1 问题
现有一土地利用分类栅格数据,如何计算出各类别所占的面积?
(来源:CSDN,中国科学数据网)
2 前言
2.1 原理
在ArcGIS中,面要素可以直接通过属性表计算几何面积,而栅格数据没办法直接计算。因为面要素是矢量数据,由具体的空间坐标构成,可以根据空间坐标来计算面积,而栅格数据由一个个像元构成,其面积只能根据相应的像元数量×像元大小来计算
2.2 前提
在计算面积前,矢量数据或栅格数据要具有投影坐标系,也就是平面坐标系,因为平面坐标系单位是米,而不是地理坐标系的单位:度,这样才能相当于在平面上计算几何面积
而且其空间参考范围要符合实际,不然在gis中计算出来的面积会与真实面积相差“十万八千里”
如果数据属性中没有投影坐标系,可以通过定义投影工具赋予其投影坐标系。由于这样会改变源数据,操作无法撤回,建议导出为一个新数据后,再为新数据定义投影
3 解决办法
3.1 重分类
利用重分类工具,输出重分类栅格后,其属性表中每个类别的count字段值就是其像元数量
3.2 属性表计算:像元数量×像元大小
在重分类栅格属性中可以查看到单位像元的大小,例如下图中所表示的像元大小=X*Y=1*1=1m²
然后打开其属性表,新建列面积字段,通过像元数量×像元大小计算其面积字段值,就得到各类别的面积啦
4 反思
如何利用ENVI快速计算出栅格数据各分类区面积?
若文章有问题或不足,欢迎批正或补充~
参考来源: