jiajialin1649 / tex Goto Github PK

\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 ] );