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 method for magnetic resonance imaging, comprising: collecting a field map information of a main magnetic field B0 using a two-dimensional gradient echo; calculating and evaluating a homogeneity of the main magnetic field B0; determining an offset value Δf of the main magnetic field B0 which is generated after a current of each channel shim coil is applied optimizing the current of each channel shim coil; determining whether a minimum standard deviation value of the offset value Δf is obtained; outputting an optimal current combination value and setting an optimum current value corresponding to each channel shim coil on a current control software when the minimum standard deviation value of the offset value Δf is obtained; and testing and evaluating the homogeneity of the main magnetic field B0 to achieve a shimming goal, wherein the optimizing the current of each channel shim coil comprises: determining an offset value Δf.sub.0 of the main magnetic field B0 which is generated without applying the current of each channel shim coil and a sensitivity field map Δf.sub.0i of each channel shim coil; linearly combining the offset value Δf.sub.0 of the main magnetic field B0 which is generated without applying the current of each channel shim coil and the sensitivity field map Δf.sub.0i of each channel shim coil according to the following equation: $\Delta f=\Delta f_0+\underset{n=1}{\overset{n}{.Math.}}\left(a_i.Math.\Delta f_{0i}\right)$ wherein i=1, 2, . . . , n, and n is the number of channels of the shim coil, and a.sub.i is a linear superposition factor of the sensitivity field map of each shim coil; and optimizing the current of each channel shim coil through minimizing the standard deviation value of the offset value Δf.

2. The shimming method of claim 1, wherein the collecting comprises collecting the field map information of the main magnetic field B0 with a pulse sequence echo number of 5, a pulse sequence repetition time of 25 ms, and a pulse flip angle of 10°.

3. The shimming method of claim 2, wherein the calculating comprises unwrapping five echo phase maps, and performing straight-line fitting of pixel points at a same position of the five echo phase maps on five corresponding echo times TE by a least square method, wherein both the offset value Δf and Δf.sub.0 of the main magnetic field B0 at this position are a slope value which is determined through:
Δϕ/2π.Math.ΔTE wherein Δϕ is a phase difference between two echoes and ΔTE is a time difference between the two echoes.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) 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:

(2) FIG. 1 is a schematic structural diagram of the local shimming system for magnetic resonance imaging according to an embodiment of the present application.

(3) FIG. 2 is a schematic perspective structural diagram of the local shim coil unit according to an embodiment of the present application.

(4) FIG. 3 is a flow chart of the shimming method for magnetic resonance imaging according to an embodiment of the present application.

(5) FIGS. 4A and 4B show the comparison of the obtained main magnetic field offset distributions after shimming and before shimming, respectively.

(6) FIGS. 5A and 5B show statistics of the magnetic field frequency distributions within the region of interest ROI after shimming and before shimming, respectively.

(7) FIGS. 6A and 6B show the temperature change graph of a phantom at normal temperature conditions after shimming and before shimming, respectively.

(8) 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

(9) 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.

(10) 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.

(11) 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.

(12) 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

(13) $\Delta B_0=\frac{\Delta \varphi }{\gamma .Math.\Delta \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:

(14) $\Delta f=\frac{\Delta \varphi }{2\pi .Math.\Delta \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.

(15) 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:

(16) $\Delta f=\Delta f_0+\underset{i=1}{\overset{n}{.Math.}}{a}_i.Math.\Delta f_{0i},$

(17) 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, a.sub.i 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.

(18) 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.

(19) 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 100 s, 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.

(20) 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.

(21) 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.