The field of the invention is magnetic resonance imaging systems and methods. More particularly, the invention relates to a system and method for magnetic resonance imaging in which two or more magnetization preparation radio frequency pulses are utilized, such as double inversion recovery magnetic resonance imaging.
When a substance such as human tissue is subjected to a uniform magnetic field, B0, applied along, for example, the z-axis of a Cartesian coordinate system, the individual magnetic moments of the spins in the tissue attempt to align with this magnetic field, B0, but precess about the field in random order at their characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic field, B1, that is applied in the x-y plane and that is near the Larmor frequency of the spins, the net aligned moment, Mz, may be rotated, or “tipped,” into the x-y plane to produce a net transverse magnetic moment, Mxy. A magnetic resonance signal is emitted by the excited spins after the excitation field, B1, is terminated. This magnetic resonance signal may be received and processed to form an image or to produce a spectrum.
The magnetic resonance signals acquired with an MRI system are signal samples in Fourier space, or what is often referred to in the art as “k-space.” Typically, a region to be imaged is scanned by a sequence of measurement cycles in which magnetic field gradients vary according to the particular localization method being used. Each measurement cycle, or pulse sequence, typically samples a portion of k-space along a sampling trajectory characteristic of that pulse sequence. This is accomplished by employing magnetic field gradients (Gx, Gy, and Gz) that have the same direction as the polarizing magnetic field, B0, but which have a gradient along the respective x, y, and z axes. By controlling the strength of these gradients during each measurement cycle, the spatial distribution of spin excitation can be controlled and the location of the resulting magnetic resonance signals can be identified. The acquisition of the magnetic resonance signal samples is referred to as sampling k-space, and a scan is completed when enough measurement cycles are performed to adequately sample k-space. The resulting set of received magnetic resonance signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
Inversion recovery (“IR”) techniques are commonly used in MRI to either increase signal contrast between tissues, or otherwise null the signals originating from particular tissue types. IR acquisitions use the application of an inversion recovery radio frequency (“RF”) pulse before the application of an excitation RF pulse and subsequent data acquisition. The duration of time between the IR RF pulse and the excitation RF pulse is referred to as the inversion time (“TI”). The contrast of IR techniques is modulated by the T1 (spin-lattice) relaxation time. The TI is selected before the MRI scan begins and determines the image contrast. Signal nulling with IR techniques is performed by selecting the TI that corresponds to the time point with roughly zero longitudinal magnetization.
As an example, IR techniques are commonly performed in brain imaging applications in order to significantly reduce, or otherwise null, the magnetic resonance signals attributable to cerebrospinal fluid (“CSF”). A double (also called dual) IR (“DIR”) technique applies two inversion recovery RF pulses in succession, such that two tissues are simultaneously nulled at a time TI2 following the second inversion recovery pulse. In the brain, DIR techniques are usually set to null signals from both white matter and CSF, producing images of gray matter. The acquisition time of whole-brain DIR images is typically on the order of 10-15 minutes, which is long in a clinical setting. Furthermore, the DIR technique is usually performed in addition to a separate single IR acquisition, which provides T1-weighted images of the brain. As such, the acquisition of both a DIR and an IR image typically requires 14-23 minutes. DIR image acquisitions, however, provide useful information, and are particularly suitable for assessing certain diseases including multiple sclerosis (“MS”), epilepsy, and plaque build-up in carotid arteries. As a result, DIR acquisitions are often included in MRI brain scans for MS patients.
It would therefore be desirable to provide systems and methods for MRI that are capable of acquiring multiple different images representative of different recovery states of magnetization in a shorter amount of scan time that can currently be achieved.
The present invention overcomes the aforementioned drawbacks by providing systems and methods for magnetic resonance imaging (“MRI”) in which pulse sequences that utilize multiple magnetization preparation radio frequency (“RF”) pulses are modified so as to acquire multiple sets of data attributable to different images of the subject. In doing so, multiple different scans of the subject can be condensed into one, single scan that requires less time to perform than the multiple different scans combined.
It is an aspect of the invention to provide a method for producing at least two images of a subject with an MRI system. The MRI system is directed to perform a pulse sequence that, during each repetition time (“TR”) interval, includes the following. A first magnetization preparation RF pulse is applied to manipulate magnetization attributable to a spin species, and a second magnetization preparation RF pulse is applied to further manipulate the magnetization attributable to the spin species. First data is acquired after applying the first magnetization preparation RF pulse and before applying the second magnetization preparation RF pulse, and second data is acquired after applying the second magnetization preparation RF pulse. The method also includes reconstructing at least one image of the subject from the first data and reconstructing at least one image of the subject from the second data.
It is another aspect of the invention to provide a method for producing a plurality of images of a subject with an MRI system. The MRI system is directed to apply a plurality of magnetization preparation RF pulses, each magnetization preparation RF pulse rotating magnetization through a flip angle. After each of the plurality of magnetization preparation RF pulses is applied, k-space data is acquired by applying at least one excitation RF pulse after each of the plurality of magnetization preparation RF pulses and by sampling magnetic resonance signals formed in response to those excitation RF pulses. A plurality of images of the subject are then reconstructed from the acquired k-space data.
The foregoing and other aspects and advantages of the invention will appear from the following description. In the description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full scope of the invention, however, and reference is made therefore to the claims and herein for interpreting the scope of the invention.
Described here are systems and methods for magnetic resonance imaging (“MRI”) using two or more magnetization preparation radio frequency (“RF”) pulses, such as two or more inversion recovery (“IR”) pulses. The systems and methods of the present invention are thus capable of acquiring multiple sets of IR data in a single scan. As one example, the systems and methods are capable of acquire both double IR (“DIR”) data and single IR data in the same scan. The total scan time achievable with the present invention is significantly reduced as compared to acquiring separate DIR and IR scans required with conventional techniques. Using existing techniques, it takes a total of 14-22 minutes to perform both a single IR scan and a DIR scan. Using the systems and methods of the present invention, both single IR data and DIR data can be obtained in a single scan that takes about 10-15 minutes to complete.
Current DIR techniques apply an IR radio frequency (“RF”) pulse, wait a specified amount of time, apply a second IR RF pulse, and then begin sampling magnetic resonance signals at a specified time after the second IR RF pulse is applied. The first and second specified time are selected before the MRI scan is performed such that data is acquired when two tissues with different T1 values are nulled, or otherwise significantly reduced. As an example, the specified times can be selected such that white matter and cerebrospinal fluid are nulled when data is acquired. The magnetization is usually allowed to regrow before the imaging process is repeated again. DIR is attractive because it is better able to depict intracortical lesions, such as those associated with multiple sclerosis, due to the high gray matter-white matter contrast achievable with DIR acquisitions.
The systems and methods of the present invention are fundamentally different than existing techniques because data is acquired after each IR RF pulse. In conventional DIR, the interval between the first and second IR RF pulses is merely a waiting period during which no data acquisition occurs. Using the systems and methods of the present invention, however, data is additionally acquired after the first IR RF pulse.
By way of example, data can be acquired after the first IR RF pulse using a technique such as those described in co-pending U.S. patent application Ser. No. 13/829,260, which is herein incorporated by reference in its entirety, to produce a spectrum of hundreds of images (e.g. 200-500) images each having a different T1 contrast, one of which has contrast similar to a typical T1-weighted exam. An example of the data acquisition technique described in U.S. patent application Ser. No. 13/829,260 is also provided in more detail below. The systems and methods of the present invention can thus employ a unique data acquisition scheme that allows for the retrospective identification of an image having an optimal tissue contrast. In these implementations, data acquisition is performed using a radial acquisition scheme such as, preferably, a vastly undersampled isotropic projection reconstruction (“VIPR”) scheme.
Alternatively, data can be acquired following each IR pulse using conventional Cartesian k-space sampling or other sampling trajectories, including radials, spirals, twirls, sinusoids, and so on. As will be described below, the sampling schemes may be different following each magnetization preparation RF pulse.
One advantage of the present invention is that multiple images with different contrast characteristics can be obtained in a single scan. For example, current brain imaging protocols require two separate scans to obtain a T1-weighted image and an image in which signals from cerebrospinal fluid are suppressed. With the systems and methods of the present invention, however, both a T1-weighted image and an image in which signals from cerebrospinal fluid are suppressed can be obtained in a single scan, thereby providing a significant time saving for the subject being imaged and increasing subject throughput.
Another advantage of the present invention is that a single T1-weighted image can be obtained in addition to a plurality of images each having different image contrasts, a so-called “spectrum” of images, in which the different image contrasts are defined by the recovery of magnetization following a double IR technique. In this manner, a subject and disease specific optimized image contrast can be achieved because, rather than trying to predict what imaging parameters will produce the best contrast for a particular subject, multiple images across a wide range of contrast settings can be acquired. This can be useful for many different applications, including pediatric imaging.
The method of the present invention can also be advantageously used for T1 quantification. The methods described here can sample data along a recovery curve that can be modeled and fitted to match the data. This allows the T1 value, a fundamental parameter describing tissue, to be accurately quantified. Being able to obtain an accurate measurement of T1 allows for new biomarkers for classifying and monitoring progression of illnesses and normal development with new tools not currently available on commercial scanners.
In general, a pulse sequence that implements the present invention includes a plurality of magnetization preparation RF pulses, each followed by a data acquisition, and all occurring within a single repetition time (“TR”) interval. For instance, a pulse sequence that implements the present invention may include a first magnetization preparation RF pulse followed by a first data acquisition, and a second magnetization preparation RF pulse followed by a second data acquisition. By way of example, the magnetization RF pulses can each be IR RF pulses. Unlike previous multiple IR methods, the method of the present invention acquires data as longitudinal magnetization is recovering following the application of each IR RF pulse. In other implementations, the magnetization preparation RF pulses may be saturation RF pulses or magnetization transfer RF pulses. The magnetization preparation RF pulses may be spatially selective, spectrally selective, or non-selective.
An example of a pulse sequence that implements the present invention is illustrated in
Unlike conventional DIR pulses sequences, after the application of the first magnetization preparation RF pulse 102 and before the application of the second magnetization preparation RF pulse 104, data is acquired from the subject being imaged, as indicated by data acquisition block 106. Data acquisition commences with the performance of a data acquisition sequence 108. By way of example, data acquisition sequence 108 includes multiple excitation RF pulses. The excitation RF pulses can each have the same flip angle, α, as illustrated by way of example in
After the application of the second magnetization preparation RF pulse 104, data is acquired from the subject being imaged, as indicated by data acquisition block 110. Data acquisition commences with the performance of a data acquisition sequence 112. By way of example, data acquisition sequence 112 includes multiple excitation RF pulses. The excitation RF pulses can each have the same flip angle, α, as illustrated by way of example in
In some configurations, one of the data acquisition sequences 108, 112 may include only a single excitation RF pulse. If a single excitation RF pulse with a flip angle at or near ninety degrees is employed in the first data acquisition sequence 108, the second inversion time, TI2, is preferably lengthened relative to a sequence that makes use of multiple excitation RF pulses with small flip angles. The reasoning for lengthening the second inversion time, TI2, when using only a single excitation RF pulse with a flip angle at or near ninety degrees in the first data acquisition sequence 108 is to allow additional time for longitudinal magnetization to regrow before the commencement of the second data acquisition sequence 112. If a single excitation RF pulse having a flip angle at or near ninety degrees is employed in the second data acquisition sequence 112, and not the first data acquisition sequence 108, then this lengthening of the second inversion time, TI2, may not be necessary.
In general, the data acquisition sequences 108, 112 can sample k-space in any one of a number of different manners. For example, when acquiring data throughout the magnetization recovery using the data acquisition techniques described in co-pending U.S. patent application Ser. No. 13/829,260, k-space is sampled along three-dimensional radial trajectories. Other k-space sampling patterns can also be used, however. As another example, Cartesian sampling can be implemented. For application where Cartesian sampling is implemented, typically data associated with only one image will be acquired during the data acquisition sequence that utilizes Cartesian sampling. Thus, when Cartesian sampling is implemented it is preferable that an inversion time be selected prior to performing the scan, as is done with conventional IR methods. By way of example, the first data acquisition sequence 108 could implement a Cartesian sampling pattern to obtain a single, T1-weighted image while the second data acquisition sequence 112 could implement the three-dimensional radial sampling technique described in co-pending U.S. patent application Ser. No. 13/829,260 to obtain a spectrum of images each having different image contrasts. As another example, the first data acquisition sequence 108 could implement the three-dimensional radial sampling technique described in co-pending U.S. patent application Ser. No. 13/829,260 to obtain a spectrum of images each having different image contrasts and the second data acquisition sequence 112 could implement a Cartesian sampling pattern to obtain a single, DIR image.
The pulse sequence is repeated until the desired amount of data has been acquired. For each TR interval, data associated with one or more images are obtained during the first data acquisition sequence 108 and data associated with one or more images are obtained during the second data acquisition sequence 112.
By way of example, the data acquired in the first data acquisition sequence 108 can be reconstructed to form a single, T1-weighted image of the subject. Alternatively, if the data acquisition techniques described in co-pending U.S. patent application Ser. No. 13/829,260 are implemented during the first data acquisition sequence 108, then the data acquired in the first data acquisition sequence 108 can be reconstructed to form a plurality of different images, each having a different image contrast determined by the recovery of longitudinal magnetization back to equilibrium following the application of the first magnetization preparation RF pulse 102.
By way of example, the data acquired in the second data acquisition sequence 112 can be reconstructed to form a single, DIR image of the subject in which signals from two different tissues are significantly reduced, or otherwise nulled. Alternatively, if the data acquisition techniques described in co-pending U.S. patent application Ser. No. 13/829,260 are implemented during the second data acquisition sequence 112, then the data acquired in the second data acquisition sequence 112 can be reconstructed to form a plurality of different images, each having a different image contrast determined by the recovery of longitudinal magnetization back to equilibrium following the application of the second magnetization preparation RF pulse 104.
It will be appreciated that the example pulse sequence illustrated in
Referring particularly now to
The pulse sequence server 310 functions in response to instructions downloaded from the workstation 302 to operate a gradient system 318 and an RF system 320. Gradient waveforms used to perform the prescribed scan are produced and applied to the gradient system 318, which excites gradient coils in an assembly 322 to produce the magnetic field gradients Gx, Gy, and Gz used for position encoding magnetic resonance signals. The gradient coil assembly 322 forms part of a magnet assembly 324 that includes a polarizing magnet 326 and a whole-body RF coil 328.
RF excitation waveforms are applied to the RF coil 328, or a separate local coil (not shown in
The RF system 320 also includes one or more RF receiver channels. Each RF receiver channel includes an RF preamplifier that amplifies the magnetic resonance signal received by the coil 328 to which it is connected, and a detector that detects and digitizes the I and Q quadrature components of the received magnetic resonance signal. The magnitude of the received magnetic resonance signal may thus be determined at any sampled point by the square root of the sum of the squares of the I and Q components:
M=√{square root over (I2+Q2)} (1);
and the phase of the received magnetic resonance signal may also be determined:
The pulse sequence server 310 also optionally receives patient data from a physiological acquisition controller 330. The controller 330 receives signals from a number of different sensors connected to the patient, such as electrocardiograph (“ECG”) signals from electrodes, or respiratory signals from a bellows or other respiratory monitoring device. Such signals are typically used by the pulse sequence server 310 to synchronize, or “gate,” the performance of the scan with the subject's heart beat or respiration. For instance, the application of the magnetization preparation RF pulses, such as the first magnetization preparation RF pulse applied in each TR interval, can be timed relative to an ECG signal so that data acquisitions will occur during specified cardiac phases.
The pulse sequence server 310 also connects to a scan room interface circuit 332 that receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 332 that a patient positioning system 334 receives commands to move the patient to desired positions during the scan.
The digitized magnetic resonance signal samples produced by the RF system 320 are received by the data acquisition server 312. The data acquisition server 312 operates in response to instructions downloaded from the workstation 302 to receive the real-time magnetic resonance data and provide buffer storage, such that no data is lost by data overrun. In some scans, the data acquisition server 312 does little more than pass the acquired magnetic resonance data to the data processor server 314. However, in scans that require information derived from acquired magnetic resonance data to control the further performance of the scan, the data acquisition server 312 is programmed to produce such information and convey it to the pulse sequence server 310. For example, during prescans, magnetic resonance data is acquired and used to calibrate the pulse sequence performed by the pulse sequence server 310. Also, navigator signals may be acquired during a scan and used to adjust the operating parameters of the RF system 320 or the gradient system 318, or to control the view order in which k-space is sampled. In all these examples, the data acquisition server 312 acquires magnetic resonance data and processes it in real-time to produce information that is used to control the scan.
The data processing server 314 receives magnetic resonance data from the data acquisition server 312 and processes it in accordance with instructions downloaded from the workstation 302. Such processing may include, for example: Fourier transformation of raw k-space data to produce two or three-dimensional images; the application of filters to a reconstructed image; the performance of a backprojection image reconstruction of acquired magnetic resonance data; the generation of functional magnetic resonance images; and the calculation of motion or flow images.
Images reconstructed by the data processing server 314 are conveyed back to the workstation 302 where they are stored. Real-time images are stored in a data base memory cache (not shown in
An example of the data acquisition technique described in co-pending U.S. patent application Ser. No. 13/829,260 is now provided in more detail. This technique provides two-dimensional or three-dimensional magnetic resonance imaging that allows improved flexibility in creating segments and in spacing TRs. In this technique, any number of consecutive projections can be combined to improve image quality. The number of consecutive projections to be combined can be chosen after the data has been acquired, thereby providing for the retrospective combination data. Furthermore, the set of a desired number of projections can be selectively adjusted by a single TR, thereby allowing for the effective inversion time to be adjusted by a given, even single, TR. Multiple sequential TRs can thus have increasing inversion times. It is noted that the transverse magnetization prior to each excitation can be spoiled using, for example, variable RF phase or gradient pulses; although, magnetization spoiling is not essential. The choice of whether to use spoiling, however, will affect image contrast.
By way of example, a clinician may choose to use eight consecutively acquired projections to reconstruct an image. If so, these projections are combined together as a data segment, S1. Because this segment includes projections P1 through P8, its average inversion time, TI, is determined by 1+2+3+4+5+6+7+8, or 36TR, divided by 8. If the clinician then selected projections P2 through P9 as the next data segment, S2, this next data segment will have a TI determined by 2+3+4+5+6+7+8+9, or 44TR, divided by 8. Notably the TI spacing of 44TR/8−36TR/8 is equal to 8TR/8, which is the TR between each data frame.
An example of this data acquisition scheme is illustrated in
Advantageously, the view ordering used during data acquisition is selected such that subsets of consecutively acquired projections substantially uniformly sample three-dimensional k-space during a time window 410 of any suitable duration. As a result, any subset of consecutively acquired projections can be selected as a data segment while still maintaining a uniform sampling of k-space. This flexibility allows for finer temporal resolution than achievable with previous projection-based IR imaging methods. For instance, temporal resolutions on the order of five milliseconds are obtainable with the present invention.
As noted above, images are reconstructed from those projections contained within a sliding temporal window that can be selected by the user after data acquisition has occurred. This process is illustrated in
Multiple inversion recovery sequences can also be performed, with the sampling between repetitions differing, as illustrated in
Also, as illustrated in
Referring now to
During the MRI scan of the subject, a magnetization preparation RF pulse, such as an IR RF pulse, associated with a given data acquisition sequence is applied, as indicated at step 706. While magnetization is recovering back to equilibrium following the application of the magnetization preparation RF pulse, k-space data is acquired from the subject in accordance with the sampling pattern and view ordering selected for the given data acquisition sequence, as indicated at step 708.
By way of example, k-space data is acquired as described above; namely, by sampling k-space along one or more radial projections during a the data acquisition sequence. A series of consecutive radial projections are sampled in this manner. Preferably, the radial projections are selected such that a selected number of consecutive radial projections substantially uniformly sample two-dimensional or three-dimensional k-space.
After all of the desired k-space data for a given data acquisition sequence have been acquired, a temporal window is selected, as indicated at step 710. This temporal window defines how many consecutively acquired projections should be included from each k-space data set when reconstructing an image. The data are grouped into data segments using the selected temporal window, as indicated at step 712. For instance, a selected number of consecutive projections are grouped into a data segment. The number of projections included in each data segment is determined by the selected temporal window, which can be defined based on the desired temporal resolution for the image reconstruction. As described above, preferably the data are grouped into data segments that are temporally shifted by only a single TR period. Thus, each data segment can be defined by an initial radial projection acquired at an initial time point and by the desired number of radial projections to be included in each data segment. After the k-space data have been grouped into data segments, images are reconstructed from the data segments, as indicated at step 714. A time series of images that depicts the subject as magnetization recovers following the application of the magnetization preparation RF pulses is thus formed. Each image in this time series has a unique image contrast that is determined by the recovery state of the magnetization at which the data in the associated data segment were acquired.
Because this time series of images can be designed to span almost the entirety of the magnetization recovery curve, it is not necessary to prospectively select an inversion time to obtain an image in which signals attributable to a specific tissue are nulled or when contrast between two different tissues is maximized. Rather, because of the high temporal resolution of this data acquisition scheme, such an image can be retrospectively identified in the time series, or spectrum, of images. This feature is highly advantageous for pediatric imaging where the precise inversion time for a child patient may not be readily known. Because children are rapidly developing, there is no consensus on what inversion time to use to achieve a particular tissue nulling or contrast. The lack of consensus is also due to the fact that children's brains do not develop homogeneously, meaning that no single inversion time may be ideal for a given child's entire head. This data acquisition scheme thus provides a solution to this problem by allowing the retrospective identification of an image with the optimal tissue contrast.
The present invention has been described in terms of one or more preferred embodiments, and it should be appreciated that many equivalents, alternatives, variations, and modifications, aside from those expressly stated, are possible and within the scope of the invention.