メルカトル図法の 1km は実際どれくらいの距離なのか

今回はメルカトル図法の距離に関するお話です。前回は緯度・経度の 1度が場所によってどれくらいの距離なのかを測りましたが、今回はメルカトル図法で示される距離の実際を検証します。

はじめに

メルカトル図法は正角図法で形は正確だが距離や面積は正確ではない

メルカトル図法は上記の特性があります。前回の記事でも触れましたが、Google Maps は 2018年8月からメルカトル図法以外に外射図法でも表示できるようになりました。これで地図の見た目は宇宙からの見た目と同じになった訳ですが、ArcGIS Online や 地理院地図など他の Web マップ サービスは依然としてメルカトル図法（Webメルカトル）が多く使用されています。

メルカトル図法を使用して、例えばこのような主題図を作成したとします。この地図には大きな問題があります。

メルカトル図法にスケールバーを配置した例

メルカトル図法では（赤道上を除いて）距離は正しくないので、スケールバーを配置するのは不適切です。スケールバーが配置されると一般的に地図の読み手は図上どこでもスケールバーと同じ距離になると理解します。上記の地図で示されるスケールバーが正確なのは赤道上のみです。ArcGIS Online のマップ ビューアーや Google Maps でもスケールバーがありますが、これは現在表示されている地図のある場所（おそらく中心）における縮尺に基づいて動的に変化します。大縮尺で表示すれば誤差の範囲内ですが、小縮尺で表示すると正確な表示にはなりません。

メルカトル図法で示すスケールバーは正確でないと説明しましたが、どれくらい正確でないのか気になります。ということで、メルカトル図法によって示された 2点間の距離と測地線距離がどれくらい異なるのか計測してみます。

前提条件

地球の大きさには様々な定義がありますが、今回も赤道半径と極半径が異なる WGS84 楕円体を使用します。

WGS84楕円体
赤道半径：6,378,137m、扁平率の逆数：298.257223563
GCS_WGS_1984 (4326)

回転楕円体における計算式で長さを求めることもできますが、難しい計算は ArcGIS にやらせようということで今回も arcpy という API を使用して計算しました。

経度0度上で北緯0度～ 90度まで1度ずつの間隔で点を置き、その緯度経度の座標をメルカトル図法の座標に換算します。換算後の座標から東西南北に 500m 延伸させ、「東西・南北で 1km の線分」を作成します。作成した線分の測地線距離を測定します。

メルカトル図法上 1000m の測地線距離と拡大率

図法上 1000m の測地線距離を求めた結果です。メルカトル図法は緯度によって実際の距離が異なり、高緯度になると顕著に拡大率が大きくなります。89度以上は計算できないので除外しました。

Spatial Reference: World_Mercator (54004)
GCS: GCS_WGS_1984 (4326)

メルカトル図法の誤解

メルカトル図法の緯度に対する拡大率は余弦からこのように求められます。

緯度(φ) 上の拡大率 = 1/cos(φ)

0°：1
60°：2
75°：2.9238
80°：5.7588
85°：11.4737

仮に「メルカトル図法は地球に円筒を被せて中心から光を当てた影」でできるとすると、そんなにゆるやかな拡大率はとはなりません。これについては以下の書籍で解説されています。

地図投影法

政春尋志

4,200円(12/03 18:38時点)

Amazon楽天市場Yahoo

Amazonの情報を掲載しています

『地図投影法』P45 column4 地図投影法についてよくある間違い
地図投影法に関して広まっている誤解の筆頭は、メルカトル図法の原理に関するものであろう。球に円筒を巻きつけて球の中心から投影すればメルカトル図法が得られるという解説や図解は、学生用参考書や一般向けの地図に関する書籍によく見られるものであるが、間違いである。こうして得られる投影は心射円筒図法であってメルカトル図法ではない。

私もこの書籍を読んではじめて誤りに気づきました。メルカトル図法の計算式はこちらに記載があります。

まとめ

計測した結果、緯度25度くらいまでは1割程度、36度あたりで2割の拡大率で、60度になるとほぼ2倍となり、高緯度になるほど顕著に拡大率が大きくなりました。そのため、メルカトル図法は距離や面積の測定には不向きだということが分かります。また、同一地点で東西と南北の拡大比率はほとんど同じです。このことから、メルカトル図法は歪みが少ない（形状は正確）ということが分かります。

2年ほど前「意外と日本って大きい。」というタイトルで The True Size Of... サイトを紹介し、ネットで盛り上がりましたがあのタイトルはナンセンスです。メルカトル図法は北にいくほど拡大されるのだから日本が大きくなるのは当たり前です。比較するなら赤道に移動させ「ロシアって意外と小さい。」というタイトルにするべきでしょう（ロシアは面積世界一なので赤道に移動しても大きいんだけど）。

ソースコード

Python ウィンドウにそのまま貼り付けてご利用ください。ArcMap 10.4、AcGIS Pro 1.3、ArcGIS Pro 2.2 で動作確認しております。arcpy.SpatialReference の引数である EPSG 番号を変更して一般的に使用される球体としての地球で計算することもできます。

#coding:UTF-8

# メルカトル図法の指定距離における測地線距離
# 経度0度、北緯0度～90度を中心に東西・南北の指定距離を
# 測地線距離に換算して比較

# ArcGIS Pro 1.3,ArcGIS Pro 2.2, ArcMap 10.4.1 で動作確認

import arcpy
import copy

# WGS84楕円体
# 赤道半径：6,378,137m、扁平率の逆数：298.257223563
# GCS_WGS_1984 (4326)
gcs = arcpy.SpatialReference(4326)
#gcs = arcpy.SpatialReference(4035)	# 地球半径 Authalic Shpere GCS_Shpere GCS_Sphere (4035)

# メルカトル図法
pcs = arcpy.SpatialReference(54004)

# 測定係数
ratio = 1000 # (m)

print("Spatial Reference: {0} ({1})".format(pcs.PCSName, pcs.PCSCode))
print("              GCS: {0} ({1})\n".format(gcs.GCSName, gcs.GCSCode))
print("メルカトル図法上 {}m の測地線距離\n".format(ratio))
print("{}\t{}\t{}\t{}\t{}".format("緯度", "東西距離 (m)".ljust(10), "南北距離 (m)".ljust(10), "東西誇張率 (%)".ljust(10), "南北誇張率 (%)".ljust(10)))

for lat in range(88, -1, -1):
    # 中心点を投影変換
    point = arcpy.Point(0, lat)
    point_gcs = arcpy.PointGeometry(point, gcs)
    point_pcs = point_gcs.projectAs(pcs)

    # 東西長
    point_width_from = copy.copy(point_pcs.trueCentroid)
    point_width_to = copy.copy(point_pcs.trueCentroid)
    point_width_to.X += 1 * ratio
    polyline_width = arcpy.Polyline(arcpy.Array([point_width_from, point_width_to]), pcs)
    length_width = polyline_width.getLength("GEODESIC")
    
    # 南北長
    point_height_from = copy.copy(point_pcs.trueCentroid)
    point_height_from.Y -= 0.5 * ratio
    point_height_to = copy.copy(point_pcs.trueCentroid)
    point_height_to.Y += 0.5 * ratio
    polyline_height = arcpy.Polyline(arcpy.Array([point_height_from, point_height_to]), pcs)
    length_height = polyline_height.getLength("GEODESIC")  
    
    print("{:0=2}\t{:0=13.7f}\t{:0=13.7f}\t{:0=13.7f}\t{:0=13.7f}".format(lat, length_width, length_height, ratio / length_width * 100, ratio / length_height * 100))