方位角と子午線収差角 — 座標北と真北の角度を計算する

⚠ 本記事は、姉妹サイトGeoPrism JPの「3つの『北』を学ぶ」（真北・磁北・座標北の違いを可視化で学ぶ概念編）の数式・実装版です。概念の説明はそちらに譲り、ここでは「座標北と真北のズレ（子午線収差角γ）」を実際にどう計算するかを扱います。数式連載（#59楕円体→#52子午線弧長→#51投影式）の第4回にあたります。

3つの北のおさらい（1段落で）

平面直角座標系のような投影座標系では、地図上の「上」（Y軸のプラス方向、座標北）は、その地点の経線が指す「真北」とは一般に一致しません。両者のズレを角度で表したものが子午線収差角γ（gamma）です。真北・磁北・座標北という3種類の「北」が存在する理由や、それぞれの実務上の使い分けは、前掲のGeoPrism JP記事で図解しています。

子午線収差角γの式と直感的な意味

平面直角座標系（ガウス・クリューゲル投影）における子午線収差角の近似式は、投影原点からの経度差Δλと緯度φを使って次のように表せます。

γ ≈ Δλ・sinφ + (三次以降の高次項)

一次近似だけで直感をつかむなら、「原点の経線から東西に離れるほど、また緯度が高くなるほど、座標北は真北から傾く」という関係です。原点の真上（Δλ=0）ではγ=0で座標北と真北は一致し、東西に離れるほど傾きが大きくなります。日本の平面直角座標系19系は原点付近を通る狭い経度帯に系を分割しているため、γの絶対値は実務上小さく収まるよう設計されています。

平面直角座標系19系の区分地図

上図は平面直角座標系1〜19系の区分地図です（系の早見表記事で紹介）。系の境界設計そのものが、各系内でのΔλを抑え、結果としてγを小さく保つ仕掛けになっています。

数値例 — 各系の端でγは何度になるか

系の原点経線から離れるほどγが大きくなる性質上、同じ系内でも「東端・西端」でγの大きさが変わります。目安として、原点から東西に数十km離れた地点では、γはおおむね数分（1分=1/60度）〜十数分程度になることが多く、原点直下ではほぼ0に近づきます。系の幅が狭い（原点付近に収まる）系ほどγの最大値は小さく、逆に南北に長い系や原点から離れた広域をカバーする系では相対的に大きくなる傾向があります。正確な値は系ごとの原点座標と対象地点の経度差から個別に計算する必要があり、本記事では構造の理解を優先し、係数の桁数までは踏み込みません。

座標方位角⇄真方位角の変換

現場で使う頻度が高いのは、次の変換です。

真方位角 = 座標方位角 + γ
座標方位角 = 真方位角 − γ

トータルステーションや測量成果は座標系（平面直角座標系）を基準にした座標方位角で扱うのが一般的ですが、天文観測や磁北からの換算（磁気偏角による補正）を絡める場合は真方位角との相互変換が必要になります。符号の向き（東偏か西偏か、γの正負の定義）は資料によって流儀が異なるため、実装時は原点に対する対象地点の東西どちら側にあるかを明示的に扱うのが安全です。

測量計算・図面での実務注意

図面の「北」表記を確認する：測量図面に描かれる矢印が座標北か真北かは、図面ごとに凡例で確認する必要があります。
系をまたぐ計算では要注意：異なる系（例えば7系と8系）の境界付近を扱う案件では、系ごとにγが異なるため、単純な角度の引き算では済まないケースがあります。
GNSS方位角との整合：RTK-GNSS受信機が出力する方位角は多くの場合、座標系や測地系の設定に依存します。設定と実際に使う成果の基準が一致しているかを確認する習慣が、この種の食い違いを防ぎます。

GeoCoreJP実装での位置づけ

GeoCoreJPでは、前回記事で紹介した投影式（緯度経度→平面直角座標）の周辺関数として、子午線収差角を計算する関数を用意する設計にしています。投影計算自体には収差角は不要ですが、「座標北で描かれた図面の角度を真北基準に直す」といった後続処理でこの値が必要になるため、投影パラメータ（系原点の緯度経度）を共有する形で実装するのが自然な設計です。