(57)【特許請求の範囲】
【請求項１】
ＭＲ装置の検査体積内に配置された対象のＭＲ撮像の方法において、
２以上の空間的に離れた画像スライスを同時に励起するマルチスライスＲＦパルスを有する撮像シーケンスを前記対象に行うステップと、
ＭＲ信号を取得するステップであって、前記ＭＲ信号が、前記検査体積内で異なる空間感度プロファイルを持つＲＦコイルのセットにより並列に受信される、当該取得するステップと、
前記取得されたＭＲ信号から各画像スライスに対するＭＲ画像を再構成するステップであって、異なる画像スライスからのＭＲ信号寄与が、前記ＲＦコイルの前記空間感度プロファイルに基づいて分離され、前記画像スライスからのＭＲ信号寄与が、サイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起された領域からのＭＲ信号寄与から分離され、前記分離されたサイドバンドアーチファクトが、前記ＲＦコイルの前記空間感度プロファイルに基づいて前記再構成されたＭＲ画像において抑制される、当該再構成するステップと、
を有する方法。
【請求項２】
前記画像スライスからのＭＲ信号寄与が、前記マルチスライスＲＦパルスの励起スペクトルに関する事前情報を考慮に入れることなしに前記サイドバンドアーチファクトから分離される、請求項１に記載の方法。
【請求項３】
前記画像スライスからのＭＲ信号寄与が、前記取得されたＭＲ信号の信号モデルを使用することにより前記サイドバンドアーチファクトから分離され、前記信号モデルが、（ｉ）前記画像スライス及び（ｉｉ）前記マルチスライスＲＦパルスの前記１以上のサイドバンドにより励起された領域からの信号寄与を有する、請求項１又は２に記載の方法。
【請求項４】
前記サイドバンドアーチファクトが、一組の線形方程式を解くことにより決定され、前記取得されたＭＲ信号に対する信号寄与（ｉ）及び（ｉｉ）の比が、反復的に調整される、請求項３に記載の方法。
【請求項５】
前記マルチスライスＲＦパルスの前記サイドバンドにより励起された領域は、サイドバンド周波数、すなわち、前記マルチスライスＲＦパルスの基本周波数のより高次の高調波が、前記撮像シーケンスのスライス選択磁場勾配の存在下で共鳴する場所である、請求項１乃至４のいずれか一項に記載の方法。
【請求項６】
前記マルチスライスＲＦパルスが、位相変調され、位相変調スキームは、位相サイクルが各画像スライスのＭＲ信号に適用されるように変化する位相シフトを有する、請求項１乃至５のいずれか一項に記載の方法。
【請求項７】
前記位相シフトが、位相符号化ステップごとに線形にインクリメントされる、請求項６に記載の方法。
【請求項８】
前記ＭＲ信号が、前記画像スライスの面内方向においてアンダーサンプリングで取得され、前記ＭＲ画像が、ＳＥＮＳＥ、ＳＭＡＳＨ又はＧＲＡＰＰＡのような並列画像再構成アルゴリズムを使用して再構成される、請求項１乃至７のいずれか一項に記載の方法。
【請求項９】
請求項１乃至８に記載の方法を実行するＭＲ装置において、前記ＭＲ装置が、検査体積内に一様な静止磁場を生成する少なくとも１つの主磁石コイルと、前記検査体積内で異なる空間方向において切り替え磁場勾配を生成する複数の勾配コイルと、異なる空間感度プロファイルを持つＲＦコイルのセットと、ＲＦパルス及び切り替え磁場勾配の時間的並びを制御する制御ユニットと、再構成ユニットとを含み、前記ＭＲ装置が、
２以上の空間的に離れた画像スライスを同時に励起するマルチスライスＲＦパルスを有する撮像シーケンスを前記対象に行うステップと、
ＭＲ信号を取得するステップであって、前記ＭＲ信号が、前記検査体積内で異なる空間感度プロファイルを持つＲＦコイルのセットにより並列に受信される、当該取得するステップと、
前記取得されたＭＲ信号から各画像スライスに対するＭＲ画像を再構成するステップであって、異なる画像スライスからのＭＲ信号寄与が、前記ＲＦコイルの前記空間感度プロファイルに基づいて分離され、前記画像スライスからのＭＲ信号寄与が、サイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起された領域からのＭＲ信号寄与から分離され、前記分離されたサイドバンドアーチファクトが、前記ＲＦコイルの前記空間感度プロファイルに基づいて前記再構成されたＭＲ画像において抑制される、当該再構成するステップと、
を実行するように構成される、ＭＲ装置。
【請求項１０】
ＭＲ装置上で実行されるコンピュータプログラムにおいて、
２以上の空間的に離れた画像スライスを同時に励起するマルチスライスＲＦパルスを有する撮像シーケンスを生成する命令と、
ＭＲ信号を取得する命令と、
前記取得されたＭＲ信号から各画像スライスに対するＭＲ画像を再構成する命令であって、異なる画像スライスからのＭＲ信号寄与が、ＲＦコイルのセットの空間感度プロファイルに基づいて分離され、前記画像スライスからのＭＲ信号寄与が、サイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起された領域からのＭＲ信号寄与から分離され、前記分離されたサイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起された領域からのＭＲ信号寄与が、前記ＲＦコイルの前記空間感度プロファイルに基づいて前記再構成されたＭＲ画像において抑制される、当該再構成する命令と、
を有する、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴（ＭＲ）撮像の分野に関する。これは、対象のＭＲ撮像の方法に関係する。本発明は、ＭＲ装置及びＭＲ装置上で実行されるコンピュータプログラムにも関する。
【背景技術】
【０００２】
二次元又は三次元画像を形成するために磁場と核スピンとの間の相互作用を使用する画像形成ＭＲ方法は、軟組織の撮像に対して、多くの点で他の撮像方法より優れており、電離放射線を必要とせず、通常は非侵襲的であるので、特に医療診断の分野で、現在幅広く使用されている。
【０００３】
一般的なＭＲ方法によると、検査されるべき患者の体は、測定が基づく座標系の軸（通常はｚ軸）を同時に規定する方向を持つ強力で一様な磁場（Ｂ0磁場）内に配置される。磁場は、磁場強度に依存して個別の核スピンに対して異なるエネルギレベルを生成する。これらのエネルギレベル間の遷移は、規定された周波数（いわゆるラーモア周波数又はＭＲ周波数）の交流電磁場（ＲＦ場、Ｂ1場とも称される）の印加により励起（スピン共鳴）されることができる。巨視的な視点から、個別の核スピンの分布は、磁化がｚ軸について歳差運動を実行するように、適切な周波数の電磁パルス（ＲＦパルス）の印加により平衡状態から偏向されることができる全体的な磁化を生成する。歳差運動は、フリップ角と称される開口の角度を持つ円錐の表面を描く。フリップ角の大きさは、印加される電磁パルスの強度及び持続時間に依存する。
【０００４】
ＲＦパルスの終了後に、磁化は、元の平衡状態まで緩和し、ｚ方向の磁化が、第１の時定数Ｔ1（スピン格子又は縦緩和時間）で再び増大され、ｚ方向に垂直な方向における磁化が、第２の時定数Ｔ2（スピン‐スピン又は横緩和時間）で緩和する。磁化の変化は、磁化の変化がｚ軸に垂直な方向において測定されるような形でＭＲ装置の検査体積内に配置及び配向される１以上の受信ＲＦコイルを用いて検出されることができる。
【０００５】
体内の空間解像度を実現するために、３つの主軸に沿って延びる線形磁場勾配が、一様な磁場に重畳され、スピン共鳴周波数の線形空間依存性を生じる。受信コイルにおいて取得される信号は、体内の異なる場所に関連付けられることができる異なる周波数の成分を含む。ＲＦコイルをにより取得されたＭＲ信号データは、空間周波数領域に対応し、ｋ空間データと称される。ｋ空間データのセットは、フーリエ変換又は他の適切な再構成アルゴリズムを用いてＭＲ画像に変換される。
【０００６】
ＭＲ取得を加速するための並列取得技術は、当技術分野において何年も前から既知である。このカテゴリにおける方法は、ＳＥＮＳＥ（感度符号化）、ＳＭＡＳＨ（空間高調波の同時取得）及びＧＲＡＰＰＡ（一般化自動較正部分並列取得）である。ＳＥＮＳＥ、ＳＭＡＳＨ及びＧＲＡＰＰＡ並びに他の並列取得技術は、並列で複数のＲＦ受信コイルから得られたアンダーサンプリングされたｋ空間データ取得を使用する。これらの方法において、複数のコイルからの（複雑な）信号データは、最終的に再構成されるＭＲ画像においてアンダーサンプリングアーチファクト（エイリアシング）を抑制するような形で複雑な重み付けで結合される。このタイプの複雑なアレイ信号結合は、時々、空間フィルタリングと称され、（ＳＭＡＳＨ及びＧＲＡＰＰＡのように）ｋ空間領域又は（ＳＥＮＳＥのように）画像領域において実行される結合、及びハイブリッドである方法を含む。
【０００７】
Larkman他（Journal of Magnetic Resonance Imaging, 13, 313-317, 2001）は、スキャン効率を増大するようにマルチスライス撮像の場合にスライス方向においても感度符号化を適用することを提案する。Breuer他（Magnetic Resonance in Medicine, 53, 684-691, 2005）は、この基本的なアイデアを改良し、「高加速における並列撮像結果における制御されたエイリアシング（controlled aliasing in parallel imaging results in higher acceleration）」（ＣＡＩＰＩＲＩＮＨＡ）と称されるアプローチを提案する。この技術は、マルチスライス取得中に各個別のスライスのエイリアシングアーチファクトの外形を修正し、後の並列画像再構成プロシージャを改良する。したがって、ＣＡＩＰＩＲＩＮＨＡは、純粋な後処理アプローチのみを使用する他のマルチスライス並列撮像と比べて効率的な並列マルチスライス撮像技術である。ＣＡＩＰＩＲＩＮＨＡにおいて、任意の厚さ及び距離の複数のスライスが、位相変調マルチスライスＲＦパルスの使用で同時に励起される。取得されるＭＲ信号データは、同時にサンプリングされ、互いに対してシフトされたように見える重畳されたスライス画像を生じる。エイリアスされたスライス画像のシフトは、フーリエシフト定理によってＲＦパルスの位相変調スキームにより制御される。位相符号化ステップごとに、マルチスライスＲＦパルスは、各スライスのＭＲ信号に対して個別の位相シフトを加える。関与するスライスの個別の信号寄与を分離する逆再構成問題の数値的条件は、このシフトを使用して改良される。ＣＡＩＰＩＲＩＮＨＡは、使用されるＲＦ受信コイルのコイル感度が撮像されるべき個別のスライスにおいて劇的に異なることのないように、スライスが互いにかなり近い場合でも、重畳されたスライス画像の分離を改良する可能性がある。
【０００８】
しかしながら、従来の並列マルチスライス撮像アプローチは、限界を持つ。複数の周波数におけるＭＲ信号が、マルチスライス（又はマルチ周波数）ＲＦパルスにより同時に励起される場合、いわゆるサイドバンドアーチファクトが、再構成画像において生じる。これらのアーチファクトは、マルチスライスＲＦパルスの１以上のサイドバンドにより意図せずに励起された領域からのＭＲ信号により引き起こされる。サイドバンド周波数は、それぞれのＲＦパルスの基本（メインバンド）周波数のより高次の高調波でありうる。マルチスライスＲＦパルスのこのようなサイドバンドは、使用されるＭＲ装置のハードウェア制約、例えばＲＦ増幅器の非線形性により、実際には回避不可能である。サイドバンドアーチファクトの特性は、ＲＦコイル構成の個別のロード、検査体積内のＢ1分布、及びマルチバンド励起に関与する基本周波数に依存する。
【発明の概要】
【発明が解決しようとする課題】
【０００９】
先行する記載から、改良された並列マルチスライスＭＲ撮像技術に対する必要性が存在することが、容易に理解される。本発明の目的は、サイドバンドアーチファクトの効率的な抑制を持つマルチスライスＭＲ撮像を可能にすることである。
【課題を解決するための手段】
【００１０】
本発明によると、ＭＲ装置の検査体積内に配置された対象のＭＲ撮像の方法が、開示される。前記方法は、
‐２以上の空間的に離れた画像スライスを同時に励起するマルチスライスＲＦパルスを有する撮像シーケンスを前記対象に行うステップと、
‐ＭＲ信号を取得するステップであって、前記ＭＲ信号は、前記検査体積内で異なる空間感度プロファイルを持つＲＦコイルのセットにより並列に受信される、当該取得するステップと、
‐前記取得されたＭＲ信号から各画像スライスに対するＭＲ画像を再構成するステップであって、異なる画像スライスからのＭＲ信号寄与が、前記ＲＦコイルの空間感度プロファイルに基づいて分離され、サイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起される領域からのＭＲ信号寄与が、前記ＲＦコイルの空間感度プロファイルに基づいて前記再構成されたＭＲ画像において抑制される、当該再構成するステップと、
を有する。
【００１１】
本発明の主旨は、純粋に（例えば、既知のＳＥＮＳＥアルゴリズムのような）並列画像再構成アルゴリズムを使用することにより最終的に再構成されたスライス画像においてサイドバンドアーチファクトを抑制することである。このために、前記画像スライスからのＭＲ信号寄与が、前記マルチスライスＲＦパルスの励起スペクトルに関する事前情報を考慮に入れることなしに本発明によって前記サイドバンドアーチファクトから分離される。換言すると、本発明は、本質的にサイドバンドアーチファクトを持たないＭＲ画像を再構成することができるために（メインバンド周波数に関するサイドバンド成分の振幅のような）前記使用されたＲＦパルスのサイドバンドスペクトルの詳細に関して仮定する又は知ることを要求しない。なされる必要のある唯一の仮定は、前記マルチスライスＲＦパルスのサイドバンドにより励起される領域の場所に関する。これらは、前記サイドバンド周波数（典型的には、それぞれのＲＦパルスの基本周波数のより高次の高調波）が、それぞれ印加されたスライス選択磁場勾配の存在下で共鳴する場所である。本発明により達成される抑制は、前記サイドバンドアーチファクトの正確な先験的な場所に対してあまり敏感ではない。これは、反復的アプローチにおいて、少数の反復のみが、収束を達成するのに必要とされるように見える。前記サイドバンドアーチファクトのアンフォールディングは、適切な場所における実際のサイドバンド寄与の分離に対する反復的収束を生じる内的一貫性を含む。このアンフォールディングは、サイドバンド信号寄与に対して支配的である画像スライスのＭＲ信号寄与を使用する。
【００１２】
本発明の好適な実施例において、前記画像スライスからのＭＲ信号寄与は、前記取得されたＭＲ信号の信号モデルを使用することにより前記サイドバンドアーチファクトから分離され、前記信号モデルは、（ｉ）前記画像スライス及び（ｉｉ）前記マルチスライスＲＦパルスの前記１以上のサイドバンドにより（潜在的に）励起された画像スライスの外の領域からの信号寄与を有する。可能な実施例において、前記信号モデルのパラメータは、前記マルチスライスＲＦパルスのサイドバンドエネルギがメインバンドエネルギより大幅に小さいと仮定されることができるので、一般に１より大幅に小さいこれら２つの信号寄与の比であってもよい。前記サイドバンドアーチファクト、すなわち、前記画像スライスの外の領域からのＭＲ信号寄与は、例えば、前記２つの信号寄与の比が反復的に調整される従来のＳＥＮＳＥ再構成スキームのように、一組の線形方程式を解くことにより再構成される（及びこれにより前記最終的に再構成されたＭＲ画像において抑制／減算される）ことができる。実際的な場合に、２乃至５の反復が、収束を達成するのに十分である。スペクトルの粗い特性のみが、取得された磁気共鳴信号をモデル化するのに必要とされてもよい。このモデルは、少なくとも反復的アプローチにおいて自動的に決定されることができる最初のパラメータとして、メインバンド及びサイドバンド寄与の比を含む。
【００１３】
本発明の好適な実施例において、前記マルチスライスＲＦパルスは、位相変調され、位相変調スキームは、位相サイクルが各画像スライスのＭＲ信号に加えられるような、変化する位相シフトを有する。このようにして、本発明の技術は、（上記の）既知のＣＡＩＰＩＲＩＮＨＡスキームと組み合わせられる。好ましくは、前記位相シフトは、位相符号化ステップごとに線形にインクリメントされ、個別の位相インクリメントは、各画像スライスに加えられる。このようにして、各スライス画像の個別のシフトは、フーリエシフト定理によって前記ＲＦパルスの位相変調スキームにより制御される。
【００１４】
本発明の他の好適な実施例において、前記ＭＲ信号は、前記画像スライスの面内方向においてアンダーサンプリングで取得される。前記画像スライスのＭＲ画像は、この場合、ＳＥＮＳＥ、ＳＭＡＳＨ又はＧＲＡＰＰＡのような、それ自体が既知である並列画像再構成アルゴリズムにより再構成されることができる。
【００１５】
ここまでに記載された本発明の方法は、検査体積内に一様な静止磁場を生成する少なくとも１つの主磁石コイルと、前記検査体積内で異なる空間方向において切り替え磁場勾配を生成する複数の勾配コイルと、並列に体からのＭＲ信号を受信し、異なる空間感度プロファイルを持つＲＦコイルのセットと、ＲＦパルス及び切り替え磁場勾配の時間的な並びを制御する制御ユニットと、再構成ユニットとを含むＭＲ装置を用いて実行されることができる。本発明の方法は、例えば、前記ＭＲ装置の前記再構成ユニット及び／又は前記制御ユニットの対応するプログラミングにより実施されることができる。
【００１６】
本発明の方法は、現在、臨床的に使用されているほとんどのＭＲ装置で有利に実行されることができる。このために、単に、前記ＭＲ装置が上で説明された本発明の方法ステップを実行するように制御されるコンピュータプログラムを使用することが、必要である。前記コンピュータプログラムは、データ担体、又は前記ＭＲ装置の前記制御ユニットにおけるインストールのためにダウンロードされるようにデータネットワークのいずれに存在してもよい。
【００１７】
添付の図面は、本発明の好適な実施例を開示する。理解されるべきは、しかしながら、図面が、図示の目的であり、本発明の限定の既定として設計されているわけではないことである。
【図面の簡単な説明】
【００１８】
【図１】本発明の方法を実行するＭＲ装置を示す。
【図２ａ】従来のマルチスライス取得スキームを使用して取得されたＭＲスライス画像を示す。
【図２ｂ】本発明によるサイドバンドアーチファクトの抑制を用いてＭＲスライス画像を示す。
【図２ｃ】図２ａのＭＲスライス画像に対するサイドバンドアーチファクトの寄与を示す。
【発明を実施するための形態】
【００１９】
図１を参照すると、ＭＲ装置１が示される。前記装置は、実質的に一様な、時間的に一定の主磁場が検査体積を通るｚ軸に沿って作成されるような超伝導又は常伝導主磁石コイル２を有する。
【００２０】
磁気共鳴生成及び操作システムは、ＭＲ撮像を実行するように、核磁気スピンを反転又は励起し、磁気共鳴を誘起し、磁気共鳴をリフォーカスし、磁気共鳴を操作し、前記磁気共鳴を空間的に及び他の形で符号化し、スピンを飽和させるように等、一連のＲＦパルス及び切り替え磁場勾配を印加する。
【００２１】
より具体的には、勾配パルス増幅器３は、前記検査体積のｘ、ｙ及びｚ軸に沿った全身勾配コイル４、５及び６の選択されたものに電流パルスを印加する。デジタルＲＦ周波数送信器７は、前記検査体積内にＲＦパルスを送信するように、送受信スイッチ８を介して、全身体積ＲＦコイル９にＲＦパルス又はパルスパケットを送信する。典型的なＭＲ撮像シーケンスは、互いと一緒に取られる短い持続時間のＲＦパルスセグメントのパケットからなり、いかなる印加された磁場勾配も、核磁気共鳴の選択された操作を達成する。前記ＲＦパルスは、飽和させる、共鳴を励起する、磁化を反転させる、共鳴をリフォーカスする、又は共鳴を操作し、前記検査体積内に配置された体１０の一部を選択するのに使用される。
【００２２】
並列撮像を用いる体１０の限定領域のＭＲ画像の生成に対して、局所ＲＦコイル１１、１２、１３のセットは、撮像に対して選択された領域に隣接するように配置される。
【００２３】
結果として生じるＭＲ信号は、ＲＦコイル１１、１２、１３により取得され、好適には前置増幅器（図示されない）を含む受信器１４により復調される。受信器１４は、送受信スイッチ８を介してＲＦコイル９、１１、１２及び１３に接続される。
【００２４】
ホストコンピュータ１５は、エコープラナー撮像（ＥＰＩ）、エコー体積撮像、勾配及びスピンエコー撮像、高速スピンエコー撮像等のような、複数のＭＲ撮像シーケンスのいずれかを生成するように勾配パルス増幅器３及び送信器７を制御する。選択されたシーケンスに対して、受信器１４は、各ＲＦ励起パルスに続いて急速に連続して単一の又は複数のＭＲデータラインを受信する。データ取得システム１６は、受信された信号のアナログ・デジタル変換を実行し、各ＭＲデータラインを更なる処理に適したデジタルフォーマットに変換する。近年のＭＲ装置において、データ取得システム１６は、生画像データの取得に特化される別個のコンピュータである。
【００２５】
最終的に、デジタル生画像データが、フーリエ変換又は他の適切な再構成アルゴリズムを適用する再構成プロセッサ１７により画像表現に再構成される。前記ＭＲ画像は、患者を通る平面スライス、平行な平面スライスのアレイ、又は三次元体積等を表してもよい。前記画像は、次いで、例えば結果として生じるＭＲ画像の人間可読ディスプレイを提供するビデオモニタ１８により、画像表現のスライス、投影又は他の部分を視覚化に対する適切なフォーマットに変換するためにアクセスされてもよい画像メモリに記憶される。
【００２６】
図１を参照し続け、更に図２を参照すると、本発明の撮像アプローチの一実施例が説明される。
【００２７】
本発明によると、患者の体１０は、マルチスライスＲＦパルスを有する撮像シーケンスを受け、これにより２以上の空間的に離れた画像スライス内の核スピンが同時に励起される。前記撮像シーケンスにより生成された前記ＭＲ信号は、異なる空間感度プロファイルを持つＲＦコイル１１、１２、１３により並列に取得される。従来のマルチスライス技術と同様に、ＭＲ画像は、前記取得されたＭＲ信号から各画像スライスに対して再構成され、異なる画像スライスからのＭＲ信号寄与は、ＲＦコイル１１、１２、１３の（既知の）空間感度プロファイルに基づいて分離される。実際には従来のＳＥＮＳＥアンフォールディングアルゴリズムに対応する、前記画像スライスの分離に対して適用されるアルゴリズムは、以下に詳細に記載される。
【００２８】
まず、Ｎの異なる画像スライスにわたって、Ｍの受信コイルの各々の取得されたＭＲ信号ｍにおいて１つの場所ｘに寄与する全ての画像場所ｘiを考慮する。これは、行列ベクトル表記法において、
Ｓｐ＝ｍ
と記載されることができる。
【００２９】
ここで、ベクトルｍは、前記Ｎの異なる画像スライスの感度で重み付けされた信号寄与ｐi(ｘi)の線形結合としてＭのＲＦコイル１１、１２、１３の各々における取得されたＭＲ信号ｍj(ｘ)を示し、行列Ｓは、位置ｘiにおけるスライスＩ及びコイルｊに対するコイル感度であるＳijを持つ（Ｎ×Ｍ）の感度行列を示す。符号化行列の逆を含むこの連立方程式を解くことは、対応するＮのスライス固有ＭＲ信号を含むベクトルｐを生じる。
(ＳHＳ)-1ＳHｍ＝ｐ
【００３０】
行列(ＳHＳ)-1ＳHは、Ｓの疑似逆であり、そのノルムは、ＭＲ信号取得から最終画像までの誤差伝搬を記述する。このノルムは、良好な条件の場合に小さい。
【００３１】
本発明によると、サイドバンドアーチファクト、すなわち、前記マルチスライスＲＦパルスの１以上のサイドバンドにより励起された領域からのＭＲ信号寄与は、ＲＦコイル１１、１２、１３の空間感度プロファイルに基づいて前記再構成されたＭＲ画像において抑制される。これを達成するために、場所xiにおけるＮのスライス固有メインバンドＭＲ信号寄与を含むベクトルｐ及び追加的にＬのサイドバンドＭＲ信号寄与、すなわち、前記マルチスライスＲＦパルスの前記サイドバンド周波数成分により潜在的に励起された前記画像スライスの外の領域からのＭＲ信号寄与を含むベクトルｐ'を有する信号モデルが、使用される。このモデルを用いて、前記Ｍの受信コイルの各々による前記取得されたＭＲ信号は、
0006332891.tif 000002
として書かれる行列であることができる。
【００３２】
ここで、行列Ｓは、（Ｎ＋Ｌ）×Ｍの感度行列を示し、Ｓijは、メインバンド寄与（ｉ＝１．．Ｎ）及びサイドバンド寄与（ｉ＝Ｎ＋１．．Ｎ＋Ｌ）並びにコイルｊに対するコイル感度である。この連立方程式は、主に既知の正則化ＳＥＮＳＥフレームワークを使用することにより解かれることができる。
0006332891.tif 000003
【００３３】
ここで、Ｒ／Ｒ'は、正則化行列であり、σは、メインバンド及びサイドバンド寄与の比を表す。
0006332891.tif 000004
【００３４】
パラメータσは、ユーザパラメータとして得られることができ、又は上の式を反復的に解くことにより自動的に決定されることができ、σは以下のように更新される。
σ＝average(ｐ'√Ｒ')／average(ｐ√Ｒ)
【００３５】
ここで、"average"は、全ての画像ボクセル又は所定の位置の周りの所定の領域内の画像ボクセルにわたる平均と理解されるべきである。収束は、実際には、少数の２乃至５の反復後に達成されるべきである。より一般的なモデルにおいて、パラメータσは、パラメータセットσ1，σ2，...σLが使用されることができるように、各サイドバンドに対して異なっていてもよい。ベクトルｐの解は、サイドバンドアーチファクトを持たないＮのＭＲスライス画像を表す。
【００３６】
これは、図２に示される。図２ａは、マルチスライス励起を追使用してＲＦコイルのセットにより並列にファントムから取得されたＭＲ信号データから従来の形式で再構成されたＭＲスライス画像を示す。白い矢印は、ＭＲ画像の中心における強力なサイドバンドアーチファクトを示す。図２ｂは、本発明によって前記サイドバンドアーチファクトを抑制する上記の修正ＳＥＮＳＥフレームワークを使用して同じＭＲ信号データから再構成されたＭＲスライス画像を示す。図に見られるように、前記サイドバンドアーチファクトは、図２ｂにおいておおよそ完全に取り除かれている。図２ｃは、上のフレームワークを使用してメインバンド信号（ｐ）から分離されたサイドバンドアーチファクト（ｐ'）を示す。
【００３７】
本発明の一実施例において、ＳＥＮＳＥは、適切な減少係数を使用して面内位相符号化方向において追加的に使用されてもよい。
【図１】
0006332891.tif 000005
【図２ａ）】
0006332891.tif 000006
【図２ｂ）】
0006332891.tif 000007
【図２ｃ）】
0006332891.tif 000008