経緯度

経緯度（けいいど、テンプレート:Lang-en）とは、経度（テンプレート:Lang）および緯度（テンプレート:Lang）を指し、地球を含む天体表面上で位置（点）を示すための座標表現である。本稿では地理座標系で用いられる経緯度を説明する。

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

経度（$\lambda$）、緯度（$\varphi$）、および垂直線（赤）。

ECEF直交座標・地理座標・局所座標の関係（回転楕円体面上）。$(X,Y,Z)$および方位角$\theta$の取り方は右手系。

地理経緯度

経緯度は基本的にその地表点の垂直ベクトルに基づき、そのベクトルの方向を球面座標で角度表現したものである。

{経度$\lambda$、緯度$\varphi$｝⇔｛局所垂直ベクトル$(\cos \varphi \cos \lambda ,\cos \varphi \sin \lambda ,\sin \varphi )$｝。

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

経緯度の歴史

天文経緯度

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

地理経緯度

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

さらに近年は全地球的な準拠楕円体に基づく方式の採用が増えている（全地球的測地系）。

地理経緯度の変換式

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

$\{\begin{array}{c}x=(N(\varphi )+h)\cos \varphi \cos \lambda ,\\ y=(N(\varphi )+h)\cos \varphi \sin \lambda ,\\ z=(N(\varphi )(1-e^2)+h)\sin \varphi ,\end{array}$

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

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

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

回転楕円体面に沿う最短距離の式

微小量

回転楕円体面に沿う最短距離（測地線距離）$s$ の微小量は上記から得られる（微分幾何学）。$h=0$（$U=0$） の下で、

$ds=\sqrt{d{E}^2+d{N}^2}=\sqrt{{(N\left(\varphi \right)\cos \varphi d\lambda )}^2+{\left(M\left(\varphi \right)d\varphi \right)}^2}.$

ただし、両極が特異点となる。

短距離近似式

二点間測地線距離$\mathrm{\Delta }s$は、短距離の場合には、簡素な近似形を導出できる。$\mathrm{\Delta }\lambda ={\lambda }_1-{\lambda }_2,\mathrm{\Delta }\varphi ={\varphi }_1-{\varphi }_2,$ ${\varphi }_{\text{m}}=\frac{{\varphi }_1+{\varphi }_2}{2}$とおいて、短距離条件は、$|\mathrm{\Delta }\varphi |\ll 1$かつ$|\cos {\varphi }_{\text{m}}\mathrm{\Delta }\lambda |\ll 1$と表される。

これに従うと$\mathrm{\Delta }s$ の近似式が導出される（誤差は $|\mathrm{\Delta }{s}_{\text{approx}}-\mathrm{\Delta }{s}_{\text{true}}|⪅\frac{\mathrm{\Delta }{s}^3}{24a^2}$。したがって30kmでは1ppmの精度）。

$\mathrm{\Delta }s=\sqrt{{(2N\left({\varphi }_{\text{m}}\right)\cos {\varphi }_{\text{m}}\sin \frac{\mathrm{\Delta }\lambda }{2})}^2+{(M\left({\varphi }_{\text{m}}\right)\mathrm{\Delta }\varphi \cos \frac{\mathrm{\Delta }\lambda }{2})}^2}$.

平面法

中緯度・低緯度限定の簡易な短距離近似計算式として、下記の平面法 (flat-surface method) が普及している。

$\mathrm{\Delta }s=\sqrt{{(N\left({\varphi }_{\text{m}}\right)\cos {\varphi }_{\text{m}}\mathrm{\Delta }\lambda )}^2+{\left(M\left({\varphi }_{\text{m}}\right)\mathrm{\Delta }\varphi \right)}^2}.$

導出は、$|\mathrm{\Delta }\varphi |\ll 1$かつ$|\mathrm{\Delta }\lambda |\ll 1$と仮定し、$|\mathrm{\Delta }\lambda |$の三角関数を一次近似して得られる。すなわち上記の微小量式を率直に一次式（$d$→$\mathrm{\Delta }$）と見なすことに相当する。

• しかしながら、この $|\mathrm{\Delta }\lambda |\ll 1$ については高緯度（および極近傍）では必ずしも適切な短距離条件とは言えず、それによる三角関数の近似を行ったことから両極に特異性を生じさせるなど難点を持つが、高緯度（および極近傍）を除けば短距離近似として妥当であり多用される。

ガウスの平均緯度法（中間緯度法）

二点間測地線計算の短距離条件の球面近似の一種で 、上記よりも近似精度が改善される （ Rapp (1991) §6.4 ）（誤差は $|\mathrm{\Delta }{s}_{\text{approx}}-\mathrm{\Delta }{s}_{\text{true}}|⪅\frac{\mathrm{\Delta }{s}^3}{400a^2}$。したがって100kmでは1ppmの精度）。

$\mathrm{\Delta }s=2N\left({\varphi }_{\text{m}}\right)\arcsin \sqrt{{(\sin \frac{\mathrm{\Delta }\lambda }{2}\cos {\varphi }_{\text{m}})}^2+{(\cos \frac{\mathrm{\Delta }\lambda }{2}\sin \frac{M\left({\varphi }_{\text{m}}\right)\mathrm{\Delta }\varphi }{2N\left({\varphi }_{\text{m}}\right)})}^2}.$

さらにより近似精度を改善した計算法も歴史的に多くの研究者によって開発されている（例えば、Bowring's_method_for_short_lines）。また高次の級数計算もしくは反復を含んでいる計算法も多い。

経度・緯度を並べる順序

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

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

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

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

• 右手系：$x$方向を右横方向、$y$方向を上縦方向
• 左手系：$x$方向を上縦方向、$y$方向を右横方向

方位角との対応関係

テンプレート:Main 方位角は上記と対応した関係が存在する：

• 右手系（反時計回り）：（東→北→西→南）
• 左手系（時計回り）：（北→東→南→西）

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

右手系経緯度の採用

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

• OpenLayers
• MapboxGL
• KML
• GeoJSON
• Well-known text
• MySQL
• MongoDB
• Redis
• Oracle Spatial
• Solr
• Elasticsearch
• Geocouch
• OSRM

polygonの頂点配列が時計周り順

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

• Shapefile
• ArcGIS API for JavaScript
• D3.js
• PostGIS
• SpatiaLite
• TopoJSON

左手系経緯度の採用

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

• Leaflet
• Google Maps API
• Apple MapKit
• ArangoDB
• GeoRSS
• Open Geospatial Consortium (OGC)の Spatial Reference System (SRS)

左手系地図投影法の採用

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

• Mapbox Vector Tiles

関連項目

• 経緯線
• 日本経緯度原点
• IERS基準子午線
• 地理座標系
• メルカトル図法
• 平面直角座標系（日本の規格。左手系）
• ISO 6709 (座標の文字列表現の国際標準。原則として左手系）
• 国家座標
• Vincenty法