好文转载收藏，原文来自：http://hi.baidu.com/dengshubin/blog/item/d4e04e4ce8e5d2f0d62afc68.html

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

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