[go: up one dir, main page]

JPH07163537A - NMR imaging method - Google Patents

NMR imaging method

Info

Publication number
JPH07163537A
JPH07163537A JP6208457A JP20845794A JPH07163537A JP H07163537 A JPH07163537 A JP H07163537A JP 6208457 A JP6208457 A JP 6208457A JP 20845794 A JP20845794 A JP 20845794A JP H07163537 A JPH07163537 A JP H07163537A
Authority
JP
Japan
Prior art keywords
signal
magnetic field
nmr
imaging method
image signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP6208457A
Other languages
Japanese (ja)
Inventor
Koichi Sano
耕一 佐野
Tetsuo Yokoyama
哲夫 横山
Hideaki Koizumi
英明 小泉
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.)
Hitachi Ltd
Original Assignee
Hitachi 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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP6208457A priority Critical patent/JPH07163537A/en
Publication of JPH07163537A publication Critical patent/JPH07163537A/en
Pending legal-status Critical Current

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 【目的】 特殊なハードウェアを用いることなく、NM
R検査装置本体のみで呼吸の動きを検出し、また、その
影響を補正可能としたNMRイメージング方法を提供す
ること。 【構成】 静磁場,傾斜磁場,高周波磁場の発生手段
と、検査対象物からのNMR信号を取出す検出手段と、
検出された信号に対し画像再構成を含む各種演算を行う
手段とを有するNMR検査装置におけるイメージング方
法において、画像信号計測時に呼吸の動きに関する情報
を画像信号とともに計測し、計測した動きに関する情報
に基づいて前記画像信号を補正することを特徴とするN
MRイメージング方法。
(57) [Abstract] [Purpose] NM without using special hardware
To provide an NMR imaging method capable of detecting respiratory movements only by the R inspection apparatus main body and correcting the influence thereof. [Structure] Static magnetic field, gradient magnetic field, high-frequency magnetic field generation means, detection means for extracting an NMR signal from an inspection object,
In an imaging method in an NMR inspection apparatus having means for performing various calculations including image reconstruction on a detected signal, information on respiratory movements is measured together with the image signal at the time of measuring the image signal, and based on the information on the measured movement. N to correct the image signal
MR imaging method.

Description

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

【0001】[0001]

【産業上の利用分野】本発明はNMRイメージング方
法、すなわち、NMR現象を利用した体内断層撮影方式
に関し、特にハードウェアを用いることなしに、呼吸性
の動きに基づく影響を補正するようにしたイメージング
方法に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an NMR imaging method, that is, an internal tomography method using an NMR phenomenon, and an imaging method for correcting an effect based on respiratory motion without using hardware. Regarding the method.

【0002】[0002]

【従来の技術】NMR検査装置におけるイメージングに
おいては、撮影に2〜20分程度の時間を要するため、
胸部や腹部を撮影する場合、呼吸の影響により、画像の
ぼけやアーチファクトを生ずるという問題がある。この
ため、通常は、何等かの手段により、呼吸の動きを検出
して同期をとったり、計測データを補正する手法が用い
られている。例えば、従来の装置では、下記の如き検出
方式が用いられていた。 (1)エアバッグ方式:腹部に固定したエアバッグ内の圧
力が、呼吸に応じて変化することを利用して検出する方
式。 (2)バンド方式:硫酸亜鉛溶液が満たされたチューブを
腹部に巻き付け、呼吸に応じてこのチューブが伸展し、
電気抵抗が変化することを利用して検出する方式。 (3)サーミスタ方式:サーミスタを鼻孔近くに取り付
け、呼吸によって変化する温度を検出する方式。 なお、この種の装置と関連するものとしては、「呼吸同
期NMRの開発」(NMR医学,vol.5,No.1(1985))に
記載されたものがある。
2. Description of the Related Art Imaging in an NMR inspection apparatus requires a time of about 2 to 20 minutes, so that
When photographing the chest and abdomen, there is a problem that blurring of the image and artifacts occur due to the influence of breathing. Therefore, usually, a method of detecting the movement of the breath for synchronization by some means or correcting the measurement data is used. For example, in the conventional device, the following detection method was used. (1) Airbag method: A method in which the pressure inside the airbag fixed to the abdomen is detected by utilizing the fact that it changes according to breathing. (2) Band method: A tube filled with zinc sulfate solution is wrapped around the abdomen, and this tube expands according to breathing,
A method of detecting by utilizing the change in electrical resistance. (3) Thermistor method: A method in which a thermistor is attached near the nostril to detect the temperature that changes with breathing. As a device related to this type of device, there is one described in "Development of Respiratory Synchronous NMR" (NMR Medicine, vol. 5, No. 1 (1985)).

【0003】[0003]

【発明が解決しようとする課題】上記従来技術は、いず
れも、専用ハードウェアを新たに付加する必要があると
同時に、撮影時の操作性にも問題があった。すなわち、
呼吸と同期をとって撮影するため、撮影時間が通常の6
〜8倍かかるという問題もあった。本発明は上記事情に
鑑みてなされたもので、その目的とするところは、従来
の技術における上述の如き問題を解消し、特殊なハード
ウェアを用いることなく、NMR検査装置本体のみで呼
吸の動きを検出し、また、その影響を補正可能としたN
MRイメージング方法を提供することにある。
In each of the above-mentioned conventional techniques, it is necessary to newly add dedicated hardware, and at the same time, there is a problem in operability at the time of photographing. That is,
Since the shooting is performed in synchronization with breathing, the shooting time is normally 6
There was also a problem that it took ~ 8 times. The present invention has been made in view of the above circumstances, and an object of the present invention is to solve the above-described problems in the related art and to perform breathing movements only by the NMR inspection apparatus main body without using special hardware. Is detected, and its influence can be corrected.
An object is to provide an MR imaging method.

【0004】[0004]

【課題を解決するための手段】本発明の上記目的は、静
磁場,傾斜磁場,高周波磁場の発生手段と、検査対象物
からのNMR信号を取出す検出手段と、検出された信号
に対し画像再構成を含む各種演算を行う手段とを有する
NMR検査装置におけるイメージング方法において、画
像信号計測時に呼吸の動きに関する情報を画像信号とと
もに計測し、呼吸の動きに関する情報を位相情報として
得て、これに基づいて前記画像信号を補正することを特
徴とするNMRイメージング方法,画像信号を計測した
直後に、静止物体には影響を与えず、動く物体にのみ影
響を与える傾斜磁場を印加して、呼吸の動きに関する情
報を、速度に応じて異なる信号になるように計測し、こ
れに基づいて前記画像信号を補正することを特徴とする
NMRイメージング方法、または、画像信号を計測する
際に、呼吸の動きに関する情報を、呼吸で動く方向が位
相エンコード方向と一致するようにして、同時期に計測
し、該計測データに基づいて画像信号を補正することを
特徴とするNMRイメージング方法によって達成され
る。
The above object of the present invention is to provide a static magnetic field, a gradient magnetic field, a high frequency magnetic field generating means, a detecting means for extracting an NMR signal from an object to be inspected, and an image reproduction for the detected signal. In an imaging method in an NMR inspection apparatus having a means for performing various calculations including the configuration, information on respiratory movements is measured together with the image signal during image signal measurement, and information on respiratory movements is obtained as phase information, and based on this. An NMR imaging method characterized by correcting the image signal by using a gradient magnetic field that does not affect a stationary object but only a moving object immediately after measuring the image signal, thereby performing respiratory motion. NMR imager, characterized in that information relating to Method or when measuring an image signal, information about the movement of the breath is measured at the same time so that the moving direction of the breath matches the phase encoding direction, and the image signal is corrected based on the measured data. Is achieved by an NMR imaging method.

【0005】[0005]

【作用】本発明に係るNMRイメージング方法において
は、撮影時、画像信号を計測する際、または、画像信号
を計測した直後に、呼吸に関する情報を含む信号を計測
するシーケンスを付加するようにしたものである。図4
〜図6に、その一例としての、撮影時、画像信号を計測
した直後に、呼吸に関する情報を含む信号を計測するた
めに付加するシーケンス例を示す。なお、図4〜図6中
の矢印で示した部分が、新しく付加された部分である。
また、ここでは、y方向に動きがあるものとしている。
図4においては、位相エンコードパルス58で同57を
キャンセルし、59で示す読出し用傾斜磁場をy方向に
印加する。このとき観測される信号63をフーリエ変換
すると、y軸に投影した画像信号が得られる。この信号
の端点位置から呼吸の位置を知ることができる。図5に
おいては、図4に示した例と同様に、位相エンコードパ
ルス78で同77をキャンセルし、読出し用傾斜磁場7
9をy方向に印加する。更に、上記位相位相エンコード
パルス78と、傾斜磁場79との間に、x方向の位相回
りを押さえるため、x方向傾斜磁場82を印加する。こ
のとき観測される信号84をフーリエ変換した信号の位
相を求めると、y方向の各位置における速度を求めるこ
とができる。また、各計測時に得られた上記データを積
分することにより、y軸に沿った呼吸の位置を知ること
ができる。
In the NMR imaging method according to the present invention, a sequence for measuring a signal containing respiration information is added at the time of photographing, when measuring an image signal, or immediately after measuring an image signal. Is. Figure 4
6 shows an example of a sequence added to measure a signal including information about respiration immediately after measuring an image signal at the time of photographing, as an example thereof. The parts indicated by arrows in FIGS. 4 to 6 are newly added parts.
In addition, here, it is assumed that there is movement in the y direction.
In FIG. 4, the phase encoding pulse 58 cancels 57, and a read gradient magnetic field indicated by 59 is applied in the y direction. When the signal 63 observed at this time is subjected to Fourier transform, an image signal projected on the y axis is obtained. The position of respiration can be known from the end point position of this signal. In FIG. 5, as in the example shown in FIG. 4, the phase encode pulse 78 cancels the pulse 77 and the read gradient magnetic field 7 is output.
9 is applied in the y direction. Further, an x-direction gradient magnetic field 82 is applied between the phase / phase encode pulse 78 and the gradient magnetic field 79 in order to suppress the phase rotation in the x direction. If the phase of the signal obtained by Fourier-transforming the signal 84 observed at this time is obtained, the velocity at each position in the y direction can be obtained. Further, by integrating the above data obtained at each measurement, the position of breathing along the y axis can be known.

【0006】図6においては、位相エンコードパルス9
8で同97をキャンセルし、y方向にフローエンコード
パルス99を印加し、読出し用傾斜磁場102をx方向
に印加する。このとき観測される信号104をフーリエ
変換した信号の位相を求めると、x方向の各位置におけ
る速度を求めることができる。また、図5に示した例と
同様に、各計測時に得られたデータを積分することによ
り、x軸に沿った呼吸の位置を知ることができる。な
お、位相エンコードパルス自身、1800前後で2回印加さ
れているためフローエンコードパルスになっており、図
6に示すパルス99を印加しなくても、動き(呼吸等)の
検出は可能であるが、傾斜磁場強度が計測毎に異なるた
め、呼吸による動きの速度に適した感度に設定するため
には、上述のフローエンコードパルス99があった方が
良い。ここで、位相エンコードによるフローエンコード
の量の補正は、ソフトウェアで行っても良いし、パルス
シーケンスの中に組込んでも良い。
In FIG. 6, the phase encode pulse 9
At 8, the same 97 is canceled, the flow encode pulse 99 is applied in the y direction, and the read gradient magnetic field 102 is applied in the x direction. When the phase of the signal obtained by Fourier-transforming the signal 104 observed at this time is obtained, the velocity at each position in the x direction can be obtained. Further, similarly to the example shown in FIG. 5, the position of breathing along the x-axis can be known by integrating the data obtained at each measurement. The phase encoding pulse itself, 180 0 has become flow encode pulse for being applied twice before and after, even without applying a pulse 99 shown in FIG. 6, it is possible detection of a motion (respiration, etc.) However, since the gradient magnetic field strength is different for each measurement, the flow encode pulse 99 described above is preferable in order to set the sensitivity suitable for the speed of movement due to respiration. Here, the correction of the amount of flow encoding by phase encoding may be performed by software or may be incorporated in the pulse sequence.

【0007】[0007]

【実施例】以下、本発明の実施例を図面に基づいて詳細
に説明する。図2は、本発明の一実施例に係るNMRイ
メージング装置の構成を示すブロック図である。図2に
おいて、21は被検体からNMR信号を出力させるため
に発生させる各種パルスおよび磁場をコントロールする
機能を有するシーケンス制御部、22は被検体の特定の
核種を共鳴させるために、高周波パルスを発生させる記
能を有する高周波パルス送信器、23は後述する磁場駆
動部24を制御して、NMR信号の共鳴周波数を決定す
る静磁場と、強さおよび方向を任意にコントロールでき
る傾斜磁場を発生させる機能を有する磁場制御部を示し
ている。また、24は上述の磁場制御部23から出力さ
れるコントロール信号に基づいて、計測に必要な磁場を
発生させる磁場制御部、25は被検体から発生するNM
R信号を、検波後、計測を行う機能を有する受信器、2
6は該受信器25から取込んだ計測信号を基に、画像再
構成および各種演算を行い、再構成された画像をCRT
ディスプレイ27に表示する機能を有する処理装置、2
7はCRTディスプレイを示している。
Embodiments of the present invention will now be described in detail with reference to the drawings. FIG. 2 is a block diagram showing the configuration of the NMR imaging apparatus according to the embodiment of the present invention. In FIG. 2, reference numeral 21 is a sequence control unit having a function of controlling various pulses and magnetic fields generated for outputting an NMR signal from the subject, and 22 is a high-frequency pulse for resonating a specific nuclide of the subject. A high-frequency pulse transmitter 23 having the function of controlling the magnetic field driving unit 24 described later to generate a static magnetic field that determines the resonance frequency of the NMR signal and a gradient magnetic field whose strength and direction can be arbitrarily controlled. 2 shows a magnetic field control unit having. Further, 24 is a magnetic field control unit that generates a magnetic field necessary for measurement based on the control signal output from the magnetic field control unit 23 described above, and 25 is an NM generated from the subject.
A receiver having a function of measuring the R signal and then measuring it, 2
Reference numeral 6 performs image reconstruction and various calculations based on the measurement signal fetched from the receiver 25, and displays the reconstructed image on the CRT.
A processing device having a function of displaying on the display 27, 2
Reference numeral 7 indicates a CRT display.

【0008】上述の如く構成された本実施例の動作を、
以下、図7に示した処理フローチャートに従って説明す
る。 ステップ111:図4のパルスシーケンスにおいて、位
相エンコードパルス57を順次変化させて、画像再生に
必要な回数分、以下のシーケンスで位置検出用信号63
を計測する。 ステップ112:位相エンコードパルス57と同一の大
きさを有するパルス58を印加して、位相エンコードパ
ルス57の影響をキャンセルする。ここでは、1800パル
ス52が印加されているため、パルス58でパルス57
の影響をキャンセルすることができるものである。 ステップ113:y方向読出し用傾斜磁場59を印加
し、位置検出用信号63を計測する。 ステップ114:計測信号63をフーリエ変換し、y軸
に投影したデータを得る。 ステップ115:上記投影データの上端点を検出する。
各データの上端点から呼吸の位置が検出できる。
The operation of this embodiment constructed as described above is
Hereinafter, description will be given according to the processing flowchart shown in FIG. 7. Step 111: In the pulse sequence of FIG. 4, the phase encode pulse 57 is sequentially changed, and the position detection signal 63 is generated by the number of times required for image reproduction in the following sequence.
To measure. Step 112: A pulse 58 having the same magnitude as the phase encode pulse 57 is applied to cancel the influence of the phase encode pulse 57. Here, since the 180 0 pulse 52 is applied, the pulse 58 is replaced by the pulse 57.
The effect of can be canceled. Step 113: The gradient magnetic field 59 for reading in the y direction is applied, and the signal 63 for position detection is measured. Step 114: Fourier transform the measurement signal 63 to obtain data projected on the y-axis. Step 115: Detect the upper end point of the projection data.
The position of respiration can be detected from the upper end point of each data.

【0009】次に、図5に示したシーケンスを用いる方
法について説明する。図4に示したシーケンスと異なる
のは、x方向傾斜磁場82が印加される点である。これ
は、傾斜磁場80と81の影響で、x方向の動きによる
位相変化が生ずるため、それをキャンセルするために行
うものである。このとき計測される信号84をフーリエ
変換して位相を求めると、各y方向の位置における速度
を得ることができることは前述の通りである。これは、
位相と速度とが比例関係にあるためである。なお、この
点に関しては、本出願人が先に提案した特願昭60−1
50194号「NMR血流イメージング方式」の明細書の
記載を参照されたい。この速度を各計測データについて
積分すると位置を求めることができ、呼吸の位置を検出
することができる。また、フーリエ変換後のデータの絶
対値は、図4に示したシーケンスで得られる投影データ
と等しいので、図7に示した方法を併用し、精度を高め
ることができる。
Next, a method using the sequence shown in FIG. 5 will be described. The difference from the sequence shown in FIG. 4 is that the x-direction gradient magnetic field 82 is applied. This is done to cancel the phase change due to the movement in the x direction due to the influence of the gradient magnetic fields 80 and 81. As described above, the velocity at each position in the y direction can be obtained by Fourier-transforming the signal 84 measured at this time to obtain the phase. this is,
This is because the phase and the speed have a proportional relationship. Regarding this point, Japanese Patent Application No. 60-1 previously proposed by the applicant.
See the description in the specification of 50194 "NMR Blood Flow Imaging Method". The position can be obtained by integrating this velocity for each measurement data, and the position of respiration can be detected. Further, since the absolute value of the data after the Fourier transform is equal to the projection data obtained by the sequence shown in FIG. 4, the method shown in FIG. 7 can be used together to improve the accuracy.

【0010】次に、図6に示したシーケンスについて、
図8の処理フローチャートに従って説明する。 ステップ121:図6のパルスシーケンスにおいて、位
相エンコードパルス97を順次変化させて、画像再生に
必要な回数分、以下のシーケンスで位置検出用信号11
4を計測する。 ステップ122:位相エンコードパルス97と同一の大
きさを有するパルス98を印加し、位相エンコードパル
ス97の影響をキャンセルする。 ステップ123:y方向フローエンコードパルス99を
印加し、y方向の動きの速度に応じた位相変化を与え
る。 ステップ124:x方向読出し用傾斜磁場102を印加
し、速度検出用信号104を計測する。 ステップ125:計測信号104のフーリエ変換を行
う。 ステップ126:フーリエ変換後データの位相成分を求
め、各x方向の位置における速度を得る。以上、位置の
検出から呼吸の状態を知ることができる。
Next, regarding the sequence shown in FIG.
A description will be given according to the processing flowchart of FIG. Step 121: In the pulse sequence of FIG. 6, the phase encode pulse 97 is sequentially changed, and the position detection signal 11 is generated by the following sequence for the number of times necessary for image reproduction.
Measure 4. Step 122: A pulse 98 having the same magnitude as the phase encode pulse 97 is applied to cancel the influence of the phase encode pulse 97. Step 123: The y-direction flow encode pulse 99 is applied to give a phase change according to the speed of movement in the y-direction. Step 124: The x-direction reading gradient magnetic field 102 is applied, and the velocity detection signal 104 is measured. Step 125: Fourier transform the measurement signal 104. Step 126: Obtain the phase component of the Fourier transformed data to obtain the velocity at each position in the x direction. As described above, the state of breathing can be known from the detection of the position.

【0011】図9(a)〜(c)は、それぞれ、図4〜図6
に示したシーケンスに対応するものである。図9(a)
は、図4に示したシーケンスにおいて、観測される信号
63をフーリエ変換すると、y軸に投影した画像信号1
33が得られる。この信号の端点位置から、呼吸の位置
を知ることができる(矢印参照)。図9(b)は、図5に示
したシーケンスにおいて、観測される信号84をフーリ
エ変換した信号の位相を求めると、y方向の各位置にお
ける速度を求めることができる(矢印参照)。図9(c)
は、図6に示したシーケンスにおいて、観測される信号
104をフーリエ変換した信号の位相を求めると、x方
向の各位置における速度を求めることができる(矢印参
照)。上記実施例によれば、特殊なハードウェアを用い
ることなく、NMRイメージング装置本体のみで信号計
測時の呼吸の動きを検出し、呼吸と同期をとることがで
き、経済的に、高画質の撮影を行うことができるという
効果が得られる。
FIGS. 9A to 9C show FIGS. 4 to 6, respectively.
It corresponds to the sequence shown in. Figure 9 (a)
Is the image signal 1 projected on the y-axis when the observed signal 63 is Fourier transformed in the sequence shown in FIG.
33 is obtained. The position of the breath can be known from the end point position of this signal (see arrow). In FIG. 9B, when the phase of the signal obtained by Fourier-transforming the observed signal 84 in the sequence shown in FIG. 5 is obtained, the velocity at each position in the y direction can be obtained (see arrow). Figure 9 (c)
In the sequence shown in FIG. 6, if the phase of the signal obtained by Fourier-transforming the observed signal 104 is obtained, the velocity at each position in the x direction can be obtained (see arrow). According to the above-described embodiment, it is possible to detect the movement of respiration at the time of signal measurement only with the main body of the NMR imaging device without using special hardware, and to synchronize with the respiration, which allows economical and high-quality imaging. The effect of being able to do is obtained.

【0012】以下、本発明の他の実施例として、画像信
号を計測した直後に、動きに関する情報を含む信号を計
測するシーケンスを付加するとともに、該計測情報に基
づいて計測画像信号を補正する方式について説明する。
信号の計測方式としては、先に図4〜図6に示したシー
ケンスを用いるものとする。なお、図4に示したシーケ
ンスはy軸に投影した画像信号を得る方式、他はy方向
およびx方向の速度を得る方式であり、速度からは積分
することにより位置を知ることが可能である。まず、図
10に従って、検出信号の補正の原理を説明する。今、
図10に示す如く、y方向だけに動きがある場合を考え
る。y方向だけなので、以後、一次元データとして説明
する。息を吸って、最も腹部が伸びたときの状態が斜線
で示されているとすると、このときの投影データは、図
10(b)の斜線部のようになる。このときの、上端と下
端との間隔をL、画像中心と最下端との距離をhとす
る。また、任意の時刻tの、息を吸ったときの状態を破
線で示し、それをL/a(t)で表わす。
As another embodiment of the present invention, a method of adding a sequence for measuring a signal including information about motion immediately after measuring an image signal and correcting the measured image signal based on the measurement information Will be described.
As the signal measuring method, the sequence shown in FIGS. 4 to 6 is used. The sequence shown in FIG. 4 is a system for obtaining an image signal projected on the y-axis, and the other is a system for obtaining the velocities in the y-direction and the x-direction, and the position can be known by integrating from the velocity. . First, the principle of correction of the detection signal will be described with reference to FIG. now,
Consider a case where there is a movement only in the y direction as shown in FIG. Since it is only in the y direction, it will be described as one-dimensional data hereinafter. Assuming that the state in which the abdomen is extended most by inhalation is indicated by the diagonal lines, the projection data at this time is as shown by the hatched portions in FIG. At this time, the distance between the upper end and the lower end is L, and the distance between the image center and the lowermost end is h. Further, a state when breathing in at an arbitrary time t is shown by a broken line and is represented by L / a (t).

【0013】このとき、斜線の投影データに相当する計
測信号、すなわち、投影データのフーリエスペクトルを
F(ω)とすると、破線部の計測信号G(ω)は、
At this time, if the measurement signal corresponding to the projection data of the diagonal line, that is, the Fourier spectrum of the projection data is F (ω), the measurement signal G (ω) in the broken line portion is

【外1】 と表わされる。従って式(1)中の変数a(t),hの値が
わかっていれば、計測信号G(ω)からF(ω)を求めるこ
とができる。この処理は、上述の式(1)より、位相変化
[Outer 1] Is represented. Therefore, if the values of the variables a (t) and h in equation (1) are known, F (ω) can be obtained from the measurement signal G (ω). This process is based on the above equation (1)

【外2】 の補正と、[Outside 2] Correction of

【外3】 からF(ω)を求めるリサンプリング処理の2つから成る
ことがわかる。つまり、hとa(t)の値がわかれば、補
正ができる。
[Outside 3] Therefore, it can be seen that it consists of two resampling processes for obtaining F (ω). That is, if the values of h and a (t) are known, the correction can be performed.

【0014】ここで、hの値は通常の画像再生を行うこ
とで簡単に求められ、a(t)の値は先に示した図4〜図
6のシーケンスで計測できる。以上の処理を、実際の画
像の計測信号レベルで示すと、図11のようになる。y
方向のみに動きがあるものとすれば、計測信号は
Here, the value of h can be easily obtained by performing normal image reproduction, and the value of a (t) can be measured by the sequences of FIGS. 4 to 6 shown above. The above processing is shown in FIG. 11 in terms of actual image measurement signal levels. y
If there is movement only in the direction, the measurement signal is

【外4】 となる。x方向の動きがないので、横方向にフーリエ変
換すると、
[Outside 4] Becomes Since there is no movement in the x direction, Fourier transform in the horizontal direction

【外5】 となる。従来の方法では、このデータをそのまま、縦方
向にフーリエ変換し、画像I(x,y)’が得られる。こ
の画像には、動きの影響が含まれており、画質が劣化し
ている。
[Outside 5] Becomes In the conventional method, this data is directly subjected to Fourier transform in the vertical direction to obtain an image I (x, y) '. This image contains the influence of motion, and the image quality is degraded.

【0015】これに対して、本実施例においては、横方
向にフーリエ変換した後、位相補正とリサンプリング処
理から成る体動補正処理を行い、動きを含まない計測信
号F(x,ωy)を求めた後、縦方向フーリエ変換を行っ
て、動きの影響を補正した画像I(x,y)を求める。以
下、具体例により説明を続ける。本実施例においては、
図2に示した如き装置を用いて、図12に示す処理フロ
ーチャートに従って処理を行うものとする。 ステップ151:図4のパルスシーケンスにおいて、位
相エンコードパルス57,58を順次変化させて、画像
再生に必要な回数分、画像計測信号62および位置検出
用信号63を計測する。 ステップ152:前ステップで得た画像計測信号から、
そのまま再生した画像を得て、図10(b)のパラメータ
hの値を求める。また、位置検出用信号63をフーリエ
変換し、各信号計測時のy軸への投影データを求め、息
を最も吸った状態の長さLを求める。そのLから、各計
測時点a(t)を求める。以上で、補正に必要なパラメー
タh,L,a(t)が得られる。
On the other hand, in the present embodiment, after the Fourier transform in the horizontal direction, the body motion correction process including the phase correction and the resampling process is performed, and the measurement signal F (x, ω y ) containing no motion is obtained. Then, a vertical direction Fourier transform is performed to obtain an image I (x, y) in which the influence of motion is corrected. Hereinafter, the description will be continued by using a specific example. In this embodiment,
It is assumed that processing is performed according to the processing flow chart shown in FIG. 12 by using the device as shown in FIG. Step 151: In the pulse sequence of FIG. 4, the phase encode pulses 57 and 58 are sequentially changed to measure the image measurement signal 62 and the position detection signal 63 by the number of times required for image reproduction. Step 152: From the image measurement signal obtained in the previous step,
The reproduced image is obtained as it is, and the value of the parameter h in FIG. Further, the position detection signal 63 is Fourier-transformed to obtain projection data on the y-axis at the time of measuring each signal, and the length L of the most inhaled state is obtained. Each measurement time point a (t) is obtained from the L. As described above, the parameters h, L, a (t) necessary for the correction can be obtained.

【0016】ステップ153:計測信号を横方向にフー
リエ変換して、前記式(4)の信号を求める。 ステップ154:各計測信号毎に、a(t),hの値から
次式の値を求め、上記式(4)の値に掛ける。
Step 153: Fourier-transform the measurement signal in the horizontal direction to obtain the signal of the above equation (4). Step 154: For each measurement signal, the value of the following equation is obtained from the values of a (t) and h, and is multiplied by the value of the above equation (4).

【外6】 このとき、次の位相補正された信号が得られる。[Outside 6] At this time, the next phase-corrected signal is obtained.

【外7】 ステップ155:上記式(6)の値から、縦方向であるω
y軸方向に リサンプリング処理により、(x,ωy)を求
める。ここで リサンプリング処理としては、どのよう
な補間方法を用いても良く、例えば、線形補間あるいは
スプライン補間が適用できる。
[Outside 7] Step 155: From the value of the above formula (6), it is ω which is the vertical direction.
(x, ω y ) is obtained by resampling processing in the y- axis direction. Here, any interpolation method may be used as the resampling processing, and for example, linear interpolation or spline interpolation can be applied.

【0017】ステップ156:位相補正とリサンプリン
グ処理を行った(x,ωy)に対して、縦方向にフーリエ
変換を行うことにより、呼吸による動きの影響を除去し
た画像を得る。上記実施例によれば、特殊なハードウェ
アを用いることなく、NMRイメージング装置本体のみ
で信号計測時の呼吸の動きを検出し、呼吸と同期をとる
ことなしに、経済的に、高画質の撮影を行うことができ
るという効果がある。なお、上記各実施例は本発明の一
例を示したものであり、本発明はこれに限定されるべき
ものではないことは言うまでもないことである。図3
は、本発明に関連する技術の例としての、撮影に先立っ
て短時間内に呼吸を検出するシーケンス例を示すもので
ある。ここでは、y方向に動きがあるものとする。図の
破線で囲まれた部分が、この技術の特徴部分である。
Step 156: An image in which the influence of motion due to respiration is removed is obtained by performing Fourier transform in the vertical direction on (x, ω y ) for which phase correction and resampling processing have been performed. According to the above-described embodiment, the motion of the respiration at the time of signal measurement is detected only by the NMR imaging apparatus main body without using any special hardware, and the high-quality imaging can be economically performed without synchronizing with the respiration. There is an effect that can be done. It is needless to say that each of the above embodiments is an example of the present invention and the present invention should not be limited to this. Figure 3
FIG. 4 shows an example of a sequence for detecting breathing within a short time prior to imaging, as an example of a technique related to the present invention. Here, it is assumed that there is movement in the y direction. The part surrounded by the broken line in the figure is a characteristic part of this technique.

【0018】この技術では、短時間内に何回も繰り返し
て、信号を計測する。計測信号をフーリエ変換すること
によって、呼吸による動きの位置を振幅情報から、速度
を位相情報から、それぞれ、得ることができる。RFパ
ルス31で、スピンをα0(αは90〜180)倒す。1800に近
い程、緩和が早くなるので、繰り返し時間を短縮できる
が、計測信号が小さくなる。その後、1800パルス32を
印加し、y方向の傾斜磁場Gy40を印加しながら、信
号44を計測する。この信号から、呼吸による動きの位
置や速度を検出し、呼吸と同期して撮影するタイミング
を調べる。このタイミングとずれている場合には、該当
するタイミングまで、このシーケンスを繰り返し計測し
て、待つ。上述の呼吸の状態の検知は、次に述べる如き
原理による。
In this technique, the signal is measured repeatedly within a short time. By Fourier transforming the measurement signal, it is possible to obtain the position of the motion due to respiration from the amplitude information and the velocity from the phase information. The RF pulse 31 defeats the spin α 0 (α is 90 to 180). Closer to 180 0, since relaxation is faster and can shorten the repetition time, the measurement signal is reduced. Thereafter, 180 0 pulse 32 is applied, while applying a gradient magnetic field G y 40 in the y-direction, to measure the signal 44. From this signal, the position and speed of movement due to respiration are detected, and the timing of photographing in synchronization with respiration is examined. If there is a deviation from this timing, this sequence is repeatedly measured and waits until the corresponding timing. The detection of the above-mentioned breathing state is based on the principle as described below.

【0019】上記計測信号44をフーリエ変換したデー
タは、画像信号を、y軸上に投影したデータとなってい
る。振幅は投影値であり、位相はその各投影データのy
方向の動きの速度に比例した値となっている。従って、
位置は投影データの端点位置を検出することにより、ま
た、速度は位相値から知ることができる。例えば、図3
の破線内のシーケンスを100msec単位で繰り返せば、最
大100msecの遅れで検出できる。上述の技術を、図2に
示したNMRイメージング装置に適用した場合の動作に
ついて、図1に示したフローチャートにより説明してお
く。なお、ここでは、y方向に呼吸性の動きがある場合
について説明するが、他の方向の場合でも、同様に考え
て良い。 ステップ11:図3のパルスシーケンスにおいて、位相
エンコードパルス41を順次変化させて、画像再生に必
要な回数分、ステップ12〜同17まで、以下のシーケ
ンスで、画像信号を計測する。
The Fourier-transformed data of the measurement signal 44 is data obtained by projecting the image signal on the y-axis. The amplitude is the projection value and the phase is the y of each projection data.
The value is proportional to the speed of movement in the direction. Therefore,
The position can be known by detecting the end position of the projection data, and the velocity can be known from the phase value. For example, in FIG.
If the sequence within the broken line is repeated every 100 msec, it can be detected with a maximum delay of 100 msec. The operation when the above-described technique is applied to the NMR imaging apparatus shown in FIG. 2 will be described with reference to the flowchart shown in FIG. Note that, here, a case where there is a respiratory movement in the y direction will be described, but the same may be applied to the case of other directions. Step 11: In the pulse sequence of FIG. 3, the phase encode pulse 41 is sequentially changed, and the image signal is measured by the following sequence from step 12 to step 17 for the number of times required for image reproduction.

【0020】ステップ12:y方向の動きを伴なう呼吸
を検出するシーケンスで、信号44を計測する。すなわ
ち、α0パルス(90≦α<180)31と1800パルス32と
を、傾斜磁場Gz35,36とともに印加する。信号計測
時には 斜磁場Gy40を同時に印加する。 ステップ13:計測信号44をフーリエ変換し、フーリ
エ変換後のデータの振幅と位相を求める。 ステップ14:振幅データは、画像信号のy軸への投影
したデータとなっているので、その値が“0”になる端
点の位置を検出する。 ステップ15:位相データは、投影データの各y座標点
上の位相を示しているので、速度と位相との関係から、
速度を算出する。
Step 12: The signal 44 is measured in a sequence for detecting respiration accompanied by movement in the y direction. That is, the α 0 pulse (90 ≦ α <180) 31 and the 180 0 pulse 32 are applied together with the gradient magnetic fields G z 35 and 36. The gradient magnetic field G y 40 is applied simultaneously during signal measurement. Step 13: Fourier transform the measurement signal 44 to obtain the amplitude and phase of the data after the Fourier transform. Step 14: Since the amplitude data is the data obtained by projecting the image signal onto the y-axis, the position of the end point where the value becomes "0" is detected. Step 15: Since the phase data indicates the phase on each y coordinate point of the projection data, from the relationship between the velocity and the phase,
Calculate speed.

【0021】ステップ16:予め定めた撮影すべき呼吸
の位置か否か、速度と位置から調べる。もし、該当位置
でなければ、ステップ12に戻って、繰り返す。但し、
緩和時間の関係で、直ちに繰り返しても信号が出ないの
で、間を少しあける。また、もし、該当位置であれば、
ステップ17に進む。 ステップ17:通常のパルスシーケンスで、撮影を行
い、画像信号45を計測する。 ステップ18:呼吸に同期して、繰り返し得られた画像
信号に基づき、画像再生を行う。 以上の処理手順を経て得られた画像は、呼吸に同期して
いるため、患者の動きに伴なうムービング・アーチファ
クト等の、画像劣化を生ずることがなく、高画質の画像
が得られる。
Step 16: It is checked from the speed and position whether or not the position of the breath to be imaged is determined in advance. If it is not the corresponding position, the process returns to step 12 and is repeated. However,
Due to the relaxation time, there is no signal even if it is repeated immediately, so leave some space. Also, if the corresponding position,
Proceed to step 17. Step 17: Photographing is performed by the normal pulse sequence, and the image signal 45 is measured. Step 18: Image reproduction is performed in synchronization with respiration based on the image signal obtained repeatedly. The image obtained through the above processing procedure is synchronized with respiration, so that a high-quality image can be obtained without causing image deterioration such as moving artifact due to the movement of the patient.

【0022】[0022]

【発明の効果】以上、詳細に説明した如く、本発明によ
れば、NMRイメージング装置本体のみで呼吸の動きを
検出し、その影響を補正可能とするNMRイメージング
方法を実現できるという顕著な効果を奏するものであ
る。
As described above in detail, according to the present invention, it is possible to realize a remarkable effect that an NMR imaging method capable of detecting respiratory movements only by the main body of the NMR imaging apparatus and correcting the influence thereof can be realized. It plays.

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

【図1】本発明に関連する技術の例としての、撮影に先
立って短時間内に呼吸を検出する処理のフローチャート
である。
FIG. 1 is a flowchart of a process of detecting respiration within a short time before capturing an image, as an example of a technique related to the present invention.

【図2】本発明の一実施例に係る、NMRイメージング
装置の構成を示すブロック図である。
FIG. 2 is a block diagram showing a configuration of an NMR imaging apparatus according to an embodiment of the present invention.

【図3】本発明に関連する技術の例としての、撮影に先
立って短時間内に呼吸を検出するシーケンス例を示す図
である。
FIG. 3 is a diagram showing an example of a sequence for detecting breathing within a short time prior to imaging, as an example of a technique related to the present invention.

【図4】一実施例に係る、撮影時、画像信号を計測した
直後に、呼吸に関する情報を含む信号を計測するために
付加するシーケンス例を示す図(その1)である。
FIG. 4 is a diagram (No. 1) showing an example of a sequence added to measure a signal including information related to respiration immediately after measuring an image signal during photographing according to an embodiment.

【図5】一実施例に係る、撮影時、画像信号を計測した
直後に、呼吸に関する情報を含む信号を計測するために
付加するシーケンス例を示す図(その2)である。
FIG. 5 is a diagram (No. 2) showing an example of a sequence added to measure a signal including information about respiration immediately after measuring an image signal during photographing according to an embodiment.

【図6】一実施例に係る、撮影時、画像信号を計測した
直後に、呼吸に関する情報を含む信号を計測するために
付加するシーケンス例を示す図(その3)である。
FIG. 6 is a diagram (No. 3) showing an example of a sequence added to measure a signal including information about respiration immediately after measuring an image signal during photographing according to an embodiment.

【図7】実施例の動作を示す処理フローチャート(その
1)である。
FIG. 7 is a processing flowchart (part 1) showing the operation of the embodiment.

【図8】実施例の動作を示す処理フローチャート(その
2)である。
FIG. 8 is a processing flowchart (part 2) showing the operation of the embodiment.

【図9】図4〜図6に示したシーケンスの概念を説明す
る図である。
FIG. 9 is a diagram for explaining the concept of the sequences shown in FIGS.

【図10】動きの補正に必要なパラメータの説明図であ
る。
FIG. 10 is an explanatory diagram of parameters required for motion correction.

【図11】他の実施例に係る、処理の概念を示す図であ
る。
FIG. 11 is a diagram showing a concept of processing according to another embodiment.

【図12】図11に示した処理の詳細を示すフローチャ
ートである。
FIG. 12 is a flowchart showing details of the processing shown in FIG.

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

21 シーケンス制御部 22 高周波パルス送信部 23 磁場制御部 24 磁場駆動部 25 受信器 26 処理装置 27 CRTディスプレイ 57,58,77,78,97,98 位相エンコードパル
ス 63 画像信号 84,104 速度信号
21 sequence control unit 22 high frequency pulse transmission unit 23 magnetic field control unit 24 magnetic field drive unit 25 receiver 26 processing device 27 CRT display 57,58,77,78,97,98 phase encode pulse 63 image signal 84,104 velocity signal

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】 静磁場,傾斜磁場,高周波磁場の発生手
段と、検査対象物からのNMR信号を取出す検出手段
と、検出された信号に対し画像再構成を含む各種演算を
行う手段とを有するNMR検査装置におけるイメージン
グ方法において、画像信号計測時に、呼吸の動きに関す
る情報を画像信号とともに計測し、呼吸の動きに関する
情報を位相情報として得て、これに基づいて前記画像信
号を補正することを特徴とするNMRイメージング方
法。
1. A static magnetic field, a gradient magnetic field, a high-frequency magnetic field generating means, a detecting means for extracting an NMR signal from an inspection object, and a means for performing various calculations including image reconstruction on the detected signal. In the imaging method of the NMR inspection apparatus, when measuring an image signal, information relating to respiratory movement is measured together with the image signal, information relating to respiratory movement is obtained as phase information, and the image signal is corrected based on this. And an NMR imaging method.
【請求項2】 静磁場,傾斜磁場,高周波磁場の発生手
段と、検査対象物からのNMR信号を取出す検出手段
と、検出された信号に対し画像再構成を含む各種演算を
行う手段とを有するNMR検査装置におけるイメージン
グ方法において、画像信号を計測した直後に、静止物体
には影響を与えず、動く物体にのみ影響を与える傾斜磁
場を印加して、呼吸の動きに関する情報を、速度に応じ
て異なる信号になるように計測し、これに基づいて前記
画像信号を補正することを特徴とするNMRイメージン
グ方法。
2. A static magnetic field, a gradient magnetic field, a high frequency magnetic field generating means, a detecting means for extracting an NMR signal from an inspection object, and a means for performing various calculations including image reconstruction on the detected signal. In an imaging method in an NMR inspection apparatus, immediately after measuring an image signal, a gradient magnetic field that does not affect a stationary object but only a moving object is applied to obtain information on respiratory motion according to a velocity. An NMR imaging method characterized by measuring so as to obtain different signals and correcting the image signal based on the measurement.
【請求項3】 前記計測信号を観測する際、読出し用傾
斜磁場を、呼吸の動きを検出したい方向に対して垂直な
方向に印加することを特徴とする、請求項2記載のNM
Rイメージング方法。
3. The NM according to claim 2, wherein when observing the measurement signal, a read gradient magnetic field is applied in a direction perpendicular to a direction in which respiratory movement is desired to be detected.
R imaging method.
【請求項4】 静磁場,傾斜磁場,高周波磁場の発生手
段と、検査対象物からのNMR信号を取出す検出手段
と、検出された信号に対し画像再構成を含む各種演算を
行う手段とを有するNMR検査装置におけるイメージン
グ方法において、画像信号を計測する際に、呼吸の動き
に関する情報を、呼吸で動く方向が位相エンコード方向
と一致するようにして同時期に計測し、該計測データに
基づいて画像信号を補正することを特徴とするNMRイ
メージング方法。
4. A static magnetic field, a gradient magnetic field, a high-frequency magnetic field generating means, a detecting means for extracting an NMR signal from an inspection object, and a means for performing various calculations including image reconstruction on the detected signal. In the imaging method in the NMR inspection apparatus, when measuring an image signal, information about the movement of the breath is measured at the same time so that the moving direction of the breath matches the phase encoding direction, and an image is obtained based on the measured data. A method for NMR imaging, which comprises correcting a signal.
【請求項5】 前記画像信号の補正を、呼吸の動きによ
る影響で生じた計測信号上の位相成分の変動に関して行
うことを特徴とする、請求項4記載のNMRイメージン
グ方法。
5. The NMR imaging method according to claim 4, wherein the correction of the image signal is performed with respect to a fluctuation of a phase component on the measurement signal caused by the influence of respiratory movement.
【請求項6】 前記画像信号の呼吸の動きによる影響で
生じた計測位置のずれを、リサンプリング処理により補
正することを特徴とする、請求項4または5記載のNM
Rイメージング方法。
6. The NM according to claim 4, wherein the displacement of the measurement position caused by the influence of the respiratory movement of the image signal is corrected by resampling processing.
R imaging method.
JP6208457A 1994-09-01 1994-09-01 NMR imaging method Pending JPH07163537A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP6208457A JPH07163537A (en) 1994-09-01 1994-09-01 NMR imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP6208457A JPH07163537A (en) 1994-09-01 1994-09-01 NMR imaging method

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP61207055A Division JPS63164943A (en) 1986-09-03 1986-09-03 NMR imaging method

Publications (1)

Publication Number Publication Date
JPH07163537A true JPH07163537A (en) 1995-06-27

Family

ID=16556514

Family Applications (1)

Application Number Title Priority Date Filing Date
JP6208457A Pending JPH07163537A (en) 1994-09-01 1994-09-01 NMR imaging method

Country Status (1)

Country Link
JP (1) JPH07163537A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009515586A (en) * 2005-11-10 2009-04-16 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Signal processing unit for generating images
JP2009160378A (en) * 2007-12-10 2009-07-23 Toshiba Corp Magnetic resonance imaging system
KR20160018989A (en) * 2014-08-08 2016-02-18 경북대학교 산학협력단 Apparatus and method for predicting respiratory motion
JP2018047160A (en) * 2016-09-23 2018-03-29 キヤノンメディカルシステムズ株式会社 Magnetic resonance imaging apparatus

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6120541A (en) * 1984-07-05 1986-01-29 株式会社島津製作所 Method for correcting distortion of nuclear magnetic resonance images due to respiratory body movements
JPS61133850A (en) * 1984-11-29 1986-06-21 エヌ・ベー・フイリツプス・フルーイランペンフアブリケン Nuclear magnetic resonance distribution determining method and device
JPS61143035A (en) * 1984-11-21 1986-06-30 ゼネラル・エレクトリツク・カンパニイ Method and apparatus for determining breathing cycle of patient

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6120541A (en) * 1984-07-05 1986-01-29 株式会社島津製作所 Method for correcting distortion of nuclear magnetic resonance images due to respiratory body movements
JPS61143035A (en) * 1984-11-21 1986-06-30 ゼネラル・エレクトリツク・カンパニイ Method and apparatus for determining breathing cycle of patient
JPS61133850A (en) * 1984-11-29 1986-06-21 エヌ・ベー・フイリツプス・フルーイランペンフアブリケン Nuclear magnetic resonance distribution determining method and device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009515586A (en) * 2005-11-10 2009-04-16 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Signal processing unit for generating images
JP2009160378A (en) * 2007-12-10 2009-07-23 Toshiba Corp Magnetic resonance imaging system
KR20160018989A (en) * 2014-08-08 2016-02-18 경북대학교 산학협력단 Apparatus and method for predicting respiratory motion
JP2018047160A (en) * 2016-09-23 2018-03-29 キヤノンメディカルシステムズ株式会社 Magnetic resonance imaging apparatus

Similar Documents

Publication Publication Date Title
EP2628016B1 (en) Mr data acquisition using physiological monitoring
US7221973B2 (en) Magnetic resonance tomography method and apparatus with motion correction
US9684979B2 (en) MRI 3D cine imaging based on intersecting source and anchor slice data
JPH11113878A (en) Magnetic resonance imaging method
US6889071B2 (en) Acquisition of high-temporal free-breathing MR images
US4905699A (en) Method and apparatus for NMR imaging
WO2006064582A1 (en) Magnetic resonance imaging system
JPWO2004080301A1 (en) Magnetic resonance imaging system
JPH11113877A (en) Magnetic resonance imaging method
JP3993799B2 (en) Magnetic resonance imaging system
JP4138258B2 (en) Magnetic resonance imaging system
JP4343317B2 (en) Magnetic resonance imaging system
JPH064066B2 (en) Magnetic resonance imaging equipment
JP2641486B2 (en) NMR imaging equipment
JP5454846B2 (en) Magnetic resonance imaging system
JP2005040416A (en) Magnetic resonance imaging apparatus
CN1312631C (en) Evaluation of diagnosis process according to time changing
JPH07163537A (en) NMR imaging method
JPH10201737A (en) Mr method for monitoring physical displacement with image assistance and mr system for executing this method
JPWO2000072752A1 (en) Magnetic resonance imaging method and apparatus
JP2595006B2 (en) MR imaging method
JPH06327649A (en) Image reconstruction method in nuclear magnetic resonance imaging apparatus
US20100292562A1 (en) Method for generating MR (Magnetic resonance) images of a moving partial area of an object
JPH09248286A (en) Inspection equipment using nuclear magnetic resonance
JP2003325477A (en) Magnetic resonance imaging apparatus