土地利用类型转换矩阵可以从两幅栅格图中计算得到,也可以从两个矢量文件中计算获得。下面介绍在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,执行土地利用类型转换矩阵计算过程。
![点击查看大图](http://www.upgis.com/bbs/attachments/month_0910/09101320064b1e1b7c8c2b0148.jpg)
图1 Define Equivalent Classes对话框
![点击查看大图](http://www.upgis.com/bbs/attachments/month_0910/091013200625c0fb050f59114f.jpg)
图2 选择数据参数
3、查看结果
(1)如图3为得到的土地利用类型转换矩阵结果。横字段表示前一时间段(Initial State)的土地利用类别,纵字段为后一时间段(Final State)的土地利用类别。横字段和纵字段交叉处表示变化值,如有2520900平方米林地用地变化为草地。
![点击查看大图](http://www.upgis.com/bbs/attachments/month_0910/09101320065dccad25cdb2b688.jpg)
图3 土地利用类型转换矩阵
(2) 还可以为每一个地类生成一个变换掩膜图像,图4所示为其中一个地类的掩膜图像。掩膜图像的灰度值表示变化类型,如这里的2{草地}表示林地变化为草地的像元。
![点击查看大图](http://www.upgis.com/bbs/attachments/month_0910/0910132006db8be3027972a772.jpg)
图4 变化掩膜图像