Samuli Siltanen

Zenith Purisha, Sakari S. Karhula, Juuso Ketola et al.

X-ray tomography is a reliable tool for determining the inner structure of 3D object with penetrating X-rays. However, traditional reconstruction methods such as 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 available memory. However, computing 3D reconstruction from such a sparsely sampled dataset is very sensitive to measurement noise and modelling errors. An automatic regularization method is proposed for robust reconstruction, based on enforcing sparsity in the three-dimensional shearlet transform domain. The inputs of the algorithm are the projection data and {\it a priori} known expected degree of sparsity, denoted $0<{\mathcal C}_{pr}\leq 1$. The number ${\mathcal C}_{pr}$ 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 while retaining the quality of the reconstruction images and of the morphometric paramaters.

Sarah Jane Hamilton, Samuli Siltanen

Electrical impedance tomography (EIT) is a non-invasive imaging method in which an unknown physical body is probed with electric currents applied on the boundary, and the internal conductivity distribution is recovered from the measured boundary voltage data. The reconstruction task is a nonlinear and ill-posed inverse problem, whose solution calls for special regularized algorithms, such as D-bar methods which are based on complex geometrical optics solutions (CGOs). In many applications of EIT, such as monitoring the heart and lungs of unconscious intensive care patients or locating the focus of an epileptic seizure, data acquisition on the entire boundary of the body is impractical, restricting the boundary area available for EIT measurements. An extension of the D-bar method to the case when data is collected only on a subset of the boundary is studied by computational simulation. The approach is based on solving a boundary integral equation for the traces of the CGOs using localized basis functions (Haar wavelets). The numerical evidence suggests that the D-bar method can be applied to partial-boundary data in dimension two and that the traces of the partial data CGOs approximate the full data CGO solutions on the available portion of the boundary, for the necessary small $k$ frequencies.

Samuli Siltanen, Janne P. Tamminen

The aim of electrical impedance tomography is to form an image of the conductivity distribution inside an unknown body using electric boundary measurements. The computation of the image from measurement data is a non-linear ill-posed inverse problem and calls for a special regularized algorithm. One such algorithm, the so-called D-bar method, is improved in this work by introducing new computational steps that remove the so far necessary requirement that the conductivity should be constant near the boundary. The numerical experiments presented suggest two conclusions. First, for most conductivities arising in medical imaging, it seems the previous approach of using a best possible constant near the boundary is sufficient. Second, for conductivities that have high contrast features at the boundary, the new approach produces reconstructions with smaller quantitative error and with better visual quality.

Janne P. Tamminen, Tanja Tarvainen, Samuli Siltanen

The D-bar method at negative energy is numerically implemented. Using the method we are able to numerically reconstruct potentials and investigate exceptional points at negative energy. Subsequently, applying the method to Diffusive Optical Tomography, a new way of reconstructing the diffusion coefficient from the associated Complex Geometrics Optics solution is suggested and numerically validated.

Saara Isoranta, Emilia L. K. Blåsten, Lílian Ferreira de Freitas et al.

We approach image denoising from a perception-driven perspective: how can we select the parameters that are best suited for human visual perception? We combine research methods in mathematics and psychology to develop a mathematical framework for measuring perceived similarity. We construct a sample set of differently denoised photographs by using the same base image as input data and by tuning the parameter value in a total variation denoising algorithm. A comparison test is conducted with human participants to survey perceived differences between the images. Analyzing the results with psychometric scaling provides us with a HaarPSI value to use as a threshold in discretizing parameter grids. As a result, we obtain psychometrically scaled, openly available image sets that are ready to use in further experiments in perception-driven imaging, as well as a framework for ensuing experiments involving comparison tests.

Sara Sippola, Siiri Rautio, Andreas Hauptmann et al.

Electrical impedance tomography (EIT) is a non-invasive imaging method with diverse applications, including medical imaging and non-destructive testing. The inverse problem of reconstructing internal electrical conductivity from boundary measurements is nonlinear and highly ill-posed, making it difficult to solve accurately. In recent years, there has been growing interest in combining analytical methods with machine learning to solve inverse problems. In this paper, we propose a method for estimating the convex hull of inclusions from boundary measurements by combining the enclosure method proposed by Ikehata with neural networks. We demonstrate its performance using experimental data. Compared to the classical enclosure method with least squares fitting, the learned convex hull achieves superior performance on both simulated and experimental data.

Juan Pablo Agnelli, Fernando S. Moura, Siiri Rautio et al.

Electrical impedance tomography (EIT) is a non-invasive imaging method for recovering the internal conductivity of a physical body from electric boundary measurements. EIT combined with machine learning has shown promise for the classification of strokes. However, most previous works have used raw EIT voltage data as network inputs. We build upon a recent development which suggested the use of special noise-robust Virtual Hybrid Edge Detection (VHED) functions as network inputs, although that work used only highly simplified and mathematically ideal models. In this work we strengthen the case for the use of EIT, and VHED functions especially, for stroke classification. We design models with high detail and mathematical realism to test the use of VHED functions as inputs. Virtual patients are created using a physically detailed 2D head model which includes features known to create challenges in real-world imaging scenarios. Conductivity values are drawn from statistically realistic distributions, and phantoms are afflicted with either hemorrhagic or ischemic strokes of various shapes and sizes. Simulated noisy EIT electrode data, generated using the realistic Complete Electrode Model (CEM) as opposed to the mathematically ideal continuum model, is processed to obtain VHED functions. We compare the use of VHED functions as inputs against the alternative paradigm of using raw EIT voltages. Our results show that (i) stroke classification can be performed with high accuracy using 2D EIT data from physically detailed and mathematically realistic models, and (ii) in the presence of noise, VHED functions outperform raw data as network inputs.

Markus Juvonen, Samuli Siltanen, Fernando Silva de Moura

The photographic dataset collected for the Helsinki Deblur Challenge 2021 (HDC2021) contains pairs of images taken by two identical cameras of the same target but with different conditions. One camera is always in focus and produces sharp and low-noise images the other camera produces blurred and noisy images as it is gradually more and more out of focus and has a higher ISO setting. Even though the dataset was designed and captured with the HDC2021 in mind it can be used for any testing and benchmarking of image deblurring algorithms. The data is available here: https://doi.org/10.5281/zenodo.477228

Tatiana A. Bubba, Mathilde Galinier, Matti Lassas et al.

We propose a novel convolutional neural network (CNN), called $Ψ$DONet, designed for learning pseudodifferential operators ($Ψ$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 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 $Ψ$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 $Ψ$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 $Ψ$DOs or Fourier integral operators.

David Villacis, Santeri Kaupinmäki, Samuli Siltanen et al.

This is a photographic dataset collected for testing image processing algorithms. The idea is to have images that can exploit the properties of total variation, therefore a set of playing cards was distributed on the scene. The dataset is made available at www.fips.fi/photographic_dataset2.php

Teemu Helenius, Samuli Siltanen

This is a photographic dataset collected for testing image processing algorithms. The idea is to have sets of different but statistically similar images. In this work the images show randomly distributed peppercorns. The dataset is made available at www.fips.fi/photographic_dataset.php .

Kati Niinimäki, Matti Lassas, Keijo Hämäläinen et al.

A computational method is introduced for choosing the regularization parameter for total variation (TV) regularization. The 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 X-ray tomography data measured from a walnut and compared to the S-curve method. The proposed method seems to automatically adapt to the desired resolution and noise level, and it yields useful results in the tests. The results are comparable to those of the S-curve method; however, the S-curve method needs a priori information about the sparsity of the unknown, while the proposed method does not need any a priori information (apart from the choice of a desired resolution). Mathematical analysis is presented for (partial) understanding of the properties of the proposed parameter choice method. It is rigorously proven that the TV norms of the reconstructions converge with any choice of regularization parameter.

Sarah Hamilton, Juan Manuel Reyes, Samuli Siltanen et al.

The Regularized D-bar method for Electrical Impedance Tomography 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 approach that combines the rigor of the D-bar approach with the edge-preserving nature of Total Variation 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 while still solving the full nonlinear problem. 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.