Optimal photometry of point sources: Joint source flux and background determination on array detectors -- from theory to practical implementation

Mario L. Vicuña,Jorge F. Silva,Rene A. Mendez,Marcos E. Orchard,Sebastian Espinosa,Jeremy Tregloan-Reed

2023-11-21

Abstract:In this paper we study the joint determination of source and background flux for point sources as observed by digital array detectors. We explicitly compute the two-dimensional Cramér-Rao absolute lower bound (CRLB) as well as the performance bounds for high-dimensional implicit estimators from a generalized Taylor expansion. This later approach allows us to obtain computable prescriptions for the bias and variance of the joint estimators. We compare these prescriptions with empirical results from numerical simulations in the case of the weighted least squares estimator (introducing an improved version, denoted stochastic weighted least-squares) as well as with the maximum likelihood estimator, finding excellent agreement. We demonstrate that these estimators provide quasi-unbiased joint estimations of the flux and background, with a variance that approaches the CRLB very tightly and are, hence, optimal, unlike the case of sequential estimation used commonly in astronomical photometry which is sub-optimal. We compare our predictions with numerical simulations of realistic observations, as well as with observations of a bona-fide non-variable stellar source observed with TESS, and compare it to the results from the sequential estimation of background and flux, confirming our theoretical expectations. Our practical estimators can be used as benchmarks for general photometric pipelines, or for applications that require maximum precision and accuracy in absolute photometry.

Instrumentation and Methods for Astrophysics,Data Analysis, Statistics and Probability

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to simultaneously optimize the estimation of source flux and background flux for point sources (such as stars) in astronomical photometry. Traditional methods usually adopt a sequential estimation approach, that is, estimating the background first and then the source flux. This approach has limitations when dealing with high - precision photometry. The author proposes a joint - estimation method. By considering the source flux and background flux simultaneously, the accuracy and precision of the estimation can be improved.

Specifically, the main contributions of the paper include:
1. **Theoretical derivation**: The two - dimensional Cramér - Rao absolute lower bound (CRLB) is clearly calculated, and the performance bounds of high - performance implicit estimators are derived from the perspective of the generalized Taylor expansion.
2. **Practical application**: An improved weighted least - squares method (stochastic weighted least - squares, SWLS) is introduced and compared with the maximum - likelihood estimation method. It is verified that these estimators provide almost unbiased joint estimations when close to the CRLB.
3. **Numerical simulation and empirical analysis**: The theoretical predictions are verified through numerical simulation and real observational data (such as TESS satellite data), and the superior performance of the joint - estimation method under different conditions is demonstrated.

### Formula presentation

#### 1. Point - source luminance distribution
The luminance distribution of a point source can be expressed by the following formula:
\[
\tilde{F}_{xc,\tilde{F}}(x)=\tilde{F}\cdot\phi(x - x_c,\sigma)
\]
where $\phi(x - x_c,\sigma)$ is a one - dimensional point - spread function (PSF), which describes the light distribution of the point source on the detector.

#### 2. Total expected flux per pixel
The total expected flux for each pixel is:
\[
\lambda_i(x_c,\tilde{F}) = E\{I_i\}=\tilde{F}\cdot g_i(x_c)+\tilde{B},\quad\forall i = 1,\ldots,n
\]
where $g_i(x_c)=\int_{x_i-\Delta x/2}^{x_i+\Delta x/2}\phi(x - x_c,\sigma)dx$ represents the contribution of the point source on the $i$ - th pixel.

#### 3. Cramér - Rao lower bound (CRLB)
For the joint - estimation problem, the Fisher information matrix is:
\[
I_\theta=\begin{bmatrix}
I_{1,1}&I_{1,2}\\
I_{2,1}&I_{2,2}
\end{bmatrix}
=\begin{bmatrix}
\sum_{i = 1}^n\frac{g_i(x_c)^2}{\lambda_i(x_c,\tilde{F})}&\sum_{i = 1}^n\frac{g_i(x_c)}{\lambda_i(x_c,\tilde{F})}\\
\sum_{i = 1}^n\frac{g_i(x_c)}{\lambda_i(x_c,\tilde{F})}&\sum_{i = 1}^n\frac{1}{\lambda_i(x_c,\tilde{F})}
\end{bmatrix}
\]

The corresponding CRLB is:
\[
\text{Var}(\hat{\tilde{F}})\geq\frac{I_{2,2}}{|I_\theta|},\quad\text{Var}(\hat{\tilde{B}})\geq\frac{I_{1,1}}{|I_\theta|}
\]

#### 4. WLS estimator
The cost function of the WLS estimator is:
\[
J_{WLS}(\alpha,I)=\sum_{i = 1}^n w_i(I_i-\alpha_1\cdot g_i(x_c)-\alpha_2)^2
\]