JP2820169B2 - 心電図を用いた呼吸波形推定方法 - Google Patents
心電図を用いた呼吸波形推定方法Info
- Publication number
- JP2820169B2 JP2820169B2 JP3107073A JP10707391A JP2820169B2 JP 2820169 B2 JP2820169 B2 JP 2820169B2 JP 3107073 A JP3107073 A JP 3107073A JP 10707391 A JP10707391 A JP 10707391A JP 2820169 B2 JP2820169 B2 JP 2820169B2
- Authority
- JP
- Japan
- Prior art keywords
- electrocardiogram
- vector
- qrs
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
【０００１】
【産業上の利用分野】この発明は、呼吸を計測すること
なく、呼吸波形を心電図から推定する呼吸推定方法に関
するものである。
なく、呼吸波形を心電図から推定する呼吸推定方法に関
するものである。
【０００２】
【従来の技術】呼吸の計測は、伸縮により電気抵抗が変
化するセンサを被測定者の腹部または胸部に装着して行
う方法（例えば、特開昭６２−３６２号公報、わかりや
すいＭＥ［産報出版（株）、１９７１年４月１０日初版
発行］Ｐ．５５〜５６）が簡便であり、従来より良く用
いられている。この方法では前記被測定者の呼吸による
腹部または胸部の動きに比例して、前記センサの抵抗が
変化することを利用し、呼吸を電気信号として検出す
る。前記電気信号は、高域通過フィルタを用いて不要な
直流成分と低周波成分を除去した後、増幅器で所定の信
号レベルに増幅される。その後ブラウン管オシロスコー
プ等の表示装置により表示したり、データレコーダ等の
記録装置により記録する。このほかの例としては、検出
手段に温度センサを用い鼻孔付近に装着し、呼気と吸気
の温度差を検出することで呼吸を計測する方法（実開昭
６２−１４５６０３号公報、実開昭６４−３６０３号公
報、及び実開平１−６５００９号公報）、前記温度セン
サと胸部インピーダンス検出を併用した呼吸モニタ装置
（実公平１−３７９４９号公報）、呼気と吸気の気流に
よって音を発生させて呼吸を検出する方法（特開昭６１
−７９４４５号公報）等がある。
化するセンサを被測定者の腹部または胸部に装着して行
う方法（例えば、特開昭６２−３６２号公報、わかりや
すいＭＥ［産報出版（株）、１９７１年４月１０日初版
発行］Ｐ．５５〜５６）が簡便であり、従来より良く用
いられている。この方法では前記被測定者の呼吸による
腹部または胸部の動きに比例して、前記センサの抵抗が
変化することを利用し、呼吸を電気信号として検出す
る。前記電気信号は、高域通過フィルタを用いて不要な
直流成分と低周波成分を除去した後、増幅器で所定の信
号レベルに増幅される。その後ブラウン管オシロスコー
プ等の表示装置により表示したり、データレコーダ等の
記録装置により記録する。このほかの例としては、検出
手段に温度センサを用い鼻孔付近に装着し、呼気と吸気
の温度差を検出することで呼吸を計測する方法（実開昭
６２−１４５６０３号公報、実開昭６４−３６０３号公
報、及び実開平１−６５００９号公報）、前記温度セン
サと胸部インピーダンス検出を併用した呼吸モニタ装置
（実公平１−３７９４９号公報）、呼気と吸気の気流に
よって音を発生させて呼吸を検出する方法（特開昭６１
−７９４４５号公報）等がある。
【０００３】一方、心電図より呼吸を推定する方法とし
ては、第５回生体・生理工学シンポジウムで発表された
「心電図からの１回換気量推定」（同シンポジウム論文
集、Ｐ．２７７〜２８０）がある。この方法では図６に
示すように、心電図標準肢誘導の第Ｉ誘導及び第ＩＩＩ
誘導の心電図を計測して、各々の心電図のＱＲＳの大き
さを求め、６１及び６２にプロットしたものをベクトル
合成して得られる２１のＱＲＳベクトルを求め、このＱ
ＲＳベクトルと第Ｉ誘導軸との角度θの変化から呼吸波
形を推定するというものである。
ては、第５回生体・生理工学シンポジウムで発表された
「心電図からの１回換気量推定」（同シンポジウム論文
集、Ｐ．２７７〜２８０）がある。この方法では図６に
示すように、心電図標準肢誘導の第Ｉ誘導及び第ＩＩＩ
誘導の心電図を計測して、各々の心電図のＱＲＳの大き
さを求め、６１及び６２にプロットしたものをベクトル
合成して得られる２１のＱＲＳベクトルを求め、このＱ
ＲＳベクトルと第Ｉ誘導軸との角度θの変化から呼吸波
形を推定するというものである。
【０００４】次に上記心電図から呼吸を推定する場合の
動作について説明する。図７はこの方法によるハードウ
ェアの構成図である。図７において、７１は標準肢誘導
の第Ｉ誘導と第ＩＩＩ誘導の心電図を増幅するための心
電図アンプ、７２はＡ／Ｄ変換器、７３はコンピュータ
である。７１の増幅器で増幅された前記２種類の心電図
は各々１ミリ秒で７２のＡ／Ｄ変換器でディジタル値に
変換され、７３のコンピュータに取り込まれる。図８は
前記サンプルデータより呼吸波形を推定するためのアル
ゴリズムを示すフローチャートである。図８において、
まず８１で前記サンプルデータに含まれる商用周波数成
分を除去するための前後１０点の移動平均を行う。続い
て８２のＱＲＳ決定部で、前記心電図サンプルデータの
１心拍毎のＱ波、Ｒ波、Ｓ波のピーク時点を自動検出す
る。そして８３のＱＲＳベクトル算出部で、検出した前
記各ピーク時点における前記第Ｉ、第ＩＩＩ誘導心電図
波形の大きさを求め、各々を第Ｉ誘導座標成分及び第Ｉ
ＩＩ誘導座標成分ベクトルとして、両ベクトルを合成し
て、ＱＲＳベクトルの角度を算出する。算出した角度デ
ータは８４のデータ補間部で３次スプライン補間法にて
補間し連続曲線を求め、８５でこの連続曲線を推定呼吸
曲線として出力するものである。
動作について説明する。図７はこの方法によるハードウ
ェアの構成図である。図７において、７１は標準肢誘導
の第Ｉ誘導と第ＩＩＩ誘導の心電図を増幅するための心
電図アンプ、７２はＡ／Ｄ変換器、７３はコンピュータ
である。７１の増幅器で増幅された前記２種類の心電図
は各々１ミリ秒で７２のＡ／Ｄ変換器でディジタル値に
変換され、７３のコンピュータに取り込まれる。図８は
前記サンプルデータより呼吸波形を推定するためのアル
ゴリズムを示すフローチャートである。図８において、
まず８１で前記サンプルデータに含まれる商用周波数成
分を除去するための前後１０点の移動平均を行う。続い
て８２のＱＲＳ決定部で、前記心電図サンプルデータの
１心拍毎のＱ波、Ｒ波、Ｓ波のピーク時点を自動検出す
る。そして８３のＱＲＳベクトル算出部で、検出した前
記各ピーク時点における前記第Ｉ、第ＩＩＩ誘導心電図
波形の大きさを求め、各々を第Ｉ誘導座標成分及び第Ｉ
ＩＩ誘導座標成分ベクトルとして、両ベクトルを合成し
て、ＱＲＳベクトルの角度を算出する。算出した角度デ
ータは８４のデータ補間部で３次スプライン補間法にて
補間し連続曲線を求め、８５でこの連続曲線を推定呼吸
曲線として出力するものである。
【０００５】
【発明が解決しようとする課題】上記呼吸波形推定方法
では、心拍数と呼吸数の関係が（呼吸数÷心拍数）≧１
／２となるような場合、ＱＲＳベクトルの本来の連続変
動が再現されないため、呼吸波形推定に誤差を生じると
いう問題点がある。
では、心拍数と呼吸数の関係が（呼吸数÷心拍数）≧１
／２となるような場合、ＱＲＳベクトルの本来の連続変
動が再現されないため、呼吸波形推定に誤差を生じると
いう問題点がある。
【０００６】また、従来から良く行われている呼吸計測
では、前述のごとく別途呼吸センサを被測定者に装着す
る必要があり被測定者の負荷が大きくなるという問題が
あった。従って心電図から呼吸波形を推定すれば、前記
呼吸センサの装着は不用となる。しかし上記の心電図か
ら呼吸波形を推定する方法においては、被測定者の四肢
に電極を装着するため、被測定者の動きは大きく制限さ
れるという問題点は解決されない。さらに標準肢誘導に
対応するように、胸部等で心電図を誘導して呼吸を推定
するようにして被測定者の四肢を解放しても、心電図を
２チャンネル計測しなければならず、心電図増幅器が二
つ必要となる。
では、前述のごとく別途呼吸センサを被測定者に装着す
る必要があり被測定者の負荷が大きくなるという問題が
あった。従って心電図から呼吸波形を推定すれば、前記
呼吸センサの装着は不用となる。しかし上記の心電図か
ら呼吸波形を推定する方法においては、被測定者の四肢
に電極を装着するため、被測定者の動きは大きく制限さ
れるという問題点は解決されない。さらに標準肢誘導に
対応するように、胸部等で心電図を誘導して呼吸を推定
するようにして被測定者の四肢を解放しても、心電図を
２チャンネル計測しなければならず、心電図増幅器が二
つ必要となる。
【０００７】この発明は上記のような問題点を解決する
ためになされたもので、心拍数と呼吸数の関係が（呼吸
数÷心拍数）≧１／２となるような場合でも精度の高い
推定ができる呼吸推定方法を得ることを目的とする。
ためになされたもので、心拍数と呼吸数の関係が（呼吸
数÷心拍数）≧１／２となるような場合でも精度の高い
推定ができる呼吸推定方法を得ることを目的とする。
【０００８】また、呼吸波形を心電図から推定すること
により被測定者に呼吸センサを装着することを不用に
し、さらに被測定者の四肢の動きの制約を取り除くとと
もに、計測する心電図を１チャンネルとしハードウェア
の簡略化を目的とする。
により被測定者に呼吸センサを装着することを不用に
し、さらに被測定者の四肢の動きの制約を取り除くとと
もに、計測する心電図を１チャンネルとしハードウェア
の簡略化を目的とする。
【０００９】
【課題を解決するための手段】この発明に係る呼吸推定
方法は、ＱＲＳベクトルの角度だけでなくＴベクトルの
角度も考慮して、前記両ベクトルの角度データを一つの
データファイルに合成して、補間するデータ数を２倍に
して呼吸曲線を推定することにより、前述の問題点であ
る（呼吸数÷心拍数）≧１／２となるような場合でも精
度の高い推定が出来るようにしたものである。
方法は、ＱＲＳベクトルの角度だけでなくＴベクトルの
角度も考慮して、前記両ベクトルの角度データを一つの
データファイルに合成して、補間するデータ数を２倍に
して呼吸曲線を推定することにより、前述の問題点であ
る（呼吸数÷心拍数）≧１／２となるような場合でも精
度の高い推定が出来るようにしたものである。
【００１０】また、被測定者の四肢の動きの制約を取り
除くとともに計測する心電図を１チャンネルとするた
め、前記ＱＲＳベクトルの角度変化が振幅変化となる射
影軸を決め、この射影軸に沿って前記被測定者の胸部に
装着する電極位置を決定し、胸部誘導心電図１チャンネ
ルだけを計測し、この心電図のＱＲＳの大きさから呼吸
波形を推定するようにしたものである。
除くとともに計測する心電図を１チャンネルとするた
め、前記ＱＲＳベクトルの角度変化が振幅変化となる射
影軸を決め、この射影軸に沿って前記被測定者の胸部に
装着する電極位置を決定し、胸部誘導心電図１チャンネ
ルだけを計測し、この心電図のＱＲＳの大きさから呼吸
波形を推定するようにしたものである。
【００１１】さらに、上記呼吸波形推定の精度を向上さ
せるため、ＱＲＳベクトルの角度変化及びＴベクトルの
角度変化が各々ＱＲＳ波、Ｔ波の振幅変化となる１つの
射影軸を決め、この射影軸から胸部に装着する電極装着
位置を決定して導出された心電図のＱＲＳ波の大きさと
Ｔ波の大きさから呼吸波形を推定するようにしたもので
ある。
せるため、ＱＲＳベクトルの角度変化及びＴベクトルの
角度変化が各々ＱＲＳ波、Ｔ波の振幅変化となる１つの
射影軸を決め、この射影軸から胸部に装着する電極装着
位置を決定して導出された心電図のＱＲＳ波の大きさと
Ｔ波の大きさから呼吸波形を推定するようにしたもので
ある。
【００１２】
【作用】この発明における請求項１記載の呼吸推定方法
は、ＱＲＳベクトル及びＴベクトルの両角度変化を用い
て呼吸波形を推定することにより、推定精度が向上す
る。
は、ＱＲＳベクトル及びＴベクトルの両角度変化を用い
て呼吸波形を推定することにより、推定精度が向上す
る。
【００１３】また、請求項２及び請求項３記載の呼吸推
定方法は、四肢に電極を装着する必要がないので、両腕
を動かすことができ被測定者の負荷軽減となるとともに
計測する心電図は１チャンネルだけなので、測定系の簡
略化が行える。
定方法は、四肢に電極を装着する必要がないので、両腕
を動かすことができ被測定者の負荷軽減となるとともに
計測する心電図は１チャンネルだけなので、測定系の簡
略化が行える。
【００１４】
実施例１．図１は請求項１記載の一実施例を示すフロー
チャートである。図１において１１は従来例と同様に標
準肢誘導によって計測された２チャンネルの心電図デー
タである。１２〜１５においては従来例の図８の８２及
び８３のＱＲＳベクトル角度の計算に加えＴベクトルの
角度を計算するように構成されている。１６では前記１
心拍毎の角度データを、ＱＲＳ、Ｔ、ＱＲＳ、Ｔ・・・
というように同一時間軸上に各々のピーク時点の時間関
係はそのままで生起順にデータを並べ、１つのデータ列
にまとめる。このときＱＲＳベクトルとＴベクトルの角
度変化の大きさが異なっている場合は、予め求めておい
た補正係数を用いてＴベクトルの角度データを補正す
る。１７および１８では、前記データ列に補間を行い連
続曲線に変換し、推定呼吸波形として出力する。
チャートである。図１において１１は従来例と同様に標
準肢誘導によって計測された２チャンネルの心電図デー
タである。１２〜１５においては従来例の図８の８２及
び８３のＱＲＳベクトル角度の計算に加えＴベクトルの
角度を計算するように構成されている。１６では前記１
心拍毎の角度データを、ＱＲＳ、Ｔ、ＱＲＳ、Ｔ・・・
というように同一時間軸上に各々のピーク時点の時間関
係はそのままで生起順にデータを並べ、１つのデータ列
にまとめる。このときＱＲＳベクトルとＴベクトルの角
度変化の大きさが異なっている場合は、予め求めておい
た補正係数を用いてＴベクトルの角度データを補正す
る。１７および１８では、前記データ列に補間を行い連
続曲線に変換し、推定呼吸波形として出力する。
【００１５】前記のような呼吸推定方法においては１心
拍について得られる角度データが２つであり、補間に使
われるデータはＱＲＳベクトルのみの場合の２倍であ
る。従って、従来例に見られる心拍数と呼吸数の関係が
（呼吸数÷心拍数）≧１／２となるような場合でも精度
の高い推定ができることになる。
拍について得られる角度データが２つであり、補間に使
われるデータはＱＲＳベクトルのみの場合の２倍であ
る。従って、従来例に見られる心拍数と呼吸数の関係が
（呼吸数÷心拍数）≧１／２となるような場合でも精度
の高い推定ができることになる。
【００１６】実施例２． 図２はＱＲＳベクトルの呼吸による角度変化をＱＲＳの
大きさの単調変化として検出するための胸部に装着する
電極位置の決め方の説明図である。図において２１は呼
吸によって変化するＱＲＳベクトルを表わし、２２は前
記ＱＲＳベクトルの呼吸変化の角度差を表わす。このと
き心臓位置を通る２３の破線で示される方向に、前記Ｑ
ＲＳベクトルを射影すると２４で示されるように、前記
角度差に対応する大きさの差となる。従って、２５で示
される位置に電極を装着して胸部より心電図を検出すれ
ば、ＱＲＳの大きさの変化から呼吸波形の推定ができ
る。
大きさの単調変化として検出するための胸部に装着する
電極位置の決め方の説明図である。図において２１は呼
吸によって変化するＱＲＳベクトルを表わし、２２は前
記ＱＲＳベクトルの呼吸変化の角度差を表わす。このと
き心臓位置を通る２３の破線で示される方向に、前記Ｑ
ＲＳベクトルを射影すると２４で示されるように、前記
角度差に対応する大きさの差となる。従って、２５で示
される位置に電極を装着して胸部より心電図を検出すれ
ば、ＱＲＳの大きさの変化から呼吸波形の推定ができ
る。
【００１７】図３は上記の如く計測された心電図から呼
吸波形を推定するための一実施例のフローチャートであ
る。図において３１は、計測した１チャンネルの心電図
をサンプルしたデータであり、３２はこのデータから前
記実施例１同様にＱＲＳピーク時点を検出する。３３で
は３２で求めたピーク時点のＱＲＳ波の大きさを求め
る。ＱＲＳの大きさとしては図４に示すように、４１の
基線Ｂを基準としたＱ、Ｒ、Ｓ振幅の代数和（基線Ｂを
基準とした符号つき加算）即ち（−ＢＱ）＋（＋ＢＲ）
＋（−ＢＳ）、または４２〜４４で示される各々の面積
Ｓ1、Ｓ2、Ｓ3 で算出されるＳ 2 −Ｓ 1 −Ｓ 3 、または
この面積差から計算される平均電位等を用いる。このよ
うにして求めたＱＲＳ波の大きさのデータ列を１７にお
いて補間し連続波形に変換し、１８で呼吸波形として出
力する。
吸波形を推定するための一実施例のフローチャートであ
る。図において３１は、計測した１チャンネルの心電図
をサンプルしたデータであり、３２はこのデータから前
記実施例１同様にＱＲＳピーク時点を検出する。３３で
は３２で求めたピーク時点のＱＲＳ波の大きさを求め
る。ＱＲＳの大きさとしては図４に示すように、４１の
基線Ｂを基準としたＱ、Ｒ、Ｓ振幅の代数和（基線Ｂを
基準とした符号つき加算）即ち（−ＢＱ）＋（＋ＢＲ）
＋（−ＢＳ）、または４２〜４４で示される各々の面積
Ｓ1、Ｓ2、Ｓ3 で算出されるＳ 2 −Ｓ 1 −Ｓ 3 、または
この面積差から計算される平均電位等を用いる。このよ
うにして求めたＱＲＳ波の大きさのデータ列を１７にお
いて補間し連続波形に変換し、１８で呼吸波形として出
力する。
【００１８】前記の呼吸推定方法では、従来例の呼吸推
定方法と推定精度は同じであるが、胸部より１チャンネ
ルの心電図を導出するだけで済むので測定系の簡略化が
できるとともに、被測定者の負担軽減が行える。
定方法と推定精度は同じであるが、胸部より１チャンネ
ルの心電図を導出するだけで済むので測定系の簡略化が
できるとともに、被測定者の負担軽減が行える。
【００１９】実施例３．図５は上記のＱＲＳ波の大きさ
から推定した場合より推定精度を上げるための一実施例
のフローチャートである。図において、５１の心電図デ
ータは図３にＴベクトルの方向も加えて決定した電極位
置による胸部誘導心電図のサンプルデータである。３２
〜５３はＱＲＳ波とＴ波の大きさを求める。Ｔ波の大き
さもＱＲＳ波同様、ピーク値、面積、平均値等が用いら
れる。５４では前記実施例１と同様にＴ波の大きさの補
正を行い同一時間軸上に生起時間順にデータを並べて、
データ数２倍の補間データを作成する。１７では前記Ｑ
ＲＳ波の大きさの補間と同様にデータを補間し連続曲線
を求める。１８では求めた前記連続曲線を推定呼吸波形
として出力する。
から推定した場合より推定精度を上げるための一実施例
のフローチャートである。図において、５１の心電図デ
ータは図３にＴベクトルの方向も加えて決定した電極位
置による胸部誘導心電図のサンプルデータである。３２
〜５３はＱＲＳ波とＴ波の大きさを求める。Ｔ波の大き
さもＱＲＳ波同様、ピーク値、面積、平均値等が用いら
れる。５４では前記実施例１と同様にＴ波の大きさの補
正を行い同一時間軸上に生起時間順にデータを並べて、
データ数２倍の補間データを作成する。１７では前記Ｑ
ＲＳ波の大きさの補間と同様にデータを補間し連続曲線
を求める。１８では求めた前記連続曲線を推定呼吸波形
として出力する。
【００２０】
【発明の効果】この発明は、以上説明したような方法で
呼吸波形を推定するので、以下に記載される効果を奏す
る。
呼吸波形を推定するので、以下に記載される効果を奏す
る。
【００２１】１心拍について、ＱＲＳベクトルとＴベク
トルの２つの角度データを用いることにより従来法の２
倍のデータが得られ、推定呼吸波形の精度が向上すると
ともに心拍数に対して呼吸速度が速い場合にも推定可能
となる。
トルの２つの角度データを用いることにより従来法の２
倍のデータが得られ、推定呼吸波形の精度が向上すると
ともに心拍数に対して呼吸速度が速い場合にも推定可能
となる。
【００２２】また、心ベクトルに対し異なる角度をなす
２種類の誘導による心電図を計測して得られるＱＲＳベ
クトルの方向を基にして胸部に装着する電極位置を決定
し、この電極位置で計測される心電図のＱＲＳ波の大き
さの変化から呼吸波形を推定することにより、計測すべ
き心電図は１チャンネルで済む。さらに被測定者の四肢
に電極を付けなくても良いため、被測定者の負担軽減が
できる。
２種類の誘導による心電図を計測して得られるＱＲＳベ
クトルの方向を基にして胸部に装着する電極位置を決定
し、この電極位置で計測される心電図のＱＲＳ波の大き
さの変化から呼吸波形を推定することにより、計測すべ
き心電図は１チャンネルで済む。さらに被測定者の四肢
に電極を付けなくても良いため、被測定者の負担軽減が
できる。
【００２３】また、胸部誘導心電図１チャンネルのデー
タについても、Ｔ波の大きさを推定データに使用するこ
とにより推定精度の向上が図れる。
タについても、Ｔ波の大きさを推定データに使用するこ
とにより推定精度の向上が図れる。
【００２４】また、ホルター心電計のように長時間にわ
たって心電図を記録し、検査や診断を行う場合がよくあ
るが、本発明を応用すれば、従来より行われている心電
図の解析に加え、推定呼吸波形の検査及び診断も行える
という利点がある。
たって心電図を記録し、検査や診断を行う場合がよくあ
るが、本発明を応用すれば、従来より行われている心電
図の解析に加え、推定呼吸波形の検査及び診断も行える
という利点がある。
【図１】この発明の実施例１を示すフローチャートであ
る。
る。
【図２】電極装着位置を決める説明図である。
【図３】この発明の実施例２を示すフローチャートであ
る。
る。
【図４】ＱＲＳ波の大きさの説明図である。
【図５】この発明の実施例３を示すフローチャートであ
る。
る。
【図６】ＱＲＳベクトルの求め方の説明図である。
【図７】従来の呼吸推定法のハードウェアの構成図であ
る。
る。
【図８】従来の呼吸推定法を示すフローチャートであ
る。
る。
１１ 標準肢誘導心電図２チャンネルのサンプルデータ １２ Ｑ波、Ｒ波、Ｓ波の各ピーク時点の検出 １３ Ｔ波のピーク時点の検出 １４ ＱＲＳベクトルの角度の計算 １５ Ｔベクトルの角度の計算 １６ ＱＲＳベクトル、Ｔベクトルの両角度データの１
つのデータ列に変換 １７ データ補間 １８ 推定呼吸波形の出力 ２１ ＱＲＳベクトル ２３ 呼吸によるＱＲＳベクトルの変化角度 ２３ 射影軸 ２４ ＱＲＳの大きさ ２５ 電極 ３１ ＱＲＳベクトルを基に決めた電極位置で計測した
心電図データ ３２ ＱＲＳ波のピーク検出 ３３ ＱＲＳ波の大きさ算出 ４１ 心電図基線 ４２ Ｑ波の面積 ４３ Ｒ波の面積 ４４ Ｓ波の面積 ５１ ＱＲＳとＴベクトルから決めた電極位置で計測し
た心電図データ ５２ Ｔ波ピークの検出 ５３ Ｔ波の大きさの算出 ５４ 補間データの作成 ６１ 第Ｉ誘導心電図のＱＲＳ波の大きさ ６２ 第ＩＩＩ誘導心電図のＱＲＳ波の大きさ ７１ 心電図増幅器 ７２ Ａ／Ｄ変換器 ７３ コンピュータ ８１ ２チャンネル心電図データの移動平均 ８２ ＱＲＳ波決定部 ８３ ＱＲＳベクトル算出部 ８４ ３次スプライン補間 ８５ 推定呼吸波形出力
つのデータ列に変換 １７ データ補間 １８ 推定呼吸波形の出力 ２１ ＱＲＳベクトル ２３ 呼吸によるＱＲＳベクトルの変化角度 ２３ 射影軸 ２４ ＱＲＳの大きさ ２５ 電極 ３１ ＱＲＳベクトルを基に決めた電極位置で計測した
心電図データ ３２ ＱＲＳ波のピーク検出 ３３ ＱＲＳ波の大きさ算出 ４１ 心電図基線 ４２ Ｑ波の面積 ４３ Ｒ波の面積 ４４ Ｓ波の面積 ５１ ＱＲＳとＴベクトルから決めた電極位置で計測し
た心電図データ ５２ Ｔ波ピークの検出 ５３ Ｔ波の大きさの算出 ５４ 補間データの作成 ６１ 第Ｉ誘導心電図のＱＲＳ波の大きさ ６２ 第ＩＩＩ誘導心電図のＱＲＳ波の大きさ ７１ 心電図増幅器 ７２ Ａ／Ｄ変換器 ７３ コンピュータ ８１ ２チャンネル心電図データの移動平均 ８２ ＱＲＳ波決定部 ８３ ＱＲＳベクトル算出部 ８４ ３次スプライン補間 ８５ 推定呼吸波形出力
───────────────────────────────────────────────────── フロントページの続き (72)発明者 寺下 裕美 尼崎市塚口本町８丁目１番１号 三菱電 機株式会社 中央研究所内 (56)参考文献 特開 昭52−135587（ＪＰ，Ａ) 特開 昭50−97181（ＪＰ，Ａ) 特開 平１−99530（ＪＰ，Ａ) 実開 昭50−91392（ＪＰ，Ｕ) 特公 昭49−39796（ＪＰ，Ｂ１) (58)調査した分野(Int.Cl.6，ＤＢ名) A61B 5/0472 A61B 5/08
Claims (3)
- 【請求項１】 心ベクトルに対し異なる角度をなす２種
類の誘導による心電図を計測し、前記２種類の心電図の
ＱＲＳ波よりＱＲＳベクトルを、Ｔ波よりＴベクトルを
各々求め、前記両ベクトルの各々の角度変化を用いて呼
吸波形を推定することを特徴とする心電図を用いた呼吸
波形推定方法。 - 【請求項２】 請求項１記載の２種類の心電図によるＱ
ＲＳベクトルの方向を求め、前記ベクトルの方向を基に
電極の胸部装着位置を決定し、前記電極位置で計測した
１チャンネルの心電図のＱＲＳ波の大きさより呼吸波形
を推定することを特徴とする心電図を用いた呼吸波形推
定方法。 - 【請求項３】 胸部誘導電極の装着位置をＴベクトルの
方向も考慮して決定し、呼吸波形の推定もＴ波の大きさ
を考慮して行うことを特徴とする請求項２記載の心電図
を用いた呼吸波形推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3107073A JP2820169B2 (ja) | 1991-05-13 | 1991-05-13 | 心電図を用いた呼吸波形推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3107073A JP2820169B2 (ja) | 1991-05-13 | 1991-05-13 | 心電図を用いた呼吸波形推定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04336035A JPH04336035A (ja) | 1992-11-24 |
| JP2820169B2 true JP2820169B2 (ja) | 1998-11-05 |
Family
ID=14449801
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3107073A Expired - Fee Related JP2820169B2 (ja) | 1991-05-13 | 1991-05-13 | 心電図を用いた呼吸波形推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2820169B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9980664B2 (en) | 2013-03-22 | 2018-05-29 | Panasonic Intellectual Property Management Co., Ltd. | Biological signal measurement system, apparatus, method, and computer program thereof |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6413397B2 (ja) * | 2014-06-30 | 2018-10-31 | Tdk株式会社 | 呼吸状態推定装置、呼吸状態推定方法及びプログラム |
| JP6788656B2 (ja) * | 2015-07-21 | 2020-11-25 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | Ｅｃｇデータにおける活動の自動ラベリングのための方法及びシステム |
| CN111603166A (zh) * | 2020-06-01 | 2020-09-01 | 河南华南医电科技有限公司 | 一种基于心电向量图的呼吸信号检测方法 |
-
1991
- 1991-05-13 JP JP3107073A patent/JP2820169B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9980664B2 (en) | 2013-03-22 | 2018-05-29 | Panasonic Intellectual Property Management Co., Ltd. | Biological signal measurement system, apparatus, method, and computer program thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH04336035A (ja) | 1992-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7343198B2 (en) | System, software, and method for detection of sleep-disordered breathing using an electrocardiogram | |
| US10617325B2 (en) | Respiratory rate measurement using a combination of respiration signals | |
| US6415174B1 (en) | ECG derived respiratory rhythms for improved diagnosis of sleep apnea | |
| US5025784A (en) | Apparatus and method for detecting and processing impedance rheogram | |
| EP0606301A1 (en) | Methods of impedance cardiography and heartbeat determination as well as an apparatus for implementing said methods | |
| Tayler et al. | Signal distortion in the electrocardiogram due to inadequate phase response | |
| JPH1170089A (ja) | 心臓波形特性表示方法 | |
| US5913308A (en) | Apparatus and method for determining respiratory effort from muscle tremor information in ECG signals | |
| US20030069512A1 (en) | Method and system for measuring T-wave alternans by alignment of alternating median beats to a cubic spline | |
| Zhou et al. | Philips QT interval measurement algorithms for diagnostic, ambulatory, and patient monitoring ECG applications | |
| Mason et al. | Quantitative assessment of respiratory derivation algorithms | |
| US7962201B2 (en) | Methods of generating a virtual lead associated with a physiological recording | |
| US4269195A (en) | Apparatus for measuring a pulmonary function | |
| CENGIZ et al. | The effect of ventilation on the accuracy of pulmonary artery and wedge pressure measurements | |
| Tarassenko et al. | Multi-sensor fusion for robust computation of breathing rate | |
| JP2820169B2 (ja) | 心電図を用いた呼吸波形推定方法 | |
| Bowers et al. | Respiratory rate derived from principal component analysis of single lead electrocardiogram | |
| JP2509828B2 (ja) | 呼吸波測定装置 | |
| JP4323277B2 (ja) | 睡眠時無呼吸型判別方法 | |
| Weltman et al. | Impedance pneumograph recording across the arms | |
| Govreen-Segal et al. | Real-time PC-based system for dynamic beat-to-beat QT-RR analysis | |
| Negi et al. | Development of a real-time breathing-rate monitor using difference operation method and adaptive windowing on dry-electrode ECG signal | |
| EP4085826B1 (en) | Smartphone heart rate and breathing rate determination using accuracy measurement weighting | |
| JP2829389B2 (ja) | 呼吸検出装置 | |
| EP2037800B1 (en) | Method and analytical device for the analysis of respiration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |