経緯度 （けいいど、英語 : longitude and latitude ）とは、経度 （longitude ）および緯度 （latitude ）を指し、地球 を含む天体 表面上で位置（点）を示すための座標表現である。本稿では地理座標系 で用いられる経緯度を説明する。
基本的に、その天体 の表面点の垂直ベクトルを考え、その向きを球面座標 （角度 ）で表現する[ 1]
経度（
λ
{\displaystyle \lambda }
ϕ
{\displaystyle \phi }
ECEF直交座標・地理座標・局所座標の関係（回転楕円体面上）。
(
X
,
Y
,
Z
)
{\displaystyle (X,Y,Z)}
方位角
θ
{\displaystyle \theta }
右手系 。 経緯度は基本的にその地表点の垂直ベクトルに基づき、そのベクトルの方向を球面座標 で角度 表現したものである。
{経度
λ
{\displaystyle \lambda }
ϕ
{\displaystyle \phi }
(
cos
ϕ
cos
λ
,
cos
ϕ
sin
λ
,
sin
ϕ
)
{\displaystyle (\cos \phi \cos \lambda ,\,\cos \phi \sin \lambda ,\,\sin \phi )}
地理座標系で用いられる地理経緯度（geographic longitude and latitude）[ 2] 回転楕円体 と見なし、その面の法線 ベクトル方向に基づく[ 3]
歴史的には、地表の鉛直線 に基づく垂直方向（天頂 ）が天球 のどこを指すかによって決めた天文経緯度（astronomical longitude and latitude）が使われてきた。これは地球の重力の鉛直線偏差の影響（加えて地球の極運動 の影響）を被っている。従って、距離・面積との関係も簡素にならない。
地理学・測地学の発展とともに、経緯度原点を国内に設け、その地点の天文経緯度を原点として位置づけ、接する準拠楕円体 に基づく地理経緯度を用いる方式が行われた（地域的測地系 ）。
さらに近年は全地球的な準拠楕円体 に基づく方式の採用が増えている（全地球的測地系 ）。
地理座標（経度
λ
{\displaystyle \lambda }
ϕ
{\displaystyle \phi }
楕円体高 ）
h
{\displaystyle h}
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
地球楕円体 の長半径
a
{\displaystyle a}
離心率
e
=
f
(
2
−
f
)
{\displaystyle e={\sqrt {f(2-f)}}}
{
x
=
(
N
(
ϕ
)
+
h
)
cos
ϕ
cos
λ
,
y
=
(
N
(
ϕ
)
+
h
)
cos
ϕ
sin
λ
,
z
=
(
N
(
ϕ
)
(
1
−
e
2
)
+
h
)
sin
ϕ
,
{\displaystyle {\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}}}
(
d
x
d
y
d
z
)
=
(
−
sin
λ
−
sin
ϕ
cos
λ
cos
ϕ
cos
λ
cos
λ
−
sin
ϕ
sin
λ
cos
ϕ
sin
λ
0
cos
ϕ
sin
ϕ
)
(
d
E
d
N
d
U
)
,
(
d
E
d
N
d
U
)
=
(
(
N
(
ϕ
)
+
h
)
cos
ϕ
0
0
0
M
(
ϕ
)
+
h
0
0
0
1
)
(
d
λ
d
ϕ
d
h
)
,
N
(
ϕ
)
≜
a
1
−
e
2
sin
2
ϕ
,
M
(
ϕ
)
≜
a
(
1
−
e
2
)
(
1
−
e
2
sin
2
ϕ
)
3
/
2
=
{
N
(
ϕ
)
}
3
(
1
−
e
2
)
a
2
.
{\displaystyle {\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}}}={\frac {\{N(\phi )\}^{3}(1-e^{2})}{a^{2}}}.\end{aligned}}}
微小量三成分はどれも互いに直交方向となる。
h
=
0
{\displaystyle h=0}
子午線弧 （経線 弧）の曲率半径 は
M
(
ϕ
)
{\displaystyle M(\phi )}
卯酉線 弧は
N
(
ϕ
)
{\displaystyle N(\phi )}
緯線 弧は
N
(
ϕ
)
cos
ϕ
{\displaystyle N(\phi )\cos \phi }
[ 4] [ 5]
(
x
,
y
,
z
)
{\displaystyle (x,y,z)}
(
λ
,
ϕ
,
h
)
{\displaystyle (\lambda ,\,\phi ,\,h)}
ϕ
{\displaystyle \phi }
[ 6]
回転楕円体面に沿う最短距離（測地線 距離）
s
{\displaystyle s}
微分幾何学 ）。
h
=
0
{\displaystyle h=0}
U
=
0
{\displaystyle U=0}
d
s
=
d
E
2
+
d
N
2
=
(
N
(
ϕ
)
cos
ϕ
d
λ
)
2
+
(
M
(
ϕ
)
d
ϕ
)
2
.
{\displaystyle ds={\sqrt {dE^{2}+dN^{2}}}={\sqrt {\left(N\left(\phi \right)\cos \phi \,d\lambda \right)^{2}+\left(M\left(\phi \right)d\phi \right)^{2}}}.}
ただし、両極が特異点となる。
二点間測地線 距離
Δ
s
{\displaystyle \Delta s}
Δ
λ
=
λ
1
−
λ
2
,
Δ
ϕ
=
ϕ
1
−
ϕ
2
,
{\displaystyle \Delta \lambda =\lambda _{1}-\lambda _{2},\ \Delta \phi =\phi _{1}-\phi _{2},}
ϕ
m
=
ϕ
1
+
ϕ
2
2
{\displaystyle \phi _{\textrm {m}}={\frac {\phi _{1}+\phi _{2}}{2}}}
|
Δ
ϕ
|
≪
1
{\displaystyle |\Delta \phi |\ll 1}
|
cos
ϕ
m
Δ
λ
|
≪
1
{\displaystyle |\cos \phi _{\textrm {m}}\Delta \lambda |\ll 1}
これに従うと
Δ
s
{\displaystyle \Delta s}
|
Δ
s
approx
−
Δ
s
true
|
⪅
Δ
s
3
24
a
2
{\displaystyle |\Delta s_{\text{approx}}-\Delta s_{\text{true}}|\lessapprox {\frac {\Delta s^{3}}{24a^{2}}}}
Δ
s
=
(
2
N
(
ϕ
m
)
cos
ϕ
m
sin
Δ
λ
2
)
2
+
(
M
(
ϕ
m
)
Δ
ϕ
cos
Δ
λ
2
)
2
{\displaystyle \Delta s={\sqrt {\left(2N\left(\phi _{\textrm {m}}\right)\cos \phi _{\textrm {m}}\sin {\frac {\Delta \lambda }{2}}\right)^{2}+\left(M\left(\phi _{\textrm {m}}\right)\Delta \phi \,\cos {\frac {\Delta \lambda }{2}}\right)^{2}}}}
他の計算式としては、
|
Δ
ϕ
|
≪
1
{\displaystyle |\Delta \phi |\ll 1}
|
Δ
λ
|
≪
1
{\displaystyle |\Delta \lambda |\ll 1}
d
{\displaystyle d}
Δ
{\displaystyle \Delta }
[ 7] [ 8]
Δ
s
=
(
N
(
ϕ
m
)
cos
ϕ
m
Δ
λ
)
2
+
(
M
(
ϕ
m
)
Δ
ϕ
)
2
.
{\displaystyle \Delta s={\sqrt {\left(N\left(\phi _{\textrm {m}}\right)\cos \phi _{\textrm {m}}\Delta \lambda \right)^{2}+\left(M\left(\phi _{\textrm {m}}\right)\Delta \phi \right)^{2}}}.}
しかしながら、この
|
Δ
λ
|
≪
1
{\displaystyle |\Delta \lambda |\ll 1}
[ 9]
さらに中長距離へ近似精度を改善した計算法も歴史的に多くの研究者によって開発されている。それらは高次の級数計算もしくは反復を含んでいることが多い[ 10]
二点間測地線計算の球面近似の一種で[ 11] [ 12]
|
Δ
s
approx
−
Δ
s
true
|
⪅
Δ
s
3
400
a
2
{\displaystyle |\Delta s_{\text{approx}}-\Delta s_{\text{true}}|\lessapprox {\frac {\Delta s^{3}}{400a^{2}}}}
Δ
s
=
2
N
(
ϕ
m
)
arcsin
(
sin
Δ
λ
2
cos
ϕ
m
)
2
+
(
cos
Δ
λ
2
sin
M
(
ϕ
m
)
Δ
ϕ
2
N
(
ϕ
m
)
)
2
.
{\displaystyle \Delta s=2N\left(\phi _{\textrm {m}}\right)\arcsin {\sqrt {\left(\sin {\frac {\Delta \lambda }{2}}\cos \phi _{\textrm {m}}\right)^{2}+\left(\cos {\frac {\Delta \lambda }{2}}\sin {\frac {M\left(\phi _{\textrm {m}}\right)\Delta \phi }{2N\left(\phi _{\textrm {m}}\right)}}\right)^{2}}}.}
並べる順序には、異なる慣行が存在する。正負については、東経 を正の経度
λ
{\displaystyle \lambda }
北緯 を正の緯度
ϕ
{\displaystyle \phi }
南緯 向きを正の余緯度 とする。
右手系 では：（経度 、緯度 、及び高度 ）の順とする[ 13] [ 14] これに対して左手系 [ 15]
x
{\displaystyle x}
y
{\displaystyle y}
x
,
y
{\displaystyle x,\ y}
[ 編集 ] 地図学 における地図投影法 の表式で
x
,
y
{\displaystyle x,\ y}
右手系 で表されることが多い。
右手系 ：
x
{\displaystyle x}
y
{\displaystyle y}
左手系 ：
x
{\displaystyle x}
y
{\displaystyle y}
[ 16] [ 17] 方位角 は上記と対応した関係が存在する：
方位角を
θ
{\displaystyle \theta }
(
x
,
y
,
z
)
=
(
cos
θ
,
sin
θ
,
0
)
{\displaystyle (x,y,z)=(\cos \theta ,\sin \theta ,0)}
下記では右手系経緯度が採用されている。
右手系経緯度を採用しているもののうち、polygon の頂点配列順については時計周り 順（左手系）を採用しているものがある：
下記では左手系経緯度（緯度、経度の順）が採用されている。
下記では左手系の地図投影法を採用し、平面座標の
x
{\displaystyle x}
y
{\displaystyle y}
[ 20]
^ 天体 が球体であれば、球面上の垂直ベクトルは中心を通るので、地理経緯度は地心経緯度 に等しい。^ 地理経緯度は測地経緯度、測地学的経緯度（geodetic longitude and latitude）とも呼ばれる。
^ 扁長もしくは扁平楕円体座標系 とは異なる。
^ ムーニエの定理 も参照。^ 微分関係式は、
d
[
N
(
ϕ
)
cos
ϕ
]
d
ϕ
=
−
M
(
ϕ
)
sin
ϕ
{\displaystyle {\frac {d[N(\phi )\cos \phi ]}{d\phi }}=-M(\phi )\sin \phi }
^ 解くべき
ϕ
{\displaystyle \phi }
p
cos
ϕ
−
z
sin
ϕ
−
e
2
N
(
ϕ
)
=
0
,
p
=
x
2
+
y
2
{\displaystyle {\begin{aligned}&{\frac {p}{\cos \phi }}-{\frac {z}{\sin \phi }}-e^{2}N(\phi )=0,\\&p={\sqrt {x^{2}+y^{2}}}\end{aligned}}}
κ
=
p
z
tan
ϕ
{\displaystyle \kappa ={\frac {p}{z}}\tan \phi }
方程式 に帰着できる：
κ
−
1
−
e
2
a
κ
p
2
+
(
1
−
e
2
)
z
2
κ
2
=
0
{\displaystyle \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
(
κ
−
1
−
(
1
−
e
2
)
)
p
2
+
z
2
κ
2
{\displaystyle h=e^{-2}(\kappa ^{-1}-(1-e^{2})){\sqrt {p^{2}+z^{2}\kappa ^{2}}}}
^ Williams, E. (2013年). “Aviation Formulary. ”. 2024年6月23日 閲覧。 ^ 日本では「Hubeny の（簡易）式」などと呼ばれることもある（ただしその名称は適切ではない）。
^ 180度経線に対しても特異性を持つが、対処は容易である。
^ 例えば「ガウスの平均（中間）緯度法」の式を級数展開したものとして、 Hubeny, K. (1954). Entwicklung der Gauss'schen Mittelbreitenformeln, Österreichische Zeitschrift für Vermessungswesen, Hubeny, K. (1959). Weiterentwicklung der Gauss'schen Mittelbreitenformeln. Zeitschrift für Vermessungswesen.
^ したがって「haversine関数 を用いる大円距離計算 」（円の弦長に基づき弧長を求める）を回転楕円体 （
e
≠
0
{\displaystyle e\neq 0}
^
Rapp, R, H (1991). Geometric Geodesy, Part I (Report). Ohio Start Univ. hdl :1811/24333 。
^ 和漢 の用例でも、この（経度 ・緯度 ）の順である「経緯度」である（例えば「日本経緯度原点 」、「経緯線 」）。^ 右手系 の別慣行の変数及び順序は：（余緯度 、経度 、及び高度 ）。数学・物理学における球面座標系 の標準はこれに当たる。^ a b この左手系 の使用は一般的には非推奨とされている。ただし測量 、航海術 や地理学 などの分野はこの左手系の使用は極めて標準的である。
^ 左手系の別慣行では、
x
{\displaystyle x}
y
{\displaystyle y}
^ 平面直角座標系 （日本の規格）では左手系である。^ 右手系 の別慣行では：（南 →東→北→西）^ OGCによるSRS/CRS の定義では大多数の測地系は axis order を左手系経緯度と定義する。
^ 他にSVG フォーマットでは左手系座標が採用されている。
![{\displaystyle {\frac {d[N(\phi )\cos \phi ]}{d\phi }}=-M(\phi )\sin \phi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e8f66bdbf4b6bf637e02ec238a05b9aad8afbb67)
