JP2001112735A - Magnetic resonance imaging apparatus - Google Patents
Magnetic resonance imaging apparatusInfo
- Publication number
- JP2001112735A JP2001112735A JP29657299A JP29657299A JP2001112735A JP 2001112735 A JP2001112735 A JP 2001112735A JP 29657299 A JP29657299 A JP 29657299A JP 29657299 A JP29657299 A JP 29657299A JP 2001112735 A JP2001112735 A JP 2001112735A
- Authority
- JP
- Japan
- Prior art keywords
- phase
- echo signal
- measurement
- magnetic field
- echo
- 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
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 28
- 238000005259 measurement Methods 0.000 claims abstract description 128
- 238000012937 correction Methods 0.000 claims description 51
- 238000000034 method Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 16
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 15
- 230000003068 static effect Effects 0.000 claims description 8
- 230000001678 irradiating effect Effects 0.000 claims description 5
- 238000002592 echocardiography Methods 0.000 abstract description 26
- 238000005070 sampling Methods 0.000 abstract description 22
- 238000003384 imaging method Methods 0.000 abstract description 6
- 238000012935 Averaging Methods 0.000 abstract description 2
- 230000008030 elimination Effects 0.000 abstract 1
- 238000003379 elimination reaction Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 4
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、被検体中の水素や
燐等からの核磁気共鳴(以下、「NMR」という)信号を測
定し、核の密度分布や緩和時間分布等を映像化する磁気
共鳴イメージング(MRI)装置に関し、特にエコー信号間
の位相変化により生じるアーチファクトを効果的に抑制
することを可能としたMRI装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention measures a nuclear magnetic resonance (hereinafter, referred to as "NMR") signal from hydrogen, phosphorus, or the like in a subject, and visualizes a nuclear density distribution, a relaxation time distribution, and the like. The present invention relates to a magnetic resonance imaging (MRI) apparatus, and more particularly to an MRI apparatus capable of effectively suppressing an artifact caused by a phase change between echo signals.
【0002】[0002]
【従来の技術】MRIを用いた撮像方法では、一般に読み
出し傾斜磁場の反転を利用してNMR信号をエコーとして
計測する方法が用いられる。エコー信号は一定のサンプ
リング時間計測して所定サンプリング数の時系列データ
とする。ところでエコー信号は、例えば正極性の読み出
し傾斜磁場の印加時間と印加強度の積(積分値)と負極
性の読み出し傾斜磁場の印加時間と印加強度の積とが一
致するところで発生する、即ち両者が一致した時点でエ
コーはピークとなるので、サンプリング時間はこのピー
クとなる時点を中心として設定される。2. Description of the Related Art In an imaging method using MRI, a method of measuring an NMR signal as an echo by using inversion of a readout gradient magnetic field is generally used. The echo signal is measured for a fixed sampling time to obtain time-series data of a predetermined sampling number. An echo signal is generated when, for example, the product (integral value) of the application time and the applied intensity of the read gradient magnetic field of the positive polarity and the product of the applied time and the applied intensity of the read gradient magnetic field of the negative polarity coincide with each other. Since the echo becomes a peak at the time of the coincidence, the sampling time is set around the time of the peak.
【0003】しかし、MRI装置ではこのように計測した
エコー信号が計測空間の中央からずれる場合がある。こ
れは永久磁石や超電導磁石等の静磁場発生装置によって
高度の静磁場均一性が保たれていても、測定対象である
被検体の磁化率によって磁場のオフセットが変化するこ
とや、傾斜磁場パルスの不完全性や出力応答の誤差等に
よる。[0003] However, in the MRI apparatus, the echo signal thus measured may be shifted from the center of the measurement space. This is because even if a high static magnetic field uniformity is maintained by a static magnetic field generator such as a permanent magnet or a superconducting magnet, the magnetic field offset changes due to the magnetic susceptibility of the object to be measured, This is due to imperfections or errors in output response.
【0004】撮像方法として1回或いは数回のRF照射で
1枚の画像再構成に必要なエコー信号を取得するシング
ルショット或いはマルチショットのエコープレナーイメ
ージング(EPI)法がある。この撮像方法では、読み出
し傾斜磁場を反転しながら時系列的に複数のエコー信号
を取得するので、これら複数のエコーを計測空間に配置
した場合、偶数エコーと奇数エコーではエコーのピーク
がそれぞれ異なった位置となる。[0004] One or several times of RF irradiation as an imaging method
There is a single-shot or multi-shot echo planar imaging (EPI) method for acquiring an echo signal necessary for reconstructing one image. In this imaging method, a plurality of echo signals are acquired in time series while reversing the readout gradient magnetic field. Therefore, when these plurality of echoes are arranged in the measurement space, the peaks of the echoes differ between the even-numbered echo and the odd-numbered echo. Position.
【0005】図9(a)はこの様子を模式的に示したもの
であり、例えば、奇数エコーの時系列データを図示する
計測空間901の左から順に配置するとすると、偶数エコ
ーは計測空間の右から配置することになる。従って奇数
エコーのピーク位置9021,9022・・・が計測空間中央904か
ら左側にずれた場合、偶数エコーのピーク位置9031,903
2・・・は計測空間中央904に対し右側にずれる。FIG. 9 (a) schematically shows this state. For example, if time series data of odd echoes are arranged in order from the left of the illustrated measurement space 901, even-numbered echoes are located on the right of the measurement space. Will be arranged from. Therefore, when the peak positions 9021, 9022... Of the odd echoes are shifted to the left from the center 904 of the measurement space, the peak positions 9031, 903 of the even echoes are shifted.
2 ... are shifted to the right with respect to the center 904 of the measurement space.
【0006】このようなデータを読み出し方向のフーリ
エ変換すると、偶数エコー、奇数エコーで位相の変化が
逆になるため、図9(b)のようなアーチファクト905が画
像領域906の両側に生じる。これは一般にN/2アーチファ
クトと呼ばれる。When such data is Fourier-transformed in the readout direction, the phase change is reversed between the even-numbered echo and the odd-numbered echo, so that an artifact 905 as shown in FIG. This is commonly referred to as N / 2 artifact.
【0007】このようなN/2アーチファクトを信号処理
によって除去する手法が提案されている(特開平8-21517
4号、特開平5-68674号など)。これらの方法では、被検
体の画像形成のための撮影(本計測)に先立って、位相エ
ンコード傾斜磁場を加えない状態で、あらかじめ静磁場
補正用のデータを取得しておき、この補正用データを用
いて計測したエコー信号(本計測データ)を補正する。A method for removing such N / 2 artifacts by signal processing has been proposed (JP-A-8-21517).
4, JP-A-5-68674). In these methods, prior to imaging (main measurement) for forming an image of a subject, static magnetic field correction data is acquired in advance without applying a phase encoding gradient magnetic field, and this correction data is obtained. The echo signal (main measurement data) measured by using this is corrected.
【0008】例えば、特開平8-215174号公報に記載され
た方法では、補正用データとして一組の偶数エコー、奇
数エコーを取得し、これらエコーをフーリエ変換した後
に、エコー間の位相差を求め、本計測データを補正す
る。また、特開平5-68674号公報に記載された方法で
は、補正用データ及び本計測データをそれぞれ読み出し
方向に1次元フーリエ変換し、本計測データから補正用
データの位相を減算することにより補正する。For example, in the method described in Japanese Patent Application Laid-Open No. 8-215174, a set of even-numbered echoes and odd-numbered echoes is acquired as correction data, and these echoes are subjected to Fourier transform, and then the phase difference between the echoes is determined. Then, the main measurement data is corrected. In the method described in Japanese Patent Application Laid-Open No. 5-68674, correction data and main measurement data are each subjected to one-dimensional Fourier transform in the readout direction, and correction is performed by subtracting the phase of the correction data from the main measurement data. .
【0009】[0009]
【発明が解決しようとする課題】しかし特開平8-215174
号公報に記載された補正方法では、計測空間における位
相エンコード方向の位相変化までは考慮されていなかっ
た。即ち、エコートレイン(1回の励起内で連続して計
測される複数のエコー)内で傾斜磁場の渦電流などによ
り位相差が異なる場合があるが、このような変化は上述
した補正方法では考慮されていない。特開平5-68674号
公報に記載された方法では、位相エンコード方向のオフ
セットが最適に調整されていない場合、図10(a)で示
すように補正用エコー信号を一次元フーリエ変換した補
正用データ1001,1002・・・の位相にノイズ的な変化1011,1
012が混入することがある。このような補正用データを
用いて補正した画像1020には、(b)に示すように線状の
アーチファクト(ストリークアーチファクト)1031,1032
が生じる。また、補正用データの位相に2次的な変化も
混入していることがあるため、補正後の画像が図中、10
41,1042で示すような歪みが生じる場合もある。SUMMARY OF THE INVENTION However, Japanese Patent Application Laid-Open No. 8-215174
In the correction method described in Japanese Patent Application Laid-Open Publication No. H10-265, the phase change in the phase encoding direction in the measurement space is not considered. In other words, the phase difference may differ due to the eddy current of the gradient magnetic field in the echo train (a plurality of echoes continuously measured in one excitation), and such a change is considered in the above-described correction method. It has not been. In the method described in JP-A-5-68674, when the offset in the phase encoding direction is not optimally adjusted, the correction data obtained by performing one-dimensional Fourier transform on the correction echo signal as shown in FIG. Noise-like change in the phase of 1001,1002 ...
012 may be mixed. The image 1020 corrected using such correction data has linear artifacts (streak artifacts) 1031 and 1032 as shown in (b).
Occurs. Also, since a secondary change may be mixed in the phase of the correction data, the corrected image is shown in FIG.
In some cases, distortion as shown by 41 and 1042 occurs.
【0010】そこで本発明は、位相エンコード方向のオ
フセットや2次的に混入位相ノイズも含めた補正が可能
であり、高精度でN/2アーチファクトを低減したMR画像
を得ることを目的とする。Accordingly, an object of the present invention is to obtain an MR image which can be corrected including an offset in a phase encoding direction and a secondary mixed phase noise, and has high accuracy and reduced N / 2 artifacts.
【0011】[0011]
【課題を解決するための手段】上記目的を達成するた
め、本発明では少なくとも1回の繰り返し時間内で計測
する本計測エコーに対応する数の前計測エコーを得て、
時間的に隣接する前計測エコーについてそれぞれ位相変
化を算出し、本計測エコーの位相エンコード方向の位相
変化を反映した補正用位相マップを作成する。また得ら
れた位相差(位相変化)について、その精度を向上する
処理を施す。このような補正用位相マップを用いて本計
測エコーを補正あるいはパルスシーケンスを調整するこ
とにより、N/2アーチファクトを高精度で除去すること
ができる。In order to achieve the above object, according to the present invention, a number of pre-measurement echoes corresponding to the main measurement echo measured within at least one repetition time is obtained.
The phase change is calculated for each of the pre-measurement echoes that are temporally adjacent to each other, and a correction phase map reflecting the phase change of the main measurement echo in the phase encoding direction is created. Further, the obtained phase difference (phase change) is subjected to a process for improving its accuracy. By correcting the main measurement echo or adjusting the pulse sequence using such a correction phase map, N / 2 artifacts can be removed with high accuracy.
【0012】即ち、本発明のMRI装置は、被検体に高周
波パルスを照射する手段と、読み出し傾斜磁場パルスを
印加する手段と、位相エンコード傾斜磁場パルスを印加
する手段と、少なくとも1回の繰り返し時間内で発生す
るエコー信号を時系列的に検出する手段と、エコー信号
から画像を作成し、表示する手段と、所定のパルスシー
ケンスに従い前記各手段を制御する制御手段とを有する
MRI装置において、前記制御手段は、画像取得のための
本計測エコー信号とは別に、位相エンコード傾斜磁場パ
ルスを印加せずに少なくとも1回の繰り返し時間内で計
測した前計測エコー信号を用いて、各前計測エコー信号
間の位相変化をそれぞれ算出し、それぞれの位相変化を
用いて補正用位相マップを作成するとともに、前記補正
用位相マップを用いて本計測エコー信号を取得するため
のパルスシーケンスのエコー計測時間を調整する手段を
有する。That is, the MRI apparatus of the present invention comprises a means for irradiating a subject with a high-frequency pulse, a means for applying a readout gradient magnetic field pulse, a means for applying a phase encoding gradient magnetic field pulse, and at least one repetition time. Means for detecting an echo signal generated in the apparatus in time series, means for creating and displaying an image from the echo signal, and control means for controlling the respective means according to a predetermined pulse sequence.
In the MRI apparatus, the control means, apart from the main measurement echo signal for image acquisition, using a pre-measurement echo signal measured within at least one repetition time without applying a phase encoding gradient magnetic field pulse, A phase change between each pre-measurement echo signal is calculated, and a phase map for correction is created using each phase change, and a pulse sequence for acquiring the main measurement echo signal using the phase map for correction is calculated. It has means for adjusting the echo measurement time.
【0013】また本発明のMRI装置は、前記制御手段
が、画像取得のための本計測エコー信号とは別に、位相
エンコード傾斜磁場パルスを印加せずに少なくとも1回
の繰り返し時間内で計測した前計測エコー信号を用い
て、各前計測エコー信号間の位相変化をそれぞれ算出
し、それぞれの位相変化を用いて補正用位相マップを作
成するとともに、少なくとも1回の繰り返し時間内で計
測した本計測エコー信号をそれぞれ前記補正用位相マッ
プの対応する補正データを用いて補正する手段を有す
る。In the MRI apparatus according to the present invention, the control means may perform a measurement before and after the measurement within at least one repetition time without applying a phase encoding gradient magnetic field pulse separately from the main measurement echo signal for acquiring an image. Using the measurement echo signal, calculate the phase change between each pre-measurement echo signal, create a correction phase map using each phase change, and measure the main measurement echo within at least one repetition time. Means for correcting the signals using the corresponding correction data of the correction phase map.
【0014】本発明のMRI装置では、連続して計測され
る複数のエコー信号のそれぞれについて位相変化を算出
し、補正用位相マップを作成するので、位相エンコード
方向に位相変化がある場合でも精度よく補正することが
できる。In the MRI apparatus of the present invention, a phase change is calculated for each of a plurality of continuously measured echo signals, and a correction phase map is created. Can be corrected.
【0015】また本発明のMRI装置は、均一磁場空間を
形成するための静磁場を発生する静磁場発生手段と、均
一磁場空間に置かれた被検体に対し高周波パルスを照射
する高周波パルス発生手段と、前記被検体に対しスライ
ス方向、位相エンコード方向および読み出し方向の各傾
斜磁場を印加する傾斜磁場発生手段と、前記被検体から
のエコー信号を検出する検出手段と、前記エコー信号か
ら画像を再構成する画像処理手段と、得られた画像を表
示する画像表示手段と、所定のシーケンスに従い高周波
パルス発生手段、傾斜磁場発生手段の発生タイミングを
制御するとともに検出手段、画像処理手段、画像表示手
段をそれぞれ制御する制御手段とを備えたMRI装置にお
いて、前記制御手段は、位相エンコードを印加しない前
計測シーケンスの後に位相エンコードを付加した本計測
シーケンスを実行するよう制御し、前記シーケンスにお
いて、少なくとも1回の繰り返し時間内で連続的に計測
した前計測エコー信号を用いて、各前計測エコー信号間
の位相変化をそれぞれ算出するとともに、計測空間上に
時系列に配置される本計測エコー信号を対応する前記位
相変化の情報に基づいて、偶数番目取得の本計測エコー
信号の位相を奇数番目取得の本計測エコー信号の位相に
揃える、奇数番目取得の本計測エコー信号の位相を偶数
番目取得の本計測エコー信号の位相に揃える、もしくは
計測中心の位相に揃える、のいずれかの処理を実行する
位相補正手段を備えたことを特徴とする。Further, the MRI apparatus of the present invention comprises a static magnetic field generating means for generating a static magnetic field for forming a uniform magnetic field space, and a high frequency pulse generating means for irradiating a subject placed in the uniform magnetic field space with a high frequency pulse. Gradient magnetic field generating means for applying a gradient magnetic field in the slice direction, phase encoding direction, and readout direction to the subject; detecting means for detecting an echo signal from the subject; The image processing means to be constituted, the image display means for displaying the obtained image, the high-frequency pulse generation means according to a predetermined sequence, controlling the generation timing of the gradient magnetic field generation means and the detection means, the image processing means, the image display means In an MRI apparatus having control means for controlling each, the control means is provided after a pre-measurement sequence in which phase encoding is not applied. The main measurement sequence to which phase encoding is added is controlled to be executed, and in the sequence, a phase change between each pre-measurement echo signal is measured using pre-measurement echo signals continuously measured within at least one repetition time. Based on the phase change information corresponding to the main measurement echo signals arranged in time series on the measurement space, the phase of the even measurement acquisition main measurement echo signal is changed to the odd measurement acquisition measurement signal. The phase of the main measurement echo signal of the odd-number acquisition is adjusted to the phase of the main measurement echo signal of the even-number acquisition, or the phase of the measurement center is adjusted. It is characterized by having.
【0016】本発明の好適な態様のMRI装置は、補正用
位相マップの作成において、算出した位相変化の精度を
向上する処理を含む。これら処理は、具体的には、主値
回り補正、位相差の平均処理を含む。The MRI apparatus according to a preferred embodiment of the present invention includes a process for improving the accuracy of the calculated phase change in creating a correction phase map. More specifically, these processes include main value rotation correction and phase difference averaging.
【0017】[0017]
【発明の実施の形態】以下、本発明のMRI装置につい
て、図面を参照して詳述する。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, an MRI apparatus according to the present invention will be described in detail with reference to the drawings.
【0018】図8は本発明が適用されるMRI装置の全体
構成を示す図で、このMRI装置は、被検体801を囲む所定
の空間に静磁場を発生する磁石802と、該空間に傾斜磁
場を発生する傾斜磁場コイル803と、この空間に高周波
磁場を発生するRFコイル804と、被検体801が発生するMR
信号を検出するRFプローブ805とを備えている。ベッド8
12は被検体が横たわるためのものである。FIG. 8 is a diagram showing the overall configuration of an MRI apparatus to which the present invention is applied. This MRI apparatus includes a magnet 802 for generating a static magnetic field in a predetermined space surrounding a subject 801 and a gradient magnetic field in the space. , A RF coil 804 for generating a high-frequency magnetic field in this space, and an MR coil for generating a subject 801.
An RF probe 805 for detecting a signal is provided. Bed 8
Numeral 12 is for the subject to lie down.
【0019】傾斜磁場コイル803は、X、Y、Zの3方向の
傾斜磁場コイルで構成され、傾斜磁場電源809からの信
号に応じてそれぞれ傾斜磁場を発生する。傾斜磁場に
は、被検体の所定のスライスを選択するために印加され
るスライス選択傾斜磁場、エコー信号を位相エンコード
するために印加される位相エンコード傾斜磁場およびエ
コー信号に周波数エンコードするとともにグラディエン
トエコーを発生させるための読み出し傾斜磁場があり、
本発明においては極性の反転する読み出し傾斜磁場を使
用する。The gradient magnetic field coil 803 is composed of gradient magnetic field coils in three directions of X, Y and Z, and generates a gradient magnetic field in accordance with a signal from a gradient magnetic field power supply 809. The gradient magnetic field includes a slice selection gradient magnetic field applied to select a predetermined slice of the subject, a phase encoding gradient magnetic field applied to phase-encode the echo signal, and frequency encoding of the echo signal and gradient echo. There is a read-out gradient magnetic field to generate
In the present invention, a read gradient magnetic field whose polarity is inverted is used.
【0020】RFコイル804はRF送信部810の信号に応じて
高周波磁場を発生する。RFプローブ805の信号は、信号
検出部806で検出され、信号処理部807で信号処理され、
また計算により画像信号に変換される。画像は表示部80
8で表示される。The RF coil 804 generates a high-frequency magnetic field according to a signal from the RF transmitting unit 810. The signal of the RF probe 805 is detected by the signal detection unit 806, and the signal is processed by the signal processing unit 807.
It is converted into an image signal by calculation. The image is displayed on the display unit 80
Displayed as 8.
【0021】傾斜磁場電源809、RF送信部810、信号検出
部806は、一般にパルスシーケンスと呼ばれるタイムチ
ャートに従い制御部811で制御される。The gradient magnetic field power supply 809, the RF transmitter 810, and the signal detector 806 are controlled by a controller 811 according to a time chart generally called a pulse sequence.
【0022】制御部811は、上述した信号処理部807で信
号処理されたデジタル信号であるエコー信号に各種演算
を行うとともに、反転する読み出し傾斜磁場の各周期内
で発生するエコー信号を用いて位相変化の算出、補正用
位相マップの作成、補正用位相マップを用いたパルスシ
ーケンスの調整および/または補正用位相マップを用い
た本計測エコーの補正などを行う機能を備えている。こ
れら制御部811の各機能は、制御部811を構成するコンピ
ュータのプログラムとして構築される。The control unit 811 performs various operations on the echo signal, which is a digital signal processed by the signal processing unit 807, and uses the echo signal generated in each cycle of the inverted read-out gradient magnetic field to perform phase calculation. It has a function of calculating a change, creating a correction phase map, adjusting a pulse sequence using the correction phase map, and / or correcting a main measurement echo using the correction phase map. Each function of the control unit 811 is constructed as a computer program constituting the control unit 811.
【0023】次に、このような構成における本発明の第
1の実施形態を説明する。この実施形態では、図1に示
すように本計測200に先だって前計測シーケンス100を実
行し、この前計測シーケンスで得た補正用エコー101を
用いて本計測シーケンスのサンプリング時間を調整す
る。Next, a first embodiment of the present invention having such a configuration will be described. In this embodiment, as shown in FIG. 1, the pre-measurement sequence 100 is executed prior to the main measurement 200, and the sampling time of the main measurement sequence is adjusted using the correction echo 101 obtained in the pre-measurement sequence.
【0024】本計測シーケンスは、シングルショット或
いはマルチショットのEPIシーケンスであり、高周波パ
ルスにより被検体の所定領域を励起した後、連続して反
転する読み出し傾斜磁場を印加しながら、これと同時に
位相エンコード傾斜磁場を印加して、位相エンコードの
異なる複数のエコ―信号を計測する。位相エンコードの
数は通常1枚の画像あたり64、128、256、512等の値が
選ばれる。シングルショットEPIでは、これら位相エン
コードの異なる信号を一度の励起で計測し、マルチショ
ットEPIでは分割して計測する。各エコー信号は通常12
8、256、512、1024個のサンプリングデータからなる時
系列信号として得られる。This measurement sequence is a single-shot or multi-shot EPI sequence in which a predetermined area of the subject is excited by a high-frequency pulse, and then a read-out gradient magnetic field that is continuously inverted is applied while simultaneously performing phase encoding. A gradient magnetic field is applied to measure a plurality of echo signals having different phase encodings. As the number of phase encodings, values such as 64, 128, 256, and 512 are usually selected for one image. In the single shot EPI, these signals having different phase encodings are measured by one excitation, and in the multi shot EPI, the signals are divided and measured. Each echo signal is usually 12
It is obtained as a time-series signal composed of 8, 256, 512, and 1024 sampling data.
【0025】以下、本計測シーケンスがマルチショット
EPIである場合を例にして説明する。図2は補正前の本
計測シーケンスを示す図で、このパルスシーケンスでは
高周波パルス201とスライス選択傾斜磁場パルス202を印
加した後、位相エンコードのオフセットを与えるパルス
203と読み出し傾斜磁場のオフセットを与えるパルス204
を印加し、位相エンコード傾斜磁場パルス205を離散的
に印加しながら反転する読み出し傾斜磁場206の各周期
内で各位相エンコードのエコー信号102が発生するの
で、これをサンプリング時間207の間おのおのサンプリ
ングし時系列データ102を得る。サンプリング時間207は
典型的にはそれぞれ1ms程度である。In the following, this measurement sequence is a multi-shot
The case of EPI will be described as an example. FIG. 2 is a diagram showing a main measurement sequence before correction. In this pulse sequence, a pulse for applying a phase encoding offset after applying a high-frequency pulse 201 and a slice selection gradient magnetic field pulse 202.
203 and pulse 204 giving read gradient magnetic field offset
The echo signal 102 of each phase encoding is generated within each period of the readout gradient magnetic field 206 which is inverted while discretely applying the phase encoding gradient magnetic field pulse 205, and this is sampled for each sampling time 207. Time series data 102 is obtained. The sampling time 207 is typically about 1 ms each.
【0026】このようなシーケンス208を複数回(N回)繰
り返し、画像再構成に必要な全ての本計測エコ-信号h
(n,m,t)102を取得する。なお、nは繰り返し番号(1≦n≦
N)、mはエコー番号(1≦m≦M)、t は時間を表す。一例と
して、1回の高周波パルス印加で計測する本計測エコー
信号の数が16(M=16)であり、読み出し方向のデータ数が
128(x=128)、位相エンコード方向のエンコード数が128
であるとすると、1枚の画像を得るのに必要な繰り返し
数は8(N=8)となる。This sequence 208 is repeated a plurality of times (N times), and all the actual measurement echo signals h required for image reconstruction are obtained.
(n, m, t) 102 is acquired. Note that n is a repetition number (1 ≦ n ≦
N) and m are echo numbers (1 ≦ m ≦ M), and t is time. As an example, the number of main measurement echo signals measured by one high-frequency pulse application is 16 (M = 16), and the number of data in the readout direction is
128 (x = 128), the number of encodes in the phase encode direction is 128
, The number of repetitions necessary to obtain one image is 8 (N = 8).
【0027】既に述べたようにこのように計測された本
計測エコーは、エコーのピーク位置がサンプリング時間
207の中心と一致しないために、計測空間において奇数
番目(mが奇数)のエコーのピーク位置と、偶数番目
(mが偶数)のエコーのピーク位置が計測空間中央に対
し、それぞれ逆方向にずれる(図9(a))。As described above, the main measurement echo measured in this manner has a peak position of the echo corresponding to the sampling time.
Since it does not coincide with the center of 207, the peak position of the odd-numbered echo (m is an odd number) and the peak position of the even-numbered echo (m is an even number) in the measurement space are respectively shifted in opposite directions with respect to the center of the measurement space. (FIG. 9 (a)).
【0028】この計測空間におけるエコー信号のピーク
ずれは、フーリエ変換後の空間での位相変化に対応す
る。そこで、まず前計測によって奇数エコーと偶数エコ
ーの位相変化を求め、この位相変化からピーク位置のず
れを求め、奇数エコーおよび偶数エコーのピーク位置を
一致するように本計測シーケンスにおけるサンプリング
時間を調整する。The peak shift of the echo signal in the measurement space corresponds to a phase change in the space after the Fourier transform. Therefore, first, the phase change between the odd-numbered echo and the even-numbered echo is obtained by the previous measurement, the shift of the peak position is obtained from this phase change, and the sampling time in the main measurement sequence is adjusted so that the peak positions of the odd-numbered echo and the even-numbered echo coincide. .
【0029】図3はこのような位相検出ための前計測パ
ルスシーケンスの一例を示す図であり、本計測シーケン
スに対応してマルチショットEPIシーケンスを基本とし
ている。但し、図3のパルスシーケンスは位相エンコー
ド傾斜磁場を加えない点が本計測のEPIシーケンスとは
異なる。FIG. 3 is a diagram showing an example of a pre-measurement pulse sequence for detecting such a phase, which is based on a multi-shot EPI sequence corresponding to the main measurement sequence. However, the pulse sequence of FIG. 3 is different from the EPI sequence of this measurement in that the phase encoding gradient magnetic field is not applied.
【0030】即ち、まず検知する磁化を含む被検体に高
周波パルス201を照射すると同時にスライスを選択する
傾斜磁場パルス202を印加し、画像化するスライスを選
択する。次いで読み出し傾斜磁場のオフセットを与える
パルス204を印加した後、連続して反転する読み出し傾
斜磁場パルス206を印加し、反転する読み出し傾斜磁場2
06の各周期内でエコー信号101が時系列的に発生するの
で、これを時間範囲207の間おのおのサンプリングし時
系列データ101を得る。反転する読み出し傾斜磁場206の
強度、印加タイミング、時間範囲207は本計測シーケン
スと同じである。That is, first, a high-frequency pulse 201 is applied to a subject including a magnetization to be detected, and at the same time, a gradient magnetic field pulse 202 for selecting a slice is applied to select a slice to be imaged. Next, after applying a pulse 204 for giving an offset of the read gradient magnetic field, a read gradient magnetic field pulse 206 for continuously inverting is applied, and a read gradient magnetic field
Since the echo signal 101 is generated in time series in each cycle of 06, the echo signal 101 is sampled for each time range 207 to obtain the time series data 101. The intensity, application timing, and time range 207 of the read gradient magnetic field 206 to be reversed are the same as those in the main measurement sequence.
【0031】このようにシーケンス208内で取得された
時系列データ101、即ち前計測エコ-信号p(m,t)101を用
いて位相差マップを作成する。なお、mはエコー番号(1
≦m≦M)、tは時間を表す。このため、まず前計測エコー
信号p(m,t)101を読み出し方向にフーリエ変換し、補正
用データp(m,x)105を取得する(ステップ103)。ここ
で、xは読み出し方向の位置を表し、1≦x≦Xである。A phase difference map is created using the time series data 101 acquired in the sequence 208, that is, the previously measured echo signal p (m, t) 101. M is the echo number (1
≦ m ≦ M), and t represents time. For this reason, first, the pre-measurement echo signal p (m, t) 101 is Fourier-transformed in the readout direction to obtain correction data p (m, x) 105 (step 103). Here, x represents a position in the reading direction, and 1 ≦ x ≦ X.
【0032】次いで偶数エコーの補正用データと奇数エ
コーの補正用データとから両者間の位相差s(k,x)108 を
計算する(ステップ107)。具体的には、下記計算を行
う。Next, a phase difference s (k, x) 108 between the correction data for the even-numbered echo and the correction data for the odd-numbered echo is calculated (step 107). Specifically, the following calculation is performed.
【0033】[0033]
【数1】 ここで、kは1≦k≦M/2であり、re[]、im[]はそれぞれ信
号の実部と虚部を表し、||は信号の絶対値を表す。こ
の計算では、奇数エコー(2k-1番目のエコー)とそれと
時系列的に隣接する2つの偶数エコー(2k-2番目と2k番
目)との位相差をそれぞれ求め、その平均をとり、位相
差s(k,x)としている。(Equation 1) Here, k is 1 ≦ k ≦ M / 2, re [] and im [] represent the real and imaginary parts of the signal, respectively, and || represents the absolute value of the signal. In this calculation, the phase difference between the odd-numbered echo (2k-1th echo) and two even-numbered echoes (2k-2nd and 2kth) that are adjacent in time series is calculated, and the average is calculated. s (k, x).
【0034】次に位相差s(k,x)の精度を向上するための
処理(109,110)を行う。上式からもわかるようにこの
計算は、-π〜+πの範囲でなされるため、-πを超える
ものは折り返して+πの方へ、+πを超えるものは折り近
して-πの方へそれぞれ値が離散的に変化する(これを
主値回りといい、局所的に位相が主値回りするものを主
値飛びという)。上述のように求めた位相差s(k,x)に
は、この主値飛びノイズが混入しているため、これを下
記計算により除去をする(ステップ109)。Next, processing (109, 110) for improving the accuracy of the phase difference s (k, x) is performed. As can be seen from the above equation, since this calculation is performed in the range of -π to + π, those exceeding -π are folded back toward + π, and those exceeding + π are folded close to -π The values each change discretely (this is called around the main value, and when the phase locally turns around the main value, it is called the main value jump). Since the main value jump noise is mixed in the phase difference s (k, x) obtained as described above, it is removed by the following calculation (step 109).
【0035】[0035]
【数2】 ここで、i、j はフィルタ範囲の点を表し、それぞれ-X/
2≦i≦X/2、-X/2≦j≦X/2、である。(Equation 2) Here, i and j represent points in the filter range, and -X /
2 ≦ i ≦ X / 2 and −X / 2 ≦ j ≦ X / 2.
【0036】さらに主値飛びノイズ以外に位相差に混入
するノイズを除去するため、下式により位相差を平均す
る(ステップ110)。Further, in order to remove noise mixed into the phase difference other than the main value jump noise, the phase difference is averaged by the following equation (step 110).
【0037】[0037]
【数3】 ここで、i、jはフィルタ範囲の点を表し、それぞれX/2
≦i≦X/2、X/2≦j≦X/2、である。(Equation 3) Here, i and j represent points in the filter range, and X / 2
≦ i ≦ X / 2, X / 2 ≦ j ≦ X / 2.
【0038】これらの処理109、110によりピークずれの
検出の精度を向上することができる。The accuracy of peak shift detection can be improved by these processes 109 and 110.
【0039】次にこのように求めた位相差マップを用い
て本計測シーケンスにおけるサンプリング時間207を設
定する。図4は、エコー信号のピークの計測空間中央か
らのずれと位相変化との関係を示したものであり、図示
するように、エコー信号402のピークが計測空間の中央4
01にあるときには、エコーを一次元フーリエ変換した空
間での位相の変化は、直線405で示すように0であるが、
エコー信号403のピークが中央401からずれたときには
(ずれ量404)、このエコー信号を一次元フーリエ変換し
た空間では、位相は一次関数的に直線406のように変化
する。従って位相変化の傾きを求めることによって、ず
れ量(時間)を求めることができる。Next, the sampling time 207 in the main measurement sequence is set using the phase difference map obtained in this manner. FIG. 4 shows the relationship between the shift of the peak of the echo signal from the center of the measurement space and the phase change. As shown in FIG.
When it is at 01, the change in phase in the space obtained by one-dimensional Fourier transform of the echo is 0 as shown by a straight line 405,
When the peak of the echo signal 403 deviates from the center 401
(Shift amount 404), in a space obtained by one-dimensional Fourier transform of the echo signal, the phase changes like a straight line 406 in a linear function. Therefore, the shift amount (time) can be obtained by obtaining the inclination of the phase change.
【0040】このため、ステップ110により作成した位
相差s'(k,x)をフィッティング処理して位相変化の傾き
(1次成分)を計算する(ステップ111)。フィッティン
グとして、例えば最小2乗法を採用できる。このような
フィッティングを位相差s'(k,x)の全てのkについて行
う。For this reason, the phase difference s ′ (k, x) created in step 110 is subjected to fitting processing to calculate the slope of the phase change.
(First order component) is calculated (step 111). As the fitting, for example, the least square method can be adopted. Such fitting is performed for all k of the phase difference s ′ (k, x).
【0041】ステップ111でs'(k,x)の各kについて求め
た傾きa(k)から次式によりサンプリング時間207のずれ
Δt(k)112を計算する。In step 111, a deviation Δt (k) 112 of the sampling time 207 is calculated from the gradient a (k) obtained for each k of s ′ (k, x) by the following equation.
【0042】[0042]
【数4】 ここで、Spはサンプリング207のピッチ(1点のサンプリ
ング時間)、Atはサンプリング207の時間である。1ピク
セルの位相差は2πであるから、この計算により、サブ
ピクセル単位で高精度にピークずれを検出できることが
わかる。例えば、0.1ピクセル単位のピーク差が、エコ
ー毎に分かる。(Equation 4) Here, Sp is the pitch of sampling 207 (sampling time of one point), and At is the time of sampling 207. Since the phase difference of one pixel is 2π, it can be seen from this calculation that the peak shift can be detected with high accuracy in subpixel units. For example, a peak difference of 0.1 pixel unit is known for each echo.
【0043】次に、求めたΔt(k)から、本計測シーケン
スのエコー計測時間を計算する(ステップ113)。具体
的には、図2の本計測シーケンスにおいて奇数番めのエ
コーのサンプリング時間はそのままにして、偶数番目の
エコーのサンプリング開始時間をΔt(k)ずらす。このよ
うに調整された本計測シーケンスを図5に示す。ここで
は本計測シーケンス中、読み出し傾斜磁場(Gr)、サンプ
リング(A/D)、発生したエコー信号(echo)についてのみ
示した。Next, the echo measurement time of the main measurement sequence is calculated from the obtained Δt (k) (step 113). Specifically, the sampling time of the even-numbered echo is shifted by Δt (k) while the sampling time of the odd-numbered echo is kept as it is in the main measurement sequence of FIG. FIG. 5 shows the main measurement sequence adjusted in this manner. Here, only the readout gradient magnetic field (Gr), sampling (A / D), and the generated echo signal (echo) are shown during the main measurement sequence.
【0044】図示するように、エコー信号102はそのピ
ーク位置が計測空間中心から 501(5011,5012・・・)だけ
ずれている。このずれ501は、上述の計算によって求め
たΔt(k)/2に対応する。そこで、偶数番目のエコーのサ
ンプリング時間207を、Δt(k)(502)だけずらす。これ
により計測空間のエコー信号のピーク位置がそろうこと
になる。As shown, the peak position of the echo signal 102 is shifted from the center of the measurement space by 501 (5011, 5012...). This shift 501 corresponds to Δt (k) / 2 obtained by the above calculation. Therefore, the sampling time 207 of the even-numbered echo is shifted by Δt (k) (502). As a result, the peak positions of the echo signals in the measurement space are aligned.
【0045】このようにステップ103から113までの前計
測処理114で計算した本計測シーケンスを実行し、本計
測エコー信号102 を取得する。図6に取得した本計測エ
コー102の計測空間601の配列を示す。図示するように、
偶数エコー602のピーク6021,6022・・・は、サンプリング
時間をずらすことにより、奇数エコー603のピーク6031,
6032…と同じ位置になる。これらは全体として計測空間
601の中心604からずれているが、各エコーのピークがそ
ろっているため、これを2次元フーリエ変換して画像再
構成119を行い作成した画像118ではN/2アーチファクト
が大幅に低減される。As described above, the main measurement sequence calculated in the pre-measurement processing 114 from steps 103 to 113 is executed, and the main measurement echo signal 102 is obtained. FIG. 6 shows an arrangement of the measurement space 601 of the main measurement echo 102 acquired. As shown
The peaks 6021, 6022,... Of the even-numbered echoes 602 are changed by shifting the sampling time.
The same position as 6032 ... These are the measurement space as a whole
Although it is shifted from the center 604 of the 601, the peaks of the respective echoes are aligned. Therefore, in the image 118 created by performing the two-dimensional Fourier transform and performing the image reconstruction 119, the N / 2 artifact is greatly reduced.
【0046】尚、図6では、計測空間の中心からのずれ
がエコー番号に関わらず同一である場合を示している
が、上述したように本発明では位相差はエコー番号
(k)の関数として求められるので、エコートレイン内
で位相差が変化する場合でも正確にピークを揃えること
ができる。FIG. 6 shows a case where the deviation from the center of the measurement space is the same regardless of the echo number. However, as described above, in the present invention, the phase difference is a function of the echo number (k). Since it is required, even if the phase difference changes within the echo train, the peaks can be accurately aligned.
【0047】以上、本発明の第1の実施形態を説明し
た。この実施形態では偶数エコーのサンプリング時間を
調整する場合を説明したが、奇数エコーを調整しても良
いし、両者を調整することも可能である。またシーケン
スについてもマルチショットのみならず、ワンショット
型のシーケンスや、マルチスライス、3D計測等に本発明
の処理を適用することも可能である。As described above, the first embodiment of the present invention has been described. In this embodiment, the case where the sampling time of the even-numbered echo is adjusted has been described. However, the odd-numbered echo may be adjusted, or both may be adjusted. The processing of the present invention can be applied not only to multi-shot sequences but also to one-shot type sequences, multi-slices, 3D measurement, and the like.
【0048】第1の実施形態によるMRI装置によれば、補
正用データからエコー間の位相変化を算出し本計測シー
ケンスをエコー毎にサブピクセル単位で最適に調整する
ので、画像歪みや線状のアーチファクトを生じることな
く、画質の良い撮影を実現できる。According to the MRI apparatus of the first embodiment, the phase change between echoes is calculated from the correction data, and the main measurement sequence is optimally adjusted in subpixel units for each echo. High quality image capturing can be realized without causing artifacts.
【0049】また補正用データを用いて画像再構成時に
補正処理を行うのではなく、本計測シーケンスを調整し
ているので、再構成時間を延長することなく本計測を行
うことができる。Further, since the main measurement sequence is adjusted without performing the correction processing at the time of image reconstruction using the correction data, the main measurement can be performed without extending the reconstruction time.
【0050】次に、本発明の第2の実施形態を説明す
る。この実施形態でも、本計測に先だって前計測シーケ
ンスを実行し、エコー毎の補正用データを作成すること
は第1の実施形態と同様であるが、第2の実施形態で
は、この補正用データp(m,x)101を用いて本計測シーケ
ンスで計測した各エコーの位相を補正する。Next, a second embodiment of the present invention will be described. In this embodiment as well, the pre-measurement sequence is executed prior to the main measurement to create correction data for each echo as in the first embodiment, but in the second embodiment, the correction data p The phase of each echo measured in the main measurement sequence is corrected using (m, x) 101.
【0051】図7にこの実施形態による処理の手順を示
す。図1と同一の処理については図1と同じ符号で示し
た。ここでもマルチショットEPIを例にして説明する。
まず図2に示す前計測シーケンスを実行し、補正用エコ
ーp(m,t)を計測する(ステップ101)。ここでmはエコー
番号、tは時間であり、本計測シーケンスの1繰り返し
時間で計測するエコーと同じ数のエコーを計測する。こ
の補正用エコーをエコー毎に読み出し方向にフーリエ変
換し、補正用データp(m,x)を得る(ステップ103)。FIG. 7 shows a processing procedure according to this embodiment. The same processes as those in FIG. 1 are denoted by the same reference numerals as those in FIG. Here, the multi-shot EPI will be described as an example.
First, the pre-measurement sequence shown in FIG. 2 is executed to measure the correction echo p (m, t) (step 101). Here, m is an echo number, and t is time, and the same number of echoes as those measured in one repetition time of the main measurement sequence are measured. This correction echo is Fourier-transformed in the readout direction for each echo to obtain correction data p (m, x) (step 103).
【0052】この補正用データから奇数エコーと偶数エ
コーの間の位相差s(k,x)を計算する(ステップ107)。
この計算は、第1の実施形態において位相差s(k,x)を求
めた計算と同じである。従って、この場合にも主値飛び
ノイズの除去(ステップ109)と位相差の平均を求める
処理(ステップ110)を行う。このように精度を向上し
た位相差s'(k,x)は、補正用位相マップとして制御部の
メモリに格納される。The phase difference s (k, x) between the odd echo and the even echo is calculated from the correction data (step 107).
This calculation is the same as the calculation for obtaining the phase difference s (k, x) in the first embodiment. Therefore, also in this case, the process of removing the main value jump noise (step 109) and the process of obtaining the average of the phase difference (step 110) are performed. The phase difference s ′ (k, x) with improved accuracy is stored in the memory of the control unit as a correction phase map.
【0053】一方、図2に示す本計測シーケンスを実行
し、本計測エコーh(n,m,t)を計測する(ステップ10
2)。ここでnは繰り返し番号、mはエコー番号、tは時間
である。次いで本計測エコーh(n,m,t)を読み出し方向に
フーリエ変換し、本計測データh(n,m,x)を得る(ステッ
プ106)。On the other hand, the main measurement sequence shown in FIG. 2 is executed, and the main measurement echo h (n, m, t) is measured (step 10).
2). Here, n is a repetition number, m is an echo number, and t is time. Next, the main measurement echo h (n, m, t) is Fourier-transformed in the readout direction to obtain main measurement data h (n, m, x) (step 106).
【0054】この本計測データを前計測で求めた補正用
位相マップs'(k,x)を用いて補正する(ステップ111)。
補正の計算は、例えば次式により行う。The main measurement data is corrected using the correction phase map s' (k, x) obtained in the previous measurement (step 111).
The calculation of the correction is performed by, for example, the following equation.
【0055】[0055]
【数5】 上式からわかるように、この例では奇数エコーは、計測
したデータをそのまま用い、偶数エコーのみを補正して
いる。これによって図6に示したのと同様に計測空間に
おいて偶数エコー602のピーク6011,6012・・・は、奇数エ
コー603のピーク6031,6032…と同じ位置になる。この場
合、奇数エコーを補正し、偶数エコーに揃えるようにし
てもよいし、両者を補正し、そのピークが計測空間中心
になるように揃えてもよい。(Equation 5) As can be seen from the above equation, in this example, the odd-numbered echo uses the measured data as it is, and corrects only the even-numbered echo. .., The peaks 6011, 6012,... Of the even-numbered echoes 602 are in the same positions as the peaks 6031, 6032,. In this case, odd-numbered echoes may be corrected and aligned with even-numbered echoes, or both may be corrected and their peaks may be aligned with the center of the measurement space.
【0056】このように奇数エコーと偶数エコーの位相
を揃えた後、本計測データを位相エンコード方向にフー
リエ変換することにより(ステップ112)画像再構成さ
れる。このように再構成した画像113ではN/2アーチファ
クトが大幅に低減される。After the phases of the odd-numbered echo and the even-numbered echo are aligned in this manner, the image is reconstructed by subjecting the main measurement data to Fourier transform in the phase encoding direction (step 112). In the reconstructed image 113, N / 2 artifacts are significantly reduced.
【0057】この実施形態でも位相差はエコー番号
(k)の関数として求められ、エコー番号mが異なる本
計測エコーについてそれぞれ対応する位相差データを用
いて補正するので、位相エンコード方向に位相差の変動
がある場合にも、即ち、エコートレイン内で位相差が変
化する場合でも、正確に補正することができる。これに
より画像歪みや線状のアーチファクトを生じることな
く、画質の良い撮影を実現することができる。Also in this embodiment, the phase difference is obtained as a function of the echo number (k), and the main measurement echoes having different echo numbers m are corrected using the corresponding phase difference data. Even when there is a fluctuation, that is, even when the phase difference changes within the echo train, it is possible to correct accurately. As a result, it is possible to realize high quality image capturing without causing image distortion and linear artifacts.
【0058】以上、本発明の第2の実施形態を説明した
が、この実施形態も第1の実施形態と同様にワンショッ
トシーケンスやマルチスライス計測、3D計測にも適用す
ることができる。Although the second embodiment of the present invention has been described above, this embodiment can also be applied to one-shot sequence, multi-slice measurement, and 3D measurement as in the first embodiment.
【0059】[0059]
【発明の効果】本発明は以上のように構成されたので、
下記の効果が得られる。 全エコーについて、偶数エコー、奇数エコー信号の位
相差をそれぞれ計算するので、 エコートレイン内で傾
斜磁場の渦電流などにより、上記差がエコー毎に異なっ
ていても、正確に本計測エコーを補正できる。或いは本
計測シーケンスを調整できる。 読み出し方向フーリエ変換後の補正用データについ
て、精度を向上する処理を行うので、従来の補正では除
去できなかった脳低部などのストリークアーチファクト
を除去することができ、画質のよいMR画像を得ることが
できる。The present invention has been configured as described above.
The following effects are obtained. For all echoes, the phase difference between the even-numbered echo and the odd-numbered echo signal is calculated, so even if the difference is different for each echo due to the eddy current of the gradient magnetic field in the echo train, the actual measurement echo can be accurately corrected. . Alternatively, the main measurement sequence can be adjusted. Performs processing to improve the accuracy of the correction data after the Fourier transform in the readout direction, so that streak artifacts such as the lower brain that could not be removed by conventional correction can be removed, and MR images with good image quality can be obtained. Can be.
【図1】本発明のMRI装置の第1の実施形態における処理
手順を示す図。FIG. 1 is a diagram showing a processing procedure in a first embodiment of the MRI apparatus of the present invention.
【図2】本発明のMRI装置で実行される本計測パルスシ
ーケンスの一例を示す図。FIG. 2 is a diagram showing an example of a main measurement pulse sequence executed by the MRI apparatus of the present invention.
【図3】本発明のMRI装置で実行される前計測のパルス
シーケンスの一例を示す図。FIG. 3 is a diagram showing an example of a pulse sequence of pre-measurement executed by the MRI apparatus of the present invention.
【図4】エコーのピーク位置のずれと位相変化の関係を
説明する図。FIG. 4 is a view for explaining a relationship between a shift of a peak position of an echo and a phase change.
【図5】第1の実施形態のMRI装置で実行される、調整後
の本計測パルスシーケンスの一例を示す図。FIG. 5 is a view showing an example of a main measurement pulse sequence after adjustment, which is executed by the MRI apparatus of the first embodiment.
【図6】図5のパルスシーケンスで計測したエコーの計
測空間配置を示す図。FIG. 6 is a diagram showing a measurement space arrangement of echoes measured by the pulse sequence of FIG. 5;
【図7】本発明のMRI装置の第2の実施形態における処
理手順を示す図。FIG. 7 is a diagram showing a processing procedure in a second embodiment of the MRI apparatus of the present invention.
【図8】本発明が適用されるMRI装置の全体構成を示す
図。FIG. 8 is a diagram showing an overall configuration of an MRI apparatus to which the present invention is applied.
【図9】N/2アーチファクトの発生原理を説明する図。FIG. 9 is a view for explaining the principle of generation of N / 2 artifacts.
【図10】その他のアーチファクトを説明する図。FIG. 10 is a view for explaining other artifacts.
803・・・読み出し傾斜磁場を印加する手段、位相エンコー
ド傾斜磁場を印加する手段 804・・・エコー信号を検出する手段 805・・・高周波パルスを照射する手段 811・・・制御手段803 means for applying a read gradient magnetic field, means for applying a phase encoding gradient magnetic field 804 means for detecting an echo signal 805 means for applying a high-frequency pulse 811 means for controlling
Claims (3)
読み出し傾斜磁場パルスを印加する手段と、位相エンコ
ード傾斜磁場パルスを印加する手段と、少なくとも1回
の繰り返し時間内で発生するエコー信号を時系列的に検
出する手段と、エコー信号から画像を作成し、表示する
手段と、所定のパルスシーケンスに従い前記各手段を制
御する制御手段とを有する磁気共鳴イメージング装置に
おいて、 前記制御手段は、画像取得のための本計測エコー信号と
は別に、位相エンコード傾斜磁場パルスを印加せずに少
なくとも1回の繰り返し時間内で計測した前計測エコー
信号を用いて、各前計測エコー信号間の位相変化をそれ
ぞれ算出し、それぞれの位相変化を用いて補正用位相マ
ップを作成するとともに、前記補正用位相マップを用い
て本計測エコー信号を取得するためのパルスシーケンス
のエコー計測時間を調整する手段を備えたことを特徴と
する磁気共鳴イメージング装置。A means for irradiating a subject with a high-frequency pulse;
Means for applying a read gradient magnetic field pulse, means for applying a phase encoding gradient magnetic field pulse, means for detecting an echo signal generated within at least one repetition time in time series, and creating an image from the echo signal. A magnetic resonance imaging apparatus comprising: a display unit; and a control unit that controls each of the units in accordance with a predetermined pulse sequence. The control unit, apart from the main measurement echo signal for acquiring an image, Using a pre-measurement echo signal measured within at least one repetition time without applying a pulse, a phase change between each pre-measurement echo signal is calculated, and a phase map for correction is calculated using each phase change. A pulse sequence echometer for creating and obtaining a main measurement echo signal using the correction phase map Magnetic resonance imaging apparatus characterized by comprising means for adjusting the time.
読み出し傾斜磁場パルスを印加する手段と、位相エンコ
ード傾斜磁場パルスを印加する手段と、少なくとも1回
の繰り返し時間内で発生するエコー信号を時系列的に検
出する手段と、エコー信号から画像を作成し、表示する
手段と、所定のパルスシーケンスに従い前記各手段を制
御する制御手段とを有する磁気共鳴イメージング装置に
おいて、 前記制御手段は、画像取得のための本計測エコー信号と
は別に、位相エンコード傾斜磁場パルスを印加せずに少
なくとも1回の繰り返し時間内で計測した前計測エコー
信号を用いて、各前計測エコー信号間の位相変化をそれ
ぞれ算出し、それぞれの位相変化を用いて補正用位相マ
ップを作成するとともに、少なくとも1回の繰り返し時
間内で計測した本計測エコー信号をそれぞれ前記補正用
位相マップの対応する補正データを用いて補正する手段
を有する磁気共鳴イメージング装置。2. A means for irradiating a subject with a high-frequency pulse,
Means for applying a read gradient magnetic field pulse, means for applying a phase encoding gradient magnetic field pulse, means for detecting an echo signal generated within at least one repetition time in time series, and creating an image from the echo signal. A magnetic resonance imaging apparatus comprising: a display unit; and a control unit that controls each of the units in accordance with a predetermined pulse sequence. The control unit, apart from the main measurement echo signal for acquiring an image, Using a pre-measurement echo signal measured within at least one repetition time without applying a pulse, a phase change between each pre-measurement echo signal is calculated, and a phase map for correction is calculated using each phase change. While creating, the main measurement echo signal measured within at least one repetition time of each of the correction phase map A magnetic resonance imaging apparatus having means for performing correction using corresponding correction data.
生する静磁場発生手段と、均一磁場空間に置かれた被検
体に対し高周波パルスを照射する高周波パルス発生手段
と、前記被検体に対しスライス方向、位相エンコード方
向および読み出し方向の各傾斜磁場を印加する傾斜磁場
発生手段と、前記被検体からのエコー信号を検出する検
出手段と、前記エコー信号から画像を再構成する画像処
理手段と、得られた画像を表示する画像表示手段と、所
定のシーケンスに従い高周波パルス発生手段、傾斜磁場
発生手段の発生タイミングを制御するとともに検出手
段、画像処理手段、画像表示手段をそれぞれ制御する制
御手段とを備えた磁気共鳴イメージング装置において、 前記制御手段は、位相エンコードを印加しない前計測シ
ーケンスの後に位相エンコードを付加した本計測シーケ
ンスを実行するよう制御し、 前記シーケンスにおいて、少なくとも1回の繰り返し時
間内で連続的に計測した前計測エコー信号を用いて、各
前計測エコー信号間の位相変化をそれぞれ算出するとと
もに、計測空間上に時系列に配置される本計測エコー信
号を対応する前記位相変化の情報に基づいて、偶数番目
取得の本計測エコー信号の位相を奇数番目取得の本計測
エコー信号の位相に揃える、奇数番目取得の本計測エコ
ー信号の位相を偶数番目取得の本計測エコー信号の位相
に揃える、もしくは計測中心の位相に揃える、のいずれ
かの処理を実行する位相補正手段を備えたことを特徴と
する磁気共鳴イメージング装置。3. A static magnetic field generating means for generating a static magnetic field for forming a uniform magnetic field space; a high frequency pulse generating means for irradiating a subject placed in the uniform magnetic field space with a high frequency pulse; On the other hand, gradient magnetic field generating means for applying each gradient magnetic field in the slice direction, phase encoding direction and readout direction, detecting means for detecting an echo signal from the subject, and image processing means for reconstructing an image from the echo signal Image display means for displaying the obtained image, high-frequency pulse generation means according to a predetermined sequence, control means for controlling the generation timing of the gradient magnetic field generation means and control means for detecting means, image processing means, image display means, respectively In the magnetic resonance imaging apparatus provided with, the control means, the phase after the pre-measurement sequence without applying the phase encode Controlling to execute the main measurement sequence to which encoding is added, and in the sequence, using the pre-measurement echo signals continuously measured within at least one repetition time, the phase change between each pre-measurement echo signal is respectively While calculating, based on the information on the phase change corresponding to the main measurement echo signal arranged in time series on the measurement space, the phase of the even measurement acquisition main measurement echo signal is changed to the odd measurement acquisition main measurement echo signal. Phase correction means for performing one of the following processes: aligning the phase of the main measurement echo signal of the odd-number acquisition with the phase of the main measurement echo signal of the even-number acquisition, or aligning the phase with the phase of the measurement center. A magnetic resonance imaging apparatus, characterized in that:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29657299A JP4390328B2 (en) | 1999-10-19 | 1999-10-19 | Magnetic resonance imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29657299A JP4390328B2 (en) | 1999-10-19 | 1999-10-19 | Magnetic resonance imaging system |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2001112735A true JP2001112735A (en) | 2001-04-24 |
| JP2001112735A5 JP2001112735A5 (en) | 2006-12-07 |
| JP4390328B2 JP4390328B2 (en) | 2009-12-24 |
Family
ID=17835286
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP29657299A Expired - Fee Related JP4390328B2 (en) | 1999-10-19 | 1999-10-19 | Magnetic resonance imaging system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4390328B2 (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005007181A (en) * | 2003-06-18 | 2005-01-13 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for improving signal separation of metabolite in mr spectrometry |
| JP2008006117A (en) * | 2006-06-30 | 2008-01-17 | Hitachi Medical Corp | Magnetic resonance imaging system |
| JP2009273530A (en) * | 2008-05-13 | 2009-11-26 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
| US7877129B2 (en) | 2001-09-13 | 2011-01-25 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and RF reception coil apparatus |
| JP2015047507A (en) * | 2013-09-04 | 2015-03-16 | 株式会社東芝 | Magnetic resonance imaging system |
| JP2017100014A (en) * | 2017-03-07 | 2017-06-08 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging apparatus and image processing device |
| JP2018011659A (en) * | 2016-07-19 | 2018-01-25 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging apparatus and installation method |
| WO2021197365A1 (en) * | 2020-03-31 | 2021-10-07 | 浙江大学 | Method and apparatus for frequency drift correction of magnetic resonance cest imaging, and medium and imaging device |
-
1999
- 1999-10-19 JP JP29657299A patent/JP4390328B2/en not_active Expired - Fee Related
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7877129B2 (en) | 2001-09-13 | 2011-01-25 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and RF reception coil apparatus |
| JP2005007181A (en) * | 2003-06-18 | 2005-01-13 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for improving signal separation of metabolite in mr spectrometry |
| JP2008006117A (en) * | 2006-06-30 | 2008-01-17 | Hitachi Medical Corp | Magnetic resonance imaging system |
| JP2009273530A (en) * | 2008-05-13 | 2009-11-26 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
| JP2015047507A (en) * | 2013-09-04 | 2015-03-16 | 株式会社東芝 | Magnetic resonance imaging system |
| JP2018011659A (en) * | 2016-07-19 | 2018-01-25 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging apparatus and installation method |
| JP2017100014A (en) * | 2017-03-07 | 2017-06-08 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging apparatus and image processing device |
| WO2021197365A1 (en) * | 2020-03-31 | 2021-10-07 | 浙江大学 | Method and apparatus for frequency drift correction of magnetic resonance cest imaging, and medium and imaging device |
| JP2023513191A (en) * | 2020-03-31 | 2023-03-30 | 浙江大学 | Magnetic Resonance CEST Imaging Frequency Drift Correction Method, Apparatus, Medium and Imaging Device |
| JP7408184B2 (en) | 2020-03-31 | 2024-01-05 | 浙江大学 | Magnetic Resonance CEST Imaging Frequency Drift Correction Method, Apparatus, Media and Imaging Device |
| US11927659B2 (en) | 2020-03-31 | 2024-03-12 | Zhejiang University | Method and apparatus for frequency drift correction of magnetic resonance CEST imaging, and medium and imaging device |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4390328B2 (en) | 2009-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3544782B2 (en) | Magnetic resonance diagnostic equipment | |
| JP6081273B2 (en) | Slice-specific phase correction in slice multiplexing | |
| US7847546B2 (en) | Magnetic resonance imaging apparatus | |
| JP3529446B2 (en) | Correction method of read gradient magnetic field polarity in EPI and GRASE MRI | |
| JP5607235B2 (en) | Magnetic resonance imaging system | |
| JP4991689B2 (en) | Magnetic resonance imaging apparatus and method | |
| US7418286B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| JPH11113878A (en) | Magnetic resonance imaging method | |
| US20110112393A1 (en) | Magnetic resonance imaging device | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| US6728568B1 (en) | Magnetic resonance imaging method and device | |
| JP5272184B2 (en) | Magnetic resonance imaging system | |
| JP4610611B2 (en) | Magnetic resonance imaging device | |
| US6885885B1 (en) | Magnetic resonance imaging method and device | |
| JP4390328B2 (en) | Magnetic resonance imaging system | |
| JP5285244B2 (en) | Magnetic resonance imaging system | |
| JP4330247B2 (en) | Nuclear magnetic resonance imaging system | |
| JP2004089275A (en) | Phase correction method in magnetic resonance imaging device | |
| JP2002085376A (en) | Nuclear magnetic resonance imaging device and method | |
| JP4678916B2 (en) | Magnetic resonance imaging system | |
| WO2001024695A1 (en) | Nuclear magnetic resonance imaging device and method | |
| JP4832510B2 (en) | Magnetic resonance imaging device | |
| CN111435156B (en) | Method, storage medium and magnetic resonance apparatus for obtaining operating parameters | |
| JPH09238912A (en) | Magnetic resonance imaging equipment | |
| JP4068114B2 (en) | Magnetic resonance imaging system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061019 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061019 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080818 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080826 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081015 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090407 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090603 |
|
| 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: 20091006 |
|
| 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: 20091006 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121016 Year of fee payment: 3 |
|
| 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: 20121016 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131016 Year of fee payment: 4 |
|
| LAPS | Cancellation because of no payment of annual fees |