JP3836424B2 - Magnetic resonance imaging device - Google Patents

Magnetic resonance imaging device Download PDF

Info

Publication number
JP3836424B2
JP3836424B2 JP2002376451A JP2002376451A JP3836424B2 JP 3836424 B2 JP3836424 B2 JP 3836424B2 JP 2002376451 A JP2002376451 A JP 2002376451A JP 2002376451 A JP2002376451 A JP 2002376451A JP 3836424 B2 JP3836424 B2 JP 3836424B2
Authority
JP
Japan
Prior art keywords
phase
magnetic resonance
echo
resonance imaging
imaging apparatus
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.)
Expired - Fee Related
Application number
JP2002376451A
Other languages
Japanese (ja)
Other versions
JP2004202043A (en
Inventor
勝也 平野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2002376451A priority Critical patent/JP3836424B2/en
Publication of JP2004202043A publication Critical patent/JP2004202043A/en
Application granted granted Critical
Publication of JP3836424B2 publication Critical patent/JP3836424B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、磁気共鳴撮影装置に関し、とくに、ナビゲータエコー(navigator echo)を利用する磁気共鳴撮影装置に関する。
【0002】
【従来の技術】
従来、磁気共鳴撮影装置では、撮影遂行の過程でナビゲータエコーを獲得し、それを用いて撮影結果を補正することが行われる(例えば、特許文献1参照)。
【0003】
【特許文献1】
特開2001−276017号公報(第7−10頁、図3−7)
【0004】
【発明が解決しようとする課題】
ナビゲータエコーに基づいて撮影の対象のプロファイル(profile)を作成し、所定の判定点におけるプロファイルの信号強度に基づいて対象の呼吸の位相を判定して所定の呼吸位相の像を撮影する場合、呼吸状態の経時的変化等により、当初設定した判定点が適当でなくなり、撮影が不正確になることがある。
【0005】
そこで、本発明の課題は、呼吸状態の変化にかかわらず正しい撮影を行う磁気共鳴撮影装置を実現することである。
【0006】
【課題を解決するための手段】
(1)上記の課題を解決するためのひとつの観点での本発明は、撮影の対象からナビゲータエコーを繰り返し獲得するナビゲータエコー獲得手段と、前記ナビゲータエコーに基づいて対象のプロファイルを作成する作成手段と、対象の呼吸位相が予め定められた位相に一致するか否かを予め定められた個所における前記プロファイルの信号強度に基づいて判定する判定手段と、呼吸位相が予め定められた位相に一致しない頻度が予め定められた限度を超えたとき前記判定手段の判定個所を修正する修正手段と、前記判定に基づいて撮影の対象から予め定められた呼吸位相におけるイメージングエコーを獲得するイメージングエコー獲得手段と、前記イメージングエコーに基づいて画像を生成する生成手段と、を具備することを特徴とする磁気共鳴撮影装置である。
【0007】
(2)上記の課題を解決するための他の観点での本発明は、ナビゲータエコーおよびイメージングエコーを繰り返し獲得する獲得手段と、前記ナビゲータエコーに基づいて対象のプロファイルを作成する作成手段と、対象の呼吸位相が予め定められた位相に一致するか否かを予め定められた個所における前記プロファイルの信号強度に基づいて判定する判定手段と、呼吸位相が予め定められた位相に一致しない頻度が予め定められた限度を超えたとき前記判定手段の判定個所を修正する修正手段と、前記判定に基づいて前記イメージングエコーのうち予め定められた呼吸位相におけるイメージングエコーを採用する採用手段と、前記採用されたイメージングエコーに基づいて画像を生成する生成手段と、を具備することを特徴とする磁気共鳴撮影装置である。
【0008】
上記各観点での発明では、修正手段により、呼吸位相が予め定められた位相に一致しない頻度が予め定められた限度を超えたとき判定手段の判定個所を修正するので、呼吸状態の変化にかかわらず正しい撮影を行うことができる。
【0009】
前記ナビゲータエコーはスピンエコーであることが、ナビゲータエコーとして適切な信号を得る点で好ましい。前記イメージングエコーはグラディエントエコーであることが、イメージングエコーとして適切な信号を得る点で好ましい。前記プロファイルは体軸方向の1次元プロファイルであることが、呼吸位相の判定が容易な点で好ましい。前記予め定められた位相は呼気位相であることが、呼気位相の像を撮影する点で好ましい。前記予め定められた位相は吸気位相であることが、吸気位相の像を撮影する点で好ましい。前記頻度は呼吸位相が予め定められた位相に一致しない回数または時間であることが、呼吸状態の変化に適切に対応する点で好ましい。
【0010】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、本発明は実施の形態に限定されるものではない。図1に磁気共鳴撮影装置のブロック(block)図を示す。本装置の構成によって、本発明の装置に関する実施の形態の一例が示される。
【0011】
同図に示すように、本装置はマグネットシステム(magnet system)100を有する。マグネットシステム100は主磁場コイル(coil)部102、勾配コイル部106およびRF(radio frequency)コイル部108を有する。これら各コイル部は概ね円筒状の形状を有し、互いに同軸的に配置されている。マグネットシステム100の概ね円柱状の内部空間(ボア:bore)に、撮影の対象1がクレードル(cradle)500に搭載されて図示しない搬送手段により搬入および搬出される。
【0012】
主磁場コイル部102はマグネットシステム100の内部空間に静磁場を形成する。静磁場の方向は概ね対象1の体軸の方向に平行である。すなわちいわゆる水平磁場を形成する。主磁場コイル部102は例えば超伝導コイルを用いて構成される。なお、超伝導コイルに限らず常伝導コイル等を用いて構成してもよい。
【0013】
勾配コイル部106は、互いに垂直な3軸すなわちスライス(slice)軸、位相軸および周波数軸の方向において、それぞれ静磁場強度に勾配を持たせるための3つの勾配磁場を生じる。
【0014】
静磁場空間における互いに垂直な座標軸をx,y,zとしたとき、いずれの軸もスライス軸とすることができる。その場合、残り2軸のうちの一方を位相軸とし、他方を周波数軸とする。また、スライス軸、位相軸および周波数軸は、相互間の垂直性を保ったままx,y,z軸に関して任意の傾きを持たせることも可能である。これはオブリーク(oblique)とも呼ばれる。なお、本装置では対象1の体軸の方向をz軸方向とする。
【0015】
スライス軸方向の勾配磁場をスライス勾配磁場ともいう。位相軸方向の勾配磁場を位相エンコード(encode)勾配磁場またはフェーズエンコード(phase encode)勾配磁場ともいう。周波数軸方向の勾配磁場をリードアウト(read out)勾配磁場ともいう。リードアウト勾配磁場は周波数エンコード勾配磁場と同義である。このような勾配磁場の発生を可能にするために、勾配コイル部106は図示しない3系統の勾配コイルを有する。以下、勾配磁場を単に勾配ともいう。
【0016】
RFコイル部108は静磁場空間に対象1の体内のスピン(spin)を励起するための高周波磁場を形成する。以下、高周波磁場を形成することをRF励起信号の送信ともいう。また、RF励起信号をRFパルス(pulse)ともいう。励起されたスピンが生じる電磁波すなわち磁気共鳴信号は、RFコイル部108によって受信される。
【0017】
磁気共鳴信号は、周波数ドメイン(domain)すなわちフーリエ(Fourier)空間の信号となる。位相軸方向および周波数軸方向の勾配により、磁気共鳴信号のエンコードを2軸で行うので、磁気共鳴信号は2次元フーリエ空間における信号として得られる。フェーズエンコード勾配およびリードアウト勾配は、2次元フーリエ空間における信号のサンプリング位置を決定する。以下、2次元フーリエ空間をkスペース(k−space)ともいう。
【0018】
勾配コイル部106には勾配駆動部130が接続されている。勾配駆動部130は勾配コイル部106に駆動信号を与えて勾配磁場を発生させる。勾配駆動部130は、勾配コイル部106における3系統の勾配コイルに対応して、図示しない3系統の駆動回路を有する。
【0019】
RFコイル部108にはRF駆動部140が接続されている。RF駆動部140はRFコイル部108に駆動信号を与えてRFパルスを送信し、対象1の体内のスピンを励起する。
【0020】
RFコイル部108にはデータ収集部150が接続されている。データ収集部150は、RFコイル部108が受信した受信信号をディジタルデータ(digital data)として収集する。
【0021】
勾配駆動部130、RF駆動部140およびデータ収集部150にはシーケンス(sequence)制御部160が接続されている。シーケンス制御部160は、勾配駆動部130ないしデータ収集部150をそれぞれ制御して磁気共鳴信号の収集を遂行する。
【0022】
シーケンス制御部160は、例えばコンピュータ(computer)等を用いて構成される。シーケンス制御部160は図示しないメモリ(memory)を有する。メモリはシーケンス制御部160用のプログラムおよび各種のデータを記憶している。シーケンス制御部160の機能は、コンピュータがメモリに記憶されたプログラムを実行することにより実現される。
【0023】
データ収集部150の出力側はデータ処理部170に接続されている。データ収集部150が収集したデータがデータ処理部170に入力される。データ処理部170は、例えばコンピュータ等を用いて構成される。データ処理部170は図示しないメモリを有する。メモリはデータ処理部170用のプログラムおよび各種のデータを記憶している。
【0024】
データ処理部170はシーケンス制御部160に接続されている。データ処理部170はシーケンス制御部160の上位にあってそれを統括する。本装置の機能は、データ処理部170がメモリに記憶されたプログラムを実行することによりを実現される。
【0025】
データ処理部170は、データ収集部150が収集したデータをメモリに記憶する。メモリ内にはデータ空間が形成される。このデータ空間はkスペースに対応する。データ処理部170は、kスペースのデータを2次元逆フ−リエ変換することにより画像を再構成する。
【0026】
データ処理部170には表示部180および操作部190が接続されている。表示部180は、グラフィックディスプレー(graphic display)等で構成される。操作部190はポインティングデバイス(pointing device)を備えたキーボード(keyboard)等で構成される。
【0027】
表示部180は、データ処理部170から出力される再構成画像および各種の情報を表示する。操作部190は、使用者によって操作され、各種の指令や情報等をデータ処理部170に入力する。使用者は表示部180および操作部190を通じてインタラクティブ(interactive)に本装置を操作する。
【0028】
図2に、他の方式の磁気共鳴撮影装置のブロック図を示す。同図に示す磁気共鳴撮影装置は、本発明の実施の形態の一例である。本装置の構成によって、本発明の装置に関する実施の形態の一例が示される。
【0029】
本装置は、図1に示した装置とは方式を異にするマグネットシステム100’を有する。マグネットシステム100’以外は図1に示した装置と同様な構成になっており、同様な部分は同一の符号を付して説明を省略する。
【0030】
マグネットシステム100’は主磁場マグネット部102’、勾配コイル部106’およびRFコイル部108’を有する。これら主磁場マグネット部102’および各コイル部は、いずれも空間を挟んで互いに対向する1対のものからなる。また、いずれも概ね円盤状の形状を有し中心軸を共有して配置されている。マグネットシステム100’の内部空間(ボア)に、対象1がクレードル500に搭載されて図示しない搬送手段により搬入および搬出される。
【0031】
主磁場マグネット部102’はマグネットシステム100’の内部空間に静磁場を形成する。静磁場の方向は概ね対象1の体軸方向と直交する。すなわちいわゆる垂直磁場を形成する。主磁場マグネット部102’は例えば永久磁石等を用いて構成される。なお、永久磁石に限らず超伝導電磁石あるいは常伝導電磁石等を用いて構成してもよい。
【0032】
勾配コイル部106’は、互いに垂直な3軸すなわちスライス軸、位相軸および周波数軸の方向において、それぞれ静磁場強度に勾配を持たせるための3つの勾配磁場を生じる。
【0033】
静磁場空間における互いに垂直な座標軸をx,y,zとしたとき、いずれの軸もスライス軸とすることができる。その場合、残り2軸のうちの一方を位相軸とし、他方を周波数軸とする。また、スライス軸、位相軸および周波数軸は、相互間の垂直性を保ったままx,y,z軸に関して任意の傾きを持たせること、すなわちオブリークが可能である。本装置でも対象1の体軸の方向をz軸方向とする。3軸方向の勾配磁場の発生を可能にするために、勾配コイル部106’は図示しない3系統の勾配コイルを有する。
【0034】
RFコイル部108’は静磁場空間に対象1の体内のスピンを励起するためのRFパルスを送信する 励起されたスピンが生じる電磁波すなわち磁気共鳴信号は、RFコイル部108’によって受信される。RFコイル部108’の受信信号がデータ収集部150に入力される。
【0035】
図3に、磁気共鳴撮影に用いるパルスシーケンス(pulse sequence)の一例を示す。このパルスシーケンスは、スピンエコー(SE:Spin Echo)法のパルスシーケンスである。
【0036】
すなわち、(1)はSE法におけるRF励起用の90°パルスおよび180°パルスのシーケンスであり、(2)、(3)、(4)および(5)は、同じくそれぞれ、スライス勾配Gs、リードアウト勾配Gr、フェーズエンコード勾配GpおよびスピンエコーMRのシーケンスである。なお、90°パルスおよび180°パルスはそれぞれ中心信号で代表する。パルスシーケンスは時間軸tに沿って左から右に進行する。
【0037】
同図に示すように、90°パルスによりスピンの90°励起が行われる。このときスライス勾配Gsが印加され所定のスライスについての選択励起が行われる。90°励起から所定の時間後に、180°パルスによる180°励起すなわちスピン反転が行われる。このときもスライス勾配Gsが印加され、同じスライスについての選択的反転が行われる。
【0038】
90°励起とスピン反転の間の期間に、リードアウト勾配Grおよびフェーズエンコード勾配Gpが印加される。リードアウト勾配Grによりスピンのディフェーズ(dephase)が行われる。フェーズエンコード勾配Gpによりスピンのフェーズエンコードが行われる。
【0039】
スピン反転後、リードアウト勾配Grでスピンをリフェーズ(rephase)してスピンエコーMRを発生させる。スピンエコーMRはデータ収集部150によりビューデータ(view data)として収集される。このようなパルスシーケンスが周期TR(repetition time)で64〜512回繰り返される。繰り返しのたびにフェーズエンコード勾配Gpを変更し、毎回異なるフェーズエンコードを行う。これによって、64〜512ビューのビューデータが得られる。
【0040】
磁気共鳴撮影用パルスシーケンスの他の例を図4に示す。このパルスシーケンスは、グラディエントエコー(GRE:Gradient Echo)法のパルスシーケンスである。
【0041】
すなわち、(1)はGRE法におけるRF励起用のα°パルスのシーケンスであり、(2)、(3)、(4)および(5)は、同じくそれぞれ、スライス勾配Gs、リードアウト勾配Gr、フェーズエンコード勾配GpおよびグラディエントエコーMRのシーケンスである。なお、α°パルスは中心信号で代表する。パルスシーケンスは時間軸tに沿って左から右に進行する。
【0042】
同図に示すように、α°パルスによりスピンのα°励起が行われる。αは90以下である。このときスライス勾配Gsが印加され所定のスライスについての選択励起が行われる。
【0043】
α°励起後、フェーズエンコード勾配Gpによりスピンのフェーズエンコードが行われる。次に、リードアウト勾配Grによりまずスピンをディフェーズし、次いでスピンをリフェーズして、グラディエントエコーMRを発生させる。グラディエントエコーMRはデータ収集部150によりビューデータとして収集される。このようなパルスシーケンスが周期TRで64〜512回繰り返される。繰り返しのたびにフェーズエンコード勾配Gpを変更し、毎回異なるフェーズエンコードを行う。これによって、64〜512ビューのビューデータが得られる。
【0044】
図3または図4のパルスシーケンスによって得られたビューデータが、データ処理部170のメモリに収集される。なお、パルスシーケンスはSE法またはGRE法に限るものではなく、例えばファーストスピンエコー(FSE:FastSpin Echo)法等、他の適宜の技法のものであってよい。データ処理部170は、メモリに収集したビューデータに基づいて画像を再構成する。
【0045】
このような撮影が、対象1の呼吸の所定の位相について行われる。すなわち呼吸同期の撮影が行われる。呼吸同期撮影を行う場合の、本装置の動作のフロー(flow)図を図5に示す。
【0046】
同図に示すように、ステップ(step)501で、位相指定が行われる。位相指定のフロー図を図6に示す。同図に示すように、ステップ601で、ナビゲータエコー獲得が行われる。ナビゲータエコー獲得には図7に示すようなパルスシーケンスが用いられる。このパルスシーケンスは、図3に示したスピンエコー法のパルスシーケンスにおいて、位相エンコード勾配Gpを0としたものに相当する。
【0047】
このようなパルスシーケンスを用いて、図8に示すように、対象1のスライスAについて90°励起を行い、スライスBについて180°励起を行い、両スライスが交差する直線部Cから生じるスピンエコーをナビゲータエコーとして読み出す。スライスAおよびBは、直線部Cが体軸と概ね平行でかつ体軸からずれ、肺野から腹部にかけて延在するように設定される。
【0048】
なお、90°励起をスライスBについて行い、180°励起をスライスAについて行うようにしてもよい。スライスA,Bは、基本的には、後述のイメージングエコーを獲得するスライスとは無関係に、ナビゲータエコーを適切に獲得できる位置に設定される。
【0049】
次に、ステップ603で、プロファイル作成が行われる。プロファイル作成は、ナビゲータエコーを1次元逆フーリエ変換することによって行われる。これによって、直線C上のスピンの信号強度分布すなわち1次元のプロファイルが得られる。
【0050】
次に、ステップ605で、プロファイル表示が行われる。これによって表示部180にプロファイルが表示される。直線Cが肺野から腹部にかけて設定されているので、例えば、図9に示すようなプロファイルが得られる。このプロファイルは、中間部を境にして信号強度が大きく異なる。信号強度が小さい部分は肺野に相当し、信号強度が大きい部分は腹部に相当し、中間部が横隔膜に相当する。
【0051】
横隔膜に相当する部分は呼吸とともに体軸方向に往復運動する。以下、横隔膜に相当する部分を単に横隔膜ともいう。呼吸に伴う横隔膜の変位を図10に示す。同図に示すように、横隔膜は吸気時に腹部側に移動し、呼気時に胸部側に移動する。
【0052】
次に、ステップ607で、プロファイル上に位相判定点を設定することが行われる。位相判定点の設定は、使用者により操作部190を用いて行われる。これによって、例えば、体軸上の一点が位相判定点Pとして設定される。位相判定点Pは、横隔膜の吸気時の最大変位位置と呼気時の最大変位位置と間の所望の位置に設定される。
【0053】
図11〜13に、位相判定点Pの設定の3つの例をそれぞれ示す。図11は呼気側寄りに設定した例である。このような設定により、撮影位相として呼気位相が指定される。図12は吸気側寄りに設定した例である。このような設定により吸気位相が撮影位相として指定される。図13は呼気と吸気の中間に設定した例である。これにより呼気と吸気の中間の位相が指定される。
【0054】
呼吸位相は横隔膜が位相判定点Pを中心とする所定距離±αの範囲に入っているか否かに基づいて判定される。以下、この範囲を判定領域ともいう。図11に示したような設定では、破線で示すように、横隔膜が判定領域に入っているとき呼気位相と判定される。
【0055】
横隔膜が判定領域に入っているか否かの検出は判定領域におけるプロファイルの形状に基づいて行われる。プロファイルは横隔膜のところで信号強度が大きく変化するので検出は容易である。図12に示した設定では、横隔膜が判定領域に入ったことにより呼気位相と判定される。図13に示した設定では呼気と吸気の中間位相と判定される。
【0056】
図5に戻って、ステップ505で、ナビゲータエコーを獲得することが行われる。ナビゲータエコーの獲得は、図6のステップ601と同様にして行われる。すなわちスピンエコーとして獲得される。スピンエコーは90°励起と180°励起によって発生するので、プロファイルを求めるためのナビゲータエコーとして好適である。
【0057】
次に、ステップ507で、プロファイル作成が行われる。プロファイル作成は、図6のステップ601と同様に、ナビゲータエコーを用いて行われる。これによって、図9に示したプロファイルと同様なプロファイルが得られる。
【0058】
次に、ステップ511で、指定位相か否かが判定される。指定位相か否かの判定は、前述したように、判定領域に横隔膜が入っているか否かに基づいて行われる。
【0059】
指定位相のときは、ステップ523で、イメージングエコー獲得が行われる。イメージングエコーの獲得は、例えば、図3に示したパルスシーケンスまたは図4に示したパルスシーケンスによって行われる。なお、図4に示したグラディエントエコー法のパルスシーケンスの法のほうが、一回のRF励起でイメージングエコーすなわちグラディエントエコーが得られる点で好ましい。イメージングエコー獲得は1ビューないし数ビューの適宜のビュー数について行われる。
【0060】
次に、ステップ525で、全ビュー獲得済みか否かが判定され、そうでないときはステップ505に戻る。以下、呼吸位相が指定位相に一致しかつ全ビュー採用済みでないときは、ステップ505〜525の動作が繰り返される。このような動作によって、指定位相に一致する呼吸位相におけるイメージングエコーが逐次獲得される。
【0061】
ステップ511で指定位相でないと判定されたときは、ステップ535で、不一致頻度が限度を超えたか否かが判定され、限度を超えていないときはステップ505に戻る。以下、呼吸位相が指定位相に不一致でかつその頻度が限度を超えないときは、ステップ505〜535の動作が繰り返される。これによって、位相不一致の間は、ナビゲータエコーだけが獲得され、イメージングエコーは獲得されない。
【0062】
呼吸の周期性により、位相の不一致は周期的に発生する。そのため、不一致の頻度には許容限度が設けられている。この限度は、例えば、一定時間内の不一致回数とされる。あるいは、不一致の継続時間としてもよい。
【0063】
ステップ505〜525の動作およびステップ505〜535の動作の繰り返しによって、呼吸位相が指定位相に一致する期間のイメージングエコーが獲得される。なお、不一致回数または時間の計測は、位相が一致するたびにリセット(reset)される。
【0064】
このような動作の過程で、時間の経過とともに対象の入眠等により呼吸状態が変化し、プロファイルが例えば図14または図15に示すようになることがあり得る。図14は、横隔膜の移動範囲が吸気側に移動し、相対的に、判定領域が横隔膜の移動範囲から外れてしまった状態を示す。図15は、横隔膜の移動範囲が呼気側に移動し、相対的に、判定領域が横隔膜の移動範囲から外れてしまった状態を示す。
【0065】
このようになった場合は、呼吸が継続しているにもかかわらず横隔膜が判定領域に入らなくなるので、常に位相不一致と判定される。このため、不一致回数または不一致時間がやがて限度を超えるに至る。
【0066】
限度を超えた場合は、ステップ537で、位相判定点修正が行われる。位相判定点の修正は、データ処理部170により自動的に行われる。すなわち、例えば図16に示すように、腹部の信号強度と肺野の信号強度の中間の値を持つ閾値THに基づいて、横隔膜の吸気時の最大変位位置p1および呼気時の最大変位位置p2をそれぞれ求め、それら両位置の間に新たな位相判定点P’を定める。これによって、例えば呼気位相の判定を再開することが可能となる。吸気位相の判定および呼気と吸気の中間位相の判定用には、図17および図18に示すように新たな位相判定点P’をそれぞれ定める。
【0067】
このような位相判定点修正により、呼吸状態の変化にかかわらずその位相を正しく判定することができるので、指定位相に一致する呼吸位相におけるイメージングエコーを全ビュー分獲得することができる。
【0068】
全ビュー獲得済みとなったときは、ステップ527で、画像再構成が行われる。イメージングエコーは全て指定位相において獲得されたものであるから、再構成画像は高品質のものとなる。再構成された画像は、ステップ529で、表示および記憶が行われる。
【0069】
図19に、本装置の機能ブロック図を示す。同図に示すように、本装置は、ナビゲータエコー獲得部702、プロファイル作成部704、位相判定部706、イメージングエコー獲得部708、判定点修正部710および画像生成部712を有する。
【0070】
ナビゲータエコー獲得部702は、ナビゲータエコーを獲得してプロファイル作成部704に入力する。プロファイル作成部704は、ナビゲータエコーからプロファイルを作成して位相判定部706に入力する。位相判定部706は、プロファイルから呼吸位相を判定し、判定結果をイメージングエコー獲得部708および判定点修正部710に入力する。判定点修正部710は、位相不一致の頻度が限度を超えたとき、位相判定部706の位相判定点を修正する。
【0071】
イメージングエコー獲得部708は、位相判定結果に基づき所定位相のイメージングエコー獲得して画像生成部712に入力する。画像生成部712は、入力されたイメージングエコーに基づいて画像を生成する。
【0072】
ナビゲータエコー獲得部702は、図5のフロー図におけるステップ505の機能に相当する。ナビゲータエコー獲得部702は、本発明におけるナビゲータエコー獲得手段の実施の形態の一例である。プロファイル作成部704は、図5のフロー図におけるステップ507の機能に相当する。プロファイル作成部704は、本発明における作成手段の実施の形態の一例である。
【0073】
位相判定部706は、図5のフロー図におけるステップ511の機能に相当する。位相判定部706は、本発明における判定手段の実施の形態の一例である。イメージングエコー獲得部708は、図5のフロー図におけるステップ523の機能に相当する。イメージングエコー獲得部708は、本発明におけるイメージングエコー獲得手段の実施の形態の一例である。
【0074】
判定点修正部710は、図5のフロー図におけるステップ535,537の機能に相当する。判定点修正部710は、本発明における修正手段の実施の形態の一例である。画像生成部712は、図5のフロー図におけるステップ527の機能に相当する。画像生成部712は、本発明における生成手段の実施の形態の一例である。
【0075】
図20に、本装置の他の動作のフロー図を示す。同図において図5に示したものと同様の動作は同一の符号を付して説明を省略する。同図に示すように、ステップ505’で、ナビゲータエコーとイメージングエコー(imaging echo)を獲得することが行われる。ナビゲータエコーの獲得は、図6のステップ601と同様にして行われる。イメージングエコーの獲得は、例えば、図3に示したパルスシーケンスまたは図4に示したパルスシーケンスによって行われる。
【0076】
図21に、ナビゲータエコーとイメージングエコーを獲得するタイムチャート(time chart)を示す。同図に示すように、まずナビゲータエコーを獲得し、次いでイメージングエコーを獲得する。このようなエコー獲得を周期的に繰り返す。これによって、例えば64〜512ビューのイメージングエコーが逐次獲得される。なお、ナビゲータエコーを1つ獲得するたびにイメージングエコーを複数個獲得するようにしてもよい。
【0077】
このようにして、ナビゲータエコーと対にして獲得されたイメージングエコーは、ステップ511の位相判定に応じて、指定位相に一致するときはステップ523’でイメージングエコー採用が行われ、不一致のときはステップ533でイメージングエコー不採用が行われる。
【0078】
このような動作によって、ステップ505’で獲得したイメージングエコーのうち、指定位相に一致する呼吸位相で獲得したものが、画像再構成用の信号として採用され、指定位相に一致しない呼吸位相で獲得したものは採用されない。なお、採用されなかったビューのイメージングエコーはナビゲータエコーとともに取り直しが行われる。
【0079】
不採用すなわち位相不一致のときは、ステップ525で、不一致頻度が限度を超えたか否かが判定され、限度を超えていないときはステップ505’に戻るが、限度を超えた場合は、ステップ537で、位相判定点修正が行われる。位相判定点修正は前述と同様にして行われる。これによって呼吸状態の変化にかかわらず所定の位相での撮影の継続が可能となり、不採用のビューのイメージングエコーの取り直しと相まって、全ビューのデータを指定の位相で獲得することができる。
【0080】
図22に、本装置の機能ブロック図を示す。同図に示すように、本装置は、信号獲得部702’、プロファイル作成部704、位相判定部706、信号採用部708’、判定点修正部710’および画像生成部712を有する。同図において、図19に示したものと同様の部分は同一の符号を付して説明を省略する。
【0081】
信号獲得部702’は、ナビゲータエコーおよびイメージングエコーを獲得し、ナビゲータエコーをプロファイル作成部704に入力し、イメージングエコーを信号採用部708’に入力する。
【0082】
プロファイル作成部704は、ナビゲータエコーからプロファイルを作成して位相判定部706に入力する。位相判定部706は、プロファイルから呼吸位相を判定し、判定結果を信号採用部708’および判定点修正部710に入力する。
【0083】
信号採用部708’は、呼吸位相に応じてイメージングエコーの採用/不採用を決定し、採用したイメージングエコーを画像生成部712に入力する。判定点修正部710’は、位相不一致の頻度が限度を超えたとき、位相判定部706の位相判定点を修正する。画像生成部712は、採用されたイメージングエコーに基づいて画像を生成する。
【0084】
信号獲得部702’は、図20のフロー図におけるステップ505’の機能に相当する。信号獲得部702’は、本発明における獲得手段の実施の形態の一例である。プロファイル作成部704は、図20のフロー図におけるステップ507の機能に相当する。プロファイル作成部704は、本発明における作成手段の実施の形態の一例である。
【0085】
位相判定部706は、図20のフロー図におけるステップ511の機能に相当する。位相判定部706は、本発明における判定手段の実施の形態の一例である。信号採用部708’は、図20のフロー図におけるステップ523’,533の機能に相当する。信号採用部708’は、本発明における採用手段の実施の形態の一例である。
【0086】
判定点修正部710’は、図20のフロー図におけるステップ525,537の機能に相当する。判定点修正部710’は、本発明における修正手段の実施の形態の一例である。画像生成部712は、図20のフロー図におけるステップ527の機能に相当する。画像生成部712は、本発明における生成手段の実施の形態の一例である。
【0087】
【発明の効果】
以上詳細に説明したように、本発明によれば、呼吸状態の変化にかかわらず正しい撮影を行う磁気共鳴撮影装置を実現することができる。
【図面の簡単な説明】
【図1】本発明の実施の形態の一例の装置のブロック図である。
【図2】本発明の実施の形態の一例の装置のブロック図である。
【図3】本発明の実施の形態の一例の装置が実行するパルスシーケンスの一例を示す図である。
【図4】本発明の実施の形態の一例の装置が実行するパルスシーケンスの一例を示す図である。
【図5】本発明の実施の形態の一例の装置の動作のフロー図である。
【図6】図5のフロー図の一部の詳細図である。
【図7】ナビゲータエコー獲得用のパルスシーケンスの一例を示す図である。
【図8】ナビゲータエコー獲得のためのスライス選択の一例を示す図である。
【図9】プロファイルを示す図である。
【図10】プロファイルを示す図である。
【図11】位相判定点設定の一例を示す図である。
【図12】位相判定点設定の一例を示す図である。
【図13】位相判定点設定の一例を示す図である。
【図14】位相判定点の相対的なずれを示す図である。
【図15】位相判定点の相対的なずれを示す図である。
【図16】位相判定点の修正を示す図である。
【図17】位相判定点の修正を示す図である。
【図18】位相判定点の修正を示す図である。
【図19】本発明の実施の形態の一例の装置の機能ブロック図である。
【図20】本発明の実施の形態の一例の装置の動作のフロー図である。
【図21】ナビゲータエコーとイメージングエコーの獲得のタイムチャートを示す図である。
【図22】本発明の実施の形態の一例の装置の機能ブロック図である。
【符号の説明】
1 対象
100,100’ マグネットシステム
102 主磁場コイル部
102’ 主磁場マグネット部
106,106’ 勾配コイル部
108,108’ RFコイル部
130 勾配駆動部
140 RF駆動部
150 データ収集部
160 シーケンス制御部
170 データ処理部
180 表示部
190 操作部
500 クレードル
702 ナビゲータエコー獲得部
704 プロファイル作成部
706 位相判定部
708 イメージングエコー獲得部
710,710’ 判定点修正部
712 画像生成部
702’ 信号獲得部
708’ 信号採用部
[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a magnetic resonance imaging apparatus, and more particularly, to a magnetic resonance imaging apparatus using a navigator echo.
[0002]
[Prior art]
Conventionally, in a magnetic resonance imaging apparatus, navigator echoes are acquired in the course of performing imaging, and the imaging results are corrected using the navigator echo (see, for example, Patent Document 1).
[0003]
[Patent Document 1]
JP 2001-276017 A (page 7-10, FIG. 3-7)
[0004]
[Problems to be solved by the invention]
When a profile of a subject to be imaged is created based on the navigator echo, and the respiratory phase of the target is determined based on the signal intensity of the profile at a predetermined determination point, an image of the predetermined respiratory phase is captured. Due to changes in the state over time, the initially set determination point may not be appropriate, and imaging may be inaccurate.
[0005]
Accordingly, an object of the present invention is to realize a magnetic resonance imaging apparatus that performs correct imaging regardless of changes in the respiratory state.
[0006]
[Means for Solving the Problems]
(1) One aspect of the present invention for solving the above-described problem is that a navigator echo acquisition unit that repeatedly acquires a navigator echo from an object to be imaged, and a creation unit that generates a target profile based on the navigator echo Determining means for determining whether the target respiratory phase matches a predetermined phase based on the signal intensity of the profile at a predetermined location, and the respiratory phase does not match the predetermined phase Correction means for correcting the determination part of the determination means when the frequency exceeds a predetermined limit; and imaging echo acquisition means for acquiring an imaging echo in a predetermined respiratory phase from an imaging target based on the determination; And generating means for generating an image based on the imaging echo, It is a ringing imaging apparatus.
[0007]
(2) According to another aspect of the present invention for solving the above-described problem, an acquisition unit that repeatedly acquires navigator echoes and imaging echoes, a creation unit that creates a target profile based on the navigator echoes, and an object Determining means for determining whether or not the respiratory phase matches a predetermined phase based on the signal intensity of the profile at a predetermined location, and a frequency at which the respiratory phase does not match the predetermined phase A correction unit that corrects a determination portion of the determination unit when a predetermined limit is exceeded; an adoption unit that employs an imaging echo in a predetermined respiratory phase among the imaging echoes based on the determination; Generating means for generating an image based on the obtained imaging echo, and magnetic resonance It is a shadow system.
[0008]
In the inventions according to the above aspects, the correction means corrects the determination point of the determination means when the frequency at which the respiratory phase does not match the predetermined phase exceeds a predetermined limit. The correct shooting can be performed.
[0009]
The navigator echo is preferably a spin echo from the viewpoint of obtaining an appropriate signal as the navigator echo. The imaging echo is preferably a gradient echo from the viewpoint of obtaining an appropriate signal as the imaging echo. The profile is preferably a one-dimensional profile in the body axis direction from the viewpoint of easy determination of the respiratory phase. The predetermined phase is preferably an expiratory phase from the viewpoint of photographing an expiratory phase image. The predetermined phase is preferably an intake phase from the viewpoint of capturing an image of the intake phase. The frequency is preferably the number of times or the time when the respiratory phase does not match a predetermined phase, from the viewpoint of appropriately responding to a change in the respiratory state.
[0010]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. The present invention is not limited to the embodiment. FIG. 1 shows a block diagram of the magnetic resonance imaging apparatus. An example of an embodiment relating to the apparatus of the present invention is shown by the configuration of the apparatus.
[0011]
As shown in the figure, this apparatus has a magnet system 100. The magnet system 100 includes a main magnetic field coil unit 102, a gradient coil unit 106, and an RF (radio frequency) coil unit 108. Each of these coil portions has a substantially cylindrical shape and is arranged coaxially with each other. A subject 1 to be imaged is mounted on a cradle 500 in a generally cylindrical inner space (bore) of the magnet system 100 and is carried in and out by a conveying means (not shown).
[0012]
The main magnetic field coil unit 102 forms a static magnetic field in the internal space of the magnet system 100. The direction of the static magnetic field is generally parallel to the direction of the body axis of the object 1. That is, a so-called horizontal magnetic field is formed. The main magnetic field coil unit 102 is configured using, for example, a superconducting coil. In addition, you may comprise using not only a superconductive coil but a normal conductive coil.
[0013]
The gradient coil section 106 generates three gradient magnetic fields for giving gradients to the static magnetic field strength in the directions of three axes perpendicular to each other, that is, the slice axis, the phase axis, and the frequency axis.
[0014]
When the coordinate axes perpendicular to each other in the static magnetic field space are x, y, and z, any of the axes can be a slice axis. In that case, one of the remaining two axes is a phase axis, and the other is a frequency axis. In addition, the slice axis, the phase axis, and the frequency axis can have arbitrary inclinations with respect to the x, y, and z axes while maintaining the perpendicularity therebetween. This is also referred to as an oblique. In this device, the direction of the body axis of the target 1 is the z-axis direction.
[0015]
The gradient magnetic field in the slice axis direction is also called a slice gradient magnetic field. The gradient magnetic field in the phase axis direction is also referred to as a phase encode gradient magnetic field or a phase encode gradient magnetic field. The gradient magnetic field in the frequency axis direction is also referred to as a read out gradient magnetic field. The readout gradient magnetic field is synonymous with the frequency encoding gradient magnetic field. In order to make it possible to generate such a gradient magnetic field, the gradient coil unit 106 has three gradient coils (not shown). Hereinafter, the gradient magnetic field is also simply referred to as a gradient.
[0016]
The RF coil unit 108 forms a high-frequency magnetic field for exciting spins in the body of the target 1 in the static magnetic field space. Hereinafter, the formation of a high-frequency magnetic field is also referred to as transmission of an RF excitation signal. The RF excitation signal is also referred to as an RF pulse. An electromagnetic wave generated by the excited spin, that is, a magnetic resonance signal is received by the RF coil unit 108.
[0017]
The magnetic resonance signal is a signal in the frequency domain (Fourier) space. Since the magnetic resonance signal is encoded in two axes by the gradient in the phase axis direction and the frequency axis direction, the magnetic resonance signal is obtained as a signal in a two-dimensional Fourier space. The phase encoding gradient and readout gradient determine the sampling position of the signal in two-dimensional Fourier space. Hereinafter, the two-dimensional Fourier space is also referred to as k-space.
[0018]
A gradient driving unit 130 is connected to the gradient coil unit 106. The gradient driving unit 130 gives a driving signal to the gradient coil unit 106 to generate a gradient magnetic field. The gradient drive unit 130 has three systems of drive circuits (not shown) corresponding to the three systems of gradient coils in the gradient coil unit 106.
[0019]
An RF drive unit 140 is connected to the RF coil unit 108. The RF drive unit 140 gives a drive signal to the RF coil unit 108 and transmits an RF pulse to excite spins in the body of the subject 1.
[0020]
A data collection unit 150 is connected to the RF coil unit 108. The data collection unit 150 collects reception signals received by the RF coil unit 108 as digital data.
[0021]
A sequence control unit 160 is connected to the gradient driving unit 130, the RF driving unit 140, and the data collecting unit 150. The sequence controller 160 controls the gradient driver 130 or the data collector 150 to collect magnetic resonance signals.
[0022]
The sequence control unit 160 is configured using, for example, a computer. The sequence control unit 160 has a memory (not shown). The memory stores a program for the sequence control unit 160 and various data. The function of the sequence control unit 160 is realized by the computer executing a program stored in the memory.
[0023]
The output side of the data collection unit 150 is connected to the data processing unit 170. Data collected by the data collection unit 150 is input to the data processing unit 170. The data processing unit 170 is configured using, for example, a computer. The data processing unit 170 has a memory (not shown). The memory stores a program for the data processing unit 170 and various data.
[0024]
The data processing unit 170 is connected to the sequence control unit 160. The data processing unit 170 is above the sequence control unit 160 and controls it. The function of this apparatus is realized by the data processing unit 170 executing a program stored in the memory.
[0025]
The data processing unit 170 stores the data collected by the data collection unit 150 in a memory. A data space is formed in the memory. This data space corresponds to k-space. The data processing unit 170 reconstructs an image by performing two-dimensional inverse Fourier transform on k-space data.
[0026]
A display unit 180 and an operation unit 190 are connected to the data processing unit 170. The display unit 180 is configured by a graphic display or the like. The operation unit 190 includes a keyboard having a pointing device.
[0027]
The display unit 180 displays the reconstructed image and various information output from the data processing unit 170. The operation unit 190 is operated by the user and inputs various commands and information to the data processing unit 170. The user operates the apparatus interactively through the display unit 180 and the operation unit 190.
[0028]
FIG. 2 shows a block diagram of another type of magnetic resonance imaging apparatus. The magnetic resonance imaging apparatus shown in the figure is an example of an embodiment of the present invention. An example of an embodiment relating to the apparatus of the present invention is shown by the configuration of the apparatus.
[0029]
This apparatus has a magnet system 100 ′ having a different system from the apparatus shown in FIG. Except for the magnet system 100 ′, the configuration is the same as that of the apparatus shown in FIG.
[0030]
The magnet system 100 ′ includes a main magnetic field magnet unit 102 ′, a gradient coil unit 106 ′, and an RF coil unit 108 ′. Each of the main magnetic field magnet section 102 ′ and each coil section is composed of a pair of facing each other across a space. Moreover, all have a substantially disk shape and are arranged sharing the central axis. The object 1 is mounted on the cradle 500 and carried into and out of the internal space (bore) of the magnet system 100 ′ by a conveying means (not shown).
[0031]
The main magnetic field magnet unit 102 ′ forms a static magnetic field in the internal space of the magnet system 100 ′. The direction of the static magnetic field is substantially orthogonal to the body axis direction of the target 1. That is, a so-called vertical magnetic field is formed. The main magnetic field magnet section 102 ′ is configured using, for example, a permanent magnet. In addition, you may comprise using not only a permanent magnet but a superconductive electromagnet or a normal electromagnet.
[0032]
The gradient coil section 106 ′ generates three gradient magnetic fields for giving a gradient to the static magnetic field strength in the directions of three axes perpendicular to each other, that is, the slice axis, the phase axis, and the frequency axis.
[0033]
When the coordinate axes perpendicular to each other in the static magnetic field space are x, y, and z, any of the axes can be a slice axis. In that case, one of the remaining two axes is a phase axis, and the other is a frequency axis. In addition, the slice axis, the phase axis, and the frequency axis can have an arbitrary inclination with respect to the x, y, and z axes while maintaining the perpendicularity therebetween, that is, oblique. Also in this apparatus, the direction of the body axis of the target 1 is the z-axis direction. In order to enable generation of gradient magnetic fields in the three-axis directions, the gradient coil unit 106 'has three gradient coils (not shown).
[0034]
The RF coil unit 108 ′ transmits an RF pulse for exciting the spin in the body of the subject 1 to the static magnetic field space. The electromagnetic wave generated by the excited spin, that is, the magnetic resonance signal is received by the RF coil unit 108 ′. A reception signal of the RF coil unit 108 ′ is input to the data collection unit 150.
[0035]
FIG. 3 shows an example of a pulse sequence used for magnetic resonance imaging. This pulse sequence is a pulse sequence of a spin echo (SE: Spin Echo) method.
[0036]
That is, (1) is a sequence of 90 ° pulses and 180 ° pulses for RF excitation in the SE method, and (2), (3), (4) and (5) are respectively the slice gradient Gs and the lead. This is a sequence of an out gradient Gr, a phase encode gradient Gp, and a spin echo MR. The 90 ° pulse and the 180 ° pulse are represented by center signals. The pulse sequence proceeds from left to right along the time axis t.
[0037]
As shown in the figure, 90 ° excitation of spin is performed by a 90 ° pulse. At this time, the slice gradient Gs is applied, and selective excitation for a predetermined slice is performed. After a predetermined time from the 90 ° excitation, 180 ° excitation by a 180 ° pulse, that is, spin inversion is performed. At this time, the slice gradient Gs is applied, and selective inversion is performed for the same slice.
[0038]
In the period between 90 ° excitation and spin reversal, a readout gradient Gr and a phase encode gradient Gp are applied. Spin dephase is performed by the lead-out gradient Gr. Spin phase encoding is performed by the phase encoding gradient Gp.
[0039]
After the spin inversion, the spin is rephased at the readout gradient Gr to generate the spin echo MR. The spin echo MR is collected as view data by the data collecting unit 150. Such a pulse sequence is repeated 64 to 512 times with a period TR (repetition time). The phase encoding gradient Gp is changed every time it is repeated, and a different phase encoding is performed each time. Thereby, view data of 64 to 512 views is obtained.
[0040]
Another example of a pulse sequence for magnetic resonance imaging is shown in FIG. This pulse sequence is a pulse sequence of a gradient echo (GRE) method.
[0041]
That is, (1) is a sequence of α ° pulses for RF excitation in the GRE method, and (2), (3), (4) and (5) are respectively slice gradient Gs, readout gradient Gr, It is a sequence of a phase encoding gradient Gp and a gradient echo MR. The α ° pulse is represented by a center signal. The pulse sequence proceeds from left to right along the time axis t.
[0042]
As shown in the figure, the α ° excitation of the spin is performed by the α ° pulse. α is 90 or less. At this time, the slice gradient Gs is applied, and selective excitation for a predetermined slice is performed.
[0043]
After the α ° excitation, spin phase encoding is performed by the phase encoding gradient Gp. Next, the spin is first dephased by the readout gradient Gr, and then the spin is rephased to generate a gradient echo MR. The gradient echo MR is collected as view data by the data collection unit 150. Such a pulse sequence is repeated 64 to 512 times with a period TR. The phase encoding gradient Gp is changed every time it is repeated, and a different phase encoding is performed each time. Thereby, view data of 64 to 512 views is obtained.
[0044]
View data obtained by the pulse sequence of FIG. 3 or 4 is collected in the memory of the data processing unit 170. Note that the pulse sequence is not limited to the SE method or the GRE method, and may be another appropriate technique such as a Fast Spin Echo (FSE) method. The data processing unit 170 reconstructs an image based on the view data collected in the memory.
[0045]
Such imaging is performed for a predetermined phase of respiration of the subject 1. That is, breathing synchronous imaging is performed. FIG. 5 shows a flow diagram of the operation of this apparatus when performing breathing synchronous imaging.
[0046]
As shown in the figure, phase designation is performed in step 501. A flow diagram for phase designation is shown in FIG. As shown in the figure, in step 601, navigator echo acquisition is performed. A pulse sequence as shown in FIG. 7 is used for navigator echo acquisition. This pulse sequence corresponds to the spin echo method pulse sequence shown in FIG. 3 with the phase encode gradient Gp set to zero.
[0047]
Using such a pulse sequence, as shown in FIG. 8, 90 ° excitation is performed on slice A of object 1, 180 ° excitation is performed on slice B, and a spin echo generated from a straight line portion C where both slices intersect is detected. Read as navigator echo. The slices A and B are set so that the straight line portion C is substantially parallel to the body axis and is offset from the body axis and extends from the lung field to the abdomen.
[0048]
Note that 90 ° excitation may be performed for slice B and 180 ° excitation may be performed for slice A. The slices A and B are basically set at positions where navigator echoes can be appropriately acquired irrespective of slices for acquiring imaging echoes described later.
[0049]
Next, in step 603, profile creation is performed. Profile creation is performed by one-dimensional inverse Fourier transform of navigator echoes. As a result, a signal intensity distribution of spins on the straight line C, that is, a one-dimensional profile is obtained.
[0050]
Next, in step 605, profile display is performed. As a result, the profile is displayed on the display unit 180. Since the straight line C is set from the lung field to the abdomen, for example, a profile as shown in FIG. 9 is obtained. This profile is greatly different in signal intensity at the middle part. The portion with low signal intensity corresponds to the lung field, the portion with high signal intensity corresponds to the abdomen, and the intermediate portion corresponds to the diaphragm.
[0051]
A portion corresponding to the diaphragm reciprocates in the body axis direction with breathing. Hereinafter, a portion corresponding to the diaphragm is also simply referred to as a diaphragm. FIG. 10 shows the displacement of the diaphragm accompanying respiration. As shown in the figure, the diaphragm moves to the abdomen during inhalation and moves to the chest during exhalation.
[0052]
Next, in step 607, a phase determination point is set on the profile. The phase determination point is set by the user using the operation unit 190. Thereby, for example, one point on the body axis is set as the phase determination point P. The phase determination point P is set to a desired position between the maximum displacement position during inhalation of the diaphragm and the maximum displacement position during expiration.
[0053]
11 to 13 show three examples of setting the phase determination point P, respectively. FIG. 11 shows an example set closer to the expiration side. With this setting, the expiration phase is designated as the imaging phase. FIG. 12 shows an example set closer to the intake side. With this setting, the intake phase is designated as the imaging phase. FIG. 13 shows an example of setting between the expiration and the inspiration. This designates an intermediate phase between expiration and inspiration.
[0054]
The respiratory phase is determined based on whether or not the diaphragm is within a predetermined distance ± α centered on the phase determination point P. Hereinafter, this range is also referred to as a determination region. In the setting as shown in FIG. 11, as shown by the broken line, when the diaphragm is in the determination region, the expiration phase is determined.
[0055]
Whether the diaphragm is in the determination region is detected based on the shape of the profile in the determination region. The profile is easy to detect because the signal intensity varies greatly at the diaphragm. In the setting shown in FIG. 12, the expiration phase is determined when the diaphragm enters the determination region. In the setting shown in FIG. 13, it is determined as an intermediate phase between expiration and inspiration.
[0056]
Returning to FIG. 5, in step 505, navigator echo is acquired. Acquisition of the navigator echo is performed in the same manner as in step 601 in FIG. That is, it is acquired as a spin echo. Since the spin echo is generated by 90 ° excitation and 180 ° excitation, it is suitable as a navigator echo for obtaining a profile.
[0057]
Next, in step 507, profile creation is performed. The profile creation is performed using navigator echoes as in step 601 in FIG. As a result, a profile similar to the profile shown in FIG. 9 is obtained.
[0058]
Next, in step 511, it is determined whether or not the designated phase is reached. As described above, whether or not the phase is the designated phase is determined based on whether or not the diaphragm is included in the determination region.
[0059]
When the phase is the designated phase, imaging echo acquisition is performed in step 523. The acquisition of the imaging echo is performed by, for example, the pulse sequence shown in FIG. 3 or the pulse sequence shown in FIG. Note that the pulse sequence method of the gradient echo method shown in FIG. 4 is preferable in that an imaging echo, that is, a gradient echo can be obtained by one RF excitation. Imaging echo acquisition is performed for an appropriate number of views from one view to several views.
[0060]
Next, in step 525, it is determined whether or not all views have been acquired. If not, the process returns to step 505. Hereinafter, when the breathing phase matches the designated phase and not all views have been adopted, the operations in steps 505 to 525 are repeated. By such an operation, imaging echoes in the respiratory phase that matches the specified phase are sequentially acquired.
[0061]
If it is determined in step 511 that the phase is not the designated phase, it is determined in step 535 whether or not the mismatch frequency exceeds the limit. If it does not exceed the limit, the process returns to step 505. Thereafter, when the respiratory phase does not match the designated phase and the frequency does not exceed the limit, the operations in steps 505 to 535 are repeated. Thereby, during the phase mismatch, only navigator echoes are acquired and imaging echoes are not acquired.
[0062]
Due to the periodicity of breathing, phase mismatch occurs periodically. Therefore, an allowable limit is set for the frequency of mismatch. This limit is, for example, the number of mismatches within a certain time. Or it is good also as continuation time of mismatch.
[0063]
By repeating the operations of Steps 505 to 525 and the operations of Steps 505 to 535, an imaging echo in a period in which the respiratory phase matches the designated phase is acquired. Note that the number of mismatches or time measurement is reset every time the phases match.
[0064]
In the course of such an operation, the respiratory state may change due to sleep on the subject as time passes, and the profile may become as shown in FIG. 14 or FIG. 15, for example. FIG. 14 shows a state where the moving range of the diaphragm has moved to the intake side, and the determination region has relatively deviated from the moving range of the diaphragm. FIG. 15 shows a state in which the movement range of the diaphragm has moved to the expiration side, and the determination region has relatively deviated from the movement range of the diaphragm.
[0065]
In such a case, the diaphragm does not enter the determination region even though breathing continues, and therefore it is always determined that the phases do not match. For this reason, the number of mismatches or the mismatch time eventually exceeds the limit.
[0066]
If the limit is exceeded, phase decision point correction is performed at step 537. The correction of the phase determination point is automatically performed by the data processing unit 170. That is, for example, as shown in FIG. 16, the maximum displacement position p1 at the time of inhalation of the diaphragm and the maximum displacement position p2 at the time of expiration are based on a threshold value TH having an intermediate value between the signal intensity of the abdomen and the signal intensity of the lung field. Each is obtained, and a new phase determination point P ′ is determined between these two positions. Thereby, for example, the determination of the expiration phase can be resumed. For the determination of the inspiration phase and the determination of the intermediate phase between expiration and inspiration, a new phase determination point P ′ is determined as shown in FIGS.
[0067]
By correcting the phase determination point in this way, it is possible to correctly determine the phase regardless of the change in the respiratory state. Therefore, it is possible to acquire imaging echoes for all views in the respiratory phase that matches the specified phase.
[0068]
When all views have been acquired, image reconstruction is performed in step 527. Since all the imaging echoes are acquired at the designated phase, the reconstructed image is of high quality. The reconstructed image is displayed and stored at step 529.
[0069]
FIG. 19 shows a functional block diagram of this apparatus. As shown in the figure, the apparatus includes a navigator echo acquisition unit 702, a profile creation unit 704, a phase determination unit 706, an imaging echo acquisition unit 708, a determination point correction unit 710, and an image generation unit 712.
[0070]
The navigator echo acquisition unit 702 acquires navigator echoes and inputs them to the profile creation unit 704. The profile creation unit 704 creates a profile from the navigator echo and inputs it to the phase determination unit 706. The phase determination unit 706 determines the respiratory phase from the profile, and inputs the determination result to the imaging echo acquisition unit 708 and the determination point correction unit 710. The determination point correction unit 710 corrects the phase determination point of the phase determination unit 706 when the frequency of phase mismatch exceeds the limit.
[0071]
The imaging echo acquisition unit 708 acquires an imaging echo having a predetermined phase based on the phase determination result and inputs the acquired imaging echo to the image generation unit 712. The image generation unit 712 generates an image based on the input imaging echo.
[0072]
The navigator echo acquisition unit 702 corresponds to the function of step 505 in the flowchart of FIG. The navigator echo acquisition unit 702 is an example of an embodiment of navigator echo acquisition means in the present invention. The profile creation unit 704 corresponds to the function of step 507 in the flowchart of FIG. The profile creation unit 704 is an example of an embodiment of creation means in the present invention.
[0073]
The phase determination unit 706 corresponds to the function of step 511 in the flowchart of FIG. The phase determination unit 706 is an example of an embodiment of determination means in the present invention. The imaging echo acquisition unit 708 corresponds to the function of step 523 in the flowchart of FIG. The imaging echo acquisition unit 708 is an example of an embodiment of the imaging echo acquisition means in the present invention.
[0074]
The decision point correction unit 710 corresponds to the functions of steps 535 and 537 in the flowchart of FIG. The determination point correction unit 710 is an example of an embodiment of the correction means in the present invention. The image generation unit 712 corresponds to the function of step 527 in the flowchart of FIG. The image generation unit 712 is an example of an embodiment of generation means in the present invention.
[0075]
FIG. 20 shows a flowchart of another operation of the present apparatus. In this figure, the same operations as those shown in FIG. As shown in the figure, navigator echoes and imaging echoes are acquired in step 505 ′. Acquisition of the navigator echo is performed in the same manner as in step 601 in FIG. The acquisition of the imaging echo is performed by, for example, the pulse sequence shown in FIG. 3 or the pulse sequence shown in FIG.
[0076]
FIG. 21 shows a time chart for acquiring navigator echoes and imaging echoes. As shown in the figure, a navigator echo is first acquired, and then an imaging echo is acquired. Such echo acquisition is repeated periodically. Thereby, for example, imaging echoes of 64 to 512 views are sequentially acquired. A plurality of imaging echoes may be acquired each time one navigator echo is acquired.
[0077]
In this way, when the imaging echo acquired as a pair with the navigator echo matches the designated phase according to the phase determination at step 511, the imaging echo is adopted at step 523 ′, and when it does not match, the step is adopted. At 533, the imaging echo is not adopted.
[0078]
By such an operation, among the imaging echoes acquired in step 505 ′, the imaging echo acquired at the respiratory phase that matches the specified phase is adopted as a signal for image reconstruction and acquired at the respiratory phase that does not match the specified phase. Things are not adopted. Note that the imaging echoes of the views that have not been adopted are reacquired together with the navigator echoes.
[0079]
If not adopted, i.e., phase mismatch, it is determined in step 525 whether or not the mismatch frequency has exceeded the limit, and if not exceeded, the process returns to step 505 ′. The phase determination point is corrected. The phase determination point correction is performed in the same manner as described above. As a result, it is possible to continue imaging at a predetermined phase regardless of the change in the respiratory state, and it is possible to acquire data of all views at a specified phase in combination with re-acquisition of imaging echoes of views that have not been adopted.
[0080]
FIG. 22 shows a functional block diagram of this apparatus. As shown in the figure, the apparatus includes a signal acquisition unit 702 ′, a profile creation unit 704, a phase determination unit 706, a signal adoption unit 708 ′, a determination point correction unit 710 ′, and an image generation unit 712. In the figure, the same parts as those shown in FIG.
[0081]
The signal acquisition unit 702 ′ acquires navigator echoes and imaging echoes, inputs the navigator echoes to the profile creation unit 704, and inputs the imaging echoes to the signal adoption unit 708 ′.
[0082]
The profile creation unit 704 creates a profile from the navigator echo and inputs it to the phase determination unit 706. The phase determination unit 706 determines the respiratory phase from the profile, and inputs the determination result to the signal adoption unit 708 ′ and the determination point correction unit 710.
[0083]
The signal adopting unit 708 ′ determines whether or not to adopt the imaging echo according to the respiratory phase, and inputs the adopted imaging echo to the image generating unit 712. The determination point correction unit 710 ′ corrects the phase determination point of the phase determination unit 706 when the frequency of phase mismatch exceeds the limit. The image generation unit 712 generates an image based on the adopted imaging echo.
[0084]
The signal acquisition unit 702 ′ corresponds to the function of step 505 ′ in the flowchart of FIG. The signal acquisition unit 702 ′ is an example of an embodiment of acquisition means in the present invention. The profile creation unit 704 corresponds to the function of step 507 in the flowchart of FIG. The profile creation unit 704 is an example of an embodiment of creation means in the present invention.
[0085]
The phase determination unit 706 corresponds to the function of step 511 in the flowchart of FIG. The phase determination unit 706 is an example of an embodiment of determination means in the present invention. The signal adoption unit 708 ′ corresponds to the functions of steps 523 ′ and 533 in the flowchart of FIG. The signal adopting unit 708 ′ is an example of an embodiment of the adopting means in the present invention.
[0086]
The determination point correction unit 710 ′ corresponds to the function of steps 525 and 537 in the flowchart of FIG. The determination point correction unit 710 ′ is an example of an embodiment of correction means in the present invention. The image generation unit 712 corresponds to the function of step 527 in the flowchart of FIG. The image generation unit 712 is an example of an embodiment of generation means in the present invention.
[0087]
【The invention's effect】
As described above in detail, according to the present invention, it is possible to realize a magnetic resonance imaging apparatus that performs correct imaging regardless of changes in the respiratory state.
[Brief description of the drawings]
FIG. 1 is a block diagram of an exemplary apparatus according to an embodiment of the present invention.
FIG. 2 is a block diagram of an exemplary apparatus according to an embodiment of the present invention.
FIG. 3 is a diagram showing an example of a pulse sequence executed by an example apparatus according to an embodiment of the present invention.
FIG. 4 is a diagram showing an example of a pulse sequence executed by an example apparatus according to an embodiment of the present invention.
FIG. 5 is a flowchart of the operation of the apparatus according to the embodiment of the present invention.
6 is a detailed view of a portion of the flow diagram of FIG.
FIG. 7 is a diagram showing an example of a pulse sequence for navigator echo acquisition.
FIG. 8 is a diagram illustrating an example of slice selection for navigator echo acquisition;
FIG. 9 is a diagram showing a profile.
FIG. 10 is a diagram showing a profile.
FIG. 11 is a diagram illustrating an example of phase determination point setting.
FIG. 12 is a diagram illustrating an example of phase determination point setting.
FIG. 13 is a diagram illustrating an example of phase determination point setting.
FIG. 14 is a diagram showing a relative shift of phase determination points.
FIG. 15 is a diagram showing a relative shift of phase determination points.
FIG. 16 is a diagram illustrating correction of a phase determination point.
FIG. 17 is a diagram illustrating correction of a phase determination point.
FIG. 18 is a diagram illustrating correction of a phase determination point.
FIG. 19 is a functional block diagram of an exemplary apparatus according to an embodiment of the present invention.
FIG. 20 is a flowchart of the operation of the apparatus according to the example of the embodiment of the present invention.
FIG. 21 is a diagram illustrating a time chart of acquisition of navigator echoes and imaging echoes.
FIG. 22 is a functional block diagram of an exemplary apparatus according to an embodiment of the present invention.
[Explanation of symbols]
1 Target 100, 100 'Magnet system 102 Main magnetic field coil unit 102' Main magnetic field magnet unit 106, 106 'Gradient coil unit 108, 108' RF coil unit 130 Gradient drive unit 140 RF drive unit 150 Data collection unit 160 Sequence control unit 170 Data processing unit 180 Display unit 190 Operation unit 500 Cradle 702 Navigator echo acquisition unit 704 Profile creation unit 706 Phase determination unit 708 Imaging echo acquisition unit 710, 710 ′ Determination point correction unit 712 Image generation unit 702 ′ Signal acquisition unit 708 ′ Signal adoption Part

Claims (9)

撮影の対象からナビゲータエコーを繰り返し獲得するナビゲータエコー獲得手段と、
前記ナビゲータエコーに基づいて対象のプロファイルを作成する作成手段と、
対象の呼吸位相が予め定められた位相に一致するか否かを予め定められた個所における前記プロファイルの信号強度に基づいて判定する判定手段と、
呼吸位相が予め定められた位相に一致しない頻度が予め定められた限度を超えたとき前記判定手段の判定個所を修正する修正手段と、
前記判定に基づいて撮影の対象から予め定められた呼吸位相におけるイメージングエコーを獲得するイメージングエコー獲得手段と、
前記イメージングエコーに基づいて画像を生成する生成手段と、
を具備することを特徴とする磁気共鳴撮影装置。
Navigator echo acquisition means for repeatedly acquiring navigator echo from the subject of shooting;
Creating means for creating a target profile based on the navigator echo;
Determining means for determining whether or not the breathing phase of the subject matches a predetermined phase based on the signal intensity of the profile at a predetermined location;
Correction means for correcting the determination location of the determination means when the frequency at which the respiratory phase does not match a predetermined phase exceeds a predetermined limit;
An imaging echo acquisition means for acquiring an imaging echo in a predetermined respiratory phase from an imaging target based on the determination;
Generating means for generating an image based on the imaging echo;
A magnetic resonance imaging apparatus comprising:
ナビゲータエコーおよびイメージングエコーを繰り返し獲得する獲得手段と、
前記ナビゲータエコーに基づいて対象のプロファイルを作成する作成手段と、
対象の呼吸位相が予め定められた位相に一致するか否かを予め定められた個所における前記プロファイルの信号強度に基づいて判定する判定手段と、
呼吸位相が予め定められた位相に一致しない頻度が予め定められた限度を超えたとき前記判定手段の判定個所を修正する修正手段と、
前記判定に基づいて前記イメージングエコーのうち予め定められた呼吸位相におけるイメージングエコーを採用する採用手段と、
前記採用されたイメージングエコーに基づいて画像を生成する生成手段と、
を具備することを特徴とする磁気共鳴撮影装置。
An acquisition means for repeatedly acquiring navigator echoes and imaging echoes;
Creating means for creating a target profile based on the navigator echo;
Determining means for determining whether or not the breathing phase of the subject matches a predetermined phase based on the signal intensity of the profile at a predetermined location;
Correction means for correcting the determination location of the determination means when the frequency at which the respiratory phase does not match a predetermined phase exceeds a predetermined limit;
Adopting means for adopting an imaging echo in a predetermined respiratory phase among the imaging echoes based on the determination;
Generating means for generating an image based on the employed imaging echo;
A magnetic resonance imaging apparatus comprising:
前記ナビゲータエコーはスピンエコーである、
ことを特徴とする請求項1または請求項2に記載の磁気共鳴撮影装置。
The navigator echo is a spin echo,
3. The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
前記イメージングエコーはグラディエントエコーである、
ことを特徴とする請求項1ないし請求項3のうちのいずれか1つに記載の磁気共鳴撮影装置。
The imaging echo is a gradient echo,
The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
前記プロファイルは体軸方向の1次元プロファイルである、ことを特徴とする請求項1ないし請求項4のうちのいずれか1つに記載の磁気共鳴撮影装置。The magnetic resonance imaging apparatus according to claim 1, wherein the profile is a one-dimensional profile in a body axis direction. 前記予め定められた位相は呼気位相である、
ことを特徴とする請求項1ないし請求項5のうちのいずれか1つに記載の磁気共鳴撮影装置。
The predetermined phase is an expiration phase;
6. The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
前記予め定められた位相は吸気位相である、
ことを特徴とする請求項1ないし請求項5のうちのいずれか1つに記載の磁気共鳴撮影装置。
The predetermined phase is an intake phase;
6. The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
前記頻度は呼吸位相が予め定められた位相に一致しない回数である、
ことを特徴とする請求項1ないし請求項7のうちのいずれか1つに記載の磁気共鳴撮影装置。
The frequency is the number of times the breathing phase does not match a predetermined phase.
The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
前記頻度は呼吸位相が予め定められた位相に一致しない時間である、
ことを特徴とする請求項1ないし請求項7のうちのいずれか1つに記載の磁気共鳴撮影装置。
The frequency is the time when the respiratory phase does not match a predetermined phase,
The magnetic resonance imaging apparatus according to claim 1, wherein the magnetic resonance imaging apparatus is a magnetic resonance imaging apparatus.
JP2002376451A 2002-12-26 2002-12-26 Magnetic resonance imaging device Expired - Fee Related JP3836424B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002376451A JP3836424B2 (en) 2002-12-26 2002-12-26 Magnetic resonance imaging device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002376451A JP3836424B2 (en) 2002-12-26 2002-12-26 Magnetic resonance imaging device

Publications (2)

Publication Number Publication Date
JP2004202043A JP2004202043A (en) 2004-07-22
JP3836424B2 true JP3836424B2 (en) 2006-10-25

Family

ID=32813919

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002376451A Expired - Fee Related JP3836424B2 (en) 2002-12-26 2002-12-26 Magnetic resonance imaging device

Country Status (1)

Country Link
JP (1) JP3836424B2 (en)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4711732B2 (en) * 2005-05-12 2011-06-29 株式会社日立メディコ Magnetic resonance imaging device
JP4896526B2 (en) 2006-01-11 2012-03-14 株式会社東芝 Magnetic resonance imaging system
JP4912808B2 (en) * 2006-09-22 2012-04-11 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Magnetic resonance imaging system
JP2008093098A (en) * 2006-10-10 2008-04-24 Ge Medical Systems Global Technology Co Llc Imaging apparatus
JP4956170B2 (en) * 2006-12-15 2012-06-20 株式会社日立メディコ Magnetic resonance imaging system
JP2008148918A (en) * 2006-12-18 2008-07-03 Ge Medical Systems Global Technology Co Llc Mri system and its control method
JP5097406B2 (en) * 2007-01-22 2012-12-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー MRI equipment
JP5454846B2 (en) * 2008-10-15 2014-03-26 株式会社東芝 Magnetic resonance imaging system
JP5498060B2 (en) * 2009-06-04 2014-05-21 株式会社東芝 Magnetic resonance imaging system
CN101721210B (en) 2008-10-15 2013-08-21 株式会社东芝 Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN101966081B (en) * 2010-08-05 2012-02-01 华东师范大学 Method for correcting breathing information and influencing resting state function magnetic resonance brain imaging
US9579041B2 (en) * 2012-05-24 2017-02-28 Toshiba Medical Systems Corporation Semi-automated non-contrast magnetic resonance angiography (MRA)
JP5555287B2 (en) * 2012-07-23 2014-07-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー MRI equipment
JP2013198795A (en) * 2013-07-08 2013-10-03 Toshiba Corp Magnetic resonance imaging apparatus
JP5889841B2 (en) * 2013-07-08 2016-03-22 株式会社東芝 Magnetic resonance imaging system
JP6105703B1 (en) 2015-09-29 2017-03-29 ミネベアミツミ株式会社 Biological information monitoring system
JP7539255B2 (en) * 2020-05-21 2024-08-23 富士フイルムヘルスケア株式会社 Magnetic Resonance Imaging

Also Published As

Publication number Publication date
JP2004202043A (en) 2004-07-22

Similar Documents

Publication Publication Date Title
JP3836424B2 (en) Magnetic resonance imaging device
US7570050B2 (en) Magnetic resonance imaging apparatus and method
US7209777B2 (en) Method and apparatus for automated tracking of non-linear vessel movement using MR imaging
EP1113288B1 (en) Respiratory displacement and velocity measurement using navigator MRI echo signals
US9301706B2 (en) Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same
JP3869337B2 (en) Magnetic resonance imaging device
JP3668816B1 (en) Magnetic resonance imaging system
US9684979B2 (en) MRI 3D cine imaging based on intersecting source and anchor slice data
JP5171373B2 (en) Magnetic resonance imaging apparatus and RF pulse applying method
JP2001000417A (en) Magnetic resonance imaging method for heart using multiple slabs and multiple windows
JP2010075573A (en) Magnetic resonance imaging instrument
JP4180936B2 (en) Magnetic resonance imaging device
KR20010090734A (en) Slice ordering method for breath-hold abdominal mr imaging
US20080064951A1 (en) Magnetic resonance imaging system
JP3701616B2 (en) Magnetic resonance imaging device
JP3930439B2 (en) Eddy current correction method and magnetic resonance imaging apparatus
WO2006109472A1 (en) Nuclear magnetic resonance imaging device and method
JP2004024669A (en) Magnetic resonance imaging apparatus
JP4349647B2 (en) Magnetic resonance imaging system
JP2006158762A (en) Mri apparatus
JP3753668B2 (en) RF pulse tuning device
US7595639B2 (en) System and method of angular elliptic centric view ordering for 3D MR acquisitions
JP3688782B2 (en) MRI equipment
JP4067938B2 (en) Magnetic resonance imaging device
JP2008148918A (en) Mri system and its control method

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20050126

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060612

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060620

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060726

R150 Certificate of patent or registration of utility model

Ref document number: 3836424

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090804

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100804

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100804

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100804

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110804

Year of fee payment: 5

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110804

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120804

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120804

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120804

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130804

Year of fee payment: 7

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130804

Year of fee payment: 7

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees