(57)【要約】【課題】取得したＢ１マップからＢ１不均一性を簡便に評価すると共にＢ１不均一性を容易に調整可能な評価指標を算出することができる磁気共鳴イメージング装置を提供する。【解決手段】本実施形態の磁気共鳴イメージング装置は、人体を模擬したファントムに対して高周波磁場を印加して得られる前記ファントムのＲＦ受信信号の強度分布に基づいてＢ１マップを生成するＢ１マップ生成部と、前記Ｂ１マップの長軸又は短軸を抽出し、抽出した前記長軸の水平方向に対する長軸傾斜角、又は抽出した前記短軸の垂直方向に対する短軸傾斜角を前記Ｂ１マップの評価指標として算出する指標算出部と、入力される高周波パルスの位相および振幅の少なくとも一方を、算出された前記評価指標を使用して調整することによって高周波磁場の均一性を調整できるＱＤコイルと、を備えたことを特徴とする。【選択図】  図３
2014236868.tif 000002
【特許請求の範囲】
【請求項１】
人体を模擬したファントムに対して高周波磁場を印加して得られる前記ファントムのＲＦ受信信号の強度分布に基づいてＢ１マップを生成するＢ１マップ生成部と、
前記Ｂ１マップの長軸又は短軸を抽出し、抽出した前記長軸の水平方向に対する長軸傾斜角、又は抽出した前記短軸の垂直方向に対する短軸傾斜角を前記Ｂ１マップの評価指標として算出する指標算出部と、
入力される高周波パルスの位相および振幅の少なくとも一方を、算出された前記評価指標を使用して調整することによって高周波磁場の均一性を調整できるＱＤコイルと、
を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記Ｂ１マップ生成部は、
人体の一部を模したファントムの所定の断面を撮像して前記Ｂ１マップを生成し、
前記指標算出部は、
前記Ｂ１マップのコンタ形状が略楕円形状であるとして、その前記長軸又は短軸を抽出する、
ことを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記指標算出部は、
前記所定の断面の水平方向をＸ軸とし、前記Ｘ軸に直交する軸をＹ軸とするとき、前記Ｙ軸方向において画素値が最大となるＹ座標の位置Ｙmaxを、前記Ｘ軸に沿った各位置Ｘに対して求め、求めた各位置（Ｘ、Ｙmax）から近似される直線を前記長軸として抽出する、
ことを特徴とする請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記指標算出部は、
前記Ｂ１マップ内の最大画素値を抽出し、抽出した最大画素値のＸ座標を中心とした所定の範囲の各Ｘに対して、前記Ｙ軸方向において画素値が最大となるＹ座標の位置Ｙmaxを求め、求めた各位置（Ｘ、Ｙmax）を最小二乗法で近似する直線を前記長軸として抽出する、
ことを特徴とする請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
前記指標算出部は、
前記所定の断面の水平方向をＸ軸とし、前記Ｘ軸に直交する軸をＹ軸とするとき、前記Ｘ軸方向において画素値が最大となるＸ座標の位置Ｘmaxを、前記Ｙ軸に沿った各位置Ｙに対して求め、求めた各位置（Ｘmax、Ｙ）から近似される直線を前記短軸として抽出する、
ことを特徴とする請求項２に記載の磁気共鳴イメージング装置。
【請求項６】
前記指標算出部は、
前記Ｂ１マップ内の最大画素値を抽出し、抽出した最大画素値のＹ座標を中心とした所定の範囲の各Ｙに対して、前記Ｘ軸方向において画素値が最大となるＸ座標の位置Ｘmaxを求め、求めた各位置（Ｘmax、Ｙ）を最小二乗法で近似する直線を前記短軸として抽出する、
ことを特徴とする請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記指標算出部は、
前記短軸傾斜角と前記長軸傾斜角とを重み付け加算して合成傾斜角を算出し、前記合成傾斜角を前記Ｂ１マップの評価指標とする、
ことを特徴とする請求項１乃至６のいずれか１項に記載の磁気共鳴イメージング装置。
【請求項８】
前記Ｂ１マップを表示する表示部、をさらに備え、
前記表示部は、算出した前記評価指標をさらに表示する、
ことを特徴とする請求項１乃至７のいずれか１項に記載の磁気共鳴イメージング装置。
【請求項９】
前記所定の断面はアキシャル断面である、
ことを特徴とする請求項２乃至６のいずれか１項に記載の磁気共鳴イメージング装置。
【請求項１０】
人体を模擬したファントムに対して高周波磁場を印加して得られる前記ファントムのＲＦ受信信号の強度分布に基づいてＢ１マップを生成し、
前記Ｂ１マップの長軸又は短軸を抽出し、抽出した前記長軸の水平方向に対する長軸傾斜角、又は抽出した前記短軸の垂直方向に対する短軸傾斜角を前記Ｂ１マップの評価指標として算出し、
ＱＤコイルに入力される高周波パルスの位相および振幅の少なくとも一方を、算出された前記評価指標を使用して調整することによって高周波磁場の均一性を調整する、
ことを特徴とする磁気共鳴イメージング装置の高周波磁場均一性調整方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は磁気共鳴イメージング装置及びその高周波磁場均一性調整方法に関するに関する。
【背景技術】
【０００２】
磁気共鳴イメージングは、静磁場中に置かれた被検体の原子核スピンをラーモア周波数のＲＦ信号で磁気的に励起し、この励起に伴って発生する磁気共鳴信号から画像を再構成する撮像法である。
【０００３】
この磁気共鳴イメージングにおいて、送信されるパルスの周波数が増大すると、被検体内におけるＲＦ磁場が不均一となる。このＲＦ磁場の不均一性は、Ｂ１不均一性とも呼ばれる。そこで、Ｂ１不均一性を改善することが重要となる。
【０００４】
例えば、アスペクト比に応じてＲＦ送信パルスの振幅および位相を決定するＢ１補正方法が考案されている(例えば特許文献１参照)。また、ＲＦ送信パルスの振幅および位相を変えながらＢ１マップを収集し、ＲＦ送信パルスの振幅および位相を決定する方法も知られている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２００８−５１４２５９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、近年、磁場強度が高い磁気共鳴イメージング装置、例えば３Ｔ（テスラ）のＭＲＩ装置が開発され、実用化されてきている。このため、ＲＦパルスの周波数が一層大きくなり、従来のようにアスペクト比に応じて求めた振幅および位相でRFパルスを送信しても、Ｂ１不均一性が解消できない場合がある。一方、従来のようにＲＦパルスの振幅および位相を変化させながらＢ１マップを収集する場合には、撮像時間が長くなるという問題やＢ１不均一性の判定が難しいという問題点がある。特に、Ｂ１マップを取得したとしても、そのＢ１マップの評価指標は必ずしも確立されているとは言えず、取得したＢ１マップから経験的に、或いは試行錯誤的にＢ１シミングを行っているのが実態である。
【０００７】
そこで、取得したＢ１マップからＢ１不均一性を簡便に評価すると共に、Ｂ１不均一性を容易に調整可能な評価指標を算出することができる磁気共鳴イメージング装置、及びそのＢ１不均一性調整方法が望まれている。
【課題を解決するための手段】
【０００８】
本実施形態の磁気共鳴イメージング装置は、人体を模擬したファントムに対して高周波磁場を印加して得られる前記ファントムのＲＦ受信信号の強度分布に基づいてＢ１マップを生成するＢ１マップ生成部と、前記Ｂ１マップの長軸又は短軸を抽出し、抽出した前記長軸の水平方向に対する長軸傾斜角、又は抽出した前記短軸の垂直方向に対する短軸傾斜角を前記Ｂ１マップの評価指標として算出する指標算出部と、入力される高周波パルスの位相および振幅の少なくとも一方を、算出された前記評価指標を使用して調整することによって高周波磁場の均一性を調整できるＱＤコイルと、を備えたことを特徴とする。
【図面の簡単な説明】
【０００９】
【図１】実施形態の磁気共鳴イメージング装置の全体構成例を示す構成図。
【図２】ＲＦコイル（ＱＤコイル）および送信器の詳細構成例を示す図。
【図３】プロセッサで実現する各種機能のうち、実施形態のＢ１シミングに関連する機能を抽出した機能ブロック図。
【図４】実施形態の磁気共鳴イメージング装置の動作のうち、Ｂ１シミングに関連する動作例を示すフローチャート。
【図５】ファントムの一例を示す図。
【図６】ファントムのアキシャル断面のＢ１マップを模式的に例示した第１の図。
【図７】ファントムのアキシャル断面のＢ１マップに長軸及び短軸を付加した図。
【図８】横長略楕円形状のＢ１マップから長軸を抽出する処理の説明図。
【図９】Ｂ１マップから長軸を抽出する処理例を示すフローチャート。
【図１０】横長略楕円形状のＢ１マップから短軸を抽出する処理の説明図。
【図１１】Ｂ１マップから短軸を抽出する処理例を示すフローチャート。
【発明を実施するための形態】
【００１０】
以下、本発明の実施形態を添付図面に基づいて説明する。
【００１１】
（１）全体構成
図１は、本実施形態における磁気共鳴イメージング装置１の全体構成を示すブロック図である。実施形態の磁気共鳴イメージング装置１は、磁石架台１００、寝台２００、制御キャビネット構成品３００、コンソール４００等を備えて構成される。
【００１２】
磁石架台１００は、静磁場磁石１０、傾斜磁場コイル１１、ＲＦコイル１２等を有しており、これらの構成品は円筒状の筐体に収納されている。寝台２００は、寝台本体２０と天板２１を有している。
【００１３】
一方、制御キャビネット構成品３００は、静磁場用電源３０、傾斜磁場電源３１（Ｘ軸用３１ｘ、Ｙ軸用３１ｙ、Ｚ軸用３１ｚ）、ＲＦ受信器３２、ＲＦ送信器３３、シーケンスコントローラ３４等を備えている。また、コンソール４００は、プロセッサ４０、記憶部４１、入力部４３、表示部４２等を有するコンピュータとして構成されている。
【００１４】
磁石架台１００の静磁場磁石１０は、概略円筒形状をなしており、被検体（患者）の撮像領域であるボア（静磁場磁石１０の円筒内部の空間）内に静磁場を発生させる。静磁場磁石１０は超電導コイルを内蔵し、液体ヘリウムによって超電導コイルが極低温に冷却されている。静磁場磁石１０は、励磁モードにおいて静磁場用電源３０から供給される電流を超電導コイルに印加することで静磁場を発生し、その後、永久電流モードに移行すると、静磁場用電源３０は切り離される。一旦永久電流モードに移行すると、静磁場磁石１０は長時間、例えば１年以上に亘って、大きな静磁場を発生し続ける。なお、静磁場磁石１０を永久磁石として構成しても良い。
【００１５】
傾斜磁場コイル１１も概略円筒形状をなし、静磁場磁石１０の内側に固定されている。この傾斜磁場コイル１１は、傾斜磁場電源（３１ｘ、３１ｙ、３１ｚ）から供給される電流によりＸ軸，Ｙ軸，Ｚ軸の方向に傾斜磁場を被検体に印加する。
【００１６】
寝台２００の寝台本体２０は天板２１を上下方向に移動可能であり、撮像前に天板２１に載った被検体を所定の高さまで移動させる。その後、撮影時には天板２１を水平方向に移動させて被検体をボア内に移動させる。
【００１７】
ＲＦコイル１２は全身用コイルとも呼ばれ、傾斜磁場コイル１１の内側に被検体を挟んで対向するように固定されている。ＲＦコイル１２は、ＲＦ送信器３３から伝送されるＲＦパルスを被検体に向けて送信する一方、また、水素原子核の励起によって被検体から放出される磁気共鳴信号を受信する。
【００１８】
ＲＦ送信器３３は、シーケンスコントローラ３４からの指示に基づいて、ＲＦコイル１２にＲＦパルスを送信する。一方、ＲＦ受信器３２は、ＲＦコイル１２によって受信された磁気共鳴信号を検出し、検出した磁気共鳴信号をデジタル化して得られる生データをシーケンスコントローラ３４に対して送信する。
【００１９】
シーケンスコントローラ３４は、コンソール４００による制御のもと、傾斜磁場電源３１、ＲＦ送信器３３およびＲＦ受信器３２をそれぞれ駆動することによって被検体のスキャンを行う。そして、シーケンスコントローラ３４は、スキャンを行ってＲＦ受信器３２から生データを受信すると、その生データをコンソール４００に送信する。
【００２０】
コンソール４００は、ＭＲＩ装置１全体を制御する。具体的には、検査技師等のマウスやキーボード等（入力部４２）の操作によって撮像条件その他の各種情報や指示を受け付ける。そして、プロセッサ４０は、入力された撮像条件に基づいてシーケンスコントローラ３４にスキャンを実行させる一方、シーケンスコントローラ３４から送信された生データに基づいて画像を再構成する。再構成された画像は表示部４２に表示され、或いは記憶部４１に保存される。
【００２１】
なお、図１において人体に付した黒丸は、磁場中心（アイソセンタ）を表している。
【００２２】
図２は、図１に示すＲＦコイル１２および送信器３３の詳細構成例を示す図である。
【００２３】
ＲＦコイル１２はＱＤコイルとして構成される。ＱＤコイルには、通常２つの送信チャンネルｃｈ１、ｃｈ２から送信信号が給電され、第１の送信チャンネルｃｈ１と第２の送信チャンネルｃｈ２によって、空間的に互いに直交する磁界を発生する。ＱＤコイルには、円筒状のバードケージタイプ型のものや、独立した２つのコイルで構成されるタイプのものがある。バードケージタイプのＱＤコイルでは、図２に示すように、円筒上で互いに９０度離れた位置にある２つの給電点に、チャンネルｃｈ１とチャンネルｃｈ２の送信信号が夫々給電される。一方、独立２コイルタイプのＱＤでは、それぞれのコイルにチャンネルｃｈ１とチャンネルｃｈ２の送信信号が別々に給電される。
【００２４】
ＲＦ送信器３３では、ＲＦ送信信号を、９０度分配器３３１によって電気的に９０度位相の異なる第１及び第２の送信チャンネルｃｈ１、ｃｈ２の信号に分配した後、各送信ｃｈ１， ｃｈ２のＲＦ送信信号をＲＦアンプ（２）３３０Ａ、ＲＦアンプ（１）３３０Ｂでそれぞれ増幅するように構成されている。
【００２５】
ＱＤコイルは、送信時においては、電気的位相が９０度離れた２つのＲＦ信号を空間的に直交する方向にそれぞれ印加することになり、理想的な条件下においては、撮像領域のＸＹ面に均一な高周波磁場を形成することができる。
【００２６】
しかしながら、人体を撮像する場合、誘電効果により高周波磁場の不均一性が発生することが知られている。この高周波磁場の不均一性はＢ１不均一性とも呼ばれている。特に３Ｔ等の強い静磁場を発生する磁気共鳴イメージング装置１では、Ｂ１不均一性がより大きくなる。そして、このＢ１不均一性は、撮像対象部位の形状、大きさ等によって変化する。
【００２７】
そこで、高周波磁場の均一性を調整する手段として、第１及び第２の送信チャンネルｃｈ１、ｃｈ２には、位相調整部３３２Ａ、３３２Ｂ、及び利得調整部３３３Ａ、３３３Ｂが夫々設けられている。なお、９０度分配器３３１、位相調整部３３２Ａ、３３２Ｂ及び利得調整部３３３Ａ、３３３Ｂで、振幅／位相調整部３４０を構成する。
【００２８】
位相調整部３３２Ａ、３３２Ｂと利得調整部３３３Ａ、３３３Ｂによって、チャンネル間の相対的な位相（位相バランス）と相対的な利得（利得バランス）を調整することができる。位相調整部３３２Ａ、３３２Ｂと利得調整部３３３Ａ、３３３Ｂは、夫々のチャンネルに別々に設けても良いし、いずれか一方のチャンネルにのみ設けても良い。
【００２９】
以下の説明では、図２に示すように、第１の送信チャンネルｃｈ１の位相φ１と利得Ｇ１を固定し、第２の送信チャンネルｃｈ２の位相φ２と利得Ｇ２を可変にしてチャンネル間の位相差と利得差を調整するものとする。具体的には、第２の送信チャンネルｃｈ２の位相φ２を、φ２＝φ１＋α、とする一方、利得Ｇ２を、Ｇ２＝Ｇ１＋βとし、αとβを調整することにより、チャンネル間の位相差と利得差を調整するものとしている。なお、位相φ２と利得Ｇ２のうち、いずれか一方のみを調整可能な構成としても良い。例えば、位相φ２のみを調整する構成としてもよい。
【００３０】
なお、位相調整部３３２Ａ、３３２Ｂや利得調整部３３３Ａ、３３３Ｂを調整することによって高周波磁場を均一化する調整を以下ではＢ１シミングと呼ぶものとする。
【００３１】
（２）Ｂ１シミング
図３は、図１に示すプロセッサ４０で実現する各種機能のうち、実施形態のＢ１シミングに関連する機能を抽出した機能ブロック図である。図３に示す撮像部５００は、図１の磁石架台１００と制御キャビネット構成品３００で構成されるが、一部を省略している。
【００３２】
実施形態の磁気共鳴イメージング装置１は、前述したようにＱＤコイルとして構成されるＲＦコイル１２、Ｂ１マップ生成部６０、評価指標算出部７０を少なくとも有して構成される。このうち、Ｂ１マップ生成部６０及び評価指標算出部７０は、コンソール４００のプロセッサ４０が所定のプログラムを実行することによって実現される機能構成である。
【００３３】
図４は、実施形態の磁気共鳴イメージング装置１の動作のうち、Ｂ１シミングに関連する動作例を示すフローチャートである。以下、図４のフローチャートに沿って、図３の機能ブロック図を参照しつつ実施形態の磁気共鳴イメージング装置１の動作を説明する。
【００３４】
まず、ステップＳＴ１０で、人体を模擬したファントム６００に対して撮像（Ｂ１マップ撮像）を行う。ここで用いるファントム６００は、人体の誘電効果とほぼ同一の誘電効果を得るような物質を人体の形状を模擬した容器に入れたものである。
【００３５】
例えば、導電率、誘電率、Ｔ１値、Ｔ２値等が人体と近似するような、硫酸銅、塩化ニッケルなどの水溶液を容器に入れてファントム６００を作成する。容器の形状は必ずしも人体全体の形状を模したものである必要はない。例えば、仰臥位又は腹臥位の状態で天板２１に載った人体のうち、頭部、胸部、腹部等の所定の撮像部位に対して、その左右方向（ＲＬ方向）及び背腹方向（ＡＰ方向）の寸法や形状を模擬した容器に前記の水溶液を収納したものを撮像対象のファントム６００とする。
【００３６】
図５は、ＸＹ断面が略楕円形状をなすファントム６００の一例を示す図である。ファントム６００のＸ方向の長さは、仰臥位又は腹臥位の人体の所定の撮像部位のＲＬ方向の長さに対応し、Ｙ方向の長さはＡＰ方向に対応する。Ｚ方向の長さは必ずしも厳密に規定する必要はなく、アイソセンタ近傍において、所定の厚みのスラブ、或いはスライス（アキシャル断面（ＸＹ断面））を撮像するのに十分な長さであれば良い。
【００３７】
図４のステップＳＴ１０では、上記のようなファントム６００を天板２１上のアイソセンタ近傍に設置し、そのファントム６００に対してＢ１マップ撮像を行う。
【００３８】
Ｂ１マップ撮像は、例えば、スピンエコー法に基づく公知のダブルアングル法などを用いる。具体的には、励起パルスとリフォーカスパルスのフリップ角を９０°と１８０°とする第１の撮像と、励起パルスとリフォーカスパルスのフリップ角を４５°と９０°とする第２の撮像をファントム６００に対して行う。第１の撮像と第２の撮像では、上記のフリップ角以外の撮像条件は共通にして行われる。
【００３９】
上記のダブルアングル法などで収集されたデータは、シーケンスコントローラ３４を介してＢ１マップ生成部６０に送られ、ここでＢ１マップが生成される（図４のステップＳＴ１１）。
【００４０】
Ｂ１マップとはＲＦ出力（Ｂ１値）の強度分布のことであり、通常、Ｂ１マップは、ＲＦ受信信号（ＭＲ信号）の強度分布に基づいて作成される。ＲＦ受信信号の強度分布として、実際のＲＦ信号の強度（絶対値）を用いて表しても良いし、或いは、相対値を用いて表してもよい。以下では、Ｂ１値は、本来出力されるべきＲＦ強度に対して、実際に観測されたＲＦ強度の比（相対値）であるものとして説明する。２次元分布としてＢ１マップを生成する場合、Ｂ１マップの断面は任意の方向の断面を取りうるが、一般に、アキシャル断面の２次元分布（アキシャル断面のＢ１マップ）を生成することが多い。そこで、以下ではアキシャル断面のＢ１マップを例にして説明する。
【００４１】
図６は、図５に示したファントム６００のアキシャル断面のＢ１マップを模式的に例示した図である。一般に、Ｂ１マップのＢ１値は中心部で高く、Ｂ１値のコンタ形状（Ｂ１値の等高線形状）は横長の楕円形に近い形状となる。図６では、Ｂ１値の高い領域が濃く、Ｂ１値の低い領域が薄くなるようにＢ１値の高低を濃淡で表現している。
【００４２】
ファントム６００自体は水平に設置されているため、Ｂ１均一性が確保されていれば、本来、Ｂ１マップの横長楕円形状はほぼ水平となるはずである。しかしながら、Ｂ１不均一性の影響によって横長楕円形状は図６に例示するように水平方向に対してある傾きをもつ。このため、Ｂ１マップのＢ１値を水平方向に見た場合、Ｘ方向の中心線Ｘ’−Ｘ″に対して対称とはならない。例えば、Ｂ１マップの上側の領域にあるＹａ−Ｙａ’線上では、中心線Ｘ’−Ｘ″の右側の領域の方がＢ１値が高く、左側の領域ではＢ１値が低くなっている。また、Ｂ１マップの下側の領域にあるＹｂ−Ｙｂ’線上では、逆に中心線Ｘ’−Ｘ″の右側の領域の方がＢ１値が低く、左側の領域ではＢ１値が高くなっている。このような左右のアンバランスは、人体を撮像した場合、画像診断上好ましくなく、場合によっては診断を誤る可能性もある。
【００４３】
そこで、本実施形態の磁気共鳴イメージング装置１では、図７に示すように、Ｂ１マップから長軸または短軸を抽出し、Ｂ１不均一性に起因するＢ１マップのアンバランスを表す評価指標として、長軸又は短軸の傾斜角を採用するものとしている。つまり、Ｂ１マップのコンタ形状が略楕円形状であると仮定してその長軸又は短軸を抽出した上で（図４のステップＳＴ１２）、長軸の水平方向に対する傾斜角の絶対値（長軸傾斜角）や、短軸の垂直方向に対する傾斜角の絶対値（短軸傾斜角）を算出し（ステップＳＴ１３）、算出した傾斜角をＢ１マップの評価指標としている。
【００４４】
以下、より具体的な算出方法について図８乃至図１１を参照しつつ説明する。図８は、横長略楕円形状のＢ１マップから長軸を抽出する処理の説明図であり、図９はその処理のフローチャートである。
【００４５】
図８中の中央が盛り上がった太い曲線は、Ｙ方向におけるＢ１値（Ｂ１マップの画素値）の最大値を、Ｘ軸に沿った各位置Ｘ（Ｘ方向のピクセル位置）対して求め、求めたＢ１値の最大値をＸ座標に対してプロットしたものである。また、図８中の太い破線は、Ｂ１値が最大となるＹ座標の位置Ｙmaxを、前記Ｘ軸に沿った各位置Ｘに対して求め、座標（Ｘ、Ｙmax）で表される連続点を、横軸Ｘ座標、縦軸Ｙ座標としてプロットしたものである。上記の、Ｙ方向におけるＢ１値の最大値をＸ軸に沿った各位置Ｘ対して求める処理、及び、Ｂ１値が最大となるＹ座標の位置ＹmaxをＸ軸に沿った各位置Ｘに対して求める処理は、例えば、図９のフローチャートのステップＳＴ２０乃至ステップＳＴ２４の処理によって行われる。
【００４６】
ステップＳＴ２５では、全領域での最大Ｂ１値（図８の黒丸で示されるＢ１値）のＸ座標を中心とする所定の範囲内（図８の２つの細い矢印で示される範囲）で、連続点（Ｘ、Ｙmax）を最小二乗近似法で直線近似する。そして、このようにして得られた近似直線を、Ｂ１マップの近似楕円形の長軸とする（ステップＳＴ２６）。この長軸は、図８では、中央部の太い実線で表現されている。
【００４７】
上記の処理による長軸の抽出は、図３の機能ブロック図においては、評価指標算出部７０の長軸／短軸抽出部７１が行う。また、この長軸の抽出処理は、図４のフローチャートにおいては、ステップＳＴ１２の処理で行う。
【００４８】
長軸が抽出されると、評価指標算出部７０の傾斜角算出部７２は長軸傾斜角θＬを算出する（ステップＳＴ１３）。具体的には、図８に示すように、求めた長軸の水平方向（Ｘ軸方向）に対する角度を算出し、この角度（長軸傾斜角θＬ）を、Ｂ１不均一性の評価指標とする。なお、長軸の水平方向に対する角度は正負両方の値を取り得るが、この角度の絶対値を長軸傾斜角θＬと呼ぶものとする。
【００４９】
実施形態の磁気共鳴イメージング装置１は、Ｂ１マップを表示部４２に表示するとき、求めた長軸を図７のようにＢ１マップに重ねて表示してもよく、さらに、算出した評価指標としての長軸傾斜角θＬを表示部４２の適宜の場所に表示しても良い（ステップＳＴ１４）。
【００５０】
このようにして、取得したＢ１マップから、その評価指標としての長軸傾斜角θＬが自動的に算出され、そして表示部４２に表示される、ユーザは長軸傾斜角θＬの値に基づいて、Ｂ１シミングを実施する（図４のステップＳＴ１５）。具体的には、ユーザインタフェース（入力部４３等）を介して、振幅／位相調整部３４０の位相調整部３３２Ａや利得調整部３３３Ａに、長軸傾斜角θＬに基づく調整値を設定することにより、Ｂ１均一化の調整、即ち、Ｂ１シミングを実施する。より具体的には、長軸傾斜角θＬがゼロに近づくように、位相差の調整値αや、振幅差の調整値βを調整する。
【００５１】
なお、どの程度の長軸傾斜角θＬがあったときに、どの程度位相調整部３３２Ａの位相差をシフトさせれば（或いは、どの程度利得調整部３３３Ａの利得差を増減させれば）Ｂ１マップの楕円の傾斜を水平にすることができるか（即ち、長軸傾斜角θＬをゼロにすることができるか）、といった情報は、過去の実績等から既に得られている。例えば、長軸傾斜角θＬが１５度あった場合、位相調整部３３２Ａの位相差を、長軸傾斜角θＬの１／１０程度（即ち、１．５度程度）シフトさせれば、Ｂ１マップの楕円形状の傾きをほぼ水平にし、Ｂ１シミングを行うことができる、といった情報或いは知見が、過去の実績から得られている。
【００５２】
したがって、取得したＢ１マップから長軸傾斜角θＬが算出されさえすれば、ユーザは長軸傾斜角θＬからＢ１シミングを実施することが可能である。
【００５３】
また、長軸傾斜角θＬと、位相調整部３３２Ａの位相差のシフト量（或いは、利得調整部３３３Ａの利得差を増減量）との関係を、従来の実績値に基づいてデータベース化しておけば、装置がこのデータベースを参照することにより、長軸傾斜角θＬから位相調整部３３２Ａの位相調整量αや利得調整部３３３Ａの利得調整量βを装置が求めることができ、Ｂ１マップの取得からＢ１シミングの実施までを自動化することも可能である。
【００５４】
ここまでは、Ｂ１マップの評価指標として長軸傾斜角θＬを用いることを説明したが、長軸傾斜角θＬに換えて短軸傾斜角θＳを評価指標としてもよい。この場合、図７に示す短軸をＢ１マップから抽出し、抽出した短軸の垂直方向（Ｙ軸方向）に対する傾斜角の絶対値が、評価指標としての短軸傾斜角θＳとなる。
【００５５】
図１０は横長略楕円形状のＢ１マップから短軸を抽出する処理の説明図であり、図１１はその処理のフローチャートである。図１０は、図８のＸ軸とＹ軸とを入れ替えたものであり、それ以外の説明は図８と重複するため、具体的な説明は省略する。また、図１１のステップＳＴ３０乃至ステップＳＴ３６の処理も、図９のステップＳＴ２０乃至ステップＳＴ２６の処理に対してＸ軸とＹ軸とを入れ替えたものであるため、具体的な説明を省略する。
【００５６】
図１１のステップＳＴ３６においてＢ１マップの短軸が求まると、評価指標算出部７０の傾斜角算出部７２は短軸傾斜角θＳを算出する（図４のステップＳＴ１３）。具体的には、図１０に示すように、求めた短軸の垂直方向（Ｙ軸方向）に対する角度を算出し、この角度（短軸傾斜角θＳ）を、Ｂ１不均一性の評価指標とする。
【００５７】
そして、図４のステップＳＴ１５において、評価指標としての短軸傾斜角θＳに基づいて、評価指標を長軸傾斜角θＬとした場合と同様に、位相調整部３３２Ａと利得調整部３３３Ａに対して位相／振幅調整を行ってＢ１シミングを実施する。
【００５８】
上述した実施形態の説明では、長軸傾斜角θＬ或いは短軸傾斜角θＳのいずれか一方をＢ１マップの評価指標としている。
【００５９】
これに対して、長軸傾斜角θＬと短軸傾斜角θＳの両方を組み合わせて１つの合成評価指標を作成し、この合成評価指標を用いてＢ１シミングをするようにしても良い。
【００６０】
前述したように長軸や短軸の抽出の際には、Ｂ１値の大きな中央部分の領域のデータに対して最小二乗近似法を用いて近似直線を求めており、Ｂ１マップの外側近傍における歪んだ領域（楕円形状に近似できない領域）のデータの使用を避けている。しかしながら、このようにしても、上記手順で算出した長軸傾斜角θＬや短軸傾斜角θＳの値は、Ｂ１マップの中心から離れた領域における多少歪んだ楕円形状の影響を受ける可能性がある。
【００６１】
そこで、歪んだ楕円形状の影響をさらに低減する手法として、長軸傾斜角θＬと短軸傾斜角θＳを重み付け加算して合成傾斜角θを算出し、この合成傾斜角θを評価指標として使用する実施形態が有効である。具体的には、合成傾斜角θを次式で算出する。
【００６２】
θ＝ｗ・θＬ＋（１−ｗ）・θＳ    （０≦ｗ≦１）
上式において、ｗは重み係数であり、ｗが１に近いと長軸の重みが増し、ゼロに近いと短軸の重みが増す。重み係数ｗはユーザインターフェース（入力部４３）を介して設定することができ、Ｂ１マップのコンタ形状（楕円形状への近似の程度）等に応じて適宜変更することができる。
【００６３】
以上説明してきたように、実施形態の磁気共鳴イメージング装置１、及びそのＢ１不均一性調整方法によれば、取得したＢ１マップからＢ１不均一性を簡便に評価すると共に、Ｂ１不均一性を容易に調整可能な評価指標を算出することができる。
【００６４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００６５】
１  磁気共鳴イメージング装置
１２  ＲＦコイル
３４  シーケンスコントローラ
４２  表示部
４３  ユーザインタフェース（入力部）
６０  Ｂ１マップ生成部
７０  評価指標算出部
７１  長軸／短軸抽出部
７２  傾斜角算出部
３４０  振幅／位相調整部
４００  コンソール
５００  撮像部
【図１】
2014236868.tif 000003
【図２】
2014236868.tif 000004
【図３】
2014236868.tif 000005
【図４】
2014236868.tif 000006
【図５】
2014236868.tif 000007
【図８】
2014236868.tif 000008
【図９】
2014236868.tif 000009
【図１０】
2014236868.tif 000010
【図１１】
2014236868.tif 000011
【図６】
2014236868.tif 000012
【図７】
2014236868.tif 000013