JP4625677B2 - Magnetic resonance imaging apparatus and image correction evaluation method - Google Patents
Magnetic resonance imaging apparatus and image correction evaluation method Download PDFInfo
- Publication number
- JP4625677B2 JP4625677B2 JP2004309971A JP2004309971A JP4625677B2 JP 4625677 B2 JP4625677 B2 JP 4625677B2 JP 2004309971 A JP2004309971 A JP 2004309971A JP 2004309971 A JP2004309971 A JP 2004309971A JP 4625677 B2 JP4625677 B2 JP 4625677B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- image
- image data
- correction
- static magnetic
- 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
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明は、被検体の原子核スピンをラーモア周波数のＲＦ信号で磁気的に励起し、この励起に伴って発生する磁気共鳴信号から画像を再構成する磁気共鳴イメージング装置および画像補正評価方法に係り、特に磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することが可能な磁気共鳴イメージング装置および画像補正評価方法に関する。 The present invention relates to a magnetic resonance imaging apparatus and an image correction evaluation method for magnetically exciting a nuclear spin of a subject with an RF signal having a Larmor frequency and reconstructing an image from a magnetic resonance signal generated along with the excitation. In particular, the present invention relates to a magnetic resonance imaging apparatus and an image correction evaluation method capable of evaluating whether or not image distortion correction due to magnetic field inhomogeneity has been appropriately performed.
従来、医療現場におけるモニタリング装置として、磁気共鳴イメージング（ＭＲＩ：Magnetic Resonance Imaging）装置が利用される。 Conventionally, a magnetic resonance imaging (MRI) apparatus is used as a monitoring apparatus in a medical field.
ＭＲＩ装置は、静磁場を形成する筒状の静磁場用磁石内部にセットされた被検体の撮像領域に傾斜磁場コイルでＸ軸、Ｙ軸、Ｚ軸方向の傾斜磁場を形成するとともにＲＦ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ）コイルから高周波（ＲＦ）信号を送信することにより被検体内の原子核スピンを磁気的に共鳴させ、励起により生じた核磁気共鳴（ＮＭＲ：Nuclear Magnetic Resonance）信号を利用して被検体の画像を再構成する装置である。 The MRI apparatus forms gradient magnetic fields in the X-axis, Y-axis, and Z-axis directions in an imaging region of a subject set inside a cylindrical static magnetic field magnet that forms a static magnetic field by using a gradient magnetic field coil and RF (Radio). By transmitting a radio frequency (RF) signal from a frequency coil, the nuclear spin in the subject is magnetically resonated, and an image of the subject is obtained using a nuclear magnetic resonance (NMR) signal generated by excitation. Is a device for reconfiguring.
このようなＭＲＩ装置では、静磁場を形成する静磁場用磁石の性能限界や磁場中に被検体が入ることによる磁化率の変動によって、局所的に静磁場強度の均一性が乱されるという問題がある。静磁場強度の均一性が乱されると、再構成された画像に歪みが生じることが知られている。 In such an MRI apparatus, there is a problem that the uniformity of the static magnetic field strength is locally disturbed due to the performance limit of the static magnetic field magnet that forms the static magnetic field and the change in magnetic susceptibility due to the subject entering the magnetic field. There is. It is known that when the uniformity of the static magnetic field strength is disturbed, the reconstructed image is distorted.
静磁場強度の不均一性の影響を特に強く受けるケースとしては、ＥＰＩ（echo planar imaging）法による撮像を行う場合が挙げられる。ＥＰＩ法は、一回の励起パルスで画像再構成のためのすべてのデータを取得する超高速スキャン法であり、シングルショットＥＰＩ法とも呼ばれる。ＥＰＩ法による撮像は、例えば被検体の頭部の検査における脳梗塞などの描出の際に、被検体の体動や拍動の影響を避けるために行われるが、ＥＰＩ法によって撮像された拡散強調画像は静磁場強度の不均一性の影響を強く受け、静磁場強度が十分に均一でないと再構成された画像に歪みが生じる。 As a case that is particularly strongly affected by the non-uniformity of the static magnetic field strength, there is a case where imaging is performed by an EPI (echo planar imaging) method. The EPI method is an ultra-high-speed scanning method that acquires all data for image reconstruction with a single excitation pulse, and is also called a single shot EPI method. Imaging by the EPI method is performed in order to avoid the influence of the body motion and pulsation of the subject at the time of rendering a cerebral infarction in the examination of the head of the subject, for example. The image is strongly affected by the non-uniformity of the static magnetic field strength, and the reconstructed image is distorted if the static magnetic field strength is not sufficiently uniform.
そこで、静磁場強度の不均一性による問題を解決する方法として、従来、静磁場強度の分布を推定し、推定して得られた静磁場強度の分布をもとに再構成された画像の位置および画像値を補正するという方法がとられている（例えば非特許文献１参照）。 Therefore, as a method of solving the problem due to the non-uniformity of the static magnetic field strength, the position of the image reconstructed based on the static magnetic field strength distribution obtained by estimating the static magnetic field strength distribution and estimating the distribution of the static magnetic field strength has been conventionally used. And a method of correcting the image value is used (see, for example, Non-Patent Document 1).
静磁場強度の均一性を乱し、再構成された画像に歪みをもたらす要因としては、傾斜磁場の誤差、静磁場用磁石単体による静磁場誤差および被検体が静磁場内に入ることによる静磁場の誤差の３つの要因がある。このため、静磁場の誤差は、傾斜磁場の誤差、静磁場用磁石単体による静磁場誤差および被検体が静磁場内に入ることによる静磁場の誤差の総和で定義することができる。 Factors that disturb the uniformity of the static magnetic field intensity and cause distortion in the reconstructed image include errors in the gradient magnetic field, static magnetic field errors due to the single static magnetic field magnet, and static magnetic field due to the subject entering the static magnetic field. There are three factors of error. Therefore, the static magnetic field error can be defined as the sum of the gradient magnetic field error, the static magnetic field error due to the single static magnetic field magnet, and the static magnetic field error due to the subject entering the static magnetic field.
さらに、静磁場の誤差から画像データの歪み量が画像空間の位置ごとに求められ、求められた画像データの歪み量に基づいて画像空間の位置ごとに画像値が補正される。すなわち、まず画像空間において、画像値が得られている３次元状の位置が、歪みがなかった場合の位置にシフトされる。しかし、シフトされた位置は、必ずしも本来あるべき格子点上とならないため、シフト後の各位置における複数の画像値から補間によって格子点上の画像値が求められる。
従来行われる静磁場強度の不均一性による画像歪みの補正においては、静磁場の各誤差要因に応じて静磁場の歪み量の分布を推定することが精度上重要である。特に静磁場内に入ることによる静磁場の誤差については、被検体が動くことによって変動するので、静磁場の誤差分布の推定値が誤差を持つ可能性がある。 In correction of image distortion due to static magnetic field intensity non-uniformity performed in the past, it is important in terms of accuracy to estimate the distribution of the amount of distortion of the static magnetic field according to each error factor of the static magnetic field. In particular, the error of the static magnetic field due to entering the static magnetic field varies depending on the movement of the subject, so there is a possibility that the estimated value of the error distribution of the static magnetic field has an error.
そして、静磁場の誤差分布が精度よく推定されないまま、静磁場の誤差分布に基づいて画像の補正が行われると、適切に画像歪みが改善されない恐れがある。このため、静磁場の誤差分布が精度よく推定されているか否か、あるいは静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価する技術の開発が望まれる。 If the image correction is performed based on the static magnetic field error distribution without accurately estimating the static magnetic field error distribution, the image distortion may not be appropriately improved. Therefore, it is desired to develop a technique for evaluating whether or not the error distribution of the static magnetic field is accurately estimated, or whether or not the image distortion correction due to the non-uniformity of the static magnetic field has been appropriately performed.
本発明はかかる従来の事情に対処するためになされたものであり、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することが可能な磁気共鳴イメージング装置および画像補正評価方法を提供することを目的とする。 The present invention has been made in order to cope with such a conventional situation. The error distribution of the static magnetic field formed in the imaging region is appropriately estimated, and the image distortion due to the static magnetic field inhomogeneity is appropriately corrected. It is an object of the present invention to provide a magnetic resonance imaging apparatus and an image correction evaluation method capable of evaluating whether or not a failure has occurred.
本発明に係る磁気共鳴イメージング装置は、上述の目的を達成するために、請求項１に記載したように、位相エンコード方向、リードアウト方向およびスライスエンコード方向の少なくとも１方向の傾斜磁場の極性を反転させてスキャンが実行されるように複数通りの撮影条件を設定する撮影条件設定手段と、前記撮影条件設定手段により設定された前記各撮影条件に従って撮影領域に高周波信号を送信することによる原子核の核磁気共鳴に伴って発生した磁気共鳴信号に基づいてデータを生成するデータ収集手段と、前記データに画像再構成処理を施すことにより前記撮影条件ごとに画像データをそれぞれ再構成する画像再構成手段と、前記撮影領域の静磁場分布に基づいて前記各撮影条件により得られた前記各画像データの歪みをそれぞれ補正する画像歪み補正手段と、前記画像歪み補正手段による補正後における前記各画像データに基づいて前記各画像データの補正が適切に行われたか否かを評価する画像補正評価手段とを有することを特徴とするものである。
In order to achieve the above object, a magnetic resonance imaging apparatus according to the present invention reverses the polarity of a gradient magnetic field in at least one of a phase encoding direction, a readout direction, and a slice encoding direction as described in
また、本発明に係る画像補正評価方法は、上述の目的を達成するために、請求項４に記載したように、位相エンコード方向、リードアウト方向およびスライスエンコード方向の少なくとも１方向の傾斜磁場の極性を反転させた複数通りの撮影条件に従うスキャンの実行により収集されたデータに画像再構成処理を施すことにより前記撮影条件ごとに画像データをそれぞれ再構成するステップと、前記画像データの撮影領域における静磁場分布に基づいて前記各撮影条件により得られた前記各画像データの歪みをそれぞれ補正するステップと、補正後における前記各画像データに基づいて前記各画像データの補正が適切に行われたか否かを評価するステップとを有することを特徴とするものである。
In order to achieve the above object, the image correction evaluation method according to the present invention provides the polarity of the gradient magnetic field in at least one of the phase encoding direction, the readout direction, and the slice encoding direction as described in
本発明に係る磁気共鳴イメージング装置および画像補正評価方法においては、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することができる。 In the magnetic resonance imaging apparatus and the image correction evaluation method according to the present invention, whether the error distribution of the static magnetic field formed in the imaging region is appropriately estimated and the image distortion due to the static magnetic field inhomogeneity is appropriately corrected. You can evaluate whether or not.
本発明に係る磁気共鳴イメージング装置および画像補正評価方法の実施の形態について添付図面を参照して説明する。 Embodiments of a magnetic resonance imaging apparatus and an image correction evaluation method according to the present invention will be described with reference to the accompanying drawings.
図１は本発明に係る磁気共鳴イメージング装置の実施の形態を示す構成図である。 FIG. 1 is a configuration diagram showing an embodiment of a magnetic resonance imaging apparatus according to the present invention.
磁気共鳴イメージング装置２０は、静磁場を形成する筒状の静磁場用磁石２１、この静磁場用磁石２１の内部に設けられたシムコイル２２、傾斜磁場コイルユニット２３およびＲＦコイル２４を図示しないガントリに内蔵した構成である。
The magnetic
また、磁気共鳴イメージング装置２０には、制御系２５が備えられる。制御系２５は、静磁場電源２６、傾斜磁場電源２７、シムコイル電源２８、送信器２９、受信器３０、シーケンスコントローラ３１およびコンピュータ３２を具備している。制御系２５の傾斜磁場電源２７は、Ｘ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚで構成される。また、コンピュータ３２には、入力装置３３、表示装置３４、演算装置３５および記憶装置３６が備えられる。
In addition, the magnetic
静磁場用磁石２１は静磁場電源２６と接続され、静磁場電源２６から供給された電流により撮像領域に静磁場を形成させる機能を有する。尚、静磁場用磁石２１は超伝導コイルで構成される場合が多く、励磁の際に静磁場電源２６と接続されて電流が供給されるが、一旦励磁された後は非接続状態とされるのが一般的である。また、静磁場用磁石２１を永久磁石で構成し、静磁場電源２６が設けられない場合もある。
The static
また、静磁場用磁石２１の内側には、同軸上に筒状のシムコイル２２が設けられる。シムコイル２２はシムコイル電源２８と接続され、シムコイル電源２８からシムコイル２２に電流が供給されて静磁場が均一化されるように構成される。
A cylindrical shim coil 22 is coaxially provided inside the static
傾斜磁場コイルユニット２３は、Ｘ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚで構成され、静磁場用磁石２１の内部において筒状に形成される。傾斜磁場コイルユニット２３の内側には寝台３７が設けられて撮像領域とされ、寝台３７には被検体Ｐがセットされる。ＲＦコイル２４はガントリに内蔵されず、寝台３７や被検体Ｐ近傍に設けられる場合もある。
The gradient magnetic
また、傾斜磁場コイルユニット２３は、傾斜磁場電源２７と接続される。傾斜磁場コイルユニット２３のＸ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚはそれぞれ、傾斜磁場電源２７のＸ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚと接続される。
The gradient magnetic
そして、Ｘ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚからそれぞれＸ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚに供給された電流により、撮像領域にそれぞれＸ軸方向の傾斜磁場Ｇｘ、Ｙ軸方向の傾斜磁場Ｇｙ、Ｚ軸方向の傾斜磁場Ｇｚを形成することができるように構成される。
The X-axis gradient magnetic
ＲＦコイル２４は、送信器２９および受信器３０と接続される。ＲＦコイル２４は、送信器２９からＲＦ信号を受けて被検体Ｐに送信する機能と、被検体Ｐ内部の原子核スピンのＲＦ信号による励起に伴って発生したＭＲ信号を受信して受信器３０に与える機能を有する。
The RF coil 24 is connected to the
一方、制御系２５のシーケンスコントローラ３１は、傾斜磁場電源２７、送信器２９および受信器３０と接続される。シーケンスコントローラ３１は傾斜磁場電源２７、送信器２９および受信器３０を駆動させるために必要な制御情報、例えば傾斜磁場電源２７に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したパルスシーケンスを記憶する機能と、記憶した所定のパルスシーケンスで規定される撮影条件に従って傾斜磁場電源２７、送信器２９および受信器３０を駆動させることによりＸ軸傾斜磁場Ｇｘ、Ｙ軸傾斜磁場Ｇｙ，Ｚ軸傾斜磁場ＧｚおよびＲＦ信号を発生させる機能を有する。
On the other hand, the
また、シーケンスコントローラ３１は、受信器３０におけるＭＲ信号の検波およびＡ／Ｄ変換により得られた複素データである生データ（raw data）を受けてコンピュータ３２に与えるように構成される。
The
このため、送信器２９には、シーケンスコントローラ３１から受けた制御情報に基づいてＲＦ信号をＲＦコイル２４に与える機能が備えられる一方、受信器３０には、ＲＦコイル２４から受けたＭＲ信号を検波して所要の信号処理を実行するとともにＡ／Ｄ変換することにより、デジタル化された複素データである生データを生成する機能と生成した生データをシーケンスコントローラ３１に与える機能とが備えられる。
For this reason, the
すなわち、静磁場用磁石２１、シムコイル２２、傾斜磁場コイルユニット２３、ＲＦコイル２４並びに制御系２５の各構成要素により、磁気共鳴イメージング装置２０は、パルスシーケンスとして設定された各撮影条件に従って、静磁場中の被検体Ｐに対して傾斜磁場の印加およびＲＦ信号の送信を行なう一方、被検体Ｐ内部におけるＲＦ信号による原子核の核磁気共鳴に伴って発生したＭＲ信号を受信し、デジタル化することにより生データを生成する。
That is, the magnetic
また、コンピュータ３２の記憶装置３６に保存されたプログラムを演算装置３５で実行することにより、コンピュータ３２には各種機能が備えられる。ただし、プログラムによらず、特定の回路を設けてコンピュータ３２を構成してもよい。
Further, the
図２は、図１に示す磁気共鳴イメージング装置２０におけるコンピュータ３２の機能ブロック図である。
FIG. 2 is a functional block diagram of the
コンピュータ３２は、画像補正評価プログラムによりシーケンスコントローラ制御装置４０、生データデータベース４１、画像再構成装置４２、画像データデータベース４３、静磁場分布データベース４４、画像歪み補正装置４５、撮影条件設定装置４６、画像補正評価装置の一例としての相互相関計算装置４７、画像補正評価装置の一例としての信号強度比較装置４８および磁場分布計算装置４９として機能する。
The
シーケンスコントローラ制御装置４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に所要のパルスシーケンスを与えることにより駆動制御させる機能と、シーケンスコントローラ３１から生データを受けて生データデータベース４１に形成されたｋ空間（フーリエ空間）に配置する機能とを有する。
The sequence controller control device 40 receives the raw data from the
撮影条件設定装置４６は、パルスシーケンスを生成することにより撮影条件を設定する機能と、生成したパルスシーケンスをシーケンスコントローラ制御装置４０に与える機能とを有する。ここで、撮影条件設定装置４６には、静磁場分布測定用のプリスキャンの実行のためのパルスシーケンスを生成する他、位相エンコード（ＰＥ）方向、リードアウト（ＲＯ）方向およびスライスエンコード（ＳＥ）方向の少なくとも１方向の傾斜磁場の極性を反転させて、複数通りの撮影条件によりスキャンが実行されるように複数通りのパルスシーケンスを生成する機能が備えられる。
The imaging
このため、生データデータベース４１には、撮影条件設定装置４６により設定された複数通りの撮影条件のスキャンの実行によって受信器３０において生成された各生データが保存され、生データデータベース４１に形成されたｋ空間に各生データがそれぞれ撮影条件ごとに配置される。また、静磁場分布測定用のプリスキャンが実行された場合には、静磁場分布測定用の生データが保存される。
Therefore, the
画像再構成装置４２は、生データデータベース４１から生データを取り込んでフーリエ変換処理や最大値投影処理等の所定の画像再構成処理を施すことにより、被検体Ｐの画像データを再構成する機能と、得られた画像データを画像データデータベース４３に書き込む機能とを有する。ただし、診断時以外には、被検体Ｐ以外の例えばファントムの画像データを再構成するようにしてもよい。
The
このため、画像データデータベース４３には、撮影条件ごとの被検体Ｐやファントムの画像データが保存される。
Therefore, the
静磁場分布データベース４４には、予め静磁場分布が計測あるいは求められて保存される。静磁場分布データベース４４には、磁場分布計算装置４９によって求められた静磁場分布を保存することもできる。
In the static magnetic
画像歪み補正装置４５は、画像データデータベース４３に保存された撮影条件ごとの被検体Ｐやファントムの画像データの歪みを静磁場分布データベース４４から読み込んだ静磁場分布に基づいてそれぞれ補正する機能と、補正後の各画像データを相互相関計算装置４７および信号強度比較装置４８にそれぞれ与える機能とを有する。
The image
相互相関計算装置４７は、複数通りの撮影条件によって得られた補正後の画像データを画像歪み補正装置４５から受けて、複数の画像データ間の相互相関係数を求める機能と、求められた相互相関係数に基づいて画像データの補正が適切に行われたか否かを評価する画像補正評価装置としての機能とを有する。また、必要に応じて相互相関計算装置４７には、求められた複数の画像データ間の相互相関係数や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えることにより表示させる機能が備えられる。
The
信号強度比較装置４８は、複数通りの撮影条件によって得られた補正後の画像データを画像歪み補正装置４５から受けて、複数の画像データ間でピクセルごとに信号強度を比較する機能と、信号強度の比較結果にに基づいて画像データの補正が適切に行われたか否かを評価する画像補正評価装置としての機能とを有する。また、必要に応じて信号強度比較装置４８には、複数の画像データ間における信号強度の比較結果や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えることにより表示させる機能が備えられる。
The signal
磁場分布計算装置４９は、任意の方法により静磁場の分布を推定する機能と、推定して得られた静磁場分布データを静磁場分布データベース４４に書き込む機能とを有する。静磁場の分布の推定には、必要に応じて生データデータベース４１に保存された静磁場分布測定用の生データが参照される。さらに、傾斜磁場を測定するための傾斜磁場測定コイル５０を撮影領域近傍に設けて、傾斜磁場測定コイル５０による傾斜磁場の測定結果を読み込んで参照するように構成することもできる。
The magnetic field
次に磁気共鳴イメージング装置２０の動作について説明する。
Next, the operation of the magnetic
図３は、図１に示す磁気共鳴イメージング装置２０により、静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価する際の流れの一例を示すフローチャートであり、図中Ｓに数字を付した符号はフローチャートの各ステップを示す。
FIG. 3 is a flowchart showing an example of a flow when evaluating whether or not the image distortion correction due to the non-uniformity of the static magnetic field has been appropriately performed by the magnetic
まずステップＳ１において、予め任意の方法により静磁場の誤差分布が推定されて静磁場分布データベース４４に保存される。静磁場の誤差分布は計測してもよい。
First, in
静磁場の誤差分布の推定手法の一例としては、静磁場の誤差の発生要因を傾斜磁場の誤差、静磁場用磁石２１単体による静磁場誤差および被検体Ｐが静磁場内に入ることによる静磁場の誤差の３つの要因に分類し、各要因による静磁場の誤差をそれぞれ推定した後、加算する方法がある。
As an example of the estimation method of the static magnetic field error distribution, the cause of the static magnetic field error is the gradient magnetic field error, the static magnetic field error due to the static
第１の要因である傾斜磁場の誤差は、理想的には線形であるべき傾斜磁場の線形からの誤差により生じるものである。傾斜磁場の誤差に基づく静磁場の歪み量は、傾斜磁場の強度とは独立して決定される。 The error of the gradient magnetic field, which is the first factor, is caused by an error from the linearity of the gradient magnetic field that should ideally be linear. The amount of distortion of the static magnetic field based on the gradient magnetic field error is determined independently of the strength of the gradient magnetic field.
撮影領域におけるガントリ座標系の任意の点の座標をｘ、ｙ、ｚとすると、傾斜磁場の誤差は、ｘ，ｙ，ｚ方向の各チャンネルにおいて勾配を有する３次元状の分布となる。そして、理想的な線形の傾斜磁場の強度をＧｘ，Ｇｙ，Ｇｚとして、ｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差の割合ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）、ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）、ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）を、式（１−１）、式（１−２）、式（１−３）のように表すと、ΔＧｘ（ｘ，ｙ，ｚ）、ΔＧｙ（ｘ，ｙ，ｚ）、ΔＧｚ（ｘ，ｙ，ｚ）は、ｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差量を意味することになる。 If the coordinates of an arbitrary point of the gantry coordinate system in the imaging region are x, y, and z, the error of the gradient magnetic field becomes a three-dimensional distribution having a gradient in each channel in the x, y, and z directions. Then, assuming that the intensity of the ideal linear gradient magnetic field is Gx, Gy, Gz, the ratio of the gradient magnetic field error in each channel in the x, y, z directions GxErrorRatio (x, y, z), GyErrorRatio (x, y, z), GzErrorRatio (x, y, z) can be expressed as Equation 1-1, Equation 1-2, and Equation 1-3, and ΔGx (x, y, z), ΔGy ( x, y, z) and ΔGz (x, y, z) mean an error amount of the gradient magnetic field in each channel in the x, y, and z directions.
［数１］
ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｘ（ｘ，ｙ，ｚ）／Ｇｘ…（１−１）
ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｙ（ｘ，ｙ，ｚ）／Ｇｙ…（１−２）
ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｚ（ｘ，ｙ，ｚ）／Ｇｚ…（１−３）
[Equation 1]
GxErrorRatio (x, y, z) = ΔGx (x, y, z) / Gx (1-1)
GyErrorRatio (x, y, z) = ΔGy (x, y, z) / Gy (1-2)
GzErrorRatio (x, y, z) = ΔGz (x, y, z) / Gz (1-3)
さらに、傾斜磁場の誤差量ΔＧｘ（ｘ，ｙ，ｚ）、ΔＧｙ（ｘ，ｙ，ｚ）、ΔＧｚ（ｘ，ｙ，ｚ）と、ガントリ座標系の座標ｘ、ｙ、ｚとからｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）を式（２−１）、式（２−２）、式（２−３）のように推定することができる。 Further, the gradient magnetic field errors ΔGx (x, y, z), ΔGy (x, y, z), ΔGz (x, y, z) and the coordinates x, y, z of the gantry coordinate system are x, y. , Z-direction static magnetic field error distribution ΔBgradx (x, y, z), ΔBgrady (x, y, z), ΔBgradz (x, y, z) due to errors in the gradient magnetic field in each channel in the z and z directions are expressed by the formula (2- 1), equation (2-2), and equation (2-3).
［数２］
ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）＝ΔＧｘ（ｘ，ｙ，ｚ）×ｘ …（２−１）
ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）＝ΔＧｙ（ｘ，ｙ，ｚ）×ｙ …（２−２）
ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）＝ΔＧｚ（ｘ，ｙ，ｚ）×ｚ …（２−３）
[Equation 2]
ΔBgradx (x, y, z) = ΔGx (x, y, z) × x (2-1)
ΔBgrady (x, y, z) = ΔGy (x, y, z) × y (2-2)
ΔBgradz (x, y, z) = ΔGz (x, y, z) × z (2-3)
また、第２の要因である静磁場用磁石２１単体による静磁場誤差は、理想的な静磁場の強度からの実際の静磁場の強度の誤差により生じるものである。静磁場用磁石２１単体による静磁場誤差は、撮像視野（ＦＯＶ：ｆｉｅｌｄ ｏｆ ｖｉｅｗ）が大きくなると急激に悪化するという性質がある。また、静磁場用磁石２１単体による静磁場誤差は、傾斜磁場の強度に依存する一方、静磁場中における被検体Ｐの形状や有無からは独立している。
Further, the static magnetic field error due to the static
また、第３の要因である被検体Ｐが静磁場内に入ることによる静磁場の誤差は、被検体Ｐが静磁場内に入ることで磁化率に違いが生じた結果、生じる静磁場の誤差分である。被検体Ｐが静磁場内に入ることによる静磁場の誤差は、傾斜磁場の強度に依存し、かつ静磁場中における被検体Ｐの形状や有無にも依存する。 In addition, the static magnetic field error caused by the subject P entering the static magnetic field, which is the third factor, is a static magnetic field error resulting from the difference in magnetic susceptibility caused by the subject P entering the static magnetic field. Minutes. The error of the static magnetic field due to the subject P entering the static magnetic field depends on the strength of the gradient magnetic field and also depends on the shape and presence of the subject P in the static magnetic field.
そして、静磁場の誤差分布は、第１、第２、第３の各要因による静磁場の誤差成分の和として定義することができる。すなわち、ガントリ座標系の座標ｘ、ｙ、ｚにおける静磁場の誤差をΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）、傾斜磁場に起因する誤差成分をΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）、静磁場用磁石２１単体による静磁場誤差をΔＢｍａｇ（ｘ，ｙ，ｚ）、被検体Ｐが静磁場内に入ることによる静磁場の誤差をΔＢｏｂｊ（ｘ，ｙ，ｚ）とすると、式（３−１）、式（３−２）、式（３−３）のように表すことができる。
The static magnetic field error distribution can be defined as the sum of static magnetic field error components due to the first, second, and third factors. That is, the error of the static magnetic field at the coordinates x, y, z in the gantry coordinate system is caused by ΔBox (x, y, z), Δ Boy (x, y, z), ΔBoz (x, y, z), and the gradient magnetic field. The error components are ΔBgradx (x, y, z), ΔBgrady (x, y, z), ΔBgradz (x, y, z), the static magnetic field error due to the single static
［数３］
ΔＢｏｘ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−１）
ΔＢｏｙ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−２）
ΔＢｏｚ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−３）
[Equation 3]
ΔBox (x, y, z) = ΔBgradx (x, y, z) + ΔBmag (x, y, z) + ΔBobj (x, y, z) (3-1)
ΔBoy (x, y, z) = ΔBgrady (x, y, z) + ΔBmag (x, y, z) + ΔBobj (x, y, z) (3-2)
ΔBoz (x, y, z) = ΔBgradz (x, y, z) + ΔBmag (x, y, z) + ΔBobj (x, y, z) (3-3)
従って、式（３−１）、式（３−２）、式（３−３）からガントリ座標系のｘ、ｙ、ｚ方向の静磁場の歪み量ΔＸ、ΔＹ，ΔＺは、式（４−１）、式（４−２）、式（４−３）のように求めることができる。 Therefore, the distortion amounts ΔX, ΔY, and ΔZ of the static magnetic field in the x, y, and z directions of the gantry coordinate system are expressed by the equation (4-) from the equations (3-1), (3-2), and (3-3). 1), equation (4-2), and equation (4-3).
［数４］
ΔＸ＝ΔＢｏｘ（ｘ，ｙ，ｚ）／Ｇｘ
＝ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｘ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｘ …（４−１）
ΔＹ＝ΔＢｏｙ（ｘ，ｙ，ｚ）／Ｇｙ
＝ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｙ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｙ …（４−２）
ΔＺ＝ΔＢｏｚ（ｘ，ｙ，ｚ）／Ｇｚ
＝ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｚ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｚ …（４−３）
[Equation 4]
ΔX = ΔBox (x, y, z) / Gx
= GxErrorRatio (x, y, z) × x + {ΔBmag (x, y, z) + ΔBobj (x, y, z)} / Gx (4-1)
ΔY = ΔBoy (x, y, z) / Gy
= GyErrorRatio (x, y, z) × y + {ΔBmag (x, y, z) + ΔBobj (x, y, z)} / Gy (4-2)
ΔZ = ΔBoz (x, y, z) / Gz
= GzErrorRatio (x, y, z) × z + {ΔBmag (x, y, z) + ΔBobj (x, y, z)} / Gz (4-3)
式（４−１）、式（４−２）、式（４−３）において、｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝の値は、イメージング用のスキャンに先立って静磁場分布測定用のプリスキャンを実行することにより求めることができる。イメージング用のスキャンがＥＰＩシーケンスによるスキャンである場合には、静磁場分布測定用のプリスキャンをＥＰＩシーケンスの受信ゲイン調整や傾斜磁場の遅れ時間調整のためのプリスキャンとともに実行することができる。 In equations (4-1), (4-2), and (4-3), the value of {ΔBmag (x, y, z) + ΔBobj (x, y, z)} is prior to the scan for imaging. This can be obtained by executing a pre-scan for measuring the static magnetic field distribution. When the imaging scan is a scan using an EPI sequence, the pre-scan for measuring the static magnetic field distribution can be executed together with the pre-scan for adjusting the reception gain of the EPI sequence and adjusting the delay time of the gradient magnetic field.
すなわち、静磁場分布測定用のプリスキャンの実行により設定されたスライス面内の静磁場分布（中心周波数からのずれを位相で表した分布図）を測定するためのデータを収集することができる。この場合、静磁場分布測定用のプリスキャンの実行の際における面内分解能はイメージング用のスキャンの実行の際における面内分解能よりも粗くすることができる。 That is, it is possible to collect data for measuring the static magnetic field distribution (distribution map in which the deviation from the center frequency is represented by the phase) in the slice plane set by executing the pre-scan for measuring the static magnetic field distribution. In this case, the in-plane resolution at the time of executing the pre-scan for the static magnetic field distribution measurement can be made coarser than the in-plane resolution at the time of executing the imaging scan.
静磁場分布測定用のプリスキャンの一例としては、ＦｉｅｌｄＥｃｈｏ法によりエコー時間（ＴＥ：ｅｃｈｏ ｔｉｍｅ）を変えて２つの画像データを収集するものが挙げられる。そして、式（５）に示すように、静磁場分布測定用のプリスキャンにより収集された各画像データのピクセルごとの位相差Δθ（ｘ，ｙ，ｚ）から｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝の値を求めることができる。 As an example of the pre-scan for measuring the static magnetic field distribution, one that collects two image data by changing the echo time (TE: echo time) by the FieldEcho method can be cited. Then, as shown in equation (5), {ΔBmag (x, y, z) from the phase difference Δθ (x, y, z) for each pixel of each image data collected by the pre-scan for measuring the static magnetic field distribution. The value of + ΔBobj (x, y, z)} can be obtained.
［数５］
Δθ（ｘ，ｙ，ｚ）＝（ＴＥ１−ＴＥ２）×｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝ …（５）
但し、式（５）において、ＴＥ１、ＴＥ２は、それぞれ互いに異なるＴＥである。
[Equation 5]
Δθ (x, y, z) = (TE1−TE2) × {ΔBmag (x, y, z) + ΔBobj (x, y, z)} (5)
However, in Formula (5), TE1 and TE2 are TEs different from each other.
つまり、ＴＥを変えて収集された画像データを差分することで、静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）のうちΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）の項がキャンセルされる。 In other words, ΔBgradx out of the static magnetic field errors ΔBox (x, y, z), ΔBoy (x, y, z), ΔBoz (x, y, z) by differentiating the image data collected by changing TE. The terms (x, y, z), ΔBgrady (x, y, z), and ΔBgradz (x, y, z) are canceled.
従って、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）を求めることができれば、式（３−１），式（３−２）、式（３−３）および式（４−１），式（４−２）、式（４−３）から静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）および静磁場の歪み量ΔＸ、ΔＹ，ΔＺを求めることができる。 Accordingly, the static magnetic field error distribution ΔBgrad. x (x, y, z), ΔBgrad. y (x, y, z), ΔBgrad. If z (x, y, z) can be obtained, formula (3-1), formula (3-2), formula (3-3), formula (4-1), formula (4-2), formula From (4-3), the static magnetic field errors ΔBox (x, y, z), ΔBoy (x, y, z), ΔBoz (x, y, z) and the static magnetic field distortion amounts ΔX, ΔY, ΔZ are obtained. Can do.
傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）は、例えば傾斜磁場測定コイル５０により測定することができる。
Error distribution ΔBgrad. Of the static magnetic field due to the error of the gradient magnetic field. x (x, y, z), ΔBgrad. y (x, y, z), ΔBgrad. z (x, y, z) can be measured by the gradient magnetic
そこで、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）が傾斜磁場測定コイル５０により測定されて磁場分布計算装置４９に与えられる。
Therefore, the static magnetic field error distribution ΔBgrad. x (x, y, z), ΔBgrad. y (x, y, z), ΔBgrad. z (x, y, z) is measured by the gradient magnetic
一方、撮影条件設定装置４６により静磁場分布測定用のプリスキャンのためのパルスシーケンスが生成されてシーケンスコントローラ制御装置４０に与えられる。このため、シーケンスコントローラ制御装置４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に静磁場分布測定用のパルスシーケンスを与えることにより静磁場分布測定用のプリスキャンを実行させる。
On the other hand, a pulse sequence for pre-scan for static magnetic field distribution measurement is generated by the imaging
すなわち、予め寝台３７に被検体Ｐあるいはファントムがセットされ、静磁場電源２６により励磁された静磁場用磁石２１（超伝導磁石）の撮像領域に静磁場が形成される。また、シムコイル電源２８からシムコイル２２に電流が供給されて撮像領域に形成された静磁場が均一化される。
That is, the subject P or phantom is set on the
さらに、シーケンスコントローラ３１は、シーケンスコントローラ制御装置４０から受けた静磁場分布測定用のパルスシーケンスに従って傾斜磁場電源２７、送信器２９および受信器３０を駆動させることにより被検体Ｐあるいはファントムがセットされた撮像領域にＸ軸傾斜磁場Ｇｘ、Ｙ軸傾斜磁場Ｇｙ，Ｚ軸傾斜磁場Ｇｚを形成させるとともに、ＲＦ信号を発生させる。
Furthermore, the
このため、送信器２９からパルスシーケンスに応じてＲＦコイル２４にＲＦ信号が与えられ、ＲＦコイル２４から被検体ＰにＲＦ信号が送信される。ＲＦ信号が被検体Ｐに送信されると、被検体Ｐの内部において原子核の核磁気共鳴により生じたＭＲ信号が、ＲＦコイル２４により受信されて受信器３０に与えられる。受信器３０は、ＲＦコイル２４からＭＲ信号を受けて、前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリング等の各種信号処理を実行する。さらに、受信器３０は、ＭＲ信号をＡ／Ｄ変換することにより、デジタルデータのＭＲ信号である生データを生成する。そして、受信器３０は、生成した生データをシーケンスコントローラ３１に与える。
Therefore, an RF signal is given from the
シーケンスコントローラ３１は、受信器３０から受けた生データをシーケンスコントローラ制御装置４０に与え、シーケンスコントローラ制御装置４０は生データデータベース４２に生データを入力する。このようにして、例えばＴＥを変えたＦｉｅｌｄＥｃｈｏ法による２回に亘る静磁場分布測定用のプリスキャンが実行され、生データデータベース４２には、静磁場分布測定用のプリスキャンにより収集された生データが蓄積される。
The
そして、磁場分布計算装置４９は、静磁場分布測定用のプリスキャンの実行により収集された生データを生データデータベース４１から読み込む一方、傾斜磁場測定コイル５０から傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）を受けて、例えば式（３−１），式（３−２）、式（３−３）、式（４−１），式（４−２）、式（４−３）および式（５）により静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）および静磁場の歪み量ΔＸ、ΔＹ，ΔＺを求める。
The magnetic field
さらに、磁場分布計算装置４９は、求められた静磁場の歪み量ΔＸ、ΔＹ，ΔＺを静磁場分布データとして静磁場分布データベース４４に書き込む。
Further, the magnetic field
尚、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）は、傾斜磁場測定コイル５０による測定の他、例えば系を模擬した磁場モデルを用いてシミュレーションによる計算で求めることもできる。
Note that the static magnetic field error distribution ΔBgrad. x (x, y, z), ΔBgrad. y (x, y, z), ΔBgrad. z (x, y, z) can be obtained by calculation by simulation using a magnetic field model simulating the system, for example, in addition to measurement by the gradient magnetic
次に、ステップＳ２において、撮影条件設定装置４６により被検体Ｐの画像収集用のパルスシーケンスが生成されることにより、撮影条件が設定される。このとき、生成される画像収集用のパルスシーケンスは、ＰＥ方向、ＲＯ方向およびＳＥ方向の少なくとも１方向の傾斜磁場の極性を反転させて、例えば２通りの撮影条件によりスキャンが実行されるような２通りのパルスシーケンスとされる。
Next, in step S <b> 2, the imaging condition is set by generating a pulse sequence for collecting images of the subject P by the imaging
図４は、図１に示す磁気共鳴イメージング装置２０により静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価するために生成されるパルスシーケンスの一例を示す図である。
FIG. 4 is a diagram illustrating an example of a pulse sequence generated for evaluating whether or not the image distortion correction due to the static magnetic field inhomogeneity has been appropriately performed by the magnetic
図４において、ＳＥ１、ＲＯ１、ＰＥ１は、それぞれＳＥ方向、ＲＯ方向およびＰＥ方向の傾斜磁場を規定する第１のパルスシーケンスであり、ＳＥ２、ＲＯ２、ＰＥ２は、それぞれＳＥ方向、ＲＯ方向およびＰＥ方向の傾斜磁場を規定する第２のパルスシーケンスである。図４に示すように、例えばＰＥ方向およびＲＯ方向の２方向の傾斜磁場の極性を互いに反転させて、２通りの撮影条件によりスキャンが実行されるような第１、第２の２通りのパルスシーケンスが撮影条件設定装置４６により生成される。
In FIG. 4, SE1, RO1, and PE1 are first pulse sequences that define gradient magnetic fields in the SE direction, RO direction, and PE direction, respectively. SE2, RO2, and PE2 are SE direction, RO direction, and PE direction, respectively. It is the 2nd pulse sequence which prescribes | regulates the gradient magnetic field. As shown in FIG. 4, for example, first and second two kinds of pulses in which the polarities of the gradient magnetic fields in two directions of the PE direction and the RO direction are reversed with each other and the scan is executed under two kinds of imaging conditions. A sequence is generated by the imaging
そして、撮影条件設定装置４６により生成された２通りのパルスシーケンスは、シーケンスコントローラ制御装置４０に与える。
Then, the two pulse sequences generated by the imaging
次に、ステップＳ３において、撮影条件設定装置４６により設定された２つの撮影条件に従ってスキャンが実行され、画像データが再構成される。すなわち、シーケンスコントローラ制御装置４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に画像収集用の２つのパルスシーケンスを与えることにより画像収集用のスキャンを実行させる。
Next, in step S3, scanning is executed in accordance with the two shooting conditions set by the shooting
この結果、傾斜磁場の極性を互いに反転させて得られた各生データが、生データデータベース４１に形成されたｋ空間に撮影条件ごとにそれぞれ配置される。
As a result, each raw data obtained by reversing the polarities of the gradient magnetic fields is arranged in the k space formed in the
さらに、画像再構成装置４２は、生データデータベース４１から各生データを取り込んでフーリエ変換処理等の所定の画像再構成処理をそれぞれ施すことにより、２つの画像データを再構成する。このようにして得られた各画像データは、傾斜磁場を反転させて収集された生データから再構成されたため、互いに歪み方が異なる画像データとなる。そして、各画像データは画像データデータベース４３に書き込まれて保存される。
Further, the
次に、ステップＳ４において、画像データデータベース４３に保存された歪み方の異なる各画像データが画像歪み補正装置４５に読み込まれて、それぞれ静磁場分布データベース４４に保存された静磁場分布に基づく各画像データの歪みの補正が行われる。
Next, in step S <b> 4, each image data stored in the
図５は、図１に示す磁気共鳴イメージング装置２０により静磁場分布に基づいて画像データを補正する際の方法の一例を示す説明図である。
FIG. 5 is an explanatory diagram showing an example of a method for correcting image data based on the static magnetic field distribution by the magnetic
まず、静磁場分布から各画像データの歪み量が画像空間の位置ごとに求められる。例えばｘ方向のみに無視できない歪が存在する場合には、ｘ方向の歪み量Δｘが求められる。そして、歪み量Δｘに基づいて、例えば歪み後の２点の座標（Ｘ１，Ｙ１，Ｚ１）、（Ｘ２，Ｙ２，Ｚ２）におけるそれぞれの画像値Ｐ１（Ｘ１，Ｙ１，Ｚ１）、Ｐ２（Ｘ２，Ｙ２，Ｚ２）が歪み前における各座標（Ｘ１’，Ｙ１’，Ｚ１’）、（Ｘ２’，Ｙ２’，Ｚ２’）における画像値Ｐ１（Ｘ１’，Ｙ１’，Ｚ１’）、Ｐ２（Ｘ２’，Ｙ２’，Ｚ２’）とされる。つまり、ｘ方向に歪み量Δｘだけ画像値Ｐ１（Ｘ１，Ｙ１，Ｚ１）、Ｐ２（Ｘ２，Ｙ２，Ｚ２）の位置（Ｘ１，Ｙ１，Ｚ１）、（Ｘ２，Ｙ２，Ｚ２）がシフトされる。
First, the amount of distortion of each image data is obtained for each position in the image space from the static magnetic field distribution. For example, when there is a distortion that cannot be ignored only in the x direction, the distortion amount Δx in the x direction is obtained. Based on the distortion amount Δx, for example, the respective image values P 1 (X 1 , Y 2 ) at the coordinates (X 1 , Y 1 , Z 1 ) and (X 2 , Y 2 , Z 2 ) of the two points after distortion. 1 , Z 1 ), P 2 (X 2 , Y 2 , Z 2 ) are coordinates (X 1 ′, Y 1 ′, Z 1 ′), (X 2 ′, Y 2 ′, Z 2 ′) before distortion. ) Image values P 1 (X 1 ′, Y 1 ′, Z 1 ′) and P 2 (X 2 ′, Y 2 ′, Z 2 ′). That is, the image value by distortion amount Δx in the x-direction P 1 (X 1, Y 1 , Z 1), the position of P 2 (X 2, Y 2 , Z 2) (
さらに、本来あるべき格子点の座標（ｘ，ｙ，ｚ）における画像値Ｐ（ｘ，ｙ，ｚ）が、格子点の座標（ｘ，ｙ，ｚ）近傍の位置（Ｘ１’，Ｙ１’，Ｚ１’）、（Ｘ２’，Ｙ２’，Ｚ２’）における画像値Ｐ１（Ｘ１’，Ｙ１’，Ｚ１’）、Ｐ２（Ｘ２’，Ｙ２’，Ｚ２’）から補間によって画像空間の位置ごとに求められる。このとき通常、補間にはｃｕｂｉｃ ｓｐｌｉｎｅ以上の関数が用いられ、Ｌｉｎｅａｒではボケが大きく不十分であると考えられる。 Further, the image value P (x, y, z) at the coordinates (x, y, z) of the lattice point that should be originally exists at the position (X 1 ′, Y 1 ) in the vicinity of the coordinates (x, y, z) of the lattice point. ', Z 1 '), (X 2 ', Y 2 ', Z 2 ') image values P 1 (X 1 ', Y 1 ', Z 1 '), P 2 (X 2 ', Y 2 ', Z 2 ′) is obtained for each position in the image space by interpolation. At this time, normally, a function higher than the cubic spline is used for interpolation, and the linear is considered to be largely unsatisfactory.
そして、このようにして画像歪み補正装置４５により得られた補正後の２つの各画像データは相互相関計算装置４７および信号強度比較装置４８にそれぞれ与えられる。
The two pieces of corrected image data obtained by the image
次に、ステップＳ５において、相互相関計算装置４７は、傾斜磁場を反転させた２つの撮影条件によって得られた補正後の画像データ間の相互相関係数を求め、求められた相互相関係数に基づいて画像データの補正が適切に行われたか否かを評価する。
Next, in step S5, the
つまり、２つの画像データの収集条件は、傾斜磁場の極性の向きが互いに反転していることを除けば実質的に同一であるため、画像データの歪みの補正が適切に行われている程、より各画像データ上の物体の形状が一致し、２つの画像データ間における相互相関が強くなって相互相関係数が大きくなると考えられる。逆に、画像データの歪みの補正が適切に行われていない場合には、２つの画像データ間における相互相関が弱くなって相互相関係数が小さくなると考えられる。 In other words, the two image data collection conditions are substantially the same except that the polar directions of the gradient magnetic fields are reversed from each other. Therefore, the more appropriately the distortion of the image data is corrected, It is considered that the shape of the object on each image data is more consistent, the cross-correlation between the two image data becomes stronger, and the cross-correlation coefficient becomes larger. On the contrary, if the distortion of the image data is not properly corrected, it is considered that the cross-correlation coefficient between the two image data becomes weak and the cross-correlation coefficient becomes small.
そこで、相互相関計算装置４７は、２つの画像データ間の相互相関係数を予め設定された閾値と比較し、相互相関係数が閾値以上または閾値を超える値である場合には補正が適切に行われたと評価する一方、相互相関係数が閾値未満または閾値以下である場合には補正が適切に行われていないと評価する。
Therefore, the
また、必要に応じて相互相関計算装置４７は、求められた２つの画像データ間の相互相関係数や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えて表示させる。このため、ユーザは、画像データの補正が適切に行われたか否かを容易に把握することができる。
Further, if necessary, the
また、ステップＳ６において、信号強度比較装置４８は、傾斜磁場を反転させた２つの撮影条件によって得られた補正後の画像データ間でピクセルごとに信号強度を比較し、信号強度の比較結果にに基づいて画像データの補正が適切に行われたか否かを評価する。
In step S6, the signal
つまり、前述のように２つの画像データの収集条件は、傾斜磁場の極性の向きが互いに反転していることを除けば実質的に同一であるため、画像データの歪みの補正が適切に行われている程、より各画像データ上の物体の形状が一致し、２つの画像データ間における信号強度の差分が小さくなると考えられる。逆に、画像データの歪みの補正が適切に行われていない場合には、２つの画像データ間における信号強度の差分が大きくなると考えられる。 That is, as described above, the collection conditions of the two image data are substantially the same except that the directions of the polarities of the gradient magnetic fields are reversed, so that the distortion of the image data is appropriately corrected. The closer the shape of the object on each image data, the smaller the difference in signal intensity between the two image data. Conversely, if the distortion of the image data is not properly corrected, the difference in signal strength between the two image data is considered to be large.
そこで、信号強度比較装置４８は、２つの画像データ間の信号強度の差分値を予め設定された閾値と比較し、信号強度の差分値が閾値以上または閾値を超える値である場合には補正が適切に行われていないと評価する一方、信号強度の差分値が閾値未満または閾値以下である場合には補正が適切に行われたと評価する。
Therefore, the signal
また、必要に応じて信号強度比較装置４８は、求められた２つの画像データ間の信号強度の差分値や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えて表示させる。このため、ユーザは、画像データの補正が適切に行われたか否かを容易に把握することができる。
Further, the signal
つまり以上のような磁気共鳴イメージング装置２０は、傾斜磁場を反転させて２つの画像データを作成し、歪み補正後の２つの画像データの相関係数や信号強度の比較結果に基づいて補正の適否を評価するものである。
That is, the magnetic
このため、磁気共鳴イメージング装置２０によれば、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することができる。
Therefore, according to the magnetic
尚、磁気共鳴イメージング装置２０の一部の機能や処理を省略してもよく、処理順序を変更してもよい。また、撮影条件は２通りに限らず、例えば、基準となる撮影条件、ＳＥ方向のみの傾斜磁場を反転させた撮影条件、ＲＯ方向のみの傾斜磁場を反転させた撮影条件の３通りの撮影条件とする場合のように、３つ以上の撮影条件としてもよい。
Note that some functions and processing of the magnetic
２０ 磁気共鳴イメージング装置
２１ 静磁場用磁石
２２ シムコイル
２３ 傾斜磁場コイルユニット
２４ ＲＦコイル
２５ 制御系
２６ 静磁場電源
２７ 傾斜磁場電源
２８ シムコイル電源
２９ 送信器
３０ 受信器
３１ シーケンスコントローラ
３２ コンピュータ
３３ 入力装置
３４ 表示装置
３５ 演算装置
３６ 記憶装置
３７ 寝台
４０ シーケンスコントローラ制御装置
４１ 生データデータベース
４２ 画像再構成装置
４３ 画像データデータベース
４４ 静磁場分布データベース
４５ 画像歪み補正装置
４６ 撮影条件設定装置
４７ 相互相関計算装置
４８ 信号強度比較装置
４９ 磁場分布計算装置
５０ 傾斜磁場測定コイル
Ｐ 被検体
20 Magnetic
Claims (4)
前記撮影条件設定手段により設定された前記各撮影条件に従って撮影領域に高周波信号を送信することによる原子核の核磁気共鳴に伴って発生した磁気共鳴信号に基づいてデータを生成するデータ収集手段と、
前記データに画像再構成処理を施すことにより前記撮影条件ごとに画像データをそれぞれ再構成する画像再構成手段と、
前記撮影領域の静磁場分布に基づいて前記各撮影条件により得られた前記各画像データの歪みをそれぞれ補正する画像歪み補正手段と、
前記画像歪み補正手段による補正後における前記各画像データに基づいて前記各画像データの補正が適切に行われたか否かを評価する画像補正評価手段と、
を有することを特徴とする磁気共鳴イメージング装置。 An imaging condition setting means for setting a plurality of imaging conditions so that the scan is executed by inverting the polarity of the gradient magnetic field in at least one direction of the phase encoding direction, the readout direction, and the slice encoding direction;
Data collecting means for generating data based on magnetic resonance signals generated along with nuclear magnetic resonance of the nuclei by transmitting a high-frequency signal to the imaging region according to each imaging condition set by the imaging condition setting means;
Image reconstruction means for reconstructing image data for each of the photographing conditions by performing image reconstruction processing on the data;
Image distortion correction means for correcting distortion of each image data obtained by each imaging condition based on the static magnetic field distribution of the imaging region;
Image correction evaluation means for evaluating whether or not the correction of each image data is appropriately performed based on each image data after correction by the image distortion correction means;
A magnetic resonance imaging apparatus comprising:
前記画像データの撮影領域における静磁場分布に基づいて前記各撮影条件により得られた前記各画像データの歪みをそれぞれ補正するステップと、
補正後における前記各画像データに基づいて前記各画像データの補正が適切に行われたか否かを評価するステップと、
を有することを特徴とする画像補正評価方法。 The imaging condition is obtained by performing image reconstruction processing on data collected by executing a scan according to a plurality of imaging conditions in which the polarity of the gradient magnetic field in at least one of the phase encoding direction, the readout direction, and the slice encoding direction is reversed Each reconstructing the image data for each,
Correcting each distortion of the image data obtained by the imaging conditions based on the static magnetic field distribution in the imaging area of the image data;
Evaluating whether each of the image data is appropriately corrected based on each of the image data after correction; and
An image correction evaluation method characterized by comprising:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004309971A JP4625677B2 (en) | 2004-10-25 | 2004-10-25 | Magnetic resonance imaging apparatus and image correction evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004309971A JP4625677B2 (en) | 2004-10-25 | 2004-10-25 | Magnetic resonance imaging apparatus and image correction evaluation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2006116216A JP2006116216A (en) | 2006-05-11 |
| JP4625677B2 true JP4625677B2 (en) | 2011-02-02 |
Family
ID=36534644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2004309971A Expired - Fee Related JP4625677B2 (en) | 2004-10-25 | 2004-10-25 | Magnetic resonance imaging apparatus and image correction evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4625677B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5181235B2 (en) * | 2008-02-25 | 2013-04-10 | 株式会社日立メディコ | Image processing apparatus, medical image photographing apparatus, program, and image processing method |
| WO2013161910A1 (en) | 2012-04-24 | 2013-10-31 | 株式会社東芝 | Pet-mri device |
| CN115079073B (en) * | 2022-03-10 | 2023-04-11 | 杭州永川科技有限公司 | Frequency difference quasi-static magnetic induction imaging method, system, equipment and medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61198043A (en) * | 1985-02-28 | 1986-09-02 | Shimadzu Corp | Static magnetic field variation detecting method for mr-ct device |
| JPS62176443A (en) * | 1986-01-29 | 1987-08-03 | 横河メディカルシステム株式会社 | Method for correcting position and concentration strain of nmr image |
| JPH09220213A (en) * | 1996-02-19 | 1997-08-26 | Siemens Ag | Functional image generating method and nuclear spin tomography system |
| JP2001276014A (en) * | 2000-03-30 | 2001-10-09 | Toshiba Corp | Mri instrument and functional mr imaging method |
| JP2003116815A (en) * | 2001-10-15 | 2003-04-22 | Hitachi Medical Corp | Magnetic resonance imaging system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6456042A (en) * | 1987-08-26 | 1989-03-02 | Toshiba Corp | Magnetic resonance imaging apparatus |
-
2004
- 2004-10-25 JP JP2004309971A patent/JP4625677B2/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61198043A (en) * | 1985-02-28 | 1986-09-02 | Shimadzu Corp | Static magnetic field variation detecting method for mr-ct device |
| JPS62176443A (en) * | 1986-01-29 | 1987-08-03 | 横河メディカルシステム株式会社 | Method for correcting position and concentration strain of nmr image |
| JPH09220213A (en) * | 1996-02-19 | 1997-08-26 | Siemens Ag | Functional image generating method and nuclear spin tomography system |
| JP2001276014A (en) * | 2000-03-30 | 2001-10-09 | Toshiba Corp | Mri instrument and functional mr imaging method |
| JP2003116815A (en) * | 2001-10-15 | 2003-04-22 | Hitachi Medical Corp | Magnetic resonance imaging system |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2006116216A (en) | 2006-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4812420B2 (en) | Magnetic resonance imaging apparatus and image correction evaluation method | |
| US7535227B1 (en) | Method and apparatus for correcting distortion in MR images caused by metallic implants | |
| CN102711602B (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US8274286B2 (en) | System and method for multi-spectral MR imaging near metal | |
| JP4619674B2 (en) | Magnetic resonance imaging system | |
| US8093894B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method for improving uniformity in sensitivity map | |
| US20120074940A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US7689015B2 (en) | Magnetic resonance imaging apparatus and image correction estimating method | |
| CN104323775A (en) | Magnetic resonance imaging apparatus, image processing apparatus, and image processing method | |
| WO2010035569A1 (en) | Magnetic resonance imaging device | |
| JP2004073538A (en) | Magnetic resonance imaging apparatus | |
| JP2012040362A (en) | Magnetic resonance imaging method, magnetic resonance imaging apparatus, and control device of magnetic resonance imaging apparatus | |
| JP5666779B2 (en) | Magnetic resonance imaging system | |
| US6469505B1 (en) | Method and apparatus to reduce perturbation field effects in MR images by restricting the region of interest | |
| US9846216B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US20150168521A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| JP5259177B2 (en) | Magnetic resonance imaging system | |
| US11675038B2 (en) | Diffusion MR imaging with fat suppression | |
| JP4625677B2 (en) | Magnetic resonance imaging apparatus and image correction evaluation method | |
| US20240183924A1 (en) | Distortion-free diffusion and quantitative magnetic resonance imaging with blip up-down acquisition of spin- and gradient-echoes | |
| Andersson | Geometric distortions in diffusion MRI | |
| JP3705995B2 (en) | Gradient coil manufacturing method, gradient coil, and magnetic resonance imaging apparatus | |
| US9916648B1 (en) | Method for correcting a diffusion image having an artifact | |
| US11402454B2 (en) | Correction of distorted diffusion-weighted magnetic resonance image data | |
| JP4040745B2 (en) | MR device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071025 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100401 |
|
| 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: 20101012 |
|
| 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: 20101108 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4625677 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131112 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| 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 |