1. Field of the Invention
The present invention relates to a method and an apparatus for magnetic resonance (MR) system, particularly for simultaneously acquiring multi-slice/slab magnetic resonance imaging (MRI) signals.
2. Description of the Prior Art
a. Principle of Two-Dimensional Magnetic Resonance Imaging
The principle of two-dimensional (2D) MRI procedure is described as follows: As a subject is placed in a static magnetic field, a region of the subject can be excited by using a radio-frequency (RF) coil and giving signals with respect to all the excitation and relaxation of nucleus excitations and relaxations in the region. With a (magnetic) gradient applied, the RF coil can receive those signals, which can be processed to a MR image. If the change in the structure or functionality of the region is to be realized, the gradient may be adjusted so that slices can be acquired from various locations in the region.
In the following, 2D spatial encoding is described. In the beginning, a slice-selection gradient GZ is turned on and a slice normal to the Z-direction is excited by RF pulses with suitable frequencies. A phase-encoding gradient GY is then turned on for a period of time and then is turned off, so that those nuclei have a certain phase difference in the Y-direction. As a result, a frequency-encoding gradient GX is turned on, while starting to receive signals. Due to this phase difference, the sum of signals with different frequencies (in the X-direction) is received, giving a line in a space spanned by two variables, i.e., the phase difference and the frequency. This space is termed the k-space. A similar procedure may be repeated, where the phase-encoding gradient GY is changed and thus those nuclei have another certain phase difference in the Y-direction. At this phase difference, the sum of signals with different frequencies (in the X-direction) is received, giving another line in the k-space. The whole sampling in the k-space is completed as the sums of signals with different frequencies are received at different phase differences. When the phase difference is changed, the slice-selection gradient GZ is turned on to excite the selected slice.
A 2D Fourier transform is performed after the k-space sampling is completed. That is, the phase differences and the frequencies are transformed into signal intensities at locations on the XY-plane, forming an image of a horizontal slice or a slab. Thus, a 2D MR image is made.
b. Principle of Three-Dimensional Magnetic Resonance Imaging
The principle of the 3D MRI is similar to that of the 2D MRI, except for a difference in the spatial encoding. A phase-encoding gradient GY and a slab-selection gradient GZ are turned on for a period of time and then are turned off, so that those nuclei have certain phase differences in the Y-direction and Z-direction. A similar procedure may be repeated, where the phase-encoding gradient GY and the slab-selection gradient GZ are changed. A 3D Fourier transform is performed after the k-space sampling is completed, giving a 3D MR image.
c. Common Drawback of Two or Three-Dimensional Magnetic Resonance Images
Referring to
Time required for acquiring the images of all the slices
=NEX×Npe×TR×Nslice, (Eq. 1)
where NEX denotes the average number of signaling and Npe is the whole number of encoding. For the 2D MR image, Npe denotes the number of phase encoding Np, TR denotes the time required for acquiring a line in the k-space, and Nslice denotes the number of the slices. For example, if there are total 256 images to be acquired, then Nslice=256, NEX=1, Npe=128, TR=0.1 second, and the time required for acquiring the images of all the slices is about 54 minutes. This is therefore a time-consuming procedure.
Referring to
d. Other Related Techniques
MRI is a useful tool for biomedical applications to obtain real-time images. Any possible method to accelerate the MRI scan time is highly attractive. Thus, a great deal of manpower and resources have been invested in this research field, resulting in the development of various ways of acceleration such as simultaneous excitation-time division multiple acquisition, phased array coil acceleration, and reduction in data reception.
The following list some of the approaches developed to date:
i. Simultaneous multi-slice acquisition of MR images by Hadamard encoded excitation (called SIMA sequence for short). Multiple slices are excited simultaneously, but a composite image of the multiple slices is received. The characteristic of the RF electromagnetic wave excited from each slice is required for solving a certain polynomial to give the complete information of an individual image. This technique has the disadvantage that extra time for computation may be needed and that N excitations and N receptions are required for computing the information of N slices.
ii. Simultaneous parallel inclined readout image technique (called SMA sequence for short). In this approach, multiple slices are excited simultaneously, and in receiving, different magnetic field intensities are applied to different slices by means of the gradient coil. This technique has the disadvantage that the additional gradient coil may result in inevitable blur in the image and may thus reduce the image quality.
iii. Use of multi-coil array for the separation of signals from multiple slices simultaneously excited (called SENSE sequence for short). In this approach, the images of different slices are received simultaneously using the sensitivity difference of multiple coils at different locations. The real images with respect to different locations are then computed. This technique has the disadvantage that extra coils are needed; for example, four times acceleration requires at least four coils. The acceleration effect is not proportional to the amount of coils.
iv. MAMBA. One stepped gradient coil is added, besides common linear gradient coils, so that the images of different slices have different resonant frequencies. This technique has the disadvantage that extra coils are needed (MAMBA gradient coils) and the acceleration factor and the multi-slice pitch are unchanged due to the fact that the coils cannot be arbitrarily adjusted.
Therefore, whether simultaneous excitation-time division multiple acquisition or phased array coil acceleration is used, conventional multi-slice image processing techniques cannot, at several times acceleration, simultaneously excite and acquire multi-slice, real-time images under the same pulse time-sequence without adding coils or computation algorithms, or reducing the acquired data. Moreover, these conventional multi-slice/slab MR image processing techniques require additional costs for adding coils or computation algorithms.
A method for simultaneously acquiring multi-slice/slab MRI signals from a subject is describer herein.
In one aspect, a method for simultaneously acquiring multi-slice/slab MRI signals comprises steps of: (a) applying one or more than one RF pulse, which carries at least two frequency components, and a slice/slab selection gradient so that at least two slices/slabs of the subject respectively corresponding to the at least two frequency components are excited simultaneously; (b) applying a spatial encoding gradient; and (c) applying a slice/slab separation gradient for separating at least two slices/slabs.
In another aspect, an apparatus for simultaneously acquiring multi-slice/slab MRI signals from a subject is capable of imaging by generating a MRI signal from the subject and by reducing the MRI signal to spatial encoding data with respect to each slice/slab of the subject. The apparatus comprises an RF excitation module, an RF receiving module, a gradient output module, and a sequence controller for controlling the activation of those modules.
The sequence controller may be used for controlling the RF excitation module, the gradient output module and the RF receiving module, so as to perform the method for simultaneously acquiring multi-slice/slab MRI signals. It comprises steps of controlling the subject.
The RF excitation module may be controlled and used for applying to the subject one or more than one RF pulse, which carries at least two frequency components, so that at least two slices/slabs of the subject respectively corresponding to at least two frequency components are excited simultaneously. The gradient output module may be controlled and used for applying to the subject of a spatial encoding gradient, a slice/slab selection gradient, and a slice/slab separation gradient for separating the at least two slices/slabs. The RF receiving module may be controlled and used for receiving the MRI signal excited from the subject.
In still another aspect, a method for simultaneously acquiring multi-slice/slab MRI signals from a subject, the method being capable of imaging a slice/slab of a subject, comprises the steps of: (a) applying one or more than one RF pulse, which carries at least two frequency components, and a slice/slab selection gradient so that at least two slices/slabs of the subject respectively corresponding to at least two frequency components are excited simultaneously; (b) applying a spatial encoding gradient; and (c) applying a slice/slab separation gradient, which may be used for separating at least two slices/slabs, and receiving the MRI signal excited from the subject; and (d) reconstructing the MRI signal by performing spatial encoding and 2D Fourier transform, so as to give the real-time image of each slice/slab.
These and other features, aspects, and embodiments are described below in the section entitled “Detailed Description.”
Features, aspects, and embodiments are described in conjunction with the attached drawings, in which:
In the following description, at first it will be noted that similar elements are denoted by the same numbers and, for the purpose of convenience, the symbol “/” means “of”. In addition, the coordinate (X, Y, Z) is used with respect to the coordinate of images but is not an absolute spatial coordinate about the MRI system.
The methods and apparatus described herein make it possible to simultaneously excite and acquire MR images of different locations in the subject. These embodiments are compatible with a variety of existing MRI systems. Moreover, they do not require extra coils/RF channels and extra time for computation of image information and extra computer equipment. The embodiments described herein can be applied to MRI systems by means of echo planar imaging, perfusion, image flow, angiogram, image temperature, T1 imaging (lattice-spin relaxation time constant), T2 imaging (spin-spin relaxation time constant), diffusion and the like.
1. Apparatus:
Referring to
A subject 6 can be positioned within a measurement space 30. In the measurement space 30, there are a uniform magnetic field generated by the static magnetic field output module 3 and a gradient generated by the gradient coils 42, which is under the control of the gradient controller 41. The uniform magnetic field cooperating with the gradient is used for the subject 6 to generate the magnetization as a source of the MRI signal.
2. Methods for Imaging Using the Above Apparatus:
a. 2D MRI:
Referring to
For example, the RF pulse 101 carrying frequencies f1 and f2 for exciting two slices has f1 and f2 designed to have a frequency difference fsep such that
f
sep
=d
sep
×γ×G
ss. (Eq. 2)
dsep denotes the absolute distance (in centimeters) between the two adjacent slices that correspond to the two frequencies, wherein two adjacent slices mean two neighboring slices excited by a plurality of frequencies carried by an RF pulse. γ denotes the atomic gyromagnetic ratio. Gss denotes the intensity of the slice-selection gradient (in Gauss/centimeter).
In step 202, the gradient output module 4 applies to the subject 6 a spatial encoding gradient Gspen and at least one slice separation gradient Gsep for separating at least two slices, wherein the spatial encoding gradient Gspen comprises a phase-encoding gradient GY and a frequency-encoding gradient GX.
The ratio of Gsep to Gspen, which may be the phase-encoding gradient GY or the frequency-encoding gradient GX and cooperates with Gsep to receive MRI signals, should meet a relation as follows:
G
sep
/G
spen≧FOVspen/dsep. (Eq. 3)
FOVspen denotes a field (e.g. width) (in centimeters) of view along the direction of the spatial encoding gradient Gspen. dsep denotes the absolute distance (in centimeters) between the two adjacent slices that correspond to the two frequencies. The two adjacent slices can be separated completely as held by Eq. 3.
In step 203, the gradient output module 4, in the course of applying the spatial encoding gradient (phase-encoding gradient GY and a frequency-encoding gradient GX), cooperates with the RF receiving module 22 to receive the MRI signals excited from the subject 6.
In step 204, the MRI signal is reconstructed by performing spatial encoding and 2D Fourier transform, so as to give the real-time image of each slice.
b. 3D MRI:
The difference between 3D MRI and 2D MRI is that in the latter, one slice is excited in one scan and then the image information is given through 2D(?) spatial encoding, whereas in the former one slab is excited in one scan and then the image information is given through 3D spatial encoding.
Referring to
For example, the RF pulse 101 carrying frequencies f1 and f2 for exciting two slabs has f1 and f2 designed to have a frequency difference fsep as held by Eq. 2.
In step 302, the gradient output module 4 applies to the subject 6 a spatial encoding gradient Gspen and at least one slab separation gradient Gsep for separating the at least two slabs, wherein the spatial encoding gradient Gspen comprises a phase-encoding gradient GY, a frequency-encoding gradient GX, and a slab-selection gradient GZ.
The ratio of Gsep to Gspen should obey Eq. 3 which may be the phase-encoding gradient GY, the frequency-encoding gradient GX or the slab-selection gradient GZ and cooperates with Gsep to receive MRI signals.
In step 303, the gradient output module 4, in the course of applying the spatial encoding gradient, cooperates with the RF receiving module 22 to receive the MRI signals excited from the subject 6.
In step 304, the MRI signal is reconstructed by performing spatial encoding and 3D Fourier transform, so as to give the real-time image of each slab.
3. Control Procedure:
a. 2D MRI:
Referring to
The sequence controller 1 can be configured to output driving signals to drive the activation of the RF excitation module 21, the RF receiving module 22 and the gradient output module 4, wherein the steps of controlling, in accordance with one embodiment, are as follows.
i. The sequence controller 1 drives the RF excitation module 21 to generate to the subject 6 one or more than one RF pulse 101 that carries at least two frequency components. For example, the RF pulse 101 carries frequencies f1 and f2, which have a difference fsep meeting Eq. 2. Yet, the sequence controller 1 drives the gradient controller 41 to control the gradient coils 42 for generating the slice-selection gradient Gss 104.
ii. The sequence controller 1 applies a spatial encoding gradient Gspen to the subject 6, with respect to each encoding direction for each slice.
iii. The sequence controller 1 drives the gradient controller 41 to control the plurality of gradient coils 42 for generating to the subject 6 at least one slice separation gradient Gsep 105.
In this example, the spatial encoding gradient Gspen comprises a phase-encoding gradient 102 and a frequency-encoding gradient 103 and, in the course of applying the spatial encoding gradient Gspen, at least one slice separation gradient Gsep 105 is applied while the MRI signals excited from the subject 6 are received. The ratio of Gsep 105 to Gspen, which cooperates with Gsep to receive the MRI signals, should meet Eq. 3 so that two adjacent slices can be separated completely.
Depending on the embodiment, the receiving module 22 can receive the MRI signals excited from the subject 6 while the frequency-encoding gradient 103 and the at least one slice separation gradient 105 are applied. The MRI signals can then be reconstructed by the image processing module 73 performing transformation such as spatial encoding and 2D Fourier transform, so as to give the real-time image data of the slices. The reconstructed data, which are the images of the separated slices, can then be output to be displayed on a display device 52.
b. 3D MRI:
Referring to
i. The sequence controller 1 drives the RF excitation module 21 to generate to the subject 6 one or more than one RF pulse 101 that carries at least two frequency components; for example, the RF pulse 101 carries frequencies f1 and f2, which have a difference fsep meeting Eq. 2. Yet, the sequence controller 1 drives the gradient controller 41 to control the gradient coils 42 for generating the slab-selection gradient Gss 104.
ii. The sequence controller 1 applies a spatial encoding gradient Gspen to the subject 6, with respect to each encoding direction for each slab.
iii. The sequence controller 1 drives the gradient controller 41 to control the plurality of gradient coils 42 for generating to the subject 6 at least one slab separation gradient Gsep 105.
In this example, the spatial encoding gradient Gspen comprises a phase-encoding gradient 102, a frequency-encoding gradient 103 and a slab-selection gradient 106. In the course of applying the spatial encoding gradient Gspen, at least one slab separation gradient Gsep 105 is applied while the MRI excited signals from the subject 6 are received. The ratio of Gsep 105 to Gspen, which cooperates with Gsep to receive the MRI signals, should meet Eq. 3 so that two adjacent slabs can be separated completely.
Depending on the embodiment, the receiving module 22 can receive the MRI signals excited from the subject 6 while the frequency-encoding gradient 103 and at least one slab separation gradient 105 are applied. The MRI signals can then be reconstructed by the image processing module 73 performing spatial encoding and transformation such as 3D Fourier transform, so as to give the real-time image data of the slabs. The reconstructed data, which is the images of the separated slabs, can then be output to be displayed on a display device 52.
4. Reduction of Blur to Give a Sharp Image:
a. Cause of Image Blur:
Referring to
b. Method for Overcoming Image Blur:
General equations for calculating the image blur are as follows.
blur(mm)=resZ×GZ/GX, (Eq. 4-1)
blur(pixel)=blur(mm)/resX. (Eq. 4-2)
When Eq. 4-1 and Eq. 4-2 are applied to 2D MRI, resZ denotes the thickness of one slice, resX denotes the resolution in the X-direction, GZ denotes the slice-selection gradient, and GX frequency-encoding gradient.
When Eq. 4-1 and Eq. 4-2 are applied to 3D MRI, the only difference lies in that resZ denotes the ratio of the thickness of the slab to the number of spatial encoding in the Z-direction. The other parameters are the same as those for 2D MRI.
Now set forth Eq. 3 again, which describes the criterion for complete separation of two adjacent slices/slabs.
G
sep
/G
spen≧FOVspen/dsep (Eq. 3)
With the appropriate substitutions, Eq. 3 can be combined with Eq. 4-1 and Eq. 4-2 so that;
resZ=[blur(pixel)×dsep×resX]/FOVX, (Eq. 5)
where resZ denotes the resolution of the image in the Z-direction, blur is the degree of image blur (in unit of pixel), dsep denotes the absolute distance between the two adjacent slices/slabs, resX denotes the resolution in the X-direction, and FOVspen denotes a field (e.g. width) of view along the direction of the spatial encoding gradient.
According to Eq. 5, the blur is required to be set in a preset range so that a sharp, unblurred, image may be given. The following is an example, which illustrates how to set the imaging conditions into practice. Assume the imaging conditions to be as follows: blur=3 pixels; dsep=1.67 cm; resX=1 mm; FOVspen=10 mm. Replacing these parameter values in Eq. 5, one has that resZ=0.5 mm. In other words, resZ must be determined to be less than 0.5 mm so that a sharp image can be given.
The preset blur can be set according to practical needs. For example, there is no change in the quality of image as the blur is less than 1 pixel, and therefore the image is of very good quality. The change is difficult to view from the naked eye as the blur is in the range of 1-3 pixels; making an image of good quality. When the blur is 3 pixels up to 6 pixels, the quality of image is of poor quality. When the blur is more than 6 pixels, the quality of the image is bad and many features become hard to determine. Nevertheless, adjustment may be made according to various needs since different blurs may find their uses in different situations.
For example, with the imaging conditions set forth above, the resultant images of a SMA series gives inevitable blur due to the additional applied gradients, whereas the techniques described herein reduce the blur and therefore a sharp image of the slice/slab is acquired simultaneously.
5. Experimental Results:
Referring to
Referring to
Similarly, referring to
Simultaneous acquisition of a MRI signal from 8 slices is possible, as well as 2 slices and 4 slices described above. As shown in
Thus, the systems and methods described herein can provide the following:
1. High imaging speed: The systems and methods described herein can be used to simultaneously excite and acquire multiple MR images of different locations in the subject. Thus, they can save the time for imaging of MRI systems.
2. High compatibility: The systems and methods described herein are compatible with a variety of existing MRI systems. Moreover, they do not require extra coils/RF channels and extra time for computation of image information and extra computer equipment. Thus, the systems and methods described herein can promote the efficiency of MRI systems.
3. High image sharpness: The systems and methods described herein can be used to simultaneously give multiple sharp images by utilizing techniques of high-resolution slice/slab acquisition.
While certain embodiments have been described above, it will be understood that the embodiments described are by way of example only. Accordingly, the systems and methods described herein should not be limited based on the described embodiments. Rather, the systems and methods described herein should only be limited in light of the claims that follow when taken in conjunction with the above description and accompanying drawings.
Number | Date | Country | Kind |
---|---|---|---|
097116790 | May 2008 | TW | national |
097130053 | Aug 2008 | TW | national |