(57)【特許請求の範囲】
【請求項１】
第１の範囲内に並んだ複数のｋ空間データである複数の第１のデータに対して、前記第１の範囲と隣接する第２の範囲内に並んだ複数のｋ空間データである複数の第２のデータを付加するデータ付加部であって、当該複数の第１のデータ及び当該複数の第２のデータからなる複数の第３のデータが前記第１の範囲と前記第２の範囲とが並ぶ方向に並べられた場合の当該第３のデータに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有するデータ付加部と、
前記複数の第３のデータに対して再構成処理を行い、複数のＭＲ画像を生成する生成部と、
前記生成部により生成された複数のＭＲ画像から、前記第２の範囲に対応するＭＲ画像を削除する削除部と、
を備える、磁気共鳴イメージング装置。
【請求項２】
前記複数の第１のデータを記憶する記憶部を更に備え、
前記データ付加部は、前記記憶部に記憶された前記複数の第１のデータに基づいて前記複数の第２のデータを生成し、生成した当該複数の第２のデータを当該複数の第１のデータに付加する、請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記データ付加部は、生成対象である前記複数の第２のデータに含まれる複数の第２のフレームのうち、前記第１の範囲に最も近い第２のフレームについて、前記複数の第１のデータに含まれる複数の第１のフレームのうち、当該第２のフレームに最も近い第１のフレームを設定し、当該第２のフレームに対応するサンプリングパターンで収集された前記第１のフレームの中から、対応するＭＲ画像の信号強度について、設定した第１のフレームと類似する第１のフレームを、当該第２のフレームとして生成する、請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項４】
前記データ付加部は、前記複数の第２のフレームのうち、前記第１の範囲から最も遠い第２のフレームについて、前記複数の第１のフレームのうち、当該第２のフレームから最も遠い第１のフレームを設定し、当該第２のフレームに対応するサンプリングパターンで収集された前記第１のフレームの中から、対応するＭＲ画像の信号強度について、設定した第１のフレームと類似する第１のフレームを、当該第２のフレームとして生成する、請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
前記データ付加部は、前記複数の第２のフレームのうち、前記第１の範囲に最も近い第２のフレーム及び前記第１の範囲から最も遠い第２のフレーム以外の第２のフレームについて、前記複数の第１のデータに含まれる複数の第１のフレームのうち、前記第１の範囲に最も近い第２のフレームに最も近い第１のフレームに対応するＭＲ画像の信号強度と、前記第１の範囲から最も遠い第２のフレームから最も遠い第１のフレームに対応するＭＲ画像の信号強度とに基づいて信号強度を算出することにより、当該第２のフレームを生成する、請求項４に記載の磁気共鳴イメージング装置。
【請求項６】
前記データ付加部は、設定した第１のフレームに対してフーリエ変換を適用することにより得られた折り返しを含む第１の画像空間と、前記第１の範囲に最も近い第２のフレーム又は前記第１の範囲から最も遠い第２のフレームに対応するサンプリングパターンで収集された複数の前記第１のフレームのそれぞれに対して前記フーリエ変換を適用することにより得られた折り返しを含む第２の画像空間のそれぞれとの類似度を算出し、算出した類似度が示す類似性が最も高い第２の画像空間に対応する前記第１のフレームを、当該第２のフレームとして生成する、請求項４又は５に記載の磁気共鳴イメージング装置。
【請求項７】
前記第１のデータを収集する収集部を更に備え、
前記データ付加部は、前記収集部による前記第１のデータの収集に加えて、前記第２のデータを収集し、収集した前記第２のデータを前記第１のデータに付加する、請求項１に記載の磁気共鳴イメージング装置。
【請求項８】
第１の範囲内に並んだ複数のｋ空間データである複数の第１のデータ、及び、前記第１の範囲と隣接する第２の範囲内に並んだ複数のｋ空間データである複数の第２のデータからなる複数の第３のデータに対して、当該複数の第３のデータが前記第１の範囲と前記第２の範囲とが並ぶ方向に並べられた場合の当該第３のデータに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有するように窓関数を適用する窓関数適用部と、
前記窓関数適用部により窓関数が適用された前記複数の第３のデータに対して再構成処理を行い、複数のＭＲ画像を生成する生成部と、
前記生成部により生成された複数のＭＲ画像から、前記第２の範囲に対応するＭＲ画像を削除する削除部と、
を備える、磁気共鳴イメージング装置。
【請求項９】
第１の範囲内に並んだ複数のｋ空間データである複数の第１のデータに対して、前記第１の範囲と隣接する第２の範囲内に並んだ複数のｋ空間データである複数の第２のデータを付加するデータ付加部であって、当該複数の第１のデータ及び当該複数の第２のデータからなる複数の第３のデータが前記第１の範囲と前記第２の範囲とが並ぶ方向に並べられた場合の当該第３のデータに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有するように付加するデータ付加部と、
前記複数の第３のデータに対して再構成処理を行い、複数のＭＲ画像を生成する生成部と、
前記生成部により生成された複数のＭＲ画像から、前記第２の範囲に対応するＭＲ画像を削除する削除部と、
を備える、画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージング装置及び画像処理装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置（ＭＲＩ（Magnetic  Resonance  Imaging）装置）は、核磁気共鳴現象を利用して被検体内部の情報を画像化するための装置である。ＭＲＩ装置は、コイルを利用して、被検体内部にある特定の原子核（例えば、水素原子の原子核）からの核磁気共鳴信号をサンプリングすることで、ｋ空間データと呼ばれるデータを収集し、ｋ空間データにフーリエ変換を適用することで、ＭＲ（Magnetic  Resonance）画像を得る。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】J.  Tsao,  P.  Boesiger,  K.  P.  Pruessmann,  ``k-t  BLAST  and  k-t  SENSE:  Dynamic  MRI  With  High  Frame  Rate  Exploiting  Spatiotemporal  Correlations,’’  Magnetic  Resonance  in  Medicine  50:1031-1042  (2003).
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、変換対象のデータが周期的に並んでいると仮定した上で変換を施す処理を行う際に発生し得る特有のアーティファクトの発生を抑制することができる磁気共鳴イメージング装置及び画像処理装置を提供することである。
【課題を解決するための手段】
【０００５】
実施形態の磁気共鳴イメージング装置は、データ付加部と、生成部と、削除部とを備える。データ付加部は、第１の範囲内に並んだ複数のｋ空間データである複数の第１のデータに対して、第１の範囲と隣接する第２の範囲内に並んだ複数のｋ空間データである複数の第２のデータを付加する。複数の第１のデータ及び複数の第２のデータからなる複数の第３のデータが第１の範囲と第２の範囲とが並ぶ方向に並べられた場合の当該第３のデータに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有する。生成部は、複数の第３のデータに対して再構成処理を行い、複数のＭＲ画像を生成する。削除部は、生成部により生成された複数のＭＲ画像から、第２の範囲に対応するＭＲ画像を削除する。
【図面の簡単な説明】
【０００６】
【図１】図１は、第１の実施形態に係るＭＲＩ装置を示すブロック図。
【図２】図２は、第１の実施形態における処理手順を示すフローチャート。
【図３】図３は、第１の実施形態に係る収集位置の一例を示す図。
【図４】図４は、第１の実施形態における画像生成ステップＳ１０４の処理手順を示すフローチャート。
【図５】図５は、第１の実施形態に係る再構成用ｋ空間データの一例を示す図。
【図６】図６は、第１の実施形態に係るステップＳ２０１の詳細な処理手順を示すフローチャート。
【図７】図７は、第１の実施形態に係る再構成用ｋ空間データに対応するＭＲ画像の信号強度の一例を示す図。
【図８】図８は、第２の実施形態に係るＭＲＩ装置の構成を示す機能ブロック図。
【図９】図９は、第２の実施形態における処理手順を示すフローチャート。
【図１０】図１０は、第２の実施形態におけるステップＳ４０５の処理手順を示すフローチャート。
【図１１】図１１は、第３の実施形態に係るＭＲＩ装置の構成を示す機能ブロック図。
【図１２】図１２は、第３の実施形態における処理手順を示すフローチャート。
【図１３】図１３は、ステップＳ６０５の処理の流れの一例を示すフローチャート。
【図１４】図１４は、第４の実施形態に係る画像処理装置の構成を示す機能ブロック図。
【発明を実施するための形態】
【０００７】
以下、各実施形態に係る磁気共鳴イメージング装置及び画像処理装置を詳細に説明する。
【０００８】
（第１の実施形態）
図１は、第１の実施形態に係るＭＲＩ装置１００を示す機能ブロック図である。図１に示すように、ＭＲＩ装置１００は、静磁場磁石１０１と、傾斜磁場コイル１０２と、傾斜磁場電源１０３と、寝台１０４と、寝台制御回路１０５と、送信コイル１０６と、送信回路１０７と、受信コイルアレイ１０８と、受信回路１０９と、シーケンス制御回路１１０と、計算機システム１２０とを備える。なお、ＭＲＩ装置１００に被検体Ｐ（例えば、人体）は含まれない。
【０００９】
静磁場磁石１０１は、中空の円筒形状（円筒の軸に直交する断面が楕円状となるものを含む）に形成された磁石であり、内部の空間に一様な静磁場を発生する。静磁場磁石１０１は、例えば、永久磁石、超伝導磁石などである。
【００１０】
傾斜磁場コイル１０２は、中空の円筒形状（円筒の軸に直交する断面が楕円状となるものを含む）に形成されたコイルであり、静磁場磁石１０１の内側に配置される。傾斜磁場コイル１０２は、互いに直交するＸ，Ｙ，Ｚの各軸に対応する３つのコイルが組み合わされて形成されており、これら３つのコイルは、傾斜磁場電源１０３から個別に電流の供給を受けて、Ｘ、Ｙ、Ｚの各軸に沿って磁場強度が変化する傾斜磁場を発生する。ここで、傾斜磁場コイル１０２によって発生するＸ，Ｙ，Ｚ各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ及びリードアウト用傾斜磁場Ｇｒにそれぞれ対応する。スライス選択用傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Ｇｅは、空間的位置に応じてＭＲ信号の位相を変化させるために利用される。リードアウト用傾斜磁場Ｇｒは、空間的位置に応じてＭＲ信号の周波数を変化させるために利用される。
【００１１】
傾斜磁場電源１０３は、傾斜磁場コイル１０２に電流を供給する。例えば、傾斜磁場電源１０３は、傾斜磁場コイル１０２を形成する３つのコイルのそれぞれに、個別に電流を供給する。
【００１２】
寝台１０４は、被検体Ｐが載置される天板１０４ａを備え、寝台制御回路１０５による制御のもと、天板１０４ａを、被検体Ｐが載置された状態で傾斜磁場コイル１０２の空洞（撮像口）内へ挿入する。通常、寝台１０４は、長手方向が静磁場磁石１０１の中心軸と平行になるように設置される。
【００１３】
寝台制御回路１０５は、計算機システム１２０による制御のもと、寝台１０４を駆動して天板１０４ａを長手方向及び上下方向へ移動する。
【００１４】
送信コイル１０６は、傾斜磁場コイル１０２の内側に配置され、送信回路１０７からＲＦパルスの供給を受けて、高周波磁場を発生する。
【００１５】
送信回路１０７は、対象とする原子の種類及び磁場の強度で決まるラーモア周波数に対応するＲＦ（Radio  Frequency）パルスを送信コイル１０６に供給する。
【００１６】
受信コイルアレイ１０８は、傾斜磁場コイル１０２の内側に配置され、高周波磁場の影響によって被検体Ｐから発せられる磁気共鳴信号（以下、ＭＲ信号と称する）を受信する。受信コイルアレイ１０８は、ＭＲ信号を受信すると、受信したＭＲ信号を受信回路１０９へ出力する。なお、第１の実施形態において、受信コイルアレイ１０８は、１以上、典型的には複数の受信コイルを有するコイルアレイである。
【００１７】
受信回路１０９は、受信コイルアレイ１０８から出力されるＭＲ信号に基づいてＭＲデータを生成する。例えば、受信回路１０９は、受信コイルアレイ１０８から出力されるＭＲ信号をデジタル変換することによってＭＲデータを生成する。また、受信回路１０９は、生成したＭＲデータをシーケンス制御回路１１０へ送信する。
【００１８】
なお、受信回路１０９は、静磁場磁石１０１や傾斜磁場コイル１０２などを備える架台装置側に備えられていてもよい。ここで、第１の実施形態において、受信コイルアレイ１０８の各コイルエレメント（各受信コイル）から出力されるＭＲ信号は、適宜分配・合成されることで、チャネルなどと呼ばれる単位で受信回路１０９に出力される。このため、ＭＲデータは、受信回路１０９以降の後段の処理においてチャネル毎に取り扱われる。コイルエレメントの総数とチャネルの総数との関係は、同一の場合もあれば、コイルエレメントの総数に対してチャネルの総数が少ない場合、あるいは反対に、コイルエレメントの総数に対してチャネルの総数が多い場合もある。以下において、「チャネル毎」のように表記する場合、その処理が、コイルエレメント毎に行われてもよいし、あるいは、コイルエレメントが分配・合成されたチャネル毎に行われてもよいことを示す。なお、分配・合成のタイミングは、上述したタイミングに限られるものではない。ＭＲ信号若しくはＭＲデータは、後述する再構成処理の前までに、チャネル単位に分配・合成されればよい。
【００１９】
シーケンス制御回路１１０は、計算機システム１２０から送信されるシーケンス情報に基づいて、傾斜磁場電源１０３、送信回路１０７及び受信回路１０９を駆動することによって、被検体Ｐの撮像を行う。例えば、シーケンス制御回路１１０は、プロセッサにより実現される。ここで、シーケンス情報は、撮像を行うための手順を定義した情報である。シーケンス情報には、傾斜磁場電源１０３が傾斜磁場コイル１０２に供給する電源の強さや電源を供給するタイミング、送信回路１０７が送信コイル１０６に送信するＲＦパルスの強さやＲＦパルスを印加するタイミング、受信回路１０９がＭＲ信号を検出するタイミングなどが定義される。シーケンス制御回路１１０は、プロセッサにより実現される。
【００２０】
なお、シーケンス制御回路１１０は、傾斜磁場電源１０３、送信回路１０７及び受信回路１０９を駆動して被検体Ｐを撮像した結果、受信回路１０９からＭＲデータを受信すると、受信したＭＲデータを計算機システム１２０へ転送する。
【００２１】
計算機システム１２０は、ＭＲＩ装置１００の全体制御や、データ収集、画像生成などを行う。計算機システム１２０は、インタフェース回路１２１、画像生成回路１２２、記憶回路１２３、入力回路１２４、ディスプレイ１２５及び処理回路１２６を有する。
【００２２】
インタフェース回路１２１は、シーケンス情報をシーケンス制御回路１１０へ送信し、シーケンス制御回路１１０からＭＲデータを受信する。また、インタフェース回路１２１は、ＭＲデータを受信すると、受信したＭＲデータを記憶回路１２３に格納する。記憶回路１２３に格納されたＭＲデータは、処理回路１２６によってｋ空間に配置される。この結果、記憶回路１２３は、複数チャネル分のｋ空間データを記憶する。
【００２３】
処理回路１２６は、ＭＲＩ装置１００の全体制御を行う。具体的には、処理回路１２６は、入力回路１２４を介して操作者から入力される撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御回路１１０に送信することによって撮像を制御する。また、処理回路１２６は、撮像の結果としてシーケンス制御回路１１０から送られるＭＲデータに基づいて行われるＭＲ画像の生成を制御したり、ディスプレイ１２５による表示を制御したりする。処理回路１２６は、プロセッサにより実現される。処理回路１２６は、収集機能１２６ａを有する。
【００２４】
収集機能１２６ａは、入力回路１２４を介して操作者から入力される撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御回路１１０に送信することによって撮像を制御する。また、収集機能１２６ａは、空間方向のサンプリング位置を変化させながら、複数のチャネルについて、時系列のｋ空間データを収集する。なお、時系列のｋ空間データは、複数のｋ空間データである。そして、収集機能１２６ａは、収集した時系列のｋ空間データを記憶回路１２３に格納する。
【００２５】
ここで、例えば、処理回路１２６の構成要素である収集機能１２６ａは、コンピュータによって実行可能なプログラムの形態で記憶回路１２３に記憶されている。処理回路１２６は、収集機能１２６ａに対応するプログラムを記憶回路１２３から読み出し、読み出したプログラムを実行することで、収集機能１２６ａを実現する。換言すると、収集機能１２６ａに対応するプログラムを読み出した状態の処理回路１２６は、図１の処理回路１２６内に示された収集機能１２６ａを有することとなる。
【００２６】
なお、収集機能１２６ａは、特許請求の範囲に記載された収集部の一例である。
【００２７】
画像生成回路１２２は、収集機能１２６ａで収集されて記憶回路１２３に記憶された時系列のｋ空間データを用いて、ＭＲ画像の生成を行う。画像生成回路１２２は、プロセッサにより実現される。画像生成回路１２２は、データ付加機能１２２ａ、生成機能１２２ｂ、削除機能１２２ｃを有する。
【００２８】
データ付加機能１２２ａは、収集機能１２６ａにより収集されて記憶回路１２３に記憶されたｋ空間データから、複数の付加ｋ空間データを生成し、収集機能１２６ａによって収集されたｋ空間データに付加する。
【００２９】
生成機能１２２ｂは、データ付加機能１２２ａにより付加ｋ空間データが付加された各チャネルの時系列のｋ空間データに対して、時空間方向に逆フーリエ変換を適用することで、各チャネルの変換後時空間データを得る。そして、生成機能１２２ｂは、各チャネルの変換後時空間データと、各チャネルの感度分布情報とを用いて、アンフォールドを行ってアンフォールドデータを得る。そして、生成機能１２２ｂは、アンフォールドデータに対して、時間軸方向にフーリエ変換を適用する。
【００３０】
削除機能１２２ｃは、データ付加機能１２２ａで付加された付加ｋ空間データに対応する時間に該当するアンフォールドデータを削除する。
【００３１】
なお、画像生成回路１２２によって得られたＭＲ画像は、必要に応じてディスプレイ１２５で表示されたり、記憶回路１２３に送られて記憶されたりする。画像生成回路１２２は、プロセッサにより実現される。
【００３２】
ここで、例えば、画像生成回路１２２の構成要素であるデータ付加機能１２２ａ、生成機能１２２ｂ及び削除機能１２２ｃの各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶回路１２３に記憶されている。画像生成回路１２２は、各プログラムを記憶回路１２３から読み出し、読み出した各プログラムを実行することで、各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の画像生成回路１２２は、図１の画像生成回路１２２内に示された各機能を有することとなる。なお、図１においては、単一の画像生成回路１２２にて、データ付加機能１２２ａ、生成機能１２２ｂ及び削除機能１２２ｃの各処理機能が実現されるものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路１２６を構成し、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。
【００３３】
なお、データ付加機能１２２ａは、特許請求の範囲に記載されたデータ付加部の一例である。また、生成機能１２２ｂは、特許請求の範囲に記載された生成部の一例である。また、削除機能１２２ｃは、特許請求の範囲に記載された削除部の一例である。
【００３４】
記憶回路１２３は、インタフェース回路１２１によって受信されたＭＲデータや、収集機能１２６ａによってｋ空間に配置された時系列データ（ｋ−ｔ空間データ）、生成機能１２２ｂによって生成されたＭＲ画像の画像データなどを記憶する。また、記憶回路１２３は、各種のプログラムを記憶する。記憶回路１２３は、例えば、ＲＡＭ（Random  Access  Memory）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。
【００３５】
入力回路１２４は、医師や診療放射線技師等の操作者からの各種指示や情報入力を受け付ける。入力回路１２４は、例えば、トラックボール、スイッチボタン、マウス、キーボード等によって実現される。入力回路１２４は、処理回路１２６に接続されており、操作者から受け取った入力操作を電気信号に変換して処理回路１２６へと出力する。
【００３６】
ディスプレイ１２５は、処理回路１２６による制御のもと、各種ＧＵＩ（Graphical  User  Interface）や、生成機能１２２ｂによって生成されたＭＲ画像等を表示する。
【００３７】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central  Preprocess  Unit）、ＧＰＵ(Graphics  Processing  Unit)、或いは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ））、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））等の回路を意味する。なお、記憶回路１３２にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。
【００３８】
ところで、ＭＲＩ装置において核磁気共鳴信号（ＭＲ信号）は、１次元データとしてサンプリングされる。そこでＭＲＩ装置では、２次元あるいは３次元のＭＲ画像を得るために、ｋ空間上で１次元のサンプリングを繰り返し行い、ＭＲ画像の生成に必要なデータを収集する。出力対象のＭＲ画像と同じ解像度（フルサンプリング）でｋ空間データがサンプリングされれば、得られたｋ空間データにフーリエ変換を適用することで、ＭＲ画像の生成が可能である。
【００３９】
ＭＲＩ装置のサンプリングには時間がかかることが知られている。また、ＭＲＩ装置を用いた時系列データの撮像においてサンプリングに時間がかかると、撮像速度が低下し、対象物の動きに追従できなくなる。これに対し、高速な撮像を実現するための様々な技術が研究、開発されており、その１つとして、パラレルイメージング（ＰＩ（Parallel  Imaging））と呼ばれる技術がある。パラレルイメージングでは、フルサンプリングに対して間引いたサンプリングパターンを用い、複数のコイルを用いて複数のｋ空間データを収集する。間引きされた複数のｋ空間データにそのままフーリエ変換を適用するとエイリアシングが生じるが、パラレルイメージングでは、コイルの物理的な配置によって生じる感度の違いを利用することで、間引きによるエイリアシングのないＭＲ画像を生成する。
【００４０】
パラレルイメージングの１種であるＳＥＮＳＥ（sensitivity  encoding）の手順は次のとおりである。予め、準備スキャン等によって各コイルの感度分布情報を収集しておく。次に、ＭＲＩ装置は、イメージングスキャンで得られた各コイルの間引きｋ空間データを用いてフーリエ変換を行い、各コイルのＭＲ画像を得る。その後、ＭＲＩ装置は、各コイルのＭＲ画像を、「真のＭＲ画像とコイルの感度分布との積を折り返し位置の分だけ加算して得られた画像である」と考えて、真のＭＲ画像を推定する。
【００４１】
ＳＥＮＳＥを拡張した手法の１種であるｋ−ｔ  ＳＥＮＳＥでは、時系列のＭＲ画像に対して時間軸方向にフーリエ変換を適用した空間上で、ＳＥＮＳＥによるエイリアシングの除去処理を実行して、エイリアシングが除去されたスペクトルデータを生成する。そして、エイリアシングが除去されたスペクトルデータに対して逆フーリエ変換を適用して、時系列のＭＲ画像を生成する。なお、フーリエ変換は、変換対象のデータが周期的に並んでいると仮定した上で、変換対象のデータに変換を施す処理の一例である。
【００４２】
なお、パラレルイメージングにおいて、フルサンプリングのサンプル数を間引きサンプリング時のサンプル数で除した値は、Ｒ（Reduction  Factor）と呼ばれる。
【００４３】
ｋ−ｔ  ＳＥＮＳＥのように時間軸方向にフーリエ変換を適用する手法（以降、時空間再構成法）は、時間軸方向に時系列のＭＲ画像が周期的に並んでいることを仮定した上で、変換を施す処理を行う。そのため、時間軸方向の始めと終わりとで、対応するＭＲ画像の信号強度が大きく異なるような時系列のｋ空間データを収集し、収集した時系列のｋ空間データに対応する複数のＭＲ画像（時系列のＭＲ画像）が周期的に並んでいると仮定した場合には、周期的に並んでいると仮定された場合の複数のＭＲ画像において、信号強度が大きく変化する箇所が存在する。このような場合には、周期的に並んでいると仮定された場合の複数のＭＲ画像の信号強度を表す関数は、連続性を有さない。そのため、生成されるＭＲ画像に特有のアーティファクトが発生し、再構成精度が低下する場合がある。ここで、時間軸方向の始めと終わりとで、対応するＭＲ画像の信号強度が大きく異なるような時系列のｋ空間データとしては、例えば、パフュージョンイメージングやインバージョンパルス等により収集された時系列のｋ空間データが挙げられる。また、ｋ空間データに対応するＭＲ画像とは、ｋ空間データに対して再構成処理が行われて生成されたＭＲ画像を指す。また、上述の信号強度を表す関数は、例えば、複数のｋ空間データの収集時刻又は収集順序を入力値とし、その入力値に対応するｋ空間データに対して再構成処理が行われて生成されたＭＲ画像内の注目位置における信号強度を出力値とする関数である。
【００４４】
このようなことから、本実施形態では、周期的に並んでいると仮定された場合に、対応する関数が連続性を有しないような複数のｋ空間データに対して、連続性を有するようになるｋ空間データを付加することで、時空間再構成法を適用した際に、アーティファクトの発生を抑制して、再構成精度の低下を抑制するようにする。
【００４５】
以下の実施形態において、「ｘ空間」とは、水素分布画像空間（便宜上、本実施形態では水素原子以外を撮像の対象とする場合も含めて水素分布画像と呼ぶ）のことであり、「ｘ空間データ」とは、ｘ空間内における信号点の集合のことである。ｘ空間上の異なる信号点は、ｘ空間上の異なる位置の信号点と対応付けられる。例えば、脳の３次元の水素分布画像は、３次元のｘ空間データであり、心臓のある断面を撮像した２次元のｘ空間画像は、２次元のｘ空間データである。
【００４６】
また、「ｋ−ｔ空間データ」とは、時系列のｋ空間データである。ｋ−ｔ空間データに対して空間方向にフーリエ変換を適用すると、時系列のｘ空間データ（ｘ−ｔ空間データ）が得られ、この時系列のｘ空間データに対して時間軸方向にフーリエ変換を適用すると、ｘ−ｆ空間データが得られる。即ち、ｘ−ｆ空間データは、ｘ−ｔ空間データがフーリエ変換後の係数値で表現されたものである。ｘ−ｆ空間データ内における各信号点は、ｘ軸位置とｆ軸位置とが定まれば、その位置が定まる。なお、「ｆ」は、周波数（frequency）に対応する。
【００４７】
以下の実施形態において、生成機能１２２ｂは、各チャネルのｋ−ｔ空間データに対して空間方向にフーリエ変換を適用するとともに時間軸方向にフーリエ変換を適用することで、各チャネルのｘ−ｆ空間データを得る。また、生成機能１２２ｂは、チャネル数分のｘ−ｆ空間データとチャネル数分の感度分布情報とを用いてアンフォールドを行うことで、エイリアシングの無い、１つのｘ−ｆ空間データを得る。また、生成機能１２２ｂは、ｘ−ｆ空間データに対して時間軸方向に逆フーリエ変換を適用し、ｘ−ｔ空間データ（水素分布画像）を生成する。なお、以下では、チャネル数分のｘ−ｆ空間データを、「アンフォールド前のｘ−ｆ空間データ」と呼び、アンフォールド後に得られた１つのｘ−ｆ空間データを、「アンフォールド後のｘ−ｆ空間データ」と呼び、区別する場合がある。
【００４８】
次に、ＳＥＮＳＥにおけるアンフォールド処理を説明する。以下では、アンフォールドされた画像で対応する信号点の位置をｐ＝１,…,ｎｐで表し、ｋ空間をｎｐ信号点間隔でサンプリングすることを考える。このとき、各コイルで得られたｋ空間データから、フーリエ変換（あるいは逆フーリエ変換）によって生成される画像には、折り返し（ｆｏｌｄ）信号が重ねられる。ＳＥＮＳＥでは、各コイルの感度が既知の場合に、得られた折り返し画像から、折り返されていない（アンフォールド）画像ｚ（ｐ）を復元することを考える。なお、以下の説明では、フーリエ変換、逆フーリエ変換いずれについても、フーリエ変換と呼ぶ。
【００４９】
ｎｐ信号点間隔でサンプリングしたｋ空間データにフーリエ変換を適用すると、ｎｐ個の信号点を重ね合わせた折り返し画像が得られる。ＳＥＮＳＥでは、折り返し画像内で特定の１信号点が生成された過程に着目する。折り返し画像内で注目する１信号点に対し、コイルをｃ＝１,…,ｎｃ、各コイルでの折り返し画像の信号点値をρ（ｃ）、アンフォールドされた画像で対応する信号点値をｚ（ｐ）、アンフォールドされた画像に対する各コイルの感度をＳ（ｃ，ｐ）で表すと、各コイルに対して、次の（１）式が成り立つ。
【数１】
0006548531.tif 000002
【００５０】
この関係式から、コイルの感度が既知であるなら、未知の変数ｚ（ｐ）がｎｐ個である式がコイル数ｎｃ本得られる。この関係式をコイルの本数並べた線形システムを、行列表記を用いてρ＝Ｓzと表す。このとき、コイルのノイズ共分散行列Ψを考慮すると、ｚは、以下のように推定することができる。
【数２】
0006548531.tif 000003
【００５１】
以下の実施形態では、ρからｚを推定する処理を「アンフォールド」と呼ぶ。ここで、ｎｐ信号点間隔のサンプリングにより収集されたｋ空間データに対して、フーリエ変換及びアンフォールドを行うことで得られた画像のＳＮＲ（Signal  to  Noise  Ratio）は、間引きを行わないフルサンプリングにより収集されたｋ空間データから生成された画像のＳＮＲ、信号点間隔ｎｐ、及びｇ−ｆａｃｔｏｒ（ｇｐ）を用いて、次の（３）式で表せることが知られている。
【数３】
0006548531.tif 000004
【００５２】
なお、ＳＮＲは、信号の分散を雑音の分散で除した値である。ＳＮＲの値が大きいほどノイズが少なく、ＳＮＲの値が小さいほどノイズが大きい。また、ｇ−ｆａｃｔｏｒは、以下の（４）式を満たす、１以上の数値である。
【数４】
0006548531.tif 000005
【００５３】
この場合、信号点間隔ｎｐは、フルサンプリングデータに対するサンプリングデータの間引き率Ｒ（Reduction  factor）に等しくなる。
【００５４】
時系列のｋ空間データからＭＲ画像を再構成する場合には、ｋ空間データに対して空間方向にフーリエ変換を適用するだけでなく時間軸方向にもフーリエ変換を適用してｘ−ｆ空間データを算出し、ｘ−ｆ空間上でアンフォールドを行うｋ−ｔ  ＳＥＮＳＥを利用することで、ＳＥＮＳＥに比較して高い画質で画像を再構成できることが知られている。なお、ｘ−ｆ空間のｘは、空間方向を表し、ｆは、変換後の空間の時間変換軸方向を表す。
【００５５】
続いて、第１の実施形態における処理手順を説明する。
【００５６】
図２は、第１の実施形態における処理手順を示すフローチャートである。まず、収集機能１２６ａが、操作者から撮像条件の入力を受け付ける（ステップＳ１０１）。続いて、収集機能１２６ａが、操作者から入力された撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御回路１１０に送信することで、準備スキャンの実行を制御する（ステップＳ１０２）。ここで、準備スキャンには、例えば、位置決め用の画像を収集するスキャンや、静磁場の不均一性を補正するシミングスキャン、感度分布情報を収集するスキャン等が含まれる。
【００５７】
準備スキャンの終了後、続いて、収集機能１２６ａは、操作者から入力された撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御回路１１０に送信することで、出力画像（例えば、診断用に出力される画像）を収集するイメージングスキャンの実行を制御する（ステップＳ１０３）。第１の実施形態において、収集機能１２６ａは、イメージングスキャンにおいて、単一又は複数チャネルについて、時系列のｋ空間データを収集する。
【００５８】
第１の実施形態においては、シーケンス制御回路１１０は、Ｒ（Reduction  factor；フルサンプリングデータに対するサンプリングデータの間引き率）に従って、データを間引きながら時系列ｋ空間データを収集する。すなわち、シーケンス制御回路１１０が収集する時系列ｋ空間データの個数は、フルサンプリングに比べて１／Ｒとなる。具体的には、シーケンス制御回路１１０は、ｋ空間データとして、「ＲＯ（Read  Out）方向の信号点の個数×ＰＥ（Phase  Encode）方向の信号点の個数×時間方向のフレーム数÷Ｒ×チャネル数」で求められる個数の信号点を収集する。図３は、第１の実施形態に係る収集位置の一例を示す図である。シーケンス制御回路１１０は、例えば、図３に示すように、ＰＥ方向の収集位置をずらしながら時系列ｋ空間データを収集する。図３において、黒丸は、サンプリングが行われる位置を表し、この例ではＲ＝３となる。
【００５９】
そして、画像生成回路１２２は、ステップＳ１０３において収集され記憶部１２３に格納された複数チャネル分の時系列のｋ空間データを用いて、ＭＲ画像の生成を行う（ステップＳ１０４）。生成されたＭＲ画像は、必要に応じて、記憶部１２３に格納されたり、表示部１２５に表示されたりする。ここで、ステップＳ１０１〜Ｓ１０３は、処理回路１２６が記憶回路１２３から収集機能１２６ａに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【００６０】
図４は、第１の実施形態におけるステップＳ１０４の処理手順を示すフローチャートである。
【００６１】
ステップＳ２０１でデータ付加機能１２２ａが実行する処理について説明する。まず、データ付加機能１２２ａは、収集機能１２６ａによって収集された複数のｋ空間データである複数の収集ｋ空間データに基づいて、収集ｋ空間データに付加するデータである付加ｋ空間データを生成する。そして、データ付加機能１２２ａは、複数の収集ｋ空間データに複数の付加ｋ空間データを付加することにより、複数の付加ｋ空間データが付加された複数の収集ｋ空間データを、再構成に用いられる複数のｋ空間データである複数の再構成用ｋ空間データとして生成する。そして、データ付加機能１２２ａは、生成した複数の再構成用ｋ空間データを、生成機能１２２ｂに受け渡す。ステップＳ２０１は、画像生成回路１２２が記憶回路１２３からデータ付加機能１２２ａに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【００６２】
例えば、複数の収集ｋ空間データは、ｋ空間のＲＯ方向、ｋ空間のＰＥ方向、及び、時間がそれぞれ軸となる時系列データである。収集ｋ空間データは、チャネルの個数分だけ存在する。複数の付加ｋ空間データは、時系列データである複数の収集ｋ空間データの前後の時間に該当する複数のｋ空間データである。図５は、第１の実施形態に係る再構成用ｋ空間データの一例を示す図である。以降の説明では、図５の例に示すように、データ付加機能１２２ａが、複数の収集ｋ空間データの時間軸方向の後ろのみに複数の付加ｋ空間データを付加する場合について説明するが、前のみに複数の付加ｋ空間データを付加してもよいし、前後に複数の付加ｋ空間データを付加してもよい。また、複数の付加ｋ空間データの生成フレーム数はＲフレーム以上のＲの倍数とする。付加ｋ空間データは収集ｋ空間データと同一の間引き方でサンプリングされたデータである。
【００６３】
データ付加機能１２２ａは、時間軸方向に、ある範囲内に並んでいる複数の収集ｋ空間データに対して、複数の収集ｋ空間データが並んでいる範囲と時間軸方向において隣接するある範囲内に並んだ複数の付加ｋ空間データを、複数の収集ｋ空間データ及び複数の付加ｋ空間データからなる複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有するように付加する。すなわち、データ付加機能１２２ａは、複数の収集ｋ空間データが並んでいる範囲と複数の収集ｋ空間データが並んでいる範囲とが並ぶ方向に並べられた場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が、連続性を有するように、複数の収集ｋ空間データに対して複数の付加ｋ空間データを付加する。
【００６４】
なお、収集ｋ空間データは、特許請求の範囲に記載された第１のデータの一例である。また、付加ｋ空間データは、特許請求の範囲に記載された第２のデータの一例である。また、再構成用ｋ空間データは、特許請求の範囲に記載された第３のデータの一例である。また、複数の収集ｋ空間データが並んでいる範囲は、特許請求の範囲に記載された第１の範囲の一例である。また、複数の付加ｋ空間データが並んでいる範囲は、特許請求の範囲に記載された第２の範囲の一例である。
【００６５】
例えば、図５の例において、データ付加機能１２２ａは、時間軸方向において時刻０〜１１の範囲内に並んだ複数の収集ｋ空間データに対して、時間軸方向において時刻０〜１１の範囲に隣接する時刻１２〜１４の範囲内に並んだ複数の付加ｋ空間データを、時刻０〜１４の範囲の複数の再構成用ｋ空間データが時間軸方向において周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が連続性を有するように付加する。なお、この例では、時刻０，１，・・・，１４のそれぞれの時刻が一つのフレームに対応する。すなわち、各時刻においてサンプリングされたｋ空間データに基づいてＭＲ画像が生成される。
【００６６】
図６は、第１の実施形態に係るステップＳ２０１の詳細な処理手順を示すフローチャートである。付加ｋ空間データの各フレームについてステップＳ３０１からステップＳ３０３が実行される。
【００６７】
データ付加機能１２２ａは、まず、目標フレームを設定する（ステップＳ３０１）。データ付加機能１２２ａは、上述した関数が上述した連続性を有するように、複数の付加ｋ空間データのうち時間的に最も早い付加ｋ空間データのフレーム（時間的に最初のフレーム）に対する目標フレームを、連続性を考慮したときに時間的に最も近い、すなわち、複数の収集ｋ空間データのうち最も遅い収集ｋ空間データのフレーム（時間的に最後のフレーム）に対応するＭＲ画像の信号強度とする。例えば、図５の例において、データ付加機能１２２ａは、複数の付加ｋ空間データのうち時間的に最も早い付加ｋ空間データのフレームである時刻１２に対する目標フレームとして、時刻１１の収集ｋ空間データのフレームの信号強度を設定する。
【００６８】
また、データ付加機能１２２ａは、複数の付加ｋ空間データのうち時間的に最も遅い付加ｋ空間データのフレームに対する目標フレームを、収集ｋ空間データのうち時間的に最も早い収集ｋ空間データのフレームの信号強度とする。例えば、図５の例において、データ付加機能１２２ａは、複数の付加ｋ空間データのうち時間的に最後の付加ｋ空間データのフレームである時刻１４に対する目標フレームとして、時刻０の収集ｋ空間データのフレームの信号強度を設定する。
【００６９】
また、データ付加機能１２２ａは、複数の付加ｋ空間データのうち、時間的に最初と最後との間に位置する付加ｋ空間データのフレームの目標フレームとして、複数の付加ｋ空間データのうち時間的に最初の付加ｋ空間データのフレームの目標フレームに対応するＭＲ画像の信号強度と最後の付加ｋ空間データのフレームの目標フレームに対応するＭＲ画像の信号強度との所定の重み付け加算により算出した信号強度を設定する。例えば、データ付加機能１２２ａは、複数の付加ｋ空間データのうち時間的に最初の付加ｋ空間データと最後の付加ｋ空間データとの間に位置する付加ｋ空間データのフレームと、最初の付加ｋ空間データのフレーム及び最後の付加ｋ空間データのフレームのそれぞれとの時間的な距離に応じて重みを算出する。例えば、データ付加機能１２２ａは、距離が長くなるほど、重みを小さくし、距離が短くなるほど、重みを大きくする。そして、データ付加機能１２２ａは、算出した重みを用いた上述の所定の重み付け加算により、複数の付加ｋ空間データのうち時間的に最初の付加ｋ空間データと最後の付加ｋ空間データとの間に位置する付加ｋ空間データのフレームに対応するＭＲ画像の信号強度を算出し、算出した信号強度を目標フレームとして設定する。
【００７０】
図５の例において、データ付加機能１２２ａは、時刻１２の付加ｋ空間データのフレームの目標フレームである時刻１１の収集ｋ空間データのフレームの信号強度と、時刻１４の付加ｋ空間データのフレームの目標フレームである時刻０の収集ｋ空間データのフレームの信号強度との所定の重み付け加算により、時刻１３に対する目標フレームを算出する。
【００７１】
続いて、データ付加機能１２２ａは、複数の収集ｋ空間データの複数のフレームの中から、ステップＳ３０１で設定した目標フレームに類似する信号強度をもつフレームを類似フレームとして選択する（ステップＳ３０２）。
【００７２】
図５の例に示すように、複数の付加ｋ空間データの各フレームは、時刻に応じてサンプリングパターン（間引きパターンとも称される）が一意に定まる。そこで、データ付加機能１２２ａは、複数の付加ｋ空間データの各時刻に対応するサンプリングパターンで収集された複数の収集ｋ空間データの複数のフレームの中から、各時刻の目標フレームと類似するフレームを探索する。この際、異なるサンプリングパターンで間引きサンプリングされたｋ空間データ間の類似度を算出することは困難であるため、データ付加機能１２２ａは、各フレーム（目標フレームも含む）の全チャネルデータを逆フーリエ変換してｘ空間に変換し、変換後の係数値の絶対値で類似度を算出する。データ付加機能１２２ａは、間引きデータを逆フーリエ変換した場合、ｘ空間データには折り返しが発生するため、折り返しが発生しているｘ空間データで類似度を算出する。例えば、データ付加機能１２２ａは、ＳＡＤ（Sum  of  Absolute  Difference）やＳＳＤ（Sum  of  Squared  Difference）等の類似度の算出手法を用いて、類似度を算出する。そして、データ付加機能１２２ａは、例えば、類似度が最大（ＳＡＤやＳＳＤが最小）となったフレームを類似フレームとして選択する。
【００７３】
続いて、データ付加機能１２２ａは、ステップＳ３０２で選択した類似フレームの収集ｋ空間データを付加ｋ空間データのフレームにコピーすることにより、収集ｋ空間データに付加ｋ空間データを付加する（ステップＳ３０３）。
【００７４】
データ付加機能１２２ａは、全ての付加ｋ空間データのフレームのデータ生成が完了したか否かを判定する（ステップＳ３０４）。完了していない場合（ステップＳ３０４；Ｎｏ）には、データ付加機能１２２ａは、ステップＳ３０１に戻る。一方、完了した場合（ステップＳ３０４；Ｙｅｓ）には、データ付加機能１２２ａは、複数の付加ｋ空間データが付加された複数の収集ｋ空間データを複数の再構成用ｋ空間データとして生成機能１２２ｂに受け渡し、Ｓ２０２に進む。ここで、ステップＳ３０１〜Ｓ３０４は、画像生成回路１２２が記憶回路１２３からデータ付加機能１２２ａに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【００７５】
図７は、第１の実施形態に係る再構成用ｋ空間データに対応するＭＲ画像の信号強度の一例を示す図である。図７の例に示すように、ステップＳ３０１〜Ｓ３０３での処理によって、複数の収集ｋ空間データに複数の付加ｋ空間データを付加した結果、時間軸方向において複数の再構成用ｋ空間データに対応する複数のＭＲ画像の両端の信号強度の差が小さくなっている。これにより、複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が、時間軸方向において連続性を有するようになる。
【００７６】
ステップＳ３０１〜Ｓ３０３において、データ付加機能１２２ａは、生成対象である複数の付加ｋ空間データの複数のフレームのうち、複数の収集ｋ空間データが並んでいる範囲に最も近いフレーム（以下、近接フレームと称する）について、複数の収集ｋ空間データの複数のフレームのうち、近接フレームに最も近いフレーム（より具体的には、フレームに対応するＭＲ画像の信号強度）を目標フレームとして設定する。そして、データ付加機能１２２ａは、近接フレームに対応するサンプリングパターンで収集された複数の収集ｋ空間データの複数のフレームの中から、対応するＭＲ画像の信号強度について、目標フレームと類似するフレームを、近接フレームとして生成する。なお、収集ｋ空間データのフレームは、特許請求の範囲に記載された第１のフレームの一例である。また、付加ｋ空間データのフレームは、特許請求の範囲に記載された第２のフレームの一例である。
【００７７】
また、ステップＳ３０１〜Ｓ３０３において、データ付加機能１２２ａは、生成対象である複数の付加ｋ空間データの複数のフレームのうち、複数の収集ｋ空間データが並んでいる範囲から最も遠いフレーム（以下、遠方フレームと称する）について、複数の収集ｋ空間データの複数のフレームのうち、遠方フレームから最も遠いフレーム（より具体的には、フレームに対応するＭＲ画像の信号強度）を目標フレームとして設定する。そして、データ付加機能１２２ａは、遠方フレームに対応するサンプリングパターンで収集された複数の収集ｋ空間データの複数のフレームの中から、対応するＭＲ画像の信号強度について、目標フレームと類似するフレームを、遠方フレームとして生成する。
【００７８】
また、ステップＳ３０１〜Ｓ３０３において、データ付加機能１２２ａは、複数の付加ｋ空間データの複数のフレームのうち、近接フレーム及び遠方フレーム以外のフレームについて、近接フレームの目標フレームに対応するＭＲ画像の信号強度と、遠方フレームの目標フレームに対応するＭＲ画像の信号強度とに基づいて信号強度を算出することにより、近接フレーム及び遠方フレーム以外のフレームを生成する。
【００７９】
また、ステップＳ３０２において、データ付加機能１２２ａは、目標フレームに対してフーリエ変換を適用することにより得られた画像空間（目標画像空間とも称される）と、近接フレーム、又は、遠方フレームに対応するサンプリングパターンで収集された複数の収集ｋ空間データの複数のフレームのそれぞれに対してフーリエ変換を適用することにより得られた画像空間（以下、候補画像空間と称する）のそれぞれとの類似度を算出する。そして、データ付加機能１２２ａは、算出した類似度が示す類似性が最も高い候補画像空間に対応する収集ｋ空間データのフレームを、近接フレーム、又は、遠方フレームとして生成する。なお、目標画像空間は、特許請求の範囲に記載された第１の画像空間の一例である。また、候補画像空間は、特許請求の範囲に記載された第２の画像空間の一例である。
【００８０】
図４の説明に戻る。ステップＳ２０２で生成機能１２２ｂが実行する処理について説明する。生成機能１２２ｂは、複数の再構成用ｋ空間データに対して再構成処理を行い、複数のＭＲ画像を生成する。例えば、生成機能１２２ｂは、データ付加機能１２２ａから複数の再構成用ｋ空間データが受け渡されると、複数の再構成用ｋ空間データを時空間方向に逆フーリエ変換してアンフォールド前のｘ−ｆ空間データを取得する。そして、生成機能１２２ｂは、アンフォールド前のｘ−ｆ空間データをアンフォールドし、アンフォールド後のｘ−ｆ空間データを生成する。そして、生成機能１２２ｂは、アンフォールド後のｘ−ｆ空間データに対して、時間軸方向にフーリエ変換を適用し、ｘ−ｔ空間データを取得する。そして、生成機能１２２ｂは、ｘ−ｔ空間データを削除機能１２２ｃに受け渡す。このｘ−ｔ空間データは、時系列に沿った複数のＭＲ画像を示すデータである。
【００８１】
なお、フーリエ変換及び逆フーリエ変換は、変換対象のデータが周期的に並んでいると仮定した上で、変換対象のデータに対して変換を行う処理である。すなわち、ステップＳ２０２では、生成機能１２２ｂは、変換対象のデータが周期的に並んでいると仮定した上で変換対象のデータに対して変換を行う処理により、再構成用ｋ空間データから、時系列のＭＲ画像を生成する。ここで、再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応するＭＲ画像の信号強度を表す関数は、連続性を有するので、ステップＳ２０２では、生成したＭＲ画像において、逆フーリエ変換特有のアーティファクトの発生を抑制することができる。
【００８２】
削除機能１２２ｃは、ステップＳ２０１で複数の収集ｋ空間データに付加された複数の付加ｋ空間データの時刻に該当するデータをｘ−ｔ空間データから削除し、ｘ−ｔ空間データのフレーム数を複数の収集ｋ空間データのフレーム数と同一にし（ステップＳ２０３）、処理を終了する。例えば、図５の例において、削除機能１２２ｃは、複数の付加ｋ空間データの時刻１２〜１４に該当するデータをｘ−ｔ空間データから削除し、ｘ−ｔ空間データのフレーム数を複数の収集ｋ空間データのフレーム数と同一にする。
【００８３】
すなわち、ステップＳ２０３では、削除機能１２２ｃは、生成機能１２２ｂにより生成されたｘ−ｔ空間データが示す複数のＭＲ画像から、付加ｋ空間データが並んでいる範囲に対応するＭＲ画像を削除する。これにより、収集ｋ空間データのみから生成されるフレーム数のＭＲ画像を得ることができる。
【００８４】
第１の実施形態によれば、複数の収集ｋ空間データが周期的に並んでいる場合の収集ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が連続性を有しない場合にも、複数の付加ｋ空間データを付加して、複数の再構成用ｋ空間データが周期的に並んでいる場合の複数の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が連続性を有するようにする。したがって、第１の実施形態によれば、変換対象のデータが周期的に並んでいると仮定した上で変換を行う処理、例えば、時空間再構成法を実行した際に発生し得る特有のアーティファクトの発生を抑制することができる。また、この結果、再構成精度の低下を抑制することができる。
【００８５】
（第１の実施形態の変形例１：類似フレーム探索領域の限定）
なお、データ付加機能１２２ａは、ステップＳ３０２で類似フレームを選択する際に、複数の収集ｋ空間データの全フレームを探索対象とするのではなく、探索対象を一部のフレームに限定してもよい。例えば、データ付加機能１２２ａは、近接フレームに類似する類似フレームを選択する際には、探索対象を、全フレーム中、時間軸方向において後ろの１／４のフレームにする。また、データ付加機能１２２ａは、遠方フレームに類似する類似フレームを選択する際には、探索対象を、全フレーム中、時間軸方向において前の１／４のフレームにする。また、データ付加機能１２２ａは、近接フレームと遠方フレームとの間のフレームに類似する類似フレームを選択する際には、探索対象を、全フレーム中、時間軸方向において後ろの１／４のフレームと前の１／４のフレームを除いたフレームにする。これにより、探索時間を短縮することができる。
【００８６】
（第１の実施形態の変形例２：類似度をｋ空間データで算出）
また、データ付加機能１２２ａは、ステップ３０２で類似度を算出する際に、折り返しを含むｘ空間で類似度を算出するのではなく、間引きｋ空間データを補間して生成した間引きのないｋ空間データで類似度を算出してもよい。
【００８７】
（第１の実施形態の変形例３：探索処理を行わない）
また、データ付加機能１２２ａは、ステップＳ３０１で目標フレームを設定する際に、付加ｋ空間データのフレームのサンプリングパターンと同一のサンプリングパターンを有する複数の収集ｋ空間データの中から目標フレームに対して最も時刻が近いフレームを選択してもよい。この場合、データ付加機能１２２ａは、ステップＳ３０２では類似フレームの選択を行わず、ステップＳ３０３では、ステップＳ３０１で設定した目標フレームのｋ空間データをコピーする。これにより、探索処理を行わないため処理時間を短縮することができる。
【００８８】
また、データ付加機能１２２ａは、ステップＳ３０２で、１枚の類似フレームを選択するのではなく、付加ｋ空間データのフレームのサンプリングパターンと同一のサンプリングパターンで収集された複数の収集ｋ空間データのフレーム（以下、同一パターンフレームと称する）の重み付き加算により付加ｋ空間データを生成してもよい。例えば、データ付加機能１２２ａは、目標フレームと、全ての同一パターンフレームのそれぞれから、対応するＭＲ画像の信号強度を表す特徴量を算出し、目標フレームの特徴量との差の絶対値を用いて重みを算出する。例えば、データ付加機能１２２ａは、目標フレームの特徴量との差が小さければ重みを大きく、差が大きければ重みを小さくする。信号強度を表す特徴量は、例えば、中心周波数付近のｋ空間データの加算値、又は、逆フーリエ変換して得られるｘ空間データの信号の加算値でもよい。また、データ付加機能１２２ａは、目標フレームから算出した特徴量が、同一パターンフレームから算出した特徴量の最大値よりも大きい、又は最小値よりも小さい場合には、最大値、又は最小値をとる同一パターンフレームをコピーしてもよい。これにより、探索処理を行わないため処理時間を短縮することができる。
【００８９】
（第２の実施形態）
次に、第２の実施形態について説明する。第２の実施形態では、処理回路と画像生成回路の処理内容、及び、付加データを処理回路で収集する点が第１の実施形態と異なる。第２の実施形態において、第１の実施形態と同様の構成については同一の符号を付し、説明を省略する場合がある。
【００９０】
第１の実施形態では、データ付加機能１２２ａが付加ｋ空間データを生成するが、第２の実施形態では、データ付加機能２２６ｂが付加ｋ空間データを収集する。
【００９１】
図８は、第２の実施形態に係るＭＲＩ装置２００の構成を示す機能ブロック図である。図８の例に示すように、ＭＲＩ装置２００は、図１に示すＭＲＩ装置１００と比較して、データ付加機能２２６ｂを新たに追加した処理回路２２６を備え、データ付加機能１２２ａを除いた画像生成回路２２２を保持する点が異なる。
【００９２】
処理回路２２６は、収集機能１２６ａで空間方向のサンプリング位置を変化させながら、複数のチャネルについて、時系列のｋ空間データを収集し、データ付加機能２２６ｂで収集機能１２６ａと同一のサンプリング方法で時系列のｋ空間データを収集する。
【００９３】
画像生成回路２２２は、処理回路２２６で収集されたデータや、処理回路２２６で収集され記憶部１２３で記憶されたデータを用いて、画像の生成を行う。生成機能１２２ｂの処理内容は第１の実施形態と同じである。削除機能１２２ｃは、データ付加機能２２６ｂで複数の収集ｋ空間データに付加された複数の付加ｋ空間データに対応する時間に該当するアンフォールドデータを削除する。
【００９４】
図９は、第２の実施形態における処理手順を示すフローチャートである。ステップＳ１０１〜Ｓ１０３は第１の実施形態と同一の内容なので説明を省略する。
【００９５】
データ付加機能２２６ｂは、複数の付加ｋ空間データを収集し、ステップＳ１０３で収集した複数の収集ｋ空間データに複数の付加ｋ空間データを付加して、複数の再構成用ｋ空間データを得る（ステップＳ４０４）。なお、ステップＳ４０４では、データ付加機能２２６ｂは、更に、得られた複数の再構成用ｋ空間データを、画像生成回路２２２に受け渡すとともに記憶回路１２３に格納する。ここで、ステップＳ４０４は、画像生成回路２２２が記憶回路１２３からデータ付加機能１２２ａに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【００９６】
例えば、データ付加機能２２６ｂが収集する複数の付加ｋ空間データは、例えば、ｋ空間のＲＯ方向、ｋ空間のＰＥ方向、及び、時間がそれぞれ軸となる時系列データである。付加ｋ空間データは、チャネルの個数分だけ存在し、収集ｋ空間データと同一の間引き方でサンプリングされたデータである。データ付加機能２２６ｂは、収集機能１２６ａと同時に実行されてもよいし、前後の別の時間に実行されてもよい。データ付加機能２２６ｂは、複数の付加ｋ空間データを、複数の収集ｋ空間データの時間軸で前、又は後ろ、又は前後に付加する。複数の付加ｋ空間データの収集フレーム数はＲフレーム以上のＲの倍数とする。データ付加機能２２６ｂは、信号変化の様子を任意のタイミングで収集可能な場合には、複数の収集ｋ空間データの時間軸で後ろのフレームと同等のタイミングで付加ｋ空間データの時間軸で前半のデータを収集し、複数の収集ｋ空間データの時間軸で前のフレームと同等のタイミングで付加ｋ空間データの後半のデータを収集する。このようにして、複数の付加ｋ空間データを収集することにより、複数の収集ｋ空間データに複数の付加ｋ空間データを付加した複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応するＭＲ画像の信号強度を表す関数は、先の図７の例に示すように、時間軸方向に連続性を有する。
【００９７】
ここで、任意のタイミングで信号収集が可能なシーケンスとしては、例えば、反転回復パルスシーケンスにおけるＴＩ（inversion  time）を複数回収集し、且つ各ＴＩを任意に設定可能なシーケンス等が挙げられる。
【００９８】
生成機能１２２ｂは、複数の再構成用ｋ空間データから複数のＭＲ画像を生成し（ステップＳ４０５）、処理を終了する。ここで、ステップＳ４０５は、画像生成回路２２２が記憶回路１２３から生成機能１２２ｂに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【００９９】
図１０は、第２の実施形態におけるステップＳ４０５の処理手順を示すフローチャートである。ステップＳ５０２で生成機能１２２ｂが実行する処理について説明する。生成機能１２２ｂは、処理回路２２６から受け渡された複数の再構成用ｋ空間データを時空間方向に逆フーリエ変換してアンフォールド前のｘ−ｆ空間データを取得する。そして、生成機能１２２ｂは、アンフォールド前のｘ−ｆ空間データをアンフォールドし、アンフォールド後のｘ−ｆ空間データを生成する。そして、生成機能１２２ｂは、アンフォールド後のｘ−ｆ空間データに対して、時間軸方向にフーリエ変換を適用し、ｘ−ｔ空間データを得る。そして、生成機能１２２ｂは、ｘ−ｔ空間データを削除機能１２２ｃに受け渡す。このｘ−ｔ空間データは、時系列に沿った複数のＭＲ画像を示すデータである。
【０１００】
削除機能１２２ｃは、受け渡されたｘ−ｔ空間データの一部のデータを削除する。例えば、削除機能１２２ｃは、ステップＳ４０４で複数の収集ｋ空間データに付加された複数の付加ｋ空間データの時刻に該当するデータをｘ−ｔ空間データから削除し、収集ｋ空間データと同一のフレーム数にする（ステップＳ５０３）。ここで、ステップＳ５０２は、画像生成回路２２２が記憶回路１２３から生成機能１２２ｂに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。また、ステップＳ５０３は、画像生成回路２２２が記憶回路１２３から削除機能１２２ｃに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【０１０１】
以上、第２の実施形態について説明した。第２の実施形態に係るＭＲＩ装置２００によれば、第１の実施形態に係るＭＲＩ装置１００と同様に、変換対象のデータが周期的に並んでいると仮定した上で変換を行う処理を行う際に発生し得る特有のアーティファクトの発生を抑制することができる。
【０１０２】
（第２の実施形態の変形例：リアルタイム収集）
なお、第２の実施形態において、ｋ空間データの収集の態様が、信号変化の様子を逐次収集するリアルタイム収集の場合には、任意のタイミングでの収集ができない。この場合には、ステップＳ４０４で、データ付加機能２２６ｂは、複数の収集ｋ空間データの前後の時刻に該当する信号を複数の付加ｋ空間データとして収集する。これにより、複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数は時間軸方向に連続性を有しない可能性があるが、再構成精度が低下するフレームを収集ｋ空間データからずらすことができる。
【０１０３】
（第３の実施形態）
次に、第３の実施形態について説明する。第３の実施形態では、画像生成回路の処理内容、及び、再構成用ｋ空間データに窓関数を適用する点が第１の実施形態及び第２の実施形態と異なる。第３の実施形態において、第１の実施形態及び第２の実施形態と同様の構成については同一の符号を付し、説明を省略する場合がある。
【０１０４】
図１１は、第３の実施形態に係るＭＲＩ装置の構成を示す機能ブロック図である。図８に示すＭＲＩ装置２００と比較し、窓関数適用機能３２２ａを新たに備える点が異なる。
【０１０５】
窓関数適用機能３２２ａは、複数の再構成用ｋ空間データに対して窓関数を適用する。この窓関数は、例えば、複数の再構成用ｋ空間データに対応する複数のＭＲ画像のうち時間軸方向における両端のＭＲ画像の信号強度の差を小さくするフィルタである。このようなフィルタを複数の再構成用ｋ空間データに対して適用することにより、複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が連続性を有することとなる。なお、窓関数適用機能３２２ａは、特許請求の範囲に記載された窓関数適用部の一例である。
【０１０６】
図１２は、第３の実施形態における処理手順を示すフローチャートである。第３の実施形態では、第２の実施形態（図９）と比較して、ステップＳ４０５に替えてステップＳ６０５の処理を実行する点が第２の実施形態と異なる。
【０１０７】
図１３は、ステップＳ６０５の処理の流れの一例を示すフローチャートである。図１３の例に示すように、窓関数適用機能３２２ａは、複数の再構成用ｋ空間データの時間軸に対して窓関数を適用する。窓関数適用機能３２２ａは、複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が、時間軸方向において連続性を有するようになる窓関数を適用する（ステップＳ７０１）。例えば、窓関数適用機能３２２ａは、矩形窓などで複数の付加ｋ空間データに対応する時間帯域の信号を下げる窓関数を適用すればよい。ステップＳ７０１は、画像生成回路３２２が記憶回路１２３から窓関数適用機能３２２ａに対応するプログラムを読み出し、読み出したプログラムを実行することにより実現されるステップである。
【０１０８】
なお、複数の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度の変化は、実行するプロトコルの種類に応じて変わると考えられる。そこで、各プロトコルにおける信号変化に対応させた窓関数をプロトコルの種類毎に記憶回路１２３に予め記憶させておき、窓関数適用機能３２２ａが、ステップＳ６０５において、実行中のプロトコルの種類に対応する窓関数を取得し、取得した窓関数を複数の再構成用ｋ空間データの時間軸に対して適用してもよい。
【０１０９】
第３の実施形態によれば、データ付加機能２２６ｂによっても、所望の連続性が得られない場合であっても、複数の再構成用ｋ空間データが周期的に並んでいる場合の再構成用ｋ空間データに対応する複数のＭＲ画像の信号強度を表す関数が、時間軸方向の連続性を有するようにさせることができる。よって、第３の実施形態によれば、第１の実施形態及び第２の実施形態と同様に、変換対象のデータが周期的に並んでいると仮定した上で変換を行う処理を行う際に発生し得る特有のアーティファクトの発生を抑制することができる。
【０１１０】
（第３の実施形態の変形例１：画像空間に窓を適用する）
窓関数適用機能３２２ａは、再構成用ｋ空間データの各チャネルのｋ−ｔ空間データに対して空間方向にフーリエ変換を適用し、各チャネルのｘ−ｔ空間データを得たのちに、各チャネルのｘ−ｔ空間データに対して窓関数を適用してもよい。この場合、生成機能１２２ｂは、時間軸方向のみにフーリエ変換を適用することで、各チャネルのｘ−ｆ空間データを得る。
【０１１１】
（第３の実施形態の変形例２：第１の実施形態の変形）
図１１に示すＭＲＩ装置３００は、第２の実施形態に係るＭＲＩ装置２００の変形として説明したが、第１の実施形態に係るＭＲＩ装置１００の変形でもよい。この際には、窓関数適用部は、データ付加機能１２２ａで複数の付加ｋ空間データを生成し、生成した複数の付加ｋ空間データを複数の収集ｋ空間データに付加した後の複数の再構成用ｋ空間データに対して窓関数を適用する。この際、変形例１と同様に画像空間に変換してから窓関数を適用してもよい。
【０１１２】
なお、上述の各実施形態では、収集機能１２６ａにより収集されるデータが、時系列のｋ空間データ、すなわち、ｋ−ｔ空間データである場合について説明したが、収集されるデータはこれに限られない。例えば、ＭＲＩ装置は、ｋ空間データを任意の順序で収集してもよい。
【０１１３】
例えば、ＭＲＩ装置は、ｋ空間データをスライス方向の軸に沿った空間的な順序で収集してもよい。なお、この場合には、生成機能１２２ｂは、ＣＡＩＰＩ（Controlled  Aliasing  In  Parallel  Imaging  Results  IN  Higher  Acceleration）というＭＲ画像の再構成の手法を用いて、ＭＲ画像の再構成を行う。
【０１１４】
また、例えば、ＭＲＩ装置は、ｋ空間データを時間ｔではなくパラメータｐの軸に沿った順序で収集してもよい。ここでいうｐとは、任意の軸を指す。なお、ｐには、例えば、ｂ値などが含まれる。ｂ値とは、拡散強調画像の撮像において被検体に印加するＭＰＧ（Motion  Probing  Gradient）パルスの強さを表す。拡散強調画像の撮像では、ｂ値を変えて複数回ＭＰＧパルスが印加される。
【０１１５】
（第４の実施形態）
ここで、第１の実施形態に係るＭＲＩ装置１００の機能、第２の実施形態に係るＭＲＩ装置２００の機能、又は、第３の実施形態に係るＭＲＩ装置３００の機能を、ＭＲＩ装置とネットワークを介して接続された画像処理装置に持たせることができる。このような実施形態を第４の実施形態として、図１４を用いて説明する。
【０１１６】
図１４は、第４の実施形態に係る画像処理装置の構成を示す機能ブロック図である。図１４の例に示すシステムは、ＭＲＩ装置６００と、画像保管装置７００と、画像表示装置８００と、画像処理装置９００とを有する。ＭＲＩ装置６００と、画像保管装置７００と、画像表示装置８００と、画像処理装置９００とは、例えば、病院内に設置された院内ＬＡＮ（Local  Area  Network）５００により、直接的、又は間接的に相互に通信可能な状態となっている。例えば、ＰＡＣＳ（Picture  Archiving  and  Communication  System）が導入されている場合、各装置６００〜９００は、ＤＩＣＯＭ（Digital  Imaging  and  Communications  in  Medicine）規格に則って、画像等を相互に送受信する。
【０１１７】
ＭＲＩ装置６００は、第１の実施形態に係るＭＲＩ装置１００、第２の実施形態に係るＭＲＩ装置２００、又は、第３の実施形態に係るＭＲＩ装置３００である。ＭＲＩ装置６００は、例えば、収集ｋ空間データや再構成用ｋ空間データを画像処理装置９００に送信する。
【０１１８】
画像保管装置７００は、ＭＲＩ装置６００および画像処理装置９００により再構成されたＭＲ画像を保管するデータベースである。
【０１１９】
画像処理装置９００は、ワークステーションであり、第１の実施形態に係るＭＲＩ装置の機能、第２の実施形態に係るＭＲＩ装置の機能、又は、第３の実施形態に係るＭＲＩ装置の機能と同等の機能を有する。画像処理装置９００は、ＭＲＩ装置６００から送信された収集ｋ空間データや再構成用ｋ空間データを用いて、第１の実施形態に係るＭＲＩ装置が実行する処理、第２の実施形態に係るＭＲＩ装置が実行する処理（収集ｋ空間データを収集する処理は除く）、又は、第３の実施形態に係るＭＲＩ装置が実行する処理と同様の処理を行う。
【０１２０】
画像処理装置９００は、インタフェース回路９０１、記憶回路９０２、入力回路９０４、ディスプレイ９０５及び処理回路９０６を有する。
【０１２１】
インタフェース回路９０１は、ＭＲＩ装置６００、画像保管装置７００及び画像表示装置８００と、画像処理装置９００との通信を行うためのインタフェースである。例えば、インタフェース回路９０１は、ＭＲＩ装置６００から送信された収集ｋ空間データや再構成用ｋ空間データを受信すると、受信した収集ｋ空間データや再構成用ｋ空間データを処理回路９０６に送信する。インタフェース回路９０１は、例えば、プロセッサにより実現される。
【０１２２】
記憶回路９０２は、収集ｋ空間データや再構成用ｋ空間データや、後述の生成機能９０６ｂによって生成されたＭＲ画像などを記憶する。また、記憶回路９０２は、各種のプログラムを記憶する。記憶回路９０２は、例えば、ＲＡＭ（Random  Access  Memory）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。
【０１２３】
入力回路９０４は、操作者からの各種指示や情報入力を受け付ける。入力回路９０４は、例えば、トラックボール、スイッチボタン、マウス、キーボード等によって実現される。入力回路９０４は、処理回路９０６に接続されており、操作者から受け取った入力操作を電気信号に変換して処理回路９０６へと出力する。
【０１２４】
ディスプレイ９０５は、処理回路９０６による制御のもと、各種ＧＵＩ（Graphical  User  Interface）や、生成機能９０６ｂによって生成されたＭＲ画像等を表示する。
【０１２５】
処理回路９０６は、画像処理装置９００の全体制御を行う。具体的には、処理回路９０６は、インタフェース回路９０１から送信された収集ｋ空間データや再構成用ｋ空間データを受信すると、受信した収集ｋ空間データや再構成用ｋ空間データを用いて、第１の実施形態に係るＭＲＩ装置が実行する処理、第２の実施形態に係るＭＲＩ装置が実行する処理（収集ｋ空間データを収集する処理は除く）、又は、第３の実施形態に係るＭＲＩ装置が実行する処理と同様の処理を行う。処理回路９０６は、プロセッサにより実現される。処理回路９０６は、データ付加機能９０６ａ、生成機能９０６ｂ及び削除機能９０６ｃを有する。
【０１２６】
データ付加機能９０６ａは、上述したデータ付加機能１２２ａ又はデータ付加機能２２６ｂが発揮する機能（収集ｋ空間データを収集する機能は除く）を有する。生成機能９０６ｂは、上述した生成機能１２２ｂと同様の機能を有する。削除機能９０６ｃは、上述した削除機能１２２ｃと同様の機能を有する。
【０１２７】
ここで、例えば、処理回路９０６の構成要素であるデータ付加機能９０６ａ、生成機能９０６ｂ及び削除機能９０６ｃの各機能は、コンピュータによって実行可能なプログラムの形態で記憶回路９０２に記憶されている。処理回路９０６は、各プログラムを記憶回路９０２から読み出し、読み出した各プログラムを実行することで、各機能を実現する。換言すると、各プログラムを読み出した状態の処理回路９０６は、図１４の処理回路９０６内に示された各機能を有することとなる。
【０１２８】
以上、第４の実施形態に係る画像処理装置９００について説明した。画像処理装置９００によれば、第１の実施形態に係るＭＲＩ装置の機能、第２の実施形態に係るＭＲＩ装置の機能（収集ｋ空間データを収集する機能は除く）、又は、第３の実施形態に係るＭＲＩ装置の機能と同等の機能を有するので、変換対象のデータが周期的に並んでいると仮定した上で変換を行う処理を行う際に発生し得る特有のアーティファクトの発生を抑制することができる。
【０１２９】
（プログラム）
上述した実施形態の中で示した処理手順に示された指示は、ソフトウェアであるプログラムに基づいて実行されることが可能である。汎用の計算機システムが、このプログラムを予め記憶しておき、このプログラムを読み込むことにより、上述した実施形態のＭＲＩ装置や画像処理装置による効果と同様な効果を得ることも可能である。上述した実施形態で記述された指示は、コンピュータに実行させることのできるプログラムとして、磁気ディスク（フレキシブルディスク、ハードディスク等）、光ディスク（ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＯＭ、ＤＶＤ±Ｒ、ＤＶＤ±ＲＷ等）、半導体メモリ、又はこれに類する記録媒体に記録される。コンピュータ又は組み込みシステムが読み取り可能な記憶媒体であれば、その記憶形式は何れの形態であってもよい。コンピュータは、この記録媒体からプログラムを読み込み、このプログラムに基づいてプログラムに記述されている指示をＣＰＵで実行させれば、上述した実施形態のＭＲＩ装置や画像処理装置と同様な動作を実現することができる。もちろん、コンピュータがプログラムを取得する場合又は読み込む場合はネットワークを通じて取得又は読み込んでもよい。
【０１３０】
また、記憶媒体からコンピュータや組み込みシステムにインストールされたプログラムの指示に基づきコンピュータ上で稼働しているＯＳ（オペレーティングシステム）や、データベース管理ソフト、ネットワーク等のＭＷ（ミドルウェア）等が、上述した実施形態を実現するための各処理の一部を実行してもよい。
【０１３１】
更に、記憶媒体は、コンピュータあるいは組み込みシステムと独立した媒体に限らず、ＬＡＮ（Local  Area  Network）やインターネット等により伝達されたプログラムをダウンロードして記憶又は一時記憶した記憶媒体も含まれる。
【０１３２】
また、記憶媒体は１つに限られず、複数の媒体から、上述した実施形態における処理が実行される場合も、実施形態における記憶媒体に含まれ、媒体の構成は何れの構成であってもよい。
【０１３３】
なお、実施形態におけるコンピュータ又は組み込みシステムは、記憶媒体に記憶されたプログラムに基づき、上述した実施形態における各処理を実行するためのものであって、パソコン、マイコン等の１つからなる装置、複数の装置がネットワーク接続されたシステム等の何れの構成であってもよい。
【０１３４】
また、実施形態におけるコンピュータとは、パソコンに限らず、情報処理機器に含まれる演算処理装置、マイコン等も含み、プログラムによって実施形態における機能を実現することが可能な機器、装置を総称している。
【０１３５】
以上述べた少なくとも一つの実施形態の磁気共鳴イメージング装置及び画像処理装置によれば、変換対象のデータが周期的に並んでいると仮定した上で変換を行う処理を行う際に発生し得る特有のアーティファクトの発生を抑制することができる。
【０１３６】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１３７】
１００      ＭＲＩ装置
１２２ａ    データ付加機能
１２２ｂ    生成機能
１２２ｃ    削除機能
１２６ａ    収集機能
【図１】
0006548531.tif 000006
【図２】
0006548531.tif 000007
【図３】
0006548531.tif 000008
【図４】
0006548531.tif 000009
【図５】
0006548531.tif 000010
【図６】
0006548531.tif 000011
【図７】
0006548531.tif 000012
【図８】
0006548531.tif 000013
【図９】
0006548531.tif 000014
【図１０】
0006548531.tif 000015
【図１１】
0006548531.tif 000016
【図１２】
0006548531.tif 000017
【図１３】
0006548531.tif 000018
【図１４】
0006548531.tif 000019