Superresolution Reconstruction of Severely Undersampled Point-spread Functions Using Point-source Stacking and Deconvolution

Abstract

Point-spread function (PSF) estimation in spatially undersampled images is challenging because large pixels average fine-scale spatial information. This is problematic when fine-resolution details are necessary, as in optimal photometry where knowledge of the illumination pattern beyond the native spatial resolution of the image may be required. Here, we introduce a method of PSF reconstruction where point sources are artificially sampled beyond the native resolution of an image and combined together via stacking to return a finely sampled estimate of the PSF. This estimate is then deconvolved from the pixel-gridding function to return a superresolution kernel that can be used for optimally weighted photometry. We benchmark against the <1% photometric error requirement of the upcoming SPHEREx mission to assess performance in a concrete example. We find that standard methods like Richardson-Lucy deconvolution are not sufficient to achieve this stringent requirement. We investigate a more advanced method with significant heritage in image analysis called iterative back-projection (IBP) and demonstrate it using idealized Gaussian cases and simulated SPHEREx images. In testing this method on real images recorded by the LORRI instrument on New Horizons, we are able to identify systematic pointing drift. Our IBP-derived PSF kernels allow photometric accuracy significantly better than the requirement in individual SPHEREx exposures. This PSF reconstruction method is broadly applicable to a variety of problems and combines computationally simple techniques in a way that is robust to complicating factors such as severe undersampling, spatially complex PSFs, noise, crowded fields, or limited source numbers.

Estimating the Feasibility of 21cm-Ly$α$ Synergies using the Hydrogen Epoch of Reionization Array

“…For these optimistic values, we take 𝜎 N = 3 × 10 −20 erg s −1 cm −2 sr −1 Hz −17 for SPHEREx, and 𝜎 N = 1.5 × 10 −21 erg s −1 cm −2 sr −1 Hz −1 for CDIM Heneka et al (2017). We do expect the nominal values to be the more realistic values measured, but better thermal noise values are possible given certain instrument assumptions (Doré et al 2018;Symons et al 2021). Eq.…”

Probing Intra-Halo Light with Galaxy Stacking in CIBER Images

“…However, given external source catalogs with high astrometric accuracy, we can stack on a sub-pixel basis and reconstruct the average source profile at scales finer than the native pixel size. This "sub-pixel stacking" technique has been used in previous CIBER imager analyses (Bock et al 2013;Zemcov et al 2014), and further investigated recently in the context of optimal photometry (Symons et al 2021). We summarize the sub-pixel stacking procedure as follows:…”

“…The pixel function is a matrix with each element proportional to the counts where the subpixel and the center sub-pixel that contains the source are within the same native pixel. The position of the center sub-pixel within the native pixel a uniform probability distribution, and therefore when stacking on a large number of sources, the pixel function converges to the analytic form (Symons et al 2021):…”

Measurements of the Zodiacal Light Absolute Intensity through Fraunhofer Absorption Line Spectroscopy with CIBER

“…The average PSF is measured in each flight independently by stacking on DSS star positions with 7 < M AB < 9. As the NBS design under-samples the PSF significantly, the stack is done on a sub-pixel grid, applying the technique implemented in Symons et al (2021).…”