JPH0856914A - System for analyzing electrocardiogram - Google Patents
System for analyzing electrocardiogramInfo
- Publication number
- JPH0856914A JPH0856914A JP6200382A JP20038294A JPH0856914A JP H0856914 A JPH0856914 A JP H0856914A JP 6200382 A JP6200382 A JP 6200382A JP 20038294 A JP20038294 A JP 20038294A JP H0856914 A JPH0856914 A JP H0856914A
- Authority
- JP
- Japan
- Prior art keywords
- vpc
- artifact
- electrocardiogram
- wave
- ring buffer
- 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
Abstract
Description
【0001】[0001]
【産業上の利用分野】本発明は心電図解析方法に関し、
特にアーチファクトによるVPC偽陽性率の悪化を防止
する心電図解析方法に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an electrocardiogram analyzing method,
In particular, the present invention relates to an electrocardiogram analysis method for preventing deterioration of the VPC false positive rate due to artifacts.
【0002】[0002]
【従来の技術】解析処理を行う前に帯域フィルタ等でド
リフト等のノイズの除去を行い、その後ノイズ成分が減
少した心電図に関してVPC検出を行なっていた。2. Description of the Related Art Before analysis, noise such as drift is removed by a band-pass filter or the like, and thereafter, VPC detection is performed on an electrocardiogram in which noise components are reduced.
【0003】[0003]
【発明が解決しようとする課題】しかし、従来の方法で
は100%ノイズが除去された心電図を生成することは
できず、解析される波形の中には何らかの形でアーチフ
ァクトが存在しているために、アーチファクトを誤まっ
てVPCと認めて、全体的なVPC検出の偽陽性率を低
下せしめるという欠点があった。However, the conventional method cannot generate an electrocardiogram from which 100% noise has been removed, and some artifacts exist in the analyzed waveform. However, there is a drawback that an artifact is mistakenly recognized as VPC and the overall false positive rate of VPC detection is reduced.
【0004】本発明の目的は、前述の解析される波形の
中にアーチファクトが存在しても、アーチファクトを誤
まってVPCとすることなく、全体的なVPC検出の偽
陽性率を低下させないという効果のある心電図解析方法
を提供することにある。The object of the present invention is to prevent the false positive rate of VPC detection as a whole from being lowered even if there is an artifact in the waveform to be analyzed, without causing the artifact to be erroneously regarded as VPC. It is to provide a certain electrocardiogram analysis method.
【0005】[0005]
【課題を解決するための手段】前記目的を達成するた
め、本発明に係る心電図解析方法は、解析対象となるリ
アルタイムな一心拍に対して、モフォロジパラメータを
計算する第一の計算部と、過去10心拍分のモフォロジ
パラメータを記憶するリングバッファ記憶部と、過去1
0心拍のモフォロジパラメータの平均値を計算する第二
の計算部と、平均値との変化率を判定する第一の判定部
とを有している。In order to achieve the above object, an electrocardiogram analysis method according to the present invention includes a first calculation unit for calculating a morphology parameter for one real-time heartbeat to be analyzed, and a past calculation unit. A ring buffer storage unit that stores morphological parameters for 10 heartbeats, and a past 1
It has a second calculation unit that calculates the average value of the morphology parameter of zero heartbeat, and a first determination unit that determines the rate of change from the average value.
【0006】また、更に解析対象となる一心拍のT波の
形状を認識するモフォロジパラメータを算出する第三の
計算部と、そのモフォロジパラメータによりT波の存在
を認識する第二の判定部とを有している。Further, a third calculation unit for calculating a morphology parameter for recognizing the shape of the T wave of one heartbeat to be analyzed and a second determination unit for recognizing the existence of the T wave by the morphology parameter. Have
【0007】[0007]
【実施例】以下、本発明の実施例を図面により説明す
る。Embodiments of the present invention will be described below with reference to the drawings.
【0008】図1は、本発明の心電図解析方法を示す流
れ図である。FIG. 1 is a flow chart showing the electrocardiogram analysis method of the present invention.
【0009】第一の計算部1は、各心拍に関するモフォ
ロジパラメータを計算する。モフォロジパラメータの算
出法は、R波のピークから±100mSの内部のデータ
に対し、周長と面積を求め面積を周長の2乗で割った商
をモフォロジパラメータをする。The first calculation unit 1 calculates the morphology parameter for each heartbeat. The calculation method of the morphology parameter is to calculate the perimeter and the area with respect to the data inside ± 100 mS from the peak of the R wave, and use the quotient obtained by dividing the area by the square of the perimeter as the morphology parameter.
【0010】リングバンファマネージャ部2は、第一の
計算部1で求められたモフォロジパラメータmi を、最
新のリングバッファ記憶部8へ登録するとともに、ま
た、10個分の心拍に対するモフォロジパラメータmi
の内、最も古いデータを消去して、最新の心拍のモフォ
ロジパラメータにデータを追加する。The ring banfer manager unit 2 registers the morphology parameter m i obtained by the first calculation unit 1 in the latest ring buffer storage unit 8 and also the morphology parameter m for 10 heartbeats. i
Delete the oldest data and add the data to the latest heart beat morphology parameters.
【0011】第二の計算部3は、過去10心拍分のモフ
ォロジパラメータの平均値Mave を算出する。The second calculator 3 calculates the average value M ave of the morphology parameters for the past 10 heartbeats.
【0012】第一の判定部4は、解析対象となっている
心拍と、過去10心拍分のモフォロジパラメータの平均
値とを比較して、両者の変化率がある一定の値を越える
か否かを判定する。The first judging unit 4 compares the heartbeat to be analyzed with the average value of the morphological parameters for the past 10 heartbeats, and determines whether the rate of change between the two exceeds a certain value. To judge.
【0013】第三の計算部5は、解析対象となっている
心拍のT波部分のモフォロジパラメータを計算してい
る。T波のモフォロジパラメータの計算法は、R波のピ
ークの後方200mSより後ろ向きにT波のピークをサ
ーチして、そのピークの前後±100mSの範囲での波
形形状のモフォロジパラメータを求める。The third calculating section 5 calculates morphological parameters of the T wave portion of the heartbeat to be analyzed. The calculation method of the morphology parameter of the T wave is to search the peak of the T wave backward 200 mS behind the peak of the R wave, and obtain the morphology parameter of the waveform shape within a range of ± 100 mS before and after the peak.
【0014】第三の計算部3で求められたT波のモフォ
ロジパラメータが、実際のT波に近いか否かを判定する
のが第2の判定部6である。更に解析部7は、最終的に
解析対象となっている心拍が、正常波形かVPCかを判
定する。The second determining unit 6 determines whether or not the morphology parameter of the T wave obtained by the third calculating unit 3 is close to the actual T wave. Further, the analysis unit 7 determines whether the heartbeat finally analyzed is a normal waveform or a VPC.
【0015】8は、リングバッファ記憶部である。Reference numeral 8 is a ring buffer storage unit.
【0016】以下、本発明の実施例を図面により説明す
る。Hereinafter, embodiments of the present invention will be described with reference to the drawings.
【0017】図1は、本発明の心電図解析方法を示す流
れ図である。FIG. 1 is a flowchart showing an electrocardiogram analyzing method according to the present invention.
【0018】1は、各心拍に関するモフォロジパラメー
タMi を計算する第一の計算部である。ここで云うモフ
ォロジパラメータは、以下の方法により計算される。Reference numeral 1 is a first calculation unit for calculating a morphological parameter M i for each heartbeat. The morphology parameter referred to here is calculated by the following method.
【0019】まず、入力された一心拍の心電図波形に対
して、R波の頂点を基準として、前100mSのデータ
を開始点をstartとし、start後100mSの
データの終了点をendとする。その方法は、star
t点とend点とを直線で結び、その直線の傾きを求め
αとする。次に、start点からのデータ数をiとし
て、start点からend点の間のデータD0i に対
し、以下の演算を施してドリフトの補正後のデータDi
とする(ここではstart点とend点を含めて補正
前はD00 〜D08 ,補正後はD0 〜D8 とする)。First, with respect to the input electrocardiogram waveform of one heartbeat, with the vertex of the R wave as a reference, the start point is the data of 100 mS before and the end point of the data of 100 mS after the start. The method is star
The point t and the end point are connected by a straight line, and the inclination of the straight line is obtained and set as α. Next, assuming that the number of data from the start point is i, the following calculation is performed on the data D0 i between the start point and the end point to obtain the data D i after the drift correction.
To (here before correction, including start point and end point is D0 0 ~D0 8, after the correction to D 0 ~D 8).
【0020】Di =D0i −tanα×i−D00 次に、補正後のデータに対し、その面積Aを以下の式に
より求める。D i = D0 i −tan α × i −D0 0 Next, the area A of the corrected data is calculated by the following equation.
【0021】 [0021]
【0022】次に、補正後のデータに対し、その周長L
を以下の式により求める。Next, with respect to the corrected data, the circumference L
Is calculated by the following formula.
【0023】 [0023]
【0024】次に、A/L2 をモフォロジパラメータM
i とする。一般的に、このパラメータは2次元のパター
ン形状を定量的に表すことができることでよく知られて
いるパラメータである。Next, A / L 2 is set to the morphology parameter M
i. Generally, this parameter is a well-known parameter that can quantitatively represent a two-dimensional pattern shape.
【0025】以上の手順により求められたモフォロジパ
ラメータMi を、リングバッファマネージャ2によりリ
ングバッファ記憶部8へ登録される。ここで云うリング
バッファは、最新の10心拍分のモフォロジパラメータ
を格納することが可能で、リングバッファマネージャ2
により、もっとも古いパラメータの削除と最新のパラメ
ータの追加が同時に行われる。The morphology parameter M i obtained by the above procedure is registered in the ring buffer storage unit 8 by the ring buffer manager 2. The ring buffer referred to here is capable of storing the latest 10 morphological parameters for 10 heartbeats.
Will delete the oldest parameter and add the latest one at the same time.
【0026】次に、第二の計算部3で、上記リングバッ
ファ内の最新10心拍分のモフォロジパラメータの平均
値を算出し、平均値Mave とする。Next, the second calculator 3 calculates the average value of the morphological parameters for the latest 10 heartbeats in the ring buffer, and sets it as the average value M ave .
【0027】次に、第一の判定部4で解析対象となって
いる心拍のモフォロジパラメータMi を、上記平均値M
ave と比較して、両者の変化率が大きく、ある一定の閾
値Xを越える場合には第三の計算部5へ、また変化率が
小さく、ある一定の閾値Xを越えない場合には、解析部
7へとそれぞれ分岐される。Next, the morphological parameter M i of the heartbeat to be analyzed by the first determination unit 4 is calculated by the average value M
Compared to ave , if the rate of change of both is greater and exceeds a certain threshold X, it goes to the third calculation unit 5. If the rate of change is small and does not exceed a certain threshold X, it is analyzed. Each branch is made to a unit 7.
【0028】第三の計算部5は、第一の計算部2でR波
に関するモフォロジパラメータを求めたのと同様な方法
で、T波に関するモフォロジパラメータMi を計算す
る。そして、R波の頂点の後方200mSより後ろ向き
にT波の頂点を検出し、その頂点を基準にして前後10
0mSの波形データに関してモフォロジパラメータを計
算する。The third calculation unit 5 calculates the morphology parameter M i for the T wave by the same method as that used for obtaining the morphology parameter for the R wave in the first calculation unit 2. Then, the apex of the T wave is detected backward 200 mS behind the apex of the R wave, and a 10
A morphology parameter is calculated for the 0 ms waveform data.
【0029】第二の判定部6はモフォロジパラメータM
i がある一定の閾値Zを越えてT波らしい形状であると
判断された場合には、心電図のVPC波形として解析部
7に、また一定の閾値Zを越えずT波らしい形状でない
と判断された場合には、アーチファクトとして解析部7
には送られずに処理が終了される。The second determination unit 6 determines the morphological parameter M
If i exceeds a certain threshold Z and is determined to be a T-wave-like shape, the analysis unit 7 determines that the waveform does not exceed a certain threshold Z and is not a T-wave-like shape as a VPC waveform of an electrocardiogram. If an error occurs, the analysis unit 7
The processing is terminated without being sent to.
【0030】このような一連の処理によって、アーチフ
ァクトが解析部7に入力されなくなることにより、アー
チファクトをVPCと誤って判定しシステム全体の診断
精度の低下を抑制することができる。By such a series of processing, artifacts are not input to the analysis unit 7, so that artifacts are erroneously determined as VPCs, and a decrease in diagnostic accuracy of the entire system can be suppressed.
【0031】[0031]
【発明の効果】以上説明したように本発明は、過去10
心拍分の波形に対するモフォロジパラメータのリングバ
ッファを有し、その平均値と比較し、更に変化率が大と
判定された心拍に対し、そのT波の存在を認識すること
でVPCを誤まってアーチファクトと認識することが非
常に少ない確率で、入力波形から効果的にアーチファク
トだけを検出することができるという効果を有する。As described above, the present invention has
It has a ring buffer of morphological parameters for the waveform of the heartbeat, compares it with its average value, and recognizes the presence of the T-wave for the heartbeat whose rate of change is determined to be large. It is possible to effectively detect only the artifact from the input waveform with a very low probability of recognizing that the input signal is an artifact.
【図1】本発明に係る心電図解析方法を示す流れ図FIG. 1 is a flowchart showing an electrocardiogram analysis method according to the present invention.
【図2】補正前の心電波形を示す波形図FIG. 2 is a waveform diagram showing an electrocardiographic waveform before correction.
【図3】補正後の心電波形を示す波形図FIG. 3 is a waveform chart showing a corrected electrocardiogram waveform.
1 第一の計算部 2 リングバッファマネージャ 3 第二の計算部 4 第一の判定部 5 第三の計算部 6 第二の判定部 7 解析部 8 リングバッファ記憶部 DESCRIPTION OF SYMBOLS 1 First calculation part 2 Ring buffer manager 3 Second calculation part 4 First judgment part 5 Third calculation part 6 Second judgment part 7 Analysis part 8 Ring buffer storage part
Claims (2)
収縮(以下「VPC」という)を検出する際に、波形の
形状を表すモフォロジパラメータを過去の心拍の平均値
と比較することによりアーチファクトを認識し、これに
よりVPC検出の誤りを抑制して成ることを特徴とする
心電図解析方法。1. When detecting a ventricular extrasystole (hereinafter, referred to as "VPC") in an electrocardiographic analysis process, an artifact is detected by comparing a morphological parameter representing a waveform shape with an average value of past heartbeats. A method for analyzing an electrocardiogram, which comprises recognizing and thereby suppressing an error in VPC detection.
出する際に、波形の形状を表すモフォロジパラメータを
過去の心拍の平均値と比較することによりアーチファク
トを認識し、前記VPCを誤まってアーチファクトにす
ることを回避して成るT波形状を認識し、T波が存在し
ないと判断された場合にのみ最終的にアーチファクトを
認識しVPC検出の偽陽性率を抑圧して成ることを特徴
とする心電図解析方法。2. In the electrocardiogram analysis process, when a VPC is detected, an artifact is recognized by comparing a morphological parameter representing a waveform shape with an average value of past heartbeats, and the VPC is erroneously identified as an artifact. An electrocardiogram characterized by recognizing a T-wave shape by avoiding the above, finally recognizing an artifact only when it is determined that a T-wave does not exist, and suppressing a false positive rate of VPC detection. analysis method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6200382A JP2665161B2 (en) | 1994-08-25 | 1994-08-25 | ECG analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6200382A JP2665161B2 (en) | 1994-08-25 | 1994-08-25 | ECG analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0856914A true JPH0856914A (en) | 1996-03-05 |
JP2665161B2 JP2665161B2 (en) | 1997-10-22 |
Family
ID=16423396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP6200382A Expired - Lifetime JP2665161B2 (en) | 1994-08-25 | 1994-08-25 | ECG analysis method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2665161B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5865760A (en) * | 1996-11-25 | 1999-02-02 | Pacesetter Ab | System for detecting changes in body posture |
EP1323376A2 (en) | 2001-12-27 | 2003-07-02 | Kiyoshi Takizawa | Method and apparatus for ECG diagnosis |
JP2008539988A (en) * | 2005-05-13 | 2008-11-20 | カーディオコア ラブ、インコーポレイテッド | Method and apparatus for high-speed ECG waveform analysis |
US8036733B2 (en) | 2004-05-20 | 2011-10-11 | Kiyoshi Takizawa | Diagnostic parameter calculation method, system for diagnosis and diagnostic program |
-
1994
- 1994-08-25 JP JP6200382A patent/JP2665161B2/en not_active Expired - Lifetime
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5865760A (en) * | 1996-11-25 | 1999-02-02 | Pacesetter Ab | System for detecting changes in body posture |
EP1323376A2 (en) | 2001-12-27 | 2003-07-02 | Kiyoshi Takizawa | Method and apparatus for ECG diagnosis |
US8036733B2 (en) | 2004-05-20 | 2011-10-11 | Kiyoshi Takizawa | Diagnostic parameter calculation method, system for diagnosis and diagnostic program |
JP2008539988A (en) * | 2005-05-13 | 2008-11-20 | カーディオコア ラブ、インコーポレイテッド | Method and apparatus for high-speed ECG waveform analysis |
Also Published As
Publication number | Publication date |
---|---|
JP2665161B2 (en) | 1997-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5628321A (en) | Processing velocity information in an ultrasonic system | |
CN110013247B (en) | Detection, differentiation and positioning method for electrocardiogram P wave and T wave | |
Zhu et al. | An R-peak detection method based on peaks of Shannon energy envelope | |
US6597943B2 (en) | Method of using spectral measures to distinguish among atrialfibrillation, atrial-flutter and other cardiac rhythms | |
US4721114A (en) | Method of detecting P-waves in ECG recordings | |
JP3787336B2 (en) | Heart rate detection apparatus using PPG and method thereof | |
EP0776630B1 (en) | Apparatus for obtaining heartbeat measurements in an ECG waveform | |
US20090275849A1 (en) | Methods for Detection of Cardiac Arrhythmias | |
US9314180B2 (en) | Heart electrophysiological signal analysis system | |
US20070066907A1 (en) | Method and device for analyzing a periodic or semi-periodic signal | |
CN110742599B (en) | Electrocardiosignal feature extraction and classification method and system | |
CN111481192B (en) | Electrocardiosignal R wave detection method based on improved U-Net | |
JP2635079B2 (en) | Method of removing pace pulse signal from ECG wave signal and apparatus for removing pace pulse signal from ECG wave signal | |
US6668189B2 (en) | Method and system for measuring T-wave alternans by alignment of alternating median beats to a cubic spline | |
CN111513706A (en) | Method and device for detecting electrocardiosignals containing abnormal R waves | |
US20230309899A1 (en) | Determination method and determination apparatus for beginning of t-wave, storage medium and computer program product | |
JP2665161B2 (en) | ECG analysis method | |
CN102302373A (en) | Method for detecting heart sound signal feature points based on multiplexing multi-resolution wavelet transformation | |
CN114795241B (en) | Electrocardiosignal noise detection method and system | |
JP2004243131A (en) | Method and system for improved measurement of t-wave alternant | |
CN112244867B (en) | Electrocardiosignal analysis method and system for atrial fibrillation detection | |
KR102367639B1 (en) | Processing method for electrocardiogram signal | |
Fedotov et al. | A QRS-complex detector of the electrocardiogram signal for the long-term monitoring of the patient’s condition | |
CN110477904B (en) | House flapping signal detection method and device | |
KR101293248B1 (en) | System for discriminating premature ventricular contraction using RR interval, method for discriminating premature ventricular contraction and the recording media storing the program performing the said method |
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: 19970603 |