基于Matlab的中国地面气候资料日值数据集（V3.0）快速处理
基于Matlab的中国地⾯⽓候资料⽇值数据集（V3.0）快速处理
本次主要是对中国⽓象数据⽹的V3.0地⾯⽓象观测数据集进⾏快速处理，将原先分散的txt各个要素整合到⼀个列表中，⼀⾏代表⼀个站点。直接
从全国的数据集中来提取出来⾃⼰想要的站点数据集，该过程也包含了⼀年缺测数据低于30天的站点数据线性内插，同时参考数据说明对异常值
进⾏了处理。本次以新疆省内的⽓象站点为例⼦，从全国V3.0数据集中提取新疆站点。
%authoremail:**********************
station=xlsread('D: j\新疆.xlsx');% ⼀列新疆的站点数据
%
for year=2000:2017
datasum=[];
关于感恩的故事for mon=1:12
if mon<10
filename_pre=['F:\Day1951-2012yuanshi\DAY_PRE195101201612\SURF_CLI_CHN_MUL_DAY-PRE-13011-',int2str(year),'0',int2str(mon),'.TXT'];
filename_temp=['F:\Day1951-2012yuanshi\DAY_TEM195101201612\SURF_CLI_CHN_MUL_DAY-TEM-12001-',int2str(year),'0',int2str(mon),'.TXT'];
filename_rhu=['F:\Day1951-2012yuanshi\DAY_RHU195101201612\SURF_CLI_CHN_MUL_DAY-RHU-13003-',int2str(year),'0',int2str(mon),'.TXT'];
filename_SSD=['F:\Day1951-2012yuanshi\DAY_SSD195101201612\SURF_CLI_CHN_MUL_DAY-SSD-14032-',int2str(year),'0',int2str(mon),'.TXT'];
filename_WIN=['F:\Day1951-2012yuanshi\DAY_WIN195101201612\SURF_CLI_CHN_MUL_DAY-WIN-11002-',int2str(year),'0',int2str(mon),'.TXT'];
送给老师的祝福语
filename_prs=['F:\Day1951-2012yuanshi\DAY_PRS195101201612\SURF_CLI_CHN_MUL_DAY-PRS-10004-',int2str(year),'0',int2str(mon),'.TXT'];
else
filename_pre=['F:\Day1951-2012yuanshi\DAY_PRE195101201612\SURF_CLI_CHN_MUL_DAY-PRE-13011-',int2str(year),int2str(mon),'.TXT'];
filename_temp=['F:\Day1951-2012yuanshi\DAY_TEM195101201612\SURF_CLI_CHN_MUL_DAY-TEM-12001-',int2str(year),int2str(mon),'.TXT'];
filename_rhu=['F:\Day1951-2012yuanshi\DAY_RHU195101201612\SURF_CLI_CHN_MUL_DAY-RHU-13003-',int2str(year),int2str(mon),'.TXT'];格式工厂怎么用
filename_SSD=['F:\Day1951-2012yuanshi\DAY_SSD195101201612\SURF_CLI_CHN_MUL_DAY-SSD-14032-',int2str(year),int2str(mon),'.TXT'];
filename_WIN=['F:\Day1951-2012yuanshi\DAY_WIN195101201612\SURF_CLI_CHN_MUL_DAY-WIN-11002-',int2str(year),int2str(mon),'.TXT'];
filename_prs=['F:\Day1951-2012yuanshi\DAY_PRS195101201612\SURF_CLI_CHN_MUL_DAY-PRS-10004-',int2str(year),int2str(mon),'.TXT'];
end
data_prec=importdata(filename_pre);
data_temp=importdata(filename_temp);
data_rhu=importdata(filename_rhu);
data_ssd=importdata(filename_SSD);
data_win=importdata(filename_WIN);
data_prs=importdata(filename_prs);
totalsum=[];
for i=1:size(station,1)
sta=station(i);
syz_prec=find(data_prec(:,1)==sta);
syz_temp=find(data_temp(:,1)==sta);
syz_rhu=find(data_rhu(:,1)==sta);
字体大小怎么改
syz_ssd=find(data_ssd(:,1)==sta);
syz_win=find(data_win(:,1)==sta);
syz_prs=find(data_prs(:,1)==sta);
if length(syz_prec)==length(syz_temp) && length(syz_temp)==length(syz_rhu) && length(syz_rhu)==length(syz_ssd) && length(syz_ssd)==length(syz_win) && length(s                dates1=data_prec(syz_prec,1:7);
precz=data_prec(syz_prec,8:10);
tempz=data_temp(syz_temp,8:10);
rhuz=data_rhu(syz_rhu,8:9);
ssdz=data_ssd(syz_ssd,8);
winz=data_win(syz_win,8:9);
prsz=data_win(syz_prs,8);
total1=[dates1,precz,tempz,rhuz,ssdz,winz,prsz]; %1-7 8-10 11-13 14-15 16 17-18 ,19
totalsum=[totalsum;total1];
end
end
datasum=[datasum;totalsum];
end
谷歌网页翻译
%⾸先对异常数据进⾏处理
责任的作文
%台站海拔⾼度 +100000 当台站海拔⾼度为估测值时，在估测数据基础上加100000
dem=datasum(:,4);
dem(dem>100000)=dem(dem>100000)-100000;
datasum(:,4)=dem;
%各要素项  32766  数据缺测或⽆观测任务
datasum(datasum==32766)=NaN; %待插补数据
% 风速    +1000  当风速超过仪器测量上限时，在上限数据基础上加1000
winz=datasum(:,17:18);
winz(winz>1000)=winz(winz>1000)-1000;
datasum(:,17:18)=winz;

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。