JPH04336035A - Respiratory waveform estimating method employing electrocardiogram - Google Patents

Respiratory waveform estimating method employing electrocardiogram

Info

Publication number
JPH04336035A
JPH04336035A JP3107073A JP10707391A JPH04336035A JP H04336035 A JPH04336035 A JP H04336035A JP 3107073 A JP3107073 A JP 3107073A JP 10707391 A JP10707391 A JP 10707391A JP H04336035 A JPH04336035 A JP H04336035A
Authority
JP
Japan
Prior art keywords
electrocardiogram
qrs
vector
wave
respiratory waveform
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
JP3107073A
Other languages
Japanese (ja)
Other versions
JP2820169B2 (en
Inventor
Hiromi Shimono
下野 太海
Mieko Osuga
美恵子 大須賀
Chie Akashi
明石 千恵
Hiromi Terashita
寺下 裕美
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP3107073A priority Critical patent/JP2820169B2/en
Publication of JPH04336035A publication Critical patent/JPH04336035A/en
Application granted granted Critical
Publication of JP2820169B2 publication Critical patent/JP2820169B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PURPOSE:To improve precision during estimation of a respiratory waveform from an electrocardiogram and to relieve a burden on a person to be examined and to simplify a measuring system by changing a measuring electrocardiogram from two channels to one channel. CONSTITUTION:Respiration is estimated by taking a change in an angle of a T-vector into consideration. A QRS vector and a T-vector are projected on one and the same axis to determine a projection axis for converting a change of an angle into a change of amplitude. The position of a breast part mounting electrode is decided along the projection axis and from the inside of the magnitude of the QRS wave and the T-wave of an electrocardiogram and the magnitude of an QRS wave or the magnitudes of the QRS wave and the T-wave, a respiratory curve is estimated.

Description

【発明の詳細な説明】[Detailed description of the invention]

【0001】0001

【産業上の利用分野】この発明は、呼吸を計測すること
なく、呼吸波形を心電図から推定する呼吸推定方法に関
するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a respiration estimation method for estimating a respiration waveform from an electrocardiogram without measuring respiration.

【0002】0002

【従来の技術】呼吸の計測は、伸縮により電気抵抗が変
化するセンサを被測定者の腹部または胸部に装着して行
う方法(例えば、特開昭62−362号公報、わかりや
すいME[産報出版(株)、1971年4月10日初版
発行]P.55〜56)が簡便であり、従来より良く用
いられている。この方法では前記被測定者の呼吸による
腹部または胸部の動きに比例して、前記センサの抵抗が
変化することを利用し、呼吸を電気信号として検出する
。前記電気信号は、高域通過フィルタを用いて不要な直
流成分と低周波成分を除去した後、増幅器で所定の信号
レベルに増幅される。その後ブラウン管オシロスコープ
等の表示装置により表示したり、データレコーダ等の記
録装置により記録する。このほかの例としては、検出手
段に温度センサを用い鼻孔付近に装着し、呼気と吸気の
温度差を検出することで呼吸を計測する方法(実開昭6
2−145603号公報、実開昭64−3603号公報
、及び実開平1−65009号公報)、前記温度センサ
と胸部インピーダンス検出を併用した呼吸モニタ装置(
実公平1−37949号公報)、呼気と吸気の気流によ
って音を発生させて呼吸を検出する方法(特開昭61−
79445号公報)等がある。
[Prior Art] Respiration is measured by attaching a sensor whose electrical resistance changes due to expansion and contraction to the abdomen or chest of the person to be measured (for example, Japanese Patent Application Laid-Open No. 62-362, easy-to-understand ME [Sanpo Publishing Co., Ltd., first edition published April 10, 1971] P.55-56) is simple and has been used more often than before. In this method, respiration is detected as an electrical signal by utilizing the fact that the resistance of the sensor changes in proportion to the movement of the abdomen or chest due to respiration of the subject. The electrical signal is amplified to a predetermined signal level by an amplifier after removing unnecessary DC components and low frequency components using a high-pass filter. Thereafter, it is displayed on a display device such as a cathode ray tube oscilloscope or recorded using a recording device such as a data recorder. Another example is a method of measuring respiration by using a temperature sensor as a detection means and attaching it near the nostrils to detect the temperature difference between exhaled and inhaled air.
No. 2-145603, Japanese Utility Model Application No. 64-3603, and Japanese Utility Model Application No. 1-65009), a respiration monitor device that uses the temperature sensor and chest impedance detection in combination (
Japanese Utility Model Publication No. 1-37949), a method for detecting breathing by generating sound by exhalation and inhalation airflow (Japanese Unexamined Patent Publication No. 1-37949),
79445), etc.

【0003】一方、心電図より呼吸を推定する方法とし
ては、第5回生体・生理工学シンポジウムで発表された
「心電図からの1回換気量推定」(同シンポジウム論文
集、P.277〜280)がある。この方法では図6に
示すように、心電図標準肢誘導の第I誘導及び第III
誘導の心電図を計測して、各々の心電図のQRSの大き
さを求め、61及び62にプロットしたものをベクトル
合成して得られる21のQRSベクトルを求め、このQ
RSベクトルと第I誘導軸との角度θの変化から呼吸波
形を推定するというものである。
[0003] On the other hand, as a method for estimating respiration from an electrocardiogram, "Estimation of tidal volume from an electrocardiogram" presented at the 5th Symposium on Biological and Physiological Engineering (Proceedings of the same symposium, pp. 277-280) was proposed. be. In this method, as shown in FIG.
Measure the electrocardiogram of the leads, find the QRS magnitude of each electrocardiogram, and combine the vectors plotted in 61 and 62 to obtain 21 QRS vectors.
The respiration waveform is estimated from the change in the angle θ between the RS vector and the I-th lead axis.

【0004】次に上記心電図から呼吸を推定する場合の
動作について説明する。図7はこの方法によるハードウ
ェアの構成図である。図7において、71は標準肢誘導
の第I誘導と第III誘導の心電図を増幅するための心
電図アンプ、72はA/D変換器、73はコンピュータ
である。71の増幅器で増幅された前記2種類の心電図
は各々1ミリ秒で72のA/D変換器でディジタル値に
変換され、73のコンピュータに取り込まれる。図8は
前記サンプルデータより呼吸波形を推定するためのアル
ゴリズムを示すフローチャートである。図8において、
まず81で前記サンプルデータに含まれる商用周波数成
分を除去するための前後10点の移動平均を行う。続い
て82のQRS決定部で、前記心電図サンプルデータの
1心拍毎のQ波、R波、S波のピーク時点を自動検出す
る。そして83のQRSベクトル算出部で、検出した前
記各ピーク時点における前記第I、第III誘導心電図
波形の大きさを求め、各々を第I誘導座標成分及び第I
II誘導座標成分ベクトルとして、両ベクトルを合成し
て、QRSベクトルの角度を算出する。算出した角度デ
ータは84のデータ補間部で3次スプライン補間法にて
補間し連続曲線を求め、85でこの連続曲線を推定呼吸
曲線として出力するものである。
Next, the operation for estimating respiration from the electrocardiogram will be explained. FIG. 7 is a block diagram of hardware according to this method. In FIG. 7, 71 is an electrocardiogram amplifier for amplifying the electrocardiograms of lead I and III of the standard limb leads, 72 is an A/D converter, and 73 is a computer. The two types of electrocardiograms amplified by the amplifier 71 are each converted into digital values by the A/D converter 72 in 1 millisecond, and then input into the computer 73. FIG. 8 is a flowchart showing an algorithm for estimating a respiratory waveform from the sample data. In Figure 8,
First, in step 81, a moving average of 10 points before and after is performed in order to remove commercial frequency components included in the sample data. Subsequently, the QRS determination unit 82 automatically detects the peak points of the Q wave, R wave, and S wave for each heartbeat of the electrocardiogram sample data. Then, the QRS vector calculation unit 83 calculates the magnitudes of the lead I and III electrocardiogram waveforms at each of the detected peak points, and calculates the magnitudes of the lead I and III electrocardiogram waveforms at each of the detected peak points,
Both vectors are combined as a II lead coordinate component vector to calculate the angle of the QRS vector. The calculated angle data is interpolated by the cubic spline interpolation method in a data interpolation unit 84 to obtain a continuous curve, and in 85 this continuous curve is outputted as an estimated respiration curve.

【0005】[0005]

【発明が解決しようとする課題】上記呼吸波形推定方法
では、心拍数と呼吸数の関係が(呼吸数÷心拍数)≧1
/2となるような場合、QRSベクトルの本来の連続変
動が再現されないため、呼吸波形推定に誤差を生じると
いう問題点がある。
[Problem to be solved by the invention] In the above respiratory waveform estimation method, the relationship between heart rate and breathing rate is (breathing rate ÷ heart rate) ≧1.
/2, the problem is that the original continuous fluctuation of the QRS vector is not reproduced, resulting in an error in respiratory waveform estimation.

【0006】また、従来から良く行われている呼吸計測
では、前述のごとく別途呼吸センサを被測定者に装着す
る必要があり被測定者の負荷が大きくなるという問題が
あった。従って心電図から呼吸波形を推定すれば、前記
呼吸センサの装着は不用となる。しかし上記の心電図か
ら呼吸波形を推定する方法においては、被測定者の四肢
に電極を装着するため、被測定者の動きは大きく制限さ
れるという問題点は解決されない。さらに標準肢誘導に
対応するように、胸部等で心電図を誘導して呼吸を推定
するようにして被測定者の四肢を解放しても、心電図を
2チャンネル計測しなければならず、心電図増幅器が二
つ必要となる。
[0006] Furthermore, in the respiration measurement that has been commonly carried out in the past, there has been a problem in that, as mentioned above, it is necessary to separately attach a respiration sensor to the person to be measured, which increases the burden on the person to be measured. Therefore, if the respiratory waveform is estimated from the electrocardiogram, there is no need to wear the respiratory sensor. However, in the above-mentioned method of estimating a respiratory waveform from an electrocardiogram, electrodes are attached to the limbs of the subject, so the problem that the movement of the subject is greatly restricted cannot be solved. Furthermore, even if the patient's limbs are released by guiding the electrocardiogram on the chest or the like to estimate respiration in order to correspond to standard limb guidance, the electrocardiogram must be measured in two channels, and the electrocardiogram amplifier is Two are required.

【0007】この発明は上記のような問題点を解決する
ためになされたもので、心拍数と呼吸数の関係が(呼吸
数÷心拍数)≧1/2となるような場合でも精度の高い
推定ができる呼吸推定方法を得ることを目的とする。
[0007] This invention was made to solve the above-mentioned problems, and is highly accurate even when the relationship between heart rate and breathing rate is (breathing rate ÷ heart rate) ≧ 1/2. The purpose is to obtain a respiration estimation method that can perform estimation.

【0008】また、呼吸波形を心電図から推定すること
により被測定者に呼吸センサを装着することを不用にし
、さらに被測定者の四肢の動きの制約を取り除くととも
に、計測する心電図を1チャンネルとしハードウェアの
簡略化を目的とする。
[0008] Furthermore, by estimating the respiratory waveform from the electrocardiogram, it is not necessary to attach a respiratory sensor to the subject, and furthermore, the restriction on the movement of the limbs of the subject is removed, and the electrocardiogram to be measured is made into one channel. The purpose is to simplify the software.

【0009】[0009]

【課題を解決するための手段】この発明に係る呼吸推定
方法は、QRSベクトルの角度だけでなくTベクトルの
角度も考慮して、前記両ベクトルの角度データを一つの
データファイルに合成して、補間するデータ数を2倍に
して呼吸曲線を推定することにより、前述の問題点であ
る(呼吸数÷心拍数)≧1/2となるような場合でも精
度の高い推定が出来るようにしたものである。
[Means for Solving the Problems] The respiration estimation method according to the present invention takes into account not only the angle of the QRS vector but also the angle of the T vector, and combines the angle data of both vectors into one data file. By estimating the respiration curve by doubling the number of data to be interpolated, highly accurate estimation can be made even in the case of (respiration rate ÷ heart rate) ≧ 1/2, which is the problem mentioned above. It is.

【0010】また、被測定者の四肢の動きの制約を取り
除くとともに計測する心電図を1チャンネルとするため
、前記QRSベクトルの角度変化が振幅変化となる射影
軸を決め、この射影軸に沿って前記被測定者の胸部に装
着する電極位置を決定し、胸部誘導心電図1チャンネル
だけを計測し、この心電図のQRSの大きさから呼吸波
形を推定するようにしたものである。
[0010] Furthermore, in order to remove restrictions on the movement of the extremities of the subject and to measure the electrocardiogram in one channel, a projection axis is determined in which the angular change in the QRS vector results in an amplitude change, and the projection axis is determined along this projection axis. The electrode position to be attached to the chest of the subject is determined, only one channel of the chest lead electrocardiogram is measured, and the respiratory waveform is estimated from the magnitude of the QRS of this electrocardiogram.

【0011】さらに、上記呼吸波形推定の精度を向上さ
せるため、QRSベクトルの角度変化及びTベクトルの
角度変化が各々QRS波、T波の振幅変化となる1つの
射影軸を決め、この射影軸から胸部に装着する電極装着
位置を決定して導出された心電図のQRS波の大きさと
T波の大きさから呼吸波形を推定するようにしたもので
ある。
Furthermore, in order to improve the accuracy of the above-mentioned respiratory waveform estimation, one projection axis is determined in which the angular change of the QRS vector and the angular change of the T vector correspond to the amplitude changes of the QRS wave and the T wave, respectively. The respiratory waveform is estimated from the QRS wave size and the T wave size of the electrocardiogram derived by determining the electrode attachment position to be attached to the chest.

【0012】0012

【作用】この発明における請求項1記載の呼吸推定方法
は、QRSベクトル及びTベクトルの両角度変化を用い
て呼吸波形を推定することにより、推定精度が向上する
In the respiration estimation method according to claim 1 of the present invention, estimation accuracy is improved by estimating the respiration waveform using both angular changes of the QRS vector and the T vector.

【0013】また、請求項2及び請求項3記載の呼吸推
定方法は、四肢に電極を装着する必要がないので、両腕
を動かすことができ被測定者の負荷軽減となるとともに
計測する心電図は1チャンネルだけなので、測定系の簡
略化が行える。
[0013] Furthermore, in the respiration estimation method according to claims 2 and 3, there is no need to attach electrodes to the limbs, so both arms can be moved, reducing the load on the person to be measured, and the electrocardiogram to be measured is Since there is only one channel, the measurement system can be simplified.

【0014】[0014]

【実施例】【Example】

実施例1.図1は請求項1記載の一実施例を示すフロー
チャートである。図1において11は従来例と同様に標
準肢誘導によって計測された2チャンネルの心電図デー
タである。12〜15においては従来例の図8の82及
び83のQRSベクトル角度の計算に加えTベクトルの
角度を計算するように構成されている。16では前記1
心拍毎の角度データを、QRS、T、QRS、T・・・
というように同一時間軸上に各々のピーク時点の時間関
係はそのままで生起順にデータを並べ、1つのデータ列
にまとめる。このときQRSベクトルとTベクトルの角
度変化の大きさが異なっている場合は、予め求めておい
た補正係数を用いてTベクトルの角度データを補正する
。17および18では、前記データ列に補間を行い連続
曲線に変換し、推定呼吸波形として出力する。
Example 1. FIG. 1 is a flowchart showing an embodiment according to claim 1. In FIG. 1, reference numeral 11 indicates two-channel electrocardiogram data measured by standard limb leads as in the conventional example. 12 to 15 are configured to calculate the angle of the T vector in addition to calculating the QRS vector angles 82 and 83 in FIG. 8 in the conventional example. In 16, the above 1
Angle data for each heartbeat, QRS, T, QRS, T...
In this way, the data are arranged on the same time axis in the order of occurrence, with the time relationships of the peak points unchanged, and are combined into one data string. At this time, if the magnitude of the angular change between the QRS vector and the T vector is different, the angle data of the T vector is corrected using a correction coefficient determined in advance. In steps 17 and 18, the data string is interpolated and converted into a continuous curve, which is output as an estimated respiratory waveform.

【0015】前記のような呼吸推定方法においては1心
拍について得られる角度データが2つであり、補間に使
われるデータはQRSベクトルのみの場合の2倍である
。従って、従来例に見られる心拍数と呼吸数の関係が(
呼吸数÷心拍数)≧1/2となるような場合でも精度の
高い推定ができることになる。
[0015] In the respiration estimation method as described above, two pieces of angle data are obtained for one heartbeat, and the data used for interpolation is twice that of the case where only the QRS vector is used. Therefore, the relationship between heart rate and respiration rate seen in the conventional example is (
Even in the case where (respiration rate ÷ heart rate)≧1/2, highly accurate estimation can be performed.

【0016】実施例2.図2はQRSベクトルの呼吸に
よる角度変化をQRSの大きさの単調変化として検出す
るための胸部に装着する電極位置の決め方の説明図であ
る。図において21は呼吸によって変化するQRSベク
トルを表わし、22は前記QRSベクトルのなす角度を
表わす。このとき心臓位置を通る23の破線で示される
方向に、前記QRSベクトルを射影すると24で示され
るように、前記角度差に対応する大きさの差となる。従
って、25で示される位置に電極を装着して胸部より心
電図を検出すれば、QRSの大きさの変化から呼吸波形
の推定ができる。
Example 2. FIG. 2 is an explanatory diagram of how to determine the position of an electrode to be attached to the chest in order to detect the angular change of the QRS vector due to respiration as a monotonous change in the QRS magnitude. In the figure, 21 represents the QRS vector that changes with respiration, and 22 represents the angle formed by the QRS vector. At this time, when the QRS vector is projected in the direction shown by the broken line 23 passing through the heart position, a difference in magnitude corresponding to the angular difference is obtained, as shown by 24. Therefore, by attaching an electrode to the position indicated by 25 and detecting the electrocardiogram from the chest, the respiratory waveform can be estimated from the change in QRS magnitude.

【0017】図3は上記の如く計測された心電図から呼
吸波形を推定するための一実施例のフローチャートであ
る。図において31は、計測した1チャンネルの心電図
をサンプルしたデータであり、32はこのデータから前
記実施例1同様にQRSピーク時点を検出する。33で
は32で求めたピーク時点のQRS波の大きさを求める
。QRSの大きさとしては図4に示すように、41の基
線Bを基準としたQ、R、S振幅の代数和即ちBQ+B
R+BS、または42〜44で示される各々の面積S1
 、S2 、S3 で算出されるS1 −S2 −S3
 、または積分して得られる平均値等を用いる。このよ
うにして求めたQRS波の大きさのデータ列を17にお
いて補間し連続波形に変換し、18で呼吸波形として出
力する。
FIG. 3 is a flowchart of an embodiment for estimating a respiratory waveform from an electrocardiogram measured as described above. In the figure, numeral 31 is data obtained by sampling the measured one-channel electrocardiogram, and numeral 32 is the data from which the QRS peak time is detected in the same manner as in the first embodiment. At step 33, the magnitude of the QRS wave at the peak point determined at step 32 is determined. As shown in Figure 4, the magnitude of the QRS is the algebraic sum of the Q, R, and S amplitudes based on the baseline B of 41, that is, BQ+B.
R+BS or each area S1 indicated by 42 to 44
, S2 , S3 calculated by S1 −S2 −S3
, or the average value obtained by integration. The data string of the QRS wave size obtained in this manner is interpolated at 17 and converted into a continuous waveform, and outputted at 18 as a respiratory waveform.

【0018】前記の呼吸推定方法では、従来例の呼吸推
定方法と推定精度は同じであるが、胸部より1チャンネ
ルの心電図を導出するだけで済むので測定系の簡略化が
できるとともに、被測定者の負担軽減が行える。
The respiration estimation method described above has the same estimation accuracy as the conventional respiration estimation method, but since it is only necessary to derive one channel of electrocardiogram from the chest, the measurement system can be simplified, and the measurement This can reduce the burden of

【0019】実施例3.図5は上記のQRS波の大きさ
から推定した場合より推定精度を上げるための一実施例
のフローチャートである。図において、51の心電図デ
ータは図3にTベクトルの方向も加えて決定した電極位
置による胸部誘導心電図のサンプルデータである。32
〜53はQRS波とT波の大きさを求める。T波の大き
さもQRS波同様、ピーク値、面積、平均値等が用いら
れる。54では前記実施例1と同様にT波の大きさの補
正を行い同一時間軸上に生起時間順にデータを並べて、
データ数2倍の補間データを作成する。17では前記Q
RS波の大きさの補間と同様にデータを補間し連続曲線
を求める。18では求めた前記連続曲線を推定呼吸波形
として出力する。
Example 3. FIG. 5 is a flowchart of an embodiment for improving the estimation accuracy from the above estimation based on the size of the QRS wave. In the figure, electrocardiogram data 51 is sample data of a thoracic lead electrocardiogram based on the electrode position determined by adding the direction of the T vector to FIG. 32
~53 determines the magnitude of the QRS wave and T wave. As with the QRS wave, the peak value, area, average value, etc. are used for the magnitude of the T wave. In step 54, the magnitude of the T wave is corrected in the same manner as in Example 1, and the data are arranged on the same time axis in order of occurrence time.
Create interpolated data with twice the number of data. 17, the above Q
A continuous curve is obtained by interpolating the data in the same way as interpolating the magnitude of the RS wave. In step 18, the obtained continuous curve is output as an estimated respiratory waveform.

【0020】[0020]

【発明の効果】この発明は、以上説明したような方法で
呼吸波形を推定するので、以下に記載される効果を奏す
る。
[Effects of the Invention] Since the present invention estimates a respiratory waveform using the method described above, it produces the effects described below.

【0021】1心拍について、QRSベクトルとTベク
トルの2つの角度データを用いることにより従来法の2
倍のデータが得られ、推定呼吸波形の精度が向上すると
ともに心拍数に対して呼吸速度が速い場合にも推定可能
となる。
[0021] For one heartbeat, by using two angle data of the QRS vector and the T vector, the conventional method 2
Double the data can be obtained, improving the accuracy of the estimated respiratory waveform and making it possible to estimate even when the respiratory rate is faster than the heart rate.

【0022】また、心ベクトルに対し異なる角度をなす
2種類の誘導による心電図を計測して得られるQRSベ
クトルの方向を基にして胸部に装着する電極位置を決定
し、この電極位置で計測される心電図のQRS波の大き
さの変化から呼吸波形を推定することにより、計測すべ
き心電図は1チャンネルで済む。さらに被測定者の四肢
に電極を付けなくても良いため、被測定者の負担軽減が
できる。
[0022] Furthermore, the electrode position to be attached to the chest is determined based on the direction of the QRS vector obtained by measuring electrocardiograms from two types of leads that are at different angles to the heart vector, and measurements are taken at this electrode position. By estimating the respiratory waveform from the change in the size of the QRS wave of the electrocardiogram, only one channel of the electrocardiogram is required. Furthermore, since there is no need to attach electrodes to the limbs of the person to be measured, the burden on the person to be measured can be reduced.

【0023】また、胸部誘導心電図1チャンネルのデー
タについても、T波の大きさを推定データに使用するこ
とにより推定精度の向上が図れる。
[0023] Furthermore, the estimation accuracy can be improved by using the magnitude of the T wave as the estimation data for data of one channel of the chest lead electrocardiogram.

【0024】また、ホルター心電計のように長時間にわ
たって心電図を記録し、検査や診断を行う場合がよくあ
るが、本発明を応用すれば、従来より行われている心電
図の解析に加え、推定呼吸波形の検査及び診断も行える
という利点がある。
[0024] Furthermore, electrocardiograms are often recorded over a long period of time using a Holter electrocardiograph for testing and diagnosis, but if the present invention is applied, in addition to the conventional electrocardiogram analysis, There is an advantage that the estimated respiratory waveform can also be inspected and diagnosed.

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

【図1】この発明の実施例1を示すフローチャートであ
る。
FIG. 1 is a flowchart showing a first embodiment of the present invention.

【図2】電極装着位置を決める説明図である。FIG. 2 is an explanatory diagram for determining the electrode attachment position.

【図3】この発明の実施例2を示すフローチャートであ
る。
FIG. 3 is a flowchart showing a second embodiment of the invention.

【図4】QRS波の大きさの説明図である。FIG. 4 is an explanatory diagram of the size of a QRS wave.

【図5】この発明の実施例3を示すフローチャートであ
る。
FIG. 5 is a flowchart showing a third embodiment of the present invention.

【図6】QRSベクトルの求め方の説明図である。FIG. 6 is an explanatory diagram of how to obtain a QRS vector.

【図7】従来の呼吸推定法のハードウェアの構成図であ
る。
FIG. 7 is a hardware configuration diagram of a conventional respiration estimation method.

【図8】従来の呼吸推定法を示すフローチャートである
FIG. 8 is a flowchart showing a conventional respiration estimation method.

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

11  標準肢誘導心電図2チャンネルのサンプルデー
タ12  Q波、R波、S波の各ピーク時点の検出13
  T波のピーク時点の検出 14  QRSベクトルの角度の計算 15  Tベクトルの角度の計算 16  QRSベクトル、Tベクトルの両角度データの
1つのデータ列に変換 17  データ補間 18  推定呼吸波形の出力 21  QRSベクトル 23  呼吸によるQRSベクトルの変化角度23  
射影軸 24  QRSの大きさ 25  電極 31  QRSベクトルを基に決めた電極位置で計測し
た心電図データ 32  QRS波のピーク検出 33  QRS波の大きさ算出 41  心電図基線 42  Q波の面積 43  R波の面積 44  S波の面積 51  QRSとTベクトルから決めた電極位置で計測
した心電図データ 52  T波ピークの検出 53  T波の大きさの算出 54  補間データの作成 61  第I誘導心電図のQRS波の大きさ62  第
III誘導心電図のQRS波の大きさ71  心電図増
幅器 72  A/D変換器 73  コンピュータ 81  2チャンネル心電図データの移動平均82  
QRS波決定部 83  QRSベクトル算出部 84  3次スプライン補間 85  推定呼吸波形出力
11 Sample data of standard limb lead electrocardiogram 2 channels 12 Detection of each peak point of Q wave, R wave, and S wave 13
Detection of peak time of T wave 14 Calculation of angle of QRS vector 15 Calculation of angle of T vector 16 Converting both angle data of QRS vector and T vector into one data string 17 Data interpolation 18 Output of estimated respiratory waveform 21 QRS vector 23 Angle of change in QRS vector due to breathing 23
Projection axis 24 QRS size 25 Electrode 31 Electrocardiogram data measured at electrode positions determined based on the QRS vector 32 QRS wave peak detection 33 QRS wave size calculation 41 ECG baseline 42 Q wave area 43 R wave area 44 Area of S wave 51 Electrocardiogram data measured at the electrode position determined from QRS and T vector 52 Detection of T wave peak 53 Calculation of T wave size 54 Creation of interpolated data 61 Size of QRS wave in lead I electrocardiogram 62 Size of QRS wave in lead III electrocardiogram 71 Electrocardiogram amplifier 72 A/D converter 73 Computer 81 Moving average of 2-channel electrocardiogram data 82
QRS wave determination unit 83 QRS vector calculation unit 84 Cubic spline interpolation 85 Estimated respiratory waveform output

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】  心ベクトルに対し異なる角度をなす2
種類の誘導による心電図を計測し、前記2種類の心電図
のQRS波よりQRSベクトルを、T波よりTベクトル
を各々求め、前記両ベクトルの各々の角度変化を用いて
呼吸波形を推定することを特徴とする心電図を用いた呼
吸波形推定方法。
[Claim 1] 2 that make different angles with respect to the center vector
The method is characterized in that electrocardiograms from different types of leads are measured, a QRS vector is obtained from the QRS waves of the two types of electrocardiograms, and a T vector is obtained from the T waves of the two types of electrocardiograms, and a respiratory waveform is estimated using the angular change of each of the two vectors. A method for estimating respiratory waveforms using electrocardiograms.
【請求項2】  請求項1記載の2種類の心電図による
QRSベクトルの方向を求め、前記ベクトルの方向を基
に電極の胸部装着位置を決定し、前記電極位置で計測し
た1チャンネルの心電図のQRS波の大きさより呼吸波
形を推定することを特徴とする心電図を用いた呼吸波形
推定方法。
2. Determine the direction of the QRS vector from the two types of electrocardiograms according to claim 1, determine the chest mounting position of the electrode based on the direction of the vector, and calculate the QRS of the 1-channel electrocardiogram measured at the electrode position. A method for estimating a respiratory waveform using an electrocardiogram, which is characterized by estimating a respiratory waveform based on the size of the wave.
【請求項3】  胸部誘導電極の装着位置をTベクトル
の方向も考慮して決定し、呼吸波形の推定もT波の大き
さを考慮して行うことを特徴とする請求項2記載の心電
図を用いた呼吸波形推定方法。
3. The electrocardiogram according to claim 2, wherein the mounting position of the chest lead electrode is determined by taking into account the direction of the T vector, and the estimation of the respiratory waveform is also performed by taking into account the size of the T wave. The respiratory waveform estimation method used.
JP3107073A 1991-05-13 1991-05-13 Respiratory waveform estimation method using electrocardiogram Expired - Fee Related JP2820169B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP3107073A JP2820169B2 (en) 1991-05-13 1991-05-13 Respiratory waveform estimation method using electrocardiogram

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3107073A JP2820169B2 (en) 1991-05-13 1991-05-13 Respiratory waveform estimation method using electrocardiogram

Publications (2)

Publication Number Publication Date
JPH04336035A true JPH04336035A (en) 1992-11-24
JP2820169B2 JP2820169B2 (en) 1998-11-05

Family

ID=14449801

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3107073A Expired - Fee Related JP2820169B2 (en) 1991-05-13 1991-05-13 Respiratory waveform estimation method using electrocardiogram

Country Status (1)

Country Link
JP (1) JP2820169B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016010616A (en) * 2014-06-30 2016-01-21 株式会社東芝 Respiration state estimation apparatus, respiration state estimation method, and program
JP2018519967A (en) * 2015-07-21 2018-07-26 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and system for automatic labeling of activities in ECG data
CN111603166A (en) * 2020-06-01 2020-09-01 河南华南医电科技有限公司 Respiratory signal detection method based on vector cardiogram

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5632570B1 (en) 2013-03-22 2014-11-26 パナソニック株式会社 Biological signal measurement system, apparatus, method, and program thereof

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016010616A (en) * 2014-06-30 2016-01-21 株式会社東芝 Respiration state estimation apparatus, respiration state estimation method, and program
JP2018519967A (en) * 2015-07-21 2018-07-26 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and system for automatic labeling of activities in ECG data
CN111603166A (en) * 2020-06-01 2020-09-01 河南华南医电科技有限公司 Respiratory signal detection method based on vector cardiogram

Also Published As

Publication number Publication date
JP2820169B2 (en) 1998-11-05

Similar Documents

Publication Publication Date Title
US10617325B2 (en) Respiratory rate measurement using a combination of respiration signals
US5025784A (en) Apparatus and method for detecting and processing impedance rheogram
US4945917A (en) Method and apparatus for indicating repetition intervals of a specified component of a composite electrical signal, particularly useful for displaying fetal r-waves
Fenton et al. Automated spectral characterization of wheezing in asthmatic children
US6241683B1 (en) Phonospirometry for non-invasive monitoring of respiration
US7343198B2 (en) System, software, and method for detection of sleep-disordered breathing using an electrocardiogram
JP2000508194A (en) Non-invasive monitoring of hemodynamic parameters using impedance cardiography
JP2003523248A (en) Non-invasive determination of cardiac output, pulmonary blood flow and blood gas volume
CA2262236C (en) Phonospirometry for non-invasive monitoring of respiration
Wang et al. An advanced signal processing technique for impedance cardiography
US5913308A (en) Apparatus and method for determining respiratory effort from muscle tremor information in ECG signals
CN104523262A (en) Sleep quality detection method based on electrocardiosignals
US20220330837A1 (en) Process and device for determining a respiratory and/or cardiogenic signal
Zhou et al. Philips QT interval measurement algorithms for diagnostic, ambulatory, and patient monitoring ECG applications
CENGIZ et al. The effect of ventilation on the accuracy of pulmonary artery and wedge pressure measurements
Bloch et al. Effect of mouthpiece breathing on cardiorespiratory response to intense exercise.
Tarassenko et al. Multi-sensor fusion for robust computation of breathing rate
JP2820169B2 (en) Respiratory waveform estimation method using electrocardiogram
Caggiano et al. Respiration derived from the electrocardiogram: a quantitative comparison of three different methods
JPH0323867B2 (en)
JP2509828B2 (en) Respiratory wave measuring device
US3480006A (en) Physiological volume-pressure diagram recording device
Weltman et al. Impedance pneumograph recording across the arms
Mitchell et al. Accurate, automated, continuously displayed pulmonary artery pressure measurement
JP4323277B2 (en) Sleep apnea type discrimination method

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees