High-density functional diffuse optical tomography based on frequency-domain measurements improves image quality and spatial resolution

Mapping cortical activations underlying covert and overt language production using high-density diffuse optical tomography

Gold standard neuroimaging modalities such as functional magnetic resonance imaging (fMRI), positron emission tomography (PET), and more recently electrocorticography (ECoG) have provided profound insights regarding the neural mechanisms underlying the processing of language, but they are limited in applications involving naturalistic language production especially in developing brains, during face-to-face dialogues, or as a brain-computer interface. High-density diffuse optical tomography (HD-DOT) provides high-fidelity mapping of human brain function with comparable spatial resolution to that of fMRI but in a silent and open scanning environment similar to real-life social scenarios. Therefore, HD-DOT has potential to be used in naturalistic settings where other neuroimaging modalities are limited. While HD-DOT has been previously validated against fMRI for mapping the neural correlates underlying language comprehension and covert (i.e., "silent") language production, HD-DOT has not yet been established for mapping the cortical responses to overt (i.e., "out loud") language production. In this study, we assessed the brain regions supporting a simple hierarchy of language tasks: silent reading of single words, covert production of verbs, and overt production of verbs in normal hearing right-handed native English speakers (n = 33). First, we found that HD-DOT brain mapping is resilient to movement associated with overt speaking. Second, we observed that HD-DOT is sensitive to key activations and deactivations in brain function underlying the perception and naturalistic production of language. Specifically, statistically significant results were observed that show recruitment of regions in occipital, temporal, motor, and prefrontal cortices across all three tasks after performing stringent cluster-extent based thresholding. Our findings lay the foundation for future HD-DOT studies of imaging naturalistic language comprehension and production during real-life social interactions and for broader applications such as presurgical language assessment and brain-machine interfaces.

Effects of atlas-based anatomy on modelled light transport in the neonatal head

Objective.Diffuse optical tomography (DOT) provides a relatively convenient method for imaging haemodynamic changes related to neuronal activity on the cerebral cortex. Due to practical challenges in obtaining anatomical images of neonates, an anatomical framework is often created from an age-appropriate atlas model, which is individualized to the subject based on measurements of the head geometry. This work studies the approximation error arising from using an atlas instead of the neonate's own anatomical model.Approach.We consider numerical simulations of frequency-domain (FD) DOT using two approaches, Monte Carlo simulations and diffusion approximation via finite element method, and observe the variation in (1) the logarithm of amplitude and phase shift measurements, and (2) the corresponding inner head sensitivities (Jacobians), due to varying segmented anatomy. Varying segmentations are sampled by registering 165 atlas models from a neonatal database to the head geometry of one individual selected as the reference model. Prior to the registration, we refine the segmentation of the cerebrospinal fluid (CSF) by separating the CSF into two physiologically plausible layers.Main results.In absolute measurements, a considerable change in the grey matter or extracerebral tissue absorption coefficient was found detectable over the anatomical variation. In difference measurements, a small local 10%-increase in brain absorption was clearly detectable in the simulated measurements over the approximation error in the Jacobians, despite the wide range of brain maturation among the registered models.Significance.Individual-level atlas models could potentially be selected within several weeks in gestational age in DOT difference imaging, if an exactly age-appropriate atlas is not available. The approximation error method could potentially be implemented to improve the accuracy of atlas-based imaging. The presented CSF segmentation algorithm could be useful also in other model-based imaging modalities. The computation of FD Jacobians is now available in the widely-used Monte Carlo eXtreme software.

Review of recent advances in frequency-domain near-infrared spectroscopy technologies [Invited]

Over the past several decades, near-infrared spectroscopy (NIRS) has become a popular research and clinical tool for non-invasively measuring the oxygenation of biological tissues, with particular emphasis on applications to the human brain. In most cases, NIRS studies are performed using continuous-wave NIRS (CW-NIRS), which can only provide information on relative changes in chromophore concentrations, such as oxygenated and deoxygenated hemoglobin, as well as estimates of tissue oxygen saturation. Another type of NIRS known as frequency-domain NIRS (FD-NIRS) has significant advantages: it can directly measure optical pathlength and thus quantify the scattering and absorption coefficients of sampled tissues and provide direct measurements of absolute chromophore concentrations. This review describes the current status of FD-NIRS technologies, their performance, their advantages, and their limitations as compared to other NIRS methods. Significant landmarks of technological progress include the development of both benchtop and portable/wearable FD-NIRS technologies, sensitive front-end photonic components, and high-frequency phase measurements. Clinical applications of FD-NIRS technologies are discussed to provide context on current applications and needed areas of improvement. The review concludes by providing a roadmap toward the next generation of fully wearable, low-cost FD-NIRS systems.

Novel data types for frequency-domain diffuse optical spectroscopy and imaging of tissues: characterization of sensitivity and contrast-to-noise ratio for absorption perturbations

In frequency-domain (FD) diffuse optics it is known that the phase of photon-density waves (ϕ) has a stronger deep-to-superficial sensitivity ratio to absorption perturbations than the alternate current (AC) amplitude, or the direct current intensity (DC). This work is an attempt to find FD data types that feature similar or even better sensitivity and/or contrast-to-noise for deeper absorption perturbations than phase. One way is to start from the definition of characteristic function (Xt(ω)) of the photon's arrival time (t) and combining the real (ℜ ( X t ( ω ) ) = A C D C c o s ( ϕ ) ) and imaginary parts (ℑ [ X t ( ω ) ] = A C D C s i n ( ϕ ) ) with phase to yield new data types. These new data types enhance the role of higher order moments of the probability distribution of the photon's arrival time t. We study the contrast-to-noise and sensitivity features of these new data types not only in the single-distance arrangement (traditionally used in diffuse optics), but we also consider the spatial gradients, which we named dual-slope arrangements. We have identified six data types that for typical values of the optical properties of tissues and depths of interest, have better sensitivity or contrast-to-noise features than phase data and that can be used to enhance the limits of imaging of tissue in FD near infrared spectroscopy (NIRS). For example, one promising data type is ϕ - ℑ [ X t ( ω ) ] which shows, in the single-distance source-detector arrangement, an increase of deep-to-superficial sensitivity ratio with respect to phase by 41% and 27% at a source-detector separation of 25 and 35 mm, respectively. The same data type also shows an increase of contrast-to noise up to 35% with respect to phase when the spatial gradients of the data are considered.

System Derived Spatial-Temporal CNN for High-Density fNIRS BCI

An intuitive and generalisable approach to spatial-temporal feature extraction for high-density (HD) functional Near-Infrared Spectroscopy (fNIRS) brain-computer interface (BCI) is proposed, demonstrated here using Frequency-Domain (FD) fNIRS for motor-task classification. Enabled by the HD probe design, layered topographical maps of Oxy/deOxy Haemoglobin changes are used to train a 3D convolutional neural network (CNN), enabling simultaneous extraction of spatial and temporal features. The proposed spatial-temporal CNN is shown to effectively exploit the spatial relationships in HD fNIRS measurements to improve the classification of the functional haemodynamic response, achieving an average F1 score of 0.69 across seven subjects in a mixed subjects training scheme, and improving subject-independent classification as compared to a standard temporal CNN.

Dual-slope imaging of cerebral hemodynamics with frequency-domain near-infrared spectroscopy

Significance

This work targets the contamination of optical signals by superficial hemodynamics, which is one of the chief hurdles in non-invasive optical measurements of the human brain.

Aim

To identify optimal source-detector distances for dual-slope (DS) measurements in frequency-domain (FD) near-infrared spectroscopy (NIRS) and demonstrate preferential sensitivity of DS imaging to deeper tissue (brain) versus superficial tissue (scalp).

Approach

Theoretical studies (in-silico) based on diffusion theory in two-layered and in homogeneous scattering media. In-vivo demonstrations of DS imaging of the human brain during visual stimulation and during systemic blood pressure oscillations.

Results

The mean distance (between the two source-detector distances needed for DS) is the key factor for depth sensitivity. In-vivo imaging of the human occipital lobe with FD NIRS and a mean distance of 31 mm indicated: (1) greater hemodynamic response to visual stimulation from FD phase versus intensity, and from DS versus single-distance (SD); (2) hemodynamics from FD phase and DS mainly driven by blood flow, and hemodynamics from SD intensity mainly driven by blood volume.

Conclusions

DS imaging with FD NIRS may suppress confounding contributions from superficial hemodynamics without relying on data at short source-detector distances. This capability can have significant implications for non-invasive optical measurements of the human brain.

Multi-modulated frequency domain high density diffuse optical tomography

Frequency domain (FD) high density diffuse optical tomography (HD-DOT) utilising varying or combined modulation frequencies (mFD) has shown to theoretically improve the imaging accuracy as compared to conventional continuous wave (CW) measurements. Using intensity and phase data from a solid inhomogeneous phantom (NEUROPT) with three insertable rods containing different contrast anomalies, at modulation frequencies of 78 MHz, 141 MHz and 203 MHz, HD-DOT is applied and quantitatively evaluated, showing that mFD outperforms FD and CW for both absolute (iterative) and temporal (linear) tomographic imaging. The localization error (LOCA), full width half maximum (FWHM) and effective resolution (ERES) were evaluated. Across all rods, the LOCA of mFD was 61.3% better than FD and 106.1% better than CW. For FWHM, CW was 6.0% better than FD and mFD and for ERES, mFD was 1.20% better than FD and 9.83% better than CW. Using mFD data is shown to minimize the effect of inherently noisier FD phase data whilst maximising its strengths through improved contrast.

Lawrence K, Tachtsidis I, Tong Y, Torricelli A, Urner T, Wabnitz H, Wolf M, Wolf U, Xu S, Yang C, Yodh AG, Yücel MA, Zhou W. Optical imaging and spectroscopy for the study of the human brain: status report

This report is the second part of a comprehensive two-part series aimed at reviewing an extensive and diverse toolkit of novel methods to explore brain health and function. While the first report focused on neurophotonic tools mostly applicable to animal studies, here, we highlight optical spectroscopy and imaging methods relevant to noninvasive human brain studies. We outline current state-of-the-art technologies and software advances, explore the most recent impact of these technologies on neuroscience and clinical applications, identify the areas where innovation is needed, and provide an outlook for the future directions.

Lock-in functional near-infrared spectroscopy for measurement of the haemodynamic brain response

Here we show a method of the lock-in amplifying near-infrared signals originating within a human brain. It implies using two 90-degree rotated source-detector pairs fixed on a head surface. Both pairs have a joint sensitivity region located towards the brain. A direct application of the lock-in technique on both signals results in amplifying common frequency components, e.g. related to brain cortex stimulation and attenuating the rest, including all components not related to the stimulation: e.g. pulse, instrumental and biological noise or movement artefacts. This is a self-driven method as no prior assumptions are needed and the noise model is provided by the interfering signals themselves. We show the theory (classical modified Beer-Lambert law and diffuse optical tomography approaches), the algorithm implementation and tests on a finite element mathematical model and in-vivo on healthy volunteers during visual cortex stimulation. The proposed hardware and algorithm complexity suit the entire spectrum of (continuous wave, frequency domain, time-resolved) near-infrared spectroscopy systems featuring real-time, direct, robust and low-noise brain activity registration tool. As such, this can be of special interest in optical brain computer interfaces and high reliability/stability monitors of tissue oxygenation.

A Data Self-Calibration Method Based on High-Density Parallel Plate Diffuse Optical Tomography for Breast Cancer Imaging

When performing the diffuse optical tomography (DOT) of the breast, the mismatch between the forward model and the experimental conditions will significantly hinder the reconstruction accuracy. Therefore, the reference measurement is commonly used to calibrate the measured data before the reconstruction. However, it is complicated to customize corresponding reference phantoms based on the breast shape and background optical parameters of different subjects in clinical trials. Furthermore, although high-density (HD) DOT configuration has been proven to improve imaging quality, a large number of source-detector (SD) pairs also increase the difficulty of multi-channel correction. To enhance the applicability of the breast DOT, a data self-calibration method based on an HD parallel-plate DOT system is proposed in this paper to replace the conventional relative measurement on a reference phantom. The reference predicted data can be constructed directly from the measurement data with the support of the HD-DOT system, which has nearly a hundred sets of measurements at each SD distance. The proposed scheme has been validated by Monte Carlo (MC) simulation, breast-size phantom experiments, and clinical trials, exhibiting the feasibility in ensuring the quality of the DOT reconstruction while effectively reducing the complexity associated with relative measurements on reference phantoms.

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.

A scalable, multi-wavelength, broad bandwidth frequency-domain near-infrared spectroscopy platform for real-time quantitative tissue optical imaging

Frequency-domain near-infrared spectroscopy (FD-NIRS) provides quantitative noninvasive measurements of tissue optical absorption and scattering, as well as a safe and accurate method for characterizing tissue composition and metabolism. However, the poor scalability and high complexity of most FD-NIRS systems assembled to date have contributed to its limited clinical impact. To address these shortcomings, we present a scalable, digital-based FD-NIRS platform capable of measuring optical properties and tissue chromophore concentrations in real-time. The system provides single-channel FD-NIRS amplitude/phase, optical property, and chromophore data at a maximum display rate of 36.6 kHz, 17.9 kHz, and 10.2 kHz, respectively, and can be scaled to multiple channels as well as integrated into a handheld format. The entire system is enabled by several innovations including an ultra-high-speed k-nearest neighbor lookup table method (maximum of 250,000 inversions/s for a large 2500x700 table of absorption and reduced scattering coefficients), embedded FPGA and CPU high-speed co-processing, and high-speed data transfer (due to on-board processing). We show that our 6-wavelength, broad modulation bandwidth (1-400 MHz) system can be used to perform 2D high-density spatial mapping of optical properties and high speed quantification of hemodynamics.

Investigation of effect of modulation frequency on high-density diffuse optical tomography image quality

Significance: By incorporating multiple overlapping functional near-infrared spectroscopy (fNIRS) measurements, high-density diffuse optical tomography (HD-DOT) images human brain function with fidelity comparable to functional magnetic resonance imaging (fMRI). Previous work has shown that frequency domain high-density diffuse optical tomography (FD-HD-DOT) may further improve image quality over more traditional continuous wave (CW) HD-DOT. Aim: The effects of modulation frequency on image quality as obtainable with FD-HD-DOT is investigated through simulations with a realistic noise model of functional activations in human head models, arising from 11 source modulation frequencies between CW and 1000 MHz. Approach: Simulations were performed using five representative head models with an HD regular grid of 158 light sources and 166 detectors and an empirically derived noise model. Functional reconstructions were quantitatively assessed with multiple image quality metrics including the localization error (LE), success rate, full width at half maximum, and full volume at half maximum (FVHM). All metrics were evaluated against CW-based models. Results: Compared to CW, localization accuracy is improved by >40% throughout brain depths of 13 to 25 mm below the surface with 300 to 500 MHz modulation frequencies. Additionally, the reliable field of view in brain tissue is enlarged by 35% to 48% within an optimal frequency of 300 MHz after considering realistic noise, depending on the dynamic range of the system. Conclusions: These results point to the tremendous opportunities in further development of high bandwidth FD-HD-DOT system hardware for applications in human brain mapping.

The Validation of a Portable Functional NIRS System for Assessing Mental Workload

Portable functional near-infrared spectroscopy (fNIRS) systems have the potential to image the brain in naturalistic settings. Experimental studies are essential to validate such fNIRS systems. Working memory (WM) is a short-term active memory that is associated with the temporary storage and manipulation of information. The prefrontal cortex (PFC) brain area is involved in the processing of WM. We assessed the PFC brain during n-back WM tasks in a group of 25 college students using our laboratory-developed portable fNIRS system, WearLight. We designed an experimental protocol with 32 n-back WM task blocks with four different pseudo-randomized task difficulty levels. The hemodynamic response of the brain was computed from the experimental data and the evaluated brain responses due to these tasks. We observed the incremental mean hemodynamic activation induced by the increasing WM load. The left-PFC area was more activated in the WM task compared to the right-PFC. The task performance was seen to be related to the hemodynamic responses. The experimental results proved the functioning of the WearLight system in cognitive load imaging. Since the portable fNIRS system was wearable and operated wirelessly, it was possible to measure the cognitive load in the naturalistic environment, which could also lead to the development of a user-friendly brain–computer interface system.

Quantitative evaluation of frequency domain measurements in high density diffuse optical tomography

SIGNIFICANCE

High density diffuse optical tomography (HD-DOT) as applied in functional near-infrared spectroscopy (fNIRS) is largely limited to continuous wave (CW) data. Using a single modulation frequency, frequency domain (FD) HD-DOT has recently demonstrated better localization of focal activation as compared to CW data. We show that combining CW and FD measurements and multiple modulation frequencies increases imaging performance in fNIRS.

AIM

We evaluate the benefits of multiple modulation frequencies, combining different frequencies as well as CW data in fNIRS HD-DOT.

APPROACH

A layered model was used, with activation occurring within a cortex layer. CW and FD measurements were simulated at 78, 141, and 203 MHz with and without noise. The localization error, full width half maximum, and effective resolution were evaluated.

RESULTS

Across the average of the three metrics, at 141 MHz, FD performed 8.4% better than CW, and the combination of CW and FD was 21.7% better than CW. FD measurements at 203 MHz performed 5% better than 78 MHz. Moreover, the three combined modulation frequencies of FD and CW performed up to 3.92% better than 141 MHz alone.

CONCLUSIONS

We show that combining CW and FD measurements offers better performance than FD alone, with higher modulation frequencies increasing accuracy. Combining CW and FD measurements at multiple modulation frequencies yields the best overall performance.

Diffuse optical tomography in the human brain: A briefly review from the neurophysiology to its applications

The present work describes the use of noninvasive diffuse optical tomography (DOT) technology to measure hemodynamic changes, providing relevant information which helps to understand the basis of neurophysiology in the human brain. Advantages such as portability, direct measurements of hemoglobin state, temporal resolution, non‐restricted movements as occurs in magnetic resonance imaging (MRI) devices mean that DOT technology can be used in research and clinical fields. In this review we covered the neurophysiology, physical principles underlying optical imaging during tissue‐light interactions, and technology commonly used during the construction of a DOT device including the source‐detector requirements to improve the image quality. DOT provides 3D cerebral activation images due to complex mathematical models which describe the light propagation inside the tissue head. Moreover, we describe briefly the use of Bayesian methods for raw DOT data filtering as an alternative to linear filters widely used in signal processing, avoiding common problems such as the filter selection or a false interpretation of the results which is sometimes due to the interference of background physiological noise with neural activity.

Sensitivity of frequency-domain optical measurements to brain hemodynamics: simulations and human study of cerebral blood flow during hypercapnia

This study characterizes the sensitivity of noninvasive measurements of cerebral blood flow (CBF) by using frequency-domain near-infrared spectroscopy (FD-NIRS) and coherent hemodynamics spectroscopy (CHS). We considered six FD-NIRS methods: single-distance intensity and phase (SDI and SDϕ), single-slope intensity and phase (SSI and SSϕ), and dual-slope intensity and phase (DSI and DSϕ). Cerebrovascular reactivity (CVR) was obtained from the relative change in measured CBF during a step hypercapnic challenge. Greater measured values of CVR are assigned to a greater sensitivity to cerebral hemodynamics. In a first experiment with eight subjects, CVRSDϕ was greater than CVRSDI (p < 0.01), whereas CVRDSI and CVRDSϕ showed no significant difference (p > 0.5). In a second experiment with four subjects, a 5 mm scattering layer was added between the optical probe and the scalp tissue to increase the extracerebral layer thickness (L ec ), which caused CVRDSϕ to become significantly greater than CVRDSI (p < 0.05). CVRSS measurements yielded similar results as CVRDS measurements but with a greater variability, possibly resulting from instrumental artifacts in SS measurements. Theoretical simulations with two-layered media confirmed that, if the top (extracerebral) layer is more scattering than the bottom (brain) layer, the relative values of CVRDSI and CVRDSϕ depend on L ec . Specifically, the sensitivity to the brain is greater for DSI than DSϕ for a thin extracerebral layer (L ec < 13 mm), whereas it is greater for DSϕ than DSI for a thicker extracerebral layer.

Signal regression in frequency-domain diffuse optical tomography to remove superficial signal contamination

Significance: Signal contamination is a major hurdle in functional near-infrared spectroscopy (fNIRS) of the human head as the NIR signal is contaminated with the changes corresponding to superficial tissue, therefore occluding the functional information originating from the cerebral region. For continuous wave, this is generally handled through linear regression of the shortest source-detector (SD) distance intensity measurement from all of the signals. Although phase measurements utilizing frequency domain (FD) provide deeper tissue sampling, the use of the shortest SD distance phase measurement for regression of superficial signal contamination can lead to misleading results, therefore suppressing cortical signals. Aim: An approach for FD fNIRS that utilizes a short-separation intensity signal directly to regress both intensity and phase measurements, providing a better regression of superficial signal contamination from both data-types, is proposed. Approach: Simulated data from realistic models of the human head are used, and signal regression using both intensity and phase-based components of the FD fNIRS is evaluated. Results: Intensity-based phase regression achieves a suppression of superficial signal contamination by 68% whereas phase-based phase regression is only by 13%. Phase-based phase regression is also shown to generate false-positive signals from the cortex, which are not desirable. Conclusions: Intensity-based phase regression provides a better methodology for minimizing superficial signal contamination in FD fNIRS.

Enhanced spatiotemporal resolution imaging of neuronal activity using joint electroencephalography and diffuse optical tomography

Significance: Electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS) are both commonly used methodologies for neuronal source reconstruction. While EEG has high temporal resolution (millisecond-scale), its spatial resolution is on the order of centimeters. On the other hand, in comparison to EEG, fNIRS, or diffuse optical tomography (DOT), when used for source reconstruction, can achieve relatively high spatial resolution (millimeter-scale), but its temporal resolution is poor because the hemodynamics that it measures evolve on the order of several seconds. This has important neuroscientific implications: e.g., if two spatially close neuronal sources are activated sequentially with only a small temporal separation, single-modal measurements using either EEG or DOT alone would fail to resolve them correctly. Aim: We attempt to address this issue by performing joint EEG and DOT neuronal source reconstruction. Approach: We propose an algorithm that utilizes DOT reconstruction as the spatial prior of EEG reconstruction, and demonstrate the improvements using simulations based on the ICBM152 brain atlas. Results: We show that neuronal sources can be reconstructed with higher spatiotemporal resolution using our algorithm than using either modality individually. Further, we study how the performance of the proposed algorithm can be affected by the locations of the neuronal sources, and how the performance can be enhanced by improving the placement of EEG electrodes and DOT optodes. Conclusions: We demonstrate using simulations that two sources separated by 2.3-3.3 cm and 50 ms can be recovered accurately using the proposed algorithm by suitably combining EEG and DOT, but not by either in isolation. We also show that the performance can be enhanced by optimizing the electrode and optode placement according to the locations of the neuronal sources.

Cap-Based Transcranial Optical Tomography in an Awake Infant

Although Blood Oxygenation Level Dependent (BOLD) functional MRI (fMRI) is widely used to examine brain function in adults, the need for general anesthesia limits its practical utility in infants and small children. Functional Near-Infrared Spectroscopy - Diffuse Optical Tomography (fNIRS-DOT) imaging promises to be an alternative brain network imaging technique. Yet current versions of continuous-wave fNIRS-DOT systems are restricted to the cortical surface measurements and do not probe deep structures that are frequently injured especially in premature infants. Herein we report a transcranial near infrared optical imaging system, called Cap-based Transcranial Optical Tomography (CTOT) able to image whole brain hemodynamic activity with 3 seconds of data acquisition time. We show the system is capable of whole brain oxygenation mapping in an awake child, and that tomographically reconstructed static CTOT-derived oxy- and deoxygenated blood volumes are spatially correlated with the time-averaged BOLD fMRI volumes. By removing time bottlenecks in the current system, dynamic CTOT mapping should be possible, which would then enable evaluation of functional connectivity in awake infants.

Dynamic time domain near-infrared optical tomography based on a SPAD camera

In many clinical applications it is relevant to observe dynamic changes in oxygenation. Therefore the ability of dynamic imaging with time domain (TD) near-infrared optical tomography (NIROT) will be important. But fast imaging is a challenge. The data acquisition of our handheld TD NIROT system based on single photon avalanche diode (SPAD) camera and 11 light sources was consequently accelerated. We tested the system on a diffusive medium simulating tissue with a moving object embedded. With 3D image reconstruction, the moving object was correctly located using only 0.2 s exposure time per source.

Recent Developments in Instrumentation of Functional Near-Infrared Spectroscopy Systems

In the last three decades, the development and steady improvement of various optical technologies at the near-infrared region of the electromagnetic spectrum has inspired a large number of scientists around the world to design and develop functional near-infrared spectroscopy (fNIRS) systems for various medical applications. This has been driven further by the availability of new sources and detectors that support very compact and wearable system designs. In this article, we review fNIRS systems from the instrumentation point of view, discussing the associated challenges and state-of-the-art approaches. In the beginning, the fundamentals of fNIRS systems as well as light-tissue interaction at NIR are briefly introduced. After that, we present the basics of NIR systems instrumentation. Next, the recent development of continuous-wave, frequency-domain, and time-domain fNIRS systems are discussed. Finally, we provide a summary of these three modalities and an outlook into the future of fNIRS technology.

Design of a source-detector array for dual-slope diffuse optical imaging

We recently proposed a dual-slope technique for diffuse optical spectroscopy and imaging of scattering media. This technique requires a special configuration of light sources and optical detectors to create dual-slope sets. Here, we present methods for designing, optimizing, and building an optical imaging array that features m dual-slope sets to image n voxels. After defining the m × n matrix (S) that describes the sensitivity of the m dual-slope measurements to absorption perturbations in each of the n voxels, we formulate the inverse imaging problem in terms of the Moore-Penrose pseudoinverse matrix of S (S+). This approach allows us to introduce several measures of imaging performance: reconstruction accuracy (correct spatial mapping), crosstalk (incorrect spatial mapping), resolution (point spread function), and localization (offset between actual and reconstructed point perturbations). Furthermore, by considering the singular value decomposition formulation, we show the significance of visualizing the first m right singular vectors of S, whose linear combination generates the reconstructed map. We also describe methods to build a physical array using a three-layer mesh structure (two polyethylene films and polypropylene hook-and-loop fabric) embedded in silicone (PDMS). Finally, we apply these methods to design two arrays and choose one to construct. The chosen array consists of 16 illumination fibers, 10 detection fibers, and 27 dual-slope sets for dual-slope imaging optimized for the size of field of view and localization of absorption perturbations. This particular array is aimed at functional near-infrared spectroscopy of the human brain, but the methods presented here are of general applicability to a variety of devices and imaging scenarios.

Optimizing sensitivity and dynamic range of silicon photomultipliers for frequency-domain near infrared spectroscopy

Diffuse optical imaging and tomography based upon frequency-domain near-infrared spectroscopy (fdNIRS) is used to noninvasively measure tissue structure and function through quantitative absolute measurements of tissue optical absorption and scattering. Here we describe how utilizing a silicon photomultiplier (SiPM) detector for fdNIRS improves performance. We discuss the operation of SiPMs, how they differ from other fdNIRS photodetectors, and show theoretically that SiPMs offer similar sensitivity to photomultiplier tube (PMT) detectors while having a higher dynamic range and lower cost, size, and operating voltage. With respect to avalanche photodiode (APD) detectors, theoretical and experimental data shows drastically increased signal to noise ratio performance, up to 25dB on human breast, head, and muscle tissue. Finally, we extend the dynamic range (∼10dB) of the SiPM through a nonlinear calibration technique which reduced absorption error by a mean 16 percentage points.

The Various Oximetric Techniques Used for the Evaluation of Blood Oxygenation

Adequate oxygen delivery to a tissue depends on sufficient oxygen content in arterial blood and blood flow to the tissue. Oximetry is a technique for the assessment of blood oxygenation by measurements of light transmission through the blood, which is based on the different absorption spectra of oxygenated and deoxygenated hemoglobin. Oxygen saturation in arterial blood provides information on the adequacy of respiration and is routinely measured in clinical settings, utilizing pulse oximetry. Oxygen saturation, in venous blood (SvO₂) and in the entire blood in a tissue (StO₂), is related to the blood supply to the tissue, and several oximetric techniques have been developed for their assessment. SvO₂ can be measured non-invasively in the fingers, making use of modified pulse oximetry, and in the retina, using the modified Beer–Lambert Law. StO₂ is measured in peripheral muscle and cerebral tissue by means of various modes of near infrared spectroscopy (NIRS), utilizing the relative transparency of infrared light in muscle and cerebral tissue. The primary problem of oximetry is the discrimination between absorption by hemoglobin and scattering by tissue elements in the attenuation measurement, and the various techniques developed for isolating the absorption effect are presented in the current review, with their limitations.

Dual-slope imaging in highly scattering media with frequency-domain near-infrared spectroscopy

We present theoretical and experimental demonstrations of a novel, to the best of our knowledge, diffuse optical imaging method that is based on the concept of dual slopes (DS) in frequency-domain near-infrared spectroscopy. We consider a special array of sources and detectors that collects intensity (I) and phase (ϕ) data with multiple DS sets. We have recently shown that DSϕ reflectance data features a deeper sensitivity with respect to DSI reflectance data. Here, for the first time, we describe a DS imaging approach based on the Moore-Penrose inverse of the sensitivity matrix for multiple DS data sets. Using a circular 8-source/9-detector array that generates 16 DS data sets at source-detector distances in the range 20-40 mm, we show that DSI images are more sensitive to superficial (<5mm) perturbations, whereas DSϕ images are more sensitive to deeper (>10mm) perturbations in highly scattering media.

Performance of measurands in time-domain optical brain imaging: depth selectivity versus contrast-to-noise ratio

Time-domain optical brain imaging techniques introduce a number of different measurands for analyzing absorption changes located deep in the tissue, complicated by superficial absorption changes and noise. We implement a method that allows analysis, quantitative comparison and performance ranking of measurands under various conditions - including different values of reduced scattering coefficient, thickness of the superficial layer, and source-detector separation. Liquid phantom measurements and Monte Carlo simulations were carried out in two-layered geometry to acquire distributions of times of flight of photons and to calculate the total photon count, mean time of flight, variance, photon counts in time windows and ratios of photon counts in different time windows. Quantitative comparison of performance was based on objective metrics: relative contrast, contrast-to-noise ratio (CNR) and depth selectivity. Moreover, the product of CNR and depth selectivity was used to rank the overall performance and to determine the optimal source-detector separation for each measurand. Variance ranks the highest under all considered conditions.

Frequency-Domain Techniques for Cerebral and Functional Near-Infrared Spectroscopy

This article reviews the basic principles of frequency-domain near-infrared spectroscopy (FD-NIRS), which relies on intensity-modulated light sources and phase-sensitive optical detection, and its non-invasive applications to the brain. The simpler instrumentation and more straightforward data analysis of continuous-wave NIRS (CW-NIRS) accounts for the fact that almost all the current commercial instruments for cerebral NIRS have embraced the CW technique. However, FD-NIRS provides data with richer information content, which complements or exceeds the capabilities of CW-NIRS. One example is the ability of FD-NIRS to measure the absolute optical properties (absorption and reduced scattering coefficients) of tissue, and thus the absolute concentrations of oxyhemoglobin and deoxyhemoglobin in brain tissue. This article reviews the measured values of such optical properties and hemoglobin concentrations reported in the literature for animal models and for the human brain in newborns, infants, children, and adults. We also review the application of FD-NIRS to functional brain studies that focused on slower hemodynamic responses to brain activity (time scale of seconds) and faster optical signals that have been linked to neuronal activation (time scale of 100 ms). Another example of the power of FD-NIRS data is related to the different regions of sensitivity featured by intensity and phase data. We report recent developments that take advantage of this feature to maximize the sensitivity of non-invasive optical signals to brain tissue relative to more superficial extracerebral tissue (scalp, skull, etc.). We contend that this latter capability is a highly appealing quality of FD-NIRS, which complements absolute optical measurements and may result in significant advances in the field of non-invasive optical sensing of the brain.