4. 矿产资源空间数据坐标系与投影

在ArcGIS中做矿产空间插值与等值线绘制，第一步不是选插值模型，而是把“位置”这件事说清楚。矿点坐标、矿区边界、地质图、DEM、遥感影像往往来自不同部门、不同年代、不同软件，坐标系和投影各不相同。如果空间参考不统一，即使插值算法完全正确，生成的等值线也会整体偏移、拉伸或与地质体错位，后续圈定靶区、估算资源量都会失去意义。本节围绕坐标系选择、投影变换、数据对齐与空间参考统一四个环节展开，目标是让所有参与插值的图层落在同一个可度量的平面坐标框架内。

一、先分清地理坐标系与投影坐标系

地理坐标系（GCS）用经纬度表示位置，单位是度，基于椭球体和基准面，例如 WGS 1984、CGCS2000。它适合全球定位与数据交换，但不适合直接量算面积和距离——因为一度经度在赤道约111 km，在纬度60°处只有约55 km。投影坐标系（PCS）则是把椭球面按某种投影方法展开到平面上，单位为米，例如 UTM 分带坐标、高斯-克吕格3度带或6度带坐标、Albers 等积投影。插值生成等值线、计算矿体面积、统计品位变异函数，都应在投影坐标系下进行。

我国矿产资源数据常见两套基准：一是 CGCS2000 国家大地坐标系，二是历史数据中的北京54、西安80。CGCS2000 与 WGS 1984 在一般矿产制图精度下差异很小，但北京54、西安80 与 CGCS2000 之间存在几十米到上百米的基准差异，必须通过七参数或三参数转换，不能简单当作同一坐标系处理。判断数据基准，优先查看数据附带的 .prj 文件、元数据或采集说明，而不是凭经验猜测。

二、投影选择：按用途和范围定参数

投影选择的核心是“保什么”。资源量估算、面积统计优先选等积投影，如 Albers 等积圆锥投影；勘查线布置、距离量算优先选等距或等角投影，如高斯-克吕格、UTM；全国或大区级编图常用 Albers，中央经线取制图区中央，双标准纬线取范围南北约1/6和5/6处。对于省级或矿区级数据，高斯-克吕格3度带更贴合我国成图习惯，需明确带号与中央经线，例如某矿区位于东经113°附近，3度带带号为38，中央经线114°。

需要强调的是：投影参数必须与数据范围匹配。把跨两个3度带的矿带强行投到单一分带，边缘变形会明显增大；此时应改用自定义的局部横轴墨卡托或Albers，并在元数据中记录中央经线、原点纬度、比例因子等参数，保证他人可复现。

三、投影变换：Project 与 Define Projection 不能混用

ArcGIS 中两个工具极易混淆。Define Projection（定义投影）只写入或修改空间参考描述，不改变坐标数值，适用于“数据本身坐标正确但缺少 .prj”的情况；Project（投影）则按数学变换真正重算坐标，适用于坐标系不一致需要转换的情况。若对一批经纬度数据误用 Define Projection 强行标成投影坐标，数值仍是度却被当作米，图层会瞬间“缩”到原点附近，这是最常见的错误之一。

矢量数据用 Project 工具，需指定输出坐标系与地理变换（Geographic Transformation）。当源为西安80、目标为CGCS2000时，应选择对应的七参数转换；若研究区较小且无精确参数，可用三参数并记录残差。栅格数据用 Project Raster，重采样方法要按数据性质选：连续型 DEM、品位插值面用双线性（BILINEAR）或三次卷积（CUBIC），分类图、土地利用、地质界线用最近邻（NEAREST），避免产生原数据中不存在的新类别值。

四、数据对齐与空间参考统一

统一空间参考后，还要保证栅格之间“格网对齐”。插值输出、DEM、遥感影像若像元大小、原点、行列数不一致，叠加分析会出现半像元错位。可在环境设置中固定 Snap Raster（捕捉栅格）、Cell Size（像元大小）和 Output Extent（输出范围），让所有栅格共用同一格网。矢量与栅格叠加前，用“按掩膜提取”或“裁剪”把范围收敛到矿区边界外扩一定距离的缓冲区，既减少计算量，也避免边界外插值外推失真。

对齐效果要通过具体检查判断，而不是“看起来差不多”。常用做法：叠加矿点与地质图，放大到1:2000查看点位是否落在已知采坑或露头中心；选取2～3个有实测CGCS2000坐标的控制点，与转换后坐标比较，平面残差一般应小于图上0.1 mm对应的地面距离；检查栅格属性中的像元大小、左上角坐标、行列数是否与基准栅格完全一致。若残差偏大，优先排查基准转换参数和中央经线设置，而不是反复调整插值参数。

完成上述步骤后，矿点、边界、地质图、DEM 才真正处于同一平面坐标框架，后续的变异函数拟合、克里金插值和等值线追踪才有可靠的位置基础。空间参考统一不是一次性的格式转换，而是贯穿数据准备全过程的约束条件。