JPH0856914A - System for analyzing electrocardiogram - Google Patents

System for analyzing electrocardiogram

Info

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
Application number
JP6200382A
Other languages
Japanese (ja)
Other versions
JP2665161B2 (en
Inventor
Masaki Kanazawa
正樹 金澤
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.)
TOCHIGI NIPPON DENKI KK
Original Assignee
TOCHIGI NIPPON DENKI KK
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 TOCHIGI NIPPON DENKI KK filed Critical TOCHIGI NIPPON DENKI KK
Priority to JP6200382A priority Critical patent/JP2665161B2/en
Publication of JPH0856914A publication Critical patent/JPH0856914A/en
Application granted granted Critical
Publication of JP2665161B2 publication Critical patent/JP2665161B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

PURPOSE: To effectively detect an artifact, to reduce a pseudo positive ratio in VPC detection and to improve global analysis accuracy in the VPC detecting processing of electrocardiogram analysis. CONSTITUTION: This system is provided with a first calculation part 1 which calculates a morphological parameter for each heartbeat, a ring buffer manager 2 which registers past ten heartbeats on a ring buffer 8, a second calculation part 3 which calculates the mean value of the heartbeats in the ring buffer 8, a first decision part 4 which decides the possibility of the artifact from a change rate with the mean value, a third calculation part 5 which calculates the morphological parameter of a T wave, a second decision part 6 which decides he existence of the T wave from a calculated numeric value, and an analysis part 7 which analyzes an electrocardiogram from which the artifact is removed.

Description

【発明の詳細な説明】Detailed Description of the Invention

【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.

【図面の簡単な説明】[Brief description of drawings]

【図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.

【符号の説明】[Explanation of symbols]

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)

【特許請求の範囲】[Claims] 【請求項1】 心電図の解析処理において、心室性期外
収縮(以下「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.
【請求項2】 心電図の解析処理において、VPCを検
出する際に、波形の形状を表すモフォロジパラメータを
過去の心拍の平均値と比較することによりアーチファク
トを認識し、前記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.
JP6200382A 1994-08-25 1994-08-25 ECG analysis method Expired - Lifetime JP2665161B2 (en)

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)

* Cited by examiner, † Cited by third party
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

Cited By (4)

* Cited by examiner, † Cited by third party
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