Gauss-KrügerからWGS84・UTMへの変換ツール

なぜドイツの座標変換ツールが必要なのか

地理空間データを15年扱ってきた中で、ドイツの公共インフラ地図や地籍抜粋、古い測量図を開いてみたら、普通の世界地図ではまったく意味の通らない座標に出くわしたことが何度あったか、もはや数えきれません。数百万台の数値で、4,593,000や5,821,000といった値の周辺に集中しています。これはデータが壊れているのではありません。これはGauss-Krügerであり、ドイツが公式に別の座標系へ移行してから数十年が経った今でも、ドイツの地理空間アーカイブのいたるところに残っています。

これを仕事として扱ってきた者としての率直な意見を言うと、ネット上にある一般的な「座標変換」ツールの多くは、次の二つのどちらかの誤りを犯しています。ひとつは、Gauss-Krügerを任意の横メルカトル図法と互換性があるかのように扱い、誤った回転楕円体を使ってしまうこと。もうひとつは、全国一律の単一の測地系変換の近似を適用し、東西ドイツの歴史的な測量網が再統一以前は互いに統合されていなかったという事実を無視してしまうことです。どちらのミスもコード上はわずかな違いですが、結果への影響は大きく、数十メートルから数百メートルにおよぶ静かな誤差を生みます。これは、実際の地図に結果を重ねてみて、区画の境界線がまったく違う畑に着地して初めて気づくような種類の誤差です。

このツールは、そうした問題を真剣に受け止めて作られています。正しい回転楕円体、正しいゾーンの規約、測地系変換のための東西地域の選択、そして単一の座標だけでなくGeoJSONファイル全体の変換にも完全に対応しています。実際のところ、座標を1つだけ変換したい人はほとんどおらず、座標がぎっしり詰まったファイルを変換する必要があるものだからです。

2つの座標系、1つの分断された歴史

Gauss-KrügerとWGS84/UTMは、同じ時代に、同じ目的のために、同じ種類の機関によって作られたものではありません。その違いを理解することで、実際に直面する癖のほとんどが説明できます。

Gauss-Krüger:19世紀の測量から生まれた座標系

物語はカール・フリードリヒ・ガウス(Carl Friedrich Gauss)から始まります。彼は1821年から1825年にかけての大半を、ハノーファー王国の三角測量の指揮に費やしました。三角測量では角度の測定が中心的な役割を果たすため、Gaussは角度を歪めることなく地球の曲面を平らな紙に投影する方法を必要としていました。彼は横メルカトル図法を応用し、選んだ子午線に沿って南北に細い帯状の範囲を設定することで、その帯の中では歪みが小さくとどまるようにしました。Gaussはこの手法の完全な一般公式を発表することはありませんでした。彼の研究の大部分は、公開された再利用可能な手順としてではなく、私的なノートとハノーファー測量に紐づいた計算例として残されました。

Gaussが素描した数学は、その後1世紀近くほとんど使われないままでした。1871年以前のドイツは数十の独立した邦国が寄せ集まった状態で、各邦国が独自の基準点と独自の局所座標系を用いて、それぞれ独自の土地測量を行っていました。統一後も、ドイツの測地学者たちは国全体をカバーする単一の共有図法をまだ持っていませんでした。この状況を大きく変えたのが、ヨハン・ハインリヒ・ルイ・クリューガー(Johann Heinrich Louis Krüger)の研究です。彼は1912年に、Gaussの図法の厳密で一般的な導出を発表し、緯度経度と投影座標の間の変換に今日も使われている級数展開の公式を完成させました。Krügerの論文は、Gaussの私的な素描を、どの測量局でも適用できる実用的で再現可能な手順に変えました。この図法がGauss-Krügerという2つの名前を持つのは、元のアイデアとそれを使えるものにした人物の両方に敬意を表しているためです。

1920年代以降、ドイツ各州はGauss-Krügerを標準の測量図法として採用し、国土を幅3度のゾーンに分割しました。各ゾーンは東経6度、9度、12度、15度のいずれかの子午線を中心とし、地籍・工学業務に許容できる範囲まで歪みを抑えるものでした。基準となる回転楕円体はBessel 1841で、Krügerの公式が採用された時点ですでに約70年が経過していましたが、既存の三角測量網に十分適合していたため、置き換えることの混乱に見合うだけの利点はありませんでした。

そこに、第二次世界大戦後のドイツ分断が重なります。西ドイツと東ドイツはそれぞれ独自のGauss-Krüger図法を運用し続けましたが、両者の三角測量網は40年以上にわたって互いに調整されることなく、それぞれ独立に整備されました。今日でも見て取れる実際上の結果として、古いGauss-Krügerデータと現代のWGS84との間の最適な変換は、その地点が旧西ドイツ領内にあるか旧東ドイツ領内にあるかによって実際に異なり、この点を多くの変換ツールがひそかに無視しています。これは、ドイツの地図データの中に今なお組み込まれている、冷戦の歴史の小さくも具体的な痕跡です。

1990年の再統一によって、Gauss-Krügerが一夜にして消え去ったわけではありません。ドイツが公式に新しい法的基準系であるETRS89/UTMへ移行したのは、1990年代から2000年代にかけて州ごとに段階的に進められたものであり、地籍局によっては2010年代に入ってようやく移行を終えたところもありました。何百万件もの古い測量記録、公共インフラ地図、地籍区画は再投影されることなく、今なおGauss-Krügerの座標のままで存在しています。古いドイツの地図、公共インフラ、あるいは歴史的な土地記録をデジタル化しようとする人は、遅かれ早かれGauss-Krügerに突き当たります。これこそ、このような変換ツールが今なお必要とされる理由です。

WGS84とUTM:衛星の時代のために作られた座標系

WGS84は、Gauss-Krügerとはまったく異なる出自を持っています。1つの王国を測量して作られたものではありません。アメリカ国防総省が、そもそも衛星測位を実現可能にするために構築したものです。1960年代から1970年代にかけての初期の衛星測地学によって、既存のどの国家的測地系も、地球全体をカバーする単一の測位システムを支えるほどの精度と一貫性を持っていないことが明らかになりました。どの国の測量網も、その国だけに有用な、局所的に最適化されたわずかな誤差を抱えていました。求められていたのは、地域ごとの寄せ集めを内部に持たず、地球上どこでも同じように機能する単一のグローバルな基準座標系でした。

WGS84は1984年に発表され、その後も静かに何度か改良が重ねられ、直近では2013年に更新されました。それぞれの改訂では、衛星追跡データの精度向上に伴い、基準回転楕円体と地球の実際の形状との適合度が高められています。その間も名称や日常的に使う座標そのものは安定していました。これは、Gauss-Krügerが実質的にそうであるような、別々に調整された地域網をつなぎ合わせたものではなく、正真正銘の単一のグローバルシステムです。

ユニバーサル横メルカトル図法(UTM)は、WGS84よりも数十年早く存在していました。アメリカ陸軍が1940年代に標準化しましたが、その根底にあるのは、GaussとKrügerがすでに練り上げていたのと同じ発想です。すなわち、地球を細い横メルカトル図法のゾーンに分割して歪みを小さく抑えるという考え方です。UTMは、ドイツのより狭い幅3度のGauss-Krügerゾーンの代わりに、幅6度の均一なゾーンを1から60まで世界規模で用いているだけです。「WGS84上のUTM」と呼ばれる座標は、1820年代にGaussが編み出したのと同じ横メルカトル図法の数学を、Bessel 1841の代わりに現代のグローバルな回転楕円体の上で、より広いゾーンを使って実行しているにすぎません。

ヨーロッパは1990年代に、ほとんどのGIS用途においてWGS84と機能的に同一といえる、ヨーロッパに固定された基準座標系ETRS89を標準として採用しました。今日のドイツの公式座標はETRS89/UTMであり、西部はゾーン32N、東部はゾーン33Nです。

数学的背景:この図法が実際にどう機能するか

ここでは、ギリシャ文字がずらりと並ぶだけでその物理的な意味がまったく説明されない、よくあるやり方ではなく、GISを学ぶ学生に説明するときのやり方でこれを解説していきます。

ステップ1:回転楕円体を知る

どの図法も、地球の形状を表す数学的モデル、すなわち両極でわずかに扁平になった回転楕円体から出発します。この回転楕円体は、長半径(赤道半径)であるaと、扁平率であるfという2つの数値だけで完全に定義されます。Gauss-Krügerのデータが基づくBessel 1841回転楕円体の場合:

• a = 6,377,397.155メートル
• 1/f = 299.1528128(したがってf ≈ 0.0033427731)

aとfから、以下の式に登場する他のすべての量が導かれます。とりわけ重要なのは、第一離心率の二乗e² = 2f − f²と、第二離心率の二乗e'² = e²/(1−e²)です。この2つの数値は、この先ずっと繰り返し登場します。

ステップ2:子午線弧長(この図法全体を支える背骨)

最も重要な構成要素は、子午線弧長M(φ)です。これは文字通り、赤道から真北に向かって地表を歩いていったとき、緯度φに到達するまでにどれだけの距離を歩く必要があるかを表します。地球は球ではなく回転楕円体であるため、この距離はφに対して単純に線形にはスケールせず、簡単な閉じた形の式では表せません。そのため、e²についての収束級数として計算されます。

M(φ) = a · [ (1 − e²/4 − 3e⁴/64 − 5e⁶/256) · φ
            − (3e²/8 + 3e⁴/32 + 45e⁶/1024) · sin(2φ)
            + (15e⁴/256 + 45e⁶/1024) · sin(4φ)
            − (35e⁶/3072) · sin(6φ) ]

最初の項以降のすべての項は、回転楕円体の扁平に対する補正であり、後続の各項はおおよそe²(Bessel 1841では約0.0067)の倍率で前の項より小さくなっていきます。そのため、4つの項だけでも、ドイツの緯度の全範囲にわたってミリメートル未満の精度が得られます。これが北距(Northing)の土台であり、これに加えられるその他すべての項は、中央子午線からどれだけ東西に離れているかに対する補正です。

ステップ3:東西方向の補正項

さらに3つの量が、いずれも緯度φにおいて計算され、中央子午線からの経度距離に対する補正を決定します。

• N = a / √(1 − e²sin²φ)(卯酉線曲率半径。この緯度において経度1度が実際にどれだけの距離に相当するかを表す)
• T = tan²φ
• C = e'²cos²φ

中央子午線からの経度差をラジアンで表したlを用いると、完全な正方向の式(Snyder, 1987)は次のとおりです。

Easting  = N·cosφ·l + (N/6)·cos³φ·(1−T+C)·l³
                     + (N/120)·cos⁵φ·(5−18T+T²+14C−58TC)·l⁵

Northing = M(φ) + (N/2)·sinφ·cosφ·l²
                 + (N/24)·sinφ·cos³φ·(5−T+9C+4C²)·l⁴
                 + (N/720)·sinφ·cos⁵φ·(61−58T+T²+270C−330TC)·l⁶

この一見すると怖じ気づくような式に見えますが、構造自体は見た目より単純です。これはlのべき乗によるテイラー級数展開であり、単純な平面近似(N·cosφ·l)を、中央子午線から離れるにつれて回転楕円体の曲率がどれだけものを歪めるかに応じて補正するものです。Gauss-Krügerの幅1.5度という狭い半ゾーンの範囲内では、l⁵やl⁶の項はすでに1ミリメートル未満しか寄与しません。この規模では、級数は極めて速く収束します。

ステップ4:投影の前後で行う測地系変換

ここまでの内容は、すべて1つの回転楕円体の内部だけで完結しています。Bessel 1841(Gauss-Krüger)とGRS80(WGS84)の間を行き来するには、別のステップが必要です。それが、緯度経度に直接ではなく3次元の地心XYZ空間で適用される、7パラメータのヘルマート変換です。この変換では、地理座標を地心XYZ座標に変換し、3つの並進、3つの微小な回転、1つの縮尺補正を適用してから、再び元に戻します。

X' = dx + (1+s)·( X − rz·Y + ry·Z )
Y' = dy + (1+s)·( rz·X + Y − rx·Z )
Z' = dz + (1+s)·( −ry·X + rx·Y + Z )

実務者としての意見を言うと、ほとんどの人がここでつまずくのは、式が難しいからではなく、このステップをまるごと省略して、Gauss-Krügerの東距・北距をあたかもWGS84の回転楕円体上の値であるかのようにそのまま使おうとするからです。それは決して正しくありません。図法の数学と測地系変換は互いに独立した2つの操作であり、どちらも軽視できません。

どのゾーンと地域を選べばよいか

• 旧西ドイツ: バーデン=ヴュルテンベルク州、バイエルン州、ブレーメン州、ハンブルク州、ヘッセン州、ニーダーザクセン州、ノルトライン=ヴェストファーレン州、ラインラント=プファルツ州、ザールラント州、シュレースヴィヒ=ホルシュタイン州
• 旧東ドイツ: ベルリン、ブランデンブルク州、メクレンブルク=フォアポンメルン州、ザクセン州、ザクセン=アンハルト州、テューリンゲン州

どちらに該当するか分からない場合は、まずどちらかの地域を選んで座標を変換し、「この場所に正しいGK/UTMゾーンを使用する」をクリックしてください。このツールは簡単な地理的な目安(おおよそ東経10.5度より東、北緯50.8度より北であれば旧東ドイツと判定)に基づいて、正しいゾーンと地域を自動的に選択します。

手計算で実践:ブランデンブルク門を変換する

仮の数値ではなく、実際の数値を使います。これは、よく知られたブランデンブルク門のWGS84座標を出発点として、このツールが実際に実行している計算の連鎖を、1ステップずつ示したものです。

出発点(WGS84):北緯52.516266度、東経13.377775度。

ステップ1:地理座標 → 地心XYZ座標(GRS80/WGS84上)

地球を中心とした3次元直交座標に変換します。

X = 3,783,993.8378 m
Y =   899,923.5184 m
Z = 5,037,966.2582 m

ステップ2:逆ヘルマート変換を適用する(西ドイツのパラメータ)

符号を反転させた西ドイツ用パラメータセット(dx=−598.1、dy=−73.7、dz=−418.2、rx=−0.202″、ry=−0.045″、rz=+2.455″、s=−6.7 ppm)を用いて、WGS84からDHDNへと戻します。

X' = 3,783,358.5749 m
Y' =   899,893.7601 m
Z' = 5,037,514.2481 m

ステップ3:地心座標 → 地理座標(Bessel 1841上)

X'、Y'、Z'を、今度はBessel 1841回転楕円体を用いて緯度経度に戻します。

DHDN latitude  = 52.51767297°
DHDN longitude = 13.37951396°

これはWGS84の座標とは一致しない点に注意してください。DHDNの緯度は約0.0014度北寄りに、経度は約0.0017度東寄りになります。地上距離にしておよそ150から190メートルにあたるこのわずかなずれは、すべて異なる回転楕円体と異なる歴史的基準座標系が積み重なった結果です。だからこそ、WGS84の座標にDHDNというラベルを貼り替えるだけでは済まないのです。

ステップ4:Gauss-Krügerゾーン4(中央子午線東経12度)へ投影する

φ = 52.51767297度(0.9166063ラジアン)、中央子午線からの経度差l = 1.379514度(0.0240771ラジアン)とすると、中間量は次のようになります。

N (radius of curvature) = 6,390,841.3354 m
T = tan²φ               = 1.700567
C = e'²cos²φ             = 0.002488
M (meridian arc length)  = 5,820,346.7101 m

これらの値を、前述の数学的背景の項で示した東距・北距の級数に代入すると:

Easting (before false easting)  =    93,631.7509 m
Northing                        = 5,821,241.2561 m

ゾーン4の加数(false easting)は4,500,000(先頭の「4」がゾーンを示す)であり、これを上で求めた素の東距に加えると:

Rechtswert = 4,500,000 + 93,631.75 = 4,593,631.75
Hochwert   =                          5,821,241.26

比較:同じ地点をUTMゾーン33Nで表すと

同じ元のWGS84座標をもう一度横メルカトル図法の式に通します。今回はGRS80回転楕円体上で直接計算するためデータム変換は不要で、中央子午線は東経15度、縮尺係数は0.9996です。結果は次のとおりです。

UTM 33N Easting  =   389,922.84 m
UTM 33N Northing = 5,819,698.02 m

どちらの結果も、地球上のまったく同じ物理的地点を表しています。2つの北距の値の間にある約370万メートルもの差は、実際の距離とは何の関係もありません。これは純粋に、Gauss-Krügerの加数・ゾーンの規約とUTMのそれとの違い、それにBessel 1841とGRS80のどちらに投影するかによるわずかなずれが合わさったものです。これはまさに、この変換ツールが取り除こうとしている種類の混乱です。

手順:変換ツールの使い方

ステップ1

「変換元」で元の形式を、「変換先」で目的の形式を選びます。Gauss-Krügerを選んだ場合は、データが測量された中央子午線に対応する正しいゾーン(2から5)と、その地点が旧西ドイツと旧東ドイツのどちらに該当するかも合わせて選択してください。

ステップ2

単一の地点を変換する場合は、「単一座標」タブに切り替えて、選んだ変換元の形式に応じたRechtswert/Hochwert、緯度/経度、または東距/北距の2つの値を入力します。

ステップ3

「変換」をクリックすると、計算の中間段階として使われたWGS84座標とともに、結果がすぐ下に表示され、地図上にもマーカーが表示されるので、位置が正しいかどうかをひと目で確認できます。ゾーンの選択を誤って入力してしまった場合も、たいていこの段階ですぐに気づきます。

ステップ4

ファイル全体を一括で変換するには、「GeoJSONファイル」タブに切り替えます。ファイルをアップロードするか内容を直接貼り付け、同じように変換元・変換先の形式を設定してから、「変換してプレビュー」をクリックします。

ステップ5

変換された図形が地図上に描画されるので、正しい位置にあるかどうかを確認できます。パネルには、変換された座標ペアの数も正確に表示されます。「converted.geojsonをダウンロード」をクリックすると結果を保存できます。すべてのフィーチャーの元のプロパティは、そのまま一切変更されずに保持されます。

このツールを実際に必要としているのは誰か

順不同で、実際にこの種のツールがどのように使われているのを見てきたかを挙げます。

1. 古いインフラ網の地図をデジタル化する公益事業会社。 2000年代より前のガス・水道・電気の配管図の多くは、今でもGauss-Krügerのままです。資産管理用のGISを刷新するということは、何十年分もの配管やケーブルの形状データを、精度を落とさずに変換することを意味します。
2. 古い地籍記録と現代のGNSS測量を照合する土地家屋調査士。 1975年に測量された境界標と、先週RTK-GPS受信機で測定された境界標を比較するには、両者が一致しているかどうかを判断する前に、まず同じ座標系に揃える必要があります。
3. 歴史的なドイツの地図をデジタル化するアーキビストや歴史家。 戦前や冷戦時代の地形図・地籍図にはGauss-Krügerのグリッドが記載されており、スキャンした地図の地物情報を現代のGIS形式に変換するには、まさにこの変換が必要です。
4. 古い観測点網を扱う環境コンサルタント。 数十年前に整備された地下水・土壌・大気質の観測網では、元の資料に記載された観測点の位置が今でもGauss-Krügerで表されていることが少なくありません。
5. ドイツ政府のオープンデータポータルと連携するソフトウェア開発者。 古いものや地域限定のドイツのオープンデータのGeoJSONエクスポートの中には、仕様上はWGS84専用のはずの形式の中に、今でもGauss-Krügerの座標がそのまま入っているものがあります。変換しないままにしておくと、Webマップは何のエラーも出さずに、たいていギニア湾付近といった見当違いの場所に地物を描画してしまいます。
6. 古い変電所や系統連系のデータを評価する風力・太陽光発電の計画担当者。 再生可能エネルギーの整備が本格化する以前からある送電インフラの記録は、今でも古い図法のままであることがよくあります。
7. 測地学を学ぶ大学の地理学・測量学専攻の学生。 Gauss-Krügerは、その歴史的な重要性に加え、上記の計算例のように変換を実際に手で計算してみることで概念が本当の意味で身につくという理由から、ドイツのVermessungskunde(測量学)の授業で今も教えられています。
8. クライアントの既存の空間データベースのデータム整合性を監査するGISコンサルタント。 一部のレイヤーは気づかれないままDHDNのままで、別のレイヤーはETRS89になっており、どちらがどちらなのか誰もラベル付けしていない、といった混在データムのデータベースは非常によくある問題であり、発見が遅れるほど代償が大きくなります。このツールは、そうした問題を素早く診断する助けになります。

技術的なFAQ

回答はあえて一般論ではなく、具体的な内容にとどめています。これらの質問をするだけの知識をすでにお持ちの方には、それにふさわしい実のある答えが必要だからです。