Matlab处理气象数据（十四）分段的逐点变化的DTR

最新推荐文章于 2021-03-19 14:12:57 发布
最新推荐文章于 2021-03-19 14:12:57 发布
阅读量669 收藏 4
点赞数
分类专栏： Matlab 气象 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ymyz1229/article/details/106062042
版权
气象 同时被 2 个专栏收录
23 篇文章 63 订阅
订阅专栏
Matlab
18 篇文章 70 订阅
订阅专栏 
% 两套数据逐点DTR，2000年前后对比
% 导入Tmax1,72*128*396
% 导入Tmax2,72*128*396
% 导入Tmin1,72*128*396
% 导入Tmin2,72*128*396
 
Ymax1=reshape(Tmax1,[72,128,12,33]);
Ymax1=nanmean(Ymax1,3);
Ymax1=squeeze(Ymax1); %<72x128x33 double>
 
Ymax2=reshape(Tmax2,[72,128,12,33]);
Ymax2=nanmean(Ymax2,3);
Ymax2=squeeze(Ymax2); 
 
Ymin1=reshape(Tmin1,[72,128,12,33]);
Ymin1=nanmean(Ymin1,3);
Ymin1=squeeze(Ymin1); 
 
Ymin2=reshape(Tmin2,[72,128,12,33]);
Ymin2=nanmean(Ymin2,3);
Ymin2=squeeze(Ymin2); 
 
DTR1=Ymax1-Ymin1;
DTR2=Ymax2-Ymin2;
 
DTR1=DTR1(:,:,3:33);
DTR1a=DTR1(:,:,1:19);
DTR1b=DTR1(:,:,20:31);
 
DTR2=DTR2(:,:,3:33);
DTR2a=DTR2(:,:,1:19);
DTR2b=DTR2(:,:,20:31);

%画NCEP数据2000年前的
for j=1:128
	for i=1:72
		y=DTR1a(i,j,:);
		n = 19;
		dt=1;
		% 计算统计量
		s = 0;
		for p = 1:n-1
		    for q = p+1:n
		        s = s + sign( y(q) - y(p) );
		    end
		end
		% 方差( assuming no tied groups )
		v = ( n * ( n - 1 ) * ( 2 * n + 5 ) ) / 18;
		% 检验统计量
		if s == 0
			z = 0;
			elseif s > 0,
			z = ( s - 1 ) / sqrt( v );
		else
			z = ( s + 1 ) / sqrt( v );
		end
		% should calculate Normal value here
		nor = 1.64;
		if abs( z ) < nor,
			sya(i,j) = 0;
		else
			% 计算斜率
			ndash = n * ( n - 1 ) / 2;
			s = zeros( ndash,1 );
			s=nan;
			r = 1;
			for p = 1:n-1
				for q = p+1:n
					s(r) = ( y(q) - y(p) ) / ( q - p ) / dt;
					r = r + 1;
				end
			end
			sDTR1a(i,j) = median( s );
		end
	end
end 
[x y]=meshgrid(72:0.5:135.5,18:0.5:53.5);
pcolor(x,y,sDTR1a);
shading flat;
load mycolor;
colormap(mycolor);
hold on
map=shaperead('E:\数据\边界\china_map1\maps\bou2_4l.shp');%加载省界带南海的边界线
bou2_4lx=[map(:).X];%提取经度
bou2_4ly=[map(:).Y];%提取纬度
provence=[bou2_4lx',bou2_4ly']; 
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2);%绘国界
axis([72 136 18 54]);%设置显示的经纬度范围
hold off
%最后调整区间为-0.2~0.2

%画NCEP数据2000年后的
for j=1:128
	for i=1:72
		y=DTR1b(i,j,:);
		n = 12;
		dt=1;
		% 计算统计量
		s = 0;
		for p = 1:n-1
		    for q = p+1:n
		        s = s + sign( y(q) - y(p) );
		    end
		end
		% 方差( assuming no tied groups )
		v = ( n * ( n - 1 ) * ( 2 * n + 5 ) ) / 18;
		% 检验统计量
		if s == 0
			z = 0;
			elseif s > 0,
			z = ( s - 1 ) / sqrt( v );
		else
			z = ( s + 1 ) / sqrt( v );
		end
		% should calculate Normal value here
		nor = 1.64;
		 
		if abs( z ) < nor,
			sya(i,j) = 0;
		else
			% 计算斜率
			ndash = n * ( n - 1 ) / 2;
			s = zeros( ndash,1 );
			s=nan;
			r = 1;
			for p = 1:n-1
				for q = p+1:n
					s(r) = ( y(q) - y(p) ) / ( q - p ) / dt;
					r = r + 1;
				end
			end
			sDTR1b(i,j) = median( s );
		end
	end
end 
[x y]=meshgrid(72:0.5:135.5,18:0.5:53.5);
pcolor(x,y,sDTR1b);
shading flat;
load mycolor;
colormap(mycolor);
hold on
map=shaperead('E:\数据\边界\china_map1\maps\bou2_4l.shp');%加载省界带南海的边界线
bou2_4lx=[map(:).X];%提取经度
bou2_4ly=[map(:).Y];%提取纬度
provence=[bou2_4lx',bou2_4ly']; 
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2);%绘国界
axis([72 136 18 54]);%设置显示的经纬度范围
hold off
%最后调整区间为-0.2~0.2

%画观测数据2000年前的
for j=1:128
	for i=1:72
		y=DTR2a(i,j,:);
		n = 19;
		dt=1;
		% 计算统计量
		s = 0;
		for p = 1:n-1
		    for q = p+1:n
		        s = s + sign( y(q) - y(p) );
		    end
		end
		% 方差( assuming no tied groups )
		v = ( n * ( n - 1 ) * ( 2 * n + 5 ) ) / 18;
		% 检验统计量
		if s == 0
			z = 0;
			elseif s > 0,
			z = ( s - 1 ) / sqrt( v );
		else
			z = ( s + 1 ) / sqrt( v );
		end
		% should calculate Normal value here
		nor = 1.64;
		if abs( z ) < nor,
			sya(i,j) = 0;
		else
			% 计算斜率
			ndash = n * ( n - 1 ) / 2;
			s = zeros( ndash,1 );
			s=nan;
			r = 1;
			for p = 1:n-1
				for q = p+1:n
					s(r) = ( y(q) - y(p) ) / ( q - p ) / dt;
					r = r + 1;
				end
			end
			sDTR2a(i,j) = median( s );
		end
	end
end
[x y]=meshgrid(72:0.5:135.5,18:0.5:53.5);
pcolor(x,y,sDTR2a);
shading flat;
load mycolor;
colormap(mycolor);
hold on
map=shaperead('E:\数据\边界\china_map1\maps\bou2_4l.shp');%加载省界带南海的边界线
bou2_4lx=[map(:).X];%提取经度
bou2_4ly=[map(:).Y];%提取纬度
provence=[bou2_4lx',bou2_4ly']; 
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2);%绘国界
axis([72 136 18 54]);%设置显示的经纬度范围
hold off
%最后调整区间为-0.2~0.2

%画观测数据2000年后的
for j=1:128
	for i=1:72
		y=DTR2b(i,j,:);
		n = 12;
		dt=1;
		% 计算统计量
		s = 0;
		for p = 1:n-1
		    for q = p+1:n
		        s = s + sign( y(q) - y(p) );
		    end
		end
		% 方差( assuming no tied groups )
		v = ( n * ( n - 1 ) * ( 2 * n + 5 ) ) / 18;
		% 检验统计量
		if s == 0
			z = 0;
			elseif s > 0,
			z = ( s - 1 ) / sqrt( v );
		else
			z = ( s + 1 ) / sqrt( v );
		end
		% should calculate Normal value here
		nor = 1.64;
		if abs( z ) < nor,
			sya(i,j) = 0;
		else
			% 计算斜率
			ndash = n * ( n - 1 ) / 2;
			s = zeros( ndash,1 );
			s=nan;
			r = 1;
			for p = 1:n-1
				for q = p+1:n
					s(r) = ( y(q) - y(p) ) / ( q - p ) / dt;
					r = r + 1;
				end
			end
			sDTR2b(i,j) = median( s );
		end
	end
end
[x y]=meshgrid(72:0.5:135.5,18:0.5:53.5);
pcolor(x,y,sDTR2b);
shading flat;
load mycolor;
colormap(mycolor);
hold on
map=shaperead('E:\数据\边界\china_map1\maps\bou2_4l.shp');%加载省界带南海的边界线
bou2_4lx=[map(:).X];%提取经度
bou2_4ly=[map(:).Y];%提取纬度
provence=[bou2_4lx',bou2_4ly']; 
plot(bou2_4lx,bou2_4ly,'-k','LineWidth',1.2);%绘国界
axis([72 136 18 54]);%设置显示的经纬度范围
hold off
%最后调整区间为-0.2~0.2

得到如下结果：

在这里插入图片描述

NCEP数据1981-1999年（左）和2000-2011年（右）的DTR逐点趋势分析

在这里插入图片描述

观测数据1981-1999年（左）和2000-2011年（右）的DTR逐点趋势分析

相关链接：
Matlab处理气象数据——目录

岩鸣杨子
关注
  • 0
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
离散系数的计算公式_数据的离散程度分析1
weixin_29424607的博客
01-28 9811
在统计学中，把反映现象总体中各个体的变量值之间差异程度的指标称为离散程度，也称为离中趋势。描述一组数据离散程度常用极差、四分位差、方差和标准差、变异系数等。极差极差(Range)也叫全距，是一组数据中最大值与最小值之差，即：R=max(Xi)-min(Xi)某班级40名同学统计学的考试成绩原始资料如表求考试成绩极差。在A11单元格输入表格名称“极差”。单击B11单元格，在编辑栏中输入“=MAX(A...
Matlab处理气象数据（九）逐点变化的显著性分析
设计の诗
05-11 3686
利用上面M-K趋势分析的结果，只让置信度大于95%的数据在地图上显示。 这里，可以先自定义一个colorbar。 mycolor = colormap; %mycolor名称任意 save mycolor mycolor; 下次画图时，只要先load mycolor; 再输入colormap(mycolor) 就可以直接用这个colorbar了。 % 逐点变化显著性 % 导入数据Y1、Y2 Y=Y2-Y1; Y=Y(:,:,3:35); YA=Y(:,:,1:19); YB=Y(:,:,20:33);
Matlab处理气象数据（十二）DTR及其异常的计算
设计の诗
05-11 2227
DTR（diurnal temperature range）是日较差，即最高温和最低温度的差值。一段时期平均DTR的增加或降低可以反映出一个区域内的气候变化状况。 load('Temmax1.mat'); load('Temmax2.mat'); load('Temmin1.mat'); load('Temmin2.mat'); DTR1=Temmax1-Temmin1; DTR2=Temmax2-Temmin2; %画两套数据的DTR plot(DTR1,'r.-','linewidth',2); ho
MK检验程序出错有请大侠-MK.m
08-13
MK检验程序出错有请大侠-MK.m % Time Series Trend Detection Tests % %     [ z, sl, lcl, ucl ] = trend % % where   z   = Mann-Kendall Statistic %         sl  = Sen's Slope Estimate %         lcl = Lower Confidence Limit of sl %         ucl = Upper Confidence Limit of sl %         y   = Time Series of Data %         dt  = Time Interval of Data % % Bob Newell, February 1996 % %-------------------------------------------------- % load unnamed1.mat; y=unnamed1; function [ z, sl, lcl, ucl ] = trendMK % n = length; %-------------------------------------------------- % Mann-Kendall Test for N > 40 % disp; if n < 41,   disp; end; % calculate statistic s = 0; for k = 1:n-1,   for j = k 1:n,     s = s sign - y );   end; end; % variance v = * ) / 18; % test statistic if s == 0,   z = 0; elseif s > 0,   z = / sqrt; else   z = / sqrt; end;             % should calculate Normal value here nor = 1.96; % results disp ] ); disp ) ] ); disp ] ); if abs < nor,   disp;   z = 0; elseif z > 0,   disp; else   disp; end; %---------------------------------------------------- %             disp; % calculate slopes ndash = n * / 2; s = zeros; i = 1; for k = 1:n-1,   for j = k 1:n,     s = - y ) / / dt; i = i 1;   end; end; % the estimate sl = median; disp ] ); % variance v = * ) / 18; m1 = fix ) / 2 ); m2 = fix ) / 2 ); s = sort; lcl = s; ucl = s; disp( [ '     Lower Confidence Limit = ' ...                              num2str ] ); disp( [ '     Upper Confidence Limit = ' ...                            num2str ] ); %---------------------------------------------------- % the end
CRU2GDL:Matlab代码计算每个行政单位的平均温度-matlab开发
05-31
Matlab代码，用于计算每个管理单位的平均温度。 可以轻松调整其他多边形（如果您有形状文件）或其他气候变量。
Matlab处理气象数据（十三）逐点变化的DTR及其异常
设计の诗
05-11 805
% 两套数据逐点DTR % 导入Tmax1,72*128*396 % 导入Tmax2,72*128*396 % 导入Tmin1,72*128*396 % 导入Tmin2,72*128*396 Ymax1=reshape(Tmax1,[72,128,12,33]); Ymax1=nanmean(Ymax1,3); Ymax1=squeeze(Ymax1); %<72x128x33 double> Ymax2=reshape(Tmax2,[72,128,12,33]); Ymax2=nanm
DTR KML:DTR-matlab开发
05-30
总的来说，"DTR KML:DTR-matlab开发"可能是一个用于处理地理空间数据、构建决策树模型并以KML格式输出的MATLAB项目。其提供的出色输出表明，该项目在数据处理、模型构建和地图可视化方面具有高效和高质量的特点。...
数字信号处理基本知识点Matlab实现.docx
07-10
《数字信号处理基本知识点Matlab实现》 在数字信号处理领域，MATLAB是一个不可或缺的工具。MATLAB由美国MathWorks公司开发，自1984年发布以来，因其强大的数值计算和可视化功能，逐渐成为科研和工程计算的重要平台...
DAQ Through Parallel Port：使用并口采集和发送数据的简单软件-matlab开发
05-31
在使用MATLAB进行并行端口操作时，确保对硬件兼容性、数据格式转换和错误处理有深入理解，这样可以避免潜在的问题。同时，如果可能，应考虑使用更现代的DAQ设备和接口，如PCI、GPIB或USB，它们通常提供更高的速度和...
基于matlab实现编写的串口调试工具，数据接收部分采用中断方式，保证了实时的数据显示
最新发布
05-01
该工具的设计者选择了中断方式来处理数据接收，这是为了确保数据能够实时地被处理和显示。中断方式下，每当串口接收到新的数据，MATLAB就会触发一个中断事件，从而立即执行相应的处理函数，而不必持续轮询，提高了...
OS85650：I／O端口扩展／DTCP协处理器.pdf
09-25
所有器件都具备USB暂停、恢复和远程唤醒功能，具有128字节TX和384字节RX的FIFO，优化了整体数据吞吐量，并具备自动硬件流控制（RTS/CTS或DTR/DSR）和软件流控制（XON/XOFF）。器件无需外部晶体，节省了成本，其睡眠...
matlab指数转换,基于Matlab的干旱指数PDSI及CRU全球气象nc数据的处理过程-由nc文件转变为tif文件...
weixin_39913628的博客
03-16 523
image.png从该网站上可直接下载全球1901-2016年的逐月PDSI数据，文件为nc格式，为进一步和其他栅格数据进行计算，需要将nc文件转变为tif文件，因此本文提供一种能够批量转换的处理方式。首先准备个样例数据ncdisp('H:\Global\PDSI\scPDSI.cru.3.25.bams2017.GLOBAL.1901.2016.nc');data1=ncread('H:\Glo...
Matlab处理气象数据（十五）分省份的趋势分析
设计の诗
05-11 3010
用省份边界数据，制作掩膜（mask），从而提取出各省的DTR数据。并且求P和R2。 %求P值 t=1:35; t=t'; t=[ones(35,1),t]; [b,bint,r,rint,stats]= regress(NCEP,t,0.05); %结果： alpha为显著性水平(缺省时设定为0.05) stats是用于检验回归模型的统计量有四个数值 R2，其中R是相关系数 F统计量值 概率P，与统计量F对应 an estimate of the error variance（一个错误的方差估计） %回归
RTS与CTS的含义
weixin_30556161的博客
10-26 2100
====================================我是分割线首先介绍下网上看到的======================================================================================== 一、RS232标准中的RTS与CTS RTS,CTS------请求发送/清除发送，用于半双工时的收发切换，属于辅助...
Matlab处理气象数据（十七）季节变化
设计の诗
05-11 5687
17.1 各省冬季和夏季的DTR趋势 以江西为例，其他省份只需将下面代码中的“江西”替换为其他省份的名称即可。 %% clc,clear load('E:\数据\201507\各省MASK\江西');%MASK是各省掩膜 load('E:\数据\201501\Tmax1');%导入NCEP数据最高温度 72*128*396 load('E:\数据\201501\Tmax2');%导入观测数据最高温度 72*128*396 load('E:\数据\201501\Tmin1');%导入NCEP数据最低温度 72
matlab写错怎么删除,这个程序是否写错？
weixin_30615567的博客
03-19 502
%还是直接上代码% Time Series Trend Detection Tests% [ z, sl, lcl, ucl ] = trend( y, dt )% where z = Mann-Kendall Statistic% sl = Sen's Slope Estimate% lcl = Lower Confidence Limit of sl% ucl = Upper Confiden...
Matlab处理气象数据——目录
热门推荐
设计の诗
05-10 1万+
整理了一下研究生期间做的一些气象数据处理方面的资料，讲研究过程和matlab代码一并贴出来，以便存档和查阅。 论文原文见知网 目录 1.读取NCEP数据 2.调用NCEP数据，进行预处理 3.NCEP数据的简单计算 4.观测数据的预处理和计算 5.NCEP数据的剪裁、插值和两套数据的对比 6.数据的面积加权 7.分段趋势线的做法 8.数据的逐点变化趋势 8.1 分别对NCEP和观测数据进行逐点趋势分析 8.2 平均、最高、最低温度差值的趋势 8.3 退耕还林前后的温度差值对比 9
Matlab处理气象数据（一）读取NCEP数据
设计の诗
05-10 1万+
数据获取 NCEP数据air.2m.mon.mean.nc，下载地址 读取 使用matlab命令读取NCEP数据。 常用的读取NC头文件的命令是 nc_dump 和 nc_dump 。 执行nc_dump命令： nc_dump('E:\数据\空间数据\NECP\air.2m.mon.mean.nc ') 得到如下结果： NetCDF-3 Classic E:\数据\空间数据\NECP\air.2m.mon.mean.nc { dimensions: lon = 192 ; lat = 94 ; le
Matlab处理气象数据（四）观测数据的预处理和计算
设计の诗
05-10 1万+
观测数据为中国气象数据网上获得的中国地面温度月值0.5°×0.5°格点数据集，时间范围是1961年1月至2013年12月。这套数据为txt格式，包含头文件。头文件信息为： NCOLS 128 NROWS 72 XLLCORNER 72.000000000000 YLLCORNER 18.000000000000 CELLSIZE 0.50000000000000 NODATA_VALUE -99.00000 从中可以得知，有128
串口DTR怎么通过matlab管理
03-31
您可以使用MATLAB的serial函数来管理串口DTR。以下是一些可能有用的示例代码： 1. 设置串口对象： ```matlab s = serial('COM1'); % 替换成您的串口号 set(s, 'BaudRate', 9600); % 设置波特率 set(s, 'DataBits',...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值