JP2665161B2 - ECG analysis method - Google Patents
ECG analysis methodInfo
- Publication number
- JP2665161B2 JP2665161B2 JP6200382A JP20038294A JP2665161B2 JP 2665161 B2 JP2665161 B2 JP 2665161B2 JP 6200382 A JP6200382 A JP 6200382A JP 20038294 A JP20038294 A JP 20038294A JP 2665161 B2 JP2665161 B2 JP 2665161B2
- Authority
- JP
- Japan
- Prior art keywords
- artifacts
- wave
- morphological
- parameter
- heartbeat
- 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
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
【０００１】[0001]
【産業上の利用分野】本発明は心電図解析方法に関し、
特にアーチファクトによるＶＰＣ偽陽性率（検出すべき
ＶＰＣがないのにあるように表現される）の悪化を防止
する心電図解析方法に関する。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 a false positive rate of VPC (expressed as if there is no VPC to be detected) due to an artifact.
【０００２】[0002]
【従来の技術】解析処理を行う前に帯域フィルタ等でド
リフト等のノイズの除去を行い、その後ノイズ成分が減
少した心電図に関してＶＰＣ検出を行なっていた。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]
【発明が解決しようとする課題】しかし、従来の方法で
は１００％ノイズが除去された心電図を生成することは
できず、解析される波形の中には何らかの形でアーチフ
ァクトが存在しているために、アーチファクトを誤まっ
てＶＰＣと認めて、全体的なＶＰＣ検出の偽陽性率を低
下せしめるという欠点があった。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] It is an object of the present invention to provide an effect that even if an artifact is present in the waveform to be analyzed, the false positive rate of VPC detection is not reduced overall without erroneous artifacts being used as VPC. It is to provide an electrocardiogram analysis method with a certain method.
【０００５】[0005]
【課題を解決するための手段】前記目的を達成するた
め、本発明に係る心電図の解析方法は、解析対象となる
リアルタイム一心拍に対して、モフォロジパラメータを
計算する第一の計算部と、過去１０心拍分のモフォロジ
パラメータを記憶するリングバッファ記憶部と、過去１
０心拍のモフォロジパラメータの平均値を計算する第二
の計算部と、平均値との変化率を判定する第一の判定部
とを有している。In order to achieve the above object, an electrocardiogram analyzing method according to the present invention comprises a first calculating unit for calculating a morphological parameter for one real-time heartbeat to be analyzed; A ring buffer storage unit for storing morphological parameters for 10 heartbeats;
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.
【０００６】また、更に解析対象となる一心拍のＴ波の
形状を認識するモフォロジパラメータを算出する第三の
計算部と、そのモフォロジパラメータによりＴ波の存在
を認識する第二の判定部とを有している。Further, a third calculating unit for calculating a morphological parameter for recognizing the shape of a T-wave of one heartbeat to be analyzed and a second determining unit for recognizing the existence of a T-wave based on the morphological parameter. Have.
【０００７】[0007]
【実施例】以下、本発明の実施例を図面により説明す
る。Embodiments of the present invention will be described below with reference to the drawings.
【０００８】図１は、本発明の心電図解析方法を示す流
れ図である。FIG. 1 is a flowchart showing an electrocardiogram analyzing method according to the present invention.
【０００９】第一の計算部１は、各心拍に関するモフォ
ロジパラメータを計算する。モフォロジパラメータの算
出法は、Ｒ波のピークから±１００ｍＳの内部のデータ
に対し、周長と面積を求め面積を周長の２乗で割った商
をモフォロジパラメータとする。The first calculation unit 1 calculates morphological parameters for each heartbeat. The morphology parameter is calculated by calculating the perimeter and area of the data within ± 100 mS from the peak of the R wave and dividing the area by the square of the perimeter as the morphology parameter.
【００１０】リングバッファマネージャ部２は、第一の
計算部１で求められたモフォロジパラメータｍi を、最
新のリングバッファ記憶部８へ登録するとともに、ま
た、１０個分の心拍に対するモフォロジパラメータｍi
の内、最も古いデータを消去して、最新の心拍のモフォ
ロジパラメータにデータを追加する。[0010] ring buffer manager unit 2, the Mofo LOGIS parameters m i obtained by the first calculating portion 1, and registers to the latest ring buffer memory 8, also Mofo LOGIS parameters m i for 10 pieces of the heart
, The oldest data is deleted, and the data is added to the morphology parameter of the latest heartbeat.
【００１１】第二の計算部３は、過去１０心拍分のモフ
ォロジパラメータの平均値Ｍave を算出する。The second calculator 3 calculates an average value M ave of morphological parameters for the past 10 heartbeats.
【００１２】第一の判定部４は、解析対象となっている
心拍と、過去１０心拍分のモフォロジパラメータの平均
値とを比較して、両者の変化率がある一定の値を越える
か否かを判定する。The first determination unit 4 compares the heartbeat to be analyzed with the average value of the morphological parameters for the past ten heartbeats, and determines whether or not the rate of change between the two exceeds a certain value. Is determined.
【００１３】第三の計算部５は、解析対象となっている
心拍のＴ波部分のモフォロジパラメータを計算してい
る。Ｔ波のモフォロジパラメータの計算法は、先述した
Ｒ波のピークの後方２００ｍＳより後ろ向きにＴ波のピ
ークをサーチして、そのピークの前後±１００ｍＳの範
囲での波形形状のモフォロジパラメータを求める。The third calculating section 5 calculates morphological parameters of the T wave portion of the heartbeat to be analyzed. The method of calculating the morphological parameter of the T wave is to search for the peak of the T wave backward 200 mS after the peak of the R wave described above, and obtain the morphological parameter of the waveform shape within a range of ± 100 mS before and after the peak.
【００１４】第三の計算部５で求められたＴ波のモフォ
ロジパラメータが、実際のＴ波に近いか否かを判定する
のが第２の判定部６である。更に解析部７は、最終的に
解析対象となっている心拍が、正常波形かＶＰＣかを判
定する。The second determination unit 6 determines whether the morphological parameter of the T wave obtained by the third calculation unit 5 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.
【００１５】８は、リングバッファ記憶部である。Reference numeral 8 denotes a ring buffer storage unit.
【００１６】以下、本発明の実施例を図面により説明す
る。Hereinafter, embodiments of the present invention will be described with reference to the drawings.
【００１７】図１は、本発明の心電図解析方法を示す流
れ図である。FIG. 1 is a flowchart showing an electrocardiogram analyzing method according to the present invention.
【００１８】１は、各心拍に関するモフォロジパラメー
タＭi を計算する第一の計算部である。ここで云うモフ
ォロジパラメータは、以下の方法により計算される。[0018] 1 is the first calculation unit for calculating the Mofo Logistics parameters M i for each heart beat. The morphology parameter referred to here is calculated by the following method.
【００１９】まず、入力された一心拍の心電図波形に対
して、Ｒ波の頂点を基準として、前１００ｍＳのデータ
の開始点をｓｔａｒｔとし、ｓｔａｒｔ後１００ｍＳの
データの終了点をｅｎｄとする。その方法は、ｓｔａｒ
ｔ点とｅｎｄ点とを直線で結び、その直線の傾きを求め
αとする。次に、ｓｔａｒｔ点からのデータ数をｉとし
て、ｓｔａｒｔ点からｅｎｄ点の間のデータＤ０i に対
し、以下の演算を施してドリフトの補正後のデータＤi
とする（ここではｓｔａｒｔ点とｅｎｄ点を含めて補正
前はＤ０0 〜Ｄ０8 ，補正後はＤ0 〜Ｄ8 とする）。First, with respect to an input electrocardiogram waveform of one heartbeat, the start point of the data 100 ms before is set to start, and the end point of the data 100 mS after start is set to end based on the peak of the R wave. The method is star
The point t and the end point are connected by a straight line, and the slope of the straight line is determined as α. Next, assuming that the number of data from the start point is i, the data D0 i between the start point and the end point is subjected to the following operation to correct the data D i after 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).
【００２０】Ｄi ＝Ｄ０i −ｔａｎα×ｉ−Ｄ０0 次に、補正後のデータに対し、その面積Ａを以下の式に
より求める。D i = D 0 i −tan α × i−D 0 0 Next, the area A of the corrected data is determined by the following equation.
【００２１】 [0021]
【００２２】次に、補正後のデータに対し、その周長Ｌ
を以下の式により求める。Next, with respect to the corrected data, the circumference L
Is determined by the following equation.
【００２３】 [0023]
【００２４】次に、Ａ／Ｌ2 をモフォロジパラメータＭ
i とする。一般的に、このパラメータは２次元のパター
ン形状を定量的に表すことができることでよく知られて
いるパラメータである。Next, A / L 2 is represented by the morphological parameter M
i. Generally, this parameter is a well-known parameter that can quantitatively represent a two-dimensional pattern shape.
【００２５】以上の手順により求められたモフォロジパ
ラメータＭi を、リングバッファマネージャ２によりリ
ングバッファ記憶部８へ登録される。ここで云うリング
バッファは、最新の１０心拍分のモフォロジパラメータ
を格納することが可能で、リングバッファマネージャ２
により、もっとも古いパラメータの削除と最新のパラメ
ータの追加が同時に行われる。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. The ring buffer mentioned here can store the morphological parameters of the latest 10 heartbeats, and the ring buffer manager 2
As a result, the oldest parameter is deleted and the latest parameter is added at the same time.
【００２６】次に、第二の計算部３で、上記リングバッ
ファ内の最新１０心拍分のモフォロジパラメータの平均
値を算出し、平均値Ｍ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.
【００２７】次に、第一の判定部４で解析対象となって
いる心拍のモフォロジパラメータＭi を、上記平均値Ｍ
ave と比較して、両者の変化率が大きく、ある一定の閾
値Ｘを越える場合には第三の計算部５へ、また変化率が
小さく、ある一定の閾値Ｘを越えない場合には、解析部
７へとそれぞれ分岐される。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.
【００２８】第三の計算部５は、第一の計算部２でＲ波
に関するモフォロジパラメータを求めたのと同様な方法
で、Ｔ波に関するモフォロジパラメータＭ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. 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.
【００２９】第二の判定部６はモフォロジパラメータＭ
i がある一定の閾値Ｚを越えてＴ波らしい形状であると
判断された場合には、心電図のＶＰＣ波形として解析部
７に、また一定の閾値Ｚを越えずＴ波らしい形状でない
と判断された場合には、アーチファクトとして解析部７
には送られずに処理が終了される。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.
【００３０】このような一連の処理によって、アーチフ
ァクトが解析部７に入力されなくなることにより、アー
チファクトをＶＰＣと誤って判定しシステム全体の診断
精度の低下を抑制することができる。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]
【発明の効果】以上説明したように本発明は、過去１０
心拍分の波形に対するモフォロジパラメータのリングバ
ッファを有し、その平均値と比較し、更に変化率が大と
判定された心拍に対し、そのＴ波の存在を認識すること
でＶＰＣを誤まってアーチファクトと認識することが非
常に少ない確率で、入力波形から効果的にアーチファク
トだけを検出することができるという効果を有する。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.
【図１】本発明に係る心電図解析方法を示す流れ図FIG. 1 is a flowchart showing an electrocardiogram analysis method according to the present invention.
【図２】補正前の心電波形を示す波形図FIG. 2 is a waveform chart showing an electrocardiogram waveform before correction.
【図３】補正後の心電波形を示す波形図FIG. 3 is a waveform chart showing a corrected electrocardiogram waveform.
１ 第一の計算部 ２ リングバッファマネージャ ３ 第二の計算部 ４ 第一の判定部 ５ 第三の計算部 ６ 第二の判定部 ７ 解析部 ８ リングバッファ記憶部 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)
応する波形の形状を表すモフォロジパラメータを求める
第一の計算手段と、前記モフォロジパラメータを過去の
心拍の平均値と比較しアーチファクトを認識する第二の
計算手段と、アーチファクトを誤って前記心室性期外収
縮に判定することを回避する判断手段とを有することを
特徴とする心電図解析装置。 The present invention relates to detection of premature ventricular contractions.
Find morphological parameters representing the shape of the corresponding waveform
First calculating means, and the morphological parameter
Second to recognize artifacts by comparing with heart rate average
Calculation means and the artifact
Having judgment means for avoiding judgment
Characteristic ECG analyzer.
応する波形の形状を表すモフォロジパラメータを求める
第一の計算手段と、前記モフォロジパラメータを過去の
心拍の平均値と比較しアーチファクトを認識する第二の
計算手段と、認識されたアーチファクトからさらにモフ
ォロジパラメータによりＴ波の形状を認識する第三の計
算手段と、前記Ｔ波が存在しないと判断されたアーチフ
ァクトのみを最終的なアーチファクトとして認識し、ア
ーチファクトを誤って前記心室性期外収縮に判定するこ
とを回避する判断手段とを有することを特徴とする心電
図解析装置。 2. The method according to claim 1, wherein the detection of a premature ventricular contraction is performed with respect to the heartbeat.
Find morphological parameters representing the shape of the corresponding waveform
First calculating means, and the morphological parameter
Second to recognize artifacts by comparing with heart rate average
More morphs from the calculation means and the recognized artifacts
A third method for recognizing the shape of the T wave using
Calculating means, and an artifact determined to be absent of the T wave.
Only the artifacts as final artifacts,
Erroneous judgment of ventricular premature contraction
And a determination means for avoiding the
Diagram analyzer.
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 JPH0856914A (en) | 1996-03-05 |
| JP2665161B2 true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SE9604319D0 (en) * | 1996-11-25 | 1996-11-25 | Pacesetter Ab | Medical detecting system |
| CN1428130A (en) | 2001-12-27 | 2003-07-09 | 泷泽清 | Diagnostic method and device, diagnostic program |
| DK1790285T3 (en) | 2004-05-20 | 2013-05-21 | Kiyoshi Takizawa | Diagnostic method and apparatus for this and diagnostic program |
| US8543193B2 (en) * | 2005-05-13 | 2013-09-24 | Cardiocore Lab., Inc. | Method and apparatus for rapid interpretive analysis of electrocardiographic waveforms |
-
1994
- 1994-08-25 JP JP6200382A patent/JP2665161B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0856914A (en) | 1996-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6597943B2 (en) | Method of using spectral measures to distinguish among atrialfibrillation, atrial-flutter and other cardiac rhythms | |
| US8774908B2 (en) | Systems for detecting cardiac arrhythmias | |
| US9955890B2 (en) | Heart electrophysiological signal analysis system | |
| Zhu et al. | An R-peak detection method based on peaks of Shannon energy envelope | |
| US5628321A (en) | Processing velocity information in an ultrasonic system | |
| Hamilton et al. | Quantitative investigation of QRS detection rules using the MIT/BIH arrhythmia database | |
| Sadhukhan et al. | R-peak detection algorithm for ECG using double difference and RR interval processing | |
| WO2021169296A1 (en) | Method and apparatus for processing electrocardiogram data, computer device, and storage medium | |
| CN108324271B (en) | Electrocardiosignal identification method and system and electrocardiosignal monitoring equipment | |
| CN103110417A (en) | Automatic electrocardiogram recognition system | |
| US6668189B2 (en) | Method and system for measuring T-wave alternans by alignment of alternating median beats to a cubic spline | |
| US20230015562A1 (en) | Device and process for ecg measurements | |
| CN111358459A (en) | Arrhythmia identification method, device, equipment and storage medium | |
| CN111839494A (en) | Heart rate monitoring method and system | |
| Ghodrati et al. | RR interval analysis for detection of atrial fibrillation in ECG monitors | |
| JP2665161B2 (en) | ECG analysis method | |
| US20230309899A1 (en) | Determination method and determination apparatus for beginning of t-wave, storage medium and computer program product | |
| CN102302373A (en) | Method for detecting heart sound signal feature points based on multiplexing multi-resolution wavelet transformation | |
| CN109199361B (en) | Method, storage medium and device for extracting R peak time in electrocardiosignal data | |
| CN114795241B (en) | Electrocardiosignal noise detection method and system | |
| CN115590524A (en) | QT interval prolonging identification method and system based on convolutional neural network | |
| CN110477904B (en) | House flapping signal detection method and device | |
| CN112244867A (en) | Electrocardiosignal analysis algorithm and system for atrial fibrillation detection | |
| JP3224569B2 (en) | Cardiac arrhythmia detection device based on ECG signal | |
| JPS6156633A (en) | Apparatus for detecting abnormality of electrocardiograph |
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 |