経緯度 (けいいど、英語 : 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} )とECEF直交座標系 ( 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 ( ϕ ) 1 − e 2 1 − e 2 sin 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}}}=N(\phi ){\frac {1-e^{2}}{1-e^{2}\sin ^{2}\phi }}.\end{aligned}}} 微小量三成分はどれも互いに直交方向となる。 h = 0 {\displaystyle h=0} では回転楕円体となり、また子午線弧 (経線 弧)の曲率半径 は M ( ϕ ) {\displaystyle M(\phi )} 、卯酉線 弧は N ( ϕ ) {\displaystyle N(\phi )} (緯線 弧は N ( ϕ ) cos ϕ {\displaystyle N(\phi )\cos \phi } )[4] となる。
( x , y , z ) {\displaystyle (x,y,z)} から ( λ , ϕ , h ) {\displaystyle (\lambda ,\,\phi ,\,h)} を求める変換計算については上記から導かれる ϕ {\displaystyle \phi } の方程式を解く必要がある[5] 。
回転楕円体面に沿う最短距離の式 詳細は「地理上の距離(英語版)」を参照
回転楕円体面に沿う最短距離(測地線 距離) 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}}}.} ただし、両極が特異点となる。
近距離近似式 上記を率直に一次式( d {\displaystyle d} → Δ {\displaystyle \Delta } )と見なせば、二点間測地線 距離 Δ s {\displaystyle \Delta s} の短距離近似計算式が導出される[6] (平面法と呼ばれることがある[7] )。
Δ 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}}},} ϕ m ≜ ϕ 1 + ϕ 2 2 . {\displaystyle \phi _{\textrm {m}}\triangleq {\frac {\phi _{1}+\phi _{2}}{2}}.} 両極に特異性を持つが[8] 、短距離近似計算としては精度が良く多用される( Δ λ ≃ 0 {\displaystyle \Delta \lambda \simeq 0} もしくは ϕ m ≃ 0 and Δ ϕ ≃ 0 {\displaystyle \phi _{\textrm {m}}\simeq 0\ {\textrm {and}}\ \Delta \phi \simeq 0} の条件下)。
さらに中長距離へ近似精度を改善した計算法も歴史的に多くの研究者によって開発されている。それらは高次の級数計算もしくは反復を含んでいることが多い[9] 。
ガウスの平均緯度法(中間緯度法) 二点間測地線計算の球面近似の一種で[10] 、中距離への近似精度が改善される ( Rapp (1991)[11] §6.4 )。
Δ s = 2 N ( ϕ m ) arcsin ( sin Δ λ 2 cos ϕ m ) 2 + ( cos Δ λ 2 sin Δ 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 {\Delta m}{2N\left(\phi _{\textrm {m}}\right)}}\right)^{2}}}.} ここで子午線弧長 Δ m ≃ M ( ϕ m ) Δ ϕ {\displaystyle \Delta m\simeq M\left(\phi _{\textrm {m}}\right)\Delta \phi } と近似して良い。上記と同様に与えた二点の緯度の中間値 ϕ m {\displaystyle \phi _{\textrm {m}}} を用いる。
式中の Δ λ {\displaystyle \Delta \lambda } 、 Δ ϕ {\displaystyle \Delta \phi } を含む三角関数および arcsin {\displaystyle \arcsin } を一次近似すれば、上記の近距離近似式が得られる。
経度・緯度を並べる順序 並べる順序には、異なる慣行が存在する。正負については、東経 を正の経度 λ {\displaystyle \lambda } 、北緯 を正の緯度 ϕ {\displaystyle \phi } 、南緯 向きを正の余緯度 とする。
右手系 では:(経度 、緯度 、及び高度 )の順とする[12] [13] 。 これに対して左手系 [14] では:(緯度、経度、及び高度)の順とする。局所座標系(地平面)の 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} 方向を右横方向[15] [16] 方位角との対応関係 方位角 は上記と対応した関係が存在する:
方位角を θ {\displaystyle \theta } として、局所座標系(地平面)の単位円は ( x , y , z ) = ( cos θ , sin θ , 0 ) {\displaystyle (x,y,z)=(\cos \theta ,\sin \theta ,0)} となる。
右手系経緯度の採用 下記では右手系経緯度が採用されている。
右手系経緯度を採用しているもののうち、polygon の頂点配列順については時計周り 順(左手系)を採用しているものがある:
左手系経緯度の採用 下記では左手系経緯度(緯度、経度の順)が採用されている。
左手系地図投影法の採用 下記では左手系の地図投影法を採用し、平面座標の x {\displaystyle x} 軸は右横方向が正、 y {\displaystyle y} 軸は下縦方向が正としている[19] 。
脚注 ^ 天体 が球体であれば、球面上の垂直ベクトルは中心を通るので、地理経緯度は地心経緯度に等しい。 ^ 地理経緯度は測地経緯度、測地学的経緯度(geodetic longitude and latitude)とも呼ばれる。 ^ 扁長もしくは扁平楕円体座標系とは異なる。 ^ ムーニエの定理 も参照。 ^ 解くべき ϕ {\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}}}} ^ 球面の場合と同様に、与えた二点の緯度の中間値 ϕ m {\displaystyle \phi _{\textrm {m}}} を用いて経線・緯線弧曲率半径の計算を行う。 ^ 日本では「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 フォーマットでは左手系座標が採用されている。 関連項目