Optimality of the Maximum Likelihood estimator in Astrometry

Abstract

The problem of astrometry is revisited from the perspective of analyzing the attainability of well-known performance limits (the Cramer-Rao bound) for the estimation of the relative position of light-emitting (usually point-like) sources on a CCD-like detector using commonly adopted estimators such as the weighted least squares and the maximum likelihood. Novel technical results are presented to determine the performance of an estimator that corresponds to the solution of an optimization problem in the context of astrometry. Using these results we are able to place stringent bounds on the bias and the variance of the estimators in close form as a function of the data. We confirm these results through comparisons to numerical simulations under a broad range of realistic observing conditions. The maximum likelihood and the weighted least square estimators are analyzed. We confirm the sub-optimality of the weighted least squares scheme from medium to high signal-to-noise found in an earlier study for the (unweighted) least squares method. We find that the maximum likelihood estimator achieves optimal performance limits across a wide range of relevant observational conditions. Furthermore, from our results, we provide concrete insights for adopting an adaptive weighted least square estimator that can be regarded as a computationally efficient alternative to the optimal maximum likelihood solution. We provide, for the first time, close-form analytical expressions that bound the bias and the variance of the weighted least square and maximum likelihood implicit estimators for astrometry using a Poisson-driven detector. These expressions can be used to formally assess the precision attainable by these estimators in comparison with the minimum variance bound.

arXiv comment24 pages, 7 figures, 2 tables, 3 appendices. Accepted by Astronomy & Astrophysics

References in corpus (15)

• The Gaia mission
• A method for optimal image subtraction
• Astrometry with Hubble Space Telescope Fine Guidance Sensors - A Review
• Gravitationally Lensed Quasars in Gaia: I. Resolving Small-Separation Lenses
• Analysis and interpretation of the Cramer-Rao lower-bound in astrometry: One dimensional case
• Gaia astrometry for stars with too few observations - a Bayesian approach
• Gaia
• Proper motions of Local Group dwarf spheroidal galaxies I: First ground-based results for Fornax
• GBOT: ground based optical tracking of the Gaia satellite
• Quasars can be used to verify the parallax zero-point of the Tycho-Gaia Astrometric Solution
• Performance of an Algorithm for Estimation of Flux, Background and Location on One-Dimensional Signals
• Performance analysis of the Least-Squares estimator in Astrometry
• Analysis of the Bayesian Cramer-Rao lower bound in astrometry: Studying the impact of prior information in the location of an object
• Characterizing the astrometric precision limit for moving targets observed with digital-array detectors
• General constraints on influential error sources for super-high accuracy star tracker

Cited by in corpus (3)

• C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: An analysis of the precision achievable when modeling spectral line within the Local Thermodynamic Equilibrium approximation
• Investigation and Application of Fitting Models for Centering Algorithms in Astrometry
• A review on fundamental bounds and estimators for photometry and astrometry of celestial point sources using array detectors, from first principles