経緯度

経緯度（けいいど、英語: longitude and latitude）とは、経度（longitude）および緯度（latitude）を指し、地球を含む天体表面上で位置（点）を示すための座標表現である。本稿では地理座標系で用いられる経緯度を説明する。

基本的に、その天体の表面点の垂直ベクトルを考え、その向きを球面座標（角度）で表現する^[1]。

経度（

\lambda

）、緯度（

\phi

）、および垂直線（赤）。

ECEF直交座標・地理座標・局所座標の関係（回転楕円体面上）。

(X,Y,Z)

および方位角

\theta

の取り方は右手系。

地理経緯度と天文経緯度[編集]

経緯度は基本的にその地表点の垂直ベクトルに基づき、そのベクトルの方向を球面座標で角度表現したものである。すなわち｛経度 $\lambda$ 、緯度 $\phi$ ｝⇔｛局所垂直ベクトル $(\cos \phi \cos \lambda ,\,\cos \phi \sin \lambda ,\,\sin \phi )$ ｝。

地理座標系で用いられる地理経緯度（geographic longitude and latitude）^[2]は、地球を回転楕円体と見なし、その面の法線ベクトル方向に基づく。

ただし歴史的には、地表の鉛直線に基づく垂直方向（天頂）が天球のどこを指すかによって決めた天文経緯度（astronomical longitude and latitude）が使われてきた。これは地球の重力の鉛直線偏差の影響（加えて地球の極運動の影響）を被っている。従って、距離・面積との関係も簡素にならない。

地理学・測地学の発展とともに、経緯度原点を国内に設け、その地点の天文経緯度を原点として位置づけ、接する準拠楕円体に基づく地理経緯度を用いる方式が行われた（地域的測地系）。

地理経緯度の変換式[編集]

地理座標（経度 $\lambda$ 、緯度 $\phi$ 、高度（楕円体高） $h$ ）とECEF直交座標系 $(x,y,z)$ との変換、および微小量の式は下記となる（地球楕円体の長半径 $a$ 、離心率 $e={\sqrt {f(2-f)}}$ ）。

{\begin{cases}x=\left(N(\phi )+h\right)\cos {\phi }\cos {\lambda },\\y=\left(N(\phi )+h\right)\cos {\phi }\sin {\lambda },\\z=\left(N(\phi )(1-e^{2})+h\right)\sin {\phi },\end{cases}}

{\begin{aligned}{\begin{pmatrix}dx\\dy\\dz\\\end{pmatrix}}&={\begin{pmatrix}-\sin \lambda &-\sin \phi \cos \lambda &\cos \phi \cos \lambda \\\cos \lambda &-\sin \phi \sin \lambda &\cos \phi \sin \lambda \\0&\cos \phi &\sin \phi \\\end{pmatrix}}{\begin{pmatrix}dE\\dN\\dU\\\end{pmatrix}},\\{\begin{pmatrix}dE\\dN\\dU\\\end{pmatrix}}&={\begin{pmatrix}\left(N(\phi )+h\right)\cos \phi &0&0\\0&M(\phi )+h&0\\0&0&1\\\end{pmatrix}}{\begin{pmatrix}d\lambda \\d\phi \\dh\\\end{pmatrix}},\\N(\phi )&\triangleq {\frac {a}{\sqrt {1-e^{2}\sin ^{2}\phi }}},\\M(\phi )&\triangleq {\frac {a(1-e^{2})}{\left(1-e^{2}\sin ^{2}\phi \right)^{3/2}}}.\end{aligned}}

微小量三成分はどれも互いに直交方向となる。 $h=0$ では回転楕円体となり、また子午線弧（経線弧）の曲率半径は $M(\phi )$ 、卯酉線弧は $N(\phi )$ （緯線弧は $N(\phi )\cos \phi$ ）^[3]となる。

$(x,y,z)$ から $(\lambda ,\,\phi ,\,h)$ を求める変換計算については上記から導かれる $\phi$ の方程式を解く必要がある^[4]。

経度・緯度を並べる順序[編集]

並べる順序には、異なる慣行が存在する。正負については、東経を正の経度 $\lambda$ 、北緯を正の緯度 $\phi$ 、南緯向きを正の余緯度とする。

右手系では：（経度、緯度、及び高度）の順とする^[5]^[6]。
これに対して左手系^[7]では：（緯度、経度、及び高度）の順とする。局所座標系（地平面）の $x$ 方向が北・緯度座標、 $y$ 方向が東・経度座標となる。

地図投影法の表式における $x,\ y$ 平面座標の取り方[編集]

地図学における地図投影法の表式で $x,\ y$ 平面座標の取り方は右手系で表されることが多い。

右手系： $x$ 方向を右横方向、 $y$ 方向を上縦方向
左手系： $x$ 方向を上縦方向、 $y$ 方向を右横方向^[8]^[9]

方位角との対応関係[編集]

詳細は「方位角」を参照

方位角は上記と対応した関係が存在する：

右手系（反時計回り）：（東→北→西→南）^[10]
左手系^[7]（時計回り）：（北→東→南→西）

方位角を $\theta$ として、局所座標系（地平面）の単位円は $(x,y,z)=(\cos \theta ,\sin \theta ,0)$ となる。

右手系経緯度の採用[編集]

下記では右手系経緯度が採用されている。

polygonの頂点配列が時計周り順[編集]

右手系経緯度を採用しているもののうち、polygonの頂点配列順については時計周り順（左手系）を採用しているものがある：

左手系経緯度の採用[編集]

下記では左手系経緯度（緯度、経度の順）が採用されている。

Leaflet
Google Maps API
Apple MapKit
ArangoDB
GeoRSS
Open Geospatial Consortium (OGC)の Spatial Reference System (SRS)^[11]

左手系地図投影法の採用[編集]

下記では左手系の地図投影法を採用し、平面座標の $x$ 軸は右横方向が正、 $y$ 軸は下縦方向が正としている^[12]。

Mapbox Vector Tiles

脚注[編集]

^ 天体が球体であれば、球面上の垂直ベクトルは中心を通るので、地理経緯度は地心経緯度に等しい。
^ 地理経緯度は測地経緯度、測地学的経緯度（geodetic longitude and latitude）とも呼ばれる。
^ ムーニエの定理も参照。
^ 解くべき $\phi$ の方程式は
${\begin{aligned}&{\frac {p}{\cos \phi }}-{\frac {z}{\sin \phi }}-e^{2}N(\phi )=0,\\&p={\sqrt {x^{2}+y^{2}}}\end{aligned}}$
で、またこれは変数 $\kappa ={\frac {p}{z}}\tan \phi$ についての方程式に帰着できる：
$\kappa -1-{\frac {e^{2}a\kappa }{\sqrt {p^{2}+(1-e^{2})z^{2}\kappa ^{2}}}}=0$
解き方はGeographic_coordinate_conversion#From_ECEF_to_geodetic_coordinates等を参照のこと。また $h=e^{-2}(\kappa ^{-1}-(1-e^{2})){\sqrt {p^{2}+z^{2}\kappa ^{2}}}$
^ 和漢の用例でも、この（経度・緯度）の順である「経緯度」である（例えば「日本経緯度原点」、「経緯線」）。
^ 右手系の別慣行の変数及び順序は：（余緯度、経度、及び高度）。数学・物理学における球面座標系の標準はこれに当たる。
^ ^a ^b この左手系の使用は一般的には非推奨とされている。ただし測量、航海術や地理学などの分野はこの左手系の使用は極めて標準的である。
^ 左手系の別慣行では、 $x$ 方向を右横方向、 $y$ 方向を下縦方向にとる。
^ 平面直角座標系（日本の規格）では左手系である。
^ 右手系の別慣行では：（南→東→北→西）
^ OGCによるSRS/CRS の定義では大多数の測地系は axis order を左手系経緯度と定義する。
^ 他にSVGフォーマットでは左手系座標が採用されている。