JP2001061797A - Fall-asleep judgment device - Google Patents
Fall-asleep judgment deviceInfo
- Publication number
- JP2001061797A JP2001061797A JP23904499A JP23904499A JP2001061797A JP 2001061797 A JP2001061797 A JP 2001061797A JP 23904499 A JP23904499 A JP 23904499A JP 23904499 A JP23904499 A JP 23904499A JP 2001061797 A JP2001061797 A JP 2001061797A
- Authority
- JP
- Japan
- Prior art keywords
- time
- value
- biological information
- information value
- threshold
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は心拍数や脈拍数のよ
うに比較的容易に得られる生体の活動情報に基づいて生
体の入眠状態の変化を検出する入眠状態判定装置に関す
るものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a sleep state determining apparatus for detecting a change in the state of sleep of a living body based on information on the activity of the living body relatively easily obtained, such as a heart rate or a pulse rate.
【0002】[0002]
【従来の技術】人体の精神的なリラックス度合いや睡眠
の状態変化を検出すれば、その状態に応じて入眠を促進
したり、より深い睡眠状態に誘うように適宜の刺激を与
えることができる。一方、人間の状態変化を検出する装
置として、脳波、眼球運動、筋電などを含む睡眠ポリグ
ラフがあるが、これは装置が大がかりであるために研究
室や病院などの計測設備を備えた場所でしか利用でき
ず、健康機器のように日常的に使用する用途には不向き
であり、睡眠ポリグラフに代わる手軽な手段によって睡
眠の状態変化を精度良く検出することが望まれている。2. Description of the Related Art By detecting the degree of mental relaxation of a human body and a change in sleep state, it is possible to provide appropriate stimulation to promote sleep onset or induce a deeper sleep state according to the state. On the other hand, as a device that detects changes in human condition, there is a polysomnogram that includes brain waves, eye movements, myoelectricity, and the like.This is a large-scale device that is used in laboratories, hospitals, and other places equipped with measurement equipment. However, it is not suitable for daily use such as health equipment, and it is desired to accurately detect a change in sleep state by a simple alternative to polysomnography.
【0003】ここにおいて、一般に心拍信号はその人間
の状態、睡眠・覚醒、安静・運動、姿勢の変化、摂食、
さらには情動や精神的なリラックスやストレス等によっ
ても絶えず変化しており、覚醒時においてはリラックス
度合いが増大するに従って単位時間あたりの心拍数や脈
拍数は減少し、入眠するとさらに大きく減少することが
知られている。[0003] Here, generally, the heartbeat signal is a human state, sleep / wake, rest / exercise, posture change, eating,
Furthermore, it is constantly changing due to emotion, mental relaxation, stress, etc., and when awake, the heart rate and pulse rate per unit time decrease as the degree of relaxation increases, and further decrease when falling asleep. Are known.
【0004】このような知見に基づいて、心拍数や脈拍
数を測定して入眠時期を検出することが提案されてお
り、特開昭63−82673号公報には就床以降の脈拍
数の増減を指標とすることが示されており、特開昭63
−150047号公報には脈波レベルの積分値を脈拍数
で除算することにより脈波1個あたりの積分値を算出し
てこの積分値の増減を指標とすることが示されており、
さらに特開平3−41926号公報には就床時の脈拍数
または呼吸数を安静覚醒時における生体情報値とみなし
て基準値とし、生体情報値がその基準値に基づいて設定
した所定の閾値以下になった時刻を入眠時刻とすること
が示されている。[0004] Based on such knowledge, it has been proposed to detect the sleep onset by measuring the heart rate and the pulse rate. Japanese Patent Application Laid-Open No. 63-82673 discloses an increase or decrease in the pulse rate after going to bed. Is indicated as an index.
Japanese Patent Publication No. -150047 discloses that an integral value per pulse wave is calculated by dividing an integral value of a pulse wave level by a pulse rate, and an increase or decrease of the integral value is used as an index.
Further, Japanese Patent Application Laid-Open No. 3-41926 discloses that a pulse rate or a respiratory rate at bedtime is regarded as a biological information value at the time of rest and awakening and is used as a reference value, and the biological information value is equal to or less than a predetermined threshold set based on the reference value. It is shown that the time when the sleep time is reached is taken to fall asleep.
【0005】[0005]
【発明が解決しようとする課題】しかし、脈拍数の減少
度合いを指標として入眠時期を判別している上記従来構
成は、実際の睡眠の状態変化との一致率に個人差がかな
りあり、睡眠ポリグラフに比較して精度がかなり悪いと
いう問題を有していた。However, in the above-described conventional configuration in which the sleep onset is determined using the degree of decrease in the pulse rate as an index, there is a considerable individual difference in the coincidence rate with the actual change in sleep state, and the sleep polygraph However, there is a problem that the accuracy is considerably lower than that of
【0006】本発明は上記問題点の解決を目的とするも
ので、心拍数または脈拍数の変化をもとに入眠時刻を精
度良く検出できるようにした入眠判定装置を提供しよう
とするものである。SUMMARY OF THE INVENTION An object of the present invention is to solve the above-mentioned problem, and an object of the present invention is to provide a sleep onset judging device which can accurately detect a sleep onset time based on a change in a heart rate or a pulse rate. .
【0007】[0007]
【課題を解決するための手段】しかして本発明の請求項
1の発明は、単位時間あたりの心拍数または脈拍数を計
測して各時刻の生体情報値とする計測部と、計測初期に
おける生体情報値に基づいて安静覚醒時における生体情
報値とみなせる基準値を設定する基準値設定部と、上記
基準値に基づいて比較的小さな第1閾値と比較的大きな
第2閾値とを設定する閾値設定部とを備えているととも
に、、生体情報値が上記第1閾値未満になった時刻もし
くは上記第2閾値未満の値を予め設定しておいた回数を
連続でとった時刻の内、早い時刻を入眠時刻と推定する
入眠時刻推定部とを具備していることに特徴を有してい
る。生体情報値の減少傾向に時間的変動を織り込むこと
から、入眠を精度良く検出することができるものであ
る。According to a first aspect of the present invention, there is provided a measuring unit for measuring a heart rate or a pulse rate per unit time and obtaining a biological information value at each time, A reference value setting unit that sets a reference value that can be regarded as a biological information value at the time of rest and awakening based on the information value, and a threshold setting that sets a relatively small first threshold and a relatively large second threshold based on the reference value And a time when the biological information value is less than the first threshold value or a time when the number of times that the biological information value is less than the second threshold value is set in advance is set to an earlier time. It is characterized by having a sleep onset time estimating unit for estimating sleep onset time. Since the temporal variation is incorporated into the decreasing tendency of the biological information value, it is possible to accurately detect falling asleep.
【0008】また、請求項2の発明は、単位時間あたり
の心拍数または脈拍数を計測して各時刻の生体情報値と
する計測部と、生体情報値の時系列の時間変化の傾向を
表すように設定したトレンド曲線を算出するトレンド曲
線算出部と、上記各測定時刻での生体情報値の増加度合
いを算出する増分算出部と、上記各測定時刻の前の一定
時間内での上記生体情報値のばらつき度合いを算出する
ばらつき度算出部と、生体情報値の増分とばらつき度の
組み合わせが予め設定してあった閾値以下になった時刻
を入眠時刻と推定する入眠時刻推定部とを具備している
ことに特徴を有している。生体情報値の安定性を織り込
んだものとするために、入眠を精度良く検出することが
できるものである。According to a second aspect of the present invention, there is provided a measuring unit which measures a heart rate or a pulse rate per unit time and obtains a biological information value at each time, and represents a time-series trend of the biological information value. A trend curve calculation unit that calculates a trend curve set as described above, an increment calculation unit that calculates the degree of increase in the biological information value at each of the measurement times, and the biological information within a certain time before each of the measurement times. A sleep time estimating unit for estimating a time when a combination of the increment of the biological information value and the degree of the fluctuation becomes equal to or less than a preset threshold as a sleep time. The feature is that it is. In order to incorporate the stability of the biological information value, it is possible to accurately detect falling asleep.
【0009】さらに請求項3の発明は、単位時間あたり
の心拍数または脈拍数を計測して各時刻の生体情報値と
する計測部と、計測初期における生体情報値に基づいて
安静覚醒時における生体情報値とみなせる基準値を設定
する基準値設定部と、上記基準値に基づいて比較的小さ
な第1閾値と比較的大きな第2閾値とを設定する閾値設
定部と、生体情報値の時系列の平均的な時間変化の傾向
を表すように設定したトレンド曲線を算出するトレンド
曲線算出部と、上記各測定時刻での生体情報値の増分を
算出する増分算出部と、上記各測定時刻の前の一定時間
内での上記生体情報値のばらつき度合いを算出するばら
つき度算出部と、生体情報値が上記第1閾値未満になっ
た時刻もしくは上記第2閾値未満の値を予め設定してお
いた回数を連続でとった時刻もしくは生体情報値の増分
とばらつき度の組み合わせが予め設定してあった閾値以
下になった時刻の内、早い時刻を入眠時刻と推定する入
眠時刻推定部とを具備していることに特徴を有してお
り、生体情報値の減少傾向に時間的変動を織り込むとと
もに生体情報値の安定性も織り込むことから、入眠をさ
らに精度良く検出することができる。Further, the invention according to claim 3 is a measuring section which measures a heart rate or a pulse rate per unit time and obtains a biological information value at each time, and a living body at rest awakening based on the biological information value at the initial stage of measurement. A reference value setting unit that sets a reference value that can be regarded as an information value; a threshold setting unit that sets a relatively small first threshold value and a relatively large second threshold value based on the reference value; A trend curve calculation unit that calculates a trend curve set to represent the tendency of an average time change, an increment calculation unit that calculates an increase in the biological information value at each of the measurement times, A variability calculator for calculating the variability of the biometric information value within a certain period of time, and the number of times that the biometric information value has become less than the first threshold or a value less than the second threshold set in advance Continuously The sleep time estimating unit for estimating the earliest time as the sleep time during the time when the combination of the measured time or the increment of the biological information value and the degree of variation become equal to or less than the preset threshold value. Since it has a characteristic and incorporates temporal fluctuation into the decreasing tendency of the biological information value and also incorporates stability of the biological information value, falling asleep can be detected with higher accuracy.
【0010】計測初期における生体情報値を基に算出し
た基準値について、その後の生体情報値が基準値を上回
る値であった場合に基準値を書き換える基準値再設定部
を合わせて具備しているものとするのも好ましい。計測
の初期状態に拘わらず入眠を精度良く検出することがで
きる。With respect to the reference value calculated based on the biological information value at the initial stage of measurement, a reference value resetting unit for rewriting the reference value when the subsequent biological information value exceeds the reference value is provided. It is also preferable to use It is possible to accurately detect falling asleep regardless of the initial state of measurement.
【0011】[0011]
【発明の実施の形態】以下本発明を実施の形態の一例に
基づいて詳述すると、図1は実施の形態の一例の基本構
成を示しており、計測部1と基準値設定部2、閾値設定
部3、そして入眠時刻推定部4で構成されている。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, the present invention will be described in detail with reference to an example of an embodiment. FIG. 1 shows a basic configuration of an example of the embodiment, wherein a measuring unit 1, a reference value setting unit 2, a threshold value It comprises a setting unit 3 and a sleep onset time estimation unit 4.
【0012】計測部1はたとえば心拍センサからの心拍
信号が入力されるもので、波形整形を行うことでパルス
状の信号を得た後、単位時間毎のパルス数を計数して、
計数値を単位時間(たとえば1分)あたりの心拍数H
(t)として出力する。The measuring section 1 receives, for example, a heart rate signal from a heart rate sensor. After obtaining a pulse-like signal by performing waveform shaping, the number of pulses per unit time is counted.
Heart rate H per unit time (for example, 1 minute)
Output as (t).
【0013】計測部1から出力される心拍数H(t)は、
安静覚醒時における生体情報値とみなせる基準値を設定
する基準値設定部2と入眠時刻推定部4に入力される。The heart rate H (t) output from the measuring unit 1 is
The values are input to a reference value setting unit 2 and a sleep onset time estimation unit 4 for setting a reference value that can be regarded as a biological information value when the user is awake.
【0014】基準値設定部2は、安静覚醒時における心
拍数とみなせる基準心拍数Hrを基準値として算出する
もので、図2に示すように、計測部1から出力された心
拍数H(t)の値の最初の3個を平均して、その結果を基
準値Hrとして採用する。ただし、前回心拍数H(t-1)
より今回心拍数H(t)が5以上大きい時は体動などによ
る雑音が入ったものとして、今回心拍数H(t)の値とし
て前回心拍数H(t-1)の値と同じものを用いて基準値H
rを算出している。The reference value setting unit 2 calculates, as a reference value, a reference heart rate Hr that can be regarded as a heart rate at the time of rest and awakening. As shown in FIG. 2, the heart rate H (t) output from the measurement unit 1 is calculated. ) Are averaged, and the result is adopted as the reference value Hr. However, the previous heart rate H (t-1)
If the heart rate H (t) is greater than 5 this time, it is assumed that noise due to body motion etc. has entered, and the value of the current heart rate H (t) is the same as the value of the previous heart rate H (t-1). Using the reference value H
r is calculated.
【0015】基準値設定部2から出力された基準値Hr
は閾値設定部3に入力される。閾値設定部3は、入力さ
れた基準値Hrに対して80%〜95%の値を第1閾値
ならびに第2閾値として設定するもので、この時、第1
閾値より第2閾値のほうが大きな値をとるものとする
(たとえば、93%を第1閾値,95%を第2閾値)。Reference value Hr output from reference value setting unit 2
Is input to the threshold setting unit 3. The threshold setting unit 3 sets a value of 80% to 95% with respect to the input reference value Hr as the first threshold and the second threshold.
It is assumed that the second threshold value is larger than the threshold value (for example, 93% is the first threshold value and 95% is the second threshold value).
【0016】入眠の判定は入眠時刻推定部4において、
計測部1から1分ごとに入力されるパルス数と閾値設定
部3から入力される2つの閾値とを用いて行われる。す
なわち、入眠時刻推定部4は、図3に示すように、入力
された心拍数H(t)を第1閾値と比較して、第1閾値よ
り小さければ入眠と推定する。また、第1閾値よりも大
きければ心拍数H(t)を第2閾値と比較し、第2閾値よ
り小さければ変数Hsを1増やす。次に、変数Hsが予
め設定しておいた数(例えば4)以上かどうかを判定
し、予め設定しておいた数以上であれば入眠と判定す
る。予め設定しておいた数未満であれば次に入力される
心拍数H(t+1)をもとに同じ解析を繰り返し、心拍数H
(x)が第1閾値より小さくなるか、変数Hsが予め設定
しておいた数以上になった時点で入眠と判断する。The sleep onset determination is performed by the sleep onset time estimating unit 4.
This is performed using the number of pulses input from the measuring unit 1 every minute and two thresholds input from the threshold setting unit 3. That is, as shown in FIG. 3, the falling-on time estimating unit 4 compares the input heart rate H (t) with the first threshold value, and if it is smaller than the first threshold value, estimates that it is falling asleep. If the value is larger than the first threshold value, the heart rate H (t) is compared with the second threshold value. If the value is smaller than the second threshold value, the variable Hs is increased by one. Next, it is determined whether or not the variable Hs is equal to or greater than a preset number (for example, 4). If the variable Hs is equal to or greater than the preset number, it is determined that the user falls asleep. If the number is less than the preset number, the same analysis is repeated based on the next input heart rate H (t + 1), and the heart rate H
When (x) becomes smaller than the first threshold value or when the variable Hs becomes equal to or more than a preset number, it is determined that the user falls asleep.
【0017】図4に他例を示す。計測部1から出力され
た心拍数H(t)は、生体情報値の時系列の変化の傾向を
表すように設定したトレンド曲線を算出するトレンド曲
線算出部5と、各測定時刻での生体情報値の増加度合い
を算出する増分算出部6と、各測定時刻の前の一定時間
内での生体情報値のばらつき度合いを算出するばらつき
度算出部7に入力される。FIG. 4 shows another example. The heart rate H (t) output from the measuring unit 1 is calculated by using a trend curve calculating unit 5 that calculates a trend curve set so as to indicate the tendency of the time series change of the biological information value, and the biological information at each measurement time. The values are input to an increment calculator 6 that calculates the degree of increase in the value and a degree-of-variation calculator 7 that calculates the degree of fluctuation of the biological information value within a certain period of time before each measurement time.
【0018】トレンド曲線算出部5では、生体情報値の
時系列の平均的な時間変化の傾向を表すように設定した
トレンド曲線を算出する。すなわち、トレンド曲線算出
部5は図5に示すように入力された心拍数H(t)からト
レンド曲線FT(t)を次式に基づいて求める。The trend curve calculating section 5 calculates a trend curve set so as to represent the average time-varying trend of the biological information value. That is, the trend curve calculation unit 5 calculates the trend curve FT (t) from the input heart rate H (t) as shown in FIG.
【0019】 t=1の時 FT(1)=H(1) t>1で且つ H(t)≧FT(t-1)の時 FT(t)=FT(t-1) t>1で且つ H(t)<FT(t-1)の時 FT(t)=H(t) このようにして求めたトレンド曲線FT(t)の一例を図
6に示す。図中、細線イが心拍数を、太線ロがトレンド
曲線FT(t)を示している。When t = 1 FT (1) = H (1) When t> 1 and when H (t) ≧ FT (t−1), FT (t) = FT (t−1) When t> 1 And when H (t) <FT (t-1) FT (t) = H (t) An example of the trend curve FT (t) thus obtained is shown in FIG. In the figure, the thin line A indicates the heart rate, and the thick line B indicates the trend curve FT (t).
【0020】トレンド曲線算出部5から出力されたトレ
ンド曲線FT(t)は各測定時刻での生体情報値の増加度
合いを算出する増分算出部6に入力される。増分算出部
6では、トレンド曲線FT(t)からの心拍数H(t)の増加
度合いを表す指標である増分ZOU(t)(ZOU(t)=H(t)−
FT(t))を算出する。増分ZOU(t)の一例を図7に示
す。The trend curve FT (t) output from the trend curve calculation section 5 is input to an increment calculation section 6 for calculating the degree of increase of the biological information value at each measurement time. In the increment calculation unit 6, an increment ZOU (t) (ZOU (t) = H (t)-), which is an index indicating the degree of increase of the heart rate H (t) from the trend curve FT (t).
FT (t)) is calculated. FIG. 7 shows an example of the increment ZOU (t).
【0021】ばらつき度算出部7は計測部1より入力さ
れた心拍数H(t)のばらつき度合いを算出する。まず、
心拍数H(t)と1分前の心拍数H(t-1)の差の絶対値(差
分SABUN(t)とする)を算出する。なお、差分SABUN(t)は |(H(t)−H(t-1)|<6の場合 SABUN(t)=|(H(t)−H(t-1)| |(H(t)−H(t-1)|≧6の場合 SABUN(t)=6 とする。The variability calculator 7 calculates the variability of the heart rate H (t) input from the measuring unit 1. First,
The absolute value of the difference between the heart rate H (t) and the heart rate H (t-1) one minute ago (referred to as a difference SABUN (t)) is calculated. The difference SABUN (t) is | (H (t) −H (t−1) | <6. SABUN (t) = | (H (t) −H (t−1) || (H (t ) −H (t−1) | ≧ 6 SABUN (t) = 6.
【0022】さらに、差分SABUN(t)を用いてばらつき度
BARA(t)を t=1の場合 BARA(1)=SABUN(1) t=2の場合 BARA(2)=SABUN(1)+SABUN(2) t=3の場合 BARA(3)=SABUN(1)+SABUN(2)+SABUN(3) t=4の場合 BARA(4)=SABUN(1)+SABUN(2)+SABUN(3)+SABUN(4) t≧5の場合 BARA(t)=SABUN(t-4)+SABUN(t-3)+SABUN(t-2)+SABUN
(t-1)+SABUN(t) として求める。このようにして求めたばらつき度BARA
(t)の一例を図8に示す。Further, using the difference SABUN (t),
When BARA (t) is t = 1 BARA (1) = SABUN (1) When t = 2 BARA (2) = SABUN (1) + SABUN (2) When t = 3 BARA (3) = SABUN (1 ) + SABUN (2) + SABUN (3) When t = 4 BARA (4) = SABUN (1) + SABUN (2) + SABUN (3) + SABUN (4) When t ≧ 5 BARA (t) = SABUN (t-4 ) + SABUN (t-3) + SABUN (t-2) + SABUN
Calculate as (t-1) + SABUN (t). Variation degree BARA obtained in this way
FIG. 8 shows an example of (t).
【0023】入眠時刻推定部4による入眠の推定は、増
分算出部6とばらつき度算出部7から1分毎に入力され
てくる増分ZOU(t)とばらつき度BARA(t)及び予め設定し
ている閾値との比較によって行う。すなわち、入眠時刻
推定部5では、増分ZOU(t)とばらつき度BARA(t)の両方
が予め設定してある閾値(たとえば、増分ZOU(t)が3で
ばらつき度BARA(t)が7)より小さくなれば入眠と推定す
る。The sleep onset estimation by the sleep onset time estimating section 4 is performed by setting the increment ZOU (t) and the variation BARA (t) inputted every minute from the increment calculating section 6 and the variation calculating section 7 and presetting them. This is performed by comparing with a threshold value. That is, in the sleep onset time estimation unit 5, both the increment ZOU (t) and the variation BARA (t) are preset thresholds (for example, the increment ZOU (t) is 3 and the variation BARA (t) is 7). If it becomes smaller, it is estimated that it falls asleep.
【0024】図9に更に他例を示す。図中の計測部1、
基準値設定部2、閾値設定部3、トレンド曲線算出部
5、増分算出部6、ばらつき度算出部7は上記の2例に
おけるものと同じである。FIG. 9 shows still another example. The measurement unit 1 in the figure,
The reference value setting section 2, threshold value setting section 3, trend curve calculation section 5, increment calculation section 6, and variation degree calculation section 7 are the same as those in the above two examples.
【0025】入眠の推定は、入眠時刻推定部4におい
て、計測部1から1分毎に入力される心拍数H(t)と、
閾値設定部3から入力される第1閾値及び第2閾値とを
用いてまずなされる。すなわち入眠時刻推定部4では図
10に示すように入力された心拍数H(t)を第1閾値と
比較し、第1閾値より小さければ入眠と推定する。ま
た、第1閾値よりも大きければ心拍数H(t)を第2閾値
と比較し、第2閾値より小さければ変数Hsを1増やし
て変数Hsが予め設定しておいた数(例えば4)以上かど
うかを判定し、予め設定しておいた数以上であれば入眠
と判定する。予め設定しておいた数未満であれば次に入
力される心拍数H(t+1)をもとに同じ解析を繰り返し、
心拍数H(x)が第1閾値より小さくなるか、変数Hsが
予め設定しておいた数以上になるか、増分ZOU(t)とばら
つき度BARA(t)の両方が予め設定されていた閾値(たとえ
ば、増分ZOU(t)が3でばらつき度(BARA(t)が7)より小
さくなれば入眠と推定する。The sleep onset estimation is performed by a sleep onset time estimating unit 4 in which the heart rate H (t) inputted from the measuring unit 1 every minute is represented by:
This is first performed using the first threshold and the second threshold input from the threshold setting unit 3. That is, the sleep onset time estimation unit 4 compares the input heart rate H (t) with the first threshold as shown in FIG. If the heart rate H (t) is larger than the first threshold value, the heart rate H (t) is compared with the second threshold value. If the heart rate H (t) is smaller than the second threshold value, the variable Hs is increased by 1 and the variable Hs is equal to or more than a preset number (for example, 4). It is determined whether or not it falls asleep if the number is equal to or greater than a preset number. If it is less than the preset number, repeat the same analysis based on the next input heart rate H (t + 1),
Either the heart rate H (x) is smaller than the first threshold value, the variable Hs is equal to or greater than a preset number, or both the increment ZOU (t) and the variation BARA (t) are preset. If the threshold value (for example, the increment ZOU (t) is 3 and the degree of variation (BARA (t)) is smaller than 7), it is estimated that the user falls asleep.
【0026】図11に別の例を示す。これは基準値設定
部2で求めた基準値Hrを、より適切な値に補正する基
準値再設定部8を設けたものを示しており、計測部1か
ら出力される心拍数H(t)は、安静覚醒時における生体
情報値とみなせる基準値を設定する基準値設定部2と、
生体の初期状態によっては基準値を再設定する基準値再
設定部8に入力する。基準値設定部2では、前述のよう
に心拍数H(t)の最初の3個の値を平均して基準値Hr
を設定する。FIG. 11 shows another example. This shows that a reference value resetting unit 8 for correcting the reference value Hr obtained by the reference value setting unit 2 to a more appropriate value is provided, and the heart rate H (t) output from the measuring unit 1 is shown. Is a reference value setting unit 2 that sets a reference value that can be regarded as a biological information value at the time of rest and awakening,
The reference value is input to the reference value resetting unit 8 for resetting the reference value depending on the initial state of the living body. As described above, the reference value setting unit 2 averages the first three values of the heart rate H (t) to obtain the reference value Hr.
Set.
【0027】一方、基準値再設定部8は、図12に示す
ように、計測部1より出力された心拍数H(t)からの連
続3個の値を平均した値Hpと基準値設定部2において
設定された基準値Hrとを比較し、新平均値Hpのほう
が基準値Hrより大きければ、新平均値Hpの値を基準
値Hrに代入する。新たな基準値Hrは閾値設定部4に
入力され、新たな第1閾値および第2閾値が設定され
て、これらの値に基づいて入眠時刻推定部4において入
眠の判定が行われる。この場合、計測開始時から3分間
の心拍数の平均に比べて、その後の心拍数が増加した場
合にも精度良く入眠時刻を推定することができる。On the other hand, as shown in FIG. 12, a reference value resetting unit 8 calculates a value Hp obtained by averaging three consecutive values from the heart rate H (t) output from the measuring unit 1 and a reference value setting unit. Then, the new average value Hp is compared with the reference value Hr, and if the new average value Hp is larger than the reference value Hr, the value of the new average value Hp is substituted for the reference value Hr. The new reference value Hr is input to the threshold setting unit 4, a new first threshold and a new second threshold are set, and the sleep onset time estimation unit 4 determines sleep on the basis of these values. In this case, it is possible to accurately estimate the sleep onset time even when the heart rate increases after that as compared to the average of the heart rate for three minutes from the start of the measurement.
【0028】以上の各例では、生体の活動状況の指標と
なる生体情報値として心拍数を用いたが、脈拍数でも同
様の処理を行うことができるのはもちろんである。In each of the above examples, the heart rate is used as the biological information value as an index of the activity state of the living body. However, the same processing can be performed with the pulse rate.
【0029】[0029]
【発明の効果】以上のように請求項1の発明において
は、計測初期における生体情報値に基づいて安静覚醒時
における生体情報値とみなせる基準値を設定する基準値
設定部と、上記基準値に基づいて比較的小さな第1閾値
と比較的大きな第2閾値を設定する閾値設定部とを備え
るとともに、生体情報値が上記第1閾値未満になった時
刻もしくは上記第2閾値未満の値を予め設定しておいた
回数を連続でとった時刻の内、早い時刻を入眠時刻と推
定する入眠時刻推定部を具備しているものであり、生体
情報値の減少傾向に時間的変動を織り込んで入眠を判定
することから、入眠を精度良く検出することができる。As described above, according to the first aspect of the present invention, a reference value setting unit for setting a reference value that can be regarded as a biological information value at the time of rest and awakening based on the biological information value at the initial stage of measurement; A threshold setting unit that sets a relatively small first threshold and a relatively large second threshold based on the threshold value, and presets a time at which the biological information value falls below the first threshold or a value less than the second threshold. It is provided with a sleep onset time estimating unit for estimating an earlier time as a sleep onset time among the times when the number of times set in succession is taken. From the determination, falling asleep can be detected with high accuracy.
【0030】また、請求項2の発明では、生体情報値の
時系列の平均的な時間変化の傾向を表すように設定した
トレンド曲線を算出するトレンド曲線算出部と、上記各
測定時刻での生体情報値の増分を算出する増分算出部
と、上記各測定時刻の前の一定時間内での上記生体情報
値のばらつき度合いを算出するばらつき度算出部と、生
体情報値の増分とばらつき度の組み合わせが予め設定し
てあった閾値以下になった時刻を入眠時刻と推定する入
眠時刻推定部を具備しているものであり、この場合、生
体情報値の安定性が織り込まれていることになり、従っ
て、入眠を精度良く検出することができる。Further, according to the invention of claim 2, a trend curve calculating section for calculating a trend curve set so as to represent a tendency of an average time change of the time series of the biological information value; An increment calculator for calculating the increment of the information value, a variability calculator for calculating the variance of the biological information value within a certain time period before each of the measurement times, and a combination of the increment and the variability of the biological information value Is provided with a sleep onset time estimating unit for estimating the time at which the time becomes equal to or less than a preset threshold as a sleep onset time, in this case, the stability of the biological information value is incorporated, Therefore, falling asleep can be accurately detected.
【0031】また、請求項3の発明では、生体情報値の
減少傾向と時間的変動に安定性も織り込むことから、入
眠を精度良く検出することができる。According to the third aspect of the present invention, since the stability is incorporated into the decreasing tendency and the temporal variation of the biological information value, it is possible to accurately detect falling asleep.
【0032】そして請求項4の発明によれば、計測初期
における生体情報値を基に算出した基準値を決定した後
で、生体情報値が基準値を上回る値であらた場合に基準
値を書き換えるようにしているために、生体の初期状態
に拘わらず入眠を精度良く検出することができる。According to the fourth aspect of the present invention, after the reference value calculated based on the biological information value at the initial stage of the measurement is determined, the reference value is rewritten when the biological information value exceeds the reference value. Therefore, falling asleep can be accurately detected regardless of the initial state of the living body.
【図面の簡単な説明】[Brief description of the drawings]
【図1】本発明の実施の形態の一例のブロック回路図で
ある。FIG. 1 is a block circuit diagram of an example of an embodiment of the present invention.
【図2】同上の基準値設定部の処理内容を示すフローチ
ャートである。FIG. 2 is a flowchart showing processing contents of a reference value setting unit according to the first embodiment;
【図3】同上の入眠時刻推定部の処理内容を示すフロー
チャートである。FIG. 3 is a flowchart showing processing contents of a sleep onset time estimating unit according to the embodiment.
【図4】他例のブロック回路図である。FIG. 4 is a block circuit diagram of another example.
【図5】同上のトレンド曲線算出部の処理内容を示すフ
ローチャートである。FIG. 5 is a flowchart showing processing contents of a trend curve calculation unit according to the embodiment.
【図6】同上のトレンド曲線の一例を示すグラフであ
る。FIG. 6 is a graph showing an example of the above trend curve.
【図7】同上の増分の一例を示すグラフである。FIG. 7 is a graph showing an example of the above increment.
【図8】同上のばらつき度の一例を示すグラフである。FIG. 8 is a graph showing an example of a variation degree according to the embodiment.
【図9】更に他例のブロック回路図である。FIG. 9 is a block circuit diagram of still another example.
【図10】同上の入眠時刻推定部の処理内容を示すフロ
ーチャートである。FIG. 10 is a flowchart showing a processing content of a sleep onset time estimating unit according to the embodiment.
【図11】別の例のブロック回路図である。FIG. 11 is a block circuit diagram of another example.
【図12】同上の基準値再設定部の処理内容を示すフロ
ーチャートである。FIG. 12 is a flowchart showing a processing content of a reference value resetting unit according to the embodiment.
1 計測部 2 基準値設定部 3 閾値設定部 4 入眠時刻推定部 DESCRIPTION OF SYMBOLS 1 Measurement part 2 Reference value setting part 3 Threshold setting part 4 Sleeping time estimation part
───────────────────────────────────────────────────── フロントページの続き (72)発明者 道盛 章弘 大阪府門真市大字門真1048番地松下電工株 式会社内 (72)発明者 中野 紀夫 大阪府門真市大字門真1048番地松下電工株 式会社内 (72)発明者 萩原 啓 大阪府門真市大字門真1048番地松下電工株 式会社内 Fターム(参考) 4C017 AA02 AA10 4C038 PP05 PQ00 PS00 ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Akihiro Michimori 1048 Kadoma Kadoma, Osaka Prefecture Matsushita Electric Works, Ltd. (72) Inventor Kei Hagiwara 1048 Kadoma Kadoma, Osaka Prefecture F-term in Matsushita Electric Works Co., Ltd. (reference) 4C017 AA02 AA10 4C038 PP05 PQ00 PS00
Claims (4)
計測して各時刻の生体情報値とする計測部と、計測初期
における生体情報値に基づいて安静覚醒時における生体
情報値とみなせる基準値を設定する基準値設定部と、上
記基準値に基づいて比較的小さな第1閾値と比較的大き
な第2閾値とを設定する閾値設定部とを備えるととも
に、生体情報値が上記第1閾値未満になった時刻もしく
は上記第2閾値未満の値を予め設定しておいた回数を連
続で示した時刻の内、早い時刻を入眠時刻と推定する入
眠時刻推定部を具備していることを特徴とする入眠判定
装置。1. A measuring unit that measures a heart rate or a pulse rate per unit time and obtains a biological information value at each time, and a reference value that can be regarded as a biological information value at the time of rest and awakening based on the biological information value at the beginning of measurement. And a threshold value setting unit that sets a relatively small first threshold value and a relatively large second threshold value based on the reference value, and the biological information value is less than the first threshold value. A sleep onset time estimating unit for estimating an earlier time as a sleep onset time among times in which the number of consecutive times or a value less than the second threshold value is preset. Sleep detection device.
計測して各時刻の生体情報値とする計測部と、生体情報
値の時系列の時間変化の傾向を表すように設定したトレ
ンド曲線を算出するトレンド曲線算出部と、上記各測定
時刻での生体情報値の増加度合いを算出する増分算出部
と、上記各測定時刻の前の一定時間内での上記生体情報
値のばらつき度合いを算出するばらつき度算出部と、生
体情報値の増分とばらつき度の組み合わせが予め設定し
てあった閾値以下になった時刻を入眠時刻と推定する入
眠時刻推定部とを具備していることを特徴とする入眠判
定装置。2. A measuring unit that measures a heart rate or a pulse rate per unit time and obtains a biological information value at each time, and a trend curve set so as to represent a tendency of a time series of the biological information value over time. A trend curve calculating unit for calculating, an increment calculating unit for calculating an increasing degree of the biological information value at each of the measurement times, and calculating a variation degree of the biological information value within a fixed time before each of the measuring times. A degree-of-variation calculating section, and a sleep-onset time estimating section for estimating a sleep-onset time at a time when a combination of the increment of the biometric information value and the degree of variance becomes equal to or less than a preset threshold. Sleep detection device.
計測して各時刻の生体情報値とする計測部と、計測初期
における生体情報値に基づいて安静覚醒時における生体
情報値とみなせる基準値を設定する基準値設定部と、上
記基準値に基づいて比較的小さな第1閾値と比較的大き
な第2閾値とを設定する閾値設定部と、生体情報値の時
系列の平均的な時間変化の傾向を表すように設定したト
レンド曲線を算出するトレンド曲線算出部と、上記各測
定時刻での生体情報値の増分を算出する増分算出部と、
上記各測定時刻の前の一定時間内での上記生体情報値の
ばらつき度合いを算出するばらつき度算出部と、生体情
報値が上記第1閾値未満になった時刻もしくは上記第2
閾値未満の値を予め設定しておいた回数を連続でとった
時刻もしくは生体情報値の増分とばらつき度の組み合わ
せが予め設定してあった閾値以下になった時刻の内、早
い時刻を入眠時刻と推定する入眠時刻推定部とを具備し
ていることを特徴とする入眠判定装置。3. A measuring unit that measures a heart rate or a pulse rate per unit time and obtains a biological information value at each time, and a reference value that can be regarded as a biological information value at the time of rest and awakening based on the biological information value at the initial stage of measurement. A threshold value setting unit that sets a relatively small first threshold value and a relatively large second threshold value based on the reference value, and a time-series average time change of the biological information value. A trend curve calculation unit that calculates a trend curve set to represent the trend, and an increment calculation unit that calculates an increase in the biological information value at each of the measurement times,
A variation calculation unit that calculates the variation of the biological information value within a certain time period before each of the measurement times; and a time when the biological information value falls below the first threshold or the second
Of the times when the number of times less than the threshold value is set in advance, or the time when the combination of the increment of the biometric information value and the degree of variation becomes equal to or less than the preset threshold value, the earliest time is the sleep time. A sleep onset determining unit for estimating sleep onset.
した基準値について、その後の生体情報値が基準値を上
回る値であった場合に基準値を書き換える基準値再設定
部を合わせて具備していることを特徴とする請求項1ま
たは3記載の入眠判定装置。4. A reference value resetting unit for rewriting a reference value calculated based on a biological information value at an initial stage of measurement when the subsequent biological information value exceeds the reference value. The sleep onset judging device according to claim 1 or 3, wherein:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP23904499A JP3829544B2 (en) | 1999-08-26 | 1999-08-26 | Sleep detection device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP23904499A JP3829544B2 (en) | 1999-08-26 | 1999-08-26 | Sleep detection device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001061797A true JP2001061797A (en) | 2001-03-13 |
| JP3829544B2 JP3829544B2 (en) | 2006-10-04 |
Family
ID=17039050
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP23904499A Expired - Fee Related JP3829544B2 (en) | 1999-08-26 | 1999-08-26 | Sleep detection device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3829544B2 (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003052649A (en) * | 2001-08-10 | 2003-02-25 | Kawasaki Heavy Ind Ltd | Care receiver monitoring device |
| JP2003310584A (en) * | 2002-04-26 | 2003-11-05 | Yoshida Dental Mfg Co Ltd | Dental care instrument and dental consulting method |
| CN106913313A (en) * | 2015-12-28 | 2017-07-04 | 中国移动通信集团公司 | A kind of sleep monitor method and system |
| JP2018019879A (en) * | 2016-08-03 | 2018-02-08 | ヤマハ株式会社 | Sleep detection device, sleep detection method and program |
| CN109893091A (en) * | 2017-12-07 | 2019-06-18 | 八乐梦医用床有限公司 | Sleep waking state determination method, device, equipment and its storage medium |
| CN110101392A (en) * | 2018-02-01 | 2019-08-09 | 德国福维克控股公司 | System with vital sign data sensor and analytical unit |
| JP2021003381A (en) * | 2019-06-26 | 2021-01-14 | パラマウントベッド株式会社 | Sleep evaluation apparatus and determination device |
| JP2021159571A (en) * | 2020-04-02 | 2021-10-11 | 株式会社ミラクルプランニング | Sleep analysis method and equipment |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5104815B2 (en) * | 2009-05-19 | 2012-12-19 | ダイキン工業株式会社 | Heart rate detector |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0341926A (en) * | 1989-07-07 | 1991-02-22 | Matsushita Electric Works Ltd | Detector for change in sleeping state and sleeping state controller |
| JPH03151938A (en) * | 1989-11-07 | 1991-06-28 | Nippondenso Co Ltd | Awakening degree detector |
| JPH0349686Y2 (en) * | 1988-02-18 | 1991-10-23 | ||
| JPH05184558A (en) * | 1992-01-17 | 1993-07-27 | Mitsubishi Motors Corp | Driver abnormality detection device |
| JP2001061820A (en) * | 1999-08-26 | 2001-03-13 | Matsushita Electric Works Ltd | Sleeping state judging method and sleeping state judging device |
-
1999
- 1999-08-26 JP JP23904499A patent/JP3829544B2/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0349686Y2 (en) * | 1988-02-18 | 1991-10-23 | ||
| JPH0341926A (en) * | 1989-07-07 | 1991-02-22 | Matsushita Electric Works Ltd | Detector for change in sleeping state and sleeping state controller |
| JPH03151938A (en) * | 1989-11-07 | 1991-06-28 | Nippondenso Co Ltd | Awakening degree detector |
| JPH05184558A (en) * | 1992-01-17 | 1993-07-27 | Mitsubishi Motors Corp | Driver abnormality detection device |
| JP2001061820A (en) * | 1999-08-26 | 2001-03-13 | Matsushita Electric Works Ltd | Sleeping state judging method and sleeping state judging device |
Non-Patent Citations (1)
| Title |
|---|
| 井辺浩行、道盛章弘: "「心拍数の変動と睡眠状態の関係について その3」", ヒューマン・インターフェース・シンポジウム論文集, vol. 13, JPN4006003837, 1997, JP, pages 189 - 192, ISSN: 0000717156 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003052649A (en) * | 2001-08-10 | 2003-02-25 | Kawasaki Heavy Ind Ltd | Care receiver monitoring device |
| JP2003310584A (en) * | 2002-04-26 | 2003-11-05 | Yoshida Dental Mfg Co Ltd | Dental care instrument and dental consulting method |
| CN106913313A (en) * | 2015-12-28 | 2017-07-04 | 中国移动通信集团公司 | A kind of sleep monitor method and system |
| JP2018019879A (en) * | 2016-08-03 | 2018-02-08 | ヤマハ株式会社 | Sleep detection device, sleep detection method and program |
| CN109893091B (en) * | 2017-12-07 | 2022-08-19 | 八乐梦医用床有限公司 | Sleep-awake state determination method, apparatus, device, and storage medium therefor |
| JP2019098068A (en) * | 2017-12-07 | 2019-06-24 | パラマウントベッド株式会社 | Sleep state determination device and program |
| CN109893091A (en) * | 2017-12-07 | 2019-06-18 | 八乐梦医用床有限公司 | Sleep waking state determination method, device, equipment and its storage medium |
| US11647945B2 (en) | 2017-12-07 | 2023-05-16 | Paramount Bed Co., Ltd. | Evaluating apparatus and non-transitory computer readable medium storing program |
| CN110101392A (en) * | 2018-02-01 | 2019-08-09 | 德国福维克控股公司 | System with vital sign data sensor and analytical unit |
| CN110101392B (en) * | 2018-02-01 | 2022-12-13 | 德国福维克控股公司 | System with vital sign data sensor and analysis unit |
| JP2021003381A (en) * | 2019-06-26 | 2021-01-14 | パラマウントベッド株式会社 | Sleep evaluation apparatus and determination device |
| JP7198165B2 (en) | 2019-06-26 | 2022-12-28 | パラマウントベッド株式会社 | Sleep evaluation device and determination device |
| JP2021159571A (en) * | 2020-04-02 | 2021-10-11 | 株式会社ミラクルプランニング | Sleep analysis method and equipment |
| JP7432233B2 (en) | 2020-04-02 | 2024-02-16 | 株式会社人間と科学の研究所 | Sleep analysis method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3829544B2 (en) | 2006-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3788126B2 (en) | Sleep state determination device | |
| EP2914171B1 (en) | Measuring psychological stress from cardiovascular and activity signals | |
| CN100464703C (en) | Sleep state estimation device and method thereof | |
| JP3733133B2 (en) | Sleep state estimation device | |
| JP4582642B2 (en) | Sleep stage determination device | |
| US7351206B2 (en) | Apparatus for and method of biotic sleep state determining | |
| EP2120701B1 (en) | Computer program product, device and method for estimating a change in level of arousal of a user | |
| JPH0341926A (en) | Detector for change in sleeping state and sleeping state controller | |
| WO2005082252A1 (en) | Method for judging stage of sleep | |
| CN101095612A (en) | Device and method for monitoring biological information | |
| JP2016107095A (en) | Sleep stage determination device, sleep stage determination method, sleep stage determination program, and stress degree determination method | |
| JP2013013542A (en) | Apparatus and program for estimating state of consciousness | |
| JP2004089267A (en) | Sleep depth estimation device and futon having the same | |
| JP3829544B2 (en) | Sleep detection device | |
| JPH07143972A (en) | Method and apparatus for judging sleeping condition | |
| JP2950038B2 (en) | Sleep evaluation device | |
| US20180146915A1 (en) | Sleep stage determination apparatus and sleep stage determination method | |
| JP2007000280A (en) | Arousal level drop determination device | |
| JP2006271474A (en) | Sleep state estimation device | |
| JP2718303B2 (en) | Sleep state determination device | |
| JP2022159825A (en) | Biological state determination device, biological state determination method, and biological state determination program | |
| US20240164678A1 (en) | Fatigue Evaluation Based On Sleep Quality And Physical Activity | |
| JPH04272744A (en) | sleep detection device | |
| KR20220143599A (en) | Method and detecting device for determining tiredness of user | |
| CN116235009A (en) | Temperature estimation device, air-conditioning control device, air-conditioning control system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051019 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051025 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051226 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060228 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060501 |
|
| 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: 20060620 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060703 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090721 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090721 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090721 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100721 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100721 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110721 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120721 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120721 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130721 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130721 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140721 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |