地图云
当前位置:地图知识 / GIS知识 / 坐标投影 / 详情

高斯-克吕格投影(等角横切椭圆柱投影)

收藏

高斯-克吕格投影是由德国数学家、物理学家、天文学家高斯于19 世纪20 年代拟定,后经德国大地测量学家克吕格于1912 年对投影公式加以补充,故称为高斯-克吕格投影,又名"等角横切椭圆柱投影”,是地球椭球面和平面间正形投影的一种。

定义

高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”,又名"等角横切椭圆柱投影”,地球椭球面和平面间正形投影的一种。

德国数学家、物理学家、天文学家高斯(Carl Friedrich Gauss,1777一1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(Johannes Kruger,1857~1928)于1912年对投影公式加以补充,故名。该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线的条件,确定函数的形式,从而得到高斯一克吕格投影公式。投影后,除中央子午线和赤道为直线外,其他子午线均为对称于中央子午线的曲线。设想用一个椭圆柱横切于椭球面上投影带的中央子午线,按上述投影条件,将中央子午线两侧一定经差范围内的椭球面正形投影于椭圆柱面。将椭圆柱面沿过南北极的母线剪开展平,即为高斯投影平面。取中央子午线与赤道交点的投影为原点,中央子午线的投影为纵坐标x轴,赤道的投影为横坐标y轴,构成高斯克吕格平面直角坐标系。

高斯-克吕格投影在长度和面积上变形很小,中央经线无变形,自中央经线向投影带边缘,变形逐渐增加,变形最大之处在投影带内赤道的两端。由于其投影精度高,变形小,而且计算简便(各投影带坐标一致,只要算出一个带的数据,其他各带都能应用),因此在大比例尺地形图中应用,可以满足军事上各种需要,能在图上进行精确的量测计算。

它是假设一个椭圆柱面与地球椭球体面横切于某一条经线上,按照等角条件将中央经线东、西各3°或1.5°经线范围内的经纬线投影到椭圆柱面上,然后将椭圆柱面展开成平面而成的。

这种投影,将中央经线投影为直线,其长度没有变形,与球面实际长度相等,其余经线为向极点收敛的弧线,距中央经线愈远,变形愈大。 赤道线投影后是直线,但有长度变形。除赤道外的其余纬线,投影后为凸向赤道的曲线,并以赤道为对称轴。经线和纬线投影后仍然保持正交。所有长度变形的线段,其长度变形比均大于1. 随远离中央经线,面积变形也愈大。若采用分带投影的方法,可使投影边缘的变形不致过大。

我国各种大、中比例尺地形图采用了不同的高斯-克吕格投影带。其中大于1:1万的地形图采用3°带;1:2.5万至1:50万的地形图采用6°带。

高斯克吕格投影这一投影的几何概念是,假想有一个椭圆柱与地球椭球体上某一经线相切,其椭圆柱的中心轴与赤道平面重合,将地球椭球体面有条件地投影到椭球圆柱面上高斯克吕格投影条件:

a) 中央经线和赤道投影为互相垂直的直线,且为投影的对称轴;

b) 具有等角投影的性质;

c) 中央经线投影后保持长度不变;

假想有一个椭圆柱面横套在地球椭球体外面,并与某一条子午线(此子午线称为中央子午线或轴子午线)相切,椭圆柱的中心轴通过椭球体中心,然后用一定投影方法,将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上,再将此柱面展开即成为投影面,此投影为高斯投影。高斯投影是正形投影的一种。

bd315c6034a85edf3606297449540923dc5475db.jpg

 

分带投影

按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自0度子午线起每隔经差6度自西向东分带,带号依次编为第 1、2…60带。三度带是在六度带的基础上分成的,它的中央子午线与六度带的中央子午线和分带子午线重合,即自 1.5度子午线起每隔经差3度自西向东分带,带号依次编为三度带第1、2…120带。我国的经度范围西起 73°东至135°,可分成六度带十一个,各带中央经线依次为75°、81°、87°、……、117°、123°、129°、135°,或三度带二十二个。六度带可用于中小比例尺(如 1:250000)测图,三度带可用于大比例尺(如 1:10000)测图,城建坐标多采用三度带的高斯投影。

高斯投影6度带

自0子午线起每隔经差6自西向东分带,依次编号1,2,3,…。我国6度带中央子午线的经度,由69°起每隔6°而至135°,共计12带(12~23带),带号用N表示,中央子午线的经度用L0表示,它们的关系是:L0=6N-3。

 
d967f3a0-006b-4afe-a193-a029852195f2.jpeg
 

高斯投影3度带

它的中央子午线一部分同6度带中央子午线重合,一部分同6度带的分界子午线重合,如用n表示3度带的带号,表示带中央子午线经度,它们的关系图所示。我国带共计22带(24~45带)。

高斯平面直角坐标系

在投影面上,中央子午线和赤道的投影都是直线,并且以中央子午线和赤道的交点0作为坐标原点,以中央子午线的投影为纵坐标轴(X),以赤道的投影为横坐标轴(Y),这样便形成了高斯平面直角坐标系。 

在我国x坐标都是正的,y坐标的最大值(在赤道上,6°带)约为330km。为了避免出现负的横坐标,可在横坐标上加上500 000m。此外还应在坐标前面再冠以带号。这种坐标称为国家统一坐标。

例如,有一点y=19 623 456.789m,该点位在19带内,位于中央子午线以东,其相对于中央子午线而言的横坐标则是:首先去掉带号,再减去500 000m,最后得=123 456.789m。

高斯平面投影的特点

①中央子午线无变形;

②无角度变形,图形保持相似;

③离中央子午线越远,变形越大。

椭球面三角系归算高斯投影面

将椭球面三角系归算到高斯投影面的主要内容是:

(1)将起始点p的大地坐标(L,B)归算为高斯平面直角坐标(X,Y);为了检核还应进行反算,亦即根据X,Y反算L,B。

(2)通过计算该点的子午线收敛角γ及方向δ改正,将椭球面上起算边大地方位角A归算到高斯平面上相应边PK的坐标方位角α。

(3)通过计算各方向的曲率改正和方向改正,将椭球面上各三角形内角归算到高斯平面上的由相应直线组成的三角形内角。

(4)通过计算距离改正Δs,将椭球面上起算边PK的长度S归算到高斯平面上的直线长度s。

(5)当控制网跨越两个相邻投影带,需要进行平面坐标的邻带换算。

高斯-克吕格投影的变形分析

①中央经线上无变形,满足投影后长度比不变的条件;

②除中央经线上长度比为1以外,其它任何点长度比均大于1;

③在同一条纬线上,离中央经线越远则变形越大,最大值位于投影带边缘。

④在同一条经线上,纬度越低变形越大,最大值位于赤道上。

⑤等角投影无角度变形,面积比为长度比的平方。

⑥长度比的等变形线平行于中央轴子午线。

高斯-克吕格投影的特点

  • 投影后中央子午线为直线,长度不变形,其余经线投影对称并且凹向中央子午线,离中央子午线越远,变形越大
  • 赤道的投影也是直线,并与中央子午线正交,其余纬线投影为凸向赤道的对称曲线
  • 投影后经纬线仍然保持相会垂直,角度无变形
  • 长度和面积变形很小

高斯-克吕格投影坐标

高斯-克吕格投影是按分带方法各自进行投影,故各带坐标成独立系统。以中央经线投影为纵轴(x),赤道投影为横轴(y),两轴交点即为各带的坐标原点。纵坐标以赤道为零起算,赤道以北为正,以南为负。我国位于北半球,纵坐标均为正值。横坐标如以中央经线为零起算,中央经线以东为正,以西为负,横坐标出现负值,使用不便,故规定将坐标纵轴西移500公里当作起始轴,凡是带内的横坐标值均加500公里。由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。   

高斯-克吕格投影与UTM投影

某些国外的软件如ARC/INFO或国外仪器的配套软件如多波束的数据处理软件等,往往不支持高斯-克吕格投影,但支持UTM投影,因此常有把UTM投影坐标当作高斯-克吕格投影坐标提交的现象。

UTM投影全称为“通用横轴墨卡托投影”,是等角横轴割圆柱投影(高斯-克吕格为等角横轴切圆柱投影),圆柱割地球于南纬80度、北纬84度两条等高圈,该投影将地球划分为60个投影带,每带经差为6度,已被许多国家作为地形图的数学基础。UTM投影与高斯投影的主要区别在南北格网线的比例系数上,高斯-克吕格投影的中央经线投影后保持长度不变,即比例系数为1,而UTM投影的比例系数为0.9996。UTM投影沿每一条南北格网线比例系数为常数,在东西方向则为变数,中心格网线的比例系数为0.9996,在南北纵行最宽部分的边缘上距离中心点大约 363公里,比例系数为 1.00158。

高斯-克吕格投影与UTM投影可近似采用 Xutm=0.9996 * X高斯,Yutm=0.9996 * Y高斯进行坐标转换。以下举例说明(基准面为WGS84):

 

输入坐标(度)

高斯投影(米)

UTM投影(米)

Xutm=0.9996 * X高斯, Yutm=0.9996 * Y高斯

纬度值(X

32

3543600.9

3542183.5

3543600.9*0.9996 ≈ 3542183.5

经度值(Y

121

21310996.8

311072.4

(310996.8-500000)*0.9996+500000 ≈ 311072.4

 

注:坐标点(32,121)位于高斯投影的21带,高斯投影Y值21310996.8中前两位“21”为带号;坐标点(32,121)位于UTM投影的51带,上表中UTM投影的Y值没加带号。因坐标纵轴西移了500000米,转换时必须将Y值减去500000乘上比例因子后再加500000。