(57)【要約】【課題】パラレルイメージングに用いるコイルの感度分布の精度を向上可能であり、画質を向上可能なコイルの感度分布作成方法、磁気共鳴撮影方法および磁気共鳴撮影装置を提供する。【解決手段】複数の受信コイルＣ１，Ｃ２のそれぞれについて、各コイルが検出した磁気共鳴信号をフーリエ変換して位相のデータθＣ１，θＣ２を算出して出力し、位相のデータθＣ１，θＣ２と所定の一つのコイルから得られる位相のデータとの差分を計算し位相差分データθ’’Ｃ１，θ’’Ｃ２を出力し、位相差分データθ’’Ｃ１，θ’’Ｃ２に基づいてフィッティング処理を施し、それぞれのコイルＣ１，Ｃ２について、被検体の領域と被検体の領域以外の領域とを含んで感度領域の位相データ分布θ’Ｃ１，θ’Ｃ２を算出して出力し、磁気共鳴信号の大きさのデータ分布Ｍ’Ｃ１，Ｍ’Ｃ２とそれぞれ合成して感度分布Ｓ１，Ｓ２を生成する。【選択図】  図５
2004344327.tif 000002
【特許請求の範囲】
【請求項１】
感度領域内の被検体からの磁気共鳴信号を検出するコイルの感度分布を作成するコイルの感度分布作成方法であって、
複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを入手するステップと、
前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分である位相差分データをそれぞれ算出するステップと、
前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出するステップと
を含むコイルの感度分布作成方法。
【請求項２】
前記磁気共鳴信号の分析において、前記位相のデータに加えて大きさのデータを入手し、
前記フィッティング処理において、前記位相差分データのフィッティングに加えて前記大きさのデータに基づく大きさ処理データをフィッティングし、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の大きさ処理データ分布を算出し、
前記位相データ分布と前記大きさ処理データ分布とに基づいて、前記複数のコイルのそれぞれについて前記感度分布を作成する
請求項１に記載のコイルの感度分布作成方法。
【請求項３】
前記複数のコイルのそれぞれの前記位相差分データに、所定のオフセット値を加える
請求項１または２に記載のコイルの感度分布作成方法。
【請求項４】
前記所定の一つのコイルを、全ての前記コイルのうちから検出する前記磁気共鳴信号の強度に基づいて選択する
請求項１〜３のいずれかに記載のコイルの感度分布作成方法。
【請求項５】
前記所定の一つのコイルとして、全ての前記コイルのうち前記感度領域が最も広いコイルを用いる
請求項１〜３のいずれかに記載のコイルの感度分布作成方法。
【請求項６】
コイルの感度領域内の被検体からの磁気共鳴信号の位相方向のエンコードステップを間引いて得られる折返し偽像が存在する画像から、複数のコイルの前記感度領域における感度分布の差に基づいて前記折返し偽像を除去する磁気共鳴撮影方法であって、
前記感度分布を、
前記複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを入手するステップと、
前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分である位相差分データをそれぞれ算出するステップと、
前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出するステップと
を含んで生成する
磁気共鳴撮影方法。
【請求項７】
前記磁気共鳴信号の分析において、前記位相のデータに加えて大きさのデータを入手し、
前記フィッティング処理において、前記位相差分データのフィッティングに加えて前記大きさのデータに基づく大きさ処理データをフィッティングし、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の大きさ処理データ分布を算出し、
前記位相データ分布と前記大きさ処理データ分布とに基づいて、前記複数のコイルのそれぞれについて前記感度分布を作成する
請求項６に記載の磁気共鳴撮影方法。
【請求項８】
前記複数のコイルのそれぞれの前記位相差分データに、所定のオフセット値を加える
請求項６または７に記載の磁気共鳴撮影方法。
【請求項９】
前記所定の一つのコイルを、全ての前記コイルのうちから検出する前記磁気共鳴信号の強度に基づいて選択する
請求項６〜８のいずれかに記載の磁気共鳴撮影方法。
【請求項１０】
前記所定の一つのコイルとして、全ての前記コイルのうち前記感度領域が最も広いコイルを用いる
請求項６〜８のいずれかに記載の磁気共鳴撮影方法。
【請求項１１】
コイルの感度領域内の被検体からの磁気共鳴信号を位相方向のエンコードステップを間引いて検出して生成した折返し偽像が存在する画像から、複数のコイルの前記感度領域における感度分布の差に基づいて前記折返し偽像を除去する磁気共鳴撮影装置であって、
前記複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを算出して出力する分析部と、
前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分を計算し位相差分データを出力する位相差分計算部と、
前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出して出力するフィッティング処理部と
を有する磁気共鳴撮影装置。
【請求項１２】
前記分析部は、前記磁気共鳴信号の分析において前記位相のデータに加えて大きさのデータを算出して出力し、
前記フィッティング処理部は、前記位相差分データのフィッティングに加えて前記分析部からの前記大きさのデータに基づく大きさ処理データをフィッティングし、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の大きさ処理データ分布を算出し、
前記フィッティング処理部からの前記位相データ分布と前記大きさ処理データ分布とを合成して複数の前記コイルのそれぞれについて前記感度分布を作成する合成部をさらに有する
請求項１１に記載の磁気共鳴撮影装置。
【請求項１３】
前記位相差分計算部は、前記複数のコイルのそれぞれの前記位相差分データに、所定のオフセット値を加える
請求項１１または１２に記載の磁気共鳴撮影装置。
【請求項１４】
前記位相差分計算部は、全ての前記コイルのうちから検出する前記磁気共鳴信号の強度に基づいて選択された前記所定の一つのコイルからの前記位相のデータを用いて前記位相差分データを計算する
請求項１１〜１３のいずれかに記載の磁気共鳴撮影装置。
【請求項１５】
前記位相差分計算部は、全ての前記コイルのうち前記感度領域が最も広いコイルを前記所定の一つのコイルとして前記位相差分データを計算する
請求項１１〜１３のいずれかに記載の磁気共鳴撮影装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、コイルによって検出した被検体の被検部位からの磁気共鳴信号に基づいて被検部位の画像を生成する磁気共鳴撮影に関する。特定的には、本発明は磁気共鳴信号を間引いて検出して生成した折返し偽像の存在する画像から折返し偽像を取り除くために用いるコイルの感度分布の作成方法、この作成方法により作成した感度分布を用いた磁気共鳴撮影方法および磁気共鳴撮影装置に関する。
【０００２】
【従来の技術】
磁気共鳴撮影（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ：ＭＲＩ）において、パラレルイメージング（Ｐａｒａｌｌｅｌ  Ｉｍａｇｉｎｇ）と呼ばれる撮影方法が知られている（たとえば、非特許文献１参照）。
パラレルイメージング法においては、位相方向のエンコードステップを間引いて被検体からの磁気共鳴信号を検出することにより撮影時間を短縮して撮影を行ない、撮影視野が狭まり折返り偽像（アーチファクト）が存在する画像を生成する。そして、複数のコイルの感度分布差に基づいて折返り偽像を取り除き、撮影視野が広がった画像を最終的に得る。
また、パラレルイメージング法を用いた磁気共鳴撮影装置も知られている（たとえば、特許文献１参照）。
【０００３】
【特許文献１】
特開２００２−２４８０８９号公報
【非特許文献１】
ケー・ピー・プリュースマン（Ｋ．Ｐ．Ｐｒｕｅｓｓｍａｎｎ），  エム・ヴァイガー（Ｍ．Ｗｅｉｇｅｒ），  エム・ビー・シャイデガー（Ｍ．Ｂ．Ｓｃｈｅｉｄｅｇｇｅｒ），  ピー・ベージガー（Ｐ．Ｂｏｅｓｉｇｅｒ），  マグネティック・レゾナンス・イン・メディシン（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  ｉｎ  Ｍｅｄｉｃｉｎｅ），１９９９年，４２，  ９５２
【０００４】
【発明が解決しようとする課題】
パラレルイメージング法においては磁気共鳴信号を検出する複数のコイルの検出感度分布の差を用いて折返り偽像を除去する。検出感度分布を単に感度分布と言うこともあり、感度分布の差を単に感度分布差と言うこともある。しかしながら、コイルの感度分布を正確に、またコイルの感度領域において全体的に得ることは困難であった。このため、折返し偽像を取り除いて最終的に得られる画像の画質は、理論上得られる画質よりも低かった。
【０００５】
本発明の目的は、パラレルイメージングに用いるコイルの感度分布の精度を向上可能なコイルの感度分布作成方法を提供することにある。
また、コイルの感度分布の精度向上に伴い、画質を向上させることが可能な磁気共鳴撮影方法および磁気共鳴撮影装置を提供することも本発明の目的である。
【０００６】
【課題を解決するための手段】
本発明に係るコイルの感度分布作成方法は、感度領域内の被検体からの磁気共鳴信号を検出するコイルの感度分布を作成するコイルの感度分布作成方法であって、複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを入手するステップと、前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分である位相差分データをそれぞれ算出するステップと、前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出するステップとを含む。
【０００７】
本発明に係る磁気共鳴撮影方法は、コイルの感度領域内の被検体からの磁気共鳴信号の位相方向のエンコードステップを間引いて得られる折返し偽像が存在する画像から、複数のコイルの前記感度領域における感度分布の差に基づいて前記折返し偽像を除去する磁気共鳴撮影方法であって、前記感度分布を、前記複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを入手するステップと、前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分である位相差分データをそれぞれ算出するステップと、前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出するステップとを含んで生成する。
【０００８】
本発明に係る磁気共鳴撮影装置は、コイルの感度領域内の被検体からの磁気共鳴信号を位相方向のエンコードステップを間引いて検出して生成した折返し偽像が存在する画像から、複数のコイルの前記感度領域における感度分布の差に基づいて前記折返し偽像を除去する磁気共鳴撮影装置であって、前記複数のコイルのそれぞれについて、各コイルが検出した前記磁気共鳴信号を分析して位相のデータを算出して出力する分析部と、前記複数のコイルのそれぞれについて、当該コイルから得られる前記位相のデータと所定の一つのコイルから得られる前記位相のデータとの差分を計算し位相差分データを出力する位相差分計算部と、前記位相差分データに基づいてフィッティング処理を施し、それぞれの前記コイルについて、前記被検体の領域と当該被検体の領域以外の領域とを含んで前記感度領域の位相データ分布を算出して出力するフィッティング処理部とを有する。
【０００９】
本発明においては、折返し偽像の除去に用いるコイルの感度分布を算出するために、コイルが被検体から検出した磁気共鳴信号を分析部により分析する。分析部は、複数のコイルのそれぞれについて、各コイルが検出した磁気共鳴信号を分析して磁気共鳴信号に関する位相のデータを算出する。位相差分計算部は、分析部が算出した各コイルに関する位相のデータと所定の一つのコイルに関する位相のデータとの差分を計算し、位相差分データを出力する。フィッティング処理部は、各コイルに関する位相差分データに基づいてフィッティング処理を施し、各コイルについて位相データ分布を算出する。フィッティング処理部は、各コイルの感度領域における被検体の領域と被検体の領域以外の領域とを含んで、各コイルの位相データ分布を計算する。
以上のようにして求められた位相データ分布を利用して、各コイルの感度分布が最終的に生成される。複数のコイルの感度領域における感度分布の差に基づいて、折返し偽像が除去される。
【００１０】
【発明の実施の形態】
以下、本発明の実施の形態について、添付図面を参照しながら述べる。
【００１１】
図１は本発明の一実施の形態に係る磁気共鳴撮影（ＭＲＩ）装置の概略構成を示すブロック図であり、図２は図１に示すＭＲＩ装置に用いるマグネットアセンブリの一例の外観を示す斜視図である。
【００１２】
図１に示すように、本実施の形態に係るＭＲＩ装置１００は、マグネットアセンブリ１、主磁場電源２、勾配磁場駆動回路３、ＲＦ電力増幅器４、前置増幅器５、計算機７、シーケンス記憶回路８、ゲート変調回路９、ＲＦ発振回路１０、Ａ／Ｄ変換器１１、位相検波器１２、操作コンソール１３、表示装置１４および記憶部３０を有する。
【００１３】
マグネットアセンブリ１とＭＲＩ装置１００のその他の構成要素とは、適宜離して配置される。
本実施の形態に係るマグネットアセンブリ１はいわゆるシリンドリカルタイプと呼ばれるものである。シリンドリカルタイプのマグネットアセンブリ１は、図２に示すように内部にボア（Ｂｏｒｅ）１３１ａを有しており、このボア１３１ａを取巻くように図示しない主磁場発生部、勾配磁場発生部および送信コイルを配置した構成を有している。
【００１４】
本実施の形態においては、図２に示すようにボア１３１ａの形状を円筒形にしている。ボア１３１ａは、その内部に撮影のために被検体４を収容し、ボア１３１ａの内部が撮影空間となる。
被検体４は、テーブル１３３によって支持されたクレードル１３２上に載置される。クレードル１３２がテーブル１３３上においてボア１３１ａの長手方向に沿ってボア１３１ａ側へ移動することによって、被検体４がボア１３１ａの内部に収容される。
以下では、円筒形のボア１３１ａの長手方向に沿う方向をＺ方向とする。また、Ｚ方向に直交する水平方向をＸ方向、鉛直方向をＹ方向とする。
【００１５】
主磁場発生部は、たとえば、常電導磁石を用いて構成する。主磁場発生部は主磁場電源２に接続され、主磁場電源２からの供給電流に応じた磁界強度の均一な主磁場をボア１３１ａ内に形成する。
主磁場発生部が形成する主磁場は、被検体４の撮影中常に磁界強度が一定でボア１３１内において均一となるように保たれるため、静磁場とも呼ばれる。なお、マグネットアセンブリ１における静磁場の方向はＺ方向に沿う方向である。
【００１６】
主磁場発生部は、常電導磁石の他に超伝導磁石や永久磁石等の他の磁石を用いて構成することもできる。永久磁石を用いる場合には主磁場電源２は必要ない。
【００１７】
勾配磁場発生部は、後述する受信コイルＣ１，Ｃ２が検出する磁気共鳴信号に３次元の位置情報を持たせるために３系統のコイルを有する。
勾配磁場発生部は勾配磁場駆動回路３に接続され、勾配磁場駆動回路３からの駆動信号に応じた勾配磁場を発生する。
【００１８】
送信コイルは、静磁場および勾配磁場が印加される被検体４の被検部位に対して回転磁場を印加するためのものである。回転磁場は、被検部位の原子核のスピンをその中心軸を傾けて回転させるためのＲＦ（Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ）帯の磁場である。送信コイルはＲＦ電力増幅器４に接続され、ＲＦ電力増幅器４から供給される電力に応じた回転磁場を被検部位に印加する。
【００１９】
また、ＭＲＩ装置１００は、静磁場および勾配磁場中の被検部位への回転磁場の印加によるエコーとして発生する被検部位からの磁気共鳴信号を検出するための受信コイルを有する。受信コイルとしては、たとえば、被検体４の頭部や腹部や肩等の被検部位を部分的に覆うような専用のコイルを用いる。
図２においては、受信コイルの一例として、ＸＺ平面に平行に対向配置してフェーズドアレイ化した矩形状のコイルＣ１，Ｃ２を挙げている。コイルＣ１とコイルＣ２との間に被検体４が挟まれる。
受信コイルＣ１，Ｃ２はそれぞれ前置増幅器５に接続され、受信コイルＣ１，Ｃ２が検出した磁気共鳴信号が、前置増幅器５によって増幅される。
【００２０】
シーケンス記憶回路８は、勾配磁場駆動回路３と、ゲート変調回路９と、Ａ／Ｄ変換器１１と、計算機７の後述する制御部２１とに接続されている。
ゲート変調回路９は、ＲＦ電力増幅器４とＲＦ発振回路１０とにさらに接続されている。ＲＦ発振回路１０は、前記増幅器５およびＡ／Ｄ変換器１１に接続される位相検波器１２にさらに接続される。
【００２１】
シーケンス記憶回路８は、制御部２１が生成した磁気共鳴画像取得のためのパルスシーケンスを記憶し、制御部２１から撮影指令信号が入力されたときには記憶したパルスシーケンスに応じて勾配磁場駆動回路３、ゲート変調回路９およびＡ／Ｄ変換器１１を駆動制御する。
【００２２】
勾配磁場駆動回路３は、シーケンス記憶回路８からの制御に応じて、所定のパルス波形で勾配磁場発生部を駆動する駆動信号を勾配磁場発生部に出力する。
ＲＦ発振回路１０は、所定周波数のＲＦ信号を生成してゲート変調回路９および位相検波器１２に出力する。ゲート変調回路９は、シーケンス記憶回路８からの制御に基づいて所定のタイミングでＲＦ発振回路１０からのＲＦ信号を変調し、所定の波形のＲＦ励起パルス信号を、ＲＦ電力増幅器４を介して送信コイルに印加する。
【００２３】
位相検波器１２は、ＲＦ発振回路１０から出力されるＲＦ信号を参照信号とし、前置増幅器５からの出力信号（受信コイルＣ１，Ｃ２によって検出され、前置増幅器５により増幅された磁気共鳴信号）を位相検波したのちにＡ／Ｄ変換器１１に出力する。Ａ／Ｄ変換器１１は、シーケンス記憶回路８からの制御に基づいて、位相検波器１２による位相検波後のアナログ信号をデジタル信号に変換し、計算機７の画像処理部２３に出力する。Ａ／Ｄ変換器１１から画像処理部２３に出力されるデジタル信号のデータを、ローデータ（Ｒａｗ  Ｄａｔａ）と呼ぶ。
【００２４】
計算機７は、図１に示すように構成要素として制御部２１と画像処理部２３と記憶部３０とを有する。制御部２１および画像処理部２３はそれぞれ異なるハードウェアを用いて実現してもよいし、１つのハードウェアを兼用させてもよい。
Ａ／Ｄ変換器１１に接続される画像処理部２３が制御部２１に接続される。
また、制御部２１および画像処理部２３は、たとえばＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）やハードディスクドライブによって実現される記憶部３０に適宜アクセス可能になっている。
さらに、制御部２１には、たとえばキーボードやマウスによって実現され、オペレータの操作に応じた信号を制御部２１に出力するための操作コンソール１３と、テレビモニター等の表示装置１４とが接続される。
【００２５】
画像処理部２３は、Ａ／Ｄ変換器１１から送信されたデジタルデータに対して以下に詳述するパラレルイメージング法のアンフォールディング（Ｕｎｆｏｌｄｉｎｇ）処理を施し、折返し偽像（アーチファクト）の無い磁気共鳴画像を生成する。画像処理部２３が生成した磁気共鳴画像は記憶部３０に記憶することが可能である。
【００２６】
制御部２１は、操作コンソール１３を介してオペレータから指令が入力されたときには、パラレルイメージング法によって磁気共鳴画像を取得するように勾配磁場駆動回路３、ゲート変調回路９およびＡ／Ｄ変換器１１の駆動状態を規定するパルスシーケンスを生成し、シーケンス記憶回路８に送信する。
なお、パルスシーケンスは、勾配磁場駆動回路３、ゲート変調回路９およびＡ／Ｄ変換器１１に対する駆動パルス信号の波形のデータを意味する。
また、制御部２１は、画像処理部２３が生成した画像および記憶部３０に記憶させていた画像を適宜表示装置１４に表示させる。
【００２７】
ここで、パラレルイメージング法について、図３を参照しながら簡単に述べる。なお、パラレルイメージング法はＳＥＮＳＥ（Ｓｅｎｓｉｔｉｖｉｔｙ  Ｅｎｃｏｄｉｎｇ）法とも呼ばれ、その詳細は、たとえば、前出の非特許文献１に記載されている。
図３には、フェーズドアレイ化した２個のコイルＣ１，Ｃ２をそれぞれ用いて２枚の断層像を撮影し、この２枚の断層像とコイルＣ１，Ｃ２のそれぞれの感度分布を用いて１枚の最終的な断層像を生成する例を挙げている。
【００２８】
１）パラレルイメージング法を適用した撮影においては、まず、最終的に得たい画像ＤＶの撮影視野（Ｆｉｅｌｄ  ｏｆ  Ｖｉｅｗ：ＦＯＶ）の大きさを決める。このとき定めた大きさのＦＯＶをＦｕｌｌ  ＦＯＶと呼ぶ。図３においては、Ｙ方向においてＦＯＶの大きさを決めた例を挙げている。図３に示す矩形状の画像ＤＶはＹ方向におけるＦｕｌｌ  ＦＯＶ画像となる。ＦＯＶの大きさは操作コンソール１３を介して制御部２１に入力する。なお、図３に示すＸ，Ｙ，Ｚの各方向は、図２におけるＸ，Ｙ，Ｚ方向に対応しており、以下ではＸＹ平面に平行な断層像の撮影を例として挙げる。また、Ｘ方向が磁気共鳴信号の位置情報を規定する方向のうちの周波数エンコーディング方向、Ｙ方向が位相エンコーディング方向、Ｚ方向がスライス方向とする。これらの３方向は、前述の勾配磁場発生部の３系統のコイルの各コイルが発生する勾配磁場の方向に対応している。
【００２９】
２）ＦＯＶを規定した後には、２個の受信コイルＣ１，Ｃ２をそれぞれ用いて、各コイルＣ１，Ｃ２の感度分布を入手するための撮影を行なう。受信コイルの感度分布を得るための撮影を、キャリブレーションスキャンと言う。
キャリブレーションスキャンにおいては、位相エンコーディング方向においてエンコードステップを間引かずにスキャンを行なう。
制御部は、キャリブレーションスキャンのためのパルスシーケンスを生成してシーケンス記憶回路８に送信し、キャリブレーションスキャンを実行させる。
【００３０】
図３に示す画像ＤＳ１およびＤＳ２は、コイルＣ１，Ｃ２の感度分布を示すために、ＸＹ平面に平行な断面において楕円形のファントムを撮影したときの画像である。キャリブレーションスキャンにおいてはエンコードステップを間引かずに撮影するため、画像ＤＳ１およびＤＳ２はＦｕｌｌ  ＦＯＶ画像となる。画像ＤＳ１およびＤＳ２において、白く表示されている領域ほどコイルＣ１，Ｃ２の感度が高く、検出した磁気共鳴信号の強度が高いことを示している。
矩形状のコイルＣ１，Ｃ２はＸＺ平面に平行に配置される。図３に示すＹ１側に配置されるコイルＣ１を用いて撮影された画像ＤＳ１においてはＹ１側の信号強度が高く、Ｙ２側に配置されるコイルＣ２を用いて撮影された画像ＤＳ２においてはＹ２側の信号強度が高いことが分かる。
このように、２個のコイルＣ１，Ｃ２の感度分布には差が存在する。
【００３１】
３）キャリブレーションスキャンの終了後に、コイルＣ１，Ｃ２を用いて位相エンコーディング方向においてエンコードステップを間引きながら、１）のステップの断面と同じ断面において撮影を行なう。この撮影を本スキャンという。
本スキャンのためのパルスシーケンスとしては、たとえば、エコープラナーイメージング（Ｅｃｈｏ  Ｐｌａｎａｒ  Ｉｍａｇｉｎｇ：ＥＰＩ）法やスピンエコー法等の各種パルスシーケンスを適宜用いることができる。制御部２１は、オペレータから指定された本スキャンのためのパルスシーケンスを生成してシーケンス記憶回路８に送信し、実行させる。
図３に示す画像ＤＡ１がコイルＣ１を用いて得られた本スキャンの画像であり、画像ＤＡ２がコイルＣ２を用いて得られた本スキャンの画像である。本スキャンにおいて得られる画像ＤＡ１，ＤＡ２のＦＯＶは、エンコードステップを間引くことに起因して位相エンコーディング方向に小さくなっている。このＦＯＶを、Ｓｍａｌｌ  ＦＯＶと呼ぶ。
【００３２】
図３においては、一例として位相エンコードステップを１／２に間引く例を挙げる。位相エンコードステップ数を１／２にしたときには、Ｓｍａｌｌ  ＦＯＶの大きさはＦｕｌｌ  ＦＯＶの大きさの１／２になる。また、本スキャンのスキャン時間は、エンコードステップを間引かない時の１／２になる。
Ｓｍａｌｌ  ＦＯＶの画像には、エンコードステップを間引くことに起因する折返し偽像が存在する。たとえば、画像ＤＡ１，ＤＡ２に示すように、被検体の領域のうちＳｍａｌｌ  ＦＯＶの領域をはみ出た部分が、折返された偽像Ａｆとなり被検体の画像ＩＢ上に重なって表示される。
【００３３】
４）Ｓｍａｌｌ  ＦＯＶの画像ＤＡ１，ＤＡ２に対して２）のステップにおいて得られたコイルＣ１，Ｃ２のそれぞれの感度分布を用いてアンフォールディング処理を施し、折返し偽像の無いＦｕｌｌ  ＦＯＶの画像ＤＶを得る。これにより、エンコードステップを１／２に間引いたときには、半分のスキャン時間（撮影）でＦｕｌｌ  ＦＯＶ画像ＤＶを得ることができる。
画像処理部２３によって実行されるアンフォールディング処理の内容については後述する。
【００３４】
上述のように、パラレルイメージング法においてはコイルＣ１，Ｃ２の感度分布差を利用してＳｍａｌｌ  ＦＯＶ画像から折返し偽像を取り除く。
図４は、Ｓｍａｌｌ  ＦＯＶ画像ＤＡ１，ＤＡ２とＦｕｌｌ  ＦＯＶ画像ＤＶとを用いて、コイルの感度分布について述べるための図であり、（ａ）がＳｍａｌｌ  ＦＯＶ画像ＤＡ１，ＤＡ２を、（ｂ）がＦｕｌｌ  ＦＯＶ画像ＤＶをそれぞれ示している。
【００３５】
たとえば、図４（ａ）の画像ＤＡ１，ＤＡ２におけるピクセルＰＸ１をアンフォールディング処理することを考える。画像ＤＡ１，ＤＡ２にはＳｍａｌｌ  ＦＯＶの領域外の被検体の画像が折返されて存在しているため、ピクセルＰＸ１をアンフォールディング処理するためには、Ｆｕｌｌ  ＦＯＶ画像におけるピクセルＰＸ２とピクセルＰＸ３の情報が必要となる。
コイルＣ１，Ｃ２が画像を生成可能な感度領域ＶＳのうち、被検体が存在する領域に相当する画像ＩＢ内のピクセルＰＸ３については、被検体４から磁気共鳴信号を入手することが可能なためコイルＣ１，Ｃ２の感度を特定することができる。しかしながら、画像ＩＢの領域以外の領域ＶＴ、すなわちキャリブレーションスキャンにおいて被検体４が存在していなかった領域のピクセルＰＸ２については、被検体４が存在していなかったためピクセルＰＸ２に対応する位置からの磁気共鳴信号を検出できなかったことになる。このため、ピクセルＰＸ２におけるコイルＣ１，Ｃ２の感度を直接的に特定することはできず、ノイズを用いてコイルＣ１，Ｃ２の感度分布を作成することになる。
【００３６】
したがって、本実施の形態においては、コイルＣ１，Ｃ２の検出信号に基づいてフィッティング（ｆｉｔｔｉｎｇ）処理を行ない、感度領域ＶＳ全体における感度分布を作成する。
図５は、本実施の形態における感度分布作成、および作成した感度分布を用いた磁気共鳴撮影における処理の流れを示すブロック図である。
また、図６は、図１に示すＭＲＩ装置１００における計算機７の画像処理部２３の機能ブロック図である。
【００３７】
図６に示すように、画像処理部２３は、フーリエ変換部２４と、規格化部２５と、位相差分計算部２６と、フィッティング処理部２７と、合成部２８と、アンフォールド処理部２９とを有する。
本発明における分析部の一実施態様がフーリエ変換部２４に相当する。
【００３８】
本実施の形態においては、受信コイルＣ１，Ｃ２によるキャリブレーションスキャンの検出信号に基づくローデータを分析して、ローデータから検出信号の位相のデータと大きさ（ｍａｇｎｉｔｕｄｅ）のデータとを取出す。ローデータの分析の一手法として、フーリエ変換部２４はＡ／Ｄ変換器１１から入力されたローデータをフーリエ変換し、位相のデータと大きさのデータを算出する。
フーリエ変換部２４は、算出した位相のデータは位相差分計算部２６に出力し、大きさのデータは規格化部２５に出力する。
また、フーリエ変換部２４は、本スキャンによるローデータにフーリエ変換を施し、図３に示すようなＳｍａｌｌ  ＦＯＶの画像ＤＡ１，ＤＡ２の画像データを生成し、アンフォールド処理部２９に出力する。
【００３９】
位相差分計算部２６は、フーリエ変換部２４から入力された、受信コイル毎の位相のデータに対し、所定の一つの受信コイルから得られた位相データとの差分を計算する。
たとえば、コイルＣ１，Ｃ２のそれぞれの位相のデータに対して、コイルＣ１の位相のデータとの差分を計算する。
位相差分計算部２６は、計算により得られた位相差分データをフィッティング処理部２７に出力する。
【００４０】
規格化部２５は、フーリエ変換部２４から入力された受信コイル毎の検出信号の大きさのデータに対し、コイル毎の特性の差を解消するための規格化（ｎｏｒｍａｌｉｚｉｎｇ）を施す。規格化部２５は、規格化後の規格化データをフィッティング処理部２７に出力する。
本発明における大きさ処理データの一実施態様が規格化データである。
【００４１】
フィッティング処理部２７は、位相差分計算部２６および規格化部２５からそれぞれ入力された受信コイル毎の位相差分データおよび規格化データにそれぞれフィッティング処理を施し、各受信コイルの感度領域全体における位相差分データ分布および規格化データ分布を算出する。
フィッティング処理部２７は、算出した位相差分データ分布および規格化データ分布の情報を合成部２８に出力する。
【００４２】
合成部２８は、フィッティング処理部２７から入力された位相差分データ分布および規格化データ分布の情報を合成し、コイルＣ１，Ｃ２ごとに図４（ｂ）に示すような感度領域ＶＳ全域における感度分布を作成する。
合成部２８は、作成した各コイルの感度分布の情報を、アンフォールド処理部２９に出力する。
【００４３】
アンフォールド処理部２９は、フーリエ変換部２４から送信されたＳｍａｌｌ  ＦＯＶ画像ＤＡ１，ＤＡ２の画像データに対して合成部２８から送信されたコイルＣ１，Ｃ２の感度分布の情報を用いてアンフォールディング処理を施し、Ｆｕｌｌ  ＦＯＶの画像ＤＶを生成する。
【００４４】
以下に、図５に沿ってさらに詳細に述べる。
まず、フーリエ変換部２４には、ブロック５１，５３に示すようなＳｍａｌｌ  ＦＯＶ画像ＤＡ１，ＤＡ２の元となるローデータＲＡ１，ＲＡ２が入力される。また、フーリエ変換部２４にはブロック５５，５７に示すようなコイルＣ１，Ｃ２のそれぞれに対応したキャリブレーションスキャンによるローデータＲＣ１，ＲＣ２も入力される。
【００４５】
ブロック６０に示すように、フーリエ変換部２４は、入力されたローデータＲＡ１，ＲＡ２，ＲＣ１，ＲＣ２に対してそれぞれ２次元フーリエ変換（２ＤＦＴ）を施す。
ブロック６１，６２に示すように、ローデータＲＡ１，ＲＡ２からはＳｍａｌｌ  ＦＯＶ画像ＤＡ１，ＤＡ２の画像データＡ１，Ａ２がそれぞれ生成される。また、ブロック６５，６７に示すように、ローデータＲＣ１，ＲＣ２からはコイルＣ１，Ｃ２に対応したフィッティングを施す前の感度分布を表わす画像ＩＣ１，ＩＣ２の画像データがそれぞれ生成される。
前述のように、フーリエ変換部２４はフーリエ変換によりローデータから位相のデータおよび大きさのデータを取出す。したがって、ブロック６１，６２，６５，６７における画像ＤＡ１，ＤＡ２，ＩＣ１，ＩＣ２を生成するデータは、それぞれ位相のデータと大きさのデータを有する複素数のデータとなる。
【００４６】
ブロック７１，７５にそれぞれ示すように、画像ＩＣ１の複素数のデータのうち、大きさのデータをＭＣ１、位相のデータをθＣ１とする。
また、ブロック７３，７７にそれぞれ示すように、画像ＩＣ２の複素数のデータのうち、大きさのデータをＭＣ２、位相のデータをθＣ２とする。
【００４７】
図７は、位相のデータθＣ１およびθＣ２に対するフィッティング処理の概念をイメージとして捉えるための図である。
図７（ａ），（ｄ）は、位相データθＣ１，θＣ２をそれぞれ用いて生成した、図３，４と同じＸＹ平面に平行な断面における画像ＩθＣ１，ＩθＣ２をそれぞれ示している。
画像ＩθＣ１，ＩθＣ２におけるパターンＤＰＣ１，ＤＰＣ２はそれぞれ、位相データθＣ１，θＣ２の分布をイメージとして表わすためのものである。
前述のように磁気共鳴信号は被検体４が存在する領域からのみ検出されるため、パターンＤＰＣ１，ＤＰＣ２は被検体の画像ＩＢに相当する領域にのみ現われる。
【００４８】
また、コイルＣ１，Ｃ２による検出信号には、ＭＲＩ装置１００全体としての位相情報やボア１３１ａにおける磁場の不均一等の外乱成分が混入する。画像ＩθＣ１，ＩθＣ２において、この外乱成分のイメージを、外乱成分パターンＵＭＧとして描いている。外乱成分はコイルＣ１，Ｃ２の両方にほぼ等しく混入するため、外乱成分パターンＵＭＧは画像ＩθＣ１，ＩθＣ２においてほぼ同じ位置に現われる。
【００４９】
本実施の形態においては、ブロック８３に示すように、コイルＣ１，Ｃ２の位相データθＣ１，θＣ２に対し、所定の一つの受信コイルの位相データとの差分を位相差分計算部２６によって計算し、ブロック９５，９７に示すように位相差分データθ’’Ｃ１，θ’’Ｃ２をそれぞれ算出する。
このとき、本実施の形態に係る位相差分計算部２６は、所定のオフセット値Ｃｓｔを加えて位相差分データθ’’Ｃ１，θ’’Ｃ２を計算する。オフセット値Ｃｓｔは、たとえば、所定の大きさの定数である。オフセット値Ｃｓｔを加えることにより、重み付け等の処理が可能になる。
【００５０】
たとえば、コイルＣ１の位相データθＣ１を用いて位相差分データを計算する場合には、θ’’Ｃ１＝θＣ１−θＣ１＋Ｃｓｔ＝Ｃｓｔとなり、θ’’Ｃ２＝θＣ２−θＣ１＋Ｃｓｔとなる。
図７（ｂ），（ｅ）は、このようにして求めた位相差分データθ’’Ｃ１，θ’’Ｃ２を用いて生成した画像Ｉθ’’Ｃ１，Ｉθ’’Ｃ２をそれぞれ示している。画像Ｉθ’’Ｃ１は、オフセット値Ｃｓｔに基づく均一な画像が被検体画像ＩＢの領域に分布しているように表示される。
また、画像Ｉθ’’Ｃ２における被検体の画像ＩＢの領域には、位相差分データθ’’Ｃ２に基づくパターンＤＤＰＣ２が現われる。
【００５１】
受信コイルがＮ個存在する場合には、Ｎ個の受信コイルのｎ番目（１≦ｎ≦Ｎ）の受信コイルの位相差分データθ’’ｎは、ｎ番目のコイルの位相データθｎと差分計算に用いるｌ番目のコイルの位相データθｌとを用いて、θ’’ｎ＝θｎ−θｌ＋Ｃｓｔと表わされる。ｌは位相の差分をＮ個の受信コイルについて計算するときには固定であるが、１からＮの間で任意に選択することができる。
【００５２】
フィッティング処理部２７が、ブロック１００に示すように位相差分データθ’’Ｃ１，θ’’Ｃ２に対してそれぞれフィッティング処理を施す。
フィッティング処理部２７は、たとえば、以下に示す２次元多項式によって位相差分データθ’’Ｃ１，θ’’Ｃ２をフィッティングし、ブロック１０５，１０７に示すように位相データ分布θ’Ｃ１，θ’Ｃ２を得る。
【００５３】
【数１】
θ’ｎ  ＝Ａｎ・ｘ２＋Ｂｎ・ｙ２＋Ｃｎ・ｘ＋Ｄｎ・ｙ＋Ｅｎ・ｘｙ＋Ｆｎ…（１）
【００５４】
上式（１）において、引き数ｎ（１≦ｎ≦Ｎ）はコイルＣ１，Ｃ２等のコイルを表わし、変数ｘ，ｙは図７に示すような、画像のｘ座標、ｙ座標をそれぞれ表わす。変数ｘ，ｙの各次数の係数が、位相差分データによって一意に決まる。
【００５５】
図７（ｃ）がコイルＣ１の位相データ分布θ’Ｃ１によって得られる画像Ｉθ’Ｃ１を表わしており、図７（ｆ）がコイルＣ２の位相データ分布θ’Ｃ２によって得られる画像Ｉθ’Ｃ２を表わしている。
フィッティングにより、感度領域ＶＳ全体において位相データの分布を得ることができる。画像Ｉθ’Ｃ１においては位相差分データθ’’Ｃ１が被検体画像ＩＢ以外の領域ＶＴにも外挿され、感度領域ＶＳ全体にオフセット値Ｃｓｔによる画像が均一に分布する画像となる。画像Ｉθ’Ｃ２は、図７（ｆ）に示すように、パターンＤＤＰＣ２が被検体画像ＩＢ以外の領域ＶＴを含む感度領域ＶＳ全体に広がった画像となる。
【００５６】
以上のような手順により、位相のデータについてはフィッティング処理が終了し、コイルＣ１，Ｃ２の感度領域ＶＳ全体における分布を得ることができる。
大きさのデータについても同様に、以下に示すようにフィッティングして感度領域ＶＳ全体における分布を計算する。
【００５７】
コイルＣ１，Ｃ２がそれぞれ検出した磁気共鳴信号の大きさのデータＭＣ１，ＭＣ２は、図５のブロック８１に示すように、規格化部２５によって規格化される。規格化後の大きさのデータＭＣ１，ＭＣ２を、ブロック９１，９３に示すようにそれぞれ規格化データＭ’’Ｃ１，Ｍ’’Ｃ２とする。
【００５８】
図８は、規格化データＭ’’Ｃ１，Ｍ’’Ｃ２に対するフィッティング処理の概念をイメージとして捉えるための図である。
図８（ａ），（ｃ）は、規格化データＭ’’Ｃ１，Ｍ’’Ｃ２をそれぞれ用いて生成した、図３，４，７と同じＸＹ平面に平行な断面における画像ＩＭＣ１，ＩＭＣ２をそれぞれ示している。
画像ＩＭＣ１，ＩＭＣ２においては、規格化データＭ’’Ｃ１，Ｍ’’Ｃ２の分布を表わすイメージであるパターンＤＭＣ１，ＤＭＣ２が、それぞれ被検体画像ＩＢの領域にのみ現われる。
【００５９】
ブロック１００に示すように、フィッティング処理部２７は規格化データＭ’’Ｃ１，Ｍ’’Ｃ２にフィッティング処理を施す。
フィッティング処理部２７は、たとえば、以下に示す２次元多項式によって規格化データＭ’’Ｃ１，Ｍ’’Ｃ２をフィッティングし、ブロック１０１，１０３に示すように規格化データ分布Ｍ’Ｃ１，Ｍ’Ｃ２を得る。
【００６０】
【数２】
Ｍ’ｎ  ＝ａｎ・ｘ２＋ｂｎ・ｙ２＋ｃｎ・ｘ＋ｄｎ・ｙ＋ｅｎ・ｘｙ＋ｆｎ…（２）
【００６１】
上式（２）における引き数ｎおよび変数ｘ，ｙの意味は式（１）の場合と同じである。変数ｘ，ｙの各次数の係数が、規格化データＭ’’Ｃ１，Ｍ’’Ｃ２によって一意に決まる。
【００６２】
図８（ｂ）がコイルＣ１の規格化データ分布Ｍ’Ｃ１によって得られる画像ＩＭ’Ｃ１を表わしており、図８（ｄ）がコイルＣ２の規格化データ分布Ｍ’Ｃ２によって得られる画像ＩＭ’Ｃ２を表わしている。
被検体画像ＩＢにおける規格化データＭ’’Ｃ１，Ｍ’’Ｃ２の分布と規格化データ分布Ｍ’Ｃ１，Ｍ’Ｃ２は完全には一致しない場合もあるが、基本的には、画像ＩＭ’Ｃ１，ＩＭ’Ｃ２は図８（ｂ），（ｄ）に示すようにパターンＤＭＣ１，ＤＭＣ２がそれぞれ感度領域ＶＳ全体に広がった画像となる。
【００６３】
以上により、コイルＣ１，Ｃ２のそれぞれについて感度領域ＶＳ全体における規格化データ分布と位相データ分布を得ることができる。
合成部２８は、ブロック１１１，１１３に示すように、下記式（３）により規格化データ分布Ｍ’Ｃ１と位相データ分布θ’Ｃ１とを合成してコイルＣ１の感度領域ＶＳ全体についての感度分布Ｓ１を生成し、規格化データ分布Ｍ’Ｃ２と位相データ分布θ’Ｃ２とを合成してコイルＣ２の感度分布Ｓ２を生成する。
【００６４】
【数３】
Ｓｎ＝Ｍ’ｎ  ｅ−ｉθ’ｎ…（３）
【００６５】
式（３）におけるθ’ｎは式（１）を表わし、Ｍ’ｎは式（２）を表わす。また、引き数ｎはコイルＣ１，Ｃ２等のコイルを表わす。ここでは、ｎ＝１のときコイルＣ１を表わし、ｎ＝２のときコイルＣ２を表わす。受信コイルがたとえばＮ個存在する場合には、１≦ｎ≦Ｎとなる。
また、図９（ａ），（ｂ）は、感度分布Ｓ１，Ｓ２をイメージングした画像ＤＳ１，ＤＳ２をそれぞれ示している。画像ＤＳ１，ＤＳ２におけるパターンＰＳ１，ＰＳ２は、コイルＣ１，Ｃ２の感度分布をそれぞれ示すパターンである。図９（ａ），（ｂ）に示すように、フィッティングにより感度領域ＶＳ全体における各コイルの感度分布を入手することができる。
【００６６】
ブロック１２０に示すように、アンフォールド処理部２９は、画像データＡ１，Ａ２および感度分布Ｓ１，Ｓ２を用いて、下記式（４）によりブロック１３０に示すようにＦｕｌｌ  ＦＯＶ画像ＤＶの画像データＶを計算して生成する。
【００６７】
【数４】
Ｖ（ｐｉｘｅｌ）＝（Ｓ（ｃｏｉｌ，ｐｉｘｅｌ）＊  Ｓ（ｃｏｉｌ，ｐｉｘｅｌ））−１  Ｓ（ｃｏｉｌ，ｐｉｘｅｌ）＊  Ａ（ｃｏｉｌ）…（４）
【００６８】
上式（４）における引き数ｃｏｉｌは複数の受信コイルのチャンネルを表わす。たとえば、本実施の形態においてはコイルＣ１またはコイルＣ２となる。したがって、たとえば、感度分布Ｓには上述の感度分布Ｓ１またはＳ２を、画像データＡには画像データＡ１またはＡ２をそれぞれ当てはめる。
また、引き数ｐｉｘｅｌは折返りの引き数であり、リダクション要素Ｒに相当する。Ｒ＝２であればエンコーディング数を１／２に間引いてスキャンし、引き数ｐｉｘｅｌは２までとなる。
なお、式（４）におけるＳ＊  は感度分布Ｓの行列の随伴（ｃｏｎｊｕｇａｔｅ  ｔｒａｎｓｐｏｓｅ）行列である。
画像データＶを用いることにより、図３および図４（ｂ）に示すようなＦｕｌｌ  ＦＯＶの画像ＤＶを生成することができる。
【００６９】
以上のように、本実施の形態においては、各受信コイルの受信信号の位相データについて、所定の一つのコイルの位相データとの差分を計算する。これにより、磁場不均一による位相特性やＭＲＩ装置全体の位相特性等の、各受信コイル特有の位相特性以外の外乱としての位相特性が相殺され、各受信コイルの位相特性を正確に入手することが可能になる。キャリブレーションスキャンにおいてグラディエントエコー法を用いる場合には磁場不均一の影響が大きいため、差分を計算することは特に効果的である。
また、位相データの差分についてフィッティングを施すため、フィッティングを容易に実行することができ、たとえば２次多項式等の比較的簡単な関数の形に位相データの分布を回帰させることができる。フィッティングにより、コイルの感度領域全体について位相特性を得ることができる。差分を計算したときにもコイル間における相対的位相関係は保持されるため、フォールディング処理に対する影響はない。
【００７０】
さらに、本実施の形態においては、位相データだけでなく大きさのデータについてもフィッティングを施し、フィッティング後の位相データと大きさのデータを用いて各受信コイルの感度分布を生成する。これにより、コイルの感度領域全体についての感度分布を得ることができる。位相については差分を計算することにより外乱の影響が取り除かれているため、感度分布の精度が向上する。
キャリブレーションスキャンと本スキャンとにおいて被検体４の位置がずれていたとしても、感度領域全体において感度分布を得ることができるため、アンフォールディング処理を確実に、かつ、ある程度の精度で実行することができる。その結果、被検体の位置ずれによるＦｕｌｌ  ＦＯＶ画像の画質劣化を抑制可能になる。また、フィッティングにより滑らかな感度分布が得られるため、感度分布を表わす画像のノイズが減少し、感度分布画像からＦｕｌｌ  ＦＯＶ画像へ伝播されるノイズが減少する。その結果、Ｆｕｌｌ  ＦＯＶ画像のノイズが減少してＳ／Ｎ比が向上し、Ｆｕｌｌ  ＦＯＶ画像の画質が向上する。
感度分布は画像処理部２３における比較的簡単な計算によって生成することができるため、本実施の形態に係るような感度領域全体において精度を向上可能な感度分布を得るためにＭＲＩ装置１００等のハードウェアに特別な変更を施す必要はない。
【００７１】
変形形態
位相差分計算部２６における差分計算に用いるｌ番目のコイルは、各コイルが受信する磁気共鳴信号の強度に基づいて選択するようにしてもよい。たとえば、最も強い磁気共鳴信号を検出したコイルを選択して使用する。または、１スキャンにおける検出信号の２次元積分を計算し、１スキャンにおける平均的な検出信号強度が最も高いコイルを選択して使用してもよい。
コイルの選択は、たとえば、入手したデータに基づいて位相差分計算部２６が自動的に判断する形態にすることによって実現することができる。もしくは、各コイルの受信信号の強度データを表示装置１４に表示させるようにし、表示されたデータに基づいてオペレータが指定するようにしてもよい。
信号強度が相対的に高いコイルを選択することによって、信号とノイズとの判別が明確となり、感度分布の精度を向上させることができる。
【００７２】
磁気共鳴信号の強度の代わりに、感度領域が最も広いコイルを選択するようにしてもよい。基本的に、感度領域はコイルにおける磁気共鳴信号を検出可能な検出面の面積に比例する。このため、差分に用いる受信コイルは、各受信コイルの形状に基づいて位相差分計算部２６に対して予め指定しておくことができる。
感度領域が最も広いコイルを基準とすることにより、位相データの差分の値を必ず計算することが可能になり、フィッティング処理およびアンフォールディング処理を確実に実行することができるようになる。
【００７３】
なお、上記実施の形態および図面に記載の内容は本発明を説明するための一例であり、上記実施の形態および図面に記載の形状や数値等の条件は特許請求の範囲内において適宜変更可能である。たとえば、上記実施の形態においては矩形状のコイルＣ１，Ｃ２を例に挙げて述べたが、受信コイルにはバードケージコイルやサドルコイル等の他のコイルを用いることも可能であり、受信コイルの数は３個以上でもよい。また、本発明は、図２に示すようなシリンドリカルタイプのマグネットアセンブリを有するＭＲＩ装置だけでなく、撮影空間の大部分が開放されている、いわゆるオープンタイプのマグネットアセンブリを有するＭＲＩ装置にも適用可能である。
さらに、たとえば、感度領域ＶＳの全領域ではなく、所定の領域までフィッティングにより感度分布を計算するようにしてもよい。
【００７４】
【発明の効果】
以上のように、本発明によれば、パラレルイメージングに用いるコイルの感度分布の精度を向上可能なコイルの感度分布作成方法を提供することができる。
また、本発明によれば、コイルの感度分布の精度向上に伴い、画質を向上させることが可能な磁気共鳴撮影方法および磁気共鳴撮影装置を提供することもできる。
【図面の簡単な説明】
【図１】本発明の一実施の形態に係るＭＲＩ装置の概略構成を示すブロック図である。
【図２】図１に示すＭＲＩ装置に用いるマグネットアセンブリの一例の外観を示す斜視図である。
【図３】パラレルイメージング法の概略について述べるための図である。
【図４】コイルの感度分布について述べるための図である。
【図５】本発明の一実施の形態における感度分布作成、および作成した感度分布を用いた磁気共鳴撮影における処理の流れを示すブロック図である。
【図６】図１に示すＭＲＩ装置における計算機の画像処理部の機能ブロック図である。
【図７】位相のデータに対するフィッティング処理の概念について述べるための図である。
【図８】大きさのデータに対するフィッティング処理の概念について述べるための図である。
【図９】複数のコイルについて最終的に得られる感度分布を表わす図である。
【符号の説明】
１…マグネットアセンブリ
４…被検体
７…計算機
２１…制御部
２３…画像処理部
２４…フーリエ変換部
２５…規格化部
２６…位相差分計算部
２７…フィッティング処理部
２８…合成部
２９…アンフォールド処理部
１００…ＭＲＩ（磁気共鳴撮影）装置
Ｃ１，Ｃ２…受信コイル
【図１】
2004344327.tif 000003
【図２】
2004344327.tif 000004
【図３】
2004344327.tif 000005
【図４】
2004344327.tif 000006
【図５】
2004344327.tif 000007
【図６】
2004344327.tif 000008
【図７】
2004344327.tif 000009
【図８】
2004344327.tif 000010
【図９】
2004344327.tif 000011