Spatiotemporal Coherence Weighting for In Vivo Cardiac Photoacoustic Image Beamformation

Improvement of LED-based photoacoustic imaging using lag-coherence factor (LCF) beamforming

BACKGROUND

Owing to its portability, affordability, and energy-efficiency, LED-based photoacoustic (PA) imaging is increasingly becoming popular when compared to its laser-based alternative, mainly for superficial vascular imaging applications. However, this technique suffers from low SNR and thereby limited imaging depth. As a result, visual image quality of LED-based PA imaging is not optimal, especially in sub-surface vascular imaging applications.

PURPOSE

Combination of linear ultrasound (US) probes and LED arrays are the most common implementation in LED-based PA imaging, which is currently being explored for different clinical imaging applications. Traditional delay-and-sum (DAS) is the most common beamforming algorithm in linear array-based PA detection. Side-lobes and reconstruction-related artifacts make the DAS performance unsatisfactory and poor for a clinical-implementation. In this work, we explored a new weighting-based image processing technique for LED-based PAs to yield improved image quality when compared to the traditional methods.

METHODS

We are proposing a lag-coherence factor (LCF), which is fundamentally based on the combination of the spatial auto-correlation of the detected PA signals. In LCF, the numerator contains lag-delay-multiply-and-sum (DMAS) beamformer instead of a conventional DAS beamformer. A spatial auto-correlation operation is performed between the detected US array signals before using DMAS beamformer. We evaluated the new method on both tissue-mimicking phantom (2D) and human volunteer imaging (3D) data acquired using a commercial LED-based PA imaging system.

RESULTS

Our novel correlation-based weighting technique showed LED-based PA image quality improvement when it is combined with conventional DAS beamformer. Both phantom and human volunteer imaging results gave a direct confirmation that by introducing LCF, image quality was improved and this method could reduce side-lobes and artifacts when compared to the DAS and coherence-factor (CF) approaches. Signal-to-noise ratio, generalized contrast-to-noise ratio, contrast ratio and spatial resolution were evaluated and compared with conventional beamformers to assess the reconstruction performance in a quantitative way. Results show that our approach offered image quality enhancement with an average signal-to-noise ratio and spatial resolution improvement of around 20% and 25% respectively, when compared with conventional CF based DAS algorithm.

CONCLUSIONS

Our results demonstrate that the proposed LCF based algorithm performs better than the conventional DAS and CF algorithms by improving signal-to-noise ratio and spatial resolution. Therefore, our new weighting technique could be a promising tool to improve the performance of LED-based PA imaging and thus accelerate its clinical translation.

Signal restoration algorithm for photoacoustic imaging systems

In a photoacoustic (PA) imaging system, the detectors are bandwidth-limited. Therefore, they capture PA signals with some unwanted ripples. This limitation degrades the resolution/contrast and induces sidelobes and artifacts in the reconstructed images along the axial direction. To compensate for the limited bandwidth effect, we present a PA signal restoration algorithm, where a mask is designed to extract the signals at the absorber positions and remove the unwanted ripples. This restoration improves the axial resolution and contrast in the reconstructed image. The restored PA signals can be considered as the input of the conventional reconstruction algorithms (e.g., Delay-and-sum (DAS) and Delay-multiply-and-sum (DMAS)). To compare the performance of the proposed method, DAS and DMAS reconstruction algorithms were performed with both the initial and restored PA signals on numerical and experimental studies (numerical targets, tungsten wires, and human forearm). The results show that, compared with the initial PA signals, the restored PA signals can improve the axial resolution and contrast by 45% and 16.1 dB, respectively, and suppress background artifacts by 80%.

Improving Minimum Variance Beamforming with Sub-Aperture Processing for Photoacoustic Imaging

Minimum variance (MV) beamforming improves resolution and reduces sidelobes when compared to delay-and-sum (DAS) beamforming for photoacoustic imaging (PAI). However, some level of sidelobe signal and incoherent clutter persist degrading MV PAI quality. Here, an adaptive beamforming algorithm (PSAPMV) combining MV formulation and sub-aperture processing is proposed. In PSAPMV, the received channel data are split into two complementary nonoverlapping sub-apertures and beamformed using MV. A weighting matrix based on similarity between sub-aperture beamformed images was derived and multiplied with the full aperture MV image resulting in suppression of sidelobe and incoherent clutter in the PA image. Numerical simulation experiments with point targets, diffuse inclusions and microvasculature networks are used to validate PSAPMV. Quantitative evaluation was done in terms of main-lobe-to-side-lobe ratio, full width at half maximum (FWHM), contrast ratio (CR) and generalized contrast-to-noise ratio (gCNR). PSAPMV demonstrated improved beamforming performance both qualitatively and quantitatively. PSAPMV had higher resolution (FWHM =0.19 mm) than MV (0.21 mm) and DAS (0.22mm) in point target simulations, better target detectability (gCNR =0.99) than MV (0.89) and DAS (0.84) for diffuse inclusions and improved contrast (CR in microvasculature simulation, DAS = 15.38, MV = 22.42, PSAPMV = 51.74 dB).

LightSpeed: A Compact, High-Speed Optical-Link-Based 3D Optoacoustic Imager

Wide-scale adoption of optoacoustic imaging in biology and medicine critically depends on availability of affordable scanners combining ease of operation with optimal imaging performance. Here we introduce LightSpeed: a low-cost real-time volumetric handheld optoacoustic imager based on a new compact software-defined ultrasound digital acquisition platform and a pulsed laser diode. It supports the simultaneous signal acquisition from up to 192 ultrasound channels and provides a hig-bandwidth direct optical link (2x 100G Ethernet) to the host-PC for ultra-high frame rate image acquisitions. We demonstrate use of the system for ultrafast (500Hz) 3D human angiography with a rapidly moving handheld probe. LightSpeed attained image quality comparable with a conventional optoacoustic imaging systems employing bulky acquisition electronics and a Q-switched pulsed laser. Our results thus pave the way towards a new generation of compact, affordable and high-performance optoacoustic scanners.

Subaperture Processing-Based Adaptive Beamforming for Photoacoustic Imaging

Delay-and-sum (DAS) beamformers, when applied to photoacoustic (PA) image reconstruction, produce strong sidelobes due to the absence of transmit focusing. Consequently, DAS PA images are often severely degraded by strong off-axis clutter. For preclinical in vivo cardiac PA imaging, the presence of these noise artifacts hampers the detectability and interpretation of PA signals from the myocardial wall, crucial for studying blood-dominated cardiac pathological information and to complement functional information derived from ultrasound imaging. In this article, we present PA subaperture processing (PSAP), an adaptive beamforming method, to mitigate these image degrading effects. In PSAP, a pair of DAS reconstructed images is formed by splitting the received channel data into two complementary nonoverlapping subapertures. Then, a weighting matrix is derived by analyzing the correlation between subaperture beamformed images and multiplied with the full-aperture DAS PA image to reduce sidelobes and incoherent clutter. We validated PSAP using numerical simulation studies using point target, diffuse inclusion and microvasculature imaging, and in vivo feasibility studies on five healthy murine models. Qualitative and quantitative analysis demonstrate improvements in PAI image quality with PSAP compared to DAS and coherence factor weighted DAS (DAS _CF ). PSAP demonstrated improved target detectability with a higher generalized contrast-to-noise (gCNR) ratio in vasculature simulations where PSAP produces 19.61% and 19.53% higher gCNRs than DAS and DAS _CF , respectively. Furthermore, PSAP provided higher image contrast quantified using contrast ratio (CR) (e.g., PSAP produces 89.26% and 11.90% higher CR than DAS and DAS _CF in vasculature simulations) and improved clutter suppression.

Enhancement of in vivo cardiac photoacoustic signal specificity using spatiotemporal singular value decomposition

SIGNIFICANCE

Photoacoustic imaging (PAI) can be used to infer molecular information about myocardial health non-invasively in vivo using optical excitation at ultrasonic spatial resolution. For clinical and preclinical linear array imaging systems, conventional delay-and-sum (DAS) beamforming is typically used. However, DAS cardiac PA images are prone to artifacts such as diffuse quasi-static clutter with temporally varying noise-reducing myocardial signal specificity. Typically, multiple frame averaging schemes are utilized to improve the quality of cardiac PAI, which affects the spatial and temporal resolution and reduces sensitivity to subtle PA signal variation. Furthermore, frame averaging might corrupt myocardial oxygen saturation quantification due to the presence of natural cardiac wall motion. In this paper, a spatiotemporal singular value decomposition (SVD) processing algorithm is proposed to reduce DAS PAI artifacts and subsequent enhancement of myocardial signal specificity.

AIM

Demonstrate enhancement of PA signals from myocardial tissue compared to surrounding tissues and blood inside the left-ventricular (LV) chamber using spatiotemporal SVD processing with electrocardiogram (ECG) and respiratory signal (ECG-R) gated in vivo murine cardiac PAI.

APPROACH

In vivo murine cardiac PAI was performed by collecting single wavelength (850 nm) photoacoustic channel data on eight healthy mice. A three-dimensional (3D) volume of complex PAI data over a cardiac cycle was reconstructed using a custom ECG-R gating algorithm and DAS beamforming. Spatiotemporal SVD was applied on a two-dimensional Casorati matrix generated using the 3D volume of PAI data. The singular value spectrum (SVS) was then filtered to remove contributions from diffuse quasi-static clutter and random noise. Finally, SVD processed beamformed images were derived using filtered SVS and inverse SVD computations.

RESULTS

Qualitative comparison with DAS and minimum variance (MV) beamforming shows that SVD processed images had better myocardial signal specificity, contrast, and target detectability. DAS, MV, and SVD images were quantitatively evaluated by calculating contrast ratio (CR), generalized contrast-to-noise ratio (gCNR), and signal-to-noise ratio (SNR). Quantitative evaluations were done at three cardiac time points (during systole, at end-systole (ES), and during diastole) identified from co-registered ultrasound M-Mode image. Mean CR, gCNR, and SNR values of SVD images at ES were 245, 115.15, and 258.17 times higher than DAS images with statistical significance evaluated with one-way analysis of variance.

CONCLUSIONS

Our results suggest that significantly better-quality images can be realized using spatiotemporal SVD processing for in vivo murine cardiac PAI.

Photoacoustic tomography reconstruction using lag-based delay multiply and sum with a coherence factor improves in vivo ovarian cancer diagnosis

Ovarian cancer is the fifth most common cause of death due to cancer, and it is the deadliest of all gynecological cancers. Diagnosing ovarian cancer via conventional photoacoustic delay-and-sum beamforming (DAS) presents several challenges, such as poor image resolution and low lesion to background tissue contrast. To address these concerns, we propose an improved beamformer named lag-based delay multiply and sum combined with coherence factor (DMAS-LAG-CF). Simulations and phantom experiments demonstrate that compared with the conventional DAS, the proposed algorithm can provide 1.39 times better resolution and 10.77 dB higher contrast. For patient data, similar performance on contrast ratios has been observed. However, since the diagnostic accuracy between cancer and benign/normal groups is a significant measure, we have extracted photoacoustic histogram features of mean, kurtosis and skewness. DMAS-LAG-CF can improve cancer diagnosis with an AUC of 0.91 for distinguishing malignant vs. benign ovarian lesions when mean and skewness are used as features.