Helle Majander

Antti Hannukainen, Nuutti Hyvönen, Helle Majander et al.

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating the distribution of optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. Photoacoustic images are three-dimensional and high-resolution. Furthermore, high-resolution reconstructions of the optical parameters are targeted. Therefore, in order to provide a practical method for quantitative photoacoustic tomography, the inversion algorithm needs to be able to perform successfully with problems of prominent size. In this work, an efficient approach for the inverse problem of quantitative photoacoustic tomography is proposed, assuming an edge-preferring prior for the optical parameters. The method is based on iteratively combining priorconditioned LSQR with a lagged diffusivity step and a linearisation of the measurement model, with the needed multiplications by Jacobians performed in a matrix-free manner. The algorithm is tested with three-dimensional numerical simulations. The results show that the approach can be used to produce accurate and high quality estimates of absorption and diffusion in complex three-dimensional geometries with moderate computation time and cost.

Nuutti Hyvönen, Helle Majander, Stratos Staboulis

Electrical impedance tomography aims at reconstructing the conductivity inside a physical body from boundary measurements of current and voltage at a finite number of contact electrodes. In many practical applications, the shape of the imaged object is subject to considerable uncertainties that render reconstructing the internal conductivity impossible if they are not taken into account. This work numerically demonstrates that one can compensate for inaccurate modeling of the object boundary in two spatial dimensions by estimating the locations and sizes of the electrodes as a part of a reconstruction algorithm. The numerical studies, which are based on both simulated and experimental data, are complemented by proving that the employed complete electrode model is approximately conformally invariant, which suggests that the obtained reconstructions in mismodeled domains reflect conformal images of the true targets. The numerical experiments also confirm that a similar approach does not, in general, lead to a functional algorithm in three dimensions.

Antti Hannukainen, Lauri Harhanen, Nuutti Hyvönen et al.

In optical tomography a physical body is illuminated with near-infrared light and the resulting outward photon flux is measured at the object boundary. The goal is to reconstruct internal optical properties of the body, such as absorption and diffusivity. In this work, it is assumed that the imaged object is composed of an approximately homogeneous background with clearly distinguishable embedded inhomogeneities. An algorithm for finding the maximum a posteriori estimate for the absorption and diffusion coefficients is introduced assuming an edge-preferring prior and an additive Gaussian measurement noise model. The method is based on iteratively combining a lagged diffusivity step and a linearization of the measurement model of diffuse optical tomography with priorconditioned LSQR. The performance of the reconstruction technique is tested via three-dimensional numerical experiments with simulated measurement data.