JPH07241280A - Data collecting method and imaging method of mri device and mri device - Google Patents

Data collecting method and imaging method of mri device and mri device

Info

Publication number
JPH07241280A
JPH07241280A JP6033745A JP3374594A JPH07241280A JP H07241280 A JPH07241280 A JP H07241280A JP 6033745 A JP6033745 A JP 6033745A JP 3374594 A JP3374594 A JP 3374594A JP H07241280 A JPH07241280 A JP H07241280A
Authority
JP
Japan
Prior art keywords
data
sampling
echo
time
center
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
JP6033745A
Other languages
Japanese (ja)
Other versions
JP3447099B2 (en
Inventor
Yoshikazu Ikezaki
吉和 池崎
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 Healthcare Japan Corp
Original Assignee
GE Yokogawa Medical System Ltd
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 Yokogawa Medical System Ltd filed Critical GE Yokogawa Medical System Ltd
Priority to JP03374594A priority Critical patent/JP3447099B2/en
Publication of JPH07241280A publication Critical patent/JPH07241280A/en
Application granted granted Critical
Publication of JP3447099B2 publication Critical patent/JP3447099B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

PURPOSE:To decrease the artifacts and ringing artifacts by phase distortions even if phase correction is not executed by setting the number of sampling before an echo center in time at a specific value and sampling echo signals, thereby collecting data. CONSTITUTION:Data are collected by sampling the echo signals at the sampling point of N-pieces of the number (m) of sampling before the echo center in time and the number (p) of sampling after the echo center in time at the time of collecting the data with an MRI device 100. At this time, m<p, NX20%<=mNX35%, m+p=N are set. A real number component Re and a virtual number component Im are obtd. by subjecting the prescribed data to two-dimensional Fourier transform and an image of the absolute value Z=(Re<2>+ Im<2>)<1/2> is formed in image reconstitution with a calculator 7. As a result, the artifacts by the phase distortion are decreased without making phase correction. In addition, dealing with the short echo time is made possible without losing the advantage of an asymmetrical echo method.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は、MRI(Magnetic R
esonance Imaging)装置におけるデータ収集方法および
イメージング方法およびMRI装置に関し、さらに詳し
くは、画像再構成時間を短縮できるようにしたMRI装
置におけるデータ収集方法およびイメージング方法およ
びMRI装置に関する。
BACKGROUND OF THE INVENTION This invention is applied to MRI (Magnetic R
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a data acquisition method, an imaging method, and an MRI apparatus in an esonance imaging apparatus, and more specifically, to a data acquisition method, an imaging method, and an MRI apparatus in an MRI apparatus that can reduce image reconstruction time.

【0002】[0002]

【従来の技術】図8は、MRI装置における通常のフル
エコー法によるデータ収集方法の説明図である。すなわ
ち、k−スペースの周波数軸方向のデータ数をNとする
とき、エコー中心(プロトン励起用のパルス系列SGに
おける励起パルスからエコー時間TE後の時間位置)よ
り時間的に前のサンプリング数N/2およびエコー中心
より時間的に後のサンプリング数N/2のN個のサンプ
リング点でエコー信号をサンプリングし、サンプリング
データを収集する。従って、図9に示すように、収集し
たサンプリングデータは、k−スペースSを埋めつく
す。画像再構成では、前記k−スペースSのデータに対
し2次元フーリエ変換を行い、実数成分Reと虚数成分
Imを得て、絶対値Z=√{Re2+Im2}の画像を生
成する。
2. Description of the Related Art FIG. 8 is an explanatory view of a data acquisition method by a normal full echo method in an MRI apparatus. That is, when the number of data in the frequency axis direction of the k-space is N, the number of samplings N / time preceding the echo center (the time position after the echo time TE from the excitation pulse in the pulse sequence SG for proton excitation) 2 and the echo signal is sampled at N sampling points of the sampling number N / 2, which is later than the echo center in time, and sampling data is collected. Therefore, as shown in FIG. 9, the collected sampling data fills the k-space S. In image reconstruction, two-dimensional Fourier transform is performed on the data in the k-space S to obtain a real number component Re and an imaginary number component Im, and an image having an absolute value Z = √ {Re 2 + Im 2 } is generated.

【0003】図10は、MRI装置における従来の非対
称エコー法によるデータ収集方法の説明図である。すな
わち、エコー中心より時間的に前のサンプリング数m
(通常、mはN×10%程度である。)およびエコー中
心より時間的に後のサンプリング数N/2の(m+N/
2)個のサンプリング点でエコー信号をサンプリング
し、サンプリングデータを収集する。従って、図11に
示すように、収集したサンプリングデータはk−スペー
スを埋めつくさず、欠けている部分Eを生じる。画像再
構成では、まず、周波数軸方向の中心軸Gyがエコー中
心に対応するk−スペースSを想定し、前記中心軸Gy
に対称な(m+m)の範囲のサンプリングデータを用い
て位相補正を行う。この位相補正を行うのは、主磁場不
均一,生体内の磁化率の違いによる局所的磁場不均一,
血流による位相シフトなどに起因する位相歪をサンプリ
ングデータが含むので、この位相歪を除去するためであ
る。そして、k−スペースの半分のデータH(図11)
に対して2次元フーリエ変換を行い、実数成分Reと虚
数成分Imを得て、その実数成分Reの画像を生成す
る。通常、虚数成分Imを使用すると望ましい画像が得
られないため、虚数成分Imは使用しない。
FIG. 10 is an explanatory view of a conventional data collection method by the asymmetric echo method in the MRI apparatus. That is, the number of samplings m which is before the echo center in time
(Normally, m is about N × 10%.) And (m + N /) of the sampling number N / 2 after the echo center in time.
2) Echo signals are sampled at the sampling points to collect sampling data. Therefore, as shown in FIG. 11, the collected sampling data does not fill the k-space, resulting in a missing portion E. In the image reconstruction, first, the central axis Gy in the frequency axis direction is assumed to be the k-space S corresponding to the echo center, and the central axis Gy is used.
Phase correction is performed using sampling data in the range (m + m) that is symmetrical with respect to. This phase correction is performed due to inhomogeneity of the main magnetic field, local inhomogeneity due to the difference in magnetic susceptibility in the living body,
This is because the sampling data includes the phase distortion caused by the phase shift due to the blood flow, and the phase distortion is removed. Then, the data H of half the k-space (FIG. 11)
Then, a two-dimensional Fourier transform is performed to obtain a real number component Re and an imaginary number component Im, and an image of the real number component Re is generated. Usually, since the desired image cannot be obtained by using the imaginary number component Im, the imaginary number component Im is not used.

【0004】この特殊な画像再構成は数学的には後述の
(数1)でt0=0とおいた場合に相当し、第2項がゼ
ロとなり、実数成分より強度が1/2の所望の画像が得
られる。
This special image reconstruction mathematically corresponds to the case where t0 = 0 in (Expression 1) described later, the second term becomes zero, and the desired image whose intensity is 1/2 that of the real number component is obtained. Is obtained.

【0005】[0005]

【発明が解決しようとする課題】上記フルエコー法で
は、位相補正は全く要らないが、エコー中心より前にサ
ンプリング数N/2のデータを収集するため、エコー時
間TEを短くすることが出来ない。このため、エコー時
間TEを短くしたい場合は、上記非対称エコー法が用い
られている。しかし、上記従来の非対称エコー法では、
位相補正が必要であるため、画像再構成時間が長くなっ
てしまう問題点がある。なお、位相補正を行わないと、
位相歪によって、画像にアーチファクトを生じてしま
う。また、従来例に示すような特殊な画像再構成を用い
ないとリンギングアーチファクトを生じてしまう。そこ
で、この発明の目的は、非対称エコー法を改良して、上
記位相補正を不要とし、画像再構成時間を短縮できるよ
うにしたMRI装置におけるデータ収集方法およびイメ
ージング方法およびMRI装置を提供することにある。
In the above-mentioned full echo method, no phase correction is required, but since the sampling number N / 2 of data is collected before the echo center, the echo time TE cannot be shortened. Therefore, when it is desired to shorten the echo time TE, the asymmetric echo method is used. However, in the above conventional asymmetric echo method,
Since the phase correction is necessary, there is a problem that the image reconstruction time becomes long. If you do not correct the phase,
The phase distortion causes artifacts in the image. In addition, ringing artifacts will occur unless a special image reconstruction as in the conventional example is used. Therefore, an object of the present invention is to provide a data acquisition method, an imaging method, and an MRI apparatus in an MRI apparatus, which is an improvement of the asymmetric echo method so that the above-mentioned phase correction is unnecessary and the image reconstruction time can be shortened. is there.

【0006】[0006]

【課題を解決するための手段】第1の観点では、この発
明は、MRI装置において、エコー中心より時間的に前
のサンプリング数をm、エコー中心より時間的に後のサ
ンプリング数をp、k−スペースの周波数軸方向の全デ
ータ数をNとするとき、m<p、且つ、N×20%≦m
≦N×35%、且つ、m+p≦Nとしてエコー信号をサ
ンプリングしデータを収集することを特徴とするMRI
装置におけるデータ収集方法を提供する。
According to a first aspect of the present invention, in an MRI apparatus, the sampling number before the echo center is m, and the sampling number after the echo center is p, k. -Where N is the total number of data in the frequency axis direction of the space, m <p and N × 20% ≦ m
MRI characterized in that the echo signal is sampled and data is collected with ≦ N × 35% and m + p ≦ N.
A data collection method in an apparatus is provided.

【0007】第2の観点では、この発明は、MRI装置
において、エコー中心より時間的に前のサンプリング数
をm、エコー中心より時間的に後のサンプリング数を
p、k−スペースの周波数軸方向の全データ数をNとす
るとき、m<p、且つ、N×20%≦m≦N×35%、
且つ、m+p=Nとしてエコー信号をサンプリングしサ
ンプリングデータを収集し、そのサンプリングデータに
対して2次元フーリエ変換を施し、実数成分Reおよび
虚数成分Imを得て、絶対値Z=√{Re2+Im2}の
画像を生成することを特徴とするMRI装置におけるイ
メージング方法を提供する。
According to a second aspect of the present invention, in the MRI apparatus, the number of samplings before the echo center in time is m, the number of samplings after the echo center in time is p, and the number of samplings after the echo center is in the frequency axis direction of the k-space. , Where N is the total number of data in N, and N × 20% ≦ m ≦ N × 35%,
In addition, the echo signal is sampled with m + p = N, sampling data is collected, two-dimensional Fourier transform is applied to the sampling data, the real number component Re and the imaginary number component Im are obtained, and the absolute value Z = √ {Re 2 + Im 2 } is provided to provide an imaging method in an MRI apparatus.

【0008】第3の観点では、この発明は、MRI装置
において、エコー中心より時間的に前のサンプリング数
をm、エコー中心より時間的に後のサンプリング数を
p、k−スペースの周波数軸方向の全データ数をNとす
るとき、m<p、且つ、N×20%≦m≦N×35%、
且つ、m+p<Nとしてエコー信号をサンプリングしサ
ンプリングデータを収集し、周波数軸方向の中心がエコ
ー中心に対応するk−スペース中で前記サンプリングデ
ータの欠けた部分に0を代入し、そのk−スペースのデ
ータに対して2次元フーリエ変換を施し、実数成分Re
および虚数成分Imを得て、絶対値Z=√{Re2+I
2}の画像を生成することを特徴とするMRI装置に
おけるイメージング方法を提供する。
According to a third aspect of the present invention, in the MRI apparatus, the number of samplings before the echo center in time is m, the number of samplings after the echo center in time is p, and the number of samplings after the echo center is in the frequency axis direction of the k-space. , Where N is the total number of data in N, and N × 20% ≦ m ≦ N × 35%,
Also, the echo signal is sampled by collecting m + p <N, sampling data is collected, 0 is substituted in the missing portion of the sampling data in the k-space whose center in the frequency axis direction corresponds to the echo center, and the k-space 2D Fourier transform is applied to the data of
And the imaginary number component Im, the absolute value Z = √ {Re 2 + I
There is provided an imaging method in an MRI apparatus characterized by generating an image of m 2 }.

【0009】第4の観点では、この発明は、エコー中心
より時間的に前のサンプリング数mがエコー中心より時
間的に後のサンプリング数pより少ない非対称のサンプ
リング点でデータを収集するMRI装置において、k−
スペースの周波数軸方向の全データ数をNとするとき、
N×20%≦m≦N×35%、且つ、m+p≦Nとして
エコー信号をサンプリングしデータを収集するデータ収
集手段を具備したことを特徴とするMRI装置を提供す
る。
In a fourth aspect, the present invention is an MRI apparatus for collecting data at an asymmetrical sampling point in which the sampling number m before the echo center is smaller than the sampling number p after the echo center in time. , K-
When the total number of data in the space frequency axis direction is N,
There is provided an MRI apparatus comprising a data collection means for sampling an echo signal and collecting data by setting N × 20% ≦ m ≦ N × 35% and m + p ≦ N.

【0010】第5の観点では、この発明は、エコー中心
より時間的に前のサンプリング数mがエコー中心より時
間的に後のサンプリング数pより少ない非対称のサンプ
リング点でデータを取得するMRI装置において、k−
スペースの周波数軸方向の全データ数をNとするとき、
N×20%≦m≦N×35%、且つ、m+p=Nとして
エコー信号をサンプリングしサンプリングをデータを収
集するデータ収集手段と、そのサンプリングデータに対
して2次元フーリエ変換を施し実数成分Reおよび虚数
成分Imを得て絶対値Z=√{Re2+Im2}の画像を
生成する画像再構成手段とを具備したことを特徴とする
MRI装置を提供する。
In a fifth aspect, the present invention is an MRI apparatus for acquiring data at an asymmetric sampling point in which the sampling number m before the echo center in time is smaller than the sampling number p after the echo center in time. , K-
When the total number of data in the space frequency axis direction is N,
N × 20% ≦ m ≦ N × 35%, and m + p = N, data collection means for sampling an echo signal and collecting data by sampling, and two-dimensional Fourier transform for the sampling data, and a real number component Re and An MRI apparatus comprising: an image reconstructing unit that obtains an imaginary number component Im and generates an image having an absolute value Z = √ {Re 2 + Im 2 }.

【0011】第6の観点では、この発明は、エコー中心
より時間的に前のサンプリング数mがエコー中心より時
間的に後のサンプリング数pより少ない非対称のサンプ
リング点でデータを取得するMRI装置において、k−
スペースの周波数軸方向の全データ数をNとするとき、
N×20%≦m≦N×35%、且つ、m+p<Nとして
エコー信号をサンプリングしサンプリングをデータを収
集するデータ収集手段と、周波数軸方向の中心がエコー
中心に対応するk−スペース中で前記サンプリングデー
タの欠けた部分に“0”を代入しそのk−スペースのデ
ータに対して2次元フーリエ変換を施し実数成分Reお
よび虚数成分Imを得て絶対値Z=√{Re2+Im2
の画像を生成する画像再構成手段とを具備したことを特
徴とするMRI装置を提供する。
In a sixth aspect, the present invention is an MRI apparatus for acquiring data at an asymmetric sampling point in which the sampling number m before the echo center is smaller than the sampling number p after the echo center in time. , K-
When the total number of data in the space frequency axis direction is N,
In N × 20% ≦ m ≦ N × 35%, and m + p <N, data collecting means for sampling the echo signal and collecting the sampling data, and in the k-space whose center in the frequency axis direction corresponds to the echo center "0" is substituted for the missing portion of the sampling data, the two-dimensional Fourier transform is applied to the data in the k-space to obtain the real number component Re and the imaginary number component Im, and the absolute value Z = √ {Re 2 + Im 2 }
And an image reconstructing means for generating an image of the above.

【0012】[0012]

【作用】図12に示す非対称エコー信号x’(t)は、
非対称エコー法によるサンプリングデータに相当する。
エコー中心を時刻t=0とし、また、非対称エコー法で
のサンプリング開始時刻を−t0とし、 h(t)=1 t≧−t0 h(t)=0 t<−t0 なる窓関数h(t)を想定すると、前記非対称エコー信
号x’(t)は、窓関数h(t)を対称エコー信号x
(t)に乗算したものと等価であるから、 x’(t)=x(t)・h(t) と表せる。そこで、非対称エコー信号x’(t)のフー
リエ変換F[x’(t)]は、 F[x’(t)]=F[x(t)・h(t)] =F[x(t)]*F[h(t)] …(a) となる。ここで、t0=0とすると、 F[h(t)]=δ(f)/2−i/(2πf) (但し、δ(f)はデルタ関数) であるから、 F[h(t+t0)] =exp{i・2π・t0}・(δ(f)/2−i/(2πf)) …(b) となる。式(a)に式(b)を代入して整理し、 F[x’(t)]=I(f) F[x(t)] =X(f) とおけば、次の(数1)が得られる。
The asymmetric echo signal x '(t) shown in FIG.
It corresponds to sampling data by the asymmetric echo method.
The center of echo is set to time t = 0, the sampling start time in the asymmetric echo method is set to -t0, and h (t) = 1 t ≧ -t0 h (t) = 0 t <-t0 ) Is assumed, the asymmetric echo signal x ′ (t) has a window function h (t) as the symmetric echo signal x
Since it is equivalent to the product of (t) multiplied, it can be expressed as x ′ (t) = x (t) · h (t). Therefore, the Fourier transform F [x ′ (t)] of the asymmetric echo signal x ′ (t) is F [x ′ (t)] = F [x (t) · h (t)] = F [x (t) )] * F [h (t)] (a). Here, if t0 = 0, then F [h (t)] = δ (f) / 2−i / (2πf) (where δ (f) is a delta function), so F [h (t + t0) ] = Exp {i · 2π · t0} · (δ (f) / 2−i / (2πf)) (b) By substituting the expression (b) into the expression (a) and rearranging, F [x ′ (t)] = I (f) F [x (t)] = X (f) ) Is obtained.

【0013】[0013]

【数1】 [Equation 1]

【0014】(数1)の実数成分は、対称エコー信号x
(t)の1/2パワーの画像と,t0により決まる周波
数成分のシンク(sinc)関数と対称エコー信号のた
たみ込みの和になっている。このことは、t0が小さい
場合はリンギングアーチファクトを無視できないが、t
0が大きい場合はリンギングアーチファクトを無視でき
ることを示している。従来のようにt0が小さいデータ
で従来例のような特殊な画像再構成を用いないとリンギ
ングアーチファクトが問題となる。
The real number component of (Equation 1) is a symmetric echo signal x
It is the sum of the half power image of (t), the sinc function of the frequency component determined by t0, and the convolution of the symmetrical echo signal. This means that ringing artifacts cannot be ignored if t0 is small, but
A large 0 indicates that ringing artifacts can be ignored. If t0 is small as in the conventional case and the special image reconstruction as in the conventional example is not used, the ringing artifact becomes a problem.

【0015】さて、この発明のMRI装置におけるデー
タ収集方法およびイメージング方法およびMRI装置で
は、エコー中心より時間的に前のサンプリング数mをk
−スペースの周波数軸方向の全データ数Nの20%以上
とするから、従来の非対称エコー法(通常、10%程
度)に比べて、上記t0が大きくなっている。従って、
位相補正を行わないが、位相歪によるアーチファクトを
低減できることになる。また、従来例のような特殊な画
像再構成を行わなくてもリンギングアーチファクトが無
視できるぐらい小さい。また、エコー中心より時間的に
前のサンプリング数mをk−スペースの周波数軸方向の
全データ数Nの35%以下とするから、非対称エコー法
の利点を失わず、短いエコー時間TEに対応できる。そ
して、位相補正を行わないため、画像再構成時間を短縮
できるようになる。
In the data collecting method, the imaging method and the MRI apparatus in the MRI apparatus of the present invention, the sampling number m before the echo center in time is k.
Since the space is 20% or more of the total number N of data in the frequency axis direction, the above-mentioned t0 is larger than that of the conventional asymmetric echo method (usually about 10%). Therefore,
Although no phase correction is performed, the artifacts due to the phase distortion can be reduced. Also, ringing artifacts are so small that they can be ignored without performing special image reconstruction as in the conventional example. Moreover, since the sampling number m temporally before the echo center is set to 35% or less of the total number N of data in the frequency axis direction of the k-space, the advantage of the asymmetric echo method is not lost and a short echo time TE can be supported. . Since the phase correction is not performed, the image reconstruction time can be shortened.

【0016】[0016]

【実施例】以下、図に示す実施例によりこの発明をさら
に詳しく説明する。なお、これによりこの発明が限定さ
れるものではない。図1は、この発明のMRI装置の一
実施例のブロック図である。このMRI装置100にお
いて、マグネットアセンブリ1は、内部に被検体を挿入
するための空間部分(孔)を有し、この空間部分を取り
まくようにして、被検体に一定の静磁場を印加する静磁
場コイルと、勾配磁場を発生するための勾配磁場コイル
(勾配コイルは、読み出し,位相エンコード,スライス
選択の各軸のコイルを備えている)と、被検体内の原子
核のスピンを励起するためのRFパルスを与える送信コ
イルと、被検体からのNMR(Nuclear Magnetic Reson
ance)信号を検出する受信コイル等が配置されている。
静磁場コイル,勾配磁場コイル,送信コイルおよび受信
コイルは、それぞれ主磁場電源2,勾配磁場駆動回路
3,RF電力増幅器4および前置増幅器5に接続されて
いる。
The present invention will be described in more detail with reference to the embodiments shown in the drawings. The present invention is not limited to this. FIG. 1 is a block diagram of an embodiment of the MRI apparatus of the present invention. In this MRI apparatus 100, the magnet assembly 1 has a space portion (hole) for inserting a subject therein, and a static magnetic field for applying a constant static magnetic field to the subject so as to surround this space portion. A coil, a gradient magnetic field coil for generating a gradient magnetic field (the gradient coil includes coils for reading, phase encoding, and slice selection axes), and an RF for exciting spins of nuclei in a subject. A transmitter coil that gives a pulse and an NMR (Nuclear Magnetic Resonance) from the subject.
A receiving coil or the like for detecting the signal is arranged.
The static magnetic field coil, the gradient magnetic field coil, the transmitting coil and the receiving coil are connected to the main magnetic field power source 2, the gradient magnetic field driving circuit 3, the RF power amplifier 4 and the preamplifier 5, respectively.

【0017】シーケンス記憶回路8は、計算機7からの
指令に従い、スピンエコー法やグラジエントエコー法な
どのパルス系列に基づいて、勾配磁場駆動回路3を操作
し、前記マグネットアセンブリ1の勾配磁場コイルから
勾配磁場を発生させると共に、ゲート変調回路9を操作
し、RF発振回路10の高周波出力信号を所定タイミン
グ・所定包絡線のパルス状信号に変調し、それをRFパ
ルスとしてRF電力増幅器4に加え、RF電力増幅器4
でパワー増幅した後、前記マグネットアセンブリ1の送
信コイルに印加し、被検体の目的のスライスを選択励起
する。
The sequence storage circuit 8 operates the gradient magnetic field drive circuit 3 in accordance with a command from the computer 7 on the basis of a pulse sequence such as a spin echo method or a gradient echo method, so that the gradient magnetic field coil of the magnet assembly 1 receives a gradient. While generating a magnetic field, the gate modulation circuit 9 is operated to modulate the high frequency output signal of the RF oscillation circuit 10 into a pulsed signal having a predetermined timing and a predetermined envelope, which is applied as an RF pulse to the RF power amplifier 4 to generate an RF signal. Power amplifier 4
After the power is amplified by (1), it is applied to the transmission coil of the magnet assembly 1 to selectively excite the target slice of the subject.

【0018】前置増幅器5は、マグネットアセンブリ1
の受信コイルで検出された被検体からのNMR信号を増
幅し、位相検波器12に入力する。位相検波器12は、
RF発振回路10の出力を参照信号とし、前置増幅器5
からのNMR信号を位相検波して、A/D変換器11に
与える。A/D変換器11は、位相検波後のアナログ信
号をディジタル信号に変換して、計算機7に入力する。
計算機7は、A/D変換器11からのデジタル信号に対
する画像再構成演算を行い、目的のスライスのイメージ
(プロトン密度像)を生成する。このイメージは、表示
装置6にて表示される。また、計算機7は、操作卓13
から入力された情報を受け取るなどの全体的な制御を受
け持つ。
The preamplifier 5 includes a magnet assembly 1
The NMR signal from the subject detected by the receiving coil is amplified and input to the phase detector 12. The phase detector 12 is
The output of the RF oscillation circuit 10 is used as a reference signal, and the preamplifier 5
The NMR signal from is phase-detected and given to the A / D converter 11. The A / D converter 11 converts the analog signal after phase detection into a digital signal and inputs it to the computer 7.
The computer 7 performs an image reconstruction operation on the digital signal from the A / D converter 11 to generate an image of a target slice (proton density image). This image is displayed on the display device 6. In addition, the computer 7 is a console 13
Responsible for overall control such as receiving information input from.

【0019】図2は、上記MRI装置100における非
対称エコー法によるデータ収集方法の第1例の説明図で
ある。すなわち、エコー中心より時間的に前のサンプリ
ング数m(第1例では、mはN×25%とする。)およ
びエコー中心より時間的に後のサンプリング数p(第1
例では、pはN×75%とする。)のN個のサンプリン
グ点でエコー信号をサンプリングし、サンプリングデー
タを収集する。従って、図3に示すように、収集したサ
ンプリングデータはk−スペースを埋めつくす。画像再
構成では、前記k−スペースSのデータに対し2次元フ
ーリエ変換を行い、実数成分Reと虚数成分Imを得
て、絶対値Z=√{Re2+Im2}の画像を生成する。
以上の第1例のデータ収集方法およびイメージング方法
では、エコー中心より時間的に前のサンプリング数mを
k−スペースの周波数軸方向の全データ数Nの25%と
するから、t0が大きく、位相補正を行わないが位相歪
によるアーチファクトを低減できる。また、従来例のよ
うな特殊な画像再構成を行わなくてもリンギングアーチ
ファクトが無視できるぐらい小さい。また、非対称エコ
ー法の利点を失わず、短いエコー時間TEに対応でき
る。また、位相補正を行わないため、画像再構成時間を
短縮することが出来る。さらに、従来は収集していなか
った高周波成分(図3に示すQの部分)を含むため、画
像分解能を上げることが出来る。
FIG. 2 is an explanatory diagram of a first example of a data collection method by the asymmetric echo method in the MRI apparatus 100. That is, the sampling number m temporally before the echo center (m is N × 25% in the first example) and the sampling number p temporally after the echo center (first number)
In the example, p is N × 75%. ), The echo signal is sampled at N sampling points to collect sampling data. Therefore, as shown in FIG. 3, the collected sampling data fills the k-space. In image reconstruction, two-dimensional Fourier transform is performed on the data in the k-space S to obtain a real number component Re and an imaginary number component Im, and an image having an absolute value Z = √ {Re 2 + Im 2 } is generated.
In the above-described data acquisition method and imaging method of the first example, the sampling number m temporally before the echo center is set to 25% of the total number N of data in the frequency axis direction of the k-space, so that t0 is large and the phase is large. Although not corrected, the artifacts due to phase distortion can be reduced. Also, ringing artifacts are so small that they can be ignored without performing special image reconstruction as in the conventional example. Moreover, the advantage of the asymmetrical echo method is not lost, and a short echo time TE can be dealt with. Further, since the phase correction is not performed, the image reconstruction time can be shortened. Further, since the high frequency component (Q portion shown in FIG. 3) which has not been collected conventionally is included, the image resolution can be improved.

【0020】図4は、上記MRI装置100における非
対称エコー法によるデータ収集方法の第2例の説明図で
ある。すなわち、エコー中心より時間的に前のサンプリ
ング数m(第2例では、mはN×35%とする。)およ
びエコー中心より時間的に後のサンプリング数p(第2
例では、pはN×50%とする。)の(m+p)個のサ
ンプリング点でエコー信号をサンプリングし、サンプリ
ングデータを収集する。従って、図5に示すように、収
集したサンプリングデータはk−スペースSを埋めつく
さず、欠けている部分F1を生じる。画像再構成では、
まず、周波数軸方向の中心軸Gyがエコー中心に対応す
るk−スペースSを想定し、サンプリングデータの欠け
ている部分F1に“0”を代入する。そして、そのk−
スペースSのデータに対して2次元フーリエ変換を施
し、実数成分Reおよび虚数成分Imを得て、絶対値Z
=√{Re2+Im2}の画像を生成する。
FIG. 4 is an explanatory diagram of a second example of the data collection method by the asymmetric echo method in the MRI apparatus 100. That is, the sampling number m temporally before the echo center (m is N × 35% in the second example) and the sampling number p temporally after the echo center (second)
In the example, p is N × 50%. ), The echo signal is sampled at (m + p) sampling points to collect sampling data. Therefore, as shown in FIG. 5, the collected sampling data does not fill the k-space S, resulting in a missing portion F1. In image reconstruction,
First, assuming that the center axis Gy in the frequency axis direction is the k-space S corresponding to the center of the echo, "0" is substituted for the missing portion F1 of the sampling data. And that k-
Two-dimensional Fourier transform is applied to the data of the space S to obtain the real number component Re and the imaginary number component Im, and the absolute value Z
An image of = √ {Re 2 + Im 2 } is generated.

【0021】以上の第2例のデータ収集方法およびイメ
ージング方法では、エコー中心より時間的に前のサンプ
リング数mをk−スペースの周波数軸方向の全データ数
Nの35%とするから、t0が大きく、位相補正を行わ
ないがリンギングアーチファクトを低減できる。また、
非対称エコー法の利点を失わず、短いエコー時間TEに
対応できる。また、位相補正を行わないため、画像再構
成時間を短縮することが出来る。
In the above-described data acquisition method and imaging method of the second example, since the sampling number m temporally before the echo center is set to 35% of the total number N of data in the frequency axis direction of k-space, t0 is Large, no phase correction is performed, but ringing artifacts can be reduced. Also,
A short echo time TE can be supported without losing the advantage of the asymmetric echo method. Further, since the phase correction is not performed, the image reconstruction time can be shortened.

【0022】図6は、上記MRI装置100における非
対称エコー法によるデータ収集方法の第3例の説明図で
ある。すなわち、エコー中心より時間的に前のサンプリ
ング数m(第3例では、mはN×20%とする。)およ
びエコー中心より時間的に後のサンプリング数p(第3
例では、pはN×30%とする。)の(m+p)個のサ
ンプリング点でエコー信号をサンプリングし、サンプリ
ングデータを収集する。従って、図7に示すように、収
集したサンプリングデータはk−スペースSを埋めつく
さず、欠けている部分F1,F2を生じる。画像再構成
では、まず、周波数軸方向の中心軸Gyがエコー中心に
対応するk−スペースSを想定し、サンプリングデータ
の欠けている部分F1,F2に“0”を代入し、そのk
−スペースSのデータに対して2次元フーリエ変換を施
し、実数成分Reおよび虚数成分Imを得て、絶対値Z
=√{Re2+Im2}の画像を生成する。
FIG. 6 is an explanatory diagram of a third example of the data collection method by the asymmetric echo method in the MRI apparatus 100. That is, the sampling number m temporally before the echo center (m is N × 20% in the third example) and the sampling number p temporally after the echo center (third example).
In the example, p is N × 30%. ), The echo signal is sampled at (m + p) sampling points to collect sampling data. Therefore, as shown in FIG. 7, the collected sampling data does not completely fill the k-space S and produces missing portions F1 and F2. In image reconstruction, first, assuming that the central axis Gy in the frequency axis direction corresponds to the echo center, a k-space S is assumed, and "0" is substituted for the missing portions F1 and F2 of the sampling data, and k
Two-dimensional Fourier transform is performed on the data in the space S to obtain the real number component Re and the imaginary number component Im, and the absolute value Z
An image of = √ {Re 2 + Im 2 } is generated.

【0023】以上の第3例のデータ収集方法およびイメ
ージング方法では、エコー中心より時間的に前のサンプ
リング数mをk−スペースの周波数軸方向の全データ数
Nの20%とするから、t0が大きく、位相補正を行わ
ないが位相歪によるアーチファクトおよびリンギングア
ーチファクトを低減できる。また、非対称エコー法の利
点を失わず、短いエコー時間TEに対応できる。また、
位相補正を行わないため、画像再構成時間を短縮するこ
とが出来る。
In the above-described data acquisition method and imaging method of the third example, since the sampling number m temporally before the echo center is set to 20% of the total number N of data in the frequency axis direction of k-space, t0 is Although it is large and phase correction is not performed, it is possible to reduce artifacts and ringing artifacts due to phase distortion. Moreover, the advantage of the asymmetrical echo method is not lost, and a short echo time TE can be dealt with. Also,
Since no phase correction is performed, the image reconstruction time can be shortened.

【0024】[0024]

【発明の効果】この発明のMRI装置におけるデータ収
集方法およびイメージング方法およびMRI装置によれ
ば、エコー中心より時間的に前のサンプリング数mをk
−スペースの周波数軸方向の全データ数Nの20%以上
35%以下とするから、位相補正を行わなくても位相歪
によるアーチファクトおよびリンギングアーチファクト
を低減できる。また、非対称エコー法の利点を失わず、
短いエコー時間TEに対応できる。そして、位相補正を
行わないため、画像再構成時間を短縮することが出来
る。特に、タイムオブフライト・アンジオグラフィにお
いて有用である。
According to the data collecting method, the imaging method, and the MRI apparatus in the MRI apparatus of the present invention, the sampling number m before the echo center in time is k.
-Because 20% or more and 35% or less of the total number N of data in the frequency axis direction of the space, it is possible to reduce the artifacts and ringing artifacts due to the phase distortion without performing the phase correction. Also, without losing the advantages of the asymmetric echo method,
It can cope with a short echo time TE. Since the phase correction is not performed, the image reconstruction time can be shortened. It is especially useful in time-of-flight angiography.

【図面の簡単な説明】[Brief description of drawings]

【図1】この発明のMRI装置の一実施例のブロック図
である
FIG. 1 is a block diagram of an embodiment of an MRI apparatus of the present invention.

【図2】図1のMRI装置における非対称エコー法によ
るデータ収集方法の第1例の説明図である。
FIG. 2 is an explanatory diagram of a first example of a data collection method by the asymmetric echo method in the MRI apparatus of FIG.

【図3】図2のデータ収集方法にかかるk−スペースの
説明図である。
FIG. 3 is an explanatory diagram of k-space according to the data collection method of FIG.

【図4】図1のMRI装置における非対称エコー法によ
るデータ収集方法の第2例の説明図である。
4 is an explanatory diagram of a second example of a data collection method by the asymmetric echo method in the MRI apparatus of FIG.

【図5】図4のデータ収集方法にかかるk−スペースの
説明図である。
5 is an explanatory diagram of k-space according to the data collection method of FIG.

【図6】図1のMRI装置における非対称エコー法によ
るデータ収集方法の第3例の説明図である。
6 is an explanatory diagram of a third example of a data collection method by the asymmetric echo method in the MRI apparatus of FIG.

【図7】図6のデータ収集方法にかかるk−スペースの
説明図である。
7 is an explanatory diagram of k-space according to the data collection method of FIG.

【図8】フルエコー法によるデータ収集方法の説明図で
ある。
FIG. 8 is an explanatory diagram of a data collection method by the full echo method.

【図9】図8のデータ収集方法にかかるk−スペースの
説明図である。
9 is an explanatory diagram of k-space according to the data collection method of FIG.

【図10】従来の非対称エコー方法によるデータ収集方
法の説明図である。
FIG. 10 is an explanatory diagram of a data collection method by a conventional asymmetric echo method.

【図11】図10のデータ収集方法にかかるk−スペー
スの説明図である。
11 is an explanatory diagram of k-space according to the data collection method of FIG.

【図12】非対称エコー信号と対称エコー信号と窓関数
の説明図である。
FIG. 12 is an explanatory diagram of an asymmetric echo signal, a symmetric echo signal, and a window function.

【符号の説明】[Explanation of symbols]

100 MR装置 1 マグネットアセンブリ 7 計算機 8 シーケンスコントローラ 9 ゲート変調回路 10 RF発振回路 13 操作卓 N k−スペースの周波数軸方向の全デー
タ数 m エコー中心より時間的に前のサンプリ
ング数 p エコー中心より時間的に後のサンプリ
ング数 TE エコー時間 SG プロトン励起用のパルス系列
100 MR device 1 Magnet assembly 7 Computer 8 Sequence controller 9 Gate modulation circuit 10 RF oscillation circuit 13 Operator's console N k-Number of all data in the frequency axis direction of k-space m Number of samplings temporally before echo center p Time from echo center Sampling number after TE TE Echo time SG Pulse sequence for proton excitation

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】 MRI装置において、エコー中心より時
間的に前のサンプリング数をm、エコー中心より時間的
に後のサンプリング数をp、k−スペースの周波数軸方
向の全データ数をNとするとき、 m<p、且つ、N×20%≦m≦N×35%、且つ、m
+p≦Nとしてエコー信号をサンプリングしデータを収
集することを特徴とするMRI装置におけるデータ収集
方法。
1. In the MRI apparatus, the number of samplings before the echo center in time is m, the number of samplings after the echo center in time is p, and the total number of data in the frequency axis direction of the k-space is N. When, m <p, N × 20% ≦ m ≦ N × 35%, and m
A data collection method in an MRI apparatus, characterized in that an echo signal is sampled as + p ≦ N to collect data.
【請求項2】 MRI装置において、エコー中心より時
間的に前のサンプリング数をm、エコー中心より時間的
に後のサンプリング数をp、k−スペースの周波数軸方
向の全データ数をNとするとき、 m<p、且つ、N×20%≦m≦N×35%、且つ、m
+p=Nとしてエコー信号をサンプリングしサンプリン
グデータを収集し、そのサンプリングデータに対して2
次元フーリエ変換を施し、実数成分Reおよび虚数成分
Imを得て、絶対値Z=√{Re2+Im2}の画像を生
成することを特徴とするMRI装置におけるイメージン
グ方法。
2. In the MRI apparatus, the number of samples before the echo center in time is m, the number of samples after the echo center in time is p, and the total number of data in the frequency axis direction of k-space is N. When, m <p, N × 20% ≦ m ≦ N × 35%, and m
+ P = N, the echo signal is sampled, sampling data is collected, and 2 is applied to the sampling data.
A method for imaging in an MRI apparatus, characterized in that a two -dimensional Fourier transform is performed to obtain a real number component Re and an imaginary number component Im to generate an image having an absolute value Z = √ {Re 2 + Im 2 }.
【請求項3】 MRI装置において、エコー中心より時
間的に前のサンプリング数をm、エコー中心より時間的
に後のサンプリング数をp、k−スペースの周波数軸方
向の全データ数をNとするとき、 m<p、且つ、N×20%≦m≦N×35%、且つ、m
+p<Nとしてエコー信号をサンプリングしサンプリン
グデータを収集し、周波数軸方向の中心がエコー中心に
対応するk−スペース中で前記サンプリングデータの欠
けた部分に“0”を代入し、そのk−スペースのデータ
に対して2次元フーリエ変換を施し、実数成分Reおよ
び虚数成分Imを得て、絶対値Z=√{Re2+Im2
の画像を生成することを特徴とするMRI装置における
イメージング方法。
3. In the MRI apparatus, the number of samplings before the echo center in time is m, the number of samplings after the echo center in time is p, and the total number of data in the frequency axis direction of k-space is N. When, m <p, N × 20% ≦ m ≦ N × 35%, and m
The echo signal is sampled by collecting + p <N, sampling data is collected, "0" is substituted for the missing portion of the sampling data in the k-space whose center in the frequency axis direction corresponds to the echo center, and the k-space 2D Fourier transform is performed on the data of to obtain the real number component Re and the imaginary number component Im, and the absolute value Z = √ {Re 2 + Im 2 }
Imaging method in an MRI apparatus, which is characterized by generating an image.
【請求項4】 エコー中心より時間的に前のサンプリン
グ数mがエコー中心より時間的に後のサンプリング数p
より少ない非対称のサンプリング点でデータを収集する
MRI装置において、 k−スペースの周波数軸方向の全データ数をNとすると
き、 N×20%≦m≦N×35%、且つ、m+p≦Nとして
エコー信号をサンプリングしデータを収集するデータ収
集手段を具備したことを特徴とするMRI装置。
4. A sampling number m preceding the echo center in time is a sampling number p following the echo center in time.
In an MRI apparatus that collects data at less asymmetrical sampling points, when N is the total number of data in the frequency axis direction of k-space, N × 20% ≦ m ≦ N × 35% and m + p ≦ N An MRI apparatus comprising a data collecting means for sampling an echo signal and collecting data.
【請求項5】 エコー中心より時間的に前のサンプリン
グ数mがエコー中心より時間的に後のサンプリング数p
より少ない非対称のサンプリング点でデータを取得する
MRI装置において、 k−スペースの周波数軸方向の全データ数をNとすると
き、 N×20%≦m≦N×35%、且つ、m+p=Nとして
エコー信号をサンプリングしサンプリングをデータを収
集するデータ収集手段と、そのサンプリングデータに対
して2次元フーリエ変換を施し実数成分Reおよび虚数
成分Imを得て絶対値Z=√{Re2+Im2}の画像を
生成する画像再構成手段とを具備したことを特徴とする
MRI装置。
5. A sampling number m preceding the echo center in time is a sampling number p following the echo center in time.
In an MRI apparatus that acquires data at a smaller number of asymmetrical sampling points, when N is the total number of data in the frequency axis direction of k-space, N × 20% ≦ m ≦ N × 35%, and m + p = N A data collecting means for sampling an echo signal and collecting data by sampling, and a two-dimensional Fourier transform on the sampling data to obtain a real number component Re and an imaginary number component Im to obtain an absolute value Z = √ {Re 2 + Im 2 }. An MRI apparatus comprising: an image reconstructing unit that generates an image.
【請求項6】 エコー中心より時間的に前のサンプリン
グ数mがエコー中心より時間的に後のサンプリング数p
より少ない非対称のサンプリング点でデータを取得する
MRI装置において、 k−スペースの周波数軸方向の全データ数をNとすると
き、 N×20%≦m≦N×35%、且つ、m+p<Nとして
エコー信号をサンプリングしサンプリングをデータを収
集するデータ収集手段と、周波数軸方向の中心がエコー
中心に対応するk−スペース中で前記サンプリングデー
タの欠けた部分に“0”を代入しそのk−スペースのデ
ータに対して2次元フーリエ変換を施し実数成分Reお
よび虚数成分Imを得て絶対値Z=√{Re2+Im2
の画像を生成する画像再構成手段とを具備したことを特
徴とするMRI装置。
6. A sampling number m preceding the echo center in time is a sampling number p following the echo center in time.
In an MRI apparatus that acquires data at less asymmetrical sampling points, when N is the total number of data in the frequency axis direction of k-space, N × 20% ≦ m ≦ N × 35% and m + p <N Data collecting means for sampling an echo signal and collecting data by sampling, and "0" is assigned to the missing portion of the sampling data in the k-space whose center in the frequency axis direction corresponds to the echo center, and the k-space 2D Fourier transform is applied to the data of to obtain the real number component Re and the imaginary number component Im, and the absolute value Z = √ {Re 2 + Im 2 }
And an image reconstructing means for generating an image of the MRI apparatus.
JP03374594A 1994-03-03 1994-03-03 MRI equipment Expired - Lifetime JP3447099B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP03374594A JP3447099B2 (en) 1994-03-03 1994-03-03 MRI equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP03374594A JP3447099B2 (en) 1994-03-03 1994-03-03 MRI equipment

Publications (2)

Publication Number Publication Date
JPH07241280A true JPH07241280A (en) 1995-09-19
JP3447099B2 JP3447099B2 (en) 2003-09-16

Family

ID=12394959

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03374594A Expired - Lifetime JP3447099B2 (en) 1994-03-03 1994-03-03 MRI equipment

Country Status (1)

Country Link
JP (1) JP3447099B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003061928A (en) * 2001-08-14 2003-03-04 Ge Medical Systems Global Technology Co Llc Magnetic resonance signal acquisition device and magnetic resonance imaging apparatus
US9797969B2 (en) 2013-09-12 2017-10-24 Siemens Aktiengesellschaft Method and magnetic resonance apparatus to generate a spatially selective excitation

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003061928A (en) * 2001-08-14 2003-03-04 Ge Medical Systems Global Technology Co Llc Magnetic resonance signal acquisition device and magnetic resonance imaging apparatus
US9797969B2 (en) 2013-09-12 2017-10-24 Siemens Aktiengesellschaft Method and magnetic resonance apparatus to generate a spatially selective excitation

Also Published As

Publication number Publication date
JP3447099B2 (en) 2003-09-16

Similar Documents

Publication Publication Date Title
US7015696B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP3796446B2 (en) MRI equipment
JP3796455B2 (en) MRI equipment
US6528996B1 (en) Diffusion-weighted imaging method and apparatus for fast pulse sequence with MPG pulses
JP2001299720A (en) Method for collecting data of correction of magnetic field drift, method for correcting magnetic field drift and mri device
JPH05300895A (en) Selective excitation method of nuclear spin in mri apparatus
US4713616A (en) Nuclear magnetic resonance imaging apparatus
JP2007090001A (en) Mr scan method and mri apparatus
US5685304A (en) Diffusion sensitizing imaging method and MRI apparatus
US20020050816A1 (en) MR imaging method, phase error measuring method, and MRI apparatus
JP3447099B2 (en) MRI equipment
US20220057467A1 (en) Epi mr imaging with distortion correction
JP3576641B2 (en) MRI equipment
JP3748661B2 (en) MRI equipment
JP2006061235A (en) Magnetic resonance imaging device
JP3091203B2 (en) Magnetic resonance imaging equipment
US20240094320A1 (en) Dixon-type water/fat separation mr imaging
JP3205061B2 (en) Magnetic resonance imaging equipment
KR20010060324A (en) Method and apparatus to compensate for image artifacts caused by magnet vibration in an mr imaging system
JP3615614B2 (en) MRI equipment
JPH10248824A (en) Diffusion accentuated imaging method and mri system
JP3439845B2 (en) Method for generating ramped RF pulse for MRI and MRI apparatus
JPH0824235A (en) Method for collecting data in mri appparatus, and mri apparatus
JPS628046A (en) Method for measuring nuclear magnetic resonance signal
JPH06133944A (en) Mr imaging system

Legal Events

Date Code Title Description
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: 20070704

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20080704

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20080704

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20090704

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20090704

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20090704

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20100704

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20100704

Year of fee payment: 7

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

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

Free format text: PAYMENT UNTIL: 20100704

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

Free format text: PAYMENT UNTIL: 20100704

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20110704

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20110704

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20110704

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20120704

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20120704

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20120704

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20120704

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20130704

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20130704

Year of fee payment: 10

EXPY Cancellation because of completion of term