Quantitative photoacoustic computed tomography (qPACT) holds great promise to advance a variety of important clinical applications with its potential to estimate vital physiological properties such as oxygen saturation. However, the qPACT reconstruction problem is highly nonlinear and ill-posed. Conventional spectral unmixing methods often oversimplify the problem, resulting in suboptimal accuracy. Alternatively, more principled image reconstruction approaches that comprehensively model the imaging physics are computationally burdensome and require the design of effective regularization strategies. To overcome these limitations, learning-based methods have been proposed. To date, however, the effectiveness of such methods on full-scale problems in which clinically relevant variability in anatomy and physiological parameters is considered has not been established. To address this, we investigated the use of a convolutional neural network with spatial and channel attention modules to estimate the three-dimensional (3D) distribution of tissue oxygenation within vessels and lesions in the female breast. The network was provided with input data comprising noise-corrupted 3D initial pressure distributions corresponding to three wavelengths (757, 800, 850 nm). An additional novel aspect of our study was the use of realistic 3D numerical breast phantoms that described stochastic variations in breast anatomy and functional properties, which enabled a meaningful, quantitative, and systematic evaluation of the proposed method. This study represents an important contribution to the field of qPACT and will guide the exploration of learning-based methods to help translate this important technology by delineating potential prospects and limitations.
KEYWORDS: Image restoration, Photoacoustic tomography, Tunable filters, Data modeling, Education and training, Data acquisition, Signal filtering, Linear filtering, Image filtering, 3D modeling
Photoacoustic computed tomography (PACT) is being actively developed for breast cancer imaging. In 3D PACT imagers for breast imaging, a hemispherical measurement geometry that encloses the breast has been employed. Such measurement data are referred to as “half-scan” data. Existing closed-form reconstruction methods assume a closed measurement aperture; however, the direct application of these methods to half-scan data results in inaccurate images with artifacts. Previous studies have demonstrated that half-scan data are “complete” in the sense that they contain sufficient information for accurate and stable reconstruction of an object contained within a hemispherical measurement aperture. However, direct closed-form methods for use with half-scan data have not been reported. Although optimization-based iterative image reconstruction methods are applicable, they are computationally intensive. In this work, for the first time, a semi-analytic image reconstruction method of filtered backprojection (FBP) form was proposed for use with half-scan PACT data. To accomplish this, the unknown data filtering operation is learned in a data-driven way by use of a linear U-Net neural network. To investigate the method, stochastic 3D numerical breast phantoms (NBPs) were used for model training and testing. As a result of the completeness of the half-scan data, we demonstrate that the learned FBP method can be widely applied, even when the to-be-reconstructed object differs considerably from those that were used to learn the data filtering. This is a key feature of the method that will enable it to have an important practical impact on PACT.
SignificanceDynamic photoacoustic computed tomography (PACT) is a valuable imaging technique for monitoring physiological processes. However, current dynamic PACT imaging techniques are often limited to two-dimensional spatial imaging. Although volumetric PACT imagers are commercially available, these systems typically employ a rotating measurement gantry in which the tomographic data are sequentially acquired as opposed to being acquired simultaneously at all views. Because the dynamic object varies during the data-acquisition process, the sequential data-acquisition process poses substantial challenges to image reconstruction associated with data incompleteness. The proposed image reconstruction method is highly significant in that it will address these challenges and enable volumetric dynamic PACT imaging with existing preclinical imagers.AimThe aim of this study is to develop a spatiotemporal image reconstruction (STIR) method for dynamic PACT that can be applied to commercially available volumetric PACT imagers that employ a sequential scanning strategy. The proposed reconstruction method aims to overcome the challenges caused by the limited number of tomographic measurements acquired per frame.ApproachA low-rank matrix estimation-based STIR (LRME-STIR) method is proposed to enable dynamic volumetric PACT. The LRME-STIR method leverages the spatiotemporal redundancies in the dynamic object to accurately reconstruct a four-dimensional (4D) spatiotemporal image.ResultsThe conducted numerical studies substantiate the LRME-STIR method’s efficacy in reconstructing 4D dynamic images from tomographic measurements acquired with a rotating measurement gantry. The experimental study demonstrates the method’s ability to faithfully recover the flow of a contrast agent with a frame rate of 10 frames per second, even when only a single tomographic measurement per frame is available.ConclusionsThe proposed LRME-STIR method offers a promising solution to the challenges faced by enabling 4D dynamic imaging using commercially available volumetric PACT imagers. By enabling accurate STIRs, this method has the potential to significantly advance preclinical research and facilitate the monitoring of critical physiological biomarkers.
SignificanceWhen developing a new quantitative optoacoustic computed tomography (OAT) system for diagnostic imaging of breast cancer, objective assessments of various system designs through human trials are infeasible due to cost and ethical concerns. In prototype stages, however, different system designs can be cost-efficiently assessed via virtual imaging trials (VITs) employing ensembles of digital breast phantoms, i.e., numerical breast phantoms (NBPs), that convey clinically relevant variability in anatomy and optoacoustic tissue properties.AimThe aim is to develop a framework for generating ensembles of realistic three-dimensional (3D) anatomical, functional, optical, and acoustic NBPs and numerical lesion phantoms (NLPs) for use in VITs of OAT applications in the diagnostic imaging of breast cancer.ApproachThe generation of the anatomical NBPs was accomplished by extending existing NBPs developed by the U.S. Food and Drug Administration. As these were designed for use in mammography applications, substantial modifications were made to improve blood vasculature modeling for use in OAT. The NLPs were modeled to include viable tumor cells only or a combination of viable tumor cells, necrotic core, and peripheral angiogenesis region. Realistic optoacoustic tissue properties were stochastically assigned in the NBPs and NLPs.ResultsTo advance optoacoustic and optical imaging research, 84 datasets have been released; these consist of anatomical, functional, optical, and acoustic NBPs and the corresponding simulated multi-wavelength optical fluence, initial pressure, and OAT measurements. The generated NBPs were compared with clinical data with respect to the volume of breast blood vessels and spatially averaged effective optical attenuation. The usefulness of the proposed framework was demonstrated through a case study to investigate the impact of acoustic heterogeneity on OAT images of the breast.ConclusionsThe proposed framework will enhance the authenticity of virtual OAT studies and can be widely employed for the investigation and development of advanced image reconstruction and machine learning-based methods, as well as the objective evaluation and optimization of the OAT system designs.
Previous studies of dynamic photoacoustic computed tomography (PACT) consider the case where complete data can be rapidly acquired and employed to directly reconstruct a sequence of images. However, such frame-by-frame methods do not apply to commercially available volumetric PACT imaging systems with rotating gantries because the object varies during data acquisition. Furthermore, the rotation speed and the laser repetition rate limit the number of tomographic views per frame. In this study, a low-rank matrix estimation-based spatiotemporal image reconstruction method attuned to rotating-gantry volumetric PACT systems is proposed, and its accuracy is shown by numerical and experimental studies.
The ability to perform dynamic imaging of time-varying physiological processes in small animal models is critically needed to understand the progression of human diseases and develop new therapies. Photoacoustic computed tomography (PACT) has been recognized as a promising tool for small animal imaging because of its relatively low expense, high resolution, and good signal-to-noise ratio. By exploiting the optical absorption of hemoglobin or exogenous contrast agents, dynamic PACT holds excellent potential for measuring important time-varying biomarkers like tumor vascular perfusion. Nonetheless, current dynamic PACT technologies possess several limitations. Most three-dimensional (3D) PACT imagers employ a tomographic measurement process in which a gantry containing acoustic transducers is rotated about the animal. Such a rotating gantry is advantageous for limiting the cost of the system due to the decreased number of acoustic transducers and associated electronics and for enabling convenient delivery of the light to the object. However, this presents significant challenges for dynamic image reconstruction because only a few tomographic views are available to reconstruct each temporal frame. This work presents an efficient and accurate dynamic image reconstruction method that can be deployed with widely available 3D imagers using rotating gantries. In particular, a low-rank matrix estimation based spatiotemporal image reconstruction (LRME-STIR) algorithm is proposed. In a stylized virtual dynamic contrast-enhanced imaging study, the proposed LRME-STIR algorithm is shown to accurately recover a well characterized dynamic numerical murine phantom in which tumor vascular perfusion and breathing motion are modeled.
The circular Radon transform (CRT) is widely employed as an imaging model for wave-based tomographic bioimaging modalities like ultrasound reflectivity tomography. A complete set of CRT data function is known to have redundancies. However, no explicit non-iterative image reconstruction method is known for inverting temporally-truncated data. To address this, a learning-based approach is proposed to establish a filtered backprojection (FBP) method for use with the half-time CRT data function. The proposed method approximates a mapping that is known to exist in theory; therefore, it is fundamentally different than many deep-learning based reconstruction methods that seek to establish a non-existent mapping. Thus, the proposed method performs well on unforeseen data. The learned half-time FBP achieves image quality comparable to a conventional full-time FBP method although it uses half of the complete data.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.