磁気共鳴イメージング装置およびプログラム
【課題】リンギングやブラーリングはできるだけ低減することができ、更に、コントラストを調整できるようにする。
【解決手段】エコー時間ＴＥおよびエコー数n_maxなどの値に基づいて、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ＝n_center＝３９を決定する。そして、ｋ空間の原点Ｃに、エコー番号n_center＝３９が割り当てられるように、係数ratio_yzの値を設定し、関数func_n（ky，kz）を決定する。決定した関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てていき、データ収集領域Ｒａｃｑを複数の曲線でn_max個の領域に分割する。
【解決手段】エコー時間ＴＥおよびエコー数n_maxなどの値に基づいて、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ＝n_center＝３９を決定する。そして、ｋ空間の原点Ｃに、エコー番号n_center＝３９が割り当てられるように、係数ratio_yzの値を設定し、関数func_n（ky，kz）を決定する。決定した関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てていき、データ収集領域Ｒａｃｑを複数の曲線でn_max個の領域に分割する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置、およびプログラムに関する。
【背景技術】
【０００２】
３Ｄ（Dimension）撮像法として、３Ｄ ＦＳＥ（Fast Spin
Echo）法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開2005-288195号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
図１５は、３Ｄ ＦＳＥ法で使用されるシーケンスの一例を示す図である。
尚、図１５の説明に当たっては、記号「iter_max」、「iter」、「n_max」、「n」を用いる。各記号の意味は、以下の通りである。
iter_max：データ収集シーケンスＡＣＱｉｔｅｒが繰り返される回数
iter：データ収集シーケンスＡＣＱｉｔｅｒの繰返し番号
n_max：１回のデータ収集シーケンスＡＣＱｉｔｅｒで収集されるエコー数
n：エコー番号
【０００５】
データ収集シーケンスＡＣＱｉｔｅｒ（iter＝１〜iter_max）は、iter_max回繰返し実行される。図１５では、繰返し番号iter＝１、２、ｍ、ｍ＋１、iter_maxのデータ収集シーケンスＡＣＱ１、ＡＣＱ２、ＡＣＱｍ、ＡＣＱｍ＋１、ＡＣＱｉｔｅｒ＿ｍａｘが示されている。例えば、繰返し回数iter_max＝５００の場合、５００回のデータ収集シーケンスＡＣＱ１〜ＡＣＱ５００が実行される。
【０００６】
データ収集シーケンスＡＣＱｉｔｅｒを１回実行することにより、n_max個のエコーＥ１〜Ｅｎ＿ｍａｘが収集される。図１５では、繰返し番号iter＝１のデータ収集シーケンスＡＣＱ１を実行することにより収集されるエコーＥ１〜Ｅｎ＿ｍａｘが示されている。しかし、他の繰返し番号iterのデータ収集シーケンスＡＣＱｉｔｅｒであっても、同様に、エコーＥ１〜Ｅｎ＿ｍａｘが収集される。
【０００７】
図１６〜図１８は、図１５に示す３ＤＦＳＥを実行するときのビュー・オーダリングの一例の説明図である。
【０００８】
図１６は、ｋｙ−ｋｚ面を示す図である。ｋｙ−ｋｚ面は、データの収集が行われるデータ収集領域Ｒａｃｑと、データの収集が行われないデータ非収集領域Ｒｎｏｎとが規定されている。データ収集領域Ｒａｃｑは、ｋｙ−ｋｚ面に内接する円の領域として規定されている。
【０００９】
データ収集領域Ｒａｃｑは、n_max個の領域に分割される。つまり、データ収集領域Ｒａｃｑは、エコー数n_maxと同数の領域に分割される（図１７参照）。
【００１０】
図１７は、データ収集領域Ｒａｃｑがn_max個の領域に分割されたときのｋｙ−ｋｚ面を示す図である。
【００１１】
図１７では、データ収集領域Ｒａｃｑは、ｋｚ軸に平行な直線Ｌによって、面積の等しいn_max個の領域Ｒ１〜Ｒｎ＿ｍａｘに分割されている。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点には、それぞれ、エコー番号n＝１〜n_maxが割り当てられる。例えば、領域Ｒ１のサンプリング点には、エコー番号n＝１が割り当てられており、領域Ｒｎ＿ｍａｘのサンプリング点には、エコー番号n＝n_maxが割り当てられている。尚、図１７では、ｋ空間の原点Ｃを含む領域は、「Ｒｎ＿ｃｅｎｔｅｒ」で表されており、領域Ｒｎ＿ｃｅｎｔｅｒのサンプリング点に割り当てられたエコー番号ｎは、「n_center」で表されている。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、エコー番号n＝１〜n_maxを割り当てた後、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに、サンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。
【００１２】
図１８は、繰返し番号iter＝１〜iter_maxを割り当てる方法の説明図である。
先ず、領域Ｒ１内のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。図１８では、説明の便宜上、繰返し番号iter＝１、２、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられる。以下、同様に、他の領域Ｒ２〜Ｒｎ＿ｍａｘ内のサンプリング点についても、繰返し番号iter＝１〜iter_maxを割り当てる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘには、繰返し番号iter＝１〜iter_maxが割り当てられたサンプリング点が存在する。図１８では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝２が付されたサンプリング点は、ラインＪ２で繋いで示されており、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。
【００１３】
上記のようにして、各領域Ｒ１〜Ｒｎ＿ｍａｘにおいて、サンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。サンプリング点に繰返し番号iterを割り当てた後、データ収集シーケンスＡＣＱ１〜ＡＣＱｉｔｅｒ＿ｍａｘを実行し、収集したエコーをデータ収集領域に配置する（図１９参照）。
【００１４】
図１９は、収集したエコーをデータ収集領域に配置するときの説明図である。
繰返し番号iter＝１のデータ収集シーケンスＡＣＱ１で収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝１のサンプリング点（ラインＪ１のサンプリング点）に配置される。例えば、エコーＥ１は、領域Ｒ１のiter＝１のサンプリング点に配置され、エコーＥ２は、領域Ｒ２のiter＝１のサンプリング点に配置され、エコーＥｎ＿ｍａｘは、領域Ｒｎ＿ｍａｘのiter＝１のサンプリング点に配置される。尚、図１９では、ｋ空間の原点Ｃを含む領域Ｒｎ＿ｃｅｎｔｅｒに配置されるエコーは、記号「Ｅｎ＿ｃｅｎｔｅｒ」で表されている。
【００１５】
以下同様に、繰返し番号iter＝ｉのデータ収集シーケンスＡＣＱｉで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝ｉのサンプリング点に配置される。例えば、繰返し番号iter＝２のデータ収集シーケンスＡＣＱ２で収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝２のサンプリング点（ラインＪ２のサンプリング点）に配置される。また、繰返し番号iter＝ｍのデータ収集シーケンスＡＣＱｍで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝ｍのサンプリング点（ラインＪｍのサンプリング点）に配置される。同様に、繰返し番号iter＝iter_maxのデータ収集シーケンスＡＣＱｉｔｅｒ＿ｍａｘで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝iter_maxのサンプリング点（ラインＪｉｔｅｒ＿ｍａｘ）に配置される。したがって、同じエコー番号ｎのエコーは、同じ領域に配置される。例えば、エコー番号ｎ＝n_centerのエコーＥｎ＿ｃｅｎｔｅｒは、領域Ｒｎ＿ｃｅｎｔｅｒに配置される。
【００１６】
上記のようにして、データ収集領域Ｒａｃｑにデータが配置される。
しかし、上記のビュー・オーダリングでは、データ収集領域Ｒａｃｑをｋｚ方向に延在する直線Ｌで分割しているので（図１７参照）、エコー番号ｎが同一のエコーは、ｋｚ方向に配置される。したがって、ｋｙ方向にリンギング（ringing）やブラーリング（blurring）が発生しやすいという問題がある。
【００１７】
また、上述した３Ｄ ＦＳＥの他に、３Ｄ ＳＳＦＰが知られている。
図２０は、３Ｄ ＳＳＦＰを実行するときのビュー・オーダリングの一例を示す図である。
【００１８】
図２０に示すビュー・オーダリングでは、データ収集領域Ｒａｃｑを、ｋ空間の原点Ｃを中心とした半径の異なる同心円によってｎ_maｘ個の領域Ｒ１〜Ｒｎ＿ｍａｘに分割する。そして、領域Ｒ１内のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。図２０では、説明の便宜上、繰返し番号iter＝１、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられる。以下、同様に、他の領域Ｒ２〜Ｒｎ＿ｍａｘ内のサンプリング点についても、繰返し番号iter＝１〜iter_maxを割り当てられる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘには、繰返し番号iter＝１〜iter_maxが割り当てられたサンプリング点が存在する。図２０では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。このようにして繰返し番号iterが割り当てられ、データが収集される。しかし、図２０に示す方法では、ｋ空間の原点Ｃのエコー番号ｎは、ｎ＝１に限られるので、コントラストの調整を行うことができないという問題がある。
【００１９】
したがって、リンギングやブラーリングはできるだけ低減することができ、更に、コントラストを調整できることが望まれている。
【課題を解決するための手段】
【００２０】
本発明の第１の態様は、ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置であって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割手段、
を有する磁気共鳴イメージング装置である。
【００２１】
本発明の第２の態様は、ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置のプログラムであって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割処理、
を計算機に実行させるためのプログラムである。
【発明の効果】
【００２２】
上記の複数の曲線でデータ収集領域を分割することによって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の一形態の磁気共鳴イメージング装置の概略図である。
【図２】データ収集領域Ｒａｃｑを分割するときに設定される基準点Ｐ（ky0，kz0）を示す図である。
【図３】選び出されたサンプリング点を示す図である。
【図４】エコー番号ｎ＝２が割り当てられるサンプリング点を選び出すときの説明図である。
【図５】エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出すときの説明図である。
【図６】データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割した様子を示す図である。
【図７】領域Ｒ１の拡大図である。
【図８】繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
【図９】各領域Ｒ１〜Ｒｎ＿ｍａｘごとに割り当てられた繰返し番号iter＝１〜iter_maxを示す図である。
【図１０】式（１Ｂ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【図１１】式（１Ｃ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【図１２】本形態のビュー・オーダリングに従って被検体を撮影するときのフローを示す図である。
【図１３】別の基準点を基準にした曲線でデータ収集領域Ｒａｃｑを分割したときの一例を示す図である。
【図１４】円又は楕円の曲線とは別の曲線でデータ収集領域を分割する一例を示す図である。
【図１５】３ＤＦＳＥ法で使用されるシーケンスの一例を示す図である。
【図１６】ｋｙ−ｋｚ面を示す図である。
【図１７】データ収集領域Ｒａｃｑがn_max個の領域に分割されたときのｋｙ−ｋｚ面を示す図である。
【図１８】繰返し番号iter＝１〜iter_maxを割り当てる方法の説明図である。
【図１９】収集したエコーをデータ収集領域に配置するときの説明図である。
【図２０】３Ｄ ＳＳＦＰを実行するときのビュー・オーダリングの一例を示す図である。
【発明を実施するための形態】
【００２４】
以下、発明を実施するための形態を詳細に説明するが、本発明は、以下の形態に限定されることはない。
【００２５】
図１は、本発明の一形態の磁気共鳴イメージング装置の概略図である。
磁気共鳴イメージング装置（以下、ＭＲＩ装置と呼ぶ。ＭＲＩ（Magnetic Resonance Imaging））１００は、磁場発生装置２、テーブル３、受信コイル４などを有している。
【００２６】
磁場発生装置２は、被検体１２が収容されるボア２１、超伝導コイル２２、勾配コイル２３、および送信コイル２４などを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。
【００２７】
テーブル３は、クレードル３１を有している。クレードル３１は、ボア２１に移動できるように構成されている。クレードル３１によって、被検体１２はボア２１に搬送される。
【００２８】
受信コイル４は、被検体１２の腹部から胸部に渡って取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００２９】
ＭＲＩ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、操作部１０、および表示部１１を有している。
【００３０】
シーケンサ５は、中央処理装置９の制御を受けて、被検体１２のスキャンを実行するための情報を送信器６および勾配磁場電源７に送る。
【００３１】
送信器６は、シーケンサ５から送られた情報に基づいて、ＲＦコイル２４を駆動する駆動信号を出力する。
【００３２】
勾配磁場電源７は、シーケンサ５から送られた情報に基づいて、勾配コイル２３を駆動する駆動信号を出力する。
【００３３】
受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、中央処理装置９に伝送する。
【００３４】
中央処理装置９は、シーケンサ５および表示部１１に必要な情報を伝送したり、受信器８から受け取った信号に基づいて画像を再構成するなど、ＭＲＩ装置１００の各種の動作を実現するように、ＭＲＩ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、エコー番号決定手段９１〜分割手段９３などを有している。
【００３５】
エコー番号決定手段９１は、ｋ空間の原点に割り当てられるエコー番号を決定する。
関数決定手段９２は、後述する関数func_n（ky，kz）およびfunc_iter（ky，kz）を決定する。
分割手段９３は、ｋ空間のデータ収集領域を複数の領域に分割する。
【００３６】
中央処理装置９は、エコー番号決定手段９１〜分割手段９３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。中央処理装置９は、課題を解決するための手段に記載された計算機の一例である。
【００３７】
操作部１０は、オペレータ１３の操作に応じて、種々の命令を中央処理装置９に入力する。表示部１１は種々の情報を表示する。
【００３８】
ＭＲＩ装置１００は、上記のように構成されている。
次に、図１５に示すシーケンスを用いてデータを収集するときの本形態におけるｋ空間のビュー・オーダリングの原理について説明する。
【００３９】
尚、以下の説明では、図１６に示すように、ｋｙ−ｋｚ面には、データの収集が行われるデータ収集領域Ｒａｃｑと、データの収集が行われないデータ非収集領域Ｒｎｏｎとが規定されているとする。ただし、データ収集領域Ｒａｃｑは、円以外の形状の領域としてもよい。また、ｋｙ−ｋｚ面の全領域をデータ収集領域Ｒａｃｑとしてもよい。
【００４０】
本形態では、データ収集領域Ｒａｃｑを、複数の領域に分割する。以下に、データ収集領域Ｒａｃｑを複数の領域に分割する方法について、図２〜図６を参照しながら説明する。
【００４１】
図２は、データ収集領域Ｒａｃｑを分割するときに用いられる基準点Ｐ（ky0，kz0）を示す図である。尚、図２では、サンプリング点は黒丸「●」で示されている。
【００４２】
基準点Ｐ（ky0，kz0）は、ｋ空間の原点Ｃとは異なる位置に位置決めされている。本形態では、基準点Ｐ（ky0，kz0）を基準にして、データ収集領域Ｒａｃｑのサンプリング点にエコー番号ｎ＝１〜n_maxを割り当てる。エコー番号ｎの割当ては、以下の式（１）で表される関数func_n（ky，kz）を用いて行われる。
【数１】
式（１）に含まれる係数ratio_yzは、調整可能な値である。ここでは、係数ratio_yz＝１．０であるとする。係数ratio_yz＝１．０のときの関数func_n（ky，kz）は、式（１Ａ）で表される。
【数２】
【００４３】
式（１Ａ）の関数func_n（ky，kz）は、基準点Ｐ（ky0，kz0）とサンプリング点Ｐ（ky，kz）との間の距離を表している。したがって、式（１Ａ）の関数func_n（ky，kz）の値が小さいほど、サンプリング点Ｐ（ky，kz）は基準点Ｐ（ky0，kz0）に近いことを意味し、一方、式（１Ａ）の関数func_n（ky，kz）の値が大きいほど、サンプリング点Ｐ（ky，kz）は基準点Ｐ（ky0，kz0）から離れていることを意味している。図２では、代表して、サンプリング点Ｐαの関数func_n（ky，kz）の値ｒαと、サンプリング点Ｐβの関数func_n（ky，kz）の値ｒβとが示されている。ｒα＜ｒβであるので、サンプリング点Ｐαは、サンプリング点Ｐβよりも、基準点Ｐ（ky0，kz0）に近いことがわかる。
【００４４】
エコー番号ｎ＝１〜n_maxを割り当てる場合、先ず、データ収集領域Ｒａｃｑのサンプリング点の中から、関数func_n（ky，kz）の値の小さい順に、エコー番号ｎ＝１が割り当てられるサンプリング点をiter_max個選び出す。尚、iter_maxとは、図１５に示すデータ収集シーケンスＡＣＱｉｔｅｒが繰り返される回数であるので、例えばiter_max＝５００の場合、エコー番号ｎ＝１が割り当てられるサンプリング点は、５００個選び出されることになる。
【００４５】
図３は、選び出されたサンプリング点を示す図である。
尚、図３では、選び出されたiter_max個のサンプリング点Ｐ１１〜Ｐ１ｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。また、選び出されたサンプリング点のうち、２つのサンプリング点について、符号「Ｐ１１」および符号「Ｐ１ｚ」で示してある。
【００４６】
選び出されたiter_max個のサンプリング点Ｐ１１〜Ｐ１ｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐ１ｚの関数func_n（ky，kz）の値ｒ１ｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒ１ｚの円の曲線ＣＬ１を、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。この曲線ＣＬ１を用いて、サンプリング点Ｐ１１〜Ｐ１ｚの領域Ｒ１を決めることができる。領域Ｒ１に存在するサンプリング点Ｐ１１〜Ｐ１ｚに、エコー番号ｎ＝１が割り当てられる。
【００４７】
サンプリング点Ｐ１１〜Ｐ１ｚにエコー番号ｎ＝１を割り当てたら、次に、エコー番号ｎ＝２が割り当てられるサンプリング点を選び出す（図４参照）。
【００４８】
図４は、エコー番号ｎ＝２が割り当てられるサンプリング点を選び出すときの説明図である。
【００４９】
エコー番号ｎ＝２が割り当てられるサンプリング点を選び出す場合、領域Ｒ１のサンプリング点を除いて、関数func_n（ky，kz）の値の小さい順に、サンプリング点をiter_max個選び出す。ここでは、サンプリング点Ｐ２１〜Ｐ２ｚが、エコー番号ｎ＝２が割り当てられるサンプリング点として選び出されたとする。尚、図４では、選択されたサンプリング点Ｐ２１〜Ｐ２ｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。
【００５０】
選び出されたサンプリング点Ｐ２１〜Ｐ２ｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐ２ｚの関数func_n（ky，kz）の値ｒ２ｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒ２ｚの円の曲線ＣＬ２を、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。サンプリング点Ｐ２１〜Ｐ２ｚの領域Ｒ２は、２本の曲線ＣＬ１およびＣＬ２を用いて決めることができる。領域Ｒ２に存在するサンプリング点Ｐ２１〜Ｐ２ｚに、エコー番号ｎ＝２が割り当てられる。
【００５１】
以下同様に、iter_max個のサンプリング点を選び出しながら、データ収集領域Ｒａｃｑの残りのサンプリング点にもエコー番号ｎを割り当てていく（図５参照）。
【００５２】
図５は、エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出すときの説明図である。
【００５３】
エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出す場合、領域Ｒ１〜Ｒｋ−１のサンプリング点を除いて、関数func_n（ky，kz）の値の小さい順に、サンプリング点をiter_max個選び出す。ここでは、サンプリング点Ｐｋ１〜Ｐｋｚが、エコー番号ｎ＝ｋが割り当てられるサンプリング点として選び出されたとする。尚、図５では、選択されたサンプリング点Ｐｋ１〜Ｐｋｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。
【００５４】
選び出されたサンプリング点Ｐｋ１〜Ｐｋｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐｋｚの関数func_n（ky，kz）の値ｒｋｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒｋｚの円の曲線ＣＬｋを、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。サンプリング点Ｐｋ１〜Ｐｋｚの領域Ｒｋは、２本の曲線ＣＬｋ−１およびＣＬｋを用いて決めることができる。領域Ｒｋに存在するサンプリング点Ｐｋ１〜Ｐｋｚに、エコー番号ｎ＝ｋが割り当てられる。
【００５５】
エコー番号ｎ＝ｋを割り当てた後、以下同様に、関数func_n（ky，kz）の値の小さい順に、iter_max個のサンプリング点を選び出しながら、残りのサンプリング点にもエコー番号ｎを割り当てる。このように、サンプリング点にエコー番号ｎ＝１〜n_maxを割り当てていき、データ収集領域Ｒａｃｑを分割するための曲線を決定することにより、データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割することができる（図６参照）。
【００５６】
図６は、データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割した様子を示す図である。尚、図６では、サンプリング点は、図示省略されている。
【００５７】
図６に示すように、データ収集領域Ｒａｃｑは、曲線ＣＬ１〜ＣＬｗによって複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割される。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点には、エコー番号n＝１〜n_maxが割り当てられる。尚、ｋ空間の原点Ｃを含む領域は、「Ｒｎ＿ｃｅｎｔｅｒ」で表されており、領域Ｒｎ＿ｃｅｎｔｅｒのサンプリング点に割り当てられたエコー番号ｎは、「n_center」で表されている。式（１Ａ）の関数func_n（ky，kz）でデータ収集領域Ｒａｃｑを分割した場合、エコー番号n_centerとエコー数n_maxとの関係は、以下の式で表すことができる。
【数３】
【００５８】
例えば、エコー数n_max＝１００の場合、n_center＝３９となる。したがって、エコー数n_max＝１００の場合、領域Ｒｎ＿ｃｅｎｔｅｒには、エコー番号n_center＝３９が割り当てられる。
【００５９】
図６に示すようにデータ収集領域Ｒａｃｑを分割したら、各領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。先ず、領域Ｒ１のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる（図７および図８参照）。
【００６０】
図７および図８は、領域Ｒ１のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
【００６１】
図７は、領域Ｒ１の拡大図である。領域Ｒ１のサンプリング点は、黒丸「●」で示されている。尚、説明の便宜上、領域Ｒ１には、９個のサンプリング点のみが示されており、９個のサンプリング点のうち、４個のサンプリング点に、符号Ｐ１１、Ｐ１２、Ｐ１３、Ｐ１４が付されている。
【００６２】
本形態では、以下の関数func_iter（ky，kz）を用いて、各サンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。
【数４】
【００６３】
ここでは、式（３）に含まれる係数ratio_yzは、ratio_yz＝１．０であるとする。係数ratio_yz＝１．０のときの関数func_iter（ky，kz）は、式（３Ａ）で表される。
【数５】
【００６４】
以下に、式（３Ａ）の関数func_iter（ky，kz）を用いて繰返し番号iter＝１〜iter_maxを割り当てる方法について説明する（図８参照）。
【００６５】
図８は、繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
関数func_iter（ky，kz）は、基準点Ｐ（ky0，kz0）に対するサンプリング点の角度θを表す関数であり、角度θは、基準線Ｌｒｅｆと線分ＬＳとによって規定される。ここで、基準線Ｌｒｅｆは、基準点Ｐ（ky0，kz0）からｋｚ軸方向に延在する線を表しており、線分ＬＳは、基準点Ｐ（ky0，kz0）と各サンプリング点とを結ぶ線を表している。本形態では、関数func_iter（ky，kz）を用いて、基準線Ｌｒｅｆと線分ＬＳとの成す角度θを求め、角度θの小さい順に、サンプリング点に繰返し番号iterを割り当てる。例えば、図８では、領域Ｒ１に含まれるサンプリング点の中で、サンプリング点Ｐ１２の角度θ＝θ１が最小角となるので、サンプリング点Ｐ１２に繰返し番号iter＝１が割り当てられる。そして、サンプリング点Ｐ１２の次に角度θが小さいのは、サンプリング点Ｐ１３の角度θ＝θ２であるので、サンプリング点Ｐ１３に繰返し番号iter＝２が割り当てられる。尚、図８は、サンプリング点Ｐ１３他に、サンプリング点Ｐ１４も、角度θ＝θ２であるので、角度θが同じになるサンプリング点が複数個存在している。この場合は、基準点Ｐ（ky0，kz0）からの距離が近い順に、小さい繰返し番号iterを割り当てる。したがって、サンプリング点Ｐ１３に繰返し番号iter＝２が割り当てられ、サンプリング点Ｐ１４には返し番号iter＝３が割り当てられる。以下同様に、角度θの小さい順に、繰返し番号iterを割り当てていき、角度θが最大値θmaxとなるサンプリング点Ｐ１１に、繰返し番号iter＝iter_maxを割り当てる。このようにして、領域Ｒ１のサンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。
【００６６】
上記の説明では、領域Ｒ１のサンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる方法について説明しているが、他の領域Ｒ２〜Ｒｎ＿ｍａｘのサンプリング点についても、関数func_iter（ky，kz）を用いて、角度θの小さい順に、繰返し番号iterが割り当てられる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに、繰返し番号iter＝１〜iter_maxが割り当てられる（図９参照）。
【００６７】
図９は、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに割り当てられた繰返し番号iter＝１〜iter_maxを示す図である。
【００６８】
図９では、説明の便宜上、繰返し番号iter＝１、ｋ、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられている。図９では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝ｋが付されたサンプリング点は、ラインＪｋで繋いで示されており、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。
【００６９】
このようにして、各領域Ｒ１〜Ｒｎ＿ｍａｘにおいて、サンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。
【００７０】
尚、上記の説明では、式（１）の係数ratio_yzが、ratio_yz＝１．０の場合について説明されている。しかし、係数ratio_yzは１．０以外の値（つまり、係数ratio_yz≠１．０）であってもよい。例えば、係数ratio_yz＝１．５の場合、関数func_n（ky，kz）は、式（１Ｂ）で表され、係数ratio_yz＝３．０の場合、関数func_n（ky，kz）は、式（１Ｃ）で表される。
【数６】
【００７１】
図１０は、式（１Ｂ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図、図１１は、式（１Ｃ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【００７２】
係数ratio_yz≠１．０の場合、図１０および図１１に示すように、データ収集領域Ｒａｃｑは、基準点Ｐ（ky0，kz0）を基準にした楕円の曲線ＣＬによって、複数の領域に分割される。
【００７３】
図１０の場合、領域Ｒｎ＿ｃｅｎｔｅｒに割り当てられるエコー番号ｎ_centerは、以下の式（４）で表すことができる。また、図１１の場合、領域Ｒｎ＿ｃｅｎｔｅｒに割り当てられるエコー番号ｎ_centerは、以下の式（５）で表すことができる。
ｎ_center＝０．３０＊n_max ・・・（４）
ｎ_center＝０．１６＊n_max ・・・（５）
【００７４】
式（４）では、エコー数n_max＝１００の場合、n_center＝３０となるので、図１０では、ｋ空間の原点Ｃには、エコー番号n＝３０が割り当てられる。一方、式（５）では、エコー数n_max＝１００の場合、n_center＝１６となるので、図１１では、ｋ空間の原点Ｃには、エコー番号n＝１６が割り当てられる。したがって、本形態では、関数func_n（ky，kz）の係数ratio_yzの値によって、ｋ空間の原点Ｃに割り当てられるエコー番号n_centerを調整できることが分かる。関数func_n（ky，kz）の係数ratio_yzの値を設定する方法については、後述する。
【００７５】
データ収集領域Ｒａｃｑを、図１０又は図１１に示すように分割したら、式（３）の関数func_iter（ky，kz）を用いて、各領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。図１０の場合は、係数ratio_yz＝１．５であるので、式（３）のratio_yzに１．５を代入すると、以下の式（３Ｂ）が得られる。一方、図１１の場合は、係数ratio_yz＝３．０であるので、式（３）のratio_yzに３．０を代入すると、以下の式（３Ｃ）が得られる。
【数７】
【００７６】
したがって、図１０では、式（３Ｂ）で表される関数func_iter（ky，kz）を用いて、繰返し番号iter＝１〜iter_maxが割り当てられる。一方、図１１では、式（３Ｃ）で表される関数func_iter（ky，kz）を用いて、繰返し番号iter＝１〜iter_maxが割り当てられる。
以上のようにして、ビュー・オーダリングが行われる。
【００７７】
本形態では、データ収集領域Ｒａｃｑを分割するための曲線は、ｋｚ方向だけでなく、ｋｙ方向にも変化する。したがって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【００７８】
次に、本形態のビュー・オーダリングに従って被検体を撮影するときのフローについて説明する。尚、以下の説明では、図１５に示す３Ｄ
ＦＳＥのシーケンスを用いてデータを収集する場合のフローについて説明するが、本発明は、３Ｄ ＦＳＥのシーケンスの他に、３Ｄ ＳＳＦＰのシーケンス等の他の３Ｄシーケンスを用いてデータを収集する場合にも適用することができる。
【００７９】
図１２は、本形態のビュー・オーダリングに従って被検体を撮影するときのフローを示す図である。
【００８０】
ステップＳＴ１では、オペレータ１３が、図１５に示すシーケンスのエコー時間ＴＥおよびエコー数ｎ_maxを設定する。ここでは、オペレータ１３は、エコー数n_max＝１００に設定したとする。
【００８１】
ステップＳＴ２では、エコー番号決定手段９１（図１参照）が、エコー時間ＴＥおよびエコー数n_maxなどの値に基づいて、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ＝n_centerを決定する。ここでは、ｎ_center＝３９と決定されたとする。エコー番号ｎ_centerを決定したら、ステップＳＴ３に進む。
【００８２】
ステップＳＴ３では、関数決定手段９２（図１参照）が、先ず、ｋ空間の原点Ｃに、ステップＳＴ２で決定されたエコー番号n_center＝３９が割り当てられるように、式（１）の係数ratio_yzの値を設定する。係数ratio_yzの値を設定する方法の一例としては、エコー数n_maxの値と、エコー番号n_centerの値と、係数ratio_yzの値との対応関係を事前に作成しておく方法がある。本形態では、ステップＳＴ１において、オペレータ１３によってエコー数n_max＝１００が設定されており、ステップＳＴ２において、エコー番号n_center＝３９が決定されているので、上記の対応関係を事前に作成しておくことにより、係数ratio_yzの値を設定することができる。ここでは、係数ratio_yz＝１．０に設定されたとする。関数決定手段９２は、設定した係数ratio_yz＝１．０を式（１）に代入する。係数ratio_yz＝１．０を式（１）に代入することにより、データ収集領域を分割するために用いられる関数func_n（ky，kz）が決定される。係数ratio_yz＝１．０の場合、決定された関数func_n（ky，kz）は、式（１Ａ）で表される。関数func_n（ky，kz）を決定したら、ステップＳＴ４に進む。
【００８３】
ステップＳＴ４では、分割手段９３（図１参照）が、ステップＳＴ３で決定された関数func_n（ky，kz）を用いてデータ収集領域Ｒａｃｑのサンプリング点にエコー番号を割り当てる。そして、割り当てられたエコー番号に基づいて、図６に示すように、データ収集領域Ｒａｃｑに曲線ＣＬ１〜ＣＬｚを規定し、複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割する。分割したら、ステップＳＴ５に進む。
【００８４】
ステップＳＴ５では、関数決定手段９２が、ステップＳＴ３で設定した係数ratio_yz＝１．０を式（３）に代入する。係数ratio_yz＝１．０を式（３）に代入することにより、領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に繰返し番号を割り当てるために用いられる関数func_iter（ky，kz）が決定される。係数ratio_yz＝１．０の場合、決定された関数func_iter（ky，kz）は、式（３Ａ）で表される。関数func_iter（ky，kz）を決定したら、ステップＳＴ６に進む。
【００８５】
ステップＳＴ６では、ステップＳＴ５で決定された関数fucn_iterを用いて、領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、番号iterが割り当てられる（図９参照）。番号iterを割り当てたら、ステップＳＴ７に進み、データを収集するためのスキャンを実行し、フローを終了する。
【００８６】
本形態では、データ収集領域Ｒａｃｑを分割するための曲線は、ｋｚ方向だけでなく、ｋｙ方向にも変化する。したがって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【００８７】
また、関数func_n（ky，kz）の基準点Ｐ（ky0，kz0）は、ｋ空間の原点Ｃからずれている。したがって、関数func_n（ky，kz）に含まれる係数ratio_yzの値を変えるだけで、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ_centerを調整することができるので、コントラストを容易に調整することができる。
【００８８】
本形態では、データ収集領域Ｒａｃｑを分割する曲線は、基準点Ｐ（ky0，kz0）を基準にして規定されている（例えば、図６参照）。しかし、別の基準点を基準にして曲線を規定してもよい（図１３参照）。
【００８９】
図１３は、別の基準点を基準にした曲線でデータ収集領域Ｒａｃｑを分割したときの一例を示す図である。
【００９０】
図１３では、基準点Ｐ′を基準にして、データ収集領域Ｒａｃｑを分割する曲線ＣＬ１〜ＣＬｗを規定している。基準点Ｐ′は、図６に示す基準点Ｐ（ky0，kz0）よりも、ｋ空間の原点Ｃに近い位置に設定されている。このように、データ収集領域Ｒａｃｑを分割する曲線の基準点は、ｋ空間の原点Ｃとは異なる位置であれば、様々な位置に設定することが可能である。
【００９１】
また、本形態では、データ収集領域Ｒａｃｑを円又は楕円の曲線で分割しているが、別の曲線で分割してもよい（図１４参照）。
【００９２】
図１４は、円又は楕円の曲線とは別の曲線でデータ収集領域を分割する一例を示す図である。
【００９３】
図１４では、基準点Ｐ（ky0，kz0）を基準にした菱形の曲線でデータ収集領域Ｒａｃｑを分割した場合が示されている。このように、データ収集領域Ｒａｃｑを分割する曲線は、円、楕円、双曲線、放物線などのような曲線であってもよいし、菱形のように複数の直線の組合せで構成される曲線であってもよい。
【００９４】
また、本形態では、関数func_iter（ky，kz）は、サンプリング点の位置を角度で表しているが、角度以外の別のパラメータを用いてサンプリング点の位置を表してもよい。
【符号の説明】
【００９５】
２ 磁場発生装置
３ テーブル
４ 受信コイル
５ シーケンサ
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 中央処理装置
１０ 操作部
１１ 表示部
１２ 被検体
１３ オペレータ
２１ ボア
２２ 超伝導コイル
２３ 勾配コイル
２４ 送信コイル
３１ クレードル
９１ エコー番号決定手段
９２ 関数決定手段
９３ 分割手段
９４ 番号割当て手段
１００ ＭＲＩ装置
【技術分野】
【０００１】
本発明は、エコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置、およびプログラムに関する。
【背景技術】
【０００２】
３Ｄ（Dimension）撮像法として、３Ｄ ＦＳＥ（Fast Spin
Echo）法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開2005-288195号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
図１５は、３Ｄ ＦＳＥ法で使用されるシーケンスの一例を示す図である。
尚、図１５の説明に当たっては、記号「iter_max」、「iter」、「n_max」、「n」を用いる。各記号の意味は、以下の通りである。
iter_max：データ収集シーケンスＡＣＱｉｔｅｒが繰り返される回数
iter：データ収集シーケンスＡＣＱｉｔｅｒの繰返し番号
n_max：１回のデータ収集シーケンスＡＣＱｉｔｅｒで収集されるエコー数
n：エコー番号
【０００５】
データ収集シーケンスＡＣＱｉｔｅｒ（iter＝１〜iter_max）は、iter_max回繰返し実行される。図１５では、繰返し番号iter＝１、２、ｍ、ｍ＋１、iter_maxのデータ収集シーケンスＡＣＱ１、ＡＣＱ２、ＡＣＱｍ、ＡＣＱｍ＋１、ＡＣＱｉｔｅｒ＿ｍａｘが示されている。例えば、繰返し回数iter_max＝５００の場合、５００回のデータ収集シーケンスＡＣＱ１〜ＡＣＱ５００が実行される。
【０００６】
データ収集シーケンスＡＣＱｉｔｅｒを１回実行することにより、n_max個のエコーＥ１〜Ｅｎ＿ｍａｘが収集される。図１５では、繰返し番号iter＝１のデータ収集シーケンスＡＣＱ１を実行することにより収集されるエコーＥ１〜Ｅｎ＿ｍａｘが示されている。しかし、他の繰返し番号iterのデータ収集シーケンスＡＣＱｉｔｅｒであっても、同様に、エコーＥ１〜Ｅｎ＿ｍａｘが収集される。
【０００７】
図１６〜図１８は、図１５に示す３ＤＦＳＥを実行するときのビュー・オーダリングの一例の説明図である。
【０００８】
図１６は、ｋｙ−ｋｚ面を示す図である。ｋｙ−ｋｚ面は、データの収集が行われるデータ収集領域Ｒａｃｑと、データの収集が行われないデータ非収集領域Ｒｎｏｎとが規定されている。データ収集領域Ｒａｃｑは、ｋｙ−ｋｚ面に内接する円の領域として規定されている。
【０００９】
データ収集領域Ｒａｃｑは、n_max個の領域に分割される。つまり、データ収集領域Ｒａｃｑは、エコー数n_maxと同数の領域に分割される（図１７参照）。
【００１０】
図１７は、データ収集領域Ｒａｃｑがn_max個の領域に分割されたときのｋｙ−ｋｚ面を示す図である。
【００１１】
図１７では、データ収集領域Ｒａｃｑは、ｋｚ軸に平行な直線Ｌによって、面積の等しいn_max個の領域Ｒ１〜Ｒｎ＿ｍａｘに分割されている。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点には、それぞれ、エコー番号n＝１〜n_maxが割り当てられる。例えば、領域Ｒ１のサンプリング点には、エコー番号n＝１が割り当てられており、領域Ｒｎ＿ｍａｘのサンプリング点には、エコー番号n＝n_maxが割り当てられている。尚、図１７では、ｋ空間の原点Ｃを含む領域は、「Ｒｎ＿ｃｅｎｔｅｒ」で表されており、領域Ｒｎ＿ｃｅｎｔｅｒのサンプリング点に割り当てられたエコー番号ｎは、「n_center」で表されている。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、エコー番号n＝１〜n_maxを割り当てた後、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに、サンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。
【００１２】
図１８は、繰返し番号iter＝１〜iter_maxを割り当てる方法の説明図である。
先ず、領域Ｒ１内のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。図１８では、説明の便宜上、繰返し番号iter＝１、２、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられる。以下、同様に、他の領域Ｒ２〜Ｒｎ＿ｍａｘ内のサンプリング点についても、繰返し番号iter＝１〜iter_maxを割り当てる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘには、繰返し番号iter＝１〜iter_maxが割り当てられたサンプリング点が存在する。図１８では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝２が付されたサンプリング点は、ラインＪ２で繋いで示されており、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。
【００１３】
上記のようにして、各領域Ｒ１〜Ｒｎ＿ｍａｘにおいて、サンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。サンプリング点に繰返し番号iterを割り当てた後、データ収集シーケンスＡＣＱ１〜ＡＣＱｉｔｅｒ＿ｍａｘを実行し、収集したエコーをデータ収集領域に配置する（図１９参照）。
【００１４】
図１９は、収集したエコーをデータ収集領域に配置するときの説明図である。
繰返し番号iter＝１のデータ収集シーケンスＡＣＱ１で収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝１のサンプリング点（ラインＪ１のサンプリング点）に配置される。例えば、エコーＥ１は、領域Ｒ１のiter＝１のサンプリング点に配置され、エコーＥ２は、領域Ｒ２のiter＝１のサンプリング点に配置され、エコーＥｎ＿ｍａｘは、領域Ｒｎ＿ｍａｘのiter＝１のサンプリング点に配置される。尚、図１９では、ｋ空間の原点Ｃを含む領域Ｒｎ＿ｃｅｎｔｅｒに配置されるエコーは、記号「Ｅｎ＿ｃｅｎｔｅｒ」で表されている。
【００１５】
以下同様に、繰返し番号iter＝ｉのデータ収集シーケンスＡＣＱｉで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝ｉのサンプリング点に配置される。例えば、繰返し番号iter＝２のデータ収集シーケンスＡＣＱ２で収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝２のサンプリング点（ラインＪ２のサンプリング点）に配置される。また、繰返し番号iter＝ｍのデータ収集シーケンスＡＣＱｍで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝ｍのサンプリング点（ラインＪｍのサンプリング点）に配置される。同様に、繰返し番号iter＝iter_maxのデータ収集シーケンスＡＣＱｉｔｅｒ＿ｍａｘで収集されたエコーＥ１〜Ｅｎ＿ｍａｘは、それぞれ、領域Ｒ１〜Ｒｎ＿ｍａｘのiter＝iter_maxのサンプリング点（ラインＪｉｔｅｒ＿ｍａｘ）に配置される。したがって、同じエコー番号ｎのエコーは、同じ領域に配置される。例えば、エコー番号ｎ＝n_centerのエコーＥｎ＿ｃｅｎｔｅｒは、領域Ｒｎ＿ｃｅｎｔｅｒに配置される。
【００１６】
上記のようにして、データ収集領域Ｒａｃｑにデータが配置される。
しかし、上記のビュー・オーダリングでは、データ収集領域Ｒａｃｑをｋｚ方向に延在する直線Ｌで分割しているので（図１７参照）、エコー番号ｎが同一のエコーは、ｋｚ方向に配置される。したがって、ｋｙ方向にリンギング（ringing）やブラーリング（blurring）が発生しやすいという問題がある。
【００１７】
また、上述した３Ｄ ＦＳＥの他に、３Ｄ ＳＳＦＰが知られている。
図２０は、３Ｄ ＳＳＦＰを実行するときのビュー・オーダリングの一例を示す図である。
【００１８】
図２０に示すビュー・オーダリングでは、データ収集領域Ｒａｃｑを、ｋ空間の原点Ｃを中心とした半径の異なる同心円によってｎ_maｘ個の領域Ｒ１〜Ｒｎ＿ｍａｘに分割する。そして、領域Ｒ１内のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。図２０では、説明の便宜上、繰返し番号iter＝１、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられる。以下、同様に、他の領域Ｒ２〜Ｒｎ＿ｍａｘ内のサンプリング点についても、繰返し番号iter＝１〜iter_maxを割り当てられる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘには、繰返し番号iter＝１〜iter_maxが割り当てられたサンプリング点が存在する。図２０では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。このようにして繰返し番号iterが割り当てられ、データが収集される。しかし、図２０に示す方法では、ｋ空間の原点Ｃのエコー番号ｎは、ｎ＝１に限られるので、コントラストの調整を行うことができないという問題がある。
【００１９】
したがって、リンギングやブラーリングはできるだけ低減することができ、更に、コントラストを調整できることが望まれている。
【課題を解決するための手段】
【００２０】
本発明の第１の態様は、ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置であって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割手段、
を有する磁気共鳴イメージング装置である。
【００２１】
本発明の第２の態様は、ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置のプログラムであって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割処理、
を計算機に実行させるためのプログラムである。
【発明の効果】
【００２２】
上記の複数の曲線でデータ収集領域を分割することによって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の一形態の磁気共鳴イメージング装置の概略図である。
【図２】データ収集領域Ｒａｃｑを分割するときに設定される基準点Ｐ（ky0，kz0）を示す図である。
【図３】選び出されたサンプリング点を示す図である。
【図４】エコー番号ｎ＝２が割り当てられるサンプリング点を選び出すときの説明図である。
【図５】エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出すときの説明図である。
【図６】データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割した様子を示す図である。
【図７】領域Ｒ１の拡大図である。
【図８】繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
【図９】各領域Ｒ１〜Ｒｎ＿ｍａｘごとに割り当てられた繰返し番号iter＝１〜iter_maxを示す図である。
【図１０】式（１Ｂ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【図１１】式（１Ｃ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【図１２】本形態のビュー・オーダリングに従って被検体を撮影するときのフローを示す図である。
【図１３】別の基準点を基準にした曲線でデータ収集領域Ｒａｃｑを分割したときの一例を示す図である。
【図１４】円又は楕円の曲線とは別の曲線でデータ収集領域を分割する一例を示す図である。
【図１５】３ＤＦＳＥ法で使用されるシーケンスの一例を示す図である。
【図１６】ｋｙ−ｋｚ面を示す図である。
【図１７】データ収集領域Ｒａｃｑがn_max個の領域に分割されたときのｋｙ−ｋｚ面を示す図である。
【図１８】繰返し番号iter＝１〜iter_maxを割り当てる方法の説明図である。
【図１９】収集したエコーをデータ収集領域に配置するときの説明図である。
【図２０】３Ｄ ＳＳＦＰを実行するときのビュー・オーダリングの一例を示す図である。
【発明を実施するための形態】
【００２４】
以下、発明を実施するための形態を詳細に説明するが、本発明は、以下の形態に限定されることはない。
【００２５】
図１は、本発明の一形態の磁気共鳴イメージング装置の概略図である。
磁気共鳴イメージング装置（以下、ＭＲＩ装置と呼ぶ。ＭＲＩ（Magnetic Resonance Imaging））１００は、磁場発生装置２、テーブル３、受信コイル４などを有している。
【００２６】
磁場発生装置２は、被検体１２が収容されるボア２１、超伝導コイル２２、勾配コイル２３、および送信コイル２４などを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。
【００２７】
テーブル３は、クレードル３１を有している。クレードル３１は、ボア２１に移動できるように構成されている。クレードル３１によって、被検体１２はボア２１に搬送される。
【００２８】
受信コイル４は、被検体１２の腹部から胸部に渡って取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００２９】
ＭＲＩ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、操作部１０、および表示部１１を有している。
【００３０】
シーケンサ５は、中央処理装置９の制御を受けて、被検体１２のスキャンを実行するための情報を送信器６および勾配磁場電源７に送る。
【００３１】
送信器６は、シーケンサ５から送られた情報に基づいて、ＲＦコイル２４を駆動する駆動信号を出力する。
【００３２】
勾配磁場電源７は、シーケンサ５から送られた情報に基づいて、勾配コイル２３を駆動する駆動信号を出力する。
【００３３】
受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、中央処理装置９に伝送する。
【００３４】
中央処理装置９は、シーケンサ５および表示部１１に必要な情報を伝送したり、受信器８から受け取った信号に基づいて画像を再構成するなど、ＭＲＩ装置１００の各種の動作を実現するように、ＭＲＩ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、エコー番号決定手段９１〜分割手段９３などを有している。
【００３５】
エコー番号決定手段９１は、ｋ空間の原点に割り当てられるエコー番号を決定する。
関数決定手段９２は、後述する関数func_n（ky，kz）およびfunc_iter（ky，kz）を決定する。
分割手段９３は、ｋ空間のデータ収集領域を複数の領域に分割する。
【００３６】
中央処理装置９は、エコー番号決定手段９１〜分割手段９３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。中央処理装置９は、課題を解決するための手段に記載された計算機の一例である。
【００３７】
操作部１０は、オペレータ１３の操作に応じて、種々の命令を中央処理装置９に入力する。表示部１１は種々の情報を表示する。
【００３８】
ＭＲＩ装置１００は、上記のように構成されている。
次に、図１５に示すシーケンスを用いてデータを収集するときの本形態におけるｋ空間のビュー・オーダリングの原理について説明する。
【００３９】
尚、以下の説明では、図１６に示すように、ｋｙ−ｋｚ面には、データの収集が行われるデータ収集領域Ｒａｃｑと、データの収集が行われないデータ非収集領域Ｒｎｏｎとが規定されているとする。ただし、データ収集領域Ｒａｃｑは、円以外の形状の領域としてもよい。また、ｋｙ−ｋｚ面の全領域をデータ収集領域Ｒａｃｑとしてもよい。
【００４０】
本形態では、データ収集領域Ｒａｃｑを、複数の領域に分割する。以下に、データ収集領域Ｒａｃｑを複数の領域に分割する方法について、図２〜図６を参照しながら説明する。
【００４１】
図２は、データ収集領域Ｒａｃｑを分割するときに用いられる基準点Ｐ（ky0，kz0）を示す図である。尚、図２では、サンプリング点は黒丸「●」で示されている。
【００４２】
基準点Ｐ（ky0，kz0）は、ｋ空間の原点Ｃとは異なる位置に位置決めされている。本形態では、基準点Ｐ（ky0，kz0）を基準にして、データ収集領域Ｒａｃｑのサンプリング点にエコー番号ｎ＝１〜n_maxを割り当てる。エコー番号ｎの割当ては、以下の式（１）で表される関数func_n（ky，kz）を用いて行われる。
【数１】
式（１）に含まれる係数ratio_yzは、調整可能な値である。ここでは、係数ratio_yz＝１．０であるとする。係数ratio_yz＝１．０のときの関数func_n（ky，kz）は、式（１Ａ）で表される。
【数２】
【００４３】
式（１Ａ）の関数func_n（ky，kz）は、基準点Ｐ（ky0，kz0）とサンプリング点Ｐ（ky，kz）との間の距離を表している。したがって、式（１Ａ）の関数func_n（ky，kz）の値が小さいほど、サンプリング点Ｐ（ky，kz）は基準点Ｐ（ky0，kz0）に近いことを意味し、一方、式（１Ａ）の関数func_n（ky，kz）の値が大きいほど、サンプリング点Ｐ（ky，kz）は基準点Ｐ（ky0，kz0）から離れていることを意味している。図２では、代表して、サンプリング点Ｐαの関数func_n（ky，kz）の値ｒαと、サンプリング点Ｐβの関数func_n（ky，kz）の値ｒβとが示されている。ｒα＜ｒβであるので、サンプリング点Ｐαは、サンプリング点Ｐβよりも、基準点Ｐ（ky0，kz0）に近いことがわかる。
【００４４】
エコー番号ｎ＝１〜n_maxを割り当てる場合、先ず、データ収集領域Ｒａｃｑのサンプリング点の中から、関数func_n（ky，kz）の値の小さい順に、エコー番号ｎ＝１が割り当てられるサンプリング点をiter_max個選び出す。尚、iter_maxとは、図１５に示すデータ収集シーケンスＡＣＱｉｔｅｒが繰り返される回数であるので、例えばiter_max＝５００の場合、エコー番号ｎ＝１が割り当てられるサンプリング点は、５００個選び出されることになる。
【００４５】
図３は、選び出されたサンプリング点を示す図である。
尚、図３では、選び出されたiter_max個のサンプリング点Ｐ１１〜Ｐ１ｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。また、選び出されたサンプリング点のうち、２つのサンプリング点について、符号「Ｐ１１」および符号「Ｐ１ｚ」で示してある。
【００４６】
選び出されたiter_max個のサンプリング点Ｐ１１〜Ｐ１ｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐ１ｚの関数func_n（ky，kz）の値ｒ１ｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒ１ｚの円の曲線ＣＬ１を、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。この曲線ＣＬ１を用いて、サンプリング点Ｐ１１〜Ｐ１ｚの領域Ｒ１を決めることができる。領域Ｒ１に存在するサンプリング点Ｐ１１〜Ｐ１ｚに、エコー番号ｎ＝１が割り当てられる。
【００４７】
サンプリング点Ｐ１１〜Ｐ１ｚにエコー番号ｎ＝１を割り当てたら、次に、エコー番号ｎ＝２が割り当てられるサンプリング点を選び出す（図４参照）。
【００４８】
図４は、エコー番号ｎ＝２が割り当てられるサンプリング点を選び出すときの説明図である。
【００４９】
エコー番号ｎ＝２が割り当てられるサンプリング点を選び出す場合、領域Ｒ１のサンプリング点を除いて、関数func_n（ky，kz）の値の小さい順に、サンプリング点をiter_max個選び出す。ここでは、サンプリング点Ｐ２１〜Ｐ２ｚが、エコー番号ｎ＝２が割り当てられるサンプリング点として選び出されたとする。尚、図４では、選択されたサンプリング点Ｐ２１〜Ｐ２ｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。
【００５０】
選び出されたサンプリング点Ｐ２１〜Ｐ２ｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐ２ｚの関数func_n（ky，kz）の値ｒ２ｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒ２ｚの円の曲線ＣＬ２を、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。サンプリング点Ｐ２１〜Ｐ２ｚの領域Ｒ２は、２本の曲線ＣＬ１およびＣＬ２を用いて決めることができる。領域Ｒ２に存在するサンプリング点Ｐ２１〜Ｐ２ｚに、エコー番号ｎ＝２が割り当てられる。
【００５１】
以下同様に、iter_max個のサンプリング点を選び出しながら、データ収集領域Ｒａｃｑの残りのサンプリング点にもエコー番号ｎを割り当てていく（図５参照）。
【００５２】
図５は、エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出すときの説明図である。
【００５３】
エコー番号ｎ＝ｋが割り当てられるサンプリング点を選び出す場合、領域Ｒ１〜Ｒｋ−１のサンプリング点を除いて、関数func_n（ky，kz）の値の小さい順に、サンプリング点をiter_max個選び出す。ここでは、サンプリング点Ｐｋ１〜Ｐｋｚが、エコー番号ｎ＝ｋが割り当てられるサンプリング点として選び出されたとする。尚、図５では、選択されたサンプリング点Ｐｋ１〜Ｐｋｚのみを黒丸で示しており、他のサンプリング点は図示省略されている。
【００５４】
選び出されたサンプリング点Ｐｋ１〜Ｐｋｚの関数func_n（ky，kz）の値で、最も大きい値は、サンプリング点Ｐｋｚの関数func_n（ky，kz）の値ｒｋｚである。したがって、基準点Ｐ（ky0，kz0）を中心とした半径ｒｋｚの円の曲線ＣＬｋを、データ収集領域Ｒａｃｑを分割するための曲線として用いることができる。サンプリング点Ｐｋ１〜Ｐｋｚの領域Ｒｋは、２本の曲線ＣＬｋ−１およびＣＬｋを用いて決めることができる。領域Ｒｋに存在するサンプリング点Ｐｋ１〜Ｐｋｚに、エコー番号ｎ＝ｋが割り当てられる。
【００５５】
エコー番号ｎ＝ｋを割り当てた後、以下同様に、関数func_n（ky，kz）の値の小さい順に、iter_max個のサンプリング点を選び出しながら、残りのサンプリング点にもエコー番号ｎを割り当てる。このように、サンプリング点にエコー番号ｎ＝１〜n_maxを割り当てていき、データ収集領域Ｒａｃｑを分割するための曲線を決定することにより、データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割することができる（図６参照）。
【００５６】
図６は、データ収集領域Ｒａｃｑを複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割した様子を示す図である。尚、図６では、サンプリング点は、図示省略されている。
【００５７】
図６に示すように、データ収集領域Ｒａｃｑは、曲線ＣＬ１〜ＣＬｗによって複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割される。領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点には、エコー番号n＝１〜n_maxが割り当てられる。尚、ｋ空間の原点Ｃを含む領域は、「Ｒｎ＿ｃｅｎｔｅｒ」で表されており、領域Ｒｎ＿ｃｅｎｔｅｒのサンプリング点に割り当てられたエコー番号ｎは、「n_center」で表されている。式（１Ａ）の関数func_n（ky，kz）でデータ収集領域Ｒａｃｑを分割した場合、エコー番号n_centerとエコー数n_maxとの関係は、以下の式で表すことができる。
【数３】
【００５８】
例えば、エコー数n_max＝１００の場合、n_center＝３９となる。したがって、エコー数n_max＝１００の場合、領域Ｒｎ＿ｃｅｎｔｅｒには、エコー番号n_center＝３９が割り当てられる。
【００５９】
図６に示すようにデータ収集領域Ｒａｃｑを分割したら、各領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる。先ず、領域Ｒ１のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てる（図７および図８参照）。
【００６０】
図７および図８は、領域Ｒ１のサンプリング点に、繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
【００６１】
図７は、領域Ｒ１の拡大図である。領域Ｒ１のサンプリング点は、黒丸「●」で示されている。尚、説明の便宜上、領域Ｒ１には、９個のサンプリング点のみが示されており、９個のサンプリング点のうち、４個のサンプリング点に、符号Ｐ１１、Ｐ１２、Ｐ１３、Ｐ１４が付されている。
【００６２】
本形態では、以下の関数func_iter（ky，kz）を用いて、各サンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。
【数４】
【００６３】
ここでは、式（３）に含まれる係数ratio_yzは、ratio_yz＝１．０であるとする。係数ratio_yz＝１．０のときの関数func_iter（ky，kz）は、式（３Ａ）で表される。
【数５】
【００６４】
以下に、式（３Ａ）の関数func_iter（ky，kz）を用いて繰返し番号iter＝１〜iter_maxを割り当てる方法について説明する（図８参照）。
【００６５】
図８は、繰返し番号iter＝１〜iter_maxを割り当てるときの説明図である。
関数func_iter（ky，kz）は、基準点Ｐ（ky0，kz0）に対するサンプリング点の角度θを表す関数であり、角度θは、基準線Ｌｒｅｆと線分ＬＳとによって規定される。ここで、基準線Ｌｒｅｆは、基準点Ｐ（ky0，kz0）からｋｚ軸方向に延在する線を表しており、線分ＬＳは、基準点Ｐ（ky0，kz0）と各サンプリング点とを結ぶ線を表している。本形態では、関数func_iter（ky，kz）を用いて、基準線Ｌｒｅｆと線分ＬＳとの成す角度θを求め、角度θの小さい順に、サンプリング点に繰返し番号iterを割り当てる。例えば、図８では、領域Ｒ１に含まれるサンプリング点の中で、サンプリング点Ｐ１２の角度θ＝θ１が最小角となるので、サンプリング点Ｐ１２に繰返し番号iter＝１が割り当てられる。そして、サンプリング点Ｐ１２の次に角度θが小さいのは、サンプリング点Ｐ１３の角度θ＝θ２であるので、サンプリング点Ｐ１３に繰返し番号iter＝２が割り当てられる。尚、図８は、サンプリング点Ｐ１３他に、サンプリング点Ｐ１４も、角度θ＝θ２であるので、角度θが同じになるサンプリング点が複数個存在している。この場合は、基準点Ｐ（ky0，kz0）からの距離が近い順に、小さい繰返し番号iterを割り当てる。したがって、サンプリング点Ｐ１３に繰返し番号iter＝２が割り当てられ、サンプリング点Ｐ１４には返し番号iter＝３が割り当てられる。以下同様に、角度θの小さい順に、繰返し番号iterを割り当てていき、角度θが最大値θmaxとなるサンプリング点Ｐ１１に、繰返し番号iter＝iter_maxを割り当てる。このようにして、領域Ｒ１のサンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。
【００６６】
上記の説明では、領域Ｒ１のサンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる方法について説明しているが、他の領域Ｒ２〜Ｒｎ＿ｍａｘのサンプリング点についても、関数func_iter（ky，kz）を用いて、角度θの小さい順に、繰返し番号iterが割り当てられる。したがって、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに、繰返し番号iter＝１〜iter_maxが割り当てられる（図９参照）。
【００６７】
図９は、各領域Ｒ１〜Ｒｎ＿ｍａｘごとに割り当てられた繰返し番号iter＝１〜iter_maxを示す図である。
【００６８】
図９では、説明の便宜上、繰返し番号iter＝１、ｋ、ｍ、およびiter_maxが割り当てられたサンプリング点のみを黒丸「●」で示してあるが、他のサンプリング点も、繰返し番号iterが割り当てられている。図９では、同じ繰返し番号iterが割り当てられたサンプリング点は、１本のライン上に示されている。例えば、繰返し番号iter＝１が付されたサンプリング点は、ラインＪ１で繋いで示されている。同様に、繰返し番号iter＝ｋが付されたサンプリング点は、ラインＪｋで繋いで示されており、繰返し番号iter＝ｍが付されたサンプリング点は、ラインＪｍで繋いで示されており、繰返し番号iter＝iter_maxが付されたサンプリング点は、ラインＪｉｔｅｒ＿ｍａｘで繋いで示されている。
【００６９】
このようにして、各領域Ｒ１〜Ｒｎ＿ｍａｘにおいて、サンプリング点に繰返し番号iter＝１〜iter_maxが割り当てられる。
【００７０】
尚、上記の説明では、式（１）の係数ratio_yzが、ratio_yz＝１．０の場合について説明されている。しかし、係数ratio_yzは１．０以外の値（つまり、係数ratio_yz≠１．０）であってもよい。例えば、係数ratio_yz＝１．５の場合、関数func_n（ky，kz）は、式（１Ｂ）で表され、係数ratio_yz＝３．０の場合、関数func_n（ky，kz）は、式（１Ｃ）で表される。
【数６】
【００７１】
図１０は、式（１Ｂ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図、図１１は、式（１Ｃ）の関数func_n（ky，kz）を用いてサンプリング点にエコー番号ｎを割り当てたときのデータ収集領域Ｒａｃｑの分割の様子を示す図である。
【００７２】
係数ratio_yz≠１．０の場合、図１０および図１１に示すように、データ収集領域Ｒａｃｑは、基準点Ｐ（ky0，kz0）を基準にした楕円の曲線ＣＬによって、複数の領域に分割される。
【００７３】
図１０の場合、領域Ｒｎ＿ｃｅｎｔｅｒに割り当てられるエコー番号ｎ_centerは、以下の式（４）で表すことができる。また、図１１の場合、領域Ｒｎ＿ｃｅｎｔｅｒに割り当てられるエコー番号ｎ_centerは、以下の式（５）で表すことができる。
ｎ_center＝０．３０＊n_max ・・・（４）
ｎ_center＝０．１６＊n_max ・・・（５）
【００７４】
式（４）では、エコー数n_max＝１００の場合、n_center＝３０となるので、図１０では、ｋ空間の原点Ｃには、エコー番号n＝３０が割り当てられる。一方、式（５）では、エコー数n_max＝１００の場合、n_center＝１６となるので、図１１では、ｋ空間の原点Ｃには、エコー番号n＝１６が割り当てられる。したがって、本形態では、関数func_n（ky，kz）の係数ratio_yzの値によって、ｋ空間の原点Ｃに割り当てられるエコー番号n_centerを調整できることが分かる。関数func_n（ky，kz）の係数ratio_yzの値を設定する方法については、後述する。
【００７５】
データ収集領域Ｒａｃｑを、図１０又は図１１に示すように分割したら、式（３）の関数func_iter（ky，kz）を用いて、各領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に繰返し番号iter＝１〜iter_maxを割り当てる。図１０の場合は、係数ratio_yz＝１．５であるので、式（３）のratio_yzに１．５を代入すると、以下の式（３Ｂ）が得られる。一方、図１１の場合は、係数ratio_yz＝３．０であるので、式（３）のratio_yzに３．０を代入すると、以下の式（３Ｃ）が得られる。
【数７】
【００７６】
したがって、図１０では、式（３Ｂ）で表される関数func_iter（ky，kz）を用いて、繰返し番号iter＝１〜iter_maxが割り当てられる。一方、図１１では、式（３Ｃ）で表される関数func_iter（ky，kz）を用いて、繰返し番号iter＝１〜iter_maxが割り当てられる。
以上のようにして、ビュー・オーダリングが行われる。
【００７７】
本形態では、データ収集領域Ｒａｃｑを分割するための曲線は、ｋｚ方向だけでなく、ｋｙ方向にも変化する。したがって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【００７８】
次に、本形態のビュー・オーダリングに従って被検体を撮影するときのフローについて説明する。尚、以下の説明では、図１５に示す３Ｄ
ＦＳＥのシーケンスを用いてデータを収集する場合のフローについて説明するが、本発明は、３Ｄ ＦＳＥのシーケンスの他に、３Ｄ ＳＳＦＰのシーケンス等の他の３Ｄシーケンスを用いてデータを収集する場合にも適用することができる。
【００７９】
図１２は、本形態のビュー・オーダリングに従って被検体を撮影するときのフローを示す図である。
【００８０】
ステップＳＴ１では、オペレータ１３が、図１５に示すシーケンスのエコー時間ＴＥおよびエコー数ｎ_maxを設定する。ここでは、オペレータ１３は、エコー数n_max＝１００に設定したとする。
【００８１】
ステップＳＴ２では、エコー番号決定手段９１（図１参照）が、エコー時間ＴＥおよびエコー数n_maxなどの値に基づいて、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ＝n_centerを決定する。ここでは、ｎ_center＝３９と決定されたとする。エコー番号ｎ_centerを決定したら、ステップＳＴ３に進む。
【００８２】
ステップＳＴ３では、関数決定手段９２（図１参照）が、先ず、ｋ空間の原点Ｃに、ステップＳＴ２で決定されたエコー番号n_center＝３９が割り当てられるように、式（１）の係数ratio_yzの値を設定する。係数ratio_yzの値を設定する方法の一例としては、エコー数n_maxの値と、エコー番号n_centerの値と、係数ratio_yzの値との対応関係を事前に作成しておく方法がある。本形態では、ステップＳＴ１において、オペレータ１３によってエコー数n_max＝１００が設定されており、ステップＳＴ２において、エコー番号n_center＝３９が決定されているので、上記の対応関係を事前に作成しておくことにより、係数ratio_yzの値を設定することができる。ここでは、係数ratio_yz＝１．０に設定されたとする。関数決定手段９２は、設定した係数ratio_yz＝１．０を式（１）に代入する。係数ratio_yz＝１．０を式（１）に代入することにより、データ収集領域を分割するために用いられる関数func_n（ky，kz）が決定される。係数ratio_yz＝１．０の場合、決定された関数func_n（ky，kz）は、式（１Ａ）で表される。関数func_n（ky，kz）を決定したら、ステップＳＴ４に進む。
【００８３】
ステップＳＴ４では、分割手段９３（図１参照）が、ステップＳＴ３で決定された関数func_n（ky，kz）を用いてデータ収集領域Ｒａｃｑのサンプリング点にエコー番号を割り当てる。そして、割り当てられたエコー番号に基づいて、図６に示すように、データ収集領域Ｒａｃｑに曲線ＣＬ１〜ＣＬｚを規定し、複数の領域Ｒ１〜Ｒｎ＿ｍａｘに分割する。分割したら、ステップＳＴ５に進む。
【００８４】
ステップＳＴ５では、関数決定手段９２が、ステップＳＴ３で設定した係数ratio_yz＝１．０を式（３）に代入する。係数ratio_yz＝１．０を式（３）に代入することにより、領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に繰返し番号を割り当てるために用いられる関数func_iter（ky，kz）が決定される。係数ratio_yz＝１．０の場合、決定された関数func_iter（ky，kz）は、式（３Ａ）で表される。関数func_iter（ky，kz）を決定したら、ステップＳＴ６に進む。
【００８５】
ステップＳＴ６では、ステップＳＴ５で決定された関数fucn_iterを用いて、領域Ｒ１〜Ｒｎ＿ｍａｘのサンプリング点に、番号iterが割り当てられる（図９参照）。番号iterを割り当てたら、ステップＳＴ７に進み、データを収集するためのスキャンを実行し、フローを終了する。
【００８６】
本形態では、データ収集領域Ｒａｃｑを分割するための曲線は、ｋｚ方向だけでなく、ｋｙ方向にも変化する。したがって、リンギングやブラーリングをｋｙ方向およびｋｚ方向の両方向に分散させることができる。
【００８７】
また、関数func_n（ky，kz）の基準点Ｐ（ky0，kz0）は、ｋ空間の原点Ｃからずれている。したがって、関数func_n（ky，kz）に含まれる係数ratio_yzの値を変えるだけで、ｋ空間の原点Ｃに割り当てられるエコー番号ｎ_centerを調整することができるので、コントラストを容易に調整することができる。
【００８８】
本形態では、データ収集領域Ｒａｃｑを分割する曲線は、基準点Ｐ（ky0，kz0）を基準にして規定されている（例えば、図６参照）。しかし、別の基準点を基準にして曲線を規定してもよい（図１３参照）。
【００８９】
図１３は、別の基準点を基準にした曲線でデータ収集領域Ｒａｃｑを分割したときの一例を示す図である。
【００９０】
図１３では、基準点Ｐ′を基準にして、データ収集領域Ｒａｃｑを分割する曲線ＣＬ１〜ＣＬｗを規定している。基準点Ｐ′は、図６に示す基準点Ｐ（ky0，kz0）よりも、ｋ空間の原点Ｃに近い位置に設定されている。このように、データ収集領域Ｒａｃｑを分割する曲線の基準点は、ｋ空間の原点Ｃとは異なる位置であれば、様々な位置に設定することが可能である。
【００９１】
また、本形態では、データ収集領域Ｒａｃｑを円又は楕円の曲線で分割しているが、別の曲線で分割してもよい（図１４参照）。
【００９２】
図１４は、円又は楕円の曲線とは別の曲線でデータ収集領域を分割する一例を示す図である。
【００９３】
図１４では、基準点Ｐ（ky0，kz0）を基準にした菱形の曲線でデータ収集領域Ｒａｃｑを分割した場合が示されている。このように、データ収集領域Ｒａｃｑを分割する曲線は、円、楕円、双曲線、放物線などのような曲線であってもよいし、菱形のように複数の直線の組合せで構成される曲線であってもよい。
【００９４】
また、本形態では、関数func_iter（ky，kz）は、サンプリング点の位置を角度で表しているが、角度以外の別のパラメータを用いてサンプリング点の位置を表してもよい。
【符号の説明】
【００９５】
２ 磁場発生装置
３ テーブル
４ 受信コイル
５ シーケンサ
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 中央処理装置
１０ 操作部
１１ 表示部
１２ 被検体
１３ オペレータ
２１ ボア
２２ 超伝導コイル
２３ 勾配コイル
２４ 送信コイル
３１ クレードル
９１ エコー番号決定手段
９２ 関数決定手段
９３ 分割手段
９４ 番号割当て手段
１００ ＭＲＩ装置
【特許請求の範囲】
【請求項１】
ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置であって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割手段、
を有する磁気共鳴イメージング装置。
【請求項２】
サンプリング点にエコー番号を割り当てるための第１の関数を決定する関数決定手段を有する、請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
ｋ空間の原点に割り当てられるエコー番号を決定するエコー番号決定手段を有し、
前記関数決定手段は、
１回のデータ収集シーケンスで収集されるエコー数と、前記エコー番号決定手段により決定されたエコー番号とに基づいて、前記第１の関数を決定する、請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記分割手段は、
前記第１の関数によってサンプリング点に割り当てられたエコー番号に基づいて、前記曲線を規定する、請求項２又は３に記載の磁気共鳴イメージング装置。
【請求項５】
前記関数決定手段は、
前記複数の領域のサンプリング点に前記データ収集シーケンスの繰返し番号を割り当てるための第２の関数を決定する、請求項２〜４のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項６】
前記第２の関数は、
前記基準点に対する前記サンプリング点の角度を表す関数である、請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記曲線は、円又は楕円である、請求項１〜６のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項８】
ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置のプログラムであって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割処理、
を計算機に実行させるためのプログラム。
【請求項１】
ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置であって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割手段、
を有する磁気共鳴イメージング装置。
【請求項２】
サンプリング点にエコー番号を割り当てるための第１の関数を決定する関数決定手段を有する、請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
ｋ空間の原点に割り当てられるエコー番号を決定するエコー番号決定手段を有し、
前記関数決定手段は、
１回のデータ収集シーケンスで収集されるエコー数と、前記エコー番号決定手段により決定されたエコー番号とに基づいて、前記第１の関数を決定する、請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記分割手段は、
前記第１の関数によってサンプリング点に割り当てられたエコー番号に基づいて、前記曲線を規定する、請求項２又は３に記載の磁気共鳴イメージング装置。
【請求項５】
前記関数決定手段は、
前記複数の領域のサンプリング点に前記データ収集シーケンスの繰返し番号を割り当てるための第２の関数を決定する、請求項２〜４のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項６】
前記第２の関数は、
前記基準点に対する前記サンプリング点の角度を表す関数である、請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記曲線は、円又は楕円である、請求項１〜６のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項８】
ｋｙ−ｋｚ面に規定されたデータ収集領域を複数の領域に分割し、前記複数の領域に配置されるエコーを収集するためのデータ収集シーケンスを繰り返し実行する磁気共鳴イメージング装置のプログラムであって、
ｋｙ−ｋｚ面の原点とは異なる点を基準にして規定される複数の曲線によって、前記データ収集領域を複数の領域に分割する分割処理、
を計算機に実行させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０】
【公開番号】特開２０１２−１１５４７８（Ｐ２０１２−１１５４７８Ａ）
【公開日】平成２４年６月２１日（２０１２．６．２１）
【国際特許分類】
【出願番号】特願２０１０−２６７８７３（Ｐ２０１０−２６７８７３）
【出願日】平成２２年１１月３０日（２０１０．１１．３０）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
【公開日】平成２４年６月２１日（２０１２．６．２１）
【国際特許分類】
【出願日】平成２２年１１月３０日（２０１０．１１．３０）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
[ Back to top ]