(57)【特許請求の範囲】
【請求項１】
被検体の心拍を表す生体信号のトリガポイントに同期してパルスシーケンスを繰り返し実行し、前記被検体からデータを収集する磁気共鳴イメージング装置であって、
前記生体信号のトリガポイントと前記パルスシーケンスとの間の遅延時間を決定する遅延時間決定手段であって、前記被検体の心拍が変動した場合、前記遅延時間の値を更新する遅延時間決定手段と、
前記被検体の心位相に基づいて、前記被検体から収集されたデータを取り直すか否かを決定するデータ取直し決定手段と、
を有し、
前記パルスシーケンスは、
前記生体信号のトリガポイントから前記遅延時間が経過した時点で送信されるＲＦパルスと、
前記ＲＦパルスが送信された後に前記被検体からデータを収集するためのデータ収集シーケンスであって、ｋ空間の低周波領域のデータを、データ収集シーケンスの始端側で収集するためのデータ収集シーケンスと、
を有し、
前記データ取直し決定手段は、
前記データ収集シーケンスの開始前における前記被検体の心位相に基づいて、前記被検体から収集されたデータを取り直すか否かを決定する、磁気共鳴イメージング装置。
【請求項２】
被検体の心拍を表す生体信号のトリガポイントに同期してパルスシーケンスを繰り返し実行し、前記被検体からデータを収集する磁気共鳴イメージング装置であって、
前記生体信号のトリガポイントと前記パルスシーケンスとの間の遅延時間を決定する遅延時間決定手段であって、前記被検体の心拍が変動した場合、前記遅延時間の値を更新する遅延時間決定手段と、
前記被検体の心位相に基づいて、前記被検体から収集されたデータを取り直すか否かを決定するデータ取直し決定手段と、
を有し、
前記パルスシーケンスは、
前記生体信号のトリガポイントから前記遅延時間が経過した時点で送信されるＲＦパルスと、
前記ＲＦパルスが送信された後に前記被検体からデータを収集するためのデータ収集シーケンスであって、ｋ空間の低周波領域のデータを、データ収集シーケンスの終端側で収集するためのデータ収集シーケンスと
を有し、
前記データ取直し決定手段は、
前記データ収集シーケンスの終了後における前記被検体の心位相に基づいて、前記被検体から収集されたデータを取り直すか否かを決定する、磁気共鳴イメージング装置。
【請求項３】
前記被検体の心拍が変動したか否かを判断する判断手段を有し、
前記遅延時間決定手段は、
前記判断手段が前記被検体の心拍が変動したと判断した場合、前記遅延時間の値を更新する、請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項４】
前記判断手段は、
前記被検体の心拍数に基づいて、前記被検体の心拍が変動したか否かを判断する、請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
前記ＲＦパルスは反転パルスである、請求項４に記載の磁気共鳴イメージング装置。
【請求項６】
前記遅延時間決定手段は、前記データ収集シーケンスが心拡張期で実行されるように、前記遅延時間を決定する、請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記遅延時間決定手段は、以下の式を用いて前記遅延時間を決定する、請求項６に記載の磁気共鳴イメージング装置。
ＴＤ１＝ｋ・（６０／ＨＲ）＋ｕ・（６０／ＨＲ）−ＴＩ
ここで、ＴＤ１：前記遅延時間
ＨＲ：前記被検体の心拍数
ＴＩ：前記反転パルスの反転時間
ｋ、ｕ：係数
ただし、係数ｋは、以下の式が成り立つときの整数値である。
ＴＩ−ｕ・（60/HR）＜ｋ・（60/HR）＜ＴＩ−（ｕ−１）・（60/HR）
【請求項８】
前記遅延時間決定手段は、以下の式を用いて前記遅延時間を決定する、請求項６に記載の磁気共鳴イメージング装置。
ＴＤ１＝ｕ・（６０／ＨＲ）−ＴＩ
ここで、ＴＤ１：前記遅延時間
ＨＲ：前記被検体の心拍数
ＴＩ：前記ＲＦパルスの反転時間
ｕ：係数
【請求項９】
前記生体信号に基づいて前記被検体の心拍数を算出する心拍数算出手段を有する、請求項４〜８のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項１０】
前記パルスシーケンスは、前記ＲＦパルスと前記データ収集シーケンスとの間に、一つ以上の別のＲＦパルスを有する、請求項１〜９のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項１１】
前記生体信号のトリガポイントは、前記被検体の脈波信号のピーク、又は前記被検体の心電信号のＲ波である、請求項１〜１０のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体からデータを収集する磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
動脈血を撮影する方法として、ＦＳＥ（Fast Spin Echo）を用いた非造影ＭＲＡ（Magnetic
Resonance Angiography）が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開2008-086748号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＦＳＥを用いた非造影ＭＲＡでは、心収縮期のときにデータを収集してしまうと、動脈血を強調して描出することができない場合がある。したがって、データ収集は、できるだけ心拡張期のときの行う必要がある。そこで、心拡張期のときにデータ収集が行われるように、心拍に同期して撮影することが行われている。しかし、撮影中に被検体の心拍が変動してしまうと、心拡張期にデータを収集することができない場合がある。
【０００５】
本発明は、上記の事情に鑑み、被検体の心拍が変動しても所望の心位相のときにデータを収集することが可能な磁気共鳴イメージング装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の問題を解決する本発明の磁気共鳴イメージング装置は、
被検体の心拍を表す生体信号のトリガポイントに同期してパルスシーケンスを繰り返し実行し、前記被検体からデータを収集する磁気共鳴イメージング装置であって、
前記生体信号のトリガポイントと前記パルスシーケンスとの間の遅延時間を決定する遅延時間決定手段を有し、
前記遅延時間決定手段は、前記被検体の心拍が変動した場合、前記遅延時間の値を更新する。
【発明の効果】
【０００７】
本発明では、被検体の心拍が変動した場合、被検体の心拍とパルスシーケンスとの間の遅延時間の値を更新している。したがって、被検体の心拍が変動しても所望の心位相のときにデータを収集することが可能となる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１の実施形態の磁気共鳴イメージング装置の概略図である。
【図２】被検体１４の撮影部位の一例を示す図である。
【図３】被検体１４を撮影するときに使用されるパルスシーケンスＰＳの一例を示す図である。
【図４】心拍同期法の説明図である。
【図５】ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝６０の場合の係数ｋの値について説明する図である。
【図６】ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝４０の場合の係数ｋの値について説明する図である。
【図７】係数ｋ＝０の場合のＰＧ信号とパルスシーケンスＰＳとの関係について説明する図である。
【図８】撮影中に被検体１４の心拍が変動した様子を示す図である。
【図９】ＭＲＩ装置１の処理フローを示す図である。
【図１０】図９の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【図１１】図９の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【図１２】図９の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【図１３】心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃが、心拍数ａのときのＰＧ信号のピーク間隔ＰＰint＝ｗａより短くなった場合の説明図である。
【図１４】被検体１４に不整脈が発生した場合のＰＧ信号とパルスシーケンスＰＳとの関係を示す図である。
【図１５】データを取り直すときのフローを説明する図である。
【図１６】図１５の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【図１７】データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’を求めることが好ましい場合の説明図である。
【発明を実施するための形態】
【０００９】
以下、図面を参照しながら、発明を実施するための形態を詳細に説明するが、本発明は、以下の実施形態に限定されることはない。
【００１０】
（１）第１の実施形態
図１は、本発明の第１の実施形態の磁気共鳴イメージング装置の概略図である。
【００１１】
磁気共鳴イメージング装置（以下、ＭＲＩ（Magnetic Resonance Imaging）装置と呼ぶ）１は、磁場発生装置２、テーブル３、脈波センサ４、受信コイル５などを有している。
【００１２】
磁場発生装置２は、被検体１４が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は勾配パルスを印加し、送信コイル２４はＲＦパルスを送信する。
【００１３】
テーブル３は、被検体１４を搬送するためのクレードル３１を有している。クレードル３１によって、被検体１４はボア２１に搬送される。
【００１４】
脈波センサ４は、被検体１４の脈波を検出する。
【００１５】
受信コイル５は、被検体１４の頭部１４ａに取り付けられており、頭部１４ａからの磁気共鳴信号を受信する。
【００１６】
ＭＲＩ装置１は、更に、シーケンサ６、送信器７、勾配磁場電源８、ＰＧユニット９、受信器１０、中央処理装置１１、入力装置１２、および表示装置１３を有している。
【００１７】
シーケンサ６は、中央処理装置１１の制御を受けて、ＲＦパルスの情報（中心周波数、バンド幅など）を送信器７に送り、勾配磁場の情報（勾配磁場の強度など）を勾配磁場電源８に送る。
【００１８】
送信器７は、シーケンサ６から送られた情報に基づいて、送信コイル２４を駆動する。
【００１９】
勾配磁場電源８は、シーケンサ６から送られた情報に基づいて、勾配コイル２３を駆動する。
【００２０】
ＰＧユニット９は、脈波センサ４からのＰＧ（Peripheral Gating）信号にデジタル化処理を含む各種の処理を施し、中央処理装置１１に伝送する。
【００２１】
受信器１０は、受信コイル５で受信された磁気共鳴信号を信号処理し、中央処理装置１１に伝送する。
【００２２】
中央処理装置１１は、シーケンサ６および表示装置１３に必要な情報を伝送したり、受信器１０から受け取った信号に基づいて画像を再構成するなど、ＭＲＩ装置１の各種の動作を実現するように、ＭＲＩ装置１の各部の動作を総括する。また、中央処理装置１１は、ＰＧユニット９からの信号に基づいて、パルスシーケンスを実行するときの遅延時間ＴＤ１（後述する図４参照）を決定する。中央処理装置１１は、例えばコンピュータ（computer）によって構成される。尚、中央処理装置１１は、特許請求の範囲における遅延時間決定手段、判断手段、データ取直し決定手段、および心拍数算出手段の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２３】
入力装置１２は、オペレータ１５の操作に応じて、種々の命令を中央処理装置１１に入力する。
【００２４】
表示装置１３は種々の情報を表示する。
【００２５】
上記のように構成されたＭＲＩ装置１を用いて、被検体１４を撮影する。
【００２６】
図２は、被検体１４の撮影部位の一例を示す図、図３は、被検体１４を撮影するときに使用されるパルスシーケンスＰＳの一例を示す図である。
【００２７】
本実施形態では、被検体１４の頭部１４ａおよび頸部１４ｂを撮影する。ＭＲＩ装置１は、心臓１４ｃからの動脈血１４ｄの流入効果を利用して、撮像領域Ｒにおける動脈血１４ｄを描出するためのパルスシーケンスＰＳを実行する。
【００２８】
パルスシーケンスＰＳは、選択反転パルスＳＩＲ（Selective Inversion Recovery）およびデータ収集シーケンスＤＡＱを有している。
【００２９】
選択反転パルスＳＩＲは、被検体１４の撮像領域Ｒ（図２参照）の血液（動脈血および静脈血）の縦磁化を反転させるパルスである。選択反転パルスＳＩＲの反転時間ＴＩは、選択反転パルスＳＩＲによって血液の縦磁化を反転させてから、反転した血液の縦磁化がヌルポイントに到達するまでの時間に設定されている。反転時間ＴＩは、例えば１（sec）〜１．５（sec）程度の値である。反転時間ＴＩの間に、心臓１４ｃからの動脈血１４ｄが撮像領域Ｒに流入するので、反転時間ＴＩが経過した時点でデータ収集シーケンスＤＡＱを実行することによって、動脈血が強調して描出されるとともに静脈血が抑制されたＭＲ画像を得ることができる。データ収集シーケンスＤＡＱは、例えばＦＳＥ（Fast Spin Echo）や、ＳＳＦＰ（Steady
State Free Precession）である。
【００３０】
動脈血１４ｄの流入効果を利用して動脈血１４ｄを描出する場合、心臓１４ｅの心収縮期にデータを収集してしまうと、フローボイド（Flow Void）によって動脈血１４ｄを高コントラストで描出することができなくなるので、データ収集シーケンスＤＡＱはできるだけ心拡張期に実行することが望まれる。そこで、本実施形態では、心拍同期法によってパルスシーケンスＰＳを実行している。以下に、心拍同期法について説明する。
【００３１】
図４は、心拍同期法の説明図である。
図４（ａ）は、ＰＧ信号を示す図であり、図４（ｂ）は、パルスシーケンスＰＳを示す図である。
【００３２】
上記のように、データ収集シーケンスＤＡＱは、心拡張期ＤＰに実行されることが望ましい。心拡張期ＤＰの開始時点は、ＰＧ信号のピークＰｎ＋ｋからの時間Ｔｓによって規定することができ、心拡張期ＤＰの終了時点は、ＰＧ信号のピークＰｎ＋ｋからの時間Ｔｅによって規定することができる。時間ＴｓおよびＴｅは、ＰＧ信号のピーク間隔ＰＰintの値を用いて、例えば、以下の式で表すことができる。
Ｔｓ＝ｖ１×ＰＰint  ・・・（１）
Ｔｅ＝ｖ２×ＰＰint  ・・・（２）
【００３３】
式（１）の係数ｖ１および式（２）の係数ｖ２は、予め決められる値であり、例えば、ｖ１＝０．１５、ｖ２＝０．８５である。ｖ１およびｖ２の値は、経験的に決めてもよく、ＰＧ信号のピーク間隔ＰＰintの値に基づいて決めてもよい。
【００３４】
データ収集シーケンスＤＡＱを心拡張期ＤＰに実行するためには、データ収集シーケンスＤＡＱは、ＰＧ信号のピークＰｎ＋ｋから遅延時間ＴＤ２だけ遅延させて実行すればよい。遅延時間ＴＤ２は、ＰＧ信号のピーク間隔ＰＰintの値に基づいて決めており、以下の式で表される。
ＴＤ２＝ｕ・ＰＰint  ・・・（３）
ただし、ｖ１≦ｕ＜ｖ２
【００３５】
係数ｕは、予め決められる値であり、例えば、ｕ＝０．２である。ｕの値は、経験的に決めてもよく、ＰＧ信号のピーク間隔ＰＰintの値に基づいて決めてもよい。ｕ＝０．２の場合、ＴＤ２＝０．２ＰＰintとなる。ピークＰｎ＋ｋの直前および直後は心収縮期ＳＰであるが、ピークＰｎ＋ｋから０．２ＰＰintだけ遅れた時刻では、心拡張期ＤＰに十分に移行していると考えられるので、データ収集シーケンスＤＡＱを心拡張期ＤＰに開始することができる。ただし、係数ｕの値を係数ｖ２の値に近づけすぎると、ピーク間隔ＰＰintの値によっては、データ収集シーケンスＤＡＱが心拡張期ＤＰからはみ出てしまうので、係数ｕの値は、できるだけ係数ｖ１の値に近いことが望ましい（尚、ここでは、データ収集シーケンスＤＡＱに掛かる時間ＴＡは、心拡張期ＤＰよりもある程度短いとする）。したがって、ＰＧ信号のピークＰｎ＋ｋを、データ収集シーケンスＤＡＱを実行するためのトリガポイントにすることができれば、ピークＰｎ＋ｋから遅延時間ＴＤ２だけ遅れてデータ収集シーケンスＤＡＱを開始させることによって、データを心拡張期ＤＰに収集できる。
【００３６】
しかし、パルスシーケンスＰＳは、データ収集シーケンスＤＡＱの前に、選択反転パルスＳＩＲを有している。したがって、データ収集シーケンスＤＡＱの直前に発生するＰＧ信号のピークＰｎ＋ｋをトリガポイントとすることはできず、データ収集シーケンスＤＡＱの直前に発生するＰＧ信号のピークＰｎ＋ｋをトリガポイントにしなければならない。このため、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行するには、データ収集シーケンスＤＡＱがＰＧ信号のピークＰｎ＋ｋから遅延時間ＴＤ２だけ遅れて実行されるように、パルスシーケンスＰＳの遅延時間ＴＤ１を決定する必要がある。以下に、パルスシーケンスＰＳの遅延時間ＴＤ１の決定方法について説明する。
【００３７】
図４から、遅延時間ＴＤ１、反転時間ＴＩ、ピークＰｎとピークＰｎ＋ｋとの間の間隔Ｌ、および遅延時間ＴＤ２の間には、以下の関係式が成り立つ。
ＴＤ１＋ＴＩ＝Ｌ＋ＴＤ２  ・・・（４）
【００３８】
ここで、ピーク間隔ＰＰintは、パルスシーケンスＰＳが実行されている間は変化しない一定値であると仮定する。この場合、ピークＰｎとピークＰｎ＋ｋとの間の間隔Ｌは、以下の式のように、ピーク間隔ＰＰintの整数倍で表される。
Ｌ＝ｋ・ＰＰint  ・・・（５）
ここで、ｋ：整数
【００３９】
式（５）を式（４）に代入すると、以下の式が得られる。
ＴＤ１＋ＴＩ＝ｋ・ＰＰint＋ＴＤ２  ・・・（６）
【００４０】
式（６）から、遅延時間ＴＤ１は、以下の式で表される。
ＴＤ１＝ｋ・ＰＰint＋ＴＤ２−ＴＩ  ・・・（７）
【００４１】
式（３）を式（７）に代入すると、以下の式が得られる。
ＴＤ１＝ｋ・ＰＰint＋ｕ・ＰＰint−ＴＩ  ・・・（８）
【００４２】
また、ピーク間隔ＰＰintは、心拍数ＨＲを用いて、以下の式で表される。
ＰＰint＝６０／ＨＲ  ・・・（９）
【００４３】
したがって、式（９）を式（８）に代入すると、以下の式が得られる。
ＴＤ１＝ｋ・（６０／ＨＲ）＋ｕ・（６０／ＨＲ）−ＴＩ  ・・・（１０）
【００４４】
式（１０）において、反転時間ＴＩは、パルスシーケンスＰＳによって決まる固定値であり、心拍数ＨＲは、ＰＧ信号から求められる値である。また、係数ｕは、式（３）を参照しながら説明したように、予め決められる値である。したがって、係数ｋの値を求めることができれば、遅延時間ＴＤ１を決定することができる。以下に、係数ｋの算出方法について説明する。
【００４５】
係数ｋを算出する場合、先ず、反転時間ＴＩの値と遅延時間ＴＤ２の値とを比較する。反転時間ＴＩはパルスシーケンスＰＳによって決まる固定値であり、遅延時間ＴＤ２は心拍数ＨＲから求められる値であるので（式（３）および（９）参照）、反転時間ＴＩの値および遅延時間ＴＤ２の値は知ることができる。反転時間ＴＩの値と遅延時間ＴＤ２の値とを比較すると、以下の場合に分けられる。
ＴＩ≧ＴＤ２、又はＴＩ＜ＴＤ２  ・・・（１１）
【００４６】
以下に、係数ｋを算出する手順について、ＴＩ≧ＴＤ２の場合と、ＴＩ＜ＴＤ２の場合とに分けて考察する。
（１）  ＴＩ≧ＴＤ２の場合
先ず、図４から、以下の式が成り立つことがわかる。
ＴＩ−ＴＤ２＜Ｌ＜ＴＩ−ＴＤ２＋ＰＰint  ・・・（１２）
【００４７】
ここで、遅延時間ＴＤ２は式（３）で表され、間隔Ｌは式（５）で表され、ピーク間隔ＰＰintは式（９）で表されるので、式（１２）は、以下の式に変形できる。
ＴＩ−ｕ・（60/HR）＜ｋ・（60/HR）＜ＴＩ−（ｕ−１）・（60/HR）  ・・・（１３）
【００４８】
式（１３）において、反転時間ＴＩはパルスシーケンスＰＳによって決まる固定値であり、係数ｕは予め決められる値であり、心拍数ＨＲはＰＧ信号から計測される値である。したがって、式（１３）が成り立つ整数ｋを求めることができる。例えば、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝６０の場合、係数ｋは、以下のようになる。
【００４９】
図５は、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝６０の場合の係数ｋの値について説明する図である。
【００５０】
ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝６０の場合、式（１３）は、以下のようになる。
１．３＜ｋ＜２．３
【００５１】
したがって、式（１３）を満たす整数ｋは、ｋ＝２となることがわかる。ｋ＝２の場合、ＰＧ信号のピークＰｎとＰｎ＋２との間の間隔Ｌは、Ｌ＝２・ＰＰintとなる。
【００５２】
ｋ＝２、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝６０を式（１０）に代入すると、ＴＤ１＝０．７（sec）となる。したがって、ＰＧ信号のピークＰｎから０．７（sec）だけ遅延させてパルスシーケンスＰＳを実行させることによって、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。
【００５３】
また、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝４０の場合、係数ｋは、以下のようになる。
【００５４】
図６は、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝４０の場合の係数ｋの値について説明する図である。
【００５５】
ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝４０の場合、式（１３）は、以下のようになる。
１．２＜１．５ｋ＜２．７
【００５６】
したがって、式（１３）を満たす整数ｋは、ｋ＝１となることがわかる。ｋ＝１の場合、ＰＧ信号のピークＰｎとＰｎ＋１との間の間隔Ｌは、Ｌ＝ＰＰintとなる。
【００５７】
ｋ＝１、ＴＩ＝１．５（sec）、ｕ＝０．２、ＨＲ＝４０を式（１０）に代入すると、ＴＤ１＝０．３（sec）となる。したがって、ＰＧ信号のピークＰｎから０．３（sec）だけ遅延させてパルスシーケンスＰＳを実行させることによって、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。
【００５８】
（２）  ＴＩ＜ＴＤ２の場合
ＴＩ＜ＴＤ２の場合は、図４において、ＰＧ信号のピークＰｎとＰｎ＋ｋが一致していると見なすことができる。したがって、ピークＰｎとＰｎ＋ｋとの間の間隔Ｌ＝０となるので、ｋ＝０となる。係数ｋ＝０の場合、ＰＧ信号とパルスシーケンスＰＳとの関係は、以下のようになる。
【００５９】
図７は、係数ｋ＝０の場合のＰＧ信号とパルスシーケンスＰＳとの関係について説明する図である。
【００６０】
ｋ＝０の場合、ＰＧ信号のピーク間隔ＰＰintの間に、パルスシーケンスＰＳが実行されることになる。ｋ＝０を式（１０）に代入すると、以下の式が得られる。
ＴＤ１＝ｕ・（６０／ＨＲ）−ＴＩ  ・・・（１０）’
【００６１】
上記のように、反転時間ＴＩは、パルスシーケンスＰＳによって決まる固定値であり、心拍数ＨＲはＰＧ信号から求められる値であり、係数ｕは予め決められる値である。したがって、式（１０）’を用いて、遅延時間ＴＤ１を決定することができる。
【００６２】
以上説明したように、係数ｋを算出することができるので、式（１２）を用いることいよって、遅延時間ＴＤ１が求められることがわかる。
【００６３】
しかし、撮影中に被検体１４の心拍が変動することがある。
【００６４】
図８は、撮影中に被検体１４の心拍が変動した様子を示す図である。
図８（ａ）は、心拍数ＨＲ＝ａのときのＰＧ信号を示す図、図８（ｂ）は、心拍数ＨＲがａからｂに変化した場合のＰＧ信号を示す図、図８（ｃ）は、ｎ回目のパルスシーケンスＰＳと、ｎ＋１回目のパルスシーケンスＰＳとを示す図である。
【００６５】
図８を参照すると、心拍数ＨＲ＝ａの場合（図８（ａ）参照）、遅延時間ＴＤ１＝ｔｄ１１でパルスシーケンスＰＳを実行することによって、データ収集シーケンスＤＡＱは心拡張期ＤＰａに実行される。しかし、心拍数ＨＲがＨＲ＝ｂに変化すると（図８（ｂ）参照）、データ収集シーケンスＤＡＱの大部分が心収縮期ＳＰｂに実行されてしまうことになる。したがって、心拍が変動したにもかかわらず、遅延時間ＴＤ１＝ｔｄ１１でパルスシーケンスＰＳを実行し続けてしまうと、動脈を十分に強調して描出することができない恐れがある。そこで、第１の実施形態のＭＲＩ装置１は、撮影中に心拍が変動した場合、遅延時間ＴＤ１の値を更新し、パルスシーケンスＰＳを実行している。
【００６６】
以下に、パルスシーケンスＰＳを実行するためのＭＲＩ装置１の処理フローについて、図９〜図１２を参照しながら説明する。
【００６７】
図９は、ＭＲＩ装置１の処理フローを示す図、図１０〜図１２は、図９の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【００６８】
図９の処理フローの説明に当たっては、必要に応じて図１０〜図１２を参照しながら説明する。
【００６９】
先ず、ステップＳ１において、被検体１４の心拍数ＨＲを算出する。ここでは、時刻ｔ１における被検体１４の心拍数ＨＲを算出する（図１０参照）。算出された時刻ｔ１における心拍数ＨＲは、ＨＲ＝ａであるとする。心拍数ＨＲ＝ａを算出した後、ステップＳ２に進む。
【００７０】
ステップＳ２では、図４〜図７を参照しながら説明した手順で、心拍数ＨＲ＝ａのときの係数ｋを算出する。ここでは、係数ｋ＝１であるとする（したがって、ＰＧ信号とパルスシーケンスＰＳとの関係は、図６で表される）。係数ｋを算出した後、ステップＳ３に進む。
【００７１】
ステップＳ３では、算出した心拍数ＨＲ＝ａおよび係数ｋ＝１を式（１０）に代入し、遅延時間ＴＤ１を計算する。式（１０）において、ｕは予め決められている値であり、反転時間ＴＩもパルスシーケンスＰＳに応じて決まっている固定値であるので、算出した心拍数ＨＲ＝ａおよび係数ｋ＝１を式（１０）に代入することにより、遅延時間ＴＤ１を決定することができる。ここでは、決定された遅延時間ＴＤ１＝ｔｄ１１とする。遅延時間ＴＤ１を決定した後、ステップＳ４に進む。
【００７２】
ステップＳ４では、パルスシーケンスＰＳを実行する。ここでは、先ず、１回目のパルスシーケンスＰＳが実行される。１回目のパルスシーケンスＰＳは、ＰＧ信号のピークＰ１１から遅延時間ｔｄ１１だけ遅れて実行される。遅延時間ｔｄ１１は、ステップＳ３において式（１０）を用いて決定された値である。したがって、データ収集シーケンスＤＡＱを、ＰＧ信号のピークＰ１２から遅延時間ＴＤ２だけ遅れて開始させることができるので、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。１回目のパルスシーケンスＰＳを実行した後、ステップＳ５に進む。
【００７３】
ステップＳ５では、ｋ空間を埋めるために必要な全データが収集されたか否かを判断する。全データが収集された場合、フローを終了する。一方、全データが収集されていない場合は、ステップＳ６に進む。ここでは、全データが収集されていないと判断され、ステップＳ６に進む。
【００７４】
ステップＳ６では、時刻ｔｎにおける心拍数ＨＲを算出する。ここでは、時刻ｔ２における心拍数ＨＲを算出する。算出された時刻ｔ２における心拍数ＨＲは、ＨＲ＝ａであるとする。心拍数ＨＲ＝ａを算出した後、ステップＳ７に進む。
【００７５】
ステップＳ７では、時刻ｔｎにおける心拍数ＨＲと、時刻ｔｎ−１における心拍数ＨＲとが同じであるか否かを求めることによって、心拍が変動したか否かを判断する。時刻ｔｎにおける心拍数ＨＲが、時刻ｔｎ−１に算出された心拍数ＨＲと同じである場合、心拍は変動していないと判断される。この場合、遅延時間ＴＤ１を変更することなく、ステップＳ４に戻る。一方、時刻ｔｎにおける心拍数ＨＲが、時刻ｔｎ−１に算出された心拍数ＨＲとは異なる場合、心拍が変動したと判断され、ステップＳ８に進む。ここでは、時刻ｔ２における心拍数ＨＲが、時刻ｔ１における心拍数ＨＲと同じであるか否かが判断される。時刻ｔ２における心拍数ＨＲ、および時刻ｔ１における心拍数ＨＲは、いずれもＨＲ＝ａであるので、心拍数ＨＲは同じである。したがって、心拍は変動していないと判断され、遅延時間ＴＤ１を変更することなくステップＳ４に戻る。
【００７６】
ステップＳ４では、パルスシーケンスＰＳを実行する。ここでは、２回目のパルスシーケンスＰＳが実行される。遅延時間ＴＤ１は変更されていないので、２回目のパルスシーケンスＰＳは、ＰＧ信号のピークＰ２１から遅延時間ｔｄ１１だけ遅れて実行される。したがって、２回目のパルスシーケンスＰＳにおいても、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。２回目のパルスシーケンスＰＳを実行した後、ステップＳ５に進む。
【００７７】
ステップＳ５では、ｋ空間を埋めるために必要な全データが収集されたか否かを判断する。ここでは、全データが収集されていないと判断され、ステップＳ６に進む。
【００７８】
ステップＳ６では、時刻ｔｎにおける心拍数ＨＲを算出する。ここでは、時刻ｔ３における心拍数ＨＲを算出する（図１１参照）。算出された時刻ｔ３における心拍数ＨＲは、ＨＲ＝ａであるとする。心拍数ＨＲ＝ａを算出した後、ステップＳ７に進む。
【００７９】
ステップＳ７では、時刻ｔｎにおける心拍数ＨＲと、時刻ｔｎ−１における心拍数ＨＲとが同じであるか否かを求めることによって、心拍が変動したか否かを判断する。ここでは、時刻ｔ３における心拍数ＨＲが、時刻ｔ２における心拍数ＨＲ（図１０参照）と同じであるか否かが判断される。時刻ｔ３における心拍数ＨＲ、および時刻ｔ２における心拍数ＨＲは、いずれもＨＲ＝ａであるので、心拍数ＨＲは同じである。したがって、心拍は変動していないと判断され、遅延時間ＴＤ１を変更することなくステップＳ４に戻る。
【００８０】
ステップＳ４では、パルスシーケンスＰＳを実行する。ここでは、３回目のパルスシーケンスＰＳが実行される（図１１参照）。遅延時間ＴＤ１は変更されていないので、３回目のパルスシーケンスＰＳは、ＰＧ信号のピークＰ３１から遅延時間ｔｄ１１だけ遅れて実行される。したがって、３回目のパルスシーケンスＰＳにおいても、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。３回目のパルスシーケンスＰＳを実行した後、ステップＳ５に進む。
【００８１】
ステップＳ５では、ｋ空間を埋めるために必要な全データが収集されたか否かを判断する。ここでは、全データが収集されていないと判断され、ステップＳ６に進む。
【００８２】
ステップＳ６では、時刻ｔｎにおける心拍数ＨＲを算出する。ここでは、時刻ｔ４における心拍数ＨＲを算出する。図１１では、時刻ｔ４の前の時刻ｔｊにおいて、心拍数ＨＲが、ＨＲ＝ａからＨＲ＝ｂに変化したとする。したがって、時刻ｔ４における心拍数ＨＲは、ＨＲ＝ｂと算出される。心拍数ＨＲ＝ｂを算出した後、ステップＳ７に進む。
【００８３】
ステップＳ７では、時刻ｔｎにおける心拍数ＨＲと、時刻ｔｎ−１における心拍数ＨＲとが同じであるか否かを求めることによって、心拍が変動したか否かを判断する。ここでは、時刻ｔ４における心拍数ＨＲが、時刻ｔ３における心拍数ＨＲと同じであるか否かが判断される。時刻ｔ４における心拍数ＨＲはＨＲ＝ｂであるが、時刻ｔ３における心拍数ＨＲはＨＲ＝ａである。したがって、心拍数ＨＲは異なるので、心拍は変動したと判断され、ステップＳ８に進む。
【００８４】
ステップＳ８では、図４〜図７を参照しながら説明した手順で、心拍数ＨＲ＝ｂのときの係数ｋを算出する。ここでは、係数ｋ＝１であるとする。係数ｋを算出した後、ステップＳ９に進む。
【００８５】
ステップＳ９では、ステップＳ６で算出した心拍数ＨＲ＝ｂと、ステップＳ８で算出した係数ｋ＝１を式（１０）に代入し、遅延時間ＴＤ１を計算する。これにより、新たな遅延時間ＴＤ１が決定される。ここでは、決定された新たな遅延時間ＴＤ１＝ｔｄ１２とする。新たな遅延時間ＴＤ１を決定した後、ステップＳ１０に進む。
【００８６】
ステップＳ１０では、遅延時間ＴＤ１を、ステップＳ９において決定された新たな遅延時間ＴＤ１に更新する。ステップＳ９では、新たな遅延時間ＴＤ１＝ｔｄ１２であるので、遅延時間ＴＤ１はＴＤ１＝ｔｄ１２に更新される。遅延時間ＴＤ１を更新したら、ステップＳ４に戻る。
【００８７】
ステップＳ４では、パルスシーケンスＰＳを実行する。ここでは、４回目のパルスシーケンスＰＳが実行される。遅延時間ＴＤ１は更新されているので、４回目のパルスシーケンスＰＳは、ＰＧ信号のピークＰ４１から新たな遅延時間ｔｄ１２だけ遅れて実行される。遅延時間ｔｄ１２は、ステップＳ１０において式（１０）を用いて決定された値である。したがって、データ収集シーケンスＤＡＱを、ＰＧ信号のピークＰ４２から遅延時間ＴＤ２だけ遅れて開始させることができるので、データ収集シーケンスＤＡＱを心拡張期ＤＰに実行することができる。
【００８８】
以下同様に、ステップＳ５において全データが収集されたと判断されるまで、図９に示すフローが繰り返し実行される。
【００８９】
図１２は、ＰＧ信号と、ｎ回目およびｎ＋１回目のパルスシーケンスＰＳとを示す図である。
【００９０】
図１２では、ｎ回目のパルスシーケンスＰＳが実行される前は、遅延時間ＴＤ１＝ｔｄ１１、心拍数ＨＲ＝ａであるが、ｎ＋１回目のパルスシーケンスＰＳが実行される前に、心拍数ＨＲ＝ｃに変化している。しかし、時刻ｔｎ＋１において、心拍数ＨＲ＝ｃと算出されるので、心拍数ＨＲ＝ｃに対応した遅延時間ＴＤ１＝ｔｄ１３が決定される。したがって、ｎ＋１回目のパルスシーケンスＰＳは、心拡張期ＤＰに実行される。
【００９１】
以上説明したように、第１の実施形態では、パルスシーケンスＰＳを実行する前に心拍数ＨＲを算出し、心拍が変動した場合、遅延時間ＴＤ１を更新する。したがって、撮影中に被検体１４の心拍が変動しても、心拡張期でデータを収集することができ、高品質なＭＲ画像を得ることができる。
【００９２】
尚、図１２では、心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃが、心拍数ｂのときのＰＧ信号のピーク間隔ＰＰint＝ｗａよりも長い場合について示されている。しかし、心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃは、心拍数ａのときのＰＧ信号のピーク間隔ＰＰint＝ｗａより短くなる場合もある。以下に、心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃが、心拍数ａのときのＰＧ信号のピーク間隔ＰＰint＝ｗａより短くなった場合について説明する。
【００９３】
図１３は、心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃが、心拍数ａのときのＰＧ信号のピーク間隔ＰＰint＝ｗａより短くなった場合の説明図である。
【００９４】
心拍数ＨＲは変動値であるが、パルスシーケンスＰＳのデータ収集シーケンスＤＡＱの時間ＴＡは固定値である。したがって、心拍数ｃのときのＰＧ信号のピーク間隔ＰＰint＝ｗｃが狭くなりすぎると、データ収集シーケンスＤＡＱの一部ＤＡＱｓは心拡張期ＤＰに含まれているが、残りの部分ＤＡＱｒは心収縮期ＳＰに含まれてしまうことがある。しかし、このような場合であっても、データ収集シーケンスＤＡＱの一部ＤＡＱｓにおいて、ｋ空間の低周波領域のデータを収集しておけば、十分なコントラストを有するＭＲ画像を得ることができる。したがって、ピーク間隔ＰＰintが狭くなっても、式（１０）を用いて遅延時間ＴＤ１を決定することによって、高品質なＭＲ画像を得ることができる。
【００９５】
尚、第１の実施形態では、式（１０）を用いて遅延時間ＴＤ１を決定している。しかし、式（１０）の代わりに、式（８）を用いて遅延時間ＴＤ１を決定してもよい。
【００９６】
第１の実施形態では、脈波信号を用いて被検体１４の心拍数ＨＲを算出している。しかし、脈波信号の代わりに、心電信号を用いて被検体１４の心拍数ＨＲを算出してもよい。
【００９７】
第１の実施形態では、パルスシーケンスＰＳは、データ収集シーケンスＤＡＱの前に、選択反転パルスＳＩＲを有しているが、本発明は、パルスシーケンスＰＳが選択反転パルスＳＩＲを備えておらず、データ収集シーケンスＤＡＱしか備えていない場合にも適用できる。パルスシーケンスＰＳがデータ収集シーケンスＤＡＱしか備えていない場合、遅延時間ＴＤ１は、ＰＧ信号のピークと、データ収集シーケンスＤＡＱの一番最初に送信されるＲＦパルスとの間の時間になる。また、パルスシーケンスＰＳは、必要に応じて、データ収集シーケンスＤＡＱの後に、一つ以上のＲＦパルスを有していてもよい。
【００９８】
また、第１の実施形態では、パルスシーケンスＰＳは、データ収集シーケンスＤＡＱの前に、選択反転パルスＳＩＲを有している。しかし、パルスシーケンスＰＳは、必要に応じて、選択反転パルスＳＩＲの代わりに、別のＲＦパルス（例えば、非選択反転パルス、フリップ角αがα≠９０°のＲＦパルス）を有していてもよい。また、パルスシーケンスＰＳは、必要に応じて、選択反転パルスＳＩＲとデータ収集シーケンスＤＡＱとの間に、一つ以上の別のＲＦパルスを有していてもよい。
【００９９】
更に、第１の実施形態では、パルスシーケンスＰＳは、データ収集シーケンスＤＡＱの前に、選択反転パルスＳＩＲのみを有しているが、パルスシーケンスＰＳは、必要に応じて、データ収集シーケンスＤＡＱの前に、複数のＲＦパルスを有していてもよい。この場合、遅延時間ＴＤ１は、ＰＧ信号のピークと、複数のＲＦパルスのうちの一番最初に送信されるＲＦパルスとの間の時間となる。
（２）第２の実施形態
第２の実施形態のＭＲＩ装置のハードウェア構成は、第１の実施形態のＭＲＩ装置と同じであるので、ハードウェア構成の説明は省略し、相違点を主に説明する。
【０１００】
図１４は、被検体１４に不整脈が発生した場合のＰＧ信号とパルスシーケンスＰＳとの関係を示す図である。
【０１０１】
図１４（ａ）は、被検体１４に不整脈が発生していない場合のＰＧ信号を示す図、図１４（ｂ）は、被検体１４に不整脈が発生した場合のＰＧ信号を示す図、図１４（ｃ）は、ｎ回目のパルスシーケンスＰＳおよびｎ＋１回目のパルスシーケンスを示す図である。
【０１０２】
図１４（ａ）に示すように、心拍数ＨＲ＝ａの場合、被検体１４に不整脈が発生しなければ、ＰＧ信号のピーク間隔ＰＰintは、ＰＰint＝ｗａで一定である。したがって、遅延時間ＴＤ１＝ｔｄ１１でｎ回目のパルスシーケンスＰＳを実行した場合、データ収集シーケンスＤＡＱは、心拡張期ＤＰに実行される。
【０１０３】
しかし、図１４（ｂ）に示すように、心拍数ＨＲ＝ａであっても、被検体１４に不整脈が発生した場合、ＰＧ信号のピーク間隔ＰＰintが、瞬間的にＰＰint＝ｗａ’に変化する。したがって、ＰＧ信号のピーク間隔ＰＰintがｗａ’からｗａに戻っても、データ収集シーケンスＤＡＱの一部ＤＡＱｓが心収縮期ＳＰに実行されてしまうことになる。そこで、第２の実施形態では、データ収集シーケンスＤＡＱの一部ＤＡＱｓが心収縮期ＳＰに実行される可能性が高い場合、次のパルスシーケンスＰＳでデータを取り直している。以下に、データを取り直すときのフローについて説明する。
【０１０４】
図１５は、データを取り直すときのフローを説明する図、図１６は、図１５の処理フローが実行されているときのＰＧ信号およびパルスシーケンスＰＳを示す図である。
【０１０５】
尚、図１５の説明に当たっては、図１６を参照しながら説明する。
【０１０６】
先ず、ステップＳ１において、ｎ回目のパルスシーケンスＰＳのデータ収集シーケンスＤＡＱの直前の時刻ｔｘにおける心位相Ｈが、心拡張期ＤＰであるのか心収縮期ＳＰであるのかを求める。心拡張期ＤＰの範囲は、図４を参照しながら説明したように、時間Ｔｓ（式（１）参照）と、時間Ｔｅ（式（２）参照）で求めることができる。したがって、ＰＧ信号のピークＰｎ３から時刻ｔｘまでの時間差をΔｔとすると、以下の式（１４）が成り立つ場合、時刻ｔｘにおける心位相Ｈは、心拡張期ＤＰであると考えられる。
Ｔｓ≦Δｔ≦Ｔｅ  ・・・（１４）
【０１０７】
一方、以下の式（１５）が成り立つ場合、時刻ｔｘにおける心位相Ｈは、心収縮期ＳＰであると考えられる。
Δｔ＜Ｔｓ  又は  Δｔ＞Ｔｅ  ・・・（１５）
【０１０８】
図１６を参照すると、Δｔ＞Ｔｅである。したがって、式（１５）が成り立つので、心位相Ｈは、心収縮期ＳＰであると判断される。心位相Ｈを求めた後、ステップＳ２に進む。
【０１０９】
ステップＳ２では、心位相Ｈに基づいてデータを取り直すか否かを決定する。時刻ｔｘにおける心位相Ｈが心収縮期ＳＰであるということは、時刻ｔｘの直後に開始されるデータ収集シーケンスＤＡＱの始端側の部分ＤＡＱｓが心収縮期ＳＰの間に実行される可能性が高いことを意味する。したがって、心位相Ｈが心収縮期ＳＰの場合、ｎ回目のパルスシーケンスＰＳによって収集されたデータは取り直す必要があると決定し、フローを終了する。
【０１１０】
ｎ回目のパルスシーケンスＰＳによって収集されたデータは取り直す必要があると決定した場合、ｎ＋１回目のパルスシーケンスＰＳは、以下のように行われる。
【０１１１】
先ず、時刻ｔｎ＋１において、心拍数ＨＲが算出される。ここでは、時刻ｔｎ＋１における心拍数ＨＲ＝ｃであるとする。したがって、心拍数ＨＲは、ａからｃに変化しているので、新たな遅延時間ＴＤ１＝ｔｄ１３が決定される。ｎ＋１回目のパルスシーケンスＰＳは、新たな遅延時間ｔｄ１３で実行される。ただし、ｎ回目のパルスシーケンスＰＳによって収集されたデータは取り直す必要があると決定されているので、ｎ＋１回目のパルスシーケンスＰＳでは、ｎ回目のパルスシーケンスＰＳによって収集されたデータが取り直される。
【０１１２】
ｎ＋１回目のパルスシーケンスＰＳにおいても、図１５に示すフローが実行される。ｎ＋１回目のパルスシーケンスＰＳにおいては、ステップＳ１において、データ収集シーケンスＤＡＱの直前の時刻ｔｙにおける心位相Ｈが、心拡張期ＤＰであるのか心収縮期ＳＰであるのかが求められる。ｎ＋１回目のパルスシーケンスＰＳにおいては、Ｔｓ≦Δｔ≦Ｔｅである。したがって、式（１４）が成り立つので、心位相Ｈは、心拡張期ＤＰであると判断される。心位相Ｈを求めた後、ステップＳ２に進む。
【０１１３】
ステップＳ２では、心位相Ｈに基づいてデータを取り直すか否かを決定する。時刻ｔｙにおける心位相Ｈが心拡張期ＤＰであるということは、時刻ｔｙの直後に開始されるデータ収集シーケンスＤＡＱの大部分は、心拡張期ＤＰの間に実行される可能性が高いことを意味する。したがって、心位相Ｈが心拡張期ＤＰの場合、ｎ＋１回目のパルスシーケンスＰＳによって収集されたデータは取り直す必要がないと決定し、フローを終了する。
【０１１４】
以下同様に、図１５に示すフローを実行しながら、パルスシーケンスＰＳが繰り返し実行される。
【０１１５】
第２の実施形態では、データ収集シーケンスＤＡＱの開始直前の心位相Ｈが心収縮期ＳＰの場合は、データを取り直すことができるので、より高画質なＭＲ画像が得られる。
【０１１６】
尚、上記の説明では、データ収集シーケンスＤＡＱの開始直前の心位相Ｈが心収縮期ＳＰであるのか心拡張期ＤＰであるのかを求めている。しかし、データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’が心収縮期ＳＰであるのか心拡張期ＤＰであるのかを求めてもよい。ただし、データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’は、データ収集シーケンスＤＡＱの開始直前における心位相Ｈと異なる場合がある。例えば、図１６のｎ回目のパルスシーケンスＰＳでは、データ収集シーケンスＤＡＱの開始直前の時刻ｔｘにおける心位相Ｈは心収縮期ＳＰと判断され、一方、データ収集シーケンスＤＡＱの終了直後の時刻ｔｘ’における心位相Ｈ’は心拡張期ＤＰと判断される。しかし、心位相が心収縮期ＳＰであると判断されるのか、それとも心拡張期ＤＰであると判断されるのかに応じて、データが取り直されるか否かが決まるので、心位相の判断を、データ収集シーケンスＤＡＱの開始直前に行うか終了直後に行うかは重要である。以下に、心位相の判断を、データ収集シーケンスＤＡＱの開始直前に行う方が好ましい場合と、終了直後に行う方が好ましい場合の一例について説明する。
【０１１７】
ｋ空間の高周波領域を埋めるデータは、低周波領域を埋めるデータに比べて、画像のコントラストに与える影響は小さい。したがって、ｋ空間の高周波領域を埋めるデータが、心収縮期ＳＰの間に収集されても、画像のコントラストが極端に低下することはない。しかし、ｋ空間の低周波領域を埋めるデータが、心収縮期ＳＰの間に収集されてしまうと、画像のコントラストが大きく低下することがある。したがって、ｋ空間の低周波領域を埋めるデータが、心収縮期ＳＰの間に収集された場合は、データの取直しが実行されるようにすることが望ましい。そこで、ｋ空間の低周波領域を埋めるデータが、データ収集シーケンスＤＡＱの最初の方で収集される場合は、データ収集シーケンスＤＡＱの開始直前における心位相Ｈを求めることが好ましい。例えば、ｎ回目のパルスシーケンスＰＳを参照すると、データ収集シーケンスＤＡＱの始端側の部分ＤＡＱｓは、心収縮期ＳＰの間に実行されている。したがって、ｋ空間の低周波領域を埋めるデータが、データ収集シーケンスＤＡＱの最初の方で収集される場合、ｎ回目のパルスシーケンスＰＳでは、ｋ空間の低周波領域を埋めるデータは、心収縮期ＳＰに収集されることになる。しかし、データ収集シーケンスＤＡＱの開始直前における心位相Ｈを求めていれば、心位相Ｈは心収縮期ＳＰと判断されるので、ｎ＋１回目のパルスシーケンスＰＳで、データが取り直される。したがって、コントラストの低下が抑制された高品質な画像を得ることができる。
【０１１８】
一方、ｋ空間の低周波領域を埋めるデータが、データ収集シーケンスＤＡＱの最後の方で収集される場合は、データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’を求めることが好ましい（図１７参照）。
【０１１９】
図１７は、データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’を求めることが好ましい場合の説明図である。
【０１２０】
図１７において、ｎ回目のパルスシーケンスＰＳを参照すると、データ収集シーケンスＤＡＱの終端側の部分ＤＡＱｅは、心収縮期ＳＰの間に実行されている。したがって、ｋ空間の低周波領域を埋めるデータが、データ収集シーケンスＤＡＱの最後の方で収集される場合、ｎ回目のパルスシーケンスＰＳでは、ｋ空間の低周波領域を埋めるデータは、心収縮期ＳＰに収集されることになる。しかし、データ収集シーケンスＤＡＱの終了直後における心位相Ｈ’を求めていれば、心位相Ｈ’は心収縮期ＳＰと判断されるので、ｎ＋１回目のパルスシーケンスＰＳ（図示せず）で、データが取り直される。したがって、コントラストの低下が抑制された高品質な画像を得ることができる。
【符号の説明】
【０１２１】
１  ＭＲＩ装置
２  磁場発生装置
３  テーブル
４  脈波センサ
５  受信コイル
６  シーケンサ
７  送信器
８  勾配磁場電源
９  ＰＧユニット
１０  受信器
１１  中央処理装置
１２  入力装置
１３  表示装置
１４  被検体
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  送信コイル
３１  クレードル
【図１】
0005601678.tif 000002
【図２】
0005601678.tif 000003
【図３】
0005601678.tif 000004
【図４】
0005601678.tif 000005
【図５】
0005601678.tif 000006
【図６】
0005601678.tif 000007
【図７】
0005601678.tif 000008
【図８】
0005601678.tif 000009
【図９】
0005601678.tif 000010
【図１０】
0005601678.tif 000011
【図１１】
0005601678.tif 000012
【図１２】
0005601678.tif 000013
【図１３】
0005601678.tif 000014
【図１４】
0005601678.tif 000015
【図１５】
0005601678.tif 000016
【図１６】
0005601678.tif 000017
【図１７】
0005601678.tif 000018