15 June 2016 Reproducibility of parameters of postocclusive reactive hyperemia measured by diffuse optical tomography
Author Affiliations +
J. of Biomedical Optics, 21(6), 066012 (2016). doi:10.1117/1.JBO.21.6.066012
The application of near-infrared spectroscopy (NIRS) to assess microvascular function has shown promising results. An important limitation when using a single source-detector pair, however, is the lack of depth sensitivity. Diffuse optical tomography (DOT) overcomes this limitation using an array of sources and detectors that allow the reconstruction of volumetric hemodynamic changes. This study compares the key parameters of postocclusive reactive hyperemia measured in the forearm using standard NIRS and DOT. We show that while the mean parameter values are similar for the two techniques, DOT achieves much better reproducibility, as measured by the intraclass correlation coefficient (ICC). We show that DOT achieves high reproducibility for muscle oxygen consumption (ICC: 0.99), time to maximal HbO2 (ICC: 0.94), maximal HbO2 (ICC: 0.99), and time to maximal HbT (ICC: 0.99). Absolute reproducibility as measured by the standard error of measurement is consistently smaller and close to zero (ideal value) across all parameters measured by DOT compared to NIRS. We conclude that DOT provides a more robust characterization of the reactive hyperemic response and show how the availability of volumetric hemodynamic changes allows the identification of areas of temporal consistency, which could help characterize more precisely the microvasculature.
Vidal-Rosas, Billings, Chico, and Coca: Reproducibility of parameters of postocclusive reactive hyperemia measured by diffuse optical tomography



Atherosclerosis and resulting cardiovascular diseases such as stroke and myocardial infarction are a major cause of death in developed countries. These account for more than 32% of mortality worldwide,1 and in England and Wales cardiovascular disease was responsible for almost 30% of deaths in 2011.2 Many noninvasive methods have been developed to assess the peripheral vascular system and identify signs of atherosclerosis at an early stage. Optical methods, in particular, have received attention due to their capability of measuring tissue oxygenation and blood perfusion3,4 and because they offer several attractive features, such as portability, compactness, fast data acquisition, and noninvasiveness.

Near-infrared spectroscopy (NIRS) can determine changes in tissue hemodynamics and oxygenation5,6 by measuring tissue absorbance at several wavelengths in the near-infrared range of the electromagnetic spectrum (650 to 950 nm). Typically, NIRS employs a few source-detector pairs to carry out measurements. As a consequence, the spatial resolution of NIRS, which is dictated by the optode separation, is relatively low.7 Diffuse optical tomography (DOT) overcomes this limitation by employing a larger number of sources and detectors to enable three-dimensional (3-D) volumetric reconstruction8 of changes in tissue hemodynamics.

NIRS and DOT have been used in a wide range of applications including functional imaging of the brain,9 assessment of muscle oxygenation,10 and cancer detection. It is widely accepted that DOT outperforms NIRS.8,11 For example, localized changes in hemodynamics due to motor tasks in adults and motor-sensory brain activation in neonates were accurately measured using DOT, while NIRS could not even detect relative changes11 because of low spatial sampling. In another study, DOT could discriminate the somatosensory activation of two fingers, while in the same experiment, a 12-channel NIRS setting failed to resolve the activation.12 Among the factors affecting the accuracy of NIRS concentration calculations are the differences in the pathlength factor: location, spatial extent, and heterogeneous distribution of absorption changes, e.g., multiple absorption foci. Although these sources of error could be minimized,13 DOT accounts for these problems implicitly.7

DOT requires solving two distinct problems: the forward problem and the inverse problem. The forward problem requires solving of the equation governing the photon transport in tissue to predict the detector measurements; typically this involves solving a diffusion equation over a 3-D spatial domain using the finite element method. The inverse problem involves estimating the optical properties of tissue to minimize the difference between experimental and model-predicted measurements; this is achieved by solving a nonlinear optimization problem, which can take hours to complete even on a high-end workstation. However, as we have recently demonstrated,14 it is possible to significantly speed up reconstruction of hemodynamic changes in complex tissue structures by using reduced order models of photon transport in tissue. This makes it possible to perform real-time monitoring of hemodynamic responses even with relatively modest computing resources.

NIRS, like other methods of measuring microvascular function, can detect significant differences in microvascular function between groups of healthy controls and patients with peripheral vascular diseases,1516.17 peripheral artery disease (PAD), and coronary heart disease.1718.19 However, these techniques are highly variable and have not been shown to be helpful in predicting an individual’s risk of future cardiovascular disease. Therefore, from a clinical point of view, refinements that improve reproducibility and reduce variability are highly desirable.

Typically, the evaluation of microvascular function by NIRS relies solely on a few measuring channels, and although repeatability and accuracy of results are promising,18,20,21 there has not been a direct comparison of the parameters obtained during and after arterial occlusion by NIRS versus DOT. Furthermore, while multichannel systems for analysis of microvascular function are available,2223.24.25 studies on repeatability are lacking.

This study aims to evaluate and compare, using experimental data collected from healthy volunteers, the intrasubject reproducibility of key parameters of the hemodynamic response during postocclusive reactive hyperemia, obtained using DOT and NIRS, and to highlight the potential advantages of DOT in assessing endothelial function.





This study was approved by the Research Ethics Committee of the University of Sheffield. A group of 17 subjects was recruited, after giving informed consent. The group consisted of 11 men and 6 women. Baseline characteristics are listed in Table 1. Smokers or those with a history of cardiovascular disease were excluded.

Table 1

Baseline characteristic of study subjects (n=17).

Age (year)33±7.5
Weight (kg)70.5±14.1
Height (cm)172.1±8.3
BMI (kg/m2)23.6±3.4

It is known that measurement of vascular function is greatly influenced by external factors such as recent activity, diet, time of day, and so on. These influences greatly reduce the applicability of vascular function measurement to clinical practice since in the real world it is difficult or impossible to control these factors. For this reason, we did not ask participants to fast or refrain from physical activity or caffeine for either of their visits.



DOT and NIRS measurements were obtained using a dynamic near-infrared optical tomography (DYNOT) instrument (NIRx). The system illuminates the tissue with four laser diodes at wavelengths λ=725, 760, 810, and 830 nm.

The diodes are modulated at distinct frequencies and then coupled with 30, 1 mm multimode optical fibers—optodes—acting both as sources and detectors. Synchronous detection allows parallel measurement at a sampling frequency of 1.8 Hz.26 The optodes were organized in a hexagonal pattern with an interoptode spacing of 8 mm as shown in Fig. 1(a) and placed in a solid plastic holder to avoid movement artifacts [Fig. 1(b)]. For simplicity, hemoglobin concentration was calculated using only the wavelengths at 760 and 830 nm, chosen based on their symmetry with respect to the isobestic point of the extinction coefficients of hemoglobin.

Fig. 1

Experimental setup for NIRS and DOT measurement of reactive hyperemia. (a) Array with 30 optodes placed 8 mm apart in a hexagonal pattern. Optodes in red were used for NIRS measurements. All optodes were used for DOT reconstructions. (b) Placement of optode array on volar forearm of subject. (c) Rectangular finite element mesh used to model forearm tissue with dimensions 92  mm×35  mm×24  mm (W×L×D). Blue circles at z=0  mm indicate the location of the 30 optodes.



Near-Infrared Spectroscopy

NIRS employs the modified Beer–Lambert law (MBLL)


OD=  logII0=ϵCLB+G
to convert from changes in absorption to changes of de/oxyhemoglobin.5,27 In Eq. (1) OD is the optical density and I0 and I are incident and detected light intensities, respectively. ϵ represents the extinction coefficient of the tissue, and C is the concentration of the chromophore. L denotes the mean path length of detected photons. B is the path length factor, which accounts for the compensation of the increase of path lengths at various wavelengths caused by the scattering phenomena. G is defined as a geometric factor used to compensate the objective with different geometrical shapes. Typically, L, B, and G are constants with monochromatic illumination in a turbid media with unchanging geometry.

A change in the concentration of the chromophore causes a change in the intensity measured. The parameters ϵ and L remain constant, and it is assumed that B and G remain constant. Under this assumption, Eq. (1) can be written as


where ΔOD is the change in optical density. ODFinal and ODInitial are the detected optical density and the optical density of incident light. If and Ii are the measured intensities before and after the change in concentration; ΔC is the change in concentration. In Eq. (2), B is often referred to as the differential pathlength factor28 (DPF). DPF was determined experimentally for a number of tissues including forearm, calf, adult, and infant head.29 Changes in detected light are dominated mainly by oxygenated (HbO2) and deoxygenated hemoglobin (HbR) such that


where ϵHbO2λ, ϵHbRλ, Bλ are the extinction coefficients for oxyhemoglobin and deoxyhemoglobin and DPF at a given wavelength. By measuring the change in intensity at two wavelengths, it is possible to determine the concentration changes in HbO2 and HbR


Δ[HbR]=ϵHbO2λ2ΔODλ1Bλ1ϵHbO2λ1ΔODλ2Bλ2(ϵHbRλ1ϵHbO2λ2ϵHbRλ2εHbO2λ1)L,Δ[HbO2]=  ϵHbRλ1ΔODλ2Bλ2ϵHbRλ2ΔODλ1Bλ1(ϵHbRλ1ϵHbO2λ2ϵHbRλ2ϵHbO2λ1)L.
Typical source-detector separations in NIRS studies in reflectance mode are in the range 20 to 50 mm. NIRS measurements taken with the DYNOT equipment for separations larger than 45 mm were noisy (CV>5%), therefore this study was limited to a source-detector separation of L=32  mm obtained by selecting the central fibers indicated with red circles in Fig. 1(a). Although DPF depends on inter-optode spacing and wavelength, experiments carried out in the forearm suggest that this parameter becomes constant for distances larger than 25 mm.2930.31 For this study, DPF=4.0 was used for the calculation of concentration changes.


Diffuse Optical Tomography

Photon transport in tissue was modeled using the diffusion approximation of the radiative transport equation.32,33 This is a more accurate description of photon transport than the MBLL because it takes into account the random scattering of light produced by tissue. Consider the medium ΩR3 with boundary Ω, the diffusion equation in the steady-state domain is


where ϕi(r) is the spatially varying photon fluence at r due to source qi, μa is the absorption coefficient, D=[3(μa+μs)]1 is the diffusion coefficient, and μs is the reduced scattering coefficient. The source term represents an isotropic point source qi(r)=δ(rri) located at a depth of one scattering length inside the medium (d=1/μs). The boundary condition is usually of Robin type34


where the term A accounts for the refractive index boundary mismatch at the interface. The quantity measured by a detector located at ξjΩ, given the point source qi(r), is the outward flux Γi(ξj), and it is calculated from Fick’s law


where n(ξj) denotes the direction of the normal vector to the boundary at the detector location ξj. Equations (5)–(7) constitute the forward problem in DOT, which can also be represented by a scalar operator mapping35 between the space of optical parameters of interest, μa in this case, and the space of measurements as


where yi,j is the output of the j’th detector given the source i.


Image reconstruction

3-D volumetric reconstruction of the absorption coefficient was carried out using a modified version of the normalized difference approach proposed by Pei et al.36 At each sampling time tk, the 3-D absorption coefficient map μa(r,tk)=μaref(r)+Δμa(r,tk), discretized over the 3-D mesh, is determined by minimizing, using a nonlinear conjugated-gradient algorithm,37 the following cost function38


The first summation in the right side of Eq. (9) is over the volume elements of the 3-D mesh. In Eq. (9), Ns is the total number of sources, Nd,i is the number of detectors for the i’th source, yi,j(tk) is the flux measured by the j’th optode given the source i at time tk, yi,j0 is a reference state defined as the mean of the baseline measurements, and P^i,j[μaref(r)] and P^i,j[μa(r,tk)] are the model predicted measurements for the reference and estimated absorption coefficient at time tk. Essentially, for each time sample, the initial guess in the optimization process is the reference state   μaref(r). The conjugate gradient descent algorithm is applied iteratively for n steps to compute the estimate μa(r,tk). Note that the optimized image from a previous time sample can be used also as the initial guess. However, in this paper, concentration calculations were carried out independently of previous or subsequent samples.


Calculation of hemoglobin concentration

The absorption coefficient is related to the extinction coefficient and concentration as μa=εC. Assuming that the primary source of absorption changes is a combination of hemoglobin chromophores


where ϵHbO2λ and ϵHbRλ are the extinction coefficients for oxyhemoglobin and deoxyhemoglobin at a given wavelength. Reconstruction of absorption changes at two different wavelengths provides two independent equations, which can be solved simultaneously to calculate deoxyhemoglobin and oxyhemoglobin concentration changes as




Finite element modeling

A simplified finite element model of the skin and muscle tissue was used to implement the DOT reconstruction algorithms. The tissue was modeled as a rectangular cuboid with two layers [Fig. 1(c)] with dimensions 92  mm×35  mm×24  mm (length×width×depth). A tetrahedral mesh with 6900 elements and 1848 nodes [Fig. 1(c)] was generated for this domain. Each layer was assigned realistic optical parameters within the range of clinically normal tissue: μa=0.01  mm1 and μs=1.00  mm1 for skin39 and μa=0.02  mm1 and μs=0.50  mm1 for muscle.40

A simulation study was carried out to calculate the photon measurement density functions for different source-detector configurations. PDMFs characterize the regions of the tissue that contribute to the measurement signal and can be used to determine a suitable combination of measurements to reconstruct optical properties of tissue within a region of interest (ROI).

It is well known that the sensitivity in the reflectance mode is higher at the surface and diminishes as function of the depth.41,42 The study confirmed that for the chosen optode arrangement, depths larger than 12  mm contributed less than 5% to the sensitivity function. As a result, the ROI used in the analysis was a rectangular cuboid with dimensions 80  mm×20  mm×12  mm (length×width×depth) as shown in Fig. 2.

Fig. 2

The ROI is indicated with the volume in green. The location of the optodes in relation to the ROI is indicated with the blue circles. The dimensions of the ROI are 80  mm×20  mm×12  mm (length×width×depth). The ROI is located directly under the array of 30 optodes.


Yu et al.43 showed through experiments and simulation studies that for source–detector separations larger than 1 cm, the hyperemic response is mainly influenced by the autoregulation of muscle tissue, while smaller separations are primarily affected by subcutaneous tissue layer. For this reason, first- and second-nearest neighbor source–detector pairs were not used in the reconstruction.

3-D maps of absorption changes were obtained by minimizing the cost function in Eq. (9). To minimize the occurrence of artifacts, a two-step sign constraint algorithm was used during the reconstruction process.44 Positivity/negativity constraints were imposed after each iteration in separate reconstructions, and then the final absorption value was calculated as the sum of the two partial solutions. Positivity/negativity constraints force the algorithm to seek only positive/ negative changes in order to minimize the cost function given in Eq. (9). In each case, convergence was usually achieved after 5 to 10 iterations. Pei et al.44 demonstrated theoretically and experimentally that by using sign constraints, the quality and specificity of the recovered images improved significantly.

Finally, Eq. (11) was used to compute HbO2, HbR, and HbT (HbO2+HbR). Near-infrared parameters of postocclusive reactive hyperemia (PORH) were calculated on a nodal basis and then averaged on the ROI.


Study Protocol

The weight, height, and age for each subject were recorded. The baseline characteristics of the subjects are summarized in Table 1. Microvascular function during forearm PORH was assessed on two different occasions 1 to 2 weeks apart. In order to assess “real-world” performance, subjects were not asked to fast or abstain from exercise, and the time of day of assessment was not standardized. Subjects were requested to sit comfortably on a chair and to locate the right arm on a platform such that the dorsal-ventral axis of the elbow was held parallel to the surface of the platform [Fig. 1(b)]. The hand was placed such that the palm was facing upward. A manual pneumatic blood pressure cuff was placed above the elbow, and the tomographic probe was placed on the ventral side of the forearm on the brachioradialis muscle and attached with tape [Fig. 1(a)].

The test consisted of 3 min baseline, followed by 5 min of arterial occlusion (induced by inflating the cuff to 180 mmHg). After 5 min, the cuff was rapidly deflated and the response measured for a further 5 min (13 min total measurement time). Light attenuation changes were recorded at the four wavelengths, but only the data at 760 and 830 nm were used for processing. Hemodynamics changes were calculated using NIRS and DOT techniques. The former uses the MBLL to convert absorption changes into de/oxyhemoglobin5,27 [Eq. (4)], while the latter is based in the diffusion process of light [Eqs. (5) and (6)].35

The parameters measured during and after the arterial occlusion test are (Fig. 3):

  • (i) mVO2 (mlO2/min/100  g), muscle oxygen consumption: calculated from the gradient of the HbO2 at the beginning of the arterial occlusion and converted into ml O2/min/100  g.10

  • (ii) 1/2 THbO2 (s), 1/2 time of recovery of HbO2: time needed for half recovery of HbO2 from maximum deoxygenation at the end of the occlusion period to maximum reoxygenation during reactive hyperemia.

  • (iii) THbO2 (s), time to maximal HbO2: time interval after cuff release until maximum HbO2 value is reached.

  • (iv) MaxHbO2 (μM), maximal amplitude of HbO2: maximal hyperemic response reached by HbO2.

  • (v) MaxHbT (μM), maximal amplitude of HbT: maximal hyperemic response of HbT.

  • (vi) THbT (s), time to maximal HbT: time to peak value of HbT after cuff release.

  • (vii) Increase rate to maximal HbO2 (μM/s), reoxygenation rate of HbO2: calculated by dividing Max HbO2 by the time to maximal HbO2.

  • (viii) Increase rate to maximal HbT (μM/s), reoxygenation rate of HbT: calculated by dividing Max HbT by the time and maximal HbT.

  • (ix) AUC HbO2 (a.u.), area under the curve of HbO2 after cuff release.

  • (x) AUC HbR (a.u.), postocclusion area under the curve of HbT.

  • (xi) AUC HbT (a.u.), area under the curve of HbT after cuff release.

Fig. 3

Typical hemodynamic response before, during, and after 5 min brachial arterial occlusion measured by NIRS. Some characteristic parameters are indicated: oxygen consumption (mVO2), maximal amplitude of HbO2 (Max. HbO2), time to maximal of HbO2 (THbO2), and area under the curve of HbO2 (AUC).


The reproducibility was assessed using the intraclass correlation coefficient (ICC) and the standard error of measurement (SEM)45 based on NIRS parameters measured on two separate occasions. ICC provides the degree to which individuals maintain their rank order across repeated tests, and SEM is a measure of absolute repeatability. To interpret repeatability based on ICC values, the following guidelines were considered:46 ICC0.8 indicates “excellent,” 0.6ICC<0.8 is defined as “good,” 0.4ICC<0.6 indicates “moderate,” 0.2ICC<0.4 indicates “fair,” and ICC<0.2 is defined as “poor.” On the other hand, SEM is expressed as a percentage and the ideal value is zero. In general, the smaller the SEM, the smaller the variability between individual repeated measurements. ICC was calculated for near-infrared parameters derived using DOT and NIRS, and this is denoted by ICCDOT and ICCNIRS, respectively (similarly for SEM: SEMDOT and SEMNIRS).




Microvascular Function Reproducibility

Table 2 lists the mean and standard deviation of near-infrared parameters, calculated using DOT, during and after the arterial occlusion test, together with ICC and SEM repeatability measures. Similarly, Table 3 displays mean and standard deviation of near-infrared parameters measured using NIRS. Mean values in both tables are very similar; however, the standard deviation for near-infrared parameters derived using DOT is smaller in most cases. More importantly, the ICCDOT is higher in all but one case, and the repeatability is extremely high (1) for mVO2, THbO2, Max. HbO2, and THbT. In comparison, only one parameter (mVO2) measured with NIRS reached the “excellent” ICC category. The least repeatable parameter is AUC HbR, and this is consistent in the two measuring approaches.

Table 2

Intrasubject reproducibility of PORH parameters obtained using DOT.

ParameterTest 1Test 2ICCSEM
Inc. rate to max HbO21.53±0.511.46±0.550.7019.60
Inc. rate to max HbT0.72±0.390.62±0.470.5145.16
AUC HbO21461.79±572.511625.32±900.880.5432.22
AUC HbR708.47±586.48611.91±551.930.6352.11
AUC HbT1313.68±963.061741.28±775.310.2245.11

Table 3

Intrasubject reproducibility of PORH parameters obtained using NIRS spectroscopy.

ParameterTest 1Test 2ICCSEM
Inc. rate to max HbO21.06±0.520.85±0.380.2141.94
Inc. rate to max HbT0.73±0.400.60±0.450.4050.17
AUC HbO21053.21±740.531131.97±838.850.7833.89
AUC HbR810.67±657.291016.41±840.480.4561.01
AUC HbT1385.63±929.661897.82±836.570.1948.47

SEM calculated for tomographic imaging showed consistently smaller values across all parameters. This suggests that region-wise averaged measurements provide more robust and repeatable measurements than point-wise measurements, as typically measured in NIRS.



Our results are in agreement with similar studies using NIRS;20,47 however, the main finding is that tomographic imaging provides more robust and repeatable results. Figure 4 shows the mean and SEM of de/oxyhemoglobin and total hemoglobin across all subjects computed using NIRS and DOT. NIRS calculations are noisier than DOT reconstructions, despite using the same filtering techniques in both methods. One reason might be that single-channel measurements are more susceptible to systematic errors, such as position-dependent differences in muscle oxygenation or probe localization and interobserver reproducibility. On the other hand, DOT samples a larger tissue volume. Averaging over a large number of voxels also smooths the response. Furthermore, the image reconstruction method for DOT implements regularization that further improves the SNR.8

Fig. 4

Average (solid lines) and standard deviation (vertical bars) of HbO2, HbR, and HbT responses across all subjects in the first experiment obtained using (a) NIRS spectroscopy and (b) DOT.


Figure 5(a) shows the volume averaged HbO2, HbR, and HbT for a representative subject while Fig. 5(b) shows a tomographic reconstruction for HbO2 at its maximum value (t=516  s). The heterogeneity of tissue is evident, and this has been reported in similar studies involving two-dimensional images of venous occlusion tests.48 Furthermore, spatial dependence of the hemodynamic responses in the forearm has also been reported in nerve stimulation studies49 and in exercising hand in the reflectance50 and transmittance51 modes. This variability that occurs in all three spatial directions may contribute to diminished repeatability if the probe location is not consistent across examinations.

Fig. 5

(a) Typical hemodynamic responses of oxygenated (HbO2), de-oxygenated (HbR) and total hemoglobin (HbT) and (b) HbO2 (x,y,z,t) at t=516  s, i.e., the time average HbO2 reaches maximum.


NIRS and DOT calculations have clear methodological differences, both prone to systematic errors. For NIRS, the key parameter is the DPF, which corrects for the effect of scattering in tissue and therefore is wavelength and tissue dependent. This variable has been the subject of several investigations, and there is an agreement on its main characteristics at population levels.2930.31,52 However, in most of the cases, DPF is not calculated as part of the arterial/venous occlusion protocol, and calculations are based on data and tables available in the literature.

The techniques used in DOT are more complex and also susceptible to systematic errors such as parameterization of optical parameters, initial guess, or reconstruction convergence criteria. In this study, a differential approach was followed, which has demonstrated to be insensitive to boundary effects and the medium’s initial guess.36 On the other hand, the major disadvantage is that it is not possible to determine the absolute distribution of optical parameters, but only the change of absorption or diffusion from a given baseline.

In general, our results show time-based parameters have excellent repeatability. These parameters have shown great promise in distinguishing between healthy volunteers and patients with PAD and diabetes.16,17 For these patient groups, the reactive hyperemic response showed delayed times of recovery. Parameters derived from HbO2 are also the most repeatable, and this has been observed by different groups.16,20,21


Regions of Hemodynamic Consistence

The availability of 3-D maps allows the selection of more specific areas of interest. Selection techniques, either manual or automatic, are available and have been shown to enhance the response by improving the signal-to-noise ratio. Wang53 used a binarized segmentation technique to select the regions with the dominant response, while in another study by Cuccia et al.,54 regions with large vessels were avoided to focus on regions containing only microvasculature. Recently, Khalil et al.55 employed correlation analysis to define regions of temporal consistency that correlate significantly with the averaged HbT. Their results showed some differences in the evolution of HbT between a healthy volunteer and a patient with PAD, which were not evident from the average signals alone.

The best selection strategy is likely to be disease-dependent, and for this reason, we did not attempt to develop one based only on healthy volunteer subjects. However, to illustrate the potential advantages of DOT compared with NIRS, we carried out the following simple analysis:

  • Step 1. Compute the average total hemoglobin time series over the ROI, denoted by HbTROI(t), and determine the time point satisfying

    This point is located at tmax=467  s [Fig. 6(a)].

  • Step 2. Determine the node (xM,yM,zM) satisfying

    for t=467  s [Fig. 6(b)].

  • Step 3. For each node (xi,yi,zi) within the ROI, compute the cross-correlation between HbT(xM,yM,zM,t)=HbTmax(t) and HbT(xi,yi,zi,t)=HbTi(t) using


    ρ(xi,yi,zi)=  k=1N[HbTi(k)HbT¯i][HbTmax(k)HbT¯max]k=1N[HbTi(k)HbT¯i]2t=0N1[HbTref(k)HbT¯max]2,
    where N is the total number of samples and the bar over each variable denotes the time average over the duration of the experiment. The analysis was done on the entire experiment to illustrate the different dynamic signatures that can be obtained when different reference signals are used. Note also that HbTmax(t) is not necessarily the time series with the highest HbT across the experiment.

The resulting 3-D correlation map ρ(x,y,z), shown in Fig. 6(c), indicates very strong correlation (ρi>0.9) with the reactive hyperemic response HbTmax(t) for nodes located at a depth >2  mm, which largely correspond to muscle tissue.

Fig. 6

(a) Total hemoglobin curve (HbT) for a healthy subject; the point at maximal HbT (t=467  s) is marked with a point. (b) 3-D map of HbT(x,y,z) at the time of maximal HbT. The node with maximum contribution to the HbT response, HbT(xM,yM,zM), is indicated by the arrow and its time response is denoted with HbTmax(t). (c) 3-D map of correlation coefficients ρ(x,y,z) computed between HbTi(t) at each node location (xi, yi, zi) and HbTmax(t). (d) 3-D map of correlation coefficients ρ(x,y,z) computed between HbTi(t) at each node location (xi, yi, zi) and HbTmean(t).


In contrast, the correlation coefficients for nodes located in the superficial layer (skin) are very low (ρi0.4). A different correlation map was computed using as a reference the averaged HbT signal over the ROI, denoted HbTROI(t).

The corresponding 3-D correlation map, shown in Fig. 6(d), also exhibits strong correlation for nodes (ρi>0.9) located bellow the skin layer. The correlation maps can be used to define tissue compartments with similar hemodynamic profiles. For example, regions of hemodynamic consistency (RHC) can be defined based on the correlation maps shown in Fig. 6(c) and 6(d) by selecting all nodes with correlation coefficient ρ0.9.

The resulting RHCs are shown in Fig. 7(a); the red and green volumes indicate the regions with high correlation (ρ0.9) with the HbTmean and HbTmax(t), respectively. The blue volume represents the intersection between the two RHC. The average HbT time series for each RHC are shown in Fig. 7(b) together with the averaged HbT over the entire ROI [HbTROI(t)].

Fig. 7

(a) Compartments defined by the two RHC: the A and B volumes indicate the regions with high correlation (ρ0.9) with the HbTmean(t) and HbTmax(t), respectively. The AB represent the intersection between the the A and B regions. The boundaries of the ROI are denoted with the black box. (b) Total hemoglobin averaged over the entire ROI [HbTmean(t)] is indicated with the solid line. The dashed and dotted lines denote the average HbT response of all the nodes with high correlation (ρ0.9) with HbTmean(t) and HbTmax(t), respectively.


Interestingly, while the three signals are different before the end of the occlusion period, after the cuff is released only the region correlated to HbTmax(t) is still distinct compared with the average signal over the entire ROI. This tissue “compartment” identified by our analysis [Fig. 7(a)] appears to exhibit extensive reactive hyperemia both in terms of significantly increased blood flow relative to the baseline as well as duration of hyperemia.

The DOT-estimated MaxHbT and MaxHbO2 also have good reproducibility (ICC=0.69 and ICC=0.99, respectively), suggesting that an RHC identified based on HbTmax(t) could be a good candidate to focus on in future studies that aim to distinguish between healthy volunteers and patients with vascular disorder.

In contrast, the relatively poor reproducibility of MaxHbT measured using NIRS (ICC=0.37), which clearly limits its diagnostic potential, highlights the challenges posed by the spatial and temporal heterogeneity of tissue oxygenation and hemodynamics when assessing endothelial function.

To further illustrate the advantage of using DOT to resolve the spatiotemporal properties of the hemodynamic response, we created a short movie (Fig. 8) showing the reconstructed HbT changes in different slice planes along the x axis.

Fig. 8

(a) Dynamical changes of HbO2 (dotted line), HbR (dashed line), and HbT (solid line) for a healthy volunteer. (Video 1, top panel) The vertical line indicates the specific time point where a tomographic map is displayed in the lower panel. (b) Visualization of volumetric changes of HbT (video 1, bottom panel) using y-z slice planes located at x=30, 15, 0, 15, and 30 mm. (Video 1, MPEG-4, 941 KB) [URL: http://dx.doi.org/10.1117/1.JBO.21.6.066012.1].




This study shows that DOT achieves excellent reproducibility of key PORH parameters. These parameters include muscle oxygen consumption (mVO2), time to maximal HbO2 (THbO2), maximal HbO2 (MaxHbO2), and time to maximal HbT (THbT). Although a direct comparison of these parameters may be enough to distinguish between healthy volunteers and patients with vascular disease, our analysis suggests that obtaining reliable signatures of vascular disease may require first the identification of an appropriate RHC. This is only possible if a full volumetric reconstruction of hemodynamics, as that provided by DOT, is available for the particular ROI. The choice of an appropriate RHC is beyond the scope of this study, and it is the subject of future research.

The availability of volumetric hemodynamic parameters clearly offers more opportunities for the analysis and characterization of PORH and warrants further efforts to evaluate DOT’s potential to measure endothelial function in a clinical environment.

Despite the increased complexity of the instrumentation and reconstruction algorithms used to implement DOT, there are many examples of inexpensive portable and wearable DOT instruments.56 The availability of freely available software such as NIRFAST57 and TOAST58 also encourages the development of DOT-based technology. The use of reduced-order models14 to speed up the reconstruction process makes it now possible to perform analyses in real-time using mobile devices with modest computational resources.

Overall, the inherent advantages of DOT compared with other imaging modalities, combined with the availability of algorithms and portability of the instrumentation, makes DOT an ideal method for routinely and noninvasively assessing the cardiovascular function, inside and outside the hospital.


The authors gratefully acknowledge that this work was supported by the University of Sheffield as part of the “Engineering for Life” 2022 Futures Project and by EPSRC grants EP/H00453X/1 (S.A.B. and D.C.) and GR/T01006/01 (D.C.), by BBSRC grant BB/K010123/1 (S.A.B. and D.C.) and by an ERC Advanced Investigator Grant (S.A.B.). E.E.V.R. gratefully acknowledges the support from a grant from the Mexican National Research Council for Science and Technology (CONACYT).



G. A. Roth et al., “Global and regional patterns in cardiovascular mortality from 1990 to 2013,” Circulation 132(17), 1667–1678 (2015).JBOPFO1083-3668http://dx.doi.org/10.1161/CIRCULATIONAHA.114.008720CIRCAZ0009-7322Google Scholar


Cardiovascular-Disease-Team, Cardiovascular Disease Outcomes Strategy, Department of Health, London, p. 89 (2013).Google Scholar


D. A. Benaron et al., “Noninvasive methods for estimating in vivo oxygenation,” Clin. Pediatr. 31(5), 258–273 (1992).http://dx.doi.org/10.1177/000992289203100501Google Scholar


M. Ferrari, M. Muthalib and V. Quaresima, “The use of near-infrared spectroscopy in understanding skeletal muscle physiology: recent developments,” Philos. Trans. R Soc. A 369(1955), 4577–4590 (2011).http://dx.doi.org/10.1098/rsta.2011.0230Google Scholar


M. Cope and D. T. Delpy, “System for long-term measurement of cerebral blood and tissue oxygenation on newborn-infants by near-infrared trans-illumination,” Med. Biol. Eng. Comput. 26(3), 289–294 (1988).MBECDY0140-0118http://dx.doi.org/10.1007/BF02447083Google Scholar


F. F. Jobsis, “Noninvasive, infrared monitoring of cerebral and myocardial oxygen sufficiency and circulatory parameters,” Science 198(4323), 1264–1267 (1977).SCIEAS0036-8075http://dx.doi.org/10.1126/science.929199Google Scholar


D. A. Boas, A. M. Dale and M. A. Franceschini, “Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy,” Neuroimage 23, S275–S288 (2004).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2004.07.011Google Scholar


A. P. Gibson, J. C. Hebden and S. R. Arridge, “Recent advances in diffuse optical imaging,” Phys. Med. Biol. 50(4), R1–R43 (2005).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/50/4/R01Google Scholar


A. Villringer and B. Chance, “Non-invasive optical spectroscopy and imaging of human brain function,” Trends Neurosci. 20(10), 435–442 (1997).TNSCDR0166-2236http://dx.doi.org/10.1016/S0166-2236(97)01132-6Google Scholar


M. C. P. Van Beekvelt et al., “Performance of near-infrared spectroscopy in measuring local O 2 consumption and blood flow in skeletal muscle,” J. Appl. Physiol. 90(2), 511–519 (2001).Google Scholar


D. A. Boas et al., “The accuracy of near infrared spectroscopy and imaging during focal changes in cerebral hemodynamics,” Neuroimage 13(1), 76–90 (2001).NEIMEF1053-8119http://dx.doi.org/10.1006/nimg.2000.0674Google Scholar


C. Habermehl et al., “Somatosensory activation of two fingers can be discriminated with ultrahigh-density diffuse optical tomography,” NeuroImage 59(4), 3201–3211 (2012).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2011.11.062Google Scholar


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 18(4), 865–879 (2003).NEIMEF1053-8119http://dx.doi.org/10.1016/S1053-8119(03)00021-1Google Scholar


E. E. Vidal-Rosas et al., “Reduced-order modeling of light transport in tissue for real-time monitoring of brain hemodynamics using diffuse optical tomography,” J. Biomed. Opt. 19(2), 026008 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.2.026008Google Scholar


M. Vardi and A. Nini, “Near-infrared spectroscopy for evaluation of peripheral vascular disease. a systematic review of literature,” Eur. J. Vasc. Endovascular Surg. 35(1), 68–74 (2008).http://dx.doi.org/10.1016/j.ejvs.2007.07.015Google Scholar


T. Jarm et al., “Postocclusive reactive hyperemia in healthy volunteers and patients with peripheral vascular disease measured by three noninvasive methods,” in Advances in Experimental Medicine and Biology, , J. F. Dunn and H. M. Swartz, Eds., pp. 661–669, Springer, Boston, Massachusetts (2003).Google Scholar


R. Kragelj et al., “Parameters of postocclusive reactive hyperemia measured by near infrared spectroscopy in patients with peripheral vascular disease and in healthy volunteers,” Ann. Biomed. Eng. 29(4), 311–320 (2001).ABMECF0090-6964http://dx.doi.org/10.1114/1.1359451Google Scholar


M. Gayda et al., “Muscle VO2 and forearm blood flow repeatability during venous and arterial occlusions in healthy and coronary heart disease subjects,” Clin. Hemorheol. Microcirculation 59(2), 177–183 (2015).http://dx.doi.org/10.3233/CH-141836Google Scholar


R. C. Mesquita et al., “Diffuse optical characterization of an exercising patient group with peripheral artery disease,” J. Biomed. Opt. 18(5), 057007 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.5.057007Google Scholar


S. Lacroix et al., “Reproducibility of near-infrared spectroscopy parameters measured during brachial artery occlusion and reactive hyperemia in healthy men,” J. Biomed. Opt. 17(7), 077010 (2012).http://dx.doi.org/10.1117/1.JBO.17.7.077010Google Scholar


R. Kragelj, T. Jarm and D. Miklavcic, “Reproducibility of parameters of postocclusive reactive hyperemia measured by near infrared spectroscopy and transcutaneous oximetry,” Ann. Biomed. Eng. 28(2), 168–173 (2000).INPEEY0266-5611http://dx.doi.org/10.1114/1.241ABMECF0090-6964Google Scholar


M. A. Khalil et al., “Detection of peripheral arterial disease within the foot using vascular optical tomographic imaging: a clinical pilot study,” Eur. J. Vasc. Endovascular Surg. 49(1), 83–89 (2015).http://dx.doi.org/10.1016/j.ejvs.2014.10.010Google Scholar


A. Torricelli et al., “Functional brain imaging by multi-wavelength time-resolved near infrared spectroscopy,” Opto-electron. Rev. 16(2), 131–135 (2008).OELREM1230-3402http://dx.doi.org/10.2478/s11772-008-0011-6Google Scholar


T. Muehlemann, D. Haensse and M. Wolf, “Wireless miniaturized in-vivo near infrared imaging,” Opt. Express 16(14), 10323–10330 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.010323Google Scholar


K. J. Kek et al., “Optical imaging instrument for muscle oxygenation based on spatially resolved spectroscopy,” Opt. Express 16(22), 18173–18187 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.018173Google Scholar


C. H. Schmitz et al., “Dynamic studies of small animals with a four-color diffuse optical tomography imager,” Rev. Sci. Instrum. 76(9), 094302 (2005).RSINAK0034-6748http://dx.doi.org/10.1063/1.2038467Google Scholar


M. Cope, The Development of a Near-Infrared Spectroscopy System and Its Application for Noninvasive Monitoring of Cerebral Blood and Tissue Oxygenation, University College London, London (1991).Google Scholar


D. T. Delpy et al., “Estimation of optical pathlength through tissue from direct time of flight measurement,” Phys. Med. Biol. 33(12), 1433–1442 (1988).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/33/12/008Google Scholar


P. Van der Zee et al., “Experimentally measured optical pathlengths for the adult head, calf and forearm and the head of the newborn infant as a function of inter optode spacing,” Adv. Exp. Med. Biol. 316, 143–153 (1992).AEMBAP0065-2598http://dx.doi.org/10.1007/978-1-4615-3404-4Google Scholar


A. Duncan et al., “Optical pathlength measurements on adult head, calf and forearm and the head of the newborn-infant using phase-resolved optical spectroscopy,” Phys. Med. Biol. 40(2), 295–304 (1995).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/40/2/007Google Scholar


M. Ferrari et al., “Time-resolved spectroscopy of the human forearm,” J. Photochem. Photobiol. B 16(2), 141–153 (1992).JPPBEG1011-1344http://dx.doi.org/10.1016/1011-1344(92)80005-GGoogle Scholar


S. R. Arridge, “Methods in diffuse optical imaging,” Philos. Trans. R. Soc. A 369(1955), 4558–4576 (2011).PTRMAD1364-503Xhttp://dx.doi.org/10.1098/rsta.2011.0311Google Scholar


M. S. Patterson, B. Chance and B. C. Wilson, “Time resolved reflectance and transmittance for the noninvasive measurement of tissue optical-properties,” Appl. Opt. 28(12), 2331–2336 (1989).APOPAI0003-6935http://dx.doi.org/10.1364/AO.28.002331Google Scholar


M. Schweiger et al., “The finite element method for the propagation of light in scattering media: boundary and source conditions,” Med. Phys. 22(11), 1779–1792 (1995).MPHYA60094-2405http://dx.doi.org/10.1118/1.597634Google Scholar


S. R. Arridge, “Optical tomography in medical imaging,” Inverse Prob. 15(2) R41–R93 (1999).http://dx.doi.org/10.1088/0266-5611/15/2/022Google Scholar


Y. L. Pei, H. L. Graber and R. L. Barbour, “Influence of systematic errors in reference states on image quality and on stability of derived information for dc optical imaging,” Appl. Opt. 40(31), 5755–5769 (2001).JBOPFO1083-3668http://dx.doi.org/10.1364/AO.40.005755APOPAI0003-6935Google Scholar


J. Nocedal and S. Wright, Numerical Optimization, 2nd Ed., Springer, New York (2006).Google Scholar


A. Y. Bluestone et al., “Three-dimensional optical tomographic brain imaging in small animals, part 1: hypercapnia,” J. Biomed. Opt. 9(5), 1046–1062 (2004).http://dx.doi.org/10.1117/1.1784471Google Scholar


T. Lister, P. A. Wright and P. H. Chappell, “Optical properties of human skin,” J. Biomed. Opt. 17(9), 090901 (2012).http://dx.doi.org/10.1117/1.JBO.17.9.090901Google Scholar


S. L. Jacques, “Optical properties of biological tissues: a review,” Phys. Med. Biol. 58(11), R37–R61 (2013).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/58/11/R37Google Scholar


S. R. Arridge and M. Schweiger, “Photon-measurement density-functions. 2. Finite-element-method calculations,” Appl. Opt. 34(34), 8026–8037 (1995).APOPAI0003-6935http://dx.doi.org/10.1364/AO.34.008026Google Scholar


S. R. Arridge, “Photon-measurement density functions. Part I: analytical forms,” Appl. Opt. 34(31), 7395–7409 (1995).APOPAI0003-6935http://dx.doi.org/10.1364/AO.34.007395Google Scholar


G. Q. Yu et al., “Time-dependent blood flow and oxygenation in human skeletal muscles measured with noninvasive near-infrared diffuse optical spectroscopies,” J. Biomed. Opt. 10(2), 024027 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.1884603Google Scholar


Y. L. Pei, H. L. Graber and R. L. Barbour, “Normalized-constraint algorithm for minimizing inter-parameter crosstalk in DC optical tomography,” Opt. Express 9(2), 97–109 (2001).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.9.000097Google Scholar


T. A. Baumgartner, “Statistical tools in evaluation,” Chapter 2 in Statistical Tools in Evaluation, 3rd Ed., and A. S. Jackson, Eds., Houghton Mifflin, London (1987).Google Scholar


D. G. Altman, Practical Statistics for Medical Research, p. 623, Chapman & Hall, London (1999).Google Scholar


M. C. P. van Beekvelt et al., “In vivo quantitative near-infrared spectroscopy in skeletal muscle during incremental isometric handgrip exercise,” Clin. Physiol. Funct. Imaging 22(3), 210–217 (2002).http://dx.doi.org/10.1046/j.1475-097X.2002.00420.xGoogle Scholar


C. Casavola et al., “Blood flow and oxygen consumption with near-infrared spectroscopy and venous occlusion: spatial maps and the effect of time and pressure of inflation,” J. Biomed. Opt. 5(3), 269–276 (2000).JBOPFO1083-3668http://dx.doi.org/10.1117/1.429995Google Scholar


D. K. Chen et al., “Spectral and spatial dependence of diffuse optical signals in response to peripheral nerve stimulation,” Biomed. Opt. Express 1(3), 923–942 (2010).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.1.000923Google Scholar


M. Maris et al., “Functional near-infrared imaging of deoxygenated hemoglobin during exercise of the finger extensor muscles using the frequency-domain technique,” Bioimaging 2(4), 174–183 (1994).BOIMEL0966-9051http://dx.doi.org/10.1002/1361-6374(199412)2:4<174::AID-BIO2>3.3.CO;2-HGoogle Scholar


E. M. C. Hillman et al., “Time resolved optical tomography of the human forearm,” Phys. Med. Biol. 46(4), 1117–1130 (2001).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/46/4/315Google Scholar


S. J. Matcher, M. Cope and D.T. Delpy, “In vivo measurements of the wavelength dependence of tissue-scattering coefficients between 760 and 900 nm measured with time-resolved spectroscopy,” Appl. Opt. 36(1), 386–396 (1997).APOPAI0003-6935http://dx.doi.org/10.1364/AO.36.000386Google Scholar


C. Y. Wang et al., “Diffuse optical multipatch technique for tissue oxygenation monitoring: clinical study in intensive care unit,” IEEE Trans. Biomed. Eng. 59(1), 87–94 (2012).IEBEAX0018-9294http://dx.doi.org/10.1109/TBME.2011.2147315Google Scholar


D. J. Cuccia et al., “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt. 14(2), 024012 (2009).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3088140Google Scholar


M. A. Khalil et al., “Dynamic diffuse optical tomography imaging of peripheral arterial disease,” Biomed. Opt. Express 3(9), 2288–2298 (2012).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.002288Google Scholar


F. Scholkmann et al., “A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology,” NeuroImage 85, 6–27 (2014).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2013.05.004Google Scholar


H. Dehghani et al., “Near infrared optical tomography using NIRFAST: algorithm for numerical model and image reconstruction,” Commun. Numer. Methods Eng. 25(6), 711–732 (2009).CANMER0748-8025http://dx.doi.org/10.1002/cnm.v25:6Google Scholar


M. Schweiger and S. Arridge, “The Toast++ software suite for forward and inverse modeling in optical tomography,” J. Biomed. Opt. 19(4), 040801 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.4.040801Google Scholar


Ernesto E. Vidal-Rosas is a research associate at the University of Sheffield. He received his BEng and MSc degrees in electronic engineering from the Orizaba Institute of Technology and the National Centre for Research and Technology Development Cuernavaca in 2003 and 2006, respectively, and his PhD degree from the University of Sheffield in 2011. His main research interests include the development and application of novel methods for image reconstruction and analysis for diffuse optical tomography.

Stephen A. Billings is a professor in the Department of Automatic Control and Systems Engineering at the University of Sheffield and leads the Signal Processing and Complex Systems research group. His research interests include system identification and information processing for nonlinear systems, narmax methods, model validation, prediction, spectral analysis, adaptive systems, nonlinear systems analysis, wavelets, machine vision, cellular automata, spatiotemporal systems, EEG, fMRI and optical imagery of the brain, synthetic biology, crystal growth, and related fields.

Timothy Chico is reader in cardiovascular medicine at the University of Sheffield. He qualified in medicine from Leeds Medical School before undertaking an MD at the University of Sheffield. He previously held a GSK clinician scientist fellowship and works in the Bateson Centre lot the University of Sheffield. His research interests are in the effect of blood flow on vascular remodelling.

Daniel Coca is a professor at the University of Sheffield. He received his MEng degree in electrical engineering from Transilvania University of Brasov in 1993 and his PhD degree in automatic control and systems engineering from the University of Sheffield in 1997. His research interests include modeling, analysis, and control of complex systems, inverse problems, diffuse optical tomography and biological data analysis using reconfigurable computers.

Ernesto E. Vidal-Rosas, Stephen A. Billings, Timothy J. Chico, Daniel Coca, "Reproducibility of parameters of postocclusive reactive hyperemia measured by diffuse optical tomography," Journal of Biomedical Optics 21(6), 066012 (15 June 2016). http://dx.doi.org/10.1117/1.JBO.21.6.066012
Submission: Received ; Accepted

Near infrared spectroscopy



Tissue optics


Scanning electron microscopy


Back to Top