过两点的大圆轴与球面交点的名称及推导公式咨询

一、名称定义 #

这两个对跖点被称为大圆的极点（Great Circle Poles）。从极点到大圆上任意一点的球面距离恒为90°（π/2弧度），且极点所在轴线与大圆所在平面垂直。

二、交点位置 #

它们是垂直于该大圆所在平面的轴线与球面的两个交点，两点互为对跖点（即两点关于球心对称，球面距离为180°）。

三、公式推导步骤 #

假设球面为单位球（半径不影响经纬度结果），已知两点经纬度：点1(φ₁, λ₁)、点2(φ₂, λ₂)，其中φ为纬度（北纬取正、南纬取负），λ为经度（东经取正、西经取负）。

1. 转换为三维笛卡尔坐标 #

将经纬度转换为单位球面上的三维向量：

点1坐标：
x₁ = cosφ₁ · cosλ₁
y₁ = cosφ₁ · sinλ₁
z₁ = sinφ₁

点2坐标：
x₂ = cosφ₂ · cosλ₂
y₂ = cosφ₂ · sinλ₂
z₂ = sinφ₂

2. 计算极点方向向量 #

两个点向量的叉乘结果就是垂直于大圆平面的向量，对应其中一个极点的方向：

V = P₁ × P₂ = (
  y₁z₂ - y₂z₁,
  z₁x₂ - z₂x₁,
  x₁y₂ - x₂y₁
)

其反方向向量-V对应另一个对跖极点的方向。

3. 转换回经纬度 #

将极点的笛卡尔坐标(x, y, z)转换为经纬度：

纬度φ_p = arcsin(z)
经度λ_p = arctan2(y, x)

• arctan2(y, x)的结果范围是[-π, π]，可转换为[-180°, 180°]或[0°, 360°]的经度格式；
• 若x=0且y=0，说明极点在南北极，经度可任意取值。

特殊情况说明 #

若两点互为对跖点，则过它们的大圆有无穷多个，此时极点也有无穷多个，无法通过上述方法得到唯一解。

内容的提问来源于stack exchange，提问作者user258279