- 积分
- 1019
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-7-15
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我们通常下载的为逐日的气象数据,需要处理为不同站点的文本,程序如下:clc;
clear;
datadir='F:\ChinaStations\'
files= dir(fullfile(datadir,'*.txt'));
n=length(files);
StaName_one=ones(0,1)*nan;
StaName_two=ones(0,1)*nan;
stationsnames=[50136;50468;50527;50742;50834;50873;50953;51076;51133;51431;51463;51567;51573;51628;51709;51777;51828;52203;52267;52418;52533;52681;52754;52818;52866;52889;52983;53068;53336;53487;53543;53614;53772;53817;53845;53923;53963;54102;54135;54161;54292;54324;54342;54511;54527;54539;54662;54764;54823;54936;55228;55299;55591;56029;56043;56137;56146;56173;56187;56196;56294;56385;56651;56666;56739;56778;56959;56985;57036;57083;57131;57178;57245;57461;57494;57516;57602;57604;57649;57687;57816;57874;57957;57993;58141;58144;58208;58238;58265;58321;58362;58457;58531;58606;58665;58737;58847;59287;59316;59431;59644;59758;59948;59981];
num=length(stationsnames);
for j=1:num
MerData=ones(0,50)*nan;
data_temp=ones(0,50)*nan;
for i=1:n
filename=files(i).name;
data=textread([datadir,filename]);
idx=find(data(:,1)==stationsnames(j));
MerData=data(idx,:);
data_temp=[data_temp;MerData];
end
xlswrite(['F:\AlbedoAllStations\',num2str(stationsnames(j)),'.xlsx'],data_temp);
end
|
|