(57)【特許請求の範囲】
【請求項１】
ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する収集部と、
前記ｋ空間データに基づいて磁気共鳴画像を生成する生成部とを備え、
前記複数のストリップのそれぞれは、前記ｋ空間における一方の周縁部から他方の周辺部にわたる範囲をサンプリングするように配置され、
前記複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、前記ｋ空間の中心に対する前記アンダーサンプリング領域の反対側の位置がサンプリングされるように配置される、
磁気共鳴イメージング装置。
【請求項２】
前記複数のストリップは、複数のアンダーサンプリング領域が前記ｋ空間内に形成されるように配置され、かつ、当該複数のアンダーサンプリング領域それぞれについて前記ｋ空間の中心に対する反対側の位置がサンプリングされるように配置される、
請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記生成部は、前記複数のストリップそれぞれから収集されたｋ空間データを統合し、前記ｋ空間の中心に対する反対側の位置からサンプリングされたｋ空間データを用いた画像再構成を行うことで、前記磁気共鳴画像を生成する、
請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項４】
前記生成部は、前記磁気共鳴画像を生成する前に、前記ストリップごとに動き推定量を設定し、設定された動き推定量を用いて、前記ｋ空間データにおける位相及び回転のうち少なくとも一方を補正し、補正されたｋ空間データを統合する、
請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する収集部と、
前記ｋ空間データに基づいて磁気共鳴画像を生成する生成部とを備え、
前記複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、前記ｋ空間の中心に対する前記アンダーサンプリング領域の反対側の位置がサンプリングされるように配置され、
前記生成部は、前記磁気共鳴画像を生成する前に、オーバーサンプリングされた領域において、各ストリップに対して、ストリップ間のデータの相関関係の質、又は、面内の動きのペナルティに関する重みを付与し、付与された重みを用いて、前記複数のストリップそれぞれから収集されたｋ空間データを統合し、前記ｋ空間の中心に対する反対側の位置からサンプリングされたｋ空間データを用いた画像再構成を行うことで、前記磁気共鳴画像を生成する、
磁気共鳴イメージング装置。
【請求項６】
ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する収集部と、
前記ｋ空間データに基づいて磁気共鳴画像を生成する生成部とを備え、
前記複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、前記ｋ空間の中心に対する前記アンダーサンプリング領域の反対側の位置がサンプリングされるように配置され、
前記生成部は、前記磁気共鳴画像を生成する前に、オーバーサンプリングされた領域において、各ストリップに含まれる各データ収集ラインに対して、データ収集ライン間の相対的なエコーの質に関する重みを付与し、付与された重みを用いて、前記複数のストリップそれぞれから収集されたｋ空間データを統合し、前記ｋ空間の中心に対する反対側の位置からサンプリングされたｋ空間データを用いた画像再構成を行うことで、前記磁気共鳴画像を生成する、
磁気共鳴イメージング装置。
【請求項７】
前記各データ収集ラインに付与される重みは、前記オーバーサンプリングされた領域における２本以上のライン間のエコーのコントラストの変化に従って定められる、
請求項６に記載の磁気共鳴イメージング装置。
【請求項８】
ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する収集部と、
前記ｋ空間データに基づいて磁気共鳴画像を生成する生成部とを備え、
前記複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、前記ｋ空間の中心に対する前記アンダーサンプリング領域の反対側の位置がサンプリングされるように配置され、
前記複数のストリップのうちの少なくとも１つのストリップが、前記ｋ空間内にある所定の領域の周縁において、他のストリップと重なるように配置される、
磁気共鳴イメージング装置。
【請求項９】
前記他のストリップとの重なりの量は、動きの補正における許容範囲として設定される、請求項８に記載の磁気共鳴イメージング装置。
【請求項１０】
複数のストリップは、隣り合うストリップ間の角度が不均一な状態で前記ｋ空間内に放射状に配置される、
請求項１に記載の磁気共鳴イメージング装置。
【請求項１１】
前記複数のストリップのうちの少なくとも１つのストリップが、前記ｋ空間の中心に対して位置がずれるように配置される、
請求項１〜１０のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項１２】
前記複数のストリップそれぞれが、各ストリップの少なくとも短軸方向における前記中心からのオフセット分だけ位置がずれるように配置される、
請求項１１に記載の磁気共鳴イメージング装置。
【請求項１３】
前記アンダーサンプリング領域は、前記ｋ空間内にある略円形領域の周縁上に略均一に分散配置される、
請求項１〜１２のいずれか一つに記載の磁気共鳴イメージング装置。
【請求項１４】
前記収集部は、高速スピンエコー法のパルスシーケンスを実行することによって、前記ｋ空間データを収集する、
請求項１〜１３のいずれか一つに記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
従来、磁気共鳴イメージング装置によるデータ収集の方法として、ｋ空間内において放射状のストリップに沿ってデータを収集する放射状ストリップスキャン手法が知られている。例えば、ここでいう放射状ストリップスキャン手法の一つとして、ＰＲＯＰＥＬＬＥＲ（Periodically  Rotated  Overlapping  Parallel  Lines  with  Enhanced  Reconstruction）が知られている。また、このような放射状ストリップスキャン手法において、ｋ空間内をアンダーサンプリングすることによって、データ収集を高速化する方法が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】米国特許第７０２３２０７号明細書
【特許文献２】米国特許第７１２２３４８号明細書
【特許文献３】米国特許第７３８２１２７号明細書
【特許文献４】米国特許第７３７２２６９号明細書
【特許文献５】米国特許第７４８２８０６号明細書
【特許文献６】米国特許第７６１９４１０号明細書
【特許文献７】米国特許出願公開第２００９／０１２９６４８号明細書
【特許文献８】米国特許第７５３５２２２号明細書
【特許文献９】米国特許出願公開第２０１４／０１１１２０３号明細書
【特許文献１０】米国特許第７５３８５４９号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、アンダーサンプリングによって生じるアーチファクトを低減することができる磁気共鳴イメージング装置を提供することである。
【課題を解決するための手段】
【０００５】
実施形態に係る磁気共鳴イメージング（Magnetic  Resonance  Imaging：ＭＲＩ）装置は、収集部と、生成部とを備える。収集部は、ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する。生成部は、前記ｋ空間データに基づいて磁気共鳴画像を生成する。前記複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、前記ｋ空間の中心に対する前記アンダーサンプリング領域の反対側の位置がサンプリングされるように配置される。
【図面の簡単な説明】
【０００６】
【図１】図１は、１又は複数の実施形態に係る、アンダーサンプリングされたｋ空間を伴う、放射状ストリップスキャン手法（ＰＲＯＰＥＬＬＥＲ等）の高速化に適したＭＲＩ装置の大まかな概略ブロック図である。
【図２】図２は、従来のＰＲＯＰＥＬＬＥＲ手法によるブレードの配置パターンを示す図である。
【図３】図３は、従来の高速化されたＰＲＯＰＥＬＬＥＲ手法によるブレードの配置パターン及びアンダーサンプリング領域の分布を示す図である。
【図４】図４は、一部の実施形態に係るアンダーサンプリング領域のｋ空間の分布を示す図である。
【図５】図５は、一部の実施形態に係る、スキャンブレードが９個である構成例であり、アンダーサンプリング領域が図４に示すパターンに対応する例を示す図である。
【図６】図６は、一部の実施形態に係る、スキャンブレードが５個である別の構成例であり、アンダーサンプリング領域が図４に示すパターンに対応する例を示す図である。
【図７】図７は、一部の実施形態に係る、ブレードを重ね合わせてイメージング対象の回転によるアンダーサンプリングを抑制する手法を示す図である。
【図８】図８は、１又は複数の実施形態に係る、アンダーサンプリングされたｋ空間を伴う、放射状ストリップスキャン手法（例えば、ＰＲＯＰＥＬＬＥＲ等）を高速化するための処理のフローチャートを示す図である。
【図９】図９は、一部の実施形態に係る、図８の処理で使用されてよい、アンダーサンプリングされたｋ空間データを利用した、ＭＲＩ再構成のための処理を示すフローチャートである。
【発明を実施するための形態】
【０００７】
本発明の実施形態は、磁気共鳴イメージングにおいて、アンダーサンプリングされた（undersampled）ｋ空間を利用して、例えば、ＰＲＯＰＥＬＬＥＲ（Periodically  Rotated  Overlapping  Parallel  Lines  with  Enhanced  Reconstruction）等の放射状ストリップスキャン手法を高速化するための装置及び方法に関する。
【０００８】
本明細書に開示する本発明の実施形態例は、アンダーサンプリングされたｋ空間を利用して、放射状ストリップスキャン手法の高速化を実現する。実施形態において、ｋ空間データは、ｋ空間におけるアンダーサンプリング領域を最適化するパターンに従った高速化スキャンにより収集される。この時、アーチファクトがより少ない状態の画像再構成が行われるように収集される。
【０００９】
本明細書に開示する本発明の実施形態例は、デカルト軌道内ではなくｋ空間内において放射状ストリップに沿ってデータを収集する放射状ストリップスキャン手法（限定はされないが、ＰＲＯＰＥＬＬＥＲ又は別の手法等）を用いて収集した、アンダーサンプリングされたｋ空間データから得られた磁気共鳴（Magnetic  Resonance：ＭＲ）画像のスキャン速度や質を改善するように機能する。ＰＲＯＰＥＬＬＥＲは、多種類の動きアーチファクトを軽減するのに有利である。ＰＲＯＰＥＬＬＥＲ手法は、他の様々な名称でも知られており、例えば、ＪＥＴ（登録商標）（例えば、東芝製のＭＲＩ機器関連）、ＢＬＡＤＥ（登録商標）（例えば、シーメンス製のＭＲＩ機器関連）、Ｍｕｌｔｉ−Ｖａｎｅ（登録商標）（例えば、フィリップス製のＭＲＩ機器関連）、ＲＡＤＡＲ（登録商標）（例えば、日立製のＭＲＩ機器関連）がある。
【００１０】
ＰＲＯＰＥＬＬＥＲでは、放射状に配置された１組のストリップを回転させることでサンプリングｋ空間を構成する。各放射状ストリップは、所定の高周波（Radio  Frequency：ＲＦ）パルスエコーシーケンスに従って読み出される複数の並列位相符号化ラインを備える。ＰＲＯＰＥＬＬＥＲは、第１の撮影で、所定数の放射状ストリップについてｋ空間データを収集するステップと、所定の回転角度で上記１組の放射状ストリップを回転させるステップと、回転させた１組の放射状ストリップを用いて第２の撮影を行うステップと、完全な（又は略完全な）ｋ空間が得られるまでこの処理を繰り返すステップとによって機能する。この場合、回転する１組の放射状のブレードで画定される中心（例えば、ｋ空間の中心）を持つ円の形をとる領域がサンプリングされる。
【００１１】
ＰＲＯＰＥＬＬＥＲでの軌道は完全な放射状ではない。それは、その軌道ライン同士が、角度及び中心点からの直角に対するオフセットによって違うからである。ＰＲＯＰＥＬＬＥＲシーケンスの具体的な態様、その収集、及びその再構成は、角度によるステップ、又は角度を持った１つのブレード内のデカルト位相符号の変位に強く関連する。したがって、ＰＲＯＰＥＬＬＥＲを「放射状」と呼ぶことは、ＰＲＯＰＥＬＬＥＲの詳細についての記述として理想的であるとは限らない。しかし、ほぼ円形のｋ空間領域の周縁を取り巻く十分な対象範囲を与えるのに必要な相異なるライン本数を考えると、本明細書においては、用語「放射状」が、そのような対象範囲の一般的な記述子として、特にｋ空間領域の周縁において使用される。「放射状」という用語が本明細書で使用されているのは、以下を強調するためである。ｋ空間のラインは、一般的な「デカルト位相で符号化されたｋ空間」とは対照的に、収集されたｋ空間の周縁に沿った位置ならどこで終わってもよい。「デカルト位相で符号化されたｋ空間」では、収集されたデータのラインが終わるのは、矩形領域の２つの対向する辺であって他の２つの辺ではないというのが一般的である。本明細書においては、軌道が中心点回りの角度位置による回転という意味でしか変化しないという「放射状」という用語は、他の意味では使用されない。
【００１２】
ＰＲＯＰＥＬＬＥＲは、多くの場合、高速スピンエコー法（Fast  Spin  Echo：ＦＳＥ）のシーケンスを基本として行われる。この技術が画像品質に関し極めて有利である臨床的な状況としては、非協力的な患者（例えば、子どもや認知症の成人）、動きの激しい解剖学的部位（例えば、肝臓や胎児イメージング）、本来的な動き感度が高いシーケンス（例えば、拡散強調イメージング）、一般的に動きアーチファクトが顕著な部位（例えば、頸椎）、又は磁化率の歪みが許容不可能な部位（例えば、高分解能の拡散強調エコープラナーイメージング法）等が挙げられる。ＰＲＯＰＥＬＬＥＲは、スキャン対象のある種の動きには実質的に影響されない。なぜなら、各繰り返し時間（Repetition  Time：ＴＲ）を、特にその時間についての位置ずれの量に従って、おそらくは個々のブレードを補正して平均的な位置へ可能な限り近づけることにより、収集を長くするように補正できるからである。
【００１３】
ＰＲＯＰＥＬＬＥＲシーケンスの中で、エコーシーケンスは、少なくとも部分的にはＲＦリフォーカスされたスピンエコーの生成に依存する、ＦＳＥのエコートレイン又はＧＲＡＳＥ（Gradient  And  Spin  Echo）法のエコートレインを使用する場合がほとんどである。しかしながら、ＰＲＯＰＥＬＬＥＲシーケンスの例では、放射状に方向づけられたブレード内のラインが、勾配リフォーカスされたエコーを用いて収集される例も知られている。
【００１４】
ＰＲＯＰＥＬＬＥＲの最大のデメリットは、空間分解能が同じ場合に他の方法より遅いことである。これは、ＰＲＯＰＥＬＬＥＲでは、本来、ｋ空間の中心がオーバーサンプリング（oversampling）されるという事実に関係する。このオーバーサンプリングが、ＰＲＯＰＥＬＬＥＲの動き特性が良好である主な理由である。それにより、少なくともｋ空間中央部において、ある種の動きの平均化が可能になる。さらに、連続するＴＲから得られるｋ空間中央部におけるデータの比較が可能になり、それにより、動きの検出及び数種類の動き（例えば、何らかの回転又は平行移動）の補正が可能になる。
【００１５】
ＰＲＯＰＥＬＬＥＲは、解剖学上の位置ずれの総計（特に２次元マルチスライス収集について平面内の動き）を補正するのに使用できる。ＰＲＯＰＥＬＬＥＲは、高分解能のために多数の撮影が必要な場合に、データの一貫性の改善にも寄与できる。例えば、通常、拡散強調イメージングでは、撮影ごとに信号位相に大きな変動が生じるが、その場合に、ＰＲＯＰＥＬＬＥＲの利用により一貫性の程度が高くなる可能性がある。ＰＲＯＰＥＬＬＥＲには、他にも、オフレゾナンスによる比較的無害なアーチファクトをはじめとする利点がある。上記のようにＰＲＯＰＥＬＬＥＲには臨床的に利用価値のある利点が多くあるが、放射状の収集では、患者の快適さ等の要因により収集時間を短縮するために何らかの形で角度上のアンダーサンプリングを伴うことが多い。ＰＲＯＰＥＬＬＥＲは、これを呈示できる。ＰＲＯＰＥＬＬＥＲでは、ｋ空間領域の周縁部分はアンダーサンプリングされる可能性がある一方で、収集されたｋ空間のそれ以外の内側部分は、オーバーサンプリングされるか、又は臨界サンプリングされるかのいずれかである可能性がある。
【００１６】
ＭＲＩシーケンスでは、ｋ空間内のサンプルは、対象画像の画像再構成領域（Field  Of  View：ＦＯＶ）について、少なくともある方向に沿って、ナイキスト（Nyquist）基準に対して相対的に定まるサンプル密度に従ってアンダーサンプリング、オーバーサンプリング、又は臨界サンプリングされる。しかしながら、ＰＲＯＰＥＬＬＥＲでは、このサンプリングは、別の有用な方策として説明できる。通常通りに、１つのブレードにより、そのブレードに対応するｋ空間内領域について、ｋ空間内の臨界サンプリングが提供される。そこで、ｋ空間内の任意の位置は、ｋ空間におけるその近傍のデータに寄与するブレードの数によって、アンダーサンプリングされる、オーバーサンプリングされる、又は臨界サンプリングされる、のいずれかに類型化することができる。そのような基準により、ｋ空間内の小領域は、その小領域に関わるデータを提供するブレードが２つ以上ある場合にオーバーサンプリングされる。小領域は、１ブレードのみがそのデータに寄与する場合に、臨界サンプリングされる。データの小領域は、その領域についてのｋ空間データを収集するブレードがない場合に、アンダーサンプリングされる（又は「サンプリングされない／非サンプリングとなる（unsampled）」）。このように、オーバーサンプリングは、ｋ空間内の複数のブレードによるブレード対応領域の重なり等に関連する。
【００１７】
図１に示すＭＲＩ装置は、ガントリ１０（概略断面で示す）と、これに接続された各種の関連システム構成要素２０とを有する。少なくともガントリ１０は、通常はシールドルーム内に配置される。図１に示すＭＲＩ装置のジオメトリは、実質的に同軸の円筒形に配置された静磁場Ｂ０磁石１２と、Ｇｘ、Ｇｙ、及びＧｚの傾斜磁場コイルセット１４と、大型の全身用ＲＦコイル（Whole  Body  RF  Coil：ＷＢＣ）アセンブリ１６とを有する。この円筒形に配置される要素の横軸に沿って、患者寝台１１により支持される患者９の関心解剖部位を実質的に取り囲むものとして、イメージングボリューム１８が示される。
【００１８】
実施形態によっては、より小型のアレイＲＦコイル（図示せず）を患者の一部分に接続し、イメージングボリュームが、この小型のアレイＲＦコイルで取り囲まれる領域を含むようにしてもよい。当業者には明らかなように、ＷＢＣと比較して、複数の、表面コイル等のような小型コイルやアレイは、特定の身体部分（例えば、腕、肩、肘、手首、膝、脚、胸、背骨等）に合わせて設計されることが多い。本明細書では、そのような小型のＲＦコイルを、アレイコイル（Array  Coil：ＡＣ）又はフェーズドアレイコイル（Phased  Array  Coil：ＰＡＣ）と呼ぶ。それらの複数のコイルには、イメージングボリューム内にＲＦ信号を送信するように構成される少なくとも１個のコイルと、イメージングボリュームからＲＦ信号を受信するように構成される１個以上の受信コイルとを含んでもよい。
【００１９】
ＭＲＩシステムコントローラ２２は、ディスプレイ２４、キーボード２６、及びプリンタ２８に接続された入出力ポートを有する。当然のことながら、ディスプレイ２４は、制御入力もできるようにタッチスクリーンタイプのものでもよい。
【００２０】
ＭＲＩシステムコントローラ２２又はＭＲＩシステムコントローラ２２に連結されたコンピュータを操作して、パルスシーケンスに関する情報をＭＲＩシーケンスコントローラ３０に供給したり、インストールされたソフトウェアプログラムに従って、システム全体の処理を管理したりしてもよい。また、ＭＲＩシステムコントローラ２２は、自動音声合成技術を用いて生成された音声メッセージにより、例えば息止め等のタスクを行うように患者に指示する要素として機能してもよい。
【００２１】
ＭＲＩシステムコントローラ２２はＭＲＩシーケンスコントローラ３０に接続され、ＭＲＩシーケンスコントローラ３０は、Ｇｘ、Ｇｙ、及びＧｚの傾斜磁場コイルドライバ３２、並びにＲＦ送信機３４及び送受信スイッチ３６（同じＲＦコイルが送信と受信の両方に使用される場合）を制御する。ＭＲＩシーケンスコントローラ３０は、ＭＲＩイメージング（核磁気共鳴（Nuclear  Magnetic  Resonance：ＮＭＲ）イメージングとしても知られている）技術を実装するための適切なプログラムコード構造３８を含む。ＭＲＩイメージング技術として、例えば、パラレルイメージング等のイメージングシーケンスがある。
【００２２】
パルスシーケンス情報は、パルスシーケンスに従ってＧｘ、Ｇｙ、及びＧｚの傾斜コイルドライバ３２、並びにＲＦ送信器３４及び送受信スイッチ３６（同じＲＦコイルが送信と受信の両方に使用される場合）を操作するのに必要とされる情報を含む。そのような情報は、ｘコイル、ｙコイル、及びｚコイルに印加されるパルス電流の強度、持続時間、及び印加タイミングを含む。また、デューティサイクルすなわちＴＲは、ＭＲＩシステムコントローラ２２とＭＲＩシーケンスコントローラ３０のうち少なくとも一方により制御される。パルスシーケンスとデューティサイクル等の送信は、例えば、ＰＲＯＰＥＬＬＥＲに限定はされないが、ＰＲＯＰＥＬＬＥＲ等の放射状収集等の信号収集パターンに関するプログラムコードに基づいて、ＭＲＩシステムコントローラ２２、ＭＲＩシーケンスコントローラ３０、及びＭＲＩデータプロセッサ４２のうち少なくとも１つにより制御してもよい。
【００２３】
ＭＲＩ装置の関連システム構成要素２０は、ＭＲＩデータプロセッサ４２に入力を与えるＲＦ受信機４０を含む。ＭＲＩデータプロセッサ４２はその入力を使用して、ディスプレイ２４に送られる処理画像データを生成する。また、ＭＲＩデータプロセッサ４２は、システム構成パラメータを含むメモリ４６と、プログラムコード構造４４と、プログラム／データ記憶装置５０とにアクセスするように構成される。また、ＭＲＩプログラムコード構造４４及び５０は、磁気共鳴（Magnetic  Resonance：ＭＲ）画像を再構成するための制御ロジックに加え、全身用ＲＦコイル１６やその他の受信コイルからＭＲデータを収集するための制御ロジックを含んでもよい。例えば、図８及び図９に関連して後述する処理８００及び処理９００等の処理は、ＭＲＩシステムコントローラ２２、ＭＲＩシーケンスコントローラ３０、及びＭＲＩデータプロセッサ４２のうち少なくとも１つにより実行できる。
【００２４】
図１では、ＲＦ送信機３４、送受信スイッチ３６、及びＲＦ受信機４０がＲＦコイル１６から離れて位置するように図示するが、実施形態によっては、それらのいずれも、ＲＦコイル１６やその他のＲＦコイルのごく近く又は表面に位置してもよい。
【００２５】
また図１に、ＭＲＩ装置のプログラム／データ記憶装置５０の一般的な説明を示す。ここで、プログラム／データ記憶装置５０においては、格納されたプログラムコード構造、及び構成済み又は所定のデータが、ＭＲＩ装置の種々のデータ処理構成要素にアクセスできる非一時的コンピュータ可読記憶媒体に格納される。（上記の格納されたプログラムコード構造は、例えば、放射状ストリップスキャン手法に従ったデータ収集を制御するため、アンダーサンプリングされたｋ空間データ収集時の画像再構成改善のため、及び、操作者の入力を受け付けるグラフィカルユーザーインタフェース（Graphical  User  Interface：ＧＵＩ）を定義するための構造である。また、上記の構成済み又は所定のデータは、例えば、プログラムの実行を制御するための特定の閾値設定、及び、圧縮センシングによるコイル感度及び画像再構成の同時推定を制御するためのパラメータである。）当業者には明らかなように、プログラム／データ記憶装置５０をセグメント化して、少なくとも一部分を、関連システム構成要素２０の処理コンピュータのうち、通常操作においてその格納されたプログラムコード構造を最優先で必要とする別のコンピュータに直接接続してもよい（すなわち、共通の場所に格納してＭＲＩシステムコントローラ２２に直接接続するのではない）。
【００２６】
実際に、当業者には明らかなように、図１は、後述する例示的実施形態を実現するために変更された典型的なＭＲＩ装置の、非常に大まかな概略図を示したものである。システム構成要素は様々な論理集合の「ボックス」に分割することができ、通常、多数のデジタル信号プロセッサ（Digital  Signal  Processor：ＤＳＰ）と、マイクロプロセッサと、専用処理回路（例えば、高速アナログディジタル（Analog-Digital：ＡＤ）変換用、高速フーリエ変換用、アレイ処理用等）とを含む。これらのプロセッサの各々は、通常、クロック制御された「状態マシン」であり、物理データ処理回路は、クロックサイクル（又は、所定数のクロックサイクル）毎に、ある物理状態から別の物理状態に移る。
【００２７】
処理回路（例えば、中央演算処理装置（Central  Processing  Unit：ＣＰＵ）、レジスタ、バッファ、演算装置等）の物理状態が、操作過程において、あるクロックサイクルから別のクロックサイクルに徐々に変化するだけでなく、関連データ記憶媒体（例えば、磁気記憶媒体内のビット記憶場所）の物理状態も、このような装置の操作過程において、ある状態から別の状態に変換される。例えば、画像再構成処理や時としてコイル感度マップ生成処理の終わりに、物理的記憶媒体内のコンピュータ可読でアクセス可能なデータ値の記憶場所の配列は、ある先行状態（例えば、全て一様に「０」値、又は全て「１」値）から新しい状態に変換され、そのような配列における物理的場所の物理状態は、最小値と最大値との間で変化して、実世界の物理的事象及び物理的条件（例えば、イメージングボリューム空間内の患者の内部物理構造）を表す。当業者には明らかなように、そのような格納データ値の配列は、ＭＲＩ装置内で特定のシーケンスの動作状態を引き起こし遷移させる特定構造のコンピュータ制御プログラムコードと同様に、命令レジスタに順次読み込まれＭＲＩ装置の関連システム構成要素２０の１つ以上のＣＰＵによって実行されると物理的構造を表し構成する。
【００２８】
図２は、従来のＰＲＯＰＥＬＬＥＲスキャンで使用される例示的なデータ収集パターン２００を示す図である。この図は、患者の脳の軸方向スライスのスキャンで収集されたデータが存在するｋ空間に対応してもよい、ｋ空間の略円形領域２０２を概略的に示す。図２に示す従来の手法では、ＰＲＯＰＥＬＬＥＲストリップ（本明細書において、「ブレード」、「放射状ストリップ」、又は「ストリップ」と呼ぶこともある）は、画像に変換されるｋ空間内の略円形領域が、隣り合うブレードの間に有角形状の（angular）（例えば、楔形の）アンダーサンプリング領域を含まないように配置される。例えば、ストリップは、１個のストリップの上部右側端部が、上記円形領域の周縁上又は非常に近くにある別の１個のストリップの上部左側端部と同一平面上にあるように配置される。それにより、直接隣り合うストリップのどのペアの間にも有角形状の隙間部（例えば、「有角形状のアンダーサンプリング領域」と呼んだ上記領域）が存在しないようになる。長方形のストリップの各短端と略円形領域２０２の周縁との間に非サンプリング領域があるが、それらの領域は、比較的小さな領域であること、及び位置がｋ空間中心から遠いこと、さらに、略円形領域２０２の形成自体がどちらかというと人工的なもので非本質的であることから、重要ではないと考えられる。正多角形等の比較的整ったｋ空間の形におけるデータ収集が、根本的に間違っているというわけではない。図示の例では、ｋ空間の略円形領域２０２のイメージングに合計４個のストリップが使用されている。ストリップ２０４、２０６、２０８、及び２１０のそれぞれは、四角形である。各ストリップには、そのストリップの長軸に平行な複数のデータ収集ライン（例えば、ライン２１２）が含まれる。例えば、ストリップ２０６に対応するｋ空間データを収集するために、ストリップ２０６の長さ方向に沿ったライン２１２等の複数の互いに平行なラインに沿ってデータが収集される。典型的な実施態様では、ＦＳＥの第１の「撮影」すなわちＴＲで、ストリップ２０４のデータが収集され、第２の「撮影」すなわちＴＲで、ストリップ２０６のデータが収集され、さらに同様に続く。ストリップについてのデータ収集は、そのストリップにおける互いに平行なデータ収集ラインのそれぞれについてデータを収集するステップを含む。
【００２９】
従来技術では、図３に示すようなストリップ配置を用いて図２に示すＰＲＯＰＥＬＬＥＲ手法を高速化するようにしていた。パターン３００では、４個のストリップ３０４、３０６、３０８、及び３１０が配置されており、直接隣り合うストリップの間のそれぞれに、同一の有角形状の隙間部が存在するようになっている。有角形状の隙間部、すなわち有角形状の非サンプリング領域はそれぞれ、ｋ空間内にあって、ストリップに沿って領域３０２がスキャンされる時にサンプリングされない領域である。例えば、図３に示すパターン３００では、データが収集される必要のあるｋ空間内領域をこのように縮小してデータを収集するために、有角形状の非サンプリング領域はサンプリングされないので、データ収集速度は改善される（例えば、速くなる）。
【００３０】
図示の例では、４個のストリップ３０４、３０６、３０８、及び３１０が配置されて、８個の有角形状の非サンプリング領域３１２〜３２６が形成される。違う個数のストリップを用いる場合は、有角形状の非サンプリング領域の個数も異なる。例えば、高速化方法でｎ個のストリップが略円形領域３０２のスキャンに使用される場合、偶数２ｎ個の有角形状のアンダーサンプリング領域が形成される。ｎ個のストリップを用いて２ｎ個の有角形状の非サンプリング領域が角度３６０度にわたって形成されて、これらのストリップがｋ空間円の中心を中心として対称になっている場合、有角形状の非サンプリング領域が角度位置αにあるのであれば、別の有角形状の非サンプリング領域が角度位置α＋１８０度に必ずある。例えば、楔形３１２及び３２０は、互いに直径上反対側である（すなわち、スキャンされる円形領域の直径の両側にある）。言い方を変えると、楔形３１２及び３２０は、１８０度離れている。
【００３１】
しかしながら、そのようなＰＲＯＰＥＬＬＥＲを高速化する従来技術では、ＭＲ画像では許容不可能なことが多いストリーク（streak）アーチファクト（例えば、ストリーキングと呼ばれる）が生じる。ストリークアーチファクトが生じる理由は、有角形状の非サンプリング領域が別の楔形の非サンプリング領域に対して正反対（例えば、直径上反対側）の位置にあることである。ある非サンプリング領域に対して直径上反対の位置に別の非サンプリング領域がある場合、これらの非サンプリング領域は、ｋ空間の対称性に頼って正確な位置に入れ込むことができず、結果として、最終的な画像にストリークアーチファクトが形成される。
【００３２】
ＭＲＩでは、欠損データに対処するための様々な再構成アルゴリズムが広く知られている。最も簡易な方法であって解析が最も容易な方法は、対応する測定されたサンプルがない、あるいは、少なくとも、スキャンのためのナイキストの臨界サンプリングによる空間形成の程度において近傍にサンプルがないようなｋ空間内位置のデータにゼロ値を付与するものである。この方法は、「ゼロフィリング（zero-filling）再構成」と呼ばれる場合が多い。そのような方法は、簡易ではあるが、大きくなることもある点拡がりのアーチファクトを取り込むことになる。
【００３３】
より優れた、より効果的な再構成方法は、概して、ｋ空間の直径上反対の領域から収集されたデータを利用するものである。その収集データは、直接的に利用される場合も間接的に利用される場合もある。ただし、欠損データが反対側のデータでは補強できない場合、アーチファクトが生じる。これらのアーチファクトは、点拡がり関数（point  spread  function）の性質を帯びることが多く、その点拡がり関数は、根本的には、利用可能なｋ空間データ（収集されたものと数学的に合成されたものを含むデータ）の領域（又は形状）のフーリエ変換である。円のような非常に単純な定型の形状と関連する点拡がり関数であれば、ほとんど無害なアーチファクト、あるいは、まったくアーチファクトであると判断されないような点の広がりを与える。それに比べて、より複雑な形状と関連する点拡がり関数は、概して、許容不可能な度合いが相当高い。ｋ空間内の欠損した楔形が原因で、対応する再構成画像中に、離間した、多くの場合強烈であるストリーク様の構造がいかに生成されるか等である。
【００３４】
したがって、様々な従来の高速化技術によっても、ＰＲＯＰＥＬＬＥＲは速度が上がらないままで、少なくとも、一般的に利用されている２次元フーリエ変換（Two-Dimensional  Fourier  Transform：２ＤＦＴ）デカルト位相エンコーディング方法等の他のｋ空間軌道の多くと比較すると速度が遅いままである。この相対的な遅さの理由の一部は、位相エンコーディングの低減及び非対称の分解能、あるいはＦＯＶよりも外側の信号によるエイリアスの除去に関するシーケンスの柔軟性の低下、もしくは、矩形のＦＯＶ画像を得る手法の効率性又は互いに平行なイメージングコイルパターンへの最適な空間方向感度を最大限に利用する能力の低下に帰することができる。また、この相対的な遅さの一因は、提案されているＰＲＯＰＥＬＬＥＲ高速化の多くにおいては、その高速化がデカルト法に適用でき、デカルト法の場合に同程度に、又はより好ましく機能することである。
【００３５】
さらに、この遅さの理由の一部は、スキャンシーケンスにおけるある種の動きの補正能力の維持に必要な実務的配慮を行うせいである可能性もある（例えば、任意のストリップ上の部分フーリエ次数は、拡散強調スキャンにおける回転性の動きに対処するために、非効率さという点で漸進的にコストを増す、等）。いくつかの高速化手法には、追加的な処理及び推定が取り入れられているが、それは、おそらくは追加的なアーチファクト、動きエラーの除去性能の低下、ＳＮ比（Signal  to  Noise  Ratio：ＳＮＲ）の損失等という結果になりうる。一般的には、多くのＭＲ物理学者の考えによると、放射状収集では、そのｋ空間の同等の正方形をカバーするデカルト位相エンコードによるスキャンにおける平行収集ライン本数に比べて、ｋ空間の円の周縁をカバーするにはπ／２（約１．５７）倍の角度によるステップが必要とされる。この点は、正確に同一条件での比較ではないかもしれないが、放射状あるいは放射状に近いｋ空間パターンの課題を示す。
【００３６】
したがって、ＰＲＯＰＥＬＬＥＲは多種類の動き補正に好適である一方で、様々な従来の高速化手法によって改良された場合でも、非ＰＲＯＰＥＬＬＥＲスキャン（多くの場合、同様の方法、又はさらに多くの方法で高速化された非ＰＲＯＰＥＬＬＥＲスキャン）に対してスキャン時間の点では好適でないという課題が残る。同じ分解能に対してスキャン時間が長いことは、依然として、ＰＲＯＰＥＬＬＥＲの臨床使用に向けての障害である。
【００３７】
本実施形態に係るＭＲＩ装置は、収集部と、生成部とを備える。収集部は、ｋ空間に配置された複数のストリップで表されるデータ収集パターンに沿ってｋ空間データを収集する。生成部は、収集されたｋ空間データに基づいて磁気共鳴画像を生成する。複数のストリップは、少なくとも１つのアンダーサンプリング領域を間に挟むように配置され、かつ、ｋ空間の中心に対するアンダーサンプリング領域の反対側の位置がサンプリングされるように配置される。例えば、収集部は、ＭＲＩシーケンスコントローラ３０によって実現され、生成部は、ＭＲＩデータプロセッサ４２によって実現される。
【００３８】
ここで、複数のストリップは、複数のアンダーサンプリング領域がｋ空間内に形成されるように配置され、かつ、当該複数のアンダーサンプリング領域それぞれについてｋ空間の中心に対する反対側の位置がサンプリングされるように配置される。
【００３９】
また、生成部は、複数のストリップそれぞれから収集されたｋ空間データを統合し、ｋ空間の中心に対する反対側の位置からサンプリングされたｋ空間データを用いた画像再構成を行うことで、磁気共鳴画像を生成する。
【００４０】
本明細書において開示される実施形態では、スキャン領域の放射状ストリップ間の非サンプリング部分の存在を許しながらも、アンダーサンプリングによるアーチファクトを除去するための機能改善を、特にＰＯＣＳ（Projection  Onto  Convex  Sets）再構成等の部分フーリエ手法の使用について提供する。それらの例示的実施形態では、データがサンプリングされない有角形状の隙間部、すなわち「楔形状」を伴う放射状の収集領域によってデータを収集するが、これらの「楔形状」は、ｋ空間中心を挟んで正反対に位置する欠損データから成る、直径上反対側にある対を作らないように配置される。データがｋ空間の一方側で欠損していても、ｋ空間中心を挟んで反映されるｋ空間の対応する位置で対応するデータが収集可能であれば、欠損しているｋ空間データを正確に「充填する」機能が大幅に改善されるので、再構成補正の実行可能性が向上する。
【００４１】
ＰＲＯＰＥＬＬＥＲの性能を改善するために、様々な従来の技術において、各ブレードのサイズや構成（部分フーリエ、自動キャリブレーションを伴うパラレルイメージング等）、それぞれのブレードからのデータを収集する好ましい順序について研究されてきており、成功具合はまちまちである。しかしながら、従来の技術では、本明細書の例示的実施形態にあるようなデータ収集、すなわち、それぞれのブレードの移動が違ったり、それぞれのブレードの角形成が不均一であったり、又は隙間の位置が最適になっていたりするようなデータ収集は提供されていない。
【００４２】
それらの例示的実施形態によれば、患者のスキャン対象部位に対応するｋ空間内の円形領域の周縁に、奇数個の有角形状の（例えば、楔型の）非サンプリング領域が配置されるデータ収集パターンが設定される。例えば、アンダーサンプリング領域は、ｋ空間内にある略円形領域の周縁上に略均一に分散配置される。図４に示す実施形態のパターン４００では、９個の有角形状の非サンプリング領域４１２〜４２８が、ｋ空間内の円形領域４０２の周縁上に均一に分散配置される。図４は、非サンプリング（本明細書では「アンダーサンプリング」とも記載される）領域が、有角形状のアンダーサンプリング領域と別の有角形状のアンダーサンプリング領域とが直径上反対の場所に位置することがないように円形領域４０２の周縁に配されるようにする場合の、一部の実施形態に係る一般的な配置のしかたを概略的に示す。有角形状のアンダーサンプリング領域のそれぞれは、隣接する放射状ストリップの間の隙間部である。
【００４３】
なお、図４はやや様式化されている。もし、例えば、各ストリップが「直径」のように円形のｋ空間全体にわたる範囲をカバーする場合、あるいは、ストリップが単純な矩形形状の場合は、サンプリングパターンを図示通りに集めるのは実行不可能かもしれない。この点は、次に、図５及び図６にあるように現実的で実行可能なサンプリングパターンを用いて検討する。図４は、楔形の隙間部の基本的な位置を示すための図であり、楔形の正確な形状あるいはストリップの正しい配置を示すためのものではない。
【００４４】
一部の実施形態に従って配置される場合、少なくとも１個の有角形状非サンプリング領域に対しては、別の１個の有角形状非サンプリング領域は、これら２個の有角形状非サンプリング領域の位置がｋ空間中心を対称点として鏡映対称になるような位置にあることはない。一部の実施形態の有角形状非サンプリング領域には、別の１個の有角形状非サンプリング領域が、任意の２個の有角形状アンダーサンプリング領域の位置がｋ空間中心を対称点として鏡映対称になるような位置にあるものはない。複数の有角形状の非サンプリング領域を患者の１スキャンスライスに対応するｋ空間領域の周縁に沿って配置するに際し、どの有角形状の非サンプリング領域も別の有角形状の非サンプリング領域に対して直径上反対の位置にはないようにした場合、ｋ空間データの欠損が基本的にないようにｋ空間情報が導出されて、必要に応じて充填される。それにより、そのスライスの、実質的にアーチファクトのない画像の再構成が可能になる。この実質的に完全なｋ空間データ再構成は、ｋ空間のエルミート対称の性質によって可能になる。画像域でほとんどの複素位相を除去するために多数の方法があるが、必要に応じて、再構成にそのいずれかを利用してもよく、当該技術分野で広く知られている通り、それによりデータは位相除去後に本質的にエルミートとなる。
【００４５】
実施形態では、エコートレイン長（Echo  Train  Length：ＥＴＬ）を減らすために意図的にアンダーサンプリングの区域を導入しながら各放射状ストリップにおいて使用するライン本数を少なくすることや、使用する放射状ストリップ数を少なくすることで、ＰＲＯＰＥＬＬＥＲ又はＰＲＯＰＥＬＬＥＲ類似の手法でスキャンを高速化する。
【００４６】
図５は、一部の実施形態に係る、データ収集パターンを決定する手法を示す図である。アンダーサンプリング領域の望ましい位置を設定した後に、データ収集の放射状ストリップが、アンダーサンプリング領域間のスペースを充填するように設定される。図５の手法によれば、望ましいアンダーサンプリング領域の数に等しい数の放射状ストリップ５０４、５０６、５０８、５１０、５３４、５３６、５３８、５４０、５４２が設定され、その読み出し方向は長さ方向である。各放射状ストリップはより細くしてもよい。図示の通り、アンダーサンプリング領域５１２、５１４、５１６、５１８、５２０、５２２、５２４、５２６、５２８は互いに直径上反対の位置にはない一方で、有角形状のアンダーサンプリング領域は、この例では、略円形のｋ空間領域５０２の周縁上に均一に分散配置される。特にこの一組の制約は、もし選択されれば、ストリップ数は必ず奇数となることを意味する。
【００４７】
例えば、収集部は、ＦＳＥ法のパルスシーケンスを実行することによって、ｋ空間データを収集する。
【００４８】
ある典型的な実施形態では、放射状ストリップ５０４に関するデータは第１のＴＲで収集され、放射状ストリップ５０６に関するデータは第２のＴＲで収集され、というように続き、１ＴＲにつき１ストリップが対応する。しかしながら、複数のＴＲからのデータが１個のストリップに関与することも可能であり、反対に、１つのＴＲに複数のストリップについてのデータが提供されることも可能であると予期される。
【００４９】
一部の実施形態では放射状ストリップ（例えば、図５に図示するような）は大きさが同じであるが、他の実施形態では放射状ストリップは大きさが不均一である。放射状ストリップの大きさは、各ストリップにおける収集データの、ひいてはエコー数に関係するライン本数（例えば、Ｔ２−ＦＳＥ法のスキャン用データの、９〜２５本のライン）によって決定してもよい。ライン間の間隔の取り方は、所望のＦＯＶに数学的に関係する。矩形画像ＦＯＶであることが理由で、ストリップの向きに応じて幾何形状にばらつきを持たせることも考えられる。同様に、アレイコイルが様々な方向においてパラレルイメージング性能に違いがあることも、角度の関数としてのストリップの幾何形状にばらつきを持たせる理由としてもよい。
【００５０】
一般的には、少なくとも実施態様を簡単にするという理由で、１回のスキャンにおけるすべての放射状ストリップは同じ大きさであるのが望ましい。
【００５１】
図６は、一部の実施形態に係る、データ収集パターンを決定する別の手法を示す図である。図６に図示する手法では、図５の手法で使用されたものより幅広のストリップを用いて、やや多めにデータを収集することもできる。図６に図示する通り、一例示的実施形態では、９個の有角形状アンダーサンプリング領域に対応するために５個の大きなブレードを必要とするデータ収集パターンを使用する。図６に図示する大きなブレードのそれぞれは、ｋ空間のうち、隣り合う有角形状アンダーサンプリング領域の１つから別の１つまでの領域全体を範囲とする。
【００５２】
例えば、複数のストリップのうちの少なくとも１つのストリップが、ｋ空間内にある所定の領域の周縁において、他のストリップと重なるように配置される。このとき、例えば、他のストリップとの重なりの量は、動きの補正における許容範囲として設定される。
【００５３】
図６に図示するように、円形領域の周縁における一部のｋ空間領域は、複数の放射状ストリップでカバーしてもよい。例えば、放射状ストリップ６１２及び６１４は、有角形状アンダーサンプリング領域６２２及び６２４の間において、放射状ストリップパターン６００に関連付けられた円形領域の周縁において重なり、放射状ストリップ６１２及び６１４のそれぞれについて１回ずつ収集が行われると、その部分の収集は２回行われることになる。このように、図６に図示する実施形態によると、データ収集の対象範囲は一様ではない。また、（２ｎ＋１）個の有角形状のアンダーサンプリング領域に対して、この手法では、（ｎ＋１）個のサンプリング対象ストリップで収集を行う。
【００５４】
一部の実施形態では、ＥＴＬを短縮するために各ストリップで使用するスキャン対象ラインを少なくすることにより、アンダーサンプリング区域を取り入れてＰＲＯＰＥＬＬＥＲ又はＰＲＯＰＥＬＬＥＲ類似の手法を高速化する。例えば、分解能が２３０×２３０のスキャンが望ましいとする。その場合、ナイキスト特性（例えば、ラインの合計本数＝分解能×π／２）に従って、サンプリング対象データを収集するのに約３６０本のラインが必要になる。
【００５５】
従来の手法によれば、９個のストリップでサンプリングされない有角形状隙間部がないスキャンは、３６０本のラインを得るのに１ストリップにつき４０本のラインを必要とする。（実際にすべてのラインを収集してもよいが、パラレルイメージングの場合は、一部分のみ実際に収集し、残りの部分は数学的に合成することができるであろう。）なお、従来のＰＲＯＰＥＬＬＥＲスキャンは、典型的には、サンプリングされない楔形隙間部がない状態で行われる。それは、そのような隙間部によってストリークアーチファクトが生じるからである。
【００５６】
図５に図示する実施形態を利用して、ｋ空間が９個の有角形状アンダーサンプリング領域でカバーされる。ｋ空間の周縁部分における、サンプリングされない有角形状の割合（すなわち、隙間部の割合）は、例えば、２５％欠損又は３５％欠損として設定可能としてよい。２５％が欠損している場合、幅の狭い放射状ストリップの各対は、３０本のライン（４０本ではなく）を収集する。あるいは、ラインを３０本持つ比較的幅の広いブレードではなく、それぞれラインを１５本持つ幅の狭い２つのブレードを使用してもよい。欠損部分が３５％である場合、同様の判断により、幅の広いブレード１個あたり２６本のライン、又は幅の狭いブレード１個あたり１３本のラインが算出される。
【００５７】
したがって、実施形態においては、ＥＴＬは、ＴＲの短縮のため、あるいは１ＴＲあたりにより多くのスライスについてのデータを収集するため、又は１つのエコートレインにおけるＴ２減衰の低減を実現するために、短縮される。あるいは、エコー間スペーシング（Inter-Echo  Spacing：ＩＥＳ）を長くすることができ、それは、おそらくは、ＳＮＲの向上のため、最終的な再構成、又は、動き補正推定等の中間的な計算のためである。以上やその他のトレードオフは様々なアーチファクトを抑制するために行われる。図６に関連して説明した実施形態についても、同様の計算が行われる。
【００５８】
一部の実施形態では、使用するストリップ数を少なくすることにより、アンダーサンプリング領域を取り入れてＰＲＯＰＥＬＬＥＲ又はＰＲＯＰＥＬＬＥＲ類似の手法でのスキャンを高速化する。例えば、短縮したＥＴＬについて上述した例では、図５に従う実施形態により、少ない数のブレードを使用してｋ空間をカバーできる。
【００５９】
説明のための実施形態として、ストリップ数を９から７に減らすことが考えられる。その場合、３６０／７本、すなわち５１．５本のラインが、放射状ストリップ１個あたりに望ましいであろう（この数は、実際に収集されるライン、及び楔形状の隙間部を充填するのに必要な追加の仮想的なラインの両方を含む）。したがって、その場合も４０本のラインが収集されるのであれば（例えば、ＥＴＬを同様に維持して）、各ブレードの隙間部の幅は約１１．５ライン、すなわち１１．５本の欠損ラインに対応する。この場合は、ＴＲ及びＥＴＬは同様のまま維持される。スキャン時間は９×ＴＲから７×ＴＲに短縮され、これはスキャン時間を２２％短縮したことになる（プレスキャンは含まず）。図６に関連して説明した実施形態についても、同様の判断が行われる。
【００６０】
当然ながら、実施形態において、設計上のトレードオフを行うことは可能である。例えば、有角形状アンダーサンプリング領域は、高速化要因を高めるために拡大することができる。しかしながら、アンダーサンプリング領域の拡大により、ｋ空間の中心が少なくとも一部の放射状ストリップの端に近づく可能性があるので、ｋ空間中央領域のオーバーサンプリングが減ることにより、結局、動きエラーに対する補正能力が低下する。
【００６１】
図５に関連して説明した実施形態におけるある程度の角度上の重なり（例えば、アンダーサンプリング領域とは対立するものとして）を組み入れるために、任意選択される機能により、節約されたスキャン時間がある程度もとに戻ることになる。放射状ストリップパターン５００では、対象が放射状ストリップ５３４を収集するＴＲの間、時計回りのパターンで回転する場合、動き補正では、放射状ストリップ５３４を反時計回りに回転させてよく、新たな（例えば、図５に図示しない）アンダーサンプリング領域が、放射状ストリップ５３４及び５０４の間に（例えば、略円形のｋ空間領域５０２の上部右側に）出現する可能性がある。放射状ストリップの重なりが一定量であり、ストリップの動きエラーがその重なり量より小さい場合、そのような回転では新たなアンダーサンプリング領域は出現しない。重なり量より小さい動きや回転は、単に、放射状ストリップを、重なりが小さくなるように移動させることになるのみである。言い換えると、一部の位置にある有角形状の非サンプリング領域は、より広範囲の動きへのより強い対抗力を得ることと引き換えに、画像品質上の実質的なペナルティなしに有利に取り入れることができる。図７に示すパターン７００は、図５に示す例示的配置においてストリップ対のうち１つが一定量重なった状態を示す。このように、すべての放射状ストリップスキャン手法において中心領域で重なりが発生しているが、一部の実施形態では、さらに、収集したｋ空間領域の周縁部において、あるストリップの短手方向に沿った端部の片方において２個以上のストリップが重なることもある。
【００６２】
例えば、図７は、ストリップ７３４及び７０４が短手方向に沿った端部の片方において距離７１２だけ重なるように、（図５のストリップ５３４と比べて）幅方向が広くなったストリップ７３４を示す図である。そのような重なりが発生すると、それらのブレードのそれぞれについて、重なった領域を収集するスキャンが行われるため、余計な時間がかかることがある。しかしながら、スキャン対象被検体を距離７１２未満の量だけ回転すると、２つのストリップ７３４及び７０４の間には新たなアンダーサンプリング領域が出現することはない。なぜなら、動き補正では、単純に、ストリップ７３４（あるいは、ストリップ７３４及びストリップ７０４の少なくとも一方）を、図５のストリップ５３４及び５０４によって示されるように、最上部においてそれらのストリップが同一平面でぴったりと重なるまで（周縁で接する、臨界ナイキストサンプリング距離以上離間しないようにする）を最大限として回転させてよいからである。重なり量（例えば、距離７１２）は、動きの補正における望ましい誤差範囲として設定されてよい。このように、実施形態では、対象の動きに対してある程度の排除性（ロバスト性）を確保するために選択的に重なりを取り込むことが図られている。
【００６３】
上述した実施形態では、ｋ空間内において、有角形状非サンプリング領域のうち少なくとも１個は、その直径上反対側の位置に別の有角形状アンダーサンプリング領域がない。しかしながら、１又は複数個の有角形状の非サンプリング領域の直径上反対側の位置に別の有角形状の非サンプリング領域があってよいが、そのような互いに直径上反対側の位置にある有角形状の非サンプリング領域の対は、その対の一方のアンダーサンプリング領域が、もう一方の領域よりも大幅に大きいという実施形態も考えられる。一例として、片方の楔形のｋ空間周縁部の角度距離が、反対側の位置にある楔型の角度の大きさの３倍であってもよい。大きさの違いは、ｋ空間内の大きい方の非サンプリング領域のデータが、ｋ空間の対称性に基づく手法（例えば、エルミート対称、複素共役）でどの程度完全に埋められるかに影響し、さらに、結果としては、欠損したｋ空間データにより再構成画像内で生じるアーチファクトの量に影響する可能性がある。
【００６４】
例えば、複数のストリップは、隣り合うストリップ間の角度が不均一な状態で前記ｋ空間内に放射状に配置されてもよい。
【００６５】
また、例えば、複数のストリップのうちの少なくとも１つのストリップが、ｋ空間の中心に対して位置がずれるように配置されてもよい。例えば、複数のストリップそれぞれが、各ストリップの少なくとも短軸方向におけるｋ空間の中心からのオフセット分だけ位置がずれるように配置されてもよい。
【００６６】
図８は、１又は複数の実施形態に係る、アンダーサンプリングされたｋ空間を伴う、放射状ストリップスキャン手法（例えば、ＰＲＯＰＥＬＬＥＲ等）の高速化のための処理８００を定義するコンピュータプログラムコード構造を示す図である。動作Ｓ８０２〜Ｓ８１０は、図示の順序で実行されても別の順序で実行されてもよい。一部の実施形態では、動作Ｓ８０２〜Ｓ８１０のうち少なくとも１つは実行されなくてもよく、処理８００に少なくとも１つの他の動作を含んでもよい。処理８００は、図１に示すＭＲＩデータプロセッサ４２、ＭＲＩシステムコントローラ２２、ＭＲＩシーケンスコントローラ３０のうち少なくとも１つにより実行や制御がなされてもよい。
【００６７】
処理８００が開始されると、動作Ｓ８０２においてアンダーサンプリング領域の最適化パターンが決定される（事前に設定済みでない場合）。アンダーサンプリング領域のパターン例は、図４に図示し上述したものである。
【００６８】
動作Ｓ８０４において、決定されたアンダーサンプリング領域パターンに対応する放射状ストリップ配置が決定される（ここでも、事前に設定済みでない場合）。図５〜７は、動作Ｓ８０２で決定される通りにアンダーサンプリング領域の最適な配置が実現されるような、放射状ストリップの配置例を示す。
【００６９】
一部の実施形態では、どの程度の隙間部（例えば、隣接する放射状ストリップ間の有角形状アンダーサンプリング領域の大きさ）が生成されるかと、どの程度の重なりが組み込まれるかとの間で、柔軟なトレードオフが行われる。例えば、図５に示す、各ストリップの終端に重なりのない放射状ストリップの配置は、図７に示すように変更して、少なくとも一部のストリップの幅方向に沿った終端にいくらかの重なりが生じるようにすることもできる。ｋ空間の一方側における角度の重なりは、再構成が、例えば、面内の回転に対してどの程度の排除性（例えば、実行可能な補正の範囲）を持つようになるかに対応する。
【００７０】
動作Ｓ８０６において、ｋ空間のアンダーサンプリングデータを収集するためにＭＲＩスキャンが行われる。ＭＲＩスキャンは、ここまでのステップで決定された放射状ストリップの配置を用いて、ＰＲＯＰＥＬＬＥＲ又はＰＲＯＰＥＬＬＥＲ類似の放射状ストリップキャンに従って行われる。
【００７１】
動作Ｓ８０８において、収集されたｋ空間のアンダーサンプリングデータからＭＲ画像が再構成される。処理９００は、図９に関連して後述するが、無視できる程度のアーチファクトしか含まない最終的な再構成ＭＲ画像を再構成するのに用いられる。例えば、ｋ空間内でアンダーサンプリング領域が互いに直径上反対の位置にあることで従来の高速化ＰＲＯＰＥＬＬＥＲ手法の場合に生じるストリーキングアーチファクトは、実施形態では軽減又は除去される。動作Ｓ８１０において、再構成ＭＲ画像が出力される。再構成画像は、ディスプレイ２４等のディスプレイ上に表示、プリンタ２８等のプリンタに出力、例えばメモリ４６又はプログラム／データ記憶装置５０等の記憶装置に格納、又はさらに処理を行うためにネットワーク経由で他の装置へ送信される。
【００７２】
例えば、生成部は、ＭＲ画像を生成する前に、ストリップごとに動き推定量を設定し、設定された動き推定量を用いて、ｋ空間データにおける位相及び回転のうち少なくとも一方を補正し、補正されたｋ空間データを統合する。
【００７３】
図９は、動作Ｓ８０８においてアンダーサンプリングされたｋ空間データから画像再構成を行うための処理９００を示すフローチャートである。動作Ｓ９０２〜Ｓ９０８は、図示の順序で実行されても別の順序で実行されてもよい。一部の実施形態では、動作Ｓ９０２〜Ｓ９０８のうち少なくとも１つは実行されなくてもよく、処理９００に少なくとも１つの他の動作を含んでもよい。処理９００は、図１に示すＭＲＩデータプロセッサ４２及びＭＲＩシステムコントローラ２２のうち少なくとも一方により実行や制御がなされてもよい。
【００７４】
動作Ｓ９０２において、各放射状ストリップについて、収集されたｋ空間データの動き補正が行われる。上述したように、ｋ空間の中央領域のオーバーサンプリングは、ＰＲＯＰＥＬＬＥＲ等の放射状ストリップスキャン手法では本来備わった性質である。中心領域のオーバーサンプリングにより、各ＴＲ又は各ストリップの標準位置への動き補正が可能になる。標準位置は、例えば、すべて又は複数個の放射状ストリップの平均値に基づき、又は放射状ストリップの任意の１つに基づき選択されてよい。
【００７５】
動き補正は、スキャン領域に対する被写体の回転や平行移動によって必要とされてもよい。上述したように、そのような被写体の動きにより、アンダーサンプリング領域がどこに位置するかが変化する。拡散強調スキャンでは局所的な平行移動の動きの方が大きな問題となるが、非拡散強調スキャンでは回転性の動きの方が影響度の高い問題である。どの放射状ストリップもｋ空間中心に対して所定の位置を持つので、任意の放射状ストリップについて収集されたｋ空間データは、低分解能フーリエ変換データが移動してどの程度中心から外れたかを決定することでその所定の位置に平行移動できる。質量中心の計算のような簡単な距離計算で十分であるが、画像域の最大相関の算出が使用されてもよい。同様に、画像域における平行移動は、ｋ空間内の位置の関数である線形位相移動に対応するので、例えば、ストリップ対についてｋ空間の複素位相差への線形適合を推定することも合理的である。回転は、ｋ空間内又は画像空間内のいずれかにおいて、少なくとも１本の放射状軸に沿って中心を超えて、収集されたデータにおける特徴により検出できる。回転は、例えば、ｋ空間内における特定の複数の円形状に沿った角度移動を変数とする関数として、データの最良適合を求めることで推定してよい。ここでも、適合化は、ストリップの対からのデータ間で行うこともできるが、何らかの集約データに対して各ストリップについて行ってもよく、その場合は、おそらく反復的（すなわち、集約データの漸次改良をしながら）となるであろう。
【００７６】
又は、一部の実施形態では、ｋ空間内又は画像空間内でｋ空間全体を分析することによって相互相関を検出し、それにより動きを補正してもよい。
【００７７】
動作Ｓ９０４において、任意選択として、各ストリップが、別の形の位相回転エラーについて補正される。例えば、拡散イメージングでは、パルスシーケンスにおける拡散の大きい傾斜ローブ（gradient  lobe）の存在下で、小さな動きにより、大きな複素位相偏移がＭＲ信号に与えられることがある。この動きに起因したエラーは、硬直的な身体の平行移動及び回転の効果についての数学的モデルに従うことになるであろうが、そのような重大な位相エラーは、単一の画像画素よりもかなり小さな、画素の変位によって起こりうる。画像域での信号位相の検出及び補正は、多くの場合、画像画素位置の位置移動なしに、画像域において補正できる。この補正は、高分解能のマルチショットの拡散イメージングを行う場合に非常に重要である。そのような補正なしでは、多数の撮影について複素数データを正しく集約できない。画像域画素の複素位相回転についての補正なしでは、複数のＴＲ撮影（すなわち、複数のブレード）の集約において、量に変動のある疑似的な位相の打ち消し、及び従来のマグニチュードイメージ（magnitude  image）における多様な信号損失又は暗化が生じる。他の場合のように、位相エラーの推定は、ｋ空間域又は画像域のいずれにおいても行うことができる。同様に、位相エラーの補正は、ｋ空間域又は画像域のいずれにおいても行うことができる。動作Ｓ９０４は、図９に示すようにｋ空間に対して実行されるが、顕著に硬直的な身体の動きであるような動きには好適な選択である。このような様々な代替的補正は、現在、当業者には広く知られている。
【００７８】
動作Ｓ９０６において、任意選択として、各ストリップに対応するデータを合成する際にストリップ内における重み付け又はストリップ間の重み付けが反映されるように、放射状ストリップへの重み付けが行われる。重みの変更は、ｋ空間合成の前かつ密度補償の前に、（ａ）相関関係の質の低さ、及び、面内の動きのペナルティ等のうちの少なくとも１つ、又は、（ｂ）ＴＥ（Echo  Time）が望ましいＴＥからどの程度離れているか等、すなわち、コントラストがどの程度準最適かによって決まるエコーの質の意味を与えるように行われてよい。
【００７９】
例えば、生成部は、ＭＲ画像を生成する前に、オーバーサンプリングされた領域において、各ストリップに対して、ストリップ間のデータの相関関係の質、又は、面内の動きのペナルティに関する重みを付与し、付与された重みを用いてｋ空間データを統合する。
【００８０】
また、例えば、生成部は、ＭＲ画像を生成する前に、オーバーサンプリングされた領域において、各ストリップに含まれる各データ収集ラインに対して、データ収集ライン間の相対的なエコーの質に関する重みを付与し、付与された重みを用いてｋ空間データを統合する。このとき、例えば、各データ収集ラインに付与される重みは、オーバーサンプリングされた領域における２本以上のライン間のエコーのコントラストの変化に従って定められる。
【００８１】
図５及び図６等の中央領域においてオーバーサンプリングされたデータを検討する。効率的で広く理解されている高速フーリエ変換（Fast  Fourier  Transform：ＦＦＴ）及び離散フーリエ変換（Discrete  Fourier  Transform：ＤＦＴ）のアルゴリズムを用いて再構成するには、規則的に間をおいた矩形のデカルト格子上の各格子交点において使用される１つのサンプルが必要になるであろう。しかし、オーバーサンプリング領域では、理想的な各サンプル位置のごく近傍で、複数のデータ収集位置が利用可能である。１つのｋ空間位置のサンプルを推定するために近傍の収集データを集約する場合、その集約はどう行われるべきであろうか。単純にデータの平均をとるのは、ランダムな熱雑音の信号対ノイズに関しては適切な対処法である。しかし、ＰＲＯＰＥＬＬＥＲイメージングでは、もっとよい方法がある可能性がある。ＰＲＯＰＥＬＬＥＲでは、多くの場合、動きによって変更あるいは「破損」されたデータを扱う。動き破損の少ないと思われるデータを優先的に使用することにより改良を図るのは合理的である。簡単な対処法は、理想的なデカルト格子上のサンプル位置の近傍における利用可能な各収集データ点に重みを関連付けることである。同様に、単一のＦＳＥエコートレイン内の異なる複数のエコーが、理想に近いか遠いかを判定してもよい。より理想的なエコーからのデータに、より高い重みが付けられる。収集されたｋ空間領域全体の複数の部分において、データは、データの１つのストリップからのみ供給されている。冗長データのない領域では、重み付けは非現実的となる可能性がある。あるデータが使用するに十分近い場合は、想定されるデータ品質においてグラデーション（gradation）があったとしても、単純にそれを使用すべきである。
【００８２】
再構成に格子を設定するのは、非デカルト軌道を用いて収集されたデータを再構成するための標準的なＭＲ方法であり、そこでは重み付け関数を活用する。理想的なｋ空間格子点の近傍内で収集された複数のサンプルは、その格子点に足し込まれる際、各サンプルは収集されたサンプル位置と望ましい理想的な（ｋ空間内の）格子点位置間距離の関数によって機械的に乗算、すなわち重み付けされる。この重み付け関数は、ｓｉｎｃ関数又はｓｉｎｃ＊ｗｉｎｄｏｗ関数の関数形式をとり、このｗｉｎｄｏｗ関数としては、信号処理におけるカイザーベッセルウィンドウ（Ｋａｉｓｅｒ−Ｂｅｓｓｅｌ  ｗｉｎｄｏｗ）のようなものを使用できる。ＰＲＯＰＥＬＬＥＲ再構成では、この格子化重み付け値をさらに他の重み付け情報と組み合わせることが妥当と思われる。特に、ある１画素の格子化重み付け値は、撮影間のデータ品質の一貫性に関連する別の重み付け関数、又はより理想的でないエコー時間からのデータにペナルティを与える別の重み付け関数、あるいは追加的なこれらの関数の両方で乗算することができる。
【００８３】
簡単であるが実用的なデータ品質の指標として、ｋ空間の中央位置における、ｋ空間データの値、すなわち、ｋ空間データの値の絶対的な大きさを使用できる。この指標は、面内平行移動の動きについてｋ空間を補正した後に評価できるであろう。そのようなｋ空間中央部の値の変動は、拡散強調スキャンの間の、影響度の高い非硬直的な身体の動きによる信号キャンセルの指標として使用できる。同様に、そのような変動は、イメージング面の垂直方向に沿った平行移動、又はイメージング面の回転による傾きからの、イメージングスライスの変化と関連している可能性がある。なお、当然ながら、そのようなイメージングスライスの変化は、ＭＲＩのハードウェアシステムの機械的な座標と相対的に関連するものではなく、むしろ、意図されたスライス及び関心被検体の座標（おそらくは人間である被写体の脳内のスライス）と相対的に関連する。また、例えば、「スライスの面内回転についての補正」について言及するのはよくあることではあるが、この数学的補正は、「関心解剖部位の面内回転」（解剖部位が反対方向に回転したこと）から物質界で発生する可能性がある。
【００８４】
そのような重み付け関数の例は、この場合、ストリップ全体の何らかのばらつきに関する指標であるが、次のように表せる。
【００８５】
diff（strip）＝ABS（k0（strip）−k0（aggregate））            ・・・（１）
relativediff（strip）＝diff（strip）／ABS（k0（aggregate））    ・・・（２）
weight（strip）＝SQR（MAX（relativediff（strip），0））        ・・・（３）
【００８６】
このような重みは、０から１の範囲の値を返す。上記の等式（１）〜（３）は、あるストリップの差違値を、そのストリップのｋ空間中心領域値と全ストリップのｋ空間中心領域集約値との差の絶対値として決定し、全ストリップのｋ空間中心領域集約値の絶対値に対する、そのストリップの差違値の比に基づいて、そのストリップの相対的差違値を決定することを示す。そうすると、そのストリップの重みは、そのストリップの相対的差違値の二乗として算出できる。
【００８７】
放射状ストリップにおける各データ収集ラインは各エコーに対応する。一般的には、任意の特定の放射状ストリップについて、幅方向の中央にあるデータライン（例えば、長軸）は、画像において比較的大きい支配的な特徴についての特色をよく示しているＴ２コントラストを呈示するという理由で最も正確なデータを含んでいると考えられる。そのため、ラインは、中央のラインが最も大きな重みを持ち、ラインが幅方向中央のデータラインから離れるにつれ各重みが次第に小さくなるように重み付けされてよい。これにより、理想から遠い緩和（relaxation）コントラストにペナルティが課される。そのような重み付け関数の簡単な例が、以下のように表せる。
【００８８】
Wecho＝ABS（TE−TEcenter）／（a＊TEmax）      ・・・（４）
【００８９】
上記式において、「ａ」は、１．０より大きい値になるように選択される自由設計パラメータである。等式（４）は、エコーの重みが最大ＴＥに対する、そのエコーのＴＥとｋ空間中心領域ＴＥとの絶対差の比として決められることを示す。
【００９０】
オーバーサンプリングされた領域においてデータライン又はストリップに重みを付与する場合に、各ｋ空間位置の値は、その位置に重なっているストリップそれぞれからの対応する値の重み付き平均値によって決められてもよい。
【００９１】
一部の実施形態によると、重み付けをしたとしても、データを何らかの標準に従わせることは有益である場合もある。例えば、放射状ストリップは、特定のＴＲにおいて、その時の意図したスライスとは異なるスライスのデータを収集してしまう場合（イメージングスライスに垂直な方向での大きな患者の動きによる）、動き補正が行われた後でも、値がまだ標準的なものとはかなり違っている場合もある。そのような相違が検出された場合には重みを調整してもよいし、そのようなデータは無視することもできる。データの無視は重み付けのステップにおいて、重み付け係数を０．０にすることで実現できる。収集されたデータが標準に十分近いかどうかを判断するデータ閾値は、いくつかある手法のいずれかを用いて決定でできる。一例として、あるストリップの最大ｋ空間データ幅が、最初の平均化された集約後のｋ空間の最大データ空間幅の０．５倍未満である場合に、その特定のストリップの重み付けがゼロ値になってもよい。
【００９２】
動作Ｓ９０８において、動作Ｓ９０２〜Ｓ９０６の少なくとも１つによって処理された後の収集データを使用して部分フーリエ再構成が行われる。いくつかある部分フーリエ法のいずれか（例えば、ホモダイン法、Margosian法、又はＰＯＣＳ法）によって実行される。ホモダイン法は、通常は画像空間において、画像のローパスフィルタを利用するものである。あるいは、集約後のｋ空間において窓関数を利用する、Margosian法も適用できる。あるいは、部分フーリエ再構成用のＰＯＣＳ法を、すべての放射状ストリップについての統合されたデータセットに対して行ってもよい（各放射状ストリップが別々に動きについて補正された後）。
【００９３】
ＰＯＣＳ法の利用は、アンダーサンプリング領域が充填される前の収集データから始めてもよい。動き補正をはじめとする補正は、各放射状ストリップについて既に行われていてもよい。ＰＯＣＳ法では、まず、収集データを、アンダーサンプリング領域のゼロ値化によって画像空間（例えば、実空間）に変換する。この段階で得られる画像には、例えば、ｋ空間の未充填領域が原因で、不具合がある。ローパス位相フィルタのような、複素数値化された（Ｒ＋ｊＩ）データのフィルタリングが行われてもよい（画素信号の非負性を強制する等の、追加的処理も可能である）。上記のフィルタリングされていない画像から、この際立って位相が整った信号Ｒを引く減算が行われる。これらの減算した画像は、実際のところ、アンダーサンプリングに由来する画像アーチファクト部分を多く含むと同時に、少ないが望ましい正しい信号部分を含む。
【００９４】
減算した画像からの生成物は、ｋ空間に変換され、最初に収集されたｋ空間データと統合されて「完全に充填された」ｋ空間（その後、これを元に好適なフーリエ変換再構成によって再構成画像が得られる）に逆変換される。
【００９５】
ある条件のもとでは、欠損データのあるｋ空間に対してフーリエ再構成を行うというＰＯＣＳアルゴリズムの概念の利用は、旧来の離散的なＧｅｒｃｈｂｅｒｇ−Ｐａｐｏｕｌｉｓの信号処理方法として簡略になる。これらの部分フーリエアルゴリズムは、ＭＲＩ技術文献で広く知られている。
【００９６】
なお、少なくとも一部の実施形態によると、部分フーリエ再構成は、複数のストリップからの収集データを統合したデータに対して行われる。それとは対照的に、一部の従来の手法では、同等のフーリエ再構成を、データの統合前に各ストリップについて行う。部分フーリエ再構成を行う利点には、（１）集約後データの信号対ノイズが利用可能になる、（２）位相フィルタリング改善のために、集約後データ内から利用可能なデータの範囲が広い、（３）オーバーサンプリング部分での選択的重み付けにより集約後データにおいて高品質化されたデータになっている、ということが含まれうる。これら３つの利点を持つ方法は、画像品質を高める場合がある。
【００９７】
例示的実施形態により、従来のＰＲＯＰＥＬＬＥＲ及びＰＲＯＰＥＬＬＥＲ類似の手法に優る利点がいくつか提供される。例えば、近い又は同じ画像品質であれば、スキャン速度が速い。さらに、結果的なサンプリングパターン及びストリップの重なりパターンは、面内の回転する動きに起因する、アンダーサンプリングによるアーチファクトに対して良好な排除性を持つ。また、例示的実施形態によれば、統合データに対する部分フーリエ再構成により、各ブレードに対する部分フーリエ再構成を行う手法に比べて、十分にサンプリングされたデータを持つ中央領域が大きくなる。十分に収集されたデータを持つ中央領域が大きいほど、より高い空間周波数成分が正しい信号位相により呈示されるような、画像領域における再構成を実現することができる。
【００９８】
特定の実施形態を説明したが、これらの実施形態は、単なる例示として提示したものであり、発明の範囲を限定するものではない。実際、本明細書で説明した新規の実施形態は、様々な別の態様で具体化が可能である。さらには、本発明の要旨を逸脱することなく、本明細書で説明した実施形態の態様において様々な省略、置換、及び変更が可能である。添付の特許請求の範囲及びそれらの均等物は、本発明の範囲や要旨に収まると考えられる態様又は変更を包含するものである。
【符号の説明】
【００９９】
１０  ガントリ
２０  関連システム構成要素
３０  ＭＲＩシーケンスコントローラ
４２  ＭＲＩデータプロセッサ
【図１】
0006691797.tif 000002
【図２】
0006691797.tif 000003
【図３】
0006691797.tif 000004
【図４】
0006691797.tif 000005
【図５】
0006691797.tif 000006
【図６】
0006691797.tif 000007
【図７】
0006691797.tif 000008
【図８】
0006691797.tif 000009
【図９】
0006691797.tif 000010