The present invention relates generally to an improved method of medical imaging over large areas, and more particularly, to a method and apparatus of acquiring magnetic resonance (MR) images over an area that is greater than the optimal imaging area of an MR scanner using continuous table movement through the MR scanner without incurring slab-boundary artifacts.
When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field B0), the individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but process about it in random order at their characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic field (excitation field B1) which is in the x-y plane and which is near the Larmor frequency, the net aligned moment, or “longitudinal magnetization”, Mz, may be rotated, or “tipped”, into the x-y plane to produce a net transverse magnetic moment Mt. A signal is emitted by the excited spins after the excitation signal B, is terminated and this signal may be received and processed to form an image.
When utilizing these signals to produce images, magnetic field gradients (Gx Gy and Gz) are employed. Typically, the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used. The resulting set of received MR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.
In such MRI systems, the volume for acquiring MR data with optimal gradient linearity, having a uniform magnetic field Bo, and uniform radio frequency (RF) homogeneity is of limited extent. Desired fields-of-view (FOV) that exceed this limited volume are traditionally acquired in sections, with table motion between scans. The resulting concatenated images often exhibit discontinuities at the slab junctions. These slab-boundary artifacts result in non-ideal images. When these artifacts are either severe or occur in a critical region-of-interest, complete re-acquisition of data may be needed for a thorough analysis. In addition, the time for table motion extends scan time beyond that required for data acquisition.
Known methods designed to eliminate slab boundary artifacts in angiograms and 3D fast spin echo acquisitions include Sliding Interleaved ky Acquisition (SLINKY) and Shifted Interleaved Multi-Volume Acquisition (SIMVA). These methods however, move the slab-selective excitation while keeping the table stationary rather than moving the table and keeping the slab position stationary. As a result, these methods are limited by the inherent optimal imaging volume of the MR scanner. These methods employ phase encoding in the z-direction and do not account for those situations where the z matrix is not equal to the number of kx-ky subsets. Similarly, other known methods designed to eliminate slab boundary artifacts have likewise moved the imaging slab rather than the patient table. Other known techniques implement a continuously moving table in both the phase encoding and frequency direction. However, prior techniques that move the patient table in the frequency-encode direction encode the data as if they were acquiring an image of the entire FOV and data is combined prior to any Fourier transform.
Other known systems employ stepped table and/or moving table approach with an array of receiver coils that move with the imaged object. Data is collected from each coil independently as it moves through the homogeneous volume of the scanner. None of these known systems however, collect data from a slab thickness that is fixed relative to the magnet, place the frequency encoding axis in the direction of table motion, and combine the data after Fourier transforming in the direction of table motion.
It would therefore be desirable to have a new method and apparatus that allows coverage of large FOV without slab-boundary artifacts in the resulting concatenated images.
The present invention relates to a system and method of acquiring large FOV MR images using continuous table motion for increased volume coverage that results in reconstructed images without discontinuities.
Slab-boundary artifacts are eliminated and scan time is reduced in an acquisition sequence by continuously moving an imaging object with respect to the optimal volume of the imaging apparatus, or vice versa. The thickness of the slab, which is smaller than the desired FOV, is selected to remain within the optimal volume of the MR system. The selected slab position remains fixed relative to the magnet of the MR scanner during the scan, and the table is moved continuously during scanning of the entire FOV. MR data is acquired by applying an excitation that excites spins and applying magnetic field gradient waveforms to encode the volume of interest. The volume of interest is restricted in the direction of table motion. The magnetic field gradients traverse k-space following a trajectory that is uniform in the k-space dimension that is in the direction of table motion. The magnetic field gradient waveforms that encode the k-space directions perpendicular to table motion are divided into subsets. During each acquisition, all the k-space data in the direction of table motion are acquired for a subset of at least one other k-space dimension. After acquisition, data is transformed in the direction of table motion, sorted, and aligned to match anatomical locations in the direction of table motion. This procedure is repeated to fill the entire matrix. A final image is reconstructed by transforming the data in the remaining dimension(s) perpendicular to table motion. This approach provides reconstructed images absent of slab-boundary artifacts over a large FOV and reduces the time needed for image acquisition.
In accordance with one aspect of the invention, a method of imaging large volumes without resulting slab-boundary artifacts includes defining a desired FOV larger than an optimal imaging volume of an MR scanner and selecting a slab thickness in a first direction that is smaller than the desired FOV but that is within the optimal imaging volume of the MR scanner. MR data is then acquired by exciting and encoding spins to acquire data that is restricted to the selected slab thickness. The imaging object is then moved continuously with respect to the imaging area, or vice versa. This process is repeated until the desired FOV is fully encoded using a series of cyclically repeated magnetic field gradient waveforms. The selected trajectory need only be uniform in the direction of motion.
In accordance with another aspect of the invention, an MRI apparatus is disclosed to acquire multiple sets of MR data with a moving table and reconstruct MR images without slab-boundary artifacts that includes a magnetic resonance imaging system having an RF transceiver system and a plurality of gradient coils positioned about a bore of a magnet to impress a polarizing magnetic field. An RF switch is controlled by a pulse module to transmit and receive RF signals to and from an RF coil assembly to acquire MR data A patient table is moveable fore and aft in the MR system about the magnetic bore to translate the patient so that an FOV larger than the optimal scanning area of the MRI system can be scanned. A computer is programmed to receive input defining a desired FOV larger than the optimal imaging volume of the MR system. The computer is also capable of defining a fixed slab with respect to the magnet to acquire the MR data therein. The computer then acquires full MR data using frequency encoding in the direction of table motion for a subset of the data in the direction(s) perpendicular to table motion while repositioning the patient table and maintaining the position of the fixed slab. The algorithm is repeated, collecting the necessary MR data across the defined FOV. The patient table can be tracked using the acquired MR data. To reconstruct the image, the MR data is first transformed in the direction of table motion, and then interpolated and aligned to match anatomic locations across slab boundaries. Thereafter, the MR data is transformed with respect to the remaining dimension(s) to reconstruct an MR image.
Yet another aspect of the invention includes a computer program having a set of instructions executed by a computer to control a medical imaging scanner and create large FOV images without significant boundary artifacts. The computer is caused to select an FOV spanning an area greater than a predefined optimal image area of the medical imaging scanner, apply an RF pulse to excite a region of interest, and apply magnetic field gradients to encode the region of interest in a first direction. The volume of interest is limited in the direction of table motion by acquiring the data in such a way that an acquisition filter can be used to restrict the spatial extent. Image data can then be acquired in the first direction as a subset of remaining direction(s) for each acquisition. The computer causes a patient table to move continually with respect to the medical imaging scanner, and repeat the image data acquisition with encoding in the first direction for the encoding in remaining direction(s) until sufficient data is acquired across the entire FOV. An image can be reconstructed without slab-boundary artifacts after aligning anatomical data in the first direction.
Various other features, objects and advantages of the present invention will be made apparent from the following detailed description and the drawings.
The drawings illustrate embodiments presently contemplated for carrying out the invention.
In the drawings:
Referring to
The system control 32 includes a set of modules connected together by a backplane 32a. These include a CPU module 36 and a pulse generator module 38 which connects to the operator console 12 through a serial link 40. It is through link 40 that the system control 32 receives commands from the operator to indicate the scan sequence that is to be performed. The pulse generator module 38 operates the system components to carry out the desired scan sequence and produces data which indicates the timing, strength and shape of the RF pulses produced, and the timing and length of the data acquisition window. The pulse generator module 38 connects to a set of gradient amplifiers 42, to indicate the timing and shape of the gradient pulses that are produced during the scan. The pulse generator module 38 can also receive patient data from a physiological acquisition controller 44 that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes attached to the patient. And finally, the pulse generator module 38 connects to a scan room interface circuit 46 which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 46 that a patient positioning system 48 receives commands to move the patient to the desired position for the scan.
The gradient waveforms produced by the pulse generator module 38 are applied to the gradient amplifier system 42 having Gx, Gy, and Gz amplifiers. Each gradient amplifier excites a corresponding physical gradient coil in a gradient coil assembly generally designated 50 to produce the magnetic field gradients used for spatially encoding acquired signals. The gradient coil assembly 50 forms part of a magnet assembly 52 which includes a polarizing magnet 54 and a whole-body RF coil 56. A transceiver module 58 in the system control 32 produces pulses which are amplified by an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive switch 62. The resulting signals emitted by the excited nuclei in the patient may be sensed by the same RF coil 56 and coupled through the transmit/receive switch 62 to a preamplifier 64. The amplified MR signals are demodulated, filtered, and digitized in the receiver section of the transceiver 58. The transmit/receive switch 62 is controlled by a signal from the pulse generator module 38 to electrically connect the RF amplifier 60 to the coil 56 during the transmit mode and to connect the preamplifier 64 to the coil 56 during the receive mode. The transmit/receive switch 62 can also enable a separate RF coil (for example, a surface coil) to be used in either the transmit or receive mode.
The MR signals picked up by the RF coil 56 are digitized by the transceiver module 58 and transferred to a memory module 66 in the system control 32. A scan is complete when an array of raw k-space data has been acquired in the memory module 66. This raw k-space data is processed as necessary and rearranged into separate k-space data arrays for each image to be reconstructed, and each of these is input to an array processor 68 which operates to Fourier transform the data into an array of image data. This image data is conveyed through the serial link 34 to the computer system 20 where it is stored in memory, such as disk storage 28. In response to commands received from the operator console 12, this image data may be archived in long term storage, such as on the tape drive 30 or any mass storage device, or it may be further processed by the image processor 22 and conveyed to the operator console 12 and presented on the display 16.
The present invention includes a method and system suitable for use with the above-referenced MRI system, or any similar or equivalent system for obtaining MR images. When used with a trajectory that is uniform in the direction of table motion, such as z, data acquired in sections can be reconstructed to form large FOV images without significant slab-boundary artifacts.
As is known, in any MRI system, there is a limited spatial volume having optimal gradient linearity, uniform magnetic polarizing field B0 and uniform RF homogeneity. Larger volumes are traditionally imaged in sections with table motion between scans that result in concatenated images that exhibit discontinuities at slab junctions.
One embodiment to generate large FOV images without slab-boundary artifacts includes a 3D imaging technique that steps the table position in increments during the scanning procedure. In general, after defining the desired FOV, a slab thickness is selected to remain within the MRI system's optimal imaging volume. As a result, the slab thickness is also smaller than the desired spatial coverage. MR data is acquired by repeatedly applying an excitation that excites spins and applying magnetic field gradient waveforms to encode the volume of interest. Preferably, a z slab selective RF pulse is used to excite the volume of interest in the z-direction and a 3D k-space trajectory encodes the volume selected. Some example of 3D k-space trajectories include a 3D EPI k-space trajectory, a cylindrical-stack of EPI k-space trajectory, a stack-of-spirals k-space trajectory, a stack-of-TWIRL k-space trajectory, a stack-of-projection-reconstruction k-space trajectory, and a 3DFT k-space trajectory. More generally though, the spatial extent of the acquired data is restricted in the direction of motion using a method such as using a slab-selective RF pulse or by using an acquisition strategy that allows the acquisition filter to restrict the volume of interest. Therefore, in the examples shown and in a preferred embodiment, the direction of table motion is in the z-direction. However, as one skilled in the art will readily recognize, the x-y-z orientation is arbitrary and the invention is not so limited. In addition, in the example shown and in a preferred embodiment the table moves relative to the system. However, as one skilled in the art will readily recognize it is equivalent to have the system move relative to the subject and the invention is not so limited.
The table is moved to cover the desired FOV while the slab position remains fixed relative to the magnet in the MRI system. At each table position, full k-space data in the direction of table motion is acquired for a subset of the data in the other two dimension(s) by using one or a series of magnetic field gradient waveforms. In this embodiment, at each table position, full k, data is acquired for a subset of kx-ky by using one or a series of magnetic field gradient waveforms. Preferably, RF excitations and magnetic field gradient waveforms are continued during table motion to preserve a steady-state condition.
Referring to
It is noted that when Nz=Nwf, each of the table increments throughout the desired FOV 110 are of equal distance. Accordingly, any of the aforementioned parameters can be adjusted as desired. That is, the slab thickness 114 may be made larger or smaller, or the number of table positions and data acquisition sets 106 can be increased or decreased, above or below the 16 that are shown. The minimum number of table positions desired to reconstruct an image is equal to the number of waveform subsets, Nwf. Additionally, the number of z-pixels, Nz, retained in each excited slab, as well as the number of kx-ky subsets can be modified as desired. Preferably however, Nz is kept greater than or equal to Nwf. One such modification will be described with reference to FIG. 3.
Still referring to
After Fourier transforming with respect to z, or the direction of table motion, data acquired at different table positions are then sorted, interpolated as required, and aligned to match anatomic z locations, thereby filling the z-kx-ky space. In general, the table step distances are multiples of the z-resolution and are selected to ensure complete sampling of the z-kx-ky matrix. As previously indicated, after the acquired data are Fourier transformed with respect to z, data at the slab edges can be dropped as required by imperfections in the slab-select profile. A final reconstructed image 130 is formed by gridding, if necessary, and Fourier transforming the fully sampled data array 120 with respect to x and y.
For complete sampling, the number of z-pixels retained, Nwf must be at least equal to the number of kx-ky waveform subsets, Nwf. However, referring to
Still referring to
To generalize, the table step distances are multiples of the pixel size in the direction of table motion and are selected to ensure complete sampling of the 3D matrix. The number of table steps required depends on the relative number of pixels retained in the excited slab, Nz, for the case of motion in the z direction, and the number of magnetic field gradient waveform subsets required to fully encode the dimension(s) perpendicular to table motion, Nwf. As previously indicated, for complete sampling, Nz, should be at least equal to Nwf. However, if Nz is greater than Nwf, faster overall scan times are achievable. It is noted that the subsets of magnetic field gradient waveforms, (gx gy) are defined by one or a series of such waveforms that differ between subsets. This set, or series, of magnetic field gradient waveforms, that encode the kx, ky subsets, are then repeated in a cyclic manner to obtain the data sets 106,
The present technique can be used with any 3D k-space trajectory that is uniform in the direction of table motion. In one embodiment implemented on a system with table motion in the z direction, a cylindrical stack of EPI trajectory (CSEPI) sequence was utilized that phase encodes in k, between planes of interleaved EPI waveforms that are clipped to a circle in kx-ky. The data is gridded before Fourier transformation in x and y. The waveforms were designed with Gmax=22 mT/m and a slew rate of 120 mT/m/ms. Since k-space data for a single anatomic location is being compiled from data collected at different table positions, any variation in phase or amplitude in the direction of table motion will result in image artifacts. To minimize such artifacts, one embodiment uses a linear phase slab-select RF pulse with sharp transitions and low in-slice ripple, for example 0.5%, and drops one slice at each slab edge. Frequency variations in the direction of table motion that result from magnetic field inhomogeneity are demodulated during reconstruction. Further, decreasing the slab thickness also decreases the sensitivity to such variations.
(Nz−Nwf+1)×(z-resolution), Eqn. 1.
and data is again acquired for the new slab 144, 146, until all the slabs have been scanned across the FOV 148, 168 which completes the scan at 170. Again, the direction of table motion is arbitrarily chosen and represented as a z-direction. However, the invention is not so limited. As one will readily recognize, the assignment of any character set can be used to represent a three dimensional axis.
It is noted that at each table position such as the acquisition at table position 1 in
Referring now to
In another embodiment, the present invention includes a process and system suitable for use with the MRI system of
Moreover, the present invention allows for acquisition of slightly over half of the full k-space data in the transverse dimension(s) with reconstruction using homodyne processing. Further, slightly over half of the full k-space data may be acquired in the readout direction, i.e., frequency encoded direction, which is in the direction of table motion. In addition, phase encodes in the transverse k-space dimension(s) may be acquired in any temporal order and/or with different temporal frequencies. For example, phase encodes closest to the k-space origin could be acquired with a higher temporal frequency than the, phase encodes further from the k-space origin. Furthermore, a subset of phase encodes may be acquired during movement of the table over the FOV during forward movement of the patient table and then acquiring a different subset of phase encodes in the same FOV during reverse table direction. For example, low spatial frequency k-space lines may be acquired during a first pass over the FOV to acquire one time point of time-varying contrast and then acquire remaining phase encoding lines while the patient table is continuously moved back over the FOV.
Additionally, the present invention allows for acquisition of MR data for a larger slab thickness in the direction of continuous table motion than is used for image reconstruction. The data corresponding to the extra slab thickness may be selected so that data on subsequent repetitions of the same phase encoding would have overlapping coverage in the direction of table motion. This overlapping data may be used to estimate table velocity and position and/or other correction factors. Further, this overlapping data may be acquired using the same polarity of readout gradient or an opposite polarity. It may also have the same or different echo time. Data for the extra slab thickness may be acquired for all phase encodes or only the phase encode at the origin of the k-space in the transverse dimension(s). The acquired MR data may also be used for creating an image as well as for estimating correction factors, or may be used as an extra acquisition that is used to estimate the correction factors. For example, MR data corresponding to the kx=ky=0 line may be acquired with double the imaging slab thickness for each repetition of the full set of waveforms. From that, change in table position may be determined by performing a cross-correlation.
The present invention also enables the use of less data than acquired in the direction of table motion thereby allowing the dropping of data at a leading and a trailing edge of acquisition to combat imperfections in filtering. Estimations of imperfections in the magnetic field and center frequency induced by the magnet alone may be determined before starting table motion. Despite continuous table motion, these imperfections will be fixed relative to the magnet. Imperfections induced by object-dependent susceptibility will be superimposed on these fixed errors. The present invention is applicable with any type of image contrast including proton density, T1 weighted, T2 weighted, phase contrast, and SFFP type sequences. Saturation bands may be applied superior and/or inferior to the imaging slab to suppress inflowing blood from one or both directions or to create a different unique contrast effect. Further, the present invention allows for application of echo-planar type readouts in the direction of table motion to acquire different lines in k-space or, depending upon the application, reacquire the same k-space line. Such data may be used for image formation or it may be used to determine correction factors and/or table positions/velocity. If different k-space lines are acquired, the acquisition order may be linear or interleaved. The present invention also allows for correcting gradient non-linearities in a direction of table motion on each acquired line and for transverse gradient non-linearities by correcting each kx-ky location individually based on the gradient non-linearity at the table location at the time of acquisition.
As is known and previously stated, in any MRI system, there is a limited spatial volume having optimal gradient linearity, uniform magnetic polarizing field B0, and uniform RF homogeneity. Larger volumes are traditionally imaged in sections with table motion between scans that result in concatenated images that exhibit discontinuities at slab junctions. The present invention includes an imaging technique that moves the table position continuously during the scanning procedure to create large FOV images without slab boundary artifacts and without adding scan time for table motion between fixed-location acquisitions. Acquired images can be 2D single-slice, 2D multi-slice, 2D projection, or 3D.
In general, after defining the desired FOV, a slab thickness is selected to remain within the MRI system's optimal imaging volume. As a result, the slab thickness is also smaller than the desired spatial coverage. MR data is acquired by repeatedly applying an excitation that excites spins and applying magnetic field gradient waveforms to encode the volume of interest. Preferably, the spatial extent of the acquired data is restricted in the direction of motion using applied magnet field gradient waveforms that allow the acquisition filter to restrict the slab thickness. The FOV can be restricted using spatially-selective RF pulses in directions perpendicular to the direction of table motion. In the examples shown and in one preferred embodiment, the direction of continuous table motion is in the z-direction. However, as one skilled in the art will readily recognize, the x-y-z orientation is somewhat arbitrary and the invention is therefore not so limited. In addition, in the example shown and in a preferred embodiment the table moves relative to the system. However, as one skilled in the art will readily recognize it is equivalent to have the system move relative to the subject and the invention is accordingly not so limited.
The table is moved to cover the desired FOV while the slab position remains fixed relative to the magnet in the MRI system. During each acquisition, full k-space data representing the selected slab in the direction of table motion is acquired for a subset of the data in the other dimension(s) by using one or a series of magnetic field gradient waveforms.
Referring to
Any of the aforementioned parameters can be adjusted as desired. That is, the slab thickness 314 may be made larger or smaller and the number of data acquisitions 306 can be increased or decreased, above or below the 16 that are shown. The maximum table velocity is determined by the slab thickness divided by the time required to collect all the data corresponding to one cycle through the set of gradient waveforms. Additionally, the number of z-pixels, Nz, retained in each excited slab, as well as the number of gradient waveforms in a set, Nwf, can be modified as desired.
Still referring to
After Fourier transforming with respect to z, or the direction of table motion, data acquired at different table positions are then sorted and aligned to match anatomic z locations, thereby filling the z-ktransverse space. Data can be interpolated as necessary to improve the accuracy of the alignment. As previously indicated, after the acquired data are Fourier transformed with respect to z, data at the slab edges can be dropped as required by imperfections in filter performance or for other reasons as described below. A final reconstructed image 330 is formed by gridding and/or performing other processing as necessary, and Fourier transforming the fully sampled data array 320 in remaining dimension(s).
It is noted that the subsets of magnetic field gradient waveforms, (gx gy, gz) are defined by one or a series of such waveforms that differ between subsets. This set, or series, of magnetic field gradient waveforms, that encode the ktransverse subsets, are then repeated in a cyclic manner to obtain the data sets 306, FIG. 8.
The present technique can be used with any 2D or 3D k-space trajectory that is uniform in the direction of table motion. For example, a 2D or 3D acquisition could be performed with the frequency-encoding axis in the direction of table motion and phase encoding gradients positioned on a Cartesian grid in the k-space dimension(s) that are perpendicular to the direction of table motion. Phase encoding gradients for a 3D acquisition could also be positioned on concentric circles, in the shape of a spiral, in rays from the center of k-space, or in any other pattern. In a 2D or 3D acquisition, phase encoding gradients could be placed in a pattern with density that varies with k-space position. If a non-Cartesian and/or non-uniform placement of ktransverse locations is utilized, the data would be gridded in the transverse dimension(s) during reconstruction.
Partial k-space techniques could be employed. For partial k-space in the transverse direction(s), slightly more than half of the data is acquired in the transverse direction. Homodyne processing would then be utilized in the reconstruction. Partial k-space techniques could also be used in the direction of table motion. Homodyne processing could also be used in this case.
Echo-planar type acquisitions could be employed with readout in the direction of table motion. This would allow full k-space data in the direction of table motion to be acquired for multiple transverse k-space positions following each RF excitation. Alternatively, full k-space data in the direction of table motion could be acquired multiple times for the same transverse k-space position following a single RF excitation. Such data acquisition could allow for the calculation of image correction factors.
Different k-space locations in the transverse dimension(s) may also be acquired with different slab thicknesses in the direction of table motion and therefore could be acquired with different temporal frequency during table motion.
As another variation, k-space locations in the transverse dimension(s) can be acquired in a variety of orders to achieve different modulations of k-space or different contrast properties. K-space locations in the transverse dimension(s) may be acquired sequentially from one edge of k-space to the other, centrically from the center out, in an interleaved fashion, or in any other pattern that is desired.
K-space locations in the transverse dimension(s) could also be divided into an arbitrary number of subsets that could be acquired in multiple table passages. For example, you could acquire the ktransverse data in two stages. First, you could acquire the ktransverse locations within a certain radius of the center of k-space. Since less time would be required to acquire this subset of the required data, the table could move faster over the selected FOV. As an example application, this would allow the center of k-space, which has the largest impact on image contrast, to be acquired at the point of optimal contrast following administration of an extrinsic MR contrast agent. The table would then be moved back over the same FOV again to acquire higher spatial frequency k-space data.
As another extension, data for a larger selected slab thickness could be acquired in the direction of table motion. The extra data could be acquired for all transverse k-space positions or only some positions, for example the kx=ky=0 line. The extra slab thickness could be selected such that the acquired data on subsequent repetitions of the same gradient waveform set would acquire data from overlapping spatial regions. This overlapping data could be used to estimate the distance of table motion between the two acquisitions or could be used to estimate other image correction factors. A portion of the data from these acquisitions with the increased slab thickness could be used to reconstruct the image in the normal fashion described in
Since k-space data for a single anatomic location in the direction of table motion is compiled from data collected at different table positions, any variation in phase or amplitude in the direction of table motion will result in image artifacts. To minimize such artifacts, frequency variations in the direction of table motion that result from magnetic field inhomogeneity are demodulated during reconstruction. Frequency variations induced by the applied external magnetic field may be estimated before table motion begins. Gradient non-linearities can be corrected by assuming that these imperfections are separable between the direction of table motion and the transverse dimension(s). Gradient non-linearities in the direction of table motion are corrected for each data acquisition. Transverse non-linearities are corrected using a different factor for each transverse k-space location, or group of locations, that is/are determined using knowledge of the table position at the time of data acquisition. Signal phase induced by table motion during acquisition may be corrected during reconstruction.
It is noted that at each data acquisition position such as the acquisition at position 1 in
Referring now to
Referring to
If at 412, the process determines that there is sufficient data acquired for the selected subset to reconstruct the desired FOV 426, the process determines if the last subset has been encoded at 428. If not 428, 430, the process increments the subset and resets the waveform, i.e., J=1 at 432. At 434, the table is decelerated to zero velocity and subsequently reversed and accelerated until a desired constant velocity is reached. An RF and Jth gradient waveform of the I+1 subset on each gradient axis are then applied at 410 while the table moves in a second or reverse direction. If at 428 the last subset has been encoded and sufficient data has been acquired to reconstruct a final image for the desired FOV 435, the scan ends at 436.
Methods such as those described here could be used to acquire images with a variety of contrast properties. As examples, images could display proton-density weighting, T1 weighting, T2 weighting, or any other weighting available in stationary-table imaging. Image contrast could be manipulated by interleaving the repetition of imaging RF and gradient waveforms with RF and gradient waveforms designed to saturate, invert, or otherwise manipulate spins in regions superior or inferior to the selected imaging slab or with specific spectral properties. Such saturation regions could suppress in-flowing materials from the direction saturated, could saturate lipid spins, or could achieve different contrast effects due to the movement of the patient table.
There are a number of advantages associated with the present invention. For example, the present invention allows for 3D imaging with continuous motion and frequency encoding in the direction of the patient and/or patient table motion. Two-dimensional single-slice or multi-slice imaging with continuous motion and frequency encoding in the direction of table motion is also possible. Further, 2D projection with continuous motion and frequency encoding in the direction of motion may be achieved with the present invention. The present invention further allows for different spatial distributions of phase encode positions such as spiral, concentric rings, variable density, etc. and using gridding reconstruction techniques. The present invention further allows for partial k-space imaging methods with continuous table motion as well as different phase encode orderings with continuous table motion. The present invention also enables for tracking of table position using acquired larger slab thickness in a direction of table movement thereby allowing for the overlapping of data. Additionally, the table may be moved several times over a desired anatomy to achieve greater image quality. The present invention also allows for performing off-resonance and gradient linearity correction with moving table as well as using saturation bands with a moving table. EPI readout in a direction of table motion to either reacquire the same k-space line or a different line is possible with the present invention. Further, the present invention allows for different temporal update rates for different phase encode positions.
Accordingly, the present invention includes a method of imaging large volumes without resulting slab-boundary artifacts that includes defining a desired FOV larger than an optimal imaging volume of an MR scanner, selecting a slab thickness in a first direction that is smaller than the desired FOV, but within the optimal imaging volume of the MR scanner, then exciting and encoding spins to acquire data that is restricted to the selected slab thickness. MR data is acquired that includes full encoding data in the direction of table motion for a subset of the remaining direction(s) and continual moving one of the optimal imaging volume and an imaging object to acquire another set of MR data. This process of acquiring sets of MR data and continual moving the table is repeated until the desired FOV is covered. Acquired MR data can be used and analyzed to track motion of the table. In a preferred embodiment, each MR data acquisition includes acquiring all k-space data in the direction of table motion for a selected subset of the remaining transverse dimension(s), and further includes defining a set or series of magnetic field gradient waveforms that encode a k-space trajectory that is uniform in a k-space dimension(s) along the direction of movement.
The invention also includes an MRI apparatus, such as that disclosed with reference to
Additionally, the present invention also includes a computer program to control a medical image scanner and create images across scanning boundaries without boundary artifacts. The computer program has a set of instructions to control a computer to select an FOV spanning an area greater than a predefined optimal imaging area of the medical image scanner and acquire MR data by repeatedly applying an RF pulse to excite a region in the selected FOV and magnetic field gradients to encode the region in a direction. The instructions further control the computer to apply a k-space trajectory to encode the region in the direction of table motion and acquire data for a subset of the remaining direction(s). The computer program includes instructions to continuously reposition the imaging object with respect to the predefined optimal imaging area or vice versa and repeat the image data acquisition and the continuous imaging area repositioning until complete image data are acquired across the entire FOV to reconstruct an image of the FOV. Positioning of the optimal imaging area can be tracked using acquired MR data, or other known processes.
The present invention results in images that are free of slab-boundary artifacts, reduces the time required for data acquisition, and allows spatial coverage that is limited only by the range of table motion.
The present invention has been described in terms of the preferred embodiment, and it is recognized that equivalents, alternatives, and modifications aside from those expressly stated, are possible and within the scope of the appending claims.
The present application is a continuation of U.S. patent application Ser. No. 09/682,699 filed Oct. 5, 2001, now U.S. Pat. No. 6,794,869 which is a continuation-in-part of U.S. patent application Ser. No. 09/681,420 filed Mar. 30, 2001 now U.S. Pat. No. 6,801,034.
Number | Name | Date | Kind |
---|---|---|---|
5928148 | Wang et al. | Jul 1999 | A |
6385478 | Hajnal | May 2002 | B1 |
6445181 | Pelc et al. | Sep 2002 | B1 |
6584337 | Dumoulin et al. | Jun 2003 | B2 |
6617850 | Welch et al. | Sep 2003 | B2 |
20020021128 | Kuhara | Feb 2002 | A1 |
20020115929 | Machida | Aug 2002 | A1 |
20020140423 | Brittain | Oct 2002 | A1 |
20020173715 | Kruger et al. | Nov 2002 | A1 |
20030011369 | Brittain et al. | Jan 2003 | A1 |
20030100825 | Demoulin et al. | May 2003 | A1 |
20030102864 | Welch et al. | Jun 2003 | A1 |
20040155654 | Brittain | Aug 2004 | A1 |
Number | Date | Country |
---|---|---|
H6-304153 | Nov 1994 | JP |
H6-311977 | Nov 1994 | JP |
H8-173396 | Jul 1996 | JP |
Number | Date | Country | |
---|---|---|---|
20040155654 A1 | Aug 2004 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 09682699 | Oct 2001 | US |
Child | 10749167 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 09681420 | Mar 2001 | US |
Child | 09682699 | US |