ArcGIS土地利用转移矩阵.doc

(word完整版)ArcGIS土地利用转移矩阵
一、数据准备(图1）
　　准备两幅不同时相的土地利用现状图（shp格式)，每幅图的属性表都要有一个表示土地利用类型的字段，并且要使用不同的名称加以区分，如Type1995,Type2000.土地利用类型名称必须统一，并且完整，如都使用“城镇用地”、“有林地”等。

二、数据融合（图2）
　　在ArcMap里分别打开两个时相的图层，打开ArcToolbox，选择Data Management Tools | Generalization | Dissolve工具。Input Feature选择要融合的图层，Output Feature Class选择输出结果存储的位置及名称，Dissolve Field(s）选择土地利用类型字段(如Type1995），然后勾选Creat multipart features选项，点击OK完成。重复此过程,对另一时相数据进行融合。此步骤使相同利用类型的记录融合为一个记录,以提高后面步骤的计算速度。

三、叠置分析（图3)
　　在ArcMap中打开两个时相融合后的数据，在ArcToolbox中选择Analysis Tools | Overlay | Intersect工具，Input Features选择两个时相的图层，Output Feature Class选择叠加结果存储的位置及名称，其余选项可以忽略，单击【OK】完成。

四、计算面积并导出属性表（图4—6）
　　在ArcMap中打开叠加后的图层数据，在该图层上右键打开属性表,选择Option ｜ Add field… 新建一个字段，命名为NewArea。

　　在Editer工具条中选择Editer | Start Editing，然后在属性表中NewArea字段上单击右键选择Calculate Geometry… ，在打开的Calculate Geometry对话框中，Property选择Area，Units选择要使用的面积单位,单击【OK】完成图斑面积计算.依次选择Editer | Save Edits / End Editing保存和退出编辑状态.

　　在属性表中选择Option ｜ Export… 将属性表保存为dbf文件。

五、制作转移矩阵（图7—10）(以Excel2007为例)
　　在Excel中打开上一步保存的dbf,另存为Excel格式并打开。在Excel中选中所有数据（不要点左上角，只选择有效数据），点击【插入】选项卡，选择【数据透视表】｜【数据透视表】,点击【确定】.

　　在打开的数据透视表中按图示将字段拖入相应区域。

　　Excel自动计算矩阵，将该表稍事整饰就得到美观的土地利用转移矩阵.矩阵中r(I， j）就表示i类型向j类型转移的土地面积，空值表示i类型向j类型没有转移。

阅读全文(131) / 评论 / 扔小纸条 / 给曾经善良过留言
收藏： QQ书签 del.icio.us / 订阅： Google 抓虾

【分享】利用ARCGIS求算土地转移矩阵

作者： gaochangju （站内联系TA）    收录： 2009—06-19    发布： 2009—06-19

-
利用ARCGIS生成土地利用转移矩阵
转移矩阵在土地利用变更调查时经常用到.arcgis里没有直接的工具来计算转移矩阵，但我们可以利用已有的工具间接来求。下面介绍一下利用arcgis求转移矩阵的方法。
现有2005、2008两期土地利用现状数据：

在求转移矩阵前先做好下列工作:
1、备份数据。
2、分别将两期数据用dissolve工具按各自的地类编码合并要素，注意选中最下面的“create multipat features”选项

处理好数据后,用intersect工具求两期数据的交集

交集的结果：

将属性表增加面积字段,并计算面积

结果：

从上图中已经可以看出土地转移的情况：dlbm2005是2005年时该块地的地类代码，dlbm2008是2008年时该地块的地类代码，最后一列是相应地块的面积。
接下来将属性表导出到dbf文件,然后在excel中打开，执行数据菜单中的“数据透视表和数据透视图"命令

点完成，然后将dlbm2005作为行字段，dlbm2008作为列字段，area作为数据拖到相应的区域

生成的数据透视表

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

欢迎您的光临，Word文档下载后可修改编辑.双击可删除页眉页脚.谢谢！希望您提出您宝贵的意见,你的意见是我进步的动力。赠语; 1、如果我们做与不做都会有人笑，如果做不好与做得好还会有人笑，那么我们索性就做得更好，来给人笑吧！ 2、现在你不玩命的学，以后命玩你。3、我不知道年少轻狂,我只知道胜者为王。4、不要做金钱、权利的奴隶；应学会做“金钱、权利”的主人。5、什么时候离光明最近？那就是你觉得黑暗太黑的时候。6、最值得欣赏的风景,是自己奋斗的足迹。 7、压力不是有人比你努力,而是那些比你牛×几倍的人依然比你努力.

word版本