4. 轨迹坐标系统与投影转换

野外GPS轨迹整理中,坐标问题是最容易出错、也最容易被忽视的环节。一条轨迹在QGIS里看起来“位置不对”,十有八九不是采集设备的问题,而是坐标系定义或投影转换没有处理正确。本节围绕地理坐标系与投影坐标系的区别、野外常用坐标系的适用场景,以及QGIS中轨迹投影变换的实际操作展开。

地理坐标系与投影坐标系:两种描述位置的方式

地理坐标系(Geographic Coordinate System,GCS)用经纬度表示位置,单位是度,本质是把地球近似为旋转椭球体后,用角度描述地表点。它适合存储和交换数据,因为不依赖具体投影,全球通用。但经纬度不能直接用于量算距离和面积——一度经度在赤道约111 km,在纬度40°处只有约85 km,直接用度数做长度统计会得到错误结果。

投影坐标系(Projected Coordinate System,PCS)则是把椭球面按某种数学规则展开到平面上,单位是米。它由“地理坐标系 + 投影方法 + 投影参数”三部分共同定义。只有投影坐标系下的轨迹才能直接计算长度、面积、缓冲距离。因此野外轨迹整理的通行做法是:采集与存储用地理坐标系,量算与分析用投影坐标系。

需要特别强调一个概念:坐标转换(不同基准面之间,如WGS84到CGCS2000)与投影变换(同一基准面下经纬度到平面米制)是两件事。前者涉及基准转换参数,可能带来米级甚至更大的位移;后者只是数学展开,不改变点的实际地理位置。混淆二者是轨迹错位的常见根源。

野外常用坐标系及其选择

野外地质调查中常见的坐标系主要有三类:

  • WGS84:GPS原始输出的坐标系,手持GPS、手机记录、无人机航迹默认多为WGS84。它是全球通用的地心坐标系,适合作为轨迹的“原始存档”坐标系。
  • CGCS2000:我国现行国家大地坐标系,与WGS84在厘米级上基本一致(差异通常在10 cm以内),在一般地质填图精度下可近似等同处理,但正式成果提交时应明确标注为CGCS2000。
  • 地方坐标系:如北京54、西安80,或各省市的独立坐标系(如“××市城建坐标系”),常带有特定的中央子午线和投影面高程。这类坐标系与WGS84之间存在明确的基准转换参数,必须使用测区提供的转换参数(三参数或七参数)进行转换,不能随意套用。
判断原则:如果轨迹只用于本队内部查看和粗略量算,统一到WGS84或CGCS2000的地理坐标系即可;如果要与已有地形图、地质图、征地红线叠加,则必须转换到该图件所用的投影坐标系,并确认基准面一致。

QGIS中的坐标转换工具

QGIS的坐标处理能力集中在几个位置。图层属性中的“源”选项卡用于定义数据本身的坐标系(CRS),这是“告诉QGIS数据是什么”,不是转换。项目属性中的“CRS”选项卡决定显示时使用哪个坐标系,QGIS会自动做实时投影变换。真正把数据另存为新坐标系,则通过右键图层 → 导出 → 要素另存为,在“CRS”下拉框中选择目标坐标系。

对于需要基准转换的情况,QGIS在转换时会弹出或使用默认的转换路径。可在“设置 → 选项 → 坐标转换”中查看和指定转换参数,例如WGS84到西安80的七参数。若测区提供了地方转换参数,应在此处或通过“矢量 → 数据管理工具 → 重新投影图层”指定,避免使用默认的近似转换。

推荐操作顺序(以一条WGS84轨迹为例):
1. 导入轨迹 → 图层属性确认源CRS为 EPSG:4326(WGS84)
2. 项目CRS设为 EPSG:4547(CGCS2000 / 3度带 114°E)
3. 右键图层 → 导出 → 要素另存为 → CRS选 EPSG:4547
4. 对新图层用“字段计算器”或“导出添加几何属性”计算长度
5. 与地形图叠加检查套合情况

轨迹投影变换实战

假设一条野外轨迹以GPX格式导入,包含约1200个点,原始坐标系为WGS84。目标是转换到CGCS2000的3度带投影(中央子午线114°E,EPSG:4547),用于计算路线总长并与1:1万地形图套合。

第一步,导入后立即检查图层CRS是否为EPSG:4326。若GPX导入后显示为未知坐标系,需手动指定,否则后续转换全部失效。第二步,将项目CRS设为EPSG:4547,此时轨迹在画布上会“跳”到正确位置,但这只是显示变换,数据本身仍是经纬度。第三步,执行“要素另存为”,目标CRS选EPSG:4547,得到一份米制坐标的新图层。第四步,打开新图层的属性表,用“导出 → 添加几何属性”生成length字段,单位为米,即可统计总长。第五步,叠加测区地形图,观察轨迹是否落在正确的沟谷、山脊位置上;若整体偏移一个固定量,多半是基准转换参数缺失或选错。

轨迹从采集到量算的坐标处理链 野外采集 GPS / 手机 / 无人机 输出 WGS84 经纬度 定义源坐标系 EPSG:4326 确认而非转换 投影变换 EPSG:4547 CGCS2000 3度带 米制量算 长度 / 面积 与地形图套合 基准转换(另一条路径) WGS84 → 西安80 / 北京54 / 地方坐标 需三参数或七参数,不可省略 地方投影坐标 与既有图件叠加 判断要点:整体固定偏移 → 查基准转换参数;长度明显偏大或偏小 → 查是否误用经纬度直接量算; 与图件套合不上但形状正确 → 查投影带号与中央子午线是否选错。

实际作业中还有两个细节值得注意。其一,投影带号选择错误会导致轨迹整体平移数百公里,表现为“形状对、位置错”,此时应核对测区经度范围与3度带、6度带的对应关系。其二,若轨迹需要与扫描地质图配准,图件本身的坐标系往往需要先确认,再决定轨迹转换到哪个投影,而不是反过来。

不要用“设置项目CRS”代替“导出重投影”。项目CRS只影响显示,导出的数据仍是原坐标系;若直接把显示结果截图或量算,得到的数值可能基于临时变换,无法复现。正式成果必须落到一份明确标注CRS的矢量文件上。

完成投影变换后,建议保留两份数据:一份WGS84原始轨迹作为存档,一份目标投影坐标系的轨迹用于分析出图,并在文件名或属性中写明坐标系与EPSG代码。这样在后续汇交、复核或与其他人协作时,坐标来源清晰可查,避免反复返工。