Abstract:
A method of inducing tag lines in an MRI image is provided which allows adjustment of tag line width and spacing. The method includes determining a desired tag line width ( 505 ) and then selecting a parameter of a SINC modulating function, such as main lobe bandwidth, in accordance with the desired tag line width. The method also includes determining a desired tag line separation ( 51 ) and selecting a parameter of an ARF pulse train based upon the desired tag line separation. A gradiant field is applied to the specimen being imaged. The RF pulse train is modulated in accordance with the SINC modulating function ( 520 ) and the modulated RF pulse train is applied to the specimen being imaged prior to an image acquisition operation. In one embodiment of the invention, the gradiant is activated between pulses of the RF pulses.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates generally to magnetic resonance imaging and more particularly relates to magnetic resonance imaging using spatial modulation of magnetization to establish a frame of reference in the form of tag lines on the resulting image.  
       BACKGROUND OF THE INVENTION  
       [0002]     Magnetic resonance imaging (MRI) is a well established tool for medical imaging. In magnetic resonance imaging (MRI), the proton density of physical tissue characteristics is sensed to discriminate between healthy and diseased tissue. In the case of tissue that is in motion, such as the human heart, diseased tissue of the heart can be detected by analysis of wall motion of a heart chamber and wall thickness during the cardiac cycle.  
         [0003]     U.S. Pat. No. 4,953,554 describes a system and method for MRI imaging which uses a modulated RF signal to impose reference void lines on the resulting image. The &#39;554 patent recognizes that such non-invasive MRI tagging is useful in establishing a frame of reference in the image data which is important in the analysis of moving bodies, such as the heart. The method described in the &#39;554 patent includes modulating an. RF pulse that is applied to the body being imaged prior to acquiring image data. More particularly, the &#39;554 patent discloses modulating an RF carrier pulse in the form (1+2 cos ω 0 t) with a function in the form of (sin(t)/t) just prior to conducting imaging acquisition. However, this form of modulated signal results in a non-uniform gradient which can result in aliasing.  
         [0004]     There have been a number of improvements developed to spatially modulate the magnetization field, thereby providing a reference frame on the image in a non-invasive manner. Most notable among these improvements are the SPAMM and DANTE tagging techniques.  
         [0005]     The SPAMM method, which is described for example in the article “A cine-SPAMM Sequence for NMR Imaging of Pulsatile Cerbrospinal Fluid Flow,” by Wayte et. al., Phys. Med. Biol., pp. 455-463, 1993, produces a sinusoidal modulation of saturated magnetization. An improvement to SPAMM, know as high-order SPAMM, improves the resulting modulation profile by employing binomial modulation of the RF pulse train. While these systems have proven effective in a number of applications, neither SPAMM nor high order SPAMM provide the user with the flexibility to alter the ratio of the tag width (width of the void line) to tag separation. Thus, the user is limited in the form of the resulting image tagging that can be achieved.  
         [0006]     DANTE tagging employs a delta function to effect tagging. By effecting spatial modulation with a series of narrow RF pulses of uniform amplitude, the resulting saturation tags exhibit a SINC profile. The DANTE system is further described, for example, in the article entitled “Myocardial Tagging with B1 insensitive Adiabatic DANTE Inversion Sequences” by Tesekos, et al., 1994, pp 395-401, which is hereby incorporated by reference in its entirety.  
       SUMMARY OF THE INVENTION  
       [0007]     It is an object of the present invention to provide MRI tagging using a modulated RF pulse train which yields a rectangular profile of saturated magnetization.  
         [0008]     It is a further object of the present invention to provide a method of MRI tagging which provides flexibility in selecting the resulting tag width and tag separation.  
         [0009]     In accordance with the present invention, a first method for inducing tag lines in an MRI image begins with the steps of determining a desired tag line width and determining a desired tag line separation. One or more parameters of a SINC modulating function are then selected in accordance with the desired tag line width. One or more parameters of an RF pulse train based are selected based upon the desired tag line separation. A gradient field is applied to the specimen being imaged. The RF pulse train is modulated in accordance with the SINC modulating function and the modulated RF pulse train is applied to the specimen being imaged prior to an image acquisition operation.  
         [0010]     Preferably, the SINC modulating function is defined at least in part by the bandwidth of a main lobe, and by altering the bandwidth of the main lobe the width of the resulting tag lines is altered accordingly. It is also preferable for the RF pulses of the RF pulse train to be separated by a delay time between pulses which can be altered to determine the resulting tag line separation. Because these parameters of the modulating waveform and RF pulse train can be varied substantially independently of each other, a great deal of flexibility is provided.  
         [0011]     Each pulse in the RF pulse train has a flip angle associated with a pulse amplitude and pulse duration. Preferably, the pulse amplitude and pulse duration can be varied, in an inverse relationship, to maintain the flip angle associated with the pulse while optimizing or reducing the peak power requirement. In one embodiment, the pulses in the RF pulse train have a common amplitude and the duty cycle of the pulses is modulated to deliver the appropriate flip angle per pulse in accordance with the SINC modulating function.  
         [0012]     In the methods set forth above, the step of applying the gradient preferably occurs in a pulsed manner such that the gradient is active only between pulses of the RF pulse train. 
     
    
     BRIEF DESCRIPTION OF THE DRAWING  
       [0013]     Further objects, features and advantages of the invention will become apparent from the following detailed description taken in conjunction with the accompanying figures showing illustrative embodiments of the invention, in which  
         [0014]      FIG. 1  is a simplified block diagram of an MRI apparatus;  
         [0015]      FIG. 2A  is a timing diagram illustrating an RF magnetization pulse train in accordance with the present invention in the presence of a constant magnetic gradient Gx;  
         [0016]      FIG. 2B  is a graph of the spatial modulation of magnetization, having a tag width to separation ratio of 1:2, which results from the RF pulse train of  FIG. 2A ;  
         [0017]      FIG. 3A  is a timing diagram illustrating an RF magnetization pulse train in accordance with the present invention in the presence of a constant magnetic gradient Gx;  
         [0018]      FIG. 3B  is a graph of the spatial modulation of magnetization, having a tag width to separation ratio of 1:4, which results from the RF pulse train of  FIG. 3A ;  
         [0019]      FIG. 4A  is a timing diagram illustrating an RF magnetization pulse train in accordance with the present invention in the presence of a constant magnetic gradient Gx;  
         [0020]      FIG. 4B  is a graph of the spatial modulation of magnetization, having a tag width to separation ratio of 1:8, which results from the RF pulse train of  FIG. 4A ;  
         [0021]      FIG. 5  is a flow chart illustrating a method of applying a SINC modulated RF pulse train to establish tag lines in a specimen being imaged in accordance with the present invention;  
         [0022]      FIGS. 6A and 6B  are graphs illustrating the timing relationship between the SINC modulated RF pulses and a pulsed gradient signal wherein the gradient pulses are active between the RF pulses;  
         [0023]      FIGS. 7A and 7B  are graphs illustrating the timing relationship between the SINC modulated RF pulses and a pulsed gradient signal in an equivalent relationship to that depicted in  FIG. 6 , wherein the RF pulses have a longer duty cycle and lower peak amplitude;  
         [0024]      FIGS. 8A and 8B  are graphs illustrating the timing relationship between the SINC modulated RF pulses and a pulsed gradient signal in an equivalent relationship to that depicted in  FIG. 6 , wherein the RF pulses have a constant magnitude but exhibit a varying duty cycle determined by the SINC modulation;  
         [0025]      FIGS. 9A, 9B  and  9 C are examples of images acquired during a phantom study using the waveforms of  FIGS. 2A, 3A  and  4 A, respectively;  
         [0026]      FIG. 9D  is an example of an image acquired during a phantom study illustrating an extension of the present invention into two dimensional tag lines;  
         [0027]      FIG. 10  is an example of an image acquired during a phantom study using the prior art DANTE technique for tagging the image;  
         [0028]      FIG. 11  is an example of an image of a mouse heart acquired with the present invention; and  
         [0029]      FIG. 12  is an example of an image of a mouse heart acquired with the prior art DANTE technique for tagging the image. 
     
    
       [0030]     Throughout the figures, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components or portions of the illustrated embodiments. Moreover, while the subject invention will now be described in detail with reference to the figures, it is done so in connection with the illustrative embodiments. It is intended that changes and modifications can be made to the described embodiments without departing from the true scope and spirit of the subject invention as defined by the appended claims.  
       DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS  
       [0031]     The method of the present invention may be used with standard MRI apparatus if such apparatus has control capabilities where parameters of the radio frequency (RF) pulse train and gradient may be easily varied. Referring to  FIG. 1 , a simplified block diagram of an MRI device is illustrated. The MRI device includes a main magnet  10  which provides a magnetic field B 0  that generates a steady magnetic field realizing a polarization of the nuclei of the protons of the specimen or subject for which an image is desired. Within magnet  10  there is a cavity or space in which the specimen or human to be examined is placed.  
         [0032]     The apparatus also includes a gradient system for producing spatial linear field gradients. These gradient fields are generally established by a set of three orthogonal direct current coils  11 ,  12  and  13 , which generate the three principal gradients G y , G x , and G z . These coils are driven by gradient generator  14 , which in turn is controlled by a controller  16  which communicates with the host computer  20 . Typical gradients used in MRI image processing are the well known slice select, readout, and phase encoding gradients.  
         [0033]     Typical MRI systems also generally include a radio frequency (RF) coil  17  which generates a radio frequency field in the specimen being analyzed and senses a free induction decay or spin echo signal which is generated after termination of the radio frequency pulse. RF pulse unit  18  excites RF coil  17 .  
         [0034]     The signal processor  19  receives the small microvoltage level spin echo signals which are reconstructed by computer  20  to form an image. The image is digitized and stored in the memory section of computer  20  for later display on display unit  21 .  
         [0035]     Protons with their magnetic axis aligned in the transverse plane are termed “saturated,” as will be discussed below. If the imaging step is conducted with a magnetic axis of a proton already in the transverse plane, then due to the subsequent gradient crushing or dephasing, the magnetic axes will be aligned orthogonally to the slice plane (along the Z axis) which produces no signal, thereby resulting in a perceptible void in the image.  
         [0036]     Referring to  FIGS. 2A and 2B , the present invention uses a SINC function to modulate the RF pulses  200  from generator  18  in a manner which yields relatively uniform rectangular magnetization profile tags  202 . As will be set forth below, the present invention allows for control of both the width of the tags and the separation of the tags.  
         [0037]     A first order approximation of the presaturated magnetization profile produced by an RF pulse train in the presence of a constant gradient G x  is the Fourier transform of the RF pulse train waveform. Assume the width of individual RF pulses is infinitesimal, the SINC modulated RF pulse train can be written as:  
               rf   ⁡     (   t   )       =     sin   ⁢           ⁢     c   ⁡     (       π   ⁢           ⁢   t       Δ   ⁢           ⁢     t   1         )       ×     comb   ⁡     (     t   ,     Δ   ⁢           ⁢     t   2         )                 [   1   ]             
 
 where  
       sin   ⁢           ⁢     c   ⁡     (       π   ⁢           ⁢   t       Δ   ⁢           ⁢     t   1         )           
 
 is the SINC function with the first zero crossing at Δt 1  and comb(t, Δt 2 ) is the comb function of spacing Δt 2 . Thus the spatial modulation of magnetization can be approximated as:  
                 M   sat     ⁡     (   x   )       =       rect   ⁡     (       γG   x     ⁢   Δ   ⁢           ⁢     t   1     ⁢   x     )       ⊕     comb   ⁡     (     x   ,     1       γG   x     ⁢   Δ   ⁢           ⁢     t   2           )                 [   2   ]             
 
 which provides rectangular tags of width 1/(γG x Δt 1 ) and separation 1/(γG x Δt 2 ). As is evident from this equation, both the tag width and separation are influenced by the gradient G x . Thus, in addition to influencing the tag line properties with the characteristics of the SINC function and comb function, the magnitude of the gradient can also be used to define these properties. 
 
         [0038]      FIGS. 2A, 2B ,  3 A,  3 B and  4 A,  4 B illustrate examples of SINC modulated RF pulse trains and corresponding saturated magnetization profiles of varying tag width to separation ratios, i.e., 1:2, 1:4 and 1:8, respectively. The ratio of tag width to tag separation is controlled by controlling the ratio of Δt 2  to Δt 1 . This approach can also be extended to 2D and 3D tagging without causing signal loss in untagged areas.  
         [0039]     While expressed above as an ideal comb function, In practical implementations, the width of individual RF pulses is finite. This finite width may cause undesirable shading across the image in the presence of the constant gradient G x . Assuming the duration of individual rectangular RF pulses is Δt 3 , the SINC modulated RF pulse train can be expressed as:  
               rf   ⁡     (   t   )       =       (     sin   ⁢           ⁢     c   ⁡     (       π   ⁢           ⁢   t       Δ   ⁢           ⁢     t   1         )       ×     comb   ⁡     (     t   ,     Δ   ⁢           ⁢     t   2         )         )     ⊕     rect   ⁡     (     t     Δ   ⁢           ⁢     t   3         )                 [   3   ]             
 
 and the saturated magnetization profile then approximated by:  
                 M   sat     ⁡     (   x   )       =       (       rect   ⁡     (       γG   x     ⁢   Δ   ⁢           ⁢     t   1     ⁢   x     )       ⊕     comb   ⁡     (     x   ,     1       γG   x     ⁢   Δ   ⁢           ⁢     t   2           )         )     ×   sin   ⁢           ⁢       c   ⁡     (     π   ⁢           ⁢   x   ⁢           ⁢   γ   ⁢           ⁢     G   x     ⁢   Δ   ⁢           ⁢     t   3       )       .               [   4   ]             
 
         [0040]     Thus the shading is characterized by a SINC modulation, sinc(πχγG x Δt 3 ). This shading can be visible and detrimental when tagging a large field of view with very small tag separation. Reducing the duration of individual RF pulses, Δt 3 , can overcome this problem. However, doing so may impose high RF peak power requirements in some applications, such as human cardiac studies.  
         [0041]     To address the problem of shading in a manner which does not put undo limitations on the peak power requirements of the RF generator, it has been discovered that the constant gradient G x  can be replaced by gradient segments between RF pulses without affecting the saturation profile as long as the gradient segment integral (i.e., the phase accumulation) between adjacent RF pulses is preserved.  FIGS. 6A and 6B  illustrate, by way of timing diagrams, a first example where gradient segments are active only between the RF pulses of the SINC modulated RF waveform.  
         [0042]      FIG. 5  is a simplified flowchart describing the generation of the waveforms depicted in the timing diagrams of  FIGS. 6A and 6B . As described above, the present invention provides for the selection of both tag width and separation by controlling parameters of the gradient, the SINC modulation function and the RF pulses of the comb function being modulated. In step  505 , the desired tag width is set by adjusting the lobe width of the SINC function (Δt 1 ). As is evident in equation 2, set forth above, both tag width and separation are also simultaneously affected by the gradient, G x . The tag separation can be adjusted by altering gradient and the time between RF pulses in the comb function (Δt 2 ), as illustrated in step  510 . While the gradient alters both the tag width and separation, it effects each of these parameters simultaneously. The tag width and separation are altered independently of each other by varying the SINC function and comb function as described above.  
         [0043]     Prior to activation of the first RF pulse, the gradient, G x , is turned off (step  515 ) and the RF pulse is then applied (step  520 ). At the end of the RF pulse, the gradient Gx is again turned on (step  525 ). This process is completed for each pulse in the SINC modulated RF pulse train. If there are additional pulses in the pulse train to be applied (step  530 ), the process idles with the gradient active until the start time of the next RF pulse (step  540 ). If there are no more pulses in the pulse train, then the process advances to conventional MRI image acquisition (step  535 ). The gradient, G x , is then deactivated (step  545 ) and the next RF pulse, whose amplitude and duration are determined by the SINC and comb function parameters, is applied (step  550 ). At the end of the pulse, the process returns, in a loop, to step  525 .  
         [0044]      FIGS. 7A and 7B  are timing diagrams which illustrate an equivalent example to that depicted in  FIG. 6 . As compared to the example in  FIG. 6A , in the case of  FIG. 7A , the amplitude of the individual RF pulses is reduced, but the duty cycle of such pulses have been extended to preserve the power of each pulse and maintain an equivalent resulting spatial magnetization tag. This alteration permits the use of longer pulse duration to reduce peak power requirement, and allows equalized RF pulse amplitudes to minimize the total pulse train duration, reducing delays to the initialization of imaging after the QRS complex trigger. This modification of the RF pulse could occur, for example, in steps  505  and/or  510  of the method depicted in  FIG. 5 .  
         [0045]      FIGS. 8A and 8B  illustrate a further equivalent timing relationship to that illustrated in  FIG. 6 . In  FIG. 8A , the amplitude of each of the RF pulse segments is now constant rather than varying in accordance with the SINC modulation. However, the duty cycle of the RF pulses, which was previously constant, is now varied to provide the SINC modulation of the RF signal. The gradient pulses illustrated in  FIG. 8B  are now pulse width modulated as well to correspond to the varying duty cycle of the RF signal. Note that SINC modulated RF pulse trains can be further iteratively optimized by simulating the Bloch equation numerically.  
         [0046]     The present methods are generally incorporated into conventional MRI apparatus by way of programming the host computer  20  and/or controller  16  in order to generate the desired waveforms from the RF pulse generator  18  and gradient generator  14 . The software is generally written in any number of conventional programming languages and can be stored and transported on conventional computer readable media, such as magnetic storage disks (floppy diskettes, hard disks and the like), optical disks (CD-ROMS) and the like.  
       Experimental Results  
       [0047]     The present invention has been implemented using a Bruker AVANCE 400WB spectrometer (Bruker NMR, Inc., Billerica, Mass.) with an 89 mm vertical bore magnet of 9.4 T (Oxford Instruments Ltd., UK) using a 30-mm-i.d. quadrature RF probe and a shielded gradient system up to 100 G/cm. The SINC modulated RF pulse trains in conjunction with a constant gradient illustrated in  FIGS. 2A, 2B ,  3 A,  3 B, and  4 A,  4 B, were implemented for both a phantom and an in vivo animal study. The RF pulse trains were 6 ms long and modulated by a 3-lobe SINC function. The duration of individual RF pulses was 20 μs, 40 μs and 80 μs using waveforms illustrated in  FIGS. 2A, 3A  and  4 A, respectively.  
         [0048]     The DANTE tagging method, consisting of an RF pulse train of uniform amplitude and constant gradient, was also implemented for comparison. To achieve a tag width to separation ratio equivalent to that in  FIG. 2B  for comparison (1:4), the number of RF pulses in DANTE pulse train was chosen to be four (4). In all experiments, the RF transmitter attenuation was calibrated for the tagging RF pulse train to produce a 90° flip angle. The tagging sequence was followed by a 1.5 ms crusher gradient.  
         [0049]     For the phantom study, a 20 mm tube of water doped with copper sulfate was used. Images were obtained by RF pulse train tagging followed by a 2D gradient echo (GE) encoding with the following sequence parameters: FOV=22 mm, acquisition matrix=256×256, slice thickness=2 mm, TR/TE=300/3 ms, flip angle=30°, number of averages=1. The RF pulse trains in  FIGS. 2A, 3A , and  4 A were employed to demonstrate various spatial modulation of magnetization.  
         [0050]     The in vivo mouse heart study was performed using a normal wild-type adult mouse (C57BL/6, 25 g). During the imaging experiment the mouse was anesthetized with isoflurane gas (1.5 vol. % at 2 L/min air flow) via a nose cone. The ECG was recorded from the front limbs using subcutaneous silver electrodes. During imaging heart rate was approximately 500 beats per minute. Imaging was performed with an ECG-gated fast 2D GE sequence using TR/TE=110/1.8 ms, FOV=26 mm, acquisition matrix=168×168, slice thickness=1.5 mm, flip angle=20°, number of averages=8. The total acquisition time was approximately 2.5 min. The tagging waveform in  FIG. 2B  was employed. There was a  3  ms delay between the end of the tagging sequence and the start of the GE imaging sequence, i.e., slice selective RF excitation. Delay between the QRS complex trigger and the start of the GE sequence was approximately 10 ms, therefore, images were acquired during ventricular systole.  
         [0051]      FIGS. 9A, 9B  and  9 C are images which illustrate the results from the phantom study employing the waveforms of  FIGS. 2A, 3A  and  4 A, respectively. As predicted, tags with sharp edges were obtained with SINC modulated RF pulse trains.  FIGS. 9A, 9B  and  9 C demonstrate the flexibility to alter the ratio of tag width to tag separation, which cannot be easily provided by the DANTE and SPAMM techniques.  FIG. 10  is an image acquired using the prior art DANTE method during the phantom study. The improvement of tag contrast over the conventional DANTE tagging technique is apparent.  FIGS. 11 and 12  illustrate ID tagging of a wild-type mouse heart by the SINC modulated RF pulse train of the present invention and the conventional DANTE technique, respectively. The nominal tag width and separation were 0.16 mm and 0.65 mm, respectively. The sharper-edged tagging by the SINC modulated RF train is apparent in the myocardium and the entire field of view of  FIG. 9A .  
         [0052]      FIG. 9D  illustrates an extension of the present invention into two dimensional tag lines. The tag lines are added in a direction orthogonal to the initial tag lines, thereby forming a grid, by following the same procedure described above for G x  with either the G y  or G z  gradients.  
         [0053]     Although the present invention has been described in connection with specific exemplary embodiments, it should be understood that various changes, substitutions and alterations can be made to the disclosed embodiments without departing from the spirit and scope of the invention as set forth in the appended claims.