JPH11267111A - Mriシステムで拡散重み付き画像を成形する方法 - Google Patents

Mriシステムで拡散重み付き画像を成形する方法

Info

Publication number
JPH11267111A
JPH11267111A JP11026798A JP2679899A JPH11267111A JP H11267111 A JPH11267111 A JP H11267111A JP 11026798 A JP11026798 A JP 11026798A JP 2679899 A JP2679899 A JP 2679899A JP H11267111 A JPH11267111 A JP H11267111A
Authority
JP
Japan
Prior art keywords
pulse
gradient
pulses
echo
nmr
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.)
Granted
Application number
JP11026798A
Other languages
English (en)
Other versions
JP4236320B2 (ja
Inventor
C Mckinnon Gream
グリーム・シー・マッキンノン
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of JPH11267111A publication Critical patent/JPH11267111A/ja
Application granted granted Critical
Publication of JP4236320B2 publication Critical patent/JP4236320B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341Diffusion imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Vascular Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 【課題】 核磁気共鳴システムにおいて拡散重み付き画
像の画像アーティファクトを抑制する。 【解決手段】 核磁気共鳴システムにおいて拡散重み付
け高速スピン・エコー(FSE)パルス・シーケンス
(図4)を採用して、画像を形成するためのデータを取
得する。各々のFSEパルス・シーケンスの前に拡散重
み付けパルス・シーケンスが設けられ、画像アーティフ
ァクトを発生させるような拡散重み付きスピン磁化成分
を、FSEパルス・シーケンスが実行される前に印加さ
れる勾配パルス(370、376)とRFパルス(37
2)との組み合わせによって抑制する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明の分野は、核磁気共鳴
イメージング方法及びシステムである。より具体的に
は、本発明は、拡散重み付き(diffusion weighted)M
Rイメージングに用いられる改良された高速スピン・エ
コー・パルス・シーケンスに関する。
【0002】
【従来の技術】磁気モーメントを有しているあらゆる核
は、核が置かれている磁場の方向に沿って整列しようと
する。しかしながら、整列しようとするときに、核は磁
場の強度及び特定の核種についての特性(核の磁気回転
比γ)に依存する固有の角周波数(ラーモア周波数)で
磁場の方向の周りを歳差運動する。この現象を示す核
を、ここでは「スピン」と呼ぶ。
【0003】人体組織のような物体が一様の磁場(分極
磁場B0 )にさらされると、組織内のスピンの個々の磁
気モーメントは、この分極磁場に沿って整列しようとす
るが、各スピン固有のラーモア周波数においてランダム
な秩序で分極磁場の周りを歳差運動する。分極磁場の方
向には正味の磁気モーメントMz が発生するが、直交方
向の平面(x−y平面)すなわち横方向平面内にランダ
ムに配向した磁気成分は互いに打ち消し合う。しかしな
がら、x−y平面内に存在し且つラーモア周波数に近い
周波数を持つ磁場(励起磁場B1 )に物体すなわち組織
がさらされると、整列後の正味のモーメントMz はx−
y平面に向かって回転すなわち「傾いて」、ラーモア周
波数でx−y平面内で回転すなわちスピンする正味の横
磁気モーメントMt を発生する。この現象を表す実用的
な値は、励起信号B1 の停止後に励起されたスピンによ
って放出される信号中に存在する。この核磁気共鳴(N
MR)現象を利用する測定シーケンスには多様な種類の
ものがある。
【0004】NMRを利用して画像を形成するときに、
被検体の特定の位置からNMR信号を取得する技術が用
いられる。典型的には、イメージングしたい領域(関心
領域)が一連のNMR測定サイクルによって走査される
が、このNMR測定サイクルは、用いられる特定の局在
化方法によって異なる。結果として得られるNMR受信
信号のセットはディジタル化され処理されて、多くの周
知の再構成手法のうちの1つを用いることにより画像を
再構成する。このような走査(スキャン)を実行するた
めには、言うまでもなく、被検体内の特定の位置からの
NMR信号を引き出すことが必要である。これを達成す
るために、分極磁場B0 と同じ方向を有しているが、x
軸、y軸及びz軸に沿った勾配をそれぞれ有している磁
場(Gx、Gy 及びGz )を用いる。各々のNMRサイ
クル中のこれらの勾配の強度を制御することにより、ス
ピン励起の空間分布を制御すると共に、結果として得ら
れるNMR信号の位置を同定(または識別)することが
できる。
【0005】医用画像を形成するのに現状で用いられて
いる殆どのNMR走査は、必要なデータを取得するのに
何分もの時間を要求する。この走査時間を減少させる
と、患者のスループットが増大し、患者の安心感を高め
ると共に、モーション・アーティファクトの減少により
画質が向上するので、走査時間を減少させることは重要
な課題である。
【0006】短時間でNMR画像データを取得するとい
う概念は、1977年にピーター・マンスフィールド
(Peter Mansfield )によってエコー・プラナ・パルス
・シーケンスが提案されて以来、知られている(J.P
hys.誌、C.10:L55L58、1977年参
照)。標準的なパルス・シーケンスとは対照的に、エコ
ー・プラナ・パルス・シーケンスは、各々のRF励起パ
ルスについて1セットのNMR信号を発生させる。これ
らのNMR信号を別個に位相エンコーディングして、持
続時間にして20ミリ秒〜100ミリ秒の単一のパルス
・シーケンスで64のビューから成る走査の全体が取得
されるようにすることができる。エコー・プラナ・イメ
ージング(EPI)の利点は周知であり、このパルス・
シーケンスの多くの変形が、米国特許第4,678,9
96号、同第4,733,188号、同第4,716,
369号、同第4,355,282号、同第4,58
8,948号及び同第4,752,735号に開示され
ている。残念なことに、最新の高速勾配システムを用い
た場合でも、EPIパルス・シーケンスは、渦電流に依
存し且つ磁化率に誘起される画像の歪みを伴うという難
点を有している。
【0007】エコー・プラナ・イメージング法の一つの
変形に、高速取得緩和重み付け(Rapid Acquisition Re
laxation Enhanced ;RARE)シーケンスがあり、こ
のシーケンスは、Magnetic Resonanc
e in Medicine誌、第3巻、第823頁〜
第833頁(1986年)に所載のジェー・ヘニング
(J. Hennig )等による論文「RAREイメージング:
臨床MRのための高速イメージング法(RARE Imaging:
A Fast Imaging Method for Clinical MR )」に記載さ
れている。RAREシーケンスとEPIシーケンスとの
本質的な相違点は、エコー信号を形成する方式にある。
RAREシーケンスは、カー・パーセル・メイブーム・
ジル(Carr-Purcell-Meiboom-Gill )シーケンスによっ
て発生されるRF再収束(refocusing)され
たエコーを利用しており、他方、EPI法は、勾配呼び
戻し(gradient recalling)されたエコーを用いてい
る。この「高速スピン・エコー(FSE)」パルス・シ
ーケンスは一般的には、多数のk空間線を1回の励起で
取得するための問題のない手法と考えられている。例え
ば、この手法は、磁場の不均一性及び勾配のタイミング
誤差に対する感受性がエコー・プラナ・イメージングよ
りも遥かに小さい。更に、エコー・プラナ・イメージン
グの場合とは異なり、読み出し勾配が常に正であるの
で、勾配忠実度についての問題が少ない。
【0008】しかしながら、高速スピン・エコー・パル
ス・シーケンスには、正しく設定しないとかなりの画像
アーティファクトを形成する可能性のあるいくつかの重
要なパラメータが存在している。これらのパラメータに
は、無線周波数(RF)パルスの間隔及び位相関係、並
びに読み出し勾配パルスの面積がある。第一に、RF励
起パルスの中心と第1のRF再集束(refocusing)パル
スの中心との間の時間を、隣り合うRF再集束パルスの
中心同士の間の時間の2分の1にしなければならない。
第二に、エコー及びRF再集束パルスは、同じ位相角を
有していなければならない。これは通常、RF励起パル
スの位相をRF再集束パルスの位相に対して90°に設
定することにより達成されている。これらの要求条件に
関連して、RF励起パルスと第1のRF再集束パルスと
の間における読み出し勾配パルスの面積を、後続の再集
束パルスの各々の間における読み出し勾配パルスの面積
の2分の1に等しくしなければならないという事実があ
る。
【0009】従来の高速スピン・エコー(FSE)イメ
ージングの場合には、上述の重要なパラメータは、比較
的簡単な方式で制御することができる。しかしながら、
RFパルスとエコーとの間の位相制御の要求水準が達成
困難であるような多くのイメージング状況が存在してい
る。これらの状況の1つが、拡散重み付きイメージング
であり、この場合には、大きな勾配パルスが採用されて
おり、結果的に生ずる渦電流がより波及する。
【0010】拡散重み付きイメージングは、パルス・シ
ーケンスの開始時に1対の大きな勾配パルスを用いて、
取得されるNMR信号がスピンの動きに感受性を有する
ようにしている。一般的には、このようなイメージング
は、シングル・ショットEPIパルス・シーケンスを用
いて行われているが、深刻な画像の歪みが発生する可能
性がある。バッツ(Butts )等によってMagn.Re
son.Med.誌、第38巻、第741頁〜第749
頁(1997年)に提案された1つの解決法は、位相誤
差を補正するためにナヴィゲータ(navigator)信号を
用いたマルチ・ショットEPI取得を利用するものであ
る。
【0011】拡散重み付けFSEが試みられてはいる
が、パルス・シーケンスの開始時の大きなモーション・
エンコーディング勾配が渦電流を誘起し、これが励起パ
ルスと再集束パルスとの間の位相関係に干渉する。ノリ
ス(Norris)等は、Magn.Reson.Med.
誌、第27巻、第142頁〜第164頁(1992年)
に、2つのコヒーレンス経路を分離すると共にコヒーレ
ンス信号のうちの1つのみを用いるようなFSEパルス
・シーケンスにおいて位相を制御する方法を提案してい
る。このアプローチの1つの問題点は、補正しないと画
像に深刻なゴーストを発生させるエコー信号の振幅の強
い振動である。同様の着想がシック(Schick)によって
Magn.Reson.Med.誌、第38巻、第63
8頁〜第644頁(1997年)に提案されており、こ
こでは、エコー信号の振幅が増大している。アルソプ
(Alsop )は、Magn.Reson.Med.誌、第
38巻、第527頁〜第533頁(1997年)にこれ
らのエコー信号の振幅の振動を減少させる方法を開示し
ている。
【0012】
【発明が解決しようとする課題】本発明の一般的な目的
は、拡散重み付けFSE取得における画像アーティファ
クトを抑制することにある。所望の拡散重み付きスピン
磁化が、スピン磁化を縦軸に向かって回転させることに
より「保存」される。FSEパルス・シーケンスが、所
望の拡散重み付きスピン磁化を回転させて横方向平面に
戻して、一連のNMRエコー信号を発生させる。
【0013】
【課題を解決するための手段】本発明は、改良された高
速スピン・エコー・シーケンスに関し、より具体的に
は、拡散重み付きイメージングを行うための高速スピン
・エコー・パルス・シーケンスに関する。本発明は、従
来の高速スピン・エコー・パルス・シーケンスに対する
改良であり、この改良においては、拡散準備(diffusio
n preparation )パルス・シーケンスを実行することに
より拡散重み付きスピン横磁化を発生させ、第1の位相
外し(dephasing )勾配パルスを印加ずることによりス
ピン横磁化を位相外しし、RFパルスを印加することに
より位相外しされたスピン横磁化の1つの成分を縦軸に
沿って回転させ、第2の位相外し勾配パルスを印加する
ことによりスピン横磁化の残りの成分の位相外しさせ、
そして高速スピン・エコー・パルス・シーケンスを実行
することにより縦軸に沿って回転した上述の1つの成分
から一連のNMRエコー信号を発生させる。各々のNM
Rエコー信号取得に先立って拡散重み付きスピン磁化を
位相戻し(rephasing )するために、高速スピン・エコ
ー・パルス・シーケンスにおける各々のRF再集束パル
スの後に勾配パルスが形成される。
【0014】
【発明の実施の形態】先ず、図1について説明すると、
図1には、本発明を組み込んだ好ましいMRIシステム
の主要な構成要素が示されている。このシステムの動作
は、キーボード及び制御パネル102と表示装置104
とを含んでいるオペレータ・コンソール100から制御
される。コンソール100はリンク116を介して別個
の計算機システム107と連絡しており、計算機システ
ム107によりオペレータはスクリーン104上の画像
の形成及び表示を制御することが可能になる。計算機シ
ステム107は、バックプレーンを介して互いに連絡し
ている多数のモジュールを含んでいる。これらのモジュ
ールは、画像プロセッサ・モジュール106と、CPU
モジュール108と、画像データを記憶するフレーム・
バッファとして当業界で知られているメモリ・モジュー
ル113とを含んでいる。計算機システム107は、画
像データ及びプログラムを記憶するためのディスク記憶
装置111及びテープ・ドライブ112に結合されてお
り、また高速シリアル・リンク115を介して別個のシ
ステム制御部122と連絡する。
【0015】システム制御部122は、バックプレーン
118によってまとめて接続された一組のモジュールを
含んでいる。これらのモジュールは、CPUモジュール
119及びパルス発生器モジュール121を含んでい
る。パルス発生器モジュール121は、シリアル・リン
ク125を介してオペレータ・コンソール100に接続
されている。リンク125を介して、システム制御部1
22は実行されるべき走査シーケンスを指示する命令
(コマンド)をオペレータから受け取る。パルス発生器
モジュール121は、システムの構成要素を動作させ
て、所望の走査シーケンスを実行する。パルス発生器モ
ジュール121は、発生すべきRFパルスのタイミン
グ、振幅及び形状、並びにデータ取得ウィンドウのタイ
ミング及び長さを指示するデータを発生する。パルス発
生器モジュール121は、一組の勾配増幅器127に接
続されており、走査中に発生される勾配パルスのタイミ
ング及び形状を指示する。パルス発生器モジュール12
1はまた、患者に接続された多数の異なるセンサからの
信号、例えば電極からの心電図(ECG)信号またはベ
ローズからの呼吸信号を受信する生理学データ取得制御
装置129から、患者のデータを受信する。最後に、パ
ルス発生器モジュール121は走査室インタフェイス回
路133にも接続されており、走査室インタフェイス回
路133は、患者及び磁石システムの状態に関連した様
々なセンサからの信号を受信する。走査室インタフェイ
ス回路133を介して、患者位置決めシステム134も
また、走査のための望ましい位置に患者を移動させるた
めの命令を受信する。
【0016】パルス発生器モジュール121によって発
生された勾配波形は、Gx 増幅器とGy 増幅器とGz
幅器とで構成されている勾配増幅器システム127に印
加される。各々の勾配増幅器は、全体的に参照番号13
9で示す勾配コイル・アセンブリ内の対応する勾配コイ
ルを励起して、取得される信号を位置エンコーディング
するのに用いられる磁場勾配を発生する。勾配コイル・
アセンブリ139は、分極磁石140と全身用RFコイ
ル152とを含んでいる磁石アセンブリ141の一部を
形成している。システム制御部122内のトランシーバ
・モジュール150がパルスを発生し、これらのパルス
はRF増幅器151によって増幅されて、送受信切換え
(T/R)スイッチ154を介してRFコイル152に
結合される。患者の内部の励起した核によって放出され
る結果として生ずる信号は、同じRFコイル152によ
って検知して、T/Rスイッチ154を介して前置増幅
器153に結合することが出来る。増幅されたNMR信
号は、トランシーバ150の受信器部において復調され
濾波されて、ディジタル化される。T/Rスイッチ15
4は、パルス発生器モジュール121からの信号によっ
て制御されて、送信モード時にはRF増幅器151をコ
イル152に電気的に接続し、受信モード時には前置増
幅器153をコイル152に電気的に接続する。T/R
スイッチ154はまた、送信モードまたは受信モードの
いずれの場合にも、分離型RFコイル(例えば、頭部コ
イルまたは表面コイル)を用いることを可能にする。
【0017】RFコイル152によって検出されたNM
R信号は、トランシーバ・モジュール150によってデ
ィジタル化されて、システム制御部122内のメモリ・
モジュール160へ転送される。走査が完了してデータ
配列全体がメモリ・モジュール160内に取得されたと
きに、アレイ・プロセッサ161が動作して、このデー
タをフーリエ変換して画像データの配列とする。この画
像データは、シリアル・リンク115を介して計算機シ
ステム107へ伝送され、ここで、ディスク・メモリ1
11に記憶される。オペレータ・コンソール100から
受信された命令に応答して、この画像データをテープ・
ドライブ112に保管してもよいし、または画像プロセ
ッサ106によって更に処理してオペレータ・コンソー
ル100へ伝送すると共に表示装置104に表示しても
よい。
【0018】図1及び図2について具体的に説明する。
トランシーバ150は、電力増幅器151を介してコイ
ル152AにRF励起磁場B1 を発生させ、またその結
果得られるコイル152B内に誘導された信号を受信す
る。上述のように、コイル152A及びコイル152B
は、図2に示すような分離型であってもよいし、または
図1に示すような単一の全身用コイルであってもよい。
RF励起磁場の基本周波数すなわち搬送周波数は、周波
数合成器200の制御下で発生される。周波数合成器2
00は、CPUモジュール119及びパルス発生器モジ
ュール121から一組のディジタル信号(CF)を受信
する。これらのディジタル信号は、出力201に発生さ
れるRF搬送波信号の周波数及び位相を表す。命令に従
って発生されたRF搬送波は、変調器及びアップ・コン
バータ202に印加され、ここで、その振幅は、やはり
パルス発生器モジュール121から受信された信号R
(t)に応答して変調される。信号R(t)は、発生す
べきRF励起パルスの包絡線を画定しており、記憶され
た一連のディジタル値を順次読み出すことによりモジュ
ール121内で発生される。これらの記憶されたディジ
タル値はまた、オペレータ・コンソール100から変更
可能であり、任意の所望のRFパルス包絡線を発生させ
ることができる。
【0019】出力205に発生されたRF励起パルスの
振幅は、バックプレーン118からディジタル命令TA
を受信する励起信号減衰器回路206によって減衰され
る。減衰されたRF励起パルスは、RFコイル152A
を駆動する電力増幅器151へ印加される。トランシー
バのこの部分に関する更なる詳細については、ここに参
考として引用する米国特許第4,952,877号に記
載されている。
【0020】図1及び図2について説明を続ける。被検
体によって発生された信号は、受信器コイル152Bに
よって検知されて、前置増幅器153を介して受信信号
減衰器207の入力へ印加される。受信信号減衰器20
7は、バックプレーン118から受信されたディジタル
減衰信号(RA)によって決定されている量だけ信号を
更に増幅する。
【0021】受信信号はラーモア周波数またはそれに近
い周波数であり、この高周波信号はダウン・コンバータ
208によって次の2段階の処理で下降変換(ダウン・
コンバート)される。すなわち、ダウン・コンバータ2
08は、先ず、NMR信号を線201の搬送波信号と混
合し、次いでその結果の差信号を線204の2.5MH
zの基準信号と混合する。下降変換されたNMR信号
は、アナログ−ディジタル(A/D)変換器209の入
力へ印加される。A/D変換器209は、アナログ信号
をサンプリングしディジタル化して、これをディジタル
検出器及び信号プロセッサ210へ印加する。ディジタ
ル検出器及び信号プロセッサ210は、受信した信号に
対応する16ビットの同相(I)値及び16ビットの直
角位相(Q)値を発生する。その結果得られる、受信信
号のディジタル化されたI値及びQ値のストリームは、
バックプレーン118を介してメモリ・モジュール16
0へ出力され、ここで画像を再構成するのに用いられ
る。
【0022】2.5MHzの基準信号、250kHzの
サンプリング信号、並びに5MHz、10MHz及び6
0MHzの基準信号が、共通の20MHzマスタ・クロ
ック信号から基準周波数発生器203によって発生され
る。受信器に関する更なる詳細については、ここに参考
として引用する米国特許第4,992,736号に記載
されている。
【0023】次に、図3について説明する。図3には、
2DFT−RAREシーケンスとして知られている従来
の高速スピン・エコーNMRパルス・シーケンスが示さ
れている。簡明にするために、図3には4つのエコー信
号301〜304のみが示されているが、各々の「ショ
ット」毎に、より多くのエコー信号が発生され取得され
ることが理解されよう。これらのNMRエコー信号は、
z スライス選択勾配パルス306の存在下で印加され
る90°RF励起パルス305によって発生され、患者
のスライス内の縦磁化を回転させて横磁化を形成する。
この横磁化は、各々の選択的180°RF再集束パルス
307によって再収束されて、NMRスピン・エコー信
号301〜304を発生させる。これらのNMRスピン
・エコー信号301〜304は、Gx 読み出し勾配パル
ス308の存在下で取得される。各々のNMRスピン・
エコー信号301〜304は、それぞれのGy 位相エン
コーディング・パルス309〜313によって別々に位
相エンコーディングされる。各々の位相エンコーディン
グ・パルスの大きさは異なっていて、256個の値に段
階的に変化して、1回の完全な走査中に256個の別々
のビューが取得される。これにより、y方向に256個
の別々の画素(ピクセル)を有している画像を再構成す
ることが可能になる。各々のNMRスピン・エコー信号
は、各々の信号の256個のサンプルをディジタル化す
ることにより取得される。その結果、1つの画像の走査
の完了時には、16ショット分(256/16=16)
の図3のパルス・シーケンスが実行されており、256
×256要素の複素数の配列が取得されている。この高
速スピン・エコー・パルス・シーケンスの変形は、クラ
ッシャ(crusher )勾配パルス316を含んでいる。こ
れらのクラッシャ勾配パルス316は、等しい面積を有
しており、各々の再集束RFパルス307の直前及び直
後に印加される。
【0024】この画像データ配列に2次元(2D)フー
リエ変換を施し、次いで、結果として得られる各々の複
素要素の絶対値を算出することにより、1枚の画像が再
構成される。このようにして、各々の画素(ピクセル)
の輝度が変換後の配列の対応する要素の大きさによって
決定されるような256×256画素の画像が形成され
る。
【0025】図3の従来の高速スピン・エコー・パルス
・シーケンスでは、RF励起パルス305は90°のフ
リップ角を有しており、RF再集束パルスは名目上18
0°のフリップ角である。しかしながら、実用上は、9
0°から180°にわたる他のフリップ角が用いられて
いる。加えて、「90°x 」及び「180°y 」という
表記は、RFパルス305と再集束パルス307との間
のRF位相角が90°であることを示している。この直
交関係は、この位相関係が維持されているときに発生さ
れるエコー信号がすべて同じ符号を有するので、重要で
ある。対照的に、RFパルスがすべて同相にあるときに
は、発生されるエコー信号は符号が交番する。従って、
RF励起パルス305と再集束パルス307との間の位
相角が0°である場合、180°の再集束パルスが用い
られているときにはエコー信号の振幅が振動(oscillat
e )する。このように、従来の高速スピン・エコー・パ
ルス・シーケンスは、RF励起パルス305と再集束R
Fパルス307との間の直交位相関係の適正な維持に敏
感である。
【0026】従来の高速スピン・エコー・パルス・シー
ケンスのこの位相感受性を図5に示す。図5には、「同
相」のRF励起パルスがエコー信号の振幅に及ぼす影響
が示されている。曲線330は、180°RF再集束パ
ルスが用いられるときのエコー信号の振幅の符号の交番
の様子を示している。曲線332及び334は、RF再
集束パルスのフリップ角をそれぞれ135°及び90°
に減少させたときのエコー信号の振幅を示している。エ
コー信号の振幅は、エコー番号が増大するにつれて急速
に0に向かって減衰している。これは、スピン・エコー
と誘導されたエコーとの間の相殺に起因している。ま
た、信号の振幅の振動の時間が、360°を集束パルス
のフリップ角で除算したものに等しくなることに留意さ
れたい。
【0027】次に、図4について説明する。本発明を実
施するのに用いられるパルス・シーケンスの好ましい実
施例は、標準的なFSEシーケンスの前に縦磁化の拡散
準備(diffusion preparation )を含む。この拡散準備
は、z勾配軸、y勾配軸及びx勾配軸にそれぞれ沿った
大きなモーション・エンコーディング勾配パルス35
0、352及び354によって行われる。これらのモー
ション・エンコーディング勾配パルス350、352及
び354は、選択性90°RF励起パルス356と選択
性180°RF再集束パルス358とで構成されている
スピン・エコー・シーケンスの一部である。RF励起パ
ルス356は、スライス選択勾配パルス360によって
決定されるスライス位置において横磁化を発生させ、R
Fパルス358は第2のスライス選択勾配パルス362
の存在下でこの横磁化を再収束させる。このようにし
て、動きに感受性のある横磁化は、破線364によって
示される時刻に再収束する。
【0028】高速スピン・エコー(FSE)パルス・シ
ーケンスは、スライス選択勾配パルス306の存在下で
発生される90°RF励起パルス305によって開始す
る。前述のように、180°RF再集束パルス307及
び位相エンコーディング勾配パルス309が適正な時刻
に印加されてNMRエコー信号(図4には示されていな
い)を発生させ、これらのNMRエコー信号は読み出し
勾配パルス308の印加中に取得される。図4には、F
SEパルス・シーケンスの一部しか示していないが、当
業者には、ショット内に所望されるだけの数のNMRエ
コー信号についてこのパルス・シーケンスが続行される
ことが理解されよう。拡散重み付きイメージングでは、
ETL=40〜72のFSEショットが典型的である。
【0029】本発明は、上述のパルス・シーケンスに対
する改良であり、読み出し勾配Gxに沿って位相外し勾
配パルス370を印加して、拡散重み付け過程中に発生
される横磁化を位相外しする。この意味する所は、各々
のボクセル(voxel )内で、横磁化を有すスピンの位相
が、横方向平面(すなわち、x−y平面)内で360°
にわたって広げられることである。次いで、スライス選
択勾配374の存在下で90°パルス372を印加し
て、y軸に沿って配向していたスピン磁化を回転させて
縦方向のz軸に戻す。次いで、Gy 位相エンコーディン
グ軸に沿って第2の位相外し勾配パルス376を印加し
て、横方向平面に残存する残留横磁化を位相外しする。
90°RF励起パルス305は、RFパルス372によ
ってz軸に沿って「保存」されているアーティファクト
無しのスピン磁化を回転させて横方向平面に戻し、この
平面において、カー・パーセル・メイブーム・ジル(Ca
rr-Purcell-Meiboom-Gill )シーケンスでNMRエコー
信号を発生させる。
【0030】横磁化の拡散重み付き用成分のみを再収束
させるために、各々のRF再集束パルス307の直後に
x 読み出し軸に沿って位相戻し(rephasing )勾配パ
ルス380が印加される。この位相戻し勾配パルス38
0は、拡散準備シーケンス中に回復した可能性のあるT
1重み付き信号成分を抑制する。加えて、位相戻し勾配
パルス380は、位相外し勾配パルス370と同じ寸法
とする。各々の拡散重み付きエコー信号が取得された後
に、同じ寸法の巻き戻し(rewinder)勾配パルス382
が印加される。巻き戻し勾配382は、隣り合う再集束
パルス307の間での勾配面積を励起パルス305と第
1の再集束パルス307との間での勾配面積の2倍とす
ることを要求するCPMG条件が満たされることを保証
するためにこの位置に存在する。
【0031】本発明から逸脱することなく上記の好まし
い実施例の多くの変形が可能である。位相外し勾配パル
ス370の極性を反転させることもできるし、第2の位
相外し勾配パルス376を任意の方向に印加することも
できる。また、勾配パルス370、380及び382は
これらの3つの勾配のすべてが同じ方向に印加される限
りにおいて任意の方向に印加することもできる。180
°よりも小さいフリップ角を有する再集束パルスを用い
る等、当業者に周知の他の変形もまた可能である。
【図面の簡単な説明】
【図1】本発明を採用しているMRIシステムのブロッ
ク図である。
【図2】図1のMRIシステムの一部を形成しているト
ランシーバの電気ブロック図である。
【図3】従来の高速スピン・エコー・パルス・シーケン
スを示すグラフである。
【図4】本発明を採用している改良された高速スピン・
エコー・パルス・シーケンスを示すグラフである。
【図5】従来の高速スピン・エコー・パルス・シーケン
スにおけるエコー信号の振幅の変化を示すグラフであ
る。
【符号の説明】
100 オペレータ・コンソール 102 キーボード及び制御パネル 104 表示装置 106 画像プロセッサ 107 計算機システム 108、119 CPUモジュール 111 ディスク記憶装置 112 テープ・ドライブ 113、160 メモリ・モジュール 115 高速シリアル・リンク 116 リンク 118 バックプレーン 121 パルス発生器モジュール 122 システム制御部 125 シリアル・リンク 127 勾配増幅器 129 生理学データ取得制御装置 133 走査室インタフェイス回路 134 患者位置決めシステム 139 勾配コイル・アセンブリ 140 分極磁石 141 磁石アセンブリ 150 送受信器 151 RF増幅器 152 全身型RFコイル 152A、152B 分離型コイル 153 前置増幅器 154 送受信切換えスイッチ 161 アレイ・プロセッサ 200 周波数合成器 201、204、205、212 出力線 202 変調器及びアップ・コンバータ 203 基準周波数発生器 206 励起信号減衰器回路 207 受信信号減衰器 208 ダウン・コンバータ 209 アナログ−ディジタル変換器 210 ディジタル検出器及び信号プロセッサ

Claims (11)

    【特許請求の範囲】
  1. 【請求項1】 磁気共鳴イメージング(MRI)システ
    ムを用いて拡散重み付き画像を形成する方法であって、 a)RF励起パルスにより横磁化を生じさせ、バイポー
    ラ勾配パルスの印加により該横磁化を拡散重み付けする
    拡散重み付けパルス・シーケンスを実行する工程と、 b)第1の勾配パルスの印加により前記の拡散重み付き
    横磁化を位相外しする工程と、 c)前記拡散重み付き横磁化の一成分を縦軸に向かって
    回転させるRFパルスを印加することにより、該成分を
    保存する工程と、 d)第2の勾配パルスを印加することにより、横磁化を
    位相外しする工程と、 e)RF励起パルスにより前記の保存された拡散重み付
    き磁化を横方向平面に向かって回転させ、一連の核磁気
    共鳴エコー信号を対応する一連のRF再集束パルスによ
    り発生させる高速スピン・エコー・パルス・シーケンス
    を実行する工程と、 f)取得された核磁気共鳴(NMR)エコー信号から画
    像を再構成する工程とを含んでいる前記方法。
  2. 【請求項2】 前記バイポーラ勾配パルスが3つの軸に
    沿って印加される請求項1に記載の方法。
  3. 【請求項3】 前記高速スピン・エコー・パルス・シー
    ケンスの実行中に各々のRF再集束パルスの後で且つ対
    応するNMRエコー信号の取得の前に第3の勾配パルス
    が印加され、該第3の勾配パルスの各々の寸法は前記第
    1の勾配パルスの寸法と実質的に同じである請求項1に
    記載の方法。
  4. 【請求項4】 前記高速スピン・エコー・パルス・シー
    ケンスの実行中に各々のNMRエコー信号の取得の後に
    第4の勾配パルスが印加され、該第4の勾配パルスの各
    々の寸法は前記第3の勾配パルスの寸法と実質的に同じ
    であるが、極性が前記第3の勾配パルスの極性と反対で
    ある請求項3に記載の方法。
  5. 【請求項5】 各々のNMRエコー信号は、1つの勾配
    軸に沿って方向付けされた読み出し勾配パルスの印加中
    に取得され、第2の勾配軸に沿って方向付けされた位相
    エンコーディング勾配パルスが、各々のNMRエコー信
    号の前に発生されている請求項4に記載の方法。
  6. 【請求項6】 前記バイポーラ勾配パルスが前記勾配軸
    の各々に沿って印加される請求項5に記載の方法。
  7. 【請求項7】 前記第1、第3及び第4の勾配パルス
    は、読み出し勾配軸に沿って方向付けされている請求項
    5に方法。
  8. 【請求項8】 前記第2の勾配パルスは、位相エンコー
    ディング勾配軸に沿って方向付けされている請求項7に
    方法。
  9. 【請求項9】 前記RF励起パルスの各々及び前記RF
    再集束パルスの各々の印加中にスライス選択勾配パルス
    が印加される請求項8に方法。
  10. 【請求項10】 前記バイポーラ勾配パルスは同じ極性
    の2つの大きな勾配パルスで構成されており、該2つの
    大きな勾配パルスの間にRFパルスが発生される請求項
    1に記載の方法。
  11. 【請求項11】 画像を再構成するのに用いられる更な
    るNMRエコー信号を取得するために前記の工程(a)
    〜工程(e)が繰り返される請求項1に記載の方法。
JP02679899A 1998-02-13 1999-02-04 Mriシステムで拡散重み付き画像を形成する方法 Expired - Fee Related JP4236320B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/023572 1998-02-13
US09/023,572 US6078176A (en) 1996-11-08 1998-02-13 Fast spin echo pulse sequence for diffusion weighted imaging

Publications (2)

Publication Number Publication Date
JPH11267111A true JPH11267111A (ja) 1999-10-05
JP4236320B2 JP4236320B2 (ja) 2009-03-11

Family

ID=21815929

Family Applications (1)

Application Number Title Priority Date Filing Date
JP02679899A Expired - Fee Related JP4236320B2 (ja) 1998-02-13 1999-02-04 Mriシステムで拡散重み付き画像を形成する方法

Country Status (5)

Country Link
US (1) US6078176A (ja)
JP (1) JP4236320B2 (ja)
CN (1) CN1138156C (ja)
DE (1) DE19905720B4 (ja)
IL (1) IL128354A0 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010527740A (ja) * 2007-05-31 2010-08-19 シーアール ディベロップメント アーベー 拡散磁気共鳴画像化のための方法とシステム
KR20130086321A (ko) * 2012-01-24 2013-08-01 지멘스 메디컬 솔루션즈 유에스에이, 인크. 모션 정정 mr 확산 이미징을 위한 시스템
CN113811783A (zh) * 2019-04-26 2021-12-17 Cr发展公司 执行弥散加权磁共振测量的方法

Families Citing this family (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100331449B1 (ko) * 1999-11-12 2002-04-09 윤종용 패스트 스핀 에코 영상법에서 위상 부호화 경사 자계펄스의 위상 오류 보정방법
DE10021496C2 (de) 2000-05-03 2002-12-05 Siemens Ag Verfahren und Kernspintomographiegerät zur Durchführung dieses Verfahrens zur Erzeugung eines Kernspinbildes eines Objektes mit verschiedenen Pulssequenzen im k-Raum
JP3875479B2 (ja) * 2000-10-20 2007-01-31 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴撮影装置
JP2004513720A (ja) * 2000-11-14 2004-05-13 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 傾斜磁場の非線形性の補正を行う磁気共鳴撮像方法
US6546275B2 (en) 2001-06-25 2003-04-08 Wisconsin Alumni Research Foundation Determination of the arterial input function in dynamic contrast-enhanced MRI
EP1459090A1 (en) * 2001-12-19 2004-09-22 Koninklijke Philips Electronics N.V. Diffusion weighted multiple spin echo (rare) sequence with periodically amplitude-modulated crusher gradients
US6859203B2 (en) * 2002-05-15 2005-02-22 Koninklijke Philips Electronics N.V. Sweeping real-time single point fiber
US6750651B2 (en) * 2002-07-03 2004-06-15 The Board Of Trustees Of The Leland Stanford Junior University Fat suppression in MRI using oscillating steady-state free precession
DE10250922B4 (de) * 2002-10-31 2008-01-31 Siemens Ag Verfahren zur Ermittlung des ADC-Koeffizienten in der Diffusionsgewichteten Magnetresonanz-Bildgebung bei Verwendung von Steady-State-Sequenzen
CN101273276A (zh) * 2005-09-22 2008-09-24 威斯康星校友研究基金会 利用高度约束的背投来重构跳动着的心脏的图像
JP5037075B2 (ja) * 2005-12-22 2012-09-26 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP4249215B2 (ja) * 2006-10-06 2009-04-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP5022696B2 (ja) * 2006-12-22 2012-09-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
US7804299B2 (en) * 2007-05-04 2010-09-28 Wisconsin Alumni Research Foundation Diffusion weighted preparatory sequence for magnetic resonance imaging pulse sequence
WO2008147923A1 (en) * 2007-05-22 2008-12-04 Imaging Biometrics Method for detecting tumor cell invasion using short diffusion times
US8320647B2 (en) * 2007-11-20 2012-11-27 Olea Medical Method and system for processing multiple series of biological images obtained from a patient
CN101236241B (zh) * 2008-02-28 2010-12-08 上海师范大学 一种对扩散张量磁共振图像进行恢复处理的方法
JP5280089B2 (ja) * 2008-04-23 2013-09-04 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
US8875261B2 (en) * 2008-10-22 2014-10-28 International Business Machines Corporation Rules driven multiple passwords
DE102009014498B4 (de) * 2009-03-23 2011-03-24 Siemens Aktiengesellschaft Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquisitionstechnik
JP5633899B2 (ja) * 2009-06-30 2014-12-03 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
EP2585845A1 (en) * 2010-06-24 2013-05-01 Ramot at Tel Aviv University Ltd. Magnetic resonance analysis using a plurality of pairs of bipolar gradient pulses
US9612301B2 (en) * 2011-02-14 2017-04-04 Duke University High-throughput and motion insensitive MRI accelerated with multi-echo planar acquisition and related systems
CN108318843B (zh) * 2011-09-22 2021-11-16 东芝医疗系统株式会社 磁共振成像装置以及磁共振成像方法
JP6139119B2 (ja) * 2012-01-13 2017-05-31 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
DE102013201671B3 (de) 2013-02-01 2014-07-31 Siemens Aktiengesellschaft Verfahren zum Erfassen von MR-Daten und zur Bestimmung eines B1-Magnetfelds sowie entsprechend ausgestaltete Magnetresonanzanlage
DE102013219754B4 (de) * 2013-09-30 2018-11-22 Siemens Healthcare Gmbh Optimierung der Geräuschentwicklung einer 3D-Gradientenchosequenz in einer Magnetresonanzanlage
US9594144B2 (en) 2014-04-23 2017-03-14 General Electric Company Low-noise magnetic resonance imaging using low harmonic pulse sequences
US10921410B2 (en) * 2015-02-03 2021-02-16 Koninklijke Philips N.V. Method and system for susceptibility weighted magnetic resonance imaging
DE102015223658B4 (de) * 2015-11-30 2017-08-17 Siemens Healthcare Gmbh Verfahren zum Erfassen von Magnetresonanz-Signalen eines Untersuchungsobjekts
CN108020796B (zh) * 2016-10-31 2021-10-29 西门子(深圳)磁共振有限公司 一种磁共振扩散加权成像方法和装置
EP3364207A1 (en) * 2017-02-17 2018-08-22 Koninklijke Philips N.V. Phase error correction for bipolar read out gradients
CN107245041A (zh) * 2017-05-25 2017-10-13 北京万鹏朗格医药科技有限公司 一种精氨酸布洛芬的制备方法
CN109085522B (zh) * 2018-07-02 2021-07-09 上海东软医疗科技有限公司 一种磁共振弥散加权成像和波谱信号的采集方法和装置
CN112684393B (zh) * 2020-12-28 2021-11-23 电子科技大学 一种增强磁共振谱仪动态范围的方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3193467B2 (ja) * 1992-07-31 2001-07-30 ジーイー横河メディカルシステム株式会社 Mri装置
US5317262A (en) * 1992-08-05 1994-05-31 The United States Of America As Represented By The Department Of Health And Human Services Single shot magnetic resonance method to measure diffusion, flow and/or motion
JPH06133944A (ja) * 1992-10-26 1994-05-17 Shimadzu Corp Mrイメージング装置
US5492123A (en) * 1993-08-13 1996-02-20 Siemens Medical Systems, Inc. Diffusion weighted magnetic resonance imaging
US5684400A (en) * 1995-06-12 1997-11-04 Ge Yokogawa Medical Systems, Limited Diffusion sensitizing imaging method, dynamic imaging method, and MRI apparatus
US5864233A (en) * 1997-02-22 1999-01-26 General Electric Company Method to reduce eddy current effects in diffusion-weighted echo planar imaging

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010527740A (ja) * 2007-05-31 2010-08-19 シーアール ディベロップメント アーベー 拡散磁気共鳴画像化のための方法とシステム
US8565854B2 (en) 2007-05-31 2013-10-22 Cr Development Ab Method and system for diffusion magnetic resonance imaging
KR20130086321A (ko) * 2012-01-24 2013-08-01 지멘스 메디컬 솔루션즈 유에스에이, 인크. 모션 정정 mr 확산 이미징을 위한 시스템
CN113811783A (zh) * 2019-04-26 2021-12-17 Cr发展公司 执行弥散加权磁共振测量的方法
CN113811783B (zh) * 2019-04-26 2024-02-13 兰德姆沃克成像有限公司 执行弥散加权磁共振测量的方法

Also Published As

Publication number Publication date
IL128354A0 (en) 2000-01-31
DE19905720B4 (de) 2010-10-14
JP4236320B2 (ja) 2009-03-11
CN1138156C (zh) 2004-02-11
CN1234508A (zh) 1999-11-10
DE19905720A1 (de) 1999-09-30
US6078176A (en) 2000-06-20

Similar Documents

Publication Publication Date Title
JP4236320B2 (ja) Mriシステムで拡散重み付き画像を形成する方法
US7495439B2 (en) MRI method for reducing artifacts using RF pulse at offset frequency
US5256967A (en) Fast NMR image acquisition with spectrally selective inversion pulses
US6265873B1 (en) Non-CPMG fast spin echo MRI method
US5378985A (en) Fast spin echo prescan for MRI system
JP4229487B2 (ja) マクスウェル項誤差を補正する方法
JP3381178B2 (ja) 磁気共鳴画像システム
US5422576A (en) Magnetic resonance angiography using fast spin echo pulse sequence
US8072215B2 (en) Magnetic resonance device and method
US5347216A (en) Fast NMR image acquisition with spectrally selective inversion pulse
JP3445300B2 (ja) 高速スピンエコーnmrパルス系列での勾配モーメントナリング法
JP4395516B2 (ja) 脂肪抑制を用いる血管の等方性画像形成
JPH0654828A (ja) Nmrシステム
US9211082B2 (en) Method for magnetic resonance imaging using saturation harmonic induced rotary saturation
US5315249A (en) Stabilized fast spin echo NMR pulse sequence
US5339035A (en) MR imaging with rectangular magnetization transfer pulse
US6380738B1 (en) Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system
US6329821B1 (en) Method and apparatus to compensate for image artifacts caused by magnet vibration in an MR imaging system
US20110115487A1 (en) Method and magnetic resonance system for imaging particles
AU754136B2 (en) Recovery of signal void arising from field inhomogeneities in magnetic resonance echo planar imaging
US5291891A (en) Monitoring body functions using fast NMR pulse sequences
US5541513A (en) MRI center point artifact elimination using realtime receiver phase control
US6239597B1 (en) Method and apparatus for rapid T2 weighted MR image acquisition
US5233302A (en) Masking motion ghost artifacts in NMR images
JPH0263435A (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20080402

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080624

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080804

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081007

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20081010

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: 20081118

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20081216

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20111226

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20121226

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20121226

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20131226

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees