図1を参照して、磁気共鳴(MR)スキャナ10は、静磁場(B0)を生成する主磁気や磁場勾配コイルセットなどのコンポーネント(図示せず)を収容又は保持するハウジング12と、図示されたMRスキャナ10のケースにおいて、MRスキャナ10のボア16内にあるイメージング領域の内外に平行移動可能な被検者ベッドなどのMR被検者ロードシステム14とを有する。図示されたMRスキャナ10は、Koninklijke Philips Electronics N.V.(オランダのEindhoven)から入手可能なAchievaTM MRスキャナである。しかしながら、実質的に任意のMRスキャナが利用可能である。
N個(Nは2以上の整数)のRF送信チャネル20が図式的に示される図1に示されるように、複数のRF送信チャネル20が設けられる。複数のRF送信チャネル20は、B1送信フィールドと示されることもあるRF送信フィールドを生成するため、マルチチャネル送信モードにおいて動作可能である。B1送信フィールドのRF周波数は、好ましくは、磁気共鳴周波数であるか、又はその近傍である。所与の静磁場(B0)について、磁気共鳴周波数は、静磁場強度(|B0|)と、核磁気共鳴を受けることが意図される核の性質である磁気回転定数(γ)との積により与えられる。
複数のRF送信チャネル20は、様々に実現可能である。例えば、ある実施例では、複数のRF送信チャネル20は、RF送信チャネルの個数Nが2になるように、独立に駆動されるI及びQポートを有する単一のバードケージタイプの容積RFコイルとして実現される。他の実施例では、複数のRF送信チャネル20は、N個の独立した表面コイル、変性した全身用RFコイルのN個の分離されたロッド又はラングなどのN個の独立したコイル要素のセットとして実現される。これらの実施例では、N個の独立したコイル要素は、例えば、別々に収容されたコイル要素、電気的に絶縁されるが、共通のハウジングに物理的に収容されるコイル要素(例えば、専用のN要素コイルアレイアセンブリなど)などとして様々に構成されてもよい。
さらに、1以上の磁気共鳴受信コイルが設けられる。ある実施例では、複数のRF送信チャネル20の1つ、一部又はすべてが、磁気共鳴を受けるための受信モードに適切にスイッチされる送受信コイルとして構成される。他の実施例では、複数のRF送信チャネル20から分離した1以上の磁気共鳴受信コイル20(図示せず)が、磁気共鳴受信処理を実行するため設けられる。
図1を続けて参照して、MRシステムはさらに、放射線技師や他のユーザがMRスキャナ10にMRイメージングデータを取得させ、MR被検者ロードシステム14を介しイメージング被検者の自動的なロード及びアンロードなどの他の機構を実行させるためにMRスキャナ10とインタフェースをとることが可能なMRシステムコントローラユーザインタフェースモジュール22を有する。
典型的なイメージングシーケンスでは、イメージング対象の被検者は、ロードシステム14を用いてボア16のイメージング領域にロードされ、複数のRF送信チャネル20のRF送信チャネルが、被検者の磁気共鳴を励起するためマルチチャネル送信モードにおいて起動され、磁界勾配コイルが、磁気共鳴を空間的に制限及び/又は符号化するか、または操作するため、磁気共鳴励起の前、期間中及び/又は後に作動され、磁気共鳴がMR受信コイルを介し受信され、所得されたMRデータストレージ24に格納される。取得されたMRデータは、再構成されたMRイメージストレージ28に格納される1以上の再構成されたMRイメージを生成するため、MRイメージ再構成モジュール26により適切に再構成される。再構成モジュール26は、MRイメージングデータの取得中に使用される空間符号化により動作する再構成アルゴリズムを利用する。例えば、MRイメージングデータがデカルト符号化を用いてkスペースサンプルとして取得される場合、フーリエ変換ベースの再構成アルゴリズムは、再構成モジュール26によって適切に利用されてもよい。
この例示的なイメージングシーケンスでは、複数のRF送信チャネルのRF送信チャネルは、被検者において磁気共鳴を励起するため、マルチチャネル送信モードにおいて起動される。マルチチャネル送信モードでは、各RF送信チャネルは、RF励起振幅及び位相に関して独立に制御される。従って、N個のRFチャネルについて、2N個の独立に調整可能なパラメータがある。実質的に(空間的に)均一なB1送信フィールドを提供し、所望のRF送信パワーのB1送信フィールドを提供するため、これら2N個のパラメータを調整することが所望される。実質的に均一のB1送信フィールドを提供するため、RFチャネルを調整することは、RFシムとして知られる。所望のRF送信パワーを提供するためのRFチャネルの調整は、典型的には、ターゲットの90度のフリップ角度などの被検者において所望されるフリップ角度を提供するため、又は比吸収率(SAR)又は他の被検者の安全性の指標などを制限するなどのために実行される。2N個のマルチチャネル送信パラメータの所与のセットに対するB1送信フィールドの均一性は、被検者が受けるイメージングの電気的及び/又は磁気的感受性性質によって実質的に影響を受ける可能性があり、これにより、“最適”な送信パラメータは一般に被検者に固有のものとなる。静的磁界(B0)が増加するに従って、B1送信フィールドに対する被検者の影響は増加する傾向がある。
図1を参照し続けて、MRシステムはさらに、個々のRF送信チャネルについて取得したB1マップに基づき、複数のRF送信チャネル20のRF送信チャネルのRF振幅及び位相を最適化するRFシムRF送信パワー最適化モジュール30を有する。利用されるB1マップは、好ましくは、B1送信フィールドに対する上述された被検者のロード効果を説明するため、必ずしも必要ではないが、ロードされた被検者によって取得される。最適化された振幅及び位相は、被検者のイメージング中にMRシステムコントローラユーザインタフェースモジュール22による収集及び利用のため、RF送信チャネル振幅及び位相パラメータストレージ32に格納される。
処理モジュール22,26,30は、図1の例示的な実施例では、コンピュータ42のプロセッサであるデジタルプロセッサ40により適切に実現される。デジタルプロセッサ40は、マルチコアマイクロプロセッサ、マイクロプロセッサ及び連動するグラフィカル処理ユニット(GPU)又は数値演算コプロセッサなどのケースなど、複数のプロセッサであってもよいことが理解されるべきである。さらに、デジタルプロセッサ40は、コンピュータの一部でない専用プロセッサとして構成されてもよい。さらに、各種処理モジュール22,26,30は、異なるプロセッサにより実現されてもよく、及び/又は非デジタルプロセッサコンポーネントを含むものであってもよい。例えば、再構成モジュール26は、アナログパイプラインコンポーネントを含むものであってもよい。MRシステムコントローラユーザインタフェースモジュール22のユーザインタフェースコンポーネントは、MRスキャナ構成の再構成されたイメージを表示し、又は他のユーザにより知覚可能な出力を提供するためのコンピュータ42の図示されたディスプレイ44などの適切なユーザインタフェースハードウェア、又はユーザ入力を受信するためのマウス、トラックボール、タッチセンシティブスクリーンなどの他のユーザ入力装置にアクセスする。各種データストレージコンポーネント24,28,32は、ハードディスクドライブ、RAM(Random Access Memory)などのコンピュータ42の1以上の記憶媒体として適切に実現される。データストレージコンポーネント24,28,32はまた、ネットワークにアクセス可能なPACS(Picture Archiving and Communications System)、外部のハードドライブ、光ディスクなどの他の記憶媒体により実現されてもよい。
また、各種処理モジュール22,26,30は、取得したB1マップを用いて(i)マルチチャネルRF送信フィールドをシム処理し、(ii)RF送信パワーを最適化する複数のRF送信チャネル20について最適化された振幅及び位相パラメータの計算を含むモジュール30により実行される処理を含む、ここに開示される処理を実行するため、コンピュータ42の図示されたプロセッサ又は他のプロセッサにより実行可能な命令を格納する記憶媒体により実現可能である。このような命令を格納する記憶媒体は、例えば、ハードディスクドライブ、他の磁気記憶媒体、光ディスク又は他の光記憶媒体、RAM、ROM、フラッシュメモリ、他の電子記憶媒体などであってもよい。
図2及び3を参照して、RFシムRF送信パワー最適化モジュール30により適切に実行される最適化された振幅及び位相パラメータの計算の例示的な具体例が説明される。ここに開示されるアプローチは、RF送信チャネルについて取得されたB1マップを用いてシム処理とRF送信パワー最適化との双方を実行する。これは、RF送信パワーを測定及び調整するために追加的なMRデータの取得を実行することを回避し、パワー最適化に用いられるRF送信パワーメトリックの選択に関してフレキシビリティを提供する。例えば、RF送信パワーメトリックは、関心領域における平均RF送信パワーとすることができ(例えば、心臓イメージングのケースにおいて心臓を包囲するなど)、又は関心スライスの平均RF送信パワーとすることができ、又は関心のある空間上のポイントにおけるRF送信パワーとすることが可能である。
図2及び3の例示的な具体例は、各RF送信チャネルについて(コンプレックス)B1を取得することによって開始される。このため、処理60において、マッピング対象のRF送信チャネルが選択される。処理62において、選択されたRF送信チャネルについて、振幅スケールが1.0に設定され、相対位相が0°に設定され、パワーレベルがPcalibとして示される計測パワーレベルに設定される。より一般には、処理62において、これらのパラメータは選択された計測又はリファレンスレベルに設定され、例えば、0°以外のリファレンス相対位相を利用することが想定される。処理64では、選択されたRF送信チャネル以外のすべてのRF送信チャネルについて、振幅スケールは0.0に設定され、パワーレベルはゼロに設定される。処理68において、B1マップは、振幅スケール=1.0、相対位相=0°、パワーレベル=Pcalibとなるパラメータを有する選択されたチャネルのみからの送信を用いて取得される。ループ又は繰り返し処理70では、複数のRF送信チャネル20の(コンプレックス)B1マップセット72を生成するため、複数のRF送信チャネル20の各RF送信チャネルを選択及びマッピングするため、処理60,62,64,68が繰り返される。
B1マッピング処理68に適したアプローチでは、関心スライス又はボリュームの2次元又は3次元B1マップ(好ましくは、ロードされたイメージング対象の内部又は同時)が取得される。B1マッピングは、所定のターゲットB1振幅(振幅スケール1.0など)及びRFパワー(パワーPcalibなど)のRFパルスを適切に利用してもよい。パワーレベルPcalibは、固定的であって、典型的には低いパワーレベルとすることが可能であり、任意的には従来のRFドライブスケール決定から導出される。B1マップは、コンプレックスB1値(すなわち、位相情報を含むB1値)をマッピングし、ターゲット又は名目的なB1値に対する実際のB1値又は相対的なB1値を表すものであってもよい。所与のRF送信チャネルのB1マップは、当該RF送信チャネルの実際の送信感度を表す。
図2を参照し続けて、複数のRF送信チャネル20のB1マップセット72が取得されると、計算処理80において、取得したB1マップ72を用いて(i)マルチチャネルRF送信フィールドをシム処理すると共に、(ii)RF送信パワーを最適化することによって、複数のRF送信チャネル20について最適化された振幅及び位相パラメータが計算される。
図3を参照して、計算処理80を実現する例示的な適切な処理が説明される。例示的なアプローチはまず、マルチチャネルRF送信フィールドの空間均一性を最適化するためシム処理を計算し、その後に所望のRF送信パワーメトリックを実現するため、シム処理されたRF送信チャネルの振幅を調整する。図3において実現されるシム処理は繰り返し的なものであり、初期的な振幅(又は振幅スケール)及び相対的位相が複数のRF送信チャネル20の各RF送信チャネルについて選択される処理82からスタートする。初期的な振幅及び位相は、B1送信フィールドの均一性を繰り返し向上させるため繰り返し調整されるべきであり、従って、最終的に最適化された値に近い初期値を有することが繰り返しの計算時間を低減するが、初期値は一般には重大ではない。いくつかの実施例では、振幅スケール=1.0及び相対位相=0°が、すべてのRF送信チャネルの初期値として使用される。あるいは、事前情報が利用可能である場合、処理82において、それは初期値を設定するため利用可能である。例えば、以前の同様の被検者(例えば、体重が類似している、体の大きさが類似しているなど)について決定された最適化された振幅及び位相が、初期値として利用されてもよい。処理84において、B1マップ72は、これら初期的な振幅及び位相値に基づき調整される。これは、コンプレックスB1値と初期的な振幅スケール値とを乗算し、B1位相を初期的な相対位相値だけシフトすることによってピクセル単位で実行可能である。このように調整されたB1マップは、その後に処理84において合成され、処理82において選択された初期的なパラメータにより作動する複数のRF送信チャネル20を用いてマルチチャネル送信モードにより取得されるB1マップを生成する。
処理82において選択された初期的なパラメータにより作動する複数のRF送信チャネル20を用いてマルチチャネル送信モードにより取得されるこのB1マップは、処理88において、空間均一性に関して解析される。処理88は、RF送信フィールドの均一性の指標を有する性能指数を適切に利用する。いくつかの実施例では、分散係数がRF送信フィールド均一性を測定する性能指数として利用される。しかしながら、他の均一性の性能指数が利用可能である。処理88が均一性が満たされていないことを検出した場合(例えば、計算された分散の性能指数が許容される最大分散閾値より大きい)、処理90において、振幅(又は振幅スケール)及び位相は、性能指数を向上させるため調整される。処理90は、各種振幅及び位相パラメータに関する分散の偏微分を計算し、勾配降下改良ステップを利用するものなど、何れか適切な繰り返し調整アルゴリズムを利用可能である。その後、処理は、調整処理90により調整される振幅及び位相パラメータにより作動する複数のRF送信チャネル20を用いてマルチチャネル送信モードにおいて取得される調整されたB1マップを生成するため処理84に戻り、新たな性能指数が、処理88において繰り返し調整されたパラメータが良好な空間均一性のマルチチャネル送信モードB1マップを生成していると判断されるまで、処理88における最大分散閾値又は他の良好な均一性基準と比較される処理86において計算される。この最終的なマップは、シム処理されたB1マップ92として適切にみなされる。
処理82,84,86,88,90により実現される繰り返しのシム処理は例示的な具体例であり、他のシム処理が利用されてもよい。一般に、最小分散係数について(又は他の均一性最適化基準によって測定される)個々のB1マップ72を合成するため最適な相対振幅及び位相パラメータを決定する何れかの適合方法が利用されてもよい。合成されたB1マップの均一性をテストしながら、逐次繰り返される位相及び振幅係数に関する力ずくのアプローチがまた想定される。
シム処理されたB1マップ92は、シム処理されたマルチチャネルRF励起の複数のRF送信チャネル20によって適用されるとイメージング対象の内部にあるシム処理されたB1フィールドを表す。シム処理82,84,86,88,90により最適化される振幅は、最適化された相対振幅である。これは、マルチチャネル送信モードにおいてB1送信フィールドの均一性を決定するのは、互いに関して最適化された振幅の値であるためである。従って、シム処理82,84,86,88,90により出力された最適化された相対振幅は、(一般に)何れか特定のRF送信パワーレベルを提供しない。しかしながら、シム処理されたB1マップ92の効果的な性質は、これらの値が所望のB1振幅(又は等価的には所望のRF送信パワーレベル)を実現するため、個々のチャネルパワー及び位相に直接関連しうることである。
従って、シム処理されたB1マップ92は、各チャネルのB1マップを取得するため用いられる既知のパワーレベルを、シム解析(処理82,84,86,88,90)から導出されたシム係数を利用した訂正に続に従って取得されたB1フィールド分布及び振幅に関連付けることによって、RFパワーレベル(すなわち、ドライブスケール)を導出するのに利用される。これは、最も均一的な励起を提供するため決定された位相及び振幅係数により各RFチャネルを導出する際にターゲットB1フィールドが正確に取得されることを保障する。このため、RF送信パワーメトリックは、例えば、(i)関心領域における平均RF送信パワー、(ii)関心スライスにおける平均RF送信パワー、(iii)関心空間ポイントにおけるRF送信パワーなどとすることができる。完全なシム処理されたB1マップ92が、処理94による処理のため利用可能であるため、関心イメージングタスクに適したRF送信パワーメトリックを選択するのに有意なフレキシビリティがある。例えば、イメージの中心において90°のフリップ角度を有することが重要である場合、RF送信パワーメトリックは、イメージングボリュームの中心におけるRF送信パワーとすることが可能である。スライスをイメージングするため、RF送信パワーメトリックの選択は、当該スライスにおける平均RF送信パワーであってもよい。
処理94により決定されるRF送信パワーメトリックは、処理96においてパワースケーリングファクタを決定するため、RF送信パワーメトリックの所望値と比較され、RF送信チャネルのシム処理された振幅が、RFシム処理と所望のRF送信パワーとの双方を実現するため、最適化された振幅及び位相98に到達するようパワースケーリングファクタによりスケーリングされる。例えば、処理94により決定されるRF送信パワーメトリックがB1meansとして示され(振幅単位により)、RF送信パワーメトリックの所望値がB1targetとして示される場合(再び振幅単位により)、スケーリングファクタは、B1target/B1meansとなる。その後、これらの振幅は、当該スケーリングファクタにより適切にスケーリングされる。この調整の実行において、RF送信パワーメトリックの選択は振幅単位であり、この振幅がスケーリングファクタ(B1target/B1means)によりスケーリングされることによって、対応するRF送信パワーがファクタ(B1target/B1means)2によりスケーリングされることになることに留意すべきである。RF送信パワーメトリックの選択は、振幅単位又はパワー単位とすることができる。パワー単位の具体例を利用して、処理94により決定されるRF送信パワーメトリックがP1meansとして示され(パワー単位により)、RF送信パワーメトリックの所望値がP1targetとして示される場合(パワー単位により)、振幅のスケーリングファクタは(P1target/P1means)1/2となり、対応するRF送信パワーは(P1target/P1means)によりスケーリングされる。
図3の実施例では、シム処理は、まず例示的な処理82,84,86,88,90により実行され、その後に処理94,96,98によって、取得したB1マップ72を用いてシム処理とRF送信パワー最適化との双方により実行される。
他の実施例では、シム処理及びRF送信パワー最適化は、取得したB1マップを再び用いて単一の処理により同時に実行可能である。例えば、このような実施例では、判定ブロック88において用いられる性能指数は、(i)RF送信フィールド均一性の指標(分散係数など)及び(ii)RF送信フィールドパワーの指標(関心スライス又は領域における平均B1フィールドなど)を合成する性能指数となるよう修正される。このような実施例では、例えば、性能指数は、(i)分散係数と、(ii)RF送信フィールドパワー(B1means)の指標とターゲットRF送信フィールドパワー(B1target)とを比較する項(B1target−B1means)2との加重和であってもよい。この修正された性能指数によると、繰り返しの処理82,84,86,88,90が、何れの側面が最適を支配するか選択する2つの項の間の重み付けによって、シム処理(分散係数項を最適化することによって)及びRF送信パワー((B1target−B1means)2の項を最適化することによって)を同時に実行可能である。本実施例では、処理94,96,98は、修正された性能指数が最適化処理82,84,86,88,90がRF送信パワーメトリックを最適化することを保障するため、適切に省略される。
図2のB1マッピングアプローチでは、各RF送信チャネルのB1マップは、当該チャネルのみをB1マッピングシーケンスにおいて作動することによって取得される。しかしながら、他のB1マッピングアプローチは、B1マップ72のセットを生成するため利用可能である。例えば、各B1マッピングの取得において、1つを除くすべてのチャネルが起動され、B1マッピング取得が複数回繰り返され(RF送信チャネル20の個数Nに等しい)、異なるチャネルは各回において起動されない、1つを除くすべてのマッピングアプローチが利用可能である。1つを除くすべてのアプローチでは、各チャネルの相対位相は、初期的には方形励起に対して固定され、以降のB1マップ取得は異なるチャネル振幅をゼロに設定してもよい。固定的な関係を用いて異なるRF送信チャネルグループが起動され、独立したRF送信チャネルと同じ数のB1マップが取得されるまで、B1マップが取得される毎に当該関係が置換される、このアプローチに対する変形がまた適している。B1マッピングデータをN個のチャネルのB1マップセット72に変換するため、物理チャネルが仮想チャネル(要素の組み合わせから構成された)にマッピングされる。このような1つを除くすべて又は他の組み合わせ的なマッピング処理は、B1マッピング処理のロウバスト性を向上させ、適合処理を促進することが可能である。
本出願は、1以上の好適な実施例を説明した。上述した詳細な説明を参照及び理解した者には改良及び変形が相当するであろう。本出願は、添付した請求項又はその均等の範囲内にある限り、このようなすべての改良及び変形を含むものとして解釈されることが意図される。