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

《土地利用转移矩阵生成的几种方法.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、土地利用转移矩阵生成的几种方法土地利用转移矩阵生成的几种方法根据你的数据类型选用不同的数据生成方法若你的数据是Raster格式：则有如下方法1ErdasImagine----Interpreter---GisAnalysis---Matrix，输入两个时相的Raster数据即可做这一步之前记得先对两时相的数据进行重编码（nterpreter---GisAnalysis---Recode）一般运行如果出现错误肯定是重编码没做好，请继续查证。2先在Erdas中利用Modeler计算如下公式NC(I,J)=NC(I)*10+NC(J),(J>I)其中：NC(I,J)表示i，j两年份的土地利用变

2、化图；NC（i）表示i年份遥感分类影像；NC（j）表示j年份的遥感分类影像。在此计算的基础上，将以上变化影像图转化为BIL格式，再利用ARC/INFOGRID模块将影像转为GRID格式，然后利用GRID模块中的属性表（vat）查看命令对影像灰度值进行统计，最后得出土地利用转化举证。（注:此方法本人尚未实现过，不知可行否）。若数据是Vector格式1ErdasImagine----Interpreter---GisAnalysis---Matrix，输入两个时相的Vector数据即可此时注意输出栅格大小不应设的太小要不一运行就会提示你的空间不足做这一步之前，请做好前期的地理编码。2ArcV

3、iew3.3加载spatialanalysis模块把两时相的Vector图转成grid格式（当然中间有一些单位的设置根据你做的图的分辨率来设置即可）analysis---mapcaculate直接计算即可。3把两期解译完的Vector文件在arctoolbox——overlay——union中叠加，注意：两个文件不能用同一个字段名，比如一个用93Type，另一个时相则用00Type叠加后的文件在Arcmap中打开，选中文件，然后点右键——Property——空间查询，输入条件语句，比如：93Type=‘1’And00Type=‘2’；查询结果即为第一种类型转化为第二种类型的图形，可以另建

4、一图层比如：12，把查询结果复制到12图层上。统计出面积，依进行，就可以得到土地利用类型转移矩阵。