4. 地层数据获取与预处理

在ArcGIS Pro中做地层与断层的空间叠加分析，第一步并不是打开叠加工具，而是把“地层”这件事变成一份可被空间运算的数据。地层信息在现实中以三种典型形态存在：公开地质图（面状地质单元）、钻孔数据（点状分层记录）、剖面图（线状或栅格化的地层界线）。它们各自携带不同的空间语义，必须先经过获取、导入、格式转换与属性整理，才能进入统一的地层框架。本节围绕这三类来源，讲清楚每一步实际要做什么、参数如何设定、结果如何判断。

一、三类来源的获取与空间语义

公开地质图是最常用的底图来源，例如1:20万或1:5万区域地质图，通常以Shapefile、GeoPackage或PDF形式发布。若拿到的是矢量格式，重点检查其属性表是否包含“地层代号”“岩性”“时代”等字段；若只有PDF或扫描图，则需要先做地理配准（Georeferencing），用至少4个已知控制点把图像纠正到目标坐标系，再通过屏幕数字化把地层界线描成面要素。钻孔数据一般以Excel或CSV表格提供，每行是一个“钻孔—层位”记录，字段包括孔号、X坐标、Y坐标、层顶深度、层底深度、地层代号；导入时用“XY表转点”生成点要素，再按孔号聚合。剖面图则常以图片或CAD文件存在，需要沿剖面线建立线要素，并把地层界线转成带高程属性的线或点。

三类来源的空间语义差异决定了后续处理方式：地质图是面，可直接参与叠加；钻孔是点，需要先做空间插值或作为约束点；剖面是线，需要投影到平面或转成剖面栅格。理解这一点，才能避免把点、线、面混在一起直接做相交运算。

二、导入与格式转换的关键操作

导入阶段最容易出问题的是坐标系。假设目标工程使用CGCS2000 3度带高斯投影（中央经线117°E），而地质图是WGS84经纬度，钻孔表里的X/Y是北京54坐标。此时不能直接“投影”了事，而应先用“投影”工具把矢量数据转到目标坐标系，对钻孔表则用“XY表转点”时指定其原始坐标系，再投影。若原始坐标系未知，可用“定义投影”先赋一个假定坐标系，但必须在属性中记录这一假设，避免后续叠加出现系统性偏移。

格式转换方面，Shapefile字段名限制10个字符，中文属性容易截断，建议转成文件地理数据库（File Geodatabase）要素类，字段名可保留完整中文。PDF地质图经配准后，用“创建要素类”建面图层，编辑时开启捕捉，沿地层界线逐段数字化，闭合后写入地层代号。钻孔Excel导入前，先检查层顶深度是否从0起算、层底是否大于层顶、是否存在重叠层；用字段计算器生成“层顶高程 = 孔口高程 − 层顶深度”，同理得层底高程，这是后续三维或厚度计算的基础。

三、属性整理与质量检查

属性整理的目标是让三类数据说同一种“地层语言”。核心是建立统一的地层代号字典：例如把地质图的“J3”与钻孔的“上侏罗统”映射到同一代号“J3”，并补充岩性、时代两个辅助字段。可用“添加字段”+“计算字段”批量赋值，也可用“连接”把外部对照表挂到要素类上。厚度字段由层底高程减层顶高程得到，出现负值说明层序录入颠倒，需回查原始记录。

质量检查建议至少做三项：一是坐标系一致性检查，用“要素类属性”核对每个图层的空间参考是否与工程一致；二是拓扑检查，对面图层建“不能有重叠”“不能有空隙”规则，对线图层检查悬挂点；三是属性完整性检查，用“按属性选择”筛出地层代号为空或厚度为负的记录，逐条修正。完成这些后，数据才具备进入叠加分析的条件。

把这一步做扎实，后续的地层与断层叠加才有可靠的空间与属性基础；反之，坐标系错位或代号不统一，会让相交、擦除等运算产生大量伪结果，且难以追溯。