JP6936696B2 - 磁気共鳴イメージング装置及びその制御方法 - Google Patents

Description

本発明は磁気共鳴イメージング装置（以下、ＭＲＩ装置という）に関し、特に体内にインプラント等の金属を含む被検体に対し、金属周辺に生じる画像の歪みを低減する技術に関する。

ＭＲＩ装置では、強い磁場内に被検体を置き、核磁気共鳴を起こさせる高周波磁場を印加し、印加終了後に被検体から発生する核磁気共鳴信号を検出する。この際、核磁気共鳴信号に位置情報を付加するために、静磁場に磁場勾配を与える。この傾斜磁場は、直交する３軸方向にパルス状に印加され、スライス傾斜磁場、位相エンコード傾斜磁場、周波数エンコード傾斜磁場などと呼ばれる。

このような強い磁場内で検査を行うＭＲＩでは、発熱や誘導電流を引き起こす金属は禁忌であったが、近年、インプラント（体内埋め込み型デバイス）を持った被検者に対するＭＲＩ撮像の要望が増加しており、条件付きＭＲＩ対応のインプラントも開発されている。しかし、条件付きＭＲＩ対応デバイスの適用により人体への影響を極力抑えることができるものの、依然として金属周辺の画質の劣化の問題がある。具体的には、静磁場により金属が磁化されることで磁場が発生し、その金属周辺の磁場が歪む。この金属周辺の磁場の歪みは、ＭＲＩで得られる画像において金属周辺で周波数方向及びスライス方向に強い歪みとなって現れる。

この歪みに対し、特許文献１に記載された技術では、スライス位置を選択する傾斜磁場を印加しないことでスライス方向の歪みを低減するとともに、高周波磁場の照射（送信）及び受信の周波数帯域を複数回シフトして複数の３次元画像を取得し、これらを合成することで周波数方向の歪みを低減している。一つの周波数帯域（Ｂｉｎ）の３次元画像をＢｉｎ画像ということとする。

米国特許第７８２１２６４号明細書

特許文献１に記載された技術では、照射の周波数帯域をシフトする際に、隣り合う周波数帯域の間隔を狭くし、隣り合う周波数帯域を重ねることで、合成後の画像に周波数帯域の谷間を発生させないようにしている。しかしながら、ある共鳴周波数をもつ1つのスピンは隣り合う複数のＢｉｎ画像に描出される。受信周波数中心はＢｉｎごとに異なるため、共鳴周波数と受信周波数中心の違いにより、1つのスピンはＢｉｎ画像ごとに周波数エンコード方向に異なる位置に描出される。これを合成すると周波数エンコード方向にボケた画像になる。特に、撮像範囲に対して金属が小さく、金属の影響を受けない静磁場中心周波数付近のスピンが画像の大部分を占める場合、画像全体が周波数エンコード方向にボケることになる。

そこで本発明は、周波数帯域の異なるＢｉｎ画像を合成して画像を作成する際に、周波数方向のボケを回避すること、特に異なるＢｉｎ画像間で静磁場中心周波数付近のスピンが周波数エンコード方向にずれるのを防止することを課題とする。

上記課題を解決するため、本発明は周波数帯域の異なるＢｉｎ画像を取得する際に、送信周波数の制御とは別に受信周波数を制御し、画像で支配的なスピンがＢｉｎ画像毎に周波数方向にずれるのを低減する。

具体的には、本発明のＭＲＩ装置は、静磁場発生磁石と、前記静磁場発生磁石が形成する静磁場空間に、傾斜磁場を発生する傾斜磁場コイルと、前記静磁場空間に置かれた被検体に、高周波磁場を印加する送信部と、前記被検体から発生される核磁気共鳴信号を受信する受信部と、前記高周波磁場の送信周波数を予め定めた周波数差で変化させて、複数の計測を行う制御部と、前記複数の計測で取得した複数の前記被検体の画像を合成する画像演算部と、を備える。前記制御部は、前記複数の計測のうち、基準周波数の計測に対し、周波数方向に隣接する周波数の計測において、前記基準周波数の計測との受信周波数差を０又は最小とするように、前記受信周波数を制御する。受信周波数の制御に代えて前記複数の計測で得た画像を補正することも可能である。

また本発明のＭＲＩ装置の制御方法は、複数の異なる送信周波数を用いて複数の計測を行い、当該複数の計測により得た複数の三次元画像を合成して１つの画像を取得する。複数の撮像の際の受信周波数を前記送信周波数の変化に応じて変化させる制御を行い、その際、基準周波数の計測に対し周波数方向に隣接する計測の受信周波数を、前記基準周波数との差が０又は最小となるように各計測における前記受信周波数を制御する。

本発明によれば、最も支配的なスピンの周波数（基準周波数という）、例えば静磁場中心周波数を中心周波数とするＢｉｎの計測に対し、予め定めた周波数差ずらしたＢｉｎの計測において、その信号受信の際の受信周波数を、基準周波数側にずらし、基準周波数近傍の受信周波数を基準周波数に一致させることで、最も支配的なスピンについて周波数方向のずれをなくすことができる。これにより複数のＢｉｎ画像を合成したときのボケを低減し画質の向上を図ることができる。

ＭＲＩ装置の一実施形態の全体構成を示すブロック図。 ＭＲＩ装置で実行されるパルスシーケンスの一例を示す図。 第一実施形態の制御部（ＣＰＵ）の構成を示す機能ブロック図。 第一実施形態の撮像の手順を示すフロー図。 撮像条件設定のためのＧＵＩの一例。 送信周波数（中心周波数）と受信周波数（中心周波数）との関係を示す図で、（ａ）はボケ補正なしの場合、（ｂ）は第一実施形態の場合、（ｃ）はその変形例１の場合、（ｄ）はその変形例２の場合である。 受信周波数制御を行わない場合の、（ａ）は金属が存在するときの周波数分布、（ｂ）は各Ｂｉｎ画像の信号、（ｃ）は各Ｂｉｎ画像の合成後の信号、を示す。 受信周波数制御を行った場合の、（ａ）は金属が存在するとき周波数分布、（ｂ）は各Ｂｉｎ画像の信号、（ｃ）は各Ｂｉｎ画像の合成後の信号、を示す。 ファントム撮像の画像を示す図である。 第二実施形態の制御部（ＣＰＵ）の構成を示す機能ブロック図。 第二実施形態の撮像の手順を示すフロー図。 撮像条件設定のためのＧＵＩの一例。 第三実施形態の信号処理部（ＣＰＵ）の構成を示す機能ブロック図。 第三実施形態の補正処理を説明する図で、（ａ）は計測毎の受信周波数を示す図、（ｂ）は第一実施形態の制御をする場合の受信周波数（目標受信周波数）を示す図である。

以下、添付図面に従って本発明のＭＲＩ装置の好ましい実施形態について詳述する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。

最初に、本発明が適用されるＭＲＩ装置の一例の全体概要を図１に基づいて説明する。図１は、本発明に係るＭＲＩ装置の一実施形態の全体構成を示すブロック図である。このＭＲＩ装置は、ＮＭＲ現象を利用して被検体の断層画像を得るもので、静磁場発生装置２と、傾斜磁場発生部３と、送信部５と、受信部６と、信号処理部７と、シーケンサ４と、中央処理装置（ＣＰＵ）８とを備えている。

静磁場発生装置２は、永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源（静磁場発生コイル）を備え、発生する静磁場の方向によって、被検体１のその体軸と直交する方向に磁場を発生する垂直磁場方式と体軸方向に磁場を発生する水平磁場方式がある。

傾斜磁場発生部３は、ＭＲＩ装置の座標系（静止座標系）であるＸ、Ｙ、Ｚの３軸方向に傾斜磁場を印加する傾斜磁場コイル９と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源１０とから成り、後述のシ−ケンサ４からの命令に従ってそれぞれのコイルの傾斜磁場電源１０を駆動することにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを印加する。これら３軸方向の傾斜磁場を組み合わせることで、任意の方向の傾斜磁場を印加することができる。例えば、スライス面（撮影断面）に直交する方向にスライス方向傾斜磁場パルス（Ｇｓ）を印加して被検体１に対するスライス面（または３次元計測時はスラブ面）を設定し、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード方向傾斜磁場パルス（Ｇｐ）と周波数エンコード方向傾斜磁場パルス（Ｇｆ）とを印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。３次元計測時はスライス方向に傾斜磁場パルス（Ｇｓ）を印加してスラブ面を設定するとともに、エンコード傾斜磁場としてスライスエンコード傾斜磁場パルスを追加する。

シーケンサ４は、高周波磁場パルス（以下、「ＲＦパルス」という）と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、ＣＰＵ８の制御で動作し、被検体１の断層画像のデータ収集に必要な種々の命令を送信部５、傾斜磁場発生部３、および受信部６に送る。

送信部５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体２にＲＦパルスを照射するもので、高周波発振器１１と変調器１２と高周波増幅器１３と送信側の高周波コイル（送信コイル）１４ａとを備える。高周波発振器１１から出力されたＲＦパルスをシーケンサ４からの指令によるタイミングで変調器１２により振幅変調し、この振幅変調されたＲＦパルスを高周波増幅器１３で増幅した後に被検体１に近接して配置された高周波コイル１４ａに供給することにより、ＲＦパルスが被検体１に照射される。

受信部６は、被検体１の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出するもので、受信側の高周波コイル（受信コイル）１４ｂと信号増幅器１５と直交位相検波器１６とＡ／Ｄ変換器１７とを備える。送信コイル１４ａから照射された電磁波によって誘起された被検体1の応答のＮＭＲ信号が被検体１に近接して配置された受信コイル１４ｂで検出され、信号増幅器１５で増幅された後、シーケンサ４からの指令によるタイミングで直交位相検波器１６により直交する二系統の信号に分割され、それぞれがとＡ／Ｄ変換器１７でディジタル量に変換されて、信号処理部７に送られる。なお直交位相検波器１６は、高周波発振器１１から送られる周波数信号を参照周波数として信号を検出する。即ち参照周波数により受信周波数が決まる。

信号処理部７は、各種データ処理と処理結果の表示及び保存等を行うもので、その機能の一部はＣＰＵ８により実行される。信号処理部７は、さらに光ディスク１９、磁気ディスク１８等の外部記憶装置と、ディスプレイ２０とを有する。受信部６からのデータがＣＰＵ８に入力されると、ＣＰＵ８が信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層画像をディスプレイ２０に表示すると共に、外部記憶装置の磁気ディスク１８等に記録する。

ＣＰＵ８は、シーケンサ４を介したＭＲＩ装置の各部を制御する制御部として機能するとともに上述した信号処理部７としての機能を実行する。このため、各種制御情報や上記信号処理部７で行う処理の制御情報を入力するための操作部２５を備える。操作部２５は、トラックボール又はマウス２３、及び、キーボード２４から成る。この操作部２５はディスプレイ２０に近接して配置され、操作者がディスプレイ２０を見ながら操作部２５を通してインタラクティブにＭＲＩ装置の各種処理を制御する。以下、シーケンサ４及びＣＰＵ８の制御機能を併せて制御部という。

制御部は、ＲＦパルスや傾斜磁場パルス及びエコー信号の計測のタイミング等を予め設定したパルスシーケンスと、操作部２５を介して操作者が設定する撮像パラメータとにより、撮像に用いるパルスシーケンス（撮像シーケンス）を計算する。パルスシーケンスは、撮像方法や検査プロトコル等によって予め定められたものが記憶装置等に格納されており、制御部は記憶装置等に格納されたパルスシーケンスから、操作者が設定したパルスシーケンスと撮像パラメータを読出し、撮像シーケンスを決定する。

本実施形態では、撮像シーケンスとして、スライス選択傾斜磁場を印加しない３次元計測のパルスシーケンスを実行する。３次元計測は、ＲＦパルスの周波数（中心周波数）を異ならせて複数回実行し、各計測で得られる３次元画像（Ｂｉｎ画像）を合成し、１つの画像を得る。３Ｄ−ＳＳＦＰ（Ｓｔｅａｄｙ−ｓｔａｔｅ ｆｒｅｅ ｐｒｅｃｅｓｓｉｏｎ）のパルスシーケンスの一例を図２に示す。ただし、パルスシーケンスは、図示するものに限定されず、スピンエコー系のパルスシーケンスでもよいし、また１回の繰り返し（ＴＲ）で１つのエコーを計測するパルスシーケンスでも、複数のエコーを計測するマルチエコーシーケンスでもよい。

このような撮像シーケンスを用いた計測において、制御部は、計測毎にＲＦパルスの送信周波数と受信部６の受信周波数を設定する。具体的には、送信周波数は、中心周波数が所定の周波数ｆ、例えば静磁場中心周波数、であるＢｉｎ（周波数帯域）の計測を中心としたとき、その両側に所定の周波数差Δｆずつ、ずれた複数のＢｉｎを設定する。受信周波数は、送信周波数が所定の周波数ｆである計測の受信周波数と、送信周波数が所定の周波数ｆと所定の周波数差Δｆである計測の受信周波数との差が０または最小となるように制御する。

所定の周波数ｆは、例えば画像において最も支配的なスピンの共鳴周波数であり、この周波数を送信及び受信周波数とする計測に対し、近傍の周波数の計測の受信周波数差を最小（例えば０）にすることにより、各計測で得た画像を合成する際に、周波数方向のずれを低減することができる。

以上説明したＭＲＩ装置の構成と、ＭＲＩ装置が採用する撮像シーケンスの概要を踏まえ、送信及び受信周波数制御の具体的な実施形態を説明する。

＜第一実施形態＞
本実施形態は、所定の送信周波数を中心に、それより高い周波数側及び低い周波数側に複数の周波数帯域（Ｂｉｎ）の計測を行い、各計測の送信周波数の周波数差を一定とする。所定の送信周波数は、静磁場中心周波数とする。受信周波数は、所定の送信周波数の計測では所定周波数とし、それより高い或いは低いＢｉｎの計測では、いずれも所定の周波数差分、静磁場中心周波数側にずらした周波数とする。

まず本実施形態を実現するための制御部（ＣＰＵ８）の構成例を図３に示す。図示するように、本実施形態の制御部８０は、計測制御部８１、送受信制御部８３及び表示制御部８５を含む。計測制御部８１は撮像シーケンスの計算を行い、シーケンサ４を介して送信部５、傾斜磁場発生部３及び受信部６の動作を制御する。送受信制御部８３は、高周波発振器１１を制御し、送信部５の高周波コイル１４ａが照射するＲＦパルスの中心周波数を、計測毎に所定の周波数差Δｆで変化させる。また受信時に直交位相検波器１６が参照する高周波発振器１１の周波数、即ち受信周波数を計測毎に制御する。表示制御部８５は、画像やＧＵＩをディスプレイ２０に表示させるための制御を行うとともにＧＵＩを介して操作者が入力した情報を受け付ける。
また信号処理部７は、画像演算に関わるＣＰＵの機能（画像演算部７０）として、画像再構成部７１、画像合成部７３及びその他の演算部を備える。

次に本実施形態における撮像及び制御の手順を、図４を参照して説明する。
撮像が開始されると、まず、制御部８０（表示制御部８５）は操作者による撮像条件や撮像パラメータを設定するためのＧＵＩをディスプレイ２０に表示し、撮像条件の設定を受け付ける（Ｓ４０１）。

撮像条件設定のためのＧＵＩの一例を図５に示す。図中、左側の列５１とその右側のボックス５２は撮像パラメータ名称とその値を入力するボックスで、ボタン５３で撮像を開始し、ボタン５４で撮像を中止する。パラメータのうち、「ＴＲ」はパルスシーケンスの繰り返し時間、「Ｓｌｉｃｅ＃」はスライスエンコード数を、「Ｐｈａｓｅ＃」は位相エンコード数を、「Ｆｒｅｑ＃」は周波数エンコード数を、「Ｓｐｅｃｔｒａｌ ｂｉｎ＃」はＢｉｎ数（計測数）を、「ＭｕｌｔｉＡｃｑ＃」は計測の分割数を、「ＢａｎｄＷｉｄｔｈ」は受信帯域幅（単位ｋＨｚ）を示す。「ＢｌｕｒＣｏｒｒｅｃｔ」は、ボケ補正のＯＮ／ＯＦＦの選択を受け付ける。上述した本実施形態による受信周波数制御は「ＢｌｕｒＣｏｒｒｅｃｔ」（ボケ補正）の「ＯＮ」が選択されたときに実行される。操作者は、例えば、観察したい部位や組織と被検体内に埋め込まれた金属との関係を考慮して、ボケ補正のＯＮ／ＯＦＦを選択する。

なお図５中の四角内に示す各パラメータの数値は単なる例示にすぎないが、以下の説明において各パラメータは図５に記載の値で計測するものとして説明に用いる。また図５に示されていないが、本実施形態の撮像、即ち複数のＢｉｎの計測、において、各Ｂｉｎの周波数差（周波数方向の間隔）Δｆはデフォルトで設定されているものとする。周波数差Δｆは、合成後の画像に周波数帯域の谷間を発生させない程度に狭いことが好ましく、例えば、１ｋＨｚ（１０００Ｈｚ）とする。

撮像パラメータが設定され、「Ｓｔａｒｔ」ボタン５３により撮像が開始されると、計測制御部８１はシーケンサ４を起動し、予め設定された撮像シーケンスにより撮像を開始する。この撮像シーケンスは、前述したようにスライス傾斜磁場を印加せずに、位相エンコード傾斜磁場及びスライスエンコード傾斜磁場を用いた３次元計測のパルスシーケンスである。即ち、シーケンサ４は傾斜磁場発生部３にスライス選択傾斜磁場（Ｇｓ）を印加しない命令を送る。また送受信制御部８３は、送信部５にＲＦパルスの周波数ＲＦ_Ｆｒｅｑを式（１）で求まる値とする命令を送る（Ｓ４０２）。

ただし、式（１）においてＣｅｎｔｅｒＦは静磁場中心周波数とする。またｉはＢｉｎ画像の番号（計測の番号と同じ）を示し、図５に示すようにＢｉｎ数が２０の場合、ｉは−９から１０までの値とする。即ち、ｉ＝０の送信周波数ＲＦ_ＦｒｅｑはＣｅｎｔｅｒＦと等しく、送信周波数をＢｉｎごとに１０００Ｈｚ間隔とする。

次いで送受信制御部８３は受信周波数を設定する。ここで送受信制御部８３は、分岐ステップＳ４０３において、「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮかＯＦＦかにより条件分岐し、「ＯＦＦ」の場合には、ステップＳ４０４の処理を、「ＯＮ」の場合にはステップＳ４０５の処理を実行する。

ステップＳ４０４の処理は、各計測において受信周波数を送信周波数に合わせる処理であり、受信部６に受信周波数帯域の中心周波数ＡＤ Ｃｅｎｔｅｒを式（２）で求まる値とする命令を送る。すなわち、送信周波数と同様に、受信周波数をＢｉｎごとに１０００Ｈｚ間隔（周波数差）とする。

一方、ボケ補正ＯＮの場合には、送受信制御部８３は受信部６に中心周波数ＡＤ Ｃｅｎｔｅｒを式（３）で求まる値とする命令を送る。すなわち、受信周波数の中心をＢｉｎごとに１０００Ｈｚ間隔から、中心に向けて１０００Ｈｚずつ近づけた値とする（Ｓ４０５）。

この場合の送信周波数（中心周波数）と受信周波数（中心周波数）との関係を図６（ｂ）に示す。図６（ａ）は受信周波数の制御を行わない場合（「ＢｌｕｒＣｏｒｒｅｃｔ」ＯＦＦ）の場合である。図６では、「ＣｅｎｔｅｒＦ」を「ｆ０」と示している。図示するように、この処理では、送信周波数を静磁場中心周波数に合わせたＢｉｎの計測（ｉ＝０）に対し、周波数方向両側に位置するＢｉｎの計測（ｉ＝−１、ｉ＝１）では、その受信周波数は、中央の計測の受信周波数と同じに制御される。

シーケンサ４はステップＳ４０４又はステップＳ４０５の受信周波数条件で撮像を行い、エコー信号を受信し、信号処理部７に送る（Ｓ４０６）。

設定した全ての計測（２０の計測）でエコー信号の受信が完了するまで、上述したステップを繰り返す（分岐ステップＳ４０７）。この繰り返しの繰り返し時間ＧｒｏｕｐＴＲは、次式（４）で求められる時間となる。

式（４）は、「ＧｒｏｕｐＴＲ」はパルスシーケンスの繰り返し時間ＴＲ内に、Ｂｉｎ画像の計測を「ＭｕｌｔｉＡｃｑ＃」で設定された数（図５の例では４）だけ分割して計測することを示している。

計測毎にエコー信号の受信が完了すると、画像再構成部７１がフーリエ変換等の演算により画像再構成し３Ｄ画像（Ｂｉｎ画像）を得る。全ての計測でエコー信号受信が完了した後、画像合成部７３が全てのＢｉｎ画像を合成し、撮像が完了する（Ｓ４０８）。なおＢｉｎの計測順序は、特に限定されず、ｉ＝−９から順にｉ＝１０まで計測してもよいし、各Ｂｉｎについてｋ空間の低周波領域のエコー信号を取得した後、各Ｂｉｎの高周波領域のエコー信号を取得するなど任意に変更することが可能である。また必要に応じて、体動モニタリングやナビゲーションエコーなどを採用し、合成前のＢｉｎ画像に体動補正などを行ってもよい。

本実施形態によれば、静磁場中心周波数を中心周波数とするＢｉｎ画像に対し周波数方向に隣接するＢｉｎ画像の受信周波数を静磁場中心周波数側に揃えることにより、画像において支配的なスピンの画像上の位置が周波数方向にずれるのを防止し、ボケの発生が抑制された合成画像を得ることができる。

以下、本実施形態により、周波数エンコード方向の位置ずれが解消されることを詳述する。

本実施形態の撮像において、周波数エンコード方向の位置ｘにあり共鳴周波数ＳｐｉｎＦｒｅｑをもつスピンは、Ｂｉｎ画像において式（５）で求まる位置Ｘに描出される。

ΔＸは、共鳴周波数ＳｐｉｎＦｒｅｑと受信の中心周波数ＡＤ Ｃｅｎｔｅｒとの差で決まる周波数方向のずれである。

ここで、ＲＦパルスの周波数帯域幅（ＢａｎｄＷｉｄｔｈ）の半値幅を２２５０Ｈｚとすると、１つのスピンは隣り合う３つのＢｉｎ画像に描出される。それらのＢｉｎ画像においてΔＸが１より小さい場合は同じ位置に描出されるが、１より大きい場合は異なる位置に描出される。静磁場中心周波数付近のスピンに着目し、ＳｐｉｎＦｒｅｑを静磁場中心周波数ＣｅｎｔｅｒＦ、周波数エンコード数Ｆｒｅｑ＃を４００、帯域幅ＢａｎｄＷｉｄｔｈを２００とすると、「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＦＦの場合は、
ΔＸ［０］＝０
ΔＸ［１］＝２
ΔＸ［−１］＝２
となるため、これを合成するとボケが発生する。

一方、「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮの場合は、
ΔＸ［０］＝ΔＸ［１］＝ΔＸ［−１］＝０
となるため、これを合成するとボケが発生しない。

また一様な水の中に金属がある状態を例に、本実施形態による金属周辺の描出について、図７および図８を用いて説明する。これら図面は、横軸は周波数エンコード方向の位置を示し、（ａ）は金属から生じる磁場による共鳴周波数変化の分布を、（ｂ）は各Ｂｉｎ画像にて描出される信号分布を、（ｃ）は合成後の信号分布を示す。両図に示すように、金属から生じる磁場が金属からの距離の二乗に反比例し、静磁場中心周波数付近の領域が広いことを示している。

図７は「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＦＦの場合、図８は「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮの場合であり、ＯＦＦの場合を破線で示している。図８からわかるように、「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮの場合は、ＯＦＦの場合に描出される位置より周波数エンコード方向にシフトする。しかし、Ｂｉｎ画像の信号の周波数は、ｉ＝０およびｉ＝１では広い範囲で分布するため、周波数エンコード方向のシフト量に対して、合成後の信号形状の変化が少なく、またｉ＝０以外のＢｉｎ画像において一様にシフトするため、合成後の信号は形状が保たれたままシフトする。つまり、計測の目的である金属アーチファクト軽減の効果が「ＢｌｕｒＣｏｒｒｅｃｔ」のＯＮ、ＯＦＦに関わらず同様であることを示す。

本実施形態によるファントム撮像の画像を図９に示す。図９中、画像９１は金属とプラスチック格子をプラスチック容器内に設置した写真であり、撮像時には容器に水を満たして撮像した。画像９２は「ＢｌｕｒＣｏｒｒｅｃｔ」＝ＯＦＦの場合の撮像結果であり、画像９３は「ＢｌｕｒＣｏｒｒｅｃｔ」＝ＯＮの場合の撮像結果である。尚、周波数エンコード方向はどちらも画像の上下方向である。画像９３は、画像９２よりもプラスチック格子のボケが改善されていることがわかる。金属周辺が周波数方向にシフトしているが、金属周辺のアーチファクトの描出は同程度であった。

以上、説明したように、本実施形態によれば、操作者に金属に起因するボケの防止を行うか否かの選択を行わせる手段が提供される。また本実施形態によれば、送信周波数とは別に受信周波数を制御することで、静磁場中心周波数付近のスピンにボケが発生せず、且つ金属アーチファクト軽減の効果が従来と同様の画像が得られる撮像が可能となる。

＜変形例１＞
上述した実施形態では、中心周波数が静磁場中心周波数であるＢｉｎ（基準）に対し周波数方向に両側の複数のＢｉｎの計測において、順次、受信周波数を静磁場中心周波数側にシフトする制御を行ったが、受信周波数のシフトは静磁場中心周波数であるＢｉｎの両側の波数全てのＢｉｎの計測で行う必要はない。本変形例では、隣接するＢｉｎの受信周波数のみを制御する。その他の構成は第一実施形態と同様であり、以下、異なる点を中心に説明する。

制御部８０による制御の手順は、第一実施形態と同様に図４に示すフローに従い行われる。但し、ステップＳ４０５において、送受信制御部８３は受信部６に、受信周波数ＡＤＣｅｎｔｅｒを式（６）で求まる値とする命令を送る。

すなわち、受信周波数の中心を、中心（ｉ＝０）の近傍（ｉ＝1、又は−１）の２つＢｉｎだけ中心に向けて１０００Ｈｚずつ近づけた値（ここではｉ＝０のときの中心周波数と同じ値）とする。本変形例における送信周波数と受信周波数との関係を、図６（ｃ）に示す。

本変形例においても、第一実施形態と同様に、ボケの発生が抑制できる。すなわち、本変形例で制御される受信周波数を、周波数エンコード方向のずれ量ΔＸを求める前掲の式（５）に当てはめると、静磁場中心周波数付近のスピンに着目して、ＳｐｉｎＦｒｅｑを中心周波数ＣｅｎｔｅｒＦ、周波数エンコード数Ｆｒｅｑ＃を４００、帯域幅を２００とすると、「ＢｌｕｒＣｏｒｒｅｃｔ」ＯＮのときの受信周波数では、
ΔＸ［０］＝ΔＸ［１］＝ΔＸ［−１］＝０
となるため、これら計測で得たＢｉｎ画像を合成するとボケが発生しない。

金属周辺の描出については、第一実施形態では、中央のＢｉｎ以外では中心周波数が静磁場中心周波数側に１０００Ｈｚずつずれているため金属周辺が周波数方向にシフトしていたが、本変形例では、中心近傍の２つのＢｉｎ画像以外は、送信周波数と同一に保たれるため、金属周辺における周波数エンコード方向のシフトは発生しない。但し、中心近傍の２つのＢｉｎでは画像がシフトするため、１０００Ｈｚから２０００Ｈｚの位置および−１０００Ｈｚから−２０００Ｈｚの位置において合成後の信号形状に段差が発生する。

＜変形例２＞
変形例１では、第一実施形態に対し、金属周辺の周波数エンコード方向へのシフトが発生しないように受信周波数に変更を加えたものであるが、本変形例では、さらに変形例１で発生した合成後の信号形状の段差を抑制するように受信周波数を制御する。

具体的には、中心のＢｉｎ以外のＢｉｎについて、周波数方向に隣接するＢｉｎとの周波数差が、中心（ｉ＝０）に近いＢｉｎほど１０００Ｈｚに近く、遠いＢｉｎほど０Ｈｚに近い値とする。

このため、図４に示す制御フローのステップＳ４０５において、送受信制御部８３は、受信部６に、受信周波数ＡＤＣｅｎｔｅｒを式（７）で求まる値とする命令を送る。

本変形例における送信周波数と受信周波数との関係を、図６（ｄ）に示す。

本変形例で制御される受信周波数を、周波数エンコード方向のずれ量ΔＸを求める前掲の式（５）に当てはめると、静磁場中心周波数付近のスピンに着目して、ＳｐｉｎＦｒｅｑを中心周波数ＣｅｎｔｅｒＦ、周波数エンコード数Ｆｒｅｑ＃を４００、帯域幅を２００とすると、「ＢｌｕｒＣｏｒｒｅｃｔ」ＯＮのときの受信周波数では、
ΔＸ［０］＝０、ΔＸ［１］＝−０．０２、ΔＸ［−１］＝０．０２
となるため、これを合成するとほぼボケが発生しない。

また、金属周辺の描出について、「ＢｌｕｒＣｏｒｒｅｃｔ」ＯＮの場合は、ＯＦＦの場合に描出される位置より周波数エンコード方向にシフトするが、金属に近づくほどシフト量が小さくなるため、金属の位置や大きさが正しく描出される。ただし、Ｂｉｎごとにシフト量が異なるため、金属に近づくほど合成後の信号形状に段差が発生する。

以上、第一実施形態の２つの変形例を説明したが、これら変形例においても、第一実施形態と同様に、静磁場中心周波数付近のスピンにボケが発生せず、且つ且つ金属アーチファクトの軽減効果が得られる。また変形例１では金属周辺の周波数エンコード方向へのシフトが発生しない画像が得られる。さらに変形例１、変形例２では、Ｂｉｎ間の受信周波数の間隔の相違によって、変形例１では中心周波数近傍に、変形例２では金属に近づくほど信号形状に段差が発生する。従って、「ＢｌｕｒＣｏｒｒｅｃｔ」ＯＮの際に、観察しようとする部位に応じて、第一実施形態、変形例１及び変形例２の制御のいずれかを操作者が選択できるようにＧＵＩを設定してもよい。

＜第二実施形態＞
第一実施形態では、送信周波数の中心周波数が静磁場中心周波数であるＢｉｎの計測を中心として、その両側に複数のＢｉｎを設定し、それらの受信周波数を制御したが、本実施形態は、被検体を構成するスピンのうち最も信号強度の高いスピンの周波数を中心周波数とするＢｉｎを中心として、その両側の設定した複数のＢｉｎの受信周波数を制御する。信号強度が最も高い周波数は、例えば、プリスキャンを実行することで求めることができる。

本実施形態の制御部８０及び信号処理部７（画像演算部７０）の構成を図１０に示す。図１０において、図３と同様の要素は同じ符号で示す。本実施形態では、図３の構成に対し、信号処理部７に最大信号強度の周波数を算出する周波数算出部７５が追加されていることが異なる。また計測制御部８１は、第一実施形態の撮像シーケンスの他に、プリスキャンのシーケンスを実行するようにプログラムされている。

上記構成を踏まえ、本実施形態における撮像及び制御の手順を、図１１を参照して説明する。

撮像が開始されると、まず、制御部８０（表示制御部８５）は操作者による撮像条件や撮像パラメータを設定するためのＧＵＩをディスプレイ２０に表示し、撮像条件の設定を受け付ける（Ｓ１４０１）。

本実施形態では、図１２に示すようなＧＵＩが表示される。図１２に示す各パラメータは図５に示すパラメータと同じであるが、図１２ではさらに「ＢｌｕｒＣｏｒｒｅｃｔ」が選択されると、「Ｐｒｅｓｃａｎ」を実行するか否か（ＹＥＳ、ＮＯ）の選択ブロック５５が表示される。「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮであって「Ｐｒｅｓｃａｎ」が実行しない場合は、第一実施形態及びその変形例のいずれかの手順となる。

ステップＳ１４０１で「Ｐｒｅｓｃａｎ」が選択されている場合には（Ｓ１４０２）、撮像制御部８１はシーケンサ４に指令を送り、プリスキャンを実行する（Ｓ１４０３）。プリスキャンでは、スライス方向傾斜磁場パルス（Ｇｓ）、位相エンコード方向傾斜磁場パルス（Ｇｐ）、周波数エンコード方向傾斜磁場パルス（Ｇｆ）およびスライスエンコード傾斜磁場パルスを印加せず、式（１）で示す照射周波数にて、設定されたＢｉｎ数繰り返しＮＭＲ信号を取得する。

信号処理部７はＮＭＲ信号を受け取ると、周波数算出部７５はＮＭＲ信号をフーリエ変換し、周波数スペクトルを取得する。そして、得られた周波数スペクトルから最も信号の強い周波数（最大強度周波数という）を求める（Ｓ１４０４）。最大強度周波数ＦｒｅｑＭａｘ（静磁場中心周波数との差）は、受信周波数のＢｉｎ間の周波数差（間隔）の算出に用いられる。

送受信制御部８３は、第一実施形態のステップＳ４０２と同様に、中心となるＢｉｎの送信周波数とＢｉｎ間の周波数差を決定する（Ｓ１４０５）。次に、分岐処理Ｓ１４０６において、ステップＳ１０１で「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮとＯＦＦのいずれに設定されたかを判断し、ＯＦＦの場合は第一実施形態のステップＳ４０４と同様で、所定の周波数差、即ち送信周波数の周波数差と同じで且つ各Ｂｉｎの送信周波数＝受信周波数となるように受信周波数を設定する（Ｓ１４０７）。

一方「ＢｌｕｒＣｏｒｒｅｃｔ」がＯＮの場合は、受信周波数ＡＤＣｅｎｔｅｒを次式（８）で求まる値とする。

式（８）において、Ｍｏｄは、引数１つ目を引数２つ目で割った余りを返す関数である。すなわち、受信周波数の中心を、最大強度周波数ＦｒｅｑＭａｘの余りだけずらし、且つＢｉｎごとに最大強度周波数ＦｒｅｑＭａｘに向けて１０００Ｈｚずつ近づけた値とする。

本実施形態で制御される受信周波数を、周波数エンコード方向のずれ量ΔＸを求める前掲の式（５）に当てはめるとＦｒｅｑＭａｘ付近のスピンに着目して、ＳｐｉｎＦｒｅｑをＦｒｅｑＭａｘ、Ｆｒｅｑ＃を４００、ＢａｎｄＷｉｄｔｈを２００とすると、ＢｌｕｒＣｏｒｒｅｃｔがＯＮの場合は、受信周波数の中心がＦｒｅｑＭａｘとなるＢｉｎをｊ−1、ｊ、ｊ＋１とすると
ΔＸ［ｊ］＝ΔＸ［ｊ＋１］＝ΔＸ［ｊ−１］＝０
となるため、これを合成するとＦｒｅｑＭａｘ付近にボケが発生しない。即ち、最も画像において支配的なスピンが合成前の各Ｂｉｎ画像において周波数方向にシフトしないため、合成後の画像のボケが抑制される。

ただし、ＢｌｕｒＣｏｒｒｅｃｔがＯＮの場合は、ＯＦＦの場合に描出される位置より周波数エンコード方向にシフトする。またＢｉｎによってシフト方向が正負反転するため、合成後の信号形状が変化する。

本実施形態は、特に、被検者の体格に対して大きな金属がある場合に有効である。このような場合、ほぼ全てのスピンが金属から発生する磁場の影響を受けて静磁場中心周波数付近から外れる状況がある。このような状況では、静磁場中心周波数付近のＢｉｎ画像を合成することにより発生するボケよりも、静磁場中心周波数から離れた周波数帯のＢｉｎ画像を合成することにより発生するボケの方が、画像全体に与える影響が大きくなる。本実施形態では画像全体に与える影響が最も大きい周波数ＦｒｅｑＭａｘを中心に受信周波数を制御するので、画像全体に対する影響の大きいボケを解消することができ、画質を向上することができる。

＜第三実施形態＞
第一実施形態とその変形例及び第二実施形態では、いずれも、画像に影響を与える周波数を中心に受信周波数の中心周波数を制御したが、本実施形態では、受信周波数を制御するのではなく、Ｂｉｎ画像を受信周波数の間隔に応じて周波数方向に補正する手段を備えることが特徴である。

本実施形態の制御部及び信号処理部の構成を図１３に示す。図示するように、本実施形態の信号処理部７（画像演算部７０）は、合成前のＢｉｎ画像を補正する補正部７７を備える。補正部７７は、画像に影響を与える周波数（例えば静磁場中心周波数）を中心周波数とするＢｉｎに対し、その近傍のＢｉｎ画像の周波数エンコード方向の位置をシフトする補正を行う。

送受信制御部８３は、各Ｂｉｎ画像の計測において、受信周波数を送信周波数の周波数差（間隔）と同じに制御する。即ち静磁場中心周波数に対し、受信周波数の中心周波数を１０００Ｈｚずつシフトさせて、各Ｂｉｎの計測を行う。この様子を図１４に示す。

こうして取得したＢｉｎ画像は、第一実施形態により制御した受信周波数で撮像する場合と比較すると、図１４（ｂ）に示すように、ＣｅｎｔｅｒＦのＢｉｎ画像の両側のＢｉｎでは、１０００Ｈｚの周波数の差がある。この周波数差Δｆは、Ｂｉｎ画像では周波数方向の位置ずれΔＸと、式（５）に示す関係がある。即ち
ΔＸ（ｉ）＝Δｆ（ｉ）×Ｆｒｅｑ＃／（Ｂａｎｄｗｉｄｔｈ×１０００）

そこで補正部７７は、ｉ番目のＢｉｎ画像に対し、各画素のＸ座標を
Ｘ（ｉ）＝ｘ−ΔＸ（ｉ）
とする補正を行う。

これにより第一実施形態における受信周波数制御と同様の効果を得ることができる。なお以上の説明では、第一実施形態の受信周波数制御と同様の補正をする場合を説明したが、これは一例であり、第一実施形態の変形例１，２の受信周波数制御についても、制御すべき目標周波数と実際に設定した受信周波数との差Δｆから同様の補正を行うことができる。また以上の説明では、送受信制御部８３が、受信周波数を送信周波数と同じ周波数差（間隔、１０００Ｈｚ）で制御する場合を説明したが、受信周波数は静磁場中心周波数に固定しておいもよく、その場合にも、制御すべき目標周波数との差を用いて同様の補正が可能である。

以上、本発明のＭＲＩ装置の各実施形態を説明したが、本発明は少なくとも画像に影響を与える周波数の近傍において、周波数方向の位置ずれを防止する受信周波数の制御を行うことが特徴であり、上記各実施形態に限定されることなく、種々の変更が可能である。例えば、また受信周波数の制御を行う手段として、ＭＲＩ装置の制御系であるＣＰＵを例示したが、このＣＰＵとは別に制御部を設けてもよいし、ＣＰＵ以外の構成要素で実現することも可能である。また実施形態で説明したＣＰＵで行われる各種機能や演算の一部は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）やＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）などのハードウェアで実現することも可能である。さらに、ＭＲＩ装置を構成する要素の一部を省略したり、或いは別の要素を追加したものも本発明に包含される。

２：静磁場発生装置、３：傾斜磁場発生部、４：シーケンサ、５：送信部、６：受信部、７：信号処理部（画像演算部）、８：ＣＰＵ、７１：画像再構成部、７３：画像合成部、７５：周波数算出部、７７：補正部、８１：計測制御部、８３：送受信制御部、８５：表示制御部。

Claims (10)

1. 静磁場発生磁石と、
   前記静磁場発生磁石が形成する静磁場空間に、傾斜磁場を発生する傾斜磁場コイルと、
   前記静磁場空間に置かれた被検体に、高周波磁場を印加する送信部と、
   前記被検体から発生される核磁気共鳴信号を受信する受信部と、
   前記高周波磁場の送信周波数を予め定めた周波数差で変化させて、複数の計測を行う制御部と、
   前記複数の計測で取得した複数の前記被検体の画像を合成する画像演算部と、を備え、
   前記複数の計測のうち、基準周波数の計測に対し、周波数方向に隣接する周波数の計測において、前記基準周波数の計測との受信周波数差を０又は最小とするように、前記受信周波数を制御する送受信制御部又は前記複数の計測で得た画像を補正する補正部をさらに備えることを特徴とする磁気共鳴イメージング装置。
2. 請求項１に記載の磁気共鳴イメージング装置であって、
   前記制御部は、前記受信部の受信周波数を制御する送受信制御部を備えることを特徴とする磁気共鳴イメージング装置。
3. 請求項２に記載の磁気共鳴イメージング装置であって、
   前記送受信制御部は、前記基準周波数の計測に対し、周波数方向に隣接する計測の受信周波数の中心周波数を前記基準周波数と同じ周波数とし、その他の計測は、所定の周波数差で前記受信周波数を変化させることを特徴とする磁気共鳴イメージング装置。
4. 請求項２に記載の磁気共鳴イメージング装置であって、
   前記送受信制御部は、前記基準周波数の計測に対し、周波数方向に隣接する計測の受信周波数の中心周波数を前記基準周波数と同じ周波数とし、その他の計測は、受信周波数の中心周波数を送信周波数の中心周波数に一致させることを特徴とする磁気共鳴イメージング装置。
5. 請求項２に記載の磁気共鳴イメージング装置であって、
   前記送受信制御部は、前記基準周波数の計測に対し、その他の計測は、送信周波数と同一の周波数を上限として、前記基準周波数から周波数方向に離れるにつれて、計測間の周波数差が大きくなるように、各計測の受信周波数を制御することを特徴とする磁気共鳴イメージング装置。
6. 請求項１に記載の磁気共鳴イメージング装置であって、
   前記基準周波数は、静磁場中心における核磁気共鳴周波数であることを特徴とする磁気共鳴イメージング装置。
7. 請求項１に記載の磁気共鳴イメージング装置であって、
   前記基準周波数は、前記核磁気共鳴信号の周波数スペクトルにおいて信号値が最大となる周波数であることを特徴とする磁気共鳴イメージング装置。
8. 請求項７に記載の磁気共鳴イメージング装置であって、
   前記制御部は、前記周波数スペクトルを、送信周波数を異ならせた複数の予備計測でそれぞれ取得した核磁気共鳴信号から取得することを特徴とする磁気共鳴イメージング装置。
9. 請求項１に記載の磁気共鳴イメージング装置であって
   前記画像演算部は、画像の合成に用いる画像を補正する補正部を備え、
   前記補正部は、計測毎に、前記基準周波数と受信周波数との周波数差を用いて、前記画像の周波数エンコード方向の補正量を算出することを特徴とする磁気共鳴イメージング装置。
10. 磁気共鳴イメージング装置を制御し、複数の異なる送信周波数を用いて複数の計測を行い、当該複数の計測により得た複数の三次元画像を合成して１つの画像を取得する磁気共鳴イメージング装置の制御方法であって、
    複数の撮像の際の受信周波数を前記送信周波数の変化に応じて変化させる制御を行い、その際、基準周波数の計測に対し周波数方向に隣接する計測の受信周波数を、前記基準周波数との差が０又は最小となるように各計測における前記受信周波数を制御することを特徴とする磁気共鳴イメージング装置の制御方法。

Images