Quantitative MR imaging of water and fat using a quadruple-echo sequence

An MRI apparatus for generating MRI data and quantitating fractions of a first species of nuclei having a first Larmor frequency and a second species of nuclei having a second Larmor frequency different from said first Larmor frequency, said apparatus comprising:

• a) means for magnetically orienting nuclei in a body to be imaged;
• b) means for imposing a slice-selection gradient on the nuclei to sensitize a planar slice of nuclei for imaging;
• c) means for nutating the selected nuclei by a RF transmission pulse;
• d) means for imposing a phase distribution on the nutated nuclei by a phase-encoding and a readout gradient pulses;
• e) means for reversing the readout gradient to generate and record in a first raw data matrix a first field-echo, with the nuclei of said first species and said second species having a phase difference of substantially 180.degree. at the center of said first field echo;
• f) means for reversing the readout gradient after said first field-echo to generate and record in a second raw data matrix a second field-echo, with the nuclei of said first species and said second species having a phase difference of substantially zero at the center of said second field-echo;
• g) means for again reversing the readout gradient after said second field-echo to generate and record in a third data matrix a third field-echo, with the nuclei of said first species and said second species having a phase difference of substantially 180.degree. at the center of said third field-echo;
• h) means for again reversing the readout gradient after said third field-echo to generate and record in a fourth raw data matrix a fourth field-echo, with the nuclei of said first species and said second species having a substantially zero degree phase difference at the center of said fourth field-echo;
• i) means for varying said phase-encoding gradient pulses and repeating means (b) through (h) to complete said raw data matrices for said phase-encoding directions;
• j) means for controlling RF and field gradient pulses so that the time intervals between the center of said RF and said first field-echo and between the centers of said successive field-echoes are same T;
• k) means for controlling the nutation angle of said RF pulse and of repeating said acquisition steps for said new data matrices with different TR for generating two different sets of said raw data matrices with TR.sup.2 =2 TR.sup.1 and with said RF pulse having a nutation angle of 90.degree.;
• and
• l) means for processing said raw data matrices to quantitate fractions of said first species and second species.

A method for quantitating fractions of a first species of nuclei having a first Larmor frequency and a second species of nuclei having a second Larmor frequency different from said first Larmor frequency, said method comprising the steps of:

• a) generating and acquiring four successive and equally distanced field-echoes, wherein the phase difference between said first species and said second species are alternately 180.degree., 0 and 180.degree., and 0 at the corresponding echo centers;
• b) repeating step (a) with various phase-encoding to complete four raw data matrices for all phase-encoding and frequency-encoding directions, wherein said first field-echo is in a first raw data matrix, said second field-echo in a second raw data matrix, said third field-echo in a third raw data matrix, and said fourth field-echo in a fourth raw data matrix;
• c) time-reversing said second and fourth field-echoes of said raw data matrices;
• d) Fourier transforming said raw data matrices to yield images;
• e) evaluating phase factors of field inhomogeneities with a guided region-growing phase unwrapping algorithm;
• f) phase-correcting said images in accordance to said phase factors;
• g) determining whether a region contains both said species of nuclei or only one of the two said species;
• h) correcting for T.sub.2 * relaxation attenuation based on the outcome of step (g);
• i) correcting for T.sub.1 relaxation attenuation;
• j) calculating equilibrium magnetization of said first species and of said second species;
• k) determining fractions of said first and second species;
• l) repeating steps (g) through (k) for all regions of interest;
• and
• m) generating fraction images of said first and second species.

A propagated MRI signal sequence embodied in a plurality of magnetic field generating and detection coils comprising a MRI apparatus, comprising:

• an RF nuclei nutation pulse;
• an imaging region selection magnetic gradient producing pulse;
• a phase encoding magnetic gradient producing pulse;
• and
• a series of read-out magnetic gradient producing pulses having alternating polarity, the read-out pulses being arranged to occur sequentially with respect to the center of the nutation pulse such that four field-echo NMR signals are produced by an object being imaged by the MRI apparatus, the field-echo signals having an inter-echo time, T, that is determined according to a known chemical-shift difference between water and fat such that an angular phase difference of .pi. radians evolves between water and fat NMR signal components during the inter-echo time.

In an MRI system, a method of obtaining NMR imaging data for use in generating a fat fraction image, comprising the steps of:

• generating a RF excitation pulse;
• generating a phase-encoding magnetic gradient pulse after generating the RF excitation pulse;
• generating a first read-out magnetic gradient pulse after generating the RF excitation pulse and the phase-encoding magnetic gradient pulse;
• generating a second read-out magnetic gradient pulse after generating the first read-out magnetic gradient pulse;
• generating a third read-out magnetic gradient pulse after generating
• the second read-out magnetic gradient pulse;
• generating a fourth read-out magnetic gradient pulse after generating the third read-out magnetic gradient pulse;
• wherein field-echo signals separated by an equal inter-echo time period (T) are timely elicited by controlling a timing and polarity of said first through fourth read-out magnetic gradients with respect to said RF excitation pulse such that NMR signals produced by a water species of nuclei develop an angular phase difference of .pi. radians (180.degree.) with respect to the NMR signals produced by a fat species of nuclei during the inter-echo time period;
• and
• acquiring first, second, third and fourth NMR field-echo signals emitted by excited fat and water species of nuclei during said first through fourth read-out magnetic gradient pulses.

In an MRI system, a method of generating a fat fraction magnetic resonance image of a region of an object comprised of a water species of nuclei having a first Larmor frequency and a fat species of nuclei having a second Larmor frequency different from said first Larmor frequency, said method comprising:

• a) generating a first plurality of quadruple field-echo NMR imaging sequences corresponding to a first MRI scan, said first plurality of imaging sequences having first, second, third and fourth read-out magnetic gradient pulses and a sequence repetition period (TR) of a first predetermined duration;
• wherein first, second, third and fourth NMR field-echo signals emitted by excited fat and water species of nuclei during each imaging sequence are acquired;
• b) generating a second plurality of quadruple field-echo NMR imaging sequences corresponding to a second MRI scan, said second plurality of imaging sequences having first, second, third and fourth read-out magnetic gradient pulses and a sequence repetition period (TR) of a second predetermined duration;
• wherein first, second, third and fourth NMR field-echo signals emitted by excited fat and water species of nuclei during each imaging sequence are acquired;
• and
• c) processing the NMR field-echo signals acquired in steps (a) and (b) to generate fat fraction magnetic resonance image data.

In an MRI system, a method of processing NMR field-echo signals used to produce a species fraction image constituting either a fat-fraction image or a water-fraction image, comprising the steps of:

• removing effects of magnetic field inhomogeneities from field-echo signal data by applying a region growing phase unwrapping algorithm to determine signal phase;
• determining if NMR signal components of fat species of nuclei and water species of nuclei are present in a field-echo signal corresponding to a pixel area of the fraction image;
• determining if image data for a pixel region of a fraction image is comprised of only water species of NMR signal components or only fat species of NMR signal components;
• depending on species present, determining values for water and or fat species of NMR signal components;
• compensating for T.sub.2 * relaxation effects;
• determining equilibrium magnetization values for water and fat species signals which will act to compensate for the effects of T.sub.1 relaxation;
• and
• for each pixel area of the fraction image, determining a fraction value corresponding to a percentage of signal due to either a fat species or a water species of nuclei out of the total signal due to both fat and water species of nuclei, wherein said fraction value for each pixel area is used in producing either a fat-fraction image or a water-fraction image.

A computer program for developing water and fat species fraction image data from NMR signal data acquired via a quadruple field-echo MRI sequence, said program comprising instructions for controlling a computer to:

• remove effects of magnetic field inhomogeneities by a applying a region growing phase unwrapping algorithm to determine signal phase;
• determine if NMR signal components of fat species of nuclei and water species of nuclei are present in a field-echo signal corresponding to a pixel area of the fat-fraction image;
• determine if image data for a pixel region of an image is comprised of only water species of NMR signal components or only fat species of NMR signal components;
• determine T.sub.2 * factor values for water and/or fat signal components, if present;
• compensate for T.sub.2 * relaxation effects;
• determine equilibrium magnetization values for water and fat species signals which compensate for the effects of T.sub.1 relaxation;
• determine a fat fraction value, .function., for each pixel area of the fat fraction image according to a percentage of signal due to a fat species of nuclei out of the total signal produced by both fat and water species of nuclei for each pixel area;
• and
• determine a water fraction value, w, for each pixel area of the water fraction image according to a percentage of signal due to a water species of nuclei out of the total signal produced by both fat and water species of nuclei for each pixel area.

An apparatus for developing either water or fat species-fraction image data from MRI quadruple field-echo sequence data, said apparatus comprising a computer including a storage memory and I/O devices, said memory having stored therein rules for performing NMR signal analysis, said computer programmed to:

• remove effects of magnetic field inhomogeneities by applying a region growing phase unwrapping algorithm to determine signal phase;
• determine if NMR signal components of fat species of nuclei and water species of nuclei are present in a field-echo signal corresponding to a pixel area of the fat-fraction image;
• determine if image data for a pixel region of an image is comprised of only water species of NMR signal components or only fat species of NMR signal components;
• determine values for water and/or fat species signal components, if present;
• compensate for T.sub.2 * relaxation effects;
• determine equilibrium magnetization values for water and fat species signals which will act to compensate for the effects of T.sub.1 relaxation;
• and
• determine either a fat fraction value or a water fraction value for each pixel area of the species-fraction image according to a percentage of signal due to either a fat or a water species of nuclei out of the total signal due to both fat and water species of nuclei.