Overview of diffuse optical tomography and its clinical applications

Deep orthogonal multi-wavelength fusion for tomogram-free diagnosis in diffuse optical imaging

Novel portable diffuse optical tomography (DOT) devices for breast cancer lesions hold great promise for non-invasive, non-ionizing breast cancer screening. Critical to this capability is not just the identification of lesions but rather the complex problem of discriminating between malignant and benign lesions. To accurately reconstruct the highly heterogeneous tissue of a cancer lesion in healthy breast tissue using DOT, multiple wavelengths can be leveraged to maximize signal penetration while minimizing sensitivity to noise. However, these wavelength responses can overlap, capture common information, and correlate, potentially confounding reconstruction and downstream end tasks. We show that an orthogonal fusion loss regularizes multi-wavelength DOT leading to improved reconstruction and accuracy of end-to-end discrimination of malignant versus benign lesions. We further show that our raw-to-task model significantly reduces computational complexity without sacrificing accuracy, making it ideal for real-time throughput, desired in medical settings where handheld devices have severely restricted power budgets. Furthermore, our results indicate that image reconstruction is not necessary for unbiased classification of lesions with a balanced accuracy of 77% and 66% on the synthetic dataset and clinical dataset, respectively, using the raw-to-task model. Code is available at https://github.com/sfu-mial/FuseNet.

Functional ultrasound imaging of human brain activity through an acoustically transparent cranial window

Visualization of human brain activity is crucial for understanding normal and aberrant brain function. Currently available neural activity recording methods are highly invasive, have low sensitivity, and cannot be conducted outside of an operating room. Functional ultrasound imaging (fUSI) is an emerging technique that offers sensitive, large-scale, high-resolution neural imaging; however, fUSI cannot be performed through the adult human skull. Here, we used a polymeric skull replacement material to create an acoustic window compatible with fUSI to monitor adult human brain activity in a single individual. Using an in vitro cerebrovascular phantom to mimic brain vasculature and an in vivo rodent cranial defect model, first, we evaluated the fUSI signal intensity and signal-to-noise ratio through polymethyl methacrylate (PMMA) cranial implants of different thicknesses or a titanium mesh implant. We found that rat brain neural activity could be recorded with high sensitivity through a PMMA implant using a dedicated fUSI pulse sequence. We then designed a custom ultrasound-transparent cranial window implant for an adult patient undergoing reconstructive skull surgery after traumatic brain injury. We showed that fUSI could record brain activity in an awake human outside of the operating room. In a video game "connect the dots" task, we demonstrated mapping and decoding of task-modulated cortical activity in this individual. In a guitar-strumming task, we mapped additional task-specific cortical responses. Our proof-of-principle study shows that fUSI can be used as a high-resolution (200 μm) functional imaging modality for measuring adult human brain activity through an acoustically transparent cranial window.

Enhanced model iteration algorithm with graph neural network for diffuse optical tomography

Diffuse optical tomography (DOT) employs near-infrared light to reveal the optical parameters of biological tissues. Due to the strong scattering of photons in tissues and the limited surface measurements, DOT reconstruction is severely ill-posed. The Levenberg-Marquardt (LM) is a popular iteration method for DOT, however, it is computationally expensive and its reconstruction accuracy needs improvement. In this study, we propose a neural model based iteration algorithm which combines the graph neural network with Levenberg-Marquardt (GNNLM), which utilizes a graph data structure to represent the finite element mesh. In order to verify the performance of the graph neural network, two GNN variants, namely graph convolutional neural network (GCN) and graph attention neural network (GAT) were employed in the experiments. The results showed that GCNLM performs best in the simulation experiments within the training data distribution. However, GATLM exhibits superior performance in the simulation experiments outside the training data distribution and real experiments with breast-like phantoms. It demonstrated that the GATLM trained with simulation data can generalize well to situations outside the training data distribution without transfer training. This offers the possibility to provide more accurate absorption coefficient distributions in clinical practice.

Agarose Gel Characterization for the Fabrication of Brain Tissue Phantoms for Infrared Multispectral Vision Systems

Synthetic phantoms that recreate the characteristics of biological tissues are valuable tools for systematically studying and comprehending physiologies, pathologies, and biological processes related to tissues. The reproduction of mechanical and optical properties allows for the development and evaluation of novel systems and applications in areas such as imaging, optics, ultrasound, or dosimetry, among others. This paper proposes a methodology for manufacturing agarose-based phantoms that mimics the optical properties of healthy brain tissue within the wavelength infrared range of 800 to 820 nm. The fabrication of such phantoms enables the possibility of testing and experimentation in controlled and safe environments toward the design of new near-infrared multispectral imaging systems in neurosurgery. The results of an experimental optical characterization study indicate the validity and reliability of the proposed method for fabricating brain tissue phantoms in a cost-effective and straightforward fashion.

Imaging inside highly scattering media using hybrid deep learning and analytical algorithm

Imaging through highly scattering media is a challenging problem with numerous applications in biomedical and remote-sensing fields. Existing methods that use analytical or deep learning tools are limited by simplified forward models or a requirement for prior physical knowledge, resulting in blurry images or a need for large training databases. To address these limitations, we propose a hybrid scheme called Hybrid-DOT that combines analytically derived image estimates with a deep learning network. Our analysis demonstrates that Hybrid-DOT outperforms a state-of-the-art ToF-DOT algorithm by improving the PSNR ratio by 4.6 dB and reducing the resolution by a factor of 2.5. Furthermore, when compared to a deep learning stand-alone model, Hybrid-DOT achieves a 0.8 dB increase in PSNR, 1.5 times the resolution, and a significant reduction in the required dataset size (factor of 1.6-3). The proposed model remains effective at higher depths, providing similar improvements for up to 160 mean-free paths.

Depth-sensitive diffuse speckle contrast topography for high-density mapping of cerebral blood flow in rodents

Significance

Frequent assessment of cerebral blood flow (CBF) is crucial for the diagnosis and management of cerebral vascular diseases. In contrast to large and expensive imaging modalities, such as nuclear medicine and magnetic resonance imaging, optical imaging techniques are portable and inexpensive tools for continuous measurements of cerebral hemodynamics. The recent development of an innovative noncontact speckle contrast diffuse correlation tomography (scDCT) enables three-dimensional (3D) imaging of CBF distributions. However, scDCT requires complex and time-consuming 3D reconstruction, which limits its ability to achieve high spatial resolution without sacrificing temporal resolution and computational efficiency.

Aim

We investigate a new diffuse speckle contrast topography (DSCT) method with parallel computation for analyzing scDCT data to achieve fast and high-density two-dimensional (2D) mapping of CBF distributions at different depths without the need for 3D reconstruction.

Approach

A new moving window method was adapted to improve the sampling rate of DSCT. A fast computation method utilizing MATLAB functions in the Image Processing Toolbox™ and Parallel Computing Toolbox™ was developed to rapidly generate high-density CBF maps. The new DSCT method was tested for spatial resolution and depth sensitivity in head-simulating layered phantoms and in-vivo rodent models.

Results

DSCT enables 2D mapping of the particle flow in the phantom at different depths through the top layer with varied thicknesses. Both DSCT and scDCT enable the detection of global and regional CBF changes in deep brains of adult rats. However, DSCT achieves fast and high-density 2D mapping of CBF distributions at different depths without the need for complex and time-consuming 3D reconstruction.

Conclusions

The depth-sensitive DSCT method has the potential to be used as a noninvasive, noncontact, fast, high resolution, portable, and inexpensive brain imager for basic neuroscience research in small animal models and for translational studies in human neonates.

Efficient algorithm to calculate the optical properties of breast tumors by high-order perturbation theory

An efficient algorithm to obtain the solutions for n-th order terms of perturbation expansions in absorption, scattering, and cross-coupling for light propagating in human tissue is presented. The proposed solution is free of any approximations and makes possible fast and efficient estimates of mammographic, optical tomographic, and fluorescent images, applying a perturbation order of 30 and more. The presented analysis sets the general limits for the applicability of the perturbation approach as a function of tumor size and optical properties of the human tissue. The convergence tests of the efficient calculations for large absorbing objects show excellent agreement with the reference data from finite element method calculations. The applicability of the theory is demonstrated in experiments on breast-like phantoms with high absorbing and low-scattering lesions.

One-dimensional convolutional neural network for Jacobian in Diffuse Optical Tomography

Non-linear least square minimization algorithms are often employed to solve Diffuse Optical Tomography (DOT) inverse problem. However, it is time-consuming to calculate the Jacobian matrix. This work has proposed a data-driven neural network method to improve computational efficiency. The singular value decomposition is employed to compute the updated Jacobian and a mapping from boundary measurements to the singular values based on a convolutional neural network (CNN) is learned to obtain the singular values. The method is validated with 3D numerical simulation data. We have demonstrated that the approach can save computation time compared to Adjoint method, and reconstructed absorption coefficient close to Adjoint method.Clinical Relevance- These results are not focused on clinical relevance currently, but in the future may be helpful to accelerant DOT reconstruction in clinic.

An Optical Sensory System for Assessment of Residual Cancer Burden in Breast Cancer Patients Undergoing Neoadjuvant Chemotherapy

Breast cancer patients undergoing neoadjuvant chemotherapy (NAC) require precise and accurate evaluation of treatment response. Residual cancer burden (RCB) is a prognostic tool widely used to estimate survival outcomes in breast cancer. In this study, we introduced a machine-learning-based optical biosensor called the Opti-scan probe to assess residual cancer burden in breast cancer patients undergoing NAC. The Opti-scan probe data were acquired from 15 patients (mean age: 61.8 years) before and after each cycle of NAC. Using regression analysis with k-fold cross-validation, we calculated the optical properties of healthy and unhealthy breast tissues. The ML predictive model was trained on the optical parameter values and breast cancer imaging features obtained from the Opti-scan probe data to calculate RCB values. The results show that the ML model achieved a high accuracy of 0.98 in predicting RCB number/class based on the changes in optical properties measured by the Opti-scan probe. These findings suggest that our ML-based Opti-scan probe has considerable potential as a valuable tool for the assessment of breast cancer response after NAC and to guide treatment decisions. Therefore, it could be a promising, non-invasive, and accurate method for monitoring breast cancer patient's response to NAC.

Functional photoacoustic imaging: from nano- and micro- to macro-scale

Functional photoacoustic imaging is a promising biological imaging technique that offers such unique benefits as scalable resolution and imaging depth, as well as the ability to provide functional information. At nanoscale, photoacoustic imaging has provided super-resolution images of the surface light absorption characteristics of materials and of single organelles in cells. At the microscopic and macroscopic scales. photoacoustic imaging techniques have precisely measured and quantified various physiological parameters, such as oxygen saturation, vessel morphology, blood flow, and the metabolic rate of oxygen, in both human and animal subjects. This comprehensive review provides an overview of functional photoacoustic imaging across multiple scales, from nano to macro, and highlights recent advances in technology developments and applications. Finally, the review surveys the future prospects of functional photoacoustic imaging in the biomedical field.

Artificial Intelligence-Aided Optical Imaging for Cancer Theranostics

The use of artificial intelligence (AI) to assist biomedical imaging have demonstrated its high accuracy and high efficiency in medical decision-making for individualized cancer medicine. In particular, optical imaging methods are able to visualize both the structural and functional information of tumors tissues with high contrast, low cost, and noninvasive property. However, no systematic work has been performed to inspect the recent advances on AI-aided optical imaging for cancer theranostics. In this review, we demonstrated how AI can guide optical imaging methods to improve the accuracy on tumor detection, automated analysis and prediction of its histopathological section, its monitoring during treatment, and its prognosis by using computer vision, deep learning and natural language processing. By contrast, the optical imaging techniques involved mainly consisted of various tomography and microscopy imaging methods such as optical endoscopy imaging, optical coherence tomography, photoacoustic imaging, diffuse optical tomography, optical microscopy imaging, Raman imaging, and fluorescent imaging. Meanwhile, existing problems, possible challenges and future prospects for AI-aided optical imaging protocol for cancer theranostics were also discussed. It is expected that the present work can open a new avenue for precision oncology by using AI and optical imaging tools.

Testicular torsion diagnosis and injury assessment using photoacoustic oxygenation imaging

Testicular torsion (TT) is a medical emergency that requires immediate diagnostic evaluation. Photoacoustic imaging (PAI) has the potential to provide spatially resolved oxygen saturation (sO₂), which can serve as a valuable marker in TT diagnosis. We investigated the potential of PAI as an alternative method for TT diagnosis and testicular injury assessment. We measured sO₂ levels in different degrees of TT models using PAI at various time points. Based on histopathological results, we found that the averaged sO₂ per pixel (sO2®) and reduction of sO2® (rsO₂) in twisted testicles had significant correlations with hypoxic conditions. Both sO2® and rsO₂ exhibited excellent diagnostic abilities in detecting TT and identifying ischemia/hypoxia injury following TT. Furthermore, PAI-measured sO₂ demonstrated favorable diagnostic capabilities in discriminating if the testicle had suffered irreversible injury. In summary, PAI presents a potentially promising novel approach in evaluating TT and warrants further clinical investigation.

Compact breast shape acquisition system for improving diffuse optical tomography image reconstructions

Diffuse optical tomography (DOT) has been investigated for diagnosing malignant breast lesions, but its accuracy relies on model-based image reconstructions, which in turn depends on the accuracy of breast shape acquisition. In this work, we have developed a dual-camera structured light imaging (SLI) breast shape acquisition system tailored for a mammography-like compression setting. Illumination pattern intensity is dynamically adjusted to account for skin tone differences, while thickness-informed pattern masking reduces artifacts due to specular reflections. This compact system is affixed to a rigid mount that can be installed into existing mammography or parallel-plate DOT systems without the need for camera-projector re-calibration. Our SLI system produces sub-millimeter resolution with a mean surface error of 0.26 mm. This breast shape acquisition system results in more accurate surface recovery, with an average 1.6-fold reduction in surface estimation errors over a reference method via contour extrusion. Such improvement translates to 25% to 50% reduction in mean squared error in the recovered absorption coefficient for a series of simulated tumors 1-2 cm below the skin.

Functional Near-Infrared Spectroscopy as a Personalized Digital Healthcare Tool for Brain Monitoring

The sustained growth of digital healthcare in the field of neurology relies on portable and cost-effective brain monitoring tools that can accurately monitor brain function in real time. Functional near-infrared spectroscopy (fNIRS) is one such tool that has become popular among researchers and clinicians as a practical alternative to functional magnetic resonance imaging, and as a complementary tool to modalities such as electroencephalography. This review covers the contribution of fNIRS to the personalized goals of digital healthcare in neurology by identifying two major trends that drive current fNIRS research. The first major trend is multimodal monitoring using fNIRS, which allows clinicians to access more data that will help them to understand the interconnection between the cerebral hemodynamics and other physiological phenomena in patients. This allows clinicians to make an overall assessment of physical health to obtain a more-detailed and individualized diagnosis. The second major trend is that fNIRS research is being conducted with naturalistic experimental paradigms that involve multisensory stimulation in familiar settings. Cerebral monitoring of multisensory stimulation during dynamic activities or within virtual reality helps to understand the complex brain activities that occur in everyday life. Finally, the scope of future fNIRS studies is discussed to facilitate more-accurate assessments of brain activation and the wider clinical acceptance of fNIRS as a medical device for digital healthcare.

Influence of Tumor Volume on the Fluence Rate Within Human Breast Model Using Continuous-Wave Diffuse Optical Imaging: A Simulation Study

Objective: This article investigates the effect of varying breast tumor size on the fluence rate distribution within a breast model during the diffuse optical imaging procedure. Background: Early detection of breast cancer is of significant importance owing to its wide spread among women worldwide. Mastectomy surgery became very common due to the late detection of breast cancers by the conventional diagnostic methods such as X-ray mammography and magnetic resonance imaging. On the contrary, optical imaging techniques provide a safe and more sensitive methodology, which is suitable for the early detection criteria. Methods: The implementation was performed based on simulating multiple detectors placed on the outer surface of a human breast model to compute the optical fluence rate after probing the breast (normal and different tumor sizes) with laser irradiation. Different laser wavelengths ranging from the red to near-infrared rays spectral range were examined to determine the optimum fluence rate that shows the highest capability to differentiate between normal and cancerous breasts. A three-dimensional breast model was created using the COMSOL multiphysics package where the optical fluence rate was estimated based on the finite-element solution of the diffusion equation. Results: To evaluate the efficiency of the suggested technique for identifying cancers and discriminate them from normal breast at various wavelengths (600-1000 nm) and several tumor sizes. Conclusions: The obtained results reveal different fluence rate distributions in the breast with different radius tumors, especially at 600 nm due to the significant differences in the scattering coefficient between malignancies and healthy tissue.

Impact of Anatomical Variability on Sensitivity Profile in fNIRS-MRI Integration

Functional near-infrared spectroscopy (fNIRS) is an important non-invasive technique used to monitor cortical activity. However, a varying sensitivity of surface channels vs. cortical structures may suggest integrating the fNIRS with the subject-specific anatomy (SSA) obtained from routine MRI. Actual processing tools permit the computation of the SSA forward problem (i.e., cortex to channel sensitivity) and next, a regularized solution of the inverse problem to map the fNIRS signals onto the cortex. The focus of this study is on the analysis of the forward problem to quantify the effect of inter-subject variability. Thirteen young adults (six males, seven females, age 29.3 ± 4.3) underwent both an MRI scan and a motor grasping task with a continuous wave fNIRS system of 102 measurement channels with optodes placed according to a 10/5 system. The fNIRS sensitivity profile was estimated using Monte Carlo simulations on each SSA and on three major atlases (i.e., Colin27, ICBM152 and FSAverage) for comparison. In each SSA, the average sensitivity curves were obtained by aligning the 102 channels and segmenting them by depth quartiles. The first quartile (depth < 11.8 (0.7) mm, median (IQR)) covered 0.391 (0.087)% of the total sensitivity profile, while the second one (depth < 13.6 (0.7) mm) covered 0.292 (0.009)%, hence indicating that about 70% of the signal was from the gyri. The sensitivity bell-shape was broad in the source-detector direction (20.953 (5.379) mm FWHM, first depth quartile) and steeper in the transversal one (6.082 (2.086) mm). The sensitivity of channels vs. different cortical areas based on SSA were analyzed finding high dispersions among subjects and large differences with atlas-based evaluations. Moreover, the inverse cortical mapping for the grasping task showed differences between SSA and atlas based solutions. In conclusion, integration with MRI SSA can significantly improve fNIRS interpretation.

Early cancer detection using the fluorescent Ashwagandha chitosan nanoparticles combined with near-infrared light diffusion characterization: in vitro study

Early cancer diagnosis through characterizing light propagation and nanotechnology increases the survival rate. The present research is aimed at evaluating the consequence of using natural nanoparticles in cancer therapy and diagnosis. Colon cancer cells were differentiated from the normal cells via investigating light diffusion combined with the fluorescence effect of the Ashwagandha chitosan nanoparticles (Ash C NPs). Ionic gelation technique synthesized the Ash C NPs. High-resolution transmission electron microscope, dynamic light scattering, and zeta potential characterized Ash C NPs. Fourier transform infrared spectroscopy analyzed Ash C NPs, chitosan, and Ashwagandha root water extract. Moreover, the MTT assay evaluated the cytotoxicity of Ash C NPs under the action of near-infrared light (NIR) irradiation. The MTT assay outcomes were statistically analyzed by Bonferroni post hoc multiple two-group comparisons using one-way variance analysis (ANOVA). Based on the Monte-Carlo simulation technique, the spatially resolved steady-state diffusely reflected light from the cancerous and healthy cells is acquired. The diffuse equation reconstructed the optical fluence rate using the finite element technique. The fluorescent effect of the nanoparticles was observed when the cells were irradiated with NIR. The MTT assay revealed a decrease in the cell viability under the action of Ash C NPs with and without laser irradiation. Colon cancer and normal cells were differentiated based on the optical characterization after laser irradiation. The light diffusion equation was successfully resolved for the fluence rate on cells' surfaces showing different normal and cancer cells values. Ash C NPs appeared its fluorescent effect in the presence of NIR laser.

Sub-mm resolution tomographic imaging in turbid media by an ultra-high density multichannel approach

We demonstrate an ultra-high-density source-detector (SD) diffuse optical tomography system scalable to thousands of combinatorial SD pairs per cm3 of total voxel volume. We demonstrate the imaging of dynamic targets (including phantom arteries) with 100 um resolution at over 10 Hz frame rate within turbid media (> 60 MFP). Further, as a step toward a wearable mobile imager, we introduce monolithic mm-size dense semiconductor laser array chips as sources for potential unobtrusive epidermal tomographic use.

Optical neuroimaging: advancing transcranial magnetic stimulation treatments of psychiatric disorders

Transcranial magnetic stimulation (TMS) has been established as an important and effective treatment for various psychiatric disorders. However, its effectiveness has likely been limited due to the dearth of neuronavigational tools for targeting purposes, unclear ideal stimulation parameters, and a lack of knowledge regarding the physiological response of the brain to TMS in each psychiatric condition. Modern optical imaging modalities, such as functional near-infrared spectroscopy and diffuse optical tomography, are promising tools for the study of TMS optimization and functional targeting in psychiatric disorders. They possess a unique combination of high spatial and temporal resolutions, portability, real-time capability, and relatively low costs. In this mini-review, we discuss the advent of optical imaging techniques and their innovative use in several psychiatric conditions including depression, panic disorder, phobias, and eating disorders. With further investment and research in the development of these optical imaging approaches, their potential will be paramount for the advancement of TMS treatment protocols in psychiatry.

Anthropomorphic Polydimethylsiloxane silicone-based phantom for Diffuse Optical Imaging

This work presents a method for constructing phantoms suitable for diffuse optical mammography. They are based on Polydimethylsiloxane silicones, with the characteristic of being anthropomorphic, and having similar mechanical and optical properties as a real breast. These phantoms are useful for testing the performance of diffuse optical imaging devices in the near infrared, both in transmittance and reflectance geometries, since they can be constructed containing inclusions, to simulate breast tumors. An alternative component to be used as scattering agent, that is easier to handle than traditional scattering agents, is also studied.

The optical properties of the phantoms were tested varying the concentration of scattering and absorbing agents, while their mechanical properties were modified by adding a silicone fluid to the basic mixture.

Finally, the phantoms were tested by Diffuse Optical Imaging experiments, and these images were compared to the ones obtained by conventional ultrasound techniques.

Results show that the constructed anthropomorphic phantoms properly reproduce the optical and mechanical characteristics of human breasts, and are suitable to be used in Diffuse Optical Imaging.

•

We constructed anthropomorphic phantoms for Diffuse Optical Imaging.

•

They simulate the optical and mechanical characteristics of a human breast.

•

A new scattering agent was successfully introduced.

•

Results of Diffuse Optical images are compared to Ultrasound images.

Noninvasive hemoglobin sensing and imaging: optical tools for disease diagnosis

SIGNIFICANCE

Measurement and imaging of hemoglobin oxygenation are used extensively in the detection and diagnosis of disease; however, the applied instruments vary widely in their depth of imaging, spatiotemporal resolution, sensitivity, accuracy, complexity, physical size, and cost. The wide variation in available instrumentation can make it challenging for end users to select the appropriate tools for their application and to understand the relative limitations of different methods.

AIM

We aim to provide a systematic overview of the field of hemoglobin imaging and sensing.

APPROACH

We reviewed the sensing and imaging methods used to analyze hemoglobin oxygenation, including pulse oximetry, spectral reflectance imaging, diffuse optical imaging, spectroscopic optical coherence tomography, photoacoustic imaging, and diffuse correlation spectroscopy.

RESULTS

We compared and contrasted the ability of different methods to determine hemoglobin biomarkers such as oxygenation while considering factors that influence their practical application.

CONCLUSIONS

We highlight key limitations in the current state-of-the-art and make suggestions for routes to advance the clinical use and interpretation of hemoglobin oxygenation information.

Optical Property Measurement and Temperature Monitoring in High-Intensity Focused Ultrasound Therapy by Diffuse Optical Tomography: A Correlation Study

In this article, we propose a new approach utilizing diffuse optical tomography (DOT) to monitoring the changes in tissues’ optical properties and temperature in high-intensity focused ultrasound (HIFU) therapy. By correlating the tissue reduced scattering coefficient (μs’) reconstructed by DOT and the temperature measured by a thermocouple, the quantitative relationship between μs’ and temperature in HIFU treatment was explored. The experiments were conducted using porcine and chicken breast muscle tissues during HIFU; the temperature of each tissue sample was recorded using a thermocouple. To incorporate the temperature dependency of tissue optical properties, both polynomial and exponential models were utilized to fit the experimental data. The results show that the change of μs’ during HIFU treatment could be detected in real-time using DOT and that this change of μs’ is quantitatively correlated with tissue temperature. Furthermore, while the tissue-type-dependent relationship between μs’ and temperature is non-linear in nature, it is stable and repeatable. Therefore, our approach has the potential to be used to predict temperature of tissue during HIFU treatment.

Different brain activation patterns in the prefrontal area between self-paced and high-speed driving tasks

The purpose of this study was to investigate the effects on prefrontal cortex brain activity when participants attempted to stop a car accurately at a stop line when driving at different speeds using functional near-infrared spectroscopy (fNIRS). Twenty healthy subjects with driving experience drove their own cars for a distance of 60 m five times each at their own pace or as fast as possible. The variation in the distance between the stop line and the car was not significantly different between the self-paced and high-speed tasks. However, oxygenated hemoglobin concentration in the prefrontal cortex was significantly higher in the high-speed task than in the self-paced task. These findings suggest that driving at high speed requires more divided attention than driving at self-paced speed, even though the participants were able to stop the car at the same distance from the target. This study shows the advantages and usefulness of fNIRS .

Challenges and advances in optical 3D mesoscale imaging

Optical mesoscale imaging is a rapidly developing field that allows the visualisation of larger samples than is possible with standard light microscopy, and fills a gap between cell and organism resolution. It spans from advanced fluorescence imaging of micrometric cell clusters to centimetre-size complete organisms. However, with larger volume specimens, new problems arise. Imaging deeper into tissues at high resolution poses challenges ranging from optical distortions to shadowing from opaque structures. This manuscript discusses the latest developments in mesoscale imaging and highlights limitations, namely labelling, clearing, absorption, scattering, and also sample handling. We then focus on approaches that seek to turn mesoscale imaging into a more quantitative technique, analogous to quantitative tomography in medical imaging, highlighting a future role for digital and physical phantoms as well as artificial intelligence.

A Model-Based Iterative Learning Approach for Diffuse Optical Tomography

Diffuse optical tomography (DOT) utilises near-infrared light for imaging spatially distributed optical parameters, typically the absorption and scattering coefficients. The image reconstruction problem of DOT is an ill-posed inverse problem, due to the non-linear light propagation in tissues and limited boundary measurements. The ill-posedness means that the image reconstruction is sensitive to measurement and modelling errors. The Bayesian approach for the inverse problem of DOT offers the possibility of incorporating prior information about the unknowns, rendering the problem less ill-posed. It also allows marginalisation of modelling errors utilising the so-called Bayesian approximation error method. A more recent trend in image reconstruction techniques is the use of deep learning, which has shown promising results in various applications from image processing to tomographic reconstructions. In this work, we study the non-linear DOT inverse problem of estimating the (absolute) absorption and scattering coefficients utilising a 'model-based' learning approach, essentially intertwining learned components with the model equations of DOT. The proposed approach was validated with 2D simulations and 3D experimental data. We demonstrated improved absorption and scattering estimates for targets with a mix of smooth and sharp image features, implying that the proposed approach could learn image features that are difficult to model using standard Gaussian priors. Furthermore, it was shown that the approach can be utilised in compensating for modelling errors due to coarse discretisation enabling computationally efficient solutions. Overall, the approach provided improved computation times compared to a standard Gauss-Newton iteration.

Haemodynamic Signatures of Temporal Integration of Visual Mirror Symmetry

EEG, fMRI and TMS studies have implicated the extra-striate cortex, including the Lateral Occipital Cortex (LOC), in the processing of visual mirror symmetries. Recent research has found that the sustained posterior negativity (SPN), a symmetry specific electrophysiological response identified in the region of the LOC, is generated when temporally displaced asymmetric components are integrated into a symmetric whole. We aim to expand on this finding using dynamic dot-patterns with systematically increased intra-pair temporal delay to map the limits of temporal integration of visual mirror symmetry. To achieve this, we used functional near-infrared spectroscopy (fNIRS) which measures the changes in the haemodynamic response to stimulation using near infrared light. We show that a symmetry specific haemodynamic response can be identified following temporal integration of otherwise meaningless dot-patterns, and the magnitude of this response scales with the duration of temporal delay. These results contribute to our understanding of when and where mirror symmetry is processed in the visual system. Furthermore, we highlight fNIRS as a promising but so far underutilised method of studying the haemodynamics of mid-level visual processes in the brain.

MCDataset: a public reference dataset of Monte Carlo simulated quantities for multilayered and voxelated tissues computed by massively parallel PyXOpto Python package

SIGNIFICANCE

Current open-source Monte Carlo (MC) method implementations for light propagation modeling are many times tedious to build and require third-party licensed software that can often discourage prospective researchers in the biomedical optics community from fully utilizing the light propagation tools. Furthermore, the same drawback also limits rigorous cross-validation of physical quantities estimated by various MC codes.

AIM

Proposal of an open-source tool for light propagation modeling and an easily accessible dataset to encourage fruitful communications amongst researchers and pave the way to a more consistent comparison between the available implementations of the MC method.

APPROACH

The PyXOpto implementation of the MC method for multilayered and voxelated tissues based on the Python programming language and PyOpenCL extension enables massively parallel computation on numerous OpenCL-enabled devices. The proposed implementation is used to compute a large dataset of reflectance, transmittance, energy deposition, and sampling volume for various source, detector, and tissue configurations.

RESULTS

The proposed PyXOpto agrees well with the original MC implementation. However, further validation reveals a noticeable bias introduced by the random number generator used in the original MC implementation.

CONCLUSIONS

Establishing a common dataset is highly important for the validation of existing and development of MC codes for light propagation in turbid media.

Speckle-resolved optical coherence tomography for mesoscopic imaging within scattering media

Light scattering poses a challenge for imaging deep in scattering media as the ballistic light exponentially attenuates with depth. In contrast to the ballistic light, the multiply scattered light penetrates deeper and also contains information about the sample. One technique to image deeper is to selectively detect only a subset of the multiply scattered light, namely the 'snake' photons, which are predominantly forward scattered and retain more direct information than the more strongly scattered light. In this work, we develop a technique, termed speckle-resolved optical coherence tomography (srOCT), for efficiently detecting these 'snake' photons to enable imaging deeper in scattering media. The system couples spatio-angular filtering with speckle-resolved interferometric detection to preferentially and efficiently detect the weakly scattered 'snake' photons. With our proof-of-concept system, we demonstrate depth-resolved imaging beyond the ballistic limit, up to a depth of 90 round-trip MFPs in a scattering phantom and a depth of 4.5 mm of chicken tissue at 0.4 mm axial and lateral resolution.

BlenderPhotonics: an integrated open-source software environment for three-dimensional meshing and photon simulations in complex tissues

SIGNIFICANCE

Rapid advances in biophotonics techniques require quantitative, model-based computational approaches to obtain functional and structural information from increasingly complex and multiscaled anatomies. The lack of efficient tools to accurately model tissue structures and subsequently perform quantitative multiphysics modeling greatly impedes the clinical translation of these modalities.

AIM

Although the mesh-based Monte Carlo (MMC) method expands our capabilities in simulating complex tissues using tetrahedral meshes, the generation of such domains often requires specialized meshing tools, such as Iso2Mesh. Creating a simplified and intuitive interface for tissue anatomical modeling and optical simulations is essential toward making these advanced modeling techniques broadly accessible to the user community.

APPROACH

We responded to the above challenge by combining the powerful, open-source three-dimensional (3D) modeling software, Blender, with state-of-the-art 3D mesh generation and MC simulation tools, utilizing the interactive graphical user interface in Blender as the front-end to allow users to create complex tissue mesh models and subsequently launch MMC light simulations.

RESULTS

Here, we present a tutorial to our Python-based Blender add-on-BlenderPhotonics-to interface with Iso2Mesh and MMC, which allows users to create, configure and refine complex simulation domains and run hardware-accelerated 3D light simulations with only a few clicks. We provide a comprehensive introduction to this tool and walk readers through five examples, ranging from simple shapes to sophisticated realistic tissue models.

CONCLUSIONS

BlenderPhotonics is user friendly and open source, and it leverages the vastly rich ecosystem of Blender. It wraps advanced modeling capabilities within an easy-to-use and interactive interface. The latest software can be downloaded at http://mcx.space/bp.

A modified source-detector configuration for the discrimination between normal and diseased human breast based on the continuous-wave diffuse optical imaging approach: a simulation study

Breast tumors are among the most common types of tumors that can affect both genders. It may spread to the whole breast without any symptoms. Therefore, the early detection and accurate diagnosis of breast tumors are significantly important. Current approaches for breast cancer screening such as positron emission tomography (PET) and magnetic resonance imaging (MRI) have some limitations of being time and money-consuming. In addition, mammography screening is not recommended for women under forty. Consequently, optical techniques have been introduced as safe and functional alternatives. Diffuse optical imaging is a non-invasive imaging technique that utilizes near-infrared light to examine biological tissues based on measuring the optical transmission and/or reflection at various locations on the tissue surface. In this paper, we propose a modified arrangement between the laser source and the detectors for distinguishing tumors from normal breast tissue. A three-dimensional model of the normal human breast with three types of tumors is developed using a COMSOL simulation software based on the finite element solution of Helmholtz equation to estimate optical fluence distribution. The breast model consists of four layers: skin, fat, glandular, and muscle, where the tumor is included in the glandular layer. Different wavelengths were used to determine the most proper wavelength for the discrimination between the normal tissue and tumor. The obtained results were verified using the receiver operating characteristic (ROC) method. The resultant fluence images show different features between normal breast and breast with tumor especially using 600-nm incident laser as demonstrated by the obtained ROC curves.

Multitask Deep Learning Reconstruction and Localization of Lesions in Limited Angle Diffuse Optical Tomography

Diffuse optical tomography (DOT) leverages near-infrared light propagation through tissue to assess its optical properties and identify abnormalities. DOT image reconstruction is an ill-posed problem due to the highly scattered photons in the medium and the smaller number of measurements compared to the number of unknowns. Limited-angle DOT reduces probe complexity at the cost of increased reconstruction complexity. Reconstructions are thus commonly marred by artifacts and, as a result, it is difficult to obtain an accurate reconstruction of target objects, e.g., malignant lesions. Reconstruction does not always ensure good localization of small lesions. Furthermore, conventional optimization-based reconstruction methods are computationally expensive, rendering them too slow for real-time imaging applications. Our goal is to develop a fast and accurate image reconstruction method using deep learning, where multitask learning ensures accurate lesion localization in addition to improved reconstruction. We apply spatial-wise attention and a distance transform based loss function in a novel multitask learning formulation to improve localization and reconstruction compared to single-task optimized methods. Given the scarcity of real-world sensor-image pairs required for training supervised deep learning models, we leverage physics-based simulation to generate synthetic datasets and use a transfer learning module to align the sensor domain distribution between in silico and real-world data, while taking advantage of cross-domain learning. Applying our method, we find that we can reconstruct and localize lesions faithfully while allowing real-time reconstruction. We also demonstrate that the present algorithm can reconstruct multiple cancer lesions. The results demonstrate that multitask learning provides sharper and more accurate reconstruction.

Preclinical evaluation of a clinical prototype transrectal diffuse optical tomography system for monitoring photothermal therapy of focal prostate cancer

SIGNIFICANCE

Our work demonstrates in preclinical models that continuous-wave transrectal diffuse optical tomography (TRDOT) can be used to accurately monitor photothermal therapy (PTT) and, in particular, the progression of the photocoagulation boundary toward the rectum. When used in patients, this should prevent rectal damage during PTT, thereby achieving maximum treatment efficacy while ensuring safety, using a technology platform suitable for wide dissemination.

AIM

We aim to validate that TRDOT measurements analyzed using a shape-based image-reconstruction algorithm (SBDOT) allow localization of the photocoagulation boundary during PTT within ±1  mm toward the rectum in the transverse plane.

APPROACH

TRDOT measurements were performed in tissue-simulating phantoms, ex vivo tissues, and an in vivo canine prostate model. The accuracy and sensitivity of reconstructing the size and location of the coagulation zone were determined, based on changes in the tissue absorption and reduced scattering coefficients upon photocoagulation. The reconstruction also yields the native and coagulated tissue optical properties.

RESULTS

The TRDOT measurements and SBDOT reconstruction algorithm were confirmed to perform sufficiently well for clinical translation in PTT monitoring, recovering the location of the coagulation boundary within ±1  mm compared to the true value as determined by direct visualization postexcision and/or MRI.

CONCLUSIONS

Implementing previously described TRDOT instrumentation and SBDOT image reconstruction in different tissue models confirms the potential for clinincal translation, including required refinements of the system and reconstruction algorithm.

Model equations of light scattering properties and a characteristic time of light propagation for polydisperse colloidal suspensions at different volume fractions

We developed model equations of light scattering properties and a characteristic time of light propagation for polydisperse colloidal suspensions at different volume fractions. By the model equations, we examined numerical results using the first-order (dependent) scattering theory (FST) and radiative transfer theory in 600-980 nm wavelength. The model equations efficiently treat the interference of electric fields scattered from colloidal particles by a single effective coefficient, providing fast computation. Meanwhile, the FST provides accurate but complicated treatment. We found the interference effects on the scattering properties and characteristic time depend linearly on wavelength. Dimensionless analysis showed a simple mechanism of the interference effects, independently of wavelength and source-detector distance.

Numerical Study of Near-Infrared Light Propagation in Aqueous Alumina Suspensions Using the Steady-State Radiative Transfer Equation and Dependent Scattering Theory

Understanding light propagation in liquid phantoms, such as colloidal suspensions, involves fundamental research of near-infrared optical imaging and spectroscopy for biological tissues. Our objective is to numerically investigate light propagation in the alumina colloidal suspensions with the mean alumina particle diameter of 55 nm at the volume fraction range 1–20%. We calculated the light scattering properties using the dependent scattering theory (DST) on a length scale comparable to the optical wavelength. We calculated the steady-state radiative transfer and photon diffusion equations (RTE and PDE) using the DST results based on the finite difference method in a length scale of the mean free path. The DST calculations showed that the scattering and reduced scattering coefficients become more prominent at a higher volume fraction. The anisotropy factor is almost zero at all the volume fractions, meaning the scattering is isotropic. The comparative study of the RTE with the PDE showed that the diffusion approximation holds at the internal region with all the volume fractions and the boundary region with the volume fraction higher than 1%. Our findings suggest the usefulness of the PDE as a light propagation model for the alumina suspensions rather than the RTE, which provides accurate but complicated computation.

Accuracy of retrieving optical properties from liquid tissue phantoms using a single integrating sphere

Using integrating spheres (ISs) in conjunction with the inverse adding-doubling algorithm (IAD) offers a well-established, rigorous protocol for determining optical absorption (μ_a) and reduced scattering (μs') coefficients of thin, optically homogeneous, turbid media. Here, we report the performance and use of a single IS system for experimentally retrieving optical properties in phantom media whose optical properties were well controlled. The IS system was used to measure the total reflectance and transmittance between 500 and 800 nm in liquid phantoms that were prepared to span a wide range of optical scattering and absorption coefficients. Measurements on phantoms were sequentially made using one of two broadband light sources-a halogen lamp or a supercontinuum laser. We report on the accuracy of IAD-derived optical coefficients using IS measurements made on phantoms-directly or by employing one of two previously reported correction methods. The first (sample-substitution error) correction was experimentally achieved while the second used Monte Carlo-based corrections with IAD. When experimentally calibrated reflectance and transmittance values were directly used as inputs to the IAD, mean absolute errors in recovered optical coefficients were larger than 0.4cm^-1 for absorption and more than 6cm^-1 for scattering across all phantoms and wavelengths measured. These errors reduced to 0.06-0.17cm^-1 and 0.7-2cm^-1 for μ_a and μs', respectively, with the use of corrections. Choice of light sources used, sample geometry (relative to optical coefficients), signal-to-noise of measurements, and the selection of correction methods are discussed.

Objective Numerical Evaluation of Diffuse, Optically Reconstructed Images Using Structural Similarity Index

Diffuse optical tomography is emerging as a non-invasive optical modality used to evaluate tissue information by obtaining the optical properties' distribution. Two procedures are performed to produce reconstructed absorption and reduced scattering images, which provide structural information that can be used to locate inclusions within tissues with the assistance of a known light intensity around the boundary. These methods are referred to as a forward problem and an inverse solution. Once the reconstructed image is obtained, a subjective measurement is used as the conventional way to assess the image. Hence, in this study, we developed an algorithm designed to numerically assess reconstructed images to identify inclusions using the structural similarity (SSIM) index. We compared four SSIM algorithms with 168 simulated reconstructed images involving the same inclusion position with different contrast ratios and inclusion sizes. A multiscale, improved SSIM containing a sharpness parameter (MS-ISSIM-S) was proposed to represent the potential evaluation compared with the human visible perception. The results indicated that the proposed MS-ISSIM-S is suitable for human visual perception by demonstrating a reduction of similarity score related to various contrasts with a similar size of inclusion; thus, this metric is promising for the objective numerical assessment of diffuse, optically reconstructed images.

Fast non-line-of-sight imaging with high-resolution and wide field of view using synthetic wavelength holography

The presence of a scattering medium in the imaging path between an object and an observer is known to severely limit the visual acuity of the imaging system. We present an approach to circumvent the deleterious effects of scattering, by exploiting spectral correlations in scattered wavefronts. Our Synthetic Wavelength Holography (SWH) method is able to recover a holographic representation of hidden targets with sub-mm resolution over a nearly hemispheric angular field of view. The complete object field is recorded within 46 ms, by monitoring the scattered light return in a probe area smaller than 6 cm × 6 cm. This unique combination of attributes opens up a plethora of new Non-Line-of-Sight imaging applications ranging from medical imaging and forensics, to early-warning navigation systems and reconnaissance. Adapting the findings of this work to other wave phenomena will help unlock a wider gamut of applications beyond those envisioned in this paper.

Recent Technical Progression in Photoacoustic Imaging—Towards Using Contrast Agents and Multimodal Techniques

For combining optical and ultrasonic imaging methodologies, photoacoustic imaging (PAI) is the most important and successful hybrid technique, which has greatly contributed to biomedical research and applications. Its theoretical background is based on the photoacoustic effect, whereby a modulated or pulsed light is emitted into tissue, which selectively absorbs the optical energy of the light at optical wavelengths. This energy produces a fast thermal expansion in the illuminated tissue, generating pressure waves (or photoacoustic waves) that can be detected by ultrasonic transducers. Research has shown that optical absorption spectroscopy offers high optical sensitivity and contrast for ingredient determination, for example, while ultrasound has demonstrated good spatial resolution in biomedical imaging. Photoacoustic imaging combines these advantages, i.e., high contrast through optical absorption and high spatial resolution due to the low scattering of ultrasound in tissue. In this review, we focus on advances made in PAI in the last five years and present categories and key devices used in PAI techniques. In particular, we highlight the continuously increasing imaging depth achieved by PAI, particularly when using exogenous reagents. Finally, we discuss the potential of combining PAI with other imaging techniques.

Grey-box modeling and hypothesis testing of functional near-infrared spectroscopy-based cerebrovascular reactivity to anodal high-definition tDCS in healthy humans

Transcranial direct current stimulation (tDCS) has been shown to evoke hemodynamics response; however, the mechanisms have not been investigated systematically using systems biology approaches. Our study presents a grey-box linear model that was developed from a physiologically detailed multi-compartmental neurovascular unit model consisting of the vascular smooth muscle, perivascular space, synaptic space, and astrocyte glial cell. Then, model linearization was performed on the physiologically detailed nonlinear model to find appropriate complexity (Akaike information criterion) to fit functional near-infrared spectroscopy (fNIRS) based measure of blood volume changes, called cerebrovascular reactivity (CVR), to high-definition (HD) tDCS. The grey-box linear model was applied on the fNIRS-based CVR during the first 150 seconds of anodal HD-tDCS in eleven healthy humans. The grey-box linear models for each of the four nested pathways starting from tDCS scalp current density that perturbed synaptic potassium released from active neurons for Pathway 1, astrocytic transmembrane current for Pathway 2, perivascular potassium concentration for Pathway 3, and voltage-gated ion channel current on the smooth muscle cell for Pathway 4 were fitted to the total hemoglobin concentration (tHb) changes from optodes in the vicinity of 4x1 HD-tDCS electrodes as well as on the contralateral sensorimotor cortex. We found that the tDCS perturbation Pathway 3 presented the least mean square error (MSE, median <2.5%) and the lowest Akaike information criterion (AIC, median -1.726) from the individual grey-box linear model fitting at the targeted-region. Then, minimal realization transfer function with reduced-order approximations of the grey-box model pathways was fitted to the ensemble average tHb time series. Again, Pathway 3 with nine poles and two zeros (all free parameters), provided the best Goodness of Fit of 0.0078 for Chi-Square difference test of nested pathways. Therefore, our study provided a systems biology approach to investigate the initial transient hemodynamic response to tDCS based on fNIRS tHb data. Future studies need to investigate the steady-state responses, including steady-state oscillations found to be driven by calcium dynamics, where transcranial alternating current stimulation may provide frequency-dependent physiological entrainment for system identification. We postulate that such a mechanistic understanding from system identification of the hemodynamics response to transcranial electrical stimulation can facilitate adequate delivery of the current density to the neurovascular tissue under simultaneous portable imaging in various cerebrovascular diseases.

Optical breast atlas as a testbed for image reconstruction in optical mammography

We present two optical breast atlases for optical mammography, aiming to advance the image reconstruction research by providing a common platform to test advanced image reconstruction algorithms. Each atlas consists of five individual breast models. The first atlas provides breast vasculature surface models, which are derived from human breast dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) data using image segmentation. A finite element-based method is used to deform the breast vasculature models from their natural shapes to generate the second atlas, compressed breast models. Breast compression is typically done in X-ray mammography but also necessary for some optical mammography systems. Technical validation is presented to demonstrate how the atlases can be used to study the image reconstruction algorithms. Optical measurements are generated numerically with compressed breast models and a predefined configuration of light sources and photodetectors. The simulated data is fed into three standard image reconstruction algorithms to reconstruct optical images of the vasculature, which can then be compared with the ground truth to evaluate their performance.

Indocyanine Green-Camptothecin Co-Loaded Perfluorocarbon Double-Layer Nanocomposite: A Versatile Nanotheranostics for Photochemotherapy and FDOT Diagnosis of Breast Cancer

Breast cancer remains the most frequently diagnosed cancer and is the leading cause of neoplastic disease burden for females worldwide, suggesting that effective therapeutic and/or diagnostic strategies are still urgently needed. In this study, a type of indocyanine green (ICG) and camptothecin (CPT) co-loaded perfluorocarbon double-layer nanocomposite named ICPNC was developed for detection and photochemotherapy of breast cancer. The ICPNCs were designed to be surface modifiable for on-demand cell targeting and can serve as contrast agents for fluorescence diffuse optical tomography (FDOT). Upon near infrared (NIR) irradiation, the ICPNCs can generate a significantly increased production of singlet oxygen compared to free ICG, and offer a comparable cytotoxicity with reduced chemo-drug dosage. Based on the results of animal study, we further demonstrated that the ICPNCs ([ICG]/[CPT] = 40-/7.5-μM) in association with 1-min NIR irradiation (808 nm, 6 W/cm²) can provide an exceptional anticancer effect to the MDA-MB-231 tumor-bearing mice whereby the tumor size was significantly reduced by 80% with neither organ damage nor systemic toxicity after a 21-day treatment. Given a number of aforementioned merits, we anticipate that the developed ICPNC is a versatile theranostic nanoagent which is highly promising to be used in the clinic.

An Adjustable Dark-Field Acoustic-Resolution Photoacoustic Imaging System with Fiber Bundle-Based Illumination

Photoacoustic (PA) imaging has become one of the major imaging methods because of its ability to record structural information and its high spatial resolution in biological tissues. Current commercialized PA imaging instruments are limited to varying degrees by their bulky size (i.e., the laser or scanning stage) or their use of complex optical components for light delivery. Here, we present a robust acoustic-resolution PA imaging system that consists of four adjustable optical fibers placed 90° apart around a 50 MHz high-frequency ultrasound (US) transducer. In the compact design concept of the PA probe, the relative illumination parameters (i.e., angles and fiber size) can be adjusted to fit different imaging applications in a single setting. Moreover, this design concept involves a user interface built in MATLAB. We first assessed the performance of our imaging system using in vitro phantom experiments. We further demonstrated the in vivo performance of the developed system in imaging (1) rat ear vasculature, (2) real-time cortical hemodynamic changes in the superior sagittal sinus (SSS) during left-forepaw electrical stimulation, and (3) real-time cerebral indocyanine green (ICG) dynamics in rats. Collectively, this alignment-free design concept of a compact PA probe without bulky optical lens systems is intended to satisfy the diverse needs in preclinical PA imaging studies.

Photoacoustic imaging for the monitoring of local changes in oxygen saturation following an adrenaline injection in human forearm skin

Clinical monitoring of blood oxygen saturation (sO₂) is traditionally performed using optical techniques, such as pulse oximetry and diffuse reflectance spectroscopy (DRS), which lack spatial resolution. Photoacoustic imaging (PAI) is a rapidly developing biomedical imaging technique that is superior to previous techniques in that it combines optical excitation and acoustic detection, providing a map of chromophore distribution in the tissue. Hitherto, PAI has primarily been used in preclinical studies, and only a few studies have been performed in patients. Its ability to measure sO₂ with spatial resolution during local vasoconstriction after adrenaline injection has not yet been investigated. Using PAI and spectral unmixing we characterize the heterogeneous change in sO₂ after injecting a local anesthetic containing adrenaline into the dermis on the forearm of seven healthy subjects. In comparison to results obtained using DRS, we highlight contrasting results obtained between the two methods arising due to the so-called 'window effect' caused by a reduced blood flow in the superficial vascular plexus. The results demonstrate the importance of spatially resolving sO₂ and the ability of PAI to assess the tissue composition in different layers of the skin.

Spatial light interference microscopy: principle and applications to biomedicine

In this paper, we review spatial light interference microscopy (SLIM), a common-path, phase-shifting interferometer, built onto a phase-contrast microscope, with white-light illumination. As one of the most sensitive quantitative phase imaging (QPI) methods, SLIM allows for speckle-free phase reconstruction with sub-nanometer path-length stability. We first review image formation in QPI, scattering, and full-field methods. Then, we outline SLIM imaging from theory and instrumentation to diffraction tomography. Zernike's phase-contrast microscopy, phase retrieval in SLIM, and halo removal algorithms are discussed. Next, we discuss the requirements for operation, with a focus on software developed in-house for SLIM that enables high-throughput acquisition, whole slide scanning, mosaic tile registration, and imaging with a color camera. We introduce two methods for solving the inverse problem using SLIM, white-light tomography, and Wolf phase tomography. Lastly, we review the applications of SLIM in basic science and clinical studies. SLIM can study cell dynamics, cell growth and proliferation, cell migration, mass transport, etc. In clinical settings, SLIM can assist with cancer studies, reproductive technology, blood testing, etc. Finally, we review an emerging trend, where SLIM imaging in conjunction with artificial intelligence brings computational specificity and, in turn, offers new solutions to outstanding challenges in cell biology and pathology.

Fundamentals of Optical Imaging

Optical imaging, which possesses noninvasive and high-resolution features for biomedical imaging, has been used to study various biological samples, from in vitro cells, ex vivo tissue, to in vivo imaging of living organism. Furthermore, optical imaging also covers a very wide scope of spatial scale, from submicron sized organelles to macro-scale live biological samples, enabling it a powerful tool for biomedical studies. Before introducing these superior optical imaging methods to researchers, first of all, it is necessary to present the basic concept of light-matter interactions such as absorption, scattering, and fluorescence, which can be used as the imaging contrast and also affect the imaging quality. And then the working mechanism of various imaging modalities including fluorescence microscopy, confocal microscopy, multiphoton microscopy, super-resolution microscopy, optical coherence tomography (OCT), diffuse optical tomography (DOT), etc. will be presented. Meanwhile, the main features and typical bioimaging applications of these optical imaging technologies are discussed. Finally, the perspective of future optical imaging methods is presented. The aim of this chapter is to introduce the background and principle of optical imaging for grasping the mechanism of advanced optical imaging modalities introduced in the following chapters.

The Use of Supercontinuum Laser Sources in Biomedical Diffuse Optics: Unlocking the Power of Multispectral Imaging

Optical techniques based on diffuse optics have been around for decades now and are making their way into the day-to-day medical applications. Even though the physics foundations of these techniques have been known for many years, practical implementation of these technique were hindered by technological limitations, mainly from the light sources and/or detection electronics. In the past 20 years, the developments of supercontinuum laser (SCL) enabled to unlock some of these limitations, enabling the development of system and methodologies relevant for medical use, notably in terms of spectral monitoring. In this review, we focus on the use of SCL in biomedical diffuse optics, from instrumentation and methods developments to their use for medical applications. A total of 95 publications were identified, from 1993 to 2021. We discuss the advantages of the SCL to cover a large spectral bandwidth with a high spectral power and fast switching against the disadvantages of cost, bulkiness, and long warm up times. Finally, we summarize the utility of using such light sources in the development and application of diffuse optics in biomedical sciences and clinical applications.

Real-Time Handheld Probe Tracking and Image Formation Using Digital Frequency-Domain Diffuse Optical Spectroscopy

OBJECTIVE

Frequency-domain diffuse optical spectroscopic imaging (FD-DOS) is a non-invasive method for measuring absolute concentrations of tissue chromophores such as oxy- and deoxy-hemoglobin in vivo. The utility of FD-DOS for clinical applications such as monitoring chemotherapy response in breast cancer has previously been demonstrated, but challenges for further clinical translation, such as slow acquisition speed and lack of user feedback, remain. Here, we propose a new high speed FD-DOS instrument that allows users to freely acquire measurements over the tissue surface, and is capable of rapidly imaging large volumes of tissue.

METHODS

We utilize 3D monocular probe tracking combined with custom digital FD-DOS hardware and a high-speed data processing pipeline for the instrument. Results are displayed during scanning over the surface of the sample using a probabilistic Monte Carlo light propagation model.

RESULTS

We show this instrument can measure absorption and scattering coefficients with an error of 7% and 1% respectively, with 0.7 mm positional accuracy. We demonstrate the equivalence of our visualization methodology with a standard interpolation approach, and demonstrate two proof-of-concept in vivo results showing superficial vasculature in the human forearm and surface contrast in a healthy human breast.

CONCLUSION

Our new FD-DOS system is able to compute chromophore concentrations in real-time (1.5 Hz) in vivo.

SIGNIFICANCE

This method has the potential to improve the quality of FD-DOS image scans while reducing measurement times for a variety of clinical applications.

Neuroimaging of depression with diffuse optical tomography during repetitive transcranial magnetic stimulation

Repetitive transcranial magnetic stimulation (rTMS) is an effective and safe treatment for depression; however, its potential has likely been hindered due to non-optimized targeting, unclear ideal stimulation parameters, and lack of information regarding how the brain is physiologically responding during and after stimulation. While neuroimaging is ideal for obtaining such critical information, existing modalities have been limited due to poor resolutions, along with significant noise interference from the electromagnetic spectrum. In this study, we used a novel diffuse optical tomography (DOT) device in order to advance our understanding of the neurophysiological effects of rTMS in depression. Healthy and depressed subjects aged 18–70 were recruited. Treatment parameters were standardized with targeting of the left dorsolateral prefrontal cortex with a magnetic field intensity of 100% of motor threshold, pulse frequency of 10 per second, a 4 s stimulation time and a 26 s rest time. DOT imaging was simultaneously acquired from the contralateral dorsolateral prefrontal cortex. Six healthy and seven depressed subjects were included for final analysis. Hemoglobin changes and volumetric three-dimensional activation patterns were successfully captured. Depressed subjects were observed to have a delayed and less robust response to rTMS with a decreased volume of activation compared to healthy subjects. In this first-in-human study, we demonstrated the ability of DOT to safely and reliably capture and compare cortical response patterns to rTMS in depressed and healthy subjects. We introduced this emerging optical functional imaging modality as a novel approach to investigating targeting, new treatment parameters, and physiological effects of rTMS in depression.

Deconvolution of hemodynamic responses along the cortical surface using personalized functional near infrared spectroscopy

In functional near infrared spectroscopy (fNIRS), deconvolution analysis of oxy and deoxy-hemoglobin concentration changes allows estimating specific hemodynamic response functions (HRF) elicited by neuronal activity, taking advantage of the fNIRS excellent temporal resolution. Diffuse optical tomography (DOT) is also becoming the new standard reconstruction procedure as it is more accurate than the modified Beer Lambert law approach at the sensor level. The objective of this study was to assess the relevance of HRF deconvolution after DOT constrained along the cortical surface. We used local personalized fNIRS montages which consists in optimizing the position of fNIRS optodes to ensure maximal sensitivity to subject specific target brain regions. We carefully evaluated the accuracy of deconvolution when applied after DOT, using realistic simulations involving several HRF models at different signal to noise ratio (SNR) levels and on real data related to motor and visual tasks in healthy subjects and from spontaneous pathological activity in one patient with epilepsy. We demonstrated that DOT followed by deconvolution was able to accurately recover a large variability of HRFs over a large range of SNRs. We found good performances of deconvolution analysis for SNR levels usually encountered in our applications and we were able to reconstruct accurately the temporal dynamics of HRFs in real conditions.

In vivo photoacoustic assessment of the oxygen saturation changes in the human radial artery: a preliminary study associated with age

SIGNIFICANCE

We demonstrate the potential of probing the sO2 change under blood flow in vivo using photoacoustic (PA) imaging and sheds light on the complex relationship between RBC aggregation and oxygen delivery.

AIM

To conduct in vivo assessments of the sO2 in the radial artery of healthy volunteers and simultaneously probe the relation between the sO2 and hemodynamic behavior such as red blood cell (RBC) aggregation.

APPROACH

The effects of PA-based measurements of blood hemodynamics were studied as a function of the subjects' age (20s, 30s, and 40s). The pulsatile blood flow in the human radial artery of 12 healthy subjects was imaged in the 700 to 900 nm optical wavelength range using a linear array-based PA system.

RESULTS

The PA power when blood velocity is minimum (Pamax) was larger than the one attained at maximum blood velocity (Pamin), consistent with predictions based on the cyclical variation of RBC aggregation during pulsatile flow. The difference between Pamin and Pamax at 800 nm (ΔPa800) increased with age (1.7, 2.2, and 2.6 dB for age group of 20s, 30s, and 40s, respectively). The sO2 computed from Pamax was larger than the one from Pamin.

CONCLUSIONS

The ΔPa800 increased with participant age. The ΔPa800 metric could be a surrogate of noninvasively monitoring the age-induced changes in RBC aggregation. The sO2 change during a cycle of pulsatile blood flow also increased with age, demonstrating that RBC aggregation can affect the sO2 change.