土地利用转移矩阵图怎么做_土地利用转移矩阵生成的几种方法

土地利用转移矩阵生成的几种方法

土地利用转移矩阵生成的几种方法

根据你的数据类型选用不同的数据生成方法

若你的数据是

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

图层上。统计出面积，依进行，就可以得到土

地利用类型转移矩阵

。