当我们遇到要进行分区统计时如统计每个市的各类土地利用面积,或各个市的气象要素时,通常会用到Arcgis中的分区统计功能,但当要统计的指标特别多时,又是许多年的时候,并且需要直接一步整理成我们最后需要的格式时,Arcgis的分区统计功能可能实现起来并不方便。而采用Matlab进行分区统计可以完美的解决该问题,并且能够分区统计指标的各种指标,包括最大值,平均值等等,能够统计的范围比Arcgis更加广泛。
前期的数据准备
既然是分区统计,首先我们应该将分区的矢量文件转换了栅格数据,通常每个矢量文件的属性表中都有fid这个字段,本文假设fid是0到6,采用这个字段来进行转换为tif文件,见下图,在转栅格数据的时候,注意在环境设置时,将分辨率,范围和投影设置的和要统计的栅格数据一致。此外,所有要统计的栅格文件均需要变为tif格式的才行,Arcgis中提供将栅格数据批量转变为tif的功能。
var cpro_id = "u6885494";