在前面的一系列教程中压固，我們可以得到全國每個氣象站點(diǎn)完整的ANUSPLIN插值格式的數(shù)據(jù)怠缸，當(dāng)我們的研究范圍不在全國而在某個地理單元，如黃土高原時苞笨，如何快速的提取出來自己想要的數(shù)據(jù)呢族阅。首先選擇好自己想要的站點(diǎn)妻柒，整理成excel格式，如下所示：
image.png
然后通過代碼進(jìn)行匹配自動獲取研究區(qū)的ANUSPLIN插值格式的數(shù)據(jù)耘分，具體代碼如下：
sta=xlsread('D:\DATA\example.xlsx');
sta1=sta(:,1);
for i=1980:2014 %研究時段的開始和結(jié)束年份
filename=strcat('H:\Day1951-2012yuanshi\插值格式\1980-2016\缺測修正的全國日氣象要素\修正中國PRE',int2str(i),'.xlsx');
data=xlsread(filename);data1=data(:,1);datasum=[];
for j=1:length(sta1)
sy=find(data1==sta1(j));
if length(sy)==1
data2=data(sy,:);
datasum=[datasum;data2];
end
end
datasum=[data(1,:);datasum];
filename=strcat('H:\Day1951-2012yuanshi\插值格式\1980-2016\缺測修正的全國日氣象要素\黃土高原降水',int2str(i),'.xlsx');
xlswrite(filename,datasum)
end
輸出的結(jié)果如下圖所示，紅色部分需要自己添加表頭绑警，分別是station求泰，x，y和dem
image.png
image.png
更多需求计盒，請查看個人介紹
