JP2001500777A - Ecgによってトリガされるmr撮像方法及び装置 - Google Patents

Ecgによってトリガされるmr撮像方法及び装置

Info

Publication number
JP2001500777A
JP2001500777A JP11509576A JP50957699A JP2001500777A JP 2001500777 A JP2001500777 A JP 2001500777A JP 11509576 A JP11509576 A JP 11509576A JP 50957699 A JP50957699 A JP 50957699A JP 2001500777 A JP2001500777 A JP 2001500777A
Authority
JP
Japan
Prior art keywords
magnetic field
human body
measuring
signal
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.)
Granted
Application number
JP11509576A
Other languages
English (en)
Other versions
JP4104674B2 (ja
Inventor
イー フィッシャー,ステファン
エイ ウィックライン,サミュエル
エイチ ロレンズ,クリスティン
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.)
Koninklijke Philips NV
Original Assignee
Philips Electronics NV
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 Philips Electronics NV filed Critical Philips Electronics NV
Publication of JP2001500777A publication Critical patent/JP2001500777A/ja
Application granted granted Critical
Publication of JP4104674B2 publication Critical patent/JP4104674B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5673Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker
    • 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/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Physiology (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Power Engineering (AREA)
  • Pulmonology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Cardiology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 本発明は、人体のECGデータを測定する段階と、MR信号の測定を人体の周期的運動と同期するためにECGデータからの捕捉周期を決定する段階とを含む磁気共鳴(MR)によって静磁場の中に配置された人体の一部を撮像する方法であって、ECGデータを測定する段階はベクトルカルジオグラムを測定する段階を含む方法に関する。ベクトルカルジオグラムの方向性情報から、捕捉周期はMR信号の測定と同期するよう決定される。本発明法の利点は、ベクトルカルジオグラムの外乱に対して減少された感度を有する。ECGトリガ動作が得られることである。

Description

【発明の詳細な説明】 ECGによってトリガされるMR撮像方法及び装置 本発明は、人体のECGデータを測定する段階と、 MR信号の測定を人体の周期的運動と同期させるためにECGデータから捕捉 周期を決定する段階と、捕捉周期中にRF場及び傾斜磁場を含むMRシーケンス を発生する段階とを含み、MR信号から人体の一部の画像を再構成する段階を更 に含む、静磁場の中に配置された人体の一部を磁気共鳴(MR)によって撮像す る方法に関する。 本発明は更に、係る方法を実行するためのMR装置に関する。 本発明は更に、MR装置における使用に適したECG装置に関する。 係る方法及び装置は、米国特許第5,526,813号より既知である。 本願の文脈において、k空間は、MR信号が軌道に沿って測定され、測定され た値が人体の画像の逆フーリエ変換された値を生ずる空間周波数ドメインである と理解されることに注意すべきである。k空間中の軌道は、核スピンの励起から MR信号の測定の実際の時点までの時間間隔に亘る一時傾斜磁場の時間積分によ って決定される。更に、傾斜磁場は、定常磁場の上に重畳され、3つの夫々直交 方向に静磁場の傾斜を生ずる一時磁場であると理解される。概して、第1の方向 の傾斜磁場は読取り傾斜磁場と称され、第2の方向の傾斜磁場は位相符号化傾斜 磁場と称され、第3の方向の傾斜磁場はスライス符号化傾斜磁場と称される。 既知の方法では、MR画像の質はMR信号の捕捉を人体の心拍と同期させるこ とによって改善される。ここでMR信号の捕捉周期は測定されたECGデータ、 例えばR波から決定される。捕捉周期の 誤った決定を減少するために、ECGデータはフィルタされる。フィルタリング によって、人体から獲得されたECGデータから雑音が減少される。人体の領域 のMR画像は次に連続する同期された捕捉周期中に測定された人体の一部からの 1組のMR信号から再構成される。既知の方法の欠点は、ECGデータのフィル タリングは信頼性がないことがあり、結果としてフィルタされたECGデータか ら誤った捕捉周期が決定されうることである。 本発明は特に、誤った捕捉周期の決定を減少することを目的とする。このため 、本発明による方法は、上記ECGデータを測定する段階はベクトルカルジオグ ラムを測定する段階を含むことを特徴とする。本願の文脈において、ベクトルカ ルジオグラムは、心房心室脱分極及び心室脱分極の間の平均瞬間空間心臓ベクト ルの時間経過の登録として理解されるべきである。測定空間の中で形成される瞬 間ベクトルの特徴ループはP、QRS及びT軌道又はループとしてラベル付けさ れる。 本発明は、測定されたECGデータは捕捉周期中に発生されたRFパルスによ って妨害されるだけでなく、ECGデータは傾斜を切り替えることによる人体上 の電極からECG装置まで導線の中の誘導、及び静磁場が存在する場合は血流に よって誘導されるECG信号自体の中の雑音の誘導によって妨害されるという概 念に基づく。更に、ベクトルカルジオグラムを測定することにより、捕捉周期の 決定においてベクトルカルジオグラムの方向性情報を組み込むことが可能である 。測定されたベクトルはある平面の中の所定のループを記述し、一方アーティフ ァクトは異なる平面及び方向の中で優性である。このようにして、例えば心拍の R波を決定することが可能である。捕捉周期は次に、例えばR波上でトリガされ うる。更に、ベクトルカルジオグラムの性質は、カルジオグラムの生理パラメー タを傾斜切替雑音から弁別するために適用されうる。例えば、静磁場上に適用さ れた切り替えられた傾斜磁場は全ての導線の中に同様 の電圧差を誘導する。このように、測定されたベクトルは傾斜磁場を変化させる ことによって変更され、一方測定されたベクトルの方向は変化されない。この特 性は、例えば傾斜切替雑音からQRS複合体を弁別するために使用されうる。本 発明による方法の更なる利点は、MR検査中、人体のECGの診断値が更に改善 されうることである。 本発明による方法の特定の変形は、上記ベクトルカルジオグラムを測定する段 階は、直交する導線の組によって人体に複数の電極の配置を貼り付ける段階を含 むことを特徴とする。直交導線系は、それ自体としてC.Lentner,Geigy Scienti fic Tables,Heart and Circulation,45,1990 より既知である。直交導線系は 、相互に直交し、人体の横軸、サジタル軸及び垂直軸に対して直交する導線を含 み、導線の振幅はベクトルの観点からは等しくなり、導線は、単一の点に対して だけでなく、起電力が発生される心臓の中の全ての点に対して同一の強さ及び方 向を有する。 直交導線の組を使用することにより、心臓の分極及び脱分極の記述のために単 純なモデルが適用されうる。結果として、測定されたベクトルカルジオグラムは 直交ベクトル成分へ分解されうる。ECGの生理パラメータ、例えばP波、QR S複合体及びT波は、投影面に特徴ループを形成する。3つの直交ベクトル成分 のうちの2つによって生じた平面上のこれらの3次元のループの投影は、パラメ ータ化されうる。獲得されたパラメータは、静磁場の存在なしで測定された人体 のベクトルカルジオグラムのパラメータと比較されうる。更に、血流の相互作用 によって生じた静磁場の中の測定されたベクトルカルジオグラムのアーティファ クト及び静磁場は、生理パラメータのループのパラメータと比較して異なるパラ メータによって記述されるループによって表わされる。 本発明による方法の更なる変形は、上記複数の電極の配置は、補正された直交 する導線の組を含むことを特徴とする。 本発明による方法の更なる変形は、上記ベクトルカルジオグラムの第1の弁別 関数の第1の実際値を決定する段階を更に含み、上記捕捉周期は第1の実際値及 び第1の所定の基準関数から決定されることを特徴とする。これらの段階の結果 として、捕捉周期は、例えば心臓の電気軸とベクトルカルジオグラムの瞬間ベク トルとの間の角度から獲得されうる。他の可能性は、捕捉周期がベクトルカルジ オグラムの瞬間ベクトルと心臓の電気軸とのスカラー積から獲得されることであ る。この方法は、R波と望ましくない血流に関するアーティファクトの比率を増 加させるための簡単な方法を提供する。捕捉周期がベクトルカルジオグラムの2 つの隣接する点によって記述される平面の法線ベクトルの絶対値から決定される こともまた可能である。この決定は、ECGのR波の間は、ベクトルカルジオグ ラムの習慣ベクトルの速度は、ECGの他の周期中のベクトルカルジオグラムの 瞬間ベクトルの速度よりも大きいことに基づく。 基準関数は、第1の実際値の最大値の割合、例えば75%に基づく一定の閾値 でありうる。他の可能性は、基準関数の閾値の代わりに、基準関数がファジー論 理を含むことである。ファジー論理は、一定の閾値の代わりに確率によって動作 し、即ち、ベクトルカルジオグラムの瞬間ベクトルが心臓の電気軸の回りの球状 の領域の中であれば、R波が検出される確率は高い。この規準は、他の規準、例 えば局部最大検出と組み合わせられ得る。即ち、ベクトルカルジオグラムの隣接 する点間の差が大きい局部最大では、これらの点がR波を表わすという大きな確 度がある。他の規準は、次のR波が低い可能性は、第1のR波の検出の直ぐ後、 及び進む時と共に増加するということに基づく。例えば、2つの連続するR波の 間の標準時間間隔の80%の後では、次のRのための確率は1となる。 本発明による変形例は、基準ベクトルカルジオグラムを獲得する段階と、基準 ベクトルカルジオグラムから第1の弁別関数の基準値を決定する段階とを更に含 み、上記捕捉周期は第1の実際値、基準 値及び第2の所定の基準関数の組み合わせから決定されることを特徴とする。こ れらの段階の結果、第1の弁別関数の実際値は、静磁場の外側に配置される人体 について決定された第1の弁別関数の規準値によって重みづけされる。 本発明による方法の更なる実施例は、上記ベクトルカルジオグラムの第2の弁 別関数の第2の実際値を決定する段階を更に含み、上記捕捉周期は上記第1の実 際値、第2の実際値及び第2の所定の基準関数の組み合わせから決定されること を特徴とする。 結果として捕捉周期は、異なる弁別関数、例えばベクトルカルジオグラムの瞬 間ベクトルと心臓の電気軸との間のスカラー積の弁別関数及びベクトル図の瞬間 的なベクトルと心臓の電気軸との間の確度の弁別関数等の実際値の加重された組 み合わせによって決定されうる。 本発明の方法の更なる実施例は、静磁場のアイソセンタに関して複数の電極か ら選択された1つの電極の方向を決定するために、2つの直交する傾斜磁場を周 期的に印加することによって回転磁場を発生する段階と、磁場を回転させること によって選択された電極及び人体からなるコイルの中に誘導された方向信号を測 定し、測定された方向信号から上記選択された電極の静磁場のアイソセンタに関 する方向を決定する段階とを更に有することを特徴とする。これは、電極の方向 を静磁場のアイソセンタに対して決定するための容易な方法である。アイソセン タは、印加される傾斜磁場の影響が略ゼロの静磁場の中の位置である。電極位置 は、傾斜切替情報を使用する適合フィルタリング及び心臓の双極子ベクトルモデ ルの計算のために重要である。回転する磁場の周波数は、ECGのフィルタリン グ及び抑制が容易であるよう選択されうる。測定された信号の位相は、アイソセ ンタに関する電極の方向を決定するために使用されうる。 本発明による方法の更なる変形は、上記複数の電極から選択され た電極の傾斜磁場の方向に沿った位置を決定するために、傾斜磁場の方向に沿っ た位置から線形に依存する振幅を有する第1の傾斜磁場を発生し、上記選択され た電極及び人体からなるコイルの中の第1の誘導電圧を測定する段階と、傾斜磁 場の方向に沿った位置の2次又はより高次の関数によって変化する振幅を有する 第2の傾斜磁場を発生し、コイルの中の第2の誘導電圧を測定する段階と、第1 の誘導電圧と第2の誘導電圧との比から傾斜磁場の方向に沿った電極の位置を決 定する段階とを更に有することを特徴とする。 本発明はまた、静磁場の中に配置されたRF場及び傾斜磁場を含むMRシーケ ンスを受ける人体の一部を磁気共鳴(MR)によって撮像する方法であって、人 体のベクトルカルジオグラムを測定する段階と、ベクトルカルジオグラムの第1 の弁別関数の第1の実際値を決定する段階と、第1の実際値及び第1の所定の基 準関数から生理パラメータを決定する段階とを含む方法に関する。結果として、 発生されたRF場、印加された傾斜及び存在する静磁場による望ましくないアー ティファクトを実質的に含まずに、診断のための人体のECGが獲得されうる。 本発明によるMR装置の実施例は、上記ECG信号を測定するための手段は、 ベクトルカルジオグラムを測定するための手段を含むことを特徴とする。 本発明によるMR装置における使用に適したECG装置の実施例は、人体のベ クトルカルジオグラムを測定する手段と、ベクトルカルジオグラムの第1の弁別 関数の第1の実際値を決定する手段とを含むことを特徴とする。 本発明の上述及び他の面は以下説明される実施例を参照して明らかとなろう。 図中、 図1はMR装置の実施例を示す図であり、 図2は既知の方法によって測定された2つのエコーカルジオグラ ムを示す図であり、 図3はフランク導線組の配置を概略的に示す図であり、 図4は静磁場が存在しないときの健康なボランティアの第1のベクトルカルジ オグラムを示す図であり、 図5は静磁場が存在するときの健康なボランティアの第2のベクトルカルジオ グラムを示す図であり、 図6は傾斜座標系x,y,z及びコイル60の位置を示す図であり、 図7はフィルタされた電極信号及び回転する傾斜場を時間の関数として示す図 である。 図1は、MR装置の実施例を示す図である。MR装置1は、定常磁場を発生す るための第1の磁石系2を含む。図示される座標系のz方向は、磁石系2の中の 定常磁場の方向に対応する。MR装置はまた、z方向に向けられた一時磁場を発 生し、x,y及びz方向に夫々傾きを発生するための第2の磁石系3を含む。こ れらの方向は座標系の主方向と一致する必要はない。更に、本願では、x方向、 y方向及びz方向に傾きを有する一時傾斜磁場は、夫々、読取り傾斜磁場、位相 符号化傾斜磁場及びスライス選択傾斜磁場と称される。電力供給手段4は、第2 の磁石系3へ電力を供給する。磁石系2は例えば人体の一部といった検査される べき対象7の一部を収容するために十分大きい検査空間を囲む。RF送信器コイ ル5はRF磁場を発生するために使用され、送信器/受信器回路9を通じてRF 源及び変調器6へ接続される。RF送信器コイル5は検査空間の中の人体7の一 部の回りに配置される。MR装置はまた、送信器/受信器回路9を通じて信号増 幅及び変調ユニット10へ接続される。受信器コイル及びRF送信器コイル5は 同一のコイルでありうる。制御ユニット11は、RFパルス及び一時傾斜磁場を 含むRF撮像シーケンスを発生するために、変調器6及び電力供給手段4を制御 する。検査空間の中に配置された人体の一部の中に核スピンを励起 した後、受信器コイル5はMR信号を受信する。そこから導出される位相及び振 幅はサンプリングされ、更に増幅及び変調ユニット10の中で処理される。画像 再構成ユニット12は画像を形成するよう与えられたMR信号を処理する。画像 は例えばモニタ13上に表示される。制御ユニット11はまた画像再構成ユニッ ト12を制御する。更に、MR装置は、MR画像シーケンスの発生を、人体の周 期的な運動、例えばMR装置の中に配置された人体の心拍と同期させるためにE CG装置14を含む。ECG信号を測定するために、ECG装置は既知の方法で 、人体に取り付けられた電極15を通じて接続される。例えばDas Elektrokardi ogram by T.Klinge,Thieme Verlag, 1987に記載される。 ECG装置は、例えば測定されたECGにR波が存在するときに制御ユニット のためのトリガ信号を発生する。しかしながら、静磁場の存在はECGの中にア ーティファクトを発生する。静磁場が存在する場合に測定されたECGと、静磁 場が存在しない場合に測定されたECGとの差は、図2を参照して明らかとなろ う。 図2は、健康なボランティアの第1及び第2のECGを示す図である。第1の ECGはMR装置の静磁場の存在なしに記録される。第2のECGはMR装置の 例えば1.5Tの静磁場が存在するときに記録される。静磁場が存在するときに 記録された第2のECG21は、第1のECG20と比較して幾つかの追加的な 信号ピークを含む。追加的な信号ピークから誤った捕捉周期が決定されうる。 誤った捕捉周期を減少するため、本発明による方法では、捕捉周期は人体のE CGの代わりに人体から獲得されるベクトルカルジオグラムから導出される。そ れによりECG装置は第1の基準ベクトルカルジオグラムと、人体がMR装置の 中に配置される前の心臓の電気軸とを測定する。ベクトルカルジオグラムは、補 正された直交導線組によって人体に複数の電極を配置した後に測定される。例え ば、補正された直交導線組は4つの電極を有する。4つの電極の信 号は3つの信号Xf,Yf,Zfへ変換される。3つの信号Xf,Yf,Zfは、導線 が相互に直交し、人体の横軸、サジタル軸及び垂直軸に直交する直交導線系を記 述する。座標系Xf,Yf,Zfの向きは、MR装置の座標系x,y,zとは異な ることに注意すべきである。 図3は、フランク導線組による電極の配置を示す図である。図3は、A,C, E,F,H,I及びMによって示される7つの電極位置を示す。フランク導線組 は、1956年Circulation 13号、第737頁で発表されたE.Frankによ る文献から既知である。或いは、4電極配置が使用されうる。更に、ベクトルカ ルジオグラムは、3つの標準化された平面、即ち、前面、縦面及び左又は右サジ タル面で記録される。これらの平面は、直交導線X,Y及びZによって定義され る。XY平面は前面を、XZ平面は縦面を、ZY平面はサジタル平面を記述する 。時間経過を視覚化するために、ベクトルループは例えば1又は2秒毎に中断さ れうる。本発明による方法を説明するために、図4及び図5中、静磁場が存在す る場合及び静磁場が存在しない場合の人体のベクトルカルジオグラムの例を図示 する。 図4は、静磁場が存在しない場合の健康なボランティアのQRSループの前面 即ちXY面中の第1の投影40を示す図である。図5は、例えば体軸に平行な1 .5Tの静磁場が存在する場合の同一の健康なボランティアのQRSループの前 面中の第2の投影50を示す図である。QRSループの第2の投影50はよく保 存され、一方血流の相互作用及び静磁場によって生じた他のループ、例えばルー プ51はQRSループの投影50の平面及び方向とは異なる平面及び方向で観察 されうる。 本発明による方法では誤った捕捉周期の数を減少させるために、QRS波と望 ましくない血流に関連するアーティファクトとの間の差を拡大する弁別関数が導 入される。QRS波と望ましくない血流 に関連するアーティファクトとの間の決定は、弁別関数の実際値及び所定の基準 関数に基づく。本発明の第1の方法では、第1の弁別関数は、例えばベクトルカ ルジオグラムの瞬間ベクトルと心臓の電気軸とのスカラー積、即ち式で表わすと 、となる。式において、ψは瞬間ベクトルと心臓の電気軸を表わすベクトルRとの 間の角度を表わす。このようにして、例えばR波とベクトルカルジオグラムの中 の望ましくない血流に関連するアーティファクトとの比率は、1から1.5へ改 善されうる。基準関数の例は、例えば弁別関数s1(t)の最大振幅の75%の 一定の閾値である。 R波と血流に関連するアーティファクトとの間の分化を改善するために、R波 と望ましくない血流アーティファクトとの間の比率は、スケーリングされた関数 s4を使用することによって増加されえ、s1(t)はECGとベクトルRとのス カラー積を表わし、Rは磁石の外側で獲得された心臓のRベクトルを表わすとす ると、 によって定義される。信号s4(t)は3倍に増加されたR波を示す。例えば振 幅、値4の閾値及び1次導関数に基づく単純なR波検出アルゴリズムは、s4( t)に対して適用されうる。 本発明によれば上記の方法は、例えばMR装置1のECG装置14において実 施されうる。人体の測定されたベクトルカルジオグラムの弁別関数の実際値及び 一定の閾値に基づいて、ECG装置は捕捉周期を開始させるために制御ユニット 11へ送信されるトリガ信号を発生する。 更に、本発明の第2の方法では、第2の弁別関数は、例えば、ベクトルカルジ オグラムの瞬間ベクトルの振幅Vecg (t)を瞬間ベクトルと心臓の電気軸との 間の角度αの関数fwの重み付けとして 定義される。式で表わすと、 s2(t)=Vecg(t).fw(α)(2) となる。係る関数fwの例は、によって与えられる。第2の基準関数は、再び例えば弁別関数s2(t)の最大 振幅の75%の一定の閾値に基づく。 本発明による第3のより進んだ方法では、第3の弁別関数は、例えばベクトル カルジオグラムの3つの連続する点によって定義される平面の法線ベクトルとし て定義される。式で表わすと、 v(t)=vx,vy,vz=Vecg(t)-Vecg(t-dt)x(Vecg(t)-Vecg(t+dt))(3) となる。定義される法線ベクトルは、ベクトルカルジオグラム軌道の速度v(t )としてみなされうる。するとQRS波と望ましくない血流アーティファクトと の間の決定は、R波の間、ベクトルカルジオグラム軌道の速度v(t)は他の周 期の間の速度よりも大きいという観察に基づく。R波の持続時間は約60msで ある。第3の基準関数は、再び例えばベクトルカルジオグラム軌道の最大速度の 90%の一定の閾値に基づきうる。 本発明の更なる方法では、v(t)の絶対値は、静磁場が存在するときに人体 から獲得されたベクトル図の軌道の速度と、静磁場が存在しないときの同一の人 体からのベクトルカルジオグラム軌道の基準速度の速度とのスカラー積によって 重み付けされうる。 誤った捕捉の数を更に減少させるために、本発明の第4の方法では、第4の弁 別関数s3は、第1、第2及び第3の弁別関数及びそれらの夫々の閾値に基づい て例えば決定D1,D2,D3の重み付けされた組み合わせに基づいて定義される 。例えば、 s3=c1D1(s1)+c2D2(s2)+c3D3(v(t)) において、例えばc1=c2=c3=0.5である。決定は、値s3及び一定の閾値 1に基づいて行われうる。 或いは、本発明による他の方法では、ベクトルカルジオグラムのR波と望まし くない血流に関するアーティファクトとの間の決定は、ファジー論理によって定 義される。ファジー論理は、例えば、3つの確率関数P1 (t),P2 (t), P3 (t)の組み合わせに基づき、ここで第1の確率関数P1 (t)はベクトル カルジオグラムの瞬間ベクトルの終端点がベクトルカルジオグラムの電気軸Rの 半径dの球の中にある確率P1 (t)として定義される。式で表わすと、 となる。第2の確率関数P2 (t)は局部極値、例えば、ベクトルカルジオグラ ムの2つの連続する点の間の局部極値として定義され、確率P2 (t)はR波の 間は高いと観察される。式で表わすと、 となる。第3の確率関数P3 (t)は最近検出されたR波の後の次の続くR波の 検出の確率であり、最近検出されたR波の直後の観察では、確率P3 (t)は低 く、心臓の平均周期時間Tの80パーセントの間隔tが経過した後は、確率P3 (t)は高い。式で表わすと、 となる。R波と望ましくない血流アーティファクトとの間で決定を行う他の可能 性は、例えば、ベクトルカルジオグラムの中の所定のパターンを認識することを 学習しうるニューラルネットに基づきうる。 電極場所に関する知識は、傾斜切替情報を使用する適応フィルタリング及びベ クトル双極子モデルのために重要である。 MR装置の中で人体に配置された補正された直交導線組の4つの電極の場所を 決定するために、本発明によれば静磁場上に回転傾斜磁場を印加するために傾斜 増幅器へ送信される制御信号を含む専用 シーケンスを発生するために、制御ユニットが更に配置されている。電極の場所 を決定するための方法を図6を参照して説明する。 図6は、傾斜座標系x,y,z及びコイル60の位置を示す図である。コイル 60は、導線、人体及びECG装置によって形成される。制御ユニットは、傾斜 座標系が例えばx,y平面に回転磁場を発生するよう傾斜増幅器への制御信号を 発生する。回転磁場B(t)は、B0は静磁場を表わし、Gは傾斜磁場の振幅( mT/m)を表わし、ωは回転磁場の周波数を表わすとすると、 によって表わされうる。回転傾斜磁場は、コイル60の中に信号u(t)を誘導 し、Sはコイルの表面を表わすとすると、信号u(t)は、 によって与えられる。誘導された信号は、コイル60の接続61,62上で測定 されうる。更に、回転磁場の周波数ωは、測定された信号の中のECG情報が容 易に例えば100Hzへ抑制されうるようにされる。測定された信号の中のEC G情報は、例えば回転傾斜磁場によって誘導された電極信号を通過させる高域通 過フィルタによって抑制されうる。フィルタリングされた信号の位相ψは、MR 装置の中の傾斜磁場のアイソセンタに関する測定された電極の方向を決定するた めに使用されうる。本発明のこの面を図7を参照して説明する。 図7は、相対的な位置x=y=1における電極のフィルタリングされた電極信 号を示す図である。細実線72は測定されフィルタリングされた電極信号を表わ し、破線71は時間に関する第1の正弦関数としてy方向に傾斜磁場の振幅を表 わし、太実線70は、第1の正弦関数に対してπ/2の位相シフトを有する時間 に関する第2の正弦関数としてx方向に傾斜磁場の振幅を表わす。測定された信 号は夫々、x方向及びy方向の両方の誘導された電圧を表わす。例 えばx方向の傾斜磁場と誘導された電圧との間の位相差ψから、アイソセンタに 関する電極の場所の方向は、 によって決定されうる。 更に、フィルタリングされた信号の振幅は電極の場所からアイソセンタまでの 距離に比例するため、正弦信号の振幅は、電極の場所からアイソセンタまでの相 対的な距離を決定するために使用されうる。ECG導線及び人体によって画成さ れるコイルの寸法は正確に知られていないため、座標系x,y,zの中心の絶対 距離は受信された正弦信号からは決定されえない。 本発明による傾斜座標系の中の絶対位置r(rx,ry,rz)を決定するため に、線形及びより高次の傾斜磁場についてx,y,z実験が行われた。例えば図 6に示される座標系x,y,zのコイル60の距離rを決定するために、傾斜座 標系は、第1の実験においてx方向に線形傾斜磁場を発生する。最大傾斜強さの 立上り時間dtに亘って、U1は接続61,62上の誘導された電圧を表わし、 Sは磁場の中のコイル表面を表わし、Gは傾斜磁場の振幅を表わし、rxはx方 向に沿った位置の成分を表わすとすると、コイル60の電圧は、 として与えられる。更に、傾斜座標系は、第2の実験では高次の傾斜磁場、例え ば2次傾斜磁場を発生する。2次磁場の立上り時間dtに亘って、コイル60の 電圧U2は、 として与えられる。電圧U1の電圧U2に対する比率は、傾斜座標系におけるコイ ル60の距離rxを生ずる。すると距離rxは、 として与えられる。座標系x,y,zにおいてコイル60の位置(rx,ry,rz )を決定するために、上述の実験はy方向及びz方向に繰り返される必要があ る。

Claims (1)

  1. 【特許請求の範囲】 1. 人体のECGデータを測定する段階と、 MR信号の測定を人体の周期的運動と同期させるためにECGデータから捕捉 周期を決定する段階と、該捕捉周期中にRF場及び傾斜磁場を含むMRシーケン スを発生する段階とを含み、該MRシーケンスに亘って獲得されたMR信号から 人体の一部の画像を再構成する段階を更に含む、静磁場の中に配置された人体の 一部を磁気共鳴(MR)によって撮像する方法であって、 上記ECGデータを測定する段階はベクトルカルジオグラムを測定する段階を 含むことを特徴とする方法。 2. 上記ベクトルカルジオグラムを測定する段階は、直交する導線の組によっ て人体に複数の電極の配置を付ける段階を含むことを特徴とする請求項1記載の 方法。 3. 上記複数の電極の配置は、補正された直交する導線の組を含むことを特徴 とする請求項2記載の方法。 4. 上記ベクトルカルジオグラムの第1の弁別関数の第1の実際値を決定する 段階を更に含み、上記捕捉周期は該第1の実際値及び第1の所定の基準関数から 決定されることを特徴とする請求項1乃至3のうちいずれか1項記載の方法。 5. 基準ベクトルカルジオグラムを獲得する段階と、該基準ベクトルカルジオ グラムから第1の弁別関数の基準値を決定する段階とを更に含み、上記捕捉周期 は第1の実際値、基準値及び第2の所定の基準関数の組み合わせから決定される ことを特徴とする請求項4記載の方法。 6. 上記ベクトルカルジオグラムの第2の弁別関数の第2の実際値を決定する 段階を更に含み、上記捕捉周期は上記第1の実際値、該第2の実際値及び第2の 所定の基準関数の組み合わせから決定されることを特徴とする請求項4記載の方 法。 7. 上記捕捉周期は、第1の実際値、第2の実際値、基準値及び第3の所定の 基準関数の組み合わせから決定されることを特徴とする請求項5及び6記載の方 法。 8. 静磁場のアイソセンタに関して複数の電極から選択された1つの電極の方 向を決定するために、 2つの直交する傾斜磁場を周期的に印加することによって回転磁場を発生する 段階と、 磁場を回転させることによって選択された電極及び人体からなるコイルの中に 誘導された方向信号を測定し、該測定された方向信号から上記選択された電極の 静磁場のアイソセンタに関する方向を決定する段階とを更に有することを特徴と する請求項2乃至7のうちいずれか1項記載の方法。 9. 上記複数の電極から選択された電極の傾斜磁場の方向に沿った位置を決定 するために、 傾斜磁場の方向に沿った位置から線形に依存する振幅を有する第1の傾斜磁場 を発生し、上記選択された電極及び人体からなるコイルの中の第1の誘導電圧を 測定する段階と、 傾斜磁場の方向に沿った位置の2次又はより高次の関数によって変化する振幅 を有する第2の傾斜磁場を発生し、コイルの中の第2の誘導電圧を測定する段階 と、 第1の誘導電圧と第2の誘導電圧との比から傾斜磁場の方向に 沿った電極の位置を決定する段階とを更に有することを特徴とする、請求項2乃 至7のうちいずれか1項記載の方法。 10. 略均一な定常磁場の中に配置された人体の領域の画像の捕捉用のMR装 置であって、 (a)定常磁場を発生するための磁石と、 (b)一時傾斜磁場を発生する手段と、 (c)RFパルスを発生する手段と、 (d)MR信号を発生する手段と、 (e)ECG信号を測定する手段と、 (f)人体の領域の画像を決定するよう測定されたMR信号を処理する処理手段 と、 (h)傾斜磁場を発生するための手段及びRFパルスを発生するための手段のた めの制御信号を発生するためにトリガ信号によって同期され、夫々がRF場及び 傾斜磁場を含む複数のシーケンスが発生されるよう配置された制御ユニットと、 (i)MR信号の測定を心臓の周期的運動に同期させるために測定されたECG 信号から上記制御ユニットのためのトリガ信号を決定する同期手段とを含むMR 装置であって、 上記ECG信号を測定するための手段は、ベクトルカルジオグラムを測定する ための手段を含むことを特徴とするMR装置。 11. 請求項10記載のMR装置における使用に適したECG装置であって、 人体のベクトルカルジオグラムを測定する手段と、 ベクトルカルジオグラムの第1の弁別関数の第1の実際値を決定する手段と、 該第1の実際値及び第1の所定の基準関数からMR装置用のトリガ信号を発生 するための発生手段とを含むことを特徴とするECG 装置。
JP50957699A 1997-07-23 1998-07-13 Ecgによってトリガされるmr撮像方法及び装置 Expired - Lifetime JP4104674B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US08/898,830 US5987348A (en) 1997-07-23 1997-07-23 ECG triggered MR imaging method and apparatus
US08/898,830 1997-07-23
PCT/IB1998/001062 WO1999004688A1 (en) 1997-07-23 1998-07-13 Ecg triggered mr imaging method and apparatus

Publications (2)

Publication Number Publication Date
JP2001500777A true JP2001500777A (ja) 2001-01-23
JP4104674B2 JP4104674B2 (ja) 2008-06-18

Family

ID=25410093

Family Applications (1)

Application Number Title Priority Date Filing Date
JP50957699A Expired - Lifetime JP4104674B2 (ja) 1997-07-23 1998-07-13 Ecgによってトリガされるmr撮像方法及び装置

Country Status (5)

Country Link
US (1) US5987348A (ja)
EP (1) EP0925019B1 (ja)
JP (1) JP4104674B2 (ja)
DE (1) DE69828372T2 (ja)
WO (1) WO1999004688A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011524770A (ja) * 2008-06-18 2011-09-08 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁気共鳴イメージングのための心電計及びそのための電極パッチ
JP2013000590A (ja) * 2011-06-17 2013-01-07 General Electric Co <Ge> 心電図検査におけるノイズ検出システム及び方法
JP2016093488A (ja) * 2014-11-06 2016-05-26 株式会社東芝 信号処理装置、撮像装置及び信号処理方法
US11058361B2 (en) 2014-11-06 2021-07-13 Canon Medical Systems Corporation Signal processing apparatus, imaging apparatus, and signal processing method

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6144201A (en) * 1997-12-26 2000-11-07 Kabushiki Kaisha Toshiba MR imaging utilizing ECG gating technique
US6198959B1 (en) * 1998-03-27 2001-03-06 Cornell Research Foundation Inc. Coronary magnetic resonance angiography using motion matched acquisition
WO2001046709A1 (en) * 1999-12-20 2001-06-28 Koninklijke Philips Electronics N.V. Mri apparatus with flexible gradient control
US7572231B2 (en) * 2000-01-31 2009-08-11 Pearlman Justin D Method of and system for signal separation during multivariate physiological monitoring
AU2001231265A1 (en) * 2000-01-31 2001-08-07 Justin D. Pearlman Multivariate cardiac monitor
US6501979B1 (en) 2000-03-09 2002-12-31 Koninklijke Philips Electronics N.V. Methods and devices for combined ECG and PPU controlled magnetic resonance imaging
WO2001067960A2 (en) 2000-03-16 2001-09-20 Cedara Software Corp. System and method for processing an image sequence
US6845261B2 (en) * 2001-04-23 2005-01-18 General Electric Company System for correlation of MR images with physiological data
US7039455B1 (en) 2001-10-09 2006-05-02 Medrad, Inc. Apparatus and method for removing magnetic resonance imaging-induced noise from ECG signals
AU2003274980A1 (en) * 2002-09-17 2004-04-08 U.S. Government As Represented By The Secretary Of The Army Needle with fiberoptic capability
WO2004080302A1 (ja) * 2003-03-14 2004-09-23 Hitachi Medical Corporation 磁気共鳴イメージング装置
US20050059880A1 (en) * 2003-09-11 2005-03-17 Mathias Sanjay George ECG driven image reconstruction for cardiac imaging
WO2005103750A1 (en) 2004-04-27 2005-11-03 Koninklijke Philips Electronics N.V. Magnetic resonance imaging
US7738943B2 (en) * 2004-11-19 2010-06-15 General Electric Company System and method for identifying periodic motion complexes for MR imaging triggering
US20060173249A1 (en) * 2005-01-11 2006-08-03 Invivo Corporation Method of producing a synthesized bipolar ECG waveform from a scalar ECG waveform
DE102005027438B4 (de) * 2005-06-14 2011-12-22 Siemens Ag Verfahren zur EKG-Triggerung einer Messsequenz eines Magnetresonanzgeräts
US8626266B1 (en) 2006-06-01 2014-01-07 Perinatronics Medical Systems, Inc. ECG triggered heart and arterial magnetic resonance imaging
DE102006053613A1 (de) * 2006-11-14 2008-05-29 Siemens Ag Vorrichtung und Verfahren zum Synchronisieren einer Magnetresonanzanlage mit dem Herzrhytmus eines Patienten
DE102006058332A1 (de) * 2006-12-11 2008-06-19 Siemens Ag Verfahren zur Steuerung einer Messung eines Magnetresonanzgeräts anhand eines Elektrokardiogramm-Signals
US8554311B2 (en) 2011-06-17 2013-10-08 General Electric Company System and method of noise reduction in an electrocardiology study
WO2013048667A1 (en) * 2011-09-29 2013-04-04 General Electric Company System and method of noise reduction in an electrocardiology study
US8659297B2 (en) 2012-02-27 2014-02-25 Perinatronics Medical Systems, Inc. Reducing noise in magnetic resonance imaging using conductive loops
JP2016528982A (ja) * 2013-08-07 2016-09-23 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 磁気共鳴イメージングのための改善されたecgに基づくトリガ
US10751001B2 (en) 2013-12-03 2020-08-25 General Electric Company Systems and methods for tracking and analysis of electrical-physiological interference
WO2015163369A1 (ja) 2014-04-25 2015-10-29 株式会社東芝 心電波形検出装置、及び撮像装置
DE102016220473A1 (de) * 2016-10-19 2018-04-19 Siemens Healthcare Gmbh Detektieren und Unterdrücken von Störsignalen eines differenziellen Spannungsmesssystems
JP6943634B2 (ja) * 2017-06-02 2021-10-06 キヤノンメディカルシステムズ株式会社 信号処理装置、撮像装置及び信号処理方法
US11457834B2 (en) 2018-07-20 2022-10-04 General Electric Company System and method for generating ECG reference data for MR imaging triggering

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4136690A (en) * 1977-10-31 1979-01-30 Del Mar Avionics Method and apparatus for vector analysis of ECG arrhythmias
IL67815A (en) * 1982-02-12 1988-01-31 Sanz Ernst Method and apparatus for cardiogonometry
CN1009330B (zh) * 1987-03-30 1990-08-29 创建基立有限公司 一种电脑电信号检测处理装置
US4850370A (en) * 1987-07-22 1989-07-25 Dower Gordon E Method and apparatus for sensing and analyzing electrical activity of the human heart
US5000182A (en) * 1989-08-11 1991-03-19 Picker International, Inc. Cardiac synchronization magnetic resonance imaging
EP0443120A1 (de) * 1990-01-22 1991-08-28 Studer Revox Ag Verfahren zur Ermittlung von Messwerten am menschlichen und tierischen Körper
JPH03228745A (ja) * 1990-02-05 1991-10-09 Toshiba Corp Mri装置
US5477144A (en) * 1993-02-18 1995-12-19 Johns Hopkins University Method for reducing image artifacts in magnetic resonance images acquired with synchronization to physiological cycles
JP3496889B2 (ja) * 1993-06-24 2004-02-16 株式会社東芝 磁気共鳴イメージング装置
US5458116A (en) * 1994-03-30 1995-10-17 Siemens Medical Systems, Inc. Display system and method for three-dimensional vectorcardiograms
US5692508A (en) * 1996-04-15 1997-12-02 Siemens Medical Systems, Inc. Cardiac-gated 3-dimensional MR angiography
AU2414497A (en) * 1996-04-25 1997-11-12 Per Samuel Bjorgaas Method and instrument for examination of heart/arteries
US5830143A (en) * 1997-01-21 1998-11-03 Wisconsin Alumnin Research Foundation Gated time-resolved contrast-enhanced 3D MR angiography

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011524770A (ja) * 2008-06-18 2011-09-08 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁気共鳴イメージングのための心電計及びそのための電極パッチ
JP2013000590A (ja) * 2011-06-17 2013-01-07 General Electric Co <Ge> 心電図検査におけるノイズ検出システム及び方法
JP2016093488A (ja) * 2014-11-06 2016-05-26 株式会社東芝 信号処理装置、撮像装置及び信号処理方法
US11058361B2 (en) 2014-11-06 2021-07-13 Canon Medical Systems Corporation Signal processing apparatus, imaging apparatus, and signal processing method

Also Published As

Publication number Publication date
DE69828372D1 (de) 2005-02-03
JP4104674B2 (ja) 2008-06-18
DE69828372T2 (de) 2005-12-08
EP0925019A1 (en) 1999-06-30
WO1999004688A1 (en) 1999-02-04
EP0925019B1 (en) 2004-12-29
US5987348A (en) 1999-11-16

Similar Documents

Publication Publication Date Title
JP4104674B2 (ja) Ecgによってトリガされるmr撮像方法及び装置
US11026636B2 (en) Method for generating a medical data set of a moving body part
Fischer et al. Novel real‐time R‐wave detection algorithm based on the vectorcardiogram for accurate gated magnetic resonance acquisitions
KR890001910B1 (ko) 동맥 흐름에 대한 비침투성 mr 영상 형성 향상장치
US6597937B2 (en) Self-adaptive tracking and phase encoding during data collection for contrast-enhanced MRA and dynamic agent uptake studies
EP3413075B1 (en) Mri involving the generation of a physiological motion signal using a pilot tone navigator
US5603322A (en) Time course MRI imaging of brain functions
US5303705A (en) Evoked 23NA MR imaging of sodium currents in the brain
US5000182A (en) Cardiac synchronization magnetic resonance imaging
KR20040071643A (ko) 신체 동작에 자기 공명 영상화 데이터를 동기화시키는 방법
JPH0616771B2 (ja) 磁気共鳴イメージングのための生体電位信号プロセッサ
US5653233A (en) Method and apparatus for improved temporal resolution in dynamic MRI
JPH10512482A (ja) 画像アーティファクトを抑制したデジタル消去磁気共鳴血管造影法
KR20140031820A (ko) 생체 검사 대상자를 자극하는 것에 의한 검사 대상자의 사전 결정된 체적 세그먼트의 mr 영상화 방법 및 자기 공명 시스템
JP4152630B2 (ja) 可変サンプリングを用いた高速/息止め3dmrデータ獲得方法と装置
US20120281901A1 (en) Magnetic resonance imaging apparatus and blood vessel image capturing method
JP2004508859A (ja) 磁気共鳴方法
JP5037866B2 (ja) 磁気共鳴イメージング装置
US4895157A (en) Magnetic resonance imaging system
US20140303482A1 (en) Magnetic resonance imaging method for imaging components with short transverse relaxation times (t2) in a human or an animal heart
WO2005022185A1 (en) Method for cardiac magnetic resonance imaging
Tsukada et al. Noninvasive visualization of multiple simultaneously activated regions on torso magnetocardiographic maps during ventricular depolarization
EP4100756B1 (en) Method to shim a magnetic field and magnetic resonance measurement arrangement
EP2491858A1 (en) Magnetic resonance tomography (MRT) apparatus and method of operating a magnetic resonance (MR) apparatus
JPH0622926A (ja) 磁気共鳴映像装置及び磁気共鳴映像方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050707

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070306

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20070605

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20070723

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070906

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20080226

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080326

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110404

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120404

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120404

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130404

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140404

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term