1、ArcGIS土地利用转移矩阵一、数据准备(图 1)准备两幅不同时相的土地利用现状图 (shp 格式) ,每幅图的属性表都要有一个表 示土地利用类型的字段, 并且要使用不同的名称加以区分, 如 Type1995 ,Type2000 土地利用类型名称必须统一,并且完整,如都使用“城镇用地”、“有林地”等。二、数据融合(图 2)在 ArcMap 里分别打开两个时相的图层,打开 ArcToolbox ,选择 Data Management Tools | Generalization | Dissolve 工具。 Input Feature 选择要融合的图层, Output Feature Class
2、 选择输出结果存储的位置及名称, Dissolve Field(s) 选择 土地利用类型字段(如 Type1995 ),然后勾选 Creat multipart features 选项,点击 OK 完成。重复此过程,对另一时相数据进行融合。此步骤使相同利用类型的记录融合 为一个记录,以提高后面步骤的计算速度。三、叠置分析(图 3)在 ArcMap 中打开两个时相融合后的数据, 在 ArcToolbox 中选择 Analysis Tools | Overlay | Intersect 工具, Input Features 选择两个时相的图层, Output Feature Class 选择叠加结
3、果存储的位置及名称,其余选项可以忽略,单击【 OK 】完成。四、计算面积并导出属性表(图 4-6 )在 ArcMap 中打开叠加后的图层数据, 在该图层上右键打开属性表, 选择 Option| Add field 新建一个字段,命名为 NewArea在 Editer 工具条中选择 Editer | Start Editing ,然后在属性表中 NewArea 字段上 单击右键选择 Calculate Geometry ,在打开的 Calculate Geometry 对话框中, Property 选择 Area ,Units 选择要使用的面积单位,单击【 OK 】完成图斑面积计算。 依次选择
4、Editer | Save Edits / End Editing 保存和退出编辑状态。在属性表中选择 Option | Export 将属性表保存为 dbf 文件五、制作转移矩阵(图 7-10 )(以 Excel2007 为例)在 Excel 中打开上一步保存的 dbf ,另存为 Excel 格式并打开。在 Excel 中选中所有数据(不要点左上角, 只选择有效数据),点击【插入】选项卡,选择【数据透视表】 |【数据透视表】,点击【确定】。在打开的数据透视表中按图示将字段拖入相应区域Excel 自动计算矩阵, 将该表稍事整饰就得到美观的土地利用转移矩阵。 矩阵中 r(I,j)就表示 i 类型
5、向 j 类型转移的土地面积,空值表示 i 类型向 j 类型没有转移。阅读全文 (131) / 评论 / 扔小纸条 / 给曾经善良过留言 收藏: QQ 书签 del.icio.us / 订阅: Google 抓虾分享】利用 ARCGIS 求算土地转移矩阵作者 : gaochangju (站内联系 TA) 收录: 2009-06-19 发布: 2009-06-19利用 ARCGIS 生成土地利用转移矩阵但我转移矩阵在土地利用变更调查时经常用到。 arcgis 里没有直接的工具来计算转移矩阵, 们可以利用已有的工具间接来求。下面介绍一下利用 arcgis 求转移矩阵的方法。 现有 2005 、 20
6、08 两期土地利用现状数据:在求转移矩阵前先做好下列工作:create1、备份数据。2 、分别将两期数据用 dissolve 工具按各自的地类编码合并要素, 注意选中最下面的multipat features ”选项处理好数据后,用intersect 工具求两期数据的交集交集的结果:将属性表增加面积字段,并计算面积结果:从上图中已经可以看出土地转移的情况: dlbm2005 是 2005 年时该块地的地类代码,dlbm2008 是 2008 年时该地块的地类代码,最后一列是相应地块的面积。接下来将属性表导出到 dbf 文件, 然后在 excel 中打开, 执行数据菜单中的 “数据透视表和 数据
7、透视图”命令点完成,然后将 dlbm2005 作为行字段, dlbm2008 作为列字段, area 作为数据拖到相应的区域生成的数据透视表若数据是 Vector 格式1Erdas Imagine Interpreter-Gis Analysis-Matrix ,输入两个时相的Vector 数据即可此时注意 输出栅格大小不应设的太小 要不一运行就会提示你的空间不足 做这一步之前,请做好前期的地理编码。2ArcView3.3 加载 spatial analysis 模块把两时相的 Vector 图转成 grid 格式(当然中间有一些单位的设置根据你做的图 的分辨率来设置即可) analysis-mapcaculate 直接计算即可。3把两期解译完的 Vector 文件在 arctoolbox overlay union 中叠加, 注意:两个文件不能用同一个字段名,比如一个用 93Type ,另一个时相则用 00Type叠加后的文件在 Arcmap 中打开,选中文件,然后点右键 Property 空 间查询,输入条件语句,比如: 93Type= 1And 00Type= 2 ;查询结果 即为第一种类型转化为第二种类型的图形,可以另建一图层比如: 12 ,把查询 结果复制到 12 图层上。统计出面积,依进行,就可以得到土地利用类型转移矩 阵。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1