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

arXiv: https://arxiv.org/abs/2203.04908 · HTML: https://mlsanity.com/paper/2203.04908/

## Key facts

- Cited by 4 papers in the map (data through 14 Oct 2026).
- Builds on 2 of its 49 references that are in the map.
- First posted on arXiv 9 Mar 2022.

## Builds on (most cited first, 2 of 2)

- [On the Variance of the Adaptive Learning Rate and Beyond](https://mlsanity.com/paper/1908.03265/) (2019)
- [Rethinking the modeling of the instrumental response of telescopes with a differentiable optical model](https://mlsanity.com/paper/2111.12541/) (2021)

## Built on by (most cited first, 4 of 4)

- [A Physics-Informed Blur Learning Framework for Imaging Systems](https://mlsanity.com/paper/2502.11382/) (2025)
- [Point spread function modelling for astronomical telescopes: a review focused on weak gravitational lensing studies](https://mlsanity.com/paper/2306.07996/) (2023)
- [On Uncertainty Calibration for Equivariant Functions](https://mlsanity.com/paper/2510.21691/) (2025)
- [An Image Quality Evaluation and Masking Algorithm Based On Pre-trained Deep Neural Networks](https://mlsanity.com/paper/2405.03408/) (2024)

## Abstract (arXiv)

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 degraded observations that can infer a super-resolved PSF at any wavelength and position in the FOV. 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 complexity from the instrumental response into the forward model. The proposed model relies on stochastic gradient descent to estimate its parameters. 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 6-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. Code available at https://github.com/tobias-liaudat/wf-psf.
