[论文解读] The astrometric core solution for the Gaia mission. Overview of models, algorithms and software implementation

L. Lindegren, U. Lammers+4 more|Dec 18, 2011

Geophysics and Gravity Measurements参考文献 41被引用 125

一句话总结

本文提出了盖亚天体测量核心解的数学与算法框架，这是一种对约5亿个参数（包括位置、视差、自行、卫星姿态及仪器校准）的全局最小二乘估计。AGIS软件系统在模拟数据上成功验证了其可行性，实现了亚微角秒级精度，并在10 Tflop/s的计算设备上，用约60天完成真实任务规模的计算。

ABSTRACT

The Gaia satellite will observe about one billion stars and other point-like sources. The astrometric core solution will determine the astrometric parameters (position, parallax, and proper motion) for a subset of these sources, using a global solution approach which must also include a large number of parameters for the satellite attitude and optical instrument. The accurate and efficient implementation of this solution is an extremely demanding task, but crucial for the outcome of the mission. We provide a comprehensive overview of the mathematical and physical models applicable to this solution, as well as its numerical and algorithmic framework. The astrometric core solution is a simultaneous least-squares estimation of about half a billion parameters, including the astrometric parameters for some 100 million well-behaved so-called primary sources. The global nature of the solution requires an iterative approach, which can be broken down into a small number of distinct processing blocks (source, attitude, calibration and global updating) and auxiliary processes (including the frame rotator and selection of primary sources). We describe each of these processes in some detail, formulate the underlying models, from which the observation equations are derived, and outline the adopted numerical solution methods with due consideration of robustness and the structure of the resulting system of equations. Appendices provide brief introductions to some important mathematical tools (quaternions and B-splines for the attitude representation, and a modified Cholesky algorithm for positive semidefinite problems) and discuss some complications expected in the real mission data.

研究动机与目标

• 解决从盖亚任务处理约10亿条天体测量观测数据以生成高精度天球参考系的挑战。
• 开发一种全局迭代解法，同时估计天体测量参数、卫星姿态及仪器校准参数。
• 确保在处理高达1亿颗主星时，具备数值鲁棒性与计算可行性。
• 设计一种可扩展的软件系统，能够处理盖亚数据中庞大的数据量与复杂的相互依赖关系。
• 为天球参考系与仪器校准提供基础，这对盖亚任务后续所有数据处理至关重要。

提出的方法

• 将天体测量核心解表述为对约5亿个未知数的全局最小二乘平差，包括源位置、视差、自行、姿态及校准参数。
• 实施一种分块迭代求解过程，包括源、姿态、校准与全局更新等独立模块，并辅以坐标系旋转与主星选择等辅助过程。
• 采用四元数表示卫星姿态，以实现高效且稳定的计算；使用B样条对姿态与校准参数进行平滑的时间插值。
• 应用改进的Cholesky算法求解由最小二乘问题产生的正半定法方程组。
• 设计解法以支持逐步优化，从接近最终解的先验估计出发，从而减少所需迭代次数。
• 利用200万颗主星的模拟数据验证方法，证明其具备数值稳定性与收敛性。

实验结果

研究问题

• RQ1 如何高效且稳健地实现对5亿个天体测量与仪器参数的全局最小二乘解？
• RQ2 为确保对1亿颗主星实现亚微角秒级精度的天体测量解，需要哪些数学模型与数值算法？
• RQ3 如何在可行的时间与资源约束下管理盖亚海量数据流的计算复杂度？
• RQ4 四元数与B样条在实现卫星姿态与校准随时间稳定且精确表示中起到何种作用？
• RQ5 在盖亚任务规模下，从近似初始估计出发，该迭代解法在多大程度上能可靠收敛？

主要发现

• AGIS软件系统成功实现了全局迭代解法，在200万颗主星的模拟数据上验证了其数值有效性与鲁棒性。
• 姿态参数γ的RMS不确定性约为3×10⁻⁵，与理论误差估计及真实参数值一致。
• 在14节点IBM集群（0.65 Tflop/s）上运行135次迭代耗时近6天，对应约3×10¹⁷次浮点运算。
• 按全任务规模（10⁸颗主星）估算，总计算工作量约为1.5×10¹⁹ flop，需在10 Tflop/s机器上运行约60天。
• 由于从接近最终解的初始估计出发，解法具备可行性，可能显著减少所需迭代次数。
• AGIS的所有核心组件均已就绪，系统预计能应对真实数据的复杂性，实现盖亚任务雄心勃勃的天体测量精度目标。