3パラメータと7パラメータ変換 — Helmert変換のしくみ

,

3パラメータと7パラメータ変換 — Helmert変換のしくみ

3パラメータと7パラメータ変換 — Helmert変換のしくみ

異なる測地系(データム)の座標を変換する方法にはいくつかの流儀がありますが、世界的に古くから使われてきた古典的な方法がヘルマート変換(Helmert transformation)です。「3パラメータ変換」「7パラメータ変換」という言葉を目にしたことがあるかもしれません。今回は、この2つの違いと、なぜ日本の測地系変換ではこの方式ではなく別の方法が主役になったのかを整理します。


ヘルマート変換の基本的な考え方

ヘルマート変換は、緯度・経度・楕円体高を一度3次元の直交座標(X, Y, Z)に変換したうえで、

  1. 平行移動(シフト)
  2. 回転
  3. スケール(拡大縮小)

という3種類の操作を組み合わせて別の測地系のXYZに移し、最後にまた緯度経度に戻す、という流れで計算します。剛体を「動かして・回して・伸縮させる」だけの単純なモデルなので、実装も比較的シンプルです。

3パラメータ変換 — 平行移動だけのシンプル版

3パラメータ変換は、このうち平行移動(ΔX, ΔY, ΔZ)だけを使い、回転とスケールは考慮しない最も単純な方式です。2つの測地系の原点位置のズレをそのままXYZ方向のオフセットとして補正するイメージで、計算がシンプルな分、広い範囲をひとつのパラメータで近似すると精度は数m程度に留まることが一般的です。

旧日本測地系から世界測地系への座標のズレ

日本でも、旧日本測地系(Tokyo Datum)から世界測地系への簡易的な変換に、この種の平行移動ベースの手法が使われていた時期があります。国土地理院はTokyo DatumとITRF系(GRS80楕円体)との間の平行移動パラメータを公表しており、大まかな変換であれば少ないパラメータでも実用になります。ただしTokyo Datum自体が明治期からの測量網の蓄積誤差を内包しているため、単純な平行移動だけでは場所によって残差(ズレ)にばらつきが出ることが知られています。

7パラメータ変換 — 回転とスケールも含むフル版

7パラメータ変換(ブルサ・ヴォルフ変換とも呼ばれます)は、平行移動3つに加えて回転3つ・スケール1つを加えた、合計7つのパラメータを使う方式です。

平行移動:ΔX, ΔY, ΔZ(3つ)
回転  :Rx, Ry, Rz(3つ)
スケール:s(1つ)
合計  :7パラメータ

回転とスケールを考慮できる分、3パラメータ変換より広い範囲でより高い精度を実現できます。国際的な測地基準系(ITRFの各エポック間の変換など)でも、この7パラメータ形式が標準的に使われています。

それでも「剛体変換」の限界がある

7パラメータ変換であっても、対象の領域全体をひとつの剛体として扱っている点は変わりません。つまり「全体を一律に動かして・回して・伸縮させる」モデルなので、地震による局所的な地殻変動や、明治期以来の測量誤差の空間的な不均一性のように、場所によって歪み方が違うケースをうまく表現できません。

なぜ日本の高精度変換はヘルマート変換を使わないのか

GeoCoreJPが実装しているTKY2JGD・PatchJGD・SemiDynaEXE・POS2JGDは、いずれもヘルマート変換ではなくグリッド(格子)補正方式です。日本全国を細かいメッシュに区切り、メッシュごとに個別の補正量(緯度方向・経度方向のずれ量)をあらかじめ用意しておき、対象地点は周囲のメッシュ点からバイリニア内挿で補正量を求めます。

この方式なら、メッシュ単位で補正量を変えられるため、「関東は測量誤差の蓄積が大きいが、九州は小さい」「2016年の地震で熊本周辺だけが動いた」といった空間的に不均一な歪みを表現できます。剛体変換であるヘルマート変換では原理的に難しいことが、格子方式なら自然に扱えるわけです。

姉妹サイトGeoPrism JPの「測地系を変換する2つの流儀」では、このパラメータ変換とグリッド変換の違いを可視化で学べます。あわせて読むと理解が深まります。

GeoCoreJPでの位置づけ

GeoCoreJPは日本国内の高精度変換ではグリッド方式を採用していますが、ヘルマート変換(3〜7パラメータ)の考え方そのものは、他国のデータムとの大まかな変換や、精度要件がゆるい場面での基礎知識として今も有効です。剛体変換で足りるか、メッシュ単位の歪み補正が要るか——精度要件に応じて手法を選ぶ、という視点が重要です。


まとめ

  • ヘルマート変換は「平行移動・回転・スケール」を組み合わせる古典的なデータム変換手法
  • 3パラメータ変換は平行移動のみ(シンプルだが広域では精度が数m程度に留まりやすい)、7パラメータ変換は回転・スケールも含む(より高精度だが、それでも剛体変換)
  • どちらも「領域全体を一律に動かす」モデルであり、地震による局所変動や測量誤差の空間的な不均一性は表現できない
  • 日本の高精度変換(TKY2JGD・PatchJGD・SemiDynaEXE・POS2JGD)がグリッド補正方式を採用しているのは、この限界を超えるため

関連記事

出典:
– 国土地理院「測地成果の変換」関連資料 https://www.gsi.go.jp/
– 一般的な測地学教科書におけるヘルマート変換(ブルサ・ヴォルフ変換)の定義


開発者より: アプリ・Kindle本・オープンソースの一覧は GitHub: amru195704 にまとめています。


お願い
本記事の情報は参考目的で掲載しており、正確性・完全性を保証するものではありません。誤記・不正確な情報がございましたら、コメント欄よりご指摘いただければ、確認のうえ修正いたします。


アプリを入手(App Store)GeoConverterPro(座標変換) | GeoPrism JP(測地系の可視化・学習)


コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

Back to top