遥感专题讲座——影像信息提取(六、土地利用类型转换混淆矩阵)

土地利用类型转换（变化）混淆矩阵

查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。

土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI 下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。

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 变化掩膜图像