Search

Purpose: Conventional non-Cartesian compressed sensing requires multiple nonuniform Fourier transforms every iteration, which is computationally expensive. Accordingly, time-consuming reconstructions have slowed the adoption of undersampled 3D non-Cartesian acquisitions into clinical protocols. In this work we investigate several approaches to minimize reconstruction times without sacrificing accuracy., Methods: The reconstruction problem can be reformatted to exploit the Toeplitz structure of matrices that are evaluated every iteration, but it requires larger oversampling than what is strictly required by nonuniform Fourier transforms. Accordingly, we investigate relative speeds of the two approaches for various nonuniform Fourier transform kernel sizes and oversampling for both GPU and CPU implementations. Second, we introduce a method to minimize matrix sizes by estimating the image support. Finally, density compensation weights have been used as a preconditioning matrix to improve convergence, but this increases noise. We propose a more general approach to preconditioning that allows a trade-off between accuracy and convergence speed., Results: When using a GPU, the Toeplitz approach was faster for all practical parameters. Second, it was found that properly accounting for image support can prevent aliasing errors with minimal impact on reconstruction time. Third, the proposed preconditioning scheme improved convergence rates by an order of magnitude with negligible impact on noise., Conclusion: With the proposed methods, 3D non-Cartesian compressed sensing with clinically relevant reconstruction times (<2 min) is feasible using practical computer resources. Magn Reson Med 79:2685-2692, 2018. © 2017 International Society for Magnetic Resonance in Medicine., (© 2017 International Society for Magnetic Resonance in Medicine.)

Purpose: To create a B 0 map and correct for off-resonance with minimal scan time increase for two-dimensional (2D) or 3D non-Cartesian acquisitions., Methods: Rewinding trajectories that bring the zeroth gradient moment to zero every repetition time (TR) were used to estimate the off-resonance with a center-out 3D cones trajectory, which required an increase in the minimum TR by 5%. The off-resonance estimation and correction was implemented using an algorithm based on binning and object-domain phase correction. B 0 maps using BMART (B 0 mapping using rewinding trajectories) were compared to maps obtained using separate scans with multiple echo time (TE) in a phantom and human brain., Results: Excellent agreement between BMART and the multiple-TE method were observed, and images corrected with BMART were deblurred., Conclusion: BMART can correct for off-resonance without requiring an additional scan, and can be easily applied to center-out or projection trajectories (2D or 3D). Magn Reson Med 78:664-669, 2017. © 2016 International Society for Magnetic Resonance in Medicine., (© 2016 International Society for Magnetic Resonance in Medicine.)

Purpose: We propose a method to acquire B1 distribution plots by encoding in B1 instead of image space. Using this method, B1 data is acquired in a different way from traditional spatial B1 mapping, and allows for quick measurement of high dynamic range B1 data., Methods: To encode in B1, we acquire multiple projections of a slice, each along the same direction, but using a different phase sensitivity to B1. Using a convex optimization formulation, we reconstruct histograms of the B1 distribution estimates of the slice., Results: We verify in vivo B1 distribution measurements by comparing measured distributions to distributions calculated from reference spatial B1 maps using the Earth Mover's Distance. Phantom measurements using a surface coil show that for increased spatial B1 variations, measured B1 distributions using the proposed method more accurately estimate the distribution than a low-resolution spatial B1 map, resulting in a 37% Earth Mover's Distance decrease while using fewer measurements., Conclusion: We propose and validate the performance of a method to acquire B1 distribution information directly without acquiring a spatial B1 map. The method may provide faster estimates of a B1 field for applications that do not require spatial B1 localization, such as the transmit gain calibration of the scanner, particularly for high dynamic B1 ranges. Magn Reson Med 77:229-236, 2017. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

Purpose: To develop a robust motion estimation method for free-breathing body MRI using dense coil arrays., Methods: Self-navigating pulse sequences can measure subject motion without using external motion monitoring devices. With dense coil arrays, individual coil elements can provide localized motion estimates. An averaged motion estimate over all coils is often used for motion compensation. However, this motion estimate may not accurately represent the dominant motion within the imaging volume. In this work, a coil clustering method is proposed to automatically determine the dominant motion for dense coil arrays. The feasibility of the proposed method is investigated in free-breathing abdominal MRI and cardiac MRI, and compared with manual motion estimate selection for respiratory motion estimation and electrocardiography for cardiac motion estimation., Results: Automated motion estimation achieved similar respiratory motion estimation compared to manual selection (averaged correlation coefficient 0.989 and 0.988 for abdominal MRI and cardiac MRI, respectively), and accurate cardiac triggering compared to electrocardiography (averaged temporal variability 17.5 ms)., Conclusion: The proposed method can provide accurate automated motion estimation for body MRI using dense coil arrays. It can enable self-navigated free-breathing abdominal and cardiac MRI without the need for external motion monitoring devices. Magn Reson Med 76:197-205, 2016. © 2015 Wiley Periodicals, Inc., (© 2015 Wiley Periodicals, Inc.)

Purpose: Accurate measurement of the nonuniform transmit radiofrequency field is necessary for magnetic resonance imaging applications. The radiofrequency field excitation amplitude (B1) is often obtained by acquiring a B1 map. We modify the B1 estimation using adiabatic refocusing (BEAR) method to extend its range to lower B1 magnitudes., Theory and Methods: The BEAR method is a phase-based B1 mapping method, wherein hyperbolic secant pulses induce a phase sensitivity to B1. The measurable B1 range is limited due to the adiabatic threshold of the pulses. We redesign the method to use flattened hyperbolic secant pulses, which have lower adiabatic thresholds. We optimize the flattened hyperbolic secant parameters to minimize phase sensitivity to frequency variations., Results: We validate the performance of the new method via simulation and in vivo at 3T, and show that for n ≤ 8, accurate B1 maps can be acquired using reduced nominal peak B1 values., Conclusion: The adiabatic threshold for the BEAR method is reduced with flattened hyperbolic secant pulses, which are optimized for accurate phase-to-B1 mapping over a frequency range, and allow for lower nominal B1 values. At 3T, the nominal B1 is decreased by 52% and the sensitivity to B1 is increased by a factor of 3.8. This can improve the method's applicability for measurement of low B1., (© 2015 Wiley Periodicals, Inc.)

Purpose: Accurate measurement of the nonuniform transmit radiofrequency field is useful for many applications in magnetic resonance imaging, such as calibrating the scanner transmit system, evaluating coil performance, and improving image quality and quantitation. The radiofrequency field excitation amplitude (B(1)) is often obtained by acquiring a B(1) map. In this study, a new B(1) mapping method is proposed., Theory and Methods: The use of two adiabatic full passage pulses with different magnitudes applied as successive refocusing pulses results in a linear relationship between phase and B(1) field strength that is insensitive to the repetition time, off-resonance effects, T(1), and T(2). Using this method, B(1) mapping can be localized to a slice or three-dimensional (3D) volume, with a spin-echo acquisition that is appropriate for fast projection measurements., Results: This new method is shown to agree well with the Bloch-Siegert B(1) mapping method for both phantom and in vivo B(1) measurements at 1.5T, 3T, and 7T. The method's ability to acquire accurate projection B(1) measurements is also demonstrated., Conclusion: This method's high dynamic range, ability to make fast projection measurements, and linear quantitative relationship between phase and B1 make it an ideal candidate for use in robust transmitter gain calibration., (Copyright © 2013 Wiley Periodicals, Inc.)

Purpose: The balanced steady-state free precession (bSSFP) pulse sequence has shown to be of great interest due to its high signal-to-noise ratio efficiency. However, bSSFP images often suffer from banding artifacts due to off-resonance effects, which we aim to minimize in this article., Methods: We present a general and fast two-step algorithm for 1) estimating the unknowns in the bSSFP signal model from multiple phase-cycled acquisitions, and 2) reconstructing band-free images. The first step, linearization for off-resonance estimation (LORE), solves the nonlinear problem approximately by a robust linear approach. The second step applies a Gauss-Newton algorithm, initialized by LORE, to minimize the nonlinear least squares criterion. We name the full algorithm LORE-GN., Results: We derive the Cramér-Rao bound, a theoretical lower bound of the variance for any unbiased estimator, and show that LORE-GN is statistically efficient. Furthermore, we show that simultaneous estimation of T1 and T2 from phase-cycled bSSFP is difficult, since the Cramér-Rao bound is high at common signal-to-noise ratio. Using simulated, phantom, and in vivo data, we illustrate the band-reduction capabilities of LORE-GN compared to other techniques, such as sum-of-squares., Conclusion: Using LORE-GN we can successfully minimize banding artifacts in bSSFP., (Copyright © 2013 Wiley Periodicals, Inc.)

Purpose: To develop new magnetization-prepared imaging schemes based on a three-dimensional (3D) concentric cylinders trajectory., Methods: The 3D concentric cylinders trajectory, which is robust to off-resonance effects and timing delays while requiring fewer excitations than a comparable 3D Cartesian (3DFT) sequence, is used as the readout for magnetization-prepared sequences exploiting its inherently centric-ordered structure. Two applications: (i) T1 -weighted brain imaging with an inversion-recovery-prepared radiofrequency-spoiled gradient-echo (IR-SPGR) sequence, (ii) non-contrast-enhanced (NCE) peripheral angiography with a magnetization-prepared balanced steady-state free precession (bSSFP) sequence are presented to demonstrate the effectiveness of the proposed method. For peripheral angiography, the scan efficiency is further improved by interleaving different preparations at different rates and by carefully designing the sampling geometry for an efficient parallel imaging method., Results: In vivo brain scans with an IR-SPGR sequence and lower extremity scans with a magnetization-prepared bSSFP sequence for NCE peripheral angiography both demonstrate that the proposed sequences with concentric cylinders effectively capture the transient magnetization-prepared contrast with faster scan times than a corresponding 3DFT sequence. The application of peripheral angiography also shows the feasibility of the proposed interleaving schemes and parallel imaging method., Conclusion: The 3D concentric cylinders trajectory is a robust and efficient readout that is well-suited for magnetization-prepared imaging., (Copyright © 2013 Wiley Periodicals, Inc.)

Purpose: To develop a new velocity-selective (VS) excitation pulse sequence which is robust to field inhomogeneity, and demonstrate its application to non-contrast-enhanced peripheral MR angiography (MRA)., Methods: The off-resonance-robust VS saturation pulse is designed by incorporating 180° refocusing pulses into the k-space-based reference design and tailoring sequence parameters in a velocity region of interest. The VS saturation pulse is used as magnetization preparation for non-contrast-enhanced peripheral MRA to suppress background tissues but not arterial blood based on their velocities. Non-contrast-enhanced peripheral MRA using the proposed VS preparation was tested in healthy volunteers and a patient with arterial stenosis., Results: Calf angiograms obtained using the new VS preparation show more uniform background suppression than the reference VS preparation, as demonstrated by larger mean values and smaller standard deviations of artery-to-vein and artery-to-muscle contrast-to-noise ratios (71.0 ± 11.4 and 75.3 ± 12.1 versus 61.7 ± 22.7 and 58.5 ± 27.8). Two-station peripheral MRA using the new VS preparation identifies stenosis of the femoral and popliteal arteries in the patient, as validated by digital subtraction angiography., Conclusion: Non-contrast-enhanced MRA using the new VS magnetization preparation can reliably provide high angiographic contrast in the lower extremities with significantly improved immunity to field inhomogeneity., (Copyright © 2012 Wiley Periodicals, Inc.)

Fast imaging trajectories are used in MRI to speed up the acquisition process, but imperfections in the gradient system create artifacts in the reconstructed images. Artifacts result from the deviation between k-space trajectories achieved on the scanner and their original prescription. Measuring or approximating actual k-space trajectories with predetermined gradient timing delays reduces the artifacts, but are generally based on a specific trajectory and scan orientation. A single linear time-invariant characterization of the gradient system provides a method to predict k-space trajectories scanned in arbitrary orientations through convolution. This is done efficiently, by comparing the Fourier transforms of the input and measured waveforms of a single high-bandwidth test gradient waveform. This new method is tested for spiral, interleaved echo-planar, and three-dimensional cones imaging, demonstrating its ability to reduce reconstructed image artifacts for various k-space trajectories., (Copyright © 2011 Wiley-Liss, Inc.)

Small perturbations of steady-state sequence parameters can induce very large spectral profile deviations that are localized to specific off-resonant frequencies, denoted critical frequencies. Although, a small number of studies have previously considered the use of these highly specific modulations for MR angiography and elastography, many potential applications still remain to be explored. An analysis of this phenomenon using a linear systems technique and a geometric magnetization trajectory technique shows that the critical frequencies correspond to singularities in the steady-state signal equation. An interleaved acquisition combined with a complex difference technique yields a spectral profile containing sharp peaks interleaved with wide stopbands, while a complex sum technique yields a spectral profile similar to that of balanced steady-state free precession. Simulations and phantom experiments are used to demonstrate a novel application of this technique for positive-contrast imaging of superparamagnetic iron-oxide nanoparticles. The technique is shown to yield images with high levels of positive contrast and good water and fat background suppression. The technique can also simultaneously yield images with contrast similar to balanced steady-state free precession., (Copyright © 2011 Wiley Periodicals, Inc.)

The article "A k-space analysis of small-tip-angle excitation" introduced a spatial frequency interpretation of the effect of RF excitation pulses. This introduction describes where the initial ideas for this paper came from, and traces out some of the applications that have been developed using this perspective., (Copyright © 2011. Published by Elsevier Inc.)

Motion--both rigid-body and nonrigid--is the main limitation to in vivo, high-resolution larynx imaging. In this work, a new real-time motion compensation algorithm is introduced. Navigator data are processed in real time to compute the displacement information, and projections are corrected using phase modulation in k-space. Upon automatic feedback, the system immediately reacquires the data most heavily corrupted by nonrigid motion, i.e., the data whose corresponding projections could not be properly corrected. This algorithm overcomes the shortcomings of the so-called diminishing variance algorithm by combining it with navigator-based rigid-body motion correction. Because rigid-body motion correction is performed first, continual bulk motion no longer impedes nor prevents the convergence of the algorithm. Phantom experiments show that the algorithm properly corrects for translations and reacquires data corrupted by nonrigid motion. Larynx imaging was performed on healthy volunteers, and substantial reduction of motion artifacts caused by bulk shift, swallowing, and coughing was achieved., (Copyright © 2011 Wiley-Liss, Inc.)

Magnetization-prepared acquisitions offer a trade-off between image contrast and scan efficiency for magnetic resonance imaging. Because the prepared signals gradually decay, the contrast can be improved by frequently repeating the preparation, which in turn significantly increases the scan time. A common solution is to perform the data collection progressing from low- to high-spatial-frequency samples following each preparation. Unfortunately, this leads to loss of spatial resolution, and thereby image blurring. In this work, a new technique is proposed that first corrects the signal decay in high-frequency data to mitigate the resolution loss and improve the image contrast without reducing the scan efficiency. The proposed technique then employs a sparsity-based nonlinear reconstruction to further improve the image quality. In addition to reducing the amplified high-frequency noise, this reconstruction extrapolates missing k-space samples in the case of undersampled compressed-sensing acquisitions. The technique is successfully demonstrated for noncontrast-enhanced flow-independent angiography of the lower extremities, an application that substantially benefits from both the signal compensation and the nonlinear reconstruction.

Background: Coronary wall cardiovascular magnetic resonance (CMR) is a promising noninvasive approach to assess subclinical atherosclerosis, but data are limited in subjects over 60 years old, who are at increased risk. The purpose of the study was to evaluate coronary wall CMR in an asymptomatic older cohort., Results: Cross-sectional images of the proximal right coronary artery (RCA) were acquired using spiral black-blood coronary CMR (0.7 mm resolution) in 223 older, community-based patients without a history of cardiovascular disease (age 60-72 years old, 38% female). Coronary measurements (total vessel area, lumen area, wall area, and wall thickness) had small intra- and inter-observer variabilities (r = 0.93~0.99, all p < 0.0001), though one-third of these older subjects had suboptimal image quality. Increased coronary wall thickness correlated with increased coronary vessel area (p < 0.0001), consistent with positive remodeling. On multivariate analysis, type 2 diabetes was the only risk factor associated with increased coronary wall area and thickness (p = 0.03 and p = 0.007, respectively). Coronary wall CMR measures were also associated with coronary calcification (p = 0.01-0.03)., Conclusions: Right coronary wall CMR in asymptomatic older subjects showed increased coronary atherosclerosis in subjects with type 2 diabetes as well as coronary calcification. Coronary wall CMR may contribute to the noninvasive assessment of subclinical coronary atherosclerosis in older, at-risk patient groups.

In this article, a robust methodology for in vivo T(1) mapping is presented. The approach combines a gold standard scanning procedure with a novel fitting procedure. Fitting complex data to a five-parameter model ensures accuracy and precision of the T(1) estimation. A reduced-dimension nonlinear least squares method is proposed. This method turns the complicated multi-parameter minimization into a straightforward one-dimensional search. As the range of possible T(1) values is known, a global grid search can be used, ensuring that a global optimal solution is found. When only magnitude data are available, the algorithm is adapted to concurrently restore polarity. The performance of the new algorithm is demonstrated in simulations and phantom experiments. The new algorithm is as accurate and precise as the conventionally used Levenberg-Marquardt algorithm but much faster. This gain in speed makes the use of the five-parameter model viable. In addition, the new algorithm does not require initialization of the search parameters. Finally, the methodology is applied in vivo to conventional brain imaging and to skin imaging. T(1) values are estimated for white matter and gray matter at 1.5 T and for dermis, hypodermis, and muscle at 1.5 T, 3 T, and 7 T.

Efficient acquisition strategies for magnetization-prepared imaging based on the three-dimensional (3D) stack-of-rings k-space trajectory are presented in this work. The 3D stack-of-rings can be acquired with centric ordering in all three dimensions for greater efficiency in capturing the desired contrast. In addition, the 3D stack-of-rings naturally supports spherical coverage in k-space for shorter scan times while achieving isotropic spatial resolution. While non-Cartesian trajectories generally suffer from greater sensitivity to system imperfections, the 3D stack-of-rings can enhance magnetization-prepared imaging with a high degree of robustness to timing delays and off-resonance effects. As demonstrated with phantom scans, timing errors and gradient delays only cause a bulk rotation of the 3D stack-of-rings reconstruction. Furthermore, each ring can be acquired with a time-efficient retracing design to resolve field inhomogeneities and enable fat/water separation. To demonstrate its effectiveness, the 3D stack-of-rings are considered for the case of inversion-recovery-prepared structural brain imaging. Experimental results show that the 3D stack-of-rings can achieve higher signal-to-noise ratio and higher contrast-to-noise ratio within a shorter scan time when compared to the standard inversion-recovery-prepared sequence based on 3D Cartesian encoding. The design principles used for this specific case of inversion-recovery-prepared brain imaging can be applied to other magnetization-prepared imaging applications., ((c) 2010 Wiley-Liss, Inc.)

Partially parallel imaging with localized sensitivities is a fast parallel image reconstruction method for both Cartesian and non-Cartesian trajectories, but suffers from aliasing artifacts when there are deviations from the assumption of perfect localization. Such reconstructions would normally crop the individual coil images to remove the artifacts prior to combination. However, the sampling densities in variable-density k-space trajectories support different field-of-views for separate regions in k -space. In fact, the higher sampling density of low frequencies can be used to reconstruct a bigger field-of-view without introducing aliasing artifacts and the resulting image signal-to-noise ratio (SNR) can be improved. A novel, fast variable-density parallel imaging method is presented, which reconstructs different field-of-views from separate frequencies according to the local sampling density in k-space. Aliasing-suppressed images can be produced with high SNR-efficiency without the need for accurate estimation of coil sensitivities and complex or iterative computations.

As a noninvasive modality, MR is attractive for in vivo skin imaging. Its unique soft tissue contrast makes it an ideal imaging modality to study the skin water content and to resolve the different skin layers. In this work, the challenges of in vivo high-resolution skin imaging are addressed. Three 3D Cartesian sequences are customized to achieve high-resolution imaging and their respective performance is evaluated. The balanced steady-state free precession (bSSFP) and gradient echo (GRE) sequences are fast but can be sensitive to off-resonance artifacts. The fast large-angle spin echo (FLASE) sequence provides a sharp depiction of the hypodermis structures but results in more specific absorption rate (SAR). The effect of increasing the field strength is assessed. As compared to 1.5 T, signal-to-noise ratio at 3 T slightly increases in the hypodermis and almost doubles in the dermis. The need for fat/water separation is acknowledged and a solution using an interleaved three-point Dixon method and an iterative reconstruction is shown to be effective. The effects of motion are analyzed and two techniques to prevent motion and correct for it are evaluated. Images with 117 x 117 x 500 microm(3) resolution are obtained in imaging times under 6 min., ((c) 2010 Wiley-Liss, Inc.)

There has been recent interest in positive-contrast MRI methods for noninvasive tracking of cells labeled with superparamagnetic iron-oxide nanoparticles. Low-tip-angle balanced steady-state free precession sequences have been used for fast, high-resolution, and flow-insensitive positive-contrast imaging; however, the contrast can be compromised by the limited suppression of the on-resonant and fat signals. In this work, a new technique that produces positive contrast with alternating repetition time steady-state free precession is proposed to achieve robust background suppression for a broad range of tissue parameters. In vitro and in vivo experiments demonstrate the reliability of the generated positive contrast. The results indicate that the proposed method can enhance the suppression level by up to 18 dB compared with conventional balanced steady-state free precession.

Although balanced steady-state free precession (bSSFP) imaging yields high signal-to-noise ratio (SNR) efficiency, the bright lipid signal is often undesirable. The bSSFP spectrum can be shaped to suppress the fat signal with scan-efficient alternating repetition time (ATR) bSSFP. However, the level of suppression is limited, and the pass-band is narrow due to its nonuniform shape. A multiple repetition time (TR) bSSFP scheme is proposed that creates a broad stop-band with a scan efficiency comparable with ATR-SSFP. Furthermore, the pass-band signal uniformity is improved, resulting in fewer shading/banding artifacts. When data acquisition occurs in more than a single TR within the multiple-TR period, the echoes can be combined to significantly improve the level of suppression. The signal characteristics of the proposed technique were compared with bSSFP and ATR-SSFP. The multiple-TR method generates identical contrast to bSSFP, and achieves up to an order of magnitude higher stop-band suppression than ATR-SSFP. In vivo studies at 1.5 T and 3 T demonstrate the superior fat-suppression performance of multiple-TR bSSFP., ((c) 2009 Wiley-Liss, Inc.)

The Periodically Rotated Overlapping ParallEL Lines with Enhanced Reconstruction (PROPELLER) method for magnetic resonance imaging data acquisition and reconstruction has the highly desirable property of being able to correct for motion during the scan, making it especially useful for imaging pediatric or uncooperative patients and diffusion imaging. This method nominally supports a circular field of view (FOV), but tailoring the FOV for noncircular shapes results in more efficient, shorter scans. This article presents new algorithms for tailoring PROPELLER acquisitions to the desired FOV shape and size that are flexible and precise. The FOV design also allows for rotational motion which provides better motion correction and reduced aliasing artifacts. Some possible FOV shapes demonstrated are ellipses, ovals and rectangles, and any convex, pi-symmetric shape can be designed. Standard PROPELLER reconstruction is used with minor modifications, and results with simulated motion presented confirm the effectiveness of the motion correction with these modified FOV shapes. These new acquisition design algorithms are simple and fast enough to be computed for each individual scan. Also presented are algorithms for further scan time reductions in PROPELLER echo-planar imaging (EPI) acquisitions by varying the sample spacing in two directions within each blade.

The concentric rings two-dimensional (2D) k-space trajectory enables flexible trade-offs between image contrast, signal-to-noise ratio (SNR), spatial resolution, and scan time. However, to realize these benefits for in vivo imaging applications, a robust method is desired to deal with fat signal in the acquired data. Multipoint Dixon techniques have been shown to achieve uniform fat suppression with high SNR-efficiency for Cartesian imaging, but application of these methods for non-Cartesian imaging is complicated by the fact that fat off-resonance creates significant blurring artifacts in the reconstruction. In this work, two fat-water separation algorithms are developed for the concentric rings. A retracing design is used to sample rings near the center of k-space through multiple revolutions to characterize the fat-water phase evolution difference at multiple time points. This acquisition design is first used for multipoint Dixon reconstruction, and then extended to a spectroscopic approach to account for the trajectory's full evolution through 3D k-t space. As the trajectory is resolved in time, off-resonance effects cause shifts in frequency instead of spatial blurring in 2D k-space. The spectral information can be used to assess field variation and perform robust fat-water separation. In vivo experimental results demonstrate the effectiveness of both algorithms.

Balanced SSFP achieves high SNR efficiency, but suffers from bright fat signal. In this work, a multiple-acquisition fat-water separation technique using alternating repetition time (ATR) balanced SSFP is proposed. The SSFP profile can be modified using alternating repetition times and appropriate phase cycling to yield two spectra where fat and water are in-phase and out-of-phase, respectively. The signal homogeneity and the broad width of the created in-phase and out-of-phase profiles lead to signal cancellation over a broad stop-band. The stop-band suppression is achieved for a wide range of flip angles and tissue parameters. This property, coupled with the inherent flexibility of ATR SSFP in repetition time selection, makes the method a good candidate for fat-suppressed SSFP imaging. The proposed method can be tailored to achieve a smaller residual stop-band signal or a decreased sensitivity to field inhomogeneity depending on application-specific needs., ((c) 2008 Wiley-Liss, Inc.)

Passband balanced-steady-state free precession (b-SSFP) fMRI is a recently developed method that utilizes the passband (flat portion) of the b-SSFP off-resonance response to measure MR signal changes elicited by changes in tissue oxygenation following increases in neuronal activity. Rapid refocusing and short readout durations of b-SSFP, combined with the relatively large flat portion of the b-SSFP off-resonance spectrum allows distortion-free full-brain coverage with only two acquisitions. This allows for high-resolution functional imaging, without the spatial distortion frequently encountered in conventional high-resolution functional images. Finally, the 3D imaging compatibility of the b-SSFP acquisitions permits isotropic-voxel-size high-resolution acquisitions. In this study we address some of the major technical issues involved in obtaining passband b-SSFP-based functional brain images with practical imaging parameters and demonstrate the advantages through breath-holding and visual field mapping experiments., ((c) 2008 Wiley-Liss, Inc.)

Nonuniform, non-Cartesian k-space trajectories enable fast scanning with reduced motion and flow artifacts. In such cases, the data are usually convolved with a kernel and resampled onto a Cartesian grid before reconstruction. For trajectories such as undersampled variable-density spirals, the mainlobe width of the kernel for undersampled high spatial frequencies has to be larger to limit the amount of aliasing energy. Continuously varying the kernel extent is time consuming. By dividing k-space into several annuli and using appropriate mainlobe widths for each, the aliasing energy and noise can be reduced at the expense of lower resolution towards the edge of the field of view (FOV). Resolution can instead be preserved at the center of the FOV, which is expected to be free of artifacts, without any artifact reduction. The image reconstructed from each annulus can be deapodized separately. The method can be applied to most k-space trajectories used in MRI., (2007 Wiley-Liss, Inc)

The concentric rings two-dimensional (2D) k-space trajectory provides an alternative way to sample polar data. By collecting 2D k-space data in a series of rings, many unique properties are observed. The concentric rings are inherently centric-ordered, provide a smooth weighting in k-space, and enable shorter total scan times. Due to these properties, the concentric rings are well-suited as a readout trajectory for magnetization-prepared studies. When non-Cartesian trajectories are used for MRI, off-resonance effects can cause blurring and degrade the image quality. For the concentric rings, off-resonance blur can be corrected by retracing rings near the center of k-space to obtain a field map with no extra excitations, and then employing multifrequency reconstruction. Simulations show that the concentric rings exhibit minimal effects due to T(2) (*) modulation, enable shorter scan times for a Nyquist-sampled dataset than projection-reconstruction imaging or Cartesian 2D Fourier transform (2DFT) imaging, and have more spatially distributed flow and motion properties than Cartesian sampling. Experimental results show that off-resonance blurring can be successfully corrected to obtain high-resolution images. Results also show that concentric rings effectively capture the intended contrast in a magnetization-prepared sequence., (2007 Wiley-Liss, Inc)

Balanced steady-state free precession (SSFP) is hindered by the inherent off-resonance sensitivity and unwanted bright fat signal. Multiple-acquisition SSFP combination methods, where multiple datasets with different fixed RF phase increments are acquired, have been used for shaping the SSFP spectrum to solve both problems. We present a new combination method (weighted-combination SSFP or WC-SSFP) that preserves SSFP contrast and enables banding-reduction and fat-water separation. Methods addressing the banding artifact have focused on either getting robust banding-reduction (complex-sum SSFP) or improved SNR efficiency (sum-of-squares SSFP). The proposed method achieves both robust banding-reduction and an SNR efficiency close to that of the sum-of-squares method. A drawback of fat suppression methods that create a broad stop-band around the fat resonance is the wedge shape of the stop-band leading to imperfect suppression. WC-SSFP improves the suppression of the stop-band without affecting the pass-band performance, and prevents fat signal from obscuring the tissues of interest in the presence of considerable resonant frequency variations. The method further facilitates the use of SSFP imaging by providing a control parameter to adjust the level of banding-reduction or fat suppression to application-specific needs., ((c) 2007 Wiley-Liss, Inc.)

Purpose: To describe and evaluate a fast, fluid-suppressed 2D multislice steady-state free precession (SSFP) neuroimaging sequence., Materials and Methods: We developed a fast fluid-attenuated inversion-recovery SSFP sequence for use in neuroimaging. The inversion time (TI) was optimized to yield good cerebrospinal fluid (CSF) suppression while conserving white matter (WM)/lesion contrast across a broad range of flip angles. Multiple SSFP acquisitions were combined using the sum-of-squares (SOS) method to maximize SNR efficiency while minimizing SSFP banding artifacts. We compared our fluid-attenuated inversion-recovery (FLAIR) SSFP sequence with FLAIR fast spin-echo (FSE) in both normal subjects and a volunteer with multiple sclerosis. SNR measurements were performed to ascertain the SNR efficiency of each sequence., Results: Our FLAIR SSFP sequence demonstrated excellent CSF suppression and good gray matter (GM)/WM contrast. Coverage of the entire brain (5-mm slices, 24-cm FOV, 256 x 192 matrix) was achieved with FLAIR SSFP in less than half the scan time of a corresponding FLAIR FSE sequence with similar SNR, yielding improvements of more than 50% in SNR efficiency. Axial scans of a volunteer with multiple sclerosis show clearly visible plaques and very good visualization of brain parenchyma., Conclusion: We have demonstrated the feasibility of a very fast fluid-suppressed neuroimaging technique using SSFP., ((c) 2006 Wiley-Liss, Inc.)

Nanocrystals with advanced magnetic or optical properties have been actively pursued for potential biological applications, including integrated imaging, diagnosis and therapy. Among various magnetic nanocrystals, FeCo has superior magnetic properties, but it has yet to be explored owing to the problems of easy oxidation and potential toxicity. Previously, FeCo nanocrystals with multilayered graphitic carbon, pyrolytic carbon or inert metals have been obtained, but not in the single-shelled, discrete, chemically functionalized and water-soluble forms desired for biological applications. Here, we present a scalable chemical vapour deposition method to synthesize FeCo/single-graphitic-shell nanocrystals that are soluble and stable in water solutions. We explore the multiple functionalities of these core-shell materials by characterizing the magnetic properties of the FeCo core and near-infrared optical absorbance of the single-layered graphitic shell. The nanocrystals exhibit ultra-high saturation magnetization, r1 and r2 relaxivities and high optical absorbance in the near-infrared region. Mesenchymal stem cells are able to internalize these nanoparticles, showing high negative-contrast enhancement in magnetic-resonance imaging (MRI). Preliminary in vivo experiments achieve long-lasting positive-contrast enhancement for vascular MRI in rabbits. These results point to the potential of using these nanocrystals for integrated diagnosis and therapeutic (photothermal-ablation) applications.

3D (kx, ky, t)-space analysis is invoked to provide insights into dynamic magnetic resonance imaging (DMRI) with arbitrary k-space sampling trajectories. The effects of 3D sampling in (kx, ky, t) are analyzed theoretically and verified with computer simulation. The analyses show that a 3D sampling pattern that is more isotropic in (kx, ky, t)-space and denser around the (kx=0, ky=0, t)-line likely results in improved temporal qualities of DMRI. The isotropy of 3D sampling is quantified using Voronoi 3D cells and the isoperimetric theorem. This 3D sampling perspective provides a theoretical framework for understanding the temporal qualities of various DMRI methods.

An ungated spiral phase-contrast (USPC) method was used to measure cardiac output (CO) rapidly and conveniently. The USPC method, which was originally designed for small peripheral vessels, was modified to assess CO by measuring flow in the ascending aorta (AA). The modified USPC used a 12-interleaf spiral trajectory to acquire full-image data every 283 ms with 2-mm spatial resolution. The total scan time was 5 s. For comparison, a triggered real-time (TRT) method was used to indirectly calculate CO by measuring left-ventricular (LV) volume. The USPC and TRT measurements from all normal volunteers agreed. In a patient with patent ductus arteriosus (PDA), high CO was measured with USPC, which agreed well with the invasive cardiac-catheterized measurement. In normal volunteers, CO dropped about 20-30% with Valsalva maneuvering, and increased about 100% after exercise. Continuous 28-s cycling between Valsalva maneuvering and free-breathing showed that USPC can temporally resolve physiological CO changes., (Copyright 2006 Wiley-Liss, Inc.)

The 3D Cones k-space trajectory has many desirable properties for rapid and ultra-short echo time magnetic resonance imaging. An algorithm is presented that generates the 3D Cones gradient waveforms given a desired field of view and resolution. The algorithm enables a favorable trade-off between increases in readout time and decreases in the total number of required readouts. The resulting trajectory is very signal-to-noise ratio (SNR) efficient and has excellent aliasing properties. A rapid high-resolution ultra-short echo time imaging sequence is used to compare the 3D Cones trajectory to 3D projection reconstruction (3DPR) sampling schemes. For equivalent scan times, the 3D Cones trajectory has better SNR performance and fewer aliasing artifacts as compared to the 3DPR trajectory., (Magn Reson Med, 2006. (c) 2006 Wiley-Liss, Inc.)

Balanced steady-state free precession (SSFP) sequences use fully re-focussed gradient waveforms to achieve a high signal and useful image contrast in short scan times. Despite these strengths, the clinical feasibility of balanced SSFP is still limited both by bright fat signal and by the signal voids that result from off-resonance effects such as field or susceptibility variations. A new method, dual-acquisition phase-sensitive SSFP, combines the signals from two standard balanced SSFP acquisitions to separate fat and water while simultaneously reducing the signal voids. The acquisitions are added in quadrature and then phase corrected using a simple algorithm before fat and water can be identified simply by the sign of the signal. This method is especially useful for applications at high field, where the RF power deposition, spatial resolution requirements and gradient strength limit the minimum repetition times. Finally, dual-acquisition phase-sensitive SSFP can be combined with other magnetization preparation schemes to produce specific image contrast in addition to separating fat and water signals.

Multiple receiver-coil data collection is an effective approach to reduce scan time. There are many parallel imaging techniques that reduce scan time using multiple receiver coils. One of these methods, partially parallel imaging with localized sensitivities (PILS), utilizes the localized sensitivity of each coil. The advantages of PILS over other parallel imaging methods include the simplicity of the algorithm, good signal-to-noise ratio (SNR) properties, and the fact that there is no additional complexity involved in applying the algorithm to arbitrary k-space trajectories. This PILS method can be further improved to provide truly parallel broadband imaging with the use of multiple-demodulation hardware. By customizing the demodulation based on each coil's location, the k-space sampling rate can be chosen based on each coil's localized sensitivity region along the readout direction. A simulated demodulation of data from 2D Fourier transform (FT) and spiral trajectories is shown to demonstrate the method's feasibility., (Copyright (c) 2005 Wiley-Liss, Inc.)

Reconstruction of magnetic resonance images from data not falling on a Cartesian grid is a Fourier inversion problem typically solved using convolution interpolation, also known as gridding. Gridding is simple and robust and has parameters, the grid oversampling ratio and the kernel width, that can be used to trade accuracy for computational memory and time reductions. We have found that significant reductions in computation memory and time can be obtained while maintaining high accuracy by using a minimal oversampling ratio, from 1.125 to 1.375, instead of the typically employed grid oversampling ratio of two. When using a minimal oversampling ratio, appropriate design of the convolution kernel is important for maintaining high accuracy. We derive a simple equation for choosing the optimal Kaiser-Bessel convolution kernel for a given oversampling ratio and kernel width. As well, we evaluate the effect of presampling the kernel, a common technique used to reduce the computation time, and find that using linear interpolation between samples adds negligible error with far less samples than is necessary with nearest-neighbor interpolation. We also develop a new method for choosing the optimal presampled kernel. Using a minimal oversampling ratio and presampled kernel, we are able to perform a three-dimensional (3-D) reconstruction in one-eighth the time and requiring one-third the computer memory versus using an oversampling ratio of two and a Kaiser-Bessel convolution kernel, while maintaining the same level of accuracy.

Objective: The aim of this work was to show the potential utility of a novel rapid 3D fat-suppressed MRI method for joint imaging., Conclusion: Phase-sensitive steady-state free precession provides rapid 3D joint imaging with robust fat suppression and excellent cartilage delineation.

Purpose: To verify the potential of ungated spiral phase-contrast (USPC), which has been shown to provide accurate and reproducible time-averaged measurements of pulsatile flow, for rapid measurement of renal artery blood flow (RABF) in vivo., Materials and Methods: The RABF rates of 11 normal human subjects and one patient with renal failure were measured with USPC within six seconds., Results: Rapid USPC scans produced reproducible RABF measurements (SD < or = 9%) that agreed with the normal RABF rates known from the literature. The RABF rates of the patient with renal failure were substantially less (<50-65%) than the normal RABF rates., Conclusion: The results demonstrate that it is now possible to obtain rapid and consistent RABF measurements within six seconds with USPC.

Balanced steady-state free precession (SSFP) imaging sequences require short repetition times (TRs) to avoid off-resonance artifacts. The use of slab-selective excitations is common, as this can improve imaging speed by limiting the field of view (FOV). However, the necessarily short-duration excitations have poor slab profiles. This results in unusable slices at the slab edge due to significant flip-angle variations or aliasing in the slab direction. Variable-rate selective excitation (VERSE) is a technique by which a time-varying gradient waveform is combined with a modified RF waveform to provide the same excitation profile with different RF power and duration characteristics. With the use of VERSE, it is possible to design short-duration pulses with dramatically improved slab profiles. These pulses achieve high flip angles with only minor off-resonance sensitivity, while meeting SAR limits at 1.5 T. The improved slab profiles will enable more rapid 3D imaging of limited volumes, with more consistent image contrast across the excited slab., (Copyright 2004 Wiley-Liss, Inc.)

Refocused steady-state free precession (SSFP) is limited by its high sensitivity to local field variation, particularly at high field strengths or the long repetition times (TRs) necessary for high resolution. Several methods have been proposed to reduce SSFP banding artifact by combining multiple phase-cycled SSFP acquisitions, each differing in how individual signal magnitudes and phases are combined. These include maximum-intensity SSFP (MI-SSFP) and complex-sum SSFP (CS-SSFP). The reduction in SSFP banding is accompanied by a loss in signal-to-noise ratio (SNR) efficiency. In this work a general framework for analyzing banding artifact reduction, contrast, and SNR of any multiple-acquisition SSFP combination method is presented. A new sum-of-squares method is proposed, and a comparison is performed between each of the combination schemes. The sum-of-squares SSFP technique (SOS-SSFP) delivers both robust banding artifact reduction and higher SNR efficiency than other multiple-acquisition techniques, while preserving SSFP contrast., (Copyright 2004 Wiley-Liss, Inc.)

Magnetic resonance imaging (MRI) is limited in many cases by long scan times and low spatial resolution. Recent advances in gradient systems hardware allow very rapid imaging sequences, such as steady-state free precession (SSFP), which has repetition times (TRs) of 2-5 ms. The design of these rapid sequences demands time-optimal preparatory gradient waveforms to provide maximum readout duty-cycle, and preserve spatial resolution and SNR while keeping TRs low. Time-optimal gradient waveforms can be synthesized analytically for certain simple cases. However, certain problems, such as time-optimal 2D and 3D gradient design with moment constraints, either may not have a solution or must be solved numerically. We show that time-optimal gradient design is a convex-optimization problem, for which very efficient solution methods exist. These methods can be applied to solve gradient design problems for oblique gradient design, spiral imaging, and flow-encoding using either a constant slew rate or the more exact voltage-limited gradient models. Ultimately, these methods provide a time-optimal solution to many 2D and 3D gradient design problems in a sufficiently short time for interactive imaging., (Copyright 2003 Wiley-Liss, Inc.)

Purpose: To test the feasibility of using real-time phase contrast (PC) magnetic resonance imaging (MRI) to track velocities (1-20 cm/second) of skeletal muscle motion., Materials and Methods: To do this we modified a fast real-time spiral PC pulse sequence to accommodate through-plane velocity encoding in the range of -20 to +20 cm/second. We successfully imaged motion of the biceps brachii and triceps brachii muscles during elbow flexion and extension in seven unimpaired adult subjects using real-time PC MRI., Results: The velocity data demonstrate that the biceps brachii and the triceps brachii, antagonistic muscles, move in opposite directions during elbow flexion and extension with velocity values in the muscle tissue ranging from -10 to +10 cm/second., Conclusion: With further development, real-time PC MRI may provide a means to analyze muscle function in individuals with neurologic or movement disorders who cannot actively complete the repeated motions required for dynamic MRI techniques, such as cine PC MRI, that are more commonly used in musculoskeletal biomechanics applications., (Copyright 2003 Wiley-Liss, Inc.)

A new synthesis algorithm, based on the Shinnar-Le Roux (SLR) transform, can be used to generate fully refocused steady-state pulse sequences with arbitrary magnetization profiles as a function of off-resonant precession. This is accomplished by appropriate periodic oscillation of the RF excitation magnitude and phase from echo to echo. The technique is applied to the design of refocused steady-state free precession (SSFP) sequences with flat profiles, providing the opportunity for banding-artifact-free imaging with steady-state contrast. The algorithm is also used to generate refocused-SSFP sequences with an arbitrarily broad region of attenuated signal. These sequences are implemented and applied to the problem of steady-state fat suppression. Preliminary results show signal levels that agree well with theory, and a broad region of suppressed signal at each echo. Total imaging time is kept identical to that of a standard refocused-SSFP experiment through echo equalization and interleaving. 3D images from the leg of a normal volunteer acquired in 44 s demonstrate the applicability of the technique to fat-suppressed imaging., (Copyright 2003 Wiley-Liss, Inc.)

Flow jets containing velocities up to 5-7 m/s are common in patients with congenital defects and patients with valvular disease (stenosis and regurgitation). The quantitation of peak velocity and flow volume in these jets is clinically significant but requires specialized imaging sequences. Conventional 2DFT phase contrast sequences require lengthy acquisitions on the order of several minutes. Conventional spiral phase contrast sequences are faster, but are highly corrupted by flow artifacts at these high velocities due to phase dispersion and motion during the excitation and readout. A new prospectively gated method based on spiral phase contrast is presented, which has a sufficiently short measurement interval (<4 ms) to minimize flow artifacts, while achieving high spatial resolution (2 x 2 x 4 mm(3)) to minimize partial volume effects, all within a single breathhold. A complete single-slice phase contrast movie loop with 22 ms true temporal resolution is acquired in one 10-heartbeat breathhold. Simulations indicate that this technique is capable of imaging through-plane jets with velocities up to 10 m/s, and initial studies in aortic stenosis patients show accurate in vivo measurement of peak velocities up to 4.2 m/s (using echocardiography as a reference)., (Copyright 2003 Wiley-Liss, Inc.)

A novel ungated spiral phase-contrast (USPC) imaging method was developed for rapid measurement of time-averaged blood-flow rates in the presence of pulsatility. The spatial point-spread function was analyzed to provide an intuitive understanding of how spiral trajectories, which sample the k-space origin at every excitation, can mitigate the effects of pulsatility. Pulsatile flow phantom experiments were performed to validate the accuracy and repeatability of the USPC method. The measurement of flow in the renal and femoral arteries of normal volunteers were also performed. The phantom results (error < or = +9%, SD(phantom) < or = 2%, time-averaged pulsatile-flow rates = 3-15 ml/s) and in vivo results (SD(renal) < or = 8%, SD(femoral) < or = 14%) demonstrate the potential of the USPC method for rapidly and repeatedly measuring accurate time-averaged blood flow even in relatively small arteries and in the presence of strong pulsatility., (Copyright 2003 Wiley-Liss, Inc.)

A new method of encoding flow velocity as image phase in a refocused steady-state free precession (SSFP) sequence, called steady-state phase contrast (SSPC), can be used to generate velocity images rapidly while retaining high signal. Magnitude images with refocused-SSFP contrast are simultaneously acquired. This technique is compared with the standard method of RF-spoiled phase contrast (PC), and is found to have more than double the phase-signal to phase-noise ratio (PNR) when compared with standard PC at reasonable repetition intervals (TRs). As TR decreases, this advantage increases exponentially, facilitating rapid scans with high PNR efficiency. Rapid switching between the two necessary steady states can be accomplished by the insertion of a single TR interval with no flow-encoding gradient. The technique is implemented in a 2DFT sequence and validated in a phantom study. Preliminary results indicate that further TR reduction may be necessary for high-quality cardiac images; however, images in more stationary structures, such as the descending aorta and carotid bifurcation, exhibit good signal-to-noise ratio (SNR) and PNR. Comparisons with standard-PC images verify the PNR advantage predicted by theory., (Copyright 2002 Wiley-Liss, Inc.)

Objective: Fluctuating equilibrium magnetic resonance (FEMR) is a rapid three-dimensional (3D) imaging sequence with high signal-to-noise ratio (SNR). FEMR may be useful for detecting cartilage defects in the knee. At 1.5 T, FEMR uses a TR with odd multiples of 2.2 ms for fat/water separation. With a TR of 6.6 ms, high-resolution 3D imaging of cartilage is possible., Design and Patients: The knees of 10 volunteers and two patients were imaged on a GE Signa 1.5 T scanner using an extremity coil. Scans were preceded by a shimming sequence optimizing linear terms. Subjects were imaged with FEMR, proton-density fast spin-echo (PD-FSE), T2-weighted fast spin-echo (T2-FSE), and 3D fat-suppressed spoiled-gradient-recalled echo (3D-SPGR)., Results: SNR and contrast-to-noise efficiency measurements for cartilage using FEMR were superior to those using PD-FSE, T2-FSE, and 3D-FS-SPGR. FSE images showed bright synovial fluid with limited cartilage detail. 3D-SPGR had comparable resolution to FEMR but suboptimal cartilage/fluid contrast and longer scan times (8 min versus 2 min). Cartilage surface detail, outlined by bright synovial fluid, was best seen on the FEMR images., Discussion: FEMR obtains high-resolution 3D images of the entire knee in 2 min with excellent cartilage/fluid contrast. FEMR is sensitive to field inhomogeneity and requires shimming. Surface defects are outlined by bright synovial fluid, and cartilage has higher signal-to-noise efficiency compared with PD-FSE, T2-FSE, and 3D-SPGR techniques.

Purpose: The Magnetic resonance imaging (MRI) is an emerging technology for catheter-based imaging and interventions. Real-time MRI is a promising methodfor overcoming catheter and physiologic motion for intravascular imaging., Methods: All imaging was performed on a 1.5 T Signa MRI scanner with high-speed gradients. Multiple catheter coils were designed and constructed, including low-profile, stub-matched coils. Coil sensitivity patterns and SNR measurements were compared. Real-time imaging was performed with an interleaved spiral sequence using a dedicated workstation, providing real-time data acquisition, image reconstruction and interactive control and display. Real-time "black-blood" imaging was achieved through incorporation of off-slice saturation pulses. The imaging sequence was tested in a continuous flow phantom and then in vivo in the rabbit aorta using a 2 mm catheter coil., Results: The real-time intravascular imaging sequence achieved 120-440 micron resolution at up to 16 frames per second. Low-profile stub-tuned catheter coils achieved similar SNR to larger traditional coil designs. In the phantom experiments, addition of real-time black-blood saturation pulses effectively suppressed the flow signal and allowed visualization of the phantom wall. In vivo experiments clearly showed real-time intravascular imaging of the rabbit aortic wall with minimal motion artifacts and effective blood signal suppression., Conclusions: Real-time imaging with low-profile coil designs provides significant enhancements to intravascular MRI.