在前面的一系列教程中，我们可以得到全国每个气象站点完整的ANUSPLIN插值格式的数据，当我们的研究范围不在全国而在某个地理单元，如黄土高原时，如何快速的提取出来自己想要的数据呢。首先选择好自己想要的站点，整理成excel格式，如下所示：
image.png
然后通过代码进行匹配自动获取研究区的ANUSPLIN插值格式的数据，具体代码如下：
sta=xlsread('D:\DATA\example.xlsx');
sta1=sta(:,1);
for i=1980:2014 %研究时段的开始和结束年份
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
输出的结果如下图所示，红色部分需要自己添加表头，分别是station，x，y和dem
image.png
image.png
更多需求，请查看个人介绍
