US7417430B2 - Continuous moving-table MRI contrast manipulation and/or update of scanning parameters - Google Patents
Continuous moving-table MRI contrast manipulation and/or update of scanning parameters Download PDFInfo
- Publication number
- US7417430B2 US7417430B2 US11/568,321 US56832105A US7417430B2 US 7417430 B2 US7417430 B2 US 7417430B2 US 56832105 A US56832105 A US 56832105A US 7417430 B2 US7417430 B2 US 7417430B2
- Authority
- US
- United States
- Prior art keywords
- magnetic resonance
- resonance imaging
- data
- subject
- acquired
- 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.)
- Expired - Fee Related, 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/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/56375—Intentional motion of the sample during MR, e.g. moving table imaging
-
- 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/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
-
- 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/5605—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by transferring coherence or polarization from a spin species to another, e.g. creating magnetization transfer contrast [MTC], polarization transfer using nuclear Overhauser enhancement [NOE]
-
- 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/5607—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood imaging
-
- 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/56375—Intentional motion of the sample during MR, e.g. moving table imaging
- G01R33/56383—Intentional motion of the sample during MR, e.g. moving table imaging involving motion of the sample as a whole, e.g. multistation MR or MR with continuous table motion
Definitions
- the following relates to the magnetic resonance imaging arts. It finds particular application in cancer screening, angiography, and other diagnostic imaging procedures advantageously performed over the whole body or large portions thereof, and will be described with particular reference thereto. However, it also finds application in continuous moving-table magnetic resonance imaging generally.
- Magnetic resonance scanners have a limited field of view that is typically substantially smaller than the average person. Thus, magnetic resonance imaging is not readily used to perform “whole body” or other extended region-of-interest imaging in which the region of interest is larger than the scanner field of view. This limits the usefulness of magnetic resonance imaging for diagnostic tasks such as cancer screening, angiography, or the like.
- One approach for magnetic resonance imaging of extended region of interest larger than the scanner field of view is the multi-station method.
- imaging progresses in discrete steps. During each step, one scanner field of view is imaged. Subsequently, the subject is advanced through the scanner bore by a distance equal to the length of the scanner field of view in the axial direction, and another scanner field of view is imaged. The successively imaged fields of view are sewn together to form an image of the extended region of interest.
- the multi-station method has certain disadvantages.
- the start-and-stop motion repeated between each successive station can be disturbing to the patient.
- the imaging must stop between each discrete stepwise advancement of the imaging subject through the scanner bore, which increases the total imaging session time and generates substantial discontinuities at the boundaries between successive scanner fields of view.
- Another approach for imaging an extended region of interest larger than the scanner field of view is continuous moving table magnetic resonance imaging.
- the patient is continuously axially advanced through the scanner bore, typically at a constant table velocity, and magnetic resonance imaging data is acquired during the continuous table advancement.
- This approach may be less disturbing to the patient because it eliminates the discrete start-and-stop table movements of the multi-station method.
- discontinuities between successively imaged fields of view are decreased by omission of the temporally intervening discrete stepwise advancement.
- Discontinuities can be still further decreased in the continuous moving table approach by using a reduced field of view in the axial direction, that is, by imaging a thick slab field of view which is substantially thinner in the axial direction than the total scanner field of view. It will be appreciated that using a reduced field of view in the multi-station method is generally not advisable since it substantially increases the number of discrete advancement steps or imaged stations.
- the continuously moving imaging subject corresponds to an imaging subject coordinate system which is moving with reference to the magnetic fields and other aspects of the stationary magnetic resonance imaging scanner. This creates difficult timing issues in ensuring collection of a complete, but substantially non-redundant, data set.
- patient-dependent tuned imaging parameters such as coil loading, resonance frequency, patient-induced susceptibility and flip angle variations typically change during the continuous patient advancement.
- each voxel defined with respect to the patient that is, defined in the imaging subject coordinate system
- continuous moving table imaging has been applied mainly to contrast-enhanced angiography.
- the progression of an injected magnetic contrast-enhancing bolus is followed through the circulatory system.
- Continuous table motion facilitates tracking the dynamic bolus advancement.
- the contrast agent reduces the T 1 value of blood to increase the blood/tissue contrast. Because the contrast is provided principally by the contrast agent, imaging artifacts due to the continuous patient advancement are largely suppressed.
- a continuous motion magnetic resonance imaging method in which an imaged subject moves continuously during data acquisition. Values are determined of scan parameters including at least a velocity, a field of view length in a direction of subject movement, and a preparation sequence time. The imaged subject is continuously moved at the determined velocity. During the continuous moving, magnetic resonance imaging data of the imaged subject are acquired using the determined scan parameters. The acquiring includes performing one or more preparatory operations during one or more time intervals. The acquired data are reconstructed into an image of the field of view. The acquiring and reconstructing are continued during the continuous moving to generate an elongated image spanning a plurality of fields of view.
- a continuous motion magnetic resonance imaging system is disclosed.
- a processor is provided for determining values of scan parameters including at least a subject velocity, a field of view length in the direction of the table velocity, and a preparation sequence time.
- a table is provided for continuously moving a supported imaging subject at the determined velocity.
- a magnetic resonance imaging scanner through which the table continuously moves the subject acquires magnetic resonance imaging data of the imaged subject during the continuous moving using the determined scan parameters.
- the acquiring includes performing one or more preparatory operations during one or more time intervals.
- a reconstruction processor reconstructs the acquired data into an image of the field of view.
- One advantage resides in facilitating the use of continuous moving table magnetic resonance imaging for non-contrast-enhanced imaging, phase-sensitive imaging, and other types of magnetic resonance imaging.
- Another advantage resides in reduced image artifacts due to accumulated phase errors, magnetic field inhomogeneity, and radio frequency field non-uniformity.
- Another advantage is that seamless whole body and large area images are generated.
- Yet another advantage resides in reduced total imaging session time.
- the invention may take form in various components and arrangements of components, and in various process operations and arrangements of process operations.
- the drawings are only for the purpose of illustrating preferred embodiments and are not to be construed as limiting the invention.
- FIG. 1 diagrammatically shows an example continuous moving table magnetic resonance imaging system.
- FIG. 2 diagrammatically shows an example k-space scan performed by the example continuous moving table magnetic resonance imaging system of FIG. 1 .
- FIG. 3 diagrammatically shows a preferred k-space sampling in successive segments of a k-space scan that enhances the effectiveness of an accompanying contrast-manipulating preparatory sequence.
- a magnetic resonance imaging scanner 10 performs magnetic resonance imaging in an examination region 12 .
- the magnetic resonance imaging scanner typically includes main magnetic field coils that produce a substantially temporally and spatially constant Bo magnetic field in the examination region 12 , magnetic field gradient coils for selectively producing magnetic field gradients in the examination region 12 , and one or more radio frequency coils for selectively exciting magnetic resonances in an imaging subject disposed in the examination region 12 .
- a continuous linear translation mechanism 14 moves a pallet or other subject support 16 supporting an associated imaging subject (not shown) in a linear direction, designated the z-direction in FIG. 1 , through the examination region 12 .
- the associated imaging subject is a patient lying prone on the table 16 , and the z-direction of continuous motion corresponds to the axial direction of the patient.
- other imaging subjects and orientations of imaging subjects can be employed.
- the scanner 10 produces a substantially spatially and temporally constant B 0 magnetic field parallel to the z-direction. Although every effort is made to make the B 0 field spatially and temporally constant, there are inhomogeneities in the B 0 field through which the subject passes as the subject moves through the examination region 12 .
- a magnetic resonance imaging controller 20 controls energizing of the magnetic field gradient coils and the radio frequency coil or coils of the scanner 10 to generate and spatially encode magnetic resonance in a portion of the imaging subject disposed in the examination region 12 . This portion disposed in the examination region 12 continuously changes due to the continuous motion of the supporting table 16 .
- a selected k-space trajectory is traversed within the examination region 12 during acquisition of magnetic resonance signals, such as a Cartesian trajectory, a plurality of radial trajectories, or a spiral trajectory.
- a radio frequency coil or coils of the scanner 10 is used for receiving magnetic resonance k-space samples that are stored in a k-space data memory 24 .
- the magnetic resonance imaging scanner is also used, in some embodiments, to acquire one or more spatial maps or calibrations, for example mapping B 0 magnetic field inhomogeneity, radio frequency excitation non-uniformity, a flip angle map, or the like, which are stored in a map/calibration data memory 26 .
- the k-space samples in these embodiments are corrected by a correction processor 30 that corrects the acquired k-space samples for phase errors produced by the continuous motion of the imaging subject through inhomogeneity in the B 0 field, through non-uniformities of the applied radio frequency field, and through applied magnetic field gradients.
- the correction processor 30 receives mapping or calibration data from the map/calibration data memory 26 and also receives the table position as a function of time (denoted “z t ” in FIG. 1 ), and uses this information to associate each acquired k-space sample with the position of the continuously moving associated imaging subject in the inhomogenous magnetic fields, applied magnetic field gradients, and/or non-uniform radio frequency fields at the time the k-space sample was acquired.
- An anatomy alignment processor 34 also receives the table position as a function of time (denoted “z t ” in FIG. 1 ) and uses this information to transform the k-space data into imaging subject coordinates of the continuously moving associated imaging subject.
- a reconstruction processor 40 reconstructs the transformed k-space samples into a reconstructed image that is stored in an images memory 42 , for example by using a Fourier transform based reconstruction algorithm.
- the alignment processor 34 aligns the k-space data with the imaging subject coordinate system by a spatial shift in the z-direction performed after Fourier transform reconstruction, or by a demodulation of the k-space data.
- a user interface 44 enables a radiologist, technician, or other operator to view the reconstructed images, transmit the reconstructed images over a hospital network, the Internet, or another computer network, print the images, or otherwise manipulate the reconstructed images.
- the user interface 44 preferably also interfaces the radiologist, technician, or other operator with the magnetic resonance imaging controller 20 to enable control of the operation of the magnetic resonance imaging scanner 10 .
- a k-space scan 48 images a field of view having a selected length (denoted FOV) in the direction of the table velocity, i.e., the z-direction.
- the length of the field of view in the z-direction generally corresponds to a thick slab, and is less than or equal to the length of the scanner field of view in the z-direction.
- k-space is broken up into a number of segments (the number of segments being denoted M) such as two or more segments.
- segments 1 , 2 , 3 , 4 , . . . M are illustrated.
- Each segment is acquired using a sequence including a preparatory sequence 54 which is performed over a preparatory sequence time (T p ), and a plurality (denoted N M ) of k-space phase encode/sampling cycles 58 each having a duration of a repetition time (TR).
- one or more (denoted N d ) gradient switching startup cycles 62 each of duration TR are performed prior to acquisition of the phase encode/sampling cycles 58 to establish a steady state for the subsequent imaging. (These features are only expressly shown for segment 3 in FIG. 2 .)
- the preparatory sequence 54 has duration T p
- the startup cycles 62 have a duration N d ⁇ TR (where N d can be equal to zero, indicating no startup cycles are employed)
- the k-space phase encode/sampling cycles 58 have a duration N M ⁇ TR.
- N N 1 ⁇ N 2 where N 1 and N 2 are the number of phase encoding lines in the two phase encoding directions, respectively.
- the acquisition time for non-Cartesian k-space trajectories are also readily computed by those skilled in the art.
- the z-direction length of the field of view of the k-space scan can be smaller than the z-direction length of the scanner field of view.
- Using a smaller FOV than the scanner field of view can advantageously reduce the magnitude of discontinuities in the final reconstructed image.
- the startup cycles 62 are advantageous for certain types of imaging techniques, such as fast field echo (FFE) imaging, in which a steady-state should be achieved before acquiring k-space sampling image data.
- FFE fast field echo
- the preparatory operations may be suitably performed only once per field of view (i.e., once per k-space scan), as opposed to being performed for each segment of the k-space scan.
- it may be sufficient to perform field mapping, coil loading, resonance frequency measurements, or the like once in the k-space scan, with the result used to correct data in all the segments of the field of view.
- a contrast-modification preparatory operation such as a slab saturation preparation
- the contrast modification may persist long enough so that a single contrast enhancement preparation is operative for the entire k-space scan.
- the k-space scan 48 diagrammed in FIG. 2 places the preparatory sequence 54 in a single time interval disposed prior to the phase encode/sampling cycles 58 .
- Certain preparatory operations may be or are preferably interleaved with the phase encode steps.
- Some contrast-manipulating preparation pulses, for example, may have short persistence and need to be applied before each phase encode/sampling cycle, once every eight phase encode/sampling cycles, or so forth. Interleaving also advantageously reduces the discontinuity between segments that is present in the example k-space scan 48 due to the preparatory sequence 54 being grouped in a single time interval equal to T p .
- the preparation sequence time (denoted T p ) is computed as the sum of the time intervals over which the preparatory operations are performed.
- Equation (3) typically, many of the parameters of Equation (3) are fixed for a given imaging session.
- the radiologist, technician, or other operator may select a specific FFE imaging sequence having a specific repetition time (TR), a specific number of phase encode steps (N), a specific number of segments (M), and an operator-selected length (FOV) of the field of view in the z-direction.
- the preparatory sequence 54 may be a stored preparation sequence of known duration (T p ) selected by the operator.
- T p the table velocity
- Equation (3) is readily evaluated to determine a suitable table velocity (v) which ensures that each successive acquired field of view is adjacent to and non-overlapping with the next field of view in the imaging subject coordinates system. This value of the table velocity ensures that the continuous moving table imaging provides spatially continuous and non-redundant k-space sampling of the imaging subject over a series of successive fields of view.
- some overlap of the fields of view is chosen, for example by reducing the velocity computed using Equation (3) by a fixed percentage value.
- Voxel values in the overlapping region of adjacent fields of view can be averaged to reduce artifacts due to discontinuities at the boundary between fields of view.
- a parameters processor 70 evaluates Equation (3) to determine the missing parameter value.
- the parameters processor 70 computes the table velocity (v) based on Equation (3) such that the acquired field of view has the desired length (FOV) in the direction of the table velocity.
- Other parameters besides table velocity can be similarly computed, depending upon which values are specified by the user or by recalled stored sequences, and which values remain unknown.
- the preparatory sequence 54 may include applying a contrast-enhancing preparation sequence such as a fat saturation preparation, an outer volume pre-saturation preparation, a T2 preparation, a magnetization transfer contrast enhancement, or the like.
- a contrast-enhancing preparation sequence such as a fat saturation preparation, an outer volume pre-saturation preparation, a T2 preparation, a magnetization transfer contrast enhancement, or the like.
- the duration of each segment is T scan /M; typically, the number of segments (M) is selected such that the segment duration T scan /M is short enough so that the generated contrast enhancement substantially endures throughout the segment duration.
- the effect of the contrast enhancement is optionally further improved by selecting the acquisition of k-space such that the central portion of k-space is acquired first.
- SPIR spectral partial inversion recovery
- the preparatory sequence 54 may additionally or alternatively include other preparatory operations, such as: acquiring Bo magnetic field mapping data; acquiring flip angle mapping data; or acquiring radio frequency mapping data.
- the mapping data can also be used to adjust the scan parameters to adjust for these changes in real time.
- patient-dependent tuned parameters such as: spatially varying subject-induced susceptibility (B 0 , resonance frequency); properties of the send and/or receive radio frequency coils used in the acquisition of magnetic resonance imaging data (such as coil loading, amplitude, phase, flip angle, or so forth); anatomy dependent parameters that may influence the scanning parameters; or so forth. Since the imaging subject is continuously moving during the imaging, these patient-dependent tuned parameters or inhomogeneity maps continuously change during the imaging. Accordingly, the preparatory operations for each k-space scan 48 optionally include mapping or other determination of patient-dependent parameters.
- phase errors caused by main B 0 magnetic field inhomogeneities are corrected by Fourier transforming the data into the image domain (for example, by accessing the Fourier algorithm of the reconstruction processor 40 ) and multiplying the transformed data by a B 0 map acquired either before the start of the scan or during the preparatory sequence 54 .
- phase errors caused by non-uniformities in the radio frequency excitation can be corrected by Fourier transforming the data into the image domain and multiplying the transformed data by a radio frequency map acquired during the preparatory sequence 54 .
- Distortions caused by non-linear magnetic field gradients can also be corrected by the correction processor 30 using techniques such as those described in Glover et al, U.S. Pat. No. 4,591,789.
- the correction of acquired k-space samples can be performed for each acquired k-space line.
- table velocity denoted v
- several lines of k-space can be corrected together. This latter approach recognizes that while the patient-dependent tuned parameters continuously vary as the table 16 continuously moves, those changes can be expected to be slowly varying if the table velocity is small.
- correction of acquired k-space samples is advantageous for substantially any imaging sequence, it is particularly useful in phase-sensitive imaging such as FFE sequences, where accumulated phase errors can be particularly problematic.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Vascular Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
A continuous table motion magnetic resonance imaging system includes a processor (70) for determining values of scan parameters including at least a subject velocity, a field of view length in the direction of the table velocity, and a preparation sequence time. The acquisition includes performing one or more preparatory operations (54) to manipulate image contrast or to update scanning parameters during table motion. A table (16) is provided for continuously moving a supported imaging subject at the determined velocity. A magnetic resonance imaging scanner (10) through which the table continuously moves the subject acquires magnetic resonance imaging data of the imaged subject during the continuous moving using the determined scan parameters. A reconstruction processor (30, 34, 40) reconstructs the acquired data into an image of the field of view.
Description
This application claims the benefit of U.S. provisional application Ser. No. 60/565,881 filed Apr. 28, 2004, which is incorporated herein by reference.
The following relates to the magnetic resonance imaging arts. It finds particular application in cancer screening, angiography, and other diagnostic imaging procedures advantageously performed over the whole body or large portions thereof, and will be described with particular reference thereto. However, it also finds application in continuous moving-table magnetic resonance imaging generally.
Magnetic resonance scanners have a limited field of view that is typically substantially smaller than the average person. Thus, magnetic resonance imaging is not readily used to perform “whole body” or other extended region-of-interest imaging in which the region of interest is larger than the scanner field of view. This limits the usefulness of magnetic resonance imaging for diagnostic tasks such as cancer screening, angiography, or the like.
One approach for magnetic resonance imaging of extended region of interest larger than the scanner field of view is the multi-station method. In this approach, imaging progresses in discrete steps. During each step, one scanner field of view is imaged. Subsequently, the subject is advanced through the scanner bore by a distance equal to the length of the scanner field of view in the axial direction, and another scanner field of view is imaged. The successively imaged fields of view are sewn together to form an image of the extended region of interest.
The multi-station method has certain disadvantages. The start-and-stop motion repeated between each successive station can be disturbing to the patient. Moreover, the imaging must stop between each discrete stepwise advancement of the imaging subject through the scanner bore, which increases the total imaging session time and generates substantial discontinuities at the boundaries between successive scanner fields of view.
Another approach for imaging an extended region of interest larger than the scanner field of view is continuous moving table magnetic resonance imaging. In this approach, the patient is continuously axially advanced through the scanner bore, typically at a constant table velocity, and magnetic resonance imaging data is acquired during the continuous table advancement. This approach may be less disturbing to the patient because it eliminates the discrete start-and-stop table movements of the multi-station method. Moreover, discontinuities between successively imaged fields of view are decreased by omission of the temporally intervening discrete stepwise advancement. Discontinuities can be still further decreased in the continuous moving table approach by using a reduced field of view in the axial direction, that is, by imaging a thick slab field of view which is substantially thinner in the axial direction than the total scanner field of view. It will be appreciated that using a reduced field of view in the multi-station method is generally not advisable since it substantially increases the number of discrete advancement steps or imaged stations.
These advantages of continuous moving table imaging are, however, accompanied by significant disadvantages. The continuously moving imaging subject corresponds to an imaging subject coordinate system which is moving with reference to the magnetic fields and other aspects of the stationary magnetic resonance imaging scanner. This creates difficult timing issues in ensuring collection of a complete, but substantially non-redundant, data set. Moreover, patient-dependent tuned imaging parameters, such as coil loading, resonance frequency, patient-induced susceptibility and flip angle variations typically change during the continuous patient advancement. Still further, each voxel defined with respect to the patient (that is, defined in the imaging subject coordinate system) is continuously moving through magnetic field gradients, the B0 and B1 fields with their inherent inhomogeneities, and the like during the imaging. Such voxel motion can introduce accumulated phase errors or other errors in the collected k-space samples that translate into image artifacts in the reconstructed image.
Heretofore, continuous moving table imaging has been applied mainly to contrast-enhanced angiography. In this application, the progression of an injected magnetic contrast-enhancing bolus is followed through the circulatory system. Continuous table motion facilitates tracking the dynamic bolus advancement. In one particular approach, the contrast agent reduces the T1 value of blood to increase the blood/tissue contrast. Because the contrast is provided principally by the contrast agent, imaging artifacts due to the continuous patient advancement are largely suppressed.
There remains an unfulfilled need in the art for improved apparatuses and methods enabling continuous moving table magnetic resonance imaging to be more generally applied in areas such as non-contrast-enhanced imaging, phase-sensitive imaging, and the like. The present invention contemplates improved apparatuses and methods that overcome the aforementioned limitations and others.
According to one aspect, a continuous motion magnetic resonance imaging method is provided in which an imaged subject moves continuously during data acquisition. Values are determined of scan parameters including at least a velocity, a field of view length in a direction of subject movement, and a preparation sequence time. The imaged subject is continuously moved at the determined velocity. During the continuous moving, magnetic resonance imaging data of the imaged subject are acquired using the determined scan parameters. The acquiring includes performing one or more preparatory operations during one or more time intervals. The acquired data are reconstructed into an image of the field of view. The acquiring and reconstructing are continued during the continuous moving to generate an elongated image spanning a plurality of fields of view.
According to another aspect, a continuous motion magnetic resonance imaging system is disclosed. A processor is provided for determining values of scan parameters including at least a subject velocity, a field of view length in the direction of the table velocity, and a preparation sequence time. A table is provided for continuously moving a supported imaging subject at the determined velocity. A magnetic resonance imaging scanner through which the table continuously moves the subject acquires magnetic resonance imaging data of the imaged subject during the continuous moving using the determined scan parameters. The acquiring includes performing one or more preparatory operations during one or more time intervals. A reconstruction processor reconstructs the acquired data into an image of the field of view.
One advantage resides in facilitating the use of continuous moving table magnetic resonance imaging for non-contrast-enhanced imaging, phase-sensitive imaging, and other types of magnetic resonance imaging.
Another advantage resides in reduced image artifacts due to accumulated phase errors, magnetic field inhomogeneity, and radio frequency field non-uniformity.
Another advantage is that seamless whole body and large area images are generated.
Yet another advantage resides in reduced total imaging session time.
Numerous additional advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments.
The invention may take form in various components and arrangements of components, and in various process operations and arrangements of process operations. The drawings are only for the purpose of illustrating preferred embodiments and are not to be construed as limiting the invention.
With reference to FIG. 1 , a magnetic resonance imaging scanner 10 performs magnetic resonance imaging in an examination region 12. The magnetic resonance imaging scanner typically includes main magnetic field coils that produce a substantially temporally and spatially constant Bo magnetic field in the examination region 12, magnetic field gradient coils for selectively producing magnetic field gradients in the examination region 12, and one or more radio frequency coils for selectively exciting magnetic resonances in an imaging subject disposed in the examination region 12. A continuous linear translation mechanism 14 moves a pallet or other subject support 16 supporting an associated imaging subject (not shown) in a linear direction, designated the z-direction in FIG. 1 , through the examination region 12. In many diagnostic imaging applications, the associated imaging subject is a patient lying prone on the table 16, and the z-direction of continuous motion corresponds to the axial direction of the patient. However, other imaging subjects and orientations of imaging subjects can be employed.
During imaging, the scanner 10 produces a substantially spatially and temporally constant B0 magnetic field parallel to the z-direction. Although every effort is made to make the B0 field spatially and temporally constant, there are inhomogeneities in the B0 field through which the subject passes as the subject moves through the examination region 12. A magnetic resonance imaging controller 20 controls energizing of the magnetic field gradient coils and the radio frequency coil or coils of the scanner 10 to generate and spatially encode magnetic resonance in a portion of the imaging subject disposed in the examination region 12. This portion disposed in the examination region 12 continuously changes due to the continuous motion of the supporting table 16. By applying selected magnetic field gradients, a selected k-space trajectory is traversed within the examination region 12 during acquisition of magnetic resonance signals, such as a Cartesian trajectory, a plurality of radial trajectories, or a spiral trajectory. A radio frequency coil or coils of the scanner 10 is used for receiving magnetic resonance k-space samples that are stored in a k-space data memory 24.
The magnetic resonance imaging scanner is also used, in some embodiments, to acquire one or more spatial maps or calibrations, for example mapping B0 magnetic field inhomogeneity, radio frequency excitation non-uniformity, a flip angle map, or the like, which are stored in a map/calibration data memory 26. The k-space samples in these embodiments are corrected by a correction processor 30 that corrects the acquired k-space samples for phase errors produced by the continuous motion of the imaging subject through inhomogeneity in the B0 field, through non-uniformities of the applied radio frequency field, and through applied magnetic field gradients. The correction processor 30 receives mapping or calibration data from the map/calibration data memory 26 and also receives the table position as a function of time (denoted “zt” in FIG. 1 ), and uses this information to associate each acquired k-space sample with the position of the continuously moving associated imaging subject in the inhomogenous magnetic fields, applied magnetic field gradients, and/or non-uniform radio frequency fields at the time the k-space sample was acquired.
An anatomy alignment processor 34 also receives the table position as a function of time (denoted “zt” in FIG. 1 ) and uses this information to transform the k-space data into imaging subject coordinates of the continuously moving associated imaging subject. A reconstruction processor 40 reconstructs the transformed k-space samples into a reconstructed image that is stored in an images memory 42, for example by using a Fourier transform based reconstruction algorithm. The alignment processor 34 aligns the k-space data with the imaging subject coordinate system by a spatial shift in the z-direction performed after Fourier transform reconstruction, or by a demodulation of the k-space data.
A user interface 44 enables a radiologist, technician, or other operator to view the reconstructed images, transmit the reconstructed images over a hospital network, the Internet, or another computer network, print the images, or otherwise manipulate the reconstructed images. The user interface 44 preferably also interfaces the radiologist, technician, or other operator with the magnetic resonance imaging controller 20 to enable control of the operation of the magnetic resonance imaging scanner 10.
With continuing reference to FIG. 1 and with further reference to FIG. 2 , a k-space scan 48 images a field of view having a selected length (denoted FOV) in the direction of the table velocity, i.e., the z-direction. The length of the field of view in the z-direction generally corresponds to a thick slab, and is less than or equal to the length of the scanner field of view in the z-direction.
Optionally, k-space is broken up into a number of segments (the number of segments being denoted M) such as two or more segments. A single segment (M=1) can also be used (that is, the k-space may optionally not be broken up into a plurality of segments). In example FIG. 2 , segments 1, 2, 3, 4, . . . M are illustrated. Each segment is acquired using a sequence including a preparatory sequence 54 which is performed over a preparatory sequence time (Tp), and a plurality (denoted NM) of k-space phase encode/sampling cycles 58 each having a duration of a repetition time (TR). Optionally, one or more (denoted Nd) gradient switching startup cycles 62 each of duration TR are performed prior to acquisition of the phase encode/sampling cycles 58 to establish a steady state for the subsequent imaging. (These features are only expressly shown for segment 3 in FIG. 2 .)
As labeled in FIG. 2 , the preparatory sequence 54 has duration Tp, the startup cycles 62 have a duration Nd×TR (where Nd can be equal to zero, indicating no startup cycles are employed), and the k-space phase encode/sampling cycles 58 have a duration NM×TR. Assuming no time gaps between the preparatory sequence 54, the startup cycles 62, and the k-space phase encode/sampling cycles 58, the scan time (denoted Tscan) for the k-space scan 48 is given by:
T scan =M×(T p +N d ×TR+N M ×TR) (1).
The total number of phase encode steps (denoted N) is given by the number of phase encode steps in each segment (NM) multiplied by the number of segments (M); that is, N=NM×M. Using this expression, Equation (1) can be rewritten as:
T scan =M×T p+(N+M×N d)×TR (2).
If the field of view is two-dimensional, then the total number of phase encoding steps (denoted N) corresponds to the number of lines of k-space in the two dimensional slice. On the other hand, in three-dimensional imaging, the total number of phase encoding steps is the product of the number of lines in the two phase encoding directions. That is, N=N1×N2 where N1 and N2 are the number of phase encoding lines in the two phase encoding directions, respectively. The acquisition time for non-Cartesian k-space trajectories are also readily computed by those skilled in the art.
T scan =M×(T p +N d ×TR+N M ×TR) (1).
The total number of phase encode steps (denoted N) is given by the number of phase encode steps in each segment (NM) multiplied by the number of segments (M); that is, N=NM×M. Using this expression, Equation (1) can be rewritten as:
T scan =M×T p+(N+M×N d)×TR (2).
If the field of view is two-dimensional, then the total number of phase encoding steps (denoted N) corresponds to the number of lines of k-space in the two dimensional slice. On the other hand, in three-dimensional imaging, the total number of phase encoding steps is the product of the number of lines in the two phase encoding directions. That is, N=N1×N2 where N1 and N2 are the number of phase encoding lines in the two phase encoding directions, respectively. The acquisition time for non-Cartesian k-space trajectories are also readily computed by those skilled in the art.
During the course of the k-space scan 48, the table 16 is continuously moving at a fixed linear velocity (denoted v). The length of the field of view (FOV) of the imaged volume in the z-direction of the table velocity is therefore given by v×Tscan, or equivalently, Tscan=FOV/v. Substituting this expression into Equation (2) and rearranging yields:
FOV=v×[(N+M×N d)×TR+M×T p] (3).
It is to be appreciated that the z-direction length of the field of view of the k-space scan (denoted FOV) can be the same as the z-direction length of the scanner field of view. Alternatively, the z-direction length of the field of view of the k-space scan (denoted FOV) can be smaller than the z-direction length of the scanner field of view. Using a smaller FOV than the scanner field of view can advantageously reduce the magnitude of discontinuities in the final reconstructed image.
FOV=v×[(N+M×N d)×TR+M×T p] (3).
It is to be appreciated that the z-direction length of the field of view of the k-space scan (denoted FOV) can be the same as the z-direction length of the scanner field of view. Alternatively, the z-direction length of the field of view of the k-space scan (denoted FOV) can be smaller than the z-direction length of the scanner field of view. Using a smaller FOV than the scanner field of view can advantageously reduce the magnitude of discontinuities in the final reconstructed image.
The startup cycles 62 are advantageous for certain types of imaging techniques, such as fast field echo (FFE) imaging, in which a steady-state should be achieved before acquiring k-space sampling image data. In other types of imaging in which achievement of a steady state or quasi-steady state condition is not advantageous, the startup cycles are suitably omitted (Nd=0), in which case Equation (3) simplifies to:
FOV=v×[N×TR+M×T p] (4).
FOV=v×[N×TR+M×T p] (4).
In some embodiments, the preparatory operations, or some of the preparatory operations, may be suitably performed only once per field of view (i.e., once per k-space scan), as opposed to being performed for each segment of the k-space scan. For example, it may be sufficient to perform field mapping, coil loading, resonance frequency measurements, or the like once in the k-space scan, with the result used to correct data in all the segments of the field of view. Similarly, in a contrast-modification preparatory operation such as a slab saturation preparation, the contrast modification may persist long enough so that a single contrast enhancement preparation is operative for the entire k-space scan.
In the most general case, there are one or more preparatory operations performed once per field of view, and one or more other preparatory operations performed for each segment. Denoting the time per field of view for the once-per-field of view preparatory operation or operations by the time Tp,1, and the time per segment for the once-per-segment operation or operations by the time Tp,s, the average preparation time (denoted Tp and used in Equations (1)-(4)) is given by:
T p=(T p,1 /M)+T p,s (5).
If all preparatory operations are performed on a per-segment basis, Tp,1=0 and Tp=Tp,s. On the other hand, if all preparatory operations are performed on a per-field of view basis, Tp,s=0 and Tp=Tp,1/M.
T p=(T p,1 /M)+T p,s (5).
If all preparatory operations are performed on a per-segment basis, Tp,1=0 and Tp=Tp,s. On the other hand, if all preparatory operations are performed on a per-field of view basis, Tp,s=0 and Tp=Tp,1/M.
The k-space scan 48 diagrammed in FIG. 2 places the preparatory sequence 54 in a single time interval disposed prior to the phase encode/sampling cycles 58. Certain preparatory operations, however, may be or are preferably interleaved with the phase encode steps. Some contrast-manipulating preparation pulses, for example, may have short persistence and need to be applied before each phase encode/sampling cycle, once every eight phase encode/sampling cycles, or so forth. Interleaving also advantageously reduces the discontinuity between segments that is present in the example k-space scan 48 due to the preparatory sequence 54 being grouped in a single time interval equal to Tp. When the preparatory sequence is partially or wholly distributed or interleaved amongst the phase encode/sampling steps, the preparation sequence time (denoted Tp) is computed as the sum of the time intervals over which the preparatory operations are performed.
Typically, many of the parameters of Equation (3) are fixed for a given imaging session. For example, the radiologist, technician, or other operator may select a specific FFE imaging sequence having a specific repetition time (TR), a specific number of phase encode steps (N), a specific number of segments (M), and an operator-selected length (FOV) of the field of view in the z-direction. Similarly, the preparatory sequence 54 may be a stored preparation sequence of known duration (Tp) selected by the operator. In this example, only the table velocity is left unspecified. Accordingly, Equation (3) is readily evaluated to determine a suitable table velocity (v) which ensures that each successive acquired field of view is adjacent to and non-overlapping with the next field of view in the imaging subject coordinates system. This value of the table velocity ensures that the continuous moving table imaging provides spatially continuous and non-redundant k-space sampling of the imaging subject over a series of successive fields of view.
Optionally, some overlap of the fields of view is chosen, for example by reducing the velocity computed using Equation (3) by a fixed percentage value. Voxel values in the overlapping region of adjacent fields of view can be averaged to reduce artifacts due to discontinuities at the boundary between fields of view.
With reference to FIG. 1 , more generally a parameters processor 70 evaluates Equation (3) to determine the missing parameter value. Thus, for example, if all scan parameters except the table velocity (v) are known based on the type of scan to be performed and other information, then the parameters processor 70 computes the table velocity (v) based on Equation (3) such that the acquired field of view has the desired length (FOV) in the direction of the table velocity. Other parameters besides table velocity can be similarly computed, depending upon which values are specified by the user or by recalled stored sequences, and which values remain unknown.
Moreover, if more than one scan parameter remains unknown, Equation (3) can be evaluated to determine suitable combinations of the plurality of unknown scan parameters. The radiologist, technician, or other operator can then select from the available combinations. For example, if the table velocity and field of view are both unknown, Equation (3) can be solved for the ratio FOV/v according to:
FOV/v=(N+M×N d)×TR+M×T p (6).
The user can then be presented with a set of proposed pairs of values for FOV and v each of which satisfy the computed FOV/v ratio. The user can then select one of the proposed pairs which provides a reasonable table velocity and field of view.
FOV/v=(N+M×N d)×TR+M×T p (6).
The user can then be presented with a set of proposed pairs of values for FOV and v each of which satisfy the computed FOV/v ratio. The user can then select one of the proposed pairs which provides a reasonable table velocity and field of view.
The preparatory sequence 54 may include applying a contrast-enhancing preparation sequence such as a fat saturation preparation, an outer volume pre-saturation preparation, a T2 preparation, a magnetization transfer contrast enhancement, or the like. The duration of each segment is Tscan/M; typically, the number of segments (M) is selected such that the segment duration Tscan/M is short enough so that the generated contrast enhancement substantially endures throughout the segment duration.
With reference to FIG. 3 , the effect of the contrast enhancement is optionally further improved by selecting the acquisition of k-space such that the central portion of k-space is acquired first. FIG. 3 plots the values of k-space samples acquired in two successive segments. In the first segment, positive k-space values are acquired starting at k=0 and increasing until k=kmax. In the second segment, negative k-space values are acquired starting at k=0 and decreasing until k=−kmax. The preparatory sequences (denoted by “(1)” in FIG. 3 ) immediately precede the k-space sampling; hence the data around k=0 is acquired immediately after each preparatory sequence. In one example, the preparatory sequences (denoted “(1)” in FIG. 3 ) each include a T2 preparation followed by a fat suppression pulse using spectral partial inversion recovery (SPIR), followed by two startup cycles.
The preparatory sequence 54 may additionally or alternatively include other preparatory operations, such as: acquiring Bo magnetic field mapping data; acquiring flip angle mapping data; or acquiring radio frequency mapping data. Such optional mapping data update or are stored in the map/calibration data memory 26 (see FIG. 1 ), and if acquired are used by the correction processor 30 to correct the k-space samples for various scanner inhomogeneities or phase errors. The mapping data can also be used to adjust the scan parameters to adjust for these changes in real time. Many of these inhomogeneities are reflected in patient-dependent tuned parameters, such as: spatially varying subject-induced susceptibility (B0, resonance frequency); properties of the send and/or receive radio frequency coils used in the acquisition of magnetic resonance imaging data (such as coil loading, amplitude, phase, flip angle, or so forth); anatomy dependent parameters that may influence the scanning parameters; or so forth. Since the imaging subject is continuously moving during the imaging, these patient-dependent tuned parameters or inhomogeneity maps continuously change during the imaging. Accordingly, the preparatory operations for each k-space scan 48 optionally include mapping or other determination of patient-dependent parameters.
In one example of a correction suitably performed by the correction processor 30, phase errors caused by main B0 magnetic field inhomogeneities are corrected by Fourier transforming the data into the image domain (for example, by accessing the Fourier algorithm of the reconstruction processor 40) and multiplying the transformed data by a B0 map acquired either before the start of the scan or during the preparatory sequence 54. Similarly, phase errors caused by non-uniformities in the radio frequency excitation can be corrected by Fourier transforming the data into the image domain and multiplying the transformed data by a radio frequency map acquired during the preparatory sequence 54. Distortions caused by non-linear magnetic field gradients can also be corrected by the correction processor 30 using techniques such as those described in Glover et al, U.S. Pat. No. 4,591,789.
The correction of acquired k-space samples can be performed for each acquired k-space line. Alternatively, if the table velocity (denoted v) is small enough, several lines of k-space can be corrected together. This latter approach recognizes that while the patient-dependent tuned parameters continuously vary as the table 16 continuously moves, those changes can be expected to be slowly varying if the table velocity is small. While correction of acquired k-space samples is advantageous for substantially any imaging sequence, it is particularly useful in phase-sensitive imaging such as FFE sequences, where accumulated phase errors can be particularly problematic.
The invention has been described with reference to the preferred embodiments. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (20)
1. A continuous motion magnetic resonance imaging method in which an imaged subject moves continuously during data acquisition, the method comprising:
determining values of scan parameters including at least a velocity, a field of view length in a direction of subject movement, and a preparation sequence time;
continuously moving the imaged subject at the determined velocity;
during the continuous moving, acquiring magnetic resonance imaging data of the imaged subject using the determined scan parameters, the acquiring including performing one or more preparatory operations during one or more time intervals;
correcting the acquired magnetic resonance imaging data for a spatial non-uniformity based on a non-uniformity mapping performed as one of the one or more preparatory operations;
reconstructing the acquired and corrected data into an image of the field of view; and
continuing the acquiring, correcting, and reconstructing during the continuous moving to generate an elongated image spanning a plurality of fields of view.
2. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the correcting comprises:
correcting the acquired magnetic resonance imaging data for variation of a tuned parameter due to the continuous moving.
3. The continuous motion magnetic resonance imaging method as set forth in claim 2 , wherein the tuned parameter is selected from a group consisting of: (i) spatial variation of the B0 magnetic field, (ii) spatial variation of the resonance frequency, (iii) a property of a radio frequency coil used in the acquisition of magnetic resonance imaging data, and (iv) anatomy dependent parameters influencing one or more of the scan parameters.
4. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the non-uniformity mapping includes measuring at least one of magnetic field and radio frequency inhomogeneity, and the correcting includes correcting the acquired data in accordance with the measured inhomogeneities.
5. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the determining of values of the scan parameters comprises:
computing a value for the velocity based on values of the other scan parameters.
6. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the determining of values of the scan parameters comprises:
computing a value of a selected one of the scan parameters based on values of the other scan parameters.
7. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the one or more preparatory operations performed during the one or more time intervals include at least one of:
applying a contrast enhancing preparation sequence,
acquiring a B0 magnetic field map,
acquiring a flip angle map, and
acquiring a radio frequency map.
8. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the one or more preparatory operations performed during the one or more time intervals include:
applying a contrast manipulating preparation sequence, the acquiring of magnetic resonance imaging data not involving an administered contrast agent.
9. The continuous motion magnetic resonance imaging method as set forth in claim 1 , wherein the reconstructing comprises:
transforming the acquired data into subject coordinates based on position of the continuously moving imaged subject as a function of time during the acquiring; and
reconstructing the transformed acquired data.
10. The continuous motion magnetic resonance imaging method as set forth in claim 9 , wherein the determining of values of the scan parameters comprises:
determining values of the scan parameters such that successive reconstructed images of fields of view are spatially continuous in the direction of subject movement.
11. The continuous motion magnetic resonance imaging method as set forth in claim 9 , wherein the determining of values of the scan parameters comprises:
determining values of the scan parameters such that successive reconstructed images of fields of view are spatially adjacent and non-overlapping in the direction of subject movement.
12. The continuous motion magnetic resonance imaging method as set forth in claim 9 , wherein the one or more time intervals during which the one or more preparatory operations are performed sum to the determined preparation time.
13. A continuous motion magnetic resonance imaging method in which an image subject moves continuously during data acquisition, the method comprising:
continuously moving the imaged subject at a determined velocity;
during the continuous moving, acquiring magnetic resonance imaging data of the imaged subject;
correcting the acquired magnetic resonance imaging data for phase error introduced by traversal of the imaged subject through an inhomogeneous magnetic field during the continuous moving;
reconstructing the acquired and corrected data into an image of the field of view; and
continuing the acquiring and reconstructing during the continuous moving to generate an elongated image spanning a plurality of fields of view.
14. The continuous motion magnetic resonance imaging method as set forth in claim 13 , further comprising:
performing an initial pre-scan including continuously moving the imaged subject and acquiring mapping data, the pre-scan being performed prior to the acquiring of magnetic resonance imaging data;
computing an imaged subject position dependent calibration based on the acquired mapping data; and
wherein the correcting of the acquired magnetic resonance imaging data for phase error based on the imaged subject position dependent calibration.
15. The continuous motion magnetic resonance imaging method as set forth in claim 13 , wherein the correcting comprises:
transforming the acquired magnetic resonance imaging data into the image domain; and
multiplying the transformed data by a B0 map.
16. The continuous motion magnetic resonance imaging method as set forth in claim 13 , wherein the correcting comprises:
transforming the acquired magnetic resonance imaging data into the image domain; and
multiplying the transformed data by a radio frequency map.
17. A continuous motion magnetic resonance imaging system comprising:
a processor for determining values of scan parameters including at least a subject velocity, a field of view length in the direction of the table velocity, and a preparation sequence time;
a table for continuously moving a supported imaging subject at the determined velocity;
a magnetic resonance imaging scanner through which the table continuously moves the subject, the scanner acquiring magnetic resonance imaging data of the imaged subject during the continuous moving using the determined scan parameters, the acquiring including performing one or more preparatory operations during one or more time intervals; and
a reconstruction processor for reconstructing the acquired data into an image of the field of view, the reconstruction processor including a correction processor for correcting the acquired data for radio frequency excitation inhomogeneity, the reconstruction processor reconstructing the corrected acquired data.
18. A continuous motion magnetic resonance imaging system comprising:
a processor for determining values of scan parameters including at least a subject velocity, a field of view length in the direction of the table velocity, and a preparation sequence time;
a table for continuously moving a supported imaging subject at the determined velocity;
a magnetic resonance imaging scanner through which the table continuously moves the subject, the scanner acquiring magnetic resonance imaging data of the imaged subject during the continuous moving using the determined scan parameters, the acquiring including performing one or more preparatory operations during one or more time intervals; and
a reconstruction processor for reconstructing the acquired data into an image of the field of view, the reconstruction processor including a correction processor for correcting the acquired data for magnetic field inhomogeneity, the reconstruction processor reconstructing the corrected acquired data.
19. The continuous motion magnetic resonance imaging system as set forth in claim 18 , wherein the reconstruction processor further comprises:
an anatomy alignment processor for transforming the acquired data into imaging subject coordinates based on position of the continuously moving table as a function of time during the acquiring.
20. The continuous motion magnetic resonance imaging system as set forth in claim 18 , wherein the correction processor corrects the acquired data for motion induced accumulated phase error due to a magnetic field gradient.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US56588104P | 2004-04-28 | 2004-04-28 | |
PCT/IB2005/051071 WO2005106522A1 (en) | 2004-04-28 | 2005-03-30 | Continuous moving-table mri involving contrast manipulation and/or update of scanning parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
US20070222442A1 US20070222442A1 (en) | 2007-09-27 |
US7417430B2 true US7417430B2 (en) | 2008-08-26 |
Family
ID=34963825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/568,321 Expired - Fee Related US7417430B2 (en) | 2004-04-28 | 2005-03-30 | Continuous moving-table MRI contrast manipulation and/or update of scanning parameters |
Country Status (5)
Country | Link |
---|---|
US (1) | US7417430B2 (en) |
EP (1) | EP1743188A1 (en) |
JP (1) | JP2007534414A (en) |
CN (1) | CN1977182A (en) |
WO (1) | WO2005106522A1 (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060173277A1 (en) * | 2005-02-03 | 2006-08-03 | Daniel Elgort | Adaptive imaging parameters with MRI |
US20090066328A1 (en) * | 2005-04-11 | 2009-03-12 | Tetsuhiko Takahashi | Nuclear Magnetic Resonance Imaging Apparatus and Method |
US20090219020A1 (en) * | 2005-04-28 | 2009-09-03 | Shinji Kurokawa | Magnetic resonance imaging apparatus |
US20100001727A1 (en) * | 2005-12-21 | 2010-01-07 | Lucio Frydman | Method and apparatus for acquiring high resolution spectral data or high definition images in inhomogeneous environments |
US20100079141A1 (en) * | 2008-09-08 | 2010-04-01 | Alto Stemmer | Magnetic resonance system and method for correction of distortion due to continuous movement of patient table |
US20100156413A1 (en) * | 2008-12-19 | 2010-06-24 | Sloan Kettering Institute For Cancer Research | Corrected nuclear magnetic resonance imaging using magnetization transfer |
US20100237863A1 (en) * | 2009-03-19 | 2010-09-23 | Alto Stemmer | Method and device for controlling acquisition of magnetic resonance data in a magnetic resonance apparatus |
US20100244823A1 (en) * | 2006-03-31 | 2010-09-30 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and method |
DE102010038777A1 (en) * | 2010-08-02 | 2012-02-02 | Siemens Aktiengesellschaft | Method for creating MR images and corresponding magnetic resonance system |
US20120032676A1 (en) * | 2010-08-06 | 2012-02-09 | Toshiba Medical Systems Corporation | Spatial intensity correction for RF shading non-uniformities in MRI |
US20120133360A1 (en) * | 2010-11-26 | 2012-05-31 | Kensuke Shinoda | Magnetic resonance diagnosis apparatus and data acquisition method of magnetic resonance spectroscopy |
US20130320979A1 (en) * | 2011-02-15 | 2013-12-05 | Koninklijke Philips N.V. | Magnetic resonance spectroscopy with automatic phase and b0 correction using interleaved water reference scan |
US8810242B2 (en) | 2010-08-06 | 2014-08-19 | Kabushiki Kaisha Toshiba | Spatial intensity correction for RF shading non-uniformities in MRI |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006028015A1 (en) * | 2004-09-06 | 2006-03-16 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and method |
US20060184027A1 (en) * | 2005-01-06 | 2006-08-17 | Kabushiki Kaisha Toshiba | Diagnostic imaging system, magnetic resonance imaging apparatus, and method of diagnostic imaging |
WO2006106448A2 (en) * | 2005-04-06 | 2006-10-12 | Koninklijke Philips Electronics, N.V. | Sense mr parallel imaging with continuously moving bed |
US8320647B2 (en) | 2007-11-20 | 2012-11-27 | Olea Medical | Method and system for processing multiple series of biological images obtained from a patient |
WO2009074917A1 (en) * | 2007-12-11 | 2009-06-18 | Koninklijke Philips Electronics N.V. | Reducing motion artefacts in mri |
DE102008028919B4 (en) * | 2008-06-18 | 2012-03-29 | Siemens Aktiengesellschaft | Contacting system and method for contacting magnetic resonance local coils with a signal processing unit of a magnetic resonance tomograph and magnetic resonance tomograph |
CN102084264B (en) * | 2008-07-04 | 2016-01-20 | 皇家飞利浦电子股份有限公司 | MR imaging with extended field of view |
DE102009014903B4 (en) * | 2009-03-25 | 2012-01-12 | Siemens Aktiengesellschaft | Magnetic resonance projection angiography with continuous table displacement |
DE102009017778A1 (en) * | 2009-04-20 | 2010-11-18 | Siemens Aktiengesellschaft | Method, magnetic resonance apparatus and computer program for acquisition of measurement data of an examination area of an examination subject during a continuous movement of the examination area through a magnetic resonance apparatus |
US8600476B2 (en) | 2011-04-21 | 2013-12-03 | Siemens Aktiengesellschaft | Patient support table control system for use in MR imaging |
US9295406B2 (en) * | 2011-05-05 | 2016-03-29 | Siemens Medical Solutions Usa, Inc. | Automatic or semi-automatic whole body MR scanning system |
RU2015110990A (en) * | 2012-08-27 | 2016-10-20 | Конинклейке Филипс Н.В. | MOVEMENT TRACKING BASED ON QUICK IMAGE RECEIPT |
CN104161517B (en) * | 2014-06-12 | 2016-06-15 | 中国科学院苏州生物医学工程技术研究所 | A kind of formation method strengthened for MRI contrast |
US10456103B2 (en) | 2016-06-14 | 2019-10-29 | Shanghai United Imaging Healthcare Co., Ltd. | Unified trajectory generation process and system |
US11360175B2 (en) * | 2017-06-28 | 2022-06-14 | Canon Medical Systems Corporation | Magnetic resonance imaging apparatus |
WO2019127443A1 (en) * | 2017-12-29 | 2019-07-04 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for synchronous motion optimization of device with moving components |
CN112345990B (en) * | 2020-10-22 | 2022-02-08 | 上海交通大学 | Three-echo balance steady-state free precession pulse sequence imaging method |
CN113219390A (en) * | 2021-05-21 | 2021-08-06 | 电子科技大学 | Method for dynamically updating scanning parameters through magnetic resonance |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4591789A (en) | 1983-12-23 | 1986-05-27 | General Electric Company | Method for correcting image distortion due to gradient nonuniformity |
US6385478B1 (en) | 1998-12-24 | 2002-05-07 | Picker International, Inc. | Magnetic resonance imaging apparatus |
US6408201B1 (en) * | 2000-06-09 | 2002-06-18 | General Electric Company | Method and apparatus for efficient stenosis identification in peripheral arterial vasculature using MR imaging |
US20020140423A1 (en) | 2001-03-30 | 2002-10-03 | Brittain Jean Helen | Moving table MRI with frequency-encoding in the z-direction |
US20020173715A1 (en) | 2001-04-09 | 2002-11-21 | Kruger David G. | Method for acquiring MRI data from a large field of view using continuous table motion |
US20030011369A1 (en) | 2001-03-30 | 2003-01-16 | Brittain Jean H. | Moving table MRI with frequency-encoding in the z-direction |
US6671536B2 (en) | 2001-09-25 | 2003-12-30 | Wisconsin Alumni Research Foundation | Magnetic resonance angiography using floating table projection imaging |
EP1398642A2 (en) | 2002-09-12 | 2004-03-17 | General Electric Company | Extended volume imaging using MRI with parallel reception |
US6963768B2 (en) * | 2002-05-16 | 2005-11-08 | General Electric Company | Whole body MRI scanning with moving table and interactive control |
US7251520B2 (en) * | 2003-07-08 | 2007-07-31 | General Electric Company | Method and apparatus of slice selective magnetization preparation for moving table MRI |
-
2005
- 2005-03-30 WO PCT/IB2005/051071 patent/WO2005106522A1/en not_active Application Discontinuation
- 2005-03-30 US US11/568,321 patent/US7417430B2/en not_active Expired - Fee Related
- 2005-03-30 CN CNA2005800135474A patent/CN1977182A/en active Pending
- 2005-03-30 EP EP05718598A patent/EP1743188A1/en not_active Withdrawn
- 2005-03-30 JP JP2007510168A patent/JP2007534414A/en not_active Withdrawn
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4591789A (en) | 1983-12-23 | 1986-05-27 | General Electric Company | Method for correcting image distortion due to gradient nonuniformity |
US6385478B1 (en) | 1998-12-24 | 2002-05-07 | Picker International, Inc. | Magnetic resonance imaging apparatus |
US6408201B1 (en) * | 2000-06-09 | 2002-06-18 | General Electric Company | Method and apparatus for efficient stenosis identification in peripheral arterial vasculature using MR imaging |
US20020140423A1 (en) | 2001-03-30 | 2002-10-03 | Brittain Jean Helen | Moving table MRI with frequency-encoding in the z-direction |
US20030011369A1 (en) | 2001-03-30 | 2003-01-16 | Brittain Jean H. | Moving table MRI with frequency-encoding in the z-direction |
US20020173715A1 (en) | 2001-04-09 | 2002-11-21 | Kruger David G. | Method for acquiring MRI data from a large field of view using continuous table motion |
US6671536B2 (en) | 2001-09-25 | 2003-12-30 | Wisconsin Alumni Research Foundation | Magnetic resonance angiography using floating table projection imaging |
US6963768B2 (en) * | 2002-05-16 | 2005-11-08 | General Electric Company | Whole body MRI scanning with moving table and interactive control |
EP1398642A2 (en) | 2002-09-12 | 2004-03-17 | General Electric Company | Extended volume imaging using MRI with parallel reception |
US7251520B2 (en) * | 2003-07-08 | 2007-07-31 | General Electric Company | Method and apparatus of slice selective magnetization preparation for moving table MRI |
Non-Patent Citations (3)
Title |
---|
Aldefeld, B., et al.; Sensitivity to RF and Bo Field Imperfections in Continuous Moving-Table MR Imaging; 2004; Proc. Int'l. Soc. Mag. Reson. Med.; 11; p. 102. |
Shankaranayanan, A., et al., Helical MR: Continuously Moving Table Axial Imaging with Radical Acquisitions; 2003; MRM; 50:1053-1060. |
Zhu, Y., et al., Extended Field-of-View Imaging with Table Translation and Frequency Sweeping; 2003; MRM; 49:1106-1112. |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7643864B2 (en) * | 2005-02-03 | 2010-01-05 | Case Western Reserve University | Adaptive imaging parameters with MRI |
US20060173277A1 (en) * | 2005-02-03 | 2006-08-03 | Daniel Elgort | Adaptive imaging parameters with MRI |
US20090066328A1 (en) * | 2005-04-11 | 2009-03-12 | Tetsuhiko Takahashi | Nuclear Magnetic Resonance Imaging Apparatus and Method |
US7635978B2 (en) * | 2005-04-11 | 2009-12-22 | Hitachi Medical Corporation | Nuclear magnetic resonance imaging apparatus and method |
US20090219020A1 (en) * | 2005-04-28 | 2009-09-03 | Shinji Kurokawa | Magnetic resonance imaging apparatus |
US7701214B2 (en) * | 2005-04-28 | 2010-04-20 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus |
US20100001727A1 (en) * | 2005-12-21 | 2010-01-07 | Lucio Frydman | Method and apparatus for acquiring high resolution spectral data or high definition images in inhomogeneous environments |
US7944206B2 (en) * | 2005-12-21 | 2011-05-17 | Yeda Research And Development Co. Ltd. | Method and apparatus for acquiring high resolution spectral data or high definition images in inhomogeneous environments |
US20100244823A1 (en) * | 2006-03-31 | 2010-09-30 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and method |
US7944209B2 (en) * | 2006-03-31 | 2011-05-17 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and method |
US20100079141A1 (en) * | 2008-09-08 | 2010-04-01 | Alto Stemmer | Magnetic resonance system and method for correction of distortion due to continuous movement of patient table |
US8294462B2 (en) * | 2008-09-08 | 2012-10-23 | Siemens Aktiengesellschaft | Magnetic resonance system and method for correction of distortion due to continuous movement of patient table |
US8405392B2 (en) * | 2008-12-19 | 2013-03-26 | Sloan Kettering Institute For Cancer Research | Corrected nuclear magnetic resonance imaging using magnetization transfer |
US20100156413A1 (en) * | 2008-12-19 | 2010-06-24 | Sloan Kettering Institute For Cancer Research | Corrected nuclear magnetic resonance imaging using magnetization transfer |
US20100237863A1 (en) * | 2009-03-19 | 2010-09-23 | Alto Stemmer | Method and device for controlling acquisition of magnetic resonance data in a magnetic resonance apparatus |
US8604786B2 (en) * | 2009-03-19 | 2013-12-10 | Siemens Aktiengesellschaft | Method and device for controlling acquisition of magnetic resonance data in a magnetic resonance apparatus |
DE102010038777A1 (en) * | 2010-08-02 | 2012-02-02 | Siemens Aktiengesellschaft | Method for creating MR images and corresponding magnetic resonance system |
US8643365B2 (en) | 2010-08-02 | 2014-02-04 | Siemens Aktiengesellschaft | Method and magnetic resonance system to generate magnetic resonance images |
DE102010038777B4 (en) * | 2010-08-02 | 2012-12-13 | Siemens Aktiengesellschaft | Creation of MR images of a predetermined volume section within an examination object with continuous table displacement |
US20120032676A1 (en) * | 2010-08-06 | 2012-02-09 | Toshiba Medical Systems Corporation | Spatial intensity correction for RF shading non-uniformities in MRI |
US8217652B2 (en) * | 2010-08-06 | 2012-07-10 | Kabushiki Kaisha Toshiba | Spatial intensity correction for RF shading non-uniformities in MRI |
US8810242B2 (en) | 2010-08-06 | 2014-08-19 | Kabushiki Kaisha Toshiba | Spatial intensity correction for RF shading non-uniformities in MRI |
US20120133360A1 (en) * | 2010-11-26 | 2012-05-31 | Kensuke Shinoda | Magnetic resonance diagnosis apparatus and data acquisition method of magnetic resonance spectroscopy |
US8890521B2 (en) * | 2010-11-26 | 2014-11-18 | Kabushiki Kaisha Toshiba | Magnetic resonance diagnosis apparatus and data acquisition method of magnetic resonance spectroscopy |
US20130320979A1 (en) * | 2011-02-15 | 2013-12-05 | Koninklijke Philips N.V. | Magnetic resonance spectroscopy with automatic phase and b0 correction using interleaved water reference scan |
US9470770B2 (en) * | 2011-02-15 | 2016-10-18 | Koninklijke Philips N.V. | Magnetic resonance spectroscopy with automatic phase and B0 correction using interleaved water reference scan |
Also Published As
Publication number | Publication date |
---|---|
EP1743188A1 (en) | 2007-01-17 |
JP2007534414A (en) | 2007-11-29 |
US20070222442A1 (en) | 2007-09-27 |
WO2005106522A1 (en) | 2005-11-10 |
WO2005106522A8 (en) | 2006-12-14 |
CN1977182A (en) | 2007-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7417430B2 (en) | Continuous moving-table MRI contrast manipulation and/or update of scanning parameters | |
JP4646015B2 (en) | Method of operating a magnetic resonance imaging (MRI) system | |
JP4166068B2 (en) | Table-moving MRI with frequency encoding in the z direction | |
US9687172B2 (en) | System for motion corrected MR diffusion imaging | |
JP4610611B2 (en) | Magnetic resonance imaging device | |
US6037771A (en) | Sliding thin-slab acquisition of three-dimensional MRA data | |
US8185188B2 (en) | Method and apparatus for magnetic resonance imaging on the basis of a gradient echo sequence | |
EP1472978B1 (en) | Magnetic resonance imaging method and apparatus | |
JP2004527301A (en) | Method of acquiring MRI data from a large field of view using continuous table movement | |
EP1615044A1 (en) | Continuous table motion MRI involving variable fields of view during the scan | |
Scott et al. | Beat-to-beat respiratory motion correction with near 100% efficiency: a quantitative assessment using high-resolution coronary artery imaging | |
JPS63300750A (en) | Method for reducing artificial effect of image | |
EP1693680B1 (en) | Continuous table motion MRI involving phase correction | |
JP2002530173A (en) | Magnetic resonance method and apparatus | |
JP4060459B2 (en) | Magnetic resonance imaging system | |
US11543482B2 (en) | Magnetic resonance imaging using motion-compensated image reconstruction | |
JP4330247B2 (en) | Nuclear magnetic resonance imaging system | |
JP2004089275A (en) | Phase correction method in magnetic resonance imaging device | |
JP3907944B2 (en) | Magnetic resonance imaging method and apparatus | |
JP4678916B2 (en) | Magnetic resonance imaging system | |
EP0905525A1 (en) | Sliding thin-slab aquisition of three-dimensional MRA data | |
JP4832510B2 (en) | Magnetic resonance imaging device | |
JP2004329269A (en) | Magnetic resonance imaging apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: KONINKLIJKE PHILIPS ELECTRONICS N.V., NETHERLANDS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ALDEFELD, BERND;BOERNERT, PETER;KEUPP, JOCHEN;REEL/FRAME:018437/0812 Effective date: 20040427 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20120826 |