The present invention relates to a core body temperature estimation device, a core body temperature estimation method, and a core body temperature estimation program.
In recent years, the tendency toward high temperature has become remarkable and the risk of heatstroke has increased. Accordingly, various techniques for heatstroke countermeasures have been disclosed (see, for example, PTL 1 and PTL 2). The technique described in PTL 1 relates to air-conditioned clothes having a body temperature monitoring function for the purpose of comprehensively managing the physical conditions of a worker, in addition to preventing an increase in body temperature in a hot environment.
The air-conditioned clothes having the body temperature monitoring function include a clothing portion worn by the worker, an air blowing device that blows air to the space between the clothing portion and the body of the worker, a body temperature measurement device that measures the body temperature of the worker wearing the clothing portion, and an air blowing control device that controls the driving of the air blowing device, based on the measurement result of the body temperature by the body temperature measurement device. When the measured body temperature is lower than a predetermined temperature, the air blowing control device stops air blowing by the air blowing device or reduces the amount of air blown. When the measured body temperature is higher than the predetermined temperature, the air blowing control device increases the amount of air blown by the air blowing device. A heart rate measurement device for measuring the heart rate of the worker is further included.
The technique described in PTL 2 detects the risk of heatstroke of a worker before an abnormality occurs in the core body temperature of the worker. This technique includes workload information input means for inputting the workload information of a worker, a respiration sensor for detecting the respiration of the worker in real time, an electrocardiograph for detecting the heart rate of the worker in real time, information recording means for recording respiration curve information acquired by the respiration sensor and heart rate information acquired by the electrocardiograph, with the respiration curve information being in sync with the heart rate information, RSA calculation means for deriving respiratory sinus arrhythmia (RSA), based on the respiration curve information and the heart rate information recorded in the information recording means, and display means for displaying an RSA calculation value calculate by the RSA calculation means and the workload information.
[PTL 1] Japanese Patent Application Publication No. 2017-166075
[PTL 2] Japanese Patent Application Publication No. 2017-27123
Heatstroke is considered to be associated with core body temperature. Accordingly, the core body temperature of a worker is desirably measured to determine the risk of heatstroke of the worker. However, it is difficult to measure the core body temperature of the worker during labor work. Therefore, the core body temperature needs to be estimated based on physiological indexes that can be measured more easily than the core body temperature.
An object of the present invention is to provide a technique for estimating a core body temperature, based on physiologicalindexes.
The following means are adopted to solve the problems described above.
That is, according to a first aspect, there is provided a core body temperature estimation device including: a core body temperature acquisition unit that acquires a core body temperature, at a first time, of an estimation target subject targeted for core body temperature estimation; a beat acquisition unit that acquires a beat, in a period including the first time, of the estimation target subject; and an estimation unit that estimates the core body temperature of the estimation target subject based on the core body temperature at the first time acquired by the core body temperature acquisition unit and the beat in the period including the first time acquired by the beat acquisition unit, by using a core body temperature estimation model that estimates a core body temperature at a predetermined time based on a core body temperature at an initial time and a beat from the initial time to the predetermined time.
The disclosed aspect may be achieved by causing an information processing device to execute a program. That is, the disclosed structure can be identified as a program that causes an information processing device to execute the processes performed by individual means of the aspect described above or as a computer readable recording medium that records the program. In addition, the disclosed structure may be identified as a method of executing processes executed by the individual means described above. The disclosed structure may be identified as a system including the information processing device that executes the processes executed by the individual means described above.
According to the present invention, it is possible to provide a technique for estimating core body temperature based on physiological indexes.
An embodiment will be described below with reference to the drawings. The structure of the embodiment is an example and the structure of the invention is not limited to the specific structure of the disclosed embodiment. In carrying out the invention, the specific structure according to an embodiment may be adopted as appropriate.
The core body temperature estimation model building device 100 acquires the core body temperature and heart rate information when a plurality of model building subjects are working under predetermined stresses and builds a core body temperature estimation model for estimating the core body temperature.
The core body temperature acquisition unit 102 acquires the core body temperature of each of the model building subjects measured by thermometer or the like for measuring the core body temperature. The core body temperature is, for example, rectal temperature.
The heart rate information acquisition unit 104 acquires the heart rate information of each of the model building subjects during work measured by an electrocardiograph or the like for measuring the heart rate information. The heart rate information is, for example, electrocardiogram data or heart rate data. The electrocardiograph is an example of a pulsometer that measures beats.
The estimation model building unit 106 builds a core body temperature estimation model that estimates the core body temperature by a regression analysis based on the core body temperature after the initial time and the heart rate information after the initial time of each of the model building subjects.
The storage unit 108 stores programs, models, data, and the like used in the core body temperature estimation model building device 100. The storage unit 108 stores the built core body temperature estimation model, the core body temperature, the heart rate information, and the like.
The core body temperature estimation device 200 acquires the core body temperature, at the initial time, of the estimation target subject targeted for core body temperature estimation and the heart rate information of the estimation target subject during work, and estimates the core body temperature of the estimation target subject using the core body temperature estimation model built by the core body temperature estimation model building device 100. The model building subjects and the estimation target subject are also collectively referred to simply as the subjects.
The initial core body temperature acquisition unit 202 acquires the core body temperature of the estimation target subject before work measured by a thermometer or the like for measuring the core body temperature.
The heart rate information acquisition unit 204 acquires the heart rate information of the estimation target subject during work measured by an electrocardiograph or the like for measuring the heart rate information. The heart rate information is, for example, electrocardiogram data or heart beat data. The heart rate information acquisition unit 204 is an example of the beat acquisition unit.
The estimation unit 206 estimates the change amount of the core body temperature in each of the reference periods using the core body temperature estimation model built by the core body temperature estimation model building device 100 based on the core body temperature at the first time (arbitrary time) of the estimation target subject and the heart rate information in the period including the first time.
The storage unit 208 stores programs, models, data, and the like used in the core body temperature estimation device 200. The storage unit 208 stores the core body temperature estimation model built by the core body temperature estimation model building device 100, the estimated change amount of the core body temperature, the core body temperature, the heart rate information, and the like.
Although heart rate information is used here, other beat information may be used instead of heart rate information. A heart rate is the beat of a heart. A heart rate is one type of a beat. A beat is a motion that occurs when an internal organ is repeatedly contracted and relaxed. Other examples of a beat include a fingertip beat wave and an earlobe beat wave. A heart rate is the number of heart beats in a certain period. A pulse beat is the beat of an artery. A beat interval is one cycle of a beat.
The core body temperature estimation model building device 100 and the core body temperature estimation device 200 can be achieved by use of a dedicated or general-purpose computer such as a personal computer (PC), a work station (WS), a smartphone, a mobile phone, a tablet terminal, a car navigation device, and a personal digital assistant (PDA), or an electronic device having a computer.
The information processing device 90 can achieve the function that meets a predetermined purpose by causing the processor 91 to load a program stored in the recording media into the work area of the memory 92 and execute the program and the components and the like to be controlled through the execution of the program.
The processor 91 is, for example, a central processing unit (CPU) or a digital signal processor (DSP).
The memory 92 includes, for example, a random access memory (RAM) and a read only memory (ROM). The memory 92 is also referred to as a main storage device.
The storage unit 93 is, for example, an erasable programmable ROM (EPROM) or a hard disk drive (HDD). In addition, the storage unit 93 may include a removable medium, that is, a portable recording medium. The removable medium is, for example, a universal serial bus (USB) memory or a disc recording medium such as a compact disc (CD) or a digital versatile disc (DVD). The storage unit 93 is also referred to as a secondary storage device.
The storage unit 93 stores various types of programs, various types of data, and various types of tables in the recording medium in a readable and writable manner. The storage unit 93 stores an operating system (OS), various types of programs, various types of tables, and the like. The information stored in the storage unit 93 may be stored in the memory 92. In addition, the information stored in the memory 92 may be stored in the storage unit 93.
The operating system is the software that mediates between software and hardware, manages memory space, manages files, manages processes and tasks, and so on. The operating system includes a communication interface. The communication interface is the program that exchanges data with other external devices and the like connected via the communication control unit 96. The external devices and the like include, for example, other computers, external storage devices, and the like.
The input unit 94 includes a keyboard, a pointing device, a wireless remote controller, a touch panel, and the like. In addition, the input unit 94 may include an input device for video or images such as a camera or an input device for audio such as a microphone.
The output unit 95 includes a display device such as a liquid crystal display (LCD), an EL (electroluminescence) panel, a cathode ray tube (CRT) display, a plasma display panel (PDP), and an output device such as a printer. In addition, the output unit 95 may include an output device for audio such as a speaker.
The communication control unit 96 is connected to another device and controls communication between the computer 90 and the other device. The communication control unit 96 is, for example, a local area network (LAN) interface board, a wireless communication circuit for wireless communication, and a communication circuit for wired communication. The LAN interface board and the wireless communication circuit are connected to a network such as the Internet.
The computer that achieves the core body temperature estimation model building device 100 achieves the functions of the core body temperature acquisition unit 102, the heart rate information acquisition unit 104, and the estimation model building unit 106 by causing the processor to load programs stored in an auxiliary storage device into the main storage device and execute the programs. On the other hand, the storage unit 108 is provided in the storage area of the main storage device or the auxiliary storage device.
The computer that achieves the core body temperature estimation device 200 achieves the functions of the initial core body temperature acquisition unit 202, the heart rate information acquisition unit 204, and the estimation unit 206 by causing the processor to load programs stored in the auxiliary storage device into the main storage device and execute the programs. On the other hand, the storage unit 208 is provided in the storage area of the main storage device or the auxiliary storage device.
In S101, the core body temperature acquisition unit 102 of the core body temperature estimation model building device 100 acquires the core body temperature of each of the model building subjects measured by a thermometer or the like for measuring the core body temperature. The core body temperature acquisition unit 102 acquires the core body temperature from the start time to the end time of the exercise stress test of each of the plurality of model building subjects through a clinical thermometer or the like connected to the core body temperature estimation model building device 100 directly or via a network or the like. In addition, the core body temperature acquisition unit 102 may cause the model building subjects or other managers to input the core body temperature through input means such as a keyboard so as to acquire the core body temperature from the start time to the end time of the exercise stress test. The core body temperature acquisition unit 102 stores the acquired core body temperatures in the storage unit 108.
In S102, the heart rate information acquisition unit 104 acquires the heart rate information from the start time to the end time of the exercise stress test of each of the plurality of model building subjects measured by an electrocardiograph or the like for measuring the heart rate information. The electrocardiograph is, for example, a device including electrodes attached to the skin of a model building subject and a transmission unit for storing and transmitting heart rate information measured by the electrodes. In addition, the electrocardiograph may be a wearable terminal including electrodes provided on a wear (such as a T-shirt) and a transmission unit for transmitting the electrocardiogram data measured by the electrodes. In addition, the electrocardiograph may also be a wearable terminal worn on a wrist for measurement such as a wristwatch or a wristband that acquires heart rate information, or a clip-type wearable terminal worn on a fingertip or an ear. The heart rate information acquisition unit 104 acquires the heart rate information from the start time to the end time of the exercise stress test of each of the model building subjects through an electrocardiograph or the like connected to the core body temperature estimation model building device 100 directly or via a network or the like. In addition, the heart rate information acquisition unit 104 may acquire the heart rate information via a network or the like from another information processing device or the like that has acquired the heart rate information from the electrocardiograph or the like. The heart rate information acquisition unit 104 acquires the heart rate information after the initial time. The heart rate information acquisition unit 104 stores the acquired heart rate information of the model building subjects in the storage unit 108.
The processing steps S101 and S102 may be exchanged with each other. In addition, the processing steps S101 and S102 may be performed in parallel. The core body temperature acquisition unit 102 and the heart rate information acquisition unit 104 may acquire the core body temperature and the heart rate information in parallel with the measurement of the core body temperature and the heart rate information during the exercise stress test of the model building subjects. The core body temperature and the heart rate information are stored in the storage unit 108 together with the time information indicating the measurement time.
In S103, the estimation model building unit 106 builds a core body temperature estimation model that estimates the core body temperature by a regression analysis based on the acquired core body temperature and heart rate information after the initial time of each of the model building subjects. For example, the core body temperature estimation model estimates the change amount of the core body temperature in each of the predetermined reference period at predetermined reference periods from the initial time. The core body temperature estimation model is built using, for example, the heart rate information and the core body temperature at the start of one reference period as explanatory variables and the change amount of the core body temperature in the reference period as an objective variable. The objective variable may be an estimated value of core body temperature. The core body temperature estimation model may be, for example, a model that estimates the core body temperature for each of the reference periods from the initial time.
The estimation model building unit 106 calculates the R-R interval (RRI) based on the acquired heart rate information. The RRI is the time difference between the time of occurrence of an R wave and the time of occurrence of the preceding R wave on an electrocardiogram. The RRI is has a length of one heart rate (one cycle). The RRI is an example of a beat interval. The variability in the RRI is referred to as heart rate variability (HRV). The estimation unit 106 calculates the RRI based on the acquired heart rate information. The estimation model building unit 106 divides the calculated RRI of each heart rate for each of the reference periods (for example, every three minutes) from the initial time and creates a Poincaré plot for each of the reference periods. The length of one heart rate calculated by another well-known method may be used instead of the RRI. The reference period is not limited to three minutes and may be longer or shorter than three minutes. The reference period is preferably one minute or more from a statistical viewpoint for calculating an index that reflects the heart rate variability, which will be described later. In addition, the reference period is preferably 15 minutes or shorter from the viewpoint of the accuracy of estimation of the core body temperature and the operation.
The Poincaré plot in
Where, N is the total number of points of one Poincaré plot. The indexes SD1 and SD2 reflect the heart rate variability. It should be noted that SD1 and SD2 are 0 or positive values. SD1 and SD2 are examples of indicators acquired by the analysis of the Poincaré plot of the RRI.
The estimation model building unit 106 calculates SD2 for each of the reference periods based on the RRI. SD2(x) represents SD2 of the x-th reference period (x-th reference period). The estimation unit 106 builds a core body temperature estimation model for estimating the core body temperature for each of reference periods (for example, every three minutes) from the initial time. Here, when the reference period is P, the time after a lapse of the first reference period from the initial time (t=0) is t=1×P and the time after a lapse of the j-th reference period is t=jP. That is, the first reference period begins with the initial time (t=0) and ends with the time (t=1×P) after a lapse of the first reference period. The j-th reference period begins with the time (t=(j−1)×P) after a lapse of the (j−1)-th reference period and ends with the time (t=j×P) after a lapse of the j-th reference period. The estimation model building unit 106 builds a core body temperature estimation model by performing a regression analysis that uses, as the explanatory variables, SD2(j) and SD2(j−1) acquired by analysis of the Poincaré plot of the RRIs of the j-th reference period and the (j−1)-th reference period and the estimated value (which is assumed to be TMP(j−1)) of the core body temperature at time t=(j−1)P after a lapse of the (j−1)-th reference period (at the start of the j-th reference period) and uses, as the objective variable, the change amount (change amount of the core body temperature in the j-th reference period) of the core body temperature from time t=(j−1)P to time t=jP. A statistical analysis method such as a multi-regression analysis or a logistic regression analysis can be used as the regression analysis. When the change amount (change amount of the core body temperature in the j-th reference period) of the core body temperature from time t=(j−1)P to time t=jP is ΔTMP(j), ΔTMP(j) is expressed as follows.
ΔTMP(j)=ƒ(SD2(j),SD2(j−1),TMP(j−1)) [Math. 2]
Where, f(SD2(j), SD2(j−1), TMP(j−1)) is a function of SD2(j), SD2(j−1), and TMP(j−1). In addition, TMP(0) is the core body temperature at the initial time. For example, f(SD2(j), SD2(j−1), TMP(j−1)) is expressed as follows.
ƒ(SD2(j),SD2(j−1),TMP(j−1))=A1+A2×log SD2(j)+A3×SD2(j)/SD2(j−1)+A4×TMP(j−1) [Math. 3]
Where, A1, A2, A3, and A4 are coefficients acquired by a regression analysis. In addition, log is common logarithm. The formula acquired by a regression analysis is not limited to the one illustrated here. In addition, the estimated value TMP(j) of the core body temperature at time t=jP is acquired as follows.
Where, TMP(0) is the core body temperature at the initial time (t=0) acquired by measurement. The estimated value TMP(j) of the core body temperature at time t=jP is acquired by adding the change amount of the core body temperature in each of the reference periods from time t=0 to time t=jP to the core body temperature TMP(0) at the initial time. The estimation model building unit 106 stores the built core body temperature estimation model (ΔTMP, TMP) in the storage unit 108. As a result, the core body temperature estimation model is built by the core body temperature estimation model building device 100. In the core body temperature estimation model, (ΔTMP, TMP) is expressed as a function of TMP(0) and SD2. Since the core body temperature estimation model includes SD2, the core body temperature can be calculated based on the heart rate variability.
The estimation model building unit 106 may use other physiological indexes as the explanatory variables in building the core body temperature estimation model. Other physiological indexes may be, for example, the gender, the age, the skin temperature, the blood pressure, the exhaled gas, the pulse wave, the blood oxygen concentration, the blood flow amount, the exercise volume, and the respiratory rate. Furthermore, the temperature inside clothes, the wet bulb globe temperature (WBGT), meteorological data, and the like may be added to the explanatory variables. More accurate estimation can be performed by building a core body temperature estimation model in consideration of these explanatory variables.
In S201, the initial core body temperature acquisition unit 202 of the core body temperature estimation device 200 acquires the core body temperature, at the initial time, of the estimation target subject measured by a thermometer or the like for measuring the core body temperature. The initial time is the start time of the period in which the core body temperature is estimated. The initial core body temperature acquisition unit 202 acquires the core body temperature, at the initial time, of the estimation target subject from a thermometer or the like connected to the core body temperature estimation device 200 directly or via a network or the like. In addition, the initial core body temperature acquisition unit 202 may acquire the core body temperature by causing the estimation target subject, another manager, or the like to input the core body temperature with input means such as a keyboard. The initial core body temperature acquisition unit 202 stores the acquired core body temperature in the storage unit 208. The initial core body temperature acquisition unit 202 acquires the core body temperature, at the initial time, of the estimation target subject and does not acquire the core body temperature after the initial time. When it is difficult to measure the core body temperature, at the initial time, of the estimation target subject, the core body temperature estimated based on the eardrum temperature or the skin surface temperature, at the initial time, of the estimation target subject may be used as the core body temperature at the initial time.
In S202, the heart rate information acquisition unit 204 acquires the heart rate information of the estimation target subject measured by an electrocardiograph or the like for measuring the heart rate information. The electrocardiograph may be similar to the electrocardiograph used in S102 in
The core body temperature and the heart rate information are stored in the storage unit 208 together with the time information indicating the measured time.
In S203, the estimation unit 206 estimates the core body temperature using the core body temperature estimation model built by the operation flow in
The estimation unit 206 extracts the core body temperature estimation model stored in the storage unit 208. In addition, the estimation unit 206 calculates the RRI of each heart rate based on the acquired heart rate information. The estimation unit 206 divides the calculated RRI of each heart rate for each of the reference periods (for example, every three minutes) from the time one reference period before the initial time and creates a Poincaré plot for each of the reference periods. The estimation unit 206 calculates the index SD2 for each of the reference periods based on the Poincaré plot of RRI. The estimation unit 206 the change amount (ΔTMP) of the core body temperature in each of the reference periods using the core body temperature estimation model based on the core body temperature at the initial time (arbitrary time) of the estimation target subject and the index SD2 for each of the reference periods. In addition, the estimation unit 206 estimates the core body temperature (TMP) in each of the reference periods (end time of each of the reference periods) of the estimation target subject by adding the change amount of the core body temperature in each of the reference periods to the core body temperature at the initial time. The estimated core body temperature and the like are stored in the storage unit 208.
In addition, the estimation unit 206 may warn of the risk of heatstroke or the like when the estimated core body temperature becomes a predetermined threshold or higher. The estimation unit 206 may output a warning of the risk of heatstroke or the like by using output means such as a display or a speaker. This enables the estimation target subject or the like to grasp the risk of heatstroke and the like. The estimation unit 206 may simply output the estimated core body temperature using the output means. This enables the estimation target subject or the like to grasp the estimated value of the core body temperature based on the heart rate information.
The estimation unit 206 may estimate the core body temperature for each of the reference periods by calculating SD2 based on the RRI. Furthermore, the estimation unit 206 may output the core body temperature estimated for each of the reference periods using the output means. This enables the estimation target subject or the like to grasp the core body temperature in real time by the output means of the core body temperature estimation device 200.
The core body temperature estimation device 200 may be integrated with the core body temperature estimation model building device 100 that builds the core body temperature estimation model of the operation flow in
In the core body temperature estimation device 200, the calculation of ΔTMP(1) or the like requires SD2(0), which is calculated based on the heart rate information from the time (t=−P) time P before the initial time to the initial time (time t=0). For example, the core body temperature TMP(1) at time t=P is measured with a thermometer or the like (or the core body temperature TMP(1) at time t=P is estimated based on the eardrum temperature and the skin surface temperature at time t=P) and acquired, and then ΔTMP(j) and TMP(j) after the time t=2×P(j=2) may be calculated. In addition, by considering the core body temperature TMP(1) at time t=P to be the same as the core body temperature TMP(0) at time t=0, TMP(1) may be the same as TMP(0). That is, the heart rate information from time t=−P to time t=0 is not necessary.
The core body temperature estimation model building device 100 acquires the core body temperatures (such as, for example, the rectal temperatures) and the heart rate information (such as, for example, the electrocardiogram data) of the plurality of model building subjects after the initial time. The core body temperature estimation model building device 100 calculates the RRI (or other beat intervals) based on the heart rate information and creates a Poincaré plot of the RRI for each of the reference periods. The core body temperature estimation model building device 100 calculates SD2 for each of the reference periods based on the RRI. The core body temperature estimation model building device 100 performs a regression analysis that uses SD2(j) and SD2(j−1) acquired from the analysis of the Poincare plot of the RRI and the estimated value TMP(j−1) of the estimated core body temperature as the explanatory variables and the estimated change amount ΔTMP(j) of the core body temperature as the objective variable and builds a core body temperature estimation model. The estimated value TMP(j) of the core body temperature at time t=jP is expressed as TMP(j−1)+ΔTMP(j). In addition, TMP(0) is the core body temperature at the initial time. The core body temperature estimation model building device 100 can build a core body temperature estimation model that estimates the core body temperature based on the core body temperature and the heart rate information of the model building subject after the initial time.
The core body temperature estimation device 200 acquires the core body temperature of the model building subject at the initial time. The core body temperature estimation device 200 acquires the heart rate information in the period including the initial time of the estimation target subject. The core body temperature estimation device 200 estimates the core body temperature of the estimation target subject based on the core body temperature at the initial time and the heart rate information in the period including the initial time using the core body temperature estimation model. The core body temperature estimation device 200 can easily estimate the core body temperature, which is difficult to constantly measure during work or the like, by using a physiological index such as the heart rate information with a core body temperature estimation model.
A core body temperature estimation model has been built by the core body temperature estimation model building device 100 using the core body temperature and heart rate information measured when the model building subject has performed the exercise stress test described above. The estimated value of the core body temperature has been calculated by the core body temperature estimation device 200 based on the built core body temperature estimation model using the core body temperature at the initial time and the measured heart rate information. The error between the measured value of the core body temperature and the estimated value of the core body temperature based on the core body temperature estimation model was −0.007° C. on average (−0.50° C. at maximum) and the average error rate is −0.02% (−1.30% at maximum). The core body temperature can be estimated accurately using the core body temperature and the heart rate information at the initial time based on the core body temperature estimation model.
A program that enables a computer or other machine or device (referred to below as a computer etc.) to achieve any of the above functions can be recorded in a computer etc. readable recording medium. Then, the function can be provided by causing the computer etc. to read and execute the program in this recording medium.
Here, a computer etc. readable recording medium is a recording medium that can store information such as data and programs by electrical, magnetic, optical, mechanical, or chemical action and can be read by a computer etc. Such a recording medium may be provided with components constituting a computer, such as a CPU and a memory so that the CPU executes programs.
In addition, such recording media that can be removed from a computer etc. include, for example, a flexible disk, a magneto-optical disk, a CD-ROM, a CD-R/W, a DVD, a DAT, an 8-mm tape, a memory card, and the like.
In addition, recording media that are fixed to a computer etc. include a hard disk drive, a ROM, an SSD, and the like.
Although an embodiment of the present invention has been described above, the embodiment is only an example, the present invention is not limited the embodiment, and various changes such as combinations of components can be made based on the knowledge of those skilled in the art without departing from the concept of the claims.
Number | Date | Country | Kind |
---|---|---|---|
2019-164711 | Sep 2019 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2020/034305 | 9/10/2020 | WO |