土地利用转移矩阵生成的几种方法
土地利用转移矩阵生成的几种方法
根据你的数据类型选用不同的数据生成方法
若你的数据是
Raster
格式:则有如下方法
1
Erdas
Imagine----Interpreter---Gis
Analysis---Matrix
,输入两个时相的
Raster
数据即可
做这一步之前记得先对两时相的数据进行重编码(
nterpreter---Gis
Analysis---Recode
)
一般运行如果出现错误
肯定是重编码没做好,请继续查证。
2
先在
Erdas
中利用
Modeler
计算
如下公式
NC(I,J)=NC(I)*10+NC(J),(J>I)
其中:
NC(I,J)
表示
i
,
j
两年份的土地利用变化图;
NC
(
i
)表示
i
年份遥感
分类影像;
NC
(
j
)表示
j
年份的遥感分类影像。
在此计算的基础上,将以上变化影像图转化为
BIL
格式,再利用
ARC/INFO
GRID
模块将影像转为
GRID
格式,然后利用
GRID
模块中的属性
表(
vat
)查看命令对影像灰度值进行统计,最后得出土地利用转化举
证。(注
:
此方法本人尚未实现过,不知可行否)。
若数据是
Vector
格式
1
Erdas
Imagine----Interpreter---Gis
Analysis---Matrix
,输入两个时相的
Vector
数据即可
此时注意
输出栅格大小不应设的太小
要不一运行就会提示你的空间不
足
做这一步之前,请做好前期的地理编码。
2 ArcView3.3
加载
spatial analysis
模块
把两时相的
Vector
图转成
grid
格式(当然中间有一些单位的设置根据你
做的图的分辨率来设置即可)
analysis---mapcaculate
直接计算即可。
3
把两期解译完的
Vector
文件在
arctoolbox——overlay——union
中叠加,
注意:两个文件不能用同一个字段名,比如一个用
93Type
,另一个时相
则用
00Type
叠加后的文件在
Arcmap
中打开,选中文件,然后点右键
——Property
——
空间查询,输入条件语句,比如:
93Type=‘1’And
00Type=‘2’
;查询
结果即为第一种类型转化为第二种类型的图形,可以另建一图层比如:
12
,把查询结果复制到
12
图层上。统计出面积,依进行,就可以得到土
地利用类型转移矩阵
。