Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images -> Monitor Keywords

Site News |

Monitor Keywords |

Monitor Archive |

Organizer |

Account Info |

10/05/06 - USPTO Class 324 | 61 views | #20060220645 | Prev - Next | About this Page

Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images

Title: Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images

Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images description/claims

The Patent Description & Claims data below is from USPTO Patent Application 20060220645, Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images.

Brief Patent Description - Full Patent Description - Patent Application Claims

BACKGROUND OF THE INVENTION

[0001] Conventional magnetic resonance imaging (MRI) methods rely on a linear relationship between magnetic fields and spatial positions of the nuclear spins within the object being imaged. This relationship is realized by using a combination of homogeneous magnets and linear gradient coils. Deviations from this linear relationship can cause undesired image distortions and sometimes artifacts.

[0002] Field anomalies are deviations of magnetic fields from prescribed values. Anomalies can be static as well as sequence-induced. Static anomalies may occur from, for example, field inhomogeneities produced by the magnet used to generate the main static magnetic field. Sequence-induced anomalies may occur from, for example, eddy-currents or the Maxwell fields.

[0003] Two remedies often used for correcting anomalies in magnetic field are: (i) minimization of inhomogeneities by shimming the magnetic field using special shimming hardware, and (ii) post data acquisition image correction to compensate for the effects of the field anomalies. These remedies require accurate information about the spatial distribution of the magnetic field at the time when the image data is acquired.

[0004] Conventionally, maps of the magnetic field are generated using dedicated field-mapping sequences. Since images are acquired using a different sequence, such field maps contain no information about sequence-induced field anomalies relevant to the imaging sequence. In cases when the sequence-induced changes are significant, such as during balanced steady-state free precession (SSFP) sequences, it is preferable that the field maps be generated using the exact same sequence that is used for imaging.

[0005] In addition, field maps are often generated from the phases of the acquired images. As anyone skilled in the art understands, phases of acquired magnetic resonance (MR) signals can be affected by factors other than the static magnetic field. Such factors include, but not limited to, chemical shifts and radio-frequency (RF) field orientation. For example, in the single quadrature Dixon (SVD) method recently introduced for water-fat separated SSFP, maps of magnetic field are obtained from the phases of the isolated echo images. Since the phases of the echo signals are influenced by field inhomogeneities as well as chemical shifts, the magnetic field maps so generated are susceptible to interference from chemical shifts.

[0006] There is therefore a need for improved methods to generate magnetic field maps that are relevant to imaging sequences, wherein the maps are not degraded by interferences such as that from chemical shifts.

BRIEF DESCRIPTION OF THE INVENTION

[0007] A system has been developed for mapping the magnetic field during a balanced steady-state free precession (SSFP) sequence. A field map is generated using the data acquired during the same phase-incremented or frequency-shifted SSFP sequences used to generate the images. The system maps the magnetic field by analyzing the magnitude images pixel-by-pixel. The system need not rely on phase information.

[0008] A field map is generated by analyzing the magnitude images acquired using a phase-incremented or frequency-shifted balanced SSFP sequence. Pixel intensities are surveyed as a function of the phase increment or frequency shift used during the data acquisition. The precession angles of the nuclear spins during each time of repetition (TR) period are determined for each pixel from the radio frequency (RF) phase increment or frequency shift that yields a minimum pixel intensity. The map of the precession angles is unwrapped to form a field map of the magnetic field.

[0009] The magnetic field map so generated includes information regarding the static magnetic field inhomogeneities and any effects of sequence-induced field changes. The magnetic field map can be used to automatically shim the magnetic field so as to compensate for the static field inhomogeneities as well as to compensate for the sequence-induced field anomalies during SSFP. The magnetic field map can also be used during post-data acquisition processing to correct for the effects of the field anomalies, such as in the processes for separation of water and fat signals.

BRIEF DESCRIPTION OF THE DRAWINGS

[0010] FIG. 1 is a schematic diagram of an MRI system.

[0011] FIG. 2 is an exemplary balanced SSFP imaging sequence.

[0012] FIG. 3 is a graph of the theoretical spectrum of the magnitude of the steady-state transverse magnetization as a function of the precession angle acquired in a balanced SSFP sequence.

[0013] FIG. 4 is a flow chart of steps for acquiring phase-incremented or frequency-shifted SSFP data and generating maps of the magnetic field.

[0014] FIG. 5 is a schematic diagram of non-uniformly sampled k-space data that is sorted into data for imaging and data for field mapping.

[0015] FIG. 6 is a flow chart of steps for acquiring and processing phase-incremented or frequency-shifted SSFP data for water-fat separation.

DETAILED DESCRIPTION OF THE INVENTION

[0016] As shown in FIG. 1, a magnetic resonance imaging (MRI) system typically includes a large magnet 10 to impose a static magnetic field (B.sub.0), gradient coils 14 for imposing spatially distributed gradient magnetic fields (G.sub.x, G.sub.y, and G.sub.z) along three orthogonal coordinates, and RF coils 15 and 16 to transmit and receive RF signals to and from the selected nuclei of the object being imaged. The object 13 lies on a movable table 12 such that a portion of the object to be imaged is moved, in three-dimensions, in an "imaging volume" 11 between the magnet and coils, which defines a field of view (FOV) of the MRI system.

[0017] To acquire MRI data, the MRI system generates magnetic gradient and RF nutation pulses via MRI pulse sequence controllers 17 and 18 that operate under the control of a programmable processor 19, e.g., a workstation computer 24. In addition, the processor 19 controls a gradient pulse amplifier 20, and RF source and amplifier circuits 21 and 22. The MR signal circuits (RF detector) 25 are suitably interfaced with MR signal RF coils 15, 16 located within the shielded MRI system gantry. The received MR RF echo signal responses are digitized by a digitizer 23 and passed to the processor 19, which may include an array processors or the like for image processing and suitable computer program storage media (not shown) wherein programs are stored and selectively utilized so as to control the acquisition and processing of MR signal data and to produce images on a display of control terminal 24. The MRI system control terminal 24 may include a suitable keyboard, touch screen or other input devices for exerting operator control. Images may also be recorded directly on film, stored electronically or printed on a suitable media by a printing device.

[0018] FIG. 2 schematically shows a three-dimensional (3D) magnetic resonance imaging sequence with fully refocused steady-state free precession (SSFP). A radio frequency (RF) excitation pulse 30 is applied to the object being scanned. The RF pulse 30 is repeatedly applied with a time of repetition (TR). An RF echo signal (S) is received during a data acquisition (ADC) period 32. The received echo signal occurs after a time to echo (TE) period that begins at the center of the RF pulse 30. Gradient magnetic fields are applied in the x, y and z directions to select and encode a set of slices (Gslice) of the object for imaging and to spatially phase encode (Gphase) the excited nuclei of the selected slices. A readout magnetic field gradient (Gread) is applied during the ADC period to frequency-encode the received echo signal.

[0019] Steady-state free precession (SSFP) is a technique used to generate MRI signals from precessing hydrogen nuclei that do not completely return to their thermal equilibrium state. The SSFP sequence uses a series of RF excitation pulses and magnetic gradient pulses applied at repetition times (TR) significantly shorter than the spin-lattice (T1) and the spin-spin (T2) relaxation times of hydrogen nuclei within the object being imaged. In a fully refocused SSFP sequence, the magnetic gradient pulses are fully balanced, i.e., the total area (gradient moments) of all gradient pulses in each TR period is zero for each gradient channel.

[0020] Theory of Magnetic Field Mapping During SSFP:

[0021] During a TR period of the SSFP sequence, the nuclear spins precess to accumulate a precession angle (o) that depends on the field offset: .phi.=.gamma..DELTA.B.sub.0TR

[0022] The magnitude of the SSFP signal (S) depends on this spin precession angle: S = a .times. .times. e I .times. .times. .PHI. + b c .times. .times. cos .function. ( .PHI. ) + d a = - ( 1 - E 1 ) .times. E 2 .times. sin .function. ( .alpha. ) b = ( 1 - E 1 ) .times. sin .function. ( .alpha. ) c = E 2 .function. ( E 1 - 1 ) .times. ( 1 + cos .function. ( .alpha. ) ) d = 1 - E 1 .times. cos .function. ( .alpha. ) - ( E 1 - cos .function. ( .alpha. ) ) .times. E 2 2 E 1 = e - TR / T 1 E 2 = e - TR / T 2

[0023] where S is the steady-state signal generated in a SSFP sequence; .DELTA.B.sub.0 is the field offset that includes static field inhomogeneities as well as sequence-induced field anomalies; .gamma. is the gyromagnetic ratio of the nuclei imaged; .alpha. is the flip angle used by the SSFP sequence; TR is the time-of-repetition of the SSFP sequence; T.sub.1 and T.sub.2 are spin-lattice and spin-spin relaxation times of the nuclear spins, respectively.

[0024] FIG. 3 is a graph of the magnitude of the steady-state transverse magnetization acquired in a balanced SSFP sequence as a ratio to the equilibrium magnetization. FIG. 3 shows a theoretical spectrum of magnitudes of the SSFP signal assuming TR=10 ms, a 55.degree., T.sub.1=500 ms and T.sub.2=80 ms.

[0025] As is evident from FIG. 3, the SSFP signal is smallest when the spin precession angle is zero, i.e. when the spin precession frequency matches that of the RF carrier of the MRI system. In a phase-incremented SSFP sequence, the phase of the RF excitation pulse is incremented from one TR period to next in accordance with .theta..sub.i+1=.theta..sub.i+.delta. where .theta. represents the RF phase and .delta. represents the RF phase increment. This phase increment scheme has the same effect on the SSFP signals as that of a frequency shift with an equivalent frequency offset of .DELTA. .times. .times. f = .delta. 2 .times. .pi. .times. .times. TR .

[0026] The phase-incremented SSFP signal is minimal when the spin precession angle matches that of the phase increment of the RF pulse. By finding the RF phase increment at which the minimum RF signal occurs, the spin precession angle of the nuclear spins for each pixel of the images can be determined.

[0027] Because the spectrum of SSFP signal has a period of 2.pi., the measured spin precession angles are wrapped into between -.pi. and +.pi.. Field offset is related to the measured spin precession angle by: .phi..sub.meas=.gamma..DELTA.B.sub.0TR.+-.n2.pi.

[0028] Unwrapping of the .phi..sub.meas map yields a magnetic field map (.DELTA.B.sub.0) .DELTA. .times. .times. B 0 = unwrap .times. { .PHI. meas } .gamma. .times. .times. TR

[0029] FIG. 4 shows a flow chart of an exemplary method for generating a magnetic field map from phase-incremented or frequency-shifted SSFP acquisition. In step 100, multiple scans of SSFP are performed with phase cycling or frequency-shifting. Each scan is performed with a prescribed value of phase-increment between the repetitions or an equivalent value of frequency-shift. The scans are taken of the same object that is being imaged using the exact same imaging sequence.

[0030] For phase-incremented SSFP acquisition, the phase of the RF pulse is incremented between scan repetitions. The phase increment (step 101) is varied from scan to scan to cover a range of 0 to 2.pi., as follows: .phi..sub.i=.phi..sub.i-1+.delta..sub.j .delta..sub.j=j2.pi./N

[0031] where i is the index for repetition in a range of 0 to M-1; j is the index for scan in a range of 0 to N-1, and M and N are the total number of repetitions per scan and the total number of scans, respectively.

[0032] For frequency-shifted SSFP acquisition, the RF carrier frequency is incremented (step 101) from scan to scan according to: f.sub.j=f.sub.0+.DELTA.f.sub.j .DELTA.f.sub.j=j2.pi./(NTR)

[0033] where f.sub.0 is the RF carrier frequency without offset; f.sub.j is the RF carrier frequency used for scan j; .DELTA.f.sub.j is the RF frequency offset used for scan j, and TR is the repetition time of the SSFP sequence.

[0034] A single set of scans can be used both for field mapping and imaging. Alternatively, a non-uniform scheme of phase increment or frequency shift can be used in step 100. In such non-uniform scheme, the center of the k-space is acquired more frequently than the outer portions of the k-space data.

[0035] FIG. 5 illustrates an example of such non-uniformly acquired k-space data. In this example, the center portion of k-space is sampled twice as frequently as the outer portion of the k-space. For field mapping, only the center portion is used so as to cover a larger number of steps of phase increments or frequency shifts. The full k-space data, that is less frequently acquired but has higher spatial resolution, is used to obtain images. Such non-uniform schemes allow for field mapping with increased steps of phase increments or frequency shifts without significantly increasing the total scan time. Increased steps of phase increments or frequency shifts improve the accuracy in maps of the magnetic field.

[0036] In the case of non-uniform acquisition, step 102 sorts the acquired k-space data into sets for imaging and sets for field mapping. This step is skipped in the case of uniform acquisition wherein the same data is used for imaging as well as for field mapping.

[0037] In step 103, the field mapping data are transformed from k-space to image space. Algorithms for transforming k-space data to image data are conventional and well known to persons of ordinary skill in the art of MRI. The exact algorithm for the transformation can be selected based on the method used for acquiring the k-space data. For a conventional Fourier acquisition, a Fast Fourier Transform (FFT) may be used to convert the k-space data to images.

[0038] In step 104, pixel intensities are analyzed between the scans. For each pixel, the scan index that has the minimum intensity is identified as j.sub.min. In step 105, the left and right neighbors of j.sub.min are first identified as follows: j.sub.left=j.sub.min-1 and j.sub.right=j.sub.min+1 if 0<j.sub.min<N-1 j.sub.left=j.sub.N-1 and j.sub.right=j.sub.min+1 if j.sub.min=0 j.sub.left=j.sub.min-1 and j.sub.right=j.sub.0 if j.sub.min=N-1

[0039] The adjacent three data points, (I.sub.left, J.sub.left), (I.sub.min, J.sub.min), and (I.sub.right, J.sub.right) are fitted to a quadratic function for improving the precision of the minima location: I.sub.j=I.sub.0+a(j-d.sub.0).sup.2

[0040] where I.sub.j and j are the pixel intensity and scan index, respectively, of the three above-identified pairs; and I.sub.0 and d.sub.0 are the minimum intensity and its scan index, respectively, that the quadratic fitting produces. The scan index (d.sub.0) need not be an integer.

[0041] The scan index (d.sub.0) is used to determine the phase increment responsible for generating minimum pixel intensity as follows: .delta..sub.0=d.sub.02.pi./N

[0042] In step 106, steps 104 and 105 are repeated for all pixels to generate a precession map that represents the spin precession angles pixel-by-pixel during each TR that yield minimum pixel intensities.

[0043] In step 107, the precession map is unwrapped to generate the field map as follows: .DELTA.B.sub.0=unwrap{.delta..sub.0}/(.gamma.TR)

[0044] where .DELTA.B.sub.0 is the field offset map; unwrap{ } represents the phase-unwrapping process such as that described in U.S. Pat. No. 5,909,119; .gamma. is the gyromagnetic ratio of the nuclei imaged, and TR is the repetition time of the SSFP sequence.

[0045] In step 107, for cases where both water and fat contribute to the final images, the field map may be produced from a precession map that is scaled such that water and fat signals have a difference of 2.pi. in precession angle that is then removed by the unwrapping process: .DELTA. .times. .times. B 0 = unwrap .times. { k .times. .times. .delta. 0 } / ( .gamma. k TR ) k = 1 TR .DELTA. .times. .times. f water - fat

[0046] where .DELTA.f.sub.water-fat is the difference in resonance frequency in Hertz between the water and fat nuclear spins.

[0047] The field map (.DELTA.B.sub.0) represents the true magnetic fields during the actual SSFP imaging sequence. This map may be used to set the shimming hardware in the MRI system to compensate the field anomalies that are reflected in the field map. Determining the settings of a shimming system according to field maps is well known to persons of ordinary skill in the art. An example of shimming an MRI system using a magnetic field map is described in the article by Kim et al., "Regularized Higher Order in Vivo Shimming", Journal of Magnetic Resonance in Medicine (MRM) 48:715-722 (2002).

[0048] The field map may also be used to separate water and fat signals in SSFP imaging. FIG. 6 is a flow chart showing steps for acquiring and processing phase-incremented or frequency-shifted SSFP data with water-fat separation. In step 200, phase-incremented or frequency-shifted SSFP data are acquired in accordance with the above description for field mapping. The TR and TE of the sequence are selected according to: TR=2TE=(2n+1)/(2.DELTA.f.sub.water-fat), n=0,1,2 . . .

[0049] The non-uniform scheme of data acquisition can be used to increase the number of phase increments or frequency shifts for field mapping. In step 201, the acquired k-space data is sorted to data set for field mapping and data set for water-fat imaging. In step 202, the above-described method for field mapping using magnitude SSFP images is implemented to yield maps of the magnetic field during the imaging sequence.

[0050] In step 203, the scans in the data set for imaging are combined to separate the individual echo signals: S 0 = 1 N .times. j = 0 N - 1 .times. S j S - 1 = 1 N .times. j = 0 N - 1 .times. S j .times. e I .function. ( j2.pi. ) / N

[0051] where i= {square root over (-1)}, is the scan index, N is the total number of scans in the data set, and S.sub.0 and S.sub.-1 are echo signals with distinctive echo formation pathways.

[0052] In step 204, S.sub.0 and S.sub.-1 are processed to transform from k-space to image space, using algorithms dependent on the method by which the k-space is filled, generating I.sub.0 and I.sub.-1 respectively.

[0053] In step 205, the isolated image data are phase corrected using the field map generated in step 202 as follows: I.sub.0'=I.sub.0e.sup.-i.gamma..DELTA.B0.sup.TE I.sub.-1'I.sub.-1e.sup.o.gamma..DELTA.B.sup.0.sup.(TR-TE)

[0054] Water and fat images are generated, in step 206, according to: I.sub.water=Real(I'.sub.0)+Real(I'.sub.-1) I.sub.fat=Imag(I.sub.0')-Imag(I.sub.-1')

[0055] While the invention has been described in connection with what is presently considered to be the most practical and preferred embodiment, it is to be understood that the invention is not to be limited to the disclosed embodiment, but on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.

Brief Patent Description - Full Patent Description - Patent Application Claims

Click on the above for other options relating to this Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images patent application.
###

How KEYWORD MONITOR works... a FREE service from FreshPatents
1. Sign up (takes 30 seconds). 2. Fill in the keywords to be monitored.
3. Each week you receive an email with patent applications related to your keywords.
Start now! - Receive info on patent apps like Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images or other areas of interest.
###

Previous Patent Application:
Rf transmission method, mri apparatus, and auxiliary coil
Next Patent Application:
Magnetic resonance imaging device
Industry Class:
Electricity: measuring and testing

###

Design/code © 2009 FreshContext LLC/Freshpatents.com. Website Terms and Conditions - Also read: About this Page
Patent data source: patents published by the United States Patent and Trademark Office (USPTO)
Information published here is for research/educational purposes only (and in conjunction with our Keyword Monitor) and is not meant to be used in place of the full USPTO patent document/images or a comprehensive patent archive search. Complete official applications are on file at the USPTO and often contain additional data/images (Freshpatents is currently text-only). FreshPatents.com is not affiliated with or endorsed by the USPTO or firms/individuals or products/designs/ideas related to listed patents.

FreshPatents.com Support
Thank you for viewing the Magnetic field mapping during ssfp using phase-incremented or frequency-shifted magnitude images patent info.
IP-related news and info

Results in 0.09232 seconds

Other interesting Feshpatents.com categories:
Accenture , Agouron Pharmaceuticals , Amgen , AT&T , Bausch & Lomb , Callaway Golf 174

* Protect your Inventions
* US Patent Office filing

Provisional Patent
Utility Patent

PATENT INFO