Image improvement in linear-array photoacoustic imaging using high resolution coherence factor weighting technique

Evaluation of 10 current image reconstruction algorithms for linear array photoacoustic imaging

Various reconstruction algorithms have been implemented for linear array photoacoustic imaging systems with the goal of accurately reconstructing the strength absorbers within the tissue being imaged. Since the existing algorithms have been introduced by different research groups and the context of performance evaluation was not consistent, it is difficult to make a fair comparison between them. In this study, we systematically compared the performance of 10 published image reconstruction algorithms (DAS, UBP, pDAS, DMAS, MV, EIGMV, SLSC, GSC, TR, and FD) using in-vitro phantom data. Evaluations were conducted based on lateral resolution of the reconstructed images, computational time, target detectability, and noise sensitivity. We anticipate the outcome of this study will assist researchers in selecting appropriate algorithms for their linear array PA imaging applications.

Restoring the imaging quality of circular transducer array-based PACT using synthetic aperture focusing technique integrated with 2nd-derivative-based back projection scheme

Circular-array-based photoacoustic computed tomography (CA-PACT) is a promising imaging tool owing to its broad acoustic detection coverage and fidelity. However, CA-PACT suffers from poor image quality outside the focal zone along both elevational and lateral dimensions. To address this challenge, we proposed a novel reconstruction strategy by integrating the synthetic aperture focusing technique (SAFT) with the 2nd derivative-based back projection (2nd D-BP) algorithm to restore the image quality outside the focal zone along both the elevational and lateral axes. The proposed solution is a two-phase reconstruction scheme. In the first phase, with the assistance of an acoustic lens, we designed a circular array-based SAFT algorithm to restore the resolution and SNR along the elevational axis. The acoustic lens pushes the boundary of the upper limit of the SAFT scheme to achieve enhanced elevational resolution. In the second phase, we proposed a 2nd D-BP scheme to improve the lateral resolution and suppress noises in 3D imaging results. The 2nd D-BP strategy enhances the image quality along the lateral dimension by up-converting the high spatial frequencies of the object's absorption pattern. We validated the effectiveness of the proposed strategy using both phantoms and in vivo human experiments. The experimental results demonstrated superior image quality (7-fold enhancement in elevational resolution, 3-fold enhancement in lateral resolution, and an 11-dB increase in SNR). This strategy provides a new paradigm in the PACT system as it significantly enhances the spatial resolution and imaging contrast in both the elevational and lateral dimensions while maintaining a large focal zone.

A Robust Lamb Wave Imaging Approach to Plate-Like Structural Health Monitoring of Materials With Transducer Array Position Errors

The Lamb-wave-based damage imaging via beamforming techniques, which can visualize the location of damage in the structure intuitively, is one of the most promising methods in the field of structural health monitoring (SHM). However, transducer array position errors are inevitable in practical application, which may lead to serious degradation in imaging performance. In this study, it is shown that the uncertainty of the steering vectors led by the imprecise position of transducers in an array can be suppressed by the doubly constrained robust Capon beamformer (DCRCB). After the unwanted side lobes are restrained by the DCRCB-based coherence factor (CF) weighting, an effective adaptive beamforming damage imaging method robust to transducer position errors is proposed. The numerical simulation and imaging experiment of damage on an aluminum plate are carried out to verify the effectiveness of the proposed algorithm. The results show that the proposed Lamb wave damage imaging method performs better than the reported beamforming ones in literature in terms of resolution, contrast, and robustness to transducer position errors.

Achieving depth-independent lateral resolution in AR-PAM using the synthetic-aperture focusing technique

Acoustic-resolution photoacoustic microscopy (AR-PAM) is a promising imaging modality that renders images with ultrasound resolution and extends the imaging depth beyond the optical ballistic regime. To achieve a high lateral resolution, a large numerical aperture (NA) of a focused transducer is usually applied for AR-PAM. However, AR-PAM fails to hold its performance in the out-of-focus region. The lateral resolution and signal-to-noise ratio (SNR) degrade substantially, thereby leading to a significantly deteriorated image quality outside the focal area. Based on the concept of the synthetic-aperture focusing technique (SAFT), various strategies have been developed to address this challenge. These include 1D-SAFT, 2D-SAFT, adaptive-SAFT, spatial impulse response (SIR)-based schemes, and delay-multiply-and-sum (DMAS) strategies. These techniques have shown progress in achieving depth-independent lateral resolution, while several challenges remain. This review aims to introduce these developments in SAFT-based approaches, highlight their fundamental mechanisms, underline the advantages and limitations of each approach, and discuss the outlook of the remaining challenges for future advances.

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.

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.

Spatiotemporal Coherence Weighting for In Vivo Cardiac Photoacoustic Image Beamformation

Photoacoustic (PA) image reconstruction generally utilizes delay-and-sum (DAS) beamforming of received acoustic waves from tissue irradiated with optical illumination. However, nonadaptive DAS reconstructed cardiac PA images exhibit temporally varying noise which causes reduced myocardial PA signal specificity, making image interpretation difficult. Adaptive beamforming algorithms such as minimum variance (MV) with coherence factor (CF) weighting have been previously reported to improve the DAS image quality. In this article, we report on an adaptive beamforming algorithm by extending CF weighting to the temporal domain for preclinical cardiac PA imaging (PAI). The proposed spatiotemporal coherence factor (STCF) considers multiple temporally adjacent image acquisition events during beamforming and cancels out signals with low spatial coherence and temporal coherence, resulting in higher background noise cancellation while preserving the main features of interest (myocardial wall) in the resultant PA images. STCF has been validated using the numerical simulations and in vivo ECG and respiratory-signal-gated cardiac PAI in healthy murine hearts. The numerical simulation results demonstrate that STCF weighting outperforms DAS and MV beamforming with and without CF weighting under different levels of inherent contrast, acoustic attenuation, optical scattering, and signal-to-noise (SNR) of channel data. Performance improvement is attributed to higher sidelobe reduction (at least 5 dB) and SNR improvement (at least 10 dB). Improved myocardial signal specificity and higher signal rejection in the left ventricular chamber and acoustic gel region are observed with STCF in cardiac PAI.

A dynamic generalized coherence factor for side lobe suppression in ultrasound imaging

Coherence-based weighting techniques have been widely studied to weight beamsummed data to improve image quality in ultrasound imaging. Although generalized coherence factor (GCF) enhances the robustness of coherence factor (CF) with preserved speckle pattern by including some incoherent components, the side lobe suppression performance is insufficient due to constant cut-off frequency M₀. To address this problem, we introduced in this paper a dynamic GCF method, referred to as DGCF-C, based on the amplitude standard deviation and the convolution output of aperture data. The cut-off frequency is adaptively selected for GCF at each imaging point using the amplitude standard deviation of aperture data. Moreover, the convolution output of aperture data is used to calculate the dynamic GCF. The proposed method is evaluated in simulation and tissue-mimicking phantom studies. The image quality was analyzed in terms of resolution, contrast ratio (CR), generalized contrast-to-noise ratio (GCNR), speckle signal-to-noise ratio (sSNR), and signal-to-noise ratio (SNR). The results demonstrate that DGCF-C (M_max=2) achieves mean resolution improvements of 35.1% in simulation, and 32.6% in experiment, compared with GCF (M₀=1). Moreover, DGCF-C (M_max=4) outperforms GCF (M₀=2) with an average GCNR improvement of 13.5% and an average sSNR improvement of 15.2%, which indicates the better-preservation of speckle.

Miniaturized phased-array ultrasound and photoacoustic endoscopic imaging system

Visualization and detection of early-stage gynecological malignancies represents a challenge for imaging due to limiting factors including tissue accessibility, device ease of use, and accuracy of imaging modalities. In this work, we introduce a miniaturized phased-array ultrasound and photoacoustic endoscopic probe which is capable of providing structural, functional, and molecular data for the characterization of gynecologic disease. The proposed probe consists of a 64-element ultrasound phased-array transducer coupled to a fiber-optic light delivery system for co-registered ultrasound and photoacoustic imaging. The fabricated US and PA imaging endoscope's diameter is 7.5 mm, allowing for potential passage through the cervical canal and thus an intimate contact with gynecological tissues such as the cervical canal and uterus. The developed endoscopic probe was tested and characterized in a set of tissue-mimicking phantoms. US and PA resolutions were measured experimentally using 200 μm diameter wires, resulting in apparent axial and lateral diameters of 289 μm and 299 μm for US, and 308 μm and 378 μm for PA, respectively. The probe's abilities to operate in both discrete and integrated illumination/acquisition were tested in gelatin phantoms with embedded optical absorbers with the results demonstrating the ability to acquire volumetric dual-modal US and PA images.

Efficient nonlinear beamformer based on P'th root of detected signals for linear-array photoacoustic tomography: application to sentinel lymph node imaging

In linear-array transducer-based photoacoustic (PA) imaging, B-scan PA images are formed using the raw channel PA signals. Delay-and-sum (DAS) is the most prevalent algorithm due to its simple implementation, but it leads to low-quality images. Delay-multiply-and-sum (DMAS) provides a higher image quality in comparison with DAS while it imposes a computational burden of O  (  M2  )  . We introduce a nonlinear (NL) beamformer for linear-array PA imaging, which uses the p'th root of the detected signals and imposes the complexity of DAS [O  (  M  )  ]. The proposed algorithm is evaluated numerically and experimentally [wire-target and in-vivo sentinel lymph node (SLN) imaging], and the effects of the parameter p are investigated. The results show that the NL algorithm, using a root of p (NL_p), leads to lower sidelobes and higher signal-to-noise ratio compared with DAS and DMAS, for (p  >  2). The sidelobes level (for the wire-target phantom), at the depth of 11.4 mm, are about -31, -52, -52, -67, -88, and -109  dB, for DAS, DMAS, NL_2, NL_3, NL_4, and NL_5, respectively, indicating the superiority of the NL_p algorithm. In addition, the best value of p for SLN imaging is reported to be 12.