経緯度
経緯度(けいいど、英語: longitude and latitude)とは、経度(longitude)および緯度(latitude)を指し、地球を含む天体表面上で位置(点)を示すための座標表現である。本稿では地理座標系で用いられる経緯度を説明する。
基本的に、その天体の表面点の垂直ベクトルを考え、その向きを球面座標(角度)で表現する[1] 。
地理経緯度
[編集 ]経緯度は基本的にその地表点の垂直ベクトルに基づき、そのベクトルの方向を球面座標で角度表現したものである。
- {経度{\displaystyle \lambda }、緯度{\displaystyle \phi }}⇔{局所垂直ベクトル{\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 }、高度(楕円体高){\displaystyle h})とECEF直交座標系{\displaystyle (x,y,z)}との変換、および微小量の式は下記となる(地球楕円体の長半径 {\displaystyle a}、離心率 {\displaystyle e={\sqrt {f(2-f)}}})。
- {\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}}}
- {\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}}}
微小量三成分はどれも互いに直交方向となる。{\displaystyle h=0}では回転楕円体となり、また子午線弧(経線弧)の曲率半径は{\displaystyle M(\phi )}、卯酉線弧は{\displaystyle N(\phi )}(緯線弧は{\displaystyle N(\phi )\cos \phi })となる[4] [5] 。
{\displaystyle (x,y,z)}から{\displaystyle (\lambda ,,円\phi ,,円h)}を求める変換計算については上記から導かれる {\displaystyle \phi } の方程式を解く必要がある[6] 。
回転楕円体面に沿う最短距離の式
[編集 ]微小量
[編集 ]回転楕円体面に沿う最短距離(測地線距離){\displaystyle s} の微小量は上記から得られる(微分幾何学)。{\displaystyle h=0}({\displaystyle U=0}) の下で、
- {\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}}}.}
ただし、両極が特異点となる。
短距離近似式
[編集 ]二点間測地線距離{\displaystyle \Delta s}は、短距離の場合には、簡素な近似形を導出できる。{\displaystyle \Delta \lambda =\lambda _{1}-\lambda _{2},\ \Delta \phi =\phi _{1}-\phi _{2},} {\displaystyle \phi _{\textrm {m}}={\frac {\phi _{1}+\phi _{2}}{2}}}とおいて、短距離条件は、{\displaystyle |\Delta \phi |\ll 1}かつ{\displaystyle |\cos \phi _{\textrm {m}}\Delta \lambda |\ll 1}と表される。
これに従うと{\displaystyle \Delta s} の近似式が導出される(誤差は {\displaystyle |\Delta s_{\text{approx}}-\Delta s_{\text{true}}|\lessapprox {\frac {\Delta s^{3}}{24a^{2}}}}; 30kmでは1ppmの精度)。
- {\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}}}}.
他の計算式としては、{\displaystyle |\Delta \phi |\ll 1}かつ{\displaystyle |\Delta \lambda |\ll 1}と仮定すると(上記の微小量式を率直に一次式({\displaystyle d}→{\displaystyle \Delta })と見なすことに相当)、より簡素な下記の近似計算式が導出される [7] (平面法などと呼ばれることがある[8] )。
- {\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}}}.}
しかしながら、この {\displaystyle |\Delta \lambda |\ll 1} は高緯度(および極近傍)では必ずしも適切な短距離条件とは言えず、それによる三角関数の近似を行ったことから両極に特異性を生じさせるなど難点を持つが[9] 、高緯度(および極近傍)を除けば短距離近似として妥当であり多用される。
さらに中長距離へ近似精度を改善した計算法も歴史的に多くの研究者によって開発されている。それらは高次の級数計算もしくは反復を含んでいることが多い[10] 。
ガウスの平均緯度法(中間緯度法)
[編集 ]二点間測地線計算の球面近似の一種で[11] 、近似精度が改善される ( Rapp (1991)[12] §6.4 )(誤差は {\displaystyle |\Delta s_{\text{approx}}-\Delta s_{\text{true}}|\lessapprox {\frac {\Delta s^{3}}{400a^{2}}}}; 100kmでは1ppmの精度)。
- {\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] では:(緯度、経度、及び高度)の順とする。局所座標系(地平面)の {\displaystyle x} 方向が北・緯度座標、{\displaystyle y} 方向が東・経度座標となる。
地図投影法の表式における {\displaystyle x,\ y} 平面座標の取り方
[編集 ]地図学における地図投影法の表式で {\displaystyle x,\ y} 平面座標の取り方は右手系で表されることが多い。
- 右手系:{\displaystyle x}方向を右横方向、{\displaystyle y}方向を上縦方向
- 左手系:{\displaystyle x}方向を上縦方向、{\displaystyle y}方向を右横方向[16] [17]
方位角との対応関係
[編集 ]方位角は上記と対応した関係が存在する:
方位角を {\displaystyle \theta }として、局所座標系(地平面)の単位円は {\displaystyle (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の頂点配列順については時計周り順(左手系)を採用しているものがある:
左手系経緯度の採用
[編集 ]下記では左手系経緯度(緯度、経度の順)が採用されている。
- Leaflet
- Google Maps API
- Apple MapKit
- ArangoDB
- GeoRSS
- Open Geospatial Consortium (OGC)の Spatial Reference System (SRS)[19]
左手系地図投影法の採用
[編集 ]下記では左手系の地図投影法を採用し、平面座標の{\displaystyle x}軸は右横方向が正、{\displaystyle y}軸は下縦方向が正としている[20] 。
脚注
[編集 ]- ^ 天体が球体であれば、球面上の垂直ベクトルは中心を通るので、地理経緯度は地心経緯度に等しい。
- ^ 地理経緯度は測地経緯度、測地学的経緯度(geodetic longitude and latitude)とも呼ばれる。
- ^ 扁長もしくは扁平楕円体座標系とは異なる。
- ^ ムーニエの定理も参照。
- ^ 微分関係式は、{\displaystyle {\frac {d[N(\phi )\cos \phi ]}{d\phi }}=-M(\phi )\sin \phi }
- ^ 解くべき {\displaystyle \phi } の方程式は
- {\displaystyle {\begin{aligned}&{\frac {p}{\cos \phi }}-{\frac {z}{\sin \phi }}-e^{2}N(\phi )=0,\\&p={\sqrt {x^{2}+y^{2}}}\end{aligned}}}
- {\displaystyle \kappa -1-{\frac {e^{2}a\kappa }{\sqrt {p^{2}+(1-e^{2})z^{2}\kappa ^{2}}}}=0}
- ^ 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関数を用いる大円距離計算」(円の弦長に基づき弧長を求める)を回転楕円体({\displaystyle e\neq 0})へ拡張した形となっている。
- ^ Rapp, R, H (1991). Geometric Geodesy, Part I (Report). Ohio Start Univ. hdl:1811/24333。
- ^ 和漢の用例でも、この(経度・緯度)の順である「経緯度」である(例えば「日本経緯度原点」、「経緯線」)。
- ^ 右手系の別慣行の変数及び順序は:(余緯度、経度、及び高度)。数学・物理学における球面座標系の標準はこれに当たる。
- ^ a b この左手系の使用は一般的には非推奨とされている。ただし測量、航海術や地理学などの分野はこの左手系の使用は極めて標準的である。
- ^ 左手系の別慣行では、{\displaystyle x}方向を右横方向、{\displaystyle y}方向を下縦方向にとる。
- ^ 平面直角座標系(日本の規格)では左手系である。
- ^ 右手系の別慣行では:(南→東→北→西)
- ^ OGCによるSRS/CRS の定義では大多数の測地系は axis order を左手系経緯度と定義する。
- ^ 他にSVGフォーマットでは左手系座標が採用されている。
関連項目
[編集 ]- 経緯線
- 日本経緯度原点
- IERS基準子午線
- 地理座標系
- メルカトル図法
- 平面直角座標系(日本の規格。左手系)
- ISO 6709 (座標の文字列表現の国際標準。原則として左手系)
- 国家座標
- Vincenty法
この項目は、地理用語に関連した書きかけの項目 です。この項目を加筆・訂正などしてくださる協力者を求めています(プロジェクト:地理/Portal:地理学・Portal:地理)。