Research Papers: General

Visualization of light propagation in visible Chinese human head for functional near-infrared spectroscopy

[+] Author Affiliations
Ting Li, Hui Gong, Qingming Luo

Huazhong University of Science and Technology, Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics, 1037 Luoyu Road, Wuhan, Hubei 430074 China

J. Biomed. Opt. 16(4), 045001 (April 04, 2011). doi:10.1117/1.3567085
History: Received March 16, 2010; Revised January 17, 2011; Accepted February 24, 2011; Published April 04, 2011; Online April 04, 2011
Text Size: A A A

Open Access Open Access

* Address all correspondence to: Qingming Luo, Huazhong University of Science and Technology, Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics, 1037 Luoyu Road, Wuhan, Hubei 430074 China. Tel: 86 27 8779 2033; Fax: 86 27 8779 2034; E-mail: qluo@mail.hust.edu.cn

Using the visible Chinese human data set, which faithfully represents human anatomy, we visualize the light propagation in the head in detail based on Monte Carlo simulation. The simulation is verified to agree with published experimental results in terms of a differential path-length factor. The spatial sensitivity profile turns out to seem like a fat tropical fish with strong distortion along the folding cerebral surface. The sensitive brain region covers the gray matter and extends to the superficial white matter, leading to a large penetration depth (>3 cm). Finally, the optimal source-detector separation is suggested to be narrowed down to 3–3.5 cm, while the sensitivity of the detected signal to brain activation reaches the peak of 8%. These results indicate that the cerebral cortex folding geometry actually has substantial effects on light propagation, which should be necessarily considered for applications of functional near-infrared spectroscopy.

Figures in this Article

Functional near-infrared spectroscopy (fNIRS) and imaging have been used to measure functional brain activities.15 The sensitivity of the measured signal to brain activation and the brain volume sampled are crucial for clinical application; however, these factors are not feasible to be measured by experiment.67 Comparatively, simulations can give us reliable information of light propagation in the head67 by which we can obtain quantitative information of those factors in terms of a partial path-length factor (PPF) in brain tissue and the spatial sensitivity distribution of the fNIRS signal.69 Recently, simulation studies based on head models containing a cerebrospinal fluid (CSF) layer have demonstrated that light propagation is highly affected by the presence of CSF.7,1011 Because the cerebral surface is folded as gyri and sulci filling with CSF, it is probable that the cerebral cortex folding geometry is also important for light propagation in the head. Possibly, the gyri and sulci might change the shape of spatial sensitivity profiles.

Two investigations have been concerned with the effect of cerebral cortex folding geometry on light propagation. Okada et al. used a layered slab, with slots to imitate the sulci, and found little effect of the sulci on the spatial sensitivity profile for 3- or 4-cm source-detector separation (LSD).10 In the other study, using a two-dimensional (2-D) head model based on a single magnetic resonance imaging (MRI) scan, Fukui et al. showed that the spatial sensitivity profile formed like a banana, but distorted around the structure of cerebral surface, with LSD of 2–5 cm.6 Accordingly, the influence of cerebral cortex folding geometry on light propagation in a real human head is still unknown. A 3-D head model, most realistically representing the reliable and detailed anatomic structure of human head, is essential to examining or exploring the effect of cerebral cortex folding geometry on light propagation.

The visible Chinese human (VCH) data set was newly developed according to high-resolution cryosectional color photographs of a reference adult male.1214 Following a fine segmentation supervised by anatomists, a voxelized model (voxel size: 0.01×0.01 ×0.02 cm3) that most faithfully represented the human anatomical structures was established.1215 The VCH head is thus regarded as a most realistic head model containing the precise cerebral cortex folding geometry. To simulate light propagation in the VCH head, the simulation method applicable for 3-D voxelized media is needed, because the VCH head is a voxelized model. Actually, a software program targeted for Monte Carlo simulation of light propagation in voxelized media (MCVM), which is applicable for the VCH, has been developed.16 By using MCVM, our interest here is to evaluate the light propagation in the VCH with concern on the effect of cerebral cortex folding geometry.

In this study, near-infrared light propagation in the forehead of VCH was simulated. The values of the differential optical path-length factor (DPF) in cases of different LSD agreed well with published experimental data on real adult males. The spatial sensitivity profile was shown to have a fat tropical fish shape, with a larger sampled brain region covering the surface of the white matter, showing strong effect of cerebral cortex folding geometry. Lastly, the spatial sensitivity profile and PPF were presented under different LSD. The optimal LSD for fNIRS was suggested to be limited to 3–3.5 cm, while the sensitivity of the detected signal to brain activation reached the peak of 8%.

Anatomical Head Model of Visible Chinese Human

Details on developing the VCH phantom, including the image acquisition and segmentation, have been described in Refs. 1214. Briefly, the VCH specimen, a male adult, was frozen in the standing posture and horizontally sectioned at 0.02-cm intervals. Digital color photographs of the top surface after each sectioning were obtained with the resolution of 0.01 cm/pixel, which is higher than computed tomography (CT) or magnatic resonance imaging (MRI) data (such as Zubal model1718). The segmentation of VCH, unlike CT or MRI data, was performed based on these color photographs, with little dependence on pixel gray values. Under the direction of experienced anatomists, distinct tissues in each slice were visibly distinguished, manually segmented, and labeled with a specified identification number. The advantages, including high resolution and precise segmentation, make the VCH most realistically representative human anatomical structures. Particularly, the upstanding posture of VCH, instead of lying down in MRI or CT, was able to present a more realistic geometry of brain surface for light propagation prediction.

The color photographs and the segmented data sets of the VCH forehead were extracted for use. The head model has been segmented into eight types of tissue (see Fig. 1 for one slice), including scalp, skull, CSF, muscle, visible blood vessels, and gray and white matter. The optical properties of these tissues for near-infrared light at 800 nm are listed in Table 1.1921 In particular, these optical properties were equal to or close to the related previous studies, such as Ref. 6. For the segmented data set, every four consecutive pixels along both anterior-posterior and left-right directions in every two consecutive slices were chosen to form voxels each of 0.04×0.04×0.04 cm3. This was just the voxel size in our simulation on a common personal computer to balance the computation time and precision. We finally obtained a 3-D matrix of 450×400× 150 voxels containing the indices for different head tissue types. The indices were used to assign tissue-specific optical properties (Table 1) for each voxel.

Graphic Jump LocationF1 :

Visible Chinese human (VCH) head model: (a) Digital color photograph of one head slice in the VCH data set and (b) segmentation of the slice in (a).

Table Grahic Jump Location
Optical properties of head tissues for 800-nm light.
Monte Carlo Simulation

Monte Carlo simulation was used to predict light propagation in the VCH forehead. The Monte Carlo software used here was targeted for 3-D voxelized media (MCVM), the algorithm of which has been described in Ref. 16. Photon propagation in the head was determined by scattering coefficient, tiny absorption coefficient, phase function, etc. Notably, the scattering length was segmented into a series of substeps, each of which should be smaller than the voxel side length. Once the photon traveled a substep, it was checked whether the photon had crossed a boundary between different tissues to track each reaction rigorously and precisely. Particularly, when a photon crossed a refractive index mismatch boundary, reflection and transmission would be considered, including the surface of the head and interfaces between different tissue types. Additionally, the scattering length would be corrected when a photon crossed the boundary of voxels with different scattering coefficients.

The light sources and detectors were placed on the frontal head with LSD ranging from 1.5 to 5.1 cm. They were located in the slice (Fig. 1), ∼1 cm above the eyebrows. Each simulation was performed with 2×107 photons. The path lengths of detected photons in the head were weighted by the detected photon intensity and accumulated to calculate the DPF in the head and PPF only in brain tissues (the gray matter and white matter). Accordingly, the trajectories of the detected photons weighted were accumulated to obtain the spatial sensitivity distribution.10 Among these variables, DPF is the only one feasible to be measured6,10 and, thus, the predicted DPFs would be compared to the published experimental data on real adult males for validation of the simulation. Considering that all the calculations of these variables are based on photon trajectories-statistics, we believe that the accuracy of DPF would support the validation of the other variables predicted. The volume sampled by fNIRS was calculated as the spatial sensitivity distribution.68 The brain region sampled by fNIRS at varied LSD was analyzed in terms of spatial sensitivity profile only in brain tissues.

Differential Optical Path Length Factor

As mentioned above, calculated DPFs by Monte Carlo simulation with the VCH head model were compared to the experimental data measured on a real adult male forehead, for validation of the simulation. In the simulation, the light source was located on the middle forehead in the slice shown in Fig. 1, while the detectors were placed on the right forehead in the same slice with varied LSD. The experimental data were collected from results reported by Van der Zee et al.,22 and other groups.9,2324 In addition, the DPFs on the layered head model in Refs. 11 and 21 were calculated to analyze the effect of cerebral cortex folding geometry in human head. As shown in Fig. 2, in cases of different LSD, the DPFs calculated from the VCH head model (open circles with a solid line) agreed well with the mean DPFs measured by experiment from real male heads (filled circles). By contrast, the DPFs calculated from the layered head model (dotted line) were higher than all measured DPFs when LSD > 1.5 cm. Accordingly, the cerebral cortex folding geometry was important in improving simulation precision of DPFs.

Graphic Jump LocationF2 :

Comparison of calculated DPFs with experimentally measured values.

Spatial Sensitivity Profile

Video 1 was presented to show a full view of 3-D spatial sensitivity distribution in the VCH head slices from top to bottom. Particularly, the spatial sensitivity profile in the slice where the source and detector located was shown in Figs. 3(a1) and 3(a2). Here, the source (S) and detector (D) was placed on the left and right frontal head, respectively, with LSD of 3.9 cm. For comparison, simplified head models employed by Okada10 were also tested by our Monte Carlo simulation [Figs. 3(b1) and 3(b2)]. The gray lines in Fig. 3(b1) denoted the interface between different tissues in the layered head model including the outer layer, CSF, white matter, and gray matter. The models related to Figs. 3(b2) and 3(b3) added slots to the above layered model, the white region represented the scalp and skull, and the light gray represented CSF, while the dark gray and black represented the gray matter and white matter, respectively. The contours in all figures in this paper were drawn for 10%, 1% (semi-intense), and 0.1% of the maximum sensitivity.

Graphic Jump LocationF3 :

Spatial sensitivity profiles in the (a1, a2) VCH head model and (b1, b2, b3) simplified head models. The simplified head models in (b1, b2) were derived from Ref. 10. The width of slots in (b3) is triple of that in (b2).

The results showed that the spatial sensitivity profile in the VCH head [Figs. 3(a1) and 3(a2)] formed unlike the well-known “banana” shape [Fig. 3(b1)], which covered uniformly between the source and the detector. Instead, it looked like a fat tropical fish, and covered the gray matter and even the surface of the white matter. Significant distortions were observed around the cerebral cortex folding [Figs. 3(a1) and 3(a2)]. The spatial sensitivity profile penetrated deeper to the brain in the case of thicker CSF, no matter the light source and detector located in the opposite hemispheres [Fig. 3(a2)] or not [Fig. 3(a1)]. Particularly, it extended deeper by 30% under the central sulcus than regions nearby in Fig. 3(a2). The penetration depth and extra penetration depth under the sulcus over the nearby region were a little lower in Fig. 3(a1) than those in Fig. 3(a2). Accordingly, the cerebral cortex folding geometry was suggested to significantly affect the spatial sensitivity profile in human head. This was inconsistent with the previous finding based on the simplified head model [Fig. 3(b2)].10

In order to further examine the effect of cerebral cortex folding geometry, we modified the simplified head model in Fig. 3(b2). Note that the size of sulcus in the VCH head model was quite larger than 1 mm, which was used in the simplified model.10 Then we increased this size to 3 mm, and performed the simulation. As shown in Fig. 3(b3), the semi-intense (1%) sensitivity profile remarkably extended into the sulcus, whereas in the original model [Fig. 3(b2)], it did not. In addition, we simulated on Zubal MRI head model, of which voxel size equaled 1.1×1.1×1.4 mm, and observed the distortion of spatial sensitivity profile by cerebral cortex folding geometry (Fig. 4). However, the distortion on the VCH model was significantly stronger than that on the Zubal model. Therefore, the underestimation of sulcus size might be a crucial reason for the unobserved effect of cerebral cortex folding geometry in the previous study.10

Graphic Jump LocationF4 :

Spatial sensitivity profile of near-infrared light in Zubal head model with LSD sized of 3.5 cm.

Effect of Source-Detector Separation

Because LSD is a widely concerned factor, which affects both spatial sensitivity profile and PPF, we next explored the light propagation in the VCH frontal head in cases of different LSD (Figs. 56). The light source (S) was placed on the middle forehead, while the detector (D) was on the right forehead with varied LSD (Fig. 5).

Graphic Jump LocationF5 :

Spatial sensitivity profiles varied with LSD.

Graphic Jump LocationF6 :

PPF as a function of LSD. (a) PPFs for gray matter and white matter. (b) Percentage of PPFs for brain tissue in DPFs.

The spatial sensitivity profile can penetrate to the surface white matter, although the semisensitive brain region is distributed along the cerebral surface and is still confined to the gray matter. When LSD increased from 1.5 to 5.1 cm, the brain volume interrogated keeps nearly constant; however, the brain region with semi-intense sensitivity gets broader and the penetration depth of the semi-intense spatial sensitivity profile decreases. Given that the broadening of the sampling region would reduce the spatial resolution of fNIRS signal, smaller LSDs for adult head would be preferred.

Then, PPF for brain tissues (the gray matter and white matter, PPFbrain), indicating the sensitivity of detected light intensity signal to brain activation, and the ratio of PPFbrain to DPF, indicating the contribution of brain tissue to the detected signal were calculated. As shown in Fig. 6, photons that penetrated into the brain tissues could barely be detected in the case of LSD <1 cm. PPF for the gray matter increased sharply as LSD was >1 cm, kept nearly constant as long as LSD reached 2.5 cm, and decreased slightly as LSD was >3.5 cm [Fig. 6]. In contrast, the PPF for the white matter maintained small, even in the case of larger LSD. However, the positive PPF value indicated the existence of brain activation signals from the white matter in the detected signals. On the other hand, the PPFbrain/DPF ratio increases as LSD ranged 1–3 cm, and was sustained at 8% in the case of larger LSD. Hence, to obtain a larger PPFbrain component in the detected signal, LSD should be not <3 cm. As above, considering both the spatial sensitivity profile and PPFs in brain tissues, we suggested LSD be within 3–3.5 cm for adult heads.

On the basis of the VCH head model, the spatial sensitivity distribution and PPF in brain tissues for the adult male heads are predicted by Monte Carlo simulation. The coincidence of calculated DPFs of the VCH head with the experimentally measured DPFs of a real adult head proved our simulation, and our results showed that the cerebral cortex folding geometry indeed dramatically affected light propagation in the head: (i) The spatial sensitivity profile formed like a fat tropical fish. The sensitive brain region was distributed along the folding cerebral surface. (ii) The DPFs were closer to the mean values of experimentally measured DPFs and (iii) larger penetration depth. The sensitive brain region extended inward the cerebral cortex, covering the gray matter and even the surface of the white matter. (4) Larger PPFbrain proportion.

First, the previous study, showing little effect of cerebral cortex folding on light propagation, was performed on a simplified head model with slots mimicking the sulcus.10 Our simulation on that model also resulted in the common banana shape of the spatial sensitivity profile. However, increasing the sulcus width, the effect of the sulcus on the spatial sensitivity profile was obvious. The semi-intense (1%) sensitivity profile extended into the sulcus [Figs. 3(b2) and 3(b3)]. Therefore, in those studies using the 2-D or 3-D head models based on MRI slices of the adult head,6,2528 the sulcus width in the models was suggested to be not large enough to affect the spatial sensitivity profile (Fig. 4). And in this study based on the VCH head model, the sulcus width was apparently larger, leading to the transformation of the spatial sensitivity profile from the banana shape to the tropical fish shape. Take the central sulcus in the VCH forehead as an example, the width of which was particularly large. The spatial sensitivity profile in this region penetrated deeper than the nearby region by 30% in the case of LSD ≤ 4.2 cm. Notably, because of the high spatial resolution, tissue contrast, and segmentation precision of the VCH data set, the VCH head model is considered to most realistically represent the cerebral cortex folding geometry of the real adult head. Hence, the tropical fish shape of spatial sensitivity profile is believed to be more realistic.

Second, the DPFs in the VCH head model agreed well with the mean values of the experimentally measured DPFs on real human heads when LSD increase from 1 to 5.5 cm, whereas the DPFs in the layered head model11,21 almost exceeded the higher ends of experimentally measured DPFs as LSD ranged from 1.5 to 5.5 cm (Fig. 2). The published DPFs in a 2-D head model base on a single MRI slice6 were a little lower than those in the layered head model, but were still clearly higher than the mean values of the experimentally measured DPFs. Thus, the folding geometry of the cerebral cortex played an important role in improving the prediction precision of DPFs in adult heads.

Third, the penetration depth in the VCH head model exceeded previous expectation. The sampled brain volume has been commonly supposed to be confined to the shallow region of the gray matter, and the white matter is nearly impossible to be detected by the fNIRS signal.6 The penetration depth measured right beneath the source and detector in the case of commonly used LSD in fNIRS (2.5–4.5 cm) was previously predicted to be ∼2 cm.29 In contrast, this value was >3 cm in our VCH head model [Figs. 3(a1) and 3(a2)]. And, our results showed that the sensitive brain region covered the gray matter and penetrated into the superficial white matter (Figs. 3,5). The small but positive PPFs for the white matter also proved the existence of signal from white matter in the total detected signal [Fig. 6]. We consider this to be strong evidence to support that brain activation signals can be detected by fNIRS.

Lastly, while LSD > 3 cm, the PPFbrain proportion in the total DPF was demonstrated to be maximally 8%, which was higher than those (≤5%) in reports using simplified head models or MRI-based head models.6,10,25,29 Because this parameter denotes the contribution of brain tissue to the measurement of absolute absorption, the higher ratio value might be good news for detecting absolute concentration of oxy- and deoxy-hemoglobin, to interpret brain activities based on the neurovascular coupling hypothesis. Additionally, regarding the relative measurement in practical use, there would be nearly no changes in the signals from the skull and CSF (∼61% in our prediction), and a large portion (>15%) of blood in the scalp would be squeezed out by the probe pressure. Accordingly, <24% of the detected signals were involved in the relative measurement, and the proportion of that for brain activation was, therefore, >35%.

Additionally, there was another interesting result to support the effect of cerebral cortex folding geometry. We discovered that PPFs for gray matter were substantially deviated from its trend line as LSD ≤ 4.0 cm. If the cerebral surface were parallel to the skull, without cerebral cortex folding, then the PPFs for the gray matter should be right on the trend line, as in another published paper based on layered head models.11 However, in head models with cerebral cortex folding geometry, the geometry was heterogeneous along the light path, resulting in the deviation of PPFs from the trend line in cases of different LSD. Actually, a slight deviation was also observed in the report using the single-MRI slice–based 2-D head model.6

Furthermore, in this study, the optimization of LSD, which is an important concern for fNIRS detection, was explored. Previously, there was limited information on the optimal LSD in fNIRS measurement on adult heads.30 Commonly, the range of 2.5–4.5 cm was acceptable.6,28 One recent research suggested that 2.25–3.25 cm was the optimal range of LSD in detecting fast optical signal of brain activity.31 Comparatively, our results indicated a narrower optimal range of LSD. We discovered that with regard to the spatial sensitivity profile, larger LSD was actually not effective in increasing the sampled brain volume or penetration depth. Nevertheless, it produced broadening of the sampled brain region, which would reduce the spatial resolution of fNIRS detection. On the other hand, PPFs in brain tissues and its ratio to DPF both reached the peak at 3–3.5 cm of LSD (Fig. 6). Thus, the range of 3–3.5 cm was suggested to be optimal. Notably, this range actually agreed with the LSD (3 cm) in typical use, which was empirically determined in fNIRS multichannel detection of brain activation in adults.910,29

In summary, this study visualized the light propagation in the VCH head model based on Monte Carlo simulation. Cerebral cortex folding geometry was demonstrated to have strong effects on light propagation in the head. The spatial sensitivity profile formed like a fat tropical fish, showing strong effect of cerebral cortex folding geometry. The sensitive brain region penetrated into the surface layer of the white matter, leading to a much larger penetration depth compared to previous studies. The optimal LSD was proposed to be in the range of 3–3.5 cm, while the proportion of brain activation signal was maximally 8%.

This work was supported by the Program for Changjiang Scholars and Innovative Research Team in University, the 111 projects, and the Ph.D. Programs Foundation of Ministry of Education of China (Grant No. UNSPECIFIED 20070487058 ).

Izzetoglu  K., , Bunce  S., , Onaral  B., , Pourrezaei  K., , and Chance  B., “ Functional optical brain imaging using near-infrared during cognitive tasks. ,” Int. J. Hum.-Comput. Int.. 17, , 211–227  ((2004)).
Villringer  A., , and Chance  B., “ Non-invasive optical spectroscopy and imaging of human brain function. ,” Trends. Neurosci.. 20, , 435–442  ((1997)).
Li  T., , Luo  Q., , and Gong  H., “ Gender-specific hemodynamics in prefrontal cortex during a visual-verbal working memory task by near-infrared spectroscopy. ,” Behav. Brain. Res.. 209, , 148–153  ((2010)).
Li  T., , Li  L., , Luo  Q., , and Gong  H., “ Assessing working memory in real-life situations with functional near-infrared spectroscopy. ,” J. Innovat. Opt. Health. Sci.. 2, , 423–430  ((2009)).
Irani  F., , Platek  S. M., , Bunce  S., , Ruocco  A. C., , and Chute  D., “ Functional near infrared spectroscopy (fNIRS): An emerging neuroimaging technology with important applications for the study of brain disorders. ,” Clin. Neuropsychol.. 21, , 9–37  ((2007)).
Fukui  Y., , Ajichi  Y., , and Okada  E., “ Monte Carlo prediction of near-infrared light propagation in realistic adult and neonatal head models. ,” Appl. optics. 42, , 2881–2887  ((2003)).
Hoshi  Y., “ Functional near-infrared spectroscopy: current status and future prospects. ,” J. Biomed. Opt.. 12, , 062106  ((2007)).
Okada  E., , Firbank  M., , and Delpy  D. T., “ The effect of overlying tissue on the spatial sensitivity profile of near-infrared spectroscopy. ,” Phys. Med. Boil.. 40, , 2093–2108  ((1995)).
Hoshi  Y., , Shimada  M., , Sato  C., , and Iguchi  Y., “ Reevaluation of near-infrared light propagation in the adult human head: implications for functional near-infrared spectroscopy. ,” J. Biomed. Opt.. 10, , 064032  ((2005)).
Okada  E., , Firbank  M., , Schweiger  M., , Arridge  S. R., , Cope  M., , and Delpy  D. T., “ Theoretical and experimental investigation of near-infrared light propagation in a model of the adult head. ,” Appl. Opt.. 36, , 21–31  ((1997)).
Okada  E., , and Delpy  D. T., “ Near-infrared light propagation in an adult head model. I. modeling of low-level scattering in the cerebrospinal fluid layer. ,” Appl. Opt.. 42, , 2906–2914  ((2003)).
Zhang  S. X., , Heng  P. A., , Liu  Z. J., , Tan  L. W., , Qiu  M. G., , Li  Q. Y., , Liao  R. X., , Li  K., , Cui  G. Y., , Guo  Y. L., , Yang  X. P., , Liu  G. J., , Shan  J. L., , Liu  J. J., , Zhang  W. G., , Chen  X. H., , Chen  J. H., , Wang  J., , Chen  W., , Lu  M., , You  J., , Pang  X. L., , Xiao  H., , Xie  Y. M., , and Cheng  J. C., “ The Chinese visible human (CVH) datasets incorporate technical and imaging advances on earlier digital humans. ,” J. Anat.. 204, , 165–173  ((2004)).
Li  A., , Liu  Q., , Zeng  S., , Tang  L., , Zhong  S., , and Luo  Q., “ Construction and visualization of high-resolution three-dimensional anatomical structure datasets for Chinese digital human. ,” Chin. Sci. Bull.. 53, , 1848–1854  ((2008)).
Zhang  G., , Luo  Q., , Zeng  S., , and Liu  Q., “ The development and application of the visible Chinese human model for Monte Carlo dose calculations. ,” Health Phys.. 94, , 118–125  ((2008)).
Zhang  G. Z., , Liu  Q., , and Luo  Q. M., “ Monte Carlo simulations for external neutron dosimetry based on the visible Chinese human phantom. ,” Phys. Med. Boil.. 52, , 7367–7383  ((2007)).
Li  T., , Luo  Q., , and Gong  H., “ MCVM: Monte Carlo modeling of photon migration in voxelized media. ,” J. Innovat. Opt. Health Sci.. 3, , 91–102  ((2010)).
Zubal  I. G., , Harrell  C. R., , Smith  E. O., , Rattner  Z., , Gindi  G., , and Hoffer  P. B., “ Computerized three-dimensional segmented human anatomy. ,” Med. Phys.. 21, , 299–302  ((1994)).
Toronov  V., , D’Amico  E., , Hueber  D., , Gratton  E., , Barbieri  B., , and Webb  A., “ Optimization of the signal-to-noise ratio of frequency-domain instrumentation for near-infrared spectro-imaging of the human brain. ,” Opt. Express.. 11, , 2717–2729  ((2003))
Cheong  W. F., , Prahl  S. A., , and Welch  A. J., “ A review of the optical properties of biological tissues. ,” IEEE J. Quantum Electron.. 26, , 2166–2185  ((1990)).
Yaroslavsky  A. N., , Schulze  P. C., , and Yaroslavsky  I. V., “ Optical properties of selected native and coagulated human brain tissues in vitro in the visible and near infrared spectral range. ,” Phys. Med. Biol. , 47, , 2059–2073  ((2002)).
Okada  E., and Delpy  D. T., “ 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, , 2915–2922  ((2003)).
van der Zee  P., , Cope  M., , Arridge  S. R., , Essenpreis  M., , Potter  L. A., , Edwards  A. D., , Wyatt  J. S., , McCormick  D. C., , Roth  S. C., , Reynolds  E. O. R., , and Deply  D. T., “ 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)).
Duncan  A., , Meek  J. H., , Clemence  M., , Elwell  C. E., , Tyszczuk  L., , Cope  M., , and Delpy  D. T., “ 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, , 295–304  ((1995)).
Zhao  H. J., , Tanikawa  Y., , Gao  F., , Onodera  Y., , Sassaroli  A., , Tanaka  K., , and Yamada  Y., “ Maps of optical differential pathlength factor of human adult forehead, somatosensory motor and occipital regions at multi-wavelengths in NIR. ,” Phys. Med. Biol. , 47, , 2075–2093  ((2002)).
Strangman  G., , Franceschini  M. A., , and Boas  D. A., “ Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters. ,” NeuroImage. 18, , 865–879  ((2003)).
Franceschini  M. A., and Boas  D. A., “ Noninvasive measurement of neuronal activity with near-infrared optical imaging. ,” NeuroImage. 21, , 372–386  ((2004)).
Boas  D. A., and Dale  A. M., “ Simulation study of magnetic resonance imaging-guided cortically constrained diffuse optical tomography of human brain function. ,” Appl. Opt.. 44, , 1957–1968  ((2005)).
Mudra  R., , Nadler  A., , Keller  E., , and Niederer  P., “ Analysis of near-infrared spectroscopy and indocyanine green dye dilution with Monte Carlo simulation of light propagation in the adult brain. ,” J. Biomed. Opt.. 11, , 044009  ((2006)).
Firbank  M., , Okada  E., , and Delpy  D. T., “ A theoretical study of the signal contribution of regions of the adult head to near-infrared spectroscopy studies of visual evoked responses. ,” NeuroImage. 8, , 69–78  ((1998)).
Taga  G., , Homae  F., , and Watanabe  H., “ Effects of source-detector distance of near infrared spectroscopy on the measurement of the cortical hemodynamic response in infants. ,” NeuroImage. 38, , 452–460  ((2007)).
Gratton  G., , Brumback  C. R., , Gordon  B. A., , Pearson  M. A., , Low  K. A., , and Fabiani  M., “ Effects of measurement method, wavelength, and source-detector distance on the fast optical signal. ,” NeuroImage. 32, , 1576–1590  ((2006)).
© 2011 Society of Photo-Optical Instrumentation Engineers (SPIE)

Citation

Ting Li ; Hui Gong and Qingming Luo
"Visualization of light propagation in visible Chinese human head for functional near-infrared spectroscopy", J. Biomed. Opt. 16(4), 045001 (April 04, 2011). ; http://dx.doi.org/10.1117/1.3567085


Figures

Graphic Jump LocationF1 :

Visible Chinese human (VCH) head model: (a) Digital color photograph of one head slice in the VCH data set and (b) segmentation of the slice in (a).

Graphic Jump LocationF2 :

Comparison of calculated DPFs with experimentally measured values.

Graphic Jump LocationF3 :

Spatial sensitivity profiles in the (a1, a2) VCH head model and (b1, b2, b3) simplified head models. The simplified head models in (b1, b2) were derived from Ref. 10. The width of slots in (b3) is triple of that in (b2).

Graphic Jump LocationF4 :

Spatial sensitivity profile of near-infrared light in Zubal head model with LSD sized of 3.5 cm.

Graphic Jump LocationF5 :

Spatial sensitivity profiles varied with LSD.

Graphic Jump LocationF6 :

PPF as a function of LSD. (a) PPFs for gray matter and white matter. (b) Percentage of PPFs for brain tissue in DPFs.

Tables

Table Grahic Jump Location
Optical properties of head tissues for 800-nm light.

References

Izzetoglu  K., , Bunce  S., , Onaral  B., , Pourrezaei  K., , and Chance  B., “ Functional optical brain imaging using near-infrared during cognitive tasks. ,” Int. J. Hum.-Comput. Int.. 17, , 211–227  ((2004)).
Villringer  A., , and Chance  B., “ Non-invasive optical spectroscopy and imaging of human brain function. ,” Trends. Neurosci.. 20, , 435–442  ((1997)).
Li  T., , Luo  Q., , and Gong  H., “ Gender-specific hemodynamics in prefrontal cortex during a visual-verbal working memory task by near-infrared spectroscopy. ,” Behav. Brain. Res.. 209, , 148–153  ((2010)).
Li  T., , Li  L., , Luo  Q., , and Gong  H., “ Assessing working memory in real-life situations with functional near-infrared spectroscopy. ,” J. Innovat. Opt. Health. Sci.. 2, , 423–430  ((2009)).
Irani  F., , Platek  S. M., , Bunce  S., , Ruocco  A. C., , and Chute  D., “ Functional near infrared spectroscopy (fNIRS): An emerging neuroimaging technology with important applications for the study of brain disorders. ,” Clin. Neuropsychol.. 21, , 9–37  ((2007)).
Fukui  Y., , Ajichi  Y., , and Okada  E., “ Monte Carlo prediction of near-infrared light propagation in realistic adult and neonatal head models. ,” Appl. optics. 42, , 2881–2887  ((2003)).
Hoshi  Y., “ Functional near-infrared spectroscopy: current status and future prospects. ,” J. Biomed. Opt.. 12, , 062106  ((2007)).
Okada  E., , Firbank  M., , and Delpy  D. T., “ The effect of overlying tissue on the spatial sensitivity profile of near-infrared spectroscopy. ,” Phys. Med. Boil.. 40, , 2093–2108  ((1995)).
Hoshi  Y., , Shimada  M., , Sato  C., , and Iguchi  Y., “ Reevaluation of near-infrared light propagation in the adult human head: implications for functional near-infrared spectroscopy. ,” J. Biomed. Opt.. 10, , 064032  ((2005)).
Okada  E., , Firbank  M., , Schweiger  M., , Arridge  S. R., , Cope  M., , and Delpy  D. T., “ Theoretical and experimental investigation of near-infrared light propagation in a model of the adult head. ,” Appl. Opt.. 36, , 21–31  ((1997)).
Okada  E., , and Delpy  D. T., “ Near-infrared light propagation in an adult head model. I. modeling of low-level scattering in the cerebrospinal fluid layer. ,” Appl. Opt.. 42, , 2906–2914  ((2003)).
Zhang  S. X., , Heng  P. A., , Liu  Z. J., , Tan  L. W., , Qiu  M. G., , Li  Q. Y., , Liao  R. X., , Li  K., , Cui  G. Y., , Guo  Y. L., , Yang  X. P., , Liu  G. J., , Shan  J. L., , Liu  J. J., , Zhang  W. G., , Chen  X. H., , Chen  J. H., , Wang  J., , Chen  W., , Lu  M., , You  J., , Pang  X. L., , Xiao  H., , Xie  Y. M., , and Cheng  J. C., “ The Chinese visible human (CVH) datasets incorporate technical and imaging advances on earlier digital humans. ,” J. Anat.. 204, , 165–173  ((2004)).
Li  A., , Liu  Q., , Zeng  S., , Tang  L., , Zhong  S., , and Luo  Q., “ Construction and visualization of high-resolution three-dimensional anatomical structure datasets for Chinese digital human. ,” Chin. Sci. Bull.. 53, , 1848–1854  ((2008)).
Zhang  G., , Luo  Q., , Zeng  S., , and Liu  Q., “ The development and application of the visible Chinese human model for Monte Carlo dose calculations. ,” Health Phys.. 94, , 118–125  ((2008)).
Zhang  G. Z., , Liu  Q., , and Luo  Q. M., “ Monte Carlo simulations for external neutron dosimetry based on the visible Chinese human phantom. ,” Phys. Med. Boil.. 52, , 7367–7383  ((2007)).
Li  T., , Luo  Q., , and Gong  H., “ MCVM: Monte Carlo modeling of photon migration in voxelized media. ,” J. Innovat. Opt. Health Sci.. 3, , 91–102  ((2010)).
Zubal  I. G., , Harrell  C. R., , Smith  E. O., , Rattner  Z., , Gindi  G., , and Hoffer  P. B., “ Computerized three-dimensional segmented human anatomy. ,” Med. Phys.. 21, , 299–302  ((1994)).
Toronov  V., , D’Amico  E., , Hueber  D., , Gratton  E., , Barbieri  B., , and Webb  A., “ Optimization of the signal-to-noise ratio of frequency-domain instrumentation for near-infrared spectro-imaging of the human brain. ,” Opt. Express.. 11, , 2717–2729  ((2003))
Cheong  W. F., , Prahl  S. A., , and Welch  A. J., “ A review of the optical properties of biological tissues. ,” IEEE J. Quantum Electron.. 26, , 2166–2185  ((1990)).
Yaroslavsky  A. N., , Schulze  P. C., , and Yaroslavsky  I. V., “ Optical properties of selected native and coagulated human brain tissues in vitro in the visible and near infrared spectral range. ,” Phys. Med. Biol. , 47, , 2059–2073  ((2002)).
Okada  E., and Delpy  D. T., “ 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, , 2915–2922  ((2003)).
van der Zee  P., , Cope  M., , Arridge  S. R., , Essenpreis  M., , Potter  L. A., , Edwards  A. D., , Wyatt  J. S., , McCormick  D. C., , Roth  S. C., , Reynolds  E. O. R., , and Deply  D. T., “ 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)).
Duncan  A., , Meek  J. H., , Clemence  M., , Elwell  C. E., , Tyszczuk  L., , Cope  M., , and Delpy  D. T., “ 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, , 295–304  ((1995)).
Zhao  H. J., , Tanikawa  Y., , Gao  F., , Onodera  Y., , Sassaroli  A., , Tanaka  K., , and Yamada  Y., “ Maps of optical differential pathlength factor of human adult forehead, somatosensory motor and occipital regions at multi-wavelengths in NIR. ,” Phys. Med. Biol. , 47, , 2075–2093  ((2002)).
Strangman  G., , Franceschini  M. A., , and Boas  D. A., “ Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters. ,” NeuroImage. 18, , 865–879  ((2003)).
Franceschini  M. A., and Boas  D. A., “ Noninvasive measurement of neuronal activity with near-infrared optical imaging. ,” NeuroImage. 21, , 372–386  ((2004)).
Boas  D. A., and Dale  A. M., “ Simulation study of magnetic resonance imaging-guided cortically constrained diffuse optical tomography of human brain function. ,” Appl. Opt.. 44, , 1957–1968  ((2005)).
Mudra  R., , Nadler  A., , Keller  E., , and Niederer  P., “ Analysis of near-infrared spectroscopy and indocyanine green dye dilution with Monte Carlo simulation of light propagation in the adult brain. ,” J. Biomed. Opt.. 11, , 044009  ((2006)).
Firbank  M., , Okada  E., , and Delpy  D. T., “ A theoretical study of the signal contribution of regions of the adult head to near-infrared spectroscopy studies of visual evoked responses. ,” NeuroImage. 8, , 69–78  ((1998)).
Taga  G., , Homae  F., , and Watanabe  H., “ Effects of source-detector distance of near infrared spectroscopy on the measurement of the cortical hemodynamic response in infants. ,” NeuroImage. 38, , 452–460  ((2007)).
Gratton  G., , Brumback  C. R., , Gordon  B. A., , Pearson  M. A., , Low  K. A., , and Fabiani  M., “ Effects of measurement method, wavelength, and source-detector distance on the fast optical signal. ,” NeuroImage. 32, , 1576–1590  ((2006)).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.