EP4594767A1 - Physiological navigation in mri using inherent slab-selective free-induction-decay signal - Google Patents
Physiological navigation in mri using inherent slab-selective free-induction-decay signalInfo
- Publication number
 - EP4594767A1 EP4594767A1 EP23793572.1A EP23793572A EP4594767A1 EP 4594767 A1 EP4594767 A1 EP 4594767A1 EP 23793572 A EP23793572 A EP 23793572A EP 4594767 A1 EP4594767 A1 EP 4594767A1
 - Authority
 - EP
 - European Patent Office
 - Prior art keywords
 - mri
 - navigator signal
 - sequence
 - data
 - navigator
 - 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.)
 - Pending
 
Links
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/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/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
 
 
Definitions
- MRI magnetic resonance imaging
 - the navigation data can be acquired using the MRI system separate and distinct from capturing the image data.
 - these techniques typically capture the navigation data in lieu of image data during normal acquisition timing cycles, which can increase the overall time to acquire the necessary image data.
 - one or more MRI sequence cycles are used to acquire the navigator signal and, when the navigator signal indicates that no respiratory motion is present, the MRI sequence cycles capturing the image data of the structure under observation commence. Reducing the timing required to capture the necessary MRI image data can improve the comfort of a patient, who does not have to be confined in an MRI machine for as long.
 - a technique for efficiently capturing a navigator signal during the same MRI sequence cycle as the image data is desired.
 - a navigator signal may be captured during a portion of the MRI pulse sequence after slice (or slab) selection and prior to k-space encoding for capturing image data using well-known MRI pulse sequence trajectories.
 - the navigator signal enables a subject being imaged to continue breathing normally during image acquisition as motion compensation can be performed retroactively using post-processing techniques (e.g., binning) based on the captured navigator signal.
 - a method applied using a magnetic resonance imaging (MRI) apparatus includes: capturing data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, processing the navigator signal to identify timing of a respiratory cycle of a subject being measured, and performing data reconstruction using the MRI apparatus based at least on the processed navigator signal.
 - Performing data reconstruction can include performing retrospective MRI data binning based on the processed navigator signal.
 - the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus.
 - RF radio frequency
 - the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence.
 - the slice select ramp down portion of the MRI sequence is defined as a period immediately subsequent to the RF pulse generation during which current in the at least one gradient coil is reduced to zero.
 - the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through the at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition.
 - processing the navigator signal comprises filtering the navigator signal via an angular k-space filter.
 - the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to a common respiration rate.
 - performing data acquisition using the MRI apparatus comprises gating image data according to a phase of the respiratory cycle identified using the processed navigator signal.
 - the MRI sequence comprises a two dimensional (2D) or three-dimensional (3D) imaging sequence using a slice-selective excitation.
 - the RF pulse comprises a slice- selective RF pulse.
 - the slice-selective RF pulse can include a sine or a ultra-short echo time (UTE) pulse.
 - the MRI apparatus includes a plurality of RF receivers coupled to a plurality of RF coils.
 - capturing the data for the navigator signal further comprises selecting an optimal RF coil from the plurality of coils to sample.
 - a magnetic resonance imaging (MRI) system includes: an MRI apparatus; and a control system connected to the MRI apparatus.
 - the control system includes at least one processor and a memory storing a set of instructions that, in response to being executed by the at least one processor, cause the MRI apparatus to: capture data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, process the navigator signal to identify timing of a respiratory cycle of a subject being measured, and perform data reconstruction using the MRI apparatus based at least on the processed navigator signal.
 - Performing data reconstruction can include performing data sorting and image formation based at least on the processed navigator signal.
 - the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus
 - the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence.
 - the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition.
 - processing the navigator signal comprises filtering the navigator signal via an angular k-space filter.
 - the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to common respiration rates.
 - performing data acquisition using the MRI apparatus comprises binning image data according to a phase of the respiratory cycle identified using the processed navigator signal.
 - performing data acquisition using the MRI apparatus comprises gating image data according to a phase of the respiratory cycle identified using the processed navigator signal.
 - the MRI sequence comprises a two-dimensional (2D) or three-dimensional (3D) imaging sequence using slice-selective excitation.
 - the RF pulse comprises a slice-selective RF pulse.
 - a non-transitory computer- readable storage medium storing instructions that, in response to being executed by one or more processors, cause a magnetic resonance imaging (MRI) apparatus to perform the method of the first aspect.
 - MRI magnetic resonance imaging
 - FIG. 1 illustrates a magnetic resonance imaging (MRI) system, in accordance with some embodiments of the present disclosure.
 - MRI magnetic resonance imaging
 - FIG. 2 illustrates a control system for the MRI system of FIG. 1, in accordance with some embodiments of the present disclosure.
 - FIG. 3A shows a timing diagram for an MRI sequence, in accordance with some embodiments of the present disclosure.
 - FIG. 3B shows the timing associated with three separate navigator signals, in accordance with some embodiments of the present disclosure.
 - FIG. 3C illustrates an MRI sequence corresponding with a 2D Cartesian trajectory, in accordance with some embodiments of the present disclosure.
 - FIG. 3D shows various signal processing steps applied to a captured raw navigator signal, in accordance with some embodiments of the present disclosure.
 - FIG. 3E illustrates the various navigator signals described herein, pre- and postprocessing, in accordance with some embodiments of the present disclosure.
 - FIG. 4 shows a set of short-axis slice images from a 3D stack-of-spirals ultra-short echo-time acquisition, in accordance with some embodiments.
 - FIG. 5 illustrates a comparison between the four different navigator signals conducted during an experiment using ten subjects, in accordance with some embodiments.
 - FIG. 6 is a flowchart of a method for acquiring MRI images using navigator signals, in accordance with some embodiments of the present disclosure.
 - Physiological information is used for MRI to guide acquisition timing for motion- resolved or motion-compensated imaging.
 - Existing means for generating signals to track a patient’s respiration cycle include use of devices on patients, such as electrodes for ECG (electrocardiogram) or pressure-sensitive respiratory bellows.
 - Hardware modifications are also possible for radio frequency (RF) sensing, but this in not universally available.
 - modem advanced techniques for acquisition of free-breathing MRI images have two means to detect physiological motion in the data collection itself: 1) use of a navigation-only signal acquisition to determine the liver position, and 2) designing the image acquisition such that there is frequent sampling of a "direct-current" signal, at a timing where there is no spatial encoding and high MRI signal, that fluctuates with respiration.
 - navigation segments can disrupt image acquisitions that rely on steady-state or fast imaging regimes, and implementing frequent sampling of "direct-current" signals can be difficult as well as not being possible to implement for every application.
 - the present disclosure describes acquisition of navigator signals by frequently sampling the signal emitted after the RF pulse and the slice select gradient refocus portion of the sequence, which has some unspecified spatial encoding - the "unintended" encoding present following a slice-selective pulse - used in the majority of MRI imaging techniques for choosing volumes of interest.
 - This unwanted signal is acquired during image acquisition "dead-time"
 - it is therefore an opportunity to acquire navigation information without disrupting the timing of the MRI sequence.
 - the technique described below can be implemented in parallel with image acquisition without disrupting the image data.
 - MRI applications used for determining respiratory motion, cardiac motion, and bulk subject motion, applied to cardiac imaging, pulmonary imaging, abdominal imaging, neuroimaging, musculoskeletal imaging are all within the scope of the present disclosure.
 - FIG. 1 illustrates an MRI system 100, in accordance with some embodiments of the present disclosure.
 - Magnetic resonance (MR) measurements can be obtained using an MRI apparatus 100 as illustrated in FIG. 1.
 - the MRI apparatus 100 includes a controller/interface 102 that is configured to apply selected magnetic fields such as constant or pulsed field gradients to an object (e.g., a patient) located proximate to the MRI apparatus 100.
 - An axial controller 104 is in communication with an axial magnet 106 that is generally configured to produce a substantially constant magnetic field Bo proximate the MRI apparatus 100.
 - the axial magnet 106 may comprise a coil that includes superconducting material such as niobium-tin (Nb3Sn) or Magnesium Diboride (MgB2).
 - Nb3Sn niobium-tin
 - MgB2 Magnesium Diboride
 - a gradient controller 108 is configure to apply a constant or timevarying magnetic field gradient in a selected direction using a number of gradient coils 110.
 - the gradient controller 108 may cause electrical current to be applied to the gradient coils 110 to apply diffusion or slice selection gradients to the object.
 - the magnetic field gradient vector components may be referred to as Gx, Gy, and Gz corresponding to an x-axis, y-axis, and z-axis, respectively.
 - each axis is associated with a separate gradient coil, e.g., a first gradient coil 110-1 is associated with an x-axis, a second gradient coil 110-2 is associated with a y-axis, and a third gradient coil 110-3 is associated with a z-axis.
 - a magnetic field gradient may be oriented in any direction, meaning the magnitude of the magnetic field changes along an axis corresponding to that direction, such that two or more magnetic field gradient vector components are non-zero.
 - An RF generator 114 is configured to generate one or more RF pulses of the electromagnetic field to the object being imaged.
 - the RF generator 114 is connected to one or more transmitter coils 115, and delivers an electrical current to the transmitter coils 115 to generate the RF pulse.
 - One or more RF receivers 116 are coupled to corresponding receiver coils 118 and configured to detect or measure RF signals emitted by the object in response to the excitation from the RF pulse.
 - the RF receiver 116 includes at least an amplifier and an analog-to- digital converter (ADC) to convert the received signal strength to a digital value to be processed by a computing device or other processing system.
 - ADC analog-to- digital converter
 - the gradient controller 108 can be configured to produce pulses or other gradient fields, using the gradient coils 110, along one or more directions relative to the orientation of the MRI apparatus or the object. By selection of such gradients and via other applied RF pulses using the transmitter coils 115, various imaging or MRI sequences can be applied.
 - An MRI sequence refers to a specific control scheme for generating pulses and or currents in the various coils to effectively acquire data to form at least one MRI image.
 - the object or specimen may be divided into volume elements (voxels) and MR signals for a plurality of gradient directions are acquired.
 - MR signals for a plurality of gradient directions are acquired.
 - components e.g., atoms, molecules
 - image data can be acquired by varying the gradient and RF pulses to measure a signal response corresponding to that slice and or location.
 - a processing device 124 such as a personal computer, a printed circuit board with one or more processors, memory, and/or other chipsets attached thereto, or the like, is provided for acquisition control and/or processing and analysis of collected data.
 - the processing device 124 may include volatile (e.g., dynamic random access memory (DRAM)) or non-volatile (e.g., hard disk drive (HDD), solid state drive (SSD), flash memory, or the like) memory.
 - the processing device 124 may also include a network interface controller (NIC) or other communications interface to communicate with other computing devices over, e.g., a network such as a local area network (LAN) or wide area network (WAN) such as the Internet.
 - NIC network interface controller
 - data can be transmitted from the processing device 124 to a service (e.g., a service hosted on a cloud platform) for remote processing or storage.
 - a service e.g., a service hosted on a cloud platform
 - Signal acquisition, instrument control, and signal analysis can be performed locally on the processing device 124, or in a distributed computing environment using a plurality of processing devices connected via a network.
 - signal acquisition and signal processing can be performed on different devices.
 - the processing device 124 can also be configured to implement an interface (e.g., a graphical user interface (GUI)) that enables a technician or operator of the MRI apparatus to control the MRI sequence scheme to select magnetic field magnitudes, gradient field direction and magnitude, RF pulse frequency, duration, and/or shape, and the like.
 - GUI graphical user interface
 - FIG. 2 illustrates a control system 200 for the MRI system 100 of FIG. 1, in accordance with some embodiments of the present disclosure.
 - the control system 200 is used to implement the functions of the hardware described above for the MRI system 100.
 - implementation may be in the form of a set of computer-readable instructions executed by a computing device.
 - any system that implements said functionality via any combination of hardware, firmware, or software is within the scope of the present disclosure.
 - the following description may be within the context of a single computing device, these functions can be implemented on a distributed system over multiple computing devices.
 - the control system 200 includes one or more processors 202, a memory 204, and a communications bus 206.
 - the communications bus 206 communicatively connects various system components, such as the processor(s) 202 and/or the memory 204.
 - the memory 204 may include any combination of read only memory (ROM), random access memory (RAM), or non-volatile memory (e.g., HDD, SSD, Flash, etc.).
 - the communications bus 206 may include any number of bus architecture types well-known in the art including but not limited to a peripheral component interconnect (PCI) bus, PCI Express (PCIe) bus, Universal Serial Bus (USB), or a dedicated chip-to-chip bus.
 - PCI peripheral component interconnect
 - PCIe PCI Express
 - USB Universal Serial Bus
 - the control system 200 may also include a number of input/output (VO) devices and/or interfaces such as a communication interface 212 (e.g., network interface controller (NIC)), keyboard/mouse 214, a video interface 216 (e.g., HDMI, etc.), audio interface 218 (e.g., speakers, microphone, etc.), or the like.
 - a communication interface 212 e.g., network interface controller (NIC)
 - keyboard/mouse 214 e.g., keyboard/mouse 214
 - video interface 216 e.g., HDMI, etc.
 - audio interface 218 e.g., speakers, microphone, etc.
 - the control system 200 may include software stored in non-volatile storage devices 220 (e.g., HDD, SSD, Flash, etc.), which may then be loaded into the memory 204.
 - the software may include an operating system and one or more applications.
 - the memory 204 may also store program modules and/or data that are loaded and/or used by at least one of the applications.
 - the program modules may include software add-ons that add customized functionality to an application.
 - the data may include fdes or other data structures storing information such as parameters that define an MRI sequence.
 - a GUI may be displayed on a display device 240.
 - the GUI may allow a technician or operator to enter commands or values to control the MRI apparatus 100, start or stop an MRI image acquisition sequence, view images generated in response to measured signals, and the like. Commands or values can be entered manually via the one or more I/O devices.
 - the display device 240 may include a touch-sensitive screen that allows a technician or operator to enter commands directly using the display device 240.
 - the control system 200 may include a communications interface 212 to enable the control system 200 to communicate with other remote devices, such as the MRI apparatus 100 or a server device 260, over a network.
 - the network interface may be wired (e.g., IEE 802.3 Ethernet) or wireless (e.g., 802.11 WiFi).
 - the network interface may be an optical interface (e.g., SONET).
 - the control system 200 may be directly connected to the MRI apparatus 100 via a wired interface.
 - Timing the acquisition of image data to match up with a particular phase of the subject’s respiratory cycle can help to improve image quality, although this can also extend the time required to capture all of the necessary image data.
 - a novel method for capturing the state of the person’s respiratory cycle during image acquisition while not unduly extending the total period required to acquire all of the necessary image data can help ease the comfort of the subject and save time to allow for more subjects to be imaged over the course of a day.
 - a navigator signal has been captured using, e.g., an ECG based on electrodes placed on the subject’s body or using bellows that capture the chest motion during inspiration and expiration.
 - ECG ECG
 - bellows that capture the chest motion during inspiration and expiration.
 - Another technique for generating a navigator signal uses the MRI system itself to capture images of a particular structure in the body (e.g., the person’s diaphragm, liver, etc.) in order to capture the person’s respiratory cycle due to motion of these structures. This typically involves capturing at least one image at a slice position corresponding to these structures and then tracking motion by analyzing a region of interest in the image.
 - these navigator signals were captured during normal image acquisition phase of the MRI sequence and may have added significant delays to capturing of the intended images.
 - Free-breathing image acquisition is an ideal solution for robustness and patient comfort.
 - a DC signal can be sampled each repetition time (TR) as a trajectory crosses k-space center, and this signal can be used as a free-breathing navigator signal to determine physiological motion.
 - TR repetition time
 - 3D “stack” trajectories the DC signal is measured less frequently, and may not have sufficient sampling frequency to resolve the respiratory cycle.
 - traditional trajectories especially 3D Cartesian, 3D stack-of-spirals, and 3D stack-of-stars trajectories
 - 3D Cartesian, 3D stack-of-spirals, and 3D stack-of-stars trajectories may need to be modified, thereby extending the required image acquisition time, in order to use a DC navigator signal.
 - a potentially “free” selfnavigator signal can be captured by measuring the signal immediately following an RF pulse for imaging acquisition.
 - This navigator signal can be implemented with minimal or no penalty in acquisition efficiency.
 - the navigator signal can be used for respiratory binning of a 3D stack-of-spiral image acquisition trajectory.
 - FIG. 3A shows a timing diagram for an MRI sequence, in accordance with some embodiments of the present disclosure.
 - the MRI sequence is implemented by generating a number of signals in a particular order in time.
 - the signals shown in FIG. 3A include an ADC signal that indicates when the receiver coil 118 is being sampled by the ADC in the RF receiver 116; an RF signal that indicates when the RF pulse is generated by the RF generator 114; and Gradient signals (e.g., Gx, Gy, Gz) that indicate the magnetic field gradients generated by the gradient controller 108 and the gradient coils 110.
 - Gradient signals e.g., Gx, Gy, Gz
 - 3A corresponds to a 3D stack-of-spirals acquisition trajectory within which the navigator signal (SSnav) is sampled during the slice-select ramp-down portion of the MRI sequence and before the image-encoding gradients begin.
 - Image encoding refers to the portion of the MRI sequence where spatial information is encoded by varying the current in one or more gradient coils to implement the k-space trajectory for image acquisition. This follows the slice select and/or slice rewinder portion of the MRI sequence.
 - the MRI sequence begins with slice select by generating a gradient signal in the Gz coil 110-3. Although the gradient shown here is axial along the z-axis, the gradient may be generated in any particular orientation using one or more of the gradient coils 110.
 - the RF pulse is generated in the RF generator 114.
 - slice select ramp down begins by ramping down the current in the Gz coil 110-3. It is during this stage of the MRI sequence that a “free” slice-select navigator (SSnav) signal can be sampled.
 - SSnav slice-select navigator
 - FIG. 3B shows the timing associated with three separate navigator signals in accordance with some embodiments of the present disclosure.
 - the navigator signals are referred to as: SSnav; PS Snav, and a MRI reference DCnav.
 - the S Snav signal is generated by sampling the ADC signal during a slice select ramp down portion of the MRI sequence before the current in the Gz coil is reduced to zero.
 - the PS Snav signal is generated while no current is passing through the Gz coil to reduce any signal variations due to the changing magnetic field gradient.
 - the S Snav and PSSnav signals may be referred to herein as “free” navigator signals because they can be collected in a time in the MRI sequence that requires no or minimal delay to the normal image acquisition timing of most traditional MRI sequences.
 - the slice-select zeroth gradient moment Moz (the gradient time integral) demonstrated in the detailed view in FIG. 3B shows different amounts of spatial encoding between the different navigator signals.
 - the proposed navigator signals, S Snav and PSSnav are collected before the image-encoding gradients are applied during normal MRI acquisition. There is no delay or minimal delay to image acquisition using the “free” navigator signals unlike significant delays added in other common methods for collecting navigator signals such as those that rely on image processing to identify the liver location or motion or an intermittent superiorinferior signal acquisition.
 - some MRI sequences sometimes referred to as speed-optimal imaging sequences using Cartesian trajectories, where image-encoding may begin during slice select ramp down, there may be a minimal penalty incurred due to delay of starting the imageencoding gradients until the navigator signal has been sampled.
 - 3B has been purposefully modified to require slice rewinder during each repetition time TR in order to generate a period for sampling the reference DCnav signal.
 - many MRI sequences do not include this slice rewinder each TR, and so this navigator signal may not be captured without modifying the current MRI sequences.
 - FIG. 3C illustrates an MRI sequence corresponding with a 2D Cartesian trajectory, in accordance with some embodiments of the present disclosure.
 - the SSnav signal is demonstrated within the timing of the 2D Cartesian MRI sequence acquisition.
 - the SSnav signal is placed during the slice-select ramp down portion of the MRI sequence, which requires that the imageencoding gradients (G x and G y ) be delayed until the start of the slice-select rewinding gradient after the navigator signal is acquired.
 - FIG. 3D shows various signal processing steps applied to a captured raw navigator signal, in accordance with some embodiments of the present disclosure.
 - the raw navigator signal can be filtered by an angular filter for k-space modulation, such as described in Di Sopra el al., “An automated approach to fully self-gated free-running cardiac and respiratory motion-resolved 5D whole- heart MRI,” J. Magnetic Resonance in Medicine, Vol. 82:6 (December 2019), which is herein incorporated by reference in its entirety.
 - temporal filtering can also be applied to remove any non-respiratory based temporal variations.
 - a band-pass filter can be applied to filter out variations in the signal that don’t appear to be cyclic at the same rate of a person’s respiration cycle.
 - any frequency response or temporal filter band may be selected that is appropriate for a given subject for common respiratory rates.
 - FIG. 3E illustrates the various navigator signals described herein, pre- and postprocessing, in accordance with some embodiments of the present disclosure.
 - the collected “free” navigator signals are compared to a conventional navigator signal generated based on a bellows placed on a patient’s chest.
 - the DCnav signal is also shown for reference. It can be seen that the “free” navigator signals appear to capture the patient’s respiratory cycle with a similar effectiveness as compared to the traditional bellows navigator signal and/or the DC nav signal.
 - FIG. 4 shows a set of short-axis slice images from a 3D stack-of- spirals ultra-short echo-time (UTE) acquisition, in accordance with some embodiments.
 - UTE ultra-short echo-time
 - the navigator signals described above may be incorporated into a number of common MRI sequences.
 - these navigator signals may be particularly suited for capturing a 2D trajectory or 3D stack trajectories associated with a slice selective pulse such as a sine pulse or UTE pulse.
 - Such MRI sequences are not particularly suited for conventional DCbased navigator signals as the trajectories in those pulse sequences do not return to the center of k-space at sufficient frequency to capture the respiratory cycle.
 - the SSnav and PSSnav signals can be captured without first returning to the center of k-space and, therefore, are particularly suited to 2D or 3D stack trajectories. As shown in FIG.
 - the slice images have been resolved for respiratory phase (highlighted by temporal line profiles depicting the lung-liver interfact in FIG. 4b, drawn from the dashed line in FIG. 4a), as determined by the SSnav, PSSnav and reference DCnav signal acquired every repetition time.
 - FIG. 5 illustrates a comparison between the four different navigator signals conducted during an experiment using ten subjects, in accordance with some embodiments.
 - the charts shown in FIG. 5 indicate the number of peaks detected (i.e., respiratory cycles) and a mean interval between peaks (i.e., average duration of the respiratory cycle), in seconds, over a 3 minute image acquisition interval for the four different subjects.
 - peaks detected i.e., respiratory cycles
 - mean interval between peaks i.e., average duration of the respiratory cycle
 - FIG. 6 is a flowchart of a method 600 for acquiring MRI images using navigator signals, in accordance with some embodiments of the present disclosure.
 - the method 600 can be implemented one or more processors of the control system 200. However, those of skill in the art will recognize that any system suitable of controlling an MRI apparatus is contemplated as being within the scope of the present disclosure.
 - the method 600 can be implemented as a set of instructions that are executed by one or more processors of a computing device that is communicatively connected to an MRI apparatus.
 - data for a navigator signal is captured using an MRI apparatus.
 - the navigator signal is captured by sampling (e.g., using an ADC) an RF receiver 116 connected to an RF coil 118.
 - the sample can include a number of data points captured every TR (sampling cycle) to capture MRI signal that includes some undesigned slice-encoding post slice-excitation, which can reduce the signal amplitude.
 - this signal will be proportional to the image volume being captured, the navigator signal over time will include motion within the imaging volume.
 - the frequent sampling every TR drives the filtering for the low signal navigator such that the physiological motion can be extracted from the filtering process.
 - the MRI apparatus may include multiple RF receiver coils 118 (which may each be coupled to a separate RF receiver 116 circuit) disposed at different locations proximate to the subject being imaged.
 - capturing the data for the navigator signal may further comprise selecting an optimal RF coil from the plurality of coils 118 to sample. Selecting the optimal RF coil can include comparing a power output of the signal received by each RF coil and selecting the RF coil corresponding with the highest power output as the optimal coil or selecting the coil that exhibits the most dominant motion within the imaging volume, such as described in Zhang et al., “Robust Self-Navigated Body MRI using Dense Coil Arrays,” J. Magnetic Resonance in Medicine, Vol. 76:1 (July 2016), which is incorporated by reference herein in its entirety.
 - the navigator signal is captured at a time in the MRI sequence between a slice select portion of the MRI sequence and a gradient refocus portion of the MRI sequence.
 - the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus.
 - RF radio frequency
 - the navigator signal can be captured after slice select ramp down is complete (i.e., after the current in the gradient coil reaches zero amps) but prior to a gradient refocus portion of the MRI sequence.
 - This period of the MRI sequence can be referred to as a free induction decay period, and the navigator captured at this time is referred to as the PSSnav signal.
 - the gradient refocus portion refers to a brief period after slice select encoding where the current in the gradient coil is reversed to return the magnetic field to a DC state.
 - the navigator signal is processed to identify timing of a respiratory cycle of a subject being measured.
 - the navigator signal can be filtered via an angular k- space filter.
 - the navigator signal can also be filtered by a temporal bandpass filter.
 - other types of filters can be applied, such as a low-pass or high-pass temporal filter.
 - data reconstruction is performed using the MRI apparatus based at least one the processed navigator signal.
 - Image signals are collected in a conventional manner after the navigator signal has been sampled.
 - the navigator signal can be utilized retroactively to generate MRI images for display or storage.
 - the navigator signal can be used to identify a phase of a respiratory cycle corresponding to the image data captured during each repetition time (TR).
 - the image data for a plurality of RF pulses over a plurality of TR periods can be binned according to the different phases of the respiratory cycle.
 - the image data for a particular bin i .e., corresponding to a particular phase of the respiratory cycle
 - a plurality of images are generated corresponding to a plurality of different phases of the respiratory cycle, thereby increasing the chances that motion artifacts are subdued or resolved.
 - the navigator signals may be used with a number of different types of trajectories, including, e.g., Cartesian, radial and spiral trajectories.
 - the navigator-based techniques described above are particularly suited for use with two-dimensional (2D) or three-dimensional (3D) imaging sequence using a stack trajectories in k-space.
 - the navigator-based technique may also be suitable for use with slice-selective RF pulse types such as sine or ultrashort echo time (UTE).
 
Landscapes
- Physics & Mathematics (AREA)
 - Health & Medical Sciences (AREA)
 - Engineering & Computer Science (AREA)
 - Cardiology (AREA)
 - Radiology & Medical Imaging (AREA)
 - Physiology (AREA)
 - Power Engineering (AREA)
 - Pulmonology (AREA)
 - General Health & Medical Sciences (AREA)
 - Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
 - Biophysics (AREA)
 - Signal Processing (AREA)
 - Life Sciences & Earth Sciences (AREA)
 - High Energy & Nuclear Physics (AREA)
 - Condensed Matter Physics & Semiconductors (AREA)
 - General Physics & Mathematics (AREA)
 - Magnetic Resonance Imaging Apparatus (AREA)
 
Abstract
A method and system for acquiring free-breathing magnetic resonance (MR) images are provided. The method utilizes a "free" navigator signal acquired by measuring MR signals during or following a slice select ramp down portion of an MRI pulse sequence. After processing the captured navigator signal with sampling-based and temporal filters, the resulting processed navigator signal may include information that correlates strongly with a respiratory cycle of a subject being imaged. Such navigator signal can be used retroactively to process the MRI image data to improve the image quality or resolve respiratory motion. This type of navigator signal can be used to allow fee-breathing of the subject during image acquisition, improving the comfort of the subject, while not significantly increasing the overall image capture time of many common MRI sequence.
  Description
 PHYSIOLOGICAL NAVIGATION IN MRI USING INHERENT SLAB -SELECTIVE FREE-
    INDUCTION-DECAY SIGNAL 
    CROSS-REFERENCE TO RELATED APPLICATIONS 
    [0001] This application claims the benefit of U.S. Provisional Application No. 63/411,055, filed September 28, 2022, which is herein incorporated by reference in its entirety. 
    STATEMENT REGARDING
    FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
    [0002] This invention was made with Government support under Grant Numbers Z01- HL006257 and Z01-HL006213 awarded by the National Institutes of Health. The Government has certain rights in this invention. 
    BACKGROUND 
    [0003] In magnetic resonance imaging (MRI), physiological "navigation" is required for high-quality imaging which accounts for respiratory motion (and other motion), to avoid image artifacts and blurring. "Free-breathing" acquisitions are ideal for patient comfort and are more robust when patients are less compliant to holding their breath (i.e. children, or patients under anesthesia). 
    [0004] Recently, techniques for synchronizing image acquisition with timing of the patient’s respiratory cycle using a navigator signal have been proposed. The navigation data can be acquired using the MRI system separate and distinct from capturing the image data. However, these techniques typically capture the navigation data in lieu of image data during normal acquisition timing cycles, which can increase the overall time to acquire the necessary image data. Typically, one or more MRI sequence cycles are used to acquire the navigator signal and, when the navigator signal indicates that no respiratory motion is present, the MRI sequence cycles capturing the image data of the structure under observation commence. Reducing the timing required to capture the necessary MRI image data can improve the comfort of a patient, 
who does not have to be confined in an MRI machine for as long. Thus, a technique for efficiently capturing a navigator signal during the same MRI sequence cycle as the image data is desired. 
    SUMMARY 
    [0005] A navigator signal may be captured during a portion of the MRI pulse sequence after slice (or slab) selection and prior to k-space encoding for capturing image data using well-known MRI pulse sequence trajectories. The navigator signal enables a subject being imaged to continue breathing normally during image acquisition as motion compensation can be performed retroactively using post-processing techniques (e.g., binning) based on the captured navigator signal. 
    [0006] In accordance with a first aspect of the present disclosure, a method applied using a magnetic resonance imaging (MRI) apparatus is provided. The method includes: capturing data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, processing the navigator signal to identify timing of a respiratory cycle of a subject being measured, and performing data reconstruction using the MRI apparatus based at least on the processed navigator signal. Performing data reconstruction can include performing retrospective MRI data binning based on the processed navigator signal. The slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus. 
    [0007] In accordance with an embodiment of the first aspect, the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence. The slice select ramp down portion of the MRI sequence is defined as a period immediately subsequent to the RF pulse generation during which current in the at least one gradient coil is reduced to zero. 
    [0008] In accordance with an embodiment of the first aspect, the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through the at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition. 
[0009] In accordance with an embodiment of the first aspect, processing the navigator signal comprises filtering the navigator signal via an angular k-space filter. 
    [0010] In accordance with an embodiment of the first aspect, the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to a common respiration rate. 
    [0011] In accordance with an embodiment of the first aspect, performing data acquisition using the MRI apparatus comprises gating image data according to a phase of the respiratory cycle identified using the processed navigator signal. 
    [0012] In accordance with an embodiment of the first aspect, the MRI sequence comprises a two dimensional (2D) or three-dimensional (3D) imaging sequence using a slice-selective excitation. 
    [0013] In accordance with an embodiment of the first aspect, the RF pulse comprises a slice- selective RF pulse. The slice-selective RF pulse can include a sine or a ultra-short echo time (UTE) pulse. 
    [0014] In accordance with an embodiment of the first aspect, the MRI apparatus includes a plurality of RF receivers coupled to a plurality of RF coils. In addition, capturing the data for the navigator signal further comprises selecting an optimal RF coil from the plurality of coils to sample. 
    [0015] In accordance with a second aspect of the present disclosure, a magnetic resonance imaging (MRI) system is provided. The MRI system includes: an MRI apparatus; and a control system connected to the MRI apparatus. The control system includes at least one processor and a memory storing a set of instructions that, in response to being executed by the at least one processor, cause the MRI apparatus to: capture data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, process the navigator signal to identify timing of a respiratory cycle of a subject being measured, and perform data reconstruction using the MRI apparatus based at least on the processed navigator signal. Performing data reconstruction can include performing data sorting and image formation based at least on the processed navigator signal. The slice select portion of the MRI sequence is a 
period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus 
    [0016] In accordance with an embodiment of the second aspect, the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence. 
    [0017] In accordance with an embodiment of the second aspect, wherein the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition. 
    [0018] In accordance with an embodiment of the second aspect, processing the navigator signal comprises filtering the navigator signal via an angular k-space filter. 
    [0019] In accordance with an embodiment of the second aspect, the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to common respiration rates. 
    [0020] In accordance with an embodiment of the second aspect, performing data acquisition using the MRI apparatus comprises binning image data according to a phase of the respiratory cycle identified using the processed navigator signal. 
    [0021] In accordance with an embodiment of the second aspect, performing data acquisition using the MRI apparatus comprises gating image data according to a phase of the respiratory cycle identified using the processed navigator signal. 
    [0022] In accordance with an embodiment of the second aspect, the MRI sequence comprises a two-dimensional (2D) or three-dimensional (3D) imaging sequence using slice-selective excitation. 
    [0023] In accordance with an embodiment of the second aspect, the RF pulse comprises a slice-selective RF pulse. 
    [0024] In accordance with a third aspect of the present disclosure, a non-transitory computer- readable storage medium storing instructions is provided that, in response to being executed by one or more processors, cause a magnetic resonance imaging (MRI) apparatus to perform the method of the first aspect. 
 BRIEF DESCRIPTION OF THE FIGURES 
    [0025] FIG. 1 illustrates a magnetic resonance imaging (MRI) system, in accordance with some embodiments of the present disclosure. 
    [0026] FIG. 2 illustrates a control system for the MRI system of FIG. 1, in accordance with some embodiments of the present disclosure. 
    [0027] FIG. 3A shows a timing diagram for an MRI sequence, in accordance with some embodiments of the present disclosure. 
    [0028] FIG. 3B shows the timing associated with three separate navigator signals, in accordance with some embodiments of the present disclosure. 
    [0029] FIG. 3C illustrates an MRI sequence corresponding with a 2D Cartesian trajectory, in accordance with some embodiments of the present disclosure. 
    [0030] FIG. 3D shows various signal processing steps applied to a captured raw navigator signal, in accordance with some embodiments of the present disclosure. 
    [0031] FIG. 3E illustrates the various navigator signals described herein, pre- and postprocessing, in accordance with some embodiments of the present disclosure. 
    [0032] FIG. 4 shows a set of short-axis slice images from a 3D stack-of-spirals ultra-short echo-time acquisition, in accordance with some embodiments. 
    [0033] FIG. 5 illustrates a comparison between the four different navigator signals conducted during an experiment using ten subjects, in accordance with some embodiments. 
    [0034] FIG. 6 is a flowchart of a method for acquiring MRI images using navigator signals, in accordance with some embodiments of the present disclosure. 
    DETAILED DESCRIPTION 
    [0035] Physiological information is used for MRI to guide acquisition timing for motion- resolved or motion-compensated imaging. Existing means for generating signals to track a patient’s respiration cycle include use of devices on patients, such as electrodes for ECG (electrocardiogram) or pressure-sensitive respiratory bellows. Hardware modifications are also possible for radio frequency (RF) sensing, but this in not universally available. 
    [0036] Alternatively, modem advanced techniques for acquisition of free-breathing MRI images have two means to detect physiological motion in the data collection itself: 1) use of a 
navigation-only signal acquisition to determine the liver position, and 2) designing the image acquisition such that there is frequent sampling of a "direct-current" signal, at a timing where there is no spatial encoding and high MRI signal, that fluctuates with respiration. However, navigation segments can disrupt image acquisitions that rely on steady-state or fast imaging regimes, and implementing frequent sampling of "direct-current" signals can be difficult as well as not being possible to implement for every application. 
    [0037] The present disclosure describes acquisition of navigator signals by frequently sampling the signal emitted after the RF pulse and the slice select gradient refocus portion of the sequence, which has some unspecified spatial encoding - the "unintended" encoding present following a slice-selective pulse - used in the majority of MRI imaging techniques for choosing volumes of interest. As this unwanted signal is acquired during image acquisition "dead-time", it is therefore an opportunity to acquire navigation information without disrupting the timing of the MRI sequence. In other words, while other existing techniques may require delay of the image acquisition during the MRI sequence, the technique described below can be implemented in parallel with image acquisition without disrupting the image data. It has been determined that there is physiological data present in the dead-time signal, and processing to use this signal for respiratory gating has been proven to be effective. As this signal can be acquired during the image data acquisition dead-time, this can be added to a number of MRI imaging sequences with minimal impact. 
    [0038] MRI applications used for determining respiratory motion, cardiac motion, and bulk subject motion, applied to cardiac imaging, pulmonary imaging, abdominal imaging, neuroimaging, musculoskeletal imaging are all within the scope of the present disclosure. 
    [0039] FIG. 1 illustrates an MRI system 100, in accordance with some embodiments of the present disclosure. Magnetic resonance (MR) measurements can be obtained using an MRI apparatus 100 as illustrated in FIG. 1. The MRI apparatus 100 includes a controller/interface 102 that is configured to apply selected magnetic fields such as constant or pulsed field gradients to an object (e.g., a patient) located proximate to the MRI apparatus 100. An axial controller 104 is in communication with an axial magnet 106 that is generally configured to produce a substantially constant magnetic field Bo proximate the MRI apparatus 100. The axial magnet 106 may comprise a coil that includes superconducting material such as niobium-tin (Nb3Sn) or 
Magnesium Diboride (MgB2). A gradient controller 108 is configure to apply a constant or timevarying magnetic field gradient in a selected direction using a number of gradient coils 110. The gradient controller 108 may cause electrical current to be applied to the gradient coils 110 to apply diffusion or slice selection gradients to the object. The magnetic field gradient vector components may be referred to as Gx, Gy, and Gz corresponding to an x-axis, y-axis, and z-axis, respectively. In an embodiment, each axis is associated with a separate gradient coil, e.g., a first gradient coil 110-1 is associated with an x-axis, a second gradient coil 110-2 is associated with a y-axis, and a third gradient coil 110-3 is associated with a z-axis. A magnetic field gradient may be oriented in any direction, meaning the magnitude of the magnetic field changes along an axis corresponding to that direction, such that two or more magnetic field gradient vector components are non-zero. 
    [0040] An RF generator 114 is configured to generate one or more RF pulses of the electromagnetic field to the object being imaged. The RF generator 114 is connected to one or more transmitter coils 115, and delivers an electrical current to the transmitter coils 115 to generate the RF pulse. One or more RF receivers 116 are coupled to corresponding receiver coils 118 and configured to detect or measure RF signals emitted by the object in response to the excitation from the RF pulse. Typically, the RF receiver 116 includes at least an amplifier and an analog-to- digital converter (ADC) to convert the received signal strength to a digital value to be processed by a computing device or other processing system. 
    [0041] The gradient controller 108 can be configured to produce pulses or other gradient fields, using the gradient coils 110, along one or more directions relative to the orientation of the MRI apparatus or the object. By selection of such gradients and via other applied RF pulses using the transmitter coils 115, various imaging or MRI sequences can be applied. An MRI sequence refers to a specific control scheme for generating pulses and or currents in the various coils to effectively acquire data to form at least one MRI image. 
    [0042] During imaging, the object or specimen may be divided into volume elements (voxels) and MR signals for a plurality of gradient directions are acquired. It will be appreciated that by selecting a magnetic field gradient along a specific direction and then exciting the components (e.g., atoms, molecules) of the object at a specific frequency of RF pulse, certain components of the object in a plane normal to the direction of the gradient may be excited more 
strongly than other components of the object further away from the plane. Thus, image data can be acquired by varying the gradient and RF pulses to measure a signal response corresponding to that slice and or location. 
    [0043] A processing device 124 such as a personal computer, a printed circuit board with one or more processors, memory, and/or other chipsets attached thereto, or the like, is provided for acquisition control and/or processing and analysis of collected data. The processing device 124 may include volatile (e.g., dynamic random access memory (DRAM)) or non-volatile (e.g., hard disk drive (HDD), solid state drive (SSD), flash memory, or the like) memory. The processing device 124 may also include a network interface controller (NIC) or other communications interface to communicate with other computing devices over, e.g., a network such as a local area network (LAN) or wide area network (WAN) such as the Internet. In some embodiments, data can be transmitted from the processing device 124 to a service (e.g., a service hosted on a cloud platform) for remote processing or storage. Signal acquisition, instrument control, and signal analysis can be performed locally on the processing device 124, or in a distributed computing environment using a plurality of processing devices connected via a network. For example, in some embodiments, signal acquisition and signal processing can be performed on different devices. The processing device 124 can also be configured to implement an interface (e.g., a graphical user interface (GUI)) that enables a technician or operator of the MRI apparatus to control the MRI sequence scheme to select magnetic field magnitudes, gradient field direction and magnitude, RF pulse frequency, duration, and/or shape, and the like. The processing device 124 can implement the interface and any other functions in computer-executable instructions capable of being executed by the one or more processors of the processing device 124, the computer-executable instructions stored in a memory, such as the HDD and loaded in the DRAM for execution. In addition, parameters for the various MRI sequence schema can be loaded into the memory as well. In general, control and data acquisition of the MRI apparatus 100 can be controlled locally or remotely via instructions and/or data transmitted remotely to the processing device 124 via a network connection. 
    [0044] FIG. 2 illustrates a control system 200 for the MRI system 100 of FIG. 1, in accordance with some embodiments of the present disclosure. The control system 200 is used to implement the functions of the hardware described above for the MRI system 100. In some 
embodiments, implementation may be in the form of a set of computer-readable instructions executed by a computing device. However, any system that implements said functionality via any combination of hardware, firmware, or software is within the scope of the present disclosure. In addition, although the following description may be within the context of a single computing device, these functions can be implemented on a distributed system over multiple computing devices. 
    [0045] As shown in FIG. 2, the control system 200 includes one or more processors 202, a memory 204, and a communications bus 206. The communications bus 206 communicatively connects various system components, such as the processor(s) 202 and/or the memory 204. The memory 204 may include any combination of read only memory (ROM), random access memory (RAM), or non-volatile memory (e.g., HDD, SSD, Flash, etc.). The communications bus 206 may include any number of bus architecture types well-known in the art including but not limited to a peripheral component interconnect (PCI) bus, PCI Express (PCIe) bus, Universal Serial Bus (USB), or a dedicated chip-to-chip bus. The control system 200 may also include a number of input/output (VO) devices and/or interfaces such as a communication interface 212 (e.g., network interface controller (NIC)), keyboard/mouse 214, a video interface 216 (e.g., HDMI, etc.), audio interface 218 (e.g., speakers, microphone, etc.), or the like. 
    [0046] The control system 200 may include software stored in non-volatile storage devices 220 (e.g., HDD, SSD, Flash, etc.), which may then be loaded into the memory 204. The software may include an operating system and one or more applications. 
    [0047] In some embodiments, the memory 204 may also store program modules and/or data that are loaded and/or used by at least one of the applications. The program modules may include software add-ons that add customized functionality to an application. The data may include fdes or other data structures storing information such as parameters that define an MRI sequence. 
    [0048] A GUI may be displayed on a display device 240. The GUI may allow a technician or operator to enter commands or values to control the MRI apparatus 100, start or stop an MRI image acquisition sequence, view images generated in response to measured signals, and the like. Commands or values can be entered manually via the one or more I/O devices. In some embodiments, the display device 240 may include a touch-sensitive screen that allows a technician or operator to enter commands directly using the display device 240. 
[0049] Again, the control system 200 may include a communications interface 212 to enable the control system 200 to communicate with other remote devices, such as the MRI apparatus 100 or a server device 260, over a network. The network interface may be wired (e.g., IEE 802.3 Ethernet) or wireless (e.g., 802.11 WiFi). In some embodiments, the network interface may be an optical interface (e.g., SONET). In other embodiments, the control system 200 may be directly connected to the MRI apparatus 100 via a wired interface. 
    MRI Image Acquisition using Navigator Signals 
    [0050] When capturing images using the MRI apparatus 100, it may be necessary to keep the patient or subject being imaged relatively still in order to get better quality images. However, while a person breathes, there will necessarily be motion of the subject’s body and internal organs as the diaphragm moves and the lungs inhale or exhale air. Traditionally, images with reduced motion artifacts were acquired by requesting the patient to hold their breath while the image data was acquired. Newer techniques use navigator signals that detect motion of the patient due to respiration and time image acquisition or inform image post-processing accordingly, allowing the subject to continue breathing normally while image data is acquired. Timing the acquisition of image data to match up with a particular phase of the subject’s respiratory cycle can help to improve image quality, although this can also extend the time required to capture all of the necessary image data. A novel method for capturing the state of the person’s respiratory cycle during image acquisition while not unduly extending the total period required to acquire all of the necessary image data can help ease the comfort of the subject and save time to allow for more subjects to be imaged over the course of a day. 
    [0051] Conventionally, a navigator signal has been captured using, e.g., an ECG based on electrodes placed on the subject’s body or using bellows that capture the chest motion during inspiration and expiration. However, these techniques are not as comfortable to the patient and take additional steps to set up by a technician. Another technique for generating a navigator signal uses the MRI system itself to capture images of a particular structure in the body (e.g., the person’s diaphragm, liver, etc.) in order to capture the person’s respiratory cycle due to motion of these structures. This typically involves capturing at least one image at a slice position corresponding to these structures and then tracking motion by analyzing a region of interest in the image. However, these navigator signals were captured during normal image acquisition 
phase of the MRI sequence and may have added significant delays to capturing of the intended images. 
    [0052] Free-breathing image acquisition is an ideal solution for robustness and patient comfort. A DC signal can be sampled each repetition time (TR) as a trajectory crosses k-space center, and this signal can be used as a free-breathing navigator signal to determine physiological motion. However, with 2D or 3D “stack” trajectories, the DC signal is measured less frequently, and may not have sufficient sampling frequency to resolve the respiratory cycle. As such, traditional trajectories (especially 3D Cartesian, 3D stack-of-spirals, and 3D stack-of-stars trajectories) may need to be modified, thereby extending the required image acquisition time, in order to use a DC navigator signal. In contrast, as proposed herein, a potentially “free” selfnavigator signal can be captured by measuring the signal immediately following an RF pulse for imaging acquisition. This navigator signal can be implemented with minimal or no penalty in acquisition efficiency. In an example application, the navigator signal can be used for respiratory binning of a 3D stack-of-spiral image acquisition trajectory. 
    [0053] FIG. 3A shows a timing diagram for an MRI sequence, in accordance with some embodiments of the present disclosure. The MRI sequence is implemented by generating a number of signals in a particular order in time. The signals shown in FIG. 3A include an ADC signal that indicates when the receiver coil 118 is being sampled by the ADC in the RF receiver 116; an RF signal that indicates when the RF pulse is generated by the RF generator 114; and Gradient signals (e.g., Gx, Gy, Gz) that indicate the magnetic field gradients generated by the gradient controller 108 and the gradient coils 110. The MRI sequence shown in FIG. 3A corresponds to a 3D stack-of-spirals acquisition trajectory within which the navigator signal (SSnav) is sampled during the slice-select ramp-down portion of the MRI sequence and before the image-encoding gradients begin. Image encoding refers to the portion of the MRI sequence where spatial information is encoded by varying the current in one or more gradient coils to implement the k-space trajectory for image acquisition. This follows the slice select and/or slice rewinder portion of the MRI sequence. 
    [0054] The MRI sequence begins with slice select by generating a gradient signal in the Gz coil 110-3. Although the gradient shown here is axial along the z-axis, the gradient may be generated in any particular orientation using one or more of the gradient coils 110. Once the 
magnetic field gradient has been set, the RF pulse is generated in the RF generator 114. After the RF pulse is complete, slice select ramp down begins by ramping down the current in the Gz coil 110-3. It is during this stage of the MRI sequence that a “free” slice-select navigator (SSnav) signal can be sampled. 
    [0055] FIG. 3B shows the timing associated with three separate navigator signals in accordance with some embodiments of the present disclosure. The navigator signals are referred to as: SSnav; PS Snav, and a MRI reference DCnav. The S Snav signal is generated by sampling the ADC signal during a slice select ramp down portion of the MRI sequence before the current in the Gz coil is reduced to zero. The PS Snav signal is generated while no current is passing through the Gz coil to reduce any signal variations due to the changing magnetic field gradient. The S Snav and PSSnav signals may be referred to herein as “free” navigator signals because they can be collected in a time in the MRI sequence that requires no or minimal delay to the normal image acquisition timing of most traditional MRI sequences. The slice-select zeroth gradient moment Moz (the gradient time integral) demonstrated in the detailed view in FIG. 3B shows different amounts of spatial encoding between the different navigator signals. The DCnav signal is generated after the current in the Gz coil is briefly reversed to refocus the zeroth gradient moment at Moz = 0. While the DCnav signal is shown here for comparison to the “free” navigator signals, it may not be an ideal navigator signal because some MRI sequences do not refocus the magnetic field at sufficient frequency to be able to capture enough data to use as a navigator signal. 
    [0056] It will be appreciated that the proposed navigator signals, S Snav and PSSnav are collected before the image-encoding gradients are applied during normal MRI acquisition. There is no delay or minimal delay to image acquisition using the “free” navigator signals unlike significant delays added in other common methods for collecting navigator signals such as those that rely on image processing to identify the liver location or motion or an intermittent superiorinferior signal acquisition. In some MRI sequences, sometimes referred to as speed-optimal imaging sequences using Cartesian trajectories, where image-encoding may begin during slice select ramp down, there may be a minimal penalty incurred due to delay of starting the imageencoding gradients until the navigator signal has been sampled. 
[0057] It will be appreciated that the MRI sequence shown in FIG. 3B has been purposefully modified to require slice rewinder during each repetition time TR in order to generate a period for sampling the reference DCnav signal. In practice, many MRI sequences do not include this slice rewinder each TR, and so this navigator signal may not be captured without modifying the current MRI sequences. 
    [0058] FIG. 3C illustrates an MRI sequence corresponding with a 2D Cartesian trajectory, in accordance with some embodiments of the present disclosure. The SSnav signal is demonstrated within the timing of the 2D Cartesian MRI sequence acquisition. The SSnav signal is placed during the slice-select ramp down portion of the MRI sequence, which requires that the imageencoding gradients (Gx and Gy) be delayed until the start of the slice-select rewinding gradient after the navigator signal is acquired. 
    [0059] FIG. 3D shows various signal processing steps applied to a captured raw navigator signal, in accordance with some embodiments of the present disclosure. As shown in FIG. 3D, the raw navigator signal can be filtered by an angular filter for k-space modulation, such as described in Di Sopra el al., “An automated approach to fully self-gated free-running cardiac and respiratory motion-resolved 5D whole- heart MRI,” J. Magnetic Resonance in Medicine, Vol. 82:6 (December 2019), which is herein incorporated by reference in its entirety. In addition, temporal filtering can also be applied to remove any non-respiratory based temporal variations. For example, a band-pass filter can be applied to filter out variations in the signal that don’t appear to be cyclic at the same rate of a person’s respiration cycle. Alternatively, any frequency response or temporal filter band may be selected that is appropriate for a given subject for common respiratory rates. 
    [0060] FIG. 3E illustrates the various navigator signals described herein, pre- and postprocessing, in accordance with some embodiments of the present disclosure. As shown in FIG. 3E, the collected “free” navigator signals are compared to a conventional navigator signal generated based on a bellows placed on a patient’s chest. The DCnav signal is also shown for reference. It can be seen that the “free” navigator signals appear to capture the patient’s respiratory cycle with a similar effectiveness as compared to the traditional bellows navigator signal and/or the DC nav signal. 
[0061] FIG. 4 shows a set of short-axis slice images from a 3D stack-of- spirals ultra-short echo-time (UTE) acquisition, in accordance with some embodiments. It will be appreciated that the navigator signals described above (e.g., SSnav or PSSnav) may be incorporated into a number of common MRI sequences. In particular, these navigator signals may be particularly suited for capturing a 2D trajectory or 3D stack trajectories associated with a slice selective pulse such as a sine pulse or UTE pulse. Such MRI sequences are not particularly suited for conventional DCbased navigator signals as the trajectories in those pulse sequences do not return to the center of k-space at sufficient frequency to capture the respiratory cycle. In contrast, the SSnav and PSSnav signals can be captured without first returning to the center of k-space and, therefore, are particularly suited to 2D or 3D stack trajectories. As shown in FIG. 4, the slice images have been resolved for respiratory phase (highlighted by temporal line profiles depicting the lung-liver interfact in FIG. 4b, drawn from the dashed line in FIG. 4a), as determined by the SSnav, PSSnav and reference DCnav signal acquired every repetition time. 
    [0062] FIG. 5 illustrates a comparison between the four different navigator signals conducted during an experiment using ten subjects, in accordance with some embodiments. The charts shown in FIG. 5 indicate the number of peaks detected (i.e., respiratory cycles) and a mean interval between peaks (i.e., average duration of the respiratory cycle), in seconds, over a 3 minute image acquisition interval for the four different subjects. When comparing the different navigator signals, no significant difference was found between the number of peaks detected (p > 0.43, ANOVA), and no significant differences were found between mean respiratory interval (p > 0.31, ANO V A). 
    [0063] FIG. 6 is a flowchart of a method 600 for acquiring MRI images using navigator signals, in accordance with some embodiments of the present disclosure. The method 600 can be implemented one or more processors of the control system 200. However, those of skill in the art will recognize that any system suitable of controlling an MRI apparatus is contemplated as being within the scope of the present disclosure. In an embodiment, the method 600 can be implemented as a set of instructions that are executed by one or more processors of a computing device that is communicatively connected to an MRI apparatus. 
    [0064] At 602, data for a navigator signal is captured using an MRI apparatus. In an embodiment, the navigator signal is captured by sampling (e.g., using an ADC) an RF receiver 
116 connected to an RF coil 118. The sample can include a number of data points captured every TR (sampling cycle) to capture MRI signal that includes some undesigned slice-encoding post slice-excitation, which can reduce the signal amplitude. As this signal will be proportional to the image volume being captured, the navigator signal over time will include motion within the imaging volume. The frequent sampling every TR drives the filtering for the low signal navigator such that the physiological motion can be extracted from the filtering process. 
    [0065] In an embodiment, the MRI apparatus may include multiple RF receiver coils 118 (which may each be coupled to a separate RF receiver 116 circuit) disposed at different locations proximate to the subject being imaged. In such cases, capturing the data for the navigator signal may further comprise selecting an optimal RF coil from the plurality of coils 118 to sample. Selecting the optimal RF coil can include comparing a power output of the signal received by each RF coil and selecting the RF coil corresponding with the highest power output as the optimal coil or selecting the coil that exhibits the most dominant motion within the imaging volume, such as described in Zhang et al., “Robust Self-Navigated Body MRI using Dense Coil Arrays,” J. Magnetic Resonance in Medicine, Vol. 76:1 (July 2016), which is incorporated by reference herein in its entirety. 
    [0066] In an embodiment, the navigator signal is captured at a time in the MRI sequence between a slice select portion of the MRI sequence and a gradient refocus portion of the MRI sequence. The slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus. Once the RF pulse is complete, and as the current in the gradient coil is ramping down to zero current, a SSnav signal can be captured by the MRI apparatus. Alternatively, in another embodiment, the navigator signal can be captured after slice select ramp down is complete (i.e., after the current in the gradient coil reaches zero amps) but prior to a gradient refocus portion of the MRI sequence. This period of the MRI sequence can be referred to as a free induction decay period, and the navigator captured at this time is referred to as the PSSnav signal. The gradient refocus portion refers to a brief period after slice select encoding where the current in the gradient coil is reversed to return the magnetic field to a DC state. 
[0067] At 604, the navigator signal is processed to identify timing of a respiratory cycle of a subject being measured. In an embodiment, the navigator signal can be filtered via an angular k- space filter. In some embodiments, the navigator signal can also be filtered by a temporal bandpass filter. Alternatively, other types of filters can be applied, such as a low-pass or high-pass temporal filter. 
    [0068] At 606, data reconstruction is performed using the MRI apparatus based at least one the processed navigator signal. Image signals are collected in a conventional manner after the navigator signal has been sampled. In an embodiment, the navigator signal can be utilized retroactively to generate MRI images for display or storage. The navigator signal can be used to identify a phase of a respiratory cycle corresponding to the image data captured during each repetition time (TR). In one case, the image data for a plurality of RF pulses over a plurality of TR periods can be binned according to the different phases of the respiratory cycle. The image data for a particular bin (i .e., corresponding to a particular phase of the respiratory cycle) is then used to reconstruct an image or images corresponding to that particular phase of the respiratory cycle. As such, a plurality of images are generated corresponding to a plurality of different phases of the respiratory cycle, thereby increasing the chances that motion artifacts are subdued or resolved. 
    [0069] It will be appreciated that the steps of the method 600 maybe repeated for each TR of an MRI sequence in order to capture data points for the desired k-space trajectories. The navigator signals may be used with a number of different types of trajectories, including, e.g., Cartesian, radial and spiral trajectories. In particular, the navigator-based techniques described above are particularly suited for use with two-dimensional (2D) or three-dimensional (3D) imaging sequence using a stack trajectories in k-space. Furthermore, the navigator-based technique may also be suitable for use with slice-selective RF pulse types such as sine or ultrashort echo time (UTE). 
    [0070] All references, including publications, patent applications, and patents, cited herein are hereby incorporated by reference to the same extent as if each reference were individually and specifically indicated to be incorporated by reference and were set forth in its entirety herein. 
[0071] The use of the terms “a” and “an” and “the” and “at least one” and similar referents in the context of describing the invention (especially in the context of the following claims) are to be construed to cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted by context. The use of the term “at least one” followed by a list of one or more items (for example, “at least one of A and B”) is to be construed to mean one item selected from the listed items (A or B) or any combination of two or more of the listed items (A and B), unless otherwise indicated herein or clearly contradicted by context. The terms “comprising,” “having,” “including,” and “containing” are to be construed as open-ended terms (i.e., meaning “including, but not limited to,”) unless otherwise noted. Recitation of ranges of values herein are merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or exemplary language (e.g., “such as”) provided herein, is intended merely to better illuminate the invention and does not pose a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention. 
    [0072] Preferred embodiments of this invention are described herein, including the best mode known to the inventors for carrying out the invention. Variations of those preferred embodiments may become apparent to those of ordinary skill in the art upon reading the foregoing description. The inventors expect skilled artisans to employ such variations as appropriate, and the inventors intend for the invention to be practiced otherwise than as specifically described herein. Accordingly, this invention includes all modifications and equivalents of the subject matter recited in the claims appended hereto as permitted by applicable law. Moreover, any combination of the above-described elements in all possible variations thereof is encompassed by the invention unless otherwise indicated herein or otherwise clearly contradicted by context. 
  Claims
1. A method applied using a magnetic resonance imaging (MRI) apparatus, the method comprising: capturing data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, wherein the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus; processing the navigator signal to identify timing of a respiratory cycle of a subject being measured; and performing data reconstruction using the MRI apparatus based at least on the processed navigator signal. 
    2. The method of claim 1, wherein the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence. 
    3. The method of claim 1, wherein the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through the at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition. 
    4. The method of claim 1, wherein processing the navigator signal comprises filtering the navigator signal via an angular k-space filter. 
    5. The method of claim 4, wherein the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to a common respiration rate. 
    6. The method of claim 4, wherein the MRI apparatus includes a plurality of RF receivers coupled to a plurality of RF coils, and wherein capturing the data for the navigator signal further comprises selecting an optimal RF coil from the plurality of coils to sample. 
    7. The method of claim 1, wherein performing data acquisition using the MRI apparatus comprises binning image data according to a phase of the respiratory cycle identified using the processed navigator signal. 
    8. The method of claim 1, wherein the MRI sequence comprises a two-dimensional (2D) or three-dimensional (3D) imaging sequence using a stack trajectory in k-space. 
    9. The method of claim 1, wherein the RF pulse comprises a slice-selective pulse. 
    10. A magnetic resonance imaging (MRI) system, the MRI system comprising: an MRI apparatus; and a control system connected to the MRI apparatus, wherein the control system includes at least one processor and a memory storing a set of instructions that, in response to being executed by the at least one processor, cause the MRI apparatus to: capture data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, wherein the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus; process the navigator signal to identify timing of a respiratory cycle of a subject being measured; and perform data reconstruction using the MRI apparatus based at least on the processed navigator signal. 
    11. The MRI system of claim 10, wherein the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence. 
    12. The MRI system of claim 10, wherein the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through 
the at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of image-encoding gradients for data acquisition. 
    13. The MRI system of claim 10, wherein processing the navigator signal comprises filtering the navigator signal via an angular k-space filter. 
    14. The MRI system of claim 13, wherein the filtering further comprises applying a temporal band-pass filter centered on a frequency corresponding to a common respiration rate. 
    15. The MRI system of claim 10, wherein performing data reconstruction using the MRI apparatus comprises binning image data according to a phase of the respiratory cycle identified using the processed navigator signal. 
    16. The MRI system of claim 10, wherein the MRI sequence comprises a two- dimensional (2D) or three-dimensional (3D) imaging sequence using a stack trajectory in k- space. 
    17. The MRI system of claim 10, wherein the RF pulse comprises of a slice-selective pulse. 
    18. A non-transitory computer-readable storage medium storing instructions that, in response to being executed by one or more processors, cause a magnetic resonance imaging (MRI) apparatus to: capture data for a navigator signal between a slice select portion of an MRI sequence and a gradient refocus portion of the MRI sequence, wherein the slice select portion of the MRI sequence is a period where a current is passed through at least one gradient coil of the MRI apparatus while a radio frequency (RF) pulse is generated by a RF coil of the MRI apparatus; process the navigator signal to identify timing of a respiratory cycle of a subject being measured; and 
 perform data reconstruction using the MRI apparatus based at least on the processed navigator signal. 
    19. The non-transitory computer-readable storage medium of claim 18, wherein the data for the navigator signal is captured during a slice select ramp down portion of the MRI sequence. 
    20. The non-transitory computer-readable storage medium of claim 18, wherein the data for the navigator signal is captured after a slice select ramp down portion of the MRI sequence, when no current is passed through the at least one gradient coil, and before both a gradient refocus portion of the MRI sequence and a start of generation of image-encoding gradients for data acquisition. 
    Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title | 
|---|---|---|---|
| US202263411055P | 2022-09-28 | 2022-09-28 | |
| PCT/US2023/033997 WO2024072976A1 (en) | 2022-09-28 | 2023-09-28 | Physiological navigation in mri using inherent slab-selective free-induction-decay signal | 
Publications (1)
| Publication Number | Publication Date | 
|---|---|
| EP4594767A1 true EP4594767A1 (en) | 2025-08-06 | 
Family
ID=88506900
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date | 
|---|---|---|---|
| EP23793572.1A Pending EP4594767A1 (en) | 2022-09-28 | 2023-09-28 | Physiological navigation in mri using inherent slab-selective free-induction-decay signal | 
Country Status (2)
| Country | Link | 
|---|---|
| EP (1) | EP4594767A1 (en) | 
| WO (1) | WO2024072976A1 (en) | 
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title | 
|---|---|---|---|---|
| DE102009055961B4 (en) * | 2009-11-27 | 2018-10-18 | Siemens Healthcare Gmbh | Method for detecting a respiratory movement of an examination subject of corresponding signal data by means of magnetic resonance | 
| US9797974B2 (en) * | 2013-01-30 | 2017-10-24 | The Board Of Trustees Of The Leland Stanford Junior University | Nonrigid motion correction in 3D using autofocusing with localized linear translations | 
- 
        2023
        
- 2023-09-28 WO PCT/US2023/033997 patent/WO2024072976A1/en not_active Ceased
 - 2023-09-28 EP EP23793572.1A patent/EP4594767A1/en active Pending
 
 
Also Published As
| Publication number | Publication date | 
|---|---|
| WO2024072976A1 (en) | 2024-04-04 | 
Similar Documents
| Publication | Publication Date | Title | 
|---|---|---|
| CN102078196B (en) | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance | |
| US10222443B2 (en) | Generating motion information for an at least partially moving examination region | |
| CN103654781B (en) | Check that object carries out method and the magnetic resonance equipment of nuclear magnetic resonance by excitation | |
| JP6141883B2 (en) | MRI with motion compensation using navigator acquired using Dixon method | |
| JP3668816B1 (en) | Magnetic resonance imaging system | |
| CN102596026B (en) | MR imaging apparatus and MR imaging method | |
| EP3408679B1 (en) | System and method for optimized diffusion-weighted imaging | |
| JP2007229443A (en) | Magnetic resonance imaging apparatus and imaging condition setting method in magnetic resonance imaging apparatus | |
| JP6549612B2 (en) | Improved multiphase dynamic contrast magnetic resonance imaging method | |
| CN102743172A (en) | Magnetic resonance method and apparatus for triggered acquisition of magnetic resonance measurement data | |
| US20130307542A1 (en) | High-throughput and motion insensitive mri accelerated with multi-echo echo planar acquisition and related systems | |
| US10470685B2 (en) | Method and apparatus for capturing magnetic resonance image | |
| JP6708525B2 (en) | Magnetic resonance measuring apparatus and image processing method | |
| US20180271400A1 (en) | Magnetic resonance 2d navigator technique | |
| JPWO2011034004A1 (en) | Magnetic resonance imaging apparatus and gradient magnetic field application method | |
| JP2015213744A (en) | Magnetic resonance imaging system | |
| Grimm et al. | Optimal channel selection for respiratory self-gating signals | |
| JP2007185300A (en) | Magnetic resonance imaging apparatus | |
| CN102949190A (en) | Method and apparatus for simultaneously generating multi-type magnetic resonance images | |
| KR101787369B1 (en) | Medical imaging apparatus, magnetic reasonance imaging apparatus and method of controlling the same | |
| US10962617B2 (en) | Methods and apparatus for scan time reductions in magnetic resonance imaging using outer volume supression | |
| JP2004209084A (en) | Inspection apparatus using nuclear magnetic resonance | |
| JP6742758B2 (en) | Magnetic resonance device and program | |
| EP4594767A1 (en) | Physiological navigation in mri using inherent slab-selective free-induction-decay signal | |
| CN108261199A (en) | For acquiring the MRI system and method for time resolution MR imagings to perform object using grouped data | 
Legal Events
| Date | Code | Title | Description | 
|---|---|---|---|
| STAA | Information on the status of an ep patent application or granted ep patent | 
             Free format text: STATUS: UNKNOWN  | 
        |
| STAA | Information on the status of an ep patent application or granted ep patent | 
             Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE  | 
        |
| PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase | 
             Free format text: ORIGINAL CODE: 0009012  | 
        |
| STAA | Information on the status of an ep patent application or granted ep patent | 
             Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE  | 
        |
| 17P | Request for examination filed | 
             Effective date: 20250425  | 
        |
| AK | Designated contracting states | 
             Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC ME MK MT NL NO PL PT RO RS SE SI SK SM TR  |