根据你的数据类型选用不同的数据生成方法
若你的数据是
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
图层上。统计出面积,依进行,就可以得到土地利用类型转移
矩阵
。
若你的数据是
1 Erdas
做这一步之前记得先对两时相的数据进行重编码(
一般运行如果出现错误
2
NC(I,J)=NC(I)*10+NC(J),(J>I)
其中:
在此计算的基础上,将以上变化影像图转化为
若数据是
1 Erdas
此时注意
做这一步之前,请做好前期的地理编码。
2 ArcView3.3
把两时相的
3
叠加后的文件在