jiajialin1649 / tex Goto Github PK
View Code? Open in Web Editor NEWLicense: MIT License
License: MIT License
\documentclass{article}
\usepackage{ctex}%使用该编码可以支持中文
\usepackage{setspace}
\begin{document}
\section{geostrophic current}
{\zihao{6}大尺度的海水在压强梯度力和coriolis力平衡下的流动,流动基本是近水平的,近似认为定常的}
{\zihao{7}倾斜流的定义:海水分布均匀但是还还海面倾斜造成不均匀的水压场,在这种压力产生的地转流成为地转流}
{\heiti 倾斜流公式}
\begin{equation}
0=-\frac{1}{\rho}\frac{\partial \rho}{\partial x}+fv
\end{equation}
\begin{equation}
0=-\frac{1}{\rho}\frac{\partial \rho}{\partial y}+fu
\end{equation}
\begin{equation}
0=-\frac{1}{\rho}\frac{\partial \rho}{\partial z}-g
\end{equation}
\begin{equation}
\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}=0
\end{equation}
\begin{equation}
u\frac{\partial \theta}{\partial x}+v\frac{\partial \theta}{\partial y}+w\frac{\partial \theta}{\partial z}=0
\end{equation}
\begin{equation}
\rho=\rho(s,\theta)
\end{equation}
\begin{enumerate}
\zihao{6}\item \quad 倾斜流沿着水位线流动
\item \quad 流动方向右边为高水位
\item \quad 倾斜流从海面到海底的整个水柱具有相同的速度
\item \quad 均匀海洋中的海绵坡度可作为倾斜流大小的度量
\end{enumerate}
\heiti 梯度流公式
\begin{equation}
\frac{\partial u}{\partial t}-fv=-\frac{1}{\rho}\frac{\partial p}{\partial x}
\end{equation}
\begin{equation}
\frac{\partial v}{\partial t}-fu=-\frac{1}{\rho}\frac{\partial p}{\partial y}
\end{equation}
\begin{equation}
0=-\frac{1}{\rho}\frac{\partial p}{\partial z}-g
\end{equation}
\begin{equation}
\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}=0
\end{equation}
\begin{equation}
\frac{d\theta}{dt}=K_\theta_Z\frac{\partial^2 \theta}{\partial z^2}+K_\theta _t(\frac{\partial \theta^2}{\partial x^2}+\frac{\partial \theta^2}{\partial y^2})
\end{equation}
\begin{equation}
\frac{ds}{dt}=K_S_Z\frac{\partial^2 s}{\partial z^2}+K_S_t(\frac{\partial s^2}{\partial x^2}+\frac{\partial s^2}{\partial y^2})
\end{equation}
\begin{equation}
\rho=\rho(s,\theta,p)
\end{equation}
\begin{enumerate}
\zihao{6}\item \quad 水平流速与压强梯度垂直,梯度流沿着等压线方向流动
\item \quad 在北半球,高压在运动方向的右边;相反南半球高压在运动方向的左边
\item \quad 水平无辐散,忽略f随y的变化,梯度流沿着等温度线方向流动;沿着等密面、等盐线方向流动
\end{enumerate}
其他概念
\zihao{6}利用海兰汉森公式计算流速
流速零面
位势高度图
热成风方程的应用
\section{复习思考题}
\begin{enumerate}
\zihao{6}\item 什么是海流
\item 什么是地转流及其分类
\item 给出梯度流的定义方程以及表达式
\item 推倒方程,讨论倾斜流的特征
\item 梯度流的温度、盐度、密度的分布,说明
\item 梯度流什么情况不适用
\item 如何确定零流速面
\item 说明动力高度和流速的关系
\item 推导地转流随深度变化的公式
\end{enumerate}
\clearpage
\section{Ekman Spiral}
\end{document}
fname='subset_1_of_S2B_MSI_2021_08_28_02_39_17_T51RVP_L2R.nc';
finfo=ncinfo(fname);
b1_n=ncread(fname,'rhos_442');
b2_n=ncread(fname,'rhos_492');
b3_n=ncread(fname,'rhos_559');
b4_n=ncread(fname,'rhos_665');
b5_n=ncread(fname,'rhos_704');
b6_n=ncread(fname,'rhos_739');
b7_n=ncread(fname,'rhos_780');
b8_n=ncread(fname,'rhos_833');
b9_n=ncread(fname,'rhos_864');
b10_n=ncread(fname,'rhos_1610');
b11_n=ncread(fname,'rhos_2186');
% % %calculate ndwi20
NDWI_20=(b3_n-b10_n)./(b3_n+b10_n);%水体指数的判据
NDWI20=NDWI_20;
NDWI20(NDWI20<0)=0;
NDWI20(NDWI20>0)=1;%将水体和陆地提取出来分离
% figure
% subplot(1,2,1);
% imshow(NDWI20);
% title('NDWI水体提取','FontSize',9)
%
% subplot(1,2,2);
% rgb1=zeros(size(b3_n,1),size(b3_n,2),3);%创建了三维数组,制定好
% %数组的行数、列数、以及数组的深度,代表着数组是一个rgb图像
% rgb1(:,:,1)=adapthisteq(b10_n);
% rgb1(:,:,2)=adapthisteq(b3_n);
% rgb1(:,:,3)=adapthisteq(b4_n);
% imshow(rgb1);
% title(['20NDWI真彩图'],'FontSize',9)
%
% %水体指数2
MNDWI_21 =(b2_n-b11_n)./(b2_n+b11_n);
MNDWI21 =MNDWI_21;
MNDWI21(MNDWI21<0)=0;
MNDWI21(MNDWI21>0)=1;
%
% figure
% subplot(1,2,1);
% imshow(MNDWI20);
% title('MNDWI水体提取','FontSize',9)
% subplot(1,2,2);
% rgb2=zeros(size(b2_n,1),size(b2_n,2),3);%创建了三维数组,制定好
% %数组的行数、列数、以及数组的深度,代表着数组是一个rgb图像
% rgb2(:,:,1)=adapthisteq(b11_n);
% rgb2(:,:,2)=adapthisteq(b2_n);
% rgb2(:,:,3)=adapthisteq(b3_n);
% imshow(rgb2);
% title(['20MNDWI真彩图'],'FontSize',9)
%
% %水体指数3
% W12015_20=1.7204+171b2_n+3b6_n-70b8_n-45b10_n-71b11_n;
% W1201520=W12015_20;
% W1201520(W1201520<0)=0;
% W1201520(W1201520>0)=1;
% figure
% subplot(1,2,1);
% imshow(W1201520);
% title('W2015水体提取','FontSize',9)
% subplot(1,2,2);
% rgb3=zeros(size(b2_n,1),size(b2_n,2),3);
% rgb3(:,:,1)=adapthisteq(b8_n);
% rgb3(:,:,2)=adapthisteq(b6_n);
% rgb3(:,:,3)=adapthisteq(b2_n);
% imshow(rgb3);
% title(['w2015真彩图'],'FontSize',9)
%
% %4
% AWEINsh_20=4(b2_n-b10_n)-(0.25b8_n+2.75b11_n);
% AWEINsh20 =AWEINsh_20;
% AWEINsh20(AWEINsh20<0)=0;
% AWEINsh20(AWEINsh20>0)=1;
% figure
% subplot(1,2,1);
% imshow(AWEINsh20)
% title('AWEINsh水体提取','FontSize',9)
% subplot(1,2,2);
% rgb4=zeros(size(b2_n,1),size(b2_n,2),3);
% rgb4(:,:,1)=adapthisteq(b2_n);
% rgb4(:,:,2)=adapthisteq(b6_n);
% rgb4(:,:,3)=adapthisteq(b8_n);
% imshow(rgb4);
% title('AWEINsh真彩图','FontSize',9)
%
% %5
% AWElsh_20=b1_n+2.5b2_n-1.5(b8_n+b10_n)-0.25*b11_n;
% AWElsh20 =AWElsh_20;
% AWElsh20(AWElsh20<0)=0;
% AWElsh20(AWElsh20>0)=1;
% figure
% subplot(1,2,1);
% imshow(AWElsh20)
% title('AWEIlsh水体提取','FontSize',9)
% subplot(1,2,2);
% rgb5=zeros(size(b6_n,1),size(b6_n,2),3);
% rgb5(:,:,1)=adapthisteq(b8_n);
% rgb5(:,:,2)=adapthisteq(b6_n);
% rgb5(:,:,3)=adapthisteq(b2_n);
% imshow(rgb5);
% title('AWEIlsh真彩图','FontSize',9)
clear b9_n b8_n b7_n b6_n b5_n b4_n b3_n b2_n b1_n b10_n b11_n
fname='subset_0_of_S2B_MSI_2020_08_13_02_39_23_T51RVP_L2R.nc';
finfo=ncinfo(fname);
b1_n=ncread(fname,'rhos_442');
b2_n=ncread(fname,'rhos_492');
b3_n=ncread(fname,'rhos_559');
b4_n=ncread(fname,'rhos_665');
b5_n=ncread(fname,'rhos_704');
b6_n=ncread(fname,'rhos_739');
b7_n=ncread(fname,'rhos_780');
b8_n=ncread(fname,'rhos_833');
b9_n=ncread(fname,'rhos_864');
b10_n=ncread(fname,'rhos_1610');
b11_n=ncread(fname,'rhos_2186');
% % % % %calculate ndwi20
% % % NDWI_20=(b3_n-b10_n)./(b3_n+b10_n);%水体指数的判据
% % % NDWI20=NDWI_20;
% % % NDWI20(NDWI20<0)=0;
% % % NDWI20(NDWI20>0)=1;%将水体和陆地提取出来分离
% % % figure
% % % subplot(2,5,1);
% % % imshow(NDWI20);
% % % title('NDWI水体提取','FontSize',6)
% % %
% % % subplot(2,5,6);
% % % rgb1=zeros(size(b3_n,1),size(b3_n,2),3);%创建了三维数组,制定好
% % % %数组的行数、列数、以及数组的深度,代表着数组是一个rgb图像
% % % rgb1(:,:,1)=adapthisteq(b10_n);
% % % rgb1(:,:,2)=adapthisteq(b3_n);
% % % rgb1(:,:,3)=adapthisteq(b4_n);
% % % imshow(rgb1);
% % % title(['20NDWI真彩图'],'FontSize',6)
% % %
% % %
% % % %水体指数2
MNDWI_20 =(b2_n-b11_n)./(b2_n+b11_n);
MNDWI20 =MNDWI_20;
MNDWI20(MNDWI20<0)=0;
MNDWI20(MNDWI20>0)=1;
% % %
% % % subplot(2,5,2);
% % % imshow(MNDWI20);
% % % title('MNDWI水体提取','FontSize',6)
% % % subplot(2,5,7);
% % % rgb2=zeros(size(b2_n,1),size(b2_n,2),3);%创建了三维数组,制定好
% % % %数组的行数、列数、以及数组的深度,代表着数组是一个rgb图像
% % % rgb2(:,:,1)=adapthisteq(b11_n);
% % % rgb2(:,:,2)=adapthisteq(b2_n);
% % % rgb2(:,:,3)=adapthisteq(b3_n);
% % % imshow(rgb2);
% % % title(['20MNDWI真彩图'],'FontSize',6)
% % %
% % % %水体指数3
% % % W12015_20=1.7204+171b2_n+3b6_n-70b8_n-45b10_n-71b11_n;
% % % W1201520=W12015_20;
% % % W1201520(W1201520<0)=0;
% % % W1201520(W1201520>0)=1;
% % %
% % % subplot(2,5,3);
% % % imshow(W1201520);
% % % title('W2015水体提取','FontSize',6)
% % % subplot(2,5,8);
% % % rgb3=zeros(size(b2_n,1),size(b2_n,2),3);
% % % rgb3(:,:,1)=adapthisteq(b8_n);
% % % rgb3(:,:,2)=adapthisteq(b6_n);
% % % rgb3(:,:,3)=adapthisteq(b2_n);
% % % imshow(rgb3);
% % % title(['w2015真彩图'],'FontSize',6)
% % %
% % %
% % % AWEINsh_20=4(b2_n-b10_n)-(0.25b8_n+2.75b11_n);
% % % AWEINsh20 =AWEINsh_20;
% % % AWEINsh20(AWEINsh20<0)=0;
% % % AWEINsh20(AWEINsh20>0)=1;
% % %
% % % subplot(2,5,4);
% % % imshow(AWEINsh20)
% % % title('AWEINsh水体提取','FontSize',6)
% % % subplot(2,5,9);
% % % rgb4=zeros(size(b2_n,1),size(b2_n,2),3);
% % % rgb4(:,:,1)=adapthisteq(b2_n);
% % % rgb4(:,:,2)=adapthisteq(b6_n);
% % % rgb4(:,:,3)=adapthisteq(b8_n);
% % % imshow(rgb4);
% % % title('AWEINsh真彩图','FontSize',6)
% % %
% % %
% % % AWElsh_20=b1_n+2.5b2_n-1.5(b8_n+b10_n)-0.25*b11_n;
% % % AWElsh20 =AWElsh_20;
% % % AWElsh20(AWElsh20<0)=0;
% % % AWElsh20(AWElsh20>0)=1;
% % %
% % % subplot(2,5,5);
% % % imshow(AWElsh20)
% % % title('AWEIlsh水体提取','FontSize',6)
% % % subplot(2,5,10);
% % % rgb5=zeros(size(b6_n,1),size(b6_n,2),3);
% % % rgb5(:,:,1)=adapthisteq(b8_n);
% % % rgb5(:,:,2)=adapthisteq(b6_n);
% % % rgb5(:,:,3)=adapthisteq(b2_n);
% % % imshow(rgb5);
% % % title('AWEIlsh真彩图','FontSize',6)
%这段代码用于显示两个不同年份的MNDWI(修正归一化差异水体指数)图像之间的差异。
% 它计算了两个MNDWI图像之间的差异,
% 将结果进行了缩放,并以带有颜色条的图像形式显示差异的情况。
MNDWI_21_2=MNDWI_21.*2;
NMNDWI_20_2=MNDWI_20;
dif=MNDWI_21_2-NMNDWI_20_2;
figure;
imshow(dif);
c=colorbar;
%caxis([-1 2]);
%c.Ticks = -1:1:2;
%c.TickLabels={'土地退化为水','两年均是土地','两年均是水体','水干涸为土地'};
map = [33/255,230/255,122/255;222/255,244/255,250/255;176/255,224/255,230/255;220/255,20/255,60/255];
colormap(map);
title('两年水体土地差异');
% set( gcf , 'position' , [ 0 50 800 600 ] );
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.