Image Reconstruction In Optoacoustic Tomography Research Articles

Invariant representations in deep learning for optoacoustic imaging

Image reconstruction in optoacoustic tomography (OAT) is a trending learning task highly dependent on measured physical magnitudes present at sensing time. A large number of different settings and also the presence of uncertainties or partial knowledge of parameters can lead to reconstruction algorithms that are specifically tailored and designed to a particular configuration, which could not be the one that will ultimately be faced in a final practical situation. Being able to learn reconstruction algorithms that are robust to different environments (e.g., the different OAT image reconstruction settings) or invariant to such environments is highly valuable because it allows us to focus on what truly matters for the application at hand and discard what are considered spurious features. In this work, we explore the use of deep learning algorithms based on learning invariant and robust representations for the OAT inverse problem. In particular, we consider the application of the ANDMask scheme due to its easy adaptation to the OAT problem. Numerical experiments are conducted showing that when out-of-distribution generalization (against variations in parameters such as the location of the sensors) is imposed, there is no degradation of the performance and, in some cases, it is even possible to achieve improvements with respect to standard deep learning approaches where invariance robustness is not explicitly considered.

Combining band-frequency separation and deep neural networks for optoacoustic imaging

In this paper we consider the problem of image reconstruction in optoacoustic tomography. In particular, we devise a deep neural architecture that can explicitly take into account the band-frequency information contained in the sinogram. This is accomplished by two means. First, we jointly use a linear filtered back-projection method and a fully dense UNet for the generation of the images corresponding to each one of the frequency bands considered in the separation. Secondly, in order to train the model, we introduce a special loss function consisting of three terms: (i) a separating frequency bands term; (ii) a sinogram-based consistency term and (iii) a term that directly measures the quality of image reconstruction and which takes advantage of the presence of ground-truth images present in training dataset. Numerical experiments show that the proposed model, which can be easily trainable by standard optimization methods, presents an excellent generalization performance quantified by a number of metrics commonly used in practice. Also, in the testing phase, our solution has a comparable (in some cases lower) computational complexity, which is a desirable feature for real-time implementation of optoacoustic imaging.

Impact of nonstationary optical illumination on image reconstruction in optoacoustic tomography.

Optoacoustic tomography (OAT), also known as photoacoustic tomography, is a rapidly emerging hybrid imaging technique that possesses great potential for a wide range of biomedical imaging applications. In OAT, a laser is employed to illuminate the tissue of interest and acoustic signals are produced via the photoacoustic effect. From these data, an estimate of the distribution of the absorbed optical energy density within the tissue is reconstructed, referred to as the object function. This quantity is defined, in part, by the distribution of light fluence within the tissue that is established by the laser source. When performing three-dimensional imaging of large objects, such as a female human breast, it can be difficult to achieve a relatively uniform coverage of light fluence within the volume of interest when the position of the laser source is fixed. To circumvent this, researchers have proposed illumination schemes in which the relative position of the laser source and ultrasound probe is fixed, and both are rotated together to acquire a tomographic dataset. A problem with this rotating-illumination scheme is that the tomographic data are inconsistent; namely, the acoustic data recorded at each tomographic view angle (i.e., probe position) are produced by a distinct object function. In this work, the impact of this data inconsistency on image reconstruction accuracy is investigated systematically. This is accomplished by use of computer-simulation studies and application of mathematical results from the theory of microlocal analysis. These studies specify the set of image discontinuities that can be stably reconstructed with a nonstationary optical illumination setup. The study also includes a comparison of the ability of iterative and analytic image reconstruction methods to mitigate artifacts attributable to the data inconsistency.

Real-Time Model-Based Inversion in Cross-Sectional Optoacoustic Tomography.

Analytical (closed-form) inversion schemes have been the standard approach for image reconstruction in optoacoustic tomography due to their fast reconstruction abilities and low memory requirements. Yet, the need for quantitative imaging and artifact reduction has led to the development of more accurate inversion approaches, which rely on accurate forward modeling of the optoacoustic wave generation and propagation. In this way, multiple experimental factors can be incorporated, such as the exact detection geometry, spatio-temporal response of the transducers, and acoustic heterogeneities. The model-based inversion commonly results in very large sparse matrix formulations that require computationally extensive and memory demanding regularization schemes for image reconstruction, hindering their effective implementation in real-time imaging applications. Herein, we introduce a new discretization procedure for efficient model-based reconstructions in two-dimensional optoacoustic tomography that allows for parallel implementation on a graphics processing unit (GPU) with a relatively low numerical complexity. By on-the-fly calculation of the model matrix in each iteration of the inversion procedure, the new approach results in imaging frame rates exceeding 10 Hz, thus enabling real-time image rendering using the model-based approach.

Discrete Imaging Models for Three-Dimensional Optoacoustic Tomography Using Radially Symmetric Expansion Functions

Optoacoustic tomography (OAT), also known as photoacoustic tomography, is an emerging computed biomedical imaging modality that exploits optical contrast and ultrasonic detection principles. Iterative image reconstruction algorithms that are based on discrete imaging models are actively being developed for OAT due to their ability to improve image quality by incorporating accurate models of the imaging physics, instrument response, and measurement noise. In this work, we investigate the use of discrete imaging models based on Kaiser-Bessel window functions for iterative image reconstruction in OAT. A closed-form expression for the pressure produced by a Kaiser-Bessel function is calculated, which facilitates accurate computation of the system matrix. Computer-simulation and experimental studies are employed to demonstrate the potential advantages of Kaiser-Bessel function-based iterative image reconstruction in OAT.

Image reconstruction in optoacoustic tomography for dispersive acoustic media

Conventional image reconstruction methods for optoacoustic tomography (OAT) assume an idealized, non-dispersive acoustic medium. However, the linear attenuation coefficient and the phase velocity of acoustic waves propagating in soft tissue depend on temporal frequency and satisfy a known dispersion law. These frequency-dependent effects are incorporated into an optoacoustic wave equation, and a corresponding reconstruction method for OAT is developed. The improvement in image fidelity that can be achieved over conventional reconstruction methods is demonstrated by use of computer-simulation studies.