利用ArcGIS实现SHP文件边界坐标批量导出为TXT格式
1. 为什么需要导出SHP边界坐标？
在实际工作中，我们经常会遇到需要将地理边界数据导出为文本格式的需求。比如城市规划部门需要将地块边界坐标提供给设计院做三维建模，环保机构需要将监测区域边界导入到数据分析软件中，或者测绘单位需要将测量结果以通用格式交付给客户。SHP文件作为GIS领域最常用的矢量数据格式，存储了丰富的空间信息，但有时我们需要将这些数据以更简单的文本形式呈现。
我遇到过这样一个案例：某物流公司需要将全国配送站点的位置信息导入到他们的路径规划系统中。由于该系统不支持直接读取SHP文件，我们就需要先将站点坐标导出为TXT格式。这个过程中，ArcGIS提供了完整的解决方案，能够高效准确地完成数据转换。
2. 准备工作：认识你的数据
2.1 理解SHP文件结构
SHP文件实际上是由多个文件组成的，包括.shp（主文件存储几何数据）、.shx（索引文件）和.dbf（属性表）。当我们谈论导出边界坐标时，主要处理的是.shp文件中的几何信息。对于面状要素（Polygon），其边界由一系列有序的点组成，这些点的坐标就是我们最终要导出的数据。
2.2 检查坐标系信息
在开始导出前，务必确认数据的坐标系。右键点击图层选择"属性"，在"源"选项卡中可以查看当前坐标系。如果是地理坐标系（如WGS84），坐标将以经纬度表示；如果是投影坐标系，则会显示为平面坐标（如米）。这个区别很重要，因为它决定了你最终导出的坐标值格式。
我曾经帮一个客户处理数据时，就因为没有注意坐标系，导致导出的坐标值单位错误，所有点位都偏移了几十万米。后来发现原始数据使用的是UTM投影，而我们误以为是经纬度坐标。
3. 完整操作流程详解
3.1 将面要素转换为点要素
由于SHP文件中的面要素是由闭合环线组成的，直接导出会得到完整的边界线。如果我们需要每个边界点的坐标，就需要先将面转换为点：
- 打开ArcToolbox（如果没看到，可以在菜单栏点击"Geoprocessing">"ArcToolbox"）
- 导航至"Data Management Tools">"Features"
- 双击"Feature Vertices to Points"工具
- 在弹出窗口中：
- 输入要素选择你的面图层
- 输出要素类指定保存位置和名称
- 点类型选择"ALL"（导出所有顶点）
这个步骤会创建一个新的点图层，其中每个点对应原始面要素的一个顶点。我建议在处理前先备份原始数据，因为转换过程是不可逆的。
3.2 添加XY坐标字段
有了点图层后，我们需要为每个点计算坐标值：
- 在ArcToolbox中导航至"Data Management Tools">"Features"
- 双击"Add XY Coordinates"工具
- 选择上一步生成的点图层作为输入要素
- 点击确定运行工具
运行完成后，打开属性表会看到新增的POINT_X和POINT_Y字段（如果是地理坐标系则是POINT_X和POINT_Y）。这些字段存储了每个点的坐标值。需要注意的是，这个工具会直接修改原始数据，所以最好在副本上操作。
3.3 导出属性表为文本文件
现在我们已经有了包含坐标值的属性表，接下来就是导出为TXT：
- 右键点击图层选择"Open Attribute Table"
- 在属性表窗口中点击"Table Options"按钮（右上角）
- 选择"Export"
- 在保存类型中选择"Text File(*.txt)"
- 指定保存路径和文件名
- 在"Field"选项卡中可以选择要导出的字段（建议只保留必要的坐标字段）
- 点击"OK"完成导出
导出的TXT文件默认是逗号分隔的，可以用Excel或其他文本编辑器打开。如果需要其他分隔符，可以在导出后使用文本编辑器的替换功能进行转换。
4. 高级技巧与常见问题
4.1 批量处理多个SHP文件
如果需要处理大量SHP文件，手动操作会很耗时。这时可以使用Model Builder创建自动化流程：
- 打开Model Builder（在菜单栏点击"Geoprocessing">"Model Builder"）
- 将"Feature Vertices to Points"和"Add XY Coordinates"工具拖入模型
- 连接两个工具（第一个工具输出作为第二个工具的输入）
- 右键点击第一个工具选择"Make Variable">"From Parameter">"Input Features"
- 将这个变量设为模型参数（右键变量选择"Model Parameter"）
- 保存模型后，就可以批量输入多个SHP文件了
我曾在一次项目中用这个方法处理了300多个县域边界文件，原本需要几天的工作量，现在只需要运行一次模型就能完成。
4.2 处理带洞的多边形
有些面要素可能包含"洞"（比如湖泊中的岛屿），这时直接转换会得到内外边界的点。如果只需要外边界，可以在"Feature Vertices to Points"工具中选择"DANGLE"选项，这样只会保留多边形外环的点。
4.3 坐标精度控制
默认情况下，ArcGIS会输出完整的双精度坐标值。如果不需要这么高的精度，可以在导出前使用"Calculate Field"工具对坐标值进行四舍五入：
- 在属性表中右键点击POINT_X字段
- 选择"Calculate Field"
- 在表达式框中输入"Round([POINT_X],4)"（保留4位小数）
- 对POINT_Y字段重复相同操作
这个技巧在减小文件体积方面特别有用，我曾经将一个50MB的TXT文件通过控制精度缩减到了15MB，而精度损失对项目影响微乎其微。
5. 实际应用案例
去年我参与了一个城市绿地规划项目，需要将全市200多个公园的边界坐标导出给3D建模团队。原始数据是一个包含所有公园的SHP文件，每个公园作为一个独立的面要素。通过以下步骤高效完成了任务：
- 首先使用"Multipart to Singlepart"工具将多部分要素拆分为单个公园
- 然后按上述方法将每个公园面转换为点并添加XY坐标
- 使用"Split by Attributes"工具按公园名称分别导出
- 最后用Python脚本批量重命名所有输出文件
整个过程从数据准备到最终交付只用了半天时间，客户对效率和准确性都非常满意。这个案例展示了ArcGIS在处理复杂空间数据转换时的强大能力。
6. 替代方案与工具比较
虽然ArcGIS是处理SHP文件的主流工具，但也有一些替代方案值得了解。QGIS作为开源GIS软件，同样提供类似功能。在QGIS中可以通过"Vector">"Geometry Tools">"Extract Vertices"完成点转换，然后使用"Vector">"Research Tools">"Export/Add Geometry Columns"添加坐标值。
相比之下，ArcGIS的优点是流程更标准化，工具集成度更高，特别适合企业级批量处理。而QGIS的优势在于免费开源，对于预算有限的个人或小团队是不错的选择。根据我的经验，ArcGIS在处理大型数据集时稳定性更好，特别是在执行批量操作时不易崩溃。