JP6109508B2 - Magnetic resonance imaging system - Google Patents

Magnetic resonance imaging system Download PDF

Info

Publication number
JP6109508B2
JP6109508B2 JP2012188305A JP2012188305A JP6109508B2 JP 6109508 B2 JP6109508 B2 JP 6109508B2 JP 2012188305 A JP2012188305 A JP 2012188305A JP 2012188305 A JP2012188305 A JP 2012188305A JP 6109508 B2 JP6109508 B2 JP 6109508B2
Authority
JP
Japan
Prior art keywords
magnetic field
gradient magnetic
waveform
current
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
Application number
JP2012188305A
Other languages
Japanese (ja)
Other versions
JP2014045775A (en
Inventor
豊 待井
豊 待井
博志 草原
博志 草原
由守 葛西
由守 葛西
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.)
Canon Medical Systems Corp
Original Assignee
Toshiba Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Medical Systems Corp filed Critical Toshiba Medical Systems Corp
Priority to JP2012188305A priority Critical patent/JP6109508B2/en
Publication of JP2014045775A publication Critical patent/JP2014045775A/en
Application granted granted Critical
Publication of JP6109508B2 publication Critical patent/JP6109508B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明の実施形態は、磁気共鳴イメージング(MRI: Magnetic Resonance Imaging)装置に関する。   Embodiments described herein relate generally to a magnetic resonance imaging (MRI) apparatus.

MRI装置は、静磁場中に置かれた被検体の原子核スピンをラーモア周波数の高周波(RF: radio frequency)信号で磁気的に励起し、この励起に伴って発生する磁気共鳴(MR: magnetic resonance)信号から画像を再構成する画像診断装置である。   The MRI apparatus magnetically excites the nuclear spin of the subject placed in a static magnetic field with a radio frequency (RF) signal of Larmor frequency, and generates magnetic resonance (MR) generated by this excitation. An image diagnostic apparatus that reconstructs an image from a signal.

MRI装置の傾斜磁場電源から傾斜磁場コイルに出力される電流の波形は、典型的には直線的に立ち上がり、一定の値となった後に直線的に立ち下がる台形波である。従って、傾斜磁場コイルから撮像領域に印加される傾斜磁場の波形も、傾斜磁場コイルに出力される電流の波形に応じた波形となる。   The waveform of the current output from the gradient magnetic field power supply of the MRI apparatus to the gradient coil is a trapezoidal wave that typically rises linearly and falls linearly after reaching a constant value. Therefore, the waveform of the gradient magnetic field applied from the gradient magnetic field coil to the imaging region is also a waveform corresponding to the waveform of the current output to the gradient magnetic field coil.

エコープラナーイメージング(EPI: echo planar imaging)法等の高速撮像法では、撮像時間の短縮のために読出し用傾斜磁場の立ち上がり中からMRエコー信号のサンプリングを行う場合がある。傾斜磁場の強度が一定の期間にサンプリングを行う場合には、時間的に等間隔にMR信号のサンプリングを行うことで、k空間(位相空間)上において空間的に等間隔のMRデータを収集することができる。   In a high-speed imaging method such as an echo planar imaging (EPI) method, the MR echo signal may be sampled during the rising of the readout gradient magnetic field in order to shorten the imaging time. When sampling is performed in a period in which the gradient magnetic field strength is constant, MR data is sampled at equal intervals in time to collect MR data at equal intervals in the k space (phase space). be able to.

これに対し、読出し用傾斜磁場の立ち上がり中からMRエコー信号のサンプリングを行う場合には、時間的に等間隔にサンプリングを行うとk空間上において不等間隔のMRデータが収集されることとなる。このため、MRデータの収集後において、傾斜磁場電源から傾斜磁場コイルに出力される電流の波形に基づく補間処理を含むリグリッディング(regridding)処理が実行される。具体的には、傾斜磁場コイルに出力される電流値とMRデータのサンプリング間隔の積が一定となるように時間的に不等間隔に収集されたMRデータと等価なk空間データがregridding処理によって時間的に等間隔に収集されたMRデータから生成される。これにより、k空間上において等間隔のk空間データが画像データの生成用に得られる。   On the other hand, when MR echo signals are sampled from when the gradient magnetic field for reading is rising, if sampling is performed at regular intervals in time, MR data at irregular intervals will be collected in the k space. . For this reason, after collecting MR data, a regridding process including an interpolation process based on the waveform of the current output from the gradient magnetic field power source to the gradient magnetic field coil is executed. Specifically, k-space data equivalent to MR data collected at unequal intervals in time so that the product of the current value output to the gradient coil and the sampling interval of MR data is constant is obtained by regridding processing. Generated from MR data collected at regular intervals in time. Thereby, k-space data at equal intervals on the k-space is obtained for generating image data.

特開2010−172383号公報JP 2010-172383 A

上述のような高速撮像法に限らず、MRIでは、よりイメージングに適した波形を有する傾斜磁場の印加が望まれる。   In MRI, not only the high-speed imaging method as described above, it is desired to apply a gradient magnetic field having a waveform more suitable for imaging.

そこで、本発明は、MRイメージングにより適した波形を有する傾斜磁場の印加を行うことが可能な磁気共鳴イメージング装置を提供すること目的とする。   Accordingly, an object of the present invention is to provide a magnetic resonance imaging apparatus capable of applying a gradient magnetic field having a waveform more suitable for MR imaging.

本発明の実施形態に係る磁気共鳴イメージング装置は、イメージング手段及び条件設定手段を備える。イメージング手段は、傾斜磁場コイル及び傾斜磁場電源を含む傾斜磁場発生システムを用いて被検体の磁気共鳴イメージングを行う。条件設定手段は、前記傾斜磁場電源から出力される第1の電流の波形に応じて撮像領域に印加される第1の傾斜磁場の波形に基づいて、前記磁気共鳴イメージング用に前記傾斜磁場電源から出力される第2の電流の波形及び前記撮像領域に印加される第2の傾斜磁場の波形の少なくとも一方を設定又は提示する。また、条件設定手段は、前記傾斜磁場電源から出力される電流の波形と前記撮像領域に印加される傾斜磁場の波形との間における誤差を表す指標が小さくなるように前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される。
また、本発明の実施形態に係る磁気共鳴イメージング装置は、イメージング手段及び条件設定手段を備える。イメージング手段は、傾斜磁場コイル及び傾斜磁場電源を含む傾斜磁場発生システムを用いて被検体の磁気共鳴イメージングを行う。条件設定手段は、前記傾斜磁場電源から出力される第1の電流の波形に応じて撮像領域に印加される第1の傾斜磁場の波形に基づいて、前記磁気共鳴イメージング用に前記傾斜磁場電源から出力される第2の電流の波形及び前記撮像領域に印加される第2の傾斜磁場の波形の少なくとも一方を設定又は提示する。また、条件設定手段は、前記第1の傾斜磁場の波形として、前記第1の電流の波形に基づくシミュレーションによって取得された傾斜磁場の波形を用いて、前記傾斜磁場電源から出力される電流の波形と前記撮像領域に印加される傾斜磁場の波形との間における誤差を表す指標が小さくなるように前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される。
A magnetic resonance imaging apparatus according to an embodiment of the present invention includes an imaging unit and a condition setting unit. The imaging means performs magnetic resonance imaging of a subject using a gradient magnetic field generation system including a gradient magnetic field coil and a gradient magnetic field power supply. The condition setting means uses the gradient magnetic field power supply for the magnetic resonance imaging based on the waveform of the first gradient magnetic field applied to the imaging region in accordance with the waveform of the first current output from the gradient magnetic field power supply. At least one of the waveform of the output second current and the waveform of the second gradient magnetic field applied to the imaging region is set or presented. Further, the condition setting means is configured to reduce the index indicating the error between the waveform of the current output from the gradient magnetic field power supply and the waveform of the gradient magnetic field applied to the imaging region, so that the second current waveform is reduced. And at least one of the waveforms of the second gradient magnetic field.
The magnetic resonance imaging apparatus according to the embodiment of the present invention includes an imaging unit and a condition setting unit. The imaging means performs magnetic resonance imaging of a subject using a gradient magnetic field generation system including a gradient magnetic field coil and a gradient magnetic field power supply. The condition setting means uses the gradient magnetic field power supply for the magnetic resonance imaging based on the waveform of the first gradient magnetic field applied to the imaging region in accordance with the waveform of the first current output from the gradient magnetic field power supply. At least one of the waveform of the output second current and the waveform of the second gradient magnetic field applied to the imaging region is set or presented. Further, the condition setting means uses the gradient magnetic field waveform obtained by the simulation based on the first current waveform as the first gradient magnetic field waveform, and the waveform of the current output from the gradient magnetic field power source. And at least one of the waveform of the second current and the waveform of the second gradient magnetic field is determined so that an index representing an error between the waveform and the waveform of the gradient magnetic field applied to the imaging region is reduced. The

本発明の実施形態に係る磁気共鳴イメージング装置の構成図。1 is a configuration diagram of a magnetic resonance imaging apparatus according to an embodiment of the present invention. 図1に示すコンピュータの機能ブロック図。The functional block diagram of the computer shown in FIG. 図2に示す撮像条件設定部において設定されるEPIシーケンスの一例を示す図。The figure which shows an example of the EPI sequence set in the imaging condition setting part shown in FIG. 図2に示す傾斜磁場波形設定部において第1の電流波形及び第1の傾斜磁場波形に基づいて生成された第2の電流波形及び第2の傾斜磁場波形の一例を示す図。The figure which shows an example of the 2nd current waveform and 2nd gradient magnetic field waveform which were produced | generated based on the 1st current waveform and the 1st gradient magnetic field waveform in the gradient magnetic field waveform setting part shown in FIG. 図2に示す撮像条件設定部において設定される傾斜磁場の波形の取得用のパルスシーケンスの一例を示すシーケンスチャート。The sequence chart which shows an example of the pulse sequence for acquisition of the waveform of the gradient magnetic field set in the imaging condition setting part shown in FIG. 図1に示す傾斜磁場電源及び傾斜磁場コイルを含む傾斜磁場発生システムの詳細構成例を示す回路図。The circuit diagram which shows the detailed structural example of the gradient magnetic field generation system containing the gradient magnetic field power supply and gradient magnetic field coil shown in FIG. 図6に示す傾斜磁場発生システムの等価回路モデルの一例を示す回路図。The circuit diagram which shows an example of the equivalent circuit model of the gradient magnetic field generation system shown in FIG. 図1に示す磁気共鳴イメージング装置により傾斜磁場電源からの出力電流波形と傾斜磁場波形との間における誤差を考慮してイメージングを実行する際の流れを示すフローチャート。3 is a flowchart showing a flow when imaging is performed in consideration of an error between an output current waveform from a gradient magnetic field power supply and a gradient magnetic field waveform by the magnetic resonance imaging apparatus shown in FIG. 1.

本発明の実施形態に係る磁気共鳴イメージング装置について添付図面を参照して説明する。   A magnetic resonance imaging apparatus according to an embodiment of the present invention will be described with reference to the accompanying drawings.

図1は本発明の実施形態に係る磁気共鳴イメージング装置の構成図である。   FIG. 1 is a configuration diagram of a magnetic resonance imaging apparatus according to an embodiment of the present invention.

磁気共鳴イメージング装置20は、静磁場を形成する筒状の静磁場用磁石21、この静磁場用磁石21の内部に設けられたシムコイル22、傾斜磁場コイル23及びRFコイル24を備えている。   The magnetic resonance imaging apparatus 20 includes a cylindrical static magnetic field magnet 21 that forms a static magnetic field, a shim coil 22, a gradient magnetic field coil 23, and an RF coil 24 provided inside the static magnetic field magnet 21.

また、磁気共鳴イメージング装置20には、制御系25が備えられる。制御系25は、静磁場電源26、傾斜磁場電源27、シムコイル電源28、送信器29、受信器30、シーケンスコントローラ31及びコンピュータ32を具備している。制御系25の傾斜磁場電源27は、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zで構成される。また、コンピュータ32には、入力装置33、表示装置34、演算装置35及び記憶装置36が備えられる。   In addition, the magnetic resonance imaging apparatus 20 includes a control system 25. The control system 25 includes a static magnetic field power supply 26, a gradient magnetic field power supply 27, a shim coil power supply 28, a transmitter 29, a receiver 30, a sequence controller 31, and a computer 32. The gradient magnetic field power source 27 of the control system 25 includes an X-axis gradient magnetic field power source 27x, a Y-axis gradient magnetic field power source 27y, and a Z-axis gradient magnetic field power source 27z. In addition, the computer 32 includes an input device 33, a display device 34, an arithmetic device 35, and a storage device 36.

静磁場用磁石21は静磁場電源26と接続され、静磁場電源26から供給された電流により撮像領域に静磁場を形成させる機能を有する。尚、静磁場用磁石21は超伝導コイルで構成される場合が多く、励磁の際に静磁場電源26と接続されて電流が供給されるが、一旦励磁された後は非接続状態とされるのが一般的である。また、静磁場用磁石21を永久磁石で構成し、静磁場電源26が設けられない場合もある。   The static magnetic field magnet 21 is connected to a static magnetic field power supply 26 and has a function of forming a static magnetic field in the imaging region by a current supplied from the static magnetic field power supply 26. In many cases, the static magnetic field magnet 21 is composed of a superconducting coil, and is connected to the static magnetic field power supply 26 at the time of excitation and supplied with current. It is common. In some cases, the static magnetic field magnet 21 is composed of a permanent magnet and the static magnetic field power supply 26 is not provided.

また、静磁場用磁石21の内側には、同軸上に筒状のシムコイル22が設けられる。シムコイル22はシムコイル電源28と接続され、シムコイル電源28からシムコイル22に電流が供給されて静磁場が均一化されるように構成される。   A cylindrical shim coil 22 is coaxially provided inside the static magnetic field magnet 21. The shim coil 22 is connected to the shim coil power supply 28, and is configured such that a current is supplied from the shim coil power supply 28 to the shim coil 22 to make the static magnetic field uniform.

傾斜磁場コイル23は、X軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zで構成され、静磁場用磁石21の内部において筒状に形成される。傾斜磁場コイル23の内側には寝台37が設けられて撮像領域とされ、寝台37には被検体Pがセットされる。RFコイル24にはガントリに内蔵されたRF信号の送受信用の全身用コイル(WBC: whole body coil)や寝台37や被検体P近傍に設けられるRF信号の受信用の局所コイルなどがある。   The gradient magnetic field coil 23 includes an X-axis gradient magnetic field coil 23 x, a Y-axis gradient magnetic field coil 23 y, and a Z-axis gradient magnetic field coil 23 z, and is formed in a cylindrical shape inside the static magnetic field magnet 21. A bed 37 is provided inside the gradient magnetic field coil 23 as an imaging region, and the subject P is set on the bed 37. The RF coil 24 includes a whole body coil (WBC) for transmitting and receiving an RF signal built in the gantry, a bed 37 and a local coil for receiving an RF signal provided near the subject P.

また、傾斜磁場コイル23は、傾斜磁場電源27と接続される。傾斜磁場コイル23のX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zはそれぞれ、傾斜磁場電源27のX軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zと接続される。   The gradient magnetic field coil 23 is connected to a gradient magnetic field power supply 27. The X-axis gradient magnetic field coil 23x, the Y-axis gradient magnetic field coil 23y, and the Z-axis gradient magnetic field coil 23z of the gradient magnetic field coil 23 are respectively an X-axis gradient magnetic field power supply 27x, a Y-axis gradient magnetic field power supply 27y, and a Z-axis gradient magnetic field coil 27z. It is connected to the magnetic field power supply 27z.

そして、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zからそれぞれX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zに供給された電流により、撮像領域にそれぞれX軸方向の傾斜磁場Gx、Y軸方向の傾斜磁場Gy、Z軸方向の傾斜磁場Gzを形成することができるように構成される。   The X-axis gradient magnetic field power source 27x, the Y-axis gradient magnetic field power source 27y, and the Z-axis gradient magnetic field power source 27z are supplied with currents supplied to the X-axis gradient magnetic field coil 23x, the Y-axis gradient magnetic field coil 23y, and the Z-axis gradient magnetic field coil 23z, respectively. In the imaging region, a gradient magnetic field Gx in the X-axis direction, a gradient magnetic field Gy in the Y-axis direction, and a gradient magnetic field Gz in the Z-axis direction can be formed, respectively.

RFコイル24は、送信器29及び受信器30の少なくとも一方と接続される。送信用のRFコイル24は、送信器29からRF信号を受けて被検体Pに送信する機能を有し、受信用のRFコイル24は、被検体P内部の原子核スピンのRF信号による励起に伴って発生したNMR信号を受信して受信器30に与える機能を有する。   The RF coil 24 is connected to at least one of the transmitter 29 and the receiver 30. The transmission RF coil 24 has a function of receiving an RF signal from the transmitter 29 and transmitting it to the subject P, and the reception RF coil 24 is accompanied by excitation of the nuclear spin inside the subject P by the RF signal. It has a function of receiving the NMR signal generated in this way and giving it to the receiver 30.

一方、制御系25のシーケンスコントローラ31は、傾斜磁場電源27、送信器29及び受信器30と接続される。シーケンスコントローラ31は傾斜磁場電源27、送信器29及び受信器30を駆動させるために必要な制御情報、例えば傾斜磁場電源27に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したシーケンス情報を記憶する機能と、記憶した所定のシーケンスに従って傾斜磁場電源27、送信器29及び受信器30を駆動させることによりX軸傾斜磁場Gx、Y軸傾斜磁場Gy,Z軸傾斜磁場Gz及びRF信号を発生させる機能を有する。   On the other hand, the sequence controller 31 of the control system 25 is connected to the gradient magnetic field power supply 27, the transmitter 29 and the receiver 30. The sequence controller 31 is control information necessary for driving the gradient magnetic field power supply 27, the transmitter 29, and the receiver 30, for example, operation control information such as the intensity, application time, and application timing of the pulse current to be applied to the gradient magnetic field power supply 27. And the gradient magnetic field power supply 27, the transmitter 29 and the receiver 30 are driven according to the stored predetermined sequence to drive the X-axis gradient magnetic field Gx, the Y-axis gradient magnetic field Gy, and the Z-axis gradient magnetic field. It has the function of generating Gz and RF signals.

また、シーケンスコントローラ31は、受信器30におけるNMR信号の検波及びA/D (analog to digital)変換により得られた複素データである生データ(raw data)を受けてコンピュータ32に与えるように構成される。   The sequence controller 31 is configured to receive raw data, which is complex data obtained by detection of an NMR signal and A / D (analog to digital) conversion in the receiver 30, and supply the received raw data to the computer 32. The

このため、送信器29には、シーケンスコントローラ31から受けた制御情報に基づいてRF信号をRFコイル24に与える機能が備えられる一方、受信器30には、RFコイル24から受けたNMR信号を検波して所要の信号処理を実行するとともにA/D変換することにより、デジタル化された複素データである生データを生成する機能と生成した生データをシーケンスコントローラ31に与える機能とが備えられる。   For this reason, the transmitter 29 is provided with a function of applying an RF signal to the RF coil 24 based on the control information received from the sequence controller 31, while the receiver 30 detects the NMR signal received from the RF coil 24. Then, by executing required signal processing and A / D conversion, a function of generating raw data that is digitized complex data and a function of supplying the generated raw data to the sequence controller 31 are provided.

また、コンピュータ32の記憶装置36に保存されたプログラムを演算装置35で実行することにより、コンピュータ32には各種機能が備えられる。ただし、プログラムの少なくとも一部に代えて、各種機能を有する特定の回路を磁気共鳴イメージング装置20に設けてもよい。   Further, the computer 32 is provided with various functions by executing the program stored in the storage device 36 of the computer 32 by the arithmetic unit 35. However, a specific circuit having various functions may be provided in the magnetic resonance imaging apparatus 20 instead of at least a part of the program.

図2は、図1に示すコンピュータ32の機能ブロック図である。   FIG. 2 is a functional block diagram of the computer 32 shown in FIG.

コンピュータ32の演算装置35は、記憶装置36に保存されたプログラムを実行することにより撮像条件設定部40及びデータ処理部41として機能する。撮像条件設定部40は、傾斜磁場波形設定部40A、傾斜磁場波形シミュレーション部40B及び傾斜磁場波形取得部40Cを有する。また、記憶装置36は、k空間データ記憶部42、画像データ記憶部43及び傾斜磁場波形記憶部44として機能する。   The computing device 35 of the computer 32 functions as the imaging condition setting unit 40 and the data processing unit 41 by executing a program stored in the storage device 36. The imaging condition setting unit 40 includes a gradient magnetic field waveform setting unit 40A, a gradient magnetic field waveform simulation unit 40B, and a gradient magnetic field waveform acquisition unit 40C. The storage device 36 functions as a k-space data storage unit 42, an image data storage unit 43, and a gradient magnetic field waveform storage unit 44.

撮像条件設定部40は、入力装置33からの指示情報に基づいてパルスシーケンスを含む撮像条件を設定し、設定した撮像条件をシーケンスコントローラ31に出力する機能を有する。   The imaging condition setting unit 40 has a function of setting imaging conditions including a pulse sequence based on instruction information from the input device 33 and outputting the set imaging conditions to the sequence controller 31.

図3は、図2に示す撮像条件設定部40において設定されるEPIシーケンスの一例を示す図である。   FIG. 3 is a diagram showing an example of an EPI sequence set in the imaging condition setting unit 40 shown in FIG.

図3において横軸は時間を、RFはRF送信パルスを、GSSはスライス選択(SS: slice selection)傾斜磁場パルスを、GPEは位相エンコード(PE: phase encode)傾斜磁場パルスを、GROは読み出し(RO: readout) 傾斜磁場パルスを、ECHOはMR受信エコー信号を、それぞれ示す。 In FIG. 3, the horizontal axis represents time, RF represents an RF transmission pulse, G SS represents a slice selection (SS) gradient magnetic field pulse, G PE represents a phase encode (PE) gradient magnetic field pulse, and G RO Indicates a readout (RO: readout) gradient magnetic field pulse, and ECHO indicates an MR received echo signal.

図3に示すように、EPIシーケンスは、90°RFパルス及び180°RFパルスの印加によってSS傾斜磁場パルスの印加により選択されたスライスを励起した後、PE傾斜磁場パルス及びRO傾斜磁場パルスを繰り返し印加することによって、選択されたスライスから連続的に複数のMR受信エコー信号を収集するパルスシーケンスである。90°RFパルスの印加時刻からk空間中心におけるエコー信号の収集時刻までの時間は、実効エコー時間(effective TE: effective echo time)と呼ばれる。   As shown in FIG. 3, the EPI sequence repeats a PE gradient pulse and an RO gradient pulse after exciting a selected slice by applying an SS gradient pulse by applying a 90 ° RF pulse and a 180 ° RF pulse. This is a pulse sequence in which a plurality of MR reception echo signals are collected continuously from a selected slice by applying. The time from the time when the 90 ° RF pulse is applied to the time when the echo signal is collected at the center of the k-space is called the effective echo time (effective TE).

撮像条件設定部40において設定されたパルスシーケンスは、シーケンスコントローラ31による傾斜磁場電源27や送信器29等のイメージング系の制御に使用される。すなわち、パルスシーケンスで定まる波形を有する電流の出力指示情報が制御信号としてシーケンスコントローラ31から傾斜磁場電源27及び送信器29に出力される。   The pulse sequence set in the imaging condition setting unit 40 is used for controlling the imaging system such as the gradient magnetic field power supply 27 and the transmitter 29 by the sequence controller 31. That is, current output instruction information having a waveform determined by a pulse sequence is output from the sequence controller 31 to the gradient magnetic field power supply 27 and the transmitter 29 as a control signal.

そして、傾斜磁場コイル23及び傾斜磁場電源27を含む傾斜磁場発生システムによって撮像領域にパルスシーケンスに応じた波形を有する傾斜磁場が印加される。一方、送信器29及び送信用のRFコイル24を含むRF磁場発生システムによって撮像領域にパルスシーケンスに応じた波形を有するRF磁場パルスが印加される。そして、傾斜磁場発生システム及びRF磁場発生システムを含むイメージング系を用いて被検体PのMRイメージングを行うことができる。   A gradient magnetic field having a waveform corresponding to the pulse sequence is applied to the imaging region by a gradient magnetic field generation system including the gradient magnetic field coil 23 and the gradient magnetic field power supply 27. On the other hand, an RF magnetic field pulse having a waveform corresponding to the pulse sequence is applied to the imaging region by the RF magnetic field generation system including the transmitter 29 and the transmission RF coil 24. Then, MR imaging of the subject P can be performed using an imaging system including a gradient magnetic field generation system and an RF magnetic field generation system.

このようなパルスシーケンスは従来、傾斜磁場電源27から傾斜磁場コイル23に出力される電流の波形と、傾斜磁場コイル23から撮像領域に印加される傾斜磁場の波形とが比例関係にあるという仮定で設計されている。例えば、図3に例示されるEPIシーケンスのように、傾斜磁場電源27から出力される電流の波形を台形波とすれば、撮像領域に印加される傾斜磁場の波形も台形波になるという前提でパルスシーケンスが設定される。これは、矩形波の場合においても同様である。   Such a pulse sequence is based on the assumption that the waveform of the current output from the gradient magnetic field power supply 27 to the gradient magnetic field coil 23 and the waveform of the gradient magnetic field applied from the gradient magnetic field coil 23 to the imaging region are proportional. Designed. For example, assuming that the waveform of the current output from the gradient magnetic field power supply 27 is a trapezoidal wave as in the EPI sequence illustrated in FIG. 3, the waveform of the gradient magnetic field applied to the imaging region is also a trapezoidal wave. A pulse sequence is set. The same applies to a rectangular wave.

しかしながら、撮像領域に印加される傾斜磁場の波形は、傾斜磁場電源27から出力される電流の波形に対して厳密には相似しないことがシミュレーションや傾斜磁場の測定によって判明した。すなわち、傾斜磁場電源27から出力される電流の波形と撮像領域に印加される傾斜磁場の波形とは比例関係になく、傾斜磁場の波形には電流の波形に対する歪が生じる。   However, it has been found by simulation and measurement of the gradient magnetic field that the waveform of the gradient magnetic field applied to the imaging region is not strictly similar to the waveform of the current output from the gradient magnetic field power supply 27. That is, the waveform of the current output from the gradient magnetic field power supply 27 and the waveform of the gradient magnetic field applied to the imaging region are not proportional, and the gradient magnetic field waveform is distorted with respect to the current waveform.

この原因としては、傾斜磁場発生システムの特性が正確には設計値としての理想的な特性となっていないことが考えられる。代表的な原因としては、傾斜磁場コイル23の抵抗がチャンネル毎に異なること、傾斜磁場コイル23の応答が相互誘導によって抵抗及びインダクタンスとともに電流の周波数に依存して変化することなどが挙げられる。   This may be because the characteristics of the gradient magnetic field generation system are not exactly the ideal characteristics as design values. Typical causes include that the resistance of the gradient coil 23 varies from channel to channel, and that the response of the gradient coil 23 changes depending on the frequency of the current together with the resistance and inductance due to mutual induction.

尚、EPIシーケンスでは、傾斜磁場コイル23に出力される電流の周波数は、エコー間隔(ETS: echo train space)に相当する。従って、EPIシーケンスの場合には、ETSに依存して傾斜磁場の歪量が変化する。特に、EPIシーケンスでは、傾斜磁場の立ち上りの波形を1/4サイン波形にしようとすると、傾斜磁場の歪量が顕著となる。   In the EPI sequence, the frequency of the current output to the gradient coil 23 corresponds to the echo interval (ETS: echo train space). Therefore, in the case of the EPI sequence, the strain amount of the gradient magnetic field changes depending on the ETS. In particular, in the EPI sequence, if the rising waveform of the gradient magnetic field is set to a 1/4 sine waveform, the amount of gradient magnetic field distortion becomes significant.

上述したように傾斜磁場の波形が歪み、電流波形との誤差が大きくなると、傾斜磁場電源27は無駄な電力を消費していることになる。すなわち、傾斜磁場発生システムの性能を十分に活用していないこととなる。   As described above, when the waveform of the gradient magnetic field is distorted and an error from the current waveform becomes large, the gradient magnetic field power supply 27 consumes useless power. That is, the performance of the gradient magnetic field generation system is not fully utilized.

また、特にEPIシーケンス等の高速撮像用のパルスシーケンスでは、データ収集時間の短縮を図るために台形波であるRO傾斜磁場パルスの立ち上がり期間においてもMRデータが収集される。この場合、時間的に等間隔に収集されたMRデータをk空間に配置すると空間的に不等間隔となるため、k空間において等間隔のk空間データを得るためのregridding処理が実行される。   In particular, in a pulse sequence for high-speed imaging such as an EPI sequence, MR data is collected even during the rising period of the RO gradient magnetic field pulse that is a trapezoidal wave in order to shorten the data collection time. In this case, if MR data collected at equal intervals in time is arranged in k space, it becomes spatially unequal, and therefore regridding processing for obtaining k space data at equal intervals in k space is executed.

しかしながら、regridding処理は、傾斜磁場電源27から出力される電流の波形に基づいて実行される。すなわち、regridding処理では、傾斜磁場波形の電流波形に対する歪みが考慮されていない。このため、傾斜磁場波形の電流波形に対する歪みが無視できない場合には、regridding処理によって正確なk空間データを生成することが困難となる。加えて、傾斜磁場波形の電流波形に対する誤差が大きくなり、傾斜磁場波形の面積が小さくなると、画質の劣化やアーチファクトの発生に繋がる。   However, the regridding process is executed based on the waveform of the current output from the gradient magnetic field power supply 27. That is, in the regridding process, the distortion of the gradient magnetic field waveform with respect to the current waveform is not considered. For this reason, when the distortion of the gradient magnetic field waveform with respect to the current waveform cannot be ignored, it is difficult to generate accurate k-space data by the regridding process. In addition, when the error of the gradient magnetic field waveform with respect to the current waveform increases and the area of the gradient magnetic field waveform decreases, the image quality deteriorates and artifacts are generated.

そこで、撮像条件設定部40の傾斜磁場波形設定部40Aには、傾斜磁場電源27から出力される第1の電流の波形に応じて撮像領域に印加される第1の傾斜磁場の波形に基づいて、MRイメージング用に傾斜磁場電源27から出力される第2の電流の波形及び撮像領域に印加される第2の傾斜磁場の波形の少なくとも一方を設定又は提示する機能が備えられる。具体的には、傾斜磁場波形設定部40Aは、傾斜磁場電源27から出力される電流の波形と撮像領域に印加される傾斜磁場の波形との間における誤差を表す指標が小さくなるように第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を設定又は提示するように構成される。   Therefore, the gradient magnetic field waveform setting unit 40A of the imaging condition setting unit 40 is based on the waveform of the first gradient magnetic field applied to the imaging region in accordance with the waveform of the first current output from the gradient magnetic field power supply 27. A function for setting or presenting at least one of the waveform of the second current output from the gradient magnetic field power supply 27 for MR imaging and the waveform of the second gradient magnetic field applied to the imaging region is provided. Specifically, the gradient magnetic field waveform setting unit 40A is configured to reduce the index indicating the error between the waveform of the current output from the gradient magnetic field power supply 27 and the waveform of the gradient magnetic field applied to the imaging region. Is configured to set or present at least one of the current waveform and the second gradient magnetic field waveform.

図4は、図2に示す傾斜磁場波形設定部40Aにおいて第1の電流波形及び第1の傾斜磁場波形に基づいて生成された第2の電流波形及び第2の傾斜磁場波形の一例を示す図である。   FIG. 4 is a diagram showing an example of the second current waveform and the second gradient magnetic field waveform generated based on the first current waveform and the first gradient magnetic field waveform in the gradient magnetic field waveform setting unit 40A shown in FIG. It is.

図4の各グラフにおいて、横軸は時間を示し、縦軸は電流と傾斜磁場のスケールを合わせた相対的な振幅を示す。図4(A)は、撮像条件設定部40において初期設定された第1の電流波形及び第1の電流波形に対応する第1の傾斜磁場波形を示す。また、図4(B)は、傾斜磁場波形設定部40Aにおいて第1の電流波形及び第1の傾斜磁場波形の補正波形として生成された第2の電流波形及び第2の傾斜磁場波形を示す。   In each graph of FIG. 4, the horizontal axis represents time, and the vertical axis represents the relative amplitude of the current and gradient magnetic field scales. FIG. 4A shows a first current waveform initially set in the imaging condition setting unit 40 and a first gradient magnetic field waveform corresponding to the first current waveform. FIG. 4B shows the second current waveform and the second gradient magnetic field waveform generated as the correction waveform of the first current waveform and the first gradient magnetic field waveform in the gradient magnetic field waveform setting unit 40A.

図4(A)に示すように、パルスシーケンスにおいてRO傾斜磁場パルスを台形波として定義して傾斜磁場電源27から第1の電流を出力しても、撮像領域に印加される第1の傾斜磁場の波形には、第1の電流波形に対する歪みが生じる。これは、矩形波等の他の波形についても同様であり、矩形波として電流を出力しても、傾斜磁場は正確な矩形波とはならない。   As shown in FIG. 4A, even if the RO gradient magnetic field pulse is defined as a trapezoidal wave in the pulse sequence and the first current is output from the gradient magnetic field power supply 27, the first gradient magnetic field applied to the imaging region. In this waveform, distortion occurs with respect to the first current waveform. The same applies to other waveforms such as a rectangular wave. Even if a current is output as a rectangular wave, the gradient magnetic field does not become an accurate rectangular wave.

そこで、第1の傾斜磁場波形の第1の電流波形に対する誤差を低減させる補正処理を行うことによって図4(B)に示すように、乖離量が低減された第2の電流波形及び第2の傾斜磁場波形を得ることができる。具体的には、振幅が不連続に変化する部分が滑らかに変化するように補正された第2の電流波形が得られる。従って、初期波形である第1の電流波形が台形波であれば、立ち上がりの部分を鈍らせた波形が第2の電流波形となる。   Therefore, as shown in FIG. 4B, by performing a correction process for reducing an error of the first gradient magnetic field waveform with respect to the first current waveform, as shown in FIG. A gradient magnetic field waveform can be obtained. Specifically, a second current waveform corrected to smoothly change a portion where the amplitude changes discontinuously is obtained. Therefore, if the first current waveform that is the initial waveform is a trapezoidal wave, the waveform in which the rising portion is blunted becomes the second current waveform.

尚、第2の傾斜磁場波形は、第2の電流波形に応じた波形となるため、第2の傾斜磁場波形は、第1の傾斜磁場波形とは異なる波形となる。   Since the second gradient magnetic field waveform is a waveform corresponding to the second current waveform, the second gradient magnetic field waveform is different from the first gradient magnetic field waveform.

傾斜磁場波形設定部40Aでは、第2の電流波形及び第2の傾斜磁場波形のいずれかを求めるようにしてもよいし、第2の電流波形及び第2の傾斜磁場波形の組合せを求めるようにしてもよい。そして、求められた第2の電流波形及び第2の傾斜磁場波形は、撮像条件として自動設定することができる。或いは、第2の電流波形及び第2の傾斜磁場波形を撮像条件の候補として表示装置34に一旦提示し、入力装置33の操作によるユーザの確認を促すようにしてもよい。   In the gradient magnetic field waveform setting unit 40A, either the second current waveform or the second gradient magnetic field waveform may be obtained, or a combination of the second current waveform and the second gradient magnetic field waveform may be obtained. May be. The obtained second current waveform and second gradient magnetic field waveform can be automatically set as imaging conditions. Alternatively, the second current waveform and the second gradient magnetic field waveform may be once presented on the display device 34 as candidates for imaging conditions, and the user may be prompted to confirm the operation by operating the input device 33.

図4に示す第1及び第2の傾斜磁場波形は、それぞれ第1及び第2の電流波形に基づくシミュレーションによって計算された波形である。図4に例示されるように、電流波形に対応する傾斜磁場の波形は、傾斜磁場コイル23及び傾斜磁場電源27を含む傾斜磁場発生システムを適切にモデル化した傾斜磁場波形シミュレーションによって求めることができる。従って、第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を設定又は提示するための元データとなる第1の傾斜磁場の波形として、第1の電流の波形に基づくシミュレーションによって取得された傾斜磁場の波形を用いることができる。   The first and second gradient magnetic field waveforms shown in FIG. 4 are waveforms calculated by simulation based on the first and second current waveforms, respectively. As illustrated in FIG. 4, the gradient magnetic field waveform corresponding to the current waveform can be obtained by a gradient magnetic field waveform simulation that appropriately models a gradient magnetic field generation system including the gradient coil 23 and the gradient magnetic field power supply 27. . Therefore, the waveform of the first gradient magnetic field that is the original data for setting or presenting at least one of the waveform of the second current and the waveform of the second gradient magnetic field is acquired by simulation based on the waveform of the first current. A gradient magnetic field waveform can be used.

但し、シミュレーションに限らず、他の方法によって傾斜磁場電源27から出力される電流の波形に対応する傾斜磁場の波形を求めることもできる。具体例として、パルスシーケンスをシーケンスコントローラ31に出力することによって傾斜磁場電源27から実際に対象となる波形の電流を出力させ、撮像領域に印加された傾斜磁場を測定する方法が挙げられる。また、別の具体例として、パルスシーケンスをシーケンスコントローラ31に出力することによって傾斜磁場電源27から実際に対象となる波形の電流を出力させ、傾斜磁場が印加された状態で収集されたMR信号に基づいて傾斜磁場を計算する方法が挙げられる。   However, the waveform of the gradient magnetic field corresponding to the waveform of the current output from the gradient magnetic field power supply 27 can be obtained by other methods without being limited to the simulation. As a specific example, there is a method of outputting a pulse sequence to the sequence controller 31 to output a current having a waveform actually targeted from the gradient magnetic field power supply 27 and measuring a gradient magnetic field applied to the imaging region. As another specific example, by outputting a pulse sequence to the sequence controller 31, the current of the waveform actually targeted is output from the gradient magnetic field power supply 27, and the MR signal collected in a state where the gradient magnetic field is applied is output. A method of calculating a gradient magnetic field based on the above is mentioned.

つまり、傾斜磁場発生システムにより実際に傾斜磁場電源27から第1の電流の波形を有する電流を出力させることによって第1の傾斜磁場の波形を有する傾斜磁場を撮像領域に印加し、撮像領域に印加された第1の傾斜磁場の波形に基づいて第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を設定又は提示することができる。   That is, the gradient magnetic field generating system applies a gradient magnetic field having the first gradient magnetic field waveform to the imaging region by actually outputting a current having the first current waveform from the gradient magnetic field power supply 27, and applies it to the imaging region. Based on the waveform of the first gradient magnetic field, at least one of the waveform of the second current and the waveform of the second gradient magnetic field can be set or presented.

従って、第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を求めるために必要となる第1の傾斜磁場の波形としては、傾斜磁場波形シミュレーションによって取得された傾斜磁場の波形、第1の電流の波形を傾斜磁場電源27から実際に出力させることによって測定された傾斜磁場の波形或いは第1の電流の波形を傾斜磁場電源27から実際に出力させることによって収集されたMR信号に基づいて求められた傾斜磁場の波形など、様々な方法によって取得した波形を用いることができる。   Therefore, the waveform of the first gradient magnetic field required for obtaining at least one of the waveform of the second current and the waveform of the second gradient magnetic field includes the waveform of the gradient magnetic field acquired by the gradient magnetic field waveform simulation, Based on the MR signal collected by actually outputting the waveform of the gradient magnetic field measured by actually outputting the waveform of one current from the gradient magnetic field power supply 27 or the waveform of the first current from the gradient magnetic field power supply 27. Waveforms acquired by various methods such as the waveform of the gradient magnetic field obtained in this way can be used.

そこで、傾斜磁場波形シミュレーション部40Bには、傾斜磁場波形シミュレーションによって傾斜磁場電源27から出力される電流の波形に対応する傾斜磁場の波形を求める機能が備えられる。また、傾斜磁場波形取得部40Cには、傾斜磁場波形シミュレーション以外の任意の方法によって傾斜磁場電源27から出力される電流の波形に対応する傾斜磁場の波形を求める機能が備えられる。そして、傾斜磁場波形シミュレーション部40B又は傾斜磁場波形取得部40Cにおいて求められた傾斜磁場の波形は、傾斜磁場波形設定部40Aに与えられるように構成される。   Therefore, the gradient magnetic field waveform simulation unit 40B has a function of obtaining a gradient magnetic field waveform corresponding to the waveform of the current output from the gradient magnetic field power supply 27 by the gradient magnetic field waveform simulation. The gradient magnetic field waveform acquisition unit 40C has a function of obtaining a gradient magnetic field waveform corresponding to the waveform of the current output from the gradient magnetic field power supply 27 by an arbitrary method other than the gradient magnetic field waveform simulation. The gradient magnetic field waveform obtained by the gradient magnetic field waveform simulation unit 40B or the gradient magnetic field waveform acquisition unit 40C is configured to be provided to the gradient magnetic field waveform setting unit 40A.

傾斜磁場波形の測定は、サーチコイル及び積分器を含む傾斜磁場測定系を用いて行うことができる。その場合には、図2に示すように、サーチコイル50、積分器51、増幅器52及びA/D変換器53が磁気共鳴イメージング装置20に備えられる。   The gradient magnetic field waveform can be measured using a gradient magnetic field measurement system including a search coil and an integrator. In that case, as shown in FIG. 2, the magnetic resonance imaging apparatus 20 includes a search coil 50, an integrator 51, an amplifier 52, and an A / D converter 53.

サーチコイル50は、傾斜磁場が印加される撮像領域に設置される。そして、サーチコイル50によって傾斜磁場の微分成分が検出される。サーチコイル50によって検出された傾斜磁場の微分成分は積分器51に出力される。そして、積分器51における傾斜磁場の微分成分の積分処理によって傾斜磁場の波形が再現される。   The search coil 50 is installed in an imaging region to which a gradient magnetic field is applied. Then, the differential component of the gradient magnetic field is detected by the search coil 50. The differential component of the gradient magnetic field detected by the search coil 50 is output to the integrator 51. Then, the waveform of the gradient magnetic field is reproduced by integration processing of the differential component of the gradient magnetic field in the integrator 51.

積分器51において取得された傾斜磁場の波形は、増幅器52で増幅された後、A/D変換器53においてデジタル情報に変換される。そして、デジタル信号化された傾斜磁場の波形がA/D変換器53からコンピュータ32の傾斜磁場波形取得部40Cに出力されるように構成される。   The waveform of the gradient magnetic field acquired by the integrator 51 is amplified by the amplifier 52 and then converted into digital information by the A / D converter 53. The digital magnetic field gradient magnetic field waveform is output from the A / D converter 53 to the gradient magnetic field waveform acquisition unit 40C of the computer 32.

一方、傾斜磁場の波形は、傾斜磁場の波形の取得用のパルスシーケンスに従って収集されたMR信号に基づいて計算することもできる。そこで、撮像条件設定部40には、傾斜磁場の波形の取得用のパルスシーケンスをデータ収集条件として設定し、設定したパルスシーケンスを含むデータ収集条件をシーケンスコントローラ31に出力する機能が備えられる。   On the other hand, the gradient magnetic field waveform can also be calculated based on the MR signal collected according to the pulse sequence for acquiring the gradient magnetic field waveform. Therefore, the imaging condition setting unit 40 has a function of setting a pulse sequence for acquiring a gradient magnetic field waveform as a data acquisition condition and outputting the data acquisition condition including the set pulse sequence to the sequence controller 31.

また、傾斜磁場波形取得部40Cには、傾斜磁場の波形の取得用のパルスシーケンスに従って収集されたMR信号をシーケンスコントローラ31からデータ処理部41を通じて取得し、取得したMR信号に基づいて傾斜磁場の波形を求める機能が備えられる。   The gradient magnetic field waveform acquisition unit 40C acquires MR signals collected according to the pulse sequence for acquiring the gradient magnetic field waveform from the sequence controller 31 through the data processing unit 41, and based on the acquired MR signals, A function for obtaining a waveform is provided.

図5は、図2に示す撮像条件設定部40において設定される傾斜磁場の波形の取得用のパルスシーケンスの一例を示すシーケンスチャートである。   FIG. 5 is a sequence chart showing an example of a pulse sequence for acquiring the gradient magnetic field waveform set in the imaging condition setting unit 40 shown in FIG.

図5において横軸は時間を、RFはRF送信パルス及びデータ収集期間を、Gは傾斜磁場パルスを、それぞれ示す。図5に示すように、RF励起パルスとともにSS傾斜磁場パルスGssを印加した後に、エンコード傾斜磁場パルスGek及び計測対象となる傾斜磁場Gmeasuredを順に印加するパルスシーケンスを、エンコード傾斜磁場パルスGekの強度をステップ状に変化させながら繰返すことによって傾斜磁場波形の取得用のMR信号を収集することができる。   In FIG. 5, the horizontal axis represents time, RF represents an RF transmission pulse and a data acquisition period, and G represents a gradient magnetic field pulse. As shown in FIG. 5, after applying the SS gradient magnetic field pulse Gss together with the RF excitation pulse, a pulse sequence in which the encode gradient magnetic field pulse Gek and the gradient magnetic field Gmeasured to be measured are sequentially applied, and the intensity of the encode gradient magnetic field pulse Gek is set. The MR signal for acquiring the gradient magnetic field waveform can be collected by repeating the process while changing it stepwise.

MR信号の収集期間は計測対象となる傾斜磁場Gmeasuredの印加期間中とされ、時刻t (t=t1, t2, t3, ..., tm)における時系列の自由誘導減衰(FID: free induction decay)信号が収集される。従って、エンコード傾斜磁場パルスGekの強度をGek=Ge1, Ge2, Ge3, ..., Genまでステップ状に変化させると、m×n個のMRデータS(t, Gek)が収集される。   The MR signal acquisition period is the period during which the gradient magnetic field Gmeasured to be measured is applied, and time-series free induction decay (FID) at time t (t = t1, t2, t3, ..., tm) ) A signal is collected. Therefore, when the intensity of the encode gradient magnetic field pulse Gek is changed stepwise from Gek = Ge1, Ge2, Ge3,..., Gen, m × n MR data S (t, Gek) are collected.

そして、傾斜磁場波形取得部40Cでは、MRデータS(t, Gek)に基づいて計測対象となる傾斜磁場Gmeasuredを求めることができる。具体的には、式(1)の演算によって時間的に隣接するMRデータS(t, Gek)間の位相差を位相角とするk空間データD(t, Gek)が求められる。
D(t, Gek)=S(t, Gek)*S*(t+Δt, Gek) (1)
但し、式(1)においてデータS*(t, Gek)は、データS(t, Gek)の共役複素数であり、Δtは隣接する時刻間における時間差である。
Then, the gradient magnetic field waveform acquisition unit 40C can obtain the gradient magnetic field Gmeasured to be measured based on the MR data S (t, Gek). Specifically, k-space data D (t, Gek) whose phase angle is a phase difference between MR data S (t, Gek) that are temporally adjacent is obtained by the calculation of equation (1).
D (t, Gek) = S (t, Gek) * S * (t + Δt, Gek) (1)
However, in equation (1), data S * (t, Gek) is a conjugate complex number of data S (t, Gek), and Δt is a time difference between adjacent times.

次に、エンコード傾斜磁場パルスGekの強度が同一のk空間データD(t, Gek)を加算することによって時刻tごとのk空間データd(t)が求められる。そうすると、式(2)により隣接する時刻間における傾斜磁場差ΔG(t)を求めることができる。
ΔG(t)=arctan{d(t)}/(2πγzΔt) (2)
但し、式(2)においてarctan()は逆正接関数、γは磁気回転比、zは傾斜磁場Gmeasuredの印加軸方向におけるスライス位置、Δtは隣接する時刻間における時間差である。
Next, k-space data d (t) for each time t is obtained by adding the k-space data D (t, Gek) having the same intensity of the encode gradient magnetic field pulse Gek. Then, the gradient magnetic field difference ΔG (t) between adjacent times can be obtained from the equation (2).
ΔG (t) = arctan {d (t)} / (2πγzΔt) (2)
In Equation (2), arctan () is an arctangent function, γ is a magnetic rotation ratio, z is a slice position in the direction of the application axis of the gradient magnetic field Gmeasured, and Δt is a time difference between adjacent times.

そして、式(2)により得られた傾斜磁場差ΔG(t)を任意の時刻Tまで積算することによって時刻Tにおける傾斜磁場の強度G(T)、つまり測定対象となる傾斜磁場Gmeasuredの波形を求めることができる。   Then, by integrating the gradient magnetic field difference ΔG (t) obtained by the equation (2) until an arbitrary time T, the gradient G of the gradient magnetic field at time T, that is, the waveform of the gradient magnetic field Gmeasured to be measured is obtained. Can be sought.

尚、図5に示されたパルスシーケンスは一例であり、他のパルスシーケンスを用いて傾斜磁場の波形を求めることもできる。例えば、図5に示すエンコード傾斜磁場パルスGekに代えてリフェーズパルスの印加を伴う時系列のFID信号S(t)の収集を1回実行し、FID信号S(t)に基づいて傾斜磁場の波形を求めることもできる。この場合には、エンコード傾斜磁場パルスGekの強度ごとのk空間データの加算処理を省略した計算によって同様に傾斜磁場Gmeasuredの波形を求めることができる。   Note that the pulse sequence shown in FIG. 5 is an example, and the waveform of the gradient magnetic field can be obtained using another pulse sequence. For example, instead of the encode gradient magnetic field pulse Gek shown in FIG. 5, the collection of the time-series FID signal S (t) accompanied by the application of a rephase pulse is executed once, and the gradient magnetic field is generated based on the FID signal S (t). A waveform can also be obtained. In this case, the waveform of the gradient magnetic field Gmeasured can be similarly obtained by a calculation in which the k-space data addition process for each intensity of the encode gradient magnetic field pulse Gek is omitted.

一方、上述したように傾斜磁場波形シミュレーション部40Bにおける傾斜磁場波形シミュレーションによれば、実際に傾斜磁場電源27から電流を出力させずに電流波形に対応する傾斜磁場波形を計算によって求めることができる。傾斜磁場波形シミュレーションは、傾斜磁場電源27、傾斜磁場コイル26及びシーケンスコントローラ31のように傾斜磁場の生成に寄与する構成要素を可能な限り含む傾斜磁場発生システムの等価回路モデルを用いて実行することができる。   On the other hand, as described above, according to the gradient magnetic field waveform simulation in the gradient magnetic field waveform simulation unit 40B, the gradient magnetic field waveform corresponding to the current waveform can be obtained by calculation without actually outputting the current from the gradient magnetic field power supply 27. The gradient magnetic field waveform simulation is executed by using an equivalent circuit model of a gradient magnetic field generation system including as much as possible the components contributing to the generation of the gradient magnetic field, such as the gradient magnetic field power supply 27, the gradient magnetic field coil 26, and the sequence controller 31. Can do.

図6は、図1に示す傾斜磁場電源27及び傾斜磁場コイル23を含む傾斜磁場発生システムの詳細構成例を示す回路図である。   FIG. 6 is a circuit diagram illustrating a detailed configuration example of the gradient magnetic field generation system including the gradient magnetic field power supply 27 and the gradient magnetic field coil 23 illustrated in FIG. 1.

図6に示すように、傾斜磁場電源27は、ブレーカ122と、整流器123と、直流電源124と、電解コンデンサ126、126'、126"と、傾斜磁場アンプ128、128'、128"と、電流検出器130、130'、130"とを有する。即ち、図1に示すX、Y、Z軸傾斜磁場電源27x、27y、27zは、ブレーカ122と、整流器123と、定電圧(CV: Constant Voltage)/定電流(CC: Constant Current)特性を持つ直流電源124とを共有する。   As shown in FIG. 6, the gradient magnetic field power supply 27 includes a breaker 122, a rectifier 123, a DC power supply 124, electrolytic capacitors 126, 126 ′, 126 ″, gradient magnetic field amplifiers 128, 128 ′, 128 ″, and a current. In other words, the X, Y, and Z axis gradient magnetic field power supplies 27x, 27y, and 27z shown in FIG. 1 include a breaker 122, a rectifier 123, and a constant voltage (CV). ) / DC power supply 124 having a constant current (CC) characteristic.

尚、CV/CC特性を持つ直流電源は、負荷が軽いときには定電圧を出力して負荷が重くなり、ある電流以上を流す必要が生じた際には、一定電流を負荷に供給する制御を行う電源である。   Note that a DC power supply with CV / CC characteristics outputs a constant voltage when the load is light and the load becomes heavy, and when it is necessary to flow more than a certain current, control is performed to supply a constant current to the load. Power supply.

ブレーカ122は、外部の交流電源120からの出力電流が定格電流値を超えた場合に、交流電源120と、整流器123との間を電気的に遮断する。整流器123は、交流電源120からの交流の供給電力を直流電力に変換して直流電源124に供給する。   Breaker 122 electrically disconnects between AC power supply 120 and rectifier 123 when the output current from external AC power supply 120 exceeds the rated current value. The rectifier 123 converts AC supply power from the AC power source 120 into DC power and supplies it to the DC power source 124.

直流電源124は、整流器123を介して供給される直流電流で電解コンデンサ126、126'、126"を充電し、傾斜磁場アンプ128、128'、128"に直流電流を供給する。また、直流電源124は、傾斜磁場アンプ128、128'、128"側の負荷が軽い場合には定電圧源として動作し、負荷が重い場合には定電流源として動作する。   The DC power supply 124 charges the electrolytic capacitors 126, 126 ′, 126 ″ with a DC current supplied via the rectifier 123, and supplies the DC current to the gradient magnetic field amplifiers 128, 128 ′, 128 ″. The DC power source 124 operates as a constant voltage source when the load on the gradient magnetic field amplifiers 128, 128 ′, and 128 ″ is light, and operates as a constant current source when the load is heavy.

各傾斜磁場アンプ128、128'、128"は、それぞれ+側入力端子(+in)、−側入力端子(-in)、+側出力端子(+out)及び−側出力端子(-out)を有する。また、傾斜磁場アンプ128、128'、128"は、直流電源124から電力供給をそれぞれ受けると共に、シーケンスコントローラ31からの制御信号を+側入力端子においてそれぞれ受ける。シーケンスコントローラ31から傾斜磁場アンプ128、128'128"に入力される各制御信号は、パルスシーケンスに従ってX、Y、Z軸傾斜磁場コイル23x、23y、23zによりそれぞれ撮像領域に印加すべき傾斜磁場波形に相似する波形を有する電圧信号である。   Each gradient magnetic field amplifier 128, 128 ', 128 "has a + side input terminal (+ in), a-side input terminal (-in), a + side output terminal (+ out), and a-side output terminal (-out), respectively. The gradient magnetic field amplifiers 128, 128 ′, and 128 ″ receive power supply from the DC power supply 124 and receive control signals from the sequence controller 31 at the + side input terminals. The control signals input from the sequence controller 31 to the gradient magnetic field amplifiers 128 and 128'128 "are gradient magnetic field waveforms to be applied to the imaging region by the X, Y, and Z axis gradient magnetic field coils 23x, 23y, and 23z, respectively, according to the pulse sequence. Is a voltage signal having a waveform similar to.

電流検出器130、130'、130"はそれぞれ、傾斜磁場アンプ128、128'、128"の−側出力端子に流入する電流の電流値を検出する。電流検出器130、130'、130"において検出される電流の大きさは、傾斜磁場アンプ128、128'、128"の+側出力端子から出力される電流に等しい。これは、各傾斜磁場アンプ128、128'、128"の+側出力端子から出力される電流はそれぞれ、傾斜磁場コイル23x、23y、23zを流れて傾斜磁場アンプ128、128'、128"の−側出力端子に戻るからである。   The current detectors 130, 130 ′, and 130 ″ detect current values of currents flowing into the negative side output terminals of the gradient magnetic field amplifiers 128, 128 ′, and 128 ″, respectively. The magnitude of the current detected by the current detectors 130, 130 ′, 130 ″ is equal to the current output from the + side output terminal of the gradient magnetic field amplifier 128, 128 ′, 128 ″. This is because currents output from the + side output terminals of the gradient magnetic field amplifiers 128, 128 ′, and 128 ″ flow through the gradient magnetic field coils 23x, 23y, and 23z, respectively, and − of the gradient magnetic field amplifiers 128, 128 ′, and 128 ″. This is because it returns to the side output terminal.

電流検出器130、130'、130"はそれぞれ、検出した電流の電流値を示す電圧信号を生成して、生成した電圧信号を傾斜磁場アンプ128、128'、128"の−側入力端子に入力する。   Each of the current detectors 130, 130 ′, and 130 ″ generates a voltage signal indicating the current value of the detected current, and inputs the generated voltage signal to the − side input terminal of the gradient magnetic field amplifier 128, 128 ′, and 128 ″. To do.

傾斜磁場アンプ128、128'、128"は、+/−側入力端子間の誤差信号が0になるように電流を出力する電流源として動作する。ここで、傾斜磁場アンプ128、128'、128"の出力電流は電流検出器130、130'、130"によって負帰還となっている。このため、+側入力端子への入力電圧に比例する電流を+側出力端子から出力するように、フィードバック制御が行われる。   The gradient magnetic field amplifiers 128, 128 ′, and 128 ″ operate as current sources that output current so that the error signal between the +/− side input terminals becomes 0. Here, the gradient magnetic field amplifiers 128, 128 ′, and 128 The output current of "is negatively fed back by the current detectors 130, 130 ', 130". Therefore, feedback is performed so that a current proportional to the input voltage to the + side input terminal is output from the + side output terminal. Control is performed.

図7は、図6に示す傾斜磁場発生システムの等価回路モデルの一例を示す回路図である。   FIG. 7 is a circuit diagram showing an example of an equivalent circuit model of the gradient magnetic field generation system shown in FIG.

図7において、X軸傾斜磁場電源27xは、図6に示すブレーカ122、整流器123、直流電源124、電解コンデンサ126、傾斜磁場アンプ128及び電流検出器130とに対応する。同様に、Y軸傾斜磁場電源27yは、ブレーカ122、整流器123、直流電源124、電解コンデンサ126'、傾斜磁場アンプ128'及び電流検出器130'に対応し、Z軸傾斜磁場電源27zは、ブレーカ122、整流器123、直流電源124、電解コンデンサ126"、傾斜磁場アンプ128"及び電流検出器130"に対応する。そこで、ここでは、X軸傾斜磁場の傾斜磁場発生システムの等価回路モデルについて説明する。   In FIG. 7, the X-axis gradient magnetic field power supply 27x corresponds to the breaker 122, rectifier 123, DC power supply 124, electrolytic capacitor 126, gradient magnetic field amplifier 128, and current detector 130 shown in FIG. Similarly, the Y-axis gradient magnetic field power supply 27y corresponds to the breaker 122, the rectifier 123, the DC power supply 124, the electrolytic capacitor 126 ′, the gradient magnetic field amplifier 128 ′, and the current detector 130 ′, and the Z-axis gradient magnetic field power supply 27z corresponds to the breaker. 122, a rectifier 123, a DC power supply 124, an electrolytic capacitor 126 ", a gradient magnetic field amplifier 128", and a current detector 130 ". Therefore, here, an equivalent circuit model of an X axis gradient magnetic field gradient magnetic field generation system will be described. .

図7に示すように等価回路モデル140xは、1次側として、X軸傾斜磁場電源27x、X軸傾斜磁場コイル23xの抵抗成分に相当する抵抗23xR及びX軸傾斜磁場コイル23xのインダクタンス成分に相当するコイル23xLを直列に接続した構成である。   As shown in FIG. 7, the equivalent circuit model 140x is equivalent to the X-axis gradient magnetic field power source 27x, the resistance component 23xR corresponding to the resistance component of the X-axis gradient magnetic field coil 23x, and the inductance component of the X-axis gradient magnetic field coil 23x as the primary side. The coil 23xL to be connected is connected in series.

また、等価回路モデル140xは、抵抗141R及び仮想コイル141Lの直列回路を第1の2次側回路として有する。さらに、等価回路モデル140xは、抵抗142R及び仮想コイル142Lの直列回路を第2の2次側回路として有する。コイル23xLと、仮想コイル141Lとが互いに磁気的に結合している。また、コイル23xLと、仮想コイル142Lとが互いに磁気的に結合している。   The equivalent circuit model 140x has a series circuit of a resistor 141R and a virtual coil 141L as a first secondary circuit. Furthermore, the equivalent circuit model 140x has a series circuit of a resistor 142R and a virtual coil 142L as a second secondary circuit. The coil 23xL and the virtual coil 141L are magnetically coupled to each other. The coil 23xL and the virtual coil 142L are magnetically coupled to each other.

周波数が高くなると、X、Y、Z軸傾斜磁場コイル23x、23y、23zの各インピーダンスは、1つの抵抗成分及び1つのインダクタンス成分の和で表される単純なモデルのように単純に増加しない。例えば、高周波電流が導体を流れる時、実際には、電流密度が導体の表面で高く、表面から離れると低くなる。すなわち、表皮効果により、周波数が高くなるほど電流が表面へ集中する。この結果、導体の交流抵抗は高くなる。   As the frequency increases, the impedances of the X, Y, and Z axis gradient magnetic field coils 23x, 23y, and 23z do not simply increase as in a simple model represented by the sum of one resistance component and one inductance component. For example, when high-frequency current flows through a conductor, the current density is actually high on the surface of the conductor and decreases as it moves away from the surface. That is, due to the skin effect, the current concentrates on the surface as the frequency increases. As a result, the AC resistance of the conductor is increased.

従って、表皮効果等を考慮すると、傾斜磁場発生システムのインピーダンスを表す多項式において、X軸傾斜磁場コイル23xの抵抗23xRの抵抗値が含まれる項も、周波数に依存して変化する項とすることが望ましい。そこで、例えば、抵抗23xRの抵抗値にも角周波数ωが乗じられるような等価回路モデル、つまりインピーダンスの虚部のみならず実部も周波数依存性を有する等価回路モデルを傾斜磁場波形シミュレーションに用いることが望ましい。   Therefore, in consideration of the skin effect and the like, in the polynomial representing the impedance of the gradient magnetic field generating system, the term including the resistance value of the resistor 23xR of the X-axis gradient magnetic field coil 23x may be a term that varies depending on the frequency. desirable. Therefore, for example, an equivalent circuit model in which the resistance value of the resistor 23xR is multiplied by the angular frequency ω, that is, an equivalent circuit model in which not only the imaginary part of impedance but also the real part has frequency dependence is used for the gradient magnetic field waveform simulation. Is desirable.

また、実際には、X、Y、Z軸傾斜磁場コイル23x、23y、23zにパルス電流を供給すると渦電流が発生する。このため、渦電流による渦磁場がX、Y、Z軸方向の各傾斜磁場Gx、Gy、Gzに加わって傾斜磁場分布の歪みが生じる。渦磁場を考慮すると、相互インダクタンスも含まれた等価回路モデルを傾斜磁場波形シミュレーションに用いることが望ましい。   In practice, an eddy current is generated when a pulse current is supplied to the X, Y, and Z axis gradient magnetic field coils 23x, 23y, and 23z. For this reason, the eddy magnetic field due to the eddy current is added to the gradient magnetic fields Gx, Gy, and Gz in the X, Y, and Z axis directions, and the gradient magnetic field distribution is distorted. Considering the eddy magnetic field, it is desirable to use an equivalent circuit model including the mutual inductance for the gradient magnetic field waveform simulation.

また、傾斜磁場発生システムには、所定の周波数以上の高周波電流を遮断するためのチョークコイルが設けられる場合がある。この場合、複数の相互インダクタンスが含まれた等価回路モデルを傾斜磁場波形シミュレーションに用いることが望ましい。   Moreover, the gradient magnetic field generation system may be provided with a choke coil for cutting off a high-frequency current having a predetermined frequency or higher. In this case, it is desirable to use an equivalent circuit model including a plurality of mutual inductances for the gradient magnetic field waveform simulation.

図7に示す等価回路モデル140xは、このような点を考慮して作成されたモデルの一例である。従って、図7に例示される等価回路モデル140xに限らず、傾斜磁場発生システムの回路構成や特性に応じた等価回路モデルを作成することができる。   An equivalent circuit model 140x shown in FIG. 7 is an example of a model created in consideration of such points. Therefore, not only the equivalent circuit model 140x illustrated in FIG. 7 but also an equivalent circuit model according to the circuit configuration and characteristics of the gradient magnetic field generation system can be created.

図7に示す等価回路モデル140xの場合には、X軸方向における傾斜磁場の波形Gx(t)を式(3)によって計算することができる。
Gx(t)= aIout(t)+bI1(t)+cI2(t) (3)
但し、式(3) において、aは図7に示すコイル23xLの電流感度、bは仮想コイル141Lの電流感度、cは仮想コイル142Lの電流感度、Iout(t)はX軸傾斜磁場電源27xからの出力電流波形、I1(t)は第1の2次側回路を流れる電流、I2(t)は第2の2次側回路を流れる電流である。尚、電流感度は、コイルに電流を流すことで発生する傾斜磁場の強度[T/m]を、コイルに流す電流値[A]で除算して得られる定数である。
In the case of the equivalent circuit model 140x shown in FIG. 7, the gradient magnetic field waveform Gx (t) in the X-axis direction can be calculated by Equation (3).
Gx (t) = aIout (t) + bI 1 (t) + cI 2 (t) (3)
In equation (3), a is the current sensitivity of the coil 23xL shown in FIG. 7, b is the current sensitivity of the virtual coil 141L, c is the current sensitivity of the virtual coil 142L, and Iout (t) is from the X-axis gradient magnetic field power supply 27x. , I 1 (t) is a current flowing through the first secondary circuit, and I 2 (t) is a current flowing through the second secondary circuit. The current sensitivity is a constant obtained by dividing the strength [T / m] of the gradient magnetic field generated by flowing a current through the coil by the current value [A] flowing through the coil.

式(3)によって傾斜磁場の波形Gx(t)を計算すれば、X軸傾斜磁場電源27xからの出力電流Iout(t)に比例しない渦磁場等の磁場が加算された正確な傾斜磁場の波形を求めることができる。すなわち、仮想コイル141L、142Lによって生じる仮想的な傾斜磁場波形bI1(t)、cI2(t)をX軸傾斜磁場電源27xからの出力電流Iout(t)に比例する傾斜磁場波形aIout(t)に加算した正確な傾斜磁場の波形Gx(t)を求めることができる。 If the gradient magnetic field waveform Gx (t) is calculated by Equation (3), an accurate gradient magnetic field waveform in which a magnetic field such as an eddy magnetic field that is not proportional to the output current Iout (t) from the X-axis gradient magnetic field power supply 27x is added. Can be requested. That is, the gradient magnetic field waveforms aIout (t) in which the virtual gradient magnetic field waveforms bI 1 (t) and cI 2 (t) generated by the virtual coils 141L and 142L are proportional to the output current Iout (t) from the X-axis gradient magnetic field power supply 27x. It is possible to obtain an accurate gradient magnetic field waveform Gx (t) added to ().

式(3)においてX軸傾斜磁場電源27xからの出力電流Iout(t)は、パルスシーケンスで定められるため既知である。一方、仮想コイル141Lを流れる電流I1(t)及び仮想コイル142Lを流れる電流I2(t)は、それぞれ式(4-1)及び式(4-2)により計算することができる。
R1I1(t)+L1{dI1(t)/dt}+M1{dIout(t)/dt}=0 (4-1)
R2I2(t)+L2{dI2(t)/dt}+M2{dIout(t)/dt}=0 (4-2)
但し、式(4-1)及び式(4-2)において、R1は抵抗141Rの抵抗値、R2は抵抗142Rの抵抗値、L1は仮想コイル141Lの自己インダクタンス値、L2は仮想コイル142Lの自己インダクタンス値、M1はコイル26xLと仮想コイル141Lとの相互インダクタンス値、M2はコイル26xLと仮想コイル142Lとの相互インダクタンス値である。
In Expression (3), the output current Iout (t) from the X-axis gradient magnetic field power supply 27x is known because it is determined by a pulse sequence. On the other hand, the current I 1 (t) flowing through the virtual coil 141L and the current I 2 (t) flowing through the virtual coil 142L can be calculated by equations (4-1) and (4-2), respectively.
R 1 I 1 (t) + L 1 {dI 1 (t) / dt} + M 1 {dIout (t) / dt} = 0 (4-1)
R 2 I 2 (t) + L 2 {dI 2 (t) / dt} + M 2 {dIout (t) / dt} = 0 (4-2)
However, in Equation (4-1) and Equation (4-2), R 1 is the resistance value of the resistor 141R, R 2 is the resistance value of the resistor 142R, L 1 is the self-inductance value of the virtual coil 141L, and L 2 is the virtual value. The self-inductance value of the coil 142L, M 1 is the mutual inductance value of the coil 26xL and the virtual coil 141L, and M 2 is the mutual inductance value of the coil 26xL and the virtual coil 142L.

更に、等価回路モデル140xにおいて、X軸傾斜磁場電源44xから見たX軸傾斜磁場コイル26xのインピーダンスは、X軸傾斜磁場コイル26xにおける抵抗23xRの抵抗値Rload及びコイル23xLの自己インダクタンス値Lload並びに式(4-1)及び式(4-2)の係数となっている抵抗141R及び抵抗142Rの各抵抗値R1、R2、仮想コイル141L及び仮想コイル142Lの自己インダクタンス値L1、L2、コイル26xLと仮想コイル141L及び仮想コイル142Lとの相互インダクタンス値M1、M2を含む式で表すことができる。具体的には、実部及び虚部についての2つの式でX軸傾斜磁場コイル26xのインピーダンスを表すことができる。 Further, in the equivalent circuit model 140x, the impedance of the X-axis gradient magnetic field coil 26x as viewed from the X-axis gradient magnetic field power supply 44x includes the resistance value Rload of the resistor 23xR and the self-inductance value Lload of the coil 23xL in the X-axis gradient magnetic field coil 26x. (4-1) and the resistance values R 1 and R 2 of the resistor 141R and the resistor 142R that are coefficients of the equation (4-2), the self-inductance values L 1 and L 2 of the virtual coil 141L and the virtual coil 142L, It can be expressed by an equation including mutual inductance values M 1 and M 2 of the coil 26xL, the virtual coil 141L, and the virtual coil 142L. Specifically, the impedance of the X-axis gradient magnetic field coil 26x can be expressed by two equations for the real part and the imaginary part.

X軸傾斜磁場コイル26xのインピーダンスを表す式において、X軸傾斜磁場コイル26xにおける抵抗23xRの抵抗値Rload及びコイル23xLの自己インダクタンス値Lloadは、設計値又は測定値を用いることができる。そうすると、X軸傾斜磁場コイル26xのインピーダンスの周波数特性を表す式と、X軸傾斜磁場コイル26xのインピーダンスの周波数特性の測定値とのカーブフィッティングによって式(4-1)及び式(4-2)の係数を求めることができる。   In the equation representing the impedance of the X-axis gradient coil 26x, a design value or a measured value can be used as the resistance value Rload of the resistor 23xR and the self-inductance value Lload of the coil 23xL in the X-axis gradient coil 26x. Then, equations (4-1) and (4-2) are obtained by curve fitting between an equation representing the frequency characteristic of the impedance of the X axis gradient magnetic field coil 26x and a measured value of the frequency characteristic of the impedance of the X axis gradient magnetic field coil 26x. Can be obtained.

尚、コイル141Lを流れる電流I1(t)の初期値及びコイル142Lを流れる電流I2(t)の初期値は、過去の傾斜磁場の印加から十分な時間が経過していれば、双方ともゼロとすることができる。この結果、式(3)における未知数がなくなり、傾斜磁場の波形Gx(t)を求めることができる。 Note that the initial value of the current I 1 (t) flowing through the coil 141L and the initial value of the current I 2 (t) flowing through the coil 142L are both set if a sufficient time has elapsed since the past application of the gradient magnetic field. Can be zero. As a result, the unknown in equation (3) disappears, and the gradient magnetic field waveform Gx (t) can be obtained.

従って、パルスシーケンスにおいて定義された傾斜磁場電源27からの第1の電流Iout(t)に基づいて、第1の電流Iout(t)に対応する第1の傾斜磁場の波形Gx(t)を傾斜磁場波形シミュレーションによって計算することができる。すなわち、傾斜磁場発生システムの等価回路における仮想コイルによって生成される傾斜磁場の計算を含むシミュレーションによって取得された傾斜磁場の波形を、第1の傾斜磁場波形として用いることができる。   Therefore, the waveform Gx (t) of the first gradient magnetic field corresponding to the first current Iout (t) is tilted based on the first current Iout (t) from the gradient magnetic field power supply 27 defined in the pulse sequence. It can be calculated by magnetic field waveform simulation. In other words, the gradient magnetic field waveform obtained by the simulation including the calculation of the gradient magnetic field generated by the virtual coil in the equivalent circuit of the gradient magnetic field generation system can be used as the first gradient magnetic field waveform.

尚、非線形成分の項を含む式によって電流波形に対応する傾斜磁場波形を計算することが可能であれば、任意のモデルを用いてシミュレーションを行うことができる。換言すれば、傾斜磁場の非線形成分の項が生じるように任意の方法で傾斜磁場発生システムのモデル化を行うことができる。   In addition, if it is possible to calculate the gradient magnetic field waveform corresponding to the current waveform by an equation including the term of the nonlinear component, the simulation can be performed using an arbitrary model. In other words, the gradient magnetic field generation system can be modeled by an arbitrary method so that the term of the nonlinear component of the gradient magnetic field is generated.

傾斜磁場波形設定部40Aにおいて初期波形としての第1の電流波形及び第1の傾斜磁場波形が定まると、第1の電流波形と第1の傾斜磁場波形との間における誤差を表す指標を低減させる処理によってイメージング用の最終波形となる第2の電流波形及び第2の傾斜磁場波形を生成することが可能となる。   When the first current waveform and the first gradient magnetic field waveform as the initial waveform are determined in the gradient magnetic field waveform setting unit 40A, an index representing an error between the first current waveform and the first gradient magnetic field waveform is reduced. By processing, it is possible to generate the second current waveform and the second gradient magnetic field waveform that are the final waveform for imaging.

第1の電流波形と第1の傾斜磁場波形との間における誤差を表す指標Dindexとしては、例えば互いにスケールを合わせた第1の電流波形の面積Siと第1の傾斜磁場波形の面積Sgとの差を用いることができる。この場合における誤差の指標Dindexは、式(5)で定義することができる。
Dindex=Si-Sg (5)
As an index Dindex representing an error between the first current waveform and the first gradient magnetic field waveform, for example, an area Si of the first current waveform and an area Sg of the first gradient magnetic field waveform that are scaled with each other. Differences can be used. In this case, the error index Dindex can be defined by Equation (5).
Dindex = Si-Sg (5)

別の例として、互いにスケールを合わせた第1の電流波形の面積Siに対する第1の電流波形の面積Siと第1の傾斜磁場波形の面積Sgとの差の割合を誤差の指標Dindexとすることもできる。この場合における誤差の指標Dindexは、式(6)で定義することができる。
Dindex=(Si-Sg)/Si (6)
As another example, the ratio of the difference between the area Si of the first current waveform and the area Sg of the first gradient magnetic field waveform to the area Si of the first current waveform scaled with each other is used as the error index Dindex. You can also. In this case, the error index Dindex can be defined by Equation (6).
Dindex = (Si-Sg) / Si (6)

更に別の例として、第1の電流波形と第1の傾斜磁場波形との差として得られる誤差波形を周波数分解し、EPI法において位相方向に出現するN/2アーチファクトの原因となる高周波成分の量を誤差の指標Dindexとすることもできる。   As yet another example, the error waveform obtained as the difference between the first current waveform and the first gradient magnetic field waveform is frequency-resolved, and high frequency components that cause N / 2 artifacts appearing in the phase direction in the EPI method are obtained. The quantity can also be used as an error index Dindex.

このように任意の方法で電流波形と傾斜磁場波形との間における誤差を表す指標Dindexを定義することができる。誤差を表す指標Dindexが定義されると、指標Dindexを低減化又は最小化する計算によって第2の電流波形及び第2の傾斜磁場波形を生成することが可能となる。   In this way, the index Dindex representing the error between the current waveform and the gradient magnetic field waveform can be defined by an arbitrary method. When the index Dindex representing the error is defined, it is possible to generate the second current waveform and the second gradient magnetic field waveform by calculation that reduces or minimizes the index Dindex.

第2の電流波形及び第2の傾斜磁場波形の計算方法としては、例えば、遺伝的アルゴリズム等の最適化アルゴリズムを用いた最適化計算によって誤差の指標を最小化する方法が挙げられる。すなわち、第1の電流波形と第1の傾斜磁場波形と間における誤差の指標を最小にする最適化アルゴリズムによって第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を求めることができる。この方法は、電流に対応する傾斜磁場波形の計算を繰返す処理となるため、傾斜磁場波形シミュレーションによる傾斜磁場波形の計算が繰り返されることとなる。   Examples of the calculation method of the second current waveform and the second gradient magnetic field waveform include a method of minimizing an error index by optimization calculation using an optimization algorithm such as a genetic algorithm. That is, at least one of the second current waveform and the second gradient magnetic field waveform can be obtained by an optimization algorithm that minimizes an error index between the first current waveform and the first gradient magnetic field waveform. . Since this method repeats the calculation of the gradient magnetic field waveform corresponding to the current, the calculation of the gradient magnetic field waveform by the gradient magnetic field waveform simulation is repeated.

別の計算方法として第1の傾斜磁場波形に対する誤差を最小にするカーブフィッティングによって第2の電流の波形を求める方法が挙げられる。具体的には、指数関数の係数をパラメータとして立ち上がり期間を鈍らせた電流波形を生成し、指数関数の係数の最適化によって第1の傾斜磁場波形に近づけた電流波形として第2の電流波形を求めることができる。   As another calculation method, there is a method of obtaining the second current waveform by curve fitting that minimizes an error with respect to the first gradient magnetic field waveform. Specifically, a current waveform with a dull rise period is generated using the exponential function coefficient as a parameter, and the second current waveform is approximated to the first gradient magnetic field waveform by optimizing the exponential function coefficient. Can be sought.

更に別の計算方法として第1の傾斜磁場波形を初期の傾斜磁場波形とし、傾斜磁場波形を傾斜磁場波形シミュレーションに電流波形として入力させて傾斜磁場波形を再計算するループ計算を繰返す方法が挙げられる。このループ計算を所定回数実行すると、誤差が低減された第2の電流波形及び第2の傾斜磁場波形を得ることができる。すなわち、傾斜磁場コイル23から撮像領域に印加される傾斜磁場の波形を傾斜磁場電源27から出力される電流の波形とする計算を繰返すことによって第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を求めることができる。   As another calculation method, there is a method of repeating a loop calculation in which the first gradient magnetic field waveform is set as an initial gradient magnetic field waveform, the gradient magnetic field waveform is input as a current waveform to the gradient magnetic field waveform simulation, and the gradient magnetic field waveform is recalculated. . When this loop calculation is executed a predetermined number of times, the second current waveform and the second gradient magnetic field waveform with reduced errors can be obtained. That is, by repeating the calculation using the gradient magnetic field waveform applied from the gradient coil 23 to the imaging region as the waveform of the current output from the gradient magnetic field power supply 27, the second current waveform and the second gradient magnetic field waveform are repeated. At least one of them can be obtained.

第2の電流の波形及び第2の傾斜磁場の波形の計算に当たって必要となる誤差の許容値は、要求される画質や撮像目的に応じて任意に決定することができる。また、傾斜磁場の強度等の撮像条件によっても傾斜磁場波形の歪み量が変化するため、誤差の許容値を変えることができる。   An allowable value of an error necessary for calculating the waveform of the second current and the waveform of the second gradient magnetic field can be arbitrarily determined according to the required image quality and imaging purpose. Further, since the amount of distortion of the gradient magnetic field waveform also changes depending on the imaging conditions such as the strength of the gradient magnetic field, the allowable value of error can be changed.

一例として、N/2アーチファクト等のアーチファクトの発生を回避することが可能な電流波形と傾斜磁場波形との間における誤差の上限を予め計測する方法が挙げられる。そして、信号領域において電流波形と傾斜磁場波形との誤差がアーチファクトの回避に要求される誤差の上限を超える部分が信号領域全体の所定の割合以下となるように誤差の許容値を決定してもよい。或いは、傾斜磁場波形シミュレーションによる傾斜磁場波形の計算を繰返すことによって電流波形と傾斜磁場波形との間における誤差が収束する場合には、収束値に基づいて誤差の許容値を決定することができる。   As an example, there is a method of measuring in advance an upper limit of an error between a current waveform and a gradient magnetic field waveform that can avoid occurrence of artifacts such as N / 2 artifacts. Even if the allowable value of the error is determined so that the portion where the error between the current waveform and the gradient magnetic field waveform exceeds the upper limit of the error required for avoiding the artifact in the signal region is less than a predetermined ratio of the entire signal region. Good. Alternatively, when the error between the current waveform and the gradient magnetic field waveform converges by repeating the calculation of the gradient magnetic field waveform by the gradient magnetic field waveform simulation, the allowable value of the error can be determined based on the convergence value.

また、複数の誤差の許容値を設定することもできる。この場合には、誤差の許容値ごとに第2の電流波形と第2の傾斜磁場波形の組合せが求められることとなる。そこで、複数の誤差の許容値に対応する複数の第2の電流波形と第2の傾斜磁場波形の組合せをイメージング用の候補として表示装置34に提示し、入力装置33の操作によってユーザが選択できるようにしてもよい。   Also, a plurality of error tolerance values can be set. In this case, a combination of the second current waveform and the second gradient magnetic field waveform is obtained for each error tolerance. Therefore, a combination of a plurality of second current waveforms and second gradient magnetic field waveforms corresponding to a plurality of allowable error values is presented as a candidate for imaging on the display device 34, and can be selected by the user by operating the input device 33. You may do it.

尚、第2の電流波形及び第2の傾斜磁場波形をシミュレーション等によって求める際は、対応するパルスシーケンスを用いた撮像条件のうち最も傾斜磁場波形が歪む条件として求めることが望ましい。例えば、撮影視野(FOV: field of view)を最小とし、かつ最大の電流レベルで電流を出力させる条件が最も傾斜磁場波形が歪む条件に該当する。   When the second current waveform and the second gradient magnetic field waveform are obtained by simulation or the like, it is desirable to obtain the most distorted gradient magnetic field waveform among the imaging conditions using the corresponding pulse sequence. For example, the condition for minimizing the field of view (FOV) and outputting the current at the maximum current level corresponds to the most distorted gradient magnetic field waveform.

また、第2の電流波形の積分値は、誤差の低減処理の結果、初期波形である第1の電流波形の積分値よりも小さくなる。従って、画像データの分解能が劣化しないように、実際には、第2の電流波形の積分値が第1の電流波形の積分値と同等になるように、第2の電流波形の波高値を調整してイメージング用の撮像条件とすることが望ましい。そこで、第2の電流波形のスケーリングを行うことによって簡易に第2の電流波形の波高値及び面積を調整することができる。これは、第1の傾斜磁場波形と第2の傾斜磁場波形との関係においても同様である。   Further, the integrated value of the second current waveform becomes smaller than the integrated value of the first current waveform that is the initial waveform as a result of the error reduction process. Therefore, the peak value of the second current waveform is actually adjusted so that the integrated value of the second current waveform is equal to the integrated value of the first current waveform so that the resolution of the image data does not deteriorate. Thus, it is desirable to set the imaging conditions for imaging. Therefore, the peak value and area of the second current waveform can be easily adjusted by scaling the second current waveform. The same applies to the relationship between the first gradient magnetic field waveform and the second gradient magnetic field waveform.

第2の電流波形及び第2の傾斜磁場波形の算出は、イメージング毎に行うことができるが、磁気共鳴イメージング装置20の設計時や据付時に行うようにしてもよい。そこで、一旦求められた撮像条件や装置ごとに適切な第2の電流波形及び第2の傾斜磁場波形をライブラリ化することができる。また、傾斜磁場波形シミュレーションの結果や傾斜磁場波形の測定結果についてもライブラリ化することができる。   The calculation of the second current waveform and the second gradient magnetic field waveform can be performed for each imaging, but may be performed when the magnetic resonance imaging apparatus 20 is designed or installed. Therefore, it is possible to create a library of the second current waveform and the second gradient magnetic field waveform that are appropriate for each imaging condition and apparatus that have been obtained. Also, the results of gradient magnetic field waveform simulations and the measurement results of gradient magnetic field waveforms can be made into a library.

具体的には、傾斜磁場電源27から出力され得る互いに異なる複数の電流の波形に対応する複数の傾斜磁場の波形を傾斜磁場波形記憶部44に記憶しておくことができる。同様に、初期波形として設定された第1の電流の波形及び第1の傾斜磁場の波形に対応する第2の電流の波形及び第2の傾斜磁場の波形についても傾斜磁場波形記憶部44に保存してライブラリ化することが可能である。   Specifically, a plurality of gradient magnetic field waveforms corresponding to a plurality of different current waveforms that can be output from the gradient magnetic field power supply 27 can be stored in the gradient magnetic field waveform storage unit 44. Similarly, the waveform of the second current and the waveform of the second gradient magnetic field corresponding to the waveform of the first current and the waveform of the first gradient magnetic field set as the initial waveform are also stored in the gradient magnetic field waveform storage unit 44. And can be made into a library.

これにより、傾斜磁場波形設定部40Aは、傾斜磁場波形記憶部44から取得した第1の電流の波形に対応する第1の傾斜磁場の波形に基づいて第2の電流の波形及び第2の傾斜磁場の波形の少なくとも一方を設定又は提示することが可能となる。このため、シミュレーションに要する時間を短縮したり、傾斜磁場波形の取得のための傾斜磁場電源27からの電流波形の出力を省略することが可能となる。更に、第2の電流の波形及び第2の傾斜磁場の波形の生成のための処理自体も省略することが可能である。   As a result, the gradient magnetic field waveform setting unit 40 </ b> A generates the second current waveform and the second gradient based on the first gradient magnetic field waveform corresponding to the first current waveform acquired from the gradient magnetic field waveform storage unit 44. It becomes possible to set or present at least one of the waveforms of the magnetic field. For this reason, it is possible to shorten the time required for the simulation and to omit the output of the current waveform from the gradient magnetic field power supply 27 for acquiring the gradient magnetic field waveform. Furthermore, the processing itself for generating the second current waveform and the second gradient magnetic field waveform can be omitted.

一方、コンピュータ32のデータ処理部41は、撮像条件設定部40において設定された撮像条件下におけるイメージングスキャンによって収集されたMR信号をシーケンスコントローラ31から取得してk空間データ記憶部42に形成されたk空間に配置する機能、k空間データ記憶部42からk空間データを取り込んでフーリエ変換(FT: Fourier transform)を含む画像再構成処理を施すことにより画像データを再構成する機能、再構成して得られた画像データを画像データ記憶部43に書き込む機能、画像データ記憶部43から取り込んだ画像データに必要な画像処理を施して表示装置34に表示させる機能を有する。   On the other hand, the data processing unit 41 of the computer 32 acquires the MR signals collected by the imaging scan under the imaging conditions set in the imaging condition setting unit 40 from the sequence controller 31 and is formed in the k-space data storage unit 42. a function to arrange in k-space, a function to reconstruct image data by taking k-space data from the k-space data storage unit 42 and performing image reconstruction processing including Fourier transform (FT). A function of writing the obtained image data in the image data storage unit 43 and a function of performing necessary image processing on the image data captured from the image data storage unit 43 and displaying the image data on the display device 34 are provided.

次に磁気共鳴イメージング装置20の動作及び作用について説明する。   Next, the operation and action of the magnetic resonance imaging apparatus 20 will be described.

図8は、図1に示す磁気共鳴イメージング装置20により傾斜磁場電源27からの出力電流波形と傾斜磁場波形との間における誤差を考慮してイメージングを実行する際の流れを示すフローチャートである。   FIG. 8 is a flowchart showing a flow when imaging is performed in consideration of an error between the output current waveform from the gradient magnetic field power supply 27 and the gradient magnetic field waveform by the magnetic resonance imaging apparatus 20 shown in FIG.

まずステップS1において、撮像条件設定部40においてEPIシーケンス等のパルスシーケンスを含む撮像条件が設定される。   First, in step S1, the imaging condition setting unit 40 sets imaging conditions including a pulse sequence such as an EPI sequence.

次に、ステップS11において、傾斜磁場波形設定部40Aは、傾斜磁場波形と電流波形とが相似するとの前提で設定されたパルスシーケンスにおける傾斜磁場の制御値を、傾斜磁場を生成するための初期の電流波形である第1の電流波形として設定する。   Next, in step S11, the gradient magnetic field waveform setting unit 40A uses the gradient magnetic field control value in the pulse sequence set on the assumption that the gradient magnetic field waveform and the current waveform are similar to each other as an initial value for generating the gradient magnetic field. The first current waveform that is a current waveform is set.

次に、ステップS12において、傾斜磁場波形設定部40Aは、初期設定された第1の電流波形に対応する第1の傾斜磁場波形を取得する。第1の傾斜磁場波形は、傾斜磁場波形シミュレーション部40Bにおける、第1の電流波形を入力とする傾斜磁場波形シミュレーションによって求めることができる。   Next, in step S12, the gradient magnetic field waveform setting unit 40A obtains a first gradient magnetic field waveform corresponding to the initially set first current waveform. The first gradient magnetic field waveform can be obtained by a gradient magnetic field waveform simulation in which the first current waveform is input in the gradient magnetic field waveform simulation unit 40B.

尚、撮像断面がオブリーク断面である場合のように、撮像条件や磁気共鳴イメージング装置20の状態によって傾斜磁場波形の電流波形に対する歪量が変化する場合がある。例えば、傾斜磁場の印加方向や強度に応じて傾斜磁場波形の歪量が変化する場合がある。   Note that the amount of distortion of the gradient magnetic field waveform with respect to the current waveform may change depending on the imaging conditions and the state of the magnetic resonance imaging apparatus 20 as in the case where the imaging cross section is an oblique cross section. For example, the distortion amount of the gradient magnetic field waveform may change depending on the application direction and strength of the gradient magnetic field.

そこで、磁気共鳴イメージング装置20の状態を表す情報や撮像条件をパラメータとする傾斜磁場波形シミュレーションによって磁気共鳴イメージング装置20の状態及び撮像条件に応じたより正確な傾斜磁場波形を求めることができる。   Therefore, a more accurate gradient magnetic field waveform corresponding to the state of the magnetic resonance imaging apparatus 20 and the imaging conditions can be obtained by a gradient magnetic field waveform simulation using information representing the state of the magnetic resonance imaging apparatus 20 and imaging conditions as parameters.

或いは、傾斜磁場波形シミュレーションに依らず、第1の電流波形を傾斜磁場電源27から出力させるプレスキャンを行って傾斜磁場波形を取得するようにしてもよい。具体的には、サーチコイル50及び積分器51で構成される傾斜磁場波形の検出システムを用いて傾斜磁場波形を検出したり、プレスキャンによって傾斜磁場の印加期間中に収集されたFID信号に基づく演算によって傾斜磁場波形を求めることができる。この場合には、傾斜磁場波形取得部40Cにより第1の電流波形に対応する第1の傾斜磁場波形が取得される。   Alternatively, the gradient magnetic field waveform may be acquired by performing a pre-scan for outputting the first current waveform from the gradient magnetic field power supply 27 without depending on the gradient magnetic field waveform simulation. Specifically, a gradient magnetic field waveform is detected using a gradient magnetic field waveform detection system including a search coil 50 and an integrator 51, or based on a FID signal collected during pre-scanning application of the gradient magnetic field. A gradient magnetic field waveform can be obtained by calculation. In this case, the first gradient magnetic field waveform corresponding to the first current waveform is acquired by the gradient magnetic field waveform acquisition unit 40C.

次に、ステップS13において、傾斜磁場波形設定部40Aは、第1の電流波形と第1の傾斜磁場波形との間における誤差を算出する。次に、ステップS14において、傾斜磁場波形設定部40Aは、第1の電流波形と第1の傾斜磁場波形との間における誤差が誤差の許容値として予め決定された閾値以下であるか否かを判定する。   Next, in step S13, the gradient magnetic field waveform setting unit 40A calculates an error between the first current waveform and the first gradient magnetic field waveform. Next, in step S14, the gradient magnetic field waveform setting unit 40A determines whether or not an error between the first current waveform and the first gradient magnetic field waveform is equal to or less than a threshold value determined in advance as an error tolerance. judge.

そして、第1の電流波形と第1の傾斜磁場波形との間における誤差が閾値以下ではないと判定された場合には、ステップS15において、傾斜磁場波形設定部40Aが電流波形を変更する。具体的には、電流波形をパラメータとして誤差を最小にする最適化アルゴリズムや指数関数による傾斜磁場波形へのカーブフィッティングによって電流波形が変更される。或いは、スケールを合わせた傾斜磁場波形自体を電流波形に設定するようにしてもよい。   If it is determined that the error between the first current waveform and the first gradient magnetic field waveform is not less than or equal to the threshold, the gradient magnetic field waveform setting unit 40A changes the current waveform in step S15. Specifically, the current waveform is changed by an optimization algorithm that minimizes an error using the current waveform as a parameter or curve fitting to a gradient magnetic field waveform by an exponential function. Or you may make it set the gradient magnetic field waveform itself which matched the scale to a current waveform.

そして、再びステップS12において変更された電流波形に対応する傾斜磁場波形が任意の方法で取得される。但し、傾斜磁場波形の取得が繰返される場合には、傾斜磁場波形シミュレーションにより計算で求めることが現実的である。   Then, a gradient magnetic field waveform corresponding to the current waveform changed in step S12 is acquired again by an arbitrary method. However, when the acquisition of the gradient magnetic field waveform is repeated, it is practical to obtain it by calculation using a gradient magnetic field waveform simulation.

そして、上述したような電流波形の変更及び変更された傾斜磁場波形の取得は、ステップS14において電流波形と傾斜磁場波形との間における誤差が閾値以下であると判定されるまで繰返される。この結果、乖離量が閾値以下となった第2の電流波形と第2の傾斜磁場波形の組合せを得ることができる。   Then, the change of the current waveform and the acquisition of the changed gradient magnetic field waveform as described above are repeated until it is determined in step S14 that the error between the current waveform and the gradient magnetic field waveform is equal to or less than the threshold value. As a result, it is possible to obtain a combination of the second current waveform and the second gradient magnetic field waveform in which the deviation amount is equal to or less than the threshold value.

そして、得られた第2の電流波形と第2の傾斜磁場波形の組合せは、ステップS14において撮像条件として提示又は設定される。また、撮像条件設定部40において、第2の電流波形の撮像条件としての確定処理や他の必要な撮像条件の設定を行うことができる。   Then, the obtained combination of the second current waveform and the second gradient magnetic field waveform is presented or set as an imaging condition in step S14. Further, the imaging condition setting unit 40 can perform a determination process as an imaging condition of the second current waveform and other necessary imaging conditions.

撮像条件設定部40において撮像条件の設定が完了すると、ステップS2においてシーケンスコントローラ31や静磁場用磁石21等を含むイメージング系によってMRイメージングが実行される。   When the setting of the imaging condition is completed in the imaging condition setting unit 40, MR imaging is executed by the imaging system including the sequence controller 31, the static magnetic field magnet 21 and the like in step S2.

具体的には予め寝台37に被検体Pがセットされ、静磁場電源26により励磁された静磁場用磁石21(超伝導磁石)の撮像領域に静磁場が形成される。また、シムコイル電源28からシムコイル22に電流が供給されて撮像領域に形成された静磁場が均一化される。   Specifically, the subject P is set on the bed 37 in advance, and a static magnetic field is formed in the imaging region of the static magnetic field magnet 21 (superconducting magnet) excited by the static magnetic field power supply 26. Further, a current is supplied from the shim coil power supply 28 to the shim coil 22, and the static magnetic field formed in the imaging region is made uniform.

そして、入力装置33から撮像条件設定部40にスキャン開始指示が与えられると、撮像条件設定部40はパルスシーケンスを含む撮影条件をシーケンスコントローラ31に出力する。このため、シーケンスコントローラ31は、パルスシーケンスに従って傾斜磁場電源27、送信器29及び受信器30を駆動させることにより被検体Pがセットされた撮像領域に傾斜磁場を形成させるとともに、RFコイル24からRF信号を発生させる。   Then, when a scan start instruction is given from the input device 33 to the imaging condition setting unit 40, the imaging condition setting unit 40 outputs imaging conditions including a pulse sequence to the sequence controller 31. For this reason, the sequence controller 31 drives the gradient magnetic field power source 27, the transmitter 29, and the receiver 30 according to the pulse sequence to form a gradient magnetic field in the imaging region where the subject P is set, and from the RF coil 24 to the RF Generate a signal.

このため、被検体Pの内部における磁気共鳴により生じたMR信号が、RFコイル24により受信されて受信器30に与えられる。受信器30は、RFコイル24からMR信号を受けて、所要の信号処理を実行した後、A/D変換することにより、デジタルデータのMR信号である生データを生成する。受信器30は、生成した生データをシーケンスコントローラ31に与える。シーケンスコントローラ31は、生データをコンピュータ32のデータ処理部41に出力する。   For this reason, the MR signal generated by the magnetic resonance inside the subject P is received by the RF coil 24 and given to the receiver 30. The receiver 30 receives the MR signal from the RF coil 24, performs necessary signal processing, and then performs A / D conversion to generate raw data that is an MR signal of digital data. The receiver 30 gives the generated raw data to the sequence controller 31. The sequence controller 31 outputs the raw data to the data processing unit 41 of the computer 32.

そうすると、データ処理部41は、イメージングスキャンによって収集されたMR信号をk空間データ記憶部42に形成されたk空間に配置する。続いて、データ処理部41はk空間データ記憶部42からk空間データを取り込んで画像再構成処理を施すことにより画像データを生成する。そして、画像データに必要な画像処理を施して表示装置34に表示させたり、画像データ記憶部43に書き込んで保存することができる。   Then, the data processing unit 41 arranges MR signals collected by the imaging scan in the k space formed in the k space data storage unit 42. Subsequently, the data processing unit 41 takes in k-space data from the k-space data storage unit 42 and performs image reconstruction processing to generate image data. Then, the image data can be subjected to necessary image processing and displayed on the display device 34 or can be written and stored in the image data storage unit 43.

尚、図8には、傾斜磁場波形との誤差が小さくなるような傾斜磁場電源27の出力電流波形をイメージング毎に事前に求める場合の流れを示したが、磁気共鳴イメージング装置20の設計時や据付時において傾斜磁場波形との誤差が小さくなるような電流波形を求めるようにしてもよい。   FIG. 8 shows a flow in the case where the output current waveform of the gradient magnetic field power supply 27 is obtained in advance for each imaging so that an error from the gradient magnetic field waveform becomes small. You may make it obtain | require a current waveform that an error with a gradient magnetic field waveform becomes small at the time of installation.

磁気共鳴イメージング装置20やパルスシーケンスの設計時に傾斜磁場波形との誤差が小さくなるようなイメージング用の第2の電流波形を求める場合には、傾斜磁場電源27及び傾斜磁場コイル23を含む磁気共鳴イメージング装置20の設計値に基づいて傾斜磁場発生システムの特性を表すパラメータを見積もることができる。そして、見積もられたパラメータ値を用いた傾斜磁場波形シミュレーションによって傾斜磁場波形との誤差が小さくなるようなイメージング用の第2の電流波形をパルスシーケンス毎に計算することができる。更に、計算されたパルスシーケンス毎のイメージング用の第2の電流波形を傾斜磁場波形記憶部44に記憶させ、イメージングの際にパルスシーケンスに対応する第2の電流波形及び第2の電流波形に対応する第2の傾斜磁場波形を撮像条件として用いることができる。   When obtaining the second current waveform for imaging so that an error from the gradient magnetic field waveform becomes small when designing the magnetic resonance imaging apparatus 20 or the pulse sequence, the magnetic resonance imaging including the gradient magnetic field power supply 27 and the gradient magnetic field coil 23 is used. Parameters representing the characteristics of the gradient magnetic field generation system can be estimated based on the design values of the apparatus 20. Then, a second current waveform for imaging that can reduce an error from the gradient magnetic field waveform by the gradient magnetic field waveform simulation using the estimated parameter value can be calculated for each pulse sequence. Further, the calculated second current waveform for imaging for each pulse sequence is stored in the gradient magnetic field waveform storage unit 44, and corresponds to the second current waveform and the second current waveform corresponding to the pulse sequence during imaging. The second gradient magnetic field waveform to be used can be used as the imaging condition.

一方、磁気共鳴イメージング装置20の据付時においても同様に傾斜磁場波形との誤差が小さくなるようなイメージング用の第2の電流波形を求めて傾斜磁場波形記憶部44に保存しておくことができる。但し、傾斜磁場コイル23のインダクタンスやキャパシタンス等の傾斜磁場発生システムの設計誤差のばらつきに応じた適切な第2の電流波形を求めることができる。その場合には、傾斜磁場コイル23等のハードウェアの特性のばらつきを模擬した傾斜磁場波形シミュレーションによって第2の電流波形が計算されることとなる。或いは、サーチコイル50及び積分器51で構成される傾斜磁場波形の検出システムを用いた傾斜磁場波形の検出やMR信号に基づく傾斜磁場波形の取得が必要となる。   On the other hand, when the magnetic resonance imaging apparatus 20 is installed, a second current waveform for imaging that similarly reduces the error from the gradient magnetic field waveform can be obtained and stored in the gradient magnetic field waveform storage unit 44. . However, it is possible to obtain an appropriate second current waveform corresponding to variations in the design error of the gradient magnetic field generation system such as the inductance and capacitance of the gradient magnetic field coil 23. In this case, the second current waveform is calculated by a gradient magnetic field waveform simulation that simulates the variation in the characteristics of hardware such as the gradient coil 23. Alternatively, it is necessary to detect a gradient magnetic field waveform using a gradient magnetic field waveform detection system including a search coil 50 and an integrator 51 and to acquire a gradient magnetic field waveform based on an MR signal.

傾斜磁場波形シミュレーションを実行する場合には、傾斜磁場波形シミュレーションによって計算された傾斜磁場波形と、傾斜磁場波形の検出システムによって検出された傾斜磁場波形又はMR信号に基づいて取得された傾斜磁場波形とを比較することによって傾斜磁場波形シミュレーションに用いられるパラメータを傾斜磁場発生システムの特性に適合させることができる。そして、パルスシーケンス毎に傾斜磁場波形が最も歪む条件で傾斜磁場波形シミュレーションを実行することにより、傾斜磁場発生システムの特性のばらつきに応じた第2の電流波形及び第2の電流波形に対応する第2の傾斜磁場波形を求めることができる。   When executing the gradient magnetic field waveform simulation, the gradient magnetic field waveform calculated by the gradient magnetic field waveform simulation, the gradient magnetic field waveform detected based on the gradient magnetic field waveform detected by the gradient magnetic field waveform detection system or the MR signal, and The parameters used in the gradient magnetic field waveform simulation can be adapted to the characteristics of the gradient magnetic field generation system. Then, by executing the gradient magnetic field waveform simulation under the condition that the gradient magnetic field waveform is most distorted for each pulse sequence, the second current waveform corresponding to the variation in the characteristics of the gradient magnetic field generation system and the second current waveform corresponding to the second current waveform are obtained. 2 gradient magnetic field waveforms can be obtained.

つまり以上のような磁気共鳴イメージング装置20は、初期波形として設定された傾斜磁場電源27から第1の電流波形に対応する第1の傾斜磁場波形を求め、電流波形と傾斜磁場波形との間における誤差が低減されるように補正された第2の電流波形及び第2の傾斜磁場波形をイメージング用に提示又は設定するようにしたものである。   That is, the magnetic resonance imaging apparatus 20 as described above obtains the first gradient magnetic field waveform corresponding to the first current waveform from the gradient magnetic field power supply 27 set as the initial waveform, and between the current waveform and the gradient magnetic field waveform. The second current waveform and the second gradient magnetic field waveform corrected so as to reduce the error are presented or set for imaging.

このため、磁気共鳴イメージング装置20によれば、傾斜磁場電源27から出力される電流波形により近い波形を有する傾斜磁場を撮像領域に印加することが可能となる。従って、傾斜磁場波形を台形や矩形に近づけるための急激な立ち上りを有する電流の出力を回避すすることができる。この結果、傾斜磁場アンプへの負荷が減り、傾斜磁場発生システムの性能を最大限引き出すことが可能となる。   For this reason, according to the magnetic resonance imaging apparatus 20, it is possible to apply a gradient magnetic field having a waveform closer to the current waveform output from the gradient magnetic field power supply 27 to the imaging region. Therefore, it is possible to avoid output of a current having a sharp rise for making the gradient magnetic field waveform close to a trapezoid or a rectangle. As a result, the load on the gradient magnetic field amplifier is reduced, and the performance of the gradient magnetic field generation system can be maximized.

特に、EPIにおいてETL (echo train length)に課される制限を緩和することができる。また、EPIにおけるregridding処理を、より正確に行うことが可能となる。従って、画質の劣化やアーチファクトを低減することができる。   In particular, restrictions imposed on ETL (echo train length) in EPI can be relaxed. In addition, regridding processing in EPI can be performed more accurately. Therefore, image quality degradation and artifacts can be reduced.

以上、特定の実施形態について記載したが、記載された実施形態は一例に過ぎず、発明の範囲を限定するものではない。ここに記載された新規な方法及び装置は、様々な他の様式で具現化することができる。また、ここに記載された方法及び装置の様式において、発明の要旨から逸脱しない範囲で、種々の省略、置換及び変更を行うことができる。添付された請求の範囲及びその均等物は、発明の範囲及び要旨に包含されているものとして、そのような種々の様式及び変形例を含んでいる。   Although specific embodiments have been described above, the described embodiments are merely examples, and do not limit the scope of the invention. The novel methods and apparatus described herein can be implemented in a variety of other ways. Various omissions, substitutions, and changes can be made in the method and apparatus described herein without departing from the spirit of the invention. The appended claims and their equivalents include such various forms and modifications as are encompassed by the scope and spirit of the invention.

20 磁気共鳴イメージング装置
21 静磁場用磁石
22 シムコイル
23 傾斜磁場コイル
24 RFコイル
25 制御系
26 静磁場電源
27 傾斜磁場電源
28 シムコイル電源
29 送信器
30 受信器
31 シーケンスコントローラ
32 コンピュータ
33 入力装置
34 表示装置
35 演算装置
36 記憶装置
37 寝台
40 撮像条件設定部
41 データ処理部
42 k空間データ記憶部
43 画像データ記憶部
44 傾斜磁場波形記憶部
50 サーチコイル
51 積分器
52 増幅器
53 A/D変換器
P 被検体
20 Magnetic Resonance Imaging Device 21 Magnet for Static Magnetic Field 22 Shim Coil 23 Gradient Magnetic Field Coil 24 RF Coil 25 Control System 26 Static Magnetic Field Power Supply 27 Gradient Magnetic Field Power Supply 28 Shim Coil Power Supply 29 Transmitter 30 Receiver 31 Sequence Controller 32 Computer 33 Input Device 34 Display Device 35 arithmetic unit 36 storage unit 37 bed 40 imaging condition setting unit 41 data processing unit 42 k-space data storage unit 43 image data storage unit 44 gradient magnetic field waveform storage unit 50 search coil 51 integrator 52 amplifier 53 A / D converter P Specimen

Claims (10)

傾斜磁場コイル及び傾斜磁場電源を含む傾斜磁場発生システムを用いて被検体の磁気共鳴イメージングを行うイメージング手段と、
前記傾斜磁場電源から出力される第1の電流の波形に応じて撮像領域に印加される第1の傾斜磁場の波形に基づいて、前記磁気共鳴イメージング用に前記傾斜磁場電源から出力される第2の電流の波形及び前記撮像領域に印加される第2の傾斜磁場の波形の少なくとも一方を設定又は提示する条件設定手段と、
を備え
前記条件設定手段は、前記傾斜磁場電源から出力される電流の波形と前記撮像領域に印加される傾斜磁場の波形との間における誤差を表す指標が小さくなるように前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される磁気共鳴イメージング装置。
An imaging means for performing magnetic resonance imaging of a subject using a gradient magnetic field generation system including a gradient coil and a gradient magnetic field power supply;
Based on the waveform of the first gradient magnetic field applied to the imaging region in accordance with the waveform of the first current output from the gradient magnetic field power supply, the second output from the gradient magnetic field power supply for the magnetic resonance imaging. Condition setting means for setting or presenting at least one of the waveform of the current and the waveform of the second gradient magnetic field applied to the imaging region;
Equipped with a,
The condition setting means includes the second current waveform and the second current waveform so that an index indicating an error between the waveform of the current output from the gradient magnetic field power supply and the waveform of the gradient magnetic field applied to the imaging region is small. magnetic resonance imaging apparatus that will be configured to determine at least one of said second gradient magnetic field waveform.
傾斜磁場コイル及び傾斜磁場電源を含む傾斜磁場発生システムを用いて被検体の磁気共鳴イメージングを行うイメージング手段と、
前記傾斜磁場電源から出力される第1の電流の波形に応じて撮像領域に印加される第1の傾斜磁場の波形に基づいて、前記磁気共鳴イメージング用に前記傾斜磁場電源から出力される第2の電流の波形及び前記撮像領域に印加される第2の傾斜磁場の波形の少なくとも一方を設定又は提示する条件設定手段と、
を備え、
前記条件設定手段は、前記第1の傾斜磁場の波形として、前記第1の電流の波形に基づくシミュレーションによって取得された傾斜磁場の波形を用いて、前記傾斜磁場電源から出力される電流の波形と前記撮像領域に印加される傾斜磁場の波形との間における誤差を表す指標が小さくなるように前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される磁気共鳴イメージング装置。
An imaging means for performing magnetic resonance imaging of a subject using a gradient magnetic field generation system including a gradient coil and a gradient magnetic field power supply;
Based on the waveform of the first gradient magnetic field applied to the imaging region in accordance with the waveform of the first current output from the gradient magnetic field power supply, the second output from the gradient magnetic field power supply for the magnetic resonance imaging. Condition setting means for setting or presenting at least one of the waveform of the current and the waveform of the second gradient magnetic field applied to the imaging region;
With
The condition setting means uses a gradient magnetic field waveform obtained by simulation based on the first current waveform as the first gradient magnetic field waveform, and a waveform of a current output from the gradient magnetic field power supply. It is configured to determine at least one of the waveform and the second gradient of the waveform of the second current as an index representing the error becomes smaller between the waveform of a gradient magnetic field to be applied to the imaging area Magnetic resonance imaging device.
前記条件設定手段は、前記第1の傾斜磁場の波形として、前記第1の電流の波形を前記傾斜磁場電源から出力させることによって測定された傾斜磁場の波形を用いるように構成される請求項記載の磁気共鳴イメージング装置。 Said condition setting means, wherein a first waveform of the gradient magnetic field, the first current waveform the inclination claim configured to use the measured waveform of the gradient magnetic field by the output from the magnetic field supply 1 The magnetic resonance imaging apparatus described. 前記条件設定手段は、前記第1の傾斜磁場の波形として、前記第1の電流の波形を前記傾斜磁場電源から出力させることによって収集された磁気共鳴信号に基づいて求められた傾斜磁場の波形を用いるように構成される請求項記載の磁気共鳴イメージング装置。 The condition setting means includes a gradient magnetic field waveform obtained based on a magnetic resonance signal collected by outputting the first current waveform from the gradient magnetic field power source as the first gradient magnetic field waveform. magnetic resonance imaging apparatus according to claim 1, wherein configured to use. 前記条件設定手段は、前記第1の傾斜磁場の波形として、前記傾斜磁場発生システムの等価回路における仮想コイルによって生成される傾斜磁場の計算を含むシミュレーションによって取得された傾斜磁場の波形を用いるように構成される請求項記載の磁気共鳴イメージング装置。 The condition setting means uses a gradient magnetic field waveform obtained by a simulation including a calculation of a gradient magnetic field generated by a virtual coil in an equivalent circuit of the gradient magnetic field generation system as the waveform of the first gradient magnetic field. The magnetic resonance imaging apparatus according to claim 2 configured. 前記傾斜磁場電源から前記第1の電流の波形を有する電流を出力させることによって前記第1の傾斜磁場の波形を有する傾斜磁場を前記撮像領域に印加する傾斜磁場印加手段を更に備え、
前記条件設定手段は、前記撮像領域に印加された前記第1の傾斜磁場の波形に基づいて前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を設定又は提示するように構成される請求項1乃至のいずれか1項に記載の磁気共鳴イメージング装置。
A gradient magnetic field applying means for applying a gradient magnetic field having the waveform of the first gradient magnetic field to the imaging region by outputting a current having the waveform of the first current from the gradient magnetic field power supply;
The condition setting means sets or presents at least one of the waveform of the second current and the waveform of the second gradient magnetic field based on the waveform of the first gradient magnetic field applied to the imaging region. magnetic resonance imaging apparatus according to any one of constituted claims 1 to 5.
前記傾斜磁場電源から出力され得る互いに異なる複数の電流の波形に対応する複数の傾斜磁場の波形を記憶する記憶手段を更に備え、
前記条件設定手段は、前記記憶手段から取得した前記第1の電流の波形に対応する前記第1の傾斜磁場の波形に基づいて前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を設定又は提示するように構成される請求項1乃至のいずれか1項に記載の磁気共鳴イメージング装置。
A storage means for storing a plurality of gradient magnetic field waveforms corresponding to a plurality of different current waveforms that can be output from the gradient magnetic field power supply;
The condition setting means includes a waveform of the second current and a waveform of the second gradient magnetic field based on the waveform of the first gradient magnetic field corresponding to the waveform of the first current acquired from the storage means. magnetic resonance imaging apparatus according to any one of claims 1 to 5 arranged to set or display at least one.
前記条件設定手段は、前記指標を最小にする最適化アルゴリズムによって前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される請求項1又は2記載の磁気共鳴イメージング装置。 Said condition setting means, the optimization algorithm by said second current waveform and the magnetic of the second at least one of the ask configured claim 1 or 2, wherein the waveform of a gradient magnetic field for the index to the minimum Resonance imaging device. 前記条件設定手段は、前記第1の傾斜磁場の波形に対する誤差を最小にするカーブフィッティングによって前記第2の電流の波形を求めるように構成される請求項1又は2記載の磁気共鳴イメージング装置。 Said condition setting means, the first magnetic resonance imaging apparatus according to claim 1 or 2, wherein configured by curve fitting to minimize the error with respect to the waveform of a gradient magnetic field seek waveform of the second current. 前記条件設定手段は、前記傾斜磁場コイルから前記撮像領域に印加される傾斜磁場の波形を前記傾斜磁場電源から出力される電流の波形とする計算を繰返すことによって前記第2の電流の波形及び前記第2の傾斜磁場の波形の少なくとも一方を求めるように構成される請求項1乃至のいずれか1項に記載の磁気共鳴イメージング装置。 The condition setting means repeats the calculation using the gradient magnetic field waveform applied from the gradient magnetic field coil to the imaging region as the waveform of the current output from the gradient magnetic field power source, and the second current waveform and the magnetic resonance imaging apparatus according to any one of the second gradient magnetic field of claims configured to determine at least one waveform to claim 1 to 8.
JP2012188305A 2012-08-29 2012-08-29 Magnetic resonance imaging system Expired - Fee Related JP6109508B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012188305A JP6109508B2 (en) 2012-08-29 2012-08-29 Magnetic resonance imaging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012188305A JP6109508B2 (en) 2012-08-29 2012-08-29 Magnetic resonance imaging system

Publications (2)

Publication Number Publication Date
JP2014045775A JP2014045775A (en) 2014-03-17
JP6109508B2 true JP6109508B2 (en) 2017-04-05

Family

ID=50606104

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012188305A Expired - Fee Related JP6109508B2 (en) 2012-08-29 2012-08-29 Magnetic resonance imaging system

Country Status (1)

Country Link
JP (1) JP6109508B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9971008B2 (en) * 2014-09-30 2018-05-15 Toshiba Medical Systems Corporation MRI gradient trajectory mapping
WO2016174497A1 (en) * 2015-04-27 2016-11-03 Synaptive Medical (Barbados) Inc. System and method for image warp correction for magnetic resonance imaging
JP6817775B2 (en) 2016-10-11 2021-01-20 株式会社東芝 Correction device, correction method and magnetic resonance imaging device

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6464639A (en) * 1987-09-07 1989-03-10 Hitachi Medical Corp Nuclear magnetic resonance imaging apparatus
JPH02211124A (en) * 1989-02-13 1990-08-22 Hitachi Medical Corp Magnetic resonance imaging device
JPH03173531A (en) * 1989-12-01 1991-07-26 Hitachi Ltd Method and device for eddy current compensation of magnetic resonance imaging device
JPH07148136A (en) * 1993-11-30 1995-06-13 Shimadzu Corp Mr imaging device
JP3576658B2 (en) * 1995-10-24 2004-10-13 ジーイー横河メディカルシステム株式会社 Driving waveform generation device for gradient magnetic field, eddy current estimation device, and MRI device
JPH1189817A (en) * 1997-09-22 1999-04-06 Hitachi Medical Corp Magnetic resonance imaging equipment
JP3921369B2 (en) * 2001-10-02 2007-05-30 株式会社日立製作所 Inspection apparatus using nuclear magnetic resonance and method for adjusting gradient magnetic field waveform
JP4106053B2 (en) * 2002-07-04 2008-06-25 株式会社日立メディコ Magnetic resonance imaging apparatus and eddy current compensation derivation method
JP3930439B2 (en) * 2003-02-06 2007-06-13 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Eddy current correction method and magnetic resonance imaging apparatus
US7112964B2 (en) * 2004-08-02 2006-09-26 General Electric Company Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom
JP5666779B2 (en) * 2009-01-27 2015-02-12 株式会社東芝 Magnetic resonance imaging system
US9588206B2 (en) * 2011-06-30 2017-03-07 Hitachi, Ltd. Magnetic resonance imaging apparatus and high-frequency magnetic field dertermination method

Also Published As

Publication number Publication date
JP2014045775A (en) 2014-03-17

Similar Documents

Publication Publication Date Title
WO2012173095A1 (en) Magnetic resonance imaging device and control device thereof
US20230324489A1 (en) Eddy current mitigation systems and methods
JP5916240B2 (en) Magnetic resonance imaging apparatus and gradient magnetic field waveform estimation method
JP5611661B2 (en) Magnetic resonance imaging system
US9720066B2 (en) Magnetic resonance imaging apparatus and control method thereof
US9317917B2 (en) Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
JP5886024B2 (en) Magnetic resonance imaging system
JP2012040362A (en) Magnetic resonance imaging method, magnetic resonance imaging apparatus, and control device of magnetic resonance imaging apparatus
US8362773B2 (en) System and method for modeling gradient coil operation induced magnetic field drift
JP5865785B2 (en) Magnetic resonance imaging apparatus and gradient magnetic field generation system load calculation method
JP2013000173A (en) Magnetic resonance imaging apparatus and control device of the same
EP3413071A1 (en) Mr imaging using a stack-of-stars acquisition with variable contrast
US7689015B2 (en) Magnetic resonance imaging apparatus and image correction estimating method
JP2018514277A (en) System and method for correcting image distortion in nuclear magnetic resonance imaging
JP5638393B2 (en) Magnetic resonance imaging apparatus and method
JP6552804B2 (en) Magnetic resonance imaging system
JP6109508B2 (en) Magnetic resonance imaging system
US9547062B2 (en) Method and magnetic resonance apparatus to generate an artifact-free magnetic resonance image data set
JP5846450B2 (en) Magnetic resonance imaging apparatus, measurement space coordinate correction method, and image reconstruction method
KR20130135776A (en) Method to determine a subject-specific b1 distribution of an examination subject in a measurement volume in the magnetic resonance technique, magnetic resonance system, computer program, and electronically readable data medium
JP2013017811A (en) Magnetic resonance imaging apparatus
JPH0751249A (en) Method and device for inspecting object of measurement by locally resolving magnetic resonance
JP5718148B2 (en) Magnetic resonance imaging apparatus and dual slice measurement method
JP5971689B2 (en) Magnetic resonance imaging device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150714

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20160510

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160512

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160614

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160812

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170308

R150 Certificate of patent or registration of utility model

Ref document number: 6109508

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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