Complete head cerebral sensitivity mapping for diffuse correlation spectroscopy using subject-specific magnetic resonance imaging models

Beneath the surface: revealing deep-tissue blood flow in human subjects with massively parallelized diffuse correlation spectroscopy

Significance

Diffuse correlation spectroscopy (DCS) allows label-free, non-invasive investigation of microvascular dynamics deep within tissue, such as cerebral blood flow (CBF). However, the signal-to-noise ratio (SNR) in DCS limits its effective cerebral sensitivity in adults, in which the depth to the brain, through the scalp and skull, is substantially larger than in infants.

Aim

Therefore, we aim to increase its SNR and, ultimately, its sensitivity to CBF through new DCS techniques.

Approach

We present an in vivo demonstration of parallelized DCS (PDCS) to measure cerebral and muscular blood flow in healthy adults. Our setup employs an innovative array with hundreds of thousands single photon avalanche diodes (SPAD) in a 500 × 500 grid to boost SNR by averaging all independent pixel measurements. We tested this device on different total pixel counts and frame rates. A secondary, smaller array was used for reference measurements from shallower tissue at lower source-detector-separation (SDS).

Results

The new system can measure pulsatile blood flow in cerebral and muscular tissue, at up to 4 cm SDS, while maintaining a similar measurement noise as compared with a previously published 32 × 32 PDCS system at 1.5 cm SDS. Data from a cohort of 15 adults provide strong experimental evidence for functional CBF activity during a cognitive memory task and allowed analysis of pulse markers. Additional control experiments on muscular blood flow in the forearm with a different technical configuration provide converging evidence for the efficacy of this technique.

Conclusions

Our results outline successful PDCS measurements with large SPAD arrays to enable detect CBF in human adults. The ongoing development of SPAD camera technology is expected to result in larger and faster detectors in the future. In combination with new data processing techniques, tailored for the sparse signal of binary photon detection events in SPADs, this could lead to even greater SNR increase and ultimately greater depth sensitivity of PDCS.

Heuristic absorption calculation in bilayered media from a white Monte Carlo dataset

Monte Carlo (MC) simulations can adequately describe photon migration in layered media; however, storing and querying the resulting dataset may be computationally prohibitive when detailed path data are needed for each photon trajectory. A heuristic approach that significantly reduces the stored information to the average path length traveled by the photons detected in each layer is proposed. Its accuracy is evaluated by comparing it with the exact time point spread function (TPSF) for a bilayered medium. This method, almost exact for small variations in absorption, is potentially useful to provide a small dataset for lookup tables to be used in inverse problems.

Revisiting equivalent optical properties for cerebrospinal fluid to improve diffusion-based modeling accuracy in the brain

Significance

The diffusion approximation (DA) is used in functional near-infrared spectroscopy (fNIRS) studies despite its known limitations due to the presence of cerebrospinal fluid (CSF). Many of these studies rely on a set of empirical CSF optical properties, recommended by a previous simulation study, that were not selected for the purpose of minimizing DA modeling errors.

Aim

We aim to directly quantify the accuracy of DA solutions in brain models by comparing those with the gold-standard solutions produced by the mesh-based Monte Carlo (MMC), based on which we derive updated recommendations.

Approach

For both a five-layer head and Colin27 atlas models, we obtain DA solutions by independently sweeping the CSF absorption (μ a ) and reduced scattering (μ s ' ) coefficients. Using an MMC solution with literature CSF optical properties as a reference, we compute the errors for surface fluence, total brain sensitivity, and brain energy deposition, and identify the optimized settings where such error is minimized.

Results

Our results suggest that previously recommended CSF properties can cause significant errors (8.7% to 52%) in multiple tested metrics. By simultaneously sweepingμ a andμ s ' , we can identify infinite numbers of solutions that can exactly match DA with MMC solutions for any single tested metric. Furthermore, it is also possible to simultaneously minimize multiple metrics at multiple source/detector separations, leading to our updated recommendation of settingμ s ' = 0.15    mm - 1 while maintaining physiologicalμ a for CSF in DA simulations.

Conclusions

Our updated recommendation of CSF equivalent optical properties can greatly reduce the model mismatches between DA and MMC solutions at multiple metrics without sacrificing computational speed. We also show that it is possible to eliminate such a mismatch for a single or a pair of metrics of interest.

Lawrence K. All-optics technique for monitoring absolute cerebral blood flow: validation against magnetic resonance imaging perfusion

Significance

The ability to monitor cerebral blood flow (CBF) at the bedside is essential to managing critical-care patients with neurological emergencies. Diffuse correlation spectroscopy (DCS) is ideal because it is non-invasive, portable, and inexpensive. We investigated a near-infrared spectroscopy (NIRS) approach for converting DCS measurements into physiological units of blood flow.

Aim

Using magnetic resonance imaging perfusion as a reference, we investigated the accuracy of absolute CBF measurements from a bolus-tracking NIRS method that used transient hypoxia as a flow tracer and hypercapnia-induced increases in CBF measured by DCS.

Approach

Twelve participants (7 female, 28 ± 6 years) completed a hypercapnia protocol with simultaneous CBF recordings from DCS and arterial spin labeling (ASL). Nine participants completed the transient hypoxia protocol while instrumented with time-resolved NIRS. The estimate of baseline CBF was subsequently used to calibrate hypercapnic DCS data.

Results

Moderately strong correlations at baseline ( slope = 0.79 andR 2 = 0.59 ) and during hypercapnia ( slope = 0.90 andR 2 = 0.58 ) were found between CBF values from calibrated DCS and ASL (range 34 to 85    mL / 100    g / min ).

Conclusions

Results demonstrated the feasibility of an all-optics approach that can both quantify CBF and perform continuous perfusion monitoring.

How the heterogeneity of the severely injured brain affects hybrid diffuse optical signals: case examples and guidelines

Significance

A shortcoming of the routine clinical use of diffuse optics (DO) in the injured head has been that the results from commercial near-infrared spectroscopy-based devices are not reproducible, often give physiologically invalid values, and differ among systems. Besides the limitations due to the physics of continuous-wave light sources, one culprit is the head heterogeneity and the underlying morphological and functional abnormalities of the probed tissue.

Aim

The aim is to investigate the effect that different tissue alterations in the damaged head have on DO signals and provide guidelines to avoid data misinterpretation.

Approach

DO measurements and computed tomography scans were acquired on brain-injured patients. The relationship between the signals and the underlying tissue types was classified on a case-by-case basis.

Results

Examples and suggestions to establish quality control routines were provided. The findings suggested guidelines for carrying out DO measurements and speculations toward improved devices.

Conclusions

We advocate for the standardization of the DO measurements to secure a role for DO in neurocritical care. We suggest that blind measurements are unacceptably problematic due to confounding effects and care using a priori and a posteriori quality control routines that go beyond an assessment of the signal-to-noise ratio that is typically utilized.

Modified Beer-Lambert algorithm to measure pulsatile blood flow, critical closing pressure, and intracranial hypertension

We introduce a frequency-domain modified Beer-Lambert algorithm for diffuse correlation spectroscopy to non-invasively measure flow pulsatility and thus critical closing pressure (CrCP). Using the same optical measurements, CrCP was obtained with the new algorithm and with traditional nonlinear diffusion fitting. Results were compared to invasive determination of intracranial pressure (ICP) in piglets (n = 18). The new algorithm better predicted ICP elevations; the area under curve (AUC) from logistic regression analysis was 0.85 for ICP ≥ 20 mmHg. The corresponding AUC for traditional analysis was 0.60. Improved diagnostic performance likely results from better filtering of extra-cerebral tissue contamination and measurement noise.

Time-domain diffuse correlation spectroscopy at large source detector separation for cerebral blood flow recovery

Time-domain diffuse correlation spectroscopy (td-DCS) enables the depth discrimination in tissue's blood flow recovery, considering the fraction of photons detected with higher time of flight (TOF) and longer pathlength through the tissue. However, the recovery result depends on factors such as the instrument response function (IRF), analyzed TOF gate start time, gate width and the source-detector separation (SDS). In this research we evaluate the performance of the td-DCS technique at three SDSs of 1.5, 2 and 2.5 cm to recover cerebral blood flow (CBF). To do that we presented comprehensive characterization of the td-DCS system through a series of phantom experiments. First by quality metrices such as coefficient of variation and contrast-to-noise ratios, we identified optimal time gate(s) of the TOF to extract dynamics of particles. Then using sensitivity metrices, each SDS ability to detect dynamics of particles in superficial and deeper layer was evaluated. Finally, td-DCS at each SDS was tested on healthy volunteers during cuff occlusion test and breathing tasks. According to phantom measurements, the sensitivity to estimate perfusion within the deep layer located at depth of 1.5 cm from the surface can be increased more than two times when the SDS increases from 1.5 cm to 2.5 cm.

Lawrence K, Chui J, Diop M. Assessing changes in regional cerebral hemodynamics in adults with a high-density full-head coverage time-resolved near-infrared spectroscopy device

Significance

Cerebral oximeters have the potential to detect abnormal cerebral blood oxygenation to allow for early intervention. However, current commercial systems have two major limitations: (1) spatial coverage of only the frontal region, assuming that surgery-related hemodynamic effects are global and (2) susceptibility to extracerebral signal contamination inherent to continuous-wave near-infrared spectroscopy (NIRS).

Aim

This work aimed to assess the feasibility of a high-density, time-resolved (tr) NIRS device (Kernel Flow) to monitor regional oxygenation changes across the cerebral cortex during surgery.

Approach

The Flow system was assessed using two protocols. First, digital carotid compression was applied to healthy volunteers to cause a rapid oxygenation decrease across the ipsilateral hemisphere without affecting the contralateral side. Next, the system was used on patients undergoing shoulder surgery to provide continuous monitoring of cerebral oxygenation. In both protocols, the improved depth sensitivity of trNIRS was investigated by applying moment analysis. A dynamic wavelet filtering approach was also developed to remove observed temperature-induced signal drifts.

Results

In the first protocol (28 ± 5 years; five females, five males), hair significantly impacted regional sensitivity; however, the enhanced depth sensitivity of trNIRS was able to separate brain and scalp responses in the frontal region. Regional sensitivity was improved in the clinical study given the age-related reduction in hair density of the patients (65 ± 15 years; 14 females, 13 males). In five patients who received phenylephrine to treat hypotension, different scalp and brain oxygenation responses were apparent, although no regional differences were observed.

Conclusions

The Kernel Flow has promise as an intraoperative neuromonitoring device. Although regional sensitivity was affected by hair color and density, enhanced depth sensitivity of trNIRS was able to resolve differences in scalp and brain oxygenation responses in both protocols.

Introduction to the Optics and the Brain 2023 feature issue

A feature issue is being presented by a team of guest editors containing papers based on contributed submissions including studies presented at Optics and the Brain, held April 24-27, 2023 as part of Optica Biophotonics Congress: Optics in the Life Sciences, in Vancouver, Canada.

Validation of the Openwater wearable optical system: cerebral hemodynamic monitoring during a breath-hold maneuver

Significance

Bedside cerebral blood flow (CBF) monitoring has the potential to inform and improve care for acute neurologic diseases, but technical challenges limit the use of existing techniques in clinical practice.

Aim

Here, we validate the Openwater optical system, a novel wearable headset that uses laser speckle contrast to monitor microvascular hemodynamics.

Approach

We monitored 25 healthy adults with the Openwater system and concurrent transcranial Doppler (TCD) while performing a breath-hold maneuver to increase CBF. Relative blood flow (rBF) was derived from changes in speckle contrast, and relative blood volume (rBV) was derived from changes in speckle average intensity.

Results

A strong correlation was observed between beat-to-beat optical rBF and TCD-measured cerebral blood flow velocity (CBFv), R = 0.79 ; the slope of the linear fit indicates good agreement, 0.87 (95% CI: 0.83 - 0.92 ). Beat-to-beat rBV and CBFv were also strongly correlated, R = 0.72 , but as expected the two variables were not proportional; changes in rBV were smaller than CBFv changes, with linear fit slope of 0.18 (95% CI: 0.17 to 0.19). Further, strong agreement was found between rBF and CBFv waveform morphology and related metrics.

Conclusions

This first in vivo validation of the Openwater optical system highlights its potential as a cerebral hemodynamic monitor, but additional validation is needed in disease states.

Comparing the performance potential of speckle contrast optical spectroscopy and diffuse correlation spectroscopy for cerebral blood flow monitoring using Monte Carlo simulations in realistic head geometries

Significance

The non-invasive measurement of cerebral blood flow based on diffuse optical techniques has seen increased interest as a research tool for cerebral perfusion monitoring in critical care and functional brain imaging. Diffuse correlation spectroscopy (DCS) and speckle contrast optical spectroscopy (SCOS) are two such techniques that measure complementary aspects of the fluctuating intensity signal, with DCS quantifying the temporal fluctuations of the signal and SCOS quantifying the spatial blurring of a speckle pattern. With the increasing interest in the use of these techniques, a thorough comparison would inform new adopters of the benefits of each technique.

Aim

We systematically evaluate the performance of DCS and SCOS for the measurement of cerebral blood flow.

Approach

Monte Carlo simulations of dynamic light scattering in an MRI-derived head model were performed. For both DCS and SCOS, estimates of sensitivity to cerebral blood flow changes, coefficient of variation of the measured blood flow, and the contrast-to-noise ratio of the measurement to the cerebral perfusion signal were calculated. By varying complementary aspects of data collection between the two methods, we investigated the performance benefits of different measurement strategies, including altering the number of modes per optical detector, the integration time/fitting time of the speckle measurement, and the laser source delivery strategy.

Results

Through comparison across these metrics with simulated detectors having realistic noise properties, we determine several guiding principles for the optimization of these techniques and report the performance comparison between the two over a range of measurement properties and tissue geometries. We find that SCOS outperforms DCS in terms of contrast-to-noise ratio for the cerebral blood flow signal in the ideal case simulated here but note that SCOS requires careful experimental calibrations to ensure accurate measurements of cerebral blood flow.

Conclusion

We provide design principles by which to evaluate the development of DCS and SCOS systems for their use in the measurement of cerebral blood flow.

Comparison of diffuse correlation spectroscopy analytical models for measuring cerebral blood flow in adults

Significance

Although multilayer analytical models have been proposed to enhance brain sensitivity of diffuse correlation spectroscopy (DCS) measurements of cerebral blood flow, the traditional homogeneous model remains dominant in clinical applications. Rigorous in vivo comparison of these analytical models is lacking.

Aim

We compare the performance of different analytical models to estimate a cerebral blood flow index (CBFi) with DCS in adults.

Approach

Resting-state data were obtained on a cohort of 20 adult patients with subarachnoid hemorrhage. Data at 1 and 2.5 cm source-detector separations were analyzed with the homogenous, two-layer, and three-layer models to estimate scalp blood flow index and CBFi. The performance of each model was quantified via fitting convergence, fit stability, brain-to-scalp flow ratio (BSR), and correlation with transcranial Doppler ultrasound (TCD) measurements of cerebral blood flow velocity in the middle cerebral artery (MCA).

Results

The homogeneous model has the highest pass rate (100%), lowest coefficient of variation (CV) at rest (median [IQR] at 1 Hz of 0.18 [0.13, 0.22]), and most significant correlation with MCA blood flow velocities (R s = 0.59 , p = 0.010 ) compared with both the two- and three-layer models. The multilayer model pass rate was significantly correlated with extracerebral layer thicknesses. Discarding datasets with non-physiological BSRs increased the correlation between DCS measured CBFi and TCD measured MCA velocities for all models.

Conclusions

We found that the homogeneous model has the highest pass rate, lowest CV at rest, and most significant correlation with MCA blood flow velocities. Results from the multilayer models should be taken with caution because they suffer from lower pass rates and higher coefficients of variation at rest and can converge to non-physiological values for CBFi. Future work is needed to validate these models in vivo, and novel approaches are merited to improve the performance of the multimodel models.

Correlation and asynchronization of electroencephalogram and cerebral blood flow in active and passive stimulations

Objective.Real-time brain monitoring is of importance for intraoperative surgeries and intensive care unit, in order to take timely clinical interventions. Electroencephalogram (EEG) is a conventional technique for recording neural excitations (e.g. brain waves) in the cerebral cortex, and near infrared diffuse correlation spectroscopy (DCS) is an emerging technique that can directly measure the cerebral blood flow (CBF) in microvasculature system. Currently, the relationship between the neural activities and cerebral hemodynamics that reflects the vasoconstriction features of cerebral vessels, especially under both active and passive situation, has not been elucidated thus far, which triggers the motivation of this study.Approach.We used the verbal fluency test as an active cognitive stimulus to the brain, and we manipulated blood pressure changes as a passive challenge to the brain. Under both protocols, the CBF and EEG responses were longitudinally monitored throughout the cerebral stimulus. Power spectrum approaches were applied the EEG signals and compared with CBF responses.Main results.The results show that the EEG response was significantly faster and larger in amplitude during the active cognitive task, when compared to the CBF, but with larger individual variability. By contrast, CBF is more sensitive when response to the passive task, and with better signal stability. We also found that there was a correlation (p< 0.01,r= 0.866,R2= 0.751) between CBF and EEG in initial response during the active task, but no significant correlation (p> 0.05) was found during the passive task. The similar relations were also found between regional brain waves and blood flow.Significance.The asynchronization and correlation between the two measurements indicates the necessity of monitoring both variables for comprehensive understanding of cerebral physiology. Deep exploration of their relationships provides promising implications for DCS/EEG integration in the diagnosis of various neurovascular and psychiatric diseases.

Impaired cerebrovascular reactivity in pediatric sickle cell disease using diffuse correlation spectroscopy

Cerebrovascular reactivity (CVR), defined as the ability of cerebral vasculature to dilate in response to a vasodilatory stimulus, is an integral mechanism in brain homeostasis that is thought to be impaired in sickle cell disease (SCD). This study used diffuse correlation spectroscopy and a simple breath-hold stimulus to quantify CVR non-invasively in a cohort of 12 children with SCD and 14 controls. Median [interquartile range] CVR was significantly decreased in SCD compared to controls (2.03 [1.31, 2.44] versus 3.49 [3.00, 4.11] %/mmHg, p = 0.028). These results suggest DCS may provide a feasible means to routinely monitor CVR impairments in pediatric SCD.

scatterBrains: an open database of human head models and companion optode locations for realistic Monte Carlo photon simulations

Significance

Monte Carlo (MC) simulations are currently the gold standard in the near-infrared and diffuse correlation spectroscopy (NIRS/DCS) communities for generating light transport paths through tissue. However, realistic and diverse models that capture complex tissue layers are not easily available to all; moreover, manually placing optodes on such models can be tedious and time consuming. Such limitations may hinder the adoption of representative models for basic simulations and the use of these models for large-scale simulations, e.g., for training machine learning algorithms.

Aim

We aim to provide the NIRS/DCS communities with an open-source, user-friendly database of morphologically and optically realistic head models, as well as a succinct software pipeline to prepare these models for mesh-based Monte Carlo simulations of light transport.

Approach

Sixteen anatomical models were created from segmented T1-weighted magnetic resonance imaging (MRI) head scans and converted to tetrahedral mesh volumes. Approximately 800 companion scalp surface locations were distributed on each model, comprising full head coverage. A pipeline was created to place custom source and optical detectors at each location, and guidance is provided on how to use these parameters to set up MC simulations.

Results

The models, head surface locations, and all associated code are freely available under the scatterBrains project on Github.

Conclusions

The NIRS/DCS community benefits from having shared resources for conducting MC simulations on realistic head geometries. We hope this will make MRI-based head models and virtual optode placement easily accessible to all. Contributions to the database are welcome and encouraged.

Erratum: Complete head cerebral sensitivity mapping for diffuse correlation spectroscopy using subject-specific magnetic resonance imaging models: errata

[This corrects the article on p. 1131 in vol. 13, PMID: 35414976.].

Non-invasive diffuse optical monitoring of cerebral physiology in an adult swine-model of impact traumatic brain injury

In this study, we used diffuse optics to address the need for non-invasive, continuous monitoring of cerebral physiology following traumatic brain injury (TBI). We combined frequency-domain and broadband diffuse optical spectroscopy with diffuse correlation spectroscopy to monitor cerebral oxygen metabolism, cerebral blood volume, and cerebral water content in an established adult swine-model of impact TBI. Cerebral physiology was monitored before and after TBI (up to 14 days post injury). Overall, our results suggest that non-invasive optical monitoring can assess cerebral physiologic impairments post-TBI, including an initial reduction in oxygen metabolism, development of cerebral hemorrhage/hematoma, and brain swelling.

Advanced Neuromonitoring Modalities on the Horizon: Detection and Management of Acute Brain Injury in Children

Timely detection and monitoring of acute brain injury in children is essential to mitigate causes of injury and prevent secondary insults. Increasing survival in critically ill children has emphasized the importance of neuroprotective management strategies for long-term quality of life. In emergent and critical care settings, traditional neuroimaging modalities, such as computed tomography and magnetic resonance imaging (MRI), remain frontline diagnostic techniques to detect acute brain injury. Although detection of structural and anatomical abnormalities remains crucial, advanced MRI sequences assessing functional alterations in cerebral physiology provide unique diagnostic utility. Head ultrasound has emerged as a portable neuroimaging modality for point-of-care diagnosis via assessments of anatomical and perfusion abnormalities. Application of electroencephalography and near-infrared spectroscopy provides the opportunity for real-time detection and goal-directed management of neurological abnormalities at the bedside. In this review, we describe recent technological advancements in these neurodiagnostic modalities and elaborate on their current and potential utility in the detection and management of acute brain injury.

Portable, high speed blood flow measurements enabled by long wavelength, interferometric diffuse correlation spectroscopy (LW-iDCS)

Diffuse correlation spectroscopy (DCS) is an optical technique that can be used to characterize blood flow in tissue. The measurement of cerebral hemodynamics has arisen as a promising use case for DCS, though traditional implementations of DCS exhibit suboptimal signal-to-noise ratio (SNR) and cerebral sensitivity to make robust measurements of cerebral blood flow in adults. In this work, we present long wavelength, interferometric DCS (LW-iDCS), which combines the use of a longer illumination wavelength (1064 nm), multi-speckle, and interferometric detection, to improve both cerebral sensitivity and SNR. Through direct comparison with long wavelength DCS based on superconducting nanowire single photon detectors, we demonstrate an approximate 5× improvement in SNR over a single channel of LW-DCS in the measured blood flow signals in human subjects. We show equivalence of extracted blood flow between LW-DCS and LW-iDCS, and demonstrate the feasibility of LW-iDCS measured at 100 Hz at a source-detector separation of 3.5 cm. This improvement in performance has the potential to enable robust measurement of cerebral hemodynamics and unlock novel use cases for diffuse correlation spectroscopy.

Diffuse correlation spectroscopy: current status and future outlook

Diffuse correlation spectroscopy (DCS) has emerged as a versatile, noninvasive method for deep tissue perfusion assessment using near-infrared light. A broad class of applications is being pursued in neuromonitoring and beyond. However, technical limitations of the technology as originally implemented remain as barriers to wider adoption. A wide variety of approaches to improve measurement performance and reduce cost are being explored; these include interferometric methods, camera-based multispeckle detection, and long path photon selection for improved depth sensitivity. We review here the current status of DCS technology and summarize future development directions and the challenges that remain on the path to widespread adoption.

Influence of oversimplifying the head anatomy on cerebral blood flow measurements with diffuse correlation spectroscopy

Significance

Diffuse correlation spectroscopy (DCS) is an emerging optical modality for non-invasive assessment of an index of regional cerebral blood flow. By the nature of this noninvasive measurement, light must pass through extracerebral layers (i.e., skull, scalp, and cerebral spinal fluid) before detection at the tissue surface. To minimize the contribution of these extracerebral layers to the measured signal, an analytical model has been developed that treats the head as a series of three parallel and infinitely extending slabs (mimicking scalp, skull, and brain). The three-layer model has been shown to provide a significant improvement in cerebral blood flow estimation over the typically used model that treats the head as a bulk homogenous medium. However, the three-layer model is still a gross oversimplification of the head geometry that ignores head curvature, the presence of cerebrospinal fluid (CSF), and heterogeneity in layer thickness.

Aim

Determine the influence of oversimplifying the head geometry on cerebral blood flow estimated with the three-layer model.

Approach

Data were simulated with Monte Carlo in a four-layer slab medium and a three-layer sphere medium to isolate the influence of CSF and curvature, respectively. Additionally, simulations were performed on magnetic resonance imaging (MRI) head templates spanning a wide-range of ages. Simulated data were fit to both the homogenous and three-layer model for CBF. Finally, to mitigate the errors in potential CBF estimation due to the difficulty in defining layer thickness, we investigated an approach to identify an equivalent, "optimized" thickness via a pressure modulation.

Results

Both head curvature and failing to account for CSF lead to significant errors in the estimation of CBF. However, the effect of curvature and CSF on relative changes in CBF is minimal. Further, we found that CBF was underestimated in all MRI-templates, although the magnitude of these underestimations was highly influenced by small variations in the source and detector optode positioning. The optimized thickness obtained from pressure modulation did not improve estimation accuracy of CBF, although it did significantly improve the estimation accuracy of relative changes in CBF.

Conclusions

In sum, these findings suggest that the three-layer model holds promise for improving estimation of relative changes in cerebral blood flow; however, estimations of absolute cerebral blood flow with the approach should be viewed with caution given that it is difficult to account for appreciable sources of error, such as curvature and CSF.

Numerical approach to quantify depth-dependent blood flow changes in real-time using the diffusion equation with continuous-wave and time-domain diffuse correlation spectroscopy

Diffuse correlation spectroscopy (DCS) is a non-invasive optical technique that can measure brain perfusion by quantifying temporal intensity fluctuations of multiply scattered light. A primary limitation for accurate quantitation of cerebral blood flow (CBF) is the fact that experimental measurements contain information about both extracerebral scalp blood flow (SBF) as well as CBF. Separating CBF from SBF is typically achieved using multiple source-detector channels when using continuous-wave (CW) light sources, or more recently with use of time-domain (TD) techniques. Analysis methods that account for these partial volume effects are often employed to increase CBF contrast. However, a robust, real-time analysis procedure that can separate and quantify SBF and CBF with both traditional CW and TD-DCS measurements is still needed. Here, we validate a data analysis procedure based on the diffusion equation in layered media capable of quantifying both extra- and cerebral blood flow in the CW and TD. We find that the model can quantify SBF and CBF coefficients with less than 5% error compared to Monte Carlo simulations using a 3-layered brain model in both the CW and TD. The model can accurately fit data at a rate of <10 ms for CW data and <250 ms for TD data when using a least-squares optimizer.

Laser speckle simulation tool based on stochastic differential equations for bio imaging applications

Laser speckle-based blood flow imaging is a well-accepted and widely used method for pre-clinical and clinical applications. Although it was introduced as a method to measure only superficial blood flow (< 1mm depth), several recently introduced variants resulted in measuring deep tissue blood flow (a few cm) as well. A means of simulating laser speckles is often necessary for the analysis and development of these imaging modalities, as evident from many such attempts towards developing simulation tools in the past. Such methods often employ Fourier transforms or statistical tools to simulate speckles with desired statistical properties. We present the first method to use a stochastic differential equation to generate laser speckles with a pre-determined probability density function and a temporal auto-correlation. The method allows the choice of apriori gamma distribution along with simple exponential or more complex temporal auto-correlation statistics for simulated speckles, making it suitable for different blood flow profiles. In contrast to the existing methods that often generate speckles associated with superficial flow, we simulate both superficial and diffuse speckles leading to applications in deep tissue blood flow imaging. In addition, we have also incorporated appropriate models for noise associated with the detectors to simulate realistic speckles. We have validated our model by comparing the simulated speckles with those obtained from in-vivo studies in mice and healthy human subject.

High-density diffuse correlation tomography with enhanced depth localization and minimal surface artefacts

A spatially weighted filter applied to both the measurement and the Jacobian is proposed for high-density diffuse correlation tomography (DCT) to remove unwanted extracerebral interferences and artefacts along with better depth localization in the reconstructed blood flow images. High-density DCT is implemented by appropriate modification of recently introduced Multi-speckle Diffuse Correlation Spectroscopy (M-DCS) system. Additionally, we have used autocorrelation measurements at multiple delay-times in an iterative manner to improve the reconstruction results. The proposed scheme has been validated by simulations, phantom experiments and in-vivo human experiments.

Functional Time Domain Diffuse Correlation Spectroscopy

Time-domain diffuse correlation spectroscopy (TD-DCS) offers a novel approach to high-spatial resolution functional brain imaging based on the direct quantification of cerebral blood flow (CBF) changes in response to neural activity. However, the signal-to-noise ratio (SNR) offered by previous TD-DCS instruments remains a challenge to achieving the high temporal resolution needed to resolve perfusion changes during functional measurements. Here we present a next-generation optimized functional TD-DCS system that combines a custom 1,064 nm pulse-shaped, quasi transform-limited, amplified laser source with a high-resolution time-tagging system and superconducting nanowire single-photon detectors (SNSPDs). System characterization and optimization was conducted on homogenous and two-layer intralipid phantoms before performing functional CBF measurements in six human subjects. By acquiring CBF signals at over 5 Hz for a late gate start time of the temporal point spread function (TPSF) at 15 mm source-detector separation, we demonstrate for the first time the measurement of blood flow responses to breath-holding and functional tasks using TD-DCS.