当我们遇到要进行分区统计时如统计每个市的各类土地利用面积，或各个市的气象要素时，通常会用到Arcgis中的分区统计功能，但当要统计的指标特别多时，又是许多年的时候，并且需要直接一步整理成我们最后需要的格式时，Arcgis的分区统计功能可能实现起来并不方便。而采用Matlab进行分区统计可以完美的解决该问题，并且能够分区统计指标的各种指标，包括最大值，平均值等等，能够统计的范围比Arcgis更加广泛。

前期的数据准备

既然是分区统计，首先我们应该将分区的矢量文件转换了栅格数据，通常每个矢量文件的属性表中都有fid这个字段，本文假设fid是0到6，采用这个字段来进行转换为tif文件，见下图，在转栅格数据的时候，注意在环境设置时，将分辨率，范围和投影设置的和要统计的栅格数据一致。此外，所有要统计的栅格文件均需要变为tif格式的才行，Arcgis中提供将栅格数据批量转变为tif的功能。

image.png

matlab程序

 %%统计每个流域1982-2015年每年的降水，生长季平均气温和生长季辐射
basin=importdata('D:\qixiang\重点流域范围 - 副本\majorbasin_new_new.tif');
datatotal=[];
for year=1982:2015
    filenamepre=strcat('D:\qixiang\年尺度数据\0.1度\采样成8km\',int2str(year),'year_prec.tif');
    pre=importdata(filenamepre);
   
    filenamepre=strcat('D:\qixiang\年尺度数据\0.1度\采样成8km\',int2str(year),'年生长季总辐射.tif');
    rad_grow=importdata(filenamepre);
    
    filenamepre=strcat('D:\qixiang\年尺度数据\0.1度\采样成8km\',int2str(year),'年生长季平均风速.tif');
    win_grow=importdata(filenamepre);
    
    datasum=[];
    for j=0:6
        sy=find(basin==j);
        % 
        pre_value=pre(sy);
        pre_value=sum(pre_value(pre_value>=0))./length(sy);
        
        win_grow_value=win_grow(sy);
        win_grow_value=sum(win_grow_value(win_grow_value>=0))./length(sy);
        
        rad_grow_value=rad_grow(sy);
        rad_grow_value=sum(rad_grow_value(rad_grow_value>=0))./length(sy);
         
        datanew=[pre_value,win_grow_value,rad_grow_value];
        datasum=[datasum,datanew];
    end
        datatotal=[datatotal;datasum];
end
xlswrite('各个流域全部要素.xlsx',datatotal)