Open Access
9 March 2018 Voxel-based measurement sensitivity of spatially resolved near-infrared spectroscopy in layered tissues
Author Affiliations +
Abstract
We quantitatively investigated the measurement sensitivity of spatially resolved spectroscopy (SRS) across six tissue models: cerebral tissue, a small animal brain, the forehead of a fetus, an adult brain, forearm muscle, and thigh muscle. The optical path length in the voxel of the model was analyzed using Monte Carlo simulations. It was found that the measurement sensitivity can be represented as the product of the change in the absorption coefficient and the difference in optical path length in two states with different source–detector distances. The results clarified the sensitivity ratio between the surface layer and the deep layer at each source–detector distance for each model and identified changes in the deep measurement area when one of the detectors was close to the light source. A comparison was made with the results from continuous-wave spectroscopy. The study also identified measurement challenges that arise when the surface layer is inhomogeneous. Findings on the measurement sensitivity of SRS at each voxel and in each layer can support the correct interpretation of measured values when near-infrared oximetry or functional near-infrared spectroscopy is used to investigate different tissue structures.

1.

Introduction

Tissue oximetry using near-infrared spectroscopy (NIRS) is used in a range of fields, including brain research, sports medicine, surgery, and obstetrics.1,2 Implantable devices and devices mounted on the finger of the investigator have been developed,3,4 and the range of tissue types that can be examined has widened. NIRS instruments apply four basic techniques: time-resolved spectroscopy (TRS), spatially resolved spectroscopy (SRS), frequency-domain spectroscopy (FDS), and continuous-wave spectroscopy (CWS). TRS, FDS, and CWS record changes in light intensity at a single point, whereas SRS is a multipoint approach. CWS, which uses an adaptation of the Beer–Lambert law, is the simplest to apply and the most widely used. Recent TRS5 using data processing methods and high-performance FDS6 have been developed to measure scattering and absorption coefficients, but these systems have not been widely introduced to clinical tests because it is a complex and expensive system. Both SRS and CWS are simple and low-cost approaches to light radiation and detection and benefit from a high signal-to-noise ratio. SRS is used to measure the absolute values of hemoglobin concentration, and its measurement sensitivity can be suppressed when examining superficial tissue.7 This is because analysis is based on the difference in intensity at two points. A number of previous studies have addressed the spatial sensitivity distribution of CWS811 and diffuse optical tomography (DOT).12,13 Dehghani et al.12 and Saikia and Kanhirodan13 analyzed the influence of source–detector separation and depth-related sensitivity for reconstruction of DOT. However, as its voxel-based sensitivity for SRS oximetry has not yet been quantified, there is a need for further investigation of the roles played by source–detector separation and by the various layered structures of the target tissue. In the current study, we examined the SRS sensitivity in the source–detector axis direction and depth direction and compared the difference in measurement sensitivity of SRS and CWS at each voxel and layer.

2.

Methods

On the basis of radiative transfer theory, we performed a Monte Carlo analysis of the six tissue types shown in Fig. 1. For the directly contacted cerebral tissue and the small animal brain, a source–detector separation of up to 20 mm was used. This was increased to 40 mm when tested on the fetus, adult head, and limb muscles, mimicking the separation used by commercially available instruments or reported by research studies.14 Figure 1 shows the six model tissues used in the Monte Carlo analysis: (1) cerebral tissue contacted directly by an optical probe, (2) a small animal brain accessed via the scalp, (3) the forehead of a fetus, (4) an adult brain, (5) forearm muscle, and (6) thigh muscle. The optical properties shown in Table 1 and used in the simulations were based on previous literature.1419 To allow the simulation results to be applied to NIRS data interpretation, the models were divided into a surface part and a deep part. Figure 1(a) shows a surface to deep combination of gray matter and white matter, Figs. 1(b)1(d) show the scalp plus gray matter and white matter, and Figs. 1(e) and 1(f) show skin and muscle. The thickness of the surface component was set to 1.5 mm to represent the scalp and skin thickness of 1 to 2 mm.20 The sensitivity of CWS is known to be determined by the mean optical path length L. The relationship between path length and SRS sensitivity is as follows. IA is the light intensity detected over a short distance ρA from the light source, and IB is the intensity detected over a longer distance ρB. As the spatial intensity slope S is the difference in light intensity per unit length, S is defined as ln(IA/IB). The change in S given a change in the absorption coefficient μa is then defined as the measurement sensitivity of SRS. The values of S are taken from optical measurements, and the unknown absorption coefficients are derived from the values of S and the theoretical curves.

Fig. 1

Simulation models for (a) cerebral tissue, (b) small animal brain, (c) fetal brain, (d) adult brain, (e) forearm muscle, and (f) thigh muscle.

JBO_23_3_030503_f001.png

Table 1

Optical properties for each tissue.

TissueReduced scattering coefficient μs′ (mm−1)Absorption coefficient μa (mm−1)
Scalp1.30.020
Skull2.00.010
Cerebrospinal fluid0.30.002
Gray matter1.60.035
White matter5.00.015
Skin1.30.020
Fat1.20.003
Muscle0.70.025

In Eq. (1), suffixes 1 and 2 denote values before and after the small change in absorption, respectively. The ΔS associated with a change in optical density ΔOD and the mean path lengths LA and LB is then derived as follows:

Eq. (1)

ΔS=ln(I2A/I2B)ln(I1A/I1B)=ln(I2A/I1A)ln(I2B/I1B)=ΔODA+ΔODB=Δμa(LBLA).

In CWS, sensitivity is usually based on the mean path length, whereas in SRS the sensitivity is also influenced by the difference between LA and LB. On the basis of this relationship, the sensitivity of SRS for a small voxel at point (x,y,z) is given by

Eq. (2)

ΔSx,y,z=Δμax,y,z(LBx,y,zLAx,y,z).
By preliminary analysis to investigate the slope change due to an increase in the whole μa or the voxel μa by 1%, it was verified that Eqs. (1) and (2) held when Δμa was small and the relationship between S and μa could be regarded as linear. The sensitivity for each layer is derived by integrating the voxel values in the same layer. The SRS sensitivity for the surface and deep layers is derived as follows:

Eq. (3)

ΔSsurface=Δμasurface(LBsurfaceLAsurface),

Eq. (4)

ΔSdeep=Δμadeep(LBdeepLAdeep).
For example, LAsurface denotes the intensity-weighted average path length through the surface layer that the light detected at a position A has traversed. To clarify the influence of the change of S on the oximetry, the absorption coefficient is calculated by the commonly used equation 21 as follows:

Eq. (5)

μa13μs(SρBρA2ρA+ρB2)2.

The concentrations of oxyhemoglobin [O2Hb], deoxyhemoglobin [HHb], and tissue oxygen saturation StO2 are obtained by

Eq. (6)

[O2Hb]=ϵHHbλ2μaλ1ϵHHbλ1μaλ2ϵO2Hbλ1ϵHHbλ2ϵO2Hbλ2ϵHHbλ1,

Eq. (7)

[HHb]=ϵO2Hbλ2μaλ1ϵO2Hbλ1μaλ2ϵO2Hbλ1ϵHHbλ2ϵO2Hbλ2ϵHHbλ1,

Eq. (8)

StO2=[O2Hb][O2Hb]+[HHb],
where ϵO2Hbλ1,λ2 and ϵHHbλ1,λ2 are the extinction coefficients of HHb and O2Hb, respectively, at the wavelengths λ1 and λ2.22 In this study, the influence when the position of the absorber changes in the direction of the source–detector axis was also analyzed. We assumed two wavelengths of 770 and 830 nm and placed a 10×10×3.5-mm hypoxic region (StO2: 30%, μa770=0.044  mm1, μa830=0.032  mm1) in the normal cortex (StO2: 63%, μa770=0.035  mm1, μa830=0.035  mm1).

3.

Results and Discussion

Figure 2 shows the voxel-based sensitivity ΔSx,y,z from two model detector combinations: ρA=7  mm, ρB=35  mm, and ρA=30  mm, ρB=35  mm. The voxel was 0.5 mm on each side. The values of 10 voxels near the center (2.5<y<2.5) were summed to clarify the trend. The red/yellow/green plots represent positive changes in spatial slope for a 10% increase in absorption, and the blue/cyan plots represent negative changes. In both Figs. 2(a) and 2(b), the sensitivity of the surface layer was negative close to detector A and positive close to detector B. When blood is in the negative sensitivity area, an incorrect reverse tendency will be mixed into the calculation results of the oxygen saturation and blood volume. At a ρA of 30 mm, a strong positive sensitivity was noted within a narrow range in the deep layer, but a small negative sensitivity appeared at the side near the light source. Figure 3 shows the oxygen saturation due to the change in the relative position between the optical probe and the 30%-StO2 region. When light was detected in a pair of 7 to 35 mm, StO2 decreased by 4.5% with the approach of the low oxygen tissue. When a pair of 30 to 35 mm was used, it decreased by 8%, but StO2 at x=5  mm increased by 1% despite low oxygenation. These results suggest the following: (1) when the surface layer is heterogeneous, as is the case for veins, wounds, or inflammation, the effect on the measured value is large and can be positive or negative. (2) When both detectors A and B are distant from the light source, many sites have negative sensitivity, but information can be obtained from a narrow region of the deep layer. (3) When the gap between ρA and ρB is sufficiently long, average information can be obtained across a wide area of the deep tissue.

Fig. 2

SRS sensitivity distribution from an adult head model: (a) ρA=7  mm, ρB=35  mm and (b) ρA=30  mm, ρB=35  mm.

JBO_23_3_030503_f002.png

Fig. 3

Changes in tissue oxygen saturation StO2 due to position of low oxygen tissue in the source–detector axis direction.

JBO_23_3_030503_f003.png

Figure 4 shows the measurement sensitivity of each layer from every model for both CWS and SRS. Because SRS makes use of two distances, the SRS results are plotted on the horizontal axis using ρB. The color coding corresponds to the color given to each layer in Fig. 1. In all the models, the CWS sensitivity for both layers and the SRS sensitivity for the deep layer increased as the source–detector distance ρ increased, whereas the SRS sensitivity in the surface layer decreased. The cancelation produced by the positive and negative sensitivities reduced the sensitivity in the surface layer. In direct contact cerebral measurement [Fig. 4(a)], the region of interest is often the surface layer. When ρ was <15  mm, the sensitivity of both SRS and CWS in the surface layer was greater. As shown in Figs. 4(b)4(f), hemodynamics in the surface layer may affect the measured values. In the case of the small animal model [Fig. 4(b)], deep sensitivity was dominant, especially when using SRS. This suggests that (i) cerebral tissue can be sufficiently measured only after shaving and (ii) craniotomy is not necessary. In the fetus or neonate model brain [Fig. 4(c)], SRS measurement with the detector pair at 4 to 6 mm showed almost the same sensitivity in the surface and deep layers. Even when a probe with a ρ of <10  mm was used, sensitivity in the deep layer remained relatively strong. In the model of the adult brain [Fig. 4(d)], it was necessary to set ρB to 28 mm or more. The muscle measurement model used both thin [Fig. 4(e)] and thick [Fig. 4(f)] fat layers. The influence of skin was very slight when ρB>20  mm. It was demonstrated that, across a wide range of conditions, the influence of the surface layer could be canceled when SRS measurement was used as long as the layer was homogeneous. To cancel the changes in the surface tissue, it is necessary to avoid applying force only near the detectors and to avoid contact with heterogeneous parts such as the subcutaneous veins. By comparing the SRS sensitivity reported in this study with the CWS sensitivity reported in earlier studies,3,4,711 ways of suppressing or extracting the influence of the surface layer may be suggested.

Fig. 4

Measurement sensitivity of SRS and CWS for each model, each layer, and each detector pair. Symbols (a)–(f) correspond to Fig. 1. Detector pair of (b) is the same as (a), and pairs of (d)–(f) are shown in (c).

JBO_23_3_030503_f004.png

4.

Conclusions

The analysis results of the sensitivity of the voxel sensitivity and the influence of the low- StO2 tissue position quantitatively showed one of the factors of false negative and false positive signals on NIRS that was discussed in the previous study.23,24 By modeling actual measurements, we also compared the difference in sensitivity of each voxel and in each layer of SRS and CWS quantitatively. Measurement points and the ratio between the deep layer and the surface layer in measurement were clarified. Our findings on the measurement sensitivity of SRS at each voxel or layer can be used for hemodynamic interpretation of the measured values from different tissues.

Disclosures

Authors have no financial interests associated with the results of this study.

Acknowledgments

This work was supported in part by Japan Society for the Promotion of Science (JSPS) under Grants-in-Aid for Scientific Research (25350525).

References

1. 

A. Villringer et al., “Near infrared spectroscopy (NIRS): a new tool to study hemodynamic changes during activation of brain function in human adults,” Neurosci. Lett., 154 (1–2), 101 –104 (1993). https://doi.org/10.1016/0304-3940(93)90181-J NELED5 0304-3940 Google Scholar

2. 

T. Hamaoka et al., “The use of muscle near-infrared spectroscopy in sport, health and medical sciences: recent developments,” Philos. Trans. R. Soc. A: Math., Phys. Eng. Sci., 369 (1955), 4591 –4604 (2011). https://doi.org/10.1098/rsta.2011.0298 Google Scholar

3. 

M. Niwayama and T. Yamakawa, “Implantable thin NIRS probe design and sensitivity distribution analysis,” Electron. Lett., 50 (5), 346 –348 (2014). https://doi.org/10.1049/el.2013.3921 ELLEAK 0013-5194 Google Scholar

4. 

N. Kanayama and M. Niwayama, “Examiner’s finger-mounted fetal tissue oximetry,” J. Biomed. Opt., 19 (6), 067008 (2014). https://doi.org/10.1117/1.JBO.19.6.067008 JBOPFO 1083-3668 Google Scholar

5. 

K. Mehta et al., “Spread spectrum time-resolved diffuse optical measurement system for enhanced sensitivity in detecting human brain activity,” J. Biomed. Opt., 22 (4), 045005 (2017). https://doi.org/10.1117/1.JBO.22.4.045005 JBOPFO 1083-3668 Google Scholar

6. 

B. B. Zimmermann et al., “Frequency domain near-infrared multiwavelength imager design using high-speed, direct analog-to-digital conversion,” J. Biomed. Opt., 21 (1), 016010 (2016). https://doi.org/10.1117/1.JBO.21.1.016010 JBOPFO 1083-3668 Google Scholar

7. 

L. M. L. Dix et al., “Comparing near-infrared spectroscopy devices and their sensors for monitoring regional cerebral oxygen saturation in the neonate,” Pediatr. Res., 74 (5), 557 –563 (2013). https://doi.org/10.1038/pr.2013.133 PEREBL 0031-3998 Google Scholar

8. 

E. Okada and D. T. Delpy, “Near-infrared light propagation in an adult head model. II. Effect of superficial tissue thickness on the sensitivity of the near-infrared spectroscopy signal,” Appl. Opt., 42 (16), 2915 –2922 (2003). https://doi.org/10.1364/AO.42.002915 APOPAI 0003-6935 Google Scholar

9. 

A. V. Patil et al., “Experimental investigation of NIRS spatial sensitivity,” Biomed. Opt. Express, 2 (6), 1478 –1493 (2011). https://doi.org/10.1364/BOE.2.001478 BOEICL 2156-7085 Google Scholar

10. 

G. E. Strangman, Z. Li and Q. Zhang, “Depth sensitivity and source-detector separations for near infrared spectroscopy based on the Colin27 brain template,” PLoS One, 8 (8), e66319 (2013). https://doi.org/10.1371/journal.pone.0066319 POLNCL 1932-6203 Google Scholar

11. 

A. Messere and S. Roatta, “Influence of cutaneous and muscular circulation on spatially resolved versus standard Beer–Lambert near-infrared spectroscopy,” Physiol. Rep., 1 e00179 (2013). https://doi.org/10.1002/phy2.179 Google Scholar

12. 

H. Dehghani et al., “Depth sensitivity and image reconstruction analysis of dense imaging arrays for mapping brain function with diffuse optical tomography,” Appl. Opt., 48 (10), D137 –D143 (2009). https://doi.org/10.1364/AO.48.00D137 APOPAI 0003-6935 Google Scholar

13. 

M. J. Saikia and R. Kanhirodan, “Region-of-interest diffuse optical tomography system,” Rev. Sci. Instrum., 87 (1), 013701 (2016). https://doi.org/10.1063/1.4939054 RSINAK 0034-6748 Google Scholar

14. 

S. Wan, R. R. Anderson and J. A. Parrish, “Analytical modeling for the optical properties of the skin with in vitro and in vivo applications,” Photochem. Photobiol., 34 (4), 493 –499 (1981). https://doi.org/10.1111/php.1981.34.issue-4 PHCBAP 0031-8655 Google Scholar

15. 

H. R. Eggert and V. Blazek, “Optical properties of human brain tissue, meninges, and brain tumors in the spectral range of 200 to 900 nm,” Neurosurgery, 21 (4), 459 –464 (1987). https://doi.org/10.1227/00006123-198710000-00003 NEQUEB Google Scholar

16. 

P. Vanderzee, M. Essenpreis and D. T. Delpy, “optical-properties of brain-tissue,” Proc. SPIE, 1888 454 –465 (1993). https://doi.org/10.1117/12.154665 PSISDG 0277-786X Google Scholar

17. 

G. Mitic et al., “Time-gated transillumination of biological tissues and tissuelike phantoms,” Appl. Opt., 33 (28), 6699 –6710 (1994). https://doi.org/10.1364/AO.33.006699 APOPAI 0003-6935 Google Scholar

18. 

J. L. Robertson et al., “Effect of blood in the cerebrospinal fluid on the accuracy of cerebral oxygenation measured by near infrared spectroscopy,” Adv. Exp. Med. Biol., 812 233 –240 (2014). https://doi.org/10.1007/978-1-4939-0620-8 AEMBAP 0065-2598 Google Scholar

19. 

M. Niwayama et al., “Quantitative measurement of muscle hemoglobin oxygenation using near-infrared spectroscopy with correction for the influence of a subcutaneous fat layer,” Rev. Sci. Instrum., 71 (12), 4571 –4575 (2000). https://doi.org/10.1063/1.1322578 RSINAK 0034-6748 Google Scholar

20. 

H. M. Hori, G. Rebora and A. Crovato, “The thickness of human scalp: normal and bald,” J. Invest. Dermatol., 58 (6), 396 –399 (1972). https://doi.org/10.1111/1523-1747.ep12540633 Google Scholar

21. 

S. J. Matcher et al., “Absolute quantification methods in tissue near infrared spectroscopy,” Proc. SPIE, 2389 486 –495 (1995). https://doi.org/10.1117/12.209997 PSISDG 0277-786X Google Scholar

22. 

S. J. Matcher et al., “Performance comparison of several published tissue near-infrared spectroscopy algorithms,” Anal. Biochem., 227 (1), 54 –68 (1995). https://doi.org/10.1006/abio.1995.1252 Google Scholar

23. 

I. Tachtsidis and F. Scholkmann, “False positives and false negatives in functional near-infrared spectroscopy: issues, challenges, and the way forward,” Neurophotonics, 3 (3), 031405 (2016). https://doi.org/10.1117/1.NPh.3.3.031405 Google Scholar

24. 

S. Greenberg, T. Shear and G. Murphy, “Extracranial contamination of near-infrared spectroscopy devices,” Anesth. Analg., 124 (1), 356 –358 (2017). https://doi.org/10.1213/ANE.0000000000001290 Google Scholar
CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Masatsugu Niwayama "Voxel-based measurement sensitivity of spatially resolved near-infrared spectroscopy in layered tissues," Journal of Biomedical Optics 23(3), 030503 (9 March 2018). https://doi.org/10.1117/1.JBO.23.3.030503
Received: 12 November 2017; Accepted: 21 February 2018; Published: 9 March 2018
Lens.org Logo
CITATIONS
Cited by 9 scholarly publications and 1 patent.
Advertisement
Advertisement
KEYWORDS
Tissues

Near infrared spectroscopy

Brain

Natural surfaces

Sensors

Absorption

Animal model studies

Back to Top