JP2006255091A - Magnetic resonance imaging system - Google Patents

Magnetic resonance imaging system Download PDF

Info

Publication number
JP2006255091A
JP2006255091A JP2005075292A JP2005075292A JP2006255091A JP 2006255091 A JP2006255091 A JP 2006255091A JP 2005075292 A JP2005075292 A JP 2005075292A JP 2005075292 A JP2005075292 A JP 2005075292A JP 2006255091 A JP2006255091 A JP 2006255091A
Authority
JP
Japan
Prior art keywords
space
magnetic resonance
data included
data
mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2005075292A
Other languages
Japanese (ja)
Other versions
JP4783039B2 (en
Inventor
Masashi Okawa
真史 大川
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.)
Toshiba Corp
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Priority to JP2005075292A priority Critical patent/JP4783039B2/en
Publication of JP2006255091A publication Critical patent/JP2006255091A/en
Application granted granted Critical
Publication of JP4783039B2 publication Critical patent/JP4783039B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To enable proper DC correction by a DC value with slight error. <P>SOLUTION: A control part 107 corrects each data stored in a memory part 104 so as to reduce its direct current component. The control part 107 performs correction based on a DC value which is estimated according to all data included in a K space to which each data belongs. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、磁気共鳴現象を利用して画像を生成する磁気共鳴イメージング(Magnetic Resonance Imaging:MRI)装置に関する。   The present invention relates to a magnetic resonance imaging (MRI) apparatus that generates an image using a magnetic resonance phenomenon.

MRI装置では、収集した信号に含まれる直流(DC)成分を低減するためのDC補正を行っている。   In the MRI apparatus, DC correction for reducing a direct current (DC) component included in the collected signal is performed.

これは、K空間における四隅の一部の固定の割合に相当するデータが示す値の平均値をDC値として推定し、K空間内の全てのデータからそれぞれDC値を差し引くことで行われる。
特開平1−256944号公報
This is performed by estimating an average value of values indicated by data corresponding to a fixed ratio of a part of the four corners in the K space as a DC value, and subtracting the DC value from all the data in the K space.
JP-A-1-256944

ところが、上記のようなDC値の推定は、収集された信号の状態によっては誤差が大きくなることがある。そして大きな誤差を持ったDC値によるDC補正を行うと、DCアーチファクトとなる恐れがあった。   However, the estimation of the DC value as described above may have a large error depending on the state of the collected signal. If DC correction is performed using a DC value having a large error, there is a risk of DC artifacts.

本発明はこのような事情を考慮してなされたものであり、その目的とするところは、誤差の小さなDC値による適正なDC補正を行うことができる磁気共鳴イメージング装置を提供することにある。   The present invention has been made in view of such circumstances, and an object of the present invention is to provide a magnetic resonance imaging apparatus capable of performing appropriate DC correction using a DC value with a small error.

以上の目的を達成するために本発明は、K空間に含まれる各データをそれぞれ、前記K空間に含まれる全てのデータに基づいて直流成分を低減するように補正する補正手段を備えて磁気共鳴イメージング装置を構成した。   In order to achieve the above object, the present invention comprises correction means for correcting each data included in the K space so as to reduce a direct current component based on all the data included in the K space, and magnetic resonance. An imaging device was constructed.

本発明によれば、誤差の小さなDC値による適正なDC補正を行うことが可能となる。   According to the present invention, it is possible to perform appropriate DC correction using a DC value with a small error.

以下、図面を参照して本発明の一実施形態について説明する。
図1は本実施形態に係る磁気共鳴イメージング装置(以下、MRI装置と称する)の構成を示す図である。この図1に示すMRI装置は、静磁場磁石1、傾斜磁場コイル2、傾斜磁場電源3、寝台4、寝台制御部5、RFコイルユニット6、送信部7、受信部8、ハイブリッド回路9および計算機システム10を具備する。
Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
FIG. 1 is a diagram showing a configuration of a magnetic resonance imaging apparatus (hereinafter referred to as an MRI apparatus) according to the present embodiment. The MRI apparatus shown in FIG. 1 includes a static magnetic field magnet 1, a gradient magnetic field coil 2, a gradient magnetic field power supply 3, a bed 4, a bed control unit 5, an RF coil unit 6, a transmission unit 7, a reception unit 8, a hybrid circuit 9, and a computer. A system 10 is provided.

静磁場磁石1は、中空の円筒形をなし、内部の空間に一様な静磁場を発生する。この静磁場磁石1としては、例えば永久磁石、超伝導磁石等が使用される。   The static magnetic field magnet 1 has a hollow cylindrical shape and generates a uniform static magnetic field in an internal space. As the static magnetic field magnet 1, for example, a permanent magnet, a superconducting magnet or the like is used.

傾斜磁場コイル2は、中空の円筒形をなし、静磁場磁石1の内側に配置される。傾斜磁場コイル2は、互いに直交するX,Y,Zの各軸に対応する3つのコイルが組み合わされている。傾斜磁場コイル2は、上記の3つのコイルが傾斜磁場電源3から個別に電流供給を受けて、磁場強度がX,Y,Zの各軸に沿って傾斜する傾斜磁場を発生する。なお、Z軸方向は、例えば静磁場と同方向とする。X,Y,Z各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Geおよびリードアウト用傾斜磁場Grにそれぞれ対応される。スライス選択用傾斜磁場Gsは、任意に撮影断面を決めるために利用される。位相エンコード用傾斜磁場Geは、空間的位置に応じて磁気共鳴信号の位相をエンコードするために利用される。リードアウト用傾斜磁場Grは、空間的位置に応じて磁気共鳴信号の周波数をエンコードするために利用される。   The gradient magnetic field coil 2 has a hollow cylindrical shape and is disposed inside the static magnetic field magnet 1. The gradient coil 2 is a combination of three coils corresponding to the X, Y, and Z axes orthogonal to each other. The gradient magnetic field coil 2 generates a gradient magnetic field in which the above three coils are individually supplied with electric current from the gradient magnetic field power supply 3 and the magnetic field strength is inclined along the X, Y, and Z axes. The Z-axis direction is, for example, the same direction as the static magnetic field. The gradient magnetic fields of the X, Y, and Z axes correspond to, for example, the slice selection gradient magnetic field Gs, the phase encoding gradient magnetic field Ge, and the readout gradient magnetic field Gr, respectively. The slice selection gradient magnetic field Gs is used to arbitrarily determine an imaging section. The phase encoding gradient magnetic field Ge is used to encode the phase of the magnetic resonance signal in accordance with the spatial position. The readout gradient magnetic field Gr is used to encode the frequency of the magnetic resonance signal in accordance with the spatial position.

被検体Pは、寝台4の天板41に載置された状態で傾斜磁場コイル2の空洞(撮影口)内に挿入される。寝台4は、寝台制御部5により駆動され、天板41をその長手方向(図1中における左右方向)および上下方向に移動する。通常、この長手方向が静磁場磁石1の中心軸と平行になるように寝台4が設置される。   The subject P is inserted into the cavity (imaging port) of the gradient coil 2 while being placed on the top 41 of the bed 4. The couch 4 is driven by the couch controller 5 and moves the top board 41 in the longitudinal direction (left-right direction in FIG. 1) and up-down direction. Usually, the bed 4 is installed such that the longitudinal direction thereof is parallel to the central axis of the static magnetic field magnet 1.

RFコイルユニット6は、傾斜磁場コイル2の内側に配置される。RFコイルユニット6は、送信部7から高周波パルス(RFパルス)の供給を受けて、高周波磁場を発生する。RFコイルユニット6は、被検体から放射される磁気共鳴信号を受信する。RFコイルユニット6は、送信用コイルおよび受信用コイルを個別に有していても良いし、送受信共用のコイルを有していても良い。また、送信用コイル、受信用コイルおよび送受信共用のコイルは、1つのみが設けられても良いし、送受信範囲や送受信特性が異なる複数種類が並列的に設けられていても良い。   The RF coil unit 6 is disposed inside the gradient magnetic field coil 2. The RF coil unit 6 receives a high frequency pulse (RF pulse) from the transmission unit 7 and generates a high frequency magnetic field. The RF coil unit 6 receives a magnetic resonance signal radiated from the subject. The RF coil unit 6 may have a transmission coil and a reception coil individually, or may have a coil for both transmission and reception. Further, only one transmission coil, reception coil, and transmission / reception coil may be provided, or a plurality of types having different transmission / reception ranges and transmission / reception characteristics may be provided in parallel.

送信部7は、発振部、位相選択部、周波数変換部、振幅変調部および高周波電力増幅部(いずれも図示せず)を有している。発振部は、静磁場中における対象原子核に固有の共鳴周波数の高周波信号を発生する。位相選択部は、上記高周波信号の位相を選択する。周波数変換部は、位相選択部から出力された高周波信号の周波数を変換する。振幅変調部は、周波数変調部から出力された高周波信号の振幅を例えばシンク関数に従って変調する。高周波電力増幅部は、振幅変調部から出力された高周波信号を増幅する。そしてこれらの各部の動作の結果として送信部7は、ラーモア周波数に対応するRFパルスをRFコイルユニット6に供給するべく送出する。   The transmission unit 7 includes an oscillation unit, a phase selection unit, a frequency conversion unit, an amplitude modulation unit, and a high-frequency power amplification unit (all not shown). The oscillation unit generates a high-frequency signal having a resonance frequency unique to the target nucleus in the static magnetic field. The phase selection unit selects the phase of the high-frequency signal. The frequency conversion unit converts the frequency of the high-frequency signal output from the phase selection unit. The amplitude modulation unit modulates the amplitude of the high-frequency signal output from the frequency modulation unit, for example, according to a sync function. The high frequency power amplification unit amplifies the high frequency signal output from the amplitude modulation unit. As a result of the operation of each of these units, the transmission unit 7 transmits an RF pulse corresponding to the Larmor frequency to be supplied to the RF coil unit 6.

受信部8は、前段増幅器、位相検波器およびアナログディジタル変換器(いずれも図示せず)を有している。前段増幅器は、ハイブリッド回路9から出力される磁気共鳴信号を増幅する。位相検波器は、前置増幅器から出力される磁気共鳴信号の位相を検波する。アナログディジタル変換器は、位相検波器から出力される信号をディジタル信号に変換する。   The receiving unit 8 includes a pre-stage amplifier, a phase detector, and an analog / digital converter (all not shown). The preamplifier amplifies the magnetic resonance signal output from the hybrid circuit 9. The phase detector detects the phase of the magnetic resonance signal output from the preamplifier. The analog-digital converter converts the signal output from the phase detector into a digital signal.

ハイブリッド回路9は、送信部7から送出される高周波パルスを、送信期間にてRFコイルユニット6へ供給する。ハイブリッド回路9は、RFコイルユニット6から出力される信号を、受信期間にて受信部8へ供給する。送信期間および受信期間は、計算機システム10から指示される。またハイブリッド回路9は、RFコイルユニット6に受信用のコイルが複数設けられている場合に、それらのコイルが出力する信号を選択的に、並列的に、あるいは時分割的に入力する。どのコイルから出力される信号をどのような形態で入力するかは、例えば計算機システム10から指示される。   The hybrid circuit 9 supplies the RF coil unit 6 with the high-frequency pulse transmitted from the transmission unit 7 during the transmission period. The hybrid circuit 9 supplies the signal output from the RF coil unit 6 to the reception unit 8 during the reception period. The transmission period and the reception period are instructed from the computer system 10. Further, when a plurality of receiving coils are provided in the RF coil unit 6, the hybrid circuit 9 selectively inputs signals output from these coils in parallel or in a time division manner. It is instructed from the computer system 10, for example, in what form the signal output from which coil is input.

計算機システム10は、インタフェース部101、データ収集部102、再構成部103、記憶部104、表示部105、入力部106および制御部107を有している。   The computer system 10 includes an interface unit 101, a data collection unit 102, a reconstruction unit 103, a storage unit 104, a display unit 105, an input unit 106, and a control unit 107.

インタフェース部101には、傾斜磁場電源3、寝台制御部5、送信部7、受信部8およびハイブリッド回路9等が接続される。インタフェース部101は、これらの接続された各部と計算機システム10との間で授受される信号の入出力を行う。   The interface unit 101 is connected to the gradient magnetic field power supply 3, the bed control unit 5, the transmission unit 7, the reception unit 8, the hybrid circuit 9, and the like. The interface unit 101 inputs and outputs signals exchanged between these connected units and the computer system 10.

データ収集部102は、受信部8から出力されるディジタル信号を収集する。データ収集部102は、収集したディジタル信号、すなわち磁気共鳴信号データを、記憶部104に格納する。   The data collection unit 102 collects digital signals output from the reception unit 8. The data collection unit 102 stores the collected digital signal, that is, magnetic resonance signal data, in the storage unit 104.

再構成部103は、記憶部104に記憶された磁気共鳴信号データに対して、後処理、すなわちフーリエ変換等の再構成を実行し、被検体P内の所望核スピンのスペクトラムデータあるいは画像データを求める。   The reconstruction unit 103 performs post-processing, that is, reconstruction such as Fourier transform, on the magnetic resonance signal data stored in the storage unit 104, and obtains spectrum data or image data of the desired nuclear spin in the subject P. Ask.

記憶部104は、磁気共鳴信号データと、スペクトラムデータあるいは画像データとを、患者毎に記憶する。   The storage unit 104 stores magnetic resonance signal data and spectrum data or image data for each patient.

表示部105は、スペクトラムデータあるいは画像データ等の各種の情報を制御部107の制御の下に表示する。表示部105としては、液晶表示器などの表示デバイスを利用可能である。   The display unit 105 displays various information such as spectrum data or image data under the control of the control unit 107. As the display unit 105, a display device such as a liquid crystal display can be used.

入力部106は、オペレータからの各種指令や情報入力を受け付ける。入力部106としては、マウスやトラックボールなどのポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスを適宜に利用可能である。   The input unit 106 receives various commands and information input from the operator. As the input unit 106, a pointing device such as a mouse or a trackball, a selection device such as a mode change switch, or an input device such as a keyboard can be used as appropriate.

制御部107は、図示していないCPUやメモリ等を有しており、本実施形態のMRI装置を総括的に制御する。また制御部107は、MRI装置が一般的に備える機能を実現するための手段としての他に、以下のようないくつかの機能を実現するための手段として働く。この機能の一つは、収集された信号に含まれるDC値を推定する。また上記の機能の一つは、上記の推定したDC値に基づいてDC補正を行う。   The control unit 107 includes a CPU, a memory, and the like (not shown), and comprehensively controls the MRI apparatus according to the present embodiment. Further, the control unit 107 serves as means for realizing some functions as described below in addition to means for realizing functions generally provided in the MRI apparatus. One of these functions estimates the DC value contained in the collected signal. One of the above functions performs DC correction based on the estimated DC value.

次に以上のように構成されたMRI装置の動作について説明する。なお、被検体Pの画像を得るための動作は従来通りであるのでここでは説明を省略する。そして以下では、DC補正に係わる処理について説明する。   Next, the operation of the MRI apparatus configured as described above will be described. Since the operation for obtaining the image of the subject P is the same as the conventional one, the description thereof is omitted here. In the following, processing related to DC correction will be described.

制御部107は、データ収集部102によって収集されて記憶部104に記憶された磁気共鳴信号データについて、DC補正を行う。DC補正は、1つのK空間に属するデータを1単位として行われる。実部チャネルおよび虚部チャネルのそれぞれについても個別に処理される。   The control unit 107 performs DC correction on the magnetic resonance signal data collected by the data collection unit 102 and stored in the storage unit 104. The DC correction is performed with data belonging to one K space as one unit. Each of the real part channel and the imaginary part channel is also processed individually.

図2はDC補正のための制御部107の処理手順を示すフローチャートである。なお図2は、1つのK空間に属するデータについてのDC補正を行う処理を示している。従って、実部チャネルおよび虚部チャネルのそれぞれのK空間およびスライスや時相が異なるK空間に属する磁気共鳴信号データに関しては、それぞれに対して同様な処理が行われる。   FIG. 2 is a flowchart showing a processing procedure of the control unit 107 for DC correction. FIG. 2 shows a process for performing DC correction on data belonging to one K space. Therefore, the same processing is performed on the magnetic resonance signal data belonging to the K space and the slices of the real part channel and the imaginary part channel and the K spaces having different time phases.

ステップSa1において制御部107は、DWI(Diffusion Weighted Imaging)法を使用して再構成を行うか否かを確認する。当該MRI装置がDWI法を使用するモードに設定されているならば、制御部107はステップSa1からステップSa5へ進むが、そうでなければステップSa2へ進む。   In step Sa1, the control unit 107 confirms whether or not reconstruction is performed using a DWI (Diffusion Weighted Imaging) method. If the MRI apparatus is set to a mode that uses the DWI method, the control unit 107 proceeds from step Sa1 to step Sa5, but otherwise proceeds to step Sa2.

ステップSa2において制御部107は、EPI(Echo Planar Imaging)法を使用して再構成を行うか否かを確認する。当該MRI装置がEPI法を使用するモードに設定されているならば、制御部107はステップSa2からステップSa5へ進むが、そうでなければステップSa3へ進む。   In step Sa2, the control unit 107 confirms whether or not reconstruction is performed using an EPI (Echo Planar Imaging) method. If the MRI apparatus is set to a mode that uses the EPI method, the control unit 107 proceeds from step Sa2 to step Sa5, but otherwise proceeds to step Sa3.

ステップSa3において制御部107は、信号強度が閾値以下である否かを確認する。信号強度は、処理対象となっているK空間の範囲内での磁気共鳴信号データが表す信号の強度を表す数値である。信号強度としては例えば、処理対象となっているK空間に属する複数のデータが表す信号レベルのうちのピーク値、あるいはK空間に属する全データが表す信号レベルの平均値、または実部および虚部のそれぞれの平均値を適用できる。そして信号強度が閾値以下であるならば、制御部107はステップSa3からステップSa5へ進むが、そうでなければステップSa4へ進む。   In step Sa3, the control unit 107 confirms whether the signal strength is equal to or less than a threshold value. The signal strength is a numerical value representing the strength of the signal represented by the magnetic resonance signal data within the range of the K space to be processed. As the signal intensity, for example, a peak value of signal levels represented by a plurality of data belonging to the K space to be processed, an average value of signal levels represented by all data belonging to the K space, or a real part and an imaginary part Each average value of can be applied. If the signal strength is less than or equal to the threshold value, the control unit 107 proceeds from step Sa3 to step Sa5, but otherwise proceeds to step Sa4.

ステップSa4において制御部107は、ばらつき度が閾値以上であるか否かを確認する。ばらつき度は、処理対象となっているK空間の範囲内での磁気共鳴信号データが表す強度のばらつきの大きさを表す数値である。ばらつき度としては例えば、処理対象となっているK空間に属する複数のデータが表す信号レベルのうちのピーク値以外のレベルの偏差、または実部および虚部のピーク値以外のレベルの偏差を適用できる。そしてばらつき度が閾値以上であるならば、制御部107はステップSa4からステップSa5へ進むが、そうでなければステップSa7へ進む。   In step Sa4, the control unit 107 confirms whether or not the degree of variation is equal to or greater than a threshold value. The degree of variation is a numerical value indicating the magnitude of intensity variation represented by the magnetic resonance signal data within the range of the K space to be processed. As the degree of variation, for example, a deviation of a level other than the peak value among signal levels represented by a plurality of data belonging to the K space to be processed, or a deviation of a level other than the peak value of the real part and the imaginary part is applied. it can. If the degree of variation is greater than or equal to the threshold value, the control unit 107 proceeds from step Sa4 to step Sa5, but otherwise proceeds to step Sa7.

ステップSa5において制御部107は、処理対象となっているK空間の範囲内の全データをDC値の推定用に割り当てる。こののちに制御部107は、ステップSa7へ進む。このときの推定用データの範囲とK空間との関係は図3に示すようになる。なお図3では、正方形の枠がK空間の外縁を示し、ハッチングした領域が推定用データの範囲を示す。   In step Sa5, the control unit 107 assigns all data within the range of the K space to be processed for estimation of the DC value. After this, the control unit 107 proceeds to step Sa7. The relationship between the estimation data range and the K space at this time is as shown in FIG. In FIG. 3, the square frame indicates the outer edge of the K space, and the hatched area indicates the range of the estimation data.

一方、ステップSa6において制御部107は、処理対象となっているK空間の範囲の四隅の一部に相当するデータをDC値の推定用に割り当てる。こののちに制御部107は、ステップSa7へ進む。このときの推定用データの範囲とK空間との関係は図4に示すようになる。なお図4では、外側の正方形の枠がK空間の外縁を示し、ハッチングした領域が推定用データの範囲を示す。図4の例は、従来と同様に推定用データの範囲をK空間の一部の固定の割合の部分としている。このときの推定用データの範囲は、例えば図5や図6にハッチングで示すような範囲としても良い。すなわち、少なくともK空間の中心に相当するデータは、推定用データに割り当てないようにする。また、信号強度の大きさやばらつき度の大きさに応じて推定用データの範囲を変化させても良い。   On the other hand, in step Sa6, the control unit 107 allocates data corresponding to part of the four corners of the range of the K space to be processed for estimation of the DC value. After this, the control unit 107 proceeds to step Sa7. The relationship between the estimation data range and the K space at this time is as shown in FIG. In FIG. 4, the outer square frame indicates the outer edge of the K space, and the hatched area indicates the range of the estimation data. In the example of FIG. 4, the estimation data range is set to a fixed part of the K space as in the conventional case. The range of the estimation data at this time may be, for example, a range shown by hatching in FIG. 5 or FIG. That is, at least data corresponding to the center of the K space is not assigned to the estimation data. Further, the range of the estimation data may be changed according to the magnitude of the signal intensity and the degree of variation.

ステップSa7において制御部107は、推定用に割り当てられたデータに基づいてDC値を推定する。   In step Sa7, the control unit 107 estimates the DC value based on the data assigned for estimation.

ステップSa8において制御部107は、処理対象となっているK空間に含まれる全データから上記の推定したDC値をおのおの差し引く。   In step Sa8, the control unit 107 subtracts the estimated DC value from all data included in the K space to be processed.

かくして本実施形態によれば、DWI法を使用する場合には、K空間内の全てのデータに基づいてDC値の推定がなされる。DWI法では、収集される信号の強度が全般的に小さく、かつピーク値も周囲の強度に比べてさほど大きくならないので、多数のデータを参照することでDC値の推定精度が向上する。   Thus, according to the present embodiment, when the DWI method is used, the DC value is estimated based on all the data in the K space. In the DWI method, the intensity of the collected signal is generally small, and the peak value is not so large as compared with the surrounding intensity, so that the estimation accuracy of the DC value is improved by referring to a large number of data.

また、EPI法を使用する場合には、K空間内の全てのデータに基づいてDC値の推定がなされる。EPI法では、信号強度はある程度大きくなるが、信号のばらつきが大きくなることから、多数のデータを参照することでDC値の推定精度が向上する。   When the EPI method is used, the DC value is estimated based on all data in the K space. In the EPI method, the signal intensity is increased to some extent, but the variation in the signal is increased. Therefore, the estimation accuracy of the DC value is improved by referring to a large number of data.

また、DWI法を使用しない場合でも、信号強度が閾値以下となる程に小さい場合には、K空間内の全てのデータに基づいてDC値の推定がなされる。信号強度が小さいのであるから、多数のデータを参照することでDC値の推定精度が向上する。   Even when the DWI method is not used, if the signal intensity is small enough to be equal to or less than the threshold value, the DC value is estimated based on all data in the K space. Since the signal strength is small, the estimation accuracy of the DC value is improved by referring to a large number of data.

また、EPI法を使用しない場合でも、信号のばらつき度が閾値を超える程に信号のばらつきが大きいのであれば、K空間内の全てのデータに基づいてDC値の推定がなされる。信号のばらつきが大きいのであるから、多数のデータを参照することでDC値の推定精度が向上する。   Even if the EPI method is not used, if the signal variation is so large that the signal variation exceeds the threshold value, the DC value is estimated based on all the data in the K space. Since the signal variation is large, the estimation accuracy of the DC value is improved by referring to a large number of data.

そして、上記の4つの条件のいずれにも合致しないときには、K空間の中心に相当するデータをDC値の推定のために参照しない。これにより、ピーク値が非常に大きくなっても、DC値の推定を正確に行うことが可能となる。   When none of the above four conditions is met, the data corresponding to the center of the K space is not referred to for estimating the DC value. This makes it possible to accurately estimate the DC value even when the peak value becomes very large.

以上のように本実施形態によれば、DC値を常に小さな誤差で推定することができる。そしてこのように誤差の小さなDC値に基づくことにより、適正なDC補正が行える。この結果、DCアーチファクトの発生を抑圧することができる。   As described above, according to the present embodiment, the DC value can always be estimated with a small error. Based on the DC value with a small error in this way, appropriate DC correction can be performed. As a result, the occurrence of DC artifact can be suppressed.

この実施形態は、次のような種々の変形実施が可能である。
ステップSa1乃至ステップSa4の任意の1つ乃至3つの条件判断を省略しても良い。
あるいは、無条件で常にK空間内の全てのデータに基づいてDC値の推定を行うようにしても良い。
This embodiment can be variously modified as follows.
Arbitrary one to three condition judgments in steps Sa1 to Sa4 may be omitted.
Alternatively, the DC value may be estimated unconditionally based on all data in the K space.

なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。   Note that the present invention is not limited to the above-described embodiment as it is, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment.

本発明の一実施形態に係る磁気共鳴イメージング装置の構成を示す図。The figure which shows the structure of the magnetic resonance imaging apparatus which concerns on one Embodiment of this invention. DC補正のための図1中の制御部107の処理手順を示すフローチャート。The flowchart which shows the process sequence of the control part 107 in FIG. 1 for DC correction | amendment. 図2中のステップSa5にて割り当てる推定用データの範囲とK空間との関係を示す図。The figure which shows the relationship between the range of the data for estimation allocated in step Sa5 in FIG. 2, and K space. 図2中のステップSa6にて割り当てる推定用データの範囲とK空間との関係を示す図。The figure which shows the relationship between the range of the data for estimation allocated in step Sa6 in FIG. 2, and K space. 図2中のステップSa6にて割り当てる推定用データの範囲とK空間との関係の変形例を示す図。The figure which shows the modification of the relationship between the range of the data for estimation allocated in step Sa6 in FIG. 2, and K space. 図2中のステップSa6にて割り当てる推定用データの範囲とK空間との関係の変形例を示す図。The figure which shows the modification of the relationship between the range of the data for estimation allocated in step Sa6 in FIG. 2, and K space.

符号の説明Explanation of symbols

1…静磁場磁石、2…傾斜磁場コイル、3…傾斜磁場電源、4…寝台、5…寝台制御部、6…RFコイルユニット、7…送信部、8…受信部、9…ハイブリッド回路、10…計算機システム、101…インタフェース部、102…データ収集部、103…再構成部、104…記憶部、105…表示部、106…入力部、107…制御部。   DESCRIPTION OF SYMBOLS 1 ... Static magnetic field magnet, 2 ... Gradient magnetic field coil, 3 ... Gradient magnetic field power supply, 4 ... Bed, 5 ... Bed control part, 6 ... RF coil unit, 7 ... Transmission part, 8 ... Reception part, 9 ... Hybrid circuit, 10 DESCRIPTION OF SYMBOLS Computer system 101 ... Interface part 102 ... Data collection part 103 ... Reconstruction part 104 ... Storage part 105 ... Display part 106 ... Input part 107 ... Control part

Claims (7)

K空間に含まれる各データをそれぞれ、前記K空間に含まれる全てのデータに基づいて直流成分を低減するように補正する補正手段を備えたことを特徴とする磁気共鳴イメージング装置。   A magnetic resonance imaging apparatus comprising correction means for correcting each data included in the K space so as to reduce a direct current component based on all the data included in the K space. 前記補正手段は、
前記K空間に含まれる全てのデータに基づいて前記直流成分の大きさを推定する手段と、
この推定された前記直流成分の大きさを前記K空間に含まれる各データから差し引く手段とを具備することを特徴とする請求項1に記載の磁気共鳴イメージング装置。
The correction means includes
Means for estimating the magnitude of the DC component based on all data contained in the K space;
2. The magnetic resonance imaging apparatus according to claim 1, further comprising means for subtracting the estimated magnitude of the DC component from each data included in the K space.
前記K空間に含まれるデータを取得する条件または前記K空間に含まれるデータの状態に基づいて第1のモードまたは第2のモードを選択する選択手段をさらに備え、
かつ前記補正手段は前記K空間に含まれる各データの直流成分を、前記第1のモードが選択されている時には前記K空間に含まれる全てのデータに基づいて低減し、前記第2のモードが選択されている時には前記K空間に含まれるデータの一部のみに基づいて低減することを特徴とする請求項1に記載の磁気共鳴イメージング装置。
Further comprising selection means for selecting the first mode or the second mode based on a condition for acquiring data included in the K space or a state of data included in the K space;
The correction means reduces the direct current component of each data included in the K space based on all the data included in the K space when the first mode is selected, and the second mode The magnetic resonance imaging apparatus according to claim 1, wherein when selected, the magnetic resonance imaging apparatus is reduced based on only a part of data included in the K space.
前記選択手段は、前記K空間に含まれるデータが表す磁気共鳴信号の強度に基づいて前記第1のモードまたは前記第2のモードを選択することを特徴とする請求項3に記載の磁気共鳴イメージング装置。   The magnetic resonance imaging according to claim 3, wherein the selection unit selects the first mode or the second mode based on an intensity of a magnetic resonance signal represented by data included in the K space. apparatus. 前記選択手段は、前記K空間に含まれる複数のデータが表す強度のばらつきに基づいて前記第1のモードまたは前記第2のモードを選択することを特徴とする請求項3に記載の磁気共鳴イメージング装置。   The magnetic resonance imaging according to claim 3, wherein the selection unit selects the first mode or the second mode based on variation in intensity represented by a plurality of data included in the K space. apparatus. 前記選択手段は、前記K空間に含まれるデータを取得するためにEPI(Echo Planar Imaging)法を使用する場合に前記第1のモードを選択することを特徴とする請求項3に記載の磁気共鳴イメージング装置。   The magnetic resonance according to claim 3, wherein the selection unit selects the first mode when an EPI (Echo Planar Imaging) method is used to acquire data included in the K space. Imaging device. 前記選択手段は、前記K空間に含まれるデータを取得するためにDWI(Diffusion Weighted Imaging)法を使用する場合に前記第1のモードを選択することを特徴とする請求項3に記載の磁気共鳴イメージング装置。   4. The magnetic resonance according to claim 3, wherein the selection unit selects the first mode when a DWI (Diffusion Weighted Imaging) method is used to acquire data included in the K space. 5. Imaging device.
JP2005075292A 2005-03-16 2005-03-16 Magnetic resonance imaging system Expired - Fee Related JP4783039B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005075292A JP4783039B2 (en) 2005-03-16 2005-03-16 Magnetic resonance imaging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005075292A JP4783039B2 (en) 2005-03-16 2005-03-16 Magnetic resonance imaging system

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2011120865A Division JP5439431B2 (en) 2011-05-30 2011-05-30 Magnetic resonance imaging system

Publications (2)

Publication Number Publication Date
JP2006255091A true JP2006255091A (en) 2006-09-28
JP4783039B2 JP4783039B2 (en) 2011-09-28

Family

ID=37094933

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005075292A Expired - Fee Related JP4783039B2 (en) 2005-03-16 2005-03-16 Magnetic resonance imaging system

Country Status (1)

Country Link
JP (1) JP4783039B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108535675A (en) * 2018-04-08 2018-09-14 朱高杰 A kind of magnetic resonance multichannel method for reconstructing being in harmony certainly based on deep learning and data

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55112556A (en) * 1979-02-22 1980-08-30 Jeol Ltd Correction method of level in fourier spectrum unit
JPH01195849A (en) * 1988-01-31 1989-08-07 Shimadzu Corp Direct current constituent detecting method for nuclear magnetic resonance signal
JPH01250235A (en) * 1988-03-31 1989-10-05 Shimadzu Corp Apparatus for detecting dc component of signal of mri apparatus
JPH01256944A (en) * 1988-04-07 1989-10-13 Hitachi Medical Corp Magnetic resonance induction device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55112556A (en) * 1979-02-22 1980-08-30 Jeol Ltd Correction method of level in fourier spectrum unit
JPH01195849A (en) * 1988-01-31 1989-08-07 Shimadzu Corp Direct current constituent detecting method for nuclear magnetic resonance signal
JPH01250235A (en) * 1988-03-31 1989-10-05 Shimadzu Corp Apparatus for detecting dc component of signal of mri apparatus
JPH01256944A (en) * 1988-04-07 1989-10-13 Hitachi Medical Corp Magnetic resonance induction device

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108535675A (en) * 2018-04-08 2018-09-14 朱高杰 A kind of magnetic resonance multichannel method for reconstructing being in harmony certainly based on deep learning and data
CN108535675B (en) * 2018-04-08 2020-12-04 朱高杰 Magnetic resonance multi-channel reconstruction method based on deep learning and data self-consistency

Also Published As

Publication number Publication date
JP4783039B2 (en) 2011-09-28

Similar Documents

Publication Publication Date Title
US7769425B2 (en) Magnetic resonance device and method
US10393840B2 (en) Magnetic resonance apparatus and method for acquiring measurement data during simultaneous manipulation of spatially separate sub-volumes
JP5171373B2 (en) Magnetic resonance imaging apparatus and RF pulse applying method
JP5366437B2 (en) Magnetic resonance imaging system
JP2008012172A (en) Magnetic resonance imaging system and image processor
JP5559848B2 (en) Apparatus and method for simultaneously generating multiple types of magnetic resonance images
JP2011120873A (en) Magnetic resonance imaging apparatus and method
JP2010099455A (en) Magnetic resonance diagnostic apparatus and magnetic resonance diagnostic method
KR101743163B1 (en) Method and magnetic resonance system for reconstructing an mr image taking into account the chemical shift
KR20130018565A (en) Method to acquire mr data in a predetermined region, and corresponding magnetic resonance system
JP5596215B2 (en) Magnetic resonance diagnostic equipment
JP2009261421A (en) Magnetic resonance imaging system
US11360171B2 (en) Method for obtaining magnetic resonance imaging data and magnetic resonance imaging system
JP2010075634A (en) Magnetic resonance imaging system
US20140320128A1 (en) Method and magnetic resonance apparatus to acquire image data sets of an examination subject
JP4783039B2 (en) Magnetic resonance imaging system
JP5439431B2 (en) Magnetic resonance imaging system
JP6996865B2 (en) Magnetic resonance imaging device
JP4912156B2 (en) Magnetic resonance imaging method and apparatus
JP6181374B2 (en) Magnetic resonance imaging system
JP2007202584A (en) Magnetic resonance imaging apparatus
JP2023135351A (en) Radio frequency coil apparatus and method for controlling radio frequency coil apparatus
JP5417425B2 (en) Magnetic resonance diagnostic equipment
JP5509375B2 (en) Magnetic resonance apparatus and high frequency coil unit
CN113534018A (en) Linear compensation method and device of radio frequency amplifier and magnetic resonance imaging system

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080226

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100826

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110329

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110530

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

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

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

Free format text: PAYMENT UNTIL: 20140715

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4783039

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

Free format text: JAPANESE INTERMEDIATE CODE: R313115

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees