从同学空间找的,希望对你有所帮助吧
土地利用转移矩阵生成的几种方法：
据你的数据类型选用不同的数据生成方法
若你的数据是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图层上.统计出面积,依进行,就可以得到土地利用类型转移矩阵 .
4.ENVI下操作.
ENVI下有个change detection工具,很easy...可以很方便的计算出转移矩阵.
5.Idiris软件下
此软件下具有更加完善的土地利用分析工具.而且制图比较漂亮.
而更具有意义的是,Idiris下对土地利用转移矩阵做了更深入的应用.
可以利用计算的土地利用转移矩阵,结合markov方法进行土地利用模拟预测.
just try it ...
6.嗯.another one...Arcgis下有个Tabu?什么的工具,（搜索一下）,可以直接对其进行统计计算.
