Accelerating Multiframe Blind Deconvolution via Deep Learning

Ground-based solar image restoration is a computationally expensive procedure that involves nonlinear optimization techniques. The presence of atmospheric turbulence produces perturbations in individual images that make it necessary to apply blind deconvolution techniques. These techniques rely on the observation of many short exposure frames that are used to simultaneously infer the instantaneous state of the atmosphere and the unperturbed object. We have recently explored the use of machine learning to accelerate this process, with promising results. We build upon this previous work to propose several interesting improvements that lead to better models. As well, we propose a new method to accelerate the restoration based on algorithm unrolling. In this method, the image restoration problem is solved with a gradient descent method that is unrolled and accelerated aided by a few small neural networks. The role of the neural networks is to correct the estimation of the solution at each iterative step. The model is trained to perform the optimization in a small fixed number of steps with a curated dataset. Our findings demonstrate that both methods significantly reduce the restoration time compared to the standard optimization procedure. Furthermore, we showcase that these models can be trained in an unsupervised manner using observed images from three different instruments. Remarkably, they also exhibit robust generalization capabilities when applied to new datasets. To foster further research and collaboration, we openly provide the trained models, along with the corresponding training and evaluation code, as well as the training dataset, to the scientific community.Comment: 26 pages, 9 figures, accepted for publication in Solar Physic

Stokes imaging polarimetry using image restoration at the Swedish 1-m Solar Telescope

Aims: We aim to achieve high spatial resolution as well as high polarimetric sensitivity, using an earth-based 1m-class solar telescope, for the study of magnetic fine structure on the Sun. Methods: We use a setup with 3 high-speed, low-noise cameras to construct datasets with interleaved polarimetric states, particularly suitable for Multi-Object Multi-Frame Blind Deconvolution image restorations. We discuss the polarimetric calibration routine as well as various potential sources of error in the results. Results: We obtained near diffraction limited images, with a noise level of approximately 10^(-3) I(cont). We confirm that dark-cores have a weaker magnetic field and at a lower inclination angle with respect to the solar surface than the edges of the penumbral filament. We show that the magnetic field strength in faculae-striations is significantly lower than in other nearby parts of the faculae.Comment: Accepted for publication in Astronomy & Astrophysics, 12 pages, 11 figure

High-order aberration compensation with Multi-frame Blind Deconvolution and Phase Diversity image restoration techniques

Context. For accurately measuring intensities and determining magnetic field strengths of small-scale solar (magnetic) structure, knowledge of and compensation for the point spread function is crucial. For images recorded with the Swedish 1-meter Solar Telescope, restoration with Multi-Frame Blind Deconvolution and Joint Phase Diverse Speckle methods lead to remarkable improvements in image quality but granulation contrasts that are too low, indicating additional stray light. Aims. We propose a method to compensate for stray light from high-order atmospheric aberrations not included in MFBD and JPDS processing. Methods. To compensate for uncorrected aberrations, a reformulation of the image restoration process is proposed that allows the average effect of hundreds of high-order modes to be compensated for by relying on Kolmogorov statistics for these modes. The applicability of the method requires simultaneous measurements of Fried's parameter r0. The method is tested with simulations as well as real data and extended to include compensation for conventional stray light. Results. We find that only part of the reduction of granulation contrast in SST images is due to uncompensated high-order aberrations. The remainder is still unaccounted for and attributed to stray light from the atmosphere, the telescope with its re-imaging system and to various high-altitude seeing effects. Conclusions. We conclude that statistical compensation of high-order modes is a viable method to reduce the loss of contrast occurring when a limited number of aberrations is explicitly compensated for with MFBD and JPDS processing. We show that good such compensation is possible with only 10 recorded frames. The main limitation of the method is that already MFBD and JPDS processing introduces high-order compensation that, if not taken into account, can lead to over-compensation.Comment: in press in Astronomy & Astrophysic

Improved SOT (Hinode mission) high resolution solar imaging observations

We consider the best today available observations of the Sun free of turbulent Earth atmospheric effects, taken with the Solar Optical Telescope (SOT) onboard the Hinode spacecraft. Both the instrumental smearing and the observed stray light are analyzed in order to improve the resolution. The Point Spread Function (PSF) corresponding to the blue continuum Broadband Filter Imager (BFI) near 450 nm is deduced by analyzing i/ the limb of the Sun and ii/ images taken during the transit of the planet Venus in 2012. A combination of Gaussian and Lorentzian functions is selected to construct a PSF in order to remove both smearing due to the instrumental diffraction effects (PSF core) and the large-angle stray light due to the spiders and central obscuration (wings of the PSF) that are responsible for the parasitic stray light. A Max-likelihood deconvolution procedure based on an optimum number of iterations is discussed. It is applied to several solar field images, including the granulation near the limb. The normal non-magnetic granulation is compared to the abnormal granulation which we call magnetic. A new feature appearing for the first time at the extreme- limb of the disk (the last 100 km) is discussed in the context of the definition of the solar edge and of the solar diameter. A single sunspot is considered in order to illustrate how effectively the restoration works on the sunspot core. A set of 125 consecutive deconvolved images is assembled in a 45 min long movie illustrating the complexity of the dynamical behavior inside and around the sunspot.Comment: 15 pages, 22 figures, 1 movi

Project Tech Top study of lunar, planetary and solar topography Final report

Data acquisition techniques for information on lunar, planetary, and solar topograph

Two-Dimensional Spectroscopy of Photospheric Shear Flows in a Small delta Spot

In recent high-resolution observations of complex active regions, long-lasting and well-defined regions of strong flows were identified in major flares and associated with bright kernels of visible, near-infrared, and X-ray radiation. These flows, which occurred in the proximity of the magnetic neutral line, significantly contributed to the generation of magnetic shear. Signatures of these shear flows are strongly curved penumbral filaments, which are almost tangential to sunspot umbrae rather than exhibiting the typical radial filamentary structure. Solar active region NOAA 10756 was a moderately complex, beta-delta sunspot group, which provided an opportunity to extend previous studies of such shear flows to quieter settings. We conclude that shear flows are a common phenomenon in complex active regions and delta spots. However, they are not necessarily a prerequisite condition for flaring. Indeed, in the present observations, the photospheric shear flows along the magnetic neutral line are not related to any change of the local magnetic shear. We present high-resolution observations of NOAA 10756 obtained with the 65-cm vacuum reflector at Big Bear Solar Observatory (BBSO). Time series of speckle-reconstructed white-light images and two-dimensional spectroscopic data were combined to study the temporal evolution of the three-dimensional vector flow field in the beta-delta sunspot group. An hour-long data set of consistent high quality was obtained, which had a cadence of better than 30 seconds and sub-arcsecond spatial resolution.Comment: 23 pages, 6 gray-scale figures, 4 color figures, 2 tables, submitted to Solar Physic

Modern optical astronomy: technology and impact of interferometry

The present `state of the art' and the path to future progress in high spatial resolution imaging interferometry is reviewed. The review begins with a treatment of the fundamentals of stellar optical interferometry, the origin, properties, optical effects of turbulence in the Earth's atmosphere, the passive methods that are applied on a single telescope to overcome atmospheric image degradation such as speckle interferometry, and various other techniques. These topics include differential speckle interferometry, speckle spectroscopy and polarimetry, phase diversity, wavefront shearing interferometry, phase-closure methods, dark speckle imaging, as well as the limitations imposed by the detectors on the performance of speckle imaging. A brief account is given of the technological innovation of adaptive-optics (AO) to compensate such atmospheric effects on the image in real time. A major advancement involves the transition from single-aperture to the dilute-aperture interferometry using multiple telescopes. Therefore, the review deals with recent developments involving ground-based, and space-based optical arrays. Emphasis is placed on the problems specific to delay-lines, beam recombination, polarization, dispersion, fringe-tracking, bootstrapping, coherencing and cophasing, and recovery of the visibility functions. The role of AO in enhancing visibilities is also discussed. The applications of interferometry, such as imaging, astrometry, and nulling are described. The mathematical intricacies of the various `post-detection' image-processing techniques are examined critically. The review concludes with a discussion of the astrophysical importance and the perspectives of interferometry.Comment: 65 pages LaTeX file including 23 figures. Reviews of Modern Physics, 2002, to appear in April issu