Whole-body MR angiography using variable density sampling and dual-injection bolus-chase acquisition

Conventional bolus-chase acquisition generates peripheral runoff images using a single injection of the contrast material. Low spatial resolution, small slice coverage and venous contamination are major problems especially in the distal stations. A technique is presented herein in which whole-body magnetic resonance angiography is performed using a dual-contrast-injection four-station acquisition protocol. Bolus sharing was performed between two stations: the abdomen and calf stations share the first bolus injection, while the thorax and thigh stations share the second bolus injection. The combination of variable density sampling and elliptical centric acquisition order was applied to the abdomen and thorax stations. The scan time was extended to generate high spatial resolution arterial phase images with broad slice coverage for the calf and thigh stations. The feasibility of this technique was demonstrated using phantom and in vivo human volunteer studies.     

Quiet imaging with interleaved spiral read-out

The acoustic noise generated during an MRI sequence can be effectively reduced with the help of soft gradient pulses using sinusoidal ramps. The long slope duration, however, leads to long acquisition times. The use of interleaved spiral trajectories, calculated with long gradient slopes, is proposed to reduce the acquisition time while maintaining low acoustic noise levels. The practibility of this approach is demonstrated on phantom and volunteer images.     

Fast 3D coronary artery contrast-enhanced magnetic resonance angiography with magnetization transfer contrast, fat suppression and parallel imaging as applied on an anthropomorphic moving heart phantom

A magnetic resonance sequence for high-resolution imaging of coronary arteries in a very short acquisition time is presented. The technique is based on fast low-angle shot and uses fat saturation and magnetization transfer contrast prepulses to improve image contrast. GeneRalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) is implemented to shorten acquisition time. The sequence was tested on a moving anthropomorphic silicone heart phantom where the coronary arteries were filled with a gadolinium contrast agent solution, and imaging was performed at varying heart rates using GRAPPA. The clinical relevance of the phantom was validated by comparing the myocardial relaxation times of the phantom's homogeneous silicone cardiac wall to those of humans. Signal-to-noise ratio and contrast-to-noise ratio were higher when parallel imaging was used, possibly benefiting from the acquisition of one partition per heartbeat. Another advantage of parallel imaging for visualizing the coronary arteries is that the entire heart can be imaged within a few breath-holds.     

Artifact-free black-blood cine cardiac imaging in a single breath-hold

Fast imaging using the STimulated Echo Acquisition Mode (STEAM) sequence can produce cine images of the heart with black-blood contrast. Nevertheless, correction of deformation-related artifacts is required in order to maintain myocardial signal throughout the cardiac cycle. Recent work by our group has eliminated this artifact by combining two STEAM sequences acquired with two different demodulation gradients. Unfortunately, these two STEAM sequences were acquired on two separate breath-holds; thus, scan time doubled. In this work, we present a technique to reduce the total scan time by one half, without sacrificing image quality. The technique is based on interleaving two demodulations within one acquisition in order to obtain quality cine images of the heart in a single breath-hold. The technique was tested on animal models and human subjects, and the impact of interleaved acquisition on image quality was studied using quantitative and qualitative measures.     

Simultaneous acquisition of phosphocreatine and inorganic phosphate images for Pi:PCr ratio mapping using a RARE sequence with chemically selective interleaving

The ratio of inorganic phosphate to phosphocreatine (Pi:PCr) is a validated marker of mitochondrial function in human muscle. The magnetic resonance imaging rapid acquisition with relaxation enhancement (RARE) pulse sequence can acquire phosphorus-31 (³¹P) images with higher spatial and temporal resolution than traditional spectroscopic methods, which can then be used to create Pi:PCr ratio maps of muscle regions. While the ³¹P RARE method produces images that reflect the content of the ³¹P metabolites, it has been limited to producing an image of only one chemical shift in a scan. This increases the scan time required to acquire images of multiple chemical shifts as well as the likelihood of generating inaccurate Pi:PCr maps due to gross motion. This work is a preliminary study to demonstrate the feasibility of acquiring Pi and PCr images in a single scan by interleaving Pi and PCr chemical shift acquisitions using a chemically selective radiofrequency excitation pulse. The chemical selectivity of the excitation pulse evaluated and the Pi:PCr maps generated using the interleaved Pi and PCr acquisition method with the subject at rest and during exercise are compared to those generated using separate Pi and PCr acquisition scans. A paired t test indicated that the resulting Pi:PCr ratios for the exercised forearm muscle regions were not significantly different between the separate Pi and PCr acquisition method (3.18±1.53) (mean±standard deviation) and the interleaved acquisition method (3.41±1.66). This work demonstrates the feasibility of creating Pi:PCr ratio maps in human muscle with Pi and PCr images acquired simultaneously by interleaving between the Pi and PCr resonances in a single scan.     

Self-gated Fourier velocity encoding

Self-gating is investigated to improve the velocity resolution of real-time Fourier velocity encoding measurements in the absence of a reliable electrocardiogram waveform (e.g., fetal magnetic resonance or severe arrhythmia). Real-time flow data are acquired using interleaved k-space trajectories which share a common path near the origin of k-space. These common data provide a rapid self-gating signal that can be used to combine the interleaved data. The combined interleaves cover a greater area of k-space than a single real-time acquisition, thereby providing higher velocity resolution for a given aliasing velocity and temporal resolution. For example, this approach provided velocity spectra with a temporal resolution of 19 ms and velocity resolution of 22 cm/s over an 818 cm/s field-of-view. The method was validated experimentally using a computer-controlled pulsatile flow apparatus and applied in vivo to measure aortic-valve flow in a healthy volunteer.     

Magnetic resonance imaging of the abdominal aorta and iliac vessels using combined 3-D gadolinium-enhanced MRA and gadolinium-enhanced fat-suppressed spoiled gradient echo sequences.

This study evaluates a combined protocol consisting of breath hold immediate post gadolinium 3-D gradient echo MR angiography and blood pool phase gadolinium-enhanced breath hold 2-D fat-suppressed spoiled gradient echo (SGE) sequences in the examination of diseases of the abdominal aorta and iliac vessels. Thirty-two patients with suspected disease of the abdominal aorta, major aortic branches, or iliac vessels underwent MR angiographic study from January 1996 to January 1997. Examinations were performed on a 1.5 T MR imager using 2-D axial SGE, coronal 3-D fast imaging in steady state precession (3-D FISP) following bolus administration of 40 mL of gadolinium, and axial and coronal blood pool phase gadolinium-enhanced fat-suppressed SGE. Post-processed data, including 3-D reconstructions using maximum intensity projection (MIP), targeted MIP, and multiplanar reconstruction (MPR) were evaluated. MR findings in all patients were correlated as follows: surgery (13 patients), angiography (11 patients), contrast enhanced CT (3 patients), non-contrast enhanced CT (1 patient), color doppler US (2 patients), and previous MR study (2 patients). MR findings correlated closely with findings at surgery or other imaging studies in 31 of 32 patients. One patient had renal artery occlusion that was misinterpreted as mild stenosis. The following vascular diseases were present: aneurysm disease [10 patients: aortic aneurysm (8 patients), inflammatory aneurysm (2 patients)], thoracoabdominal aortic dissection (2 patients), arteriovenous fistula (1 patient), stenoses and/or occlusion of the abdominal aorta, major aortic branches and iliac vessels [12 patients: stenoses and/or occlusion of the abdominal aorta with stenoses of the iliac vessels (9 patients), renal artery stenosis (2 patients), occlusion of the abdominal aorta (1 patient)], and occluded artery to pancreatic transplant artery (1 patient). Five patients had normal studies. The 3-D FISP technique accurately defined the luminal contours of vessels, allowing precise depiction of vessel stenosis (i.e., renal artery stenosis or common iliac artery stenosis) and clear demonstration of relationship of aortic branch vessels (i.e., renal arteries) to underlying aortic pathology (i.e., aortic aneurysm or dissection). Blood pool phase gadolinium-enhanced fat-suppressed SGE images were useful in the evaluation of the external surface of vessel walls, and providing accurate measurement of aneurysm diameter and other associated vascular entities (i.e., inflammatory aneurysm, left-sided IVC). Targeted MIP or MPR reconstruction were important for assessing stenoses of medium sized vessels such as renal arteries and branches of the iliac arteries, and for identifying accessory arteries. The combination of immediate post gadolinium 3-D FISP and blood pool phase gadolinium-enhanced fat-suppressed SGE is useful in the evaluation of the abdominal aorta, major aortic branches and iliac vessels. Immediate post gadolinium 3-D FISP images provides diagnostically useful information regarding vessel luminal contour, while blood pool phase gadolinium-enhanced fat-suppressed SGE provides ancillary information on the vessel wall and surrounding tissue.     

Effect of increased repetition time TR on precision of inversion-recovery T(1) measurements

The effect of increased repetition time, TR, on the precision of inversion-recovery measurements of the spin-lattice relaxation time, T(1), was calculated theoretically, simulated numerically, and measured experimentally. All three methods yielded similar results. With constant inversion times, the T(1) precision was independent of TR. Therefore, 1) multiple-slice inversion-recovery fast-spin-echo T(1) maps should be made one slice at a time, not with interleaved acquisition, and 2) once the longest inversion time t(i) has been set, TR should be set just enough longer than the longest t(i) to allow data acquisition.     

Motion correction properties of the shells k-space trajectory

The feasibility of a k-space trajectory that samples data on a set of 3D shells is demonstrated with phantom and volunteer experiments. Details of an interleaved multi-shot, helical spiral pulse sequence and a gridding reconstruction algorithm that uses Voronoi diagrams are provided. The motion-correction properties of the shells k-space trajectory are described. It is shown that when used in conjunction with three point markers, k-space data acquired with the shells trajectory provide a generalization of the RINGLET method, allowing for correction of arbitrary rigid-body motion with six degrees of freedom. Use of dedicated navigator echoes or redundant acquisitions of k-space data are not required. Retrospective motion correction is demonstrated with controlled phantom experiments and with seven healthy human volunteers. The motion correction is shown to improve the images, both qualitatively and quantitatively with a metric calculated from image entropy. Advantages and challenges of the shells trajectory are discussed, with particular attention to acquisition efficiency.     

Reduction of flow- and eddy-currents-induced image artifacts in coronary magnetic resonance angiography using a linear centric-encoding SSFP sequence

Coronary magnetic resonance angiography (MRA) acquired using steady-state free precession (SSFP) sequences tends to suffer from image artifacts caused by local magnetic field inhomogeneities. Flow- and gradient-switching-induced eddy currents are important sources of such phase errors, especially under off-resonant conditions. In this study, we propose to reduce these image artifacts by using a linear centric-encoding (LCE) scheme in the phase-encoding (PE) direction. Abrupt change in gradients, including magnitude and polarity between consecutive radiofrequency cycles, is minimized using the LCE scheme. Results from numeric simulations and phantom studies demonstrated that signal oscillation can be markedly reduced using LCE as compared to conventional alternating centric-encoding (ACE) scheme. The image quality of coronary arteries was improved at both 1.5 and 3.0 T using LCE compared to those acquired using ACE PE scheme (1.5 T: ACE/LCE=2.2+/-0.8/3.0+/-0.6, P=.02; 3.0 T: ACE/LCE=2.1+/-1.1/3.0+/-0.8, P=.01). In conclusion, flow- and eddy-currents-induced imaging artifacts in coronary MRA using SSFP sequence can be markedly reduced with LCE acquisition of PE lines.     

Time-resolved contrast-enhanced coronary magnetic resonance angiography with highly constrained projection reconstruction

Contrast-enhanced magnetic resonance angiography (MRA) is a promising technique for coronary artery imaging. The blood signal changes during the contrast injection will result in image artifacts, blurring and relatively low signal-to-noise ratio, when the k-space segments from different cardiac cycles are combined to reconstruct the final image as “time averaged.” Thus, it is important to acquire data during maximal blood signal enhancement for first-pass, contrast-enhanced MRA, and relatively high temporal resolution is required. This work demonstrated the feasibility of highly constrained backprojection reconstruction for time-resolved, contrast-enhanced coronary MRA. With this method, the temporal resolution can be increased. In addition, coronary artery images around blood signal enhancement peak have significantly improved contrast-to-noise ratio and suppressed artifacts compared to the composite images which were collected during a much longer acquisition time during substantial blood signal changes.     

Scanning time efficient slinky for non-contrast MRA at low field.

To eliminate slab boundary artifact (SBA) for non-contrast multi-slab three-dimensional time-of-flight magnetic resonance angiogram (3D TOF MRA), we have previously developed a novel technique, termed SLINKY (Sliding Interleaved kY) acquisition in which a thin slab continuously "walks" along the z-axis while data are acquired in an interleaved fashion along the kY-axis. It has been demonstrated in our earlier works that SLINKY can suppress the SBA without any assumption of blood flow behavior, such as velocity or direction. At the same time, SLINKY keeps the same SNR as conventional multiple overlapping thin slab acquisition (MOTSA). Yet, this method is sensitive to any phase error along the ky axis. In our earlier application of SLINKY, we used navigator echoes to measure and correct the phase errors along the kY axis. The cost of using navigator echo collection is an increase in the imaging time. We therefore propose an improved SLINKY technique which does not use navigator echo collection for correcting phase errors, reducing the imaging time while keeping the same suppression of slab boundary artifacts. The present study demonstrates that by using a specifically designed RF pulse, the navigator echo collection can be avoided without incurring any extra ghosting or SNR reduction in the reconstructed images.     

Determination of vessel cross section for flow rate quantification

This study was motivated by the interest of measuring different cardiac parameters for which changes in the flow rate during a cardiac cycle needs to be determined at different positions along a vessel segment. These measurements result in a great number of images for which automatic contour detection is very helpful. A model-based algorithm for intraluminal contour detection has been developed in order to allow an accurate quantitative image analysis. The algorithm permits to select contours automatically on all the frames and slices of an imaging study. Images obtained on a flow phantom simulating the effects of blood circulation in large arteries have been used to validate the method. They were acquired with a specially designed interleaved multi slice and phase sequence, using a standard whole-body 2 Tesla NMR scanner. A potential in vivo application of the algorithm has been demonstrated on abdominal aorta images.     

High-resolution intracranial MRA at 7T using autocalibrating parallel imaging: initial experience in vascular disease patients

Purpose

Greater spatial resolution in intracranial three-dimensional time-of-flight (TOF) magnetic resonance angiography (MRA) is possible at higher field strengths, due to the increased contrast-to-noise ratio (CNR) from the higher signal-to-noise ratio and the improved background suppression. However, at very high fields, spatial resolution is limited in practice by the acquisition time required for sequential phase encoding. In this study, we applied parallel imaging to 7T TOF MRA studies of normal volunteers and patients with vascular disease, in order to obtain very high resolution (0.12 mm³) images within a reasonable scan time.

Materials and Methods

Custom parallel imaging acquisition and reconstruction methods were developed for 7T MRA, based on generalized autocalibrating partially parallel acquisition (GRAPPA). The techniques were compared and applied to studies of seven normal volunteers and three patients with cerebrovascular disease.

Results

The technique produced high resolution studies free from discernible reconstruction artifacts in all subjects and provided excellent depiction of vascular pathology in patients.

Conclusions

7T TOF MRA with parallel imaging is a valuable noninvasive angiographic technique that can attain very high spatial resolution.     

Quantitative multi-modal functional MRI with blood oxygenation level dependent exponential decays adjusted for flow attenuated inversion recovery (BOLDED AFFAIR)

A magnetic resonance imaging (MRI) method is described that allows interleaved measurements of transverse (R(2)(*) and R(2)) and longitudinal (R(1)) relaxation rates of tissue water in conjunction with spin labeling. The image-contrasts are intrinsically blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) weighted, but each contrast is made quantitative by two echo time (TE) and inversion recovery time (TIR) acquisitions with gradient echo (GE) and spin echo (SE) weighted echo-planar imaging (EPI). The EPI data were acquired at 7 Tesla with nominal spatial resolution of 430 x 430 x 1000 microm(3) in rat brain in vivo. The method is termed as blood oxygenation level dependent exponential decays adjusted for flow attenuated inversion recovery (BOLDED AFFAIR) and allows acquisition of R(2)(*), R(2), and CBF maps in an interleaved manner within approximately 12 minute. The basic theory of the method, associated experimental/systematic errors, and temporal restrictions are discussed. The method is validated by comparison of multi-modal maps obtained by BOLDED AFFAIR (i.e., two TE and TIR values with GE and SE sequences) and conventional approach (i.e., multiple TE and TIR values with GE and SE sequences) during varied levels of whole brain activity. Preliminary functional data from a rat forepaw stimulation model demonstrate the feasibility of this method for functional MRI (fMRI) studies. It is expected that with appropriate precautions this method in conjunction with contrast agent-based MRI has great potential for quantitative fMRI studies of mammalian cortex.     

Spectrum envelope analysis for interleaved sampled gratings with phase-shift

An analytical method is proposed for calculating reflection-spectrum envelope of interleaved sampled gratings with phase-shifts. Effects of parameters on reflection-spectrum envelope of the grating are discussed, such as two grating segment lengths, phase-shift and interleaved length. The unique optimal interleaved length can be obtained to provide the broadest flat-top reflection-spectrum envelope without ripples. Accuracy of the analytical expression is verified due to the calculated results in good agreement with the reflection spectra obtained by the piecewise uniform transform matrix method.     

Time-resolved and high-resolution MRA in a rabbit model of pulmonary embolism at 7 T: preliminary results

Purpose

Evaluate feasibility of using time-resolved and high-resolution, contrast-enhanced magnetic resonance angiography (MRA) at 7 T for characterization of an animal model of pulmonary embolism.

Methods

MRAs were performed in five rabbits using a 7-T MR scanner. Preceding the MR studies, each rabbit underwent a pulmonary artery catheterization with balloon placement. Two doses of gadodiamide were injected: first during a time-resolved MRA, immediately followed by a high-resolution acquisition. Balloon was then deflated, permitting reperfusion for 5 min. A second dose was then injected and another high-resolution MRA acquired. Measurements of signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and vessel cross-sections down to fourth-order branches were made, among other parameters.

Results

Occlusion was detected in all rabbits. Despite a TE of 0.58 ms for the time-resolved MRA, regions of nonuniform enhancement attributed to susceptibility effects at the 7-T field were observed in perfused lung. Mean SNR=7.5±3.3 and 134.2±46.5 for the lung and aorta, respectively, and mean CNR=126.7±46.4 for aorta versus lung were obtained. Diameters of vessels in lung that was never occluded were not statistically different from those in reperfused lung.

Conclusion

Results show that time-resolved and high-resolution MRA of the lung are feasible at 7 T and provide high SNR, CNR and resolution, but TEs smaller than 0.58 ms are required to avoid susceptibility artifacts in time-resolved MRAs.     

Fast 3D echo planar SSFP-based 1H spectroscopic imaging: demonstration on the rat brain in vivo

A fast proton spectroscopic imaging pulse sequence based on the condition of steady-state free precession is presented. High 3D spatial and temporal resolution is achieved using simultaneous detection of both one spatial and one spectral dimension, with a time-dependent gradient cycle known from echo planar imaging. Additionally, in order to increase the spectral width of the measurement, an interleaved acquisition scheme is shown either for systems with limited gradient switching capabilities or applications with a wide chemical shift range. The pulse sequence is implemented on a standard 4.7-T nuclear magnetic resonance animal imaging system. Measurements with a total measurement time of less than 2.5 min and a nominal voxel size of 6.75 microl using a total of 64 x 32 x 16 voxels are performed on phantoms and healthy rat brain in vivo allowing the rapid detection of signals from both uncoupled and J-coupled spin systems with high signal-to-noise ratio.     

Broad-band NMR with a high spectral and spatial resolution

Spatially selective excitation sequence CARVE (completely arbitrary regional volume excitation) excites signal from an arbitrarily shaped profile (I. Sersa, S. Macura: Magn. Reson. Med.37, 920–931, 1997) by an interleaved sequence of precalculated small tip angle radio-frequency pulses and gradient pulses. Here we propose a spatially selective observation method based on the CARVE principles which is insensitive to the relaxation and the off-resonance effects. The method, CARVED (CARVE detection), excites spins uniformly across the sample and across the spectrum but achieves spatial selectivity by weighted coaddition of the signals after the data acquisition. CARVE-D is suitable for spatially selective high-resolution nuclear magnetic resonance spectroscopy in chemically and geometrically complex systems. The method is analyzed theoretically and demonstrated experimentally on model systems.     

Spectral velocity estimation in ultrasound using sparse data sets

Velocity distributions in blood vessels can be displayed using ultrasound scanners by making a Fourier transform of the received signal and then showing spectra in an M-mode display. It is desired to show a B-mode image for orientation, and data for this have to be acquired interleaved with the flow data. This either halves the effective pulse repetition frequency f(prf) or gaps appear in the spectrum from B-mode emissions. This paper presents a technique to maintain the highest possible f(prf) and at the same time show a B-mode image. The power spectrum can be calculated from the Fourier transform of the autocorrelation function, and it is shown that the autocorrelation function can be calculated for a sparse set of data where flow and B-mode emissions are interspaced. Both short deterministic sequences of emissions and full random sequences can be used. The dynamic range of the sparse sequence is reduced compared to a full sequence. Typically, a reduction of 20 dB is found when using 66% of the data compared to using all data. The theory of the method and examples from simulations of flow in arteries are presented. The audio signal can also be generated from the spectrogram.