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

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

若你的数据是 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 图层上。统计出面积，依进行，就可以得到土地利用类型转移 矩阵 。