在前面的一系列教程中，我們可以得到全國每個(gè)氣象站點(diǎn)完整的ANUSPLIN插值格式的數(shù)據(jù)，當(dāng)我們的研究范圍不在全國而在某個(gè)地理單元，如黃土高原時(shí)，如何快速的提取出來自己想要的數(shù)據(jù)呢。首先選擇好自己想要的站點(diǎn)，整理成excel格式，如下所示：

image.png

然后通過代碼進(jìn)行匹配自動(dòng)獲取研究區(qū)的ANUSPLIN插值格式的數(shù)據(jù)，具體代碼如下：

sta=xlsread('D:\DATA\example.xlsx');
sta1=sta(:,1);
for i=1980:2014 %研究時(shí)段的開始和結(jié)束年份
    filename=strcat('H:\Day1951-2012yuanshi\插值格式\1980-2016\缺測(cè)修正的全國日氣象要素\修正中國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\缺測(cè)修正的全國日氣象要素\黃土高原降水',int2str(i),'.xlsx');
    xlswrite(filename,datasum)
end

輸出的結(jié)果如下圖所示，紅色部分需要自己添加表頭，分別是station，x，y和dem

image.png

image.png

更多需求，請(qǐng)查看個(gè)人介紹