WO2023034044A1 - Correction de distorsion dynamique pour irm à l'aide de navigateurs fid - Google Patents

Correction de distorsion dynamique pour irm à l'aide de navigateurs fid Download PDF

Info

Publication number
WO2023034044A1
WO2023034044A1 PCT/US2022/040917 US2022040917W WO2023034044A1 WO 2023034044 A1 WO2023034044 A1 WO 2023034044A1 US 2022040917 W US2022040917 W US 2022040917W WO 2023034044 A1 WO2023034044 A1 WO 2023034044A1
Authority
WO
WIPO (PCT)
Prior art keywords
mri
time
scan
fid
field
Prior art date
Application number
PCT/US2022/040917
Other languages
English (en)
Inventor
Tess E. WALLACE
Simon K. Warfield
Onur Afacan
Original Assignee
Children's Medical Center Corporation
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 Children's Medical Center Corporation filed Critical Children's Medical Center Corporation
Publication of WO2023034044A1 publication Critical patent/WO2023034044A1/fr

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/243—Spatial mapping of the polarizing magnetic field
    • 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/4806—Functional imaging of brain activation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices

Definitions

  • Magnetic resonance imaging is a non-invasive and versatile technique for imaging biological systems. Generally, MRI operates by detecting magnetic resonance (MR) signals, which are electromagnetic waves emitted by atoms in response to an applied electromagnetic field.
  • MR magnetic resonance
  • the detected MR signals may then be used to generate images of tissues of a patient, usually internal to the patient and unable to be directly viewed without invasive surgery.
  • a method for controlling creation of a magnetic resonance (MR) image from a magnetic resonance imaging (MRI) scan of a subject using an MRI scanner comprises, prior to completion of the MRI scan of the subject using the MRI scanner, evaluating MR signals captured during the MRI scan at a time to determine first information indicative of a state of a magnetic field within the MRI scanner at the time, and comparing the first information indicative of the state of the magnetic field within the MRI scanner at the time to second information indicative of an expected state of the magnetic field for the time.
  • the second information indicative of the expected state of the magnetic field at the time was derived from previously-captured MR data for the MRI scan of the subject using the MRI scanner.
  • the method further comprises, in response to detecting that a result of the comparing of the first information indicative of the state of the magnetic field at the time to the second information on the expected state of the magnetic field for the time satisfies at least one criterion, adjusting the MRI scanner and/or a processing of MR signals during the MRI scan based at least in part on the result.
  • the method further comprises capturing additional MR data of the MRI scan of the subject.
  • a system configured to evaluate a magnetic resonance image (MRI) scan.
  • MRI magnetic resonance image
  • the system comprises at least one processor and at least one non- transitory computer-readable storage medium storing processor executable instructions that, when executed by the at least one processor, cause the at least one processor to perform a method for controlling creation of a magnetic resonance (MR) image from a magnetic resonance imaging (MRI) scan of a subject using an MRI scanner.
  • the method comprises, prior to completion of the MRI scan of the subject using the MRI scanner, evaluating MR signals captured during the MRI scan at a time to determine first information indicative of a state of a magnetic field within the MRI scanner at the time, and comparing the first information indicative of the state of the magnetic field within the MRI scanner at the time to second information indicative of an expected state of the magnetic field for the time.
  • the second information indicative of the expected state of the magnetic field at the time was derived from previously-captured MR data for the MRI scan of the subject using the MRI scanner.
  • the method further comprises, in response to detecting that a result of the comparing of the first information indicative of the state of the magnetic field at the time to the second information on the expected state of the magnetic field for the time satisfies at least one criterion, adjusting the MRI scanner and/or a processing of MR signals during the MRI scan based at least in part on the result.
  • the method further comprises capturing additional MR data of the MRI scan of the subject.
  • MRI magnetic resonance imaging
  • MRI magnetic resonance image
  • the MRI system comprises an MRI scanner, at least one processor, and at least one non-transitory computer-readable storage medium storing executable instructions that, when executed by the at least one processor, cause the at least one computer hardware processor to perform a method for controlling creation of a magnetic resonance (MR) image from a magnetic resonance imaging (MRI) scan of a subject using an MRI scanner.
  • the method comprises, prior to completion of the MRI scan of the subject using the MRI scanner, evaluating MR signals captured during the MRI scan at a time to determine first information indicative of a state of a magnetic field within the MRI scanner at the time, and comparing the first information indicative of the state of the magnetic field within the MRI scanner at the time to second information indicative of an expected state of the magnetic field for the time.
  • the second information indicative of the expected state of the magnetic field at the time was derived from previously-captured MR data for the MRI scan of the subject using the MRI scanner.
  • the method further comprises, in response to detecting that a result of the comparing of the first information indicative of the state of the magnetic field at the time to the second information on the expected state of the magnetic field for the time satisfies at least one criterion, adjusting the MRI scanner and/or a processing of MR signals during the MRI scan based at least in part on the result.
  • the method further comprises capturing additional MR data of the MRI scan of the subject.
  • FIG.1 is a schematic diagram of a magnetic resonance imaging (MRI) facility for conducting an MRI scan, in accordance with some embodiments of the technology described herein.
  • FIG.2 is a flowchart of an illustrative process 200 of controlling generation of an MR image from an MRI scan, in accordance with some embodiments of the technology described herein.
  • Fig.3 is a pulse timing diagram showing example MRI scans with which some embodiments may operate.
  • FIG.4 illustrates a comparison of techniques operating in accordance with some embodiments described herein with other techniques.
  • FIG.5 is a schematic of an illustrative process that may be implemented in some embodiments for on-scanner calibration.
  • FIG.6 is a schematic of an illustrative process that may be implemented in some embodiments for dynamic distortion correction.
  • FIG.7 includes graphs of first- and second-order field offsets that may be used in some experiments in connection with some embodiments.
  • FIG.8 includes graphs of zeroth-, first-, and second-order field measurements that may be used in some experiments in connection with some embodiments.
  • FIG.9 depicts differences in signal intensity in a representative slice from an experiment in connection with some embodiments.
  • FIG.10 illustrates a comparison in voxel shifts that may be performed in an experiment in connection with some embodiments.
  • FIG.11 includes boxplot graphs for residual voxel shifts that may be calculated in an experiment in connection with some embodiments.
  • FIG.12 lists values for an example experiment in connection with some embodiments.
  • FIG.13 depicts maps of temporal signal-to-noise ratio that were calculated for an example experiment in connection with some embodiments.
  • FIG.14 depicts temporal signal-to-noise ratio distribution calculated for an experiment in connection with some embodiments.
  • FIG.15 depicts a temporal standard deviation distribution calculated for an experiment in connection with some embodiments.
  • FIG.16 depicts temporal signal-to-noise ratio distributions that were calculated for an experiment in connection with some embodiments.
  • FIG.17 includes boxplot graphs for temporal signal-to-noise ratio gains that were calculated for an experiment in connection with some embodiments.
  • FIG.18 depicts first-order shim coefficients and frequency offsets that were calculated for an experiment in connection with some embodiments.
  • FIG.19 depicts first-order shim coefficients and frequency offsets that were measured in an experiment in connection with some embodiments.
  • FIG.20 depicts zeroth- and first-order B0 shim offsets measured in an experiment in connection with some embodiments.
  • FIG.21A shows measured ⁇ B0 field maps from an experiment in connection with some embodiments.
  • FIG.21B includes graphs of shim coefficients from an experiment in connection with some embodiments.
  • FIG.22 shows zeroth- and first-order shim coefficients measured in an experiment in connection with some embodiments.
  • FIG.23 includes graphs of field fluctuations from an experiment in connection with some embodiments.
  • FIGs.24-25 depict images generated in experiments in connection with some embodiments.
  • FIG.26 includes boxplot graphs for values calculated in an experiment in connection with some embodiments.
  • FIGs.27-31 depict images generated in some experiments in connection with some embodiments.
  • FIG.32 is a schematic diagram of an illustrative computing device with which aspects described herein may be implemented.
  • MR magnetic resonance
  • MRI magnetic resonance imaging
  • some embodiments described herein relate to identifying changes that occur in a magnetic field within an MRI scanner during a course of an MRI scan, such as changes in magnetic field that are triggered by movement of the subject, biological functions of the subject (e.g., respiration), variation in operation of components of an MRI scanner during use of the components (e.g., operational changes that may be caused by variations in component temperature, electrical charge, inductance, or other characteristics during operation of the component), or other causes.
  • changes in a magnetic field within an MRI scanner during a course of an MRI scan such as changes in magnetic field that are triggered by movement of the subject, biological functions of the subject (e.g., respiration), variation in operation of components of an MRI scanner during use of the components (e.g., operational changes that may be caused by variations in component temperature, electrical charge, inductance, or other characteristics during operation of the component), or other causes.
  • techniques for conducting an MRI scan and for processing MR signals to generate an MR image presume an unchanging magnetic field, which may not be accurate
  • an MRI scanner may be adjusted to change a way in which additional MR data is captured and/or a processing of captured MR signals (that are to be used in generating an MR image) may be adjusted.
  • the detection of changes in the magnetic field and the adjusting of the MRI scanner and/or the adjusting of the processing of acquired MR signals may be performed by the MRI scanner, during the course of an MRI scan and prior to completion of the MRI scan.
  • the adjusting may be done by such another device, such as in response to a communication from the MRI scanner regarding a change in the magnetic field. Examples of techniques for detecting such changes in a magnetic field are described in more detail below.
  • information indicative of a state of the magnetic field at a time during the MRI scan is compared to information indicative of an expected state of the magnetic field at the time, to determine whether there is variation and/or a degree of variation. If a result of a comparing of the information regarding the state and the information regarding the expected state satisfies one or more criteria, adjustments may be made to the MRI scanner and/or to the processing of captured MR signals.
  • the information regarding the state of the magnetic field may be obtained from processing of free-induction decay (FID) navigator signals captured during the MRI scan.
  • FID free-induction decay
  • the information regarding an expected state of the magnetic field may be information corresponding to reference data captured during the MRI scan, such as from an initial MRI scan that in some cases may be performed faster and at lower resolution than the portions of the MRI scan from which diagnostic-quality images are to be generated.
  • the reference data may be used to calculate expected values for FID navigator signals. Differences between the values for FID navigator signals captured at a time during an MRI scan and expected values for the FID navigator signals for that time may be indicative of changes in the magnetic field within the MRI scanner. Accordingly, such information on FID navigator signals may be used in some embodiments to trigger adjustments to the MRI scanner and/or to processing of MR signals, as discussed below.
  • MRI is a non-invasive, versatile technique for studying physiology and pathology of a body, such as a human body or other animal body.
  • MRI exploits the nuclear magnetic resonance (NMR) phenomenon to distinguish different structures, phenomena, or characteristics of a subject.
  • NMR nuclear magnetic resonance
  • MRI may be employed to distinguish between various tissues, organs, anatomical anomalies (e.g., tumors), and/or to image diffusion, blood flow, blood perfusion, etc.
  • MRI operates by manipulating spin characteristics of atomic nuclei of subject material.
  • MRI techniques include aligning the spin characteristics of nuclei of the material being imaged using a generally homogeneous magnetic field and perturbing the magnetic field with a sequence of radio frequency (RF) pulses.
  • RF radio frequency
  • one or more resonant coils may be provided proximate a subject positioned within the magnetic field.
  • the RF coils are adapted to generate RF pulses, generally in the form of RF pulse sequences adapted for a particular MRI application, to excite the nuclei and cause the spin to precess about a different axis (e.g., about an axis in the direction of the applied RF pulses).
  • RF pulses When an RF pulse subsides, spins gradually realign with the magnetic field, releasing energy that can be measured to capture NMR data about the internal structure of the subject. Measurement of the energy of the signals can generated data indicative of nuclei excitation and relaxation.
  • complex digital signal processing techniques may be implemented by a device (e.g., through execution of software implementing the signal processing techniques) to generate, from variations within that nuclei excitation and relaxation data, information on different biological materials that were subjected to the MRI scan and that correspond to those variations, which in the context of biological subjects may be information about tissues or other anatomy of the subject.
  • Such complex signal processing may be used to generate MR images of the subject, such as images of anatomy of a biological subject.
  • Different types of MRI techniques are used to generate images for studying specific structures and functions.
  • Structural MRI may be used to detect tissues within a body, such as different organs or other tissues within the body or tissue variations within an organ.
  • fMRI Functional MRI
  • diffusion MRI is two MRI techniques that can be used to generate image data for studying anatomical structures.
  • fMRI measures the small changes in blood flow that occur with metabolic activity
  • diffusion MRI assesses the microstructural properties of tissue based on the dispersion of water molecules. These can be used, for example, to study small structures within the brain or brain activity.
  • Generating MRI data that is sufficient for accurate analysis may depend on the stillness of the subject being scanned. Motion of a subject during an MRI scan can result in image artifacts including ghosting, blurring, and signal dropout that can significantly degrade data quality. Motion of a subject during fMRI can disrupt the blood oxygenation level dependent (BOLD) signal measurement, which is used for measuring neural activity.
  • BOLD blood oxygenation level dependent
  • Such motion effects can be subtle and not easy to identify during acquisition but can jeopardize imaging quality. If image quality from an MRI scan is sufficiently degraded, it may not be possible for a radiologist to make a diagnosis using the image. To increase the likelihood of acquiring motion- free data, the imaging technologist or radiologist may acquire additional data (e.g., additional volumes) during an MRI scan, in the hopes that at least one of the acquired volumes will be motion-free. Or if a radiologist determines following an MRI scan that a resulting image is of insufficient quality, it may be necessary to repeat the MRI scan and, in some cases, sedate the patient during the repeated MRI scan to keep them still.
  • additional data e.g., additional volumes
  • That significant time lag between a first time MRI scan and any decision that a second or subsequent MRI scan is needed can have a significant negative impact on the patient, the hospital, the imaging technologies, and others.
  • An MRI scan is performed by an imaging technologist operating an MRI system, and an image is obtained through processing of the MR data that is obtained by the MRI system during the MRI scan.
  • the imaging technologist does not read the images or make any diagnoses. Instead, the detailed images generated using MRI are interpreted by a radiologist, such as to help diagnose or monitor treatment for a variety of conditions.
  • the radiologist may guide imaging technologists in performing MRI scans and generating images for the radiologist’s evaluation, but the radiologist does not operate the MRI system and often is not present during the MRI scan. Instead, the radiologist may be located in another area (perhaps not even in the same building or geographic area) and may be reviewing many images resulting from many MRI scans to perform diagnoses. In an attempt to prevent motion from degrading image quality, conventionally, an imaging technologist may visually monitor and provide instructions to a patient during the MRI scan. Due to the architecture of many MRI scanners, it can be difficult for a technologist to see if the subject in the MRI machine is moving, preventing the technologist from providing instructions to the subject to keep still.
  • an imaging technologist (who is trained to operate the MRI system) is not trained or qualified to identify those artifacts in the image or determine whether the image can still be used by a radiologist to make a diagnosis. Furthermore, the imaging technologist is not trained to select alternative strategies for addressing the imaging needs. Providing these skills to every imaging technologist would require a substantial investment in training and education, as well as a long period of practical experience. Thus, after the MRI scan, a radiologist evaluates a resulting image to determine if it is of diagnostic quality or if another MRI scan needs to be acquired.
  • a radiologist may not be immediately available to assess image quality, as they may be interpreting or managing other cases. This may incur delays of fifteen minutes or more, in addition to the time already taken for completing the initial MRI scan or the time needed for repeat MRI scans to be completed and evaluated. During that time, in many cases, a patient may have left the hospital or, even if in the hospital, may have been kept waiting for a time the patient may find unacceptable. While subject movement may be one source of artifacts in an MR image, it is not the sole source. Artifacts may arise as a result of biological activity within the subject, such as respiration/breathing or gastrointestinal activity.
  • MRI scanner hardware may vary during operation, such as due to changes in operating condition like changes in temperature, electrical charge, inductance, or other characteristics, or from routine wear and tear on components.
  • Such hardware variations may be associated with gradient heating and with eddy currents, which can affect quality of MR images resulting from a scan.
  • the inventors have recognized and appreciated that such artifacts may be related to perturbations in a magnetic field within an MRI scanner, and thus have recognized and appreciated that if such perturbations in the magnetic field could be detected during the MRI scan, steps could be taken to adjust the MRI scanner and/or adjust processing of MR signals captured during the MRI scan.
  • Such steps could enable accounting for the perturbations and mitigating the effects of the perturbations or artifacts in an MR image resulting from the scan.
  • MR signals captured during a scan may be processed once captured, as part of generating an MR image from the MR signals.
  • One aspect of the processing of MR signals may be accounting for distortion of the magnetic field, such as by modifying the signals to offset that distortion.
  • One approach to accounting for such distortion may be presuming that the distortion remains constant during the MRI scan and performing the same processing/modifying over the course of the MRI scan.
  • Another approach includes adjusting the manner of processing MR signals to account for changes in the magnetic field and changes in the distortion of the magnetic field, such that different processing/modifying of MR signals to account for perturbations in the magnetic field may be performed on MR signals acquired at different times during the MRI scan.
  • a manner in which the MRI scanner is operated to perform the MRI scan may be adjusted.
  • an MRI scanning sequence with which the MRI scanner is configured, which defines operations to be performed for the MRI scan may be modified. Modifying the MRI scanning sequence may include performing shim updates to the MRI scanning sequence, modifying frequency and/or gradient offsets for the MRI scanning sequence, and/or updating the current through external shim coils.
  • a reference scan may be performed to obtain information regarding a state of the magnetic field within the MRI scanner at a first time. This may be done, for example, at a start of an MRI scan. This information on the state of the magnetic field at the first time may set an expectation of the state of the magnetic field at a later time, if there are no changes. Subsequently, information on a state of the magnetic field may be obtained at a later time and compared to the information on the expected state of the magnetic field at that time. If the current state differs from the expected state, there has been a change in the magnetic field.
  • a change may be made to the MRI scanner and/or a manner of processing MR signals.
  • the information regarding a current state of the magnetic field at a time during the MRI scan may be obtained from processing of free-induction decay (FID) navigator signals captured during the MRI scan.
  • FID navigator may be embedded in an MRI scanning sequence for an MRI scan prior to an RF excitation pulse, following an RF excitation pulse, or at any suitable position in an MRI sequence, as embodiments are not limited in this respect.
  • the FID navigator may be embedded in the sequence following the RF excitation pulse, to measure the free induction decay following the pulse.
  • the pulse that the FID navigator follows may be a pulse used for imaging, a pulse used for the purpose of acquiring an FID navigator signal, or other pulse at another suitable position in the MRI sequence.
  • a signal analysis facility may instruct the MRI system to perform an MRI scan with one or more sequences including FID navigator signals.
  • such instructions may specify parameters for acquiring an MRI sequence.
  • the parameters may specify a repetition time (TR), a time to echo (TE), a flip angle, receiver bandwidth, field-of- view (FOV), a spatial resolution, a generalized auto-calibrating partially parallel acquisition (GRAPPA) acceleration factor, a total acquisition time, or any other suitable parameter, as aspects of the technology described herein are not limited in this respect.
  • sequence parameters may be chosen to match those typically used in clinical studies.
  • FID navigators are sensitive to changes in a magnetic field within an MRI scanner, such as due to subject motion, due to the localized spatial sensitivities of receiver coils in a multi-channel receiver array of the MRI scanner.
  • the inventors have also recognized that FID navigator signals may be captured by MRI scanner receiver coils without any spatial gradient encoding, which can mean the FID navigator signals may be captured and processed more quickly than some other MR signals. As such, using such FID navigator signals in this way may enable faster acquisition of information regarding a magnetic field, with less impact on the magnetization or overall scan time.
  • the information regarding an expected state of the magnetic field may be information corresponding to reference data captured during the MRI scan, such as from an initial MRI scan that in some cases may be performed faster and at lower resolution than the portions of the MRI scan from which diagnostic-quality images are to be generated.
  • the reference scan may be obtained using MR signals other than FID navigator signals, in some embodiments, but may be used to calculate expected values for FID navigator signals at a later time.
  • information on an expected state of a magnetic field may include information on values expected to result from acquisition of FID navigator signals at the later time.
  • FIG.1 is a block diagram of an example system 100 for evaluating MRI signals of an MRI scan to make estimations regarding an image that would result from the MRI scan, in accordance with some embodiments of the technology described herein.
  • system 100 includes an MRI system 110 (which may be an MRI scanner), an MRI system console 120, and a remote system 130.
  • system 100 is illustrative and that a system may have one or more other components of any suitable type in addition to or instead of the components illustrated in FIG.1.
  • one or more of the MRI system 110, the MRI system console 120, and the remote system 130 may be communicatively connected by a network 140.
  • the network 140 may be or include one or more local and/or wide-area, wired and/or wireless networks, including a local-area or wide-area enterprise network and/or the Internet. Accordingly, the network 140 may be, for example, a hard-wired network (e.g., a local area network within a healthcare facility), a wireless network (e.g., connected over Wi-Fi and/or cellular networks), a cloud-based computing network, or any combination thereof.
  • a hard-wired network e.g., a local area network within a healthcare facility
  • a wireless network e.g., connected over Wi-Fi and/or cellular networks
  • cloud-based computing network e.g., a cloud-based computing network, or any combination thereof.
  • the MRI system 110 and the MRI system console 120 may be located within the same healthcare facility and connected directly to each other or connected to each other via the network 140, while the remote system 130 may be located in a remote healthcare facility and connected to the MRI system 110 and/or the MRI system console 120 through the network 140.
  • the MRI system 110 may be configured to perform MR imaging of anatomy of a patient 102.
  • the MRI system 110 may include a B 0 magnet 112, gradient coils 114, and radio frequency (RF) transmit and receive coils 116 configured to act in concert to perform MR imaging.
  • B 0 magnet 112 may be configured to generate the main static magnetic field (B 0 ) during MR imaging.
  • the B 0 magnet 112 may be any suitable type of magnet that can generate a static magnetic field for MR imaging.
  • the B0 magnet 112 may include a superconducting magnet, an electromagnet, and/or a permanent magnet.
  • the B 0 magnet 112 may be configured to generate a static magnetic field having a particular field strength.
  • the B0 magnet 112 may be a magnet that can generate a static magnetic field having a field strength of 1.5T, or, in some embodiments, a field strength greater than or equal to 1.5T and less than or equal to 3.0T, or a field strength greater than or equal to 1.5T and less than or equal to 7.0T.
  • gradient coils 114 may be arranged to provide one or more gradient magnetic fields.
  • gradient coils 114 may be arranged to provide gradient magnetic fields along three substantially orthogonal directions (e.g., x, y, and z).
  • the gradient magnetic fields may be configured to, for example, provide spatial encoding of MR signals during MR imaging.
  • Gradient coils 114 may comprise any suitable electromagnetic coils, including discrete wire windings coils and/or laminate panel coils.
  • RF transmit and receive coils 116 may be configured to generate RF pulses to induce an oscillating magnetic field, B1, and/or to receive MR signals from nuclear spins within a target region of the imaged subject during MR imaging.
  • the RF transmit coils may be configured to generate any suitable types of RF pulses useful for performing MR imaging.
  • RF transmit and receive coils 116 may comprise any suitable RF coils, including volume coils and/or surface coils.
  • the MRI system 110 may optionally include image generator 118.
  • Image generator 118 may be configured to generate images based on MR data acquired by the MRI system 110 during MR imaging of the patient 102.
  • image generator 118 may be configured to perform image reconstruction to generate images in the image domain based on MR data in the spatial frequency domain (e.g., MR data comprising data describing k-space).
  • Image generator 118 may be implemented as hardware, software, or any suitable combination of hardware and software, as embodiments are not limited in this respect.
  • image generator 118 may be integrated with signal analysis facility 122, such as by being implemented in a same software program or collection of software programs to be executed on the MRI system console 120 and/or the MRI system 110, or another device such as a remote server or servers executing operations of the signal analysis facility 112 and/or image generator 118.
  • MR signal processing performed by the image generator 118 (and/or by the signal analysis facility 122 in some embodiments) may be adjusted during an MRI scan, such as in response to detected changes in a magnetic field (e.g., B 0 ) within the MRI system 110 during an MRI scan.
  • MRI facility 100 includes MRI system console 120 communicatively coupled to the MRI system 110.
  • MRI system console 120 may be any suitable electronic device configured to send instructions and/or information to MRI system 110, to receive information from MRI system 110, and/or to process obtained MR data.
  • MRI system console 120 may be a fixed electronic device such as a desktop computer, a rack-mounted computer, or any other suitable fixed electronic device.
  • MRI system console 120 may be a portable device such as a laptop computer, a smart phone, a tablet computer, or any other portable device that may be configured to send instructions and/or information to MRI system 110, to receive information from MRI system 110, and/or to process obtained MR data.
  • Some embodiments may include a signal analysis facility 122.
  • Signal analysis facility 122 may be configured to analyze MR data obtained by MRI system 110 from an MR imaging procedure of patient 102. Signal analysis facility 122 may be configured to, for example, analyze the obtained MR data to identify potential perturbations in a magnetic field within the MRI system 110 (e.g., perturbations of B0) and responsively trigger adjustments to operation of the MRI system 110 and/or processing of MR signals by the image generator 118, as described herein. Signal analysis facility 122 may be implemented as hardware, software, or any suitable combination of hardware and software, as embodiments are not limited in this respect.
  • the signal analysis facility 122 may be implemented in the MRI system console 120, such as by being implemented in software (e.g., executable instructions) executed by one or more processors of the MRI system console 120.
  • the signal analysis facility 122 may be additionally or alternatively implemented at one or more other elements of the system 100 of FIG.1.
  • the signal analysis facility 122 may be implemented at the MRI system 110 and/or the remote system 130 discussed below.
  • the signal analysis facility 122 may be implemented at or with another device, such as a device located remote from the system 100 and receiving data via the network 140.
  • MRI system console 120 may be accessed by MRI user 124 in order to control MRI system 110 and/or to process MR data obtained by MRI system 110.
  • the MRI user 124 may be, for example, an imaging technologist.
  • MRI user 124 may implement an MR imaging process by inputting one or more instructions into MRI system console 120 (e.g., MRI user 124 may select an MR imaging process from among several options presented by MRI system console 120).
  • MRI user 124 may implement an MR data analysis procedure by inputting one or more instructions into MRI system console 120 (e.g., MRI user 124 may select MR data instances to be analyzed by MRI system console 120).
  • MRI system console 120 also interacts with remote system 130 through network 140, in some embodiments.
  • Remote system 130 may be any suitable electronic device configured to receive information (e.g., from MRI system 110 and/or MRI system console 120) and to display generated images for viewing.
  • the remote system 130 may be remote from the MRI system 110 and MRI system console 120, such as by being located in a different room, wing, or building of a facility (e.g., a healthcare facility) than the MRI system 110, or being geographically remote from the system 110 and console 120, such as being located in another part of a city, another city, another state or country, etc.
  • remote system 130 may be a fixed electronic device such as a desktop computer, a rack-mounted computer, or any other suitable fixed electronic device.
  • remote system 130 may be a portable device such as a laptop computer, a smart phone, a tablet computer, or any other portable device that may be configured to receive and view generated images and/or to send instructions and/or information to MRI system console 120.
  • remote system 130 may receive information (e.g., MR data analysis results, generated images) from MRI system console 120 and/or MRI system 110 over the network 140.
  • a remote user 132 e.g., a radiologist or other clinician, such as the patient’s medical clinician
  • FIG.2 is a flowchart of an illustrative process 200 for controlling creation of an MRI image from an MRI scan, in accordance with some embodiments of the technology described herein.
  • Process 200 may be implemented by a signal analysis facility, such as the facility 122 of FIG.1.
  • the process 200 may be performed by a computing device configured to send instructions to an MRI system and/or to receive information from an MRI system (e.g., MRI system console 120 executing signal analysis facility 122 as described in connection with FIG.1).
  • the process 200 may be performed by an MRI scanner / MRI system, such as by image generator 118, or by one or more devices/processors located remotely (e.g., as part of a cloud computing environment, as connected through a network) from the MRI system that obtained MR data from the MRI system, such as via network 140.
  • an MRI scanner may be configured to execute an MRI scanning sequence.
  • An MRI sequence may include a T1-weighted scan, a T2- weighted scan, a functional MRI scan, a diffusion weighted scan, or any suitable MRI sequence, as aspects of the embodiments described herein are not limited in this respect.
  • an MRI sequence may further include one or more embedded navigators, such as FID navigators.
  • An FID navigator may be embedded in a sequence prior to an RF excitation pulse, following an RF excitation pulse, or at any suitable position in an MRI sequence, as aspects of the embodiments described herein are not limited in this respect.
  • the FID navigator may be embedded in the sequence following the RF excitation pulse, to measure the free induction decay signal following the pulse.
  • the pulse that the FID navigator follows may be a pulse used for imaging, a pulse used for the purpose of acquiring an FID navigator signal, or other pulse at another suitable position in the MRI sequence.
  • the MRI system may be configured to perform an MRI scan with one or more FID-navigated sequences.
  • the configuration may specify parameters for acquiring an MRI sequence.
  • the parameters may specify a repetition time (TR), a time to echo (TE), a flip angle, receiver bandwidth, field-of-view (FOV), a spatial resolution, a generalized auto- calibrating partially parallel acquisition (GRAPPA) acceleration factor, a total acquisition time, or any other suitable parameter, as aspects of the technology described herein are not limited in this respect.
  • sequence parameters may be chosen to match those typically used in clinical studies.
  • an MRI scanning sequence may specify acquisition settings for capturing FID navigator signals and/or the MRI sequence using the MRI system.
  • the sequence may specify an imaging trajectory, such as a Cartesian sampling trajectory with the center of the phase encoding direction acquired halfway through the MRI scan.
  • the sequence may specify whether acquisitions are non-selective.
  • the sequence may specify an acquisition plane, such as the sagittal, coronal, or transverse plane.
  • Process 200 begins at block 202, in which the MRI scanner is operated (e.g., in accordance with an MRI scanning sequence) to perform an initial scan to capture reference MR data for the scan.
  • the reference MR data may be used to determine a state of a magnetic field within the MRI scanner with the subject positioned within the scanner and having not yet moved.
  • This reference MR data may be subsequently used, in accordance with techniques described herein, to determine whether the magnetic field has changed, which may be associated with subject movement or other causes discussed above.
  • the MRI scanner may be operated in block 202 to perform a lower-resolution or lower-quality MRI scan than the configuration of the MRI scanner indicates for obtaining MR data for diagnostic purposes.
  • Such lower-resolution or lower-quality scans may be obtained, for example, by adjusting a number of encoding lines per shot as compared to what will be performed for obtaining MR data for diagnostic purposes.
  • the signal analysis facility analyzes the reference MR data to determine information on an expected state of the magnetic field, for subsequent analysis.
  • the signal analysis facility moves to block 204 and performs iterations of several blocks of the flowchart of FIG.2, during the MRI scan.
  • the operations may be performed during the MRI scan such that they are performed prior to completion of acquisition of MR data for the scan.
  • the operations may also be performed during the MRI scan such that each iteration is performed in association with a slice (e.g., for every slice, or for some slices) and is performed before MR data for the slice is output from the MRI scanner or from a system that analyzes MR data as it is acquired by the MRI scanner.
  • the operations may be performed in real time, as MR data is acquired and prior to output of the MR data to an element that will compile MR data and generate imagery.
  • Such operations may include operations that are performed “online” in the MRI scanner environment, as opposed to offline.
  • the MRI scanner executes the scanning sequence and, as part of the scanning sequence, obtains FID navigator signals.
  • the MRI scanner analyzes the FID navigator signals to determine information on a state of the magnetic field at the time the FID navigator signals were acquired, and determines in block 208 whether the state of the magnetic field differs from the expected state of the magnetic field as determined from the reference data. If not, then the process 200 continues and the signal analysis facility loops back to block 204 to acquire more MR data as part of the MRI scanning sequence.
  • the signal analysis facility determines in block 208 that the current state of the magnetic field is different from the expected state of the magnetic field, the signal analysis facility triggers adjustments to be made in view of the difference.
  • the signal analysis facility may change the MRI scan or MRI scanner, such as by reconfiguring the MRI scanner. Such a reconfiguration may be performed through changing the MRI scanning sequence, which may include changing one or more instructions of it.
  • the reconfiguring of the MRI scanner may include changing B0 shim offsets, updating current through external shim coils of the MRI scanner, modifying frequency or gradient offsets, and/or other operations.
  • the signal analysis facility may adjust a processing of MR signals acquired by the MRI scanner during the MRI scan.
  • the processing that is adjusted may include processing prior to output from the MRI system of MR data acquired by the MRI system at a time, such as the “online” process for analyzing MR signals.
  • the adjusting of the MR signal processing may include adjusting a spatially- resolved map of magnetic field distortion to reflect changes in field distortion determined from the analysis of the FID navigator signals analyzed by the signal analysis facility.
  • the process 200 continues to block 204 for the MRI scanner to acquire more data following the adjusting, to complete the scanning sequence.
  • the process 200 ends once the MR data to be acquired (as indicated by the scanning sequence) has been acquired.
  • Example Embodiments Described below are certain non-limiting examples of techniques for evaluating MR signals such as FID navigators during an MRI scan and, during the MRI scan, adjusting the MRI scanner and/or a processing of MR signals to generate an MR image. These are intended to be illustrative of various ways in which techniques described herein may operate. It should be appreciated that embodiments are not limited to operating in accordance with these examples.
  • First Set of Example Embodiments there is provided an approach for slice-wise dynamic distortion correction for echo planar imaging (EPI) using FID navigators (FIDnavs).
  • EPI echo planar imaging
  • FIDnavs FID navigators
  • a low-resolution reference image is used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B0 field inhomogeneity variations up to second order from measured FIDnavs.
  • Single-shot EPI forms the basis of many functional MRI studies as it facilitates measurement of BOLD signal fluctuations related to neuronal activity with high temporal resolution.
  • EPI is highly sensitive to magnetic field (B 0 ) inhomogeneities that occur at interfaces between materials with different magnetic susceptibilities and due to hardware imperfections.
  • Uncompensated variations in the encoding field produce non-linear geometric distortions and regions of signal “pile-up” where image intensities are assigned to incorrect spatial locations, while intra-voxel dephasing during long readouts creates signal voids in regions with high inhomogeneity.
  • Active B 0 shimming is typically performed prior to imaging to help offset inhomogeneities induced by the subject-specific susceptibility distribution by superimposing an optimal set of low-spatial-order spherical harmonic (SH) fields on the static magnetic field.
  • SH low-spatial-order spherical harmonic
  • the total EPI readout time can also be reduced by employing ramp sampling and/or parallel imaging acceleration to minimize artifacts, but residual distortions remain a problem.
  • Static distortions may be corrected in post-processing by computing voxel displacements in the phase-encoding direction from a low-resolution map of field offsets, usually measured from a separate gradient- echo (GRE) acquisition with multiple echo times.
  • GRE gradient- echo
  • Dynamic spatially-resolved field maps may be estimated from the phase evolution of each EPI volume if a second echo is acquired using accelerated readouts (dynamic off-resonance correction with multi-echo acquisition (DOCMA)); however, acquisition of two volumes per repetition time limits the spatial resolution that can be achieved. This spatiotemporal constraint can be avoided by jittering the echo time or alternating the spatial encoding direction for consecutive time points. While these modifications facilitate dynamic field mapping, they fundamentally alter the fMRI acquisition, potentially inducing additional BOLD contrast variations, and are not rapid enough to accurately characterize respiration-induced B 0 changes. Alternatively, coil-dependent phase offsets may be measured in a separate reference scan and subtracted from single-echo EPI data to generate a field map at each time point.
  • DOCMA dynamic off-resonance correction with multi-echo acquisition
  • fitting SH functions to field probe measurements can allow correction for field fluctuations that occur due to both hardware imperfections and physiological mechanisms. Correction can be performed either retrospectively, or prospectively by dynamically updating the scanner shim currents in real time. Real-time shimming capability may compensate for respiration-induced B0 fluctuations; however, this approach requires subject-specific training to calibrate phase changes related to chest motion and cannot correct for other sources of field perturbations. FID navigator methods can offer a more efficient solution for measuring these low- spatial-order relative B0 variations by sampling a small additional amount of k-space data, while circumventing the need for additional instrumentation.
  • Double volumetric navigators (DvNavs), based on a pair of highly-accelerated, low-resolution 3D EPI acquisitions periodically interleaved with the host sequence, may also be used to generate a spatially resolved field map for each volume, recently demonstrated in diffusion tensor imaging; however, the time required for DvNav encoding would place a substantial burden on the fMRI acquisition. Measuring linear B0 inhomogeneities from the phase difference between three orthogonal projections with different TEs to enable real-time shimming has been proposed. Cloverleaf navigators combine motion and shim measurement in a single trajectory and measure B 0 inhomogeneities directly from shifts in the k- space domain.
  • FID navigators measure the signal from each receiver coil without any spatial encoding and thus can be acquired extremely rapidly.
  • Measuring the phase evolution between two FID signals may be used to compensate for global respiratory-induced frequency fluctuations and system frequency drifts in EPI time series; however, this may be insufficient to correct for non-linear geometrical distortions that arise from more localized B 0 field variations.
  • Described herein are embodiments of a method for dynamic distortion correction in MRI (e.g., EPI time series) that leverages the spatial encoding information provided by the multi-channel RF receive coil arrays to enable characterization of spatio-temporal B 0 fluctuations from FIDnavs embedded in each slice.
  • FIDnavs may measure a signal integral over an excited slice of the underlying spin density distribution, modulated by the B0 field inhomogeneities and RF receive coil sensitivities at each spatial location. This may be approximated by the discrete summation of a complex- valued multi-channel reference image with matched contrast properties. Dynamic changes in the B0 field from a reference time point may be expressed as a series of low-spatial-order basis functions: where r denotes the spatial coordinate and n indexes time.
  • a forward model can be generated by simulating the effect of changes in these inhomogeneity coefficients on the complex-valued FIDnav signals (y j,n ), with spatial encoding provided by a multi-channel reference image: where y j,n is the measured FIDnav from coil j at acquisition n; each S j,0 is a 1 x N p vector of the complex pixel intensities of a reference image with N p pixels, measured by coil j; y is the gyromagnetic ratio (in rad/s); ⁇ is the echo time of the reference image (designed to match the sampling time of the FIDnav); ⁇ is an N p x N b matrix representing the N b basis functions and b n is an N b x 1 vector representing the inhomogeneity coefficients at acquisition n.
  • the FIDnav signal model of this example uses multi-channel GRE reference image data with matched contrast properties to be acquired for each subject. Two datasets may be acquired with reversed gradient polarities to account for phase errors that could arise due to gradient delays and other system imperfections.
  • TRs for single-shot EPI acquisitions are typically on the order of 2-3 seconds.
  • reference data acquisition may be quickened by using a segmented EPI reference image with three encoding lines per shot (Fig. 3(A)), and with blip-up and blip-down encoding, to provide self-referencing compensation for the effects of gradient imperfections on the encoded signal.
  • Fig.3 is a pulse timing diagram showing segmented EPI reference scan (A) and FID-navigated EPI time series scan (B). Resolution requirements for the reference image data may be adjusted to further reduce the time required for forward model generation, such as to ⁇ 1 minute in some cases.
  • MRI e.g., fMRI
  • FIDnavs e.g., ultra-short FIDnavs
  • spatial encoding provided by multi-channel array coils.
  • an FIDnav-based method may be comparable in performance to, or in some cases outperform, a DOCMA approach with respect to both geometric distortion correction and improvement in BOLD signal stability. Reliable phase unwrapping and the invertibility of distortion are both crucial for the DOCMA approach. A shorter acquisition may be used to reduce the difference in TE values across the multiple echoes to enable more reliable field estimates.
  • FIG. 4 illustrates that FIDnav correction outperformed DOCMA in the frontal brain cortex, where field inhomogeneity is higher. More particularly, Fig.4 illustrates distortion correction results in a representative volunteer following correction with a static voxel shift map (VSM) and with VSMs obtained from FIDnav and DOCMA field measurements for changing Y-shim (A) and nose touching (B). The brain boundary extracted from T2 FSE (shown on the right for comparison) is overlaid on each corrected image. Arrows denote regions with residual unwarping errors, observed with both static and DOCMA correction.
  • VSM static voxel shift map
  • FIDnavs enable reliable geometric distortion correction for both shim changes and nose touching. Unwrapping errors led to temporal instability in field measurements in some subjects, which may help explain why DOCMA may not improve tSNR in the no motion scans. While the DOCMA approach facilitates dynamic distortion correction without a reference scan, it is subject to several other limitations, including the need for a second echo, which restricts the maximum spatial resolution that can be achieved within a given TR, making it impractical for high resolution fMRI. Like most data-driven dynamic distortion correction approaches, DOCMA also requires reconstruction and phase unwrapping of image data to perform field correction, limiting its use for real-time applications or for correcting phase inconsistencies between segments in multi-shot EPI.
  • Some embodiments of a FIDnav-based solution described herein may have several advantages over other shim navigator approaches.
  • Orthogonal one-dimensional navigator signals have previously been used to measure linear field changes from the phase difference between projections with different TEs in image space or shifts in k-space. These methods assume that inhomogeneities along each axis can be measured independently; however, dephasing effects can cause the phase evolution of the projection data to be not necessarily linearly related to B 0 inhomogeneities in that direction. Thus, projection data alone may not be sufficient to resolve field changes in multiple directions without a priori information provided by a reference scan. Also, the time required for spatial encoding may substantially increase the minimal TE or limit the frequency of shim updates.
  • Some embodiments reduce the time needed to acquire this reference image data by using a lower- resolution (7 mm in-plane) multi-shot EPI sequence, rather than a GRE sequence used in other contexts.
  • model calibration could be integrated in a single step using a multi-shot, multi-echo EPI reference scan.
  • the effect of respiration-induced field perturbations on the reference image data is assumed to be minimal at 3T.
  • some embodiments assume that the measured field does not change within the EPI readout and thus cannot correct for higher-frequency field changes.
  • field probes are the only method that can simultaneously monitor field changes during the EPI readout itself to compensate for trajectory deviations due to eddy currents and other hardware-related instabilities.
  • Accelerated, multi-echo EPI can be used in some embodiments to provide a reference voxel shift measurement for each volume using the DOCMA approach. Acceleration may limit the amount of geometric distortion by reducing the total readout time. Adjusting the scanner shim settings and nose touching in volunteers induced small, but measurable distortions in the images in examples described below. Using non-accelerated imaging may have made the effect of correction clearer but would have limited use of DOCMA as a comparison. In the examples discussed below, images were unwarped in SPM using either a static field estimate or dynamic field estimates from FIDnavs or DOCMA and aligned using rigid-body registration.
  • zeroth-order field may should also be partially compensated by rigid-body realignment.
  • differences in tSNR represent residual unwarping errors, slice-wise shifts, realignment errors (due to changing geometric distortion between volumes) and BOLD contrast changes.
  • Motion may be a major source of magnetic field perturbations due to the changing orientation of the susceptibility distribution as, for example, a subject’s head moves within the static magnetic field. Head motion may be reduced in some embodiments using foam padding.
  • FID navigators may be used to detect and/or to correct for geometric distortions in single-shot EPI at 3T, although it should be appreciated that embodiments may operate with other EPI-based sequences and field strengths, or other types of MRI scans.
  • FIDnav-based field measurements may be applicable in both 2D multi- shot and 3D EPI as some embodiments do not require formation of an image and may have a small impact on the timing of an MRI scanning sequence.
  • higher-order field perturbations may be corrected retrospectively by updating the voxel shift map used to unwarp each EPI image. Shim changes measured from FIDnavs could be used to update the transmit frequency and first-order shim fields in real time in some embodiments, while prospective correction of higher-order field changes is becoming feasible using custom-built shim array coils.
  • Real-time shimming could enable correction for the effects of magnetic field perturbations on RF excitation, fat suppression and intra-voxel dephasing, in addition to changing geometric distortions.
  • Dynamic distortion correction with FIDnavs may translate to improved BOLD sensitivity in fMRI studies, which could aid many applications including improving localization for presurgical planning in epilepsy patients, improving the reliability of fMRI in patients with Parkinson’s and Alzheimer’s disease, and in research studies utilizing high- resolution fMRI to measure neuronal activation with sub-millimeter spatial resolution.
  • This technique could also be translated to other EPI applications, such as perfusion and diffusion- weighted imaging and non-BOLD fMRI.
  • Some embodiments described herein thus include slice-wise dynamic distortion correction using FIDnavs embedded in each slice of a single-shot EPI sequence.
  • FIDnav-based dynamic distortion correction may provide improved geometric accuracy of images acquired with magnetic field perturbations. Improvement in tSNR of the EPI time series was comparable to, and in some cases outperformed, an existing DOCMA approach.
  • Some embodiments of an FIDnav-based approach may have one or more of several advantages over data-driven methods: they may not fundamentally alter the MRI scan and can be implemented without compromising spatiotemporal resolution. Furthermore, they may not require computationally expensive image reconstruction and phase unwrapping, making it more amenable to real-time shimming applications.
  • FID-navigated dynamic distortion correction is widely applicable and is expected to improve BOLD sensitivity in fMRI experiments at field strengths of 3T and above.
  • Second Set of Example Embodiments FIDnavs may be embedded in a gradient echo sequence and a subject-specific linear calibration model may be generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. Applying measured FIDnav shim updates may in some cases compensate for applied global and linear field offsets in MRI scans of a subject.
  • FID-navigated real-time shimming may lead to a reduction in field fluctuations and a consequent improvement in T2*-weighted image quality.
  • FIDnavs may facilitate rapid measurement and application of field coefficients for slice-wise B 0 shimming. This may counteract spatiotemporal field perturbations and improve T2*-weighted image quality, which may be helpful for a variety of clinical and research applications, particularly at ultra-high field.
  • MRI scanners may leverage a homogeneous and stable magnetic field for accurate signal excitation and spatial encoding.
  • Dynamic changes in the underlying susceptibility distribution due to respiration and subject motion can give rise to spatially and temporally varying magnetic field perturbations, which can create ghosting, ringing, and signal modulation artifacts in multi- shot images.
  • Hardware imperfections such as gradient heating and eddy currents, are further sources of artifacts. These effects may be particularly evident in gradient echo sequences with long echo times, such as T 2 *-weighted imaging, which may form the basis of susceptibility weighted imaging and quantitative susceptibility mapping.
  • T 2 *-weighted imaging While imaging at ultra-high field enhances sensitivity to microscopic susceptibility variations, which improves contrast in both magnitude and phase images, the impact of macroscopic perturbations in the background magnetic field may also increase with field strength, which can be detrimental to image quality.
  • NMR field probes may be placed around the subject to directly record the B 0 field at discrete points. A model can then be fit to these measurements to describe and compensate for low-spatial-order field changes within the imaged region.
  • Field probes can monitor field changes independent of the imaging sequence, and rapid measurements may be used to provide real-time field control.
  • Field changes during MRI can also be estimated by acquiring additional navigator signals with varying degrees of complexity. Navigator-based methods typically involve estimating phase changes across selected parts of k-space or low-resolution images.
  • interleaving navigators with the imaging sequence reduces scan efficiency, particularly for rapid gradient echo sequences with little intrinsic dead time. Even acquiring a small amount of additional k- space data interleaved with the imaging acquisition can increase scan overhead by up to 30%.
  • navigators may be acquired in a segmented manner; however, this limits the frequency with which field measurements can be made.
  • Free induction decay navigators FIDnavs
  • FIDnavs Free induction decay navigators
  • FIDnavs have been used to correct for global frequency fluctuations and with the introduction of multi-channel coil arrays, this additional spatial sensitivity information may be leveraged to characterize spatially-varying B0 field distributions.
  • the ability of sensitivity encoding to recover low-spatial-order field inhomogeneity coefficients from unencoded FID navigator points has been proposed, along with the use of contrast-matched reference image data to measure and retrospectively correct for the effects of spatiotemporal field changes in both T 2 *-weighted imaging and echo planar imaging.
  • retrospective correction can be limited to reversible alterations that are encoded in the perturbed data and therefore cannot compensate for inherent losses of information due to intra-voxel dephasing, inconsistent selective excitation, and disruption of steady state.
  • Retrospective correction also typically requires implementation of offline reconstruction pipelines, which can be computationally slow.
  • Some embodiments described herein include a method for real-time field control using rapid FIDnav measurements, which may update the central frequency and first-order shim terms for each slice of an MRI scan. These measurements can be calibrated directly on the MRI scanner using spatial encoding from a prior low-resolution reference scan, facilitating dynamic updates with minimal latency.
  • Some embodiments thus include accurate online field correction and improved T2*-weighted image quality in the presence of dynamic field perturbations at different field strengths, such as at 3T or 7T.
  • FIDnavs may characterize spatiotemporal field changes using a forward model that predicts the effect of low-spatial-order changes in B0 on the measured FIDnav signal from each RF receive coil in a multi-channel array.
  • each FIDnav signal By approximating each FIDnav signal as a summation of pixel intensities from the complex-valued multi-channel reference image data, modulated by spatial variations in the magnetic field, the impact of low-spatial-order field changes on the measured FIDnav signal may be written as follows: where is the measured FIDnav from coil j at acquisition vector of the reference complex-valued coil image pixel intensities is the gyromagnetic ratio (in rads -1 T -1 ); is the sampling time of the FIDnav; matrix representing the spherical harmonic basis functions for an image with pixels, and is an vector representing the shim coefficients at acquisition If field changes meet a criterion, such as by being sufficiently small , changes in the FIDnav signal may be assumed to be linear with respect to the field inhomogeneity coefficients.
  • a field change of 4 Hz at 3T may give ⁇ 1% approximation error.
  • the ability of FIDnavs to estimate first, second, or even higher- order field profiles may in some cases be linked to having a sufficiently large number of spatially independent RF receive profiles.
  • a typical 64-channel brain array may in some cases provide sufficient degrees of freedom to estimate B 0 fields in the head up to at least second order spherical harmonics.
  • a linear calibration model may be generated on the scanner for each subject by simulating the effect of B0 field offsets up to second order on the complex- valued multi-channel FIDnav signals in each slice.
  • Y is the complex multi-channel FIDnav signal vector
  • A is the linear calibration matrix
  • b is a vector of the (unknown) inhomogeneity coefficients at each time point.
  • GRE gradient echo
  • a multi-channel reference image may be computed by phase-correcting and combining these low-resolution images and calibration of FIDnav shim changes was performed directly on the scanner using the vendor image calculation environment.
  • a product GRE sequence may be modified to include an FIDnav every repetition time between the slice-selective excitation pulse and imaging readout.
  • the subject-specific calibration matrix computed for each slice may be loaded into memory at the beginning of the scan.
  • B 0 shim coefficients may be estimated as described above, via least-squares fitting and zeroth- and first-order shim terms were sent back to the sequence prior to calculation of the next event block (defined as the RF and gradient events comprising one TR).
  • continuous processing of event blocks may be interrupted and a processing margin of 25 microseconds may be specified to allow time for the sequence to prepare the subsequent event block with updated shim information during execution of the current event block.
  • Global B 0 changes may be compensated by updating the center frequency relative to the previous TR.
  • Linear shim corrections may be implemented as offsets to the imaging gradients, which allows shim updates to be switched as rapidly as regular gradient pulses in the sequence.
  • FIG.5 is a schematic showing on-scanner calibration of FIDnav B0 shim measurements of some embodiments, which is performed for each subject at the beginning of the acquisition.
  • a pair of low-resolution images with reversed gradient encoding are acquired to create a forward model of FIDnav signal changes.
  • Changes in B 0 shim offsets are simulated using the phase-corrected multi-channel reference data to calibrate a linear model for each subject.
  • B0 shim coefficients are estimated from the measured FIDnav signals in each slice using the pre-calculated subject- specific calibration and used to update the center frequency and gradient offsets in real time.
  • B0 field changes relative to the previous TR may be assumed to be small.
  • This can enable use of a linear calibration model, which may facilitate fast estimation of shim coefficients via least-squares fitting. Since the calculation may be quick with low computational cost, latency may be linked to the rate of data transfer between the vendor image calculation environment to update the hardware instructions.
  • Applying FID-navigated B 0 shim updates may successfully counteract field perturbations induced by manual shim changes in phantoms and volunteers, with high estimation fidelity in both X, Y and Z directions.
  • T 2 *-weighted imaging may be advantageous for a range of clinical applications including Alzheimer’s disease, Huntington’s disease, seizures, and stroke.
  • T2*- weighted imaging may also be used in a variety of research applications in studies investigating brain structure, iron deposition and the effects of aging.
  • Approaches described herein may be useful for subjects who are prone to move and breathe deeply and irregularly, such as children, elderly adults, and patients with Alzheimer’s disease.
  • the ability to mitigate dynamic field changes may be advantageous to realize the full potential of these techniques at ultra-high field strengths, which suffer from the increased effects of time-varying macroscopic field changes.
  • Mitigating artifacts due to dynamic field fluctuations may also be helpful for improving the accuracy and reproducibility of T2* mapping and quantitative susceptibility mapping. Since the impact of dynamic field changes can scale with magnetic field strength, phase stabilization may commonly be employed at 7T to mitigate global field fluctuations during the scan.
  • phase stabilization is limited to retrospectively compensate for global field fluctuations, while some techniques described herein may enable prospective correction for field changes up to first order. Phase stabilization may also lead to an overall increase in scan time as it acquires an additional echo following the imaging echo, which increases the minimum TR of the sequence.
  • Previous studies have demonstrated the utility of real-time shimming using NMR field probes, which facilitate sampling of magnetic field fluctuations with high temporal resolution and without any sequence overhead.
  • field probes sample point-wise measurements and use a model to estimate low-spatial-order spatiotemporal field changes within the imaged region.
  • FIDnavs may also measure a signal directly from the receiver coils, which may be aided by some minor sequence modifications, but may also enable field measurements without specialized hardware. Magnetic field measurements from k-space and image-based navigators have also been used to perform dynamic B0 field updates and may also be used to simultaneously measure and correct for head motion.
  • the additional scan time in some experiments was approximately 15 seconds per session, which may be reduced in some circumstances using parallel imaging and by applying a pre-calculated offset to compensate for gradient delays in the reference image.
  • spatiotemporal field changes up to first order may be compensated for, based on measurements from FIDnavs embedded in each slice of a two- dimensional GRE sequence.
  • a short processing latency makes the technique well suited to 3D susceptibility weighted imaging and quantitative susceptibility mapping protocols, which are typically acquired with TRs on the order of 50-60 ms. Shim updates could be calculated and applied to the subsequent TR to closely track the true field perturbations.
  • FID-navigated real-time shimming may also be applied in some embodiments in both single-shot and segmented echo planar imaging sequences to reduce dynamic signal fluctuations in functional and diffusion-weighted MRI. As the FIDnav signal may be dependent on the underlying sequence, application of this approach to sequences with an inconsistent steady state or other dynamic contrast changes may benefit from modifications to account for signal changes unrelated to field perturbations.
  • Some embodiments may include prospective motion correction using motion measurements from an external tracking system or additional navigators for simultaneous field stabilization and motion control. Joint retrospective correction for head motion and motion- induced field changes using dual-echo image-based navigators may improve T2*-weighted image quality and susceptibility weighted imaging, with applications to improved cortical lesion detection in multiple sclerosis. Accordingly, some embodiments include a technique for real-time shimming based on rapid measurement of spatiotemporal field changes from multi-channel FIDnav measurements in each slice.
  • Some embodiments may successfully counteract spatiotemporal field perturbations up to first order and improve T2*-weighted images acquired in the presence of intentional field modulations. Reducing dynamic field fluctuations may be advantageous for realizing more or full potential of sequences based on T 2 * contrast, particularly at high field strengths and in challenging patient groups.
  • Example Experiments First Set of Experiments All imaging experiments in the first set were performed at 3T (MAGNETOM Skyra; Siemens, Er Weg, Germany) using the vendor-supplied 64-channel head coil. Five volunteers were scanned after providing written informed consent in accordance with an institutional IRB- approved protocol. An FIDnav module was inserted into a gradient-echo EPI sequence between each slice-selective excitation and EPI readout (Fig.3B).
  • FIG.6 illustrates an example experimental design to compare geometric distortion correction using voxel shift maps (VSMs) computed from conventional static ⁇ B 0 field measurements from a multi-gradient-echo (MGRE) scan, and from dynamic FIDnav- and DOCMA-based field offsets.
  • VSMs voxel shift maps
  • MGRE multi-gradient-echo
  • Static ⁇ B0 field maps were calculated by taking the Hermitian inner product (HIP) over channels of the phase difference between multi echo images and scaling to obtain ⁇ B 0,0 (r).
  • HIP Hermitian inner product
  • Coefficients describing the change in B 0 were estimated from FIDnavs and used to compute maps of ⁇ B 0,n (r). These were superposed with the static ⁇ B 0 map to compute ⁇ B0,n(r) at each time point n.
  • maps of ⁇ B0,n(r’) were computed in distorted image space from the HIP of multi-echo EPI images.
  • VSMs were computed from each ⁇ B 0 measurement and used to unwarp the first echo EPI image.
  • Experiment 1 Phantom validation of FIDnav ⁇ B0 measurements To validate the accuracy of FIDnav-based ⁇ B0 measurements in an EPI sequence and assess the minimal requirements for forward model generation, a cylindrical water-bottle was scanned with manual adjustment of the scanner shim settings. A segmented EPI reference image was acquired with TE 5 ms, TR 1000 ms, EPI factor 3 and blip-up and blip-down encoding (total acquisition time 22 s). Axial FID-navigated multi-echo EPI scans were acquired with the following parameters: TE1/TE216.0/33.8 ms, TR 1000 ms, 10 repetitions, GRAPPA acceleration factor 2.
  • the FIDnav module sampled 64 points in 0.4 ms, centered on TNAV 5 ms. Both reference and FID-navigated scans were acquired with flip angle 90°, receiver bandwidth 2232 Hz/pixel, FOV 192 mm, voxel size 3 mm isotropic, 16 slices with 20% inter-slice gap, centered at the isocenter of the magnet. First and second-order shim settings were manually adjusted up to ⁇ 10 ⁇ T/m and ⁇ 100 ⁇ T/m2 in increments of 2.5 ⁇ T/m and 25 ⁇ T/m2, respectively. EPI images were reconstructed using Dual-Polarity GRAPPA (DPG), which incorporates Nyquist ghost correction into the parallel imaging reconstruction (34).
  • DPG Dual-Polarity GRAPPA
  • Experiment 2 Field mapping and geometric distortion correction in volunteers To compare the performance of FIDnav-based voxel shift maps to the reference DOCMA method, FID-navigated multi-echo EPI images were acquired in five volunteers with the following conditions: 1) resting position, standard shim settings; 2) applied field offset in the phase-encoding (Y) direction of ⁇ 10 ⁇ T/m; and 3) while subjects raised their arm to touch their nose.
  • Scan parameters for the FID-navigated multi-echo EPI scans were as follows: TNAV 5 ms, TE1/TE218/48.1 ms, TR 2000 ms, receiver bandwidth 1860 Hz/pixel, in-plane resolution 2 mm, GRAPPA acceleration factor 2 and 6/8 partial Fourier.
  • Posterior- anterior (PA) phase-encoding was used to avoid signal pileup in the frontal cortex (35).
  • Segmented EPI reference images with blip-up and blip-down encoding were acquired for forward model generation with TE 5 ms, TR 2000 ms, receiver bandwidth 2232 Hz/pixel and matrix size 32x32, yielding an in-plane resolution of 7 mm (total acquisition time 44 s). Both sets of scans were acquired with flip angle 70°, FOV 224 mm, 20 slices with 2.5 mm slice thickness and 20% inter-slice gap, and fat suppression.
  • a multi-echo GRE field map was acquired for static distortion correction (TE1/TE24.92/7.38 ms, TR 250 ms, flip angle 25°, receiver bandwidth 800 Hz/pixel, in-plane resolution 2 mm; matched FOV, slice thickness and spacing).
  • An axial T2-weighted fast-spin-echo (FSE) image (TR 7360 ms, TE 78 ms, receiver bandwidth 220 Hz/pixel, in-plane resolution 0.9 mm, with matched FOV, slice thickness and spacing) was acquired to provide an anatomical reference.
  • FSE fast-spin-echo
  • NRMSE Normalized root-mean-square error
  • Static distortion correction was applied using field measurements from a multi-echo GRE acquisition acquired prior to each EPI series. Dynamic distortion correction was applied by superimposing second-order field changes estimated from FIDnavs on the static field map for each time point (Fig.6). Dynamic correction using the DOCMA approach was also applied for comparison. Distortion-corrected images were aligned to the third image in each series using rigid-body registration in the Insight Toolkit (ITK) with a mean-squares similarity metric and Kaiser-Bessel window (sinc) resampling.
  • ITK Insight Toolkit
  • sinc Kaiser-Bessel window
  • Temporal SNR was computed as the ratio of the mean and standard deviation of each voxel time-course for each image series as a measure of BOLD signal stability for static, FID-navigated and DOCMA distortion correction. Improvement in mean tSNR associated with FIDnav and DOCMA dynamic distortion correction was characterized by comparing to images reconstructed using the static distortion correction scheme.
  • the relative tSNR gain was calculated as follows: where (tSNR(r)) denotes the mean tSNR over all voxels within the whole-brain mask. Significance of the one-sided Wilcoxon signed rank test for median tSNR gain above zero was tested.
  • Fig.7 Phantom validation of FIDnav LlB0 measurements A comparison between field coefficients measured using the proposed FIDnav-based approach and the DOCMA method for the applied first- and second-order shim changes is shown in Fig.7. More particularly, Fig.7 depicts applied first- and second-order field offsets in a phantom showing field inhomogeneity coefficients estimated from the proposed FIDnav-based approach and the reference DOCMA method. As shown in the graphs, FIDnav measurements accurately characterized field changes up to second order in the phantom study.
  • FIDnav measurements were computed using a reference image matrix size of 32x32, as this provided a good trade-off between the efficiency of both reference data acquisition and shim parameter computation and the accuracy of higher-order field estimates.
  • Fig.8 shows errors in zeroth-, first- and second order field measurements estimated from FIDnavs using a non-linear fit that minimizes error between the measured FIDnav signal and a forward model generated from multi-channel reference images of different resolutions. Time taken to acquire reference data and compute the non-linear fit with each resolution level are also shown.
  • a reference image with matrix size 32x32 provides a good trade-off between accuracy and time taken to acquire the reference data.
  • the proposed FIDnav-based approach accurately characterized shim changes, with mean absolute errors of 0.49 ⁇ 0.06 ⁇ T/m (first-order) and 1.22 ⁇ 0.37 ⁇ T/m2 (second- order) relative to the reference DOCMA approach for maximum changes of 10 ⁇ T/m and 100 ⁇ T/m2, respectively.
  • Maps of normalized signal intensity before and after correction for first- and second- order shim changes are shown in Fig.9.
  • Fig.9 depicts differences in normalized signal intensity in a representative slice through the water bottle phantom for applied shim offsets of Y ⁇ 10 ⁇ T/m (A) and X 2 -Y 2 ⁇ 100 ⁇ T/m2 (B), relative to an undistorted reference scan, before correction, and after correction with voxel shift maps calculated from FIDnav field inhomogeneity coefficients and DOCMA field maps.
  • FIDnav field estimates enable accurate geometric distortion correction in the presence of large field offsets, comparable to that achieved with the reference DOCMA method.
  • Fig.10 provides a comparison of voxel shift maps (VSMs) in an axial slice showing the residual geometric distortion in the brain induced by 10 T/ change in the Y-shim gradient (A) and nose touching (B) in a representative volunteer, measured from FIDnavs and the reference DOCMA method. Difference maps confirm that spatial field offsets estimated by the proposed FIDnav-based approach are in good agreement with pixel-wise estimates from multi-echo EPI. Relative to the spatially-resolved VSMs measured using the DOCMA approach, correction with FIDnav field estimates led to a substantial reduction in the residual voxel shift.
  • VSMs voxel shift maps
  • Fig.11 depicts a comparison of residual voxel shifts measured using DOCMA (“measured”) and after subtraction of voxel shifts modelled using the proposed FIDnav method (“corrected”) for all subjects in the brain region with applied Y-shim gradient (A) and with nose touching (B). Boxes delineate 25th and 75th quantiles, with whiskers extending to the maximum voxel shifts. Median and mean are denoted by solid lines and open circles, respectively.
  • Fig.13 depicts maps of temporal signal-to-noise ratio (tSNR) in a representative subject for the no motion, nose-touching and deep-breathing paradigms, comparing static distortion correction and dynamic correction using FIDnav and DOCMA field measurements.
  • Dynamic distortion correction results in a substantial tSNR improvement across all paradigms, particularly in the frontal brain region during nose touching.
  • a global increase in tSNR within the whole-brain volume is observed with both FID-navigated and DOCMA-based dynamic distortion correction across all scans.
  • the tSNR increase is particularly marked in the frontal brain regions for the nose-touching paradigm. Histograms of the tSNR distribution across all five subjects within the brain region are displayed in Figure 14.
  • Fig.14 depicts temporal signal-to-noise ratio (tSNR) distribution in the whole-brain volume measured across all five volunteers for no motion, nose touching and deep breathing with static, FIDnav and DOCMA distortion correction. Histograms are shown on a linear scale (top of Fig.14), as well as a logarithmic scale (bottom of Fig.14) to show the differences in high tSNR values more clearly.
  • the tSNR distribution for the “no motion” case is shown as a dashed line in Fig.14 for comparison.
  • the proposed FIDnav method enables an increase in tSNR for all paradigms, comparable to, or in the case of “no motion” – outperforming, the DOCMA method.
  • Fig.15 depicts temporal standard deviation (tSD) distribution in the wholebrain volume measured across all five volunteers for no motion, nose touching and deep breathing with static, FIDnav and DOCMA distortion correction. Histograms are shown on a linear scale (top of Fig.15), as well as a logarithmic scale (bottom of Fig.15) to show the differences in high tSD values more clearly.
  • the tSD distribution for the “no motion” case is shown as a dashed line in Fig.15 for comparison.
  • median tSNR increased from 66 to 71 with FIDnav correction and 67 with DOCMA correction.
  • tSNR distributions for static distortion correction were shifted towards lower values.
  • median tSNR increased from 55 to 64 with FID-navigated correction and 63 with DOCMA-based correction
  • tSNR increased from 52 to 61 and 60 for FIDnav and DOCMA corrections, respectively.
  • the tSNR distributions for each subject are quantified in the table of Fig.16. The table provides a summary of the mean and standard deviation of the tSNR distribution for each subject and scan. FIDnav-based correction reliably improved the tSNR distribution across all subjects and paradigms, while DOCMA results were more variable.
  • FIDnav- based correction improved the tSNR distribution across all paradigms and subjects.
  • the tSNR gains across all subjects, relative to static distortion correction, are shown in Figure 17. More particularly, Fig. 17 depicts boxplots showing the tSNR gains across all five volunteers for the no motion, nose touching and deep breathing paradigms for FIDnav- and DOCMA-based dynamic distortion correction, relative to static distortion correction.
  • the tSNR gain was evaluated in the wholebrain region of interest for each subject and scan. For the no motion case, FIDnavs result in a significantly higher tSNR gain, relative to DOCMA, while for both paradigms designed to induce magnetic field perturbations, FID-navigated dynamic distortion correction results in a comparable tSNR improvement.
  • Fitting spherical harmonic coefficients to the measured field maps was performed to confirm that the measured B0 field changes matched the applied shim settings and that the applied FIDnav shim updates successfully counteracted these field changes.
  • This experiment was repeated at 7T (MAGNETOM Terra; Siemens Healthcare, Er Weg, Germany) using an anthropomorphic head phantom and a 32-channel head coil.
  • 7T MAGNETOM Terra; Siemens Healthcare, Er Weg, Germany
  • Maps of ⁇ B0 relative to the reference shim setting were generated from the multi-echo GRE data with and without FID-navigated field control.
  • the brain region was segmented on the reference (no motion) scan acquired in each subject using FSL’s BET tool (FMRIB, Oxford, United Kingdom). Peak signal-to-noise ratio (PSNR) and structural similarity (SSIM) were computed within the masked brain region following rigid-body realignment of each image acquired with and without real-time field control to the reference scan. Differences in PSNR and SSIM were assessed using two-sided paired Wilcoxon rank sum tests with 0.05. We conducted several additional volunteer experiments to further evaluate the performance of the proposed approach. One subject was rescanned using a 32-channel head coil to test the performance of FID-navigated shimming across different coil configurations.
  • phase stabilization computes the global phase difference between two non-phase encoded echoes; one acquired at the beginning of the scan and the other acquired following each image echo.
  • subject-specific calibration data were acquired as described above.
  • the plots show zeroth- and first-order shim coefficients measured from the GRE field maps that were acquired simultaneously with the FIDnavs, with and without real- time field control.
  • the fitted shim coefficients reflect the applied values, which varied between ⁇ 20 Hz and ⁇ 20 ⁇ T/m.
  • residual absolute zeroth- and first-order changes across acquisitions were 0.20 ⁇ 0.10 Hz, 0.16 ⁇ 0.05 ⁇ T/m (X) and 0.07 ⁇ 0.05 ⁇ T/m (Y) at 3T.
  • Residual shim coefficients at 7T were 1.01 ⁇ 0.59 Hz, 0.09 ⁇ 0.04 ⁇ T/m (X) and 0.16 ⁇ 0.07 ⁇ T/m (Y).
  • the performance of FID- navigated shimming in sagittal orientation is shown in Fig.19. More particularly, Fig.19 depicts first-order shim coefficients and frequency offsets measured from sagittal GRE field maps during applied Y and Z gradient offsets and central frequency shifts across acquisitions with FID-navigated shimming switched off and on at 3T. Applying B0 shim measurements from FIDnavs using the linear calibration model accurately compensates for controlled field offsets in the sagittal plane.
  • Fig.20 depicts zeroth- and first-order B 0 shim offsets measured in a phantom during application of shim changes of 10 ⁇ T/m and 10 Hz, relative to reference shim settings and in the presence of additional second-order shim changes, with FID- navigated field control switched off and on.
  • FIG. 21A shows measured ⁇ B 0 field maps in a representative volunteer with applied gradient offsets of 10 ⁇ T/m in X and Y, and a shift in the center frequency of 10 Hz, with FID- navigated real-time field control switched off and on.
  • Figure 21B summarizes the fitted shim coefficients across all slices with and without FID-navigated field control for nine volunteers scanned with applied shim offsets.
  • the residual measured shim coefficients across all volunteers following application of FID-navigated field control were 0.26 ⁇ 0.13 ⁇ T/m (X), 0.17 ⁇ 0.39 ⁇ T/m (Y), and 0.99 ⁇ 0.85 Hz, respectively.
  • the measured latency between FIDnav measurements and application of B 0 shim updates was ⁇ 12 ms.
  • applying ⁇ B 0 shim measurements computed from FIDnavs accurately compensates for controlled linear field offsets in volunteers. Impact on T 2 *-weighted imaging at 3T Ten volunteers were scanned with the FID-navigated T 2 *-weighted GRE protocol at 3T while performing tasks designed to dynamically modulate the B0 field.
  • FIG. 22 Examples of the zeroth- and first-order shim coefficients measured from FIDnavs in a single slice are shown in Figure 22 for representative volunteers during continuous deep breathing (A) and nose touching (B), with and without application of FID-navigated real-time shimming. More particularly, Fig.22 depicts zeroth- and first-order shim coefficients measured using FIDnavs for a single slice in a representative volunteer performing deep breathing (A) and nose touching (B) at 3T with real- time shimming switched off and on. Performing real-time field control based on FIDnav measurements reduces spatiotemporal field fluctuations. Both tasks introduced strong global B0 fluctuations, which were compensated by FID-navigated shimming.
  • FIG. 23 depicts boxplots showing standard deviation of zeroth- (A) and first-order (B) B0 field fluctuations measured across all volunteers scanned at 3T during deep breathing and nose touching with FID-navigated shimming switched off and on.
  • Real-time field control substantially reduces spatiotemporal field fluctuations across all volunteers.
  • the standard deviation of zeroth- and first-order field fluctuations during deep breathing was reduced by 40.1% ⁇ 16.1% and 13.3% ⁇ 29.9%, respectively.
  • FIG. 24 depicts axial T2*-weighted GRE images acquired at 3T in two volunteers during continuous deep breathing with FID-navigated shimming switched off and on.
  • FID-navigated shimming visibly reduces ghosting and signal modulation artifacts induced by variable susceptibility changes during deep breathing. ghosting, blurring and signal dropout artifacts are visibly reduced with application of FID-navigated shimming.
  • Figure 25 shows the ability of FID-navigated real-time shimming to compensate for spatiotemporal field fluctuations induced by nose touching in two volunteers.
  • Fig.25 depicts axial T 2 *-weighted GRE images acquired at 3T in two volunteers repeatedly performing a nose touching action inside the scanner with FID-navigated shimming switched off and on.
  • Reference T2*-weighted image acquired with no motion (FID- navigated shimming OFF) and difference images relative to this reference are shown for comparison.
  • FID-navigated shimming successfully compensates for artifacts arising from the spatiotemporal field changes induced by arm motion. Nose touching induced strong ghosting and signal loss in the frontal brain regions, which were considerably improved with application of real-time shimming. Image quality improvements were found in 9 out of 10 subjects for both deep breathing and nose touching tasks.
  • FIG.26 Boxplots of PSNR and SSIM across all volunteers are shown in Figure 26. More particularly, Fig.26 depicts boxplots showing peak signal-to-noise ratio (PSNR; A) and structural similarity (SSIM; B) relative to the reference T 2 *-weighted image for deep breathing and nose touching scans across all ten volunteers. FID-navigated real- time shimming significantly improved image quality across all volunteers. Across all volunteers, PSNR increased 5.84% ⁇ 5.05% for deep breathing and 9.29% ⁇ 6.00% for nose touching (P ⁇ 0.05). SSIM increased 7.31% ⁇ 14.10% for deep breathing and 9.12% ⁇ 6.10% for nose touching (P ⁇ 0.05).
  • Figure 27 shows an example of the impact of FID-navigated shimming during deep breathing in a T 2 *-weighted image acquired in the sagittal plane. More particularly, Fig.27 depicts sagittal T2*-weighted images acquired at 3T in a volunteer performing continuous deep breathing during the scan. A reference scan acquired during normal breathing is shown for comparison. Strong ghosting artifacts are induced by respiration in the inferior brain region, which are considerably improved with first-order FID-navigated shimming. The strong ghosting artifacts are observed near the brain stem, which is expected since field fluctuations are greatest in more inferior brain regions closest to the lungs. The amount and severity of artifacts are considerably improved with FID-navigated real-time field control.
  • FIG.28 A comparison of improvement in T2*-weighted image quality at 3T using a 32-channel and a 64-channel head coil is shown in Fig.28. More particularly, Fig.28 depicts a comparison of FID-navigated shimming performance in a volunteer scanned at 3T with a 32-channel and a 64-channel head coil. The volunteer repeated the same nose touching action with FID-navigated shimming switched off and on. A reference scan with no motion is shown for comparison. Performance of the proposed FID-navigated shimming approach was similar for both coil configurations.
  • Fig.29 depicts, in the top row, axial T2*-weighted images acquired in a volunteer during nose touching at 3T without any correction (A); with FID-navigated zeroth-order shimming (B); and with the proposed FID-navigated first-order shimming (C).
  • FIG.29 also depicts, in the bottom row, sagittal T2*-weighted images acquired in the same volunteer during deep breathing without any correction (E); with zeroth-order shimming (F); with the proposed first-order shimming (G); reference image acquired during normal breathing (H).
  • FID-navigated first-order shimming yielded substantially improved image quality compared to both no correction and zeroth-order only shimming, which is expected since breathing-induced field fluctuations vary along the head-foot axis.
  • Image quality and corresponding quantitative metrics are substantially improved with first-order correction for nose touching in the axial plane, compared to both zeroth-order and no correction, which is expected since arm motion induces substantial spatiotemporal field variations, particularly in frontal brain regions.
  • first- order shimming yielded improved performance over zeroth-order correction for deep breathing in the sagittal plane, where linear field variations are expected along the head-foot axis.
  • FIG. 30 depicts measured shim coefficients in a volunteer performing the nose touching task at 7T (A) and with FID-navigated shimming (B).
  • a reference scan acquired without any intentional motion (with PS correction) is shown for comparison (F).
  • Fig.31 depicts sagittal T 2 *-weighted images acquired during deep breathing at 7T without any phase stabilization (PS) correction applied (A); with retrospective PS (B); with FID-navigated shimming (C); and a reference scan acquired during normal breathing with retrospective PS correction (D). Zoomed in sections (E-H) show that retrospective PS and prospective FID-navigated shimming exhibit comparable performance for improving image quality.
  • Example Processing Techniques Techniques operating according to the principles described herein may be implemented in any suitable manner.
  • Such computer-executable instructions may be written using any of a number of suitable programming languages and/or programming or scripting tools, and also may be compiled as executable machine language code or intermediate code that is executed on a framework or virtual machine.
  • these computer-executable instructions may be implemented in any suitable manner, including as a number of functional facilities, each providing one or more operations to complete execution of algorithms operating according to these techniques.
  • a “functional facility,” however instantiated, is a structural component of a computer system that, when integrated with and executed by one or more computers, causes the one or more computers to perform a specific operational role.
  • a functional facility may be a portion of or an entire software element.
  • a functional facility may be implemented as a function of a process, or as a discrete process, or as any other suitable unit of processing. If techniques described herein are implemented as multiple functional facilities, each functional facility may be implemented in its own way; all need not be implemented the same way. Additionally, these functional facilities may be executed in parallel and/or serially, as appropriate, and may pass information between one another using a shared memory on the computer(s) on which they are executing, using a message passing protocol, or in any other suitable way.
  • functional facilities include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types.
  • the functionality of the functional facilities may be combined or distributed as desired in the systems in which they operate.
  • one or more functional facilities carrying out techniques herein may together form a complete software package. These functional facilities may, in alternative embodiments, be adapted to interact with other, unrelated functional facilities and/or processes, to implement a software program application, for example as a software program application such as a signal analysis facility.
  • Some exemplary functional facilities have been described herein for carrying out one or more tasks. It should be appreciated, though, that the functional facilities and division of tasks described is merely illustrative of the type of functional facilities that may implement the exemplary techniques described herein, and that embodiments are not limited to being implemented in any specific number, division, or type of functional facilities. In some implementations, all functionality may be implemented in a single functional facility.
  • Computer-executable instructions implementing the techniques described herein may, in some embodiments, be encoded on one or more computer-readable media to provide functionality to the media.
  • Computer-readable media include magnetic media such as a hard disk drive, optical media such as a Compact Disk (CD) or a Digital Versatile Disk (DVD), a persistent or non- persistent solid-state memory (e.g., Flash memory, Magnetic RAM, etc.), or any other suitable storage media.
  • Such a computer-readable medium may be implemented in any suitable manner, including as computer-readable storage media 3206 of FIG.32 described below (i.e., as a portion of a computing device 3200) or as a stand-alone, separate storage medium.
  • “computer-readable media” also called “computer-readable storage media” refers to tangible storage media. Tangible storage media are non-transitory and have at least one physical, structural component. In a “computer-readable medium,” as used herein, at least one physical, structural component has at least one physical property that may be altered in some way during a process of creating the medium with embedded information, a process of recording information thereon, or any other process of encoding the medium with information.
  • a magnetization state of a portion of a physical structure of a computer-readable medium may be altered during a recording process.
  • these instructions may be executed on one or more suitable computing device(s) operating in any suitable computer system, including the exemplary computer system of FIG.1, or one or more computing devices (or one or more processors of one or more computing devices) may be programmed to execute the computer-executable instructions.
  • a computing device or processor may be programmed to execute instructions when the instructions are stored in a manner accessible to the computing device or processor, such as in a data store (e.g., an on-chip cache or instruction register, a computer-readable storage medium accessible via a bus, a computer-readable storage medium accessible via one or more networks and accessible by the device/processor, etc.).
  • a data store e.g., an on-chip cache or instruction register, a computer-readable storage medium accessible via a bus, a computer-readable storage medium accessible via one or more networks and accessible by the device/processor, etc.
  • FIG.32 illustrates one exemplary implementation of a computing device in the form of a computing device 3200 that may be used in a system implementing techniques described herein, although others are possible.
  • FIG.32 is intended neither to be a depiction of necessary components for a computing device to operate as a device for controlling generation of an image that would result from an MRI scan in accordance with the principles described herein, nor a comprehensive depiction.
  • Computing device 3200 may comprise at least one processor 3202 (e.g., a computer hardware processor), a network adapter 3204, and computer-readable storage media 3206.
  • Computing device 3200 may be, for example, a desktop or laptop personal computer, a personal digital assistant (PDA), a smart mobile phone, or any other suitable computing device.
  • Network adapter 3204 may be any suitable hardware and/or software to enable the computing device 3200 to communicate wired and/or wirelessly with any other suitable computing device over any suitable computing network.
  • the computing network may include wireless access points, switches, routers, gateways, and/or other networking equipment as well as any suitable wired and/or wireless communication medium or media for exchanging data between two or more computers, including the Internet.
  • Computer-readable media 3206 may be adapted to store data to be processed and/or instructions to be executed by processor 3202. Processor 3202 enables processing of data and execution of instructions. The data and instructions may be stored on the computer-readable storage media 3206. The data and instructions stored on computer-readable storage media 3206 may comprise computer-executable instructions implementing techniques which operate according to the principles described herein. In the example of FIG.32, computer-readable storage media 3206 stores computer-executable instructions implementing various facilities and storing various information as described above.
  • Computer-readable storage media 3206 may store signal analysis facility 3208 configured to perform techniques described herein.
  • Media 3206 may also store an MR image generator 3210 to generate one or more images from signals captured during an MRI scan.
  • a computing device may additionally have one or more components and peripherals, including input and output devices. These devices can be used, among other things, to present a user interface. Examples of output devices that can be used to provide a user interface include printers or display screens for visual presentation of output and speakers or other sound generating devices for audible presentation of output. Examples of input devices that can be used for a user interface include keyboards, and pointing devices, such as mice, touch pads, and digitizing tablets.
  • a computing device may receive input information through speech recognition or in other audible format.
  • Embodiments have been described where the techniques are implemented in circuitry and/or computer-executable instructions. It should be appreciated that some embodiments may be in the form of a method, of which at least one example has been provided. The acts performed as part of the method may be ordered in any suitable way. Accordingly, embodiments may be constructed in which acts are performed in an order different than illustrated, which may include performing some acts simultaneously, even though shown as sequential acts in illustrative embodiments.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Signal Processing (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Selon divers modes de réalisation, la présente invention concerne des techniques destinées à commander la création d'une résonance magnétique (RM) à partir d'un balayage d'imagerie par résonance magnétique (IRM) d'un sujet à l'aide d'un appareil d'IRM, en temps réel, la commande dynamique de l'appareil d'IRM pendant le balayage IRM et/ou du traitement pendant le balayage IRM de signaux RM acquis en tant que partie de génération de données RM pour la création de l'image RM. Plus particulièrement, selon certains modes de réalisation, la présente invention concerne l'identification de changements qui se produisent dans un champ magnétique à l'intérieur d'un appareil d'IRM au cours d'un balayage IRM. En réponse à des changements dans le champ magnétique, un appareil d'IRM peut être réglé pour changer une manière dans laquelle des données RM supplémentaires sont capturées et/ou un traitement de signaux RM capturés (qui doivent être utilisés pour générer une image RM) peut être réglé.
PCT/US2022/040917 2021-08-30 2022-08-19 Correction de distorsion dynamique pour irm à l'aide de navigateurs fid WO2023034044A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US202163238640P 2021-08-30 2021-08-30
US63/238,640 2021-08-30

Publications (1)

Publication Number Publication Date
WO2023034044A1 true WO2023034044A1 (fr) 2023-03-09

Family

ID=85412814

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2022/040917 WO2023034044A1 (fr) 2021-08-30 2022-08-19 Correction de distorsion dynamique pour irm à l'aide de navigateurs fid

Country Status (1)

Country Link
WO (1) WO2023034044A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120249138A1 (en) * 2011-03-28 2012-10-04 Josef Pfeuffer Magnetic resonance method and system to correct phase information in mr images
US20130076356A1 (en) * 2011-09-26 2013-03-28 Vladimir Jellus Magnetic resonance imaging method and apparatus to correct distortions due to inhomogeneities of the basic magnetic field
US20130229177A1 (en) * 2011-08-30 2013-09-05 Roland Bammer Apparatus for Real-Time Phase Correction for Diffusion-Weighted Magnetic Resonance Imaging Using Adaptive RF Pulses

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120249138A1 (en) * 2011-03-28 2012-10-04 Josef Pfeuffer Magnetic resonance method and system to correct phase information in mr images
US20130229177A1 (en) * 2011-08-30 2013-09-05 Roland Bammer Apparatus for Real-Time Phase Correction for Diffusion-Weighted Magnetic Resonance Imaging Using Adaptive RF Pulses
US20130076356A1 (en) * 2011-09-26 2013-03-28 Vladimir Jellus Magnetic resonance imaging method and apparatus to correct distortions due to inhomogeneities of the basic magnetic field

Similar Documents

Publication Publication Date Title
Huang et al. Body MR imaging: artifacts, k-space, and solutions
Finsterbusch et al. Single, slice-specific z-shim gradient pulses improve T2*-weighted imaging of the spinal cord
JP6462306B2 (ja) 医用画像処理装置及び磁気共鳴イメージング装置
US9317917B2 (en) Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
US9983284B2 (en) MRI with dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
US10031201B2 (en) Method and apparatus for magnetic resonance imaging
US9330457B2 (en) Device, method, and program for analyzing a magnetic resonance image using phase difference distribution
Xu et al. Prospective and retrospective high order eddy current mitigation for diffusion weighted echo planar imaging
US11486953B2 (en) Phase estimation for retrospective motion correction
Gretsch et al. Investigating the accuracy of FatNav‐derived estimates of temporal B0 changes and their application to retrospective correction of high‐resolution 3 D GRE of the human brain at 7 T
Wallace et al. Dynamic distortion correction for functional MRI using FID navigators
Meineke et al. Data consistency‐driven determination of‐fluctuations in gradient‐echo MRI
Brackenier et al. Data‐driven motion‐corrected brain MRI incorporating pose‐dependent B0 fields
US10690741B2 (en) Method and systems for reducing artifacts in magnetic resonance imaging
Andersson Geometric distortions in diffusion MRI
US11402454B2 (en) Correction of distorted diffusion-weighted magnetic resonance image data
WO2023034044A1 (fr) Correction de distorsion dynamique pour irm à l'aide de navigateurs fid
US20220057467A1 (en) Epi mr imaging with distortion correction
US20160124065A1 (en) Method and apparatus for correction of magnetic resonance image recordings with the use of a converted field map
US11061097B2 (en) Data driven correction of phase depending artefacts in a magnetic resonance imaging system
US11226385B2 (en) Dixon type water/fat separation MR imaging with improved fat shift correction
Wallace et al. Rapid motion estimation and correction using self‐encoded FID navigators in 3D radial MRI
US20230293039A1 (en) Methods for acquiring a magnetic resonance image dataset and for generating a motion-corrected image dataset
EP4343356A1 (fr) Imagerie par résonance magnétique avec cartographie eau/graisse/b0
Wallace et al. Navigators without gradients

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22865300

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE