Image deconvolution and PSF reconstruction with STARRED: a wavelet-based two-channel method optimized for light-curve extraction

Martin Millon,Kevin Michalewicz,Frédéric Dux,Frédéric Courbin,Philip J. Marshall

DOI: https://doi.org/10.3847/1538-3881/ad4da7

2024-07-24

Abstract:We present STARRED, a Point Spread Function (PSF) reconstruction, two-channel deconvolution, and light curve extraction method designed for high-precision photometric measurements in imaging time series. An improved resolution of the data is targeted rather than an infinite one, thereby minimizing deconvolution artifacts. In addition, STARRED performs a joint deconvolution of all available data, accounting for epoch-to-epoch variations of the PSF and decomposing the resulting deconvolved image into a point source and an extended source channel. The output is a deep sharp frame combining all data, and the photometry of all point sources in the field of view as a function of time. Of note, STARRED also provides exquisite PSF models for each data frame. We showcase three applications of STARRED in the context of the imminent LSST survey and of JWST imaging: i) the extraction of supernovae light curves and the scene representation of their host galaxy, ii) the extraction of lensed quasar light curves for time-delay cosmography, and iii) the measurement of the spectral energy distribution of globular clusters in the "Sparkler", a galaxy at redshift z=1.378 strongly lensed by the galaxy cluster SMACS J0723.3-7327. STARRED is implemented in JAX, leveraging automatic differentiation and GPU acceleration. This enables rapid processing of large time-domain datasets, positioning the method as a powerful tool for extracting light curves from the multitude of lensed or unlensed variable and transient objects in the Rubin-LSST data, even when blended with intervening objects.

Instrumentation and Methods for Astrophysics,Astrophysics of Galaxies

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to improve the resolution and reduce deconvolution artifacts through high - precision photometry in astronomical image time series. Specifically, the authors propose a new method named STARRED for point - spread - function (PSF) reconstruction, two - channel deconvolution, and light - curve extraction. The following is a specific description of this problem:

1. **Image Blurring Caused by Atmospheric Turbulence and Optical Diffraction**:
The resolution of astronomical images is affected by the Earth's atmospheric turbulence or the optical diffraction of space telescopes, resulting in image blurring. This blurring can be corrected by understanding the point - spread - function (PSF).

2. **The Impact of the Temporal Variation of PSF on Photometry**:
The PSF changes over time, especially in different observation periods, which will lead to inaccurate photometry of mixed light sources. This systematic error cannot be reduced by increasing the exposure time; on the contrary, it will worsen as the signal - to - noise ratio increases.

3. **The Mathematical Difficulties of Deconvolution**:
The deconvolution problem is mathematically ill - posed because of the influence of pixelation and noise, resulting in no unique and stable solution. To solve this problem, a regularization term is usually introduced to constrain the solution space.

4. **The Artifact Problem Caused by Complete Deconvolution**:
Completely removing the PSF effect will cause artifacts (such as the "ringing" phenomenon) to appear in the image because the deconvolved point source cannot be accurately represented by any sampling or pixel size. This violates the Nyquist - Shannon sampling theorem.

5. **The Need for Multi - Period Data Processing**:
A method that can process multi - period data is required to jointly deconvolve all available data and consider the PSF changes in each observation period.

To meet these challenges, the STARRED method adopts the following strategies:

- **Partial Deconvolution**: Instead of attempting to completely remove the PSF effect, selectively improve the resolution to ensure that the deconvolved image satisfies the Nyquist - Shannon sampling theorem.

- **Two - Channel Decomposition**: Decompose the image into a point - source channel and an extended - source channel, and process point - sources and extended - sources separately.

- **Wavelet Regularization**: Use starlets for regularization to sparsely represent astronomical objects, thereby reducing the influence of noise.

- **Automatic Differentiation and GPU Acceleration**: Utilize the JAX library to achieve automatic differentiation and GPU acceleration for rapid processing of large - scale time - domain data sets.

Through these improvements, STARRED can provide high - resolution images and accurate light - curve extraction while maintaining a high signal - to - noise ratio, and is suitable for various astronomical applications, such as supernova light - curve extraction, time - delay measurement of lensed quasars, and spectral energy distribution measurement of distant star clusters.