**) of the bottom layer in a four-layered human head model with time-domain near-infrared measurements. The difference curve in the temporal profiles of the light attenuation between an object and a reference medium, which are obtained from Monte Carlo simulations, is divided into segments along the time axis, and a slope for each segment is calculated to obtain the depth-dependent**

*μ*_{a}*(*

**μ**_{a}**) . The reduced scattering coefficient (**

*μ*^{ seg }_{a}**) of the reference is determined by curve fitting with the temporal point spread function derived from the analytical solution of the diffusion equation to the time-resolved reflectance of the object. The deviation of**

*μ*_{s}′**from the actual**

*μ*^{ seg }_{a}**is expressed by a function of the ratio of**

*μ*_{a}**in an earlier time segment to that in a later segment for mediums with different optical properties and thicknesses of the upper layers. Using this function, it is possible to determine the**

*μ*^{ seg }_{a}**of the bottom layer in a four-layered epoxy resin-based phantom. These results suggest that the method reported here has potential for determining the**

*μ*_{a}**of the cerebral tissue in humans.**

*μ*_{a}## 1.

## Introduction

In 1977, Jöbsis^{1} was the first to report near-infrared measurements of tissue oxygenation in living animals, and near-infrared spectroscopy (NIRS) has been developing into a useful tool for monitoring of cerebral hemodynamic changes. Today, NIRS is used in a wide variety of research fields and in clinical medicine, while quantitative and selective extraction of signals originating in the brain remains difficult. Recently, there have been reports that the skin blood flow influences NIRS signals, an issue that cannot be ignored in NIRS studies.^{2}^{,}^{3} Over the past 35 years, the NIRS community has focused on solving these issues and has developed a wide range of types of NIRS measurements. Among these, time-domain measurement [time-resolved spectroscopy (TRS)] is a very promising approach to these issues, as TRS provides a temporal profile of the detected light intensity [temporal point spread function (TPSF)], which carries information about depth-dependent attenuation based on the correlation of detection time and penetration depth of photons (time-domain depth sensitivity).^{4}5.6.^{–}^{7}

Propagation of photons in living tissue is described by the radiative transport equation (RTE); however, the RTE, an integro-differential equation, is not easy to solve even by numerical methods. Therefore, a diffusion approximation to the RTE is widely used instead, since light propagation can be considered isotropic at the macroscopic level, when the scattering interaction is much larger than the absorption and the point of interest is far from light sources and boundaries.^{8} In conventional time-domain measurements of the human head, the absorption coefficient (${\mu}_{a}$) of the cerebral tissue has commonly been determined by fitting the TPSF derived from an analytical solution of the photon diffusion equation (DE) to the profile of the time-resolved reflectance on the assumption that the head is semi-infinite and homogeneous.^{9}10.^{–}^{11} This approach selectively measures cerebral optical properties, but the selective and quantitative accuracies are inadequate because of limitations of the diffusion approximation to the RTE for inhomogeneous biological tissue.^{12}13.^{–}^{14} Other methods have also been proposed to estimate the optical properties of layered mediums, such as the multilayered (time-dependent) DE,^{15}16.^{–}^{17} also a method with spatially and time-resolved reflectance,^{18} and a method considering time-dependent mean partial path lengths.^{19} However, none of these methods have been applied to practical human head measurements.

Diffuse optical tomography (DOT) is a technique for reconstructing images of optical properties using multiple light sources and detectors^{20} and is thought to be the most promising technique for the quantitative detection of focal changes in cerebral hemodynamics. DOT can be performed with time-domain,^{21} frequency-domain,^{22} and continuous-wave (CW) measurement instruments.^{23} Recently, CW high-density DOT with high spatial resolution has been developed and applied to functional hemodynamic maps of the adult human visual cortex.^{24} However, this CW-DOT only provides qualitative images of the measured changes as it is based on linear image reconstruction. Nonlinear iterative reconstruction schemes are required to obtain quantitative images, and such procedures are still under development.^{25}

We have previously proposed a novel approach, “time segment analysis,” for determination of the ${\mu}_{a}$ of each layer in a two-layered turbid medium with a single-channel TRS instrument.^{26} In that approach, first a time-attenuation difference curve, derived by the subtraction of the time-resolved reflectance of a reference medium from that of an object, was divided into segments along the time axis (e.g., 500 ps time width). Then, a depth-dependent ${\mu}_{a}$ (time-segmented ${\mu}_{a}$) was estimated from the slope of each segment. The difference between the estimated ${\mu}_{a}$ and the actual ${\mu}_{a}$ in the lower layer then depends on the ratio of the actual ${\mu}_{a}$ in the upper layer to that in the lower layer, which was linearly related to the ratio of the estimated ${\mu}_{a}$ in an early time segment to that in a later segment. Using this relationship, it was possible to correct the ${\mu}_{a}$ obtained from the later time segment to derive the ${\mu}_{a}$ of the lower layer in two-layered mediums. An advantage of this approach is that it involves a single-distance measurement, which is more suitable than a multidistance measurement for measurements of the head with its heterogeneous superficial layer and curvature. A single distance measurement is also straightforward to expand into multichannel measurements. To apply our method to human head measurements, the validity and reliability of the method must be confirmed with three or more layered models and with other models with more complex structures. In the present study, the investigation focuses on the applicability of our method^{26} to a four-layered slab model with a Monte Carlo simulation and phantom experiments.

The article here first reports the results of simulations: (1) time-segmented ${\mu}_{a}$ under several conditions in which ${\mu}_{a}$ and the reduced scattering coefficient (${{\mu}_{s}}^{\prime}$) of each layer were changed, (2) the difference between the time-segmented ${\mu}_{a}$ at a later time and the actual ${\mu}_{a}$ of the bottom layer, and its relation to the ratio of the time-segmented ${\mu}_{a}$ in an earlier time segment to that of a later one, and (3) the mean partial path length of photons detected in each time segment to be able to interpret the time-segmented ${\mu}_{a}$ values. Then, the experimental results are reported: (4) the determination of ${\mu}_{a}$ of the bottom layer of a four-layered epoxy resin-based phantom with the “time segment analysis.”

## 2.

## Methods for Determining the ${\mu}_{a}$ of the Bottom Layer

The details of the method have been described in a previous paper,^{26} and the following is a brief summary.

## 2.1.

### Theoretical Basis

Two kinds of mediums, a reference medium and an object medium for the measurements, are considered. When a light impulse is irradiated on the surface of the mediums, the reflected light intensity at time $t$ for the reference and the object is expressed by ${I}_{R}(t)={S}_{R}(t)\mathrm{exp}(-{\mu}_{a,R}\xb7{c}_{R}\xb7t)$ and $I(t)=\phantom{\rule{0ex}{0ex}}S(t)\mathrm{exp}(-{\mu}_{a}\xb7c\xb7t)$,^{27}^{,}^{28} where $c$ is the light velocity in a medium and $S(t)$ is a scattering function which is dependent on the ${{\mu}_{s}}^{\prime}$ of the medium in the photon diffusion regime. The light attenuation, $A$, is defined as the logarithm of the inverse of the reflectance. According to the time-resolved Beer–Lambert law, the difference between the attenuation of an object $[A(t)]$ and of a reference $[{A}_{R}(t)]$, ${A}^{\mathrm{diff}}(t)$, can be expressed by

## (1)

$${A}^{\mathrm{diff}}(t)=A(t)-{A}_{R}(t)=\mathrm{ln}[{I}_{R}(t)/I(t)]=({\mu}_{a}\xb7c-{\mu}_{a,R}\xb7{c}_{R})t+{S}^{\mathrm{diff}}(t),$$In the case of an object that is a layered medium in which ${\mu}_{a}$ varies with layer (i.e., depth) and the reference is homogeneous, ${A}^{\mathrm{diff}}(t)$ is not a linear function of $t$ even though ${S}^{\mathrm{diff}}(t)$ can be ignored and $c={c}_{R}$, that is, ${\mu}_{a}$ is different at each detection time. The method introduces a time-dependent apparent absorption coefficient, ${\mu}_{a}(t)$, given by $\sum _{i}{\mu}_{ai}{l}_{i}(t)/L(t)$, where the subscript $i$ refers to the $i$’th layer in the object medium, ${l}_{i}(t)$ is the stochastic time-dependent partial path length in the $i$’th layer for the photons detected at time $t$, and $L(t)$ is the total path length of photons detected at time $t$, which is defined as $L(t)=ct$. Then $\mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t$ is expressed as follows:

## (2)

$$\mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t=(\mathrm{d}/\mathrm{d}t)[{\mu}_{a}(t)ct-{\mu}_{a,R}\xb7{c}_{R}\xb7t]+\mathrm{d}{S}^{\mathrm{diff}}(t)/\mathrm{d}t\phantom{\rule{0ex}{0ex}}=\sum _{i}[\mu ai\xb7\mathrm{d}li(t)/\mathrm{d}t]-{\mu}_{a,R}\xb7{c}_{R}+\mathrm{d}{S}^{\mathrm{diff}}(t)/\mathrm{d}t,\phantom{\rule{0ex}{0ex}}$$## 2.2.

### Time Segment Analysis

To simplify the process of the analysis, the ${A}^{\mathrm{diff}}-t$ curve within a short time range (e.g., a few hundred picoseconds, $\mathrm{\Delta}t$) is considered a straight line, and a term for the mean absorption differences, $\langle \mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t\rangle $, is introduced, the value of this is obtained by $\mathrm{\Delta}{A}^{\mathrm{diff}}(t)/\mathrm{\Delta}t$, rather than $\mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t$. In this study, the ${A}^{\mathrm{diff}}-t$ curve in the range 500 to 3500 ps is divided into six time segments of 500 ps, and the slope of each time segment is estimated ($\langle \mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t\rangle $) by linear regression. The mean ${\mu}_{a}(t)$ of each time segment is referred to as the time-segmented ${\mu}_{a}$, ${\mu}_{a}^{\mathrm{seg}}$. Then, assuming that a reference and an object have the same ${{\mu}_{s}}^{\prime}$ and refractive index, this gives the following equation,

## (3)

$$\langle \mathrm{d}{A}^{\text{diff}}(t)/\mathrm{d}t\rangle =({\mu}_{a}^{\mathrm{seg}}-{\mu}_{a,R})c.$$The time-segmented ${\mu}_{a}({\mu}_{a}^{\mathrm{seg}})$ was calculated by substituting $\langle \mathrm{d}{A}^{\mathrm{diff}}(t)/\mathrm{d}t\rangle $, ${\mu}_{a},R$, and $c$ in Eq. (3).

The segments of the ${A}^{\mathrm{diff}}-t$ curve were numbered as $K=\mathrm{I},\mathrm{II},\dots ,\mathrm{VI}$ in time order, so the segment of a time range from 500 to 1000 ps is referred to as segment-I and that from 2500 to 3000 ps was as segment-V, ${\mu}_{a}^{\mathrm{seg}}$ ($K$) denotes the ${\mu}_{a}^{\mathrm{seg}}$ in segment $K$.

## 3.

## Simulation and Experiment

## 3.1.

### Monte Carlo Simulations – Calculation of the Time-Resolved Reflectances

To calculate the reflectances from four-layered (object) and from homogeneous (reference) mediums, we used the Monte Carlo code developed by Wang and Jacques.^{29} The code was modified to fit our measurement system. Source photons were perpendicularly irradiated onto the surface of the semi-infinite mediums, and photons that were emitted from a detector position on the surface were all detected with a time step of 10 ps, and the source–detector distance was 30 mm. Reflectance profiles simulated in an anisotropic setting (scattering coefficient ${\mu}_{s}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$, anisotropy parameter $g=0.9$), were compared with those in an isotropic setting (${\mu}_{s}=1.0\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$) and it was confirmed that the differences between the two were small enough to be disregarded in these simulating parameter conditions. Thus, an isotropic scattering was assumed, reducing the calculation time, and the calculations were repeated until the number of detected photons reached 1,000,000. A light velocity of $0.219\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}/\mathrm{ps}$ (refractive index 1.37) was used in all the calculations.

In the previous study,^{26} we confirmed that temporal profiles of detected light intensity in homogeneous and two-layered phantoms, which consisted of two layers with optical properties the same as those of the homogeneous one, were almost in agreement. Thus, the effect of the interface between two layers on the pattern of photon propagation was considered to be negligible in the Monte Carlo simulation and phantom measurement under present study conditions.

The statistical error of the Monte Carlo simulation was estimated by fitting a theoretical TPSF for a semi-infinite homogenous medium (based on the solution of the DE derived by employing the extrapolated boundary condition of the reflectance mode^{30}31.^{–}^{32}) to the TPSF obtained from the Monte Carlo simulation (time range of 0 to 4500 ps). The photon diffusion coefficient was expressed by $D=1/3{{\mu}_{s}}^{\prime}$.^{33} In the fitting procedure, a weighted least squares fitting method based on the Levenberg–Marquardt Method,^{34}^{,}^{35} which is an iterative improved method, was used. In case of practical TRS measurements of inhomogeneous objects, such as the human head, when the chi square nu (${\chi}_{\nu}^{2}$) value is 0.7 to 1.3, it is considered that the theoretical TPSF is well fitted to the observed one.^{36}^{,}^{37} In the present study, it was assumed that the best fitting TPSF was the reference true solution, and that the ${\chi}_{\nu}^{2}$ value mainly represented statistical bias and statistical variations that were attributed to the Monte Carlo simulation. The ${\chi}_{\nu}^{2}$ values in all the mathematical models described in Sec. 3.2 (Table 1) were $1.09\pm 0.05$ ($\text{mean}\pm \mathrm{SD}$). Thus, it was concluded that the statistical error of the Monte Carlo simulation was negligible.

## Table 1

Optical properties and layer thickness of four layered models used in simulation.

Layer (#) | μs′ (mm−1) | μa (mm−1) | d (mm) |
---|---|---|---|

Group A | |||

1 | 1.3 | 0.0013–0.02 | (d1,d2)=(3,5), (3, 7), |

2 | 1.0 | 0.011–0.014 | (4, 7), (5, 7), (5, 5), (5, 9) |

3 | 0.1 | 0.0033 | 1 |

4 | 1.7 | 0.01–0.02 | 90 |

Group B | |||

1 | 1.1 | 0.0013–0.02 | (d1,d2)=(3,7), (5, 7) |

2 | 0.8 | 0.012 | |

3 | 0.1 | 0.0033 | 1 |

4 | 1.7 | 0.01–0.02 | 90 |

Group C | |||

1 | 1.8 | 0.0013–0.02 | (d1,d2)=(5,7) |

2 | 1.2 | 0.012 | |

3 | 0.1 | 0.0033 | 1 |

4 | 1.7 | 0.01–0.02 | 90 |

## 3.2.

### Four-Layered Models for Simulation

We used a four-layered semi-infinite model as a human head model. The layers are numbered $i=1,2,3$, and 4 in order of depth from the outside, and denoted by #$i$: the layers #1 to #4 correspond to the scalp, skull, cerebrospinal fluid (CSF), and brain, respectively. The ${{\mu}_{s}}^{\prime}$, ${\mu}_{a}$, and the thickness $d$ of each layer are denoted by ${{\mu}_{si}}^{\prime}$, ${\mu}_{ai}$, and ${d}_{i}$, respectively. Table 1 shows the values used in this paper. The optical properties of each layer were chosen from the published data.^{7}^{,}^{38}^{,}^{39} The depth from the scalp to the cortical surface was set to 9.0 to 15.0 mm in the forehead region.^{40}41.^{–}^{42}

The models were classified into three groups according to the ${{\mu}_{s1}}^{\prime}$ and ${{\mu}_{s2}}^{\prime}$ values: 1.3 and $1.0\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ for Group A, 1.1 and $0.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ for Group B, and 1.8 and $1.2\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ for Group C. The ${{\mu}_{s3}}^{\prime}$ and ${{\mu}_{s4}}^{\prime}$ values were set to 0.1 and $1.7\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ in all three groups. In Group A, ${\mu}_{a}$ was varied in the range of 0.013 to $0.02\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ in layer #1, 0.011 to $0.014\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ in layer #2, and 0.01 to $0.02\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ in layer #4 separately. In Groups B and C, only ${\mu}_{a1}$ and ${\mu}_{a4}$ were varied and ${\mu}_{a2}$ was set to $0.012\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$. In all the models, ${\mu}_{a1}$ was set to be larger than ${\mu}_{a2}$ as the ${\mu}_{a}$ of the scalp is generally larger than that of the skull in adult heads,^{43}44.^{–}^{45} and ${\mu}_{a3}$ was $0.0033\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$. The thicknesses of layers #1 and #2, (${d}_{1},{d}_{2}$) (mm), were set to (3, 5), (3, 7), (4, 7), (5, 7), (5, 5), or (5, 9) in Group A, (3, 7) or (5, 7) in Group B, and (5, 7) in Group C, whereas ${d}_{3}$ and ${d}_{4}$ were set as 1 and 90 mm, respectively. We mainly evaluated temporal variations in ${\mu}_{a}^{\mathrm{seg}}$ and differences between ${\mu}_{a}^{\mathrm{seg}}$ and the actual ${\mu}_{a}$ for the Group A models. The models in Groups B and C were used to investigate the effects of the ${{\mu}_{s}}^{\prime}$ variation in the upper layers on ${\mu}_{a}^{\mathrm{seg}}$. The ${{\mu}_{s}}^{\prime}$ of the homogeneous reference was changed within the range of 1.1 to $1.8\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ as described below, while the ${\mu}_{a}$ and the thickness were set to $0\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ and 100 mm, respectively.

## 3.3.

### Determination of ${{\mu}_{s}}^{\prime}$ for the Reference Mediums

In the previous study,^{26} it was found that the results of our method depend on the ${{\mu}_{s}}^{\prime}$ value of the reference, which was ideally the same as that of the upper layer of a two-layered object medium. In general, the ${{\mu}_{s}}^{\prime}$ value of a multilayered medium determined by the curve fitting procedure as described in Sec. 3.1 is closer to the ${{\mu}_{s}}^{\prime}$ of the upper layer than that of the lower layer. This can be explained by the fact that the slope of the rising phase of the TPSF is dependent on the ${{\mu}_{s}}^{\prime}$ value but independent of the ${\mu}_{a}$ value, when the optical properties are in the range of the published data for the biological tissue.^{46} In this article, the ${{\mu}_{s}}^{\prime}$ value of the reference was determined by curve fitting with the TPSF derived from an analytical solution of the DE for a semi-infinite homogenous medium to the profile of the time-resolved reflectance of the object medium in the time range of 250 to 1000 ps. The reason why we chose this time range is described in Discussion (Sec. 5.2).

## 3.4.

### Calculation of Mean Partial Path Lengths

As shown in Eqs. (2) and (3), ${\mu}_{a}^{\mathrm{seg}}$ depends on the temporal change in ${l}_{i}(t)$ in the time range of each segment in each layer. Therefore, to interpret the calculated ${\mu}_{a}^{\mathrm{seg}}$, the mean partial path lengths in each layer for a four-layered medium were also estimated by the Monte Carlo simulation. Two models in Group A [${\mu}_{a1}=0.015$, ${\mu}_{a2}=0.012$, ${\mu}_{a4}=0.02$, $({d}_{1},{d}_{2})=(5,5)$; and ${\mu}_{a1}=0.015$, ${\mu}_{a2}=0.012$, ${\mu}_{a4}=0.01$, $({d}_{1},{d}_{2})=(5,5)$], were perpendicularly partitioned into 1 mm thick laminae and the attenuation changes ($\mathrm{\Delta}A$) caused by altering ${\mu}_{a}$ in a lamina by $+10\%$ of the initial value were calculated. The time-independent mean partial path length of photons in the lamina where ${\mu}_{a}$ is changed can be estimated by the $\mathrm{\Delta}A/\mathrm{\Delta}{\mu}_{a}$ of the whole time range of the time-resolved reflectance. The present study estimated a time-dependent mean partial path length given by the time segment unit [denoted as the time-segment dependent mean partial path length (s-mPPL)] at the depth of $z$ to $z+1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ ($z=0,1,2,\dots ,19$) obtained by dividing the $\mathrm{\Delta}A$ calculated in each time segment by $\mathrm{\Delta}{\mu}_{a}$. We also calculated a time-segment dependent mean total path length (s-mTPL) in the medium, the product of the light velocity and the mean time of flight of photons in each time segment.

## 3.5.

### Epoxy Resin-Based Phantoms

Four-layered (object) and homogeneous (reference) phantoms, $140\times 140\times 63$ (50 for the reference) mm, were prepared. The base of the phantoms was epoxy resin (Beuhler, Lake Bluff, Illinois). We adjusted the ${{F}_{s}}^{\prime}$ and ${\mu}_{a}$ to 0.1 to 1.77 and 0.0033 to $0.015\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ at 760 nm, respectively, by adding titanium oxide (Tian Kogyo, Ube, Japan) as scatterers and ink (Greenish brown: Chugai Pharmaceutical, Tokyo, Japan) as an absorber. The thickness of each layer and the optical properties of the object are shown in Table 2. These optical properties are theoretical values based on the concentrations of titanium oxide and ink. Unlike the mathematical phantoms, there were no differences in the ${{\mu}_{s}}^{\prime}$ values of the objects determined in the two time ranges used in the fitting ($1.520\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$, 250 to 1000 ps; and $1.502\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$, 0 to 4500 ps). The ${\mu}_{a}$ and ${{\mu}_{s}}^{\prime}$ of the reference were 0.001 and $1.52\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$, respectively (Table 2).

## Table 2

Optical properties at 760 nm and layer thickness of epoxy resin-based phantoms, where μs′ and Fa are theoretical values based on the concentration of titanium oxide and ink.

Layer (#) | μs′ (mm−1) | μa (mm−1) | d (mm) |
---|---|---|---|

Object | |||

1 | 1.36 | 0.016 | 4 |

2 | 1.0 | 0.012 | 6 |

3 | 0.1 | 0.0033 | 3 |

4 | 1.77 | 0.015 | 50 |

Reference | |||

1 | 1.52 | 0.001 | 50 |

## 3.6.

### Instrumentation

A two-channel TRS instrument (TRS-20, Hamamatsu Photonics, Hamamatsu, Japan) was employed. The TRS-20 consists of three pulsed laser diodes with different wavelengths (760, 800, and 830 nm) with durations of around 150 ps [full width at half maximum (FWHM)] at the repetition rate of 5 MHz, a high-speed photomultiplier tube for single photon detection, and a circuit for time-resolved measurements based on a time correlated single photon counting method. Incident light is delivered to the sample through an optical light guide [GI type, 200 μm core diameter, numerical aperture $(\mathrm{NA})=\phantom{\rule{0ex}{0ex}}0.25$]. Reflected light is collected by a fiber bundle (3-mm diameter, $\mathrm{NA}=0.26$). Minimum data acquisition time is 100 ms.

## 3.7.

### Time-Resolved Measurements of a Phantom

The incident and detecting light guides were placed on the upper surface of a phantom, separated by 30 mm. Measurements with an accumulated time of 10 s were conducted at five different positions within an area 40 mm from the edge of the phantom to avoid the distortion of photon propagation due to the edge. Specular reflection and light leakage were prevented with a black light-guide holder. Instrument responses of TRS-20 were measured by placing the incident fiber opposite the detecting fiber with a neutral density filter between the two. The instrument response of was around 250 ps FWHM at all wavelengths.

The optical properties of the phantoms were determined by curve fitting with the TPSF derived from an analytical solution of the DE with the extrapolated boundary condition, convoluted by the instrument function to the measured profile of the time-resolved reflectance as described in Sec. 3.1. The light velocity in the phantoms was assumed to be $0.19\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}/\mathrm{ps}$, corresponding to a refractive index of 1.58. Before the calculation of ${A}^{\mathrm{diff}}(t)$, the profiles of the time-resolved reflectance measured by TRS were deconvoluted by the instrument function by using Bayesian deconvolution.^{47}

## 4.

## Results

## 4.1.

### Temporal Variations in ${\mu}_{a}^{\mathrm{seg}}$

Figure 1 shows the temporal variations in ${\mu}_{a}^{\mathrm{seg}}$ obtained from the nine models in Group A. Figure 1(a) shows results in the mediums where both ${\mu}_{a1}$ and ${\mu}_{a2}$ were larger than ${\mu}_{a4}$ (squares) and where ${\mu}_{a1}$ was larger but ${\mu}_{a2}$ was equal to ${\mu}_{a4}$ (triangles). Figure 1(b) shows results in the mediums where ${\mu}_{a2}$ was smaller than ${\mathrm{FFFF}}_{a4}$ and ${\mu}_{a1}$, while ${\mu}_{a1}$ were larger than (diamonds), equal to (crosses), or smaller than ${\mu}_{a4}$ (circles). The ${\mu}_{a}^{\mathrm{seg}}$ values are expressed as the ratio to ${\mu}_{a4}$. When both ${\mu}_{a1}$ and ${\mu}_{a2}$ were larger than ${\mu}_{a4}$ [squares in Fig. 1(a)], the ratios of the ${\mu}_{a}^{\mathrm{seg}}$ to the ${\mu}_{a4}$ $[{\mu}_{a}^{\mathrm{seg}}(k)/{\mu}_{a4}]$ term were above 1 in all the time segments. The values gradually decreased to become close to 1 as the time segments become later, and reached a value slightly larger than ${\mu}_{a4}$. When ${\mu}_{a2}$ was equal to ${\mu}_{a4}$ [triangles in Fig. 1(a)], ${\mu}_{a}^{\mathrm{seg}}(k)/{\mu}_{a4}$ is nearest 1 in the time segments later than IV. For all the models in Fig. 1(a), the upper layer thickness (${d}_{1},{d}_{2}$) influenced the values of ${\mu}_{a}^{\mathrm{seg}}$ in segments-I and -II but not those in the later segments.

In the models where ${\mu}_{a2}$ was smaller than ${\mu}_{a4}$ [Fig. 1(b)], the ${\mu}_{a}^{\mathrm{seg}}(k)/{\mu}_{a4}$ in all the time segments later than segment-II were smaller than 1. When ${\mu}_{a1}$ was equal to or larger than ${\mu}_{a4}$, the ${\mu}_{a}^{\mathrm{seg}}$ was almost constant in segments-II to -VI, and a little smaller than the ${\mu}_{a4}$ values. The values of ${\mu}_{a}^{\mathrm{seg}}$ in the later time segments were independent of the value of ${\mu}_{a1}$. When both ${\mu}_{a1}$ and ${\mu}_{a2}$ were smaller than ${\mu}_{a4}$, the ${\mu}_{a}^{\mathrm{seg}}(k)/{\mu}_{a4}$ values were much smaller than one in all the time segments. The ${\mu}_{a}^{\mathrm{seg}}$ decreased between segments I and II, then ${\mu}_{a}^{\mathrm{seg}}$ increased slightly in the later time segments, and remained about 20% smaller than ${\mu}_{a4}$. For all the models in Fig. 1(b), the differences between ${\mu}_{a4}$ and ${\mu}_{a}^{\mathrm{seg}}$ in the later time segments were larger in the medium with $({d}_{1},{d}_{2})=(5,9)$ than in that with (3, 7). Further, the ${\mu}_{a}^{\mathrm{seg}}$ values of the medium with (5, 9) were smaller than those of the medium with (3, 7) in all segments. The same tendency of the temporal variations of ${\mu}_{a}^{\mathrm{seg}}$ was observed in other models with different sets of ${\mu}_{ai}$ and thicknesses.

## 4.2.

### Differences Between the ${\mu}_{a}^{\mathrm{seg}}$ and the Actual ${\mu}_{a}$ in the Bottom Layer

In the previous study of two-layered models, it was found that the ratio of ${\mu}_{a}^{\mathrm{seg}}$ (V) to the actual ${\mu}_{a}$ of the bottom layer was expressed by a function of the ratio of ${\mu}_{a}^{\mathrm{seg}}$ in an earlier time segment to that in a later one. In this article we examined the relationship between the ratio of ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ to ${\mu}_{a4}[{\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}]$ and the ratio of ${\mu}_{a}^{\mathrm{seg}}$ in segment-I to that in segment-V (${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$) obtained from Group A (Fig. 2). All data obtained from the various mediums with different sets of (${d}_{1},{d}_{2}$) were plotted together in the graph. The plot shows that the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ could be expressed as one polynomial function of the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ $[{\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}=-1.70{({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})}^{2}+\phantom{\rule{0ex}{0ex}}4.82({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})-2.35]$ with an adjusted coefficient of determination ${R}^{2}=0.976$. The deviation of ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ from the regression curve was less than 6% (mean error was 1.6%).

Figure 3 is a plot of the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ versus the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ for Groups B and C, where the ${{\mu}_{s}}^{\prime}$ values of layers #1 and #2 are different from those in Group A. The solid line represents the regression curve of the relation between ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ and ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ for the Group A results, as also shown in Fig. 2. The data points from Groups B and C are in good agreement with the regression curve obtained from the data sets of Group A (an adjusted ${R}^{2}=0.975$). The deviation of ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ in Groups B and C from the regression curve was less than 8% (mean error was 2.7%). This result suggests the conclusion that the relation between ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ and the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ can be expressed by $[-1.70{({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})}^{2}+\phantom{\rule{0ex}{0ex}}4.82({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})-2.35]$ irrespective of the ${{\mu}_{s}}^{\prime}$ values of the upper layers.

## 4.3.

### Temporal Variations in the s-mPPL in Each Layer

Figures 4(a) and 4(b) show temporal variations in the s-mPPL in the four layers. In the medium where ${\mu}_{a4}$ was larger than ${\mu}_{a2}$ [Fig. 4(a)], the s-mPPLs in layers #2 and #4 markedly increased with time, and the s-mPPL in layer #2 was the longest in the latter time segments. The temporal changes in s-mPPL in layer #1, which was longer than that in layer # 4, was small compared to those in layers #2 and #4. The s-mPPL in layer #3, which was the shortest, also showed small increments with time.

In the medium where ${\mu}_{a4}$ was smaller than ${\mu}_{a2}$ [Fig. 4(b)], the s-mPPL in layer #4 markedly increased with time, with a steeper slope than that in the medium in Fig. 4(a). The s-mPPLs in layer #2 increased with time in time segments I–IV, with amplitudes very similar to those in Fig. 4(a). In time segment-V, however, the s-mPPL in layer #2 did not increase. The temporal variations in the s-mPPLs in layers #1 and #3 were relatively small, especially after time segment-III.

## 4.4.

### Phantom Experiment

Figures 5(a) and 5(b) show the profiles of the time-resolved reflectances of the object [$\mathrm{I}(t)$] and reference [${\mathrm{I}}_{R}(t)$] at 760 nm, which were deconvoluted by the incident light, and the ${A}^{\mathrm{diff}}-t$ curve, respectively. The profiles of ${A}^{\mathrm{diff}}(t)$ were reliable in the time range of 0.5 to 3 ns (segments-I to -V). From Fig. 5(b), ${\mu}_{a}^{\mathrm{seg}}(\mathrm{I})$ and ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ were estimated at each measurement position. Substituting these values into the function $[{\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}=-1.70{({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})}^{2}+4.82({\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}})-2.35]$, ${\mu}_{a4}$ was estimated to be $0.0156\pm 0.002/\mathrm{mm}$ ($\text{mean}\pm \mathrm{SD}$). This value was very close to the actual ${\mu}_{a4}$ of the object ($0.015/\mathrm{mm}$).

## 5.

## Discussion

## 5.1.

### Characterization of ${\mu}_{a}^{\mathrm{seg}}$ in View of the Change in s-mPPL in Each Layer

The temporal variations in ${\mu}_{a}^{\mathrm{seg}}$ in four-layered mediums (Fig. 1) largely agreed with the previous results for a two-layered model.^{26} The ${\mu}_{a}^{\mathrm{seg}}$ changed with the time segment, reflecting the differences in ${\mu}_{a}$ in the depth direction. The main finding is that the ${\mu}_{a}^{\mathrm{seg}}$ values in later time segments exclusively depend on the differences between ${\mu}_{a2}$ and ${\mu}_{a4}$. That is, the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ accurately represents the ${\mu}_{a4}$ in the mediums where ${\mu}_{a2}$ was larger than or equal to ${\mu}_{a4}$ [Fig. 1(a)], whereas in mediums where ${\mu}_{a2}$ was smaller than ${\mu}_{a4}$, the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ was smaller than ${\mu}_{a4}$ even in the later time segments irrespective of ${\mu}_{a1}$ values [Fig. 1(b)]. In addition, the dependence of ${\mu}_{a}^{\mathrm{seg}}$ on the upper-layer thickness varied with the difference between ${\mu}_{a2}$ and ${\mu}_{a4}$ (black versus grey symbols in Fig. 1).

As defined in Sec. 2.2, ${\mu}_{a}^{\mathrm{seg}}$ represents the mean ${\mu}_{a}(t)$ in a segment. The ${\mu}_{a}^{\mathrm{seg}}$ value depends on the mean temporal change in ${l}_{i}(t)$ for each time segment in the individual layers [Eq. (2)]. Therefore, a consideration of the contribution of the mean change in ${l}_{i}(t)$ to the change in $L(t)$ ($=c\mathrm{\Delta}t$; it is constant in all the segments here) in each time segment is helpful to interpret the findings mentioned above. From Figs. 4(a) and 4(b), we can predict the mean temporal change in ${l}_{i}(t)$ and its amplitude within the segment (500 ps) by interpolation of s-mPPL between the time segments. Here, it should be noted that the temporal patterns of changes in the s-mPPLs provide more important information than the amplitude of the changes as described below.

When ${\mu}_{a2}$ was smaller than ${\mu}_{a4}$, s-mPPL in layer #2 increased considerably in all the time ranges [Fig. 4(a)]. This implies that the contribution of the mean change in ${l}_{2}(t)$ to the change in $L(t)$ was high enough to strongly influence the ${\mu}_{a}^{\mathrm{seg}}$ even in later time segments. This dominant change in ${l}_{2}(t)$ can explain the finding that the ${\mu}_{a}^{\mathrm{seg}}$ at the later time represented values smaller than ${\mu}_{a4}$ in the mediums where ${\mu}_{a2}$ is smaller than ${\mu}_{a4}$. The dominance of the change in ${l}_{2}(t)$ also accounts for the difference in ${\mu}_{a}^{\mathrm{seg}}$ at later time segments between the mediums with $({d}_{1},{d}_{2})=(3,7)$ and with (5, 9). The reason why there was little effect of ${\mu}_{a1}$ on the temporal pattern of changes in ${\mu}_{a}^{\mathrm{seg}}$ may be ascribed to the smaller temporal change in s-mPPL in layer #1.

When ${\mu}_{a2}$ was larger than ${\mu}_{a4}$, ${\mu}_{a}^{\mathrm{seg}}$ at segments-V and -VI mainly represented the ${\mu}_{a4}$ [Fig. 1(a)]. This is explained by the finding that there was no difference in s-mPPL in layer #2 between segment-IV and segment-V, while s-mPPL in layer #4 in segment-V was longer than that in segment-IV [Fig. 4(a)], which indicates that the contribution of the mean change in ${l}_{2}(t)$ to the change in $L(t)$ in these segments was smaller than that of ${l}_{4}(t)$. Therefore, it was concluded that the size difference between the ${\mu}_{a2}$ and the ${\mu}_{a4}$ is critical for the change in ${l}_{i}(t)$ in a later time segment, which determines ${\mu}_{a}^{\mathrm{seg}}$ values.

## 5.2.

### Determination of ${\mu}_{a}$ in the Bottom Layer by Correction of ${\mu}_{a}^{\mathrm{seg}}$

Like in two-layered models, when the ${\mu}_{a}$ value of the bottom layer was larger than that of the upper layer, the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ was smaller than the ${\mu}_{a4}$. With two-layered models, the ratio of ${\mu}_{a}^{\mathrm{seg}}$(IV) or (V) to the real ${\mu}_{a}$ of the lower layer is expressed by the affine function of the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$, which enables estimation of the ${\mu}_{a}$ of the lower layer.^{26} The difference between ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})$ and ${\mu}_{a4}$ in the four-layered models also depends on the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$, though the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ was not expressed by an affine function but by a polynomial function of the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ (Figs. 2 and 3). This relationship between the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ and ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ was independent of the ${{\mu}_{s}}^{\prime}$ values and thicknesses of the layers #1 and #2. In addition to the high ${R}^{2}$ values, experimental measurements with the epoxy resin-based phantom support the validity of using this function, although there is no functional explanation why the relationship between ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ and ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ in various object mediums with different conditions of the optical properties and thicknesses of the upper layers can be expressed by a single polynomial function.

## 5.3.

### Determination of an Adequate Reference in Time Segment Analysis

We have reported that the difference in ${{\mu}_{s}}^{\prime}$ between the upper layer of the object and the reference significantly influenced the ${\mu}_{a}^{\mathrm{seg}}$ in the time segments earlier than 2500 ps.^{26} However, this difference did not influence the ${\mu}_{a}^{\mathrm{seg}}$ after 2500 ps, and the difference in ${{\mu}_{s}}^{\prime}$ between the lower layer of the object and the reference rarely affected the ${\mu}_{a}^{\mathrm{seg}}$ values in all the time segments. The present study used a curve fitting procedure to determine the ${{\mu}_{s}}^{\prime}$ value of the reference. The effect of the time range of the curve fitting on the ${\mu}_{a}^{\mathrm{seg}}$ within the range of 250 to 3000 ps was investigated. When the ${{\mu}_{s}}^{\prime}$ of the reference was determined by fitting over the wider time ranges, the plots of ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ versus ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ did not converge on a single curve. Therefore, we used the time range of 250 to 1000 ps as the most appropriate range for fitting to determine the averaged ${{\mu}_{s}}^{\prime}$ of the upper layers. However, there were no differences in ${{\mu}_{s}}^{\prime}$ values of the epoxy resin-based phantom of the two fitting ranges (250 to 1000 ps and 0 to 4500 ps). In practical measurements, an inappropriate convolution process by the incident light pulse and the time lag of the pulse launching time may induce errors when deriving the ${{\mu}_{s}}^{\prime}$, whereas our simulation did not include noise. This may account for the discrepancy between the simulations and the experimental results.

## 5.4.

### Possibility of Application of the Present Method to Human Head Measurements

As described above, practical measurements are, in general, accompanied by noises arising from various causes, such as a change in coupling between a light guide and the object’s surface. In this article, we examined the applicability of the present method to practical measurements by using the four layered slab phantom. The difference between the estimated ($0.0156\pm 0.002/\mathrm{mm}$) and actual ${\mu}_{a}$ values ($0.015/\mathrm{mm}$) was maximally about $0.008/\mathrm{mm}$ (about 5%). Considering the deviation of ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ from the regression curve (Figs. 2 and 3) that was less than 8%, the estimation error is acceptable. Further, since the relationship between the ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ and ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ was independent of the ${{\mu}_{s}}^{\prime}$ values and thicknesses of the layers #1 and #2 in the range used in this study (${{\mu}_{s}}^{\prime}$ of 0.8 to 1.8/mm, thickness of less than 15 mm), it is expected that the applicability of the present approach could be extended to *in vivo* measurements.

In the adult human head, the total thickness of scalp, skull, and CSF is 11 to 15 mm in the forehead, and about 12 mm in the temporal and occipital areas.^{32}^{,}^{33}^{,}^{42} Therefore, it is likely that the ${\mu}_{a}$ of the brain can be estimated by the present method on the assumption that the human head is a four-layered slab model. It is however difficult to simply apply this method to measurements in the parietal region, where the total thickness is 15 to 25 mm. Further, the four-layered slab would be too simple for a human head model. As the next step, thus, we will investigate the applicability of our “time segment analysis” to more realistic human head models with the total thickness of the extracerebral tissues in the range of 8 to 25 mm.

## 6.

## Conclusion

This study demonstrated that the “time segment analysis” of time resolved reflectance could be applied to four-layered slab models. The ${\mu}_{a}^{\mathrm{seg}}(\mathrm{V})/{\mu}_{a4}$ was expressed by a unique function of the ${\mu}_{a}^{\mathrm{seg}}\text{\hspace{0.17em}}{\text{ratio}}_{\mathrm{I}/\mathrm{V}}$ on the condition that the thickness of the upper layer was less than 15 mm. Since the ${\mu}_{a}$ of the bottom layer can be determined by this function, the approach reported here has the potential to selectively and quantitatively measure hemoglobin concentrations in the human brain.

## Acknowledgments

The authors very much appreciate the helpful discussion with Y. Yamashita and M. Oda in Hamamatsu Photonics KK. This research was supported by Grants-in Aid for Scientific Research (C) from the Japan Society for the Promotion of Science (No. 16500468).

## References

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

T. Takahashiet al., “Influence of skin blood flow on near-infrared spectroscopy signals measured on the forehead during a verbal fluency task,” NeuroImage 57(3), 991–1002 (2011).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2011.05.012Google Scholar

E. Kirilinaet al., “The physiological origin of task-evoked systemic artifacts in functional near infrared spectroscopy,” NeuroImage 61(1), 70–81 (2012).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2012.02.074Google Scholar

S. Del BiancoF. MartelliG. Zaccanti, “Penetration depth of light re-emitted by a diffusive medium: theoretical and experimental investigation,” Phys. Med. Biol. 47(23), 4131–4144 (2002).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/47/23/301Google Scholar

J. SelbD. K. JosephD. A. Boas, “Time-gated optical system for depth-resolved functional brain imaging,” J. Biomed. Opt. 11(4), 044008 (2006).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2337320Google Scholar

A. Torricelliet al., “Time-resolved reflectance at null source-detector separation: improving contrast and resolution in diffuse optical imaging,” Phys. Rev. Lett. 95(7), 078101 (2005).PRLTAO0031-9007http://dx.doi.org/10.1103/PhysRevLett.95.078101Google Scholar

G. Zaccantiet al., “Optical properties of biological tissues,” Proc. SPIE 2389, 513–521 (1995).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.210000Google Scholar

A. Ishimaru, “Diffusion of a pulse in densely distributed scatters,” J. Opt. Soc. Am. 68(8), 1045–1050 (1978).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSA.68.001045Google Scholar

S. Ijichiet al., “Quantification of cerebral hemoglobin as a function of oxygenation using near-infrared time-resolved spectroscopy in a piglet model of hypoxia,” J. Biomed. Opt. 10(2), 024026 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.1899184Google Scholar

V. Quaresimaet al., “Bilateral prefrontal cortex oxygenation responses to a verbal fluency task: a multichannel time-resolved near-infrared topography study,” J. Biomed. Opt. 10(1), 011012 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.1851512Google Scholar

Y. Hoshiet al., “Resting hypofrontality in schizophrenia: a study using near-infrared time-resolved spectroscopy,” Schizophr. Res. 84(2–3), 411–420 (2006).SCRSEH0920-9964http://dx.doi.org/10.1016/j.schres.2006.03.010Google Scholar

A. LiebertA. Kienle, “Light diffusion in N-layered turbid media: frequency and time domains,” J. Biomed. Opt. 15(2), 025002 (2010).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3368682Google Scholar

A. Kienleet al., “Investigation of two-layered turbid media with time-resolved reflectance,” Appl. Opt. 37(28), 6852–6862 (1998).APOPAI0003-6935http://dx.doi.org/10.1364/AO.37.006852Google Scholar

A. H. Hielscheret al., “Time-resolved photon emission from layered turbid media,” Appl. Opt. 35(4), 719–728 (1996).APOPAI0003-6935http://dx.doi.org/10.1364/AO.35.000719Google Scholar

A. Kienleet al., “Noninvasive determination of the optical properties of two-layered turbid media,” Appl. Opt. 37(5), 779–791 (1998).APOPAI0003-6935http://dx.doi.org/10.1364/AO.37.000779Google Scholar

F. Martelliet al., “Analytical approximate solutions of the time-domain diffusion equation in layered slabs,” J. Opt. Soc. Am. A 19(1), 71–80 (2002).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.19.000071Google Scholar

F. Martelliet al., “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 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.002037Google Scholar

M. ShimadaY. HoshiY. Yamada, “Simple algorithm for the measurement of absorption coefficients of a two-layered medium by spatially resolved and time-resolved reflectance,” Appl. Opt. 44(35), 7554–7563 (2005).APOPAI0003-6935http://dx.doi.org/10.1364/AO.44.007554Google Scholar

J. Steinbrinket al., “Determining changes in NIR absorption using a layered model of the human head,” Phys. Med. Biol. 46(3), 879–896 (2001).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/46/3/320Google Scholar

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

F. GaoP. PouletY. Yamada, “Simultaneous mapping of absorption and scattering coefficients from a three-dimensional model of time-resolved optical tomography,” Appl. Opt. 39(31), 5898–5910 (2000).APOPAI0003-6935http://dx.doi.org/10.1364/AO.39.005898Google Scholar

B. W. Pogueet al., “Initial assessment of a simple system for frequency domain diffuse optical tomography,” Phys. Med. Biol. 40(10), 1709–1729 (1995).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/40/10/011Google Scholar

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

B. W. Zeffet al., “Retinotopic mapping of adult human visual cortex with high-density diffuse optical tomography,” Proc. Natl. Acad. Aci. U. S. A. 104(29), 12169–12174 (2007).Google Scholar

A. CharetteJ. BoulangerH. K. Kim, “An overview on recent radiation transport algorithm development for optical tomography imaging,” JQSRT 109(17–18), 2743–2766 (2008).JQSRAE0022-4073http://dx.doi.org/10.1016/j.jqsrt.2008.06.007Google Scholar

C. Satoet al., “Extraction of depth-dependent signals from time-resolved reflectance in layered turbid media,” J. Biomed. Opt. 10(6), 064008 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2136312Google Scholar

Y. NomuraO. HazekiM. Tamura, “Exponential attenuation of light along nonlinear path through the biological model,” Adv. Exp. Med. Biol. 248, 77–80 (1989).AEMBAP0065-2598http://dx.doi.org/10.1007/978-1-4684-5643-1Google Scholar

Y. TsuchiyaT. Urakami, “Photon migration model for turbid biological medium having various shapes,” Jpn. J. Appl. Phys. 34, L79–L81 (1995).JJAPA50021-4922http://dx.doi.org/10.1143/JJAP.34.L79Google Scholar

L. WangS. L. Jacques, Monte Carlo modeling of light transport in multi-layered tissues in standard C, University of Texas/M.D. Anderson Cancer Center, Houston (1992).Google Scholar

T. J. FarrelM. S. Patterson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo.” Med. Phys. 19(4), 879–888 (1992).MPHYA60094-2405http://dx.doi.org/10.1118/1.596777Google Scholar

R. C. Haskellet al., “ Boundary conditions for the diffusion equation in radiative transfer,” J. Opt. Soc. Am. A Opt. Image Sci. Vis. 11(10), 2727–2741 (1994).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.11.002727Google Scholar

A. KienleM. S. Patterson, “Improved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid medium,” J. Opt. Soc. Am. A Opt. Image Sci. Vis. 14(1), 246–254 (1997).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.14.000246Google Scholar

K. FurutsuY. Yamada, “Diffusion approximation for a dissipative random medium and the applications,” Phys. Rev. E 50(1), 3634–3640 (1994).PLEEE81063-651Xhttp://dx.doi.org/10.1103/PhysRevE.50.3634Google Scholar

K. Levenberg, “A method for the solution of certain problems in least squares, Quart. Appl. Math. 2, 164–168 (1944).QAMAAY0033-569XGoogle Scholar

D. Marquardt, “An algorithm for least-squares estimation of nonlinear parameters,” SIAM J. Appl. Math. 11(2), 431–441 (1963).SMJMAP0036-1399http://dx.doi.org/10.1137/0111030Google Scholar

P. B. Bevington, Data Reduction and Error Analysis for the Physical Science, McGraw-Hill, New York (1969).Google Scholar

V. D. O ConnorD. Phillips, Time-Correlated Single-Photon Counting, Academic Press, London (1985).Google Scholar

F. Bevilacquaet al., “In vivo local determination of tissue optical properties: applications to human brain,” Appl. Opt. 38(22), 4939–4950 (1999).APOPAI0003-6935http://dx.doi.org/10.1364/AO.38.004939Google Scholar

J. H. Choiet al., “Noninvasive determination of the optical properties of adult brain: near-infrared spectroscopy approach,” J. Biomed. Opt. 9(1), 221–229 (2004).JBOPFO1083-3668http://dx.doi.org/10.1117/1.1628242Google Scholar

M. Okamotoet al., “Three-dimensional probabilistic anatomical cranio-cerebral correlation via the international 10–20 system oriented for transcranial functional brain mapping,” Neuroimage 21(1), 99–111 (2004).NEIMEF1053-8119http://dx.doi.org/10.1016/j.neuroimage.2003.08.026Google Scholar

H. Zhaoet al., “Maps of optical differential pathlength factor of human adult forehead, somatosensory motor, occipital regions at multi-wavelengths in NIR,” Phys. Med. Biol. 47(1), 2075–2093 (2002).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/47/12/306Google Scholar

A. Liebertet al., “Time-resolved multidistance near-infrared spectroscopy of the adult head: intracerebral and extracerebral absorption changes from moments of distribution of times of flight of photons,” Appl. Opt. 43(15), 3037–3047 (2004).APOPAI0003-6935http://dx.doi.org/10.1364/AO.43.003037Google Scholar

D. A. BoasA. M. Dale, “Simulation study of magnetic resonance imaging-guided cortically constrained diffuse optical tomography of human brain function,” Appl. Opt. 44(9), 1957–1968 (2005).APOPAI0003-6935http://dx.doi.org/10.1364/AO.44.001957Google Scholar

Y. Hoshiet al., “Reevaluation of near-infrared light propagation in the adult human head: implications for functional near-infrared spectroscopy,” J. Biomed. Opt. 10(6), 064032 (2005).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2142325Google Scholar

E. OkadaD. T. Delpy, “Near-infrared light propagation in an adult head model. I. Modeling of low-level scattering in the cerebrospinal fluid layer,” Appl. Opt. 42(16), 2906–2914 (2003).APOPAI0003-6935http://dx.doi.org/10.1364/AO.42.002906Google Scholar

J. Johanssonet al., “Time-resolved NIR/Vis spectroscopy for analysis of solids: Pharmaceutical tablets.” Appl. Spectro. 56(6), 725–731 (2002).APSPA40003-7028http://dx.doi.org/10.1366/000370202760077676Google Scholar

T. J. KennettW. V. PrestwichA. Robertson, “Baysian deconvolution II: noise properties.” Nucl. Instr. Meth. 151, 293–301 (1978).NUIMAL0029-554Xhttp://dx.doi.org/10.1016/0029-554X(78)90503-7Google Scholar