Document Detail

Assessment of the frequency-domain multi-distance method to evaluate the brain optical properties: Monte Carlo simulations from neonate to adult.
Jump to Full Text
MedLine Citation:
PMID:  21412461     Owner:  NLM     Status:  PubMed-not-MEDLINE    
Abstract/OtherAbstract:
The near infrared spectroscopy (NIRS) frequency-domain multi-distance (FD-MD) method allows for the estimation of optical properties in biological tissue using the phase and intensity of radiofrequency modulated light at different source-detector separations. In this study, we evaluated the accuracy of this method to retrieve the absorption coefficient of the brain at different ages. Synthetic measurements were generated with Monte Carlo simulations in magnetic resonance imaging (MRI)-based heterogeneous head models for four ages: newborn, 6 and 12 month old infants, and adult. For each age, we determined the optimal set of source-detector separations and estimated the corresponding errors. Errors arise from different origins: methodological (FD-MD) and anatomical (curvature, head size and contamination by extra-cerebral tissues). We found that the brain optical absorption could be retrieved with an error between 8-24% in neonates and infants, while the error increased to 19-44% in adults over all source-detector distances. The dominant contribution to the error was found to be the head curvature in neonates and infants, and the extra-cerebral tissues in adults.
Authors:
Mathieu Dehaes; P Ellen Grant; Danielle D Sliva; Nadège Roche-Labarbe; Rudolph Pienaar; David A Boas; Maria Angela Franceschini; Juliette Selb
Publication Detail:
Type:  Journal Article     Date:  2011-02-11
Journal Detail:
Title:  Biomedical optics express     Volume:  2     ISSN:  2156-7085     ISO Abbreviation:  Biomed Opt Express     Publication Date:  2011  
Date Detail:
Created Date:  2011-03-17     Completed Date:  2011-07-14     Revised Date:  2014-09-12    
Medline Journal Info:
Nlm Unique ID:  101540630     Medline TA:  Biomed Opt Express     Country:  United States    
Other Details:
Languages:  eng     Pagination:  552-67     Citation Subset:  -    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Grant Support
ID/Acronym/Agency:
R01 HD042908/HD/NICHD NIH HHS; R01 HD042908-07/HD/NICHD NIH HHS; R01 HD042908-08/HD/NICHD NIH HHS; R01 HD042908-09/HD/NICHD NIH HHS
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Biomed Opt Express
Journal ID (publisher-id): BOE
ISSN: 2156-7085
Publisher: Optical Society of America
Article Information
Download PDF
©2011 Optical Society of America
open-access:
Received Day: 2 Month: 11 Year: 2010
Revision Received Day: 3 Month: 2 Year: 2011
Accepted Day: 4 Month: 2 Year: 2011
Electronic publication date: Day: 11 Month: 2 Year: 2011
collection publication date: Day: 1 Month: 3 Year: 2011
pmc-release publication date: Day: 11 Month: 2 Year: 2011
Volume: 2 Issue: 3
First Page: 552 Last Page: 567
ID: 3047361
PubMed Id: 21412461
Publisher Id: 137543
DOI: 10.1364/BOE.2.000552

Assessment of the frequency-domain multi-distance method to evaluate the brain optical properties: Monte Carlo simulations from neonate to adult
Mathieu Dehaes1*
P. Ellen Grant12
Danielle D. Sliva1
Nadège Roche-Labarbe2
Rudolph Pienaar1
David A. Boas2
Maria Angela Franceschini2
Juliette Selb2
1Fetal-Neonatal Neuroimaging & Development Science Center, Children's Hospital Boston and Harvard Medical School, Boston, Massachusetts 02115, USA
2Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts 02129, USA
*mathieu.dehaes@childrens.harvard.edu

1. Introduction

Near-infrared spectroscopy (NIRS) measures the attenuation of diffuse near-infrared light (650–950 nm) propagating through biological tissues at two or more wavelengths, allowing for the characterization of the two main chromophores in this spectral window: oxygenated hemoglobin (HbO) and reduced hemoglobin (HbR). Therefore, NIRS is highly sensitive to tissue vascularization and oxygenation, making it an ideal tool for cerebral monitoring and imaging. Applications to functional brain imaging have spanned a diverse range of studies, including visual [1–3], sensori-motor [4], auditory [5, 6], and cognitive [7, 8] tasks. While in this context of functional imaging NIRS measures the hemodynamic changes associated with brain activation, the determination of the absolute optical absorption of the brain, directly related to baseline hemoglobin concentrations, is a more challenging application [9, 10]. Such quantitative absolute measurements are of considerable interest for longitudinal and cross-sectional studies, allowing comparison of cerebral oxygenation and blood volume between subjects, as well as monitoring within a subject over time, opening the door to a number of diagnostic and treatment monitoring applications [11, 12].

The most commonly used modality for functional NIRS, the continuous-wave (CW) approach, is not suited for absolute characterization of tissue optical properties, because it cannot distinguish between the absorption μa and scattering μ′s of a medium [13]. Some publications have nonetheless presented CW estimates of baseline optical properties of tissues [14–16]. However they either report a tissue oxygenation index (TOI) [14, 16], which is the product μ′s × μa, or they require that oxygen saturation be varied within a physiological range through such manipulations as hypoxic events [15] which ethically cannot be performed in infants.

Unlike CW measurements, NIRS in the time-domain (TD) and frequency-domain (FD) allows for the characterization of the baseline μa and μ′s of a biological tissue [17–20] without inducing changes in saturation, or making assumptions about μ′s. Pulsed (TD method) or radiofrequency modulated (FD method) light is injected into the medium and its detection provides information about mean time of flights of photons (TD method) or phase shift (FD method), both related to μ′s. The commonly used FD multi-distance (FD-MD) method takes advantage of the phase shift and attenuation of modulated light at several distances to retrieve the baseline optical properties of the medium assumed homogeneous [18]. These TD and FD methods have been validated with simulations, phantom experiments and animal studies [21–28]. In vivo applications of baseline measurements of the optical properties of the human brain include studies in premature and term neonates and infants [10–12, 29, 30], as well as adults [28, 31–33]. For example, our group has applied the FD-MD method on a large scale to the study of healthy brain development in term and pre-term infants during the first year of life [10, 11], as well as to the assessment of brain injury in a neonatal population [12].

With the exception of a few studies, most approaches rely on a homogeneous semi-infinite model of the head because of its simplicity. However this lack of realism introduces error in the recovered optical properties. Two- [22, 24, 28, 34] or three-layer [35] geometries have been suggested instead, as well as more realistic brain structure obtained from magnetic resonance imaging (MRI) segmentation [36], but practical applications remain challenging. In newborns the measurement can be influenced by the presence of a thick low-scattering CSF layer, as well as the small head size and the stronger curvature of the head's surface. In adults, the main challenge is the thickness of the extra-cerebral layers, which can reach more than 1 cm. No systematic study of these effects has been published to our knowledge. Validation studies are difficult to perform because there is no access to actual cerebral absorption values in humans.

In the present study, our aim was to use Monte Carlo simulations to systematically investigate the accuracy of the FD-MD method for retrieving brain absorption coefficients using a semi-infinite homogeneous model with a range of age groups, from the newborn to the adult.


2. Methods
2.1. Simulated measurements: Monte Carlo approach

Realistic data were generated using a 3D Monte Carlo (MC) code [37], modeling the photon propagation in head segmentations obtained from MRI scans.

2.1.1. Head models

We simulated data sets for four different ages: a term newborn of 38 weeks gestational age, 6 and 12 month (mo.) old infants, and an adult. A realistic 3D head model was obtained from each data set by segmentation of an MR T1 scan into four tissue types (see Fig. 1): white matter (WM), gray matter (GM), cerebrospinal fluid (CSF) and a superficial layer (superf) comprising scalp and skull. The voxel size was 1.0 mm in all directions for all head models. Figure 1 shows 3D views and axial slices of the segmented heads at each age. Gray and white matter were subsequently combined into a single “brain” tissue type with identical optical properties [38].

The effect of extra-cerebral inhomogeneities was investigated by simulating a blood vessel in the scalp of the adult. Located parallel to the optical probe at a depth of 4 to 6 mm from the surface of the skin, the vessel was 30 mm long and had a diameter of 2 mm.

2.1.2. Optical properties

The optical properties of the scalp, skull and brain vary widely in the literature depending on the method of measurement [28–33]. Ex-vivo and post-mortem values are bound to differ significantly from in-vivo data due to changes in blood content. To simulate realistic absorption values, we considered a wide range of hemoglobin concentrations spanning those found in the literature (between 40 and 120 μM total hemoglobin) and converted them into absorption coefficients using tabulated values for the extinction coefficients of HbO and HbR [39]. This hemoglobin range yielded brain absorption coefficient between 0.05 and 0.3cm−1, while superficial absorption was set to 0.1 and 0.2cm−1. The absorption coefficient of the scalp blood vessel described in section 2.1.1 was set to 5cm−1, a typical value for whole blood absorption in the NIR spectrum [40]. Some studies report higher scattering in the brain than in the superficial tissues [31–33], while others present the opposite result [28], accompanied with inter-individual variability [30, 32]. For simplicity, the reduced scattering coefficient μ′s was fixed at 10cm−1 in the brain. The value of μ′s in the superficial layer was set successively to 10cm−1 and 5cm−1, to evaluate the effect of a spatial variation between tissues. The CSF was either considered part of the brain, with the same absorption and scattering coefficients, or, more realistically, as a relatively clear medium with low absorption and scattering coefficients μa = 0.04cm−1 and μ′s = 0.1cm−1. The refractive index n was set to 1.3, and the anisotropy coefficient g to 0.01 in all tissue types (the MC code does not accept a null value for g; see Table 1).

2.1.3. Probe geometry

Measurements were modeled from one source seen by four detectors positioned in a straight line, with 5 mm between successive detectors, as in a common FD-MD configuration. To test the effect of different groups of source-detector separations (SDs), a total of 10 detectors located every 5 mm were used and subsequently analyzed in subsets of four detectors, yielding 7 different probe configurations. The probe was positioned over the right parietal area (right and left parietal for the newborn) where the head surface is relatively flat (Fig. 1). The 1 mm voxelization of the heads creates a degree of “stair” effect in the curved regions, leading to artefacts in the simulated data, especially at short distances. To reduce this effect, 5 different probe locations within a small neighborhood (2 mm radius) were simulated and averaged together.

2.1.4. Monte Carlo simulations

The MC code has been described in details by Boas et al. [37]. FD measurements were simulated at a modulation frequency of 110 MHz, and the phase and modulated intensity were computed at each detector. 500 million photon packets were launched for each simulation, which took between 8 and 16 hours depending on the number of voxels included in the head.

2.2. Noise modeling

The large number of photons propagated in the MC simulations yields measurements with a very high signal to noise ratio, unrepresentative of actual experimental data. Therefore, measurements generated from MC simulations were modified by the addition of a Gaussian random noise. From each MC simulation data set of phase and intensity, we created 150 noisy measurements by adding 150 different values of phase and intensity noises. The variance of the phase and intensity noises were based on realistic values obtained experimentally with a FD instrument (Imagent, ISS Inc.) in neonates [10, 11] (noise data not published). The noise to signal ratio increases as the source-detector separation increases, for both intensity (1.6% to 4.6%) and phase (0.02% to 3.7%) from 0.5 to 5 cm, respectively. The standard deviation of the retrieved absorption coefficients over the 150 values of noise gives a realistic estimate of uncertainty due to both instrumental and physiological noises.

2.3. Optical properties recovery

The estimation of brain optical properties from simulated data was conducted using the homogeneous FD-MD method described previously [9, 18]. Briefly, in a semi-infinite homogeneous medium, the phase φ and natural logarithm of the modulated intensity of light IAC vary linearly with source-detector separation d such that

[Formula ID: e01]
ln(d2IAC)=dSAC+CACandφ=dSφ+Cφ
(1) 

where SAC and Sφ are the slopes and CAC and Cφ are the y-intercepts of the linear regressions. Absorption and reduced scattering coefficients of the homogeneous model can then be estimated by

[Formula ID: e02]
μa=ω2ν(SφSAC−SACSφ)andμs′=SAC2−Sφ23μa−μa
(2) 

where ω is the angular modulation frequency of the source intensity and ν is the speed of light in the tissue.

2.4. Calibration method

The intensity and phase signals were calibrated prior to recovery. The calibration method implemented [18] is the same as that used in experimental data acquisition [10, 11]. It calculates corrective factors for the DC, AC, and phase of the light source [41], based on measurements in a phantom block of known optical properties. In this study, we modeled the measurements of the calibrated block with MC simulations in a large homogeneous slab imitating a semi-infinite medium. The dimensions of the slab were 18 × 18 × 10cm3 and its optical properties set to μa = 0.1cm−1, μ′s = 10cm−1, g = 0.01, and n = 1.3. The corrective calibration factors (IcalibAC ,φcalib) were calculated so that the slopes of intensity and phase simulated from the slab (IslabAC ,φslab) matched the theoretical ones:

[Formula ID: e03]
IACcalib=IACtheoIACslabwithIACtheo=exp(dSACtheo)d2
(3) 

and

[Formula ID: e04]
φcalib=φtheo−φslabwithφtheo=dSφtheo
(4) 

where the expressions StheoAC and Stheoφ are detailed in [42]. These calibration factors were then applied to noisy intensity and phase signals from the simulated head data. Finally, Eq. (2) were used to recover the values of optical properties.

We show in Fig. 2 the effect of the calibration procedure in the homogeneous slab itself. In particular, the calibration helps reducing large errors due to the FD-MD approximations at short SDs. However the effect of the absorption coefficient of the calibration block appears clearly. When the absorption coefficient of the measured sample is close to that of the calibration block (0.1cm−1, green set), the calibration process improves the results. When the absorption of the medium is much higher than that of the calibration block (e.g. 0.3cm−1, yellow set) the calibration process results in an underestimation of the absorption in the medium.


3. Results
3.1. Homogeneous head model

To separate the errors due to the diffusion approximation, finite geometry, and curvature from the errors due to the layered structure of the head, we first treated the head as a homogeneous medium by setting equal the superficial and brain absorptions, i.e. μa,superf = μa,brain = [0.05,0.06, . . . ,0.29,0.3]cm−1. In Fig. 3, each graph presents the retrieved brain absorption for the 7 sets of SDs. The adult model shows better estimate than the young models. This is due to the stronger head curvature for the newborn and infants. In addition, all estimations at 5–20 mm are less accurate because of the diffusion approximation utilized in the FD-MD formulation [41].

In Fig. 4, retrieved absorption coefficients and standard deviation due to experimental noise are shown for three specific cases when μa,superf = μa,brain = [0.05,0.15,0.3]cm−1. While the retrieved values improve with longer SDs, the uncertainty increases due to decreasing signal to noise ratio. In addition, the underestimation of retrieved absorption coefficient for μa = [0.15,0.3]cm−1 is due to the unmatched calibration medium (calibration μa = 0.1cm−1).

Figure 5 presents averages and standard deviations of the retrieved μa,brain over 26 combinations of μa,superf = μa,brain = [0.05,0.06, . . . ,0.29,0.3]cm−1. In newborn and infants, the errors vary between 14–21% and 7–20%, respectively, and decreases to 3–14% for the adult head model. At very short distances (5–20 mm), the error is large for all ages (15–20%) due to the non-validity of the diffusion approximation. At long distances, the error decreases down to less than 5% for the adult head, while the stronger curvature of the head is responsible for additional errors in the case of newborn and infants.

3.2. Two-tissue head model

Figure 6 depicts the evolution of the recovered brain absorption values as a function of true absorption for the 7 groups of SDs in all head models (row-wise) treated as a two-tissue structure. The superficial absorption coefficient μa,superf was set successively to 0.1 and 0.2 cm−1 (column-wise). In each sub-figure, the dotted line corresponds to the true brain absorption, while the dash-dotted line identifies the superficial absorption. Solid curves represent the 7 groups of SDs.

The effect of the superficial layer is to pull the recovered absorption value of the brain towards the superficial absorption value, especially at short SDs and for the adult case. Because of the specific sets of optical properties we chose, this selection leads to a generally weaker than true value μa,brain, especially for high brain absorption coefficients. The three young head models present similar brain absorption estimations, mostly underestimated, probably due to the head's curvature. The adult head presents a less precise estimation including over- and underestimations, due to the large thickness of the superficial layer which contaminates the signal especially at short SDs.

Figure 7 shows the average error at each age for all SDs groups, over 52 absorption combinations covering the whole range of coefficients for the brain (μa,brain between 0.05 and 0.3 cm−1 in steps of 0.01 cm−1) and the superficial layer (μa,superf =0.1 or 0.2 cm−1).

For the newborn, the optimal SDs group is 10–25 mm and shows a relative error of 12±6%. The 6 and 12 month old infants present a similar trend: the relative error is about 10±5% for all SDs groups, except for the two shortest (5–20 and 10–25 mm). For both infants, the optimal SDs group can be selected from the five longest SDs (from 15–30 to 35–50 mm), since the error is minimized and shows similar values. The adult head shows the highest relative error for almost all SDs groups (with certain combinations of optical properties yielding errors of more than 80%). The optimal SDs group is reached at 35–50 mm with a corresponding error of 19±13%.

The effect of the scattering coefficient was investigated by modifying the value in the superficial layer to μ′s,superf = 5cm−1. The relative error decreases in the newborn (ranging from 6 to 19%) and infant models (2–10%) for all SDs except the shortest group (5–20 mm), and slightly increases in the adult head (from 26 to 45%) for almost all groups (except for 10–25 and 15–30 mm). Again, this is probably due to the thicker superficial layer in the adult model that causes a stronger error in μ′s translating through cross-talk in higher error in μa.

The same analysis was performed with the adult model in which a large vessel was introduced. Relative errors and standard deviations [%] for the 7 SDs (from 5–20 to 35–50 mm) are the following: 51.6 ± 58.2; 48.3 ± 53.9; 35.4 ± 33.4; 26.3 ± 18.1; 23.7 ± 14.7; 22.0 ± 14.5; 19.2 ± 13.2. The error is higher than the case without inhomogeneity in the first 3 groups and then similar for the next 4 groups. The contamination of the shortest estimations is related to the length of the vessel (3 cm) which was located around the beginning of the probe where the light source is positioned. This type of inhomogeneity can be detected experimentally, and data automatically discarded, because of the corresponding low correlation factor r2 of the linear fitting process from Eq. (1).

3.3. Three-tissue head model

We investigated the influence of the CSF structure by setting its optical properties to low absorption and scattering values: μa = 0.04cm−1 and μ′s = 0.1cm−1. The two hemispheres (left and right parietal) of the newborn were studied, because the right hemisphere possessed a thick “pocket” of CSF located underneath the optical probe (white arrow on the axial slice of the second row in Fig. 8), while the left hemisphere possessed a “normal” amount of fluid.

Figure 8 shows retrieved μa,brain with respect to true values for the 7 groups of SDs. The first two columns correspond to μa,superf set to 0.1 and 0.2cm−1, respectively. The third column shows an axial slice of the left and right newborn (NB) hemispheres. Each graph was created in the same way as in Fig. 6. Distortion of the value of μa,brain by the superficial layer is similar to the effect observed when ignoring CSF, with the difference that the low absorption of CSF pulls the retrieved μa,brain to lower values than in the previous case.

For the newborn, the estimations computed at shorter SDs in the left hemisphere are less distorted than those estimated in the right hemisphere. This is due to the large volume of CSF in the vicinity of the optical probe. The clear pocket of CSF contributes to modification of the signal through the combined effects of very low absorption and scattering. In general, short SDs groups, from 5–20 to 15–30 mm, are poorly estimated for all head models.

Figure 9 shows the averages of the relative errors and standard deviations in percentage [%] of μa,brain for all 52 combinations of optical properties of brain and superficial tissues. The error rises to a maximum of 51% for the newborn right hemisphere, mostly due to the effect of the CSF pocket. For the newborn left hemisphere, the optimal set of SDs is 20–35 mm and the error is 19±5%. Decreasing error with increasing SDs is still true for the adult model to a slightly higher extent. The distinction between 6 and 12 month old infants is more visible than in the case without CSF. The 12 month old infant shows less accurate estimations for almost all SDs because of the thicker superficial layer. In both cases, optimal SDs are between 20–35 and 30–45 mm, where the errors is less than 16±10%. For the adult, the longest SDs group (35–50 mm) is needed to minimize the error to a value of 23 ± 15%.

The influence of the CSF was also studied for a superficial scattering coefficient of 5 cm−1. For all head models, the error at 5–20 mm increases. For the newborn, both hemispheres show lower errors (ranging from 7 to 45%) with the lower scattering coefficient (except for 25–40 mm in the right hemisphere: error increased to 13%). In the case of the 12 mo. old infant, the error increases from 15 to 61% (except for 10–25 and 15–30 mm: down to 30 and 16%, respectively) while decreases to a minimum of 6 and 25% in the 6 mo. old infant and adult, respectively (data not shown).


4. Discussion
4.1. Effect of the curvature

Generally, the estimations computed at longer distances are strongly affected by head curvature in young subjects. The fitted slopes from Eq. (1) deviate from the theoretical slopes computed with the flat calibration slab. In practice, measurements on a curved region of the head lead to a non-linear variation of the phase and intensity with source-detector separation. This error is easily visible in the raw data and such data can be objectively discarded, as performed in [10–12].

4.2. Effect of the superficial layer

When modeling the head with two layers, the adult model shows larger errors for almost all SDs groups. This is due to the thicker layer of extra-cerebral tissue in adults. Even at the optimal SD separations, the brain absorption is still estimated with 19% error, compromising the validity of longitudinal or cross-sectional studies. Implementing more complex models of the head that incorporate the layered structure more realistically can help solve this issue.

4.3. Effect of the CSF layer

In newborn, we studied one location where a thick region of CSF was located directly underneath the optical probe, hindering the recovery of brain optical properties at short distances (from 5–20 to 15–30 mm). This lateral “CSF pocket” could have been produced by the positioning of the head during the MR scan. At this age, the skull is very flexible and the CSF can be distributed according to the pressure applied on the surface of the head. Even though the skull layer is very thin in neonates, facilitating the measurement of brain absorption coefficients, the relatively thick CSF layer can affect the performance of the method. For the adult head model, results with and without CSF were very similar. The 6 and 12 month old infants often yield more accurate estimates than the newborn, with and without the optical properties of the CSF. This is probably due to the description of the CSF, the size of the superficial layer (scalp and skull) and the head size of the infant models. This suggests that the dominant error is not related to the CSF compartment.

4.4. Effect of the scattering coefficient

Surprisingly, decreasing the scattering coefficient of the extra-cerebral tissue from 10 to 5 cm−1 produced better results for the newborn and infants, which can probably be explained by compensatory errors. On the contrary, it deteriorated the estimation for the adult at almost all SDs groups. A more extended study of this effect should yield more robust conclusions.

4.5. Large vessel, hair obstruction and extra-cerebral inhomogeneities

The presence of extra-cerebral inhomogeneities such as vessel or hair is generally uncontrolled and their locations are unknown except if angiography data are accessible. We showed how such strongly absorbing local structure can affect the retrieved absorption, but also how it deteriorates the linear fit of the slope. In practice, such data are easily rejected based on objective criteria employed in experimental studies [10–12], and the probe is slightly moved to avoid the inhomogeneities.

4.6. Derived physiological variables

Computation of baseline hemoglobin concentrations from the retrieved absorption coefficients allows clinical investigations [10–12] and derivation of physiological variables such as cerebral blood volume (CBV) and oxygen saturation (StO2). For the typical wavelengths of 690 and 830 nm, the propagation of the error from retrieved absorption to hemoglobin concentration was computed [43] for the newborn and the adult using specific extinction coefficients [39]. For the newborn, the typical error of 12% in absorption propagates into an error of 17, 20 and 11% for HbO, HbR and HbT, respectively. For the adult, the minimal error of 19% corresponds to 27, 31 and 17% for HbO, HbR and HbT, respectively. The propagation of the error for CBV has the same rate than HbT while the error for StO2 propagates to 8% and 13% for the newborn and adult, respectively.

4.7. Error arising from the diffusion approximation

The FD-MD method is based on the diffusion approximation and its formulation relies on additional approximations [41]. These approximations lead to errors in the recovered optical properties, even for a large homogeneous slab and especially at short SDs (Fig. 2). Typically, the recovered μa for a homogeneous slab is overestimated by 20% at short SDs (5–20 mm), but only 5% at large SDs (from 25–40 to 35–50 mm) when no calibration is applied: this method favors long distance estimations because it neglects terms of order greater or equal to 1 in the Taylor development [41]. Adequate calibration of the data reduces this error.

4.8. Error arising from the finite geometry

The FD-MD method also considers that the light is propagating in a semi-infinite medium. This assumption is less true for the newborn model since the head size is markedly smaller than the adult (Fig. 1). The error related to the finite geometry is higher in the newborn estimations compared to the three other head models.

4.9. Limitations

All head models originated from anatomical MRI data sets and were created with a semi-automatic tissue segmentation procedure. This method enabled us to simulate light propagation with Monte Carlo methods in a voxel-wise approach. However, since each voxel has a discrete 3D resolution of 1mm3 for each head model, the boundaries of all segmented tissues are represented non-uniformly, especially curved layers. The number of photons detected on the skin depends strongly on local geometry, i.e. on the position of the detector with respect to the voxelization defining the local curvature. This local geometry can be slightly irregular across head models, which can then lead to under- or overestimated fluences.

Also, a consequence of voxelization and voxel size is imprecision of the SD separations. The distances d used in the fit procedure are the theoretical distances used with the calibration slab phantom, which correspond to the set [5,10, . . . ,50] mm. However, the “real”, or simulated, distances for all head models are those determined by the voxel difference between the source and detector positions that belong to the skin. These distances can differ by 0.1 to 0.5 mm from those used in the flat slab, where the curvature is null. These artefacts are higher for the newborn than the adult because the neonate head presents more curvature poorly sampled by the 1 mm3 voxels . However, these artefacts are due solely to the simulation process, and are not representative of any errors observable in experimental studies. We reduced these errors by averaging multiple probe locations. These errors in our method set the minimum on the accuracy with which absorption coefficients can be recovered.

4.10. Future work

In the present study, the effect of the heterogeneous anisotropy across the head [44] was not investigated. Future simulations could incorporate realistic tissue-dependent anisotropy factors.

Future simulations at intermediate ages are necessary to determine the age limit at which the FD-MD method becomes less accurate than 10 or 15%. Also, more accurate segmentation algorithms can provide improved head models and thus better estimates for the young population [45]. In addition, the study of the FD-MD method performance in premature neonates will also be crucial in order to characterize the effects of the thick CSF layer in this population.


5. Conclusion

In this study, we investigated the performance of the FD-MD method when measuring the brain optical absorption at different ages from newborn to adult. In the neonate and infants, the error is dominated by the head curvature and can be minimized around 10–15%. At these ages, the effect of the extra-cerebral layer is small and visible only at short SDs. On the contrary, it becomes dominant in the adult case because of the larger thickness of the scalp and skull. In adults, the strong contamination by the superficial layers leads to both under- and overestimations of the brain absorption coefficient, depending on the relative absorption values assigned to the intra- and extra-cerebral tissues. Not only is the mean error much larger (above 20%) than in the younger head models, but there is a large variability of this error across sets of optical properties. This argument strongly suggests the use of a two- [34] or three-layer [35] model to minimize the deterioration produced by the extra-cerebral layers. In all cases, measurements at 5 mm should be avoided because of the diffusion approximation.


Acknowledgments

The authors thank Myong-Sun Choe (CHB) and Lilla Zollei (MGH). This work was supported by the NIH/NICHD grants # R01-HD042908 and R21-HD058725, the NIH/NCRR grant # P41-RR14075 and the Fonds Québécois de la Recherche sur la Nature et les Technologies (FQRNT-B3, Canada).


References and links
1. Honda Y.,Nakato E.,Otsuka Y.,Kanazawa S.,Kojima S.,Yamaguchi M. K.,Kakigi R.. , “How do infants perceive scrambled face?: A near-infrared spectroscopic study,” Brain Res.1308, 137–146 (Year: 2010).10.1016/j.brainres.2009.10.04619874803
2. Wilcox T.,Bortfeld H.,Woods R.,Wruck E.,Boas D. A.. , “Hemodynamic response to featural changes in the occipital and inferior temporal cortex in infants: a preliminary methodological exploration,” Dev. Sci.11(3), 361–370 (Year: 2008).10.1111/j.1467-7687.2008.00681.x18466370
3. Karen T.,Morren G.,Haensse D.,Bauschatz A. S.,Bucher H. U.,Wolf M.. , “Hemodynamic response to visual stimulation in newborn infants using functional near-infrared spectroscopy,” Hum. Brain Mapp.29(4), 453–460 (Year: 2008).10.1002/hbm.2041117525986
4. Gibson A. P.,Austin T.,Everdell N. L.,Schweiger M.,Arridge S. R.,Meek J. H.,Wyatt J. S.,Delpy D. T.,Hebden J. C.. , “Three-dimensional whole-head optical tomography of passive motor evoked responses in the neonate,” Neuroimage30(2), 521–528 (Year: 2006).10.1016/j.neuroimage.2005.08.05916246586
5. Kotilahti K.,Nissilä I.,Näsi T.,Lipiäinen L.,Noponen T.,Meriläinen P.,Huotilainen M.,Fellman V.. , “Hemodynamic responses to speech and music in newborn infants,” Hum. Brain Mapp.31(4), 595–603 (Year: 2010).19790172
6. Zaramella P.,Freato F.,Amigoni A.,Salvadori S.,Marangoni P.,Suppjei A.,Schiavo B.,Chiandetti L.. , “Brain auditory activation measured by near-infrared spectroscopy (NIRS) in neonates,” Pediatr. Res.49(2), 213–219 (Year: 2001).10.1203/00006450-200102000-0001411158516
7. Schecklmann M.,Ehlis A. C.,Plichta M. M.,Fallgatter A. J.. , “Functional near-infrared spectroscopy: a long-term reliable tool for measuring brain activity during verbal fluency,” Neuroimage43(1), 147–155 (Year: 2008).10.1016/j.neuroimage.2008.06.03218657624
8. Bortfeld H.,Fava E.,Boas D. A.. , “Identifying cortical lateralization of speech processing in infants using near-infrared spectroscopy,” Dev. Neuropsychol.34(1), 52–65 (Year: 2009).10.1080/8756564080256448119142766
9. Fantini S.,Hueber D.,Franceschini M. A.,Gratton E.,Rosenfeld W.,Stubblefield P. G.,Maulik D.,Stankovic M. R.. , “Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy,” Phys. Med. Biol.44(6), 1543–1563 (Year: 1999).10.1088/0031-9155/44/6/30810498522
10. Franceschini M. A.,Thaker S.,Themelis G.,Krishnamoorthy K. K.,Bortfeld H.,Diamond S. G.,Boas D. A.,Arvin K.,Grant P. E.. , “Assessment of infant brain development with frequency-domain near-infrared spectroscopy,” Pediatr. Res.61(5 Pt 1), 546–551 (Year: 2007).17413855
11. Roche-Labarbe N.,Carp S. A.,Surova A.,Patel M.,Boas D. A.,Grant P. E.,Franceschini M. A.. , “Noninvasive optical measures of CBV, StO(2), CBF index, and rCMRO(2) in human premature neonates’ brains in the first six weeks of life,” Hum. Brain Mapp.31(3), 341–352 (Year: 2010).10.1002/hbm.2086819650140
12. Grant P. E.,Roche-Labarbe N.,Surova A.,Themelis G.,Selb J.,Warren E. K.,Krishnamoorthy K. S.,Boas D. A.,Franceschini M. A.. , “Increased cerebral blood volume and oxygen consumption in neonatal brain injury,” J. Cereb. Blood Flow Metab.29(10), 1704–1713 (Year: 2009).10.1038/jcbfm.2009.9019675563
13. Arridge S. R.,Cope M.,Delpy D. T.. , “The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis,” Phys. Med. Biol.37(7), 1531–1560 (Year: 1992).10.1088/0031-9155/37/7/0051631197
14. Naulaers G.,Morren G.,Van Huffel S.,Casaer P.,Devlieger H.. , “Cerebral tissue oxygenation index in very premature infants,” Arch. Dis. Child. Fetal Neonatal Ed.87(3), F189–F192 (Year: 2002).10.1136/fn.87.3.F18912390989
15. Leung T. S.,Tachtsidis I.,Smith M.,Delpy D. T.,Elwell C. E.. , “Measurement of the absolute optical properties and cerebral blood volume of the adult human head with hybrid differential and spatially resolved spectroscopy,” Phys. Med. Biol.51(3), 703–717 (Year: 2006).10.1088/0031-9155/51/3/01516424590
16. Pichler G.,Grossauer K.,Peichl E.,Gaster A.,Berghold A.,Schwantzer G.,Zotter H.,Müller W.,Urlesberger B.. , “Combination of different noninvasive measuring techniques: a new approach to increase accuracy of peripheral near infrared spectroscopy,” J. Biomed. Opt.14(1), 014014 (Year: 2009).10.1117/1.307619319256702
17. Patterson M. S.,Chance B.,Wilson B. C.. , “Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties,” Appl. Opt.28(12), 2331–2336 (Year: 1989).10.1364/AO.28.00233120555520
18. Fantini S.,Franceschini M.,Maier J.,Walker S.,Barbieri B.,Gratton E.. , “Frequency-domain multichannel optical detector for non-invasive tissue spectroscopy and oximetry,” Opt. Eng.34(1), 32–42 (Year: 1995).10.1117/12.183988
19. Matcher S. J.. , “Closed-form expressions for obtaining the absorption and scattering coefficients of a turbid medium with time-resolved spectroscopy,” Appl. Opt.36(31), 8298–8302 (Year: 1997).10.1364/AO.36.00829818264371
20. Liebert A.,Wabnitz H.,Grosenick D.,Möller M.,Macdonald R.,Rinneberg H.. , “Evaluation of optical properties of highly scattering media by moments of distributions of times of flight of photons,” Appl. Opt.42(28), 5785–5792 (Year: 2003).10.1364/AO.42.00578514528944
21. Ntziachristos V.,Chance B.. , “Accuracy limits in the determination of absolute optical properties using time-resolved NIR spectroscopy,” Med. Phys.28(6), 1115–1124 (Year: 2001).10.1118/1.137367411439481
22. Pifferi A.,Torricelli A.,Taroni P.,Cubeddu R.. , “Reconstruction of absorber concentrations in a two-layer structure by use of multidistance time-resolved reflectance spectroscopy,” Opt. Lett.26(24), 1963–1965 (Year: 2001).10.1364/OL.26.00196318059746
23. Swartling J.,Dam J. S.,Andersson-Engels S.. , “Comparison of spatially and temporally resolved diffuse-reflectance measurement systems for determination of biomedical optical properties,” Appl. Opt.42(22), 4612–4620 (Year: 2003).10.1364/AO.42.00461212916630
24. Martelli F.,Del Bianco S.,Zaccanti G.,Pifferi A.,Torricelli A.,Bassi A.,Taroni P.,Cubeddu R.. , “Phantom validation and in vivo application of an inversion procedure for retrieving the optical properties of diffusive layered media from time-resolved reflectance measurements,” Opt. Lett.29(17), 2037–2039 (Year: 2004).10.1364/OL.29.00203715455772
25. Abrahamsson C.,Svensson T.,Svanberg S.,Andersson-Engels S.,Johansson J.,Folestad S.. , “Time and wavelength resolved spectroscopy of turbid media using light continuum generated in a crystal fiber,” Opt. Express12(17), 4103–4112 (Year: 2004).10.1364/OPEX.12.00410319483952
26. Pifferi A.,Torricelli A.,Bassi A.,Taroni P.,Cubeddu R.,Wabnitz H.,Grosenick D.,Möller M.,Macdonald R.,Swartling J.,Svensson T.,Andersson-Engels S.,van Veen R. L.,Sterenborg H. J.,Tualle J. M.,Nghiem H. L.,Avrillier S.,Whelan M.,Stamm H.. , “Performance assessment of photon migration instruments: the MEDPHOT protocol,” Appl. Opt.44(11), 2104–2114 (Year: 2005).10.1364/AO.44.00210415838951
27. Selb J.,Joseph D. K.,Boas D. A.. , “Time-gated optical system for depth-resolved functional brain imaging,” J. Biomed. Opt.11(4), 044008 (Year: 2006).10.1117/1.233732016965165
28. Gagnon L.,Gauthier C.,Hoge R. D.,Lesage F.,Selb J.,Boas D. A.. , “Double-layer estimation of intra- and extracerebral hemoglobin concentration with a time-resolved system,” J. Biomed. Opt.13(5), 054019 (Year: 2008).10.1117/1.298252419021399
29. Ijichi S.,Kusaka T.,Isobe K.,Okubo K.,Kawada K.,Namba M.,Okada H.,Nishida T.,Imai T.,Itoh S.. , “Developmental changes of optical properties in neonates determined by near-infrared time-resolved spectroscopy,” Pediatr. Res.58(3), 568–573 (Year: 2005).10.1203/01.PDR.0000175638.98041.0E16148075
30. Zhao J.,Ding H. S.,Hou X. L.,Zhou C. L.,Chance B.. , “In vivo determination of the optical properties of infant brain using frequency-domain near-infrared spectroscopy,” J. Biomed. Opt.10(2), 024028 (Year: 2005).10.1117/1.189134515910101
31. Choi J.,Wolf M.,Toronov V.,Wolf U.,Polzonetti C.,Hueber D.,Safonova L. P.,Gupta R.,Michalos A.,Mantulin W.,Gratton E.. , “Noninvasive determination of the optical properties of adult brain: near-infrared spectroscopy approach,” J. Biomed. Opt.9(1), 221–229 (Year: 2004).10.1117/1.162824214715077
32. Comelli D.,Bassi A.,Pifferi A.,Taroni P.,Torricelli A.,Cubeddu R.,Martelli F.,Zaccanti G.. , “In vivo time-resolved reflectance spectroscopy of the human forehead,” Appl. Opt.46(10), 1717–1725 (Year: 2007).10.1364/AO.46.00171717356614
33. Ohmae E.,Ouchi Y.,Oda M.,Suzuki T.,Nobesawa S.,Kanno T.,Yoshikawa E.,Futatsubashi M.,Ueda Y.,Okada H.,Yamashita Y.. , “Cerebral hemodynamics evaluation by near-infrared time-resolved spectroscopy: correlation with simultaneous positron emission tomography measurements,” Neuroimage29(3), 697–705 (Year: 2006).10.1016/j.neuroimage.2005.08.00816165372
34. Kienle A.,Patterson M.,Dögnitz N.,Bays R.,Wagnières G.,van den Bergh H.. , “Noninvasive determination of the optical properties of two-layered turbid media,” Appl. Opt.37(4), 779–791 (Year: 1998).10.1364/AO.37.00077918268653
35. Martelli F.,Sassaroli A.,Del Bianco S.,Zaccanti G.. , “Solution of the time-dependent diffusion equation for a three-layer medium: application to study photon migration through a simplified adult head model,” Phys. Med. Biol.52(10), 2827–2843 (Year: 2007).10.1088/0031-9155/52/10/01317473354
36. Barnett A. H.,Culver J. P.,Sorensen A. G.,Dale A.,Boas D. A.. , “Robust inference of baseline optical properties of the human head with three-dimensional segmentation from magnetic resonance imaging,” Appl. Opt.42(16), 3095–3108 (Year: 2003).10.1364/AO.42.00309512790461
37. Boas D.,Culver J.,Stott J.,Dunn A.. , “Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head,” Opt. Express10(3), 159–170 (Year: 2002).19424345
38. Strangman G.,Franceschini M. A.,Boas D. A.. , “Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters,” Neuroimage18(4), 865–879 (Year: 2003).10.1016/S1053-8119(03)00021-112725763
39. S. Prahl, “Optical absorption of hemoglobin,” http://omlcogiedu/spectra/hemoglobin/summary.html (2002).
40. Yaroslavsky A.,Yaroslavsky I.,Goldbach T.,Schwarzmaier H.-J.. , “Optical properties of blood in the nearinfrared spectral range,” Proc. SPIE2678, 314–324 (Year: 1996).10.1117/12.239516
41. Fantini S.,Franceschini M.,Gratton E.. , “Semi-infinite-geometry boundary problem for light migration in highly scattering media: a frequency-domain study in the diffusion approximation,” J. Opt. Soc. Am. B11(10), 2128–2138 (Year: 1994).10.1364/JOSAB.11.002128
42. Fantini S.,Franceschini M. A.,Fishkin J. B.,Barbieri B.,Gratton E.. , “Quantitative determination of the absorption spectra of chromophores in strongly scattering media: a light-emitting-diode based technique,” Appl. Opt.33(22), 5204–5213 (Year: 1994).10.1364/AO.33.00520420935909
43. S. Meyer, Data Analysis for Scientists and Engineers (Wiley, 1975).
44. Heiskala J.,Neuvonen T.,Grant P. E.,Nissilä I.. , “Significance of tissue anisotropy in optical tomography of the infant brain,” Appl. Opt.46(10), 1633–1640 (Year: 2007).10.1364/AO.46.00163317356605
45. Fang Q.. , “Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates,” Biomed. Opt. Express1(1), 165–175 (Year: 2010).10.1364/BOE.1.00016521170299

Figures

[Figure ID: fig01]
Fig. 1. 

Segmented head models: 3D and 2D axial views. Segmentations consist of 4 tissue types: superficial layers (superf) of scalp and skull, CSF, gray matter (GM), and white matter (WM). Sources and detectors are represented by red crosses and blue dots respectively. The black bar (3 cm) in the bottom right of the image indicates the scale of all 2D slices.



[Figure ID: fig02]
Fig. 2. 

Retrieved absorption coefficients of a homogeneous slab (calibrated with μa = 0.1cm−1). Calibrated (dashed) and uncalibrated (solid) are compared to true (dash-dotted) values with respect to SDs for different homogeneous coefficients of absorption [cm−1]: 0.05 (blue); 0.10 (green); 0.15 (red); 0.20 (cyan); 0.25 (magenta); 0.30 (yellow).



[Figure ID: fig03]
Fig. 3. 

Estimated μa,brain when the head is treated as a homogeneous medium, i.e. when μa,superf = μa,brain = [0.05,0.06, . . . ,0.29,0.3]cm−1. For each sub-figure, solid lines represent the 7 groups of SDs used in the estimation (legend is shown in top left sub-figure) and dotted lines represent the 1-ratio between estimated and true brain absorption values.



[Figure ID: fig04]
Fig. 4. 

Uncertainty and retrieved absorption for three particular cases of the homogeneous head: μa,superf = μa,brain = [0.05,0.15,0.3]cm−1 for the 7 sets of SDs.



[Figure ID: fig05]
Fig. 5. 

Averages and standard deviations of the relative error on retrieved μa,brain over 26 combinations of head absorption [0.05,0.06, . . . ,0.29,0.3]cm−1, for the 7 sets of SDs.



[Figure ID: fig06]
Fig. 6. 

Brain absorption coefficients μa,brain estimated across subject age (row-wise) and for two values of superficial absorption μa,superf = [0.1,0.2]cm−1 (column-wise). For each sub-figure, solid lines represent the 7 groups of SDs used in the estimation (legend is shown in bottom left sub-figure), dotted lines represent the 1-ratio between estimated and true brain absorption values and dash-dotted lines represent the superficial absorption.



[Figure ID: fig07]
Fig. 7. 

Two-tissue head model: averages and standard deviations of the relative error in retrieved μa,brain over 52 combinations of μa,superf = [0.1,0.2]cm−1 and μa,brain = [0.05,0.06, . . . ,0.29,0.3]cm−1 for the 7 sets of SDs.



[Figure ID: fig08]
Fig. 8. 

Brain absorption coefficients μa,brain estimated across subject age (row-wise) and for two values of superficial absorption of μa = [0.1,0.2]cm−1 (column-wise) when the CSF tissue is taken into account. For each sub-figure, solid lines represent the 7 groups of SDs used in the estimation (legend is shown in bottom left sub-figure), dotted lines represent the 1-ratio between estimated and true brain absorption values and dash-dotted lines represent the superficial absorption. Note the presence of the axial slice of the newborn (NB) left and right hemispheres superimposed with the optical probe.



[Figure ID: fig09]
Fig. 9. 

Three-tissue head model: averages and standard deviations of the relative error on retrieved μa,brain over 52 combinations of μa,superf = [0.1,0.2]cm−1 and μa,brain = [0.05,0.06, . . . ,0.29,0.3]cm−1 for the 7 sets of SDs.



Tables
[TableWrap ID: tbl01] Table 1.  Optical properties of tissue types selected for the Monte Carlo simulations
Tissue n g μ′s [cm−1] μa [cm−1]
superf 1.3 0.01 5 [0.1,0.2]
1.3 0.01 10 [0.1,0.2]
brain 1.3 0.01 10 [0.05,0.06, . . . ,0.29,0.3]
CSF 1.3 0.01 0.1 0.04
1.3 0.01 10 [0.05,0.06, . . . ,0.29,0.3]


Article Categories:
  • Diffuse Optical Imaging

Keywords: OCIS codes: (170.3660) Light propagation in tissues, (110.3080) Infrared imaging, (170.5280) Photon migration.

Previous Document:  Simulating human photoreceptor optics using a liquid-filled photonic crystal fiber.
Next Document:  Diffuse optical cortical mapping using the boundary element method.