地图投影

投影的表面常见的有锥面、柱面与平面。

平面投影:两极,圆形区域;方向不变,最短距离是直线。

需要投影的理由很简单,我们看到的地图或者在计算机屏幕看到的地图都是平面的或者说是二维的,但是地球却不是平的。所以我们必须想出一种办法让地球表面上 的点跟平面上的点一一对应起来,而这种变换的结果就是把地球表面的点对应到笛卡儿坐标系统中。投影的方式主要有三种,见图6。每一种投影都会有不同程度的 变形,要么是长度变形、要么是角度变形、要么是面积变形。

为了应用于不同的目的,有些投影没有角度变形,但是其它两种变形就很变得很大,在航海应用中最为有名墨卡托投影就是一种无角度变形的投影,但是其面积变形很大,比如靠近北极地区的格陵兰岛的面积在墨卡托投影下看起来比南美洲还大。

墨卡托投影 ( Mercator Projection)

 墨卡托(Mercator)投影,是一种"等角正切圆柱投影”,荷兰地图学家墨卡托(Gerhardus Mercator 1512-1594)在1569年拟定,假设地球被围在一中空的圆柱里,其标准纬线与圆柱相切接触,然后再假想地球中心有一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。墨卡托投影没有角度变形, 由每一点向各方向的长度比相等,它的经纬线都是平行直线,且相交成直角,经线间隔相等,纬线间隔从标准纬线向两极逐渐增大。墨卡托投影的地图上长度和面积 变形明显,但标准纬线无变形,从标准纬线向两极变形逐渐增大,但因为它具有各个方向均等扩大的特性,保持了方向和相互位置关系的正确。在地图上保持方向和角度的正确是墨卡托投影的优点, 墨卡托投影地图常用作航海图和航空图,如果循着墨卡托投影图上两点间的直线航行,方向不变可以一直到达目的地,因此它对船舰在航行中定位、确定航向都具有 有利条件,给航海者带来很大方便。“海底地形图编绘规范”(GB/T 17834-1999,海军航保部起草)中规定1:25万及更小比例尺的海图采用墨卡托投影,其中基本比例尺海底地形图(1:5万,1:25 万,1:100万)采用统一基准纬线30°,非基本比例尺图以制图区域中纬为基准纬线。基准纬线取至整度或整分。

墨卡托投影坐标系取零子午线或自定义原点经线(L0)与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标X轴,赤道的投影为横坐标Y轴,构成墨卡托平面直角坐标系。

高斯-克吕格(Gauss-Kruger)投影和UTM(Universal Transverse Mercator)投影

高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”,又名"等角横切椭圆柱投影”, 地球椭球面和平面间正形投影的一种。德国数学家、物理学家、天文学家高斯(Carl FriedrichGauss,1777一 1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(Johannes Kruger,1857~1928)于 1912年对投影公式加以补充,故名。

该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线的条件,确定函数的形式,从而得到高斯一克吕格投影公式。投影后,除中央子午线和赤道为直线外, 其他子午线均为对称于中央子午线的曲线。 设想用一个椭圆柱横切于椭球面上投影带的中央子午线,按上述投影条件,将中央子午线两侧一定经差范围内的椭球面正形投影于椭圆柱面。将椭圆柱面沿过南北极 的母线剪开展平,即为高斯投影平面。取中央子午线与赤道交点的投影为原点,中央子午线的投影为纵坐标x轴,赤道的投影为横坐标y轴,构成高斯克吕格平面直 角坐标系。

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

按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差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)测图,城建坐标多采用三度带的高斯投影。

高斯-克吕格投影坐标

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

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

高斯-克吕格投影与UTM投影都是横轴墨卡托投影的变种。目前一些国外的软件或国外进口仪器的配套软件往往不支持高斯-克吕格投影,但支持UTM投影,因此常有把UTM投影当作高斯-克吕格投影的现象。


两者的区别:

  (1)投影几何方式不同。高斯-克吕格投影是“等角横切椭圆柱投影”,投影后中央经线保持长度不变,即比例系数为1;UTM投影是“等角横轴
割圆柱投影
”,圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条割线上没有变形,中央经线上长度比0.9996。

  (2)所选参考椭球不同。在我
国高斯投影采用1975年国际椭球体;UTM投影采用1866克拉克椭球。

  (3)分带方式不同。
 两者的分带起点不同,高斯-克吕格投影自0°子午线起每隔经差6度自西向东分带,第1带的中央经度为3°;UTM投影自西经180°起每隔经差6度自西
向东分带,第1带的中央经度为-177°。  投影带编号不同。UTM投影带号与高斯投影带号相差30。高斯-克吕格投影的第1带是UTM的第31带。
 此外,两投影的东伪偏移(False_Easting)都是500公里,高斯-克吕格投影北伪偏移(False_Northing)为零,UTM北半球
投影北伪偏移(False_Northing)为零,南半球则为10000km。

  (4)计算结果不同。两者主要差别在比例因子上,高斯-克吕格投影中央经
线上的比例系数为1,UTM投影为0.9996。  高斯-克吕格投影与UTM投影可近似采用
X[UTM]=0.9996*X[高斯],Y[UTM]=0.9996*Y[高斯],进行坐标转换(注意:如坐标纵轴西移了500km,转换时必须将Y值
减去500km乘上比例因子后再加500km)。

某些国外的软件如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度。
    (5)如正向投影,选择经纬度输入数据格式,有三个选项,缺省为十进制度格式。具体输入方式如下例:
   格 式     原始纬度值                  原始经度值          输入纬度值            输入经度值
十进制度     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
二、分带方法
    我国采用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万地形图)。

Reference

LongShaoAn, 关于墨卡托(Mercator)投影

http://hi.baidu.com/liongg/item/929eb7812c6b145d27ebd96b

3°、6°带高斯-克吕格投影

时间: 2024-11-04 21:18:20

地图投影的相关文章

[转]地图投影

1.地球椭球体     地球是一个表面很复杂的球体,人们以假想的平均静止的海水面形成的“大地体”为参照,推求出近似的椭球体,理论和实践证明,该椭球体近似一个以地球短轴为轴的椭园而旋转的椭球面,这个椭球面可用数学公式表达,将自然表面上的点归化到这个椭球面上,就可以计算了.下面列举了一些常用的一些椭球及参数:1)海福特椭球(1910) 我国52年以前采用的椭球 a=6378388m b=6356911.9461279m α=0.336700336702)克拉索夫斯基椭球(1940 Krassovsk

地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念

地球椭球体(Ellipsoid) 众所周知我们的地球表面是一个凸凹不平的表面,而对于地球测量而言,地表是一个无法用数学公式表达的曲面,这样的曲面不能作为测量和制图的基准面.假想一个扁率极小的椭圆,绕大地球体短轴旋转所形成的规则椭球体称之为地球椭球体.地球椭球体表面是一个规则的数学表面,可以用数学公式表达,所以在测量和制图中就用它替代地球的自然表面.因此就有了地球椭球体的概念. 地球椭球体有长半径和短半径之分,长半径(a)即赤道半径,短半径(b)即极半径.f=(a-b)/a为椭球体的扁率,表示椭球

地图投影系列介绍(三)----地图投影

4.地图投影     4.1 投影实质 将地球椭球面上的点映射到平面上的方法,称为地图投影. 为什么要进行投影? – 地理坐标为球面坐标,不方便进行距离.方位.面积等参数的量算. – 地球椭球体为不可展曲面. – 地图为平面,符合视觉心理,并易于进行距离.方位.面积等量算和各种空间分析.  投影的实质:经纬度坐标 -> 笛卡儿平面直角坐标系          建立地球椭球面上经纬线网和平面上相应经纬线网的数学基础,也就是建立地球椭球面上的点的地理坐标(λ,φ)与平面上对应点的平面坐标(x,y)之

地图投影系列介绍(四)----投影坐标系

4.4 我国常用地图投影  我国基本比例尺地形图(1:100万.1:50万.1:25万.1:10万.1:5万.1:2.5万.1:1万.1:5000)除1:100万以外均采用高斯-克吕格Gauss-Kruger投影(横轴等角切圆柱投影,又叫横轴墨卡托Transverse Mercator投影)为地理基础.  1:100万地形图采用兰伯特Lambert投影(正轴等角割圆锥投影),其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影保持一致.  海上小于50万的地形图多用墨卡托Me

MATLAB地图工具箱学习总结(一)从地图投影说起

前言 本学期地图投影课上,李连营老师建议我们使用MATLAB完成每周的作业.从大二上学期开始接触MATLAB学习数学运算和地理数据处理的我,自然不会放过这次机会,每周找了点时间好好研究了一下,把作业比较轻松地就搞定了.不过由于网上相关的教程.资料比较少,在夹杂了多种专业词汇的同时,又没有中文翻译,也是让人感到头疼.在此想总结一下相关的函数,以帮助其他需要使用工具箱的人.作为一个学生,我所接触的仅仅是皮毛,也必然会有一些错误,希望看到的人能够指出,谢谢. PS:本文基于MATLAB  R2014a

地图投影——高斯-克吕格投影、墨卡托投影和UTM投影

地图投影(Map Projection) http://baike.baidu.com/view/94066.htm 概念: 地图投影是把地球表面的任意点,利用一定数学法则,转换到地图平面上的理论和方法. 由于地球是一个赤道略宽两极略扁的不规则的梨形球体,故其表面是一个不可展平的曲面,所以运用任何数学方法进行这种转换都会产生误差和变形,为按照不同的需求缩小误差,就产生了各种投影方法. 方法: 1.几何透视法:几何透视法是利用透视的关系,将地球体面上的点投影到投影面上的一种投影方法.如假设地球按比

地图投影为什么

地图投影为什么 地图数学要素:地图投影.比例尺.控制点.坐标网.高程系.地图分幅等. 在我的印象中,比例尺从打小开始接触地图就强调其重要性,关联着距离量测.当时还有指北针等要注意的事项,主要关注于地图的使用.后来一不小心入了GIS的门,还得学会更深入的使用地图数据甚至是编制地图.这时候,大学GIS第一门课程<地图投影>就来了.当时的注意点完全被繁琐的公式迷惑,不过看着用C语言在那小黑块的屏幕画出一幅幅漂亮的投影地图,还是相当的快乐.时至今日,又多了解些相关知识,重新回顾复习整理下. 在个人的认

Google Maps地图投影全解析(3):WKT形式表示

update20090601:EPSG对该投影的编号设定为EPSG:3857,对应的WKT也发生了变化,下文不再修改,相对来说格式都是那样,可以到http://www.epsg-registry.org 网站输入SRID进行查询.         Google Maps和Virtual Earth等的流行程度不用多讲,然而他们所使用的Web Mercator或Spherical Mercator在很长一段时间内并没有被EPSG的投影数据库所接纳.EPSG认为它不能算作科学意义上的投影,所以只是给

关于科学选择地图投影类型的探讨

关于科学选择地图投影类型的探讨 作者:左正康(11级GIS 2011221108120027) 摘要:地图投影是利用一定数学法则把地球表面的经纬线转换到平面上的理论和方法.由于地球是一个赤道略宽两级略扁的不规则的梨形球体,故其表面是一个不可展平的曲面,所以运用任何数学方法进行这种转换都会产生误差和变形,为按照不同的需求缩小误差,就产生了各种投影方法.投影方法按变形方式可以分为等角投影.等积投影.等距投影:按正轴投影时经纬网的形状可以分为几何投影(平面投影.圆锥投影.圆柱投影.多圆锥投影,几何投影

D3学习之:D3.js中的12中地图投影方式

特别感谢:1.[张天旭]的D3API汉化说明,已被引用到官方网站: 2.[馒头华华]提供的ourd3js.com上提供的学习系列教程,让我们这些新人起码有了一个方向. 不得不说,学习国外的新技术真的是一个很艰苦的过程. 在学习D3绘制地图的过程中,有朋友建议看一下其中投影的说明比较好,于是,凭借我这半吊子不到的英文水平,大致给翻译了下来,仅供参考: 原文链接:https://github.com/mbostock/d3/wiki/Geo-Projections#albers D3中一共提供了12