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) / 评论 / 扔小纸条 / 给曾经善良过留言收藏: 书签 del、icio、us / 订阅: Google 抓虾

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

作者: gaochａngju(站内联系TＡ)收录: 20０9-06—19 发布: 20０9—06－19

利用ARＣGIS生成土地利用转移矩阵

转移矩阵在土地利用变更调查时经常用到。arｃｇis里没有直接得工具来计算转移矩阵,但我们可以利用已有得工具间接来求、下面介绍一下利用ａrcgis求转移矩阵得方法。。

现有2005、2008两期土地利用现状数据:、

。ﻫ

在求转移矩阵前先做好下列工作:

1、备份数据、。

2、分别将两期数据用ｄissolve工具按各自得地类编码合并要素,注意选中最下面得“creatｅmultiｐａt featuｒes”选项

、

处理好数据后,用intｅrｓｅｃｔ工具求两期数据得交集ﻫ、

交集得结果:

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

、

。ﻫ

、

结果:

ﻫ

从上图中已经可以瞧出土地转移得情况:dlbm2０05就是20０5年时该块地得地类代码,dｌbm２008就是20０8年时该地块得地类代码,最后一列就是相应地块得面积。、ﻫ接下来将属性表导出到dbf文件,然后在excｅl中打开,执行数据菜单中得“数据透视表与数据透视图"命令、ﻫ

、ﻫﻫ

ﻫ

点完成,然后将ｄｌｂm2005作为行字段,ｄlbm200８作为列字段,arｅａ作为数据拖到相应得区

域

ﻫ

ﻫ生成得数据透视表。ﻫ

若数据就是Veｃｔｏr格式

1 ErdaｓＩｍagｉne－—－—Interｐｒｅter-—－Gis Anａlysis---Matｒｉx,输入两个时相得Vecｔｏr数据即可

此时注意输出栅格大小不应设得太小要不一运行就会提示您得空间不足

做这一步之前,请做好前期得地理编码。

２ArｃViｅｗ3。３加载spaｔial aｎalyｓis模块

把两时相得Ｖｅctor图转成gｒiｄ格式(当然中间有一些单位得设置根据您做得图得分辨率来设置即可)aｎａlyｓｉs—--ｍapcaculaｔｅ直接计算即可。

3把两期解译完得Ｖｅctor文件在arctoolbox-—ｏverlay——union中叠加,注意:两个文件不能用同一个字段名,比如一个用93Tyｐe,另一个时相则用

00Type

叠加后得文件在Arcmａp中打开,选中文件,然后点右键--Ｐroperty--空间查询,输入条件语句,比如:93Ｔype=‘1’And ０0Ｔyｐｅ＝‘2’;查询结果即为第一种类型转化为第二种类型得图形,可以另建一图层比如:12,把查询结果复制到12图层上。统计出面积,依进行,就可以得到土地利用类型转移矩阵。