FIELD OF THE INVENTION

[0001]
The invention relates generally to pulse oximeters used to detect blood oxygenation. More specifically, the invention relates to a method for taking into account human variability in pulse oximeters. The invention further relates to a sensor allowing compensation for the inaccuracies caused by human variability, the sensor being an integral part of the pulse oximeter.
BACKGROUND OF THE INVENTION

[0002]
Pulse oximetry is at present the standard of care for continuous monitoring of arterial oxygen saturation (SPO_{2}). Pulse oximeters provide instantaneous invivo measurements of arterial oxygenation, and thereby an early warning of arterial hypoxemia, for example.

[0003]
A pulse oximeter comprises a computerized measuring unit and a probe attached to the patient, typically to a finger or ear lobe. The probe includes a light source for sending an optical signal through the tissue and a photo detector for receiving the signal after transmission through the tissue. On the basis of the transmitted and received signals, light absorption by the tissue can be determined. During each cardiac cycle, light absorption by the tissue varies cyclically. During the diastolic phase, absorption is caused by venous blood, tissue, bone, and pigments, whereas during the systolic phase there is an increase in absorption, which is caused by the influx of arterial blood into the tissue. Pulse oximeters focus the measurement on this arterial blood portion by determining the difference between the peak absorption during the systolic phase and the constant absorption during the diastolic phase. Pulse oximetry is thus based on the assumption that the pulsatile component of the absorption is due to arterial blood only.

[0004]
Light transmission through an ideal absorbing sample is determined by the known LambertBeer equation as follows:
I _{out} =I _{in} e ^{−εDC} (1)
where I_{in }is the light intensity entering the sample, I_{out }is the light intensity received from the sample, D is the path length through the sample, ε is the extinction coefficient of the analyte in the sample at a specific wavelength, and C is the concentration of the analyte. When I_{in}, D, and ε are known, and I_{out }is measured, the concentration C can be calculated.

[0005]
In pulse oximetry, in order to distinguish between two species of hemoglobin, oxyhemoglobin (HbO_{2}), and deoxyhemoglobin (RHb), absorption must be measured at two different wavelengths, i.e. the probe includes two different light emitting diodes (LEDs). The wavelength values widely used are 660 nm (red) and 940 nm (infrared), since the said two species of hemoglobin have substantially different absorption values at these wavelengths. Each LED is illuminated in turn at a frequency which is typically several hundred Hz.

[0006]
The accuracy of a pulse oximeter is affected by several factors. This is discussed briefly in the following.

[0007]
Firstly, the dyshemoglobins which do not participate in oxygen transport, i.e. methemoglobin (MetHb) and carboxyhemoglobin (COHb), absorb light at the wavelengths used in the measurement. Pulse oximeters are set up to measure oxygen saturation on the assumption that the patient's blood composition is the same as that of a healthy, nonsmoking individual. Therefore, if these species of hemoglobin are present in higher concentrations than normal, a pulse oximeter may display erroneous data.

[0008]
Secondly, intravenous dyes used for diagnostic purposes may cause considerable deviation in pulse oximeter readings. However, the effect of these dyes is shortlived since the liver purifies blood efficiently.

[0009]
Thirdly, coatings like nail polish may in practice impair the accuracy of a pulse oximeter, even though the absorption caused by them is constant, not pulsatile, and thus in theory it should not have an effect on the accuracy.

[0010]
Fourthly, the optical signal may be degraded by both noise and motion artifacts. One source of noise is the ambient light received by the photodetector. Many solutions have been devised with the aim of minimizing or eliminating the effect of the movement of the patient on the signal, and the ability of a pulse oximeter to function correctly in the presence of patient motion depends on the design of the pulse oximeter. One way of canceling out the motion artefact is to use an extra wavelength for this purpose.

[0011]
A further factor affecting the accuracy of a pulse oximeter is the method used to calibrate the pulse oximeter. Usually the calibration is based on extensive empirical studies in which an average calibration curve is determined based on a high number of persons. By means of this calibration curve, which relates the oxygen saturation of blood to pulse oximeter signals, the average difference between the theory and practice (i.e. invivo measurements) is taken into account. The calibration curve typically maps the measured invivo signal to a corresponding SpO_{2 }value.

[0012]
Pulse oximeters, however, can also utilize the LambertBeer model for calculating the concentrations of the different Hb species. In this method of calibration, the measurement signals must first be transformed into signals applicable to the LambertBeer model for calculation. This transformation constitutes the calibration of the pulse oximeter, since it is the step which adapts the invivo signals to the LambertBeer theory, according to which the pulse oximeter is designed to operate. Thus, the calibration curves can also be in the form of transformations used to adapt the actual invivo measurements to the LambertBeer model.

[0013]
Transformations are discussed for example in U.S. Pat. No. 6,104,938, which discloses a calibration method based on the absorption properties of each hemoglobin component, i.e. on the extinction coefficients of blood. In this method, the effective extinction coefficients are determined for each light signal via a mathematical transformation from the extinction coefficients according to the LambertBeer theory.

[0014]
Below, the solution according to the invention is discussed with reference to a pulse oximeter utilizing the abovementioned transformations and four different wavelengths. As mentioned above, U.S. Pat. No. 6,104,938 discloses a pulse oximeter utilizing the transformations.

[0015]
FIG. 1 is a block diagram of a pulse oximeter utilizing four different wavelengths. Light from four different LEDs 10 a, 10 b, 10 c, and 10 d, each operating at a respective wavelength, passes into patient tissue, such as a finger 11. The light propagated through or reflected from the tissue is received by a photodetector 12, which converts the optical signal received into an electrical signal and feeds it to an input amplifier 13. The amplified signal is then supplied to a control unit 14, which carries out calculation of the amount of the Hbderivatives in the blood. The control unit further controls the LED drive 15 to alternately activate the LEDs. As mentioned above, each LED is typically illuminated several hundred times per second.

[0016]
When each LED is illuminated at such a high rate as compared to the pulse rate of the patient, the control unit obtains a high number of samples at each wavelength for each cardiac cycle of the patient. The value of these samples (i.e. the amplitude of the received signal) varies according to the cardiac cycle of the patient, the variation being caused by the arterial blood, as mentioned above. The control unit 14 therefore utilizes four measurement signals, as shown in FIG. 2, each being received at one of the wavelengths.

[0017]
In order for variations in extrinsic factors, such as the brightness of the LEDs, sensitivity of the detector, or thickness of the finger, to have no effect on the measurement, each signal received is normalized by extracting the AC component oscillating at the cardiac rhythm of the patient, and then dividing the AC component by the DC component of the light transmission or reflection. The signal thus obtained is independent of the abovementioned extrinsic factors. Thus in this case the control unit utilizes four normalized signals, which are in the following denoted with
$d\text{\hspace{1em}}{A}_{i}=\frac{A\text{\hspace{1em}}{C}_{i}}{D\text{\hspace{1em}}{C}_{i}},$
where i is the wavelength in question (in this basic embodiment of the multiwavelength pulse oximeter i=1, 2, 3, 4), AC_{i }is the AC component at wavelength i, and DC_{i }is the DC component at wavelength i. The signals dA_{i }are also referred to below as modulation signals. The modulation signals thus indicate how absorption is affected by the arterial blood of the patient.

[0018]
The abovedescribed measurement arrangement corresponds to a conventional fourwavelength pulse oximeter. The operation of the pulse oximeter is discussed in more detail below.

[0019]
The theory of pulse oximetry is generally presented as being based on the LambertBeer Law. According to this theory, light transmission through the tissue at each wavelength is exponentially dependent on the absorbance of the tissue (Eq. 1). This theory is generally accepted and established as a good model for pulse oximetry.

[0020]
Next to be discussed is the theory and formalism on which the method of the invention is based.

[0021]
According to the LambertBeer theory and for a system of two analytes, the signals described above can be presented as follows:
dA _{1} =dA×(ε
_{1} ^{HbO} ^{ 2 } ×HbO _{2}+ε
_{1} ^{RHb} ×RHb)
dA _{2} =dA×(ε
_{2} ^{HbO} ^{ 2 } ×HbO _{2}+ε
_{2} ^{RHb} ×RHb)
dA _{3} =dA×(ε
_{3} ^{HbO} ^{ 2 } ×HbO _{2}+ε
_{3} ^{RHb} ×RHb)
dA _{4} =dA×(ε
_{4} ^{HbO} ^{ 2 } HbO _{2}+ε
_{4} ^{RHb} ×RHb)
RHb=1
−HbO _{2 }

 where dA is a common factor which depends on the absolute values, i.e. inter alia on the total amount of hemoglobin, ε_{i} ^{HbO} ^{ 2 }is the extinction coefficient of oxyhemoglobin at wavelength i (i=14), ε_{i} ^{RHb }is the extinction coefficient of deoxyhemoglobin at wavelength i, HbO_{2 }is the concentration fraction of oxyhemoglobin, and RHb is the concentration fraction of deoxyhemoglobin.

[0023]
Using a matrix notation, the above dependencies can be expressed for a system of n wavelengths and n analytes as follows:
$\begin{array}{cc}\left(\begin{array}{c}d\text{\hspace{1em}}{A}_{1}\\ d\text{\hspace{1em}}{A}_{2}\\ \dots \\ d\text{\hspace{1em}}{A}_{n}\end{array}\right)=C*\left(\begin{array}{c}{\varepsilon}_{11}\text{\hspace{1em}}\dots \text{\hspace{1em}}{\varepsilon}_{1n}\\ {\varepsilon}_{21}\text{\hspace{1em}}\dots \text{\hspace{1em}}{\varepsilon}_{2n}\\ \dots \\ {\varepsilon}_{\mathrm{n1}}\text{\hspace{1em}}\dots \text{\hspace{1em}}{\varepsilon}_{\mathrm{nn}}\end{array}\right)\xb7\left(\begin{array}{c}{\mathrm{HbX}}_{1}\\ {\mathrm{HbX}}_{2}\\ \dots \\ {\mathrm{HbX}}_{n}\end{array}\right),& \left(2\right)\end{array}$
where dA_{i }is the differential change in absorption (i.e. the modulation signal) at wavelength λi, εij is the extinction coefficient of the hemoglobin derivative HbX_{i }at wavelength λi, and the constant C accounts for the change of units to fractional percentages of the concentrations of the analytes HbX_{j}.

[0024]
FIG. 3 shows the extinction coefficients (ε^{HbO} ^{ 2 }and ε^{RHb}) of oxyhemoglobin (HbO_{2}) and deoxyhemoglobin (RHb) as a function of the wavelength. Point P shown in the figure is the isobestic point of oxyhemoglobin (HbO_{2}) and deoxyhemoglobin (RHb). The point has the special property that the modulation signal at the wavelength in question does not depend on the respective proportions (relative concentrations) of the hemoglobin species. Thus at the wavelength of point P the effect of the relative concentrations of oxyhemoglobin and deoxyhemoglobin on the result of the measurement is nil. It should be noted, however, that the modulation signal is independent of the relative concentrations only, not of the absolute concentrations. Thus, the absolute amount of the hemoglobin species has an effect on the result of the measurement.

[0025]
As is known, there is a difference between the LambertBeer theory and the practical measurements. The difference is due to the fact that the LambertBeer theory does not take into account the scattering and nonhomogeneity of the tissue, whereas the actual extinction coefficients are also dependent on the scattering of light caused by the tissue and blood, and on the combined effect of absorption and scattering. The larger the proportion of the attenuation caused by absorption and scattering, the larger is the correction needed between the actual and the theoretical (nonscatter) domains. This correction between these two domains can be represented by the transformation curves discussed above, by means of which the actual invivo measurements are mapped to the LambertBeer model.

[0026]
The transformation can be expressed, for example, as follows:
N _{kl} ^{LB} =g _{M} ^{−1}(
N _{kl} ^{invivo}) (3)

 where
${N}_{\mathrm{kl}}=\frac{d{A}_{k}}{d{A}_{l}}$
is the modulation ratio (the superscript indicating the domain) and the subscripts k and l indicating the wavelengths in question), and g is the transformation, for instance in the form of a polynomial function, transforming the LB Nvalues to the corresponding Nvalues in the invivo domain. The g^{−1 }in Eq. 3 is the inverse transformation, i.e. the inverse function, for transforming the measured invivo values to the ideal, nonscatter, values in the LB domain

[0028]
FIGS. 4 a to 4 f illustrate the average transformation curves measured for a pulse oximeter, where the two wavelengths for measuring the two species of hemoglobin are 660 nm and 900 nm and the third wavelength is either 725 nm or 805 nm. FIGS. 4 a to 4 c illustrate the transformation curves for a pulse oximeter with the third wavelength being 725 nm, and FIGS. 4 d to 4 f illustrate the transformation curves for a pulse oximeter with the third wavelength being 805 nm. Each curve shows the LambertBeer N_{k,l }as a function of the invivo N_{kl }at wavelengths k and l.

[0029]
FIG. 5 is a flow diagram describing the general measurement principle described in U.S. Pat. No. 6,104,938. In this method, the abovementioned N_{kl} ^{invivo }values are first determined from the dA_{i }values measured (step 51). The average transformations g_{kl }are then used to convert the measured invivo values to values N_{kl} ^{LB}, which can be used in the ideal LambertBeer model (step 52). Other input values needed for the LambertBeer model are also determined (step 53). In practice these input values are the ideal (nominal) extinction coefficients of the analytes to be measured, the extinction coefficients being given for the center wavelengths used in the measurement. The converted transformation values and the nominal input values (i.e. nominal extinction coefficients) are then used according to the LambertBeer model to calculate the concentrations of the desired analytes (step 54). Thus in this approach the invivo values N_{kl} ^{invivo }measured from the tissue are converted to the ideal invitro (cuvette) environment, where the ideal oximetry model (i.e. the LambertBeer model) is applied to yield the desired concentrations.

[0030]
In the standard two wavelength pulse oximetry the prior art technique is to map the modulation ratio N_{kl} ^{invivo }directly to the SpO2 percentage measured. In this simple case the transformation is not necessary, though the transformation technique together with the solution in the LambertBeer domain can be utilized as well.

[0031]
There are two basic ways to determine the average transformation, a theoretical approach and an empirical approach. In the empirical approach the measurements are made in the tissue by taking blood samples and measuring the actual proportions of the hemoglobin species and then determining the value of N_{kl} ^{LB }on the basis of the measured proportions. The transformation is then obtained as the relationship between the values based on the blood samples and the values given by empirical measurements as measured by the pulse oximeter. The theoretical approach, in turn, is based on a known tissue model, which takes into account the characteristics of the tissue as referred to above, which are ignored in the LambertBeer model. A first value is determined for invivo N_{kl }by means of the tissue model and a second value on the basis of the LambertBeer model. The tissue parameters of the model are determined so that the known 2wavelength calibration (socalled Rcurve) is reproduced. Then using these tissue parameters and the wavelength dependence of the tissue model, the relation of the invivo N_{kl }and the LambertBeer N_{kl }is extrapolated to other wavelengths in order to obtain the transformations at these new wavelengths. Thus in the theoretical approach no new empirical measurements are made.

[0032]
In practice the transformation can be a quadratic equation yielding a correction of the order of 20 percent to the measured N_{kl} ^{invivo }value, for example. As discussed below, the transformation data (i.e. the transformation curves) are preferably stored in numeric form in the pulse oximeter or the sensor. The number of transformation curves stored in the pulse oximeter can vary, depending on the number of wavelengths used, for example. Typically there is a transformation curve for each wavelength pair.

[0033]
As mentioned above, the accuracy of a pulse oximeter utilizing an average transformation is not necessarily sufficient, especially if analytes which are weak absorbers are to be measured or if two analytes absorb similarly, whereby it is difficult to distinguish the said analytes from each other.

[0034]
Further, each patient (i.e. subject of the measurement) has a calibration curve of his or her own, which deviates from the average calibration curve calculated on the basis of a high number of patients. This is due to the fact that for each patient the characteristics of the tissue through which light is transmitted deviate from those of an average patient.

[0035]
This causes one drawback of the current pulse oximeters; they are incapable of taking this human variability into account. Human variability here refers to any and all factors causing patientspecific variation in the calibration curve, including timedependent changes in the calibration curve of a single patient. As discussed in the abovementioned U.S. Patent, subjectdependent variation can also be seen as an effect of a third substance, such as a third hemoglobin species in the blood. However, the variation can also be interpreted as a subjectdependent change in the calibration curve of the pulse oximeter.

[0036]
Without compensation for human variability, the accuracy of current pulse oximeters is about ħ2% SpO2. However, in multiwavelength applications in general, and especially if weak absorbers, such as COHb, are to be measured, the human variability represents a much more serious problem. Therefore, techniques of compensation for these inaccuracies are called for.

[0037]
It is an objective of the invention to bring about a solution by means of which the effects caused by the tissue of the subject can be taken into account when a pulse oximeter is calibrated. In other words, it is an objective of the present invention to create a pulse oximeter which can take into account the differences caused by an individual subject as compared to the average calibration or transformation curve which the current pulse oximeter relies on.

[0038]
A further objective of the invention is to bring about a generalpurpose solution for the compensation of inaccuracies caused by human variability in pulse oximetry, a solution which is not limited to the particular general calibration method employed in the pulse oximeter, but which can be applied to any pulse oximeter regardless of its current builtin calibration method.
SUMMARY OF THE INVENTION

[0039]
These and other objectives of the invention are accomplished in accordance with the principles of the present invention by providing a mechanism by means of which the subjectspecific deviation in the tissueinduced effects on the accuracy of the pulse oximeter can be taken into account. Thus, the accuracy of the pulse oximeter is improved by taking into account the subjectspecific light transmission through the tissue, and changing the values input to the ideal model, i.e. the nominal transformation and the nominal extinction coefficients, on the basis of the measurement to compensate for the subjectspecific changes.

[0040]
In the method of the invention, the effect of tissue is taken into account and the inaccuracies caused by subjectspecific variation in that effect are compensated for. This is implemented by defining a nominal calibration for the apparatus and making initial characterization measurements in order to define the characteristics which describe the conditions under which the nominal calibration has been defined. Reference data indicating the characteristics are stored for subsequent invivo measurements in which light transmission through the actual tissue of the patients is measured. (Initial characterization measurements here refer to the measurements performed before the apparatus is taken into use. The term is used to refer to A the characterization measurements without tissue, i.e. mainly characterization of the optical components of the sensor, and B the characterization measurements of the tissues in volunteered or hospitalized subjects for which the nominal calibration of the oximeter is established.) Individual subjectspecific calibration is then defined based on the nominal calibration, and the reference data created in connection with the initial characterization measurements in the subject group in the nominal calibration, the invivo measurements in an individual patient and the invivo characterization measurements, defining the tissue characteristics of the individual patient in the invivo measurement. (Invivo characterization measurements here refer to the characterizations performed when the apparatus is in actual use.) The invivo characterization also includes a step in which the information from the optical properties of the particular sensor, used in the invivo measurement of the individual patient, is read into the oximeter. Thus, the inaccuracies are eliminated by means of comparing the optical properties of the sensors and the characteristics of the tissues in the calibration measurements and the invivo measurements in the individual patient. Thus the initial characterization measurements are used to create the reference data so that light transmission measured subsequently through the tissue of a subject can be used to correct the nominal calibration for that particular subject.

[0041]
Thus in one aspect the invention provides a method for compensating for subjectspecific variability in an apparatus intended for noninvasively determining the amount of at least two light absorbing substances in the blood of a subject and being provided with emitter means for emitting radiation at a minimum of two different wavelengths and with detector means for receiving the radiation emitted, the method comprising the steps of

 carrying out initial characterization measurements, said measurements to include the measuring of radiation received by the detector,
 based on the initial characterization measurements, establishing nominal characteristics describing conditions under which the nominal calibration is established,
 calibrating the apparatus using a nominal calibration,
 storing reference data indicating the nominal characteristics and nominal calibration,
 performing invivo characterization measurements on a living tissue, said measurements to include the measuring of radiation emitted through the tissue and received by the detector means is measured,
 performing simultaneously with the invivo characterization measurement measurements, wherein the pulsative light absorption is measured,
 based on the invivo measurements, establishing characteristics describing conditions under which the invivo measurement is done,
 based on the invivo characteristics and the reference data stored, determining tissueinduced changes in the nominal characteristics, and
 compensating for subjectspecific variation in the invivo measurements by correcting the nominal calibration on the basis of the tissueinduced changes.

[0051]
In a preferred embodiment of the invention the method is divided in two steps so that the first step compensates for the inaccuracies caused by tissueinduced and sensorinduced wavelength shift and the second step compensates for the inaccuracies caused by internal effects occurring in the tissue. The first step is then used to correct the extinction coefficients of the blood analytes to be measured, and the second step is used to correct the average transformation stored in the pulse oximeter.

[0052]
In a further preferred embodiment of the invention the effect of the temperature is also compensated for in connection with the first step.

[0053]
The method is not limited to pulse oximeters explicitly using the transformations, but can be applied to any pulse oximeter. However, the method is preferably applied to a pulse oximeter based on a transformation, since in a preferred embodiment the method is implemented by carrying out changes separately in the transformation and in the extinction coefficients.

[0054]
In another aspect, the invention provides an apparatus for noninvasively determining the amount of at least two light absorbing substances in the blood of a subject, the apparatus comprising

 emitter means for emitting radiation at a minimum of two different wavelengths,
 detector means for receiving said radiation at each of said wavelengths and producing at least two electrical output signals,
 first signal processing means for processing said output signals and producing a modulation signal for each wavelength, whereby each modulation signal represents the pulsating absorption caused by the arterialized blood of the subject,
 second signal processing means for applying a predetermined calibration on said modulation signals, whereby transformed modulation signals applicable in the LambertBeer model are obtained,
 memory means for storing and reading reference data indicating nominal characteristics under which said predetermined calibration has been applied,
 first compensation means, operatively connected to the memory means, for determining tissueinduced changes in the nominal characteristics,
 second compensation means, operatively connected to the first compensation means, for defining a subjectspecific calibration by correcting the predetermined calibration on the basis of the tissueinduced changes, and
 calculation means, responsive to the second compensation means, for determining said amounts, and
 display means.

[0064]
In a still further aspect, the invention provides a sensor for collecting measurement data for a pulse oximeter intended for noninvasively determining the amount of at least two light absorbing substances in the blood of a subject, the sensor comprising

 emitter means for emitting radiation at a minimum of two different wavelengths,
 detector means for receiving said radiation at each of said wavelengths and for producing at least two electrical output signals,
 storage means including nominal calibration and reference data indicating nominal characteristics describing calibration conditions of the pulse oximeter, said data allowing apparatus connected to the sensor to determine tissueinduced changes in the nominal characteristics when radiation is emitted through said tissue.

[0068]
Preferred embodiments of the invention are discussed in more detail below.
BRIEF DESCRIPTION OF THE DRAWINGS

[0069]
In the following, the invention and its preferred embodiments are described more closely by referring to the appended drawings, wherein:

[0070]
FIG. 1 illustrates the basic embodiment of a pulse oximeter according to the present invention;

[0071]
FIG. 2 illustrates the signals utilized in the pulse oximeter of FIG. 1;

[0072]
FIG. 3 shows the extinction coefficients of two different species of hemoglobin as a function of wavelength;

[0073]
FIG. 4 a to 4 f illustrate the average transformation curves for two different pulse oximeters;

[0074]
FIG. 5 is a flow diagram illustrating the prior art calibration method;

[0075]
FIG. 6 illustrates an example of the transmission curve of human tissue, the curve being employed in the compensation of the inaccuracies caused by tissueinduced wavelength shift;

[0076]
FIG. 7 is a flow diagram illustrating the general principle according to the present invention;

[0077]
FIG. 8 is a flow diagram illustrating steps of initial characterization measurement phase;

[0078]
FIG. 9 depicts stages of the phase of establishing nominal characteristics;

[0079]
FIG. 10 is a block diagram depicting different steps of the first stage in FIG. 9;

[0080]
FIG. 11 is a block diagram depicting different steps of the second stage in FIG. 9;

[0081]
FIG. 12 shows compensation of human variability in the invivo measurement;

[0082]
FIG. 13 depicts received intensity as a function of wavelength;

[0083]
FIG. 14 depicts frequency baseline fluctuations of a plethysmographic wave signal;

[0084]
FIG. 15 illustrates an embodiment of a sensor according to the invention;

[0085]
FIG. 16 illustrates main frame blocks of a pulse oximeter, and

[0086]
FIG. 17 illustrates main blocks of a sensor.
DESCRIPTION OF THE INVENTION

[0000]
Guidelines for Implementing the Invention

[0087]
The method of the present invention is implemented in the control unit of the pulse oximeter on the basis of the four modulation signals described above, i.e. the novelty of the system resides within the control unit itself. However, to be able to perform the selfcalibration in conjunction with each patient, the control unit requires some precalculated data, which is stored in the memory of the pulse oximeter. Instead of being stored in conjunction with the control unit, this data, or at least part of it, can also be stored in the sensor part of the pulse oximeter. The sensor part, including at least the LEDs and the photo detector, is connected to the signal processing part, which includes the control unit. Consequently, depending on the overall configuration, the novelty can also reside partly in the sensor.

[0088]
Human tissue can influence the accuracy of a pulse oximeter by two different mechanisms. First a direct wavelength shift is caused in the LED emission due to the filtering effect of the tissue. Namely, on one side of the LED center wavelength the absorption may be larger than on the other, whereby the center wavelength of the transmitted or reflected light is effectively shifted towards the region with smaller absorption. The second mechanism is a subtle one. It arises from the fact that the arterial blood is in interaction with the surrounding tissue, which can either increase or decrease the effective path length through the arterial blood layer. The first mechanism is in this context termed the external mechanism, since it affects factors external to the tissue (wavelength). The second mechanism is called the internal mechanism, as it is caused by internal factors in the tissue itself.

[0089]
Therefore, adding two compensating processes to the prior art mechanism preferably compensates for the subjectspecific variations in the influence of tissue. In other words, prior art nominal calibration is corrected with two compensating process. The first process attends to the subjectspecific variation in the external mechanism, and the second process attends to the subjectspecific variation in the internal mechanism. The first process preferably controls the extinction coefficients to be input to the LambertBeer model, while the second process preferably controls the value of the transformation used to transform the modulation ratios N_{kl} ^{invivo }to the LambertBeer model N_{kl} ^{LB}. The linear equations with the unknown analyte concentrations are then solved in the LambertBeer model, as in the prior art method. The degree of these compensations is determined by DC light transmission through the tissue (the measured DC signal), measured both in the initial characterization and invivo characterization conditions.

[0000]
a) Compensating for Tissue Filter Effect to Nominal Extinction Coefficient Values

[0090]
Nominal extinction coefficients determined without tissue must be corrected by measuring DC light transmission through tissue and determining the optical characteristics of the particular sensor without tissue. Based on the withouttissue sensor characteristics new nominal extinction coefficients without tissue for the particular invivo sensor are calculated. Then in the actual invivo measurement, subjectspecific extinction coefficients, i.e. individual extinction coefficients of a patient, are calculated for each patient and for each time moment at which a change of tissue properties invivo has been observed. Finally, based on the withouttissue optical properties and the invivo tissue properties, the final extinction coefficient values are found continuously in realtime, whereupon said values are input to the LambertBeer model.

[0091]
We next discuss the mechanisms by which tissue changes the extinction coefficients (external mechanism) and the transformation functions (internal mechanism). We introduce a parameter called the Functional Light Transmission (FLT)_{i }at a wavelength i, since it is used below in order to make all DC_{i }values (measured at varying LED emission powers at the four discrete wavelengths i) comparable to each other. Using DC values comparable to each other is in practice a prerequisite for unveiling the real effect of the tissue on the measurements and the characteristics of the tissue. In order to obtain the comparable units, the DC light transmission for each LED channel (wavelength) is first measured at a certain emitter drive current, and the measured DC value is then reduced in the preamplifier to a detector current, which is normalized to an emitter current value of 1 mA. When measured without the tissue in the probe, this result is called the Current Transfer Ratio (CTR) of the probe. CTR characterizes the sensor design and the efficiency of the light transmission from the emitters to the detector. It is usually of the order of a few microAmps (of detector current) per one milliAmp (of LED current). Now the tissue (e.g. a finger) is inserted into the probe and the CTR is again measured. This result is now called the Functional Current Transfer Ratio (FCTR) because it is the CTR measured under conditions of the function of the pulse oximeter, i.e. when the tissue is in place in the probe. The FLT_{i }is then calculated for each emitter (wavelength) as follows:
FLT(emitter# k)=FCTR(1 mA−emitter−current)/CTR(1 mA−emitter−current)

[0092]
Next the CTR and the FCTR concepts will be linked to the LambertBeer absorption model and to the actual measured intensities in the pulse oximeter. The CTR obviously describes how the external probe design factors, such as the color and geometry of the probe, affect the light transmission to the detector. On the other hand, the FLT can be associated with the true transmission through the tissue in units which are normalized to the emitter efficiency. Therefore, Eq. 1 can be written in a slightly different form, as it is often written in transport theory:
I=I _{0 }exp(−α*
d)=
I _{0 }exp(−α
_{int} *d)exp(−α
_{est} *d) (4),

 where d is the tissue thickness and a is an effective absorption coefficient. The above equation has been divided into two components. The attenuation factor with α_{int }accounts for all internal absorption effects, such as blood and can be associated to the FLTvalue as the FLT equals one when no tissue (no internal attenuation) is in the probe, and the factor with α_{ext }accounts for all external attenuations, such as geometrical factors and multiple surface reflections without light penetration into the tissue, and can be associated with the CTR of the probe. (The term d′ denotes the ‘phantom’ absorption thickness parameter for the external effects.) The term α_{ext }is mainly a SpO2 probe design issue which does not influence the measurement accuracy as such, and thus it need not to be compensated for by any means. The FLT at wavelength k can now be defined as:
$\begin{array}{cc}{\mathrm{FLT}}_{k}=\frac{I}{{I}_{0}\text{\hspace{1em}}\mathrm{exp}\left({\alpha}_{\mathrm{ext}}*{d}^{\prime}\right)}=\mathrm{exp}\left({\alpha}_{\mathrm{int}}*d\right)={\mathrm{FCTR}}_{k}/{\mathrm{CTR}}_{k}.& \left(5\right)\end{array}$

[0094]
The FLT thus describes light attenuation caused by the tissue, and it can be related to the DC light transmission in the pulse oximeter.

[0095]
In the following the compensations are discussed in more detail. The compensation of subjectvariability causing wavelength shift type interference (i.e. external mechanism) is discussed first.

[0096]
In the LambertBeer model (see Eq. 2) the effective extinctions _{εij} ^{effective }for broadband emitters, such as LEDs, can be calculated as follows:
$\begin{array}{cc}{\varepsilon}_{\mathrm{ij}}^{\mathrm{effective}}=\frac{1}{W}\text{\hspace{1em}}{\int}_{\Delta \text{\hspace{1em}}\lambda}^{\text{\hspace{1em}}}{\varepsilon}_{j}\left(\lambda \right)*{\mathrm{LED}}_{i}\left(\lambda \left(T\right)\right)*\mathrm{DET}\left(\lambda \right)*\mathrm{tissue}\left(\lambda \right)\text{\hspace{1em}}\partial \lambda ,& \left(6\right)\end{array}$
where the integration is over the LED emission spectrum LED_{i}(λ(T)), DET(λ) represents the spectral sensitivity of the detector, tissue(λ) is the spectral transmission of light through the tissue, εj(λ) is the spectral extinction of the analyte in question, T is the temperature, and W=┘LED*DET*tissue*δλ represents a normalization factor.

[0097]
In a preferred embodiment of the invention, the radiation emitting means are Light Emitting Diodes (LED), but lasers emitting at one single wavelength are also possible. For lasers the effective extinction values are the extinction values at the laser wavelength, which can depend, however, on the temperature of the emitter component. In the case of a laser, Eg. 6 is thus not needed to calculate the effective extinction value. In the preferred embodiment of the invention the emitter and detector means are located at the tissue site at which the radiation is transmitted through the tissue, but the radiation can also be conducted to and from the tissue site in a light conducting fiber or in equivalent conduction means. In this case Eg. 6 shall also include a term for the spectral transmission of the radiation conductor. A sensor utilizing light conducting fibers can be as shown in FIG. 5 of the abovementioned U.S. Pat. No. 6,104,938.

[0098]
The extinction coefficients can thus be calculated according to this equation by determining all the above factors, which depend on the actual wavelength values, i.e. the optical properties of the sensor components and the tissue term. However, as the task of determining the exact spectral value of Eg. 6 is not possible in connection with a realtime pulse oximeter measurement using only a few discrete wavelength bands, in practice the result of Eg. 6 has to be approximated. The compensation is based on determining nominal extinction coefficients and approximating their wavelength dependence in advance at the factory and using this information in the real measurement situation to approximate the final subjectspecific extinction coefficients.

[0099]
The compensation algorithm will now be presented for a 4wavelength pulse oximeter according to FIG. 1, having four LEDs at nominal wavelengths of 627 nm, 645 nm, 670 nm, and 870 nm. The extinction matrix for RHb (first column), HbO_{2}, HbCO, and metHb (last column) and for the above four wavelengths (627 nm on top) is then nominally in L/(mmol*cm).
$\begin{array}{cc}{E}_{\mathrm{kl}}^{0}=\left(\begin{array}{cccc}1.132& 0.1799& 0.2734& 3.575\\ 0.9182& 0.1124& 0.1337& 2.411\\ 0.7353& 0.0885& 0.0550& 0.5796\\ 0.2071& 0.2772& 0.010& 0.5754\end{array}\right)& \left(7\right)\end{array}$

[0100]
This equation (7) describes the nominal extinction matrix for the particular sensor used in the invivo measurement. Thus the changes of the optical properties of the sensor components with respect to the sensor components in the nominal calibration can thus be directly incorporated into this extinction matrix, a new nominal extinction matrix for the particular sensor.

[0101]
The above extinction coefficients have been calculated applying Eg. 6 at nominal LED drive temperature without the tissue filtering term tissue(λ). It then represents a nominal extinction matrix for a SpO2 sensor before its attachment on the tissue site. This extinction matrix is then altered on the basis of the measured filtering effect caused by tissue, when the sensor is attached on the site.

[0102]
It is now assumed that the spectral tissue transmission is as presented in FIG. 6, which shows spectral characteristics of tissue in the same units for each wavelength, i.e. FLT as a function of the wavelength, based on an empirical measurement. In a continuous realtime SpO2 measurement, the transmission is known at 4 distinct wavelength values (the FLT values derived from the DC values in the pulse oximeter) marked in the figure. At each wavelength the slope of the tissue transmission curve can be determined or approximated using the four transmission values. The slope then determines the change in the tissue transmission in a band of a predetermined width (100 nm in this example) around the center of the LED band. We denote the slopes between 627 to 645 nm and 645 to 670 nm by A and B, respectively. This definition of the slopes is expressed as:
$\begin{array}{c}A=\frac{\mathrm{FLT}\left({\lambda}_{2}\right)\mathrm{FLT}\left({\lambda}_{1}\right)}{\left({\lambda}_{2}{\lambda}_{1}\right)*\left(\mathrm{FLT}\left({\lambda}_{1}\right)+\mathrm{FLT}\left({\lambda}_{2}\right)\right)/2}*100\\ \mathrm{and}\\ B=\frac{\mathrm{FLT}\left({\lambda}_{3}\right)\mathrm{FLT}\left({\lambda}_{2}\right)}{\left({\lambda}_{3}{\lambda}_{2}\right)*\left(\mathrm{FLT}\left({\lambda}_{2}\right)+\mathrm{FLT}\left({\lambda}_{3}\right)\right)/2}*100,\end{array}$
where FLT(λi) is the measured FLT value determined at wavelength λi. The estimation of these slopes can be improved by calculating the curvature at the center LED (645 nm). This curvature (change of the slope/nm) is
$\mathrm{curv}=\frac{BA}{\left({\lambda}_{3}{\lambda}_{1}\right)/2}.$

[0103]
Finally the expressions are obtained for the slopes s at the three red wavelengths using A and B as parameters:
$\begin{array}{cc}\left(\begin{array}{c}{s}_{{\lambda}_{1}}\\ {s}_{{\lambda}_{2}}\\ {s}_{{\lambda}_{3}}\\ {s}_{{\lambda}_{4}}\end{array}\right)=\left(\begin{array}{c}A\mathrm{curv}*\left({\lambda}_{2}{\lambda}_{1}\right)/2\\ \left(A+B\right)/2\\ B+\mathrm{curv}*\left({\lambda}_{3}{\lambda}_{2}\right)/2\\ 0.5\end{array}\right),& \left(8\right)\end{array}$
where the slope at the IR wavelength has been estimated to be constant as it cannot be determined by the other LEDs. If we had had another LED, at about 8001000 nm range, for example, it could have been used for the estimation of the IR slope. Because the extinction curves are very flat at 870 nm and the transmission is usually rather high, the tissue prefilter cannot alter the effective extinction coefficient from its nominal value significantly. The approximation of a constant transmission slope is thus considered sufficient.

[0104]
In principle these slopes could be inserted in Eg.6 in order to integrate the new true values for the extinction coefficients. However, this is impractical to do in realtime, so a simpler algorithm is presented below. We first calculate offline the relative changes of extinction coefficients for each analyte of the system using Eg.6 and assuming that the value of the slope equals a predetermined value, which is one in this example. This calculation (assuming the slope is 1) results in a shift matrix of Eq. 9:
$\begin{array}{cc}{\mathrm{Tissue}}_{\mathrm{SHIFT}}^{\mathrm{SLOPE}=1}=\left(\begin{array}{cccc}0.975& 0.942& 0.940& 0.999\\ 0.984& 0.966& 0.920& 0.916\\ 0.963& 0.986& 0.896& 0.789\\ 1.01& 1.03& 0.903& 1.05\end{array}\right),& \left(9\right)\end{array}$

 where the effective extinction of HbO_{2 }at 645 nm is 0.966 times the original value, and the effective extinction of HbCO at 670 nm is 0.896 times the original value, for example. The proportional change of the extinction coefficient is thus the matrix value minus one, i.e. (Tissue^{SLOPE−1}−1). Thus, the matrix of Eq. 9 defines the relative changes caused by the tissue, assuming that the slope of the tissue transmission curve equals one. During invivo measurement, the slope is continuously estimated using the DC values. The ratio of the slopes then indicates the relative change of a coefficient. In other words, if the relative slope is s, the relative change is s*(Tissue^{SLOPE−1}−1). The relative changes are different for different analytes since the extinction coefficients of the different analytes behave differently as a function of wavelength.

[0106]
In realtime the effective extinction coefficients can thus be calculated as follows:
$\begin{array}{cc}{E}^{\mathrm{Eff}}={E}_{\mathrm{kl}}^{0}\otimes \left(1+S\xb7\left({\mathrm{Tissue}}_{\mathrm{SHIFT}}^{\mathrm{SLOPE}=1}1\right)\right),& \left(10\right)\end{array}$
where S denotes the column array in Eq. 8 and the matrix multiplications are performed element by element ({circle over (x)}) or element by row (•), respectively.
b) Compensating for Temperature Effect to Extinction Coefficient Values

[0107]
Changes that the external temperature and the LED drive power induce to the nominal extinction coefficients must also be corrected. If the LEDs are not driven at the nominal drive currents, their effective wavelength may also be shifted by the temperature change at the LED pn junction. The wavelength shift induced by temperature is typically about 0.10.2 nm/° C. which is significant if the drive currents are high, as is usually the case at wavelengths shorter than 660 nm. Thus, the extinction matrix of Eq. 10 must also be compensated for in varying LED drive conditions.

[0108]
There are many ways to find out the temperature of the LED pn junction. One alternative is to add a temperature sensor on the LED substrate and use the reading of the sensor for the compensation of all LED emission wavelengths. Though the junction temperature follows the substrate temperature according to some empirical heat conduction model, the method may be unreliable because the LED chip contact to the substrate and the internal heat conductivity may vary considerably. A better way is therefore to determine the junction temperature directly from the forward voltage drop of the LED junction. The junction has typical diode characteristics, which can be determined offline for each LED separately after assembling the LEDs on the substrate. It is even possible to measure, with an optical spectrometer the shift of the emission as a function of the LED forward voltage. Relating the wavelength shift to the forward voltage assumes that the forward voltage is measured during the operation of the pulse oximeter. The circuit board of the pulse oximeter should thus preferably have means for performing the forward voltage measurement. But if it does not, the LED emission shift can be calibrated against the temperature sensor at the substrate. The LED manufacturer specifications for the temperature shift can then be used to calculate the corresponding wavelength shift. Still another practical compensation for the emitter temperature changes is to map empirically the relationship of the emitter drive current to the observed wavelength shift and to use this information to adjust the invivo extinction coefficients for the sensor.

[0109]
A method for temperature compensation of the LED emission is now presented, assuming that the LED forward voltage is measured on the circuit board. The wavelength shifts can then be calculated as follows
$\begin{array}{cc}\left(\begin{array}{c}\Delta \text{\hspace{1em}}{\lambda}_{1}\\ \Delta \text{\hspace{1em}}{\lambda}_{2}\\ \Delta \text{\hspace{1em}}{\lambda}_{3}\\ \Delta \text{\hspace{1em}}{\lambda}_{4}\end{array}\right)=\left(\begin{array}{c}{k}_{1}\\ {k}_{2}\\ {k}_{3}\\ {k}_{4}\end{array}\right)\xb7\left(\begin{array}{c}\Delta \text{\hspace{1em}}{V}_{1}\\ \Delta \text{\hspace{1em}}{V}_{2}\\ \Delta \text{\hspace{1em}}{V}_{3}\\ \Delta \text{\hspace{1em}}{V}_{4}\end{array}\right),& \left(11\right)\end{array}$

 where the shift coefficients k_{i }are values determined empirically in advance and ΔV_{i }are the measured changes of the forward voltage drops. For the 627645670870 nm LEDs of the sensor, the kvalues are 0.06 nm/mV, 0.06 nm/mV, 0.09 nm/mV, and 0.1 nm/mV, respectively.

[0111]
As in the compensation discussed above relating to tissue filtering, it is practical to first calculate the change in the extinction coefficients offline for a certain fixed wavelength shift. In this example the relative changes of the extinction coefficients are calculated, as in Eq. 9, for a 5 nm wavelength shift for each of the four hemoglobin derivatives. The following shift matrix is then obtained:
$\begin{array}{cc}{\mathrm{Temp}}_{\mathrm{SHIFT}}^{\Delta \text{\hspace{1em}}\lambda =5\text{\hspace{1em}}\mathrm{nm}}=\left(\begin{array}{cccc}0.919& 0.820& 0.798& 0.974\\ 0.963& 0.926& 0.823& 0.794\\ 0.941& 0.983& 0.855& 0.725\\ 1.0& 1.01& 0.963& 1.02\end{array}\right).& \left(12\right)\end{array}$

[0112]
During invivo measurement, which will be later applied to a patient in a hospital or like, the relative changes are then calculated based on the measured wavelength shift. The ratio of the wavelength shifts then indicates the relative change of a coefficient caused by temperature. In other words, if the relative change calculated for a wavelength shift of Y1 is r, the relative change for the measured (invivo) wavelength shift of Y2 is r×(Y2/Y1). The relative changes are different for different analytes, since the extinction coefficients of the different analytes behave differently as a function of wavelength.

[0113]
The temperature compensated extinction coefficients are thus:
$\begin{array}{cc}{E}_{\mathrm{TEMP}}^{\mathrm{EFF}}={E}_{\mathrm{kl}}^{0}\otimes \left(1+\left(\Delta \text{\hspace{1em}}\lambda /5\text{\hspace{1em}}\mathrm{nm}\right)\xb7\left({\mathrm{Temp}}_{\mathrm{SHIFT}}^{\Delta \text{\hspace{1em}}\lambda =5\text{\hspace{1em}}\mathrm{nm}}1\right)\right),& \left(13\right)\end{array}$

 where Δλ is the array in Eq. 11. As mentioned earlier, Δλ can also be estimated by reading the temperature indicated by the temperature sensor on the LED substrate or by measuring the LED drive current and using the mapping of the current to the wavelength shift.

[0115]
The compensation of the variability causing wavelength shift type interference can now be summed up as follows:
$\begin{array}{cc}{E}^{\mathrm{Eff}}={E}_{\mathrm{kl}}^{0}\otimes \left(1+S\xb7\left({\mathrm{Tissue}}_{\mathrm{SHIFT}}^{\mathrm{SLOPE}=1}1\right)\right)\otimes \left(1+\left(\Delta \text{\hspace{1em}}\lambda /5\text{\hspace{1em}}\mathrm{nm}\right)\xb7\left({\mathrm{Temp}}_{\mathrm{SHIFT}}^{\Delta \text{\hspace{1em}}\lambda =5\text{\hspace{1em}}\mathrm{nm}}1\right)\right).& \left(14\right)\end{array}$

 c) Compensating for Tissue Effect to Transformation Functions

[0117]
The second compensation (the internal mechanism) controls the value of the transformation used to transform the modulation ratios N_{kl} ^{invivo }to the LambertBeer model N_{kl} ^{LB}. Therefore nominal transformation values are first calculated based on DC signals obtained from the sensor when the tissue properties are averaged over a large group of people. Then, in the actual invivo measurement subjectspecific transformation values, i.e. individual tissue characteristics affecting the transformation values of a patient are calculated. Finally, based on the invivo measurement, the nominal transformation values are corrected, whereupon the corrected values are input to the transformations used to transform the modulation ratios to the LambertBeer model.

[0118]
A practical implementation of the second compensating step is now discussed by introducing a new variable called “path length multiplier”, since this will provide an easy way of understanding the technique in accordance with the invention.

[0119]
As mentioned above, the purpose of the invention is to improve the accuracy of a pulse oximeter in situations in which the blood volume, the red blood cell density or the hematocrit, the total hemoglobin (g/dl), the division between the arterial and venous blood compartment volumes, and the arterialvenous saturation difference vary and produce human variability, which worsens the accuracy of the SpO2 measurement. It is also the purpose of the invention to compensate for the effect of skin pigmentation (dark skin), which in part can be considered to belong to the tissue prefilter category of compensations, but which also influences via modifying the path length multiplier. This modification is especially important for SpO2 ear sensors, which are attached to a very thin and pigmented tissue part (of about the same thickness corresponding to the diffusion constant in human tissue).

[0120]
The interdependence of the abovedescribed transformation and the path length multiplier is first illustrated by considering the photon path lengths through a single layer of artery blood and examining how the scattering and absorption affect it. It is postulated here that multiple scattering effectively increases the photon path length through the artery and that the absorption of the surrounding tissue effectively decreases it. In this way the artery and tissue are in interaction with each other. To derive a mathematical formulation of this, relationship, the known KubelkaMunk twoflux model can be used. This model defines an absorption probability K as follows:
$\begin{array}{cc}K=\langle \frac{dl}{dz}\rangle *{\Sigma}_{a},& \left(15\right)\end{array}$
where Σa is the macroscopic absorption crosssection of the media and dl is the true average photon path length through the scattering and absorbing medium of infinitesimal layer thickness dz. The term <dl/dz>=K/Σa is a path length multiplier (plm) which enhances the arterial blood absorption from that of the LambertBeer nonscatter value because of the multiple scattering in the surrounding medium.

[0121]
The idea of the path length multiplier is applied to the LambertBeer formulation of 2λ pulse oximetry. The ratio of the change in absorption at the two probe wavelengths is defined as:
$\begin{array}{cc}\frac{d{A}_{k}}{d{A}_{l}}={N}_{\mathrm{kl}}^{\mathrm{in}\mathrm{vivo}}=\frac{{\mu}_{a}^{k}*{d}_{k}}{{\mu}_{a}^{l}*{d}_{l}},& \left(16\right)\end{array}$
where μa
^{i }is the arterial (nonscatter) absorption coefficient at wavelength i and di is the effective true optical path length. The transformation is defined by substituting equation 15 with dl=d
_{i }in equation 16:
$\begin{array}{cc}{N}_{\mathrm{kl}}^{\mathrm{in}\mathrm{vivo}}=\frac{{\mu}_{a}^{k}*{\left(\frac{K}{{\Sigma}_{a}}\right)}_{k}*\mathrm{dz}}{{\mu}_{a}^{l}*{\left(\frac{K}{{\Sigma}_{a}}\right)}_{l}*\mathrm{dz}}=\frac{{\left(\frac{K}{{\Sigma}_{a}}\right)}_{k}}{{\left(\frac{K}{{\Sigma}_{a}}\right)}_{l}}*{N}_{\mathrm{kl}}^{\mathrm{ideal}},& \left(17\right)\end{array}$

 where the ideal LambertBeer model is used for N_{kl} ^{ideal}≡μa^{k}/μa^{l}, and where the layer thickness dz is the same for all wavelengths (k, l). Equation 17 now represents the transformation (gkl)^{−1 }from N_{kl} ^{invivo}, i.e. from the measured value, to N_{kl} ^{LB}, which is the ratio of differential absorptions that would be measured if the measurement system were the ideal cuvette system of the LambertBeer model. For the transformation the equation below is obtained:
$\begin{array}{cc}{g}_{\mathrm{kl}}=\frac{{\left(\frac{K}{{\Sigma}_{a}}\right)}_{k}}{{\left(\frac{K}{{\Sigma}_{a}}\right)}_{l}}=\frac{{\mathrm{plm}}_{k}}{{\mathrm{plm}}_{l}}.& \left(18\right)\end{array}$

[0123]
Thus the transforming quantity is a ratio of path length multipliers measured at two different wavelengths (k and l). The dependence of the function g_{kl }thus refers to the absorption density of the scattering tissue in the surrounding of the infinitesimal arterial layer dz including the layer itself. This essentially means that the transformation does not require knowledge of the analyte composition in the arterial blood, but refers rather to the macroscopic light absorption, i.e. transmission through the tissue part under the sensor. That is in the language of pulse oximetry the DC component of the light transmission. This is utilized in the compensation of the invention.

[0124]
Modifying Eq.1 and leaving the attenuation of the probe design factors (i.e. CTR values) out of consideration, the relationship of the DC light transmission through the tissue and the path length multiplier can be presented as follows:
I _{out} −I _{in} e ^{−εDC} =I _{in} e ^{−εplm D1 C} (19),

 where D is the actual path length through the sample, D1 is the shortest path length through the sample (i.e. the thickness of the sample), and ε is the ideal extinction coefficient of the analyte. Here the I_{out}/I_{in }can be associated with the FLT at the wavelength in question. Plm thus describes the internal attenuation factors in the tissue and, in particular, the enhancement of the absorbancy relative to the ideal cuvette absorption.

[0126]
In nominal conditions, the path length multiplier has a certain nominal value plm^{0 }(where the superscript ‘0’ refers to the nominal value). This nominal value can be determined in the factory at the manufacturing stage of the pulse oximeter. When the DC component is measured again in connection with invivo measurement, the change in the plm from the nominal value can be used to determine the change in the average transformation.

[0127]
The term α
_{int }in Eg.4 can be expressed with the help of the path length multiplier in the LambertBeer model as
α
_{int} =plm*Σ _{α},

 where Σa accounts for all internal absorption sources and is defined in the nonscatter LambertBeer domain. The FLT at wavelength k can then be written as follows:
FLT _{k}=exp(−α_{int} *d)=exp(−plm*Σ _{α} *d) (20).

[0129]
We then ratio the logarithms of the FLTs at two wavelengths k and l, which results in Eq. 21:
$\begin{array}{cc}\begin{array}{c}\frac{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}=\frac{{\mathrm{plm}}_{k}*{\Sigma}_{a}^{k}}{{\mathrm{plm}}_{l}*{\Sigma}_{a}^{l}}\\ ={g}_{\mathrm{kl}}*\frac{{f}_{a}*{\mu}_{a}^{k}+{f}_{v}*{\mu}_{v}^{k}}{{f}_{a}*{\mu}_{a}^{l}+{f}_{v}*{\mu}_{v}^{l}}\\ ={g}_{\mathrm{kl}}*\frac{{f}_{a}\left({\mu}_{a}^{k}{\mu}_{v}^{k}\right)+{\mu}_{v}^{k}}{{f}_{a}\left({\mu}_{a}^{l}{\mu}_{v}^{l}\right)+{\mu}_{v}^{l}},\end{array}& \left(21\right)\end{array}$

 where g_{kl}=plm_{k}/plm_{i }is the transformation between the LambertBeer and invivo modulation ratios according to Eq. 18, and in which the internal absorbing tissue compartments are venous and arterial blood with volume fractions f_{v }and f_{a }and with absorption coefficients μ_{v }and μ_{a }determined in the LambertBeer domain, respectively. In the last expression we have used for the venous volume fraction the relationship f_{v}=1−f_{a}. As the arterial volume fraction is always smaller than the venous volume fraction and as the arterialvenous absorption difference is always smaller than the venous absorption, the dominating factor in the last term is μ_{v} ^{k}/μ_{v} ^{l}, i.e. the venous saturation SvO2.

[0131]
Thus the changes in the FLT and SvO2 from their nominal values provide the compensation needed for estimating the correct transformation function g
_{kl}. We can then finally write for the relative change of the transformation function g
_{kl}:
$\begin{array}{cc}\frac{{g}_{\mathrm{kl}}}{{g}_{\mathrm{kl}}^{0}}=\frac{\frac{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}/\frac{{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}^{0}}{{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}^{0}}}{{F}_{\mathrm{kl}}\left(\mathrm{SvO2},\mathrm{SaO2},{f}_{a}\right)/{{F}_{\mathrm{kl}}\left(\mathrm{SvO2},\mathrm{SaO2},{f}_{a}\right)}^{0}},& \left(22\right)\end{array}$

 where the function F_{kl }represents the ratio term
$\frac{{f}_{a}\left({\mu}_{a}^{k}{\mu}_{v}^{k}\right)+{\mu}_{v}^{k}}{{f}_{a}\left({\mu}_{a}^{l}{\mu}_{v}^{l}\right)+{\mu}_{v}^{l}}$
in Eq. 21 and the superscript 0 represents the nominal values of the nominal calibration function g_{kl} ^{0}, which is on average true for a large patient population. In fact, the log(FLT) and the F_{kl }compensation terms account for quite different human variability factors in the tissue: whereas F_{kl }mainly tracks the changes of the arterial venous saturation difference, in particular SvO2, the log(FLT) reflects the changes in the total absorption of the tissue, i.e. in the total blood volume and the total hemoglobin or hematocrit, which are not seen in F_{kl }at all. In practice, the largest corrections to the transformation function are due to the log(FLT) and F_{kl }is less important.

[0133]
The FLT in Eq. 22 is easily obtained at the two wavelengths k and l, as has been described earlier in Eg.5. The function F_{kl }represents the ratio of the absorption coefficients (in the LambertBeer nonscatter model) of the whole tissue at these same two wavelengths, i.e. it represents the internal color of the tissue. This internal absorption ratio can be measured by examining the low frequency baseline fluctuations of the plethysmographic wave signal.

[0134]
FIG. 14 depicts frequency baseline fluctuations of the plethysmographic wave signal. These fluctuations are caused by the low frequency changes (usually of respiration origin) in the blood volume or in the blood volume distribution of the tissue. Similarly, since the arterial color (=Rratio) is defined as the ratio of the arterial absorption coefficients, function F
_{kl }can be calculated as:
$\begin{array}{cc}{F}_{\mathrm{kl}}=\frac{{f}_{a}\left({\mu}_{a}^{k}{\mu}_{v}^{k}\right)+{\mu}_{v}^{k}}{{f}_{a}\left({\mu}_{a}^{l}{\mu}_{v}^{l}\right)+{\mu}_{v}^{l}}={g}_{\mathrm{kl}}^{1}\text{\hspace{1em}}\frac{{\left(\mathrm{AC}/\mathrm{DC}\right)}_{k}}{{\left(\mathrm{AC}/\mathrm{DC}\right)}_{l}}={g}_{\mathrm{kl}}^{1}\left({N}_{\mathrm{kl}}^{\mathrm{baseline}}\right),& \left(23\right)\end{array}$ 
 where AC is the amplitude (or the instantaneous slope) of the low frequency baseline fluctuation, instead of the heart pulse amplitude of the plethysmographic wave, and DC is the DC light transmission at that particular wavelength. Because the effective tissue color is mainly determined by the venous blood, function F_{kl }can be approximated as the arterial modulation ratio calculated for the venous saturation, which is usually about SaO210% i.e. F_{kl}=R(SvO2=SaO210%).

[0136]
If the venous saturation is determined by venous blood samples and the arterial saturation by the arterial blood samples, the function F_{kl }can be calculated using the real blood values (with the assumption that the corresponding blood compartment volumes are f_{a}=0.25 and f_{v}=0.75).
PREFERRED EMBODIMENT OF THE INVENTION

[0137]
FIG. 7 is a flow diagram illustrating the general principle of the present invention. The method can be divided into two groups of phases. The first group 7173 comprises measures relating to the setting up that is carried out prior to actual use of a pulse oximeter for measuring analyte concentrations of a patient whereas the second group 7475 comprises method phases performed in the actual use. Previously in the present application the phases of the first group were also called offline phases.

[0138]
In the setting up phase, initial characterization measurements are first made, preferably at the calibration stage of the pulse oximeter with a nominal wavelength pulse oximeter sensor (phase 71). Based on the measurements, nominal characteristics are established describing the conditions under which the pulse oximeter has been calibrated (phases 72). As a result of these phases, reference data are stored (phase 73), which describe the calibration conditions of the pulse oximeter. In connection with subsequent invivo measurements, the same characteristics are again estimated and tissueinduced changes in the characteristics are determined based on the measured characteristics and the reference data stored (phase 74). In addition to the tissue induced changes, the changes relating to the different sensor than in the calibration stage are incorporated with the tissue changes. In invivo measurement, the N_{kl} ^{invivo }values are determined from the dA_{i }values measured. On the basis of the changes determined, the subjectspecific calibration is then determined (phase 75) for the invivo measurements to be performed by the pulse oximeter on the subject.

[0139]
It is to be noted here that phases 7173 are performed either when the pulse oximeter has been calibrated in a known manner or at the manufacturing stage of the pulse oximeter sensor when the sensor characteristics are determined. After these phases, the nominal transformation and the nominal extinction coefficients are known to the pulse oximeter.

[0140]
Next, referring to FIG. 811 contents of phases 71 and 72 will be explained in more detail.

[0141]
FIG. 8 depicts the first steps carried out in the initial characterization phase of the setting up. It should be noted that in this phase no living tissue is needed. These steps are performed prior to the actual measurements, for example in the factory at the manufacturing stage of the pulse oximeter sensor. Thus, referring to FIG. 8, the steps are as follows.

[0142]
Step 81.

[0143]
The spectral characteristics of the emitter/detector system are measured. In other words, the LEDs are characterized for their light emission (the emission as a function of wavelength) and the detector for its spectral sensitivity. This step thus includes determination of the characteristics of the curve shown in FIG. 13, i.e. the received intensity as a function of wavelength (at least around the wavelengths used). The light transmission from the light emitter to the light detector is measured without living tissue, i.e. the CTR is determined in the sensor in a fixed setup mimicking the actual use of the sensor. For cliptype sensors this is usually the Probe Off position of the sensor. The step also includes determination of the center wavelength of each LED.

[0144]
Step 82

[0145]
Using the spectral characteristics obtained at the previous step, the effective extinction coefficients for the nominal extinction matrix are determined without the tissue term. Thus in this step Eg.6 is used without the tissue term (tissue(λ)) to form the nominal extinction matrix E_{kl} ^{0 }according to Eg.7.

[0146]
Step 83

[0147]
The tissue correction to the nominal extinction matrix is estimated for artificial tissue, in which transmission slopes are 1. In other words, the relative changes in the effective extinction coefficients due to artificial tissue filter effect are determined. Here Eg.6 is used assuming that the slope of tissue(s) equals a predetermined fixed value at each wavelength. In other words, the shift matrix of Eq. 9 is determined.

[0148]
Step 84

[0149]
Temperature corrections to the nominal extinction matrix are determined, i.e. the relative changes in the effective extinction coefficients due to wavelength shift caused by changes in temperature are determined. In other words, the matrix of Eq. 12 is determined, which indicates the relative changes for a wavelength shift of a predetermined value. Said value could be 5 nm, for example.

[0150]
Step 85

[0151]
If the LED forward voltage method is used, the LED forward voltages are characterized at a typical drive current for small ambient temperature changes. In other words, the relationship between the forward voltage shift and the temperature shift is established for each emitter. Thus, the temperature coefficients k_{l }in Eq. 11 and the nominal forward voltage drops at nominal temperature are determined for each LED.

[0152]
Step 86

[0153]
All data obtained in the previous steps are saved in a memory unit in the sensor or in the control unit, or the corresponding information is otherwise made available to the pulse oximeter, for example, by using codes, such as sensor identification numbers, which indicate the values of the information.

[0154]
As is obvious from the above, steps 81 to 86 include performing initial characterization measurements for the compensation, said measurements to include measuring the light transmission of the apparatus, establishing nominal DC transmission characteristics of the apparatus on the basis of the measurements, and for subsequent invivo measurements storing reference data that indicate the transmission characteristics established.

[0155]
Thus, in these first initial characterization measurements, the value of Eg.6 is determined using nominal values, and the nominal extinction matrix is formed. The apparatus is also provided with the data needed in the subsequent invivo compensation steps for calculating the changes in the factors included in Eg.6. In the invivo steps the said changes are determined and a new extinction matrix is formed, whereby the new extinction values are such that the external effects are compensated for.

[0156]
To sum up, after the setting up steps the pulse oximeter stores the matrices according to equations 7, 9, and 12 and the values of the shift coefficients k_{i}. In addition to this, the oximeter stores the CTR values and the center wavelengths corresponding to these values.

[0157]
FIG. 9 depicts setting up stages that are carried out while establishing calibration and nominal characteristics, which describe calibration conditions (phase 72 in FIG. 7). These setting up stages consist of three measurement and calculation stages that are carried out simultaneously. Stage 91 includes characterization of calibration conditions when external factors influence the nominal extinction matrix that was formed previously. Stage 92 includes characterization of calibration conditions when internal factors influence transformation curves. Stage 93 includes establishment of nominal calibration either for multiwavelength pulse oximeter of twowave pulse oximeter. In performing the steps the information stored during initial characterisation measurements (phase 71 in FIG. 7) is utilized. In addition, it is worth noting that when doing the measurements in said three steps a sensor is attached to the pulse oximeter and there is living tissue between the emitters and the detector. As to living tissues, a large group of people is used in order to obtain enough statistical data. At least the nominal calibration and the nominal internal factors are stored. Optionally external factors are stored.

[0158]
FIG. 10 is a block diagram depicting different steps of stage 91 (FIG. 9). The content of the steps are as follows:

[0159]
Step 101

[0160]
The forward voltage drops (ΔV) are measured for each LED and for each subject in a group of people. The voltage drops are compared with the nominal drops obtained in step 85 of FIG. 8 or the temperature is compared with the nominal temperature. The change in the forward voltage relative to the sensor nominal values are calculated, the nominal values having been stored in the sensor memory unit or in the control unit memory. If the pulse oximeter does not have forward voltage measuring means, the temperature of each LED is estimated by reading the temperature indicated by a sensor on the LED substrate, and either the manufacturer specifications or empirical data for corresponding wavelength changes or a lookup table mapping the emitter drive current to the center wavelength shift is used.

[0161]
Step 102

[0162]
Then the change in the extinction matrix for the temperature compensation is calculated. This is done for each person of the group. Calculation is done according to Eq. 1113, i.e. using the wavelength shifts determined according to Eq. 11 and calculating the matrix of Eq. 13 using the matrices of Eg.7 and 12 stored in step 84 of FIG. 8.

[0163]
Step 103

[0164]
For each person of the group the tissue transmission induced changes to the nominal extinction matrix are calculated. Thus, the DC light transmission for each LED channel (wavelength) is measured, and the value measured is normalized to an emitter current value of 1 mA. The result is the FCTR of the sensor. An estimate for the FLT is then calculated for each emitter (wavelength). In this connection, the FLT values are used for calculating the slopes (A and B). In other words, all DC values are normalized in relation to the 1 mA emitter current in order to make all values comparable to one another. Equations 810 are used for calculating the change of the nominal extinction matrix.

[0165]
Step 104

[0166]
This step is an alternative to the next step. Namely it might be practical simply to average the results obtained in the previous steps for each person of the group and then just combine the average values (temperature and tissue effects) into one extinction matrix.

[0167]
Step 105

[0168]
If step 104 is omitted, then the effective LambetBeer extinction matrix for each person of the group is calculated. The effective LambertBeer extinction coefficients are determined using Eq. 14. Thus, each of those matrices characterizes the calibrations conditions for one particular person of the group.

[0169]
Step 106

[0170]
Averaging of the individual LambetBeer matrices element by element results in a group average extinction matrix that characterize the extinction coefficients during the nominal calibration process.

[0171]
Optionally all the data is stored. However, it may be more practical that in the subsequent invivo pulse oximeter measurements the steps 101106 are continuously repeated to update extinction matrix coefficients using only the information stored in the settingup phase of FIG. 8. When the data describing the external factors—mainly the characteristics relating to the sensor—is stored, it can be used for checking up the sensor condition and alert the user if considerable, abnormal deviation from the average characteristics occur during invivo measurements. The stored data can thus be used to issue a Probe Fault Condition alarm.

[0172]
Next, referring to FIG. 11 the content of the stage 92 (FIG. 9) is explained. As stated above, that stage is performed simultaneously with stages 91 and 93.

[0173]
Step 111

[0174]
Characterization of the calibration conditions with internal factors influencing transformation curves starts by first fetching the FLTs calculated in step 103 of FIG. 10 and then by calculating the ratio of the logarithms for each wavelength pair and for each person of the group. In other words, the nominal light transmission through a finger or ear lobe or its approximation at the distinct nominal wavelengths of the sensor, is determined for each person of a group (the group may be a population of patients or volunteers or even for only one single volunteer), on whom the nominal calibration was performed. This gives for each person of the group an individual curve
$\frac{{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}^{0}}{{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}^{0}}$
as a function of N_{kl} ^{invivo }at these wavelengths (k, l), i.e. essentially as a function of the correct SpO2.

[0175]
Step 112

[0176]
A regression curve is calculated from the individual curves obtained in step 111. The logarithm ratios fitted to the regression curve and the corresponding N_{KL }values are stored in the table 1, 5^{th }column. The table is presented below. This information is stored in the sensor memory unit (or in the control unit).

[0177]
Step 113

[0178]
Ffactors, i.e. function F
^{0 }can be determined in two alternative ways:
 a) In the above measurement using the baseline fluctuations, N_{kl} ^{baseline }is calculated and transformed to the LambertBeer model by the nominal transformation (g^{0} _{kl})^{−1 }for the group of persons (patients or volunteers or even for only one single volunteer, on whom the nominal calibration was performed). This is tabulated in table 1 below as a function of N_{kl} ^{invivo }at these wavelengths (k, l). These data are stored in the sensor memory unit (or in the control unit).
 b) In the above measurement venous and arterial blood samples are taken from a position close to the sensor site and analyzed for RHb, HbO_{2}, HbCO, and metHb. The absorption coefficients μ_{a }and μ_{v }are then calculated using the measured analyte fractions. The arterial volume fraction f_{a }is then estimated. Usually it is sufficient to approximate that f_{a }is equal to 0.25. The functions F_{kl} ^{0 }are calculated in the LambertBeer Model using the venous and arterial volume fractions and
${\mu}_{a,v}=\mathrm{RHb}*{\mu}_{a,v}^{\mathrm{RHb}}+\mathrm{HbO2}*{\mu}_{a,v}^{\mathrm{HbO2}}*\mathrm{HbCO}*{\mu}_{a,v}^{\mathrm{HbCO}}+\mathrm{metHb}*{\mu}_{a,v}^{\mathrm{metHb}}.$

[0181]
As a result of steps
111
113, we have set up a lookup table for each wavelength pair, in which the following nominal information is tabulated:
TABLE 1 


k  l  N_{kl} ^{invivo}  g_{kl} ^{0}  Log(FLT_{k})^{0}/log(FLT_{l})^{0}  g_{kl} ^{−1} × (N_{kl} ^{baseline})  F^{0}kl 







 and where also the nominal transformation g^{0 }is presented and N_{kl} ^{baseline}≅g_{kl} ^{0}×(F_{kl}). Only one of the two last columns is necessary, depending on the way the values of function F are determined. It is also to be noted that the ratio is not necessary in column 5, but that it is enough to store the FLT values from which the ratio of their logarithms can be calculated.

[0183]
As is obvious from the above, steps 111113 again include measuring the DC light transmission of the apparatus, establishing nominal DC transmission characteristics for the apparatus on the basis of the measurements, and storing reference data for subsequent invivo measurements, the data indicating the transmission characteristics established.

[0184]
After the above setting up steps, which can be performed at the manufacturing stage of the pulse oximeter, the pulse oximeter is ready for use in a hospital environment or like in connection with invivo measurements.

[0185]
Now, next we consider a case where the pulse oximeter is used to detect blood oxygenation of a patient in a hospital or like, i.e. invivo measurement is to be performed. The patient puts his finger of ear lobe between the emitters and the detector whereupon light transmission through the finger or ear lobe is measured. The sensor of the pulse oximeter produces data that is processed as shown in FIG. 12. It is worth noting that during the measurement pulsative signals are collected through the whole measuring period. Accordingly, the full round of steps that will be described below are performed cyclically many times, one round per each pulsation of the heart, for example.

[0186]
FIG. 12 shows four steps of the data processing in the invivo measurement. Firstly, changes that patient tissue causes to the nominal extinction matrix are calculated, step 121. Secondly, changes that patient tissue causes to the nominal transformation curves are calculated, step 122. Thirdly, the nominal values, which were calculated in the setting up phase, are compared with the individual patient values calculated. Equation 22 is used to get the individual transformation function for the patient (step 123). Finally, the individual transformation functions obtained from step 123 and the individual extinction matrix obtained from step 121 are applied to LambertBeer model that gives analytic concentrations of the patient.

[0187]
The steps are now explained in more detail.

[0188]
Step 121 is actually repetition of steps 101102 of FIG. 10. Thus, the forward voltage drops (ΔV) are measured for each LED. The voltage drops are compared with the nominal drops stored in the memory. The change in the forward voltages relative to the sensor nominal values is calculated. Then the change in the nominal extinction matrix for the temperature compensation is calculated. Calculation is done according to Eq. 1113, i.e. using the wavelength shifts determined according to Eq. 11 and calculating the matrix of Eq. 13 using the matrices of Eg.7 and 12. Calculation in step 121 therefore results in patient specific extinction matrix.

[0189]
Accordingly, step 122 is actually repetition of step 111 of FIG. 11. The step starts by first fetching the FLTs stored in the memory during the setting up process of the oximeter and then by calculating the ratio of the logarithms for each wavelength pair and for the patient. In other words, the nominal light transmission through the finger or ear lobe is determined for the patient. This gives for the patient a curve
$\frac{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}$
as a function of N_{kl} ^{invivo }at these wavelengths (k, l). Then in Table 1 in the row where the measured modulation ratio at the wavelengths k and l equals N_{kl} ^{invivo}, the nominal
$\frac{{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}^{0}}{{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}^{0}}$
is read. The correction factor
$\frac{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}/\frac{{\mathrm{log}\left({\mathrm{FLT}}_{k}\right)}^{0}}{{\mathrm{log}\left({\mathrm{FLT}}_{l}\right)}^{0}}$
is then calculated.

[0190]
Then in step
123, comparison between the nominal values stored in Table 1 and the individual patient specific values obtained from step
112 is done and equation (22) is used to get the patient specific transformation functions. This is carried out so that in Table 1 in the row, where the measured modulation ratio at the wavelengths k and l equals N
_{kl} ^{invivo}, the nominal F
_{kl} ^{0 }is read from either of the last two columns. Then the correction factor F
_{kl}/Fk
_{kl} ^{0 }is determined in one of the following two ways:

 a) Column g^{−1}×(N^{baseline}): By using the baseline fluctuations of the measured plethysmographic wave and by using Eq. 23 with g^{−1}, the function F_{kl }is determined. The correction factor F_{kl}/F_{kl} ^{0 }is then calculated. For determining N^{baseline }and its changes, the amplitudes of the signal can be used, as is normally done for a modulation ratio N.
 b) Column F_{kl}: The blood analytes RHb, HbO_{2}, HbCO, and metHb are solved in the LambertBeer model using the nominal transformation g^{0}. This is the first approximation for the analytes. The absorption coefficients μ_{a }and μ_{v }are calculated using the measured analyte fractions in the arterial absorption and approximating the absorption coefficient in the venous blood by using the measured dyshemoglobin fractions and setting HbO2^{vena}=HbO210% and RHb^{vena}=RHb+10%, where f_{a}=0.25 is assumed. The functions F_{kl }are calculated in the LambertBeer Model using the equation:
${\mu}_{a,v}=\mathrm{RHb}*{\mu}_{a,v}^{\mathrm{RHb}}+\mathrm{HbO2}*{\mu}_{a,v}^{\mathrm{HbO2}}*\mathrm{HbCO}*{\mu}_{a,v}^{\mathrm{HbCO}}+\mathrm{metHb}*{\mu}_{a,v}^{\mathrm{metHb}}.$

[0193]
A new transformation g_{kl }is calculated using Eq. 22, and the new transformation is used for solving the analyte concentrations in the LambertBeer model. Optionally, a more accurate estimate of F_{kl }can be obtained by iteration of new analyte fractions for the new corrected transformation.

[0194]
Finally, in step 124 the patient specific transformation functions and the patient specific extinction matrix are used to solve the analyte concentrations by applying them to the LambertBeer model.

[0195]
During the above steps 121124 invivo measurements are performed, wherein the DC component of the radiation emitted through the tissue and received by the detector is measured, tissueinduced changes in the transmission characteristics are determined based on the invivo DC component and the transmission characteristics stored, and on the basis of the tissueinduced changes the subjectspecific variation in the invivo measurement is compensated for.

[0196]
The above setting up and invivo steps compensate for the nonideal characteristics of the broadband emitters or for external effects on the light source emission spectra. They do also compensate for the variation in the absorption and scattering interplay in the tissue, i.e. the internal effects, which equally influence a single line laser emitter and a broadband LED emitter. Lasers also show shifts in the emission line wavelength as a function of the temperature. Therefore, the lasers are compensated for the temperature and internal tissue effects, but not for the prefilter tissueinduced spectral shifts.

[0197]
The precalculated data utilized by the pulse oximeter can be stored in the sensor part of the pulse oximeter, whereby the same sensor can be attached to different pulse oximeter housings.

[0198]
FIG. 15 illustrates the general structure of a sensor according to the invention, the detailed configuration of the sensor being dependent on which information is stored in the sensor and which in the signal processing part, and also on the amount of the calculation appropriate in the signal processing part.

[0199]
Nevertheless, a sensor according to the invention includes the light sources (10 a10 c) and the photo detector, the light sources being adapted to emit at two or more wavelengths. In addition, the sensor includes a data storage unit M2 for storing the data on the basis of which the signal processing part can perform the abovedescribed calibration. The information necessary for the above compensations is shown in the figure. For the first compensation process the pulse oximeter needs the kvalues, the abovementioned three matrices, i.e. the nominal extinction matrix (Eg.7), the shift matrices (Eqs. 9 and 12), and the CTR/wavelength pairs. For the second compensation process, in turn, the pulse oximeter needs the information stored in Table 1 and the CTR/wavelength pairs. As mentioned above, at least part of this data determined prior to the use of the device for invivo measurements can also be stored in the control unit part of the pulse oximeter. The apparatus further preferably includes means 150 for measuring the forward voltage of the pn junction of each LED, as discussed above.

[0200]
FIG. 16 illustrates the general structure of a main frame. The main frame refers here to hardware loaded with software, which are enclosed in the housing ox a pulse oximeter. The main frame comprises thee units: control and measurement unit 161, compensation unit 162 and hemoglobin fraction calculation unit 163.

[0201]
Control and measurement unit 161 includes light source controller 1611 that forms and controls supply current to the emitters of the light source. The emitters locate in the separate sensor (see FIG. 15). Detector signal receiver 1612 receives the signal from the photo detector, filters and amplifies said signal. The amplified analog signal is converted to digital domain in A/D converter 1614 whereupon AC and DC components are separated from the digital signal. Further, the control and measurement unit 161 also includes memory read/write means 1616 for storing data in and reading data from memory 1617 in the main frame and memory in the sensor. Control and measurement unit 161 per se is known in the art except the module 1613, in which the sensor CTR and the functional light transmission are measured.

[0202]
The inventive part of the main frame resides mainly in compensation unit 162 that includes software to perform the method steps of the invention. The unit 162 comprises of extinction compensation block 1621, path length multiplier compensation block 1622, and block 1623 for calculating subject specific compensation.

[0203]
Extinction compensation block 1621 includes software adapted to create a nominal extinction and to correct said matrix by compensating it for tissue prefilter effects and temperature effect. In other words, tissue prefilter compensator 1601 calculates an amount of tissue prefiltering compensation that is needed and temperature compensator 1602 calculates an amount of the needed temperature compensation. Thus, the blocks 1601 and 1602 are adapted to perform the steps described under subtitles (a) and (b).

[0204]
Path length multiplier compensation block 1622 is adapted to control the value of the transformation used to transform the modulation ratios N_{kl} ^{invivo }to the LambertBeer model N_{kl} ^{LB}. Thus, the block performs operations described previously under subtitle (c). Log (FLT) compensator 1603 fetches the FLTs, that have been calculated in block 1612, and then calculates the ratio of the logarithms for each wavelength pair and for each person of a group. Then this block calculates the relative change of the nominal transformation function. Fcompensator 1605 14 calculates function F_{kl }that can be approximated as the arterial modulation ratio calculated for the venous saturation.

[0205]
The last block of the compensation unit 1621 is a block that calculates subjectspecific compensation based on the results of a invivo measurement and the data stored in memory 1617.

[0206]
The results obtained from block 1623 are fed to hemoglobin fraction calculation unit 163 that is per se known from the prior art.

[0207]
FIG. 17 depicts some units of the sensor. The sensor comprises memory 171 for storing data, particularly the reference data indicating calibration conditions, i.e. data produced by the compensation unit 162. Optionally reference data may also be stored in memory 1617 residing in the main frame of a pulse oximeter. For calculating a correction to the nominal extinction matrix, which is required pursuant to a wavelength shift caused by a temperature change in the emitter chip, the sensor is provided with junction sensitive element 172. The element monitors the temperatures of the pnjunctions of the LEDs and produces signals to be fed to the main frame. Further, the sensor includes light emitters 173 and 174.

[0208]
As can be seen from the above, the method of the invention is based on the DC transmission of light. By means of the DC measurements in the setting up phase, reference data is first created. During subsequent invivo DC measurements, the reference data is then utilized to filter out human variability from the invivo measurement.

[0209]
Although the method in accordance with the invention has been discussed in connection with a four wavelength pulse oximeter, it can also be employed in a basic two wavelength pulse oximeter. However, the method is more beneficial in a multiwavelength pulse oximeter where the number of analytes to be measured is greater than two.

[0210]
In the case of a twowavelength pulse oximeter, the simplest way to apply the compensation is first to formulate the calibration of the twowavelength oximeter as a first step using only one transformation function g^{−1 }(e.g. at wavelengths 660 nm and 900 nm) and a second step using a twotimestwo extinction matrix ε for these wavelengths and for the two analytes RHb and HbO2. The compensation procedures are then identical to the ones presented in the above multiwavelength method. If the calibration of the twowavelength pulse oximeter is done in the normal way using a direct mapping of the invivo measured Rratio (=N_{660900}) to the SpO2 percentage, the compensation steps could for example, be as follows: The wavelength shifts from the nominal LED center wavelength values to a change in the SpO2 value can first be coded. The wavelength shifts are determined for the temperature component as described in the above multiwavelength method and for the tissue component by mapping at the two wavelengths the change in the FLT ratio from its nominal value in the calibration conditions to a change in the SpO2 value from the nominal calibration SpO2. The tissue wavelength shift cannot be estimated as accurately as in the multiwavelength oximeter, but sufficient compensation to the tissue prefilter can still be obtained and the accuracy of the pulse oximeter can be improved. The last compensation step also includes the compensation for the internal tissue variability, which is summed with the prefilter effect.

[0211]
A distinguishing feature of the invention is that compensation is made several times during the invivo measurement. In other words, the pulse oximeter measurement is compensated heartbeatbyheartbeat, i.e. the oxygenation measurement is accurate for each patient and for each time moment in the individual patient. This is not known from the prior art. Further, the same apparatus is used both in the calibration phase, the initial characterization with tissue and the invivo characterization with tissue. The method compensates continuously and dynamically the tissue induced changes and keeps the PO measurement accurate at all times. The apparatus and method is also used in the manner that the same apparatus is suitable both the invivo characterization and the tissue characterization during the calibration. The prior art technology uses different techniques in the initial characterization and invivo characterization.

[0212]
Although the invention has been described above with reference to the examples shown in the appended drawings, it is obvious that the invention is not limited to these, but may be modified by those skilled in the art without departing from the scope and spirit of the invention. For example, instead of transformation, any other quantity by which the pulse oximeter can correct the average calibration known to it can be used to eliminate human variability.

[0213]
The invention has also been described with reference to pulse oximeters for analytes which are in the blood of a subject. The invention, however, can also be applied at different wavelength ranges, e.g. at around 1.5 μm for glucose, at which similar compensation means are called for. The other substances in the tissue modify the effective extinction of the glucose because they alter the path length multiplier at this wavelength. Similarly, the tissue prefilter and temperature effects are taken into account.