JP2009284958A - Magnetic resonance imaging system and phantom - Google Patents

Magnetic resonance imaging system and phantom Download PDF

Info

Publication number
JP2009284958A
JP2009284958A JP2008138295A JP2008138295A JP2009284958A JP 2009284958 A JP2009284958 A JP 2009284958A JP 2008138295 A JP2008138295 A JP 2008138295A JP 2008138295 A JP2008138295 A JP 2008138295A JP 2009284958 A JP2009284958 A JP 2009284958A
Authority
JP
Japan
Prior art keywords
phantom
magnetic field
gradient magnetic
image
resonance imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2008138295A
Other languages
Japanese (ja)
Inventor
Naoya Sakaguchi
直哉 坂口
Hidehisa Akimaru
英久 秋丸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2008138295A priority Critical patent/JP2009284958A/en
Publication of JP2009284958A publication Critical patent/JP2009284958A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To realize a single phantom adjustable to a plurality of kinds without being accompanied by the alteration of an arranging position in an MRI system. <P>SOLUTION: The phantom 26 has a cubic base and all of the side parts of the regular hexahedron of the cube-shaped base are formed into a curved face shape while all of the corner parts of the regular hexahedron are also formed into the curved face. The signal values of the sharpened parts (corner parts) being the shape of the regular hexahedron become high but, if the corner parts are formed into the curved surface shape like the phantom 26, the signal values of them can be prevented from becoming high and the curving of the contour of an imaged can be suppressed. Since the distance between mutually parallel sides are kept constant if the concentration of the signal values is not considered at the time of processing of the image of the phantom 26 and the image is normal, the image can be processed easily and rapidly. Since the center slice images in X-, Y- and Z-axis directions are mutually equal in the phantom 26, image processing can be performed by uniform processing without altering the installation direction of the phantom 26 at every axis. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、磁気共鳴イメージング装置及び傾斜磁場調整用のファントムに関する。   The present invention relates to a magnetic resonance imaging apparatus and a phantom for adjusting a gradient magnetic field.

磁気共鳴イメージング装置(MRI装置)は、均一な静磁場中に置かれた被検体に高周波磁場を印加することによって、被検体の任意の領域に存在する原子核(通常プロトン)に核磁気共鳴現象を生じさせ、それによって発生する核磁気共鳴信号(エコー信号)からその領域の断層像を得るものである。この際、傾斜磁場のゲイン、スリューレート、渦電流補正等の特性値の調整が正確に行う必要がある。   A magnetic resonance imaging apparatus (MRI apparatus) applies a high-frequency magnetic field to a subject placed in a uniform static magnetic field, thereby causing a nuclear magnetic resonance phenomenon to occur in nuclei (usually protons) present in an arbitrary region of the subject. A tomographic image of the region is obtained from a nuclear magnetic resonance signal (echo signal) generated thereby. At this time, it is necessary to accurately adjust characteristic values such as the gradient magnetic field gain, slew rate, and eddy current correction.

MRI装置における各調整値は装置ごとに異なるため、据付終了後に正確な調整値を設定する必要がある。このため、特許文献1に記載されているようなファントムを用いて、調整が行われている。   Since each adjustment value in the MRI apparatus is different for each apparatus, it is necessary to set an accurate adjustment value after the installation is completed. For this reason, adjustment is performed using a phantom as described in Patent Document 1.

特開平4−26411号公報JP-A-4-26411

従来技術においては、傾斜磁場のゲイン、スリューレート等のそれぞれの調整に必要なファントムおよび受信コイルは、それぞれ異なった姿勢や異なったものが必要であり、調整対象物が変わる毎に、作業者によるファントム及び受信コイルの交換作業が必要である。   In the prior art, the phantom and the receiving coil required for each adjustment of the gradient magnetic field gain, slew rate, etc. require different postures and different ones. It is necessary to replace the phantom and the receiving coil.

例えば、傾斜磁場ゲイン調整では、専用のファントムを用いて調整を行うが、X、Y、Z軸毎にファントムの設置方向を変更しながら調整しなければならなかった。また、傾斜磁場スリューレート調整時には、傾斜磁場ゲイン調整で使用したファントムではないファントムを使用する必要があった。   For example, in the gradient magnetic field gain adjustment, adjustment is performed using a dedicated phantom, but the adjustment must be performed while changing the installation direction of the phantom for each of the X, Y, and Z axes. In addition, when adjusting the gradient magnetic field slew rate, it is necessary to use a phantom other than the phantom used for adjusting the gradient magnetic field gain.

傾斜磁場のゲイン等の調整時において、作業者は、ファントムの設置方向を変更したり、調整ごとに別のファントムに変更しなければならず、非効率であった。   When adjusting the gain of the gradient magnetic field, the operator must change the installation direction of the phantom or change to another phantom for each adjustment, which is inefficient.

本発明の目的は、配置位置の変更を伴うことなく、複数種類の特性値の調整が可能な単一のファントムと、そのファントムを用いて傾斜磁場の調整が可能な磁気共鳴イメージング装置を実現することである。   An object of the present invention is to realize a single phantom capable of adjusting a plurality of types of characteristic values without changing the arrangement position, and a magnetic resonance imaging apparatus capable of adjusting a gradient magnetic field using the phantom. That is.

上記目的を達成するため、本発明は次のように構成される。   In order to achieve the above object, the present invention is configured as follows.

磁気共鳴イメージング装置において、正六面体の角部及び辺が曲面状のファントムのスライス画像を撮像し、撮像したスライス画像から得られたファントムの寸法情報と、ファントムの実際の寸法情報とに基づいて、傾斜磁場発生手段の特性値を調整する。   In the magnetic resonance imaging apparatus, a sliced image of a phantom having a regular hexahedron with curved corners and sides is taken, and based on the phantom dimensional information obtained from the taken slice image and the actual phantom dimensional information, The characteristic value of the gradient magnetic field generating means is adjusted.

また、磁気共鳴イメージング装置において、正六面体の角部及び辺が曲面状のファントムの拡散傾斜磁場パルスを印加したときのスライス画像と、拡散傾斜磁場パルスを印加していないときのスライス画像とを撮像させ、撮像した2つのスライス画像から得られたファントムの寸法情報を互いに比較した結果に基づいて、傾斜磁場発生手段の特性値を調整する。   In addition, in a magnetic resonance imaging apparatus, a slice image obtained when a diffusion gradient magnetic field pulse of a phantom whose corners and sides of a regular hexahedron are curved and a slice image obtained when no diffusion gradient magnetic field pulse is applied are taken. The characteristic value of the gradient magnetic field generating means is adjusted based on the result of comparing the phantom dimensional information obtained from the two captured slice images.

配置位置の変更を伴うことなく、複数種類の特性値の調整が可能な単一のファントムと、そのファントムを用いて傾斜磁場の調整が可能な磁気共鳴イメージング装置を実現することができる。   A single phantom capable of adjusting a plurality of types of characteristic values without changing the arrangement position and a magnetic resonance imaging apparatus capable of adjusting a gradient magnetic field using the phantom can be realized.

以下、本発明の実施形態について添付図面を参照して説明する。   Embodiments of the present invention will be described below with reference to the accompanying drawings.

図1は本発明が適用されるMRI装置の概略構成図である。図1において、磁気共鳴イメージング装置は、核磁気共鳴(NMR)現象を利用して被検体の断層画像を得るもので、静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、中央処理装置(CPU)8とを備える。   FIG. 1 is a schematic configuration diagram of an MRI apparatus to which the present invention is applied. In FIG. 1, a magnetic resonance imaging apparatus obtains a tomographic image of a subject using a nuclear magnetic resonance (NMR) phenomenon, and includes a static magnetic field generation system 2, a gradient magnetic field generation system 3, a transmission system 5, A reception system 6, a signal processing system 7, a sequencer 4, and a central processing unit (CPU) 8 are provided.

静磁場発生系2は、被検体1の周りの空間にその体軸方向または体軸と直交する方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式または常電導方式あるいは超電導方式の磁場発生手段が配置されている。傾斜磁場発生系3は、X,Y,Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とを備える。傾斜磁場電源10は、後述するシーケンサ4からの命令に従ってそれぞれのコイルを駆動することにより、X,Y,Zの3軸方向の傾斜磁場Gs,Gp,Gfを被検体1に印加する。   The static magnetic field generation system 2 generates a uniform static magnetic field in a body axis direction or a direction perpendicular to the body axis in a space around the subject 1, and a permanent magnet system or a normal conduction system around the subject 1. Alternatively, a superconducting magnetic field generating means is arranged. The gradient magnetic field generation system 3 includes a gradient magnetic field coil 9 wound in three axial directions of X, Y, and Z, and a gradient magnetic field power source 10 that drives each gradient magnetic field coil. The gradient magnetic field power supply 10 applies the gradient magnetic fields Gs, Gp, and Gf in the three-axis directions of X, Y, and Z to the subject 1 by driving each coil in accordance with a command from the sequencer 4 described later.

より具体的には、X,Y,Zのいずれかの1方向にスライス方向傾斜磁場パルス(Gs)を印加して被検体1に対するスライス面を設定し、残り2つの方向に位相エンコード方向傾斜磁場パルス(Gp)と、周波数エンコード方向傾斜磁場パルス(Gf)を印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。   More specifically, a slice direction gradient magnetic field pulse (Gs) is applied in one of X, Y, and Z to set a slice plane for the subject 1, and the phase encode direction gradient magnetic field is applied in the remaining two directions. A pulse (Gp) and a frequency encoding direction gradient magnetic field pulse (Gf) are applied, and position information in each direction is encoded in the echo signal.

シーケンサ4は、高周波磁場パルス(RFパルス)と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段である。このシーケンサ4は、CP∪8の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を、送信系5、傾斜磁場発生系3、および受信系6に送る。   The sequencer 4 is a control unit that repeatedly applies a high-frequency magnetic field pulse (RF pulse) and a gradient magnetic field pulse in a predetermined pulse sequence. The sequencer 4 operates under the control of the CP 8 and sends various commands necessary for collecting tomographic image data of the subject 1 to the transmission system 5, the gradient magnetic field generation system 3, and the reception system 6.

送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるためにRFパルスを照射するもので、高周波発振器11と、変調器12と、高周波増幅器13と、送信側の高周波コイル14aとを備える。   The transmission system 5 irradiates an RF pulse to cause nuclear magnetic resonance to the atomic nuclear spin constituting the biological tissue of the subject 1, and includes a high-frequency oscillator 11, a modulator 12, a high-frequency amplifier 13, A high-frequency coil 14a on the transmission side.

高周波発振器11から出力された高周波パルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給する。これにより、電磁波(RFパルス)が被検体1に照射される。   The high-frequency pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 at a timing according to a command from the sequencer 4, and after the amplitude-modulated high-frequency pulse is amplified by the high-frequency amplifier 13, the high-frequency pulse is arranged close to the subject 1. Is supplied to the high frequency coil 14a. Thereby, the subject 1 is irradiated with electromagnetic waves (RF pulses).

受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号(NMR信号)を検出する。この受信系6は、受信側の高周波コイル14bと、増幅器15と、直交位相検波器16と、A/D変換器17とを備える。送信側の高周波コイル14aから照射された電磁波によって誘起される被検体1の応答の電磁波(NMR信号)が被検体1に近接して配置された高周波コイル14bで検出される。そして、増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。   The receiving system 6 detects an echo signal (NMR signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the subject 1. The reception system 6 includes a reception-side high frequency coil 14 b, an amplifier 15, a quadrature phase detector 16, and an A / D converter 17. The response electromagnetic wave (NMR signal) of the subject 1 induced by the electromagnetic wave irradiated from the transmission-side high frequency coil 14 a is detected by the high frequency coil 14 b disposed in the vicinity of the subject 1. Then, after being amplified by the amplifier 15, it is divided into two orthogonal signals by the quadrature phase detector 16 at a timing according to a command from the sequencer 4, and each is converted into a digital quantity by the A / D converter 17, It is sent to the signal processing system 7.

信号処理系7は、光ディスク19、磁気ディスク18等の外部記憶装置と、CRT等からなるディスプレイ20とを有する。そして、信号処理系7は、受信系6からのデータがCPU8に入力されると、CPU8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。   The signal processing system 7 includes an external storage device such as an optical disk 19 and a magnetic disk 18 and a display 20 made up of a CRT or the like. When the data from the reception system 6 is input to the CPU 8, the signal processing system 7 performs processing such as signal processing and image reconstruction, and displays the tomographic image of the subject 1 as a result of the display 20. And recorded on the magnetic disk 18 or the like of the external storage device.

操作部25は、MRI装置の各種制御情報や信号処理系7で行う処理の制御情報を入力するもので、トラックボール又はマウス23、及び、キーボード24から構成される。この操作部25はディスプレイ20に近接して配置され、操作者がディスブレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。   The operation unit 25 is used to input various control information of the MRI apparatus and control information of processing performed by the signal processing system 7 and includes a trackball or mouse 23 and a keyboard 24. The operation unit 25 is disposed in the vicinity of the display 20, and the operator controls various processes of the MRI apparatus interactively through the operation unit 25 while looking at the display 20.

なお、図1において、送信側及び受信側の高周波コイル14a、14bと傾斜磁場コイル9は、被検体1の周りの空間に配置された静磁場発生系2の静磁場空間内に設置されている。   In FIG. 1, the high-frequency coils 14 a and 14 b and the gradient magnetic field coil 9 on the transmission side and the reception side are installed in the static magnetic field space of the static magnetic field generation system 2 arranged in the space around the subject 1. .

次に、本発明の第1の実施形態であるファントムについて説明する。   Next, the phantom which is the 1st Embodiment of this invention is demonstrated.

図2は、本発明の第1の実施形態であるファントム26の概観図である。図2に示すように、ファントム26は、立方体型ベースの形状となっている。つまり、正六面体の全ての辺部分が曲面状となっていると共に、正六面体の全ての角部分も曲面状となっている。   FIG. 2 is an overview of the phantom 26 according to the first embodiment of the present invention. As shown in FIG. 2, the phantom 26 has a cubic base shape. That is, all the sides of the regular hexahedron are curved, and all the corners of the regular hexahedron are curved.

図3は、ファントム26の中心スライスの画像を示す図である。図3に示すように、ファントム26の中心スライス画像は、正方形の角を丸くした形になり、互いに平行な辺を持つ画像となる。また、中心スライスであれば、どの断面でも同じ画像となる。   FIG. 3 is a diagram showing an image of the center slice of the phantom 26. As shown in FIG. 3, the central slice image of the phantom 26 has a shape in which square corners are rounded, and is an image having sides parallel to each other. Moreover, if it is a center slice, the same image is obtained in any cross section.

正六面体のままの形状であると、尖った部分(正六面体の角部)の信号値は高くなるが、ファントム26のように、角部が丸くなっていることにより(角部が曲面状となっていることにより)、信号値が高くなることを抑制でき、画像の輪郭が湾曲することを抑えることができる。   If the shape is a regular hexahedron, the signal value of the pointed part (corner part of the regular hexahedron) becomes high, but the corner part is rounded like the phantom 26 (the corner part has a curved surface shape). As a result, the signal value can be prevented from becoming high and the contour of the image can be prevented from being curved.

これにより、ファントム26の画像処理時に、信号値の集中を考慮せず、正常な画像であれば互いに平行な辺の間の距離は一定に保たれるため、容易に、かつ迅速に処理できる。また、ファントム26は、X、Y、Z軸方向の中心スライス画像が互いに同じであることから、軸ごとにファントムの設置方向を変更することなく、一律した処理で画像処理を行うことが可能である。   As a result, when the image processing of the phantom 26 is performed, the concentration of signal values is not taken into consideration and the distance between the sides parallel to each other is kept constant if the image is a normal image. Further, since the central slice images in the X, Y, and Z axis directions of the phantom 26 are the same, it is possible to perform image processing with uniform processing without changing the installation direction of the phantom for each axis. is there.

次に、本発明の第1の実施形態において、ファントム26を用いて傾斜磁場ゲインを自動で調整する動作例を説明する。   Next, an operation example in which the gradient magnetic field gain is automatically adjusted using the phantom 26 in the first embodiment of the present invention will be described.

最適な傾斜磁場ゲイン(Goptimal)は、撮像時に使用した傾斜磁場ゲイン(Gscan)、上記ファントムの実サイズ(Sactual)、画像再構成によって出力された画像上の上記ファントムのサイズ(Simage)を用いて、次式(1)で算出される。 The optimum gradient magnetic field gain (G optimal ) includes the gradient magnetic field gain (G scan ) used at the time of imaging, the actual size of the phantom (S actual ), and the size of the phantom on the image output by image reconstruction (S image). ) And is calculated by the following equation (1).

optimal=(Sactual/Simage)・Gscan ・・・(1)
上記式(1)で表される原理を用いて、上記ファントム画像を解析し、ファントム26のサイズSimageを算出することにより、最適な傾斜磁場ゲインを自動算出する。つまり、スライス画像から得られたファントム26の寸法で、ファントム26の実際の寸法を除算し、その値と撮影時に使用した傾斜磁場ゲインを乗算して最適な傾斜磁場ゲインを算出し、傾斜磁場発生手段の特性値を調整する。
G optimal = (S actual / S image ) · G scan (1)
Using the principle expressed by the above formula (1), the phantom image is analyzed, and the size S image of the phantom 26 is calculated to automatically calculate the optimum gradient magnetic field gain. That is, the actual size of the phantom 26 is divided by the size of the phantom 26 obtained from the slice image, and the optimum gradient magnetic field gain is calculated by multiplying the actual value of the phantom 26 and the gradient magnetic field gain used at the time of photographing. Adjust the characteristic value of the means.

次に、画像解析方法例を説明する。図4に示すように、上記ファントム26の画像上の中央部に互いに平行なプロファイルライン31、プロファイルライン32、プロファイルライン33を引き、このプロファイルライン31、32、33上の信号値から、図5に示すように、ノイズ領域34とファントム領域35とに分割する。   Next, an example of an image analysis method will be described. As shown in FIG. 4, a profile line 31, a profile line 32, and a profile line 33 that are parallel to each other are drawn at the center of the image of the phantom 26, and the signal values on the profile lines 31, 32, and 33 are used as shown in FIG. As shown in FIG. 4, the noise area 34 and the phantom area 35 are divided.

上記処理により、3本のプロファイルライン31、32、33上のファントム領域の端点座標を獲得し、その座標から図6に示す三角形(A)36と三角形(B)37のデータを抽出する。三角形36は、プロファイルライン31、32、33に直交し、プロファイルライン33とファントム26の上面との交点からプロファイルライン31まで延びる直線と、ファントム26の上面とプロファイルライン33との交点から、ファントム26の上面とプロファイルライン31との交点まで延びる直線と、上記2つの直線のそれぞれがプロファイルライン31と交差する2つの交点との間の直線とにより形成される。   By the above processing, the end point coordinates of the phantom area on the three profile lines 31, 32, 33 are obtained, and the data of the triangle (A) 36 and the triangle (B) 37 shown in FIG. 6 are extracted from the coordinates. The triangle 36 is orthogonal to the profile lines 31, 32, 33, and extends from the intersection of the profile line 33 and the upper surface of the phantom 26 to the profile line 31, and from the intersection of the upper surface of the phantom 26 and the profile line 33, the phantom 26. Are formed by a straight line extending to the intersection of the upper surface of the profile line 31 and the profile line 31 and a straight line between the two intersections where each of the two straight lines intersects the profile line 31.

また、三角形37は、画像の中心Oから、ファントム26の下面とプロファイルライン32との交点まで延びる直線と、ファントム26の下面と直交し、画像の中心Oからファントム26の下面まで延びる直線と、上記2つの直線のそれぞれが、ファントム26の下面と交差する2つの交点との間の直線とにより形成される。   The triangle 37 includes a straight line extending from the center O of the image to the intersection of the lower surface of the phantom 26 and the profile line 32, a straight line orthogonal to the lower surface of the phantom 26 and extending from the center O of the image to the lower surface of the phantom 26, Each of the two straight lines is formed by a straight line between two intersections intersecting the lower surface of the phantom 26.

三角形36、37が互いに相似であることを利用し、次式(2)によりファントム26のサイズ(Simage)を算出し、式(1)に代入することにより最適な傾斜磁場ゲインを得る。つまり、三角形36の一直線の寸法であるプロファイルライン31と33との間の距離2aは、既知であり、他の直線の距離bは、得られた画像から算出でき、三角形36は直角三角形であるので、2辺の長さから残りの辺の長さを算出することができる。そして、ファントム26の下面から中心Oまでの距離(length)を画像から算出し、これらを次式(2)に代入することにより、ファントム26のサイズ(Simage)を算出することができる。 Using the fact that the triangles 36 and 37 are similar to each other, the size (S image ) of the phantom 26 is calculated by the following equation (2), and the optimum gradient magnetic field gain is obtained by substituting it into the equation (1). That is, the distance 2a between the profile lines 31 and 33, which is a dimension of a straight line of the triangle 36, is known, and the distance b of the other straight line can be calculated from the obtained image, and the triangle 36 is a right triangle. Therefore, the length of the remaining sides can be calculated from the length of the two sides. Then, the distance (length) from the lower surface of the phantom 26 to the center O is calculated from the image, and by substituting these into the following equation (2), the size (S image ) of the phantom 26 can be calculated.

image=(4a/√((2a)+b))・length ・・・(2)
上述した例では、画像の垂直方向に対するファントム26のサイズ(Simage)を算出したが、同様の方法でファントム26の水平方向のサイズ(Simage)も算出できる。なお、上記処理により、ファントム26の互いに平行な辺(上面、下面)と、プロファイルライン31及びプロファイルライン33とが交差可能な角度までファントム26の設置時の傾きを許容できる。例えば、プロファイルライン31、プロファイルライン32、プロファイルライン33のそれぞれの間隔が25mmで、ファントム26の互いに平行な辺の長さが100mmの場合、ファントム26設置時の傾きを、三本のプロファイルライン31、32、33と交差可能な最大約16度まで許容できる。
S image = (4a / √ ((2a) 2 + b 2 )) · length (2)
In the example described above, the size (S image ) of the phantom 26 with respect to the vertical direction of the image is calculated, but the horizontal size (S image ) of the phantom 26 can also be calculated in the same manner. Note that, by the above processing, the inclination when the phantom 26 is installed can be allowed to an angle at which the parallel sides (upper surface, lower surface) of the phantom 26 and the profile line 31 and the profile line 33 can intersect. For example, when the interval between the profile line 31, the profile line 32, and the profile line 33 is 25 mm and the length of the phantom 26 parallel to each other is 100 mm, the inclination when the phantom 26 is installed is set to the three profile lines 31. , 32, 33 can be allowed up to about 16 degrees.

なお、上記式(1)、(2)の演算処理は、CPU8により実行される。   Note that the arithmetic processing of the above formulas (1) and (2) is executed by the CPU 8.

次に、本発明の第1の実施形態における、傾斜磁場ゲイン調整動作について図7のフローチャートを用いて説明する。この調整動作は、CPU8の制御指令により実行される。   Next, the gradient magnetic field gain adjustment operation in the first embodiment of the present invention will be described with reference to the flowchart of FIG. This adjustment operation is executed by a control command of the CPU 8.

先ず、作業者はどの軸を調整するのかを指示し、調整開始を指示する(ステップ1)。次に、X、Y、Z軸のうちの調整する軸から、画像の計測条件を判断し、ファントム26を撮像する(ステップ2)。取得した撮像データから、先に説明した3本のプロファイルライン31、32、33上のファントム領域の端点座標を取得し、ファントム26の画像上のサイズを算出する(ステップ3)。算出されたファントム26の画像上のサイズから最適な傾斜磁場ゲインを算出する(ステップ4)。   First, the operator instructs which axis to adjust, and instructs the start of adjustment (step 1). Next, an image measurement condition is determined from the X, Y, and Z axes to be adjusted, and the phantom 26 is imaged (step 2). The end point coordinates of the phantom area on the three profile lines 31, 32, 33 described above are acquired from the acquired imaging data, and the size of the phantom 26 on the image is calculated (step 3). An optimum gradient magnetic field gain is calculated from the calculated size of the phantom 26 on the image (step 4).

先の調整値が、調整値の取りうる値の範囲内に収まっているかを判断する(ステップ5)。収まっていなければ、ステップ6へ移動し、エラー表示を行い、調整動作を終了する。この場合は、ファントム26が傾斜しすぎていることが考えられる。   It is determined whether the previous adjustment value is within the range of values that can be taken by the adjustment value (step 5). If not, the process moves to step 6 to display an error and the adjustment operation is terminated. In this case, the phantom 26 may be inclined too much.

ステップ5において、調整値が範囲内に収まっていれば、ステップ7へ移動する。そして、作業者に調整値が表示し(ステップ7)、調整値をMRI装置に設定する(ステップ8)。他の軸の調整が必要かどうかをステップ1の指示より判断し、必要であればステップ2へ移動する。必要でなければ調整を終了する(ステップ9)。   In step 5, if the adjustment value is within the range, the process moves to step 7. Then, the adjustment value is displayed to the operator (step 7), and the adjustment value is set in the MRI apparatus (step 8). It is judged from the instruction in Step 1 whether or not other axes need to be adjusted. If necessary, the process moves to Step 2. If it is not necessary, the adjustment is terminated (step 9).

以上のように、本発明の第1の実施形態によれば、ファントム26を正六面体の角部が曲面状となった形状としたので、ファントムの配置位置の変更を伴うことなく、複数種類の特性値の調整が可能な単一のファントムと、そのファントムを用いて傾斜磁場の特性値の調整が可能な磁気共鳴イメージング装置を実現することができる。   As described above, according to the first embodiment of the present invention, since the phantom 26 has a shape in which the corners of the regular hexahedron are curved, a plurality of types of phantoms can be used without changing the arrangement position of the phantom. A single phantom capable of adjusting the characteristic value and a magnetic resonance imaging apparatus capable of adjusting the characteristic value of the gradient magnetic field using the phantom can be realized.

ここで、従来技術においては、傾斜磁場の調整のための作業は、約30分程度必要であったが、本願の第1の実施形態によれば、約4分程度で調整作業を終了することが可能である。   Here, in the prior art, the work for adjusting the gradient magnetic field required about 30 minutes, but according to the first embodiment of the present application, the adjustment work is finished in about 4 minutes. Is possible.

なお、ファントム26の材質は、アクリル樹脂を使用することができる。また、ファントム26の内容物は塩化ニッケルとすることができる。   The material of the phantom 26 can be an acrylic resin. The contents of the phantom 26 can be nickel chloride.

次に、本発明の第2の実施形態について説明する。この第2の実施形態は、上記ファントム26を用いてMPG Pulse(拡散傾斜磁場パルス)印加時の画像の歪み補正値を自動で算出する例である。   Next, a second embodiment of the present invention will be described. The second embodiment is an example in which an image distortion correction value when an MPG pulse (diffusion gradient magnetic field pulse) is applied is automatically calculated using the phantom 26.

本発明以前の技術にあっては、MPGパルス印加時の画像の歪み補正値は、従来から使用されている瓶型ファントムを撮像し、MPGパルスを印加していない画像を基準に、印加時の画像のサイズ・歪みを調整するための補正値である。従来は、作業者がこの調整値を変更しながら撮像を行い、画像を評価することで決定していた。この場合、各軸毎に、瓶型ファントムの位置や種類の変更を行わなければならず、調整に長時間が必要であった。   In the technology prior to the present invention, the distortion correction value of the image at the time of applying the MPG pulse is based on an image in which a conventionally used bottle phantom is imaged and the MPG pulse is not applied. This is a correction value for adjusting the size and distortion of the image. Conventionally, it is determined by an operator taking an image while changing the adjustment value and evaluating the image. In this case, the position and type of the bottle-type phantom must be changed for each axis, and a long time is required for adjustment.

このため、本願発明のファントム26をMPGパルス印加時の画像の歪み補正値算出処理に使用することにより、補正値調整を短時間で処理可能となる。   For this reason, the correction value adjustment can be processed in a short time by using the phantom 26 of the present invention for calculating the distortion correction value of the image when the MPG pulse is applied.

つまり、本発明のファントム26を使用することにより、第1の実施形態と同様のプロファイルラインによる解析を、画像上の縦・横軸に対して実施することで、ファントム26の画像上のサイズを算出する。   In other words, by using the phantom 26 of the present invention, the same profile line analysis as in the first embodiment is performed on the vertical and horizontal axes on the image, so that the size of the phantom 26 on the image can be reduced. calculate.

ファントム26の、MPGパルスを印加していない画像のサイズと、MPGパスルを印加した画像のサイズとを比較し、サイズ誤差を算出する。上記動作を、歪み補正値を変更しながら撮像した画像に対して繰り返し実施し、最もサイズ誤差が小さな補正値を、最適な歪み補正値として採用する。   A size error is calculated by comparing the size of the image of the phantom 26 to which the MPG pulse is not applied and the size of the image to which the MPG pulse is applied. The above operation is repeated for an image captured while changing the distortion correction value, and the correction value with the smallest size error is adopted as the optimum distortion correction value.

次に、本発明の第2の実施形態において、MPGパルス歪み補正値の調整動作について、図8のフローチャートを用いて説明する。   Next, in the second embodiment of the present invention, the MPG pulse distortion correction value adjustment operation will be described with reference to the flowchart of FIG.

先ず、作業者は変更する歪み補正値の範囲と撮像回数を指示し、調整開始を指示する(ステップ11)。MRI装置のCPU8は指示された歪み補正値の範囲と回数から、撮像に使用する歪み補正値を設定し、ファントム26の撮像を行う。(ステップ12、13)。   First, the operator instructs the range of distortion correction values to be changed and the number of times of imaging, and instructs the start of adjustment (step 11). The CPU 8 of the MRI apparatus sets a distortion correction value to be used for imaging from the range and number of instructed distortion correction values, and images the phantom 26. (Steps 12, 13).

取得した撮像データから、MPGパルスを印加していない画像と印加した画像の縦・横軸に対するプロファイルライン上のファントム領域のサイズをそれぞれ取得し、誤差を取得する(ステップ14)。次に、ステップ15において、撮像回数が指示された回数に満たない場合か否か、つまり、傾斜磁場発生手段のゲイン等の調整値を変更して撮像するか否かを判断し、調整値を変更して撮像する場合は、ステップ12へ移動する。ステップ15において、撮像回数が指示された回数を満たし、調整値を変更する必要が無い場合は、ステップ16に移動する。   From the acquired imaging data, the size of the phantom area on the profile line with respect to the vertical and horizontal axes of the image to which the MPG pulse is not applied and the applied image is acquired, and the error is acquired (step 14). Next, in step 15, it is determined whether or not the number of times of imaging is less than the specified number of times, that is, whether or not imaging is performed by changing an adjustment value such as a gain of the gradient magnetic field generating means. If the image is to be changed, the process moves to step 12. If it is determined in step 15 that the number of times of imaging satisfies the specified number and there is no need to change the adjustment value, the process moves to step 16.

MPGパルスを印加していない画像のサイズと、MPGパルスを印加した画像のサイズとの誤差が最も少ない歪み補正値を調整値として決定する(ステップ16)。そして、決定した調整値が妥当か否かを判断し(ステップ17)、妥当でなければ、エラー表示を行い(ステップ18)、処理は終了する。   A distortion correction value with the smallest error between the size of the image to which the MPG pulse is not applied and the size of the image to which the MPG pulse is applied is determined as an adjustment value (step 16). Then, it is determined whether or not the determined adjustment value is valid (step 17). If it is not valid, an error is displayed (step 18), and the process ends.

ステップ17において、調整値が妥当な範囲であると判断した場合は、調整結果を作業者に通知し(ステップ19)、調整値をMRI装置適用し(ステップ20)、調整を終了する。   If it is determined in step 17 that the adjustment value is within a reasonable range, the operator is notified of the adjustment result (step 19), the adjustment value is applied to the MRI apparatus (step 20), and the adjustment is terminated.

なお、その他の構成については、本発明の第1の実施形態と同様となるので、図示及びその説明は省略する。   Since other configurations are the same as those of the first embodiment of the present invention, illustration and description thereof are omitted.

本発明の第2の実施形態によれば、MPGパルス印加時の画像歪補正値算出処理においても、短時間で処理可能となる。   According to the second embodiment of the present invention, the image distortion correction value calculation process at the time of applying the MPG pulse can be processed in a short time.

なお、本発明のファントムの他の形状例として、図9に示すものがある。図9に示したファントム27は、球体の表面に、互いに平行する2つの面の組を3組形成し、一つの組の2つの面に直交し、球体の中心を通過する直線と、他の一つの組の2つの面に直交し、球体の中心を通過する直線とが互いに直交する形状となっている。   FIG. 9 shows another example of the shape of the phantom of the present invention. The phantom 27 shown in FIG. 9 forms three sets of two surfaces parallel to each other on the surface of the sphere, a straight line passing through the center of the sphere perpendicular to the two surfaces of one set, and the other A straight line that is orthogonal to the two surfaces of one set and passes through the center of the sphere is orthogonal to each other.

このような形状のファントム27であっても、ファントム26と同様な効果を得ることができる。なお、このファントム27は、ファントム26と比較して、中心から平面部までの距離を同一とした場合、容積はファントム26に比較してファントム27の方が小となる。   Even with the phantom 27 having such a shape, the same effect as the phantom 26 can be obtained. Note that the phantom 27 has a smaller volume than the phantom 26 compared to the phantom 26 when the distance from the center to the plane portion is the same.

本発明が適用されるMRI装置の概略構成図である。1 is a schematic configuration diagram of an MRI apparatus to which the present invention is applied. 本発明の実施形態におけるファントムの一例の外観図である。It is an external view of an example of the phantom in the embodiment of the present invention. 図2に示したファントムの中心スライスの撮像画像例を示す図である。It is a figure which shows the example of a captured image of the center slice of the phantom shown in FIG. 本発明の第1の実施形態において、画像データ解析の説明図である。In the 1st Embodiment of this invention, it is explanatory drawing of image data analysis. 本発明の第1の実施形態において、プロファイルライン上の信号値の解析結果の一例を示すグラフである。In the 1st Embodiment of this invention, it is a graph which shows an example of the analysis result of the signal value on a profile line. 本発明の第1の実施形態において、ファントムが傾いている場合に画像上のファントムのサイズを算出する方法の説明図である。In the 1st Embodiment of this invention, when the phantom inclines, it is explanatory drawing of the method of calculating the size of the phantom on an image. 本発明の第1の実施形態において、傾斜磁場ゲインの調整動作のフローチャートである。5 is a flowchart of a gradient magnetic field gain adjustment operation in the first embodiment of the present invention. 本発明の第2の実施形態において、MPGパルスの歪み補正値を調整する動作フローチャートである。9 is an operation flowchart for adjusting a distortion correction value of an MPG pulse in the second embodiment of the present invention. 本発明の実施形態におけるファントムの他の例の外観図である。It is an external view of the other example of the phantom in embodiment of this invention.

符号の説明Explanation of symbols

1・・・被検体、2・・・静磁場発生系、3・・・傾斜磁場発生系、4・・・シーケンサ、5・・・送信系、6・・・受信系、7・・・信号処理系、8・・・中央処理装置(CPU)、9・・・傾斜磁場コイル、10・・・傾斜磁場電源、l1・・・高周波発振器、12・・・変調器、13・・・高周波増幅器、14a、14b・・・高周波コイル、15・・・増幅器、16・・・直交位相検波器、17・・・A/D変換器、18・・・磁気ディスク、19・・・光ディスク、20・・・ディスプレイ、23・・・トラックボール(マウス)、24・・・キーボード、26、27・・・ファントム、31〜33・・・プロファイルライン   DESCRIPTION OF SYMBOLS 1 ... Subject, 2 ... Static magnetic field generation system, 3 ... Gradient magnetic field generation system, 4 ... Sequencer, 5 ... Transmission system, 6 ... Reception system, 7 ... Signal Processing system, 8 ... Central processing unit (CPU), 9 ... Gradient magnetic field coil, 10 ... Gradient magnetic field power supply, l1 ... High frequency oscillator, 12 ... Modulator, 13 ... High frequency amplifier 14a, 14b ... high frequency coil, 15 ... amplifier, 16 ... quadrature phase detector, 17 ... A / D converter, 18 ... magnetic disk, 19 ... optical disc, 20. ..Display, 23 ... Track ball (mouse), 24 ... Keyboard, 26, 27 ... Phantom, 31-33 ... Profile line

Claims (7)

静磁場発生手段と、傾斜磁場発生手段と、高周波送受信手段と、この高周波送受信手段が受信した核磁気共鳴信号に基づいて被検体の画像を再構成し、表示する信号処理表示手段と、上記静磁場発生手段、傾斜磁場発生手段、高周波送受信手段及び信号処理表示手段を制御する制御手段とを有する磁気共鳴イメージング装置において、
上記制御手段は、静磁場発生手段、傾斜磁場発生手段、高周波送受信手段及び信号処理表示手段を制御して、正六面体の角部及び辺が曲面状のファントムのスライス画像を撮像させ、撮像した上記スライス画像から得られた上記ファントムの寸法情報と、上記ファントムの実際の寸法情報とに基づいて、上記傾斜磁場発生手段の特性値を調整することを特徴とする磁気共鳴イメージング装置。
Static magnetic field generation means, gradient magnetic field generation means, high frequency transmission / reception means, signal processing display means for reconstructing and displaying an image of a subject based on a nuclear magnetic resonance signal received by the high frequency transmission / reception means, and the static In a magnetic resonance imaging apparatus having a magnetic field generating means, a gradient magnetic field generating means, a high frequency transmitting / receiving means, and a control means for controlling a signal processing display means,
The control means controls the static magnetic field generation means, the gradient magnetic field generation means, the high frequency transmission / reception means, and the signal processing display means to pick up a slice image of a phantom having a regular hexahedron with corners and sides curved, A magnetic resonance imaging apparatus, wherein a characteristic value of the gradient magnetic field generating means is adjusted based on dimensional information of the phantom obtained from a slice image and actual dimensional information of the phantom.
請求項1記載の磁気共鳴イメージング装置において、上記特性値は、上記傾斜磁場発生手段の傾斜磁場ゲインであることを特徴とする磁気共鳴イメージング装置。   2. The magnetic resonance imaging apparatus according to claim 1, wherein the characteristic value is a gradient magnetic field gain of the gradient magnetic field generating means. 請求項2記載の磁気共鳴イメージング装置において、上記制御手段は、上記ファントムの実際の寸法を上記スライス画像から算出された寸法で除算し、得られた値を、上記スライス画像を撮影した傾斜磁場ゲインと乗算して、得られたゲインを上記傾斜磁場発生手段の最適傾斜磁場ゲインとすることを特徴とする磁気共鳴イメージング装置。   3. The magnetic resonance imaging apparatus according to claim 2, wherein the control means divides the actual dimension of the phantom by the dimension calculated from the slice image, and the obtained value is a gradient magnetic field gain obtained by photographing the slice image. The magnetic resonance imaging apparatus is characterized in that the obtained gain is used as an optimum gradient magnetic field gain of the gradient magnetic field generating means. 請求項3記載の磁気共鳴イメージング装置において、上記制御手段は、上記スライス画像の輪郭と、互いに平行な3本のプロファイルラインとの交点位置に基づいて、上記スライス画像から得られる上記ファントムの寸法情報を算出することを特徴とする磁気共鳴イメージング装置。   4. The magnetic resonance imaging apparatus according to claim 3, wherein the control means is dimensional information of the phantom obtained from the slice image based on an intersection position of the contour of the slice image and three profile lines parallel to each other. The magnetic resonance imaging apparatus characterized by calculating. 静磁場発生手段と、傾斜磁場発生手段と、高周波送受信手段と、この高周波送受信手段が受信した核磁気共鳴信号に基づいて被検体の画像を再構成し、表示する信号処理表示手段と、上記静磁場発生手段、傾斜磁場発生手段、高周波送受信手段及び信号処理表示手段を制御する制御手段とを有する磁気共鳴イメージング装置において、
上記制御手段は、静磁場発生手段、傾斜磁場発生手段、高周波送受信手段及び信号処理表示手段を制御して、正六面体の角部及び辺が曲面状のファントムの拡散傾斜磁場パルスを印加したときのスライス画像と、拡散傾斜磁場パルスを印加していないときのスライス画像とを撮像させ、撮像した上記2つのスライス画像から得られた上記ファントムの寸法情報を互いに比較した結果に基づいて、上記傾斜磁場発生手段の特性値を調整することを特徴とする磁気共鳴イメージング装置。
Static magnetic field generation means, gradient magnetic field generation means, high frequency transmission / reception means, signal processing display means for reconstructing and displaying an image of a subject based on a nuclear magnetic resonance signal received by the high frequency transmission / reception means, and the static In a magnetic resonance imaging apparatus having a magnetic field generating means, a gradient magnetic field generating means, a high frequency transmitting / receiving means, and a control means for controlling a signal processing display means,
The control means controls the static magnetic field generation means, the gradient magnetic field generation means, the high frequency transmission / reception means, and the signal processing display means, and applies a diffusion gradient magnetic field pulse of a phantom whose corners and sides of the regular hexahedron are curved. Based on the result of comparing the phantom dimensional information obtained from the two slice images obtained by capturing the slice image and the slice image when the diffusion gradient magnetic field pulse is not applied, the gradient magnetic field A magnetic resonance imaging apparatus characterized by adjusting a characteristic value of a generating means.
磁気共鳴イメージング装置の傾斜磁場発生手段の特性値を調整するためにそのスライス画像が撮影されるファントムにおいて、
正六面体の角部及び辺が曲面状の形状を有することを特徴とするファントム。
In the phantom where the slice image is taken in order to adjust the characteristic value of the gradient magnetic field generating means of the magnetic resonance imaging apparatus,
A phantom characterized in that corners and sides of a regular hexahedron have a curved shape.
磁気共鳴イメージング装置の傾斜磁場発生手段の特性値を調整するためにそのスライス画像が撮影されるファントムにおいて、
球体の表面に、互いに平行する2つの面の組が3組形成され、一つの組の2つの面に直交し、球体の中心を通過する直線と、他の一つの組の2つの面に直交し、球体の中心を通過する直線とが互いに直交する形状であることを特徴とするファントム。
In the phantom where the slice image is taken in order to adjust the characteristic value of the gradient magnetic field generating means of the magnetic resonance imaging apparatus,
Three sets of two surfaces parallel to each other are formed on the surface of the sphere, orthogonal to the two surfaces of one set, perpendicular to the center of the sphere, and orthogonal to the two surfaces of the other set The phantom is characterized in that the straight lines passing through the center of the sphere are orthogonal to each other.
JP2008138295A 2008-05-27 2008-05-27 Magnetic resonance imaging system and phantom Pending JP2009284958A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008138295A JP2009284958A (en) 2008-05-27 2008-05-27 Magnetic resonance imaging system and phantom

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008138295A JP2009284958A (en) 2008-05-27 2008-05-27 Magnetic resonance imaging system and phantom

Publications (1)

Publication Number Publication Date
JP2009284958A true JP2009284958A (en) 2009-12-10

Family

ID=41454764

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008138295A Pending JP2009284958A (en) 2008-05-27 2008-05-27 Magnetic resonance imaging system and phantom

Country Status (1)

Country Link
JP (1) JP2009284958A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101662436B1 (en) * 2015-03-25 2016-10-05 가톨릭대학교 산학협력단 Performance evaluation phantom for magnetic resonance imaging system using ultra high field
CN110301931A (en) * 2019-07-24 2019-10-08 东软医疗系统股份有限公司 Coordinate registration die body, multi-mode imaging system and its method for registering images
CN111132730A (en) * 2017-08-02 2020-05-08 维申Rt有限公司 Calibration method for a patient monitoring system for use with a radiation therapy device
WO2023109064A1 (en) * 2021-12-14 2023-06-22 武汉联影生命科学仪器有限公司 Gradient sensitivity correction method and system

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101662436B1 (en) * 2015-03-25 2016-10-05 가톨릭대학교 산학협력단 Performance evaluation phantom for magnetic resonance imaging system using ultra high field
CN111132730A (en) * 2017-08-02 2020-05-08 维申Rt有限公司 Calibration method for a patient monitoring system for use with a radiation therapy device
CN111132730B (en) * 2017-08-02 2022-07-22 维申Rt有限公司 Calibration method for a patient monitoring system for use with a radiation therapy device
CN110301931A (en) * 2019-07-24 2019-10-08 东软医疗系统股份有限公司 Coordinate registration die body, multi-mode imaging system and its method for registering images
CN110301931B (en) * 2019-07-24 2023-10-31 沈阳智核医疗科技有限公司 Coordinate registration die body, multi-mode imaging system and image registration method thereof
WO2023109064A1 (en) * 2021-12-14 2023-06-22 武汉联影生命科学仪器有限公司 Gradient sensitivity correction method and system

Similar Documents

Publication Publication Date Title
US7768263B2 (en) Magnetic resonance imaging apparatus and method
CN113498479A (en) Correcting magnetic hysteresis in magnetic resonance imaging
JP5337162B2 (en) Magnetic resonance imaging system
JPS5848839A (en) Testing device using nuclear magnetic resonance
JP2009284958A (en) Magnetic resonance imaging system and phantom
WO2009101871A1 (en) Magnetic resonance imaging device, initial state creating method, optimum intensity deciding method, and magnetic resonance imaging method
JP6618988B2 (en) Magnetic resonance imaging apparatus and RF shimming parameter setting method
US9575152B1 (en) Magnetic resonance imaging
JP5352130B2 (en) Magnetic resonance imaging system
US7358726B2 (en) System of generating a magnetic field and MRI system using the system of generating a magnetic field
JP5186698B2 (en) Magnetic resonance imaging apparatus, image reconstruction method, and program
JP2011110131A (en) Magnetic resonance imaging apparatus
JP2004283483A (en) Magnetic resonance imager
JP5280127B2 (en) Magnetic resonance imaging system
JP5758230B2 (en) Magnetic resonance imaging apparatus and inversion RF pulse phase control method
JP6118127B2 (en) Magnetic resonance imaging system
JP5484001B2 (en) Magnetic resonance imaging apparatus and image correction method
JP5942265B2 (en) Magnetic resonance imaging apparatus and RF pulse control method
JP4969933B2 (en) Magnetic resonance imaging system
JP5322614B2 (en) Magnetic resonance imaging system
JPS63109847A (en) Nuclear magnetic resonance imaging apparatus
JP5309379B2 (en) Medical imaging device
JP2016131847A (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2005040315A (en) Multiple rf receiving coil for magnetic resonance imaging device and magnetic resonance imaging device
JP5638324B2 (en) Magnetic resonance imaging apparatus and image correction method