[論文レビュー] Introducing ADAPTSMOOTH, a new code for the adaptive smoothing of astronomical images

S. Zibetti|arXiv (Cornell University)|Nov 25, 2009

Image Enhancement Techniques被引用数 9

ひとこと要約

本論文は、局所的なS/Nに応じて画素ごとのカーネルサイズを変化させることで信号対雑音比（S/N）を最適化する、新しいスタンドアロンの宇宙画像用の適応平滑化コード adaptsmooth を紹介する。この手法は、既存の平均ベースの手法と比較して、フラックス保存性を著しく向上させ、アーティファクトを低減する独自の中央値平均化モードを実装している。銀河の光度測定テストでは、表面輝度のずれが<0.01 mag、フラックス誤差が<0.001 magに抑えられている。

ABSTRACT

We introduce and publicly release a new code, ADAPTSMOOTH, which serves to smooth astronomical images in an adaptive fashion, in order to enhance the signal-to-noise ratio (S/N). The adaptive smoothing scheme allows to take full advantage of the spatially resolved photometric information contained in an image in that at any location the minimal smoothing is applied to reach the requested S/N. Support is given to match more images on the same smoothing length, such that proper estimates of local colours can be done, with a big potential impact on multi-wavelength studies of extended sources (galaxies, nebulae). Different modes to estimate local S/N are provided. In addition to classical arithmetic-mean averaging mode, the code can operate in median averaging mode, resulting in a significant enhancement of the final image quality and very accurate flux conservation. To this goal also other code options are implemented and discussed in this paper. Finally, we analyze in great detail the effect of the adaptive smoothing on galaxy photometry, in particular in terms of surface brightness (SB) profiles and aperture photometry: deviations in SB with respect to the original image can be limited to &lt;0.01 mag, with flux difference in apertures of less than 0.001 mag.

研究の動機と目的

• 光学および近赤外宇宙画像向けに、柔軟でスタンドアロンの適応平滑化コードを開発すること。
• 表面輝度動的範囲が広い拡張源において、空間分解能と信号対雑音比（S/N）のバランスをとる課題に取り組むこと。
• 特に明るい源がかすかな背景に隣接する領域において、フラックス保存性を向上させ、アーティファクトを低減すること。
• ユーザーが多様なノイズ環境に適応できるように、複数のS/N推定モードと平滑化レベルのカットオプションを提供すること。
• 平滑化マスクを用いることで、空間的に分解された色測定を可能とし、正確な多波長光度測定を実現すること。

提案手法

• 適応平滑化は、ユーザーが定めた最小S/Nに達するよう、画素ごとのカーネルサイズを変化させることで実施され、信号が強い領域では最小限の平滑化が保証される。
• コードは2つの平均化モードをサポートする：算術平均（標準）と中央値平均化。後者はフラックス保存誤差とアーティファクトを低減することが示されている。
• 局所S/Nは3つの方法で推定される：ポisson雑音＋バックグラウンド雑音、バックグラウンド雑音のみ、または画素強度のr.m.s.（局所的信号フラクチュエーション）。
• 平滑化マスクが実装されており、多バンド画像処理を可能とし、拡張源における空間的に分解された色測定を可能としている。
• 平滑化レベルカットは、後処理オプションとして導入され、特に高S/Nと低S/N領域の境界付近でのアーティファクトを抑制する。
• 可能な限り空間分解能を保持し、高S/N領域での過剰な平滑化を避けることで、情報損失を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 表面輝度の動的範囲が広い宇宙画像において、表面輝度と全フラックスを保存するように適応平滑化を最適化する方法は何か？
• RQ2 フラックス保存誤差とアーティファクトを最小限に抑えるために、中央値平均化と算術平均化の相対的利点は何か？
• RQ3 ポisson＋バックグラウンド雑音、バックグラウンド雑音のみ、局所r.m.s. の各S/N推定法が、多様な画像環境における適応平滑化の性能と頑健性に与える影響は何か？
• RQ4 平滑化レベルカットは、高S/Nと低S/N領域の過渡領域で、フラックス保存誤差とアーティファクトをどの程度低減できるか？
• RQ5 適応平滑化により、銀河や星雲の画像において、表面輝度の測定可能な動的範囲を著しく拡張できるか、かつ光度測定の正確性を維持できるか？

主な発見

• adaptsmooth の中央値平均化モードは、算術平均化と比較して、フラックス保存誤差を約50％低減し、平滑化レベル1を除いた場合、フラックス成長曲線のずれが0.001 mag未満に抑えられている。
• 高コントラスト領域でさえも、平滑化画像の表面輝度プロファイルは元の画像と0.01 mag未満のずれで保たれている（方位平均による）。
• 平滑化レベル1のカットを併用した中央値平滑化が、未平滑化領域と平滑化領域の境界で生じるバイアスを最小限に抑え、最も正確なフラックス保存が達成される。
• ポisson＋バックグラウンド、バックグラウンドのみ、局所r.m.s. の3つのS/N推定法をサポートするコードの柔軟性により、ノイズモデルの事前知識が限られた状況でも頑健な性能が得られる。
• adaptsmooth を用いた適応平滑化により、空間分解能と光度測定の正確性を保持したまま、多バンド画像処理で局所的な色測定が信頼性を持って可能になる。
• csmooth や ftools などの既存コードと比較して、adaptsmooth は中央値平均化とより柔軟な平滑化レベルカットを独自に提供しており、光学および近赤外データにおける画像品質と光度測定の正確性を顕著に向上させている。