JP2013013542A - Apparatus and program for estimating state of consciousness - Google Patents

Apparatus and program for estimating state of consciousness Download PDF

Info

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
Application number
JP2011148101A
Other languages
Japanese (ja)
Other versions
JP5929020B2 (en
Inventor
Takahiro Yamaguchi
喬弘 山口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Central R&D Labs Inc
Original Assignee
Toyota Central R&D Labs Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Central R&D Labs Inc filed Critical Toyota Central R&D Labs Inc
Priority to JP2011148101A priority Critical patent/JP5929020B2/en
Publication of JP2013013542A publication Critical patent/JP2013013542A/en
Application granted granted Critical
Publication of JP5929020B2 publication Critical patent/JP5929020B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

PROBLEM TO BE SOLVED: To estimate stably in good accuracy, the consciousness decreased condition of a person to be estimated.SOLUTION: The time-series data of the heartbeat interval of the person to be estimated is detected by a heartbeat detection device 12, and the time-series data of the biological signal related to breath is detected by a breath detection device 14. The correlation coefficient between the time-series data of the heartbeat interval and the time-series data of the biological signal related to breath is calculated by a correlation coefficient calculating unit 32. A power spectral density characteristic amount is calculated by a power spectral density characteristic amount calculation section 36 as an indicator that shows the variation degree of breath. The consciousness-decreased condition of the person is estimated by a state of consciousness estimation section 38 based on the power spectral density characteristic and the correlation coefficient that have been calculated.

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).

特開2006-271474号公報JP 2006-271474 A 特開2006−280686号公報JP 2006-280686 A 特表2008−504942号公報Special table 2008-504942 gazette 特開2008−154681号公報JP 2008-154681 A

心拍変動の低周波成分は副交感神経活動を反映し、心拍変動の高周波成分は副交感神経活動を反映していることは古くから知られており、睡眠状態の指標として用いられることが検討されている。しかしながら、心拍変動は副交感神経活動以外の要因(日内変動等)でも変動することが知られており、上記の特許文献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.

本発明の第1の実施の形態に係る意識状態推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the consciousness state estimation apparatus which concerns on the 1st Embodiment of this invention. 心拍間隔と呼吸との相関性を説明するための図である。It is a figure for demonstrating the correlation with a heartbeat interval and respiration. 呼吸のパワースペクトル密度と周波数との関係を示すグラフである。It is a graph which shows the relationship between the power spectrum density and frequency of respiration. 本発明の第1の実施の形態に係る意識状態推定装置における意識状態推定処理ルーチンの内容を示すフローチャートである。It is a flowchart which shows the content of the consciousness state estimation process routine in the consciousness state estimation apparatus which concerns on the 1st Embodiment of this invention. 本発明の第2の実施の形態に係る意識状態推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the consciousness state estimation apparatus which concerns on the 2nd Embodiment of this invention. 本発明の第3の実施の形態に係る意識状態推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the consciousness state estimation apparatus which concerns on the 3rd Embodiment of this invention. 呼吸関連生体信号の時系列データを示す図である。It is a figure which shows the time series data of a respiratory related biological signal. 呼吸パワースペクトル密度と相関係数の絶対値とを2軸とした座標系を示す図である。It is a figure which shows the coordinate system which made the respiratory power spectral density and the absolute value of the correlation coefficient 2 axes. 本発明の第6の実施の形態に係る意識状態推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the consciousness state estimation apparatus which concerns on the 6th Embodiment of this invention. 本発明の第7の実施の形態に係る意識状態推定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the consciousness state estimation apparatus which concerns on the 7th Embodiment of this invention. 本発明の第7の実施の形態に係る意識状態推定装置における意識状態推定処理ルーチンの内容を示すフローチャートである。It is a flowchart which shows the content of the consciousness state estimation process routine in the consciousness state estimation apparatus which concerns on the 7th Embodiment of this invention.

以下、図面を参照して本発明の実施の形態を詳細に説明する。   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 state estimation device 10 according to the first embodiment includes a heartbeat detection device 12 that detects a heartbeat of an estimation target person, a respiration detection device 14 that detects a respiration of the estimation target person, Based on the time-series data of the heartbeat interval output from the heartbeat detection device 12 and the time-series data of the respiratory related biological signal output from the respiration detection device 14, it is estimated whether the estimation target person is in a state of reduced consciousness. And a computer 16 that causes the output unit 18 to output the estimation result. The heartbeat detection device 12 and the respiration detection device 14 are examples of data detection means.

心拍検出装置12は、心拍センサ20、心拍間隔算出部22、及び心拍間隔リサンプリング部24を備えている。心拍センサ20は、人体より心拍に関連する生体信号を検出するセンサであり、心拍関連生体信号の時系列データを出力する。心拍センサ20として、例えば、心電計、心音計、脈波センサ、電気インピーダンス計測装置、電磁誘導的に電気インピーダンスを計測する装置、圧電素子を用いて心拍由来の生体振動を検出する装置、又は電磁波により心拍由来の生体振動を検出する装置等を用いればよい。   The heartbeat detection device 12 includes a heartbeat sensor 20, a heartbeat interval calculation unit 22, and a heartbeat interval resampling unit 24. The heart rate sensor 20 is a sensor that detects a biological signal related to the heartbeat from the human body, and outputs time-series data of the heartbeat related biological signal. As the heart rate sensor 20, for example, an electrocardiograph, a heart sound meter, a pulse wave sensor, an electrical impedance measurement device, a device that measures electrical impedance electromagnetically, a device that detects biological vibrations derived from a heart using a piezoelectric element, or A device or the like that detects biological vibrations derived from heartbeats using electromagnetic waves may be used.

心拍間隔算出部22は、心拍センサ20により出力された心拍関連生体信号の時系列データに基づいて、心拍時刻を算出し、心拍間隔の時系列データを出力する。心拍間隔リサンプリング部24は、心拍間隔算出部22から出力された心拍間隔の時系列データを、後述する呼吸センサ26のサンプリング時刻に合わせてリサンプリングして出力する。   The heartbeat interval calculator 22 calculates a heartbeat time based on the time-series data of the heartbeat-related biological signal output from the heartbeat sensor 20, and outputs the time-series data of the heartbeat interval. The heartbeat interval resampling unit 24 resamples and outputs the time series data of the heartbeat interval output from the heartbeat interval calculation unit 22 in accordance with the sampling time of the respiration sensor 26 described later.

呼吸検出装置14は、呼吸センサ26及び低域遮断フィルタ28により構成される。呼吸センサ26は、人体から呼吸に関連する生体信号(例えば、呼吸動を示す呼吸波形)を検出するセンサであり、呼吸関連生体信号の時系列データを出力する。呼吸センサ26として、例えば、胸郭周囲長を測定する装置、電気インピーダンス計測装置、電磁誘導的に電気インピーダンスを計測する装置、圧電素子を用いて呼吸由来の生体振動を検出する装置、又は電磁波により呼吸由来の生体振動を検出する装置等を用いればよい。また、心拍センサ20が心電図を出力する場合、心電図R波の振幅より呼吸波形を推定する手法を用いてもよい。なお、呼吸関連生体信号が、呼吸に応じて変化する物理量を示す呼吸情報の一例である。   The respiration detection device 14 includes a respiration sensor 26 and a low-frequency cutoff filter 28. The respiration sensor 26 is a sensor that detects a biological signal related to respiration (for example, a respiration waveform indicating respiration movement) from the human body, and outputs time-series data of the respiration related biosignal. As the respiration sensor 26, for example, a device for measuring the circumference of the thorax, an electric impedance measuring device, a device for measuring electric impedance by electromagnetic induction, a device for detecting biological vibration derived from respiration using a piezoelectric element, or breathing by electromagnetic waves A device for detecting the biological vibration derived from the origin may be used. When the heart rate sensor 20 outputs an electrocardiogram, a method of estimating a respiratory waveform from the amplitude of the electrocardiogram R wave may be used. Note that the respiration-related biological signal is an example of respiration information indicating a physical quantity that changes according to respiration.

低域遮断フィルタ28は、呼吸関連生体信号の時系列データより直流成分を除去し、直流成分を除去した呼吸関連生体信号の時系列データを出力する。   The low-frequency cutoff filter 28 removes the direct current component from the time series data of the respiratory related biological signal and outputs the time series data of the respiratory related biological signal from which the direct current component has been removed.

コンピュータ16は、CPUと、RAMと、後述する意識状態推定処理ルーチンを実行するためのプログラムを記憶したROMとを備え、機能的には次に示すように構成されている。コンピュータ16は、心拍検出装置12から出力される心拍間隔の時系列データ及び呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、心拍間隔の時系列データと呼吸関連生体信号の時系列データとの相関係数を算出する相関係数算出部32と、呼吸関連生体信号の時系列データに基づいて、呼吸関連生体信号のパワースペクトル密度を算出するパワースペクトル密度算出部34と、算出されたパワースペクトル密度に基づいて、パワースペクトル密度特徴量を算出するパワースペクトル密度特徴量算出部36と、算出された相関係数及びパワースペクトル密度特徴量に基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部38とを備えている。なお、相関係数算出部32が、相関性判定手段の一例である。パワースペクトル密度算出部34及びパワースペクトル密度特徴量算出部36が、呼吸ばらつき判定手段の一例である。   The computer 16 includes a CPU, a RAM, and a ROM that stores a program for executing a consciousness state estimation processing routine described later, and is functionally configured as follows. The computer 16 includes a data acquisition unit 30 that acquires time-series data of heartbeat intervals output from the heartbeat detection device 12 and time-series data of respiratory-related biological signals output from the respiration detection device 14, and time-series data of heartbeat intervals. Coefficient calculating unit 32 for calculating the correlation coefficient between the time series data of the respiratory related biological signal and the power spectrum for calculating the power spectrum density of the respiratory related biological signal based on the time series data of the respiratory related biological signal Based on the density calculation unit 34, the power spectrum density feature amount calculation unit 36 that calculates the power spectrum density feature amount based on the calculated power spectrum density, and the calculated correlation coefficient and power spectrum density feature amount, A consciousness state estimation unit for estimating whether or not the estimation target person is in a consciousness lowered state. The correlation coefficient calculation unit 32 is an example of a correlation determination unit. The power spectrum density calculation unit 34 and the power spectrum density feature amount calculation unit 36 are an example of respiratory variation determination means.

データ取得部30は、心拍検出装置12から入力された心拍間隔の時系列データを取得し、メモリ(図示省略)に記憶する共に、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得してメモリに記憶する。   The data acquisition unit 30 acquires time-series data of heartbeat intervals input from the heartbeat detection device 12, stores the time-series data in a memory (not shown), and time-series data of respiratory related biological signals output from the respiration detection device 14. Is stored in the memory.

ここで、本実施の形態の原理について説明する。   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 coefficient calculation unit 32 re-samples the time series data of the heartbeat interval and the time series data of the respiratory related biological signal from which the DC component has been removed, at a predetermined time length. The correlation coefficient between the time series data of the heartbeat interval and the time series data of the respiratory related biological signal is calculated as follows.

まず、ある時間窓に含まれる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).

Figure 2013013542
Figure 2013013542

また、同じ時間窓に含まれる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).

Figure 2013013542
Figure 2013013542

ただし、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).

Figure 2013013542
Figure 2013013542

以上より、求められた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).

Figure 2013013542
Figure 2013013542

ただし、 ̄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 density calculation unit 34 performs a discrete Fourier transform on the time-series data of the respiratory related biological signal from which the DC component is removed for a predetermined time length, and calculates the power spectrum density.

ここで、呼吸または呼吸間隔のばらつきとパワースペクトル密度との関係の模式図を、図3に示す。睡眠深度が深い状態での呼吸間隔は、個体差や時間帯によって差があるものの、ある程度一定の範囲内に収まると考えられる。その範囲の最長間隔から求められた呼吸周波数をf、最短間隔から求められた周波数をfとすると、周波数fからfまでのパワースペクトル密度の和が求められる。この和が全帯域のパワースペクトル密度の和に占める割合は、呼吸周波数がfからfまでの間で安定しているかどうかの指標となり、呼吸の安定性を示すパワースペクトル密度の特徴量として使用できる。 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によって、算出されたパワースペクトル密度に基づいて、パワースペクトル密度特徴量として、予め定められた周波数帯域(fからfまで)におけるパワースペクトル密度の和が、パワースペクトル密度全体に占める割合を算出する。なお、パワースペクトル密度特徴量が、パワースクペクトル密度の周波数依存性の一例である。 Therefore, the power spectral density characteristic amount calculation unit 36, based on the power spectral density calculated as the power spectral density characteristic amount, the sum of the power spectral density in the predetermined frequency band (from f 1 to f 2) The ratio of the total power spectrum density is calculated. The power spectral density feature amount is an example of the frequency dependence of the power spectrum density.

意識状態推定部38は、算出された相関係数及びパワースペクトル密度特徴量に基づいて、以下に説明するように、推定対象者の意識低下状態を推定する。   Based on the calculated correlation coefficient and power spectral density feature quantity, the consciousness state estimation unit 38 estimates the consciousness reduction state of the estimation target person as described below.

まず、相関係数と意識低下状態との関連性について述べる。心拍と呼吸との関連については、吸気時に心拍数が上昇(心拍間隔が短縮)し、呼気時に心拍数が下降(心拍間隔が延長)することがわかっている(上記図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 state estimation unit 38 calculates 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 and the power spectrum density feature amount. The product after arbitrary weighting is calculated as an index for estimating the consciousness level. This index increases as the level of consciousness decreases (towards the sleep direction) and decreases as it increases (towards the awakening direction).

意識状態推定部38は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18により推定結果を出力させる。   The consciousness state estimation unit 38 estimates that the estimation target person is in a reduced consciousness state when the calculated index is greater than or equal to the threshold value, and otherwise, the estimation target person is not in a consciousness reduction state. presume. In addition, the consciousness state estimation unit 38 causes the output unit 18 to output an estimation result.

なお、上記の閾値として、従来の医学的治験により、汎用性のある値をあらかじめ与えることができる。また、個人差に対応した閾値を設定するため、次に述べるような方法で定めることもできる。まず、対象となる人の就寝中の心拍間隔の時系列データ及び呼吸関連生体信号の時系列データを記録する。この時系列データには、安静覚醒状態、覚醒から睡眠への移行状態、睡眠状態の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 state estimation apparatus 10 according to the first embodiment will be described.

まず、心拍検出装置12によって推定対象者の心拍間隔を検出し、リサンプリングされた心拍間隔の時系列データをコンピュータ16に入力する。また、呼吸検出装置14によって、推定対象者の呼吸関連生体信号を検出し、直流成分を除去した呼吸関連生体信号の時系列データをコンピュータ16に入力する。また、コンピュータ16において、図4に示す意識状態推定処理ルーチンが繰り返し実行される。   First, the heartbeat detection device 12 detects the heartbeat interval of the estimation target person, and inputs time-series data of the resampled heartbeat interval to the computer 16. In addition, the respiratory detection device 14 detects a respiratory related biological signal of the estimation target person, and inputs time series data of the respiratory related biological signal from which the direct current component has been removed to the computer 16. In the computer 16, the consciousness state estimation process routine shown in FIG. 4 is repeatedly executed.

まず、ステップ100で、直近の所定時間分の、リサンプリングされた心拍間隔の時系列データ及び直流成分を除去した呼吸関連生体信号の時系列データを取得する。   First, in step 100, the time-series data of the resampled heartbeat interval and the time-series data of the respiratory related biological signal from which the direct current component is removed are acquired for the most recent predetermined time.

そして、ステップ102において、上記ステップ100で取得した心拍間隔の時系列データ及び呼吸関連生体信号の時系列データについて、相関係数を算出する。   In step 102, a correlation coefficient is calculated for the time series data of the heartbeat interval and the time series data of the respiratory related biological signal acquired in step 100.

次のステップ104では、上記ステップ100で取得した呼吸関連生体信号の時系列データについて、パワースペクトル密度を算出し、ステップ106において、パワースペクトル密度特徴量として、予め定められた周波数帯域(f〜f)におけるパワースペクトル密度の和の割合を算出する。 In the next step 104, the power spectrum density is calculated for the time-series data of the respiratory related biological signal acquired in step 100, and in step 106, a predetermined frequency band (f 1 to The ratio of the sum of the power spectral densities in f 2 ) is calculated.

そして、ステップ108において、上記ステップ102で算出した相関係数、及び上記ステップ106で算出したパワースペクトル密度特徴量に基づいて、意識レベルを推定するための指標を算出する。次のステップ110では、上記ステップ108で算出された指標が、閾値以上であるか否かに基づいて、推定対象者が意識低下状態であるか否かを推定し、ステップ112において、出力部18によって推定結果を出力させて、意識状態推定処理ルーチンを終了する。   In step 108, an index for estimating the consciousness level is calculated based on the correlation coefficient calculated in step 102 and the power spectrum density feature amount calculated in step 106. In the next step 110, it is estimated whether or not the estimation target person is in a state of reduced consciousness based on whether or not the index calculated in step 108 is equal to or greater than a threshold. In step 112, the output unit 18 To output the estimation result, and the consciousness state estimation processing routine is terminated.

以上説明したように、第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 heartbeat detection device 212 of the consciousness state estimation device 210 according to the second embodiment includes a heartbeat sensor 220, a low-pass filter 222, and a heartbeat waveform resampling unit 224. The heartbeat sensor 220 is a sensor that detects a biological signal (electrocardiogram waveform) related to the heartbeat from the human body, and outputs time-series data of the heartbeat-related biological signal indicating the electrocardiogram waveform. The heartbeat detection device 212 is an example of a signal detection unit, and a heartbeat-related biological signal indicating an electrocardiogram waveform is an example of a signal that changes according to the heartbeat. The respiration detection device 14 is an example of a data detection unit.

低域通過フィルタ222のカットオフ周波数は,平均的な心拍周波数及び呼吸周波数を考慮して決定されるが、心拍の周波数より高い周波数で設定することが望ましい。低域通過フィルタ222によって、呼吸に由来する成分が含まれる周波数帯域の信号が分離される。   The cut-off frequency of the low-pass filter 222 is determined in consideration of an average heartbeat frequency and a respiration frequency, but is preferably set at a frequency higher than the heartbeat frequency. The low-pass filter 222 separates a signal in a frequency band including a component derived from respiration.

心拍波形リサンプリング部224は、高域遮断された心拍関連生体信号の時系列データを、呼吸センサ26のサンプリング時刻に合わせてリサンプリングして出力する。   The heartbeat waveform resampling unit 224 resamples and outputs the time-series data of the heartbeat-related biological signal that has been cut off in the high frequency in accordance with the sampling time of the respiration sensor 26.

コンピュータ216は、心拍検出装置212から出力される心拍関連生体信号の時系列データ及び呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、心拍関連生体信号の時系列データと呼吸関連生体信号の時系列データとの間の相互情報量を算出する相互情報量算出部232と、パワースペクトル密度算出部34と、パワースペクトル密度特徴量算出部36と、算出された相互情報量及びパワースペクトル密度特徴量に基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部38とを備えている。なお、相互情報量算出部232が、相関性判定手段の一例であり、相互情報量が、相関性の一例である。   The computer 216 includes a data acquisition unit 30 that acquires time-series data of heartbeat-related biological signals output from the heartbeat detection device 212 and time-series data of respiratory-related biological signals output from the respiration detection device 14, and a heartbeat-related biological signal. A mutual information calculation unit 232, a power spectrum density calculation unit 34, a power spectrum density feature amount calculation unit 36, and a calculation And a consciousness state estimation unit 38 for estimating whether or not the estimation target person is in a consciousness reduction state based on the mutual information amount and the power spectral density feature amount. Note that the mutual information amount calculation unit 232 is an example of a correlation determination unit, and the mutual information amount is an example of a correlation.

相互情報量算出部232は、リサンプリングされた心拍関連生体信号の時系列データ及び直流成分を除去した呼吸関連生体信号の時系列データのうち、予め定められた時間長における、心拍関連生体信号の時系列データ及び呼吸関連生体信号の時系列データ間の相互情報量を、以下のように算出する。   The mutual information amount calculation unit 232 includes the time-series data of the resampled heartbeat-related biological signal and the time-series data of the respiratory-related biological signal from which the direct current component has been removed. The mutual information amount between the time series data and the time series data of the respiratory related biological signal is calculated as follows.

まず、ある時間窓に含まれる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).

Figure 2013013542
Figure 2013013542

また、同じ時間窓に含まれる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).

Figure 2013013542
Figure 2013013542

ただし、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).

Figure 2013013542
Figure 2013013542

以上より求められた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).

Figure 2013013542
Figure 2013013542

上記(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 consciousness state estimator 38 arbitrarily selects the mutual information amount of the time-series data of the respiratory-related biological signal from which the DC component is removed and the power spectrum density feature amount with respect to the time-series data of the resampled heartbeat-related biological signal. Is calculated as an index for estimating the level of consciousness.

意識状態推定部38は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。   The consciousness state estimation unit 38 estimates that the estimation target person is in a reduced consciousness state when the calculated index is greater than or equal to the threshold value, and otherwise, the estimation target person is not in a consciousness reduction state. presume.

なお、第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 computer 316 of the consciousness state estimation apparatus 310 according to the third embodiment is based on the data acquisition unit 30, the correlation coefficient calculation unit 32, and time-series data of respiratory related biological signals. , A respiratory variability calculator 334 that calculates the degree of respiratory variability, and a consciousness state estimator 338 that estimates whether or not the estimation target person is in a state of reduced consciousness based on the calculated correlation coefficient and respiratory variability. And.

呼吸ばらつき度算出部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 degree calculation unit 334 calculates the time-series data of the respiratory-related biological signal as described below. An autocorrelation function is calculated and used as the degree of respiratory variability.

まず、ある時間窓における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).

Figure 2013013542
Figure 2013013542

ただし、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).

Figure 2013013542
Figure 2013013542

ただし、−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 state estimation unit 338 calculates 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 and the calculated autocorrelation function with respect to the time-series data of the resampled heartbeat interval. The product after arbitrary weighting is calculated as an index for estimating the consciousness level. This index increases as the level of consciousness decreases (towards the sleep direction) and decreases as it increases (towards the awakening direction).

意識状態推定部338は、算出した上記の指標が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18によって、推定結果を出力させる。   The consciousness state estimation unit 338 estimates that the estimation target person is in a state of reduced consciousness when the calculated index is equal to or greater than the threshold value, and otherwise, the estimation target person is not in a consciousness reduction state. presume. In addition, the consciousness state estimation unit 38 causes the output unit 18 to output an estimation result.

なお、第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 state estimation device 310 according to the fourth embodiment, based on the time-series data of the respiratory related biological signal from which the DC component is removed by the respiratory variation degree calculation unit 334 for a predetermined time length, As described above, the coefficient of variation of the breathing interval is calculated as the degree of breathing variation.

図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).

Figure 2013013542
Figure 2013013542

上記(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 state estimation device 10 according to the fifth embodiment, the consciousness state estimation unit 338 determines, based on the position on the coordinates spanned by the two axes of the absolute value of the correlation coefficient and the power spectrum density feature amount, It is estimated whether the estimation target person is in a state of reduced consciousness. A schematic diagram of the coordinates stretched at this time is shown in FIG. Note that the boundary between the consciousness lowered state (sleep state) and the other states can be determined by the same method as the method described in the first embodiment.

なお、第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 state estimation device 610 according to the sixth embodiment includes a heartbeat detection device 12, a respiration detection device 14, a body motion detection device 612 that detects the body motion of the estimation target person, An utterance detection device 614 that detects the utterance of the estimation target person and a computer 16 are provided.

体動検出装置612としては、例えば、人間に装着するセンサとしての加速度センサを用いることができ、あるいは、人間の周辺環境に装着するセンサとしての、人間と接触する部位の圧力を計測するセンサ、カメラ、電波式距離センサ、電波式速度センサ、音波式距離センサ、音波式速度センサ等を用いることができる。また、心拍センサ20及び呼吸センサ26とセンサを共有することも可能である。   As the body motion detection device 612, for example, an acceleration sensor as a sensor worn on a human can be used, or a sensor for measuring the pressure of a part in contact with a human as a sensor worn on a human surrounding environment, A camera, a radio wave distance sensor, a radio wave speed sensor, a sound wave distance sensor, a sound wave speed sensor, or the like can be used. It is also possible to share the sensor with the heart rate sensor 20 and the respiration sensor 26.

また、発話検出装置614としては、マイクロフォンなどを用いることができる。   As the speech detection device 614, a microphone or the like can be used.

コンピュータ16のデータ取得部30は、心拍検出装置12から入力された心拍間隔の時系列データを取得し、メモリ(図示省略)に記憶する共に、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得してメモリに記憶する。また、データ取得部30は、体動検出装置612から入力された体動検出結果の時系列データを取得し、メモリ(図示省略)に記憶する共に、発話検出装置614から出力される発話検出結果の時系列データを取得してメモリに記憶する。   The data acquisition unit 30 of the computer 16 acquires time-series data of heartbeat intervals input from the heartbeat detection device 12 and stores it in a memory (not shown), as well as respiratory-related biological signals output from the respiration detection device 14. Time series data is acquired and stored in memory. In addition, the data acquisition unit 30 acquires time series data of the body motion detection result input from the body motion detection device 612, stores the time series data in a memory (not shown), and outputs the speech detection result output from the speech detection device 614. Is acquired and stored in the memory.

ここで、本実施の形態の原理について説明する。   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 state estimation unit 38, and the power Based on the spectral density feature amount, an index for estimating the consciousness level is calculated. When the calculated index is equal to or greater than the threshold value, it is estimated that the estimation target person is in a state of reduced consciousness. If the calculated index is less than the threshold, whether or not at least one of the body motion and the utterance has a certain intensity or more based on the time series data of the body motion detection result and the time series data of the speech detection result Determine whether. If at least one of body movement and utterance has a certain level of intensity (when at least one of body movement and utterance is detected), it is determined that the breathing variation is caused by awakening body movement or utterance. The estimation target person's consciousness state is estimated to be an awakening state.

また、意識状態推定部38は、算出した上記の指標が閾値未満である場合であって、かつ、体動検出結果の時系列データ及び発話検出結果の時系列データに基づいて体動及び発話の何れもが検出されない場合、推定対象者の意識状態が、覚醒から睡眠の過渡状態であると推定する。   In addition, the consciousness state estimation unit 38 is a case where the calculated index is less than the threshold, and based on the time series data of the body motion detection result and the time series data of the speech detection result, If neither is detected, it is estimated that the estimation target person's consciousness state is a transitional state from awakening to sleep.

なお、第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 state estimation device 710 according to the seventh embodiment is based on the respiratory detection device 14 and the time-series data of the respiratory related biological signal output from the respiratory detection device 14. A computer 716 that estimates whether or not the person is in a state of reduced consciousness and causes the output unit 18 to output the estimation result.

コンピュータ716は、呼吸検出装置14から出力される呼吸関連生体信号の時系列データを取得するデータ取得部30と、呼吸関連生体信号の時系列データに基づいて、呼吸のばらつき度合いを算出する呼吸ばらつき度算出部334と、算出された呼吸ばらつき度合いに基づいて、推定対象者が意識低下状態であるか否かを推定する意識状態推定部738とを備えている。   The computer 716 obtains the time series data of the respiratory related biological signal output from the respiratory detection device 14, and calculates the degree of respiratory variation based on the time series data of the respiratory related biological signal. A degree calculation unit 334 and a consciousness state estimation unit 738 that estimates whether or not the estimation target person is in a consciousness reduction state based on the calculated degree of respiratory variation.

意識状態推定部738は、呼吸ばらつき度として算出された自己相関関数R(i)のi-=0以外の最大のピーク値を、意識レベルを推定する指標として、算出した上記の指標(自己相関関数R(i)のi-=0以外の最大のピーク値)が閾値以上である場合に、推定対象者が意識低下状態であると推定し、それ以外の場合には、推定対象者は意識低下状態ではないと推定する。また、意識状態推定部38は、出力部18によって、推定結果を出力させる。   The consciousness state estimator 738 uses the maximum peak value other than i− = 0 of the autocorrelation function R (i) calculated as the degree of respiratory variability as an index for estimating the consciousness level. When the maximum peak value other than i− = 0 of the function R (i) is equal to or greater than the threshold value, the estimation target person is estimated to be in a state of reduced consciousness; otherwise, the estimation target person is conscious Estimated not in a degraded state. In addition, the consciousness state estimation unit 38 causes the output unit 18 to output an estimation result.

次に、第7の実施の形態に係る意識状態推定装置710の作用について説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。   Next, the operation of the consciousness state estimation apparatus 710 according to the seventh embodiment will be described. In addition, about the process similar to 1st Embodiment, the same code | symbol is attached | subjected and detailed description is abbreviate | omitted.

まず、呼吸検出装置14によって、推定対象者の呼吸関連生体信号を検出し、直流成分を除去した呼吸関連生体信号の時系列データをコンピュータ716に入力する。また、コンピュータ716において、図11に示す意識状態推定処理ルーチンが繰り返し実行される。   First, the respiratory detection device 14 detects a respiratory related biological signal of the estimation target person, and inputs time series data of the respiratory related biological signal from which the DC component is removed to the computer 716. In the computer 716, the consciousness state estimation processing routine shown in FIG. 11 is repeatedly executed.

まず、ステップ750で、直近の所定時間分の、直流成分を除去した呼吸関連生体信号の時系列データを取得する。   First, in step 750, time series data of respiratory related biological signals from which a direct current component has been removed for the latest predetermined time is acquired.

そして、ステップ752において、上記ステップ750で取得した呼吸関連生体信号の時系列データについて、上記(10)式に従って自己相関関数R(i)のi-=0以外の最大のピーク値を算出して、意識レベルを推定するための指標とする。   Then, in step 752, the maximum peak value other than i− = 0 of the autocorrelation function R (i) is calculated for the time series data of the respiratory related biological signal acquired in step 750 according to the above equation (10). It is an index for estimating the level of consciousness.

次のステップ754では、上記ステップ752で算出された自己相関関数R(i)のi-=0以外の最大のピーク値が、閾値以上であるか否かに基づいて、推定対象者が意識低下状態であるか否かを推定し、ステップ112において、出力部18によって推定結果を出力させて、意識状態推定処理ルーチンを終了する。   In the next step 754, the estimation target person loses consciousness based on whether or not the maximum peak value other than i− = 0 of the autocorrelation function R (i) calculated in step 752 is greater than or equal to a threshold value. In step 112, the output unit 18 outputs an estimation result, and the consciousness state estimation processing routine is terminated.

以上説明したように、第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 state estimation device 12, 212 Heart rate detection device 14 Respiration detection device 16, 216, 316, 716 Computer 20, 220 Heart rate sensor 22 Heart rate interval calculation unit 26 Respiration sensor 28 Low-frequency cutoff filter 32 Correlation coefficient calculation unit 34 Power spectrum density calculation unit 36 Power spectrum density feature amount calculation unit 38, 338, 738 Consciousness state estimation unit 222 Low-pass filter 232 Mutual information amount calculation unit 334 Respiration variation degree calculation unit 612 Body motion detection device 614 Utterance detection device

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.
前記呼吸ばらつき判定手段は、前記呼吸情報の時系列データに基づいて、前記呼吸のばらつき度合いとして、呼吸間隔の標準偏差又は前記呼吸間隔の変動係数を算出する請求項3記載の意識状態推定装置。   The consciousness state estimation device according to claim 3, wherein the breathing variation determination unit calculates 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. 前記呼吸ばらつき判定手段は、前記呼吸情報の時系列データに基づいて、前記呼吸のばらつき度合いとして、前記時系列データのパワースペクトル密度の周波数依存性、又は前記時系列データの自己相関性を判定する請求項3記載の意識状態推定装置。   The breathing variation determination means determines 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 based on the time series data of the breathing information. The consciousness state estimation apparatus according to claim 3. 前記推定対象者の体動及び発話の少なくとも一方を検出する体動発話検出手段を更に含み、
前記意識状態推定手段は、前記相関性判定手段によって判定された相関性と、前記体動発話検出手段による前記体動及び発話の少なくとも一方の検出結果とに基づいて、前記意識低下状態を推定する請求項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.
JP2011148101A 2011-07-04 2011-07-04 Consciousness state estimation device and program Expired - Fee Related JP5929020B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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&#39;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