JPH10262941A - Electroencephalogram processor - Google Patents
Electroencephalogram processorInfo
- Publication number
- JPH10262941A JPH10262941A JP9069753A JP6975397A JPH10262941A JP H10262941 A JPH10262941 A JP H10262941A JP 9069753 A JP9069753 A JP 9069753A JP 6975397 A JP6975397 A JP 6975397A JP H10262941 A JPH10262941 A JP H10262941A
- Authority
- JP
- Japan
- Prior art keywords
- waveform
- electroencephalogram
- pulse wave
- unit
- pulse
- 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
- 238000011109 contamination Methods 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 claims description 58
- 238000007781 pre-processing Methods 0.000 claims description 45
- 238000000034 method Methods 0.000 claims description 35
- 238000009499 grossing Methods 0.000 claims description 34
- 210000004556 brain Anatomy 0.000 claims description 24
- 210000000624 ear auricle Anatomy 0.000 claims description 7
- 238000002156 mixing Methods 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 25
- 230000006996 mental state Effects 0.000 description 9
- 239000012528 membrane Substances 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 241000204801 Muraenidae Species 0.000 description 1
- 229910021607 Silver chloride Inorganic materials 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 235000013351 cheese Nutrition 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- HKZLPVFGJNLROG-UHFFFAOYSA-M silver monochloride Chemical compound [Cl-].[Ag+] HKZLPVFGJNLROG-UHFFFAOYSA-M 0.000 description 1
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、脳波処理装置に関
し、特に心電波形による雑音が混在している脳波の波形
から心電波形による雑音を除去するための脳波処理装置
に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an electroencephalogram processing apparatus, and more particularly to an electroencephalogram processing apparatus for removing noise due to an electrocardiographic waveform from a waveform of an electroencephalogram in which noise due to an electrocardiographic waveform is mixed.
【0002】[0002]
【従来の技術】図9は従来の脳波処理装置の一例を示す
ブロック図である。2. Description of the Related Art FIG. 9 is a block diagram showing an example of a conventional electroencephalogram processing apparatus.
【0003】従来の脳波処理装置としては、特公平6−
7824号公報に開示されている脳波処理装置が知られ
ている。A conventional brain wave processing apparatus is disclosed in
An electroencephalogram processing device disclosed in Japanese Patent No. 7824 is known.
【0004】上記脳波処理装置は、図9に示すように、
2個の脳波計測用電極91aおよび91bと、脳波計測
用電極91aおよび91bからの信号を入力してそれら
の電位差を増幅する脳電位増幅器93と、脳電位増幅器
93からのアナログ信号をデジタル信号に変換するアナ
ログデジタル変換器94と、アナログデジタル変換器9
4の出力信号を格納する一時記憶部95と、心電波形を
検出するための2個の心電図用電極92aおよび92b
と、心電図用電極92aおよび92bからの信号を入力
してそれらの電位差を増幅する心電図用増幅器96と、
心電図用増幅器96の出力信号からトリガー信号を発生
するトリガー信号発生部99と、トリガー信号発生部9
9からの信号によって一時記憶部95に格納した脳波波
形に雑音として混入している心電波形を推定して差引
き、雑音を除去した脳波波形を抽出する演算部98とを
備えている。[0004] As shown in FIG.
Two brain wave measurement electrodes 91a and 91b, a brain potential amplifier 93 that receives signals from the brain wave measurement electrodes 91a and 91b and amplifies a potential difference therebetween, and converts an analog signal from the brain potential amplifier 93 into a digital signal. An analog-to-digital converter 94 for converting and an analog-to-digital converter 9
Temporary storage unit 95 for storing the output signal of ECG 4 and two electrocardiogram electrodes 92a and 92b for detecting an electrocardiographic waveform.
And an ECG amplifier 96 for receiving signals from the ECG electrodes 92a and 92b and amplifying a potential difference therebetween,
A trigger signal generator 99 for generating a trigger signal from an output signal of the ECG amplifier 96; and a trigger signal generator 9
And an arithmetic unit 98 for estimating and subtracting an electrocardiographic waveform mixed as noise into the electroencephalogram waveform stored in the temporary storage unit 95 based on the signal from the temporary storage unit 95 and extracting a noise-free electroencephalogram waveform.
【0005】上述のように構成した脳波処理装置は、脳
波計測用電極91aおよび91bと心電図用電極92a
および92bとによって脳波および心電波形を同時に検
出し、脳波は脳電位増幅器93によって増幅してアナロ
グデジタル変換器94によってアナログ信号をデジタル
信号に変換し、一時記憶部95に格納する。一方、心電
波形は、心電図用増幅器96によって電位差を増幅さ
れ、トリガー信号発生部99によってトリガー信号を発
生させて一時記憶部95に送られる。トリガー信号とし
ては、例えばR波のような特徴的な波形が現われたとき
の時点に関する信号を発生して一時記憶部95に送る。
演算部98は、トリガー信号発生部99からのトリガー
信号によって所定の時間だけ遡及した時点を起点とし
て、脳波の波形を個々のセグメントに分割し、各セグメ
ントを順次加算して平均値を求めることによって脳波波
形のなかに混入している心電波形の推定波形を求め、起
点を同期させながらこの推定波形を心電波形が混在して
いる元の脳波波形から差し引くことにより、心電波形を
除去した脳波波形を抽出している。The electroencephalogram processing apparatus constructed as described above comprises electroencephalogram measurement electrodes 91a and 91b and an electrocardiogram electrode 92a.
And 92b simultaneously detect an electroencephalogram and an electrocardiographic waveform. The electroencephalogram is amplified by a brain potential amplifier 93, converted from an analog signal to a digital signal by an analog-to-digital converter 94, and stored in a temporary storage unit 95. On the other hand, the potential difference of the electrocardiogram waveform is amplified by an electrocardiogram amplifier 96, a trigger signal is generated by a trigger signal generation unit 99, and the trigger signal is sent to a temporary storage unit 95. As the trigger signal, for example, a signal relating to a point in time when a characteristic waveform such as an R-wave appears is generated and sent to the temporary storage unit 95.
The arithmetic unit 98 divides the waveform of the brain wave into individual segments, starting from a point in time when the trigger signal from the trigger signal generating unit 99 traces back for a predetermined time, and sequentially adds each segment to obtain an average value. The estimated waveform of the electrocardiogram waveform mixed in the electroencephalogram waveform was obtained, and the electrocardiogram waveform was removed by subtracting this estimated waveform from the original electroencephalogram waveform containing the electrocardiogram waveform while synchronizing the starting point. EEG waveforms are extracted.
【0006】[0006]
【発明が解決しようとする課題】上述したような従来の
脳波処理装置は、脳波の検出以外に、心電波形も同時に
検出しなければならないため、心電波形の検出のための
電極の装着が必要であり、このため検査準備が面倒であ
り、また、被験者が制限されて適用範囲が限定されると
いう問題点を有している。Since the conventional electroencephalogram processing device as described above must simultaneously detect an electrocardiographic waveform in addition to detecting an electroencephalogram, it is necessary to mount electrodes for detecting the electrocardiographic waveform. It is necessary, and therefore, the preparation for the examination is troublesome, and there is a problem that the subject is limited and the applicable range is limited.
【0007】本発明の目的は、上述のような従来の脳波
処理装置の問題点を解消することにより、心電波形検出
用の電極の装着よりもセンサの装着が簡便な脈派を利用
し、心電波形の特徴的な波形に対応する脈派の特徴的な
波形の出現時点を検出することにより、心電波形の推定
波形を求めるための脳波波形のセグメント分割を容易に
し、脳波波形から心電波形を除去するための計測時間を
短縮して脳波処理を容易にすることができる脳波処理装
置を提供することにある。[0007] An object of the present invention is to solve the above-mentioned problems of the conventional electroencephalogram processing apparatus, and to utilize a pulse train in which a sensor can be mounted more easily than an electrode for detecting an electrocardiographic waveform. By detecting the time when the characteristic waveform of the pulse corresponding to the characteristic waveform of the electrocardiographic waveform is detected, the segmentation of the electroencephalogram waveform for obtaining the estimated waveform of the electrocardiographic waveform is facilitated. An object of the present invention is to provide an electroencephalogram processing device capable of shortening a measurement time for removing a radio wave shape and facilitating electroencephalogram processing.
【0008】[0008]
【課題を解決するための手段】本発明の第一の脳波処理
装置は、2極間の電位差としての脳波波形を検出するた
めの活性電極および基準電極を有する脳波センサと、耳
朶または指先から脈派波形を検出する脈派センサと、前
記脳波センサの前記活性電極と前記基準電極との電位差
および前記脈派センサの電位とを増幅する増幅器と、前
記増幅器の出力のアナログ信号をデジタル信号に変換す
るアナログデジタル変換器と、前記アナログデジタル変
換器の出力信号を格納する記憶部と、前記記憶部に格納
した信号から周期性を有する脈派波形の極値およびその
極値の出現する時点を検出する脈波極値検出部と、前記
脈波極値検出部の出力信号によって所定の時間だけ遡及
した時点を前記脳波波形の分割の起点に設定する脳波分
割起点設定部と、前記脳波分割起点設定部によって設定
された起点から前記脳波波形を個々のセグメントに分割
し前記セグメントを順次加算して平均値を求めることに
よって前記脳波波形に雑音として混入している心電波形
の推定波形を求める雑音混入推定部と、前記雑音混入推
定部で得られた前記推定波形を前記起点を同期させて前
記脳波波形から差し引くことによって雑音の混入を除去
した脳波信号を抽出する脳波信号抽出部とを備えるもの
であり、特に、前記脈波極値検出部として、前記記憶部
に格納した脈派波形を表示装置の画面に表示し、前記画
面から周期性を有する脈派波形の極値およびその極値の
出現する時点を直接検出する波形表示部を有するか、ま
たは、前記記憶部に格納した脈派波形のスムージング処
理を行う脈派波形前処理部と、前記脈派波形前処理部で
スムージング処理された信号から各時点における1階差
分値を算出して前記1階差分値が充分に小さな値の所定
の値に等しいかまたは前記所定の値よりも小さい場合に
その値を極値としてその極値の出現する時点を検出する
差分型脈波極値自動検出部とを有するか、または、前記
記憶部に格納した脈派波形のスムージング処理を行う脈
派波形前処理部と、前記脈派波形前処理部でスムージン
グ処理された信号に対してウエーブレット処理を施すこ
とによって極値を求めて前記極値の出現する時点を検出
するウエーブレット変換型脈波極値自動検出部とを有す
るか、または、前記記憶部に格納した脈派波形のスムー
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号から各時点におけ
る2階差分値を算出して前記2階差分値が充分に小さな
値の所定の値に等しいかまたは前記所定の値よりも小さ
い場合を変曲点として前記変曲点の出現する時点を検出
する脈波変曲点自動検出部とを有する前記脈波極値検出
部を備えるものである。A first electroencephalogram processing device according to the present invention comprises an electroencephalogram sensor having an active electrode and a reference electrode for detecting an electroencephalogram waveform as a potential difference between two electrodes, and a pulse from an earlobe or a fingertip. A pulse sensor for detecting a pulse waveform, an amplifier for amplifying a potential difference between the active electrode and the reference electrode of the electroencephalogram sensor and a potential of the pulse sensor, and converting an analog signal output from the amplifier into a digital signal. An analog-to-digital converter, a storage unit for storing an output signal of the analog-to-digital converter, and detecting, from the signal stored in the storage unit, an extreme value of a pulse waveform having periodicity and a time point at which the extreme value appears. A pulse wave extremum detecting unit, and an electroencephalogram division starting point setting unit that sets a point in time when the output signal of the pulse wave extremal value detecting unit goes back for a predetermined time as a starting point of the division of the electroencephalogram waveform, Estimating an electrocardiographic waveform mixed as noise into the electroencephalogram waveform by dividing the electroencephalogram waveform into individual segments from the start point set by the electroencephalogram division start point setting unit and sequentially adding the segments to obtain an average value A noise mixing estimating unit for obtaining a waveform, and an electroencephalogram signal extracting unit for extracting an electroencephalogram signal from which noise has been mixed by subtracting the estimated waveform obtained by the noise mixing estimating unit from the brain wave waveform by synchronizing the starting point. In particular, as the pulse wave extremum detection unit, a pulse wave waveform stored in the storage unit is displayed on a screen of a display device, and an extreme value of a pulse wave waveform having periodicity from the screen. Either has a waveform display unit that directly detects the time point when the extreme value appears, or a pulse waveform preprocessing unit that performs smoothing processing of the pulse waveform stored in the storage unit, Calculating a first-order difference value at each point in time from the signal subjected to the smoothing processing in the wave-form preprocessing unit, and when the first-order difference value is equal to a sufficiently small value or smaller than the predetermined value; Or a difference type pulse wave extremum automatic detection unit for detecting a time point when the extremum appears as an extremum, or before performing a smoothing process of a pulse waveform stored in the storage unit. A processing unit, a wavelet transform type pulse wave extremum for obtaining an extremum by performing a wavelet process on the signal subjected to the smoothing process in the pulse wave waveform preprocessing unit and detecting a time point at which the extremum appears. A pulse wave pre-processing unit having an automatic detection unit, or performing a smoothing process of a pulse waveform stored in the storage unit, and at each time from a signal smoothed by the pulse wave pre-processing unit. Calculating the second-order difference value and detecting the time point at which the inflection point appears as an inflection point when the second-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value. And a pulse wave inflection point automatic detection unit.
【0009】本発明の第二の脳波処理装置は、2極間の
電位差としての脳波波形を検出するための活性電極およ
び基準電極を有する脳波センサと、耳朶または指先から
膜派波形を検出する脈派センサと、前記脳波センサの前
記活性電極と前記基準電極との電位差および前記脈派セ
ンサの電位とを増幅する増幅器と、前記増幅器の出力の
アナログ信号をデジタル信号に変換するアナログデジタ
ル変換器と、前記アナログデジタル変換器の出力信号を
格納する記憶部と、前記記憶部に格納した信号から周期
性を有する膜派波形の極値およびその極値の出現する時
点を検出する脈波極値検出部と、前記脈波極値検出部に
おいて検出された隣接する極値間の差として脈波振幅値
を算出する脈波振幅値算出部と、前記脈波極値検出部に
おいて検出された前記極値のなかから前記脈波振幅値算
出部によって算出された前記脈波振幅値が所定の範囲内
にある極値とその出現時点とを選出する加算セグメント
判断部とを有する判断部と、前記判断部からの出力信号
によって所定の時間だけ遡及した時点を前記脳波波形の
分割の起点に設定する脳波分割起点設定部と、前記脳波
分割起点設定部によって設定された起点から前記脳波波
形を個々のセグメントに分割し前記セグメントを順次加
算して平均値を求めることによって前記脳波波形に雑音
として混入している心電波形の推定波形を求める雑音混
入推定部と、前記雑音混入推定部で得られた前記推定波
形を前記起点を同期させて前記脳波波形から差し引くこ
とによって雑音の混入を除去した脳波信号を抽出する脳
波信号抽出部とを備えるものであり、特に、前記脈波極
値検出部として、前記記憶部に格納した脈派波形を表示
装置の画面に表示し、前記画面から周期性を有する脈派
波形の極値およびその極値の出現する時点を直接検出す
る波形表示部を有するか、または、前記記憶部に格納し
た脈派波形のスムージング処理を行う脈派波形前処理部
と、前記脈派波形前処理部でスムージング処理された信
号から各時点における1階差分値を算出して前記1階差
分値が充分に小さな値の所定の値に等しいかまたは前記
所定の値よりも小さい場合にその値を極値としてその極
値の出現する時点を検出する差分型脈波極値自動検出部
とを有するか、または、前記記憶部に格納した脈派波形
のスムージング処理を行う脈派波形前処理部と、前記脈
派波形前処理部でスムージング処理された信号に対して
ウエーブレット処理を施すことによって極値を求めて前
記極値の出現する時点を検出するウエーブレット変換型
脈波極値自動検出部とを有するか、または、前記記憶部
に格納した脈派波形のスムージング処理を行う脈派波形
前処理部と、前記脈派波形前処理部でスムージング処理
された信号から各時点における2階差分値を算出して前
記2階差分値が充分に小さな値の所定の値に等しいかま
たは前記所定の値よりも小さい場合を変曲点として前記
変曲点の出現する時点を検出する脈波変曲点自動検出部
とを有する前記脈波極値検出部を備えるものである。A second electroencephalogram processing apparatus according to the present invention comprises an electroencephalogram sensor having an active electrode and a reference electrode for detecting an electroencephalogram waveform as a potential difference between two electrodes, and a pulse for detecting a membrane wave from an earlobe or a fingertip. A group sensor, an amplifier that amplifies the potential difference between the active electrode and the reference electrode of the brain wave sensor and the potential of the pulse sensor, and an analog-to-digital converter that converts an analog signal output from the amplifier into a digital signal. A storage unit for storing an output signal of the analog-to-digital converter, and a pulse wave extremum detection for detecting, from the signal stored in the storage unit, an extremum of a periodic membrane wave having a periodicity and a time when the extremum appears. A pulse wave amplitude value calculating unit that calculates a pulse wave amplitude value as a difference between adjacent extreme values detected by the pulse wave extreme value detecting unit, and a pulse wave amplitude value detecting unit that is detected by the pulse wave extreme value detecting unit. A determination unit having an addition segment determination unit that selects an extreme value whose pulse wave amplitude value calculated by the pulse wave amplitude value calculation unit from within the extreme value is within a predetermined range and the present time of the extreme value; An electroencephalogram division starting point setting section that sets a point in time that is retroactive by a predetermined time by the output signal from the determination section as a starting point of the electroencephalogram waveform division, and individually calculates the electroencephalogram waveform from the starting point set by the electroencephalogram division starting point setting section. And a noise mixing estimating unit for obtaining an estimated waveform of an electrocardiographic waveform mixed as noise in the brain wave waveform by sequentially adding the segments to obtain an average value. An EEG signal extraction unit for extracting an EEG signal in which noise is removed by subtracting the estimated waveform from the EEG waveform while synchronizing the starting point. In particular, as the pulse wave extremum detection unit, the pulse wave stored in the storage unit is displayed on a screen of a display device, and the extreme value of the pulse wave waveform having periodicity from the screen and the time when the extreme value appears Has a waveform display unit that directly detects, or a pulse waveform pre-processing unit that performs smoothing processing of the pulse waveform stored in the storage unit, and each of the signals that have been smoothed by the pulse waveform pre-processing unit. Calculating a first-order difference value at a point in time, and when the first-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value, the value is regarded as an extreme value and the extreme value appears Or a pulse type waveform pre-processing unit for performing a smoothing process of a pulse type waveform stored in the storage unit, and a pulse type waveform pre-processing unit for performing a smoothing process on the pulse type waveform stored in the storage unit. For the processed signal A wavelet conversion type pulse wave extremum automatic detection unit that determines the extremum by performing a wavelet process and detects the time when the extremum appears, or of the pulse waveform stored in the storage unit. A pulse waveform pre-processing unit that performs a smoothing process, and calculates a second-order difference value at each time from the signal that has been subjected to the smoothing process in the pulse waveform pre-processing unit, and the second-order difference value is a predetermined value having a sufficiently small value. A pulse wave inflection point automatic detection unit that detects a point in time when the inflection point appears as an inflection point when the pulse wave extremum value is equal to or smaller than the predetermined value. It is.
【0010】[0010]
【発明の実施の形態】次に、本発明の実施の形態につい
て図面を参照して説明する。Next, embodiments of the present invention will be described with reference to the drawings.
【0011】本発明の第一の実施形態を示すブロック図
である。図10は実測された脳波データおよび心電波形
データおよび脈波データの一例を示す波形図、図11は
図10の脈波データの一周期分を示す波形図、図12は
心理状態の影響を受ける脈波データの一例を示す波形
図、図13は心理状態の影響を受けた脈波データにおけ
る波形の変動の一例を示す波形図である。FIG. 1 is a block diagram showing a first embodiment of the present invention. FIG. 10 is a waveform diagram showing an example of actually measured brain wave data, electrocardiographic waveform data and pulse wave data, FIG. 11 is a waveform diagram showing one cycle of the pulse wave data in FIG. 10, and FIG. FIG. 13 is a waveform diagram showing an example of received pulse wave data, and FIG. 13 is a waveform diagram showing an example of a waveform change in pulse wave data affected by a mental state.
【0012】図1の実施形態は、2極間の電位差として
の脳波波形を検出するための活性電極1aおよび基準電
極1bを有する脳波センサと、指先に装着して脈波波形
(指尖脈波)を検出する脈派センサ2a、または、耳朶
に装着して脈波波形(耳朶脈波)を検出する脈派センサ
2bと、脳波センサの活性電極1aと基準電極1bとの
電位差および脈派センサ2aまたは2bの電位とを増幅
する増幅器(アンプ)3と、増幅器の出力のアナログ信
号をデジタル信号に変換するアナログデジタル変換器
(A/D)4と、アナログデジタル変換器4の出力信号
を格納する記憶部5と、記憶部5に格納した脈派波形を
表示装置の画面に表示し、その画面から周期性を有する
脈派波形の極値およびその極値の出現する時点を直接検
出する波形表示部6aを有する脈波極値検出部6と、脈
波極値検出部6の出力信号によって所定の時間だけ遡及
した時点を脳波波形の分割の起点に設定する脳波分割起
点設定部7と、脳波分割起点設定部7によって設定され
た起点から脳波波形を個々のセグメントに分割しそれら
のセグメントを順次加算して平均値を求めることによっ
て脳波波形に雑音として混入している心電波形の推定波
形を求める雑音混入推定部8と、雑音混入推定部8で得
られた推定波形を上記の起点を同期させて脳波波形から
差し引くことによって雑音の混入を除去した脳波信号を
抽出する脳波信号抽出部9とを備えている。The embodiment of FIG. 1 shows an electroencephalogram sensor having an active electrode 1a and a reference electrode 1b for detecting an electroencephalogram waveform as a potential difference between two electrodes, a pulse wave waveform (fingertip pulse wave) attached to a fingertip. ), A pulse sensor 2b attached to an earlobe and detecting a pulse waveform (earlobe pulse wave), a potential difference between the active electrode 1a and the reference electrode 1b of the electroencephalogram sensor, and a pulse sensor. An amplifier (amplifier) 3 for amplifying the potential of 2a or 2b, an analog / digital converter (A / D) 4 for converting an analog signal output from the amplifier into a digital signal, and an output signal of the analog / digital converter 4 are stored. And a waveform for displaying the pulse waveform stored in the storage unit 5 on the screen of the display device, and directly detecting the extreme value of the periodic pulse waveform and the time when the extreme value appears from the screen. Display 6 A pulse wave extremum detecting unit 6 having an EEG, an electroencephalogram division starting point setting unit 7 for setting a point in time when the output signal of the pulse wave extremal value detecting unit 6 goes back for a predetermined time as a starting point of the electroencephalogram waveform division, A noise for obtaining an estimated electrocardiographic waveform mixed as a noise in an electroencephalogram waveform by dividing the electroencephalogram waveform into individual segments from the starting point set by the setting unit 7 and sequentially adding the segments to obtain an average value. A brain wave signal extraction unit that extracts a brain wave signal from which the noise is removed by subtracting the estimated waveform obtained by the noise contamination estimation unit from the brain wave waveform by synchronizing the above-mentioned starting point; ing.
【0013】次に上述のように構成した脳波処理装置の
動作について説明する。Next, the operation of the electroencephalogram processing apparatus configured as described above will be described.
【0014】本実施形態は、図10に示した脳波データ
15と、心電データ16と、脈波データ17との3種類
のデータに関連する。これらのうちの脳波データ15
は、活性電極および基準電極(例えば Ag/AgCl
皿電極)を有する脳波センサに、増幅器(アンプ)(例
えばNECメディカルシステムズ(株)製SYNAX1
200)を接続し、心電データ16は、3個1組の電極
(胸部V誘導)にそれらの電極間の電位差を増幅するア
ンプ(例えばNECメディカルシステムズ(株)製PO
LYGRAPH360)を接続し、脈波データ17は、
脈波データを電位差として検出するための脈派センサ
(例えばNECメディカルシステムズ(株)製脈波ファ
ピックアップ)に、検出した電位差を増幅するためのア
ンプ(例えばNECメディカルシステムズ(株)製PO
LYGRAPH360)を接続して計測し、各アンプの
出力のアナログ信号を記録計(例えばNECメディカル
システムズ(株)製OMNICORDER 8M14)
に出力した波形を示す図である。これらのうちの脈波デ
ータ17は、心電データ16の影響を受けて心電データ
16と同じ周期の周期性を有している。この1周期(約
700msec)の中に、極大値17aおよび17cや
極小値17bおよび17d等のいくつかの極値を有して
いる。脈波データ17を心電データ16と比較すると、
例えば脈波データ17の極大値17aの出現時点と心電
データ16のR波16aの出現時点との間隔は一定(約
340msec)である。従って、脈波データ17の極
大の出現時点がわかれば、心電データ16を使用しなく
ても、脈波データ17の極値を利用して、その時点から
一定時間遡及した時点を起点として、脳波データ15を
個々のセグメントに分割することができる。This embodiment relates to three types of data, ie, electroencephalogram data 15, electrocardiogram data 16, and pulse wave data 17 shown in FIG. EEG data 15 of these
Are the active and reference electrodes (eg, Ag / AgCl
An amplifier (amplifier) (for example, SYNAX1 manufactured by NEC Medical Systems Co., Ltd.)
200), and an electrocardiographic data 16 is applied to a set of three electrodes (chest V lead) using an amplifier (for example, a PO medical device manufactured by NEC Medical Systems Co., Ltd.) that amplifies the potential difference between the electrodes.
LYGRAPH360), and the pulse wave data 17 is
A pulse sensor for detecting pulse wave data as a potential difference (for example, a pulse wave pick-up manufactured by NEC Medical Systems Co., Ltd.) and an amplifier for amplifying the detected potential difference (for example, PO manufactured by NEC Medical Systems Co., Ltd.)
LYGRAPH360) is connected and measured, and the analog signal of the output of each amplifier is recorded by a recorder (for example, OMNICORDER 8M14 manufactured by NEC Medical Systems Co., Ltd.).
FIG. 7 is a diagram showing a waveform output to FIG. The pulse wave data 17 among these has the same periodicity as the electrocardiogram data 16 under the influence of the electrocardiogram data 16. In this one cycle (about 700 msec), there are several extreme values such as the maximum values 17a and 17c and the minimum values 17b and 17d. When the pulse wave data 17 is compared with the electrocardiogram data 16,
For example, the interval between the present time of the peak value 17a of the pulse wave data 17 and the present time of the R wave 16a of the electrocardiographic data 16 is constant (about 340 msec). Therefore, if the present point of the maximum of the pulse wave data 17 is known, the extremum of the pulse wave data 17 is used, and the time point retroactive from the time point by a certain time is used as a starting point, without using the electrocardiogram data 16. The electroencephalogram data 15 can be divided into individual segments.
【0015】本実施形態は、活性電極1aおよび基準電
極1bを有する脳波センサによって脳波データ15を計
測し、脈派センサ2aまたは脈派センサ2bによって脈
波データ17を計測し、脳波センサの活性電極1aと基
準電極1bとの電位差および脈派センサ2aまたは2b
の電位とを増幅器(アンプ)3によって増幅し、増幅器
3の出力のアナログ信号をアナログデジタル変換器(A
/D)4によってデジタル信号に変換し、、アナログデ
ジタル変換器4の出力信号をパーソナルコンピュータの
ハードディスク等の記憶部5に格納する。In the present embodiment, the electroencephalogram data 15 is measured by an electroencephalogram sensor having an active electrode 1a and a reference electrode 1b, and the pulse wave data 17 is measured by a pulse sensor 2a or a pulse sensor 2b. Potential difference between 1a and reference electrode 1b and pulse sensor 2a or 2b
Is amplified by an amplifier (amplifier) 3 and an analog signal output from the amplifier 3 is converted into an analog-to-digital converter (A).
/ D) 4 to convert the signal into a digital signal, and store the output signal of the analog-to-digital converter 4 in a storage unit 5 such as a hard disk of a personal computer.
【0016】脈波極値検出部6は、記憶部5に格納した
脳波データおよび脈派データをアナログ波形として波形
表示部6aの画面に表示し、ユーザは、その画面から例
えば周期性を有する脈派データの極大値17aの出現す
る時点を読み取る。脈派データの極値の出現時点の読み
取りは、上記のように、脳波データおよび脈派データの
波形を記録計に出力し、その出力記録用紙から読み取る
ことも可能である。The pulse wave extremum detection unit 6 displays the brain wave data and the pulse data stored in the storage unit 5 as analog waveforms on the screen of the waveform display unit 6a. The time when the maximum value 17a of the group data appears is read. The reading of the extremum of the pulse data at the present time can be performed by outputting the waveforms of the electroencephalogram data and the pulse data to the recorder as described above, and reading from the output recording paper.
【0017】脳波分割起点設定部7は、脈波極値検出部
6の出力信号によって設定された所定の時間だけ遡及し
た時点を、脳波データ15の波形を分割する起点に設定
する。具体的には、以下のようなデータ処理を行う。An electroencephalogram division starting point setting section 7 sets a point in time when the waveform of the electroencephalogram data 15 is divided by a point in time retroactive by a predetermined time set by an output signal of the pulse wave extreme value detecting section 6. Specifically, the following data processing is performed.
【0018】すなわち、検出された極値出現時点を ti *(i=1,2,……,L)(Lは検出された極値出
現時点の個数) とすると、このときの脳波データ15の波形分割の起点
は、 ti=ti *+Tretro(i=1,2,……,L) (Tretroは遡及すべき時間)となる。That is, assuming that the detected extreme value present time is t i * (i = 1, 2,..., L) (L is the number of detected extreme value present times), the electroencephalogram data 15 The starting point of the waveform division is t i = t i * + T retro (i = 1, 2,..., L) (T retro is the time to be traced back ).
【0019】図10の例において、脈波データ17の極
大値17aの出現時点を利用する場合は、 Tretro=TA ー R+TR (TA ー R は心電データ16のR波16aの出現時点と脈
波データ17の極大値17aの出現時点との差、TR は
心電データ16のR波16aの出現時点から遡及すべき
時間)となる。TR の具体的な数値としては、200m
secに設定することが可能である。[0019] In the example of FIG. 10, when using the occurrence time of the maximum value 17a of the pulse wave data 17, the T retro = T A over R + T R (T A over R is R wave 16a of the electrocardiograph data 16 the difference between the appearance time of the maximum value 17a occurrence time and the pulse wave data 17, T R is the time to be retroactively from the appearance point of the R-wave 16a of the electrocardiograph data 16). Specific numerical values of T R, 200 meters
sec.
【0020】雑音混入推定部8は、脳波分割起点設定部
7によって設定された起点ti (i=1,2,……,
L)から脳波データ17を個々のセグメントに分割し、
それらのセグメントを順次加算して平均値を求めること
によって脳波波形に雑音として混入している心電データ
の推定波形を求める。具体的には、以下のようなデータ
処理を行う。The noise contamination estimating section 8 has a starting point t i (i = 1, 2,...) Set by the brain wave dividing starting point setting section 7.
L), the EEG data 17 is divided into individual segments,
The estimated waveform of electrocardiogram data mixed as noise into the electroencephalogram waveform is obtained by sequentially adding these segments to obtain an average value. Specifically, the following data processing is performed.
【0021】すなわち、時点tにおける脳波のデジタル
データの値を y(t) とすると、起点ti (i=
1,2,……,L)に基づいた脳波データの個々のセグ
メントは、 yj *(τ+tj )=y(τ+tj )[u(τ)−u(τ
−(tj+1 −tj ))] として求めることができる。That is, assuming that the value of the digital data of the electroencephalogram at time t is y (t), the starting point t i (i =
Each segment of the electroencephalogram data based on 1, 2,..., L) is represented by y j * (τ + t j ) = y (τ + t j ) [u (τ) −u (τ
− (T j + 1 −t j ))].
【0022】但し、u(τ)はステップ関数であり、
[u(τ)−u(τ−(tj+1 −tj ))]は、矩形波
タイプの窓関数となる。更に、 Where u (τ) is a step function,
[U (τ) −u (τ− (t j + 1 −t j ))] is a rectangular wave type window function. Furthermore,
【0023】により、時点τにおける混入雑音の推定値
を求めることができる。但し、上式においては、各起点
からの時点τにおける推定値であるので、元の脳波デー
タの時点tに戻すためには、 As a result, an estimated value of the mixed noise at the time point τ can be obtained. However, in the above equation, since it is an estimated value at time τ from each starting point, in order to return to the time t of the original brain wave data,
【0024】という演算を行う。Is performed.
【0025】脳波信号抽出部9は、雑音混入推定部8で
得られた推定波形を、雑音が混入している元の脳波デー
タから差し引いて、脳波信号を抽出する。すなわち、 The electroencephalogram signal extraction unit 9 extracts an electroencephalogram signal by subtracting the estimated waveform obtained by the noise contamination estimation unit 8 from the original electroencephalogram data containing noise. That is,
【0026】によって脳波信号を抽出する。Thus, an electroencephalogram signal is extracted.
【0027】図2は本発明の第二の実施形態を示すブロ
ック図である。FIG. 2 is a block diagram showing a second embodiment of the present invention.
【0028】図2の実施形態は、図1の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した膜派波形の
スムージング処理を行う脈派波形前処理部36aと、脈
派波形前処理部36aでスムージング処理された信号か
ら各時点における1階差分値を算出してその1階差分値
が充分に小さな値の所定の値に等しいかまたは所定の値
よりも小さい場合にその値を極値としてその極値の出現
する時点を検出する差分型脈波極値自動検出部36bと
を有する脈波極値検出部36を設けたものである。その
他の構成要素は、図1の実施形態と同じである。In the embodiment of FIG. 2, instead of the pulse wave extremum detecting unit 6 of the embodiment of FIG. 1, a pulse wave pre-processing unit 36a for performing a smoothing process of the membrane wave stored in the storage unit 5; When the first-order difference value at each time is calculated from the signal subjected to the smoothing processing by the pulse waveform pre-processing unit 36a, and the first-order difference value is equal to or smaller than a sufficiently small predetermined value. And a pulse-wave extremal-value detection unit 36 having a differential-type pulse-wave extremum-automatic detection unit 36b that detects the time when the extremum appears as an extreme value. Other components are the same as those in the embodiment of FIG.
【0029】次に上述のように構成した脳波処理装置の
動作について、第一の実施形態と同じ動作を行う部分の
動作説明は省略し、第一の実施形態と異なる部分の動作
について説明する。Next, regarding the operation of the electroencephalogram processing device configured as described above, the operation of the portion performing the same operation as the first embodiment will be omitted, and the operation of the portion different from the first embodiment will be described.
【0030】脈派波形前処理部36aは、記憶部5に格
納した脈派データのスムージング処理を行う。スムージ
ング処理には、種々の方法があるが、例えば、脈派デー
タに低域通過型のデジタルフィルタを適用したり、移動
平均法を適用したりすることによってスムージング処理
を行うことが可能である。The pulse waveform preprocessing unit 36a performs smoothing processing of pulse data stored in the storage unit 5. There are various methods for the smoothing process. For example, the smoothing process can be performed by applying a low-pass digital filter to the pulse data or applying a moving average method.
【0031】差分型脈波極値自動検出部36bは、脈派
波形前処理部36aでスムージング処理された信号から
各時点における1階差分値を算出し、その1階差分値が
充分に小さな値の所定の値に等しいかまたは所定の値よ
りも小さい場合にその値を極値(極大値または極小値)
としてその極値の出現する時点を検出する。具体的に
は、以下のようなデータ処理を行う。The differential type pulse wave extremum automatic detection unit 36b calculates the first-order difference value at each time point from the signal smoothed by the pulse waveform pre-processing unit 36a, and the first-order difference value is a sufficiently small value. If the value is less than or equal to the specified value, the value is taken to the extreme value (maximum value or minimum value)
At the time when the extreme value appears. Specifically, the following data processing is performed.
【0032】すなわち、スムージング処理された脈派の
デジタルデータを xi とし、そのi番目の時点におけ
る1階差分値を That is, the digital data of the pulse train subjected to the smoothing process is defined as x i, and the first-order difference value at the i-th time point is defined as x i.
【0033】と定義する。このとき xi が極大値であ
ることは、 The definition is as follows. At this time, the fact that x i is a local maximum value means that
【0034】によって、 xi が極小値であることは、 By the way, the fact that x i is a minimum value means that
【0035】によって、それぞれ認識することが可能で
ある。Thus, each can be recognized.
【0036】例えば、図11の脈波データ17の1周期
分の波形17nに適用すると、極大値17a・極小値1
7b・極大値17c・極小値17dを順次に検出するこ
とができる。For example, when applied to the waveform 17n of one cycle of the pulse wave data 17 in FIG. 11, the maximum value 17a and the minimum value 1
7b, the maximum value 17c, and the minimum value 17d can be sequentially detected.
【0037】図3は本発明の第三の実施形態を示すブロ
ック図である。FIG. 3 is a block diagram showing a third embodiment of the present invention.
【0038】図3の実施形態は、図1の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した脈派波形の
スムージング処理を行う脈派波形前処理部56aと、脈
派波形前処理部56aでスムージング処理された信号に
対してウエーブレット処理を施すことによって極値を求
めてその極値の出現する時点を検出するヴァヴェレット
変換型脈波極値自動検出部56bとを有する脈波極値検
出部56を設けたものである。その他の構成要素は、図
1の実施形態と同じである。In the embodiment shown in FIG. 3, instead of the pulse wave extremum detection unit 6 of the embodiment shown in FIG. 1, a pulse waveform pre-processing unit 56a for performing a smoothing process of the pulse waveform stored in the storage unit 5, A wavelet transform type pulse wave extremum automatic detection unit 56b which obtains an extremum by performing a wavelet process on the signal smoothed by the pulse waveform pre-processing unit 56a and detects a time when the extremum appears A pulse wave extremum detection unit 56 having the following is provided. Other components are the same as those in the embodiment of FIG.
【0039】次に上述のように構成した脳波処理装置の
動作について、第一の実施形態と異なる部分の動作につ
いて説明する。Next, as for the operation of the electroencephalogram processing device configured as described above, the operation of a portion different from that of the first embodiment will be described.
【0040】脈派波形前処理部56aは、図2は本発明
の第二の実施形態と同様に、記憶部5に格納した脈派波
形のスムージング処理を行う。The pulse waveform pre-processing unit 56a performs smoothing processing of the pulse waveform stored in the storage unit 5 in FIG. 2, as in the second embodiment of the present invention.
【0041】ウエーブレット(Wavelet)変換型
脈波極値自動検出部56bは、脈派波形前処理部56a
でスムージング処理された信号に対してウエーブレット
(Wavelet)変換処理を施すことにより、一層明
確化された極値と、その極値の出現する時点を検出す
る。具体的には、以下のようなデータ処理を行う。The wavelet conversion type pulse wave extremum automatic detection unit 56b includes a pulse wave pre-processing unit 56a.
By performing a wavelet conversion process on the signal subjected to the smoothing process in step (1), a further clarified extremum and a time point at which the extremum appears are detected. Specifically, the following data processing is performed.
【0042】任意の実数値関数f(t)のWavele
t変換は、一般に、 Wavelet of arbitrary real-valued function f (t)
The t-transform is generally
【0043】但し、ψはマザーウエーブレット(Mot
her Wavelet)と呼ばれ、時空間的に局在す
る「波」を表わす種々の関数の総称でり、種々の関数と
しては、ハールウエーブレット(haar Wavel
et)、ガボーアウエーブレット(Gabor Wav
elet)、マルヴアーウエーブレット(Malver
Wavelet)、モーレイウエーブレット(Morl
ey Wavelet)、メキシカンウエーブレット
(Mexican Wavelet)、フレンチウエー
ブレット(French Wavelet)、シャノン
ウエーブレット(Shannon Wavelet)、
メイヤーウエーブレット(Meyer Wavele
t)、ドーブチーズウエーブレット(Daubechi
es Wavelet)、シムレットウエーブレット
(Symlet Wavelet)、コアフレットウエ
ーブレット(Coiflet Wavelet)、スプ
ラインウエーブレット(Spline Wavele
t)等が知られている。また、aおよびbは、任意の実
数であり、2つの整数jおよびkにょって、 However, ψ indicates a mother wavelet (Mot
Her Wavelet is a general term for various functions representing “waves” that are localized in space and time, and various functions include a Haar Wavelet.
et), Gabor Wave
elet), Malver Wavelet (Malver)
Wavelet), Moray Wavelet (Moll)
eye Wavelet), Mexican Wavelet (Mexican Wavelet), French Wavelet (French Wavelet), Shannon Wavelet (Shannon Wavelet),
Meyer Wavelet
t), Dove cheese wavelet (Daubechi)
es Wavelet), Simlet Wavelet, Corefret Wavelet, Cosplet Wavelet, Spline Wavelet
t) and the like are known. A and b are arbitrary real numbers, and along with two integers j and k,
【0044】と置くことにより、関数f(t)の離散化
表現の ## EQU3 ## By putting, the discrete expression of the function f (t) is
【0045】を得ることができる。更に、 Can be obtained. Furthermore,
【0046】と表現するとし、このとき信号f(t)を
f0 (t)と見做すと、上式は、 If the signal f (t) is regarded as f 0 (t) at this time, the above equation becomes:
【0047】と書くことができる。Can be written as
【0048】ウエーブレット変換処理は、周期的な信号
の分解や、パルスの分解や、ノイズの分解や、異常波の
検出や、ピークの検出等の種々の目的に使用されるが、
本実施形態においては、例えば4階のBスプラインウエ
ーブレット(SplineWavelet)をマザーウ
エーブレット(Mother Wavelet)として
用い、脈派波形データf(t)に対してf-2(t)を求
めればよい。更に、このように変換した波形から、極値
と、その極値の出現する時点を検出することができる。The wavelet transform process is used for various purposes such as periodic signal decomposition, pulse decomposition, noise decomposition, abnormal wave detection, and peak detection.
In the present embodiment, for example, the fourth floor B-spline wavelet (SplineWavelet) may be used as a mother wavelet (Mother Wavelet), and f -2 (t) may be obtained for the pulse waveform data f (t). Further, from the waveform converted in this way, the extreme value and the time when the extreme value appears can be detected.
【0049】図4は本発明の第四の実施形態を示すブロ
ック図である。FIG. 4 is a block diagram showing a fourth embodiment of the present invention.
【0050】図4の実施形態は、図1の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した膜派波形の
スムージング処理を行う脈派波形前処理部76aと、脈
派波形前処理部76aでスムージング処理された信号か
ら各時点における2階差分値を算出し、その2階差分値
が充分に小さな値の所定の値に等しいかまたは所定の値
よりも小さい場合を変曲点として変曲点の出現する時点
を検出する脈波変曲点自動検出部76bとを有する脈波
極値検出部76を設けたものである。その他の構成要素
は、図1の実施形態と同じである。In the embodiment of FIG. 4, instead of the pulse wave extremum detection unit 6 of the embodiment of FIG. 1, a pulse wave pre-processing unit 76a for performing a smoothing process of the membrane wave stored in the storage unit 5; When the second-order difference value at each time point is calculated from the signal subjected to the smoothing process in the pulse waveform pre-processing unit 76a, and the second-order difference value is equal to or smaller than a sufficiently small predetermined value. Is provided as an inflection point, and a pulse wave inflection point automatic detection unit 76b for detecting a time point at which the inflection point appears is provided. Other components are the same as those in the embodiment of FIG.
【0051】次に上述のように構成した脳波処理装置の
動作について、第一の実施形態と異なる部分の動作につ
いて説明する。Next, the operation of the electroencephalogram processing device configured as described above will be described, focusing on the operation of the portion different from that of the first embodiment.
【0052】脈派波形前処理部76aは、第二の実施形
態と同様に、記憶部5に格納した膜派波形のスムージン
グ処理を行う。The pulse waveform pre-processing unit 76a performs the smoothing process on the membrane waveform stored in the storage unit 5, as in the second embodiment.
【0053】脈波変曲点自動検出部76bは、脈派波形
前処理部76aでスムージング処理された信号から各時
点における2階差分値を算出し、その2階差分値が充分
に小さな値の所定の値に等しいかまたは所定の値よりも
小さい場合を変曲点として、変曲点の出現する時点を検
出する。具体的には、以下のようなデータ処理を行う。The pulse wave inflection point automatic detecting section 76b calculates the second-order difference value at each time from the signal smoothed by the pulse waveform pre-processing section 76a, and the second-order difference value is a sufficiently small value. The time when the inflection point appears is detected by using the case where the value is equal to or smaller than the predetermined value as the inflection point. Specifically, the following data processing is performed.
【0054】すなわち、スムージング処理された脈派の
デジタルデータを xi とし、そのi番目の時点におけ
る2階(後退)差分値 That is, the digital data of the pulse train subjected to the smoothing process is denoted by x i, and the second-order (retreat) difference value at the i-th time point
【0055】を計算する。例えば、図11の脈波データ
17の場合は、4つの変曲点17e・17f・17g・
17hがあり、変曲点17eおよび17gは Is calculated. For example, in the case of the pulse wave data 17 in FIG. 11, four inflection points 17e, 17f, 17g,
17h, and the inflection points 17e and 17g are
【0056】によって検出可能であり、変曲点17fお
よび17hは The inflection points 17f and 17h are
【0057】によって検出可能である。Can be detected.
【0058】なお、本実施形態において、例えば変曲点
17eを利用する場合の遡及時間は、第二の実施形態に
おける極大値17aを利用する場合の遡及時間とは異っ
た時間となる。In this embodiment, for example, the retroactive time when using the inflection point 17e is different from the retroactive time when using the maximum value 17a in the second embodiment.
【0059】図5は本発明の第五の実施形態を示すブロ
ック図である。FIG. 5 is a block diagram showing a fifth embodiment of the present invention.
【0060】図5の実施形態は、図1の実施形態の脈波
極値検出部6と脳波分割起点設定部7との間に、脈波極
値検出部6において検出された隣接する極値間の差とし
て脈波振幅値を算出する脈波振幅値算出部10aと、脈
波極値検出部6において検出された極値のなかから脈波
振幅値算出部10aによって算出された脈波振幅値が所
定の範囲内にある極値とその出現時点とを選出する加算
セグメント判断部10bとを有する判断部10を設けた
ものであり、その他の構成要素は、図1の実施形態と同
じである。The embodiment shown in FIG. 5 is arranged between the pulse wave extremum detection unit 6 and the electroencephalogram division starting point setting unit 7 of the embodiment shown in FIG. A pulse wave amplitude value calculating unit 10a for calculating a pulse wave amplitude value as a difference between the pulse wave amplitude value calculating unit 10a and a pulse wave amplitude calculated by the pulse wave amplitude value calculating unit 10a from the extreme values detected by the pulse wave extreme value detecting unit 6. A judgment unit 10 having an extreme value whose value is within a predetermined range and an addition segment judging unit 10b for selecting the present time is provided. Other components are the same as those in the embodiment of FIG. is there.
【0061】次に上述のように構成した脳波処理装置の
動作について、第一の実施形態と異なる部分の動作につ
いて説明する。Next, regarding the operation of the electroencephalogram processing device configured as described above, the operation of the portion different from that of the first embodiment will be described.
【0062】脈波振幅値算出部10aは、脈波極値検出
部6において検出された隣接する極値間の差として脈波
振幅値を算出する。The pulse wave amplitude value calculating section 10a calculates a pulse wave amplitude value as a difference between adjacent extreme values detected by the pulse wave extreme value detecting section 6.
【0063】加算セグメント判断部10bは、脈波極値
検出部6において検出された極値のなかから脈波振幅値
算出部10aによって算出された脈波振幅値が所定の範
囲内にある極値とその出現時点とを選出する。具体的に
は、以下のようなデータ処理を行う。The addition segment judging section 10b determines, from among the extreme values detected by the pulse wave extreme value detecting section 6, an extreme value whose pulse wave amplitude value calculated by the pulse wave amplitude value calculating section 10a is within a predetermined range. And its current time to elect. Specifically, the following data processing is performed.
【0064】すなわち、まず始めに、脈波振幅値の算出
を以下のように行う。That is, first, the calculation of the pulse wave amplitude value is performed as follows.
【0065】図11に示したような脈波データ17の1
周期分17nの時間幅(ウインドウ幅)の中で、脈波電
位値の最大値と最小値とを探索し、それらの差を脈波振
幅値とする。 図11の場合は、 脈波振幅値=|x(tA )−x(tD )| となる。但しtA は極大値17aが出現する時点、tD
は極小値17dが出現する時点である。One of the pulse wave data 17 as shown in FIG.
The maximum value and the minimum value of the pulse wave potential value are searched within the time width (window width) of the period 17n, and the difference between them is defined as the pulse wave amplitude value. In the case of FIG. 11, the pulse wave amplitude value = | x (t A ) −x (t D ) |. However, t A is the time when the maximum value 17a appears, t D
Is the time when the minimum value 17d appears.
【0066】このようにして算出した脈波振幅値は、影
響を受ける。すなわち、図12に示すように、心理状態
が緊張しているときは、図12(a)の振幅減少区間1
8aに示すように、振幅値が減少し、心理状態が弛緩し
ているときは、図12(b)の振幅漸増区間18bに示
すように、振幅値が増加することがよく知られている。
このような心理状態による影響は、心電データにも反映
されることがあり、このような心電データの変動は、脈
波データの分割の起点の変動を引き起しかねない。従っ
て、安定した心理状態のときのデータによって脈波デー
タの分割の起点の設定を行う必要があり、安定した心理
状態を検出するために、以下のようにして脈波振幅値を
求める。The pulse wave amplitude value calculated in this way is affected. That is, as shown in FIG. 12, when the mental state is nervous, the amplitude decreasing section 1 in FIG.
As shown in FIG. 8A, when the amplitude value decreases and the mental state is relaxing, it is well known that the amplitude value increases as shown in an amplitude gradual increase section 18b of FIG. 12B.
The influence of such a psychological state may be reflected on the electrocardiogram data, and such a fluctuation of the electrocardiogram data may cause a fluctuation of the starting point of the division of the pulse wave data. Therefore, it is necessary to set the starting point of the division of the pulse wave data based on the data in the stable mental state. In order to detect the stable mental state, the pulse wave amplitude value is obtained as follows.
【0067】図13は、心理状態が変動しているときの
脈波データを示す波形図である。図13において、脈波
データ20に対して脈波極値検出部6における処理を適
用すると、例えば極大値17aとして時点t1 ,t2 ,
……,t16が検出される。これらの時点の中で、安定し
た心理状態のときには脈波振幅値が所定の値ε以上であ
ると設定すると、t1 =t1 +,t2 =t2 +,t3 =
t3 +,t4 =t4 +,t5 =t5 +,t13=t6 +,t14=t
7 +,t15=t8 +,t16=t9 + が選出できる。FIG. 13 is a waveform diagram showing pulse wave data when the mental state is changing. In FIG. 13, when the processing in the pulse wave extremum detection unit 6 is applied to the pulse wave data 20, for example, the time points t 1 , t 2 ,
.., T 16 are detected. At these moments, if the pulse wave amplitude value is set to be equal to or greater than the predetermined value ε in a stable mental state, t 1 = t 1 + , t 2 = t 2 + , t 3 =
t 3 +, t 4 = t 4 +, t 5 = t 5 +, t 13 = t 6 +, t 14 = t
7 + , t 15 = t 8 + and t 16 = t 9 + can be selected.
【0068】図6は本発明の第六の実施形態を示すブロ
ック図である。FIG. 6 is a block diagram showing a sixth embodiment of the present invention.
【0069】図6の実施形態は、図5の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した膜派波形の
スムージング処理を行う脈派波形前処理部36aと、脈
派波形前処理部36aでスムージング処理された信号か
ら各時点における1階差分値を算出してその1階差分値
が充分に小さな値の所定の値に等しいかまたは所定の値
よりも小さい場合にその値を極値としてその極値の出現
する時点を検出する差分型脈波極値自動検出部36bと
を有する脈波極値検出部36を設けたものである。その
他の構成要素は、図5の実施形態と同じである。従って
本実施形態の動作は、図5の実施形態の動作の脈波極値
検出部6の動作の代りに、図2の実施形態の脈派波形前
処理部36aと差分型脈波極値自動検出部36bとの動
作を置換したものである。In the embodiment of FIG. 6, instead of the pulse wave extremum detection unit 6 of the embodiment of FIG. 5, a pulse wave pre-processing unit 36a for performing a smoothing process of the membrane wave stored in the storage unit 5, When the first-order difference value at each time is calculated from the signal subjected to the smoothing processing by the pulse waveform pre-processing unit 36a, and the first-order difference value is equal to or smaller than a sufficiently small predetermined value. And a pulse-wave extremal-value detection unit 36 having a differential-type pulse-wave extremum-automatic detection unit 36b that detects the time when the extremum appears as an extreme value. Other components are the same as those in the embodiment of FIG. Therefore, the operation of the present embodiment is different from the operation of the pulse wave extremum detector 6 of the operation of the embodiment of FIG. 5 in that the pulse wave pre-processing unit 36a of the embodiment of FIG. This is a replacement of the operation with the detection unit 36b.
【0070】図7は本発明の第七の実施形態を示すブロ
ック図である。FIG. 7 is a block diagram showing a seventh embodiment of the present invention.
【0071】図7の実施形態は、図5の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した膜派波形の
スムージング処理を行う脈派波形前処理部56aと、脈
派波形前処理部56aでスムージング処理された信号に
対してウエーブレット処理を施すことによって極値を求
めてその極値の出現する時点を検出するウエーブレット
変換型脈波極値自動検出部56bとを有する脈波極値検
出部56を設けたものである。その他の構成要素は、図
5の実施形態と同じである。従って本実施形態の動作
は、図5の実施形態の動作の脈波極値検出部6の動作の
代りに、図3の実施形態の脈派波形前処理部56aとウ
エーブレット変換型脈波極値自動検出部56bとの動作
を置換したものである。In the embodiment of FIG. 7, instead of the pulse wave extremum detection unit 6 of the embodiment of FIG. 5, a pulse wave pre-processing unit 56a for performing smoothing processing of the membrane wave stored in the storage unit 5, Wavelet conversion type pulse wave extremum automatic detection unit 56b for obtaining an extremum by performing wavelet processing on the signal smoothed by the pulse waveform preprocessing unit 56a and detecting the time when the extremum appears And a pulse wave extremum detection unit 56 having the following. Other components are the same as those in the embodiment of FIG. Therefore, the operation of the present embodiment is different from the operation of the pulse wave extremum detection unit 6 of the operation of the embodiment of FIG. 5 in that the pulse wave pre-processing unit 56a of the embodiment of FIG. This is a replacement of the operation with the automatic value detection unit 56b.
【0072】図8は本発明の第八の実施形態を示すブロ
ック図である。FIG. 8 is a block diagram showing an eighth embodiment of the present invention.
【0073】図8の実施形態は、図5の実施形態の脈波
極値検出部6の代りに、記憶部5に格納した膜派波形の
スムージング処理を行う脈派波形前処理部76aと、脈
派波形前処理部76aでスムージング処理された信号か
ら各時点における2階差分値を算出し、その2階差分値
が充分に小さな値の所定の値に等しいかまたは所定の値
よりも小さい場合を変曲点として変曲点の出現する時点
を検出する脈波変曲点自動検出部76bとを有する脈波
極値検出部76を設けたものである。その他の構成要素
は、図5の実施形態と同じである。従って本実施形態の
動作は、図5の実施形態の動作の脈波極値検出部6の動
作の代りに、図4の実施形態の脈派波形前処理部76a
と脈波変曲点自動検出部76bとの動作を置換したもの
である。In the embodiment of FIG. 8, instead of the pulse wave extremum detection unit 6 of the embodiment of FIG. 5, a pulse wave pre-processing unit 76a for performing a smoothing process of the membrane wave stored in the storage unit 5, When the second-order difference value at each time point is calculated from the signal subjected to the smoothing process in the pulse waveform pre-processing unit 76a, and the second-order difference value is equal to or smaller than a sufficiently small predetermined value. Is provided as an inflection point, and a pulse wave inflection point automatic detection unit 76b for detecting a time point at which the inflection point appears is provided. Other components are the same as those in the embodiment of FIG. Therefore, the operation of the present embodiment is different from the operation of the pulse wave extremum detecting unit 6 of the operation of the embodiment of FIG. 5 in that the pulse waveform preprocessing unit 76a of the embodiment of FIG.
And the operation of the pulse wave inflection point automatic detection unit 76b is replaced.
【0074】[0074]
【発明の効果】以上説明したように、本発明の脳波処理
装置は、2極間の電位差としての脳波波形を検出するた
めの活性電極および基準電極を有する脳波センサと、耳
朶または指先から脈派波形を検出する脈派センサと、脳
波センサの活性電極と基準電極との電位差および脈派セ
ンサの電位とを増幅する増幅器と、この増幅器の出力の
アナログ信号をデジタル信号に変換するアナログデジタ
ル変換器と、アナログデジタル変換器の出力信号を格納
する記憶部と、記憶部に格納した信号から周期性を有す
る脈派波形の極値およびその極値の出現する時点を検出
する脈波極値検出部と、脈波極値検出部の出力信号によ
って所定の時間だけ遡及した時点を脳波波形の分割の起
点に設定する脳波分割起点設定部と、脳波分割起点設定
部によって設定された起点から脳波波形を個々のセグメ
ントに分割し、セグメントを順次加算して平均値を求め
ることによって脳波波形に雑音として混入している心電
波形の推定波形を求める雑音混入推定部と、雑音混入推
定部で得られた推定波形を起点を同期させて脳波波形か
ら差し引くことによって雑音の混入を除去した脳波信号
を抽出する脳波信号抽出部とを備えることにより、心電
波形検出用の電極の装着よりもセンサの装着が簡便な脈
派を利用し、心電波形の特徴的な波形に対応する脈派の
特徴的な波形の出現時点を検出することが可能になるた
め、心電波形の推定波形を求めるための脳波波形のセグ
メント分割を容易にし、脳波波形から心電波形を除去す
るための計測時間を短縮して脳波処理を容易にすること
ができるという効果があるAs described above, the electroencephalogram processing device according to the present invention includes an electroencephalogram sensor having an active electrode and a reference electrode for detecting an electroencephalogram waveform as a potential difference between two electrodes, and a pulse wave from an earlobe or a fingertip. A pulse sensor that detects a waveform, an amplifier that amplifies the potential difference between the active electrode and the reference electrode of the electroencephalogram sensor and the potential of the pulse sensor, and an analog-to-digital converter that converts an analog signal output from the amplifier into a digital signal A storage unit for storing an output signal of the analog-to-digital converter, and a pulse wave extremum detection unit for detecting, from the signal stored in the storage unit, an extremum of a pulse waveform having periodicity and a time point at which the extremum appears A brain wave division starting point setting unit that sets a point in time that has been traced back by a predetermined time by an output signal of the pulse wave extreme value detection unit as a starting point of the division of the electroencephalogram waveform, and a brain wave division starting point setting unit. A noise mixing estimator that divides an electroencephalogram waveform into individual segments from the starting point and sequentially adds the segments to obtain an average value, thereby obtaining an estimated waveform of an electrocardiogram waveform mixed as noise into the electroencephalogram waveform; By attaching an electrode for detecting an electrocardiographic waveform by providing an electroencephalogram signal extraction unit for extracting an electroencephalogram signal from which noise has been removed by synchronizing the starting point of the estimated waveform obtained by the estimating unit and subtracting the estimated waveform from the electroencephalogram waveform, Estimation of the electrocardiographic waveform because it is possible to use the pulse that is easier to attach the sensor than to detect the characteristic waveform of the pulse corresponding to the characteristic waveform of the electrocardiographic waveform. The effect is that the segmentation of the electroencephalogram waveform for obtaining the waveform can be easily performed, and the measurement time for removing the electrocardiogram waveform from the electroencephalogram waveform can be shortened to facilitate the electroencephalogram processing.
【図1】本発明の第一の実施形態を示すブロック図であ
る。FIG. 1 is a block diagram showing a first embodiment of the present invention.
【図2】本発明の第二の実施形態を示すブロック図であ
る。FIG. 2 is a block diagram showing a second embodiment of the present invention.
【図3】本発明の第三の実施形態を示すブロック図であ
る。FIG. 3 is a block diagram showing a third embodiment of the present invention.
【図4】本発明の第四の実施形態を示すブロック図であ
る。FIG. 4 is a block diagram showing a fourth embodiment of the present invention.
【図5】本発明の第五の実施形態を示すブロック図であ
る。FIG. 5 is a block diagram showing a fifth embodiment of the present invention.
【図6】本発明の第六の実施形態を示すブロック図であ
る。FIG. 6 is a block diagram showing a sixth embodiment of the present invention.
【図7】本発明の第七の実施形態を示すブロック図であ
る。FIG. 7 is a block diagram showing a seventh embodiment of the present invention.
【図8】本発明の第八の実施形態を示すブロック図であ
る。FIG. 8 is a block diagram showing an eighth embodiment of the present invention.
【図9】従来の脳波処理装置の一例を示すブロック図で
ある。FIG. 9 is a block diagram illustrating an example of a conventional electroencephalogram processing device.
【図10】実測された脳波および心電波形および脈波の
一例を示す波形図である。FIG. 10 is a waveform chart showing an example of actually measured brain waves, electrocardiographic waveforms, and pulse waves.
【図11】図10の脈波の一周期分を示す波形図であ
る。FIG. 11 is a waveform diagram showing one cycle of the pulse wave of FIG. 10;
【図12】心理状態の影響を受ける脈波の実測データの
一例を示す波形図である。FIG. 12 is a waveform chart showing an example of actual measurement data of a pulse wave affected by a mental state.
【図13】心理状態の影響を受けた脈波データにおける
波形の変動の一例を示す波形図である。FIG. 13 is a waveform chart showing an example of waveform fluctuation in pulse wave data affected by a mental state.
1a 活性電極 1b 基準電極 2a・2b 脈派センサ 3 増幅器(アンプ) 4 アナログデジタル変換器(A/D) 5 記憶部 6・36・56・76 脈波極値検出部 6a 波形表示部 7 脳波分割起点設定部 8 雑音混入推定部 9 脳波信号抽出部 10 判断部 10a 脈波振幅値算出部 10b 加算セグメント判断部 15 脳波データ 16 心電データ 17 脈波データ 17a・17c 極大値 17b・17d 極小値 17e〜17h 変曲点 18a 振幅減少区間 18b 振幅漸増区間 20 脈波データ 36a・56a.76a 脈派波形前処理部 36b 差分型脈波極値自動検出部 56b ウエーブレット変換型脈波極値自動検出部 76b 脈波変曲点自動検出部 91a・91b 脳波計測用電極 92a・92b 心電図用電極 93 脳電位増幅器 94 アナログデジタル変換器 95 一時記憶部 96 心電図用増幅器 98 演算部 99 トリガー信号発生部 1a Active electrode 1b Reference electrode 2a ・ 2b Pulse sensor 3 Amplifier (Amplifier) 4 Analog / Digital converter (A / D) 5 Storage unit 6.36.56.76 Pulse wave extremum detection unit 6a Waveform display unit 7 EEG division Starting point setting unit 8 Noise mixing estimation unit 9 EEG signal extraction unit 10 Judgment unit 10a Pulse wave amplitude value calculation unit 10b Addition segment judgment unit 15 EEG data 16 Electrocardiogram data 17 Pulse wave data 17a / 17c Maximum value 17b / 17d Minimum value 17e 17h Inflection point 18a Amplitude decreasing section 18b Amplitude increasing section 20 Pulse wave data 36a / 56a. 76a pulse wave pre-processing unit 36b differential type pulse wave extreme value automatic detection unit 56b wavelet conversion type pulse wave extreme value automatic detection unit 76b pulse wave inflection point automatic detection unit 91a / 91b EEG measurement electrodes 92a / 92b for electrocardiogram Electrode 93 Brain potential amplifier 94 Analog-to-digital converter 95 Temporary storage unit 96 ECG amplifier 98 Operation unit 99 Trigger signal generation unit
Claims (10)
するための活性電極および基準電極を有する脳波センサ
と、耳朶または指先から脈派波形を検出する脈派センサ
と、前記脳波センサの前記活性電極と前記基準電極との
電位差および前記脈派センサの電位とを増幅する増幅器
と、前記増幅器の出力のアナログ信号をデジタル信号に
変換するアナログデジタル変換器と、前記アナログデジ
タル変換器の出力信号を格納する記憶部と、前記記憶部
に格納した信号から周期性を有する膜派波形の極値およ
びその極値の出現する時点を検出する脈波極値検出部
と、前記脈波極値検出部の出力信号によって所定の時間
だけ遡及した時点を前記脳波波形の分割の起点に設定す
る脳波分割起点設定部と、前記脳波分割起点設定部によ
って設定された起点から前記脳波波形を個々のセグメン
トに分割し前記セグメントを順次加算して平均値を求め
ることによって前記脳波波形に雑音として混入している
心電波形の推定波形を求める雑音混入推定部と、前記雑
音混入推定部で得られた前記推定波形を前記起点を同期
させて前記脳波波形から差し引くことによって雑音の混
入を除去した脳波信号を抽出する脳波信号抽出部とを備
えることを特徴とする脳波処理装置。1. An electroencephalogram sensor having an active electrode and a reference electrode for detecting an electroencephalogram waveform as a potential difference between two poles, a pulse sensor for detecting a pulse waveform from an earlobe or a fingertip, and the electroencephalogram sensor. An amplifier for amplifying a potential difference between an active electrode and the reference electrode and a potential of the pulse sensor, an analog-to-digital converter for converting an analog signal output from the amplifier to a digital signal, and an output signal of the analog-to-digital converter A pulse wave extremum detection unit that detects an extreme value of a periodic membranous waveform having a periodicity from the signal stored in the storage unit and a time point at which the extremum appears, and a pulse wave extremum detection An electroencephalogram division starting point setting unit that sets a point in time that has been traced back by a predetermined time by the output signal of the unit as the starting point of the division of the electroencephalogram waveform; A noise mixing estimating unit for obtaining an estimated waveform of an electrocardiographic waveform mixed as noise into the brain wave waveform by dividing the brain wave waveform into individual segments and sequentially adding the segments to obtain an average value; and An electroencephalogram processing device comprising: an electroencephalogram signal extraction unit configured to extract an electroencephalogram signal from which noise is removed by synchronizing the estimated waveform obtained by the contamination estimation unit with the start point and subtracting the estimated waveform from the electroencephalogram waveform. .
置の画面に表示し、前記画面から周期性を有する脈派波
形の極値およびその極値の出現する時点を直接検出する
波形表示部を有する前記脈波極値検出部を備えることを
特徴とする請求項1記載の脳波処理装置。2. A waveform display for displaying a pulse waveform stored in the storage unit on a screen of a display device, and directly detecting an extreme value of the periodic pulse waveform and a time point at which the extreme value appears from the screen. The electroencephalogram processing device according to claim 1, further comprising the pulse wave extremum detection unit having a unit.
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号から各時点におけ
る1階差分値を算出して前記1階差分値が充分に小さな
値の所定の値に等しいかまたは前記所定の値よりも小さ
い場合にその値を極値としてその極値の出現する時点を
検出する差分型脈波極値自動検出部とを有する前記脈波
極値検出部を備えることを特徴とする請求項1記載の脳
波処理装置。3. A pulse waveform preprocessing unit for performing a smoothing process on the pulse waveform stored in the storage unit, and calculating a first-order difference value at each time from the signal smoothed by the pulse waveform preprocessing unit. When the first-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value, the difference type pulse wave pole detecting the time when the extreme value appears as the extreme value The electroencephalogram processing apparatus according to claim 1, further comprising the pulse wave extremum detection unit having an automatic value detection unit.
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号に対してウエーブ
レット処理を施すことによって極値を求めて前記極値の
出現する時点を検出するウエーブレット変換型脈波極値
自動検出部とを有する前記脈波極値検出部を備えること
を特徴とする請求項1記載の脳波処理装置。4. A pulse wave pre-processing unit for performing a smoothing process on a pulse waveform stored in the storage unit, and performing a wavelet process on the signal smoothed by the pulse wave pre-processing unit. The electroencephalogram processing according to claim 1, further comprising the pulse wave extremum detection unit having a wavelet transformation type pulse wave extremum automatic detection unit that obtains an extremum and detects a time when the extremum appears. apparatus.
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号から各時点におけ
る2階差分値を算出して前記2階差分値が充分に小さな
値の所定の値に等しいかまたは前記所定の値よりも小さ
い場合を変曲点として前記変曲点の出現する時点を検出
する脈波変曲点自動検出部とを有する前記脈波極値検出
部を備えることを特徴とする請求項1記載の脳波処理装
置。5. A pulse waveform pre-processing unit for performing a smoothing process on a pulse waveform stored in the storage unit, and calculating a second-order difference value at each time from a signal smoothed by the pulse waveform pre-processing unit. A pulse wave inflection point automatic detecting the time when the inflection point appears as an inflection point when the second-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value. The electroencephalogram processing device according to claim 1, further comprising the pulse wave extremum detection unit having a detection unit.
するための活性電極および基準電極を有する脳波センサ
と、耳朶または指先から脈派波形を検出する脈派センサ
と、前記脳波センサの前記活性電極と前記基準電極との
電位差および前記脈派センサの電位とを増幅する増幅器
と、前記増幅器の出力のアナログ信号をデジタル信号に
変換するアナログデジタル変換器と、前記アナログデジ
タル変換器の出力信号を格納する記憶部と、前記記憶部
に格納した信号から周期性を有する膜派波形の極値およ
びその極値の出現する時点を検出する脈波極値検出部
と、前記脈波極値検出部において検出された隣接する極
値間の差として脈波振幅値を算出する脈波振幅値算出部
と、前記脈波極値検出部において検出された前記極値の
なかから前記脈波振幅値算出部によって算出された前記
脈波振幅値が所定の範囲内にある極値とその出現時点と
を選出する加算セグメント判断部とを有する判断部と、
前記判断部からの出力信号によって所定の時間だけ遡及
した時点を前記脳波波形の分割の起点に設定する脳波分
割起点設定部と、前記脳波分割起点設定部によって設定
された起点から前記脳波波形を個々のセグメントに分割
し前記セグメントを順次加算して平均値を求めることに
よって前記脳波波形に雑音として混入している心電波形
の推定波形を求める雑音混入推定部と、前記雑音混入推
定部で得られた前記推定波形を前記起点を同期させて前
記脳波波形から差し引くことによって雑音の混入を除去
した脳波信号を抽出する脳波信号抽出部とを備えること
を特徴とする脳波処理装置。6. An electroencephalogram sensor having an active electrode and a reference electrode for detecting an electroencephalogram waveform as a potential difference between two poles, a pulse wave sensor for detecting a pulse wave waveform from an earlobe or a fingertip, and the electroencephalogram sensor. An amplifier for amplifying a potential difference between an active electrode and the reference electrode and a potential of the pulse sensor, an analog-to-digital converter for converting an analog signal output from the amplifier to a digital signal, and an output signal of the analog-to-digital converter A pulse wave extremum detection unit that detects an extreme value of a periodic membranous waveform having a periodicity from the signal stored in the storage unit and a time point at which the extremum appears, and a pulse wave extremum detection A pulse wave amplitude value calculating unit that calculates a pulse wave amplitude value as a difference between adjacent extreme values detected by the unit, and the pulse wave amplitude from the extreme values detected by the pulse wave extreme value detecting unit. A determination unit having an addition segment determination unit that selects the extreme value and the present time when the pulse wave amplitude value calculated by the value calculation unit is within a predetermined range,
An electroencephalogram division starting point setting section that sets a point in time that is retroactive by a predetermined time by the output signal from the determination section as a starting point of the electroencephalogram waveform division, and individually calculates the electroencephalogram waveform from the starting point set by the electroencephalogram division starting point setting section. And a noise mixing estimating unit for obtaining an estimated waveform of an electrocardiographic waveform mixed as noise in the brain wave waveform by sequentially adding the segments to obtain an average value. An electroencephalogram signal processing unit, which extracts an electroencephalogram signal from which the noise is removed by subtracting the estimated waveform from the electroencephalogram waveform by synchronizing the starting point.
置の画面に表示し、前記画面から周期性を有する脈派波
形の極値およびその極値の出現する時点を直接検出する
波形表示部を有する前記脈波極値検出部を備えることを
特徴とする請求項6記載の脳波処理装置。7. A waveform display for displaying a pulse waveform stored in the storage unit on a screen of a display device, and directly detecting an extreme value of the periodic pulse waveform and a time when the extreme value appears from the screen. The electroencephalogram processing device according to claim 6, further comprising the pulse wave extremum detection unit having a unit.
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号から各時点におけ
る1階差分値を算出して前記1階差分値が充分に小さな
値の所定の値に等しいかまたは前記所定の値よりも小さ
い場合にその値を極値としてその極値の出現する時点を
検出する差分型脈波極値自動検出部とを有する前記脈波
極値検出部を備えることを特徴とする請求項6記載の脳
波処理装置。8. A pulse waveform preprocessing unit for performing a smoothing process on the pulse waveform stored in the storage unit, and calculating a first-order difference value at each time from the signal smoothed by the pulse waveform preprocessing unit. When the first-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value, the difference type pulse wave pole detecting the time when the extreme value appears as the extreme value The electroencephalogram processing device according to claim 6, further comprising the pulse wave extremum detection unit having an automatic value detection unit.
ジング処理を行う脈派波形前処理部と、前記脈派波形前
処理部でスムージング処理された信号に対してウエーブ
レット処理を施すことによって極値を求めて前記極値の
出現する時点を検出するウエーブレット変換型脈波極値
自動検出部とを有する前記脈波極値検出部を備えること
を特徴とする請求項6記載の脳波処理装置。9. A pulse waveform pre-processing unit for performing a smoothing process on a pulse waveform stored in the storage unit, and performing a wavelet process on the signal smoothed by the pulse waveform pre-processing unit. 7. The electroencephalogram processing according to claim 6, further comprising the pulse wave extremum detection unit having a wavelet transformation type pulse wave extremum automatic detection unit for obtaining an extremum and detecting a time when the extremum appears. apparatus.
ージング処理を行う脈派波形前処理部と、前記脈派波形
前処理部でスムージング処理された信号から各時点にお
ける2階差分値を算出して前記2階差分値が充分に小さ
な値の所定の値に等しいかまたは前記所定の値よりも小
さい場合を変曲点として前記変曲点の出現する時点を検
出する脈波変曲点自動検出部とを有する前記脈波極値検
出部を備えることを特徴とする請求項6記載の脳波処理
装置。10. A pulse waveform pre-processing unit for performing a smoothing process of a pulse waveform stored in the storage unit, and calculating a second-order difference value at each time from a signal smoothed by the pulse waveform pre-processing unit. A pulse wave inflection point automatic detecting the time when the inflection point appears as an inflection point when the second-order difference value is equal to a sufficiently small predetermined value or smaller than the predetermined value. The electroencephalogram processing device according to claim 6, further comprising the pulse wave extremum detection unit having a detection unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP9069753A JP3033512B2 (en) | 1997-03-24 | 1997-03-24 | EEG processing device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP9069753A JP3033512B2 (en) | 1997-03-24 | 1997-03-24 | EEG processing device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH10262941A true JPH10262941A (en) | 1998-10-06 |
JP3033512B2 JP3033512B2 (en) | 2000-04-17 |
Family
ID=13411878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP9069753A Expired - Fee Related JP3033512B2 (en) | 1997-03-24 | 1997-03-24 | EEG processing device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3033512B2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100689987B1 (en) | 2005-09-08 | 2007-03-09 | 경희대학교 산학협력단 | Apparatus and method for removing noise included in an EEG signal in an EEG and MRI system |
JP2007111118A (en) * | 2005-10-18 | 2007-05-10 | Japan Lifeline Co Ltd | Brain lesion diagnosis apparatus and method of using the same |
JP2007128237A (en) * | 2005-11-02 | 2007-05-24 | Tokyo Denki Univ | Extreme value detection method and extreme value detection program using optimal smoothing spline |
JP2008229307A (en) * | 2007-03-21 | 2008-10-02 | Korea Advanced Inst Of Sci Technol | EEG measurement method, EEG measurement apparatus, and recording medium |
CN101933805A (en) * | 2009-06-29 | 2011-01-05 | 索尼公司 | Auricle-installed device and bio-signal measurement apparatus |
JP2011110378A (en) * | 2009-11-30 | 2011-06-09 | Qrs:Kk | Electroencephalograph and electroencephalography |
CN103654766A (en) * | 2013-11-15 | 2014-03-26 | 中国航天员科研训练中心 | Low-override amplification device for weak physiological signals |
-
1997
- 1997-03-24 JP JP9069753A patent/JP3033512B2/en not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100689987B1 (en) | 2005-09-08 | 2007-03-09 | 경희대학교 산학협력단 | Apparatus and method for removing noise included in an EEG signal in an EEG and MRI system |
JP2007111118A (en) * | 2005-10-18 | 2007-05-10 | Japan Lifeline Co Ltd | Brain lesion diagnosis apparatus and method of using the same |
JP2007128237A (en) * | 2005-11-02 | 2007-05-24 | Tokyo Denki Univ | Extreme value detection method and extreme value detection program using optimal smoothing spline |
JP2008229307A (en) * | 2007-03-21 | 2008-10-02 | Korea Advanced Inst Of Sci Technol | EEG measurement method, EEG measurement apparatus, and recording medium |
JP2012030111A (en) * | 2007-03-21 | 2012-02-16 | Samsung Electronics Co Ltd | Method and instrument for measuring brain wave, and recording medium |
CN101933805A (en) * | 2009-06-29 | 2011-01-05 | 索尼公司 | Auricle-installed device and bio-signal measurement apparatus |
CN101933805B (en) | 2009-06-29 | 2012-08-08 | 索尼公司 | Auricle-installed device and bio-signal measurement apparatus |
JP2011110378A (en) * | 2009-11-30 | 2011-06-09 | Qrs:Kk | Electroencephalograph and electroencephalography |
CN103654766A (en) * | 2013-11-15 | 2014-03-26 | 中国航天员科研训练中心 | Low-override amplification device for weak physiological signals |
Also Published As
Publication number | Publication date |
---|---|
JP3033512B2 (en) | 2000-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7035685B2 (en) | Apparatus and method for measuring electroencephalogram | |
AU2006308467B2 (en) | Method and system for high-resolution extraction of quasi-periodic signals | |
WO2007041766A1 (en) | Adaptive real-time line noise suppression for electrical or magnetic physiological signals | |
JP2015519087A (en) | Real-time QRS detection using adaptive threshold | |
JPH0626541B2 (en) | EEG spectrum analyzer | |
JPH07284482A (en) | Heart rate variability waveform analysis method and device | |
KR101536361B1 (en) | Method for analyzing ecg and ecg apparatus therefor | |
WO2013052944A1 (en) | Periodic artifact reduction from biomedical signals | |
WO2009110051A1 (en) | On-line vital signal processing method effective during fetal movement and its apparatus | |
JP3033512B2 (en) | EEG processing device | |
Kouekeu et al. | Embedded QRS complex detection based on ECG signal strength and trend | |
Luengo et al. | Sparse ECG representation with a multi-scale dictionary derived from real-world signals | |
US4716907A (en) | Apparatus for detecting electroencephalogram and evoked response with monopolar derivation method | |
JP3864390B2 (en) | EEG analysis method | |
JPH10216096A (en) | Biological signal analyzer | |
Burke et al. | ECG analysis using the Mexican-Hat wavelet | |
JP3314521B2 (en) | Heart rate variability waveform analysis method and apparatus | |
Govreen-Segal et al. | Real-time PC-based system for dynamic beat-to-beat QT-RR analysis | |
US20210052179A1 (en) | Method for determining r peaks of electrocardiogram | |
CN112826451A (en) | A method and device for evaluating the depth of anesthesia and the depth of sleep | |
RU2063167C1 (en) | Method and device for producing electric cardiogram | |
US8731647B2 (en) | Method and system for analyzing the cardiac activity of a patient and uses thereof | |
JP3307071B2 (en) | Heart rate variability waveform analysis method and apparatus | |
WO2024248049A1 (en) | Electrocardiograph | |
Jayashree et al. | Optimising Seismocardiography for Precision Cardiac Diagnosis with Advanced Signal Processing Algorithms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20000118 |
|
LAPS | Cancellation because of no payment of annual fees |