兰伯特等角圆锥投影模拟器 返回
地图学与 GIS

兰伯特等角圆锥投影模拟器

用一个与地球相割于两条标准纬线的圆锥将地表投影到平面,这就是 LCC(Lambert Conformal Conic)投影。调整标准纬线、中央子午线和参考纬度,即可实时观察圆锥常数 n、平面坐标 (x, y)、比例尺因子与畸变率。这是美国 SPCS、国际航空图 ICAO ONC 与气象预报图所使用的中纬度主力投影法。

参数设置
标准纬线 1 φ₁
°
圆锥与地球相交的低纬度标准纬线
标准纬线 2 φ₂
°
高纬度标准纬线。与 φ₁ 相等时为切圆锥
中央子午线 λ₀
°
投影中央经度。SPCS 中各州取靠近州中心的值
参考纬度 φ₀
°
投影平面原点对应的纬度(伪北距起点)
投影点纬度 φ
°
要投影到平面的地表点纬度
投影点经度 λ
°
要投影到平面的地表点经度
地球半径 R
km
球模型地球半径。WGS84 平均半径为 6371 km
计算结果
圆锥常数 n
平面坐标 X (km)
平面坐标 Y (km)
比例尺因子 k
畸变率 (%)
圆锥顶点距离 (km)
圆锥投影几何示意

地球与一个沿两条标准纬线(φ₁、φ₂)相割的圆锥相交。投影点(红色)被映射到圆锥面,将圆锥展开后即得到下方的平面坐标 (x, y)。

比例尺因子随纬度变化(畸变曲线)
投影平面上的经纬网
理论与主要公式

$$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 用按带划分的横轴圆柱,而兰伯特则用一个沿中纬度两条纬线(φ₁、φ₂)与地球相割的圆锥。美国 SPCS、国际航空图 ICAO ONC 与气象预报图都把它作为标准投影。把滑块保持在 φ₁=33°、φ₂=45° 时,n 大约是 0.6305——这就是圆锥的开口程度。
🙋
n=0.6305 是不是意味着圆锥展开后只占整圆的约 63%?
🎓
正是如此。把圆锥展开为平面时,地球 360° 的经度对应平面上 n×360° ≈ 227° 的扇形。看极限就一目了然:让 φ₁ 趋向 0,n 也趋向 0,圆锥几乎变成圆柱(墨卡托极限);让 φ₁、φ₂ 都接近 89°,n 趋近 1,圆锥几乎变成平面(极球面投影极限)。兰伯特正是介于二者之间、为中纬度优化的等角投影。
🙋
"等角"是说保持角度吧?那为什么还会有畸变呢?
🎓
观察得好。"等角(conformal)"指任意一点周围的局部角度不变,即微小圆映射后仍为圆(不会变成椭圆)。但大小会随纬度变化,比例尺因子 k 描述的正是这一点。两条标准纬线(φ₁、φ₂)之间 k<1(压缩),外侧 k>1(拉伸)。看下方的畸变曲线就能发现 φ=40° 时 k≈0.994,即压缩约 0.5%。
🙋
原来如此,把目标夹在两条标准纬线之间能减小畸变。SPCS 是怎么设计的呢?
🎓
各州把 φ₁ 和 φ₂ 设在该州南北端向内约 1/6 的位置。例如科罗拉多州中央带(约北纬 37°–41°)取 φ₁=37°50′、φ₂=39°45′。这样整个州内 k 都在 0.9999~1.0001 之间,距离误差约为万分之一。1:24,000 地形图上连 10 厘米的偏差都不会出现。ICAO ONC(1:1,000,000)也是同样思路:把标准纬线放在每张图的内侧,使 |k−1| 控制在 ±0.05% 以内。
🙋
把 φ₁ 和 φ₂ 设为相同值就成了"切圆锥",它与"割圆锥"有何不同?
🎓
φ₁=φ₂ 的切圆锥只有一条 k=1 的标准纬线,畸变沿一个方向单调增大;而 φ₁≠φ₂ 的割圆锥(secant,这也是本工具默认计算的形式)在两条标准纬线上都有 k=1,二者之间 k<1,外侧 k>1——可以把畸变"从两端向内分摊"。法国旧的 NTF 用切圆锥(φ₀=48°),现代 RGF93 Lambert-93 用割圆锥(φ₁=44°、φ₂=49°)。实际工程几乎全部采用割圆锥。

常见问题

适合东西方向延伸的中纬度地区(大约南北纬 30°–60°)。美国 SPCS(州平面坐标系统)的大多数东西向州、国际航空图 ICAO ONC(1:1,000,000)、各国气象预报图、法国 NTF/RGF93 与加拿大部分地形图都采用 LCC。两条标准纬线(φ₁, φ₂)将圆锥与地球相割,使两条纬线之间比例尺小于 1、外侧大于 1,把目标区域置于两线之间即可把畸变控制在 ±1% 以内。
n 控制投影所用圆锥的开口大小。将圆锥展开为平面后,地球 360° 经度对应平面上 2π·n 的弧角。对切圆锥(φ₁=φ₂)有 n = sin(φ₁):低纬标准纬线对应浅圆锥(n→0 接近圆柱,与墨卡托相似),高纬对应尖圆锥(n→1 接近极球面投影)。LCC 即是介于墨卡托(n=0)与立体投影(n=1)之间的广义等角投影,是中纬度的最优选择。
比例尺因子 k 是地面单位距离对应平面单位距离的比值,k=1 表示无畸变,k=1.02 表示拉伸 2%,k=0.98 表示压缩 2%。在 1:50,000 航空图上 2% 的畸变意味着 10 km 的航段会有 200 m 偏移,可能超过仪表进近的容差。SPCS 与 ONC 通过限定每张图的范围,使 |k−1| 控制在 ±0.05% 以内。本工具将 |畸变| ≤ 2% 标记为黄色警告,> 10% 标记为红色不通过。
二者都是等角投影(保持局部角度),但适用范围不同。墨卡托使圆柱与赤道相切,赤道附近畸变最小,高纬度面积爆炸式放大(格陵兰看起来与非洲一样大)。兰伯特圆锥则在中纬度与地球相割,所以中纬度畸变最小。航空航海实践中:赤道至热带(如东南亚航线)用墨卡托;中纬度(北美、欧洲航空图)用兰伯特;极地航线用极球面投影。

实际应用

国际航空图 ICAO ONC(1:1,000,000):全球约 270 幅图几乎都采用 LCC。每幅图覆盖纬度 4°、经度 6°(赤道附近)的较窄范围,将两条标准纬线放在图幅上下边界向内 1/6 处,可将图幅内畸变控制在 ±0.05% 以内。飞行员用直尺在图上量得的距离基本就是真实距离,因此即便在数字化时代,LCC 依旧是仪表进近与航路规划的世界标准基础。

美国州平面坐标系统(SPCS):东西方向细长的州(加州、得州、科罗拉多等)几乎都使用 LCC 作为州平面投影。例如科罗拉多中央带取 φ₁=37°50′、φ₂=39°45′、λ₀=-105°30′,全州各处 k 都在 0.9999~1.0001。地籍登记、道路设计、GIS 基础底图乃至自动驾驶高精地图都建立在 SPCS 之上。与 NAD83 组合是现代美国工程坐标的事实标准。

气象预报图:世界气象组织 WMO 的中纬度高空天气图、500 hPa 解析图、急流图都使用 LCC 网格。气象现象的等压线与锋面"角度"很重要,因此偏好等角投影;分析对象(西风带)正好位于中纬度,LCC 最理想。NCEP/GFS 与 ECMWF 的数值产品同样以 LCC 网格输出,再由各国气象局分发到地面台站。

法国与加拿大的国家坐标系:法国从早期的四带切圆锥 NTF 过渡到现代 RGF93/Lambert-93(割圆锥,φ₁=44°、φ₂=49°、λ₀=3°、φ₀=46.5°),单一坐标系即可覆盖整个法国本土,是 IGN 地形图、地籍、车载导航、洪水灾害图等各类国家基础设施的基础。加拿大部分航空与地形图同样基于 LCC。GIS 软件中若 .prj 文件写着 PROJCS["NTF (Paris) / Lambert ..."],对应的就是该投影。

常见误解与注意事项

首先,"等角即距离精确"是常见误区。LCC 仅保持局部角度,不保持面积与距离。比例尺因子 k 与纬度相关,例如同一张图 φ=30° 时 k=1.005、φ=50° 时 k=1.008,用平面直尺量两者之间的距离会有约 0.7% 的误差——1,000 km 的航线上偏差约 7 km。任务规划与测量必须按纬度补偿 k。本工具显示的"畸变率(%)"只是投影点处的局部尺度误差,并非长距离累积误差。

其次,不要混淆球模型与椭球模型。本工具为了简化计算采用球形地球(半径 R),但实际生产用的 SPCS、RGF93 LCC 都基于 GRS80/WGS84 椭球(长半轴 a=6378137 m、扁率 1/298.257)。差异随纬度增大,北纬 45° 处两者差约 0.3%。在 GIS 软件(QGIS、ArcGIS)或 PROJ 中实际计算时,请用 EPSG 代码(例如 EPSG:2154 = Lambert-93)指定椭球与投影参数。仅做概念演示或教学时,球模型完全足够。

最后,φ=±90°(极点)处公式发散。LCC 公式中含 tan(π/4 + φ/2),φ→±90° 时该项趋于无穷,ρ→0。本工具把纬度滑块限制在 ±89° 来避免奇点;如果你自行实现,需在极点附近设阈值防止 NaN 与 Infinity。同样,当两条标准纬线极为接近(|φ₁−φ₂|<0.1°)时 n 的分母趋于 0,数值不稳定。建议判断后切换为切圆锥公式 n=sin(φ₁)(本工具已这样实现)。

使用指南

  1. 设定两条标准纬线(phi1、phi2),通常选择投影区域的南北边界,例如美国SPCS选择38°N和40°N
  2. 输入中央子午线(lam0)和参考纬度(phi0),参考纬度一般设为标准纬线的平均值
  3. 点击计算按钮,获取圆锥常数n、投影坐标(X, Y)、比例尺因子k和畸变率,评估投影适用性

具体计算示例

以法国RGF93为例,设标准纬线phi1=44°N、phi2=49°N,中央子午线lam0=3°E,参考纬度phi0=46.5°N,地点坐标为48°N/2°E。计算结果:圆锥常数n≈0.755,投影坐标X≈923.4 km、Y≈1847.2 km,比例尺因子k≈0.9996,畸变率约-0.04%,圆锥顶点距原点约2156 km。该配置下投影畸变极小,适合大比例尺地形图制作。

实务注意事项

  1. 美国SPCS采用州级LCC,如加州选择phi1=39°20'N、phi2=38°03'N,单州畸变控制在0.01%以内,避免使用全国统一投影导致的大畸变
  2. 国际航空图ICAO ONC采用phi1=10°N、phi2=60°N的宽带设置,必须容忍±0.3%的畸变以覆盖大陆区域
  3. 气象预报图多选择phi1=phi2(退化为单标准纬线),简化计算但应补偿中纬度的形状畸变
  4. 投影坐标Y值过大(超过5000 km)表示圆锥收敛角过小,需调整标准纬线间距