› ランベルト正角円錐図法 シミュレーター 戻る
地図学・GIS

ランベルト正角円錐図法 シミュレーター

2つの標準緯線で地球を「割る」ように接する円錐に地表を投影する LCC（Lambert Conformal Conic）の挙動を、円錐定数 n、投影座標 (x, y)、スケール係数、歪みのリアルタイム計算で体験できるツールです。米国 SPCS・国際航空図 ONC・気象チャートで実際に使われている中緯度の標準投影の感覚をつかめます。

パラメータ設定
標準緯線 1 φ₁
°
円錐が地球と交差する南側（または低緯度側）の標準緯線
標準緯線 2 φ₂
°
高緯度側の標準緯線。φ₁ と一致すると接円錐になる
中心子午線 λ₀
°
投影の中央経度。SPCS では各州の中央付近を選ぶ
基準緯度 φ₀
°
投影面の原点となる緯度（false northing 起点）
投影点緯度 φ
°
平面に投影したい地表点の緯度
投影点経度 λ
°
平面に投影したい地表点の経度
地球半径 R
km
球モデルの地球半径。WGS84 平均は 6371 km
計算結果
—
円錐定数 n
—
投影座標 X (km)
—
投影座標 Y (km)
—
スケール係数 k
—
歪み (%)
—
円錐頂点距離 (km)
—
円錐投影の幾何イメージ

地球儀に2つの標準緯線（φ₁・φ₂）で交差する円錐。投影点（赤）はその円錐面に写像され、円錐を平面に展開すると下の (x, y) 平面座標になります。

スケール係数 vs 緯度（歪みカーブ）
投影平面の経緯線グリッド
理論・主要公式

$$n = \frac{\ln(\cos\phi_1/\cos\phi_2)}{\ln(\tan(\pi/4+\phi_2/2)/\tan(\pi/4+\phi_1/2))},\quad x = \rho\sin(n(\lambda-\lambda_0))$$

n=円錐定数（0<n<1）、ρ=円錐軸からの半径、(x, y)=投影平面座標（km）、(φ, λ)=測地座標、λ₀=中心子午線、φ₁・φ₂=標準緯線。

$$\rho = R\,F/\tan^{n}(\pi/4+\phi/2),\qquad F = \cos\phi_1\,\tan^{n}(\pi/4+\phi_1/2)/n$$

ρ は緯度 φ から円錐頂点までの平面上距離。F は標準緯線 φ₁ で k=1 となるように決める正規化係数。R は地球半径。

$$k = \frac{\cos\phi_1\,\tan^{n}(\pi/4+\phi_1/2)}{\cos\phi\,\tan^{n}(\pi/4+\phi/2)}$$

k はスケール係数（緯線方向）。両標準緯線間で k<1、外側で k>1、両標準緯線上で k=1。正角投影なので同じ点では子午線方向のスケールも k に等しい。

ランベルト正角円錐図法 — 航空・気象地図の標準投影

🙋
「ランベルト円錐図法」って、世界地図を作る投影法のひとつですよね？メルカトルや UTM とどう違うんですか？
🎓
ざっくり言うと「対象範囲に円錐をかぶせる」発想なんだ。メルカトルは赤道に接する円筒、UTM はゾーンごとの横軸円筒なんだけど、ランベルトは中緯度の2本の緯線（φ₁、φ₂）で地球と交差する円錐を被せる。米国の SPCS や国際航空図 ONC、気象予報チャートで標準的に使われているよ。スライダーで φ₁=33°, φ₂=45° のままで n を見ると 0.6492 になっているはず——これが円錐の開き具合だ。
🙋
n が 0.6492 ってことは、円錐は丸い円錐の「6.5割」分だけ開くってことですか？
🎓
そう、平面に展開すると 360° の地球の経度幅が n×360° = 約 234° の扇形になる、というイメージだね。極限を考えると分かりやすい。φ₁ を0に近づけると n も0に近づいて、円錐がほぼ円筒（メルカトル）になる。逆に φ₁=φ₂=89° みたいに極側に寄せると n はほぼ1で、円錐がほぼ平面（極ステレオ）になる。ランベルトはその中間を、中緯度に最適化したものなんだ。
🙋
「正角」というと角度を保つって意味ですよね？でも歪みが出てるってどういうことですか？
🎓
いいところに気づいたね。「正角（conformal）」というのは「ある点の周りでの局所的な角度」を保つ性質で、つまり微小な円が平面でも円のまま（楕円にならない）ということ。ただし大きさは緯度ごとに変わる。スケール係数 k がそれを表していて、両標準緯線（φ₁・φ₂）の間では k<1（縮小）、外側では k>1（拡大）になる。下の歪みカーブを見て、φ=40° で k≈0.997、つまり 0.3% 縮小していることが分かるはずだ。
🙋
なるほど、両標準緯線の間に挟むと歪みが小さくなるんですね。SPCS だとどう設計してるんですか？
🎓
州ごとに φ₁・φ₂ をその州の南北端から内側に約 1/6 ずつ入った位置に置くんだ。たとえばコロラド州（北緯約 37°〜41°）では φ₁=37°50′、φ₂=39°45′ という具合。こうすると k は州内で 0.9999〜1.0001 の範囲、つまり距離誤差 1万分の1 程度に収まる。1:24,000 の地形図で 10cm 単位のズレも出ないレベル。航空図 ONC（1:1,000,000）も同様で、各シートで標準緯線をその範囲の境界寄りに置いて歪みを ±0.05% 以下にしてる。
🙋
φ₁ と φ₂ を同じ値にすると「接円錐」になると書いてありますが、それと「割円錐」はどう違うんですか？
🎓
φ₁=φ₂ の接円錐は、その緯線でちょうど k=1 になる単一標準緯線。そこから離れるほど一方向に歪みが増える。一方で φ₁≠φ₂ の割円錐（secant、＝この計算式の標準）は、2本の標準緯線で k=1、その間は k<1、外側は k>1 と「両端から内側に歪みを分散」できるんだ。フランスの NTF 旧座標系は接円錐（φ₀=48°）、現代の RGF93 Lambert-93 は割円錐（φ₁=44°、φ₂=49°）。実務ではほとんど割円錐が使われるよ。

よくある質問

東西に広がる中緯度地域（おおむね北緯または南緯 30°〜60°）の地図に向いています。米国の SPCS（State Plane Coordinate System）は東西に長い州の多くで LCC を採用し、国際航空図 ICAO ONC（1:1,000,000）、気象予報チャート、フランスの NTF/RGF93、カナダの一部地理院図でも標準投影として使われます。2つの標準緯線（φ₁, φ₂）で円錐を地球と「割る」ように交差させると、両標準緯線の間ではスケールが1未満（縮小）、外側では1超（拡大）になり、対象範囲を両線で挟むだけで歪みを ±1% 以下に抑えられます。
n は投影で使う円錐の「開き具合」を決めるパラメータです。地球の中心を頂点とする無限大の半径円錐と地球面の交線が標準緯線になり、その円錐を平面に展開したときの開き角 2π·n が決まります。φ₁=φ₂（接円錐）では n = sin(φ₁) で、φ₁ が低緯度なら円錐は浅く（n が0に近づくとほぼ円筒＝メルカトル類似）、高緯度なら円錐は鋭く（n が1に近づくと立体写像＝極ステレオに近づく）なります。LCC は n=0 のメルカトル、n=1 のステレオの中間に位置する一般化された正角投影です。
スケール係数 k は「地球上での1単位距離が平面で何単位として描かれるか」の比率です。k=1 が正確、k=1.02 は2%拡大、k=0.98 は2%縮小を意味します。1:50,000 の航空図で2%の歪みがあると、10km の航路で 200m のずれになり、計器進入の誤差規定を超える恐れがあります。SPCS や ONC では 1標準緯線あたり ±0.05% 以下になるよう設計範囲を限定し、本ツールでは |歪み| ≤ 2% を黄色警告、≤ 10% を赤判定としています。
両方とも正角投影（局所的に角度を保存）ですが、対象範囲が違います。メルカトルは円筒を赤道に接する形なので赤道付近で歪み最小、高緯度では面積が爆発的に拡大します（グリーンランドがアフリカ並みに大きく見える有名な歪み）。ランベルト円錐は円錐を中緯度の標準緯線に交差させるため、中緯度で歪みが最小化されます。航空航海用途では、メルカトルは赤道付近〜熱帯（東南アジア航路など）、ランベルトは中緯度（北米・欧州航空図）、極ステレオは極域（北極ルート）と使い分けます。

実世界での応用

国際航空図 ICAO ONC（1:1,000,000）：世界を 270 枚のシートに分けて発行されている国際標準の航空図で、ほぼ全てのシートが LCC を採用しています。各シートは緯度幅 4°、経度幅 6°（赤道付近）と狭く、標準緯線をシート上下から 1/6 ずつ入った位置に置くことで、シート内の歪みを ±0.05% 以下に抑えています。これにより、航路図上で定規を当てて測った距離がほぼ実距離として読めるため、計器進入や航法計算の基準として今も世界中で使われています。

米国 State Plane Coordinate System（SPCS）：東西に細長い州（カリフォルニア・テキサス・コロラドなど）の州面座標系で標準的に採用されている投影法です。たとえばコロラド州中央ゾーンでは φ₁=37°50′、φ₂=39°45′、λ₀=-105°30′ という設定で、州内のあらゆる地点で k が 0.9999〜1.0001 の範囲に収まります。地籍・測量・道路設計・GIS の基盤座標として、不動産登記から自動運転 HD マップまで幅広く使われます。NAD83 と組み合わせるのが現代の標準です。

気象予報チャート：WMO（世界気象機関）の中緯度の高層天気図、500hPa 解析図、ジェット気流マップなどは LCC で描画されます。気象現象は等圧線や前線の「角度」が重要なので正角投影が好まれ、かつ中緯度帯（偏西風帯）が解析対象なのでランベルトが理想的。NCEP/GFS や ECMWF のチャートも内部的に LCC グリッドで出力され、各国気象庁が地上局へ配信しています。

フランス・カナダの国家座標系：フランスは旧来の NTF（4ゾーンの接円錐 Lambert）から、現代の RGF93/Lambert-93（割円錐、φ₁=44°、φ₂=49°、λ₀=3°、φ₀=46.5°）に移行しました。1つの座標系でフランス本土全域をカバーでき、IGN 地形図・地籍・カーナビ・洪水ハザードマップなどあらゆる国家インフラの基盤となっています。カナダの航空図・地形図の一部も LCC ベース。GIS で .prj ファイルに「PROJCS["NTF (Paris) / Lambert ..."]」と書かれていれば、まさにこの投影です。

よくある誤解と注意点

まず、「正角だから距離も正確」という思い込み。LCC は局所的な角度を保つだけで、面積も距離も保存されません。スケール係数 k は緯度に依存するため、たとえば φ=30° で k=1.005、φ=50° で k=1.008 という地図で「両者間の距離」を平面の定規で測ると、実距離との誤差は約 0.7% になります。1,000km の距離で 7km のずれです。航路計画や測量では、必ず緯度ごとの k を考慮して補正してください。本ツールでも「歪み (%)」が表示されますが、これはあくまで「その投影点での局所スケール誤差」であり、長距離の累積誤差ではありません。

次に、「球モデルと楕円体モデルを混同すること」。このツールは説明の単純化のため球形地球（半径 R）で計算していますが、実務で使う SPCS や RGF93 の LCC は GRS80/WGS84 楕円体ベース（半長軸 a=6378137m、扁平率 1/298.257）です。両者の違いは緯度依存で、北緯 45° で計算結果に約 0.3% の差が出ます。GIS ソフト（QGIS、ArcGIS）や PROJ ライブラリで実務計算するときは、必ず楕円体パラメータを EPSG コード（例: EPSG:2154 = Lambert-93）で指定してください。逆に概念理解や授業デモのときは、本ツールのように球モデルでも全く問題ありません。

最後に、「φ=90°（極）で公式が破綻する」こと。LCC の式には tan(π/4 + φ/2) が含まれ、φ→90° で tan が無限大、ρ→0 になります。本ツールでは緯度範囲を ±89° に制限してこの特異点を避けていますが、もしコードで自前実装する場合、極点周辺で NaN や Infinity が発生するので必ず緯度上限を設定してください。同様に、両標準緯線が極めて近い（φ₁≈φ₂、差 0.1° 未満）と n の分母が 0 に近づき不安定になるので、接円錐の場合は分岐して n=sin(φ₁) を直接代入する実装にします（このツールでもそうしています）。

使い方ガイド

  1. 第1標準緯線φ₁（北側、例：40°N）と第2標準緯線φ₂（南側、例：30°N）を度単位で入力します。この2本の緯線上ではスケール係数k=1.0となります。
  2. 中心子午線λ₀（例：105°W、米国SPCS zone 14の場合）と基準緯度φ₀（投影原点、通常φ₁とφ₂の中点）を設定します。
  3. 対象地点の緯度φ、経度λを入力すると、リアルタイムで投影座標X,Y(km単位)、円錐定数n、スケール係数k、歪み(%)、円錐頂点からの距離ρが計算され、歪みカーブグラフに表示されます。

具体的な計算例

米国SPCS(State Plane Coordinate System)ゾーン14（中部テキサス）の設定例：φ₁=32°N、φ₂=30°30'N、λ₀=100°30'W、φ₀=30°N。ダラス地点（32.78°N, 96.81°W）の座標計算では円錐定数n≈0.667、投影座標X≈268,500m、Y≈450,300m、スケール係数k≈1.0003、歪み≈0.03%となります。国際航空図ONC（Operational Navigation Chart）では φ₁=60°N、φ₂=30°N設定で中緯度帯の角度歪みを±1.5%以内に抑えます。

実務での注意点