实验二：面状要素中轴线提取

一、实验目的

掌握栅格数据和矢量数据的几何量测参数计算方法，理解地理空间数据的几何参数含义，会初步应用空间量测发散性解决实际问题。

二、实验数据

国家1:10,000面状水系数据。

三、实验步骤

3.1投影转换

依次打开Data Management Tools->Projectionsand Transformations->Feature->Project选择WGS_1984_World_Mercator投影，参数设置如下图所示：

 

分析：地图投影转换是地图投影和 地图编绘的一个重要组成部分。当GIS使用的数据取自不同地图投影的图幅时就需要将一种投影数字化数据转换为所需要的投影坐标数据。

3.2采用几何量测的方法计算面状水系要素的中轴线

1、首先消除面状水系中的岛，做法是先生成外包矩形。依次打开Data Management Tools->Features->Feature Envelope to Polygon

 

◆Input raster or featuresource data：输入要进行距离分析的源数据，在这里输入投影转换好的面状水系，点击ok,如下图所示：

◆选择AnalysisTools->Overly->Symmetrical Difference, 如图所示设置相关参数：

 

◆单击ok 得到Difference数据：

◆ 打开属性表会发现属性值唯一：

◆现在把属性打散；依次打开ArcToolbox，选择Data Management Tools->Features->Multipart To Singlepart，打散后的如下图所示：

◆打开ArcToolbox，选择Data Management Tools->Features->Delete Features，在弹出的对话框中输入“打散”图层。完成数据的删除：

◆打开ArcToolbox，选择选择Data Management Tools->General->Merge，设置相关参数，点击ok，完成图层的合并。如下图所示：

◆打开ArcToolbox，选择Data Management Tools->Generalization->Dissolve，如图所示：

◆选择Data ManagementTools->Features->Feature To Line，如图所示：

◆打开Editor Toolbar工具，选择Start Editing，编辑刚得到的线图层，选择Split Tool工具，在合适的位置，将数据打断。如图所示：

◆选择要合并的线段，然后Editor->Merge,在弹出的对话框中点击ok，完成图层的合并。如图所示：

◆右击featuretoline，打开属性表，右击字段id，选择Field Calculate，在弹出的对话框中点击FID（使字段id=FID）。如下图所示：

◆选择SpatialAnalyst Tools->Distance->Euclidean Allocation，设置相关参数如图所示：

◆点击ok，完成欧式分配：

◆选择ConversionTools->Raster To Polygon,如图所示：

◆选择Data ManagementTools->Features->Polygon To Line，右击图层“rastert_eucallo1_PolygonToLi”打开属性表，选择需要删除的线段；打开ArcToolbox，选择Data ManagementTools->Features->Delete Features，在弹出的对话框中输入“rastert_eucallo1_PolygonToLi”图层。如图所示：

◆至此中轴线提取完毕

四、结论与讨论

采用几何量测的方法计算面状水系要素的中轴线十分有用，在流域分析中是一种常见的方法。这也是泰森多边形典型使用。但是步骤比较繁琐，建议多加练习。

转载自：https://blog.csdn.net/vector5210/article/details/46461519