[論文レビュー] A 3D Drizzle Algorithm for JWST and Practical Application to the MIRI Medium Resolution Spectrometer

David R. Law, J. Morrison+5 more|arXiv (Cornell University)|Jun 8, 2023

Atmospheric Ozone and ClimateEarth and Planetary Sciences被引用数 3

ひとこと要約

本論文は、JWSTのMIRI中分解能スペクトロメータ（MIRI MRS）およびNIRSpec IFUデータからスペクトルキューブを構築するための3次元ドリズルアルゴリズムを提示する。この手法は、検出器ピクセルとキューブボクセル間の体積的重複度を計算することで、空間的および分光的重複度の計算を分離し、計算コストを低減する。このアプローチにより、わずか0.03%のスペクトロフィメトリック忠実度の損失で、近似的に理想的な性能が達成され、パイプライン実装が効率的に行える。主な発見として、4ポイントのディザーと1.5×PSF FWHM以上のアパーチャ半径が、リサンプリングアーティファクトを1%未満に抑えることが判明した。また、正確な分光的分散推定には1.5～3の共分散補正係数が必要である。

ABSTRACT

We describe an algorithm for application of the classic `drizzle' technique to produce 3d spectral cubes using data obtained from the slicer-type integral field unit (IFU) spectrometers on board the James Webb Space Telescope. This algorithm relies upon the computation of overlapping volume elements (composed of two spatial dimensions and one spectral dimension) between the 2d detector pixels and the 3d data cube voxels, and is greatly simplified by treating the spatial and spectral overlaps separately at the cost of just 0.03% in spectrophotometric fidelity. We provide a matrix-based formalism for the computation of spectral radiance, variance, and covariance from arbitrarily dithered data and comment on the performance of this algorithm for the Mid-Infrared Instrument's Medium Resolution IFU Spectrometer (MIRI MRS). We derive a series of simplified scaling relations to account for covariance between cube spaxels in spectra extracted from such cubes, finding multiplicative factors ranging from 1.5 to 3 depending on the wavelength range and kind of data cubes produced. Finally, we discuss how undersampling produces periodic amplitude modulations in the extracted spectra in addition to those naturally produced by fringing within the instrument; reducing such undersampling artifacts below 1% requires a 4-point dithering strategy and spectral extraction radii of 1.5 times the PSF FWHM or greater.

研究の動機と目的

• JWSTのスライサー型IFU分光計用に、計算効率の高い3次元ドリズルアルゴリズムを開発すること。特にMIRI MRSおよびNIRSpecを対象とする。
• ディザーデータと不規則にサンプリングされた2次元検出器データから3次元スペクトルキューブを構築する際の、スペクトロフィメトリック忠実度と計算の実行可能性という課題に取り組むこと。
• 最終データキューブにおけるスパツエル間の共分散を定量化し、補正すること。これは、抽出された1次元スペクトルにおける分散推定を歪める要因である。
• 空間的および分光的アンダーサンプリング（要因約2）がスペクトルキューブ品質に与える影響を分析し、その緩和戦略を提供すること。
• リサンプリングアーティファクトを最小限に抑える観測およびデータ解析の実用的推奨事項を提示すること。

提案手法

• アルゴリズムは、2次元検出器ピクセルと3次元データキューブボクセル間の体積的重複度を計算し、空間的および分光的重複度の計算を別々に処理することで、計算コストを低減する。
• 行列形式を用いて、ディザーオブザーバーションからの分光放射度、分散、共分散を計算し、リサンプリングプロセス全体にわたる正確な統計的伝搬を可能にする。
• パイプライン用途に適したスケーラビリティを実現するため、事前にコンパイルされたC言語実装を採用する。
• 抽出された1次元スペクトルにおけるスパツエル間の共分散を補正するスケーリング要因を導出する。その値は波長およびキューブのサンプリングに応じて1.5～3の範囲をとる。
• 本アルゴリズムは、明るい星16CygBの観測を含む実際のMIRI MRSデータを用いて検証され、アーティファクトの振幅とディザー効果の有効性を評価した。
• 周期的アンダーサンプリングアーティファクトを低減するために、4ポイントのディザー戦略と、1.5×PSF FWHM以上のアパーチャ半径のしきい値を組み込む。

実験結果

リサーチクエスチョン

• RQ1 クラシックな2次元ドリズル技法を、JWSTのMIRI MRSおよびNIRSpec機器用のIFU分光計に3次元に拡張する方法は何か？ また、スペクトロフィメトリック忠実度の損失を最小限に抑えるには？
• RQ2 空間的および分光的アンダーサンプリングがスペクトルキューブ品質に与える影響は何か？ また、観測戦略によってどのように緩和できるか？
• RQ3 最終データキューブにおけるスパツエル間の共分散が、抽出された1次元スペクトルの分散に与える影響は何か？ どのような補正が必要か？
• RQ4 リサンプリングアーティファクトを抽出スペクトルで1%未満に抑えるために、どのディザーパターンと分光的アパーチャ半径が必要か？
• RQ5 大規模なJWSTデータ処理に適した3次元ドリズルの効率的パイプライン実装を可能にする計算最適化は何か？

主な発見

• 空間的および分光的重複度の計算を分離することで、3次元ドリズルアルゴリズムはスペクトロフィメトリック忠実度の損失をわずか0.03%に抑え、最小限の精度コストで効率的なパイプライン利用が可能になった。
• MIRI MRSデータキューブにおけるスパツエル間の共分散は、抽出された1次元スペクトルの分散推定を正確に行うために、波長およびアパーチャサイズに応じて1.5～3の補正係数が必要である。
• 緩和策がない場合、重度のアンダーサンプリング（要因約2）により、抽出スペクトルに最大20%の周期的振幅変調が生じる。
• 4ポイントのディザーではリサンプリングアーティファクトが5%未満に抑えられ、1.5×PSF FWHM以上の分光的抽出アパーチャが必要で、これにより1%未満にまで低減できる。
• プレビューキューブ（calwebb_spec2）におけるリサンプリングノイズは、複数のディザー露出を統合した最終キューブ（calwebb_spec3）よりも高い。
• アルゴリズムは行列形式により正確な分散伝搬を実現できるが、実用的な用途では完全な共分散行列は扱いきれない。