JP2003525675A - マルチチャンネルecg波形の電話伝送監視 - Google Patents

マルチチャンネルecg波形の電話伝送監視

Info

Publication number
JP2003525675A
JP2003525675A JP2001564657A JP2001564657A JP2003525675A JP 2003525675 A JP2003525675 A JP 2003525675A JP 2001564657 A JP2001564657 A JP 2001564657A JP 2001564657 A JP2001564657 A JP 2001564657A JP 2003525675 A JP2003525675 A JP 2003525675A
Authority
JP
Japan
Prior art keywords
signal
frequency
composite
local
ecg
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.)
Pending
Application number
JP2001564657A
Other languages
English (en)
Inventor
アンソニー マーケシニ,
ナルジス ネイデノフ,
Original Assignee
ペースアート アソシエイツ, リミテッド パートナーシップ
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 ペースアート アソシエイツ, リミテッド パートナーシップ filed Critical ペースアート アソシエイツ, リミテッド パートナーシップ
Publication of JP2003525675A publication Critical patent/JP2003525675A/ja
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0004Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by the type of physiological signal transmitted
    • A61B5/0006ECG or EEG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/30Input circuits therefor
    • A61B5/307Input circuits therefor specially adapted for particular uses
    • A61B5/308Input circuits therefor specially adapted for particular uses for electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/30Input circuits therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/904Telephone telemetry

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Cardiology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Burglar Alarm Systems (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Selective Calling Equipment (AREA)

Abstract

(57)【要約】 ECG波形等の生物学的信号の遠隔監視であって、ECG波形のFMエンコードバージョンを患者の家庭から中央ステーションへ送信し、そこで解析を行ってFM信号のゼロ交差よりもより多くの時間点において局所的周波数を推定する。複数個のECG波形を複合FM信号にエンコードし、そのうちのいずれか1つの点は同時にECG波形の数を表わすことが可能である。該複合波形を中央ステーションにおいて解析して周波数ドメインにおける分離により個々の波形を抽出し且つ再生する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】
本発明は、例えばECG(心電図)データ等の生物医学的データの遠隔的監視
の分野に関するものである。
【0002】
【従来の技術】
ECG波形及びペースメーカー情報等の心拍記録データの電話伝送監視は長年
の間使用されている。典型的に、患者の家庭における心臓トランスデューサは患
者の身体と電気的に接触している一対のECGパッドの間における電圧の形態で
の電気的ECG信号を発生する。その結果発生するECG波形は搬送波を周波数
変調するために使用され、その結果発生するFM信号はスピーカーを駆動して電
話受話器のマウスピース内に音響FM信号を発生し、該受話器は該音響信号を電
気的FM信号へ変換する。電話回路網を介して、中央ステーションは送信されて
くる信号を受取り且つそれを処理してECG波形を再生し、表示し、且つ記録し
、又はその他の情報を抽出する。心臓信号がペースメーカーに関連している場合
には、興味のある情報はペースメーカーのパルスの期間即ちパルス間の時間であ
る場合がある。音響信号へ及び電気信号へ戻す変換は、患者が適切な機器及びそ
の目的のための技能を有している場合には回避することが可能である。心臓情報
の電話伝送監視の例は米国特許第4,938,229号及び第5,467,77
3号(各々を引用によって本明細書に取込む)、及び第5,735,285号及
びこれらの特許において引用されている文献において見出すことが可能である。
【0003】 本特許明細書の譲受人はこのような装置及びサービスを提供しており、それは
更にウエブサイト<<www.Paceart.com>>において記載してい
る。典型的に、異なるECGパッド又はパッドの組合わせ(ベクトル)から、又
はペースメーカーからの情報はFM信号内に直列的に埋め込まれ、且つ例えばI
D及びタイムスタンプ等の付加的な情報のパケットが通常挿入される。唯1つの
ECG波形、又は唯1つのペースマーカーパルス、又は唯1つの情報を識別する
項目が任意の1つの時間において搬送波を変調する。数年前においては、ウィス
コンシン州ミルウォーキーにおけるカーディアックエバリューションセンター(
Cardiac Evaluation Center)と呼ばれる事業体が2
チャンネル送信器及び私有の受信器を提供しており且つ未だに提供しているもの
と思われ、それは同時に2つのECG波形を単一のFM信号へエンコードし、該
信号は私有の受信器において分離されるものであり、その分離のためにアナログ
バンドパスフィルタを使用しているものと思われる。そのシステムがFM信号を
復調するためにどのような技術を使用しているかは不知である。
【0004】 本明細書に引用によって取込んだ2つの特許において、FM信号は受信ステー
ションにおいて復調されてFM信号のゼロ交差を見つけ出し且つこれらのゼロ交
差の間の時間を測定することによってECG波形を抽出している。特に、特許さ
れているシステムは隣接するゼロ交差の間のインターバル期間中にクロックをカ
ウントし且つそのカウントを周波数へ変換し、それによりオリジナルのECG波
形を再生している。特許第5,735,285号は別のゼロ交差検知技術を提案
しているものと理解され、FM信号のデジタル化したサンプルが正の値と負の値
との間で遷移する区域の検査を行う。このようなFM信号のゼロ交差解析は、ノ
イズ及びその他の不正確性の発生源を適宜抑圧した状態で数年にわたって使用さ
れているが、(1)オリジナルのECG波形のより正確且つ信頼性のある再生、
(2)複数個のECG波形又はこのようなオリジナルの情報のより正確且つ信頼
性のある再生と結合されたその他の情報の同時的伝送、(3)異なるフォーマッ
トの情報伝送へプログラミングによって便宜的に適合させることが可能であり且
つ私有のハードウエアよりもより廉価であり且つより許容可能である汎用コンピ
ュータを使用してこのような同時的伝送を受信ステーションにおいて復調させる
ことに対する必要性が未だに残存している。
【0005】
【課題を解決するための手段】
本特許明細書は遠隔ステーションに存在する心臓状態を中央ステーションにお
いて遠隔的にモニタ即ち監視するシステム及び方法を開示している。好適実施例
においては、3個又はそれ以上のECG波形がローカル(局所的)ステーション
における患者から派生される。これらのECG波形は夫々の異なる搬送波を周波
数変調してそれにより3個又はそれ以上の夫々のFM信号を発生する。これらの
FM信号は3個又はそれ以上のECG波形からの同時的な情報を包含する複合F
M信号に結合され、且つ中央ステーションへ送信される。中央ステーションにお
いては、受取った複合FM信号を時間ドメイン及び周波数ドメインの両方におい
て処理して、複合FM信号のゼロ交差よりも実質的に互いにより近い複合FM信
号の部分において局所的周波数を推定することを包含するような態様で3個又は
それ以上の個別的なECG波形を再生する。本プロセスは局所的位相差を計算し
且つそれらを使用して局所的周波数を推定する。局所的位相差は、複合FM信号
のデジタルサンプルを該デジタルサンプルの位相シフトさせたバージョン(もの
)と結合させることによって計算される。
【0006】
【発明の実施の形態】
心臓ECGイベントモニタ、ループレコーダー及びポストイベントレコーダー
は患者の家庭等の遠隔位置において使用されるセンサー/送信器の例である。マ
ルチチャンネルセンサー/送信器は、例えば、ECGパッドの異なる組合わせか
らECG波形の幾つかのチャンネルを発生する。1つの公知のパターンは、3対
のECGパッドからの信号の間の差異から派生される3個のECG波形を使用す
るものである(3個のベクトル)。好適実施例においては、本明細書に記載する
システムは遠隔位置において3個のECG波形を単一のFM信号へ同時的にエン
コードし、且つゼロ交差検知よりもより信頼性のある技術を使用して受信ステー
ションにおいてECG波形を再生する。好適実施例においては、ECG波形の再
生は適宜のユーティリティ及びアプリケーションプログラムを稼動するPC等の
汎用コンピュータのみによって実現される。以下の説明は同時的に単一のFM信
号へエンコードされる3個のECG波形の例を使用するが、その一般的な形態に
おいては、本明細書の開示はN個の波形(尚、N≧2)に対して、及びECG波
形に加えて生物医学的信号に対して適用可能である。更に、ゼロ交差検知を使用
する場合よりもより信頼性を持って波形を再生するための本明細書において開示
する技術は、単一の波形がFM信号内にエンコードされるに過ぎない場合に対し
ても適用可能である。
【0007】 送信端(典型的に、患者の家庭)において、患者は本特許の譲受人によって現
在供給されているものとFMエンコーディング技術において同様のセンサー/送
信器を使用するものであるが、FMは3個のECG波形の各々を夫々のチャンネ
ル内へエンコードし、次いで、それら3個のFMエンコードされたチャンネルを
単一の複合FM信号へ加算する。例えば、第一ECGチャンネルFMはECGデ
ータの第一チャンネルに対して1500−1900Hzの周波数帯域における1
700Hz搬送波を変調し、第二ECGチャンネルは2150Hz搬送波上で1
950−2350Hz帯域を使用し、且つ第三ECGチャンネルは2600Hz
搬送波上で2400−2800Hz帯域を使用する。その結果得られる3個のF
M信号は複合FM信号へ加算され、それは中央ステーションへ送信される。例え
ば、これらに制限されるものではないが、送信装置のID、ペースメーカーパル
ス測定値、及びタイムスタンプ等の付加的なデータを、例えば周波数を信号帯域
周波数の外側に数ミリ秒シフトさせ、従ってFMアナログ信号内に二進データの
存在を表わす公知の且つ長い間使用されている技術によって複合FM信号内に埋
め込むことが可能である。その送信は、最初に、複合FM信号をセンサー/送信
器におけるスピーカーによって音響信号へ変換させ且つ該音響信号を中央又は受
信ステーションへ公衆電話システムを介して接続されている電話の受信器におい
て動作させることにより行うことが可能であり、又は直接的な電気的送信を使用
することが可能であり、それはオーディオステージを介して行われるものではな
い。
【0008】 中央ステーションにおいて、電話線を介して受取られた複合FM信号はサウン
ドカードを具備するPC等の汎用デジタルコンピュータへ供給され、且つ解析さ
れて3個のECGチャンネルを定義する情報、及びその中にエンコードされてい
る場合のある付加的な情報を抽出する。原理的には、中央ステーションにおいて
実施されるプロセスが受取られた複合FM信号をデジタルサンプルx(t)へ変
換し、それらを個別的なFM信号とマッチングする周波数帯域へ分離し、各帯域
に対して隣接するデジタルサンプルの間の瞬間的な周波数における差異を見つけ
出し、且つこれらの周波数差を使用してオリジナルのECG波形及びその他の興
味のあるデータを再生させる。
【0009】 本明細書において開示する実施例を使用したシステムを例示した図1を参照す
ると、例えば患者の家庭等の遠隔位置において公知の如く使用されるECGパッ
ド10はECGアナログ波形の3個のチャンネル即ちベクトルを発生する。局所
的送信器11はFMエンコーダー12−1,12−2,12−3を有しており、
その各々はECGデータの夫々のチャンネルを夫々の周波数帯域において周波数
変調したアナログ波形へエンコードする。局所的送信器11は、又、加算装置1
4を有しており、それはこれら3個のFMチャンネルを単一の複合FM信号へ結
合させる。ペースメーカー16が使用される場合には、加算装置14は当該技術
において公知の如くペースメーカーに関連した情報を複合FM信号内に埋め込む
ことが可能である。更に当該技術において公知の如く、加算装置14は例えば局
所的送信器のID、タイムスタンプ等のその他の情報を複合FM信号内に埋め込
むことも可能である。遠隔位置におけるスピーカー18は複合FM信号をオーデ
ィオ信号へ変換し、それは、遠隔位置の電話20の受信器によって複合アナログ
電気FM信号へ変換される。このFM信号は公衆電話交換システムを介して、又
は何等かのその他の通信リンクを介して解析及び記録を行うために中央ステーシ
ョンへ送信される。特別の装置及び技能が遠隔位置において使用可能である場合
には、オーディオ信号へ及びアナログ電気信号へ戻す変換を回避することが可能
であり、且つ加算装置14からの複合FM信号は適宜の通信リンクを使用して中
央ステーションへ直接的に送信することが可能である。
【0010】 中央位置において、電話ユニット30が複合FM信号を受取り且つそれを、例
えばサウンドカードを具備するPC等の適宜プログラムされた汎用コンピュータ
へ供給する。サウンドカードをADC(アナログ・デジタル変換器)32として
使用して、該中央ステーションは受取った複合FM信号を時間ドメインデジタル
サンプルx(t)のアレイへ変換し、それは、FFT(高速フーリエ変換)解析
器34が周波数ドメインデジタルサンプルfft(t)のアレイへ変換する。こ
れらのサンプルfft(t)はバンドパスフィルタ36−1,36−2,36−
3においてECG情報の3個のチャンネルに対応する3個の帯域へ分離され、且
つ各バンドパスフィルタの出力は解析器38−1,38−2,38−3の夫々の
1つにおいてIFFT(逆高速フーリエ変換)解析が行われる。ユニット40は
これらの解析器の出力を受取り且つ3個のECG波形を再生し、記録し、且つ表
示する。例えばペースメーカー情報等の付加的な情報が複合FM信号内に埋め込
まれている場合には、ペースメーカー解析器42がそれを抽出し且つそれを表示
及び記録のためにユニット40へ供給する。中央ステーションにおける装置は、
PCをプログラミングすることによって実現することが可能であり、且つ好適実
施例においてはそうである。前述したように、従来のPCサウンドカードハード
ウエア及びPCのユーティリティを使用して複合FM信号をデジタル化する。F
FT解析器34は在庫品のFFTプログラムを使用して実現することが可能であ
る。バンドパスフィルタ36は以下に説明するようにある周波数ビンをヌル化す
ることにより実現することが可能である。IFFT解析器38は在庫品のIFF
T及びヒルベルト変換プログラムを使用することによって実現することが可能で
ある。ペースメーカー解析器42は当該技術において公知であり且つ単一チャン
ネルECGデータ用に何年もの間本譲受人に使用されているように実現すること
が可能である。最後に、ユニット40はPCの従来のデータの格納及び表示能力
を使用して実現することが可能である。
【0011】 別の実施例においては、遠隔位置電話セット20から開始し且つ中央ステーシ
ョン電話セット30で終了する経路を、スピーカー18からの音をアナログ電気
信号へ変換させるためにスピーカー18の充分近くにあるマイクロフォン31に
よって置換させることが可能であり、該アナログ信号は、次いで、ADC32へ
供給される。更なる実施例(図面には示していない)としては、加算回路14か
らのアナログ電気信号を直接的にADC32へ供給することが可能であり、それ
によりスピーカー18から開始し中央ステーション電話セット30で終了する経
路を除去している。
【0012】 本明細書に開示した例示的で且つ非制限的な例においては、ECGベクトルへ
適用される処理は以下の主要なステップを包含しており、それらはサウンドカー
ド及び適宜のプログラミングを具備するPCを使用してコンピュータによって実
現する。
【0013】 1.患者の家庭、又は別の遠隔又は送信位置において、3個のECGチャンネ
ル(ベクトル)を得、その各々は夫々のECG電気波形の形態にある。このこと
は、例えば、本譲受人から入手可能な装置等の現在市販されている装置を使用し
て行うことが可能である。
【0014】 2.各ECGベクトルを使用して夫々の搬送波を周波数変調し、それにより3
個のFM ECG信号を発生し、その各々は夫々の周波数帯域内にあり、例えば
チャンネル位置に対しては1700Hzの搬送波周波数及び1500−1900
Hzの帯域幅であり、チャンネル2に対しては2150Hzの搬送波周波数及び
1950−2350Hzの帯域幅であり、且つチャンネル3に対しては2600
Hzの搬送波周波数及び2400−2800Hzの帯域幅である。各個別的なチ
ャンネルに対するエンコーディングは、又、例えば、本譲受人から現在市販され
ている装置を使用して行うことが可能である。
【0015】 3.これら3個のFM ECG信号を単一の複合FM信号へ加算する。このこ
とは、例えば、このタイプの現在市販されているアナログ加算回路を使用して行
うことが可能である。
【0016】 4.複合FM信号をオーディオ信号へ変換する。ことのことは、例えば本譲受
人から入手可能であるもの等の現在市販されている家庭送信器等のスピーカーを
使用して行うことが可能である。
【0017】 5.オーディオ信号を複合FM信号へ戻すべく変換し且つそれを中央ステーシ
ョンへ送信する。このことは、例えば本明細書に引用によって組込んだ特許に記
載されているような患者の家庭における電話セットを使用して行うことが可能で
ある。
【0018】 6.中央ステーションにおいて受けられた複合FM信号を時間ドメインデジタ
ルサンプルx(t)のアレイへデジタル化する。このことは、従来のPCのサウ
ンドカードを使用して行うことが可能である。好適なフォーマットは、8kHz
(毎秒8,000個のサンプル)のサンプリングレートで入力してくる複合FM
信号を各々が512個のサンプルによってオーバーラップされる1024サンプ
ルのアレイへデジタル化させることであり、即ち第二アレイの最初の512個の
サンプルは第一アレイの最後の512個のサンプルと同一であり、以下同様であ
る。各サンプルは16ビットの長さであり、夫々の時間(t)における複合FM
信号の瞬間的な振幅(x)を表わしている。その結果は1024のアレイのつな
がりであり、各々が16ビットの値であり、512個のサンプルによってオーバ
ーラップされている。好適実施例における計算の便宜上、該アレイは単一精度ア
レイへ変換される。
【0019】 7.デジタルサンプルx(t)を、例えば、ハミングウインドウ(Hammi
ng Window)フィルタでフィルタしてノイズを減少させる。このことは
、PCにおけるハミングウインドウフィルタ用の在庫品ユーティリティを使用し
て行うことが可能である。
【0020】 8.デジタルサンプルx(t)のアレイをFFT(高速フーリエ変換)解析器
を介して通過させて該サンプルを周波数ドメインデジタルサンプルアレイfft
(t)ヘ変換し、尚各サンプルはx(t)アレイのフーリエ級数表現の係数の値
である。このことは、PCにおいて稼動する在庫品であるFFTプログラムを使
用して行うことが可能である。その結果は1024要素x(t)アレイの各々の
対応する1025要素複素共役対称fft(t)アレイへの変換となる。fft
(t)アレイの要素は夫々の周波数に対する係数の値に関連しており、PCメモ
リにおける夫々の周波数ビン内に格納される。この時点において、3つのECG
信号の帯域幅外側の周波数に対する係数値をヌル即ちゼロとさせるために付加的
なフィルタ処理を行うことが可能である。例えば、750−3250Hzのバン
ドパスフィルタを、0Hz、7.8125Hz、15.625Hz、...、7
42.1875Hz(即ち、各fft(t)アレイに包含される要素1−95)
に対応するビン及び3257.8125Hz、3265.625Hz、...、
4000Hz(即ち、各fft(t)アレイの要素417−513)に対応する
周波数ビンをゼロ化することによって適用することが可能である。
【0021】 9.サンプルfft(t)を、各々が遠隔位置(患者の家庭)においてFMエ
ンコードされたECGチャンネルの夫々の1つの周波数帯域と一致する夫々の別
個の帯域に分割する。このことは、(フィルタした)fft(t)アレイの3つ
のコピーを作成し且つ各々において夫々ECG信号の周波数帯域外側の周波数に
対応する要素をヌル化即ちゼロとさせることによって行うことが可能である。
【0022】 10.サンプルfft(t)をIFFT(逆高速フーリエ変換)及びヒルベル
ト変換解析を介して通過させて解析信号z(t)のデジタルサンプルからなるア
レイを得、尚各z(t)サンプルは複合FM信号の時間ドメインサンプルx(t
)に一致する実部分及びx(t)のヒルベルト変換と一致する虚部分jh(t)
とを有しており、次式で表わされる。
【0023】
【数1】
【0024】 尚、ifftは逆フーリエ変換を示しており、 i=[0,N/(2−1)]の場合にはB(i)=2、 i=[N/2,N−1]の場合にはB(i)=0、 iはN個の要素からなるfft(t)アレイの要素を示しており、
【数1a】
【0025】 は共役を示しており、 fftは高速フーリエ変換を示しており、 jは虚部分を示しており、 h(t)はアレイx(t)のヒルベルト変換を示している。
【0026】 このことは、PC上でランされる在庫品のヒルベルト変換及びIFFTプログ
ラムを使用することによって行うことが可能である。上の式から明らかなように
、ヒルベルト変換はfft(t)アレイ(即ち、アレイ要素514−1025)
の全ての負の周波数ビンにおける係数値をゼロ化させ且つfft(t)アレイ(
即ち、要素1−513)の全ての正の周波数ビンにおける係数値を2倍とさせる
ことが関与する。その結果はIFFT処理され、各fft(t)アレイ(ヒルベ
ルト変換されたもの)を1024要素複素アレイz(t)ヘ変換し、その場合に
、(1)実部分は上述したウインドウ処理及びフィルタ処理によって向上された
オリジナルのデータx(t)を包含しており、且つ(2)虚部分は同一のオリジ
ナルデータのヒルベルト変換を包含している。
【0027】 11.次式に従ってx(t)の各サンプル位置に対する瞬間的な位相角p(t
)を見つけ出す。
【0028】
【数2】
【0029】 このことは、時間サンプル(t)の各々に対してすぐ上に記載した割算及び三角
関数の計算を実施するためにPCをプログラミングすることによって行うことが
可能である。その結果は、複合FM信号がアレイx(t)を発生するためにサン
プルされた各時刻(t)に対する位相角値p(t)である。
【0030】 12.次式に従ってx(t)の各サンプル位置に対する瞬間的周波数f(t)
を見つけ出す。
【0031】
【数3】
【0032】 尚、△p(t)は瞬間的位相角p(t)に関する2つの隣接するサンプルの間に
おける値の差であり、△tはx(t)の2つの隣接するサンプルの間の時間間隔
である。
【0033】 このことは、p(t)及び(t)の各対の隣接する値に対してすぐ上に記載し
た演算を実施するためにPCをプログラミングすることによって行うことが可能
であり、実効的には、各アレイx(t)に対して瞬間的な周波数値f(t)の1
024要素アレイを発生する。
【0034】 13.瞬間的な周波数f(t)を再生したECG波形のサンプルからなる振幅
へ変換し(56個のサンプル移動平均を使用)14ビット長の平均化し再生した
ECGサンプルを得る。このことは、各アレイf(t)の最初及び最後の25%
を最初に破棄することによって行うことが可能である(何故ならば、先のハフマ
ンウインドウフィルタ処理に起因してアレイ内のこれらの部分においては減衰が
大きいからである)。先に説明したx(t)アレイの50%のオーバーラップの
ために、1つのアレイf(t)から廃棄された要素は先行する及び後続のアレイ
内に存在しており、従ってこの処理は複合FM信号がサンプルされた各時刻に対
する瞬間的周波数f(t)を派生する。f(t)の56個のサンプルの移動平均
を使用する目的は、ノイズ又はその他のアーチファクト即ち人工的効果の影響を
減少させるためである。その結果は、142.85714Hz(即ち、オリジナ
ルのサンプリングレートの8,000Hzを56で割算したものであり、平均化
において使用したサンプル数)のレートにおける周波数の平均化した値のストリ
ングである。計算の便宜上、その結果得られる値を5で乗算し且つ整数の形態へ
変換し、夫々のECG信号において夫々の1/142.85714タイムスロッ
トにおける瞬間的周波数を表わす14ビット値のストリングを発生することが可
能である。
【0035】 14.FSK/ペースメーカーパルス解析に対するエッジ検知及びエンコード
の結果2ビットエンコーダーデータとなる。このことは、例えば、本譲受人から
入手可能な市販の装置において現在実施されている如くに行うことが可能である
。原理的には、その処理は複合FM信号において高周波数高振幅エッジを検知し
、FSK(周波数シフトキー)/ペースメーカーパルス解析を実施し、且つ検知
したFSK/ペースメーカーパルスデータを相次ぐ2ビット値としてエンコーデ
ィングすることを包含している。
【0036】 15.結果的に得られるデータを142.85714Hzにおいて16ビット
サンプルへフォーマット化し、その場合に上の2ビットはFSK/ペースメーカ
ーデータであり、そのうちの20ビットは10個の16ビットサンプルにわたっ
て格納される。
【0037】 16.再生されたECG及びその他の関連するデータを表示/記録する。この
ことは、例えば、本譲受人によって現在商業的に使用されており、且つこの場合
には3チャンネルデータを表示するべく適合されている単チャンネルECGデー
タに対して引用により本明細書に取込んだ特許に記載されているような周波数対
振幅変換技術を使用して行うことが可能である。
【図面の簡単な説明】
【図1】 本開示の好適な例を具体化したシステムを示した概略図。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 ネイデノフ, ナルジス アメリカ合衆国, ニュージャージー 07470, ウエイン, ハートウィック レーン 87 Fターム(参考) 4C027 AA02 DD04 JJ01

Claims (20)

    【特許請求の範囲】
  1. 【請求項1】 遠隔ステーションに存在する心臓状態を中央ステーションに
    おいて遠隔的に監視する方法において、 ローカルステーションにおいて3個又はそれ以上のECG波形を派生させ且つ
    前記波形の各々に従って夫々の搬送波を周波数変調させて3個又はそれ以上の夫
    々のFM信号を生成し、 前記FM信号を前記3個又はそれ以上のECG波形からの同時的な情報を包含
    する複合FM信号へ結合させ、 前記複合FM信号を中央ステーションへ送信し、 前記中央ステーションにおいて受取られた前記複合FM信号を時間領域及び周
    波数領域の両方において処理して3個又はそれ以上の個別的なECG波形を再生
    する、 ことを包含している方法。
  2. 【請求項2】 請求項1において、前記処理が前記複合FM信号のゼロ交差
    以外の実質的に互いにより近い前記複合FM信号の部分における局所的周波数を
    推定することを包含している方法。
  3. 【請求項3】 請求項2において、前記処理が、前記ゼロ交差よりも実質的
    により近い前記FM信号の部分における局所的位相差を推定し且つ前記局所的周
    波数を推定するために前記局所的位相差を使用することを包含している方法。
  4. 【請求項4】 請求項3において、前記局所的位相差を推定する場合に、前
    記複合FM信号のデジタルサンプルを前記デジタルサンプルの位相シフトさせた
    ものと結合させることを包含している方法。
  5. 【請求項5】 請求項4において、前記比較を行う場合に、前記サンプルを
    前記サンプルのヒルベルト変換したものと比較することを包含している方法。
  6. 【請求項6】 請求項5において、前記処理が、前記複合FM信号を時間ド
    メインにおけるデジタルサンプルへデジタル化し、前記時間ドメインサンプルを
    周波数ドメインサンプルへ変換し且つ前記周波数ドメインサンプルをフィルタし
    且つ夫々のECG波形に関連するグループへ分離し,且つ前記周波数ドメインサ
    ンプルのヒルベルト変換したものを時間ドメインへ変換させることを包含してい
    る方法。
  7. 【請求項7】 遠隔ステーションにおいて派生されたECG波形を中央ステ
    ーションにおいて再生する方法において、 遠隔ステーションにおいてECG波形を派生し、 前記ECG波形に従って搬送波を周波数変調して前記ECG波形に関連し且つ
    ゼロ交差を具備するFM信号を発生し、 前記FM信号を中央ステーションへ送信し、 前記中央ステーションにおいて前記FM信号を処理して前記FM信号のゼロ交
    差よりも実質的により頻繁に発生する時間においてその局所的周波数を推定し、 前記推定した局所的周波数を使用して前記中央ステーションにおいて前記EC
    G波形を再生する、 ことを包含している方法。
  8. 【請求項8】 請求項7において、前記局所的周波数が少なくとも前記搬送
    波周波数に等しい割合で推定される方法。
  9. 【請求項9】 請求項7において、前記局所的周波数が相次ぐゼロ交差の間
    において少なくとも数百のタイムスロットにおいて推定される方法。
  10. 【請求項10】 請求項7において、前記処理が前記FM信号のデジタル化
    したものを時間ドメインと周波数ドメインの両方において処理することを包含し
    ている方法。
  11. 【請求項11】 請求項7において、前記処理が前記FM信号の1つのバー
    ジョンをその位相シフトさせたバージョンと比較して前記FM信号の相次ぐサン
    プルに対する局所的位相差を推定し、且つ前記推定した位相差を使用して前記F
    M信号を再生することを包含している方法。
  12. 【請求項12】 ECG波形で周波数変調された搬送波を有するFM信号か
    らECG波形を再生する方法において、 前記FM信号を前記FM信号の局所的振幅を表わす第一時間ドメインデジタル
    値へデジタル化し、 前記時間ドメインデジタル値を前記ECG波形に対する興味のある帯域内の夫
    々の周波数において前記FM信号の内容を表わす周波数ドメインデジタル値へ変
    換し、 前記周波数ドメインデジタル値の位相シフトしたバージョンを発生し、 前記周波数ドメインデジタル値及びその位相シフトをさせたバージョンを夫々
    第二及び第三時間ドメインデジタル値へ変換し、 前記第二及び第三時間ドメインデジタル値を使用して前記搬送周波数の少なく
    とも2倍に対応する前記FM信号における位置において局所的周波数を計算し、 前記局所的周波数を使用して前記ECG波形を再生する、 ことを包含している方法。
  13. 【請求項13】 請求項12において、前記周波数ドメインへの変換がFF
    T処理を使用することを包含している方法。
  14. 【請求項14】 請求項12において、前記位相シフトさせたバージョンの
    発生がデジタル形態においてヒルベルト変換を使用することを包含している方法
  15. 【請求項15】 請求項12において、前記ECG波形に対する前記興味の
    ある帯域以外の周波数をフィルタ除去することを包含している方法。
  16. 【請求項16】 構成信号のタイムスロットが複数個のECG波形の各々か
    らの情報を包含しているように複数個のFMエンコードされた搬送波をECG波
    形の夫々のものと結合させることによって形成した複合FM信号から複数個のE
    CG波形を再生する方法において、前記複合波形はゼロ交差を有しており、 前記複合信号をデジタル化し、 前記複合信号をデジタルドメインにおいて処理して前記ゼロ交差よりも実質的
    により頻繁に発生する時間においてその局所的周波数を推定し、 前記推定した局所的周波数を使用して前記ECG波形を再生する、 上記各ステップを有している方法。
  17. 【請求項17】 請求項16において、前記処理が、前記複合信号を時間ド
    メイン及び周波数ドメインの両方において処理することを包含している方法。
  18. 【請求項18】 請求項17において、前記処理が、前記複合信号の周波数
    ドメインバージョンをフィルタし,且つ夫々のECG波形の周波数帯域に対応す
    るその部分を分離する、ことを包含している方法。
  19. 【請求項19】 請求項18において、前記処理が、前記複合信号の1つの
    バージョンをその位相シフトしたバージョンと比較して局所的位相シフトを推定
    し、且つ前記推定した局所的位相シフトを前記局所的周波数の推定において使用
    する、ことを包含している方法。
  20. 【請求項20】 請求項19において、前記比較が、前記複合信号の1つの
    バージョンをそのヒルベルト変換したバージョンと比較することを包含している
    方法。
JP2001564657A 2000-03-03 2001-03-01 マルチチャンネルecg波形の電話伝送監視 Pending JP2003525675A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US09/518,546 US6377843B1 (en) 2000-03-03 2000-03-03 Transtelephonic monitoring of multi-channel ECG waveforms
US09/518,546 2000-03-03
PCT/US2001/006556 WO2001065999A2 (en) 2000-03-03 2001-03-01 Transtelephonic monitoring of multi-channel ecg waveforms

Publications (1)

Publication Number Publication Date
JP2003525675A true JP2003525675A (ja) 2003-09-02

Family

ID=24064410

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001564657A Pending JP2003525675A (ja) 2000-03-03 2001-03-01 マルチチャンネルecg波形の電話伝送監視

Country Status (9)

Country Link
US (1) US6377843B1 (ja)
EP (1) EP1259158B1 (ja)
JP (1) JP2003525675A (ja)
KR (1) KR20030011785A (ja)
AT (1) ATE360393T1 (ja)
AU (2) AU2001243352B2 (ja)
CA (1) CA2402070A1 (ja)
DE (1) DE60128072T2 (ja)
WO (1) WO2001065999A2 (ja)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040210273A1 (en) * 2003-04-18 2004-10-21 Medtronic, Inc. Presentation of patient image with implantable medical device information
US7218967B2 (en) * 2003-09-26 2007-05-15 Medtronic, Inc. System and method for real-time remote monitoring of implantable medical devices
US7996075B2 (en) * 2004-10-20 2011-08-09 Cardionet, Inc. Monitoring physiological activity using partial state space reconstruction
US7729753B2 (en) * 2006-03-14 2010-06-01 Cardionet, Inc. Automated analysis of a cardiac signal based on dynamical characteristics of the cardiac signal
US7782191B2 (en) * 2007-07-25 2010-08-24 Tomas Flores Portable alarm apparatus for warning persons
US8395498B2 (en) 2007-08-31 2013-03-12 Cardiac Pacemakers, Inc. Wireless patient communicator employing security information management
US9848058B2 (en) * 2007-08-31 2017-12-19 Cardiac Pacemakers, Inc. Medical data transport over wireless life critical network employing dynamic communication link mapping
US8200319B2 (en) * 2009-02-10 2012-06-12 Cardionet, Inc. Locating fiducial points in a physiological signal
US8319631B2 (en) 2009-03-04 2012-11-27 Cardiac Pacemakers, Inc. Modular patient portable communicator for use in life critical network
US8812841B2 (en) 2009-03-04 2014-08-19 Cardiac Pacemakers, Inc. Communications hub for use in life critical network
KR100931829B1 (ko) * 2009-06-02 2009-12-15 주식회사 두성기술 페이서 장치
US8301232B2 (en) * 2010-06-08 2012-10-30 Alivecor, Inc. Wireless, ultrasonic personal health monitoring system
US8509882B2 (en) 2010-06-08 2013-08-13 Alivecor, Inc. Heart monitoring system usable with a smartphone or computer
US9351654B2 (en) 2010-06-08 2016-05-31 Alivecor, Inc. Two electrode apparatus and methods for twelve lead ECG
US9275013B2 (en) 2012-03-16 2016-03-01 Qualcomm Incorporated System and method for analysis and reconstruction of variable pulse-width signals having low sampling rates
WO2014036436A1 (en) 2012-08-30 2014-03-06 Alivecor, Inc. Cardiac performance monitoring system for use with mobile communications devices
US20140107979A1 (en) * 2012-10-16 2014-04-17 New York University Precise ferquency-pattern analysis to decompose complex systems into functionally invariant entities
WO2014074913A1 (en) 2012-11-08 2014-05-15 Alivecor, Inc. Electrocardiogram signal detection
WO2014107700A1 (en) 2013-01-07 2014-07-10 Alivecor, Inc. Methods and systems for electrode placement
US9254092B2 (en) 2013-03-15 2016-02-09 Alivecor, Inc. Systems and methods for processing and analyzing medical data
US9247911B2 (en) 2013-07-10 2016-02-02 Alivecor, Inc. Devices and methods for real-time denoising of electrocardiograms
EP3079571A4 (en) 2013-12-12 2017-08-02 Alivecor, Inc. Methods and systems for arrhythmia tracking and scoring
US20150286779A1 (en) * 2014-04-04 2015-10-08 Xerox Corporation System and method for embedding a physiological signal into a video
US9839363B2 (en) 2015-05-13 2017-12-12 Alivecor, Inc. Discordance monitoring
US11998341B2 (en) 2019-03-11 2024-06-04 Man-Rim CHOI Apparatus, method, and computer-readable recording medium for measuring size of electrocardiography signal using electrocardiography signal using Hilbert transform

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4889134A (en) * 1988-03-25 1989-12-26 Survival Technology, Inc. Device for measuring multiple channels of heartbeat activity and encoding into a form suitable for simultaneous transmission over
GB2234411A (en) * 1989-07-03 1991-01-30 Marconi Instruments Ltd Integrated circuit for digital demodulation

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3886314A (en) * 1973-09-20 1975-05-27 Parallel Data Systems No hands voice instruction for EEG telemetering system
US3898984A (en) * 1974-02-04 1975-08-12 Us Navy Ambulatory patient monitoring system
DE3114063A1 (de) * 1981-04-07 1982-10-21 Licentia Patent-Verwaltungs-Gmbh, 6000 Frankfurt Empfangssystem
US4658831A (en) * 1984-06-18 1987-04-21 Pacific Communications, Inc. Telemetry system and method for transmission of ECG signals with heart pacer signals and loose lead detection
US4938229A (en) 1988-04-29 1990-07-03 Paceart Inc. Method and apparatus for transmitting data descriptive of electromagnetic waveform
US5333617A (en) * 1989-10-16 1994-08-02 Marquette Electronics, Inc. Cardiac monitoring method and apparatus
US5467773A (en) 1993-05-21 1995-11-21 Paceart Associates, L.P. Cardiac patient remote monitoring using multiple tone frequencies from central station to control functions of local instrument at patient's home
US5454007A (en) * 1993-09-24 1995-09-26 Rockwell International Corporation Arrangement for and method of concurrent quadrature downconversion input sampling of a bandpass signal
US5735285A (en) 1996-06-04 1998-04-07 Data Critical Corp. Method and hand-held apparatus for demodulating and viewing frequency modulated biomedical signals
US5764704A (en) * 1996-06-17 1998-06-09 Symmetricom, Inc. DSP implementation of a cellular base station receiver
WO1998024212A1 (en) * 1996-11-29 1998-06-04 Micromedical Industries Limited Telemedicine system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4889134A (en) * 1988-03-25 1989-12-26 Survival Technology, Inc. Device for measuring multiple channels of heartbeat activity and encoding into a form suitable for simultaneous transmission over
GB2234411A (en) * 1989-07-03 1991-01-30 Marconi Instruments Ltd Integrated circuit for digital demodulation

Also Published As

Publication number Publication date
EP1259158A2 (en) 2002-11-27
ATE360393T1 (de) 2007-05-15
US6377843B1 (en) 2002-04-23
EP1259158B1 (en) 2007-04-25
DE60128072T2 (de) 2007-12-27
KR20030011785A (ko) 2003-02-11
DE60128072D1 (de) 2007-06-06
CA2402070A1 (en) 2001-09-13
EP1259158A4 (en) 2004-06-23
WO2001065999A2 (en) 2001-09-13
AU4335201A (en) 2001-09-17
AU2001243352B2 (en) 2005-10-27
WO2001065999A3 (en) 2002-02-14

Similar Documents

Publication Publication Date Title
JP2003525675A (ja) マルチチャンネルecg波形の電話伝送監視
AU2001243352A1 (en) Transtelephonic monitoring of multi-channel ECG waveforms
US3810102A (en) System for transmission and analysis of biomedical data
CA2079260C (en) Method and apparatus for automatically identifying a program including a sound signal
US5617871A (en) Spread spectrum telemetry of physiological signals
US5689529A (en) Communications method and apparatus for digital information
US6480820B1 (en) Method of processing auditory data
CN101641891B (zh) 调制装置、解调装置及声信号传输方法
EP0615197A1 (en) Methods for enhancements of HRV and late potentials measurements
EP0586286A1 (fr) Système de transmission de données médicales, sans fil
GB2292506A (en) Automatically identifying a program including a sound signal
US4938229A (en) Method and apparatus for transmitting data descriptive of electromagnetic waveform
US6654633B2 (en) Mobile neurological signal data acquisition system and method
US4424815A (en) Stethoscope-type recording system with stretched-in-time playback
US3913567A (en) Electrocardiac information monitoring apparatus
JPS61126836A (ja) Pcm伝送方式およびその装置
Alexander et al. Computer diagnosis of electrocardiograms. I. Equipment
WO1998034542A2 (en) Stethoscope communication and remote diagnosis system
Taylor et al. Data compression of fetal Doppler ultrasound audio signals using zero-crossings analysis
JPH0346730Y2 (ja)
US20070041529A1 (en) Synchronous parallel acoustic transmission in transtelephonic medical monitors
EP1037440A3 (en) Null symbol detecting device
Lawrence Role of Synthetic Speech in Speech Research
Gold et al. Bandpass Compressor
JPH0433176B2 (ja)

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080227

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20101125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101207

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20110304

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20110311

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20110705