国土地理院の繰り返し計算不要の新算式と順変換反復法の比較

国土地理院が「繰り返し計算のいらない」新算式を公開 — 逆変換の実測比較と、この手法が使える範囲・使えない範囲

国土地理院の測量計算サイトに2026年4月9日、こんなお知らせが掲載されました。「緯度・経度と地心直交座標の相互換算について繰り返し計算の必要ない新しい計算式に変更しています」。地心直交座標（ECEF・XYZ）から緯度経度・楕円体高への逆変換は、長らく「反復計算で解くもの」だった計算です。それが反復なしの閉形式に置き換わりました。

前回の記事「平面直角座標→経緯度座標の変換計算式が公式に「改善」された」では、GeoConverterPro（GCVP）の変換エンジンGeoCoreJPが逆変換を「順変換の反復収束」で解いていることを紹介しました。そこで本記事では2つのことを検証します。ひとつは新算式と順変換反復の精度・速度の実測比較。もうひとつは、この「反復をなくす」手法が座標計算のどこまで通用するのか——ジオイドや地殻変動補正のようなメッシュ計算にも応用できるのか、という適用範囲の評価です。

何が変わったのか

緯度経度・楕円体高（BLH）→地心直交座標（XYZ）の順変換は、昔から1行で書ける直接計算です。

X = (N + h) cosφ cosλ
Y = (N + h) cosφ sinλ
Z = ((1-f)²N + h) sinφ      N: 卯酉線曲率半径

問題は逆変換（XYZ → BLH）です。緯度φを求めるにはNが必要で、Nを求めるにはφが必要——という循環構造があるため、従来は初期値から始めてφを繰り返し更新する反復計算が定番でした。

新算式はこの反復を解消しています。P = √(X²+Y²)、Zp = Z/P などの補助量から、

t0 = Zp × (Q + ã·Zp²) / (Q − ã·f̃²)        （近似解）
t  = Zp / (1 − ã·√( f̃Q / (1 + f̃·t0²) ))   （一段補正）
φ = arctan(t)

と、近似解t0を一度だけ補正して確定させる構成です。ループも収束判定もありません。国土地理院は「計算結果については従来のものと相違がないことを確認しています」と明記しています。

実測比較：新閉形式 vs 順変換の反復

GRS80楕円体で両方式をPython（倍精度）で実装し、緯度0〜89.99°・高さ−1km（海溝）〜2万km（GNSS衛星高度）のグリッドで往復変換誤差（BLH→XYZ→BLHで元に戻るか）を測りました。緯度誤差は地表の距離に換算しています。

手法	地表〜高さ10km	高度400km	高度2万km
国土地理院・新閉形式（反復なし）	数nm（倍精度の限界）	約5µm	約0.2mm
順変換反復 2回	約0.8mm	—	約9cm
順変換反復 3回	数µm	—	約0.15mm
順変換反復 5〜6回（収束）	数nm	数nm	数nm

速度は10万点の一括計算で、新閉形式が約2.2ms、収束までの反復法が約8.4msと、新算式が約4倍高速でした。反復なしで計算量が固定なので、条件分岐がなくベクトル化・並列化にも向いています。

読み取れること

地表付近の測量用途では、どちらも実質的に誤差ゼロ（nm級＝倍精度浮動小数点の限界）。国土地理院の「従来と相違ない」は数値実験でも裏付けられました。
反復法は回数を打ち切ると精度が素直に落ちる（2回でmm級、3回でµm級）。収束閾値をきちんと設けて回し切ることが大事です。
新閉形式にも弱点はあります。P = √(X²+Y²) で割る構造のため、極の真上（P≒0）では別処理が必要です。また高度2万kmではわずかに誤差が残ります（0.2mm、実務では無視できる量）。

評価：この手法はどこまで使えるのか

「反復をなくして閉形式にする」——優秀な手法です。では、緯度経度⇔XYZ以外の計算、たとえばジオイドや地磁気、地殻変動補正のようなメッシュ計算にも同じ手が使えるのでしょうか。結論から言うと、使える計算と使えない計算がはっきり分かれます。

使える：順変換が「解析的な式」で書けている計算

新算式のトリックは、順変換が楕円体の幾何学という数式で書けているからこそ、代数変形で逆向きの閉形式を導けた、という点にあります。同じ性質を持つのは：

計算	状況
地心直交座標⇔緯度経度	今回、閉形式化された
平面直角座標⇔緯度経度（ガウス・クリューゲル）	河瀬式などで既に閉形式化済み。測量計算サイトも採用
Webメルカトル⇔緯度経度など各種投影	多くは元から閉形式、または級数展開で対応可能

つまり「楕円体と投影の幾何学」の世界では、逆変換の閉形式化は王道であり、今回の新算式はその総仕上げに近い位置づけです。

使えない：メッシュ（グリッド）ベースの計算

一方、次の計算群には原理的に適用できません。

計算	順変換の正体
ジオイド高（ジオイド2024）	格子点データのバイリニア補間
セミダイナミック補正（SemiDynaEXE）	補正パラメータ格子の補間
PatchJGD・POS2JGD（地殻変動補正）	同上
地磁気値（偏角・伏角など）	磁気モデル格子の補間

これらの順変換は数式ではなく表引き＋補間です。代数的にひっくり返す「式」がそもそも存在しないため、閉形式の逆変換は作れません。補正後座標→補正前座標のような逆計算が必要な場面では、今も昔も順変換の反復収束が事実上唯一の汎用解法です。

なお、ジオイドと地磁気はそもそも逆計算がほぼ不要な片道計算です（ジオイドは H = h − Ng の引き算、磁気偏角は参照のみ）。逆問題が本当に効いてくるのは、セミダイナミック補正やPatchJGDのような「補正前⇔補正後」の往復が求められる系です。

棲み分けのまとめ

解析式の世界（楕円体・投影）→ 閉形式化が可能で、国土地理院の新算式が優位。
メッシュの世界（ジオイド・地殻変動補正・地磁気）→ 閉形式は原理的に不可能。順変換の反復が本命であり続ける。

GeoCoreJPの「順変換の反復」の意義は変わらない

GeoCoreJPは逆変換を専用の近似式ではなく、順変換を繰り返し呼び出して目標値へ収束させる方式（収束閾値0.1mm・最大反復20回）で解いています。今回の実験は、この方式が回し切れば新算式と同じ機械精度に到達することを示しました。

そして上の適用範囲の評価が示すとおり、GCVPの変換で精度の要となるセミダイナミック補正・PatchJGD系の逆計算はメッシュの世界の住人です。ここでは閉形式という選択肢自体が存在せず、順変換反復が唯一の汎用解法です。GeoCoreJPが解析式・メッシュを問わず同じ設計（順変換の反復）で貫いているのは、公式の算式がどの版でも、順変換さえ正しければ逆変換の精度が自動的についてくるうえ、メッシュ計算にもそのまま通用するからです。設計思想の詳細は「GeoCoreJPの高精度座標変換を支える設計思想」をご覧ください。

地心直交座標（ECEF・XYZ）そのものに馴染みがない方は、姉妹アプリGeoPrism JPの解説記事「地心直交座標（ECEF・XYZ）と緯度経度・標高の関係をビジュアルで学ぶ」が入り口としておすすめです。

まとめ

2026年4月9日、国土地理院の測量計算サイトで地心直交座標⇔緯度経度の相互換算が繰り返し計算不要の新算式に変わった。
数値実験では、新算式は反復なしで倍精度の限界精度に到達し、収束までの反復法より約4倍高速。
ただしこの手法が使えるのは順変換が解析式で書ける計算（楕円体・投影の幾何学）に限られる。
ジオイド・セミダイナミック補正・PatchJGD・地磁気などのメッシュ計算には原理的に適用できず、逆計算が必要な場面では順変換の反復が引き続き唯一の汎用解法。
GeoCoreJPの順変換反復方式は、解析式・メッシュの両世界を同じ設計で扱える点にこそ価値がある。

出典

測量計算サイトトップページ（お知らせ：2026年4月9日）| 国土地理院 https://vldb.gsi.go.jp/sokuchi/surveycalc/main.html
経緯度・高さと地心直交座標の相互換算（計算式）| 国土地理院 https://vldb.gsi.go.jp/sokuchi/surveycalc/surveycalc/algorithm/trans/trans_alg.html