3度6度带高斯投影详解.docx
《3度6度带高斯投影详解.docx》由会员分享,可在线阅读,更多相关《3度6度带高斯投影详解.docx(9页珍藏版)》请在冰豆网上搜索。
![3度6度带高斯投影详解.docx](https://file1.bdocx.com/fileroot1/2023-2/2/8c46429a-aa32-4d9e-ac18-e5e82280b7a9/8c46429a-aa32-4d9e-ac18-e5e82280b7a91.gif)
3度6度带高斯投影详解
3度6度带高斯投影
选择投影的目的在于使所选投影的性质、特点适合于地图的用途,同时考虑地图在图廓范围内变形较小而且变形分布均匀。
海域使用的地图多采用保角投影,因其能保持方位角度的正确。
我国的基本比例尺地形图(1:
5千,1:
1万,1:
2.5万,1:
5万,1:
10万,1:
25万,1:
50万,1:
100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),这是一个等角横切椭圆柱投影,又叫横轴墨卡托投影(TransverseMercator);小于50万的地形图采用等角正轴割园锥投影,又叫兰勃特投影(LambertConformalConic);海上小于50万的地形图多用等角正轴圆柱投影,又叫墨卡托投影(Mercator)。
一般应该采用与我国基本比例尺地形图系列一致的地图投影系统。
地图坐标系由大地基准面和地图投影确定,大地基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的大地基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。
我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的IAG75地球椭球体建立了我国新的大地坐标系--西安80坐标系,目前GPS定位所得出的结果都属于WGS84坐标系统,WGS84基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心的坐标系。
因此相对同一地理位置,不同的大地基准面,它们的经纬度坐标是有差异的。
采用的3个椭球体参数如下(源自“全球定位系统测量规范GB/T8314-2001”):
椭球体
长半轴
短半轴
Krassovsky
6378245
6356863.0188
IAG75
6378140
6356755.2882
WGS84
6378137
6356752.3142
椭球体与大地基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的Pulkovo1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体,但它们的大地基准面显然是不同的。
在目前的GIS商用软件中,大地基准面都通过当地基准面向WGS84的转换7参数来定义,即三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值;三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时,分别绕Xt、Yt、Zt的旋转角;最后是比例校正因子,用于调整椭球大小。
北京54、西安80相对WGS84的转换参数至今没有公开,实际工作中可利用工作区内已知的北京54或西安80坐标控制点进行与WGS84坐标值的转换,在只有一个已知控制点的情况下(往往如此),用已知点的北京54与WGS84坐标之差作为平移参数,当工作区范围不大时,如青岛市,精度也足够了。
以(32°,121°)的高斯-克吕格投影结果为例,北京54及WGS84基准面,两者投影结果在南北方向差距约63米(见下表),对于几十或几百万的地图来说,这一误差无足轻重,但在工程地图中还是应该加以考虑的。
输入坐标(度)
北京54高斯投影(米)
WGS84高斯投影(米)
纬度值(X)
32
3543664
3543601
经度值(Y)
121
21310994
21310997
高斯-克吕格投影
(1)高斯-克吕格投影性质
高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”,又名"等角横切椭圆柱投影”,地球椭球面和平面间正形投影的一种。
德国数学家、物理学家、天文学家高斯(CarlFriedrichGauss,1777一1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(JohannesKruger,1857~1928)于1912年对投影公式加以补充,故名。
该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线的条件,确定函数的形式,从而得到高斯一克吕格投影公式。
投影后,除中央子午线和赤道为直线外,其他子午线均为对称于中央子午线的曲线。
设想用一个椭圆柱横切于椭球面上投影带的中央子午线,按上述投影条件,将中央子午线两侧一定经差范围内的椭球面正形投影于椭圆柱面。
将椭圆柱面沿过南北极的母线剪开展平,即为高斯投影平面。
取中央子午线与赤道交点的投影为原点,中央子午线的投影为纵坐标x轴,赤道的投影为横坐标y轴,构成高斯克吕格平面直角坐标系。
高斯-克吕格投影在长度和面积上变形很小,中央经线无变形,自中央经线向投影带边缘,变形逐渐增加,变形最大之处在投影带内赤道的两端。
由于其投影精度高,变形小,而且计算简便(各投影带坐标一致,只要算出一个带的数据,其他各带都能应用),因此在大比例尺地形图中应用,可以满足军事上各种需要,能在图上进行精确的量测计算。
(2)高斯-克吕格投影分带
按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。
分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。
通常按经差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)测图,城建坐标多采用三度带的高斯投影。
(3)高斯-克吕格投影坐标
高斯-克吕格投影是按分带方法各自进行投影,故各带坐标成独立系统。
以中央经线投影为纵轴(x),赤道投影为横轴(y),两轴交点即为各带的坐标原点。
纵坐标以赤道为零起算,赤道以北为正,以南为负。
我国位于北半球,纵坐标均为正值。
横坐标如以中央经线为零起算,中央经线以东为正,以西为负,横坐标出现负值,使用不便,故规定将坐标纵轴西移500公里当作起始轴,凡是带内的横坐标值均加500公里。
由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。
(4)高斯-克吕格投影与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。
单点转换步骤如下:
(1)选择是高斯正转换还是反转换,缺省为经纬度转换到高斯投影坐标,投影坐标单位为米。
(2)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。
(3)选择分带,3度或6度,缺省为6度。
(4)输入中央经度,20带(114°E~120°E)中央经度为117度,21带(120°E~126°E)中央经度为123度。
如正向投影,选择经纬度输入数据格式,有三个选项,缺省为十进制度格式。
具体输入方式如下例:
格式
原始纬度值
原始经度值
输入纬度值
输入经度值
十进制度
35.445901°
122.997344°
35.445901
122.997344
度分
35°26.7541′
122°59.8406′
3526.7541
12259.8406
度分秒
35°26′45.245″
122°59′50.438″
352645.245
1225950.438
(6)正投影按选定格式在“输入”栏输入经纬度值,反投影输入以米为单位的X、Y坐标值。
(7)单击“单点转换”按钮。
(8)在“输出”栏查看计算结果。
批量转换步骤如下:
(1)准备好需要转换的输入数据文件,要求是文本文件,分两列,第一列纬度值或纵向坐标值,第二列经度值或横向坐标值,两列之间用空格分开。
正向投影时,纬度值及经度值格式可以有三种选择,缺省当作十进制度处理;反向投影时,纵向及横向坐标值必须以米为单位。
下例为度分秒格式(WGS84)的6°带正投影输入数据文件testdata.txt
352645.245 1225950.438
353800.402 1230000.378
351600.519 1225959.506
345800.101 1225959.8
343600.336 1230000.26
341400.018 1225959.897
335159.17 1225959.46
333000.08 1230000.28
(2)选择是高斯正转换还是反转换,缺省为经纬度转换到高斯投影坐标,投影坐标单位为米。
(3)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。
(4)选择分带,3度或6度,缺省为6度。
(5)输入中央经度,20带(114°E~120°E)中央经度为117度,21带(120°E~126°E)中央经度为123度。
(6)如正向投影,选择输入数据文件中的经纬度输入数据格式,有三个选项,缺省为十进制度格式。
(7)单击“批量转换”按钮。
弹出打开文件对话框,输入你的数据文件名。
(8)输入转换结果文件名,单击“保存”后,程序开始进行计算。
(9)打开输出文件查看计算结果,结果分五列,第一序号,第二列输入纬度值或纵向坐标值,第三列输入经度值或横向坐标值,第四列转换后纬度值或纵向坐标值,第五列转换后经度值或横向坐标值。
下例为度分秒格式(WGS84)的6°带正投影转换结果数据文件result.txt
1 352645.245 1225950.438 3924063.3 21499758.9
2 353800.402 1230000.378 3944871.4 21500009.5
3 351600.519 1225959.506 3904193.8 21499987.5
4 345800.101 1225959.8 3870898.1 21499994.9
5 343600.336 1230000.26 3830228.5 21500006.6
6 341400.018 1225959.897 3789544.4 21499997.4
7 335159.17 1225959.46 3748846.4 21499986.1
8 333000.08 1230000.28 3708205 21500007.2
\
[个人心得]快速区分区分3度带和6度带
也许实际工作中你能用到,哈哈
整理一个。
1.我国采用6度分带和3度分带:
1∶2.5万及1∶5万的地形图采用6度分带投影,即经差为6度,从零度子午线开始,自西向东每个经差6度为一投影带,全球共分60个带,用1,2,3,4,5,……表示.即东经0~6度为第一带,其中央经线的经度为东经3度,东经6~12度为第二带,其中央经线的经度为9度。
1∶1万的地形图采用3度分带,从东经1.5度的经线开始,每隔3度为一带,用1,2,3,……表示,全球共划分120个投影带,即东经1.5~4.5度为第1带,其中央经线的经度为东经3度,东经4.5~7.5度为第2带,其中央经线的经度为东经6度.
地形图上公里网横坐标前2位就是带号,例如:
1∶5万地形图上的横坐标为20345486,其中20即为带号,345486为横坐标值。
2.当地中央经线经度的计算
六度带中央经线经度的计算:
当地中央经线经度=6°×当地带号-3°,例如:
地形图上的横坐标为20345,其所处的六度带的中央经线经度为:
6°×20-3°=117°(适用于1∶2.5万和1∶5万地形图)。
三度带中央经线经度的计算:
中央经线经度=3°×当地带号(适用于1∶1万地形图)。
一个相当EASY的方法:
在PRC陆地范围内,坐标(Y坐标,8位数,前两位是带号)带号小于 等于23的肯定是6度带,大于等于24的肯定是3度带.
3.只知道经纬度时中央经线的计算
将当地经线的整数部分除以6,再取商的整数部分加上1.再将所得结果乘以6后减去3,就可以得到当地的中央经线值。
6度分带和3度分带是怎么回事?
1.我国采用6度分带和3度分带:
1∶2.5万及1∶5万的地形图采用6度分带投影,即经差为6度,从零度子午线开始,自西向东每个经差6度为一投影带,全球共分60个带,用1,2,3,4,5,……表示.即东经0~6度为第一带,其中央经线的经度为东经3度,东经6~12度为第二带,其中央经线的经度为9度。
1∶1万的地形图采用3度分带,从东经1.5度的经线开始,每隔3度为一带,用1,2,3,……表示,全球共划分120个投影带,即东经1.5~4.5度为第1带,其中央经线的经度为东经3度,东经4.5~7.5度为第2带,其中央经线的经度为东经6度.我省位于东经113度-东经120度之间,跨第38、39、40共计3个带,其中东经115.5度以西为第38带,其中央经线为东经114度;东经115.5~118.5度为39带,其中央经线为东经117度;东经118.5度以东到山海关为40带,其中央经线为东经120度。
地形图上公里网横坐标前2位就是带号,例如:
1∶5万地形图上的横坐标为20345486,其中20即为带号,345486为横坐标值。
2.当地中央经线经度的计算
六度带中央经线经度的计算:
当地中央经线经度=6°×当地带号-3°,例如:
地形图上的横坐标为20345,其所处的六度带的中央经线经度为:
6°×20-3°=117°(适用于1∶2.5万和1∶5万地形图)。
三度带中央经线经度的计算:
中央经线经度=3°×当地带号(适用于1∶1万地形图)。
如有侵权请联系告知删除,感谢你们的配合!