Web麥卡托投影
Web麥卡托投影(英語:Web Mercator)[1]是麥卡托投影的一種變體,被Web地圖應用業界普遍採納。2005年,隨着Google地圖的採用,這一標準開始嶄露頭角。[2] 幾乎所有主要的線上地圖提供商都使用這一標準,包括谷歌地圖、Mapbox[3]、Bing地圖、OpenStreetMap、MapQuest、Esri等等。[4]其正式的EPSG識別碼是EPSG:3857,不過歷史上也曾使用過其他識別碼。
性質
[編輯]Web麥卡托是麥卡托投影地圖的一個輕微變體,主要用於基於Web的地圖程式。對於小比例尺地圖,它與標準的麥卡托用的公式一樣。但Web麥卡托在所有比例尺下都使用球面公式,但大比例尺的麥卡托地圖通常使用投影的橢球面形式。這種差異在全球比例尺下是察覺不到的,但會導致局部地區的地圖稍微偏離同一比例尺的真正的橢球面麥卡托地圖。 雖然Web麥卡托的公式是麥卡托的球面形式,但地理坐標必須得以WGS 84橢球面基準獲得。此差異會導致投影略微不符合正形投影。人們普遍對Web麥卡托與標準麥卡托的不同認識不足,導致了相當多的混亂和誤用。[4]:87錯誤轉換造成的誤差在地面上可以達到40 km。[5][6] 出於所有這些原因,美國國防部通過美國國家地理空間情報局宣佈官方用途不可以使用這種地圖投影。[7]
公式
[編輯]Web麥卡托的公式與標準球面麥卡托的公式基本相同,但是在應用縮放之前,將「世界坐標」調整為使得左上角為(0, 0),右下角為(256, 256):[8]
其中 λ 是用弧度表示的經度,而 φ 是用弧度表示的大地緯度。
因為麥卡托將極點投影在無窮遠處,所以使用Web麥卡托投影的地圖無法顯示極點。Google地圖之類的服務會在南北緯85.051129°切斷覆蓋。這不會對街道地圖造成什麼影響,因為這些服務的主要目的不在於顯示完整的地球。85.051129°這個值是能讓整個投影地圖變成正方形的緯度,在 y = 0 的情況下計算為 φ 就會得到:
球面與橢球面的混合
[編輯]這個投影既不是嚴格的橢球面,也不是嚴格的球面,EPSG的定義說這個投影「使用橢球坐標系的球面演化」。[9] 它是建立在地球表面的WGS84橢球面模型定義的地理坐標上的,但在投影的時候卻仿佛坐標是定義在球面上的。[6] 這種做法對於小比例尺的地圖(例如整個世界)是沒有爭議的,但是在大比例尺地圖(例如城市或省份)中卻很少有先例。[10]
優點和缺點
[編輯]Web麥卡托是一個球面麥卡托投影,因此它具有與球面麥卡托相同的屬性: 不管在哪,北都在上方;子午線都是等距的垂直線,角度在局部是正確的(假設是球面坐標),而且面積會隨着離赤道越來越遠而膨脹,以至於極地區域被嚴重誇大。橢球面麥卡托具有同樣的屬性,但是它將地球模擬為一個橢球體。
然而,與橢球面麥卡托不同的是,Web麥卡托並不十分符合正形投影。這意味着地面上線條之間的角度,繪製到地圖中角度會發生變化,儘管它們不會偏離到肉眼可以察覺的程度。由於Web麥卡托指定了在WGS 84橢球面模型上給出的測量坐標,因此會有偏離。通過將測量坐標投影到球面上,角度關係會發生輕微的變化,就像在一個球面上測量一樣。這是標準球面麥卡托投影法的標準做法,但與Web麥卡托不用的是,球面麥卡托通常不用於局部區域的地圖,例如街道地圖,因此繪圖所需的位置精度通常小於使用球面公式引起的角偏差。Web麥卡托的好處在於,球面形式的計算比橢球面形式的計算簡單得多,因此只需要一小部分計算資源。[11]
辨識碼
[編輯]由於EPSG註冊通過遲緩,Web麥卡托有了好幾個不同的空間參考系辨識碼(SRID),包括EPSG:900913、EPSG:3785和EPSG:3857。
EPSG:900913
[編輯]該投影坐標參考系最初缺少一個官方的SRID,OGP的地理資訊學委員會(又名EPSG)的大地測量小組委員會拒絕提供一個SRID,並稱「我們已經審查了微軟、谷歌等公司使用的坐標參考系統,並認為它在技術上有缺陷。我們不會為了收錄這些不適當的大地測量和製圖,而降低EPSG數據集的價值。」[12] 非官方的代碼「EPSG:900913」(Google轉寫為數字)開始被使用。它最初由Christopher Schmidt在他的Technical Ramblings網誌[13]中定義,並在OpenLayers 2[14]中進行了編碼,從技術上講,它將使OpenLayers成為SRID權威。
EPSG:3785
[編輯]2008年,EPSG提供了官方辨識碼EPSG:3785,其官方名稱為「Popular Visualisation CRS / Mercator」,但指出「它不是一個正式的大地測量系統」。[6] 該定義使用了地球的球面(而非橢球面)模型。
EPSG:3857
[編輯]當年晚些時候,EPSG提供了一個更新的辨識碼EPSG:3857,其正式名稱為「WGS 84 / Pseudo-Mercator」。[6] 該定義轉換為使用WGS84橢球面(EPSG:4326),而不是使用球面。
雖然這一投影與Google密切相關,但在EPSG的標準里,把微軟列為了「資訊來源」。[15]
其他辨識碼
[編輯]其他使用過的辨識碼有ESRI:102113、ESRI:102100和OSGEO:41001。[15][13]
ESRI:102113對應於EPSG:3785,而ESRI:102100對應於EPSG:3857。[16]
WKT定義
[編輯]該投影覆蓋了經度從−180°到180°,南北緯85.05°之間的地球。使用坐標參考系的知名文字表示(well-known text representation of coordinate reference systems,簡稱WKT),EPSG:3857定義如下:[15]
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"]]
參考文獻
[編輯]- ^ WGS 84 and the Web Mercator Projection NGA Office of Geomatics (PDF). National Geospatial Intelligence Agency. 2014-05-16 [2014-08-06]. (原始內容 (PDF)存檔於2017-03-29).
- ^ Google Maps & Earth Help Forum. [2017-08-29].
- ^ Our Map Data. MapBox. [2018-06-20]. (原始內容存檔於2018-12-29).
Mapbox supports the popular Web Mercator projection, and currently does not support any other projections for display.
- ^ 4.0 4.1 Battersby, Sarah E.; Finn, Michael P.; Usery, E. Lynn; Yamamoto, Kristina H. Implications of Web Mercator and Its Use in Online Mapping (PDF). Cartographica. 2014, 49 (2): 92 [2020-12-12]. doi:10.3138/carto.49.2.2313. (原始內容存檔 (PDF)於2015-02-15).
- ^ NGA: (U) NGA Advisory Notice on "Web Mercator" (UNCLASSIFIED). earth-info.nga.mil. [2018-06-07]. (原始內容存檔於2019-05-17).
- ^ 6.0 6.1 6.2 6.3 The Google Maps / Bing Maps Spherical Mercator Projection. Alastair Aitchison. [2014-10-04]. (原始內容存檔於2014-10-06).
- ^ NGA: (U) NGA Advisory Notice on "Web Mercator" (UNCLASSIFIED). [2014-10-04]. (原始內容存檔於2016-03-04).
- ^ OpenStreetMap Slippy map tilenames. openstreetmap.org. [2018-11-14]. (原始內容存檔於2022-04-21).
- ^ GmbH (https://www.klokantech.com/), Klokan Technologies. WGS 84 / Pseudo-Mercator - Spherical Mercator, Google Maps, OpenStreetMap, Bing, ArcGIS, ESRI - EPSG:3857. epsg.io. [2020-12-12]. (原始內容存檔於2022-04-17) (英語).
- ^ Strebe, Daniel "daan". A bevy of area-preserving transforms for map projection designers. Cartography and Geographic Information Science. 2019, 46 (3). doi:10.1080/15230406.2018.1452632.
- ^ The Earth is Not Round! Utah, NAD83 and WebMercator Projections. Utah GIS Portal. [2020-06-20]. (原始內容存檔於2021-12-07) (英語).
- ^ dotMorten | Spherical/Web Mercator: EPSG code 3785. www.sharpgis.net. [2020-12-12]. (原始內容存檔於2021-12-09).
- ^ 13.0 13.1 Google Projection: 900913. 2007-08-07 [2020-12-12]. (原始內容存檔於2021-06-20).
- ^ Spherical Mercator — OpenLayers. docs.openlayers.org. [2020-12-12]. (原始內容存檔於2021-04-18).
- ^ 15.0 15.1 15.2 GmbH (https://www.klokantech.com/), Klokan Technologies. WGS 84 / Pseudo-Mercator - Spherical Mercator, Google Maps, OpenStreetMap, Bing, ArcGIS, ESRI - EPSG:3857. epsg.io. [2020-12-12]. (原始內容存檔於2022-05-06).
- ^ 102100 vs. 102113 is mixed up?. GeoNet, The Esri Community. 2012-10-03 [2020-12-12]. (原始內容存檔於2021-04-17).