特許 7505872 磁気共鳴イメージング方法及び磁気共鳴イメージング装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-06-17

(45)【発行日】2024-06-25

(54)【発明の名称】磁気共鳴イメージング方法及び磁気共鳴イメージング装置

(51)【国際特許分類】

   A61B 5/055 20060101AFI20240618BHJP

【ＦＩ】

A61B5/055 311

【請求項の数】 11

(21)【出願番号】P 2019184912

(22)【出願日】2019-10-08

(65)【公開番号】P2021058434

(43)【公開日】2021-04-15

【審査請求日】2022-07-25

(73)【特許権者】

【識別番号】594164542

【氏名又は名称】キヤノンメディカルシステムズ株式会社

(74)【代理人】

【識別番号】110001771

【氏名又は名称】弁理士法人虎ノ門知的財産事務所

(72)【発明者】

【氏名】朽名 英明

(72)【発明者】

【氏名】竹島 秀則

【審査官】田中 洋介

(56)【参考文献】

【文献】米国特許出願公開第２０１５／０２７６９０９（ＵＳ，Ａ１）

【文献】特表２０１４－５０２９１０（ＪＰ，Ａ）

【文献】国際公開第２０１９／０７２７１９（ＷＯ，Ａ１）

【文献】特表２０１３－５２１９５５（ＪＰ，Ａ）

【文献】特開２００１－０４６３５４（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１７／０２４２０８６（ＵＳ，Ａ１）

【文献】Keigo Kawaji et al.，Whole Heart Coronary Imaging with Flexible Acquisition Window and Trigger Delay，PLOS ONE，2015年02月26日，Vol.10 No.1，pp.1-14，doi.org/10.1371/journal.pone.0112020

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ５／０５５

Ｇ０１Ｒ ３３／２０－３３／６４

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

磁気共鳴イメージング装置により実行される磁気共鳴イメージング方法であって、
撮像対象の被検体から生体情報を取得し、
非カルテシアン座標系において、前記生体情報の周期に基づいて、互いに回転角度の異なる複数のトラジェクトリを設定し、前記設定された互いに回転角度の異なる複数のトラジェクトリで収集したｋ空間データを配置し、
前記トラジェクトリが異なる複数の収集パターンから、前記生体情報の周期における１周期あたりのｋ空間の分散性を評価することで、収集パターンを選択し、
前記ｋ空間データを再構成して画像を生成することを含み、
前記収集パターンを選択するステップは、
前記１周期あたりのスポーク数の値を算出し、前記値を中心とした範囲における前記分散性を評価することで前記収集パターンを選択する、磁気共鳴イメージング方法。

【請求項2】

前記生体情報は呼吸に関する情報である、請求項１に記載の磁気共鳴イメージング方法。

【請求項3】

前記生体情報は心拍に関する情報である、請求項１に記載の磁気共鳴イメージング方法。

【請求項4】

前記収集パターンを選択するステップは、
前記複数の収集パターンの中から、前記被検体への安全性にかかる指標値に応じて、収集パターンを選択する、
請求項１～３のいずれか一つに記載の磁気共鳴イメージング方法。

【請求項5】

前記収集パターンを選択するステップでは、Ｂを高周波磁場として、前記指標値としてｄＢ／ｄｔを用いる、
請求項４に記載の磁気共鳴イメージング方法。

【請求項6】

前記複数のトラジェクトリでｋ空間データを収集する前に、前記回転角度を決定するパラメータと、前記生体情報の周期とに基づいて、前記ｋ空間データの収集シーケンスのパラメータを調整する、
請求項１に記載の磁気共鳴イメージング方法。

【請求項7】

前記複数のトラジェクトリでｋ空間データを収集する前に、前記回転角度を決定するパラメータと、前記生体情報の周期とに基づいて、前記ｋ空間データの収集シーケンスのパラメータの調整候補値を提示する、
請求項１に記載の磁気共鳴イメージング方法。

【請求項8】

前記複数の収集パターンは、自然数をＮとし、ψ＝１．６１８として、前記回転角度が、１８０°／（Ｎ＋１／ψ）となる、請求項４又は５に記載の磁気共鳴イメージング方法。

【請求項9】

圧縮センシングを用いて前記ｋ空間データを再構成して画像を生成する、
請求項１～８のいずれか一つに記載の磁気共鳴イメージング方法。

【請求項10】

前記非カルテシアン座標系は、ラジアル座標系であり、
前記ラジアル座標系で収集された前記ｋ空間データの一部を用いて再構成を行い動的に前記画像を生成する、請求項１～９のいずれか一つに記載の磁気共鳴イメージング方法。

【請求項11】

撮像対象の被検体から生体情報を取得する取得部と、
非カルテシアン座標系において、前記生体情報の周期に基づいて、互いに回転角度の異なる複数のトラジェクトリを設定する制御部と、
前記設定された互いに回転角度の異なる複数のトラジェクトリで収集したｋ空間データを配置する生成部とを
備え、
前記制御部は、前記トラジェクトリが異なる複数の収集パターンから、前記生体情報の周期おける１周期あたりのｋ空間の分散性を評価することで収集パターンを選択し、
前記生成部は、前記ｋ空間データを再構成して画像を生成し、
前記制御部は、前記１周期あたりのスポーク数の値を算出し、前記値を中心とした範囲における前記分散性を評価することで前記収集パターンを選択
する、磁気共鳴イメージング装置。

【発明の詳細な説明】

【技術分野】

【0001】

本明細書等に開示の実施形態は、磁気共鳴イメージング方法及び磁気共鳴イメージング装置に関する。

【背景技術】

【0002】

ラジアル系の収集シーケンスを用いて磁気共鳴イメージングを行う場合、収集を、互いに角度を黄金角ＧＡ（約１１１．２５°）だけ回転角度をずらしながら、ｋ空間の原点を通る複数のスポークで行う方法が知られている。このような収集方法を取ることで、分散性のよいｋ空間データを得ることができる。

【0003】

しかしながら、例えば、実際の撮像状況下では、呼吸や心拍など、周期的な変動がある。ここで、これらの周期的な変動と、撮像の周期とが同期してしまうと、画質の劣化につながる場合がある。

【先行技術文献】

【特許文献】

【0004】

【文献】特開２０１９－５２８９号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

明細書等に開示の実施形態が解決しようとする課題の一つは、画質を向上させることである。ただし、本明細書等に開示の実施形態により解決される課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を、本明細書等に開示の実施形態が解決する他の課題として位置づけることもできる。

【課題を解決するための手段】

【0006】

実施形態に係る磁気共鳴イメージング方法は、磁気共鳴イメージング装置により実行される磁気共鳴イメージング方法であって、撮像対象の被検体から生体情報を取得し、非カルテシアン座標系において、生体情報の周期に基づいて設定された互いに回転角度の異なる複数のトラジェクトリで収集したｋ空間データを配置し、ｋ空間データを再構成して画像を生成する。

【図面の簡単な説明】

【0007】

【図1】図１は、実施形態に係る磁気共鳴イメージング装置を示す図である。

【図2】図２は、実施形態の背景について説明した図である。

【図3】図３は、第１の実施形態に係る磁気共鳴イメージング装置が行う処理の手順を示したフローチャートである。

【図4A】図４Ａは、比較例に係る磁気共鳴イメージング方法について説明した図である。

【図4B】図４Ｂは、第１の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図4C】図４Ｃは、第１の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図5】図５は、その他の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図6A】図６Ａは、その他の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図6B】図６Ｂは、その他の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図6C】図６Ｃは、その他の実施形態に係る磁気共鳴イメージング方法について説明した図である。

【図6D】図６Ｄは、変形例に係る磁気共鳴イメージング方法について説明した図である。

【図7】図７は、変形例に係る磁気共鳴イメージング方法について説明した図である。

【図8】図８は、第２の実施形態に係る磁気共鳴イメージング装置が行う処理の手順を示したフローチャートである。

【発明を実施するための形態】

【0008】

以下、本発明の実施形態について図面を参照しながら説明する。ここで、互いに同じ構成には共通の符号を付して、重複する説明は省略する。

【0009】

（第１の実施形態）
図１は、第１の実施形態に係る磁気共鳴イメージング装置１００を示すブロック図である。図１に示すように、磁気共鳴イメージング装置１００は、静磁場磁石１０１と、静磁場電源（図示しない）と、傾斜磁場コイル１０３と、傾斜磁場電源１０４と、寝台１０５と、寝台制御回路１０６と、送信コイル１０７と、送信回路１０８と、受信コイル１０９と、受信回路１１０と、シーケンス制御回路１２０（シーケンス制御部）と、コンピューター１３０（「画像処理装置」とも称される）とを備える。なお、磁気共鳴イメージング装置１００に、被検体Ｐ（例えば、人体）は含まれない。また、図１に示す構成は一例に過ぎない。例えば、シーケンス制御回路１２０及びコンピューター１３０内の各部は、適宜統合若しくは分離して構成されてもよい。

【0010】

静磁場磁石１０１は、中空の略円筒形状に形成された磁石であり、内部の空間に静磁場を発生する。静磁場磁石１０１は、例えば、超伝導磁石等である。別の例として、静磁場磁石１０１は、永久磁石でもよい。

【0011】

傾斜磁場コイル１０３は、中空の略円筒形状に形成されたコイルであり、静磁場磁石１０１の内側に配置される。傾斜磁場コイル１０３は、互いに直交するＸ、Ｙ、及びＺの各軸に対応する３つのコイルが組み合わされて形成されており、これら３つのコイルは、傾斜磁場電源１０４から個別に電流の供給を受けて、Ｘ、Ｙ、及びＺの各軸に沿って磁場強度が変化する傾斜磁場を発生する。傾斜磁場コイル１０３によって発生するＸ、Ｙ、及びＺの各軸の傾斜磁場は、例えば、スライス用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、及びリードアウト用傾斜磁場Ｇｒである。傾斜磁場電源１０４は、傾斜磁場コイル１０３に電流を供給する。

【0012】

寝台１０５は、被検体Ｐが載置される天板１０５ａを備え、寝台制御回路１０６による制御の下、天板１０５ａを、被検体Ｐが載置された状態で、傾斜磁場コイル１０３の空洞（撮像口）内へ挿入する。通常、寝台１０５は、長手方向が静磁場磁石１０１の中心軸と平行になるように設置される。寝台制御回路１０６は、コンピューター１３０による制御の下、寝台１０５を駆動して天板１０５ａを長手方向及び上下方向へ移動する。

【0013】

送信コイル１０７は、傾斜磁場コイル１０３の内側に配置され、送信回路１０８からＲＦパルスの供給を受けて、高周波磁場を発生する。送信回路１０８は、対象とする原子の種類及び磁場強度で定まるラーモア（Ｌａｒｍｏｒ）周波数に対応するＲＦパルスを送信コイル１０７に供給する。

【0014】

受信コイル１０９は、傾斜磁場コイル１０３の内側に配置され、高周波磁場の影響によって被検体Ｐから発せられる磁気共鳴信号（以下、必要に応じて、「ＭＲ信号」と呼ぶ）を受信する。受信コイル１０９は、磁気共鳴信号を受信すると、受信した磁気共鳴信号を受信回路１１０へ出力する。

【0015】

なお、上述した送信コイル１０７及び受信コイル１０９は一例に過ぎない。送信機能のみを備えたコイル、受信機能のみを備えたコイル、若しくは送受信機能を備えたコイルのうち、１つ若しくは複数を組み合わせることによって構成されればよい。

【0016】

受信回路１１０は、受信コイル１０９から出力される磁気共鳴信号を検出し、検出した磁気共鳴信号に基づいて磁気共鳴データを生成する。具体的には、受信回路１１０は、受信コイル１０９から出力される磁気共鳴信号をデジタル変換することによって磁気共鳴データを生成する。また、受信回路１１０は、生成した磁気共鳴データをシーケンス制御回路１２０へ送信する。なお、受信回路１１０は、静磁場磁石１０１や傾斜磁場コイル１０３等を備える架台装置側に備えられてもよい。さらには、受信回路１１０の有する機能の一部、例えば磁気共鳴信号のデジタル変換を、受信コイル１０９に設けても構わない。

【0017】

シーケンス制御回路１２０は、コンピューター１３０から送信されるシーケンス情報に基づいて、傾斜磁場電源１０４、送信回路１０８及び受信回路１１０を駆動することによって、被検体Ｐの撮像を行う。ここで、シーケンス情報は、撮像を行うための手順を定義した情報である。シーケンス情報には、傾斜磁場電源１０４が傾斜磁場コイル１０３に供給する電流の強さや電流を供給するタイミング、送信回路１０８が送信コイル１０７に供給するＲＦパルスの強さやＲＦパルスを印加するタイミング、受信回路１１０が磁気共鳴信号を検出するタイミング等が定義される。例えば、シーケンス制御回路１２０は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）等の集積回路、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＭＰＵ（Ｍｉｃｒｏ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等の電子回路である。なお、シーケンス制御回路１２０が実行するパルスシーケンスの詳細については、後述する。

【0018】

さらに、シーケンス制御回路１２０は、傾斜磁場電源１０４、送信回路１０８及び受信回路１１０を駆動して被検体Ｐを撮像した結果、受信回路１１０から磁気共鳴データを受信すると、受信した磁気共鳴データをコンピューター１３０へ転送する。

【0019】

コンピューター１３０は、磁気共鳴イメージング装置１００の全体制御や、画像の生成等を行う。コンピューター１３０は、メモリ１３２、入力装置１３４、ディスプレイ１３５、処理回路１５０を備える。処理回路１５０は、インタフェース機能１３１、制御機能１３３、及び生成機能１３６を備える。

【0020】

第１の実施形態では、インタフェース機能１３１、制御機能１３３、生成機能１３６にて行われる各処理機能は、コンピューターによって実行可能なプログラムの形態でメモリ１３２へ記憶されている。処理回路１５０はプログラムをメモリ１３２から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路１５０は、図１の処理回路１５０内に示された各機能を有することになる。なお、図１においては単一の処理回路１５０にて、インタフェース機能１３１、制御機能１３３、生成機能１３６にて行われる処理機能が実現されるものとして説明するが、複数の独立したプロセッサを組み合わせて処理回路１５０を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。換言すると、上述のそれぞれの機能がプログラムとして構成され、１つの処理回路１５０が各プログラムを実行する場合であってもよい。別の例として、特定の機能が専用の独立したプログラム実行回路に実装される場合であってもよい。なお、図１において、インタフェース機能１３１、制御機能１３３、生成機能１３６は、それぞれ受付部、制御部、生成部の一例である。また、シーケンス制御回路１２０は、シーケンス制御部の一例である。

【0021】

上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）或いは、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。プロセッサはメモリ１３２に保存されたプログラムを読み出し実行することで機能を実現する。

【0022】

また、メモリ１３２にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、寝台制御回路１０６、送信回路１０８、受信回路１１０等も同様に、上記のプロセッサ等の電子回路により構成される。

【0023】

処理回路１５０は、インタフェース機能１３１により、シーケンス情報をシーケンス制御回路１２０へ送信し、シーケンス制御回路１２０から磁気共鳴データを受信する。また、磁気共鳴データを受信すると、インタフェース機能１３１を有する処理回路１５０は、受信した磁気共鳴データをメモリ１３２に格納する。

【0024】

メモリ１３２に格納された磁気共鳴データは、制御機能１３３によってｋ空間に配置される。この結果、メモリ１３２は、ｋ空間データを記憶する。

【0025】

メモリ１３２は、インタフェース機能１３１を有する処理回路１５０によって受信された磁気共鳴データや、制御機能１３３を有する処理回路１５０によってｋ空間に配置されたｋ空間データ、生成機能１３６を有する処理回路１５０によって生成された画像データ等を記憶する。例えば、メモリ１３２は、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等である。

【0026】

入力装置１３４は、操作者からの各種指示や情報入力を受け付ける。入力装置１３４は、例えば、マウスやトラックボール等のポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスである。ディスプレイ１３５は、制御機能１３３を有する処理回路１５０による制御の下、撮像条件の入力を受け付けるためのＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）や、生成機能１３６を有する処理回路１５０によって生成された画像等を表示する。ディスプレイ１３５は、例えば、液晶表示器等の表示デバイスである。

【0027】

処理回路１５０は、制御機能１３３により、磁気共鳴イメージング装置１００の全体制御を行い、撮像や画像の生成、画像の表示等を制御する。例えば、制御機能１３３を有する処理回路１５０は、撮像条件（撮像パラメータ等）の入力をＧＵＩ上で受け付け、受け付けた撮像条件に従ってシーケンス情報を生成する。また、制御機能１３３を有する処理回路１５０は、生成したシーケンス情報をシーケンス制御回路１２０へ送信する。
処理回路１５０は、生成機能１３６により、ｋ空間データをメモリ１３２から読み出し、読み出したｋ空間データにフーリエ変換等の再構成処理を施すことで、画像を生成する。

【0028】

次に、実施形態に係る背景について簡単に説明する。

【0029】

ラジアル系の収集シーケンスを用いて磁気共鳴イメージングを行う場合、収集を、互いに角度を黄金角ＧＡ（約１１１．２５°）または黄金角ＧＡから導かれる角度である１８０°－ＧＡ，３６０°－２×ＧＡ，２×ＧＡなどだけ回転角度をずらしながら、ｋ空間の原点を通る複数のスポークで行う方法が知られている。このような収集方法をとることで、ｋ空間に充填されたスポーク同士の間隔が、最悪でも３種類以下に抑制され、結果として、ｋ空間のどの方向に対しても、まんべんなくｋ空間データを収集することができる。すなわち、分散性のよいｋ空間データを得ることができる。また、これらの収集は、どの時刻においても一定の角度でスポークが回転していくので、撮像中のデータをどの時刻で一部切り出してきても、分散性の良いｋ空間データを収集できるという利点を有する。

【0030】

一方で、黄金角ＧＡだけスポークの回転角度を変化させながら磁気共鳴イメージングを行う場合、黄金角ＧＡは固定された角度であるので、撮像の自由度には限界がある。換言すると、黄金角ＧＡ、または、黄金角ＧＡから導かれる角度である１８０°－ＧＡ，３６０°－２×ＧＡ，２×ＧＡなどのなかから選ばれる固定的なスポークの回転角度を用いただけでは、画質の高い画像を得られない場合がある。例えば、実際の撮像状況下では、呼吸や心拍などの周期的な変動があるので、このような周期的な変動と、撮像の周期とが同期してしまうと、画質の劣化につながる場合がある。すなわち、呼吸や心拍などの変動の周期によっては、黄金角ＧＡを利用した回転角度を設定した撮像であっても、ｋ空間データの分散性が十分に担保できない場合がある。

【0031】

図２に、ｋ空間データの分散性についての例を説明する図が示されている。図２において、収集１０Ａは、シーケンス制御回路１２０が１呼吸周期あたり１１スポークのｋ空間データを収集する場合において、呼気時に収集されるｋ空間データの一例を示している。これに対して、収集１０Ｂは、シーケンス制御回路１２０が１呼吸周期あたり１１スポークのｋ空間データを収集する場合において、吸気時に収集されるｋ空間データの一例を示している。また、収集１１Ａは、シーケンス制御回路１２０が１呼吸周期あたり２０スポークのｋ空間データを収集する場合において、呼気時に収集されるｋ空間データの一例を示している。また、収集１１Ｂは、シーケンス制御回路１２０が１呼吸周期あたり２０スポークのｋ空間データを収集する場合において、吸気時に収集されるｋ空間データの一例を表している。ここで、図２において、点線の一つ一つは、収集される一つ一つのスポークを示す。

【0032】

図２に示されるように、１呼吸周期当たり１１スポークの収集を行う場合、吸気時のみ、あるいは呼気時のみを考えると、収集されるｋ空間の位置が、ｋ空間の特定の方向に偏っている。また、１呼吸周期当たり２０スポークの収集を行う場合、１呼吸周期当たり１１スポークの収集を行う場合と比較すると状況は改善されているものの、吸気時のみあるいは呼気時のみを考えると、収集されるｋ空間の位置が、依然としてｋ空間の特定の方向に偏っている。ここで、１呼吸周期を通じて撮像を行う場合、得られる画像は、呼気時のｋ空間データより得られた画像と、吸気時のｋ空間データより得られた画像とを合成した画像と考えられるが、呼気時と吸気時とでは、動きなどの影響により、撮像対象の位置が変化する。従って、吸気時のみあるいは呼気時のみを考えたときに、収集されるｋ空間の位置がｋ空間の特定の方向に偏っていてｋ空間の分散性が悪い場合、画質が劣化する可能性がある。

【0033】

以上のように、呼吸などの周期的な変動が存在する場合において、呼吸などの変動の周期と、撮像の周期との関係まで考慮すると、黄金角ＧＡを用いたラジアル系のサンプリングを行ったとしても、画質が劣化する可能性がある。一方で、黄金角ＧＡを用いたサンプリングが有する、ｋ空間データの高い分散性という利点があるので、活用する余地がある。

【0034】

かかる背景に基づいて、実施形態に係る磁気共鳴イメージング方法では、処理回路１５０は、制御機能１３３により、非カルテシアン座標系において、スポークの回転角度が異なる複数の収集パターンの中から収集パターンを選択する。シーケンス制御回路１２０は、選択された収集パターンでパルスシーケンスを実行してｋ空間データを収集し、処理回路１５０は、生成機能１３６により、当該ｋ空間データを再構成して画像を生成する。ここで、処理回路１５０は、制御機能１３３により、収集パターンにおけるトラジェクトリのスポークを回転させる角度を、黄金角を一般化した以下の式に基づいて設定する。下式で定義されるスポークの回転角度を、以下便宜上「プラチナ角（Ｐｌａｔｉｎｕｍ Ａｎｇｌｅ： ＰＡ）」と称することにする。

【数1】

【0035】

ここで、Ｎは自然数であり、ＰＡ（Ｎ）は、Ｎに対応するプラチナ角であり、φは黄金比に基づく定数である。ここで、Ｎ＝１でＰＡ（Ｎ）は約１１１．２５°となり上述した黄金角になることから、プラチナ角は黄金角の一般化と考えられる。そのため、シーケンス制御回路１２０は、プラチナ角を用いたラジアルサンプリングを行うことで、黄金角を用いたラジアルサンプリングと同様、高いｋ空間データの分散性を保持することができる。一方で、プラチナ角は黄金角と異なりパラメータＮの自由度があるので、処理回路１５０がパラメータＮを適切に選択することで、シーケンス制御回路１２０が、個々の被検体の生体情報に応じた適切なパルスシーケンスを実行することでき、画質が向上する。

【0036】

かかる構成について、図３及び図４Ａ～４Ｃを用いて説明する。

【0037】

図３は、第１の実施形態に係る磁気共鳴イメージング装置が行う処理の手順を示したフローチャートである。

【0038】

はじめに、ステップＳ１００において、処理回路１５０は、制御機能１３３により、撮像対象の被検体から、例えば呼吸情報収集装置や心電計等を通じて、呼吸に関する情報や、心拍に関する情報等の、生体情報を取得する。

【0039】

続いて、ステップＳ１１０Ａにおいて、処理回路１５０は、制御機能１３３により、対象とする変動の周期を測定し、並行して、スキャン条件から、単位時間当たりに収集されるスポーク数を計算する。例えば、処理回路１５０は、制御機能１３３により、ステップＳ１００で得られた撮像対象の被検体の生体情報である呼吸に関する情報に基づいて、対象とする変動の周期である呼吸周期が、「４秒」であると測定する。また、これと並行して、処理回路１５０は、制御機能１３３により、シーケンス制御回路１２０からスキャン条件を取得し、１秒あたりに収集されるスポーク数を、例えば、「４．０」と算出する。

【0040】

続いて、ステップＳ１２０において、処理回路１５０は、制御機能１３３により、変動１周期当たりのスポーク数を計算する。例えば、対象とする変動の周期である呼吸周期が「４秒」であり、１秒あたりに収集されるスポーク数が「４．０」である場合、処理回路１５０は、制御機能１３３により、変動１周期当たりのスポーク数、すなわち１呼吸周期あたりのスポーク数を、４．０×４＝１６．０と算出する。

【0041】

続いて、ステップＳ１３０において、処理回路１５０は、制御機能１３３により、互いに回転角度の異なる複数のトラジェクトリにおける、回転角度が異なる複数の収集パターンを参照し、参照した収集パターンの中から、適切な収集パターンを選択する。一例として、処理回路１５０は、制御機能１３３により、回転角度であるプラチナ角ＰＡ（Ｎ）が異なる複数のプラチナ角サンプリングパターンを当該複数の収集パターンとして参照し、参照した収集パターンの中から、配置されるｋ空間データの分布に応じて、最適なプラチナ角ＰＡ（Ｎ）となるようなサンプリングパターンを選択する。一例として、処理回路１５０は、制御機能１３３により、プラチナ角ＰＡ（Ｎ）が異なる複数の収集パターンの中から、例えば呼気時など特定の時相のデータを抽出した場合における、ｋ空間の分散性を表す指標（分散指標値）に基づいて、最適なＮを選択することにより、最適な収集パターンを選択する。

【0042】

ここで、ｋ空間の分散性をあらわす指標（分散指標値）は、例えば、黄金角またはプラチナ角サンプリングを行った場合のＳＮＲ（Ｓｉｇｎａｌ ｔｏ Ｎｏｉｓｅ Ｒａｔｉｏ）であるＳＮＲ_{ＧＯＬＤＥＮＲＡＴＩＯ}の、隣接してｋ空間に配置されるスポークの間隔が等しくなるようにサンプリングした場合のＳＮＲであるＳＮＲ_{ｕｎｉｆｏｒｍ}に対する比であり、当該指標は、以下の式で計算できる。

【数2】

【0043】

ここで、Ｐは、データ収集が行われるスポークの本数であり、Δφ_ｉは、黄金角（プラチナ角）サンプリングを行って収集したスポークを、基準となる始線となす角の大きさの順でソートした場合における、ｉ番目のスポークと、ｉ＋１番目のスポークとの角度差である。

【0044】

上述の分散指標値は、ｋ空間に配置されるスポークが、ｋ空間にスポークが等間隔で配置されている状態に近づくほど大きくなり、分散指標値は大きいほど好適であると考えられる。黄金角（プラチナ角）サンプリングを行ってｋ空間に配置された、隣接するスポークがなす角度の分散値が小さいほど、分散指標値は大きくなる。

【0045】

ここで、一呼吸あたりのスポーク数が１８本であると想定した比較例を図４Ａ、図４Ｂ、図４Ｃを用いて説明する。呼吸周期が常に安定しているとは限らないため、例えば１８±２スポークの範囲における分散性指標値の平均が最大となるＮを最適なＮとして選択することとする。

【0046】

図４Ａにおいて、グラフ１は、比較例として、一呼吸周期のうち呼気時の時相において、黄金角サンプリングを行った場合の分散指標値の値を、一呼吸周期あたりのスポーク数の関数としてプロットしたものである。一呼吸周期のうち呼気時など、特定の時相のデータを抽出した場合、１呼吸周期あたりのスポーク数が１６～２０となる領域２において、分散指標値が減少傾向にある。すなわち、領域２では、ｋ空間データに配置されるスポークの分布が偏っており、画質の劣化が想定される。従って、処理回路１５０は、制御機能１３３により、ステップＳ１２０で計算された「変動１周期当たりのスポーク数」において、Ｎを変化させてプラチナ角ＰＡ（Ｎ）でサンプリングをした場合の分散指標値が最大となるようなＮの値を、最適なＮとして選択する。

【0047】

図４Ｂにおいて、グラフ１ａは、一呼吸周期のうち呼気時の時相のデータを抽出した場合における、Ｎ＝４におけるプラチナ角ＰＡ（Ｎ）でシーケンス制御回路１２０がサンプリングをした場合の上述の分散指標値の値を、一呼吸周期当たりのスポーク数の関数としてプロットしたものである。図４Ｂに示されるように、Ｎ＝４においては、１呼吸周期当たりのスポーク数が１６～２０となる領域では分散指標値が減少傾向にあることから、この領域ではＮ＝４は適切なサンプリングでない。

【0048】

一方、図４Ｃにおいて、グラフ１ｂは、一呼吸周期のうち呼気時の時相のデータを抽出した場合における、Ｎ＝５におけるプラチナ角ＰＡ（Ｎ）でシーケンス制御回路１２０がサンプリングをした場合の上述の分散指標値の値を、一呼吸周期当たりのスポーク数の関数としてプロットしたものである。図４Ｃに示されるように、Ｎ＝５においては、１呼吸周期当たりのスポーク数が１６～２０となる領域２では、Ｎ＝４の場合と比較して分散指標値の値が改善されていることから、この領域ではＮ＝５はＮ＝４と比較して望ましいサンプリングことがわかる。処理回路１５０は、制御機能１３３により、さまざまなＮについて同様の計算を行った上で、Ｎ＝５を、最適な収集パターンとして選択する。

【0049】

なお、図４Ａ、図４Ｂ、図４Ｃを用いた説明においては、一呼吸周期のうち、呼気時の時相データを用い、分散指標値が最大となるＮを収集パターンの決定に用いている。これは、呼気時において最適なＮの値は、吸気時のデータにおいても最適となる場合が多いという仮定に基づくものである。したがって、吸気時のデータにおいて最適となるＮを別途計算することを排除することは意図しておらず、吸気時のデータにおいて最適な収集パターンを別途計算しても良い。

【0050】

続いて、ステップＳ１４０において、シーケンス制御回路１２０は、ステップＳ１３０で選択された収集パターンに基づいて、非カルテシアン座標系において、生体情報の周期に基づいてステップＳ１３０で設定された互いに回転角度の異なる複数のトラジェクトリでデータ収集を行う。例えば、シーケンス制御回路１２０は、ステップＳ１３０で選択されたＮに係るプラチナ角ＰＡ（Ｎ）により定められる、互いに回転角度の異なる複数のトラジェクトリで、２Ｄラジアル収集のパルスシーケンスを実行してデータ収集を行う。

【0051】

続いて、ステップＳ１５０において、処理回路１５０は、生成機能１３６により、ステップＳ１４０で収集したｋ空間データを配置し、当該ｋ空間データを再構成して磁気共鳴画像を生成する。

【0052】

以上述べたように、第１の実施形態においては、処理回路１５０が制御機能１３３により最適な収集パターンを分散指標値に基づいて選択した上で、シーケンス制御回路１２０が選択されたプラチナ角でラジアルサンプリングを行った。これにより、呼吸などの周期的な変動であって、例えば被検体ごとに個人差があるような状況等において、個々の被検体について、ｋ空間データの分散性を向上することができ、画質を向上することができる。

【0053】

（変形例１）
第１の実施形態では、図３のステップＳ１００及びステップＳ１１０Ａにおいて、処理回路１５０が被検体の生体情報を、例えば呼吸情報収集装置や心電計等を通じて収集する場合について説明した。実施形態はこれに限られず、処理回路１５０は、制御機能１３３により、ユーザから生体情報の入力を受け付けても良い。図５に、そのような場合のフローチャートが示されている。

【0054】

変形例１にかかる実施形態では、図３におけるステップＳ１００及びステップＳ１１０Ａの代わりに、ステップＳ１１０Ｂにおいて、処理回路１５０は、制御機能１３３により、入力装置１３４を通じてユーザから被検体の呼吸特性等の生体情報を受け付けるとともに、スキャン条件から、単位時間当たりに収集されるスポーク数を計算する。ステップＳ１２０～ステップＳ１５０においては、第１の実施形態と同様の処理を行う。

【0055】

制御機能１３３は、例えば、呼吸の周期をユーザに入力させ、生体情報として受け付ける。また、生体情報は、「呼吸が浅い」「通常」「呼吸が深い」といった定性的な選択肢を入力装置１３４がユーザに提示し、ユーザがいずれかの選択肢を選ぶことにより、制御機能１３３が選択肢に関連付けられた呼吸の周期を、生体情報として受け付けてもよい。なお、呼吸の周期に関して入力例を示したが、他の生体情報を受け付けるようにしてもよい。例えば、心拍や脈波の情報を入力装置１３４で受け付けられるようにする場合は、事前に計測しておいた心拍数を入力できるようにしたり、「脈が早い」「通常」「脈が遅い」といった選択肢から選べるようにしたりすることが例として考えられる。

【0056】

変形例１にかかる実施形態によれば、呼吸情報収集装置や心電計といった生体情報を収集する装置を被検体に装着していなくても、ユーザの入力により最適な収集パターンを選択することができる。また、被検体にとって生体情報を収集する装置を装着する負荷がなくなるとともに、装置を被検体に装着させる手順がなくワークフローの改善にもつながる。

【0057】

（変形例２）
第１の実施形態では、ステップＳ１３０において、処理回路１５０の制御機能１３３が分散指標値の値に基づいて最適な収集パターンを選択した。実施形態はこれに限られず、ステップＳ１３０において、処理回路１５０は、制御機能１３３により、複数の収集パターンの中から、被検体への安全性に係る指標値に応じて、適切な収集パターンを選択してもよい。

【0058】

変形例２にかかる実施形態では、処理回路１５０は、制御機能１３３により、安全性に係る指標値の一例として、磁場の時間変化率ｄＢ／ｄｔに基づいて、収集パターンを選択する。例えば、スポークの回転角度が大きくなると傾斜磁場の時間変化率が大きくなるため、ｄＢ／ｄｔを考慮して収集パターンを選択することが必要となる場合がある。例えば、処理回路１５０は、制御機能１３３により、ｄＢ／ｄｔが閾値を下回る収集パターンを選択する。閾値は、例えば安全規格で定められた最大ｄＢ／ｄｔの値や任意に設定した値を用いることができる。また、別の例として、処理回路１５０は、制御機能１３３により、ｄＢ／ｄｔを最小化する収集パターンを選択してもよい。さらに、別の例として、処理回路１５０は、制御機能１３３により、分散指標値の値と、安全性に係る指標値との両方とに基づいて、適切な収集パターンを選択してもよい。一例として、処理回路１５０は、制御機能１３３により、安全性に係る指標値が正常範囲に収まる収集パターンの中で、もっとも分散指標値が良い収集パターンを、適切な収集パターンとして選択する。

【0059】

変形例２にかかる実施形態によれば、安全性に係る指標値を用いて最適な収集パターンを選択することができる。これによれば、ユーザが意識せずとも、ＭＲＩ装置に要求される安全規格に準拠した範疇で画質の高い画像を得ることができる。また、安全性に係る指標値と分散指標値とを併用して収集パターンを決定すれば、最適な収集パターンを決定するまでのワークフローをより短縮することができる。

【0060】

（変形例３）
ステップ１３０で選択される最適な収集パターン及びステップＳ１４０で実行される撮像については、第１の実施形態で説明したものに限られない。変形例３に係る実施形態では、収集パターンのバリエーションについて示す。

【0061】

プラチナ角ＰＡ（Ｎ）を用いた収集のバリエーションを図６Ａ～図６Ｄに示した。図６Ａは、ｋ番目の収集のスポーク３とｋ＋１番目の収集のスポーク４との間の回転角度５ａが、プラチナ角ＰＡ（Ｎ）となる基本パターンを示す。これに対して、図６Ｂ～図６Ｄは、図６Ａの基本パターンの拡張である。図６Ｂは、ｋ番目の収集のスポーク３とｋ＋１番目の収集のスポーク４との間の回転角度５ｂが、１８０°－ＰＡ（Ｎ）である場合のパターンを示している。図６Ｃは、ｋ番目の収集のスポーク３とｋ＋１番目の収集のスポーク４との間の回転角度５ｃが、２×ＰＡ（Ｎ）である場合のパターンを示している。図６Ｄは、ｋ番目の収集のスポーク３とｋ＋１番目の収集のスポーク４との間の回転角度５ｄが、３６０°－２×ＰＡ（Ｎ）である場合のパターンを示している。

【0062】

シーケンス制御回路１２０は、プラチナ角ＰＡ（Ｎ）そのものではなく、プラチナ角ＰＡ（Ｎ）に近似する角度である疑似プラチナ角（Ｐｓｅｕｄｏ Ｐｌａｔｉｎｕｍ Ａｎｇｌｅ）を用いて収集を行ってもよい。ここで、疑似プラチナ角は、例えばプラチナ角に近い有理数であって、分子及び分数が比較的小さい整数となる角度として構成される。また、別の例として、シーケンス制御回路１２０は、プラチナ角の回転角度ずつ回転したｎ本のスポークを、基準となる始線となす角の大きさの順でソートし、そのｉ番目のスポークの角度が、３６０°×（ｉ／ｎ）となるように等間隔に改めて配置されたスポークについて、収集シーケンスを実行してもよい。

【0063】

また、ステップ１３０で選択される最適な収集パターン及びステップＳ１４０で実行される撮像については、そのほかにも、さまざまなバリエーションが考えられる。

【0064】

例えば、ステップＳ１４０において、シーケンス制御回路１２０は、単純Ｎ分割、すなわち、ｉ番目のスポークの角度が３６０°×（ｉ／Ｎ）となるようなスポークで収集シーケンスを実行してもよい。この場合、ステップＳ１３０において、処理回路１５０は、制御機能１３３により、各Ｎについて分散指標値を算出し、算出された分散指標値に基づいて、最適な収集パターンを選択する。

【0065】

また、ステップＳ１４０において、シーケンス制御回路１２０は、Ｎ／Ｍ周期での分割、すなわち、ｉ番目のスポークの角度が３６０°×（ｉ／Ｎ＊Ｍ）となるようなスポークで収集シーケンスを実行してもよい。この場合、ステップＳ１３０において、処理回路１５０は、制御機能１３３により、各Ｎ及びＭについて分散指標値を算出し、算出された分散指標値に基づいて、最適な収集パターンを選択する。

【0066】

また、ステップＳ１４０において、シーケンス制御回路１２０は、Ｇｒｏｕｐｅｄ Ｇｏｌｄｅｎ Ａｎｇｌｅ、すなわち、１セットの複数の収集を、黄金角ずつ角度をずらしながら収集するシーケンスで収集シーケンスを実行してもよい。図７に、１セットが３スポークである場合のＧｒｏｕｐｅｄ Ｇｏｌｄｅｎ Ａｎｇｌｅの例が示されている。図７において、互いに１２０度間隔で並んだスポーク７ａ，７ｂ，７ｃ と、スポーク８ａ，８ｂ，８ｃとがそれぞれ１セットの収集となり、例えばスポーク７ａとスポーク８ａとの間の角度が黄金角となる。シーケンス制御回路１２０は、例えばスポーク７ａ，７ｂ，７ｃで１セットの収集を行い、続いて、スポーク８ａ，８ｂ，８ｃで１セットの収集を行う。ステップＳ１３０において、処理回路１５０は、制御機能１３３により、１セットがｎスポークである場合のＧｒｏｕｐｅｄ Ｇｏｌｄｅｎ Ａｎｇｌｅでの収集において、各ｎにおいて分散指標値を算出し、算出された分散指標値に基づいて、最適な収集パターンを選択する。

【0067】

また、ステップＳ１４０において、シーケンス制御回路１２０は、Ｂｉｔ Ｒｅｖｅｒｓｅｄ Ｏｒｄｅｒで収集シーケンスを実行してもよい。ここで、例えばＮ，Ｍ，Ｌを自然数としてこれらの自然数の組（Ｎ，Ｍ，Ｌ）で特徴づけられるＢｉｔ Ｒｅｖｅｒｓｅｄ Ｏｒｄｅｒの収集シーケンスとは、第１ステップとして、１周をＮ分割したＮ本のスポークで順次収集を行い、第２ステップとして、１周をＮ×Ｍ分割したスポークのうち、第１ステップで収集されなかったスポークで順次収集を行い、第３ステップとして、１周をＮ×Ｍ×Ｌ分割したスポークのうち、第２ステップまでで収集されなかったスポークで順次収集を行う収集シーケンスであって、特定の順序で行われる収集シーケンスである。ステップＳ１３０において、処理回路１５０は、制御機能１３３により、Ｂｉｔ Ｒｅｖｅｒｓｅｄ Ｏｒｄｅｒの収集シーケンスを特徴づけるパラメータを変化させたときの分散指標値を算出し、算出された分散指標値に基づいて、最適な収集パターンを選択する。

【0068】

また、ステップＳ１３０の処理の別の例として、処理回路１５０は、制御機能１３３により、黄金角とそのバリエーションごとに分散指標値を算出し、算出された分散指標値に基づいて、最適な収集パターンを選択してもよい。例えば、処理回路１５０は、制制御機能１３３により、黄金角ＧＡと、そのバリエーションである１８０°－ＧＡ， ２×ＧＡ， ３６０°－２×ＧＡのそれぞれの回転角度のパルスシーケンスごとに、分散指標値を算出し、算出された分散指標値に基づいて、最適なバリエーションを選択してもよい。

【0069】

また、実施形態では、シーケンス制御回路１２０が２Ｄラジアル収集を行う場合について説明したが、実施形態は、２Ｄラジアル収集に限られず、一定の周期でエンコード方向を変えながら非カルテシアン座標で収集を行うさまざまな収集シーケンス、例えば、２Ｄスパイラル収集などの２Ｄ系収集シーケンス、及び、Ｓｔａｃｋ Ｏｆ Ｓｔａｒｓ収集、ＵＴＥ（ｕｌｔｒａｓｈｏｒｔ ＴＥ）撮像などで用いられるＫｏｏｓｈ Ｂａｌｌ型の収集、Ｃｏｎｅ Ｔｒａｊｅｃｔｏｒｙでの収集などに例示される、３Ｄ系収集シーケンスに適用可能である。また、複数の並行したトラジェクトリが構成する帯状領域であるブレードを回転させて収集するＪＥＴやＰＲＯＰＥＬＬＥＲ（Ｐｅｒｉｏｄｉｃａｌｌｙ Ｒｏｔａｔｅｄ Ｏｖｅｒｌａｐｐｉｎｇ Ｐａｒａｌｌｅｌ Ｌｉｎｅｓ ｗｉｔｈ ｅｎｈａｎｃｅｄ ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ ）などの２Ｄ系収集シーケンスにも適用可能である。つまり、１本のスポークを収集するたびに回転角度を変化させる場合に限らず、実施形態に係る発明は適用可能である。

【0070】

変形例３にかかる実施形態によれば、黄金角、プラチナ角それぞれを利用した回転角度を拡張させた収集パターンを設定できるため、撮像シーケンスの種類に柔軟に対応することが可能となる。また、黄金角を用いた公知のいくつかの収集パターン間で、被検体の生体情報およびスキャン条件を用いて最適なパターンを選択することも可能であり、公知技術を用いた撮像法を高画質化することも可能である。

【0071】

（再構成方法・撮像方法について）
第１の実施形態に係る発明は、圧縮センシング（ＣＳ：Ｃｏｍｐｒｅｓｓｅｄ Ｓｅｎｓｉｎｇ）などの各種画像再構成アルゴリズムやＤｙｎａｍｉｃ撮像と併用してもよい。

【0072】

圧縮センシングなど、画像再構成アルゴリズムによっては、前述の分散指標値が画質の評価指標として相応しくない場合があり得る。その場合には、再構成アルゴリズムに適応した別の指標値を設定し、それに基づき最適な収集パターンを選択してもよい。

【0073】

ステップＳ１３０で選択された最適な収集パターンは、単位時間あたりの角度の増加の割合が一定であるので、どの時刻でデータを切り出してきてもｋ空間データの分散性が良好となる。従って、Ｄｙｎａｍｉｃ再構成などのようにラジアルスキャンで取得したデータを複数のフレームに分割し、別個の画像として再構成する場合に適している。すなわち、非カルテシアン座標系として、ラジアル座標系を選択する場合、処理回路１５０は、生成機能１３６により、当該ラジアル座標系で収集されたｋ空間データの一部を用いて再構成を行い動的に画像を生成してもよい。

【0074】

（第２の実施形態）
第１の実施形態では、ステップＳ１３０において、処理回路１５０は、制御機能１３３により、分散指標値の値に基づいて、最適な収集パターンを選択した。実施形態はこれに限られず、処理回路１５０は、制御機能１３３により、ステップＳ１４０で実行されるパルスシーケンスの最適化を行ってもよい。

【0075】

すなわち、第２の実施形態においては、処理回路１５０は、制御機能１３３により、プラチナ角Ｐ（Ｎ）におけるＮの値を固定したまま、ＴＲ（Ｒｅｐｅｔｉｔｉｏｎ Ｔｉｍｅ）やＴＥ（Ｅｃｈｏ Ｔｉｍｅ）などのパラメータが最適になるように、ステップＳ１４０で実行されるパルスシーケンスの最適化を行う。

【0076】

かかる処理の一例が、図８に示されている。

【0077】

第１の実施形態と同様に、ステップＳ１００において、処理回路１５０は、制御機能１３３により、撮像対象の被検体から、呼吸に関する情報や、心拍に関する情報等の、生体情報を取得する。続いて、ステップＳ１１０Ｃにおいて、処理回路１５０は、制御機能１３３により、対象とする変動の周期を測定する。続いて、ステップＳ１２０において、処理回路１５０は、制御機能１３３により、変動１周期当たりのスポーク数を計算する。

【0078】

続いて、ステップＳ１３０Ｂにおいて、処理回路１５０は、制御機能１３３により、分散指標値等の値に基づいて、実行されるパルスシーケンスのパラメータを変化させ、ステップＳ１４０で実行されるパルスシーケンスの最適化を行う。例えば、処理回路１５０は、制御機能１３３により、ステップＳ１４０において互いに回転角度の異なる複数のトラジェクトリでｋ空間データを収集する前に、回転角度を決定するパラメータと、ステップＳ１１０Ａで測定された生体情報の周期と、ステップＳ１４０で実行されるｋ空間データの収集シーケンスのパラメータごとに算出された分散指標値の値に基づいて、ＴＲ（Ｒｅｐｅｔｉｔｉｏｎ Ｔｉｍｅ）やＴＥ（Ｅｃｈｏ Ｔｉｍｅ）等のｋ空間データの収集シーケンスのパラメータを調整する。一例として、処理回路１５０は、制御機能１３３により、算出された分散指標値の値が、最も好ましい値となるようなＴＲやＴＥの値を、ステップＳ１４０において実行される収集シーケンスのパラメータとして算出する。このように、収集シーケンスのパラメータごとに分散指標値が算出されることで、シーケンス制御回路１２０は、変動１周期あたりのスポーク数が適切な値となるような収集シーケンスのパラメータを算出することができる。

【0079】

また、処理回路１５０は、制御機能１３３により、ディスプレイ１３５を通じてユーザにｋ空間データの収集シーケンスのパラメータの調整候補値を提示してもよい。

【0080】

続いて、ステップＳ１４０において、シーケンス制御回路１２０は、ステップＳ１３０で選択された最適なＴＲやＴＥに基づいて、非カルテシアン座標系においてパルスシーケンスを実行してデータ収集を行う。

【0081】

続いて、ステップＳ１５０において、処理回路１５０は、生成機能１３６により、ステップＳ１４０で収集したｋ空間データを配置し、当該ｋ空間データを再構成して磁気共鳴画像を生成する。

【0082】

このように、第２の実施形態によれば、例えば、呼吸など、周期的な変動がある場合において、実行するパルスシーケンスを最適化することができ、これにより画質を向上することができる。

【0083】

（プログラム）
また、上述した実施形態の中で示した処理手順に示された指示は、ソフトウェアであるプログラムに基づいて実行されることが可能である。汎用コンピューターが、このプログラムを予め記憶しておき、このプログラムを読み込むことにより、上述した実施形態の磁気共鳴イメージング装置１００による効果と同様の効果を得ることも可能である。上述した実施形態で記述された指示は、コンピューターに実行させることのできるプログラムとして、磁気ディスク（フレキシブルディスク、ハードディスクなど）、光ディスク（ＣＤ－ＲＯＭ、ＣＤ－Ｒ、ＣＤ－ＲＷ、ＤＶＤ－ＲＯＭ、ＤＶＤ±Ｒ、ＤＶＤ±ＲＷなど）、半導体メモリ、又はこれに類する記録媒体に記録される。コンピューター又は組み込みシステムが読み取り可能な記憶媒体であれば、その記憶形式は何れの形態であってもよい。コンピューターは、この記録媒体からプログラムを読み込み、このプログラムに基づいてプログラムに記述されている指示をＣＰＵで実行させれば、上述した実施形態の磁気共鳴イメージング装置１００と同様の動作を実現することができる。また、コンピューターがプログラムを取得する場合又は読み込む場合は、ネットワークを通じて取得又は読み込んでもよい。

【0084】

また、記憶媒体からコンピューターや組み込みシステムにインストールされたプログラムの指示に基づきコンピューター上で稼働しているＯＳ（Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）や、データベース管理ソフト、ネットワーク等のＭＷ（Ｍｉｄｄｌｅｗａｒｅ）等が、上述した実施形態を実現するための各処理の一部を実行してもよい。更に、記憶媒体は、コンピューターあるいは組み込みシステムと独立した媒体に限らず、ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）やインターネット等により伝達されたプログラムをダウンロードして記憶又は一時記憶した記憶媒体も含まれる。また、記憶媒体は１つに限られず、複数の媒体から、上述した実施形態における処理が実行される場合も、実施形態における記憶媒体に含まれ、媒体の構成は何れの構成であってもよい。

【0085】

なお、実施形態におけるコンピューター又は組み込みシステムは、記憶媒体に記憶されたプログラムに基づき、上述した実施形態における各処理を実行するためのものであって、パソコン、マイコン等の１つからなる装置、複数の装置がネットワーク接続されたシステム等の何れの構成であってもよい。また、実施形態におけるコンピューターとは、パソコンに限らず、情報処理機器に含まれる演算処理装置、マイコン等も含み、プログラムによって実施形態における機能を実現することが可能な機器、装置を総称している。

【0086】

以上述べた少なくとも一つの実施形態によれば、画質を向上させることができる。

【0087】

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

【符号の説明】

【0088】

１２０ シーケンス制御回路
１３３ 制御機能
１３６ 生成機能
１５０ 処理回路

【図1】

【図2】

【図3】

【図4A】

【図4B】

【図4C】

【図5】

【図6A】

【図6B】

【図6C】

【図6D】

【図7】

【図8】