一、数据准备(图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两期土地利用现状数据:
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_0b6e101fc13f480892f2s8hgVLYlkXXG.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_90f0537fe91dfda3df93nkE9fpooJQlw.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_08e60c477f4e3c738916wAuCofAWVnAe.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_14792a3fd38e04114770eJamOU0rvlXW.jpg" \* MERGEFORMATINET
在求转移矩阵前先做好下列工作:
1、备份数据。
2、分别将两期数据用dissolve工具按各自的地类编码合并要素,注意选中最下面的“create multipat features”选项
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_aa9028c2682005d6fd0cttBRvAl8Vgtd.jpg" \* MERGEFORMATINET
处理好数据后,用intersect工具求两期数据的交集
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_130a02cb7060768b0856uPbWHLLNCdjx.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_32dadec0eb1bcc88a47bexYGEIFJGwSX.jpg" \* MERGEFORMATINET
交集的结果:
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_2a1bc4bab7f92f918cae1WkEpeqNCcjN.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_1254d5db0a9d03d76451SNmcGLcVLFRX.jpg" \* MERGEFORMATINET
将属性表增加面积字段,并计算面积
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_5af83e411359435d77cbIECo4quIaoWp.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_a158b2edfdab43ad1cd7tkGKonz5nN0m.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_dc5f14fe6bca43d86199q1YqLAoNb58I.jpg" \* MERGEFORMATINET
结果:
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_bc597ac4fb05211eb097ybdDS9CsZLXr.jpg" \* MERGEFORMATINET
从上图中已经可以看出土地转移的情况:dlbm2005是2005年时该块地的地类代码,dlbm2008是2008年时该地块的地类代码,最后一列是相应地块的面积。
接下来将属性表导出到dbf文件,然后在excel中打开,执行数据菜单中的“数据透视表和数据透视图”命令
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_5a991f955dee4bab97424MCrAtidhswb.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_434aa9c354132c730145DOki1t4W5AfE.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_955bff9400ded3c7fbaeV6UUte1hPXrL.jpg" \* MERGEFORMATINET
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_eb2081f69554d4d0e86aZA5oH845n6iJ.jpg" \* MERGEFORMATINET
点完成,然后将dlbm2005作为行字段,dlbm2008作为列字段,area作为数据拖到相应的区域
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_8eff1f70926e8166312bsggq8V6PnDUa.jpg" \* MERGEFORMATINET
生成的数据透视表
INCLUDEPICTURE "http://bbs.esrichina-bj.cn/ESRI/attachments/forumid_24/20090201_02565947d592ab8e405fNHlagAdRrEi6.jpg" \* MERGEFORMATINET
若数据是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图层上。统计出面积,依进行,就可以得到土地利用类型转移矩阵 。