3. 地质图栅格数据加载与地理配准

扫描得到的地质图只是一张“图片”，它没有空间坐标，QGIS 无法判断图上的某个点对应地球表面的哪个位置。地理配准（Georeferencing）就是给这张图片建立像素坐标与真实地理坐标之间的映射关系，使其能够与其他空间数据叠加使用。本节围绕一幅 1:5 万地质图，从加载、控制点选取、坐标变换、残差评估到结果检查，完整走一遍配准流程。

扫描地质图的加载与预处理

配准前先确认扫描件质量。理想的地质图扫描件应保留完整图廓、四个图廓角点、坐标网（方里网）交点以及图名、比例尺、坐标系统说明等图面要素。若扫描时发生旋转、裁切或明显透视变形，应先在图像处理软件中校正或重新扫描，否则后续控制点误差会被放大。分辨率建议不低于 300 dpi，1:5 万图幅单张扫描后像素尺寸通常在数千像素量级，既能保证控制点拾取精度，又不至于让 QGIS 渲染过慢。

在 QGIS 中通过“图层 → 添加图层 → 添加栅格图层”加载扫描图，此时它只是普通图像，没有坐标参考。随后打开“栅格 → 配准（Georeferencer）”，把图像载入配准窗口。配准窗口左侧是待配准图像，右侧是控制点列表，底部显示变换参数与残差信息。需要特别注意：配准前应先在 QGIS 主界面设定好工程坐标参考系统（CRS），例如选择图幅所采用的 2000 国家大地坐标系 3 度带高斯-克吕格投影带，配准结果的目标坐标即以此为准。

配准原理与控制点选取

地理配准的本质是求解一个从图像像素坐标 (x, y) 到地图坐标 (X, Y) 的变换函数。控制点（GCP）就是已知两套坐标的对应点：在图像上点出某点的像素位置，再输入或拾取它的真实地图坐标。控制点数量与分布直接决定变换精度。多项式变换至少需要 (n+1)(n+2)/2 个控制点，一次多项式需 3 个，二次需 6 个，三次需 10 个；实际作业应多选 30%～50% 作为冗余，用于残差检验。

控制点应优先选在图廓角点、方里网交点、图面清晰且位置确定的三角点或明显地物上。避免选在等高线密集处、符号中心模糊处或图边被裁切的位置。控制点要均匀覆盖整幅图，四角与中部都要有，切忌全部集中在某一角，否则该区域外的变换会严重外推失真。对 1:5 万图幅，建议选取 12～20 个控制点，其中 4 个图廓角点作为骨架，其余分布在图内方里网交点。

坐标变换方法的选择

QGIS 配准器提供多种变换类型。对于扫描质量较好、变形以平移旋转和轻微缩放为主的地质图，一次多项式（仿射变换）通常足够，它保留直线为直线，不会产生弯曲。若扫描件存在纸张褶皱或镜头畸变导致的局部非线性变形，可选用二次或三次多项式，但阶数越高越容易在控制点稀疏区产生振荡，必须配合更多均匀分布的控制点。薄板样条（Thin Plate Spline）适合局部变形明显的情形，但外推能力差，图幅边缘需有控制点约束。工程实践中，1:5 万地质图优先尝试一次多项式，残差超限再考虑升阶。

重采样方法也需选择：地质图以线划和色块为主，建议使用“最近邻”以保持符号边界锐利，避免双线性或立方卷积造成线条发虚、颜色晕染。输出栅格时设定目标分辨率，1:5 万图幅地面分辨率 5 m 对应图上 0.1 mm，可按此换算输出像素大小，并指定与工程一致的 CRS。

残差评估与优化

配准器会为每个控制点计算残差，即该点实际坐标与变换后预测坐标之差，并给出均方根（RMS）总误差。判读要点有三：一看总 RMS 是否在可接受范围，1:5 万图幅一般要求控制在图上 0.2 mm 以内，换算到地面约 10 m，示例中 1.4 m 属于较理想情况；二看单个控制点残差是否显著偏大，若某点残差是其他点的数倍，多半是该点拾取错误或坐标输入有误，应检查后删除或重选；三看残差是否存在方向性规律，若所有点 dX 同号偏大，说明存在系统性平移，需检查坐标输入或 CRS 设置。

优化手段包括：删除粗差点、在残差大的区域补充控制点、调整变换阶数、核对控制点坐标的带号与单位。每次修改后重新计算，直至残差分布均匀且总误差达标。注意不要为了压低 RMS 而盲目删除正常点，残差均匀比数值极小更有意义。

配准结果检查

配准完成后，QGIS 会生成带坐标的栅格图层（或 world file）。检查分三步：第一，将配准图与已知矢量数据（如行政边界、道路、水系）叠加，观察图廓与地物是否套合；第二，缩放到图幅四角与中部，抽查方里网交点是否落在正确坐标位置；第三，用“识别要素”工具读取图上若干点的坐标，与图面注记比对。若发现整体偏移，回到配准器检查 CRS 与坐标单位；若局部扭曲，考虑升阶或补点。确认无误后保存配准结果，并保留控制点文件（.points）以便复现和追溯。

完成本节流程后，扫描地质图便从“无坐标图片”转变为可与地形图、矢量图层叠加的空间数据，为后续的图面要素矢量化与属性表整理奠定坐标基础。