Matlab处理气象数据（十一）数据的异常值计算

%平均温度的异常值计算
load('Tem1.mat');%导入NCEP数据的面积加权年平均
load('Tem2.mat');%导入观测数据的面积加权年平均
m1=mean(Tem1); %求Tem1的平均值
m2=mean(Tem2);
a1=Tem1-m1; %求Tem1的距平
a2=Tem2-m2;
plot(a1,'r.-','linewidth',2);%画NCEP数据年平均气温折线图，红色实线实心点
hold on
plot(a2,'b.-','linewidth',2);
axis([ -inf inf -1.5 1.5]); %设置纵坐标范围
legend({'NCEP','Observed'},'Location','Northwest');%添加图例
set(gca,'xtick',[2 7 12 17 22 27 32],'xticklabel',{'1980','1985','1990','1995','2000','2005','2010'});%在x轴特定位置上添加标注
set(gca, 'FontSize',10,'FontWeight','Bold','tickdir','out') %设置标注为10号字、加粗、标记线向外
h=xlabel('Year'); %设置x轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
h=ylabel('Temperarure(\circC)'); %设置y轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
xlim([1 35])%x轴范围锁定为1~35
box off %去掉外框
hold off

NCEP数据和观测数据平均温度的异常值

%最高温度的异常值计算
load('Temmax1.mat');
load('Temmax2.mat');
amax1=Temmax1-mean(Temmax1); %求Temmax1的距平
amax2=Temmax2-mean(Temmax2);
plot(amax1,'r.-','linewidth',2);%画NCEP数据最高气温折线图，红色实线实心点
hold on
plot(amax2,'b.-','linewidth',2);
axis([ -inf inf -1.5 1.5]); %设置纵坐标范围
legend({'NCEP','Observed'},'Location','Northwest');%添加图例
set(gca,'xtick',[2 7 12 17 22 27 32],'xticklabel',{'1980','1985','1990','1995','2000','2005','2010'});%在x轴特定位置上添加标注
set(gca, 'FontSize',10,'FontWeight','Bold','tickdir','out') %设置标注为10号字、加粗、标记线向外
h=xlabel('Year'); %设置x轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
h=ylabel('Temperarure(\circC)'); %设置y轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
xlim([1 35])%x轴范围锁定为1~35
box off %去掉外框
hold off

NCEP数据和观测数据最高温度的异常值

%最低温度的异常值计算
load('Temmin1.mat');
load('Temmin2.mat');
amin1=Temmin1-mean(Temmin1); %求Temmin1的距平
amin2=Temmin2-mean(Temmin2);
plot(amin1,'r.-','linewidth',2);%画NCEP数据最低气温折线图，红色实线实心点
hold on
plot(amin2,'b.-','linewidth',2);
axis([ -inf inf -1.5 1.5]); %设置纵坐标范围
legend({'NCEP','Observed'},'Location','Northwest');%添加图例
set(gca,'xtick',[2 7 12 17 22 27 32],'xticklabel',{'1980','1985','1990','1995','2000','2005','2010'});%在x轴特定位置上添加标注
set(gca, 'FontSize',10,'FontWeight','Bold','tickdir','out') %设置标注为10号字、加粗、标记线向外
h=xlabel('Year'); %设置x轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
h=ylabel('Temperarure(\circC)'); %设置y轴名称
set(h, 'FontSize',10,'FontWeight','Bold')
xlim([1 35])%x轴范围锁定为1~35
box off %去掉外框
hold off

NCEP数据和观测数据最低温度的异常值

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

Matlab处理气象数据（十一）数据的异常值计算相关推荐

Matlab处理气象数据（十）给地图添加边界线
边界线数据从网上下载得来: 其中,"l"结尾的为只有线没有面的边界线,"p"结尾的为实心的带边界的平面.为了让数据显现出来,我们使用"bou2_4l. ...
Matlab处理气象数据（四）观测数据的预处理和计算
观测数据为中国气象数据网上获得的中国地面温度月值0.5°×0.5°格点数据集,时间范围是1961年1月至2013年12月.这套数据为txt格式,包含头文件.头文件信息为: NCOLS 128 NROW ...
Matlab处理气象数据（七）分段趋势线的做法
退耕还林工程自1999年起开始实施,土地利用的改变会影响区域气候.为了研究区域气候的变化,需要对退耕前和退耕后作对比. 在matlab中,需要对1999年前和后分别作出趋势线,这种趋势线可以采用mat ...
Matlab处理气象数据（九）逐点变化的显著性分析
利用上面M-K趋势分析的结果,只让置信度大于95%的数据在地图上显示. 这里,可以先自定义一个colorbar. mycolor = colormap; %mycolor名称任意 save mycol ...
美赛整理之Matlab读取全球海洋温度数据并显示干货
Matlab读取全球海洋温度数据并显示干货 Matlab读取全球海洋温度数据并显示干货 Matlab读取全球海洋温度数据并显示干货一.nc文件的读取二.画出从1981到2000年的全球温度海洋变化 ...
Matlab Smooth函数/丝滑数据
在处理数据的时候,一般会进行数据拟合,这时候免不了会有噪点,从而造成图像上的"飞点".幸好Matlab自带一个平滑函数smooth,笔者也是最近才晓得,做点笔记. 在Matlab ...
Matlab/Python nc格式小时数据转日期
在处理气象nc格式数据时,时间变量往往不是已YYYYMMDD格式给出日期,而是距某一天(如1900年1月1日)XXX小时或XXX日. 为此,可由以下代码将nc的时间变量转为日期: matlab代码: ...
机器学习数据预处理之离群值/异常值：图像对比法
机器学习数据预处理之离群值/异常值:图像对比法 garbage in,garbage out. 异常值是分析师和数据科学家常用的术语,因为它需要密切注意,否则可能导致错误的估计. 简单来说,异常值是一 ...
机器学习数据预处理之离群值/异常值：箱图法（boxplot method）
机器学习数据预处理之离群值/异常值:箱图法(boxplot method) garbage in,garbage out. 异常值是分析师和数据科学家常用的术语,因为它需要密切注意,否则可能导致错误的 ...

Matlab处理气象数据（十一）数据的异常值计算

Matlab处理气象数据（十一）数据的异常值计算相关推荐

最新文章

热门文章