ENVI4.8下从两幅分类结果的栅格图中计算土地利用类型转换矩阵
1、准备数据 两个时相的土地利用分类结果,它是单波段、专题类型的伪彩色图像(ENVI Classification)
2、计算转换矩阵 打开两个土地利用分类结果
(1) 在主菜单中,选择Basic Tools → Change Detection → Change Detection Statistics
(2)分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果
(3) 在Define Equivalent Classes对话框中(图1),如果两个土地利用分类名称一致,系统自动将Initial State Class和Final State Class对应,否则手动选择,单击Add Pair按钮选择
(4) 选择对应的地物类型之后,单击OK按钮,出现图2对话框。选择生成图表表示单位(Report Type):像素(Pixels)、百分比(Percent)和面积(Area)。选择Output Classification Mask Images?为YES,输出掩膜图像,选择输入路径及文件名
(5) 单击OK,执行土地利用类型转换矩阵计算过程。
图1 Define Equivalent Classes对话框
图2 选择数据参数
3、查看结果
(1)如图3为得到的土地利用类型转换矩阵结果。横字段表示前一时间段(Initial State)的土地利用类别,纵字段为后一时间段(Final State)的土地利用类别。横字段和纵字段交叉处表示变化值,如有2520900平方米林地用地变化为草地。
图3 土地利用类型转换矩阵
(2) 还可以为每一个地类生成一个变换掩膜图像,图4所示为其中一个地类的掩膜图像。掩膜图像的灰度值表示变化类型,如这里的2{草地}表示林地变化为草地的像元。
图4 变化掩膜图像