JPWO2004080301A1 - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置 Download PDFInfo
- Publication number
- JPWO2004080301A1 JPWO2004080301A1 JP2005503502A JP2005503502A JPWO2004080301A1 JP WO2004080301 A1 JPWO2004080301 A1 JP WO2004080301A1 JP 2005503502 A JP2005503502 A JP 2005503502A JP 2005503502 A JP2005503502 A JP 2005503502A JP WO2004080301 A1 JPWO2004080301 A1 JP WO2004080301A1
- Authority
- JP
- Japan
- Prior art keywords
- navigator
- displacement
- magnetic resonance
- resonance imaging
- body motion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
- G01R33/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Radiology & Medical Imaging (AREA)
- Physiology (AREA)
- Power Engineering (AREA)
- Pulmonology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
少ないナビゲーターエコー取得回数で効果的に体動アーチファクト抑制を行う。被検体(101)の呼息期、吸息期を区別した体動アーチファクト抑制画像を作成する。時間間隔をおいて計測した2以上のナビゲーターエコー(207、208)から、ナビゲーターエコー取得時刻における被検体(101)の変位及び変位方向を求め、これら変位及び変位方向から本計測エコーを取得するタイミングを決定する。変位方向は変位の一次微分から求めることができる。また、計測に先立って体動周期を近似する体動予測関数を求めておくことにより、フィッティングにより変位及び変位方向を判定できるとともに本計測エコー取得の時間、間隔を推定する。
Description
本発明は、核磁気共鳴を利用して被検体を撮像する磁気共鳴イメージング(以下、MRIという)装置に関し、特に体動アーチファクトを抑制する技術に関する。
MRIは、静磁場中に配置された被検体に、撮影対象である原子核スピンの共鳴周波数の高周波磁場を印加し、それによって発生する核磁気共鳴信号を計測することによって被検体の断層像を得る装置である。
このMRIでは、1枚のMR画像を作成する間に被検体が動くと、画像に大きなアーチファクトが生じ、画質を劣化させる。体動アーチファクトの原因として、重要なものが被検体の呼吸による体動アーチファクトである。
呼吸動による体動アーチファクトを除去する方法の一つに、付加的なエコーであるナビゲーターエコーを用いて呼吸動をモニターする方法があり、このようなナビゲーターエコーを用いて呼吸動アーチファクトを除去する技術を心臓撮影に適応した撮影方法が提案されている(文献、”Navigator Echoes in Cardiac Magnetic Resonance,Journal of Cardiovascular Magnetic Resonance”(3),183−193(2001))。
この文献では、ナビゲーターエコーから求めた呼吸による体動位置(変位)が、予め指定された所定の範囲内であるときには、それに続いて実行される撮像シーケンスにおいて取得したデータを画像再構成用のデータとして用いる。しかし、体動の変位方向については考慮されていない。
このMRIでは、1枚のMR画像を作成する間に被検体が動くと、画像に大きなアーチファクトが生じ、画質を劣化させる。体動アーチファクトの原因として、重要なものが被検体の呼吸による体動アーチファクトである。
呼吸動による体動アーチファクトを除去する方法の一つに、付加的なエコーであるナビゲーターエコーを用いて呼吸動をモニターする方法があり、このようなナビゲーターエコーを用いて呼吸動アーチファクトを除去する技術を心臓撮影に適応した撮影方法が提案されている(文献、”Navigator Echoes in Cardiac Magnetic Resonance,Journal of Cardiovascular Magnetic Resonance”(3),183−193(2001))。
この文献では、ナビゲーターエコーから求めた呼吸による体動位置(変位)が、予め指定された所定の範囲内であるときには、それに続いて実行される撮像シーケンスにおいて取得したデータを画像再構成用のデータとして用いる。しかし、体動の変位方向については考慮されていない。
本発明の目的は、呼吸動の時相に依らずに、被検体の呼吸動による体動アーチファクトを抑制可能なMRI装置を実現することである。
上記目的を達成するため、本発明は次のように、構成される。
(1)本発明の磁気共鳴イメージング装置は、所定の撮像シーケンスに基づき、画像再構成に必要なエコー信号を取得して被検体の画像を得る手段と、前記被検体の周期的体動情報を含む核磁気共鳴信号をナビゲーターエコーとして取得するナビゲーターエコー計測手段と、前記被検体の周期的体動に基づく所定部位の変位の許容範囲を設定する許容範囲設定手段と、前記ナビゲーターエコーから前記所定部位の変位を算出して該変位が前記許容範囲内にあるか否かを判定する体動判定手段とを備え、前記画像を得る手段は、前記所定部位の変位が前記許容範囲内にあるときに取得されたエコー信号を用いて画像再構成を行う。
そして、この磁気共鳴イメージング装置において、前記体動判定手段は、時間間隔を置いて計測した2以上のナビゲーターエコーから、前記所定部位の変位方向を算出し、前記所定部位の変位及び変位方向に基づいて、前記画像再構成に必要なエコー信号を取得するタイミングを決定する撮像タイミング決定手段を備える。
(2)好ましくは、上記(1)において、前記ナビゲーターエコー計測手段は、前記被検体の一心周期内に前記2以上のナビゲーターエコーを取得する。
(3)また、好ましくは、上記(2)において、前記ナビゲーターエコー計測手段は、前記一心周期内の前記2以上のナビゲーターの内、1以上を画像再構成用エコー信号の取得期間の前に取得し、他の1以上を該画像再構成用エコー信号の取得期間の後に取得する。
(4)また、好ましくは、上記(1)において、前記ナビゲーターエコー計測手段は、一心周期毎に1以上のナビゲーターエコーを取得する。
(5)また、好ましくは、上記(1)において、3以上のナビゲーターエコーを用いて前記周期的体動を表す体動予測関数を作成する体動予測関数作成手段を備え、前記体動判定手段は、前記2以上のナビゲーターエコーから求めた前記所定部位の変位を前記体動予測関数にフィッティングして、所望の時刻における前記所定部位の変位及び変位方向を算出する。
(6)また、好ましくは、上記(1)において、前記体動判定手段は、前記所定部位の変位の変動に対応して前記許容範囲を調整する。
(7)また、好ましくは、上記(1)において、前記許容範囲設定手段は、前記2以上のナビゲーターエコーの取得タイミングと前記許容範囲を設定する手段を備える。
(8)また、好ましくは、上記(1)において、前記撮像タイミング決定手段は、前記体動判定手段によって求められた前記所定部位の変位が前記許容範囲内であって旦つ前記変位方向が同一方向のときに、前記画像再構成用エコー信号の取得を行うように撮像タイミングを決定する。
(9)また、好ましくは、上記(5)において、前記撮像タイミング決定手段は、前記2以上のナビゲーターエコーの取得時刻及び前記体動予測関数から、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定する。
(10)また、好ましくは、上記(9)において、前記撮像タイミング決定手段が、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定できるまで、前記ナビゲーターエコーの計測を繰り返す。
(11)また、好ましくは、上記(5)において、前記体動予測関数は、少なくとも、被検体の体動の1/4周期に渡るナビゲーターエコーから求めた変位を用いて、被検体毎に作成される。
(12)また、好ましくは、上記(5)において、前記体動予測関数作成手段は、体動決定のために逐次取得されるナビゲーターエコーを用いて前記体動予測関数を更新する。
(13)また、好ましくは、上記(12)において、前記画像再構成用エコー信号の取得期間終了時から次の前記画像再構成用エコー信号の取得開始時までの間に、1以上のナビゲーターエコーを取得し、該ナビゲーターエコーを用いて前記体動予測関数を更新する。
(14)また、好ましくは、上記(13)において、前記体動予測関数作成手段は、前記ナビゲーターエコーを取得する期間と、該体動予測関数の更新に使用するナビゲーターエコーの期間をそれぞれ前記体動予測関数上で指定する手段を備える。
(15)また、好ましくは、上記(2)において、前記被検体の体動が呼吸動の場合に、前記体動判定手段は、一方の前記変位方向を呼息期とし、他方の前記変位方向を吸息期と判定する。
(16)また、好ましくは、上記(1)乃至(4)において、前記体動判定手段は、前記2以上のナビゲーターエコーから算出した前記所定部位の変位の一次微分又は差分に基づいて、該所定部位の変位方向を算出する。
(17)また、好ましくは、上記(1)において、前記画像再構成用エコー信号の取得期間以外はエコー信号の取得を行わずに前記撮像シーケンスを継続する。
(18)また、好ましくは、上記(1)において、2つのナビゲーターエコーの互いの時間間隔は、約100msであり、変位の検出精度は、0.7mm〜0.8mmの分解能である。
上記目的を達成するため、本発明は次のように、構成される。
(1)本発明の磁気共鳴イメージング装置は、所定の撮像シーケンスに基づき、画像再構成に必要なエコー信号を取得して被検体の画像を得る手段と、前記被検体の周期的体動情報を含む核磁気共鳴信号をナビゲーターエコーとして取得するナビゲーターエコー計測手段と、前記被検体の周期的体動に基づく所定部位の変位の許容範囲を設定する許容範囲設定手段と、前記ナビゲーターエコーから前記所定部位の変位を算出して該変位が前記許容範囲内にあるか否かを判定する体動判定手段とを備え、前記画像を得る手段は、前記所定部位の変位が前記許容範囲内にあるときに取得されたエコー信号を用いて画像再構成を行う。
そして、この磁気共鳴イメージング装置において、前記体動判定手段は、時間間隔を置いて計測した2以上のナビゲーターエコーから、前記所定部位の変位方向を算出し、前記所定部位の変位及び変位方向に基づいて、前記画像再構成に必要なエコー信号を取得するタイミングを決定する撮像タイミング決定手段を備える。
(2)好ましくは、上記(1)において、前記ナビゲーターエコー計測手段は、前記被検体の一心周期内に前記2以上のナビゲーターエコーを取得する。
(3)また、好ましくは、上記(2)において、前記ナビゲーターエコー計測手段は、前記一心周期内の前記2以上のナビゲーターの内、1以上を画像再構成用エコー信号の取得期間の前に取得し、他の1以上を該画像再構成用エコー信号の取得期間の後に取得する。
(4)また、好ましくは、上記(1)において、前記ナビゲーターエコー計測手段は、一心周期毎に1以上のナビゲーターエコーを取得する。
(5)また、好ましくは、上記(1)において、3以上のナビゲーターエコーを用いて前記周期的体動を表す体動予測関数を作成する体動予測関数作成手段を備え、前記体動判定手段は、前記2以上のナビゲーターエコーから求めた前記所定部位の変位を前記体動予測関数にフィッティングして、所望の時刻における前記所定部位の変位及び変位方向を算出する。
(6)また、好ましくは、上記(1)において、前記体動判定手段は、前記所定部位の変位の変動に対応して前記許容範囲を調整する。
(7)また、好ましくは、上記(1)において、前記許容範囲設定手段は、前記2以上のナビゲーターエコーの取得タイミングと前記許容範囲を設定する手段を備える。
(8)また、好ましくは、上記(1)において、前記撮像タイミング決定手段は、前記体動判定手段によって求められた前記所定部位の変位が前記許容範囲内であって旦つ前記変位方向が同一方向のときに、前記画像再構成用エコー信号の取得を行うように撮像タイミングを決定する。
(9)また、好ましくは、上記(5)において、前記撮像タイミング決定手段は、前記2以上のナビゲーターエコーの取得時刻及び前記体動予測関数から、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定する。
(10)また、好ましくは、上記(9)において、前記撮像タイミング決定手段が、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定できるまで、前記ナビゲーターエコーの計測を繰り返す。
(11)また、好ましくは、上記(5)において、前記体動予測関数は、少なくとも、被検体の体動の1/4周期に渡るナビゲーターエコーから求めた変位を用いて、被検体毎に作成される。
(12)また、好ましくは、上記(5)において、前記体動予測関数作成手段は、体動決定のために逐次取得されるナビゲーターエコーを用いて前記体動予測関数を更新する。
(13)また、好ましくは、上記(12)において、前記画像再構成用エコー信号の取得期間終了時から次の前記画像再構成用エコー信号の取得開始時までの間に、1以上のナビゲーターエコーを取得し、該ナビゲーターエコーを用いて前記体動予測関数を更新する。
(14)また、好ましくは、上記(13)において、前記体動予測関数作成手段は、前記ナビゲーターエコーを取得する期間と、該体動予測関数の更新に使用するナビゲーターエコーの期間をそれぞれ前記体動予測関数上で指定する手段を備える。
(15)また、好ましくは、上記(2)において、前記被検体の体動が呼吸動の場合に、前記体動判定手段は、一方の前記変位方向を呼息期とし、他方の前記変位方向を吸息期と判定する。
(16)また、好ましくは、上記(1)乃至(4)において、前記体動判定手段は、前記2以上のナビゲーターエコーから算出した前記所定部位の変位の一次微分又は差分に基づいて、該所定部位の変位方向を算出する。
(17)また、好ましくは、上記(1)において、前記画像再構成用エコー信号の取得期間以外はエコー信号の取得を行わずに前記撮像シーケンスを継続する。
(18)また、好ましくは、上記(1)において、2つのナビゲーターエコーの互いの時間間隔は、約100msであり、変位の検出精度は、0.7mm〜0.8mmの分解能である。
図1は、本発明が適用されるMRI装置の全体概要を示す図である。
図2は、本発明の第1の態様による撮像方法の一例を示す図である。
図3は、本発明の第1の態様による撮像方法の効果を説明する図である。
図4は、本発明の第1の態様による撮像方法の他の例を示す図である。
図5は、本発明の第1の態様による撮像方法の他の例を示す図である。
図6は、本発明の第2の態様による撮像方法の手順を示す図である。
図7は、図6に示す撮像方法による体動予測関数の作成を説明する図である。
図8は、図6に示す撮像方法における本計測を説明する図である。
図9は、関数フィッティング(第2のフィッティング)を説明する図である。
図10は、本発明の第2の態様による撮像方法の一例を示す図である。
図11は、本発明の第2の態様による撮像方法の他の例を示す図である。
図12は、本発明の第2の態様が適用される撮像方法を示す図である。
図13は、図12に示す撮像方法への本発明の適用を説明する図である。
図14は、ナビゲーターエコー取得タイミングを指定するためのUIの一例を示す図である。
図2は、本発明の第1の態様による撮像方法の一例を示す図である。
図3は、本発明の第1の態様による撮像方法の効果を説明する図である。
図4は、本発明の第1の態様による撮像方法の他の例を示す図である。
図5は、本発明の第1の態様による撮像方法の他の例を示す図である。
図6は、本発明の第2の態様による撮像方法の手順を示す図である。
図7は、図6に示す撮像方法による体動予測関数の作成を説明する図である。
図8は、図6に示す撮像方法における本計測を説明する図である。
図9は、関数フィッティング(第2のフィッティング)を説明する図である。
図10は、本発明の第2の態様による撮像方法の一例を示す図である。
図11は、本発明の第2の態様による撮像方法の他の例を示す図である。
図12は、本発明の第2の態様が適用される撮像方法を示す図である。
図13は、図12に示す撮像方法への本発明の適用を説明する図である。
図14は、ナビゲーターエコー取得タイミングを指定するためのUIの一例を示す図である。
以下、本発明のMRI装置の実施形態を、図面を参照して説明する。
図1は、本発明が適用されるMRI装置の全体概略構成図である。このMRI装置は、被検体101が挿入される空間に静磁場を発生する磁石102と、この空間に傾斜磁場を発生する傾斜磁場コイル103と、被検体101の撮像領域に高周波磁場を発生するRFコイル104と、被検体101が発生する核磁気共鳴信号(MR)信号を検出するRFプローブ105と、静磁場空間に被検体101を挿入するためのベッド112とを備えている。
傾斜磁場コイル103は、互いに直交する3方向(X,Y,Z)の傾斜磁場コイルで構成され、傾斜磁場電源109からの信号に応じてそれぞれ傾斜磁場を発生する。これら傾斜磁場の印加の仕方によって、被検体の撮像断面を決定し、またMR信号に位置情報を付与することができる。
RFコイル104は、RF送信部110の信号に応じて高周波磁場を発生する。RFプローブ105の信号は、信号検出部106で検出され、信号処理部107で信号処理され、画像信号に変換される。また、画像は、表示部108で表示される。
傾斜磁場電源109、RF送信部110、信号検出部106は制御部111により制御される。制御動作のタイムチャートは一般にパルスシーケンスと呼ばれ、撮像方法によって決まる種々のパルスシーケンス(撮像シーケンス)が予めプログラムとして、制御部111が有する記憶部(図示せず)に格納されている。
制御部111には、このような記憶部のほか、撮像シーケンスの選択や撮像パラメータ等の入力のための入力装置(ユーザーインターフェイス:UI)が備えられている。
本発明のMRI装置では、体動モニターとしてのナビゲーターエコーを発生し、取得するためのシーケンス(ナビゲーターシーケンス)を、撮像シーケンスとは独立して或いは選択された撮像シーケンスとを組み合わせて実行することができる。
ナビゲーターシーケンスとは、モニターしたい着目部位(例えば横隔膜など)を高周波磁場及び選択傾斜磁場を用いて局所的に励起し、この局所的な励起領域から位相エンコード傾斜磁場を付加しないエコー(ナビゲーターエコー)を取得するシーケンスである。
ナビゲーターシーケンスは、例えば、撮像シーケンスの選択に合わせて、ユーザーが適宜実行することも可能であるが、撮像シーケンスと組み合わせた所定のシーケンスとして制御部111により実行される。
信号処理部107(制御部111)は、上述した画像再構成等の演算のほか、ナビゲーターシーケンスにより取得したナビゲーターエコーをもとに着目部位の変位を算出するとともに、変位方向を判断し、その結果に基づき、撮像シーケンスにおけるデータ取得のタイミングを決定する。
次に、上記構成のMRI装置を用いた呼吸動アーチファクト抑制を心電同期撮像に適用した例を説明する。
図2、図4及び図5は、それぞれ本発明の第1の態様によるMRI装置を用いた撮像方法の実施例を示す図で、図中、202は心周期、203は体動、209はシーケンス、210はデータを示し、その下にはデータの処理手順を示している。
この第1の態様においては、異なる時刻に取得した複数のナビゲーターエコーをもとに、各時刻における被検体の所定部位の変位を算出すると共に、そのときの呼吸動が呼息期であるか吸息期であるかを判定し、変位が指定された範囲であって且つ呼息期又は吸息期のいずれか一方であるときに画像用データを取得する。
まず、第1の実施例として図2の撮像方法を説明する。
心電波201を検出後、ディレイタイム206の経過後、本計測209に先行して2つのナビゲーターシーケンス207、208を続けて実行する。ナビゲーターシーケンス207、208によって取得された2つのナビゲーターエコーを用いて、続く本計測シーケンス209において、データを取得するか否かを判断する。
このため、まず、ナビゲーターエコーを周波数エンコード方向に1次元フーリエ変換して投影像を作成し(ステップ212、216)、この投影像からモニターする着目部位(例えば横隔膜など)の変位204、205を求める(ステップ213、217)。変位の検出方法は各投影像自身の信号強度プロファイルからエッジ(横隔膜の場合であれば、肺(低信号)と肝(高信号)で信号強度が急変する部分)を検出する手段などを用いる。
次に、これら変位204、205を用いて、変位の方向(二つの変位間の傾き)を求める(ステップ214)。具体的には、これら変位204、205の時間に対する1次微分を求める。或いは単純に、二つの変位204、205の差分を取ってもよく、これにより、傾きが正か負かを判断する。そして、二つの変位間の傾きが正であれば呼息期、負であれば吸息期と判断する(ステップ215)。
2つのナビゲーターエコーから呼息期又は吸息期を判定する場合、2つのナビゲーターシーケンスの時間間隔は変位の検出精度に応じて適宜設定することが望ましい。
実際にナビゲーターエコーを用いて横隔膜の呼吸動をモニターしたところ、100msで0.6〜0.7mmの割合で動いていることが分かった。従って、ナビゲーターシーケンス207と208を100ms程度の間隔で実行する場合、検出精度は0.7〜0.8mmの分解能が必要となる。ナビゲーターシーケンス207と208の間隔を広げれば、分解能もそれに応じて低くしてもよい。
次に、2つの時刻における変位204、205が、それぞれ予め設定した一定の範囲(ゲートウィンドウ211)内にあるか否かを判定する(ステップ218)。ゲートウインドウ211は、例えばナビゲーターエコーにより得られた着目部位の信号を画像化し、その位置を基準として、適当な範囲を制御部111の入力装置から入力することにより設定することができる。
ステップ215及びステップ218の判定結果より、好ましくは、呼息期で且つ両変位204、205が211内であるとき、または、吸息期で且つ両変位204、205がゲートウインドウ211内であるときの何れか一方の条件でのみ本計測データを取得する。
ただし、変位204、205がゲートウインドウ211内であることを条件に、一周期間で呼息期と吸息期を両方同時、または複数周期間で呼息期と吸息期とを混合して本計測データを取得することにより、全体の撮像時間を短縮することも可能である。
以上のようなナビゲーターシーケンス207、208と本計測シーケンス209を、心周期毎に繰り返し1枚の画像に必要なデータを取得する(ステップ219)。
呼吸変動203は周期運動をしているので、本計測シーケンス209開始直前の変位および呼吸期(呼息期または吸息期)が同じであれば、心周期毎の繰り返しにおいて本計測シーケンス209実行時も毎回同変位となり、呼吸動による位置ずれの影響、すなわち、体動アーチファクトの低減された画像を取得することができる。
この様子を図3に示す。図中、図2と同じ要素は同一符号で示した。また、301、302は異なる一心周期における体動である。一つの心周期で検出した変位203と他の心周期で検出した変位がそれぞれ指定範囲にあれば、その向きが同じなので、変位はほぼ同様に変化し(位置変化幅310は少なく)、心時相毎にみると位置ずれは一定範囲内に保たれる。
一方、変位の方向を考慮しない場合は、仮に一つの心周期では、変位203後にマイナスの方向に位置が変化し(301)、他の心周期ではプラスの方向に位置が変化した(303)とすると、心周期毎の位置変化幅311は図中下方に示したようになる。すなわち、ナビゲーターエコー取得からの経過時間が短い時相では、心周期毎の位置変化幅311が小さくてすむが、ナビゲーターエコー取得からの経過時間が長くなると、心周期毎の位置変化幅311が増大し、呼吸動アーチファクトの抑制効果を得ることができない。
上記のように変位の方向を考慮した結果、各心時相において、一様に体動アーチファクトの抑制された画像を得ることができる。このような効果に加えて、ナビゲーターシーケンスは一心拍につき2回しか実行されないので、従来技術に比べナビゲーターシーケンスによるデッドタイムが低減し、多くの心時相の画像取得が可能となる。
第2の実施例を図4に示す。図中、402、403は心周期、404は体動、415は本計測シーケンス、416はデータを示す。図2に示す実施例では各心周期の本計測開始前に2回続けてナビゲーターエコーを取得したが、図4に示す実施例では本計測前後に1回ずつ取得する。心電波401を検出後、ディレイタイム409が経過した後、本計測シーケンス415に先行してナビゲーターシーケンス410を実行する。
ナビゲーターシーケンス410によって取得されたナビゲーターエコーは周波数エンコード方向に1次元フーリエ変換され投影像が作成され(ステップ420)、着目部位の変位405が求められる(ステップ421)。続いて、本計測シーケンス415を実行して本計測データ416を取得する(ステップ422)。その後、2つ目のナビゲーターシーケンス411が実行され、取得されたナビゲーターエコーを再構成し(ステップ423)、変位406が求められる(ステップ424)。
取得されられた2つの変位405と406の間で一次微分もしくは差分をとり(ステップ425)、呼吸動波形404の傾きが正か負かを求め、呼息期にあるか吸息期にあるかの判断を行う(ステップ426)。
次に、2つのナビゲーターエコーから求めた変位405、406がゲートウインドウ419内にあるか否かを判定する(ステップ427)。上記2つの判定ステップ426、427の判定結果より、呼息期で且つ変位405、406がいずれもゲートウインドウ419内にあるとき、または、吸息期で且つ変位405、406がいずれもゲートウインドウ419内にあるときの何れか一方の条件を満たすときに取得された本計測データを画像再構成に使用する(ステップ428)。
図示する例では、呼息期と判定され且つ変位405、406がいずれもゲートウインドウ419内である心周期402で本計測データ416を画像再構成に使用し、2つの変位407と408から吸息期であると判定された心周期403では本計測データ418を画像再構成に使用しない。
しかも、この例では変位408がゲートウインドウ419の外にある場合を示しているが、ゲートウインドウ内であっても吸息期であるので本計測データ418を画像再構成には使用しない。
なお、吸息期と判定され且つ変位405、406がいずれもゲートウインドウ419内である心周期402で本計測データ416を画像再構成に使用する場合も、変位408がゲートウインドウ419内であっても呼息期であるときは、本計測データ418を画像再構成には使用しない。
また、本実施例においても、変位405、406がゲートウインドウ419内であることを条件に、一周期間で呼息期と吸息期を両方同時、または複数周期間で呼息期と吸息期を混合して本計測データを画像再構成に使用することにより、全体の撮像時間を短縮することも可能である。
また、図4では、425→426→427の順で処理が行われる例を示したが、呼息期と吸息期の判定処理と変位がゲートウィンドウ内に有るか否かの判定の順序を交換しても良い。即ち、427→425→426の順序で処理を進めても良い。
このように、第2の実施例においても、第1の実施例と同様に、呼吸動アーチファクトが無く多くの心時相の画像取得が可能となる。また、呼吸動404は周期運動をしていることから、本計測前後で測定した変位および呼吸期が分かれば本計測シーケンスの間の呼吸動変位も推測でき、毎回同変位でデータ取得可能である。
さらに、本実施例では、本計測シーケンスの前後の変位を計測しているので、ゲートウインドウ内での本計測データの取得をより確実にできる。
第3の実施例を図5に示す。502、503は心周期、504は体動、512、514はシーケンス、513、515はデータを示す。第1及び第2の実施例では各心周期に2回ナビゲーダーエコーを取得したが、図5に示す例では本計測前に1回のみ取得する。最初の心周期502で心電波501を検出後、ディレイタイム508が経過した後、本計測に先行してナビゲーターシーケンス509を実行する。
ナビゲーターシーケンス509によって取得されたナビゲーターエコーから、その時点における着目部位の変位505が求められることは第1及び第2の実施例と同様である(ステップ516、517)。ナビゲーターシーケンス509に続いて本計測シーケンス512が実行された後、本計測データ513が取得され(ステップ518)、次の心周期503に進む。次の心周期503においても同様にディレイタイム508の後、ナビゲーターシーケンス510が実行され、エコーは再構成され(ステップ519)、変位506が得られる(ステップ520)。
心周期502で得られた変位505と心周期503で得られた変位506を用いて両者の間で一次微分もしくは差分をとり、呼吸動波形504の傾きが正か負かを求め(ステップ521)、呼息期にあるか吸息期にあるかの判断を行う(ステップ522)。
次に、変位505、506がゲートウインドウ516内にあるか否かを判定し(ステップ523)、これらステップ522、523の判定結果より、呼息期で且つ変位505、506がゲートウインドウ516内にあるとき、又は、吸息期で且つ変位505、506がゲートウインドウ516内にあるときの何れか一方の条件で本計測データ513を画像再構成に使用する(ステップ524)。
以降、同様に、第2心周期503で検出した変位506と次の心周期のナビゲーターシーケンスにより得られた変位507を用いて、呼息期か吸息期かの判定、ゲートウインドウ内か否かの判定を行い、本計測データの取得を画像再構成に使用するか否かを判断する。
ここで、変位506と507のような場合は、呼吸波形504の傾きが逆で尚且つ変位507がゲートウインドウ516の外にあることから第2心周期503の本計測シーケンス514の本計測データ515は画像再構成に使用しない。図示する場合とは異なり、吸息期であり、かつ、変位がゲートウインドウ内にあるときのみ、データを取得する場合も同様に、傾きが逆でかつ、一方の変位がゲートウインドウ外にあるときは、その間の本計測シーケンスデータを画像再構成に使用しない。
なお、この実施例でも、変位505、506がゲートウインドウ516内であることを条件に、一周期間で呼息期と吸息期を両方同時、又は複数周期間で呼息期と吸息期とを混合して本計測データを画像再構成に使用することにより、全体の撮像時間を短縮することも可能である。また、図4の場合と同様に、呼息期と吸息期の判定処理と変位がゲートウインドウ内に有るか否かの判定の順序を交換しても良い。即ち、523→521→522の順序で処理を進めても良い。
この実施例においても、前述の実施例の場合と同様に、呼吸動504は周期運動をしていることから、変位及び呼吸期が分かれば、本計測シーケンス512の間の呼吸動変位も推測でき、毎回同変位でデータ取得可能である。また、実施例1、2と同様に、呼吸動アーチファクトが無く、多くの心時相の画像取得が可能となる。
以上、説明したように、本発明の第1の態様によれば、本計測シーケンスの前或いは前後にナビゲーターシーケンスを実行し、それによって得たナビゲーターエコーから変位と変位の方向を算出し、その結果によって本計測データを取得するか否か又は画像再構成に使用するか否かを判定する。これにより、単に変位のみを用いる場合に比べ本計測データ間の位置ずれを低減し、体動アーチファクトの抑制された画像を得ることができる。
また、変位の方向を検出することにより、呼吸周期の呼息期であるか吸息期であるかを判定できるので、これらの一方のみのデータであって且つ指定された範囲内の変位にある本計測データを用いて画像を再構成することができる。
さらに、一心周期内で複数の心時相の画像データを取得する場合、心時相毎にナビゲーターシーケンスを実行する必要がないので、計測できる心時相数を多くすることができる。
次に、本発明の第2の態様によるMRI装置を用いた撮像方法を説明する。この第2の態様においても、MRI装置の構成は、第1の態様によるMRI装置と同様であるが、第2の態様では信号処理系及び制御部が、問題とする体動を近似する関数(体動予測関数)を予め格納することが特徴である。
体動予測関数は、例えば本計測に先立って、少なくとも体動の一周期にわたってナビゲーターシーケンスを実行し、ナビゲーターシーケンスによって得られたナビゲーターエコーから求めた変位を用いて被検体毎に作成することができる(第1のフィッティング)。
本計測に際しては、ナビゲーターエコーを取得すると共に、ナビゲーターエコーから求めた変位を体動予測関数でフィッティングし(第2のフィッティング)、本計測データを取得するか否かを判定する。このため、信号処理系及び制御部は、前計測としてのナビゲーターシーケンスを備え、ナビゲーターエコーを用いた関数フィッティング機能を備えている。
第2の態様によるMRI装置を用いた撮像方法の一実施例を図6〜図8に示す。図6は撮像方法の手順を示すフローチャートである。この撮像方法は、被検体毎に体動予測関数を作成し、本計測データを取得する時間と間隔を求めるまでの手順(ステップ601〜605)と、本計測のための手順(ステップ606〜616)とを含み、図7はステップ601〜605を説明する図、図8はステップ606〜616を説明する図である。
まず、体動の少なくとも一周期702にわたって連続してナビゲーターシーケンス701を実行し、複数のナビゲーターエコーを取得する(ステップ601)。ナビゲーターエコーを周波数エンコード方向に1次元フーリエ変換して投影像を作成し、この投影像からモニターする着目部位(例えば横隔膜など)の変位704、705、706・・・を、ナビゲーターエコーの取得時刻毎に求める(ステップ602)。
位置の検出方法は、各投影像の信号強度プロファイルからエッジを検出するなどの手法を用いることができる。ステップ602で求めた変位の時間変位、即ち呼吸周期を多項関数で近似する(第1のフィッティング、ステップ603)。呼吸周期を近似する関数(体動予測関数)は、例えば最小二乗法などの手法を用いて算出することができる。例えば、sin、cosなどの周期関数を含む多項式である。また、ドリフト成分を反映した1次項または体動周期より長い周期を持つ周期関数を含むことが望ましい。具体的には、次式(1)或いは、次式(2)或いは、次式(3)に示すような関数fを用いることができる。
f=Σfn(t)=a*sin(bt+c)+d −−−(1)
f=a*sin(bt+c)+dt+e −−−(2)
f=a*sin(bt+c)+d*sin(et+f)−−−(3)
ここで、nは1から始まる整数、tは所望の基準時間からの経過時間、a、b、c、d、e、fは係数である。
図7に上記ステップ601〜603で検出された体動変位704、705・・・711と体動予測関数720とを示す。呼吸の一周期702が約4秒で、ナビゲーターシーケンス701を200ms毎に実行すると、一周期に20個のナビゲーターエコーを取得することができる。
このように、体動予測関数が決定したならば、体動予測関数で予測される体動変位が指定範囲703である時間及びその間隔を求める(ステップ604)。時間は、例えば、第1回目のナビゲーターシーケンス取得時を0としてカウントしてもよいし、1回目のナビゲーターシーケンスを心電(或いは脈波)のR波に同期して開始し、R波検知のときを0としてカウントしてもよい。
次に、指定範囲703内にある体動予測関数の一次微分を求め、その値が正であるか負であるかを判定する(ステップ605)。これにより、呼吸の呼息期と吸息期とを判定する。
図7の例では、期間714は微分値が負の吸息期、期間715は微分値が正の呼息期である。これら呼息期と吸息期とについて、それぞれ、期間の開始時712、713及び期間の間隔714、715が記憶される。この場合の時間のカウントは前述の通り、1回目のナビゲーターエコー取得時或いはR波検知時を基準とする。期間714又は期間715が本計測データの取得時となる。
次に、本計測を開始する(ステップ606)。本計測は、図8(b)に示すように、ナビゲーターシーケンス801、803・・・と本計測シーケンス811、813・・・とを含み、まず少なくとも2回のナビゲーターシーケンス801、803・・・を実行し、2以上のナビゲーターエコーを取得する(ステップ607)。
これらナビゲーターエコーを再構成し、取得時における目的部位の変位802、804・・・を検出する(ステップ608)。この変位802、804・・・を予め求めた体動予測関数720(図8(a))にフィッティングする(第2のフィッティング、ステップ609)。
体動予測関数720が、例えば、図9に示すような関数である場合、本計測中に第2のフィッティングを行う方法は次のようになる。
ナビゲーターシーケンスを実行して得られた目的部位の変位P1を関数720に代入する。そこで求まる時間はT1、T1’であり、一意に決定できないが、次のナビゲーターシーケンスを実行して得られた目的部位の変位P2を関数に代入して時間T2、T2’を求め、P1に対応する時間がP2に対応する時間より遅くなるという条件を与える。これにより、時間T2’の変位にフィットしていると判断する。
このように、第2のフィッティングは、2つの変位を用いることができるが、現実には呼吸振幅は、周期毎に全く同一にはならないので、ある程度幅をもたせるために、図示するように、数回ナビゲーターシーケンスを実行し、第2のフィッティングを行うことが好ましい。
第2のフィッティングによって、ナビゲーターエコー取得時の呼吸周期における位置が確定したならば、それをもとに本計測データ取得の開始時刻及び取得時間を推定する(ステップ610)。図8に示す例では、4回のナビゲーターシーケンス801、803、805、807の実行によって得られた変位802、804、806、807を体動予測関数720にフィッティングすることによって、最後のナビゲーターエコー807取得時(Pt)から本計測シーケンスを開始するまでの時間dTとデータ取得時間を推定することができる。
ナビゲーターエコー数が少ない或いは呼吸周期の不規則な乱れなどによって時間dTとデータ取得時間が推定できない場合には、ナビゲーターエコーの取得(ステップ607)から推定までを繰り返す(ステップ611)。本計測データ取得を開始するまでの時間dTとデータ取得時間が推定できたならば、その時間が経過するまで本計測シーケンスは実行するがデータは取得しない空打ち期間812とし(ステップ612)、その後、予め推定された取得時間の間、本計測データを取得する(ステップ613)。そして、予め推定された取得時間内であるかを判定し(ステップ614)、推定取得時間内であれば、ステップ613に戻る。ステップ614で推定取得時間外になっていれば、ステップ615に進み、データ取得終了とし、ステップ616にて撮像終了か否かを判断し、終了でなければ、ステップ612に戻る。
呼吸周期の変動がなく、ほぼ一定であるとすると、それ以降、呼吸周期T毎に本計測データの取得を所定時間繰り返すことにより、常に体動が指定範囲703内であるときに本計測データを取得することができる。
本計測データ取得と次の本計測データ取得との間は本計測シーケンスは実行するがデータは取得しない空打ちをすることにより、例えば自由歳差運動定常状態の計測において画質を向上することができる。
このように、本実施例によれば、ナビゲーターエコーから得た変位を、予め求めた呼吸周期の近似関数に第2のフィッティングを行うことにより本計測データの取得開始時及び取得時間を推定するようにしたので、呼息期及び吸息期のそれぞれについて体動が指定範囲内であるときの画像を取得することができる。
また、ナビゲーターシーケンスを本計測シーケンス中に実行する必要がないので撮像時間の短縮を図ることができ、心電同期撮影の場合には、より多くの心時相のデータを取得することができる。
上記実施例は、さらに呼吸周期の不規則な変動や段階的な変化などに対応した変更が可能である。そのような変更例として、本計測と本計測との間の期間(図8の実施例では空打ち期間)を利用してナビゲーターシーケンスを行う実施例を図10及び図11に示す。
これらの実施例でも、本計測に先立って、ナビゲーターシーケンスで取得したナビゲーターエコーから呼吸周期を近似する体動予測関数を求めることは上記実施例と同じである。また、本計測において、本計測データ取得開始時を決定するために複数のナビゲーターエコーを取得することも同様である。
しかし、図10に示す実施例では、本計測データ取得から次の本計測データ取得までの間の待ち時間を利用して、ナビゲーターシーケンス901、903を複数回、実行する。このナビゲーターシーケンスで得られたナビゲーターエコー902、904はその時点における体動位置の検出に用いられるとともに、前計測によって作成された体動予測関数の更新に利用される。
即ち、まずこれらナビゲーターエコー901から検出した変位902をその際存在する体動予測関数900にフィッティング(第2のフィッティング)して(ステップ915)、本計測データの取得開始時とデータ取得時間間隔を推定する(ステップ916)。必要に応じて本計測データ取得911に先立って空打ち期間912を設定する。一方、検出した変位902を用いて新たな体動予測関数を算出する(ステップ917)。ナビゲーターシーケンス901により得られる変位902は呼吸周期の一部であるので、従前の体動予測関数とこの一部の変位とを用いて新たな体動予測関数を算出する。新たな体動予測関数で従前の体動予測関数を更新し、これを次の回の本計測913のデータ取得開始時と取得時間の推定に用いる。
本計測911後に実行されたナビゲーターシーケンス903のナビゲーターエコーから検出された変位904を、変位902によって更新された体動予測関数にフィッティングし(第2のフィッティング、ステップ918)、本計測913のデータ取得開始時と取得時間を推定し、実行する(ステップ919)。
以後、同様に本計測前に取得したナビゲーターエコーでその本計測のデータ取得開始時とデータ取得時間を推定するとともに、それ以後に用いる体動予測関数を更新しながら計測を行う。
このように、本実施例では、本計測データ取得を待つ時間を利用してナビゲーターシーケンスを実行し、逐次体動予測関数を更新するので、一時的に呼吸周期が乱れた場合や呼吸周期が変化した場合にも、データ取得時間のずれを補正することができる。
図11に示す実施例でも、本計測データの取得と次の本計測データ取得との間でナビゲーターシーケンスを実行することは図10の実施例と同じである。但し、この実施例では、一定の第2のフィッティング区間をあらかじめ設定しておき、第2のフィッティングを行うナビゲーターエコーの数を逐次追加或いは更新する。
第2のフィッティングを行う一定の区間を設定することで、定期的に本計測時間推定のための第2のフィッティングを行うことが可能となり、より実際に近い本計測時間推定ができる。
また、先に実行された本計測時間の推定によって、本計測以外の時間を予測できれば、その時間に第2のフィッティング区間をあらかじめ設定しておくことも可能である。
即ち、区間1011で検出された変位に基づき第2のフィッティングを行って本計測1021のデータ取得時を推定し、次の区間1012では区間1011に含まれるナビゲーターエコーと新たに得られたナビゲーターエコーを用いて、第2のフィッティングを行って本計測1022のデータ取得時を推定する。
以下、同様に、順次、前の区間に含まれるナビゲーターエコーの一部と新たに得られたナビゲーターエコーとを用いて第2のフィッティングを行う。
この場合、第2のフィッティング区間の間隔は、呼吸周期よりも短く、しかも第2のフィッティング区間は、それぞれオーバーラップしているため、時間の経過に伴い、例えば、第2のフィッティング区間1013で本計測1023のデータ取得時間を推定した後に、次の第2のフィッティング区間1014でも同じ本計測1023のデータ取得時間が推定されることになる。この場合には、本計測1023に近い時間に実行されたフィッティング結果に基づき本計測1023おけるデータ取得を行う。
このように、本実施例によれば、第2のフィッティングに用いるナビゲーターエコーを逐次更新していくことにより、本計測データの取得時間の推定結果をより正確にすることができる。
なお、本実施例においても、図10の実施例と同様に、フィッティングに用いるナビゲーターエコーを用いて体動予測関数を更新(第1のフィッティング)していくことが好ましい。
即ち、本計測と本計測との間に計測したナビゲーターエコーを用いて第2のフィッティングを行うとともに、体動予測関数を逐次更新(第1のフィッティング)する。これにより、呼吸周期が変化した場合にも、データ取得時間のずれを補正することができる。
以上説明したように、本発明の第2の態様によるMRI装置によれば、本計測に先立って体動周期を近似する体動予測関数を作成し、この体動予測関数と本計測に付随して取得したナビゲーターエコーとから本計測データを取得する時刻及びデータ取得時間を推測するようにしたので、本計測時における実際の変位と予測した変位とのずれを小さくすることができ、体動のアーチファクトを確実に予測することができる。
特に、呼息期と吸息期とを確実に識別して、それぞれの画像を取得することができる。また、心電同期撮像のように、一心周期内の複数の心時相について画像を得る場合に、各時相の計測毎にナビゲーターシーケンスを実行する必要がないので、実質的な計測時間を増加することができ、これにより計測時間の短縮或いは心時相数の増加を図ることができる。
さらに、本計測と本計測との間の待ち時間を利用してナビゲーターシーケンスを実行することにより、フィッティング結果をより確実にすることができ、またそのナビゲーターシーケンス結果を用いて体動予測関数を更新することにより、体動周期が変化した場合でも、ずれが発生することなく常に変位が指定範囲にあるときに本計測データを取得することができる。
また、各実施例において、体動予測関数として式(2)のような一次項を含む体動予測関数を用いた場合には、一時的な被検***置の変化にも対応することができる。
図12は、被検体の位置が動いた場合の体動波形を表した図である。図示するように、時刻t1からt2で被検体が動いたために、波形が変形し、計測部位は本来計測しようとする体動範囲1201から外れてしまっていて、初期の指定範囲のまま撮像したのでは、有効なデータを得ることができない。
このような体動波形は、例えば前記(2)式の様な一次項を含む体動予測関数又は前記(3)式の様な体動周期とそれより長い周期を持つ周期関数の重み付け和によって表される体動予測関数によって近似することができ、図13に示すように、この一次項又は長周期項による波形の変化に合わせて撮像位置及び指定範囲(ゲートウインドウ)を更新する(1201→1202→1203と変更)ことにより、計測部位が指定範囲内にあるときに本計測データの取得を行うことができる。また、ゲートウインドウの更新は、体動予測関数を用いずに行なうことも可能である。例えば、1呼吸周期毎(図4の例の場合405、406、407、408・・・)に位置の平均値を計測し、平均値の変化からゲートウインドウを図13のように変化させる。あるいは、後述する図14の様なUIを持つ画面上で直接ゲートウインドウを設定してもよい。
なお、撮像位置の更新は、リアルタイムで撮像位置を変化させるインタラクティブスキャンコントロール(ISC)技術により実現できる。
本発明のMRI装置において、ナビゲーターエコーを取得するタイミング指定するためのUIの実施形態を図14に示す。図14はUIとして表示される画面の一例を示すものであり、心電計から得られたR波1401、1402、体動予測関数1403、ゲートウィンドウゲートウィンドウ1404が表示されるとともに、ナビゲーターエコーの取得タイミングを指定するためのメニュー1405が表示される。
メニュー1405では、例えば、本計測の前か、本計測の前後かを指定することができるとともに、ナビゲーターエコーの数を指定することができる。例えば、メニューで「Head1」を選択すると本計測シーケンスの前で1つのナビゲーターエコーを取得し、「Head2」を選択すると、本計測シーケンスの前で2つのナビゲーターエコーを取得することが選択される。
また、「Head&Tail」を選択すると本計測シーケンスの前後でそれぞれ1つのナビゲーターエコーを取得することが選択される。
また、ナビゲーターシーケンスの実行タイミングの指定は、例えば、体動予測関数1403上に表示された点を指定する。指定の内容は、模式的な図形として表示される。図示する例では、本計測シーケンス1407の前に2回のナビゲーターシーケンス1406の実行が指定されたことが表示されている。
UIは図示するものに限定されないが、このように、R波、体動予測関数、ナビゲーターシーケンス及び本計測シーケンスを表示すると共に、ナビゲーターシーケンスの数、実行のタイミングを画面上で指定可能にすることにより、本発明の実行を容易にすることができる。
図1は、本発明が適用されるMRI装置の全体概略構成図である。このMRI装置は、被検体101が挿入される空間に静磁場を発生する磁石102と、この空間に傾斜磁場を発生する傾斜磁場コイル103と、被検体101の撮像領域に高周波磁場を発生するRFコイル104と、被検体101が発生する核磁気共鳴信号(MR)信号を検出するRFプローブ105と、静磁場空間に被検体101を挿入するためのベッド112とを備えている。
傾斜磁場コイル103は、互いに直交する3方向(X,Y,Z)の傾斜磁場コイルで構成され、傾斜磁場電源109からの信号に応じてそれぞれ傾斜磁場を発生する。これら傾斜磁場の印加の仕方によって、被検体の撮像断面を決定し、またMR信号に位置情報を付与することができる。
RFコイル104は、RF送信部110の信号に応じて高周波磁場を発生する。RFプローブ105の信号は、信号検出部106で検出され、信号処理部107で信号処理され、画像信号に変換される。また、画像は、表示部108で表示される。
傾斜磁場電源109、RF送信部110、信号検出部106は制御部111により制御される。制御動作のタイムチャートは一般にパルスシーケンスと呼ばれ、撮像方法によって決まる種々のパルスシーケンス(撮像シーケンス)が予めプログラムとして、制御部111が有する記憶部(図示せず)に格納されている。
制御部111には、このような記憶部のほか、撮像シーケンスの選択や撮像パラメータ等の入力のための入力装置(ユーザーインターフェイス:UI)が備えられている。
本発明のMRI装置では、体動モニターとしてのナビゲーターエコーを発生し、取得するためのシーケンス(ナビゲーターシーケンス)を、撮像シーケンスとは独立して或いは選択された撮像シーケンスとを組み合わせて実行することができる。
ナビゲーターシーケンスとは、モニターしたい着目部位(例えば横隔膜など)を高周波磁場及び選択傾斜磁場を用いて局所的に励起し、この局所的な励起領域から位相エンコード傾斜磁場を付加しないエコー(ナビゲーターエコー)を取得するシーケンスである。
ナビゲーターシーケンスは、例えば、撮像シーケンスの選択に合わせて、ユーザーが適宜実行することも可能であるが、撮像シーケンスと組み合わせた所定のシーケンスとして制御部111により実行される。
信号処理部107(制御部111)は、上述した画像再構成等の演算のほか、ナビゲーターシーケンスにより取得したナビゲーターエコーをもとに着目部位の変位を算出するとともに、変位方向を判断し、その結果に基づき、撮像シーケンスにおけるデータ取得のタイミングを決定する。
次に、上記構成のMRI装置を用いた呼吸動アーチファクト抑制を心電同期撮像に適用した例を説明する。
図2、図4及び図5は、それぞれ本発明の第1の態様によるMRI装置を用いた撮像方法の実施例を示す図で、図中、202は心周期、203は体動、209はシーケンス、210はデータを示し、その下にはデータの処理手順を示している。
この第1の態様においては、異なる時刻に取得した複数のナビゲーターエコーをもとに、各時刻における被検体の所定部位の変位を算出すると共に、そのときの呼吸動が呼息期であるか吸息期であるかを判定し、変位が指定された範囲であって且つ呼息期又は吸息期のいずれか一方であるときに画像用データを取得する。
まず、第1の実施例として図2の撮像方法を説明する。
心電波201を検出後、ディレイタイム206の経過後、本計測209に先行して2つのナビゲーターシーケンス207、208を続けて実行する。ナビゲーターシーケンス207、208によって取得された2つのナビゲーターエコーを用いて、続く本計測シーケンス209において、データを取得するか否かを判断する。
このため、まず、ナビゲーターエコーを周波数エンコード方向に1次元フーリエ変換して投影像を作成し(ステップ212、216)、この投影像からモニターする着目部位(例えば横隔膜など)の変位204、205を求める(ステップ213、217)。変位の検出方法は各投影像自身の信号強度プロファイルからエッジ(横隔膜の場合であれば、肺(低信号)と肝(高信号)で信号強度が急変する部分)を検出する手段などを用いる。
次に、これら変位204、205を用いて、変位の方向(二つの変位間の傾き)を求める(ステップ214)。具体的には、これら変位204、205の時間に対する1次微分を求める。或いは単純に、二つの変位204、205の差分を取ってもよく、これにより、傾きが正か負かを判断する。そして、二つの変位間の傾きが正であれば呼息期、負であれば吸息期と判断する(ステップ215)。
2つのナビゲーターエコーから呼息期又は吸息期を判定する場合、2つのナビゲーターシーケンスの時間間隔は変位の検出精度に応じて適宜設定することが望ましい。
実際にナビゲーターエコーを用いて横隔膜の呼吸動をモニターしたところ、100msで0.6〜0.7mmの割合で動いていることが分かった。従って、ナビゲーターシーケンス207と208を100ms程度の間隔で実行する場合、検出精度は0.7〜0.8mmの分解能が必要となる。ナビゲーターシーケンス207と208の間隔を広げれば、分解能もそれに応じて低くしてもよい。
次に、2つの時刻における変位204、205が、それぞれ予め設定した一定の範囲(ゲートウィンドウ211)内にあるか否かを判定する(ステップ218)。ゲートウインドウ211は、例えばナビゲーターエコーにより得られた着目部位の信号を画像化し、その位置を基準として、適当な範囲を制御部111の入力装置から入力することにより設定することができる。
ステップ215及びステップ218の判定結果より、好ましくは、呼息期で且つ両変位204、205が211内であるとき、または、吸息期で且つ両変位204、205がゲートウインドウ211内であるときの何れか一方の条件でのみ本計測データを取得する。
ただし、変位204、205がゲートウインドウ211内であることを条件に、一周期間で呼息期と吸息期を両方同時、または複数周期間で呼息期と吸息期とを混合して本計測データを取得することにより、全体の撮像時間を短縮することも可能である。
以上のようなナビゲーターシーケンス207、208と本計測シーケンス209を、心周期毎に繰り返し1枚の画像に必要なデータを取得する(ステップ219)。
呼吸変動203は周期運動をしているので、本計測シーケンス209開始直前の変位および呼吸期(呼息期または吸息期)が同じであれば、心周期毎の繰り返しにおいて本計測シーケンス209実行時も毎回同変位となり、呼吸動による位置ずれの影響、すなわち、体動アーチファクトの低減された画像を取得することができる。
この様子を図3に示す。図中、図2と同じ要素は同一符号で示した。また、301、302は異なる一心周期における体動である。一つの心周期で検出した変位203と他の心周期で検出した変位がそれぞれ指定範囲にあれば、その向きが同じなので、変位はほぼ同様に変化し(位置変化幅310は少なく)、心時相毎にみると位置ずれは一定範囲内に保たれる。
一方、変位の方向を考慮しない場合は、仮に一つの心周期では、変位203後にマイナスの方向に位置が変化し(301)、他の心周期ではプラスの方向に位置が変化した(303)とすると、心周期毎の位置変化幅311は図中下方に示したようになる。すなわち、ナビゲーターエコー取得からの経過時間が短い時相では、心周期毎の位置変化幅311が小さくてすむが、ナビゲーターエコー取得からの経過時間が長くなると、心周期毎の位置変化幅311が増大し、呼吸動アーチファクトの抑制効果を得ることができない。
上記のように変位の方向を考慮した結果、各心時相において、一様に体動アーチファクトの抑制された画像を得ることができる。このような効果に加えて、ナビゲーターシーケンスは一心拍につき2回しか実行されないので、従来技術に比べナビゲーターシーケンスによるデッドタイムが低減し、多くの心時相の画像取得が可能となる。
第2の実施例を図4に示す。図中、402、403は心周期、404は体動、415は本計測シーケンス、416はデータを示す。図2に示す実施例では各心周期の本計測開始前に2回続けてナビゲーターエコーを取得したが、図4に示す実施例では本計測前後に1回ずつ取得する。心電波401を検出後、ディレイタイム409が経過した後、本計測シーケンス415に先行してナビゲーターシーケンス410を実行する。
ナビゲーターシーケンス410によって取得されたナビゲーターエコーは周波数エンコード方向に1次元フーリエ変換され投影像が作成され(ステップ420)、着目部位の変位405が求められる(ステップ421)。続いて、本計測シーケンス415を実行して本計測データ416を取得する(ステップ422)。その後、2つ目のナビゲーターシーケンス411が実行され、取得されたナビゲーターエコーを再構成し(ステップ423)、変位406が求められる(ステップ424)。
取得されられた2つの変位405と406の間で一次微分もしくは差分をとり(ステップ425)、呼吸動波形404の傾きが正か負かを求め、呼息期にあるか吸息期にあるかの判断を行う(ステップ426)。
次に、2つのナビゲーターエコーから求めた変位405、406がゲートウインドウ419内にあるか否かを判定する(ステップ427)。上記2つの判定ステップ426、427の判定結果より、呼息期で且つ変位405、406がいずれもゲートウインドウ419内にあるとき、または、吸息期で且つ変位405、406がいずれもゲートウインドウ419内にあるときの何れか一方の条件を満たすときに取得された本計測データを画像再構成に使用する(ステップ428)。
図示する例では、呼息期と判定され且つ変位405、406がいずれもゲートウインドウ419内である心周期402で本計測データ416を画像再構成に使用し、2つの変位407と408から吸息期であると判定された心周期403では本計測データ418を画像再構成に使用しない。
しかも、この例では変位408がゲートウインドウ419の外にある場合を示しているが、ゲートウインドウ内であっても吸息期であるので本計測データ418を画像再構成には使用しない。
なお、吸息期と判定され且つ変位405、406がいずれもゲートウインドウ419内である心周期402で本計測データ416を画像再構成に使用する場合も、変位408がゲートウインドウ419内であっても呼息期であるときは、本計測データ418を画像再構成には使用しない。
また、本実施例においても、変位405、406がゲートウインドウ419内であることを条件に、一周期間で呼息期と吸息期を両方同時、または複数周期間で呼息期と吸息期を混合して本計測データを画像再構成に使用することにより、全体の撮像時間を短縮することも可能である。
また、図4では、425→426→427の順で処理が行われる例を示したが、呼息期と吸息期の判定処理と変位がゲートウィンドウ内に有るか否かの判定の順序を交換しても良い。即ち、427→425→426の順序で処理を進めても良い。
このように、第2の実施例においても、第1の実施例と同様に、呼吸動アーチファクトが無く多くの心時相の画像取得が可能となる。また、呼吸動404は周期運動をしていることから、本計測前後で測定した変位および呼吸期が分かれば本計測シーケンスの間の呼吸動変位も推測でき、毎回同変位でデータ取得可能である。
さらに、本実施例では、本計測シーケンスの前後の変位を計測しているので、ゲートウインドウ内での本計測データの取得をより確実にできる。
第3の実施例を図5に示す。502、503は心周期、504は体動、512、514はシーケンス、513、515はデータを示す。第1及び第2の実施例では各心周期に2回ナビゲーダーエコーを取得したが、図5に示す例では本計測前に1回のみ取得する。最初の心周期502で心電波501を検出後、ディレイタイム508が経過した後、本計測に先行してナビゲーターシーケンス509を実行する。
ナビゲーターシーケンス509によって取得されたナビゲーターエコーから、その時点における着目部位の変位505が求められることは第1及び第2の実施例と同様である(ステップ516、517)。ナビゲーターシーケンス509に続いて本計測シーケンス512が実行された後、本計測データ513が取得され(ステップ518)、次の心周期503に進む。次の心周期503においても同様にディレイタイム508の後、ナビゲーターシーケンス510が実行され、エコーは再構成され(ステップ519)、変位506が得られる(ステップ520)。
心周期502で得られた変位505と心周期503で得られた変位506を用いて両者の間で一次微分もしくは差分をとり、呼吸動波形504の傾きが正か負かを求め(ステップ521)、呼息期にあるか吸息期にあるかの判断を行う(ステップ522)。
次に、変位505、506がゲートウインドウ516内にあるか否かを判定し(ステップ523)、これらステップ522、523の判定結果より、呼息期で且つ変位505、506がゲートウインドウ516内にあるとき、又は、吸息期で且つ変位505、506がゲートウインドウ516内にあるときの何れか一方の条件で本計測データ513を画像再構成に使用する(ステップ524)。
以降、同様に、第2心周期503で検出した変位506と次の心周期のナビゲーターシーケンスにより得られた変位507を用いて、呼息期か吸息期かの判定、ゲートウインドウ内か否かの判定を行い、本計測データの取得を画像再構成に使用するか否かを判断する。
ここで、変位506と507のような場合は、呼吸波形504の傾きが逆で尚且つ変位507がゲートウインドウ516の外にあることから第2心周期503の本計測シーケンス514の本計測データ515は画像再構成に使用しない。図示する場合とは異なり、吸息期であり、かつ、変位がゲートウインドウ内にあるときのみ、データを取得する場合も同様に、傾きが逆でかつ、一方の変位がゲートウインドウ外にあるときは、その間の本計測シーケンスデータを画像再構成に使用しない。
なお、この実施例でも、変位505、506がゲートウインドウ516内であることを条件に、一周期間で呼息期と吸息期を両方同時、又は複数周期間で呼息期と吸息期とを混合して本計測データを画像再構成に使用することにより、全体の撮像時間を短縮することも可能である。また、図4の場合と同様に、呼息期と吸息期の判定処理と変位がゲートウインドウ内に有るか否かの判定の順序を交換しても良い。即ち、523→521→522の順序で処理を進めても良い。
この実施例においても、前述の実施例の場合と同様に、呼吸動504は周期運動をしていることから、変位及び呼吸期が分かれば、本計測シーケンス512の間の呼吸動変位も推測でき、毎回同変位でデータ取得可能である。また、実施例1、2と同様に、呼吸動アーチファクトが無く、多くの心時相の画像取得が可能となる。
以上、説明したように、本発明の第1の態様によれば、本計測シーケンスの前或いは前後にナビゲーターシーケンスを実行し、それによって得たナビゲーターエコーから変位と変位の方向を算出し、その結果によって本計測データを取得するか否か又は画像再構成に使用するか否かを判定する。これにより、単に変位のみを用いる場合に比べ本計測データ間の位置ずれを低減し、体動アーチファクトの抑制された画像を得ることができる。
また、変位の方向を検出することにより、呼吸周期の呼息期であるか吸息期であるかを判定できるので、これらの一方のみのデータであって且つ指定された範囲内の変位にある本計測データを用いて画像を再構成することができる。
さらに、一心周期内で複数の心時相の画像データを取得する場合、心時相毎にナビゲーターシーケンスを実行する必要がないので、計測できる心時相数を多くすることができる。
次に、本発明の第2の態様によるMRI装置を用いた撮像方法を説明する。この第2の態様においても、MRI装置の構成は、第1の態様によるMRI装置と同様であるが、第2の態様では信号処理系及び制御部が、問題とする体動を近似する関数(体動予測関数)を予め格納することが特徴である。
体動予測関数は、例えば本計測に先立って、少なくとも体動の一周期にわたってナビゲーターシーケンスを実行し、ナビゲーターシーケンスによって得られたナビゲーターエコーから求めた変位を用いて被検体毎に作成することができる(第1のフィッティング)。
本計測に際しては、ナビゲーターエコーを取得すると共に、ナビゲーターエコーから求めた変位を体動予測関数でフィッティングし(第2のフィッティング)、本計測データを取得するか否かを判定する。このため、信号処理系及び制御部は、前計測としてのナビゲーターシーケンスを備え、ナビゲーターエコーを用いた関数フィッティング機能を備えている。
第2の態様によるMRI装置を用いた撮像方法の一実施例を図6〜図8に示す。図6は撮像方法の手順を示すフローチャートである。この撮像方法は、被検体毎に体動予測関数を作成し、本計測データを取得する時間と間隔を求めるまでの手順(ステップ601〜605)と、本計測のための手順(ステップ606〜616)とを含み、図7はステップ601〜605を説明する図、図8はステップ606〜616を説明する図である。
まず、体動の少なくとも一周期702にわたって連続してナビゲーターシーケンス701を実行し、複数のナビゲーターエコーを取得する(ステップ601)。ナビゲーターエコーを周波数エンコード方向に1次元フーリエ変換して投影像を作成し、この投影像からモニターする着目部位(例えば横隔膜など)の変位704、705、706・・・を、ナビゲーターエコーの取得時刻毎に求める(ステップ602)。
位置の検出方法は、各投影像の信号強度プロファイルからエッジを検出するなどの手法を用いることができる。ステップ602で求めた変位の時間変位、即ち呼吸周期を多項関数で近似する(第1のフィッティング、ステップ603)。呼吸周期を近似する関数(体動予測関数)は、例えば最小二乗法などの手法を用いて算出することができる。例えば、sin、cosなどの周期関数を含む多項式である。また、ドリフト成分を反映した1次項または体動周期より長い周期を持つ周期関数を含むことが望ましい。具体的には、次式(1)或いは、次式(2)或いは、次式(3)に示すような関数fを用いることができる。
f=Σfn(t)=a*sin(bt+c)+d −−−(1)
f=a*sin(bt+c)+dt+e −−−(2)
f=a*sin(bt+c)+d*sin(et+f)−−−(3)
ここで、nは1から始まる整数、tは所望の基準時間からの経過時間、a、b、c、d、e、fは係数である。
図7に上記ステップ601〜603で検出された体動変位704、705・・・711と体動予測関数720とを示す。呼吸の一周期702が約4秒で、ナビゲーターシーケンス701を200ms毎に実行すると、一周期に20個のナビゲーターエコーを取得することができる。
このように、体動予測関数が決定したならば、体動予測関数で予測される体動変位が指定範囲703である時間及びその間隔を求める(ステップ604)。時間は、例えば、第1回目のナビゲーターシーケンス取得時を0としてカウントしてもよいし、1回目のナビゲーターシーケンスを心電(或いは脈波)のR波に同期して開始し、R波検知のときを0としてカウントしてもよい。
次に、指定範囲703内にある体動予測関数の一次微分を求め、その値が正であるか負であるかを判定する(ステップ605)。これにより、呼吸の呼息期と吸息期とを判定する。
図7の例では、期間714は微分値が負の吸息期、期間715は微分値が正の呼息期である。これら呼息期と吸息期とについて、それぞれ、期間の開始時712、713及び期間の間隔714、715が記憶される。この場合の時間のカウントは前述の通り、1回目のナビゲーターエコー取得時或いはR波検知時を基準とする。期間714又は期間715が本計測データの取得時となる。
次に、本計測を開始する(ステップ606)。本計測は、図8(b)に示すように、ナビゲーターシーケンス801、803・・・と本計測シーケンス811、813・・・とを含み、まず少なくとも2回のナビゲーターシーケンス801、803・・・を実行し、2以上のナビゲーターエコーを取得する(ステップ607)。
これらナビゲーターエコーを再構成し、取得時における目的部位の変位802、804・・・を検出する(ステップ608)。この変位802、804・・・を予め求めた体動予測関数720(図8(a))にフィッティングする(第2のフィッティング、ステップ609)。
体動予測関数720が、例えば、図9に示すような関数である場合、本計測中に第2のフィッティングを行う方法は次のようになる。
ナビゲーターシーケンスを実行して得られた目的部位の変位P1を関数720に代入する。そこで求まる時間はT1、T1’であり、一意に決定できないが、次のナビゲーターシーケンスを実行して得られた目的部位の変位P2を関数に代入して時間T2、T2’を求め、P1に対応する時間がP2に対応する時間より遅くなるという条件を与える。これにより、時間T2’の変位にフィットしていると判断する。
このように、第2のフィッティングは、2つの変位を用いることができるが、現実には呼吸振幅は、周期毎に全く同一にはならないので、ある程度幅をもたせるために、図示するように、数回ナビゲーターシーケンスを実行し、第2のフィッティングを行うことが好ましい。
第2のフィッティングによって、ナビゲーターエコー取得時の呼吸周期における位置が確定したならば、それをもとに本計測データ取得の開始時刻及び取得時間を推定する(ステップ610)。図8に示す例では、4回のナビゲーターシーケンス801、803、805、807の実行によって得られた変位802、804、806、807を体動予測関数720にフィッティングすることによって、最後のナビゲーターエコー807取得時(Pt)から本計測シーケンスを開始するまでの時間dTとデータ取得時間を推定することができる。
ナビゲーターエコー数が少ない或いは呼吸周期の不規則な乱れなどによって時間dTとデータ取得時間が推定できない場合には、ナビゲーターエコーの取得(ステップ607)から推定までを繰り返す(ステップ611)。本計測データ取得を開始するまでの時間dTとデータ取得時間が推定できたならば、その時間が経過するまで本計測シーケンスは実行するがデータは取得しない空打ち期間812とし(ステップ612)、その後、予め推定された取得時間の間、本計測データを取得する(ステップ613)。そして、予め推定された取得時間内であるかを判定し(ステップ614)、推定取得時間内であれば、ステップ613に戻る。ステップ614で推定取得時間外になっていれば、ステップ615に進み、データ取得終了とし、ステップ616にて撮像終了か否かを判断し、終了でなければ、ステップ612に戻る。
呼吸周期の変動がなく、ほぼ一定であるとすると、それ以降、呼吸周期T毎に本計測データの取得を所定時間繰り返すことにより、常に体動が指定範囲703内であるときに本計測データを取得することができる。
本計測データ取得と次の本計測データ取得との間は本計測シーケンスは実行するがデータは取得しない空打ちをすることにより、例えば自由歳差運動定常状態の計測において画質を向上することができる。
このように、本実施例によれば、ナビゲーターエコーから得た変位を、予め求めた呼吸周期の近似関数に第2のフィッティングを行うことにより本計測データの取得開始時及び取得時間を推定するようにしたので、呼息期及び吸息期のそれぞれについて体動が指定範囲内であるときの画像を取得することができる。
また、ナビゲーターシーケンスを本計測シーケンス中に実行する必要がないので撮像時間の短縮を図ることができ、心電同期撮影の場合には、より多くの心時相のデータを取得することができる。
上記実施例は、さらに呼吸周期の不規則な変動や段階的な変化などに対応した変更が可能である。そのような変更例として、本計測と本計測との間の期間(図8の実施例では空打ち期間)を利用してナビゲーターシーケンスを行う実施例を図10及び図11に示す。
これらの実施例でも、本計測に先立って、ナビゲーターシーケンスで取得したナビゲーターエコーから呼吸周期を近似する体動予測関数を求めることは上記実施例と同じである。また、本計測において、本計測データ取得開始時を決定するために複数のナビゲーターエコーを取得することも同様である。
しかし、図10に示す実施例では、本計測データ取得から次の本計測データ取得までの間の待ち時間を利用して、ナビゲーターシーケンス901、903を複数回、実行する。このナビゲーターシーケンスで得られたナビゲーターエコー902、904はその時点における体動位置の検出に用いられるとともに、前計測によって作成された体動予測関数の更新に利用される。
即ち、まずこれらナビゲーターエコー901から検出した変位902をその際存在する体動予測関数900にフィッティング(第2のフィッティング)して(ステップ915)、本計測データの取得開始時とデータ取得時間間隔を推定する(ステップ916)。必要に応じて本計測データ取得911に先立って空打ち期間912を設定する。一方、検出した変位902を用いて新たな体動予測関数を算出する(ステップ917)。ナビゲーターシーケンス901により得られる変位902は呼吸周期の一部であるので、従前の体動予測関数とこの一部の変位とを用いて新たな体動予測関数を算出する。新たな体動予測関数で従前の体動予測関数を更新し、これを次の回の本計測913のデータ取得開始時と取得時間の推定に用いる。
本計測911後に実行されたナビゲーターシーケンス903のナビゲーターエコーから検出された変位904を、変位902によって更新された体動予測関数にフィッティングし(第2のフィッティング、ステップ918)、本計測913のデータ取得開始時と取得時間を推定し、実行する(ステップ919)。
以後、同様に本計測前に取得したナビゲーターエコーでその本計測のデータ取得開始時とデータ取得時間を推定するとともに、それ以後に用いる体動予測関数を更新しながら計測を行う。
このように、本実施例では、本計測データ取得を待つ時間を利用してナビゲーターシーケンスを実行し、逐次体動予測関数を更新するので、一時的に呼吸周期が乱れた場合や呼吸周期が変化した場合にも、データ取得時間のずれを補正することができる。
図11に示す実施例でも、本計測データの取得と次の本計測データ取得との間でナビゲーターシーケンスを実行することは図10の実施例と同じである。但し、この実施例では、一定の第2のフィッティング区間をあらかじめ設定しておき、第2のフィッティングを行うナビゲーターエコーの数を逐次追加或いは更新する。
第2のフィッティングを行う一定の区間を設定することで、定期的に本計測時間推定のための第2のフィッティングを行うことが可能となり、より実際に近い本計測時間推定ができる。
また、先に実行された本計測時間の推定によって、本計測以外の時間を予測できれば、その時間に第2のフィッティング区間をあらかじめ設定しておくことも可能である。
即ち、区間1011で検出された変位に基づき第2のフィッティングを行って本計測1021のデータ取得時を推定し、次の区間1012では区間1011に含まれるナビゲーターエコーと新たに得られたナビゲーターエコーを用いて、第2のフィッティングを行って本計測1022のデータ取得時を推定する。
以下、同様に、順次、前の区間に含まれるナビゲーターエコーの一部と新たに得られたナビゲーターエコーとを用いて第2のフィッティングを行う。
この場合、第2のフィッティング区間の間隔は、呼吸周期よりも短く、しかも第2のフィッティング区間は、それぞれオーバーラップしているため、時間の経過に伴い、例えば、第2のフィッティング区間1013で本計測1023のデータ取得時間を推定した後に、次の第2のフィッティング区間1014でも同じ本計測1023のデータ取得時間が推定されることになる。この場合には、本計測1023に近い時間に実行されたフィッティング結果に基づき本計測1023おけるデータ取得を行う。
このように、本実施例によれば、第2のフィッティングに用いるナビゲーターエコーを逐次更新していくことにより、本計測データの取得時間の推定結果をより正確にすることができる。
なお、本実施例においても、図10の実施例と同様に、フィッティングに用いるナビゲーターエコーを用いて体動予測関数を更新(第1のフィッティング)していくことが好ましい。
即ち、本計測と本計測との間に計測したナビゲーターエコーを用いて第2のフィッティングを行うとともに、体動予測関数を逐次更新(第1のフィッティング)する。これにより、呼吸周期が変化した場合にも、データ取得時間のずれを補正することができる。
以上説明したように、本発明の第2の態様によるMRI装置によれば、本計測に先立って体動周期を近似する体動予測関数を作成し、この体動予測関数と本計測に付随して取得したナビゲーターエコーとから本計測データを取得する時刻及びデータ取得時間を推測するようにしたので、本計測時における実際の変位と予測した変位とのずれを小さくすることができ、体動のアーチファクトを確実に予測することができる。
特に、呼息期と吸息期とを確実に識別して、それぞれの画像を取得することができる。また、心電同期撮像のように、一心周期内の複数の心時相について画像を得る場合に、各時相の計測毎にナビゲーターシーケンスを実行する必要がないので、実質的な計測時間を増加することができ、これにより計測時間の短縮或いは心時相数の増加を図ることができる。
さらに、本計測と本計測との間の待ち時間を利用してナビゲーターシーケンスを実行することにより、フィッティング結果をより確実にすることができ、またそのナビゲーターシーケンス結果を用いて体動予測関数を更新することにより、体動周期が変化した場合でも、ずれが発生することなく常に変位が指定範囲にあるときに本計測データを取得することができる。
また、各実施例において、体動予測関数として式(2)のような一次項を含む体動予測関数を用いた場合には、一時的な被検***置の変化にも対応することができる。
図12は、被検体の位置が動いた場合の体動波形を表した図である。図示するように、時刻t1からt2で被検体が動いたために、波形が変形し、計測部位は本来計測しようとする体動範囲1201から外れてしまっていて、初期の指定範囲のまま撮像したのでは、有効なデータを得ることができない。
このような体動波形は、例えば前記(2)式の様な一次項を含む体動予測関数又は前記(3)式の様な体動周期とそれより長い周期を持つ周期関数の重み付け和によって表される体動予測関数によって近似することができ、図13に示すように、この一次項又は長周期項による波形の変化に合わせて撮像位置及び指定範囲(ゲートウインドウ)を更新する(1201→1202→1203と変更)ことにより、計測部位が指定範囲内にあるときに本計測データの取得を行うことができる。また、ゲートウインドウの更新は、体動予測関数を用いずに行なうことも可能である。例えば、1呼吸周期毎(図4の例の場合405、406、407、408・・・)に位置の平均値を計測し、平均値の変化からゲートウインドウを図13のように変化させる。あるいは、後述する図14の様なUIを持つ画面上で直接ゲートウインドウを設定してもよい。
なお、撮像位置の更新は、リアルタイムで撮像位置を変化させるインタラクティブスキャンコントロール(ISC)技術により実現できる。
本発明のMRI装置において、ナビゲーターエコーを取得するタイミング指定するためのUIの実施形態を図14に示す。図14はUIとして表示される画面の一例を示すものであり、心電計から得られたR波1401、1402、体動予測関数1403、ゲートウィンドウゲートウィンドウ1404が表示されるとともに、ナビゲーターエコーの取得タイミングを指定するためのメニュー1405が表示される。
メニュー1405では、例えば、本計測の前か、本計測の前後かを指定することができるとともに、ナビゲーターエコーの数を指定することができる。例えば、メニューで「Head1」を選択すると本計測シーケンスの前で1つのナビゲーターエコーを取得し、「Head2」を選択すると、本計測シーケンスの前で2つのナビゲーターエコーを取得することが選択される。
また、「Head&Tail」を選択すると本計測シーケンスの前後でそれぞれ1つのナビゲーターエコーを取得することが選択される。
また、ナビゲーターシーケンスの実行タイミングの指定は、例えば、体動予測関数1403上に表示された点を指定する。指定の内容は、模式的な図形として表示される。図示する例では、本計測シーケンス1407の前に2回のナビゲーターシーケンス1406の実行が指定されたことが表示されている。
UIは図示するものに限定されないが、このように、R波、体動予測関数、ナビゲーターシーケンス及び本計測シーケンスを表示すると共に、ナビゲーターシーケンスの数、実行のタイミングを画面上で指定可能にすることにより、本発明の実行を容易にすることができる。
本発明のMRI装置によれば、時間間隔を置いて計測した2以上のナビゲーターエコー信号を用いて求めた被検体の変位及び変位方向から、画像再構成用のエコー信号を取得するタイミングを決定するので、少ないナビゲーターエコー信号で指定する体動範囲内での画像再構成用エコー信号の取得を実現できる。
これにより、確実に体動アーチファクトを抑制することができる。また、ナビゲーターシーケンスに要する時間を低減することができるので、心電同期撮像などにおいて本計測時間が圧縮されるのを防止するとともに、より多くの心時相数の画像を取得することができる。
また、本発明によれば、変位方向を求め、変位方向が正又は負であるかに応じて撮像タイミングを決定することにより、呼吸周期における呼息期及び吸息期を識別することができ、これらの画像をそれぞれ撮像することができる。
さらに、本発明によれば、体動周期を近似する体動予測関数を予め求めておくことができ、この体動予測関数を用いることにより、より確実に本計測データを取得することができる。また、体動予測関数の周期性を利用して、本計測データを取得すべき時刻、取得する間隔を推定することができる。
また、本計測データを取得すべき時刻、取得する間隔の推定に基づき、それ以外の時間をナビゲーターシーケンスの実行やデータは取得しない撮像シーケンス(空打ち)の実行に当てることができ、ナビゲーターシーケンスを実行した場合は、それで得られたナビゲーターエコーを用いて、変位算出に用いるナビゲーターエコー信号の更新や体動予測関数の更新を行うことができる。これにより、体動予測関数を用いた撮像タイミングの推定をより正確にすることができ、また体動の変化があった場合のずれを最小にすることができる。
これにより、確実に体動アーチファクトを抑制することができる。また、ナビゲーターシーケンスに要する時間を低減することができるので、心電同期撮像などにおいて本計測時間が圧縮されるのを防止するとともに、より多くの心時相数の画像を取得することができる。
また、本発明によれば、変位方向を求め、変位方向が正又は負であるかに応じて撮像タイミングを決定することにより、呼吸周期における呼息期及び吸息期を識別することができ、これらの画像をそれぞれ撮像することができる。
さらに、本発明によれば、体動周期を近似する体動予測関数を予め求めておくことができ、この体動予測関数を用いることにより、より確実に本計測データを取得することができる。また、体動予測関数の周期性を利用して、本計測データを取得すべき時刻、取得する間隔を推定することができる。
また、本計測データを取得すべき時刻、取得する間隔の推定に基づき、それ以外の時間をナビゲーターシーケンスの実行やデータは取得しない撮像シーケンス(空打ち)の実行に当てることができ、ナビゲーターシーケンスを実行した場合は、それで得られたナビゲーターエコーを用いて、変位算出に用いるナビゲーターエコー信号の更新や体動予測関数の更新を行うことができる。これにより、体動予測関数を用いた撮像タイミングの推定をより正確にすることができ、また体動の変化があった場合のずれを最小にすることができる。
Claims (18)
- 所定の撮像シーケンスに基づき、画像再構成に必要なエコー信号を取得して被検体(101)の画像を得る手段と、
前記被検体の周期的体動情報を含む核磁気共鳴信号をナビゲーターエコー(207、208、410、509)として取得するナビゲーターエコー計測手段(105、106)と、
前記被検体の周期的体動に基づく所定部位の変位の許容範囲を設定する許容範囲設定手段(108、111)と、
前記ナビゲーターエコーから前記所定部位の変位を算出して該変位が前記許容範囲内にあるか否かを判定する体動判定手段(107)とを備え、
前記画像を得る手段は、前記所定部位の変位が前記許容範囲内にあるときに取得されたエコー信号を用いて画像再構成を行う磁気共鳴イメージング装置において、
前記体動判定手段は、時間間隔を置いて計測した2以上のナビゲーターエコーから、前記所定部位の変位方向を算出し、
前記所定部位の変位及び変位方向に基づいて、前記画像再構成に必要なエコー信号を取得するタイミングを決定する撮像タイミング決定手段(111)を備えたことを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置において、前記ナビゲーターエコー計測手段は、前記被検体の一心周期内に前記2以上のナビゲーターエコーを取得することを特徴とする磁気共鳴イメージング装置。
- 請求項2記載の磁気共鳴イメージング装置において、前記ナビゲーターエコー計測手段は、前記一心周期内の前記2以上のナビゲーターの内、1以上を画像再構成用エコー信号の取得期間の前に取得し、他の1以上を該画像再構成用エコー信号の取得期間の後に取得することを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、前記ナビゲーターエコー計測手段は、一心周期毎に1以上のナビゲーターエコーを取得することを特徴とする磁気共鳴イメージング装置。
- 請求項1項記載の磁気共鳴イメージング装置において、3以上のナビゲーターエコーを用いて前記周期的体動を表す体動予測関数を作成する体動予測関数作成手段(107、111)を備え、前記体動判定手段は、前記2以上のナビゲーターエコーから求めた前記所定部位の変位を前記体動予測関数にフィッティングして、所望の時刻における前記所定部位の変位及び変位方向を算出することを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、前記体動判定手段は、前記所定部位の変位の変動に対応して前記許容範囲を調整することを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、前記許容範囲設定手段は、前記2以上のナビゲーターエコーの取得タイミングと前記許容範囲を設定する手段を備えたことを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、前記撮像タイミング決定手段は、前記体動判定手段によって求められた前記所定部位の変位が前記許容範囲内であって旦つ前記変位方向が同一方向のときに、前記画像再構成用エコー信号の取得を行うように撮像タイミングを決定することを特徴とする磁気共鳴イメージング装置。
- 請求項5記載の磁気共鳴イメージング装置において、前記撮像タイミング決定手段は、前記2以上のナビゲーターエコーの取得時刻及び前記体動予測関数から、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定することを特徴とする磁気共鳴イメージング装置。
- 請求項9記載の磁気共鳴イメージング装置において、前記撮像タイミング決定手段が、前記画像再構成用エコー信号の取得開始時刻及び取得期間を決定できるまで、前記ナビゲーターエコーの計測を繰り返すことを特徴とする磁気共鳴イメージング装置。
- 請求項5記載の磁気共鳴イメージング装置において、前記体動予測関数は、少なくとも、被検体の体動の1/4周期に渡るナビゲーターエコーから求めた変位を用いて、被検体毎に作成されることを特徴とする磁気共鳴イメージング装置。
- 請求項5記載の磁気共鳴イメージング装置において、前記体動予測関数作成手段は、体動決定のために逐次取得されるナビゲーターエコーを用いて前記体動予測関数を更新することを特徴とする磁気共鳴イメージング装置。
- 請求項12記載の磁気共鳴イメージング装置において、前記画像再構成用エコー信号の取得期間終了時から次の前記画像再構成用エコー信号の取得開始時までの間に、1以上のナビゲーターエコーを取得し、該ナビゲーターエコーを用いて前記体動予測関数を更新することを特徴とする磁気共鳴イメージング装置。
- 請求項13記載の磁気共鳴イメージング装置において、前記体動予測関数作成手段は、前記ナビゲーターエコーを取得する期間と、該体動予測関数の更新に使用するナビゲーターエコーの期間をそれぞれ前記体動予測関数上で指定する手段を備えたことを特徴とする磁気共鳴イメージング装置。
- 請求項2記載の磁気共鳴イメージング装置において、前記被検体の体動が呼吸動の場合に、前記体動判定手段は、一方の前記変位方向を呼息期とし、他方の前記変位方向を吸息期と判定することを特徴とする磁気共鳴イメージング装置。
- 請求項1乃至4の内いずれか一項記載の磁気共鳴イメージング装置において、前記体動判定手段は、前記2以上のナビゲーターエコーから算出した前記所定部位の変位の一次微分又は差分に基づいて、該所定部位の変位方向を算出することを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、前記画像再構成用エコー信号の取得期間以外はエコー信号の取得を行わずに前記撮像シーケンスを継続することを特徴とする磁気共鳴イメージング装置。
- 請求項1記載の磁気共鳴イメージング装置において、2つのナビゲーターエコーの互いの時間間隔は、約100msであり、変位の検出精度は、0.7mm〜0.8mmの分解能であることを特徴とする磁気共鳴イメージング装置。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003069543 | 2003-03-14 | ||
JP2003069543 | 2003-03-14 | ||
JP2003390704 | 2003-11-20 | ||
JP2003390704 | 2003-11-20 | ||
PCT/JP2004/002849 WO2004080301A1 (ja) | 2003-03-14 | 2004-03-05 | 磁気共鳴イメージング装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPWO2004080301A1 true JPWO2004080301A1 (ja) | 2006-06-08 |
Family
ID=32992974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005503502A Pending JPWO2004080301A1 (ja) | 2003-03-14 | 2004-03-05 | 磁気共鳴イメージング装置 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JPWO2004080301A1 (ja) |
WO (1) | WO2004080301A1 (ja) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4443079B2 (ja) | 2001-09-13 | 2010-03-31 | 株式会社日立メディコ | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用rf受信コイル |
JP4558397B2 (ja) * | 2004-07-15 | 2010-10-06 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
US8352013B2 (en) * | 2005-01-18 | 2013-01-08 | Siemens Medical Solutions Usa, Inc. | Method and system for motion compensation in magnetic resonance (MR) imaging |
JP4639136B2 (ja) * | 2005-10-19 | 2011-02-23 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴イメージング装置 |
JP4896526B2 (ja) | 2006-01-11 | 2012-03-14 | 株式会社東芝 | 磁気共鳴イメージング装置 |
CN100570393C (zh) * | 2006-02-06 | 2009-12-16 | 株式会社东芝 | 磁共振成像装置及磁共振成像方法 |
JP5105848B2 (ja) * | 2006-02-06 | 2012-12-26 | 株式会社東芝 | 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法 |
JP2008154887A (ja) * | 2006-12-26 | 2008-07-10 | Ge Medical Systems Global Technology Co Llc | Mri装置 |
JP5433169B2 (ja) * | 2007-07-11 | 2014-03-05 | 株式会社東芝 | 磁気共鳴イメージング装置 |
US8649846B2 (en) | 2007-07-11 | 2014-02-11 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
JP5361236B2 (ja) | 2008-03-31 | 2013-12-04 | 株式会社東芝 | 磁気共鳴イメージング装置および撮像条件設定方法 |
JP5498060B2 (ja) * | 2009-06-04 | 2014-05-21 | 株式会社東芝 | 磁気共鳴映像装置 |
JP5627903B2 (ja) * | 2010-02-26 | 2014-11-19 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴イメージング装置 |
JP5442490B2 (ja) * | 2010-02-26 | 2014-03-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴イメージング装置 |
JP5908713B2 (ja) * | 2011-12-26 | 2016-04-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴装置およびプログラム |
CN104395933B (zh) * | 2012-06-27 | 2017-05-17 | 皇家飞利浦有限公司 | 运动参数估计 |
JP6483691B2 (ja) * | 2013-12-02 | 2019-03-13 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 磁気共鳴イメージングシステム及び方法 |
KR101759083B1 (ko) | 2015-09-10 | 2017-07-18 | 삼성전자주식회사 | 자기 공명 영상 장치 및 자기 공명 영상 획득 방법 |
JP7115847B2 (ja) * | 2017-12-26 | 2022-08-09 | 富士フイルムヘルスケア株式会社 | 磁気共鳴イメージング装置 |
CN110604549B (zh) * | 2019-09-23 | 2022-06-24 | 上海联影医疗科技股份有限公司 | 磁共振扫描仪的图像采集方法、装置、设备和介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05154130A (ja) * | 1991-12-06 | 1993-06-22 | Hitachi Ltd | 体動アーチファクト除去方法 |
DE4445782C1 (de) * | 1994-12-21 | 1996-07-25 | Siemens Ag | Verfahren zur Phasenkorrektur von Kernresonanzsignalen |
US6292684B1 (en) * | 1999-12-14 | 2001-09-18 | Ge Medical Systems Global Technology Company, Llc | Respiratory displacement and velocity measurement using navigator MRI echo signals |
JP3878429B2 (ja) * | 2001-04-05 | 2007-02-07 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
-
2004
- 2004-03-05 JP JP2005503502A patent/JPWO2004080301A1/ja active Pending
- 2004-03-05 WO PCT/JP2004/002849 patent/WO2004080301A1/ja active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2004080301A1 (ja) | 2004-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPWO2004080301A1 (ja) | 磁気共鳴イメージング装置 | |
JP4807825B2 (ja) | 磁気共鳴イメージング装置 | |
US7432710B2 (en) | Apparatus and method for reducing image artifact | |
JP5105848B2 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法 | |
EP2628016B1 (en) | Mr data acquisition using physiological monitoring | |
JP4789244B2 (ja) | 磁気共鳴イメージング装置 | |
US9532728B2 (en) | Magnetic resonance imaging apparatus | |
JP6008839B2 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
JP2006198407A (ja) | 磁気共鳴(mr)イメージングにおける動き補償方法および動き補償システム | |
US20150157277A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
JP2014030723A (ja) | 生理学的周期と協調した磁気共鳴画像化方法 | |
JP2006320527A (ja) | 磁気共鳴イメージング装置 | |
JP4711732B2 (ja) | 磁気共鳴撮影装置 | |
JP4558397B2 (ja) | 磁気共鳴イメージング装置 | |
JP4133348B2 (ja) | 核磁気共鳴を用いた検査装置 | |
JP2012183431A (ja) | 磁気共鳴イメージング装置 | |
JP3967210B2 (ja) | 磁気共鳴イメージング装置 | |
JP6133926B2 (ja) | 可変コントラストを用いたダイナミックイメージング | |
JP2007061545A (ja) | 磁気共鳴イメージング装置 | |
JP2005040416A (ja) | 磁気共鳴イメージング装置 | |
JP2641486B2 (ja) | Nmrイメージング装置 | |
JP2004305454A (ja) | 磁気共鳴イメージング装置 | |
US11766188B2 (en) | Magnetic resonance imaging device | |
JP4745650B2 (ja) | 磁気共鳴イメージング装置 | |
JP5421600B2 (ja) | 核磁気共鳴イメージング装置および核磁気共鳴イメージング装置の作動方法 |