LOCAL SHIMMING SYSTEM FOR MAGNETIC RESONANCE IMAGING AND METHOD THEREOF

Abstract

A local shimming system for magnetic resonance imaging and the method thereof, wherein the shimming method comprises the following steps: collecting B0 field map information using two-dimensional gradient echo (301); calculating and evaluating the homogeneity of B0 (302); optimizing the current of each channel shim coil (303); determining whether the minimum standard deviation value of f is obtained (304); outputting an optimal current combination values and setting an optimum current value corresponding to each channel of the shim coil on the current control software (305); and testing and evaluating the homogeneity of B0 to achieve the shimming goal (306).

Claims

1. A shimming system for magnetic resonance imaging, which comprises: a multi-channel local shim coil unit, which is configured to be installed on an inspection table of the magnetic resonance imaging system, wherein the local shim coil unit comprises a local multi-channel shim coil and a radio frequency receiving coil for receiving magnetic resonance signals, and the local radio frequency coil is placed inside the local shim coil and separated by a distance; a computer control system, which is configured to install and set the software controlled by the DC power and calculate field maps and calculate optimization processes; and a DC power system, which is communicatively connected to the computer control system to control the value of current of each channel.

2. The shimming system of claim 1, wherein the multi-channel local shim coil unit further comprises a housing, and the housing is provided with a plurality of holes for tuning and matching the radio frequency receiving coil.

3. The shimming system of claim 1, wherein the channels of the local shim coil unit have a number of 5, and the local shim coils are asymmetrically distributed.

4. The shimming system of claim 1, wherein the local shim coil is a circular loop having a diameter of 4 cm.

5. The shimming system of claim 1, wherein the housing has a bottom provided with a hole as a reserved channel for the ultrasound ablation treatment.

6. The shimming system of claim 5, wherein the hole has a diameter of 6 cm.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The drawings described in the present application are intended to provide a further understanding of the present application and are not intended to limit the present application. In the drawing:

[0023] FIG. 1 is a schematic structural diagram of the local shimming system for magnetic resonance imaging according to an embodiment of the present application.

[0024] FIG. 2 is a schematic perspective structural diagram of the local shim coil unit according to an embodiment of the present application.

[0025] FIG. 3 is a flow chart of the shimming method for magnetic resonance imaging according to an embodiment of the present application.

[0026] FIGS. 4A and 4B show the comparison of the obtained main magnetic field offset distributions after shimming and before shimming, respectively.

[0027] FIGS. 5A and 5B show statistics of the magnetic field frequency distributions within the region of interest ROI after shimming and before shimming, respectively.

[0028] FIGS. 6A and 6B show the temperature change graph of a phantom at normal temperature conditions after shimming and before shimming, respectively.

[0029] In the figures: 10 represents a local shim coil unit, 11 represents a magnetic resonance imaging system, 101 represents an operation room, 102 represents an equipment room, 103 represents a magnetic resonance laboratory, 104 represents a computer control system, 105 represents a network cable, 106 represents a DC power system, 107 represents a waveguide plate and a waveguide tube, 108 represents a DC power cord, 109 represents a local multi-channel shim coil, 110 represents an RF receiving coil.

DESCRIPTION OF THE EMBODIMENTS

[0030] Specific embodiments of the present application will be described below with reference to the accompanying drawings. In the specific embodiments of the present application described hereinafter, some very specific technical features are described for a better understanding of the present application, but it will be apparent to those skilled in the art that not all of the technical features are essential technical features for implementing the present application. Some specific embodiments of the present application described below are merely some exemplary embodiments of the present application and are not to be construed as limiting the present application.

[0031] FIG. 1 schematically illustrates a local shimming system for magnetic resonance imaging according to an embodiment of the present application. In the present embodiment, the local shimming system includes a multi-channel local shim coil unit 10, a DC power system 106, and a computer control system 104. The local shim coil unit 10 can be installed on an inspection table of the magnetic resonance imaging system 11 and connected to the DC power system 106 via a DC power line 108 through a waveguide plate and a waveguide tube 107. The role of the waveguide plate and the waveguide tube is to reduce noise. The local shim coil unit 10 includes a local multi-channel shim coil 109 and a radio frequency receiving coil 110 for receiving magnetic resonance signals. The RF receiving coil 110 can be placed inside the local multi-channel shim coil 109 with a distance ranging from 1 cm to 2 cm. The DC power system 106 is communicably connected to the computer control system 104. In some embodiments, computer control system 104 may be configured to install and set the software system controlled by the DC power and calculate the field maps and calculate optimization processes. The network cable 105 is configured to connect the computer to the DC power system 106.

[0032] FIG. 2 is a schematic perspective structural diagram of the local shim coil unit according to an embodiment of the present application. In the present embodiment, the local shim coil 202 is fixed on the housing 201. The housing 201 may be a three-dimensional (3D) printed resin housing having a semi-cylindrical configuration and provided with a certain number of holes 203 in order to tune and match the RF receiving coil. As shown in FIG. 2, according to this embodiment, the number of local shim coils 202 is 5, the local shim coils are circular loops having a diameter of 4 cm, and the local shim coils are asymmetrically distributed. The bottom of the resin housing is further provided with a hole 205 having a diameter of 6 cm to facilitate ultrasonic ablation treatment in combination with magnetic resonance imaging techniques, and a positioning groove 204 for positioning the radio frequency receiving coil.

[0033] FIG. 3 is a flow chart of a shimming method for magnetic resonance imaging according to an embodiment of the present application. In the present embodiment, the optimization target of the shimming is to minimize the standard deviation value of the main magnetic field offset, and the average value of the optimized main magnetic field offset is smaller than the average value of the main magnetic field offset before the optimization. In step 301, the B0 field map information is collected using a two-dimensional (2D) multi-echo GRE (gradient echo). In some embodiments, a field map is collected for each channel shimming current at 100 mA and 0 mA, and a sensitivity mapping for each channel coil is obtained after they are subtracted. In the collecting process, the number of echoes in the pulse sequence is set to 5 to reduce the effect of eddy currents on the B0 field map. In order to quickly collect field map information, the pulse sequence has a repetition time (TR) of 25 ms and the pulse flip angle is set as 10. The homogeneity of B0 is calculated and evaluated in step 302. In some embodiments, five echo phase maps are unwrapped, and then straight-line fitting of the pixel points at the same position of the five phase maps is performed on the five corresponding echo times TE by the least square method, and the slope value is the main magnetic field offset value (in Hz) at this position. From the calculation formula

[00003]${.Math.B}_0=\frac{}{.Math..Math.\mathrm{TE}}$

and the formula .sub.0=.Math.B.sub.0 of the main magnetic field offset, the calculation formula of this step can be obtained as follows:

[00004]$.Math.f=\frac{}{2.Math..Math..Math.\mathrm{TE}},$

wherein is the phase difference between the two echoes, is the gyromagnetic ratio of the imaging nucleus, TE is the time difference between the two echoes, and .sub.0 is the nuclear magnetic resonance angle frequency.

[0034] In step 303, the current of each channel shim coil element is optimized. In some embodiments, the main magnetic field offset value f.sub.0 and the sensitivity field map f.sub.0i (i=1, 2, . . . , n, wherein n is the channel number of the shim coil) of each channel shim coil element are linearly combined, which is expressed as:

[00005]$.Math.f={.Math.f}_0+\underset{i=1}{\overset{n}{.Math.}}.Math..Math.a_i.Math.{.Math.f}_{0.Math.i},$

[0035] wherein f.sub.0 is the main magnetic field offset value without applying the shimming current, f.sub.0i is the sensitivity field map of the shim coil, ai is the linear superposition factor of the sensitivity field map, and f is the main magnetic field offset value after the shimming current is applied; through the above process, with the goal of obtaining the minimum standard deviation value of f to optimize the current value of each channel coil, when the average value of the optimized f is smaller than the average value of f.sub.0 prior to optimization, it is effective shimming.

[0036] In step 304, it is determined whether a constraint condition is met. If not, operations 303-304 are repeated until the constraint condition is met. In step 305, an optimal current combination value is output, and optimum current values corresponding to every channel shim coil are set on the current control software. Finally, in step 306, the homogeneity of the main magnetic field B0 is tested and evaluated. The average value and the standard deviation value of the main magnetic field offset value are used as the evaluation index. The smaller the value is, the more uniform the main magnetic field B0 is, and the shimming goal is achieved.

[0037] The present application is experimentally verified on a Siemens 3T magnetic resonance imaging system, and a local shimming system and a shimming method are applied in magnetic resonance temperature measurement imaging. FIGS. 4A and 4B show the comparison of the obtained main magnetic field offset distributions after shimming and before shimming (i.e., in the cases of applying the optimally combined shimming current and applying no shimming current), respectively. Before and after shimming, the average value of the main magnetic field offset values decreases from 16.04 Hz to 9.48 Hz, and the standard deviation value decreases from 5.52 Hz to 3.80 Hz. The result shows that the homogeneity of the main magnetic field is obviously improved, wherein ROI represents the region of interest. FIGS. 5A and 5B show statistics of the magnetic field frequency distributions within the region of interest ROI after shimming and before shimming, respectively, wherein the more concentrated the frequency distributions, the higher the homogeneity of the main magnetic field B0. FIGS. 6A and 6B show the temperature change graph in the phantom experiments at normal temperature conditions after shimming and before shimming, respectively. Within 100s, the average value of the temperature change before and after shimming decreases from 1.01 degrees Celsius to 0.74 degrees Celsius, and the standard deviation value decreases from 1.42 degrees Celsius to 1.32 degrees Celsius. The result shows that the temperature measurement accuracy after shimming is improved. The above experimental results show that the local shimming system and the shimming method produced by the present application can improve the homogeneity of the main magnetic field B0 and improve the accuracy of magnetic resonance temperature measurement at the same time.

[0038] The present application provides a 5-channel local shim coil for magnetic resonance temperature imaging, which compensates for a change in local magnetic field due to the difference in magnetic proton susceptibility between water and fat tissue, thereby improving the precision and accuracy of temperature measurement. The magnetic resonance shimming system and the shimming method of the embodiment of the present application achieve a better shimming effect using a shim coil with a smaller number of channels, and the system is simpler and the cost is relatively low; and the shimming method provided is relatively simple and is relatively short in time; at the same time, the system is applicable in magnetic resonance temperature imaging, which improves the accuracy of temperature measurement.

[0039] Although the present application has been described in terms of the preferred embodiments, modifications, replacements, and various alternatives are possible that fall within the scope of the present application. It should also be noted that there are many alternative ways of implementing the method and system of the present application. Accordingly, it is intended that the appended claims be interpreted as including all such modifications, replacements, and various alternatives that fall within the gist and scope of the present application.