 Research
 Open Access
 Published:
Patientoriented simulation based on Monte Carlo algorithm by using MRI data
BioMedical Engineering OnLine volume 11, Article number: 21 (2012)
Abstract
Background
Although Monte Carlo simulations of light propagation in full segmented threedimensional MRI based anatomical models of the human head have been reported in many articles. To our knowledge, there is no patientoriented simulation for individualized calibration with NIRS measurement. Thus, we offer an approach for brain modeling based on image segmentation process with in vivo MRI T1 threedimensional image to investigate the individualized calibration for NIRS measurement with Monte Carlo simulation.
Methods
In this study, an individualized brain is modeled based on in vivo MRI 3D image as five layers structure. The behavior of photon migration was studied for this individualized brain detections based on threedimensional timeresolved Monte Carlo algorithm. During the Monte Carlo iteration, all photon paths were traced with various sourcedetector separations for characterization of brain structure to provide helpful information for individualized design of NIRS system.
Results
Our results indicate that the patientoriented simulation can provide significant characteristics on the optimal choice of sourcedetector separation within 3.3 cm of individualized design in this case. Significant distortions were observed around the cerebral cortex folding. The spatial sensitivity profile penetrated deeper to the brain in the case of expanded CSF. This finding suggests that the optical method may provide not only functional signal from brain activation but also structural information of brain atrophy with the expanded CSF layer. The proposed modeling method also provides multiwavelength for NIRS simulation to approach the practical NIRS measurement.
Conclusions
In this study, the threedimensional timeresolved brain modeling method approaches the realistic human brain that provides useful information for NIRS systematic design and calibration for individualized case with prior MRI data.
Background
Nearinfrared spectroscopy (NIRS) is a promising noninvasive brain imaging technique with a higher sampling rate than positron emission tomography (PET)/functional magnetic resonance imaging (fMRI) and a more precise and localized spatial resolution than Electroencephalography (EEG)/Magnetoencephalography (MEG). The NIRS technique provides information about the slow signal (i.e., hemoglobin response) and fast signal (i.e., neuronal activation) [1–5]. This optical method permitted several benefits as noninvasive, less expensive, nonionizing radiation imaging, realtime measurement, compact implementation, long time monitoring and easy operation with high time resolution and adequate spatial resolution for continuously recording oxy and deoxyhemoglobin changes of brain. Also, NIRS offers a more comprehensive measurement of brain activity than bloodoxygenationleveldependent (BOLD) fMRI.
Functional nearinfrared brain imaging is achieved with the backscattering light detection by using sourcedetector pairs on the surface of human head [5–7]. For NIRS implementation, there are several issues that including signaltonoise ratio evaluation, optimal choice of sourcedetector separation, the brain structural effects on light propagation and the brain volume sampled remain to be fully understood well. Therefore, the simulation approach is important for characterization of photon migration in human brain with various sourcedetector separations to provide helpful information for individualized design of NIRS system [8–10].
In the most previous studies, the simulation results were generally obtained by semiinfinity fivelayer structure [11–22] or twodimensional head model with a MRI slice [13, 18]. Naturally, the threedimensional brain structure modeling by utilizing in vivo MRI data provides a realistic phenomenon of photon migration dynamics. However, there is no detail description for efficient and systematic modeling method of Monte Carlo algorithm with threedimensional anatomical MRI data [14, 15, 17, 19]. Additionally, the threedimensional model which faithfully represents the realistic human head from MRI data depends on image processing.
Therefore, we offer a systematic approach for 3D brain modeling based on image segmentation process with in vivo MRI T1 threedimensional image. For investigation of individualized difference in brain structure with NIRS, an adult volunteer was modeling to implement Monte Carlo simulated with various sourcedetector separations. According to previously studies, the light guiding effect occurred in the CSF layer of human brain. The presence of a relatively clear layer such as CSF that has both low scattering and absorption coefficients has been shown especially to alter the light propagation in the head [12, 16, 19, 23, 24]. This phenomenon cannot be portrayed by diffusion approximation method because the CSF reveals low scattering property [25] but it can be observed in the Monte Carlo simulation. Accordingly, the result indicates the advantage of the Monte Carlo method for NIRS modeling. Besides, the NIRS system typically applies multiwavelength sources to detect the concentration changes of oxy and deoxyhemoglobin such as 690 nm, 780 nm and 830 nm. Therefore, this study offers a NIRS simulation method for understanding photon migration dynamics in human brain by using threedimensional MRI data with multiwavelength illumination.
Methods
Three dimensional brain MRI T1 data processing
Figure 1 shows an in vivo MRI T1 image of human brain with five layers that assigned as scalp, skull, cerebral spinal fluid (CSF), gray matter and white matter, respectively. The threedimensional brain image contains 256 × 256 × 92 voxels and each voxel size is 1 × 1 × 1 mm^{3}.
Segmentation methods are important technique used in image processing to identify the objects in the image. For segmentation of brain layers, the image process includes two steps: 1) to segment the scalp and skull by level set method and region growing approach and 2) to segment the CSF, gray matter and white matter by using probabilistic framework segmentation. Figure 2 demonstrates the contour segmentation of scalp and skull that was achieved with level set and region growing operator. The level set method is a numerical technique for tracking interfaces and shapes. The advantage of the level set method are that it is implicit, parameter free, provides a direct way to estimate the geometric properties of the evolving structure, can change the topology and is intrinsic [26–29].
In the Figure 2(a), as could be expected, there are holes (boundary gaps) in the edges of skull layer, so simply coloring and parametric form of deformable models between the edges would not work. In contrast, level sets are designed to handle topological changes naturally. However, unlike the parametric form, they are not robust to boundary gaps and suffer from several other deficiencies as well. Thus, the level set method is suitable to segment the contour of scalp and skull from MRI data. In two dimensions, the level set method amounts to representing a closed curve Γ using a level set function ø (t, x, y). The closed curve Γ is represented as the zero level set of ø (t, x, y) by [29]
If the curve Γ(t) moves in the normal direction with a speed function F, then the level set function ø satisfies the level set equation that can be written in the following general form:
To avoid the problem of further computation highly inaccurate with develop shocks, very sharp or flat shape during the evolution. The level set function ø have to reshape (also call reinitialize) to be a signed distance function periodically during the evolution. The variational level set formulation which keeps as an approximate signed distance function has been proposed that can be easily implemented by simple finite difference scheme, without the need of reinitialization as following the formula:
where μ > 0 is a parameter controlling the effect of penalizing the deviation of ø from a signed distance function; P(ø) is considered as the internal energy, defined within the curve, are designed to keep the model smooth during the deformation process. Consider a unit circle Ω⊂ℜ^{2} that can be written as:
while the ε_{m} (ø) is considered as the external energy which are computed from the underlying image data, are defined to move the model toward an object boundary or other desired features within the image. The formula of external energy can be written as:
where λ > 0 and ν are constants; L _{ g } (ø) is the length of the zero level curve of ø and A _{ g } (ø) is introduced to speed up curve evolution of level set function that are defined as:
where δ is the univariae Dirac function, and H is the Heaviside function; g is the edge indicator function defined by
where G _{ α }is the Gaussian kernel with standard deviation α and I is an image.
Figure 2(b) shows the results obtained using level sets to segment scalp and skull from the background. To obtain the outer boundary of the scalp, we started with an initial level set at the boundary of the image. To obtain the structures of skull on the inside of the brain, we started with an initial level set that was a closed curve around a point on the inside. This curve evolved to identify the boundaries of skull inside the brain. The contours generated by the level sets are closed contours.
After image segmentation by utilizing level set method, the contours of scalp and skull were segmented as shown in Figure 2(c) and 2(d). The region growing approach was then adopted to segment the two connected layers in binary image. The basic idea of region growing was starting with seeds. The grow regions from corresponding seeds revealed similar properties with their neighboring pixels [30, 31]. According to the result of region growing segmentation, the scalp and skull layers were distinguished and marked as type 1 and 2 in simulation (Figure 2(e) and 2(f)). Figure 2(g) shows the two layers modeling of the scalp and skull.
After scalp and skull labeling, the probabilistic framework was then applied to classify CSF, gray matter and white matter layers with unified segmentation, which was performed by fitting a mixture of Gaussians (MOG) model with prior information of deformable tissue probability maps [32]. The MOG model can be described by the probability density of intensity y _{i} and k _{th} Gaussian distribution with mean μ _{k} and variance σ _{k} ^{2} as
The Gaussian function indicated the probability distribution of brain tissues. The bias correction and image registration were included within the unified segmentation approach. Figure 3 shows the probability distributions of CSF, gray matter and white matter by unified segmentation from in vivo MRI data.
According to the mapping of probability distribution, each image pixel would be sorted to CSF, gray matter, or white matter by calculated maximum probability of tissue type. The layers of CSF, gray matter and white matter were assigned as type 3, 4 and 5 in simulation. Figure 4(a) shows the fivelayer brain structure after image process method and Figure 4(b) demonstrates the reconstructed threedimensional brain model by 92 slices that corresponded to original in vivo threedimensional MRI data (Figure 4(c)).
Monte Carlo algorithm
Photon migration in human tissue can be numerically simulated by the Monte Carlo algorithm [33–40]. The photon trajectory can be computed with the parameters for propagation governing: 1) the mean free path of a scattering or absorption event, 2) the boundary conditions  refraction and specular reflection, 3) scattering event  deflection and azimuth angles, 4) absorption event  energy loss, 5) detector location. The Monte Carlo model relies on the sampling of random variables from their probability density function. The probability density function of step size s _{1} is defined as:
According to the probability density function as uniform distribution, s _{1} can be obtained:
where ξ is uniformly distributed between [0,1]. Consider a multilayer structure that the photon may experience a free path over multiple layers before a scattering event, the counterpart of Eq. (11) becomes:
where μ _{ ti } is the extinction coefficient and s _{ i } is the path length of the i _{th} voxel. The Snell's law and Fresnel reflection formulas were applied at each boundary. The probability for the occurrence of a scattering process over the distance d _{s} was given by
The probability distribution of scattering angle θ was assumed by Henyey and Greenstein function with the anisotropy factor g as
The probability distribution of azimuth angle Ψ was assumed to be isotropic as
where Ψ = 2πξ with the uniform random number ξ ∈ [0, 1]. At each scattering event, an individual photon packet dropped part of its power and the energy loss can be represented by
where w is the weight of the photon packet before the scattering event. In addition, the formal solution, Mie theory, describes absorption and/or scattering event with a sphere that has been available in previously study [41]. The photonpassed voxels were all recorded with temporal evolution for photon footprint tracing. Therefore, the dynamic behavior of photon migration in human brain can be manifested. Additionally, we recorded all the paths of the received photons in the simulations, the visited layers of each photon were marked. Accordingly, spatial sensitivity profiles (SSP) of adult head models were calculated from the accumulated trajectories of photons. The spatial sensitivity has been described theoretically by photon measurement density functions or sensitivity maps [39, 42, 43].
The source and multidetectors arrangement on the surface with transverse view and sagittal view were applied to investigate the light propagation in human brain (shown in Figure 5). The sourcedetector separations in transverse view and sagittal view are 110 cm with 1 cm step. First, all cases were simulated at typical 800 nm wavelength illumination with the fivelayer scattering/absorption coefficients 1.9/0.018 (scalp), 1.6/0.016 (skull), 0.24/0.004 (CSF), 2.2/0.036 (gray matter) and 9.1/0.014 (white matter) mm^{1} [16, 21, 23, 36]. For NIRS modeling, multiwavelength sources (690, 780, and 830 nm) were applied for illumination of the adult brain. The reduced scattering coefficient μ _{s}', absorption coefficient μ _{a}, scatters' radius, refractive indices of background and scatters of brain tissues are described in Table 1 [35].
Results
Figure 6 shows the tomograms of the adult brain structures of in vivo MRI data and processed optical models for Monte Carlo simulation. The 92 twodimensional slices were used from head top to down and then the threedimensional images were reconstructed of both structures. The depth of head model was 9.2 cm.
The light source was located on the frontal surface at the 6 cm from the head top in transverse and sagittal view (shown in Figure 5). Based on time gating approach, the temporal responses of photon migration in human brain was made. Figure 7 shows the trajectories of all 2 × 10^{7} photons in the transverse and sagittal cases of brain structure. Obviously, the difference patterns of photon migration can be observed. In the result, the light can reach to white matter layer (~3 cm depth) after 100 psec in transverse views. The result of simulation revealed the photon guiding effect in CSF layer. In our simulation, all photonpassed voxels were recorded for photon migration analysis. The CSF light guiding effect can be easily observed in the movie file. This paper has supplementary downloadable material available of Additional file 1 provided by the authors. This includes a multimedia MOV format movie file, which shows the dynamic photon migration with 800 nm light pulse illumination through the adult brain model. Two crosssectional views are demonstrated as transverse and sagittal. This material is 7.65 MB in size.
Additional file 1: The movies of the dynamics of photon migration in brain models. The movie shows the dynamic photon migration with 800 nm light pulse illumination through the adult brain models. Two crosssectional views are demonstrated as transverse and sagittal. (MOV 8 MB)
Figure 8 demonstrates the paths of detected photons via sourcedetector separation in transverse and sagittal view of individualised model. Ten detectors were placed away from the light source with 110 cm as Figure 5. The red arrow indicates the location of light source and orange one represents each detector. The power of received light was decreasing while sourcedetector separation increasing. According to the result, the optimal sourcedetector separation was chosen between 2 and 4 cm for brain detection in this individualised model. However, Figure 8(c) indicates the longer propagation distance of diffuse photon in sagittal view because of its bigger CSF volume in interhemispheric fissure. Therefore, the behavior of light propagation in human brain strongly depends on individualised structure of brain.
The proposed modeling method can offer multiwavelength illumination. Figure 9 shows the curves of intensity distribution with various sourcedetector separations in this individualised brain model at 690 nm, 780 nm and 830 nm. Figure 9 implies that the propagated photon through brain was absorbed stronger at longer wavelength. Also, the result reveals that the change rate of detected intensity via sourcedetector separation can provide a quantitative analysis to evaluate the brain structure individually.
In this study, the paths of received photons from each layer were recorded. Figure 10 shows the ratios of the backscattered intensities from different layer versus the sourcedetector separation with multiwavelength. Obviously, the signal from the surface (scalp and skull) layer and cerebral cortex layer were crossed at about 3.3 cm of sourcedetector separation. In this result of this individualised model, the backscattered light from the cerebral cortex layer is greater than 50%, while the sourcedetector separation exceeds the crosspoint. Compare this result with Figure 9, the total received intensity was decreases strongly with the sourcedetector separation increasing. Hence, the sourcedetector separation in this individualised case was optimally set as 3.3 cm for NIRS measurement.
Discussion
In this paper, the patientoriented and individualized simulation for brain monitoring by using in vivo MRI data was proposed. An adult brain was modeled in threedimensional timeresolved Monte Carlo simulation to investigate the structural characteristic of individualized brain. The result indicates that the threedimensional model which faithfully represents the realistic and individualized human head from MRI data depends on image processing. In order to detect the characteristic of individualized structural with NIRS measurement, the various sourcedetector separations on human head were simulated dynamically with transverse and sagittal views in an adult brain model. Although the previously studies indicated the light guiding effect occurred in the CSF layer of human head. We described the Monte Carlo method that is capable of performing the penetration analysis with dynamic photon migration movies to show more clearly effects of light guiding by CSF. The penetration of 800 nm light via the CSF is remarkably clear in the Figure 7.
The CSF is a region that is filled with clear low scattering fluid in the head. In other words, the effect of both low scattering and absorption coefficients in the CSF layer reveals a strong effect on light propagation in the head. Thus, the photons propagate longer distance along the CSF layer can be observed clearly in the movie file of Figure 7, especially in sagittal view. In the Figure 8, we observed that the penetration depth in sagittal view is longer than transverse with sourcedetector separation from 3 cm to 10 cm, especially from 6 cm to 10 cm. In other words, the transverse crosssection contains bigger volume of gray and white matter and smaller volume of CSF than the sagittal crosssection that can observe in Figure 8. Besides, the gray and white matters generate absorption and multiscattering that cause shallow penetration in transverse crosssection. On the contrary, the CSF layer provides low extinction of light that can help light propagation longer in sagittal crosssection. The results showed that the spatial sensitivity profile in the head formed unlike the wellknown "banana" shape when the sourcedetector separations are less than 3 cm, which covered uniformly between the source and the detector, and covered the gray matter and even the surface of the white matter. Significant distortions were observed around the cerebral cortex folding. The spatial sensitivity profile penetrated deeper to the brain in the case of expanded CSF. Accordingly, the cerebral cortex folding geometry was suggested to significantly affect the spatial sensitivity profile in human head because it filled with CSF. This was discrepant with the previous finding by loose brain models [39, 40]. However, accurate modeling of the brain structure based on image segmentation process, the effect of the sulcus (filled of CSF) on the spatial sensitivity profile was obvious. Therefore, in most studies using the brain models based on MRI data of the adult brain [13–15, 17, 38, 39], the cerebral cortex folding in the models were suggested to be not exact and large enough to affect the spatial sensitivity profile. In the sagittal section, the photons propagate longer distance along the CSF layer can be observed clearly because the expanded interhemispheric fissure. This finding suggests that the optical method may provide not only functional signal from brain activation but also structural information of brain atrophy with the expanded CSF layer.
The multiwavelength simulations at 690, 780, and 830 nm demonstrate nicely effects of increasing absorption with wavelength and light guiding effect through the cerebral cortex folding.
In previous studies, the sourcedetector separation is usually chosen between 2 and 3 cm in the previous studies. It is still a tradeoff problem between the received intensity and the useful information in NIRS measurement. To our knowledge, this is the first study to provide an individualized modeling method for patientoriented simulation with NIRS measurement. In the Figure 10, the ratio of the received intensity indicates the existence of brain activation signals from the surface of cerebral cortex (surface of gray and white matter). According to the distribution of received intensity versus sourcedetector separation (shown in Figure 9) and the ratio of the received intensity from different layers (shown in Figure 10), our results suggest that the optimal choice of sourcedetector separation for this individualized case is set as 3.3 cm. In this paper, the new contributions are stated as follows:

1.
In previous studies, although the results of Monte Carlo simulation of light propagation in full segmented 3D MRI model of the human head was presented and the code was released for use by other researchers, it was an only one regular brain model. In our study, we provided an efficient and systematic modeling method of individual brain model for patientoriented measurement and analysis. The signaltonoise ratio evaluation and optimal choice of sourcedetector separation for individualized brain may provide more helpful information for NIRS systems design.

2.
Li et. al. indicated that the significant characterization on the visible Chinese human model was significantly stronger than that on the MRI model [11]. Additionally, we clarified and proved that the threedimensional model which faithfully represents the realistic human head from MRI data depends on image processing.

3.
Currently, most Monte Carlo simulations have been suited to a single wavelength. However, the NIRS system usually applies multiwavelength. In our study, we reformed the Monte Carlo simulation for multiwavelength sources to approach the practical NIRS measurement.
Conclusions
In conclusion, the threedimensional timeresolved brain modeling method approaches the realistic human brain that provides useful information for NIRS systematic design and calibration for individualized case with prior MRI data. Besides, NIRS, with its advantages, could be a useful research tool for the diagnosis of patientoriented in the near future.
References
 1.
Strangman G, Boas DA, Sutton JP: Noninvasive neuroimaging using nearinfrared light. Biol Psychiatry 2002, 52: 679–693. 10.1016/S00063223(02)015500
 2.
Wolf M, Wolf U, Choi JH, Gupta R, Safonova LP, Paunescu LA, Michalos A, Gratton E: Functional FrequencyDomain NearInfrared Spectroscopy Detects Fast Neuronal Signal in the Motor Cortex. NeuroImage 2002, 17: 1868–1875. 10.1006/nimg.2002.1261
 3.
Gratton G, Brumback CR, Gordon BA, Pearson MA, Low KA, Fabiani M: Effects of measurement method, wavelength, and sourcedetector distance on the fast optical signal. NeuroImage 2006, 32: 1576–1590. 10.1016/j.neuroimage.2006.05.030
 4.
Medvedeva AV, Kainerstorferb J, Borisova SV, Barbourc RL, VanMetera J: Eventrelated fast optical signal in a rapid object recognition task: Improving detection by the independent component analysis. Brain Res 2008, 1236: 145–158.
 5.
Hillman EMC: Optical brain imaging in vivo: techniques and applications from animal to man. J Biomed Opt 2007, 12: 051402. 10.1117/1.2789693
 6.
Obrig H, Wenzel R, Kohl M, Horst S, Wobst P, Steinbrink J, Thomas F, Villringer A: Nearinfrared spectroscopy: does it function in functional activation studies of the adult brain? Int J Psychophysiol 2000, 35: 125–142. 10.1016/S01678760(99)000483
 7.
Boas DA, Gaudette T, Strangman G, Cheng X, Marota JJA, Mandeville JB: The accuracy of near infrared spectroscopy and imaging during focal changes in cerebral hemodynamics. NeuroImage 2001, 13: 76–90.
 8.
Gebhart SC, Lin WC, Jansen AM: In vitro determination of normal and neoplastic human brain tissue optical properties using inverse addingdoubling. Phys Med Biol 2006, 51: 2011–2027. 10.1088/00319155/51/8/004
 9.
Zhao H, Tanikawa Y, Gao F, Onodera Y, Sassaroli A, Tanaka K, Yamada Y: Maps of optical differential pathlength factor of human adult forehead, somatosensory motor and occipital regions at multiwavelengths in NIR. Phys Med Biol 2002, 47: 2075–2093. 10.1088/00319155/47/12/306
 10.
Bevilacqua F, Piguet D, Marquet P, Gross JD, Tromberg BJ, Depeursinge C: In vivo local determination of tissue optical properties: applications to human brain. Appl Optics 1999, 38: 4939–4950. 10.1364/AO.38.004939
 11.
Li T, Gong H, Luo Q: Visualization of light propagation in visible Chinese human head for functional nearinfrared spectroscopy. JBO 2011, 16(045001):1–6.
 12.
Okada E, Delpy DT: Nearinfrared light propagation in an adult head model. I. Modeling of lowlevel scattering in the cerebrospinal fluid layer. Appl Optics 2003, 42: 2906–2914. 10.1364/AO.42.002906
 13.
Fukui Y, Ajichi Y, Okada E: Monte Carlo prediction of nearinfrared light propagation in realistic adult and neonatal head models. Appl Optics 2003, 42: 2881–2887. 10.1364/AO.42.002881
 14.
Mudra R, Nadler A, Keller E, Niederer P: Analysis of nearinfrared spectroscopy and indocyanine green dye dilution with Monte Carlo simulation of light propagation in the adult brain. J Biomed Opt 2006, 11: 044009. 10.1117/1.2341652
 15.
Xu Y, Graber HL, Barbour RL: Image correction algorithm for functional threedimensional diffuse optical tomography brain imaging. Appl Optics 2007, 46: 1693–1704. 10.1364/AO.46.001693
 16.
Okada E, Delpy DT: Nearinfrared light propagation in an adult head model. II. Effect of superficial tissue thickness on the sensitivity of the nearinfrared spectroscopy signal. Appl Optics 2003, 42: 2915–2922. 10.1364/AO.42.002915
 17.
Boas DA, Dale AM: Simulation study of magnetic resonance imagingguided cortically constrained diffuse optical tomography of human brain function. Appl Optics 2005, 44: 1957–1968. 10.1364/AO.44.001957
 18.
Ogoshi Y, Okada E: Analysis of light propagation in a realistic head model by a hybrid method for optical brain function measurement. Opt Rev 2005, 12: 264–269. 10.1007/s100430050264y
 19.
Heiskala J, Nissilä I, Neuvonen T, Järvenpää S, Somersalo E: Modeling anisotropic light propagation in a realistic model of the human head. Appl Optics 2005, 44: 2049–2057. 10.1364/AO.44.002049
 20.
Hayashi T, Kashio YO, Okada E: Hybrid Monte Carlodiffusion method for light propagation in tissue with a lowscattering region. Appl Optics 2003, 42: 2888–2896. 10.1364/AO.42.002888
 21.
Hoshi Y, Shimada M, Sato C, Iguchi Y: Reevaluation of nearinfrared light propagation in the adult human head: implications for functional nearinfrared spectroscopy. J Biomed Opt 2005, 10: 064032. 10.1117/1.2142325
 22.
Diamond SG, Huppert TJ, Kolehmainen V, Franceschini MA, Kaipio JP, Arridge SR, Boas DA: Dynamic physiological modeling for functional diffuse optical tomography. NeuroImage 2006, 30: 88–101. 10.1016/j.neuroimage.2005.09.016
 23.
Firbanky M, Arridgez SR, Schweigery M, Delpy DT: An investigation of light transport through scattering bodies with nonscattering regions. Phys Med Biol 1996, 41: 767–783. 10.1088/00319155/41/4/012
 24.
Wolf M, Keel M, Dietz V, von Siebenthal K, Bucher HU, Baenziger O: The influence of a clear layer on nearinfrared spectrophotometry measurements using a liquid neonatal head phantom. Phys Med Biol 1999, 44: 1743–1753. 10.1088/00319155/44/7/313
 25.
Dehghani H, Delpy DT, Arridge SR: Photon migration in nonscattering tissue and the effects on image reconstruction. Phys Med Biol 1999, 44: 2897–2906. 10.1088/00319155/44/12/303
 26.
Oshe S: Fronts Propagating with Curvature Dependent Speed: Algorithms Based on HamiltonJacobi Formulations. J Comput Phys 1998, 79: 12–49.
 27.
Raviv TR, Kiryati N, Sochen N: Segmentation by Level Sets and Symmetry. IEEE Computer Society Conference on Computer Vision and Pattern Recognition Volume 1 CVPR06 2006, 00: 1015–1022.
 28.
Leventon ME, Faugeras O, Grimson WEL, Wells WM: Level Set Based Segmentation with Intensity and Curvature Priors. IEEE Workshop on Math Met Biomed Imag Anal 2000, 00: 4–11.
 29.
Li C, Xu C, Gui C, Fox MD: Level set evolution without reinitialization: a new variational formulation. IEEE Comp Soc Conf Comp Vis Pattern Recogn 2005, 1: 430–436.
 30.
Hojjatoleslami SA, Kittler J: Region growing: a new approach. IEEE T Image Process 1998, 7: 1079–1084. 10.1109/83.701170
 31.
Tremeau A, Borel N: A region growing and merging algorithm to color segmentation. Pattern Recognition 1997, 30: 1191–1203. 10.1016/S00313203(96)001471
 32.
Ashburner J, Friston KJ: Unified segmentation. NeuroImage 2005, 26: 839–851. 10.1016/j.neuroimage.2005.02.018
 33.
Strangman G, Franceschini MA, Boas DA: Factors affecting the accuracy of nearinfrared spectroscopy concentration calculations for focal changes in oxygenation parameters. NeuroImage 2003, 18: 865–879. 10.1016/S10538119(03)000211
 34.
Lee CK, Sun CW, Lee PL, Lee HC, Yang CC, Jiang CP, Tong YP, Yeh TC, Hsieh JC: Study of photon migration with various sourcedetector separations in nearinfrared spectroscopic brain imaging based on threedimensional Monte Carlo modelling. Opt Express 2005, 13: 8339–8348. 10.1364/OPEX.13.008339
 35.
Wang LH, Jacques SL, Zheng LQ: MCML  Monte Carlo modeling of photon transport in multilayered tissues. Comput Meth Prog Bio 1995, 47: 131–146. 10.1016/01692607(95)01640F
 36.
Kirillin M, Meglinski I, Kuzmin V, Sergeeva E, Myllylä R: Simulation of optical coherence tomography images by Monte Carlo modeling based on polarization vector approach. Opt Express 2010, 18: 21714–21724. 10.1364/OE.18.021714
 37.
Churmakov DY, Meglinski IV, Greenhalgh DA: Influence of refractive index matching on the photon diffuse reflectance. Phys Med Biol 2002, 47: 4271–4285. 10.1088/00319155/47/23/312
 38.
Meglinsky IV, Matcher SJ: Modelling the sampling volume for skin blood oxygenation measurements. Med Biol Eng Comput 2001, 39: 44–50. 10.1007/BF02345265
 39.
Mansouri C, Huillier JPL, Kashou NH, Humeau A: Depth sensitivity analysis of functional nearinfrared spectroscopy measurement using threedimensional Monte Carlo modellingbased magnetic resonance imaging. Lasers Med Sci 2010, 25: 431–438. 10.1007/s1010301007544
 40.
Boas DA, Culver JP, Stott JJ, Dunn AK: Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head. Opt Express 2001, 10: 159–170.
 41.
Bohren CF, Huffman DR: Absorption and Scattering of Light by Small Particles. John Wiley & Sons; 1983.
 42.
Okada E, Firbank M, Schweiger M, Arridge SR, Cope M, Delpy DT: Theoretical and experimental investigation of nearinfrared light propagation in a model of the adult head. Appl Opt 1997, 36: 21–31.
 43.
Meglinsky IV, Matcher SJ: Analysis of the spatial distribution of the detector sensitivity in a multilayer randomly inhomogeneous medium with strong light scattering and absorption by the Monte Carlo method. Opt Spectrosc 2001, 91: 692–697.
Acknowledgements
This research was supported by National Taiwan University YongLin Biomedical Engineering Center and under grants FB0022, National Taiwan University and under grants 10R809214, the National Science Council of Taiwan and under grants NSC 1002221E010004, NSC 1002622E010003CC3, NSC 1002627E010001, a grant from Ministry of Education, Aim for the Top University Plan in National YangMing University and Yen Tjing Ling Medical Foundation.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
CC (first author) contributed in the theoretical model, proposal of the method, and writing of the manuscript. YL, CC (third author), YH and TL contributed equally in the analysis of algorithms. CS conceived of the study, and participated in its design and coordination and helped to writing the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Chuang, C., Lee, Y., Chen, C. et al. Patientoriented simulation based on Monte Carlo algorithm by using MRI data. BioMed Eng OnLine 11, 21 (2012). https://doi.org/10.1186/1475925X1121
Received:
Accepted:
Published:
Keywords
 Patientoriented simulation
 Timeresolved Monte Carlo
 Brain modeling
 Spatial sensitivity profile