The high complexity of various inverse problems poses a significant challenge to model-based reconstruction schemes, which in such situations often reach their limits. At the same time, we witness an exceptional success of data-based methodologies such as deep learning. However, in the context of inverse problems, deep neural networks mostly act as black box routines, used for instance for a somewhat unspecified removal of artifacts in classical image reconstructions. In this paper, we will focus on the severely ill-posed inverse problem of limited angle computed tomography, in which entire boundary sections are not captured in the measurements. We will develop a hybrid reconstruction framework that fuses model-based sparse regularization with data-driven deep learning. Our method is reliable in the sense that we only learn the part that can provably not be handled by model-based methods, while applying the theoretically controllable sparse regularization technique to the remaining parts. Such a decomposition into visible and invisible segments is achieved by means of the shearlet transform that allows to resolve wavefront sets in the phase space. Furthermore, this split enables us to assign the clear task of inferring unknown shearlet coefficients to the neural network and thereby offering an interpretation of its performance in the context of limited angle computed tomography. Our numerical experiments show that our algorithm significantly surpasses both pure model- and more data-based reconstruction methods.

We propose a novel convolutional neural network (CNN), called $\Psi$DONet, designed for learning pseudodifferential operators ($\Psi$DOs) in the context of linear inverse problems. Our starting poi...

Tomographic reconstruction is an ill-posed inverse problem that calls for regularization. One possibility is to require sparsity of the unknown in an orthonormal wavelet basis. This, in turn, can be achieved by variational regularization, where the penalty term is the sum of the absolute values of the wavelet coefficients. The primal-dual fixed point algorithm showed that the minimizer of the variational regularization functional can be computed iteratively using a soft-thresholding operation. Choosing the soft-thresholding parameter μ > 0 is analogous to the notoriously difficult problem of picking the optimal regularization parameter in Tikhonov regularization. Here, a novel automatic method is introduced for choosing μ , based on a control algorithm driving the sparsity of the reconstruction to an a priori known ratio of nonzero versus zero wavelet coefficients in the unknown.

We consider a statistical inverse learning problem, where the task is to estimate a function $ f $ based on noisy point evaluations of $ Af $, where $ A $ is a linear operator. The function $ Af $ is evaluated at i.i.d. random design points $ u_n $, $ n = 1, ..., N $ generated by an unknown general probability distribution. We consider Tikhonov regularization with general convex and $ p $-homogeneous penalty functionals and derive concentration rates of the regularized solution to the ground truth measured in the symmetric Bregman distance induced by the penalty functional. We derive concrete rates for Besov norm penalties and numerically demonstrate the correspondence with the observed rates in the context of X-ray tomography.

In the context of international nuclear safeguards, the International Atomic Energy Agency (IAEA) has recently approved passive gamma emission tomography (PGET) as a method for inspecting spent nuclear fuel assemblies (SFAs). The PGET instrument is essentially a single photon emission computed tomography (SPECT) system that allows the reconstruction of axial cross-sections of the emission map of an SFA. The fuel material heavily self-attenuates its gamma-ray emissions, so that correctly accounting for the attenuation is a critical factor in producing accurate images. Due to the nature of the inspections, it is desirable to use as little a priori information as possible about the fuel, including the attenuation map, in the reconstruction process. Current reconstruction methods either do not correct for attenuation, assume a uniform attenuation throughout the fuel assembly, or assume an attenuation map based on an initial filtered back-projection reconstruction. We propose a method to simultaneously reconstruct the emission and attenuation maps by formulating the reconstruction as a constrained minimization problem with a least squares data fidelity term and regularization terms. Using simulated data, we show that our approach produces clear reconstructions which allow for a highly reliable classification of spent, missing, and fresh fuel rods.

In this paper we propose a motion-aware variational approach to reconstruct moving objects from sparse dynamic data. The motivation of this work stems from x-ray imaging of plants perfused with a liquid contrast agent, aimed at increasing the contrast of the images and studying the phloem transport in plants over time. The key idea of our approach is to deploy 3D shearlets as a space-temporal prior, treating time as the third dimension. The rationale behind this model is that a continuous evolution of a cartoon-like object is well suited for the use of 3D shearlets. We provide a basic mathematical analysis of the variational model for the image reconstruction. The numerical minimization is carried out with primal-dual scheme coupled with an automated choice of the regularization parameter. We test our model on different measurement setups: a simulated phantom especially designed to resemble a plant stem, with spreading points to simulate a spreading contrast agent; a measured agarose gel phantom to demonstrate iodide diffusion and geometry prior to imaging living sample; a measured living tree grown in vitro and perfused with a liquid sugar–iodine-mix. The results, compared against a 2D static model, show that our approach provides reconstructions that capture well the time dynamic of the contrast agent onset and are encouraging to develop microCT as a tool to study phloem transport using iodine tracer.