[論文レビュー] The correct estimate of the probability of false detection of the matched filter in the detection of weak signals. II. (Further results with application to a set of ALMA and ATCA data)

R. Vio, C. Vergès+1 more|arXiv (Cornell University)|May 9, 2017

Soil Geostatistics and Mapping参考文献 15被引用数 4

ひとこと要約

本稿では、天体画像における弱い信号検出の信頼性を評価するための、標準的な誤検出確率（PFA）よりもより正確な指標として、特定誤警報確率（SPFA）を導入する。マッチドフィルタを用いた検出において、信号位置が不明な場合にPFAが誤検出リスクを著しく低く見積もる問題を解消する。SPFAは、等方的ガウスノイズ場のピーク分布を活用し、検出された源が誤検出である確率を定量的に評価する。ALMAおよびATCAのデータにおいて、SPFAは一部の検出でPFAが非常に低いにもかかわらず、信頼度が43–50%にとどまることを示しており、PFAによる過剰な有意水準の見積もりを是正する。

ABSTRACT

The matched filter (MF) is one of the most popular and reliable techniques to the detect signals of known structure and amplitude smaller than the level of the contaminating noise. Under the assumption of stationary Gaussian noise, MF maximizes the probability of detection subject to a constant probability of false detection or false alarm (PFA). This property relies upon a priori knowledge of the position of the searched signals, which is usually not available. Recently, it has been shown that when applied in its standard form, MF may severely underestimate the PFA. As a consequence the statistical significance of features that belong to noise is overestimated and the resulting detections are actually spurious. For this reason, an alternative method of computing the PFA has been proposed that is based on the probability density function (PDF) of the peaks of an isotropic Gaussian random field. In this paper we further develop this method. In particular, we discuss the statistical meaning of the PFA and show that, although useful as a preliminary step in a detection procedure, it is not able to quantify the actual reliability of a specific detection. For this reason, a new quantity is introduced called the specific probability of false alarm (SPFA), which is able to carry out this computation. We show how this method works in targeted simulations and apply it to a few interferometric maps taken with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Australia Telescope Compact Array (ATCA). We select a few potential new point sources and assign an accurate detection reliability to these sources.

研究の動機と目的

• 信号位置が不明な状況において、標準的なPFAがマッチドフィルタ検出の誤警報リスクを著しく低く見積もるという深刻な問題に対処すること。
• 弱信号状況における個々の検出の信頼性を定量化するため、PFAよりもより正確な統計的指標を開発すること。
• 等方的ガウス過程のピーク分布に基づく新しい手法を、ALMAおよびATCAの実際の干渉計測定データに適用し、検証すること。
• 天体マップにおける新たな点源候補の検出信頼度を、堅牢で定量的な評価で提供すること。

提案手法

• 標準的なPFAを補正するため、ピークの確率密度関数（PDF）を等方的ガウスノイズ場の性質に基づいてモデル化し、ピーク振幅のガウス分布を仮定するのではなく、ピークの分布を直接扱う。
• 特定誤警報確率（SPFA）を導入し、これは与えられたピークがノイズ由来である確率として定義され、フィルタドノイズ場の最大値のPDFを用いて計算される。
• 等方的ガウス過程の確率的理論に基づき、しきい値を超えるピークの期待数を用いてSPFAを導出する。
• 空間相関と非一様ノイズを考慮し、有効自由度を推定し、実測相関関数を用いてピークの等方性と独立性を検証する。
• シミュレーテッドマップおよびALMA（2バンド）およびATCAの実際の干渉計データに、明るい源のマスキング、平均0、分散1への正規化を施して適用する。
• 所望のPFAレベルに基づいてしきい値を設定し、各候補源に対してSPFA値を計算することで、検出信頼度を評価する。

実験結果

リサーチクエスチョン

• RQ1 信号位置が不明な状況において、標準的なマッチドフィルタ手法がなぜ誤検出確率を著しく低く見積もるのか。
• RQ2 標準的なPFAを超えて、特定の検出の統計的有意水準をどのように正確に定量化できるか。
• RQ3 非一様ノイズおよびグリッド化アーティファクトを伴う実宇宙データにおいて、SPFA手法が検出信頼度の推定をどの程度改善するか。
• RQ4 SPFAは、ALMAおよびATCAの干渉計マップにおいて、真の弱い源を的確に同定し、誤検出を効果的に除外できるか。

主な発見

• 標準的なPFA手法は誤警報リスクを著しく低く見積もっており、過剰な検出信頼度の主張を引き起こしている。例えば、PFA ≈ 7.4×10⁻⁴ の源では、SPFA値が検出の信頼度を50%にまで低下させている。
• ALMAデータでは、3つのマップにわたって7つの新たな点源候補が同定され、SPFA値は0.01から0.20の範囲にあり、誤検出リスクが低いことが示された。
• ATCAマップでは、11個の候補源が検出され、PFA ≈ 1.25×10⁻⁴ であり、SPFA値は0.01から0.20の範囲にあり、誤検出が存在しない確率は60%、最大で1件の誤検出である確率は93%であった。
• SPFA値はピークの経験的分布と整合しており、ヒストограм分析および自己相関関数を用いた検証により、非一様なuvカバレッジおよびグリッド化アーティファクトが存在する状況でも、本手法の堅牢性が確認された。
• ATCAマップにおいて既知の源（例：SSTSL2源）が正しく同定されたことから、実運用環境下でも本手法の信頼性が裏付けられた。
• SPFA指標は、誤差が数パーセントの範囲に収まる精度の高い、定量的で明確な検出信頼度の測定を可能にし、後続観測計画の立案に適している。