知財求人 - 知財ポータルサイト「IP Force」
▶ バイオセンス・ウエブスター・（イスラエル）・リミテッドの特許一覧
特開2015-147047ＥＣＧシグナル中に現れるＭＲＩシークエンシングノイズの動的解除
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】特開2015-147047(P2015-147047A)
(43)【公開日】2015年8月20日
(54)【発明の名称】ＥＣＧシグナル中に現れるＭＲＩシークエンシングノイズの動的解除
(51)【国際特許分類】
A61B 5/0402 20060101AFI20150724BHJP
A61B 5/055 20060101ALI20150724BHJP
【ＦＩ】
A61B5/04 310M
A61B5/05 390
【審査請求】未請求
【請求項の数】30
【出願形態】ＯＬ
【外国語出願】
【全頁数】23
(21)【出願番号】特願2015-21917(P2015-21917)
(22)【出願日】2015年2月6日
(31)【優先権主張番号】14/174,951
(32)【優先日】2014年2月7日
(33)【優先権主張国】US
(71)【出願人】
【識別番号】511099630
【氏名又は名称】バイオセンス・ウエブスター・（イスラエル）・リミテッド
【氏名又は名称原語表記】Ｂｉｏｓｅｎｓｅ Ｗｅｂｓｔｅｒ （Ｉｓｒａｅｌ）， Ｌｔｄ．
(74)【代理人】
【識別番号】100088605
【弁理士】
【氏名又は名称】加藤 公延
(74)【代理人】
【識別番号】100130384
【弁理士】
【氏名又は名称】大島 孝文
(72)【発明者】
【氏名】アサフ・ゴバリ
(72)【発明者】
【氏名】ヤロン・エフラス
【テーマコード（参考）】
4C027
4C096
【Ｆターム（参考）】
4C027AA02
4C027BB05
4C027CC01
4C027KK03
4C096AB41
4C096AD19
4C096FC20
(57)【要約】 （修正有）
【課題】磁気共鳴装置において好適なＥＣＧシグナルを検出する方法を提供する。【解決手段】初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて患者を撮像すると同時に、患者からの心電図（ＥＣＧ）シグナルの初期セットを受信することからなる、方法。初期ＭＲＩシーケンスから発生しているＥＣＧシグナルの初期セット中の初期ノイズを同定することと、その周波数成分に関する初期のノイズを特性化することとを含む。更に、初期ＭＲＩシーケンスと、周波数成分とを定義しているパラメータ間の関係を生成することと、該関係を、続くＭＲＩシーケンスのパラメータに適用することによって、フィルタを計算することと、続くＭＲＩシーケンスを用いて患者を撮像すると同時に、患者から受信した、更なる組のＥＣＧシグナル中のノイズを減少させるために、計算したフィルタを適用することと、が含まれる。
【選択図】図１
【特許請求の範囲】
【請求項1】
初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて患者を撮像すると同時に、前記患者から心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、計算された前記フィルタを適用して、前記患者から受信した、ＥＣＧシグナルの更なるセットにおけるノイズを減少させることと、を含む方法。
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、計算された前記フィルタを適用して、前記患者から受信した、ＥＣＧシグナルの更なるセットにおけるノイズを減少させることと、を含む方法。
【請求項2】
前記初期ノイズを同定することが、前記患者を撮像する前にＥＣＧシグナルの参照セットを生成することと、前記ＥＣＧシグナルの初期セットを前記参照セットと比較することとを含む、請求項１に記載の方法。
【請求項3】
前記初期ノイズを特性化することが、前記初期ノイズ上で、フーリエ変換を実施することを含む、請求項１に記載の方法。
【請求項4】
前記初期ＭＲＩシーケンスを定義している前記パラメータが、前記初期ＭＲＩシーケンスの電圧及び磁場パルスを記述している変数を含む、請求項１に記載の方法。
【請求項5】
前記関係を生成することが、前記初期ＭＲＩシーケンスを使用すると同時に、前記変数の値を変更することと、前記値における変化を、前記周波数成分の振幅における変化と比較することと、を含む、請求項４に記載の方法。
【請求項6】
前記初期ＭＲＩシーケンスと、前記続くＭＲＩシーケンスが、共通パラメータを有する共通ＭＲＩシーケンスを含む、請求項１に記載の方法。
【請求項7】
前記共通パラメータが、前記共通ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含み、該変数が、前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスに対して異なる値を有する、請求項６に記載の方法。
【請求項8】
前記フィルタを計算することが、前記周波数成分の予想される振幅に関して、予想されるノイズを計算することを含み、計算された前記フィルタを適用することが、前記ＥＣＧシグナルの更なるセットと関連した、対応する周波数の振幅から、前記予想される振幅を差し引くことを含む、請求項１に記載の方法。
【請求項9】
前記ＥＣＧシグナルが、前記患者内に配置したプローブから獲得した、内部ＥＣＧシグナルを含む、請求項１に記載の方法。
【請求項10】
前記ＥＣＧシグナルが、前記患者の皮膚から獲得した外部ＥＣＧシグナルを含む、請求項１に記載の方法。
【請求項11】
前記初期及び続くＭＲＩシーケンスが、共通の複数の相を含み、前記関係を生成することが、前記初期ＭＲＩシーケンスの選択された相を定義しているパラメータと、前記周波数成分との間の相関係を生成することを含み、前記フィルタを計算することが、前記続くＭＲＩ配列の選択された相のパラメータに、前記相関係を適用することによって、前記フィルタを計算することを含む、請求項１に記載の方法。
【請求項12】
前記ＥＣＧシグナルの初期セットが、前記患者に連結した多数のＥＣＧリード上に発生し、前記初期ノイズを同定することが、前記多数のリード上のシグナルを合計することを含む、請求項１に記載の方法。
【請求項13】
前記初期ノイズを同定することが、前記患者を撮像する前に獲得された参照ＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、請求項１に記載の方法。
【請求項14】
前記初期ノイズを同定することが、前記患者の撮像の間に、前記初期ノイズを同定することを含む、請求項１に記載の方法。
【請求項15】
前記初期ノイズを同定することが、前記患者に接続した固定電極からのＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、請求項１に記載の方法。
【請求項16】
患者からの心電図（ＥＣＧ）シグナルの初期セットを獲得するために設定された、１つ又は複数の電極と、
プロセッサであって、
初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて前記患者を撮像すると同時に、前記患者からの前記心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、前記患者から受信した、ＥＣＧシグナルの更なるセット中のノイズを減少させるために、計算された前記フィルタを適用することと、のために設定された、プロセッサと、を含む装置。
プロセッサであって、
初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて前記患者を撮像すると同時に、前記患者からの前記心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、前記患者から受信した、ＥＣＧシグナルの更なるセット中のノイズを減少させるために、計算された前記フィルタを適用することと、のために設定された、プロセッサと、を含む装置。
【請求項17】
前記初期ノイズを同定することが、前記患者を撮像することの前に、ＥＣＧシグナルの参照セットを生成することと、前記ＥＣＧシグナルの初期セットと、前記参照セットを比較することと、を含む、請求項１６に記載の装置。
【請求項18】
前記初期ノイズを特性化することが、前記初期ノイズの上でフーリエ変換を実施することを含む、請求項１６に記載の装置。
【請求項19】
前記初期ＭＲＩシーケンスを定義している前記パラメータが、前記初期ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含む、請求項１６に記載の装置。
【請求項20】
前記関係を生成することが、前記初期ＭＲＩシーケンスを使用すると同時に前記変数の値を変更することと、前記値における変化と、前記周波数成分の振幅における変化とを比較することと、を含む、請求項１９に記載の装置。
【請求項21】
前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスが、共通パラメータを有する共通ＭＲＩシーケンスを含む、請求項１６に記載の装置。
【請求項22】
前記共通パラメータが、前記共通ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含み、該変数が、前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスに対する異なる値を有する、請求項２１に記載の装置。
【請求項23】
前記フィルタを計算することが、前記周波数成分の予想される振幅に関して、予想されるノイズを計算することを含み、計算された前記フィルタを適用することが、前記ＥＣＧシグナルの更なるセットと関連した対応する周波数の振幅から前記予想される振幅を差し引くことを含む、請求項１６に記載の装置。
【請求項24】
前記１つ又は複数の電極が前記患者内に位置し、前記ＥＣＧシグナルに、前記１つ又は複数の電極から獲得された初期ＥＣＧシグナルが含まれる、請求項１６に記載の装置。
【請求項25】
前記１つ又は複数の電極が、前記患者の皮膚上に位置し、前記ＥＣＧシグナルが、前記１つ又は複数の電極から獲得された外部ＥＣＧシグナルを含む、請求項１６に記載の装置。
【請求項26】
前記初期及び続くＭＲＩシーケンスが、共通の複数の相を含み、前記関係を生成することが、前記初期ＭＲＩシーケンスの選択された相を定義しているパラメータと、前記周波数成分との間の相関係を生成することを含み、前記フィルタを計算することが、前記続くＭＲＩシーケンスの選択された相のパラメータに対して、前記相関係を適用することによって前記フィルタを計算することを含む、請求項１６に記載の装置。
【請求項27】
前記ＥＣＧシグナルの初期セットを獲得するために設定された多数のＥＣＧリードを含み、前記初期ノイズを同定することが、前記多数のＥＣＧリード上のシグナルを合計することを含む、請求項１６に記載の装置。
【請求項28】
前記初期ノイズを同定することが、前記患者を撮像する前に獲得した参照ＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、請求項１６に記載の装置。
【請求項29】
前記初期ノイズを同定することが、前記患者の撮像の間に、前記初期ノイズを同定することを含む、請求項１６に記載の装置。
【請求項30】
前記１つ又は複数の電極が、前記患者に接続した固定電極を含み、前記初期ノイズを同定することが、前記固定電極からのＥＣＧシグナルを用いて前記初期ノイズを同定することを含む、請求項１６に記載の装置。
【発明の詳細な説明】
【技術分野】
【0001】
（関連出願の相互参照）本明細書は、本発明の付与者に付与され、参考文献によって本明細書に組み込まれている、同日付にて申請された、「Ａｎａｌｏｇ Ｃａｎｃｅｌｌａｔｉｏｎ ｏｆ ＭＲＩ Ｓｅｑｕｅｎｃｉｎｇ Ｎｏｉｓｅ Ａｐｐｅａｒｉｎｇ ｉｎ ａｎ ＥＣＧ Ｓｉｇｎａｌ」という表題の米国特許明細書に関連する。
【0002】
（発明の分野）本発明は一般に、心電図（ＥＣＧ）シグナルに関し、とりわけ磁気共鳴映像法（ＭＲＩ）手順の間、ＥＣＧシグナルを検出することに関する。
【背景技術】
【0003】
磁気共鳴映像法（ＭＲＩ）は、患者の組織、特に軟組織を可視化するための、極めて強力な技術である。この技術は、原子核、典型的には水素原子核をそれらの平衡状態から励起することと、原子核が平衡状態に緩和するときに原子核により放出される共鳴高周波信号を測定することとに依存する。
【0004】
ＭＲＩ手順において、該手順によって発生した、急速に変化する磁場中に配置された患者の心状態が査定される必要がありうる。心電図（ＥＣＧ）シグナルのモニタリングが、患者の心状態の良好な指標であり、ＥＣＧモニタリングが該手順の効率を増強しうる。
【0005】
その開示内容が参考により本明細書に組み込まれている、Ｂｒｏｓｏｖｉｃｈらに付与された米国特許第７，０３９，４５５号明細書は、ＭＲＩを受けている患者から得たＥＣＧシグナルの品質を改善するための装置を記述している。装置には、差動増幅計、プレフィルタ、シグナルリミッター回路、及び積分ローパスフィルタを有する中間増幅計が含まれる。
【0006】
その開示内容が参考により本明細書に組み込まれている、Ｃｏｈｅｎに付与された米国特許第７，２８６，８７１号明細書は、繰り返し干渉汚染の存在下記録された、電気シグナルの汚染を減少させる方法を記述している。電気シグナルはデジタル化され、デジタル化は、タイミングシグナルと共に開始する。複数のデジタル化した電気シグナルを解析し、該電気シグナルをタイミングシグナルに対して同期化して、デジタル化した電気シグナルから差し引いた推定汚染シグナルを得る。
【0007】
その開示内容が参考により本明細書に組み込まれている、Ｍｏｏｒｅに付与された、米国特許第４，９９１，５８０号明細書は、ＭＲＩを受けている患者から得たＥＣＧシグナルの品質を改善するための方法を開示している。本方法には、ＭＲＩ誘導ノイズシグナルを有するＥＣＧシグナルを、先に選択された最大立ち上がり速度を有する、立ち上がり速度リミッター（ＳＲＬ）回路の入力に対して伝導することが含まれる。ＳＲＬ回路の出力は、ローパスフィルタ回路に対して連結される。
【0008】
その開示内容が参考により本明細書に組み込まれている、Ｋｒｅｇｅｒらに付与された米国特許第６，０７０，０９７号明細書は、心ＭＲＩに対するゲーティングシグナルを生成するための方法を記述している。ＭＲＩシステムには、スキャンされている患者からのＥＣＧシグナルを受信し、ゲーティングシグナルを産出する検出器システムが含まれる。ゲーティングシグナルは、ＥＣＧシグナル内の検出されたピークが、一組のＲ−波判定基準に見合った時に産出される。
【0009】
その開示内容が参考により本明細書に組み込まれている、Ｕｕｔｅｌａ Ｋｉｍｍｏに付与された欧州特許第１，８７２，７１５号明細書は、乱流電磁環境によって生成され、回収される、バイオシグナルアーチファクトの統計学的特徴の参照データ指標を記述している。患者のバイオシグナルの多数のチャネルを、乱流電磁環境中で測定する。多数のチャンネル中のアーチファクトを、参照データを用いて検出し、パラメータを、多数のチャネルの線形組み合わせに対して誘導する。正確なバイオシグナルを、多数のチャネルの望むシグナル試料に、線形組み合わせを適用することによって得、線形組み合わせは、誘導されたパラメータによって定義される。
【0010】
その開示内容が参考により本明細書に組み込まれている、Ｓｃｈｗｅｉｔｚｅｒらに付与された国際特許第ＷＯ／２０１２／１７０１１９号は、ＭＲＩスキャンシーケンスの間に、患者の体で、カテーテル電極位置を追跡するためのシステムが記述されている。システムには、潜在的にノイズを誘導する条件の間にとられ、したがってノイズによる改変にさらされた、１つ又は複数のインピーダンス測定を同定するために設定された緩和ロジックが含まれる。該緩和ロジックを、潜在的に間違いの多いインピーダンス測定を、直前に進められている獲得サイクルからとった、先に得たインピーダンス測定に置換するように構成する。
【0011】
参考により本特許出願に組み込まれた文書は、これらの組み込まれた文書内のどんな用語でも、本明細書で明示的又は暗黙的に行われる定義と相反するように定義される場合を除き、本出願の一体部分と見なされるべきであり、本明細書における定義のみが検討されるべきである。
【発明の概要】
【発明が解決しようとする課題】
【0012】
本発明の一実施形態は、次のことを含む方法を提供する：初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて患者を撮像すると同時に、前記患者から心電図（ＥＣＧ）のシグナル初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、計算された前記フィルタを適用して、前記患者から受信した、ＥＣＧシグナルの更なるセットにおけるノイズを減少させること。
【課題を解決するための手段】
【0013】
開示された実施形態において、前記初期ノイズを同定することが、前記患者を撮像する前にＥＣＧシグナルの参照セットを生成することと、前記ＥＣＧシグナルの初期セットを前記参照セットと比較することとを含む。
【0014】
さらに開示された実施形態において、前記初期ノイズを特性化することが、前記初期ノイズ上で、フーリエ変換を実施することを含む。
【0015】
また更なる開示された実施形態において、前記初期ＭＲＩシーケンスを定義しているパラメータが、初期ＭＲＩシーケンスの電圧及び磁場パルスを記述している変数を含む。典型的に、前記関係を生成することが、初期ＭＲＩシーケンスを使用すると同時に、変数の値を変更することと、値における変化を、周波数成分の振幅における変化と比較することと、を含む。
【0016】
他の実施形態において、前記初期ＭＲＩシーケンスと、前記続くＭＲＩシーケンスが、共通のパラメータを有する共通ＭＲＩシーケンスを含む。典型的に、前記共通パラメータが、前記共通ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含み、該変数が、前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスに対して異なる値を有する。
【0017】
更なる他の実施形態において、フィルタを計算することが、前記周波数成分の予想される振幅に関して、予想されるノイズを計算することを含み、計算された前記フィルタを適用することが、前記ＥＣＧシグナルの更なるセットと関連した、対応する周波数の振幅から、前記予想される振幅を差し引くことを含む。
【0018】
ＥＣＧシグナルは、患者内に配置したプローブから獲得した、内部ＥＣＧシグナルからなってよい。あるいは、又は加えて、ＥＣＧシグナルは、患者の皮膚から獲得した外部ＥＣＧシグナルからなってよい。
【0019】
前記初期及び続くＭＲＩシーケンスが、共通の複数の相を含み、関係を生成することが、初期ＭＲＩシーケンスの選択された相を定義しているパラメータと、周波数成分との間の相関係を生成することからなってよく、フィルタを計算することが、続くＭＲＩ配列の選択された相のパラメータに、前記相関係を適用することによって、フィルタを計算することからなってよい。
【0020】
またさらに他の実施形態において、初期セットのＥＣＧシグナルが、患者に連結した多数のＥＣＧリード上に発生し、初期ノイズを同定することが、多数のリード上のシグナルを合計することを含む。
【0021】
典型的に、初期ノイズを同定することが、患者を撮像する前に獲得された参照ＥＣＧシグナルを用いて、初期ノイズを同定することを含む。
【0022】
あるいは、又は加えて、初期ノイズを同定することが、患者の撮像の間、初期のノイズを同定することを含む。
【0023】
さらにあるいは、又は加えて、初期ノイズを同定することが、患者に接続した固定電極からのＥＣＧシグナルを用いて、初期ノイズを同定することを含む。
【0024】
本発明の一実施形態により更に、患者からの心電図（ＥＣＧ）シグナルの初期セットを獲得するために設定された、１つ又は複数の電極と、
プロセッサであって、
初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて前記患者を撮像すると同時に、前記患者からの前記心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義するパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、前記患者から受信した、ＥＣＧシグナルの更なるセット中のノイズを減少させるために、計算された前記フィルタを適用することと、のために設定されたプロセッサと、を含む装置が提供される。
【0025】
本発明は、以下の詳細な実施形態の説明を、図面と併せ読むことによって、より十分に理解されるであろう。
【図面の簡単な説明】
【0026】
【図1】本発明の一実施形態にしたがった、心電図（ＥＣＧ）シグナル中に現れている磁気共鳴画像法（ＭＲＩ）ノイズの動的解除のための略図的絵図。
【図2】本発明の一実施形態にしたがった、ＭＲＩ手順の間に発生したＭＲＩシグナルのシーケンスを略図的に示している、一組の電圧対時間、及び磁場対時間グラフ。
【図3】本発明の一実施形態にしたがった、ＭＲＩ手順の間に発生したＭＲＩシグナルのシーケンスを略図的に示している、他の組の電圧対時間、及び磁場対時間グラフ。
【発明を実施するための形態】
【0027】
概論本発明の実施形態は、磁気共鳴映像法（ＭＲＩ）手順が患者において実施されていると同時に操作する、ノイズキャンセリングシステムを提供する。ＭＲＩ手順の間、心内及び／又は皮膚電極を用いて獲得された心電図（ＥＣＧ）シグナルが患者より受信される。ＭＲＩ手順は、獲得したＥＣＧシグナルにノイズを注入し、システムが本ノイズを減少又は除去する。
【0028】
典型的に、システムのトレーニングセクションにて、ＭＲＩ手順が実行されない一方で、プロセッサが、発生するＥＣＧシグナルの参照セットを保存する。ＥＣＧシグナルを次いで、選択されたＭＲＩシーケンスを実施していると同時に獲得し、プロセッサが、参照セットを使用して、注入されたノイズを同定する。プロセッサが、その周波数成分に関して、ノイズを特性化し、周波数成分と、選択されたＭＲＩシーケンスを定義しているパラメータ、及びパラメータと関連した変数との間の関係を生成する。典型的には、ノイズが特性化され、ＭＲＩシーケンスそれぞれを順に操作し、各異なるシーケンスに対してノイズを同定することによって、対応する関係が生成される。
【0029】
システムの操作セクションにおいて、プロセッサが、手順の間に使用されている所定のＭＲＩシーケンスに関する関係を選択し、該関係から、所定のＭＲＩシーケンスが患者に対して適用されている時に、使用に関してノイズフィルタを計算する。フィルタを、受信したＥＣＧシグナル中のノイズを減少、又は除去さえするように、所定のＭＲＩシーケンスの適用の間に受信したＥＣＧシグナルに適用する。
【0030】
周波数成分に関してノイズを特性化することと、これらを特定のＭＲＩシーケンスのパラメータに関連づけることによって、本発明の実施形態が、ＭＲＩ手順の間に使用した各異なるＭＲＩシークエンに対して、対応する異なるフィルタを生成可能である。
【0031】
本発明のいくつかの実施形態が、完全なＭＲＩシーケンスのためのノイズフィルタを生成して良い一方で、他の実施形態が、所定のＭＲＩシーケンスの別の個々の相に対して、対応する、異なるノイズフィルタを生成してよい。
【0032】
詳細な説明ここで、本発明の実施形態にしたがって、心電図（ＥＣＧ）シグナル中に現れている磁気共鳴映像法（ＭＲＩ）シークエンシングノイズの動的解除のための、システム２０の略図的絵図である、図１を参照する。システム２０は、ＭＲＩスキャナ２２と、カテーテル等のプローブ２４と、制御コンソール２６とを備える。プローブ２４は、プローブの遠位末端３４における１つ又は複数の電極３２を用いて、患者３０の心腔２８のチャンバーにおけるＥＣＧシグナルを獲得するために使用されてよい。電極３２によって獲得されたシグナルを本明細書では、内部ＥＣＧシグナルと呼ぶ。いくつかの実施形態では、プローブ２４は、心臓焼灼を行うなど、追加の目的のために使用され得る。更に別の方法としては、プローブ２４は、変更すべきところは変更して、心臓又はその他の身体器官において、その他の治療及び／又は診断目的のために使用され得る。電極３２によって獲得された内部ＥＣＧシグナルを、典型的には、プローブ２４中のコンダクタ及び／又は光ファイバーによって、制御コンソール２６まで転送し、そこでシグナルが解析されてよい。
【0033】
プローブ２４を用いて、内部ＥＣＧシグナルを獲得することに加えて、システム２０は典型的にまた、典型的に患者の皮膚上に多数の導電性パッチ３６を配置することによって、患者３０の皮膚からＥＣＧシグナルを獲得する。パッチ３６によって獲得されたシグナルを本明細書では、外部ＥＣＧシグナルと呼ぶ。外部ＥＣＧシグナルを、ケーブル３８を介して、シグナルを解析する制御コンソール２６に伝達する。解析から得た結果を、システム２０の操作者４２に対して、ディスプレイ４０上に提示してよい。
【0034】
操作者４２、典型的には心臓専門医が、プローブの遠位末端３４が、本明細書でそこより内部ＥＣＧシグナルが獲得される心腔であると仮定される、体腔に入るように、患者３０の血管系を通してプローブ２４を挿入する。典型的に、プローブの遠位末端を、本技術分野で公知の方法によって追跡する。磁気位置追跡技術は、例えば、米国特許第５，３９１，１９９号明細書、同５，４４３，４８９号明細書、第６，７８８，９６７号明細書、同第６，６９０，９６３号明細書、同第５，５５８，０９１号明細書、同第６，１７２，４９９号明細書、同第６，１７７，７９２号明細書に、記載されており、それらの開示内容を参照により本願に援用する。インピーダンスに基づく位置追跡技術は、例えば米国特許第５，９８３，１２６号明細書、同第６，４５６，８６４号明細書、及び同第５，９４４，０２２号明細書に説明され、その開示内容が同様に参照により本明細書に組み込まれる。
【0035】
ＭＲＩスキャナ２２は、一緒になって空間的に異なる磁場Ｂ（ｘ、ｙ、ｚ）を生成する、傾斜磁場コイルなどの磁場コイル５０を備える。空間的に異なる磁場は、スキャナ内で生成される高周波（ＲＦ）信号の空間位置確認を提供する。更に、スキャナは、送信／受信コイル５２を備える。送信モードにおいて、コイル５２は患者３０へＲＦパルス化エネルギーを放射し、エネルギーのＲＦパルスが、患者の組織の原子核スピンと相互作用し、それによって、それらの平衡位置から離れて、原子核の磁気モーメントを再調整する。受信モードにおいて、コイル５２は、組織の原子核がそれらの平衡状態に緩和するときに患者の組織から受信するＲシグナル号を検出する。所定の領域内の原子核の緩和によって生成されるシグナルの周波数（ラーモア周波数）は、その領域の磁場と正比例の関係であり、比例定数は原子の磁気回転比γで与えられる。したがって、水素原子核に関して下記式（１）：【数1】
Ｂ（ｘ、ｙ、ｚ）は、この点の磁場であり、
【数2】
【0036】
いくつかの実施形態において、スキャナ２２を、スキャナプロセッサ５６によって操作し、以上で記述したＥＣＧシグナルと、内部ＥＣＧシグナルを獲得しているプローブの追跡を、ＥＣＧプロセッサ５８によって解析及び制御する。単純化のために、本明細書の記述において、コンソール２６内の単一プロセッサ６０が、システム２０を操作すると推定され、当業者が、１つ又は複数のプロセッサがシステムを操作する事象における記述を適合可能であろう。
【0037】
したがって、以上を参照してＥＣＧシグナルを解析することと、内部ＥＣＧシグナルを生成しているプローブを追跡することとに加えて、プロセッサ６０は、必要とされた磁場勾配、並びに送信／受信コイル５２を操作するための他の回路を形成することを含む、制御コイル５０に対する回路を用いることによって、スキャナ２２を操作する。
【0038】
典型的には、プロセッサ６０は一般用途コンピュータで構成され、コンピュータには、本明細書に記載する機能を実行するソフトウェアがプログラムされている。ソフトウェアは、例えば、ネットワークを介して電子形式でプロセッサ６０にダウンロードされてもよく、或いは、光学的、磁気的、又は電子的記憶媒体などの有形の持続的な媒体上で提供されてもよい。別の方法としては、プロセッサ６０のいくつかの又は全ての機能は、専用の、又はプログラム可能なデジタルハードウェア構成要素によって、又はハードウェア要素とソフトウェア要素との組み合わせを使用することによって、実施されてもよい。
【0039】
図２は、本発明の一実施形態にしたがった、ＭＲＩ手順の間に生成されるＭＲＩシグナルのシーケンス７０を略図的に説明している、第一組の電圧（Ｖ）対時間（ｔ）及び磁場（Ｂ）対時間（ｔ）グラフである。第一の電圧対時間グラフ７２は、シーケンス７０の開始時に送信される、コイル５２によって生成された、送信ＲＦパルス７４を説明している。送信ＲＦパルス長は典型的に、より大きく、又はより小さくてもよいが、２ｍｓのオーダーである。第一磁場対時間グラフ７６中に説明された、送信ＲＦパルスを包含することは、コイル５０によって生成されたスライス選択勾配（Ｇｓｓ）磁場パルスである。スライス選択勾配場は、スキャナ２２によって撮像されるべきである、患者３０中の所望の容量を同定する。
【0040】
第二の磁場対時間グラフ８０は、所望の容量内の点の垂直位を選択する、相コード勾配（Ｇｐｅ）場パルスを説明している。第三の磁場対時間グラフ８４は、所望の容量内の点の水平位を選択する、周波数コード勾配（Ｇｆｅ）場パルスを説明している。
【0041】
第二の電圧対時間グラフ８８は、送信ＲＦパルスに対する応答において、コイル５２によって受信したデータ獲得シグナルに対応している、受信ＲＦパルス９０を説明している。
【0042】
ＲＦ送信パルスの中心と、ＲＦ受信パルスの中心間の典型的な時間は、およそ３０ｍｓであってよい。全シーケンスに対する典型的な時間は、およそ４０ｍｓであってよい。しかしながら、実際の時間は、これらの値よりも長く、又は短くてもよい。
【0043】
ＭＲＩ手順の経過の間、シーケンス７０によって説明された、ＲＦパルスのシーケンスと磁場設定は、典型的には１ｓのオーダーの反復速度で、繰り返してよい。シーケンス７０の各反復において、磁場パルスを定義している１つ又は複数の変数は、典型的に患者３６の異なる領域がスキャンされうるように、様々であってよい。いくつかの場合、シーケンス７０を繰り返すにつれ、その振幅、周波数又は相のような、ＲＦ送信パルス７４の変数が変化してよい。
【0044】
シーケンス７０は、ＭＲＩ手順にて典型的に使用される、勾配エコーシーケンスを意味する、１つの型のＭＲＩシーケンスを説明する。
【0045】
図３は、本発明の一実施形態にしたがった、ＭＲＩ手順の間に生成されたＭＲＩシグナルのシーケンス１００を略図的に説明している、第二組の電圧対時間と磁場対時間グラフである。シーケンス７０に対してのように、シーケンス１００は、送信及び受信電圧対時間グラフ１１０及び１１２と、３つの磁場対時間グラフ１２０、１２２及び１２４とを含む。グラフ中で説明するように、シーケンス１００は、（シーケンス７０の１つのＲＦパルスと１つのＧｆｅパルスと対称的に）２つのＲＦ送信パルスと２つの磁場Ｇｆｅパルスを使用する。シーケンス１００は、ＭＲＩ手順中でまた典型的に使用されてよい、スピンエコーシーケンスを説明する。
【0046】
図２及び３は、本発明の実施形態にて使用されてよい、２つの型のＭＲＩシーケンスを説明している。これらのシーケンス、並びに他の可能性あるシーケンスにおける変化が、当業者に明白であり、すべてのそのようなシーケンスは、本発明の範囲内であると推定される。以下の記述において、固有識別子「ｍ」は、各型のＭＲＩシーケンスを標識するために使用し、例示の方法によって、図２によって説明されたＭＲＩシーケンスが、固有識別子ｍ＝１を有すると推測され、図３によって説明されたＭＲＩシーケンスは、固有識別子ｍ＝２を有すると推測される。
【0047】
図２及び３はまた、各ＭＲＩシーケンスが、シーケンスに対して生成される異なるパルスに対応する、多数の相からなる。ＭＲＩシーケンスの相は、一パルス、多数のパルス、及び／又はパルスの部分に関して定義されてよい。例えば、図２のＭＲＩシーケンスの選択された相は、Ｇｐｅパルスの期間として定義されてよく、図３のＭＲＩシーケンスに対して、他の選択された相は、第二Ｇｆｅパルスの最初の半分の期間として定義されてよい。
【0048】
一般に、任意のＭＲＩシーケンスは、図２及び３でグラフ化された各パルスを記述している変数に関して定義されてよい。したがって、電圧対時間グラフ７２及び１１２を参照して、変数には、シーケンスを開始しているＲＦパルスの周波数、相及び電圧の振幅、時間にわたるこれらの変数の値における変化、任意の続くＲＦパルスの変数、及び初期ＲＦパルスに関する続くＲＦパルスのタイミングが含まれる。磁場対時間グラフ７６、８０及び８４と、グラフ１２０、１２０及び１２４を参照して、変数には、磁場パルスの形状、すなわち時間によって変化するにともなった場の振幅を記述しているもの、ＲＦパルス又はグラフ７２及び１１２におけるパルスに関する、又は互いに関する、磁場パルスのタイミングが含まれる。
【0049】
よりとりわけ、各ＭＲＩシーケンスｍは、一組のパラメータに関して広く特徴付けられてよく、各パラメータは、シーケンスのパルスの一つに対応している。以下の記述において、ＭＲＩシーケンスｍに対するパラメータの組は、｛Ｓ｝ｍとして記述される。
【0050】
例えば、図２にて説明されるＭＲＩシーケンスｍ＝１は、一組のパラメータ｛Ｓ｝１＝｛ＰＲＦ、ＰＧｓｓ、ＰＧｐｅ、ＰＧｆｅ、Ｐシグナル｝を有してよく、本組の各要素が、電圧又は磁場パルスを示している。反対に、図３にて説明されるＭＲＩシーケンスｍ＝２は、一連のパラメータ｛Ｓ｝２＝｛ＰＲＦ１、ＰＲＦ２、ＰＧｓｓ、ＰＧｐｅ、ＰＧｆｅ１、ＰＧｆｅ２、Ｐシグナル｝を有する。
【0051】
任意の所望のＭＲＩシーケンスｍに対して、組中の各パラメータの、さらにより詳細な特徴が存在する。したがって、ｍ＝１の電圧対時間グラフ７２を参照して、パラメータＰＲＦが、変数、シーケンスを開始しているＲＦパルスの電圧の、周波数ｆ、相φ及び振幅Ａ、及び時間ｔにわたるこれらの変数における変化に関して特性化される。ＰＲＦは、したがって、一組の規則的な４−タプルとして表されてよく、下記式（２）【数3】
【0052】
同様に、ｍ＝１の磁場Ｇｓｓ対時間グラフ７６を参照して、ＰＧｓｓが、一組の規則的な組の変数として表されてよいように、パラメータＰＧｓｓは、時間ｔにおける可変磁場Ｂに関して特性化され、以下式（３）【数4】
【0053】
ｍ＝２の電圧対時間グラフ１１０を参照して、パラメータＰＲＦ１及びＰＲＦ２は、式【数5】
ｆ２、φ２、Ａ２は、ｍ＝２の第二パルスの時間ｔにおける、周波数、相及び振幅である）によって表されてよい。
【0054】
他の例として、場対時間グラフ１２４を参照して、パラメータＰＧｆｅ１、ＰＧｆｅ２は、式【数6】
【0055】
スキャナ２２を操作する前に、操作者４２は、各ＭＲＩシーケンスｍに対する（以上で記述される広く、そしてより詳細な様式での）各パラメータに関連したそのパラメータ及び変数に関して、｛Ｓ｝ｍを保存するために、プロセッサ６０を使用する。スキャナ２２が操作される手順の間、プロセッサは、図４のフローチャートを参照して表されるように、｛Ｓ｝ｍのパラメータと変数を想起する。
【0056】
図４は、本発明の一実施形態にしたがった、実行システム２０中で実施された段階のフローチャート２００である。フローチャートは、２つのセクション、トレーニングセクション２０２と操作セクション２０４とに分けられる。フローチャートの記述は、フローチャートの適切な段階が、単一プロセッサ６０を用いて実行されることを想定する。システム２０が、スキャナプロセッサ５６とＥＣＧプロセッサ５８を含む場合、フローチャート段階は、スキャナプロセッサ５６からのＭＲＩシーケンスｍのタイミングと同一性に対してのような、コミュニケーションを使用して、ＥＣＧプロセッサ５８によって実行されると推定されてよい。
【0057】
トレーニング相の初期段階２１０において、操作者４２は、パッチ３６を患者３０にはり、プローブ２４を患者の心臓２８に挿入する。パッチ３６は、外部ＥＣＧシグナルを獲得し、プローブ２４の電極又は電極群が、内部ＥＣＧシグナルを獲得し、一方でスキャナ２２は動作不能である。プロセッサ６０は、獲得したＥＣＧシグナルを受信し、シグナルの参照セットとして、シグナルの選択された群を保存する。
【0058】
典型的に、シグナルの参照セットには、患者３０に関して比較的定常である、電圧からのシグナルが含まれる。例えば、シグナルの参照セットは、パッチ３６からであってよい。いくつかの実施形態において、プローブ２４は、例えば心臓２８中の内部参照プローブとして使用される場合、定常であってよく、この場合、（図１で示していない）他の移動可能なプローブが心臓２８内に挿入されてよい。プローブ２４が定常である場合、ついでその電極又は電極群からのシグナルが、ＥＣＧ参照シグナルとして使用されてよい。いくつかの実施形態において、参照セットには、パッチからのシグナルの選択された１つのような、パッチシグナルのサブセット、及び／又はプローブが定常である場合、プローブ２４の電極からのシグナルのサブセットが含まれてよい。
【0059】
特性化段階２１２において、プロセッサは、操作システム２０にて使用されるべき異なるＭＲＩシーケンスのパラメータ及び変数を保存する。言い換えれば、プロセッサは、手順中で使用されるべきｍの各値に対して、組｛Ｓ｝ｍを保存する。一旦パラメータ及び変数が、保存されたならば、プロセッサは、それによってスキャナを操作するために、ＭＲＩシーケンスｍを選択する。
【0060】
記録段階２１４において、スキャナ２２を、選択されたＭＲＩシーケンスを用いて操作し、プロセッサが、スキャナを操作する一方で生成された、パッチ３６からの、及びプローブ２４の電極からの、ＥＣＧシグナルを獲得及び記録する。記録されたＥＣＧシグナルを、段階２１０中で保存された参照シグナルと一致する、すなわち、参照シグナルとして同一のリード上に存在するように選択する。記録されたＥＣＧシグナルには、スキャナ操作のために、とりわけ、シグナルを伝達しているリード内に注入されたノイズが含まれる。必要性はないけれども、スキャナ操作は典型的に、シーケンスの多数のスキャンに対してであり、操作はシーケンスの単一スキャンに対してであってよい。
【0061】
ノイズ単離段階２１６において、段階２１４の記録されたＥＣＧシグナルを、段階２１０にて保存されたＥＣＧシグナルの参照セットと比較する。比較を用いて、プロセッサ６０が、生成されるノイズシグナルを単離する。プロセッサがついで、典型的には、ノイズシグナルに対してフーリエ変換を適用し、一組の周波数と、関連した周波数振幅を、ノイズシグナルに対して導き出すことによって、周波数にしたがって生成されたノイズを特性化する。
【0062】
ＭＲＩシーケンスｍによって生成されたノイズの周波数と関連した周波数振幅の組を、本明細書で記号｛Ｎ｝ｍによって表示し、【数7】
【0063】
関係形成段階２１８において、プロセッサは、段階２１６にて生成された組｛Ｎ｝ｍと段階２１２にて保存された組｛Ｓ｝ｍ間の関係を生成する。関係は典型的に、解析されている選択されたＭＲＩシーケンスに対して、フローチャート中矢印２２０によって表されるように、｛Ｓ｝ｍの各パラメータの１つ又は複数の変数を変更することにとって生成される。特定の変数の変化が、段階２１６の周波数の組中の、各周波数における変数の効果に関して、プロセッサに情報を提供する。
【0064】
一般化関係は式、【数8】
【0065】
決定段階２２２において、プロセッサは、全てのＭＲＩシーケンスが解析されたかどうかを確認する。シーケンスが残っている場合、フローチャートが段階２２４まで進み、他のＭＲＩシーケンスを選択し、ついで段階２１４に戻る。全てのシーケンスが解析されたならば、トレーニングセクション２０２が完了し、操作セクション２０４が開始される。
【0066】
操作セクションの第一段階２３０において、プロセッサは、使用されているＭＲＩシーケンスｍを同定する。プロセッサはまた、シーケンスのために使用されている｛Ｓ｝ｍの特定の数値を決定する。数値は本明細書で、｛Ｓｐ｝ｍと称され、添字ｐは数値を表している。
【0067】
ノイズ推定段階２３２において、段階２３０の数値を使用して、ＥＣＧシグナル上で推定されるノイズの数値｛Ｎｐ｝ｍを決定する。本決定は、以下の式にしたがって、段階２１８からのマトリックス［Ｍ］ｍを使用する。【数9】
【0068】
式（８）から、｛Ｎｐ｝ｍは、【数10】
【0069】
式（１１）は、ＭＲＩシーケンスｍが動作可能であると同時に獲得されたＥＣＧシグナル上、周波数と周波数振幅の一組の規則的な対として、想定されるノイズを表す。
【0070】
フローチャートの最終フィルトレーション段階２３４において、式（１１）によって表される想定されるノイズが、全ての獲得されたＥＣＧシグナルに対するノイズフィルタとして使用される。フィルトレーションは、パッチ３６から受信したＥＣＧシグナル、及びプローブ２４の電極それぞれにおいてフーリエ変換を実施することと、対応するシグナルの各周波数ｆｐに対して、式（１１）によって与えられたような、Ａｐの値を差し引くこととによって実施される。逆フーリエ変換をついで、減算の結果において実施し、ノイズがそれから除去された、ＥＣＧシグナルを生成する。
【0071】
フローチャート２００の以上の記述は、完全ＭＲＩシーケンスｍが考慮されていること、式（１１）によって与えられたノイズフィルタ｛Ｎｐ｝ｍが、完全シーケンスが動作可能であると同時に獲得されたＥＣＧシグナルをフィルタリングするために使用されること、を推定する。フローチャートの記述は、シーケンスｍの１つ又は複数の相が解析される場合のために、変更すべきところは変更して適合可能であり、その場合、別々のノイズフィルタが計算され、各相の間に獲得されたＥＣＧシグナルの対応する部分に適用されてよい。
【0072】
したがって、シーケンスｍの所定の相ｆに対して、ＭＲＩシーケンスの相に対するパラメータの組は、｛Ｓ｝ｍｆと記述されてよい。ノイズ単離段階２１６において、生成された相ｆのノイズの周波数及び関連した周波数振幅の組に対する式が、記号｛Ｎ｝ｍｆによって示されてよく、【数11】
【0073】
段階２１８にて生成された相ｆに対する関係はついで、【数12】
【0074】
段階２３２にて計算されたノイズフィルタは、【数13】
【0075】
段階２３４において、式（１５）によって表される予想されるノイズが、相ｆの間に獲得されるＥＣＧシグナルの部分に対するノイズフィルタとして使用される。フィルトレーションは、完全ＭＲＩシーケンスｍに対して、上記したのと実質的に同一の様式、すなわち、ＥＣＧシグナルの関連部分を変換することと、式（１５）によって与えられるようなＡｐの値を差し引くことと、ついで除去されたノイズで、ＥＣＧシグナルの部分を再生成するために、逆フーリエ変換を実施することとによって、実施してよい。
【0076】
フローチャート２００の記述は、システム２０が、トレーニングセクション２０２と、続く操作セクション２０４を用いて実行されることを仮定した。しかしながら、本発明の他の実施形態において、別個のトレーイングセクションに対する必要性はなく、これらの実施形態において、システム２０のトレーニングは、「オンザフライ」で実施される。図５及び６を参照して以下で記述した、フローチャート３００及び４００の段階は、どのようにシステム２０が、分離したトレーニングセクションなしに、患者におけるＭＲＩ手順の実施の間に実行されうるのかを例示している。
【0077】
図５は、本発明の他の実施形態にしたがった、実行システム２０にて実施された段階のフローチャート３００である。第一及び第二段階３１０及び３１２は実質的に、それぞれ参照ＥＣＧシグナル及びＭＲＩシーケンス特性を保存する、フローチャート２００の段階２１０及び２１２に関して記述したのと同様である。記録段階３１４は一般に、典型的にＥＣＧ記録は、少数のＭＲＩスキャンに対して、又は単一スキャンに対してであるけれども、段階２１４と同一である。
【0078】
ノイズ単離段階３１６は一般に、段階２１６と同一であり、段階３１０からのＥＣＧ参照シグナルを用いて、段階３１４中で記録されたＥＣＧシグナル中のＭＲＩノイズを単離する。
【0079】
関係形成段階３１８はまた一般に、プロセッサが、段階３１６中で生成されたノイズ周波数と、段階３１２にて保存されたＭＲＩ特性間の関係を生成するように、段階２１８と同一である。生成された関係は、以上の式（９）によって得られ、マトリックス［Ｍ］ｍの推定を与える。
【0080】
フローチャートはついで、マトリックス［Ｍ］ｍの推定にて、ノイズ予測段階３２０に進む。
【0081】
典型的に、他のスキャンに対してスキャナが操作される、反復段階３２２によって示唆されるように、段階３１４、３１６及び３１８が反復される。反復が連続するにつれて、プロセッサは、段階３１８にて生成されたマトリックス［Ｍ］ｍの精度を改善可能である。
【0082】
ノイズ予測段階３２０は一般に、段階２３２と同一であり、プロセッサが、段階３１８にて生成されたマトリックス［Ｍ］ｍの推定を使用して、式（１０）及び（１１）にしたがって、ＥＣＧシグナル上の予測されるノイズを計算する。
【0083】
最終フィルトレーション段階３２４は、実質的に段階２３４と同一であり、ノイズが除去されたＥＣＧシグナルを生成するように、プロセッサが、全てのＥＣＧシグナルから（段階３２０からの）予測されるノイズを差し引く。
【0084】
図６は、本発明の更なる代替実施形態にしたがった、実行システム２０にて実施される段階のフローチャート４００である。第一及び第二段階４１０及び４１２は一般に、段階４１０にて、スキャナが動作不能である一方で、ＥＣＧシグナルを獲得する必要がなく、参照ＥＣＧシグナルの保存が存在しないことを除いて、段階２１０及び２１２と同一である。
【0085】
記録段階４１４は一般に、段階２１４と同一である。典型的に、記録されたＥＣＧシグナルは、パッチ３６及び／又はプローブ２４の固定電極を含む、多数のリードからのものである。
【0086】
ノイズ同定段階４１６において、多数のリード上のシグナルが合計される。合計されたシグナルにより、プロセッサ７０が、ＭＲＩスキャンによって生成されるノイズシグナルを推定可能である。段階２１６に対して以上で記述したように、プロセッサは、式（８）によって与えられたような｛Ｎ｝ｍの推定を作製して、生成されたノイズを特性化する。
【0087】
関係形成段階４１８は一般に、プロセッサが、段階４１６で生成された｛Ｎ｝ｍの推定と、段階３１２にて保存されたＭＲＩ特性との間の関係を生成するように、段階２１８と同一である。生成された関係は、以上の式（９）によって与えられ、マトリックス［Ｍ］ｍの推定を与える。
【0088】
フローチャートがついで、マトリックス［Ｍ］ｍの推定にて、ノイズ推定段階４２０に進む。
【0089】
典型的に、スキャナが更なるスキャンのために操作される、反復段階４２２によって示唆されるように、段階４１４、４１６及び４１８が反復される。反復が連続するにつれて、プロセッサは、段階４１８にて生成されたマトリックス［Ｍ］ｍの推定の精度を改善する。
【0090】
ノイズ予測段階４２０は一般に、段階２３２と同一であり、プロセッサは、段階４１８にて生成されたマトリックス［Ｍ］ｍの推定を使用し、式（１０）及び（１１）にしたがって、ＥＣＧ上の推定されるノイズを計算する。
【0091】
最終フィルトレーション段階４２４は実質的に、段階２３４と同一であり、プロセッサは、ノイズが除去されたＥＣＧシグナルを生成するように、すべてのＥＣＧシグナルから、（段階４２０にて生成される）推定さるノイズを差し引く。
【0092】
２つの別のプロセッサ５６及び５８が動作可能であり、ＥＣＧプロセッサ５８が、フローチャート２００、３００又は４００の段階を実行するために使用される場合、スキャナプロセッサ５６が、シーケンスの同一性、シーケンスの開始及び／又は終了時間、及びシーケンスの相の開始及び／又は終了時間のような、任意の所定のＭＲＩシーケンスに関連した適切な実体を、ＥＣＧプロセッサ５８に情報提供してよい。
【0093】
フローチャート３００及び４００は、最初にトレーニングセクションを有さない、操作システム２０に対する２つの異なる方法を記述している。ＭＲＩスキャンによってＥＣＧシグナル内へ導入されたノイズを推定するために、２つのフローチャートに関して記述された方法を組み合わせうることが評価されるであろう。そのような組み合わせは、ただ一つのフローチャートの工程を適用することと比較して、誘導されたＭＲＩノイズのより迅速、及び／又はより正確な同定を導きうる。
【0094】
したがって、上述の実施形態は一例として引用したものであり、また本発明は上記に具体的に図示及び記載したものに限定されないことは認識されるであろう。むしろ、本発明の範囲には、上で説明した様々な特徴の組み合わせと部分的組合わせの両方、並びにそれらの変形形態及び修正形態が含まれ、これらは、上述の説明を読めば当業者には想到するものであり、従来技術では開示されていないものである。
【0095】
〔実施の態様〕（１） 初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて患者を撮像すると同時に、前記患者から心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、計算された前記フィルタを適用して、前記患者から受信した、ＥＣＧシグナルの更なるセットにおけるノイズを減少させることと、を含む方法。
（２） 前記初期ノイズを同定することが、前記患者を撮像する前にＥＣＧシグナルの参照セットを生成することと、前記ＥＣＧシグナルの初期セットを前記参照セットと比較することとを含む、実施態様１に記載の方法。
（３） 前記初期ノイズを特性化することが、前記初期ノイズ上で、フーリエ変換を実施することを含む、実施態様１に記載の方法。
（４） 前記初期ＭＲＩシーケンスを定義している前記パラメータが、前記初期ＭＲＩシーケンスの電圧及び磁場パルスを記述している変数を含む、実施態様１に記載の方法。
（５） 前記関係を生成することが、前記初期ＭＲＩシーケンスを使用すると同時に、前記変数の値を変更することと、前記値における変化を、前記周波数成分の振幅における変化と比較することと、を含む、実施態様４に記載の方法。
【0096】
（６） 前記初期ＭＲＩシーケンスと、前記続くＭＲＩシーケンスが、共通パラメータを有する共通ＭＲＩシーケンスを含む、実施態様１に記載の方法。（７） 前記共通パラメータが、前記共通ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含み、該変数が、前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスに対して異なる値を有する、実施態様６に記載の方法。
（８） 前記フィルタを計算することが、前記周波数成分の予想される振幅に関して、予想されるノイズを計算することを含み、計算された前記フィルタを適用することが、前記ＥＣＧシグナルの更なるセットと関連した、対応する周波数の振幅から、前記予想される振幅を差し引くことを含む、実施態様１に記載の方法。
（９） 前記ＥＣＧシグナルが、前記患者内に配置したプローブから獲得した、内部ＥＣＧシグナルを含む、実施態様１に記載の方法。
（１０） 前記ＥＣＧシグナルが、前記患者の皮膚から獲得した外部ＥＣＧシグナルを含む、実施態様１に記載の方法。
【0097】
（１１） 前記初期及び続くＭＲＩシーケンスが、共通の複数の相を含み、前記関係を生成することが、前記初期ＭＲＩシーケンスの選択された相を定義しているパラメータと、前記周波数成分との間の相関係を生成することを含み、前記フィルタを計算することが、前記続くＭＲＩ配列の選択された相のパラメータに、前記相関係を適用することによって、前記フィルタを計算することを含む、実施態様１に記載の方法。（１２） 前記ＥＣＧシグナルの初期セットが、前記患者に連結した多数のＥＣＧリード上に発生し、前記初期ノイズを同定することが、前記多数のリード上のシグナルを合計することを含む、実施態様１に記載の方法。
（１３） 前記初期ノイズを同定することが、前記患者を撮像する前に獲得された参照ＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、実施態様１に記載の方法。
（１４） 前記初期ノイズを同定することが、前記患者の撮像の間に、前記初期ノイズを同定することを含む、実施態様１に記載の方法。
（１５） 前記初期ノイズを同定することが、前記患者に接続した固定電極からのＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、実施態様１に記載の方法。
【0098】
（１６） 患者からの心電図（ＥＣＧ）シグナルの初期セットを獲得するために設定された、１つ又は複数の電極と、プロセッサであって、
初期磁気共鳴画像法（ＭＲＩ）シーケンスを用いて前記患者を撮像すると同時に、前記患者からの前記心電図（ＥＣＧ）シグナルの初期セットを受信することと、
前記初期ＭＲＩシーケンスから発生している前記ＥＣＧシグナルの初期セット中の初期ノイズを同定することと、
前記初期ノイズを、その周波数成分に関して特性化することと、
前記初期ＭＲＩシーケンスを定義しているパラメータと、前記周波数成分との間の関係を生成することと、
続くＭＲＩシーケンスのパラメータに対して前記関係を適用することによって、フィルタを計算することと、
前記続くＭＲＩシーケンスを用いて前記患者を撮像すると同時に、前記患者から受信した、ＥＣＧシグナルの更なるセット中のノイズを減少させるために、計算された前記フィルタを適用することと、のために設定された、プロセッサと、を含む装置。
（１７） 前記初期ノイズを同定することが、前記患者を撮像することの前に、ＥＣＧシグナルの参照セットを生成することと、前記ＥＣＧシグナルの初期セットと、前記参照セットを比較することと、を含む、実施態様１６に記載の装置。
（１８） 前記初期ノイズを特性化することが、前記初期ノイズの上でフーリエ変換を実施することを含む、実施態様１６に記載の装置。
（１９） 前記初期ＭＲＩシーケンスを定義している前記パラメータが、前記初期ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含む、実施態様１６に記載の装置。
（２０） 前記関係を生成することが、前記初期ＭＲＩシーケンスを使用すると同時に前記変数の値を変更することと、前記値における変化と、前記周波数成分の振幅における変化とを比較することと、を含む、実施態様１９に記載の装置。
【0099】
（２１） 前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスが、共通パラメータを有する共通ＭＲＩシーケンスを含む、実施態様１６に記載の装置。（２２） 前記共通パラメータが、前記共通ＭＲＩシーケンスの電圧と磁場パルスを記述している変数を含み、該変数が、前記初期ＭＲＩシーケンスと前記続くＭＲＩシーケンスに対する異なる値を有する、実施態様２１に記載の装置。
（２３） 前記フィルタを計算することが、前記周波数成分の予想される振幅に関して、予想されるノイズを計算することを含み、計算された前記フィルタを適用することが、前記ＥＣＧシグナルの更なるセットと関連した対応する周波数の振幅から前記予想される振幅を差し引くことを含む、実施態様１６に記載の装置。
（２４） 前記１つ又は複数の電極が前記患者内に位置し、前記ＥＣＧシグナルに、前記１つ又は複数の電極から獲得された初期ＥＣＧシグナルが含まれる、実施態様１６に記載の装置。
（２５） 前記１つ又は複数の電極が、前記患者の皮膚上に位置し、前記ＥＣＧシグナルが、前記１つ又は複数の電極から獲得された外部ＥＣＧシグナルを含む、実施態様１６に記載の装置。
【0100】
（２６） 前記初期及び続くＭＲＩシーケンスが、共通の複数の相を含み、前記関係を生成することが、前記初期ＭＲＩシーケンスの選択された相を定義しているパラメータと、前記周波数成分との間の相関係を生成することを含み、前記フィルタを計算することが、前記続くＭＲＩシーケンスの選択された相のパラメータに対して、前記相関係を適用することによって前記フィルタを計算することを含む、実施態様１６に記載の装置。（２７） 前記ＥＣＧシグナルの初期セットを獲得するために設定された多数のＥＣＧリードを含み、前記初期ノイズを同定することが、前記多数のＥＣＧリード上のシグナルを合計することを含む、実施態様１６に記載の装置。
（２８） 前記初期ノイズを同定することが、前記患者を撮像する前に獲得した参照ＥＣＧシグナルを用いて、前記初期ノイズを同定することを含む、実施態様１６に記載の装置。
（２９） 前記初期ノイズを同定することが、前記患者の撮像の間に、前記初期ノイズを同定することを含む、実施態様１６に記載の装置。
（３０） 前記１つ又は複数の電極が、前記患者に接続した固定電極を含み、前記初期ノイズを同定することが、前記固定電極からのＥＣＧシグナルを用いて前記初期ノイズを同定することを含む、実施態様１６に記載の装置。
【外国語明細書】