US10317501B2 - System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields - Google Patents
System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields Download PDFInfo
- Publication number
- US10317501B2 US10317501B2 US15/660,221 US201715660221A US10317501B2 US 10317501 B2 US10317501 B2 US 10317501B2 US 201715660221 A US201715660221 A US 201715660221A US 10317501 B2 US10317501 B2 US 10317501B2
- Authority
- US
- United States
- Prior art keywords
- recited
- parameters
- acquisition parameters
- magnetic resonance
- homogeneity
- 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.)
- Active, expires
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/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
- G01R33/3614—RF power amplifiers
-
- 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
-
- 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]
-
- 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/5659—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF 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
Definitions
- Magnetic resonance fingerprinting is an imaging technique that enables quantitative mapping of tissue or other material properties based on random or pseudorandom measurements of the subject or object being imaged. Examples of parameters that can be mapped include longitudinal relaxation time, T 1 ; transverse relaxation time, T 2 ; main magnetic field map, B 0 ; and proton density, ⁇ . MRF is generally described in U.S. Pat. No. 8,723,518, which is herein incorporated by reference in its entirety.
- the random or pseudorandom measurements obtained in MRF techniques are achieved by varying the acquisition parameters from one repetition time (“TR”) period to the next, which creates a time series of images with varying contrast.
- acquisition parameters include flip angle, radio frequency (“RF”) pulse phase, TR, echo time (“TE”), and sampling patterns, such as by modifying one or more readout encoding gradients.
- RF radio frequency
- TR repetition time
- TE echo time
- sampling patterns such as by modifying one or more readout encoding gradients.
- the data acquired with MRF techniques are compared with a dictionary of signal models, or templates, that have been generated for different acquisition parameters from magnetic resonance signal models, such as Bloch equation-based physics simulations. This comparison allows estimation of the desired physical parameters, such as those mentioned above.
- the parameters for the tissue or other material in a given voxel are estimated to be the values that provide the best signal template matching.
- the present disclosure addresses the aforementioned drawbacks by providing a method for using a magnetic resonance imaging (“MRI”) system to estimate quantitative parameters of a subject in the presence of inhomogeneous magnetic fields.
- MRI magnetic resonance imaging
- Optimized acquisition parameters are estimated.
- the acquisition parameters are optimized to direct the MRI system to generate a plurality of different signal evolutions that maximize discrimination between magnetic resonance signal patterns in a minimized number of repetition time (“TR”) periods as a function of at least one relaxation time parameter and magnetic field homogeneity.
- TR repetition time
- Data are acquired with the MRI system by directing the MRI system to perform a plurality of pulse sequences using the optimized acquisition parameters, the acquired data representing the plurality of different signal evolutions.
- Quantitative parameters of the subject are computed by comparing the acquired data with a dictionary database comprising a plurality of different signal templates, wherein each signal template includes values associated with flip angle homogeneity.
- RF radio frequency
- FIG. 1 is a flowchart setting forth the steps of an example method for using magnetic resonance fingerprinting (“MRF”) to estimate quantitative parameters from data acquired using acquisition parameters that have been optimized to reduce the number of acquisitions necessary to maximize the discrimination between signal patterns as a function of relaxation time and static magnetic field homogeneity values.
- MRF magnetic resonance fingerprinting
- FIG. 2 is an example echo-planar imaging (“EPI”) pulse sequence that can be used to acquire data using an optimized set of acquisition parameters and from which quantitative parameters can be estimated using MRF.
- EPI echo-planar imaging
- FIG. 3 is a block diagram of an example of a magnetic resonance imaging (“MRI”) system.
- MRI magnetic resonance imaging
- MRF magnetic resonance fingerprinting
- MRF techniques utilize a data acquisition scheme that causes signals from different materials or tissues to be spatially and temporally incoherent by continuously varying acquisition parameters throughout the data acquisition process.
- acquisition parameters include flip angle (“FA”), radio frequency (“RF”) pulse phase, repetition time (“TR”), echo time (“TE”), and sampling patterns, such as by modifying readout encoding gradients.
- FA flip angle
- RF radio frequency
- TR repetition time
- TE echo time
- sampling patterns such as by modifying readout encoding gradients.
- the acquisition parameters are generally varied in a pseudorandom manner.
- each material or tissue is associated with a unique signal evolution or “fingerprint,” that is a function of multiple different physical parameters, including longitudinal relaxation time, T 1 ; transverse relaxation time, T 2 ; main magnetic field map, B 0 ; and proton density, ⁇ .
- Quantitative parameter maps are then generated from the acquired signals based on a comparison of the signals to a predefined dictionary of predicted signal evolutions.
- a comparison of the signals to a predefined dictionary of predicted signal evolutions is associated with different combinations of materials and acquisition parameters.
- the comparison of the acquired signals to a dictionary can be performed using any suitable matching or pattern recognition technique. This comparison results in the selection of a signal vector, which may constitute a weighted combination of signal vectors, from the dictionary that best corresponds to the observed signal evolution.
- the selected signal vector includes values for multiple different quantitative parameters, which can be extracted from the selected signal vector and used to generate the relevant quantitative parameter maps.
- the systems and methods described here generate an optimized, ultra-short sequence of FA and TR values that amplifies signal pattern differences as a function of relaxation times (e.g., T 1 , T 2 ) and static magnetic field (e.g., B 0 ) homogeneity.
- the resulting optimized sequence of FA and TR values is robust to static magnetic field variations, and therefore offers high signal discrimination among nuclear spin species, even in the presence of inhomogeneous magnetic fields. Because the optimized sequence of FA and TR values significantly reduces the number of data points implemented in the MRF pulse sequence, Cartesian data acquisitions can be utilized.
- the optimization used to determine the sequence of FA and TR values can also be designed to maintain signal-to-noise ratio (“SNR”) in the acquired signal patterns above a minimum threshold value, which may include imposing one or more penalty terms.
- SNR signal-to-noise ratio
- the optimization can be limited to search for trajectories with flip angles above a certain minimum value (e.g., 10 or 15 degrees), or for trajectories with echo time (“TE”) or TR values within a certain range of values.
- a penalty term can be added to the optimization to maximize the average transverse magnetization of a given trajectory.
- a penalty term that enforces maximization of transverse magnetization results in changes in signal above the noise floor, which allows for discrimination between different signal patterns.
- a penalty term that can be used to maintain a minimum SNR, deviations of the on-diagonal elements from a value of 1 can be penalized. For the latter example, in a noiseless scenario the on-diagonal elements should be exactly 1, so any deviation from would be the result of noise.
- FA homogeneity is accounted for in the dictionary, thereby allowing for recovery of coil homogeneity profiles in addition to quantitative parameters. Adding FA homogeneity to the dictionary also enhances the accuracy of the matching process by adding another quantitative parameter in addition to T 1 , T 2 , ⁇ (proton density), and off-resonance.
- MRI magnetic resonance imaging
- Typical MRI scanners require highly homogeneous static magnetic fields and RF coils, which increases the total cost of the MRI system hardware.
- FIG. 1 a flowchart is illustrated as setting forth the steps of one non-limiting example method for producing a quantitative parameter map from data acquired using acquisition parameters that have been selected to reduce the number of acquisitions necessary to desirably sample the relaxation parameter space.
- the method begins by generating a vector, or schedule, of acquisition parameters that has been selected to reduce the number of acquisitions necessary to sufficiently sample the relaxation parameter space while also amplifying signal pattern differences as a function of relaxation time and static magnetic field homogeneity, as indicated at step 102 .
- the sequence of acquisition parameters (e.g., FA and TR) is selected to be symmetric, such as temporally symmetric.
- a symmetric sequence of acquisition parameters nuclear spin magnetization will be at least partially rewound back to its original state at the end of each pulse sequence.
- the FA and TR series are played back half way through the imaging sequence, which as stated above partially rewinds the spins (e.g., to get back to the initial position of longitudinal magnetization).
- acquisition sequences can follow in rapid succession without having to wait as long between sequences for magnetization recovery.
- the selection or optimization of acquisition parameters may include providing an initial, randomly-generated seed vector of the acquisition parameters to be selected or optimized.
- This seed vector may have a length, N, and may be used to simulate the signal for a range, P, of relaxation parameters and static magnetic field homogeneity values.
- a first range, P 1 , of T 1 values can be used together with a second range, P 2 , of T 2 values and a third range, P 3 , of B 0 values.
- combinations of ranges of these parameters can be used (e.g., a range of only one relaxation time parameter and a range of static magnetic field homogeneity values).
- the seed vector and simulated relaxation parameters are used to form an N ⁇ P matrix, A, where
- the diagonal elements of this dot product matrix, D indicate the closeness of a match between a schedule resulting from a given relaxation parameter or static magnetic field homogeneity value and itself.
- the diagonal elements are, therefore, equal to one.
- the off-diagonal elements of the dot product matrix, D indicate the quality of matching between every two different elements of the matrix, A. Discriminating between T 1 , T 2 , and B 0 homogeneity values in the matching process requires that the dot product of a measured trajectory with the pre-computed trajectory that is stored in the dictionary be high for the correct T 1 , T 2 , and B 0 homogeneity values and, ideally, zero for all others.
- a model can be utilized.
- One non-limiting example of a model is the following optimization problem:
- f(x) is the function to simulate the schedule and compute the dot product matrix, D, given a vector, x, of acquisition parameters.
- a penalty term, ⁇ is applied as well to avoid minimizing the on-diagonal elements.
- another penalty term can be added to Eqn. (2) to maintain an SNR value above a minimum threshold value.
- a constrained non-linear solver can be used to solve Eqns. (2) or (3).
- Using the optimized schedule reduces the total number of acquisitions needed to sufficiently sample the physical parameter space, but also samples this physical parameter space in a manner that increases the discrimination between different parameters, thereby reducing the scan time required to perform MRF. Additional time savings can be achieved by using rapid data acquisition schemes, such as EPI and segmented EPI.
- the method for estimating quantitative parameters continues by acquiring data by directing an MRI system to perform pulse sequences using the optimized acquisition parameters, as indicated at step 104 .
- the pulse sequence can be any suitable pulse sequence for obtaining MRF data, but in some embodiments the pulse sequence used can be an EPI pulse sequence, which may be a spin-echo (“SE”) EPI sequence or a gradient-recalled echo (“GRE”) EPI sequence.
- SE spin-echo
- GRE gradient-recalled echo
- EPI pulse sequence has the added benefit that k-space can be fully sampled in a rapid fashion, which reduces the total number of required acquisitions, and therefore contrast doses. Reducing the total number of acquisitions not only provides a time saving measure, but also reduces the computational burden of matching the acquired data to a pre-computed dictionary because a smaller dictionary can be used with the reduced number of acquisitions.
- the method continues by reconstructing images from the acquired data, as indicated at step 106 .
- reconstruction can be achieved using a conventional Fourier transform-based reconstruction, which introduces another reduction in the computational burden compared to conventional MRF applications that require complex reconstruction techniques to reconstruct images from data acquired using non-Cartesian trajectories.
- Quantitative parameter maps are then estimated by, for example, matching the reconstructed images to one or more pre-computed dictionaries, as indicated at step 108 .
- Conventional matching algorithms can be used; however, in some configurations, an adaptive matching algorithm, such as the one described in co-pending PCT Application No. PCT/US15/11948, which is herein incorporated by reference in its entirety, can also be used.
- quantitative parameter maps can be generated by associating the estimated quantitative parameters with spatial locations in the imaged field-of-view, as indicated at step 110 .
- the pre-computed dictionaries include additional entries associated with FA homogeneity values. Using such dictionaries not only increases the accuracy of the matching process in the presence of inhomogeneous magnetic fields and when using inhomogeneous RF transmitters, but can provide for the simultaneous estimation of coil profiles in addition to quantitative parameter maps. Thus, in some examples RF coil transmission homogeneity profiles are simultaneously estimated with the quantitative parameter maps, as indicated at step 112 .
- FIG. 2 an example pulse sequence diagram of a modified spin-echo EPI technique for use in accordance with some implementations of the systems and methods described in the present disclosure, is shown.
- a selective RF excitation pulse 202 is played out in the presence of a slice select gradient pulse 204 .
- the flip angle, a values of the RF excitations are usually 90 degrees to maximize the signal to noise ratio (“SNR”).
- SNR signal to noise ratio
- the present disclosure employs a variable FA, a, that may be varied according to the optimized acquisition parameters, unlike previous MRF approaches that utilize pseudo-randomized parameter values.
- the slice select gradient pulse 204 is followed by a rephasing pulse 206 , which operates to rephase excited spins in the selected slice. Subsequently, to facilitate the formation of spin echoes, the orientation of the transverse magnetization of excited spins may be reversed by application of a selective 180 degree RF excitation pulse 208 produced in the present of a slice select gradient pulse 210 . Readout and phase-encoding waveforms are played following RF excitation 202 to produce a series of spatially encoded echoes 220 .
- a plurality of gradient waveforms are played out for traversing k-space in accordance with a desired trajectory.
- this trajectory is a Cartesian k-space trajectory.
- a series of readout gradient pulses 216 with alternating polarity are played along the readout direction to produce a series of spin-echo signals 220 .
- the readout gradient pulses 216 are spaced apart, and phase-encoding gradient blips 218 are played out therebetween to separately phase encode each acquired echo 220 .
- multiple views may be acquired for each excitation in the TR, to encompass either k-space segments or fully sample k-space in accordance with the Nyquist criterion.
- fully sampled k-space on a Cartesian grid avoids a regridding and simplifies the reconstruction process.
- multiple slices may be acquired during each scan. For instance, delays associated with wait times necessary for magnetization recovery of successive phase encoding k-space lines may limit total scan times. Thus, in certain applications, these otherwise dead times may be utilized for acquisition of multiple slices to facilitate three-dimensional (“3D”) data acquisition.
- 3D three-dimensional
- a segmented EPI sequence may also be used, depending on the desired resolution and tissue parameter values.
- ⁇ T a wait time, ⁇ T, may be needed following each segment to ensure the initial magnetization is the same for all segments.
- additional slices can be acquired, for example, if 3D acquisition is desired. As such, after a last phase-encoding k-space line of a segment is acquired for a particular slice, additional k-space segments from other slices may be interleaved therein.
- multiple acquisitions can be performed, wherein scan parameter values for each acquisition can vary from one excitation to the next. For example, ten acquisitions can take place in a manner described above, although other values are possible.
- TR, FA, and other parameter values may be varied in accordance with optimized or desired trajectories, such as those described above; in a random manner; or in a pseudorandom manner.
- the parameters can be varied sequentially, simultaneously, or in accordance with a differential sensitivity requirement for specific parameters of interest.
- the MRI system 300 includes an operator workstation 302 that may include a display 304 , one or more input devices 306 (e.g., a keyboard, a mouse), and a processor 308 .
- the processor 308 may include a commercially available programmable machine running a commercially available operating system.
- the operator workstation 302 provides an operator interface that facilitates entering scan parameters into the MRI system 300 .
- the operator workstation 302 may be coupled to different servers, including, for example, a pulse sequence server 310 , a data acquisition server 312 , a data processing server 314 , and a data store server 316 .
- the operator workstation 302 and the servers 310 , 312 , 314 , and 316 may be connected via a communication system 340 , which may include wired or wireless network connections.
- the pulse sequence server 310 functions in response to instructions provided by the operator workstation 302 to operate a gradient system 318 and a radiofrequency (“RF”) system 320 .
- Gradient waveforms for performing a prescribed scan are produced and applied to the gradient system 318 , which then excites gradient coils in an assembly 322 to produce the magnetic field gradients G x , G y , and G z that are used for spatially encoding magnetic resonance signals.
- the gradient coil assembly 322 forms part of a magnet assembly 324 that includes a polarizing magnet 326 and a whole-body RF coil 328 .
- RF waveforms are applied by the RF system 320 to the RF coil 328 , or a separate local coil to perform the prescribed magnetic resonance pulse sequence.
- Responsive magnetic resonance signals detected by the RF coil 328 , or a separate local coil are received by the RF system 320 .
- the responsive magnetic resonance signals may be amplified, demodulated, filtered, and digitized under direction of commands produced by the pulse sequence server 310 .
- the RF system 320 includes an RF transmitter for producing a wide variety of RF pulses used in MRI pulse sequences.
- the RF transmitter is responsive to the prescribed scan and direction from the pulse sequence server 310 to produce RF pulses of the desired frequency, phase, and pulse amplitude waveform.
- the generated RF pulses may be applied to the whole-body RF coil 328 or to one or more local coils or coil arrays.
- the RF system 320 also includes one or more RF receiver channels.
- An RF receiver channel includes an RF preamplifier that amplifies the magnetic resonance signal received by the coil 328 to which it is connected, and a detector that detects and digitizes the I and Q quadrature components of the received magnetic resonance signal.
- phase of the received magnetic resonance signal may also be determined according to the following relationship:
- the pulse sequence server 310 may receive patient data from a physiological acquisition controller 330 .
- the physiological acquisition controller 330 may receive signals from a number of different sensors connected to the patient, including electrocardiograph (“ECG”) signals from electrodes, or respiratory signals from a respiratory bellows or other respiratory monitoring devices. These signals may be used by the pulse sequence server 310 to synchronize, or “gate,” the performance of the scan with the subject's heart beat or respiration.
- ECG electrocardiograph
- the pulse sequence server 310 may also connect to a scan room interface circuit 332 that receives signals from various sensors associated with the condition of the patient and the magnet system. Through the scan room interface circuit 332 , a patient positioning system 334 can receive commands to move the patient to desired positions during the scan.
- the digitized magnetic resonance signal samples produced by the RF system 320 are received by the data acquisition server 312 .
- the data acquisition server 312 operates in response to instructions downloaded from the operator workstation 302 to receive the real-time magnetic resonance data and provide buffer storage, so that data is not lost by data overrun. In some scans, the data acquisition server 312 passes the acquired magnetic resonance data to the data processor server 314 . In scans that require information derived from acquired magnetic resonance data to control the further performance of the scan, the data acquisition server 312 may be programmed to produce such information and convey it to the pulse sequence server 310 . For example, during pre-scans, magnetic resonance data may be acquired and used to calibrate the pulse sequence performed by the pulse sequence server 310 .
- navigator signals may be acquired and used to adjust the operating parameters of the RF system 320 or the gradient system 318 , or to control the view order in which k-space is sampled.
- the data acquisition server 312 may also process magnetic resonance signals used to detect the arrival of a contrast agent in a magnetic resonance angiography (“MRA”) scan.
- MRA magnetic resonance angiography
- the data acquisition server 312 may acquire magnetic resonance data and processes it in real-time to produce information that is used to control the scan.
- the data processing server 314 receives magnetic resonance data from the data acquisition server 312 and processes the magnetic resonance data in accordance with instructions provided by the operator workstation 302 .
- processing may include, for example, reconstructing two-dimensional or three-dimensional images by performing a Fourier transformation of raw k-space data, performing other image reconstruction algorithms (e.g., iterative or backprojection reconstruction algorithms), applying filters to raw k-space data or to reconstructed images, generating functional magnetic resonance images, or calculating motion or flow images.
- Images reconstructed by the data processing server 314 are conveyed back to the operator workstation 302 for storage.
- Real-time images may be stored in a data base memory cache, from which they may be output to operator display 302 or a display 336 .
- Batch mode images or selected real time images may be stored in a host database on disc storage 338 .
- the data processing server 314 may notify the data store server 316 on the operator workstation 302 .
- the operator workstation 302 may be used by an operator to archive the images, produce films, or send the images via a network to other facilities.
- the MRI system 300 may also include one or more networked workstations 342 .
- a networked workstation 342 may include a display 344 , one or more input devices 346 (e.g., a keyboard, a mouse), and a processor 348 .
- the networked workstation 342 may be located within the same facility as the operator workstation 302 , or in a different facility, such as a different healthcare institution or clinic.
- the networked workstation 342 may gain remote access to the data processing server 314 or data store server 316 via the communication system 340 . Accordingly, multiple networked workstations 342 may have access to the data processing server 314 and the data store server 316 . In this manner, magnetic resonance data, reconstructed images, or other data may be exchanged between the data processing server 314 or the data store server 316 and the networked workstations 342 , such that the data or images may be remotely processed by a networked workstation 342 .
Abstract
Description
for the M different ranges. This matrix, A, can then be used to calculate a dot product matrix,
D=A T A (1).
M=√{square root over (I 2 +Q 2)} (4);
Claims (13)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/660,221 US10317501B2 (en) | 2016-07-26 | 2017-07-26 | System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201662366812P | 2016-07-26 | 2016-07-26 | |
US15/660,221 US10317501B2 (en) | 2016-07-26 | 2017-07-26 | System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields |
Publications (2)
Publication Number | Publication Date |
---|---|
US20180031667A1 US20180031667A1 (en) | 2018-02-01 |
US10317501B2 true US10317501B2 (en) | 2019-06-11 |
Family
ID=61009467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/660,221 Active 2037-11-29 US10317501B2 (en) | 2016-07-26 | 2017-07-26 | System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields |
Country Status (1)
Country | Link |
---|---|
US (1) | US10317501B2 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10180476B2 (en) * | 2015-11-25 | 2019-01-15 | The General Hospital Corporation | Systems and methods for segmented magnetic resonance fingerprinting dictionary matching |
US10429468B2 (en) * | 2016-03-25 | 2019-10-01 | The General Hospital Corporation | Simultaneous dynamic contrast enhanced and dynamic susceptibility magnetic resonance imaging using magnetic resonance fingerprinting |
US10488481B2 (en) * | 2017-02-22 | 2019-11-26 | The General Hospital Corporation | Systems and methods for multislice magetic resonance fingerprinting |
NL2022890B1 (en) * | 2019-04-08 | 2020-10-15 | Umc Utrecht Holding Bv | Parameter map determination for time domain magnetic resonance |
CN114217255B (en) * | 2021-11-29 | 2022-09-20 | 浙江大学 | Rapid liver multi-parameter quantitative imaging method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7337205B2 (en) * | 2001-03-21 | 2008-02-26 | Apple Inc. | Matrix multiplication in a vector processing system |
US20120235678A1 (en) * | 2011-03-18 | 2012-09-20 | Case Western Reserve University | Nuclear magnetic resonance (nmr) fingerprinting |
US20160349342A1 (en) * | 2014-02-11 | 2016-12-01 | The General Hospital Corporation | Systems and methods for acceleration magnetic resonance fingerprinting |
-
2017
- 2017-07-26 US US15/660,221 patent/US10317501B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7337205B2 (en) * | 2001-03-21 | 2008-02-26 | Apple Inc. | Matrix multiplication in a vector processing system |
US20120235678A1 (en) * | 2011-03-18 | 2012-09-20 | Case Western Reserve University | Nuclear magnetic resonance (nmr) fingerprinting |
US20160349342A1 (en) * | 2014-02-11 | 2016-12-01 | The General Hospital Corporation | Systems and methods for acceleration magnetic resonance fingerprinting |
Also Published As
Publication number | Publication date |
---|---|
US20180031667A1 (en) | 2018-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10422845B2 (en) | Systems and methods for steady-state magnetic resonance fingerprinting | |
US10598750B2 (en) | 3D balanced EPI magnetic resonance fingerprinting | |
US9897675B2 (en) | Magnetic resonance fingerprinting (MRF) with simultaneous multivolume acquisition | |
US10488481B2 (en) | Systems and methods for multislice magetic resonance fingerprinting | |
US10261155B2 (en) | Systems and methods for acceleration magnetic resonance fingerprinting | |
US10317501B2 (en) | System and method for magnetic resonance fingerprinting in the presence of inhomogeneous magnetic fields | |
US9726742B2 (en) | System and method for iteratively calibrated reconstruction kernel for accelerated magnetic resonance imaging | |
US11373392B2 (en) | System and method for magnetic resonance fingerprinting using a plurality of pulse sequence types | |
US10564243B2 (en) | System and method for phase cycling magnetic resonance fingerprinting (phc-MRF) | |
US10429468B2 (en) | Simultaneous dynamic contrast enhanced and dynamic susceptibility magnetic resonance imaging using magnetic resonance fingerprinting | |
US10598747B2 (en) | System and method for simultaneous multislice magnetic resonance fingerprinting with variable radio frequency encoding | |
US10761167B2 (en) | System and method for generating a magnetic resonance fingerprinting dictionary using semi-supervised learning | |
US11875509B2 (en) | System and method for determining undersampling errors for a magnetic resonance fingerprinting pulse sequence | |
US10527695B2 (en) | Systems and methods for efficient magnetic resonance fingerprinting scheduling | |
US11313931B2 (en) | System and method for quantifying T1, T2 and resonance frequency using rosette trajectory acquisition and read segmented reconstruction | |
US10429479B2 (en) | Rapid measurement of perfusion using optimized magnetic resonance fingerprinting | |
US11131733B2 (en) | System and method for magnetic resonance fingerprinting with non-locally sequential sampling of k-space | |
US11867785B2 (en) | Dual gradient echo and spin echo magnetic resonance fingerprinting for simultaneous estimation of T1, T2, and T2* with integrated B1 correction | |
US10466321B2 (en) | Systems and methods for efficient trajectory optimization in magnetic resonance fingerprinting | |
US11079448B2 (en) | System and method for proton density mapping and receiver bias correction using magnetic resonance fingerprinting (MRF) | |
US11519987B2 (en) | Magnetic resonance fingerprinting thermometry | |
US11747421B2 (en) | System and method for quantifying perfusion using a dictionary matching | |
US10670680B2 (en) | System and method for motion insensitive magnetic resonance fingerprinting | |
US20190120921A1 (en) | Blind Source Separation in Magnetic Resonance Fingerprinting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: THE GENERAL HOSPITAL CORPORATION, MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ROSEN, MATTHEW;SARRACANIE, MATHIEU;COHEN, OURI;SIGNING DATES FROM 20161130 TO 20161205;REEL/FRAME:045728/0949 Owner name: THE UNIVERSITY COURT OF THE UNIVERSITY OF ABERDEEN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROCHE, LIONEL MARC;REEL/FRAME:045729/0054 Effective date: 20170308 Owner name: THE GENERAL HOSPITAL CORPORATION, MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:THE UNIVERSITY COURT OF THE UNIVERSITY OF ABERDEEN;REEL/FRAME:045358/0666 Effective date: 20170518 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
CC | Certificate of correction | ||
AS | Assignment |
Owner name: THE GOVERNMENT OF THE UNITED STATES AS REPRESENTED BY THE SECRETARY OF THE ARMY, MARYLAND Free format text: CONFIRMATORY LICENSE;ASSIGNOR:MASSACHUSETTS GENERAL HOSPITAL;REEL/FRAME:058551/0720 Effective date: 20200625 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |