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 PDF

Info

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
Application number
US15/301,868
Inventor
Matthew D. Tisdale
Andrew J.W. Van Der Kouwe
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Hospital Corp
Original Assignee
General Hospital Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by General Hospital Corp filed Critical General Hospital Corp
Priority to US15/301,868 priority Critical patent/US20170160366A1/en
Assigned to THE GENERAL HOSPITAL CORPORATION reassignment THE GENERAL HOSPITAL CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: TISDALL, MATTHEW D., VAN DER KOUWE, ANDRE J.W.
Publication of US20170160366A1 publication Critical patent/US20170160366A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5615Echo 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/5616Echo 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image 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/5673Gating 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

    CROSS-REFERENCE TO RELATED APPLICATIONS
  • 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.”
  • STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH
  • This invention was made with government support under HD074649 awarded by the National Institutes of Health. The government has certain rights in the invention.
  • BACKGROUND
  • 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.
  • SUMMARY
  • 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.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • 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 to FIG. 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.
  • DETAILED DESCRIPTION OF THE INVENTION
  • Referring particularly now to FIG. 1, an example of a magnetic resonance imaging (“MRI”) system 100 is illustrated. 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. In general, 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. For example, 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. As an example, 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 Gx, Gy, and Gz 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:

  • 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:
  • ϕ = tan 1 ( Q I ) . ( 2 )
  • The pulse sequence server 110 also optionally receives patient data from a physiological acquisition controller 130. By way of example, 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. 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. As another example, 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. In still another example, 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. By way of example, 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. 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 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. When such images have been reconstructed and transferred to storage, 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. By way of example, 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, whether within the same facility or in a different facility as the operator workstation 102, 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.
  • 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 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.
  • Following excitation of the nuclear spins in the volume-of-interest, 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. At the same time, a partition-encoding gradient 212 is applied to spatially encode the echo signal 210 along a second, orthogonal direction in the volume-of-interest. By way of example, the phase-encoding gradient 208 may spatially encode the echo signal 210 along the y-direction, while 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. By way of example, 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.
  • Referring now to FIG. 3, the pulse sequence descried above with respect to FIG. 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 of imaging sequences 302, such as described above with respect to FIG. 2 are illustrated. Integrated therewith are a plurality of vNav sequences 304. As will be described, 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. Furthermore, 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.
  • 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, the vNav sequence 304 may have 25 three-dimensional excitation pulses. In this example, 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. If the pulses of both sequences 302, 304 are matched, 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. Expanding on this example, as long as 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. Again, as long as the scan time of the imaging 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 as separate 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 the vNav 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 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. When motion is detected, 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.
  • 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 to FIG. 3. The prospective motion-corrected process 400 begins with user selection of the imaging pulse sequences and parameters for the sequence at process 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, at process block 404, the user is prompted to select how many imaging TRs will be played between each vNav TR. Using this information, at process block 406, the overall scan time is updated to inform users of the trade-off between tracking accuracy and scan time. As indicated by decision 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 at process 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, 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.
  • 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 at process 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 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.
  • 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 at process 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, 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.
  • 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.
US15/301,868 2014-04-18 2015-04-16 System and method for magnetic resonance i maging with prospective motion control Abandoned US20170160366A1 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Cited By (3)

* Cited by examiner, † Cited by third party
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