JPS6111017A - 生体信号処理方式 - Google Patents

生体信号処理方式

Info

Publication number
JPS6111017A
JPS6111017A JP59133705A JP13370584A JPS6111017A JP S6111017 A JPS6111017 A JP S6111017A JP 59133705 A JP59133705 A JP 59133705A JP 13370584 A JP13370584 A JP 13370584A JP S6111017 A JPS6111017 A JP S6111017A
Authority
JP
Japan
Prior art keywords
peak
signal
time
order differential
pattern
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP59133705A
Other languages
English (en)
Other versions
JPH0261250B2 (ja
Inventor
克郎 岡本
赤井 孝至
将史 伊藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
Nippon Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Electric Co Ltd filed Critical Nippon Electric Co Ltd
Priority to JP59133705A priority Critical patent/JPS6111017A/ja
Publication of JPS6111017A publication Critical patent/JPS6111017A/ja
Publication of JPH0261250B2 publication Critical patent/JPH0261250B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は心電図等の生体信号の認識、計測を行う場合の
生体信号処理方式に関する。
〔従来の技術〕
一般に生体の信号は、るる形状を持つ複数種の波の組合
せ、あるいはそれらの周期的な繰り返しで構成されてい
るが、その生体信号の性質を把握するにはそれらの検波
の基線(ゼロレベル)からの立ち上り時刻、基線へ戻る
時刻、ピーク生起時刻、あるいは変曲点を示す時刻等を
計測し、波の幅(時間間隔)や高さく電圧値)を求めて
、蜂波の形状を類推する方法が一般に行われている。
しかし生体信号は、その電圧レベルが小さくハムなどの
高周波雑音が重畳し易いためその信号の計測を正確に行
うことは難かしいとされていた。例えば心電図を例にと
ると心電図信号は第2図に示すようにP、 Q几S、 
T、 U波と称される各線波から構成されているが、こ
の心電図信号を自動解析するにらたって従来は、その周
波数成分が高く、かつ振幅値が大きい為、検出が容易で
あるQR8波を中心として行われており、P波あるいは
T波の認識に関しては、雑音や基線動揺による影響のた
め正確な波形認識、計測は困難であった。
特に、不整脈などの解析の際にはP波の存在の有無、形
状などの正確な認識が必要となるがP波は周波数成分も
低く、また振幅値も小さいため筋電図やハムなどの混入
で信号成分が雑音に埋もれてしまいその弁別が難かしい
とされていた。またT波については、心筋異常等を診断
する場合に、S−T成分の傾向やT波の形状(陰性波や
2相性の区別)を認識することが必要であるが、S−T
成分には雑音が重畳し易く、また陰性(2相生)T波の
認識には基線動揺の鉱響を受け易いため、計測誤差が大
きくなるという欠点があった。
冑、従来このP波(T波)を検出するには以下の方式が
とられている。(1)心電図アナ−ログ信号をアナログ
ディジタル変換し、そのディジタル信号に対し、帯域通
過型のディジタルフィルタリング処理を行うことにより
P波やT波の生起時刻、終端時刻、ピーク生起時刻を決
定する方式、(2)心電図アナログ信号をアナログディ
ジタル変換し、そのディジタル原信号に対し、差分処理
を加えて一次微係数を求め傾きの変化が最大となる時刻
をP波やT波の生起時刻とする方式である。方式(1)
に於て使用されるディジタルフィルタをソフトウェアで
実現すると係数として精度が要求される小数点データを
必要とし、さらに高精度の乗除算処理を繰り返し実行す
るため、処理時間がかかり、必要なメモリも大容量のも
のを必要とする欠点があった。また、P波、T波を弁別
するためには異なった周波数特性を持つディジタルフィ
ルタを必要とした。方式(2)は原石電図信号に重畳し
た筋電図、ハムなど高周波雑音による影響を受け易く誤
った認識を行う可能性が大きかった。
〔発明が解決しようとする問題点〕
本発明の目的は、従来の生体信号処理方式と比較して高
周波雑音や基線動揺の影響を受けずに波形の生起時刻、
終端時刻、ピーク生起時刻全正確にしかも同一構成で求
めることができ、さらに波形の形状を2次微分出力信号
のパターンからだけで認識できる生体信号処理方式を提
供することである。
〔問題点を解決するための手段〕
本発明によれば、処理対象信号に2次微分フィルタリン
グ処理を加える操作と、2次微分出力された時系列信号
群よりピークを検出する操作と、検出されたピーク信号
系列の生起パターンにより原信号のパターンを類推し、
分類する操作と、2次微分ピーク信号系列より原信号の
生起時刻、終端時刻およびピーク生起時刻を求めさらに
原信号のピーク電圧値を計1Qtlする操作とにより構
成される生体信号処理方式が得られる。
〔発明の原理〕
生体信号等の処理対象信号の波形の一周期の基本形を第
3図(a)と考えると、ω)は(a)及び(a)と極性
の反転する波形の連続したものであり、(C)は(a)
の波形が2つ連続した波形と考えることができる5、さ
らに一般の連続波形はこれらの波形が組み合せられて連
続的な周期を持って現われると考えてよい。そこでこの
処理対象信号に2次微分操作を加える。本発明にて利用
する2次微分フィルタの式を式(1)に示す。
+0 ・(JC&+2+3’! 2)+O・(”&+3
 +”k−3)+1 ・(”&+4+”&−4)+2 
・(”*+5+”1−s)ここで(xk )は第4図に
示すサンプリングされたアクタを示す。式(1)に示し
た通シ本フィルタは係数が簡単な整数(0,工ないしは
2)を示すため、加減算およびピットシフト演算のみで
処理が可能であり、高速かつ、良好な低減2次微分特性
を実現する特長を有している。
第3図(a)山)(C)の形状を持つ信号に上記フィル
タリング操作を加えた出力結果を各々(a’ Xb’ 
X(!’ )に示す。
得られた2次微分フィルタリング出力信号のピーク生起
時刻は各々原信号における生起時刻、終端時刻および生
起時刻ないしはその近傍の時刻を示しておシ、シかも上
記の良好な低減特性より原信号に高周波雑音等が重畳し
た場合でも、その影響を受けずに顕著なピークを示すた
め雑音から信号成分を弁別することが容易である。また
、原信号に基線動揺がある場合も、2次微分出力信号の
ピーク生起時刻に゛影響を及ぼすことはないので凍原波
形の計測が正確に行われる。さらに原信号のうち第3図
(a)(b)(C)の成分に対する2次微分出力信号成
分のピークが他の成分より顕著であることよ#)(a’
)(bつ(C′)に示しだピーク列の生起パターンを認
識することにより原波形(a)(b)(e)の形状を類
推することが可能である。ここで心電図を例にとると心
電図を構成する陣波成分のうちP波、T波は一般的に第
3図(a)(b)(e)のパターンのいずれかを示しく
a)を単峰性パターン(b)を2相性パターン(C)を
双峰性パターンと称する。各(a)(b)(e)の信号
においてP(T)波の2次微分出カバターンの一般形を
(C′)とすると(a’)(b’)は各々(C′)のP
lないしLはP5あるいはPlおよびP5が存在しない
場合の特殊波形と考えられるのでP波(T波)の2次微
分フィルタリング出力信号をそのピーク数、ピーク電圧
値ピーク時間間隔等から(a’) (b’) ((!’
)のパターンに分類することによシ心電図信号解析に必
要な区分と称する各P(T)波の生起時刻終端時刻、ビ
iり生起時刻を計測可能であり、さらにその分類された
パターンからP(T)波の形状が単峰性、2相性、双峰
性かを判断することが可能となる。
そのためまず心電図原信号に関する2次微分出力信号よ
りQ RSK波に関する成分を除去ないしはQR8の生
起時刻を明確にしてP(T)液検出の際にQR8R波の
影響を及ぼさないようにする。
これはBsn波群がその2次フィルタ出力において他の
出力よりも顕著な電圧値を示すためである。次に、2次
微分出力信号に重畳している雑音を計測し、その平均電
圧レベルを求める。原信号に混入した図やハムの電圧レ
ベルが高い場合、2次微分出力信号中においても雑音レ
ベルが大きくなるため、これがP(T)波信号成分の検
出精度に影響を及ぼさないようにピーク検出の為の閾値
を動的に設定できる電圧レベルによって決めている。
この電圧レベルに相当する雑音閾値によって雑音レベル
以上にあるピーク系列を選択しその中からさらに電圧レ
ベルの絶対値が最大を示す時刻のピークを中心として前
後に存在するピークを複数個、選択する。この選択基準
としてはまず、最大電圧レベルを示す時刻の近傍一定時
間間隔内にらるピークを選び、ピークとピークの時間間
隔が定められた時間間隔の範囲に存在しないものは雑音
ないしはノツチとみなしふるい落す。こうして残ったピ
ーク系列のパターンが第3図(a)(b)(C)に示し
た2次微分出カバターンに合致するかを判別する。もし
、ピーク系列パターンがどれかに合致する場合は2次微
分出力信号ピーク生起時刻が各々原信号の生起時刻、終
端時刻、ピーク生起時刻ないしはこの近傍を示している
。原信号の性質によって、例えばT波の終端のようにそ
の変化が急岐でない場合、ピーク生起時刻が実際の区分
点よりずれを生じるのを補正する。01ビ一ク系列パタ
ーンにおいてピークが過不足を示し第3図(a’) (
b’) (eつのいずれにも合致しない場合、これを補
正してパターンの分類を行う。例えば第5図に示すよう
な2次微分出力信号上で原波形のピークと生起時刻を示
すパターンが得られた場合、b′  に相当するピーク
を再検出する為に雑音の閾値を下げてb′を求める。も
し、モデルパターンよりもピーク数が多い場合には時間
間隔、レベル等で、ピークの再検出を行い、ピーク数を
減少させる。
〔実施例〕
次に本発明の一実施例を図面を参照して説明する。第1
図は本発明の一実施例を示す図で、1は2次微分ディジ
タルフィルタであり、心電図等生体信号にフィルタリン
グ操作を加える。2は雑音検出装置であり、ディジタル
フィルタによって出力された信号より雑音等解析に不要
な2次微分出力成分の信号の平均電圧値を検出するもの
である。
3はQR8波検出装置であり2次微分出力信号よりQR
8波に相当する成分を検出する。4は減算装置であシ、
2つの入力信号の差分を行う。5は雑音閾値レベル発生
装置であり、ピーク検出の為の閾値を決定する。6はピ
ーク検出装置であり、設定された閾値以上にあるピーク
を検出する。7は電圧レベル比較装置であり、ピーク検
出装置6で検出されたピーク電圧のレベルを比較する。
8は最大レベル検出装置であり入力されたピークよシ最
大ピークを選択する。9はピークパターン生成装置であ
シ、ピーク検出装置tt、6で検出されたピーク列に基
づきパターンを生成する。10はパターン比較装置であ
り、ピークパターン生成装置9で得られたパターンをあ
らかじめ設定されたモデルパターンと比較する。11は
ピーク補正装置であシ、ピーク生起時刻の補正を行いピ
ーク時刻とビーク電圧値を出力する。12は閾値補正装
置であシバターン比較装置10でパターン比較された結
果に基きもし、パターンがモデルパターンと異なる場合
に閾値レベルを変化させる。13は原波形の生起時刻、
終端時刻の補正装置であシ、正確な区分点の時刻を出力
する。14は信号遅延装置であり、入力1言号を一定時
間、遅延する。
次に本実施例の動作を説明する。まず心電図信号第6図
(a)は2次微分ディジタルフィルタ1でフィルタリン
グ操作が行われる。このフィルタリング操作においては
第6図中)に示すような出力信号が得られるので雑音検
出装置2において(b)の11区間に存在する雑音を検
出し雑音閾値レベル発生装置5においてピーク検出装置
6内に必要な雑音除去のだめの閾値電圧を設定する。Q
R8R8量検出装置3算装置4ではP(T)波検出の際
、不要なQR8R2O3次微分出力成分第6図区間t2
を除去する。ピーク検出装#6においては、t2算装置
4からの信号よシ、設定された雑音閾値tb+以上にお
るピークを検出する。これを第7図に示す。さらに電圧
レベル比較装置7、最大レベル検出装置8でピーク系列
のうちの最大電圧レベルにあるものを選択する。(第7
図P)パターン比較装置10はピークパターン生成装置
9において得られた第8図に示すようなピークの生成パ
ターンが第3図に示したパターンのいずれかと一致する
か否かを比較する。もし、パターンが一致する場合はピ
ーク補正装置11、時刻の補正装置13ヘトリガを発生
する。ピーク補正装置11では最大レベル検出装置8で
得られた2次微分出力最大ピークレベルの時刻に基づき
原波形の性質によシ必要あればその時刻の補正を行い、
原波形を参照し原波形のピーク生起時刻及びその電圧値
を出力する。
一方、時刻の補正装置13ではピーク検出装置6、ピー
クバター/生成装置9からのピークパターンよシ、原波
形の生起時刻、終端時刻を検出し原波形の性質により、
時刻の補正を行い、出力する。もし、パターン比較装置
においてモデルパターンと一致しない場合には、閾値補
正装置12にトリガ奮かけピーク検出用の雑音閾値を変
化させ雑音閾値レベル発生装置5に戻る。
〔発明の効果〕
以上の様な方式を採用することによプ従米の方式と比較
して次のような効果がある。(1)本方式で採用した2
次微分ディジタルフィルタは加減45および遅延回路だ
けの簡単なlr、!l路構成で実現できるだめ従来のデ
ィジタルフィルタを利用したP■波波山出方式比較して
高速な波形検出が可能である。(2)本方式で採用した
2次微分グイジタルフイルタは良好な低域通過特性を持
つため、ハム、筋電信号等の高周波雑音を含んだ信号に
対しても高い検出精度を示す。(3)2次微分フィルタ
リング出力信号のみから、原信号の生起時刻、ピーク生
起時刻、終端時刻およびその近傍が検知できるので原信
号の保持、原信号の雑音処理等のための手段が不要にな
シ、構成が簡単になる。(4)外部から与えるパラメー
タの設定を変えるのみで同一の構成でP波、T波の検出
が可能であるため、信号処理装置全体の構成が簡単にな
る。
【図面の簡単な説明】
第1図は本発明の一実施例を示す図、第2図は生体信号
の例としての心電図信号を示す図、第3図は生体信号の
処理の原理を示す図、第4図は時系列で示された信号デ
ータを示す図、第51’Jは補正された信号波形を示す
図、第6図は心電図信号と2次微分フィルタリング処理
された信号を示す図、第7図はピーク検出の様子を示す
図、第8図は生成ビークパタ゛−ンを示す図。 躬4図 茶5区

Claims (1)

    【特許請求の範囲】
  1. 入力生体信号を受け2次微分フィルタリング処理するフ
    ィルタ手段と、前記フィルタ手段の出力を受けピークを
    検出するピーク検出手段と、前記ピーク検出手段からの
    ピーク検出結果の時系列的変化により入力生体信号を所
    定の分類に類推する分類手段と、前記分類手段の出力を
    受け2次微分フィルタリング処理された信号を基に前記
    入力生体信号におけるピーク生起時刻及びピーク電圧値
    を計測する計測手段とを具備することを特徴とする生体
    信号処理方式。
JP59133705A 1984-06-28 1984-06-28 生体信号処理方式 Granted JPS6111017A (ja)

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 true JPS6111017A (ja) 1986-01-18
JPH0261250B2 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008168016A (ja) * 2007-01-15 2008-07-24 Fujifilm Corp 超音波診断装置、imt計測方法及びimt計測プログラム

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008168016A (ja) * 2007-01-15 2008-07-24 Fujifilm Corp 超音波診断装置、imt計測方法及びimt計測プログラム

Also Published As

Publication number Publication date
JPH0261250B2 (ja) 1990-12-19

Similar Documents

Publication Publication Date Title
CN106725450B (zh) Qrs波群检测方法
CN109381181B (zh) 心电信号特征波形的端点检测方法
US5772603A (en) Device for filtering ECG signals
CA1258099A (en) Egg with delta wave recognition by double
CN107303177B (zh) 一种心电图t波、p波的检测方法及***
JPS63226338A (ja) Ecg波信号からペース・パルス信号を除去する方法とecg波信号からペース・パルス信号を除去する装置
AU2019101756A4 (en) Method for analyzing heart rate variability signal based on extremum energy decomposition method
US4846568A (en) Method of analyzing an electroretinogram
CN107898443A (zh) 重搏波检测方法、装置及计算机存储介质
JPS6111017A (ja) 生体信号処理方式
EP0058050B1 (en) Measuring method
US5479933A (en) Method and apparatus for processing ECG signals
JPS644328Y2 (ja)
CN114795241A (zh) 一种心电信号的噪声检测方法及***
Kotas Projective filtering of time-aligned ECG beats for repolarization duration measurement
JPS58223927A (ja) デ−タ圧縮処理装置
JPH05212006A (ja) 心拍間隔計測装置
CN114624512A (zh) 调频信号反三角抑噪测频方法
Reddy et al. Analysis of interval measurements on CSE multilead reference ECGs
KR100458165B1 (ko) 심박동 변이 신호의 내삽 방법
JPH0238214B2 (ja)
JPS61206428A (ja) 生体信号ア−チフアクト検出除去方式
RU2636905C1 (ru) Способ аналого-цифрового измерения параметров при автоматической фрагментации электрокардиосигналов
WO2021049185A1 (ja) 情報処理システム、情報処理方法及びプログラム
CN107296598B (zh) 基于光电传感器的心率测定方法及装置