Journal of Healthcare Engineering
Volume 2019 (2019), Article ID 3750495, 11 pages
https://doi.org/10.1155/2019/3750495

Review Article

Portable Near-Infrared Technologies and Devices for Noninvasive Assessment of Tissue Hemodynamics

Lin Hou,¹ Yinqiu Liu,¹ Lixia Qian,² Yucong Zheng,² Jinnan Gao,² Wenxing Cao,³ and Yu Shang¹

¹North University of China, No. 3 Xueyuan Road, Taiyuan 030051, China
²Shanxi Dayi Hospital, No. 99 Longcheng Street, Taiyuan 030032, China
³University of Electronic Science and Technology of China, No. 2006 Xiyuan Street, Chengdu 611731, China

Correspondence should be addressed to Yu Shang; yushang@nuc.edu.cn

Received 18 July 2018; Revised 24 December 2018; Accepted 14 January 2019; Published 12 February 2019

Academic Editor: Valentina Hartwig

Copyright © 2019 Lin Hou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Tissue hemodynamics, including the blood flow, oxygenation, and oxygen metabolism, are closely associated with many diseases. As one of the portable optical technologies to explore human physiology and assist in healthcare, near-infrared diffuse optical spectroscopy (NIRS) for tissue oxygenation measurement has been developed for four decades. In recent years, a dynamic NIRS technology, namely, diffuse correlation spectroscopy (DCS), has been emerging as a portable tool for tissue blood flow measurement. In this article, we briefly describe the basic principle and algorithms for static NIRS and dynamic NIRS (i.e., DCS). Then, we elaborate on the NIRS instrumentation, either commercially available or custom-made, as well as their applications to physiological studies and clinic. The extension of NIRS/DCS from spectroscopy to imaging was depicted, followed by introductions of advanced algorithms that were recently proposed. The future prospective of the NIRS/DCS and their feasibilities for routine utilization in hospital is finally discussed.

1. Introduction

Many diseases would lead to abnormal status in microvasculature blood flow, oxygenation, and oxygen metabolism [1, 2]. For example, ischemic stroke causes reduction in both blood flow and oxygenation [3]. Malignant tumor leads to higher blood flow, higher oxygen consumption rate, and lower oxygen status [4, 5]. For the people with obstructive sleep apnea, large fluctuations in microvasculature blood flow and total hemoglobin concentration were discovered [6].

There are few technologies currently available to noninvasively monitor the oxygen status in tissue. PO₂ is a technology capable of quantifying the oxygen pressure via a tiny needle [7]. This invasive pattern of measurement, however, is unfavorable to healthy population. As a conventional technology for tissue blood oxygenation measurement including the oxy-, deoxy-, total hemoglobin concentration ([HbO₂], [Hb], and THC), and oxygen saturation (StO₂), near-infrared diffuse optical spectroscopy (NIRS) has been developed for more than four decades, since 1970s [8]. NIRS oximeter was often utilized in clinical and physiological studies [8–10] for detecting the relative deep tissue (up to several centimeters), owing to the low tissue absorption in near-infrared (NIR) range (650 to 900 nm). However, the oxygenation alone only reflects the static balance between the oxygen supply and consumption [11]. Blood flow is such a parameter that elucidates how fast and efficient the oxygen is carried in the blood and supplied to the tissue. Studies show that blood flow is highly sensitive to pathophysiological alteration [12], thus could be a potential indicator for early detection of many diseases. Additionally, the blood flow reflects the reactivity patterns of vascular vasodilatation [3], thus assisting in evaluation of vessel functions. Therefore, investigation of blood flow, as well as its correlations with pathophysiological consequences, has significant implications for diagnosis and medications. Furthermore, the combination of tissue blood flow and oxygenation allows for estimate of the oxygen metabolic rate [11, 13], a valuable parameter closely associated with tissue pathology.

At present, the low-cost modalities to noninvasively and fast measure the blood flow in deep tissue are not established. Methods currently available for assessment of blood flow include Doppler ultrasound [14], laser Doppler [15], and perfusion-weighted magnetic resonance imaging (perfusion MRI) [16]. Doppler ultrasound can only assess the blood velocity in single and large vessel [14, 17]. Many diseases, however, are relevant to the abnormal blood flow in microvasculature level rather than macrovasculature level [1, 2]. While perfusion MRI can provide the microvasculature blood flow at high spatial resolution, the high cost and low mobility of MRI instrument preclude its wide use for routine disease screening. Laser Doppler can only detect blood flow at superficial tissue (∼1 mm depth) [18], thus unable to probe the diseases mostly originated from deep tissues. In summary, the routine use of the above techniques is limited due to availability, expense, and the difficulty of making continuous measurements.

In recent years, a dynamic NIRS technology, namely, diffuse correlation spectroscopy (DCS), has gained rapid development [19–23]. DCS quantifies the speckle fluctuations of diffuse light, which are sensitive to the motion of red blood cells in the microvasculature, thus permitting fast assessment of blood flow in deep tissues (up to several centimeters).

The article will briefly describe the technology and devices for NIRS and DCS, as well as their applications to physiological studies and clinic. The extension of NIRS from spectroscopy to imaging was depicted, followed by introductions of advanced algorithms that were recently proposed. The future prospective of the NIRS and DCS and their feasibilities for routine use in hospital was finally discussed.

2. NIRS Principle and Devices

When being injected into the biological tissues that are considered as turbid medium, the near-infrared light transports through the medium in “diffusive” manner [2], due to tissue absorption (mainly red blood cells) and scattering (mainly organelles and mitochondria) [1, 24]. The aim of NIRS technology is to separate the tissue absorption from scattering, through which the hemoglobin concentration is derived. Over the years, two approaches, namely, the time-resolved spectroscopy [25, 26] and frequency-modulated spectroscopy [27, 28], have been developed to account for this separation. For the time-resolved approach, a pulsed laser (a few of picoseconds) is adopted and the time-of-flight of photons is recorded by a digital time-correlated single photon counting (TCSPC) board. For frequency-modulated approach, the light is modulated at high frequency (around 100 MHz) [29], so that the phase of modulated light escaping out from the tissue could be traced.

Both these two approaches were commercialized with numerous applications [30], particularly for neurological studies. Either the time-resolved or frequency-modulated approach, however, requires complicated hardware to manipulate the light, substantially increasing the cost of instrument.

In practice, the NIRS was often simplified by using continuous-wave (CW) light instead of pulsed or frequency-modulated light. In CW oximeter, the tissue scattering (μs) was assumed to be a constant number and only the relative change of tissue absorption coefficient needs to be determined by the modified Beer–Lambert law.

where and are the light intensity detected at time t₀ and t, respectively; is the relative change of the absorption coefficient at wavelength ; and L is the mean photon length between a specific source-detector (S-D) pair. By collecting the light at two wavelengths, the oxygenation changes can be calculated as [31–33]

where and are the relative change of oxy- and deoxy-hemoglobin concentration and and are their extinction coefficients, respectively.

Figure 1 illustrates a typical scheme of CW oximeter. Briefly, a dual-wavelength LED at NIR range was controlled by transistor-transistor logic (TTL) signals, injecting the photons at each wavelength into the tissue alternately. The escaped photons were collected in parallel by all optical detectors and transferred to the computer, wherein the modified Beer–Lambert law was applied, yielding time-course oxygenation data.

Figure 1: Graphic illustration of a NIRS tissue oximeter.

The reason for using CW light as the source light of NIRS technology is that the tissue scattering (majorly contributed from organelles and mitochondria) changes much less than the tissue absorption in most of physiological status. Thus, the tissue blood oxygenation would still been properly extracted by the CW tissue oximeter. Moreover, this low-cost device could minimize the cross-talk between the absorption and scattering.

The CW tissue oximeter was also commercialized by different manufacturers (e.g., NIRO, Hamamatsu Inc.; INVOS, Somanetics Inc.) and has various clinical applications. At CW working mode, the absolute value of tissue absorption coefficient (μa) could be extracted by the near-infrared light collected from multiple source-detector (S-D) separations, according to the theory of spatial-resolved spectroscopy [34]. This theory permits generation of a slope through multiple S-D separations of optical data, from which the tissue absorption will be determined. In general, more S-D separations perform better in calculating absorption coefficient (μa) with stable and accurate values. For example, four S-D separations (2.0 to 3.5 cm) were used in a commercial tissue oximeter (OxiplexTS, ISS Inc). On the other hand, however, more S-D separation leads to higher cost of instrumentation and increases the size of optical probe. Hence, two S-D separations, i.e., the minimal number to determine a slope, were used more frequently in commercial instruments (e.g., NIRO, Hamamatsu Inc.; INVOS, Somanetics Inc.).

Despite the wide use of CW light, it is worthy to note that the tissue absorption coefficient (μa) and reduced scattering coefficient (μs′) cannot be uniquely determined from the CW tissue oximeter alone, regardless of the number of S-D separations being used [35]. To separate these two coefficients, time-resolved spectroscopy or frequency-modulated spectroscopy must be adopted.

3. DCS Principle and Devices

DCS is originated from the technology of dynamic light scattering [36, 37], in which single scattering event is allowed within thin biological sample. To account for in vivo measurement in bulk tissues or organs, advanced mathematical procedures were proposed to model the multiple scattering event, namely, diffuse correlation spectroscopy (DCS) [20, 38], or diffusing wave spectroscopy (DWS) [39, 40]. As depicted in Figure 2, a long coherence length (>5 m) laser in near-infrared range injects photons into the tissue via optical fiber. The photons travel through the tissue and are scattered back to the detector via a single mode fiber. There are fluctuations in the intensity of escaped photons, which is due to the moving of scatters (mainly red blood cells) in the microvasculature system. A digital correlator takes the detector output and uses photon arrival times to compute the light intensity autocorrelation function. From the normalized intensity autocorrelation function, the electric field temporal autocorrelation function , which satisfies the correlation diffusion equation in turbid medium [20, 38], is calculated. Approximately, the normalized electric field temporal function at early delay time decays exponentially with τ (delay time). The blood flow index (BFI) information is extracted by fitting the curve of temporal function whose decay rate depends mainly on the motion of red blood cells [20, 21, 41].

Figure 2: Graphic illustration of a tissue flowmeter.

The current optical probes for DCS/DCT are designed with source and detector fibers confined by a flat foam pad, restricting the use of this technology for tissues with different geometries. Moreover, the requirement in temporal/spatial resolution varies across different tissues and organs.

As a relative new technology, there are commercial devices for DCS (e.g., Neuro-Monitor, HemoPhotonics Inc., Spain). Similar to tissue oximeter, this commercial DCS/DCT device is primarily utilized for neurological function evaluation, limiting its application for other diseases.

4. Integrated System and Hybrid Optical Probe

Either the NIRS or DCS, when being used alone, could provide only single hemodynamic parameter (oxygenation or blood flow). At early time, the NIRS and DCS devices were applied, respectively, to the biological tissues, in order to obtain comprehensive hemodynamic information. Under this framework, however, it is difficult to measure the oxygenation and blood flow from the same tissue because the NIRS and DCS probes cannot be physically overlapped. More importantly, since the two devices are controlled independently, the source light from one device interferences with the signal collection by another device, if the two probe are placed closely on the tissue. In recent years, attempts have been made to combine the DCS with NIRS devices to form an integrated system, for simultaneous measurement of blood flow and oxygenation at microvasculature level, also permitting estimate of oxygen metabolic parameters (e.g., oxygen consumption rate) [11, 42]. Such integrated NIRS/DCS system for tissue oxygen consumption monitoring is already commercially available (MetaOx, ISS Inc., USA). Prior to the commercial products, the integrated system has been constructed by various research groups.

Figure 3 shows a typical integrated optical system combing the tissue oximeter and flowmeter. Briefly, the dual-wavelength NIR LED (oximeter) and the long-coherence NIR laser are turned on/off by the TTL switch (Figure 3(a)), launching the photon into biological tissues alternately. The photons scattered back from the tissue are collected by either optical detectors (oximeter) or single photon detectors (flowmeter). In this pattern, the two devices work on/off in sequential pattern.

Figure 3: Graphic illustration of the integrated optical system (a) and the hybrid optical probe on human lower limb (b). The application of this system to simultaneously monitor the blood flow and oxygenation changes during arterial revascularization is exhibited (c).

To comply with the integrated system, one needs a hybrid optical probe, containing all source and detector fibers in place. This configuration of fiber setup ensures the blood flow and oxygenation information measured by DCS and NIRS to cover mostly the same volume in the tissue. As mentioned earlier, the sources of DCS and NIRS are turned on/off sequentially, in order to avoid the light interference between the two technologies.

For construction of hybrid optical probe, a foam pad is used to house both NIRS and DCS fibers (Figure 3(b)) that are originally distributed in two separated pads. Within the pad, the fibers are arranged to form multiple-distance source-detector pairs for both NIRS and DCS measurements. This all-in-one probe has many advantages, such as compact size, easy to tape, and enable detection of blood flow and oxygenation in an overlapped tissue volume, as stated earlier, making it possible to estimate the metabolic parameters (e.g., the oxygen consumption rate). Moreover, the multiple separations provide hemodynamic information from different layers of tissues (e.g., scalp, skull, and brain cortex), which can be separated by using the advanced algorithms.

The design of optical probe would be adaptable to different tissues or organs. For the brain studies, there is little deformation on the head when subject to contact pressure. Hence, the materials with good elasticity and stiffness (e.g., silicon glue and polymer foam) are often used to house the S-D fibers, which would adapt the small curvature of the human head in contact pattern. For measurement of the breast tumor in which large tissue deformation occurs, noncontact pattern of measurement is preferred in order to adapt to the large curvature of breast. For such clinical applications, a lens system is designed to focus the source and detector on the tissue surface [4], permitting noncontact collection of optical signals. To minimize the instrumentation cost, the lens system will be controlled to scan over the tissue surface, greatly reducing the number of optical components (laser and detectors) [4].

The integrated NIRS system and the hybrid optical probe, depicted above, have advantages of high-throughput, portability, and easy-to-use, which have high flexibility to be used in different tissue type, disease category, and clinical setting. The successful carrying out of the hemodynamic measurements promote the translation of DCS/DCT into clinic, making it possible to routinely use this technology for early detection and therapeutic evaluation of a variety of diseases.

Table 1 lists the physiological or clinical studies wherein the integrated NIRS/DCS instrument and hybrid probe were used. Among these, a small portion of studies reports the improvements on the optical instrumentation or probe design. For example, Shang et al. [32] combined a commercial frequency-domain NIRS oximeter and a dual-length DCS flowmeter to form an integrated system. Accordingly, a hybrid probe that houses the NIRS fibers and DCS fibers in the same foam was designed. In another study, Li et al. [62] designed a noncontact probe which focuses the source and detector light on the surface of the target tissue via lens system. This probe permits simultaneous monitoring of the blood flow and oxygenation without contact with the tissues. This noncontact pattern of measurement is suitable to be used in deformable tissues such as human breast. Both contact and noncontact systems were generally validated through the manipulation protocol of temporary cuff occlusion on the limb, which reduces the blood supply to the skeletal muscle where the blood flow and oxygenation are measured.

Table 1: Studies wherein the integrated NIRS/DCS instrument and hybrid probe were utilized.

Year
Ref. no.
Authors
Target tissue
Measured variables
Physiological manipulation/diseases

2001
[20]
Cheung et al.
Rat brain
BFI, hemoglobin, StO₂
Hypercapnia

2003
[42]
Culver et al.
Rat brain
BFI, hemoglobin, StO₂ oxygen consumption rate
Ischemia

2004
[43]
Durduran et al.
Human brain
BFI, hemoglobin change, oxygen consumption rate
Cortex activation

2005
[44]
Yu et al.
Human skeletal muscle
BFI, hemoglobin, StO₂ oxygen consumption rate
Plantar flexion exercise

2006
[45]
Sunar et al.
Human tumor
BFI, hemoglobin, StO₂
Head/neck tumor

2006
[46]
Yu et al.
Human tumor
BFI, StO₂
Prostrate tumor

2007
[47]
Sunar et al.
Mouse tumor
BFI, StO₂
Melanoma tumor

2007
[48]
Zhou et al.
Human tumor
BFI, hemoglobin
Breast tumor

2009
[32]
Shang et al.
Human skeletal muscle
BFI, hemoglobin change
Cuff occlusion on upper limb

2009
[49]
Durduran et al.
Human brain
BFI, hemoglobin, StO₂
Stroke

2009
[50]
Zhou et al.
Piglet brain
BFI, hemoglobin, StO₂
Closed head injury

2010
[51]
Durduran et al.
Human infant brain
BFI, hemoglobin oxygen consumption rate
Congenital heart defects

2010
[52]
Mesquita et al.
Mouse skeletal muscle
BFI, StO₂ oxygen consumption rate
Hindlimb ischemia

2010
[13]
Roche-Labarbe et al.
Human infant brain
BFI, hemoglobin, StO₂ oxygen consumption rate
Premature

2010
[53]
Edlow et al.
Human infant brain
BFI, hemoglobin
Posture change

2010
[54]
Kim et al.
Human brain
BFI, hemoglobin change
Critically brain-injured

2011
[3]
Shang et al
Mouse brain
BFI, hemoglobin change
Ischemia

2011
[55]
Yu et al.
Human skeletal muscle
BFI, hemoglobin change
Ischemia

2011
[33]
Shang et al.
Human brain
BFI, hemoglobin change
Ischemia

2012
[56]
Cheng et al.
Human brain
BFI, hemoglobin change
Posture change

2012
[57]
Dong et al.
Human tumor
BFI, hemoglobin change
Head/neck tumor

2012
[58]
Gurley et al.
Human skeletal muscle
BFI, hemoglobin, StO₂ oxygen consumption rate
Handgrip exercise

2012
[11]
Shang et al.
Human skeletal muscle
BFI, hemoglobin, StO₂ oxygen consumption rate
Fibromyalgia

2013
[59]
Buckley et al.
Human infant brain
BFI, StO₂ oxygen consumption rate
Cardiac surgery

2013
[60]
Mesquita et al.
Human skeletal muscle
BFI, StO₂
Peripheral artery disease

2013
[61]
Mesquita et al.
Human brain
BFI, hemoglobin oxygen consumption rate
Transcranial magnetic stimulation

2013
[62]
Li et al.
Human skeletal muscle
BFI, hemoglobin change
Cuff occlusion on upper limb

2013
[63]
Shang et al.
Human skeletal muscle
BFI, hemoglobin change
Electrical stimulation

2014
[64]
Jain et al.
Human infant brain
BFI, StO₂ oxygen consumption rate
Congenital heart disease

2014
[65]
Cheng et al.
Human brain
BFI, hemoglobin change
Posture change

2014
[6]
Hou et al.
Human brain
BFI, hemoglobin change
Sleep apnea

2015
[66]
Henry et al.
Human skeletal muscle
BFI, hemoglobin, StO₂ oxygen consumption rate
Plantar flexion exercise

2016
[12]
Dong et al.
Human tumor
BFI, hemoglobin change
Head/neck tumor

2017
[67]
Baker et al.
Human skeletal muscle
BFI, hemoglobin, StO₂ oxygen consumption rate
Peripheral artery disease

2018
[68]
Ko et al.
Swine brain
BFI, hemoglobin, StO₂ oxygen consumption rate
Deep hypothermic

In addition to the instrumentation improvements as mentioned above, the majority of studies using the integrated system focused on the physiological or clinical applications. For the physiological studies, a variety of experimental manipulation were designed to challenge the target tissues, such as the hypercapnia to rat brain [20], cortex activation to human brain [43], plantar flexion or handgrip exercise to human skeletal muscles [58, 66], posture change to human brain [53, 56, 65], and the hypothermic to swine brain [68].

These physiological manipulations alter the amount of tissue hemodynamics (blood flow oxygenation and oxygen consumption rate), making it convenient to detect the potential defects or diseases (e.g., aging, autoregulation impairment, and vasovagal syncope).

For the clinical or preclinical studies, the integrated system was primarily used to probe ischemic diseases or tumor. For these purposes, animal models were often used to create the relevant diseases, including the cerebral ischemia on rat and mouse [3, 42], hindlimb ischemia on the mouse [52], and melanoma tumor on the mouse [47]. The outcomes derived from the animal studies would assist the understanding of pathophysiological mechanism and improvements of therapeutic approaches. In addition to the animal models, the integrated system was also directly applied to clinic, due to the advantages of noninvasiveness and easy-to-use. For example, Sunar et al. [47] and Dong et al. [12, 57] longitudinally monitored the blood flow and oxygenation changes in the head/neck tumor over weeks of chemoradiation therapy and found these hemodynamic changes were associated with pathophysiological outcomes. Yu et al. [46] and Zhou et al. [50] assessed the hemodynamic changes in prostate and breast tumors during photodynamic therapy and chemotherapy, respectively. Durduran et al. [49] investigated the hemodynamic responses to posture change in the stroke patients and found that the response pattern was significantly affected by the diseases. Shang et al. [33] and Yu et al. [55], on the other hand, evaluated the integrated system for intraoperative monitoring of the carotid and femoral arterial surgeries, respectively. The experiment results verified the high sensitivity of the integrated NIRS/DCS technology for ischemic detection during the vascular surgeries.

As a graphic example, Figure 3(c) shows the muscle hemodynamics (blood flow and oxygenation) during arterial vascularization [55]. It is clearly seen that the arterial clamping causes acute ischemia and hypoxia on the muscle tissue, i.e., decrease in blood flow/Δ[HbO₂] and increase in Δ[Hb]. Following the release of cuff occlusion, there is hyperemic response, elevating the blood flow to a peak value. Thereafter, both the blood flow and oxygenation recover gradually towards their baseline. Additionally, the postrevascularization improvements in tissue hemodynamics were detected. This study demonstrates the high sensitivity of NIRS/DCS in monitoring the intraoperative interventions as well as the acute therapeutic assessments.

In addition, there are reports of hemodynamic assessment on other diseases, including the congenital heart defects [51, 64], premature [13], critically brain-injured [54], fibromyalgia [11], and peripheral artery disease [60, 67]. Those clinical studies validated the wide usage of integrated NIRS/DCS system for various diseases.

5. Imaging Modalities for NIRS and DCS

The NIRS and DCS provide time variation of oxygenation and blood flow, respectively. In order to obtain the spatial contrast of the tissue hemodynamics, there is a need to extend the NIRS and DCS from spectroscopy to tomography, termed as diffuse optical tomography (DOT) and diffuse correlation tomography (DCT), respectively. The DOT/DCT share the same physical principle of NIRS/DCS, except that far more S-D pairs are needed to perform image reconstruction. As reported in literature, a high-density DOT system consists of 96 sources and 92 detectors, providing the optical measurements from more 1200 S-D pairs [69]. Similarly, a representative DCT device consists of 2 laser diodes and 15 APD and with 15 steps of line scanning, generating a total of 450 autocorrelation curves [70]. In addition to the hardware configuration, the algorithms adopted for image reconstruction is critical in determining imaging quality of blood flow or oxygenation, which will be discussed in Section 6. Apparently, a large number of optical components (laser, detector, etc.) are required for DOT or DCT, substantially elevating the cost of instrumentation. A few efforts have been made to lower down the physical S-D measurements. For example, the line or rotated scanning could efficiently reduce the number of optical components but needs the optical probe to operate in noncontact pattern. The noncontact measurement, however, is susceptible to the motion artifacts and ambient light. The optical switch, which enables sharing the laser or detector among multilocation measurements, could be a good alternative option to save the instrument cost. Nevertheless, a number of hardware source and detector pairs are still needed for the imaging of tissue hemodynamics, making the instrument no longer portable. Thus, the hemodynamic imaging is not the focus of this article.

6. Advanced Algorithms for NIRS and DCS

Apart from instrumentation, the reconstruction algorithm of DCS/NIRS also contributes greatly to the accuracy and robustness of oxygenation and blood flow values, which is highly crucial to the detective sensitivity of diagnostic outcomes [31, 71]. For DCS, as stated earlier, the conventional algorithm is based on the analytical solution of a partial differential equation (PDE) [20, 21, 41]. However, the analytical solution of PDE is dependent on tissue geometry and usually has complicated mathematical forms. For the purpose of easy implementation, the tissue geometry is often assumed as semi-infinite. There are also analytical solutions to PDE for other regular tissue geometries, such as cylinder and sphere [41, 72]. However, those solutions are rarely adopted for practical applications, due to the complicated and tedious computations. In recent years, a finite element method (FEM) for DCT was developed based on similar approach for NIRS [70]. Although the FEM was proved to be efficient in reconstruction of flow value in tissue with complicated geometry, it is difficult to fully utilize autocorrelation data, thus susceptible to data noise.

For sufficient use of the DCS/DCT data, a mathematical approach, namely, the Nth-order linear algorithm (i.e., NL algorithm), was developed for extracting BFI value. Unlike the analytical solution or FE method, the NL algorithm does not seek for PDE solution. Instead, it was originated from a combination of a Nth-order Taylor polynomial and the integral form of . After a series of mathematical procedures, the BFI is extracted by iteration of linear regressions, ultimately reaching the following expressions [73]:

Equations (3) and (4) are called first-order and Nth-order linear (NL) algorithm, respectively. Here, is the BFI of the ith tissue component; is the normalized weight of photon package collected by jth detector; and is the path length of pth photon package traveled within ith tissue component and collected by jth detector. The coefficient of matrix has the following expression:

The key parameter contained in equations (3)–(5) is , which is determined by the tissue geometry and heterogeneity, indicating that both of these two features are fully taken into consideration in the NL algorithm. The accuracy of the NL algorithm on various geometries and tissue heterogeneity was validated in DCS, through computer simulations as well as animal experiment [73, 74].

The similar procedure for heterogeneous tissue is also applicable to NIRS, leading to the following equations to calculate the changes of chromophore concentration [31, 75]:

where ODi is the optical density from ith S-D measurement; sij is the photon path length within jth tissue measured from ith S-D pair; M and N are the numbers of tissue and S-D pair, respectively; and Δμa,j is the relative change of the absorption coefficient at jth tissue. For convenience in matrix operation, equation (6) can be abbreviated as

Here, s is the matrix of photon path length. When M = N, the changes of absorption coefficient can be calculated according to the following equation:

It is possible to extract the hemodynamic information at higher accuracy by using the optical signals from more measurements (i.e., N > M). In this situation, the least-square fitting method can be used to calculate the changes of absorption coefficient at different tissues:

Due to the correlation of optical signals among different S-D pairs, the matrix s is often ill-posed, i.e., susceptible to the data noise. The noise influence on oxygenation measurement could be partially alleviated by introducing regularization approaches. One of the widely-used approaches is Tikhonov regularization, with which the final form of hemoglobin calculation is the following:

Here, α is the constant number controlling the balance between the equation fidelity and solution robustness.

7. Clinical Applications

Over the years, the NIRS and DCS for tissue oxygenation and blood flow measurements have been extensively validated to other modalities and been translational to various physiological and clinical studies, such as assessment of autoregulation capability [56], therapeutic monitoring of ischemic stroke [33], and early detection of tumors [4]. These applications are summarized in Table 1 and briefly reviewed in Section 4. Some representative applications are introduced with details in the following subsections.

7.1. Assessment of Autoregulation Capability

Studies have shown that many diseases, such as ischemic stroke, obstructive sleep apnea, vasovagal syncope, or Alzheimer’s disease, may cause deficit in autoregulation capability, which could be assessed by alteration of both blood pressure and cerebral hemodynamics. For example, change of posture could alter both blood pressure and cerebral hemodynamics. However, the changing pattern of cerebral hemodynamics in patients, dependent on the autoregulation capability, may differ significantly from those in healthy subjects, from which we can detect the neurological deficit. The technologies of NIRS and DCS were well established to measure the cerebral hemodynamics, and the blood pressure could be longitudinally assessed by a finger plethysmograph, at the pulse frequency (∼1 Hz). The physiological protocols to induce the cerebral hemodynamics include enforced breathing at 0.1 Hz, cuff occlusion on both side of lower limbs, posture transition from squatting to standing, and the head-up-tilting (HUT). Among these, the 70-degree HUT was adopted widely in physiological studies. For most of the people, it was found that there are spontaneous low frequency oscillations (LFOs) around 0.1 Hz in some physiological parameters including the blood pressure and hemodynamics, and the LFOs could be affected by the impairment in autoregulation capability. For rapid and noninvasive assessment of the autoregulation capability, a hybrid near-infrared diffuse optical instrument was constructed to simultaneously detect LFOs of cerebral blood flow (CBF) and oxygenation ([HbO₂], [Hb], and THC) [56], and a finger plethysmograph was used to detect mean arterial pressure (MAP). The coherence outcomes of LFOs between MAP and each of the hemodynamic variables show that the CBF, [HbO₂], and THC were reliable hemodynamic parameters in detecting LFOs [56], which has clinical implication in early diagnosis of cerebral impairment caused by ischemic stroke, obstructive sleep apnea, vasovagal syncope, or Alzheimer’s disease.

7.2. Therapeutic Monitoring of Ischemic Stroke

During many head and neck surgeries, such as carotid endarterectomy, the major arteries, such as common carotid artery, internal carotid artery, and middle cerebral artery, are temporarily occluded for blooding prevention and operating convenience. On the other hand, this procedure increases the risk of ischemic stroke. Therefore, real-time monitoring of cerebral hemodynamics is critical for minimizing the intraoperative stroke and acute assessment of surgical outcomes. For this purpose, the DCS and NIRS were combined in a portable device, namely, DCS flow-oximeter [32]. As a pilot study of intraoperative assessment, two fiber-optic probes were taped on both sides of forehead for cerebral hemodynamic measurements [33]. Meanwhile, the brain wave was monitored by electroencephalogram (EEG). The cerebral blood flow and oxygenation were simultaneously monitored on twelve patients throughout entire period of carotid endarterectomy. The results verified the feasibility and sensitivity of NIRS and DCS in probing transient cerebral ischemia due to arterial clamping and postsurgical hyper perfusion [33]. More importantly, the comparison outcomes show that the CBF responses to arterial clamping, measured by DCS, were faster, larger, and more sensitive than EEG responses. This study suggests the potential of NIRS/DCS technologies for routine use in surgical rooms.

7.3. Early Detection of Breast Tumors

The blood flow, which reflects the dynamic demands of oxygen in tissue, is highly sensitive to the growth of the tumor with extraordinary high metabolic rate. Hence, monitoring of microvasculature blood flow, along with oxygenation measurement, will probe the malignant tumor earlier than the morphological imaging modality such as ultrasound, CT, or MRI. The portable NIRS technology has been used to early diagnosis a few cancers, majorly, the breast tumors.

In an earlier study, a hand-held optical probe, containing four S-D pair fibers, was used to scan over the breast horizontally and vertically [5]. The optical data were transferred into a DCS device for extracting the blood flow index (BFI) at different locations. In the women with breast cancer, a high contrast in BFI was found between malignant tumor and the normal tissue, demonstrating the high sensitivity of DCS in detecting the abnormal blood flow caused by the cancer. However, the breast tissue is soft and susceptible to contact pressure. Hence, it is difficult to obtain the true BFI values in natural posturer. Recently, a noncontact probe is constructed, through which the light emitted into (or escaped out from) the tissue is delivered by the lens system [76], rather than directly contacting with the skin. The noncontact pattern of measurement efficiently avoids the breast deformation caused by the contact pressure. Moreover, the motor-controlled line scanning system permits data acquisition from a large number of S-D pairs, making it possible to reconstruct the breast tumor imaging [70]. Similar to the contact-measurement study mentioned above, the imaging from the women with malignant tumor exhibits high spatial contrast (5.9- and 10.0-fold) between the malignant tumor and normal tissue [4].

These studies indicate the potential of NIRS/DCS systems for noninvasive and early diagnosis of the breast tumors, which is greatly beneficial to the healthcare of middle and aging women.

In addition to the clinical applications mentioned above, the NIRS/DCS has been widely translated to assess the tissue hemodynamic status for other diseases or medical interferences, including the peripheral artery disease (PAD) [60], fibromyalgia [11], radiation therapy on head/neck tumor [12], and tissue transfer flaps [77]. Due to the length limit, the readers are encouraged to see references for more details.

8. Challenges and Future Prospective

Although some NIRS commercial products have been utilized for physiological and disease studies, there are some challenges for its wide translation to clinic. For example, the absolute value of oxygenation parameters ([HbO₂], [Hb], THC, and StO₂) varies over a large range in healthy population, making it difficult to set a threshold to separate the diseased tissues and normal tissues. In practice, therefore, the oxygenation contrast, either temporally or spatially, was often used for abnormal detection. Likewise, the standard procedure to detect abnormal tissue with DCS flow measurement has not been established. Besides, a few factors, such as the ambient noise, the curvature of tissue geometry, and the contact pressure between optical probe and tissue surface, will affect the NIRS/DCS technology for hemodynamic measurement.

In summary, this article provides a review of the static and dynamic NIRS technologies in comprehensive aspects, including the physical principle, the instrumentation, the algorithms for hemodynamic calculation, and the applications in physiological and clinical measurement. Over four decade development, the NIRS technology has demonstrated excellent feasibility in exploring the microvasculature oxygen kinetics and pathophysiological mechanism. The success of hemodynamic measurement greatly promoted the translation of NIRS technology to clinic. However, the accuracy and efficiency of this technology for early disease detection and therapeutic evaluation is heavily dependent on the tissue type, disease category, and clinical environment. Thus, there are still urgent needs to develop utility-driven strategy for routine application of NIRS/DCS in clinic. Ultimately, we anticipate that the portable, inexpensive NIRS/DCS systems will be widely commercialized for routine utilization in clinical rooms and patient’s bedside.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (61771433 and 61671413), National Key Research and Development Program of China (2016YFC0101605), National Key Scientific Instrument and Equipment Development Project of China (2014YQ24044508), and OIT Program of Shanxi Province.

References

Y. Shang, T. Li, and G. Yu, “Clinical applications of near-infrared diffuse correlation spectroscopy and tomography for tissue blood flow monitoring and imaging,” Physiological Measurement, vol. 38, no. 4, pp. R1–R26, 2017.

T. Durduran, R. Choe, W. B. Baker, and A. G. Yodh, “Diffuse optics for tissue monitoring and tomography,” Reports on Progress in Physics, vol. 73, no. 7, Article ID 076701, 2010.

Y. Shang, L. Chen, M. Toborek, and G. Yu, “Diffuse optical monitoring of repeated cerebral ischemia in mice,” Optics Express, vol. 19, no. 21, pp. 20301–20315, 2011.

L. He, Y. Lin, C. Huang, D. Irwin, M. M. Szabunio, and G. Yu, “Noncontact diffuse correlation tomography of human breast tumor,” Journal of Biomedical Optics, vol. 20, no. 8, Article ID 086003, 2015.

T. Durduran, R. Choe, G. Yu et al., “Diffuse optical measurement of blood flow in breast tumors,” Optics Letters, vol. 30, no. 21, pp. 2915–2917, 2005.

Y. Hou, Y. Shang, R. Cheng et al., “Obstructive sleep apnea-hypopnea results in significant variations in cerebral hemodynamics detected by diffuse optical spectroscopies,” Physiological Measurement, vol. 35, no. 10, pp. 2135–2148, 2014.

N. Khan, B. B. Williams, H. Hou, H. Li, and H. M. Swartz, “Repetitive tissue pO2 measurements by electron paramagnetic resonance oximetry: current status and future potential for experimental and clinical studies,” Antioxidants & Redox Signaling, vol. 9, no. 8, pp. 1169–1182, 2007.

F. Jobsis, “Noninvasive, infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters,” Science, vol. 198, no. 4323, pp. 1264–1267, 1977.

M. Wolf, M. Ferrari, and V. Quaresima, “Progress of near-infrared spectroscopy and topography for brain and muscle clinical applications,” Journal of Biomedical Optics, vol. 12, no. 6, Article ID 062104, 2007.

M. Ferrari, M. Muthalib, and V. Quaresima, “The use of near-infrared spectroscopy in understanding skeletal muscle physiology: recent developments,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 369, no. 1955, pp. 4577–4590, 2011.

Y. Shang, K. Gurley, B. Symons et al., “Noninvasive optical characterization of muscle blood flow, oxygenation, and metabolism in women with fibromyalgia,” Arthritis Research & Therapy, vol. 14, no. 6, p. R236, 2012.

L. X. Dong, M. Kudrimoti, D. Irwin et al., “Diffuse optical measurements of head and neck tumor hemodynamics for early prediction of chemoradiation therapy outcomes,” Journal of Biomedical Optics, vol. 21, no. 8, Article ID 085004, 2016.

N. Roche-Labarbe, S. A. Carp, A. Surova et al., “Noninvasive optical measures of CBV, StO2, CBF index, and rCMRO2 in human premature neonates' brains in the first six weeks of life,” Human Brain Mapping, vol. 31, no. 3, pp. 341–352, 2009.

M. E. Tschakovsky, J. K. Shoemaker, and R. L. Hughson, “Beat-by-beat forearm blood flow with Doppler ultrasound and strain-gauge plethysmography,” Journal of Applied Physiology, vol. 79, no. 3, pp. 713–719, 1995.

J. L. Saumet, D. L. Kellogg, W. F. Taylor, and J. M. Johnson, “Cutaneous laser-Doppler flowmetry: influence of underlying muscle blood flow,” Journal of Applied Physiology, vol. 65, no. 1, pp. 478–481, 1988.

G. Yu, T. F. Floyd, T. Durduran et al., “Validation of diffuse correlation spectroscopy for muscle blood flow with concurrent arterial spin labeled perfusion MRI,” Optics Express, vol. 15, no. 3, pp. 1064–1075, 2007.

C. Jansen, E. M. Vriens, B. C. Eikelboom, F. E. Vermeulen, J. van Gijn, and R. G. Ackerstaff, “Carotid endarterectomy with transcranial Doppler and electroencephalographic monitoring. A prospective study in 130 operations,” Stroke, vol. 24, no. 5, pp. 665–669, 1993.

D. T. Ubbink, M. J. H. M. Jacobs, G. J. Tangelder, D. W. Slaaf, and R. S. Reneman, “The usefulness of capillary microscopy, transcutaneous oximetry and laser Doppler fluxmetry in the assessment of the severity of lower limb ischaemia,” International Journal of Microcirculation, vol. 14, no. 1-2, pp. 34–44, 2008.

G. Maret and P. E. Wolf, “Multiple light scattering from disordered media. The effect of brownian motion of scatterers,” Zeitschrift für Physik B Condensed Matter, vol. 65, no. 4, pp. 409–413, 1987.

C. Cheung, J. P. Culver, K. Takahashi, J. H. Greenberg, and A. G. Yodh, “In vivocerebrovascular measurement combining diffuse near-infrared absorption and correlation spectroscopies,” Physics in Medicine and Biology, vol. 46, no. 8, pp. 2053–2065, 2001.

D. Irwin, L. Dong, Y. Shang et al., “Influences of tissue absorption and scattering on diffuse correlation spectroscopy blood flow measurements,” Biomedical Optics Express, vol. 2, no. 7, pp. 1969–1985, 2011.

D. R. Busch, J. M. Lynch, M. E. Winters et al., “Cerebral blood flow response to hypercapnia in children with obstructive sleep apnea syndrome,” Sleep, vol. 39, no. 1, pp. 209–216, 2016.

L. He, Y. Lin, Y. Shang, B. J. Shelton, and G. Yu, “Using optical fibers with different modes to improve the signal-to-noise ratio of diffuse correlation spectroscopy flow-oximeter measurements,” Journal of Biomedical Optics, vol. 18, no. 3, Article ID 037001, 2013.

G. Yu, “Near-infrared diffuse correlation spectroscopy in cancer diagnosis and therapy monitoring,” Journal of Biomedical Optics, vol. 17, no. 1, Article ID 010901, 2012.

X. Intes, J. Chen, V. Venugopal, and F. Lesage, “Time-resolved diffuse optical tomography with patterned-light illumination and detection,” Optics Letters, vol. 35, no. 13, pp. 2121–2123, 2010.

B. Chance, J. S. Leigh, H. Miyake et al., “Comparison of time-resolved and -unresolved measurements of deoxyhemoglobin in brain,” Proceedings of the National Academy of Sciences, vol. 85, no. 14, pp. 4971–4975, 1988.

M. Dehaes, A. Aggarwal, P.-Y. Lin et al., “Cerebral oxygen metabolism in neonatal hypoxic ischemic encephalopathy during and after therapeutic hypothermia,” Journal of Cerebral Blood Flow & Metabolism, vol. 34, no. 1, pp. 87–94, 2013.

J. H. Choi, S. Lee, M. Lee, D. Koh, and B. M. Kim, “Cerebral hemodynamic responses to seizure in the mouse brain: simultaneous near-infrared spectroscopy-electroencephalography study,” Journal of Biomedical Optics, vol. 15, no. 3, Article ID 037010, 2010.

G. Yu, T. Durduran, D. Furuya, J. H. Greenberg, and A. G. Yodh, “Frequency-domain multiplexing system for in vivo diffuse light measurements of rapid cerebral hemodynamics,” Applied Optics, vol. 42, no. 16, pp. 2931–2939, 2003.

K. Yoshitani, M. Kawaguchi, K. Tatsumi, K. Kitaguchi, and H. Furuya, “A comparison of the INVOS 4100 and the NIRO 300 near-infrared spectrophotometers,” Anesthesia & Analgesia, vol. 94, no. 3, pp. 586–590, 2002.

G. Strangman, M. A. Franceschini, and D. A. Boas, “Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters,” Neuroimage, vol. 18, no. 4, pp. 865–879, 2003.

Y. Shang, Y. Zhao, R. Cheng, L. Dong, D. Irwin, and G. Yu, “Portable optical tissue flow oximeter based on diffuse correlation spectroscopy,” Optics Letters, vol. 34, no. 22, pp. 3556–3558, 2009.

Y. Shang, R. Cheng, L. Dong, S. J. Ryan, S. P. Saha, and G. Yu, “Cerebral monitoring during carotid endarterectomy using near-infrared diffuse optical spectroscopies and electroencephalogram,” Physics in Medicine and Biology, vol. 56, no. 10, pp. 3015–3032, 2011.

H. Liu, D. A. Boas, Y. Zhang, A. G. Yodh, and B. Chance, “Determination of optical properties and blood oxygenation in tissue using continuous NIR light,” Physics in Medicine and Biology, vol. 40, no. 11, pp. 1983–1993, 1999.

S. R. Arridge and W. R. B. Lionheart, “Nonuniqueness in diffusion-based optical tomography,” Optics Letters, vol. 23, no. 11, pp. 882–884, 1998.

W. Brown, Dynamic Light Scattering: The Method and Some Applications, Clarendon, New York, NY, USA, 1993.

A. B. Leung, K. I. Suh, and R. R. Ansari, “Particle-size and velocity measurements in flowing conditions using dynamic light scattering,” Applied Optics, vol. 45, no. 10, pp. 2186–2190, 2006.

D. A. Boas, L. E. Campbell, and A. G. Yodh, “Scattering and imaging with diffusing temporal field correlations,” Physical Review Letters, vol. 75, no. 9, pp. 1855–1858, 1995.

D. J. Pine, D. A. Weitz, P. M. Chaikin, and E. Herbolzheimer, “Diffusing wave spectroscopy,” Physical Review Letters, vol. 60, no. 12, pp. 1134–1137, 1988.

J. Li, M. Ninck, L. Koban, T. Elbert, J. Kissler, and T. Gisler, “Transient functional blood flow change in the human brain measured noninvasively by diffusing-wave spectroscopy,” Optics Letters, vol. 33, no. 19, pp. 2233–2235, 2008.

D. A. Boas and A. G. Yodh, “Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation,” Journal of the Optical Society of America A, vol. 14, no. 1, pp. 192–215, 1997.

J. P. Culver, T. Durduran, D. Furuya, C. Cheung, J. H. Greenberg, and A. G. Yodh, “Diffuse optical tomography of cerebral blood flow, oxygenation, and metabolism in rat during focal ischemia,” Journal of Cerebral Blood Flow & Metabolism, vol. 23, no. 8, pp. 911–924, 2016.

T. Durduran, G. Yu, M. G. Burnett et al., “Diffuse optical measurement of blood flow, blood oxygenation, and metabolism in a human brain during sensorimotor cortex activation,” Optics Letters, vol. 29, no. 15, pp. 1766–1768, 2004.

G. Yu, T. Durduran, G. Lech et al., “Time-dependent blood flow and oxygenation in human skeletal muscles measured with noninvasive near-infrared diffuse optical spectroscopies,” Journal of Biomedical Optics, vol. 10, no. 2, Article ID 024027, 2005.

U. Sunar, H. Quon, T. Durduran et al., “Noninvasive diffuse optical measurement of blood flow and blood oxygenation for monitoring radiation therapy in patients with head and neck tumors: a pilot study,” Journal of Biomedical Optics, vol. 11, no. 6, Article ID 064021, 2006.

G. Yu, T. Durduran, C. Zhou et al., “Real-time in situ monitoring of human prostate photodynamic therapy with diffuse light,” Photochemistry and Photobiology, vol. 82, no. 5, pp. 1279–1284, 2006.

U. Sunar, S. Makonnen, C. Zhou et al., “Hemodynamic responses to antivascular therapy and ionizing radiation assessed by diffuse optical spectroscopies,” Optics Express, vol. 15, no. 23, pp. 15507–15516, 2007.

C. Zhou, R. Choe, N. Shah et al., “Diffuse optical monitoring of blood flow and oxygenation in human breast cancer during early stages of neoadjuvant chemotherapy,” Journal of Biomedical Optics, vol. 12, no. 5, Article ID 051903, 2007.

T. Durduran, C. Zhou, B. L. Edlow et al., “Transcranial optical monitoring of cerebrovascular hemodynamics in acute stroke patients,” Optics Express, vol. 17, no. 5, pp. 3884–3902, 2009.

C. Zhou, S. A. Eucker, T. Durduran et al., “Diffuse optical monitoring of hemodynamic changes in piglet brain with closed head injury,” Journal of Biomedical Optics, vol. 14, no. 3, Article ID 034015, 2009.

T. Durduran, C. Zhou, E. M. Buckley et al., “Optical measurement of cerebral hemodynamics and oxygen metabolism in neonates with congenital heart defects,” Journal of Biomedical Optics, vol. 15, no. 3, Article ID 037004, 2010.

R. C. Mesquita, N. Skuli, M. N. Kim et al., “Hemodynamic and metabolic diffuse optical monitoring in a mouse model of hindlimb ischemia,” Biomedical Optics Express, vol. 1, no. 4, pp. 1173–1187, 2010.

B. L. Edlow, M. N. Kim, T. Durduran et al., “The effects of healthy aging on cerebral hemodynamic responses to posture change,” Physiological Measurement, vol. 31, no. 4, pp. 477–495, 2010.

M. N. Kim, T. Durduran, S. Frangos et al., “Noninvasive measurement of cerebral blood flow and blood oxygenation using near-infrared and diffuse correlation spectroscopies in critically brain-Injured adults,” Neurocritical Care, vol. 12, no. 2, pp. 173–180, 2009.

G. Yu, Y. Shang, Y. Zhao, R. Cheng, L. Dong, and S. P. Saha, “Intraoperative evaluation of revascularization effect on ischemic muscle hemodynamics using near-infrared diffuse optical spectroscopies,” Journal of Biomedical Optics, vol. 16, no. 2, Article ID 027004, 2011.

R. Cheng, Y. Shang, D. Hayes, S. P. Saha, and G. Yu, “Noninvasive optical evaluation of spontaneous low frequency oscillations in cerebral hemodynamics,” Neuroimage, vol. 62, no. 3, pp. 1445–1454, 2012.

L. Dong, M. Kudrimoti, R. Cheng et al., “Noninvasive diffuse optical monitoring of head and neck tumor blood flow and oxygenation during radiation delivery,” Biomedical Optics Express, vol. 3, no. 2, pp. 259–272, 2012.

K. Gurley, Y. Shang, and G. Yu, “Noninvasive optical quantification of absolute blood flow, blood oxygenation, and oxygen consumption rate in exercising skeletal muscle,” Journal of Biomedical Optics, vol. 17, no. 7, Article ID 075010, 2012.

E. M. Buckley, J. M. Lynch, D. A. Goff et al., “Early postoperative changes in cerebral oxygen metabolism following neonatal cardiac surgery: effects of surgical duration,” Journal of Thoracic and Cardiovascular Surgery, vol. 145, no. 1, pp. 196–205, 2013.

R. C. Mesquita, M. Putt, M. Chandra et al., “Diffuse optical characterization of an exercising patient group with peripheral artery disease,” Journal of Biomedical Optics, vol. 18, no. 5, Article ID 057007, 2013.

R. C. Mesquita, O. K. Faseyitan, P. E. Turkeltaub et al., “Blood flow and oxygenation changes due to low-frequency repetitive transcranial magnetic stimulation of the cerebral cortex,” Journal of Biomedical Optics, vol. 18, no. 6, Article ID 067006, 2013.

T. Li, Y. Lin, Y. Shang et al., “Simultaneous measurement of deep tissue blood flow and oxygenation using noncontact diffuse correlation spectroscopy flow-oximeter,” Scientific Reports, vol. 3, no. 1, p. 1358, 2013.

Y. Shang, Y. Lin, B. A. Henry et al., “Noninvasive evaluation of electrical stimulation impacts on muscle hemodynamics via integrating diffuse optical spectroscopies with muscle stimulator,” Journal of Biomedical Optics, vol. 18, no. 10, Article ID 105002, 2013.

V. Jain, E. M. Buckley, D. J. Licht et al., “Cerebral oxygen metabolism in neonates with congenital heart disease quantified by MRI and optics,” Journal of Cerebral Blood Flow & Metabolism, vol. 34, no. 3, pp. 380–388, 2013.

R. Cheng, Y. Shang, S. Wang et al., “Near-infrared diffuse optical monitoring of cerebral blood flow and oxygenation for the prediction of vasovagal syncope,” Journal of Biomedical Optics, vol. 19, no. 1, Article ID 017001, 2014.

B. Henry, M. Zhao, Y. Shang et al., “Hybrid diffuse optical techniques for continuous hemodynamic measurement in gastrocnemius during plantar flexion exercise,” Journal of Biomedical Optics, vol. 20, no. 12, Article ID 125006, 2015.

W. B. Baker, Z. Li, S. S. Schenkel et al., “Effects of exercise training on calf muscle oxygen extraction and blood flow in patients with peripheral artery disease,” Journal of Applied Physiology, vol. 123, no. 6, pp. 1599–1609, 2017.

T. S. Ko, C. D. Mavroudis, W. B. Baker et al., “Non-invasive optical neuromonitoring of the temperature-dependence of cerebral oxygen metabolism during deep hypothermic cardiopulmonary bypass in neonatal swine,” Journal of Cerebral Blood Flow and Metabolism, Article ID 271678X18809828, 2018, In press.

A. T. Eggebrecht, S. L. Ferradal, A. Robichaux-Viehoever et al., “Mapping distributed brain function and networks with diffuse optical tomography,” Nature Photonics, vol. 8, no. 6, pp. 448–454, 2014.

Y. Lin, C. Huang, D. Irwin, L. He, Y. Shang, and G. Yu, “Three-dimensional flow contrast imaging of deep tissue using noncontact diffuse correlation tomography,” Applied Physics Letters, vol. 104, no. 12, Article ID 121103, 2014.

Y. Shang, “Advances in reconstruction algorithms for diffuse correlation spectroscopy and tomography,” in International Conference on Sensing and Imaging 2017, M. Jiang, N. Ida, A. K. Louis, and E. T. Quinto, Eds., pp. 15–28, Springer, Chengdu, China, 2017, Lecture Notes in Electrical Engineering 506.

D. Boas, Diffuse Photon Probes of Structural and Dynamical Properties of Turbid Media: Theory and Biomedical Applications, Department of Physics, University of Pennsylvania, Philadelphia, PA, USA, 1996.

Y. Shang and G. Yu, “A Nth-order linear algorithm for extracting diffuse correlation spectroscopy blood flow indices in heterogeneous tissues,” Applied Physics Letters, vol. 105, no. 13, Article ID 133702, 2014.

Y. Shang, T. Li, L. Chen, Y. Lin, M. Toborek, and G. Yu, “Extraction of diffuse correlation spectroscopy flow index by integration of Nth-order linear model with Monte Carlo simulation,” Applied Physics Letters, vol. 104, no. 19, Article ID 193703, 2014.

A. E. Cerussi, D. Jakubowski, N. Shah et al., “Spectroscopy enhances the information content of optical mammography,” Journal of Biomedical Optics, vol. 7, no. 1, pp. 60–71, 2002.

Y. Lin, L. He, Y. Shang, and G. Yu, “Noncontact diffuse correlation spectroscopy for noninvasive deep tissue blood flow measurement,” Journal of Biomedical Optics, vol. 17, no. 1, Article ID 010502, 2012.

C. Huang, J. P. Radabaugh, R. K. Aouad et al., “Noncontact diffuse optical assessment of blood flow changes in head and neck free tissue transfer flaps,” Journal of Biomedical Optics, vol. 20, no. 7, Article ID 075008, 2015.