緯度・経度の 1度はどれくらいの長さがあるのか

地球上の緯度 1度、経度 1度は実際どれくらいの長さがあるのでしょうか。これが今回のテーマです。同じネタは他のブログや Yahoo知恵袋などで扱われているのですが、どれも地球を球として計算しているので、ここでは地球を回転楕円体として計算し、1度ずつ厳密に求めてみます。

経度1度間の距離

私は（現在のところ）年に数回 GIS ソフトウェアの実習講座を行っているのですが、GIS は空間参照を理解しなければ正しい使い方ができません。そのため測地系や投影法の説明をする単元を用意しており、その中で毎回こんな質問します。

はださん

下の図を見てください。
緯度0度、つまり赤道上で東経135度から東経136度までの経線間の距離は約111kmです。日本の関西地方付近である北緯35度上で東経135度から東経136度までの経線間の距離は約91kmです。

問題

はださん

では北緯90度上で経度135度から経度136度までの経線間の距離は何キロメートルでしょうか？

...

はださん

大人になると自分から手を挙げて回答するのは恥ずかしいですよね？では4択で手を挙げてください。
１．赤道上の距離と同じだと思う人？
２．赤道上の距離より短いと思う人？
３．赤道上の距離より長いと思う人？

というと、半分くらいが「２．」を選びますが、他の半分はまごまごします。ここで、

はださん

北緯90度ってどこですかね？

というと、半分くらいが

あっ（察し）

という顔になり、

はださん

北緯90度は北極点です。つまり北極点から北極点までの距離なので 0km です。

と答えるとほとんどの人が

なんだそんな簡単な問題か。

という顔をします。

この問題は私の受け持つ講座で延べ 300人ぐらいには質問してきましたが、正解が答えられた人はこれまで 4人しかいません。その 4人は一瞬で「0です。」と答えてくれたのですが、ほとんどの方が地図を専門に仕事をしている訳ではないためか答えられません。中学校の社会科で習ったはずなのに。世の中には地球平面協会という地球が平らだと本気で主張している人たちがいますが、彼らを笑うことはできないのです。

多くの人が答えられない原因は、緯度と経度が平面上で直交しているイメージしか持てないからだと考えます。メルカトル図法がその良い例です。私の質問ではホワイトボードにあえて緯度と経度を直交して書くのですが、このように書かれると地球が丸いことを忘れてしまうようです。

Globie

地球は丸いんだよ。僕を見て思い出して。

Google Maps の投影法が変更

そんな認識の問題もあってか、2018年8月3日に Google Maps はデフォルト設定が宇宙から見たように表現される外射図法に変更されました（従来のメルカトル図法で表示することもできます）。こちらの記事では正射図法と結論づけられています。

正射図法で世界を示すと「北緯90度上で経度135度から経度136度までの経線間の距離」は 0m であることがよく分かります。

地球は丸いことを改めて認識したところで、冒頭の緯度・経度の 1度はどれくらいの距離があるのか考えてみます。

前提条件

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

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

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

経度0度上で北緯0度～90度まで1度ずつの間隔で点を置き、その点から経度を＋1度とした2点、緯度は +-0.5度して求めた 2点で「東西・南北における 1度の線分」を作成します。作成した線分の測地線距離を測定します。

緯度・経度で1度、1分の距離と、その場所の測地線距離 1m が何度相当になるかを求めました。測地線距離の説明はこちらの記事をご覧ください。

平面・測地線・航程線・大楕円線、それぞれの距離とは？

はじめに ArcMap には距離や面積が求められる [計測] ツールというのがあります。地図上でクリックして距離や面積がぱっと求められます。 [計測] ツール このツールで [計測タイプの選択] で何...

緯度・経度1度の距離と1mの角距離

北緯 90度における緯度 1度の距離は求めることができますが作成したコードだと 0 で除算することになるので 89 度までの処理としました。

この結果を見ると、以前に書いた 1マイルの距離を求めた記事でも触れたように、厳密には緯度の場所によって緯度 1度の距離も異なることが分かります。緯度1度の長さは極付近で最大となり、赤道へ近づくにつれて微妙に短くなります。経度は赤道が最長で極が 0 となります。緯度1度と経度1度の距離がほぼ等しいのは緯度 5度～ 6度の間ということが分かります。また、以前に書いた地理座標系のデータをそのまま距離に関連する空間解析を行う際に指定するべき Z係数の値ですが、Esri のブログで提示された係数は「"経度" 間 1m 相当の角距離」であることが分かりました。より厳密に係数を求めたい場合はこちらの表が利用できます。

マイルという距離単位

マイレージの記事だと思ってやってきた人はすみません、スターアライアンスとかワンワールドといったお話ではありません。この記事はマイレージの由来である距離単位としてのマイルに関するお話です。 はじめに 以...

地理座標系における Z 係数換算表

最近空間解析を行う機会が多くなり、ジオプロセシング ツールをあれこれ使ってなんだかんだ処理をしています。 地理座標系の標高データ 地形解析でよく利用するツールに傾斜角や陰影起伏があります。 傾斜角ツー...

ソースコード

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

なお、arcpy を含め、ArcGIS で計測される線分の距離は頂点間となります。同一緯線上では赤道以外、線分の両端の最短距離（測地線距離を求めるための線分）と、緯線の長さは微妙に一致しませんが、経度 1度間では約100km に対して 1m 程度と微々たる誤差だったので無視しました。経度の距離間が長くなると誤差は大きくなるので、より正確に求めたいならば、arcpy.Polyline で生成したポリラインを densify で細かく頂点分割してください。

#coding:UTF-8

# ArcGIS での距離計測
# 地球上の 1度間の測地線距離を計算
# 経度0度、北緯0度～90度を中心とした東西1度、南北1度の測地線距離を計算

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

import arcpy

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

# 天文学における地球半径
# 赤道半径：6,371,000m、扁平率の逆数：0
# Authalic Shpere GCS_Shpere GCS_Sphere (4035)
#gcs = arcpy.SpatialReference(4035)

print("Spatial Reference: {0} ({1})".format(gcs.GCSName, gcs.GCSCode))
print("{}\t{}\t{}\t{}\t{}\t{}\t{}".format("緯度", "経度1度長 (m)".ljust(10), "緯度1度長 (m)".ljust(10), "経度1分長 (m)".ljust(10),
"緯度1分長 (m)".ljust(10), "経度1m角 (度)".ljust(10), "緯度1m角 (度)".ljust(10)))
print("{}\t{}\t{}\t{}\t{}\t{}\t{}".format("90", "0".ljust(10), "?".ljust(10), "0".ljust(10),
"?".ljust(10), "None".ljust(10), "?".ljust(10)))

for lat in range(89, -1, -1):
    # 経度
    point_width_from = arcpy.Point(0, lat)
    point_width_to = arcpy.Point(1, lat)
    polyline_width = arcpy.Polyline(arcpy.Array([point_width_from, point_width_to]), gcs)
    length_width = polyline_width.getLength("GEODESIC")

    # 緯度
    point_height_from = arcpy.Point(0, lat - 0.5)
    point_height_to = arcpy.Point(0, lat + 0.5)
    polyline_height = arcpy.Polyline(arcpy.Array([point_height_from, point_height_to]), gcs)
    length_height = polyline_height.getLength("GEODESIC")

    print("{:0=2}\t{:0=11.4f}\t{:0=11.4f}\t{:0=11.4f}\t{:0=11.4f}\t{:0=11.10f}\t{:0=11.10f}".format(
    lat, length_width, length_height, length_width / 60, length_height / 60, 1 / length_width, 1 / length_height))