Investigation of light delivery geometries for photoacoustic applications using Monte Carlo simulations with multiple wavelengths, tissue types, and species characteristics

Dual-Scan Photoacoustic Tomography for the Imaging of Vascular Structure on Foot

Chronic leg ulcers are affecting approximately 6.5 million Americans, and they are associated with significant mortality, reduced quality of life, and high treatment costs. Since many chronic ulcers have underlying vascular insufficiency, accurate assessment of tissue perfusion is critical to treatment planning and monitoring. This study introduces a dual-scan photoacoustic (PA) tomography (PAT) system that can simultaneously image the dorsal and plantar sides of the foot to reduce imaging time. To account for the unique shape of the foot, the system employs height-adjustable and articulating baseball stages that can scan along the foot's contour. In vivo results from healthy volunteers demonstrate the system's ability to acquire clear images of foot vasculature, and results from patients indicate that the system can image patients with various ulcer conditions. We also investigated various PA features and examined their correlation with the foot condition. Our preliminary results indicate that vessel sharpness, occupancy, intensity, and density could all be used to assess tissue perfusion. This research demonstrated the potential of PAT for routine clinical tissue perfusion assessment.

Whole-Body Photoacoustic Imaging Techniques for Preclinical Small Animal Studies

Photoacoustic imaging is a hybrid imaging technique that has received considerable attention in biomedical studies. In contrast to pure optical imaging techniques, photoacoustic imaging enables the visualization of optical absorption properties at deeper imaging depths. In preclinical small animal studies, photoacoustic imaging is widely used to visualize biodistribution at the molecular level. Monitoring the whole-body distribution of chromophores in small animals is a key method used in preclinical research, including drug-delivery monitoring, treatment assessment, contrast-enhanced tumor imaging, and gastrointestinal tracking. In this review, photoacoustic systems for the whole-body imaging of small animals are explored and summarized. The configurations of the systems vary with the scanning methods and geometries of the ultrasound transducers. The future direction of research is also discussed with regard to achieving a deeper imaging depth and faster imaging speed, which are the main factors that an imaging system should realize to broaden its application in biomedical studies.

Enhanced Vascular Features in Porcine Gastrointestinal Endoscopy Using Multispectral Imaging

Endoscopic investigation is a predominant stan-dard while assessing the gastrointestinal tract. Even though it has been rigorously used in diagnostics for many decades, a high miss rate has been recorded. Advanced endoscopic imaging still has not found solutions to problems like early cancer detection, polyp generality, disease classification, etc. One of the less explored techniques to study early cancer detection is spectral imaging which deals with the absorption and reflection spectra of various wavelengths of light by different layers of tissue. To study tissues under various illumination, a multi-spectral light source unit that can be used along with an endoscopy system was developed with 10 different LEDs of very narrow bandwidths. Using this light source, a feasibility study was per-formed on an animal in which the upper GI tract of a porcine model was imaged and sample images were taken for processing from five different sections. Some wavelengths showed better contrast enhancements for visualization of vascular structures. Wavelength 420 nm (violet light) showed better contrast and the gradient of the line profile histogram showed the highest intensity change between the blood vessels and the surrounding mucosa. These enhancements showed that spectral imaging can potentially help in studying tissues for early cancer detection and improved visualization of the G I tract using endoscopy.

Randomized multi-angle illumination for improved linear array photoacoustic computed tomography in brain

One of the key challenges in linear array transducer-based photoacoustic computed tomography is to image structures embedded deep within the biological tissue with limited optical energy. Here, we utilized a manually controlled multi-angle illumination technique to allow the incident photons to interact with the imaging targets for longer periods of time and diffuse further in all directions. We have developed and optimized a compact probe that enables manual changes to the angle of illumination while acquiring photoacoustic signals. The performance has been demonstrated and evaluated by imaging complex blood vessel mimicking phantoms in-vitro and sheep brain samples ex-vivo. For effective image reconstruction from the data acquired by multi-angle illumination method, we have utilized a method based on the extraction of maximum intensity. In both cases, multi-angle illumination has out-performed the conventional fixed angle illumination technique to improve the overall image quality. Specifically, extraction of the imaging targets located at greater axial depths was possible using this multi-angle illumination technique.

SIMPA: an open-source toolkit for simulation and image processing for photonics and acoustics

SIGNIFICANCE

Optical and acoustic imaging techniques enable noninvasive visualisation of structural and functional properties of tissue. The quantification of measurements, however, remains challenging due to the inverse problems that must be solved. Emerging data-driven approaches are promising, but they rely heavily on the presence of high-quality simulations across a range of wavelengths due to the lack of ground truth knowledge of tissue acoustical and optical properties in realistic settings.

AIM

To facilitate this process, we present the open-source simulation and image processing for photonics and acoustics (SIMPA) Python toolkit. SIMPA is being developed according to modern software design standards.

APPROACH

SIMPA enables the use of computational forward models, data processing algorithms, and digital device twins to simulate realistic images within a single pipeline. SIMPA's module implementations can be seamlessly exchanged as SIMPA abstracts from the concrete implementation of each forward model and builds the simulation pipeline in a modular fashion. Furthermore, SIMPA provides comprehensive libraries of biological structures, such as vessels, as well as optical and acoustic properties and other functionalities for the generation of realistic tissue models.

RESULTS

To showcase the capabilities of SIMPA, we show examples in the context of photoacoustic imaging: the diversity of creatable tissue models, the customisability of a simulation pipeline, and the degree of realism of the simulations.

CONCLUSIONS

SIMPA is an open-source toolkit that can be used to simulate optical and acoustic imaging modalities. The code is available at: https://github.com/IMSY-DKFZ/simpa, and all of the examples and experiments in this paper can be reproduced using the code available at: https://github.com/IMSY-DKFZ/simpa_paper_experiments.

Assessment of photoacoustic tomography contrast for breast tissue imaging using 3D correlative virtual histology

Current breast tumor margin detection methods are destructive, time-consuming, and result in significant reoperative rates. Dual-modality photoacoustic tomography (PAT) and ultrasound has the potential to enhance breast margin characterization by providing clinically relevant compositional information with high sensitivity and tissue penetration. However, quantitative methods that rigorously compare volumetric PAT and ultrasound images with gold-standard histology are lacking, thus limiting clinical validation and translation. Here, we present a quantitative multimodality workflow that uses inverted Selective Plane Illumination Microscopy (iSPIM) to facilitate image co-registration between volumetric PAT-ultrasound datasets with histology in human invasive ductal carcinoma breast tissue samples. Our ultrasound-PAT system consisted of a tunable Nd:YAG laser coupled with a 40 MHz central frequency ultrasound transducer. A linear stepper motor was used to acquire volumetric PAT and ultrasound breast biopsy datasets using 1100 nm light to identify hemoglobin-rich regions and 1210 nm light to identify lipid-rich regions. Our iSPIM system used 488 nm and 647 nm laser excitation combined with Eosin and DRAQ5, a cell-permeant nucleic acid binding dye, to produce high-resolution volumetric datasets comparable to histology. Image thresholding was applied to PAT and iSPIM images to extract, quantify, and topologically visualize breast biopsy lipid, stroma, hemoglobin, and nuclei distribution. Our lipid-weighted PAT and iSPIM images suggest that low lipid regions strongly correlate with malignant breast tissue. Hemoglobin-weighted PAT images, however, correlated poorly with cancerous regions determined by histology and interpreted by a board-certified pathologist. Nuclei-weighted iSPIM images revealed similar cellular content in cancerous and non-cancerous tissues, suggesting malignant cell migration from the breast ducts to the surrounding tissues. We demonstrate the utility of our nondestructive, volumetric, region-based quantitative method for comprehensive validation of 3D tomographic imaging methods suitable for bedside tumor margin detection.

Mechanical modeling and characterization of human skin: A review

This paper reviews the advances made in recent years on modeling approaches and experimental techniques to characterize the mechanical properties of human skin. The skin is the largest organ of the human body that has a complex multi-layered structure with different mechanical behaviors. The mechanical properties of human skin play an important role in distinguishing between healthy and unhealthy skin. Furthermore, knowing these mechanical properties enables computer simulation, skin research, clinical studies, as well as diagnosis and treatment monitoring of skin diseases. This paper reviews the recent efforts on modeling skin using linear, nonlinear, viscoelastic, and anisotropic materials. The work also focuses on aging effects, microstructure analysis, and non-invasive methods for skin testing. A detailed explanation of the skin structure and numerical models, such as finite element models, are discussed in this work. This work also compares different experimental methods that measure the mechanical properties of human skin. The work reviews the experimental results in the literature and shows how the mechanical properties of human skin vary with the skin sites, the layers, and the structure of human skin. The paper also discusses how state-of-the-art technology can advance skin research.

In-silico investigation towards the non-invasive optical detection of blood lactate

This paper uses Monte Carlo simulations to investigate the interaction of short-wave infrared (SWIR) light with vascular tissue as a step toward the development of a non-invasive optical sensor for measuring blood lactate in humans. The primary focus of this work was to determine the optimal source-detector separation, penetration depth of light at SWIR wavelengths in tissue, and the optimal light power required for reliable detection of lactate. The investigation also focused on determining the non-linear variations in absorbance of lactate at a few select SWIR wavelengths. SWIR photons only penetrated 1.3 mm and did not travel beyond the hypodermal fat layer. The maximum output power was only 2.51% of the input power, demonstrating the need for a highly sensitive detection system. Simulations optimized a source-detector separation of 1 mm at 1684 nm for accurate measurement of lactate in blood.

Towards Transabdominal Functional Photoacoustic Imaging of the Placenta: Improvement in Imaging Depth Through Optimization of Light Delivery

Functional photoacoustic imaging of the placenta could provide an innovative tool to diagnose preeclampsia, monitor fetal growth restriction, and determine the developmental impacts of gestational diabetes. However, transabdominal photoacoustic imaging is limited in imaging depth due to the tissue's scattering and absorption of light. The aim of this paper was to investigate the impact of geometry and wavelength on transabdominal light delivery. Our methods included the development of a multilayer model of the abdominal tissue and simulation of the light propagation using Monte Carlo methods. A bifurcated light source with varying incident angle of light, distance between light beams, and beam area was simulated to analyze the effect of light delivery geometry on the fluence distribution at depth. The impact of wavelength and the effects of variable thicknesses of adipose tissue and muscle were also studied. Our results showed that the beam area plays a major role in improving the delivery of light to deep tissue, in comparison to light incidence angle or distance between the bifurcated fibers. Longer wavelengths, with incident fluence at the maximum permissible exposure limit, also increases fluence within deeper tissue. We validated our simulations using a commercially available light delivery system and ex vivo human placental tissue. Additionally, we compared our optimized light delivery to a commercially available light delivery system, and conclude that our optimized geometry could improve imaging depth more than 1.6×, bringing the imaging depth to within the needed range for transabdominal imaging of the human placenta.

Multiangle Long-Axis Lateral Illumination Photoacoustic Imaging Using Linear Array Transducer

Photoacoustic imaging (PAI) combines optical contrast with ultrasound spatial resolution and can be obtained up to a depth of a few centimeters. Hand-held PAI systems using linear array usually operate in reflection mode using a dark-field illumination scheme, where the optical fiber output is attached to both sides of the elevation plane (short-axis) of the transducer. More recently, bright-field strategies where the optical illumination is coaxial with acoustic detection have been proposed to overcome some limitations of the standard dark-field approach. In this paper, a novel multiangle long-axis lateral illumination is proposed. Monte Carlo simulations were conducted to evaluate light delivery for three different illumination schemes: bright-field, standard dark-field, and long-axis lateral illumination. Long-axis lateral illumination showed remarkable improvement in light delivery for targets with a width smaller than the transducer lateral dimension. A prototype was developed to experimentally demonstrate the feasibility of the proposed approach. In this device, the fiber bundle terminal ends are attached to both sides of the transducer’s long-axis and the illumination angle of each fiber bundle can be independently controlled. The final PA image is obtained by the coherent sum of subframes acquired using different angles. The prototype was experimentally evaluated by taking images from a phantom, a mouse abdomen, forearm, and index finger of a volunteer. The system provided light delivery enhancement taking advantage of the geometry of the target, achieving sufficient signal-to-noise ratio at clinically relevant depths.

Optimizing Irradiation Geometry in LED-Based Photoacoustic Imaging with 3D Printed Flexible and Modular Light Delivery System

Photoacoustic (PA) imaging–a technique combining the ability of optical imaging to probe functional properties of the tissue and deep structural imaging ability of ultrasound–has gained significant popularity in the past two decades for its utility in several biomedical applications. More recently, light-emitting diodes (LED) are being explored as an alternative to bulky and expensive laser systems used in PA imaging for their portability and low-cost. Due to the large beam divergence of LEDs compared to traditional laser beams, it is imperative to quantify the angular dependence of LED-based illumination and optimize its performance for imaging superficial or deep-seated lesions. A custom-built modular 3-D printed hinge system and tissue-mimicking phantoms with various absorption and scattering properties were used in this study to quantify the angular dependence of LED-based illumination. We also experimentally calculated the source divergence of the pulsed-LED arrays to be 58° ± 8°. Our results from point sources (pencil lead phantom) in non-scattering medium obey the cotangential relationship between the angle of irradiation and maximum PA intensity obtained at various imaging depths, as expected. Strong dependence on the angle of illumination at superficial depths (−5°/mm at 10 mm) was observed that becomes weaker at intermediate depths (−2.5°/mm at 20 mm) and negligible at deeper locations (−1.1°/mm at 30 mm). The results from the tissue-mimicking phantom in scattering media indicate that angles between 30–75° could be used for imaging lesions at various depths (12 mm–28 mm) where lower LED illumination angles (closer to being parallel to the imaging plane) are preferable for deep tissue imaging and superficial lesion imaging is possible with higher LED illumination angles (closer to being perpendicular to the imaging plane). Our results can serve as a priori knowledge for the future LED-based PA system designs employed for both preclinical and clinical applications.

3D Monte Carlo simulation of light distribution in mouse brain in quantitative photoacoustic computed tomography

Background

Photoacoustic computed tomography (PACT) detects light-induced ultrasound (US) waves to reconstruct the optical absorption contrast of the biological tissues. Due to its relatively deep penetration (several centimeters in soft tissue), high spatial resolution, and inherent functional sensitivity, PACT has great potential for imaging mouse brains with endogenous and exogenous contrasts, which is of immense interest to the neuroscience community. However, conventional PACT either assumes homogenous optical fluence within the brain or uses a simplified attenuation model for optical fluence estimation. Both approaches underestimate the complexity of the fluence heterogeneity and can result in poor quantitative imaging accuracy.

Methods

To optimize the quantitative performance of PACT, we explore for the first time 3D Monte Carlo (MC) simulation to study the optical fluence distribution in a complete mouse brain model. We apply the MCX MC simulation package on a digital mouse (Digimouse) brain atlas that has complete anatomy information. To evaluate the impact of the brain vasculature on light delivery, we also incorporate the whole-brain vasculature in the Digimouse atlas. k-wave toolbox was used to investigate the effect of inhomogeneous illumination on the reconstructed images and chromophore concentration estimation.

Results

The simulation results clearly show that the optical fluence in the mouse brain is heterogeneous at the global level and can decrease by a factor of five with increasing depth. Moreover, the strong absorption and scattering of the brain vasculature also induce the fluence disturbance at the local level.

Conclusions

Both global and local fluence heterogeneity contributes to the reduced quantitative accuracy of the reconstructed PACT images of mouse brain. Correcting the optical fluence distribution can improve the quantitative accuracy of PACT.