一、准备工作
-
土地利用类型参考栅格和预测栅格(.tif)
需要保证两者的投影坐标系和像元大小一致
-
用户精度(user accuracy),生产者精度(producer accuracy)、总精度(overall accuracy)的计算原理
-
kappa系数计算
可以参考这一篇Kappa系数。则对于上文中的例子来说,有
-
软件需要用到Arcmap和Excel
二、处理步骤
2.1 Arcmap:栅格转点+多值提取到点
栅格转点
将参考栅格(.tif)转为点要素(.shp)
多值提取至点
2.2 Excel:计算混淆矩阵
创建数据透视表
设置数据透视表字段列表,列标签为参考栅格对应的列,行标签为预测栅格对应的列,数值选择"FID"和“POINTID”均可,但需要右键将值字段设置中的值汇总方式更改为计数
遇到这种预测栅格的可选值(1、3、4、5)小于参考栅格的情况,需要手动补全缺失分类项
利用Excel中的函数功能来计算混淆矩阵,可以得到下表