EP1136811B1

Method and device for measuring internal information of scattering absorber

Abstract

This record has no abstract on file.

EP1136811B1, drawing sheet 1 — Sheet 1 of 70

Term

Term ended

Expired 5 November 2019, 6.9 years ago.

12 claims: 2 independent, 10 dependent

1
A measuring method of internal information of a scattering medium, comprising:a light injecting step of injecting pulsed light of two or more predetermined wavelengths into a scattering medium at a light injection position;a light detecting step of detecting the light of said two or more predetermined wavelengths having propagated inside said scattering medium, at a photodetection position to acquire a photodetection signal;characterised by a signal processing step of acquiring waveform data indicating a temporal change of intensity of the detected light, based on said photodetection signal;a mean pathlength and variance computing step of performing an operation to compute a mean pathlength of plural photons composing said detected light, and a variance, based on said waveform data;and an absorption coefficient difference calculating step of calculating a difference between absorption coefficients at said predetermined wavelengths, based on a predetermined relation holding among said mean pathlength, said variance, and the difference between the absorption coefficients at said two or more predetermined wavelengths.
7
A measuring apparatus of internal information of a scattering medium, comprising:light injecting means (3, 4, 5) for injecting pulsed light of two or more predetermined wavelengths into a scattering medium at a light injection position;light detecting means (7) for detecting the light of said two or more predetermined wavelengths having propagated inside said scattering medium, at a photodetection position to acquire a photodetection signal;characterised by a signal processing means (8) for acquiring waveform data indicating a temporal change of intensity of the detected light, based on said photodetection signal;mean pathlength and variance computing means for performing an operation to compute a mean pathlength of plural photons composing said detected light, and a variance, based on said waveform data;and absorption coefficient difference calculating means (9) for calculating a difference between absorption coefficients at said predetermined wavelengths, based on a predetermined relation holding among said mean pathlength, said variance, and the difference between the absorption coefficients at said two or more predetermined wavelengths.

Independent claims2

112 paragraphs, as filed

Technical Field

The present invention relates to methods and apparatus for measuring internal information of scattering media, e.g., an absorption coefficient or a concentration of an absorber.

Background Art

Methods of measuring an absorption coefficient of a scattering medium being a medium to be measured, a concentration of an absorber or the like, based on the Microscopic Beer-Lambert Law (hereinafter referred to as "MBL law"), include, for example, the methods disclosed in Japanese Patent Applications Laid-Open No. H08-94517, No. H10-73481, and EPO83427 by the inventors. These methods based on the MBL law possess the significant feature of being theoretically free of the influence of ① the shape of medium, ② boundary conditions, ③ scattering, and so on, and the same analytical formulas can be applied to media having any medium shape, any boundary conditions, and various scattering characteristics as long as there is no reentry of photon into the scattering medium.

The measuring methods based on the MBL law can be roughly classified under four types at present. They are (1) Time Resolved Spectroscopy (hereinafter referred to as "TRS method") making use of the time-resolved profile of re-emission, (2) Time Integrated Spectroscopy (hereinafter referred to as "TIS method") making use of the time integral value of the time-resolved profile of re-emission and the mean pathlength, (3) Time Gating Spectroscopy (hereinafter referred to as "TGS method") making use of portions cut out of the time-resolved profile by gates, and (4) Phase Modulation Spectroscopy (hereinafter referred to as "PMS method") making use of modulated light. Among these, (2) the TIS method and (4) the PMS method utilizing all re-emission are advantageous from the viewpoint that optical attenuation due to scattering is great in living tissues and it is thus important to use as much re-emission as possible in practical use. These two measuring methods are in the relation of Fourier transform with each other.

Disclosure of the Invention

While the measuring methods based on the MBL law have the many advantages as described above, the measurement accuracy thereof, however, is not sufficient yet for utilization and application in the wide range. For example, there was the problem that, in the case of a living tissue being a measured object, various individual differences, such as color of skin, presence/absence of hair, etc., affected the absolute value of light intensity or the like, so as to be the cause of degradation of the measurement accuracy. Further, the wavelength dependence of scattering coefficient also degrades the measurement accuracy.

There was another problem that reduction was insufficient in the computation time of analysis during measurement and it was thus difficult to implement real-time measurement.

The present invention has been accomplished in view of the above problems and an object of the invention is to provide measuring methods and apparatus of internal information of scattering medium capable of measurement with higher accuracy and at higher speed, when compared with the conventional measuring methods based on the MBL law.

The inventor has conducted intensive and extensive research in order to accomplish the above object and finally found that it became feasible to implement highly accurate measurement by measuring the mean pathlength and variance, particularly, for pulsed light of a plurality of different wavelength components, by the Mean and Variance based Spectroscopy (hereinafter referred to as "MVS method") utilizing the mean pathlength and variance, or physical quantities equivalent thereto, without use of information such as absolute values of light intensity or a ratio thereof, and to readily put the wavelength dependence of scattering coefficient into analytical formulas, thus accomplishing the present invention. Likewise, the inventor found that it also became feasible to implement highly accurate measurement by measuring the group delay for modulated light of a plurality of different wavelength components and the second partial derivative of logarithm of amplitude with respect to modulation frequency, thus accomplishing the present invention.

Specifically, a first measuring method of internal information of a scattering medium according to the present invention is a method comprising (1) a light injecting step of injecting pulsed light of two or more predetermined wavelengths into a scattering medium at a light injection position, (2) a light detecting step of detecting the light of the two or more predetermined wavelengths having propagated inside the scattering medium, at a photodetection position to acquire a photodetection signal, (3) a signal processing step of acquiring waveform data indicating a temporal change of intensity of the detected light, based on the photodetection signal, (4) a mean pathlength and variance computing step of performing an operation to compute a mean pathlength of plural photons composing the detected light, and a variance, based on the waveform data, and (5) an absorption coefficient difference calculating step of calculating a difference between absorption coefficients at the predetermined wavelengths, based on a predetermined relation holding among the mean pathlength, the variance, and the difference between the absorption coefficients at the two or more predetermined wavelengths.

A first measuring apparatus of internal information of a scattering medium according to the present invention is an apparatus comprising (1) light injecting means for injecting pulsed light of two or more predetermined wavelengths into a scattering medium at a light injection position, (2) light detecting means for detecting the light of the two or more predetermined wavelengths having propagated inside the scattering medium, at a photodetection position to acquire a photodetection signal, (3) signal processing means for acquiring waveform data indicating a temporal change of intensity of the detected light, based on the photodetection signal, (4) mean pathlength and variance computing means for performing an operation to compute a mean pathlength of plural photons composing the detected light, and a variance, based on the waveform data, and (5) absorption coefficient difference calculating means for calculating a difference between absorption coefficients at the predetermined wavelengths, based on a predetermined relation holding among the mean pathlength, the variance, and the difference between the absorption coefficients at the two or more predetermined wavelengths.

The above first method and apparatus according to the present invention are based on a TIMVS method, which is an MVS method in which the analysis is carried out in the time domain by the time integrated spectroscopy (TIS method). The TIMVS method using the mean pathlength and variance for the light of plural wavelength components in this way presents very significant advantages in practical use, including ① being not influenced by the individual differences and absolute values of intensity of incident light dependent upon wavelengths and positions, ② extremely simplifying the quantitative formula involving the wavelength dependence of scattering coefficient, ③ decreasing errors in determination of the zero point (t = 0) of the time axis in the time-resolved spectroscopy, and so on, in addition to the advantages presented by the conventional measuring methods based on the MBL law.

In the method and apparatus, it also becomes feasible to greatly reduce the measurement and analysis

In the method and apparatus, it also becomes feasible to greatly reduce the measurement and analysis time and implement real-time measurement, by applying the Simple Subtraction Method (SSM) developed by the inventors, which permits the mean pathlength of detected photons (the center of gravity of a time-resolved profile) and the variance to be computed at high speed by calculating the moment of the time-resolved profile with a computer. This SSM is described, for example, in Japanese Patent Application Laid-Open No. H09-61343.

Brief Description of the Drawings

<ul id="ul0001" list-style="none" compact="compact"><li>Fig. 1 is a graph to show the relationship between mean pathlength and absorption coefficient for different scattering coefficients.</li><li>Fig. 2 is a schematic diagram to show an embodiment of the measuring apparatus of internal information of scattering medium according to the present invention.</li><li>Fig. 3 is a schematic diagram to show an example of preferred specific structure of the apparatus illustrated in Fig. 2.</li><li>Fig. 4 is a flowchart to show an embodiment of the measuring method of internal information of scattering medium according to the present invention.</li><li>Fig. 5 is a flowchart to show a modification of the measuring method of internal information of scattering medium not embodying the present invention.</li><li>Fig. 6 is a schematic diagram to show another embodiment of the measuring apparatus of internal information of scattering medium according to the present invention.</li><li>Fig. 7 is a schematic diagram to show an example of preferred specific structure of the apparatus illustrated in Fig. 6.</li><li>Fig. 8 is a graph to show absorption spectra of hemoglobin.</li><li>Fig. 9 is a graph to show the relationship between absorption coefficient and mean pathlength of a phantom.</li></ul>

Best Mode for Carrying out the Invention

Detailed below with reference to the drawings are preferred embodiments of the measuring methods and apparatus of internal information of scattering medium according to the present invention. In the description of the drawings the same elements will be denoted by the same reference symbols and redundant description will be omitted.

First, the principles of the present invention will be described.

(Principles of the invention)

Calculation of absorption coefficient difference by TIS method

The survival probability of photons propagating in a zigzag path inside a scattering medium is the exponential of the product of the zigzag pathlength 1 and the scattering coefficient µa of the medium (scattering medium), exp(-µal). Namely, the attenuation is expressed by the product µal of the zigzag photon pathlength 1 and the absorption coefficient µa. At this time, the impulse response h(t) of the scattering medium is a time-causal function, which is expressed as follows.<maths id="math0001" num=""><math display="block"><mrow><mtext mathvariant="italic">h</mtext><mtext>(</mtext><mtext mathvariant="italic">t</mtext><mtext>)=</mtext><mtext mathvariant="italic">h</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,</mtext><mtext mathvariant="italic">t</mtext><mtext>)=</mtext><mtext mathvariant="italic">s</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><mtext>,</mtext><mtext mathvariant="italic">t</mtext><msub><mrow><mtext>)exp(-µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext mathvariant="italic">ct</mtext><mtext>) (1.1</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0001.tif" /></maths><maths id="math0002" num=""><math display="block"><mrow><mfrac><mrow><mtext>∂</mtext></mrow><mrow><msub><mrow><mtext>∂µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext>ln </mtext><mtext mathvariant="italic">h</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,</mtext><mtext mathvariant="italic">t</mtext><mtext>)=-</mtext><mtext mathvariant="italic">ct</mtext><mtext>=-ℓ (1.1</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0002.tif" /></maths><maths id="math0003" num=""><img file="EP1136811B1_D0003.tif" /></maths> Here µs and µa represent an anisotropic scattering coefficient and an absorption coefficient, c the speed of photons in the medium, t the time of flight, and 1 the photon pathlength (distance of flight). The function s(µs,t) is the response at the absorption coefficient µa = 0. The time of flight t can be measured by time-resolved spectroscopy. The speed of photons c is determined by the refractive index of the scattering medium and can be deemed as a constant, for example, in the case of living tissues. The above facts are called the MBL law. A transport scattering coefficient µ's (also called an equivalent scattering coefficient) used hereinafter is expressed as µ's = (1 - g)µs, using µs and an average g of cosine of scattering angles.

The temporal integration I(µs,µa) of the impulse response h(t) of the scattering medium is given as follows.<maths id="math0004" num=""><img file="EP1136811B1_D0004.tif" /></maths><maths id="math0005" num=""><img file="EP1136811B1_D0005.tif" /></maths><maths id="math0006" num=""><img file="EP1136811B1_D0006.tif" /></maths> L(µs,µa) = c<t> in above Eq. (1.2b) represents an average of pathlengths of detected photons (which is also called a mean pathlength). On this occasion, <t> indicates the center of gravity of the impulse response profile (mean time of flight of detected photons) and can be readily calculated by an operation with a computer (to compute the moment) from the temporal profile of the impulse response.

Next, absorption coefficient dependency of the mean pathlength L(µs,µa) is derived as follows.<maths id="math0007" num=""><img file="EP1136811B1_D0007.tif" /></maths> Here σ2 represents the variance of pathlength 1, which is equivalent to a result of partially differentiating the mean pathlength L(µs,µa) with respect to µa and flipping the sign, or to the second partial derivative of the temporal integration I(µs,µa) with respect to µa. This variance σ2(µs,µa) can be readily calculated by an operation with a computer from the temporal profile of the impulse response, similarly as in the case of the above mean pathlength L(µs,µa).

Likewise, we can find the third partial derivative of the temporal integration I(µs,µa) with respect to µa, i.e., the second partial derivative of the mean pathlength L(µs,µa) with respect to µa, and this value provides information about profile distortion. From the mathematical aspect, if there exists the mth partial derivative, there always exist the (m - 1)th and lower partial derivatives.

Here, assuming that the scattering coefficient µs is a constant, the Taylor expansion of the mean pathlength L(µs,µa) = L(µa) around µa1 is expressed as follows.<maths id="math0008" num="(1.4)"><math display="block"><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext> + δ)=</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)+δ</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>'(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)+</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2!</mtext></mrow></mfrac><msup><mrow><mtext>δ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>"(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)+·····</mtext></mrow></math><img file="EP1136811B1_D0008.tif" /></maths> In this equation, L' and L" represent the first and second partial derivatives of L with respect to µa. Substituting δ = µa2 - µa1 into the above equation, we obtain the following.<maths id="math0009" num="(1.5)"><math display="block"><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)-</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msup><mrow><mtext>)=δσ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)-</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><msup><mrow><mtext>δ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>''(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)·····</mtext><mspace linebreak="newline" /><msub><mrow><mtext> ≈ (µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext><mfrac><mrow><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msup><mrow><mtext>)+σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac></mrow></math><img file="EP1136811B1_D0009.tif" /></maths> This means that with the scattering coefficient being constant, a change of absorption coefficient (absorption coefficient difference), (µa2 - µa1), can be calculated using the mean pathlength of the impulse response (the center of gravity of the time-resolved profile) and the variance. This novel finding will be applied to the measurement of concentration of absorber as described hereinafter.

Temporal widths of incident light pulses used in measurement are finite and bandwidths of amplifiers and counting circuits are also finite. Therefore, the temporal profile (observed waveform or observed values) obtained in actual measurement is the convolution between the impulse response of the scattering medium and the impulse response of the measuring system (also called an instrumental function).

There are two approaches below as means for determining the true mean pathlength and variance of the impulse response of the scattering medium by eliminating the influence of characteristics of the measuring apparatus from measured values. The first approach is the deconvolution method known well. This is a method of obtaining the impulse response by deconvolution of the measured values with the instrumental function and then determining the mean pathlength and variance from the resultant profile. The second approach is a method of separately determining the mean pathlength and variance in the instrumental function and the mean pathlength and variance in the observed profile and then finding the mean pathlength and variance of the impulse response of the scattering medium from these values. In this case, each of the mean pathlength and variance in the impulse response is given by a difference between corresponding values in the observed profile and in the instrumental function.

The observed profile o(t) is expressed as follows, using the impulse response (true re-emission profile) h(t) of the medium and the impulse response of the measuring system (instrumental function) i(t).<maths id="math0010" num="(2.1)"><math display="block"><mrow><mtext mathvariant="italic">o</mtext><mtext>(</mtext><mtext mathvariant="italic">t</mtext><mtext>) = </mtext><mtext mathvariant="italic">i</mtext><mtext>(</mtext><mtext mathvariant="italic">t</mtext><mtext>) ⊗ </mtext><mtext mathvariant="italic">h(t)</mtext></mrow></math><img file="EP1136811B1_D0010.tif" /></maths> In the equation the symbol ⊗ represents the convolution operation. Let µo, µi, and µh denote centers of gravity of the profiles o(t), i(t), and h(t), respectively, and σo2, σi2, and σh2 denote their respective variances. Let us define generating functions of the above respective profiles as follows.<maths id="math0011" num=""><img file="EP1136811B1_D0011.tif" /></maths><maths id="math0012" num=""><img file="EP1136811B1_D0012.tif" /></maths><maths id="math0013" num=""><img file="EP1136811B1_D0013.tif" /></maths> This definition permits these generating functions to be differentiated any times at s = 0.

With these generating functions, the convolution of the profiles can be expressed by the product of the generating functions, i.e., as follows.<maths id="math0014" num="(2.5)"><math display="block"><mrow><mtext mathvariant="italic">O</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>) = </mtext><mtext mathvariant="italic">I</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)</mtext><mtext mathvariant="italic">H</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0014.tif" /></maths> The center of gravity µo of o(t) is given as follows.<maths id="math0015" num="(2.6)"><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext mathvariant="italic">O</mtext><mtext>'(0)</mtext></mrow><mrow><mtext mathvariant="italic">O</mtext><mtext>(0)</mtext></mrow></mfrac></mrow></math><img file="EP1136811B1_D0015.tif" /></maths> Since the following relation is derived from Eq. (2.5):<maths id="math0016" num="(2.7),"><math display="block"><mrow><mtext mathvariant="italic">O'</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>) = </mtext><mtext mathvariant="italic">I'</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)</mtext><mtext mathvariant="italic">H</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)+</mtext><mtext mathvariant="italic">I</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)</mtext><mtext mathvariant="italic">H'</mtext><mtext>(</mtext><mtext mathvariant="italic">s</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0016.tif" /></maths> we obtain the following relation therefrom.<maths id="math0017" num="(2.8)"><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><msub><mrow><mtext> = µ</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><msub><mrow><mtext> + µ</mtext></mrow><mrow><mtext mathvariant="italic">h</mtext></mrow></msub></mrow></math><img file="EP1136811B1_D0017.tif" /></maths> Since the variance is given by the following equation:<maths id="math0018" num="(2.9)."><math display="block"><mrow><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> = </mtext><mfrac><mrow><mtext mathvariant="italic">O</mtext><mtext>''(0) </mtext></mrow><mrow><mtext mathvariant="italic">O</mtext><mtext>(0)</mtext></mrow></mfrac><mtext> - </mtext><msup><mrow><mfenced open="(" close=")"><mrow><mfrac><mrow><mtext mathvariant="italic">O</mtext><mtext>'(0)</mtext></mrow><mrow><mtext mathvariant="italic">O</mtext><mtext>(0)</mtext></mrow></mfrac></mrow></mfenced></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1136811B1_D0018.tif" /></maths> we can derive the following relation therefrom.<maths id="math0019" num="(2.10)"><math display="block"><mrow><msub><mrow><mtext>σ</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>=σ</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>+σ</mtext></mrow><mrow><mtext mathvariant="italic">h</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></math><img file="EP1136811B1_D0019.tif" /></maths> Likewise, similar relations are also attained as to the third and higher moments. Accordingly, each moment for the impulse response can be quickly calculated with a computer from the preliminarily measured instrumental function and the observed profile of the medium. The procedures described above are the approach of determining the mean pathlength and variance of the impulse response and the absorption coefficient difference from the measured profile by the TIS method.

Calculation of absorption coefficient difference by PMS method

Next described is the measuring method using intensity-modulated light. The system function H(ω) indicating the frequency response of the medium is expressed by the Fourier transform of the impulse response h(t) as follows.<maths id="math0020" num=""><img file="EP1136811B1_D0020.tif" /></maths> In this equation R and X represent the real part and imaginary part, respectively, A and φ the amplitude and phase delay, respectively, and these can be readily measured by a lock-in amplifier or the like.

Then the following relations (the Cauchy-Riemann equations) hold.<maths id="math0021" num=""><math display="block"><mrow><mfrac><mrow><mtext>∂</mtext><mtext mathvariant="italic">R</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂</mtext><mtext mathvariant="italic">c</mtext><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext>= </mtext><mfrac><mrow><mtext>∂</mtext><mtext mathvariant="italic">X</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext> (3.2</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0021.tif" /></maths><maths id="math0022" num=""><math display="block"><mrow><mfrac><mrow><mtext>∂</mtext><mtext mathvariant="italic">R</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext> = </mtext><mfrac><mrow><mtext>∂</mtext><mtext mathvariant="italic">X</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂</mtext><mtext mathvariant="italic">c</mtext><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext> (3.2</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0022.tif" /></maths> Thus it is seen that the system function H(ω) is a regular function. From Eqs. (3.2a) and (3.2b), we can further derive the following relations.<maths id="math0023" num=""><math display="block"><mrow><mfrac><mrow><mtext>∂ln</mtext><mtext mathvariant="italic">A</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂</mtext><mtext mathvariant="italic">c</mtext><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext>= - </mtext><mfrac><mrow><msub><mrow><mtext>∂φ(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext> (3.2</mtext><mtext mathvariant="italic">c</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0023.tif" /></maths><maths id="math0024" num=""><math display="block"><mrow><mfrac><mrow><mtext>∂ln</mtext><mtext mathvariant="italic">A</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext>= </mtext><mfrac><mrow><msub><mrow><mtext>∂φ(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂</mtext><mtext mathvariant="italic">c</mtext><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext> (3.2</mtext><mtext mathvariant="italic">d</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0024.tif" /></maths> Then we can derive the following equation, for example, from Eq. (3.2c).<maths id="math0025" num=""><img file="EP1136811B1_D0025.tif" /></maths> This Eq. (3.3) is similar to aforementioned Eq. (1.2c) and the left side of Eq. (3.3) is observable. The integrand in the first term of the right side is the group delay and corresponds to the mean pathlength L(µs,µa) described previously. This group delay is approximated as follows, using the phase delays φ1 and φ2 at two modulation frequencies ω1 and ω2.<maths id="math0026" num="(3.4)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext>∂φ(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext>= </mtext><mfrac><mrow><msub><mrow><mtext>φ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>-φ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext>ω</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>-ω</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac></mrow></math><img file="EP1136811B1_D0026.tif" /></maths> The second term of the right side in Eq. (3.3) is a value at the absorption coefficient µa = 0. Further, the group delay ∂φ(µs,µa,ω)/∂ω can be approximated as follows when ω « cµa.<maths id="math0027" num=""><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext>∂φ(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext></mtext><mfrac><mrow><mtext>φ</mtext></mrow><mrow><mtext>ω</mtext></mrow></mfrac></mrow></math><img file="EP1136811B1_D0027.tif" /></maths> Here φ/ω represents the phase delay.

We can obtain the absorption coefficient dependency of the group delay as follows.<maths id="math0028" num="(3.5)"><math display="block"><mrow><mfrac><mrow><mtext>∂</mtext></mrow><mrow><msub><mrow><mtext>∂µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext></mtext><mfrac><mrow><msub><mrow><mtext>∂φ(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><mtext>∂ω</mtext></mrow></mfrac><mtext> = -</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext mathvariant="italic">c</mtext></mrow></mfrac><mtext></mtext><mfrac><mrow><msup><mrow><mtext>∂</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>ln</mtext><mtext mathvariant="italic">A</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><msub><mrow><mtext>∂µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> = </mtext><mtext mathvariant="italic">c</mtext><mtext></mtext><mfrac><mrow><msup><mrow><mtext>∂</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>ln</mtext><mtext mathvariant="italic">A</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>,ω)</mtext></mrow><mrow><msup><mrow><mtext>∂ω</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP1136811B1_D0028.tif" /></maths> Here ∂2lnA/∂ω2 in the right side of Eq. (3.5) can be readily measured using three modulation frequencies. This corresponds to the variance previously defined.

In the above discussion Eq. (3.3) to Eq. (3.5) were obtained from Eq. (3.2c), but it is also possible to obtain similar relations from Eqs. (3.2a), (3.2b), and (3.2d). The PMS method normally obviates the need for the aforementioned deconvolution operation for obtaining the true profile.

Just as in the case of the TIS method, the relation between the group delay and the change of absorption coefficient (absorption coefficient difference) is obtained as follows, assuming that the scattering coefficient µs is a constant.<maths id="math0029" num=""><img file="EP1136811B1_D0029.tif" /></maths> Accordingly, it is also feasible to calculate the change of absorption coefficient (absorption coefficient difference) with the scattering coefficient being constant, from actually measured values of the group delay and the absorption coefficient dependency by the PMS method, similarly as in the case of the aforementioned TIS method.

Calculation of absorber concentration by MVS method

Next described is the mean and variance based spectroscopy (MVS method) to quantify the absorber concentration from values of the mean pathlength and variance of the impulse response. In the following, the transport scattering coefficient µ's(= (1 - g)µs), which is easier to measure in general, is used instead of the scattering coefficient µs. For easier description, we will consider dual-wavelength spectroscopy of a scattering medium containing one absorber component and the principal component of water, the wavelengths used in the measurement are λ1 and λ2, and the optical constants at the respective wavelengths will be accompanied by subscripts 1 and 2, respectively. At this time, the following relations hold between the concentration C of the absorber and the absorption coefficients µa1 and µa2 of the medium at the wavelengths λ1 and λ2.<maths id="math0030" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext></mtext><mtext mathvariant="italic">= C</mtext><msub><mrow><mtext>ε</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> + µ</mtext></mrow><mrow><mtext mathvariant="italic">w</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>, µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext> = </mtext><mtext mathvariant="italic">C</mtext><msub><mrow><mtext>ε</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext> + µ</mtext></mrow><mrow><mtext mathvariant="italic">w</mtext><mtext>2</mtext></mrow></msub><mtext> (4.1</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0030.tif" /></maths><maths id="math0031" num=""><math display="block"><mrow><mtext mathvariant="italic">C</mtext><mtext> = </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><msub><mrow><mtext>ε</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - ε</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><msub><mrow><mtext>[(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>)-(µ</mtext></mrow><mrow><mtext mathvariant="italic">w</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">w</mtext><mtext>1</mtext></mrow></msub><mtext>)] (4.1</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0031.tif" /></maths> Here ε1 and ε2 indicate extinction coefficients (or absorption coefficients) for unit concentration of the absorber at the wavelengths λ1 and λ2, which are, for example, molar extinction coefficients for example, and µw1 and µw2 represent absorption coefficients of water at the wavelengths λ1 and λ2. Therefore, it is possible to quantify the concentration C of the absorber by obtaining the difference (µa2 - µa1) between the absorption coefficients from measured profiles.

In execution of the above spectroscopy, ordinary media demonstrate different values of scattering coefficients at the two wavelengths (λ1 and λ2). Then this difference between the scattering coefficients complicates the algorithm for quantification of the absorber concentration. In the following we will first find a relation between the mean pathlength and the transport scattering coefficient and then explain a new method of quantifying the absorber concentration by use of the relation.

According to the photon diffusion equation, the mean pathlength L(µ's,µa) in the case of reflection measurement, which will be simply abbreviated as Lρ, is expressed as follows.<maths id="math0032" num=""><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext mathvariant="italic">ρ</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><msub><mrow><mtext>3µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext></mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">r</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><msub><mrow><mtext>1+ρµ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext></mrow></msub></mrow></mfrac><mtext> (4.2</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0032.tif" /></maths><maths id="math0033" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext></mrow></msub><msub><mrow><mtext>=(3µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>1/2</mtext></mrow></msup><mtext> (4.2</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0033.tif" /></maths> The symbol ρ indicates the source-detector distance, the distance between a light injection position and a photodetection position.

Now let us consider a ratio of mean pathlengths in a case wherein the absorption coefficients are equal, µa, but the transport scattering coefficients are µ's1 and µ's2 (which are equal just to the transport scattering coefficients at the wavelengths λ1 and λ2) (though this will not usually occur). The ratio is given as follows.<maths id="math0034" num="(4.3)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext mathvariant="italic">ρ</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext mathvariant="italic">ρ</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac><mtext> = </mtext><mfrac><mrow><msub><mrow><mtext>µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext>1+ρµ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext><mtext>2</mtext></mrow></msub></mrow></mfrac><mtext></mtext><mfrac><mrow><msub><mrow><mtext>1+ρµ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext><mtext>1</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext>µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub></mrow></mfrac></mrow></math><img file="EP1136811B1_D0034.tif" /></maths> Here µeff represents an effective attenuation coefficient and is given as follows.<maths id="math0035" num="(4.4)"><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext><mtext>1</mtext></mrow></msub><mtext>=(3</mtext><mtext mathvariant="italic">µ</mtext><msub><mrow><mtext>'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext> µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>1/2</mtext></mrow></msup><msub><mrow><mtext>, µ</mtext></mrow><mrow><mtext mathvariant="italic">eff2</mtext></mrow></msub><msub><mrow><mtext>=(3µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>1/2</mtext></mrow></msup></mrow></math><img file="EP1136811B1_D0035.tif" /></maths> However, a photon diffusion constant independent of the absorption coefficient was used herein. When the measurement conditions are ordinary conditions, i.e., ρ > 20 mm, µ's > 0.8 mm, and µa » 0.001 mm-1, Eq. (4.3) is approximated as follows.<maths id="math0036" num="(4.5)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext mathvariant="italic">ρ</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext mathvariant="italic">ρ</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac><mtext> =</mtext><msup><mrow><mfenced open="(" close=")"><mrow><mfrac><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub></mrow><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mtext>1/2</mtext></mrow></msup></mrow></math><img file="EP1136811B1_D0036.tif" /></maths> Namely, the ratio of mean pathlengths in the case wherein the absorption coefficients are equal but the transport scattering coefficients are different is a constant determined by a ratio of the transport scattering coefficients, independent of the absorption coefficients. Since in the measurement of actual media change of the transport scattering coefficient with change of wavelength normally also varies the absorption coefficient, it is difficult to directly measure the ratio of Lρ represented by Eq. (4.3).

The approximation of Eq. (4.5) also holds similarly for the ratio of mean pathlengths under almost equal conditions to the above, i.e., in the case of ordinary transmission measurement. Fig. 1 schematically shows by L1 and L2, an example of the relation between the mean pathlength L and the absorption coefficient µa for different transport scattering coefficients.

Eq. (4.5) above is the relation obtained based on the photon diffusion approximation, but the validity of this relation for actual media was proved by the Monte Carlo calculation to simulate scattering media and by experiments with a tissue-like phantom.

To generalize the above, the ratio of mean pathlengths at the wavelengths λ1 and λ2 is approximated as follows and defined by k'.<maths id="math0037" num="(4.6)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac><mtext> ≡ </mtext><mfrac><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac><mtext> = </mtext><msup><mrow><mfenced open="(" close=")"><mrow><mfrac><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub></mrow><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mtext>1/2</mtext></mrow></msup><mtext> ≡ </mtext><mtext mathvariant="italic">k</mtext><mtext>'</mtext></mrow></math><img file="EP1136811B1_D0037.tif" /></maths> L1(µa) and L2(µa) used herein indicate the mean pathlengths for a medium having different transport scattering coefficients but an absorption coefficient of an equal value µa. Accordingly, use of the coefficient k' defined here permits the ratio of mean pathlengths to be estimated for the medium having the absorption coefficient of the equal value µa and the different transport scattering coefficients (see Fig. 2). Moreover, since the fact that Eq. (4.6) holds means that curves representing the mean pathlengths, k' × L1(µa) and L2(µa), are laid on each other on the plane (µa,L), the following relation holds for the variances σ2(µ's1,µa) and σ2(µ's2,µa) at µa of the curves L1(µa) and L2(µa) indicating the mean pathlengths.<maths id="math0038" num="(4.7)"><math display="block"><mrow><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msub><mrow><mtext>) ≡ σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>) = </mtext><mtext mathvariant="italic">k</mtext><msup><mrow><mtext>'σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>) ≡ </mtext><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>'σ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0038.tif" /></maths>

The above knowledge is very significant in the following sense. Namely, use of the coefficient k' defined in Eq. (4.6) permits us to normalize the mean pathlength L(µ's,µa) for media with different transport scattering coefficients. As a consequence, we can unitarily handle the mean pathlengths for the media with different transport scattering coefficients.

Now let us describe how to determine the concentration C of the absorber, using this relation.

The introduction of the coefficient k' as described above reduces the spectroscopy of the media exhibiting the wavelength dependence of the transport scattering coefficient to measurement of concentration of absorber with the transport scattering coefficient being constant. As a result, we need only to consider one curve indicating the mean pathlength on the plane (µa,L) and can apply Eq. (1.5) to the measurement of absorber concentration. Now, referring to Eq. (1.5) with use of the mean pathlength L2 = k'L1 normalized by k', the following relation holds as to the difference between absorption coefficients.<maths id="math0039" num=""><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)-</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>) = (µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>-µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext><mfrac><mrow><msub><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>)+σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac></mrow></math><img file="EP1136811B1_D0039.tif" /></maths> Since there are the following relations:<maths id="math0040" num=""><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>) = L(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) = </mtext><mtext mathvariant="italic">k</mtext><mtext>'</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) = </mtext><mtext mathvariant="italic">k</mtext><mtext>'</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0040.tif" /></maths><maths id="math0041" num=""><math display="block"><mrow><msub><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msup><mrow><mtext>) = σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) = </mtext><mtext mathvariant="italic">k</mtext><msup><mrow><mtext>'σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext mathvariant="italic">L</mtext><msub><mrow><mtext>(µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>,µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) = </mtext><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>'σ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0041.tif" /></maths> the following relation finally holds.<maths id="math0042" num="(4.8)"><math display="block"><mrow><mtext mathvariant="italic">k</mtext><mtext>'</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) - </mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>) = (µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>-µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext><mfrac><mrow><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>'σ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>)+σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext>) </mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac></mrow></math><img file="EP1136811B1_D0042.tif" /></maths> This means that the absorption coefficient difference can be quantified using the mean pathlengths L1(µa1), L2(µa2) and the variances σ12(µa1), σ22(µa2), which are obtained from measured values.

From the above, the absorber concentration C is finally obtained as follows.<maths id="math0043" num=""><img file="EP1136811B1_D0043.tif" /></maths> This yields the equation for quantifying the absorber concentration C by dual-wavelength spectroscopy. Since the mean pathlength term is (k'L1 - L2) herein, the problem of error in determination of the zero point of the time axis (t = 0) in time-resolved spectroscopy is largely relaxed in comparison with (L1 + L2) in the conventional method. This effect is particularly great in the case of living tissues with µ's ≒ 1.

In the case of the phase modulation spectroscopy (PMS method), the concentration is given as follows in the same manner as above.<maths id="math0044" num=""><img file="EP1136811B1_D0044.tif" /></maths> However, the notation follows the following definition.<maths id="math0045" num=""><img file="EP1136811B1_D0045.tif" /></maths><maths id="math0046" num=""><img file="EP1136811B1_D0046.tif" /></maths> Since the group delay term is also (k'φ'1 - φ'2) in above Eq. (4.10), the problem of error in determination of the zero point of the phase (φ = 0) in the phase modulation spectroscopy is largely relaxed in comparison with (φ1 + φ2) in the conventional method. This effect is also particularly great in the case of the living tissues with µ's ≒ 1, as above.

It is apparent from the above that the spectroscopy (MVS method) according to the present invention with use of the mean pathlength and variance, or physical quantities equivalent thereto, presents tremendous advantages, including ① there will arise no problem from the individual differences and the absolute value of incident photon intensity dependent upon the wavelength and position, ② the quantitative equation is very simple with inclusion of the wavelength dependency of scattering coefficient, ③ the problem of error is relaxed in the determination of the zero point of the time axis (t = 0) in the time-resolved spectroscopy, ④ the problem of error is relaxed in the determination of the zero point of the phase in the phase modulation spectroscopy, and so on, in addition to those of the conventional measuring methods based on the MBL law.

Eq. (4.8) above can also be derived as follows. First, the following is yielded from Eq. (1.3).<maths id="math0047" num=""><img file="EP1136811B1_D0047.tif" /></maths> Then the difference between mean pathlengths at the wavelengths λ1 and λ2 is given as follows.<maths id="math0048" num=""><img file="EP1136811B1_D0048.tif" /></maths> Assuming Eq. (4.6) and Eq. (4.7), the following relation is obtained.<maths id="math0049" num=""><img file="EP1136811B1_D0049.tif" /></maths> As a result, aforementioned Eq. (4.8b), i.e., the following equation is derived.<maths id="math0050" num=""><math display="block"><mrow><mtext mathvariant="italic">k</mtext><mtext>'</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>) - </mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext>) = (µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext>)</mtext><mfrac><mrow><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>'σ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><msub><mrow><mtext>)+σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext> (4.8</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0050.tif" /></maths>

Next described is a method of quantifying the ratio of transport scattering coefficients. The mean pathlength Lρ(µ's,µa) in the reflection measurement is given by aforementioned Eq. (4.2) and is rewritten as follows by abbreviating it as L.<maths id="math0051" num=""><math display="block"><mrow><mtext mathvariant="italic">L</mtext><mtext> = </mtext><mfrac><mrow><msub><mrow><mtext>3µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext></mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">r</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><msub><mrow><mtext>1+ρµ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext></mrow></msub></mrow></mfrac><mtext> (5.1</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0051.tif" /></maths><maths id="math0052" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext></mrow></msub><msub><mrow><mtext>=(3µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext></mrow></msub><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>1/2</mtext></mrow></msup><mtext> (5.1</mtext><mtext mathvariant="italic">b</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0052.tif" /></maths> In the equation ρ represents the source-detector distance. Thus the variance σ2 is yielded as follows.<maths id="math0053" num=""><math display="block"><mrow><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> = </mtext><mfrac><mrow><msup><mrow><mtext>∂</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> ln</mtext><mtext mathvariant="italic">I</mtext></mrow><mrow><msub><mrow><mtext>∂µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> = - </mtext><mfrac><mrow><mtext>∂</mtext><mtext mathvariant="italic">L</mtext></mrow><mrow><msub><mrow><mtext>∂µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></mrow></mfrac><mtext> = </mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><msub><mrow><mtext>ρµ</mtext></mrow><mrow><mtext mathvariant="italic">eff</mtext></mrow></msub></mrow></mfrac><mtext> (5.2)</mtext></mrow></math><img file="EP1136811B1_D0053.tif" /></maths> Then we obtain the following equations for the transport scattering coefficient µ's and the absorption coefficient µa, by describing them using the mean pathlength L and variance σ2.<maths id="math0054" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">L</mtext><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext></mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><msup><mrow><mtext>σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> + </mtext><msup><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> (5.3</mtext><mtext mathvariant="italic">a</mtext><mtext>)</mtext></mrow></math><img file="EP1136811B1_D0054.tif" /></maths><maths id="math0055" num=""><img file="EP1136811B1_D0055.tif" /></maths> In the ordinary measurement, i.e., when 3µ'sρ » 4, the transport scattering coefficient can be approximated as follows.<maths id="math0056" num=""><img file="EP1136811B1_D0056.tif" /></maths>

From the above, the transport scattering coefficient µ's and the absorption coefficient µa can be quantified using the mean pathlength L and variance σ2 computed from measured values. This method has the most significant feature and advantage of being simple and attains high accuracy for relatively large media, and the measurement accuracy of the transport scattering coefficient µ's and the absorption coefficient µa at this time is approximately 10%. However, this method presents excellent quantification accuracy for the ratio of transport scattering coefficients in measurement with different media or at different wavelengths. Namely, the transport scattering coefficient ratio can be measured with high accuracy and at high speed by the above method. In this case, the aforementioned subtraction method can be utilized for obtaining the mean pathlength L and variance σ2.

Accordingly, the ratio of transport scattering coefficients obtained in measurement at the wavelengths λ1 and λ2, i.e., µ's(λ2)/µ's(λ1), is given by the following equation, using the mean pathlengths L and variances σ2 obtained in measurement at the respective wavelengths.<maths id="math0057" num="(5.5)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext>µ'(λ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><msub><mrow><mtext>µ'(λ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac><mtext> ≡ </mtext><mfrac><mrow><msub><mrow><mtext>µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext>µ'</mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub></mrow></mfrac><mtext> = </mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub></mrow><mrow><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><mtext></mtext><mfrac><mrow><mtext>1+(</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext>/σ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><mtext>1+(</mtext><msub><mrow><mtext mathvariant="italic">L</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext>/σ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP1136811B1_D0057.tif" /></maths> Therefore, aforementioned k' is given as follows.<maths id="math0058" num="(5.6)"><math display="block"><mrow><mtext mathvariant="italic">k</mtext><mtext>'=</mtext><msup><mrow><mfenced open="(" close=")"><mrow><mfrac><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>2</mtext></mrow></msub><mtext></mtext></mrow><mrow><mtext>µ'</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext mathvariant="italic">s</mtext><mtext>1</mtext></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mtext>1/2</mtext></mrow></msup><mtext> = </mtext><msup><mrow><mfenced open="(" close=")"><mrow><mfrac><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>1+(</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>/σ</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>)</mtext><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><mtext mathvariant="italic">L</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>1+(</mtext><mtext mathvariant="italic">L</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>/σ</mtext><msub><mrow><mtext></mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext><msup><mrow><mtext></mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac></mrow></mfenced></mrow><mrow><mtext>1/2</mtext></mrow></msup></mrow></math><img file="EP1136811B1_D0058.tif" /></maths> It is, however, noted that the mean pathlengths and variances are expressed by their abbreviations.

From the above, the absorber concentration C can be determined from measured values, by substituting Eq. (5.6) into Eq. (4.8). If the ratio of transport scattering coefficients is known, this known value may also be used. This ratio of transport scattering coefficients may also be gained by another method.

Preferred embodiments of the measuring methods and apparatus based on the above-stated measurement principles will be described below in detail.

(First Embodiment)

The first embodiment as a preferred embodiment of the present invention will be described with reference to Fig. 2 to Fig. 4. Fig. 2 shows a measuring system according to the present invention to quantify the concentration C of an absorber in a scattering medium 1, using light of two wavelengths. For simplicity of description, let us consider a liquid medium containing one absorber, in which the scattering and absorption characteristics are uniform. The liquid is, for example, water and it is necessary to take absorption of water in the measurement into consideration.

The apparatus illustrated in Fig. 2 is provided with a light guide 3 for injection of light and an output end of the light guide 3 is located at a predetermined position on the surface of the scattering medium 1. A light source 5 is optically coupled through a wavelength selector 4 to an input end of the light guide 3 and pulsed light emitted from the light source 5 is subjected to wavelength selection to the predetermined wavelength λ1 and/or λ2 in the wavelength selector 4. The light is then guided through the light guide 3 to be injected at the position uj into the scattering medium 1.

The time width of this pulsed light can be set to any value short enough to derive the mean pathlength of the impulse response from a photodetection signal and is normally selected in the range of approximately 10 ps to 1 ns. The wavelengths of the light are properly selected according to the scattering medium 1 being a measured object and in general, for example, for living tissues, the wavelengths in the near infrared region of approximately 700 to 900 nm are normally used from the relation between transmittance of living tissue and spectral absorption coefficient of absorber to be quantified. The light source 5 can be selected from a variety of sources including light emitting diodes, laser diodes, various pulse lasers, and so on. This light source 5 may be a combination of two or more sources each emitting light of a single wavelength or light of a narrow band, or may be one simultaneously emitting light of two or more wavelengths. The structure of the light guide 3 and wavelength selector 4 is also adequately modified or set according to the structure of the light source 5 as described. The light source may also be one time-sequentially generating light of two or more wavelengths and in this case the wavelength selector 4 can be omitted.

The apparatus illustrated in Fig. 2 is provided with a light guide 6 for photodetection and an input end of the light guide 6 is located at a predetermined position on the surface of the scattering medium 1. A photodetector 7 is optically coupled to an output end of the light guide 6, and the light having propagated with being scattered inside the scattering medium 1 is guided from the position vk through the light guide 6 to the photodetector 7. The photodetector 7 converts the received signal into a photodetection signal being an electric signal. A signal processing unit 8 is electrically connected to the photodetector 7 and to the light source 5 and this signal processing unit 8 acquires waveform data indicating temporal change of intensity of detected light, based on the photodetection signal. Further, an computation processing unit 9 is electrically connected to the signal processing unit 8 and this computation processing unit 9 performs the operation to compute the mean pathlength and variance of plural photons constituting the detected light, based on the waveform data. The absorption coefficient difference (µa2 - µa1) is quantified according to aforementioned Eq. (4.8), based on these mean pathlength, variance, and the ratio of transport scattering coefficients at two wavelengths, and the concentration C of the absorber is further quantified according to aforementioned Eq. (4.9), based on this absorption coefficient difference, or directly.

The light guide 3 for injection of light, wavelength selector 4, and light source 5 described above constitute the light injecting means according to the present invention; the light guide 6 for detection of light and the photodetector 7 the light detecting means according to the present invention; and the signal processing unit 8 the signal processing means according to the present invention, respectively. The computation processing unit 9 is configured with a plurality of functions, which constitute the mean pathlength and variance computing means (or the group delay and second-partial-derivative-of-amplitude computing means) and the absorption coefficient difference calculating means according to the present invention.

It is desirable to form the part except for the light injecting surface coupled to the light guide 3 and except for the light detecting surface coupled to the light guide 6 in the surface of the scattering medium 1, in structure in which light is absorbed inside but blocked outside. In the case wherein light of plural wavelengths simultaneously propagates with being scattered inside the scattering medium 1, a wavelength selecting filter (not illustrated) may be optionally placed between the photodetector 7 and the light guide 6 to carry out the measurement therewith.

Fig. 3 shows an example of preferred structure of the photodetector 7, signal processing unit 8, and computation processing unit 9. The setup illustrated in Fig. 3 is a configuration for carrying out a highspeed temporal profile measuring method using a method called so-called time-correlated single photon counting. In this setup example, a photomultiplier tube (PMT) is used as the photodetector 7 and the signal processing unit 8 is composed of a constant fraction discriminator (CFD) 21, a time-to-amplitude converter (TAC) 22, and an AD converter (A/D) 23. Then an output signal from the PMT 7 is guided through CFD 21 to TAC 22 to be converted into an analog voltage corresponding to time. This analog voltage is further converted into a digital signal by the AD converter 23. This digital signal is one corresponding to waveform data indicating temporal change of intensity of detected light.

In the computation processing unit 9 illustrated in Fig. 3, a CPU 30 is electrically connected to the light source 5 and to the signal processing unit 8, the CPU 30 controls the timing of photodetection synchronized with injection of light and others, and the waveform data outputted from the signal processing unit 8 is guided to the CPU 30. This CPU 30 also controls or selects the wavelength of incident light and others. Specific techniques are a technique of injecting light of different wavelengths in time sharing manner and a technique of using light simultaneously containing beams of different wavelengths. Specific wavelength selecting means are, for example, an optical beam switch using a mirror, a wavelength switch using a filter, a light switch using an optical switch, and so on.

The computation processing unit 9 illustrated in Fig. 3 is further provided with a program memory 40 storing an operating system (OS) 41 and an internal information measuring program 42 detailed hereinafter, a data file memory 50 storing various data files, a data memory 61 storing data indicating obtained internal information of the scattering medium, a work memory 62 temporarily saving working data, an input device 70 with a keyboard 71 and a mouse 72 for accepting entry of data, and an output device 80 with a display 81 and a printer 82 for output of obtained data, and these are also controlled by the CPU 30 electrically connected thereto. The above memories may be configured in the form of an internal memory (hard disk) of a computer or a flexible disk.

The data file memory 50 stores various data including the waveform data, mean pathlength, instrumental function (impulse response of the measuring system), variance, ratio of transport scattering coefficients, absorption coefficient difference, etc. obtained by execution of the internal information measuring program 42, and also stores the data of measurement conditions, known values, etc. preliminarily entered by use of the input device 70. Such input data includes the shape of the measured medium, the light injection position, the photodetection position, the source-detector distance, the wavelength of light used in the measurement, the type of measurement (e.g., reflection type or transmission type), the extinction coefficient at a predetermined wavelength of the absorber as a measured object, and so on.

The photodetector 7 can be selected from photodetectors of various kinds, including photodiodes, avalanche photodiodes, PIN photodiodes, etc., as well as the photomultiplier tubes. For selection of the photodetector 7 to be used in the measurement, the necessary condition is that it has spectral sensitivity characteristics capable of detecting re-emission of the wavelengths of the measurement light used. Further, it is preferable to use the photodetector with high sensitivity or with high gain if the photon signal is weak. The source light guide 3 and the detector light guide 6 described above may also be replaced by optical fibers, lenses, or the like.

Next, the measuring method according to the present invention will be described below in detail on the basis of the flowchart of an embodiment thereof illustrated in Fig. 4 (which is a flowchart indicating processing of the internal information measuring program 42 illustrated in Fig. 3).

In the flowchart illustrated in Fig. 4, the first step is a step of injecting the pulsed light of a predetermined wavelength generated in the light source 5, through the light guide 3 into the scattering medium 1 at the light injection position uj (S110) and the next step is a step of detecting the light having propagated with being scattered inside the scattering medium 1, by the photodetector 7 through the light guide 6 located at the photodetection position vk (S120).

Then the photodetector 7 outputs the photodetection signal corresponding to the detected light and the signal processing unit 8 converts it into the waveform data indicating the temporal change of intensity of detected light (S130). The instrumental function (the impulse response of the measuring system) is preliminarily measured (S190) and is stored in the data file memory 50. On this occasion, the instrumental function is measured in such a way that the scattering medium 1 is removed from the structure illustrated in Fig. 3 and the light output end of the light guide 3 and the light input end of the light guide 6 are directly opposed and coupled in contact to each other. Therefore, the instrumental function includes the influence of the pulse width of the light source, the bandwidth of the detecting system, and so on.

Then the mean pathlength L and the variance σ2 of plural photons constituting the impulse response are computed, based on the obtained waveform data and instrumental function (S140). The mean pathlength and variance of the impulse response are the sum of the mean pathlengths and the sum of the variances of the measured profile and the instrumental function, respectively, as represented by aforementioned Eqs. (2.8) and (2.10). Since the mean pathlength is expressed by the weighted mean of time-resolved profile as indicated in aforementioned Eq. (1.2b) and the variance by aforementioned (1.3), they can be quickly obtained by computation based on the time-resolved profile, i.e., the waveform data obtained as above with the computer (computation of moment).

Then the square root of the ratio of scattering coefficients, i.e., k' is computed, based on the mean pathlength L and variance σ2, in accordance with Eq. (5.6) (S150). Then the absorption coefficient difference of the scattering medium or the concentration of the absorber is calculated based on aforementioned Eq. (4.8) or based on Eq. (4.9), respectively (S160 or S170), and the result of the calculation is outputted (S180).

The above operation to compute the mean pathlength and variance (S140) may also be arranged to obtain the mean pathlength and variance from the impulse response obtained by the deconvolution operation of the waveform data with the instrumental function. In the above operation to compute the square root of the scattering coefficient ratio (S150), it is also possible to use a value preliminarily measured by another method, as described previously.

In the above embodiment, when the light of the predetermined wavelength is pulsed light of n + 2 (≥ 3, where n is an integer not less than 1) or more wavelengths, it is possible to determine n + 1 absorption coefficient differences and quantify concentrations of n + 1 absorbers from these values.

(Second Embodiment)

The present embodiment is an example of application of the present invention to the phase modulation spectroscopy. In this case, the setup of the measuring system is one obtained by replacing the signal processing unit 8 described above and illustrated in Fig. 3, for example, with an computation unit including a lock-in amplifier. The light source 5 emits modulated light of two predetermined wavelengths λ1 and/or λ2 including three types of modulation frequency components (ω1, ω2, ω3).

Fig. 5 is a flowchart of a method not embodying the present invention applied to the phase modulation spectroscopy. In the flowchart illustrated in Fig. 5, the first step is a step of injecting the intensity-modulated light of the predetermined wavelengths generated in the light source 5, through the light guide 3 into the scattering medium 1 at the light injection position uj (S110) and the next step is a step of detecting the light having propagated with being scattered inside the scattering medium 1, by the photodetector 7 through the light guide 6 located at the photodetection position vk (S120). Then the photodetection signal corresponding to the detected light is outputted from the photodetector 7 to be supplied to the signal processing unit 8.

The lock-in amplifier included in the signal processing unit 8 extracts signals of the three predetermined frequency components from the modulated light of the wavelengths λ1 and λ2 described above (S131) and outputs the real parts R, imaginary parts X, amplitudes A, and phase delays φ stated in Eq. (3.1), about the signals of the three predetermined frequency components. The zero point of the phase is preliminarily acquired (S191). The next step in this embodiment is a step of computing the group delay of plural photons constituting the detected light of the modulation frequency of ω2, and the second partial derivative of logarithm of amplitude with respect to ω (which is equal to the partial derivative of the group delay with respect to the absorption coefficient, see Eq. (3.5)), using the amplitudes A, phase delays φ, and three modulation frequencies (ω1, ω2, ω3) of the signals of the three predetermined frequency components against injection of the modulated light of the two predetermined wavelengths (S141).

Then the square root of the ratio of scattering coefficients, i.e., k' is computed based on this group delay and the partial derivative thereof with respect to the absorption coefficient (S150). In this case, the computation is carried out by replacing the mean pathlength L with c times the group delay and replacing the variance σ2 with the second partial derivative of logarithm of c2 times the amplitude with respect to ω, in aforementioned Eq. (5.6). Then the absorption coefficient difference of the scattering medium, or the concentration of the absorber is computed based on aforementioned Eq. (4.8) or Eq. (4.9), respectively, (S160 or S170) and the result of the computation is outputted (S180). In these cases, however, the computation is also conducted by replacing the mean pathlength L with c times the group delay and replacing the variance σ2 with the second partial derivative of logarithm of c2 times the amplitude with respect to ω.

The above process is arranged to obtain the group delay and the partial derivative thereof with respect to the absorption coefficient, but in the case of ω « cµa, the process may also be modified to compute the phase delay and the partial derivative of the phase delay with respect to the absorption coefficient, because the group delay is approximated to the phase delay in that case as described above. In the above operation to compute the square root of the ratio of scattering coefficients (S150), it is also possible to use a value preliminarily measured by another method.

In the above modification, when the light of the predetermined wavelengths is pulsed light of n + 2 (≥ 3, where n is an integer not less than 1) or more wavelengths, it becomes feasible to obtain n + 1 absorption coefficient differences and quantify concentrations of n + 1 absorbers from these values.

(Third Embodiment)

Fig. 6 shows another embodiment of the present invention to illustrate a system for measuring or monitoring the concentration of hemoglobin or oxygen saturation of hemoglobin inside a scattering medium such as a human head. This system uses light of three wavelengths, i.e., wavelengths λ1, λ2, and λ3. In this case, the operation principle and structure are the same as in the first embodiment with use of pulsed light and the same as in the above modification with use of modulated light. The system of the present embodiment, however, is different in the structure of a container housing the light injecting means and light detecting means from each of the above embodiments.

In the system illustrated in Fig. 6, the light injecting means and light detecting means are housed in the container 10 provided with a fitting band to be mounted on the head 1a like a headband and are connected through a cable 12 to an external device 11 including the signal processing unit 8, the computation processing unit 9, and so on.

Fig. 7 shows the details of the container 10. The container 10 incorporates the light source 5, wavelength selector 4, source light guide 3, detector light guide 6, and photodetector 7. The light of the predetermined wavelengths λ1, λ2, λ3 emitted from the light source 5 is subjected to the wavelength selection in the wavelength selector 4 and light of the selected wavelength is guided through the light guide 3 into the head la. On this occasion, the predetermined wavelengths λ1, λ2, and λ3 are properly selected with reference to the absorption spectrum of hemoglobin illustrated in Fig. 8.

The above system is constructed in such structure that the cable 12 connects the container 10 incorporating the light injecting means and light detecting means, via a connector 13 with the external device 11 housing the signal processing unit 8 and the computation processing unit 9, but the system may also be constructed so as to connect them with each other by another means such as wireless communication, optical signals, or the like. Such arrangement permits, for example, not only measurement on the bed side or in a state of repose, but also measurement in a state of motion. Besides the head, it is also possible to implement measurement with an object such as the thigh or the like of a marathon runner during a race, for example. Further, if the system is connected to a city cable, an optical cable, or the like, it will also become feasible to implement remote measurement for a person at home from facilities such as hospitals or the like and to apply it to centralized control or the like of patient rooms in the hospitals, for example.

Here let us denote the concentrations and optical constants of hemoglobin as follows. <ul id="ul0002" list-style="none" compact="compact"><li>Cb: molar concentration (M) of deoxygenated hemoglobin,</li><li>Co: molar concentration (M) of oxygenated hemoglobin,</li><li>p1: molar extinction coefficient (mm-1M-1) at the wavelength λ1 of deoxygenated hemoglobin,</li><li>p2: molar extinction coefficient (mm-1M-1) at the wavelength λ2 of deoxygenated hemoglobin,</li><li>p3: molar extinction coefficient (mm-1M-1) at the wavelength λ3 of deoxygenated hemoglobin,</li><li>q1: molar extinction coefficient (mm-1M-1) at the wavelength λ4 of oxygenated hemoglobin,</li><li>q2: molar extinction coefficient (mm-1M-1) at the wavelength λ5 of oxygenated hemoglobin,</li><li>q3: molar extinction coefficient (mm-1M-1) at the wavelength λ6 of oxygenated hemoglobin</li></ul> Then the absorption coefficients at the respective wavelengths are expressed as follows.<maths id="math0059" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext> = </mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">b</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><mtext> = </mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">b</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><msub><mrow><mtext> µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>3</mtext></mrow></msub><mtext> = </mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>3</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">b</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>3</mtext></mrow></msub><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>3 (</mtext><mtext mathvariant="italic">A</mtext><mtext>1.1)</mtext></mrow></msub></mrow></math><img file="EP1136811B1_D0059.tif" /></maths> In these equations, a1, a2, and a3 are absorption coefficients of water and the absorbers other than hemoglobin.

Since values actually measured by the method of the present invention are absorption coefficient differences, the above equations are rewritten as follows corresponding thereto.<maths id="math0060" num=""><math display="block"><mrow><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>2</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext> = (</mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">b</mtext></mrow></msub><mtext> + (</mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><msub><mrow><mtext> µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>3</mtext></mrow></msub><msub><mrow><mtext> - µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext><mtext>1</mtext></mrow></msub><mtext> = (</mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>3</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">p</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">b</mtext></mrow></msub><mtext> + (</mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>3</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">q</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>)</mtext><msub><mrow><mtext mathvariant="italic">C</mtext></mrow><mrow><mtext mathvariant="italic">o</mtext></mrow></msub><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>3</mtext></mrow></msub><mtext> - </mtext><msub><mrow><mtext mathvariant="italic">a</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext> (</mtext><mtext mathvariant="italic">A</mtext><mtext>1.2)</mtext></mrow></math><img file="EP1136811B1_D0060.tif" /></maths> Namely, the left sides of the above equations are quantities measured by the method of the present invention. Since these equations are simultaneous equations including two unknowns Cb and Co, the two unknowns Cb and Co can be determined by solving these simultaneous equations. In this case, as to values of a1, a2, and a3, the wavelengths are normally selected so as to make a1 ≒ a2 ≒ a3. It is also possible to use standard values of living tissues for the values of a1, a2, a3.

According to the above, this embodiment allows us to determine the concentration Cb of deoxygenated hemoglobin, the concentration Co of oxygenated hemoglobin, the amount of hemoglobin (Cb + Co), and the oxygen saturation Co/(Cb + Co).

The above described the preferred embodiments of the present invention, but it is noted that the present invention is by no means limited to the above embodiments, of course.

Specifically, the light injection position and photodetection position are fixed in the above embodiments, but the light injection position and/or the photodetection position may be moved to scan. It can also be contemplated that a plurality of light injection positions and/or photodetection positions are arranged around the scattering medium.

Further, in the above embodiments the scattering characteristics during measurement can be considered to be constant when the measurement is carried out with moving the light injection position and/or the photodetection position relative to one scattering medium, or when the measurement is carried out at a fixed position of a scattering medium before and after change in the concentration of an absorber. In this case, a concentration distribution of an absorber (differences with respect to a reference value) or a concentration change of an absorber (differences with respect to a reference value) can be measured by setting k' = 1 in aforementioned Eq. (4.9) and Eq. (4.10). Particularly, the former example is effective, for example, to the measurement in Example 2 (mammography) disclosed in Japanese Patent Application Laid-Open No. H10-73481 by the inventors of the present application.

The present invention will be described in further detail on the basis of examples.

(Example 1)

The present example provides the results of simulation performed for confirming the accuracy of the method of the present invention.

Specifically, Monte Carlo simulation was conducted as to the spectroscopy using the mean pathlength and variance according to the present invention (the MVS method) with variation of absorption coefficient for slablike scattering media (media) having the thickness of 30 mm. Although this Monte Carlo simulation does not present data directly corresponding to data measured by dual-wavelength spectroscopy, it allows simulation of validity of the dual-wavelength spectroscopy, using Monte Carlo data at different concentrations of an absorber. Namely, aforementioned Eq. (4.8) with k' = 1 is applied to quantification of the difference between absorption coefficients.

Table 1 and Table 2 present the results of measurement with the source-detector distance of 5 mm in the reflection measurement and Table 3 and Table 4 the results of measurement with the source-detector distance of 30 mm in the reflection measurement. Table 5 and Table 6 present the results of measurement with the source-detector distance of 30 mm (corresponding to the thickness of the medium) in the transmission measurement. The scattering coefficient µ's = 1mm-1 in all the cases, but the mean cosine of scattering angles is g = 0.6 in Tables 1, 3, and 5 and g = 0.9 in the Tables 2, 4, and 6. <tables id="tabl0001" num="0001"><table frame="all"><title>TABLE 1</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error (%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error(%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">26.063</entry><entry namest="col3" nameend="col3" align="char" char=".">31.90</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">22.98</entry><entry namest="col3" nameend="col3" align="char" char=".">23.15</entry><entry namest="col4" nameend="col4" align="center">0.0047</entry><entry namest="col5" nameend="col5" align="center">-6</entry><entry namest="col6" nameend="col6" align="center">0.0053</entry><entry namest="col7" nameend="col7" align="center">6</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">19.30</entry><entry namest="col3" nameend="col3" align="char" char=".">16.17</entry><entry namest="col4" nameend="col4" align="center">0.0092</entry><entry namest="col5" nameend="col5" align="center">-8</entry><entry namest="col6" nameend="col6" align="center">0.0111</entry><entry namest="col7" nameend="col7" align="center">11</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">17.23</entry><entry namest="col3" nameend="col3" align="char" char=".">12.91</entry><entry namest="col4" nameend="col4" align="center">0.0097</entry><entry namest="col5" nameend="col5" align="center">-3</entry><entry namest="col6" nameend="col6" align="center">0.0114</entry><entry namest="col7" nameend="col7" align="center">14</entry></row></tbody></tgroup></table></tables><tables id="tabl0002" num="0002"><table frame="all"><title>TABLE 2</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error (%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">26.19</entry><entry namest="col3" nameend="col3" align="char" char=".">31.58</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">22.62</entry><entry namest="col3" nameend="col3" align="char" char=".">22.91</entry><entry namest="col4" nameend="col4" align="center">0.0047</entry><entry namest="col5" nameend="col5" align="center">-6</entry><entry namest="col6" nameend="col6" align="center">0.0053</entry><entry namest="col7" nameend="col7" align="center">6</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">19.00</entry><entry namest="col3" nameend="col3" align="char" char=".">15.99</entry><entry namest="col4" nameend="col4" align="center">0.0093</entry><entry namest="col5" nameend="col5" align="center">-7</entry><entry namest="col6" nameend="col6" align="center">0.0111</entry><entry namest="col7" nameend="col7" align="center">11</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">16.98</entry><entry namest="col3" nameend="col3" align="char" char=".">12.76</entry><entry namest="col4" nameend="col4" align="center">0.0097</entry><entry namest="col5" nameend="col5" align="center">-3</entry><entry namest="col6" nameend="col6" align="center">0.0116</entry><entry namest="col7" nameend="col7" align="center">16</entry></row></tbody></tgroup></table></tables><tables id="tabl0003" num="0003"><table frame="all"><title>TABLE 3</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error (%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">283.50</entry><entry namest="col3" nameend="col3" align="char" char=".">137.56</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">219.51</entry><entry namest="col3" nameend="col3" align="char" char=".">93.43</entry><entry namest="col4" nameend="col4" align="center">0.0046</entry><entry namest="col5" nameend="col5" align="center">-8</entry><entry namest="col6" nameend="col6" align="center">0.0046</entry><entry namest="col7" nameend="col7" align="center">-8</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">165.16</entry><entry namest="col3" nameend="col3" align="char" char=".">59.20</entry><entry namest="col4" nameend="col4" align="center">0.0089</entry><entry namest="col5" nameend="col5" align="center">-11</entry><entry namest="col6" nameend="col6" align="center">0.0102</entry><entry namest="col7" nameend="col7" align="center">2</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">138.77</entry><entry namest="col3" nameend="col3" align="char" char=".">44.44</entry><entry namest="col4" nameend="col4" align="center">0.0096</entry><entry namest="col5" nameend="col5" align="center">-4</entry><entry namest="col6" nameend="col6" align="center">0.0110</entry><entry namest="col7" nameend="col7" align="center">10</entry></row></tbody></tgroup></table></tables><tables id="tabl0004" num="0004"><table frame="all"><title>TABLE 4</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error (%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">284.02</entry><entry namest="col3" nameend="col3" align="char" char=".">137.65</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">219.75</entry><entry namest="col3" nameend="col3" align="char" char=".">93.32</entry><entry namest="col4" nameend="col4" align="center">0.0046</entry><entry namest="col5" nameend="col5" align="center">- 8</entry><entry namest="col6" nameend="col6" align="center">0.0046</entry><entry namest="col7" nameend="col7" align="center">-8</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">165.44</entry><entry namest="col3" nameend="col3" align="char" char=".">59.03</entry><entry namest="col4" nameend="col4" align="center">0.0089</entry><entry namest="col5" nameend="col5" align="center">-11</entry><entry namest="col6" nameend="col6" align="center">0.0104</entry><entry namest="col7" nameend="col7" align="center">4</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">139.39</entry><entry namest="col3" nameend="col3" align="char" char=".">44.40</entry><entry namest="col4" nameend="col4" align="center">0.0095</entry><entry namest="col5" nameend="col5" align="center">- 5</entry><entry namest="col6" nameend="col6" align="center">0.0110</entry><entry namest="col7" nameend="col7" align="center">10</entry></row></tbody></tgroup></table></tables><tables id="tabl0005" num="0005"><table frame="all"><title>TABLE 5</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error (%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">231.94</entry><entry namest="col3" nameend="col3" align="char" char=".">110.03</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">190.39</entry><entry namest="col3" nameend="col3" align="char" char=".">78.01</entry><entry namest="col4" nameend="col4" align="center">0.0046</entry><entry namest="col5" nameend="col5" align="center">- 8</entry><entry namest="col6" nameend="col6" align="center">0.0056</entry><entry namest="col7" nameend="col7" align="center">12</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">150.92</entry><entry namest="col3" nameend="col3" align="char" char=".">51.83</entry><entry namest="col4" nameend="col4" align="center">0.0090</entry><entry namest="col5" nameend="col5" align="center">-10</entry><entry namest="col6" nameend="col6" align="center">0.0117</entry><entry namest="col7" nameend="col7" align="center">17</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">129.03</entry><entry namest="col3" nameend="col3" align="char" char=".">40.63</entry><entry namest="col4" nameend="col4" align="center">0.0101</entry><entry namest="col5" nameend="col5" align="center">1</entry><entry namest="col6" nameend="col6" align="center">0.0104</entry><entry namest="col7" nameend="col7" align="center">4</entry></row></tbody></tgroup></table></tables><tables id="tabl0006" num="0006"><table frame="all"><title>TABLE 6</title><tgroup cols="7" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="22.50mm" /><colspec colnum="2" colname="col2" colwidth="22.50mm" /><colspec colnum="3" colname="col3" colwidth="22.50mm" /><colspec colnum="4" colname="col4" colwidth="22.50mm" /><colspec colnum="5" colname="col5" colwidth="22.50mm" /><colspec colnum="6" colname="col6" colwidth="22.50mm" /><colspec colnum="7" colname="col7" colwidth="22.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">µa</entry><entry namest="col2" nameend="col2" align="center">L</entry><entry namest="col3" nameend="col3" align="center">σ</entry><entry namest="col4" nameend="col4" align="center">Δµa</entry><entry namest="col5" nameend="col5" align="center">Error(%)</entry><entry namest="col6" nameend="col6" align="center">Old-Δµa</entry><entry namest="col7" nameend="col7" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.005</entry><entry namest="col2" nameend="col2" align="char" char=".">232.00</entry><entry namest="col3" nameend="col3" align="char" char=".">110.23</entry><entry namest="col4" nameend="col4" align="center">-----</entry><entry namest="col5" nameend="col5" align="center">-----</entry><entry namest="col6" nameend="col6" align="center">-----</entry><entry namest="col7" nameend="col7" align="center">-----</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.01</entry><entry namest="col2" nameend="col2" align="char" char=".">190.29</entry><entry namest="col3" nameend="col3" align="char" char=".">77.74</entry><entry namest="col4" nameend="col4" align="center">0.0044</entry><entry namest="col5" nameend="col5" align="center">-12</entry><entry namest="col6" nameend="col6" align="center">0.0057</entry><entry namest="col7" nameend="col7" align="center">14</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.02</entry><entry namest="col2" nameend="col2" align="char" char=".">150.39</entry><entry namest="col3" nameend="col3" align="char" char=".">52.05</entry><entry namest="col4" nameend="col4" align="center">0.0091</entry><entry namest="col5" nameend="col5" align="center">-9</entry><entry namest="col6" nameend="col6" align="center">0.0113</entry><entry namest="col7" nameend="col7" align="center">13</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.03</entry><entry namest="col2" nameend="col2" align="char" char=".">129.18</entry><entry namest="col3" nameend="col3" align="char" char=".">40.22</entry><entry namest="col4" nameend="col4" align="center">0.0098</entry><entry namest="col5" nameend="col5" align="center">-2</entry><entry namest="col6" nameend="col6" align="center">0.0116</entry><entry namest="col7" nameend="col7" align="center">16</entry></row></tbody></tgroup></table></tables>

In each of the tables µa in the left extreme column represents the absorption coefficients of media set on the occasion of Monte Carlo calculation. In the tables µa in the central part represents the absorption coefficient differences quantified by aforementioned Eq. (4.8), using the mean pathlengths L and variances σ2 computed from a set of Monte Carlo data of upper and lower neighbors. However, k' = 1. In the tables the right part presents Old-Δµa to indicate the absorption coefficient differences obtained by another method for quantifying the ratio of transport scattering coefficients (Eq. (5.3a)) for comparison. From the tables, this Old-Δµa shows a tendency of DC biasing in the reflection measurement and a tendency of the slope becoming larger than 1 in the transmission measurement, and Δµa obtained by the method according to the present invention demonstrates smaller errors than Old-Δµa.

The smaller the absorption coefficient µa, the larger the change of slope of mean pathlength L against absorption coefficient µa becomes (see Fig. 1). Therefore, there arises the error of linear approximation (equivalent to use of the mean value theorem) used for deriving aforementioned Eq. (1.5), in the small region of the absorption coefficient µa. Of course, this error decreases with decrease in the absorption coefficient difference. This tendency is seen in the results of Table 1 to Table 6 and the quantified values in the bottom row show considerably good accuracy. In the spectroscopy of actual living tissues or the like, the difference between absorption coefficients at two wavelengths becomes small just when the absorption coefficient µa is small, as described hereinafter. Since the living tissues involve the absorption of water, absolute values of absorption coefficients also become larger by that degree. Accordingly, these two points largely relax the above error in the spectroscopy of living tissues or the like.

(Example 2)

The present example provides the results of experiments conducted with a simulated phantom in order to confirm the accuracy of the method of the present invention.

The schematic structure of the experimental setup was as illustrated in Fig. 3, and the light source was two picosecond pulse generators that generated picosecond pulses of the respective wavelengths λ1 = 782 nm and λ2 = 831 nm, the repetition rate of 5 MHz, and the pulse width of about 50 ps. These picosecond pulses are injected through an optical switch and an optical attenuator into a GI fiber having the diameter of 200 µm and light emerging from an output end being the other end of the fiber is injected into the phantom being a scattering medium. Output light from the phantom is received by a bundle fiber having the diameter of 5 mm and is measured by the system based on the time-correlated single photon counting method as illustrated in Fig. 3. The instrumental functions at the two wavelengths necessary for the computation of mean pathlength are measured with the light injection and photodetection fibers kept in a closely fitted state.

The phantom used in the experiments was prepared by putting 420 ml of a 1% intralipid solution as a scattering substance into an acrylic vessel (width 120 mm, height 120 mm, and depth 40 mm) and gradually adding greenish brown ink 0.07 ml every time as an absorbing substance to the solution until the total amount of the ink became 0.56 ml, and the transmission measurement was carried out under the condition of the source-detector distance of 40 mm. The actual absorption coefficient of the phantom is the sum of the absorption coefficient of the added ink and the absorption coefficient of water (distilled water). At the wavelengths λ1 = 782 nm and λ2 = 831 nm of the pulsed light used in the measurement, theoretical values of the transport scattering coefficient are 1.0207 and 0.9531, respectively. Therefore, a ratio of them is 1.071 and k' = 1/1.035 = 0.9633. The extinction coefficient of greenish brown ink and the absorption coefficient of distilled water were measured by a spectroscope. Table 7 presents the optical parameters of the phantom in the experiments. <tables id="tabl0007" num="0007"><table frame="all"><title>TABLE 7</title><tgroup cols="4" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="39.37mm" /><colspec colnum="2" colname="col2" colwidth="39.37mm" /><colspec colnum="3" colname="col3" colwidth="39.37mm" /><colspec colnum="4" colname="col4" colwidth="39.37mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">Wavelength (nm)</entry><entry namest="col2" nameend="col2" align="center">λ1 = 782</entry><entry namest="col3" nameend="col3" align="center">λ2 = 831</entry><entry namest="col4" nameend="col4" align="center">Difference between or ratio of coefficients at the wavelengths</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="center">Absorption coefficient of water (mm-1)</entry><entry namest="col2" nameend="col2" align="center">µw1 = 2.452×10-3</entry><entry namest="col3" nameend="col3" align="center">µw2 = 3.095×10-3</entry><entry namest="col4" nameend="col4" align="center">µw2 -µw1 = 0.643×10-3</entry></row><row><entry namest="col1" nameend="col1" align="center">Extinction coefficient (mm-1·ml-1)</entry><entry namest="col2" nameend="col2" align="center">ε1 = 4.806×10-2</entry><entry namest="col3" nameend="col3" align="center">ε2 =2.758×10-2</entry><entry namest="col4" nameend="col4" align="center">ε2-ε1 = -2.048×10-2</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="center">Transport scattering coefficient (mm-1)</entry><entry namest="col2" nameend="col2" align="center">µs1 = 1.0207</entry><entry namest="col3" nameend="col3" align="center">µs2=0.9531</entry><entry namest="col4" nameend="col4" align="center">(µs2/µs1)1/2=0.9633</entry></row></tbody></tgroup></table></tables> Table 8 presents values of the mean pathlength and variance of the impulse response calculated from experimental values obtained in the experiments with the phantom carried out under the above conditions. <tables id="tabl0008" num="0008"><table frame="all"><title>TABLE 8</title><tgroup cols="5" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="31.50mm" /><colspec colnum="2" colname="col2" colwidth="31.50mm" /><colspec colnum="3" colname="col3" colwidth="31.50mm" /><colspec colnum="4" colname="col4" colwidth="31.50mm" /><colspec colnum="5" colname="col5" colwidth="31.50mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" align="center">Ink (ml)</entry><entry namest="col2" nameend="col2" align="center">L(λ1) (×102 mm)</entry><entry namest="col3" nameend="col3" align="center">L(λ2) (×102 mm)</entry><entry namest="col4" nameend="col4" align="center">σ2(λ1) (×104 mm)</entry><entry namest="col5" nameend="col5" align="center">σ2(λ2) (×104 mm)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.00</entry><entry namest="col2" nameend="col2" align="char" char=".">4.2007</entry><entry namest="col3" nameend="col3" align="char" char=".">3.8468</entry><entry namest="col4" nameend="col4" align="char" char=".">3.5658</entry><entry namest="col5" nameend="col5" align="char" char=".">2.9339</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.07</entry><entry namest="col2" nameend="col2" align="char" char=".">3.3142</entry><entry namest="col3" nameend="col3" align="char" char=".">3.3683</entry><entry namest="col4" nameend="col4" align="char" char=".">1.6680</entry><entry namest="col5" nameend="col5" align="char" char=".">1.8673</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.14</entry><entry namest="col2" nameend="col2" align="char" char=".">2.8670</entry><entry namest="col3" nameend="col3" align="char" char=".">3.0533</entry><entry namest="col4" nameend="col4" align="char" char=".">1.0439</entry><entry namest="col5" nameend="col5" align="char" char=".">1.3639</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.21</entry><entry namest="col2" nameend="col2" align="char" char=".">2.5690</entry><entry namest="col3" nameend="col3" align="char" char=".">2.8141</entry><entry namest="col4" nameend="col4" align="char" char=".">0.7275</entry><entry namest="col5" nameend="col5" align="char" char=".">1.0329</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.28</entry><entry namest="col2" nameend="col2" align="char" char=".">2.3322</entry><entry namest="col3" nameend="col3" align="char" char=".">2.6051</entry><entry namest="col4" nameend="col4" align="char" char=".">0.5572</entry><entry namest="col5" nameend="col5" align="char" char=".">0.8342</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.35</entry><entry namest="col2" nameend="col2" align="char" char=".">2.1610</entry><entry namest="col3" nameend="col3" align="char" char=".">2.4412</entry><entry namest="col4" nameend="col4" align="char" char=".">0.4399</entry><entry namest="col5" nameend="col5" align="char" char=".">0.6946</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.42</entry><entry namest="col2" nameend="col2" align="char" char=".">2.0367</entry><entry namest="col3" nameend="col3" align="char" char=".">2.3369</entry><entry namest="col4" nameend="col4" align="char" char=".">0.3605</entry><entry namest="col5" nameend="col5" align="char" char=".">0.5910</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.49</entry><entry namest="col2" nameend="col2" align="char" char=".">1.9164</entry><entry namest="col3" nameend="col3" align="char" char=".">2.2304</entry><entry namest="col4" nameend="col4" align="char" char=".">0.2993</entry><entry namest="col5" nameend="col5" align="char" char=".">0.5125</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.56</entry><entry namest="col2" nameend="col2" align="char" char=".">1.8336</entry><entry namest="col3" nameend="col3" align="char" char=".">2.1266</entry><entry namest="col4" nameend="col4" align="char" char=".">0.2611</entry><entry namest="col5" nameend="col5" align="char" char=".">0.4546</entry></row></tbody></tgroup></table></tables>

The absorption of water is unignorable in the measurement with the living tissues and living-tissue-simulated liquid phantoms. The extinction coefficient of an absorber differs depending upon wavelengths, and the spectroscopy makes use of this difference. Fig. 9 shows the relation between absorption coefficient and mean pathlength of a phantom. In this relation the absorption coefficients of the phantom were calculated using the absorption coefficient of distilled water actually measured, the quantity of added ink, and the extinction coefficient of ink actually measured. The curves in the figure are the results of fitting in the general form of the following equation indicating the mean photon pathlength obtained by the photon diffusion approximation.<maths id="math0061" num=""><math display="block"><mrow><mtext mathvariant="italic">L</mtext><mtext>=1/</mtext><mfenced open="(" close=")"><mrow><mtext mathvariant="italic">a</mtext><mtext>+</mtext><mtext mathvariant="italic">b</mtext><msqrt><msub><mrow><mtext>µ</mtext></mrow><mrow><mtext mathvariant="italic">a</mtext></mrow></msub></msqrt></mrow></mfenced></mrow></math><img file="EP1136811B1_D0061.tif" /></maths> It is seen from this Fig. 9 that Eq. (4.6) approximately holds.

Further, it is also seen that, considering the absorption coefficient difference in a set of data used in the spectroscopy in the figure, i.e., in a set of data at the same quantity of added ink, this difference increases with increase in the absorption coefficient of the phantom and the difference also becomes smaller with decrease in the absorption coefficient of the phantom. This relation generally holds for the spectroscopy of living tissues or the like. The linear approximation was used for deriving Eq. (1.5) being the fundamental equation of the MVS method. Since the curvature of mean pathlength and the absorption coefficient difference are just in the relation of one increasing with decrease of the other, as illustrated in Fig. 9, it is seen that the linear approximation holds with small error throughout the entire region of change of the absorption coefficient.

Concentrations C of the absorber quantified by the experiments are presented in Table 9. For comparison's sake, Table 9 also presents the results by the method using a fixed value for the ratio k' (k' = 0.9663, which is a theoretical value), the results on the assumption that the transport scattering coefficient has no wavelength dependence (k' = 1), and the results by use of an average of the ratios of transport scattering coefficients obtained by the experiments (k' = 0.9681). It is apparent from this table that the errors are naturally large on the assumption that the transport scattering coefficient has no wavelength dependence (k' = 1) and that the method according to the present invention assures excellent accuracy as against the other methods.

As compared with the results by use of the average of the ratios gained by the experiments (k' = 0.9681), the results with k' obtained in each measurement show smaller errors, and this is conceivably because the influence of noise in each measurement is canceled out in the process of operation. <tables id="tabl0009" num="0009"><table frame="all"><title>TABLE 9</title><tgroup cols="10" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="15.75mm" /><colspec colnum="2" colname="col2" colwidth="15.75mm" /><colspec colnum="3" colname="col3" colwidth="15.75mm" /><colspec colnum="4" colname="col4" colwidth="15.75mm" /><colspec colnum="5" colname="col5" colwidth="15.75mm" /><colspec colnum="6" colname="col6" colwidth="15.75mm" /><colspec colnum="7" colname="col7" colwidth="15.75mm" /><colspec colnum="8" colname="col8" colwidth="15.75mm" /><colspec colnum="9" colname="col9" colwidth="15.75mm" /><colspec colnum="10" colname="col10" colwidth="15.75mm" /><thead valign="top"><row><entry namest="col1" nameend="col1" /><entry namest="col2" nameend="col4" align="center">Method of the invention</entry><entry namest="col5" nameend="col6" align="center">k' = 0.9663</entry><entry namest="col7" nameend="col8" align="center">k' = 1</entry><entry namest="col9" nameend="col10" align="center">k' = 0.9681</entry></row><row><entry namest="col1" nameend="col1" align="left">Ink (ml)</entry><entry namest="col2" nameend="col2" align="center">k'</entry><entry namest="col3" nameend="col3" align="center">C (ml)</entry><entry namest="col4" nameend="col4" align="center">Error (%)</entry><entry namest="col5" nameend="col5" align="center">Ck'2 (ml)</entry><entry namest="col6" nameend="col6" align="center">Error (%)</entry><entry namest="col7" nameend="col7" align="center">Ck'1 (ml)</entry><entry namest="col8" nameend="col8" align="center">Error (%)</entry><entry namest="col9" nameend="col9" align="center">Ck'3 (ml)</entry><entry namest="col10" nameend="col10" align="center">Error (%)</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="char" char=".">0.00</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9646</entry><entry namest="col3" nameend="col3" align="center">-0.000</entry><entry namest="col4" nameend="col4" align="center">--</entry><entry namest="col5" nameend="col5" align="center">-0,001</entry><entry namest="col6" nameend="col6" align="center">--</entry><entry namest="col7" nameend="col7" align="center">-0.022</entry><entry namest="col8" nameend="col8" align="center">--</entry><entry namest="col9" nameend="col9" align="center">-0.002</entry><entry namest="col10" nameend="col10" align="center">--</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.07</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9737</entry><entry namest="col3" nameend="col3" align="center">0.071</entry><entry namest="col4" nameend="col4" align="center">1.4</entry><entry namest="col5" nameend="col5" align="center">0.078</entry><entry namest="col6" nameend="col6" align="center">11.4</entry><entry namest="col7" nameend="col7" align="center">0.046</entry><entry namest="col8" nameend="col8" align="center">-34</entry><entry namest="col9" nameend="col9" align="center">0.076</entry><entry namest="col10" nameend="col10" align="center">8.5</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.14</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9697</entry><entry namest="col3" nameend="col3" align="center">0.144</entry><entry namest="col4" nameend="col4" align="center">2.9</entry><entry namest="col5" nameend="col5" align="center">0.148</entry><entry namest="col6" nameend="col6" align="center">5.7</entry><entry namest="col7" nameend="col7" align="center">0.107</entry><entry namest="col8" nameend="col8" align="center">-24</entry><entry namest="col9" nameend="col9" align="center">0.146</entry><entry namest="col10" nameend="col10" align="center">4.3</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.21</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9709</entry><entry namest="col3" nameend="col3" align="center">0.211</entry><entry namest="col4" nameend="col4" align="center">0.5</entry><entry namest="col5" nameend="col5" align="center">0.218</entry><entry namest="col6" nameend="col6" align="center">3.8</entry><entry namest="col7" nameend="col7" align="center">0.168</entry><entry namest="col8" nameend="col8" align="center">-20</entry><entry namest="col9" nameend="col9" align="center">0.216</entry><entry namest="col10" nameend="col10" align="center">2.9</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.28</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9738</entry><entry namest="col3" nameend="col3" align="center">0.269</entry><entry namest="col4" nameend="col4" align="center">-3.9</entry><entry namest="col5" nameend="col5" align="center">0.282</entry><entry namest="col6" nameend="col6" align="center">0.7</entry><entry namest="col7" nameend="col7" align="center">0.223</entry><entry namest="col8" nameend="col8" align="center">-20</entry><entry namest="col9" nameend="col9" align="center">0.279</entry><entry namest="col10" nameend="col10" align="center">-0.4</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.35</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9652</entry><entry namest="col3" nameend="col3" align="center">0.342</entry><entry namest="col4" nameend="col4" align="center">-2.3</entry><entry namest="col5" nameend="col5" align="center">0.340</entry><entry namest="col6" nameend="col6" align="center">-2.9</entry><entry namest="col7" nameend="col7" align="center">0.273</entry><entry namest="col8" nameend="col8" align="center">-22</entry><entry namest="col9" nameend="col9" align="center">0.337</entry><entry namest="col10" nameend="col10" align="center">-3.7</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.42</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9693</entry><entry namest="col3" nameend="col3" align="center">0.409</entry><entry namest="col4" nameend="col4" align="center">-2.6</entry><entry namest="col5" nameend="col5" align="center">0.416</entry><entry namest="col6" nameend="col6" align="center">-1.0</entry><entry namest="col7" nameend="col7" align="center">0.340</entry><entry namest="col8" nameend="col8" align="center">-19</entry><entry namest="col9" nameend="col9" align="center">0.412</entry><entry namest="col10" nameend="col10" align="center">-1.9</entry></row><row><entry namest="col1" nameend="col1" align="char" char=".">0.49</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9690</entry><entry namest="col3" nameend="col3" align="center">0.487</entry><entry namest="col4" nameend="col4" align="center">-0.6</entry><entry namest="col5" nameend="col5" align="center">0.494</entry><entry namest="col6" nameend="col6" align="center">0.8</entry><entry namest="col7" nameend="col7" align="center">0.410</entry><entry namest="col8" nameend="col8" align="center">-16</entry><entry namest="col9" nameend="col9" align="center">0.489</entry><entry namest="col10" nameend="col10" align="center">-0.2</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="char" char=".">0.56</entry><entry namest="col2" nameend="col2" align="char" char=".">0.9566</entry><entry namest="col3" nameend="col3" align="center">0.549</entry><entry namest="col4" nameend="col4" align="center">-2.0</entry><entry namest="col5" nameend="col5" align="center">0.520</entry><entry namest="col6" nameend="col6" align="center">-7.1</entry><entry namest="col7" nameend="col7" align="center">0.432</entry><entry namest="col8" nameend="col8" align="center">-23</entry><entry namest="col9" nameend="col9" align="center">0.518</entry><entry namest="col10" nameend="col10" align="center">-7.5</entry></row></tbody></tgroup></table></tables>

Industrial Applicability

The measuring methods and apparatus of internal information of scattering medium according to the present invention are applicable to measurement of the absorption coefficient of the scattering medium and the concentration of the absorber, as measuring methods and apparatus of internal information of scattering medium capable of measurement with higher accuracy and at higher speed than the conventional measuring methods based on the MBL law. Particularly, the measuring methods and apparatus according to the present invention make use of the mean pathlength and variance, or physical quantities equivalent thereto, but do not make use of information of the absolute value of light intensity or the ratio thereof. Accordingly, they are useful in that a solution is given to the difficult problem in practical use; quantifying or estimating the absolute value of quantity of light injected into the medium.

Such methods make real-time measurement feasible and also possess the great feature that the measurement can be conducted independent of the shape and boundary conditions, the size of the medium, the scattering characteristics, the source-detector distance, the measurement mode of either transmission or reflection, and so on. From the above, the measuring methods and apparatus of internal information of scattering medium according to the present invention are expected to be commonly applied to apparatus for simply measuring various physiological substances in living tissues on a noninvasive and real-time basis.

Every citation, both waysCites: the store holds 8 of 9

Citations
Document	Relation	Office
EP0826958A	Cites	European Patent Office (EPO)
EP0834277A	Cites	European Patent Office (EPO)
WO9834097A	Cites	World Intellectual Property Organization (WIPO)
JP7163571A	Cites	Japan
JP10073481A	Cites	Japan
JP10111238A	Cites	Japan
US5713352A	Cites	United States of America
US5792051A	Cites	United States of America
Japanese Journal of Applied Physics, Vol. 37, No. 5A, Part I, The Japan Society of Applied Physic, 15 May 1998, pages 2724-2727, XP002924572.	Non-patent	–

11 members in 6 offices

This recordMember, by publication datePriority claim

Priority claims9

Priority claims
Document	Office	Kind	Date
31461398	Japan	A	1998-11-05
31461398	Japan	A	1998-11-05
31461398	Japan	–	1998-11-05
9906181	Japan	W	1999-11-05
9906181	Japan	W	1999-11-05
31461398	–	–	–
JP19980314613	–	–	–
JP1999006181	–	–	–
WO1999JP06181	–	–	–

Members11

Family members, oldest first
Document	Office	Kind
WO0028302A1	World Intellectual Property Organization (WIPO)	A1	2000-05-18
JP2000146828A	Japan	A	2000-05-26
AU1078200A	Australia	A	2000-05-29
EP1136811A1	European Patent Office (EPO)	A1	2001-09-26
US2001038454A1	United States of America	A1	2001-11-08
EP1136811A4	European Patent Office (EPO)	A4	2003-01-02
US6704110B2	United States of America	B2	2004-03-09
EP1136811B1This record	European Patent Office (EPO)	B1	2005-06-22
DE69925930D1	Germany	D1	2005-07-28
DE69925930T2	Germany	T2	2006-05-11
JP3950243B2	Japan	B2	2007-07-25

24 legal events, as 3 offices reported them to INPADOC

Over the term

Point at a mark for the event

Events

Legal events, newest first
	Event	Code	Office
2020-01-31	Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25	PG25	EP
2019-11-27	Patent expired after termination of 20 yearsExpiredPE20	PE20	GB
2019-11-05	Expiry of rightR071	R071	DE
2019-02-28	Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP	PGFP	EP
2019-01-31	Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP	PGFP	EP
2006-06-14	No opposition filedOpposition26N	26N	EP
2006-05-12	No opposition filed within time limitOppositionORIGINAL CODE: 0009261PLBE	PLBE	EP
2006-05-12	Information on the status of an ep patent application or granted ep patentGrantedSTATUS: NO OPPOSITION FILED WITHIN TIME LIMITSTAA	STAA	EP
2005-07-28	Corresponds to:REF	REF	EP
2005-06-22	Designated contracting statesAK	AK	EP
2005-06-22	European patent grantedGrantedFG4D	FG4D	GB
2005-05-06	(expected) grantORIGINAL CODE: 0009210GRAA	GRAA	EP
2005-05-04	Grant fee paidORIGINAL CODE: EPIDOSNIGR3GRAS	GRAS	EP
2004-12-02	Despatch of communication of intention to grant a patentORIGINAL CODE: EPIDOSNIGR1GRAP	GRAP	EP
2004-07-21	First examination report despatched17Q	17Q	EP
2004-05-26	Designated contracting states (corrected)RBV	RBV	EP
2003-01-02	Supplementary search report drawn up and despatchedA4	A4	EP
2003-01-02	Designated contracting statesAK	AK	EP
2002-08-07	Information provided on ipc code assigned before grant7G 01N 21/17 A, 7A 61B 5/145 B, 7A 61B 5/00 BRIC1	RIC1	EP
2001-09-26	Request for examination filed17P	17P	EP
2001-09-26	Designated contracting statesAK	AK	EP
2001-09-26	Designated contracting statesAK	AK	EP
2001-09-26	Request for extension of the european patentAL;LT;LV;MK;RO;SIAX	AX	EP
2001-08-10	Public reference made under article 153(3) epc to a published international application that has entered the european phaseORIGINAL CODE: 0009012PUAI	PUAI	EP

Method and device for measuring internal information of scattering absorber

Abstract

Term

Priority claims9

Members11

Over the term

Events

Numbers

Titles3

Classification

Designated states1