[논문 리뷰] Do fast stellar centroiding methods saturate the Cramér-Rao lower bound?

Mohammadjavad Vakili, David W. Hogg|arXiv (Cornell University)|2016. 10. 19.

Adaptive optics and wavefront sensing인용 수 3

한 줄 요약

이 논문은 낮은 신호 대 잡음비, 천체 기반 조건에서 천체측위 정밀도의 크래머-라오 하한(CRLB)에 도달할 수 있는지 평가한다. 구체적으로 다항식 피팅을 통해 PSF 스무딩 또는 가우시안 스무딩을 수행하는 빠른 별 중심점 추정 방법을 검토한다. PSF가 알려져 있을 경우 PSF 매칭 다항식 중심점 추정법은 CRLB에 거의 도달하며, PSF가 알려져 있지 않을 경우 고정된 가우시안 스무딩과 다항식 피팅 조합이 매우 유사한 성능을 보이며, 중심점-빛 방법은 크게 실패함을 확인한다.

ABSTRACT

One of the most demanding tasks in astronomical image processing---in terms of precision---is the centroiding of stars. Upcoming large surveys are going to take images of billions of point sources, including many faint stars, with short exposure times. Real-time estimation of the centroids of stars is crucial for real-time PSF estimation, and maximal precision is required for measurements of proper motion. The fundamental Cramér-Rao lower bound sets a limit on the root-mean-squared-error achievable by optimal estimators. In this work, we aim to compare the performance of various centroiding methods, in terms of saturating the bound, when they are applied to relatively low signal-to-noise ratio unsaturated stars assuming zero-mean constant Gaussian noise. In order to make this comparison, we present the ratio of the root-mean-squared-errors of these estimators to their corresponding Cramér-Rao bound as a function of the signal-to-noise ratio and the full-width at half-maximum of faint stars. We discuss two general circumstances in centroiding of faint stars: (i) when we have a good estimate of the PSF, (ii) when we do not know the PSF. In the case that we know the PSF, we show that a fast polynomial centroiding after smoothing the image by the PSF can be as efficient as the maximum-likelihood estimator at saturating the bound. In the case that we do not know the PSF, we demonstrate that although polynomial centroiding is not as optimal as PSF profile fitting, it comes very close to saturating the Cramér-Rao lower bound in a wide range of conditions. We also show that the moment-based method of center-of-light never comes close to saturating the bound, and thus it does not deliver reliable estimates of centroids.

연구 동기 및 목표

• 천체측위에서 이론적 정밀도 한계인 크래머-라오 하한(CRLB)에 도달할 수 있는지 평가하기 위해 빠르고 계산 효율성이 높은 중심점 추정 방법의 성능을 평가한다.
• 실제로 낮은 신호 대 잡음비(SNR) 및 변동하는 반폭하강점(FWHM) 조건에서 최적의 최대우도 추정기와 실용적인 빠른 방법 간의 성능 격차를 조사한다.
• PSF가 알려져 있을 때와 알려져 있지 않을 때의 두 상황에서 중심점 추정 방법을 비교하며, 천체 기반, 가우시안 잡음이 있는 이미지를 대상으로 한다.
• 다양한 SNR 및 FWHM 값 범위에서 다양한 방법의 제곱근 평균 제곱오차(RMSE)가 CRLB에 비해 얼마나 되는지 정량화한다.
• 대규모 설문 조사에서 고정밀 천체측위에 널리 사용되는 중심점-빛 중심점 추정법이 신뢰할 수 있는지 판단한다.

제안 방법

• 특정 분석적 PSF 모델(Moffat 유사)을 사용하여 상호 간섭이 없는 희미한 별이 포함된 천체 기반, 상관관계가 없는 가우시안 잡음 이미지를 시뮬레이션한다.
• PSF 모델에서 유도된 피셔 정보 행렬을 사용하여 SNR 및 FWHM에 따라 중심점 추정 오차의 이론적 크래머-라오 하한(CRLB)을 계산한다.
• 네 가지 중심점 추정 방법을 적용한다: (1) 진짜 PSF에 대한 최대우도 피팅, (2) PSF 매칭 스무딩 후 3×3 패치에 대한 2차 다항식 피팅, (3) 고정된 가우시안 스무딩 후 다항식 피팅(PSF 지식 없음), (4) 중심점-빛(1차 모멘트) 추정.
• 다양한 시뮬레이션 별에 대해 각 방법의 중심점 추정치의 제곱근 평균 제곱오차(RMSE)를 측정하고, CRLB와의 비율인 RMSE / CRLB를 통해 비교한다.
• 두 가지 시뮬레이션 세트를 수행한다: 하나는 고정된 FWHM에서 SNR를 변화시키고, 다른 하나는 고정된 SNR에서 FWHM를 변화시켜 조건에 대한 강건성을 평가한다.
• 산점도와 중앙값 비율을 사용하여 각 방법이 매개변수 공간 전역에서 CRLB에 얼마나 가까이 도달하는지 시각화하고 정량화한다.

실험 결과

연구 질문

• RQ1 낮은 SNR, 천체 기반 조건에서 PSF 스무딩 또는 가우시안 스무딩 이후 다항식 피팅을 수행하는 빠르고 반복하지 않는 중심점 추정 방법이 크래머-라오 하한(CRLB)에 도달할 수 있는가?
• RQ2 PSF가 알려져 있을 경우 PSF 매칭 다항식 중심점 추정의 성능은 최대우도 추정과 비교해 어떻게 되는가?
• RQ3 PSF가 알려져 있지 않을 경우 고정된 가우시안 스무딩 후 다항식 피팅의 성능은 CRLB에 비해 얼마나 되는가?
• RQ4 중심점-빛 방법은 CRLB에 얼마나 가까이 도달하지 못하는가? 고정밀 천체측위에 적합한가?
• RQ5 이 방법들의 RMSE 대 CRLB 비율은 신호 대 잡음비와 별의 PSF 반폭하강점(FWHM)에 따라 어떻게 변화하는가?

주요 결과

• PSF가 알려져 있을 경우 PSF 매칭 다항식 중심점 추정법은 다양한 SNR 및 FWHM 값 범위에서 약 1.0의 RMSE 대 CRLB 비율을 달성하여 CRLB에 거의 완전히 도달함을 나타낸다.
• PSF가 알려져 있지 않을 경우 고정된 가우시안 스무딩 후 다항식 피팅은 대부분의 테스트 조건에서 약 1.05~1.15의 RMSE 대 CRLB 비율을 기록하여 거의 최적의 성능을 보이며, 매우 효율적임을 보여준다.
• 중심점-빛 방법은 고SNR 조건에서도 일관되게 2.0 이상의 RMSE 대 CRLB 비율을 기록하여 근본적으로 최적화되지 않으며 고정밀 천체측위에 부적합함을 나타낸다.
• 모든 방법의 성능은 FWHM 증가와 SNR 감소에 따라 떨어지지만, CRLB에 대한 상대적 격차는 PSF 매칭 및 고정된 가우시안 다항식 방법에서 가장 작다.
• 진짜 PSF에 대한 최대우도 피팅은 가장 낮은 RMSE 대 CRLB 비율(약 1.0)을 기록하여 최적 기준으로서의 정당성을 확인한다.
• PSF 매칭 및 고정된 가우시안 다항식 방법의 RMSE 대 CRLB 비율은 시험 범위 내의 SNR(5~40) 및 FWHM(2~5.6 픽셀)에서 항상 10% 이내로 유지되어 강건성과 실용적 최적성을 입증한다.