月值数据转换

月值数据处理，除了日值气象要素外，我们还经常用到月值数据，针对月值气象要素如何批量变成AUSPLIN需要的格式呢。
月值数据如下：

image.png

上图第一列为气象站点，第二列为月份，第三列为年份，第四列为气象要素
需要转换成以下形式：

image.png

代码如下所示：

a=xlsread('D:\日降水插值\中国气象站点shp\china-station1.xlsx'); %已经带有全部站点的变成以m为单位的经纬度和高程
filename='D:\temp-L1\new\数据\中国插值格式的'; %存放结果数据的位置
station=a(:,1);x=a(:,2);y=a(:,3);dem=a(:,4);%站点，经度，纬度，高程;
e=dir(fullfile('*.xlsx'));% 要转换的excel数据，见图1
for i=1:size(e,1)
    a=xlsread(e(i).name);
    st=a(:,1);year=a(:,3);month=a(:,2);termtemmax=a(:,4:end);
    sta=unique(st);
    datatotal=[];
    for j=1:size(sta,1)
        sta1=sta(j);
        sy=find(st==sta1);
        datasum=[];
        for mon=1:12
            symon=find(mon==month);
            sy_inter=intersect(symon,sy);
            if length(sy_inter)<1
                value=NaN;
                datasum=[datasum,value];
            else
                datasum=[datasum,termtemmax(sy_inter)];
            end
        end
        datatotal=[datatotal;datasum];
    end
    datatotal=[sta,datatotal];
   
    zz1=[];
    for n=1:length(sta)
        st1=sta(n);
        data2=datatotal(n,:);
        sy=find(station==st1);
        if length(sy)==1
            zz=[x(sy),y(sy),dem(sy),data2];
            zz1=[zz1;zz];
        end
    end
    %除掉存在NaN值的
    datasum=[];
    for kkk=1:size(zz1,1)
        value=zz1(kkk,:);
        if sum(abs(value)>=0)==16;
           datasum=[datasum;value];
        end
    end
    datasum=[datasum(:,4),datasum(:,1:3),datasum(:,5:end)];
    datasum=[[NaN,NaN,NaN,NaN,[1:12]];datasum];
    name1=e(i).name;
    name1=strcat(filename,name1);
    xlswrite(name1,datasum);
end

通过上述代码即可以得到图2，完成月份数据的转换