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

 

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

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

若你的数据是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)表示ij 两年份的土地利用变化图;NCi)表示i年份遥感分类影像;NCj)表示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=1And 00Type=2’;查询结果即为第一种类型转化为第二种类型的图形,可以另建一图层比如:12,把查询结果复制到12图层上。统计出面积,依进行,就可以得到土地利用类型转移矩阵

最后输出土地利用变化图 ,如下图所示:


集美大学
GISwww.jmugis.cn

Detail.lee

2008-5-30

转载于:https://www.cnblogs.com/Detail-lee/archive/2008/05/30/1210710.html

  • 0
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值