双标准纬线等角圆锥投影_坐标转换—墨卡托投影、高斯克吕格投影、UTM投影

引子：上周“勘探风暴“的一个同仁非要我帮他转15个坐标。从1.5度带转为3度带坐标，说实话我不是测量专业的，稍有接触，但并不精于此。但在手头参数有限的情况下帮他完成了工作。搞地质、搞物探的人难免要遇到坐标的问题，今天就粗略讲一下，有很多资料来自网络，请大家自己多多查证推敲，有兴趣深度学习的可以找我要以下书目的视频课程。

一、投影简介 
      将地球椭球面上的点映射到平面上的数学方法，叫做地图投影。地图投影的原因如下：
1、地理坐标为球面坐标，不方便进行距离、方位、和面积的量算；
2、地球椭球面为不可展面，无法制图；
3、地图为平面，符合视觉心理，且易于进行距离、方位、和面积的量算和各种空间分析。

  

 二、墨卡托投影  
       墨卡托(Mercator)投影，是一种”等角正切圆柱投影”，荷兰地图学家墨卡托(Gerhardus Mercator 1512－1594)在1569年拟定，假设地球被围在一中空的圆柱里，其标准纬线与圆柱相切接触，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱体上，再把圆柱体展开，这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。墨卡托投影没有角度变形，由每一点向各方向的长度比相等，它的经纬线都是平行直线，且相交成直角，经线间隔相等，纬线间隔从标准纬线向两极逐渐增大。墨卡托投影的地图上长度和面积变形明显，但标准纬线无变形，从标准纬线向两极变形逐渐增大，但因为它具有各个方向均等扩大的特性，保持了方向和相互位置关系的正确。在地图上保持方向和角度的正确是墨卡托投影的优点，墨卡托投影地图常用作航海图和航空图，如果循着墨卡托投影图上两点间的直线航行，方向不变可以一直到达目的地，因此它对船舰在航行中定位、确定航向都具有有利条件，给航海者带来很大方便。“海底地形图编绘规范”(GB/T 17834-1999，海军航保部起草)中规定1：25万及更小比例尺的海图采用墨卡托投影，其中基本比例尺海底地形图(1：5万，1：25万，1：100万)采用统一基准纬线30°，非基本比例尺图以制图区域中纬为基准纬线。基准纬线取至整度或整分。
      墨卡托投影坐标系取零子午线或自定义原点经线(L0)与赤道交点的投影为原点，零子午线或自定义原点经线的投影为纵坐标X轴，赤道的投影为横坐标Y轴，构成墨卡托平面直角坐标系。三、高斯-克吕格(Gauss-Kruger)投影 
      高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”，又名”等角横切椭圆柱投影”，地球椭球面和平面间正形投影的一种。德国数学家、物理学家、天文学家高斯(Carl FriedrichＧauss，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)测图，城建坐标多采用三度带的高斯投影。

以下红色的是干货：

      我在这里支出注意的是，我们在使用6度带或者3度带的时候都要涉及到中央子午线以及带号的计算。使用的时候拿到数据第一件事情，先看X，Y坐标的整数位数，有以下两种情况：

      1、X坐标6位，Y坐标7位(东阳何生的经验总结)

     没有加带号的坐标，坐标单位是米，假偏东500公里。

      2、X坐标8位，Y坐标7位

      加了带号的坐标，坐标单位是米。X坐标最前面两位就是添加的带号，这时就要判断是3度带还是6度带，我国幅员辽阔，经度从东经72度到135度，有经验的人一看带号就能大致知道是6度分带还是3度分带；没有经验的，就随便假设一个，然后根据下面的公式算出其中央经线，再与研究区域所在的经度对照，看是否相符，从而判断出是3度分带还是6度分带。带号与中央经线一一对应，知道两者中的任何一个，都能推算出另外一个的值，计算公式如下：

我国1:2.5—1:50万地形图均采用6度分带；1:1万及更大比例尺地形图采用3度分带，以保证必要的精度。我国的经度范围西起 73°东至135°，横跨11个六度带，73~135，对应6度带号是13~23。三度带比6度带大一倍，基本上就是24到45。当地中央经线和带号计算公式 首先我们要确定我们的坐标是3度带的还是6度带的投影，然后根据经度来计算带号、中央子午线，计算公式如下：6度带：中央子午线计算公式：中央子午线L=6 ×(N＋1)－3 。N=[当地经度/6]，N值不进行四舍五入，只取整数部分，(N＋1)即为6度带的带号。3度带：中央子午线计算公式：中央子午线L=3 ×N 。N=当地经度/3，N值进行四舍五入后即为3度带的带号。

四、UTM投影 
      某些国外的软件如ARC/INFO或国外仪器的配套软件如多波束的数据处理软件等，往往不支持高斯-克吕格投影，但支持UTM投影，因此常有把UTM投影坐标当作高斯-克吕格投影坐标提交的现象。
      UTM投影全称为“通用横轴墨卡托投影”，是等角横轴割圆柱投影(高斯-克吕格为等角横轴切圆柱投影)，圆柱割地球于南纬80度、北纬84度两条等高圈，该投影将地球划分为60个投影带，每带经差为6度，已被许多国家作为地形图的数学基础。UTM投影与高斯投影的主要区别在南北格网线的比例系数上，高斯-克吕格投影的中央经线投影后保持长度不变，即比例系数为1，而UTM投影的比例系数为0.9996。UTM投影沿每一条南北格网线比例系数为常数，在东西方向则为变数，中心格网线的比例系数为0.9996，在南北纵行最宽部分的边缘上距离中心点大约 363公里，比例系数为 1.00158。 注：坐标点(32,121)位于高斯投影的21带，高斯投影Y值21310996.8中前两位“21”为带号；坐标点(32,121)位于UTM投影的51带，上表中UTM投影的Y值没加带号。因坐标纵轴西移了500000米，转换时必须将Y值减去500000乘上比例因子后再加500000。       高斯-克吕格投影与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	Xutm=3543600.9*0.9996 ≈ 3542183.5
经度值(Y)121	21310996.8	311072.4	Yutm=(310996.8-500000)*0.9996+500000≈311072.4

五、关于1.5度带转换为3度带坐标       好多人不了解1.5度带，1.5度带属于地方坐标系，跟3、6度带一样采用高斯投影，只不过带宽变成1.5 ，中央子午线为 108、109.5、111、112.5、114等。东坐标常数不变，500000。坐标换带计算，就是选择对应中央子午线就行，例如1.5是109.5，三度带就可能是108，六度带就可能是111。

        转换的秘诀：一般坐标换带软件都可以的，只要能自定义中央子午线就行，中央子午线的计算方法道理都一样。选自定义高斯投影，填上合适的中央子午线就可以转换了。

       当然，你要十分精确的转换，可能要涉及到比较严密的是用七参数布尔莎模型，即X平移，Y平移，Z平移，X旋转(WX)，Y旋转(WY)，Z旋转(WZ)，尺度变化(DM)。要求得七参数就需要在一个地区有3个以上的已知点。

希望大家多多支持“勘探风暴”的发展。我珍惜每一次与大家交流沟通的机会，用心做好每一个问题的解答，尽我所能。