3. 高斯投影正算公式推导与解析

高斯投影正算解决的是这样一个问题：已知某点在参考椭球面上的大地坐标（大地纬度 B、大地经度 L），求它在高斯投影平面上的直角坐标（x, y）。在矿山测量中，这一计算是把国家控制点、GPS 实测点纳入矿区平面坐标系统的必经环节，也是后续换带计算的基础。理解公式结构，比记住公式本身更重要——只有弄清每一项从哪来、量级多大、精度受什么影响，才能在实际计算中判断结果是否可信。

一、公式的整体结构：从椭球面到平面

高斯投影是等角横切椭圆柱投影，其正算公式以经差 l = L − L₀ 为基本变量（L₀ 为中央子午线经度），按 l 的幂级数展开。我国现行规范采用的形式为：

x = X + N/2·sinB·cosB·l² + N/24·sinB·cos³B·(5 − t² + 9η² + 4η⁴)·l⁴
      + N/720·sinB·cos⁵B·(61 − 58t² + t⁴)·l⁶

y = N·cosB·l + N/6·cos³B·(1 − t² + η²)·l³
      + N/120·cos⁵B·(5 − 18t² + t⁴ + 14η² − 58η²t²)·l⁵

式中：X 为赤道至纬度 B 的子午线弧长；N 为卯酉圈曲率半径；t = tanB；η² = e′²cos²B（e′ 为第二偏心率）；l 以弧度为单位。可以看到，x 的展开只含 l 的偶次项，y 只含奇次项，这是投影关于中央子午线对称性的直接体现：同纬度下经差反号，x 不变、y 反号。

公式的物理含义可以这样理解：x 的主项 X 是“沿中央子午线量到该纬度的弧长”，后面的 l²、l⁴ 项是对离开中央子午线后子午线收敛、长度变化的改正；y 的主项 N·cosB·l 相当于把纬圈上一段弧长近似当作平面长度，l³、l⁵ 项则是对投影变形的高阶补偿。经差越大，高阶项越不可忽略，这正是分带限制带宽的根本原因。

二、子午线弧长 X 的计算

X 是正算中量级最大、也最容易被忽视的一项。它表示从赤道沿中央子午线到纬度 B 的弧长，在 B = 45° 附近约为 4980 km，是 x 坐标的主体。X 不能写成简单的 a·B，因为子午线曲率半径随纬度变化。工程上通常采用按纬度幂级数展开的公式：

X = A₀·B − A₂/2·sin2B + A₄/4·sin4B − A₆/6·sin6B + A₈/8·sin8B

其中 A₀、A₂、A₄、A₆、A₈ 是由椭球参数 a、e² 决定的常数，
B 以弧度计。以 CGCS2000 椭球为例，A₀ ≈ 6367449.146 m。

实际作业中，X 一般直接查《高斯投影计算表》或由软件计算，手算时只需保证 B 的弧度化算正确。需要提醒的是：不同椭球（克拉索夫斯基、CGCS2000、WGS84）的 A₀ 等常数不同，同一 B 对应的 X 相差可达百米量级，因此必须确认所用椭球与矿区坐标系统一致，绝不能混用。

三、各阶系数的含义与取舍

公式中的系数并非随意拼凑，每一项都对应一种几何效应。t = tanB 反映纬度高低：低纬度处 t 小，高纬度处 t 迅速增大，因此高纬度地区高阶项的相对影响更明显。η² = e′²cos²B 反映椭球扁率带来的曲率差异，量级约 0.0067，在 l⁴、l⁵ 项中与 t² 共同作用。分母中的 2、6、24、120、720 来自泰勒展开的阶乘，说明阶次越高、贡献衰减越快。

精度取舍要看经差。以 3° 带为例，带边缘经差约 1.5°，此时 l⁵ 项对 y 的影响约为毫米级，l⁶ 项对 x 的影响更小，因此常规矿山控制测量保留到 l⁵（y）、l⁶（x）已足够。若在 6° 带边缘（l 可达 3°），l⁵ 项可达厘米级，必须保留；若只算到 l³，则会产生明显系统偏差。工程上判断的标准是：截断项的量级应小于允许误差的 1/3～1/5。

四、一个可复算的示例

设某矿点采用 CGCS2000 椭球、3° 带第 38 带，中央子午线 L₀ = 114°，已知 B = 30°00′00″、L = 114°45′00″，则 l = 0.75° = 0.01308997 rad。取 a = 6378137 m，e² = 0.0066943800229，算得 N ≈ 6383840 m，t = tan30° ≈ 0.577350，η² ≈ 0.005024。代入主项：y 主项 = N·cosB·l ≈ 6383840 × 0.866025 × 0.01308997 ≈ 72360 m；x 主项 X（30°）≈ 3320113 m，Δx 的 l² 项约 0.5×6383840×0.5×0.866025×0.00017135 ≈ 237 m。可见即使经差不到 1°，Δx 也达数百米，绝不能省略。上述数值为按公式手算的示例结果，实际作业应以软件或投影表为准并做校核。

五、结果校核与常见问题

算完后应做两项检查：一是用反算公式把（x, y）反求（B, L），往返较差应在毫米级；二是检查 y 的自然值是否在 ±334 km（3° 带）以内，超出说明带号或 L₀ 取错。常见错误包括：l 未化弧度导致结果放大 57 倍；把 y 的自然值直接当作含带号坐标使用（通用坐标需在 y 前加带号并加 500 km）；以及在高纬度矿区仍只取低阶项。把这些检查固化为计算流程的一部分，正算结果才可靠，也才能顺利衔接后续的换带计算。