US20170160366A1 - System and method for magnetic resonance i maging with prospective motion control - Google Patents
System and method for magnetic resonance i maging with prospective motion control Download PDFInfo
- Publication number
- US20170160366A1 US20170160366A1 US15/301,868 US201515301868A US2017160366A1 US 20170160366 A1 US20170160366 A1 US 20170160366A1 US 201515301868 A US201515301868 A US 201515301868A US 2017160366 A1 US2017160366 A1 US 2017160366A1
- Authority
- US
- United States
- Prior art keywords
- navigator
- imaging
- pulse sequence
- motion
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
- G01R33/5673—Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker
Definitions
- the present disclosure relates to systems and methods for magnetic resonance imaging (“MRI”). More particularly, the disclosure relates to systems and methods for tracking and controlling artifacts caused by motion during a MRI procedure.
- MRI magnetic resonance imaging
- MRI uses the nuclear magnetic resonance (“NMR”) phenomenon to produce images.
- NMR nuclear magnetic resonance
- a substance such as human tissue is subjected to a uniform magnetic field (“main magnetic field”), B 0
- the individual magnetic moments of the nuclei in the tissue attempt to align with this magnetic field, but precess about it in random order at their characteristic Larmor frequency, ⁇ .
- B 1 an excitation magnetic field
- B 1 an excitation magnetic field
- B 1 that is in the plane transverse to the main magnetic field
- B 0 and that is near the Larmor frequency, ⁇
- the net aligned magnetic moment of the nuclei may be rotated, or “tipped,” into the transverse plane to produce a net transverse magnetic moment.
- a signal is emitted by the excited nuclei, or “spins,” after the excitation magnetic field, B 1 , is terminated. The emitted signal may be received and processed to form an image.
- magnetic field gradients (G x , G y , and G z ) are employed.
- the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used.
- the resulting set of received MR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
- the measurement cycle used to acquire each MR signal is performed under the direction of a pulse sequence produced by a pulse sequencer.
- Clinically available MRI systems store a library of such pulse sequences that can be prescribed to meet the needs of many different clinical applications.
- Research MRI systems include a library of clinically-proven pulse sequences and they also enable the development of new pulse sequences.
- one popular category of pulse sequences are so-called gradient echo sequences.
- the spoiled gradient echo and three-dimensional spoiled gradient-recalled echo (SPGR) or fast low angle shoot (FLASH) pulse sequences are often used in neuroimaging applications.
- SPGR or FLASH sequence forms the basis of many 3D neuroimaging sequences, but acquisition times often stretch to several minutes. Lengthy acquisitions in neuroimaging applications can be particularly troublesome because it is imperative that the subject remain motionless during the duration of the sequence. That is, in neuroimaging applications, motion can be particularly damaging to the resulting images because of the complexity of the structures being imaged and studied in neuroimaging applications.
- Motion-correction systems in MRI can be grouped into two general methods: prospective and retrospective.
- Retrospective methods use information about the subject's motion to estimate what k-space data would have been measured if the subject had not moved during scanning.
- Prospective methods use motion-tracking data acquired during the scan to follow the subject with the gradient axes of the sequence, measuring the desired k-space data directly. Additionally, it is possible to combine the two methods so that retrospective processing corrects residual errors in the prospective system.
- a retrospective system can access all of the k-space data while performing reconstruction; a prospective system must necessarily rely only on previous measurements to estimate the current position of the patient. However, a prospective system avoids the need to estimate missing k-space data, allowing for direct reconstruction while avoiding possible sources of estimation error in the k-space data.
- two types of motion correction problems that arise in MRI can differentiate between two types of motion correction problems that arise in MRI: between-scan motion and within-scan motion.
- between-scan motion several retrospective motion correction methods are available that register either slice-by-slice or volume-by-volume to estimate the data that would have been acquired in each volume if the subject had not moved.
- Prospective motion correction can also be employed for this problem, such as the orbital navigator system that inserts 3-plane circular k-space navigators, or the PACE system that registers each completed EPI volume back to the first time-point and so requires no navigators.
- short TR protocols are often used in neuroimaging applications.
- Such short-TR protocols have very-little dead time and, thus, force navigators to be very short and provide limited k-space coverage. That is, as the TR is reduced, the effectiveness of the navigator is reduced because there is less information gathered by the navigator to use to form an estimate of the subject's head motion.
- the present disclosure overcomes the aforementioned drawbacks by providing a system and method for directing a magnetic resonance imaging (MRI) system to reconstruct an image of a subject, in which motion has been prospectively corrected.
- the method includes acquiring image data and navigator data from a subject by directing the MRI system to perform an imaging pulse sequence that includes a navigator portion.
- the navigator data is processed to determine whether motion occurred while the image data was acquired.
- the imaging pulse sequence is updated when processing the navigator data determines motion occurred.
- the imaging pulse sequence is updated by first acquiring volume navigator data by directing the MRI system to perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject.
- EPI echo-planar imaging
- the volume navigator data is then processed to generate motion estimates of motion in the VOI.
- the imaging pulse sequence is then updated based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates.
- the imaging pulse sequence, whether updated or not, is then repeated and navigator data processed to determine whether motion occurred and whether the imaging pulse sequence should be further updated until a desired amount of image data is acquired. An image of the subject is then reconstructed from the acquired image data.
- a magnetic resonance imaging (MRI) system includes a magnet system configured to generate a polarizing magnetic field about at least a portion of a subject arranged in the MRI system, a magnetic gradient system including a plurality of magnetic gradient coils configured to apply at least one magnetic gradient field to the polarizing magnetic field, and radio frequency (RF) system configured to apply an RF field to the subject and to receive magnetic resonance signals from the subject using a coil array.
- the system also includes a computer system programmed to control the gradient system and the RF system to perform an imaging pulse sequence that includes a navigator portion to acquire image data and navigator data from a subject.
- the computer system is further programmed to process the navigator data to determine whether the subject moved while the image data was acquired.
- the computer system Upon determining that the subject moved, the computer system is configured to control the gradient system and the RF system to perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject, process the volume navigator data to generate motion estimates of motion in the VOI, and update the imaging pulse sequence based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates.
- the computer system is further configured to reconstruct an image of the subject from the acquired image data.
- EPI echo-planar imaging
- FIG. 1 is a block diagram of an example of a magnetic resonance imaging (“MRI”) system.
- MRI magnetic resonance imaging
- FIG. 2 is an example of a three-dimensional spoiled gradient-recalled echo (“3D SPGR”) pulse sequence.
- 3D SPGR three-dimensional spoiled gradient-recalled echo
- FIG. 3 is an illustrative example of an integrated pulse sequence that includes both imaging sequences, each played out during an imaging repetition time (“TR”), and volumetric navigator sequences, each played out during a volumetric navigator TR;
- TR imaging repetition time
- FIG. 4 is a flowchart setting forth the steps of an example method for performing prospective motion correction using a pulse sequence such as the one described with respect to FIG. 3 ;
- FIG. 5 is a flowchart setting forth the steps of an example method for performing prospective motion correction based on an integrated pulse sequence that uses an ultra-fast navigator in an imaging pulse sequence to detect motion and trigger the performance of one or more vNav sequences.
- the MRI system 100 includes an operator workstation 102 , which will typically include a display 104 ; one or more input devices 106 , such as a keyboard and mouse; and a processor 108 .
- the processor 108 may include a commercially available programmable machine running a commercially available operating system.
- the operator workstation 102 provides the operator interface that enables scan prescriptions to be entered into the MRI system 100 .
- the operator workstation 102 may be coupled to four servers: a pulse sequence server 110 ; a data acquisition server 112 ; a data processing server 114 ; and a data store server 116 .
- the operator workstation 102 and each server 110 , 112 , 114 , and 116 are connected to communicate with each other.
- the servers 110 , 112 , 114 , and 116 may be connected via a communication system 140 , which may include any suitable network connection, whether wired, wireless, or a combination of both.
- the communication system 140 may include both proprietary or dedicated networks, as well as open networks, such as the internet.
- the pulse sequence server 110 functions in response to instructions downloaded from the operator workstation 102 to operate a gradient system 118 and a radiofrequency (“RF”) system 120 .
- Gradient waveforms necessary to perform the prescribed scan are produced and applied to the gradient system 118 , which excites gradient coils in an assembly 122 to produce the magnetic field gradients G x , G y , and G z used for position encoding magnetic resonance signals.
- the gradient coil assembly 122 forms part of a magnet assembly 124 that includes a polarizing magnet 126 and a whole-body RF coil 128 .
- RF waveforms are applied by the RF system 120 to the RF coil 128 , or a separate local coil (not shown in FIG. 1 ), in order to perform the prescribed magnetic resonance pulse sequence.
- Responsive magnetic resonance signals detected by the RF coil 128 , or a separate local coil (not shown in FIG. 1 ) are received by the RF system 120 , where they are amplified, demodulated, filtered, and digitized under direction of commands produced by the pulse sequence server 110 .
- the RF system 120 includes an RF transmitter for producing a wide variety of RF pulses used in MRI pulse sequences.
- the RF transmitter is responsive to the scan prescription and direction from the pulse sequence server 110 to produce RF pulses of the desired frequency, phase, and pulse amplitude waveform.
- the generated RF pulses may be applied to the whole-body RF coil 128 or to one or more local coils or coil arrays (not shown in FIG. 1 ).
- the RF system 120 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 128 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, therefore, be determined at any sampled point by the square root of the sum of the squares of the I and Q components:
- phase of the received magnetic resonance signal may also be determined according to the following relationship:
- the pulse sequence server 110 also optionally receives patient data from a physiological acquisition controller 130 .
- the physiological acquisition controller 130 may receive signals from a number of different sensors connected to the patient, such as electrocardiograph (“ECG”) signals from electrodes, or respiratory signals from a respiratory bellows or other respiratory monitoring device.
- ECG electrocardiograph
- Such signals are typically used by the pulse sequence server 110 to synchronize, or “gate,” the performance of the scan with the subject's heart beat or respiration.
- the pulse sequence server 110 also connects to a scan room interface circuit 132 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 132 that a patient positioning system 134 receives commands to move the patient to desired positions during the scan.
- the digitized magnetic resonance signal samples produced by the RF system 120 are received by the data acquisition server 112 .
- the data acquisition server 112 operates in response to instructions downloaded from the operator workstation 102 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 112 does little more than pass the acquired magnetic resonance data to the data processor server 114 . However, in scans that require information derived from acquired magnetic resonance data to control the further performance of the scan, the data acquisition server 112 is programmed to produce such information and convey it to the pulse sequence server 110 . For example, during prescans, magnetic resonance data is acquired and used to calibrate the pulse sequence performed by the pulse sequence server 110 .
- navigator signals may be acquired and used to adjust the operating parameters of the RF system 120 or the gradient system 118 , or to control the view order in which k-space is sampled.
- the data acquisition server 112 may also be employed to process magnetic resonance signals used to detect the arrival of a contrast agent in a magnetic resonance angiography (“MRA”) scan.
- MRA magnetic resonance angiography
- the data acquisition server 112 acquires magnetic resonance data and processes it in real-time to produce information that is used to control the scan.
- the data processing server 114 receives magnetic resonance data from the data acquisition server 112 and processes it in accordance with instructions downloaded from the operator workstation 102 .
- processing may, for example, include one or more of the following: reconstructing two-dimensional or three-dimensional images by performing a Fourier transformation of raw k-space data; performing other image reconstruction algorithms, such as iterative or backprojection reconstruction algorithms; applying filters to raw k-space data or to reconstructed images; generating functional magnetic resonance images; calculating motion or flow images; and so on.
- Images reconstructed by the data processing server 114 are conveyed back to the operator workstation 102 where they are stored.
- Real-time images are stored in a data base memory cache (not shown in FIG. 1 ), from which they may be output to operator display 112 or a display 136 that is located near the magnet assembly 124 for use by attending physicians.
- Batch mode images or selected real time images are stored in a host database on disc storage 138 .
- the data processing server 114 notifies the data store server 116 on the operator workstation 102 .
- the operator workstation 102 may be used by an operator to archive the images, produce films, or send the images via a network to other facilities.
- the MRI system 100 may also include one or more networked workstations 142 .
- a networked workstation 142 may include a display 144 ; one or more input devices 146 , such as a keyboard and mouse; and a processor 148 .
- the networked workstation 142 may be located within the same facility as the operator workstation 102 , or in a different facility, such as a different healthcare institution or clinic.
- the networked workstation 142 may gain remote access to the data processing server 114 or data store server 116 via the communication system 140 . Accordingly, multiple networked workstations 142 may have access to the data processing server 114 and the data store server 116 . In this manner, magnetic resonance data, reconstructed images, or other data may be exchanged between the data processing server 114 or the data store server 116 and the networked workstations 142 , such that the data or images may be remotely processed by a networked workstation 142 . This data may be exchanged in any suitable format, such as in accordance with the transmission control protocol (“TCP”), the internet protocol (“IP”), or other known or suitable protocols.
- TCP transmission control protocol
- IP internet protocol
- systems such as described above with respect to FIG. 1 have been used to perform neuroimaging acquisitions with short repetition times (“TRs”), which force navigators to be very short and, therefore, provide limited k-space coverage.
- TRs short repetition times
- These navigators are often referred to as “ultra-fast navigators.”
- the effectiveness of the navigators at providing an estimate of the subject's head motion is, therefore, also limited.
- vNays volume navigators
- EPI echo-planar imaging
- vNays volume navigators
- TI inflow times
- TR gaps TR gaps in which the entire vNav could be inserted.
- the present disclosure provides a way to insert a vNav into a 3D FLASH/SPGR sequence with only marginal impact on SNR or time. By doing so, the present disclosure provides a system and method for prospective motion correction of in-scan motion using a 3D FLASH/SPGR pulse sequence with vNays that does not substantially extend acquisition times.
- the pulse sequence includes an RF excitation pulse 202 that is played out in the presence of a slab-selective gradient 204 in order to produce transverse magnetization in a volume-of-interest.
- the slab-selective gradient 204 includes a rephasing lobe 206 that acts to rephase unwanted phase dispersions introduced by the slab-selective gradient 204 such that signal losses resultant from these phase dispersions are mitigated.
- a phase-encoding gradient 208 is applied to spatially encode a nuclear magnetic resonance echo signal 210 at a given echo time (TE) representative of a gradient-recalled echo along one direction in the volume-of-interest.
- a partition-encoding gradient 212 is applied to spatially encode the echo signal 210 along a second, orthogonal direction in the volume-of-interest.
- the phase-encoding gradient 208 may spatially encode the echo signal 210 along the y-direction
- the partition-encoding gradient 212 may spatially encode the echo signal 210 along the z-direction.
- a readout gradient 214 is also applied after a dephasing gradient lobe 216 to spatially encode the echo signal 210 along a third, orthogonal direction in the volume-of-interest.
- the readout gradient 214 may spatially encode the echo signal along the x-direction.
- the echo signal 210 is sampled during a data acquisition window.
- Spoiler gradients 218 , 220 , 222 may be played out along the partition-encoding, phase-encoding, and frequency-encoding directions to dephase any residual transverse magnetization in the volume-of-interest to prevent signal contamination from one repetition time (TR) period to the next.
- FIG. 3 provides a timing diagram of an integrated vNav and a SPGR/FLASH sequence 300 , showing common pulses and TRs, but varied gradients, for interleaved imaging and vNav sequences.
- a plurality of imaging sequences 302 such as described above with respect to FIG. 2 are illustrated.
- Integrated therewith are a plurality of vNav sequences 304 .
- the scan time of the imaging sequence 302 is selected such that an imaging TR 306 allows vNav sequences 304 , each having a vNav TR 308 , to be inserted at any desired point.
- the RF pulses used in the vNav sequences 304 and the imaging sequences 302 may be the same. To this end, the steady-state of the imaging sequences 302 is maintained, even when vNav sequences 304 are inserted in the integrated sequence 300 .
- the vNav sequence 304 may be a 3D-encoded echo-planar imaging (EPI) pulse sequence with a 32 3 matrix and may be acquired with 3 ⁇ 4 partial Fourier encoding in the partition direction.
- the vNav sequence 304 may have 25 three-dimensional excitation pulses.
- the vNav TR 308 is assumed to be 11 ms and the TR 306 for each imaging sequence 302 is assumed to match the vNav TR 308 .
- a train of 25 vNav TRs (i.e., one vNav) can be played instead of a TR 306 of an imaging sequence 302 without disturbing the steady state of the imaging sequence 302 .
- the TR 306 of the imaging sequence 302 is 11 ms or longer
- matching the vNav TR 308 to the imaging TR 306 can be readily achieved by adding TRs 310 to the vNav sequence 304 after its readouts and not simply attempting to fit the vNav TR 308 within dead time.
- vNav TRs 308 can be inserted at any desired point while maintaining the steady-state of the imaging sequence.
- Inserting the vNav sequences 304 as separate TRs 308 at any time allows a great deal of flexibility, but also increases overall scan time.
- the navigators were inserted in TR or TI gaps and, thus, did not increase scan time.
- 25 additional TRs are added to the total scan time with every vNav sequence while not gaining any additional imaging signal.
- the user may set how many imaging TRs will be played between each vNav sequence, and the overall scan time is updated to inform users of the trade-off between tracking accuracy and scan time.
- an additional time savings can be realized by selectively performing a limited number of vNav sequences 304 such that the vNav sequences 304 are performed only when a predetermined amount of motion is occurring.
- motion detection can be provided by inserting an ultra-fast navigator, such as a cloverleaf navigator, into each imaging sequence 302 . Examples of cloverleaf navigators are described in U.S. Pat. Nos. 6,771,068 and 6,958,605, both of which are herein incorporated by reference in their entirety.
- the ultrafast-navigator can be used as a high-sensitivity motion detector that can be relied on to determine when motion occurs without a significant scan time cost.
- one or more vNav sequences 304 can be performed to provide a more accurate measurement of the motion, as described above in detail. In this manner, the more time consuming vNav sequences 304 can be performed only when a predetermined amount of motion is detected by an ultra-fast navigator, thereby reducing overall scan time by eliminating the performance of unnecessary vNav sequences 304 . Integrating an ultra-fast navigator into the imaging sequences 302 also has the added benefit of enabling resonance frequency correction.
- the prospective motion-corrected process 400 begins with user selection of the imaging pulse sequences and parameters for the sequence at process block 402 .
- the imaging pulse sequence may be a FLASH or SPGR or other name for such gradient-echo sequences.
- the user is prompted to select how many imaging TRs will be played between each vNav TR.
- the overall scan time is updated to inform users of the trade-off between tracking accuracy and scan time.
- the user may adjust the number and spacing of the vNav sequences to reach a desired overall scan time until the user has developed an overall integrated pulse sequence that is desired to be run.
- scanning begins using the overall integrated pulse sequence built by the user as described above.
- motion estimates are generated at process block 412 .
- the motion estimates are sent back to the scanner as they are generated on the image reconstruction computer, and are immediately applied at process block 414 to the coordinates of the next imaging sequence to provide real-time and prospective correction of subject motion that is identified.
- This process is repeated until, at decision block 416 , all data has been acquired.
- This method is applicable both to three-dimensional acquisitions and to simultaneous multi-slice acquisitions, including when such acquisitions include in-place acceleration.
- a human volunteer was scanned in a 3T TIM Trio (Siemens Healthcare, Erlangen, Germany) using a 32-channel head matrix.
- the imaging sequence was a FLASH pulse sequence that used a 15 degree flip angle, 11 ms TR, 3.43 ms TE, 200 Hz/px bandwidth, 256 mm ⁇ 256 mm ⁇ 176 mm FOV, and 1 mm isotropic resolution, and 2 ⁇ GRAPPA acceleration for a total scan time of 4:48.
- One volume was acquired with this protocol while the subject remained still.
- vNav sequence was inserted after every 136 FLASH TRs (approximately 1.5 seconds), increasing the scan time to 5:38, and acquired another volume while the subject remained still. Then, two more volumes were acquired with the vNays sequences inserted, during both of which the subject was prompted to change their head position every minute, repeating the same motion pattern in both scans. In the first of these with-motion scans, correction was applied, using the process described above, for the subject's motion in real time, while in the second with-motion scan not apply the update.
- FIG. 5 a flowchart is illustrated as setting forth the steps of an example method for performing prospective motion correction based on an integrated pulse sequence that uses an ultra-fast navigator in an imaging pulse sequence to detect motion and trigger the performance of one or more vNav sequences.
- the method begins with user selection of the imaging pulse sequences and parameters for the sequence at process block 502 .
- the imaging pulse sequence may be a FLASH, SPGR, or other such gradient-echo pulse sequence.
- Selecting the imaging pulse sequences also includes selecting an ultra-fast navigator to be included in each imaging pulse sequence.
- the ultra-fast navigator may include a cloverleaf navigator.
- scanning begins using the selected imaging pulse sequences. Motion of the subject is monitored using the data acquired from the ultra-fast navigators acquired in each performance of the imaging pulse sequence, as indicated at step 506 . A determination is made at decision block 508 whether a predetermined amount of motion has occurred and been detected in the acquired navigator data. If not, imaging continues at step 504 and motion is monitored again at step 506 .
- vNav sequences are performed, as indicated at step 510 , as described above.
- motion estimates are generated at process block 512 .
- the motion estimates are sent back to the scanner as they are generated on the image reconstruction computer, and are immediately applied at process block 514 to the coordinates of the next imaging sequence to provide real-time and prospective correction of subject motion that is identified. This process is repeated until, at decision block 516 , all data has been acquired.
- This method is applicable for both three-dimensional acquisitions and simultaneous multi-slice acquisitions, including when such acquisitions include in-place acceleration.
- a system and method for performing prospective motion correction in FLASH, SPGR, or FLASH/SPGR-like pulse sequences.
- the present disclosure provides a way to insert multiple matched vNav TRs into the overall pulse sequence train and, thus, preserves the steady state.
- the present disclosure recognizes that a navigator of substantial length, including vNays, can be inserted or interleaved in a dense sequence imaging sequence (like FLASH), by substituting FLASH TRs for vNav TRs, while preserving the magnetization steady state. For example, to preserve the magnetization steady state, it is advantageous to match the vNav and FLASH TR times.
- any gradients are used during the excitation pulse, it is advantageous to match these gradients.
- any RF spoiling between the vNav and FLASH TRs may be advantageously synchronized.
- the additional TRs to accommodate the vNav sequence in the integrated FLASH-vNav pulse sequence does increase overall scan time.
- this drawback is overcome by the high registration accuracy enabled by whole-head navigators readily achieved with the vNav sequence.
- the additional scan time is still significantly less than what would be required for an MR technologist to request patient compliance and rescan.
- the provided, integrated vNav-FLASH/SPGR sequence allows users the flexibility to trade-off scan time for tracking accuracy based on the needs of their application and the subject population involved.
- the additional scan time introduced by the vNav sequences is reduced by incorporating an ultra-fast navigator, such as a cloverleaf navigator, into each imaging pulse sequence.
- the navigator data is then relied upon to monitor motion, such that vNav sequences are performed only when a significant amount of subject motion is detected. In this manner, the benefits of the vNav sequences can be realized without adding unnecessary amounts of additional scan time.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Physiology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Power Engineering (AREA)
- Pulmonology (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
A system and method for reconstructing an image of a subject, in which motion has been prospectively corrected, using a magnetic resonance imaging (MRI] system are provided. An imaging pulse sequence is used to acquire image data in addition to navigator data using an ultra-fast navigator, such as a cloverleaf navigator. When motion is detected in the navigator data, a volume navigator (vNav] sequence is performed, from which motion estimates are generated and used to update the imaging pulse sequence to prospectively correct for the motion in subsequent repetition times.
Description
- This application is based on, claims priority to, and incorporates herein by reference in its entirety, U.S. Provisional Application Ser. No. 61/981,452, filed Apr. 18, 2014, and entitled “SYSTEM AND METHOD FOR MAGNETIC RESONANCE IMAGING WITH PROSPECTIVE MOTION CONTROL” and is based on, claims priority to, and incorporates herein by reference in its entirety, U.S. Provisional Application Ser. No. 62/005,165, filed May 30, 2014, and entitled “SYSTEM AND METHOD FOR MAGNETIC RESONANCE IMAGING WITH PROSPECTIVE MOTION CONTROL.”
- This invention was made with government support under HD074649 awarded by the National Institutes of Health. The government has certain rights in the invention.
- The present disclosure relates to systems and methods for magnetic resonance imaging (“MRI”). More particularly, the disclosure relates to systems and methods for tracking and controlling artifacts caused by motion during a MRI procedure.
- MRI uses the nuclear magnetic resonance (“NMR”) phenomenon to produce images. When a substance such as human tissue is subjected to a uniform magnetic field (“main magnetic field”), B0, the individual magnetic moments of the nuclei in the tissue attempt to align with this magnetic field, but precess about it in random order at their characteristic Larmor frequency, ω. If the substance, or tissue, is subjected to an excitation magnetic field, B1, that is in the plane transverse to the main magnetic field, B0, and that is near the Larmor frequency, ω, the net aligned magnetic moment of the nuclei may be rotated, or “tipped,” into the transverse plane to produce a net transverse magnetic moment. A signal is emitted by the excited nuclei, or “spins,” after the excitation magnetic field, B1, is terminated. The emitted signal may be received and processed to form an image.
- When utilizing these emitted “MR” signals to produce images, magnetic field gradients (Gx, Gy, and Gz) are employed. Typically, the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used. The resulting set of received MR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
- The measurement cycle used to acquire each MR signal is performed under the direction of a pulse sequence produced by a pulse sequencer. Clinically available MRI systems store a library of such pulse sequences that can be prescribed to meet the needs of many different clinical applications. Research MRI systems include a library of clinically-proven pulse sequences and they also enable the development of new pulse sequences.
- Depending on the technique used, many MR scans currently require many minutes to acquire the necessary data used to produce medical images. The reduction of this scan time is an important consideration, since reduced scan time increases patient throughout, improves patient comfort, and improves image quality by reducing motion artifacts. Many different strategies have been developed to shorten the scan time.
- For example, one popular category of pulse sequences are so-called gradient echo sequences. Within this category, the spoiled gradient echo and three-dimensional spoiled gradient-recalled echo (SPGR) or fast low angle shoot (FLASH) pulse sequences are often used in neuroimaging applications. Specifically, the SPGR or FLASH sequence forms the basis of many 3D neuroimaging sequences, but acquisition times often stretch to several minutes. Lengthy acquisitions in neuroimaging applications can be particularly troublesome because it is imperative that the subject remain motionless during the duration of the sequence. That is, in neuroimaging applications, motion can be particularly damaging to the resulting images because of the complexity of the structures being imaged and studied in neuroimaging applications.
- Motion-correction systems in MRI can be grouped into two general methods: prospective and retrospective. Retrospective methods use information about the subject's motion to estimate what k-space data would have been measured if the subject had not moved during scanning. Prospective methods use motion-tracking data acquired during the scan to follow the subject with the gradient axes of the sequence, measuring the desired k-space data directly. Additionally, it is possible to combine the two methods so that retrospective processing corrects residual errors in the prospective system. A retrospective system can access all of the k-space data while performing reconstruction; a prospective system must necessarily rely only on previous measurements to estimate the current position of the patient. However, a prospective system avoids the need to estimate missing k-space data, allowing for direct reconstruction while avoiding possible sources of estimation error in the k-space data.
- Also, one can differentiate between two types of motion correction problems that arise in MRI: between-scan motion and within-scan motion. For between-scan motion, several retrospective motion correction methods are available that register either slice-by-slice or volume-by-volume to estimate the data that would have been acquired in each volume if the subject had not moved. Prospective motion correction can also be employed for this problem, such as the orbital navigator system that inserts 3-plane circular k-space navigators, or the PACE system that registers each completed EPI volume back to the first time-point and so requires no navigators.
- For in-scan motion, several methods are available, such as PROPELLER and the like that use redundant sampling of the center of k-space during each repetition time (TR) and estimate motion-free k-space data retrospectively. Also, prospective motion correction is available, such as by using cloverleaf navigators. The use of cloverleaf navigators is useful with SPGR/FLASH sequences.
- However, in order to maximize SNR/time efficiency, short TR protocols are often used in neuroimaging applications. Such short-TR protocols, by definition, have very-little dead time and, thus, force navigators to be very short and provide limited k-space coverage. That is, as the TR is reduced, the effectiveness of the navigator is reduced because there is less information gathered by the navigator to use to form an estimate of the subject's head motion.
- It would therefore be desirable to provide a system and method for controlling the competing constraints of neuroimaging applications that desire extended acquisition times and the need to control or compensate or correct for patient motion during such acquisitions.
- The present disclosure overcomes the aforementioned drawbacks by providing a system and method for directing a magnetic resonance imaging (MRI) system to reconstruct an image of a subject, in which motion has been prospectively corrected. The method includes acquiring image data and navigator data from a subject by directing the MRI system to perform an imaging pulse sequence that includes a navigator portion. The navigator data is processed to determine whether motion occurred while the image data was acquired. The imaging pulse sequence is updated when processing the navigator data determines motion occurred. The imaging pulse sequence is updated by first acquiring volume navigator data by directing the MRI system to perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject. The volume navigator data is then processed to generate motion estimates of motion in the VOI. The imaging pulse sequence is then updated based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates. The imaging pulse sequence, whether updated or not, is then repeated and navigator data processed to determine whether motion occurred and whether the imaging pulse sequence should be further updated until a desired amount of image data is acquired. An image of the subject is then reconstructed from the acquired image data.
- In accordance with yet another aspect of the disclosure, a magnetic resonance imaging (MRI) system is disclosed that includes a magnet system configured to generate a polarizing magnetic field about at least a portion of a subject arranged in the MRI system, a magnetic gradient system including a plurality of magnetic gradient coils configured to apply at least one magnetic gradient field to the polarizing magnetic field, and radio frequency (RF) system configured to apply an RF field to the subject and to receive magnetic resonance signals from the subject using a coil array. The system also includes a computer system programmed to control the gradient system and the RF system to perform an imaging pulse sequence that includes a navigator portion to acquire image data and navigator data from a subject. The computer system is further programmed to process the navigator data to determine whether the subject moved while the image data was acquired. Upon determining that the subject moved, the computer system is configured to control the gradient system and the RF system to perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject, process the volume navigator data to generate motion estimates of motion in the VOI, and update the imaging pulse sequence based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates. The computer system is further configured to reconstruct an image of the subject from the acquired image 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.
-
FIG. 1 is a block diagram of an example of a magnetic resonance imaging (“MRI”) system. -
FIG. 2 is an example of a three-dimensional spoiled gradient-recalled echo (“3D SPGR”) pulse sequence. -
FIG. 3 is an illustrative example of an integrated pulse sequence that includes both imaging sequences, each played out during an imaging repetition time (“TR”), and volumetric navigator sequences, each played out during a volumetric navigator TR; -
FIG. 4 is a flowchart setting forth the steps of an example method for performing prospective motion correction using a pulse sequence such as the one described with respect toFIG. 3 ; and -
FIG. 5 is a flowchart setting forth the steps of an example method for performing prospective motion correction based on an integrated pulse sequence that uses an ultra-fast navigator in an imaging pulse sequence to detect motion and trigger the performance of one or more vNav sequences. - Referring particularly now to
FIG. 1 , an example of a magnetic resonance imaging (“MRI”)system 100 is illustrated. TheMRI system 100 includes anoperator workstation 102, which will typically include adisplay 104; one ormore input devices 106, such as a keyboard and mouse; and aprocessor 108. Theprocessor 108 may include a commercially available programmable machine running a commercially available operating system. Theoperator workstation 102 provides the operator interface that enables scan prescriptions to be entered into theMRI system 100. In general, theoperator workstation 102 may be coupled to four servers: apulse sequence server 110; adata acquisition server 112; adata processing server 114; and adata store server 116. Theoperator workstation 102 and eachserver servers - The
pulse sequence server 110 functions in response to instructions downloaded from theoperator workstation 102 to operate agradient system 118 and a radiofrequency (“RF”)system 120. Gradient waveforms necessary to perform the prescribed scan are produced and applied to thegradient system 118, which excites gradient coils in anassembly 122 to produce the magnetic field gradients Gx, Gy, and Gz used for position encoding magnetic resonance signals. Thegradient coil assembly 122 forms part of amagnet assembly 124 that includes apolarizing magnet 126 and a whole-body RF coil 128. - RF waveforms are applied by the
RF system 120 to theRF coil 128, or a separate local coil (not shown inFIG. 1 ), in order to perform the prescribed magnetic resonance pulse sequence. Responsive magnetic resonance signals detected by theRF coil 128, or a separate local coil (not shown inFIG. 1 ), are received by theRF system 120, where they are amplified, demodulated, filtered, and digitized under direction of commands produced by thepulse sequence server 110. TheRF system 120 includes an RF transmitter for producing a wide variety of RF pulses used in MRI pulse sequences. The RF transmitter is responsive to the scan prescription and direction from thepulse sequence server 110 to produce RF pulses of the desired frequency, phase, and pulse amplitude waveform. The generated RF pulses may be applied to the whole-body RF coil 128 or to one or more local coils or coil arrays (not shown inFIG. 1 ). - The
RF system 120 also includes one or more RF receiver channels. Each RF receiver channel includes an RF preamplifier that amplifies the magnetic resonance signal received by thecoil 128 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, therefore, 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 (I 2 +Q 2)} (1); - and the phase of the received magnetic resonance signal may also be determined according to the following relationship:
-
- The
pulse sequence server 110 also optionally receives patient data from aphysiological acquisition controller 130. By way of example, thephysiological acquisition controller 130 may receive signals from a number of different sensors connected to the patient, such as electrocardiograph (“ECG”) signals from electrodes, or respiratory signals from a respiratory bellows or other respiratory monitoring device. Such signals are typically used by thepulse sequence server 110 to synchronize, or “gate,” the performance of the scan with the subject's heart beat or respiration. - The
pulse sequence server 110 also connects to a scanroom interface circuit 132 that receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scanroom interface circuit 132 that apatient positioning system 134 receives commands to move the patient to desired positions during the scan. - The digitized magnetic resonance signal samples produced by the
RF system 120 are received by thedata acquisition server 112. Thedata acquisition server 112 operates in response to instructions downloaded from theoperator workstation 102 to receive the real-time magnetic resonance data and provide buffer storage, such that no data is lost by data overrun. In some scans, thedata acquisition server 112 does little more than pass the acquired magnetic resonance data to thedata processor server 114. However, in scans that require information derived from acquired magnetic resonance data to control the further performance of the scan, thedata acquisition server 112 is programmed to produce such information and convey it to thepulse sequence server 110. For example, during prescans, magnetic resonance data is acquired and used to calibrate the pulse sequence performed by thepulse sequence server 110. As another example, navigator signals may be acquired and used to adjust the operating parameters of theRF system 120 or thegradient system 118, or to control the view order in which k-space is sampled. In still another example, thedata acquisition server 112 may also be employed to process magnetic resonance signals used to detect the arrival of a contrast agent in a magnetic resonance angiography (“MRA”) scan. By way of example, thedata acquisition server 112 acquires magnetic resonance data and processes it in real-time to produce information that is used to control the scan. - The
data processing server 114 receives magnetic resonance data from thedata acquisition server 112 and processes it in accordance with instructions downloaded from theoperator workstation 102. Such processing may, for example, include one or more of the following: reconstructing two-dimensional or three-dimensional images by performing a Fourier transformation of raw k-space data; performing other image reconstruction algorithms, such as iterative or backprojection reconstruction algorithms; applying filters to raw k-space data or to reconstructed images; generating functional magnetic resonance images; calculating motion or flow images; and so on. - Images reconstructed by the
data processing server 114 are conveyed back to theoperator workstation 102 where they are stored. Real-time images are stored in a data base memory cache (not shown inFIG. 1 ), from which they may be output tooperator display 112 or adisplay 136 that is located near themagnet assembly 124 for use by attending physicians. Batch mode images or selected real time images are stored in a host database ondisc storage 138. When such images have been reconstructed and transferred to storage, thedata processing server 114 notifies thedata store server 116 on theoperator workstation 102. Theoperator workstation 102 may be used by an operator to archive the images, produce films, or send the images via a network to other facilities. - The
MRI system 100 may also include one or morenetworked workstations 142. By way of example, anetworked workstation 142 may include adisplay 144; one ormore input devices 146, such as a keyboard and mouse; and aprocessor 148. Thenetworked workstation 142 may be located within the same facility as theoperator workstation 102, or in a different facility, such as a different healthcare institution or clinic. - The
networked workstation 142, whether within the same facility or in a different facility as theoperator workstation 102, may gain remote access to thedata processing server 114 ordata store server 116 via the communication system 140. Accordingly, multiplenetworked workstations 142 may have access to thedata processing server 114 and thedata store server 116. In this manner, magnetic resonance data, reconstructed images, or other data may be exchanged between thedata processing server 114 or thedata store server 116 and thenetworked workstations 142, such that the data or images may be remotely processed by anetworked workstation 142. This data may be exchanged in any suitable format, such as in accordance with the transmission control protocol (“TCP”), the internet protocol (“IP”), or other known or suitable protocols. - As described, systems such as described above with respect to
FIG. 1 have been used to perform neuroimaging acquisitions with short repetition times (“TRs”), which force navigators to be very short and, therefore, provide limited k-space coverage. These navigators are often referred to as “ultra-fast navigators.” The effectiveness of the navigators at providing an estimate of the subject's head motion is, therefore, also limited. - In contrast to these ultra-fast navigators, echo-planar imaging (EPI)-based navigators, or volume navigators (“vNays”), have been used to acquire a whole-head volume in roughly 275 ms and, thereby, allow high-accuracy motion tracking. However, vNays have previously only been used in sequences with significant dead time provided by inflow times (“TI”) or TR gaps in which the entire vNav could be inserted. As will be described, however, the present disclosure provides a way to insert a vNav into a 3D FLASH/SPGR sequence with only marginal impact on SNR or time. By doing so, the present disclosure provides a system and method for prospective motion correction of in-scan motion using a 3D FLASH/SPGR pulse sequence with vNays that does not substantially extend acquisition times.
- An example of a pulse sequence employed to direct an MRI system to acquire image data in accordance with some configurations of the present invention is illustrated in
FIG. 2 . The pulse sequence includes anRF excitation pulse 202 that is played out in the presence of a slab-selective gradient 204 in order to produce transverse magnetization in a volume-of-interest. The slab-selective gradient 204 includes arephasing lobe 206 that acts to rephase unwanted phase dispersions introduced by the slab-selective gradient 204 such that signal losses resultant from these phase dispersions are mitigated. - Following excitation of the nuclear spins in the volume-of-interest, a phase-
encoding gradient 208 is applied to spatially encode a nuclear magneticresonance echo signal 210 at a given echo time (TE) representative of a gradient-recalled echo along one direction in the volume-of-interest. At the same time, a partition-encoding gradient 212 is applied to spatially encode theecho signal 210 along a second, orthogonal direction in the volume-of-interest. By way of example, the phase-encoding gradient 208 may spatially encode theecho signal 210 along the y-direction, while the partition-encoding gradient 212 may spatially encode theecho signal 210 along the z-direction. Areadout gradient 214 is also applied after adephasing gradient lobe 216 to spatially encode theecho signal 210 along a third, orthogonal direction in the volume-of-interest. By way of example, thereadout gradient 214 may spatially encode the echo signal along the x-direction. Theecho signal 210 is sampled during a data acquisition window. -
Spoiler gradients - Referring now to
FIG. 3 , the pulse sequence descried above with respect toFIG. 2 can, in accordance with the present disclosure, be combined with vNav sequences to yield a pulse sequence for performing prospective motion correction in SPGR/FLASH without extended TI or TR gaps. Specifically,FIG. 3 provides a timing diagram of an integrated vNav and a SPGR/FLASH sequence 300, showing common pulses and TRs, but varied gradients, for interleaved imaging and vNav sequences. In particular, a plurality ofimaging sequences 302, such as described above with respect toFIG. 2 are illustrated. Integrated therewith are a plurality ofvNav sequences 304. As will be described, the scan time of theimaging sequence 302 is selected such that animaging TR 306 allowsvNav sequences 304, each having avNav TR 308, to be inserted at any desired point. Furthermore, the RF pulses used in thevNav sequences 304 and theimaging sequences 302 may be the same. To this end, the steady-state of theimaging sequences 302 is maintained, even whenvNav sequences 304 are inserted in theintegrated sequence 300. - As an example, the
vNav sequence 304 may be a 3D-encoded echo-planar imaging (EPI) pulse sequence with a 323 matrix and may be acquired with ¾ partial Fourier encoding in the partition direction. As such, in this example, thevNav sequence 304 may have 25 three-dimensional excitation pulses. In this example, thevNav TR 308 is assumed to be 11 ms and theTR 306 for eachimaging sequence 302 is assumed to match thevNav TR 308. If the pulses of bothsequences TR 306 of animaging sequence 302 without disturbing the steady state of theimaging sequence 302. Expanding on this example, as long as theTR 306 of theimaging sequence 302 is 11 ms or longer, matching thevNav TR 308 to theimaging TR 306 can be readily achieved by addingTRs 310 to thevNav sequence 304 after its readouts and not simply attempting to fit thevNav TR 308 within dead time. Again, as long as the scan time of theimaging sequence 302 meets the minimum TR requirement,vNav TRs 308 can be inserted at any desired point while maintaining the steady-state of the imaging sequence. - Inserting the
vNav sequences 304 asseparate TRs 308 at any time allows a great deal of flexibility, but also increases overall scan time. In previous acquisitions using vNav sequences, the navigators were inserted in TR or TI gaps and, thus, did not increase scan time. However, continuing with the non-limiting example provided above, 25 additional TRs are added to the total scan time with every vNav sequence while not gaining any additional imaging signal. To allow grater flexibility, as will be described, the user may set how many imaging TRs will be played between each vNav sequence, and the overall scan time is updated to inform users of the trade-off between tracking accuracy and scan time. - In some embodiments, however, an additional time savings can be realized by selectively performing a limited number of
vNav sequences 304 such that thevNav sequences 304 are performed only when a predetermined amount of motion is occurring. In these instances, motion detection can be provided by inserting an ultra-fast navigator, such as a cloverleaf navigator, into eachimaging sequence 302. Examples of cloverleaf navigators are described in U.S. Pat. Nos. 6,771,068 and 6,958,605, both of which are herein incorporated by reference in their entirety. The ultrafast-navigator can be used as a high-sensitivity motion detector that can be relied on to determine when motion occurs without a significant scan time cost. When motion is detected, one or morevNav sequences 304 can be performed to provide a more accurate measurement of the motion, as described above in detail. In this manner, the more time consumingvNav sequences 304 can be performed only when a predetermined amount of motion is detected by an ultra-fast navigator, thereby reducing overall scan time by eliminating the performance of unnecessaryvNav sequences 304. Integrating an ultra-fast navigator into theimaging sequences 302 also has the added benefit of enabling resonance frequency correction. - Specifically, referring to
FIG. 4 , a flowchart is provided that sets forth an example of a process for performing prospective motion correction using the pulse sequence described with respect toFIG. 3 . The prospective motion-correctedprocess 400 begins with user selection of the imaging pulse sequences and parameters for the sequence atprocess block 402. As described above, the imaging pulse sequence may be a FLASH or SPGR or other name for such gradient-echo sequences. Once the imaging sequence and parameters are selected, atprocess block 404, the user is prompted to select how many imaging TRs will be played between each vNav TR. Using this information, atprocess block 406, the overall scan time is updated to inform users of the trade-off between tracking accuracy and scan time. As indicated bydecision block 408, the user may adjust the number and spacing of the vNav sequences to reach a desired overall scan time until the user has developed an overall integrated pulse sequence that is desired to be run. - At
process block 410, scanning begins using the overall integrated pulse sequence built by the user as described above. Using the data acquired from each performance of the vNav pulse sequence, motion estimates are generated atprocess block 412. Specifically, the motion estimates are sent back to the scanner as they are generated on the image reconstruction computer, and are immediately applied at process block 414 to the coordinates of the next imaging sequence to provide real-time and prospective correction of subject motion that is identified. This process is repeated until, atdecision block 416, all data has been acquired. This method is applicable both to three-dimensional acquisitions and to simultaneous multi-slice acquisitions, including when such acquisitions include in-place acceleration. - To test the efficacy of the above-described process, a human volunteer was scanned in a 3T TIM Trio (Siemens Healthcare, Erlangen, Germany) using a 32-channel head matrix. The imaging sequence was a FLASH pulse sequence that used a 15 degree flip angle, 11 ms TR, 3.43 ms TE, 200 Hz/px bandwidth, 256 mm×256 mm×176 mm FOV, and 1 mm isotropic resolution, and 2×GRAPPA acceleration for a total scan time of 4:48. One volume was acquired with this protocol while the subject remained still. Thereafter, a vNav sequence was inserted after every 136 FLASH TRs (approximately 1.5 seconds), increasing the scan time to 5:38, and acquired another volume while the subject remained still. Then, two more volumes were acquired with the vNays sequences inserted, during both of which the subject was prompted to change their head position every minute, repeating the same motion pattern in both scans. In the first of these with-motion scans, correction was applied, using the process described above, for the subject's motion in real time, while in the second with-motion scan not apply the update.
- Referring now to
FIG. 5 , a flowchart is illustrated as setting forth the steps of an example method for performing prospective motion correction based on an integrated pulse sequence that uses an ultra-fast navigator in an imaging pulse sequence to detect motion and trigger the performance of one or more vNav sequences. The method begins with user selection of the imaging pulse sequences and parameters for the sequence atprocess block 502. As described above, the imaging pulse sequence may be a FLASH, SPGR, or other such gradient-echo pulse sequence. Selecting the imaging pulse sequences also includes selecting an ultra-fast navigator to be included in each imaging pulse sequence. For example, the ultra-fast navigator may include a cloverleaf navigator. - At
process block 504, scanning begins using the selected imaging pulse sequences. Motion of the subject is monitored using the data acquired from the ultra-fast navigators acquired in each performance of the imaging pulse sequence, as indicated atstep 506. A determination is made atdecision block 508 whether a predetermined amount of motion has occurred and been detected in the acquired navigator data. If not, imaging continues atstep 504 and motion is monitored again atstep 506. - If, however, a predetermined amount of motion is detected, then one or more vNav sequences are performed, as indicated at
step 510, as described above. Using data acquired from the vNav pulse sequences, motion estimates are generated atprocess block 512. Specifically, the motion estimates are sent back to the scanner as they are generated on the image reconstruction computer, and are immediately applied at process block 514 to the coordinates of the next imaging sequence to provide real-time and prospective correction of subject motion that is identified. This process is repeated until, atdecision block 516, all data has been acquired. This method is applicable for both three-dimensional acquisitions and simultaneous multi-slice acquisitions, including when such acquisitions include in-place acceleration. - Thus, a system and method is provided for performing prospective motion correction in FLASH, SPGR, or FLASH/SPGR-like pulse sequences. Unlike previous navigator methods that attempted to insert ultra-fast navigators into the imaging (FLASH/SPGR) TR, the present disclosure provides a way to insert multiple matched vNav TRs into the overall pulse sequence train and, thus, preserves the steady state. The present disclosure recognizes that a navigator of substantial length, including vNays, can be inserted or interleaved in a dense sequence imaging sequence (like FLASH), by substituting FLASH TRs for vNav TRs, while preserving the magnetization steady state. For example, to preserve the magnetization steady state, it is advantageous to match the vNav and FLASH TR times. Also, it is advantageous to ensure that the waveform of the excitation pulses match between the vNav and FLASH TRs. Further, if any gradients are used during the excitation pulse, it is advantageous to match these gradients. In addition, any RF spoiling between the vNav and FLASH TRs may be advantageously synchronized.
- Unlike in previous applications of vNays, where there was no scan-time increase, the additional TRs to accommodate the vNav sequence in the integrated FLASH-vNav pulse sequence does increase overall scan time. However, this drawback is overcome by the high registration accuracy enabled by whole-head navigators readily achieved with the vNav sequence. Furthermore, the additional scan time is still significantly less than what would be required for an MR technologist to request patient compliance and rescan. Additionally, the provided, integrated vNav-FLASH/SPGR sequence allows users the flexibility to trade-off scan time for tracking accuracy based on the needs of their application and the subject population involved.
- In some configurations, the additional scan time introduced by the vNav sequences is reduced by incorporating an ultra-fast navigator, such as a cloverleaf navigator, into each imaging pulse sequence. The navigator data is then relied upon to monitor motion, such that vNav sequences are performed only when a significant amount of subject motion is detected. In this manner, the benefits of the vNav sequences can be realized without adding unnecessary amounts of additional scan time.
- 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.
Claims (12)
1. A method for reconstructing an image of a subject, in which motion has been prospectively corrected, using a magnetic resonance imaging (MRI) system, the steps of the method comprising:
(a) acquiring image data and navigator data from a subject by directing the MRI system to perform an imaging pulse sequence that includes a navigator portion;
(b) processing the navigator data to determine whether motion occurred while the image data was acquired;
(c) upon determining motion occurred in step (b), updating the imaging pulse sequence by:
(i) directing the MRI system to perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject;
(ii) processing the volume navigator data to generate motion estimates of motion in the VOI;
(iii) updating the imaging pulse sequence based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates;
(d) repeating steps (a) through (c) until a desired amount of image data is acquired; and
(e) reconstructing an image of the subject from the acquired image data.
2. The method as recited in claim 1 , wherein the imaging pulse sequence includes a gradient echo pulse sequence.
3. The method as recited in claim 1 , wherein the navigator portion of the imaging pulse sequence includes performing a cloverleaf navigator.
4. The method as recited in claim 1 , wherein the image data is acquired during an imaging repetition time (TR) and the volume navigator data is acquired during a volume navigator TR that is coordinated with the imaging TR to preserve a steady-state of magnetization in the VOI associated with the imaging pulse sequence.
5. The method as recited in claim 4 , wherein the volume navigator TR is selected to match the imaging TR.
6. The method as recited in claim 1 , wherein the VOI includes a head of the subject, and the volume navigator pulse sequence is configured to acquire the volume navigator data from the head.
7. A magnetic resonance imaging (MRI) system, comprising:
a magnet system configured to generate a polarizing magnetic field about at least a portion of a subject arranged in the MRI system;
a magnetic gradient system including a plurality of magnetic gradient coils configured to apply at least one magnetic gradient field to the polarizing magnetic field;
a radio frequency (RF) system configured to apply an RF field to the subject and to receive magnetic resonance signals from the subject using a coil array;
a computer system programmed to:
control the gradient system and the RF system to perform an imaging pulse sequence that includes a navigator portion to acquire image data and navigator data from a subject;
process the navigator data to determine whether the subject moved while the image data was acquired;
upon determining that the subject moved, control the gradient system and the RF system to:
perform a volume navigator pulse sequence that uses an echo-planar imaging (EPI) technique to acquire volume navigator data from a volume-of-interest in the subject;
process the volume navigator data to generate motion estimates of motion in the VOI;
update the imaging pulse sequence based on the motion estimates, wherein the updated imaging pulse sequence prospectively corrects for the motion in the VOI using the motion estimates; and
reconstruct an image of the subject from the acquired image data.
8. The system as recited in claim 7 , wherein the imaging pulse sequence includes a gradient echo pulse sequence.
9. The system as recited in claim 7 , wherein the navigator portion of the imaging pulse sequence includes performing a cloverleaf navigator.
10. The system as recited in claim 7 , wherein the image data is acquired during an imaging repetition time (TR) and the volume navigator data is acquired during a volume navigator TR that is coordinated with the imaging TR to preserve a steady-state of magnetization in the VOI associated with the imaging pulse sequence.
11. The method as recited in claim 10 , wherein the volume navigator TR is selected to match the imaging TR.
12. The method as recited in claim 7 , wherein the VOI includes a head of the subject, and the volume navigator pulse sequence is configured to acquire the volume navigator data from the head.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/301,868 US20170160366A1 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance i maging with prospective motion control |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201461981452P | 2014-04-18 | 2014-04-18 | |
US201462005165P | 2014-05-30 | 2014-05-30 | |
PCT/US2015/026159 WO2015161067A1 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance imaging with prospective motion control |
US15/301,868 US20170160366A1 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance i maging with prospective motion control |
Publications (1)
Publication Number | Publication Date |
---|---|
US20170160366A1 true US20170160366A1 (en) | 2017-06-08 |
Family
ID=54324573
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/301,868 Abandoned US20170160366A1 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance i maging with prospective motion control |
US15/301,860 Active 2036-01-23 US10459058B2 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance imaging with prospective motion control |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/301,860 Active 2036-01-23 US10459058B2 (en) | 2014-04-18 | 2015-04-16 | System and method for magnetic resonance imaging with prospective motion control |
Country Status (2)
Country | Link |
---|---|
US (2) | US20170160366A1 (en) |
WO (2) | WO2015161059A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10302724B2 (en) | 2016-01-19 | 2019-05-28 | Siemens Healthcare Gmbh | Multi-slice gradient echo magnetic resonance imaging |
US11474173B2 (en) * | 2020-01-03 | 2022-10-18 | Siemens Healthcare Gmbh | Magnetic resonance apparatus and method for operating a magnetic resonance apparatus, computer program and electronically readable data storage medium |
US11754654B2 (en) * | 2017-06-27 | 2023-09-12 | Koninklijke Philips N.V. | Method and device for determining a motion field from k-space data |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10304219B2 (en) * | 2014-06-13 | 2019-05-28 | Siemens Medical Solutions Usa, Inc. | Intra reconstruction motion correction |
DE102015222835B4 (en) * | 2015-11-19 | 2019-06-13 | Siemens Healthcare Gmbh | Magnetic resonance imaging method with simultaneous image acquisition of several partial volumes with a synchronous image acquisition by navigators |
US20200146554A1 (en) * | 2015-12-14 | 2020-05-14 | Danmarks Tekniske Universitet | Scanning and tracking monitoring apparatus and method |
EP3595525B1 (en) * | 2017-03-17 | 2023-09-27 | The Regents of The University of California | Systems and methods for real-time motion prediction in dynamic imaging |
US11484279B2 (en) | 2018-09-24 | 2022-11-01 | Siemens Medical Solutions Usa, Inc. | Systems to assess projection data inconsistency |
US10806370B1 (en) | 2019-04-25 | 2020-10-20 | General Electric Company | MRI system and method for detection and correction of patient motion |
CN109917315B (en) * | 2019-04-30 | 2021-09-28 | 上海联影医疗科技股份有限公司 | Magnetic resonance imaging scanning method, magnetic resonance imaging scanning device, computer equipment and storage medium |
US11519992B2 (en) | 2020-08-26 | 2022-12-06 | GE Precision Healthcare LLC | Systems and methods of automatically detecting and scoring motion artifacts in magnetic resonance imaging |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5800354A (en) * | 1994-11-23 | 1998-09-01 | U.S. Phillips Corporation | Method of and device for magnetic resonance imaging |
US20030135105A1 (en) * | 2000-04-26 | 2003-07-17 | Jack Clifford R. | Alignment of multiple MR images using navigator signals |
US6980847B2 (en) * | 2002-03-25 | 2005-12-27 | New York State Department Of Mental Hygiene Office Of Mental Health | Apparatus for dynamic angular position tracking |
US7358732B2 (en) | 2005-10-24 | 2008-04-15 | The General Hospital Corporation | System, method, software arrangement and computer-accessible medium for providing real-time motion correction by utilizing clover leaf navigators |
US8406849B2 (en) | 2006-03-31 | 2013-03-26 | University Of Utah Research Foundation | Systems and methods for magnetic resonance imaging |
JP5305925B2 (en) | 2009-01-05 | 2013-10-02 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Magnetic resonance imaging system |
US10371779B2 (en) * | 2009-03-09 | 2019-08-06 | The Johns Hopkins University | Apparatus and method for magnetic resonance imaging with high spatial temporal resolutions |
JP5388749B2 (en) * | 2009-08-11 | 2014-01-15 | 株式会社東芝 | Magnetic resonance imaging system |
US8427153B2 (en) * | 2010-01-15 | 2013-04-23 | Beth Israel Deaconess Medical Center, Inc. | Method for motion correction in magnetic resonance imaging using radio frequency coil arrays |
US8981776B2 (en) * | 2011-04-22 | 2015-03-17 | The General Hospital Corporation | Method for magnetic resonance imaging with controlled aliasing |
US9687172B2 (en) * | 2012-01-24 | 2017-06-27 | National Institute of Health (NIH), The United States of America, U.S. Dept. of Health and Human Services (DHHS) | System for motion corrected MR diffusion imaging |
EP2626718A1 (en) * | 2012-02-09 | 2013-08-14 | Koninklijke Philips Electronics N.V. | MRI with motion correction using navigators acquired using a Dixon technique |
-
2015
- 2015-04-16 US US15/301,868 patent/US20170160366A1/en not_active Abandoned
- 2015-04-16 WO PCT/US2015/026144 patent/WO2015161059A1/en active Application Filing
- 2015-04-16 WO PCT/US2015/026159 patent/WO2015161067A1/en active Application Filing
- 2015-04-16 US US15/301,860 patent/US10459058B2/en active Active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10302724B2 (en) | 2016-01-19 | 2019-05-28 | Siemens Healthcare Gmbh | Multi-slice gradient echo magnetic resonance imaging |
US11754654B2 (en) * | 2017-06-27 | 2023-09-12 | Koninklijke Philips N.V. | Method and device for determining a motion field from k-space data |
US11474173B2 (en) * | 2020-01-03 | 2022-10-18 | Siemens Healthcare Gmbh | Magnetic resonance apparatus and method for operating a magnetic resonance apparatus, computer program and electronically readable data storage medium |
Also Published As
Publication number | Publication date |
---|---|
US20170219674A1 (en) | 2017-08-03 |
WO2015161067A1 (en) | 2015-10-22 |
US10459058B2 (en) | 2019-10-29 |
WO2015161059A1 (en) | 2015-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10459058B2 (en) | System and method for magnetic resonance imaging with prospective motion control | |
JP5980126B2 (en) | Simultaneous multi-slice magnetic resonance imaging using single and multi-channel receive coils | |
US9726742B2 (en) | System and method for iteratively calibrated reconstruction kernel for accelerated magnetic resonance imaging | |
US10261155B2 (en) | Systems and methods for acceleration magnetic resonance fingerprinting | |
US8957681B2 (en) | System and method for magnetic resonance imaging water-fat separation with full dynamic range using in-phase images | |
US10520570B2 (en) | System and method for tissue characterization using multislice magnetic resonance imaging | |
US10531812B2 (en) | System and method for improved cardiac imaging of subjects with adverse cardiac conditions | |
US8975895B2 (en) | Method for correcting motion-induced phase errors in magnetic resonance imaging | |
US10191132B2 (en) | Method and apparatus for multi-slice imaging of T2-relaxation time | |
US10338180B2 (en) | System and method for gradient measurement using single-point imaging | |
EP3385743A1 (en) | System and method for phase cycling magnetic resonance fingerprinting (phc-mrf) | |
US10045741B2 (en) | System and method for real-time frequency correction for magnetic resonance imaging | |
US10330760B2 (en) | System and method for assessing T2-relaxation times with improved accuracy | |
WO2023055709A1 (en) | System and method for controlling concomitant gradient phase errors in chemical‐shift encoded imaging | |
WO2017132182A1 (en) | Systems and methods for joint reconstruction of quantitative t2 and adc maps | |
US11266324B2 (en) | System and methods for fast multi-contrast magnetic resonance imaging | |
US11079448B2 (en) | System and method for proton density mapping and receiver bias correction using magnetic resonance fingerprinting (MRF) | |
US10401459B2 (en) | Systems and methods for imaging vascular calcifications with magnetic resonance imaging | |
WO2015152957A1 (en) | Inverse imaging with magnetic resonance imaging using blipped gradient encoding | |
US20240201302A1 (en) | Systems and methods for non-selective stimulated echo multislice diffusion imaging | |
US20070238974A1 (en) | Three-dimensional magnetization prepared mri acquisition with increased acquisition window |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: THE GENERAL HOSPITAL CORPORATION, MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TISDALL, MATTHEW D.;VAN DER KOUWE, ANDRE J.W.;SIGNING DATES FROM 20150501 TO 20150504;REEL/FRAME:039942/0896 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONMENT FOR FAILURE TO CORRECT DRAWINGS/OATH/NONPUB REQUEST |