測地学 基本 測地学 地理力学(英語版) ジオマティクス(英語版) 地図学 測地学の歴史(英語版) 概念 技術 基準(歴史(英語版) ) NGVD29(英語版) 海面測地系1929年 OSGB36(英語版) イギリス陸上測量1936年 SK-42(英語版) Systema Koordinat 1942 goda ED50(英語版) 欧州座標系1950年 SAD69(英語版) 南米測地系1969年 GRS80 GRS80地球楕円体1980年 NAD83 北米測地系1983年 WGS84 世界測地系1984年 NAVD88(英語版) 北米垂直測地系1988年 ETRS89(英語版) 欧州地球基準座標系システム1989年 GCJ-02 中国の暗号化された座標系2002年 Geo URI(英語版) 地点へのインターネットリンク 2010年
国際地球基準座標系 空間参照系識別子(SRID)(英語版) ユニバーサル横メルカトル図法(UTM)
ユニバーサル横メルカトル図法 (ユニバーサルよこメルカトルずほう)は、国際的に標準化された地図投影法の一種である。略してUTM図法 (Universal Transverse Mercator) とも呼ばれる。主に中縮尺 向けの図法として採用している国が多い。日本では国土地理院 発行の縮尺1:10,000から1:200,000の地形図 に使用されている。
図法の詳細 横メルカトル図法 では通常のメルカトル図法 で歪みの大きくなる高緯度地方を比較的正確に表せるが、逆に基準経度 から東西に離れた地点での歪みが大きくなる。そこで南緯80度から北緯84度までの間を西経180度から東向きに6度ずつ1から60のゾーンに分割して、各ゾーンの範囲をそれぞれの中央経度を基準子午線 とした横メルカトル図法で投影し、60枚の地図を使って両極を除く全地球を描く(正確には地球が扁球 により近いことを考慮したガウス・クリューゲル図法 を用い、両極部分はユニバーサル極心平射図法 を用いる)。
ただし、基準子午線上の縮尺を 1 とした場合、他の部分の局所的縮尺が 1 よりも大きくなり(赤道上で基準経度から3度離れると1.00137)、全体として見ても縮尺が 1 よりも大きくなる(赤道上の基準経度-3度から+3度まで、6度の長さは1.00046倍になる)。これを調整するため、中央子午線上での縮尺係数を0.9996にして、投影範囲全体の平面距離について、その相対誤差の絶対値 を4/10,000以内に収める。
このように分割して地球全体を描けば、それぞれのゾーンの地図は基準子午線から3度以内に収まっており、比較的小さな歪みで済むので、中縮尺でも実用上大きな問題は起きない。
実際には適当な縮尺にして四角い地図に切り分けるが、同じゾーン内でつなぎ合わせれば一枚の平面地図として扱うことができる。ゾーンが異なる場合は平面としてつなぎ合わせることができず、地球の丸みを復元することになる。
なお、「ユニバーサル横メルカトル図法」と呼ぶ場合、6度ごとの基準子午線の取り方と 0.9996 の係数で標準化されたものを指す。
座標換算の簡略式 以下に掲げる座標換算式はドイツの数学者・測地学者 であるヨハン・ハインリヒ・ルイ・クリューゲル(ドイツ語版) により初めて導出され1912年に発表されたもの[1] が元となっており、展開式の初めの数項しか用いていない簡潔さでありながらも中央子午線から約3,000キロメートルの範囲内でミリメートル程度の精度を有している[2] 。日本語による導出の詳細な解説も与えられている[3] 。
地球楕円体 の長半径を a {\displaystyle a\,\!} 、扁平率 を f {\displaystyle f\,\!} とし、中央子午線の経度を λ 0 {\displaystyle \lambda _{0}} とするとき、地理緯度 φ {\displaystyle \,\varphi } 、経度 λ {\displaystyle \,\lambda } の点からUTM座標並びに縮尺係数 k {\displaystyle k\,\!} 及び子午線収差(ドイツ語版) 角 γ {\displaystyle \gamma \,\!} を計算する。便宜上北半球 においては南北座標に N 0 = 0 {\displaystyle N_{0}=0} km を、南半球 においては N 0 = 10000 {\displaystyle N_{0}=10000} km をオフセット値として加える。併せて東西座標には E 0 = 500 {\displaystyle E_{0}=500} km をオフセット値として加え、 k 0 = 0.9996 {\displaystyle k_{0}=0.9996} とする。以下では、距離の単位は km とする。
換算式の表式に先立ち、幾つかの初期値を計算しておく:
n = f 2 − f , A = a 1 + n ( 1 + n 2 4 + n 4 64 + ⋯ ) , {\displaystyle n={\frac {f}{2-f}},\quad A={\frac {a}{1+n}}\left(1+{\frac {n^{2}}{4}}+{\frac {n^{4}}{64}}+\cdots \right),} α 1 = 1 2 n − 2 3 n 2 + 5 16 n 3 , α 2 = 13 48 n 2 − 3 5 n 3 , α 3 = 61 240 n 3 , {\displaystyle \alpha _{1}={\frac {1}{2}}n-{\frac {2}{3}}n^{2}+{\frac {5}{16}}n^{3},\,\,\,\alpha _{2}={\frac {13}{48}}n^{2}-{\frac {3}{5}}n^{3},\,\,\,\alpha _{3}={\frac {61}{240}}n^{3},} β 1 = 1 2 n − 2 3 n 2 + 37 96 n 3 , β 2 = 1 48 n 2 + 1 15 n 3 , β 3 = 17 480 n 3 , {\displaystyle \beta _{1}={\frac {1}{2}}n-{\frac {2}{3}}n^{2}+{\frac {37}{96}}n^{3},\,\,\,\beta _{2}={\frac {1}{48}}n^{2}+{\frac {1}{15}}n^{3},\,\,\,\beta _{3}={\frac {17}{480}}n^{3},} δ 1 = 2 n − 2 3 n 2 − 2 n 3 , δ 2 = 7 3 n 2 − 8 5 n 3 , δ 3 = 56 15 n 3 . {\displaystyle \delta _{1}=2n-{\frac {2}{3}}n^{2}-2n^{3},\,\,\,\delta _{2}={\frac {7}{3}}n^{2}-{\frac {8}{5}}n^{3},\,\,\,\delta _{3}={\frac {56}{15}}n^{3}.}
経緯度 (φ, λ) から UTM 座標 (E, N) への換算 まず、中間変数を以下のとおり定める:
t = sinh ( tanh − 1 sin φ − 2 n 1 + n tanh − 1 ( 2 n 1 + n sin φ ) ) , {\displaystyle t=\sinh \left(\tanh ^{-1}\sin \varphi -{\frac {2{\sqrt {n}}}{1+n}}\tanh ^{-1}\left({\frac {2{\sqrt {n}}}{1+n}}\sin \varphi \right)\right),} ξ ′ = tan − 1 ( t cos ( λ − λ 0 ) ) , η ′ = tanh − 1 ( sin ( λ − λ 0 ) 1 + t 2 ) , {\displaystyle \xi '=\tan ^{-1}\left({\frac {t}{\cos(\lambda -\lambda _{0})}}\right),\,\,\,\eta '=\tanh ^{-1}\left({\frac {\sin(\lambda -\lambda _{0})}{\sqrt {1+t^{2}}}}\right),} σ = 1 + ∑ j = 1 3 2 j α j cos ( 2 j ξ ′ ) cosh ( 2 j η ′ ) , τ = ∑ j = 1 3 2 j α j sin ( 2 j ξ ′ ) sinh ( 2 j η ′ ) . {\displaystyle \sigma =1+\sum _{j=1}^{3}2j\alpha _{j}\cos \left(2j\xi '\right)\cosh \left(2j\eta '\right),\,\,\,\tau =\sum _{j=1}^{3}2j\alpha _{j}\sin \left(2j\xi '\right)\sinh \left(2j\eta '\right).} その上で、換算式は以下のとおりとなる:
E = E 0 + k 0 A ( η ′ + ∑ j = 1 3 α j cos ( 2 j ξ ′ ) sinh ( 2 j η ′ ) ) , {\displaystyle E=E_{0}+k_{0}A\left(\eta '+\sum _{j=1}^{3}\alpha _{j}\cos \left(2j\xi '\right)\sinh \left(2j\eta '\right)\right),} N = N 0 + k 0 A ( ξ ′ + ∑ j = 1 3 α j sin ( 2 j ξ ′ ) cosh ( 2 j η ′ ) ) , {\displaystyle N=N_{0}+k_{0}A\left(\xi '+\sum _{j=1}^{3}\alpha _{j}\sin \left(2j\xi '\right)\cosh \left(2j\eta '\right)\right),} k = k 0 A a { 1 + ( 1 − n 1 + n tan φ ) 2 } σ 2 + τ 2 t 2 + cos 2 ( λ − λ 0 ) , {\displaystyle k={\frac {k_{0}A}{a}}{\sqrt {\left\{1+\left({\frac {1-n}{1+n}}\tan \varphi \right)^{2}\right\}{\frac {\sigma ^{2}+\tau ^{2}}{t^{2}+\cos ^{2}(\lambda -\lambda _{0})}}}},} γ = tan − 1 ( τ 1 + t 2 + σ t tan ( λ − λ 0 ) σ 1 + t 2 − τ t tan ( λ − λ 0 ) ) . {\displaystyle \gamma =\tan ^{-1}\left({\frac {\tau {\sqrt {1+t^{2}}}+\sigma t\tan(\lambda -\lambda _{0})}{\sigma {\sqrt {1+t^{2}}}-\tau t\tan(\lambda -\lambda _{0})}}\right).}
UTM 座標 (E, N, Zone, Hemi) から経緯度 (φ, λ) への換算 北半球に対しては Hemi=+1, 南半球に対しては Hemi=-1 とおく。
まず、中間変数を以下のとおり定める:
ξ = N − N 0 k 0 A , η = E − E 0 k 0 A , {\displaystyle \xi ={\frac {N-N_{0}}{k_{0}A}},\,\,\,\eta ={\frac {E-E_{0}}{k_{0}A}},} ξ ′ = ξ − ∑ j = 1 3 β j sin ( 2 j ξ ) cosh ( 2 j η ) , η ′ = η − ∑ j = 1 3 β j cos ( 2 j ξ ) sinh ( 2 j η ) , {\displaystyle \xi '=\xi -\sum _{j=1}^{3}\beta _{j}\sin \left(2j\xi \right)\cosh \left(2j\eta \right),\,\,\,\eta '=\eta -\sum _{j=1}^{3}\beta _{j}\cos \left(2j\xi \right)\sinh \left(2j\eta \right),} σ ′ = 1 − ∑ j = 1 3 2 j β j cos ( 2 j ξ ) cosh ( 2 j η ) , τ ′ = ∑ j = 1 3 2 j β j sin ( 2 j ξ ) sinh ( 2 j η ) , {\displaystyle \sigma '=1-\sum _{j=1}^{3}2j\beta _{j}\cos \left(2j\xi \right)\cosh \left(2j\eta \right),\,\,\,\tau '=\sum _{j=1}^{3}2j\beta _{j}\sin \left(2j\xi \right)\sinh \left(2j\eta \right),} χ = sin − 1 ( sin ξ ′ cosh η ′ ) . {\displaystyle \chi =\sin ^{-1}\left({\frac {\sin \xi '}{\cosh \eta '}}\right).} その上で、換算式は以下のとおりとなる:
φ = χ + ∑ j = 1 3 δ j sin ( 2 j χ ) , {\displaystyle \varphi =\chi +\sum _{j=1}^{3}\delta _{j}\sin \left(2j\chi \right),} λ 0 = Z o n e × 6 ∘ − 183 ∘ {\displaystyle \lambda _{0}=\mathrm {Z} \mathrm {o} \mathrm {n} \mathrm {e} \times 6^{\circ }-183^{\circ }\,} λ = λ 0 + tan − 1 ( sinh η ′ cos ξ ′ ) , {\displaystyle \lambda =\lambda _{0}+\tan ^{-1}\left({\frac {\sinh \eta '}{\cos \xi '}}\right),} k = k 0 A a { 1 + ( 1 − n 1 + n tan φ ) 2 } cos 2 ξ ′ + sinh 2 η ′ σ ′ 2 + τ ′ 2 , {\displaystyle k={\frac {k_{0}A}{a}}{\sqrt {\left\{1+\left({\frac {1-n}{1+n}}\tan \varphi \right)^{2}\right\}{\frac {\cos ^{2}\xi '+\sinh ^{2}\eta '}{\sigma '^{2}+\tau '^{2}}}}},} γ = H e m i × tan − 1 ( τ ′ + σ ′ tan ξ ′ tanh η ′ σ ′ − τ ′ tan ξ ′ tanh η ′ ) . {\displaystyle \gamma =\mathrm {H} \mathrm {e} \mathrm {m} \mathrm {i} \times \tan ^{-1}\left({\frac {\tau '+\sigma '\tan \xi '\tanh \eta '}{\sigma '-\tau '\tan \xi '\tanh \eta '}}\right).}
UTM直交座標系 UTM の経度6度幅のゾーン内であれば、南緯80度から北緯84度まで、比較的広い範囲を一枚の平面地図として扱える。そこで各ゾーンに直交座標系 を入れて、経緯度 ではなくゾーンとメートル 単位の直交座標を使って地球上の位置指定を行う場合がある。直交座標系の利便性は高く、メッシュを等面積に作れることから統計 分野などでも利用価値が高い。
もっとも使われる UTM-Grid は、赤道をX軸、基準子午線をY軸、両者の交点を原点とし、座標値に負の数 が現れないようX座標値に500キロメートルを加えた座標系である。南半球 の場合はY座標にも1万キロメートルを加える。軍用や国際的用途の他、オーストラリアなど一部の国の国内用直交座標系としても利用されている。
ただし、UTM での直交座標系を基本としながらも、国内用として自国に都合良く原点を設定する例もあり、必ずしも世界中で統一されているわけではない[4] 。また、座標の入れ方が同じであっても、測地系 の選択により僅かな差が出ることがある。さらに、メートル単位で6桁であるX座標値の上の桁として、ゾーン番号や測地系・座標系の種類を区別する独自番号を付加する国もある[5] 。
The U.S. Military Grid Reference System (MGRS) アメリカ軍 やNATO の位置表示システムでは、ゾーンを緯度 8度ずつ南緯80度から72度のCから北緯72度から84度のXまで(ただし数字の1と0に紛らわしいIとOは使用しない)の20に分割して、ゾーン番号と合わせて大雑把な位置指定を行う。
これとは別に、各ゾーン内を100キロメートル四方の区画に分割して、各区画にアルファベット2文字の符号を付与する。文字数が足りないので同じゾーン内に同じ符号が何度も現れるが、上記のアルファベットによる緯度分けと併用することで100キロメートル区画を一意に特定できる。さらに座標や区画分けを加えて位置指定を行う。
日本周辺のUTMゾーン ヨーロッパのUTMゾーン(数字の部分)
採用国 1947年にアメリカ陸軍工兵司令部 が座標システムとして考案したのが始まりとされる[6] 。のちに地形図の図法としても採用する。
ヨーロッパでは各国で測地系がばらばらであったが、NATO の活動に支障が生じるため、西ヨーロッパ全体をカバーする測地系を設けることとなり、それに合わせて一部の国がUTM直交座標系や地形図図法としての UTM を導入した。
日本 では、地形図や地勢図に多面体図法を用いていたが、昭和40年から25,000分の1地形図整備に合わせてUTMに順次切り替えた。日本では、51帯~56帯(中央経線は東経123度、129度、135度、141度、147度、153度)を使用している。他の国と比べると、位置指定やメッシュの手段としてのUTM直交座標系はあまり浸透していない。アメリカ軍との共同運用が多い防衛省 ・自衛隊 を除けば[7] 、1万分の1以上の地形図に日本独自の平面直角座標系 の方眼が入っているが、他は経緯度線が見られる程度である。
脚注・参考文献 ^ Krüger, L. (1912). Konforme Abbildung des Erdellipsoids in der Ebene . Royal Prussian Geodetic Institute, New Series 52. ^ C. F. F. Karney (2011), Transverse Mercator with an accuracy of a few nanometers, J. Geodesy 85 (8), 475–485 ^ 河瀬和重 (2011): Gauss-Krüger投影における経緯度座標及び平面直角座標相互間の座標換算についてのより簡明な計算方法, 国土地理院時報, 121 , 109–124. ^ たとえばイタリア(it:Proiezione di Gauss-Boaga)の東座標系は実質的にUTM33帯北であるが、X座標に加える値は2,520キロメートルである。 ^ たとえばフィンランドのETRS-TM35FINは名前のとおりUTM35帯北であるが、旧来の直交座標系と区別するためX座標に8,000キロメートルを余分に加えることがある。 ^ ただし地球を経度6度帯にゾーン分けをするアイデアは、1891年に作成が決められた世界統一規格の国際図(日本では万国図)に見られる。日本の地形図に付されている地形図番号はUTMゾーン番号と一致している(緯度帯符号はMGRSと異なる)が、これは国際図番号として戦前から使われているものである。またガウス・クリューゲル図法の採用はドイツやソ連が早い。第二次大戦 中に連合国軍として共同作戦を取る上でソ連の図法を採用した事が、UTM考案のきっかけと言われている。ソ連のシステムとUTMの違いは、縮尺係数(ソ連では 1 )とゾーン番号の付け方である(DMA Tecnical manual 8358.1)。 ^ 数値地図ユーザーズガイド(国土地理院監修、日本地図センター編刊、1992年)、p1
関連項目 ウィキメディア・コモンズには、ユニバーサル横メルカトル図法 に関連するカテゴリがあります。
典拠管理データベース
全般 国立図書館