注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

快乐小屋

 
 
 

日志

 
 

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

2012-03-08 20:28:52|  分类: ARCGIS |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

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

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

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

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

集美大学GIS(www.jmugis.cn)

  评论这张
 
阅读(219)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018