IP Library Granted Patent US 12,189,014
Granted Patent B2
US 12,189,014 · App. 18/321,708 · Granted Jan 7, 2025

Eddy current mitigation systems and methods

Inventors: Hadrien A. Dyvorne (New York, NY); Cedric Hugon (Guilford, CT); Rafael O'Halloran (Guilford, CT); Laura Sacolick (Guilford, CT)
Assignee: Hyperfine Operations, Inc.
G01R33/56518G01R33/58G01R33/445
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,189,014
App. No.
18/321,708
Granted
Jan 7, 2025
Kind
B2
Abstract

Techniques for compensating for presence of eddy currents during the operation of a magnetic resonance imaging (MRI) system in accordance with a pulse sequence, the pulse sequence comprising a gradient waveform associated with a target gradient field. The techniques include: compensating for presence of eddy currents during operation of the MRI system at least in part by correcting the gradient waveform using a nonlinear function of a characteristic of the gradient waveform to obtain a corrected gradient waveform; and operating the MRI system in accordance with the corrected gradient waveform to generate the target gradient field.

Claims (42)

1. A method of operating a magnetic resonance imaging (MRI) system, the MRI system including at least one radio frequency (RF) coil configured to apply an RF signal in accordance with a pulse sequence, the pulse sequence comprising a gradient waveform associated with a target gradient field and an RF waveform associated with a target RF signal, the method comprising:

compensating for a presence of eddy currents during operation of the MRI system at least in part by correcting the RF waveform to obtain a corrected RF waveform, the correcting comprising:

determining a shifted B 0 magnetic field strength waveform based on a corrected gradient waveform obtained at least in part by applying a first nonlinear function of a characteristic of the gradient waveform to the gradient waveform; and

correcting the RF waveform using the shifted B 0 magnetic field strength waveform to obtain the corrected RF waveform; and

operating the MRI system in accordance with the corrected RF waveform to generate the target RF signal.

2. The method of claim 1 , wherein determining the shifted B 0 magnetic field strength waveform comprises applying a first transformation to the corrected gradient waveform to obtain the shifted B 0 magnetic field strength waveform.

3. The method of claim 2 , wherein applying the first transformation to the corrected gradient waveform comprises filtering the corrected gradient waveform with a first filter determined using the first nonlinear function.

4. The method of claim 1 , further comprising determining the first nonlinear function based on first eddy current calibration data obtained from the MRI system, the first eddy current calibration data describing a B 0 magnetic field strength as a function of the characteristic of the gradient waveform.

5. The method of claim 1 , wherein the method further comprises:

compensating for the presence of eddy currents during operation of the MRI system at least in part by correcting the gradient waveform using a second nonlinear function of a characteristic of the gradient waveform to obtain a corrected gradient waveform; and

operating the MRI system in accordance with the corrected gradient waveform to generate the target gradient field.

6. The method of claim 5 , wherein the characteristic of the gradient waveform comprises an amplitude of the gradient waveform, a direction of the target gradient field, shape of the gradient waveform, and/or slew rate of the gradient waveform.

7. The method of claim 5 , wherein the second nonlinear function is a polynomial function, an exponential function, a piecewise polynomial function, or a piecewise constant function of the characteristic of the gradient waveform.

8. The method of claim 5 , wherein the second nonlinear function is determined based on second eddy current calibration data obtained from the MRI system.

9. The method of claim 8 , further comprising:

obtaining the second eddy current calibration data using a field probe to measure eddy currents in the MRI system during its operation; and

determining the second nonlinear function using the second eddy current calibration data.

10. The method of claim 5 , wherein correcting the gradient waveform comprises:

applying the second nonlinear function to the gradient waveform to obtain a scaled gradient waveform; and

applying a second transformation to the scaled gradient waveform to obtain a transformed gradient waveform, the second transformation being determined using second eddy current calibration data obtained from the MRI system.

11. The method of claim 10 , wherein applying the second transformation to the scaled gradient waveform comprises filtering the scaled gradient waveform with a filter determined using the second eddy current calibration data.

12. The method of claim 11 , wherein filtering the scaled gradient waveform is performed in a time domain at least in part by performing a convolution.

13. The method of claim 11 , wherein filtering the scaled gradient waveform is performed in a domain other than a time domain.

14. The method of claim 11 , wherein the filter is a pre-emphasis filter, and the method further comprises determining the pre-emphasis filter using an inverse Laplace transformation.

15. The method of claim 10 , further comprising:

subtracting the scaled gradient waveform from the transformed gradient waveform to obtain a correction waveform; and

combining the correction waveform with the gradient waveform to obtain the corrected gradient waveform.

16. The method of claim 15 , wherein the first nonlinear function is a nonlinear function for compensating for first eddy currents and the second nonlinear function is a nonlinear function for compensating for second eddy currents, wherein the compensating comprises:

compensating for the presence of eddy currents using multiple nonlinear functions of the characteristic of the gradient waveform to obtain a corrected gradient waveform, the multiple nonlinear functions including the first nonlinear function and the second nonlinear function.

17. The method of claim 5 , wherein operating the MRI system in accordance with the corrected gradient waveform comprises generating one or more preparation gradient field pulses prior to generating the target gradient field.

18. The method of claim 1 , wherein the pulse sequence comprises one of a diffusion weighted imaging (DWI) pulse sequence, a diffusion weighted steady-state free precession (DW-SSFP) pulse sequence, or a fast spin echo (FSE) pulse sequence.

19. At least one non-transitory computer-readable storage medium storing processor-executable instructions that, when executed by a magnetic resonance imaging (MRI) system, cause the MRI system to perform a method of operating the MRI system by applying a plurality of RF signals in accordance with a pulse sequence, the pulse sequence comprising a gradient waveform associated with a target gradient field and an RF waveform associated with a target RF signal, the method comprising:

compensating for a presence of eddy currents during operation of the MRI system at least in part by correcting the RF waveform to obtain a corrected RF waveform, the correcting comprising:

determining a shifted B 0 magnetic field strength waveform based on a corrected gradient waveform obtained at least in part by applying a first nonlinear function of a characteristic of the gradient waveform to the gradient waveform; and

correcting the RF waveform using the shifted B 0 magnetic field strength waveform to obtain the corrected RF waveform; and

operating the MRI system in accordance with the corrected RF waveform to generate the target RF signal.

20. A magnetic resonance imaging (MRI) system, comprising:

at least one controller configured to operate one or more of a plurality of magnetics components in accordance with a pulse sequence, the pulse sequence comprising a gradient waveform associated with a target gradient field and an RF waveform associated with a target RF signal, wherein the at least one controller is configured to:

compensate for a presence of eddy currents during operation of the MRI system at least in part by correcting the RF waveform to obtain a corrected RF waveform, the correcting comprising:

determining a shifted B 0 magnetic field strength waveform based on a corrected gradient waveform obtained at least in part by applying a first nonlinear function of a characteristic of the gradient waveform to the gradient waveform; and

correcting the RF waveform using the shifted B 0 magnetic field strength waveform to obtain the corrected RF waveform; and

operate the MRI system in accordance with the corrected RF waveform to generate the target RF signal.

Assignments (3)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jun 6, 2023
From: DYVORNE, HADRIEN A.; HUGON, CEDRIC; O'HALLORAN, RAFAEL; SACOLICK, LAURA
To: HYPERFINE RESEARCH, INC.
Reel/Frame 063871/0184 →
CHANGE OF NAME Recorded Jun 6, 2023
From: HYPERFINE RESEARCH, INC.
To: HYPERFINE, INC.
Reel/Frame 063882/0083 →
CHANGE OF NAME Recorded Jun 6, 2023
From: HYPERFINE, INC.
To: HYPERFINE OPERATIONS, INC.
Reel/Frame 063882/0092 →
Continuity (3)
Continuation 16994005 · Aug 14, 2020
Provisional Application 62887212 · Aug 15, 2019
Related Publication 20230324489A1 · Oct 12, 2023
References Cited (87)
US 4980641A · Breneman et al. · 1990 [cited by applicant]
US 5126672A · Le Roux · 1992 [cited by applicant]
US 5227728A · Kaufman et al. · 1993 [cited by applicant]
US 5250901A · Kaufman et al. · 1993 [cited by applicant]
US 5332969A · Tsuruno et al. · 1994 [cited by applicant]
US 5652517A · Maki et al. · 1997 [cited by applicant]
US 5864233A · Zhou et al. · 1999 [cited by applicant]
US 6191582B1 · Zur · 2001 [cited by applicant]
US 6201987B1 · Dumoulin · 2001 [cited by examiner]
US 6291997B1 · King et al. · 2001 [cited by applicant]
US 6335620B1 · Weissenberger · 2002 [cited by applicant]
US 6392411B1 · Goto · 2002 [cited by applicant]
US 6437566B1 · Heid · 2002 [cited by applicant]
US 6483305B1 · Miyamoto · 2002 [cited by applicant]
US 9541616B2 · Rothberg et al. · 2017 [cited by applicant]
US 9547057B2 · Rearick et al. · 2017 [cited by applicant]
US 9625544B2 · Poole et al. · 2017 [cited by applicant]
US 9645210B2 · McNulty et al. · 2017 [cited by applicant]
US 9817093B2 · Rothberg et al. · 2017 [cited by applicant]
US 10145913B2 · Hugon et al. · 2018 [cited by applicant]
US 10145922B2 · Rothberg et al. · 2018 [cited by applicant]
US 10222434B2 · Poole et al. · 2019 [cited by applicant]
US 10274561B2 · Poole et al. · 2019 [cited by applicant]
US 10281540B2 · Mileski et al. · 2019 [cited by applicant]
US 10281541B2 · Poole et al. · 2019 [cited by applicant]
US 10310037B2 · McNulty et al. · 2019 [cited by applicant]
US 10416264B2 · Sofka et al. · 2019 [cited by applicant]
US 10551452B2 · Rearick et al. · 2020 [cited by applicant]
US 10591561B2 · Sacolick et al. · 2020 [cited by applicant]
US 10709387B2 · Poole et al. · 2020 [cited by applicant]
US 11156688B2 · O'Halloran et al. · 2021 [cited by applicant]
US 20010041819A1 · Goto · 2001 [cited by applicant]
US 20020097049A1 · Goto · 2002 [cited by applicant]
US 20040046554A1 · Carlini · 2004 [cited by examiner]
US 20060192558A1 · Miyawaki et al. · 2006 [cited by applicant]
US 20100148774A1 · Kamata · 2010 [cited by applicant]
US 20120098535A1 · Kaneta et al. · 2012 [cited by applicant]
US 20130154642A1 · Sueoka · 2013 [cited by applicant]
US 20130187650A1 · Pfeuffer et al. · 2013 [cited by applicant]
US 20130234708A1 · Goora · 2013 [cited by examiner]
US 20140125333A1 · Hanada · 2014 [cited by examiner]
US 20140218031A1 · Lee et al. · 2014 [cited by applicant]
US 20150008921A1 · Lee · 2015 [cited by examiner]
US 20160128592A1 · Rosen et al. · 2016 [cited by applicant]
US 20160131727A1 · Sacolick et al. · 2016 [cited by applicant]
US 20160245891A1 · Ookawa · 2016 [cited by applicant]
US 20180238978A1 · McNulty et al. · 2018 [cited by applicant]
US 20190038233A1 · Poole et al. · 2019 [cited by applicant]
US 20190324098A1 · McNulty et al. · 2019 [cited by applicant]
US 20190353723A1 · Dyvorne et al. · 2019 [cited by applicant]
US 20190353726A1 · Poole et al. · 2019 [cited by applicant]
US 20200022611A1 · Nelson et al. · 2020 [cited by applicant]
US 20200022612A1 · McNulty et al. · 2020 [cited by applicant]
US 20200034998A1 · Schlemper et al. · 2020 [cited by applicant]
US 20200041588A1 · O'Halloran et al. · 2020 [cited by applicant]
US 20200045112A1 · Sacolick et al. · 2020 [cited by applicant]
US 20200058106A1 · Lazarus et al. · 2020 [cited by applicant]
US 20200200844A1 · Boskamp et al. · 2020 [cited by applicant]
US 20200209334A1 · O'Halloran et al. · 2020 [cited by applicant]
US 20200209335A1 · O'Halloran et al. · 2020 [cited by applicant]
US 20200289019A1 · Schlemper et al. · 2020 [cited by applicant]
US 20200289022A1 · Coumans et al. · 2020 [cited by applicant]
US 20200294229A1 · Schlemper et al. · 2020 [cited by applicant]
US 20200294282A1 · Schlemper et al. · 2020 [cited by applicant]
US 20200294287A1 · Schlemper et al. · 2020 [cited by applicant]
US 20200337587A1 · Sacolick et al. · 2020 [cited by applicant]
US 20200355765A1 · Chen et al. · 2020 [cited by applicant]
US 20210048498A1 · Dyvorne et al. · 2021 [cited by applicant]
US 20220043094A1 · O'Halloran et al. · 2022 [cited by applicant]
CN 1258001A · 2000 [cited by applicant]
CN 1336558A · 2002 [cited by applicant]
CN 1374069A · 2002 [cited by applicant]
EP 1004892A1 · 2000 [cited by applicant]
EP 1197760A2 · 2002 [cited by applicant]
JP H0910186A · 1997 [cited by applicant]
Eichler, et al., “Implementation of the First Order Reversal Curve Method for Identification of Weight Function in Preisach Model for Ferromagnetics,” ELEKTRO 11th International Conference; IEEE; pp. 602-607; May 16-18,… [cited by applicant]
First Office Action and Search Report on CN App. No. 201980093021.3 dated Oct. 27, 2023 (with English translation, 29 pages). [cited by applicant]
Invitation to Pay Additional Fees for International Application No. PCT/US2019/061663 mailed Mar. 5, 2020. [cited by applicant]
International Search Report and Written Opinion for International Application No. PCT/US2019/061663 mailed Jun. 24, 2020. [cited by applicant]
International Search Report and Written Opinion for International Application No. PCT/US2020/046445 mailed Dec. 3, 2020. [cited by applicant]
International Preliminary Report on Patentability for International Application No. PCT/US2019/061663 mailed Jul. 8, 2021. [cited by applicant]
International Preliminary Report on Patentability for International Application No. PCT/US2020/046445 mailed Dec. 3, 2020. [cited by applicant]
De Zanche et al., NMR probes for measuring magnetic fields and field dynamics in MR systems. Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine. Jul. 2008… [cited by applicant]
Li et al., Finite element analysis of gradient z-coil induced eddy currents in a permanent MRI magnet. Journal of Magnetic Resonance. Jan. 1, 2011;208(1):148-55. [cited by applicant]
Novak et al., Difficulty in identification of Preisach hysteresis model weighting function using first order reversal curves method in soft magnetic materials. Applied Mathematics and Computation. Feb. 15, 2018;319:469-… [cited by applicant]
Van Vaals et al., Optimization of eddy-current compensation. Journal of Magnetic Resonance (1969). Oct. 15, 1990;90(1):52-70. [cited by applicant]
Wansapura et al., Temperature mapping of frozen tissue using eddy current compensated half excitation RF pulses. Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in… [cited by applicant]