将格点数据格式转化成适用于ArcMap…

%第一次尝试将自己写的小程序总结起来，都是方便查找和自我提高用

%原始截断数据u（7021*15），第一行为经度，第一列为纬度，45年月值数据

%结果数据矩阵d(year*mon+2,x*y);y=13;x=14;

%x、y为经度、纬度格点数

function(Data)=ReShape(u,year,x,y)

data=u(2:end,2:end);%仅存放数据，不存经纬度信息

for n=0:y-1%重设表头

d(1,x*n+1:(n+1)*x)=u(1,2:end);

d(2,x*n+1:(n+1)*x)=u(n+2,1);  

end

DATA=[];

%将数据重排，注重将大块的数据小块分化！这种思维很棒有没有！！！

for i=1:mon*year

    d2=[];

    for j=1:y

        temp=data(j,:);

        d2=[d2,temp];

    end

    DATA=[DATA;d2];%将一个月的小块数据练成大块~

end

d(3:end,:)=DATA;

Data=d;