让GIS三维可视化变得简单-投影坐标系统

2022-03-28 18:16:44 浏览数 (2)

接上文地理坐标系统,此文我们一起来了解下坐标系统中的投影坐标系统

投影坐标系统 PCS(Projection Coordinate System),它也叫非地球投影坐标系统 (notearth),或者再简单点叫平面坐标系统,也就是使用基于 X,Y 值的坐标系统来描述地球上某个点所处的位置

到了这你可能会问投影坐标系统和之前的地理坐标系统是什么关系呢?

其实,地理坐标系统是球面坐标,参考面是椭球面,坐标单位是经纬度,而投影坐标系是平面坐标系,参考平面是水平面,坐标单位是米、千米等等

「PS:」 因为地图相关太难画了,本文中所有图片皆来自网络,如有侵权,请告知,即删

什么是投影

投影指的是用一组光线将物体的形状投射到一个平面上去,称为投影,这个平面我们称之为投影面,且看下图,而我们要学习的投影,待投影的物体就是地球

能展开成平面的面有平面本身,球永远也展不开成一个平面,那么我们要把球搞成一个平面只有靠投影,而球的投影方式也有很多,我们后面会介绍到

为什么需要投影

需要投影是因为地理坐标系统中经纬度本身不带单位,度分秒仅仅是一个进制,另外,同样是1度经度,在不同的纬度带表示的弧段长是不一样的,这给面积以及定量计算等带来了困难

传统的地图(相较于google earth等3D平台)是呈现在二维的平面上的,所以需要一种能够将地球球面上三维坐标转化到二维平面的一种映射方法,这种方法就是地图投影,因此,地图投影是一种由经纬度 B,L,H 到投影坐标 X,Y,Z 的映射方式

地图投影的实质就是建立地球椭球表面上的点与地图平面上的点之间的对应关系,将建立在球体上的地理坐标系下的经纬度坐标,通过一种投影方法转为平面上的直角坐标

地图投影的用途,就是建立一种平面坐标系(称为投影坐标系),很多计算、业务在平面上计算相较与曲面上计算要简便一些

比如球面上,纬度相同,同样经度差的两点,处在不同的纬度上的距离就是不一样的,这给计算带来了不便,欧几里得平面上就不存在这个问题

很多经典的GIS软件,比如ArcMap等,基本框架,包括渲染、编辑等一些功能都是基于2D平面开发的,这也和我们对地图的要求一致

在允许变形的范围内,建立统一的平面坐标系,强调这一点,是因为在设计地图投影时,地图的范围也是需要考虑的重要因素

投影的分类

将球面上的目标展平到平面上,目标肯定会发生压缩或拉伸,根据变形的性质,地图投影可以简单分为下面几类

  • 等角投影:投影后目标在各个方向上变形一致,总体保持角度不变
  • 等距投影:投影后总体保持长度不变,面积和角度有些变形
  • 等积投影:假设目标上存在一个微分圆,如果投影后圆面积不变,此投影为等积投影
  • 任意投影:投影后面积、角度和长度均发生变化的投影

如下图

按照投影的方式,可以分为下面几类

  • 圆柱投影:投影面为圆柱
  • 圆锥投影:投影面为圆锥
  • 方位投影:投影面为平面

像后面所说的高斯克吕格投影就是一种圆柱投影,如下图所示

那么按照投影面(圆柱、圆锥,平面)主轴与地轴(连接南北极的短轴)相对位置的关系,投影又可以分为下面几种方式

  • 正轴投影:投影面中心轴与地轴相互重合
  • 斜轴投影:投影面中心轴与地轴斜向相交
  • 横轴投影:投影面中心轴与地轴相互垂直
  • 相切投影:投影面与椭球体相切
  • 相割投影:投影面与椭球体相割

如下图

常见投影方式

墨卡托投影-Mercator

墨卡托投影(Mercator),又称麦卡托投影、正轴等角圆柱投影,由荷兰地图学家墨卡托(G.Mercator)于1569年创拟

简单说就是假设地球被套在一个圆柱中,赤道与圆柱相切,然后在地球中心放一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,就形成以一幅墨卡托投影的世界地图,如下图

由于可显示任两点间的正确方位,航海用途的海图、航路图大部分都这种方式绘制,这种投影中线型比例尺在图中任意一点周围都保持不变,从而可以保持大陆轮廓投影后的角度和形状不变(即等角),但墨卡托投影会使面积产生变形,极点(也就是南北极地区)的比例甚至达到了无穷大

所以它并没完整展现了整个世界,地图上最高纬度是85.05度,为了简化计算,采用球形映射,而不是椭球体形状,采用墨卡托投影只是为了方便展示地图,我们还需要知道的是,这种映射会给Y轴方向带来 0.33% 的误差,如下图

高斯克吕格投影(Gauss - Kruger)

此投影是由德国数学家、物理学家、天文学家高斯于19世纪20年代拟定,后经德国大地测量学家克吕格于1912年对投影公式加以补充,故称为 高斯克吕格投影

你可以想象一下一个椭圆柱横着套在(相切于经线圈)地球的椭球面上,围绕地轴旋转,将地球表面投影到椭圆柱面,投影到椭圆柱面后,再沿着椭圆柱面的高剪开摊成平面,如下图

经度有东西各180度,按3度或者6度剥开一瓣用于投影这个叫分带,旋转高斯克吕格的投影面,可以获取下一个分带的投影,以合适全球所有的地方,如下图

高斯克吕格投影 又名 等角横切椭圆柱投影,在英美国家称为 横轴墨卡托投影,美国编制世界各地军用地图和地球资源卫星象片所采用的 全球/通用 横轴墨卡托投影(UTM)横轴墨卡托投影 的一种变形

通用横轴墨卡托投影(UTM)

UTM 全称 Universal Transverse Mercator 中文叫做 通用横轴墨卡托投影,和高斯克吕格投影非常相似,只不过它并不是切于经线圈,而是为椭圆柱横正轴割地球椭球体(穿地球而过),椭圆柱的中心线位于椭球体赤道面上,且通过椭球体质点,从而将椭球体上的点投影到椭圆柱上,就像下面这个样子

这种投影是美国陆军工程兵测绘局于20世纪40年代提出的,美国本土采用Clarke 1866椭球体,对全球其它地方采用 WGS84 椭球体,目前大部分的遥感影像图用的就是 UTM 投影出来的投影坐标系统

网络墨卡托投影(WebMercator)

网络墨卡托投影的英文名为 PVPM(Popular Visualization Pseudo Mercator),它不是严格意义的墨卡托投影,而是一个伪墨卡托的投影方法,由Google Map最先提出

它的不严谨性在于,在投影过程中,将表示地球的椭球面作为正球面处理

传说是因为谷歌程序员懒得用椭球面来编程计算屏幕坐标,具体不晓得

它以整个世界范围赤道作为标准纬线,本初子午线作为中央经线,两者交点为坐标原点,向东向北为正,向西向南为负

我们上面说过墨卡托投影方法,那么可以知道,当纬度接近两极,y值就趋向于无穷,于是这些 懒惰的工程师 就把 Y轴 的取值范围限定在 [-20037508.3427892,20037508.3427892] 之间

电子地图通常以经纬度为单位展示地理坐标,而墨卡托投影坐标经度范围:[-180,180] ,纬度上面已知不可能到达 90° 的,取值范围是 [-85.05112877980659,85.05112877980659] ,那其余的地区怎么办呢?没人知道,用不到就不管它

另外我们还需要知道的是,WebMercator投影 也是采用的 WGS84 基准面,基准面相关请看上文地理坐标系统

最后

常用的 高斯克吕格投影 的 PCS 族,在我国,地理坐标系统只有 北京54、西安80、国家2000,根据中央经线按需选用即可

基于 UTM投影 的 PCS 族,和 高斯克吕格投影 坐标系统族类似,只不过仅基于WGS84,也是根据中央经线按需选用即可

我国的高德、百度、腾讯地图,谷歌地图、OSM地图都用了网络墨卡托来平面化展示,不过国内高德百度腾讯和国内谷歌地图的地理坐标系统经过加密,而天地图虽然用的是国家 2000 地理坐标系,但是也提供了墨卡托投影方法来展示平面地图,ArcGIS Earth 则和 Google Earth 类似,都是直接使用3857坐标系

参考资料

地理信息系统导论-第九版[1]

聊聊GIS中的坐标系|再版[2]

你所不知的有趣投影方法[3]

Reference

[1]地理信息系统导论-第九版: https://book.douban.com/subject/31438879/

[2]聊聊GIS中的坐标系|再版: https://www.cnblogs.com/onsummer/p/12081889.html

[3]你所不知的有趣投影方法: https://blog.csdn.net/kikitamoon/article/details/50634331

0 人点赞