Tag: PSF

Rethinking data-driven point spread function modeling with a differentiable optical model

Posted on by Jean-Luc Starck

 

Authors:

Tobias Liaudat, Jean-Luc Starck, Martin Kilbinger, Pierre-Antoine Frugier
Journal: Inverse Problems
Year: 2023
DOI:  
Download: ADS | arXiv

Abstract

In astronomy, upcoming space telescopes with wide-field optical instruments have a spatially varying point spread function (PSF). Specific scientific goals require a high-fidelity estimation of the PSF at target positions where no direct measurement of the PSF is provided. Even though observations of the PSF are available at some positions of the field of view (FOV), they are undersampled, noisy, and integrated into wavelength in the instrument's passband. PSF modeling represents a challenging ill-posed problem, as it requires building a model from these observations that can infer a super-resolved PSF at any wavelength and position in the FOV. Current data-driven PSF models can tackle spatial variations and super-resolution. However, they are not capable of capturing PSF chromatic variations. Our model, coined WaveDiff, proposes a paradigm shift in the data-driven modeling of the point spread function field of telescopes. We change the data-driven modeling space from the pixels to the wavefront by adding a differentiable optical forward model into the modeling framework. This change allows the transfer of a great deal of complexity from the instrumental response into the forward model. The proposed model relies on efficient automatic differentiation technology and modern stochastic first-order optimization techniques recently developed by the thriving machine-learning community. Our framework paves the way to building powerful, physically motivated models that do not require special calibration data. This paper demonstrates the WaveDiff model in a simplified setting of a space telescope. The proposed framework represents a performance breakthrough with respect to the existing state-of-the-art data-driven approach. The pixel reconstruction errors decrease six-fold at observation resolution and 44-fold for a 3x super-resolution. The ellipticity errors are reduced at least 20 times, and the size error is reduced more than 250 times. By only using noisy broad-band in-focus observations, we successfully capture the PSF chromatic variations due to diffraction. WaveDiff source code and examples associated with this paper are available at this link .

Deep Learning-based galaxy image deconvolution

Posted on by Jean-Luc Starck

 

Authors: Utsav Akhaury, Jean-Luc Starck, Pascale Jablonka, Frédéric Courbin, Kevin Michalewicz
Journal: A&A
Year: 2022
DOI:  
Download: ADS | arXiv

Abstract

With the onset of large-scale astronomical surveys capturing millions of images, there is an increasing need to develop fast and accurate deconvolution algorithms that generalize well to different images. A powerful and accessible deconvolution method would allow for the reconstruction of a cleaner estimation of the sky. The deconvolved images would be helpful to perform photometric measurements to help make progress in the fields of galaxy formation and evolution. We propose a new deconvolution method based on the Learnlet transform. Eventually, we investigate and compare the performance of different Unet architectures and Learnlet for image deconvolution in the astrophysical domain by following a two-step approach: a Tikhonov deconvolution with a closed-form solution, followed by post-processing with a neural network. To generate our training dataset, we extract HST cutouts from the CANDELS survey in the F606W filter (V-band) and corrupt these images to simulate their blurred-noisy versions. Our numerical results based on these simulations show a detailed comparison between the considered methods for different noise levels.

UNIONS: The impact of systematic errors on weak-lensing peak counts

UNIONS: The impact of systematic errors on weak-lensing peak counts

Posted on by Samuel Farrens

 

Authors: E. Ayçoberry, V. Ajani, A. Guinot, M. Kilbinger, V. Pettorino, S. Farrens, J.-L. Starck, R. Gavazzi, M. Hudson
Journal: A&A
Year: 2022
DOI:  
Download: ADS | arXiv

Abstract

Context. The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is an ongoing deep photometric multi-band survey of the Northern sky. As part of UNIONS, the Canada-France Imaging Survey (CFIS) provides r-band data which we use to study weak-lensing peak counts for cosmological inference.
Aims. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear bias, baryonic feedback, the source galaxy redshift estimate, intrinsic alignment, and the cluster member dilution.

Methods. For each uncertainty and systematic effect, we describe our mitigation scheme and the impact on cosmological parameter constraints. We obtain constraints on cosmological parameters from MCMC using CFIS data and MassiveNuS N-body simulations as a model for peak counts statistics.
Results. Depending on the calibration (local versus global, and the inclusion of the residual multiplicative shear bias), the mean matter density parameter Ωm can shift up to −0.024 (−0.5σ). We also see that including baryonic corrections can shift Ωm by +0.027 (+0.5σ) with respect to the DM-only simulations. Reducing the impact of the intrinsic alignment and cluster member dilution through signal-to-noise cuts can lead to a shift in Ωm of +0.027 (+0.5σ). Finally, with a mean redshift uncertainty of ∆z ̄ = 0.03, we see that the shift of Ωm (+0.001 which corresponds to +0.02σ) is not significant.

Conclusions. This paper investigates for the first time with UNIONS weak-lensing data and peak counts the impact of systematic effects. The value of Ωm is the most impacted and can shift up to ∼ 0.03 which corresponds to 0.5σ depending on the choices for each systematics. We expect constraints to become more reliable with future (larger) data catalogues, for which the current pipeline will provide a starting point. The code used to obtain the results is available in the following Github repository.

ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

ShapePipe: a new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

Posted on by Samuel Farrens

 

Authors: A. Guinot, M. Kilbinger, S. Farrens, A. Peel, A. Pujol, M. Schmitz, J.-L. Starck, T. Erben, R. Gavazzi, S. Gwyn, M. Hudson,  H. Hiledebrandt, T. Liaudat , et. al
Journal: A&A
Year: 2022
DOI:  
Download: ADS | arXiv

Abstract

UNIONS is an ongoing collaboration that will provide the largest deep photometric survey of the Northern sky in four optical bands to date. As part of this collaboration, CFIS is taking r-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. We perform the first weak-lensing analysis of CFIS r-band data over an area spanning 1700 deg2 of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies using the full UNIONS survey of 4800 deg2 when completed. Here we present ShapePipe, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as Ngmix for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the Point Spread Function (PSF), and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales <100 arcmin. The additive shear bias is below 5x10−4, and the residual multiplicative shear bias is at most 10−3 as measured on image simulations. Using COSEBIs we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than 4σ.

ShapePipe: A modular weak-lensing processing and analysis pipeline

ShapePipe: A modular weak-lensing processing and analysis pipeline

Posted on by Samuel Farrens

 

Authors: S. Farrens, A. Guinot, M. Kilbinger, T. Liaudat , L. Baumont, X. Jimenez, A. Peel , A. Pujol , M. Schmitz, J.-L. Starck, and A. Z. Vitorelli
Journal: A&A
Year: 2022
DOI: 10.1051/0004-6361/202243970
Download: ADS | arXiv

Abstract

We present the first public release of ShapePipe, an open-source and modular weak-lensing measurement, analysis, and validation pipeline written in Python. We describe the design of the software and justify the choices made. We provide a brief description of all the modules currently available and summarise how the pipeline has been applied to real Ultraviolet Near-Infrared Optical Northern Survey data. Finally, we mention plans for future applications and development. The code and accompanying documentation are publicly available on GitHub.

AMICO galaxy clusters in KiDS-DR3: measurement of the halo bias and power spectrum normalization from a stacked weak lensing analysis

AMICO galaxy clusters in KiDS-DR3: measurement of the halo bias and power spectrum normalization from a stacked weak lensing analysis

Posted on by Samuel Farrens

 

Authors: L. Ingoglia, G. Covone, M. Sereno, ..., S. Farrens, et al.
Journal: MNRAS
Year: 2022
DOI: 10.1093/mnras/stac046
Download: ADS | arXiv

Abstract

Galaxy clusters are biased tracers of the underlying matter density field. At very large radii beyond about 10 Mpc/\textit{h}, the shear profile shows evidence of a second-halo term. This is related to the correlated matter distribution around galaxy clusters and proportional to the so-called halo bias. We present an observational analysis of the halo bias-mass relation based on the AMICO galaxy cluster catalog, comprising around 7000 candidates detected in the third release of the KiDS survey. We split the cluster sample into 14 redshift-richness bins and derive the halo bias and the virial mass in each bin by means of a stacked weak lensing analysis. The observed halo bias-mass relation and the theoretical predictions based on the \\Lambda\CDM standard cosmological model show an agreement within \2\sigma\. The mean measurements of bias and mass over the full catalog give \M_{200c} = (4.9 \pm 0.3) \times 10^{13} M_{\odot}/\textit{h}\ and \b_h \sigma_8^2 = 1.2 \pm 0.1\. With the additional prior of a bias-mass relation from numerical simulations, we constrain the normalization of the power spectrum with a fixed matter density \\Omega_m = 0.3\, finding \\sigma_8 = 0.63 \pm 0.10\.

Deep transfer learning for blended source identification in galaxy survey data

Deep transfer learning for blended source identification in galaxy survey data

Posted on by Samuel Farrens

 

Authors: S. Farrens, A. Lacan, A. Guinot, A. Z. Vitorelli
Journal: A&A
Year: 2022
DOI: 10.1051/0004-6361/202141166
Download: ADS | arXiv

Abstract

We present BlendHunter, a proof-of-concept deep-transfer-learning-based approach for the automated and robust identification of blended sources in galaxy survey data. We take the VGG-16 network with pre-trained convolutional layers and train the fully connected layers on parametric models of COSMOS images. We test the efficacy of the transfer learning by taking the weights learned on the parametric models and using them to identify blends in more realistic Canada-France Imaging Survey (CFIS)-like images. We compare the performance of this method to SEP (a Python implementation of SExtractor) as a function of noise level and the separation between sources. We find that BlendHunter outperforms SEP by ~ 15% in terms of classification accuracy for close blends (< 10 pixel separation between sources) regardless of the noise level used for training. Additionally, the method provides consistent results to SEP for distant blends (>10 pixel separation between sources) provided the network is trained on data with noise that has a relatively close standard deviation to that of the target images. The code and data have been made publicly available to ensure the reproducibility of the results.

Semi-supervised dictionary learning with graph regularization and active points

Semi-supervised dictionary learning with graph regularization and active points

Posted on by Jean-Luc Starck

 

Authors: Khanh-Hung Tran,  Fred-Maurice Ngole-Mboula, J-L. Starck
Journal: SIAM Journal on Imaging Sciences
Year: 2020
DOI: 10.1137/19M1285469
Download: arXiv

Abstract

Supervised Dictionary Learning has gained much interest in the recent decade and has shown significant performance improvements in image classification. However, in general, supervised learning needs a large number of labelled samples per class to achieve an acceptable result. In order to deal with databases which have just a few labelled samples per class, semi-supervised learning, which also exploits unlabelled samples in training phase is used. Indeed, unlabelled samples can help to regularize the learning model, yielding an improvement of classification accuracy. In this paper, we propose a new semi-supervised dictionary learning method based on two pillars: on one hand, we enforce manifold structure preservation from the original data into sparse code space using Locally Linear Embedding, which can be considered a regularization of sparse code; on the other hand, we train a semi-supervised classifier in sparse code space. We show that our approach provides an improvement over state-of-the-art semi-supervised dictionary learning methods
.

Deep Learning for space-variant deconvolution in galaxy surveys

Deep Learning for space-variant deconvolution in galaxy surveys

Posted on by Jean-Luc Starck

 

Authors: Florent Sureau, Alexis Lechat, J-L. Starck
Journal: Astronomy and Astrophysics
Year: 2020
DOI: 10.1051/0004-6361/201937039
Download: ADS | arXiv

Abstract

The deconvolution of large survey images with millions of galaxies requires developing a new generation of methods that can take a space-variant point spread function into account. These methods have also to be accurate and fast. We investigate how deep learning might be used to perform this task. We employed a U-net deep neural network architecture to learn parameters that were adapted for galaxy image processing in a supervised setting and studied two deconvolution strategies. The first approach is a post-processing of a mere Tikhonov deconvolution with closed-form solution, and the second approach is an iterative deconvolution framework based on the alternating direction method of multipliers (ADMM). Our numerical results based on GREAT3 simulations with realistic galaxy images and point spread functions show that our two approaches outperform standard techniques that are based on convex optimization, whether assessed in galaxy image reconstruction or shape recovery. The approach based on a Tikhonov deconvolution leads to the most accurate results, except for ellipticity errors at high signal-to-noise ratio. The ADMM approach performs slightly better in this case. Considering that the Tikhonov approach is also more computation-time efficient in processing a large number of galaxies, we recommend this approach in this scenario.

In the spirit of reproducible research, the codes will be made freely available on the CosmoStat website (http://www.cosmostat.org). The testing datasets will also be provided to repeat the experiments performed in this paper.

PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing

PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing

Posted on by Samuel Farrens

 

Authors: S. Farrens, A. Grigis, L. El Gueddari, Z. Ramzi, Chaithya G. R., S. Starck, B. Sarthou, H. Cherkaoui, P.Ciuciu, J-L. Starck
Journal: Astronomy and Computing
Year: 2020
DOI: 10.1016/j.ascom.2020.100402
Download: ADS | arXiv

Abstract

We present the open-source image processing software package PySAP (Python Sparse data Analysis Package) developed for the COmpressed Sensing for Magnetic resonance Imaging and Cosmology (COSMIC) project. This package provides a set of flexible tools that can be applied to a variety of compressed sensing and image reconstruction problems in various research domains. In particular, PySAP offers fast wavelet transforms and a range of integrated optimisation algorithms. In this paper we present the features available in PySAP and provide practical demonstrations on astrophysical and magnetic resonance imaging data.

Code

PySAP Code