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书签 / 订阅: Google 抓虾