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.

We characterize the performance of the widely-used least-squares estimator in astrometry in terms of a comparison with the Cramer-Rao lower variance bound. In this inference context the performance of the least-squares estimator does not offer a closed-form expression, but a new result is presented (Theorem 1) where both the bias and the mean-square-error of the least-squares estimator are bounded and approximated analytically, in the latter case in terms of a nominal value and an interval around it. From the predicted nominal value we analyze how efficient is the least-squares estimator in comparison with the minimum variance Cramer-Rao bound. Based on our results, we show that, for the high signal-to-noise ratio regime, the performance of the least-squares estimator is significantly poorer than the Cramer-Rao bound, and we characterize this gap analytically. On the positive side, we show that for the challenging low signal-to-noise regime (attributed to either a weak astronomical signal or a noise-dominated condition) the least-squares estimator is near optimal, as its performance asymptotically approaches the Cramer-Rao bound. However, we also demonstrate that, in general, there is no unbiased estimator for the astrometric position that can precisely reach the Cramer-Rao bound. We validate our theoretical analysis through simulated digital-detector observations under typical observing conditions. We show that the nominal value for the mean-square-error of the least-squares estimator (obtained from our theorem) can be used as a benchmark indicator of the expected statistical performance of the least-squares method under a wide range of conditions. Our results are valid for an idealized linear (one-dimensional) array detector where intra-pixel response changes are neglected, and where flat-fielding is achieved with very high accuracy.

In this paper we use the Cramer-Rao lower uncertainty bound to estimate the maximum precision that could be achieved on the joint simultaneous (or 2D) estimation of photometry and astrometry of a point source measured by a linear CCD detector array. We develop exact expressions for the Fisher matrix elements required to compute the Cramer-Rao bound in the case of a source with a Gaussian light profile. From these expressions we predict the behavior of the Cramer-Rao astrometric and photometric precision as a function of the signal and the noise of the observations, and compare them to actual observations - finding a good correspondence between them. We show that the astrometric Cramer-Rao bound goes as $(S/N)^{-1}$ (similar to the photometric bound) but, additionally, we find that this bound is quite sensitive to the value of the background - suppressing the background can greatly enhance the astrometric accuracy. We present a systematic analysis of the elements of the Fisher matrix in the case when the detector adequately samples the source (oversampling regime), leading to closed-form analytical expressions for the Cramer-Rao bound. We show that, in this regime, the joint parametric determination of photometry and astrometry for the source become decoupled from each other, and furthermore, it is possible to write down expressions (approximate to first order in the small quantities F/B or B/F) for the expected minimum uncertainty in flux and position. These expressions are shown to be quite resilient to the oversampling condition, and become thus very valuable benchmark tools to estimate the approximate behavior of the maximum photometric and astrometric precision attainable under pre-specified observing conditions and detector properties.

In this paper we explore the maximum precision attainable in the location of a point source imaged by a pixel array detector in the presence of a background, as a function of the detector properties. For this we use a well-known result from parametric estimation theory, the so-called Cramer-Rao lower bound. We develop the expressions in the 1-dimensional case of a linear array detector in which the only unknown parameter is the source position. If the object is oversampled by the detector, analytical expressions can be obtained for the Cramer-Rao limit that can be readily used to estimate the limiting precision of an imaging system, and which are very useful for experimental (detector) design, observational planning, or performance estimation of data analysis software: In particular, we demonstrate that for background-dominated sources, the maximum astrometric precision goes as $B/F^2$, where $B$ is the background in one pixel, and $F$ is the total flux of the source, while when the background is negligible, this precision goes as $F^{-1}$. We also explore the dependency of the astrometric precision on: (1) the size of the source (as imaged by the detector), (2) the pixel detector size, and (3) the effect of source de-centering. Putting these results into context, the theoretical Cramer-Rao lower bound is compared to both ground- as well as spaced-based astrometric results, indicating that current techniques approach this limit very closely. Our results indicate that we have found in the Cramer-Rao lower variance bound a very powerful astrometric "benchmark" estimator concerning the maximum expected positional precision for a point source, given a prescription for the source, the background, the detector characteristics, and the detection process.