再议ANUSPLIN之从全国ANUSPLIN插值格式站点中提取研究区站点

在前面的一系列教程中，我们可以得到全国每个气象站点完整的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

更多需求，请查看个人介绍