大地测量 · 坐标正算

坐标正算

已知起点大地坐标、大地方位角与椭球面距离,基于 Vincenty 正算公式精确求解终点坐标。支持 WGS-84、CGCS2000、GRS-80 等主流参考椭球,精度达 ±0.5 mm。

首页/ 地理图志/ 实用工具/ 坐标正算
🧭

坐标正算

Vincenty 正算公式 · 迭代法 · 椭球面大地主题解算

当前椭球:WGS-84 长半轴 a = 6378137.000 m 扁率 f = 1/298.257223563

北纬为正,南纬为负

东经为正,西经为负

0°=正北,90°=正东,180°=正南,270°=正西

两点间的大地线长度

📐 计算结果 精度 ±0.5 mm
终点纬度 B2
终点经度 L2
反方位角 A21
终点归化纬度 U2
大地方线角距 σ
迭代次数

终点纬度 DMS: 终点经度 DMS:

🗺 计算示意图

示意图:椭球面上的大地方线、起点/终点与方位角关系(墨卡托投影示意)

每行格式:B1(°), L1(°), A12(°), S(m) 使用英文逗号分隔

支持直接粘贴 CSV 格式数据。点击"执行批量正算"后,结果将以下方表格形式呈现。

算法说明

Vincenty 正算公式(1975)是大地测量学中解决椭球面大地主题解算的经典迭代算法。与贝塞尔公式相比,Vincenty 公式在短距离与长距离上均具有更高精度,收敛速度快。

本工具实现了完整的 Vincenty 正算流程,迭代收敛阈值设为 1×10⁻¹² 弧度,通常 3–6 次迭代即可收敛。

  • 适用范围:任意两点间距离(< 20 000 km)
  • 精度:±0.5 mm(相对于 GeographicLib)
  • 失效情况:两点几乎对跖(antipodal)时可能不收敛

参考椭球参数

椭球长半轴 a (m)扁率 f
WGS-846378137.01/298.257223563
CGCS20006378137.01/298.257222101
GRS-806378137.01/298.257222101
北京 546378245.01/298.3
西安 806378140.01/298.257
克拉索夫斯基6378245.01/298.3

短半轴 b = a(1 − f),第一偏心率平方 e² = 2f − f²

核心公式

// 归化纬度 U1 = arctan[(1−f)·tan φ1] // 大地方线常数 sin σ = cos U1 · sin α1 cos²σ = 1 − sin²σ // 辅助量 = cos²σ · (a²−b²)/b² A = 1 + u²/16384 · (4096 + u²·(−768 + u²·(320 − 175u²))) B = u²/1024 · (256 + u²·(−128 + u²·(74 − 47u²))) // 迭代求角距 σ σ = s/(b·A) → σ = s/(b·A) + Δσ // 终点坐标 φ2 = arctan₂( sinU1·cosσ + cosU1·sinσ·cosα1 , (1−f)·√[sin²σ + (sinU1·sinσ − cosU1·cosσ·cosα1)²] ) λ = arctan₂( sinσ·sinα1 , cosU1·cosσ − sinU1·sinσ·cosα1 ) L = λ − (1−C)·f·sinσ·[σ + C·sinσ·(cos2σₘ + C·cosσ·(−1+2cos²2σₘ))] λ2 = λ1 + L

完整公式参见 Vincenty, T. (1975). Direct and Inverse Solutions of Geodesics on the Ellipsoid with Application of Nested Equations. Survey Review, 23(176), 88–93.

使用提示

  • 方位角:从正北方向顺时针计量,范围 0°–360°
  • 距离:输入两点间的椭球面大地线长度(非水平投影距离)
  • 对跖点:当起点与终点几乎互为对跖点时,Vincenty 公式可能不收敛
  • 批量模式:支持 CSV 格式粘贴,可一次性计算数百组数据

大地测量知识

深入理解坐标正算的原理、算法与应用场景

📐 什么是坐标正算?

坐标正算(Geodetic Forward Computation)是大地测量学的基本问题之一:已知起点的大地坐标 (B1, L1)、大地方位角 A12 与椭球面距离 S,求解终点的大地坐标 (B2, L2) 及反方位角 A21

与平面几何不同,地球表面是椭球面,因此"直线"在椭球面上表现为大地线(Geodesic)—两点间最短路径。坐标正算就是在大地线上进行"里程丈量"与"方向推算"。

在测绘工程中,坐标正算广泛应用于导线测量、三角网平差、GPS 基线向量解算、航空摄影测量外业控制等领域。

🔬 Vincenty 算法详解

Vincenty 正算公式由波兰大地测量学家 Thaddeus Vincenty 于 1975 年发表,是对贝塞尔(Bessel)大地主题解算公式的重大改进。

其核心思想是:将椭球面上的大地线问题转换为辅助球面上的大圆航线问题,通过引入归化纬度(Reduced Latitude)消除椭球扁率的一阶影响,再利用嵌套方程进行高精度修正。

特性Vincenty贝塞尔
短距离精度±0.5 mm±1 mm
长距离精度±0.5 mm±数米
收敛速度3–6 次迭代5–10 次迭代
对跖点稳定性可能不收敛较稳定

🌍 参考椭球与坐标系

坐标正算的结果依赖于所采用的参考椭球(Reference Ellipsoid)。不同椭球对应不同的长半轴 a 与扁率 f,进而影响大地线的曲率与长度。

WGS-84(World Geodetic System 1984)是目前全球定位系统(GPS)使用的标准椭球,由美国国防部定义并于 1984 年启用,2014 年进行了精化(G1762)。

CGCS2000(China Geodetic Coordinate System 2000)是中国当前法定的国家大地坐标系,其定义与 WGS-84 在厘米级一致,但在框架实现与历元上有所不同,高精度测绘工程必须使用 CGCS2000。

🛰 典型应用场景

导线测量:在野外测量中,从已知控制点出发,测得边长与方位角后,通过坐标正算逐点推算未知点坐标,构成导线网。

GPS 基线解算:在 GNSS 相对定位中,基线向量(ΔX, ΔY, ΔZ)可转换为大地线距离与方位角,再通过坐标正算验证解算结果。

大地控制网平差:在国家一、二等三角网与精密导线网平差中,坐标正算是观测方程构建的核心环节。

惯性导航:航空与航海惯性导航系统通过航向角与里程进行坐标正算,实现自主定位(Dead Reckoning)。

需要坐标反算或其他大地测量工具?

地理实用工具箱持续更新中。坐标反算(已知两点求方位角与距离)、高斯投影正反算、七参数转换等工具即将上线。

坐标转换器 大圆距离计算 返回工具箱