[論文レビュー] Flat-relative optimal extraction. A quick and efficient algorithm for stabilised spectrographs

M. Zechmeister, G. Anglada‐Escudé+1 more|GoeScholar The Publication Server of the Georg-August-Universität Göttingen (Georg-August-Universität Göttingen)|Nov 20, 2013

Astronomy and Astrophysical Research参考文献 13被引用数 9

ひとこと要約

本稿では、安定化されたファイバー入射分光計向けに高速かつ効率的なアルゴリズムであるFlat-relative Optimal Extraction (FOX) を提案する。この手法は、高SNRのマスターフラットを空間プロファイルマスクとして用い、最適1次元スペクトル抽出を実現する。フラットフィールド基準でのスペクトル抽出により、機器プロファイルのモデル化や別個のフラットフィールド補正の必要性が排除され、計算コストを著しく低減しつつ、標準的な最適抽出と同等の性能を達成し、PSFの非対称性に対しても頑健である。

ABSTRACT

Optimal extraction is a key step in processing the raw images of spectra as registered by two-dimensional detector arrays to a one-dimensional format. Previously reported algorithms reconstruct models for a mean one-dimensional spatial profile to assist a properly weighted extraction. We outline a simple optimal extraction algorithm including error propagation, which is very suitable for stabilised, fibre-fed spectrographs and does not model the spatial profile shape. A high signal-to-noise, master-flat image serves as reference image and is directly used as an extraction profile mask. Each extracted spectral value is the scaling factor relative to the cross-section of the unnormalised master-flat which contains all information about the spatial profile as well as pixel-to-pixel variations, fringing, and blaze. The extracted spectrum is measured relative to the flat spectrum. Using echelle spectra of the HARPS spectrograph we demonstrate a competitive extraction performance in terms of signal-to-noise and show that extracted spectra can be used for high precision radial velocity measurement. Pre- or post-flat-fielding of the data is not necessary, since all spectrograph inefficiencies inherent to the extraction mask are automatically accounted for. Also the reconstruction of the mean spatial profile by models is not needed, thereby reducing the number of operations to extract spectra. Flat-relative optimal extraction is a simple, efficient, and robust method that can be applied easily to stabilised, fibre-fed spectrographs.

研究の動機と目的

• 高精度の径速度測定に用いられる安定化・ファイバー入射分光計向けに、より高速かつ頑健なスペクトル抽出手法の開発。
• 最適抽出における空間プロファイル（多項式やガウス関数による）のモデリングを不要にすることで、低SNR領域での誤差に起因する計算コストの増大と誤差の発生を回避。
• マスターフラットを抽出基準として用いることで、フラットフィールディング、マスキング、ブレイズ補正を1ステップで統合し、フラットフィールドによる除算に起因する数値的不安定性を回避。
• HARPS や CARMENES などの大規模な調査において、1スペクトルあたりの演算数を削減することで、計算効率を向上。
• 標準的な最適抽出と同等の高SNRおよび径速度精度を維持していることを検証すること。

提案手法

• 高SNRのマスターフラット画像を、ピクセル間の変動、フラaring、ブレイズ関数を含むすべての機器的要因を符号化した参照空間プロファイルマスクとして使用。
• 各抽出スペクトルフラックス値は、そのピクセル列における未正規化マスターフラット断面に対するスケーリング係数として計算され、プロファイルモデリングを伴わない最適抽出を実現。
• 機器点スプル関数（iPSF）のフィッティングや再構築を避ける。代わりに、分光計の時間的安定性に依存し、マスターフラットが科学露出時におけるPSFを正確に表現していることを保証。
• 科学データに対して直接、マスターフラットを空間マスクとして用いて抽出を実行し、抽出、マスキング、ブレイズ補正を1ステップで統合。
• 各抽出フラックス点の信頼性の高い不確実性推定を保証するため、誤差伝搬をアルゴリズムに組み込む。
• HARPS エッセールスペクトルにこの手法を適用し、HARPS Data Reduction System を用いた標準的な最適抽出と性能を比較。

実験結果

リサーチクエスチョン

• RQ1 機器プロファイルのモデリングを伴わずに、フラットフィールド基準を最適抽出における空間プロファイルマスクとして直接使用可能か？
• RQ2 フラット基準最適抽出は、標準的な最適抽出と同等のSNRおよび径速度精度を達成できるか？
• RQ3 プロファイルフィッティングの排除とフラットフィールドによる除算の回避により、計算コストの低減と頑健性の向上が図れるか？
• RQ4 フラットフィールドの非対称性や非一様性（例：フラリング、ブレイズ効果）が存在する状況下でも、本手法は有効か？
• RQ5 大規模な調査において高速かつ安定したスペクトル抽出を要する状況に適しているか？

主な発見

• FOXは、HARPSデータに適用した際、標準的な最適抽出と区別できないSNR性能を達成しており、スペクトル品質に顕著な劣化は認められない。
• 重複するオーダーからの抽出スペクトルが良好に一致しており、本手法に内蔵されたブレイズ補正の有効性が裏付けられている。
• FOXと標準的なHARPS Data Reduction Systemを用いて抽出した径速度時間系列は、共に同じ程度の散乱（rms）を示しており、系外惑星探査における同等の精度が得られていることが確認された。
• 空間プロファイルのモデリングや別個のフラットフィールディングステップの不要化により、計算コストが顕著に低減された。
• マスターフラットがこれらの影響を内蔵的に捉えているため、PSFの非対称性や非一様性に対しても本手法は頑健である。
• スペクトルオーダーの正確な局在やゆっくり変化する空間プロファイルの必要がないため、多様な分光計構成への柔軟性が向上している。