CHEMICAL SPECIES SUPPRESSION USING RE-TRAVERSAL AROUND THE K-SPACE CENTER IN INTERLEAVED REVERSED SPIRAL IMAGING
BACKGROUND OF THE INVENTION The present invention relates to a method of differentiating a chemical species in Magnetic Resonance Imaging (MRI), and more particularly to a method of suppressing a chemical species while reducing acquisition times in MRI.
A number of suppression techniques have been implemented to differentiate chemical species in MRI imaging. For example fat suppression techniques are used in MRI to differentiate fat from water in an imaged object.
Techniques such as Dixon methods achieve suppression by using additional out-of-phase k-space acquisitions. However acquiring the additional k-space data can create an undesirable increase in acquisition times.
Techniques such as binomial, STIR, CHESS techniques use specialized Radio
Frequency RF excitation schema to selectively excite water protons or saturate fat protons. Unfortunately, using specialized RF excitation pulses also requires additional acquisition time. Similarly, spatially spectrally selective excitation RF pulses used for suppression in both forward and reversed spiral imaging are relatively long duration pulses which also cause undesirable increases in scan time.
Longer acquisition times limit the use of these techniques for fast or real- time MR applications such as cardiac imaging. Therefore, it is desirable to provide an MR pulse sequence that enables suppression of a chemical species using relatively short acquisition times.
SUMMARY OF THE INVENTION According to the present invention, a new and improved method of generating an MRI image by selectively suppressing one of two chemically shifted species is provided.
In accordance with a first aspect of the invention, the method includes generating a first spiral trajectory in k-space which acquires the center of the spiral, generating a second spiral trajectory in k-space retraversing a portion of the first spiral trajectory to the center after a delay, wherein the delay is the time required for the first chemical species to rotate out of phase with a second chemical species in the subject. The invention further includes combining the first and second spiral trajectories to suppress the signal of the first chemical species in the portion of k-space defined by the second spiral trajectory, and repeating the steps of generating the first and second spiral trajectories to form a plurality of trajectory interleaf pairs each covering a different portion of k- space.
In accordance with a second aspect of the invention, the invention includes generating an image using the trajectory interleaf pairs in which the first chemical species is suppressed relative to the second chemical species. Other features, benefits and advantages of this invention will become apparent to those skilled in the art from the following detailed description of the preferred embodiments, when read in light of the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
The invention may take form in certain components and structures, preferred embodiments which will be illustrated in the accompanying drawings wherein:
Fig. 1 illustrates steps of the invention; Fig. 2 is a diagram illustrating the second spiral trajectory re-traversing the first spiral trajectory in k-space in accordance with the invention; and
Fig. 3 is a diagram illustrating the use of several trajectory interleaf pairs.
DETAILED DESCRIPTION OF THE INVENTION Referring now to Fig. 1 a method of suppressing a first chemical species in MR imaging of a subject having first and second chemical species is shown generally at 10. In the example described herein, the first chemical species is
fat and the second chemical species is water. Although, alternatively, the first chemical species can be water and the second chemical species can be fat, or the first and second chemical species can be any other suitable chemical species exhibiting a relative difference in precession frequencies during MRI acquisitions.
The method includes generating a first spiral trajectory in k-space which acquires the center of the spiral at 12. The first spiral trajectory, shown at 22 in Fig. 2, can be generated in any suitable known manner. The first spiral trajectory 22 is shown to acquire low spatial frequencies of k-space at the center of the spiral 24. The center of the spiral 24 is located at the center, also referred to the origin, of k-space shown at 26.
The method further includes at 14 in Fig. 1 , generating a second spiral trajectory in k-space after a delay which retraverses a portion of the first spiral trajectory 22 to the center of k-space 26. The delay is the time required for the first chemical species to rotate out of phase with a second chemical species present in the subject. The delay is preferably long enough for the first chemical species to rotate 180 degrees out of phase relative to the second chemical species. In the fat suppression example shown herein, a delay 2.2msec was used. The second spiral trajectory, shown at 28 in Fig. 2, can also be generated in any suitable known manner. The first and second trajectories 22, 28 were generated using gradient waveforms with TR/TE(first)/TE(second)= 17.0/8.5/10.7msec and a 10° tip angle on a Siemens 1.5 T Magnetom Sonata (Siemens, Erlangen Germany), although any other suitable parameters can be used.
The first and second spiral trajectories 22 and 28 are preferably a reverse spiral trajectories, although alternatively the first and second spiral trajectories can both be forward spiral trajectories. The second spiral trajectory 28 is sized smaller than the first spiral trajectory 22 so it re-traverses only a portion of the first spiral trajectory, the re-traversed portion shown at 29. For example, in the example described herein, the radius, shown at r, of the second spiral trajectory 28 is approximately 30% of kmax of the first spiral trajectory, although the second
spiral trajectory can have any other suitable radius which is smaller than the radius of the first spiral trajectory.
The method further includes at 16 in Fig. 1 , combining the first and second spiral trajectories to suppress the signal of the first chemical species in the portion of k-space defined by second spiral trajectory 28. The reversed spiral trajectories 22, 28 incorporate first moment gradient nulling at the start of the spiral and re-traversal of the center of k-space 2.2 msec after the center was initially acquired. The re-traversal goes back to a point of the original trajectory 30 and follows the same trajectory as the first spiral 22 back to the origin 26. The time difference, or delay, between the first and second traversals is the time for the first chemical species (in this example fat having a 3.5ppm chemical shift) to rotate out of phase with water by 180 degrees. Therefore, the fat signal of the first and the second traversals cancel each other in the lower frequency region of k-space when combined in step 16. The method further includes at 18 in Fig. 1 , repeating the steps of generating the first spiral trajectory 14 and generating the second spiral trajectory 16 to form a plurality of trajectory interleaf pairs each covering a different portion of k-space. The gradient waveforms of the interleaf pairs were incorporated into a 12-interleave reversed spiral sequence to achieve a sufficiently high sampling density of k-space. Since each spiral interleaf reaches the k-space origin twice, two bipolar gradients were incorporated for flow compensation. Though 12-interleaf trajectory pairs were used in the example described herein, any suitable number of interleaf trajectory pairs can be used to sample a sufficient amount of k-space. For example, Fig. 3 shows the use of 3 interleaf trajectory pairs, 22a and 28a, 22b and 28b, and 22c and 28c, each covering a different portion of k-space.
The method further includes at 20 in Fig. 1 , generating an MR image using the trajectory interleaf pairs in which the first chemical species is suppressed relative to the second chemical species. The MR image can be generated in any suitable know manner using the trajectory interleaf pairs.
All the MR images were reconstructed using a modified known Block Uniform Resampling (BURS) regridding algorithm, although any other suitable
known methods of image reconstruction can be used. The slice thickness was 10mm and the Field of View (FOV) was 210mm x 210mm, although any suitable slice thickness and FOV can be used.
A short duration re-traversal around the center of k-space after the spiral trajectory reaches the center of k-space has been shown to be an effective fat suppression technique in reversed spiral imaging. The total scan time is only minimally affected since each spiral interleaf needs only 2.2ms of additional readout time for each spiral interleaf in a 1.5T MR system. In a 12 interleave spiral sequence, this results in a 26.4 msec increase (2.2ms x 12 interleaves) in scan time for fat suppression through re-traversal, versus about 60 msec additional scan time for a 5 msec spatial spectral pulse (using a 5 msec pulse with 12 interleaves).
The invention has been described with reference to preferred embodiments. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding specification. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.