WO2013002233A1 - 磁気共鳴イメージング装置および傾斜磁場波形推定方法 - Google Patents

磁気共鳴イメージング装置および傾斜磁場波形推定方法 Download PDF

Info

Publication number
WO2013002233A1
WO2013002233A1 PCT/JP2012/066316 JP2012066316W WO2013002233A1 WO 2013002233 A1 WO2013002233 A1 WO 2013002233A1 JP 2012066316 W JP2012066316 W JP 2012066316W WO 2013002233 A1 WO2013002233 A1 WO 2013002233A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic field
gradient magnetic
waveform
response function
output
Prior art date
Application number
PCT/JP2012/066316
Other languages
English (en)
French (fr)
Inventor
光 花田
佐藤　善隆
甲亮 平井
邦治 岡
瀧澤　将宏
直哉 坂口
英久 秋丸
美由紀 河村
Original Assignee
株式会社 日立メディコ
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 株式会社 日立メディコ filed Critical 株式会社 日立メディコ
Priority to US14/128,051 priority Critical patent/US9664765B2/en
Priority to CN201280029381.5A priority patent/CN103596496B/zh
Priority to JP2013522881A priority patent/JP5916240B2/ja
Publication of WO2013002233A1 publication Critical patent/WO2013002233A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging

Definitions

  • the present invention measures a nuclear magnetic resonance signal (hereinafter referred to as an NMR signal) from hydrogen, phosphorus, etc. in a subject and images a nuclear density distribution, relaxation time distribution, etc. (hereinafter referred to as MRI).
  • NMR signal nuclear magnetic resonance signal
  • MRI nuclear density distribution, relaxation time distribution, etc.
  • the present invention relates to a technique for compensating for image quality deterioration caused by gradient magnetic field distortion.
  • the MRI device measures NMR signals generated by the spins of the subject, especially the tissues of the human body, and visualizes the form and function of the head, abdomen, limbs, etc. in two or three dimensions Device.
  • a subject is placed in a static magnetic field (polarization magnetic field B0), a high-frequency magnetic field pulse is applied together with a slice selective gradient magnetic field pulse to selectively excite a specific region, and then a phase encoding gradient magnetic field pulse is applied.
  • polarization magnetic field B0 polarization magnetic field B0
  • a high-frequency magnetic field pulse is applied together with a slice selective gradient magnetic field pulse to selectively excite a specific region, and then a phase encoding gradient magnetic field pulse is applied.
  • a readout gradient magnetic field pulse the excitation range is encoded and position information is given.
  • Gradient magnetic field pulses provide a desired gradient magnetic field by supplying a pulsed current (referred to as an input gradient magnetic field waveform) from a gradient magnetic field power source to a plurality of coils that generate linear gradient magnetic fields in three orthogonal axes.
  • a pulse (referred to as an output gradient magnetic field waveform) is output.
  • the input gradient magnetic field waveform and the gradient magnetic field waveform generated in the MRI apparatus should coincide with each other, but errors occur due to various factors.
  • One of the causes of errors is the effect of eddy currents generated in magnetic parts or electric circuits constituting the MRI apparatus due to a sudden change in magnetic field due to a gradient magnetic field pulse.
  • the eddy current gives a change in the magnetic field in the opposite direction to the change in the magnetic field due to the gradient magnetic field pulse depending on time and space, and distorts the output gradient magnetic field waveform.
  • the output gradient magnetic field waveform of the gradient coil includes distortion such as response delay depending on the characteristic (Q value) of the gradient coil.
  • a circuit that performs feedback control or the like to compensate for this includes a gradient magnetic field. It is often included in a gradient magnetic field generation system including a power source. Such a control circuit also serves to compensate for eddy currents, but can also contribute to generating an output gradient magnetic field waveform different from the input gradient magnetic field waveform.
  • the distortion of the output gradient magnetic field waveform caused by the eddy current and the influence of the control circuit of the gradient magnetic field power supply varies depending on the shape of the input gradient magnetic field waveform, and causes various problems depending on the application axis of the gradient magnetic field.
  • distortion of the slice selective gradient magnetic field pulse causes an error in the excitation profile and the excitation position.
  • the influence is large in the case of VERSE (variable-rate-selective-excitation) method in which a high-frequency magnetic field pulse is applied while changing the application intensity of the gradient magnetic field pulse.
  • the distortion of the readout gradient magnetic field pulse causes artifacts such as distortion and ghost in the image.
  • EPI Echo Planar Imaging
  • spiral measurement that scans the k-space spirally
  • ultra-short TE measurement that starts scanning from the center of the k-space, etc.
  • the influence of distortion is significant.
  • One is to apply a magnetic field that compensates for eddy currents using shim coils or gradient magnetic field coils, etc., measuring the magnetic field generated by eddy currents after applying a gradient magnetic field pulse temporally and spatially, A waveform of a gradient magnetic field pulse that cancels the output characteristics of the eddy current obtained from the measurement result is output to compensate for distortion of the output gradient magnetic field waveform (Patent Documents 1 and 2).
  • Patent Document 3 a method of changing the irradiation timing of high-frequency magnetic field pulses according to the output gradient magnetic field waveform.
  • the distortion of the gradient magnetic field pulse is measured, the delay time is calculated from the center of gravity of the pulse area, and the timing of the high frequency magnetic field pulse and the gradient magnetic field pulse is changed based on the delay time, thereby reducing the excitation profile. Suppress.
  • a method has been proposed in which a change in the output gradient magnetic field waveform corresponding to the input gradient magnetic field waveform is eliminated by using a fixed input gradient magnetic field waveform (Non-Patent Document 1). In this method, when changing the excitation width of the high frequency magnetic field pulse, the amplitude of the high frequency magnetic field pulse is adjusted instead of changing the input gradient magnetic field waveform.
  • Non-Patent Document 2 A method for calculating the position and compensating for image distortion has been proposed (Non-Patent Document 2). In this method, the coefficient of the modeled mathematical formula is determined based on the appearance of the image. Also, instead of modeling the above method, a method has been proposed in which the coordinate position of the NMR signal in k-space is calculated from the result of actually measuring the gradient magnetic field waveform after the imaging is completed, and the distortion of the image is compensated (Non-Patent Document 3). ).
  • each of the above-described compensation techniques for the output gradient magnetic field waveform has the following problems.
  • Non-Patent Document 3 In the method of changing the irradiation timing of the high-frequency magnetic field pulse according to the output gradient magnetic field waveform (Patent Document 3), only the irradiation timing of the high-frequency magnetic field pulse is changed. Not compensated.
  • adjusting the amplitude of the high-frequency magnetic field pulse using only the specific input gradient magnetic field waveform is the SAR (Specific Absorption Rate). There is a risk of increasing the reference value indicating the accumulated amount.
  • Non-patent Document 2 if the system response does not fit the model formula by the RLC circuit, the distortion of the output gradient magnetic field waveform cannot be fully expressed, and the image distortion, ghost, etc. Will leave the artifact. This can occur when a complicated command such as feedback control is given by the gradient magnetic field power source. Further, since each coefficient of the model formula is determined based on the appearance of the image, there is a problem that the approximation accuracy of the model formula varies depending on the operator who determines the coefficient.
  • Non-patent Document 3 In the method of measuring the gradient magnetic field waveform after imaging and using it for reconstruction (Non-patent Document 3), fluctuations due to noise are superimposed because the measured gradient magnetic field waveform is used to calculate k-space coordinates. When processing such as a filter is performed to reduce noise, distortion of the gradient magnetic field waveform due to the filter processing becomes a problem.
  • an object of the present invention is to accurately calculate the distortion of the output gradient magnetic field waveform that changes in accordance with the input gradient magnetic field waveform and to perform various corrections using the calculated output gradient magnetic field waveform.
  • the present invention for solving the above problems calculates an input gradient magnetic field waveform and an output gradient magnetic field waveform corresponding to the input gradient magnetic field waveform, and uses the input gradient magnetic field waveform and the output gradient magnetic field waveform to generate the output gradient magnetic field waveform.
  • Various corrections are performed using the calculated value of the output gradient magnetic field waveform.
  • the high-frequency magnetic field pulse set in the imaging sequence is corrected and / or the k-space coordinates of the nuclear magnetic resonance signal obtained in the imaging sequence are corrected.
  • a response function is calculated by taking into account a plurality of elements that affect the output waveform of the gradient magnetic field, for example, the response functions of the control circuit of the eddy current and the gradient magnetic field power supply, and thereby the output of the gradient magnetic field. Since the waveform is calculated, not only eddy current but also the gradient magnetic field distortion caused by the control circuit can be compensated, and it could not be removed by the prior art without using a high-cost compensation gradient coil etc. Artifacts such as image distortion and ghosting can be suppressed.
  • FIG. 1 is a block diagram showing an overall configuration of an embodiment of an MRI apparatus according to the present invention.
  • This MRI apparatus uses an NMR phenomenon to obtain a tomographic image of a subject, and as shown in FIG. 1, a static magnetic field generation system 2, a gradient magnetic field generation system 3, a transmission system 5, and a reception system 6, a signal processing system 7, and a sequencer 4.
  • the static magnetic field generation system 2 includes a static magnetic field generation source such as a permanent magnet, a normal conducting magnet, or a superconducting magnet, and generates a uniform static magnetic field in a space where the subject 1 is placed.
  • a static magnetic field generation source such as a permanent magnet, a normal conducting magnet, or a superconducting magnet
  • a uniform static magnetic field in a space where the subject 1 is placed there are a vertical magnetic field method, a horizontal magnetic field method, etc. If the vertical magnetic field method is used, the body axis is used in the direction perpendicular to the body axis in the space around the subject 1 and if the horizontal magnetic field method is used. Generate a uniform static magnetic field in the direction.
  • the gradient magnetic field generation system 3 includes a gradient magnetic field coil 9 wound in the three axes of X, Y, and Z, which is a coordinate system (static coordinate system) of the MRI apparatus, and a gradient magnetic field power source 10 that drives each gradient magnetic field coil.
  • the gradient magnetic fields Gx, Gy, and Gz are applied in the three axis directions of X, Y, and Z by driving the gradient magnetic field power supply 10 of each coil in accordance with a command from the sequencer 4 to be described later.
  • a slice direction gradient magnetic field pulse is applied in a direction orthogonal to the slice plane (imaging cross section) to set a slice plane for the subject 1, and in the remaining two directions orthogonal to the slice plane and orthogonal to each other
  • a phase encoding direction gradient magnetic field pulse and a frequency encoding direction gradient magnetic field pulse are applied, and position information in each direction is encoded in the echo signal.
  • the sequencer 4 is a control means that repeatedly applies a high-frequency magnetic field pulse and a gradient magnetic field pulse in a predetermined pulse sequence, operates under the control of the signal processing system 7 (digital signal processing device 8), and tomographic image data of the subject 1 Various commands necessary for collection are sent to the transmission system 5, the gradient magnetic field generation system 3, and the reception system 6. There are various pulse sequences depending on the imaging method, and these are stored in advance on the magnetic disk 18 as a program.
  • the transmission system 5 irradiates the subject 1 with a high frequency magnetic field pulse in order to cause nuclear magnetic resonance to occur in the nuclear spins of the atoms constituting the biological tissue of the subject 1, and the high frequency oscillator 11, the modulator 12, and the high frequency It comprises an amplifier 13 and a high frequency coil (transmission coil) 14a on the transmission side.
  • the high-frequency magnetic field pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 at a timing according to a command from the sequencer 4, and the amplitude-modulated high-frequency magnetic field pulse is amplified by the high-frequency amplifier 13 and then close to the subject 1.
  • the subject 1 is irradiated with a high frequency magnetic field pulse.
  • the receiving system 6 detects an echo signal (NMR signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the subject 1, and receives a high-frequency coil (receiving coil) 14b on the receiving side and a signal amplifier 15 And a quadrature phase detector 16 and an A / D converter 17.
  • the quadrature phase detector 16 divides the signal into two orthogonal signals at the timing according to the command from the sequencer 4, and each signal is converted into a digital quantity by the A / D converter 17 and sent to the signal processing system 7.
  • the signal processing system 7 performs various data processing and display and storage of processing results.
  • the digital signal processing device 8, an external storage device such as a magnetic disk 18 and an optical disk 19, a display 20, and a digital signal processing device 8 has internal memories such as ROM21 and RAM22 for storing the processing results of 8 and data used by the digital signal processing device 8, and when the data from the receiving system 6 is input to the digital signal processing device 8, the digital signal processing The apparatus 8 performs processing such as signal processing and image reconstruction, and displays the tomographic image of the subject 1 as a result on the display 20 and records it on the magnetic disk 18 of the external storage device.
  • the digital signal processing device 8 performs calculation related to device characteristics using pre-measurement data, and calculation of parameters necessary for execution of the pulse sequence stored on the magnetic disk 18. Control the operation of sequencer 4.
  • the operation unit 25 inputs various control information of the MRI apparatus and control information of processing performed by the signal processing system 7 and includes a trackball or mouse 23 and a keyboard 24.
  • the operation unit 25 is disposed close to the display 20, and the operator controls various processes of the MRI apparatus interactively through the operation unit 25 while looking at the display 20.
  • the high-frequency coil 14a and the gradient magnetic field coil 9 on the transmission side face the subject 1 in the static magnetic field space of the static magnetic field generation system 2 into which the subject 1 is inserted, in the case of the vertical magnetic field method. If the horizontal magnetic field method is used, the subject 1 is installed so as to surround it.
  • the high-frequency coil 14b on the receiving side is installed so as to face or surround the subject 1.
  • the MRI apparatus of the present embodiment further includes, as a function of the signal processing system 7, an input waveform (input gradient magnetic field waveform) to the gradient magnetic field coil 9 and a gradient magnetic field waveform (output gradient magnetic field waveform) generated thereby in the static magnetic field space. And a gradient magnetic field waveform calculating means for calculating an output gradient magnetic field waveform using the response function.
  • the details of the digital signal processing apparatus 8 provided with these means are shown in FIG.
  • the digital signal processing device 8 includes a response function calculation unit 81 that calculates a response function of the gradient magnetic field, and an output gradient magnetic field waveform of the gradient magnetic field used in the actual imaging sequence using the calculated response function.
  • Output gradient magnetic field waveform calculating unit 82, gradient magnetic field waveform calculating unit 82, sequence calculating unit 83 for correcting high-frequency pulses and the like used in the imaging sequence using the output gradient magnetic field waveform calculated by gradient magnetic field waveform calculating unit 82, and gradient magnetic field waveform calculating unit 81 Using the calculated output gradient magnetic field waveform, the k-space coordinate calculation unit 84 corrects the coordinate position of the k-space where the NMR signal is arranged, and performs image reconstruction calculation using the NMR signal to generate image data.
  • a component 85 is provided.
  • the processing performed by the digital signal processing device 8 of the MRI apparatus of this embodiment is roughly divided into four processing.
  • the relationship between the four processes will be described with reference to the block diagram of FIG.
  • the first process 101 is a process performed by the response function calculation unit 81 as a pre-measurement, and is a process for calculating the response function 1112 of the control circuit of the eddy current and gradient magnetic field power supply.
  • Processes 102 to 104 are processes performed in actual imaging.
  • the process 102 is a process for estimating each output gradient magnetic field waveform using the response function 1112 obtained in the process 101 and the input gradient magnetic field waveforms 1113 and 1114 set for each imaging sequence. 82 does.
  • the process 103 is a process of recalculating the imaging sequence high-frequency magnetic field pulse using the result of the process 102, and is performed by the high-frequency magnetic field pulse calculating unit 831 of the sequence calculating unit 83.
  • the process 104 is a process for recalculating k-space coordinates corresponding to each sampling point of the measured NMR signal 1115 using the result of the process 102, and is performed by the k-space coordinate calculation unit 84.
  • each process will be described in detail.
  • step 301 the control parameters of the gradient magnetic field power supply are invalidated, and operations such as feedback control and feedforward control are set so that the gradient magnetic field power supply does not perform.
  • step 302 the output gradient magnetic field waveform is measured for each space.
  • the measurement of the output gradient magnetic field can be performed using a known technique (for example, the technique described in Patent Document 3). For example, as shown in FIG. 4, each space is measured at intervals of 50 mm in the X, Y, and Z axis directions.
  • the output gradient magnetic field waveform 1111 and the input gradient magnetic field waveform 1110 used at the time of measurement are stored in the memory RAM 22 (FIG. 1) of the signal processing system.
  • step 303 the input gradient magnetic field waveform and the measured output gradient magnetic field waveform stored in the RAM 22 are read, and the response function of the eddy current is calculated. A method for calculating the response function of the eddy current will be described later.
  • step 304 the control parameters of the gradient magnetic field power supply are validated to bring the feedback control, the feedforward control, and the like into operation. In this state, the output gradient magnetic field waveform is measured (step 305). This measurement result is also stored in the RAM 22.
  • step 306 the response function of the control circuit of the gradient magnetic field power supply is calculated. A method for calculating the response function of the control circuit of the gradient magnetic field power supply will also be described later. The above processing is performed for each gradient magnetic field coil arranged in the three-axis directions of X, Y, and Z.
  • the magnetic field generated by the eddy current distorts the output gradient magnetic field waveform 502 with respect to the input gradient magnetic field waveform 501 as shown in FIG.
  • the relationship between the input gradient magnetic field waveform and the output gradient magnetic field waveform can be expressed by the following equation, for example.
  • s is a Laplace variable
  • Gout (s) is a Laplace transform type of the output gradient magnetic field waveform
  • Gin (s) is a Laplace transform type of the input gradient magnetic field waveform
  • H (s) is an eddy current transfer characteristic.
  • the output gradient magnetic field waveform is individual data measured for each space as shown in FIG. Formula (1) is an example, and can be replaced with another expression such as the following formula.
  • the expression (1) is used. Since the input gradient magnetic field waveform Gin (s) and the output gradient magnetic field waveform Gout (s) are known, the transfer characteristic H (s) of eddy current can be theoretically calculated from the following equation. The response function of eddy current can be obtained by inverse Laplace transform of equation (3.1).
  • ILT [] represents the inverse Laplace transform.
  • the response function of eddy current is applied to the model formula, and the coefficient in the model formula is searched for and obtained.
  • An example of a model equation for the response function of eddy current is shown in the following equation (4).
  • the response function of Equation (4) expresses the transient response of the magnetic field generated by eddy current.
  • t is time
  • ⁇ 1 to ⁇ n are amplitude gains
  • ⁇ 1 to ⁇ n are time constants.
  • the number n of exponential functions in the equation (4) is determined by the type of time constant possessed by the parts constituting the MRI apparatus. According to the inventors' experience, if n is 10, a transient response of a magnetic field generated by an eddy current can be expressed sufficiently.
  • step 601 initial values of amplitude gain and time constant are set.
  • the initial value all the coefficients of the amplitude gain ⁇ are set to zero, and a positive numerical value other than zero is set to the time constant ⁇ .
  • step 603 only ⁇ i is changed, and the optimum value of ⁇ i is searched. At this time, the optimum value is a value that minimizes the evaluation amount e1 of the following equation.
  • Equation (4) is the response function shown in Equation (4).
  • Expression (5) is a numerical value indicating how similar the result of convolution of the response function with the output gradient magnetic field waveform and the input gradient magnetic field waveform is. Accordingly, there is essentially no change even if an evaluation formula other than the formula (5) is used.
  • -1 ⁇ i ⁇ a 1 range to calculate the e1 by changing the alpha i in steps 0.001, and optimum value alpha i which becomes a value e1 is minimized.
  • a search algorithm typified by the golden section method may be used.
  • step 604 After the searched optimum value ⁇ i is set to h (t), in step 604, only ⁇ i is changed, and the value of ⁇ i that minimizes the evaluation amount e1 is searched.
  • the search method is the same as ⁇ i .
  • step 605 it is determined whether the values of ⁇ i and ⁇ i have converged as a result of continuously searching for the optimum values of ⁇ i and ⁇ i and updating the values.
  • the convergence condition of ⁇ i and ⁇ i can be defined as follows, for example.
  • ⁇ i _prev and ⁇ i _prev are optimum values of ⁇ i and ⁇ i detected during the previous search. In other words, at least two repeated searches are required.
  • ⁇ i and ⁇ i are searched again. If the convergence condition is satisfied, if i is less than n, i is incremented, and the next combination of ⁇ i and ⁇ i is searched. If i is greater than or equal to n, the process ends.
  • ⁇ can take not only positive values but also negative values. It is correct from the principle of eddy current generation that ⁇ is only a positive value, but ⁇ 1 to ⁇ n are obtained in order, and the search is performed by ignoring the components of ⁇ and ⁇ to be obtained later. Can occur. The response function of eddy current is obtained by the above procedure.
  • the control circuit for the gradient magnetic field power supply causes the output gradient magnetic field waveform 702 to include behavior such as overshoot or undershoot with respect to the input gradient magnetic field waveform 701.
  • the gradient magnetic field pulse output from the gradient magnetic field power source is obtained by adding the component from the control circuit to the input waveform and can be expressed as the following equation.
  • gamp (t) is a gradient magnetic field pulse output from the gradient magnetic field power source
  • gin (t) is an input gradient magnetic field waveform
  • u (t) is a control component of the gradient magnetic field power source.
  • This control component is also defined by a model equation as in the case of the eddy current.
  • the model equation of the control component can be appropriately defined according to the control method adopted by the control circuit.
  • intelligent control such as PID control, simpler P control, PD control, PI control, or more complicated neural network can be assumed, and it is installed in the MRI apparatus to which the present invention is applied.
  • the control theory may be selected according to the specifications of the gradient magnetic field power supply. In the present embodiment, assuming that the control circuit of the gradient magnetic field power supply performs PID control, u (t) is defined by the following model equation.
  • K P proportional gain K I is an integral gain
  • K D is the derivative gain
  • Diff (t) in equation (9) is defined by the following equation.
  • delay is a time delay when calculating the input of PID control.
  • the relationship between gamp (t) and the output gradient magnetic field waveform is as follows.
  • gout (t) is the output gradient magnetic field waveform stored in the RAM 22 in step 305 of FIG.
  • h (t) is the response function of the eddy current shown in Equation (4), and is calculated by the above-described processing and is known. Therefore, when h (t) calculated in step 303 is convolved with gamp (t), which is the sum of the input gradient magnetic field waveform and the control component u (t) of the gradient magnetic field power supply, the most output gradient magnetic field waveform gout Find u (t) giving a gradient magnetic field waveform close to (t).
  • gamp (t) is the sum of the input gradient magnetic field waveform and the control component u (t) of the gradient magnetic field power supply
  • gout Find u (t) giving a gradient magnetic field waveform close to (t).
  • step 801 delay, K P , K I , and K D are all initialized to zero.
  • step 802 only the delay is changed to search for the optimum delay.
  • the optimum value is a value that minimizes the evaluation amount e2 of the following equation.
  • the search for delay at which the evaluation amount e2 is minimized is performed in the same manner as the search processing for the coefficients ⁇ 1 to ⁇ n and ⁇ 1 to ⁇ n of the response function described above. That is, an all-point search may be performed within a specific range, or a search algorithm typified by the golden section method may be used.
  • Diff (t) is obtained.
  • step 803 the optimum K P is searched.
  • the condition under which K P is optimal is a value that minimizes the evaluation quantity e2 in equation (12).
  • the K I in step 804 searches the K D in step 805.
  • the response function of the control circuit for the gradient magnetic field power source finally obtained is expressed by the following equation.
  • the response function h (t) of the eddy current and the response function of the control circuit (control component u (t)) are calculated.
  • the magnetic field waveform can be calculated.
  • the parameters of the response function h (t) (formula (4)) and u (t) for solving Equation (11) are obtained for individual data measured for each space as shown in FIG. It is stored in the RAM 22 of the signal processing system 7 and used for pulse sequence calculation and image reconstruction in actual imaging.
  • FIG. 10 shows an example of the estimation result of the output gradient magnetic field waveform when the control component u (t) of the gradient magnetic field power supply is not added.
  • Output gradient magnetic field waveform calculation process 102 >> Next, returning to the block diagram of FIG. 2, the processing 102 performed by the output gradient magnetic field waveform calculation unit 82 will be described.
  • an imaging pulse sequence is calculated by the sequence calculation unit 83 according to the imaging purpose, and a high-frequency pulse (amplitude, envelope), slice selection gradient magnetic field pulse (input gradient magnetic field waveform) included in the calculated imaging pulse sequence is calculated.
  • a read gradient magnetic field pulse input gradient magnetic field waveform 1114) and the like are stored in the RAM 22.
  • the output gradient magnetic field waveform of the imaging pulse sequence is calculated using the response function 1112 calculated by the response function calculation means 81 and the set input gradient magnetic field waveforms 1113 and 1114.
  • FIG. 11 shows a flowchart of the calculation process of the output gradient magnetic field waveform.
  • step 1101 the visual field center position in the three-axis directions of X, Y, and Z is acquired from the coordinate position of the imaging visual field.
  • the visual field center position is calculated from the following equation.
  • Axis represents each axis of X, Y, and Z (hereinafter the same)
  • FOV (Axis, min) and FOV (Axis, max) are the minimum of the imaging field of view of each axis specified according to the imaging conditions Coordinate position and maximum coordinate position.
  • step 1102 from the response function h (t) for each space stored in the RAM 22 and the control component u (t) of the gradient magnetic field power supply, the response function h (t ) And control component u (t) are selected and acquired.
  • the acquisition of the response function h (t) and the control component u (t) of the gradient magnetic field power supply may be obtained from the result calculated at the position closest to the visual field center position, or the result calculated for each space May be interpolated to obtain the response function h (t) and the control component u (t) of the gradient magnetic field power supply.
  • interpolating follow the formula below.
  • h (Axis, t) and u (Axis, t) are the response function at the visual field center position of each axis and the control component of the gradient magnetic field power supply circuit.
  • P1 and P2 represent the calculation positions of the response function and the control component of the gradient magnetic field power supply circuit that are closest in the front-rear direction of the visual field center position.
  • step 1103 an input gradient magnetic field waveform with respect to the X, Y, and Z directions is acquired.
  • the input gradient magnetic field waveform is set according to the imaging pulse sequence and is stored in the RAM 22, and is the slice gradient magnetic field pulse 1113 or the read gradient magnetic field pulse 1114 shown in FIG. In either case, the processing contents are the same.
  • step 1104 a waveform obtained by synthesizing the three-axis input gradient magnetic field waveform represented by the following equation is calculated.
  • gin (Axis, t) is the acquired input gradient magnetic field waveform of each axis.
  • step 1105 the output gradient magnetic field waveform is estimated using the response function h (t) and the control component u (t) of the gradient magnetic field power supply for the input gradient magnetic field waveform for each axis.
  • the estimation of the output gradient magnetic field waveform is expressed by the following equation as a form in which equation (11) is expanded to three axes.
  • gest (Axis, t) is an estimation result of the output gradient magnetic field waveform of each axis obtained by calculation.
  • the waveform of each axis is synthesized according to the following equation.
  • the output gradient magnetic field waveforms of the slice gradient magnetic field and readout gradient magnetic field estimated as described above are recalculated for the high-frequency magnetic field pulse (Fig. 2: 103) and the correction of the k-space coordinate position of the NMR signal (Fig. 2: Used for processing 104).
  • High-frequency magnetic field pulse calculation process 103 the high-frequency magnetic field pulse is recalculated using the output gradient magnetic field waveform gest (t) and the input gradient magnetic field waveform gin (t) of the slice gradient magnetic field obtained by the output gradient magnetic field waveform calculation unit 82.
  • Recalculation of the high frequency magnetic field pulse calculates the envelope and frequency of the high frequency pulse according to the following formula.
  • RF env '(t) is the envelope of the high frequency magnetic field pulse after recalculation
  • RF env (t) is the envelope of the high frequency magnetic field pulse before recalculation
  • RF freq ' (t) is the high frequency magnetic field pulse after recalculation
  • RF freq (t) is the frequency of the high-frequency magnetic field pulse before recalculation.
  • the amplitude, frequency, and irradiation time of the envelope of the high-frequency magnetic field pulse are changed in proportion to the intensity of the output gradient magnetic field waveform of the slice gradient magnetic field. Based on these equations, it is possible to make changes in consideration of device restrictions and imaging conditions.
  • phase of the high-frequency magnetic field pulse may be changed with time, and as a result, frequency modulation may be realized.
  • the calculation of the phase for realizing the frequency modulation follows the following equation.
  • g_standard is the gradient magnetic field strength [T / m] that is the reference for frequency calculation of the high-frequency magnetic field pulse
  • distance is the distance from the magnetic field center to the excitation position [m]
  • ⁇ is the gyromagnetic ratio [Hz / T]
  • ⁇ t Is the phase change interval time [s] of the high-frequency magnetic field pulse.
  • the range in which the phase can be changed is 0 to 360 [deg]
  • the range of distance is also limited. Therefore, in order to widen the range in which the distance can be set, it is preferable that the phase change is small. Therefore, in order to minimize the phase change, the g_standard and the frequency of the high frequency magnetic field pulse are calculated according to the following formula.
  • Average [] is the mean value calculation process within the high frequency magnetic field pulse application section
  • RF freq ' is the frequency of the high frequency magnetic field pulse after recalculation without time change
  • RF freq is before recalculation without time change The frequency of the high frequency magnetic field pulse.
  • Another example that needs to be changed is the suppression of the irradiation time of the high-frequency magnetic field pulse.
  • the irradiation time of the high-frequency magnetic field pulse changes, but depending on the pulse sequence, the irradiation time may change, especially the irradiation time may be extended. It may not be preferable.
  • the recalculation procedure of the high frequency magnetic field pulse that does not change the irradiation time will be described with reference to the flowchart shown in FIG.
  • step 1021 and 1022 the output gradient magnetic field waveform for each axis is calculated, and the output gradient magnetic field waveform for each axis is synthesized. These processes are performed by the output gradient magnetic field waveform calculation unit 82, and are the same as the processes performed in step 1103 and step 1104 in FIG.
  • step 1031 a high frequency magnetic field pulse is calculated according to equations (19) to (21).
  • step 1032 the irradiation time increase rate ⁇ RF is calculated.
  • the increase rate ⁇ RF is calculated according to the following equation.
  • gin (t) is the input gradient magnetic field waveform
  • gest (t) is the estimation result of the output gradient magnetic field waveform obtained from the equation (18)
  • StartTime is the high frequency magnetic field pulse irradiation start time
  • EndTime is the high frequency magnetic field pulse irradiation end time Represents.
  • the increase rate beta RF determines whether exceeds 1.0. If an increase rate beta RF exceeds 1.0, so that the irradiation time of the high-frequency magnetic field pulse after the calculation is extended. In that case, the process moves to step 1034 to recreate the input gradient magnetic field waveform of the slice gradient magnetic field.
  • the reconstruction of the input gradient magnetic field waveform follows the following formula. By multiplying the input gradient magnetic field waveform by the irradiation time increase rate ⁇ RF , the intensity of the input gradient magnetic field waveform is increased, thereby shortening the application time of the high-frequency magnetic field pulse. At this time, since only the entire input gradient magnetic field waveform is increased uniformly, no special shim coil or gradient magnetic field coil is required.
  • step 1021 of processing 102 After re-creating the input gradient magnetic field waveform, return the result to step 1021 of processing 102, calculate the output gradient magnetic field waveform (steps 1021 and 1022), recalculate the high-frequency magnetic field pulse (step 1031), and increase the irradiation time
  • the rate ⁇ RF is calculated (step 1032). If the calculated increase rate ⁇ RF is 1.0 or less, the process is terminated, and if it is not 1.0 or less, the input gradient magnetic field waveform is created again. At this time, gin (t) in Equation (26) is the input gradient magnetic field waveform recreated immediately before.
  • the irradiation time of the high-frequency magnetic field pulse can be finally prevented from extending compared to the irradiation time of the high-frequency magnetic field pulse before recalculation. it can. This eliminates the need to change the pulse sequence other than the slice gradient magnetic field pulse when using the recalculated high-frequency magnetic field pulse.
  • FIG. 13 shows an image 1301 obtained using the high-frequency magnetic field pulse before recalculation and an image 1302 obtained using the high-frequency magnetic field pulse after recalculation. Both images are the imaging results of ultra-short TE measurement using the VERSE method. Compared to the image 1301 obtained using the high-frequency magnetic field pulse before the recalculation, the image 1302 obtained using the high-frequency magnetic field pulse after the recalculation is a background signal (a signal from other than the imaging section) indicated by the arrow. It can be seen that a good excitation profile is obtained because of less contamination.
  • ⁇ k-space coordinate recalculation process 104 the k-space coordinate recalculation process 104 performed by the k-space coordinate calculation unit 84 in FIG. 2 will be described.
  • the NMR signal is arranged using the output gradient magnetic field waveform gest (t) and the input gradient magnetic field waveform gin (t) of the readout gradient magnetic field obtained by the output gradient magnetic field waveform calculation unit 82.
  • FIG. 14 shows an example of a flowchart showing a recalculation procedure of k-space coordinates. This processing procedure is also the same processing flow as the high-frequency magnetic field pulse recalculation flowchart shown in FIG. 12, and the input gradient magnetic field waveform of the readout gradient magnetic field is recreated based on the recalculation result of the k-space coordinates.
  • the k-space coordinates are calculated using the input gradient magnetic field waveform gin (t) and the output gradient magnetic field waveform gest (t) of the read gradient magnetic field obtained in the process 102.
  • the recalculation of k-space coordinates using the output gradient magnetic field waveform gest (t) follows the following formula.
  • ⁇ is the magnetic rotation ratio [Hz / T].
  • This equation is obtained by replacing gin (j) with gest (j) in the equation for calculating k-space coordinates using the general input gradient magnetic field waveform gin (t) (the following equation (28)).
  • step 1042 the reduction rate ⁇ k of the k space is calculated.
  • the reduction rate ⁇ k is calculated according to the following equation.
  • EndTime represents the read gradient magnetic field pulse application end time (or NMR signal measurement end time).
  • step 1043 it is determined whether or not the reduction rate ⁇ k is less than 1.0.
  • the increase rate ⁇ k is less than 1.0
  • the k-space coordinate k ′ (t) after the recalculation is reduced with respect to the k-space coordinate k (t) before the recalculation.
  • the process proceeds to step 1044, where the input gradient magnetic field waveform of the read gradient magnetic field is recreated.
  • the reconstruction of the input gradient magnetic field waveform follows the following formula.
  • the intensity of the input gradient magnetic field waveform is increased, thereby expanding the range of the k-space coordinates.
  • the input gradient magnetic field waveform can only be increased uniformly, so that no special shim coil or gradient magnetic field coil is required.
  • FIG. 15 shows an image 1501 obtained using k-space coordinates before recalculation and an image 1502 obtained using k-space coordinates after recalculation. Both images are the imaging results of EPI measurement, and the readout gradient magnetic field is applied in the horizontal direction in the image. Compared to the image 1501 obtained using the k-space coordinates before recalculation, the image 1502 obtained using the k-space coordinates after recalculation has reduced distortion of the image indicated by the arrow. It can be seen that the compensation for the readout gradient magnetic field is properly compensated.
  • Imaging is performed using the high-frequency magnetic field pulse corrected in the process 103, and the acquired NMR signal is reconstructed using the k-space coordinates corrected in the process 104.
  • an external signal is not mixed into the selected slice, and an image in which image distortion and artifacts are suppressed can be obtained.
  • Non-Patent Document 1 In the method using a fixed input gradient magnetic field waveform (Non-Patent Document 1), it can be applied only to a specific pulse sequence, whereas in the present embodiment, the output gradient magnetic field waveform for any input gradient magnetic field waveform can be estimated. There is no restriction that a fixed input gradient magnetic field waveform must be used.
  • each coefficient of the model formula is obtained as a criterion for judging the appearance of the image, but in this embodiment, based on the output gradient magnetic field waveform. Since the coefficient of the model formula is calculated, it is not necessary to rely on human judgment, and a highly accurate model formula can be calculated.
  • the configuration of the digital signal processing unit 8 shown in FIG. 2 and the four characteristic processes 101 to 104 are performed in the digital signal processing unit 8 are the same as in the first embodiment.
  • the present embodiment is characterized in that the process 101 for calculating the response function of the control circuit for the eddy current and gradient magnetic field power supply is obtained by a simpler procedure than the first embodiment. That is, in the first embodiment, the measurement of the output gradient magnetic field waveform was performed twice each in the state in which the control of the gradient magnetic field power supply circuit was invalidated and in the valid state. Only measure the magnetic field waveform.
  • the processing flow of this embodiment will be described with reference to the flowchart of FIG.
  • step 1601 the output gradient magnetic field waveform is measured for each space.
  • the control parameters of the gradient magnetic field power supply are not changed, and the apparatus is in a normal state. That is, the control parameter of the gradient magnetic field power source is in a valid state.
  • the measurement result of the output gradient magnetic field waveform is stored in the RAM 22. This process is the same as the process performed in steps 302 and 305 of FIG. 3 showing the flow of the first embodiment.
  • step 1602 the input gradient magnetic field waveform and the measured output gradient magnetic field waveform stored in the RAM 22 are read, and the response function of the control circuit of the eddy current and gradient magnetic field power supply is calculated.
  • step 1701 an eddy current response function is calculated.
  • This process is the same as the process of step 303 in the first embodiment described above. That is, a response function h (t) that minimizes the evaluation amount e1 of the following equation (31) is searched.
  • gout (t) is the output gradient magnetic field waveform stored in the RAM 22 in step 1601
  • gin (t) is the input gradient magnetic field waveform.
  • h (t) is an eddy current response function defined by equation (4), and the coefficients ⁇ 1 to ⁇ n and ⁇ 1 to ⁇ n of equation (4) are searched in the same manner as in step 303 in the first embodiment.
  • Equation (31) has the same shape as Equation (5) of the first embodiment, but the output gradient magnetic field waveform gout (t) of Equation (5) is measured with the control circuit of the gradient magnetic field power supply disabled.
  • the output gradient magnetic field waveform gout (t) of Expression (31) includes distortion due to the control circuit of the gradient magnetic field power supply. Therefore, the calculated response function h (t) of the eddy current includes not only the eddy current but also the characteristics of the control circuit of the gradient magnetic field power supply.
  • step 1702 a response function of the control circuit of the gradient magnetic field power supply is calculated.
  • This process is the same as the process in step 306 of the first embodiment. That is, a search is made for gamp (t) that minimizes the evaluation amount e2 of the equation (32).
  • Gamp (t) in Expression (32) is defined by Expression (33).
  • u (t) is a control component of the gradient magnetic field power source defined by equation (9), but the eddy current response function h (t) obtained in step 1701 includes the control of the gradient magnetic field power source. Since the response component of the circuit is also included, the component that could not be expressed by the model of the eddy current response function (equation (4)) is applied to equation (9).
  • the response function h (t) of the eddy current is recalculated in step 1703. This is by defining the response function u (t) of the control circuit of the gradient magnetic field power supply, so that each coefficient ( ⁇ 1 to ⁇ n and ⁇ 1 to ⁇ n) of the current response function h (t) (equation (4)) is optimal. This is because the value changes. Recalculate the response function h (t) by replacing gin (t) in equation (31) with gamp (t) obtained in step 1702 and each coefficient of h (t) ⁇ ⁇ ⁇ that minimizes the evaluation quantity e1 By searching.
  • a response function taking into account the eddy current response function and the control component of the gradient magnetic field power source is obtained. The response function is obtained for each piece of data measured for each space and stored in the RAM 22 of the signal processing system 7.
  • the procedure of the process 102 for estimating the output gradient magnetic field waveform and the processes 103 and 104 for recalculating the high-frequency magnetic field pulse and k-space coordinates using the response function stored in the RAM 22 is the same as that of the first embodiment. The description is omitted.
  • the number of times of calculating the response function is increased, the number of times of measurement of the output gradient magnetic field waveform can be reduced, and the work of changing the control parameter of the gradient magnetic field power source can be omitted, with a simpler procedure.
  • the response function required to recalculate the excitation profile can be calculated.
  • Other effects are the same as in the first embodiment.
  • the configuration of the digital signal processing unit 8 shown in FIG. 2 and the four processes 101 to 104 performed by the digital signal processing unit 8 are the same as in the first embodiment.
  • the present embodiment is characterized in that the response function of the control circuit of the eddy current and the gradient magnetic field power supply is obtained without being applied to a specific model equation.
  • the processing flow of this embodiment is shown in FIG. Steps 1801 and 1802 in FIG. 18 are the same as steps 1601 and 1602 in FIG. 16 described in the second embodiment.
  • step 1801 a plurality of types of input gradient magnetic field waveforms are processed in step 1801. Processing to measure the output gradient magnetic field waveform, processing to directly derive the response function of the eddy current and the gradient power supply control circuit without using a model formula in step 1802, processing to correct the derived response function in step 1803 It is.
  • the eddy current transfer characteristic H ( s) can be calculated from equation (3.1). Reprint equation (3.1). However, when the value of Gin (s) becomes zero, the value of H (s) diverges and a correct value cannot be obtained. In order to prevent this, in the first and second embodiments, the response function that gives the optimum evaluation value is derived by applying the response function of the eddy current to the model formula.
  • the response function h (t) calculated from the equation (3.1) or the equation (3.2) obtained by inverse Laplace transform is used without using the model equation.
  • a plurality of types of input gradient magnetic field waveforms used for calculating the response function of the eddy current and the gradient magnetic field power supply control circuit are prepared.
  • the output gradient magnetic field waveform is measured using a plurality of types of input gradient magnetic field waveforms, and stored in the RAM 22.
  • An input gradient magnetic field waveform having a steep shape as much as possible is preferable. This is because a gently changing waveform does not include a high frequency range, and a high frequency response cannot be calculated.
  • it is desirable to use an impulse wave but because it is difficult to output an impulse wave due to restrictions such as the output performance of the gradient magnetic field power supply, in this embodiment, a plurality of trapezoidal waves as shown in FIG. 19 are used. Use.
  • step 1802 the transfer characteristic H (s) is calculated by equation (3.1).
  • the zero point appearing in the spectrum (Laplace transform type) of one input gradient magnetic field waveform (trapezoidal wave) is replaced with the spectrum data of another input gradient magnetic field waveform (trapezoidal wave) to avoid the zero point.
  • Gin1 (s) is a Laplace transform type of trapezoidal wave 1901
  • Gout1 (s) is a Laplace transform type of output gradient magnetic field waveform when trapezoidal wave 1901 is used as an input gradient magnetic field waveform
  • Gin2 (s) is Laplace transform type of trapezoidal wave 1902
  • Gout2 (s) is a Laplace transform type of output gradient magnetic field waveform when trapezoidal wave 1902 is used as the input gradient magnetic field waveform.
  • FIG. 20 shows an outline of the amplitude spectrum obtained by Laplace transforming the trapezoidal wave group shown in FIG.
  • the positions where the zeros appear in the spectrums 2001 and 2002 of the trapezoidal waves 1901 and 1902 are different between the spectrum 2001 and the spectrum 2002.
  • the zero point can be avoided.
  • two trapezoidal waves are shown in FIG. 19, three or more trapezoidal waves may be prepared.
  • Equation (34) among the plurality of input gradient magnetic field waveforms, the combination of the input gradient magnetic field waveform having the largest spectrum value (not zero) and the output gradient magnetic field waveform resulting therefrom is selected. What is necessary is just to use for calculation of a transfer characteristic.
  • step 1803 a process for eliminating the influence of noise from the response function h (t) calculated in step 1802 is performed. This is because Gout1 (s) and Gout2 (s) are actually measured data, and noise is superimposed, and there are cases where correct transfer characteristics cannot be obtained. In the present embodiment, such noise is removed by performing two-stage filter processing in the Laplace space.
  • the transfer characteristic obtained by the equation (34) is multiplied by the optimum filter of the following equation to suppress the influence of noise.
  • H ′ (s) is the transfer characteristic after applying the optimum filter
  • ⁇ (s) is the optimum filter.
  • ⁇ (s) a Wiener filter that is a known technique is used as the optimum filter ⁇ (s).
  • the definition formula of the Wiener filter is as follows.
  • N (s) is a noise component of Gout1 (s)
  • S (s) is a signal component excluding the noise of Gout1 (s).
  • N (s) and S (s) are calculated according to a known technique, where N (s) is the highest frequency signal value of the power spectrum of Gout1 (s), and N (s ) Is subtracted S (s).
  • a noise removal filter is executed on the transfer characteristic after the optimum filter is applied.
  • a median filter is used as a noise removal filter.
  • H ′′ (s) is the transfer characteristic after the median filter is applied.
  • the response function that combines the eddy current and the gradient power supply control circuit can be obtained by inverse Laplace transform of Equation (37).
  • ILT [] represents the inverse Laplace transform.
  • the procedure for estimating the output gradient magnetic field waveform and recalculating the high-frequency magnetic field pulse / k-space coordinates using the response function obtained as described above is the same as in the first embodiment.
  • the response function can be obtained without applying a specific model formula, and various hardware configurations ( It can be flexibly applied to a gradient magnetic field coil or the like) or a software configuration (control method of gradient magnetic field power supply).
  • the present invention is not limited to the above-described embodiment, and various modifications can be made. It is.
  • the model formula that expresses the distortion of the output gradient magnetic field waveform due to the eddy current and the distortion of the output gradient magnetic field waveform by the control circuit of the gradient magnetic field power supply is not limited to that exemplified in the embodiment, and other model formulas may be used. Is possible.
  • the control circuit of the eddy current and the gradient magnetic field power source is taken up, and the case where the response function by them is used for the calculation of the output gradient magnetic field waveform has been described.
  • Other components such as residual magnetic fields with hysteresis found in MRI apparatus using permanent magnets, for example, can be added to the calculation of the output gradient magnetic field waveform.
  • the processing procedure of the flowchart shown for description of each embodiment is an example, there is a process which can be omitted, and it is also possible to add another process as needed.
  • the present invention is used when a high-frequency magnetic field pulse is calculated from an output gradient magnetic field waveform calculated for a slice gradient magnetic field, or when k-space coordinates where an NMR signal is arranged are corrected by an output gradient magnetic field waveform calculated for a read gradient magnetic field.
  • the readout gradient magnetic field also includes a two-direction gradient magnetic field in a pulse sequence such as a radial scan or a spiral scan that does not distinguish between the readout gradient magnetic field and the phase encoding gradient magnetic field.
  • the outline of the present invention clarified by the description of each of the above embodiments is as follows. That is,
  • the MRI apparatus of the present invention is An imaging means comprising a static magnetic field generating means, a gradient magnetic field generating means, a high frequency magnetic field generating means and a high frequency magnetic field detecting means, a control means for controlling the operation of the imaging means based on an imaging sequence, and the high frequency magnetic field detecting means detected
  • the computing means uses an input waveform input to the gradient magnetic field generating means and an output waveform generated by the gradient magnetic field generating means according to the input waveform, and a response function of a plurality of elements that affect the output waveform.
  • Gradient magnetic field calculation for calculating an output waveform from an input waveform of a gradient magnetic field pulse set in the imaging sequence, using a response function calculating unit that calculates a response function that is a sum, and a response function calculated by the response function calculating unit And k-space coordinates of the nuclear magnetic resonance signal obtained by the correction of the high-frequency magnetic field pulse set in the imaging sequence and / or the nuclear magnetic resonance signal obtained by the imaging sequence using the calculated value of the output waveform calculated by the gradient magnetic field calculating unit Compensating means for correcting the above is provided.
  • the response function calculating means calculates a response function for each of a plurality of elements that change the output waveform, and the gradient magnetic field calculating means uses a response function that is a sum of the response functions of the plurality of elements to output the output waveform. Perform the calculation.
  • the response function calculating means includes a model setting means for setting a response function model for each element, a difference between a calculated value of the output waveform calculated by the response function model set by the model setting means, and an actually measured value of the output waveform.
  • Model optimization means for optimizing the response function model so as to minimize the response function model, and the response function model optimized by the model optimization means is used as the response function of each element.
  • the response function calculating means uses a plurality of input values as the input values, and uses a combination of a plurality of input values and corresponding measured values, directly from the relationship between the input values and the measured values of the output waveform, Calculate the response function.
  • the response function calculation unit includes a filter unit that filters the response function after calculation.
  • the filter means includes a Wiener filter and a median filter.
  • the gradient magnetic field generating means includes a control circuit that controls an input waveform based on an output waveform of the generated gradient magnetic field,
  • the plurality of elements that change the output waveform include distortion of the output waveform accompanying control by the control circuit.
  • the response function calculating unit calculates a response function for each of a plurality of gradient magnetic fields having different positions with respect to the center of the static magnetic field generated by the static magnetic field generating unit.
  • the gradient magnetic field pulse calculated by the gradient magnetic field calculation means is a slice selection gradient magnetic field in the imaging sequence, and the compensation means uses a high-frequency pulse set in the imaging sequence based on a calculated value of the slice selection gradient magnetic field. High frequency pulse correction means for correcting is provided.
  • the high-frequency pulse correction means includes irradiation time calculation means for calculating a change in irradiation time of the high-frequency pulse before and after correction, and when the change in irradiation time calculated by the irradiation time calculation means exceeds a predetermined increase rate, The input value of the gradient magnetic field calculating means is changed, and the high-frequency pulse is corrected using the estimated value of the output waveform estimated by the gradient magnetic field calculating means from the changed input value.
  • the gradient magnetic field pulse calculated by the gradient magnetic field calculation means is a readout gradient magnetic field in the imaging sequence, and the compensation means determines that the nuclear magnetic resonance signal obtained in the imaging sequence is based on the calculated value of the readout selective gradient magnetic field.
  • K-space coordinate correcting means for correcting the arranged k-space coordinates is provided.
  • the correction means includes k-space size calculation means for calculating a change in the size of k-space coordinates before and after correction, and the change in the size of the k-space calculated by the k-space size calculation means exceeds a predetermined reduction rate.
  • the input value of the gradient magnetic field calculating means is changed, and the k-space coordinates are corrected using the recalculated value of the output waveform calculated from the changed input value by the gradient magnetic field calculating means.
  • the method of estimating the output gradient magnetic field waveform generated by the gradient magnetic field generating means of the magnetic resonance imaging apparatus of the present invention is as follows: Inputting an input value of the gradient magnetic field generation means and an actual measurement value of an output gradient magnetic field waveform generated by the gradient magnetic field generation means; Setting a response function model for each of a plurality of factors that affect the output gradient magnetic field waveform; Using each input value and actual measurement value of the input gradient magnetic field waveform to optimize each response function model; Calculating an estimated value of an output gradient magnetic field waveform using the optimized response function model for an input value of a desired gradient magnetic field; It is characterized by including.
  • the plurality of factors include an eddy current generated by a gradient magnetic field and an output of a control circuit provided in the gradient magnetic field generating means.
  • the actual measurement value includes a first actual measurement value when the control circuit of the gradient magnetic field generation unit is invalidated, and a second actual measurement value when the control circuit is enabled,
  • the response function of the output gradient magnetic field waveform by the control circuit is obtained using the first actual measurement value and the second actual measurement value.
  • 1 subject 2 static magnetic field generation system, 3 gradient magnetic field generation system, 4 sequencer, 5 transmission system, 6 reception system, 7 signal processing system, 8 digital signal processing unit (CPU), 9 gradient magnetic field coil, 10 gradient magnetic field power supply , 11 High-frequency transmitter, 12 modulator, 13 high-frequency amplifier, 14a high-frequency coil (transmitting coil), 14b high-frequency coil (receiving coil), 15 signal amplifier, 16 quadrature detector, 17 A / D converter, 18 magnetic disk , 19 optical disk, 20 display, 21 ROM, 22 RAM, 23 trackball or mouse, 24 keyboard, 25 operation unit, 81 response function calculation unit, 82 gradient magnetic field waveform calculation unit, 83 sequence calculation unit, 84 k-space coordinate calculation unit 85 Image reconstruction unit.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Nonlinear Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

　種々の要因によって生じる出力傾斜磁場波形の歪みを高精度に補償することができ、画像の歪みやゴースト等のアーチファクトのない信頼性の高い画像を得ることを目的とするためには、入力傾斜磁場波形と当該入力傾斜磁場波形に対応する出力傾斜磁場波形を算出し、入力傾斜磁場波形と出力傾斜磁場波形とを用いて、該出力傾斜磁場波形に影響を与える複数の要素の応答関数の合算である応答関数を算出し、応答関数を用いて、撮像シーケンスに設定された傾斜磁場パルスの入力傾斜磁場波形から出力傾斜磁場波形を算出し、算出した出力傾斜磁場波形の計算値を用いて、各種補正を行う。

Description

磁気共鳴イメージング装置および傾斜磁場波形推定方法
　本発明は、被検体中の水素や燐等からの核磁気共鳴信号(以下、NMR信号と呼ぶ)を測定し、核の密度分布や緩和時間分布等を画像化する磁気共鳴イメージング(以下、MRIと呼ぶ)装置に関し、特に、傾斜磁場の歪みに起因する画質劣化を補償する技術に関する。
　MRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生するNMR信号を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮影においては、被検体を静磁場(分極磁場Ｂ0)内に配置した上で、特定の領域を選択励起するためにスライス選択傾斜磁場パルスと共に高周波磁場パルスを印加し、その後、位相エンコード傾斜磁場パルスや読出し傾斜磁場パルスを印加することにより、励起範囲内をエンコードし、位置情報を付与する。
　傾斜磁場パルスは、直交する三軸方向に線形の傾斜磁場を発生する複数のコイルに、傾斜磁場電源からパルス状に変化する電流(入力傾斜磁場波形という)を供給することにより、所望の傾斜磁場パルス(出力傾斜磁場波形という)を出力するようにしている。理想的には、入力傾斜磁場波形とそれによってMRI装置に発生する傾斜磁場波形は一致しなければならないが、種々の要因により誤差が生じる。
　誤差を生じる原因の一つは、傾斜磁場パルスによる急激な磁界の変化により、MRI装置を構成する磁性体部品あるいは電気回路に発生する渦電流による影響である。渦電流は、傾斜磁場パルスによる磁界の変化とは逆方向の磁界の変化を時間及び空間に依存して与え、出力傾斜磁場波形を歪ませる。
　また一般に、傾斜磁場コイルの出力傾斜磁場波形は、傾斜磁場コイルの特性(Ｑ値)に依存する応答遅れ等の歪みを含むが、これを補償するためにフィードバック制御などを行う回路が、傾斜磁場電源を含む傾斜磁場発生系には含まれている場合が多い。このような制御回路は、渦電流の補償の役割も果たすが、入力傾斜磁場波形と異なる出力傾斜磁場波形を生成する一因にもなり得る。
　これら渦電流及び傾斜磁場電源の制御回路の影響により生じる出力傾斜磁場波形の歪みは、入力傾斜磁場波形の形状に依存して変化し、その傾斜磁場の印加軸に応じて種々の問題を生じる。
　例えば、スライス選択傾斜磁場パルスの歪みは、励起プロファイル及び励起位置に誤差を生じさせる。特に傾斜磁場パルスの印加強度を変化させながら高周波磁場パルスを印加するVERSE(variable rate selective excitation)法などの場合には影響が大きい。
　読み出し傾斜磁場パルスの歪みは、画像に歪みやゴースト等のアーチファクトを生じさせる。特に、読み出し傾斜磁場パルスの印加極性を反転させながら計測を行うEPI(Echo Planar Imaging)計測や、k空間をらせん状に走査するスパイラル計測、k空間中心から走査を開始する超短TE計測などの場合に歪みの影響が顕著である。
　上記問題を解決するためには、入力傾斜磁場波形に依存して変化する出力傾斜磁場波形の歪みを補償する必要がある。既存の補償技術として、以下に示すような手法が提案されている。
　一つは、シムコイルや傾斜磁場コイルなど用いて渦電流を補償する磁場を印加するというものであり、傾斜磁場パルスを印加した後に発生する渦電流により生じる磁場を時間的・空間的に計測し、その計測結果から得られる渦電流の出力特性を打ち消す傾斜磁場パルスの波形を出力して出力傾斜磁場波形の歪みを補償する(特許文献1、特許文献2)。
　またスライス選択傾斜磁場パルスの歪みの補償に特化した手法として、出力傾斜磁場波形に応じて高周波磁場パルスの照射タイミングを変更する方法が提案されている(特許文献3)。この方法では、傾斜磁場パルスの歪みを測定し、パルス面積の重心から遅延時間を算出し、遅延時間をもとに高周波磁場パルスと傾斜磁場パルスのタイミングを変更することで、励起プロファイルの劣化を抑制する。また、固定の入力傾斜磁場波形を用いることで、入力傾斜磁場波形に応じた出力傾斜磁場波形の変化を排除する方法も提案されている(非特許文献1)。この方法では、高周波磁場パルスの励起幅を変更する場合は、入力傾斜磁場波形を変化させるのではなく、高周波磁場パルスの振幅を調整している。
　読出し傾斜磁場パルスの歪みの補償に特化した手法としては、システム応答をRLC回路でモデル化し、入力傾斜磁場波形に対する出力傾斜磁場波形を推定し、その結果を用いてNMR信号のk空間における座標位置を算出し、画像の歪みを補償する方法が提案されている(非特許文献2)。この方法では、モデル化した数式の係数は、画像の見た目を基準にして決定している。また上記手法のモデル化の代わりに、撮像終了後に傾斜磁場波形を実測した結果からNMR信号のk空間における座標位置を算出し、画像の歪みを補償する方法も提案されている(非特許文献3)。
特開平10-272120号公報 特開2010-42275号公報 特開2006-149564号公報
A.M.Takahashi,Reduction of Slice Select Artifacts in Half Pulse Excitations used in Ultrashort TE(UTE)Imaging,ISMRM2010-4955 S.H.Cho et al,Compensation of eddy current by an R-L-C circuit model of the gradient system,Proc.Intl.Soc.Mag.Reson.Med.16:1156(2008) Ethan K.Brodsky et al,Characterizing and Correcting Gradient Errors in Non-Cartesian Imaging:Are Gradient Errors Linear Time-Invariant(LTI)?,Magnetic Resonance in Medicine 62:1466-1476(2009)
　しかしながら、上述した出力傾斜磁場波形の補償技術にはそれぞれ以下のような問題点がある。
　まず、渦電流を補償する傾斜磁場パルスを出力する方法(特許文献1、特許文献2)では、時間的・空間的に高次の複雑な変化を示す渦電流について、渦電流により生じる磁場を完全に補正する渦電流補正コイルを実現することは困難かつ高コストであるという問題がある。また、傾斜磁場電源の制御回路の振る舞いについては考慮されておらず、出力傾斜磁場波形の歪みを完全に除去できない。
　次に、高周波磁場パルスの照射タイミングを出力傾斜磁場波形に応じて変更する方法(特許文献3)では、高周波磁場パルスの照射タイミングのみを変更しているため出力傾斜磁場波形の非線形な歪みについては補償できていない。入力傾斜磁場波形を固定する方法(非特許文献1)では、特定の入力傾斜磁場波形のみを使用して高周波磁場パルスの振幅を調整することは、SAR(Specific Absorption Rate：被検体内へのエネルギ蓄積量を示す基準値)の増加を招くおそれがある。
　さらにシステム応答をRLC回路でモデル化する方法(非特許文献2)では、システムの応答がRLC回路によるモデル式に当てはまらない場合、出力傾斜磁場波形の歪みを表現しきれず、画像の歪みやゴースト等のアーチファクトを残すことになる。これは、傾斜磁場電源によりフィードバック制御などの複雑な命令がなされている場合などに起こりえる。また、モデル式の各係数が画像の見た目を基準にして決定されるため、係数を決定する作業者に応じてモデル式の近似精度がばらつくという問題がある。
　撮像後に傾斜磁場波形を実測して再構成に用いる方法(非特許文献3)では、実測の傾斜磁場波形をk空間座標の算出に用いているため、ノイズによる揺らぎが重畳することになる。ノイズを低減するためにフィルタなどの処理を施した場合には、フィルタ処理による傾斜磁場波形の歪みが問題となる。
　そこで本発明は、入力傾斜磁場波形に応じて変化する出力傾斜磁場波形の歪みを精度よく算出し、上記算出した出力傾斜磁場波形を用いて各種補正を行うことを目的とする。
　上記課題を解決する本発明は、入力傾斜磁場波形と当該入力傾斜磁場波形に対応する出力傾斜磁場波形を算出し、入力傾斜磁場波形と出力傾斜磁場波形とを用いて、該出力傾斜磁場波形に影響を与える複数の要素の応答関数の合算である応答関数を算出し、応答関数を用いて、撮像シーケンスに設定された傾斜磁場パルスの入力傾斜磁場波形から出力傾斜磁場波形を算出し、算出した出力傾斜磁場波形の計算値を用いて、各種補正を行う。補正の例としては、撮像シーケンスに設定された高周波磁場パルスの修正、及び/又は、撮像シーケンスで得られた核磁気共鳴信号のk空間座標の修正を行う。
　本発明によれば、傾斜磁場の出力波形に影響を与える複数の要素、例えば、渦電流や傾斜磁場電源の制御回路のそれぞれの応答関数を加味した応答関数を算出し、それによって傾斜磁場の出力波形を算出するので、渦電流のみならず制御回路に起因する傾斜磁場の歪みも補償することができ、高コストな補償用傾斜磁場コイル等を用いることなく、従来技術では取り除くことができなかった画像の歪みやゴースト等のアーチファクトを抑制することができる。
MRI装置の全体構成を示す図。 ディジタル信号処理装置の機能ブロック図。 第一実施形態による全体的な処理手順を示すフローチャート。 出力傾斜磁場波形の空間別測定を説明する図。 渦電流により発生する磁場による出力傾斜磁場波形の歪みを説明する図。 第一実施形態における渦電流の応答関数を求める処理手順を示すフローチャート。 傾斜磁場電源の制御回路により出力傾斜磁場波形に加えられるオーバシュートもしくはアンダーシュートを説明する図。 第一実施形態における傾斜磁場電源の制御回路の応答関数を求める処理手順を示すフローチャート。 第一実施形態の応答関数を用いて推定した出力傾斜磁場波形と、入力磁場波形及び実際の出力傾斜磁場波形とを示すグラフ。 渦電流のみを考慮した応答関数を用いて推定した出力傾斜磁場波形と、入力磁場波形及び実際の出力傾斜磁場波形とを示すグラフ。 第一実施形態における出力傾斜磁場波形の推定処理のフローチャート。 第一実施形態における入力傾斜磁場波形及び高周波磁場パルスの再計算処理手順を示すフローチャート。 再計算前後の高周波磁場パルスをそれぞれ用いて得られた画像を示す図。 第一実施形態における入力傾斜磁場波形及びk空間座標の再計算処理手順を示すフローチャート。 k空間座標の再計算処理前後の再構成画像を示す図。 第二実施形態における全体的な処理手順を示すフローチャート。 第二実施形態における応答関数算出の処理手順を示すフローチャート。 第三実施形態における全体的な処理手順を示すフローチャート。 第三実施形態の応答関数の算出処理に用いる入力傾斜磁場波形の例を示す図。 図19に示す台形波のラプラス変換波形を示す図。
　以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
　最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施形態の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体の断層画像を得るものであり、図1に示すように、静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4とを備えて構成される。
　静磁場発生系2は、図示していないが、永久磁石、常電導磁石あるいは超電導磁石の静磁場発生源を含み、被検体1が置かれる空間に均一な静磁場を発生させるものである。静磁場の方向によって、垂直磁場方式、水平磁場方式などがあり、垂直磁場方式であれば、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させる。
　傾斜磁場発生系3は、MRI装置の座標系(静止座標系)であるX、Y、Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成り、後述のシ-ケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X、Y、Zの3軸方向に傾斜磁場Gx、Gy、Gzを印加する。撮影時には、スライス面(撮影断面)に直交する方向にスライス方向傾斜磁場パルスを印加して被検体1に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード方向傾斜磁場パルスと周波数エンコード方向傾斜磁場パルスを印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。
　シーケンサ4は、高周波磁場パルスと傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、信号処理系7(ディジタル信号処理装置8)の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。パルスシーケンスは、撮像方法によって異なる種々のものがあり、それらは、予めプログラムとして磁気ディスク18に格納されている。
　送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体1に高周波磁場パルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとから成る。高周波発振器11から出力された高周波磁場パルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調された高周波磁場パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、高周波磁場パルスが被検体1に照射される。
　受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号(NMR信号)を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射された電磁波によって誘起された被検体1の応答のNMR信号が被検体1に近接して配置された高周波コイル14bで検出され、信号増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。
　信号処理系7は、各種データ処理と処理結果の表示及び保存等を行うもので、ディジタル信号処理装置8と、磁気ディスク18、光ディスク19等の外部記憶装置と、ディスプレイ20と、ディジタル信号処理装置8による処理結果やディジタル信号処理装置8が使用するデータ等を保存するROM21、RAM22などの内部メモリを有し、受信系6からのデータがディジタル信号処理装置8に入力されると、ディジタル信号処理装置8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。ディジタル信号処理装置8は、上述した画像再構成等の処理のほかに、事前計測データを用いた装置特性に関する計算や、磁気ディスク18に格納されたパルスシーケンスの実行に必要なパラメータの計算を行い、シーケンサ4の動作を制御する。
　操作部25は、MRI装置の各種制御情報や上記信号処理系7で行う処理の制御情報を入力するもので、トラックボール又はマウス23、及び、キーボード24から成る。この操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。
　なお、図1において、送信側の高周波コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置されている。また、受信側の高周波コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。
　本実施形態のMRI装置は、さらに信号処理系7の機能として、傾斜磁場コイル9への入力波形(入力傾斜磁場波形)と、それによって静磁場空間に発生する傾斜磁場波形(出力傾斜磁場波形)との関係を表す応答関数を算出する手段と、この応答関数を用いて出力傾斜磁場波形を算出する傾斜磁場波形算出手段とを備えている。これら手段を備えたディジタル信号処理装置8の詳細を図2に示す。
　図示するように、ディジタル信号処理装置8は、傾斜磁場の応答関数を計算する応答関数算出部81、算出された応答関数を用いて、実際の撮像シーケンスで使用される傾斜磁場の出力傾斜磁場波形を算出する出力傾斜磁場波形算出部82、傾斜磁場波形算出部82が算出した出力傾斜磁場波形を用いて撮像シーケンスで使用する高周波パルス等を補正するシーケンス計算部83、傾斜磁場波形算出部81が算出した出力傾斜磁場波形を用いて、NMR信号が配置されるk空間の座標位置を補正するk空間座標算出部84、NMR信号を用いて画像再構成演算を行って画像データを作成する画像再構成部85を備えている。
　以下、上記構成における本実施形態のMRI装置の動作、特にディジタル信号処理装置8で行われる処理を中心に説明する。
<第一の実施形態>
　本実施形態のMRI装置のディジタル信号処理装置8で行われる処理は、大きく分けて、4つの処理からなる。4つの処理の関係を図2のブロック図を用いて説明する。1つめの処理101は、事前測定として応答関数算出部81で行われる処理であり、渦電流及び傾斜磁場電源の制御回路の応答関数1112を算出する処理である。処理102～104は、実際の撮像において行われる処理である。処理102は、処理101で求めた応答関数1112と、撮像シーケンス毎に設定される入力傾斜磁場波形1113、1114を用いて、それぞれの出力傾斜磁場波形を推定する処理であり、傾斜磁場波形算出部82が行う。処理103は、処理102の結果を用いて、撮像シーケンス高周波磁場パルスを再計算する処理であり、シーケンス計算部83の高周波磁場パルス算出部831が行う。処理104は、処理102の結果を用いて、計測されるNMR信号1115の各サンプリング点に対応するk空間座標を再計算する処理であり、k空間座標算出部84が行う。以下、それぞれの処理について詳説する。
　<<応答関数算出処理101>>
　処理101の渦電流及び傾斜磁場電源の制御回路の応答関数算出処理について図3に示すフローチャートを用いて全体の流れを説明する。
　ステップ301において、傾斜磁場電源の制御パラメータを無効化し、フィードバック制御やフィードフォワード制御などの動作を傾斜磁場電源が行わないように設定する。傾斜磁場電源の制御を無効にした状態で、ステップ302において、出力傾斜磁場波形を空間別に測定する。出力傾斜磁場の測定は、公知技術(例えば、特許文献3に記載される技術)を用いて行うことができる。空間別の測定は、例えば図4に示すように、X、Y、Zの各軸方向に50mm間隔で測定する。この測定結果である出力傾斜磁場波形1111及び測定時に用いた入力傾斜磁場波形1110を、信号処理系のメモリRAM22(図1)に格納する。
　ステップ303において、RAM22に格納された入力傾斜磁場波形及び計測された出力傾斜磁場波形を読み出し、渦電流の応答関数を算出する。渦電流の応答関数を算出方法については後述する。次にステップ304において、傾斜磁場電源の制御パラメータを有効化し、フィードバック制御やフィードフォワード制御などが動作する状態にする。その状態において、出力傾斜磁場波形を測定する(ステップ305)。この測定結果もまた、RAM22に格納する。ステップ306において、傾斜磁場電源の制御回路の応答関数を算出する。傾斜磁場電源の制御回路の応答関数を算出方法についても後述する。以上の処理は、X、Y、Zの3軸方向に配された傾斜磁場コイル毎に行う。
　まずステップ303の渦電流の応答関数の算出について説明する。渦電流により生じる磁場は、図5に例として示すように、入力傾斜磁場波形501に対して出力傾斜磁場波形502を歪ませる。入力傾斜磁場波形と出力傾斜磁場波形の関係は例えば、次のような式で表現できる。
Figure JPOXMLDOC01-appb-I000001
　ここで、sはラプラス変数、Gout(s)は出力傾斜磁場波形のラプラス変換形、Gin(s)は入力傾斜磁場波形のラプラス変換形、H(s)は渦電流の伝達特性である。なお、出力傾斜磁場波形は、図4に示したように空間別に測定された個々のデータである。式(1)は一例であり、例えば下記式のような別の表現に置き換えることもできる。
Figure JPOXMLDOC01-appb-I000002
　以下の説明では、式(1)の表現を用いることとする。入力傾斜磁場波形Gin(s)と出力傾斜磁場波形Gout(s)が既知であることから、理論的には渦電流の伝達特性H(s)は次式から算出できる。
Figure JPOXMLDOC01-appb-I000003
　渦電流の応答関数は、式(3.1)を逆ラプラス変換することで得られる。
Figure JPOXMLDOC01-appb-I000004
　ここで、ILT[]は逆ラプラス変換を表す。ただし、式(3.1)の計算では、Gin(s)の値が零となる場合にはH(s)の値が発散し、正しい値が求められない。このことを防ぐため、本実施形態では、渦電流の応答関数をモデル式に当てはめ、モデル式中の係数を検索して求める。渦電流の応答関数のモデル式の一例を下記式(4)に示す。この式(4)の応答関数は渦電流により発生する磁場の過渡応答を表現している。
Figure JPOXMLDOC01-appb-I000005
　ここで、tは時間、α1～αnは振幅ゲイン、τ1～τnは時定数である。式(4)中の指数関数の数nは、MRI装置を構成する部品が持つ時定数の種類で決まる。発明者らの経験では、nは10あれば十分に渦電流により発生する磁場の過渡応答を表現することができる。
　渦電流の応答関数を求めるために、入力傾斜磁場波形に式(4)のh(t)を畳みこんだ結果が、最も出力傾斜磁場波形に近い傾斜磁場波形を与えるα1～αn及びτ1～τnを検索する。検索手順を図6のフローチャートを用いて説明する。
　まずステップ601で、振幅ゲインおよび時定数の初期値を設定する。初期値として、ここでは、振幅ゲインαの各係数を全てゼロに設定し、時定数τにはゼロ以外の正の数値を設定する。ループステップ602(繰り返し処理)では処理のループカウンタとしてiを設定する。このiは、αとτの添字の番号に対応する。すなわち、i=1の場合、α1とτ1の最適値を検索することになる。ステップ603では、αiのみを変化させ、最適なαiの値を検索する。このとき、最適値は、下記式の評価量e1が最小になる値である。
Figure JPOXMLDOC01-appb-I000006
　ここで、gout(t)は出力傾斜磁場波形、gin(t)は入力傾斜磁場波形、h(t)は式(4)に示した応答関数である。式(5)は、出力傾斜磁場波形と入力傾斜磁場波形に応答関数を畳み込んだ結果がどの程度類似しているかについて数値化しているものである。従って、式(5)以外の評価式を用いても本質的に変わりはない。-1<αi<1の範囲内を0.001刻みでαiを変更してe1を計算し、e1が最小になる値になるαiを最適値とする。ここで、黄金分割法などに代表される検索アルゴリズムを用いてもよい。
　検索した最適値αiをh(t)に設定した上で、ステップ604において、今度はτiのみを変化させ、評価量e1が最小になるτiの値を検索する。検索方法についてはαiと同様である。αiとτiを交互に最適値を検索し、その値を更新し続けた結果、αiとτiの値が収束しているかについて分岐ステップ605で判断する。αiとτiの収束条件は例えば以下のように定義できる。
Figure JPOXMLDOC01-appb-I000007
　ここで、αi_prevとτi_prevは前回検索時に検出されたαiとτiの最適値である。すなわち、最低2回の繰り替えし検索が必要となる。
　上記収束条件を満たさなければ、再度αiとτiの検索を行う。収束条件を満たせば、iがn未満であればiをインクリメントし、次のαiとτiの組み合わせの検索に移る。iがn以上であれば終了する。なお、αは正だけではなく負の値も取り得る。αは正の値のみになることが渦電流の発生原理からは正しいが、α1～αnを順番に求め、後で求めるα及びτの成分を無視して検索しているが故に負の値が生じ得る。以上の手順で渦電流の応答関数を求める。
　次に図3のステップ306の傾斜磁場電源の制御回路の応答関数の算出について説明する。傾斜磁場電源の制御回路は、図7に例として示すように入力傾斜磁場波形701に対して出力傾斜磁場波形702にオーバシュートもしくはアンダーシュートなどの挙動を含ませることになる。傾斜磁場電源が出力する傾斜磁場パルスは、入力波形にこの制御回路からの成分を加えたものであり、下記式のように表現できる。
Figure JPOXMLDOC01-appb-I000008
　ここで、gamp(t)は傾斜磁場電源が出力する傾斜磁場パルス、gin(t)は入力傾斜磁場波形、u(t)は傾斜磁場電源の制御成分である。
　この制御成分についても、渦電流の場合と同様にモデル式で定義する。制御成分のモデル式は、制御回路が採用している制御手法に応じて適宜定義することができる。制御手法としては、PID制御、それより簡素なP制御やPD制御、PI制御、もしくはより複雑なニューラルネットワークなどの知的制御を想定することができ、本発明を適用するMRI装置に搭載される傾斜磁場電源の仕様に合わせて制御理論を選択すればよい。本実施形態では傾斜磁場電源の制御回路がPID制御を行っている場合を想定し、u(t)を下記のモデル式で定義する。
Figure JPOXMLDOC01-appb-I000009
　ここで、KPは比例ゲイン、KIは積分ゲイン、KDは微分ゲインである。また、式(9)のDiff(t)は下記式で定義される。
Figure JPOXMLDOC01-appb-I000010
　ここで、delayはPID制御の入力を算出する際の時間遅れである。gamp(t)と出力傾斜磁場波形の関係は以下である。
Figure JPOXMLDOC01-appb-I000011
　ここで、gout(t)は図3のステップ305でRAM22に格納された出力傾斜磁場波形である。h(t)は式(4)に示した渦電流の応答関数であり、上述した処理により算出され、既知である。そこで、入力傾斜磁場波形と傾斜磁場電源の制御成分u(t)との合計であるgamp(t)に、ステップ303で算出したh(t)を畳みこんだときに、最も出力傾斜磁場波形gout(t)に近い傾斜磁場波形を与えるu(t)を求める。u(t)に含まれるdelay、KP、KI、KDを検索する手順を図8のフローチャートを用いて説明する。ステップ801では、delay、KP、KI、KDを全てゼロに初期設定する。ステップ802では、delayのみを変化させ最適なdelayを検索する。このとき、最適値は、下記式の評価量e2が最小になる値である。
Figure JPOXMLDOC01-appb-I000012
　評価量e2が最小になる値になるdelayの検索は、上述した応答関数の係数α1～αn及びτ1～τnの検索処理と同様に行う。すなわち、特定の範囲内で全点検索を行ってもよいし、黄金分割法などに代表される検索アルゴリズムを用いてもよい。delayを求めることで、Diff(t)が求まることから、次にステップ803において最適なKPを検索する。ここでもやはりKPが最適となる条件は、式(12)の評価量e2が最小になる値である。同様の処理で、ステップ804においてKIを検索し、ステップ805においてKDを検索する。最終的に求められる傾斜磁場電源の制御回路の応答関数は下記式で表わされる。
Figure JPOXMLDOC01-appb-I000013
　以上の処理により、渦電流の応答関数h(t)と制御回路の応答関数(制御成分u(t))が算出されたので、入力傾斜磁場波形が与えられれば、式(11)より出力傾斜磁場波形を計算できることになる。式(11)を 解くための上記応答関数h(t)(式(4))およびu(t) のパラメータは、図4に示したような空間別に測定された個々のデータについてそれぞれ求められ、信号処理系7のRAM22に保存され、実際の撮像においてパルスシーケンスの計算や画像再構成に使用される。
　式(11)を用いて出力傾斜磁場波形を推測した結果例を図9に示す。図中、波形901は入力傾斜磁場波形、波形902は実測して得られた出力傾斜磁場波形、そして波形903は式(11)を用いて推定した出力傾斜磁場波形である。波形902と波形903はよく一致していることから、式(11)により出力傾斜磁場波形を精度良く推定できることが分かる。一方、図10は、傾斜磁場電源の制御成分u(t)を加えない場合の出力傾斜磁場波形の推定結果例である。
実測して得られた出力傾斜磁場波形902と、u(t)をゼロとした式(11)から推定した出力傾斜磁場波形904との間に矢印で示した部分に偏差が生じている。このことから、傾斜磁場電源の制御成分を加味しなければ、高精度に出力傾斜磁場波形を推定できないことが分かる。
　<<出力傾斜磁場波形算出処理102>>
　次に、図2のブロック図に戻り、出力傾斜磁場波形算出部82で行われる処理102について説明する。実際の撮像にあたっては、撮像目的に応じてシーケンス計算部83により撮像パルスシーケンスが計算され、計算された撮像パルスシーケンスに含まれる高周波パルス(振幅、エンベロープ)、スライス選択傾斜磁場パルス(入力傾斜磁場波形1113)、読出し傾斜磁場パルス(入力傾斜磁場波形1114)などがRAM22に保存される。
　処理102では、応答関数算出手段81が算出した応答関数1112と、設定された入力傾斜磁場波形1113、1114を用いて、撮像パルスシーケンスの出力傾斜磁場波形を算出する。図11に、出力傾斜磁場波形の算出処理のフローチャートを示す。
　まずステップ1101において、撮像視野の座標位置からX、Y、Zの3軸方向における視野中心位置を取得する。視野中心位置は、下式から算出される。
Figure JPOXMLDOC01-appb-I000014
　ここで、AxisはX、Y、Zの各軸を表し(以下、同じ)、FOV(Axis,min)とFOV(Axis,max)は撮像条件に応じて指定される各軸の撮像視野の最小座標位置と最大座標位置である。
　ステップ1102において、RAM22に保存された空間別の応答関数h(t)及び傾斜磁場電源の制御成分u(t)から、ステップ1101で取得した各軸の視野中心位置に対応する応答関数h(t)および制御成分u(t)を選択し、取得する。ここで、応答関数h(t)及び傾斜磁場電源の制御成分u(t)の取得は、視野中心位置に最も近い位置で算出された結果を取得してもよいし、空間別に算出された結果を補間して応答関数h(t)及び傾斜磁場電源の制御成分u(t)を求めてもよい。補間する場合は、下記式に従う。
Figure JPOXMLDOC01-appb-I000015
　ここで、h(Axis,t)とu(Axis,t)は各軸の視野中心位置における応答関数と傾斜磁場電源回路の制御成分である。また、P1およびP2は、視野中心位置の前後方向でそれぞれ最も近い応答関数および傾斜磁場電源回路の制御成分の算出位置を表す。
　次に、ステップ1103において、X、Y、Zの3軸方向に対する入力傾斜磁場波形を取得する。入力傾斜磁場波形は、撮像パルスシーケンスに応じて設定され、RAM22に保存されているものであり、図2に示される、スライス傾斜磁場パルス1113もしくは読み出し傾斜磁場パルス1114である。どちらの場合であっても、処理の内容は同じである。
　ステップ1104において、次式で表される3軸の入力傾斜磁場波形を合成した波形を算出する。
Figure JPOXMLDOC01-appb-I000016
　ここで、gin(Axis,t)は取得した各軸の入力傾斜磁場波形である。
　ステップ1105において、軸毎の入力傾斜磁場波形に対して応答関数h(t)及び傾斜磁場電源の制御成分u(t)を用いて出力傾斜磁場波形を推定する。出力傾斜磁場波形の推定は、式(11)を3軸に拡張した形として下記式で表わされる。
Figure JPOXMLDOC01-appb-I000017
　ここで、gest(Axis,t)は計算により求めた各軸の出力傾斜磁場波形の推定結果である。
各軸の出力傾斜磁場波形を算出した後に、ステップ1106において、下記式に従い各軸の波形を合成する。
Figure JPOXMLDOC01-appb-I000018
　以上のように推定されたスライス傾斜磁場及び読出し傾斜磁場の各出力傾斜磁場波形は、それぞれ、高周波磁場パルスの再計算(図2：103)およびNMR信号のk空間座標位置の補正(図2：処理104)に用いられる。
　<<高周波磁場パルス算出処理103>>
　処理103では、出力傾斜磁場波形算出部82で求めたスライス傾斜磁場の出力傾斜磁場波形gest(t)と入力傾斜磁場波形gin(t)を用いて、高周波磁場パルスの再計算を行う。高周波磁場パルスの再計算は下記式に従い、高周波パルスのエンベロープおよび周波数を計算する。
Figure JPOXMLDOC01-appb-I000019
　ここで、RFenv'(t)は再計算後の高周波磁場パルスのエンベロープ、RFenv(t)は再計算前の高周波磁場パルスのエンベロープ、RFfreq'(t)は再計算後の高周波磁場パルスの周波数、RFfreq(t)は再計算前の高周波磁場パルスの周波数である。
　これらの式(19)～(21)からわかるように、ここでは、スライス傾斜磁場の出力傾斜磁場波形の強度に比例して高周波磁場パルスのエンベロープの振幅、周波数、照射時間を変更している。これら式を基本として、装置の制約や撮像条件を考慮した変更を行うことができる。
　変更を要する一例として、高周波磁場パルスの周波数を時間に応じて変化させることが装置の構成上困難な場合がある。その場合には、高周波磁場パルスの位相を時間と共に変化させ、結果として周波数の変調を実現してもよい。周波数変調を実現するための位相の算出は次式に従う。
Figure JPOXMLDOC01-appb-I000020
　ここで、g_standardは高周波磁場パルスの周波数計算の基準となる傾斜磁場強度[T/m]、distanceは磁場中心から励起位置までの距離[m]、γは磁気回転比[Hz/T]、Δtは高周波磁場パルスの位相変更間隔時間[s]である。
　ただし、位相を変更できる範囲は0～360[deg]であるため、distanceの範囲もまた限られる。従って、distanceの設定可能な範囲を広くするためには、位相の変化は少ない方が好ましい。そこで、位相の変化を最小とするために、g_standard及び高周波磁場パルスの周波数を下記式に従って計算する。
Figure JPOXMLDOC01-appb-I000021
　ここで、Average[　]は高周波磁場パルス印加区間内での平均値計算処理、RFfreq'は時間的変化のない再計算後の高周波磁場パルスの周波数、RFfreqは時間的変化のない再計算前の高周波磁場パルスの周波数である。式(23)及び式(24)によりg_standardとRFfreq'を算出することで、位相の変化が最小に抑えられ、distanceの上限値を大きく取ることが可能となる。
　変更を要する他の例として、高周波磁場パルスの照射時間の抑制がある。上記式(19)～式(21)を用いた高周波磁場パルスの再計算の結果として、高周波磁場パルスの照射時間は変化するが、パルスシーケンスによっては、照射時間の変化、特に照射時間の延長が好ましくない場合がある。照射時間を変化させない高周波磁場パルスの再計算手順を図12に示すフローチャートを用いて説明する。
　ステップ1021と1022において、軸毎の出力傾斜磁場波形の算出、各軸の出力傾斜磁場波形の合成を行う。これらの処理は、出力傾斜磁場波形算出部82で行われるものであり、図11のステップ1103とステップ1104で行う処理と同じである。次に、ステップ1031において、式(19)～式(21)に従い高周波磁場パルスを計算する。その結果を受け、ステップ1032において、照射時間の増加率βRFを算出する。増加率βRFの算出は次式に従う。
Figure JPOXMLDOC01-appb-I000022
　ここで、gin(t)は入力傾斜磁場波形、gest(t)は式(18)より求めた出力傾斜磁場波形の推定結果、StartTimeは高周波磁場パルス照射開始時刻、EndTimeは高周波磁場パルス照射終了時刻を表す。
　分岐ステップ1033において、増加率βRFが1.0を越えるか否かを判断する。増加率βRFが1.0を越える場合、計算後の高周波磁場パルスの照射時間が延長していることになる。
その場合にはステップ1034に移り、スライス傾斜磁場の入力傾斜磁場波形の再作成を行う。入力傾斜磁場波形の再作成は次式に従う。
Figure JPOXMLDOC01-appb-I000023
　照射時間の増加率βRFを入力傾斜磁場波形に掛けることで、入力傾斜磁場波形の強度を増加させ、それにより高周波磁場パルスの印加時間を短縮させる。この時、入力傾斜磁場波形の全体を一様に増加させるだけであるため、特別なシムコイルや傾斜磁場コイルなどは必要としない。
　入力傾斜磁場波形の再作成後、その結果を入力として処理102のステップ1021に戻し、出力傾斜磁場波形の算出(ステップ1021、1022)、高周波磁場パルスの再計算(ステップ1031)、照射時間の増加率βRFの算出(ステップ1032)を行う。再び算出した増加率βRFが1.0以下であれば処理を終了し、1.0以下になっていなければ再度、入力傾斜磁場波形を作成する。このとき、式(26)のgin(t)は直前に再作成された入力傾斜磁場波形である。
　このようにスライス傾斜磁場の入力傾斜磁場波形の再作成を繰り返すことで、最終的に高周波磁場パルスの照射時間は、再計算前の高周波磁場パルスの照射時間に比べて延長しないようにすることができる。これにより、再計算後の高周波磁場パルスを使用する際に、スライス傾斜磁場パルス以外のパルスシーケンスを変更する必要がなくなる。
　図13は、再計算前の高周波磁場パルスを用いて得られた画像1301と、再計算後の高周波磁場パルスを用いて得られた画像1302である。どちらの画像も、VERSE法を用いた超短TE計測での撮像結果である。再計算前の高周波磁場パルスを用いて得られた画像1301に比べて、再計算後の高周波磁場パルスを用いて得られた画像1302は、矢印で示す背景信号(撮影断面以外からの信号)の混入が少ないことから、良好な励起プロファイルが得られていることが分かる。
　<<k空間座標再計算処理104>>
　続いて、図2のk空間座標算出部84が行うk空間座標再計算処理104について説明する。k空間座標再計算処理104では、出力傾斜磁場波形算出部82で求めた読出し傾斜磁場の出力傾斜磁場波形gest(t)と入力傾斜磁場波形gin(t)を用いてNMR信号が配置されたk空間座標の再計算を行う。図14に、k空間座標の再計算手順を示すフローチャートの一例を示す。この処理手順も、図12に示す高周波磁場パルスの再計算のフローチャートと同様の処理の流れで、k空間座標の再計算結果によって読出し傾斜磁場の入力傾斜磁場波形を再作成する。
　まずステップ1041において、処理102において求めた、読出し傾斜磁場の入力傾斜磁場波形gin(t)と出力傾斜磁場波形gest(t)を用いて、k空間座標の算出を行う。出力傾斜磁場波形gest(t)を用いたk空間座標の再計算は下記式に従う。
Figure JPOXMLDOC01-appb-I000024
　ここで、γは磁気回転比[Hz/T]である。この式は、一般的な入力傾斜磁場波形gin(t)を用いたk空間座標の算出の式(次式(28))におけるgin(j)をgest(j)に置き換えたものである。
Figure JPOXMLDOC01-appb-I000025
　次に、再計算によって算出した座標の範囲が縮小しているか否かを調べる。再計算したk空間座標k’(t)は、もともとのk空間座標k(t)に対して、その範囲が縮小する可能性がある。その場合、k空間内の想定された範囲内に信号が充填されないことになり、画像の分解能を低下させる。そこで、ステップ1042において、k空間の縮小率βkを算出する。
縮小率βkの算出は次式に従う。
Figure JPOXMLDOC01-appb-I000026
　ここで、EndTimeは読み出し傾斜磁場パルスの印加終了時刻(もしくはNMR信号の計測終了時刻)を表す。
　分岐処理ステップ1043において、縮小率βkが1.0未満であるか否かを判断する。増加率βkが1.0未満である場合、再計算後のk空間座標k’(t)は再計算前のk空間座標k(t)に対して縮小していることになる。その場合にはステップ1044に移り、読み出し傾斜磁場の入力傾斜磁場波形の再作成を行う。入力傾斜磁場波形の再作成は次式に従う。
Figure JPOXMLDOC01-appb-I000027
　照射時間の縮小率βkの逆数を入力傾斜磁場波形に掛けることで、入力傾斜磁場波形の強度を増加させ、それによりk空間座標の範囲を広げる。この時、式(30)からわかるように、入力傾斜磁場波形は全体を一様に増加させられるだけであるため、特別なシムコイルや傾斜磁場コイルなどは必要としない。入力傾斜磁場波形の再作成後、その結果を入力として処理102に戻し、出力傾斜磁場波形の算出(ステップ1021、1022)、k空間座標の再計算(ステップ1041)、縮小率βkの算出(ステップ1042)を行う。
　分岐処理ステップ1043において、再び算出した縮小率βkが1.0以上であれば処理を終了し、1.0以上になっていなければ再度、入力傾斜磁場波形を作成する。このとき、式(30)のgin(t)は直前に再作成された入力傾斜磁場波形である。
　このように入力傾斜磁場波形の再作成を繰り返すことで、最終的にk空間座標の範囲は縮小せず、画像の分解能を保つことができる。図15は、再計算前のk空間座標を用いて得られた画像1501と、再計算後のk空間座標を用いて得られた画像1502である。どちらの画像も、EPI計測での撮像結果であり、読み出し傾斜磁場は画像における横方向に印加されている。再計算前のk空間座標を用いて得られた画像1501に比べて、再計算後のk空間座標を用いて得られた画像1502は、矢印で示す画像の歪みが低減されていることから、読み出し傾斜磁場の歪みの補償が適切に行われていることが分かる。
　撮像は、上記処理103で修正された高周波磁場パルスを用いて行われ、また上記処理104で修正されたk空間座標を用いて、取得したNMR信号の画像再構成を行う。これにより選択したスライスへの外部信号の混入がなく、且つ画像の歪みやアーチファクトの抑制された画像を得ることができる。
　以上、説明した本実施形態の、従来技術に対する利点は、次のとおりである。
　まず、渦電流を補償する傾斜磁場パルスを出力する方法(特許文献1、特許文献2)では、時間的・空間的に高次の複雑な変化示す渦電流による磁場を完全に補正することは困難かつ高コストであるのに対し、本実施形態では特別な磁場補正用コイルを必要とせず、渦電流による磁場を含めた高周波磁場パルスとk空間座標の再計算を行うことで、励起プロファイルの劣化と画像の歪みやゴースト等のアーチファクトを抑制できる。
　また、高周波磁場パルスの照射タイミングを出力傾斜磁場波形に応じて変更する方法(特許文献3)では、高周波磁場パルスの照射タイミングのみを変更しているのに対して、本実施形態では渦電流及び傾斜磁場電源の制御回路の応答関数を用いて出力傾斜磁場波形を推定するため、出力傾斜磁場波形の非線形な歪みについても補償する事が可能である。
　固定の入力傾斜磁場波形を用いる方法(非特許文献1)では、特定のパルスシーケンスのみにしか適用できないのに対して、本実施形態ではあらゆる入力傾斜磁場波形に対する出力傾斜磁場波形を推定できるため、固定の入力傾斜磁場波形を使用しなければならないといった制限がない。
　さらに、システム応答をRLC回路でモデル化する方法(非特許文献2)では、傾斜磁場電源の制御回路の振る舞いについては考慮されていないが、本実施例では傾斜磁場電源の制御回路の振る舞いも考慮するため、より高精度に高周波磁場パルス及びk空間座標を再計算できる。加えて、システム応答をRCL回路でモデル化する方法(非特許文献2)では、画像の見た目を判断基準としてモデル式の各係数を求めているが、本実施例では出力傾斜磁場波形に基づいてモデル式の係数を算出するため、人間の判断に依る必要がなく、精度の高いモデル式を算出することができる。
　また、撮像後に傾斜磁場波形を実測して再構成に用いる方法(非特許文献3)では、ノイズによる揺らぎが重畳することになるが、本実施例ではモデル式を用意して渦電流と傾斜磁場電源回路の制御成分を表現しているため、推定する出力傾斜磁場波形にノイズが重畳しない。
　<第二の実施形態>
　本実施形態においても、図2に示すディジタル信号処理部8の構成と、ディジタル信号処理部8で4つの特徴的な処理101～104が行われることは、第一実施形態と同様である。
本実施形態は、渦電流及び傾斜磁場電源の制御回路の応答関数を算出する処理101を、第一実施形態より簡易な手続きで求めることを特徴としている。すなわち、第一実施形態では、傾斜磁場電源回路の制御を無効にした状態と有効にした状態でそれぞれ2回の出力傾斜磁場波形の測定を行ったが、本実施形態では、1回の出力傾斜磁場波形の測定のみ行う。本実施形態の処理の流れを図16のフローチャートを用いて説明する。
　まず、ステップ1601において、出力傾斜磁場波形を空間別に測定する。この時、傾斜磁場電源の制御パラメータには変更を加えず、装置の通常状態としておく。すなわち、傾斜磁場電源の制御パラメータが有効になっている状態とする。出力傾斜磁場波形の測定結果は、RAM22に格納する。この処理は、第一実施形態のフローを示す図3のステップ302、305で行われる処理と同様である。
　次に、ステップ1602において、RAM22に格納された入力傾斜磁場波形及び計測された出力傾斜磁場波形を読み出し、渦電流と傾斜磁場電源の制御回路の応答関数を算出する。
　渦電流と傾斜磁場電源の制御回路の応答関数の算出について、図17のフローチャートを用いて説明する。ステップ1701において渦電流の応答関数を算出する。この処理は上述した第一実施形態でのステップ303の処理と同様である。すなわち、次式(31)の評価量e1が最小になる応答関数h(t)を検索する。
Figure JPOXMLDOC01-appb-I000028
　ここで、gout(t)はステップ1601でRAM22に格納された出力傾斜磁場波形であり、gin(t)は入力傾斜磁場波形である。h(t)は式(4)で定義される渦電流の応答関数であり、第一実施形態でのステップ303と同様に式(4)の係数α1～αn、τ1～τnを検索する。
　式(31)は、第一実施形態の式(5)と同形状であるが、式(5)の出力傾斜磁場波形gout(t)は傾斜磁場電源の制御回路を無効にして測定したものであるのに対し、式(31)の出力傾斜磁場波形gout(t)には、傾斜磁場電源の制御回路による歪みが含まれている。そのため、算出される渦電流の応答関数h(t)には、渦電流だけではなく傾斜磁場電源の制御回路の特性も含まれることになる。
　次に、ステップ1702において傾斜磁場電源の制御回路の応答関数を算出する。この処理は第一実施形態のステップ306の処理と同じである。すなわち式(32)の評価量e2が最小になるgamp(t)を検索する。式(32)のgamp(t)は、式(33)で定義される。
Figure JPOXMLDOC01-appb-I000029
　ここでu(t)は、式(9)で定義される傾斜磁場電源の制御成分であるが、ステップ1701で求めた渦電流の応答関数h(t)の中には、傾斜磁場電源の制御回路の応答成分も含まれているため、渦電流の応答関数のモデル(式(4))では表現しきれなかった成分を式(9)に当てはめて算出することになる。
　ステップ1702で傾斜磁場電源の制御回路の応答関数u(t)を求めた後に、ステップ1703において渦電流の応答関数h(t)を再計算する。これは、傾斜磁場電源の制御回路の応答関数u(t)を定義することで、電流の応答関数h(t)(式(4))の各係数(α1～αn及びτ1～τn)の最適値が変化するためである。応答関数h(t)を再計算は、式(31)のgin(t)を、ステップ1702で求めたgamp(t)に置き換えて、その評価量e1を最小とするh(t) の各係数を検索することにより行う。以上の処理により、渦電流の応答関数と傾斜磁場電源の制御成分を加味した応答関数が求められる。応答関数は、空間別に測定された個々のデータについてそれぞれ求められ、信号処理系7のRAM22に保存される。
　撮像に際し、RAM22に保存した応答関数を用いて、出力傾斜磁場波形の推定する処理102、及び高周波磁場パルスやk空間座標を再計算する処理103、104の手順は第一実施形態と同様であり、説明を省略する。
　本実施形態によれば、応答関数の算出回数が増えるものの、出力傾斜磁場波形の測定回数を減らすことができ、傾斜磁場電源の制御パラメータを変更する作業も省くことができ、より簡素な手順で励起プロファイルを再計算するために必要な応答関数を算出できる。
その他の効果は、第一実施形態と同様である。
　<第三実施形態>
　本実施形態においても、図2に示すディジタル信号処理部8の構成と、ディジタル信号処理部8で4つの処理101～104が行われることは、第一実施形態と同様である。本実施形態は、渦電流及び傾斜磁場電源の制御回路の応答関数を特定のモデル式に当てはめずに求めることを特徴としている。本実施形態の処理の流れを図18に示す。図18のステップ1801、1802は、第二実施形態で説明した図16のステップ1601、1602と同様であるが、本実施形態の特有の処理は、ステップ1801において、複数種類の入力傾斜磁場波形を用いて出力傾斜磁場波形を測定する処理、ステップ1802において、モデル式を用いずに渦電流と傾斜磁場電源制御回路の応答関数を直接導出する処理、ステップ1803において、導出した応答関数を補正する処理である。
　既に説明したように、入力傾斜磁場波形Gin(s)(ラプラス変換形)と出力傾斜磁場波形Gout(s) (ラプラス変換形)が既知であれば、理論的には渦電流の伝達特性H(s)は式(3.1)から算出できる。式(3.1)を再掲する。
Figure JPOXMLDOC01-appb-I000030
　しかし、Gin(s)の値が零となる場合にはH(s)の値が発散し、正しい値が求められない。
これを防ぐため、第一および第二実施形態では、渦電流の応答関数をモデル式に当てはめて、最適な評価値を与える応答関数を導出した。
　本実施形態では、モデル式を用いることなく、式(3.1)或いはそれを逆ラプラス変換した式(3.2)から算出される応答関数h(t)を用いる。この際、零点の影響を排除するために、渦電流と傾斜磁場電源制御回路の応答関数の算出に用いる入力傾斜磁場波形を複数種類用意する。
　具体的には、ステップ1801において、複数種類の入力傾斜磁場波形を用いて、それぞれ出力傾斜磁場波形を測定し、RAM22に保存する。入力傾斜磁場波形なるべく急峻な形状を持ったものがよい。なだらかな変化の波形では、高周波域が含まれず、高域の応答を算出できないためである。理想的には、インパルス波を用いることが望ましいが、傾斜磁場電源の出力性能などの制約によりインパルス波を出力することが難しいため、本実施形態では、図19に示すような複数の台形波を用いる。
　ステップ1802では、式(3.1)により伝達特性H(s)を算出する。その際、一つの入力傾斜磁場波形(台形波)のスペクトル(ラプラス変換形)に現れる零点を、他の入力傾斜磁場波形(台形波)のスペクトルデータを用いて置き換えることにより、零点を回避する。具体的には、次式に従う。
Figure JPOXMLDOC01-appb-I000031
　ここで、Gin1(s)は、台形波1901のラプラス変換形、Gout1(s)は、入力傾斜磁場波形として台形波1901を用いた場合の出力傾斜磁場波形のラプラス変換形、Gin2(s)は、台形波1902のラプラス変換形、Gout2(s)は、入力傾斜磁場波形として台形波1902を用いた場合の出力傾斜磁場波形のラプラス変換形である。
　図19に示す台形波群をラプラス変換して得られる振幅スペクトルの概略を図20に図示する。図20に示すように、台形波1901、1902のスペクトル2001、2002に零点が現れる位置は、スペクトル2001とスペクトル2002とで異なるので、式(34)に示すように、値が大きい方のスペクトルを用いて計算することにより、零点を回避することができる。なお、図19では台形波を2つ示しているが、3つ以上の台形波を用意してもよい。その場合には、式(34)において、複数の入力傾斜磁場波形のうち、スペクトルの値が最も大きい(零点ではない)値となる入力傾斜磁場波形とそれによる出力傾斜磁場波形の組合せを選んで伝達特性の算出に用いればよい。
　次に、ステップ1803において、ステップ1802で算出した応答関数h(t)からノイズの影響を排除する処理を行う。これは、Gout1(s)及びGout2(s)は実測データであるため、ノイズが重畳しており、正しい伝達特性が求められないケースが存在するためである。本実施形態では、このようなノイズを、ラプラス空間で2段階のフィルタ処理を行うことにより除去する。
　まず式(34)で求めた伝達特性に、次式の最適フィルタを掛け算し、ノイズの影響を抑制する。
Figure JPOXMLDOC01-appb-I000032
　ここで、H'(s)は最適フィルタ適用後の伝達特性、Φ(s)は最適フィルタである。最適フィルタΦ(s)には、公知の技術であるウィナーフィルタを用いる。ウィナーフィルタの定義式は下記である。
Figure JPOXMLDOC01-appb-I000033
　ここで、N(s)はGout1(s)のノイズ成分、S(s)はGout1(s)のノイズを除いた信号成分である。N(s)及びS(s)の求め方は、既知の技術に従い、Gout1(s)のパワースペクトルの最も高周波の信号値をN(s)とし、Gout1(s)のパワースペクトルからN(s)を減算したものをS(s)とする。
　このような最適フィルタによりある程度ノイズの影響は低減できるが、実用的には最適フィルタ適用後にもノイズの影響が問題になることがある。そこで、最適フィルタ適用後の伝達特性に対してノイズ除去フィルタを実行する。本実施例ではノイズ除去フィルタとしてメディアンフィルタを用いる。
Figure JPOXMLDOC01-appb-I000034
　ここで、H”(s)はメディアンフィルタ適用後の伝達特性である。
　最終的に、渦電流と傾斜磁場電源制御回路を合わせた応答関数は、式(37)を逆ラプラス変換することで得られる。
Figure JPOXMLDOC01-appb-I000035
　ここで、ILT[]は逆ラプラス変換を表す。
　以上のようにして求めた応答関数を用いて、出力傾斜磁場波形を推定し、高周波磁場パルス・k空間座標を再計算する手順は第一実施形態と同様である。
　このように本実施形態では、実測データに含まれるノイズをラプラス空間における2つのフィルタ処理により低減することで、応答関数を特定のモデル式に当てはめることなく求めることができ、様々なハードウェア構成(傾斜磁場コイル等)やソフトウェア構成(傾斜磁場電源の制御方式)に対して柔軟に適用することができる。
　以上、本発明のMRI装置の動作、特に傾斜磁場の応答関数算出部の動作の各実施形態を説明したが、本発明は、上述した実施形態に限定されず、種々の変更を行うことが可能である。例えば、渦電流による出力傾斜磁場波形の歪みや傾斜磁場電源の制御回路による出力傾斜磁場波形の歪みを表現するモデル式は、実施形態で例示したものに限らず、他のモデル式を用いることも可能である。また上述した実施形態では、出力傾斜磁場波形を歪ませる原因として、渦電流と傾斜磁場電源の制御回路を取り上げ、それらによる応答関数を出力傾斜磁場波形の算出に用いる場合を説明したが、その他の要素、例えば、永久磁石を用いたMRI装置で見られるヒステリシスを伴う残留磁場といったような他の成分を出力傾斜磁場波形の算出に加えることもできる。また、各実施形態の説明のために示したフローチャートの処理手順は一例であり、省略できる処理もあるし、また必要に応じて別の処理を追加することも可能である。
　本発明は、スライス傾斜磁場について算出した出力傾斜磁場波形により高周波磁場パルスを計算する場合や、読出し傾斜磁場について算出した出力傾斜磁場波形により、NMR信号が配置されるk空間座標を補正する場合に好適に適用することができるが、位相エンコード傾斜磁場についても本発明の適用を妨げるものではない。また読出し傾斜磁場には、読出し傾斜磁場と位相エンコード傾斜磁場の区別のないラジアルスキャンやスパイラルスキャンのようなパルスシーケンスにおける2方向の傾斜磁場も含まれる。
　さらに上記実施形態では、スライス傾斜磁場及び読出し傾斜磁場の両方について出力傾斜磁場波形を推定し、高周波磁場パルスの再計算およびk空間座標の再計算をともに行う場合を説明したが、いずれか一方のみを行う場合も本発明に含まれる。
　以上の各実施形態の説明で明らかになる本願発明の概要は以下の通りである。即ち、
　本発明のMRI装置は、
　静磁場発生手段、傾斜磁場発生手段、高周波磁場発生手段および高周波磁場検出手段を備えた撮像手段と、撮像シーケンスに基き前記撮像手段の動作を制御する制御手段と、前記高周波磁場検出手段が検出した核磁気共鳴信号を用いて、画像再構成を含む演算を行う演算手段と、を備え、
　前記演算手段は、前記傾斜磁場発生手段に入力される入力波形と当該入力波形によって前記傾斜磁場発生手段が発生した出力波形とを用いて、前記出力波形に影響を与える複数の要素の応答関数の合算である応答関数を算出する応答関数算出手段と、前記応答関数計算手段が計算した応答関数を用いて、前記撮像シーケンスに設定された傾斜磁場パルスの入力波形から出力波形を算出する傾斜磁場算出手段と、前記傾斜磁場算出手段が算出した出力波形の計算値を用いて、前記撮像シーケンスに設定された高周波磁場パルスの修正および／または前記撮像シーケンスで得られた核磁気共鳴信号のk空間座標の修正を行う補償手段とを備えたことを特徴とする。
　好ましくは、
　前記応答関数算出手段は、前記出力波形を変化させる複数の要素毎に応答関数を算出し、前記傾斜磁場算出手段は、前記複数の要素の応答関数の合算である応答関数を用いて出力波形の算出を行う。
　また、好ましくは、
　前記応答関数算出手段は、前記要素毎に応答関数モデルを設定するモデル設定手段と、前記モデル設定手段が設定した応答関数モデルによって算出した出力波形の計算値と、出力波形の実測値との差が最小となるように応答関数モデルを適正化するモデル適正化手段とを備え、前記モデル適正化手段が適正化した応答関数モデルを各要素の応答関数とする。
　また、好ましくは、
　前記応答関数算出手段は、前記入力値として複数の入力値を用い、複数の入力値とそれに対応する実測値との組合せを用いて、入力値と出力波形の実測値との関係から、直接、応答関数を算出する。
　また、好ましくは、
　前記応答関数算出手段は、算出後の応答関数にフィルタをかけるフィルタ手段を備える。
　また、好ましくは、
　前記フィルタ手段は、ウィナーフィルタ、メディアンフィルタを含む。
　また、好ましくは、
前記傾斜磁場発生手段は、発生した傾斜磁場の出力波形に基き入力波形を制御する制御回路を備え、
　前記出力波形を変化させる複数の要素は、前記制御回路による制御に伴う前記出力波形の歪みを含む。
　また、好ましくは、
　前記応答関数算出手段は、前記静磁場発生手段が発生する静磁場の中心に対する位置が異なる複数の傾斜磁場について、それぞれ、応答関数を算出する。
　また、好ましくは、
　前記傾斜磁場算出手段が算出する傾斜磁場パルスは、前記撮像シーケンスにおけるスライス選択傾斜磁場であり、前記補償手段は、前記スライス選択傾斜磁場の計算値に基き、前記撮像シーケンスに設定された高周波パルスを補正する高周波パルス補正手段を備える。
　また、好ましくは、
　前記高周波パルス補正手段は、補正前後の高周波パルスの照射時間の変化を算出する照射時間計算手段を備え、前記照射時間計算手段が算出した照射時間の変化が所定の増加率を超えるときに、前記傾斜磁場算出手段の入力値を変更し、前記傾斜磁場算出手段が前記変更された入力値から推定した出力波形の推定値を用いて、前記高周波パルスの補正を行う。
　また、好ましくは、
　前記傾斜磁場算出手段が算出する傾斜磁場パルスは、前記撮像シーケンスにおける読出し傾斜磁場であり、前記補償手段は、前記読出し選択傾斜磁場の計算値に基き、前記撮像シーケンスで得た核磁気共鳴信号が配置されるk空間座標を補正するk空間座標補正手段を備る。
　また、好ましくは、
　前記補正手段は、補正前後のk空間座標の大きさの変化を算出するk空間サイズ計算手段を備え、前記k空間サイズ計算手段が計算したk空間の大きさの変化が所定の縮小率を超えるときに、前記傾斜磁場算出手段の入力値を変更し、前記傾斜磁場算出手段が変更後の入力値から算出した出力波形の再計算値を用いて、前記k空間座標の補正を行う。
　本発明の磁気共鳴イメージング装置の傾斜磁場発生手段が発生する出力傾斜磁場波形を推定する方法は、
　前記傾斜磁場発生手段の入力値およびそれによって前記傾斜磁場発生手段が発生した出力傾斜磁場波形の実測値を入力するステップと、
　前記出力傾斜磁場波形に影響を与える複数の因子について、それぞれ、応答関数モデルを設定するステップと、
　前記入力した傾斜磁場波形の入力値および実測値を用いて、各応答関数モデルを適正化するステップと、
　所望の傾斜磁場の入力値に対し、前記適正化した応答関数モデルを用いて、出力傾斜磁場波形の推定値を算出するステップと、
　を含むことを特徴とする。
　好ましくは、
　前記複数の因子は、傾斜磁場により発生する渦電流と、前記傾斜磁場発生手段に備えられた制御回路の出力とを含む。
　また、好ましくは、
　前記実測値は、前記傾斜磁場発生手段の制御回路を無効にした時の第1の実測値と、前記制御回路を有効にした第2の実測値を含み、
　前記応答関数モデルを適正化するステップは、前記第1の実測値および第2の実測値を用いて、制御回路による出力傾斜磁場波形の応答関数を求める。
　また、好ましくは、
　前記傾斜磁場発生手段の少なくとも2種の入力値およびそれらによって前記傾斜磁場発生手段が発生した出力傾斜磁場波形の実測値を入力するステップと、
　前記少なくとも2種の入力値およびそれに対応する実測値を用いて、前記出力傾斜磁場波形に影響を与える複数の因子の影響を含む前記傾斜磁場発生手段の応答関数を算出するステップと、
　所望の傾斜磁場の入力値に対し、前記応答関数を用いて、出力傾斜磁場波形の推定値を算出するステップと、を含む。
　本発明によれば、MRI装置において、多大な負担を伴うことなく、種々の要因によって生じる出力傾斜磁場波形の歪みを高精度に補償することができ、これにより、画像の歪みやゴースト等のアーチファクトのない信頼性の高い画像を提供することができる。
　1　被検体、2　静磁場発生系、3　傾斜磁場発生系、4　シーケンサ、5　送信系、6　受信系、7　信号処理系、8　ディジタル信号処理部(CPU)、9　傾斜磁場コイル、10　傾斜磁場電源、11　高周波発信器、12　変調器、13　高周波増幅器、14a　高周波コイル(送信コイル)、14b　高周波コイル(受信コイル)、15　信号増幅器、16　直交位相検波器、17　A/D変換器、18　磁気ディスク、19　光ディスク、20　ディスプレイ、21　ROM、22　RAM、23　トラックボール又はマウス、24　キーボード、25　操作部、81　応答関数算出部、82　傾斜磁場波形算出部、83　シーケンス計算部、84　k空間座標算出部、85　画像再構成部。

Claims (16)

  1. 　静磁場発生手段、傾斜磁場発生手段、高周波磁場発生手段および高周波磁場検出手段を備えた撮像手段と、撮像シーケンスに基き前記撮像手段の動作を制御する制御手段と、前記高周波磁場検出手段が検出した核磁気共鳴信号を用いて、画像再構成を含む演算を行う演算手段と、を備えた磁気共鳴イメージング装置であって、
    　前記演算手段は、前記傾斜磁場発生手段に入力される入力波形と当該入力波形によって前記傾斜磁場発生手段が発生した出力波形とを用いて、前記出力波形に影響を与える複数の要素の応答関数の合算である応答関数を算出する応答関数算出手段と、前記応答関数計算手段が計算した応答関数を用いて、前記撮像シーケンスに設定された傾斜磁場パルスの入力波形から出力波形を算出する傾斜磁場算出手段と、前記傾斜磁場算出手段が算出した出力波形の計算値を用いて、前記撮像シーケンスに設定された高周波磁場パルスの修正および/または前記撮像シーケンスで得られた核磁気共鳴信号のk空間座標の修正を行う補償手段とを備えたことを特徴とする磁気共鳴イメージング装置。
  2. 　請求項1に記載の磁気共鳴イメージング装置であって、
    　前記応答関数算出手段は、前記出力波形を変化させる複数の要素毎に応答関数を算出し、前記傾斜磁場算出手段は、前記複数の要素の応答関数の合算である応答関数を用いて出力波形の算出を行うことを特徴とする磁気共鳴イメージング装置。
  3. 　請求項2に記載の磁気共鳴イメージング装置であって、
    　前記応答関数算出手段は、前記要素毎に応答関数モデルを設定するモデル設定手段と、前記モデル設定手段が設定した応答関数モデルによって算出した出力波形の計算値と、出力波形の実測値との差が最小となるように応答関数モデルを適正化するモデル適正化手段とを備え、前記モデル適正化手段が適正化した応答関数モデルを各要素の応答関数とすることを特徴とする磁気共鳴イメージング装置。
  4. 　請求項1に記載の磁気共鳴イメージング装置であって、
    　前記応答関数算出手段は、前記入力波形として複数の入力値を用い、複数の入力値とそれに対応する実測値との組合せを用いて、入力値と出力波形の実測値との関係から、直接、応答関数を算出することを特徴とする磁気共鳴イメージング装置。
  5. 　請求項4に記載の磁気共鳴イメージング装置であって、
    　前記応答関数算出手段は、算出後の応答関数にフィルタをかけるフィルタ手段を備えることを特徴とする磁気共鳴イメージング装置。
  6. 　請求項5に記載の磁気共鳴イメージング装置であって、
    　前記フィルタ手段は、ウィナーフィルタ、メディアンフィルタを含むことを特徴とする磁気共鳴イメージング装置。
  7. 　請求項1に記載の磁気共鳴イメージング装置であって、
    前記傾斜磁場発生手段は、発生した傾斜磁場の出力波形に基き入力波形を制御する制御回路を備え、
    　前記出力波形を変化させる複数の要素は、前記制御回路による制御に伴う前記出力波形の歪みを含むことを特徴とする磁気共鳴イメージング装置。
  8. 　請求項1に記載の磁気共鳴イメージング装置であって、
    　前記応答関数算出手段は、前記静磁場発生手段が発生する静磁場の中心に対する位置が異なる複数の傾斜磁場について、それぞれ、応答関数を算出することを特徴とする磁気共鳴イメージング装置。
  9. 　請求項1に記載の磁気共鳴イメージング装置であって、
    　前記傾斜磁場算出手段が算出する傾斜磁場パルスは、前記撮像シーケンスにおけるスライス選択傾斜磁場であり、前記補償手段は、前記スライス選択傾斜磁場の計算値に基き、前記撮像シーケンスに設定された高周波パルスを補正する高周波パルス補正手段を備えたことを特徴とする磁気共鳴イメージング装置。
  10. 　請求項9に記載の磁気共鳴イメージング装置であって、
    　前記高周波パルス補正手段は、補正前後の高周波パルスの照射時間の変化を算出する照射時間計算手段を備え、前記照射時間計算手段が算出した照射時間の変化が所定の増加率を超えるときに、前記傾斜磁場算出手段の入力値を変更し、前記傾斜磁場算出手段が前記変更された入力値から推定した出力波形の推定値を用いて、前記高周波パルスの補正を行うことを特徴とする磁気共鳴イメージング装置。
  11. 　請求項1に記載の磁気共鳴イメージング装置であって、
    　前記傾斜磁場算出手段が算出する傾斜磁場パルスは、前記撮像シーケンスにおける読出し傾斜磁場であり、前記補償手段は、前記読出し選択傾斜磁場の計算値に基き、前記撮像シーケンスで得た核磁気共鳴信号が配置されるk空間座標を補正するk空間座標補正手段を備えたことを特徴とする磁気共鳴イメージング装置。
  12. 　請求項11に記載の磁気共鳴イメージング装置であって、
    　前記補正手段は、補正前後のk空間座標の大きさの変化を算出するk空間サイズ計算手段を備え、前記k空間サイズ計算手段が計算したk空間の大きさの変化が所定の縮小率を超えるときに、前記傾斜磁場算出手段の入力値を変更し、前記傾斜磁場算出手段が変更後の入力値から算出した出力波形の再計算値を用いて、前記k空間座標の補正を行うことを特徴とする磁気共鳴イメージング装置。
  13. 　磁気共鳴イメージング装置の傾斜磁場発生手段が発生する出力傾斜磁場波形を推定する方法であって、
    　前記傾斜磁場発生手段の入力値およびそれによって前記傾斜磁場発生手段が発生した出力傾斜磁場波形の実測値を入力するステップと、
    　前記出力傾斜磁場波形に影響を与える複数の因子について、それぞれ、応答関数モデルを設定するステップと、
    　前記入力した傾斜磁場波形の入力値および実測値を用いて、各応答関数モデルを適正化するステップと、
    　所望の傾斜磁場の入力値に対し、前記適正化した応答関数モデルを用いて、出力傾斜磁場波形の推定値を算出するステップと、
    　を含む傾斜磁場波形推定方法。
  14. 　請求項13に記載の傾斜磁場波形推定方法であって、
    　前記複数の因子は、傾斜磁場により発生する渦電流と、前記傾斜磁場発生手段に備えられた制御回路の出力とを含むことを特徴とする傾斜磁場波形推定方法。
  15. 　請求項14に記載の傾斜磁場波形推定方法であって、
    　前記実測値は、前記傾斜磁場発生手段の制御回路を無効にした時の第1の実測値と、前記制御回路を有効にした第2の実測値を含み、
    　前記応答関数モデルを適正化するステップは、前記第１の実測値および第２の実測値を用いて、制御回路による出力傾斜磁場波形の応答関数を求めることを特徴とする傾斜磁場波形推定方法。
  16. 　磁気共鳴イメージング装置における出力傾斜磁場波形を推定する方法であって、
    　前記傾斜磁場発生手段の少なくとも2種の入力値およびそれらによって前記傾斜磁場発生手段が発生した出力傾斜磁場波形の実測値を入力するステップと、
    　前記少なくとも２種の入力値およびそれに対応する実測値を用いて、前記出力傾斜磁場波形に影響を与える複数の因子の影響を含む前記傾斜磁場発生手段の応答関数を算出するステップと、
    　所望の傾斜磁場の入力値に対し、前記応答関数を用いて、出力傾斜磁場波形の推定値を算出するステップと、を含む傾斜磁場波形推定方法。
PCT/JP2012/066316 2011-06-30 2012-06-27 磁気共鳴イメージング装置および傾斜磁場波形推定方法 WO2013002233A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US14/128,051 US9664765B2 (en) 2011-06-30 2012-06-27 Magnetic resonance imaging apparatus and gradient magnetic field waveform estimation method
CN201280029381.5A CN103596496B (zh) 2011-06-30 2012-06-27 磁共振成像装置以及倾斜磁场波形推定方法
JP2013522881A JP5916240B2 (ja) 2011-06-30 2012-06-27 磁気共鳴イメージング装置および傾斜磁場波形推定方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011146605 2011-06-30
JP2011-146605 2011-06-30

Publications (1)

Publication Number Publication Date
WO2013002233A1 true WO2013002233A1 (ja) 2013-01-03

Family

ID=47424124

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/066316 WO2013002233A1 (ja) 2011-06-30 2012-06-27 磁気共鳴イメージング装置および傾斜磁場波形推定方法

Country Status (4)

Country Link
US (1) US9664765B2 (ja)
JP (1) JP5916240B2 (ja)
CN (1) CN103596496B (ja)
WO (1) WO2013002233A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111329479A (zh) * 2018-12-19 2020-06-26 西门子保健有限责任公司 用于提供mr图像的方法和磁共振成像设备、程序和介质
WO2021107085A1 (ja) * 2019-11-28 2021-06-03 国立大学法人神戸大学 外場応答分布可視化装置及び外場応答分布可視化方法
CN113671594A (zh) * 2021-06-21 2021-11-19 成都理工大学 基于bp神经网络的碳酸盐岩高频层序自动化识别方法

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5846450B2 (ja) * 2010-11-09 2016-01-20 株式会社日立メディコ 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法
JP6078465B2 (ja) * 2011-06-30 2017-02-08 株式会社日立製作所 磁気共鳴イメージング装置および高周波磁場決定方法
US9664765B2 (en) * 2011-06-30 2017-05-30 Hitachi, Ltd. Magnetic resonance imaging apparatus and gradient magnetic field waveform estimation method
KR101474757B1 (ko) * 2013-07-08 2014-12-19 삼성전자주식회사 자장 측정 방법 및 장치
US10761155B2 (en) * 2013-07-16 2020-09-01 Siemens Aktiengesellschaft Hybrid averaging method for specific absorption rate supervision
US10444277B2 (en) * 2014-01-03 2019-10-15 Koninklijke Philips N.V. Calculation of the probability of gradient coil amplifier failure using environment data
JP6291328B2 (ja) * 2014-04-09 2018-03-14 株式会社日立製作所 磁気共鳴イメージング装置および磁気共鳴イメージング装置に搭載されるパルスシーケンスの算出方法
JP6333078B2 (ja) * 2014-06-09 2018-05-30 株式会社日立製作所 磁気共鳴イメージング装置および傾斜磁場波形調整方法
US10712420B2 (en) 2015-05-29 2020-07-14 Mayo Foundation For Medical Education And Research Systems and methods for concomitant field correction in magnetic resonance imaging with asymmetric gradients
CN106483482B (zh) * 2015-08-25 2019-08-23 上海联影医疗科技有限公司 用于磁共振成像***的梯度涡流校正方法和装置
CN106656177B (zh) * 2016-09-07 2020-01-07 东软医疗***股份有限公司 提高dac输出波形精度的方法、装置及磁共振***
CN110050198B (zh) 2016-10-10 2021-12-31 皇家飞利浦有限公司 梯度脉冲响应函数映射
WO2018175807A1 (en) * 2017-03-22 2018-09-27 Viewray Technologies, Inc. Reduction of artifacts in magnetic resonance imaging by creating inhomogeneity in the magnetic field at gradient null position of an mri system
JP7123595B2 (ja) * 2017-03-29 2022-08-23 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および医用処理装置
EP3543722A1 (en) * 2018-03-20 2019-09-25 Koninklijke Philips N.V. Magnetic resonance imaging using corrected k-space trajectories calculated from current sensor data
EP3637099B1 (en) * 2018-10-08 2021-06-23 Ecole Polytechnique Federale De Lausanne (Epfl) Image reconstruction method based on a trained non-linear mapping
BR112021012557A2 (pt) 2018-12-28 2021-09-14 Hyperfine, Inc. Correção de histerese em formação de imagem por ressonância magnética
WO2021030722A1 (en) * 2019-08-15 2021-02-18 Hyperfine Research, Inc. Eddy current mitigation systems and methods
DE102020200786A1 (de) * 2020-01-23 2021-07-29 Siemens Healthcare Gmbh Korrektur von verzerrten diffusionsgewichteten Magnetresonanzbilddaten
US20220027685A1 (en) * 2020-07-24 2022-01-27 International Business Machines Corporation Automated generation of optimization model for system-wide plant optimization
EP4194876A1 (de) * 2021-12-09 2023-06-14 Siemens Healthcare GmbH Kompensationsschnittstelle für ein magnetresonanzbildgebungssystem

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0454614A (ja) * 1990-06-22 1992-02-21 Fuji Electric Co Ltd 定電流電源の制御装置
JPH1189817A (ja) * 1997-09-22 1999-04-06 Hitachi Medical Corp 磁気共鳴イメージング装置
WO2010047245A1 (ja) * 2008-10-20 2010-04-29 株式会社 日立メディコ 磁気共鳴イメージング装置及び方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02211124A (ja) * 1989-02-13 1990-08-22 Hitachi Medical Corp 磁気共鳴イメージング装置
JPH03173531A (ja) * 1989-12-01 1991-07-26 Hitachi Ltd 磁気共鳴イメージング装置の渦電流補償方法及び渦電流補償装置
JPH07171129A (ja) * 1993-12-21 1995-07-11 Hitachi Medical Corp 磁気共鳴イメージング装置
US6271666B1 (en) * 1999-05-21 2001-08-07 General Electric Company Filtering of oscillatory B0 behavior in NMR devices
US7286871B2 (en) * 2000-08-15 2007-10-23 The Regents Of The University Of California Method and apparatus for reducing contamination of an electrical signal
US6396268B1 (en) * 2000-10-02 2002-05-28 Ge Medical Systems Global Technology Company, Llc Magnetic resonance imaging device having magnetic field disturbance compensation
JP2002224083A (ja) * 2001-02-06 2002-08-13 Hitachi Medical Corp 磁気共鳴画像診断装置
DE10306017A1 (de) * 2003-02-13 2004-09-09 Siemens Ag Verfahren zum Ermitteln einer Kompensationseinstellung für ein Wirbelstromfeld
CN101711126B (zh) * 2007-06-14 2012-03-28 株式会社日立医药 磁共振成像装置及由倾斜磁场引起的误差修正方法
US8390377B2 (en) * 2008-12-27 2013-03-05 Sumitomo Electric Industries, Ltd. Distortion compensating circuit and radio base station
US9664765B2 (en) * 2011-06-30 2017-05-30 Hitachi, Ltd. Magnetic resonance imaging apparatus and gradient magnetic field waveform estimation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0454614A (ja) * 1990-06-22 1992-02-21 Fuji Electric Co Ltd 定電流電源の制御装置
JPH1189817A (ja) * 1997-09-22 1999-04-06 Hitachi Medical Corp 磁気共鳴イメージング装置
WO2010047245A1 (ja) * 2008-10-20 2010-04-29 株式会社 日立メディコ 磁気共鳴イメージング装置及び方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111329479A (zh) * 2018-12-19 2020-06-26 西门子保健有限责任公司 用于提供mr图像的方法和磁共振成像设备、程序和介质
CN111329479B (zh) * 2018-12-19 2023-09-29 西门子保健有限责任公司 用于提供mr图像的方法、磁共振成像设备和介质
WO2021107085A1 (ja) * 2019-11-28 2021-06-03 国立大学法人神戸大学 外場応答分布可視化装置及び外場応答分布可視化方法
JP7473114B2 (ja) 2019-11-28 2024-04-23 国立大学法人神戸大学 外場応答分布可視化装置及び外場応答分布可視化方法
CN113671594A (zh) * 2021-06-21 2021-11-19 成都理工大学 基于bp神经网络的碳酸盐岩高频层序自动化识别方法
CN113671594B (zh) * 2021-06-21 2023-06-20 成都理工大学 基于bp神经网络的碳酸盐岩高频层序自动化识别方法

Also Published As

Publication number Publication date
JPWO2013002233A1 (ja) 2015-02-23
US20140125333A1 (en) 2014-05-08
JP5916240B2 (ja) 2016-05-11
US9664765B2 (en) 2017-05-30
CN103596496B (zh) 2016-08-24
CN103596496A (zh) 2014-02-19

Similar Documents

Publication Publication Date Title
JP5916240B2 (ja) 磁気共鳴イメージング装置および傾斜磁場波形推定方法
JP5611661B2 (ja) 磁気共鳴イメージング装置
US10042027B2 (en) Magnetic resonance imaging apparatus and control device of a magnetic resonance imaging apparatus
JP5740307B2 (ja) 磁気共鳴イメージング装置及び傾斜磁場印加方法
WO2012124375A1 (ja) 磁気共鳴イメージング装置
JP5638393B2 (ja) 磁気共鳴イメージング装置及び方法
US10234530B2 (en) Magnetic resonance imaging apparatus
US11137468B2 (en) Method and magnetic resonance apparatus for acquiring a magnetic resonance dataset with correction of gradient impulse response functions
JP2010042275A (ja) 磁気共鳴イメージング装置及びその残留磁場の消磁方法
WO2012063654A1 (ja) 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法
US10162027B2 (en) Magnetic resonance imaging apparatus and irradiation magnetic field distribution measurement method
US11965945B2 (en) Magnetic resonance system and shimming method and imaging method thereof
US10261159B2 (en) Method and magnetic resonance apparatus for forming RF pulses for the compensation of inhomogeneities in the edge region of the magnetic resonance scanner
JP5650044B2 (ja) 磁気共鳴イメージング装置
JP2013017811A (ja) 磁気共鳴イメージング装置
US20230358834A1 (en) Reducing and correcting magnetic field gradient deviations
US11681002B2 (en) Method for recording a magnetic resonance image data set, data carrier, computer-program product, and magnetic resonance system
US20230358839A1 (en) Reducing and correcting magnetic field gradient deviations
JP7253433B2 (ja) 磁気共鳴イメージング装置及び誤差磁場の補償方法
JP2012239561A (ja) 磁気共鳴イメージング装置及びＤｕａｌＳｌｉｃｅ計測方法
CN112630712A (zh) 一种简便的核磁共振梯度波形失真预矫正方法
JP4699744B2 (ja) Ｍｒｉ装置およびｍｒｉ装置の画質改善方法
JP2022130030A (ja) 磁気共鳴イメージング装置、画像補正方法、および静磁場不均一補正方法
JP4424968B2 (ja) 磁気共鳴イメージング装置及び方法
JP2014200590A (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12803944

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2013522881

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14128051

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12803944

Country of ref document: EP

Kind code of ref document: A1