CA2373366C - Predictive temperature measurement system - Google Patents
Predictive temperature measurement system Download PDFInfo
- Publication number
- CA2373366C CA2373366C CA002373366A CA2373366A CA2373366C CA 2373366 C CA2373366 C CA 2373366C CA 002373366 A CA002373366 A CA 002373366A CA 2373366 A CA2373366 A CA 2373366A CA 2373366 C CA2373366 C CA 2373366C
- Authority
- CA
- Canada
- Prior art keywords
- temperature
- derivative
- processor
- thermometer
- estimate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000009529 body temperature measurement Methods 0.000 title description 10
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 31
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 2
- 239000000523 sample Substances 0.000 description 42
- 238000005259 measurement Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 13
- 230000006641 stabilisation Effects 0.000 description 11
- 238000011105 stabilization Methods 0.000 description 11
- 238000013459 approach Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 239000011521 glass Substances 0.000 description 6
- 230000036760 body temperature Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000004913 activation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 0 *CCC(CCC1)CC1C1CC(C*)CCC1 Chemical compound *CCC(CCC1)CC1C1CC(C*)CCC1 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000008571 general function Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/42—Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
Abstract
A method and device for rapidly predicting the temperature of an object through a sensor reading of the object's temperature by using a processor comprising finite impulse response (FIR) filters to determine the average value, the first derivative, and the second derivative of the sensed temperature based upon a predetermined number of temperature samples, apply a respective weighting factor to each of the average value, the first derivative, and the second derivative respectively, and calculate a temperature estimate by combining the weighted average value, the first derivative, and second derivative together with an offset factor.
Description
PREDICTIVE TEMPERATURE MEASUREMENT SYSTEM
BACKGROUND
The present invention relates generally to improvements in thermometers and, more particularly, to electronic thermometers for more rapidly obtaining accurate temperature measurements.
It is common practice in the medical field to determine the body temperature of a patient by means of a temperature sensitive device that not only measures the temperature but also displays that temperature. Such temperature measurements are taken routinely in hospitals and in doctors' offices. One such device is a glass bulb thermometer incorporating a heat responsive mercury column that expands and contracts adjacent a calibrated temperature scale.
Typically, the glass thermometer is inserted into the patient, allowed to remain inserted for a sufficient time interval to enable the temperature of the thermometer to stabilize at the body temperature of the patient, and subsequently removed for reading by medical personnel. This time interval is usually on the order of 3 to 8 minutes.
The conventional temperature measurement procedure using a glass bulb thermometer or the like is prone to a number of significant deficiencies.
Temperature measurement is rather slow and, for patients who cannot be relied upon (by virtue of age or infirmity) to properly retain the thermometer for the necessary period of insertion in the body, may necessitate the physical presence of medical personnel during the relatively long measurement cycle, thus diverting their attention from other duties. Furthermore, glass bulb thermometers are not as easy to read and, hence, measurements are more susceptible to human error, particularly when the reading is made under poor lighting conditions or when read by harried personnel.
Various attempts have been made to minimize or eliminate these deficiencies of the glass bulb thermometer by using temperature sensing probes that are designed to operate in conjunction with direct-reading electrical thermometer instrumentation. In one such approach, an electric temperature sensitive device, such as a thermistor, is mounted at the end of a probe and inserted into the patient. The change in voltage or current of the device, depending on the particular implementation, is monitored and when that output signal stabilizes, a temperature is displayed in digital format. This is commonly referred to as the "direct reading" approach and it reduces the possibility of error by misreading the measured temperature. This approach has provided a significant contribution to the technology of temperature measurement.
An inherent characteristic of electronic thermometers is that they do not instantaneously measure the temperature of the site to which they are applied.
It may take a substantial period of time before the temperature sensitive device stabilizes at the temperature of the site and the temperature indicated by the thermometer is representative of the actual temperature of the body or site measured. This lag is caused by the various components of the measurement system that impede heat flow from the surface of the body or site to the temperature sensor. Some of the components are the sensor tip, the tissue of the body, and any hygienic covering applied to the sensor tip to prevent contamination between measurement subjects.
One attempt to shorten the amount of time required to obtain a temperature reading of a subject involves the use of a temperature sensitive electronic probe coupled with prediction or estimation circuitry or programming to provide a digital display of the patient's temperature before the probe has reached equilibrium with the patient. With this approach, assuming the patient's temperature is not significantly changing during the measurement period or cycle, the temperature that will prevail upon thermal stabilization of the electronic thermometer with the patient is predicted from measured temperatures and is displayed before thermal stabilization is attained. Typically, prediction of temperature is performed by monitoring the measured temperature over a period of time as well as the rate of change thereof, and processing these two variables to predict the patient's temperature.
With an electronic thermometer that operates by predicting the final, steady state temperature, an advantage is that the temperature measurement is completed before thermal stabilization is attained, thereby reducing the time required for measurement. This would lessen the risk that the patient would not hold the probe in the correct position for the entire measurement period and requires less time of the attending medical personnel. Another advantage is that because body temperature is dynamic and may significantly change during the five minute interval typically associated with traditional mercury glass thermometer measurements, a rapid determination offers more timely diagnostic information. In addition, the accuracy with which the temperature is predicted improves markedly as the processing and analysis of the data are more accurately performed. This approach has also contributed significantly to the advancement of temperature measurement technology.
Electronic thermometers using predictive-type processing and temperature determination may include a thermistor as a temperature-responsive transducer.
The thermistor approaches its final steady state temperature asymptotically with the last increments of temperature change occurring very slowly, whereas the major portion of the temperature change occurs relatively rapidly. Prior attempts have been made to monitor that initial, more rapid temperature change, extract data from that change, and predict the final temperature at which the thermistor will stabilize and therefore, determine the actual temperature of the tissue that is contacting the thermistor long before the thermistor actually stabilizes at the tissue temperature.
A prior approach used to more rapidly predict the tissue temperature prior to the thermistor reaching equilibrium with that tissue is the sampling of data points of the thermistor early in its response and from those data points, predicting a curve shape of the thermistor's response. From that curve shape, an asymptote of that curve and thus the stabilization, or steady state, temperature can be predicted. To illustrate these concepts through an example of a simpler system, consider the heat transfer physics associated with two bodies of unequal temperature, one having a large thermal mass and the other having a small thermal mass, placed in contact with each other at time t = 0. As time progresses, the temperature of the small thermal mass and the large thermal mass equilibrate to a temperature referred to as the stabilization temperature. The equation describing this process is as follows:
r(t) = TR + (TF - TR) (1 - e -41,0) = TF - (TF - TR)e -cciz:
(Eq. 1) where: T(t) is the temperature of the smaller body as a function of time, TR is the initial temperature of the smaller body, TF is the actual, steady state temperature of the system, t is time, and i is the time constant of the system.
From this relationship, when T is known at two points in time, for example Tl at time tl and T2 at time t2, the stabilization temperature TF can be predicted through application of Equation 2 below.
tZ-tl t2 tl TZ - T1 e T2 e z - T1 e ~
TF t2-tl t2 rl (Eq. 2) 1 _e eT - e`
Further, for a simple first order heat transfer system of the type described by Equation 1, it can be shown that the natural logarithm of the first time derivative of the temperature is a straight line with slope equal to -1/c as follows:
ln dT = K - t (Eq. 3.1) dt and also:
TF = T(t) + tiT l(t) (Eq. 3.2) where:
T l(t) (Eq. 3.3) T 11(t) where: K = a constant dependent upon TR, TF, and i, T' = first derivative T" = second derivative Prior techniques have attempted to apply these simple first order relationships through the use of thermistor time constants established by the 5 thermistor manufacturer. However, these techniques have failed to recognize that the temperature response curve cannot be modeled as first order and is to a great extent affected by factors not reflected by the thermistor's time constant.
When the thermometer is placed in contact with body tissue, such as a person's mouth for example, the response curve depends on the physical placement of the probe in relation to that tissue, on the heat transfer characteristics of the particular tissue, and on the hygienic cover that separates the probe from the tissue.
All of these factors contribute to the heat flow characteristics of the measurement system and they are not represented in the factory-supplied time constant of the thermistor alone. Moreover, the factors described above impede the flow of heat in series and with different resistance characteristics, thus causing an overall time response behavior that is not that of a first order system.
While electronic thermometers and prior predictive techniques have advanced the art of electronic thermometry, a need still exists for an electronic thermometer that can predict a stabilization temperature at a7i early stage of the measurement process and thereby shorten the amount of time taken to obtain a final temperature reading. Additionally, a need exists for a thermometer having an algorithm that can be computed in readily available, relatively simple, relatively inexpensive circuitry or processors. The invention fulfills these needs and others.
SUMMARY OF THE INVENTION
Briefly and in general terms, the present invention is directed to providing an electronic thermometer and a method for determining the temperature of an object or biological subject by predicting the steady state temperature at an early stage of the measurement process. The thermometer and method of the present invention relate certain variables determined from an early portion of the temperature rise curve to the predicted steady state temperature so that the predictive process requires a reduced process of data acquisition and a shortened data processing time while yielding accurate approximations of the stabilization temperature.
Thus, briefly and in general terms, in one aspect of the present invention is directed to a thermometer incorporating a temperature sensor, a processor for predicting an object's temperature based on the average value, slope, and curvature of the initial reading of the object's temperature, and a display for displaying the predicted temperature.
In a more detailed aspect, the processor comprises a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the subject and provide an estimated final temperature signal and a display connected with the processor to receive and display the estimated final temperature signal. In yet further detail, the finite impulse response filter takes a linear combination of a plurality of samples in calculating the estimate of the temperature of the subject. In another aspect, the processor adds an offset coefficient based on ambient temperature to the estimate of the temperature provided by the finite impulse response filter in providing an estimated final temperature signal.
In another aspect, a thermometer for determining the temperature of an object is provided and comprises a sensor that provides a time varying temperature signal in response to sensing the temperature of the object, a processor that receives the temperature signal, samples the temperature signal over a time frame, determines the average value, the first derivative, and the second derivative of the signal over the time frame, combines the average value, first derivative, and second derivative, and thus calculates an estimate of the temperature of the object. In a more detailed aspect, the processor applies a weighting factor to each of the average value, the first derivative, and the second derivative of the signal, and further adds an offset factor selected in accordance with the ambient temperature, to calculate a prediction of the temperature of the object.
BACKGROUND
The present invention relates generally to improvements in thermometers and, more particularly, to electronic thermometers for more rapidly obtaining accurate temperature measurements.
It is common practice in the medical field to determine the body temperature of a patient by means of a temperature sensitive device that not only measures the temperature but also displays that temperature. Such temperature measurements are taken routinely in hospitals and in doctors' offices. One such device is a glass bulb thermometer incorporating a heat responsive mercury column that expands and contracts adjacent a calibrated temperature scale.
Typically, the glass thermometer is inserted into the patient, allowed to remain inserted for a sufficient time interval to enable the temperature of the thermometer to stabilize at the body temperature of the patient, and subsequently removed for reading by medical personnel. This time interval is usually on the order of 3 to 8 minutes.
The conventional temperature measurement procedure using a glass bulb thermometer or the like is prone to a number of significant deficiencies.
Temperature measurement is rather slow and, for patients who cannot be relied upon (by virtue of age or infirmity) to properly retain the thermometer for the necessary period of insertion in the body, may necessitate the physical presence of medical personnel during the relatively long measurement cycle, thus diverting their attention from other duties. Furthermore, glass bulb thermometers are not as easy to read and, hence, measurements are more susceptible to human error, particularly when the reading is made under poor lighting conditions or when read by harried personnel.
Various attempts have been made to minimize or eliminate these deficiencies of the glass bulb thermometer by using temperature sensing probes that are designed to operate in conjunction with direct-reading electrical thermometer instrumentation. In one such approach, an electric temperature sensitive device, such as a thermistor, is mounted at the end of a probe and inserted into the patient. The change in voltage or current of the device, depending on the particular implementation, is monitored and when that output signal stabilizes, a temperature is displayed in digital format. This is commonly referred to as the "direct reading" approach and it reduces the possibility of error by misreading the measured temperature. This approach has provided a significant contribution to the technology of temperature measurement.
An inherent characteristic of electronic thermometers is that they do not instantaneously measure the temperature of the site to which they are applied.
It may take a substantial period of time before the temperature sensitive device stabilizes at the temperature of the site and the temperature indicated by the thermometer is representative of the actual temperature of the body or site measured. This lag is caused by the various components of the measurement system that impede heat flow from the surface of the body or site to the temperature sensor. Some of the components are the sensor tip, the tissue of the body, and any hygienic covering applied to the sensor tip to prevent contamination between measurement subjects.
One attempt to shorten the amount of time required to obtain a temperature reading of a subject involves the use of a temperature sensitive electronic probe coupled with prediction or estimation circuitry or programming to provide a digital display of the patient's temperature before the probe has reached equilibrium with the patient. With this approach, assuming the patient's temperature is not significantly changing during the measurement period or cycle, the temperature that will prevail upon thermal stabilization of the electronic thermometer with the patient is predicted from measured temperatures and is displayed before thermal stabilization is attained. Typically, prediction of temperature is performed by monitoring the measured temperature over a period of time as well as the rate of change thereof, and processing these two variables to predict the patient's temperature.
With an electronic thermometer that operates by predicting the final, steady state temperature, an advantage is that the temperature measurement is completed before thermal stabilization is attained, thereby reducing the time required for measurement. This would lessen the risk that the patient would not hold the probe in the correct position for the entire measurement period and requires less time of the attending medical personnel. Another advantage is that because body temperature is dynamic and may significantly change during the five minute interval typically associated with traditional mercury glass thermometer measurements, a rapid determination offers more timely diagnostic information. In addition, the accuracy with which the temperature is predicted improves markedly as the processing and analysis of the data are more accurately performed. This approach has also contributed significantly to the advancement of temperature measurement technology.
Electronic thermometers using predictive-type processing and temperature determination may include a thermistor as a temperature-responsive transducer.
The thermistor approaches its final steady state temperature asymptotically with the last increments of temperature change occurring very slowly, whereas the major portion of the temperature change occurs relatively rapidly. Prior attempts have been made to monitor that initial, more rapid temperature change, extract data from that change, and predict the final temperature at which the thermistor will stabilize and therefore, determine the actual temperature of the tissue that is contacting the thermistor long before the thermistor actually stabilizes at the tissue temperature.
A prior approach used to more rapidly predict the tissue temperature prior to the thermistor reaching equilibrium with that tissue is the sampling of data points of the thermistor early in its response and from those data points, predicting a curve shape of the thermistor's response. From that curve shape, an asymptote of that curve and thus the stabilization, or steady state, temperature can be predicted. To illustrate these concepts through an example of a simpler system, consider the heat transfer physics associated with two bodies of unequal temperature, one having a large thermal mass and the other having a small thermal mass, placed in contact with each other at time t = 0. As time progresses, the temperature of the small thermal mass and the large thermal mass equilibrate to a temperature referred to as the stabilization temperature. The equation describing this process is as follows:
r(t) = TR + (TF - TR) (1 - e -41,0) = TF - (TF - TR)e -cciz:
(Eq. 1) where: T(t) is the temperature of the smaller body as a function of time, TR is the initial temperature of the smaller body, TF is the actual, steady state temperature of the system, t is time, and i is the time constant of the system.
From this relationship, when T is known at two points in time, for example Tl at time tl and T2 at time t2, the stabilization temperature TF can be predicted through application of Equation 2 below.
tZ-tl t2 tl TZ - T1 e T2 e z - T1 e ~
TF t2-tl t2 rl (Eq. 2) 1 _e eT - e`
Further, for a simple first order heat transfer system of the type described by Equation 1, it can be shown that the natural logarithm of the first time derivative of the temperature is a straight line with slope equal to -1/c as follows:
ln dT = K - t (Eq. 3.1) dt and also:
TF = T(t) + tiT l(t) (Eq. 3.2) where:
T l(t) (Eq. 3.3) T 11(t) where: K = a constant dependent upon TR, TF, and i, T' = first derivative T" = second derivative Prior techniques have attempted to apply these simple first order relationships through the use of thermistor time constants established by the 5 thermistor manufacturer. However, these techniques have failed to recognize that the temperature response curve cannot be modeled as first order and is to a great extent affected by factors not reflected by the thermistor's time constant.
When the thermometer is placed in contact with body tissue, such as a person's mouth for example, the response curve depends on the physical placement of the probe in relation to that tissue, on the heat transfer characteristics of the particular tissue, and on the hygienic cover that separates the probe from the tissue.
All of these factors contribute to the heat flow characteristics of the measurement system and they are not represented in the factory-supplied time constant of the thermistor alone. Moreover, the factors described above impede the flow of heat in series and with different resistance characteristics, thus causing an overall time response behavior that is not that of a first order system.
While electronic thermometers and prior predictive techniques have advanced the art of electronic thermometry, a need still exists for an electronic thermometer that can predict a stabilization temperature at a7i early stage of the measurement process and thereby shorten the amount of time taken to obtain a final temperature reading. Additionally, a need exists for a thermometer having an algorithm that can be computed in readily available, relatively simple, relatively inexpensive circuitry or processors. The invention fulfills these needs and others.
SUMMARY OF THE INVENTION
Briefly and in general terms, the present invention is directed to providing an electronic thermometer and a method for determining the temperature of an object or biological subject by predicting the steady state temperature at an early stage of the measurement process. The thermometer and method of the present invention relate certain variables determined from an early portion of the temperature rise curve to the predicted steady state temperature so that the predictive process requires a reduced process of data acquisition and a shortened data processing time while yielding accurate approximations of the stabilization temperature.
Thus, briefly and in general terms, in one aspect of the present invention is directed to a thermometer incorporating a temperature sensor, a processor for predicting an object's temperature based on the average value, slope, and curvature of the initial reading of the object's temperature, and a display for displaying the predicted temperature.
In a more detailed aspect, the processor comprises a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the subject and provide an estimated final temperature signal and a display connected with the processor to receive and display the estimated final temperature signal. In yet further detail, the finite impulse response filter takes a linear combination of a plurality of samples in calculating the estimate of the temperature of the subject. In another aspect, the processor adds an offset coefficient based on ambient temperature to the estimate of the temperature provided by the finite impulse response filter in providing an estimated final temperature signal.
In another aspect, a thermometer for determining the temperature of an object is provided and comprises a sensor that provides a time varying temperature signal in response to sensing the temperature of the object, a processor that receives the temperature signal, samples the temperature signal over a time frame, determines the average value, the first derivative, and the second derivative of the signal over the time frame, combines the average value, first derivative, and second derivative, and thus calculates an estimate of the temperature of the object. In a more detailed aspect, the processor applies a weighting factor to each of the average value, the first derivative, and the second derivative of the signal, and further adds an offset factor selected in accordance with the ambient temperature, to calculate a prediction of the temperature of the object.
In further detailed aspects, the processor further comprises finite impulse response filters to calculate the average value, slope, and curvature of the temperature data. In a more detailed aspect, the processor continues to saniple the signal and calculate a new prediction for the temperature of the object with each new temperature data value sampled.
In another detailed aspect, the processor monitors a predetertnined number of the last predicted temperatures and calculates the final predicted temperature of the object based on an. average of these last predicted temperatures. In yet another detailed aspect, the processor calculates the final predicted temperature when certain selected conditions have been met. In a still further aspect, the selected conditions include a predetermined time period that must lapse after the sensor has been in contact with the object, predetermined threshold values that the first derivative and the second derivative must reach, and a maximum difference between any two of a predetennined number of the last temperature estimates that must be less than a predetermined threshold.
In accordance with another aspect or illustrative embodiment, there is provided a thermometer for measuring the temperature of an object. The thermometer includes a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object. The thermometer further includes a processor including a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the object and provide an estimated final temperature signal.
The thermometer further includes a display connected with the processor to receive and display the estimated final temperature signal.
In accordance with another illustrative embodiment, there is provided a thermometer for measuring the temperature of an object. The thermometer includes a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object. The thermometer further includes a processor connected so as to sample the temperature signal a predetermined number of times, to determine the average value, the first derivative, and the second derivative of the temperature signal based upon the predetermined number of samples acquired, and to combine the average value, first derivative, and second derivative to calculate an estimate of the temperature of the object and provide an estimated final temperature signal. The thermometer further includes a display connected with the processor to receive and display the estimated final temperature signal. The processor is configured to further apply a separate weighting factor to each of the average value, the first derivative, and the second derivative prior to combining the average value, the first derivative, and the second derivative to calculate the temperature estimate.
The processor is configured to further add an ambient temperature offset coefficient to the sum of the weighted average value, first derivative, and second derivative to calculate the temperature estimate.
In accordance with another illustrative embodiment, there is provided a method for measuring the temperature of an object. The method includes sensing the temperature of lo the object, providing a time varying temperature signal responsive to sensing the temperature of the object, and sampling the temperature signal a predetermined number of times. The method further includes determining the average value of the temperature signal based upon the predetermined number of samples acquired, determining the first derivative of the temperature signal based upon the predetermined number of samples acquired, determining the second derivative of the temperature signal based upon the predetermined number of samples acquired and calculating an estimate of the temperature of the object by combining the average value, the first derivative, and the second derivative. The method further includes providing the sampled temperature signals to at least one finite impulse response (FIR) filter and thereby performing the steps of determining the average value of the temperature signal, determining the first derivative of the temperature signal, and determining the second derivative of the temperature signal by the FIR filter.
These and other features and advantages of illustrative embodiments of the present invention will become apparent from the following more detailed description, when taken in conjunction with the accompanying drawings which illustrate, by way of example, the principles of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a view of an electronic clinical thermometer according to an illustrative embodiment of the present invention;
FIG. 2 is a block diagram of a system in accordance with an illustrative embodiment of the present invention for determining the temperature of an object before final equilibrium of the temperature sensor with the object using the thernlometer shown in FIG. l ;
FIG. 3 presents a graph of a typical thennistor response curve to a temperature that differs from, and is higher than, its own temperature;
FIG. 4 is a diagram illustrating data flow and tasks performed by the system depicted in FIG. 2 ;
FIG. 5 is a diagram illustrating the Initialize System task of FIG. 4;
In another detailed aspect, the processor monitors a predetertnined number of the last predicted temperatures and calculates the final predicted temperature of the object based on an. average of these last predicted temperatures. In yet another detailed aspect, the processor calculates the final predicted temperature when certain selected conditions have been met. In a still further aspect, the selected conditions include a predetermined time period that must lapse after the sensor has been in contact with the object, predetermined threshold values that the first derivative and the second derivative must reach, and a maximum difference between any two of a predetennined number of the last temperature estimates that must be less than a predetermined threshold.
In accordance with another aspect or illustrative embodiment, there is provided a thermometer for measuring the temperature of an object. The thermometer includes a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object. The thermometer further includes a processor including a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the object and provide an estimated final temperature signal.
The thermometer further includes a display connected with the processor to receive and display the estimated final temperature signal.
In accordance with another illustrative embodiment, there is provided a thermometer for measuring the temperature of an object. The thermometer includes a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object. The thermometer further includes a processor connected so as to sample the temperature signal a predetermined number of times, to determine the average value, the first derivative, and the second derivative of the temperature signal based upon the predetermined number of samples acquired, and to combine the average value, first derivative, and second derivative to calculate an estimate of the temperature of the object and provide an estimated final temperature signal. The thermometer further includes a display connected with the processor to receive and display the estimated final temperature signal. The processor is configured to further apply a separate weighting factor to each of the average value, the first derivative, and the second derivative prior to combining the average value, the first derivative, and the second derivative to calculate the temperature estimate.
The processor is configured to further add an ambient temperature offset coefficient to the sum of the weighted average value, first derivative, and second derivative to calculate the temperature estimate.
In accordance with another illustrative embodiment, there is provided a method for measuring the temperature of an object. The method includes sensing the temperature of lo the object, providing a time varying temperature signal responsive to sensing the temperature of the object, and sampling the temperature signal a predetermined number of times. The method further includes determining the average value of the temperature signal based upon the predetermined number of samples acquired, determining the first derivative of the temperature signal based upon the predetermined number of samples acquired, determining the second derivative of the temperature signal based upon the predetermined number of samples acquired and calculating an estimate of the temperature of the object by combining the average value, the first derivative, and the second derivative. The method further includes providing the sampled temperature signals to at least one finite impulse response (FIR) filter and thereby performing the steps of determining the average value of the temperature signal, determining the first derivative of the temperature signal, and determining the second derivative of the temperature signal by the FIR filter.
These and other features and advantages of illustrative embodiments of the present invention will become apparent from the following more detailed description, when taken in conjunction with the accompanying drawings which illustrate, by way of example, the principles of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a view of an electronic clinical thermometer according to an illustrative embodiment of the present invention;
FIG. 2 is a block diagram of a system in accordance with an illustrative embodiment of the present invention for determining the temperature of an object before final equilibrium of the temperature sensor with the object using the thernlometer shown in FIG. l ;
FIG. 3 presents a graph of a typical thennistor response curve to a temperature that differs from, and is higher than, its own temperature;
FIG. 4 is a diagram illustrating data flow and tasks performed by the system depicted in FIG. 2 ;
FIG. 5 is a diagram illustrating the Initialize System task of FIG. 4;
FIG. 6 is a diagram illustrating the Acquire and Filter Temperature Data task of FIG. 4;
FIG. 7 is a diagram illustrating the Calculate Predicted Temperature task of FIG. 4; and FIG. 8 is a block diagram of the processor functions performed by the system in accordance with aspects of the invention as depicted in FIG. 2.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
In the following description, like reference numerals will be used to refer to like or corresponding elements in the different figures of the drawings.
Although temperatures are given in both Fahrenheit and Celsius, the parameters provided are used with only Fahrenheit. Parameters usable for Celsius have not been provided although temperatures have been provided in both Fahrenheit and Celsius.
Referring now to the drawings, and particularly to FIG. 1, there is shown one embodiment of an electronic thermometer incorporating novel features of the present invention. The electronic thermometer 10 contains a probe 12 for sensing the temperature of a selected part of the body, connected by conductors 14 to the main body 16 of the thermometer. As shown in FIG. 1, the probe 12 has been removed from a storage well 17 in the main body. The main body 16 of the thermometer contains the electrical components and power supply of the thermometer, and also has a display 18 for displaying temperature values and error or alarm messages. A second probe 20 is included with the thermometer and is shown in the stored position inserted in a well 19 of the main body 16.
Also shown is a hygienic cover 21 for placement over a probe 12 or 20 before insertion of the probe into the patient.
Referring to FIG. 2, the block diagram generally shows major electronic components of one embodiment of a thermometer 22 in accordance with the present invention. The temperature sensor 24 provides temperature signals in response to the temperature sensed during measurement. In the case where a thermistor is used as the temperature sensor 24, these signals are analog voltages or currents representative of the resistance of the thermistor and therefore representative of the sensed temperature. They are converted into digital form for further processing by an analog-to-digital converter 26. The analog-to-digital converter 26 is connected to a processor 28 that receives the digital temperature signals and processes them to determine the temperature of the subject being measured. A timer 30 provides time signals to the processor 28 used during the processing of the temperature signals, and a memory 32 stores the temperature and time signal data so that the signal data can be analyzed at a subsequent time.
The memory 32 also stores empirically-derived constants used by the processor to calculate the predicted temperature. Once the signals have been processed, the processor 28 provides a signal to the display 34 to display the predicted stabilization temperature. Activating a switch 36 enables the temperature measurement functions of the thermometer 22. This switch is preferably located within the probe storage well such that removal of the probe enables the measurement.
Referring now to FIG. 3, a graph of measured temperature 38 plotted as a function of measurement time for a measurement system is shown. Although the relationship illustrated is similar in form to that specified by Equation 1, the measurement system of the present invention does not exhibit first order heat transfer behavior and therefore the curve 38 differs from the siniple exponential function of Equation 1. As discussed above, the temperature indicated 38 by the thermistor lags the actual temperature TF of the subject being measured. This lag is shown by line 38. It can be seen that as the measurement progresses from a start time, to, the temperature rapidly increases from TR to Tl between times to to tl. The rate of increase in the indicated temperature 38 is reduced between times tl and t2, and the temperature line 38 gradually tends toward the stabilization temperature TF asymptotically as the time increases even more. As discussed above, the present invention is directed to a system capable of analyzing the temperature data gathered during an early period of the measurement, for example, between times tl and t2, and predicting the final temperature TF.
Referring now to FIG. 4, the general functions (tasks) of an embodiment of the system in accordance with aspects of the invention are shown, along with the data that flows among them. A data flow does not imply an activation sequence;
control and activation are not shown on this diagram. The subsequent flow diagrams, FIGS. 5 and 6, illustrate the sequential flow of certain key tasks.
The data flows are labeled with the data that is provided to one task by another.
FIG. 8 5 presents the temperature computation functions performed by the processor 28.
With continued reference to FIG. 4, the Initialize System task 40 is run each time the thermometer is activated. It serves to execute all logic that needs to occur only once per measurement. It activates the Acquire and Filter Temperature Data task 42, which in turn activates the Calculate Predicted Temperature task 44.
FIG. 7 is a diagram illustrating the Calculate Predicted Temperature task of FIG. 4; and FIG. 8 is a block diagram of the processor functions performed by the system in accordance with aspects of the invention as depicted in FIG. 2.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
In the following description, like reference numerals will be used to refer to like or corresponding elements in the different figures of the drawings.
Although temperatures are given in both Fahrenheit and Celsius, the parameters provided are used with only Fahrenheit. Parameters usable for Celsius have not been provided although temperatures have been provided in both Fahrenheit and Celsius.
Referring now to the drawings, and particularly to FIG. 1, there is shown one embodiment of an electronic thermometer incorporating novel features of the present invention. The electronic thermometer 10 contains a probe 12 for sensing the temperature of a selected part of the body, connected by conductors 14 to the main body 16 of the thermometer. As shown in FIG. 1, the probe 12 has been removed from a storage well 17 in the main body. The main body 16 of the thermometer contains the electrical components and power supply of the thermometer, and also has a display 18 for displaying temperature values and error or alarm messages. A second probe 20 is included with the thermometer and is shown in the stored position inserted in a well 19 of the main body 16.
Also shown is a hygienic cover 21 for placement over a probe 12 or 20 before insertion of the probe into the patient.
Referring to FIG. 2, the block diagram generally shows major electronic components of one embodiment of a thermometer 22 in accordance with the present invention. The temperature sensor 24 provides temperature signals in response to the temperature sensed during measurement. In the case where a thermistor is used as the temperature sensor 24, these signals are analog voltages or currents representative of the resistance of the thermistor and therefore representative of the sensed temperature. They are converted into digital form for further processing by an analog-to-digital converter 26. The analog-to-digital converter 26 is connected to a processor 28 that receives the digital temperature signals and processes them to determine the temperature of the subject being measured. A timer 30 provides time signals to the processor 28 used during the processing of the temperature signals, and a memory 32 stores the temperature and time signal data so that the signal data can be analyzed at a subsequent time.
The memory 32 also stores empirically-derived constants used by the processor to calculate the predicted temperature. Once the signals have been processed, the processor 28 provides a signal to the display 34 to display the predicted stabilization temperature. Activating a switch 36 enables the temperature measurement functions of the thermometer 22. This switch is preferably located within the probe storage well such that removal of the probe enables the measurement.
Referring now to FIG. 3, a graph of measured temperature 38 plotted as a function of measurement time for a measurement system is shown. Although the relationship illustrated is similar in form to that specified by Equation 1, the measurement system of the present invention does not exhibit first order heat transfer behavior and therefore the curve 38 differs from the siniple exponential function of Equation 1. As discussed above, the temperature indicated 38 by the thermistor lags the actual temperature TF of the subject being measured. This lag is shown by line 38. It can be seen that as the measurement progresses from a start time, to, the temperature rapidly increases from TR to Tl between times to to tl. The rate of increase in the indicated temperature 38 is reduced between times tl and t2, and the temperature line 38 gradually tends toward the stabilization temperature TF asymptotically as the time increases even more. As discussed above, the present invention is directed to a system capable of analyzing the temperature data gathered during an early period of the measurement, for example, between times tl and t2, and predicting the final temperature TF.
Referring now to FIG. 4, the general functions (tasks) of an embodiment of the system in accordance with aspects of the invention are shown, along with the data that flows among them. A data flow does not imply an activation sequence;
control and activation are not shown on this diagram. The subsequent flow diagrams, FIGS. 5 and 6, illustrate the sequential flow of certain key tasks.
The data flows are labeled with the data that is provided to one task by another.
FIG. 8 5 presents the temperature computation functions performed by the processor 28.
With continued reference to FIG. 4, the Initialize System task 40 is run each time the thermometer is activated. It serves to execute all logic that needs to occur only once per measurement. It activates the Acquire and Filter Temperature Data task 42, which in turn activates the Calculate Predicted Temperature task 44.
10 Once the predicted temperature has been calculated, it is displayed by the Display Predicted Temperature task 46.
FIG. 5 provides a flow diagram for the Initialize System task 40. It is initiated when a probe is removed from the well 60 and initializes, tests, and calibrates the hardware devices 62, initializes the FIR (finite impulse response) filter coefficients 63, and resets the clock counter "t =0" 64 and the running temperature estimates To, Tl, T2, T3, and T4 66 to equal zero. The task 40 then proceeds to the Calculate Offset Coefficient C task 68.
If the probe is not out of the well in step 60, the task continues to increment the probe in the well timer 61. In accordance with this step, the amount of time that the probe 12 is in the well 17 of the body 16 is monitored to determine if the probe is at ambient temperature. If the probe has not been in the well 17 for a certain time period, such as one minute, the measurement system assumes that the probe is not at ambient temperature and a previously-saved ambient temperature is used. If the probe has been in the well 17 for more than a minute, it is considered to be at ambient temperature.
In an alternative embodiment, the Initialize System task 40 may be triggered by other events, such as a rapid rise in the temperature of the probe signifying contact with the patient, or the lapse of a preselected length of time following the removal of the probe from the well, or the activation of the switch 36 (FIG. 2).
FIG. 5 provides a flow diagram for the Initialize System task 40. It is initiated when a probe is removed from the well 60 and initializes, tests, and calibrates the hardware devices 62, initializes the FIR (finite impulse response) filter coefficients 63, and resets the clock counter "t =0" 64 and the running temperature estimates To, Tl, T2, T3, and T4 66 to equal zero. The task 40 then proceeds to the Calculate Offset Coefficient C task 68.
If the probe is not out of the well in step 60, the task continues to increment the probe in the well timer 61. In accordance with this step, the amount of time that the probe 12 is in the well 17 of the body 16 is monitored to determine if the probe is at ambient temperature. If the probe has not been in the well 17 for a certain time period, such as one minute, the measurement system assumes that the probe is not at ambient temperature and a previously-saved ambient temperature is used. If the probe has been in the well 17 for more than a minute, it is considered to be at ambient temperature.
In an alternative embodiment, the Initialize System task 40 may be triggered by other events, such as a rapid rise in the temperature of the probe signifying contact with the patient, or the lapse of a preselected length of time following the removal of the probe from the well, or the activation of the switch 36 (FIG. 2).
In step 68, the offset coefficient C is calculated 68 as shown in Equation 4 below:
C= B. + (B1 X Ta) + (B2 X Ta2 ) (Eq. 4) where: Ta is the ambient temperature, and parameters Bo, B1 and BZ are constant weighting factors empirically derived through the use of well known statistical regression techniques upon a large sample of actual clinical data.
The offset coefficient C is used to factor the ambient temperature into the calculation of the predicted temperature, as detailed elsewhere in the specification. The calculation of the offset factor C, as shown in Equation 4 above, thus relies on the assumption that the initial temperature reading is equal or nearly equal to the ambient temperature Ta in the environment where the thermometer is to be used, and the algorithm determines whether to actually measure T. or to use a previously stored value.
FIG. 6 provides a flow diagram for the Acquire and Filter Temperature Data task 42. The illustrated steps, once initiated, are carried out with precise timing and without interruption to ensure no loss of data and no timing inaccuracy.
The processor 28, timer 30, and analog-to-digital converter 26 (FIG. 2) acquire and filter incoming temperature sensor data to remove line noise and other interference that may affect the quality of the temperature determination.
Various techniques well known in the area of signal processing may be used in this process.
In a preferred embodiment, with each timer interrupt 76, the system updates the clock counter 77 and acquires voltage input 78. This voltage is related to the resistance, and thus the temperature, of the thermistor at the probe tip. A typical mathematical relationship applied to this voltage converts it to a temperature value 82. The system then performs range checks 83 on the temperature for the purpose of determining if the thermometer is broken or poses a possible hazard. Thresholds are set and if the thermometer provides a reading outside these thresholds, the probe is considered to be broken. As an example, if the probe reads less than 32 degrees or greater than 500 degrees F, the probe is considered to be broken (an open circuit or short circuit for example). If the probe reads at a temperature of 115 degrees F, it is considered a hot probe and a hazard alarm is provided.
Temperature data samples are stored 84 in the memory 32. The mathematical relationship between the resistance of the thermistor and the temperature takes into account the thermistor characterization data supplied by its manufacturer, along with the particular circuit details having to do with the interconnection of the thermistor to the analog-to-digital converter, and is well-known to those skilled in the art.
After the temperature value is stored in memory 84, the system determines if tissue contact exists 86 by analyzing temperature data for a temperature above 94 degrees F. Other approaches may be used to determine tissue contact. If contact does not exist, the system waits for the next timer interrupt 76 and repeats the above process until tissue contact is determined 86. At the time that tissue contact is detected, depending on whether enough samples have been acquired or not 87 the system then either awaits the next interrupt 76 to acquire another data sample 78 (if enough samples have not been acquired) or proceeds to calculate the estimate of steady state temperature To. In a preferred embodiment the timer runs at 0.1 second intervals and the system acquires twenty-one samples over a period of time of approximately two seconds prior to commencing to calculate the current temperature estimate To. This number of samples and sampling rate are currently known to provide the best compromise between temperature prediction accuracy and acceptable measurement time. During this period of time, the temperature signals provided by the probe rise along a curve approaching the temperature of the patient asymptotically and the current temperature estimate To is calculated by Equation 5 shown below:
To = AiT+ A2T'+ A3T"+ C (Eq. 5) where: T is the average value of the temperature based upon the twenty-one data samples stored in memory, T' is the first derivative, or slope, of the temperature curve described by the twenty-one samples stored in memory, T" is the second derivative, or curvature, of the temperature curve described by the twenty-one samples stored in memory, Parameters Al, Az and A3 are empirically derived constant weighting factors, and C is the offset coefficient dependent on ambient temperature discussed previously.
Ideally, the number of actual clinical temperature data samples used in deriving parameters Al, A2, A3, Bo, B1, and B2 should be as large as possible and the measured temperatures should be uniformly spread over the entire range of temperatures of potential interest, i.e. the entire range of temperatures that may be encountered by a thermometer according to the invention. However, by necessity the number of temperature data samples obtained was relatively limited and, furthermore, the majority of the temperatures measured were in the normal body temperature range. In particular, the clinical data employed consisted of approximately 240 temperature data samples that ranged from 95.5 to 104 degrees F (35.3 to 40 degrees C) measured at an ambient temperature ranging from 60 to 92 degrees (16 to 33 degrees C).
For this reason, standard regression analysis applied to these data samples produced parameters that tended to predict too low at high actual temperatures and predict too high at low actual temperatures, with a predicted temperature error that exhibited a trend, or relationship, with the stabilization temperature TF
which is essentially equal to the actual body temperature of the patient.
Because the difference between the measured temperature T(t) and the actual temperature TF tends to be a function of the slope T'(t) and the curvature T"(t) of the measured temperature rather than of the measured temperature T(t) itself, parameter Al was artificially constrained to equal 1.0 in order to eliminate the trend exhibited by the error in the computed temperature To. Thus, the remaining five parameters A2, A3, Bo, B1, and B2 were computed by regression analysis of the clinical data to yield the most accurate predictions, in terms of the lowest mean squared prediction error, across the measured temperatures, and in a preferred embodiment their values are:
A1= 1.0 Bo = 15.8877 A2= 7.6136 B1= - 0.3605 (Eq. 5.1) A3 = 8.87 B2 = 0.002123 With continued reference to FIG. 6, once it is determined that tissue contact has been established 86 and a preselected number of samples (e.g.
five) have been acquired 87 since the last temperature estimate has been calculated, the system proceeds to compute the average value T of the temperature 88 and store it in memory 90, compute the first derivative T' of the temperature 92 and store it into memory 94, compute the second derivative T" of the temperature and store it into memory 98, and then finally compute the current temperature estimate To by Equation 5 defined above.
Referring now to FIG. 7, a flow diagram for the Calculate Predicted Temperature task 44 is shown. This task is performed immediately after the Acquire and Filter Temperature Data task 42 which, as shown in FIG. 4, provides filtered and processed data to the Calculate Predicted Temperature task 44.
The first steps entail calculating a prediction 100 and updating the running temperature estimates 110, which are the last four predicted temperatures Tl, T2, T3, and T4 in this case, although other numbers of predicted temperatures may be used in other embodiments. The temperature prediction process of the present invention is performed continuously by the system once the initial twenty-one data points have been acquired (see discussion elsewhere in the specification regarding the use of finite impulse response filters), and a new temperature prediction is computed with every additional five (in this embodiment) data points consequently acquired. Therefore, a new predicted temperature (i.e. the current temperature estimate To) is calculated with every five clock cycles, and thus the oldest of the running temperature estimates T4 (or TN, depending on how many temperature estimates N are employed) must be discarded in favor of To by shifting all values over by one from Tl to T4 and finally updating Tl with the value of To.
The purpose of tracking the last four (or N) temperature estimates is to better determine the first instant when a sufficiently accurate temperature 5 prediction can be posted to the display by determining when an optimal balance between prediction accuracy and elapsed time has been reached. Determining such an optimal balance must take into account various parameters, such as a predetermined minimum wairing time; the number of temperature data samples used to calculate the temperature mean T(t), slope T'(t), and curvature T"(t);
a 10 predetermined maximum slope T'(t); a predetermined maximum and minimum curvature T"(t); a predetermined maximum deviation among the running temperature estimates; and a predetermined number of prior predictions to consider for consistency.
As shown in FIG. 7, in the preferred embodiment of the present invention 15 four such conditions were chosen for determining when the optimal balance between prediction accuracy and elapsed time has been reached. Namely, a minimum waiting time t, a maximum slope T'(t), a maximum and minimum curvature T"(t), and a maximum deviation among the running temperature estimates are determined. Each of these conditions must be met before the Display Predicted Temperature task 46 is activated. Thus, as shown in FIG. 7, in the preferred embodiment at least 5 seconds must have passed 112 since the probe was first activated and placed in contact with the patient's body to ensure that transients have died out. Next, the first derivative T' must have a positive value no larger than 0.25 degrees F(0.14 degrees C) per second 114 to ensure that the temperature is not rising too quickly and that it is therefore "leveling off' and approaching steady state. Third, the second derivative T" must also be between minus and plus 0.05 degrees F (0.028 degrees C) per second per second 116, again to ensure proper convergence toward the steady state temperature.
Fourth, the difference between any two of the running temperature estimates Tl, T2, T3, and T4 ... TNmust be no greater than 0.3 degrees F (0.17 degrees C) in the present embodiment to ensure the accuracy of the final prediction. If any of these conditions is not met, the system initiates the Acquire and Filter Temperature Data task 42 to acquire the next temperature data sample point and recalculate the current temperature estimate To. The preferred values for these four conditions, as presented above, were determined empirically to provide an optimal compromise between the time necessary to compute the prediction and its accuracy.
When all four conditions have been met, the system proceeds to calculate the final temperature estimate Tf 120 by averaging the four running temperature estimates Tl, T2, T3, and T4. Tf is the temperature predicted by the system, and as shown in FIG. 4 the Display Predicted Temperature task 46 can now be activated to display Tf upon the display 18 of the thermometer 10 (FIG. 1). By implementing these multiple conditions, the system of the invention substantially ensures the improved accuracy of the final result by evaluating the uniformity of the temperature data obtained over time and identifying the earliest instant at which the predicted temperature offers a sufficient degree of certainty to be displayed and thus relied upon by the user of the thermometer.
In a preferred embodiment the processor 28 (FIG. 2) is provided with finite impulse response ("FIR") digital filters to perform the computations necessary to determine T, T' and T". FIR filters are typically configured to operate upon a string of K numerical values by multiplying each value with a predetermined coefficient hk. In a preferred embodiment, the filter is implemented by a computer program and the CPU 28 memory 32 is used for storage. In another embodiment, the filter may comprise K registers connected in a linear configuration and, as each new value is acquired and fed to the filter, the new value is stored in the first register and all previously acquired values are shifted one register over such that eventually the first acquired value is stored in the last register. Each register k has a coefficient hk associated with it, and as a new value is shifted into each register, that value is multiplied with the coefficient hk associated with the particular register. In the present implementation these steps are performed by software, and the processor memory 32 is utilized to store the register values.
C= B. + (B1 X Ta) + (B2 X Ta2 ) (Eq. 4) where: Ta is the ambient temperature, and parameters Bo, B1 and BZ are constant weighting factors empirically derived through the use of well known statistical regression techniques upon a large sample of actual clinical data.
The offset coefficient C is used to factor the ambient temperature into the calculation of the predicted temperature, as detailed elsewhere in the specification. The calculation of the offset factor C, as shown in Equation 4 above, thus relies on the assumption that the initial temperature reading is equal or nearly equal to the ambient temperature Ta in the environment where the thermometer is to be used, and the algorithm determines whether to actually measure T. or to use a previously stored value.
FIG. 6 provides a flow diagram for the Acquire and Filter Temperature Data task 42. The illustrated steps, once initiated, are carried out with precise timing and without interruption to ensure no loss of data and no timing inaccuracy.
The processor 28, timer 30, and analog-to-digital converter 26 (FIG. 2) acquire and filter incoming temperature sensor data to remove line noise and other interference that may affect the quality of the temperature determination.
Various techniques well known in the area of signal processing may be used in this process.
In a preferred embodiment, with each timer interrupt 76, the system updates the clock counter 77 and acquires voltage input 78. This voltage is related to the resistance, and thus the temperature, of the thermistor at the probe tip. A typical mathematical relationship applied to this voltage converts it to a temperature value 82. The system then performs range checks 83 on the temperature for the purpose of determining if the thermometer is broken or poses a possible hazard. Thresholds are set and if the thermometer provides a reading outside these thresholds, the probe is considered to be broken. As an example, if the probe reads less than 32 degrees or greater than 500 degrees F, the probe is considered to be broken (an open circuit or short circuit for example). If the probe reads at a temperature of 115 degrees F, it is considered a hot probe and a hazard alarm is provided.
Temperature data samples are stored 84 in the memory 32. The mathematical relationship between the resistance of the thermistor and the temperature takes into account the thermistor characterization data supplied by its manufacturer, along with the particular circuit details having to do with the interconnection of the thermistor to the analog-to-digital converter, and is well-known to those skilled in the art.
After the temperature value is stored in memory 84, the system determines if tissue contact exists 86 by analyzing temperature data for a temperature above 94 degrees F. Other approaches may be used to determine tissue contact. If contact does not exist, the system waits for the next timer interrupt 76 and repeats the above process until tissue contact is determined 86. At the time that tissue contact is detected, depending on whether enough samples have been acquired or not 87 the system then either awaits the next interrupt 76 to acquire another data sample 78 (if enough samples have not been acquired) or proceeds to calculate the estimate of steady state temperature To. In a preferred embodiment the timer runs at 0.1 second intervals and the system acquires twenty-one samples over a period of time of approximately two seconds prior to commencing to calculate the current temperature estimate To. This number of samples and sampling rate are currently known to provide the best compromise between temperature prediction accuracy and acceptable measurement time. During this period of time, the temperature signals provided by the probe rise along a curve approaching the temperature of the patient asymptotically and the current temperature estimate To is calculated by Equation 5 shown below:
To = AiT+ A2T'+ A3T"+ C (Eq. 5) where: T is the average value of the temperature based upon the twenty-one data samples stored in memory, T' is the first derivative, or slope, of the temperature curve described by the twenty-one samples stored in memory, T" is the second derivative, or curvature, of the temperature curve described by the twenty-one samples stored in memory, Parameters Al, Az and A3 are empirically derived constant weighting factors, and C is the offset coefficient dependent on ambient temperature discussed previously.
Ideally, the number of actual clinical temperature data samples used in deriving parameters Al, A2, A3, Bo, B1, and B2 should be as large as possible and the measured temperatures should be uniformly spread over the entire range of temperatures of potential interest, i.e. the entire range of temperatures that may be encountered by a thermometer according to the invention. However, by necessity the number of temperature data samples obtained was relatively limited and, furthermore, the majority of the temperatures measured were in the normal body temperature range. In particular, the clinical data employed consisted of approximately 240 temperature data samples that ranged from 95.5 to 104 degrees F (35.3 to 40 degrees C) measured at an ambient temperature ranging from 60 to 92 degrees (16 to 33 degrees C).
For this reason, standard regression analysis applied to these data samples produced parameters that tended to predict too low at high actual temperatures and predict too high at low actual temperatures, with a predicted temperature error that exhibited a trend, or relationship, with the stabilization temperature TF
which is essentially equal to the actual body temperature of the patient.
Because the difference between the measured temperature T(t) and the actual temperature TF tends to be a function of the slope T'(t) and the curvature T"(t) of the measured temperature rather than of the measured temperature T(t) itself, parameter Al was artificially constrained to equal 1.0 in order to eliminate the trend exhibited by the error in the computed temperature To. Thus, the remaining five parameters A2, A3, Bo, B1, and B2 were computed by regression analysis of the clinical data to yield the most accurate predictions, in terms of the lowest mean squared prediction error, across the measured temperatures, and in a preferred embodiment their values are:
A1= 1.0 Bo = 15.8877 A2= 7.6136 B1= - 0.3605 (Eq. 5.1) A3 = 8.87 B2 = 0.002123 With continued reference to FIG. 6, once it is determined that tissue contact has been established 86 and a preselected number of samples (e.g.
five) have been acquired 87 since the last temperature estimate has been calculated, the system proceeds to compute the average value T of the temperature 88 and store it in memory 90, compute the first derivative T' of the temperature 92 and store it into memory 94, compute the second derivative T" of the temperature and store it into memory 98, and then finally compute the current temperature estimate To by Equation 5 defined above.
Referring now to FIG. 7, a flow diagram for the Calculate Predicted Temperature task 44 is shown. This task is performed immediately after the Acquire and Filter Temperature Data task 42 which, as shown in FIG. 4, provides filtered and processed data to the Calculate Predicted Temperature task 44.
The first steps entail calculating a prediction 100 and updating the running temperature estimates 110, which are the last four predicted temperatures Tl, T2, T3, and T4 in this case, although other numbers of predicted temperatures may be used in other embodiments. The temperature prediction process of the present invention is performed continuously by the system once the initial twenty-one data points have been acquired (see discussion elsewhere in the specification regarding the use of finite impulse response filters), and a new temperature prediction is computed with every additional five (in this embodiment) data points consequently acquired. Therefore, a new predicted temperature (i.e. the current temperature estimate To) is calculated with every five clock cycles, and thus the oldest of the running temperature estimates T4 (or TN, depending on how many temperature estimates N are employed) must be discarded in favor of To by shifting all values over by one from Tl to T4 and finally updating Tl with the value of To.
The purpose of tracking the last four (or N) temperature estimates is to better determine the first instant when a sufficiently accurate temperature 5 prediction can be posted to the display by determining when an optimal balance between prediction accuracy and elapsed time has been reached. Determining such an optimal balance must take into account various parameters, such as a predetermined minimum wairing time; the number of temperature data samples used to calculate the temperature mean T(t), slope T'(t), and curvature T"(t);
a 10 predetermined maximum slope T'(t); a predetermined maximum and minimum curvature T"(t); a predetermined maximum deviation among the running temperature estimates; and a predetermined number of prior predictions to consider for consistency.
As shown in FIG. 7, in the preferred embodiment of the present invention 15 four such conditions were chosen for determining when the optimal balance between prediction accuracy and elapsed time has been reached. Namely, a minimum waiting time t, a maximum slope T'(t), a maximum and minimum curvature T"(t), and a maximum deviation among the running temperature estimates are determined. Each of these conditions must be met before the Display Predicted Temperature task 46 is activated. Thus, as shown in FIG. 7, in the preferred embodiment at least 5 seconds must have passed 112 since the probe was first activated and placed in contact with the patient's body to ensure that transients have died out. Next, the first derivative T' must have a positive value no larger than 0.25 degrees F(0.14 degrees C) per second 114 to ensure that the temperature is not rising too quickly and that it is therefore "leveling off' and approaching steady state. Third, the second derivative T" must also be between minus and plus 0.05 degrees F (0.028 degrees C) per second per second 116, again to ensure proper convergence toward the steady state temperature.
Fourth, the difference between any two of the running temperature estimates Tl, T2, T3, and T4 ... TNmust be no greater than 0.3 degrees F (0.17 degrees C) in the present embodiment to ensure the accuracy of the final prediction. If any of these conditions is not met, the system initiates the Acquire and Filter Temperature Data task 42 to acquire the next temperature data sample point and recalculate the current temperature estimate To. The preferred values for these four conditions, as presented above, were determined empirically to provide an optimal compromise between the time necessary to compute the prediction and its accuracy.
When all four conditions have been met, the system proceeds to calculate the final temperature estimate Tf 120 by averaging the four running temperature estimates Tl, T2, T3, and T4. Tf is the temperature predicted by the system, and as shown in FIG. 4 the Display Predicted Temperature task 46 can now be activated to display Tf upon the display 18 of the thermometer 10 (FIG. 1). By implementing these multiple conditions, the system of the invention substantially ensures the improved accuracy of the final result by evaluating the uniformity of the temperature data obtained over time and identifying the earliest instant at which the predicted temperature offers a sufficient degree of certainty to be displayed and thus relied upon by the user of the thermometer.
In a preferred embodiment the processor 28 (FIG. 2) is provided with finite impulse response ("FIR") digital filters to perform the computations necessary to determine T, T' and T". FIR filters are typically configured to operate upon a string of K numerical values by multiplying each value with a predetermined coefficient hk. In a preferred embodiment, the filter is implemented by a computer program and the CPU 28 memory 32 is used for storage. In another embodiment, the filter may comprise K registers connected in a linear configuration and, as each new value is acquired and fed to the filter, the new value is stored in the first register and all previously acquired values are shifted one register over such that eventually the first acquired value is stored in the last register. Each register k has a coefficient hk associated with it, and as a new value is shifted into each register, that value is multiplied with the coefficient hk associated with the particular register. In the present implementation these steps are performed by software, and the processor memory 32 is utilized to store the register values.
The output of a FIR filter is equal to the sum of the values in all of the filter's registers multiplied by the respective coefficients hk associated with each register. Thus, generally speaking, given an input sequence x(n), n = 0, 1, 2 ... K, a FIR filter with coefficients ho, hl, ...hK will provide a filtered output signal y(n) as follows:
y(n) y k o hk x x(n - k) (Eq. 6.1) wherein: hk (k = o... x-1) = set of coefficients for the particular FIR
filter, K = number of coefficients for the FIR filter, k = index x(n) corresponds to the value of the last temperature data sample acquired and fed to the filter.
Thus, through the use of the proper coefficients, FIR filters can be programmed to extract the average value T, the slope T', and the curvature T" of the sample data curve from the temperature data sample points acquired and stored by the system.
These coefficients can be readily derived for an input sequence of twenty-one data points using known mathematical methods and, for a system of twenty-one samples (n = -10, -9 ... 10, where n=10 is the most recent saznple) acquired at 0.1 second intervals, the coefficients for the FIR filters are as shown in the following equations:
T n=-10 21 x x(n) (Eq. 6.2) T,, _lio n x x(n) (Eq.6.3) ,I=-io 0.1 x 770 2 x n2 - 770 _ l0 21 ~'ii x x(n) (Eq. 6.4) n=-10 7702 0.12 50666 -where:
>n2 = 770 n_ _10 (Eq. 6.5) >=' =50666 _10 With reference to FIG. 8, a block diagram depicting the operation of the 5 processor, including the FIR filters, is shown. The last temperature data sample acquired x(n) is provided by the analog/digital converter (not shown) to the processor 28, where it is simultaneously fed to FIR circuits FIRl 202, FIR2 204, and FIR3 206, which are configured to calculate the average value T, slope T', and curvature T", respectively, of the last twenty-one temperature data points provided 10 by the analog/digital converter, including x(n). The average value T, slope T', and curvature T" are each multiplied by the respective weighting factor Al 212, A2 214, and A3 216, as per Equation 5 above, and then summed together 218.
Finally, the offset factor C, which as detailed previously is a function of the ambient temperature Ta, is added 220 to the sum of the weighted average value T, slope T', and curvature T" to calculate the current temperature estimate To.
It must be noted that the curvature T"(t) may also be computed by calculating the slope of the slope T'(t) by concatenating two slope calculating filters FIR2. Such an approach, however, was not selected for use in the preferred embodiment of the invention because it would require calculating the slope T(t) prior to, rather than simultaneously with, the curvature T" (t), and would thus cause delays in computing the final temperature.
By employing FIR filters to derive the average value, slope, and curvature of the temperature data, the system of the invention can utilize a rather sophisticated algorithm for predicting temperature with readily available, relatively inexpensive mathematical processors, such as commonly available eight-bit processors. As one example, the eight-bit processor having a part no. of UP078064 from NEC may be used. Lastly, the algorithm of the invention was fine tuned by applying it to actual data to empirically derive weighting factors that provide the most accurate results over the widest range of final, steady state temperatures.
In an alternative embodiment of the present invention, a single FIR filter may be programmed to extract the weighted sum of the temperature mean T(t), slope T'(t), and curvature T"(t). However, the use of a single FIR filter would not allow the individual extraction of the mean T(t), slope T(t), and curvature T"(t) from the sampled data and therefore these parameters could not be individually monitored to determined when the optimal balance between prediction accuracy and elapsed time has been reached, as previously discussed. Furthermore, the use of individual FIR filters allows for the individual adjustment of each of these parameters such as, for example, adjusting the size of each FIR filter to obtain a particular amount of smoothing for each parameter.
While one form of the invention has been illustrated and described, it will be apparent that further modifications and improvements may additionally be made to the device and method disclosed herein without departing from the scope of the invention. Accordingly, it is not intended that the invention be limited, except as by the appended claims.
y(n) y k o hk x x(n - k) (Eq. 6.1) wherein: hk (k = o... x-1) = set of coefficients for the particular FIR
filter, K = number of coefficients for the FIR filter, k = index x(n) corresponds to the value of the last temperature data sample acquired and fed to the filter.
Thus, through the use of the proper coefficients, FIR filters can be programmed to extract the average value T, the slope T', and the curvature T" of the sample data curve from the temperature data sample points acquired and stored by the system.
These coefficients can be readily derived for an input sequence of twenty-one data points using known mathematical methods and, for a system of twenty-one samples (n = -10, -9 ... 10, where n=10 is the most recent saznple) acquired at 0.1 second intervals, the coefficients for the FIR filters are as shown in the following equations:
T n=-10 21 x x(n) (Eq. 6.2) T,, _lio n x x(n) (Eq.6.3) ,I=-io 0.1 x 770 2 x n2 - 770 _ l0 21 ~'ii x x(n) (Eq. 6.4) n=-10 7702 0.12 50666 -where:
>n2 = 770 n_ _10 (Eq. 6.5) >=' =50666 _10 With reference to FIG. 8, a block diagram depicting the operation of the 5 processor, including the FIR filters, is shown. The last temperature data sample acquired x(n) is provided by the analog/digital converter (not shown) to the processor 28, where it is simultaneously fed to FIR circuits FIRl 202, FIR2 204, and FIR3 206, which are configured to calculate the average value T, slope T', and curvature T", respectively, of the last twenty-one temperature data points provided 10 by the analog/digital converter, including x(n). The average value T, slope T', and curvature T" are each multiplied by the respective weighting factor Al 212, A2 214, and A3 216, as per Equation 5 above, and then summed together 218.
Finally, the offset factor C, which as detailed previously is a function of the ambient temperature Ta, is added 220 to the sum of the weighted average value T, slope T', and curvature T" to calculate the current temperature estimate To.
It must be noted that the curvature T"(t) may also be computed by calculating the slope of the slope T'(t) by concatenating two slope calculating filters FIR2. Such an approach, however, was not selected for use in the preferred embodiment of the invention because it would require calculating the slope T(t) prior to, rather than simultaneously with, the curvature T" (t), and would thus cause delays in computing the final temperature.
By employing FIR filters to derive the average value, slope, and curvature of the temperature data, the system of the invention can utilize a rather sophisticated algorithm for predicting temperature with readily available, relatively inexpensive mathematical processors, such as commonly available eight-bit processors. As one example, the eight-bit processor having a part no. of UP078064 from NEC may be used. Lastly, the algorithm of the invention was fine tuned by applying it to actual data to empirically derive weighting factors that provide the most accurate results over the widest range of final, steady state temperatures.
In an alternative embodiment of the present invention, a single FIR filter may be programmed to extract the weighted sum of the temperature mean T(t), slope T'(t), and curvature T"(t). However, the use of a single FIR filter would not allow the individual extraction of the mean T(t), slope T(t), and curvature T"(t) from the sampled data and therefore these parameters could not be individually monitored to determined when the optimal balance between prediction accuracy and elapsed time has been reached, as previously discussed. Furthermore, the use of individual FIR filters allows for the individual adjustment of each of these parameters such as, for example, adjusting the size of each FIR filter to obtain a particular amount of smoothing for each parameter.
While one form of the invention has been illustrated and described, it will be apparent that further modifications and improvements may additionally be made to the device and method disclosed herein without departing from the scope of the invention. Accordingly, it is not intended that the invention be limited, except as by the appended claims.
Claims (31)
PROPERTY OR PRIVILEGE IS CLAIMED ARE DEFINED AS FOLLOWS:
1. A thermometer for measuring the temperature of an object, comprising:
a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object;
a processor comprising a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the object and provide an estimated final temperature signal; and a display connected with the processor to receive and display the estimated final temperature signal.
a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object;
a processor comprising a finite impulse response filter connected so as to sample the temperature signal a plurality of times to calculate an estimate of the temperature of the object and provide an estimated final temperature signal; and a display connected with the processor to receive and display the estimated final temperature signal.
2. The thermometer of claim 1 wherein the finite impulse response filter takes a linear combination of a plurality of samples in calculating the estimate of the temperature of the object.
3. The thermometer of claim 1 or claim 2 wherein the processor adds an offset coefficient based on ambient temperature to the estimate of the temperature provided by the finite impulse response filter in providing an estimated final temperature signal.
4. The thermometer of claim 3 wherein:
the sensor further provides an initial signal corresponding to the ambient temperature;
and the processor further determines the square of the ambient temperature, applies a separate weighting factor to each of the ambient temperature and the square of the ambient temperature, and combines the weighted value of the ambient temperature with the weighted square of the value of the ambient temperature to calculate the ambient temperature offset coefficient.
the sensor further provides an initial signal corresponding to the ambient temperature;
and the processor further determines the square of the ambient temperature, applies a separate weighting factor to each of the ambient temperature and the square of the ambient temperature, and combines the weighted value of the ambient temperature with the weighted square of the value of the ambient temperature to calculate the ambient temperature offset coefficient.
5. The thermometer of any one of claim 1 to claim 4 wherein the processor is further connected to determine the average value, the first derivative, and the second derivative of the temperature signal based upon the predetermined number of samples acquired and the processor is also configured to combine the average value, first derivative, and second derivative to calculate the estimate of the temperature of the object.
6. The thermometer of any one of claim 1 to claim 5 wherein the processor further applies a separate weighting factor to each of the average value, the first derivative, and the second derivative prior to combining the average value, the first derivative, and the second derivative to calculate the temperature estimate.
7. The thermometer of any one of claim 1 to claim 6 wherein the predetermined number of samples are acquired over a predetermined time frame.
8. The thermometer of claim 7 wherein the processor samples the temperature signal around twenty-one times over a time frame of approximately two seconds to determine the average value, the first derivative, and the second derivative of the signal.
9. The thermometer of any one of claim 1 to claim 8 wherein the processor continuously samples the signal after acquiring an initial predetermined number of samples and continuously recalculates the temperature estimate based on the last predetermined number of samples acquired as each new sample is acquired.
10. The thermometer of claim 9 wherein the thermometer displays on the display a final temperature estimate calculated by the processor after the sensor has been in contact with the object for a predetermined length of time.
11. The thermometer of claim 9 wherein the thermometer displays on the display a final temperature estimate calculated by the processor only if the last value of the first derivative calculated by the processor lies in a predetermined range of values.
12. The thermometer of claim 9 wherein the thermometer displays on the display a final temperature estimate calculated by the processor only if the last value of the second derivative calculated by the processor lies in a predetermined range of values.
13. The thermometer of claim 8 wherein the thermometer displays on the display a final temperature estimate only if the maximum difference between any two of a predetermined number of the last temperature estimates calculated by the processor is less than a predetermined value.
14. The thermometer of claim 9 wherein the final temperature estimate displayed by the thermometer is equal to the average of a selected number of the last temperature estimates calculated by the processor.
15. The thermometer of any one of claim 1 to claim 14 wherein the processor comprises:
a first finite impulse response filter configured to calculate the average value based on the last predetermined number of samples acquired by the processor;
a second finite impulse response filter configured to calculate the first derivative based on the last predetermined number of samples acquired by the processor;
and a third finite impulse response filter configured to calculate the second derivative based on the last predetermined number of samples acquired by the processor.
a first finite impulse response filter configured to calculate the average value based on the last predetermined number of samples acquired by the processor;
a second finite impulse response filter configured to calculate the first derivative based on the last predetermined number of samples acquired by the processor;
and a third finite impulse response filter configured to calculate the second derivative based on the last predetermined number of samples acquired by the processor.
16. A thermometer for measuring the temperature of an object, comprising:
a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object;
a processor connected so as to sample the temperature signal a predetermined number of times, to determine the average value, the first derivative, and the second derivative of the temperature signal based upon the predetermined number of samples acquired, and to combine the average value, first derivative, and second derivative to calculate an estimate of the temperature of the object and provide an estimated final temperature signal; and a display connected with the processor to receive and display the estimated final temperature signal;
wherein the processor is configured to further apply a separate weighting factor to each of the average value, the first derivative, and the second derivative prior to combining the average value, the first derivative, and the second derivative to calculate the temperature estimate; and wherein the processor is configured to further add an ambient temperature offset coefficient to the sum of the weighted average value, first derivative, and second derivative to calculate the temperature estimate.
a sensor which, when positioned to sense the temperature of the object, provides a time varying temperature signal in response to the temperature of the object;
a processor connected so as to sample the temperature signal a predetermined number of times, to determine the average value, the first derivative, and the second derivative of the temperature signal based upon the predetermined number of samples acquired, and to combine the average value, first derivative, and second derivative to calculate an estimate of the temperature of the object and provide an estimated final temperature signal; and a display connected with the processor to receive and display the estimated final temperature signal;
wherein the processor is configured to further apply a separate weighting factor to each of the average value, the first derivative, and the second derivative prior to combining the average value, the first derivative, and the second derivative to calculate the temperature estimate; and wherein the processor is configured to further add an ambient temperature offset coefficient to the sum of the weighted average value, first derivative, and second derivative to calculate the temperature estimate.
17. The thermometer of claim 16 wherein:
the sensor further provides an initial signal corresponding to the ambient temperature;
and the processor further determines the square of the ambient temperature, applies a separate weighting factor to each of the ambient temperature and the square of the ambient temperature, and combines the weighted value of the ambient temperature with the weighted square of the value of the ambient temperature to calculate the ambient temperature offset coefficient.
the sensor further provides an initial signal corresponding to the ambient temperature;
and the processor further determines the square of the ambient temperature, applies a separate weighting factor to each of the ambient temperature and the square of the ambient temperature, and combines the weighted value of the ambient temperature with the weighted square of the value of the ambient temperature to calculate the ambient temperature offset coefficient.
18. A method for measuring the temperature of an object, comprising:
sensing the temperature of the object;
providing a time varying temperature signal responsive to sensing the temperature of the object;
sampling the temperature signal a predetermined number of times;
determining the average value of the temperature signal based upon the predetermined number of samples acquired;
determining the first derivative of the temperature signal based upon the predetermined number of samples acquired;
determining the second derivative of the temperature signal based upon the predetermined number of samples acquired;
calculating an estimate of the temperature of the object by combining the average value, the first derivative, and the second derivative; and providing the sampled temperature signals to at least one finite impulse response (FIR) filter and thereby performing the steps of determining the average value of the temperature signal, determining the first derivative of the temperature signal, and determining the second derivative of the temperature signal by the FIR filter.
sensing the temperature of the object;
providing a time varying temperature signal responsive to sensing the temperature of the object;
sampling the temperature signal a predetermined number of times;
determining the average value of the temperature signal based upon the predetermined number of samples acquired;
determining the first derivative of the temperature signal based upon the predetermined number of samples acquired;
determining the second derivative of the temperature signal based upon the predetermined number of samples acquired;
calculating an estimate of the temperature of the object by combining the average value, the first derivative, and the second derivative; and providing the sampled temperature signals to at least one finite impulse response (FIR) filter and thereby performing the steps of determining the average value of the temperature signal, determining the first derivative of the temperature signal, and determining the second derivative of the temperature signal by the FIR filter.
19. The method of claim 18 comprising the further steps of:
applying a first weighting factor to the average value;
applying a second weighting factor to the first derivative; and applying a third weighting factor to the second derivative.
applying a first weighting factor to the average value;
applying a second weighting factor to the first derivative; and applying a third weighting factor to the second derivative.
20. The method of claim 18 or claim 19 wherein the step of sampling the temperature signal a predetermined number of times comprises sampling the temperature signal around twenty-one times over a time frame of approximately two seconds.
21. The method of any one of claim 18 to claim 20 comprising the further steps of:
repeatedly sampling the temperature signal; and repeatedly recalculating the temperature estimate based on the last predetermined number of samples obtained as each new sample is acquired.
repeatedly sampling the temperature signal; and repeatedly recalculating the temperature estimate based on the last predetermined number of samples obtained as each new sample is acquired.
22. The method of any one of claim 18 to claim 21 comprising the further step of displaying the temperature estimate on a display.
23. The method of claim 22 wherein the step of displaying the temperature estimate is performed after the temperature of the object has been sensed for a predetermined length of time.
24. The method of claim 22 wherein the step of displaying the temperature estimate is performed only when the last value of the first derivative calculated by the processor lies in a predetermined range of values.
25. The method of claim 22 wherein the step of displaying the temperature estimate is performed only when the last value of the second derivative calculated by the processor lies in a predetermined range of values.
26. The method of claim 22 wherein the step of displaying the temperature estimate is performed only when the maximum difference between selected temperature estimates calculated by the processor is less than a predetermined value.
27. The method of claim 22 wherein the step of calculating an estimate of the temperature of the object comprises the further step of:
averaging a predetermined number of the temperature estimates calculated by the processor to calculate a final temperature estimate; and wherein the step of displaying the temperature estimate is comprised of displaying the final temperature estimate.
averaging a predetermined number of the temperature estimates calculated by the processor to calculate a final temperature estimate; and wherein the step of displaying the temperature estimate is comprised of displaying the final temperature estimate.
28. The method of any one of claim 18 to claim 27 wherein the step of providing the sampled temperature signal to at least one FIR filter comprises the steps of:
providing the sampled temperature signals to a first FIR for calculating the average value based on a predetermined number of sample values acquired;
providing the sampled temperature signals to a second finite impulse response filter for calculating the first derivative based on a predetermined number of sample values acquired; and providing the sampled temperature signals to a third finite impulse response filter for calculating the second derivative based on a predetermined number of sample values acquired.
providing the sampled temperature signals to a first FIR for calculating the average value based on a predetermined number of sample values acquired;
providing the sampled temperature signals to a second finite impulse response filter for calculating the first derivative based on a predetermined number of sample values acquired; and providing the sampled temperature signals to a third finite impulse response filter for calculating the second derivative based on a predetermined number of sample values acquired.
29. The method of any one of claim 18 to claim 28 wherein:
the step of providing a FIR filter comprises applying weights to the average value, the first derivative, and the second derivative based on a predetermined number of sample values acquired.
the step of providing a FIR filter comprises applying weights to the average value, the first derivative, and the second derivative based on a predetermined number of sample values acquired.
30. The method of claim 29 further comprising the step of combining the weighed average value, the weighted first derivative, and the weighted second derivative to determine the estimated temperature of the object.
31. The method of any one of claim 18 to claim 30 comprising the further steps of:
sensing the ambient temperature;
determining the square of the value of the ambient temperature;
applying a fourth weighting factor to the ambient temperature; and applying a fifth weighting factor to the square of the ambient temperature;
wherein the step of calculating the estimate of the temperature of the object comprises combining the weighted average value, the weighted first derivative, and the weighted second derivative with the weighted ambient temperature and the weighted square of the ambient temperature.
sensing the ambient temperature;
determining the square of the value of the ambient temperature;
applying a fourth weighting factor to the ambient temperature; and applying a fifth weighting factor to the square of the ambient temperature;
wherein the step of calculating the estimate of the temperature of the object comprises combining the weighted average value, the weighted first derivative, and the weighted second derivative with the weighted ambient temperature and the weighted square of the ambient temperature.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/313,911 | 1999-05-18 | ||
| US09/313,911 US6270252B1 (en) | 1999-05-18 | 1999-05-18 | Predictive temperature measurement system |
| PCT/US2000/013496 WO2000070316A1 (en) | 1999-05-18 | 2000-05-17 | Predictive temperature measurement system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2373366A1 CA2373366A1 (en) | 2000-11-23 |
| CA2373366C true CA2373366C (en) | 2009-12-01 |
Family
ID=23217707
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002373366A Expired - Fee Related CA2373366C (en) | 1999-05-18 | 2000-05-17 | Predictive temperature measurement system |
Country Status (10)
| Country | Link |
|---|---|
| US (2) | US6270252B1 (en) |
| EP (1) | EP1183508B1 (en) |
| JP (3) | JP2002544512A (en) |
| AT (1) | ATE345491T1 (en) |
| AU (1) | AU781183B2 (en) |
| CA (1) | CA2373366C (en) |
| DE (1) | DE60031847T2 (en) |
| ES (1) | ES2275515T3 (en) |
| HK (1) | HK1044371B (en) |
| WO (1) | WO2000070316A1 (en) |
Families Citing this family (74)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6056435A (en) * | 1997-06-24 | 2000-05-02 | Exergen Corporation | Ambient and perfusion normalized temperature detector |
| ATE348565T1 (en) * | 1998-10-20 | 2007-01-15 | Omron Healthcare Co Ltd | INFRARED THERMOMETHER |
| US7140766B2 (en) * | 1999-08-04 | 2006-11-28 | Given Imaging Ltd. | Device, system and method for temperature sensing in an in-vivo device |
| TW393572B (en) * | 1999-11-30 | 2000-06-11 | Taidoc Technology Co Ltd | Fast clinical thermometer and its forecast method |
| EP1249691A1 (en) * | 2001-04-11 | 2002-10-16 | Omron Corporation | Electronic clinical thermometer |
| JP3900865B2 (en) * | 2001-06-04 | 2007-04-04 | オムロンヘルスケア株式会社 | Infrared thermometer, infrared thermometer temperature state estimation method, information notification method, and measurement operation management method |
| JP3812884B2 (en) * | 2001-06-27 | 2006-08-23 | 日本電信電話株式会社 | Temperature detection method and temperature detector |
| US6985870B2 (en) | 2002-01-11 | 2006-01-10 | Baxter International Inc. | Medication delivery system |
| US6971790B2 (en) * | 2002-10-11 | 2005-12-06 | Welch Allyn, Inc. | Thermometry probe calibration method |
| US7016715B2 (en) * | 2003-01-13 | 2006-03-21 | Nellcorpuritan Bennett Incorporated | Selection of preset filter parameters based on signal quality |
| WO2004074794A1 (en) * | 2003-02-20 | 2004-09-02 | Ysi Incorporated | Digitally modified resistive output for a temperature sensor |
| US8237386B2 (en) | 2003-08-15 | 2012-08-07 | Apple Inc. | Methods and apparatuses for operating a data processing system |
| US7270476B2 (en) * | 2003-08-21 | 2007-09-18 | Omron Healthcare Co., Ltd. | Electronic clinical thermometer with quick temperature estimating device |
| US7021824B2 (en) * | 2003-10-20 | 2006-04-04 | Welch Allyn, Inc. | Switch assembly for thermometry apparatus |
| JP2005209281A (en) * | 2004-01-22 | 2005-08-04 | Hitachi Global Storage Technologies Netherlands Bv | Data storage device, control method for storage device and magnetic disk driving device |
| US7905857B2 (en) | 2005-07-11 | 2011-03-15 | Covidien Ag | Needle assembly including obturator with safety reset |
| US7850650B2 (en) | 2005-07-11 | 2010-12-14 | Covidien Ag | Needle safety shield with reset |
| US7828773B2 (en) | 2005-07-11 | 2010-11-09 | Covidien Ag | Safety reset key and needle assembly |
| JP4702781B2 (en) * | 2005-03-29 | 2011-06-15 | シチズンホールディングス株式会社 | Electronic thermometer |
| AU2011236062B2 (en) * | 2005-04-01 | 2012-07-05 | Carefusion 303, Inc. | Temperature Prediction System and Method |
| US7318004B2 (en) * | 2005-04-01 | 2008-01-08 | Cardinal Health 303, Inc. | Temperature prediction system and method |
| US20060276747A1 (en) | 2005-06-06 | 2006-12-07 | Sherwood Services Ag | Needle assembly with removable depth stop |
| US20060276772A1 (en) * | 2005-06-06 | 2006-12-07 | Sherwood Services Ag | Bayonet release of safety shield for needle tip |
| US7731692B2 (en) | 2005-07-11 | 2010-06-08 | Covidien Ag | Device for shielding a sharp tip of a cannula and method of using the same |
| JP4949649B2 (en) * | 2005-07-12 | 2012-06-13 | テルモ株式会社 | Electronic thermometer and its control method |
| US8374730B2 (en) * | 2005-08-25 | 2013-02-12 | Apple Inc. | Methods and apparatuses for dynamic thermal control |
| US7562234B2 (en) * | 2005-08-25 | 2009-07-14 | Apple Inc. | Methods and apparatuses for dynamic power control |
| US7654735B2 (en) * | 2005-11-03 | 2010-02-02 | Covidien Ag | Electronic thermometer |
| US20070116089A1 (en) * | 2005-11-22 | 2007-05-24 | Sherwood Services Ag | Electronic thermometer with progress indicator |
| IL178673A0 (en) * | 2006-04-03 | 2007-02-11 | Ilan Menashe | Improved digital thermometer |
| US7314310B2 (en) * | 2006-04-13 | 2008-01-01 | The General Electric Company | Predictive temperature probe with proximity sensor |
| GB0608829D0 (en) * | 2006-05-04 | 2006-06-14 | Husheer Shamus L G | In-situ measurement of physical parameters |
| US20070268954A1 (en) * | 2006-05-19 | 2007-11-22 | Sherwood Services Ag | Portable test apparatus for radiation-sensing thermometer |
| US7507019B2 (en) * | 2006-05-19 | 2009-03-24 | Covidien Ag | Thermometer calibration |
| GB0617451D0 (en) | 2006-09-05 | 2006-10-18 | Medical Prediction Ltd | |
| TW200817658A (en) * | 2006-10-04 | 2008-04-16 | Actherm Inc | Electronic clinical thermometer |
| US7549792B2 (en) | 2006-10-06 | 2009-06-23 | Covidien Ag | Electronic thermometer with selectable modes |
| US7874724B2 (en) * | 2007-04-11 | 2011-01-25 | Trane International Inc. | Method for sensing the liquid level in a compressor |
| EP1988375A1 (en) * | 2007-05-03 | 2008-11-05 | Renishaw plc | Method and apparatus for rapid temperature measurement |
| US8357104B2 (en) | 2007-11-01 | 2013-01-22 | Coviden Lp | Active stylet safety shield |
| CN101435727B (en) * | 2007-11-12 | 2011-01-26 | 深圳迈瑞生物医疗电子股份有限公司 | Temperature predicting method and apparatus |
| US7949889B2 (en) * | 2008-01-07 | 2011-05-24 | Apple Inc. | Forced idle of a data processing system |
| US7949888B2 (en) * | 2008-01-07 | 2011-05-24 | Apple Inc. | Forced idle of a data processing system |
| US8315746B2 (en) | 2008-05-30 | 2012-11-20 | Apple Inc. | Thermal management techniques in an electronic device |
| US8306772B2 (en) * | 2008-10-13 | 2012-11-06 | Apple Inc. | Method for estimating temperature at a critical point |
| FR2939508B1 (en) * | 2008-12-09 | 2011-01-07 | Snecma | METHOD AND SYSTEM FOR CORRECTING MEASUREMENT SIGNAL OF TEMPERATURE. |
| FR2939509B1 (en) * | 2008-12-09 | 2011-03-04 | Snecma | METHOD AND SYSTEM FOR ESTIMATING A VEIN TEMPERATURE IN A TURBOKIN. |
| US8136983B2 (en) * | 2009-01-01 | 2012-03-20 | Mahmoud Razzaghi | Sensor and control system |
| WO2010103436A1 (en) | 2009-03-13 | 2010-09-16 | Koninklijke Philips Electronics N.V. | Zero heat flux temperature sensing device |
| US8794829B2 (en) | 2009-12-31 | 2014-08-05 | Welch Allyn, Inc. | Temperature-measurement probe |
| US8301408B2 (en) * | 2010-03-09 | 2012-10-30 | Invensys Systems, Inc. | Temperature prediction transmitter |
| US9055925B2 (en) | 2010-07-27 | 2015-06-16 | Carefusion 303, Inc. | System and method for reducing false alarms associated with vital-signs monitoring |
| US9585620B2 (en) | 2010-07-27 | 2017-03-07 | Carefusion 303, Inc. | Vital-signs patch having a flexible attachment to electrodes |
| US9615792B2 (en) | 2010-07-27 | 2017-04-11 | Carefusion 303, Inc. | System and method for conserving battery power in a patient monitoring system |
| US9420952B2 (en) | 2010-07-27 | 2016-08-23 | Carefusion 303, Inc. | Temperature probe suitable for axillary reading |
| US9017255B2 (en) | 2010-07-27 | 2015-04-28 | Carefusion 303, Inc. | System and method for saving battery power in a patient monitoring system |
| US9357929B2 (en) * | 2010-07-27 | 2016-06-07 | Carefusion 303, Inc. | System and method for monitoring body temperature of a person |
| US8814792B2 (en) | 2010-07-27 | 2014-08-26 | Carefusion 303, Inc. | System and method for storing and forwarding data from a vital-signs monitor |
| US8380455B1 (en) * | 2012-09-04 | 2013-02-19 | Volodymyr Butsenko | Method for prediction of a response of parameter sensor |
| US9061145B2 (en) | 2012-09-17 | 2015-06-23 | Donald A. Rhodes | Technique for determining optimum treatment parameters |
| US9500535B1 (en) * | 2013-06-27 | 2016-11-22 | Amazon Technologies, Inc. | Determining temperature using multiple sensors |
| AU2015223182B2 (en) | 2014-02-25 | 2017-09-14 | Icu Medical, Inc. | Patient monitoring system with gatekeeper signal |
| US9494470B2 (en) * | 2014-03-20 | 2016-11-15 | Varian Semiconductor Equipment Associates, Inc. | Thermistor based measurement system |
| EP3552586A1 (en) | 2014-12-19 | 2019-10-16 | Coloplast A/S | A system for anal and/or stomal irrigation and a method for controlling such a system |
| EP3488829A1 (en) * | 2014-12-19 | 2019-05-29 | Coloplast A/S | A system for anal or stomal irrigation, and a method for predicting a temperature of an irrigation liquid in a reservoir of such system |
| US10061331B2 (en) * | 2015-01-22 | 2018-08-28 | Qualcomm Incorporated | Systems and methods for detecting thermal runaway |
| JP6674553B2 (en) | 2015-10-19 | 2020-04-01 | アイシーユー・メディカル・インコーポレーテッド | Hemodynamic monitoring system with detachable display unit |
| CN106706165B (en) * | 2015-11-16 | 2019-06-25 | 中国移动通信集团公司 | A kind of method and device of temperature measurement |
| CN105286812B (en) * | 2015-12-02 | 2018-04-20 | 广东宝莱特医用科技股份有限公司 | A kind of body temperature measurement method and device |
| RU2018126877A (en) * | 2015-12-21 | 2020-01-23 | Конинклейке Филипс Н.В. | METHOD FOR FORECASTING TEMPERATURE OF STABILIZATION OF THE HEAT FLOW SENSOR |
| CN105547499B (en) * | 2015-12-29 | 2018-05-25 | 广州视源电子科技股份有限公司 | The method and its system of predicted temperature |
| CA3019170A1 (en) * | 2016-04-08 | 2017-10-12 | Dan Warkander | Temperature-based estimation of scrubbing capacity of a gas scrubber |
| EP3789743A1 (en) * | 2019-09-06 | 2021-03-10 | Basler AG | System and method for generating a prediction signal and a thermal imaging camera |
| EP4161047A1 (en) * | 2021-09-29 | 2023-04-05 | Basler AG | System and method for determining a thermal image of a mobile object |
Family Cites Families (42)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US2755999A (en) | 1952-05-17 | 1956-07-24 | Gen Motors Corp | Temperature measuring and control apparatus |
| US3111032A (en) | 1961-12-07 | 1963-11-19 | Gen Electric | Temperature measurement system |
| US3942123A (en) | 1970-06-15 | 1976-03-02 | Ivac Corporation | Electronic measurement system |
| US3878724A (en) | 1972-11-16 | 1975-04-22 | Gerald F Allen | Physically realizable optimum predictor filter technique for transducers, sensors, and other massive devices, based upon inverse transfers function synthesis |
| US3978325A (en) | 1973-09-05 | 1976-08-31 | Control Electronics Co., Inc. | Electronic thermometer |
| US4068526A (en) | 1975-07-29 | 1978-01-17 | Control Electronics Co., Inc. | Electronic thermometer |
| JPS53109654A (en) * | 1977-03-05 | 1978-09-25 | Nippon Dennetsu Kk | Measuring and weighing apparatus |
| US4158965A (en) | 1978-01-05 | 1979-06-26 | Electromedics, Inc. | Electronic thermometer with heat conditioned probe |
| JPH0238895B2 (en) * | 1979-09-21 | 1990-09-03 | Terumo Corp | DENSHITAIONKEI |
| JPS58225323A (en) * | 1982-06-24 | 1983-12-27 | Terumo Corp | Electronic clinical thermometer |
| US4541734A (en) | 1982-06-24 | 1985-09-17 | Terumo Kabushiki Kaisha | Electronic clinical thermometer, and method of measuring body temperature |
| KR900000822B1 (en) | 1983-04-13 | 1990-02-17 | 다데이시덴기 가부시기가이샤 | Electronic thermometer |
| JPS60114731A (en) * | 1983-11-25 | 1985-06-21 | Omron Tateisi Electronics Co | Electronic clinical thermometer |
| JPS60157031A (en) * | 1984-01-27 | 1985-08-17 | Terumo Corp | Electronic clinical thermometer and method for measuring bodily temperature |
| JPS60209125A (en) * | 1984-04-02 | 1985-10-21 | Omron Tateisi Electronics Co | Electronic clinical thermometer |
| US4727500A (en) | 1985-05-01 | 1988-02-23 | Sherwood Medical Company | Electronic thermometer with fixed response time |
| JPS62165132A (en) | 1986-01-16 | 1987-07-21 | Omron Tateisi Electronics Co | Electronic clinical thermometer |
| JPS62175627A (en) * | 1986-01-29 | 1987-08-01 | Omron Tateisi Electronics Co | Electronic thermometer |
| US4878184A (en) | 1986-02-10 | 1989-10-31 | Omron Tateisi Electronics Co. | Electronic thermometer with predicting means |
| JPH0625700B2 (en) | 1986-03-04 | 1994-04-06 | テルモ株式会社 | Electronic thermometer |
| JPH0792405B2 (en) | 1986-05-13 | 1995-10-09 | オムロン株式会社 | Electronic thermometer |
| JPS639829A (en) * | 1986-07-01 | 1988-01-16 | Takeshi Oshima | Electronic clinical thermometer capable of short-time measurement |
| JPS6340825A (en) | 1986-08-07 | 1988-02-22 | Terumo Corp | Electronic thermometer |
| JPH0656332B2 (en) * | 1986-11-19 | 1994-07-27 | テルモ株式会社 | Electronic thermometer |
| EP0332700A4 (en) * | 1986-11-19 | 1991-08-14 | Terumo Kabushiki Kaisha | Electronic thermometer |
| JPH0795004B2 (en) * | 1986-12-24 | 1995-10-11 | テルモ株式会社 | Body temperature measuring device |
| JPS63262532A (en) * | 1987-04-21 | 1988-10-28 | Citizen Watch Co Ltd | Predicative operation type electronic clinical thermometer |
| US4836442A (en) | 1988-12-14 | 1989-06-06 | Honeywell Inc. | Compensation circuit for sensor lag to reduce undershoot and overshoot |
| JPH07111383B2 (en) * | 1989-10-05 | 1995-11-29 | テルモ株式会社 | Equilibrium temperature detection method and electronic thermometer |
| JPH04109062A (en) * | 1990-08-28 | 1992-04-10 | Hitachi Ltd | Abnormal combustion detecting device and torque control device for internal combustion engine |
| US5392031A (en) | 1992-03-17 | 1995-02-21 | Terumo Kabushiki Kaisha | Electronic clinical thermometer |
| US5519303A (en) * | 1993-09-30 | 1996-05-21 | Motorola, Inc. | Fast battery charging method and apparatus with temperature gradient detection |
| US5428376A (en) * | 1993-10-29 | 1995-06-27 | Hewlett-Packard Company | Thermal turn on energy test for an inkjet printer |
| US5509422A (en) | 1994-04-14 | 1996-04-23 | Fukami; Tetsuji | Clinical thermometer with pulsimeter |
| FR2722313B1 (en) * | 1994-07-07 | 1997-04-25 | Ela Medical Sa | METHOD FOR COMPRESSING PHYSIOLOGICAL DATA, PARTICULARLY CARDIAC ACTIVATED, PARTICULARLY FOR HOLTER RECORDING OF ELECTROCARDIOGRAMS OR ELECTROGRAMS |
| US5632555A (en) | 1994-09-09 | 1997-05-27 | Diatek, L.P. | Medical thermometer |
| JPH08114512A (en) * | 1994-10-17 | 1996-05-07 | Konbi Kk | Prediction type electronic clinical thermometer |
| US5738441A (en) * | 1995-07-11 | 1998-04-14 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Electronic clinical predictive thermometer using logarithm for temperature prediction |
| US5853364A (en) * | 1995-08-07 | 1998-12-29 | Nellcor Puritan Bennett, Inc. | Method and apparatus for estimating physiological parameters using model-based adaptive filtering |
| US5743267A (en) * | 1995-10-19 | 1998-04-28 | Telecom Medical, Inc. | System and method to monitor the heart of a patient |
| US5758652A (en) * | 1995-10-19 | 1998-06-02 | Nikolic; Serjan D. | System and method to measure the condition of a patients heart |
| DE69724591T2 (en) * | 1997-06-02 | 2004-07-01 | European Atomic Energy Community (Euratom) | Online early warning and early detection system for incipient runaway reactions using chaos theory techniques |
-
1999
- 1999-05-18 US US09/313,911 patent/US6270252B1/en not_active Expired - Lifetime
-
2000
- 2000-05-17 EP EP00937576A patent/EP1183508B1/en not_active Expired - Lifetime
- 2000-05-17 WO PCT/US2000/013496 patent/WO2000070316A1/en active IP Right Grant
- 2000-05-17 JP JP2000618700A patent/JP2002544512A/en active Pending
- 2000-05-17 AU AU52725/00A patent/AU781183B2/en not_active Ceased
- 2000-05-17 DE DE60031847T patent/DE60031847T2/en not_active Expired - Lifetime
- 2000-05-17 AT AT00937576T patent/ATE345491T1/en not_active IP Right Cessation
- 2000-05-17 CA CA002373366A patent/CA2373366C/en not_active Expired - Fee Related
- 2000-05-17 ES ES00937576T patent/ES2275515T3/en not_active Expired - Lifetime
-
2001
- 2001-03-14 US US09/809,591 patent/US6698921B2/en not_active Expired - Lifetime
-
2002
- 2002-08-14 HK HK02105961.2A patent/HK1044371B/en not_active IP Right Cessation
-
2011
- 2011-01-17 JP JP2011006924A patent/JP2011075580A/en not_active Ceased
- 2011-06-09 JP JP2011128921A patent/JP2011203269A/en active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| US6270252B1 (en) | 2001-08-07 |
| DE60031847D1 (en) | 2006-12-28 |
| EP1183508B1 (en) | 2006-11-15 |
| US20020003832A1 (en) | 2002-01-10 |
| ES2275515T3 (en) | 2007-06-16 |
| HK1044371A1 (en) | 2002-10-18 |
| EP1183508A1 (en) | 2002-03-06 |
| AU781183B2 (en) | 2005-05-12 |
| JP2011075580A (en) | 2011-04-14 |
| CA2373366A1 (en) | 2000-11-23 |
| DE60031847T2 (en) | 2007-03-01 |
| US6698921B2 (en) | 2004-03-02 |
| JP2002544512A (en) | 2002-12-24 |
| JP2011203269A (en) | 2011-10-13 |
| ATE345491T1 (en) | 2006-12-15 |
| HK1044371B (en) | 2007-08-17 |
| WO2000070316A1 (en) | 2000-11-23 |
| AU5272500A (en) | 2000-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2373366C (en) | Predictive temperature measurement system | |
| US5738441A (en) | Electronic clinical predictive thermometer using logarithm for temperature prediction | |
| EP1864098B1 (en) | Temperature prediction system and method | |
| US5632555A (en) | Medical thermometer | |
| GB2133563A (en) | Electronic clinical thermometer and method of measuring temperature | |
| JPS62165132A (en) | Electronic clinical thermometer | |
| AU2011236062B2 (en) | Temperature Prediction System and Method | |
| SE453432B (en) | ELECTRONIC THERMOMETER | |
| JPS62192625A (en) | Electronic clinical thermometer | |
| JPS62280629A (en) | Electronic clinical thermometer | |
| JPS62175627A (en) | Electronic thermometer | |
| JPS62195530A (en) | Electronic clinical thermometer | |
| JPS62190426A (en) | Electronic clinical thermometer | |
| JPS62185133A (en) | Electronic thermometer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20170517 |