墨卡托投影与瓦片地图

最新推荐文章于 2024-05-04 20:47:30 发布

LuossD

最新推荐文章于 2024-05-04 20:47:30 发布

阅读量4.9k

点赞数 22

分类专栏： QT 文章标签：其他

本文链接：https://blog.csdn.net/weixin_42582542/article/details/121080586

版权

QT 专栏收录该内容

4 篇文章 2 订阅

订阅专栏

一、开胃小知识

先写几个小知识点：

1）既然地球是圆的，但是为什么地面却是平的呢？这主要是因为人和地球比起来，相差悬殊；人眼看到的地球表面区域弧度非常小，人类察觉不到，从而产生地面平坦的错觉，我们因此会认为“地面是平面”了。

2）纬度是指某点与地球球心的连线和地球赤道面所成的线面角，其数值在0至90度之间。经纬度是经度与纬度的合称组成一个坐标系统。又称为地理坐标系统。东经正数，西经为负数；北纬为正数，南纬为负数。

3）东西半球分界线不是本初子午线；而是东经160度和与其相对的西经20度所组成的经线圈。至于为啥可参考这篇知乎：为什么东西半球的分界线不是本初子午线，而南北半球的界线却是赤道？ - 知乎

二、墨卡托投影

1、什么是墨卡托投影？

我们平时看到的地图大多是平面的，那么如何将一个三维的地球画在平面上呢？

墨卡托投影就是将三维的地球表示在一个二维平面上的方法之一，也是应用得最广泛的方法。我们平时看到的谷歌地图，百度地图，包括我们的GeoHey地图，都是使用的（WEB）墨卡托投影。

墨卡托投影的过程其实非常简单，就是将地球展开成一个圆柱，再将圆柱展开成平面。

从球心出发射出一条直线，它与球的交点投影后的位置就是这条线与圆柱的交点。当然，中间的计算过程中会做一些取舍。

引用来源：一个小游戏让你彻底弄懂墨卡托投影 - 知乎

细说墨卡托

墨卡托(Mercator)投影，是一种"等角正切圆柱投影”，荷兰地图学家墨卡托（Gerhardus Mercator 1512－1594）在1569年拟定，假设地球被围在一中空的圆柱里，其标准纬线与圆柱相切接触，然后再假想地球中心有一盏灯，把球面上的图形投影到圆柱体上，再把圆柱体展开，这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。墨卡托投影没有角度变形，由每一点向各方向的长度比相等，它的经纬线都是平行直线，且相交成直角，经线间隔相等，纬线间隔从标准纬线向两极逐渐增大。墨卡托投影的地图上长度和面积变形明显，但标准纬线无变形，从标准纬线向两极变形逐渐增大，但因为它具有各个方向均等扩大的特性，保持了方向和相互位置关系的正确。在地图上保持方向和角度的正确是墨卡托投影的优点，墨卡托投影地图常用作航海图和航空图，如果循着墨卡托投影图上两点间的直线航行，方向不变可以一直到达目的地，因此它对船舰在航行中定位、确定航向都具有有利条件，给航海者带来很大方便。

引用来源：投影和坐标系统那点事 - 简书

2、墨卡托投影的特点

墨卡托投影有以下特点：

原来的经线（子午线）都汇聚于极点，纬度愈高经线的间隔愈小，而在墨卡托投影上，经线全部画成与赤道上间隔相等的平行直线。
地球表面上假想的各纬线间隔基本是相等的，而在墨卡托投影上纬线间隔由赤道向两极逐渐扩大，而且经纬线以同一倍数扩大。
墨卡托投影上经、纬线均为平行直线而且经纬线之间互相正交。
地球上的等角航线在墨卡托投影地图上均表现为直线。
图上能正确表示航线与子午线的夹角。例如，航线在地球上与子午线的夹角为 $15^{\circ}$ ，图上也表示为 $15^{\circ}$ 。
比例尺随纬度而变化，在一张地图上纬度不同，同一距离的比例尺也不同，不能比较面积。因而在量测距离时不能使用同一分划尺（但以海里为单位量测时，因为1海里即为纬度，可使用所量地域附近的分划尺）。

墨卡托投影是按等角条件修改透视圆筒投影而得到的投影。

等角(也称为保形）是指当地图上任何一点的各方向具有相同的比例，称为局部保形，透视圆筒投影如下图所示：

从墨卡托投影图上可以看出，经线间隔的经度如果相等，则经线是等距平行的直线，纬线也是平行的直线，而且经纬线是相互垂直的。

墨卡托投影有一个特别的特性：所有罗盘等角线，或称斜航线（就是与所经过的所有经线形成相同角度的航线，也称恒向航线）在墨卡托投影下都是直线。这使得在航海领域这个投影非常重要。

墨卡托投影对透视圆筒投影改造点：要使圆筒投影称为等角的性质，必须使由赤道向两极经线逐渐伸长的倍数与经线上各点相应的纬度扩大的倍数相同。

地球赤道的圆周长为2 $\pi$ R（R为赤道半径），而各纬线圈周长为2 $\pi$ R $\cos \Phi$ （ $\Phi$ 表示对应的纬度），因此，墨卡托投影地图上纬线长与地球上实际纬线长的比值为：

投影上的纬线长/地球上实际纬线长 = $2{\pi}R/2{\pi}Rcos\phi$ = $sec\phi$

既然各纬度的纬线扩大 $\sec \Phi$ （正割）倍，为了保持等角，各纬线通过处的经线也要相应的扩大 $\sec \Phi$ 倍。

这时，经线方向上的长度比才能与纬线方向上的长度比相等。

注意：投影地图上经纬线的伸长与纬度的正割成比例变化，随纬度增高极具拉伸，到极点成为无穷大；面积的扩大更为明显，在 $60^{\circ}$ 的地方面积要扩大四倍（因为 $\sec 60^{\circ}$ = 2，面积比是长度比的两倍，所以是四倍）。如下图所示，地理上等半径圆在高纬度面积明显扩大。(一般都只截取了85度纬度以内的部分，否则会是一个无穷高的图；与此同时由于取到85度时仍然会让地图变得很大，互联网地图用的墨卡托投影在计算上也会做一些处理，让地图没那么大~~）

Web墨卡托投影

最后还有一点，那就是程序员也特别喜欢，Google基于墨卡托投影设计了一个Web墨卡托投影，使得全世界可以在一个正方形里面刚好放下，极大简化了Web地图的涉及和计算的效率。

3、墨卡托投影的缺点

墨卡托投影的地图最大的缺点就是和现实差别太大，变形非常严重。

在墨卡托投影的地图上，变形最严重的就是格陵兰岛了。地图上非洲的大小和格陵兰岛差不多大，但如果计算一下的话：非洲面积约是3020万平方千米，格陵兰岛面积约是217万平方千米，也就是说，非洲的面积是格陵兰岛的14倍！

再比如墨卡托投影的地图上的加拿大，看起来是个瘦瘦的长方形，但实际上，加拿大是个类似正方形的形状。它不仅扭曲得不像样，而且还被放大了好几倍。在地图上加拿大的面积大约是美国的三倍，实际上加拿大比美国大不了多少。

引用来源：墨卡托投影_不睡觉的怪叔叔的博客-CSDN博客_墨卡托投影

三、瓦片地图

引用来源：瓦片地图原理 - 简书

1、GIS介绍

地理信息系统（Geographic Information System或 Geo－Information system，GIS）有时又称为“地学信息系统”。它是一种特定的十分重要的空间信息系统。它是在计算机硬、软件系统支持下，对整个或部分地球表层（包括大气层）空间中的有关地理分布数据进行采集、储存、管理、运算、分析、显示和描述的技术系统。和维护性。

普通地图分类

按照比例尺大小

大比例尺地形图：1：5千—1：2.5万比例尺地形图（1cm：5k cm — 1cm：25k cm）
中比例尺地形图：1：5万—1：25 万比例尺地形图
小比例尺地形图：1：50万-1：100万比例尺地形图
我国称1：1万、1：2.5万、1：5万、1：10万、1：25万、1：50万、1：100万七种比例尺普通地图为国家基本比例尺地形图

备注：地图比例尺越大误差越小，图上测量精度越高。大比例尺地图，内容详细，几何精度高，可用于图上测量。小比例尺地图，内容概括性强，不宜于进行图上测量

实现方式
瓦片地图：由一张一张的正方形小图片拼接成的地图。
矢量地图：根据矢量数据，绘制出来的地图。

2、瓦片地图原理

瓦片
指将一定范围内的地图按照一定的尺寸和格式，按缩放级别或者比例尺，切成若干行和列的正方形栅格图片，切片后的正方形栅格图片被形象的称为瓦片（Tile）。

瓦片.png

地图瓦片.png

层次模型
瓦片地图金字塔模型是一种多分辨率层次模型，从瓦片金字塔的底层到顶层，分辨率越来越低，但表示的地理范围不变。

地图瓦片金字塔模型.png

地图瓦片层级关系.png

地图坐标系

经纬度坐标系
像素坐标系
像素坐标系，也可以成为屏幕坐标系，像素坐标系和地图的经纬度坐标系存在对应关系，屏幕上的每一个像素都对应一个经纬度点位置。不同缩放级别下，像素坐标系和经纬度坐标系的对应关系是不同的。
投影坐标系
投影坐标系在GIS系统中非常的重要，作用是将经纬度坐标系和屏幕的像素坐标系关联起来。理解投影坐标系也是GIS系统研究中的一个重要环节。

地图投影
地图投影是利用一定数学法则把地球表面的经、纬线转换到平面上的理论和方法。由于地球是一个赤道略宽两极略扁的不规则的梨形球体，故其表面是一个不可展平的曲面，所以运用任何数学方法进行这种转换都会产生误差和变形，为按照不同的需求缩小误差，就产生了各种投影方法。

常见的投影类型

高斯-克吕格投影
斜轴等面积方位投影
双标准纬线等角圆锥投影
等差分纬线多圆锥投影
正轴方位投影
墨卡托投影（正轴等角圆柱投影）

墨卡托投影
墨卡托投影，是正轴等角圆柱投影。由荷兰地图学家墨卡托(G.Mercator)于1569年创立。假想一个与地轴方向一致的圆柱切或割于地球，按等角条件，将经纬网投影到圆柱面上，将圆柱面展为平面后，即得本投影。等角条件是使地球面上微分区域内两个方向的夹角投影到平面以后，保持角度不变的条件。
百度地图和Google Maps使用的投影方法都是墨卡托投影。

墨卡托投影结果区域图.png

投影结果

图中红色矩形区域是真实投影结果。
图中蓝色正方形区域，是数学计算区域。
正方形的边长是赤道周长。
纵轴方向，纬度的变化是非线性的。
横轴方向，经度的变化是线性的。
横轴和纵轴方向，投影区域和都比真实区域大，处于拉伸状态。

墨卡托投影坐标算法

手机与地图的位置关系

手机与地图的位置关系.png

四、瓦片地图原理---续

引用来源：https://segmentfault.com/a/1190000011276788

1、经纬度

地球是一个椭球，Datum是一组用于描述这个椭球的数据集合。最常用的一个Datum是WGS84（World Geodetic System 1984），它的主要参数有：

坐标系的原点是地球质心（center of mass）；
子午线（meridian），即零度经线，位于格林威治子午线Royal Observatory所在纬度往东102.5米所对应的的经线圈；
椭球截面长轴为a=6378137米；
椭圆截面短轴为b=6356752.3142米，可选参数；
扁平比例（flattening）f=(a−b)/a=1/298.257223563；
geoid，海平面，用于定义高度，本文从略。

通过以上参数设定，我们才能对地球上的任意一个位置用经度、纬度、高度三个变量进行描述。所以当我们获取一组经纬度信息时，首先要弄明白这组信息对应的Datum。

WGS84 Datum的信息可以用下图进行概括：

2、投影

地图是显示在平面上的，因此需要将球面坐标转换为平面坐标，这个转换过程称为投影。最常见的投影是墨卡托（Mercator）投影，它具有等角性质，即球体上的两点之间的角度方位与平面上的两点之间的角度方位保持不变，因此特别适合用于导航。

Web墨卡托投影（又称球体墨卡托投影）是墨卡托投影的变种，它接收的输入是Datum为WGS84的经纬度，但在投影时不再把地球当做椭球而当做半径为6378137米的标准球体，以简化计算。其计算公式推导请参考下图：

备注：这个计算我暂时还没弄懂，看以后能不能补充！

Web墨卡托投影有两个相关的投影标准，经常搞混：

EPSG4326：Web墨卡托投影后的平面地图，但仍然使用WGS84的经度、纬度表示坐标；
EPSG3857：Web墨卡托投影后的平面地图，坐标单位为米。

3、瓦片

经过Web墨卡托投影后，地图就变为平面的一张地图。考虑到有时候我们需要看宏观的地图信息（如世界地图里每个国家的国界），有时候又要看很微观的地图信息（如导航时道路的路况信息）。为此，我们对这张地图进行等级切分。在最高级（zoom=0），需要的信息最少，只需保留最重要的宏观信息，因此用一张256x256像素的图片表示即可；在下一级（zoom=1），信息量变多，用一张512x512像素的图片表示；以此类推，级别越低的像素越高，下一级的像素是当前级的4倍。这样从最高层级往下到最低层级就形成了一个金字塔坐标体系。

对每张图片，我们将其切分为256x256的图片，称为瓦片（Tile）。这样，在最高级（zoom=0）时，只有一个瓦片；在下一级（zoom=1）时有4个瓦片；在下一级（zoom=2）时有16个瓦片，以此类推。

4、瓦片编号

瓦片生成后，就是一堆图片。怎么对这堆图片进行编号，是目前主流互联网地图商分歧最大的地方。总结起来分为四个流派：

谷歌XYZ：Z表示缩放层级，Z=zoom；XY的原点在左上角，X从左向右，Y从上向下。
TMS：开源产品的标准，Z的定义与谷歌相同；XY的原点在左下角，X从左向右，Y从下向上。
QuadTree：微软Bing地图使用的编码规范，Z的定义与谷歌相同，同一层级的瓦片不用XY两个维度表示，而只用一个整数表示，该整数服从四叉树编码规则。
百度XYZ：Z从1开始，在最高级就把地图分为四块瓦片；XY的原点在经度为0纬度位0的位置，X从左向右，Y从下向上。

下图显示了前三个流派在zoom=1层级上的瓦片编号结果：

下表总结了中国主要地图商的瓦片编号流派，点击每个链接就可以获得一个对应编号的瓦片地图：

地图商	瓦片编码	图层	链接
高德地图	谷歌XYZ	道路	`http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x=105&y=48&z=7`
高德地图	谷歌XYZ	卫星	`http://webst04.is.autonavi.com/appmaptile?style=6&x=843&y=388&z=10`
谷歌地图	谷歌XYZ	道路	`http://mt2.google.cn/vt/lyrs=m&hl=zh-CN&gl=cn&x=105&y=48&z=7`
谷歌地图	谷歌XYZ	卫星	`http://mt2.google.cn/vt/lyrs=s&hl=zh-CN&gl=cn&x=105&y=48&z=7`
谷歌地图	谷歌XYZ	地形	`http://mt0.google.cn/vt/lyrs=t&hl=zh-CN&gl=cn&x=420&y=193&z=9`
OpenStreetMap	谷歌XYZ	道路	`http://a.tile.openstreetmap.org/7/105/48.png`
腾讯地图	TMS	道路	`http://rt1.map.gtimg.com/realtimerender?z=7&x=105&y=79&type=vector&style=0`
Bing地图	QuadTree	道路	`http://r1.tiles.ditu.live.com/tiles/r1321001.png?g=100&mkt=zh-cn`
百度地图	百度XYZ	道路	`http://online4.map.bdimg.com/tile/?qt=tile&x=98&y=36&z=9&;styles=pl&scaler=1&udt=20170406`
百度地图	百度XYZ	交通	`http://its.map.baidu.com:8002/traffic/TrafficTileService?level=19&x=99052&y=20189&time=1373790856265&label=web2D&;v=017`

补充一个新的瓦片地图资源：http://map.geoq.cn/，里面有彩色、灰色、暖色和夜色版的中国地图，坐标为GCJ02。举例：http://map.geoq.cn/ArcGIS/rest/services/ChinaOnlineStreetWarm/MapServer/tile/{z}/{y}/{x}

5、关于中国的经纬度

国内的经纬度有三套系统：

WGS84：为一种大地坐标系，也是目前广泛使用的GPS全球卫星定位系统使用的坐标系。
GCJ02：又称火星坐标系，是由中国国家测绘局制订的地理信息系统的坐标系统。由WGS84坐标系经加密后的坐标系。
BD09：为百度坐标系，在GCJ02坐标系基础上再次加密。其中bd09ll表示百度经纬度坐标，bd09mc表示百度墨卡托米制坐标。

使用OpenStreetMap的坐标为WGS84；使用高德地图、腾讯地图的坐标为GCJ02；使用百度地图的坐标为BD09；谷歌地图和Bing地图的中国部分采用了高德地图的数据，所以坐标为GCJ02。

WGS84的坐标转化为GCJ02的坐标是单向的，即WGS84的坐标能够准确地变换为GCJ02坐标；但GCJ02坐标转换为WGS84时会存在精度损失。

GCJ02的坐标和BD09的坐标转换是双向的。

LuossD

关注

22
点赞
踩
36

收藏

觉得还不错? 一键收藏
1
评论
墨卡托投影与瓦片地图

一、开胃小知识先写几个小知识点：1）既然地球是圆的，但是为什么地面却是平的呢？这主要是因为人和地球比起来，相差悬殊；人眼看到的地球表面区域弧度非常小，人类察觉不到，从而产生地面平坦的错觉，我们因此会认为“地面是平面”了。2）纬度是指某点与地球球心的连线和地球赤道面所成的线面角，其数值在0至90度之间。经纬度是经度与纬度的合称组成一个坐标系统。又称为地理坐标系统。东经正数，西经为负数；北纬为正数，南纬为负数。3）东西半球分界线不是本初子午线；而是东经160度和与其相对的西经20度所组成的经线圈
复制链接

扫一扫

专栏目录