Visual Quality Enhancement in Optoacoustic Tomography Using Active Contour Segmentation Priors

Artificial Intelligence Assisted Multi-modal Photoacoustic-Ultrasound Imaging for Studying Renal Tissue Function and Hemodynamics

Combined functional-anatomic imaging modalities, which integrate the benefits of visualizing gross anatomy along with the functional or metabolic information of tissue has revolutionized the world of medical imaging. However, such existing imaging modalities are very costly. An alternative option could be a hybrid modality combining contrast-enhanced ultrasound, doppler and photoacoustic imaging. In the current study, we propose an artificial intelligence assisted multi-modal imaging platform where we have used U-net model for segmenting the anatomical features from the ultrasound images obtained from an animal model study. The neural network has performed accurately for three different cases, each with a high dice score. The model was co-validated with doppler images. Further, blood perfusion and tissue oxygenation information from the predicted anatomical structures were also studied. The present findings confirm the feasibility of using this multimodal imaging modality facilitated by artificial intelligence for better understanding of the hemodynamics of the kidney.Clinical Relevance-A multi-modal imaging technique has been proposed which would provide anatomical and functional information to the clinicians for early detection and tracking of the disease prognosis. Unlike existing imaging modalities like PET-CT (Positron Emission Tomography- Computed Tomography), the proposed modality is much more costeffective and radiation free (non-ionizing nature).

A Signal Domain Object Segmentation Method for Ultrasound and Photoacoustic Computed Tomography

Image segmentation is important in improving the diagnostic capability of ultrasound computed tomography (USCT) and photoacoustic computed tomography (PACT), as it can be included in the image reconstruction process to improve image quality and quantification abilities. Segmenting the imaged object out of the background using image domain methods is easily complicated by low contrast, noise, and artifacts in the reconstructed image. Here, we introduce a new signal domain object segmentation method for USCT and PACT which does not require image reconstruction beforehand and is automatic, robust, computationally efficient, accurate, and straightforward. We first establish the relationship between the time-of-flight of the received first arrival waves and the object's boundary which is described by ellipse equations. Then, we show that the ellipses are tangent to the boundary. By looking for tangent points on the common tangent of neighboring ellipses, the boundary can be approximated with high fidelity. Imaging experiments of human fingers and mice cross-sections showed that our method provided equivalent or better segmentations than the optimal ones by active contours. In summary, our method greatly reduces the overall complexity of object segmentation and shows great potential in eliminating user dependency without sacrificing segmentation accuracy. The method can be further seamlessly incorporated into algorithms for other processing purposes in USCT and PACT, such as high-quality image reconstruction.

Pixel-wise reconstruction of tissue absorption coefficients in photoacoustic tomography using a non-segmentation iterative method

In Photoacoustic Tomography (PAT), the acquired image represents a light energy deposition map of the imaging object. For quantitative imaging, the PAT image is converted into an absorption coefficient ( μ a ) map by dividing the light fluence (LF). Previous methods usually assume a uniform tissue μ a distribution, and consequently degrade the LF correction results. Here, we propose a simple method to reconstruct the pixel-wise μ a map. Our method is based on a non-segmentation-based iterative algorithm, which alternately optimizes the LF distribution and the μ a map. Using simulation data, as well as phantom and animal data, we implemented our algorithm and compared it to segmentation-based correction methods. The results show that our method can obtain accurate estimation of the LF distribution and therefore improve the image quality and feature visibility of the μ a map. Our method may facilitate efficient calculation of the concentration distributions of endogenous and exogenous agents in vivo.

MRI Information-Based Correction and Restoration of Photoacoustic Tomography

As an emerging molecular imaging modality, Photoacoustic Tomography (PAT) is capable of mapping tissue physiological metabolism and exogenous contrast agent information with high specificity. Due to its ultrasonic detection mechanism, the precise localization of targeted lesions has long been a challenge for PAT imaging. The poor soft-tissue contrast of the PAT image makes this process difficult and inaccurate. To meet this challenge, in this study, we first make use of the rich and clear structural information brought about by another advanced imaging modality, Magnetic Resonance Imaging (MRI), to assist organ segmentation and correct for the light fluence attenuation of PAT. We demonstrate improved feature visibility and enhanced localization of endogenous and exogenous agents in the fluence corrected PAT images. Compared with PAT-based methods, the contrast-to-noise ratio (CNR) of our MRI-assisted method increases by 29.1% in live animal experiments. Furthermore, we show that the co-registered MRI image can also be incorporated into PAT image restoration, and achieves improved anatomical landscape and soft-tissue contrast (CNR increased by 25.36%) while preserving similar spatial resolution. This PAT-MRI combination provides excellent structural, functional and molecular images of the subject, and may enable more comprehensive analysis of various preclinical research applications.

Automatic 3-D segmentation and volumetric light fluence correction for photoacoustic tomography based on optimal 3-D graph search

Preclinical imaging with photoacoustic tomography (PAT) has attracted wide attention in recent years since it is capable of providing molecular contrast with deep imaging depth. The automatic extraction and segmentation of the animal in PAT images is crucial for improving image analysis efficiency and enabling advanced image post-processing, such as light fluence (LF) correction for quantitative PAT imaging. Previous automatic segmentation methods are mostly two-dimensional approaches, which failed to conserve the 3-D surface continuity because the image slices were processed separately. This discontinuity problem further hampers LF correction, which, ideally, should be carried out in 3-D due to spatially diffused illumination. Here, to solve these problems, we propose a volumetric auto-segmentation method for small animal PAT imaging based on the 3-D optimal graph search (3-D GS) algorithm. The 3-D GS algorithm takes into account the relation among image slices by constructing a 3-D node-weighted directed graph, and thus ensures surface continuity. In view of the characteristics of PAT images, we improve the original 3-D GS algorithm on graph construction, solution guidance and cost assignment, such that the accuracy and smoothness of the segmented animal surface were guaranteed. We tested the performance of the proposed method by conducting in vivo nude mice imaging experiments with a commercial preclinical cross-sectional PAT system. The results showed that our method successfully retained the continuous global surface structure of the whole 3-D animal body, as well as smooth local subcutaneous tumor boundaries at different development stages. Moreover, based on the 3-D segmentation result, we were able to simulate volumetric LF distribution of the entire animal body and obtained LF corrected PAT images with enhanced structural visibility and uniform image intensity.

Binary photoacoustic tomography for improved vasculature imaging

SIGNIFICANCE

The proposed binary tomography approach was able to recover the vasculature structures accurately, which could potentially enable the utilization of binary tomography algorithm in scenarios such as therapy monitoring and hemorrhage detection in different organs.

AIM

Photoacoustic tomography (PAT) involves reconstruction of vascular networks having direct implications in cancer research, cardiovascular studies, and neuroimaging. Various methods have been proposed for recovering vascular networks in photoacoustic imaging; however, most methods are two-step (image reconstruction and image segmentation) in nature. We propose a binary PAT approach wherein direct reconstruction of vascular network from the acquired photoacoustic sinogram data is plausible.

APPROACH

Binary tomography approach relies on solving a dual-optimization problem to reconstruct images with every pixel resulting in a binary outcome (i.e., either background or the absorber). Further, the binary tomography approach was compared against backprojection, Tikhonov regularization, and sparse recovery-based schemes.

RESULTS

Numerical simulations, physical phantom experiment, and in-vivo rat brain vasculature data were used to compare the performance of different algorithms. The results indicate that the binary tomography approach improved the vasculature recovery by 10% using in-silico data with respect to the Dice similarity coefficient against the other reconstruction methods.

CONCLUSION

The proposed algorithm demonstrates superior vasculature recovery with limited data both visually and based on quantitative image metrics.

Nonlinear iterative perturbation scheme with simplified spherical harmonics (SP3 ) light propagation model for quantitative photoacoustic tomography

When using quantitative photoacoustic tomography (q-PAT) reconstruction to recover the optical absorption coefficients of tissue, the commonly used diffusion equation has several limitations in the case of the objects that have small geometries and high-absorption or low-scattering areas. Furthermore, the conventional perturbation reconstruction strategy is unsatisfactory when the target tissue containing large heterogeneous features. We herein present a modified q-PAT implementation that employs the higher-order photon migration model achieving the tradeoff between mathematical rigidity and computational efficiency. Besides, a nonlinear iterative method is proposed to obtain the perturbations of optical absorption considering the updating of the sensitivity matrix in calculating the fluence perturbations. Consequently, the distribution of tissue optical properties can be recovered in a robust way even if the targets with high absorption are included. The proposed approach has been validated by simulation, phantom and in vivo experiments, exhibiting promising performances in image fidelity and quantitative feasibility for practical applications.

Deep Learning for Automatic Segmentation of Hybrid Optoacoustic Ultrasound (OPUS) Images

The highly complementary information provided by multispectral optoacoustics and pulse-echo ultrasound have recently prompted development of hybrid imaging instruments bringing together the unique contrast advantages of both modalities. In the hybrid optoacoustic ultrasound (OPUS) combination, images retrieved by one modality may further be used to improve the reconstruction accuracy of the other. In this regard, image segmentation plays a major role as it can aid improving the image quality and quantification abilities by facilitating modeling of light and sound propagation through the imaged tissues and surrounding coupling medium. Here, we propose an automated approach for surface segmentation in whole-body mouse OPUS imaging using a deep convolutional neural network (CNN). The method has shown robust performance, attaining accurate segmentation of the animal boundary in both optoacoustic and pulse-echo ultrasound images, as evinced by quantitative performance evaluation using Dice coefficient metrics.

Brilliant cresyl blue enhanced optoacoustic imaging enables non-destructive imaging of mammalian ovarian follicles for artificial reproduction

In the field of reproductive biology, there is a strong need for a suitable tool capable of non-destructive evaluation of oocyte viability and function. We studied the application of brilliant cresyl blue (BCB) as an intra-vital exogenous contrast agent using multispectral optoacoustic tomography (MSOT) for visualization of porcine ovarian follicles. The technique provided excellent molecular sensitivity, enabling the selection of competent oocytes without disrupting the follicles. We further conducted in vitro embryo culture, molecular analysis (real-time and reverse transcriptase polymerase chain reaction) and DNA fragmentation analysis to comprehensively establish the safety of BCB-enhanced MSOT imaging in monitoring oocyte viability. Overall, the experimental results suggest that the method offers a significant advance in the use of contrast agents and molecular imaging for reproductive studies. Our technique improves the accurate prediction of ovarian reserve significantly and, once standardized for in vivo imaging, could provide an effective tool for clinical infertility management.

Model-Based Reconstruction of Large Three-Dimensional Optoacoustic Datasets

Iterative model-based algorithms are known to enable more accurate and quantitative optoacoustic (photoacoustic) tomographic reconstructions than standard back-projection methods. However, three-dimensional (3D) model-based inversion is often hampered by high computational complexity and memory overhead. Parallel implementations on a graphics processing unit (GPU) have been shown to efficiently reduce the memory requirements by on-the-fly calculation of the actions of the optoacoustic model matrix, but the high complexity still makes these approaches impractical for large 3D optoacoustic datasets. Herein, we show that the computational complexity of 3D model-based iterative inversion can be significantly reduced by splitting the model matrix into two parts: one maximally sparse matrix containing only one entry per voxel-transducer pair and a second matrix corresponding to cyclic convolution. We further suggest reconstructing the images by multiplying the transpose of the model matrix calculated in this manner with the acquired signals, which is equivalent to using a very large regularization parameter in the iterative inversion method. The performance of these two approaches is compared to that of standard back-projection and a recently introduced GPU-based model-based method using datasets from in vivo experiments. The reconstruction time was accelerated by approximately an order of magnitude with the new iterative method, while multiplication with the transpose of the matrix is shown to be as fast as standard back-projection.

Optoacoustic image formation approaches-a clinical perspective

Clinical translation of optoacoustic imaging is fostered by the rapid technical advances in imaging performance as well as the growing number of clinicians recognizing the immense diagnostic potential of this technology. Clinical optoacoustic systems are available in multiple configurations, including hand-held and endoscopic probes as well as raster-scan approaches. The hardware design must be adapted to the accessible portion of the imaged region and other application-specific requirements pertaining the achievable depth, field of view or spatio-temporal resolution. Equally important is the adequate choice of the signal and image processing approach, which is largely responsible for the resulting imaging performance. Thus, new image reconstruction algorithms are constantly evolving in parallel to the newly-developed set-ups. This review focuses on recent progress on optoacoustic image formation algorithms and processing methods in the clinical setting. Major reconstruction challenges include real-time image rendering in two and three dimensions, efficient hybridization with other imaging modalitites as well as accurate interpretation and quantification of bio-markers, herein discussed in the context of ongoing progress in clinical translation.

Maximum Entropy Based Non-Negative Optoacoustic Tomographic Image Reconstruction

OBJECTIVE

Optoacoustic (photoacoustic) tomography is aimed at reconstructing maps of the initial pressure rise induced by the absorption of light pulses in tissue. In practice, due to inaccurate assumptions in the forward model, noise, and other experimental factors, the images are often afflicted by artifacts, occasionally manifested as negative values. The aim of this work is to develop an inversion method which reduces the occurrence of negative values and improves the quantitative performance of optoacoustic imaging.

METHODS

We present a novel method for optoacoustic tomography based on an entropy maximization algorithm, which uses logarithmic regularization for attaining non-negative reconstructions. The reconstruction image quality is further improved using structural prior-based fluence correction.

RESULTS

We report the performance achieved by the entropy maximization scheme on numerical simulation, experimental phantoms, and in-vivo samples.

CONCLUSION

The proposed algorithm demonstrates superior reconstruction performance by delivering non-negative pixel values with no visible distortion of anatomical structures.

SIGNIFICANCE

Our method can enable quantitative optoacoustic imaging, and has the potential to improve preclinical and translational imaging applications.

Transmission-reflection optoacoustic ultrasound (TROPUS) computed tomography of small animals

Rapid progress in the development of multispectral optoacoustic tomography techniques has enabled unprecedented insights into biological dynamics and molecular processes in vivo and noninvasively at penetration and spatiotemporal scales not covered by modern optical microscopy methods. Ultrasound imaging provides highly complementary information on elastic and functional tissue properties and further aids in enhancing optoacoustic image quality. We devised the first hybrid transmission-reflection optoacoustic ultrasound (TROPUS) small animal imaging platform that combines optoacoustic tomography with both reflection- and transmission-mode ultrasound computed tomography. The system features full-view cross-sectional tomographic imaging geometry for concomitant noninvasive mapping of the absorbed optical energy, acoustic reflectivity, speed of sound, and acoustic attenuation in whole live mice with submillimeter resolution and unrivaled image quality. Graphics-processing unit (GPU)-based algorithms employing spatial compounding and bent-ray-tracing iterative reconstruction were further developed to attain real-time rendering of ultrasound tomography images in the full-ring acquisition geometry. In vivo mouse imaging experiments revealed fine details on the organ parenchyma, vascularization, tissue reflectivity, density, and stiffness. We further used the speed of sound maps retrieved by the transmission ultrasound tomography to improve optoacoustic reconstructions via two-compartment modeling. The newly developed synergistic multimodal combination offers unmatched capabilities for imaging multiple tissue properties and biomarkers with high resolution, penetration, and contrast.

Optoacoustic image reconstruction: the full inverse problem with variable bases

Optoacoustic imaging was for a long time concerned with the reconstruction of energy density or optical properties. In this work, we present the full inverse problem with respect to optical absorption and diffusion as well as speed of sound and mass density. The inverse problem is solved by an iterative gradient-based optimization procedure. Since the ill-conditioning increases with the number of sought parameters, we propose two approaches to improve the conditioning. The first approach is based on the reduction of the size of the basis for the parameter spaces, that evolves according to the particular characteristics of the solution, while maintaining the flexibility of element-wise parameter selection. The second approach is a material identification technique that incorporates prior knowledge of expected material types and uses the acoustical gradients to identify materials uniquely. We present numerical studies to illustrate the properties and functional principle of the proposed methods. Significant convergence speed-ups are gained by the two approaches countering ill-conditioning. Additionally, we show results for the reconstruction of a mouse brain from in vivo measurements.

Segmentation of vessel structures from photoacoustic images with reliability assessment

Photoacoustic imaging enables the imaging of soft biological tissue with combined optical contrast and ultrasound resolution. One of the targets of interest is tissue vasculature. However, the photoacoustic images may not directly provide the information on, for example, vasculature structure. Therefore, the images are improved by reducing noise and artefacts, and processed for better visualisation of the target of interest. In this work, we present a new segmentation method of photoacoustic images that also straightforwardly produces assessments of its reliability. The segmentation depends on parameters which have a natural tendency to increase the reliability as the parameter values monotonically change. The reliability is assessed by counting classifications of image voxels with different parameter values. The resulting segmentation with reliability offers new ways and tools to analyse photoacoustic images and new possibilities for utilising them as anatomical priors in further computations. Our MATLAB implementation of the method is available as an open-source software package [P. Raumonen, Matlab, 2018].

Photoacoustic-based sO2 estimation through excised bovine prostate tissue with interstitial light delivery

Photoacoustic (PA) imaging is capable of probing blood oxygen saturation (sO₂), which has been shown to correlate with tissue hypoxia, a promising cancer biomarker. However, wavelength-dependent local fluence changes can compromise sO₂ estimation accuracy in tissue. This work investigates using PA imaging with interstitial irradiation and local fluence correction to assess precision and accuracy of sO₂ estimation of blood samples through ex vivo bovine prostate tissue ranging from 14% to 100% sO₂. Study results for bovine blood samples at distances up to 20 mm from the irradiation source show that local fluence correction improved average sO₂ estimation error from 16.8% to 3.2% and maintained an average precision of 2.3% when compared to matched CO-oximeter sO₂ measurements. This work demonstrates the potential for future clinical translation of using fluence-corrected and interstitially driven PA imaging to accurately and precisely assess sO₂ at depth in tissue with high resolution.

Constrained Inversion and Spectral Unmixing in Multispectral Optoacoustic Tomography

Accurate extraction of physical and biochemical parameters from optoacoustic images is often impeded due to the use of unrigorous inversion schemes, incomplete tomographic detection coverage, or other experimental factors that cannot be readily accounted for during the image acquisition and reconstruction process. For instance, inaccurate assumptions in the physical forward model may lead to negative optical absorption values in the reconstructed images. Any artifacts present in the single wavelength optoacoustic images can be significantly aggravated when performing a two-step reconstruction consisting in acoustic inversion and spectral unmixing aimed at rendering the distributions of spectrally distinct absorbers. We investigate a number of algorithmic strategies with non-negativity constraints imposed at the different phases of the reconstruction process. Performance is evaluated in cross-sectional multispectral optoacoustic tomography recordings from tissue-mimicking phantoms and in vivo mice embedded with varying concentrations of contrast agents. Additional in vivo validation is subsequently performed with molecular imaging data involving subcutaneous tumors labeled with genetically expressed iRFP proteins and organ perfusion by optical contrast agents. It is shown that constrained reconstruction is essential for reducing the critical image artifacts associated with inaccurate modeling assumptions. Furthermore, imposing the non-negativity constraint directly on the unmixed distribution of the probe of interest was found to maintain the most robust and accurate reconstruction performance in all experiments.