Photoacoustic imaging (PAI) is a powerful imaging modality that relies on the PA effect. PAI works on the principle of electromagnetic energy absorption by the exogenous contrast agents and/or endogenous molecules present in the biological tissue, consequently generating ultrasound waves. PAI combines a high optical contrast with a high acoustic spatiotemporal resolution, allowing the non-invasive visualization of absorbers in deep structures. However, due to the optical diffusion and ultrasound attenuation in heterogeneous turbid biological tissue, the quality of the PA images deteriorates. Therefore, signal and image-processing techniques are imperative in PAI to provide high-quality images with detailed structural and functional information in deep tissues. Here, we review various signal and image processing techniques that have been developed/implemented in PAI. Our goal is to highlight the importance of image computing in photoacoustic imaging.

One of the major concerns in photoacoustic computed tomography (PACT) is obtaining a high-quality image using the minimum number of ultrasound transducers/view angles. This issue is of importance when a cost-effective PACT system is needed. On the other hand, analytical reconstruction algorithms such as back projection (BP) and time reversal, when a limited number of view angles is used, cause artifacts in the reconstructed image. Iterative algorithms provide a higher image quality, compared to BP, due to a model used for image reconstruction. The performance of the model can be further improved using the sparsity concept. In this paper, we propose using a novel sparse dictionary to capture important features of the photoacoustic signal and eliminate the artifacts while few transducers is used. Our dictionary is an optimum combination of Wavelet Transform (WT), Discrete Cosine Transform (DCT), and Total Variation (TV). We utilize two quality assessment metrics including peak signal-to-noise ratio and edge preservation index to quantitatively evaluate the reconstructed images. The results show that the proposed method can generate high-quality images having fewer artifacts and preserved edges, when fewer view angles are used for reconstruction in PACT.

A low-cost Photoacoustic Computed Tomography (PACT) system consisting of 16 single-element transducers has been developed. Our design proposes a fast rotating mechanism of 360° rotation around the imaging target, generating comparable images to those produced by large-number-element (e.g., 512, 1024, etc.) ring-array PACT systems. The 2D images with a temporal resolution of 1.5 s and a spatial resolution of 240 µm were achieved. The performance of the proposed system was evaluated by imaging complex phantom. The purpose of the proposed development is to provide researchers a low-cost alternative 2D photoacoustic computed tomography system with comparable resolution to the current high performance expensive ring-array PACT systems.

Photoacoustic imaging (PAI) is an emerging label-free and non-invasive modality for imaging biological tissues. PAI has been implemented in different configurations, one of which is photoacoustic computed tomography (PACT) with a potential wide range of applications, including brain and breast imaging. Hemispherical Array PACT (HA-PACT) is a variation of PACT that has solved the limited detection-view problem. Here, we designed an HA-PACT system consisting of 50 single element transducers. For implementation, we initially performed a simulation study, with parameters close to those in practice, to determine the relationship between the number of transducers and the quality of the reconstructed image. We then used the greatest number of transducers possible on the hemisphere and imaged copper wire phantoms coated with a light absorbing material to evaluate the performance of the system. Several practical issues such as light illumination, arrangement of the transducers, and an image reconstruction algorithm have been comprehensively studied.

Galvo scanners are popular devices for fast transversal scanning. A triangular signal is usually employed to drive galvo scanners at scanning rates close to the inverse of their response time where scanning deflection becomes a nonlinear function of applied voltage. To address this, the triangular signal is synthesized from several short ramps with different slopes. An optimization algorithm similar to a simulated annealing algorithm is used for finding the optimal signal shape to drive the galvo scanners. As a result, a significant reduction in the nonlinearity of the galvo scanning is obtained.