JP2021029777A - 磁気共鳴データ収集方法及び磁気共鳴イメージング装置 - Google Patents

Abstract

【課題】データ収集時の被検体の体動をより正確に把握することが可能なデータを得ること。【解決手段】実施形態に係る磁気共鳴データ収集方法は、１回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集する。前記１ショットで収集される磁気共鳴データは、ｋ空間の中心領域に対応するデータと、前記ｋ空間の辺縁領域に対応するデータとを含む。前記１ショットの収集において、前記中心領域に対応するデータはフルサンプリングで収集され、前記辺縁領域に対応するデータは間引きサンプリングで収集される。【選択図】図２

Description

本発明の実施形態は、磁気共鳴データ収集方法及び磁気共鳴イメージング装置に関する。

従来、磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）装置において、撮像を高速化させるための技術として、ｋ空間のデータを間引いて収集し、収集されたデータから圧縮センシング等によって画像を生成する撮像法が知られている。

ここで、圧縮センシングとは、信号のスパース（疎）性を利用して、ゼロ成分を削減したデータに含まれる非ゼロ成分から信号を復元する手法であり、少ないデータから画像を生成するための方法として、近年、ＭＲＩ装置の分野でも注目されている方法である。

また、一般的に、ＭＲＩ装置では、データ収集時に被検体の体動が発生すると、その影響で画像にアーチファクトが生じることが知られているが、上述した高速化のための撮像法では、体動の影響が考慮されていないことも多い。このような体動の影響を適切に補正するためには、データ収集時の被検体の体動をより正確に把握することが求められる。

特表２０１６−５３９７２２号公報 特開２０１５−０５１２６４号公報 特開２００９−２６８９０１号公報

本発明が解決しようとする課題は、データ収集時の被検体の体動をより正確に把握することが可能なデータを得ることである。

実施形態に係る磁気共鳴データ収集方法は、１回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集する。前記１ショットで収集される磁気共鳴データは、ｋ空間の中心領域に対応するデータと、前記ｋ空間の辺縁領域に対応するデータとを含む。前記１ショットの収集において、前記中心領域に対応するデータはフルサンプリングで収集され、前記辺縁領域に対応するデータは間引きサンプリングで収集される。

図１は、第１の実施形態に係るＭＲＩ装置の構成例を示す図である。 図２は、第１の実施形態に係る収集機能によって行われるデータ収集の一例を示す図である。 図３は、第１の実施形態に係る比較例を示す図である。 図４は、第１の実施形態に係る補正機能によって行われる体動補正の一例を示す図である。 図５は、第１の実施形態に係るＭＲＩ装置によって行われる処理の処理手順を示すフローチャートである。 図６は、第１の実施形態の変形例に係る収集機能によって行われる磁気共鳴データの収集の一例を示す図である。 図７は、第２の実施形態に係る収集機能によって行われるデータ収集の一例を示す図である。 図８は、第２の実施形態に係る収集機能によって行われるデータ収集の一例を示す図である。 図９は、第２の実施形態に係る比較例を示す図である。

以下、図面を参照しながら、本願に係る磁気共鳴データ収集方法及び磁気共鳴イメージング装置の実施形態について詳細に説明する。

（第１の実施形態）
図１は、第１の実施形態に係るＭＲＩ装置の構成例を示す図である。

例えば、図１に示すように、ＭＲＩ装置１００は、静磁場磁石１、傾斜磁場コイル２、傾斜磁場電源３、全身用ＲＦコイル４、局所用ＲＦコイル５、送信回路６、受信回路７、ＲＦ（Radio Frequency）シールド８、架台９、寝台１０、インタフェース１１、ディスプレイ１２、記憶回路１３、及び処理回路１４〜１６を備える。

静磁場磁石１は、被検体Ｓが配置される撮像空間に静磁場を発生させる。具体的には、静磁場磁石１は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、その内周側に形成された撮像空間に静磁場を発生させる。例えば、静磁場磁石１は、超伝導磁石や永久磁石等である。ここでいう超伝導磁石は、例えば、液体ヘリウム等の冷却剤が充填された容器と、当該容器に浸漬された超伝導コイルとから構成される。

傾斜磁場コイル２は、静磁場磁石１の内側に配置されており、被検体Ｓが配置される撮像空間に傾斜磁場を発生させる。具体的には、傾斜磁場コイル２は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、互いに直交するＸ軸、Ｙ軸及びＺ軸それぞれに対応するＸコイル、Ｙコイル及びＺコイルを有している。Ｘコイル、Ｙコイル及びＺコイルは、傾斜磁場電源３から供給される電流に基づいて、各軸方向に沿って線形に変化する傾斜磁場を撮像空間に発生させる。ここで、Ｚ軸は、静磁場磁石１によって発生する静磁場の磁束に沿うように設定される。また、Ｘ軸は、Ｚ軸に直交する水平方向に沿うように設定され、Ｙ軸は、Ｚ軸に直交する鉛直方向に沿うように設定される。これにより、Ｘ軸、Ｙ軸及びＺ軸は、ＭＲＩ装置１００に固有の装置座標系を構成する。

傾斜磁場電源３は、傾斜磁場コイル２に電流を供給することで、撮像空間に傾斜磁場を発生させる。具体的には、傾斜磁場電源３は、傾斜磁場コイル２のＸコイル、Ｙコイル及びＺコイルに個別に電流を供給することで、互いに直交するリードアウト方向、位相エンコード方向及びスライス方向それぞれに沿って線形に変化する傾斜磁場を撮像空間に発生させる。なお、以下では、リードアウト方向に沿った傾斜磁場をリードアウト傾斜磁場と呼び、位相エンコード方向に沿った傾斜磁場を位相エンコード傾斜磁場と呼び、スライス方向に沿った傾斜磁場をスライス傾斜磁場と呼ぶ。

ここで、リードアウト傾斜磁場、位相エンコード傾斜磁場及びスライス傾斜磁場は、それぞれ静磁場磁石１によって発生する静磁場に重畳されることで、被検体Ｓから発生する磁気共鳴信号に空間的な位置情報を付与する。具体的には、リードアウト傾斜磁場は、リードアウト方向の位置に応じて磁気共鳴信号の周波数を変化させることで、リードアウト方向に沿った位置情報を磁気共鳴信号に付与する。また、位相エンコード傾斜磁場は、位相エンコード方向に沿って磁気共鳴信号の位相を変化させることで、位相エンコード方向に沿った位置情報を磁気共鳴信号に付与する。また、スライス傾斜磁場は、スライス方向に沿った位置情報を磁気共鳴信号に付与する。例えば、スライス傾斜磁場は、撮像領域がスライス領域（２Ｄ撮像）の場合には、スライス領域の方向、厚さ及び枚数を決めるために用いられ、撮像領域がボリューム領域（３Ｄ撮像）の場合には、スライス方向の位置に応じて磁気共鳴信号の位相を変化させるために用いられる。これにより、リードアウト方向に沿った軸、位相エンコード方向に沿った軸、及びスライス方向に沿った軸は、撮像の対象となるスライス領域又はボリューム領域を規定するための論理座標系を構成する。

全身用ＲＦコイル４は、傾斜磁場コイル２の内周側に配置されており、撮像空間に配置された被検体ＳにＲＦ磁場を印加し、当該ＲＦ磁場の影響によって被検体Ｓから発生する磁気共鳴信号を受信する。具体的には、全身用ＲＦコイル４は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、送信回路６から供給されるＲＦパルス信号に基づいて、その内周側に位置する撮像空間に配置された被検体ＳにＲＦ磁場を印加する。また、全身用ＲＦコイル４は、ＲＦ磁場の影響によって被検体Ｓから発生する磁気共鳴信号を受信し、受信した磁気共鳴信号を受信回路７へ出力する。

局所用ＲＦコイル５は、被検体Ｓから発生した磁気共鳴信号を受信する。具体的には、局所用ＲＦコイル５は、被検体Ｓの部位ごとに用意されており、被検体Ｓの撮像が行われる際に、撮像対象の部位の表面近傍に配置される。そして、局所用ＲＦコイル５は、全身用ＲＦコイル４によって印加されたＲＦ磁場の影響によって被検体Ｓから発生した磁気共鳴信号を受信し、受信した磁気共鳴信号を受信回路７へ出力する。なお、局所用ＲＦコイル５は、被検体ＳにＲＦ磁場を印加する機能をさらに有していてもよい。その場合には、局所用ＲＦコイル５は、送信回路６に接続され、送信回路６から供給されるＲＦパルス信号に基づいて、被検体ＳにＲＦ磁場を印加する。例えば、局所用ＲＦコイル５は、サーフェスコイルや、複数のサーフェスコイルをコイルエレメントとして組み合わせて構成されたフェーズドアレイコイルである。

送信回路６は、静磁場中に置かれた対象原子核に固有のラーモア周波数に対応するＲＦパルス信号を全身用ＲＦコイル４に出力する。具体的には、送信回路６は、パルス発生器、ＲＦ発生器、変調器、及び増幅器を有する。パルス発生器は、ＲＦパルス信号の波形を生成する。ＲＦ発生器は、共鳴周波数のＲＦ信号を発生する。変調器は、ＲＦ発生器によって発生したＲＦ信号の振幅をパルス発生器によって発生した波形で変調することで、ＲＦパルス信号を生成する。増幅器は、変調器によって生成されたＲＦパルス信号を増幅して全身用ＲＦコイル４に出力する。

受信回路７は、全身用ＲＦコイル４又は局所用ＲＦコイル５から出力される磁気共鳴信号に基づいて磁気共鳴データを生成し、生成した磁気共鳴データを処理回路１５に出力する。例えば、受信回路７は、選択器、前段増幅器、位相検波器、及び、Ａ／Ｄ（Analog/Digital）変換器を備える。選択器は、全身用ＲＦコイル４又は局所用ＲＦコイル５から出力される磁気共鳴信号を選択的に入力する。前段増幅器は、選択器から出力される磁気共鳴信号を電力増幅する。位相検波器は、前段増幅器から出力される磁気共鳴信号の位相を検波する。Ａ／Ｄ変換器は、位相検波器から出力されるアナログ信号をデジタル信号に変換することで磁気共鳴データを生成し、生成した磁気共鳴データを処理回路１５に出力する。なお、ここで、受信回路７が行うものとして説明した各処理は、必ずしも全ての処理が受信回路７で行われる必要はなく、全身用ＲＦコイル４又は局所用ＲＦコイル５で一部の処理（例えば、Ａ／Ｄ変換器による処理等）が行われてもよい。

ＲＦシールド８は、傾斜磁場コイル２と全身用ＲＦコイル４との間に配置されており、全身用ＲＦコイル４によって発生するＲＦ磁場から傾斜磁場コイル２を遮蔽する。具体的には、ＲＦシールド８は、中空の略円筒状（円筒の中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、傾斜磁場コイル２の内周側の空間に、全身用ＲＦコイル４の外周面を覆うように配置されている。

架台９は、略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成された中空のボア９ａを有し、静磁場磁石１、傾斜磁場コイル２、全身用ＲＦコイル４、及びＲＦシールド８を収容している。具体的には、架台９は、ボア９ａの外周側に全身用ＲＦコイル４を配置し、全身用ＲＦコイル４の外周側にＲＦシールド８を配置し、ＲＦシールド８の外周側に傾斜磁場コイル２を配置し、傾斜磁場コイル２の外周側に静磁場磁石１を配置した状態で、それぞれを収容している。ここで、架台９が有するボア９ａ内の空間が、撮像時に被検体Ｓが配置される撮像空間となる。

寝台１０は、被検体Ｓが載置される天板１０ａを備え、被検体Ｓの撮像が行われる際に、被検体Ｓが載置された天板１０ａを撮像空間に移動する。例えば、寝台１０は、天板１０ａの長手方向が静磁場磁石１の中心軸と平行になるように設置されている。

なお、ここでは、ＭＲＩ装置１００が、静磁場磁石１、傾斜磁場コイル２及び全身用ＲＦコイル４それぞれが略円筒状に形成された、いわゆるトンネル型の構造を有する場合の例を説明するが、実施形態はこれに限られない。例えば、ＭＲＩ装置１００は、被検体Ｓが配置される撮像空間を挟んで対向するように一対の静磁場磁石、一対の傾斜磁場コイル及び一対のＲＦコイルを配置した、いわゆるオープン型の構造を有していてもよい。このようなオープン型の構造では、一対の静磁場磁石、一対の傾斜磁場コイル及び一対のＲＦコイルによって挟まれた空間が、トンネル型の構造におけるボアに相当する。

インタフェース１１は、操作者から各種指示及び各種情報の入力操作を受け付ける。具体的には、インタフェース１１は、処理回路１７に接続されており、操作者から受け取った入力操作を電気信号へ変換して処理回路１７に出力する。例えば、インタフェース１１は、撮像条件や関心領域（Region Of Interest：ＲＯＩ）の設定等を行うためのトラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。なお、本明細書において、インタフェース１１は、マウス、キーボード等の物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路もインタフェース１１の例に含まれる。

ディスプレイ１２は、各種情報及び各種画像を表示する。具体的には、ディスプレイ１２は、処理回路１７に接続されており、処理回路１７から送られる各種情報及び各種画像のデータを表示用の電気信号に変換して出力する。例えば、ディスプレイ１２は、液晶モニタやＣＲＴモニタ、タッチパネル等によって実現される。

記憶回路１３は、各種データを記憶する。具体的には、記憶回路１３は、磁気共鳴データや画像データを記憶する。例えば、記憶回路１３は、ＲＡＭ（Random Access Memory）、フラッシュメモリ等の半導体メモリ素子やハードディスク、光ディスク等によって実現される。

処理回路１４は、寝台制御機能１４ａを有する。寝台制御機能１４ａは、制御用の電気信号を寝台１０へ出力することで、寝台１０の動作を制御する。例えば、寝台制御機能１４ａは、インタフェース１１を介して、天板１０ａを長手方向、上下方向又は左右方向へ移動させる指示を操作者から受け付け、受け付けた指示に従って天板１０ａを移動するように、寝台１０が有する天板１０ａの移動機構を動作させる。

処理回路１５は、各種のパルスシーケンスを実行することで、被検体Ｓの磁気共鳴データを収集する。具体的には、処理回路１５は、処理回路１７から出力されるシーケンス実行データに従って傾斜磁場電源３、送信回路６及び受信回路７を駆動することで、各種のパルスシーケンスを実行する。ここで、シーケンス実行データは、パルスシーケンスを表すデータであり、傾斜磁場電源３が傾斜磁場コイル２に電流を供給するタイミング及び供給する電流の強さ、送信回路６が全身用ＲＦコイル４にＲＦパルス信号を供給するタイミング及び供給するＲＦパルスの強さ、受信回路７が磁気共鳴信号をサンプリングするタイミング等を規定した情報である。そして、処理回路１５は、パルスシーケンスを実行した結果として受信回路７から出力される磁気共鳴データを受信し、記憶回路１３に記憶させる。このとき、記憶回路１３に記憶される磁気共鳴データは、前述したリードアウト傾斜磁場、位相エンコード傾斜磁場、及びスライス傾斜磁場によってリードアウト方向、フェーズアウト方向及びスライス方向の各方向に沿った位置情報が付与されることで、２次元又は３次元のｋ空間を表すデータとして記憶される。

処理回路１６は、処理回路１５によって収集された磁気共鳴データに基づいて、各種の画像を生成する。具体的には、処理回路１６は、処理回路１５によって収集された磁気共鳴データを記憶回路１３から読み出し、読み出した磁気共鳴データにフーリエ変換等の再構成処理を施すことで、２次元又は３次元の画像を生成する。そして、処理回路１６は、生成した画像を記憶回路１３に記憶させる。

処理回路１７は、ＭＲＩ装置１００が有する各構成要素を制御することで、ＭＲＩ装置１００の全体制御を行う。具体的には、処理回路１７は、操作者から各種指示及び各種情報の入力操作を受け付けるためのＧＵＩ（Graphical User Interface）をディスプレイ１２に表示し、インタフェース１１を介して受け付けられた入力操作に応じて、ＭＲＩ装置１００が有する各構成要素を制御する。例えば、処理回路１７は、操作者によって入力された撮像条件に基づいてシーケンス実行データを生成し、生成したシーケンス実行データを処理回路１５に出力することで、磁気共鳴データを収集させる。また、例えば、処理回路１７は、処理回路１６を制御することで、処理回路１５によって収集された磁気共鳴データに基づいて画像を生成させる。また、例えば、処理回路１７は、操作者からの要求に応じて、記憶回路１３に記憶された画像を読み出し、読み出した画像をディスプレイ１２に表示させる。

以上、本実施形態に係るＭＲＩ装置１００の構成例について説明した。このような構成のもと、本実施形態に係るＭＲＩ装置１００は、ｋ空間のデータを間引いて収集し、収集されたデータから圧縮センシングによって画像を生成する撮像法を行う機能を有する。

そして、本実施形態に係るＭＲＩ装置１００は、このような撮像法において、被検体の体動の影響を適切に補正するために、データ収集時の被検体の体動をより正確に把握することが可能なデータを得ることができるように構成されている。

具体的には、本実施形態では、処理回路１５が、収集機能１５ａを有する。また、処理回路１７が、補正機能１７ａを有する。また、処理回路１６が、第１の生成機能１６ａと、第２の生成機能１６ｂとを有する。ここで、収集機能１５ａは、収集部の一例である。また、補正機能１７ａは、補正部の一例である。また、第１の生成機能１６ａは、第１の生成部の一例である。また、第２の生成機能１６ｂは、第２の生成部の一例である。

収集機能１５ａは、１回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集する。ここで、１ショットで収集される磁気共鳴データは、ｋ空間の中心領域に対応するデータと、ｋ空間の辺縁領域に対応するデータとを含む。そして、１ショットの収集において、中心領域に対応するデータはフルサンプリングで収集され、辺縁領域に対応するデータは間引きサンプリングで収集される。

例えば、収集機能１５ａは、ＦＳＥ（Fast Spin Echo）系のパルスシーケンス、又は、ＦＦＥ（Fast Field Echo）系のパルスシーケンスを用いて、磁気共鳴データを収集する。

なお、以下では、このように、１回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集することを「マルチショット収集」と呼ぶ。

図２は、第１の実施形態に係る収集機能１５ａによって行われるデータ収集の一例を示す図である。また、図３は、第１の実施形態に係る比較例を示す図である。

なお、ここでは、ＭＲＩ装置１００によって２Ｄ撮像が行われる場合の例を説明する。

例えば、図２に示すように、収集機能１５ａは、２Ｄ撮像が行われる場合には、リードアウト方向の位置を示すＫ_ro軸と、位相エンコード方向の位置を示すＫ_pe軸とで表される２次元のｋ空間に対応するデータを収集する。

ここで、例えば、収集機能１５ａは、リードアウト方向及び位相エンコード方向それぞれに沿って、撮像条件として設定されたマトリクスサイズに応じた数のサンプリング点を設定する。

そして、収集機能１５ａは、マルチショット収集によって、位相エンコード方向の位置を変えながら、リードアウト方向に沿ったラインごとにデータを収集する。具体的には、収集機能１５ａは、マルチショット収集によって、位相エンコード傾斜磁場の強度を変えながら複数のエコー信号を収集し、１つのエコー信号を収集するごとに、収集したエコー信号に基づいて、リードアウト方向に沿った１ライン分のデータをサンプリングする。

例えば、図２に示すように、収集機能１５ａは、Ｎ回のショット（Ｓｈｏｔ＃１、＃２、＃３・・・＃Ｎ）によって、位相エンコード方向の複数の位置について、リードアウト方向に沿ったラインごとにデータをサンプリングする。ここで、図２の右側に示す複数の黒い点は、それぞれ、各ショットで収集される複数のエコー信号によってデータがサンプリングされるラインの位置を示している。

この場合に、例えば、収集機能１５ａは、ｋ空間を位相エンコード方向に分割することで、ｋ空間の中心を含む中心領域Ｒ１と、中心領域Ｒ１の外側の２つの辺縁領域Ｒ２とを設定する。例えば、収集機能１５ａは、撮像条件として設定されたマトリクスサイズが２５６×２５６であった場合に、リードアウト方向に沿って同一直線上に並ぶ複数のサンプリング点を１ラインとして、ｋ空間の中心を通るＫ_ro軸の付近に位置する複数ライン分の領域であって、位相エンコード方向に連続する複数ライン分の領域を中心領域Ｒ１として設定する。

そして、収集機能１５ａは、マルチショット収集の各ショットにおいて、中心領域Ｒ１に対応するデータをフルサンプリングで収集し、２つの辺縁領域Ｒ２に対応するデータを間引きサンプリングで収集する。ここで、例えば、間引きサンプリングは、ランダムにデータを間引くサンプリングである。なお、ここでいう「ランダム」には、完全なランダムだけではなく、疑似的なランダムも含まれる。

例えば、図２に示すように、収集機能１５ａは、辺縁領域Ｒ２については、各ショットにおいて、領域内に含まれるリードアウト方向に沿った複数のラインを位相エンコード方向に間引いてデータを収集する。このとき、収集機能１５ａは、データを収集するラインを、ショット間で重複しないように複数のショットに所定数ずつ割り当てて、データを収集する。

一方、中心領域Ｒ１については、収集機能１５ａは、各ショットにおいて、領域内に含まれる全てのラインのデータを収集する。

このように、マルチショット収集によって中心領域Ｒ１に対応するデータをフルサンプリングで収集する場合には、例えば、図３に示すように、中心領域Ｒ１に含まれるリードアウト方向に沿った複数のラインを、ショット間で重複しないように複数のショットに所定数ずつ割り当てて、データを収集する方法も考えられる。この場合には、複数のショットが行われることで、フルサンプリングされた中心領域Ｒ１のデータが１つ得られることになる。

これに対し、本実施形態では、収集機能５ａは、例えば、図２に示すように、全てのショットにおいて、中心領域Ｒ１に含まれる全てのラインのデータを収集する。これにより、本実施形態では、ショットごとに、フルサンプリングされた中心領域Ｒ１のデータが得られることになる。

このように、収集機能５ａが、各ショットにおいて、フルサンプリングされたｋ空間の中心領域Ｒ１のデータを収集することによって、ショットごとに、画像の全体像を表す低周波成分のデータを得ることができる。このデータによれば、データ収集時の被検体の体動をより正確に把握することが可能になる。

図１に戻って、補正機能１７ａは、収集機能１５ａによってフルサンプリングで収集されたｋ空間の中心領域に対応するデータから生成されたリファレンスデータを用いて、磁気共鳴データにおける体動の影響を補正する。具体的には、補正機能１７ａは、リファレンスデータを複数のショット間で比較することで、磁気共鳴データにおける体動の影響を補正する。ここで、体動の影響の補正は、体動の影響が大きいリファレンスデータに対応するショットのデータを破棄することを含む。

図４は、第１の実施形態に係る補正機能１７ａによって行われる体動補正の一例を示す図である。

例えば、図４に示すように、補正機能１７ａは、図２に示したようにＮ回のショットによってデータが収集されていた場合に、各ショットで収集された中心領域Ｒ１に対応するデータに基づいて、ショットごとに、リファレンスデータとして低周波画像を生成する。このとき、補正機能１７ａは、ショットごとに、中心領域Ｒ１に対応するデータにフーリエ変換等の再構成処理を施すことで、低周波画像を生成する。

その後、補正機能１７ａは、１ショット目の低周波画像を基準とし、２〜Ｎショット目の各低周波画像について、１ショット目の低周波画像との間の相関度を計算する。そして、補正機能１７ａは、計算した相関度に基づいて、収集機能１５ａによって収集された磁気共鳴データにおける体動の影響を補正する。

例えば、補正機能１７ａは、ｋショット目の低周波画像について、計算した相関度に基づいて、当該画像が１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものであるか否かを判定する。

ここで、補正機能１７ａは、ｋショット目の低周波画像が１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものであると判定した場合には、両画像の間の移動量を計算する。そして、補正機能１７ａは、計算した移動量に基づいて、両画像の位置が合うように、ｋショット目の磁気共鳴データを補正する。このとき、補正機能１７ａは、中心領域Ｒ１に対応するデータ（低周波データ）だけでなく、辺縁領域Ｒ２に対応するデータ（高周波データ）も含めて、ｋショット目の磁気共鳴データの全体を補正する。

また、補正機能１７ａは、ｋショット目の低周波画像が１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものではないと判定した場合には、当該低周波画像が１ショット目の低周波画像に対して局所的な動きによる変形を行って得られるものであるか否かを判定する。

そして、補正機能１７ａは、ｋショット目の低周波画像が１ショット目の低周波画像に対して局所的な動きによる変形を行って得られるものであると判定した場合には、体動の影響が大きいとみなし、ｋショット目の磁気共鳴データを破棄する。

一方、ｋショット目の低周波画像が１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものではなく、局所的な動きによる変形を行って得られるものでもないと判定した場合には、補正機能１７ａは、体動の影響がないものとみなし、ｋショット目の磁気共鳴データをそのままとする。

なお、ここでは、補正機能１７ａが、ｋショット目の低周波画像について、１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものであるか否かを判定したうえで、単純な平行移動又は回転を行って得られるものであると判定した場合に、計算した移動量に基づいて磁気共鳴データを補正することとしたが、実施形態はこれに限られない。

例えば、補正機能１７ａは、ｋショット目の低周波画像について、１ショット目の低周波画像に対して単純な平行移動又は回転を行って得られるものであるか否かの判定は行わずに、局所的な動きによる変形を行って得られるものであるか否かの判定のみを行う。そして、補正機能１７ａは、ｋショット目の低周波画像が１ショット目の低周波画像に対して局所的な動きによる変形を行って得られるものであると判定した場合には、上述した例と同様に、ｋショット目の磁気共鳴データを破棄する。一方、ｋショット目の低周波画像が１ショット目の低周波画像に対して局所的な動きによる変形を行って得られるものではないと判定した場合には、補正機能１７ａは、ｋショット目の磁気共鳴データに対して、回転量及び平行移動量を計算して補正する。この例では、回転量および平行移動量の計算を行うか否かを判定しないが、体動の影響がないショットの磁気共鳴データについては、補正量がゼロである場合とみなせばよい。

そして、補正機能１７ａは、２〜Ｎショット目の各低周波画像について上述した処理を行った後に、破棄されずに残ったショットの磁気共鳴データを、補正後の磁気共鳴データとして記憶回路１３に保存する。

例えば、Ｍショットの磁気共鳴データが破棄されたとすると、補正機能１７ａは、破棄されずに残った（Ｎ−Ｍ）ショット分のデータを補正後の磁気共鳴データとして記憶回路１３に保存する。このとき、例えば、補正機能１７ａは、中心領域Ｒ１に対応するデータについては、（Ｎ−Ｍ）ショット分のデータをアベレージング処理等によって加算処理する。また、補正機能１７ａは、辺縁領域Ｒ２に対応するデータについては、（Ｎ−Ｍ）ショット分のデータをそのまま補正後の磁気共鳴データとして用いる。

または、例えば、補正機能１７ａは、中心領域Ｒ１に対応するデータについて、（Ｎ−Ｍ）ショット分のデータの全てを加算処理するのではなく、（Ｎ−Ｍ）ショット分のデータのうち操作者によって指定されたショット数分のデータを加算処理してもよい。

例えば、操作者によって指定されたショット数がＬであったとすると、補正機能１７ａは、中心領域Ｒ１に対応するデータについては、Ｌショット分のデータを加算処理し、残りの（Ｎ−Ｍ−Ｌ）ショット分のデータは破棄する。この場合も、補正機能１７ａは、辺縁領域Ｒ２に対応するデータについては、（Ｎ−Ｍ）ショット分のデータをそのまま補正後の磁気共鳴データとして用いる。

または、例えば、補正機能１７ａは、中心領域Ｒ１に対応するデータについて、複数ショット部のデータを加算処理するのではなく、１ショット分のデータのみを用いてもよい。この場合も、補正機能１７ａは、辺縁領域Ｒ２に対応するデータについては、（Ｎ−Ｍ）ショット分のデータをそのまま補正後の磁気共鳴データとして用いる。

なお、マルチショット収集によってデータを収集した場合には、時間的に近いショットのデータほど、体動による移動が少ないと考えられる。そのため、例えば、補正機能１７ａは、上述したように、中心領域Ｒ１に対応するデータについて、操作者によって指定されたショット数分のデータを加算処理する場合、又は、１ショット分のデータのみを用いる場合には、相関度を計算する際の基準となる１ショット目に近いショットから順に、複数又は１つのデータを選択して用いるようにする。

図１に戻って、第１の生成機能１６ａは、収集機能１５ａによってフルサンプリングで収集されたｋ空間の中心領域に対応するデータからコイル感度マップを生成する。

具体的には、第１の生成機能１６ａは、補正機能１７ａによって記憶回路１３に保存された補正後の磁気共鳴データを読み出し、読み出した磁気共鳴データに含まれる中心領域に対応するデータに基づいて、コイル感度マップを生成する。

例えば、第１の生成機能１６ａは、局所用ＲＦコイル５としてフェーズドアレイコイルが用いられる場合に、フェーズドアレイコイルのコイルエレメントごとに、コイル感度マップを生成する。

第２の生成機能１６ｂは、補正機能１７ａによって体動の影響が補正された後の磁気共鳴データに基づいて、圧縮センシングにより画像を生成する。

具体的には、第２の生成機能１６ｂは、補正機能１７ａによって記憶回路１３に保存された補正後の磁気共鳴データを読み出し、読み出した磁気共鳴データと、第１の生成機能１６ａによって生成されたコイル感度マップとに基づいて、推定画像を生成する。

例えば、第２の生成機能１６ｂは、局所用ＲＦコイル５としてフェーズドアレイコイルが用いられる場合に、フェーズドアレイコイルのコイルエレメントごとに、対応する補正後の磁気共鳴データにフーリエ変換等の再構成処理を施すことで、画像を生成する。また、第２の生成機能１６ｂは、生成した各画像の折り返しをコイルエレメントごとに生成されたコイル感度マップに基づいて展開することによって、推定画像を生成する。

そして、第２の生成機能１６ｂは、生成した推定画像と、第１の生成機能１６ａによって生成されたコイル感度マップとを用いて、圧縮センシングにより最終画像を生成する。

例えば、第２の生成機能１６ｂは、コイル感度マップを用いて、ウェーブレット変換等によって推定画像からノイズが除去された画像を生成する処理を繰り返し行うことで、最終画像を生成する。

以上、処理回路１４〜１７が有する処理機能について説明したが、例えば、各処理回路は、プロセッサによって実現される。この場合に、各処理回路が有する処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路１３に記憶される。そして、各処理回路は、記憶回路１３から各プログラムを読み出して実行することで、各プログラムに対応する処理機能を実現する。換言すると、各プログラムを読み出した状態の各処理回路は、図１の各処理回路内に示された各機能を有することとなる。

図５は、第１の実施形態に係るＭＲＩ装置１００によって行われる処理の処理手順を示すフローチャートである。

例えば、図５に示すように、本実施形態では、収集機能１５ａが、操作者から開始の指示を受け付けた場合に（ステップＳ１０１，Ｙｅｓ）、マルチショット収集により磁気共鳴データを収集する（ステップＳ１０２）。例えば、ステップＳ１０１及びＳ１０２の処理は、処理回路１５が、収集機能１５ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。

その後、補正機能１７ａが、ショットごとに、ｋ空間の中心領域に対応するデータからリファレンスデータを生成する（ステップＳ１０３）。そして、補正機能１７ａは、リファレンスデータを複数のショット間で比較することで、磁気共鳴データにおける体動の影響を補正する（ステップＳ１０４）。例えば、ステップＳ１０３及びＳ１０４の処理は、処理回路１７が、補正機能１７ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。

その後、第１の生成機能１６ａが、補正後の中心領域のデータに基づいて、コイル感度マップを生成する（ステップＳ１０５）。例えば、ステップＳ１０５の処理は、処理回路１６が、第１の生成機能１６ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。

その後、第２の生成機能１６ｂが、補正後の磁気共鳴データ及びコイル感度マップに基づいて、推定画像を生成する（ステップＳ１０６）。そして、第２の生成機能１６ｂは、推定画像及びコイル感度マップに基づいて、圧縮センシングにより最終画像を生成する（ステップＳ１０７）。例えば、ステップＳ１０６及びＳ１０７の処理は、処理回路１６が、第２の生成機能１６ｂに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。

なお、ここでは、単一のプロセッサによって各処理回路が実現されるものとして説明したが、複数の独立したプロセッサを組み合わせて各処理回路を構成し、各プロセッサがプログラムを実行することによって各処理機能を実現するものとしてもよい。また、各処理回路が有する処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。また、図１に示す例では、単一の記憶回路１３が各処理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散して配置して、処理回路が個別の記憶回路から対応するプログラムを読み出す構成としても構わない。

上述したように、第１の実施形態では、収集機能１５ａが、各ショットにおいて、フルサンプリングされたｋ空間の中心領域Ｒ１のデータを収集する。これにより、画像の全体像を表す低周波成分のデータを得ることができ、データ収集時の被検体の体動をより正確に把握することが可能になる。

そして、第１の実施形態では、補正機能１７ａが、収集機能１５ａによってフルサンプリングで収集されたｋ空間の中心領域に対応するデータから生成されたリファレンスデータを用いて、磁気共鳴データにおける体動の影響を補正する。これにより、収集されたデータにおける体動の影響をより適切に補正することができるようになる。

また、第１の実施形態では、第２の生成機能１６ｂが、補正機能１７ａによって体動の影響が補正された後の磁気共鳴データに基づいて、圧縮センシングにより画像を生成する。これにより、被検体の体動にロバストな高速撮像を実現することができる。

また、第１の実施形態では、補正機能１７ａが、ｋ空間の中心領域に対応するデータを加算処理することによって、生成される画像のＳＮＲ（Signal-to-Noise Ratio）を向上させることができる。

また、第１の実施形態では、補正機能１７ａが、体動の影響が大きいショットのデータを破棄することによって、ｋ空間の中心領域の信号に対する体動の影響が低減されるため、圧縮センシングによる画像生成時のエラーを低減することができる。

（第１の実施形態の変形例）
なお、上述した第１の実施形態では、図２に示したように、収集機能１５ａが、マルチショット収集の各ショットにおいて、中心領域Ｒ１に対応するデータをフルサンプリングで収集し、２つの辺縁領域Ｒ２に対応するデータを間引きサンプリングで収集する場合の例を説明したが、実施形態はこれに限られない。

例えば、収集機能１５ａは、２つの辺縁領域Ｒ２について、ショットごとに、いずれか一方の辺縁領域Ｒ２に対応するデータのみを収集するようにしてもよい。すなわち、この場合には、２つの辺縁領域Ｒ２に対応するデータが、Ｋ_ro軸に対して非対称に収集されることになる。

図６は、第１の実施形態の変形例に係る収集機能１５ａによって行われる磁気共鳴データの収集の一例を示す図である。

例えば、図６に示すように、収集機能１５ａは、２つの辺縁領域Ｒ２について、ショットごとに、一方の辺縁領域Ｒ２に対応するデータと、他方の辺縁領域Ｒ２に対応するデータとを交互に収集する。

このとき、収集機能１５ａは、各辺縁領域Ｒ２について、図２に示した例と同様に、データをサンプリングするラインを複数のショットに所定数ずつ重複しないように割り当てて、データを収集する。ただし、本変形例では、収集機能１５ａは、図２に示した例のように、ショットごとに、２つの辺縁領域Ｒ２に対応するデータを収集する場合と比べて、少ない数のラインを各ショットに割り当てる。

これにより、本変形例では、図２に示した例と比べて、１ショットで収集するエコー信号の数、すなわち、ＥＴＬ（Echo Train Length）を減らすことができ、ＥＴＳ（Echo Train Space）を増加させたり、Ｔ２緩和補正（T2 decay補正）を行ったりしなくても、ＥＴＬの増加により生じる画像のボケ（T2 Blur）を低減させることができる。

なお、ここでは、収集機能１５ａが、２つの辺縁領域Ｒ２について、ショットごとに、いずれか一方の辺縁領域Ｒ２に対応するデータのみを収集することとしたが、例えば、収集していない方の辺縁領域Ｒ２をブランクのままにせず、収集した方の辺縁領域Ｒ２に対応するデータから、収集しなかった方の辺縁領域Ｒ２に対応するデータを推定して補填するようにしてもよい。

例えば、収集機能１５ａは、非収集部分が収集部分と複素共役の関係にあることを利用した画像再構成方法であるＡＦＩ（Asymmetric Fourier Imaging）を利用して、ショットごとに、収集した方の辺縁領域Ｒ２に対応するデータから、収集していない方の辺縁領域Ｒ２に対応するデータを補填してもよい。

（第２の実施形態）
また、上述した実施形態では、２Ｄ撮像が行われる場合の例を説明したが、実施形態はこれに限られない。例えば、上述したデータ収集の方法は、３Ｄ撮像が行われる場合にも適用することが可能である。

そこで、以下では、第２の実施形態として、ＭＲＩ装置１００によって３Ｄ撮像が行われる場合の例を説明する。なお、以下では、第２の実施形態について、第１の実施形態と異なる点を中心に説明することとし、第１の実施形態と共通する内容については詳細な説明を省略する。

図７及び８は、第２の実施形態に係る収集機能１５ａによって行われるデータ収集の一例を示す図である。また、図９は、第２の実施形態に係る比較例を示す図である。

例えば、図７に示すように、収集機能１５ａは、３Ｄ撮像を行う場合には、リードアウト方向の位置を示すＫ_ro軸と、位相エンコード方向の位置を示すＫ_pe軸と、スライス方向の位置を示すＫ_se軸とで表される３次元のｋ空間に対応するデータを収集する。

ここで、例えば、収集機能１５ａは、リードアウト方向に平行な軸を中心軸とした円筒形状の範囲をデータ収集領域とし、当該範囲内に、リードアウト方向、位相エンコード方向及びスライス方向それぞれに沿って、撮像条件として設定されたマトリクスサイズに応じた数のサンプリング点を設定する。

そして、収集機能１５ａは、マルチショット収集によって、位相エンコード方向及びスライス方向の位置を変えながら、リードアウト方向に沿ったラインごとにデータを収集する。具体的には、収集機能１５ａは、マルチショット収集によって、位相エンコード傾斜磁場及びスライス傾斜磁場の強度を変えながら複数のエコー信号を収集し、１つのエコー信号を収集するごとに、収集したエコー信号に基づいて、リードアウト方向に沿った１ライン分のデータをサンプリングする。

このように、マルチショット収集によって円筒形状のデータ収集領域のデータを収集する場合には、例えば、図９に示すように、データ収集領域を周方向に所定の角度ずつ分割することで、複数の領域Ｒを設定し、ショットごとに、１つの領域Ｒに対応するデータを収集する方法もあり得る。なお、図９では、複数の領域Ｒのうち、１つの領域Ｒのみを示し、他の領域Ｒについては図示を省略している。

これに対し、本実施形態では、例えば、図８に示すように、収集機能１５ａは、円筒形状のデータ収集領域に、中心軸を含む中心領域Ｒ１を設定し、さらに、中心領域Ｒ１以外の範囲を周方向に所定の角度ずつ分割することで、複数の辺縁領域Ｒ２を設定する。なお、図８でも、複数の辺縁領域Ｒ２のうち、１つの辺縁領域Ｒ２のみを示し、他の辺縁領域Ｒ２については図示を省略している。

そして、収集機能１５ａは、マルチショット収集の各ショットにおいて、中心領域Ｒ１に対応するデータをフルサンプリングで収集し、１つの辺縁領域Ｒ２に対応するデータを間引きサンプリングで収集する。

例えば、図８に示すように、収集機能１５ａは、辺縁領域Ｒ２については、各ショットにおいて、領域内に含まれるリードアウト方向に沿った複数のラインを位相エンコード方向及びスライスエンコード方向に間引いてデータを収集する。例えば、収集機能１５ａは、辺縁領域Ｒ２内で、データ収集領域の中心に近い側から遠い側に向けて、径方向に間隔を空けるようにラインを移動させながらデータを収集する。

一方、中心領域Ｒ１については、収集機能１５ａは、各ショットにおいて、領域内に含まれる全てのラインのデータを収集する。例えば、収集機能１５ａは、中心領域Ｒ１内で、データ収集領域の中心軸から径方向の外側に向けて、径方向に隙間を空けないようにラインを移動させながらデータを収集する。これにより、本実施形態では、複数のショットが行われることで、ショットごとに、フルサンプリングされた中心領域Ｒ１のデータが得られることになる。

このように、３Ｄ撮像が行われる場合でも、収集機能５ａが、各ショットにおいて、フルサンプリングされたｋ空間の中心領域Ｒ１のデータを収集することによって、ショットごとに、画像の全体像を表す低周波成分のデータを得ることができる。このデータによれば、データ収集時の被検体の体動をより正確に把握することが可能になる。

なお、本実施形態では、補正機能１７ａ、第１の生成機能１６ａ及び第２の生成機能１６ｂは、処理対象のデータや画像が２次元ではなく３次元になるが、第１の実施形態で説明した処理と同様の処理を行う。

上述したように、第２の実施形態では、３Ｄ撮像が行われる場合に、収集機能１５ａが、各ショットにおいて、フルサンプリングされたｋ空間の中心領域Ｒ１のデータを収集する。これにより、３Ｄ撮像が行われる場合でも、画像の全体像を表す低周波成分のデータを得ることができ、データ収集時の被検体の体動をより正確に把握することが可能になる。また、第２の実施形態でも、第１の実施形態と同様に、収集されたデータにおける体動の影響をより適切に補正することができるようになる。また、第２の実施形態でも、第１の実施形態と同様に、被検体の体動にロバストな高速撮像を実現することができる。

なお、上述した各実施形態では、収集機能１５ａが辺縁領域Ｒ２に対応するデータを収集する際に行う間引きサンプリングが、ランダムにデータを間引くサンプリングである場合の例を説明したが、実施形態はこれに限られない。例えば、収集機能１５ａが辺縁領域Ｒ２に対応するデータを収集する際に行う間引きサンプリングは、周期的にデータを間引くサンプリング（アンダーサンプリング）であってもよい。一般的に、圧縮センシングはランダムに間引きされたデータを用いて行われる手法であるため、この場合には、第２の生成機能１６ｂは、周期的に間引かれたデータから画像を生成することが可能な画像生成法を用いて、画像を生成する。

また、上述した各実施形態では、本明細書における収集部、補正部、第１の生成部及び第２の生成部を、それぞれ、処理回路の収集機能、補正機能、第１の生成機能及び第２の生成機能によって実現する場合の例を説明したが、実施形態はこれに限られない。例えば、本明細書における収集部、補正部、第１の生成部及び第２の生成部は、実施形態で述べた収集機能、補正機能、第１の生成機能及び第２の生成機能によって実現する他にも、ハードウェアのみ、ソフトウェアのみ、又は、ハードウェアとソフトウェアとの混合によって同機能を実現するものであっても構わない。

また、上述した説明で用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、或いは、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等の回路を意味する。プロセッサは、記憶回路に保存されたプログラムを読み出して実行することで、機能を実現する。なお、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合は、プロセッサは回路内に組み込まれたプログラムを読み出して実行することで機能を実現する。また、本実施形態のプロセッサは、単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて一つのプロセッサとして構成され、その機能を実現するようにしてもよい。

ここで、プロセッサによって実行されるプログラムは、ＲＯＭ（Read Only Memory）や記憶回路等に予め組み込まれて提供される。なお、このプログラムは、これらの装置にインストール可能な形式又は実行可能な形式のファイルでＣＤ（Compact Disk）−ＲＯＭ、ＦＤ（Flexible Disk）、ＣＤ−Ｒ（Recordable）、ＤＶＤ（Digital Versatile Disk）等のコンピュータで読み取り可能な記憶媒体に記録されて提供されてもよい。また、このプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納され、ネットワーク経由でダウンロードされることにより提供又は配布されてもよい。例えば、このプログラムは、上述した各機能部を含むモジュールで構成される。実際のハードウェアとしては、ＣＰＵが、ＲＯＭ等の記憶媒体からプログラムを読み出して実行することにより、各モジュールが主記憶装置上にロードされて、主記憶装置上に生成される。

以上説明した少なくとも１つの実施形態によれば、データ収集時の被検体の体動をより正確に把握することが可能なデータを得ることができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これらの実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これらの実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１００ ＭＲＩ装置
１５ 処理回路
１５ａ 収集機能
１６ 処理回路
１６ａ 第１の生成機能
１６ｂ 第２の生成機能
１７ 処理回路
１７ａ 補正機能

Claims (18)

1. １回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集する磁気共鳴データ収集方法であって、
   前記１ショットで収集される磁気共鳴データは、ｋ空間の中心領域に対応するデータと、前記ｋ空間の辺縁領域に対応するデータとを含み、
   前記１ショットの収集において、前記中心領域に対応するデータはフルサンプリングで収集され、前記辺縁領域に対応するデータは間引きサンプリングで収集される、
   磁気共鳴データ収集方法。
2. 前記間引きサンプリングは、ランダムにデータを間引くサンプリングである、
   請求項１に記載の磁気共鳴データ収集方法。
3. 前記間引きサンプリングは、周期的にデータを間引くサンプリングである、
   請求項１に記載の磁気共鳴データ収集方法。
4. 前記フルサンプリングで収集された前記中心領域に対応するデータからコイル感度マップを生成することをさらに含む、
   請求項１〜３のいずれか一つに記載の磁気共鳴データ収集方法。
5. 前記フルサンプリングで収集された前記中心領域に対応するデータから生成されたリファレンスデータを用いて、前記磁気共鳴データにおける体動の影響を補正することをさらに含む、
   請求項１〜４のいずれか一つに記載の磁気共鳴データ収集方法。
6. 前記リファレンスデータを複数のショット間で比較することで、前記体動の影響を補正する、
   請求項５に記載の磁気共鳴データ収集方法。
7. 前記体動の影響の補正は、体動の影響が大きいリファレンスデータに対応するショットのデータを破棄することを含む、
   請求項６に記載の磁気共鳴データ収集方法。
8. 前記体動の影響が補正された後の磁気共鳴データに基づいて、圧縮センシングにより画像を生成することをさらに含む、
   請求項５〜７のいずれか一つに記載の磁気共鳴データ収集方法。
9. ＦＳＥ（Fast Spin Echo）系のパルスシーケンス、又は、ＦＦＥ（Fast Field Echo）系のパルスシーケンスを用いて、前記磁気共鳴データを収集する、
   請求項１〜８のいずれか一つに記載の磁気共鳴データ収集方法。
10. １回のＲＦ励起で複数のエコー信号を収集する１ショットの収集を繰り返し行って磁気共鳴データを収集する収集部を備え、
    前記１ショットで収集される磁気共鳴データは、ｋ空間の中心領域に対応するデータと、前記ｋ空間の辺縁領域に対応するデータとを含み、
    前記１ショットの収集において、前記中心領域に対応するデータはフルサンプリングで収集され、前記辺縁領域に対応するデータは間引きサンプリングで収集される、
    磁気共鳴イメージング装置。
11. 前記間引きサンプリングは、ランダムにデータを間引くサンプリングである、
    請求項１０に記載の磁気共鳴イメージング装置。
12. 前記間引きサンプリングは、周期的にデータを間引くサンプリングである、
    請求項１０に記載の磁気共鳴イメージング装置。
13. 前記フルサンプリングで収集された前記中心領域に対応するデータからコイル感度マップを生成する第１の生成部をさらに備える、
    請求項１０〜１２のいずれか一つに記載の磁気共鳴イメージング装置。
14. 前記フルサンプリングで収集された前記中心領域に対応するデータから生成されたリファレンスデータを用いて、前記磁気共鳴データにおける体動の影響を補正する補正部をさらに備える、
    請求項１０〜１３のいずれか一つに記載の磁気共鳴イメージング装置。
15. 前記補正部は、前記リファレンスデータを複数のショット間で比較することで、前記体動の影響を補正する、
    請求項１４に記載の磁気共鳴イメージング装置。
16. 前記体動の影響の補正は、体動の影響が大きいリファレンスデータに対応するショットのデータを破棄することを含む、
    請求項１５に記載の磁気共鳴イメージング装置。
17. 前記体動の影響が補正された後の磁気共鳴データに基づいて、圧縮センシングにより画像を生成する第２の生成部をさらに備える、
    請求項１４〜１６のいずれか一つに記載の磁気共鳴イメージング装置。
18. 前記収集部は、ＦＳＥ（Fast Spin Echo）系のパルスシーケンス、又は、ＦＦＥ（Fast Field Echo）系のパルスシーケンスを用いて、前記磁気共鳴データを収集する、
    請求項１０〜１７のいずれか一つに記載の磁気共鳴イメージング装置。

Images