WO2022013911A1 - Temperature measurement device, method, and program - Google Patents

Temperature measurement device, method, and program Download PDF

Info

Publication number
WO2022013911A1
WO2022013911A1 PCT/JP2020/027224 JP2020027224W WO2022013911A1 WO 2022013911 A1 WO2022013911 A1 WO 2022013911A1 JP 2020027224 W JP2020027224 W JP 2020027224W WO 2022013911 A1 WO2022013911 A1 WO 2022013911A1
Authority
WO
WIPO (PCT)
Prior art keywords
temperature
time
equation
estimated temperature
estimated
Prior art date
Application number
PCT/JP2020/027224
Other languages
French (fr)
Japanese (ja)
Inventor
大地 松永
雄次郎 田中
Original Assignee
日本電信電話株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to JP2022535991A priority Critical patent/JP7464124B2/en
Priority to PCT/JP2020/027224 priority patent/WO2022013911A1/en
Publication of WO2022013911A1 publication Critical patent/WO2022013911A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements

Definitions

  • the present invention relates to a temperature measuring device, method and program for non-invasively and accurately measuring the temperature inside a living body.
  • Patent Document 1 Japanese Unexamined Patent Publication No. 2020-003291
  • sensor thermometer elements
  • FIG. 7 shows a schematic diagram of a pseudo one-dimensional model in the temperature measurement inside the living body 71.
  • the temperature Tcore inside the object to be measured (living body) 71 is near the surface (the surface in contact with the object) of the sensor 72 when an object (sensor) having thermal resistance Rs is placed on the surface of the object to be measured (living body) 71. It can be estimated from the temperature Tskin of the above and the temperature Tt near the back surface (the surface on the side in contact with the outside air) of the sensor 72 using the following equation.
  • Tcore Tskin + Rbody x Hskin
  • Rbody is the thermal resistance of the living body
  • Rs is the thermal resistance of the sensor.
  • FIG. 8 shows a comparison between the true internal temperature 81 and the estimated temperature 82 when the wind is blown by a fan.
  • the difference (error) between the true internal temperature 81 and the estimated temperature 82 is that when the wind exceeds a certain threshold, the heat transport mode changes from heat conduction to convection heat transfer, and the amount of heat transferred to the outside is Due to major changes.
  • the temperature measuring method is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is a temperature measuring method inside the living body.
  • the estimated temperature of the temperature is the first equation based on the heat flux calculated by the difference between the temperature of one place near the living body and the temperature of the other place in the sensor, and the estimated temperature by the first formula. Calculated by either of the second equations to be corrected, the first estimated temperature at the first time before the time when the time differential of the heat flux exceeds the reference value and the time differential of the heat flux are the reference.
  • the difference between the second estimated temperature according to the first equation at the second time after the time exceeding the value, the first estimated temperature, and the second according to the second equation at the second time is characterized in that the difference from the estimated temperature of 3 is calculated, and the calculation formula of the estimated temperature is determined to be either the first formula or the second formula based on the two differences. ..
  • the temperature measuring method is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is the first method in the vicinity of the surface in contact with the living body in the sensor.
  • the heat flux is calculated based on the step of measuring the temperature of the first temperature and measuring the second temperature at a position different from the position where the first temperature is measured, and the difference between the first temperature and the second temperature. Steps to detect the time when the time differential of the heat flux exceeds the reference value, the first equation for calculating the estimated temperature based on the heat flux, and the first equation for correcting the estimated temperature according to the first equation.
  • the temperature measuring device is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is the first method in the vicinity of the surface in contact with the living body in the sensor.
  • the heat flux is calculated based on the step of measuring the temperature of the first temperature and measuring the second temperature at a position different from the position where the first temperature is measured, and the difference between the first temperature and the second temperature.
  • the temperature measurement program measures the temperature inside the living body based on the temperature detected by the sensor, and obtains the estimated temperature of the temperature inside the living body in the sensor. Either the first equation based on the heat flux calculated by the difference between the temperature at one location near the living body and the temperature at another location, or the second equation that corrects the estimated temperature according to the first equation.
  • the difference between the second estimated temperature according to the first equation, the first estimated temperature, and the third estimated temperature according to the second equation at the second time is calculated.
  • the temperature measuring device is made to function by executing a process of determining the estimated temperature calculation formula as either the first formula or the second formula based on the difference between the two. ..
  • the present invention it is possible to provide a temperature measuring device, a method and a program for non-invasively and accurately measuring the temperature inside a living body.
  • FIG. 1 is a schematic diagram of a change with time of temperature for explaining the temperature measuring method according to the first embodiment of the present invention.
  • FIG. 2 is a diagram showing a change over time in the time derivative of the heat flux for explaining the temperature measuring method according to the first embodiment of the present invention.
  • FIG. 3 is a block diagram showing a configuration of a temperature measuring device according to the first embodiment of the present invention.
  • FIG. 4 is a flowchart of the temperature measuring method according to the first embodiment of the present invention.
  • FIG. 5 is a diagram showing an example of the temperature measuring method according to the first embodiment of the present invention.
  • FIG. 6 is a diagram showing a configuration example of a computer according to an embodiment of the present invention.
  • FIG. 7 is a schematic diagram of a pseudo one-dimensional model in a conventional temperature measurement inside a living body.
  • FIG. 8 is a diagram showing changes over time in the deep temperature measured by the conventional method.
  • FIG. 1 is a schematic diagram of a change over time in temperature for explaining the method according to the present embodiment.
  • FIG. 1 is based on FIG. 8 described above, and shows a temperature change when a wind is applied.
  • the solid line 11 shows the true internal temperature
  • the dotted line 12 shows the estimated temperature Tcore calculated by the equation (1).
  • Tcore Tskin + Rbody x Hskin (1)
  • Hskin is a heat flux
  • Hskin (Tskin-Tt) / Rs, and is calculated from the difference between Tskin and Tt.
  • Tskin is the temperature near the front surface (for example, the surface in contact with the living body) of the sensor
  • Tt is the temperature near the back surface (the surface in contact with the outside air) of the sensor
  • Rbody is the thermal resistance of the living body
  • Rs is the heat of the sensor. It is resistance.
  • the estimated temperature 12 has an error as compared with the internal temperature 11, it is necessary to correct the estimated temperature 12 including this error in order to calculate an accurate estimated temperature.
  • the change in the estimated temperature 12 is gradual because it depends on the heat capacity of the living body and the sensor. Therefore, it is difficult to correct the estimated temperature 12 by detecting the heat transport mode by the change of the estimated temperature 12. Therefore, in the present embodiment, the change in the heat transport mode is detected by using the time derivative dHskin of the heat flux Hskin, and the error at the estimated temperature 12 is corrected.
  • FIG. 2 shows the time derivative dHskin of the heat flux Hskin with respect to the time course of the temperature shown in FIG. Unlike the change in the estimated temperature, the time derivative dHskin changes sharply when the heat transport mode changes and expresses steep peaks 21 and 22, so that it becomes easy to specify the time.
  • the estimated temperature Tcore can be corrected and calculated by the following formula.
  • Tcore Tskin + ⁇ x Rbody x Hskin (2)
  • the estimated temperature is calculated by both the equations (1) and (2) and compared with the estimated temperature before the detection of the convection change, which is more accurate. Determine the formula for calculating the estimated temperature. The details will be described below.
  • FIG. 3 shows a block diagram of the configuration of the temperature measuring device 30 according to the present embodiment.
  • the temperature measuring device 30 includes a measuring unit (sensor) 31, a storage unit 32, a calculation unit 33, and an output unit 34.
  • the measuring unit (sensor) 31 is arranged near the surface of the object to be measured (for example, a living body) in contact with the surface, and measures the temperature Tskin (first temperature) near the surface of the measuring unit (sensor) 31.
  • the thermometer element 1 and a second thermometer element for measuring the temperature Tt (second temperature) in the vicinity of the back surface of the measuring unit 31 (for example, the surface on the side in contact with the outside air) are provided.
  • the place where the temperature Tt is measured is not limited to the back surface of the measuring unit 31, and may be a place different from the place where the temperature Tskin is measured, and is separated from the place where the temperature Tskin is measured, that is, the first thermometer element. It is desirable that the location is.
  • the second thermometer element may be arranged at a position far enough away from the first temperature and the second temperature so that the heat flux can be accurately calculated.
  • thermometer elements are used, but the number is not limited to two and may be plural. It is possible to measure the heat flux more accurately by using a plurality of thermometer elements.
  • the temperature storage unit 321 stores the measurement time and the measurement temperature measured by the measurement unit 31, and the time constant, the reference value (threshold), and the correction coefficient required for the temperature measurement according to the present embodiment. , Etc. are memorized. Further, the estimated temperature storage unit 322 stores the estimated temperature calculated by the equation (1) and the equation (2) by the calculation unit 331 of the calculation unit 33.
  • the calculation unit 331 calculates the estimated temperature by the equations (1) and (2). In addition, the difference in estimated temperature before and after the detection time of the change in heat flux is calculated. Further, the calculation unit 331 calculates and calculates numerical values necessary for temperature measurement according to the present embodiment, such as heat flux and its time derivative.
  • the comparison unit 332 compares the time derivative dHskin of the heat flux calculated by the calculation unit 331 with the reference value. Further, regarding the difference in the estimated temperature calculated by the calculation unit 331, the difference based on the equation (1) and the difference based on the equation (2) are compared. In addition, the calculation formula of the estimated temperature is determined based on the result of the comparison.
  • the output unit 34 outputs (displays) the estimated temperature calculated by the equations (1) and (2).
  • the accurate temperature cannot be output, it is possible to output (display) a blank (state in which nothing is displayed) or the inability to measure the accurate temperature.
  • FIG. 4 shows a flowchart of the temperature measuring method according to the present embodiment.
  • Tt is measured (step 401).
  • the estimated temperature is calculated by the formula (1) (step 402).
  • the time derivative dHskin of the heat flux Hskin is calculated (step 403).
  • the dHskin is calculated as the difference between the heat flux Hskins at the adjacent sampling (measurement) times.
  • the interval between adjacent sampling (measurement) times is about 1 second.
  • Hskin (ta) ⁇ Tskin (ta) -Tt (ta) ⁇ / Rs.
  • Hskin (tb) ⁇ Tskin (tb) -Tt (tb) ⁇ / Rs. It is calculated by.
  • dHskin Hskin (tb) -Hskin (ta).
  • ta and tb do not have to be adjacent sampling (measurement) times, and may have a predetermined interval.
  • the dHskin value is compared with the reference value (hereinafter referred to as “convection detection threshold”) dHskin_thres (step 404).
  • dHskin_thres When dHskin is equal to or less than the convection detection threshold dHskin_thres (
  • dHskin_thres (
  • dHskin_thres the convection detection threshold dHskin_thres is set to, for example, 0.02 ° C./sec.
  • first time a time t0- ⁇ 1 (hereinafter referred to as “first time”) that is advanced by a predetermined time ⁇ 1 from the time t0 at which the change in heat flux is detected.
  • first estimated temperature a time t0- ⁇ 1 (hereinafter referred to as “first time”) that is advanced by a predetermined time ⁇ 1 from the time t0 at which the change in heat flux is detected.
  • first estimated temperature a temperature Tcore1 (hereinafter referred to as “first estimated temperature”) (step 405).
  • ⁇ 1 can be set from 60 seconds to 300 seconds.
  • the estimated temperature is calculated based on the equation (1).
  • the estimated temperature is calculated by correcting it with equation (2).
  • the estimated temperature Tcore2 (hereinafter, referred to as “second time”) calculated by the equation (1) at the time t0 + ⁇ 2 (hereinafter referred to as “second time”) after a predetermined time ⁇ 2 from the time t0 when the change in heat flux is detected.
  • the “second estimated temperature”) and the estimated temperature Tcore3 calculated by the equation (2) (hereinafter referred to as "third estimated temperature”) are acquired (step 406).
  • ⁇ 2 can be set to about 600 seconds. It is desirable that ⁇ 2 is 300 seconds to 1200 seconds.
  • the difference between the first estimated temperature and the second estimated temperature (
  • the difference between the first estimated temperature and the second estimated temperature is compared with the difference between the first estimated temperature and the third estimated temperature (step 408).
  • the formula for calculating the smaller difference between the second estimated temperature and the third estimated temperature that is, the formula (1) or the formula (2) is estimated. It is used as a temperature calculation formula.
  • FIG. 5 shows the time-dependent change 51 of the time derivative dHskin of the heat flux and the time-dependent change 52 of the estimated temperature in the temperature measuring method of the embodiment according to the present embodiment.
  • the time-dependent change 52 of the estimated temperature indicates an estimated temperature 521 calculated by the formula (1) and an estimated temperature 522 calculated by the formula (2) in order to explain the temperature measuring method of the present embodiment.
  • the time-dependent change 53 of the estimated temperature indicates the estimated temperature 531 by the temperature measuring method of the embodiment according to the present embodiment and the estimated temperature 532 by the temperature measuring method not considering the change of the external temperature as a comparative example. Further, the temperature measured by the eardrum is shown as the true internal temperature 530.
  • the estimated temperature 52 is calculated by the equation (1) because no convection has occurred.
  • the estimated temperature Tcore1 (first estimated temperature) at time t1- ⁇ 1 is selected in the estimated temperature stored in the estimated temperature storage unit 322. (Acquired) (white circle in the figure).
  • the estimated temperature 521 according to the equation (1) and the estimated temperature 522 according to the equation (2) are calculated, and at time t1 + ⁇ 2, the second estimated temperature Tcore2 (black triangle in the figure) which is the estimated temperature 521 according to the equation (1). And the third estimated temperature Tcore3 (black circle in the figure), which is the estimated temperature 522 by the equation (2), is calculated.
  • the first estimated temperature Tcore1 is obtained. Since the difference from the third estimated temperature Tcore3 is smaller, the formula (2) for calculating the third estimated temperature Tcore3 is adopted as the formula for calculating the estimated temperature to measure the temperature and calculate the estimated temperature. continue.
  • the estimated temperature (first estimated temperature) at time t2- ⁇ 1 is selected in the estimated temperature stored in the estimated temperature storage unit 322 (first estimated temperature). Acquired) (white circle in the figure).
  • the estimated temperature 521 according to the equation (1) and the estimated temperature 522 according to the equation (2) are calculated, and at time t1 + ⁇ 2, the second estimated temperature Tcore2 (black triangle in the figure) which is the estimated temperature 521 according to the equation (1). And the third estimated temperature Tcore3 (black circle in the figure), which is the estimated temperature 522 by the equation (2), is calculated.
  • the first estimated temperature Tcore1 is obtained. Since the difference from the third estimated temperature Tcore3 is smaller, the formula (2) for calculating the third estimated temperature Tcore3 is adopted as the formula for calculating the estimated temperature to measure the temperature and calculate the estimated temperature. continue.
  • the estimated temperature 531 calculated by the temperature measuring method of the embodiment according to the present embodiment is calculated by the equation (2) including the correction by the convection change after the time t1, and is the same temperature as the internal temperature 530. Is shown. Further, at time t2, the change in dHskin51 is due to the change in the external temperature, and the convection state does not change, that is, the convection does not end. Therefore, the estimated temperature is calculated from the equation (2) as before t2. The temperature is equivalent to the internal temperature of 530.
  • the estimated temperature 532 of the comparative example after the time t1, the estimated temperature is calculated by the equation (2) including the correction by the convection change, and shows the same temperature as the internal temperature 530.
  • the change in dHskin51 is due to the change in the external temperature, it is erroneously recognized as being due to the change in the convection state.
  • the estimated temperature is calculated from the equation (1), so that a difference (error) from the internal temperature 530 occurs.
  • the temperature measuring method it is possible to measure the temperature inside the living body with high accuracy without erroneously recognizing the convection change when the temperature change occurs.
  • the generation of convection and the end of convection can be detected by measuring the change of the time derivative dHskin of the heat flux in positive and negative directions. For example, when the change in dHskin is 0 or more, the change in the morphology of heat conduction accompanying the start of convection is detected, and when the change in dHskin is less than 0, the change in the morphology of heat conduction accompanying the end of convection is detected. , The convection period can be detected.
  • the calculation unit 33 sequentially uses the temperature measured by the measurement unit (sensor) 31.
  • the calculated estimated temperature is output to the output unit 34.
  • the temperature data stored in the storage unit 32 after being measured by the measurement unit (sensor) 31 is read out (read) by the calculation unit 33.
  • the calculation formula is determined, the estimated temperature is calculated, and it is output to the output unit 34. As a result, it takes, for example, about 20 minutes from the measurement to the output of the estimated temperature.
  • the temperature measured by the measuring unit may be collectively stored in the storage unit, and then the temperature data may be read out (read) to calculate the estimated temperature.
  • the temperature measuring device according to the present embodiment may be attached to the user's body as an integral part of the wearable device.
  • the measuring unit (sensor) 31 is attached to the user's body as a wearable device, and the storage unit 32 and the calculation unit 33 are stored in a smartphone or server outside the wearable device. May be provided.
  • the temperature measuring device is provided with a transmission / reception unit for each of the wearable device and an external server, etc., and the measured temperature measured by the wearable device is transmitted to the server or the like, and is stored and calculated by the server or the like.
  • the estimated temperature and the like may be output to the server or the like, or may be transmitted to the wearable device or the like and output.
  • FIG. 6 shows a configuration example of the computer 60 in the temperature measuring device according to the embodiment of the present invention.
  • the temperature measuring device can be realized by a computer 60 including a CPU (Central Processing Unit) 63, a storage device (storage unit) 62, and an interface device 61, and a program for controlling these hardware resources.
  • the measurement unit and the output unit are connected to the interface device.
  • the CPU executes the process according to the embodiment of the present invention according to the temperature measurement program stored in the storage device. In this way, the temperature measuring program makes the temperature measuring device work.
  • a computer may be provided inside the device, or at least one part of the functions of the computer may be realized by using an external computer.
  • the storage unit may also use the storage medium 64 outside the apparatus, or may read out and execute the temperature measurement program stored in the storage medium 64.
  • the storage medium 64 includes various magnetic recording media, optical magnetic recording media, CD-ROMs, CD-Rs, and various memories. Further, the temperature measurement program may be supplied to the computer via a communication line such as the Internet.
  • the temperature inside the living body can be measured non-invasively and accurately.
  • the present invention can be applied to a deep thermometer used for body temperature control of workers, athletes and the like.
  • Temperature measuring device 31 Measuring unit (sensor) 32 Storage unit 321 Temperature storage unit 322 Estimated temperature storage unit 33 Calculation unit 331 Calculation unit 332 Comparison unit 34 Output unit

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

With a temperature measurement method according to the present invention, the temperature inside a living body is measured on the basis of the temperature detected by means of a sensor (31), wherein: an estimated temperature for the temperature inside the living body is calculated, in the sensor (31), by means of a first formula, which is based on a heat flux calculated on the basis of the difference between the temperature at a site in the vicinity of the living body and the temperature at another site or by means of a second formula, which corrects the estimated temperature calculated by means of the first formula; the difference between a first estimated temperature at a first time, which is before the time at which the time derivative of the heat flux exceeds a reference value, and a second estimated temperature calculated by means of the first formula at a second time, which is after the time at which the time derivative of the heat flux exceeds the reference value, and the difference between the first estimated temperature and a third estimated temperature calculated by means of the second formula at the second time are calculated; and the formula for calculating the estimated temperature is determined to be the first formula or the second formula on the basis of the two differences. As a result, the temperature measurement method according to the present invention is capable of precisely measuring the temperature inside a living body in a non-invasive manner.

Description

温度測定装置、方法およびプログラムTemperature measuring devices, methods and programs
 本発明は、生体内部の温度を非侵襲に精度良く測定する温度測定装置、方法およびプログラムに関する。 The present invention relates to a temperature measuring device, method and program for non-invasively and accurately measuring the temperature inside a living body.
 従来、生体の内部(深部を含む)体温を非侵襲に測定する技術が知られている。例えば、特許文献1(特開2020-003291号公報)は、生体と、複数の温度計素子からなる熱流束センサを備える温度測定センサ(以下、「センサ」という。)と、外気における疑似的な一次元モデルを仮定して、生体の深部体温を推定する技術を開示している。 Conventionally, a technique for non-invasively measuring the body temperature inside a living body (including the deep part) is known. For example, Patent Document 1 (Japanese Unexamined Patent Publication No. 2020-003291) describes a living body, a temperature measuring sensor including a heat flux sensor including a plurality of thermometer elements (hereinafter referred to as “sensor”), and a pseudo in the outside air. It discloses a technique for estimating the core body temperature of a living body assuming a one-dimensional model.
 図7に、生体71内部の温度測定における疑似的な一次元モデルの模式図を示す。被測定物(生体)71の内部の温度Tcoreは、被測定物(生体)71の表面に熱抵抗Rsを有する物体(センサ)を置いたとき、センサ72における表面(その物体と接する面)近傍の温度Tskin、センサ72における裏面(外気と接する側の面)近傍の温度Ttから以下の式を用いて推定できる。 FIG. 7 shows a schematic diagram of a pseudo one-dimensional model in the temperature measurement inside the living body 71. The temperature Tcore inside the object to be measured (living body) 71 is near the surface (the surface in contact with the object) of the sensor 72 when an object (sensor) having thermal resistance Rs is placed on the surface of the object to be measured (living body) 71. It can be estimated from the temperature Tskin of the above and the temperature Tt near the back surface (the surface on the side in contact with the outside air) of the sensor 72 using the following equation.
 Tcore = Tskin + Rbody × Hskin Tcore = Tskin + Rbody x Hskin
 ここで、Hskinは熱流束であり、Hskin=(Tskin-Tt)/Rsで表される。また、Rbodyは生体の熱抵抗、Rsはセンサの熱抵抗である。 Here, Hskin is a heat flux and is represented by Hskin = (Tskin-Tt) / Rs. Further, Rbody is the thermal resistance of the living body, and Rs is the thermal resistance of the sensor.
特開2020-003291号公報Japanese Unexamined Patent Publication No. 2020-003291
 しかしながら、この推定方法では、外気への熱の輸送形態を一定と仮定するため、扇風機などで風を当てる場合や、被測定物がランニングなどにより動く場合には推定温度に誤差が生じる。 However, in this estimation method, since the mode of heat transport to the outside air is assumed to be constant, an error occurs in the estimated temperature when the object is blown by a fan or when the object to be measured moves due to running or the like.
 図8に、扇風機で風を当てた時の真の内部温度81と推定温度82との比較を示す。この真の内部温度81と推定温度82の差(誤差)は、風がある閾値を超えて当てられると、熱の輸送形態が熱伝導から対流熱伝達に変化し、熱が外に伝わる量が大きく変化することに起因する。 FIG. 8 shows a comparison between the true internal temperature 81 and the estimated temperature 82 when the wind is blown by a fan. The difference (error) between the true internal temperature 81 and the estimated temperature 82 is that when the wind exceeds a certain threshold, the heat transport mode changes from heat conduction to convection heat transfer, and the amount of heat transferred to the outside is Due to major changes.
 また、対流が生じることにより熱流束Hskinはセンサの外部へも流れるため、対流時のRbodyが変化し、誤差が生じる。 In addition, since the heat flux Hskin also flows to the outside of the sensor due to the occurrence of convection, the Rbody at the time of convection changes, and an error occurs.
 このように、生体内部の温度の測定において、推定温度に誤差が生じることが問題となっている。また、推定温度の誤差において、誤差が長い時間継続すること、誤差の原因となる熱の対流時に生体の熱抵抗Rbodyが変化することが問題となる。 As described above, there is a problem that an error occurs in the estimated temperature in the measurement of the temperature inside the living body. Further, in the error of the estimated temperature, there are problems that the error continues for a long time and that the thermal resistance Rbody of the living body changes at the time of convection of heat which causes the error.
 上述したような課題を解決するために、本発明に係る温度測定方法は、生体の内部の温度を、センサにより検出された温度を基に測定する温度測定方法であって、前記生体の内部の温度の推定温度を、前記センサにおける前記生体近傍の一の箇所の温度と他の箇所の温度との差分により算出される熱流束に基づく第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより算出し、前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と、前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と、前記第2の時刻での前記第2の式による第3の推定温度との差分を算出して、前記推定温度の算出式を、前記2つの差分に基づき、前記第1の式と前記第2の式とのいずれかに決定することを特徴とする。 In order to solve the above-mentioned problems, the temperature measuring method according to the present invention is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is a temperature measuring method inside the living body. The estimated temperature of the temperature is the first equation based on the heat flux calculated by the difference between the temperature of one place near the living body and the temperature of the other place in the sensor, and the estimated temperature by the first formula. Calculated by either of the second equations to be corrected, the first estimated temperature at the first time before the time when the time differential of the heat flux exceeds the reference value and the time differential of the heat flux are the reference. The difference between the second estimated temperature according to the first equation at the second time after the time exceeding the value, the first estimated temperature, and the second according to the second equation at the second time. It is characterized in that the difference from the estimated temperature of 3 is calculated, and the calculation formula of the estimated temperature is determined to be either the first formula or the second formula based on the two differences. ..
 また、本発明に係る温度測定方法は、生体の内部の温度を、センサにより検出された温度を基に測定する温度測定方法であって、前記センサにおいて、前記生体に接触する面近傍の第1の温度を測定し、前記第1の温度が測定される位置と異なる位置で第2の温度を測定するステップと、前記第1の温度と前記第2の温度の差を基に熱流束を算出するステップと、前記熱流束の時間微分が基準値を超える時刻を検出するステップと、前記熱流束に基づき推定温度を算出する第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより前記推定温度を算出するステップと、前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と、前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度と、前記第2の時刻での前記第2の式による第3の推定温度とを比較して、前記推定温度の算出式を決定するステップとを備える。 Further, the temperature measuring method according to the present invention is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is the first method in the vicinity of the surface in contact with the living body in the sensor. The heat flux is calculated based on the step of measuring the temperature of the first temperature and measuring the second temperature at a position different from the position where the first temperature is measured, and the difference between the first temperature and the second temperature. Steps to detect the time when the time differential of the heat flux exceeds the reference value, the first equation for calculating the estimated temperature based on the heat flux, and the first equation for correcting the estimated temperature according to the first equation. The step of calculating the estimated temperature by any of the equations of 2, the first estimated temperature at the first time before the time when the time differential of the heat flux exceeds the reference value, and the time of the heat flux. Compare the second estimated temperature by the first equation at the second time after the time when the differentiation exceeds the reference value and the third estimated temperature by the second equation at the second time. Further, the step of determining the calculation formula of the estimated temperature is provided.
 また、本発明に係る温度測定装置は、生体の内部の温度を、センサにより検出された温度を基に測定する温度測定方法であって、前記センサにおいて、前記生体に接触する面近傍の第1の温度を測定し、前記第1の温度が測定される位置と異なる位置で第2の温度を測定するステップと、前記第1の温度と前記第2の温度の差を基に熱流束を算出するステップと、前記熱流束の時間微分が基準値を超える時刻を検出するステップと、前記熱流束に基づき推定温度を算出する第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより前記推定温度を算出するステップと、前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と前記第2の時刻での前記第2の式による第3の推定温度との差分とを比較して、前記推定温度の算出式を決定するステップとを備える Further, the temperature measuring device according to the present invention is a temperature measuring method for measuring the temperature inside the living body based on the temperature detected by the sensor, and is the first method in the vicinity of the surface in contact with the living body in the sensor. The heat flux is calculated based on the step of measuring the temperature of the first temperature and measuring the second temperature at a position different from the position where the first temperature is measured, and the difference between the first temperature and the second temperature. A step of detecting a time when the time differentiation of the heat flux exceeds a reference value, a first equation for calculating an estimated temperature based on the heat flux, and a first equation for correcting the estimated temperature according to the first equation. The step of calculating the estimated temperature by any of the equations of 2 and the time differentiation of the first estimated temperature and the heat flux at the first time before the time when the time differentiation of the heat flux exceeds the reference value. According to the difference between the second estimated temperature according to the first equation at the second time after the time when the value exceeds the reference value, and the first estimated temperature and the second equation at the second time. A step of determining the calculation formula of the estimated temperature by comparing with the difference from the third estimated temperature is provided.
 また、本発明に係る温度測定プログラムは、生体の内部の温度を、センサにより検出された温度を基に測定する温度測定装置に対し、前記生体の内部の温度の推定温度を、前記センサにおける前記生体近傍の一の箇所の温度と他の箇所の温度との差分により算出される熱流束に基づく第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより算出し、前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と、前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と、前記第2の時刻での前記第2の式による第3の推定温度との差分を算出して、前記推定温度の算出式を、前記2つの差分に基づき、前記第1の式と前記第2の式とのいずれかに決定する処理を実行させることを特徴とし、温度測定装置を機能させる。 Further, the temperature measurement program according to the present invention measures the temperature inside the living body based on the temperature detected by the sensor, and obtains the estimated temperature of the temperature inside the living body in the sensor. Either the first equation based on the heat flux calculated by the difference between the temperature at one location near the living body and the temperature at another location, or the second equation that corrects the estimated temperature according to the first equation. The first estimated temperature at the first time before the time when the time differential of the heat flux exceeds the reference value, and the second time after the time when the time differential of the heat flux exceeds the reference value. The difference between the second estimated temperature according to the first equation, the first estimated temperature, and the third estimated temperature according to the second equation at the second time is calculated. The temperature measuring device is made to function by executing a process of determining the estimated temperature calculation formula as either the first formula or the second formula based on the difference between the two. ..
 本発明によれば、生体内部の温度を非侵襲に精度良く測定する温度測定装置、方法およびプログラムを提供できる。 According to the present invention, it is possible to provide a temperature measuring device, a method and a program for non-invasively and accurately measuring the temperature inside a living body.
図1は、本発明の第1の実施の形態に係る温度測定方法を説明するための温度の経時変化の概要図である。FIG. 1 is a schematic diagram of a change with time of temperature for explaining the temperature measuring method according to the first embodiment of the present invention. 図2は、本発明の第1の実施の形態に係る温度測定方法を説明するための熱流束の時間微分の経時変化を示す図である。FIG. 2 is a diagram showing a change over time in the time derivative of the heat flux for explaining the temperature measuring method according to the first embodiment of the present invention. 図3は、本発明の第1の実施の形態に係る温度測定装置の構成を示すブロック図である。FIG. 3 is a block diagram showing a configuration of a temperature measuring device according to the first embodiment of the present invention. 図4は、本発明の第1の実施の形態に係る温度測定方法のフローチャート図である。FIG. 4 is a flowchart of the temperature measuring method according to the first embodiment of the present invention. 図5は、本発明の第1の実施の形態に係る温度測定方法の実施例を示す図である。FIG. 5 is a diagram showing an example of the temperature measuring method according to the first embodiment of the present invention. 図6は、本発明の実施の形態におけるコンピュータの構成例を示す図である。FIG. 6 is a diagram showing a configuration example of a computer according to an embodiment of the present invention. 図7は、従来の生体内部の温度測定における疑似的な一次元モデルの模式図である。FIG. 7 is a schematic diagram of a pseudo one-dimensional model in a conventional temperature measurement inside a living body. 図8は、従来法により測定される深部温度の経時変化を示す図である。FIG. 8 is a diagram showing changes over time in the deep temperature measured by the conventional method.
 図1は、本実施の形態に係る方法を説明するための温度の経時変化の概要図である。図1は、前述の図8に基づくものであり、風を当てた時の温度変化を示す。実線11は真の内部温度、点線12は式(1)により算出される推定温度Tcoreを示す。 FIG. 1 is a schematic diagram of a change over time in temperature for explaining the method according to the present embodiment. FIG. 1 is based on FIG. 8 described above, and shows a temperature change when a wind is applied. The solid line 11 shows the true internal temperature, and the dotted line 12 shows the estimated temperature Tcore calculated by the equation (1).
 Tcore = Tskin + Rbody × Hskin (1) Tcore = Tskin + Rbody x Hskin (1)
 ここで、Hskinは熱流束であり、Hskin=(Tskin-Tt)/Rsで表され、TskinとTtの差分より算出される。また、Tskinはセンサにおける表面(例えば、生体と接する面)近傍の温度、Ttはセンサにおける裏面(外気と接する側の面)近傍の温度であり、Rbodyは生体の熱抵抗、Rsはセンサの熱抵抗である。 Here, Hskin is a heat flux, is represented by Hskin = (Tskin-Tt) / Rs, and is calculated from the difference between Tskin and Tt. Further, Tskin is the temperature near the front surface (for example, the surface in contact with the living body) of the sensor, Tt is the temperature near the back surface (the surface in contact with the outside air) of the sensor, Rbody is the thermal resistance of the living body, and Rs is the heat of the sensor. It is resistance.
 風を当てることにより、対流の開始にともなう熱輸送形態の変化が生じると、推定温度12には正のピーク(上に凸)が発現する。また、風を停止して対流の終了にともなう熱輸送形態(対流)の変化が生じると、推定温度12には負のピーク(下に凸)が発現する。 When the heat transport mode changes with the start of convection by applying wind, a positive peak (convex upward) appears at the estimated temperature 12. Further, when the wind is stopped and the heat transport mode (convection) changes with the end of convection, a negative peak (convex downward) appears at the estimated temperature 12.
 このように、推定温度12には、内部温度11に比べて、誤差が生じるので、正確な推定温度を算出するためには、この誤差を含む推定温度12を補正する必要がある。 As described above, since the estimated temperature 12 has an error as compared with the internal temperature 11, it is necessary to correct the estimated temperature 12 including this error in order to calculate an accurate estimated temperature.
 推定温度の補正においては、まず、補正を要するとき、すなわち熱輸送形態(対流)の変化が生じるときを検知する必要がある。図1に示すように、推定温度12の変化は、生体やセンサの熱容量に依存するため、緩やかである。そこで、推定温度12の変化により熱輸送形態を検出して、推定温度12を補正することは難しい。そこで、本実施の形態においては、熱流束Hskinの時間微分dHskinを用いて熱輸送形態の変化を検出して、推定温度12における誤差を補正する。 In the correction of the estimated temperature, it is first necessary to detect when the correction is required, that is, when the heat transport mode (convection) changes. As shown in FIG. 1, the change in the estimated temperature 12 is gradual because it depends on the heat capacity of the living body and the sensor. Therefore, it is difficult to correct the estimated temperature 12 by detecting the heat transport mode by the change of the estimated temperature 12. Therefore, in the present embodiment, the change in the heat transport mode is detected by using the time derivative dHskin of the heat flux Hskin, and the error at the estimated temperature 12 is corrected.
 図2に、図1に示す温度の経時変化における熱流束Hskinの時間微分dHskinを示す。推定温度の変化と異なり、時間微分dHskinは熱輸送形態が変化すると鋭敏に変化し急峻なピーク21、22を発現するので、その時刻を特定することが容易になる。 FIG. 2 shows the time derivative dHskin of the heat flux Hskin with respect to the time course of the temperature shown in FIG. Unlike the change in the estimated temperature, the time derivative dHskin changes sharply when the heat transport mode changes and expresses steep peaks 21 and 22, so that it becomes easy to specify the time.
 また、推定温度Tcoreは、以下の式により、補正して算出することができる。 Further, the estimated temperature Tcore can be corrected and calculated by the following formula.
 Tcore = Tskin + α × Rbody × Hskin(2) Tcore = Tskin + α x Rbody x Hskin (2)
 ここで、αは補正係数であり、1.03以上1.15以下であることが望ましい。本実施の形形態では、α = 1.05とする。このように、式(2)は式(1)による推定温度を補正するものである。 Here, α is a correction coefficient, and it is desirable that it is 1.03 or more and 1.15 or less. In this embodiment, α = 1.05. As described above, the equation (2) corrects the estimated temperature according to the equation (1).
 以上のように、推定温度には対流変化により誤差が生じるため、対流変化の検知を契機に補正を行う。具体的には、熱流束Hskinの時間微分dHskinが基準値(閾値)を超えた場合に対流変化が生じたと判定して補正を行う。 As described above, an error occurs in the estimated temperature due to the convection change, so the correction is performed when the convection change is detected. Specifically, when the time derivative dHskin of the heat flux Hskin exceeds the reference value (threshold value), it is determined that a convection change has occurred and correction is performed.
 しかしながら、外気温や室内温度などのセンサ外部の環境温度(以下、「外部温度」という。)が変化したときも同様に、Hskinの時間微分dHskinが閾値を超える。その結果、対流変化が生じていないにもかかわらず、対流変化が生じたものと誤認識(検知)することになる。この場合には、本来用いるべき式(1)ではなく式(2)を用いて推定温度が算出されるため、誤差が生じる。 However, when the environmental temperature outside the sensor (hereinafter referred to as "external temperature") such as the outside air temperature and the room temperature changes, the time derivative dHskin of Hskin also exceeds the threshold value. As a result, even though the convection change has not occurred, it is erroneously recognized (detected) that the convection change has occurred. In this case, since the estimated temperature is calculated using the formula (2) instead of the formula (1) that should be originally used, an error occurs.
 この誤差を解消するために、本実施の形態では、式(1)と式(2)両方の式で推定温度を算出して、対流変化の検知前の推定温度と比較して、より正確な推定温度を算出する式を決定する。以下に、詳細を説明する。 In order to eliminate this error, in the present embodiment, the estimated temperature is calculated by both the equations (1) and (2) and compared with the estimated temperature before the detection of the convection change, which is more accurate. Determine the formula for calculating the estimated temperature. The details will be described below.
<第1の実施の形態>
 本発明の第1の実施の形態について図1~図6を参照して説明する。
<First Embodiment>
The first embodiment of the present invention will be described with reference to FIGS. 1 to 6.
<温度測定装置の構成>
 図3に、本実施の形態に係る温度測定装置30の構成のブロック図を示す。温度測定装置30は、測定部(センサ)31と、記憶部32と、演算部33と、出力部34とを備える。
<Configuration of temperature measuring device>
FIG. 3 shows a block diagram of the configuration of the temperature measuring device 30 according to the present embodiment. The temperature measuring device 30 includes a measuring unit (sensor) 31, a storage unit 32, a calculation unit 33, and an output unit 34.
 測定部(センサ)31は、被測定物(例えば、生体)の表面と接触する面近傍に配置され、測定部(センサ)31の表面近傍での温度Tskin(第1の温度)を測定する第1の温度計素子と、測定部31の裏面(例えば、外気と接する側の面)近傍での温度Tt(第2の温度)を測定する第2の温度計素子を備える。 The measuring unit (sensor) 31 is arranged near the surface of the object to be measured (for example, a living body) in contact with the surface, and measures the temperature Tskin (first temperature) near the surface of the measuring unit (sensor) 31. The thermometer element 1 and a second thermometer element for measuring the temperature Tt (second temperature) in the vicinity of the back surface of the measuring unit 31 (for example, the surface on the side in contact with the outside air) are provided.
 ここで、温度Ttを測定する箇所は、測定部31の裏面に限らず、温度Tskinを測定する箇所とは異なる箇所であればよく、温度Tskinを測定する箇所すなわち第1の温度計素子から離れた箇所であることが望ましい。第2の温度計素子は、第1の温度と第2の温度から正確に熱流束を算出できる程度離れた箇所に配置されればよい。 Here, the place where the temperature Tt is measured is not limited to the back surface of the measuring unit 31, and may be a place different from the place where the temperature Tskin is measured, and is separated from the place where the temperature Tskin is measured, that is, the first thermometer element. It is desirable that the location is. The second thermometer element may be arranged at a position far enough away from the first temperature and the second temperature so that the heat flux can be accurately calculated.
 また、本実施の形態では、2個の温度計素子を用いるが、2個に限らず複数であればよい。複数の温度計素子を用いた方が精度よく熱流束を測定できる。 Further, in the present embodiment, two thermometer elements are used, but the number is not limited to two and may be plural. It is possible to measure the heat flux more accurately by using a plurality of thermometer elements.
 記憶部32において、温度等記憶部321は、測定部31で測定される測定時刻と測定温度を記憶し、本実施の形態に係る温度測定に必要な時定数、基準値(閾値)、補正係数、などの数値等を記憶する。また、推定温度記憶部322は、演算部33の算出部331で式(1)と式(2)により算出される推定温度を記憶する。 In the storage unit 32, the temperature storage unit 321 stores the measurement time and the measurement temperature measured by the measurement unit 31, and the time constant, the reference value (threshold), and the correction coefficient required for the temperature measurement according to the present embodiment. , Etc. are memorized. Further, the estimated temperature storage unit 322 stores the estimated temperature calculated by the equation (1) and the equation (2) by the calculation unit 331 of the calculation unit 33.
 演算部33において、算出部331は、式(1)と式(2)により推定温度を算出する。また、熱流束の変化の検知時刻の前後での推定温度の差分を算出する。また、算出部331は、熱流束やその時間微分などの本実施の形態に係る温度測定に必要な数値等について算出、演算する。 In the calculation unit 33, the calculation unit 331 calculates the estimated temperature by the equations (1) and (2). In addition, the difference in estimated temperature before and after the detection time of the change in heat flux is calculated. Further, the calculation unit 331 calculates and calculates numerical values necessary for temperature measurement according to the present embodiment, such as heat flux and its time derivative.
 また、比較部332は、算出部331で算出される熱流束の時間微分dHskinを基準値と比較する。また、算出部331で算出される推定温度の差分について、式(1)に基づく差分と式(2)に基づく差分を比較する。また、比較の結果に基づき、推定温度の算出式を決定する。 Further, the comparison unit 332 compares the time derivative dHskin of the heat flux calculated by the calculation unit 331 with the reference value. Further, regarding the difference in the estimated temperature calculated by the calculation unit 331, the difference based on the equation (1) and the difference based on the equation (2) are compared. In addition, the calculation formula of the estimated temperature is determined based on the result of the comparison.
 出力部34は、式(1)と式(2)により算出される推定温度を出力(表示)する。正確な温度を出力できない場合には、ブランク(何も表示しない状態)や正確な温度を測定できないことを出力(表示)することもできる。 The output unit 34 outputs (displays) the estimated temperature calculated by the equations (1) and (2). When the accurate temperature cannot be output, it is possible to output (display) a blank (state in which nothing is displayed) or the inability to measure the accurate temperature.
<温度測定方法>
 本実施の形態に係る方法の概要を、図4を参照して説明する。図4に、本実施の形態に係る温度測定方法のフローチャート図を示す。
<Temperature measurement method>
An outline of the method according to the present embodiment will be described with reference to FIG. FIG. 4 shows a flowchart of the temperature measuring method according to the present embodiment.
 初めに、センサ(測定部)により、センサの表面(生体と接する面)近傍の温度(第1の温度)Tskinと、センサの裏面(外気と接する側の面)近傍の温度(第2の温度)Ttを測定する(ステップ401)。 First, by the sensor (measurement unit), the temperature near the front surface (surface in contact with the living body) of the sensor (first temperature) Tskin and the temperature near the back surface of the sensor (surface in contact with outside air) (second temperature). ) Tt is measured (step 401).
 次に、式(1)により、推定温度を算出する(ステップ402)。 Next, the estimated temperature is calculated by the formula (1) (step 402).
 次に、熱流束Hskinの時間微分dHskinを算出する(ステップ403)。dHskinは、隣接するサンプリング(測定)時刻での熱流束Hskinの差分として算出する。ここで、例えば、隣接するサンプリング(測定)時刻の間隔は1秒程度である。 Next, the time derivative dHskin of the heat flux Hskin is calculated (step 403). The dHskin is calculated as the difference between the heat flux Hskins at the adjacent sampling (measurement) times. Here, for example, the interval between adjacent sampling (measurement) times is about 1 second.
 例えば、時刻taでTskin(ta)とTt(ta)が測定されるとき、Hskin(ta)は、Hskin(ta)={Tskin(ta)-Tt(ta)}/Rsで算出される。 For example, when Tskin (ta) and Tt (ta) are measured at time ta, Hskin (ta) is calculated by Hskin (ta) = {Tskin (ta) -Tt (ta)} / Rs.
 同様に、時刻taに引き続き、時刻tbでTskin(tb)とTt(tb)が測定されるとき、Hskin(tb)は、Hskin(tb)={Tskin(tb)-Tt(tb)}/Rsで算出される。 Similarly, when Tskin (tb) and Tt (tb) are measured at time tb following time ta, Hskin (tb) is Hskin (tb) = {Tskin (tb) -Tt (tb)} / Rs. It is calculated by.
 このとき、dHskinは、dHskin=Hskin(tb)-Hskin(ta)で算出される。 At this time, dHskin is calculated by dHskin = Hskin (tb) -Hskin (ta).
 ここで、taとtbは隣接するサンプリング(測定)時刻でなくてもよく、所定の間隔を有してもよい。 Here, ta and tb do not have to be adjacent sampling (measurement) times, and may have a predetermined interval.
 次に、dHskinの値が基準値(以下、「対流検出閾値」という。)dHskin_thresを比較する(ステップ404)。 Next, the dHskin value is compared with the reference value (hereinafter referred to as “convection detection threshold”) dHskin_thres (step 404).
 dHskinが対流検出閾値dHskin_thres以下の場合には(|dHskin|≦dHskin_thres)、対流によって熱伝達の形態が変化していないと判定して、以下のステップを行わず、上述のTcoreを、この時刻での推定温度と決定して、引き続き、次の時刻で測定を実行する。 When dHskin is equal to or less than the convection detection threshold dHskin_thres (| dHskin | ≦ dHskin_thres), it is determined that the form of heat transfer has not changed due to convection. Determine the estimated temperature of and continue to measure at the next time.
 一方、dHskinが対流検出閾値dHskin_thresを上回る(|dHskin|>dHskin_thres)場合、すなわち熱流束の変化が検出される場合には、このdHskinの変化が対流の変化によるものか、外部温度の変化など他の要因によるものかを、以下の通り、判定する。ここで、対流検出閾値dHskin_thresは、例えば、0.02℃/secとする。 On the other hand, when dHskin exceeds the convection detection threshold dHskin_thres (| dHskin |> dHskin_thres), that is, when a change in heat flux is detected, whether this change in dHskin is due to a change in convection, a change in external temperature, or the like. It is determined as follows whether it is due to the factor of. Here, the convection detection threshold dHskin_thres is set to, for example, 0.02 ° C./sec.
 まず、推定温度記憶部322に記憶される推定温度において、熱流束の変化が検出される時刻t0より所定時間τ1遡った時刻t0-τ1(以下、「第1の時刻」という。)での推定温度Tcore1(以下、「第1の推定温度」という。)を選択(取得)する(ステップ405)。ここで、τ1は、60秒間~300秒間を設定できる。 First, at the estimated temperature stored in the estimated temperature storage unit 322, estimation is performed at a time t0-τ1 (hereinafter referred to as “first time”) that is advanced by a predetermined time τ1 from the time t0 at which the change in heat flux is detected. Select (acquire) the temperature Tcore1 (hereinafter referred to as “first estimated temperature”) (step 405). Here, τ1 can be set from 60 seconds to 300 seconds.
 次に、dHskinの変化が外部温度の変化などの要因によると想定して、式(1)に基づき推定温度を算出する。 Next, assuming that the change in dHskin is due to factors such as the change in external temperature, the estimated temperature is calculated based on the equation (1).
 同時に、dHskinの変化が対流の変化などの要因によると想定して、式(2)により補正して、推定温度を算出する。 At the same time, assuming that the change in dHskin is due to factors such as changes in convection, the estimated temperature is calculated by correcting it with equation (2).
 次に、熱流束の変化が検出される時刻t0より所定時間τ2後の時刻t0+τ2(以下、「第2の時刻」という。)での、式(1)により算出される推定温度Tcore2(以下、「第2の推定温度」という。)と、式(2)により算出される推定温度Tcore3(以下、「第3の推定温度」という。)を取得する(ステップ406)。ここで、τ2は、600秒間程度に設定できる。τ2は、300秒~1200秒であることが望ましい。 Next, the estimated temperature Tcore2 (hereinafter, referred to as “second time”) calculated by the equation (1) at the time t0 + τ2 (hereinafter referred to as “second time”) after a predetermined time τ2 from the time t0 when the change in heat flux is detected. The "second estimated temperature") and the estimated temperature Tcore3 calculated by the equation (2) (hereinafter referred to as "third estimated temperature") are acquired (step 406). Here, τ2 can be set to about 600 seconds. It is desirable that τ2 is 300 seconds to 1200 seconds.
 次に、第1の推定温度と第2の推定温度との差分(|Tcore1-Tcore2|)を算出する。また、第1の推定温度と第3の推定温度(|Tcore1-Tcore3|)との差分を算出する(ステップ407)。 Next, the difference between the first estimated temperature and the second estimated temperature (| Tcore1-Tcore2 |) is calculated. Further, the difference between the first estimated temperature and the third estimated temperature (| Tcore1-Tcore3 |) is calculated (step 407).
 次に、第1の推定温度と第2の推定温度との差分と、第1の推定温度と第3の推定温度との差分とを比較する(ステップ408)。 Next, the difference between the first estimated temperature and the second estimated temperature is compared with the difference between the first estimated temperature and the third estimated temperature (step 408).
 これらの差分を比較した結果、第2の推定温度と第3の推定温度とのいずれかのうち、差分の小さい方を算出する式、すなわち式(1)、式(2)のいずれかを推定温度の算出式として採用する。 As a result of comparing these differences, the formula for calculating the smaller difference between the second estimated temperature and the third estimated temperature, that is, the formula (1) or the formula (2) is estimated. It is used as a temperature calculation formula.
 採用した算出式を用いて、温度の測定、推定温度の算出を継続する(ステップ409、410)。 Continue measuring the temperature and calculating the estimated temperature using the adopted calculation formula (steps 409, 410).
<実施例>
 本実施の形態に係る温度測定方法の一例としての実施例を、図5を参照に説明する。
<Example>
An example of the temperature measuring method according to the present embodiment will be described with reference to FIG.
 図5に、本実施の形態に係る実施例の温度測定方法における、熱流束の時間微分dHskinの経時変化51と推定温度の経時変化52を示す。推定温度の経時変化52は、本実施例の温度測定方法を説明するために、式(1)により算出される推定温度521と式(2)により算出される推定温度522とを示す。 FIG. 5 shows the time-dependent change 51 of the time derivative dHskin of the heat flux and the time-dependent change 52 of the estimated temperature in the temperature measuring method of the embodiment according to the present embodiment. The time-dependent change 52 of the estimated temperature indicates an estimated temperature 521 calculated by the formula (1) and an estimated temperature 522 calculated by the formula (2) in order to explain the temperature measuring method of the present embodiment.
 また、推定温度の経時変化53は、本実施の形態に係る実施例の温度測定方法による推定温度531と、比較例として、外部温度の変化を考慮しない温度測定方法による推定温度532を示す。また、真の内部温度530として、鼓膜で測定される温度を示す。 Further, the time-dependent change 53 of the estimated temperature indicates the estimated temperature 531 by the temperature measuring method of the embodiment according to the present embodiment and the estimated temperature 532 by the temperature measuring method not considering the change of the external temperature as a comparative example. Further, the temperature measured by the eardrum is shown as the true internal temperature 530.
 対流変化が生じていない状態から温度測定を開始する場合を一例として説明する。初めに、推定温度52は、対流が生じていないので、式(1)で算出される。 The case where the temperature measurement is started from the state where the convection change does not occur will be described as an example. First, the estimated temperature 52 is calculated by the equation (1) because no convection has occurred.
 次に、時刻t1でdHskin51が基準値(閾値)511を超えたとき、推定温度記憶部322に記憶される推定温度において、時刻t1-τ1での推定温度Tcore1(第1の推定温度)が選択(取得)される(図中、白丸)。 Next, when dHskin51 exceeds the reference value (threshold value) 511 at time t1, the estimated temperature Tcore1 (first estimated temperature) at time t1-τ1 is selected in the estimated temperature stored in the estimated temperature storage unit 322. (Acquired) (white circle in the figure).
 引き続き、式(1)による推定温度521と式(2)による推定温度522が算出され、時刻t1+τ2で、式(1)による推定温度521である第2の推定温度Tcore2(図中、黒三角)と、式(2)による推定温度522である第3の推定温度Tcore3(図中、黒丸)が算出される。 Subsequently, the estimated temperature 521 according to the equation (1) and the estimated temperature 522 according to the equation (2) are calculated, and at time t1 + τ2, the second estimated temperature Tcore2 (black triangle in the figure) which is the estimated temperature 521 according to the equation (1). And the third estimated temperature Tcore3 (black circle in the figure), which is the estimated temperature 522 by the equation (2), is calculated.
 次に、第1の推定温度Tcore1と第2の推定温度Tcore2との差分と、第1の推定温度Tcore1と第3の推定温度Tcore3との差分とを比較した結果、第1の推定温度Tcore1と第3の推定温度Tcore3との差分の方が小さいので、第3の推定温度Tcore3を算出する式(2)を、推定温度を算出する式に採用して、温度の測定、推定温度の算出を継続する。 Next, as a result of comparing the difference between the first estimated temperature Tcore1 and the second estimated temperature Tcore2 and the difference between the first estimated temperature Tcore1 and the third estimated temperature Tcore3, the first estimated temperature Tcore1 is obtained. Since the difference from the third estimated temperature Tcore3 is smaller, the formula (2) for calculating the third estimated temperature Tcore3 is adopted as the formula for calculating the estimated temperature to measure the temperature and calculate the estimated temperature. continue.
 次に、時刻t2でdHskin51が基準値(閾値)511を超えたとき、推定温度記憶部322に記憶される推定温度において、時刻t2-τ1での推定温度(第1の推定温度)が選択(取得)される(図中、白丸)。 Next, when dHskin51 exceeds the reference value (threshold value) 511 at time t2, the estimated temperature (first estimated temperature) at time t2-τ1 is selected in the estimated temperature stored in the estimated temperature storage unit 322 (first estimated temperature). Acquired) (white circle in the figure).
 引き続き、式(1)による推定温度521と式(2)による推定温度522が算出され、時刻t1+τ2で、式(1)による推定温度521である第2の推定温度Tcore2(図中、黒三角)と、式(2)による推定温度522である第3の推定温度Tcore3(図中、黒丸)が算出される。 Subsequently, the estimated temperature 521 according to the equation (1) and the estimated temperature 522 according to the equation (2) are calculated, and at time t1 + τ2, the second estimated temperature Tcore2 (black triangle in the figure) which is the estimated temperature 521 according to the equation (1). And the third estimated temperature Tcore3 (black circle in the figure), which is the estimated temperature 522 by the equation (2), is calculated.
 次に、第1の推定温度Tcore1と第2の推定温度Tcore2との差分と、第1の推定温度Tcore1と第3の推定温度Tcore3との差分とを比較した結果、第1の推定温度Tcore1と第3の推定温度Tcore3との差分の方が小さいので、第3の推定温度Tcore3を算出する式(2)を、推定温度を算出する式に採用して、温度の測定、推定温度の算出を継続する。 Next, as a result of comparing the difference between the first estimated temperature Tcore1 and the second estimated temperature Tcore2 and the difference between the first estimated temperature Tcore1 and the third estimated temperature Tcore3, the first estimated temperature Tcore1 is obtained. Since the difference from the third estimated temperature Tcore3 is smaller, the formula (2) for calculating the third estimated temperature Tcore3 is adopted as the formula for calculating the estimated temperature to measure the temperature and calculate the estimated temperature. continue.
 本実施の形態に係る実施例の温度測定方法により算出される推定温度531では、時刻t1以降、推定温度は対流変化による補正を含む式(2)により算出されて、内部温度530と同等の温度を示す。また、時刻t2で、dHskin51の変化が外部温度の変化によるものであり、対流状態は変化せず、すなわち、対流は終了しないので、t2以前と同様に、推定温度は式(2)より算出されて、内部温度530と同等の温度を示す。 In the estimated temperature 531 calculated by the temperature measuring method of the embodiment according to the present embodiment, the estimated temperature is calculated by the equation (2) including the correction by the convection change after the time t1, and is the same temperature as the internal temperature 530. Is shown. Further, at time t2, the change in dHskin51 is due to the change in the external temperature, and the convection state does not change, that is, the convection does not end. Therefore, the estimated temperature is calculated from the equation (2) as before t2. The temperature is equivalent to the internal temperature of 530.
 一方、比較例の推定温度532では、時刻t1以降、推定温度は対流変化による補正を含む式(2)により算出されて、内部温度530と同等の温度を示す。しかしながら、時刻t2で、dHskin51の変化が外部温度の変化によるものであるにもかかわらず、対流状態の変化によるものと誤認識される。その結果、対流が終了したものとみなされ、推定温度は式(1)より算出されるので、内部温度530と差異(誤差)が生じる。 On the other hand, in the estimated temperature 532 of the comparative example, after the time t1, the estimated temperature is calculated by the equation (2) including the correction by the convection change, and shows the same temperature as the internal temperature 530. However, at time t2, although the change in dHskin51 is due to the change in the external temperature, it is erroneously recognized as being due to the change in the convection state. As a result, it is considered that the convection has ended, and the estimated temperature is calculated from the equation (1), so that a difference (error) from the internal temperature 530 occurs.
 このように、本実施の形態に係る温度測定方法によれば、気温変化発生時に対流変化の誤認識をすることなく、高精度で生体内部の温度を測定できる。 As described above, according to the temperature measuring method according to the present embodiment, it is possible to measure the temperature inside the living body with high accuracy without erroneously recognizing the convection change when the temperature change occurs.
 本発明に係る実施の形態では、熱流束の時間微分dHskinの変化を正負で測定することにより、対流発生と対流終了を検出できる。例えば、dHskinの変化が0以上の時に、対流開始にともなう熱伝導の形態の変化を検出して、dHskinの変化が0未満の時に、対流終了にともなう熱伝導の形態の変化を検出することにより、対流期間を検出できる。 In the embodiment of the present invention, the generation of convection and the end of convection can be detected by measuring the change of the time derivative dHskin of the heat flux in positive and negative directions. For example, when the change in dHskin is 0 or more, the change in the morphology of heat conduction accompanying the start of convection is detected, and when the change in dHskin is less than 0, the change in the morphology of heat conduction accompanying the end of convection is detected. , The convection period can be detected.
 本発明に係る実施の形態では、熱流束の時間微分の変化が検知されずに対流が生じていない場合には、測定部(センサ)31で測定された温度を基に、順次演算部33で算出される推定温度が、出力部34に出力される。 In the embodiment of the present invention, when the change in the time derivative of the heat flux is not detected and the convection does not occur, the calculation unit 33 sequentially uses the temperature measured by the measurement unit (sensor) 31. The calculated estimated temperature is output to the output unit 34.
 一方、熱流束の時間微分の変化が検知される場合には、測定部(センサ)31で測定された後に、記憶部32で記憶された温度データを、演算部33で読み出して(読み込んで)算出式が決定され、推定温度が算出され、出力部34に出力される。その結果、測定から推定温度の出力まで、例えば、20分間程度要する。 On the other hand, when the change in the time derivative of the heat flux is detected, the temperature data stored in the storage unit 32 after being measured by the measurement unit (sensor) 31 is read out (read) by the calculation unit 33. The calculation formula is determined, the estimated temperature is calculated, and it is output to the output unit 34. As a result, it takes, for example, about 20 minutes from the measurement to the output of the estimated temperature.
 また、測定部(センサ)で測定された温度をまとめて記憶部に記憶してから、温度データを読み出して(読み込んで)推定温度を算出してもよい。 Alternatively, the temperature measured by the measuring unit (sensor) may be collectively stored in the storage unit, and then the temperature data may be read out (read) to calculate the estimated temperature.
 本実施の形態に係る温度測定装置は、ウェアラブルデバイス一体としてユーザの身体に装着されてもよい。 The temperature measuring device according to the present embodiment may be attached to the user's body as an integral part of the wearable device.
 または、本発明に係る実施の形態に係る温度測定装置は、測定部(センサ)31をウェアラブルデバイスとしてユーザの身体に装着して、ウェアラブルデバイス外部のスマートフォンやサーバ等に記憶部32、演算部33を備えてもよい。この場合、温度測定装置はウェアラブルデバイスと外部のサーバ等それぞれに送受信部を備え、ウェアラブルデバイスで測定される測定温度をサーバ等に送信し、サーバ等で記憶、計算を行う。最後に、推定温度等(深部温度は測定されないことの表示等を含む)はサーバ等に出力されてもよいし、ウェアラブルデバイス等に送信されて出力されてもよい。 Alternatively, in the temperature measuring device according to the embodiment of the present invention, the measuring unit (sensor) 31 is attached to the user's body as a wearable device, and the storage unit 32 and the calculation unit 33 are stored in a smartphone or server outside the wearable device. May be provided. In this case, the temperature measuring device is provided with a transmission / reception unit for each of the wearable device and an external server, etc., and the measured temperature measured by the wearable device is transmitted to the server or the like, and is stored and calculated by the server or the like. Finally, the estimated temperature and the like (including the display that the deep temperature is not measured) may be output to the server or the like, or may be transmitted to the wearable device or the like and output.
<コンピュータの構成例>
 図6に、本発明の実施の形態に係る温度測定装置におけるコンピュータ60の構成例を示す。温度測定装置は、CPU(Central  Processing  Unit)63、記憶装置(記憶部)62およびインタフェース装置61を備えたコンピュータ60と、これらのハードウェア資源を制御するプログラムによって実現することができる。ここで、インタフェース装置に、測定部と、出力部が接続される。CPUは、記憶装置に格納された温度測定プログラムに従って本発明の実施の形態における処理を実行する。このように、温度測定プログラムは温度測定装置を機能させる。
<Computer configuration example>
FIG. 6 shows a configuration example of the computer 60 in the temperature measuring device according to the embodiment of the present invention. The temperature measuring device can be realized by a computer 60 including a CPU (Central Processing Unit) 63, a storage device (storage unit) 62, and an interface device 61, and a program for controlling these hardware resources. Here, the measurement unit and the output unit are connected to the interface device. The CPU executes the process according to the embodiment of the present invention according to the temperature measurement program stored in the storage device. In this way, the temperature measuring program makes the temperature measuring device work.
 本発明の実施の形態に係る温度測定装置では、コンピュータを装置内部に備えてもよいし、コンピュータの機能の少なくとも1部を外部コンピュータを用いて実現してもよい。また、記憶部も装置外部の記憶媒体64を用いてもよく、記憶媒体64に格納された温度測定プログラムを読み出して実行してもよい。記憶媒体64には、各種磁気記録媒体、光磁気記録媒体、CD-ROM、CD-R、各種メモリを含む。また、温度測定プログラムはインターネットなどの通信回線を介してコンピュータに供給されてもよい。 In the temperature measuring device according to the embodiment of the present invention, a computer may be provided inside the device, or at least one part of the functions of the computer may be realized by using an external computer. Further, the storage unit may also use the storage medium 64 outside the apparatus, or may read out and execute the temperature measurement program stored in the storage medium 64. The storage medium 64 includes various magnetic recording media, optical magnetic recording media, CD-ROMs, CD-Rs, and various memories. Further, the temperature measurement program may be supplied to the computer via a communication line such as the Internet.
 以上のように、本実施の形態に係る温度測定装置、方法およびプログラムによれば、生体内部の温度を非侵襲に精度良く測定できる。 As described above, according to the temperature measuring device, method and program according to the present embodiment, the temperature inside the living body can be measured non-invasively and accurately.
 本発明の実施の形態では、温度測定装置の構成および温度測定方法等において、各構成部の構造、寸法、材料等の一例を示したが、これに限らない。温度測定装置の構成および温度測定方法等の機能を発揮し効果を奏するものであればよい。 In the embodiment of the present invention, an example of the structure, dimensions, materials, etc. of each component is shown in the configuration of the temperature measuring device, the temperature measuring method, and the like, but the present invention is not limited to this. Anything that exerts functions such as the configuration of the temperature measuring device and the temperature measuring method and exerts an effect may be used.
  本発明は、作業者、競技者等の体温管理に用いる深部温度計に適用することができる。 The present invention can be applied to a deep thermometer used for body temperature control of workers, athletes and the like.
30  温度測定装置
31  測定部(センサ)
32  記憶部
321 温度等記憶部
322 推定温度記憶部
33  演算部
331 算出部
332 比較部
34  出力部 
30 Temperature measuring device 31 Measuring unit (sensor)
32 Storage unit 321 Temperature storage unit 322 Estimated temperature storage unit 33 Calculation unit 331 Calculation unit 332 Comparison unit 34 Output unit

Claims (8)

  1.  生体の内部の温度を、センサにより検出された温度を基に測定する温度測定方法であって、
     前記生体の内部の温度の推定温度を、前記センサにおける前記生体近傍の一の箇所の温度と他の箇所の温度との差分により算出される熱流束に基づく第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより算出し、
     前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と、前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と、前記第2の時刻での前記第2の式による第3の推定温度との差分を算出して、
     前記推定温度の算出式を、前記2つの差分に基づき、前記第1の式と前記第2の式とのいずれかに決定する温度測定方法。
    It is a temperature measurement method that measures the temperature inside a living body based on the temperature detected by a sensor.
    The first equation based on the heat flux calculated by the difference between the temperature of one place near the living body and the temperature of the other place in the sensor for the estimated temperature of the temperature inside the living body, and the first equation. Calculated by one of the second equations that correct the estimated temperature by the equation,
    The first estimated temperature at the first time before the time when the time derivative of the heat flux exceeds the reference value, and the second time after the time when the time derivative of the heat flux exceeds the reference value. Calculate the difference between the second estimated temperature according to the equation 1 and the difference between the first estimated temperature and the third estimated temperature according to the second equation at the second time.
    A temperature measuring method for determining the estimated temperature calculation formula as either the first formula or the second formula based on the difference between the two.
  2.  生体の内部の温度を、センサにより検出された温度を基に測定する温度測定方法であって、
     前記センサにおいて、前記生体に接触する面近傍の第1の温度を測定し、前記第1の温度が測定される位置と異なる位置で第2の温度を測定するステップと、
     前記第1の温度と前記第2の温度の差を基に熱流束を算出するステップと、
     前記熱流束の時間微分が基準値を超える時刻を検出するステップと、
     前記熱流束に基づき推定温度を算出する第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより前記推定温度を算出するステップと、
     前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と前記第2の時刻での前記第2の式による第3の推定温度との差分とを比較して、前記推定温度の算出式を決定するステップと
     を備える温度測定方法。
    It is a temperature measurement method that measures the temperature inside a living body based on the temperature detected by a sensor.
    In the sensor, a step of measuring a first temperature near a surface in contact with the living body and measuring a second temperature at a position different from the position where the first temperature is measured, and a step of measuring the second temperature.
    The step of calculating the heat flux based on the difference between the first temperature and the second temperature, and
    The step of detecting the time when the time derivative of the heat flux exceeds the reference value, and
    A step of calculating the estimated temperature by either the first equation for calculating the estimated temperature based on the heat flux and the second equation for correcting the estimated temperature by the first equation.
    The first estimated temperature at the first time before the time when the time differential of the heat flux exceeds the reference value and the first at the second time after the time when the time differential of the heat flux exceeds the reference value. The difference between the first estimated temperature and the third estimated temperature according to the second equation at the second time is compared with the difference from the second estimated temperature according to the above equation, and the estimated temperature is compared. A temperature measurement method with steps to determine the formula for.
  3.  前記第1の式と前記第2の式とのいずれかのうち、前記差分の小さい方を算出する式を、前記推定温度の算出式として決定するステップと
     を備える請求項2に記載の温度測定方法。
    The temperature measurement according to claim 2, further comprising a step of determining the formula for calculating the smaller difference between the first formula and the second formula as the calculation formula for the estimated temperature. Method.
  4.  前記第1の式が、
     Tcore = Tskin + Rbody × (Tskin-Tt)/Rs
    であり、前記第2の式が、
     Tcore = Tskin + α × Rbody × (Tskin-Tt)/Rs
    であり、Tskinが第1の温度であり、Ttが第2の温度であり、Rbodyは生体の熱抵抗であり、Rsはセンサの熱抵抗であり、αが補正係数である、
    請求項1から請求項3のいずれか一項に記載の温度測定方法。
    The first equation is
    Tcore = Tskin + Rbody × (Tskin-Tt) / Rs
    The second equation is
    Tcore = Tskin + α × Rbody × (Tskin-Tt) / Rs
    Tskin is the first temperature, Tt is the second temperature, Rbody is the thermal resistance of the living body, Rs is the thermal resistance of the sensor, and α is the correction coefficient.
    The temperature measuring method according to any one of claims 1 to 3.
  5.  前記第2の式における補正係数が1.03以上1.15以下であることを特徴とする請求項4に記載の温度測定方法。 The temperature measuring method according to claim 4, wherein the correction coefficient in the second equation is 1.03 or more and 1.15 or less.
  6.  前記熱流束の時間微分が0以上のときに対流の開始にともなう熱輸送形態の変化を判定し、前記熱流束の時間微分が0より小さいときに前記対流の終了にともなう熱輸送形態の変化を判定することを特徴とする請求項1から請求項5のいずれか一項に記載の温度測定方法。 When the time derivative of the heat flux is 0 or more, the change in the heat transport mode accompanying the start of convection is determined, and when the time derivative of the heat flux is smaller than 0, the change in the heat transport mode accompanying the end of the convection is determined. The temperature measuring method according to any one of claims 1 to 5, wherein the determination is made.
  7.  生体の内部の温度を、センサにより検出された温度を基に測定する温度測定装置であって、
     前記生体に接触する面近傍に配置され、第1の温度を測定する第1の温度計素子と、第1の温度計素子から離れた位置に配置され、第2の温度を測定する第2の温度計素子とを有する前記センサと、
     前記センサにおける前記生体近傍の一の箇所の温度と他の箇所の温度との差分により算出される熱流束に基づく第1の式と、前記第1の式による推定温度を補正する第2の式と、それぞれにより推定温度を算出する算出部と、
     前記差分の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度を記憶する推定温度記憶部と、
     前記第1の推定温度と、前記差分の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による推定温度と、前記第2の時刻での前記第2の式による推定温度とを比較する比較部と
     を備える温度測定装置。
    A temperature measuring device that measures the temperature inside a living body based on the temperature detected by a sensor.
    A first thermometer element arranged near the surface in contact with the living body and measuring the first temperature, and a second thermometer element arranged at a position away from the first thermometer element and measuring the second temperature. The sensor having a thermometer element and
    A first equation based on the heat flux calculated by the difference between the temperature at one location near the living body and the temperature at another location in the sensor, and a second equation for correcting the estimated temperature according to the first equation. And the calculation unit that calculates the estimated temperature by each,
    An estimated temperature storage unit that stores the first estimated temperature at the first time before the time when the time derivative of the difference exceeds the reference value.
    According to the first estimated temperature, the estimated temperature according to the first equation at the second time after the time when the time derivative of the difference exceeds the reference value, and the second equation at the second time. A temperature measuring device equipped with a comparison unit that compares with the estimated temperature.
  8.  生体の内部の温度を、センサにより検出された温度を基に測定する温度測定装置に対し、
     前記生体の内部の温度の推定温度を、前記センサにおける前記生体近傍の一の箇所の温度と他の箇所の温度との差分により算出される熱流束に基づく第1の式と、前記第1の式による推定温度を補正する第2の式とのいずれかにより算出し、
     前記熱流束の時間微分が基準値を超える時刻より前の第1の時刻での第1の推定温度と、前記熱流束の時間微分が基準値を超える時刻以降の第2の時刻での前記第1の式による第2の推定温度との差分と、前記第1の推定温度と、前記第2の時刻での前記第2の式による第3の推定温度との差分を算出して、
     前記推定温度の算出式を、前記2つの差分に基づき、前記第1の式と前記第2の式とのいずれかに決定する処理を実行させることを特徴とする、温度測定装置を機能させるための温度測定プログラム。
    For a temperature measuring device that measures the temperature inside a living body based on the temperature detected by a sensor.
    The first equation based on the heat flux calculated by the difference between the temperature of one place near the living body and the temperature of the other place in the sensor for the estimated temperature of the temperature inside the living body, and the first equation. Calculated by one of the second equations that correct the estimated temperature by the equation,
    The first estimated temperature at the first time before the time when the time derivative of the heat flux exceeds the reference value, and the second time after the time when the time derivative of the heat flux exceeds the reference value. Calculate the difference between the second estimated temperature according to the equation 1 and the difference between the first estimated temperature and the third estimated temperature according to the second equation at the second time.
    To function a temperature measuring device, which comprises executing a process of determining the estimated temperature calculation formula as one of the first formula and the second formula based on the difference between the two. Temperature measurement program.
PCT/JP2020/027224 2020-07-13 2020-07-13 Temperature measurement device, method, and program WO2022013911A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2022535991A JP7464124B2 (en) 2020-07-13 2020-07-13 Temperature measuring device, method and program
PCT/JP2020/027224 WO2022013911A1 (en) 2020-07-13 2020-07-13 Temperature measurement device, method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/027224 WO2022013911A1 (en) 2020-07-13 2020-07-13 Temperature measurement device, method, and program

Publications (1)

Publication Number Publication Date
WO2022013911A1 true WO2022013911A1 (en) 2022-01-20

Family

ID=79555234

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/027224 WO2022013911A1 (en) 2020-07-13 2020-07-13 Temperature measurement device, method, and program

Country Status (2)

Country Link
JP (1) JP7464124B2 (en)
WO (1) WO2022013911A1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019026439A1 (en) * 2017-07-31 2019-02-07 帝人株式会社 Deep body temperature estimation system, heat stress warning system, and deep body temperature estimation method
JP2019097819A (en) * 2017-11-30 2019-06-24 株式会社テクノ・コモンズ Biological data measuring apparatus
WO2019129469A1 (en) * 2017-12-29 2019-07-04 Medectis Ip Ltd. Non-invasive technique for body core temperature determination
JP2020008521A (en) * 2018-07-12 2020-01-16 オリンパス株式会社 Device and method for measuring temperature

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019026439A1 (en) * 2017-07-31 2019-02-07 帝人株式会社 Deep body temperature estimation system, heat stress warning system, and deep body temperature estimation method
JP2019097819A (en) * 2017-11-30 2019-06-24 株式会社テクノ・コモンズ Biological data measuring apparatus
WO2019129469A1 (en) * 2017-12-29 2019-07-04 Medectis Ip Ltd. Non-invasive technique for body core temperature determination
JP2020008521A (en) * 2018-07-12 2020-01-16 オリンパス株式会社 Device and method for measuring temperature

Also Published As

Publication number Publication date
JP7464124B2 (en) 2024-04-09
JPWO2022013911A1 (en) 2022-01-20

Similar Documents

Publication Publication Date Title
JP4658818B2 (en) Temperature estimation method and apparatus
US10890892B2 (en) Abnormality determination apparatus, non-transitory computer readable medium encoded with a program, abnormality determination system and abnormality determination method
JP2009505062A5 (en)
JP5139666B2 (en) Brightness adjusting device and brightness adjusting method
CN103234647B (en) Temperature correction method and temperature correction system of embedded system
JP2001324305A (en) Image correspondent position detector and range finder equipped with the same
US9702681B2 (en) System and method for temperature compensation of measurement machine
JPWO2021166128A5 (en) Abnormality estimation device, abnormality estimation method, and program
US12088784B2 (en) Thermal imaging camera and shutter operation state monitoring method thereof
JP2005516222A (en) Method and system for correcting and calibrating accelerometers with bias instability
WO2022013911A1 (en) Temperature measurement device, method, and program
KR101262277B1 (en) Detection methode for collision of robot
CN111537076B (en) Method and system for inhibiting temperature drift of infrared equipment in starting stage
WO2022013912A1 (en) Temperature measurement device, method, and program
KR102212113B1 (en) Measurement and Correction of Bio Environmental Information in Electronic Fabric Structure
WO2020245909A1 (en) Temperature measurement method and program
TW202122826A (en) Distance estimation device and method thereof and signal power calibration method
JP2008151739A (en) Temperature estimation method and device
RU2719310C2 (en) System and method for measurement of path length using pocket electronic device
WO2021260834A1 (en) Temperature measurement device, method, and program
JP2009248209A (en) Method of estimating thermal displacement of machine tool
JP7473074B2 (en) TEMPERATURE ESTIMATION METHOD, TEMPERATURE ESTIMATION PROGRAM, AND TEMPERATURE ESTIMATION DEVICE
JP2012083134A (en) Epicenter distance estimation method with improved immediate responsiveness
TWI484142B (en) A multi-sensing element correction system, a correction method and a recording medium
JP6945935B2 (en) Body temperature measurement processing program, body temperature measurement device and body temperature measurement processing system equipped with this program

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20945454

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022535991

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20945454

Country of ref document: EP

Kind code of ref document: A1