4. 坐标系转换原理与参数选择

在 CGCS2000 转换排错中，绝大多数“转完偏几十米”“高程整体抬升”“局部对得上、远处错得离谱”的问题，根源都不在软件操作，而在模型选错或参数用错。本节把两套最常用的转换模型讲清楚：它们各自描述什么几何关系、需要几个参数、参数从哪来、什么条件下能用、精度怎么评估。

一、四参数与七参数：两种不同的几何假设

四参数模型（二维相似变换）只处理平面关系，把源平面坐标 (x, y) 变换到目标平面坐标 (X, Y)，包含 2 个平移量 ΔX、ΔY，1 个旋转角 θ，1 个尺度因子 k（常写成 m = 1 + k）。其数学形式为：

X = ΔX + (1+k)·(x·cosθ − y·sinθ)
Y = ΔY + (1+k)·(x·sinθ + y·cosθ)

四参数不改变高程，适用于两个平面坐标系之间只存在平移、旋转和均匀缩放的情形，典型场景是同一椭球下的地方独立坐标系与 CGCS2000 投影坐标之间的互转，或城市独立坐标与 2000 国家坐标系的衔接。

七参数模型（三维布尔莎模型）处理三维空间直角坐标 (X, Y, Z)，包含 3 个平移量 (ΔX, ΔY, ΔZ)、3 个旋转量 (εx, εy, εz) 和 1 个尺度因子 m。其小角度近似形式为：

X₂ = ΔX + (1+m)·X₁ + εz·Y₁ − εy·Z₁
Y₂ = ΔY − εz·X₁ + (1+m)·Y₁ + εx·Z₁
Z₂ = ΔZ + εy·X₁ − εx·Y₁ + (1+m)·Z₁

七参数描述的是两个不同椭球（基准）之间的空间关系，例如北京54、西安80 与 CGCS2000 之间，因为椭球长半轴、扁率、地心位置都不同，必须用三维模型才能同时吸收平面与高程的系统性差异。注意旋转角的符号约定：我国常用“坐标轴旋转、逆时针为正”的约定，而部分软件（如某些国外工具）采用相反符号，混用时会出现“旋转量反号”的典型错误。

二、参数从哪里来：三条获取途径

途径一：国家或行业公布的区域参数。测绘主管部门会发布大区域七参数或分带四参数，例如某些省份的 CGCS2000 与西安80 转换参数。这类参数覆盖范围大、权威性高，但属于“区域平均”，在边缘地带可能残差偏大，使用前应确认参数文件的适用区域、中央子午线和椭球定义。

途径二：用公共点反算。这是工程中最常见的做法。收集至少 3 个（四参数）或 4 个（七参数）同时具备两套坐标的控制点，用最小二乘求解参数。四参数至少 2 个点即可解算、3 个点可检核；七参数至少 3 个点、建议 4 个以上并留 1 个检核点。公共点应分布均匀、覆盖整个测区，避免全部挤在一角——这是“局部准、远处偏”的头号原因。

途径三：从既有成果或软件参数库导入。如已有历史项目的转换参数文件（.txt、.trd 等），可直接导入，但必须核对参数单位（旋转量是弧度还是秒）、符号约定和椭球参数，否则会出现量级错误。

三、参数适用条件与精度评估

选模型前先问三个问题：是否跨椭球？跨椭球必须用七参数（或先做基准转换再做投影）；是否只关心平面？若高程另有来源且平面为同一基准，四参数更稳健；测区范围多大？四参数在几十公里内精度较好，大范围应分区或改用七参数。

精度评估不能只看解算时的残差。正确做法是预留检核点：用未参与解算的控制点代入参数，计算 ΔX、ΔY、ΔH 的差值，统计中误差。工程上常以检核点平面残差 ≤ 5 cm、高程残差 ≤ 10 cm 作为可用参考线（具体限差应按项目规范确定）。若残差呈系统性方向偏移，说明参数区域不匹配；若残差随距离增大而增大，说明尺度或旋转有问题；若个别点残差极大，应排查该点本身的坐标错误或点位变动。

最后强调一点：参数不是“一次解算、永久通用”的。控制点更新、测区扩展、中央子午线变更，都应重新评估参数适用性。把参数来源、公共点清单、检核残差记录成文档，是后续排错时最有力的依据。