2. 高斯投影基本原理
高斯-克吕格投影是矿山测量中把椭球面坐标转换为平面坐标的基本工具。它的核心思想可以概括为:设想用一个横椭圆柱面套在参考椭球外面,使椭圆柱面与某一子午线(中央子午线)相切,然后将椭球面上的点按等角条件投影到椭圆柱面上,再把椭圆柱面沿母线剪开展平,得到平面坐标。理解这一过程,关键要抓住四个问题:等角条件意味着什么、投影带如何划分、中央子午线如何选择、投影变形沿什么规律分布。
一、等角投影条件
高斯投影属于正形投影(等角投影),其数学条件是:椭球面上任意两条微小线段的夹角,投影到平面后保持不变;同时,同一点处各方向的长度比与方向无关。用公式表达,即投影函数满足柯西-黎曼条件:
x = x(B, L), y = y(B, L)
∂x/∂L = ∂y/∂B, ∂x/∂B = -∂y/∂L
其中 B 为纬度,L 为经度差(相对中央子午线)。满足该条件后,平面上图形与椭球面图形保持相似,角度不变形,这对矿山测量尤为重要——井下导线、地表控制网的角度观测值可以直接参与平面计算,不必做角度改化。但要注意:等角不等于等长,长度变形依然存在,只是同一测站各方向的长度比相同,因而图形保持相似。
二、投影带划分:6度带与3度带
高斯投影只在中央子午线附近变形较小,离中央子午线越远,长度变形越大。为把变形限制在允许范围内,需按经差把椭球面划分成若干投影带,每带单独投影、单独使用中央子午线。我国常用两种分带:
| 分带方式 | 经差 | 带号与中央子午线关系 | 带号范围(我国) | 最大长度变形(约) |
|---|---|---|---|---|
| 6度带 | 6° | L₀ = 6N − 3 | 13~23 | 1/1000 量级 |
| 3度带 | 3° | L₀ = 3n | 25~45 | 1/10000 量级 |
6度带自 0° 子午线起,每 6° 一带,带号 N 从 1 开始,第 N 带的中央子午线为 6N−3。例如 N=20,中央子午线为 117°。3度带自 1.5° 起每 3° 一带,带号 n 的中央子午线为 3n,例如 n=39,中央子午线为 117°。可以看出,3度带的中央子午线有一半与 6度带重合(如 117°),另一半落在 6度带的分界子午线上(如 118.5°)。3度带带宽更窄,变形更小,常用于大比例尺测图和矿山控制测量;6度带多用于中小比例尺测图。
三、中央子午线的选择
中央子午线是投影的对称轴,也是长度比等于 1 的基准线。选择原则是:使测区尽量靠近中央子午线,且不跨越带的分界子午线。实际作业中,若测区位于标准分带边缘、变形超限,可采用任意带(自定义中央子午线),把中央子午线选在测区中央,例如矿区经度范围为 116°40′~117°20′,可取 117° 为中央子午线。任意带投影的坐标必须注明中央子午线经度,否则坐标无法唯一还原。
四、投影变形规律及其影响
高斯投影的长度变形由两部分组成:一是投影本身引起的变形,离中央子午线越远越大,近似与经差(或横坐标 y)的平方成正比;二是高程归化引起的变形,把地面长度归算到参考椭球面时长度缩短,近似与高程成正比。综合长度变形可近似表示为:
ΔS/S ≈ y²/(2R²) − H/R
式中 y 为距中央子午线的横坐标(自然值),R 为地球平均曲率半径(约 6371 km),H 为地面点高出椭球面的高程。第一项使长度变长,第二项使长度变短,二者符号相反,在一定条件下可以相互抵消。工程上通常要求投影长度变形不大于 1/40000(即 2.5 cm/km)。以 R=6371 km 估算,当 y≈50 km 时,第一项约为 1/32000,已接近限值;当 H≈160 m 时,第二项约为 1/40000。这说明:离中央子午线过远或高程过大的测区,必须采取任意带投影或高程归化改正。
对矿山测量的实际影响体现在三方面。第一,控制网坐标必须明确所属投影带和中央子午线,跨带或换带时不能直接混用坐标。第二,井上下对照、地表沉陷监测等需要高精度长度比对的场合,应优先采用 3度带或任意带,把测区控制在中央子午线两侧约 40~50 km 以内。第三,当矿区高程较大(如山区矿山)时,投影变形与高程归化变形可能叠加或抵消,应计算综合变形,必要时建立矿区独立坐标系(选择抵偿高程面或任意中央子午线),使关键巷道和地表监测线的长度变形满足限值要求。
综上,高斯投影以等角条件保证图形相似,以分带和中央子午线选择控制变形,以变形公式评估成果可用性。掌握这四点,才能正确理解后续的换带计算与矿区独立坐标系建立。