Multi-CCD Point Spread Function Modelling
arXiv:2011.09835 · doi:10.1051/0004-6361/202039584
Abstract
Galaxy imaging surveys observe a vast number of objects that are affected by the instrument's Point Spread Function (PSF). Weak lensing missions, in particular, aim at measuring the shape of galaxies, and PSF effects represent an important source of systematic errors which must be handled appropriately. This demands a high accuracy in the modelling as well as the estimation of the PSF at galaxy positions. Sometimes referred to as non-parametric PSF estimation, the goal of this paper is to estimate a PSF at galaxy positions, starting from a set of noisy star image observations distributed over the focal plane. To accomplish this, we need our model to first of all, precisely capture the PSF field variations over the Field of View (FoV), and then to recover the PSF at the selected positions. This paper proposes a new method, coined MCCD (Multi-CCD PSF modelling), that creates, simultaneously, a PSF field model over all of the instrument's focal plane. This allows to capture global as well as local PSF features through the use of two complementary models which enforce different spatial constraints. Most existing non-parametric models build one model per Charge Coupled Device (CCD), which can lead to difficulties in capturing global ellipticity patterns. We first test our method on a realistic simulated dataset comparing it with two state-of-the-art PSF modelling methods (PSFEx and RCA). We outperform both of them with our proposed method. Then we contrast our approach with PSFEx on real data from CFIS (Canada France Imaging Survey) that uses the CFHT (Canada-France-Hawaii Telescope). We show that our PSF model is less noisy and achieves a 22% gain on pixel Root Mean Squared Error (RMSE) with respect to PSFEx. We present, and share the code of, a new PSF modelling algorithm that models the PSF field on all the focal plane that is mature enough to handle real data.
19 pages, 13 figures. Submitted
References in corpus (29)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Astropy: A Community Python Package for Astronomy
- The NumPy array: a structure for efficient numerical computation
- Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
- Cosmology with cosmic shear observations: a review
- Bayesian Galaxy Shape Measurement for Weak Lensing Surveys - III. Application to the Canada-France-Hawaii Telescope Lensing Survey
- Gravitational Lensing Analysis of the Kilo Degree Survey
- CFHTLenS: The Canada-France-Hawaii Telescope Lensing Survey - Imaging Data and Catalogue Products
- Shear calibration biases in weak lensing surveys
- Evidence for the accelerated expansion of the Universe from weak lensing tomography with COSMOS
- Systematic errors in weak lensing: application to SDSS galaxy-galaxy weak lensing
- Dark Energy Survey Year 1 Results: Weak Lensing Shape Catalogues
- The DES Science Verification Weak Lensing Shear Catalogues
- Principal Component Analysis of the Time- and Position-Dependent Point Spread Function of the Advanced Camera for Surveys
- Defining a weak lensing experiment in space
- The Canada-France Imaging Survey: First results from the u-band component
- Wasserstein Dictionary Learning: Optimal Transport-based unsupervised non-linear dictionary learning
- Dark Energy Survey Year 3 Results: Point-Spread Function Modeling
- Interpolating point spread function anisotropy
- The Sloan Digital Sky Survey extended point spread functions
- Fast Point Spread Function Modeling with Deep Learning
- PSF--NET: A Non-parametric Point Spread Function Model for Ground Based Optical Telescopes
- PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing
- PSF modelling for very wide-field CCD astronomy
- Super-resolution method using sparse regularization for point-spread function recovery
- Point Spread Function Modelling for Wide Field Small Aperture Telescopes with a Denoising Autoencoder
- Euclid: Nonparametric point spread function field recovery through interpolation on a graph Laplacian
- Pushing Point spread function reconstruction to the next level. Application to SPHERE/ZIMPOL
- Constraint matrix factorization for space variant PSFs field restoration
Cited by in corpus (8)
- ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data
- Point spread function modelling for astronomical telescopes: a review focused on weak gravitational lensing studies
- Rethinking data-driven point spread function modeling with a differentiable optical model
- ShapePipe: A modular weak-lensing processing and analysis pipeline
- UNIONS: a direct measurement of intrinsic alignment with BOSS/eBOSS spectroscopy
- Galaxy-Point Spread Function correlations as a probe of weak-lensing systematics with UNIONS data
- Rethinking the modeling of the instrumental response of telescopes with a differentiable optical model
- Breaking the degeneracy in stellar spectral classification from single wide-band images