查看“︁Web墨卡托投影”︁的源代码

{{NoteTA|G1=Geography|G2=IT}}
[[File:OpenStreetMap homepage 2018 en.png|thumb|400px|2018年OpenStreetMap的主页，OpenStreetMap的标准风格，像大多数网络地图一样，使用Web墨卡托投影]]
'''Web墨卡托投影'''（{{lang-en|Web Mercator}}）<ref>{{Cite web|url = http://earth-info.nga.mil/GandG/wgs84/web_mercator/%28U%29%20NGA_SIG_0011_1.0.0_WEBMERC.pdf|title = WGS 84 and the Web Mercator Projection NGA Office of Geomatics|date = 2014-05-16|accessdate = 2014-08-06|website = |publisher = National Geospatial Intelligence Agency|last = |first = |archive-date = 2017-03-29|archive-url = https://web.archive.org/web/20170329052456/http://earth-info.nga.mil/GandG/wgs84/web_mercator/(U)%20NGA_SIG_0011_1.0.0_WEBMERC.pdf}}</ref>是[[麥卡托投影法|墨卡托投影]]的一种变体，被[[电子地图服务|Web地图]]应用业界[[業界標準|普遍采纳]]。2005年，随着[[Google地图]]的采用，这一标准开始崭露头角。<ref>{{cite web|url=https://productforums.google.com/d/msg/maps/A2ygEJ5eG-o/KbZr_B0h2hkJ|title=Google Maps & Earth Help Forum|accessdate=2017-08-29}}</ref> 几乎所有主要的在线地图提供商都使用这一标准，包括谷歌地图、[[Mapbox]]<ref name="MapBox">{{cite web
 | url       = https://www.mapbox.com/help/how-mapbox-data-works/
 | title       = Our Map Data
 | publisher       = MapBox
 | access-date       = 2018-06-20
 | quote       = Mapbox supports the popular Web Mercator projection, and currently does not support any other projections for display.
 | archive-date       = 2018-12-29
 | archive-url       = https://web.archive.org/web/20181229025832/https://www.mapbox.com/help/how-mapbox-data-works/
 }}</ref>、[[Bing地图]]、[[開放街圖|OpenStreetMap]]、{{le|MapQuest}}、[[美國環境系統研究所公司|Esri]]等等。<ref name="Battersby">{{cite journal
 |journal      = Cartographica
 |title        = Implications of Web Mercator and Its Use in Online Mapping
 |last1        = Battersby
 |first1       = Sarah E.
 |last2        = Finn
 |first2       = Michael P.
 |last3        = Usery
 |first3       = E. Lynn
 |last4        = Yamamoto
 |first4       = Kristina H.
 |year         = 2014
 |volume       = 49
 |issue        = 2
 |page         = 92
 |doi          = 10.3138/carto.49.2.2313
 |url          = http://cegis.usgs.gov/projection/pdf/Battersby_Implications%20of%20Web%20Mecator%20and%20Its%20Use%20in%20Online%20Mapping.pdf
 |access-date  = 2020-12-12
 |archive-date = 2015-02-15
 |archive-url  = https://web.archive.org/web/20150215113220/http://cegis.usgs.gov/projection/pdf/Battersby_Implications%20of%20Web%20Mecator%20and%20Its%20Use%20in%20Online%20Mapping.pdf
 |dead-url     = no
}}</ref>其正式的EPSG标识符是EPSG:3857，不过历史上也曾使用过其他标识符。

== 性质 ==
Web墨卡托是墨卡托投影地图的一个轻微变体，主要用于基于Web的地图程序。对于[[比例尺 (地圖)|小比例尺地图]]，它与标准的墨卡托用的公式一样。但Web墨卡托在所有比例尺下都使用球面公式，但大比例尺的墨卡托地图通常使用投影的椭球面形式。这种差异在全球比例尺下是察觉不到的，但会导致局部地区的地图稍微偏离同一比例尺的真正的椭球面墨卡托地图。
虽然Web墨卡托的公式是墨卡托的球面形式，但地理坐标必须得以[[World Geodetic System|WGS 84]]椭球面基准获得。此差异会导致投影略微不符合{{le|正形投影|Conformal map projection}}。人们普遍对Web墨卡托与标准墨卡托的不同认识不足，导致了相当多的混乱和误用。<ref name="Battersby"/>{{rp|87}}错误转换造成的误差在地面上可以达到40&nbsp;km。<ref>{{Cite web|url=http://earth-info.nga.mil/GandG/wgs84/web_mercator/|title=NGA: (U) NGA Advisory Notice on "Web Mercator" (UNCLASSIFIED)|website=earth-info.nga.mil|access-date=2018-06-07|archive-date=2019-05-17|archive-url=https://web.archive.org/web/20190517111709/http://earth-info.nga.mil/GandG/wgs84/web_mercator/}}</ref><ref name="al">{{cite web|url=http://alastaira.wordpress.com/2011/01/23/the-google-maps-bing-maps-spherical-mercator-projection/|title=The Google Maps / Bing Maps Spherical Mercator Projection|work=Alastair Aitchison|accessdate=2014-10-04|archive-date=2014-10-06|archive-url=https://web.archive.org/web/20141006110906/http://alastaira.wordpress.com/2011/01/23/the-google-maps-bing-maps-spherical-mercator-projection/|dead-url=no}}</ref>
出于所有这些原因，[[美国国防部]]通过[[美國國家地理空間情報局]]宣布官方用途不可以使用这种地图投影。<ref>{{cite web|url=http://earth-info.nga.mil/GandG/wgs84/web_mercator/index.html|title=NGA: (U) NGA Advisory Notice on "Web Mercator" (UNCLASSIFIED)|publisher=|accessdate=2014-10-04|archive-date=2016-03-04|archive-url=https://web.archive.org/web/20160304080752/http://earth-info.nga.mil/GandG/wgs84/web_mercator/index.html}}</ref>

=== 公式 ===
Web墨卡托的公式与标准球面墨卡托的公式基本相同，但是在应用缩放之前，将“世界坐标”调整为使得左上角为{{nowrap|(0, 0)}}，右下角为{{nowrap|(256, 256)}}：<ref>{{cite web |url=https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames |title=OpenStreetMap Slippy map tilenames |publisher=openstreetmap.org |accessdate=2018-11-14 |archive-date=2022-04-21 |archive-url=https://web.archive.org/web/20220421120137/https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames }}</ref>
:<math>\begin{align}
  x &= \left\lfloor\frac{256}{2\pi} 2^{\text{zoom level}} (\lambda + \pi)\right\rfloor \text{ pixels}\\[5pt]
  y &= \left\lfloor\frac{256}{2\pi} 2^{\text{zoom level}} \left(\pi - \ln \left[\tan \left(\frac{\pi}{4} + \frac{\varphi}{2} \right) \right]\right)\right\rfloor \text{ pixels}
\end{align}</math>
其中 ''λ'' 是用弧度表示的[[经度]]，而 ''φ'' 是用弧度表示的[[纬度#大地纬度和地心纬度|大地纬度]]。

因为墨卡托将极点投影在无穷远处，所以使用Web墨卡托投影的地图无法显示极点。Google地图之类的服务会在南北纬85.051129°切断覆盖。这不会对街道地图造成什么影响，因为这些服务的主要目的不在于显示完整的地球。85.051129°这个值是能让整个投影地图变成正方形的纬度，在 {{nowrap|1= ''y'' = 0}} 的情况下计算为 ''φ'' 就会得到：
:<math>\begin{align}
\varphi_{\text{max}} = \left[2\arctan(e^{\pi}) - \frac{\pi}{2}\right]
\end{align}</math>

=== 球面与椭球面的混合 ===
这个投影既不是严格的椭球面，也不是严格的球面，EPSG的定义说这个投影“使用椭球坐标系的球面演化”。<ref>{{cite web |last1=GmbH (https://www.klokantech.com/) |first1=Klokan Technologies |title=WGS 84 / Pseudo-Mercator - Spherical Mercator, Google Maps, OpenStreetMap, Bing, ArcGIS, ESRI - EPSG:3857 |url=https://epsg.io/3857 |website=epsg.io |language=en |access-date=2020-12-12 |archive-date=2022-04-17 |archive-url=https://web.archive.org/web/20220417102312/https://epsg.io/3857 }}</ref> 它是建立在地球表面的[[世界大地测量系统|WGS84]]椭球面模型定义的地理坐标上的，但在投影的时候却仿佛坐标是定义在球面上的。<ref name="al"/> 这种做法对于小比例尺的地图（例如整个世界）是没有争议的，但是在大比例尺地图（例如城市或省份）中却很少有先例。<ref name ="Strebe2018">{{cite journal
| last = Strebe
| first = Daniel "daan"
| journal = Cartography and Geographic Information Science
| title = A bevy of area-preserving transforms for map projection designers
| doi = 10.1080/15230406.2018.1452632
| year = 2019
| volume = 46
| issue = 3 }}</ref>

=== 优点和缺点 ===
Web墨卡托是一个球面墨卡托投影，因此它具有与球面墨卡托相同的属性: 不管在哪，北都在上方；子午线都是等距的垂直线，角度在局部是正确的（假设是球面坐标），而且面积会随着离赤道越来越远而膨胀，以至于极地区域被严重夸大。椭球面墨卡托具有同样的属性，但是它将地球模拟为一个椭球体。

然而，与椭球面墨卡托不同的是，Web墨卡托并不十分符合正形投影。这意味着地面上线条之间的角度，绘制到地图中角度会发生变化，尽管它们不会偏离到肉眼可以察觉的程度。由于Web墨卡托指定了在WGS 84椭球面模型上给出的测量坐标，因此会有偏离。通过将测量坐标投影到球面上，角度关系会发生轻微的变化，就像在一个球面上测量一样。这是标准球面墨卡托投影法的标准做法，但与Web墨卡托不用的是，球面墨卡托通常不用于局部区域的地图，例如街道地图，因此绘图所需的位置精度通常小于使用球面公式引起的角偏差。Web墨卡托的好处在于，球面形式的计算比椭球面形式的计算简单得多，因此只需要一小部分计算资源。<ref>{{Cite web|title=The Earth is Not Round! Utah, NAD83 and WebMercator Projections|url=https://gis.utah.gov/nad83-and-webmercator-projections/|access-date=2020-06-20|website=Utah GIS Portal|language=en|archive-date=2021-12-07|archive-url=https://web.archive.org/web/20211207171136/https://gis.utah.gov/nad83-and-webmercator-projections/}}</ref>

== 识别码 ==
由于EPSG注册通过迟缓，Web墨卡托有了好几个不同的[[SRID|空间参考系识别码]]（SRID），包括EPSG:900913、EPSG:3785和EPSG:3857。

=== EPSG:900913 ===
该投影坐标参考系最初缺少一个官方的SRID，OGP的地理信息学委员会（又名EPSG）的大地测量小组委员会拒绝提供一个SRID，并称“我们已经审查了微软、谷歌等公司使用的坐标参考系统，并认为它在技术上有缺陷。我们不会为了收录这些不适当的大地测量和制图，而降低EPSG数据集的价值。”<ref>{{Cite web|url=http://www.sharpgis.net/post/SphericalWeb-Mercator-EPSG-code-3785|title=dotMorten &#124; Spherical/Web Mercator: EPSG code 3785|website=www.sharpgis.net|access-date=2020-12-12|archive-date=2021-12-09|archive-url=https://web.archive.org/web/20211209172652/http://www.sharpgis.net/post/SphericalWeb-Mercator-EPSG-code-3785}}</ref> 非官方的代码“EPSG:900913”（Google转写为数字）开始被使用。它最初由Christopher Schmidt在他的Technical Ramblings博客<ref name="Schmidt">{{Cite web|url=https://crschmidt.net/blog/archives/243/google-projection-900913/|title=Google Projection: 900913|date=2007-08-07|access-date=2020-12-12|archive-date=2021-06-20|archive-url=https://web.archive.org/web/20210620045242/https://crschmidt.net/blog/archives/243/google-projection-900913/}}</ref>中定义，并在OpenLayers 2<ref>{{Cite web|url=http://docs.openlayers.org/library/spherical_mercator.html|title=Spherical Mercator — OpenLayers|website=docs.openlayers.org|access-date=2020-12-12|archive-date=2021-04-18|archive-url=https://web.archive.org/web/20210418033655/http://docs.openlayers.org/library/spherical_mercator.html}}</ref>中进行了编码，从技术上讲，它将使OpenLayers成为SRID权威。

=== EPSG:3785 ===
2008年，EPSG提供了官方识别码EPSG:3785，其官方名称为“Popular Visualisation CRS / Mercator”，但指出“它不是一个正式的[[大地测量学|大地测量]]系统”。<ref name="al"/> 该定义使用了地球的球面（而非椭球面）模型。

=== EPSG:3857 ===
当年晚些时候，EPSG提供了一个更新的识别码EPSG:3857，其正式名称为“WGS 84 / Pseudo-Mercator”。<ref name="al" /> 该定义转换为使用[[世界大地测量系统|WGS84]]椭球面（EPSG:4326），而不是使用球面。

虽然这一投影与Google密切相关，但在EPSG的标准里，把微软列为了“信息来源”。<ref name="epsg.io">{{Cite web|url=http://epsg.io/|title=WGS 84 / Pseudo-Mercator - Spherical Mercator, Google Maps, OpenStreetMap, Bing, ArcGIS, ESRI - EPSG:3857|first=Klokan Technologies|last=GmbH (https://www.klokantech.com/)|website=epsg.io|access-date=2020-12-12|archive-date=2022-05-06|archive-url=https://web.archive.org/web/20220506214359/https://epsg.io/}}</ref>

=== 其他识别码 ===
其他使用过的识别码有ESRI:102113、ESRI:102100和OSGEO:41001。<ref name="epsg.io" /><ref name="Schmidt" />

ESRI:102113对应于EPSG:3785，而ESRI:102100对应于EPSG:3857。<ref>{{Cite web|url=https://community.esri.com/t5/arcgis-rest-api-questions/102100-vs-102113-is-mixed-up/m-p/426380#M2050|title=102100 vs. 102113 is mixed up?|date=2012-10-03|website=GeoNet, The Esri Community|access-date=2020-12-12|archive-date=2021-04-17|archive-url=https://web.archive.org/web/20210417221029/https://community.esri.com/t5/arcgis-rest-api-questions/102100-vs-102113-is-mixed-up/m-p/426380#M2050}}</ref>

== WKT定义 ==
该投影覆盖了经度从−180°到180°，南北纬85.05°之间的地球。使用'''坐标参考系的知名文本表示'''（{{lang|en|well-known text representation of coordinate reference systems}}，简称WKT），EPSG:3857定义如下：<ref name="epsg.io"/>
<pre>
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"],
    AUTHORITY["EPSG","3857"]]
</pre>

==参考文献==
{{reflist|30em}}

[[Category:地理坐标系]]
[[Category:电子地图]]
[[Category:Google地圖]]
[[Category:地圖投影法]]