Search

An interactive real-time imaging system capable of rapid coronary artery imaging is described. High-resolution spiral and circular echo planar trajectories were used to achieve 0.8 x 1.6 mm2 resolution in 135 ms (CEPI) or 1.13 x 1.13 mm2 resolution in 189 ms (spirals), over a 20-cm FOV. Using a sliding window reconstruction, display rates of up to 37 images/sec were achieved. Initial results indicate this technique can perform as a high-quality 2D coronary localizer and with SNR improvement may enable rapid screening of the coronary tree., (Copyright 2001 Wiley-Liss, Inc.)

A new spiral imaging technique incorporates the acquisition of a field map into imaging interleaves. Variable density spiral trajectories are designed to oversample the central region of k-space, and interleaves are acquired at two different echo times. A field map is extracted from this data and multifrequency reconstruction is used to form an off-resonance corrected image using the entire dataset. Simulation, phantom, and in vivo results indicate that this technique can be used to achieve higher image and/or field map spatial resolution compared to conventional techniques. Magn Reson Med 45:521-524, 2001., (Copyright 2001 Wiley-Liss, Inc.)

A connectivity algorithm combined with a new gray-level morphological filter dramatically improves the segmentation of tortuous coronary arteries from 3D MRI. Small coronary arteries are segmented from the larger ventricles with a new filter. These blood vessels are segmented from the noise background with connectivity. Coronary angiograms were computed in nine datasets acquired on volunteers with 3D stack of spirals and contrast-enhanced navigator sequences by both a maximum intensity projection and surface rendering. Surface images provided depth information needed to distinguish branching arteries from crossing veins. Magn Reson Med 43:892-895, 2000., (Copyright 2000 Wiley-Liss, Inc.)

Density-weighted k-space sampling with spiral trajectories is used to reduce spatial side lobes in chemical-shift imaging (CSI). In this method, more time is spent collecting data at the center of k space and less time at the edges of k space in order to make the sampling density proportional to a given apodization function, subject to constraints imposed by gradient performance and Nyquist sampling. The efficient k-space coverage of spiral-based trajectories enables good control over the sampling density within practical in vivo scan times. The density-weighted acquisition is compared to a conventional, nonweighted spiral sampling without the application of a window function. For a fixed voxel size and imaging time, the noise variance is observed to be the same for both cases, while spatial side lobes are greatly reduced with the variable-density sampling. This method is demonstrated on a normal volunteer by imaging of brain metabolites at 1.5 T with both single slice CSI and volumetric CSI. Magn Reson Med 42:314-323, 1999., (Copyright 1999 Wiley-Liss, Inc.)

A cardiac motion compensation method using magnetic resonance signal-based triggering is presented. The method interlaces a triggering pulse sequence with an imaging sequence. The triggering sequence is designed to measure aortic blood velocity, from which cardiac phase can be inferred. The triggering sequence is executed repeatedly and the acquired data processed after each sequence iteration. When the desired phase of the cardiac cycle is detected, data are acquired using the imaging sequence. A signal-processing unit of a conventional scanner is used to process the triggering data in real time and issue triggering commands. Alternatively, a workstation, with a bus adaptor, can access data as they are acquired, process and display the data, and issue triggering commands. With a graphical user interface, the triggering pulse sequence and data-processing techniques can be modified instantaneously to optimize triggering. The technique is demonstrated with coronary artery imaging using both conventional two-dimensional Fourier transform scans and spiral trajectories.

In this work, three-dimensional (3D) spiral imaging has been utilized for magnetic resonance coronary angiography. Spiral-based 3D techniques can dramatically reduce imaging time requirements compared with 3D Fourier Transform imaging. The method developed here utilized a "stack of spirals" trajectory, to traverse 3D k-space rapidly. Both thick-slab volumes encompassing the entire coronary tree with isotropic resolution and thin-slab volumes targeted to a particular vessel of interest were acquired. Respiratory compensation was achieved using the diminishing variance algorithm. T2-prepared contrast was also applied in some cases to improve contrast between vessel and myocardium, while off-resonance blurring was minimized by applying a linear correction to the acquired data. Images from healthy volunteers were displayed using a curved reformatting technique to view long segments of vessel in a single projection. The results demonstrate that this 3D spiral technique is capable of producing high-quality coronary magnetic resonance angiograms.

Coronary arteries are segmented from the blood pool using mathematical morphology operations from a 3D magnetic resonance spiral acquisition on a continuously breathing healthy volunteer. The segmented volume is maximal intensity projected at different views to yield coronary angiograms showing the left anterior descending artery (LAD), right coronary artery (RCA), and left circumflex artery (LCX). Magnetic resonance coronary angiography provides a retrospective rotating view of the coronary artery tree that complements oblique reformatted sections.

A real-time interactive MRI system capable of localizing coronary arteries and imaging arrhythmic hearts in real-time is described. Non-2DFT acquisition strategies such as spiral-interleaf, spiral-ring, and circular echo-planar imaging provide short scan times on a conventional scanner. Real-time gridding reconstruction at 8-20 images/s is achieved by distributing the reconstruction on general-purpose UNIX workstations. An X-windows application provides interactive control. A six-interleaf spiral sequence is used for cardiac imaging and can acquire six images/s. A sliding window reconstruction achieves display rates of 16-20 images/s. This allows cardiac images to be acquired in real-time, with minimal motion and flow artifacts, and without breath holding or cardiac gating. Abdominal images are acquired at over 2.5 images/s with spiral-ring or circular echo-planar sequences. Reconstruction rates are 8-10 images/s. Rapid localization in the abdomen is demonstrated with the spiral-ring acquisition, whereas peristaltic motion in the small bowel is well visualized using the circular echo-planar sequence.

Echo-planar imaging (EPI) is very susceptible to flow artifacts. Two ways to improve its flow properties are presented. First, "partial flyback" is proposed to reduce artifacts arising from flow in the readout direction. Near the center of k-space, only the even echoes of the EPI echo-train are used. Partial flyback is shown to improve the readout-flow properties at the expense of a slight worsening of the phase-encode flow and off-resonance properties. We recommend that the flyback region acquire 95% of the energy in k-space. Second, "inside-out" EPI is used to reduce artifacts arising from flow in the phase-encode direction. Data collection begins at the center of k-space, with separate interleaves to acquire the top and bottom, halves of k-space. Partial flyback is combined with partial-Fourier EPI and inside-out EPI. Partial-flyback inside-out EPI has worse off-resonance properties than partial-flyback partial-Fourier EPI but demonstrates better flow properties and does not require partial k-space reconstruction.

A fast and robust method for correcting magnetic resonance image distortion due to field inhomogeneity is proposed and applied to spiral k-space scanning. The method consists of acquiring a local field map, finding the best fit to a linear map, and using it to deblur the image distortions due to local frequency variations. The linear field map is determined using a maximum likelihood estimator with weights proportional to the pixel intensity. The method requires little additional computation and is robust in low signal regions and near abrupt field changes. Additionally, it can be used in combination with other deblurring methods. The application of this method is illustrated in conjunction with a multislice, T2-weighted, breath-held spiral scan of the liver.

A technique has been developed whereby motion can be detected in real time during the acquisition of data. This enables the implementation of several algorithms to reduce or eliminate motion effects from an image as it is being acquired. One such algorithm previously described is the acceptance/rejection method. This paper deals with another real-time algorithm called the diminishing variance algorithm (DVA). With this method, a complete set of preliminary data is acquired along with information about the relative motion position of each frame of data. After all the preliminary data are acquired, the position information is used to determine which data frames are most corrupted by motion. Frames of data are then reacquired, starting with the most corrupted one. The position information is continually updated in an iterative process; therefore, each subsequent reacquisition is always done on the worst frame of data. The algorithm has been implemented on several different types of sequences. Preliminary in vivo studies indicate that motion artifacts are dramatically reduced.

The imaging of dynamic processes in the body is of considerable interest in interventional examinations as well as kinematic studies, and spiral imaging is a fast magnetic resonance imaging technique ideally suited for such fluoroscopic applications. In this manuscript, magnetic resonance fluoroscopy pulse sequences in which interleaved spirals are used to continuously acquire data and reconstruct one movie frame for each repetition time interval are implemented. For many applications, not all of k-space needs to be updated each frame, and nonuniform k-space sampling can be used to exploit this rapid imaging strategy by allowing variable update rates for different spatial frequencies. Using the appropriate reconstruction algorithm, the temporal updating rate for each spatial frequency is effectively proportional to the corresponding k-space sampling density. Results from a motion phantom as well as in in vivo gadolinium diethylenetriaminopentaacetic acid (Gd-DTPA) bolus tracking studies in a rat model demonstrate the high temporal resolution achievable using these techniques as well as the tradeoffs available with nonuniform sampling densities. This paper focuses on the acquisition of real-time dynamic information, and all images presented are reconstructed retrospectively. The issues of real-time data reconstruction and display are not addressed.

A technique is described for discriminating blood-oxygen-level-dependent (BOLD) signal changes originating from large venous vessels and those that arise from the cortical parenchyma based on examining the temporal delay of each pixel's response. Photic stimulation experiments were performed with a conventional 1.5 T scanner and correlated each pixel's time-course with sine and cosine functions at the frequency of the stimulus. It was found that the signal in pixels anatomically associated with gray matter was delayed between 4 and 8 s compared with the stimulus, whereas the signal in pixels correlated with visible vessels and sulci was generally delayed from 8 to 14 s. This larger delay is consistent with the longer time required for blood to reach the larger vessels. Finally, stimulus-induced NMR phase changes were observed for the largest vessels, which are attributed to bulk susceptibility shifts.

A magnetization-prepared, T2-weighted sequence (T2 Prep) is used to suppress muscle and venous structures. When combined with lipid suppression, this technique improves the visualization of the coronary arteries. T2 Prep was designed to be rebust in the presence of flow as well as B0 and B1 inhomogeneities and may be combined with virtually any imaging technique. Here, it is implemented with both a single-slice spiral acquisition and a multi-slice spiral method that acquires up to 15 slices in a single breath-holding interval.

A velocity k-space formalism facilitates the analysis of flow effects for imaging sequences involving time-varying gradients such as echo-planar and spiral. For each sequence, the velocity k-space trajectory can be represented by kv(kr); that is, its velocity-frequency (kv) position as a function of spatial-frequency (kr) position. In an echo-planar sequence, kv is discontinuous and asymmetric. However, in a spiral sequence, kv is smoothly varying, circularly symmetric, and small near the kr origin. To compare the effects of these trajectory differences, simulated images were generated by computing the k-space values for an in-plane vessel with parabolic flow. Whereas the resulting echo-planar images demonstrate distortions and ghosting that depend on the vessel orientation, the spiral images exhibit minimal artifacts.

A technique has been developed whereby motion can be detected in real time during the acquisition of data. This enables the implementation of an algorithm to accept or reject and reacquire data during a scan. Frames of data with motion are rejected and reacquired on the fly so that by the end of the scan, a complete motion-free data set has been acquired. The algorithm has been implemented on several different types of sequences. Preliminary in vivo studies indicate that motion artifacts are dramatically reduced.

Time-resolved velocity imaging using the magnetic resonance phase contrast technique can provide clinically important quantitative flow measurements in vivo but suffers from long scan times when based on conventional spin-warp sequences. This can be particularly problematic when imaging regions of the abdomen and thorax because of respiratory motion. We present a rapid phase contrast sequence based on an interleaved spiral k-space data acquisition that permits time-resolved, three-direction velocity imaging within a breath-hold. Results of steady and pulsatile flow phantom experiments are presented, which indicate excellent agreement between our technique and through plane flow measurements made with an in-line ultrasound probe. Also shown are results of normal volunteer studies of the carotids, renal arteries, and heart.

A spectroscopic imaging sequence incorporating a two-shot lactate editing method was used in two human brain studies to image lactate and NAA. The subtractive editing method allows separate images of lactate, NAA, and lipids to be collected during a single study with no SNR penalty. The sequence uses a spectral-spatial excitation for slice selection and water suppression, and employs inversion recovery and an echo time of 136 ms for additional lipid suppression.

A flow-independent method for imaging the coronary arteries within a breath-hold on a standard whole-body MR imager was developed. The technique is based on interleaved spiral k-space scanning and forms a cardiac-gated image in 20 heartbeats. The spiral readouts have good flow properties and generate minimal flow artifacts. The oblique slices are positioned so that the arteries are in the plane and so that the chamber blood does not obscure the arteries. Fat suppression by a spectral-spatial pulse improves the visualization of the arteries.

For several non-2D Fourier transform imaging methods, off-resonant reconstruction does not just cause geometric distortion, but changes the shape of the point spread function and causes blurring. This effect is well known for projection reconstruction and spiral k-space scanning sequences. We introduce here a method that automatically removes blur introduced by magnetic field inhomogeneity and susceptibility without using a resonant frequency map, making these imaging methods more useful. In this method, the raw data are modulated to several different frequencies and reconstructed to create a series of base images. Determination of degree of blur is done by calculating a focusing measure for each point in each base image and a composite image is then constructed using only the unblurred regions from each base image. This method has been successfully applied to phantom and in vivo images using projection-reconstruction and spiral-scan sequences.

Using a k-space interpretation of small-tip excitation, a single excitation pulse has been designed that is simultaneously selective in space and resonant frequency. An analytic expression for the response of this pulse has been derived. The pulse has been implemented on a 1.5-T imaging system. The pulse has been applied to a rapid gradient-echo imaging sequence that forms both water and fat images within a breath-holding interval. These rapid images are free of the chemical shift artifacts at organ boundaries that typically afflict conventional rapid images. The pulse can be applied to a variety of other sequences, such as multislice water/fat sequences and rapid k-space scanning sequences.