Abstract. Diffuse optical tomography (DOT) is an imaging modality in which images of the optical properties of a biological tissue, specifically diffusion [Formula: see text] and absorption [Formula: see text], are estimated based on measurements of near-infrared light on the surface of the body. In practical applications, one often lacks exact knowledge of the measurement domain boundary. This poses a significant challenge, as inaccuracies in the boundary shape of the computational domain may result in substantial artifacts in the reconstructed images. In this study, the following two results are achieved in a two-dimensional setting: (i) when the measurement domain [Formula: see text] is known, it is shown that knowledge of the Robin-to-Neumann map for two modulation frequencies uniquely determines [Formula: see text] and [Formula: see text], both assumed to be isotropic. (ii) When the measurement domain [Formula: see text] is not exactly known, a method is proposed for simultaneous reconstruction of [Formula: see text] and [Formula: see text] as well as the boundary [Formula: see text]. For this, DOT measurements are needed for two modulation frequencies. This approach yields reconstructed coefficients that are a conformal deformation of the true coefficients in the exact domain. The new method is demonstrated using simulated noisy data.

Abstract. In this paper, we propose a new data assimilation approach based on iterative filtering and smoothing in the expectation-maximization fashion, incorporating optical flow prediction models for dynamic state estimation. The concept is suitable for applications where state estimation relies on two-dimensional images and where no natural physical prediction model is available. We apply the proposed approach to dynamic X-ray images and both real and synthetic satellite data, demonstrating that iterative filtering and smoothing with optical flow models yields improved results compared to using an identity model approach. The quality of the state estimates improves both visually and in terms of the root mean squared error, typically after a couple of iterations. However, in some cases, continued iterations may lead to deteriorating results. Therefore, monitoring the quality of both optical flow and state estimates is crucial in the iterative approach.

Abstract. A novel reconstruction method is introduced for the severely ill-posed inverse problem of limited-angle tomography. It is well known that, depending on the available measurement, angles specify a subset of the wavefront set of the unknown target, while some oriented singularities remain invisible in the data. Topological Interface recovery for Limited-angle Tomography, or TILT, is based on lifting the visible part of the wavefront set under a universal covering map. In the space provided, it is possible to connect the appropriate pieces of the lifted wavefront set correctly using dual-tree complex wavelets, a dedicated metric, and persistent homology. The result is not only a suggested invisible boundary but also a computational representation for all interfaces in the target.

The objective of electrical impedance tomography (EIT) is to reconstruct the internal conductivity of a physical body based on current and voltage measurements at the boundary of the body. In many medical applications the exact shape of the domain boundary and contact impedances are not available. This is problematic as even small errors in the boundary shape of the computation domain or in the contact impedance values can produce large artifacts in the reconstructed images, which results in a loss of relevant information. A method is proposed that simultaneously reconstructs the conductivity, the contact impedances, and the boundary shape from EIT data. The approach consists of three steps: first, the unknown contact impedances and an anisotropic conductivity reproducing the measured EIT data in a model domain are computed. Second, using isothermal coordinates, a deformation is constructed that makes the conductivity isotropic. The final step minimizes the error of true and reconstructed known geometric properties (like the electrode lengths) using conformal deformations. The feasibility of the method is illustrated with experimental EIT data, with robust and accurate reconstructions of both conductivity and boundary shape.

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 point is the iterative soft thresholding algorithm (ISTA), a well-known algorithm to solve sparsity-promoting minimization problems. We show that, under rather general assumptions on the forward operator, the unfolded iterations of ISTA can be interpreted as the successive layers of a CNN, which in turn provides fairly general network architectures that, for a specific choice of the parameters involved, allow us to reproduce ISTA, or a perturbation of ISTA for which we can bound the coefficients of the filters. Our case study is the limited-angle X-ray transform and its application to limited-angle computed tomography (LA-CT). In particular, we prove that, in the case of LA-CT, the operations of upscaling, downscaling, and convolution, which characterize our $\Psi$DONet and most deep learning schemes, can be exactly determined by combining the convolutional nature of the limited-angle X-ray transform and basic properties defining an orthogonal wavelet system. We test two different implementations of $\Psi$DONet on simulated data from limited-angle geometry, generated from the ellipse data set. Both implementations provide equally good and noteworthy preliminary results, showing the potential of the approach we propose and paving the way to applying the same idea to other convolutional operators which are $\Psi$DOs or Fourier integral operators.

X-ray tomography is a reliable tool for determining the inner structure of 3-D object with penetrating X-rays. However, traditional reconstruction methods, such as Feldkamp-Davis-Kress (FDK), require dense angular sampling in the data acquisition phase leading to long measurement times, especially in X-ray micro-tomography to obtain high-resolution scans. Acquiring less data using greater angular steps is an obvious way for speeding up the process and avoiding the need to save huge data sets. However, computing 3-D reconstruction from such a sparsely sampled data set is difficult because the measurement data are usually contaminated by errors, and linear measurement models do not contain sufficient information to solve the problem in practice. An automatic regularization method is proposed for robust reconstruction, based on enforcing sparsity in the 3-D shearlet-transform domain. The inputs of the algorithm are the projection data and a priori known expected degree of sparsity, denoted as 0pr≤ 1. The number Cpr can be calibrated from a few dense-angle reconstructions and fixed. Human subchondral bone samples were tested, and morphometric parameters of the bone reconstructions were then analyzed using standard metrics. The proposed method is shown to outperform the baseline algorithm (FDK) in the case of sparsely collected data. The number of X-ray projections can be reduced up to 10% of the total amount 300 projections over 180° with uniform angular step while retaining the quality of the reconstruction images and of the morphometric parameters.

The regularized D-bar method for electrical impedance tomography (EIT) provides a rigorous mathematical approach for solving the full nonlinear inverse problem directly, i.e., without iterations. It is based on a low-pass filtering in the (nonlinear) frequency domain. However, the resulting D-bar reconstructions are inherently smoothed, leading to a loss of edge distinction. In this paper, a novel method that combines a D-bar approach with the edge-preserving nature of total variation (TV) regularization is presented. The method also includes a data-driven contrast adjustment technique guided by the key functions (CGO solutions) of the D-bar method. The new TV-enhanced D-bar method produces reconstructions with sharper edges and improved contrast. This is achieved by using the TV-induced edges to increase the truncation radius of the scattering data in the nonlinear frequency domain, thereby increasing the radius of the low-pass filter. The algorithm is tested on numerically simulated noisy EIT data and demonstrates significant improvements in edge preservation and contrast which can be highly valuable for absolute EIT imaging.

A computational method is introduced for choosing the regularization parameter for total variation (TV) regularization. A partial understanding of the properties of the method is provided by rigorously proving that the TV norms of the reconstructions converge with any choice of regularization parameter. The computational approach is based on computing reconstructions at a few different resolutions and various values of regularization parameter. The chosen parameter is the smallest one resulting in approximately discretization-invariant TV norms of the reconstructions. The method is tested with simulated and experimental X-ray tomography data and compared to the S-curve method. The results are comparable to those of the S-curve method. However, the S-curve method needs quantitative a priori information about the expected sparsity (TV norm) of the unknown, while the proposed method does not require such input parameters.

This paper presents a new glottal inverse filtering (GIF) method that utilizes Markov chain Monte Carlo (MCMC) algorithm. First, initial estimates of the vocal tract and glottal flow are eval-uated by an existing GIF method, the iterative adaptive inverse filtering (IAIF). Simultaneously, the initially estimated glottal flow is synthesized using the Klatt model and filtered with the estimated vocal tract filter. In the MCMC estimation process, the first few poles of the initial vocal tract model and the Klatt parameter are refined in order to minimize the error between the original and the synthetic signals. MCMC converges to the optimal result, and the final estimate of the vocal tract is found by averaging the parameter values of the Markov chain. Ex-periments show that the MCMC-based GIF method gives more accurate results compared to the original IAIF method.

Diagnostic and operational tasks based on dental radiology often require three-dimensional (3-D) information that is not available in a single X-ray projection image. Comprehensive 3-D information about tissues can be obtained by computerized tomography (CT) imaging. However, in dental imaging a conventional CT scan may not be available or practical because of high radiation dose, low-resolution or the cost of the CT scanner equipment. In this paper, we consider a novel type of 3-D imaging modality for dental radiology. We consider situations in which projection images of the teeth are taken from a few sparsely distributed projection directions using the dentist's regular (digital) X-ray equipment and the 3-D X-ray attenuation function is reconstructed. A complication in these experiments is that the reconstruction of the 3-D structure based on a few projection images becomes an ill-posed inverse problem. Bayesian inversion is a well suited framework for reconstruction from such incomplete data. In Bayesian inversion, the ill-posed reconstruction problem is formulated in a well-posed probabilistic form in which a priori information is used to compensate for the incomplete information of the projection data. In this paper we propose a Bayesian method for 3-D reconstruction in dental radiology. The method is partially based on Kolehmainen et al. 2003. The prior model for dental structures consist of a weighted l1 and total variation (TV)-prior together with the positivity prior. The inverse problem is stated as finding the maximum a posteriori (MAP) estimate. To make the 3-D reconstruction computationally feasible, a parallelized version of an optimization algorithm is implemented for a Beowulf cluster computer. The method is tested with projection data from dental specimens and patient data. Tomosynthetic reconstructions are given as reference for the proposed method.

The aim of X-ray tomography is to reconstruct an unknown physical body from a collection of projection images. When the projection images are only available from a limited angle of view, the reconstruction problem is a severely ill-posed inverse problem. Statistical inversion allows stable solution of the limited-angle tomography problem by complementing the measurement data by a priori information. In this work, the unknown attenuation distribution inside the body is represented as a wavelet expansion, and a Besov space prior distribution together with positivity constraint is used. The wavelet expansion is thresholded before reconstruction to reduce the dimension of the computational problem. Feasibility of the method is demonstrated by numerical examples using in vitro data from mammography and dental radiology.

The problem this paper addresses is how to use the two-dimensional D-bar method for electrical impedance tomography with experimental data collected on finitely many electrodes covering a portion of the boundary of a body. This requires an approximation of the Dirichlet-to-Neumann, or voltage-to-current density map, defined on the entire boundary of the region, from a finite number of matrix elements of the current-to-voltage map. Reconstructions from experimental data collected on a saline filled tank containing agar heart and lung phantoms are presented, and the results are compared to reconstructions by the NOSER algorithm on the same data.

We study the problem of gathering information about wind velocity from continuous-wave clear-air Doppler radar measurements. The radar is assumed to be a monostatic fixed-frequency Doppler radar, and the wind velocity as well as the reflectivity are modeled as random fields with statistical parameters depending on the altitude. We seek to reconstruct the hodograph curve of the wind profile, i.e., the projection of the wind velocity profile to the ground plane. We show that under certain assumptions of the wind field, the problem reduces to a well-known problem occurring in classical X-ray tomography. Numerical simulations based on the use of X-ray inversion methods are presented.

A direct (noniterative) reconstruction algorithm for electrical impedance tomography in the two-dimensional (2-D), cross-sectional geometry is reviewed. New results of a reconstruction of a numerically simulated phantom chest are presented. The algorithm is based on the mathematical uniqueness proof by A. I. Nachman [1996] for the 2-D inverse conductivity problem. In this geometry, several of the clinical applications include monitoring heart and lung function, diagnosis of pulmonary embolus, diagnosis of pulmonary edema, monitoring for internal bleeding, and the early detection of breast cancer.