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
electrocardiogram
mofo
wave
artifacts
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
正樹 金澤
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
Application status is Expired - Lifetime legal-status Critical

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]

【産業上の利用分野】本発明は心電図解析方法に関し、 The present invention relates to an electrocardiogram analysis method,
特にアーチファクトによるVPC偽陽性率の悪化を防止する心電図解析方法に関する。 And more particularly to an electrocardiogram analysis method for preventing the deterioration of the VPC false positive rate due to artifacts.

【0002】 [0002]

【従来の技術】解析処理を行う前に帯域フィルタ等でドリフト等のノイズの除去を行い、その後ノイズ成分が減少した心電図に関してVPC検出を行なっていた。 Performs removal of noise such as drift in band filter or the like before performing the Related Art analyzing process, then the noise component had performed VPC detected for electrocardiogram decreased.

【0003】 [0003]

【発明が解決しようとする課題】しかし、従来の方法では100%ノイズが除去された心電図を生成することはできず、解析される波形の中には何らかの形でアーチファクトが存在しているために、アーチファクトを誤まってVPCと認めて、全体的なVPC検出の偽陽性率を低下せしめるという欠点があった。 [SUMMARY OF THE INVENTION] However, it is impossible to generate an electrocardiogram of 100% noise has been removed by conventional methods, some of the waveform to be analyzed in due to the presence of artifacts in some form , admits VPC waiting erroneous artifacts, has a drawback that allowed to reduce the false-positive rate of the overall VPC detection.

【0004】本発明の目的は、前述の解析される波形の中にアーチファクトが存在しても、アーチファクトを誤まってVPCとすることなく、全体的なVPC検出の偽陽性率を低下させないという効果のある心電図解析方法を提供することにある。 An object of the present invention, effect that even if there is an artifact in the waveform to be analyzed in the foregoing, without the VPC waiting erroneous artifacts, not to reduce the false positive rate of the overall VPC detection It is to provide an electrocardiogram analysis method with.

【0005】 [0005]

【課題を解決するための手段】前記目的を達成するため、本発明に係る心電図解析方法は、解析対象となるリアルタイムな一心拍に対して、モフォロジパラメータを計算する第一の計算部と、過去10心拍分のモフォロジパラメータを記憶するリングバッファ記憶部と、過去1 To achieve the above object, according to an aspect of, electrocardiogram analysis method according to the present invention is to provide real-time one heartbeat to be analyzed, a first calculation unit for calculating the Mofo Logistics parameters, past a ring buffer memory section for storing Mofo LOGIS parameters of 10 heart beats, the past 1
0心拍のモフォロジパラメータの平均値を計算する第二の計算部と、平均値との変化率を判定する第一の判定部とを有している。 0 and a second calculation unit for calculating an average value of Mofo LOGIS parameters heart rate, and a first determination unit that determines the rate of change of the average value.

【0006】また、更に解析対象となる一心拍のT波の形状を認識するモフォロジパラメータを算出する第三の計算部と、そのモフォロジパラメータによりT波の存在を認識する第二の判定部とを有している。 Moreover, a third calculator for further calculating the Mofo Logistics parameter recognize T-wave shape of one heartbeat to be analyzed, and a second determination unit recognizes the presence of the T-wave by its Mofo Logistics parameter It has.

【0007】 [0007]

【実施例】以下、本発明の実施例を図面により説明する。 BRIEF DESCRIPTION by examples of the present invention with reference to the accompanying drawings.

【0008】図1は、本発明の心電図解析方法を示す流れ図である。 [0008] Figure 1 is a flow diagram illustrating an electrocardiogram analysis method of the present invention.

【0009】第一の計算部1は、各心拍に関するモフォロジパラメータを計算する。 [0009] The first calculation unit 1 calculates the Mofo Logistics parameters for each heart beat. モフォロジパラメータの算出法は、R波のピークから±100mSの内部のデータに対し、周長と面積を求め面積を周長の2乗で割った商をモフォロジパラメータをする。 Calculation of Mofo LOGIS parameters to the internal data of ± 100 mS from the peak of R-wave, the quotient obtained by dividing by the square of the area sought perimeter and area perimeter to the Mofo Logistics parameter.

【0010】リングバンファマネージャ部2は、第一の計算部1で求められたモフォロジパラメータm iを、最新のリングバッファ記憶部8へ登録するとともに、また、10個分の心拍に対するモフォロジパラメータm i [0010] Ring vans file manager portion 2, a Mofo LOGIS parameters m i obtained by the first calculating portion 1, and registers to the latest ring buffer memory 8, also Mofo LOGIS parameter m for 10 pieces of the heart i
の内、最も古いデータを消去して、最新の心拍のモフォロジパラメータにデータを追加する。 Of, to erase the oldest data, to add data to the Mofo Logis parameters of the latest heart rate.

【0011】第二の計算部3は、過去10心拍分のモフォロジパラメータの平均値M aveを算出する。 [0011] The second calculation unit 3 calculates the average value M ave of Mofo LOGIS parameters for the past 10 heartbeats.

【0012】第一の判定部4は、解析対象となっている心拍と、過去10心拍分のモフォロジパラメータの平均値とを比較して、両者の変化率がある一定の値を越えるか否かを判定する。 [0012] The first determination unit 4, and heart that is analyzed, whether by comparing the average value of Mofo LOGIS parameters for the past 10 heart beats, exceeds a certain value both change rate the judges.

【0013】第三の計算部5は、解析対象となっている心拍のT波部分のモフォロジパラメータを計算している。 [0013] The third calculation portion 5, and calculates the Mofo LOGIS parameters of the T-wave portion of the heart rate has been analyzed. T波のモフォロジパラメータの計算法は、R波のピークの後方200mSより後ろ向きにT波のピークをサーチして、そのピークの前後±100mSの範囲での波形形状のモフォロジパラメータを求める。 Calculation of Mofo LOGIS parameters of the T-wave is to search for the peak of the T wave rearwardly from the rear 200mS peak of R-wave to determine the Mofo LOGIS parameters of the waveform shape in the range of about ± 100 mS of the peak.

【0014】第三の計算部3で求められたT波のモフォロジパラメータが、実際のT波に近いか否かを判定するのが第2の判定部6である。 [0014] Mofo Logistics parameter of the third calculation portion 3 at the obtained T-wave, to determine whether similar to the actual T wave is a second determination unit 6. 更に解析部7は、最終的に解析対象となっている心拍が、正常波形かVPCかを判定する。 Further analysis unit 7, the heart rate which is finally analyzed is determined whether normal waveform or VPC.

【0015】8は、リングバッファ記憶部である。 [0015] 8 is a ring buffer storage unit.

【0016】以下、本発明の実施例を図面により説明する。 [0016] Hereinafter, will be explained by examples of the present invention with reference to the accompanying drawings.

【0017】図1は、本発明の心電図解析方法を示す流れ図である。 [0017] Figure 1 is a flow diagram illustrating an electrocardiogram analysis method of the present invention.

【0018】1は、各心拍に関するモフォロジパラメータM iを計算する第一の計算部である。 [0018] 1 is the first calculation unit for calculating the Mofo Logistics parameters M i for each heart beat. ここで云うモフォロジパラメータは、以下の方法により計算される。 Mofo Logistics parameter referred to herein is calculated by the following method.

【0019】まず、入力された一心拍の心電図波形に対して、R波の頂点を基準として、前100mSのデータを開始点をstartとし、start後100mSのデータの終了点をendとする。 [0019] First, the input one heartbeat electrocardiogram wave, based on the apex of the R-wave, the data before 100mS the starting point and start, the end point of data of start after 100mS and end The. その方法は、star The method, star
t点とend点とを直線で結び、その直線の傾きを求めαとする。 A t point and end point connected by straight lines, and α determine the slope of the straight line. 次に、start点からのデータ数をiとして、start点からend点の間のデータD0 iに対し、以下の演算を施してドリフトの補正後のデータD i Then, the number of data from the start point as a i, for the data D0 i between start point of the end point, the data D i of the corrected drift by performing the following operations
とする(ここではstart点とend点を含めて補正前はD0 0 〜D0 8 ,補正後はD 0 〜D 8とする)。 To (here before correction, including start point and end point is D0 0 ~D0 8, after the correction to D 0 ~D 8).

【0020】D i =D0 i −tanα×i−D0 0次に、補正後のデータに対し、その面積Aを以下の式により求める。 [0020] Next D i = D0 i -tanα × i -D0 0, with respect to data after correction, obtain the area A by the following equation.

【0021】 [0021]

【0022】次に、補正後のデータに対し、その周長L [0022] Then, for the data after the correction, the circumferential length L
を以下の式により求める。 The calculated by the following equation.

【0023】 [0023]

【0024】次に、A/L 2をモフォロジパラメータM Next, Mofo the A / L 2 Logis parameters M
i とする。 And i. 一般的に、このパラメータは2次元のパターン形状を定量的に表すことができることでよく知られているパラメータである。 Typically, this parameter is a parameter well known that it can be quantitatively representing a two-dimensional pattern shape.

【0025】以上の手順により求められたモフォロジパラメータM iを、リングバッファマネージャ2によりリングバッファ記憶部8へ登録される。 The more Mofo LOGIS parameters M i obtained by the procedure, is registered by the ring buffer manager 2 to the ring buffer memory 8. ここで云うリングバッファは、最新の10心拍分のモフォロジパラメータを格納することが可能で、リングバッファマネージャ2 Ring buffer referred to herein are capable to store Mofo LOGIS parameters of the latest 10 heart beats, the ring buffer manager 2
により、もっとも古いパラメータの削除と最新のパラメータの追加が同時に行われる。 By, addition of the oldest parameters deleted and the latest parameters is carried out at the same time.

【0026】次に、第二の計算部3で、上記リングバッファ内の最新10心拍分のモフォロジパラメータの平均値を算出し、平均値M aveとする。 Next, in the second calculation portion 3 calculates an average value of Mofo LOGIS parameters of the most recent 10 heart beats in the ring buffer, and the average value M ave.

【0027】次に、第一の判定部4で解析対象となっている心拍のモフォロジパラメータM iを、上記平均値M Next, the Mofo LOGIS parameters M i heart rate that is analyzed by the first determination unit 4, the average value M
aveと比較して、両者の変化率が大きく、ある一定の閾値Xを越える場合には第三の計算部5へ、また変化率が小さく、ある一定の閾値Xを越えない場合には、解析部7へとそれぞれ分岐される。 Compared to ave, both the rate of change is large and if it exceeds a certain threshold X is the third calculation unit 5, also small rate of change, if not exceeding a certain threshold X, the analysis They are respectively branched into parts 7.

【0028】第三の計算部5は、第一の計算部2でR波に関するモフォロジパラメータを求めたのと同様な方法で、T波に関するモフォロジパラメータM iを計算する。 [0028] The third calculation portion 5, the first calculator 2 in a manner similar to that sought Mofo Logistics parameters for R-wave method, calculates the Mofo Logistics parameters M i for T-wave. そして、R波の頂点の後方200mSより後ろ向きにT波の頂点を検出し、その頂点を基準にして前後10 Then, to detect the apex of the T-wave rearwardly from the rear 200mS vertices of R wave, longitudinal with respect to the its apex 10
0mSの波形データに関してモフォロジパラメータを計算する。 Calculating the Mofo Logistics parameters for waveform data of 0 mS.

【0029】第二の判定部6はモフォロジパラメータM [0029] The second determination unit 6 Mofo Logistics parameter M
iがある一定の閾値Zを越えてT波らしい形状であると判断された場合には、心電図のVPC波形として解析部7に、また一定の閾値Zを越えずT波らしい形状でないと判断された場合には、アーチファクトとして解析部7 when i is determined to be T-wave likely shape beyond a certain threshold Z which is, it is determined that the analysis unit 7 as VPC waveform of an electrocardiogram, also not T-wave seems shape does not exceed a certain threshold Z If the, the analysis unit 7 as an artifact
には送られずに処理が終了される。 The process is terminated without sent to.

【0030】このような一連の処理によって、アーチファクトが解析部7に入力されなくなることにより、アーチファクトをVPCと誤って判定しシステム全体の診断精度の低下を抑制することができる。 [0030] Such a series of processing, by the artifact is not input to the analyzing unit 7, it is possible to suppress the reduction in the diagnostic accuracy of the overall system is determined erroneously artifacts and VPC.

【0031】 [0031]

【発明の効果】以上説明したように本発明は、過去10 The present invention described above, according to the present invention, the past 10
心拍分の波形に対するモフォロジパラメータのリングバッファを有し、その平均値と比較し、更に変化率が大と判定された心拍に対し、そのT波の存在を認識することでVPCを誤まってアーチファクトと認識することが非常に少ない確率で、入力波形から効果的にアーチファクトだけを検出することができるという効果を有する。 Has a ring buffer Mofo LOGIS parameters for heart beats of the waveform is compared with the average value, further to heart rate change rate is determined to be large, erroneous and VPC by recognizing the presence of the T wave artifact with very little probability to recognize and has the effect that it is possible to detect only effectively artifacts from the input waveform.

【図面の簡単な説明】 BRIEF DESCRIPTION OF THE DRAWINGS

【図1】本発明に係る心電図解析方法を示す流れ図 Flow diagram illustrating an electrocardiogram analysis method according to the invention, FIG

【図2】補正前の心電波形を示す波形図 FIG. 2 is a waveform diagram showing a correction before the electrocardiographic waveform

【図3】補正後の心電波形を示す波形図 FIG. 3 is a waveform diagram showing the electrocardiographic waveform after correction

【符号の説明】 DESCRIPTION OF SYMBOLS

1 第一の計算部 2 リングバッファマネージャ 3 第二の計算部 4 第一の判定部 5 第三の計算部 6 第二の判定部 7 解析部 8 リングバッファ記憶部 1 The first computing unit 2 ring buffer manager 3 the second calculator 4 first determination unit 5 third calculator 6 second determination unit 7 analyzer 8 ring buffer memory unit

Claims (2)

    【特許請求の範囲】 [The claims]
  1. 【請求項1】 心電図の解析処理において、心室性期外収縮(以下「VPC」という)を検出する際に、波形の形状を表すモフォロジパラメータを過去の心拍の平均値と比較することによりアーチファクトを認識し、これによりVPC検出の誤りを抑制して成ることを特徴とする心電図解析方法。 1. A analysis of the electrocardiogram, in detecting premature ventricular contractions (hereinafter referred to as "VPC"), the artifacts by comparing Mofo Logistics parameter representing the shape of the waveform and the average value of past heart recognized, electrocardiogram analysis method characterized by comprising thereby suppressing an error in the VPC detection.
  2. 【請求項2】 心電図の解析処理において、VPCを検出する際に、波形の形状を表すモフォロジパラメータを過去の心拍の平均値と比較することによりアーチファクトを認識し、前記VPCを誤まってアーチファクトにすることを回避して成るT波形状を認識し、T波が存在しないと判断された場合にのみ最終的にアーチファクトを認識しVPC検出の偽陽性率を抑圧して成ることを特徴とする心電図解析方法。 2. A analysis of the electrocardiogram, when detecting the VPC, the Mofo LOGIS parameter representing the shape of the waveform recognize artifacts by comparing the average value of past heart rate, the artifacts erroneous the VPC recognizes avoided to made T wave shape that, electrocardiogram, characterized in that formed by suppressing the false positive rate of only finally recognizes artifacts VPC detected when it is determined that the T-wave is absent analysis method.
JP6200382A 1994-08-25 1994-08-25 Electrocardiogram 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 Electrocardiogram analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP6200382A JP2665161B2 (en) 1994-08-25 1994-08-25 Electrocardiogram 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 Electrocardiogram 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 analyzing electrocardiogram wave speed
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 analyzing electrocardiogram wave speed

Also Published As

Publication number Publication date
JP2665161B2 (en) 1997-10-22

Similar Documents

Publication Publication Date Title
Saini et al. QRS detection using K-Nearest Neighbor algorithm (KNN) and evaluation on standard ECG databases
EP1450898B1 (en) Dual chamber method and apparatus for diagnosis and treatment of arrhythmias
JP4149598B2 (en) Method of automatically setting the collimator x-ray imaging system during image acquisition and x-ray imaging system,
JP4745346B2 (en) Signal template in implantable cardiac device
US5792065A (en) Method and apparatus for determining T-wave marker points during QT dispersion analysis
US6937887B2 (en) Atrial fibrillation detection method and apparatus
EP0853285B1 (en) Method and system for detecting acute myocardial infarction
US8433398B2 (en) Signal analysis system for heart condition determination
Gritzali et al. Detection of the P and T waves in an ECG
EP0570896B1 (en) Method for tachycardia discrimination or rate adaptive pacing with multiple sensors
JP4352315B2 (en) Signal processing method / apparatus and a pulse photometer using the same
US5817027A (en) Method and apparatus for classifying heartbeats in an ECG waveform
US7272437B2 (en) Systems for processing electrocardiac signals having superimposed complexes
US5464020A (en) Diagnosing and treating subacute cardiac dysfunction
US5271411A (en) Method and apparatus for ECG signal analysis and cardiac arrhythmia detection
US20040077941A1 (en) Method and system for image improvement with ECG gating and dose reduction in CT imaging
US20070219453A1 (en) Automated analysis of a cardiac signal based on dynamical characteristics of the cardiac signal
Afonso et al. ECG beat detection using filter banks
JP3578766B2 (en) Cardiac monitoring diagnostic system
US6684100B1 (en) Curvature based method for selecting features from an electrophysiologic signals for purpose of complex identification and classification
Caiani et al. Warped-average template technique to track on a cycle-by-cycle basis the cardiac filling phases on left ventricular volume
EP1178855B1 (en) Monitoring apparatus using wavelet transforms for the analysis of heart rhythms
US8346349B2 (en) Method and apparatus for predicting patient outcomes from a physiological segmentable patient signal
US6312382B1 (en) Method and apparatus for extracting cardiac information from acoustic information acquired with an ultrasound device
EP0776630A1 (en) Method and apparatus for obtaining heartbeat measurements in an ECG waveform

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