JP3112930B2 - Magnetic resonance equipment - Google Patents
Magnetic resonance equipmentInfo
- Publication number
- JP3112930B2 JP3112930B2 JP02340398A JP34039890A JP3112930B2 JP 3112930 B2 JP3112930 B2 JP 3112930B2 JP 02340398 A JP02340398 A JP 02340398A JP 34039890 A JP34039890 A JP 34039890A JP 3112930 B2 JP3112930 B2 JP 3112930B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- eddy current
- magnetic resonance
- resonance signal
- gradient
- 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.)
- Expired - Fee Related
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【発明の詳細な説明】 ［発明の目的］ （産業上の利用分野） 本発明は磁気共鳴装置に係り、特に勾配磁場コイル駆
動時に発生する渦電流磁場の計測方法に関する。The present invention relates to a magnetic resonance apparatus, and more particularly to a method for measuring an eddy current magnetic field generated when a gradient magnetic field coil is driven.
（従来の技術） 磁気共鳴映像法は、固有の磁気モーメントを持つ核の
集団が一様な静磁場中に置かれたときに、特定の周波数
で回転する高周波磁場のエネルギーを共鳴的に吸収する
現象を利用して、物質の化学的および物理的な微視的情
報を映像化する手法である。この方法では磁気共鳴信号
に空間情報（位置情報）を持たせるために、勾配磁場の
使用が必要不可欠である。勾配磁場はスイッチングされ
て印加されるため、そのスイッチングにより勾配コイル
外部のシールド円筒導体上に渦電流が誘起される。この
渦電流によって生ずる磁場（以下、渦電流磁場という）
により磁気共鳴信号が変調を受け、画像上ではぼけ等の
劣化が生じ、スペクトル上では歪み等の劣化が生じる要
因となる。(Prior art) Magnetic resonance imaging resonatesly absorbs the energy of a high frequency magnetic field rotating at a specific frequency when a group of nuclei having a unique magnetic moment is placed in a uniform static magnetic field. This is a method of visualizing chemical and physical microscopic information of a substance using a phenomenon. In this method, the use of a gradient magnetic field is indispensable for giving spatial information (position information) to the magnetic resonance signal. Since the gradient magnetic field is switched and applied, the switching induces an eddy current on the shield cylindrical conductor outside the gradient coil. The magnetic field generated by this eddy current (hereinafter referred to as eddy current magnetic field)
Causes a magnetic resonance signal to be modulated, which causes deterioration such as blur on an image and causes deterioration such as distortion on a spectrum.
この問題を解決する方法として、 （ａ）渦電流の時間応答の逆応答に相当する成分で勾配
コイル電流を変調し、渦電流磁場を補償する方法 （ｂ）勾配コイルの外部に、漏洩磁場遮断のためのアク
ティブシールド勾配コイル（ASGC）を設置する方法 が提案されている。As a method for solving this problem, (a) a method of modulating the gradient coil current with a component corresponding to the inverse response of the time response of the eddy current to compensate for the eddy current magnetic field, and (b) blocking the leakage magnetic field outside the gradient coil A method has been proposed to install an active shield gradient coil (ASGC).
しかしながら、（ａ）の方法では渦電流によって形成
される勾配磁場と勾配コイル電流の形成するそれとの間
に異なった空間的非線形性が存在するため、渦電流磁場
の補償可能な領域が限定され、この範囲で測定して逆応
答を求める必要が生ずる。また、この方法を有効にする
ためには、勾配コイル中心と渦電流磁場中心を一致させ
る、いわゆる渦電流磁場中心合わせを行わなければなら
ないので、勾配コイル中心から測定点の距離の精度が必
要とされる。従って、渦電流磁場中心合わせに対して、
位置精度向上のための特別な治具を必要とする。また、
実際には複数のシールド円筒導体が存在し、勾配コイル
の中心と各シールド導体から発生する渦電流磁場の中心
を合わせることはかなりの困難を伴い、補償しきれない
渦電流磁場（以下、残留渦電流磁場という）が残ってし
まう。However, in the method (a), since a different spatial nonlinearity exists between the gradient magnetic field formed by the eddy current and that generated by the gradient coil current, the area where the eddy current magnetic field can be compensated is limited, It is necessary to measure in this range to obtain a reverse response. Also, in order to make this method effective, the center of the gradient coil and the center of the eddy current magnetic field must be matched, that is, so-called eddy current magnetic field center alignment must be performed. Is done. Therefore, for centering of the eddy current magnetic field,
A special jig for improving the position accuracy is required. Also,
In practice, there are several shielded cylindrical conductors, and it is quite difficult to align the center of the gradient coil with the center of the eddy current magnetic field generated from each shield conductor. Current magnetic field) remains.
また、（ｂ）の方法では勾配コイルの中心とアクティ
ブシールド勾配コイルの中心を厳密に一致させることが
難しく、その結果やはり渦電流磁場が残ってしまう。In the method (b), it is difficult to exactly match the center of the gradient coil with the center of the active shield gradient coil. As a result, an eddy current magnetic field remains.
さらに、（ａ）（ｂ）のいずれの方法においても、渦
電流磁場対策後の勾配磁場スイッチング時の磁場の評価
を行う必要がある。Further, in any of the methods (a) and (b), it is necessary to evaluate the magnetic field at the time of the gradient magnetic field switching after the eddy current magnetic field countermeasure.
以上のように勾配磁場をスイッチングさせて空間情報
を含んだ磁気共鳴信号を取得する磁気共鳴診断装置で
は、渦電流磁場強度の空間分布計測が必要不可欠であ
る。従来、この計測にはピックアップコイル出力を時間
積分したり、あるいは小NMRファントムからの磁気共鳴
信号の位相を時間微分する方法が提案されていた。しか
しながら、これらの方法では１ポイント毎に測定して時
定数および強度を求めねばならず、空間分布測定に非常
に手間がかかるという問題がある。この問題解決のため
には、オートスキャナのような特別な治具を新たに作ら
ねばならない。As described above, in the magnetic resonance diagnostic apparatus that obtains a magnetic resonance signal including spatial information by switching the gradient magnetic field, it is essential to measure the spatial distribution of the eddy current magnetic field intensity. Conventionally, for this measurement, a method has been proposed in which the output of the pickup coil is integrated over time or the phase of the magnetic resonance signal from the small NMR phantom is differentiated over time. However, in these methods, the time constant and the intensity have to be obtained by measuring each point, and there is a problem that the spatial distribution measurement is very troublesome. To solve this problem, a special jig such as an auto scanner must be newly created.
また、（ａ）（ｂ）いずれの方法においても渦電流磁
場が残ってしまうため、その残留渦電流磁場の影響を補
正することが不可欠である。この補正方法として、一次
元の磁気共鳴スペクトルデータに対するものは提案され
ている。しかし、この方法は１ポイント毎にスペクトル
データを取得する局所スペクトロスコピーに関してのみ
補正が可能で、空間内多点のスペクトルデータを一時に
取得する多次元化学シフトイメージング法には適用でき
ない。すなわち、多次元化学シフトイメージング法を実
施する場合、スライス勾配磁場やディフェーズ勾配磁場
といった各エンコード（位置情報の位相情報への変換）
毎に強度の変化しない勾配磁場により生じる渦電流磁場
が問題となるが、従来の技術では、このような渦電流磁
場の影響を補正することはできない。Further, since the eddy current magnetic field remains in any of the methods (a) and (b), it is essential to correct the influence of the residual eddy current magnetic field. As this correction method, a method for one-dimensional magnetic resonance spectrum data has been proposed. However, this method can correct only local spectroscopy that acquires spectral data for each point, and cannot be applied to a multidimensional chemical shift imaging method that acquires spectral data at multiple points in space at a time. That is, when performing the multi-dimensional chemical shift imaging method, each encoding (conversion of position information into phase information) such as a slice gradient magnetic field and a dephase gradient magnetic field.
An eddy current magnetic field generated by a gradient magnetic field whose strength does not change every time becomes a problem, but the effect of such an eddy current magnetic field cannot be corrected by the conventional technology.
（発明が解決しようとする課題） 上述したように、従来の技術による渦電流磁場計測方
法では、 測定に非常に手間がかかり、それを避けようとすると
オートスキャナ等の特別な治具を必要とする、 勾配コイル中心から渦電流磁場測定点までの距離精度
を出すために、特別な治具を必要とする、 といった問題があった。(Problems to be Solved by the Invention) As described above, in the eddy current magnetic field measurement method according to the conventional technique, the measurement is extremely troublesome, and a special jig such as an auto scanner is required to avoid the measurement. There is a problem that a special jig is required to obtain the distance accuracy from the gradient coil center to the eddy current magnetic field measurement point.
また、従来の残留渦電流磁場の影響を補正する方法は
局所スペクトロスコピーにのみ有効で、多次元化学シフ
トイメージング法には適用できないという問題があっ
た。Further, the conventional method of correcting the influence of the residual eddy current magnetic field is effective only for local spectroscopy, and has a problem that it cannot be applied to a multidimensional chemical shift imaging method.
本発明はこのような問題点に鑑みてなされたもので、
その目的はオートスキャナのような特別の治具を必要と
しない簡便な構成により、しかも勾配コイル中心からの
距離精度向上のために特別の治具を必要とすることな
く、渦電流磁場の空間分布を計測できる磁気共鳴装置を
提供することにある。The present invention has been made in view of such problems,
Its purpose is a simple configuration that does not require a special jig such as an auto scanner, and also does not require a special jig to improve the distance accuracy from the center of the gradient coil. It is an object of the present invention to provide a magnetic resonance apparatus capable of measuring a magnetic field.
本発明の他の目的は、一次元の磁気共鳴スペクトルデ
ータに対してのみならず、多次元化学シフトイメージン
グ法に対しても渦電流磁場の影響を補正することができ
る磁気共鳴装置を提供することにある。Another object of the present invention is to provide a magnetic resonance apparatus capable of correcting the influence of an eddy current magnetic field not only for one-dimensional magnetic resonance spectrum data but also for multidimensional chemical shift imaging. It is in.
［発明の構成］ （課題を解決するための手段） 本発明に係る磁気共鳴装置は、一様な静磁場中に置か
れた被検体に高周波磁場および勾配磁場を印加すること
により磁気共鳴信号を取得することが可能な磁気共鳴装
置において、渦電流磁場を発生させるための勾配磁場を
印加した後、高周波磁場およびエンコード用勾配磁場を
印加することにより磁気共鳴信号を取得する手段と、こ
の手段により取得された磁気共鳴信号に対して、空間周
波数軸に関するフーリエ変換を施すことにより空間内の
各ポイントにおける磁気共鳴信号を求める手段と、この
手段により求められた空間内の各ポイントにおける磁気
共鳴信号に対して、渦電流磁場の時定数を空間内で固定
した線形最小二乗法を施すことより空間内の各ポイント
毎に渦電流磁場の強度を求める手段とを具備したことを
特徴とする。[Structure of the Invention] (Means for Solving the Problems) A magnetic resonance apparatus according to the present invention generates a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field to a subject placed in a uniform static magnetic field. In a magnetic resonance apparatus capable of acquiring, after applying a gradient magnetic field for generating an eddy current magnetic field, a means for acquiring a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field for encoding, Means for obtaining a magnetic resonance signal at each point in space by applying a Fourier transform to the acquired magnetic resonance signal with respect to the spatial frequency axis; and a means for obtaining a magnetic resonance signal at each point in space obtained by this means. By applying the linear least squares method with the time constant of the eddy current magnetic field fixed in space, the strength of the eddy current magnetic field is obtained for each point in space. Characterized by comprising a that means.
ここで、渦電流磁場成分の強度は、例えば空間内の各
ポイントにおける磁気共鳴信号の位相情報、または空間
内の各ポイントにおける磁気共鳴信号から求められる渦
電流磁場波形に基づいて求められる。Here, the intensity of the eddy current magnetic field component is obtained based on, for example, phase information of a magnetic resonance signal at each point in space or an eddy current magnetic field waveform obtained from a magnetic resonance signal at each point in space.
本発明に係る他の磁気共鳴装置は、一様な静磁場中に
置かれた被検体に高周波磁場および勾配磁場を印加する
ことにより磁気共鳴信号を取得することが可能な磁気共
鳴装置において、渦電流磁場を発生させるための勾配磁
場を印加し、さらに高周波磁場およびエンコード用勾配
磁場を印加した後、勾配磁場を印加しない状態で磁気共
鳴信号を取得する手段と、この手段により取得された磁
気共鳴信号に対して、空間周波数軸に関するフーリエ変
換を施すことにより空間内の各ポイントにおける磁気共
鳴信号を求める手段と、この手段により求められた磁気
共鳴信号に基づいて、空間内の各ポイント毎に渦電流磁
場情報を求める手段とを具備したことを特徴とする。Another magnetic resonance apparatus according to the present invention is a magnetic resonance apparatus capable of acquiring a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field to a subject placed in a uniform static magnetic field. A means for applying a gradient magnetic field for generating a current magnetic field, further applying a high-frequency magnetic field and a gradient magnetic field for encoding, and then acquiring a magnetic resonance signal without applying the gradient magnetic field; and a magnetic resonance signal acquired by the means. Means for obtaining a magnetic resonance signal at each point in space by performing a Fourier transform on the spatial frequency axis on the signal, and eddy for each point in space based on the magnetic resonance signal obtained by this means. Means for obtaining current magnetic field information.
ここで、渦電流磁場情報は、例えば空間内の各ポイン
トにおける磁気共鳴信号の位相情報、または空間内の各
ポイントにおける磁気共鳴信号から求められる解離磁場
波形に基づいて求められる。Here, the eddy current magnetic field information is obtained based on, for example, phase information of a magnetic resonance signal at each point in space or a dissociation magnetic field waveform obtained from a magnetic resonance signal at each point in space.
（作用） このように本発明では、被検体から取得された空間情
報を持たせた磁気共鳴信号に対して、時間軸について固
定の空間周波数軸に関するフーリエ変換を施すことによ
り求められた空間内の各ポイントにおける磁気共鳴信号
から、渦電流磁場各成分の各ポイント毎の強度さらには
時定数が求められる。(Operation) As described above, according to the present invention, the magnetic resonance signal having the spatial information acquired from the subject is subjected to the Fourier transform with respect to the spatial frequency axis fixed on the time axis, thereby obtaining the Fourier transform in the space. From the magnetic resonance signal at each point, the intensity of each component of the eddy current magnetic field component and the time constant are determined.
これにより従来のようにオートスキャナのごとき特別
な治具を必要とすることなく、かつ勾配コイル中心から
渦電流磁場測定点までの距離精度向上のために特別な治
具を必要とせずに、渦電流磁場の空間分布が高精度に求
まる。This eliminates the need for a special jig such as an auto scanner as in the past, and without the need for a special jig to improve the distance accuracy from the gradient coil center to the eddy current magnetic field measurement point. The spatial distribution of the current magnetic field can be determined with high accuracy.
また、本発明では上記と同様にして求められた空間内
の各ポイントにおける磁気共鳴信号に対して、各ポイン
トにおける各エンコード毎に変化しない勾配磁場により
生じるスピンの位相回りの逆数を乗じ、さらに空間軸に
ついて固定の時間軸に関するフーリエ変換を施すこと
で、例えば多次元化学シフトイメージングにおいてスラ
イス用勾配磁場やディフェーズ勾配磁場のようにエンコ
ード毎に変化しないに勾配磁場によって生じる渦電流磁
場の影響が補正された、各ポイントにおける磁気共鳴ス
ペクトルデータが得られる。Further, in the present invention, the magnetic resonance signal at each point in the space obtained in the same manner as described above is multiplied by the reciprocal of the phase around the spin generated by the gradient magnetic field that does not change for each encoding at each point, and By applying a Fourier transform on the fixed time axis, the effect of the eddy current magnetic field caused by the gradient magnetic field is not changed for each encoding such as the slice gradient magnetic field and the dephase gradient magnetic field in multidimensional chemical shift imaging. The obtained magnetic resonance spectrum data at each point is obtained.
（実施例） 以下、図面を参照して本発明の一実施例を説明する。Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
第１図は、本発明の一実施例に係る磁気共鳴装置の構
成を示すブロック図である。FIG. 1 is a block diagram showing a configuration of a magnetic resonance apparatus according to one embodiment of the present invention.
同図において、静磁場磁石１とその内側に設けられた
勾配コイル２およびシムコイル４により、図示しない被
検体に一様な静磁場とそれと同一方向で互いに直交する
x,y,z三方向に線形傾斜磁場分布を持つ勾配磁場が印加
される。勾配コイル２は勾配コイル電源５により駆動さ
れ、シムコイル４はシムコイル電源６により駆動され
る。勾配コイル２の内側に設けられたプローブ３は、送
信部７から高周波信号が供給されることによって被検体
に高周波磁場を印加し、被検体からの磁気共鳴信号を受
信する。プローブ３は送受両用でも、送受別々に設けて
もよい。プローブ３で受信された磁気共鳴信号は受信部
８で検波された後、データ収集部９に転送され、ここで
A/D変換されてから計算機システム10に送られ、データ
処理がなされる。In the figure, a static magnetic field and a gradient coil 2 and a shim coil 4 provided inside the static magnetic field magnet 1 make a subject (not shown) uniform in a uniform static magnetic field and orthogonal to each other in the same direction.
A gradient magnetic field having a linear gradient magnetic field distribution in three directions of x, y, and z is applied. The gradient coil 2 is driven by a gradient coil power supply 5, and the shim coil 4 is driven by a shim coil power supply 6. The probe 3 provided inside the gradient coil 2 applies a high-frequency magnetic field to the subject by receiving a high-frequency signal from the transmission unit 7, and receives a magnetic resonance signal from the subject. The probe 3 may be used for both transmission and reception, or may be provided separately for transmission and reception. The magnetic resonance signal received by the probe 3 is detected by the receiving unit 8 and then transferred to the data collecting unit 9 where it is detected.
After the A / D conversion, the data is sent to the computer system 10 where the data is processed.
以上の勾配コイル用電源５、シムコイル用電源６、送
信部７、受信部８およびデータ収集部９は、全てシーケ
ンス制御部12によって制御され、またシーケンス制御部
12は計算機システム10によって制御される。計算機シス
テム10はコンソール11からの指令により制御される。デ
ータ収集部９から計算機システム10に入力された磁気共
鳴信号は、フーリエ変換、位相計算等が行われ、それに
基づいて被検体内の所望原子核の密度分布、渦電流磁場
分布等の画像データが再構成される。この画像データは
画像ディスプレイ13に送られ、画像として表示される。The power supply 5 for the gradient coil, the power supply 6 for the shim coil, the transmission unit 7, the reception unit 8, and the data collection unit 9 are all controlled by the sequence control unit 12, and the sequence control unit
12 is controlled by the computer system 10. The computer system 10 is controlled by a command from the console 11. The magnetic resonance signal input from the data acquisition unit 9 to the computer system 10 is subjected to Fourier transform, phase calculation, and the like, and based on this, image data such as the density distribution of desired nuclei in the subject and the eddy current magnetic field distribution are regenerated. Be composed. This image data is sent to the image display 13 and displayed as an image.
次に、第２図に示すフローチャートを参照して、本実
施例における渦電流磁場の空間分布を求める手順につい
て説明する。Next, a procedure for obtaining the spatial distribution of the eddy current magnetic field in the present embodiment will be described with reference to the flowchart shown in FIG.
渦電流磁場の空間分布を求める場合、まずステップ20
において予め渦電流磁場各成分の時定数が求められる。
但し、この段階では渦電流磁場各成分の強度は求める必
要はない。また、後に第３図および第４図で説明する渦
電流磁場波形を求める場合は、この処理の流れの中で時
定数を求めることができるので、この渦電流磁場各成分
の時定数を求めるステップは必要ない。When calculating the spatial distribution of the eddy current magnetic field, first, step 20
, The time constant of each component of the eddy current magnetic field is obtained in advance.
However, at this stage, it is not necessary to determine the intensity of each component of the eddy current magnetic field. Further, when an eddy current magnetic field waveform described later with reference to FIGS. 3 and 4 is obtained, a time constant can be obtained in the flow of this processing. Is not required.
渦電流磁場の各成分の時定数を求めるには、第５図
（ａ）（ｂ）に示すようなシーケンスにより磁気共鳴信
号を取得してその位相を計算し、さらに位相を時間微分
することにより実際の渦電流磁場波形を求めてから、最
小二乗法により時定数各成分のみを算出するという、既
に良く知られている方法を用いればよい。但し、第５図
は勾配コイルであるGiコイル（Gi:Gx,Gy,Gz）駆動時の
渦電流磁場各成分の強度を求めるシーケンスである。こ
こで、ガリイ・ハロルド・グローバー氏等の発明（特開
昭62−189057号公報参照）のように、横緩和時間T2 ＊の
減衰が各時定数より短い場合は、第６図（ａ）（ｂ）に
示すように勾配磁場の印加時間をずらし、これから求ま
る渦電流磁場と第５図（ａ）（ｂ）のシーケンスから求
まる渦電流磁場とをつなげることにより、長い期間の渦
電流磁場を求めることができる。また、勾配磁場の印加
時間を更にずらすことにより、より長い期間の渦電流磁
場を求めることができる。渦電流磁場を求める他の方法
として、これも良く知られているようにピックアップコ
イルの出力を時間積分するというものがあるが、本発明
ではこの方法を採用しても全く問題はない。To obtain the time constant of each component of the eddy current magnetic field, a magnetic resonance signal is obtained by a sequence as shown in FIGS. 5 (a) and 5 (b), its phase is calculated, and the phase is differentiated with time. It is sufficient to use a well-known method of calculating only the components of the time constant by the least square method after obtaining the actual eddy current magnetic field waveform. However, FIG. 5 shows a sequence for obtaining the strength of each component of the eddy current magnetic field when driving a Gi coil (Gi: Gx, Gy, Gz) as a gradient coil. Here, when the decay of the transverse relaxation time T 2 * is shorter than each time constant, as in the invention of Galli-Harold Glover (see JP-A-62-189057), FIG. As shown in (b), the application time of the gradient magnetic field is shifted, and the eddy current magnetic field obtained from this is connected to the eddy current magnetic field obtained from the sequence of FIGS. You can ask. Further, by further shifting the application time of the gradient magnetic field, an eddy current magnetic field for a longer period can be obtained. As another method of obtaining the eddy current magnetic field, as is well known, there is a method of integrating the output of the pickup coil with time. However, in the present invention, there is no problem even if this method is employed.
ここで、求められた時定数をτijと定義しておく。添
字のｉは各勾配コイルによる勾配磁場の勾配方向の種
類、即ちｉ＝x,y,zを表し、ｊは各勾配コイル駆動時に
発生する渦電流磁場の各成分を表している。また、以下
で説明する渦電流磁場各成分の強度をAij（ｘ）で定義
しておく。添字の意味は、時定数τijの場合と同様であ
る。Here, the obtained time constant is defined as τij. The subscript i represents the type of gradient direction of the gradient magnetic field by each gradient coil, i.e., i = x, y, z, and j represents each component of the eddy current magnetic field generated when each gradient coil is driven. Further, the intensity of each component of the eddy current magnetic field described below is defined by A ij (x). The meaning of the subscript is the same as that of the time constant τij.
上記のように第２図のステップ20で渦電流磁場各成分
の時定数を求めた後、ステップ21のシーケンス印加によ
り、ステップ22で空間情報を持った磁気共鳴信号Ｓ（t,
k）を取得する。ここで、ｋ＝（kx,ky,kz）は空間周波
数軸であり、以下の次式（１）により定義される。After obtaining the time constants of the respective components of the eddy current magnetic field in step 20 in FIG. 2 as described above, the sequence application in step 21 causes the magnetic resonance signal S (t,
get k). Here, k = (kx, ky, kz) is a spatial frequency axis and is defined by the following equation (1).
ki＝γGiti …（１） ｉ＝x,y,z γ：磁気回転比 Gi:勾配磁場強度 ti:勾配磁場印加時間 次に、ステップ23において、磁気共鳴信号Ｓ（t,k）
に対し、時間ｔについて固定の空間周波数軸ｋに関する
フーリエ変換を施すことにより、空間内の各ポイントの
磁気共鳴信号Ｓ（t,x）を算出する。続いて、ステップ2
4で磁気共鳴信号Ｓ（t,x）の位相Φｉ（t,x）を計算す
る。この位相Φｉ（t,x）は、次式（２）により求める
ことができる。ki = γG i t i (1) i = x, y, z γ: magnetic rotation ratio G i : gradient magnetic field intensity t i : gradient magnetic field application time Next, in step 23, the magnetic resonance signal S (t, k )
By performing a Fourier transform on a fixed spatial frequency axis k with respect to time t, a magnetic resonance signal S (t, x) at each point in space is calculated. Then, step 2
In step 4, the phase Φi (t, x) of the magnetic resonance signal S (t, x) is calculated. This phase Φi (t, x) can be obtained by the following equation (2).
ここで、Re（t,x）,Im（t,x）はそれぞれ磁気共鳴信
号Ｓ（t,k）の実部および虚部である。また、位相Φｉ
（t,x）の時間微分が設定値φｓを越えた時は、式
（２）に２πを加える。 Here, Re (t, x) and Im (t, x) are the real and imaginary parts of the magnetic resonance signal S (t, k), respectively. Also, the phase Φi
When the time derivative of (t, x) exceeds the set value φs, 2π is added to equation (2).
次に、ステップ25で渦電流磁場Ei（t,x）の強度Aij
（ｘ）を算出する。一般に、渦電流磁場Ei（t,x）は次
式（３）により表わせる。Next, in step 25, the intensity Aij of the eddy current magnetic field Ei (t, x)
(X) is calculated. Generally, the eddy current magnetic field Ei (t, x) can be expressed by the following equation (3).
ところで、静磁場磁石１により形成される静磁場には
不均一性ΔＨ（ｘ）があり、磁気共鳴信号の観測時間に
おいて、この磁場不均一性ΔＨ（ｘ）の影響を無視する
ことはできない。磁場不均一性ΔＨ（ｘ）は、ボクセル
内の不均一磁場の平均値で定義される。このとき、時刻
ｔにおける磁気共鳴信号の位相Φｉ（t,x）は次式
（４）により求められる。 By the way, the static magnetic field formed by the static magnetic field magnet 1 has inhomogeneity ΔH (x), and the influence of the magnetic field inhomogeneity ΔH (x) cannot be ignored during the observation time of the magnetic resonance signal. The magnetic field inhomogeneity ΔH (x) is defined by the average value of the inhomogeneous magnetic field in the voxel. At this time, the phase Φi (t, x) of the magnetic resonance signal at the time t is obtained by the following equation (4).
式（４）は、さらに式（５）のように変形される。 Equation (4) is further transformed to equation (5).
式（５）において、渦電流磁場各成分の時定数τijは
前述のように既知であるから、fij（ｔ）を式（６）の
ように定義すると、位相Φｉ（t,x）はさらに式（７）
のように表わせる。よって、線形最小二乗法により渦電
流磁場強度Aij（ｘ）、磁場不均一性ΔＨ（ｘ）を求め
ることができる。 In equation (5), since the time constant τij of each component of the eddy current magnetic field is known as described above, if f ij (t) is defined as in equation (6), the phase Φi (t, x) is further Equation (7)
Can be expressed as Therefore, the eddy current magnetic field strength A ij (x) and the magnetic field inhomogeneity ΔH (x) can be obtained by the linear least square method.
最後に、ステップ26で渦電流磁場各成分の時定数τij
に対して、それぞれ渦電流磁場強度Aij（ｘ）をグレイ
スケール等で表示することにより、第１図の画像ディス
プレイ13上で渦電流磁場の強度分布を画像として表示す
ることができる。 Finally, in step 26, the time constant τ ij of each component of the eddy current magnetic field
On the other hand, by displaying the eddy current magnetic field strength A ij (x) in a gray scale or the like, the eddy current magnetic field strength distribution can be displayed as an image on the image display 13 in FIG.
次に、第２図のステップ21に示した空間情報を持つ磁
気共鳴信号Ｓ（t,x）の取得のためのシーケンスについ
て説明しておく。このシーケンスとしては、第７図に示
す4DCSI−FIDあるいは第８図に示す4DCSI−SEが用いら
れる。なお、第７図および第８図の例では、ｚ方向の勾
配磁場Gzを発生する勾配コイル駆動時に発生する渦電流
磁場計測用のシーケンスを示した。渦電流磁場計測用勾
配磁場に対しエンコード用勾配磁場は非常に小さくとれ
るため、第８図においてこれによる影響は無視できる。Next, the sequence for acquiring the magnetic resonance signal S (t, x) having the spatial information shown in step 21 of FIG. 2 will be described. As this sequence, 4DCSI-FID shown in FIG. 7 or 4DCSI-SE shown in FIG. 8 is used. 7 and 8 show the sequence for measuring the eddy current magnetic field generated at the time of driving the gradient coil for generating the gradient magnetic field Gz in the z direction. Since the encoding gradient magnetic field can be made very small compared to the eddy current magnetic field measuring gradient magnetic field, the influence of this can be ignored in FIG.
また、勾配磁場の立上がり、立下がりを考慮しなけれ
ばならないので、実際に観測される渦電流磁場強度ij
（t,x）は、FIDの場合は式（８）のようになり、またSE
の場合は式（９）のようになる。但し、SEの場合は勾配
磁場強度EGiBにおける渦電流磁場強度をij（ｘ）と定
義した。Also, since the rise and fall of the gradient magnetic field must be considered, the eddy current magnetic field strength ij
(T, x) is as shown in equation (8) for FID, and SE
Equation (9) is obtained in the case of. However, in the case of SE, the eddy current magnetic field strength at the gradient magnetic field strength E GiB was defined as ij (x).
また、FIDにおいてはRFパルスおよびエンコード用勾
配磁場によるデッドタイムがあるため、式（８）は式
（10）のように変形される。 Further, in the FID, since there is a dead time due to the RF pulse and the gradient magnetic field for encoding, Expression (8) is transformed into Expression (10).
従って、得られた渦電流磁場強度Aij（ｘ）をそれぞ
れ式（11）、式（12）により補正して、真の渦電流磁場
強度Aij（ｘ）を求める必要がある。 Therefore, it is necessary to correct the obtained eddy current magnetic field strength A ij (x) by the equations (11) and (12) to obtain the true eddy current magnetic field strength A ij (x).
また、シーケンス設定時に各々のAij（ｘ）が小さく
なり過ぎないように、勾配磁場強度およびパルス幅を調
整してやらなければならない。さらに、第８図のSEの場
合、渦電流磁場計測用の勾配磁場強度およびパルス幅
は、次式（13）を満たすように設定する必要があること
を付加えておく。 Also, the gradient magnetic field strength and the pulse width must be adjusted so that each A ij (x) does not become too small when setting the sequence. Furthermore, in the case of SE in FIG. 8, it is added that the gradient magnetic field intensity and pulse width for eddy current magnetic field measurement need to be set so as to satisfy the following expression (13).
EGiA（taGi−ta）＝EGiB（tbGi−tb） …（13） 第３図および第４図は本発明の他の実施例における渦
電流磁場の空間分布を求める手順を示すフローチャート
であり、T2 ＊が時定数τijと比較して長い場合に適した
例である。第３図におけるステップ31〜34は第２図にお
けるステップ21〜24と同じであり、ステップ34において
位相Φ（t,x）を算出した後、各ポイントの渦電流磁場
波形Ei（t,x）を算出し、これから渦電流磁場各成分の
時定数τij（ｘ）および強度Aij（ｘ）を求め、画像と
して表示する（ステップ35〜37）。E GiA (ta Gi −ta) = E GiB (tb Gi −tb) (13) FIGS. 3 and 4 are flowcharts showing a procedure for obtaining the spatial distribution of the eddy current magnetic field in another embodiment of the present invention. This is an example suitable for the case where T 2 * is longer than the time constant τij. Steps 31 to 34 in FIG. 3 are the same as steps 21 to 24 in FIG. 2. After calculating the phase Φ (t, x) in step 34, the eddy current magnetic field waveform Ei (t, x) at each point is calculated. Is calculated, and the time constant τ ij (x) and the intensity A ij (x) of each component of the eddy current magnetic field are obtained from the calculated values and displayed as images (steps 35 to 37).
第４図におけるステップ41〜43は、第２図におけるス
テップ21〜23と同じであり、ステップ43において各ポイ
ントにおける磁気共鳴信号Ｓ（t,x）を算出した後、各
ポイントの渦電流磁場波形Ei（t,x）を算出し、これか
ら渦電流磁場各成分の時定数τij（ｘ）および強度Aij
（ｘ）を求め、画像として表示する（ステップ44〜4
6）。Steps 41 to 43 in FIG. 4 are the same as steps 21 to 23 in FIG. 2. After calculating the magnetic resonance signal S (t, x) at each point in step 43, the eddy current magnetic field waveform at each point is calculated. Ei (t, x) is calculated, and from this, the time constant τ ij (x) and the intensity A ij of each component of the eddy current magnetic field are calculated.
(X) is obtained and displayed as an image (steps 44 to 4).
6).
これら第３図および第４図の場合は、渦電流磁場の時
定数τij（ｘ）を予め求めておくための第２図に示した
ステエップ20は必要ない。第３図のステップ35において
は、次式（14）のように位相Φｉ（t,x）を時間微分す
ることにより渦電流磁場波形Ei（t,x）が得られ、また
第４図のステップ44においては、磁気共鳴信号Ｓ（t,
x）の実部Re（t,x）、虚部Im（t,x）から、式（15）に
より渦電流磁場波形Ei（t,x）が得られる。In the case of FIGS. 3 and 4, the step 20 shown in FIG. 2 for obtaining the time constant τ ij (x) of the eddy current magnetic field in advance is not necessary. In step 35 of FIG. 3, an eddy current magnetic field waveform Ei (t, x) is obtained by time-differentiating the phase Φi (t, x) as in the following equation (14). At 44, the magnetic resonance signal S (t,
From the real part Re (t, x) and the imaginary part Im (t, x) of x), an eddy current magnetic field waveform Ei (t, x) is obtained by Expression (15).
次に、ステップ36,45で渦電流磁場波形Ei（t,x）に最
小二乗法を適用することにより、時定数τij（ｘ）と強
度Aij（ｘ）を求めることができる。なお、ステップ36,
45の他の方法として数ポイントの渦電流磁場波形から最
小二乗法により時定数τij（ｘ）を求め、誤差を取除く
ために平均等の処理を起った後に、線形最小二乗法によ
り強度Aij（ｘ）を求めることもできる。 Next, in steps 36 and 45, the time constant τ ij (x) and the intensity A ij (x) can be obtained by applying the least square method to the eddy current magnetic field waveform Ei (t, x). Step 36,
45 As another method, the time constant τ ij (x) is obtained from the eddy current magnetic field waveforms of several points by the least square method, and processing such as averaging is performed to remove errors, and the intensity is then obtained by the linear least square method. A ij (x) can also be determined.
次に、ステップ37,46では第２図のステップ26と同様
に、渦電流磁場強度Aij（ｘ）の画像化を行うと同時
に、同じ手法で渦電流磁場時定数τij（ｘ）の空間分布
の画像化を行う。Next, in steps 37 and 46, similarly to step 26 in FIG. 2, imaging of the eddy current magnetic field strength A ij (x) is performed, and at the same time, the space of the eddy current magnetic field time constant τ ij (x) is Image the distribution.
また、T2 ＊が時定数τijに比べて短い場合でも、第９
図や第10図に示したように、渦電流磁場計測用の勾配磁
場の印加時間をずらして取得した磁気共鳴信号から求ま
る渦電流磁場を、勾配磁場の印加時間をずらさずに取得
した磁気共鳴信号から求まる渦電流磁場とつなげること
により、長い期間の渦電流磁場を求めることができ、印
加時間のずらし幅を変えることにより、より長い期間の
渦電流磁場を求めることができる。Even when T 2 * is shorter than the time constant τ ij , the ninth
As shown in Fig. 10 and Fig. 10, the eddy current magnetic field obtained from the magnetic resonance signal obtained by shifting the application time of the gradient magnetic field for measuring the eddy current magnetic field was obtained without changing the application time of the gradient magnetic field. By connecting to the eddy current magnetic field obtained from the signal, an eddy current magnetic field for a long period can be obtained. By changing the shift width of the application time, an eddy current magnetic field for a longer period can be obtained.
また、第２図の方法においてステップ25で位相Φｉ
（t,x）の時間微分 から渦電流磁場強度Aij（ｘ）を求めたり、第３図の方
法において位相Φｉ（t,x）から時定数τij（ｘ）と強
度Aij（ｘ）を求めたり、あるいは第３図の方法で位相
Φｉ（t,x）を使って時定数固定展開で強度Aij（ｘ）を
求める、といった上記方法を組合わせた方法で行っても
上記実施例と同様の効果が得られる。In addition, in the method of FIG.
Time derivative of (t, x) Asking for an eddy current field strength A ij (x) from, or seek constant τ ij (x) and the intensity A ij (x) when the phase Φi (t, x) in the method of FIG. 3 or FIG. 3, The same effect as in the above embodiment can be obtained by a method combining the above methods, such as finding the intensity A ij (x) by time constant fixed expansion using the phase Φi (t, x) in the method described above.
なお、以上の実施例では渦電流磁場計測についてのみ
説明したが、同様の方法により勾配コイルとシムコイル
のとカップリング磁場の計測が可能であり、またシーケ
ンスから渦電流磁場計測用の勾配磁場を除くことによ
り、磁気共鳴信号を取得する空間内各ポイントでの変動
磁場を計測することも可能である。Although only the eddy current magnetic field measurement has been described in the above embodiment, it is possible to measure the coupling magnetic field with the gradient coil and the shim coil by the same method, and to exclude the gradient magnetic field for the eddy current magnetic field measurement from the sequence. This makes it possible to measure the fluctuating magnetic field at each point in the space where the magnetic resonance signal is acquired.
次に、第11図に示すフローチャートを参照して本実施
例における渦電流磁場の影響を補正する手順について説
明する。ここでは、各エンコード毎に変化しない勾配磁
場を使用したシーケンスの一例として、スライス勾配磁
場を使用した３次元化学シフトイメージングのシーケン
ス（第12図参照）を例にとり、スライス勾配磁場により
発生する渦電流磁場の影響を補正する方法について述べ
る。なお、エンコード勾配磁場により生じる渦電流磁場
は非常に小さいため、一般には補正の必要はない。Next, a procedure for correcting the influence of the eddy current magnetic field in the present embodiment will be described with reference to the flowchart shown in FIG. Here, as an example of a sequence using a gradient magnetic field that does not change for each encoding, a sequence of three-dimensional chemical shift imaging using a slice gradient magnetic field (see FIG. 12) is taken as an example, and an eddy current generated by the slice gradient magnetic field is taken as an example. A method for correcting the influence of the magnetic field will be described. Since the eddy current magnetic field generated by the encoding gradient magnetic field is very small, there is generally no need for correction.
まず、第11図のステップ41においてパルスシーケンス
の印加により磁気共鳴信号を取得する。密度分布をρ
（x,ω）、スライス勾配磁場により発生する各ポイント
の渦電流磁場によるスピンの位相回りを により定義した時、磁気共鳴信号ｓ（t,k）は次式（1
6）により定義される。First, in step 41 of FIG. 11, a magnetic resonance signal is obtained by applying a pulse sequence. Let the density distribution be ρ
(X, ω), the phase around the spin caused by the eddy current magnetic field at each point generated by the slice gradient magnetic field When the magnetic resonance signal s (t, k) is defined by the following equation,
6).
次に、ステップ52においてｓ（t,k）に対して、時間
について固定で、空間周波数軸ｋに関するフーリエ変換
を施すことにより、各ポイントの磁気共鳴信号Ｓ（t,
x）を算出する。 Next, in step 52, the magnetic resonance signal S (t, k) at each point is subjected to a Fourier transform on s (t, k) with respect to the spatial frequency axis k, which is fixed with respect to time.
x) is calculated.
次に、ステップ53において、予め求めておいた、スラ
イス勾配磁場により発生する渦電流磁場によるスピンの
位相回りの逆数、つまり をステップ52で得られた磁気共鳴信号Ｓ（t,x）に各ポ
イント毎に乗じることによって、渦電流磁場による影響
を補正した、次式（18）に示す磁気共鳴信号Ｓ（t,x）
を求める。 Next, in step 53, the reciprocal of the phase around the spin due to the eddy current magnetic field generated by the slice gradient magnetic field, which is obtained in advance, that is, Is multiplied for each point by the magnetic resonance signal S (t, x) obtained in step 52, thereby correcting the influence of the eddy current magnetic field. The magnetic resonance signal S (t, x) shown in the following equation (18)
Ask for.
最後に、ステップ53で得られた磁気共鳴信号Ｓ（t,
x）に対して、ステップ54で時間ｔに関してフーリエ変
換を施すことにより、渦電流磁場による影響が補正され
た空間多点の磁気共鳴スペクトルＳ（ω,x）が求められ
る。 Finally, the magnetic resonance signal S (t,
By performing a Fourier transform on time t in step 54 with respect to x), a magnetic resonance spectrum S (ω, x) at multiple points in space in which the influence of the eddy current magnetic field is corrected is obtained.
次に、第11図におけるステップ53で用いる位相回りの
求め方について説明する。第13図は、この位相回りを求
める手順を示すフローチャートである。Next, a method of obtaining the phase rotation used in step 53 in FIG. 11 will be described. FIG. 13 is a flowchart showing a procedure for obtaining the phase rotation.
まず、ステップ61において先と同様にして渦電流磁場
の空間分布を計測し、それに基づいて時定数τij（ｘ）
および強度Aij（ｘ）を算出する。この際、計測される
渦電流磁場の空間分布は離散的に求められるものである
ため、ステップ62で渦電流磁場の空間分布を連続関数に
変換し、時定数τij（ｘ）および強度Aij（ｘ）を空間
座標ｘの連続関数に変換する。この変換には、例えば最
小二乗法を用いてルジャンドル関数等の直交関数展開な
どの方法を用いればよい。First, in step 61, the spatial distribution of the eddy current magnetic field is measured in the same manner as described above, and the time constant τ ij (x)
And intensity A ij (x). At this time, since the measured spatial distribution of the eddy current magnetic field is discretely obtained, the spatial distribution of the eddy current magnetic field is converted into a continuous function in step 62, and the time constant τ ij (x) and the intensity A ij (X) is converted into a continuous function of the spatial coordinates x. For this conversion, for example, a method such as expansion of orthogonal functions such as Legendre functions using the least squares method may be used.
次に、ステップ63でスペクトル計測ボクセルの代表値
における位相φed（t,x）を計算する。第12図に示した
３次元化学シフトイメージングにおけるスライス勾配磁
場により生じる渦電流磁場について述べる。勾配磁場G0
に対して渦電流磁場Azj（ｘ）が生じる時、第12図のパ
ルスシーケンスにおいて磁気共鳴信号の取得時に混入す
る渦電流磁場の大きさBed（t,x）は、次式（19）で表さ
れる。Next, in step 63, the phase φ ed (t, x) at the representative value of the spectrum measurement voxel is calculated. The eddy current magnetic field generated by the slice gradient magnetic field in the three-dimensional chemical shift imaging shown in FIG. 12 will be described. Gradient magnetic field G 0
When the eddy current magnetic field A zj (x) is generated, the magnitude B ed (t, x) of the eddy current magnetic field mixed at the time of acquiring the magnetic resonance signal in the pulse sequence of FIG. It is represented by
よって、位相φed（t,x）は次式（20）で表わされ
る。 Therefore, the phase φ ed (t, x) is expressed by the following equation (20).
次に、ステップ64で位相φed（t,x）から、渦電流磁
場による位相回り を算出する。 Next, in step 64, from the phase φ ed (t, x), the phase Is calculated.
なお、以上の実施例ではスライス勾配磁場により発生
する渦電流磁場の影響を補正する例について述べたが、
他の各エンコード毎に変化しない勾配磁場、例えばディ
フェーズ用の勾配磁場によって発生する渦電流磁場の影
響を補正する場合にも本発明を適用できることはいうま
でもない。In the above embodiment, the example in which the influence of the eddy current magnetic field generated by the slice gradient magnetic field is corrected has been described.
It is needless to say that the present invention can be applied to a case where the influence of an eddy current magnetic field generated by a gradient magnetic field which does not change for each other encoding, for example, a gradient magnetic field for dephase is corrected.
［発明の効果］ 本発明によれば、高周波磁場およびエンコード用勾配
磁場を印加することにより空間情報を持たせた磁気共鳴
信号を取得し、これに時間軸固定の空間周波数軸に関す
るフーリエ変換を施して空間内各ポイントにおける磁気
共鳴信号を求め、この磁気共鳴信号を用いて渦電流磁場
の各ポイント毎の強度やさらに時定数を求めることによ
り、オートスキャナ等の特別な治具を必要とせずに、し
かも勾配コイル中心から渦電流磁場測定点までの距離精
度向上のために特別な治具を必要とせずに、渦電流磁場
の空間分布計測を行うことができる。[Effects of the Invention] According to the present invention, a magnetic resonance signal having spatial information is obtained by applying a high-frequency magnetic field and a gradient magnetic field for encoding, and a Fourier transform is performed on the time-fixed spatial frequency axis. By obtaining the magnetic resonance signal at each point in the space, and using this magnetic resonance signal to obtain the intensity and time constant of each point of the eddy current magnetic field, there is no need for a special jig such as an auto scanner. In addition, the spatial distribution of the eddy current magnetic field can be measured without requiring a special jig to improve the distance accuracy from the gradient coil center to the eddy current magnetic field measurement point.
また、本発明によれば各ポイントにおける磁気共鳴信
号に対して、各ポイントにおける各エンコード毎に変化
しない勾配磁場により生じるスピンの位相回りの逆数を
乗じ、さらに空間軸について固定の時間軸に関するフー
リエ変換を施すことで、一次元スペクトルデータに対し
てのみでなく、例えば多次元化学シフトイメージングに
対しても、スライス用勾配磁場やディフェーズ勾配磁場
といったエンコード毎に変化しないに勾配磁場によって
生じる渦電流磁場の影響を補正することが可能となる。According to the present invention, the magnetic resonance signal at each point is multiplied by the reciprocal of the phase around the spin generated by the gradient magnetic field that does not change for each encoding at each point, and the Fourier transform on the time axis fixed on the spatial axis is further performed. The eddy current magnetic field generated by the gradient magnetic field that does not change for each encoding, such as the slice gradient magnetic field and the dephase gradient magnetic field, not only for one-dimensional spectral data but also for multidimensional chemical shift imaging, for example, Can be corrected.
第１図は本発明の一実施例に係る磁気共鳴装置の構成を
示すブロック図、第２図、第３図および第４図は本発明
の実施例における渦電流磁場計測の流れを示すフローチ
ャート、第５図および第６図は渦電流磁場の時定数を磁
気共鳴信号から求めるためのシーケンスを示す図、第７
図および第８図は渦電流磁場の強度の空間分布を求める
ためのシーケンスを示す図、第９図および第10図はT2 ＊
が短い場合に渦電流磁場波形を空間的に求めるためのシ
ーケンスを示す図、第11図は同実施例における渦電流磁
場による影響を補正するための手順を示すフローチャー
ト、第12図はスライス勾配磁場を使用した３次元化学シ
フトイメージングのシーケンスを示す図、第13図は同実
施例における位相回りを求める手順を示すフローチャー
トである。 １……静磁場磁石、２……勾配コイル、３……プロー
ブ、４……シムコイル、５……勾配コイル、６……シム
コイル電源、７……送信部、８……受信部、９……デー
タ収集部、10……計算機システム、11……コンソール、
12……シーケンス制御部、13……画像ディスプレイ。FIG. 1 is a block diagram showing a configuration of a magnetic resonance apparatus according to one embodiment of the present invention, FIGS. 2, 3 and 4 are flowcharts showing a flow of eddy current magnetic field measurement in an embodiment of the present invention, 5 and 6 show a sequence for obtaining the time constant of the eddy current magnetic field from the magnetic resonance signal.
FIG. 8 and FIG. 8 show a sequence for obtaining the spatial distribution of the intensity of the eddy current magnetic field, and FIGS. 9 and 10 show T 2 *.
FIG. 11 is a diagram showing a sequence for spatially obtaining an eddy current magnetic field waveform in the case where is short, FIG. 11 is a flowchart showing a procedure for correcting the influence of the eddy current magnetic field in the embodiment, and FIG. 12 is a slice gradient magnetic field. FIG. 13 is a diagram showing a sequence of three-dimensional chemical shift imaging using, and FIG. 13 is a flowchart showing a procedure for obtaining a phase rotation in the embodiment. 1 ... static magnetic field magnet, 2 ... gradient coil, 3 ... probe, 4 ... shim coil, 5 ... gradient coil, 6 ... shim coil power supply, 7 ... transmission unit, 8 ... reception unit, 9 ... Data collection unit, 10 Computer system, 11 Console
12: Sequence control unit, 13: Image display.
───────────────────────────────────────────────────── フロントページの続き (72)発明者 守 清巳 神奈川県川崎市幸区小向東芝町１番地 株式会社東芝総合研究所内 (72)発明者 石原 康利 神奈川県川崎市幸区小向東芝町１番地 株式会社東芝総合研究所内 (72)発明者 佐藤 幸三 神奈川県川崎市幸区小向東芝町１番地 株式会社東芝総合研究所内 (56)参考文献 特開 昭62−189057（ＪＰ，Ａ) 特開 平３−16554（ＪＰ，Ａ) 米国特許4910460（ＵＳ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 5/055 ＪＩＣＳＴファイル（ＪＯＩＳ)──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Kiyomi Mori 1 Toshiba-cho, Komukai, Saiwai-ku, Kawasaki-shi, Kanagawa Prefecture (72) Inventor Yasutoshi Ishihara 1 Toshiba-cho, Komukai-shiku, Kawasaki-shi, Kanagawa Address: Toshiba Research Institute Co., Ltd. (72) Kozo Sato, Inventor: Kobayashi-ku, Kawasaki City, Kanagawa Prefecture 1, Komukai Toshiba Town, Toshiba Research Institute Co., Ltd. (56) References JP 62-188907 (JP, A) JP Hei 3-16554 (JP, A) U.S. Pat. No. 4,910,460 (US, A) (58) Fields investigated (Int. Cl. 7 , DB name) A61B 5/055 JICST file (JOIS)
Claims (6)
磁場および勾配磁場を印加することにより磁気共鳴信号
を取得することが可能な磁気共鳴装置において、 渦電流磁場を発生させるための勾配磁場を印加した後、
高周波磁場およびエンコード用勾配磁場を印加すること
により磁気共鳴信号を取得する手段と、 この手段により取得された磁気共鳴信号に対して、空間
周波数軸に関するフーリエ変換を施すことにより空間内
の各ポイントにおける磁気共鳴信号を求める手段と、 この手段により求められた空間内の各ポイントにおける
磁気共鳴信号に対して、前記渦電流磁場の時定数を空間
内で固定した線形最小二乗法を施すことより空間内の各
ポイント毎に前記渦電流磁場の強度を求める手段と を具備したことを特徴とする磁気共鳴装置。1. A magnetic resonance apparatus capable of acquiring a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field to a subject placed in a uniform static magnetic field, for generating an eddy current magnetic field. After applying the gradient magnetic field of
Means for acquiring a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field for encoding; and applying a Fourier transform to a spatial frequency axis to the magnetic resonance signal acquired by this means, thereby obtaining a magnetic resonance signal at each point in space. Means for obtaining a magnetic resonance signal; and applying a linear least squares method in which the time constant of the eddy current magnetic field is fixed in space to the magnetic resonance signal at each point in space obtained by this means. Means for calculating the intensity of the eddy current magnetic field at each point of the magnetic resonance apparatus.
は、前記空間内の各ポイントにおける磁気共鳴信号の位
相情報に基づいて該渦電流磁場成分の強度を求めるもの
であることを特徴とする請求項１記載の磁気共鳴装置。2. The method according to claim 1, wherein said means for determining the intensity of the eddy current magnetic field component is for obtaining the intensity of the eddy current magnetic field component based on phase information of a magnetic resonance signal at each point in the space. The magnetic resonance apparatus according to claim 1.
は、前記空間内の各ポイントにおける磁気共鳴信号から
求められる渦電流磁場波形に基づいて該渦電流磁場成分
の強度を求めるものであることを特徴とする請求項１記
載の磁気共鳴装置。3. The means for determining the intensity of an eddy current magnetic field component is based on an eddy current magnetic field component obtained from a magnetic resonance signal at each point in the space. The magnetic resonance apparatus according to claim 1, wherein:
磁場および勾配磁場を印加することにより磁気共鳴信号
を取得することが可能な磁気共鳴装置において、 渦電流磁場を発生させるための勾配磁場を印加し、さら
に高周波磁場およびエンコード用勾配磁場を印加した
後、勾配磁場を印加しない状態で磁気共鳴信号を取得す
る手段と、 この手段により取得された磁気共鳴信号に対して、空間
周波数軸に関するフーリエ変換を施すことにより空間内
の各ポイントにおける磁気共鳴信号を求める手段と、 この手段により求められた磁気共鳴信号に基づいて、空
間内の各ポイント毎に前記渦電流磁場情報を求める手段
と を具備したことを特徴とする磁気共鳴装置。4. A magnetic resonance apparatus capable of acquiring a magnetic resonance signal by applying a high-frequency magnetic field and a gradient magnetic field to a subject placed in a uniform static magnetic field. Means for applying a gradient magnetic field, applying a high-frequency magnetic field and a gradient magnetic field for encoding, and then acquiring a magnetic resonance signal without applying the gradient magnetic field; Means for obtaining a magnetic resonance signal at each point in space by performing a Fourier transform on a frequency axis; and obtaining the eddy current magnetic field information for each point in space based on the magnetic resonance signal obtained by this means. A magnetic resonance apparatus comprising:
空間内の各ポイントにおける磁気共鳴信号の位相情報に
基づいて該渦電流磁場情報を求めるものであることを特
徴とする請求項４記載の磁気共鳴装置。5. The apparatus according to claim 4, wherein said means for obtaining eddy current magnetic field information obtains the eddy current magnetic field information based on phase information of a magnetic resonance signal at each point in the space. Magnetic resonance equipment.
空間内の各ポイントにおける磁気共鳴信号から求められ
る渦電流磁場波形に基づいて該渦電流磁場情報を求める
ものであることを特徴とする請求項４記載の磁気共鳴装
置。6. The eddy current magnetic field information obtaining means obtains the eddy current magnetic field information based on an eddy current magnetic field waveform obtained from a magnetic resonance signal at each point in the space. The magnetic resonance apparatus according to claim 4.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2-259065 | 1990-09-28 | ||
| JP25906590 | 1990-09-28 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04189344A JPH04189344A (en) | 1992-07-07 |
| JP3112930B2 true JP3112930B2 (en) | 2000-11-27 |
Family
ID=17328837
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP02340398A Expired - Fee Related JP3112930B2 (en) | 1990-09-28 | 1990-11-30 | Magnetic resonance equipment |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3112930B2 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5883207B2 (en) * | 2007-12-28 | 2016-03-09 | 株式会社東芝 | Magnetic resonance imaging apparatus and control program for magnetic resonance imaging apparatus |
| JP5121773B2 (en) * | 2009-05-11 | 2013-01-16 | 株式会社東芝 | Magnetic resonance imaging system |
| JP2012183233A (en) | 2011-03-07 | 2012-09-27 | Toshiba Corp | Magnetic resonance imaging apparatus |
| JP5971683B2 (en) * | 2011-03-15 | 2016-08-17 | 東芝メディカルシステムズ株式会社 | Magnetic resonance imaging system |
-
1990
- 1990-11-30 JP JP02340398A patent/JP3112930B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH04189344A (en) | 1992-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5151656A (en) | Correction of nmr data acquired by an echo-planar technique | |
| Morgan et al. | Correction of spatial distortion in EPI due to inhomogeneous static magnetic fields using the reversed gradient method | |
| US4698591A (en) | Method for magnetic field gradient eddy current compensation | |
| US5351006A (en) | Method and apparatus for correcting spatial distortion in magnetic resonance images due to magnetic field inhomogeneity including inhomogeneity due to susceptibility variations | |
| KR100481740B1 (en) | A method for measuring and compensating for spatially and temporally varying magnetic fields induced by eddy currents | |
| DE19804823B4 (en) | Correct artifacts caused by Maxwell terms in Magnetic Resonance Echo Planar images | |
| JP2009034514A (en) | Determination method of imaging parameter, optimization device of imaging parameter, computer program product and electronically readable data medium | |
| US11585884B2 (en) | Continual trajectory correction in magnetic resonance imaging | |
| DE19631915A1 (en) | Magnetic resonance data acquisition from two sites during single scan | |
| US6933720B2 (en) | Sequence preconditioning for ultra-fast magnetic resonance imaging | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| JPH0244220B2 (en) | ||
| US5079504A (en) | Magnetic resonance imaging system | |
| US6157192A (en) | Recovery of signal void arising from field inhomogeneities in echo planar imaging | |
| JPS62254744A (en) | Method and apparatus for measuring nuclear magnetism distribution of part of matter | |
| JPH0315455B2 (en) | ||
| US6636756B2 (en) | K-space synthesis for MR imaging in the presence of gradient field nonlinearity | |
| JP3112930B2 (en) | Magnetic resonance equipment | |
| JP3538249B2 (en) | MRIS device | |
| JP2805405B2 (en) | Magnetic resonance imaging equipment | |
| US4746865A (en) | Method and system for magnetic resonance imaging | |
| US4724387A (en) | Method and apparatus for high resolution nuclear magnetic resonance imaging | |
| JPH01141656A (en) | Method for correcting magnetic field drift of nmr imaging apparatus | |
| JP3907944B2 (en) | Magnetic resonance imaging method and apparatus | |
| JP2607466B2 (en) | Inspection equipment using nuclear magnetic resonance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |