Neural blind deconvolution to reconstruct high-resolution ground-based solar observations

TL;DR

This paper introduces a physics-informed neural network approach for high-resolution solar image reconstruction from short exposure bursts, overcoming atmospheric limitations and surpassing existing methods in resolving small-scale solar features.

Contribution

The paper presents a novel neural network method that jointly estimates the PSF and high-resolution images, improving solar observation reconstructions beyond current techniques.

Findings

Successfully reconstructs small-scale solar features

Outperforms state-of-the-art reconstruction methods

Demonstrated on synthetic and real data from GREGOR and DKIST

Abstract

Ground-based solar observations enable unprecedented spatial, spectral, and temporal resolution of the lower solar atmosphere, yet Earths turbulent atmosphere imposes significant limitations, requiring advanced post-facto image reconstruction. State-of-the-art reconstruction methods are based on restoring a burst of short exposure frames to a single observation. Limitations of these techniques arise due to the sparse information about the atmospheric point spread function (PSF) that degrade the observations and consequently the quality of reconstructions. We develop a novel image reconstruction method to achieve unprecedented spatial resolution from short exposure image bursts. This can provide high-quality reconstructions and therefore advance the study of the smallest spatial scales from the solar photosphere to the chromosphere. In this study, we present a novel approach for high-resolution solar image reconstruction based on physics-informed neural networks. In the training process, the neural network maps coordinate points directly to their corresponding intensity values while simultaneously updating the PSF parameters. The method convolves the true object from the neural network with the estimated PSFs and optimizes the network by minimizing the loss between the synthesized and real short-exposure image burst. This approach enables the simultaneous estimation of both the degrading PSF and the real high-resolution intensity distribution. We demonstrate the method on synthetic intensity data derived from a radiative MHD simulation and apply it to high-resolution observations from GREGOR and DKIST. Our results demonstrate the ability to reconstruct small-scale solar features that exceed the reconstruction performance of state-of-the-art reconstruction methods. With this approach we lay the foundation for future spatially varying PSFs.