磁気共鳴装置およびプログラム
【課題】送信コイルのチャンネル数が増加しても、短時間でシミングの計算を行うことを提供する。
【解決手段】送信コイルに供給される信号のパターンを変更し、信号のパターンごとに送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求める。そして、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布Ｂ１２π／８と、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布Ｂ１１０π／８とを選択する。次に、送信磁場分布Ｂ１２π／８と送信磁場分布Ｂ１１０π／８との合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を求め、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを決定する。
【解決手段】送信コイルに供給される信号のパターンを変更し、信号のパターンごとに送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求める。そして、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布Ｂ１２π／８と、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布Ｂ１１０π／８とを選択する。次に、送信磁場分布Ｂ１２π／８と送信磁場分布Ｂ１１０π／８との合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を求め、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを決定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁場の不均一性を低減するための磁気共鳴装置およびプログラムに関する。
【背景技術】
【０００２】
磁気共鳴装置では、被検体の電気的特性によって、送信磁場分布が不均一になることが知られている。そこで、送信磁場分布を均一化するためのＲＦシミングが提案されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−０２９６４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＲＦシミングの一つとして、送信コイルの各給電点に供給する信号の位相および振幅を最適化する静的ＲＦシミング（static RF shimming）が知られている。この方法では、送信コイルのチャンネル数が少ない場合には、短時間で最適化を行うことができるが、送信コイルのチャンネル数が多くなるにつれて、最適化に必要な計算時間が増大するという問題がある。近年、送信コイルの多チャンネル化が進みつつあるので、送信コイルのチャンネル数が増加しても、短時間で最適化できることが望まれている。
【課題を解決するための手段】
【０００５】
本発明の第１の態様は、ｎチャンネルの送信コイルを有する磁気共鳴装置であって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成手段と、
前記送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択手段と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出手段と、
を有する磁気共鳴装置である。
【０００６】
本発明の第２の態様は、ｎチャンネルの送信コイルを有する磁気共鳴装置のプログラムであって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成処理と、
前記送信磁場分布作成処理により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択処理と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出処理と、
を計算機に実行させるためのプログラムである。
【発明の効果】
【０００７】
送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択し、前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出している。したがって、信号を算出するために必要な変数の数を少なくすることができ、前記送信コイルの各チャンネルに供給する信号を算出するために必要な計算時間を短縮することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】送信コイル２４の構造の説明図である。
【図３】ＲＦシミングの手順を示すフローである。
【図４】ステップＳＴ２で実行されるスキャンを示す図である。
【図５】各スキャンを実行するときに送信コイル２４に供給される信号のパターンを示す図である。
【図６】スキャンＳＣ１を実行したときの送信磁場分布Ｂ１Δθの求め方の説明図である。
【図７】スキャンＳＣ５を実行したときの送信磁場分布Ｂ１１０π／８を示す図である。
【図８】スキャンＳＣ１〜ＳＣ８ごとに求められた送信磁場分布を示す図である。
【図９】相関を示す図である。
【図１０】チャンネルＣＨ１の信号の算出方法の説明図である。
【図１１】第２の形態におけるＲＦシミングの手順を示すフローである。
【図１２】ステップＳＴ１１で実行されるスキャンを示す図である。
【図１３】各スキャンを実行するときに送信コイル２４の各チャンネルに供給される信号を示す図である。
【発明を実施するための形態】
【０００９】
以下、発明を実施するための形態を説明するが、本発明は、以下の形態に限定されることはない。
【００１０】
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１１】
マグネット２は、被検体１２が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。
【００１２】
図２は、送信コイル２４の構造の説明図である。
送信コイル２４は、ｎチャンネルのコイルである。第１の形態では、ｎ＝８、つまり、送信コイル２４は８チャンネルＣＨ１〜ＣＨ８のコイルである。尚、ｎの値は、８に限定されることはなく、８より小さくてもよく、大きくてもよい。各チャンネルＣＨ１〜ＣＨ８は、例えば、マイクロストリップコイルである。各チャンネルＣＨ１〜ＣＨ８には、後述する送信器６から、ＲＦパルスを送信するための信号が供給される。
【００１３】
図１に戻って説明を続ける。
テーブル３は、被検体１２を支持するためのクレードル３ａを有している。クレードル３ａが移動することによって、被検体１２がボア２１に搬入される。
受信コイル４は被検体１２の頭部に取り付けられている。
【００１４】
ＭＲ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、操作部１０、および表示部１１を有している。
【００１５】
シーケンサ５は、中央処理装置９の制御を受けて、パルスシーケンスの情報を送信器６および勾配磁場電源７に送る。
【００１６】
送信器６は、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に信号を供給する。
勾配磁場電源７は勾配コイル２３に信号を供給する。
受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、中央処理装置９に出力する。
【００１７】
中央処理装置９は、シーケンサ５および表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、送信磁場分布作成手段９１、送信磁場分布選択手段９２、算出手段９３などを有している。
【００１８】
送信磁場分布作成手段９１が、スキャンによって取得されたデータに基づいて、送信磁場分布Ｂ１Δθ（図８参照）を作成する。
【００１９】
送信磁場分布選択手段９２は、送信磁場分布作成手段９１により作成された送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い第１の送信磁場分布と、第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する。
【００２０】
算出手段９３は、第１の送信磁場分布と第２の送信磁場分布とに基づいて、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に供給する信号を算出する。
【００２１】
中央処理装置９は、送信磁場分布作成手段９１、送信磁場分布選択手段９２、算出手段９３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２２】
操作部１０は、オペレータにより操作され、種々の情報を中央処理装置９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２３】
次に、第１の形態において、送信コイルの送信磁場分布を均一にするためのＲＦシミングを行う手順について説明する。
【００２４】
図３は、ＲＦシミングの手順を示すフローである。
ステップＳＴ１では、図１に示すように、被検体１２をボア２１に搬入する。被検体１２をボア２１に搬入したら、ステップＳＴ２に進む。
【００２５】
ステップＳＴ２では、送信コイル２４の送信磁場分布Ｂ１Δθを求める。第１の形態では、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に供給する信号のパターンを変更しながらスキャンを行い、信号のパターンごとに、送信コイル２４の送信磁場分布Ｂ１Δθを求める。以下に、ステップＳＴ２で実行されるスキャンについて説明する。
【００２６】
図４は、ステップＳＴ２で実行されるスキャンを示す図、図５は、各スキャンを実行するときに送信コイル２４に供給される信号のパターンを示す図である。
【００２７】
ステップＳＴ２では、図４に示すように、スキャンＳＣ１〜ＳＣ８が実行される。各スキャンＳＣ１〜ＳＣ８を実行する場合、各チャンネルＣＨ１〜ＣＨ８には、図５に示すような信号が供給される。具体的には、以下の通りである。
【００２８】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：ＩΔθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋Δθ）｝
チャンネルＣＨ３：Ｉ２Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋２・Δθ）｝
チャンネルＣＨ４：Ｉ３Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋３・Δθ）｝
チャンネルＣＨ５：Ｉ４Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋４・Δθ）｝
チャンネルＣＨ６：Ｉ５Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋５・Δθ）｝
チャンネルＣＨ７：Ｉ６Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋６・Δθ）｝
チャンネルＣＨ８：Ｉ７Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋７・Δθ）｝
ここで、Ａ０：振幅
θ０＋ｋ・Δθ（ｋ＝０〜７）：位相
ただし、Δθ＝２π・（Ｍ／８） （Ｍ＝１〜８）
【００２９】
例えば、スキャンＳＣ１を実行する場合、図５に示すように、Ｍ＝１、つまり、Δθ＝２π・（１／８）＝２π／８となる。したがって、スキャンＳＣ１を実行する場合、各チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθは、Δθに２π／８を代入して、以下のようになる。
【００３０】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：Ｉ２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２π／８）｝
チャンネルＣＨ３：Ｉ４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋４π／８）｝
チャンネルＣＨ４：Ｉ６π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋６π／８）｝
チャンネルＣＨ５：Ｉ８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋８π／８）｝
チャンネルＣＨ６：Ｉ１０π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１０π／８）｝
チャンネルＣＨ７：Ｉ１２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１２π／８）｝
チャンネルＣＨ８：Ｉ１４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１４π／８）｝
【００３１】
また、スキャンＳＣ２を実行する場合、Ｍ＝２、つまり、Δθ＝２π・（２／８）＝４π／８となる。したがって、スキャンＳＣ２を実行する場合、各チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθは、Δθに４π／８を代入して、以下のようになる。
【００３２】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：Ｉ４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋４π／８）｝
チャンネルＣＨ３：Ｉ８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋８π／８）｝
チャンネルＣＨ４：Ｉ１２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１２π／８）｝
チャンネルＣＨ５：Ｉ１６π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１６π／８）｝
チャンネルＣＨ６：Ｉ２０π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２０π／８）｝
チャンネルＣＨ７：Ｉ２４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２４π／８）｝
チャンネルＣＨ８：Ｉ２８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２８π／８）｝
【００３３】
以下同様に、スキャンＳＣ３〜ＳＣ８を実行するときのΔθは、図５に示すように、それぞれ６π／８〜１６π／８となる。したがって、Δθの値を変更することによって、チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθのパターンを変更することができる。このように、チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθのパターンを変更しながら、スキャンＳＣ１〜ＳＣ８を実行し、送信磁場分布Ｂ１Δθを求める。
【００３４】
図６は、スキャンＳＣ１を実行したときの送信磁場分布Ｂ１Δθの求め方の説明図である。
スキャンＳＣ１では、チャンネルＣＨ１〜ＣＨ８に、Δθ＝２π／８の信号Ｉ０〜Ｉ１４π／８が供給される。スキャンＳＣ１によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１（図１参照）が、スキャンＳＣ１によって取得されたデータに基づいて、チャンネルＣＨ１〜ＣＨ８に信号Ｉ０〜Ｉ１４π／８が供給されたときの送信磁場分布Ｂ１Δθ（＝Ｂ１２π／８）を作成する。したがって、スキャンＳＣ１を実行することによって、送信磁場分布Ｂ１２π／８を作成することができる。
【００３５】
以下同様に、スキャンＳＣ２〜ＳＣ８を実行し、送信磁場分布を求める。例えば、図７には、スキャンＳＣ５を実行したときの送信磁場分布Ｂ１１０π／８が示されている。スキャンＳＣ５では、チャンネルＣＨ１〜ＣＨ８に、Δθ＝１０π／８の信号Ｉ０〜Ｉ７０π／８が供給される。スキャンＳＣ５によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１が、スキャンＳＣ５によって取得されたデータに基づいて、チャンネルＣＨ１〜ＣＨ８に信号Ｉ０〜Ｉ７０π／８が供給されたときの送信磁場分布Ｂ１Δθ（＝Ｂ１１０π／８）を作成する。したがって、スキャンＳＣ５を実行することによって、送信磁場分布Ｂ１１０π／８を作成することができる。
【００３６】
図８に、スキャンＳＣ１〜ＳＣ８ごとに求められた送信磁場分布Ｂ１Δθを、符号「Ｂ１２π／８」〜「Ｂ１１６π／８」で示す。送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求めたら、ステップＳＴ３に進む。
【００３７】
ステップＳＴ３では、送信磁場分布選択手段９２（図１参照）が、ステップＳＴ２で求められた送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布を選択する。磁場の均一性が最も高い送信磁場分布を求める方法の一例としては、送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の各々について、ピクセル間の信号強度のばらつきを算出し、算出した信号強度のばらつきを比較する方法がある。第１の形態では、磁場の均一性が最も高い送信磁場分布として、送信磁場分布Ｂ１２π／８が選択されたとする。送信磁場分布Ｂ１２π／８を選択したら、ステップＳＴ４に進む。
【００３８】
ステップＳＴ４では、送信磁場分布選択手段９２が、送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布を選択する。以下に、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布を選択する方法について説明する。
【００３９】
送信磁場分布選択手段９２は、先ず、送信磁場分布Ｂ１２π／８と、その他の送信磁場分布Ｂ１４π／８〜Ｂ１１６π／８の各々との相関を計算する。相関は、一般的な統計的手法によって算出することができる。図９に、算出した相関を概略的に示す。図９では、算出した相関を記号「Ｃ１」〜「Ｃ７」で示してある。送信磁場分布選択手段９２は、算出したＣ１〜Ｃ７の値を比較し、相関Ｃ１〜Ｃ７の中から、最小値となる相関を特定する。ここでは、相関Ｃ４が最小値になったとする。したがって、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布は、送信磁場分布Ｂ１１０π／８となる。送信磁場分布Ｂ１１０π／８を選択したら、ステップＳＴ５に進む。
【００４０】
ステップＳＴ５では、算出手段９３（図１参照）が、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８と、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８とに基づいて、均一性の高い送信磁場分布を実現するための各チャンネルの信号を算出する。ステップＳＴ５は、ステップＳＴ５１〜ステップＳＴ５３を有している。以下に、各ステップＳＴ５１〜ＳＴ５３について説明する。
【００４１】
ステップＳＴ５１では、算出手段９３は、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８と、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８とを合成し、合成送信磁場分布を作成する。合成送信磁場分布は、以下の式で表される。
Ｂ１ｃｏｍｂ（ａｍｐ，θ）
＝Ｂ１２π／８＋ａｍｐ・ｅｘｐ（ｉθ）・Ｂ１１０π／８ ・・・（１）
ここで、Ｂ１ｃｏｍｂ（ａｍｐ，θ）：合成送信磁場分布
ａｍｐ：振幅
θ：位相
尚、振幅ａｍｐおよび位相θは、変数である。
【００４２】
合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を作成したら、ステップＳＴ５２に進む。
ステップＳＴ５２では、算出手段９３は、式（１）に基づいて、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを求める。この組合せは、例えば最適化手法によって求めることができる。ここでは、ａｍｐ＝ａｍｐｕｎｉ、θ＝θｕｎｉと算出されたとする。振幅ａｍｐおよび位相θの組合せを算出したら、ステップＳＴ５３に進む。
【００４３】
ステップＳＴ５３では、算出手段９３は、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するために各チャンネルＣＨ１〜ＣＨ８に供給する信号Ｊ１〜Ｊ８を算出する。以下に、信号Ｊ１〜Ｊ８の求め方について説明する。
【００４４】
図１０は、チャンネルＣＨ１の信号の算出方法の説明図である。
図１０（ａ）は、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８（図６参照）が得られるときの各チャンネルＣＨ１〜ＣＨ８の信号Ｉ０〜Ｉ１４π／８を示す。また、図１０（ｂ）は、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８（図７参照）が得られるときの各チャンネルＣＨ１〜ＣＨ８の信号Ｉ０〜Ｉ７０π／８を示す。
【００４５】
先ず、図１０（ｂ）に示されている信号Ｉ０〜Ｉ７０π／８に対して、ステップＳＴ５２で算出された振幅ａｍｐｕｎｉおよび位相θｕｎｉを付加する。具体的には、信号Ｉ０〜Ｉ７０π／８の振幅Ａ０にａｍｐｕｎｉを乗算し、信号Ｉ０〜Ｉ７０π／８の位相にθｕｎｉを加算する。図１０（ｃ）に、信号Ｉ０〜Ｉ７０π／８に振幅ａｍｐｕｎｉおよび位相θｕｎｉを付加することにより得られた信号を、符号「Ｉ０′」〜「Ｉ７０π／８′」で示す。信号Ｉ０′〜Ｉ７０π／８′は、具体的には、以下の式で表される。
【００４６】
Ｉ０′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋θｕｎｉ）｝
Ｉ１０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋１０π／８＋θｕｎｉ）｝
Ｉ２０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋２０π／８＋θｕｎｉ）｝
Ｉ３０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋３０π／８＋θｕｎｉ）｝
Ｉ４０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋４０π／８＋θｕｎｉ）｝
Ｉ５０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋５０π／８＋θｕｎｉ）｝
Ｉ６０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋６０π／８＋θｕｎｉ）｝
Ｉ７０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋７０π／８＋θｕｎｉ）｝
【００４７】
次に、図１０（ａ）の信号Ｉ０〜Ｉ１４π／８と、図１０（ｃ）の信号Ｉ０′〜Ｉ７０π／８′とを、各チャンネルごとに加算する。図１０（ｄ）に、加算により得られた各チャンネルの信号ＪＩ〜Ｊ８を示す。このようにして、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための各チャンネルの信号ＪＩ〜Ｊ８を算出することができる。信号Ｊ１〜Ｊ８を算出したら、ＲＦシミングのフローを終了する。
【００４８】
第１の形態では、送信コイル２４に供給される信号のパターンを変更し、信号のパターンごとに送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求めている（ステップＳＴ２）。そして、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布Ｂ１２π／８を選択している（ステップＳＴ３）。ただし、被検体１２の電気的特性などが原因で、送信磁場分布Ｂ１２π／８の中にも、感度の低い部分と、感度の高い部分が現れるので、このような感度のばらつきは、できるだけ小さくすることが望ましい。そこで、第１の形態では、送信磁場分布Ｂ１２π／８だけでなく、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布Ｂ１１０π／８も求め（ステップＳＴ４）、送信磁場分布Ｂ１２π／８と送信磁場分布Ｂ１１０π／８との合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を求めている（ステップＳＴ５１）。送信磁場分布Ｂ１１０π／８は、送信磁場分布Ｂ１２π／８との相関が最も弱いので、送信磁場分布Ｂ１１０π／８を用いることによって、送信磁場分布Ｂ１２π／８の中の感度の低い部分を、感度が高くなるように補正することができ、更に、送信磁場分布Ｂ１２π／８の中の感度の高い部分を、感度が低くなるように補正することができる。したがって、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを決定することによって、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための信号ＪＩ〜Ｊ８を算出することができる。
【００４９】
また、第１の形態では、ステップＳＴ２で求めた８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、２つの送信磁場分布Ｂ１２π／８およびＢ１１０π／８だけを選択して、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）の式を決定している。したがって、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を全て用いなくても、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式を決定することができるので、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数の数を少なくすることができる。第１の形態では、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数は２つ（振幅ａｍｐおよび位相θ）だけで済む。したがって、送信コイル２４に供給する信号ＪＩ〜Ｊ８を決定するために必要な計算時間を短縮することができる。
【００５０】
尚、第１の形態では、ステップＳＴ２において、チャンネルＣＨ１〜ＣＨ８に供給される信号の位相を変化させることによって、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させている。しかし、振幅Ａ０を変化させることによって、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させてもよい。また、振幅と位相の両方を変化させて、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させてもよい。
【００５１】
また、第１の形態では、マグネット２に内蔵された送信コイル２４のシミングを行う場合について説明されている。しかし、本発明は、マグネット２とは別に備えられた送信コイル（例えば、人体の特定の部位の撮影に特化した専用コイル）のシミングを行う場合にも適用することができる。
【００５２】
（２）第２の形態
第１の形態は、送信磁場分布Ｂ１Δθ（図８参照）を実測により求めている。しかし、送信磁場分布Ｂ１Δθを計算により求めてもよい。第２の形態では、送信磁場分布Ｂ１Δθを計算により求めてＲＦシミングを行う例について説明する。
【００５３】
図１１は、第２の形態におけるＲＦシミングの手順を示すフローである。
ステップＳＴ１では、図１に示すように、被検体１２をボア２１に搬入する。被検体１２をボア２１に搬入したら、ステップＳＴ１１に進む。
【００５４】
ステップＳＴ１１では、送信コイル２４の各チャンネルごとに送信磁場分布を求める。
図１２は、ステップＳＴ１１で実行されるスキャンを示す図、図１３は、各スキャンを実行するときに送信コイル２４の各チャンネルに供給される信号を示す図である。
【００５５】
ステップＳＴ１１では、図１３に示すように、スキャンＳＣ１′〜ＳＣ８′が実行される。各スキャンＳＣ１′〜ＳＣ８′を実行する場合、各チャンネルＣＨ１〜ＣＨ８には、図１３に示すような信号が供給される。
【００５６】
例えば、スキャンＳＣ１′を実行する場合、図１３（ａ）に示すように、チャンネルＣＨ１にのみ信号Ｉ０が供給される。スキャンＳＣ１′によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１（図１参照）が、スキャンＳＣ１′によって取得されたデータに基づいて、チャンネルＣＨ１の送信磁場分布Ｂ１１を作成する。したがって、スキャンＳＣ１′を実行することによって、チャンネルＣＨ１による送信磁場分布Ｂ１１を作成することができる。
【００５７】
また、スキャンＳＣ２′を実行する場合、図１３（ｂ）に示すように、チャンネルＣＨ２にのみ信号Ｉ０が供給される。スキャンＳＣ２′によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１が、スキャンＳＣ２′によって取得されたデータに基づいて、チャンネルＣＨ２の送信磁場分布Ｂ１２を作成する。したがって、スキャンＳＣ２′を実行することによって、チャンネルＣＨ２による送信磁場分布Ｂ１２を作成することができる。
【００５８】
以下同様に、スキャンＳＣ３′〜ＳＣ８′を実行する場合、図１４（ｃ）〜（ｈ）に示すように、各チャンネルＣＨ３〜ＣＨ８にのみ信号Ｉ０を供給することにより、チャンネルＣＨ３〜ＣＨ８の送信磁場分布Ｂ１３〜Ｂ１８を求めることができる。したがって、スキャンＳＣ１′〜ＳＣ８′を実行することによって、チャンネルＣＨ１〜ＣＨ８ごとに、送信磁場分布Ｂ１１〜Ｂ１８を求めることができる。送信磁場分布１１〜Ｂ１８を求めた後、ステップＳＴ２に進む。
【００５９】
ステップＳＴ２では、送信磁場分布作成手段９１が、ステップＳＴ１１で求めたチャンネルＣＨ１〜ＣＨ８ごとの送信磁場分布Ｂ１１〜Ｂ１８に基づいて、送信コイル２４に信号Ｉ０〜Ｉ７Δθを供給したときに得られる送信磁場分布Ｂ１Δθ（図８参照）を算出する。送信磁場分布Ｂ１Δθは、例えば以下の式を用いて算出することができる。
【００６０】
【数１】
送信磁場分布Ｂ１ｎ（ｎ＝１〜８）はステップＳＴ１１で求められているので、求めた送信磁場分布Ｂ１ｎを上記の式に代入することによって、送信磁場分布Ｂ１Δθ（Ｂ１２π／８〜Ｂ１１６π／８）を算出することができる。送信磁場分布Ｂ１Δθを算出したら、ステップＳＴ３に進む。
ステップＳＴ３〜ＳＴ５は、第１の形態と同じであるので、説明は省略する。
【００６１】
第２の形態のように、送信磁場分布Ｂ１Δθを計算により求めても、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための信号ＪＩ〜Ｊ８を算出することができる。
【００６２】
また、計算により８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を算出した後は、第１の形態と同様に、２つの送信磁場分布Ｂ１２π／８およびＢ１１０π／８だけを選択して、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式を決定している。したがって、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数は２つ（振幅ａｍｐおよび位相θ）だけで済み、送信コイル２４に供給する信号ＪＩ〜Ｊ８を決定するために必要な計算時間を短縮することができる。
【符号の説明】
【００６３】
２ マグネット
３ テーブル
３ａ クレードル
５ シーケンサ
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 中央処理装置
１０ 操作部
１１ 表示部
１２ 被検体
９１ 送信磁場分布作成手段
９２ 送信磁場分布選択手段
９３ 算出手段
１００ ＭＲ装置
【技術分野】
【０００１】
本発明は、磁場の不均一性を低減するための磁気共鳴装置およびプログラムに関する。
【背景技術】
【０００２】
磁気共鳴装置では、被検体の電気的特性によって、送信磁場分布が不均一になることが知られている。そこで、送信磁場分布を均一化するためのＲＦシミングが提案されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−０２９６４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＲＦシミングの一つとして、送信コイルの各給電点に供給する信号の位相および振幅を最適化する静的ＲＦシミング（static RF shimming）が知られている。この方法では、送信コイルのチャンネル数が少ない場合には、短時間で最適化を行うことができるが、送信コイルのチャンネル数が多くなるにつれて、最適化に必要な計算時間が増大するという問題がある。近年、送信コイルの多チャンネル化が進みつつあるので、送信コイルのチャンネル数が増加しても、短時間で最適化できることが望まれている。
【課題を解決するための手段】
【０００５】
本発明の第１の態様は、ｎチャンネルの送信コイルを有する磁気共鳴装置であって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成手段と、
前記送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択手段と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出手段と、
を有する磁気共鳴装置である。
【０００６】
本発明の第２の態様は、ｎチャンネルの送信コイルを有する磁気共鳴装置のプログラムであって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成処理と、
前記送信磁場分布作成処理により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択処理と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出処理と、
を計算機に実行させるためのプログラムである。
【発明の効果】
【０００７】
送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択し、前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出している。したがって、信号を算出するために必要な変数の数を少なくすることができ、前記送信コイルの各チャンネルに供給する信号を算出するために必要な計算時間を短縮することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】送信コイル２４の構造の説明図である。
【図３】ＲＦシミングの手順を示すフローである。
【図４】ステップＳＴ２で実行されるスキャンを示す図である。
【図５】各スキャンを実行するときに送信コイル２４に供給される信号のパターンを示す図である。
【図６】スキャンＳＣ１を実行したときの送信磁場分布Ｂ１Δθの求め方の説明図である。
【図７】スキャンＳＣ５を実行したときの送信磁場分布Ｂ１１０π／８を示す図である。
【図８】スキャンＳＣ１〜ＳＣ８ごとに求められた送信磁場分布を示す図である。
【図９】相関を示す図である。
【図１０】チャンネルＣＨ１の信号の算出方法の説明図である。
【図１１】第２の形態におけるＲＦシミングの手順を示すフローである。
【図１２】ステップＳＴ１１で実行されるスキャンを示す図である。
【図１３】各スキャンを実行するときに送信コイル２４の各チャンネルに供給される信号を示す図である。
【発明を実施するための形態】
【０００９】
以下、発明を実施するための形態を説明するが、本発明は、以下の形態に限定されることはない。
【００１０】
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１１】
マグネット２は、被検体１２が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。
【００１２】
図２は、送信コイル２４の構造の説明図である。
送信コイル２４は、ｎチャンネルのコイルである。第１の形態では、ｎ＝８、つまり、送信コイル２４は８チャンネルＣＨ１〜ＣＨ８のコイルである。尚、ｎの値は、８に限定されることはなく、８より小さくてもよく、大きくてもよい。各チャンネルＣＨ１〜ＣＨ８は、例えば、マイクロストリップコイルである。各チャンネルＣＨ１〜ＣＨ８には、後述する送信器６から、ＲＦパルスを送信するための信号が供給される。
【００１３】
図１に戻って説明を続ける。
テーブル３は、被検体１２を支持するためのクレードル３ａを有している。クレードル３ａが移動することによって、被検体１２がボア２１に搬入される。
受信コイル４は被検体１２の頭部に取り付けられている。
【００１４】
ＭＲ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、操作部１０、および表示部１１を有している。
【００１５】
シーケンサ５は、中央処理装置９の制御を受けて、パルスシーケンスの情報を送信器６および勾配磁場電源７に送る。
【００１６】
送信器６は、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に信号を供給する。
勾配磁場電源７は勾配コイル２３に信号を供給する。
受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、中央処理装置９に出力する。
【００１７】
中央処理装置９は、シーケンサ５および表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、送信磁場分布作成手段９１、送信磁場分布選択手段９２、算出手段９３などを有している。
【００１８】
送信磁場分布作成手段９１が、スキャンによって取得されたデータに基づいて、送信磁場分布Ｂ１Δθ（図８参照）を作成する。
【００１９】
送信磁場分布選択手段９２は、送信磁場分布作成手段９１により作成された送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い第１の送信磁場分布と、第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する。
【００２０】
算出手段９３は、第１の送信磁場分布と第２の送信磁場分布とに基づいて、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に供給する信号を算出する。
【００２１】
中央処理装置９は、送信磁場分布作成手段９１、送信磁場分布選択手段９２、算出手段９３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２２】
操作部１０は、オペレータにより操作され、種々の情報を中央処理装置９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２３】
次に、第１の形態において、送信コイルの送信磁場分布を均一にするためのＲＦシミングを行う手順について説明する。
【００２４】
図３は、ＲＦシミングの手順を示すフローである。
ステップＳＴ１では、図１に示すように、被検体１２をボア２１に搬入する。被検体１２をボア２１に搬入したら、ステップＳＴ２に進む。
【００２５】
ステップＳＴ２では、送信コイル２４の送信磁場分布Ｂ１Δθを求める。第１の形態では、送信コイル２４の各チャンネルＣＨ１〜ＣＨ８に供給する信号のパターンを変更しながらスキャンを行い、信号のパターンごとに、送信コイル２４の送信磁場分布Ｂ１Δθを求める。以下に、ステップＳＴ２で実行されるスキャンについて説明する。
【００２６】
図４は、ステップＳＴ２で実行されるスキャンを示す図、図５は、各スキャンを実行するときに送信コイル２４に供給される信号のパターンを示す図である。
【００２７】
ステップＳＴ２では、図４に示すように、スキャンＳＣ１〜ＳＣ８が実行される。各スキャンＳＣ１〜ＳＣ８を実行する場合、各チャンネルＣＨ１〜ＣＨ８には、図５に示すような信号が供給される。具体的には、以下の通りである。
【００２８】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：ＩΔθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋Δθ）｝
チャンネルＣＨ３：Ｉ２Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋２・Δθ）｝
チャンネルＣＨ４：Ｉ３Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋３・Δθ）｝
チャンネルＣＨ５：Ｉ４Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋４・Δθ）｝
チャンネルＣＨ６：Ｉ５Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋５・Δθ）｝
チャンネルＣＨ７：Ｉ６Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋６・Δθ）｝
チャンネルＣＨ８：Ｉ７Δθ＝Ａ０・ｅｘｐ｛ｉ（θ０＋７・Δθ）｝
ここで、Ａ０：振幅
θ０＋ｋ・Δθ（ｋ＝０〜７）：位相
ただし、Δθ＝２π・（Ｍ／８） （Ｍ＝１〜８）
【００２９】
例えば、スキャンＳＣ１を実行する場合、図５に示すように、Ｍ＝１、つまり、Δθ＝２π・（１／８）＝２π／８となる。したがって、スキャンＳＣ１を実行する場合、各チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθは、Δθに２π／８を代入して、以下のようになる。
【００３０】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：Ｉ２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２π／８）｝
チャンネルＣＨ３：Ｉ４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋４π／８）｝
チャンネルＣＨ４：Ｉ６π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋６π／８）｝
チャンネルＣＨ５：Ｉ８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋８π／８）｝
チャンネルＣＨ６：Ｉ１０π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１０π／８）｝
チャンネルＣＨ７：Ｉ１２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１２π／８）｝
チャンネルＣＨ８：Ｉ１４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１４π／８）｝
【００３１】
また、スキャンＳＣ２を実行する場合、Ｍ＝２、つまり、Δθ＝２π・（２／８）＝４π／８となる。したがって、スキャンＳＣ２を実行する場合、各チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθは、Δθに４π／８を代入して、以下のようになる。
【００３２】
チャンネルＣＨ１：Ｉ０＝Ａ０・ｅｘｐ（ｉθ０）
チャンネルＣＨ２：Ｉ４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋４π／８）｝
チャンネルＣＨ３：Ｉ８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋８π／８）｝
チャンネルＣＨ４：Ｉ１２π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１２π／８）｝
チャンネルＣＨ５：Ｉ１６π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋１６π／８）｝
チャンネルＣＨ６：Ｉ２０π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２０π／８）｝
チャンネルＣＨ７：Ｉ２４π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２４π／８）｝
チャンネルＣＨ８：Ｉ２８π／８＝Ａ０・ｅｘｐ｛ｉ（θ０＋２８π／８）｝
【００３３】
以下同様に、スキャンＳＣ３〜ＳＣ８を実行するときのΔθは、図５に示すように、それぞれ６π／８〜１６π／８となる。したがって、Δθの値を変更することによって、チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθのパターンを変更することができる。このように、チャンネルＣＨ１〜ＣＨ８に供給される信号Ｉ０〜Ｉ７Δθのパターンを変更しながら、スキャンＳＣ１〜ＳＣ８を実行し、送信磁場分布Ｂ１Δθを求める。
【００３４】
図６は、スキャンＳＣ１を実行したときの送信磁場分布Ｂ１Δθの求め方の説明図である。
スキャンＳＣ１では、チャンネルＣＨ１〜ＣＨ８に、Δθ＝２π／８の信号Ｉ０〜Ｉ１４π／８が供給される。スキャンＳＣ１によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１（図１参照）が、スキャンＳＣ１によって取得されたデータに基づいて、チャンネルＣＨ１〜ＣＨ８に信号Ｉ０〜Ｉ１４π／８が供給されたときの送信磁場分布Ｂ１Δθ（＝Ｂ１２π／８）を作成する。したがって、スキャンＳＣ１を実行することによって、送信磁場分布Ｂ１２π／８を作成することができる。
【００３５】
以下同様に、スキャンＳＣ２〜ＳＣ８を実行し、送信磁場分布を求める。例えば、図７には、スキャンＳＣ５を実行したときの送信磁場分布Ｂ１１０π／８が示されている。スキャンＳＣ５では、チャンネルＣＨ１〜ＣＨ８に、Δθ＝１０π／８の信号Ｉ０〜Ｉ７０π／８が供給される。スキャンＳＣ５によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１が、スキャンＳＣ５によって取得されたデータに基づいて、チャンネルＣＨ１〜ＣＨ８に信号Ｉ０〜Ｉ７０π／８が供給されたときの送信磁場分布Ｂ１Δθ（＝Ｂ１１０π／８）を作成する。したがって、スキャンＳＣ５を実行することによって、送信磁場分布Ｂ１１０π／８を作成することができる。
【００３６】
図８に、スキャンＳＣ１〜ＳＣ８ごとに求められた送信磁場分布Ｂ１Δθを、符号「Ｂ１２π／８」〜「Ｂ１１６π／８」で示す。送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求めたら、ステップＳＴ３に進む。
【００３７】
ステップＳＴ３では、送信磁場分布選択手段９２（図１参照）が、ステップＳＴ２で求められた送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布を選択する。磁場の均一性が最も高い送信磁場分布を求める方法の一例としては、送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の各々について、ピクセル間の信号強度のばらつきを算出し、算出した信号強度のばらつきを比較する方法がある。第１の形態では、磁場の均一性が最も高い送信磁場分布として、送信磁場分布Ｂ１２π／８が選択されたとする。送信磁場分布Ｂ１２π／８を選択したら、ステップＳＴ４に進む。
【００３８】
ステップＳＴ４では、送信磁場分布選択手段９２が、送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布を選択する。以下に、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布を選択する方法について説明する。
【００３９】
送信磁場分布選択手段９２は、先ず、送信磁場分布Ｂ１２π／８と、その他の送信磁場分布Ｂ１４π／８〜Ｂ１１６π／８の各々との相関を計算する。相関は、一般的な統計的手法によって算出することができる。図９に、算出した相関を概略的に示す。図９では、算出した相関を記号「Ｃ１」〜「Ｃ７」で示してある。送信磁場分布選択手段９２は、算出したＣ１〜Ｃ７の値を比較し、相関Ｃ１〜Ｃ７の中から、最小値となる相関を特定する。ここでは、相関Ｃ４が最小値になったとする。したがって、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布は、送信磁場分布Ｂ１１０π／８となる。送信磁場分布Ｂ１１０π／８を選択したら、ステップＳＴ５に進む。
【００４０】
ステップＳＴ５では、算出手段９３（図１参照）が、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８と、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８とに基づいて、均一性の高い送信磁場分布を実現するための各チャンネルの信号を算出する。ステップＳＴ５は、ステップＳＴ５１〜ステップＳＴ５３を有している。以下に、各ステップＳＴ５１〜ＳＴ５３について説明する。
【００４１】
ステップＳＴ５１では、算出手段９３は、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８と、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８とを合成し、合成送信磁場分布を作成する。合成送信磁場分布は、以下の式で表される。
Ｂ１ｃｏｍｂ（ａｍｐ，θ）
＝Ｂ１２π／８＋ａｍｐ・ｅｘｐ（ｉθ）・Ｂ１１０π／８ ・・・（１）
ここで、Ｂ１ｃｏｍｂ（ａｍｐ，θ）：合成送信磁場分布
ａｍｐ：振幅
θ：位相
尚、振幅ａｍｐおよび位相θは、変数である。
【００４２】
合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を作成したら、ステップＳＴ５２に進む。
ステップＳＴ５２では、算出手段９３は、式（１）に基づいて、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを求める。この組合せは、例えば最適化手法によって求めることができる。ここでは、ａｍｐ＝ａｍｐｕｎｉ、θ＝θｕｎｉと算出されたとする。振幅ａｍｐおよび位相θの組合せを算出したら、ステップＳＴ５３に進む。
【００４３】
ステップＳＴ５３では、算出手段９３は、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するために各チャンネルＣＨ１〜ＣＨ８に供給する信号Ｊ１〜Ｊ８を算出する。以下に、信号Ｊ１〜Ｊ８の求め方について説明する。
【００４４】
図１０は、チャンネルＣＨ１の信号の算出方法の説明図である。
図１０（ａ）は、ステップＳＴ３で選択された送信磁場分布Ｂ１２π／８（図６参照）が得られるときの各チャンネルＣＨ１〜ＣＨ８の信号Ｉ０〜Ｉ１４π／８を示す。また、図１０（ｂ）は、ステップＳＴ４で選択された送信磁場分布Ｂ１１０π／８（図７参照）が得られるときの各チャンネルＣＨ１〜ＣＨ８の信号Ｉ０〜Ｉ７０π／８を示す。
【００４５】
先ず、図１０（ｂ）に示されている信号Ｉ０〜Ｉ７０π／８に対して、ステップＳＴ５２で算出された振幅ａｍｐｕｎｉおよび位相θｕｎｉを付加する。具体的には、信号Ｉ０〜Ｉ７０π／８の振幅Ａ０にａｍｐｕｎｉを乗算し、信号Ｉ０〜Ｉ７０π／８の位相にθｕｎｉを加算する。図１０（ｃ）に、信号Ｉ０〜Ｉ７０π／８に振幅ａｍｐｕｎｉおよび位相θｕｎｉを付加することにより得られた信号を、符号「Ｉ０′」〜「Ｉ７０π／８′」で示す。信号Ｉ０′〜Ｉ７０π／８′は、具体的には、以下の式で表される。
【００４６】
Ｉ０′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋θｕｎｉ）｝
Ｉ１０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋１０π／８＋θｕｎｉ）｝
Ｉ２０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋２０π／８＋θｕｎｉ）｝
Ｉ３０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋３０π／８＋θｕｎｉ）｝
Ｉ４０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋４０π／８＋θｕｎｉ）｝
Ｉ５０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋５０π／８＋θｕｎｉ）｝
Ｉ６０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋６０π／８＋θｕｎｉ）｝
Ｉ７０π／８′＝Ａ０・ａｍｐｕｎｉ・ｅｘｐ｛ｉ（θ０＋７０π／８＋θｕｎｉ）｝
【００４７】
次に、図１０（ａ）の信号Ｉ０〜Ｉ１４π／８と、図１０（ｃ）の信号Ｉ０′〜Ｉ７０π／８′とを、各チャンネルごとに加算する。図１０（ｄ）に、加算により得られた各チャンネルの信号ＪＩ〜Ｊ８を示す。このようにして、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための各チャンネルの信号ＪＩ〜Ｊ８を算出することができる。信号Ｊ１〜Ｊ８を算出したら、ＲＦシミングのフローを終了する。
【００４８】
第１の形態では、送信コイル２４に供給される信号のパターンを変更し、信号のパターンごとに送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を求めている（ステップＳＴ２）。そして、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、磁場の均一性が最も高い送信磁場分布Ｂ１２π／８を選択している（ステップＳＴ３）。ただし、被検体１２の電気的特性などが原因で、送信磁場分布Ｂ１２π／８の中にも、感度の低い部分と、感度の高い部分が現れるので、このような感度のばらつきは、できるだけ小さくすることが望ましい。そこで、第１の形態では、送信磁場分布Ｂ１２π／８だけでなく、送信磁場分布Ｂ１２π／８との相関が最も弱い送信磁場分布Ｂ１１０π／８も求め（ステップＳＴ４）、送信磁場分布Ｂ１２π／８と送信磁場分布Ｂ１１０π／８との合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）を求めている（ステップＳＴ５１）。送信磁場分布Ｂ１１０π／８は、送信磁場分布Ｂ１２π／８との相関が最も弱いので、送信磁場分布Ｂ１１０π／８を用いることによって、送信磁場分布Ｂ１２π／８の中の感度の低い部分を、感度が高くなるように補正することができ、更に、送信磁場分布Ｂ１２π／８の中の感度の高い部分を、感度が低くなるように補正することができる。したがって、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）が最も均一になるときの振幅ａｍｐおよび位相θの組合せを決定することによって、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための信号ＪＩ〜Ｊ８を算出することができる。
【００４９】
また、第１の形態では、ステップＳＴ２で求めた８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８の中から、２つの送信磁場分布Ｂ１２π／８およびＢ１１０π／８だけを選択して、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）の式を決定している。したがって、８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を全て用いなくても、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式を決定することができるので、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数の数を少なくすることができる。第１の形態では、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数は２つ（振幅ａｍｐおよび位相θ）だけで済む。したがって、送信コイル２４に供給する信号ＪＩ〜Ｊ８を決定するために必要な計算時間を短縮することができる。
【００５０】
尚、第１の形態では、ステップＳＴ２において、チャンネルＣＨ１〜ＣＨ８に供給される信号の位相を変化させることによって、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させている。しかし、振幅Ａ０を変化させることによって、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させてもよい。また、振幅と位相の両方を変化させて、チャンネルＣＨ１〜ＣＨ８に供給される信号のパターンを変化させてもよい。
【００５１】
また、第１の形態では、マグネット２に内蔵された送信コイル２４のシミングを行う場合について説明されている。しかし、本発明は、マグネット２とは別に備えられた送信コイル（例えば、人体の特定の部位の撮影に特化した専用コイル）のシミングを行う場合にも適用することができる。
【００５２】
（２）第２の形態
第１の形態は、送信磁場分布Ｂ１Δθ（図８参照）を実測により求めている。しかし、送信磁場分布Ｂ１Δθを計算により求めてもよい。第２の形態では、送信磁場分布Ｂ１Δθを計算により求めてＲＦシミングを行う例について説明する。
【００５３】
図１１は、第２の形態におけるＲＦシミングの手順を示すフローである。
ステップＳＴ１では、図１に示すように、被検体１２をボア２１に搬入する。被検体１２をボア２１に搬入したら、ステップＳＴ１１に進む。
【００５４】
ステップＳＴ１１では、送信コイル２４の各チャンネルごとに送信磁場分布を求める。
図１２は、ステップＳＴ１１で実行されるスキャンを示す図、図１３は、各スキャンを実行するときに送信コイル２４の各チャンネルに供給される信号を示す図である。
【００５５】
ステップＳＴ１１では、図１３に示すように、スキャンＳＣ１′〜ＳＣ８′が実行される。各スキャンＳＣ１′〜ＳＣ８′を実行する場合、各チャンネルＣＨ１〜ＣＨ８には、図１３に示すような信号が供給される。
【００５６】
例えば、スキャンＳＣ１′を実行する場合、図１３（ａ）に示すように、チャンネルＣＨ１にのみ信号Ｉ０が供給される。スキャンＳＣ１′によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１（図１参照）が、スキャンＳＣ１′によって取得されたデータに基づいて、チャンネルＣＨ１の送信磁場分布Ｂ１１を作成する。したがって、スキャンＳＣ１′を実行することによって、チャンネルＣＨ１による送信磁場分布Ｂ１１を作成することができる。
【００５７】
また、スキャンＳＣ２′を実行する場合、図１３（ｂ）に示すように、チャンネルＣＨ２にのみ信号Ｉ０が供給される。スキャンＳＣ２′によって取得されたデータは、中央処理装置９に伝送される。中央処理装置９では、送信磁場分布作成手段９１が、スキャンＳＣ２′によって取得されたデータに基づいて、チャンネルＣＨ２の送信磁場分布Ｂ１２を作成する。したがって、スキャンＳＣ２′を実行することによって、チャンネルＣＨ２による送信磁場分布Ｂ１２を作成することができる。
【００５８】
以下同様に、スキャンＳＣ３′〜ＳＣ８′を実行する場合、図１４（ｃ）〜（ｈ）に示すように、各チャンネルＣＨ３〜ＣＨ８にのみ信号Ｉ０を供給することにより、チャンネルＣＨ３〜ＣＨ８の送信磁場分布Ｂ１３〜Ｂ１８を求めることができる。したがって、スキャンＳＣ１′〜ＳＣ８′を実行することによって、チャンネルＣＨ１〜ＣＨ８ごとに、送信磁場分布Ｂ１１〜Ｂ１８を求めることができる。送信磁場分布１１〜Ｂ１８を求めた後、ステップＳＴ２に進む。
【００５９】
ステップＳＴ２では、送信磁場分布作成手段９１が、ステップＳＴ１１で求めたチャンネルＣＨ１〜ＣＨ８ごとの送信磁場分布Ｂ１１〜Ｂ１８に基づいて、送信コイル２４に信号Ｉ０〜Ｉ７Δθを供給したときに得られる送信磁場分布Ｂ１Δθ（図８参照）を算出する。送信磁場分布Ｂ１Δθは、例えば以下の式を用いて算出することができる。
【００６０】
【数１】
送信磁場分布Ｂ１ｎ（ｎ＝１〜８）はステップＳＴ１１で求められているので、求めた送信磁場分布Ｂ１ｎを上記の式に代入することによって、送信磁場分布Ｂ１Δθ（Ｂ１２π／８〜Ｂ１１６π／８）を算出することができる。送信磁場分布Ｂ１Δθを算出したら、ステップＳＴ３に進む。
ステップＳＴ３〜ＳＴ５は、第１の形態と同じであるので、説明は省略する。
【００６１】
第２の形態のように、送信磁場分布Ｂ１Δθを計算により求めても、最も均一性の高い合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐｕｎｉ，θｕｎｉ）を実現するための信号ＪＩ〜Ｊ８を算出することができる。
【００６２】
また、計算により８個の送信磁場分布Ｂ１２π／８〜Ｂ１１６π／８を算出した後は、第１の形態と同様に、２つの送信磁場分布Ｂ１２π／８およびＢ１１０π／８だけを選択して、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式を決定している。したがって、合成送信磁場分布Ｂ１ｃｏｍｂ（ａｍｐ，θ）の式に必要な変数は２つ（振幅ａｍｐおよび位相θ）だけで済み、送信コイル２４に供給する信号ＪＩ〜Ｊ８を決定するために必要な計算時間を短縮することができる。
【符号の説明】
【００６３】
２ マグネット
３ テーブル
３ａ クレードル
５ シーケンサ
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 中央処理装置
１０ 操作部
１１ 表示部
１２ 被検体
９１ 送信磁場分布作成手段
９２ 送信磁場分布選択手段
９３ 算出手段
１００ ＭＲ装置
【特許請求の範囲】
【請求項１】
ｎチャンネルの送信コイルを有する磁気共鳴装置であって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成手段と、
前記送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択手段と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出手段と、
を有する磁気共鳴装置。
【請求項２】
前記送信コイルの各チャンネルに供給する信号のパターンを変更し、前記被検体をスキャンするスキャン手段を有し、
前記送信磁場分布作成手段は、前記スキャンにより得られたデータに基づいて、前記信号のパターンごとに、前記送信磁場分布を作成する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の位相を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記位相は、以下の式で表される、請求項３に記載の磁気共鳴装置。
θ＝θ０＋ｋ・Δθ （ｋ＝０〜ｎ−１）
ただし、Δθ＝２π・（Ｍ／８） （Ｍ＝１〜ｎ）
【請求項５】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の振幅を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項６】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の位相および振幅を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項７】
前記送信磁場作成手段は、
前記送信コイルの各チャンネルごとに送信磁場分布を作成し、前記各チャンネルごとに作成した前記送信磁場分布に基づいて、前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を作成する、請求項１に記載の磁気共鳴装置。
【請求項８】
前記送信コイルの各チャンネルごとの送信磁場分布のデータを取得するためのスキャンを実行するスキャン手段を有する、請求項７に記載の磁気共鳴装置。
【請求項９】
前記算出手段は、
前記第１の送信磁場分布と前記第２の送信磁場分布との合成送信磁場分布を作成し、前記合成送信磁場分布が最も均一になるときの振幅および位相を算出し、前記算出した振幅および位相に基づいて、前記各チャンネルに供給する電流を算出する、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１０】
ｎチャンネルの送信コイルを有する磁気共鳴装置のプログラムであって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成処理と、
前記送信磁場分布作成処理により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択処理と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出処理と、
を計算機に実行させるためのプログラム。
【請求項１】
ｎチャンネルの送信コイルを有する磁気共鳴装置であって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成手段と、
前記送信磁場分布作成手段により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択手段と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出手段と、
を有する磁気共鳴装置。
【請求項２】
前記送信コイルの各チャンネルに供給する信号のパターンを変更し、前記被検体をスキャンするスキャン手段を有し、
前記送信磁場分布作成手段は、前記スキャンにより得られたデータに基づいて、前記信号のパターンごとに、前記送信磁場分布を作成する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の位相を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記位相は、以下の式で表される、請求項３に記載の磁気共鳴装置。
θ＝θ０＋ｋ・Δθ （ｋ＝０〜ｎ−１）
ただし、Δθ＝２π・（Ｍ／８） （Ｍ＝１〜ｎ）
【請求項５】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の振幅を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項６】
前記スキャン手段は、前記送信コイルに信号を供給する送信器を有し、
前記送信器は、前記信号の位相および振幅を変更することによって、前記信号のパターンを変更する、請求項２に記載の磁気共鳴装置。
【請求項７】
前記送信磁場作成手段は、
前記送信コイルの各チャンネルごとに送信磁場分布を作成し、前記各チャンネルごとに作成した前記送信磁場分布に基づいて、前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を作成する、請求項１に記載の磁気共鳴装置。
【請求項８】
前記送信コイルの各チャンネルごとの送信磁場分布のデータを取得するためのスキャンを実行するスキャン手段を有する、請求項７に記載の磁気共鳴装置。
【請求項９】
前記算出手段は、
前記第１の送信磁場分布と前記第２の送信磁場分布との合成送信磁場分布を作成し、前記合成送信磁場分布が最も均一になるときの振幅および位相を算出し、前記算出した振幅および位相に基づいて、前記各チャンネルに供給する電流を算出する、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１０】
ｎチャンネルの送信コイルを有する磁気共鳴装置のプログラムであって、
前記送信コイルの各チャンネルに供給する信号のパターンを変更したときの前記送信コイルの送信磁場分布を、前記信号のパターンごとに作成する送信磁場分布作成処理と、
前記送信磁場分布作成処理により得られた複数の送信磁場分布の中から、送信磁場の均一性が最も高い第１の送信磁場分布と、前記第１の送信磁場分布との相関が最も弱い第２の送信磁場分布とを選択する送信磁場分布選択処理と、
前記第１の送信磁場分布と前記第２の送信磁場分布とに基づいて、前記送信コイルの各チャンネルに供給する信号を算出する算出処理と、
を計算機に実行させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【公開番号】特開２０１３−２７５１８（Ｐ２０１３−２７５１８Ａ）
【公開日】平成２５年２月７日（２０１３．２．７）
【国際特許分類】
【出願番号】特願２０１１−１６５１５１（Ｐ２０１１−１６５１５１）
【出願日】平成２３年７月２８日（２０１１．７．２８）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
【公開日】平成２５年２月７日（２０１３．２．７）
【国際特許分類】
【出願日】平成２３年７月２８日（２０１１．７．２８）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
[ Back to top ]