区域栅格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