[論文レビュー] 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ソフトウェアシステムは、シミュレートされたデータ上で実現可能性を示し、1マイクロアーセコンド未満の精度を達成し、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億の未知数に対するグローバル最小二乗調整として天測コアソリューションを定式化する。
• 反復的ソリューションプロセスを、源、姿勢、キャリブレーション、グローバル更新という明確なブロックに分割し、フレーム回転やプライマリーストア選択といった補助プロセスを含む。
• 効率的かつ安定した衛星姿勢表現のためのクaternionを用い、姿勢およびキャリブレーションパラメータの滑らかな時系列補間にBスプラインを適用する。
• 最小二乗問題から生じる正定値半定値の正規方程式系を解くために、修正されたコレスキー法を適用する。
• 初期推定値が最終解に近く、反復回数を削減できるように、段階的改善を可能にする構造をソリューションに組み込む。
• 200万個のプライマリーストアを想定したシミュレートデータを用いて、数値的安定性と収束性を検証する。

実験結果

リサーチクエスチョン

• RQ1 約5億の天測および機器パラメータを同時に推定するグローバル最小二乗ソリューションを、効率的かつ堅牢に実装する方法は何か？
• RQ2 1億個のプライマリーストアに対して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のすべての基本的コンponentsが整っており、実世界のデータの複雑さを考慮しても、ガイヤーの野心的な天測精度目標を達成できると予想される。