分区 matlab,基于matlab的批量分区统计栅格数据集

当我们遇到要进行分区统计时如统计每个市的各类土地利用面积，或各个市的气象要素时，通常会用到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)

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