如何在AltAz坐标系下基于天球观测点求解拟合圆的圆心

完全可以不经过平面投影，直接基于AltAz坐标完成天球上共圆点的圆心求解，且球面域原生拟合的精度显著高于投影到平面后做最小二乘的结果，不会引入投影畸变误差。

• 首先要避开一个常见误区：不要直接把方位角当X轴、高度角当Y轴，套用平面圆的最小二乘公式。AltAz是球面曲面上的曲线坐标，天球上的小圆展开到A-h二维坐标网格上不是正圆，当小圆的角半径较大时，形变会非常明显，直接套平面公式得到的结果完全没有参考价值。

具体实现方法 #

最简洁稳定的直接拟合方案不需要写复杂的球面三角函数方程，按以下步骤操作即可：

1. 坐标转换：把所有观测点的AltAz坐标转成单位天球对应的三维单位向量，转换时注意统一角度单位为弧度，方位角从北点向东计量、高度角从地平向上计量，转换逻辑为：

   # 单点点坐标转换逻辑
   x = cos(高度角) * sin(方位角)
   y = cos(高度角) * cos(方位角)
   z = sin(高度角)
   

2. 构造约束方程：天球上同一个小圆上的所有点，到圆心的角距为固定值（即小圆的球面角半径θ）。设待求圆心的单位向量为c = [cx, cy, cz]，任意观测点的单位向量为v_i = [xi, yi, zi]，二者的点积恒等于角距的余弦值：

   xi*cx + yi*cy + zi*cz = cosθ
   

3. 最小二乘求解：把cx, cy, cz, cosθ作为待求参数，给6个观测点分别列方程组成超定方程组，加上单位向量约束cx² + cy² + cz² = 1，用带约束的线性最小二乘即可直接解出圆心向量，再反算回AltAz坐标就是你要的小圆天球坐标。
   你提到所有点位都在圆的同一侧，这个先验信息刚好可以排除球面拟合常见的对径点二义性问题，得到的解是唯一稳定的。

方案对比说明 #

• 如果你的6个点覆盖的天区范围很小（角半径小于10度），把点投影到切平面做XY平面圆拟合的误差在可接受范围内；如果覆盖天区更大，或者需要做望远镜指向校准、天文定位这类高精度场景，球面原生拟合的误差会比平面投影拟合低1~2个数量级。
• 如果你坚持要完全在AltAz二维曲线坐标下做拟合，也可以把球面小圆方程展开为方位角、高度角的三角函数形式，构造非线性最小二乘目标函数迭代求解，但这种方法数值稳定性差、推导繁琐，实际工程中没人会用，和上述单位向量法的数学本质完全等价，没有任何优势。

注意：求解时不要忘了给圆心向量加单位模长约束，不加约束解出来的向量不在天球面上，结果是错的。

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