JPH0261250B2 - - Google Patents
Info
- Publication number
- JPH0261250B2 JPH0261250B2 JP59133705A JP13370584A JPH0261250B2 JP H0261250 B2 JPH0261250 B2 JP H0261250B2 JP 59133705 A JP59133705 A JP 59133705A JP 13370584 A JP13370584 A JP 13370584A JP H0261250 B2 JPH0261250 B2 JP H0261250B2
- Authority
- JP
- Japan
- Prior art keywords
- peak
- signal
- time
- pattern
- order differential
- 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.)
- Expired - Lifetime
Links
- 238000001914 filtration Methods 0.000 claims description 12
- 238000000034 method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- 238000003672 processing method Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 description 20
- 238000012937 correction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000002051 biphasic effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000002902 bimodal effect Effects 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は心電図等の生体信号の認識、計測を行
う場合の生体信号処理方式に関する。
う場合の生体信号処理方式に関する。
一般に生体の信号は、ある形状を持つ複数種の
波の組合せ、あるいはそれらの周期的な繰り返し
で構成されているが、その生体信号の性質を把握
するにはそれらの棘波の基線（ゼロレベル）から
の立ち上り時刻、基線へ戻る時刻、ピーク生起時
刻、あるいは変曲点を示す時刻等を計測し、波の
幅（時間間隔）や高さ（電圧値）を求めて、棘波
の形状を類推する方法が一般に行われている。し
かし生体信号は、その電圧レベルが小さくハムな
どの高周波雑音が重畳し易いためその信号の計測
を正確に行うことは難かしいとされていた。例え
ば心電図を例にとると心電図信号は第２図に示す
ようにP.QRS.T.U波と称される各棘波から構成
されているが、この心電図信号を自動解析するあ
たつて従来は、その周波数成分が高く、かつ振幅
値が大きい為、検出が容易であるQRS波を中心
として行われており、Ｐ波あるいはＴ波の認識に
関しては、雑音や基線動揺による影響のため正確
な波形認識、計測は困難であつた。
波の組合せ、あるいはそれらの周期的な繰り返し
で構成されているが、その生体信号の性質を把握
するにはそれらの棘波の基線（ゼロレベル）から
の立ち上り時刻、基線へ戻る時刻、ピーク生起時
刻、あるいは変曲点を示す時刻等を計測し、波の
幅（時間間隔）や高さ（電圧値）を求めて、棘波
の形状を類推する方法が一般に行われている。し
かし生体信号は、その電圧レベルが小さくハムな
どの高周波雑音が重畳し易いためその信号の計測
を正確に行うことは難かしいとされていた。例え
ば心電図を例にとると心電図信号は第２図に示す
ようにP.QRS.T.U波と称される各棘波から構成
されているが、この心電図信号を自動解析するあ
たつて従来は、その周波数成分が高く、かつ振幅
値が大きい為、検出が容易であるQRS波を中心
として行われており、Ｐ波あるいはＴ波の認識に
関しては、雑音や基線動揺による影響のため正確
な波形認識、計測は困難であつた。
特に、不整脈などの解析の際にはＰ波の存在の
有無、形状などの正確な認識が必要となるがＰ波
は周波数成分も低く、また振幅値も小さいため筋
電図やハムなどの混入で信号成分が雑音に埋もれ
てしまいその弁別が難かしいとされていた。また
Ｔ波については、心筋異常等を診断する場合に、
Ｓ−Ｔ成分の傾向やＴ波の形状（陰性波や２相性
の区別）を認識することが必要であるが、Ｓ−Ｔ
成分には雑音が重畳し易く、また陰性（２相性）
Ｔ波の認識には基線動揺の鉱響を受け易いため、
計測誤差が大きくなるという欠点があつた。
有無、形状などの正確な認識が必要となるがＰ波
は周波数成分も低く、また振幅値も小さいため筋
電図やハムなどの混入で信号成分が雑音に埋もれ
てしまいその弁別が難かしいとされていた。また
Ｔ波については、心筋異常等を診断する場合に、
Ｓ−Ｔ成分の傾向やＴ波の形状（陰性波や２相性
の区別）を認識することが必要であるが、Ｓ−Ｔ
成分には雑音が重畳し易く、また陰性（２相性）
Ｔ波の認識には基線動揺の鉱響を受け易いため、
計測誤差が大きくなるという欠点があつた。
尚、従来このＰ波（Ｔ波）を検出するには以下
の方式がとられている。(1)心電図アナログ信号を
アナログデイジタル変換し、そのデイジタル信号
に対し、帯域通過型のデイジタルフイルタリング
処理を行うことによりＰ波やＴ波の生起時刻、終
端時刻、ピーク生起時刻を決定する方式、(2)心電
図アナログ信号をアナログデイジタル変換し、そ
のデイジタル原信号に対し、差分処理を加えて一
次微係数を求め傾きの変化が最大となる時刻をＰ
波やＴ波の生起時刻とする方式である。方式(1)に
於て使用されるデイジタルフイルタをソフトウエ
アで実現すると係数として精度が要求される小数
点データを必要とし、さらに高精度の乗除算処理
を繰り返し実行するため、処理時間がかかり、必
要なメモリも大容量のものを必要とする欠点があ
つた。また、Ｐ波、Ｔ波を弁別するためには異な
つた周波数特性を持つデイジタルフイルタを必要
とした。方(2)は原心電図信号に重畳した筋電図、
ハムなど高周波雑音による影響を受け易く誤つた
認識を行う可能性が大きかつた。
の方式がとられている。(1)心電図アナログ信号を
アナログデイジタル変換し、そのデイジタル信号
に対し、帯域通過型のデイジタルフイルタリング
処理を行うことによりＰ波やＴ波の生起時刻、終
端時刻、ピーク生起時刻を決定する方式、(2)心電
図アナログ信号をアナログデイジタル変換し、そ
のデイジタル原信号に対し、差分処理を加えて一
次微係数を求め傾きの変化が最大となる時刻をＰ
波やＴ波の生起時刻とする方式である。方式(1)に
於て使用されるデイジタルフイルタをソフトウエ
アで実現すると係数として精度が要求される小数
点データを必要とし、さらに高精度の乗除算処理
を繰り返し実行するため、処理時間がかかり、必
要なメモリも大容量のものを必要とする欠点があ
つた。また、Ｐ波、Ｔ波を弁別するためには異な
つた周波数特性を持つデイジタルフイルタを必要
とした。方(2)は原心電図信号に重畳した筋電図、
ハムなど高周波雑音による影響を受け易く誤つた
認識を行う可能性が大きかつた。
本発明の目的は、従来の生体信号処理方式と比
較して高周波雑音や基線動揺の影響を受けずに波
形の生起時刻、終端時刻、ピーク生起時刻を正確
にしかも同一構成で求めることができ、さらに波
形の形状を２次微分出力信号のパターンからだけ
で認識できる生体信号処理方式を提供することで
ある。
較して高周波雑音や基線動揺の影響を受けずに波
形の生起時刻、終端時刻、ピーク生起時刻を正確
にしかも同一構成で求めることができ、さらに波
形の形状を２次微分出力信号のパターンからだけ
で認識できる生体信号処理方式を提供することで
ある。
本発明によれば、処理対象信号に２次微分フイ
ルタリング処理を加える操作と、２次微分出力さ
れた時系列信号群よりピークを検出する操作と、
検出されたピーク信号系列の生起パターンにより
原信号のパターンを類推し、分類する操作と、２
次微分ピーク信号系列より原信号の生起時刻、終
端時刻およびピーク生起時刻を求めさらに原信号
のピーク電圧値を計測する操作とにより構成され
る生体信号処理方式が得られる。
ルタリング処理を加える操作と、２次微分出力さ
れた時系列信号群よりピークを検出する操作と、
検出されたピーク信号系列の生起パターンにより
原信号のパターンを類推し、分類する操作と、２
次微分ピーク信号系列より原信号の生起時刻、終
端時刻およびピーク生起時刻を求めさらに原信号
のピーク電圧値を計測する操作とにより構成され
る生体信号処理方式が得られる。
生体信号等の処理対象信号の波形の一周期の基
本形を第３図ａと考えると、ｂはａ及びａと極性
の反転する波形の連続したものであり、ｃはａの
波形が２つ連続した波形と考えることができる。
本形を第３図ａと考えると、ｂはａ及びａと極性
の反転する波形の連続したものであり、ｃはａの
波形が２つ連続した波形と考えることができる。
さらに一般の連続波形はこれらの波形が組み合
せられて連続的な周期を持つて現われると考えて
よい。そこでこの処理対象信号に２次微分操作を
加える。本発明にて利用する２次微分フイルタの
式を式(1)に示す。
せられて連続的な周期を持つて現われると考えて
よい。そこでこの処理対象信号に２次微分操作を
加える。本発明にて利用する２次微分フイルタの
式を式(1)に示す。
Ｆ（xk）＝１／100T2｛−２・（xk＋x-k）−２・
（xk+1＋xk-1）
＋０・（xk+2＋xk-2）＋０・（xk+3＋xk-2）
＋１・（xk+4＋xk-4）＋２・（xk+5＋xk-5）
＋１・（xk+6＋xk-6）｝ ………(1)
ここで｛xk｝は第４図に示すサンプリングされ
た各時系列データである。又１／100T2はスケール フアクタを示す。式(1)に示した通り本フイルタは
係数が簡単な整数（０，１ないしは２）を示すた
め、加減算およびビツトシフト演算のみで処理が
可能であり、高速かつ、良好な低減２次微分特性
を実現する特長を有している。
た各時系列データである。又１／100T2はスケール フアクタを示す。式(1)に示した通り本フイルタは
係数が簡単な整数（０，１ないしは２）を示すた
め、加減算およびビツトシフト演算のみで処理が
可能であり、高速かつ、良好な低減２次微分特性
を実現する特長を有している。
第３図ａ，ｂ，ｃの形状を持つ信号に上記フイ
ルタリング操作を加えた出力結果を各々a′，b′，
c′に示す。
ルタリング操作を加えた出力結果を各々a′，b′，
c′に示す。
得られた２次微分フイルタリング出力信号のピ
ーク生起時刻は各々原信号における生起時刻、終
端時刻および生起時刻ないしはその近傍の時刻を
示しており、しかも上記の良好な低減特性より原
信号に高周波雑音等が重畳した場合でも、その影
響を受けずに顕著なピークを示すため雑音から信
号成分を弁別することが容易である。また、原信
号に基線動揺がある場合も、２次微分出力信号の
ピーク生起時刻に影響を及ぼすことはないので原
原波形の計測が正確に行われる。さらに原信号の
うち第３図ａ，ｂ，ｃの成分に対する２次微分出
力信号成分のピークが他の成分より顕著であるこ
とよりa′，b′，c′に示したピーク列の生起パター
ンを認識することにより原波形ａ，ｂ，ｃの形状
を類推することが可能である。ここで心電図を例
にとると心電図を構成する棘波成分のうちＰ波、
Ｔ波は一般的に第３図ａ，ｂ，ｃのパターンのい
ずれかを示しａを単峰性パターンｂを２相性パタ
ーンｃを双峰性パターンと称する。各ａ，ｂ，ｃ
の信号においてＰ（Ｔ）波の２次微分出力パター
ンの一般形をc′とするとa′，b′は各々c′のP1ない
しＬはP5あるいはP1およびP5が存在しない場合
の特殊波形と考えられるのでＰ波（Ｔ波）の２次
微分フイルタリング出力信号をそのピーク数、ピ
ーク電圧値ピーク時間間隔等からa′，b′，c′のパ
ターンに分類することにより心電図信号解析に必
要な区分と称する各Ｐ（Ｔ）波の生起時刻終端時
刻、ピーク生起時刻を計測可能であり、さらにそ
の分類されたパターンからＰ（Ｔ）波の形状が単
峰性、２相性、双峰性かを判断することが可能と
なる。
ーク生起時刻は各々原信号における生起時刻、終
端時刻および生起時刻ないしはその近傍の時刻を
示しており、しかも上記の良好な低減特性より原
信号に高周波雑音等が重畳した場合でも、その影
響を受けずに顕著なピークを示すため雑音から信
号成分を弁別することが容易である。また、原信
号に基線動揺がある場合も、２次微分出力信号の
ピーク生起時刻に影響を及ぼすことはないので原
原波形の計測が正確に行われる。さらに原信号の
うち第３図ａ，ｂ，ｃの成分に対する２次微分出
力信号成分のピークが他の成分より顕著であるこ
とよりa′，b′，c′に示したピーク列の生起パター
ンを認識することにより原波形ａ，ｂ，ｃの形状
を類推することが可能である。ここで心電図を例
にとると心電図を構成する棘波成分のうちＰ波、
Ｔ波は一般的に第３図ａ，ｂ，ｃのパターンのい
ずれかを示しａを単峰性パターンｂを２相性パタ
ーンｃを双峰性パターンと称する。各ａ，ｂ，ｃ
の信号においてＰ（Ｔ）波の２次微分出力パター
ンの一般形をc′とするとa′，b′は各々c′のP1ない
しＬはP5あるいはP1およびP5が存在しない場合
の特殊波形と考えられるのでＰ波（Ｔ波）の２次
微分フイルタリング出力信号をそのピーク数、ピ
ーク電圧値ピーク時間間隔等からa′，b′，c′のパ
ターンに分類することにより心電図信号解析に必
要な区分と称する各Ｐ（Ｔ）波の生起時刻終端時
刻、ピーク生起時刻を計測可能であり、さらにそ
の分類されたパターンからＰ（Ｔ）波の形状が単
峰性、２相性、双峰性かを判断することが可能と
なる。
そのためまず心電図原信号に関する２次微分出
力信号よりQRS棘波に関する成分を除去ないし
はQRSの生起時刻を明確にしてＰ（Ｔ）波検出の
際にQRS棘波の影響を及ぼさないようにする。
これはQRS棘波群がその２次フイルタ出力にお
いて他の出力よりも顕著な電圧値を示すためであ
る。次に、２次微分出力信号に重畳している雑音
を計測し、その平均電圧レベルを求める。原信号
に混入した図やハムの電圧レベルが高い場合、２
次微分出力信号中においても雑音レベルが大きく
なるため、これがＰ（Ｔ）波信号成分の検出精度
に影響を及ぼさないようにピーク検出の為の閾値
を動的に設定できる電圧レベルによつて決めてい
る。この電圧レベルに相当する雑音閾値によつて
雑音レベル以上にあるピーク系列を選択しその中
からさらに電圧レベルの絶対値が最大を示す時刻
のピークを中心として前後に存在するピークを複
数個、選択する。この選択基準としてはまず、最
大電圧レベルを示す時刻の近傍一定時間間隔内に
あるピークを選び、ピークとピークの時間間隔が
定められた時間間隔の範囲に存在しないものは雑
音ないしはノツチとみなしふるい落す。こうして
残つたピーク系列のパターンが第３図ａ，ｂ，ｃ
に示した２次微分出力パターンに合致するかを判
別する。もし、ピーク系列パターンがどれかに合
致する場合は２次微分出力信号ピーク生起時刻が
各々原信号の生起時刻、終端時刻、ピーク生起時
刻ないしはこの近傍を示している。原信号の性質
によつて、例えばＴ波の終端のようにその変化が
急峻でない場合、ピーク生起時刻が実際の区分点
よりずれを生じるのを補正する。０１ピーク系列
パターンにおいてピークが過不足を示し第３図
a′，b′，c′のいずれにも合致しない場合、これを
補正してパターンの分類を行う。例えば第５図に
示すような２次微分出力信号上で原波形のピーク
と生起時刻を示すパターンが得られた場合、b′に
相当するピークを再検出する為に雑音の閾値を下
げてb′を求める。もし、モデルパターンよりもピ
ーク数が多い場合には時間間隔、レベル等で、ピ
ークの再検出を行い、ピーク数を減少させる。
力信号よりQRS棘波に関する成分を除去ないし
はQRSの生起時刻を明確にしてＰ（Ｔ）波検出の
際にQRS棘波の影響を及ぼさないようにする。
これはQRS棘波群がその２次フイルタ出力にお
いて他の出力よりも顕著な電圧値を示すためであ
る。次に、２次微分出力信号に重畳している雑音
を計測し、その平均電圧レベルを求める。原信号
に混入した図やハムの電圧レベルが高い場合、２
次微分出力信号中においても雑音レベルが大きく
なるため、これがＰ（Ｔ）波信号成分の検出精度
に影響を及ぼさないようにピーク検出の為の閾値
を動的に設定できる電圧レベルによつて決めてい
る。この電圧レベルに相当する雑音閾値によつて
雑音レベル以上にあるピーク系列を選択しその中
からさらに電圧レベルの絶対値が最大を示す時刻
のピークを中心として前後に存在するピークを複
数個、選択する。この選択基準としてはまず、最
大電圧レベルを示す時刻の近傍一定時間間隔内に
あるピークを選び、ピークとピークの時間間隔が
定められた時間間隔の範囲に存在しないものは雑
音ないしはノツチとみなしふるい落す。こうして
残つたピーク系列のパターンが第３図ａ，ｂ，ｃ
に示した２次微分出力パターンに合致するかを判
別する。もし、ピーク系列パターンがどれかに合
致する場合は２次微分出力信号ピーク生起時刻が
各々原信号の生起時刻、終端時刻、ピーク生起時
刻ないしはこの近傍を示している。原信号の性質
によつて、例えばＴ波の終端のようにその変化が
急峻でない場合、ピーク生起時刻が実際の区分点
よりずれを生じるのを補正する。０１ピーク系列
パターンにおいてピークが過不足を示し第３図
a′，b′，c′のいずれにも合致しない場合、これを
補正してパターンの分類を行う。例えば第５図に
示すような２次微分出力信号上で原波形のピーク
と生起時刻を示すパターンが得られた場合、b′に
相当するピークを再検出する為に雑音の閾値を下
げてb′を求める。もし、モデルパターンよりもピ
ーク数が多い場合には時間間隔、レベル等で、ピ
ークの再検出を行い、ピーク数を減少させる。
次に本発明の一実施例を図面を参照して説明す
る。第１図は本発明の一実施例を示す図で、１は
２次微分デイジタルフイルタであり、心電図等生
体信号にフイルタリング操作を加える。２は雑音
検出装置であり、デイジタルフイルタによつて出
力された信号より雑音等解析に不要な２次微分出
力成分の信号の平均電圧値を検出するものであ
る。３はQRS波検出装置であり２次微分出力信
号よりQRS波に相当する成分を検出する。４は
減算装置であり、２つの入力信号の差分を行う。
５は雑音閾値レベル発生装置であり、ピーク検出
の為の閾値を決定する。６はピーク検出装置であ
り、設定された閾値以上にあるピークを検出す
る。７は電圧レベル比較装置であり、ピーク検出
装置６で検出されたピーク電圧のレベルを比較す
る。８は最大レベル検出装置であり入力されたピ
ークより最大ピークを選択する。９はピークパタ
ーン生成装置であり、ピーク検出装置６で検出さ
れたピーク列に基づきパターンを生成する。１０
はパターン比較装置であり、ピークパターン生成
装置９で得られたパターンをあらかじめ設定され
たモデルパターンと比較する。１１はピーク補正
装置であり、ピーク生起時刻の補正を行いピーク
時刻とピーク電圧値を出力する。１２は閾値補正
装置でありパターン比較装置１０でパターン比較
された結果に基きもし、パターンがモデルパター
ンと異なる場合に閾値レベルを変化させる。１３
は原波形の生起時刻、終端時刻の補装置であり、
正確な区分点の時刻を出力する。１４は信号遅延
装置であり、入力信号を一定時間、遅延する。
る。第１図は本発明の一実施例を示す図で、１は
２次微分デイジタルフイルタであり、心電図等生
体信号にフイルタリング操作を加える。２は雑音
検出装置であり、デイジタルフイルタによつて出
力された信号より雑音等解析に不要な２次微分出
力成分の信号の平均電圧値を検出するものであ
る。３はQRS波検出装置であり２次微分出力信
号よりQRS波に相当する成分を検出する。４は
減算装置であり、２つの入力信号の差分を行う。
５は雑音閾値レベル発生装置であり、ピーク検出
の為の閾値を決定する。６はピーク検出装置であ
り、設定された閾値以上にあるピークを検出す
る。７は電圧レベル比較装置であり、ピーク検出
装置６で検出されたピーク電圧のレベルを比較す
る。８は最大レベル検出装置であり入力されたピ
ークより最大ピークを選択する。９はピークパタ
ーン生成装置であり、ピーク検出装置６で検出さ
れたピーク列に基づきパターンを生成する。１０
はパターン比較装置であり、ピークパターン生成
装置９で得られたパターンをあらかじめ設定され
たモデルパターンと比較する。１１はピーク補正
装置であり、ピーク生起時刻の補正を行いピーク
時刻とピーク電圧値を出力する。１２は閾値補正
装置でありパターン比較装置１０でパターン比較
された結果に基きもし、パターンがモデルパター
ンと異なる場合に閾値レベルを変化させる。１３
は原波形の生起時刻、終端時刻の補装置であり、
正確な区分点の時刻を出力する。１４は信号遅延
装置であり、入力信号を一定時間、遅延する。
次に本実施例の動作を説明する。まず心電図信
号第６図ａは２次微分デイジタルフイルタ１でフ
イルタリング操作が行われる。このフイルタリン
グ操作においては第６図ｂに示すような出力信号
が得られるので雑音検出装置２においてｂのt1区
間に存在する雑音を検出し雑音閾値レベル発生装
置５においてピーク検出装置６内に必要な雑音除
去のための閾値電圧を設定する。QRS波検出装
置３、減算装置４ではＰ（Ｔ）波検出の際、不要
なQRS棘波の２次微分出力成分第６図区間t2を除
去する。ピーク検出装置６においては、減算装置
４からの信号より、設定された雑音閾値th1以上
にあるピークを検出する。これを第７図に示す。
さらに電圧レベル比較装置７、最大レベル検出装
置８でピーク系列のうちの最大電圧レベルにある
ものを選択する。（第７図Ｐ）パターン比較装置
１０はピークパターン生成装置９において得られ
た第８図に示すようなピークの生成パターンが第
３図に示したパターンのいずれかと一致するか否
かを比較する。もし、パターンが一致する場合は
ピーク補正装置１１、時刻の補正装置１３へトリ
ガを発生する。ピーク補正装置１１では最大レベ
ル検出装置８で得られた２次微分出力最大ピーク
レベルの時刻に基づき原波形の性質により必要あ
ればその時刻の補正を行い、原波形を参照し原波
形のピーク生起時刻及びその電圧値を出力する。
号第６図ａは２次微分デイジタルフイルタ１でフ
イルタリング操作が行われる。このフイルタリン
グ操作においては第６図ｂに示すような出力信号
が得られるので雑音検出装置２においてｂのt1区
間に存在する雑音を検出し雑音閾値レベル発生装
置５においてピーク検出装置６内に必要な雑音除
去のための閾値電圧を設定する。QRS波検出装
置３、減算装置４ではＰ（Ｔ）波検出の際、不要
なQRS棘波の２次微分出力成分第６図区間t2を除
去する。ピーク検出装置６においては、減算装置
４からの信号より、設定された雑音閾値th1以上
にあるピークを検出する。これを第７図に示す。
さらに電圧レベル比較装置７、最大レベル検出装
置８でピーク系列のうちの最大電圧レベルにある
ものを選択する。（第７図Ｐ）パターン比較装置
１０はピークパターン生成装置９において得られ
た第８図に示すようなピークの生成パターンが第
３図に示したパターンのいずれかと一致するか否
かを比較する。もし、パターンが一致する場合は
ピーク補正装置１１、時刻の補正装置１３へトリ
ガを発生する。ピーク補正装置１１では最大レベ
ル検出装置８で得られた２次微分出力最大ピーク
レベルの時刻に基づき原波形の性質により必要あ
ればその時刻の補正を行い、原波形を参照し原波
形のピーク生起時刻及びその電圧値を出力する。
一方、時刻の補正装置１３ではピーク検出装置
６、ピークパターン生成装置９からのピークパタ
ーンより、原波形の生起時刻、終端時刻を検出し
原波形の性質により、時刻の補正を行い、出力す
る。もし、パターン比較装置においてモデルパタ
ーンと一致しない場合には、閾値補正装置１２に
トリガをかけピーク検出用の雑音閾値を変化させ
雑音閾値レベル発生装置５に戻る。
６、ピークパターン生成装置９からのピークパタ
ーンより、原波形の生起時刻、終端時刻を検出し
原波形の性質により、時刻の補正を行い、出力す
る。もし、パターン比較装置においてモデルパタ
ーンと一致しない場合には、閾値補正装置１２に
トリガをかけピーク検出用の雑音閾値を変化させ
雑音閾値レベル発生装置５に戻る。
以上の様な方式を採用することにより従来の方
式と比較して次のような効果がある。(1)本方式で
採用した２次微分デイジタルフイルタは加減算器
および遅延回路だけの簡単な回路構成で実現でき
るため従来のデイジタルフイルタを利用したＰ
（Ｔ）波検出方式に比較して高速な波形検出が可
能である。(2)本方式で採用した２次微分デイジタ
ルフイルタは良好な低域通過特性を持つため、ハ
ム、筋電信号等の高周波雑音を含んだ信号に対し
ても高い検出精度を示す。(3)２次微分フイルタリ
ング出力信号のみから、原信号の生起時刻、ピー
ク生起時刻、終端時刻およびその近傍が検知でき
るので原信号の保持、原信号の雑音処理等のため
の手段が不要になり、構成が簡単になる。(4)外部
から与えるパラメータの設定を変えるのみで同一
の構成でＰ波、Ｔ波の検出が可能であるため、信
号処理装置全体の構成が簡単になる。
式と比較して次のような効果がある。(1)本方式で
採用した２次微分デイジタルフイルタは加減算器
および遅延回路だけの簡単な回路構成で実現でき
るため従来のデイジタルフイルタを利用したＰ
（Ｔ）波検出方式に比較して高速な波形検出が可
能である。(2)本方式で採用した２次微分デイジタ
ルフイルタは良好な低域通過特性を持つため、ハ
ム、筋電信号等の高周波雑音を含んだ信号に対し
ても高い検出精度を示す。(3)２次微分フイルタリ
ング出力信号のみから、原信号の生起時刻、ピー
ク生起時刻、終端時刻およびその近傍が検知でき
るので原信号の保持、原信号の雑音処理等のため
の手段が不要になり、構成が簡単になる。(4)外部
から与えるパラメータの設定を変えるのみで同一
の構成でＰ波、Ｔ波の検出が可能であるため、信
号処理装置全体の構成が簡単になる。
第１図は本発明の一実施例を示す図、第２図は
生体信号の例としての心電図信号を示す図、第３
図は生体信号の処理の原理を示す図、第４図は時
系列で示された信号データを示す図、第５図は補
正された信号波形を示す図、第６図は心電図信号
と２次微分フイルタリング処理された信号を示す
図、第７図はピーク検出の様子を示す図、第８図
は生成ピークパターンを示す図。
生体信号の例としての心電図信号を示す図、第３
図は生体信号の処理の原理を示す図、第４図は時
系列で示された信号データを示す図、第５図は補
正された信号波形を示す図、第６図は心電図信号
と２次微分フイルタリング処理された信号を示す
図、第７図はピーク検出の様子を示す図、第８図
は生成ピークパターンを示す図。
Claims (1)
- １ 入力生体信号を受け２次微分フイルタリング
処理するフイルタ手段と、前記フイルタ手段の出
力を受けピークを検出するピーク検出手段と、前
記ピーク検出手段からのピーク検出結果の時系列
的変化により入力生体信号を所定の分類に類推す
る分類手段と、前記分類手段の出力を受け２次微
分フイルタリング処理された信号を基に前記入力
生体信号におけるピーク生起時刻及びピーク電圧
値を計測する計測手段とを具備することを特徴と
する生体信号処理方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59133705A JPS6111017A (ja) | 1984-06-28 | 1984-06-28 | 生体信号処理方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59133705A JPS6111017A (ja) | 1984-06-28 | 1984-06-28 | 生体信号処理方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6111017A JPS6111017A (ja) | 1986-01-18 |
| JPH0261250B2 true JPH0261250B2 (ja) | 1990-12-19 |
Family
ID=15110961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59133705A Granted JPS6111017A (ja) | 1984-06-28 | 1984-06-28 | 生体信号処理方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6111017A (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5230106B2 (ja) * | 2007-01-15 | 2013-07-10 | 富士フイルム株式会社 | 超音波診断装置、ｉｍｔ計測方法及びｉｍｔ計測プログラム |
-
1984
- 1984-06-28 JP JP59133705A patent/JPS6111017A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6111017A (ja) | 1986-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US3940692A (en) | Apparatus for monitoring recurrent waveforms | |
| US5772603A (en) | Device for filtering ECG signals | |
| CA1258099A (en) | Egg with delta wave recognition by double | |
| JPS63226338A (ja) | Ｅｃｇ波信号からペース・パルス信号を除去する方法とｅｃｇ波信号からペース・パルス信号を除去する装置 | |
| US4846568A (en) | Method of analyzing an electroretinogram | |
| US4718097A (en) | Method and apparatus for determining the endpoints of a speech utterance | |
| US6820017B1 (en) | Method for determining the amplitude and phase angle of a measuring signal corresponding to a current or voltage of an electrical power supply network | |
| JPS60246737A (ja) | 生体信号用信号調整器 | |
| Martínez et al. | A robust T wave alternans detector based on the GLRT for Laplacian noise distribution | |
| EP0058050B1 (en) | Measuring method | |
| US4321932A (en) | Electrode impedance monitoring method apparatus for electrocardiography | |
| JPH0261250B2 (ja) | ||
| US3991748A (en) | Circuit arrangement for the processing of physiological measuring signals | |
| JPS644328Y2 (ja) | ||
| US3581219A (en) | System for dynamically adjusting clipping in editing electrocardiogram waves | |
| Froning et al. | Problems and limitations of ECG baseline estimation and removal using a cubic spline technique during exercise ECG testing: recommendations for proper implementation | |
| JPH05212006A (ja) | 心拍間隔計測装置 | |
| JPH0238214B2 (ja) | ||
| Reddy et al. | Analysis of interval measurements on CSE multilead reference ECGs | |
| JPS58223927A (ja) | デ−タ圧縮処理装置 | |
| JPH03203412A (ja) | 信号波形の動揺の補正方法 | |
| Dotsinsky et al. | Continuously Tested and Used QRS Detection Algorithm: Free Access to the MATLAB Code | |
| RU2122349C1 (ru) | Способ обнаружения qrs-комплекса и определения момента его появления | |
| RU2242164C2 (ru) | Способ выявления информативных параметров st-сегмента электрокардиосигнала и устройство для его осуществления | |
| Pushkareva et al. | The Research of Noise Stability of the Prony's Method While Processing of Cardiological Time Series |