WO2013105384A1 - 磁気共鳴イメージング装置及び領域撮像方法 - Google Patents
磁気共鳴イメージング装置及び領域撮像方法 Download PDFInfo
- Publication number
- WO2013105384A1 WO2013105384A1 PCT/JP2012/082139 JP2012082139W WO2013105384A1 WO 2013105384 A1 WO2013105384 A1 WO 2013105384A1 JP 2012082139 W JP2012082139 W JP 2012082139W WO 2013105384 A1 WO2013105384 A1 WO 2013105384A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- magnetic field
- static magnetic
- phase
- excitation frequency
- shimming
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 99
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 57
- 230000005284 excitation Effects 0.000 claims abstract description 222
- 230000003068 static effect Effects 0.000 claims abstract description 192
- 238000012545 processing Methods 0.000 claims description 103
- 238000004364 calculation method Methods 0.000 claims description 65
- 238000005259 measurement Methods 0.000 claims description 58
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 31
- 230000036278 prepulse Effects 0.000 claims description 26
- 238000009828 non-uniform distribution Methods 0.000 claims description 9
- 230000006866 deterioration Effects 0.000 abstract description 4
- 238000000034 method Methods 0.000 description 44
- 238000009826 distribution Methods 0.000 description 42
- 230000006870 function Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 14
- 230000005540 biological transmission Effects 0.000 description 9
- 238000005481 NMR spectroscopy Methods 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 5
- 230000015556 catabolic process Effects 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 125000004429 atom Chemical group 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/387—Compensation of inhomogeneities
- G01R33/3875—Compensation of inhomogeneities using correction coil assemblies, e.g. active shimming
-
- 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
-
- 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/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
Definitions
- the present invention relates to a magnetic resonance imaging apparatus (hereinafter referred to as “MRI apparatus”) using a nuclear magnetic resonance phenomenon, and more particularly to a technique for imaging a selected region in response to static magnetic field inhomogeneity.
- MRI apparatus magnetic resonance imaging apparatus
- 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.
- the NMR signal is given different phase encoding depending on the gradient magnetic field, frequency-encoded, and measured as time series data.
- the measured NMR signal is reconstructed into an image by two-dimensional or three-dimensional Fourier transform.
- Ultra-high magnetic field MRI equipment of 3T or more increases the static magnetic field non-uniformity, which is the spatial non-uniformity of the static magnetic field.
- SSFP steady-state free precession
- EPI Echo-Planar-Imaging sequences that measure echo signals are sensitive to static magnetic field inhomogeneities. Therefore, an artifact generated based on the static magnetic field inhomogeneity becomes prominent in an image obtained by using these pulse sequences with an ultrahigh magnetic field MRI apparatus of 3T or more.
- the excitation profile may be deteriorated by static magnetic field inhomogeneity.
- Bo shimming As a standard method for reducing the static magnetic field inhomogeneity, Bo shimming is used in which an appropriate current is supplied to the shim coil to adjust the static magnetic field uniformity of the imaging target region.
- a method called local Bo shimming that generates a compensation magnetic field that reduces local static magnetic field inhomogeneity is used.
- the static magnetic field inhomogeneity unique to that apparatus is measured, and the static magnetic field inhomogeneity of the imaging space including the imaging target is measured during imaging. Then, a compensation magnetic field for reducing the static magnetic field inhomogeneity in the imaging target region is generated in the shim coil from both data.
- SSPF sequence band artifacts cannot be resolved by local Bo shimming alone. Therefore, in the SSFP sequence, in addition to local Bo shimming, a method is employed in which the excitation frequency of the RF pulse is slightly shifted so that band artifacts are not superimposed on the heart region. To determine the excitation frequency, create a static magnetic field inhomogeneity estimation map after local Bo shimming from static magnetic field inhomogeneity, and calculate the excitation frequency to reduce the average static magnetic field inhomogeneity in the short axis direction of the heart. The method is used.
- Patent Document 1 a method in which a plurality of images with slightly shifted excitation frequencies are taken, and an operator visually obtains an appropriate excitation frequency so that artifacts are not superimposed on the short axis of the heart.
- the position of the band artifact is unintentionally shifted from the appropriate position. End up. Further, if the position of the band artifact is not appropriate even after adjusting the local Bo shimming and the excitation frequency, it is difficult for the operator to adjust.
- the present invention has been made in view of the above problems, and an MRI apparatus capable of reducing image quality degradation due to non-uniform static magnetic fields in response to imaging conditions without increasing the burden on the operator and An object is to provide a region imaging method.
- the present invention obtains a shimming current so that the static magnetic field inhomogeneity in the selected region is not reduced, shims the selected region using the obtained local Bo shimming current, and generates a static magnetic field in the selected region.
- An irradiation phase increment amount (RF-Phase) or an adjusted excitation frequency (f0 ′) which is an excitation frequency of the RF pulse for exciting the selected region in a state where non-uniformity is reduced is obtained.
- the irradiation phase increment amount (RF-Phase) and the adjusted excitation frequency (f0 ′), which is the excitation frequency correspond to each other.
- imaging is performed with an RF pulse using at least one of the increment amount (RF-Phase) of the irradiation phase and the adjusted excitation frequency (f0 ′) that is the excitation frequency.
- the MRI apparatus and the region imaging method of the present invention it is possible to reduce image quality deterioration due to non-uniform static magnetic field corresponding to imaging conditions without increasing the burden on the operator.
- FIG. 1 A flowchart showing a processing flow for selecting a control method according to an imaging parameter
- (a) Figure is a functional block diagram of the first embodiment
- (b) Figure is a functional block diagram of the second embodiment.
- FIG. 1 is a block diagram showing the overall configuration of an embodiment of an MRI apparatus according to the present invention.
- This MRI apparatus uses a NMR phenomenon to obtain a tomographic image of a subject 101.
- a static magnetic field generating magnet 102, a gradient magnetic field coil 103, a gradient magnetic field power supply 109, and an RF transmission coil 104 and RF transmitter 110, RF receiver coil 105 and signal processor 107, measurement control unit 111, overall control unit 112, display / operation unit 118, and top plate on which the subject 101 is mounted generates a static magnetic field.
- a bed 106 to be taken in and out of the magnet 102.
- the static magnetic field generating magnet 102 generates a uniform static magnetic field in the direction perpendicular to the body axis of the subject 101 in the vertical magnetic field method and in the body axis direction in the horizontal magnetic field method.
- a permanent magnet type, normal conducting type or superconducting type static magnetic field generating source is arranged around the.
- the gradient magnetic field coil 103 is a coil wound in the three-axis directions of X, Y, and Z that are the real space coordinate system (stationary coordinate system) of the MRI apparatus, and each gradient magnetic field coil is a gradient magnetic field that drives it.
- a current is supplied to the power source 109.
- the gradient magnetic field power supply 109 of each gradient coil is driven according to a command from the measurement control unit 111 described later, and supplies a current to each gradient coil.
- gradient magnetic fields Gx, Gy, and Gz are generated in the three-axis directions of X, Y, and Z.
- a slice gradient magnetic field pulse (Gs) is applied in a direction orthogonal to the slice plane (imaging cross section) to set a slice plane for the subject 101, orthogonal to the slice plane and orthogonal to each other.
- Phase encoding gradient magnetic field pulse (Gp) and frequency encoding (leadout) gradient magnetic field pulse (Gf) are applied in the remaining two directions, and position information in each direction is encoded in the NMR signal (echo signal). .
- the gradient coil 103 is also provided with a shim coil that generates a compensation magnetic field that is supplied with shimming current to reduce non-uniformity of the static magnetic field.
- Each shim coil has a component coil that generates a compensation magnetic field of each order.
- a secondary component x ⁇ 2, y ⁇ 2, xy, yz, zx, (x ⁇ 2-y ⁇ 2) component, etc.
- a further higher-order component may be included.
- the 0th-order (Bo component) component is compensated by the excitation frequency f0 of the RF pulse, and the 1st-order component is also used as the gradient magnetic field coil.
- the RF transmission coil 104 is a coil that irradiates the subject 101 with an irradiation RF magnetic field pulse (hereinafter abbreviated as an RF pulse), and is connected to the RF transmission unit 110 and supplied with a high-frequency pulse current.
- an NMR phenomenon is induced in the spins of atoms constituting the living tissue of the subject 101.
- the RF transmission unit 110 is driven in accordance with a command from the measurement control unit 111 (to be described later), amplitude-modulates and amplifies the high-frequency pulse, and then the RF transmission unit 104 is placed near the subject 101 after being amplified. By supplying, the subject 101 is irradiated with the RF pulse.
- the RF receiving coil 105 is a coil that receives an echo signal emitted by the NMR phenomenon of spin that constitutes the living tissue of the subject 101, and is connected to the signal processing unit 107 so that the received echo signal is sent to the signal processing unit 107. Sent.
- the signal processing unit 107 performs processing for detecting an echo signal received by the RF receiving coil 105. Specifically, in accordance with a command from the measurement control unit 111 described later, the signal processing unit 107 amplifies the received echo signal and divides it into two orthogonal signals by quadrature detection, For example, 128, 256, 512, etc.) are sampled, and each sampling signal is A / D converted into a digital quantity. Therefore, the echo signal is obtained as time-series digital data (hereinafter referred to as echo data) composed of a predetermined number of sampling data. Then, the signal processing unit 107 performs various processes on the echo data, and sends the processed echo data to the measurement control unit 111.
- echo data time-series digital data
- the measurement control unit 111 mainly transmits various commands for collecting echo data necessary for reconstruction of the tomographic image of the subject 101 to the gradient magnetic field power source 109, the RF transmission unit 110, and the signal processing unit 107. And a control unit for controlling them. Specifically, the measurement control unit 111 operates under the control of the overall control unit 112 described later, and controls the gradient magnetic field power source 109, the RF transmission unit 110, and the signal processing unit 107 based on a predetermined sequence of control data. Then, it is necessary to reconstruct the image of the imaging region of the subject 101 by repeatedly performing the irradiation of the RF pulse and the application of the gradient magnetic field pulse to the subject 101 and the detection of the echo signal from the subject 101. Control the collection of accurate echo data.
- the application amount of the phase encoding gradient magnetic field is changed in the case of two-dimensional imaging, and the application amount of the slice encoding gradient magnetic field is further changed in the case of three-dimensional imaging.
- Values such as 128, 256, and 512 are normally selected as the number of phase encodings, and values such as 16, 32, and 64 are normally selected as the number of slice encodings.
- the overall control unit 112 controls the measurement control unit 111 and controls various data processing and processing result display and storage, and includes an arithmetic processing unit (CPU) 114, a memory 113, and a magnetic disk. And the like, and a network IF 116 that interfaces with an external network. Further, an external storage unit 117 such as an optical disk may be connected to the overall control unit 112. Specifically, the measurement control unit 111 is controlled to execute the collection of echo data, and when the echo data is input from the measurement control unit 111, the arithmetic processing unit 114 converts the encoded information applied to the echo data. Based on this, it is stored in an area corresponding to the k space in the memory 113.
- the statement that the echo data is arranged in the k space means that the echo data is stored in an area corresponding to the k space in the memory 113.
- a group of echo data stored in an area corresponding to the k space in the memory 113 is also referred to as k space data.
- the arithmetic processing unit 114 performs processing such as signal processing and image reconstruction by Fourier transform on the k-space data, and displays the resulting image of the subject 101 on the display / operation unit 118 described later.
- the data is recorded in the internal storage unit 115 or the external storage unit 117, or transferred to an external device via the network IF 116.
- the display / operation unit 118 includes a display unit for displaying the reconstructed image of the subject 101, a trackball or a mouse and a keyboard for inputting various control information of the MRI apparatus and control information for processing performed by the overall control unit 112. Etc., and an operation unit.
- the operation unit is disposed in the vicinity of the display unit, and an operator interactively controls various processes of the MRI apparatus through the operation unit while looking at the display unit.
- the radionuclide to be imaged by the MRI apparatus is a hydrogen nucleus (proton) which is the main constituent material of the subject, as is widely used in clinical practice.
- proton the main constituent material of the subject
- the form or function of the human head, abdomen, limbs, etc. is imaged two-dimensionally or three-dimensionally.
- the area imaging method according to the present invention includes the following steps. That is, A static magnetic field inhomogeneous distribution representing a spatial distribution of the static magnetic field inhomogeneity in a region including a region selected through the ROI setting (hereinafter abbreviated as a selected region) is measured.
- the local Bo shimming current (C) is obtained so that the static magnetic field inhomogeneity in the selected region is not reduced. Then, shimming the selected region using the determined local Bo shimming current (C), the adjusted excitation frequency (f0 ′) for exciting the selected region in a state where the static magnetic field inhomogeneity of the selected region is reduced Ask for.
- the obtained local Bo shimming current (C) is passed through the shim coil, and the RF pulse of the adjusted excitation frequency (f0 ′) is applied in a state where static magnetic field inhomogeneity in the excitation region is reduced.
- the increment amount (RF-Phase) may be directly obtained from the frequency corresponding to the static magnetic field inhomogeneity without being obtained from the adjusted excitation frequency (f0 ′).
- phase control refers to obtaining the increment amount (RF-Phase) of the irradiation phase of the RF pulse used in the pulse sequence based on the static magnetic field inhomogeneity, and supplying a shimming current to the shimming unit to In the state where the uniformity is reduced, the measurement of the echo signal is controlled by changing the irradiation phase of the RF pulse by the increment amount (RF-Phase).
- frequency control is the calculation of shimming current based on the static magnetic field inhomogeneity, and the post-adjustment excitation of RF pulses in the pulse sequence in a state where the shimming current is supplied to the shimming part and the static magnetic field inhomogeneity is reduced
- shimming current and post-adjustment excitation frequency are calculated for each selection area and adjusted to the RF pulse of the pulse sequence according to how the multiple selection areas overlap. Controlling the measurement of the echo signal using the post-excitation frequency.
- Examples 1 and 2 are examples of phase control
- Examples 3 to 6 are examples of frequency control
- Examples 7 and 8 are examples of using both phase control and frequency control.
- the static magnetic field inhomogeneous distribution for example, a dual-echo-Gradient-Echo sequence with an excitation frequency f0 can be used. If the phase image obtained from the first echo signal is ⁇ 1, the phase image obtained from the second echo signal is ⁇ 2, and the time interval between the peak of the first echo signal and the peak of the second echo signal is ⁇ , the static magnetic field is inhomogeneous
- the distribution is as shown in equation (1). In this specification, the static magnetic field inhomogeneity is expressed by frequency.
- Gx, Gy, and Gz are primary compensation magnetic field components generated by primary shim coils x, y, and z, respectively. Since the primary shim coil is usually used also as the gradient magnetic field coil, the current for generating the primary compensation magnetic field is supplied to the gradient magnetic field coil.
- Gxy is a secondary compensation magnetic field component generated by the secondary shim coil xy. Although only xy is described as a second or higher order term, in general, higher order terms such as x ⁇ 2, y ⁇ 2, yz, zx, (x ⁇ 2-y ⁇ 2) may be included.
- the post-adjustment excitation frequency (f0 ′) is determined by equation (3).
- f0 ' f0 + ⁇ f0 + ⁇ f0> (3)
- f0 is the excitation frequency of dual echo Gradient Echo used to measure the static magnetic field inhomogeneous distribution
- ⁇ f0 is the frequency offset obtained by equation (2)
- ⁇ f0> is the static magnetic field non-uniformity after shimming. This is a frequency offset obtained based on the uniform estimation map.
- the adjusted excitation frequency (f0 ') of the RF pulse that excites the selected region is calculated based on the state where the static magnetic field inhomogeneity is reduced, and is used when measuring the static magnetic field inhomogeneity distribution.
- the frequency corresponding to the static magnetic field inhomogeneity ( ⁇ f0 + ⁇ f0>) is added to the excitation frequency of the RF pulse.
- the post-shimming static magnetic field inhomogeneity estimation map is obtained by adding the magnetic field generated when the shim coil current values [Gx, Gy, Gz, Gxy] are passed to ⁇ B0.
- ⁇ f0> for example, the average of the maximum and minimum values of the post-shimming static magnetic field inhomogeneity estimation map in the selected region, the arithmetic average of the static magnetic field inhomogeneity estimation map in the selected region, or the selected region The average of the histogram of the static magnetic field inhomogeneity estimation map is calculated.
- the SSFP sequence irradiates RF pulses with a constant flip angle continuously for a certain repetition time TR. At that time, the irradiation phase is increased by 180 degrees.
- the excitation frequency of the SSFP sequence is f0
- the spin excited at the adjusted excitation frequency (f0 ') derived from equation (3) rotates by (f0'-f0) * TR * 360 [deg.] Between TRs. Therefore, the behavior is the same as when RF-Phase is 180+ (f0'-f0) * TR * 360.
- % represents a remainder operator
- the operation of A% B means a remainder obtained by dividing A by B. That is, RF-Phase is a phase amount that gives a phase rotation corresponding to the frequency corresponding to non-uniform static magnetic field ( ⁇ f0 + ⁇ f0>) to the irradiation phase of the RF pulse.
- RF-Phase can also be obtained through calculation of the adjusted excitation frequency (f0 '), or directly from the frequency component ( ⁇ f0 + ⁇ f0>) corresponding to the static magnetic field inhomogeneity.
- f0 which is the actual excitation frequency of the SSFP sequence, measures the spectrum of the imaging slice cross section as pre-measurement (pre-scan), and may be automatically set to the peak frequency of the water spectrum.
- pre-scan pre-measurement
- the irradiation phase is changed in increments (RF-Phase) determined by equation (4), and the actual excitation frequency Even if is f0, the excitation frequency has the same effect as f0 ′. Furthermore, since the position of the band artifact is controlled by the RF-Phase, even if the actual excitation frequency (f0) is changed, the position of the band artifact is fixed or moved to the desired position by controlling the RF-Phase. It becomes possible to do.
- Imaging parameters related to the selection of the control method include the number of selected areas and the sequence type of the imaging sequence.
- the sequence type means a sequence type of the main sequence.
- step 200 the arithmetic processing unit 114 determines whether the selected region set by the operator is a single region or a plurality of regions.
- the case of a single region is a case where the main sequence and its pre-pulse sequence or other main sequence both excite only the same region.
- the case of a plurality of regions is a case where the main sequence and its pre-pulse sequence or another main sequence excite different regions. If the area is a single area 201, the process proceeds to step 202. If the area is a plurality of areas 205, the process proceeds to step 206.
- the arithmetic processing unit 114 determines whether the sequence type is SSFP or other.
- the imaging sequence is generated by calculating the shimming current value and the excitation frequency using the known technique described in Non-Patent Document 1 (210).
- the arithmetic processing unit 114 determines whether or not the pre-pulse sequence is used together.
- the imaging sequence is generated by calculating the local Bo shimming current value and the RF-Phase of the SSFP sequence using Example 1 or 2 below (211).
- Example 7 When using the pre-pulse sequence together, use Example 7 to calculate the local Bo shimming current value and adjusted excitation frequency of the pre-pulse sequence, calculate the local Bo shimming current value and RF-Phase of the SSFP sequence, and image A sequence is generated (212).
- the local Bo shimming current value and adjusted excitation frequency of the pre-pulse sequence are calculated according to how the regions to be excited overlap, using Example 8 below. Then, the local Bo shimming current value and RF-Phase of the SSFP sequence are calculated to generate an imaging sequence (213).
- Example 1 A first embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- the local Bo shimming current flowing through the shim coil and the RF pulse implementation amount (RF-Phase) are calculated according to the static magnetic field inhomogeneity, and the calculated local Bo shimming current and RF-Phase are calculated. Use to image.
- the first embodiment will be described in detail with reference to FIGS.
- the arithmetic processing unit 114 includes a region selection unit 301, a static magnetic field nonuniform distribution measurement unit 302, a local Bo shimming current calculation unit 303, an excitation frequency calculation unit 304, and an RF-Phase calculation unit 305. And a setting unit 306.
- the region selection unit 301 displays a positioning image on the operation / display unit 118 and accepts an ROI setting input for selecting the region of the operator. Then, the selection area selected by setting the ROI is set as the imaging area. Then, the position information of the selected region is notified to the static magnetic field nonuniform distribution measuring unit 302 and the local Bo shimming current calculating unit 303.
- the static magnetic field inhomogeneous distribution measuring unit 302 generates control data of a measurement sequence for measuring the static magnetic field inhomogeneous distribution in the region including the selected region and notifies the measurement control unit 111 of the measurement sequence control data. Control the measurement of static magnetic field inhomogeneous distribution. Then, the measurement control unit 111 activates the static magnetic field inhomogeneous distribution measurement sequence, and notifies the local Bo shimming current calculation unit 303 of the measured echo data. As described above, the dual echo Gradient Echo sequence can be used for the static magnetic field inhomogeneous distribution measurement sequence.
- the local Bo shimming current calculation unit 303 calculates a local Bo shimming current (C) such that the static magnetic field inhomogeneity in the selected region is reduced. Therefore, first, in order to extract the static magnetic field inhomogeneous distribution in the selected region, mask processing is performed on the static magnetic field inhomogeneous distribution in the region including the selected region. Next, the local Bo shimming current (C) is calculated so that the static magnetic field inhomogeneity in the selected region is reduced. The details of the calculation of the local Bo shimming current (C) are as described above. Then, the calculated local Bo shimming current (C) is notified to the excitation frequency calculation unit 304 and the setting unit 306.
- the excitation frequency calculation unit 304 uses the local Bo shimming current (C) to calculate an adjusted excitation frequency (f0 ′) for exciting the selected region so that the static magnetic field inhomogeneity in the selected region is reduced.
- C Bo shimming current
- the RF-Phase calculation unit 305 calculates the RF-Phase as the phase rotation amount corresponding to the adjusted excitation frequency (f0 ′) calculated by the excitation frequency calculation unit 304.
- RF-Phase is calculated as a phase rotation amount corresponding to the frequency ( ⁇ f0 + ⁇ f0>) corresponding to the static magnetic field inhomogeneity. The details of the RF-Phase calculation are as described above. Then, the RF-Phase is notified to the setting unit 306.
- the setting unit 306 includes the local Bo shimming current (C) calculated by the local Bo shimming current calculation unit 303, the RF-Phase calculated by the RF-Phase calculation unit 305, and the actual excitation frequency (f0).
- the local Bo shimming current, RF-Phase, and excitation frequency in the selected area are set. Then, using the set local Bo shimming current (C), RF-Phase, and excitation frequency (f0), control data of the imaging sequence is generated and notified to the measurement control unit 111.
- step 400 the region selection unit 301 displays a positioning image on the operation / display unit 118, and accepts an ROI setting input for region selection by the operator. Then, the static magnetic field inhomogeneity measurement unit 302 is notified of the position information of the selected region selected by the ROI.
- step 401 the static magnetic field inhomogeneity measuring unit 302 measures the static magnetic field inhomogeneous distribution in the area including the selected area. Details are as described above.
- step 402 the local Bo shimming current calculation unit 303 performs a mask process on the static magnetic field inhomogeneity distribution measured in step 401 so as to select only the selected region. Then, a local Bo shimming current (C) that reduces the static magnetic field inhomogeneity in the selected region is calculated. Details are as described above.
- the excitation frequency calculation unit 304 calculates the adjusted excitation frequency (f0 ′) so that the static magnetic field inhomogeneity in the selected region is reduced. Details are as described above.
- step 404 the RF-Phase calculation unit 305 calculates an RF-Phase corresponding to the adjusted excitation frequency (f0 ′) calculated in step 403. Details are as described above.
- step 405 the setting unit 306 determines the local Bo shimming current (C) calculated in step 402 and the RF calculated in step 404 as the local Bo shimming current, RF-Phase, and excitation frequency when exciting the selected region.
- -Phase and the actual excitation frequency (f0) used in the calculation of step 404 are set, respectively, and control data of an imaging sequence described later is generated and notified to the measurement control unit 111.
- step 406 the measurement control unit 111 performs phase control to change the irradiation phase of the RF pulse by increments using the local Bo shimming current (C), RF-Phase and excitation frequency (f0) set in step 405. Apply and control imaging.
- C Bo shimming current
- f0 excitation frequency
- RF-Phase may be calculated as a phase rotation amount corresponding to the frequency ( ⁇ f0 + ⁇ f0>) corresponding to the static magnetic field inhomogeneity derived from the above. In this case, the calculation of the adjusted excitation frequency (f0 ′) in step 403 is omitted.
- a plurality of selected areas are set by a plurality of ROI settings in step 400, and a static magnetic field inhomogeneous distribution of an area including the plurality of selected areas is obtained in 401. Steps 402 to 406 are repeated for each region. Alternatively, Example 8 described later is executed. The above is the description of the processing flow of the first embodiment.
- the MRI apparatus and the region imaging method of the first embodiment obtain the increment amount (RF-Phase) of the irradiation phase of the RF pulse used for the pulse sequence based on the static magnetic field inhomogeneity, and the shimming unit In the state in which the shimming current is supplied and the static magnetic field inhomogeneity is reduced, the phase control for changing the irradiation phase of the RF pulse by the increment amount is applied to control the echo signal measurement. As a result, it is possible to reduce image quality degradation caused by the static magnetic field non-uniform distribution. Further, since the position of the band artifact is controlled by RF-Phase, even if the operator changes the excitation frequency (f0), it is possible to fix the band artifact at the same position without changing the position of the band artifact.
- Example 2 Next, a second embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- an appropriate RF-Phase is automatically calculated.
- the static magnetic field inhomogeneity in the ROI is 1 / TR [Hz] or more, band artifacts are mixed in the ROI.
- the operator can easily adjust the position of the band artifact.
- the outline of the processing flow for setting the local Bo shimming current and the adjusted excitation frequency (f0 ′) is the same as that in the first embodiment, but the increment amount is obtained by adding the phase set by the operator.
- the position of the band artifact can be adjusted.
- GUI graphical user interface
- this GUI includes an RF-Phase display unit 501 that displays the value of RF-Phase calculated by the equation (4) described in the first embodiment, and an RF frequency by the operator.
- An adjustment value input unit 502 to which an adjustment value ⁇ of -Phase is input.
- a display that clearly indicates that the automatically adjusted value is used on the RF-Phase display unit 501 (e.g., ⁇ Auto '' )
- the operator manually inputs and sets the RF-Phase adjustment value ⁇ in the adjustment value input unit 502 of the GUI.
- a local Bo shimming current (C) for the selected region and an adjusted excitation frequency (f0 ′) that reduces band artifacts in the selected region are calculated.
- the calculation method is the same as in the first embodiment.
- each variable is the same as the equation (4).
- the frequency component corresponding to the static magnetic field inhomogeneity ( ⁇ f0 + ⁇ f0>) is used without using the adjusted excitation frequency (f0 ′) as in the above-described first embodiment.
- the adjustment value ⁇ set by the operator may be added to the phase rotation amount corresponding to) to obtain the RF-Phase.
- the subsequent processing is the same as in the first embodiment.
- the RF-Phase calculation unit 305 calculates the RF-Phase using the adjustment value ⁇ of the RF-Phase. Specifically, the adjusted excitation frequency (f0 ′) calculated by the excitation frequency calculation unit 304 and the RF ⁇ set and input by the operator via the adjustment value input unit 502 displayed on the display / operation unit 118 Based on the phase adjustment value ⁇ , RF-Phase is calculated using equation (5). In other words, the RF-Phase is calculated based on the RF-Phase corresponding to the adjusted excitation frequency (f0 ′) calculated by the excitation frequency calculation unit 304 and the RF-Phase adjustment value ⁇ .
- the equation (5) is To calculate the RF-Phase. Then, the calculated RF-Phase is notified to the setting unit 306.
- the RF-Phase calculation unit 305 displays the GUI shown in FIG. 5 on the display / operation unit 118, and accepts the input setting of the adjustment value ⁇ of the RF-Phase by the operator. Then, the RF-Phase is calculated based on the adjusted excitation frequency (f0 ′) calculated by the excitation frequency calculation unit 304 and the adjustment value ⁇ of the RF-Phase. Alternatively, the RF-Phase is calculated based on the phase rotation amount corresponding to the frequency ( ⁇ f0 + ⁇ f0>) corresponding to the static magnetic field inhomogeneity and the RF-Phase adjustment value ⁇ set and input by the operator. . Details are as described above.
- a plurality of selected areas are set by a plurality of ROI settings in step 400, and a static magnetic field defect of an area including the plurality of selected areas is set in 401.
- a uniform distribution is obtained, and steps 402 to 406 are repeated for each selected region.
- Example 8 described later is executed. The above is the description of the processing flow of the second embodiment.
- the operator sets the adjustment value of the RF-Phase with respect to the increment amount (RF-Phase) calculated according to the static magnetic field inhomogeneous distribution.
- the operator can easily adjust the position of the band artifact.
- the position of the band artifact is controlled by RF-Phase, so even if the operator changes the excitation frequency (f0), the position of the band artifact is not changed. It can be fixed.
- Example 3 A third embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- a local Bo shimming current flowing through the shim coil and an adjusted excitation frequency are calculated for each region, and images are sequentially captured for each region. Therefore, in the third embodiment, when the selection area 1 and the selection area 2 in FIG. 6 (a) do not have a common area, or when the selection area 1 and the selection area 2 in FIG. 6 (b) have a common area. It is suitable for imaging.
- the third embodiment is suitable when at least part of the two regions do not overlap each other.
- the common area where the two areas overlap is imaged twice.
- the third embodiment will be described in detail with reference to FIGS.
- each function of the arithmetic processing unit 114 for realizing the area imaging method of the third embodiment will be described based on the functional block diagram shown in FIG. Compared with the block diagram shown in FIG. 3 (a) described in the first embodiment, the functional block of the arithmetic processing unit 114 according to the third embodiment has no RF-Phase calculation unit 305 and includes a determination unit 705. Have.
- the functional block of the arithmetic processing unit 114 according to the third embodiment has no RF-Phase calculation unit 305 and includes a determination unit 705. Have.
- only parts different from the functional block diagram of the first embodiment will be described, and description of the same functional contents will be omitted.
- the area selection unit 301 is the same as that in the first embodiment.
- the static magnetic field nonuniform distribution measuring unit 302 is the same as that of the first embodiment.
- the local Bo shimming current calculation unit 303 calculates a local Bo shimming current (C) that reduces the static magnetic field inhomogeneity in the selected region. Details of the calculation of the local Bo shimming current are the same as in the first embodiment.
- the excitation frequency calculation unit 304 sets the adjusted excitation frequency (f0 ′) for exciting the selected region so that the static magnetic field inhomogeneity in the selected region is reduced. Details of the calculation of the adjusted excitation frequency (f0 ′) are the same as in the first embodiment.
- the determination unit 705 uses the local Bo shimming current (C) calculated for the selected region and the adjusted excitation frequency (f0 ′) to generate a static magnetic field inhomogeneous distribution (that is, a post-shimming static magnetic field inhomogeneity estimation map). It is determined whether or not the variance in the selected area satisfies a predetermined criterion. As a determination criterion, for the static magnetic field inhomogeneous distribution in the selected region in the post-shimming static magnetic field inhomogeneity estimation map, the variance is below a certain threshold (for example, 100 Hz), or the difference between the maximum and minimum is below the threshold. And so on. Then, the local Bo shimming current (C) that satisfies the determination criterion and the adjusted excitation frequency (f0 ′) are notified to the setting unit 306.
- a certain threshold for example, 100 Hz
- the setting unit 306 determines the local Bo shimming current (C) and the adjusted excitation frequency (f0 ′) that the static magnetic field inhomogeneous distribution in the selected region satisfies a predetermined criterion, and the local Bo shimming current and the excitation frequency in the selected region. Set as. Then, using the set local Bo shimming current (C) and the adjusted excitation frequency (f0 ′), control data of the imaging sequence is generated and notified to the measurement control unit 111.
- step 800 the region selection unit 301 displays a positioning image on the operation / display unit 118, and accepts setting input of each ROI for selection of the selection region 1 and the selection region 2 by the operator. Then, the position information of the selected selection region 1 and selection region 2 is notified to the static magnetic field inhomogeneity measurement unit 302.
- step 801 the static magnetic field inhomogeneity measuring unit 302 measures the static magnetic field inhomogeneous distribution in a region including the selection region 1 and the selection region 2. Details are as described above.
- step 802 the local Bo shimming current calculation unit 303 performs mask processing on the static magnetic field inhomogeneity distribution measured in step 801 so as to select only the selected region 1. Then, a local Bo shimming current (C 1 ) is calculated so that the static magnetic field inhomogeneity of the selected selection region 1 is reduced. Details are as described above.
- step 803 the excitation frequency calculation unit 304 sets the adjusted excitation frequency (f0 1 ′) of the selection region 1 so that the static magnetic field inhomogeneity of the selection region 1 is reduced. Details are as described above.
- step 804 the local Bo shimming current calculation unit 303 performs a mask process on the static magnetic field inhomogeneity distribution measured in step 801 so as to select only the selected region 2. Then, a local Bo shimming current (C 2 ) that reduces the static magnetic field inhomogeneity in the selected selection region 2 is calculated. Details are as described above.
- step 805 the excitation frequency calculation unit 304 sets the adjusted excitation frequency (f0 2 ′) of the selection region 2 so that the static magnetic field inhomogeneity of the selection region 2 is reduced. Details are as described above.
- step 806 the determination unit 305 determines the local Bo shimming current (C 1 , C 2 ) calculated in steps 802 and 804 and the adjusted excitation frequency calculated in steps 803 and 805 (selected regions 1 and 2 ).
- f0 1 ', f0 2 ') is used to determine whether or not the static magnetic field inhomogeneous distribution (that is, the post-shimming static magnetic field inhomogeneity estimation map) in the selected region 1 and selected region 2 satisfies a predetermined criterion. .
- the evaluation criteria are as described above. When the evaluation criteria are satisfied (Yes), the process proceeds to Step 807, and when not satisfied (No), the process proceeds to Step 808.
- the setting unit 306 uses the local Bo shimming current (C 1 ) calculated in step 802 and the adjusted excitation calculated in step 803 as the local Bo shimming current and excitation frequency when exciting the selected region 1.
- Each of the frequencies (f0 1 ′) is set, control data for an imaging sequence to be described later is generated, and notified to the measurement control unit 111.
- the local Bo shimming current and excitation frequency when exciting the selected region 2 the local Bo shimming current (C 2 ) calculated in step 804 and the adjusted excitation frequency (f0 2 ′) calculated in step 805 are respectively set, and control data of an imaging sequence to be described later is generated and notified to the measurement control unit 111.
- step 808 the determination unit 305 gives a warning display or a warning sound to the operation / display unit 118 that the imaging condition needs to be changed because the determination condition is not satisfied.
- step 809 the operator resets the imaging parameters and the like so that the static magnetic field inhomogeneous distribution in the selection area 1 and the selection area 2 satisfies a predetermined determination criterion. For example, readjustment of selected areas 1 and 2 is performed, TR is extended in the case of SSFP sequence, and excitation position is readjusted in the case of pencil beam. Thereafter, the process proceeds to step 801, where the static magnetic field inhomogeneous distribution is measured again, and the processes of the above steps are repeated.
- the measurement control unit 111 calculates the local Bo shimming currents (C 1 and C 2 ) and the adjusted excitation frequencies (f0 1 ′ and f0 2 ′) in the selection region 1 and the selection region 2 set in step 807. To control the excitation of the selection region 1 and the selection region 2.
- the adjusted excitation frequency (f0 1 ′) and the local Bo shimming current (C 1 ) are set according to the imaging parameters of the selection region 1 and the selection region 1 Start excitation.
- the actual excitation frequency of each RF pulse such as slice selection is calculated using the set adjusted excitation frequency (f0 1 ′).
- the set local Bo shimming current (C 1 ) is passed through each shim coil (903), and the calculated RF pulse (901) of each excitation frequency is irradiated.
- the adjusted excitation frequency (f0 2 ′) and the local Bo shimming current (C 2 ) are set according to the imaging parameters of the selection region 2, and the excitation of the selection region 2 is started.
- the actual excitation frequency of each RF pulse such as slice selection is calculated using the set post-adjustment excitation frequency (f0 2 ′).
- the set local Bo shimming current (C 2 ) is passed through each shim coil (904), and the calculated RF pulse (902) of each excitation frequency is irradiated.
- the MRI apparatus and the region imaging method of the third embodiment calculate the local Bo shimming current and the adjusted excitation frequency according to the static magnetic field inhomogeneity for each selected region, and the local Bo for each selected region.
- Each selected region is excited by changing the shimming current and the adjusted excitation frequency.
- the optimal local Bo shimming current and adjusted excitation frequency are used according to the static magnetic field inhomogeneity for each selection area. It is possible to reduce image quality deterioration caused by the static magnetic field non-uniform distribution.
- Example 4 a fourth embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- the response time of a local Bo shimming current of a higher-order component equal to or higher than the secondary component is extended as compared with the local Bo shimming current of a primary component (x, y, z component).
- the fourth embodiment only the primary component of the local Bo shimming current is changed in the selection regions 1 and 2, and the high-order component is shared between the selection region 1 and the selection region 2.
- the outline of the processing flow for setting the local Bo shimming current and the adjusted excitation frequency is the same as that of the third embodiment described above, but the details of the algorithm for calculating the adjusted excitation frequency for the selection region 2 are different.
- the fourth embodiment will be described in detail mainly focusing on the points different from the third embodiment.
- the local Bo shimming current (C 1 ) and the adjusted excitation frequency (f0 1 ′) for the selected region 1 are calculated.
- the calculation method is the same as in the third embodiment.
- the local Bo shimming current (C) 2 [Gx 2 , Gy 2 , Gz 2 , Gxy 2 ] for the selected region 2 is calculated.
- the adjusted excitation frequency (f0 2 ′) for the selection region 2 is calculated.
- the post-shimming static magnetic field inhomogeneity estimation map includes the primary component [Gx 2 , Gy 2 , Gz 2 ] of the local Bo shimming current for the selected region 2 and the local Bo shimming current (C 1 ) for the selected region 1 It is assumed that a higher-order component [Gxy 1 ] of the second or higher order is used. That is, (2) using the formula, when determining the frequency offset?
- the primary component [Gx 2 , Gy 2 , Gz 2 ] of the local Bo shimming current for the selected region 2 and the second or higher order of the local Bo shimming current for the selected region 1 Is estimated using a higher-order component [Gxy 1 ].
- the subsequent processing is the same as in the third embodiment.
- the excitation frequency calculator 304 sets the adjusted excitation frequency for exciting the selected region so that the static magnetic field inhomogeneity in the selected region is reduced, as in the third embodiment.
- C 2 [ ⁇ f0 2 , Gx 2 , Gy 2 , Gz 2 , Gxy 1 ]) in the expression (2).
- the post-shimming static magnetic field inhomogeneity estimation map of the selection region 2 when obtaining the frequency offset ⁇ f0 2 > of the selection region 2 of the expression (3) is expressed by using the primary component [Gx 1 , Gy 1 , Gzv] and a higher-order component [Gxy 1 ] of the second or higher order of the local Bo shimming current for the selected region 1 is estimated.
- step 805 the excitation frequency calculation unit 304 sets the adjusted excitation frequency (f0 2 ′) of the selection region 2 so that the static magnetic field inhomogeneity of the selection region 2 is reduced.
- the second-order or higher-order component [Gxy 2 ] of the local Bo shimming current for the selected region 2 is used instead of the second-order or higher-order component [Gxy 1 ] of the local Bo-shimming current for the selected region 1 is used. Details are as described above.
- step 807 the setting unit 306 uses the local Bo shimming current and the adjusted excitation frequency at the time of exciting the selection region 1 as the primary and secondary or higher of the local Bo shimming current (C 1 ) calculated in step 802.
- the high-order component and the adjusted excitation frequency (f0 1 ′) calculated in step 803 are set, respectively, and control data of an imaging sequence described later is generated and notified to the measurement control unit 111.
- the measurement control unit 111 calculates the local Bo shimming currents (C 1 and C 2 ) and the adjusted excitation frequencies (f0 1 ′ and f0 2 ′) in the selection region 1 and the selection region 2 set in step 807. Used to control the imaging of the selection area 1 and the selection area 2.
- the primary and secondary higher-order components of the adjusted excitation frequency (f0 1 ′) and the local Bo shimming current (C 1 ) of the selected region 1 To start excitation of the selected region 1.
- the RF pulse (901) is irradiated while the set local Bo shimming current (C 1 ) flows through each shim coil (903).
- the second-order and higher components of the local Bo shimming current have the same value between the selected regions.
- Example 5 uses a local Bo shimming current common to each region. Therefore, the fifth embodiment is suitable when the selection area 2 in FIG. 2 (c) is included in the selection area 1. Hereinafter, Example 5 will be described in detail.
- the excitation frequency calculator 304 sets the adjusted excitation frequency for exciting the selected region so that the static magnetic field inhomogeneity in the selected region is reduced, as in the third embodiment.
- used local Bo shimming current for the selected region 1 (C) 1 [Gx 1, Gy 1, Gz 1, Gxy 1].
- C 2 [ ⁇ f0 2 , Gx 1 , Gy 1 , Gz 1 , Gxy 1 ]) in the expression (2).
- the post-shim static magnetic field inhomogeneity estimation map of the selection region 2 when the frequency offset ⁇ f0 2 > of the selection region 2 of the equation (3) is obtained is expressed as a local Bo shimming current [Gx 1 , Gy 1 , Gz 1 , Gxy 1 ] is used for estimation.
- step 805 the excitation frequency calculation unit 304 sets the adjusted excitation frequency (f0 2 ′) of the selection region 2 so that the static magnetic field inhomogeneity of the selection region 2 is reduced.
- the local Bo shimming current (C) 1 [Gx 1 , Gy 1 , Gz 1 , Gxy 1 ] for the selected region 1 is used as the local Bo shimming current for the selected region 2. Details are as described above.
- the excitation of the selected region 1 is started by setting the adjusted excitation frequency (f0 1 ') and the local Bo shimming current (C 1 ) of the selected region 1. To do. At that time, as in FIG. 9, the RF pulse (901) is irradiated while the set local Bo shimming current (C 1 ) flows through each shim coil (903).
- the local Bo shimming current (C 2) as a selection region 1 local Bo shimming current selection region 2 (C 1 ) Is used as it is, and excitation of the selected region 2 is started.
- the RF pulse (1102) is irradiated while the set local Bo shimming current (C 1 ) flows through each shim coil (1104).
- Example 6 a sixth embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- the sixth embodiment is an example in which any one of the above-described third to fifth embodiments is carried out for the purpose of ensuring the minimum required static magnetic field uniformity.
- Example 5 As described above, from the viewpoint of static magnetic field inhomogeneity, the method of Example 5 has the lowest uniformity, then becomes the method of Example 4, and the method of Example 3 has the highest uniformity. Therefore, in the present embodiment 6, the method of each embodiment is determined in this order, and the method of the embodiment in which a predetermined determination criterion is satisfied is adopted, and the local Bo shimming current of each selected region is adjusted. Set the excitation frequency.
- Example 6 will be described in detail.
- step 1200 the selection area 1 and the selection area 2 are set on the positioning image as in the third embodiment.
- step 1201 the static magnetic field inhomogeneous distribution of the area including the selection area 1 and the selection area 2 is measured as in the third embodiment.
- step 1202 the local Bo shimming current (C 1 ) and the adjusted excitation frequency (f0 1 ′) for the selected region 1 are calculated according to the third embodiment described above. Details are as described in Example 3 above.
- step 1203 the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f 0 2 ′) for the selected region 2 are calculated according to Example 3 described above.
- the determination unit 305 determines that the selection field 1 and 2 are static magnetic field inhomogeneity (that is, post-shimming static magnetic field non-uniformity) using the local Bo shimming current and the adjusted excitation frequency set in steps 1202 and 1203. It is determined whether or not the uniform estimation map) satisfies a predetermined determination criterion.
- the static magnetic field inhomogeneous distribution in the post-shimming static magnetic field inhomogeneity estimation map has a variance that is less than a certain threshold (for example, 100 Hz), or the difference between the maximum and minimum is less than the threshold And so on. If the determination condition is satisfied (Yes), the process proceeds to step 1209. If the determination condition is not satisfied (No), the process proceeds to step 1205.
- step 1205 the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f 0 2 ′) of the selected region 2 are calculated again according to the method of Example 4 described above. Details are as described in Example 4 above.
- step 1206 the determination unit 305 determines that the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f0 2 ′) set in step 1205 are used for the selected region 2 and that the static magnetic field is not uniform (that is, after shimming). It is determined whether or not the static magnetic field inhomogeneity estimation map) satisfies a predetermined criterion.
- the determination criteria are the same as in step 1204. When the determination criterion is satisfied (Yes), the process proceeds to Step 1209, and when the determination criterion is not satisfied (No), the process proceeds to Step 1207.
- step 1207 the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f 0 2 ′) of the selected region 2 are calculated again according to the method of Example 3 described above.
- step 1208 the determination unit 305 determines that the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f0 2 ′) set in step 1207 are used for the selected region 2 and that the static magnetic field is not uniform (that is, after shimming). It is determined whether or not the static magnetic field inhomogeneity estimation map) satisfies a predetermined criterion.
- the determination criteria are the same as in step 1204. When the determination criterion is satisfied (Yes), the process proceeds to Step 1209, and when the determination criterion is not satisfied (No), the process proceeds to Step 1210.
- the setting unit 306 sets the value calculated in step 1202 to the imaging parameter for setting the local Bo shimming current (C 1 ) and the adjusted excitation frequency (f0 1 ′) when exciting the selected region 1. To do. Similarly, the setting unit 306 satisfies the determination condition in the imaging parameters for setting the local Bo shimming current (C 2 ) and the adjusted excitation frequency (f0 2 ′) when exciting the selection region 2, according to the third to third embodiments. Set the value calculated by any of 5 methods.
- the measurement control unit 111 calculates the local Bo shimming currents (C 1 and C 2 ) and the adjusted excitation frequencies (f0 1 ′ and f0 2 ′) in the selection region 1 and the selection region 2 set in step 1209. Used to control the imaging of the selection area 1 and the selection area 2. It is assumed that the post-adjustment excitation frequency and local Bo shimming current changing processes are in accordance with Examples 3, 4, and 5.
- step 1211 the determination unit 305 gives a warning display or a warning sound to the operation / display unit 118 that the imaging condition needs to be changed because the determination criterion is not satisfied.
- step 1212 the imaging parameters and the like are reset so that the static magnetic field inhomogeneous distribution in the selection area 1 and the selection area 2 satisfies a predetermined determination criterion. For example, readjustment of selected areas 1 and 2 is performed, TR is extended in the case of SSFP sequence, and excitation position is readjusted in the case of pencil beam. Thereafter, the process proceeds to step 1201, the static magnetic field inhomogeneous distribution is measured again, and the processes in the above steps are repeated.
- a predetermined determination criterion For example, readjustment of selected areas 1 and 2 is performed, TR is extended in the case of SSFP sequence, and excitation position is readjusted in the case of pencil beam.
- the MRI apparatus and the area imaging method according to the sixth embodiment obtain shimming currents for one selected area by a plurality of methods, and satisfy the predetermined determination criteria.
- the shimming current with the least change in shimming current between the two is set for the one selected region.
- the simplest shimming current can be obtained and set according to the static magnetic field inhomogeneous distribution of the selected area to be excited, so that the imaging for each selected area is optimal according to the static magnetic field inhomogeneous distribution of the selected area Can be
- Example 3 to 6 the frequency control for each selected region, that is, an example in which the adjusted excitation frequency (f0 ′) of the RF pulse corresponding to the static magnetic field inhomogeneity for each selected region is obtained and applied is described. did. However, if the main sequence that excites the selected region is the SSFP sequence, the phase rotation amount corresponding to the frequency corresponding to the static magnetic field inhomogeneity ( ⁇ f0 + ⁇ f0>) from the adjusted excitation frequency (f0 ') From this, the increment amount (RF-Phase) may be obtained, and phase control may be applied to the RF pulse in the selected region.
- Example 7 a seventh embodiment of the MRI apparatus and the region imaging method of the present invention will be described.
- this Example 7 when a single region is a selected region, an SSFP sequence is a main sequence, and a prepulse sequence preceding this main sequence is used together, frequency control is applied to the prepulse sequence, and the main (SSFP) Phase control is applied to the sequence. Therefore, this Example 7 is basically the same as Example 1 or 2 described above, but the adjusted excitation frequency (f0 ′) obtained so as to reduce the static magnetic field inhomogeneity in the selected region is prepulsed. The excitation frequency of the sequence is different. In the following, portions different from the first or second embodiment will be mainly described, and description of the same contents will be omitted.
- Each function of the arithmetic processing unit 114 for realizing the region imaging method of the seventh embodiment is the same as the functional block of the first embodiment or the second embodiment, and the description thereof is omitted.
- processing flow of Example 7 is basically the same as the processing flow shown in FIG. 4 described in the first embodiment or the second embodiment, but the processing contents of steps 405 and 406 are different.
- steps 405 and 406 are different.
- Only the parts different from the above-described Example 1 or Example 2 will be described, and description of the same parts will be omitted.
- step 405 the setting unit 306 sets the adjusted excitation frequency (f0 ′) calculated in step 403 as the excitation frequency of the pre-pulse sequence in addition to the processing content of step 404 of the first embodiment. Then, control data of the prepulse sequence and the main sequence is generated and notified to the measurement control unit 111.
- step 406 the measurement control unit 111 controls the execution of the pre-pulse sequence using the local Bo shimming current (C) set in step 405 and the adjusted excitation frequency (f0 ′), and the same local Bo shimming current ( The execution of the main sequence is controlled using C), RF-Phase, and excitation frequency (f0).
- C local Bo shimming current
- f0 ′ adjusted excitation frequency
- f0 excitation frequency
- the MRI apparatus and the region imaging method of the seventh embodiment change the static magnetic field inhomogeneous distribution of the selected region when the selected region is excited using the imaging sequence composed of the pre-pulse sequence and the SSFP sequence. Accordingly, the shimming current and the adjusted excitation frequency (f0 ′) are calculated, and the RF-Phase corresponding to the adjusted excitation frequency (f0 ′) is calculated. Then, the shimming current and the adjusted excitation frequency (f0 ′) are applied to the prepulse sequence, and the same shimming current and the increment amount (RF-Phase) are applied to the SSFP sequence.
- the position of the band artifact is controlled by RF-Phase, so even if the operator changes the excitation frequency (f0), the position of the band artifact is not changed. It can be fixed.
- the SSFP sequence is a main sequence, and when a pre-pulse sequence preceding this main sequence is used in combination, the pre-pulse sequence is selected using frequency control.
- the selected region 2 is excited using phase control for the SSFP sequence. Therefore, in the eighth embodiment, the pre-pulse sequence is basically the same as those in the third to sixth embodiments, and the SSFP sequence is the same as in the first or second embodiment. However, the control of shimming currents in both selected regions is the same as in the third to sixth embodiments. In the following, portions different from the above-described embodiments will be described, and description of the same contents will be omitted.
- Each function of the arithmetic processing unit 114 for realizing the area imaging method of the eighth embodiment will be described with reference to FIGS. 13 (a) and 13 (b).
- Each function of the present embodiment is basically the same as the functional block described with reference to FIGS. 3 (a) and 3 (b) in the first embodiment or the second embodiment, but is selected by a pre-pulse sequence for frequency control.
- the frequency calculation unit 304 In order to excite the region 1, the frequency calculation unit 304 notifies the setting unit 306 of the adjusted excitation frequency (f0 ′) calculated.
- the RF-Phase calculation unit 305 calculates an RF-Phase corresponding to the adjusted excitation frequency (f0 2 ′) calculated in Step 805 as the RF-Phase when exciting the selected region 2.
- the setting unit 306 uses the local Bo shimming current and the excitation frequency at the time of exciting the selected region 1 as the local Bo shimming current (C 1 ) calculated in Step 802 and the adjusted excitation frequency (Step 803).
- f0 1 ′) are respectively set, control data of a pre-pulse sequence to be described later is generated and notified to the measurement control unit 111.
- Each RF-Phase is set to generate control data for a main (SSFP) sequence, which will be described later, and notify the measurement control unit 111 of the control data.
- step 810 the measurement control unit 111 controls excitation of the selected region 1 by the prepulse sequence using the local Bo shimming current (C 1 ) and the adjusted excitation frequency (f0 1 ′) set in step 807. Thereafter, excitation of the selected region 2 by the main (SSFP) sequence is controlled using the local Bo shimming current (C 1 ), RF-Phase, and excitation frequency (f 0) set in step 807.
- C 1 the local Bo shimming current
- RF-Phase RF-Phase
- the MRI apparatus and the region imaging method of the eighth embodiment calculate the local Bo shimming current and the adjusted excitation frequency or increment amount (RF-Phase) according to the static magnetic field inhomogeneous distribution for each selected region.
- the local Bo shimming current and the adjusted excitation frequency or RF-Phase are changed for each selected region to excite each selected region.
- the optimum local Bo shimming current and adjusted excitation frequency or RF-Phase are used according to the static magnetic field inhomogeneous distribution for each selected region. It is possible to reduce image quality degradation caused by the static magnetic field non-uniform distribution for each.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
本発明に係る領域撮像法は、以下のステップを有する。即ち、
ROI設定を介して選択された領域(以下、選択領域と略記する)を含む領域についての静磁場不均一の空間的分布を表す静磁場不均一分布を測定する。
そして、求めた局所Boシミング電流(C)を用いて選択領域をシミングして、該選択領域の静磁場不均一を低減した状態での、該選択領域を励起する調整後励起周波数(f0')を求める。
静磁場不均一を考慮した調整後励起周波数の算出法の概要について説明する。
次に、(1)式で求めた静磁場不均一分布に対して、励起する選択領域のみを対象とするようにマスク処理を施し、選択領域の静磁場不均一が小さくなるような局所Boシミング電流(C)を算出する。具体的には、選択領域に関して、MRI装置据付時に測定したシムコイルの特性が反映された特性行列をAと、シムコイル電流値をC(=[Δf0,Gx,Gy,Gz,Gxy])としたときの(2)式の行列をCに関して解く。
ここで、Δf0は周波数オフセットである。Gx、Gy、Gzは、それぞれ1次のシムコイルx、y、zが発生する1次の補償磁場成分である。なお、1次のシムコイルは傾斜磁場コイルと兼用されるのが通常なので、この1次の補償磁場を発生するための電流は傾斜磁場コイルに供給される。一方、Gxyは2次のシムコイルxyが発生する2次の補償磁場成分である。2次以上の項としてxyのみを記載したが、一般にはx^2、y^2、yz、zx、(x^2-y^2)などの高次項を含んでも良い。
ここで、f0は静磁場不均一分布を測定する際に用いたdual echo Gradient Echoの励起周波数であり、Δf0は(2)式により求めた周波数オフセットであり、<f0>はシミング後静磁場不均一推定マップを基に求めた周波数オフセットである。つまり、選択領域を励起するRFパルスの調整後励起周波数(f0')は、静磁場不均一が低減された状態に基づいて算出されるものであり、静磁場不均一分布を測定する際に用いたRFパルスの励起周波数に、静磁場不均一に対応する周波数分(Δf0+<f0>)を加算したものとなる。シミング後静磁場不均一推定マップは、シムコイル電流値[Gx,Gy,Gz,Gxy]を流したときに発生する磁場をΔB0に加算して求める。<f0>の算出には、たとえば、選択領域内のシミング後静磁場不均一推定マップの最大値と最小値の平均や、選択領域内の静磁場不均一推定マップの相加平均や、選択領域内の静磁場不均一推定マップのヒストグラムの平均などにより算出する。
次に、(3)式の調整後励起周波数(f0')に相当する、RFパルス照射位相の照射毎のインクリメント量である、RF-Phaseの算出方法について説明する。
RF-Phase = [ 180 + 360*TR*(f0'-f0) ] % 360
= [180 + 360*TR*(Δf0+<f0>)] % 360 (4)
ただし、%は剰余演算子を表し、A%Bの演算はAをBで割った余りを意味する。つまり、RF-Phaseは、静磁場不均一に対応する周波数分(Δf0+<f0>)に相当する位相回転をRFパルスの照射位相に与える位相量である。(4)式から、RF-Phaseは調整後励起周波数(f0')の算出を介して求めることも可能あり、静磁場不均一に対応する周波数分(Δf0+<f0>)から直接求めることも可能であることが理解される。ここで、SSFPシーケンスの実際の励起周波数であるf0は、前計測(プリスキャン)として撮像スライス断面のスペクトルを計測し、水のスペクトルのピーク周波数に自動設定してもよいし、操作者が上記スペクトルの形状からでマニュアル調整して設定してもよいとする。f0'は(3)式により算出した励起周波数である。
本発明のMRI装置及び領域撮像方法は、撮像条件としての撮像パラメータに応じて、上記位相制御又は周波数制御のいずれかを行う。制御方法の選択に係る撮像パラメータとして、選択領域の数と撮像シーケンスのシーケンス種とがある。ここで撮像シーケンスが画像用のエコー信号の計測を行うメインシーケンスとそのメインシーケンスに先行するプリパルスシーケンスとからなる場合には、シーケンス種とはメインシーケンスのシーケンス種を意味する。それらの撮像パラメータに応じた制御方法の選択を行う処理フローを、図2のフローチャートに基づいて説明する。
本発明のMRI装置及び領域撮像方法の実施例1について説明する。本実施例1は、静磁場不均一に応じて、シムコイルに流す局所Boシミング電流とRFパルスのインプリメント量(RF-Phase)とを算出して、算出した局所Boシミング電流とRF-Phaseとを用いて撮像する。以下、図3(a),4に基づいて本実施例1を詳細に説明する。
最初に、本実施例1の領域撮像方法を実現するための演算処理部114の各機能を、図3(a)に示す機能ブロック図に基づいて説明する。本実施例1に係る演算処理部114は、領域選択部301と、静磁場不均一分布測定部302と、局所Boシミング電流算出部303と、励起周波数算出部304と、RF-Phase算出部305と、設定部306を有して成る。
次に、上記各機能部が連携して行なう本実施例1の処理フローを図4に示すフローチャートに基づいて説明する。本処理フローは、予めプログラムとして内部記憶装置115に記憶されており、演算処理部114が内部記憶装置115からそのプログラムを読み込んで実行することにより実施される。以下、各処理ステップの処理内容を詳細に説明する。
以上までが本実施例1の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例2を説明する。前述の実施例1では、自動的に適切なRF-Phaseを算出したが、ROI内の静磁場不均一が1/TR [Hz]以上の場合にはROI内にバンドアーチファクトが混入してしまう。実施例2では、このような場合に、操作者がバンドアーチファクトの位置を容易に調整できるようにする。具体的には、局所Boシミング電流および調整後励起周波数(f0')を設定する処理フローの概要は前述の実施例1と同じであるが、操作者が設定入力する位相を追加してインクリメント量(RF-Phase)を算出することで、バンドアーチファクトの位置を調整できるようにする。
初めに、本実施例2の、操作者がバンドアーチファクトの位置を調整する場合のインクリメント量(RF-Phase)算出方法の概要を説明する。
= [180 +Φ+ 360*TR*(Δf0+<f0>)] % 360 (5)
ここで、各変数は、(4)式と同様である。この場合にも、(5)式に示すように、前述の実施例1と同様に、調整後励起周波数(f0')を用いることなく、静磁場不均一に対応する周波数分(Δf0+<f0>)に相当する位相回転量に操作者により設定された調整値Φを加算してRF-Phaseを求めてもよい。
以降の処理は実施例1と同様である。
次に、本実施例2の領域撮像方法を実現するための演算処理部114の各機能を図3(b)に基づいて説明する。本実施例2に係る演算処理部114の機能ブロックは、前述の実施例1と基本的には同様であるが、RF-Phase算出部305の機能が異なる。以下、前述の実施例1と異なる箇所のみを説明し、同一箇所の説明は省略する。
次に、上記各機能部が連携して行なう本実施例2の処理フローを図4に示すフローチャートに基づいて説明する。本実施例2の処理フローは、前述の実施例1で説明した図4の処理フローと基本的には同様であるが、ステップ404の処理内容が異なる。以下、前述の実施例1と異なる箇所のみを説明し、同一箇所の説明は省略する。
以上までが本実施例2の処理フローの説明である。
本発明のMRI装置及び領域撮像方法の実施例3について説明する。本実施例3は、領域毎に、シムコイルに流す局所Boシミング電流と調整後励起周波数とを算出して、領域毎に順次撮像する。したがって、本実施例3は、図6(a)の選択領域1と選択領域2が共通領域を持たない場合、又は、図6(b)の選択領域1と選択領域2が共通領域を持つ場合の撮像に好適である。換言すれば、2つの領域が互いに少なくとも一部が重ならない場合に本実施例3は好適となる。なお、本実施例3においては、2つの領域が重なる共通領域は2回撮像されることになる。以下、図7,8に基づいて本実施例3を詳細に説明する。
最初に、本実施例3の領域撮像方法を実現するための演算処理部114の各機能を、図7に示す機能ブロック図に基づいて説明する。本実施例3に係る演算処理部114の機能ブロックは、前述の実施例1で説明した図3(a)に示すブロック図と比較して、RF-Phase算出部305がなく、判定部705を有して成る。以下、前述の実施例1の機能ブロック図と異なる部分のみ説明し、同じ機能内容については説明を省略する。
静磁場不均一分布測定部302は、実施例1と同様である。
局所Boシミング電流算出部303は、選択領域の静磁場不均一が小さくなるような局所Boシミング電流(C)を算出する。局所Boシミング電流の算出の詳細は、実施例1と同様である。
次に、上記各機能部が連携して行なう本実施例3の処理フローを図8に示すフローチャートに基づいて説明する。本処理フローは、予めプログラムとして内部記憶装置115に記憶されており、演算処理部114が内部記憶装置115からそのプログラムを読み込んで実行することにより実施される。以下、各処理ステップの処理内容を詳細に説明する。
以上までが本実施例3の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例4を説明する。一般的に、2次成分以上の高次成分の局所Boシミング電流の応答時間は、1次成分(x、y、z成分)の局所Boシミング電流と比較して延長する。そこで、実施例4は、局所Boシミング電流の1次成分のみ選択領域1と2で変更し、高次成分を選択領域1と選択領域2とで共有する。その結果、局所Boシミング電流および調整後励起周波数を設定する処理フローの概要は前述の実施例3と同じであるが、選択領域2に対する調整後励起周波数の算出アルゴリズムの詳細が異なる。以下、主に前述の実施例3と異なる箇所を中心にして、本実施例4を詳細に説明する。
最初に、本実施例4の各領域の静磁場不均一分布に応じた調整後励起周波数の算出法の概要を説明する。
以降の処理は実施例3と同様である。
次に、本実施例4の領域撮像方法を実現するための演算処理部114の各機能を説明する。本実施例4に係る演算処理部114の機能ブロックは、前述の実施例3と同様であるが、励起周波数算出部304の機能が異なる。以下、前述の実施例3と異なる箇所のみを説明し、同一箇所の説明は省略する。
次に、上記各機能部が連携して行なう本実施例4の処理フローを説明する。本実施例4の処理フローは、前述の実施例3で説明した図8の処理フローと同様であるが、ステップ805、807、810の処理内容が異なる。以下、前述の実施例3と異なるステップのみ説明し、同じ処理内容のステップの説明を省略する。
以上までが本実施例4の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例5について説明する。実施例5は、各領域で共通の局所Boシミング電流を用いる。したがって、本実施例5は、図2(c)の選択領域2が選択領域1に含まれている場合に好適である。以下、本実施例5を詳細に説明する。
最初に、本実施例5の領域撮像方法を実現するための演算処理部114の各機能を説明する。本実施例5に係る演算処理部114の機能ブロックは、前述の実施例3と同様であるが、励起周波数算出部304の機能が異なる。以下、前述の実施例3と異なる箇所のみを説明し、同一箇所の説明は省略する。
次に、上記各機能部が連携して行なう本実施例5の処理フローを説明する。本実施例5の処理フローは、前述の実施例3と比較し、ステップ804の処理が無く、ステップ805、810の処理内容が異なる。以下前述の実施例3と異なるステップのみ説明し、同じ処理内容のステップの説明を省略する。
以上までが本実施例5の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例6について説明する。本実施例6は、必要最低限の静磁場均一度を確保する目的で、前述の各実施例3~5のいずれかを実施する例である。
本実施例6の処理フローを図12に示すフローチャートに基づいて説明する。本処理フローは、予めプログラムとして内部記憶装置115に記憶されており、演算処理部114が内部記憶装置115からそのプログラムを読み込んで実行することにより実施される。以下、各処理ステップの処理内容を詳細に説明する。
ステップ1203で、前述の実施例3に従い、選択領域2に対する局所Boシミング電流(C2)と調整後励起周波数(f02')が算出される。
以上までが本実施例6の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例7について説明する。本実施例7は、単一領域を選択領域とし、SSFPシーケンスをメインシーケンスとして、このメインシーケンスに先行するプリパルスシーケンスを併用する場合に、プリパルスシーケンスには、周波数制御を適用し、メイン(SSFP)シーケンスには位相制御を適用する。そのため、本実施例7は、基本的には前述の実施例1又は2と同様であるが、選択領域の静磁場不均一が小さくなるようなに求めた調整後励起周波数(f0')をプリパルスシーケンスの励起周波数とすることが異なる。以下、主に前述の実施例1又は2と異なる箇所を説明し、同じ内容については説明を省略する。
本実施例7の領域撮像方法を実現するための演算処理部114の各機能は、前述の実施例1又は実施例2の機能ブロックと同様であるので説明を省略する。
本実施例7の処理フローは、基本的には前述の実施例1又は実施例2で説明した図4に示した処理フローと同様であるが、ステップ405、406の処理内容が異なる。以下、前述の実施例1又は実施例2と異なる箇所のみを説明し、同一箇所の説明は省略する。
以上までが本実施例7の処理フローの説明である。
次に、本発明のMRI装置及び領域撮像方法の実施例8について説明する。本実施例8は、複数領域を選択領域とする場合において、SSFPシーケンスをメインシーケンスとし、このメインシーケンスに先行するプリパルスシーケンスを併用する場合に、プリパルスシーケンスには周波数制御を用いて選択領域1を励起し、SSFPシーケンスには位相制御を用いて選択領域2を励起する。そのため、本実施例8は、基本的には、プリパルスシーケンスについては前述の実施例3~6と同様であり、SSFPシーケンスについては、前述の実施例1又は2と同様である。ただし、両選択領域のシミング電流の制御については、前述の実施例3~6と同様である。以下、前述の各実施例と異なる箇所を説明し、同じ内容については説明を省略する。
本実施例8の領域撮像方法を実現するための演算処理部114の各機能を図13(a)(b)に基づいて説明する。本実施例の各機能は、前述の実施例1又は実施例2で図3(a)(b)を用いて説明した機能ブロックと基本的には同じであるが、周波数制御のプリパルスシーケンスで選択領域1を励起するために、周波数算出部304から設定部306に算出した調整後励起周波数(f0')を通知する。
上記各機能部が連携して行なう本実施例8の処理フローを説明する。本実施例8の処理フローは、前述の実施例3で説明した図8の処理フローと基本的には同様であるが、ステップ807,810の処理内容が異なる。以下、前述の実施例3と異なる箇所のみを説明し、同一箇所の説明は省略する。
以上までが本実施例8の処理フローの説明である。
Claims (15)
- 被検体が配置される撮像領域に静磁場を発生する静磁場発生部と、
シミング電流が供給されて、前記静磁場の空間的な不均一である静磁場不均一を低減する補償磁場を発生するシミング部と、
所定のパルスシーケンスに基づいて前記被検体からのエコー信号の計測を制御する計測制御部と、
前記エコー信号を用いて前記静磁場不均一を求める静磁場不均一測定部と、
前記静磁場不均一に基づいて、前記シミング電流を求めるシミング電流算出部と、
前記静磁場不均一に基づいて、パルスシーケンスに用いるRFパルスの照射位相のインクリメント量(RF-Phase)を求めるRF-Phase算出部と、
前記シミング部に前記シミング電流が供給されて前記静磁場不均一が低減された状態で、前記RFパルスの照射位相を前記インクリメント量ずつ変えて、前記被検体からのエコー信号の計測を制御する計測制御部と、
を備えたことを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置において、
前記静磁場不均一を低減する選択領域の設定を受け付ける領域選択部を備え、
前記シミング電流算出部は、前記選択領域の静磁場不均一に基づいて前記シミング電流を求め、
前記RF-Phase算出部は、前記選択領域の静磁場不均一に基づいて前記インクリメント量(RF-Phase)を求めることを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置において、
前記静磁場不均一測定部が前記静磁場不均一を求める際に用いたRFパルスの励起周波数に、該静磁場不均一に対応する周波数分を加算して、前記静磁場不均一が低減された状態での前記RFパルスの励起周波数である調整後励起周波数を算出する励起周波数算出部を備え、
前記RF-Phase算出部は、前記調整後励起周波数に基づいて、前記インクリメント量(RF-Phase)を算出することを特徴とする磁気共鳴イメージング装置。 - 請求項3記載の磁気共鳴イメージング装置において、
前記励起周波数算出部は、前記静磁場不均一が低減された状態に基づいて、前記調整後励起周波数を算出することを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置において、
前記RF-Phase算出部は、前記静磁場不均一に対応する周波数分に相当する位相回転量として前記インクリメント量(RF-Phase)を算出することを特徴とする磁気共鳴イメージング装置。 - 請求項5記載の磁気共鳴イメージング装置において、
前記インクリメント量(RF-Phase)の調整値(Φ)が入力される調整値入力部を備え、
前記RF-Phase算出部は、前記調整値(Φ)を用いて前記インクリメント量(RF-Phase)を算出することを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至6のいずれか一項に記載の磁気共鳴イメージング装置において、
前記パルスシーケンスはSSFPシーケンスであり、
前記計測制御部は、前記SSFPシーケンスのRFパルスの照射位相に前記インクリメントずつ変える位相制御を適用することを特徴とする磁気共鳴イメージング装置。 - 請求項3又は4に記載の磁気共鳴イメージング装置において、
前記パルスシーケンスは、プリパルスシーケンスとSSFPシーケンスとを有して成り、
前記計測制御部は、前記プリパルスシーケンスのRFパルスには、その励起周波数を前記調整後励起周波数とする周波数制御を適用し、前記SSFPシーケンスのRFパルスには、その照射位相に前記位相制御を適用することを特徴とする磁気共鳴イメージング装置。 - 請求項8記載の磁気共鳴イメージング装置において、
前記プリパルスシーケンスと前記SSFPシーケンスとが異なる選択領域を励起する場合に、
前記シミング電流算出部は、前記選択領域毎の静磁場不均一に基づいて、前記選択領域毎のシミング電流を求め、
前記励起周波数算出部は、前記選択領域毎に、前記静磁場不均一が低減された状態での前記RFパルスの調整後励起周波数を算出し、
前記RF-Phase算出部は、前記SSFPシーケンスにより励起される選択領域での前記インクリメント量(RF-Phase)を、該選択領域での調整後励起周波数に基づいて算出することを特徴とする磁気共鳴イメージング装置。 - 請求項2記載の磁気共鳴イメージング装置において、
前記領域選択部は複数の選択領域の設定を受け付け、
前記シミング電流算出部は、前記選択領域毎の静磁場不均一に基づいて、前記選択領域毎のシミング電流を求め、
前記励起周波数算出部は、前記選択領域毎に、前記静磁場不均一が低減された状態での前記RFパルスの調整後励起周波数を算出し、
前記RF-Phase算出部は、前記選択領域毎の調整後励起周波数に基づいて、前記選択領域毎の前記インクリメント量(RF-Phase)を算出することを特徴とする磁気共鳴イメージング装置。 - 被検体が配置される撮像領域に静磁場を発生する静磁場発生部と、
シミング電流が供給されて、前記静磁場の空間的な不均一である静磁場不均一を低減する補償磁場を発生するシミング部と、
所定のパルスシーケンスに基づいて前記被検体からのエコー信号の計測を制御する計測制御部と、
前記エコー信号を用いて前記静磁場不均一を求める静磁場不均一測定部と、
前記静磁場不均一に基づいて前記シミング電流を算出し、前記シミング部に前記シミング電流が供給されて前記静磁場不均一が低減された状態での、前記パルスシーケンスのRFパルスの調整後励起周波数を算出する演算処理部と、
複数の選択領域の設定を受け付ける領域選択部と、
を備え、
前記演算処理部は、前記複数の選択領域の重なり方に応じて、前記選択領域毎に、前記シミング電流と前記調整後励起周波数を算出し、
前記計測制御部は、前記シミング部に前記シミング電流が供給されて前記静磁場不均一が低減された状態で、前記パルスシーケンスのRFパルスに前記調整後励起周波数を用いて前記エコー信号の計測を制御する
ことを特徴とする磁気共鳴イメージング装置。 - 請求項11記載の磁気共鳴イメージング装置において、2つの領域が互いに少なくとも一部が重ならない場合に、
前記演算処理部は、前記選択領域毎にその前記静磁場の不均一分布に応じて、前記シミング電流と前記励起周波数を算出し、
前記計測制御部は、前記選択領域毎に前記シミング電流と前記励起周波数を変えて、それぞれの選択領域からエコー信号の計測を制御することを特徴とする磁気共鳴イメージング装置。 - 請求項11記載の磁気共鳴イメージング装置において、
前記演算処理部は、前記選択領域毎の前記シミング電流の2次以上の成分を同じ値とすることを特徴とする磁気共鳴イメージング装置。 - 請求項11記載の磁気共鳴イメージング装置において、1つの領域が他の一つの領域を含む場合に、
前記演算処理部は、前記選択領域毎の前記シミング電流の全成分をそれぞれ同じ値とすることを特徴とする磁気共鳴イメージング装置。 - シミング電流が供給されて静磁場不均一を低減するシミング部を備えた磁気共鳴イメージング装置における領域撮像方法であって、
前記静磁場不均一を計測する静磁場不均一測定ステップと、
前記静磁場不均一に基づいて、前記シミング電流を求めるシミング電流算出ステップと、
前記静磁場不均一に基づいて、前記パルスシーケンスに用いるRFパルスの照射位相のインクリメント量(RF-Phase)を求めるRF-Phase算出ステップと、
前記シミング部に前記シミング電流が供給されて前記静磁場不均一が低減された状態で、前記RFパルスの照射位相を前記インクリメント量ずつ変えて、被検体からエコー信号を計測する計測ステップと、
を備えたことを特徴とする領域撮像方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013553217A JP6071905B2 (ja) | 2012-01-11 | 2012-12-12 | 磁気共鳴イメージング装置及び領域撮像方法 |
US14/365,278 US9678184B2 (en) | 2012-01-11 | 2012-12-12 | Method for increment of RF-phase based on static magnetic field inhomogeneity in MRI apparatus |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-002644 | 2012-01-11 | ||
JP2012002644 | 2012-01-11 | ||
JP2012119153 | 2012-05-25 | ||
JP2012-119153 | 2012-05-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013105384A1 true WO2013105384A1 (ja) | 2013-07-18 |
Family
ID=48781344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2012/082139 WO2013105384A1 (ja) | 2012-01-11 | 2012-12-12 | 磁気共鳴イメージング装置及び領域撮像方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9678184B2 (ja) |
JP (1) | JP6071905B2 (ja) |
WO (1) | WO2013105384A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019195590A (ja) * | 2018-05-11 | 2019-11-14 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置およびマルチスライス撮像方法 |
US11237235B2 (en) | 2015-04-10 | 2022-02-01 | Synaptive Medical Inc. | Shimming coils for magnetic resonance imaging |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013216529B4 (de) * | 2013-08-21 | 2019-05-23 | Siemens Healthcare Gmbh | Verfahren insbesondere zur patientenadaptiven B0-Homogenisierung von MR-Systemen unter Verwendung unterschiedlicher Typen von Shim-Spulen |
US10335049B2 (en) * | 2014-04-22 | 2019-07-02 | Case Western Reserve University | Magnetic resonance fingerprinting with steady state precession (MRF-FISP) |
DE102014216612A1 (de) * | 2014-08-21 | 2016-02-25 | Siemens Aktiengesellschaft | Verfahren zur Durchführung mindestens einer Justierungsmessung in einem Magnetresonanzgerät, sowie entsprechend ausgestaltetes Magnetresonanzgerät, Computerprogramm und computerlesbares Speichermedium |
CN105676155B (zh) * | 2014-11-19 | 2019-02-22 | 西门子(深圳)磁共振有限公司 | 磁共振成像***的磁场不均匀性值获取方法和装置和失真校正方法和装置 |
JP5943159B2 (ja) * | 2016-01-04 | 2016-06-29 | 株式会社日立製作所 | 磁気共鳴イメージング装置 |
CN109073721B (zh) * | 2016-04-28 | 2022-02-08 | 皇家飞利浦有限公司 | 处置计划评价工具 |
CA3030629A1 (en) * | 2016-07-11 | 2018-01-18 | Synaptive Medical (Barbados) Inc. | Adaptive shim coils for mr imaging |
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 |
EP3591418A1 (en) * | 2018-07-03 | 2020-01-08 | Koninklijke Philips N.V. | Mri method for b0-mapping |
US11243287B2 (en) * | 2019-10-02 | 2022-02-08 | Synaptive Medical Inc. | Real-time compensation of high-order concomitant magnetic fields |
CN114325520B (zh) * | 2020-09-30 | 2024-07-19 | 西门子(深圳)磁共振有限公司 | 磁体升场方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003210431A (ja) * | 2002-01-23 | 2003-07-29 | Ge Medical Systems Global Technology Co Llc | Fid信号の位相ずれ測定方法、se・ste信号の位相ずれ測定方法、mrイメージング方法およびmri装置 |
JP2005152175A (ja) * | 2003-11-25 | 2005-06-16 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴撮像装置および方法 |
JP2009050605A (ja) * | 2007-08-29 | 2009-03-12 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置,プログラム |
JP2011229632A (ja) * | 2010-04-27 | 2011-11-17 | Hitachi Medical Corp | 磁気共鳴イメージング装置及びシミング方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5670159B2 (ja) * | 2009-11-26 | 2015-02-18 | 株式会社東芝 | 磁気共鳴イメージング装置 |
EP2506026A1 (en) * | 2011-03-29 | 2012-10-03 | Universitätsklinikum Freiburg | Method of dynamically compensating for magnetic field heterogeneity in magnetic resonance imaging |
-
2012
- 2012-12-12 JP JP2013553217A patent/JP6071905B2/ja active Active
- 2012-12-12 WO PCT/JP2012/082139 patent/WO2013105384A1/ja active Application Filing
- 2012-12-12 US US14/365,278 patent/US9678184B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003210431A (ja) * | 2002-01-23 | 2003-07-29 | Ge Medical Systems Global Technology Co Llc | Fid信号の位相ずれ測定方法、se・ste信号の位相ずれ測定方法、mrイメージング方法およびmri装置 |
JP2005152175A (ja) * | 2003-11-25 | 2005-06-16 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴撮像装置および方法 |
JP2009050605A (ja) * | 2007-08-29 | 2009-03-12 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置,プログラム |
JP2011229632A (ja) * | 2010-04-27 | 2011-11-17 | Hitachi Medical Corp | 磁気共鳴イメージング装置及びシミング方法 |
Non-Patent Citations (1)
Title |
---|
N. K. BANGERTER: "SNR Analysis of Multiple Acquisition SSFP", PROCEEDINGS OF INTERNATIONAL SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE, 18 May 2002 (2002-05-18) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11237235B2 (en) | 2015-04-10 | 2022-02-01 | Synaptive Medical Inc. | Shimming coils for magnetic resonance imaging |
JP2019195590A (ja) * | 2018-05-11 | 2019-11-14 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置およびマルチスライス撮像方法 |
JP7139143B2 (ja) | 2018-05-11 | 2022-09-20 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置およびマルチスライス撮像方法 |
Also Published As
Publication number | Publication date |
---|---|
US9678184B2 (en) | 2017-06-13 |
JP6071905B2 (ja) | 2017-02-01 |
US20140327440A1 (en) | 2014-11-06 |
JPWO2013105384A1 (ja) | 2015-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6071905B2 (ja) | 磁気共鳴イメージング装置及び領域撮像方法 | |
JP5523564B2 (ja) | 磁気共鳴イメージング装置および送信感度分布算出方法 | |
US9606205B1 (en) | Magnetic resonance imaging apparatus, RF shimming method, and magnetic resonance imaging method | |
JP6371554B2 (ja) | 磁気共鳴イメージング装置 | |
US9594140B2 (en) | Magnetic resonance imaging apparatus and method for calculating correction value as application amount of refocusing pulse for UTE sequence | |
US9846215B2 (en) | MRI embodiments for controlling an arrangement order of multiple echoes in a k-space | |
WO2016021603A1 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
WO2012138902A1 (en) | B0-based modulation of b1 excitation in mri | |
US20140077805A1 (en) | Magnetic resonance imaging apparatus and flip angle determination method | |
JP5317571B2 (ja) | 磁気共鳴撮影装置 | |
JP6017443B2 (ja) | 磁気共鳴イメージング装置および照射磁場分布計測方法 | |
JP5808659B2 (ja) | 磁気共鳴イメージング装置及びT1ρイメージング法 | |
JP2008307303A (ja) | 磁気共鳴イメージング装置 | |
JP6084430B2 (ja) | 再収束rfパルスのフリップ角制御法及び磁気共鳴イメージング装置 | |
JP5595759B2 (ja) | 磁気共鳴イメージング装置 | |
JP6487554B2 (ja) | 磁気共鳴イメージング装置 | |
US10156616B2 (en) | Nuclear magnetic resonance imaging apparatus and RF shimming method | |
JP4685496B2 (ja) | 磁気共鳴イメージング装置 | |
JP2014087442A5 (ja) | ||
JP2016131847A (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
US9594136B2 (en) | Magnetic resonance imaging apparatus and RF pulse control method | |
JP6103965B2 (ja) | 磁気共鳴イメージング装置及びrfシミング方法 | |
JP2012239561A (ja) | 磁気共鳴イメージング装置及びDualSlice計測方法 | |
JP5283213B2 (ja) | 磁気共鳴イメージング装置 | |
JP2014200571A (ja) | 磁気共鳴イメージング装置及び2項パルス制御法 |
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: 12865401 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013553217 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14365278 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: 12865401 Country of ref document: EP Kind code of ref document: A1 |