US20120271583A1 - System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging - Google Patents
System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging Download PDFInfo
- Publication number
- US20120271583A1 US20120271583A1 US13/397,813 US201213397813A US2012271583A1 US 20120271583 A1 US20120271583 A1 US 20120271583A1 US 201213397813 A US201213397813 A US 201213397813A US 2012271583 A1 US2012271583 A1 US 2012271583A1
- Authority
- US
- United States
- Prior art keywords
- high order
- computer
- current
- basis
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient 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/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56341—Diffusion imaging
Definitions
- Embodiments of the invention relate generally to magnetic resonance (MR) imaging and, more particularly, to correcting high order eddy-current-induced distortion in diffusion-weighted echo planar imaging.
- MR magnetic resonance
- polarizing field B 0 When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field B 0 ), the individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but precess about it in random order at their characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic field (excitation field B i ) which is in the x-y plane and which is near the Larmor frequency, the net aligned moment, or “longitudinal magnetization”, M z , may be rotated, or “tipped”, into the x-y plane to produce a net transverse magnetic moment M t . A signal is emitted by the excited spins after the excitation signal B 1 is terminated and this signal may be received and processed to form an image.
- excitation field B i which is in the x-y plane and which is near the Larmor frequency
- magnetic field gradients (G x , G y , and G z ) are employed.
- the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used.
- the resulting set of received NMR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
- DW-EPI Diffusion-Weighted Echo Planar Imaging
- Dual spin echo (also called twice refocused) DW-EPI has been proposed to provide a certain level of inherent eddy current cancellation, but with a significant increase in echo time and decrease in signal-to-noise ratio (SNR).
- SNR signal-to-noise ratio
- a typical dual spin echo protocol may generate about half as much SNR as the corresponding single spin echo (also called Stejkal-Tanner sequence) protocol on liver imaging at 3T.
- increasing NEX is not an option to increase SNR because of the associated increase in scan time. Therefore, it is desirable to keep single spin-echo while reducing the resulting distortion in practice.
- an MRI apparatus comprises a magnetic resonance imaging (MRI) system having a plurality of gradient coils positioned about a bore of a magnet, and an RF transceiver system and an RF switch controlled by a pulse module to transmit RF signals to an RF coil assembly to acquire MR images.
- the MRI apparatus also comprises a computer programmed to acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy-current-generated magnetic field error of an imaging system.
- the computer is also programmed to process the calibration data to generate a plurality of basis coefficients and a plurality of time constants and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence.
- the computer is further programmed to execute a diffusion-weighted imaging scan that comprises application of a DW-EPI pulse sequence to acquire MR data from an imaging subject and reconstruction of an image based on the acquired MR data.
- the computer is also programmed to apply high order eddy-current-generated magnetic field error correction during image reconstruction configured to reduce high order eddy-current-induced distortion in the reconstructed image.
- a method for correcting high order eddy-current-induced distortions in diffusion-weighted echo planar imaging comprises acquiring calibration data from a calibration scan, the calibration data configured to characterize high order eddy currents of an imaging system, processing the calibration data to generate a plurality of basis coefficients and a plurality of time constants, and calculating a plurality of basis correction coefficients based on the plurality of basis coefficients and based on the plurality of time constants.
- the method also comprises applying a DW-EPI pulse sequence to acquire MR data from an imaging subject, reconstructing an image based on the acquired MR data, and wherein the reconstructing the image comprises applying high order eddy-current-generated magnetic field error correction configured to reduce high order eddy-current-induced distortion in the image.
- a non-transitory computer readable medium having stored thereon a computer program comprising a set of instructions, which, when executed by a computer, causes the computer to acquire calibration data from a calibration scan configured to characterize high order eddy current generated magnetic field error of an imaging system and to process the calibration data.
- the set of instructions also causes the computer to generate a plurality of basis coefficients and a plurality of time constants based on the processed calibration data and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a DW-EPI pulse sequence.
- the set of instructions also causes of the computer to apply the DW-EPI pulse sequence to acquire MR data from an imaging subject, to reconstruct an image based on the acquired MR data, and to apply high order eddy-current-generated magnetic field error correction during reconstruction of the image configured to reduce high order eddy current induced distortion in the reconstructed image.
- FIG. 1 is a schematic block diagram of an MR imaging system for use with embodiments of the invention.
- FIG. 2 is a pulse sequence diagram showing ideal gradient and RF waveforms for a single spin echo diffusion-weighted EPI scan.
- FIG. 3 is a flowchart showing an HOEC correction technique according to an embodiment of the invention.
- FIG. 4 is a flowchart showing an HOEC calibration data processing algorithm that may be used in the HOEC correction technique of FIG. 3 according to an embodiment of the invention.
- FIG. 5 is a flowchart showing another HOEC data processing algorithm that may be used in the HOEC correction technique of FIG. 3 according to an embodiment of the invention.
- FIG. 6 is a flowchart showing an algorithm for calculating HOEC terms that may be used in the HOEC correction technique of FIG. 3 according to an embodiment of the invention.
- FIG. 7 is a flowchart showing an algorithm for retrospective compensation of HOEC terms that may be used in the HOEC correction technique of FIG. 3 according to an embodiment of the invention.
- FIG. 1 the major components of a magnetic resonance imaging (MRI) system 10 incorporating an embodiment of the invention are shown.
- the operation of the system is controlled for certain functions from an operator console 12 which in this example includes a keyboard or other input device 13 , a control panel 14 , and a display screen 16 .
- the console 12 communicates through a link 18 with a separate computer system 20 that enables an operator to control the production and display of images on the display screen 16 .
- the computer system 20 includes a number of modules which communicate with each other through a backplane 20 a. These modules include an image processor module 22 , a CPU module 24 and a memory module 26 , known in the art as a frame buffer for storing image data arrays.
- the computer system 20 communicates with a separate system control 32 through a high speed serial link 34 .
- the input device 13 can include a mouse, joystick, keyboard, track ball, touch activated screen, light wand, voice control, card reader, push-button, or any similar or equivalent input device, and may be used for interactive geometry prescription.
- the system control 32 includes a set of modules connected together by a backplane 32 a. These include a CPU module 36 and a pulse generator module 38 which connects to the operator console 12 through a serial link 40 . It is through link 40 that the system control 32 receives commands from the operator to indicate the scan sequence that is to be performed.
- the pulse generator module 38 operates the system components to carry out the desired scan sequence and produces data which indicates the timing, strength and shape of the RF pulses produced, and the timing and length of the data acquisition window.
- the pulse generator module 38 connects to a set of gradient amplifiers 42 , to indicate the timing and shape of the gradient pulses that are produced during the scan.
- the pulse generator module 38 can also receive patient data from a physiological acquisition controller 44 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes attached to the patient. And finally, the pulse generator module 38 connects to a scan room interface circuit 46 which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 46 that a patient positioning system 48 receives commands to move the patient to the desired position for the scan.
- the gradient waveforms produced by the pulse generator module 38 are applied to the gradient amplifier system 42 having Gx, Gy, and Gz amplifiers.
- Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly generally designated 50 to produce the magnetic field gradients used for spatially encoding acquired signals.
- the gradient coil assembly 50 forms part of a resonance assembly 52 which includes a polarizing magnet 54 and a whole-body RF coil 56 .
- a transceiver module 58 in the system control 32 produces pulses which are amplified by an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive switch 62 .
- the resulting signals emitted by the excited nuclei in the patient may be sensed by the same RF coil 56 and coupled through the transmit/receive switch 62 to a preamplifier 64 .
- the amplified MR signals are demodulated, filtered, and digitized in the receiver section of the transceiver 58 .
- the transmit/receive switch 62 is controlled by a signal from the pulse generator module 38 to electrically connect the RF amplifier 60 to the coil 56 during the transmit mode and to connect the preamplifier 64 to the coil 56 during the receive mode.
- the transmit/receive switch 62 can also enable a separate RF coil (for example, a surface coil) to be used in either the transmit or receive mode.
- the MR signals picked up by the RF coil 56 are digitized by the transceiver module 58 and transferred to a memory module 66 in the system control 32 .
- a scan is complete when an array of raw k-space data has been acquired in the memory module 66 .
- This raw k-space data is rearranged into separate k-space data arrays for each image to be reconstructed, and each of these is input to an array processor 68 which operates to Fourier transform the data into an array of image data.
- This image data is conveyed through the serial link 34 to the computer system 20 where it is stored in memory.
- this image data may be archived in long term storage or it may be further processed by the image processor 22 and conveyed to the operator console 12 and presented on the display 16 .
- Embodiments of the invention correct high order eddy current (HOEC) induced diffusion gradient direction dependent distortion in Diffusion-Weighted Echo Planar Imaging (DW-EPI).
- “high order” means spatial order higher than or equal to 2, as compared to linear and constant orders, which mean order 1 and 0, respectively.
- embodiments of the invention can also be applied to linear and constant orders as well.
- the distortions due to HOEC may strongly affect the acquired data, especially for body applications where large FOV is used and large slice coverage is desired.
- a system calibration may be used to characterize eddy currents of spatial orders less than or equal to P, where P is usually 3 to 5 , but can be any order in general.
- DW-EPI pulse sequence with correction amplitudes for gradients in readout, phase encoding, and slice axes, and receiver frequency on a per-slice basis may be played to prospectively remove the effect of the quasi-linear HOEC terms.
- a DW-EPI pulse sequence 70 is shown including 90° and 180° RF pulses 72 and 74 , respectively.
- RF pulses 72 , 74 can be transmitted by RF coil 56 to generate an echo signal 76 , which can be encoded with spatial information.
- Echo signal 76 can also be received by coil 56 or by another coil, such as a surface coil, for use in reconstructing an image.
- the sequence shown in FIG. 2 further includes read-out, phase-encoding, and slice-selection gradients G ro , G pe , and G sl , respectively.
- Readout gradient G ro comprises a pre-phasing pulse 78 and read-out pulses 80 .
- phase-encoding gradient G pe comprises a pre-phasing pulse 82 and phase-encoding pulses 84 .
- Slice-selection gradient G sl comprises slice selection pulses 86 for the 90° RF pulse 72 and 88 for the 180° RF pulse 74 , as well as 86 a for slice refocusing.
- Diffusion-weighting gradient G d used in a diffusion-weighted EPI scan.
- Diffusion-weighting gradient G d comprises two equivalent trapezoidal pulses 90 and 92 , placed at either side of the 180° RF pulse 74 .
- G d can contain components on all three logical axes.
- the read-out, phase-encoding, and slice axes i.e., the logical axes
- the components of G d in u, v, and w axes are denoted as G u , G v , and G w , respectively.
- the physical axes are denoted as x, y, and z.
- FIG. 3 shows an HOEC correction technique 94 according to an embodiment of the invention.
- Technique 94 begins at block 96 by performing an HOEC calibration to characterize the eddy-current-generated magnetic field error of a particular imaging system such as MRI system 10 of FIG. 1 .
- the calibration can either be phantom-based or local-pickup-coil-based.
- phantom-based method gradient echo images are often collected at multiple time points following an eddy current generating gradient.
- the local-pickup-coil-based method a number of local coils, each with a small sample, are used to obtain free induction decay signals at their respective spatial locations.
- HOEC calibration can be done as frequently as needed, but is in general only needed once per system installation or when there are significant system hardware (e.g., gradient coil) changes.
- a 4-dimensional (3D in space and 1D in time) eddy current field data set is generated after HOEC calibration scan.
- FIGS. 4 and 5 illustrate HOEC data processing algorithms 116 , 118 that may be used in block 98 of the HOEC correction technique 94 of FIG. 3 according to an embodiment of the invention.
- HOEC data processing algorithm 116 is used when the HOEC calibration of block 96 of technique 94 is performed using the phantom-based method as described above.
- Data from the HOEC calibration scan are acquired at block 120 .
- a 3D phase unwrapping is applied to the phase angle of the data, and the phase angle is scaled at block 124 by a factor proportional to the echo time to yield a magnetic field data set.
- Each time point of the magnetic field data set is then spatially fitted at block 126 to polynomial or spherical harmonic bases of order up to P to generate basis coefficients, where P is usually 3 to 5, but can be any order in general.
- magnitude weights or masks can be optionally used during the spatial fitting.
- the resulting basis coefficients are then temporal fitted along the time axis at block 128 using a single-exponential or multi-exponential model.
- B n (x, y, z) are assumed to be polynomial bases for convenience of discussion.
- HOEC data processing algorithm 118 is used when the HOEC calibration of block 96 of technique 94 is performed using the local-pickup-coil-based method as described above.
- Data from the HOEC calibration scan are acquired at block 132 .
- time derivatives are taken on the phase angle of the data to obtain the magnetic field offset at a coil location.
- Each time point of the magnetic field offset is then spatially fitted at block 136 to polynomial or spherical harmonic bases of order up to P to generate basis coefficients, where P is usually 3 to 5, but can be any order in general.
- magnitude weights or masks can be optionally used during the spatial fitting.
- the resulting basis coefficients are then temporal fitted along the time axis at block 138 , and a set of ( ⁇ mn , ⁇ mn ) pairs is generated.
- the ( ⁇ mn , ⁇ mn ) pairs are saved or stored on the host computer of the scanner for future use at block 140 .
- HOEC data processing performed in HOEC data processing algorithms 116 and 118 only needs to be done once per system installation. However, algorithms 116 and 118 can be performed as frequently as needed.
- DW-EPI protocol dependent HOEC terms are calculated at block 100 .
- block 100 can handle arbitrary imaging planes.
- an arbitrary imaging plane means straight axial, coronal, or sagittal scan plane, as well as any oblique plane.
- an algorithm 142 for calculating the HOEC terms for block 100 of FIG. 3 is shown.
- the diffusion gradient components G u , G v , G w of the DW-EPI pulse sequence to be used are obtained.
- logical gradients G u , G v , G w are converted to physical components G x , G y , G z by applying the 3 ⁇ 3 axis rotation matrix R:
- the ( ⁇ mn , ⁇ mn ) pairs determined via block 98 from the HOEC calibration scan at block 96 of technique 94 are obtained, and the pulse sequence type and sequence timing dependent constant, ⁇ nm , is calculated at block 148 .
- Derivation of ⁇ mn can either be analytical or using convolution. Although all gradient waveforms can be included to obtain ⁇ nm , contributions from the diffusion gradients are often dominant, which allow for simplified analysis to obtain ⁇ mn . For example, when single spin echo DW-EPI is used such as that shown in FIG. 2 , it can be derived that:
- ⁇ mn ( 1 - ⁇ t 1 ⁇ mn ) ⁇ ( 1 - ⁇ t 2 ⁇ mn ) ⁇ ( 1 + ⁇ t 3 ⁇ mn ) t 1 , ( Eqn . ⁇ 1 )
- t 1 , t 2 , and t 3 are sequence timing related constants shown in FIG. 2 . Note that ⁇ nm for other pulse sequences such as the dual spin echo or stimulated echo DW-EPI can also be determined analytically.
- the HOEC physical basis coefficients d n (t) for the nth basis function ⁇ n (x, y, z) at time t after the last diffusion gradient are calculated based on the equation:
- G m is the x, y, or z component of the diffusion gradient amplitude.
- the transpose of another rotation matrix is applied to d 1 (t), d 2 (t), . . . , d N (t) to convert them into HOEC logical basis coefficients c 1 (t), c 2 (t), . . . , c N (t):
- F is an N ⁇ N matrix that transforms polynomial bases from logical to physical coordinates. The actual form of F depends on the polynomial order and how the basis functions are numbered. Without loss of generality, the bases are in the following order: 1, x, y, z, x 2 ,xy, xz, y 2 , yz, z 2 , x 3 , x 2 y, x 2 z, xy 2 , xyz, xz 2 , y 3 , y 2 z, yz 2 , z 3 , . . .
- F can be determined by the relationship between B n (x, y, z) and B n (u, v, w), where
- the protocol dependent HOEC terms can be compensated for in technique 94 retrospectively.
- Retrospective compensation includes applying a DW-EPI pulse sequence to acquire data from the imaging subject at block 102 . Then, during image reconstruction at block 104 , the HOEC terms are retrospectively compensated for.
- FIG. 7 illustrates a retrospective compensation algorithm 162 for block 104 according to an embodiment of the invention. Note that due to the application of F, bases are now in logical axes.
- Image distortion due to HOEC fields can be corrected using a field map based method in image reconstruction according to this embodiment.
- Algorithm 162 does distortion correction in the image domain on a per-slice basis.
- the w dependency in the notation is dropped.
- the eddy current field f(u, v).
- the value of the HOEC related field map f(x, y) can be calculated at block 164 by the equation:
- t 0 is the time point at which the HOEC field is approximated.
- the choice of t 0 can be arbitrary, but it is preferable to be chosen at the time when the k-space center data is acquired (i.e., at the echo time).
- the eddy current field is then normalized by the effective gradient in the phase encoding axis G PE to obtain the pixel shift map h(u, v) at block 166 by the equation:
- G PE 1/( ⁇ TL) with ⁇ being the gyromagnetic ratio of the nucleus of interest, T being the EPI echo spacing and L being the field of view in the phase encoding axis.
- I corrected ⁇ ( u , v ) ( 1 + ⁇ h ⁇ ( u , v ) ⁇ v ) ⁇ I distorted ⁇ ( u , v + h ⁇ ( u , v ) ) , ( Eqn . ⁇ 5 )
- the image reconstructed at block 104 may be displayed to a user or stored on an image storage database for future use at block 106 .
- compensation of HOEC-induced distortions for DW-EPI can make single spin echo DW-EPI more practical. While, single spin echo has SNR and scan time benefits over dual spin echo DWI, embodiments of the invention can also be applied to dual spin echo and other variants of DW-EPI sequences to reduce distortion. Embodiments of the invention can produce a significant leverage of whole body DWI, which may have a long scan time, low SNR, and large image distortions. Brain DWI may benefit as well, especially for high b-value, large parallel imaging factor cases where SNR can be a problem.
- a technical contribution for the disclosed method and apparatus is that it provides for a computer implemented correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging.
- the computer readable storage medium includes a plurality of components such as one or more of electronic components, hardware components, and/or computer software components. These components may include one or more computer readable storage media that generally stores instructions such as software, firmware and/or assembly language for performing one or more portions of one or more implementations or embodiments of a sequence. These computer readable storage media are generally non-transitory and/or tangible. Examples of such a computer readable storage medium include a recordable data storage medium of a computer and/or storage device.
- the computer readable storage media may employ, for example, one or more of a magnetic, electrical, optical, biological, and/or atomic data storage medium. Further, such media may take the form of, for example, floppy disks, magnetic tapes, CD-ROMs, DVD-ROMs, hard disk drives, and/or electronic memory. Other forms of non-transitory and/or tangible computer readable storage media not list may be employed with embodiments of the invention.
- Such components can be combined or divided in an implementation of a system. Further, such components may include a set and/or series of computer instructions written in or implemented with any of a number of programming languages, as will be appreciated by those skilled in the art.
- other forms of computer readable media such as a carrier wave may be employed to embody a computer data signal representing a sequence of instructions that when executed by one or more computers causes the one or more computers to perform one or more portions of one or more implementations or embodiments of a sequence.
- an MRI apparatus comprises a magnetic resonance imaging (MRI) system having a plurality of gradient coils positioned about a bore of a magnet, and an RF transceiver system and an RF switch controlled by a pulse module to transmit RF signals to an RF coil assembly to acquire MR images.
- the MRI apparatus also comprises a computer programmed to acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy-current-generated magnetic field error of an imaging system.
- the computer is also programmed to process the calibration data to generate a plurality of basis coefficients and a plurality of time constants and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence.
- the computer is further programmed to execute a diffusion-weighted imaging scan that comprises application of a DW-EPI pulse sequence to acquire MR data from an imaging subject and reconstruction of an image based on the acquired MR data.
- the computer is also programmed to apply high order eddy-current-generated magnetic field error correction during image reconstruction configured to reduce high order eddy-current-induced distortion in the reconstructed image.
- a method for correcting high order eddy-current-induced distortions in diffusion-weighted echo planar imaging comprises acquiring calibration data from a calibration scan, the calibration data configured to characterize high order eddy currents of an imaging system, processing the calibration data to generate a plurality of basis coefficients and a plurality of time constants, and calculating a plurality of basis correction coefficients based on the plurality of basis coefficients and based on the plurality of time constants.
- the method also comprises applying a DW-EPI pulse sequence to acquire MR data from an imaging subject, reconstructing an image based on the acquired MR data, and wherein the reconstructing the image comprises applying high order eddy-current-generated magnetic field error correction configured to reduce high order eddy-current-induced distortion in the image.
- a non-transitory computer readable medium having stored thereon a computer program comprising a set of instructions, which, when executed by a computer, causes the computer to acquire calibration data from a calibration scan configured to characterize high order eddy current generated magnetic field error of an imaging system and to process the calibration data.
- the set of instructions also causes the computer to generate a plurality of basis coefficients and a plurality of time constants based on the processed calibration data and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a DW-EPI pulse sequence.
- the set of instructions also causes of the computer to apply the DW-EPI pulse sequence to acquire MR data from an imaging subject, to reconstruct an image based on the acquired MR data, and to apply high order eddy-current-generated magnetic field error correction during reconstruction of the image configured to reduce high order eddy current induced distortion in the reconstructed image.
Abstract
A computer is programmed to acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy current (HOEC) generated magnetic field error of an imaging system. The computer is also programmed to process the calibration data to generate a plurality of basis coefficients and a plurality of time constants and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence. The computer is further programmed to execute a diffusion-weighted imaging scan that comprises application of a DW-EPI pulse sequence to acquire MR data from an imaging subject and reconstruction of an image based on the acquired MR data. The computer is also programmed to apply HOEC-generated magnetic field error correction during image reconstruction configured to reduce HOEC-induced distortion in the reconstructed image.
Description
- The present application is a non-provisional of, and claims priority to, U.S. Provisional Patent Application Ser. No. 61/476,936, filed Apr. 19, 2011, the disclosure of which is incorporated herein by reference.
- Embodiments of the invention relate generally to magnetic resonance (MR) imaging and, more particularly, to correcting high order eddy-current-induced distortion in diffusion-weighted echo planar imaging.
- When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field B0), the individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but precess about it in random order at their characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic field (excitation field Bi) which is in the x-y plane and which is near the Larmor frequency, the net aligned moment, or “longitudinal magnetization”, Mz, may be rotated, or “tipped”, into the x-y plane to produce a net transverse magnetic moment Mt. A signal is emitted by the excited spins after the excitation signal B1 is terminated and this signal may be received and processed to form an image.
- When utilizing these signals to produce images, magnetic field gradients (Gx, Gy, and Gz) are employed. Typically, the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used. The resulting set of received NMR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
- It is well known that Diffusion-Weighted Echo Planar Imaging (DW-EPI) often suffers from diffusion encoding direction dependent distortions due to diffusion gradient generated eddy current field. These distortions, if not corrected, can lead to mis-registration among DW images of different directions and inaccuracies in any post processing operations involving DW image combination. Dual spin echo (also called twice refocused) DW-EPI has been proposed to provide a certain level of inherent eddy current cancellation, but with a significant increase in echo time and decrease in signal-to-noise ratio (SNR). For example, a typical dual spin echo protocol may generate about half as much SNR as the corresponding single spin echo (also called Stejkal-Tanner sequence) protocol on liver imaging at 3T. In many cases (e.g., whole body DW-EPI), increasing NEX is not an option to increase SNR because of the associated increase in scan time. Therefore, it is desirable to keep single spin-echo while reducing the resulting distortion in practice.
- Conventional distortion correction methods have focused on correcting only the linear and constant eddy currents (also called B0 eddy currents), either by pre-emphasis or by explicitly modifying gradient waveforms and receive frequency. However, uncompensated eddy currents of high spatial order due to gradient coil leakage field, or simply high order eddy currents (HOEC), can also be significant with the desire for increased b values and the increase of gradient amplitude and slew rate in modern MR scanners. Because of the high spatial order, distortions generated by the magnetic fields created by these eddy currents are not only diffusion gradient direction dependent, but also slice dependent.
- It would therefore be desirable to have a system and method capable of correcting distortion due to HOEC in DW-EPI.
- According to an aspect of the invention, an MRI apparatus comprises a magnetic resonance imaging (MRI) system having a plurality of gradient coils positioned about a bore of a magnet, and an RF transceiver system and an RF switch controlled by a pulse module to transmit RF signals to an RF coil assembly to acquire MR images. The MRI apparatus also comprises a computer programmed to acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy-current-generated magnetic field error of an imaging system. The computer is also programmed to process the calibration data to generate a plurality of basis coefficients and a plurality of time constants and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence. The computer is further programmed to execute a diffusion-weighted imaging scan that comprises application of a DW-EPI pulse sequence to acquire MR data from an imaging subject and reconstruction of an image based on the acquired MR data. The computer is also programmed to apply high order eddy-current-generated magnetic field error correction during image reconstruction configured to reduce high order eddy-current-induced distortion in the reconstructed image.
- According to another aspect of the invention, a method for correcting high order eddy-current-induced distortions in diffusion-weighted echo planar imaging (DW-EPI) comprises acquiring calibration data from a calibration scan, the calibration data configured to characterize high order eddy currents of an imaging system, processing the calibration data to generate a plurality of basis coefficients and a plurality of time constants, and calculating a plurality of basis correction coefficients based on the plurality of basis coefficients and based on the plurality of time constants. The method also comprises applying a DW-EPI pulse sequence to acquire MR data from an imaging subject, reconstructing an image based on the acquired MR data, and wherein the reconstructing the image comprises applying high order eddy-current-generated magnetic field error correction configured to reduce high order eddy-current-induced distortion in the image.
- According to yet another aspect of the invention, a non-transitory computer readable medium having stored thereon a computer program comprising a set of instructions, which, when executed by a computer, causes the computer to acquire calibration data from a calibration scan configured to characterize high order eddy current generated magnetic field error of an imaging system and to process the calibration data. The set of instructions also causes the computer to generate a plurality of basis coefficients and a plurality of time constants based on the processed calibration data and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a DW-EPI pulse sequence. The set of instructions also causes of the computer to apply the DW-EPI pulse sequence to acquire MR data from an imaging subject, to reconstruct an image based on the acquired MR data, and to apply high order eddy-current-generated magnetic field error correction during reconstruction of the image configured to reduce high order eddy current induced distortion in the reconstructed image.
- Various other features and advantages will be made apparent from the following detailed description and the drawings.
- The drawings illustrate embodiments presently contemplated for carrying out embodiments of the invention.
- In the drawings:
-
FIG. 1 is a schematic block diagram of an MR imaging system for use with embodiments of the invention. -
FIG. 2 is a pulse sequence diagram showing ideal gradient and RF waveforms for a single spin echo diffusion-weighted EPI scan. -
FIG. 3 is a flowchart showing an HOEC correction technique according to an embodiment of the invention. -
FIG. 4 is a flowchart showing an HOEC calibration data processing algorithm that may be used in the HOEC correction technique ofFIG. 3 according to an embodiment of the invention. -
FIG. 5 is a flowchart showing another HOEC data processing algorithm that may be used in the HOEC correction technique ofFIG. 3 according to an embodiment of the invention. -
FIG. 6 is a flowchart showing an algorithm for calculating HOEC terms that may be used in the HOEC correction technique ofFIG. 3 according to an embodiment of the invention. -
FIG. 7 is a flowchart showing an algorithm for retrospective compensation of HOEC terms that may be used in the HOEC correction technique ofFIG. 3 according to an embodiment of the invention. - Referring to
FIG. 1 , the major components of a magnetic resonance imaging (MRI)system 10 incorporating an embodiment of the invention are shown. The operation of the system is controlled for certain functions from anoperator console 12 which in this example includes a keyboard orother input device 13, acontrol panel 14, and adisplay screen 16. Theconsole 12 communicates through alink 18 with aseparate computer system 20 that enables an operator to control the production and display of images on thedisplay screen 16. Thecomputer system 20 includes a number of modules which communicate with each other through abackplane 20 a. These modules include animage processor module 22, aCPU module 24 and amemory module 26, known in the art as a frame buffer for storing image data arrays. Thecomputer system 20 communicates with aseparate system control 32 through a highspeed serial link 34. Theinput device 13 can include a mouse, joystick, keyboard, track ball, touch activated screen, light wand, voice control, card reader, push-button, or any similar or equivalent input device, and may be used for interactive geometry prescription. - The
system control 32 includes a set of modules connected together by abackplane 32 a. These include aCPU module 36 and apulse generator module 38 which connects to theoperator console 12 through aserial link 40. It is throughlink 40 that thesystem control 32 receives commands from the operator to indicate the scan sequence that is to be performed. Thepulse generator module 38 operates the system components to carry out the desired scan sequence and produces data which indicates the timing, strength and shape of the RF pulses produced, and the timing and length of the data acquisition window. Thepulse generator module 38 connects to a set ofgradient amplifiers 42, to indicate the timing and shape of the gradient pulses that are produced during the scan. Thepulse generator module 38 can also receive patient data from aphysiological acquisition controller 44 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes attached to the patient. And finally, thepulse generator module 38 connects to a scanroom interface circuit 46 which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scanroom interface circuit 46 that apatient positioning system 48 receives commands to move the patient to the desired position for the scan. - The gradient waveforms produced by the
pulse generator module 38 are applied to thegradient amplifier system 42 having Gx, Gy, and Gz amplifiers. Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly generally designated 50 to produce the magnetic field gradients used for spatially encoding acquired signals. Thegradient coil assembly 50 forms part of aresonance assembly 52 which includes a polarizingmagnet 54 and a whole-body RF coil 56. Atransceiver module 58 in thesystem control 32 produces pulses which are amplified by anRF amplifier 60 and coupled to theRF coil 56 by a transmit/receiveswitch 62. The resulting signals emitted by the excited nuclei in the patient may be sensed by thesame RF coil 56 and coupled through the transmit/receiveswitch 62 to apreamplifier 64. The amplified MR signals are demodulated, filtered, and digitized in the receiver section of thetransceiver 58. The transmit/receive switch 62 is controlled by a signal from thepulse generator module 38 to electrically connect theRF amplifier 60 to thecoil 56 during the transmit mode and to connect thepreamplifier 64 to thecoil 56 during the receive mode. The transmit/receiveswitch 62 can also enable a separate RF coil (for example, a surface coil) to be used in either the transmit or receive mode. - The MR signals picked up by the
RF coil 56 are digitized by thetransceiver module 58 and transferred to amemory module 66 in thesystem control 32. A scan is complete when an array of raw k-space data has been acquired in thememory module 66. This raw k-space data is rearranged into separate k-space data arrays for each image to be reconstructed, and each of these is input to anarray processor 68 which operates to Fourier transform the data into an array of image data. This image data is conveyed through theserial link 34 to thecomputer system 20 where it is stored in memory. In response to commands received from theoperator console 12 or as otherwise directed by the system software, this image data may be archived in long term storage or it may be further processed by theimage processor 22 and conveyed to theoperator console 12 and presented on thedisplay 16. - Embodiments of the invention correct high order eddy current (HOEC) induced diffusion gradient direction dependent distortion in Diffusion-Weighted Echo Planar Imaging (DW-EPI). As used herein, “high order” means spatial order higher than or equal to 2, as compared to linear and constant orders, which mean
order 1 and 0, respectively. However, embodiments of the invention can also be applied to linear and constant orders as well. The distortions due to HOEC may strongly affect the acquired data, especially for body applications where large FOV is used and large slice coverage is desired. As discussed herein, a system calibration may be used to characterize eddy currents of spatial orders less than or equal to P, where P is usually 3 to 5, but can be any order in general. DW-EPI pulse sequence with correction amplitudes for gradients in readout, phase encoding, and slice axes, and receiver frequency on a per-slice basis may be played to prospectively remove the effect of the quasi-linear HOEC terms. - Referring to
FIG. 2 , a DW-EPI pulse sequence 70 is shown including 90° and 180°RF pulses RF pulses RF coil 56 to generate anecho signal 76, which can be encoded with spatial information.Echo signal 76 can also be received bycoil 56 or by another coil, such as a surface coil, for use in reconstructing an image. - To spatially encode
echo signal 76 in accordance with echo planar imaging, the sequence shown inFIG. 2 further includes read-out, phase-encoding, and slice-selection gradients Gro, Gpe, and Gsl, respectively. Readout gradient Gro comprises apre-phasing pulse 78 and read-outpulses 80. Similarly, phase-encoding gradient Gpe comprises apre-phasing pulse 82 and phase-encodingpulses 84. Slice-selection gradient Gsl comprisesslice selection pulses 86 for the 90°RF pulse RF pulse 74, as well as 86 a for slice refocusing. - Still referring to
FIG. 2 , there is shown diffusion-weighting gradient Gd, used in a diffusion-weighted EPI scan. Diffusion-weighting gradient Gd comprises two equivalenttrapezoidal pulses RF pulse 74. Note that in general, Gd can contain components on all three logical axes. In the following, the read-out, phase-encoding, and slice axes (i.e., the logical axes) are denoted as u, v, and w respectively, and the components of Gd in u, v, and w axes are denoted as Gu, Gv, and Gw, respectively. The physical axes are denoted as x, y, and z. -
FIG. 3 shows anHOEC correction technique 94 according to an embodiment of the invention.Technique 94 begins atblock 96 by performing an HOEC calibration to characterize the eddy-current-generated magnetic field error of a particular imaging system such asMRI system 10 ofFIG. 1 . The calibration can either be phantom-based or local-pickup-coil-based. In the phantom-based method, gradient echo images are often collected at multiple time points following an eddy current generating gradient. In the local-pickup-coil-based method, a number of local coils, each with a small sample, are used to obtain free induction decay signals at their respective spatial locations. To use the local pickup coil method for HOEC measurement, multiple data acquisitions are often needed where the coil fixture is repositioned at each acquisition so that sufficient data are obtained for HOEC characterization. The HOEC calibration can be done as frequently as needed, but is in general only needed once per system installation or when there are significant system hardware (e.g., gradient coil) changes. A 4-dimensional (3D in space and 1D in time) eddy current field data set is generated after HOEC calibration scan. - At
block 98, data from the HOEC calibration ofblock 96 are first preprocessed and then fitted to mathematical models to characterize the underlying HOEC according to an algorithm based on which of the abovementioned calibration methods is used inblock 96.FIGS. 4 and 5 illustrate HOECdata processing algorithms block 98 of theHOEC correction technique 94 ofFIG. 3 according to an embodiment of the invention. - Referring to
FIG. 4 , HOECdata processing algorithm 116 is used when the HOEC calibration ofblock 96 oftechnique 94 is performed using the phantom-based method as described above. Data from the HOEC calibration scan are acquired atblock 120. Atblock 122, a 3D phase unwrapping is applied to the phase angle of the data, and the phase angle is scaled atblock 124 by a factor proportional to the echo time to yield a magnetic field data set. Each time point of the magnetic field data set is then spatially fitted atblock 126 to polynomial or spherical harmonic bases of order up to P to generate basis coefficients, where P is usually 3 to 5, but can be any order in general. Note that magnitude weights or masks can be optionally used during the spatial fitting. - The resulting basis coefficients are then temporal fitted along the time axis at
block 128 using a single-exponential or multi-exponential model. The end results of HOEC data processing is a set of (αmn, τmn) pairs, where αmn are basis coefficients and τmn are time constants (for notational simplicity, a single exponential is assumed) of the nth spatial basis function Bn(x, y, z), n=1,2, . . . , N, with diffusion donor axis m, where m is the x, y, or z axis. Bn(x, y, z) are assumed to be polynomial bases for convenience of discussion. Note that this is without any loss of generalization because spherical harmonic bases are linear combinations of polynomial bases and can be easily converted into polynomials. Note also that the total number of bases N=(P+1)(P+2)(P+3)/6. The (αmn, τmn) pairs are saved or stored on the host computer of the scanner for future use atblock 130. - Referring to
FIG. 5 , HOECdata processing algorithm 118 is used when the HOEC calibration ofblock 96 oftechnique 94 is performed using the local-pickup-coil-based method as described above. Data from the HOEC calibration scan are acquired atblock 132. Atblock 134, time derivatives are taken on the phase angle of the data to obtain the magnetic field offset at a coil location. Each time point of the magnetic field offset is then spatially fitted atblock 136 to polynomial or spherical harmonic bases of order up to P to generate basis coefficients, where P is usually 3 to 5, but can be any order in general. As noted above, magnitude weights or masks can be optionally used during the spatial fitting. - Similar to that described above with respect to HOEC
data processing algorithm 116, the resulting basis coefficients are then temporal fitted along the time axis atblock 138, and a set of (αmn, τmn) pairs is generated. The (αmn, τmn) pairs are saved or stored on the host computer of the scanner for future use atblock 140. - Similar to the HOEC calibration scan of
block 96 oftechnique 94, the HOEC data processing performed in HOECdata processing algorithms algorithms - Referring back to
FIG. 3 , DW-EPI protocol dependent HOEC terms are calculated atblock 100. Note thatblock 100 can handle arbitrary imaging planes. As used herein, an arbitrary imaging plane means straight axial, coronal, or sagittal scan plane, as well as any oblique plane. As shown inFIG. 6 , analgorithm 142 for calculating the HOEC terms forblock 100 ofFIG. 3 is shown. Atblock 144, the diffusion gradient components Gu, Gv, Gw of the DW-EPI pulse sequence to be used are obtained. Atblock 145, logical gradients Gu, Gv, Gw are converted to physical components Gx, Gy, Gz by applying the 3×3 axis rotation matrix R: -
- Note that R is a unitary matrix (i.e., R−1=RT).
- At
block 146, the (αmn, τmn) pairs determined viablock 98 from the HOEC calibration scan atblock 96 oftechnique 94 are obtained, and the pulse sequence type and sequence timing dependent constant, βnm, is calculated atblock 148. - Derivation of βmn can either be analytical or using convolution. Although all gradient waveforms can be included to obtain βnm, contributions from the diffusion gradients are often dominant, which allow for simplified analysis to obtain βmn. For example, when single spin echo DW-EPI is used such as that shown in
FIG. 2 , it can be derived that: -
- where t1, t2, and t3 are sequence timing related constants shown in
FIG. 2 . Note that βnm for other pulse sequences such as the dual spin echo or stimulated echo DW-EPI can also be determined analytically. - At
block 149, the HOEC physical basis coefficients dn(t) for the nth basis function βn(x, y, z) at time t after the last diffusion gradient are calculated based on the equation: -
- where Gm is the x, y, or z component of the diffusion gradient amplitude.
- At
block 150, the transpose of another rotation matrix, the basis rotation matrix F, is applied to d1(t), d2(t), . . . , dN(t) to convert them into HOEC logical basis coefficients c1(t), c2(t), . . . , cN(t): -
- where “T” denotes matrix transpose. F is an N×N matrix that transforms polynomial bases from logical to physical coordinates. The actual form of F depends on the polynomial order and how the basis functions are numbered. Without loss of generality, the bases are in the following order: 1, x, y, z, x2,xy, xz, y2, yz, z2, x3, x2y, x2z, xy2, xyz, xz2, y3, y2z, yz2, z3, . . . , where lower order bases lead higher order bases, and for bases that have the same polynomial order, the ones that have higher x exponent lead, or, in case of the same x exponent, the ones that have higher y exponent lead. F can be determined by the relationship between Bn(x, y, z) and Bn(u, v, w), where
-
- For example, for up to 3 rd order polynomials,
-
- At
block 151, all HOEC logical basis coefficients are saved for future use. - Referring back to
FIG. 3 , the protocol dependent HOEC terms can be compensated for intechnique 94 retrospectively. Retrospective compensation includes applying a DW-EPI pulse sequence to acquire data from the imaging subject atblock 102. Then, during image reconstruction atblock 104, the HOEC terms are retrospectively compensated for.FIG. 7 illustrates aretrospective compensation algorithm 162 forblock 104 according to an embodiment of the invention. Note that due to the application of F, bases are now in logical axes. - Image distortion due to HOEC fields can be corrected using a field map based method in image reconstruction according to this embodiment.
Algorithm 162 does distortion correction in the image domain on a per-slice basis. For simplicity, the w dependency in the notation is dropped. Denote the eddy current field as f(u, v). At w=w0, the value of the HOEC related field map f(x, y) can be calculated atblock 164 by the equation: -
- where the summation is done over all the spatial basis functions and t0 is the time point at which the HOEC field is approximated. The choice of t0 can be arbitrary, but it is preferable to be chosen at the time when the k-space center data is acquired (i.e., at the echo time).
- The eddy current field is then normalized by the effective gradient in the phase encoding axis GPE to obtain the pixel shift map h(u, v) at
block 166 by the equation: -
- where GPE=1/(γTL) with γ being the gyromagnetic ratio of the nucleus of interest, T being the EPI echo spacing and L being the field of view in the phase encoding axis.
- Denote the image under normal image reconstruction (i.e., without retrospective correction) as Idistorted(u, v). The image after retrospective correction Icorrected(u, v) can be obtained at
block 168 by the equation: -
- where the first term is the intensity correction and the second term is the geometry correction. Note that both Idistorted(u, v) and Icorrected(u, v) could be magnitude or complex images. According to Eqns. 3 and 4, h(u, v) are linear combinations of Bn(u, v, w0). Because Bn(u, v, w0) are polynomial bases, it is straightforward to obtain the analytical form of ∂h(u, v)/∂v, which makes it easy to obtain the intensity correction factor in Eqn. 5. Evaluation of the geometry correction term can be accomplished by image domain interpolation. HOEC-corrected image results from performing
algorithm 162. - One can see that this technique can be applied to multishot EPI (top-down or bottom-up interleaved acquisition) as well. For S-shot EPI, the effective echo spacing becomes T/S and GPE in Eqn. 4 becomes S/(γTL), which scales the overall distortion by a factor of 1/S. Note that multishot EPI can reduce TE versus single shot EPI, leading to a different t0 in Eqn. 3.
- Referring again to
FIG. 3 , the image reconstructed atblock 104 may be displayed to a user or stored on an image storage database for future use atblock 106. - According to embodiments of the invention, compensation of HOEC-induced distortions for DW-EPI can make single spin echo DW-EPI more practical. While, single spin echo has SNR and scan time benefits over dual spin echo DWI, embodiments of the invention can also be applied to dual spin echo and other variants of DW-EPI sequences to reduce distortion. Embodiments of the invention can produce a significant leverage of whole body DWI, which may have a long scan time, low SNR, and large image distortions. Brain DWI may benefit as well, especially for high b-value, large parallel imaging factor cases where SNR can be a problem.
- A technical contribution for the disclosed method and apparatus is that it provides for a computer implemented correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging.
- One skilled in the art will appreciate that embodiments of the invention may be interfaced to and controlled by a computer readable storage medium having stored thereon a computer program. The computer readable storage medium includes a plurality of components such as one or more of electronic components, hardware components, and/or computer software components. These components may include one or more computer readable storage media that generally stores instructions such as software, firmware and/or assembly language for performing one or more portions of one or more implementations or embodiments of a sequence. These computer readable storage media are generally non-transitory and/or tangible. Examples of such a computer readable storage medium include a recordable data storage medium of a computer and/or storage device. The computer readable storage media may employ, for example, one or more of a magnetic, electrical, optical, biological, and/or atomic data storage medium. Further, such media may take the form of, for example, floppy disks, magnetic tapes, CD-ROMs, DVD-ROMs, hard disk drives, and/or electronic memory. Other forms of non-transitory and/or tangible computer readable storage media not list may be employed with embodiments of the invention.
- A number of such components can be combined or divided in an implementation of a system. Further, such components may include a set and/or series of computer instructions written in or implemented with any of a number of programming languages, as will be appreciated by those skilled in the art. In addition, other forms of computer readable media such as a carrier wave may be employed to embody a computer data signal representing a sequence of instructions that when executed by one or more computers causes the one or more computers to perform one or more portions of one or more implementations or embodiments of a sequence.
- Therefore, according to an embodiment of the invention, an MRI apparatus comprises a magnetic resonance imaging (MRI) system having a plurality of gradient coils positioned about a bore of a magnet, and an RF transceiver system and an RF switch controlled by a pulse module to transmit RF signals to an RF coil assembly to acquire MR images. The MRI apparatus also comprises a computer programmed to acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy-current-generated magnetic field error of an imaging system. The computer is also programmed to process the calibration data to generate a plurality of basis coefficients and a plurality of time constants and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence. The computer is further programmed to execute a diffusion-weighted imaging scan that comprises application of a DW-EPI pulse sequence to acquire MR data from an imaging subject and reconstruction of an image based on the acquired MR data. The computer is also programmed to apply high order eddy-current-generated magnetic field error correction during image reconstruction configured to reduce high order eddy-current-induced distortion in the reconstructed image.
- According to another embodiment of the invention, a method for correcting high order eddy-current-induced distortions in diffusion-weighted echo planar imaging (DW-EPI) comprises acquiring calibration data from a calibration scan, the calibration data configured to characterize high order eddy currents of an imaging system, processing the calibration data to generate a plurality of basis coefficients and a plurality of time constants, and calculating a plurality of basis correction coefficients based on the plurality of basis coefficients and based on the plurality of time constants. The method also comprises applying a DW-EPI pulse sequence to acquire MR data from an imaging subject, reconstructing an image based on the acquired MR data, and wherein the reconstructing the image comprises applying high order eddy-current-generated magnetic field error correction configured to reduce high order eddy-current-induced distortion in the image.
- According to yet another embodiment of the invention, a non-transitory computer readable medium having stored thereon a computer program comprising a set of instructions, which, when executed by a computer, causes the computer to acquire calibration data from a calibration scan configured to characterize high order eddy current generated magnetic field error of an imaging system and to process the calibration data. The set of instructions also causes the computer to generate a plurality of basis coefficients and a plurality of time constants based on the processed calibration data and to calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a DW-EPI pulse sequence. The set of instructions also causes of the computer to apply the DW-EPI pulse sequence to acquire MR data from an imaging subject, to reconstruct an image based on the acquired MR data, and to apply high order eddy-current-generated magnetic field error correction during reconstruction of the image configured to reduce high order eddy current induced distortion in the reconstructed image.
- This written description uses examples to disclose embodiments of the invention, including the best mode, and also to enable any person skilled in the art to practice the embodiments of the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of embodiments of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.
Claims (20)
1. An MRI apparatus comprising:
a magnetic resonance imaging (MRI) system having a plurality of gradient coils positioned about a bore of a magnet, and an RF transceiver system and an RF switch controlled by a pulse module to transmit RF signals to an RF coil assembly to acquire MR images; and
a computer programmed to:
acquire calibration data from a calibration scan, the calibration data configured to characterize high order eddy current generated magnetic field error of an imaging system;
process the calibration data to generate a plurality of basis coefficients and a plurality of time constants;
calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a given pulse sequence;
execute a diffusion-weighted imaging scan comprising:
applying a DW-EPI pulse sequence to acquire MR data from an imaging subject; and
reconstructing an image based on the acquired MR data; and
apply high order eddy-current-generated magnetic field error correction during image reconstruction configured to reduce high order eddy-current-induced distortion in the reconstructed image.
2. The MRI apparatus of claim 1 wherein the computer is further programmed to calculate a high order eddy-current-related field map.
3. The MRI apparatus of claim 2 wherein the computer in being programmed to calculate the high order eddy-current-related field map, is programmed to calculate the high order eddy-current-related field map based on the equation:
where cn(t0) are basis correction coefficients of the plurality of basis correction coefficients for a basis function at a time point t0 at which the high order eddy current field is approximated, and Bn(u, v, w0) are polynomial bases.
4. The MRI apparatus of claim 2 wherein the computer is further programmed to apply one of an intensity correction and a geometry correction to data based on the high order eddy-current-related field map.
5. The MRI apparatus of claim 4 wherein the computer in being programmed to apply the one of the intensity correction and the geometry correction, is programmed to apply the one of the intensity correction and the geometry correction based on the equation:
where h(u, v) is a pixel shift map.
6. The MRI apparatus of claim 4 wherein the computer is further programmed to calculate a pixel shift map, h(u, v), based on the high order eddy-current-related field map.
7. The MRI apparatus of claim 1 wherein the computer is programmed to apply the high order eddy-current-generated magnetic field error correction for an arbitrary imaging plane.
8. The MRI apparatus of claim 1 wherein the computer, in being programmed to process the calibration data, is programmed to:
apply a 3D phase unwrapping to a phase angle of the calibration data;
scale the unwrapped calibration data to generate a magnetic field data set;
spatially fit the magnetic field data set to a harmonic basis to generate basis coefficients; and
temporally fit the basis coefficients along a time axis using one of a single-exponential model and a multi-exponential model.
9. The MRI apparatus of claim 1 wherein the computer, in being programmed to process the calibration data, is programmed to:
take a time derivative on a phase angle of the calibration data to obtain a magnetic field offset at a coil location;
spatially fit each time point of the magnetic field offset to a harmonic basis to generate basis coefficients; and
temporally fit the basis coefficients along a time axis using one of a single-exponential model and a multi-exponential model.
10. The MRI apparatus of claim 1 wherein the computer, in being programmed to calculate the plurality of basis correction coefficients, is programmed to calculate the plurality of basis correction coefficients based on the equation:
where Gm is the X, Y, or Z component of the diffusion gradient amplitude, βmn is a pulse sequence type and sequence timing related constant, αmn are basis coefficients, and τmn are time constants.
11. The MRI apparatus of claim 1 wherein the computer is further programmed to display the reconstructed image to a user.
12. A method for correcting high order eddy-current-induced distortions in diffusion-weighted echo planar imaging (DW-EPI) comprising:
acquiring calibration data from a calibration scan, the calibration data configured to characterize high order eddy currents of an imaging system;
processing the calibration data to generate a plurality of basis coefficients and a plurality of time constants;
calculating a plurality of basis correction coefficients based on the plurality of basis coefficients and based on the plurality of time constants; and
applying a DW-EPI pulse sequence to acquire MR data from an imaging subject;
reconstructing an image based on the acquired MR data; and
wherein the reconstructing the image comprises applying high order eddy-current-generated magnetic field error correction configured to reduce high order eddy-current-induced distortion in the image.
13. The method of claim 12 further comprising calculating a high order eddy-current-related field map.
14. The method of claim 13 further comprising applying one of an intensity correction and a geometry correction to data based on the high order eddy-current-related field map.
15. The method of claim 14 further comprising generating the data based on the high order eddy-current-related field map by calculating a pixel shift map, h(u, v), based on the high order eddy-current-related field map.
16. A non-transitory computer readable medium having stored thereon a computer program comprising a set of instructions, which, when executed by a computer, causes the computer to:
acquire calibration data from a calibration scan configured to characterize high order eddy current generated magnetic field error of an imaging system;
process the calibration data;
generate a plurality of basis coefficients and a plurality of time constants based on the processed calibration data;
calculate a plurality of basis correction coefficients based on the plurality of basis coefficients, the plurality of time constants, and gradient waveforms in a DW-EPI pulse sequence;
apply the DW-EPI pulse sequence to acquire MR data from an imaging subject;
reconstruct an image based on the acquired MR data; and
apply high order eddy-current-generated magnetic field error correction during reconstruction of the image configured to reduce high order eddy current induced distortion in the reconstructed image.
17. The computer readable medium of claim 16 wherein the set of instructions further causes the computer to:
calculate a high order eddy-current-related field map;
calculate a pixel shift map, h(u, v), based on the high order eddy-current-related field map; and
apply an intensity and geometry correction to the pixel shift map.
18. The computer readable medium of claim 17 wherein the set of instructions that causes the computer to calculate the high order eddy-current-related field map causes the computer to calculate the high order eddy-current-related field map based on the equation:
where cn(t0) are basis correction coefficients of the plurality of basis correction coefficients for a basis function at a time point t0 at which the high order eddy current field is approximated, and Bn(u, v, w0) are polynomial bases.
19. The method of claim 18 wherein the set of instructions that causes the computer to calculate a pixel shift map causes the computer to calculate the pixel shift map based on the equation:
where GPE=1/(γTL) with γ being the gyromagnetic ratio of the nucleus of interest, T being the EPI echo spacing, and L being the field of view in the phase encoding axis.
20. The method of claim 19 wherein the set of instructions that causes the computer to apply the intensity and geometry correction causes the computer to apply the intensity and geometry correction based on the equation:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/397,813 US20120271583A1 (en) | 2011-04-19 | 2012-02-16 | System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161476936P | 2011-04-19 | 2011-04-19 | |
US13/397,813 US20120271583A1 (en) | 2011-04-19 | 2012-02-16 | System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
Publications (1)
Publication Number | Publication Date |
---|---|
US20120271583A1 true US20120271583A1 (en) | 2012-10-25 |
Family
ID=47021995
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/364,003 Active 2033-05-02 US8942945B2 (en) | 2011-04-19 | 2012-02-01 | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
US13/397,813 Abandoned US20120271583A1 (en) | 2011-04-19 | 2012-02-16 | System and method for retrospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/364,003 Active 2033-05-02 US8942945B2 (en) | 2011-04-19 | 2012-02-01 | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
Country Status (3)
Country | Link |
---|---|
US (2) | US8942945B2 (en) |
JP (1) | JP6018401B2 (en) |
CN (1) | CN102778662B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120271584A1 (en) * | 2011-04-19 | 2012-10-25 | Dan Xu | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
CN104007407A (en) * | 2014-06-13 | 2014-08-27 | 河海大学 | Gradient coil design method with combination of fonctionelle and simulated annealing algorithm |
CN104007406A (en) * | 2014-06-13 | 2014-08-27 | 河海大学 | Method for designing elliptic cylindrical surface transverse gradient coil based on spatial alternation |
US9632162B2 (en) | 2013-12-06 | 2017-04-25 | Toshiba Medical Systems Corporation | Method of, and apparatus for, correcting distortion in medical images |
US10761172B2 (en) | 2015-08-25 | 2020-09-01 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for an eddy-current field compensation in magnetic resonance imaging |
CN113366330A (en) * | 2018-07-31 | 2021-09-07 | 海珀菲纳股份有限公司 | Low field diffusion weighted magnetic resonance imaging |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8587310B2 (en) * | 2008-09-25 | 2013-11-19 | Hitachi Medical Corporation | Magnetic resonance imaging device |
EP2402780A1 (en) * | 2010-06-23 | 2012-01-04 | Koninklijke Philips Electronics N.V. | Method of reconstructing a magnetic resonance image of an object considering higher-order dynamic fields |
CN104799859B (en) * | 2014-01-27 | 2019-01-08 | 西门子(深圳)磁共振有限公司 | A kind of bearing calibration of Diffusion-weighted imaging image and device |
JP6591208B2 (en) * | 2015-06-09 | 2019-10-16 | キヤノンメディカルシステムズ株式会社 | Magnetic resonance imaging system |
CN105891755B (en) * | 2016-02-25 | 2018-07-03 | 吉林大学 | The bearing calibration of aircraft hanging fluxgate magnetic gradient tensor instrument |
EP3589197A4 (en) * | 2017-03-01 | 2020-12-09 | The Regents of The University of California | System and method for improved diffusion-weighted imaging |
US10809339B2 (en) * | 2017-03-24 | 2020-10-20 | GE Precision Healthcare LLC | System and method for performing magnetic resonance diffusion weighted imaging of an object |
JP7253433B2 (en) | 2019-04-19 | 2023-04-06 | 富士フイルムヘルスケア株式会社 | MAGNETIC RESONANCE IMAGING APPARATUS AND ERROR MAGNETIC COMPENSATION METHOD |
CN112545484B (en) * | 2020-12-01 | 2023-08-08 | 上海联影医疗科技股份有限公司 | Medical imaging system calibration time determination method, device and computer equipment |
TWI771141B (en) * | 2021-08-06 | 2022-07-11 | 長庚大學 | Brain Extraction and Coil Shim Correction Model Establishment in Brain Imaging |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5351006A (en) * | 1992-02-07 | 1994-09-27 | Board Of Trustees Of The Leland Stanford Junior University | Method and apparatus for correcting spatial distortion in magnetic resonance images due to magnetic field inhomogeneity including inhomogeneity due to susceptibility variations |
US6335620B1 (en) * | 1998-12-22 | 2002-01-01 | Siemens Aktiengesellschaft | Method of acquiring eddy currents that are caused by switched magnetic field gradients of a magnetic resonance apparatus and that contain cross-terms |
US20020039449A1 (en) * | 2000-08-28 | 2002-04-04 | Shigetoshi Nouda | Image compression apparatus |
US20110225325A1 (en) * | 2010-03-10 | 2011-09-15 | Renesas Electronics Corporation | Data transfer circuit and method |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4698591A (en) | 1986-01-03 | 1987-10-06 | General Electric Company | Method for magnetic field gradient eddy current compensation |
JPH0654820A (en) * | 1991-09-30 | 1994-03-01 | Toshiba Corp | Magnetic resonance diagnostic device |
US5770943A (en) * | 1996-12-30 | 1998-06-23 | General Electric Company | Method for measuring and compensating for spatially and temporally varying magnetic fields induced by eddy currents |
US5864233A (en) | 1997-02-22 | 1999-01-26 | General Electric Company | Method to reduce eddy current effects in diffusion-weighted echo planar imaging |
US6211675B1 (en) | 1998-11-12 | 2001-04-03 | General Electric Company | Automatic measurement of gradient field distortion |
JP2002224083A (en) * | 2001-02-06 | 2002-08-13 | Hitachi Medical Corp | Magnetic resonance imaging diagnostic apparatus |
JP3734086B2 (en) * | 2003-03-12 | 2006-01-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Nuclear magnetic resonance imaging system |
US6815952B1 (en) * | 2003-05-12 | 2004-11-09 | The University Of Queensland | Magnetic resonance diffusion imaging with eddy-current compensation |
US7218110B2 (en) * | 2005-09-02 | 2007-05-15 | Toshiba America Mri, Inc. | Method and apparatus for diffusion magnetic resonance imaging with the effects of eddy currents compensated |
WO2008146551A1 (en) * | 2007-05-09 | 2008-12-04 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and image classification method |
CN100580476C (en) * | 2007-11-27 | 2010-01-13 | 华东师范大学 | Correcting method for quick-speed spin echo pulse sequence and application thereof |
US8154293B2 (en) * | 2009-07-31 | 2012-04-10 | General Electric Company | Apparatus and method for acquisition and reconstruction of non-CPMG propeller sequences |
US8942945B2 (en) * | 2011-04-19 | 2015-01-27 | General Electric Company | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
-
2012
- 2012-02-01 US US13/364,003 patent/US8942945B2/en active Active
- 2012-02-16 US US13/397,813 patent/US20120271583A1/en not_active Abandoned
- 2012-04-17 JP JP2012093447A patent/JP6018401B2/en active Active
- 2012-04-19 CN CN201210129813.1A patent/CN102778662B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5351006A (en) * | 1992-02-07 | 1994-09-27 | Board Of Trustees Of The Leland Stanford Junior University | Method and apparatus for correcting spatial distortion in magnetic resonance images due to magnetic field inhomogeneity including inhomogeneity due to susceptibility variations |
US6335620B1 (en) * | 1998-12-22 | 2002-01-01 | Siemens Aktiengesellschaft | Method of acquiring eddy currents that are caused by switched magnetic field gradients of a magnetic resonance apparatus and that contain cross-terms |
US20020039449A1 (en) * | 2000-08-28 | 2002-04-04 | Shigetoshi Nouda | Image compression apparatus |
US20110225325A1 (en) * | 2010-03-10 | 2011-09-15 | Renesas Electronics Corporation | Data transfer circuit and method |
Non-Patent Citations (1)
Title |
---|
Shen et al., "Correction of High-Order Eddy Current Induced Geometric Distortion in Diffusion-Weighted Echo-Planar Images", Magnetic Resonance in Medicine 52:1184-1189 (2004). * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120271584A1 (en) * | 2011-04-19 | 2012-10-25 | Dan Xu | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
US8942945B2 (en) * | 2011-04-19 | 2015-01-27 | General Electric Company | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging |
US9632162B2 (en) | 2013-12-06 | 2017-04-25 | Toshiba Medical Systems Corporation | Method of, and apparatus for, correcting distortion in medical images |
CN104007407A (en) * | 2014-06-13 | 2014-08-27 | 河海大学 | Gradient coil design method with combination of fonctionelle and simulated annealing algorithm |
CN104007406A (en) * | 2014-06-13 | 2014-08-27 | 河海大学 | Method for designing elliptic cylindrical surface transverse gradient coil based on spatial alternation |
US10761172B2 (en) | 2015-08-25 | 2020-09-01 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for an eddy-current field compensation in magnetic resonance imaging |
CN113366330A (en) * | 2018-07-31 | 2021-09-07 | 海珀菲纳股份有限公司 | Low field diffusion weighted magnetic resonance imaging |
Also Published As
Publication number | Publication date |
---|---|
JP6018401B2 (en) | 2016-11-02 |
CN102778662B (en) | 2016-08-03 |
US8942945B2 (en) | 2015-01-27 |
CN102778662A (en) | 2012-11-14 |
US20120271584A1 (en) | 2012-10-25 |
JP2012245350A (en) | 2012-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8942945B2 (en) | System and method for prospective correction of high order eddy-current-induced distortion in diffusion-weighted echo planar imaging | |
JP4114989B2 (en) | Method for compensating the magnetic field of a magnetic resonance system | |
US9778336B2 (en) | System and method for rapid, multi-shot segmented magnetic resonance imaging | |
US7495439B2 (en) | MRI method for reducing artifacts using RF pulse at offset frequency | |
US6603989B1 (en) | T2 contrast in magnetic resonance imaging with gradient echoes | |
US6369568B1 (en) | Fast spin echo phase correction for MRI system | |
US9316707B2 (en) | System and method of receive sensitivity correction in MR imaging | |
US7492153B2 (en) | System and method of parallel imaging with calibration to a separate coil | |
US8085041B2 (en) | Three-point method and system for fast and robust field mapping for EPI geometric distortion correction | |
US7847545B2 (en) | System and method for correcting flow velocity measurements in phase contrast imaging using magnetic field monitoring | |
US20100308827A1 (en) | System and method for multi-spectral mr imaging near metal | |
US20120169338A1 (en) | System And Method For Generating A Magnetic Resonance Image Using Compressed Sensing And Parallel Imaging | |
US20110103670A1 (en) | System and method for reducing blurring artifacts in multi-spectral mr imaging near metal | |
US8483457B2 (en) | System and method of image artifact reduction using self-navigated real-time phase correction in echo planar imaging | |
US8704518B2 (en) | System and method of high signal-to-noise ratio magnetic resonance imaging screening | |
US6469505B1 (en) | Method and apparatus to reduce perturbation field effects in MR images by restricting the region of interest | |
EP1372110B1 (en) | Method and system for image reconstruction | |
US6380738B1 (en) | Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system | |
US6528998B1 (en) | Method and apparatus to reduce the effects of maxwell terms and other perturbation magnetic fields in MR images | |
US7239137B2 (en) | Method and apparatus for fast spin echo (FSE) prescan phase correction | |
US10054657B2 (en) | Apparatus and method for multishot diffusion weighted imaging with array spatial pseudo-sensitivity encoding technique | |
US8680860B2 (en) | System and method for reducing localized signal fluctuation | |
US20110241672A1 (en) | Apparatus and method for parallel transmission of rf pulses in a spin echo sequence | |
US20220057467A1 (en) | Epi mr imaging with distortion correction | |
US8854039B2 (en) | System and method for quantifying MR data acquired during a UTE sequence |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GENERAL ELECTRIC COMPANY, NEW YORK Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:XU, DAN;MAIER, JOSEPH K.;KING, KEVIN F.;AND OTHERS;REEL/FRAME:027715/0089 Effective date: 20120208 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |