JP2013013542A - Apparatus and program for estimating state of consciousness - Google Patents
Apparatus and program for estimating state of consciousness Download PDFInfo
- Publication number
- JP2013013542A JP2013013542A JP2011148101A JP2011148101A JP2013013542A JP 2013013542 A JP2013013542 A JP 2013013542A JP 2011148101 A JP2011148101 A JP 2011148101A JP 2011148101 A JP2011148101 A JP 2011148101A JP 2013013542 A JP2013013542 A JP 2013013542A
- Authority
- JP
- Japan
- Prior art keywords
- consciousness
- series data
- state
- breathing
- estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 claims abstract description 93
- 230000003595 spectral effect Effects 0.000 claims abstract description 30
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 193
- 230000000241 respiratory effect Effects 0.000 claims description 105
- 238000001228 spectrum Methods 0.000 claims description 47
- 230000007423 decrease Effects 0.000 claims description 31
- 230000009467 reduction Effects 0.000 claims description 27
- 230000008859 change Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 27
- 206010012373 Depressed level of consciousness Diseases 0.000 abstract 2
- 238000000034 method Methods 0.000 description 19
- 238000005311 autocorrelation function Methods 0.000 description 15
- 230000000694 effects Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 230000009471 action Effects 0.000 description 6
- 238000012952 Resampling Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 230000008035 nerve activity Effects 0.000 description 3
- 210000001002 parasympathetic nervous system Anatomy 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 206010062519 Poor quality sleep Diseases 0.000 description 2
- 230000037007 arousal Effects 0.000 description 2
- 210000000467 autonomic pathway Anatomy 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000001121 heart beat frequency Effects 0.000 description 2
- 230000001734 parasympathetic effect Effects 0.000 description 2
- 230000036387 respiratory rate Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000008667 sleep stage Effects 0.000 description 2
- 208000003443 Unconsciousness Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000035565 breathing frequency Effects 0.000 description 1
- 230000002612 cardiopulmonary effect Effects 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 210000003016 hypothalamus Anatomy 0.000 description 1
- 238000002847 impedance measurement Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 210000005037 parasympathetic nerve Anatomy 0.000 description 1
- 210000001034 respiratory center Anatomy 0.000 description 1
- 238000009531 respiratory rate measurement Methods 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000020509 sex determination Effects 0.000 description 1
- 230000003860 sleep quality Effects 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
本発明は、意識状態推定装置及びプログラムに係り、特に、推定対象者の意識低下状態を推定する意識状態推定装置及びプログラムに関する。 The present invention relates to a consciousness state estimation device and a program, and more particularly, to a consciousness state estimation device and a program that estimate a consciousness reduction state of an estimation target person.
従来より、心拍間隔変動の低周波成分(0.04−0.15Hz)及び高周波成分(0.15−0.40Hz)と関連が深い、心拍間隔変動の長期標準偏差及び短期標準偏差を指標として、睡眠状態を推定する睡眠状態推定装置が知られている(特許文献1)。 Conventionally, the long-term standard deviation and short-term standard deviation of heartbeat interval variation, which are closely related to the low-frequency component (0.04-0.15 Hz) and high-frequency component (0.15-0.40 Hz) of heartbeat interval variation, are used as indicators. A sleep state estimation apparatus that estimates a sleep state is known (Patent Document 1).
また、呼吸数、呼吸数のばらつき、又は呼吸周期のばらつきを判定基準のひとつとする、睡眠段階判定装置が知られている(特許文献2)。 Further, a sleep stage determination device is known that uses a respiratory rate, a variation in the respiratory rate, or a variation in the respiratory cycle as one of the determination criteria (Patent Document 2).
また、心電図から求められたR−R間隔(≒心拍間隔)の時系列データと心電図上のR波の振幅から算出した代理呼吸曲線(≒呼吸波形)の時系列データについて、両データのリサンプリングデータのコヒーレンスとクロススペクトルを算出し、それらの指標より心肺カップリング状態を推定し、睡眠の質及び睡眠呼吸障害の評価を行う技術が知られている(特許文献3)。 In addition, resampling of both data is performed on time series data of RR interval (≈ heartbeat interval) obtained from the electrocardiogram and proxy respiratory curve (≈ breath waveform) calculated from the amplitude of the R wave on the electrocardiogram. There is known a technique for calculating coherence and cross spectrum of data, estimating a cardiopulmonary coupling state from these indexes, and evaluating sleep quality and sleep respiratory disorder (Patent Document 3).
また、呼吸変動に対する心拍変動の位相遅れを算出し、位相遅れが大きくなるほど睡眠深度が深いと判定する睡眠深度判定装置が知られている(特許文献4)。 Further, a sleep depth determination device that calculates a phase delay of heart rate variability with respect to respiratory variation and determines that the sleep depth is deeper as the phase lag increases is known (Patent Document 4).
心拍変動の低周波成分は副交感神経活動を反映し、心拍変動の高周波成分は副交感神経活動を反映していることは古くから知られており、睡眠状態の指標として用いられることが検討されている。しかしながら、心拍変動は副交感神経活動以外の要因(日内変動等)でも変動することが知られており、上記の特許文献1に記載の技術では、睡眠状態を正確に計測することができない、という問題がある。 It has long been known that the low-frequency component of heart rate variability reflects parasympathetic activity, and the high-frequency component of heart rate variability reflects parasympathetic activity, and it has been studied to be used as an indicator of sleep state . However, it is known that heart rate variability varies due to factors other than parasympathetic nerve activity (such as diurnal variation), and the technique described in Patent Document 1 cannot accurately measure the sleep state. There is.
また、上記の特許文献2に記載の技術では、呼吸数又は呼吸周期を検出するために呼吸ピークを明確化する必要があり、形状が不安定になりやすい呼吸波形では検出のロバスト性が問題となる。この問題は、不安定な非接触の呼吸計測系ではより大きな問題となる。 In the technique described in Patent Document 2, it is necessary to clarify the respiration peak in order to detect the respiration rate or the respiration cycle, and the robustness of detection is a problem in the respiration waveform that tends to be unstable in shape. Become. This problem becomes more serious in an unstable non-contact respiratory measurement system.
上記の特許文献3に記載の技術では、クロススペクトルを評価するとき、従来の心拍変動と同様に、心拍変動の超低周波成分、低周波成分、及び高周波成分と分けて評価する等、その主眼はあくまでも心拍変動解析にある。そのため、上記の特許文献3に記載の手法では、周波数分解能に対する要求が高いため、解析に使用する時間窓を非常に長く取る必要がある。したがって、睡眠判定の精度が、非常に時間分解能の低いものにならざるを得ず、リアルタイムに判定する必要がある用途には適さない、という問題がある。 In the technique described in the above-mentioned patent document 3, when evaluating the cross spectrum, as in the case of the conventional heart rate variability, the main focus is to separately evaluate the ultra-low frequency component, the low frequency component, and the high frequency component of the heart rate variability. Is in heart rate variability analysis. For this reason, in the method described in Patent Document 3, since there is a high demand for frequency resolution, it is necessary to take a very long time window for analysis. Therefore, there is a problem that the accuracy of sleep determination must be extremely low in time resolution, and is not suitable for applications that need to be determined in real time.
上記の特許文献4に記載の技術では、呼吸変動に対する心拍変動の位相遅れに着目しているが、位相遅れは呼吸及び心拍の検出手段によって異なるため、必ずしも精度良く計測することができない、という問題がある。 The technique described in Patent Document 4 focuses on the phase lag of heartbeat variability with respect to respiratory variability, but the problem is that phase lag varies depending on the breathing and heartbeat detection means, and cannot always be measured with high accuracy. There is.
本発明は、上記の問題点を解決するためになされたもので、推定対象者の意識低下状態を安定して精度良く推定することができる意識状態推定装置及びプログラムを提供することを目的とする。 The present invention has been made to solve the above-described problems, and an object of the present invention is to provide a consciousness state estimation device and a program that can stably and accurately estimate a consciousness reduction state of an estimation target person. .
上記の目的を達成するために第1の発明に係る意識状態推定装置は、推定対象者の心拍間隔の時系列データ及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段と、前記データ検出手段によって検出された前記心拍間隔の時系列データ及び前記呼吸情報の時系列データの相関性を判定する相関性判定手段と、前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、を含んで構成されている。 In order to achieve the above object, a consciousness state estimation device according to a first aspect of the present invention is a time series of respiratory information indicating time series data of heart rate intervals of an estimation target person and physical quantities that change according to the breath of the estimation target person. Data detection means for detecting data, correlation determination means for determining the correlation between the time series data of the heartbeat interval detected by the data detection means and the time series data of the respiratory information, and the correlation determination means Consciousness state estimation means for estimating the state of reduced consciousness of the person to be estimated is configured so that the consciousness state decreases as the determined correlation is higher.
第2の発明に係るプログラムは、コンピュータを、推定対象者の心拍間隔の時系列データ及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記心拍間隔の時系列データ及び前記呼吸情報の時系列データの相関性を判定する相関性判定手段、及び前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段として機能させるためのプログラムである。 According to a second aspect of the present invention, there is provided a program that causes a computer to detect time series data of a heart rate interval of an estimation target person and time series data of respiratory information indicating a physical quantity that changes in accordance with the estimation target person's breathing. Correlation determination means for determining the correlation between the detected time-series data of the heartbeat interval and the time-series data of the respiratory information, and the higher the correlation determined by the correlation determination means, the lower the consciousness state. And a program for functioning as a consciousness state estimating means for estimating a consciousness reduction state of the estimation target person.
第1の発明及び第2の発明によれば、データ検出手段によって、推定対象者の心拍間隔の時系列データ及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出する。相関性判定手段によって、前記データ検出手段によって検出された前記心拍間隔の時系列データ及び前記呼吸情報の時系列データの相関性を判定する。 According to the first and second aspects of the invention, the data detection means converts the time series data of the heart rate interval of the estimation target person and the time series data of the respiratory information indicating the physical quantity that changes according to the breathing of the estimation target person. To detect. Correlation determining means determines the correlation between the time series data of the heartbeat interval and the time series data of the respiratory information detected by the data detecting means.
そして、意識状態推定手段によって、前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する。 Then, the state of consciousness estimation of the estimation target person is estimated so that the state of consciousness decreases as the correlation determined by the correlation determination unit increases.
このように、検出された心拍間隔の時系列データ及び呼吸情報の時系列データの相関性を判定し、相関性が高いほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を安定して精度良く推定することができる。 In this way, the correlation between the detected time-series data of the heartbeat interval and the time-series data of the respiratory information is determined, and the consciousness state of the estimation target person is estimated so that the consciousness state decreases as the correlation increases. As a result, it is possible to stably and accurately estimate the state of reduced consciousness of the estimation target person.
第3の発明に係る意識状態推定装置は、推定対象者の心臓の拍動に応じて変化する信号を検出すると共に、前記検出された信号から、呼吸に由来する成分が含まれる予め定められた周波数帯域の信号を分離する信号検出手段と、前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段と、前記信号検出手段によって分離された前記周波数帯域の信号、及び前記データ検出手段によって検出された前記呼吸情報の時系列データの相関性を判定する相関性判定手段と、前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、を含んで構成されている。 A consciousness state estimation device according to a third aspect of the present invention detects a signal that changes in accordance with the heartbeat of the estimation target person, and includes a predetermined component that includes a component derived from respiration from the detected signal. A signal detection means for separating a signal in a frequency band; a data detection means for detecting time-series data of respiration information indicating a physical quantity that changes in accordance with the respiration of the person to be estimated; and the frequency separated by the signal detection means Correlation determination means for determining the correlation between the band signal and the time-series data of the respiratory information detected by the data detection means, and the higher the correlation determined by the correlation determination means, the lower the state of consciousness. As mentioned above, it includes consciousness state estimating means for estimating the state of consciousness reduction of the estimation target person.
第4の発明に係るプログラムは、コンピュータを、推定対象者の心臓の拍動に応じて変化する信号を検出すると共に、前記検出された信号から、呼吸に由来する成分が含まれる予め定められた周波数帯域の信号を分離する信号検出手段によって分離された前記周波数帯域の信号、及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記呼吸情報の時系列データの相関性を判定する相関性判定手段、及び前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段として機能させるためのプログラムである。 According to a fourth aspect of the invention, a computer detects a signal that changes in accordance with the heartbeat of the estimation target person, and includes a predetermined component that includes a component derived from respiration from the detected signal. Detected by data detection means for detecting time-series data of the respiratory information indicating the signal of the frequency band separated by the signal detection means for separating the signal of the frequency band and the physical quantity changing according to the breathing of the estimation target person. A correlation determination unit that determines the correlation of the time series data of the respiratory information, and a state of reduced consciousness of the estimation target person such that the higher the correlation determined by the correlation determination unit, the lower the consciousness state It is a program for functioning as a consciousness state estimation means for estimating
第3の発明及び第4の発明によれば、信号検出手段によって、推定対象者の心臓の拍動に応じて変化する信号を検出すると共に、前記検出された信号から、呼吸に由来する成分が含まれる予め定められた周波数帯域の信号を分離する。データ検出手段によって、前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出する。 According to 3rd invention and 4th invention, while the signal detection means detects the signal which changes according to the heartbeat of an estimation object person, the component derived from respiration is detected from the detected signal. A signal of a predetermined frequency band included is separated. The data detection means detects time-series data of respiration information indicating a physical quantity that changes in accordance with the respiration of the estimation target person.
そして、相関性判定手段によって、前記信号検出手段によって分離された前記周波数帯域の信号、及び前記データ検出手段によって検出された前記呼吸情報の時系列データの相関性を判定する。意識状態推定手段によって、前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する。 Then, the correlation determination unit determines the correlation between the signal in the frequency band separated by the signal detection unit and the time series data of the respiratory information detected by the data detection unit. The consciousness state estimation means estimates the consciousness reduction state of the estimation target person so that the consciousness state decreases as the correlation determined by the correlation determination means increases.
このように、推定対象者の心臓の拍動に応じて変化する信号から分離された、呼吸に由来する成分が含まれる周波数帯域の信号、及び呼吸情報の時系列データの相関性を判定し、相関性が高いほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を安定して精度良く推定することができる。 In this way, the correlation between the signal in the frequency band including the component derived from respiration, separated from the signal that changes according to the heart beat of the estimation target person, and the time-series data of the respiration information is determined, By estimating the consciousness reduction state of the estimation target person so that the consciousness state decreases as the correlation increases, the estimation target person's consciousness reduction state can be stably and accurately estimated.
第1の発明及び第3の発明に係る意識状態推定装置は、前記呼吸情報の時系列データに基づいて呼吸のばらつき度合いを判定する呼吸ばらつき判定手段を更に含み、前記意識状態推定手段は、前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いと、前記相関性判定手段によって判定された相関性とに基づいて、前記推定対象者の意識低下状態を推定するようにすることができる。これによって、呼吸のばらつき度合いを更に考慮して、推定対象者の意識低下状態をより精度良く推定することができる。 The consciousness state estimation device according to the first and third aspects of the present invention further includes a breathing variation determination unit that determines a degree of breathing variation based on the time series data of the breathing information, and the consciousness state estimation unit includes: Based on the degree of breathing variation determined by the breathing variation determination unit and the correlation determined by the correlation determination unit, it is possible to estimate the state of reduced consciousness of the estimation target person. Thereby, the degree of consciousness of the estimation target person can be estimated with higher accuracy in consideration of the degree of variation in breathing.
また、上記の呼吸ばらつき判定手段は、前記呼吸情報の時系列データに基づいて、前記呼吸のばらつき度合いとして、呼吸間隔の標準偏差又は前記呼吸間隔の変動係数を算出するようにすることができる。 Further, the breathing variation determining means can calculate a standard deviation of the breathing interval or a coefficient of variation of the breathing interval as the degree of breathing variation based on the time series data of the breathing information.
また、上記の呼吸ばらつき判定手段は、前記呼吸情報の時系列データに基づいて、前記呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定するようにすることができる。 In addition, the breathing variation determination unit may be configured such that, based on the time series data of the breathing information, the degree of breathing variation, the frequency dependence of the power spectrum density of the time series data, or the autocorrelation of the time series data. Can be determined.
第1の発明及び第3の発明に係る意識状態推定装置は、前記推定対象者の体動及び発話の少なくとも一方を検出する体動発話検出手段を更に含み、前記意識状態推定手段は、前記相関性判定手段によって判定された相関性と、前記体動発話検出手段による前記体動及び発話の少なくとも一方の検出結果とに基づいて、前記意識低下状態を推定するようにすることができる。 The consciousness state estimation device according to the first and third inventions further includes body movement utterance detection means for detecting at least one of the body movement and utterance of the person to be estimated, wherein the consciousness state estimation means includes the correlation The reduced consciousness state can be estimated based on the correlation determined by the sex determination unit and the detection result of at least one of the body movement and the utterance by the body movement utterance detection unit.
第5の発明に係る意識状態推定装置は、推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段と、前記データ検出手段によって検出された前記呼吸情報の時系列データに基づいて、呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する呼吸ばらつき判定手段と、前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いが低いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、を含んで構成されている。 A consciousness state estimation device according to a fifth aspect of the present invention is a data detection unit that detects time-series data of respiratory information indicating a physical quantity that changes in accordance with the breathing of the estimation target person, and the respiratory information detected by the data detection unit. Based on the time series data, the breathing variation determining means for determining the frequency dependence of the power spectrum density of the time series data or the autocorrelation of the time series data as the degree of breathing variation, and the breathing variation determining means Consciousness state estimation means for estimating a state of reduced consciousness of the person to be estimated so that the consciousness state is lowered as the degree of variation in breathing determined by the above is lower.
第6の発明に係るプログラムは、コンピュータを、推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記呼吸情報の時系列データに基づいて、呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する呼吸ばらつき判定手段、及び前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いが低いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段として機能させるためのプログラムである。 According to a sixth aspect of the invention, there is provided a program based on the time-series data of the respiratory information detected by the data detection means for detecting the time-series data of the respiratory information indicating the physical quantity that changes in accordance with the respiration of the estimation target person. Thus, as the degree of breathing variation, the breathing variation judging means for judging the frequency dependence of the power spectrum density of the time series data or the autocorrelation of the time series data, and the breathing judged by the breathing variation judging means This is a program for functioning as a consciousness state estimating means for estimating the consciousness-reduced state of the person to be estimated so that the consciousness state is lowered as the degree of variation in the level of the subject is reduced.
第5の発明及び第6の発明によれば、データ検出手段によって、推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出する。呼吸ばらつき判定手段によって、前記データ検出手段によって検出された前記呼吸情報の時系列データに基づいて、呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する。 According to the fifth and sixth inventions, the data detection means detects time-series data of respiration information indicating physical quantities that change in accordance with the respiration of the estimation target person. Based on the time series data of the respiratory information detected by the data detection means by the breathing variation determination means, the frequency dependence of the power spectrum density of the time series data, or the time series data Determine autocorrelation.
そして、意識状態推定手段によって、前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いが低いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する。 Then, the consciousness state estimation unit estimates the consciousness reduction state of the estimation target person so that the consciousness state decreases as the degree of variation in the breathing determined by the breathing variation determination unit decreases.
このように、呼吸のばらつき度合いとして、呼吸情報の時系列データのパワースペクトル密度の周波数依存性又は自己相関性を判定し、
呼吸のばらつき度合いが低いほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を安定して精度良く推定することができる。
Thus, as the degree of variation in breathing, determine the frequency dependence or autocorrelation of the power spectrum density of the time-series data of breathing information,
By estimating the consciousness reduction state of the estimation target person so that the consciousness state decreases as the degree of breathing variation decreases, the estimation target person's consciousness reduction state can be estimated stably and accurately.
第5の発明に係る意識状態推定装置は、前記推定対象者の体動及び発話の少なくとも一方を検出する体動発話検出手段を更に含み、前記意識状態推定手段は、前記呼吸ばらつき判定手段によって判定された前記呼吸のはらつき度合いと、前記体動発話検出手段による前記体動及び発話の少なくとも一方の検出結果とに基づいて、前記意識低下状態を推定するようにすることができる。 The consciousness state estimation device according to a fifth aspect of the present invention further includes body movement utterance detection means for detecting at least one of body motion and utterance of the person to be estimated, wherein the consciousness state estimation means is determined by the breathing variation determination means. The state of reduced consciousness can be estimated based on the degree of blurring of the breath and the detection result of at least one of the body motion and speech by the body motion speech detection means.
以上説明したように、本発明の意識状態推定装置及びプログラムによれば、推定対象者の意識低下状態を安定して精度良く推定することができる、という効果が得られる。 As described above, according to the consciousness state estimation apparatus and program of the present invention, there is an effect that it is possible to stably and accurately estimate the consciousness reduction state of the estimation target person.
以下、図面を参照して本発明の実施の形態を詳細に説明する。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
図1に示すように、第1の実施の形態に係る意識状態推定装置10は、推定対象者の心拍を検出する心拍検出装置12と、推定対象者の呼吸を検出する呼吸検出装置14と、心拍検出装置12から出力される心拍間隔の時系列データ及び呼吸検出装置14から出力される呼吸関連生体信号の時系列データに基づいて、推定対象者が意識低下状態であるか否かを推定し、推定結果を出力部18により出力させるコンピュータ16とを備えている。なお、心拍検出装置12及び呼吸検出装置14が、データ検出手段の一例である。
As shown in FIG. 1, the consciousness
心拍検出装置12は、心拍センサ20、心拍間隔算出部22、及び心拍間隔リサンプリング部24を備えている。心拍センサ20は、人体より心拍に関連する生体信号を検出するセンサであり、心拍関連生体信号の時系列データを出力する。心拍センサ20として、例えば、心電計、心音計、脈波センサ、電気インピーダンス計測装置、電磁誘導的に電気インピーダンスを計測する装置、圧電素子を用いて心拍由来の生体振動を検出する装置、又は電磁波により心拍由来の生体振動を検出する装置等を用いればよい。
The
心拍間隔算出部22は、心拍センサ20により出力された心拍関連生体信号の時系列データに基づいて、心拍時刻を算出し、心拍間隔の時系列データを出力する。心拍間隔リサンプリング部24は、心拍間隔算出部22から出力された心拍間隔の時系列データを、後述する呼吸センサ26のサンプリング時刻に合わせてリサンプリングして出力する。
The
呼吸検出装置14は、呼吸センサ26及び低域遮断フィルタ28により構成される。呼吸センサ26は、人体から呼吸に関連する生体信号(例えば、呼吸動を示す呼吸波形)を検出するセンサであり、呼吸関連生体信号の時系列データを出力する。呼吸センサ26として、例えば、胸郭周囲長を測定する装置、電気インピーダンス計測装置、電磁誘導的に電気インピーダンスを計測する装置、圧電素子を用いて呼吸由来の生体振動を検出する装置、又は電磁波により呼吸由来の生体振動を検出する装置等を用いればよい。また、心拍センサ20が心電図を出力する場合、心電図R波の振幅より呼吸波形を推定する手法を用いてもよい。なお、呼吸関連生体信号が、呼吸に応じて変化する物理量を示す呼吸情報の一例である。
The
低域遮断フィルタ28は、呼吸関連生体信号の時系列データより直流成分を除去し、直流成分を除去した呼吸関連生体信号の時系列データを出力する。
The low-
コンピュータ16は、CPUと、RAMと、後述する意識状態推定処理ルーチンを実行するためのプログラムを記憶したROMとを備え、機能的には次に示すように構成されている。コンピュータ16は、心拍検出装置12から出力される心拍間隔の時系列データ及び呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、心拍間隔の時系列データと呼吸関連生体信号の時系列データとの相関係数を算出する相関係数算出部32と、呼吸関連生体信号の時系列データに基づいて、呼吸関連生体信号のパワースペクトル密度を算出するパワースペクトル密度算出部34と、算出されたパワースペクトル密度に基づいて、パワースペクトル密度特徴量を算出するパワースペクトル密度特徴量算出部36と、算出された相関係数及びパワースペクトル密度特徴量に基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部38とを備えている。なお、相関係数算出部32が、相関性判定手段の一例である。パワースペクトル密度算出部34及びパワースペクトル密度特徴量算出部36が、呼吸ばらつき判定手段の一例である。
The
データ取得部30は、心拍検出装置12から入力された心拍間隔の時系列データを取得し、メモリ(図示省略)に記憶する共に、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得してメモリに記憶する。
The
ここで、本実施の形態の原理について説明する。 Here, the principle of the present embodiment will be described.
まず、人間の呼吸は、視床下部・橋などの呼吸中枢による無意識の制御を受けているが、大脳皮質より意識的に制御することが可能である。これにより、呼吸は意識的に深さや速さを変えたり、止めたりすることが可能となっている。ここで、意識状態が低下すると、意識による制御が弱まるため、呼吸は無意識の制御のみを受け、呼吸間隔が一定となると考えられる。 First, human breathing is unconsciously controlled by the respiratory center such as the hypothalamus and bridge, but can be controlled more consciously than the cerebral cortex. This makes it possible to consciously change or stop the depth and speed of breathing. Here, if the state of consciousness is lowered, control by consciousness is weakened, so that respiration is only subjected to unconscious control, and the breathing interval is considered to be constant.
また、呼吸は副交感神経系を介して心臓の制御にも影響を与えている。具体的には、肺での換気効率を上げるため、吸気時には心拍数を上げ、呼気時には心拍数を下げるように心臓を制御している。このように心拍間隔には呼吸との相関性が存在するが、心拍間隔は他の要因の影響も受けているため、呼吸によって心拍間隔が一意に決まることはない。しかし、図2に示すように、意識状態が低下している時には、人体に対する外部刺激の多くが遮断されるため、呼吸以外の心拍変動要因が減少し、呼吸波形と心拍間隔の相関性が強まると考えられる。 Respiration also affects the control of the heart via the parasympathetic nervous system. Specifically, in order to increase ventilation efficiency in the lungs, the heart is controlled so that the heart rate is increased during inspiration and the heart rate is decreased during expiration. As described above, there is a correlation between the heartbeat interval and respiration, but since the heartbeat interval is also influenced by other factors, the heartbeat interval is not uniquely determined by respiration. However, as shown in FIG. 2, when the state of consciousness is reduced, since many external stimuli to the human body are blocked, heart rate variability factors other than respiration are reduced, and the correlation between the respiration waveform and the heart rate interval is strengthened. it is conceivable that.
そこで、本実施の形態では、相関係数算出部32によって、リサンプリングされた心拍間隔の時系列データ及び直流成分を除去した呼吸関連生体信号の時系列データのうち、予め定められた時間長における、心拍間隔の時系列データ及び呼吸関連生体信号の時系列データ間の相関係数を、以下のように算出する。
Therefore, in the present embodiment, the correlation
まず、ある時間窓に含まれるN点の心拍間隔の時系列データをX(n)とすると、X(n)は、以下の(1)式のように表される。 First, when the time series data of the N-th heartbeat interval included in a certain time window is X (n), X (n) is expressed as the following equation (1).
また、同じ時間窓に含まれるM点の呼吸波形データをY(m)とすると、Y(m)は以下の(2)式のように表される。 Further, if the respiratory waveform data at point M included in the same time window is Y (m), Y (m) is expressed as the following equation (2).
ただし、n=1,2,…N,m=1,2,…,Mとする。この時、X(n)をt2 mに合わせてスプライン補間したものをX'(m)とすると、X'(m)は、以下の(3)式のように表される。 Here, n = 1, 2,... N, m = 1, 2,. At this time, assuming that X ′ (m) is the result of spline interpolation with X (n) adjusted to t 2 m , X ′ (m) is expressed by the following equation (3).
以上より、求められたX'(m)とY(m)の相関係数γは、以下の(4)式で求められる。 From the above, the obtained correlation coefficient γ between X ′ (m) and Y (m) is obtained by the following equation (4).
ただし、 ̄x'は、心拍間隔の時系列データのスプライン補間X'(m)の平均値、 ̄yは,呼吸波形データY(m)の平均値を示す。上記(4)式によれば、心拍間隔の時系列データと呼吸波形データの類似度が高いほど、相関係数は大きくなる。 Here,  ̄x ′ represents an average value of spline interpolation X ′ (m) of time series data of heartbeat intervals, and  ̄y represents an average value of respiratory waveform data Y (m). According to the above equation (4), the correlation coefficient increases as the similarity between the time series data of the heartbeat interval and the respiratory waveform data increases.
また、パワースペクトル密度算出部34は、予め定められた時間長分の、直流成分を除去した呼吸関連生体信号の時系列データに対して離散的フーリエ変換を行い、パワースペクトル密度を算出する。
The power spectrum
ここで、呼吸または呼吸間隔のばらつきとパワースペクトル密度との関係の模式図を、図3に示す。睡眠深度が深い状態での呼吸間隔は、個体差や時間帯によって差があるものの、ある程度一定の範囲内に収まると考えられる。その範囲の最長間隔から求められた呼吸周波数をf1、最短間隔から求められた周波数をf2とすると、周波数f1からf2までのパワースペクトル密度の和が求められる。この和が全帯域のパワースペクトル密度の和に占める割合は、呼吸周波数がf1からf2までの間で安定しているかどうかの指標となり、呼吸の安定性を示すパワースペクトル密度の特徴量として使用できる。 Here, FIG. 3 shows a schematic diagram of the relationship between respiration or breathing interval variation and power spectral density. The breathing interval when the sleep depth is deep is considered to be within a certain range to some extent, although there are differences depending on individual differences and time zones. Assuming that the respiration frequency obtained from the longest interval in the range is f 1 and the frequency obtained from the shortest interval is f 2 , the sum of power spectral densities from frequencies f 1 to f 2 is obtained. Percentage this sum to the sum of the power spectral density of the entire band becomes the indication of whether the breathing frequency is stable between from f 1 to f 2, as the feature amount of the power spectral density showing the stability of the respiratory Can be used.
そこで、パワースペクトル密度特徴量算出部36によって、算出されたパワースペクトル密度に基づいて、パワースペクトル密度特徴量として、予め定められた周波数帯域(f1からf2まで)におけるパワースペクトル密度の和が、パワースペクトル密度全体に占める割合を算出する。なお、パワースペクトル密度特徴量が、パワースクペクトル密度の周波数依存性の一例である。
Therefore, the power spectral density characteristic
意識状態推定部38は、算出された相関係数及びパワースペクトル密度特徴量に基づいて、以下に説明するように、推定対象者の意識低下状態を推定する。
Based on the calculated correlation coefficient and power spectral density feature quantity, the consciousness
まず、相関係数と意識低下状態との関連性について述べる。心拍と呼吸との関連については、吸気時に心拍数が上昇(心拍間隔が短縮)し、呼気時に心拍数が下降(心拍間隔が延長)することがわかっている(上記図2参照)。ここで、人間が覚醒状態にある時は、外部からの刺激、体動や発話等の能動的な活動の影響のため、心拍に対する副交感神経系を介した呼吸の寄与は低下すると考えられる。これに対して、人間の意識が低下している状態にある時(睡眠等)は、外部からの刺激は遮断され、体動や発話等の能動的な活動も消失するため、心拍に対する副交感神経系を介した呼吸の寄与は相対的に上昇すると考えられる。以上より、リサンプリングされた心拍間隔の時系列データに対する、直流成分を除去した呼吸関連生体信号の時系列データの、相関係数の絶対値が大きいほど、推定対象者となる人間の意識レベルは低下(睡眠に近い状態)し、相関係数の絶対値が小さいほど推定対象者となる人間の意識レベルは上昇(覚醒に近い状態)していると考えられる。 First, the relationship between the correlation coefficient and the state of reduced consciousness is described. Regarding the relationship between heart rate and respiration, it is known that the heart rate increases (increases the heart rate interval) during inspiration, and decreases (expands the heart rate interval) during expiration (see FIG. 2 above). Here, when a human is in an arousal state, it is considered that the contribution of breathing via the parasympathetic nervous system to the heartbeat decreases due to the influence of active activities such as external stimulation, body movement and speech. On the other hand, when human consciousness is reduced (sleeping, etc.), external stimuli are blocked and active activities such as body movements and speech are lost. The contribution of breathing through the system is thought to increase relatively. From the above, the higher the absolute value of the correlation coefficient of the time-series data of the respiratory related biological signal from which the direct current component is removed with respect to the time-series data of the resampled heartbeat interval, the higher the consciousness level of the human being who is the estimation target It is considered that the consciousness level of the human being who is the estimation target increases (the state is close to awakening) as the absolute value of the correlation coefficient decreases (the state is close to sleep).
次に、呼吸のパワースペクトル密度特徴量と人間の意識状態との関係について述べる。人間が覚醒状態で安静にしている時、呼吸はある程度一定となるが、そこは意識の作用によりある程度の揺らぎが生じている。また、覚醒状態であっても活動状態(動作や発話等をしている状態)では、呼吸のばらつきは大きくなると考えられる。一方、覚醒状態から睡眠状態に至る過渡期においては、呼吸間隔が顕著にばらつくことが知られており、安静覚醒時に比べて呼吸がばらつく方向に変化すると考えられる。さらに、睡眠状態に至ると、意識による作用がなくなるため、呼吸間隔は一定となり、呼吸が安定化する方向に変化する。なお、パワースペクトル密度特徴量は、呼吸が安定するほど増加する。 Next, the relationship between the respiratory power spectral density feature and the human consciousness state will be described. When a person is awake and rests, the breathing is constant to some extent, but there is some fluctuation due to the action of consciousness. In addition, even in an awake state, in an active state (a state in which an operation or speech is performed), it is considered that the variation in breathing becomes large. On the other hand, in the transition period from the awake state to the sleep state, it is known that the respiration interval varies remarkably, and it is considered that the respiration varies in a direction in which the respiration varies compared to the rest awake state. Furthermore, when the sleep state is reached, the action due to consciousness is lost, so the breathing interval becomes constant, and the breathing changes in a stabilizing direction. Note that the power spectral density feature amount increases as the respiration stabilizes.
そこで、意識状態推定部38は、リサンプリングされた心拍間隔の時系列データに対する直流成分を除去した呼吸関連生体信号の時系列データの相関係数の絶対値と、パワースペクトル密度特徴量とに対し、任意の重み付けをした後の積を、意識レベルを推定する指標として算出する。本指標は意識レベルが低下する(睡眠方向に向う)ほど増加し、上昇する(覚醒方向に向う)ほど減少する。
Therefore, the consciousness
意識状態推定部38は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18により推定結果を出力させる。
The consciousness
なお、上記の閾値として、従来の医学的治験により、汎用性のある値をあらかじめ与えることができる。また、個人差に対応した閾値を設定するため、次に述べるような方法で定めることもできる。まず、対象となる人の就寝中の心拍間隔の時系列データ及び呼吸関連生体信号の時系列データを記録する。この時系列データには、安静覚醒状態、覚醒から睡眠への移行状態、睡眠状態の3状態の時系列データが少なくとも含まれている。そこで、就寝時のリサンプリングされた心拍間隔の時系列データに対する直流成分を除去した呼吸関連生体信号の時系列データの相関係数の絶対値、及びパワースペクトル密度特徴量を算出し、それを睡眠状態と睡眠状態以外の状態でクラスタリング等をすることにより、個人差に対応した閾値を定めることが可能である。 As the above threshold value, a versatile value can be given in advance by a conventional medical trial. Moreover, in order to set the threshold value corresponding to individual differences, it can be determined by the following method. First, time-series data of heartbeat intervals during sleep of a subject person and time-series data of respiratory related biological signals are recorded. This time series data includes at least time series data of three states of a resting awake state, a transition state from awakening to sleep, and a sleep state. Therefore, the absolute value of the correlation coefficient of the time series data of the respiratory related biological signal from which the DC component is removed from the time series data of the resampled heartbeat interval at bedtime and the power spectrum density feature amount are calculated, By performing clustering or the like in a state other than the state and the sleep state, it is possible to determine a threshold corresponding to individual differences.
次に、第1の実施の形態に係る意識状態推定装置10の作用について説明する。
Next, the operation of the consciousness
まず、心拍検出装置12によって推定対象者の心拍間隔を検出し、リサンプリングされた心拍間隔の時系列データをコンピュータ16に入力する。また、呼吸検出装置14によって、推定対象者の呼吸関連生体信号を検出し、直流成分を除去した呼吸関連生体信号の時系列データをコンピュータ16に入力する。また、コンピュータ16において、図4に示す意識状態推定処理ルーチンが繰り返し実行される。
First, the
まず、ステップ100で、直近の所定時間分の、リサンプリングされた心拍間隔の時系列データ及び直流成分を除去した呼吸関連生体信号の時系列データを取得する。
First, in
そして、ステップ102において、上記ステップ100で取得した心拍間隔の時系列データ及び呼吸関連生体信号の時系列データについて、相関係数を算出する。
In
次のステップ104では、上記ステップ100で取得した呼吸関連生体信号の時系列データについて、パワースペクトル密度を算出し、ステップ106において、パワースペクトル密度特徴量として、予め定められた周波数帯域(f1〜f2)におけるパワースペクトル密度の和の割合を算出する。
In the
そして、ステップ108において、上記ステップ102で算出した相関係数、及び上記ステップ106で算出したパワースペクトル密度特徴量に基づいて、意識レベルを推定するための指標を算出する。次のステップ110では、上記ステップ108で算出された指標が、閾値以上であるか否かに基づいて、推定対象者が意識低下状態であるか否かを推定し、ステップ112において、出力部18によって推定結果を出力させて、意識状態推定処理ルーチンを終了する。
In
以上説明したように、第1の実施の形態に係る意識状態推定装置によれば、検出された心拍間隔の時系列データ及び呼吸関連生体信号の時系列データの間の相関係数を算出し、相関係数が高いほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を精度良く推定することができる。また、呼吸のばらつき度合いを示すパワースペクトル密度特徴量を算出して、パワースペクトル特徴量が大きいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態をより精度良く推定することができる。 As described above, according to the consciousness state estimation device according to the first embodiment, the correlation coefficient between the time-series data of the detected heartbeat interval and the time-series data of the respiratory related biological signal is calculated, By estimating the consciousness reduction state of the estimation target person so that the consciousness state decreases as the correlation coefficient increases, the estimation target person's consciousness reduction state can be accurately estimated. In addition, by calculating the power spectrum density feature amount indicating the degree of variation in breathing, and estimating the consciousness state of the estimation subject so that the consciousness state decreases as the power spectrum feature amount increases, The state of reduced consciousness can be estimated with higher accuracy.
本実施の形態では、呼吸関連生体信号の時系列データのパワースペクトル密度の周波数依存性を求め、呼吸のばらつき指標として計測している。呼吸が安定している時(≒睡眠時)は、パワースペクトル密度の周波数依存性のピークが呼吸の特徴帯域周辺に出現するが、呼吸にばらつきが多い時(≒覚醒時や居眠り時)は、より低周波帯域に周波数依存性のピークが移行し、またはピークが存在しなくなると考えられる。これにより、呼吸のピーク等を検出する必要がなくなり、検出のロバスト性が向上するため、推定対象者の意識低下状態を安定して精度良く推定することができる。また、電気インピーダンス法や電波式呼吸センサ、圧力式呼吸センサ等、いかなるセンサで呼吸を検出した場合であっても、本実施の形態の技術を容易に適用することができる。 In the present embodiment, the frequency dependence of the power spectral density of the time-series data of the respiratory related biological signal is obtained and measured as a breathing variation index. When breathing is stable (≒ sleeping), a frequency-dependent peak of power spectral density appears around the feature band of breathing, but when there is a lot of variation in breathing (≈ when awake or dozing) It is considered that a frequency-dependent peak shifts to a lower frequency band or no peak exists. Thereby, it is not necessary to detect the peak of respiration and the like, and the robustness of the detection is improved. Therefore, it is possible to stably and accurately estimate the state of reduced consciousness of the estimation target person. Further, the technique of the present embodiment can be easily applied even when respiration is detected by any sensor such as an electric impedance method, a radio wave respiration sensor, or a pressure respiration sensor.
従来、睡眠による人の意識状態の変化(あるいは睡眠段階等の睡眠の程度)の評価手法としては、心拍間隔の変動の低周波数成分(0.05〜0.15Hz程度)及び高周波数成分(0.15〜0.4Hz程度)のパワースペクトル密度により、自律神経活動を評価する手法が広く用いられていた。しかし、様々な要因で変動するため、評価の精度が低く、実用上の問題となっている。本実施の形態では、意識によって制御可能な呼吸、及び呼吸の影響を受けるが意識による制御は不可能な心拍と呼吸の比較により、より精度の高い意識状態推定装置を実現する。 Conventionally, as a method for evaluating a change in a person's consciousness state due to sleep (or the degree of sleep such as a sleep stage), a low frequency component (about 0.05 to 0.15 Hz) and a high frequency component (0 A method for evaluating autonomic nerve activity based on a power spectral density of about .15 to 0.4 Hz has been widely used. However, since it fluctuates due to various factors, the accuracy of evaluation is low, which is a practical problem. In the present embodiment, a more accurate consciousness state estimation device is realized by comparing respiration that can be controlled by consciousness and heartbeat and respiration that are affected by respiration but cannot be controlled by consciousness.
心拍間隔の変動のうちの高周波数成分を用いた自律神経活動評価は、副交感神経系を介した呼吸による心拍制御を、心拍間隔を用いて間接的に評価しているものであり、間接的であるため、他の要因の影響を受けやすいのは自明である。本実施の形態では、呼吸及び呼吸と心拍の関係を直接的に評価するものであり、従来手法に比べて評価精度を向上させることができる。 Autonomic nerve activity evaluation using high-frequency components of heart rate interval fluctuations is an indirect evaluation of heart rate control by breathing via the parasympathetic nervous system using heart rate intervals. It is obvious that it is susceptible to other factors. In this embodiment, the relationship between respiration and respiration and heartbeat is directly evaluated, and the evaluation accuracy can be improved as compared with the conventional method.
本実施の形態では、時間軸上の相関性を見るため、必要とする時間窓が短く、意識低下状態を推定するときの時間分析能が短くなる。 In the present embodiment, since the correlation on the time axis is seen, the necessary time window is short, and the time analysis ability when estimating the state of consciousness is shortened.
次に、第2の実施の形態について説明する。なお、第1の実施の形態と同様となる部分については、同一符号を付して説明を省略する。 Next, a second embodiment will be described. In addition, about the part similar to 1st Embodiment, the same code | symbol is attached | subjected and description is abbreviate | omitted.
第2の実施の形態では、心拍センサによって、心電図波形を示す心拍関連生体信号を検出している点と、心拍関連生体信号の時系列データ及び呼吸関連生体信号の時系列データの間の相互情報量を算出している点とが、第1の実施の形態と主に異なっている。 In the second embodiment, the mutual information between the point where the heartbeat sensor detects the heartbeat-related biological signal indicating the electrocardiogram waveform and the time-series data of the heartbeat-related biological signal and the time-series data of the respiratory-related biological signal The point of calculating the amount is mainly different from the first embodiment.
図5に示すように、第2の実施の形態に係る意識状態推定装置210の心拍検出装置212は、心拍センサ220、低域通過フィルタ222、及び心拍波形リサンプリング部224を備えている。心拍センサ220は、人体より心臓の拍動に関連する生体信号(心電図波形)を検出するセンサであり、心電図波形を示す心拍関連生体信号の時系列データを出力する。なお、心拍検出装置212が、信号検出手段の一例であり、心電図波形を示す心拍関連生体信号が、心臓の拍動に応じて変化する信号の一例である。呼吸検出装置14が、データ検出手段の一例である。
As shown in FIG. 5, the
低域通過フィルタ222のカットオフ周波数は,平均的な心拍周波数及び呼吸周波数を考慮して決定されるが、心拍の周波数より高い周波数で設定することが望ましい。低域通過フィルタ222によって、呼吸に由来する成分が含まれる周波数帯域の信号が分離される。
The cut-off frequency of the low-
心拍波形リサンプリング部224は、高域遮断された心拍関連生体信号の時系列データを、呼吸センサ26のサンプリング時刻に合わせてリサンプリングして出力する。
The heartbeat
コンピュータ216は、心拍検出装置212から出力される心拍関連生体信号の時系列データ及び呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、心拍関連生体信号の時系列データと呼吸関連生体信号の時系列データとの間の相互情報量を算出する相互情報量算出部232と、パワースペクトル密度算出部34と、パワースペクトル密度特徴量算出部36と、算出された相互情報量及びパワースペクトル密度特徴量に基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部38とを備えている。なお、相互情報量算出部232が、相関性判定手段の一例であり、相互情報量が、相関性の一例である。
The
相互情報量算出部232は、リサンプリングされた心拍関連生体信号の時系列データ及び直流成分を除去した呼吸関連生体信号の時系列データのうち、予め定められた時間長における、心拍関連生体信号の時系列データ及び呼吸関連生体信号の時系列データ間の相互情報量を、以下のように算出する。
The mutual information
まず、ある時間窓に含まれるN点の心拍関連生体信号の時系列データをX(n)とすると、X(n)は、以下の(5)式のように表される。 First, assuming that time series data of N points of heartbeat-related biological signals included in a certain time window is X (n), X (n) is expressed by the following equation (5).
また、同じ時間窓に含まれるM点の呼吸関連生体信号の時系列データをY(m)とすると、Y(m)は以下の(6)式のように表される。 If time series data of M-point respiratory-related biological signals included in the same time window is Y (m), Y (m) is expressed as the following equation (6).
ただし、n=1,2,…N、m=1,2,…,Mとする。このとき,X(n)をt2 mに合わせてスプライン補間したものをX'(m)とすると、X'(m)は以下の(7)式のように表される。 Here, n = 1, 2,... N, m = 1, 2,. At this time, if X ′ (m) is obtained by performing spline interpolation with X (n) adjusted to t 2 m , X ′ (m) is expressed by the following equation (7).
以上より求められたX'(m)とY(m)について、以下の処理を行う。 The following processing is performed for X ′ (m) and Y (m) obtained from the above.
まず、X'(m)とY(m)をそれぞれ平均=0、標準偏差=1となるよう正規化した時系列データ ̄X'(m)と ̄Y(m)を求める。次に、 ̄X'(m)をP個の区間に分けた時にp番目の区間Xbin(p)に属する ̄X'(m)の点数をLpとする。 ̄Y(m)をQ個の区間に分けた時にq番目の区間Ybin(q)に属する ̄Y(m)の点数をLqとする。 ̄X'(m)がXbin(p)に属し、かつ ̄Y(m)がYbin(q)に属する点数をLp,qとする。このとき、X'とYの相互情報量I(X,Y')は、以下の(8)式のように求められる。 First, time series data  ̄X ′ (m) and  ̄Y (m) are obtained by normalizing X ′ (m) and Y (m) so that the average = 0 and the standard deviation = 1, respectively. Next, when  ̄X ′ (m) is divided into P sections, the score of  ̄X ′ (m) belonging to the p-th section Xbin (p) is L p . When  ̄Y (m) is divided into Q sections, the score of  ̄Y (m) belonging to the qth section Ybin (q) is L q . Let L p, q be the number of points where  ̄X ′ (m) belongs to Xbin (p) and  ̄Y (m) belongs to Ybin (q). At this time, the mutual information I (X, Y ′) between X ′ and Y is obtained as in the following equation (8).
上記(8)式によれば、心拍関連生体信号の時系列データをスプライン補完したものと呼吸関連生体信号の時系列データとの相関性(線形・非線形含む)が完全に独立な場合、相互情報量は0となり、相関性が高くなるほど、相互情報量は増加する。 According to the above equation (8), if the correlation (including linear and non-linear) between the time-series data of the heartbeat-related biological signal and the time-series data of the respiratory-related biological signal is completely independent, the mutual information The amount is 0, and the mutual information amount increases as the correlation increases.
意識状態推定部38は、リサンプリングされた心拍関連生体信号の時系列データに対する、直流成分を除去した呼吸関連生体信号の時系列データの相互情報量と、パワースペクトル密度特徴量とに対し、任意の重み付けをした後の積を、意識レベルを推定する指標として算出する。
The
意識状態推定部38は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。
The consciousness
なお、第2の実施の形態に係る意識状態推定装置の他の構成は、第1の実施の形態と同様であるため、説明を省略する。 In addition, since the other structure of the consciousness state estimation apparatus which concerns on 2nd Embodiment is the same as that of 1st Embodiment, description is abbreviate | omitted.
以上説明したように、第2の実施の形態に係る意識状態推定装置によれば、推定対象者の心臓の拍動に応じて変化する心拍関連生体信号から分離された、呼吸に由来する成分が含まれる周波数帯域の信号、及び呼吸関連生体信号の時系列データの相互情報量を算出し、相互情報量が大きいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を安定して精度良く推定することができる。また、呼吸のばらつき度合いを示すパワースペクトル密度特徴量を算出して、パワースペクトル特徴量が大きいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態をより精度良く推定することができる。 As described above, according to the consciousness state estimation device according to the second embodiment, the component derived from respiration separated from the heartbeat-related biological signal that changes in accordance with the heartbeat of the estimation target person is included. By calculating the mutual information amount of the time-series data of the signal in the frequency band included and the respiratory related biological signal, and estimating the consciousness lowering state of the estimation target person so that the consciousness state decreases as the mutual information amount increases. Thus, it is possible to stably and accurately estimate the state of consciousness reduction of the estimation target person. In addition, by calculating the power spectrum density feature amount indicating the degree of variation in breathing, and estimating the consciousness state of the estimation subject so that the consciousness state decreases as the power spectrum feature amount increases, The state of reduced consciousness can be estimated with higher accuracy.
また、呼吸と心拍の相関性に注目しており、その指標の一つである相互情報量は心拍及び呼吸間の非線形の相関性も反映できる指標である。本指標を用いることにより、心拍及び呼吸がいかなる検出手段で検出されたものであっても、その相関性を数値化することが可能になる。 Further, attention is focused on the correlation between respiration and heartbeat, and the mutual information amount, which is one of the indices, is an index that can also reflect the non-linear correlation between heartbeat and respiration. By using this index, it is possible to quantify the correlation between heart rate and respiration detected by any detection means.
なお、上記の第1の実施の形態及び第2の実施の形態では、パワースペクトル密度特徴量として、睡眠状態でパワースペクトル密度が大きくなる周波数帯域のパワースペクトル密度の和が、パワースペクトル密度全体に占める割合を算出する場合を例に説明したが、これに限定されるものではない。 In the first embodiment and the second embodiment described above, the sum of the power spectrum density in the frequency band in which the power spectrum density increases in the sleep state is the power spectrum density feature amount in the entire power spectrum density. Although the case of calculating the occupying ratio has been described as an example, the present invention is not limited to this.
例えば、パワースペクトル密度特徴量として、パワースペクトル密度のピーク周波数を求めてもよい。上記図3に示すように、呼吸にばらつきがある状態では呼吸のパワースペクトル密度がピークとなる周波数は直流成分に極めて近いが、呼吸が安定している状態では呼吸ピークはより高周波側に存在することがわかる。以上より、呼吸パワースペクトル密度のピーク周波数を、呼吸の安定性を示すパワースペクトル密度特徴量として用いることができる。 For example, the peak frequency of the power spectral density may be obtained as the power spectral density feature quantity. As shown in FIG. 3, the frequency at which the breathing power spectral density reaches a peak is extremely close to a direct current component when there is variation in breathing, but the breathing peak is on the higher frequency side when breathing is stable. I understand that. From the above, the peak frequency of the respiratory power spectral density can be used as a power spectral density feature amount indicating the stability of respiration.
また、上記のピーク周波数が周波数f1及びf2の間にあるとき、上記のピーク周波数におけるパワースペクトル密度が、パワースペクトル密度全体に占める割合を、パワースペクトル特徴量として算出するようにしてもよい。呼吸にばらつきがある状態では、全帯域のパワースペクトル密度に占める、ピーク周波数におけるパワースペクトル密度の割合が、呼吸が安定している状態に比べて低下すると考えられる。以上より、上記のピーク周波数におけるパワースペクトル密度が全帯域のパワースペクトル密度に占める割合を、呼吸の安定を示すパワースペクトル密度特徴量として用いることができる。 Further, when the peak frequency is between the frequencies f1 and f2, the ratio of the power spectrum density at the peak frequency to the entire power spectrum density may be calculated as the power spectrum feature amount. In a state where there is variation in respiration, it is considered that the ratio of the power spectral density at the peak frequency to the power spectral density in the entire band is lower than in a state where the respiration is stable. From the above, the ratio of the power spectrum density at the peak frequency to the power spectrum density of the entire band can be used as the power spectrum density feature amount indicating the stability of respiration.
また、上記の実施の形態では、心拍と呼吸との相関性として、自己相関関数を算出する場合を例に説明したが、これに限定されるものではなく、上記の第1の実施の形態と同様に、心拍と呼吸との相関性として、相互情報量を算出しても良い。 In the above embodiment, the case where the autocorrelation function is calculated as the correlation between the heartbeat and the respiration is described as an example. However, the present invention is not limited to this, and the first embodiment and the above embodiment are not limited thereto. Similarly, mutual information may be calculated as the correlation between heartbeat and respiration.
次に、第3の実施の形態について説明する。なお、第1の実施の形態と同様の構成となる部分については、同一符号を付して説明を省略する。 Next, a third embodiment will be described. In addition, about the part which becomes the structure similar to 1st Embodiment, the same code | symbol is attached | subjected and description is abbreviate | omitted.
第3の実施の形態では、呼吸関連生体信号の時系列データから、呼吸のばらつき度を算出している点が、第1の実施の形態と異なっている。 The third embodiment is different from the first embodiment in that the degree of variation in breathing is calculated from time-series data of breathing related biological signals.
図6に示すように、第3の実施の形態に係る意識状態推定装置310のコンピュータ316は、データ取得部30と、相関係数算出部32と、呼吸関連生体信号の時系列データに基づいて、呼吸のばらつき度を算出する呼吸ばらつき度算出部334と、算出された相関係数及び呼吸ばらつき度に基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部338とを備えている。
As shown in FIG. 6, the
呼吸ばらつき度算出部334は、予め定められた時間長分の、直流成分を除去した呼吸関連生体信号の時系列データに基づいて、以下に説明するように、呼吸関連生体信号の時系列データの自己相関関数を算出し、呼吸ばらつき度とする。
Based on the time-series data of the respiratory-related biological signal from which the DC component has been removed for a predetermined time length, the respiratory variation
まず、ある時間窓におけるM点の呼吸関連生体信号の時系列データY(m)を以下の(9)式のように定義する。 First, time series data Y (m) of M-point respiratory-related biological signals in a certain time window is defined as the following equation (9).
ただし、m=1,2,…,Mとする。この時、自己相関関数R(i)は、以下の(10)式に従って算出される。 Here, m = 1, 2,..., M. At this time, the autocorrelation function R (i) is calculated according to the following equation (10).
ただし、−M<i<Mとする。上記(10)式によれば、呼吸が均一である時は、自己相関関数R(i)のi-=0以外の最大のピーク値は大きくなり、呼吸が不均一な時は、自己相関関数R(i)のi-=0以外の最大のピーク値は小さくなる。また、自己相関関数R(i)のi-=0以外の最大のピーク値が大きいほど、呼吸ばらつき度が小さくなる。 However, −M <i <M. According to the above equation (10), when respiration is uniform, the maximum peak value other than i− = 0 of autocorrelation function R (i) is large, and when respiration is nonuniform, autocorrelation function The maximum peak value of R (i) other than i− = 0 is small. In addition, the greater the maximum peak value other than i− = 0 of the autocorrelation function R (i), the smaller the degree of respiratory variability.
意識状態推定部338は、リサンプリングされた心拍間隔の時系列データに対する、直流成分を除去した呼吸関連生体信号の時系列データの相関係数の絶対値と、算出された自己相関関数とに対し、任意の重み付けをした後の積を、意識レベルを推定する指標として算出する。本指標は意識レベルが低下する(睡眠方向に向う)ほど増加し、上昇する(覚醒方向に向う)ほど減少する。
The consciousness
意識状態推定部338は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18によって、推定結果を出力させる。
The consciousness
なお、第3の実施の形態に係る意識状態推定の他の構成及び作用については、第1の実施の形態と同様であるため、説明を省略する。 In addition, about the other structure and effect | action of consciousness state estimation which concern on 3rd Embodiment, since it is the same as that of 1st Embodiment, description is abbreviate | omitted.
以上説明したように、第3の実施の形態に係る意識状態推定装置によれば、呼吸のばらつき度合いを示す指標として、呼吸関連生体信号の時系列データの自己相関関数R(i)のi-=0以外の最大のピーク値を算出して、自己相関関数R(i)のi-=0以外の最大のピーク値が大きいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を精度良く推定することができる。 As described above, according to the consciousness state estimation apparatus according to the third embodiment, i − of the autocorrelation function R (i) of the time-series data of the respiratory related biological signal is used as an index indicating the degree of variation in respiration. The maximum peak value other than = 0 is calculated, and the consciousness state of the estimation target person is reduced so that the consciousness state decreases as the maximum peak value other than i− = 0 of the autocorrelation function R (i) increases. By estimating, the estimation target person's consciousness fall state can be estimated accurately.
次に、第4の実施の形態について説明する。なお、第4の実施の形態に係る意識状態推定装置は、第3の実施の形態と同様の構成となっているため、同一符号を付して説明を省略する。 Next, a fourth embodiment will be described. In addition, since the consciousness state estimation apparatus according to the fourth embodiment has the same configuration as that of the third embodiment, the same reference numerals are given and description thereof is omitted.
第4の実施の形態では、呼吸のばらつき度合いとして、変動係数を算出している点が、第3の実施の形態と主に異なっている。 The fourth embodiment is mainly different from the third embodiment in that a variation coefficient is calculated as the degree of breathing variation.
第4の実施の形態に係る意識状態推定装置310では、呼吸ばらつき度算出部334によって、予め定められた時間長分の、直流成分を除去した呼吸関連生体信号の時系列データに基づいて、以下に説明するように、呼吸間隔の変動係数を、呼吸ばらつき度として算出する。
In the consciousness
図7に示すような呼吸波形(呼吸関連生体信号の時系列データ)を考えると、ある時間窓の中にn+1点の呼吸ピークがあるとき、n個の呼吸点間の間隔(呼吸間隔)が得られる。ここで、n個の呼吸間隔の標準偏差をSDnとし、平均をMEANnとすると、変動係数CVnは、以下の(11)式に従って算出される。 Considering the respiration waveform (time-series data of respiration-related biological signals) as shown in FIG. 7, when there are n + 1 respiration peaks in a certain time window, the interval between the n respiration points (respiration interval) is can get. Here, assuming that the standard deviation of n breath intervals is SD n and the average is MEAN n , the variation coefficient CV n is calculated according to the following equation (11).
上記(11)式によれば、呼吸間隔がばらつくとSDnが増加するため、変動係数CVnは増加する。すなわち、変動係数CVnが大きいほど、呼吸ばらつき度が大きくなる。 According to the above equation (11), when the breathing interval varies, SD n increases, so the variation coefficient CV n increases. That is, the larger the coefficient of variation CV n, respiratory variation degree becomes large.
なお、第4の実施の形態に係る意識状態推定の他の構成及び作用については、第1の実施の形態及び第3の実施の形態と同様であるため、説明を省略する。 In addition, about the other structure and effect | action of consciousness state estimation which concern on 4th Embodiment, since it is the same as that of 1st Embodiment and 3rd Embodiment, description is abbreviate | omitted.
以上説明したように、第4の実施の形態に係る意識状態推定装置によれば、呼吸のばらつき度合いを示す指標として、変動係数を算出して、変動係数の値が小さいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を精度良く推定することができる。 As described above, according to the consciousness state estimation device according to the fourth embodiment, the variation coefficient is calculated as an index indicating the degree of variation in breathing, and the consciousness state decreases as the variation coefficient value decreases. As described above, by estimating the consciousness reduction state of the estimation target person, it is possible to accurately estimate the consciousness reduction state of the estimation target person.
なお、呼吸のばらつき度合いを示す指標として、変動係数を算出する場合を例に説明したが、これに限定されるものではなく、呼吸のばらつき度合いを示す指標として、呼吸関連生体信号の時系列データの標準偏差を算出してもよい。この場合には、標準偏差の値が小さいほど意識状態が低下するように、推定対象者の意識低下状態を推定すればよい。 Note that the case where the coefficient of variation is calculated as an index indicating the degree of breathing variation has been described as an example. However, the present invention is not limited to this, and time series data of respiratory related biological signals is used as an index indicating the degree of breathing variation. The standard deviation may be calculated. In this case, the state of consciousness reduction of the estimation target person may be estimated so that the state of consciousness decreases as the standard deviation value decreases.
次に、第5の実施の形態について説明する。なお、第5の実施の形態に係る意識状態推定装置の構成は、第1の実施の形態と同様の構成となるため、同一符号を付して説明を省略する。 Next, a fifth embodiment will be described. In addition, since the structure of the consciousness state estimation apparatus which concerns on 5th Embodiment becomes a structure similar to 1st Embodiment, it attaches | subjects the same code | symbol and abbreviate | omits description.
第5の実施の形態では、意識低下状態の推定方法が、第1の実施の形態と異なっている。 In the fifth embodiment, the method for estimating the state of consciousness is different from that in the first embodiment.
第5の実施の形態に係る意識状態推定装置10では、意識状態推定部338によって、相関係数の絶対値と、パワースペクトル密度特徴量との2軸で張られる座標上の位置に基づいて、推定対象者が意識低下状態であるか否かを推定する。この時に張られる座標の模式図を図8に示す。なお、意識低下状態(睡眠状態)とそれ以外の状態との境界は、上記の第1の実施の形態で説明した方法と同様の方法により、定めることができる。
In the consciousness
なお、第5の実施の形態に係る意識状態推定の他の構成及び作用については、第1の実施の形態と同様であるため、説明を省略する。 In addition, about the other structure and effect | action of consciousness state estimation which concern on 5th Embodiment, since it is the same as that of 1st Embodiment, description is abbreviate | omitted.
次に、第6の実施の形態について説明する。なお、第1の実施の形態と同様の構成となる部分については、同一符号を付して説明を省略する。 Next, a sixth embodiment will be described. In addition, about the part which becomes the structure similar to 1st Embodiment, the same code | symbol is attached | subjected and description is abbreviate | omitted.
第6の実施の形態では、推定対象者の体動及び発話を検出し、検出結果を更に考慮して、意識状態を推定している点が、第1の実施の形態と主に異なっている。 The sixth embodiment is mainly different from the first embodiment in that the body movement and utterance of the estimation target person are detected, and the consciousness state is estimated further considering the detection result. .
図9に示すように、第6の実施の形態に係る意識状態推定装置610は、心拍検出装置12と、呼吸検出装置14と、推定対象者の体動を検出する体動検出装置612と、推定対象者の発話を検出する発話検出装置614と、コンピュータ16とを備えている。
As shown in FIG. 9, the consciousness
体動検出装置612としては、例えば、人間に装着するセンサとしての加速度センサを用いることができ、あるいは、人間の周辺環境に装着するセンサとしての、人間と接触する部位の圧力を計測するセンサ、カメラ、電波式距離センサ、電波式速度センサ、音波式距離センサ、音波式速度センサ等を用いることができる。また、心拍センサ20及び呼吸センサ26とセンサを共有することも可能である。
As the body
また、発話検出装置614としては、マイクロフォンなどを用いることができる。
As the
コンピュータ16のデータ取得部30は、心拍検出装置12から入力された心拍間隔の時系列データを取得し、メモリ(図示省略)に記憶する共に、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得してメモリに記憶する。また、データ取得部30は、体動検出装置612から入力された体動検出結果の時系列データを取得し、メモリ(図示省略)に記憶する共に、発話検出装置614から出力される発話検出結果の時系列データを取得してメモリに記憶する。
The
ここで、本実施の形態の原理について説明する。 Here, the principle of the present embodiment will be described.
上記の実施の形態では、睡眠状態と睡眠以外の状態のみを判別する場合を例に示したが、睡眠以外の状態には覚醒状態と覚醒から睡眠に移行する過渡状態が含まれており、これを判別することができればさらに有用である。しかし、呼吸のばらつきには、体動や発話等の人間の能動的活動に起因するものと、覚醒から睡眠に移行する過渡状態において意図せず生じるものの2種類があり、睡眠以外の状態を覚醒状態と覚醒から睡眠に移行する過渡状態に分けるためには、呼吸のばらつきが能動的活動によるものかどうかを判別する必要がある。 In the above embodiment, the case where only the sleep state and the state other than the sleep are discriminated is shown as an example. However, the state other than the sleep includes the awake state and the transition state from the awake state to the sleep state. It is more useful if it can be determined. However, there are two types of variability in breathing: those caused by active human activities such as body movements and speech, and those that occur unintentionally during the transition from wakefulness to sleep. In order to divide the state into a transitional state from wakefulness to sleep, it is necessary to determine whether or not the variation in breathing is due to active activity.
そこで、本実施の形態では、意識状態推定部38によって、リサンプリングされた心拍間隔の時系列データに対する、直流成分を除去した呼吸関連生体信号の時系列データの相関係数の絶対値と、パワースペクトル密度特徴量とに基づいて、意識レベルを推定する指標を算出し、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定する。算出した上記の指標が閾値未満である場合には、体動検出結果の時系列データ及び発話検出結果の時系列データに基づいて、体動及び発話の少なくとも一方が一定以上の強度であるか否かを判定する。体動及び発話の少なくとも一方が一定以上の強度である場合(体動及び発話の少なくとも一方が検出された場合)、呼吸のばらつきは覚醒状態の体動または発話に起因するものであると判断し、推定対象者の意識状態が覚醒状態であると推定する。
Therefore, in the present embodiment, the absolute value of the correlation coefficient of the time-series data of the respiratory related biological signal from which the direct current component has been removed with respect to the time-series data of the heartbeat interval resampled by the consciousness
また、意識状態推定部38は、算出した上記の指標が閾値未満である場合であって、かつ、体動検出結果の時系列データ及び発話検出結果の時系列データに基づいて体動及び発話の何れもが検出されない場合、推定対象者の意識状態が、覚醒から睡眠の過渡状態であると推定する。
In addition, the consciousness
なお、第6の実施の形態に係る意識状態推定の他の構成及び作用については、第1の実施の形態と同様であるため、説明を省略する。 In addition, about the other structure and effect | action of consciousness state estimation which concern on 6th Embodiment, since it is the same as that of 1st Embodiment, description is abbreviate | omitted.
このように、体動又は発話の検出結果を考慮することにより、睡眠以外の状態として、覚醒状態と、覚醒から睡眠に移行する過渡状態とを区別して推定することができる。 Thus, by considering the detection result of body movement or utterance, it is possible to distinguish and estimate an arousal state and a transient state from awakening to sleep as a state other than sleep.
次に、第7の実施の形態について説明する。なお、第1の実施の形態及び第3の実施の形態と同様の構成となる部分については、同一符号を付して説明を省略する。 Next, a seventh embodiment will be described. In addition, about the part which becomes the structure similar to 1st Embodiment and 3rd Embodiment, the same code | symbol is attached | subjected and description is abbreviate | omitted.
第7の実施の形態では、呼吸関連生体信号の時系列データに基づいて、推定対象者が意識低下状態であるか否かを推定している点が、第1の実施の形態と異なっている。 The seventh embodiment is different from the first embodiment in that it is estimated whether or not the estimation target person is in a state of reduced consciousness based on time-series data of respiratory related biological signals. .
図10に示すように、第7の実施の形態に係る意識状態推定装置710は、呼吸検出装置14と、呼吸検出装置14から出力される呼吸関連生体信号の時系列データに基づいて、推定対象者が意識低下状態であるか否かを推定し、推定結果を出力部18により出力させるコンピュータ716とを備えている。
As shown in FIG. 10, the consciousness
コンピュータ716は、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、呼吸関連生体信号の時系列データに基づいて、呼吸のばらつき度合いを算出する呼吸ばらつき度算出部334と、算出された呼吸ばらつき度合いに基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部738とを備えている。
The
意識状態推定部738は、呼吸ばらつき度として算出された自己相関関数R(i)のi-=0以外の最大のピーク値を、意識レベルを推定する指標として、算出した上記の指標(自己相関関数R(i)のi-=0以外の最大のピーク値)が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18によって、推定結果を出力させる。
The
次に、第7の実施の形態に係る意識状態推定装置710の作用について説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
Next, the operation of the consciousness
まず、呼吸検出装置14によって、推定対象者の呼吸関連生体信号を検出し、直流成分を除去した呼吸関連生体信号の時系列データをコンピュータ716に入力する。また、コンピュータ716において、図11に示す意識状態推定処理ルーチンが繰り返し実行される。
First, the
まず、ステップ750で、直近の所定時間分の、直流成分を除去した呼吸関連生体信号の時系列データを取得する。
First, in
そして、ステップ752において、上記ステップ750で取得した呼吸関連生体信号の時系列データについて、上記(10)式に従って自己相関関数R(i)のi-=0以外の最大のピーク値を算出して、意識レベルを推定するための指標とする。
Then, in
次のステップ754では、上記ステップ752で算出された自己相関関数R(i)のi-=0以外の最大のピーク値が、閾値以上であるか否かに基づいて、推定対象者が意識低下状態であるか否かを推定し、ステップ112において、出力部18によって推定結果を出力させて、意識状態推定処理ルーチンを終了する。
In the
以上説明したように、第7の実施の形態に係る意識状態推定装置によれば、呼吸のばらつき度合いを示す指標として、呼吸関連生体信号の時系列データの自己相関関数R(i)のi-=0以外の最大のピーク値を算出して、自己相関関数R(i)のi-=0以外の最大のピーク値が大きいほど意識状態が低下するように、推定対象者の意識低下状態を推定することにより、推定対象者の意識低下状態を精度良く推定することができる。 As described above, according to the consciousness state estimation apparatus according to the seventh embodiment, i − of the autocorrelation function R (i) of the time-series data of the respiratory related biological signal is used as an index indicating the degree of variation in respiration. The maximum peak value other than = 0 is calculated, and the consciousness state of the estimation target person is reduced so that the consciousness state decreases as the maximum peak value other than i− = 0 of the autocorrelation function R (i) increases. By estimating, the estimation target person's consciousness fall state can be estimated accurately.
なお、上記の第7の実施の形態において、呼吸のばらつき度合いを示す指標として、呼吸関連生体信号の時系列データの自己相関関数R(i)のi-=0以外の最大のピーク値を算出する場合を例に説明したが、上記の第1の実施の形態と同様に、パワースペクトル密度の周波数依存性を示す特徴量を算出するようにしてもよい。 In the seventh embodiment, the maximum peak value other than i− = 0 of the autocorrelation function R (i) of the time-series data of the respiratory related biological signal is calculated as an index indicating the degree of variation in respiration. Although the case of doing this has been described as an example, the feature quantity indicating the frequency dependence of the power spectrum density may be calculated as in the first embodiment.
また、上記の第6の実施の形態と同様に、推定対象者の体動又は発話を検出し、検出結果を更に考慮して、推定対象者の意識低下状態を推定するようにしてもよい。 Further, similarly to the sixth embodiment described above, body movement or speech of the estimation target person may be detected, and the detection target's consciousness lowered state may be estimated in consideration of the detection result.
また、上記の第1の実施の形態〜第7の実施の形態では、推定対象者が人間である場合を例に説明したが、これに限定されるものではなく、動物などを推定対象として、動物の意識低下状態を推定するようにしてもよい。 Moreover, in said 1st Embodiment-7th Embodiment, although the case where an estimation object person was a human was demonstrated to the example, it is not limited to this, An animal etc. are made into an estimation object, You may make it estimate the consciousness fall state of an animal.
また、上記の第1の実施の形態〜第7の実施の形態で説明した、心拍と呼吸との相関性を示す指標、呼吸のばらつき度合いを示す指標を組み合わせて、推定対象者の意識低下状態を推定するようにしてもよい。 Further, by combining the index indicating the correlation between heartbeat and respiration and the index indicating the degree of variation in respiration described in the first to seventh embodiments, the state of reduced consciousness of the estimation target person May be estimated.
また、相関関数の値として、リサンプリングされた心拍間隔の時系列データに対する、直流成分を除去した呼吸関連生体信号の時系列データの位相遅れが任意値を取る時の相関係数の絶対値を使用するようにしてもよい。 In addition, as the value of the correlation function, the absolute value of the correlation coefficient when the phase delay of the time series data of the respiratory related biological signal from which the DC component is removed with respect to the time series data of the resampled heartbeat interval takes an arbitrary value. It may be used.
なお、本発明のプログラムは、記録媒体に格納して提供することができる。 The program of the present invention can be provided by being stored in a recording medium.
10、210、310、610、710 意識状態推定装置
12、212 心拍検出装置
14 呼吸検出装置
16、216、316、716コンピュータ
20、220 心拍センサ
22 心拍間隔算出部
26 呼吸センサ
28 低域遮断フィルタ
32 相関係数算出部
34 パワースペクトル密度算出部
36 パワースペクトル密度特徴量算出部
38、338、738 意識状態推定部
222 低域通過フィルタ
232 相互情報量算出部
334 呼吸ばらつき度算出部
612 体動検出装置
614 発話検出装置
10, 210, 310, 610, 710 Consciousness
Claims (11)
前記データ検出手段によって検出された前記心拍間隔の時系列データ及び前記呼吸情報の時系列データの相関性を判定する相関性判定手段と、
前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、
を含む意識状態推定装置。 Data detection means for detecting time series data of a heartbeat interval of the estimation target person and time series data of respiratory information indicating a physical quantity that changes in accordance with the breathing of the estimation target person;
Correlation determining means for determining the correlation between the time series data of the heartbeat interval and the time series data of the respiratory information detected by the data detecting means;
Consciousness state estimating means for estimating a state of reduced consciousness of the estimation subject so that the higher the correlation determined by the correlation determining means, the lower the consciousness state;
Consciousness state estimation device.
前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段と、
前記信号検出手段によって分離された前記周波数帯域の信号、及び前記データ検出手段によって検出された前記呼吸情報の時系列データの相関性を判定する相関性判定手段と、
前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、
を含む意識状態推定装置。 A signal detecting means for detecting a signal that changes in accordance with the heart beat of the estimation target person, and separating a signal in a predetermined frequency band including a component derived from respiration from the detected signal;
Data detection means for detecting time-series data of respiration information indicating a physical quantity that changes in accordance with the respiration of the estimation target person;
Correlation determining means for determining the correlation between the signal of the frequency band separated by the signal detecting means and the time series data of the respiratory information detected by the data detecting means;
Consciousness state estimating means for estimating a state of reduced consciousness of the estimation subject so that the higher the correlation determined by the correlation determining means, the lower the consciousness state;
Consciousness state estimation device.
前記意識状態推定手段は、前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いと、前記相関性判定手段によって判定された相関性とに基づいて、前記推定対象者の意識低下状態を推定する請求項1又は2記載の意識状態推定装置。 Breathing variation determination means for determining the degree of variation in breathing based on the time series data of the breathing information,
The consciousness state estimation unit estimates a consciousness reduction state of the estimation target person based on the degree of breathing variation determined by the breathing variation determination unit and the correlation determined by the correlation determination unit. The consciousness state estimation apparatus according to claim 1 or 2.
前記意識状態推定手段は、前記相関性判定手段によって判定された相関性と、前記体動発話検出手段による前記体動及び発話の少なくとも一方の検出結果とに基づいて、前記意識低下状態を推定する請求項1〜請求項5の何れか1項記載の意識状態推定装置。 Further comprising body movement utterance detection means for detecting at least one of the body movement and utterance of the estimation target person,
The consciousness state estimation unit estimates the consciousness reduction state based on the correlation determined by the correlation determination unit and the detection result of at least one of the body movement and the utterance by the body movement utterance detection unit. The consciousness state estimation apparatus according to any one of claims 1 to 5.
前記データ検出手段によって検出された前記呼吸情報の時系列データに基づいて、呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する呼吸ばらつき判定手段と、
前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いが低いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段と、
を含む意識状態推定装置。 Data detection means for detecting time-series data of respiratory information indicating physical quantities that change in accordance with the estimation target person's breath;
Based on the time-series data of the respiratory information detected by the data detection means, the frequency dependence of the power spectrum density of the time-series data or the autocorrelation of the time-series data is determined as the degree of breathing variation. Breathing variation determination means;
Consciousness state estimating means for estimating a state of consciousness reduction of the estimation subject so that the consciousness state decreases as the degree of variation in breathing determined by the breathing variation determination means decreases;
Consciousness state estimation device.
前記意識状態推定手段は、前記呼吸ばらつき判定手段によって判定された前記呼吸のはらつき度合いと、前記体動発話検出手段による前記体動及び発話の少なくとも一方の検出結果とに基づいて、前記意識低下状態を推定する請求項7記載の意識状態推定装置。 Further comprising body movement utterance detection means for detecting at least one of the body movement and utterance of the estimation target person,
The consciousness state estimation unit is configured to reduce the consciousness based on the degree of fluctuation of the breath determined by the breathing variation determination unit and the detection result of at least one of the body movement and the utterance by the body movement utterance detection unit. The consciousness state estimation apparatus according to claim 7, wherein the state is estimated.
推定対象者の心拍間隔の時系列データ及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記心拍間隔の時系列データ及び前記呼吸情報の時系列データの相関性を判定する相関性判定手段、及び
前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段
として機能させるためのプログラム。 Computer
Time series data of the heart rate interval detected by data detecting means for detecting time series data of the heart rate interval of the estimation target person and time series data of respiration information indicating a physical quantity that changes according to the breathing of the estimation target person, and Correlation determining means for determining the correlation of respiratory information time-series data, and estimating the state of reduced consciousness of the person to be estimated such that the higher the correlation determined by the correlation determining means, the lower the consciousness state Program to function as a consciousness state estimation means.
推定対象者の心臓の拍動に応じて変化する信号を検出すると共に、前記検出された信号から、呼吸に由来する成分が含まれる予め定められた周波数帯域の信号を分離する信号検出手段によって分離された前記周波数帯域の信号、及び前記推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記呼吸情報の時系列データの相関性を判定する相関性判定手段、及び
前記相関性判定手段によって判定された相関性が高いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段
として機能させるためのプログラム。 Computer
A signal that changes in accordance with the heartbeat of the person to be estimated is detected, and separated from the detected signal by a signal detection unit that separates a signal in a predetermined frequency band including a component derived from respiration. Correlation of the time-series data of the respiration information detected by the data detection means for detecting the time-series data of the respiration information indicating the physical quantity that changes in accordance with the respiration of the estimated person and the estimation target person's respiration. Correlation determination means for determining, and a consciousness state estimation means for estimating a state of reduced consciousness of the estimation target person so that the consciousness state decreases as the correlation determined by the correlation determination means increases. program.
推定対象者の呼吸に応じて変化する物理量を示す呼吸情報の時系列データを検出するデータ検出手段によって検出された前記呼吸情報の時系列データに基づいて、呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する呼吸ばらつき判定手段、及び
前記呼吸ばらつき判定手段によって判定された前記呼吸のばらつき度合いが低いほど意識状態が低下するように、前記推定対象者の意識低下状態を推定する意識状態推定手段
として機能させるためのプログラム。 Computer
Based on the time series data of the respiratory information detected by the data detection means for detecting the time series data of the respiratory information indicating the physical quantity that changes in accordance with the breath of the estimation target person, the time series data Breathing variation determining means for determining the frequency dependence of the power spectral density of the device or the autocorrelation of the time series data, and the state of consciousness decreases as the degree of variation in breathing determined by the breathing variation determining means decreases. And a program for functioning as a consciousness state estimating means for estimating a consciousness reduction state of the estimation target person.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011148101A JP5929020B2 (en) | 2011-07-04 | 2011-07-04 | Consciousness state estimation device and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011148101A JP5929020B2 (en) | 2011-07-04 | 2011-07-04 | Consciousness state estimation device and program |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015232956A Division JP2016047305A (en) | 2015-11-30 | 2015-11-30 | Consciousness state estimation device and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013013542A true JP2013013542A (en) | 2013-01-24 |
JP5929020B2 JP5929020B2 (en) | 2016-06-01 |
Family
ID=47686827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011148101A Expired - Fee Related JP5929020B2 (en) | 2011-07-04 | 2011-07-04 | Consciousness state estimation device and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5929020B2 (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014193263A (en) * | 2013-02-27 | 2014-10-09 | Delta Tooling Co Ltd | Biological state estimation device, and computer program |
JP2014236332A (en) * | 2013-05-31 | 2014-12-15 | 日本放送協会 | Viewing state estimation device and program therefor |
JP2016510239A (en) * | 2013-02-05 | 2016-04-07 | セントル ホスピタリエ リージョナル ユニヴェルシテール ドゥ リールCentre Hospitalier Regional Universitaire de Lille | Method for measuring physiological parameters such as biological rhythm based on a measuring device associated with at least two sensors |
JP2016107095A (en) * | 2014-12-05 | 2016-06-20 | 公立大学法人首都大学東京 | Sleep stage determination device, sleep stage determination method, sleep stage determination program, and stress degree determination method |
JP2017121304A (en) * | 2016-01-06 | 2017-07-13 | 日本電信電話株式会社 | Information processing device, information processing method, and program |
JP2017537710A (en) * | 2014-12-11 | 2017-12-21 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | System and method for determining spectral boundaries of sleep stage classification |
KR101853102B1 (en) | 2016-05-16 | 2018-05-02 | 연세대학교 산학협력단 | Accelerometer based sleep sensing device |
WO2021187307A1 (en) * | 2020-03-17 | 2021-09-23 | データソリューションズ株式会社 | Living body abnormality detection device, living body abnormality detection method, and program |
JP2021527551A (en) * | 2018-06-19 | 2021-10-14 | マイカーディオ エルエルシー | Systems and methods for assessing health status or health status |
US11523760B2 (en) | 2019-05-21 | 2022-12-13 | Honda Motor Co., Ltd. | Arousal state estimation apparatus and arousal state estimation method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005319256A (en) * | 2004-05-07 | 2005-11-17 | Toyama Prefecture | Estimate method of parasympathetic activity |
JP2006026302A (en) * | 2004-07-21 | 2006-02-02 | Jepico Corp | Sleep analyzing apparatus |
-
2011
- 2011-07-04 JP JP2011148101A patent/JP5929020B2/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005319256A (en) * | 2004-05-07 | 2005-11-17 | Toyama Prefecture | Estimate method of parasympathetic activity |
JP2006026302A (en) * | 2004-07-21 | 2006-02-02 | Jepico Corp | Sleep analyzing apparatus |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016510239A (en) * | 2013-02-05 | 2016-04-07 | セントル ホスピタリエ リージョナル ユニヴェルシテール ドゥ リールCentre Hospitalier Regional Universitaire de Lille | Method for measuring physiological parameters such as biological rhythm based on a measuring device associated with at least two sensors |
JP2014193263A (en) * | 2013-02-27 | 2014-10-09 | Delta Tooling Co Ltd | Biological state estimation device, and computer program |
JP2014236332A (en) * | 2013-05-31 | 2014-12-15 | 日本放送協会 | Viewing state estimation device and program therefor |
JP2016107095A (en) * | 2014-12-05 | 2016-06-20 | 公立大学法人首都大学東京 | Sleep stage determination device, sleep stage determination method, sleep stage determination program, and stress degree determination method |
JP2017537710A (en) * | 2014-12-11 | 2017-12-21 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | System and method for determining spectral boundaries of sleep stage classification |
JP2017121304A (en) * | 2016-01-06 | 2017-07-13 | 日本電信電話株式会社 | Information processing device, information processing method, and program |
KR101853102B1 (en) | 2016-05-16 | 2018-05-02 | 연세대학교 산학협력단 | Accelerometer based sleep sensing device |
JP2021527551A (en) * | 2018-06-19 | 2021-10-14 | マイカーディオ エルエルシー | Systems and methods for assessing health status or health status |
JP7253047B2 (en) | 2018-06-19 | 2023-04-05 | マイカーディオ エルエルシー | Systems and methods for assessment of health conditions or health conditions |
US11523760B2 (en) | 2019-05-21 | 2022-12-13 | Honda Motor Co., Ltd. | Arousal state estimation apparatus and arousal state estimation method |
WO2021187307A1 (en) * | 2020-03-17 | 2021-09-23 | データソリューションズ株式会社 | Living body abnormality detection device, living body abnormality detection method, and program |
Also Published As
Publication number | Publication date |
---|---|
JP5929020B2 (en) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5929020B2 (en) | Consciousness state estimation device and program | |
EP3334337B1 (en) | Monitoring of sleep phenomena | |
CN105142515B (en) | The method for determining the mankind's sleep stage for being conducive to waking up | |
JP3733133B2 (en) | Sleep state estimation device | |
JP6190466B2 (en) | Biological signal measuring instrument and contact state estimation method | |
US20200121207A1 (en) | Method of processing a signal representing a physiological rhythm | |
US20140088373A1 (en) | System and method for determining sleep stage | |
JP6602657B2 (en) | Biological information acquisition apparatus, biological information acquisition program, and biological information acquisition method | |
JP2006020810A (en) | Device and program for estimating state of sleeping | |
JP6627112B2 (en) | Biological function testing device, method of operating biological function testing device, and program | |
EP3387988A1 (en) | Biological state estimation device, biological state estimation method, computer program, and recording medium | |
JP2006271474A (en) | Apparatus for estimating sleeping state | |
JP2018082931A (en) | Arousal level processing method and arousal level processing apparatus | |
JP2016047305A (en) | Consciousness state estimation device and program | |
WO2016057806A1 (en) | Weaning readiness indicator, sleeping status recording device, and air providing system applying nonlinear time-frequency analysis | |
KR101696787B1 (en) | Sleep-disordered breathing detection method and apparatus with continuous blood pressure signal measured during sleep | |
JP6372824B2 (en) | CVHR shape measuring device | |
JP2017124014A (en) | Drowsiness detection program, drowsiness detection method, and drowsiness detection device | |
JP6588035B2 (en) | Biological condition analyzer and computer program | |
JP6557489B2 (en) | Biological state estimation device and computer program | |
JP2020073108A (en) | Sleep stage determination method, sleep stage determination device and sleep stage determination program | |
JP4868514B2 (en) | Apparatus and method for outputting result of estimation of biological state | |
JP6553397B2 (en) | Sleep feeling determination method and sleep feeling determination apparatus | |
WO2023048158A1 (en) | Sleep apnea syndrome determination device, sleep apnea syndrome determination method, and program | |
JP7158641B1 (en) | Apnea hypopnea index estimation device, method and program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140513 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150121 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150224 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150427 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150929 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20151130 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160405 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160418 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5929020 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |