Search

Purpose Flow quantification by phase‐contrast MRI is hampered by spatially varying background phase offsets. Correction performance by polynomial regression on stationary tissue may be affected by outliers such as wrap‐around or constant flow. Therefore, we propose an alternative, M‐estimate SAmple Consensus (MSAC) to reject outliers, and improve and fully automate background phase correction. Methods The MSAC technique fits polynomials to randomly drawn small samples from the image. Over several trials, it aims to find the best consensus set of valid pixels by rejecting outliers to the fit and minimizing the residuals of the remaining pixels. The robustness of MSAC to its few parameters was investigated and verified using third‐order polynomial correction fits on a total of 118 2D flow (97 with wrap‐around) and 18 4D flow data sets (14 with wrap‐around), acquired at 1.5 T and 3 T. Background phase was compared with standard stationary correction and phantom correction. Pulmonary/systemic flow ratios in 2D flow were derived, and exemplary 4D flow analysis was performed. Results The MSAC technique is robust over a range of parameter choices, and a unique set of parameters is suitable for both 2D and 4D flow. In 2D flow, phase errors were significantly reduced by MSAC compared with stationary correction (p = 0.005), and stationary correction shows larger errors in pulmonary/systemic flow ratios compared with MSAC. In 4D flow, MSAC shows similar performance as stationary correction. Conclusions The MSAC method provides fully automated background phase correction to 2D and 4D flow data and shows improved robustness over stationary correction, especially with outliers present.

PURPOSE To introduce and investigate a method for free-breathing three-dimensional (3D) B1+ mapping of the human body at ultrahigh field (UHF), which can be used to generate homogenous flip angle (FA) distributions in the human body at UHF. METHODS A 3D relative B1+ mapping sequence with a radial phase-encoding (RPE) k-space trajectory was developed and applied in 11 healthy subjects at 7T. An RPE-based actual flip angle mapping method was applied with a dedicated B1+ shim setting to calibrate the relative B1+ maps yielding absolute B1+ maps of the individual transmit channels. The method was evaluated in a motion phantom and by multidimensional in vivo measurements. Additionally, 3D gradient echo scans with and without static phase-only B1+ shims were used to qualitatively validate B1+ shim predictions. RESULTS The phantom validation revealed good agreement for B1+ maps between dynamic measurement and static reference acquisition. The proposed 3D method was successfully validated in vivo by comparing magnitude and phase distributions with a 2D Cartesian reference. 3D B1+ maps free from visible motion artifacts were successfully acquired for 11 subjects with body mass indexes ranging from 19 kg/m2 to 34 kg/m2 . 3D respiration-resolved absolute B1+ maps indicated FA differences between inhalation and exhalation up to 15% for one channel and up to 24% for combined channels for shallow breathing. CONCLUSION The proposed method provides respiration-resolved absolute 3D B1+ maps of the human body at UHF, which enables the investigation and development of 3D B1+ shimming and parallel transmission methods to further enhance body imaging at UHF.

Purpose: Cardiovascular magnetic resonance first-pass perfusion for the pixel-wise detection of coronary artery disease is rapidly becoming the clinical standard, yet no widely available method exists for its assessment and validation. This study introduces a novel phantom capable of generating spatially dependent flow values to enable assessment of new perfusion imaging methods at the pixel level. Methods: A synthetic multicapillary myocardial phantom mimicking transmural myocardial perfusion gradients was designed and manufactured with high-precision 3D printing. The phantom was used in a stationary flow setup providing reference myocardial perfusion rates and was scanned on a 3T system. Repeated first-pass perfusion MRI for physiological perfusion rates between 1 and 4 mL/g/min was performed using a clinical dual-sequence technique. Fermi function-constrained deconvolution was used to estimate pixel-wise perfusion rate maps. Phase contrast (PC)-MRI was used to obtain velocity measurements that were converted to perfusion rates for validation of reference values and cross-method comparison. The accuracy of pixel-wise maps was assessed against simulated reference maps. Results: PC-MRI indicated excellent reproducibility in perfusion rate (coefficient of variation [CoV] 2.4-3.5%) and correlation with reference values (R2 = 0.985) across the full physiological range. Similar results were found for first-pass perfusion MRI (CoV 3.7-6.2%, R2 = 0.987). Pixel-wise maps indicated a transmural perfusion difference of 28.8-33.7% for PC-MRI and 23.8-37.7% for first-pass perfusion, matching the reference values (30.2-31.4%). Conclusion: The unique transmural perfusion pattern in the phantom allows effective pixel-wise assessment of first-pass perfusion acquisition protocols and quantification algorithms before their introduction into routine clinical use.

PURPOSE Two-dimensional selective excitation (2DRF) allows shortening 4D flow scan times by reducing the FOV, but the longer 2DRF pulse duration decreases the temporal resolution, yielding underestimated peak flow values. Multiple k-space lines per cardiac phase, nl ≥ 2, are commonly applied in 4D flow MRI to shorten the inherent long scan times. We demonstrate that 2DRF 4D flow with nl ≥ 2 can be easily combined with UNFOLD (UNaliasing by Fourier-encoding the Overlaps using the temporaL Dimension), a technique that allows regaining nominally the temporal resolution of the respective acquisition with nl = 1, to assure peak flow quantification. METHODS Two different 2DRF pulses with spiral k-space trajectories were designed and integrated into a 4D flow sequence. Flow phantom experiments and 7 healthy control 4D flow in vivo measurements, with and without UNFOLD reconstructions, were compared with conventional reconstruction and 1D slab-selective excitation (1DRF) by evaluating time-resolved flow curves, peak flow, peak velocity, blood flow volume per cardiac cycle, and spatial aliasing. RESULTS Applying UNFOLD to 4D flow imaging with 2DRF and reduced FOV increased the quantified in vivo peak flow values significantly by 3.7% ± 2.3% to 5.2% ± 2.4% (P < .05). Accordingly, the peak flow underestimation of 2DRF scans compared with conventional 1DRF scans decreased with UNFOLD. Finally, 2DRF combined with UNFOLD accelerated the 4D flow acquisition 3.5 ± 1.4 fold by reducing the FOV and increasing the effective temporal resolution by 6.7% compared with conventional 1D selective excitation, with 2 k-space lines per cardiac phase. CONCLUSION Two-dimensional selective excitation combined with UNFOLD allows limiting the FOV to shorten 4D flow scan times and compensates for the loss in temporal resolution with 2DRF (Δt = 64.8 ms) compared with 1DRF (Δt = 43.2 ms), yielding an effective resolution of Δteff = 40.5 ms to enhance peak flow quantification.

PURPOSE To improve the efficiency of native and postcontrast high-resolution cardiac T1 mapping by utilizing cardiac motion correction. METHODS Common cardiac T1 mapping techniques only acquire data in a small part of the cardiac cycle, leading to inefficient data sampling. Here, we present an approach in which 80% of each cardiac cycle is used for T1 mapping by integration of cardiac motion correction. Golden angle radial data was acquired continuously for 8 s with in-plane resolution of 1.3 × 1.3 mm2 . Cine images were reconstructed for nonrigid cardiac motion estimation. Images at different TIs were reconstructed from the same data, and motion correction was performed prior to T1 mapping. Native T1 mapping was evaluated in healthy subjects. Furthermore, the technique was applied for postcontrast T1 mapping in 5 patients with suspected fibrosis. RESULTS Cine images with high contrast were obtained, leading to robust cardiac motion estimation. Motion-corrected T1 maps showed myocardial T1 times similar to cardiac-triggered T1 maps obtained from the same data (1288 ± 49 ms and 1259 ± 55 ms, respectively) but with a 34% improved precision (spatial variation: 57.0 ± 12.5 ms and 94.8 ± 15.4 ms, respectively, P < 0.0001) due to the increased amount of data. In postcontrast T1 maps, focal fibrosis could be confirmed with late contrast-enhancement images. CONCLUSION The proposed approach provides high-resolution T1 maps within 8 s. Data acquisition efficiency for T1 mapping was improved by a factor of 5 by integration of cardiac motion correction, resulting in precise T1 maps.

PURPOSE: To provide high-resolution cardiac T(1) mapping of various cardiac phases and cine imaging within a single breath-hold using continuous golden ratio-based radial acquisition and model-based iterative image reconstruction. METHODS: Data acquisition was performed continuously using golden ratio-based radial sampling and multiple inversion pulses were applied independent of the heart rate. Native T(1) maps of diastole and systole were reconstructed with in-plane resolution of 1.3 × 1.3 mm(2) using model-based iterative image reconstruction. Cine images with 30 cardiac phases were reconstructed from the same data using kt-SENSE. The method was evaluated in a commercially available T(1) phantom and 10 healthy subjects. In vivo T(1) assessment was carried out segment-wise. RESULTS: Evaluation in the phantom demonstrated accurate T(1) times (R(2) > 0.99) and insensitivity to the heart rate. In vivo T(1) values did not differ between systole and diastole, and T(1) times assessed by the proposed approach were longer than measured with a modified Look-Locker inversion recovery (MOLLI) sequence, except for lateral segments. Cine images had a consistent dark-blood contrast and functional assessment was in agreement with assessment based on Cartesian cine scans (difference in ejection fraction: 0.26 ± 2.65%, P = 0.65). CONCLUSION: The proposed approach provides native T(1) maps of diastole and systole with high spatial resolution and cine images simultaneously within 16 s, which could strongly improve the scan efficiency.

Purpose To provide 3D multicontrast anatomical MR with high isotropic resolution and metabolic positron emission tomography (PET) images using a respiratory motion-compensated simultaneous PET-MR examination with high scan efficiency. Theory and Methods Standard abdominal PET-MR examinations combine MR data obtained during multiple breath-holds with free-breathing PET acquisitions, limiting the achievable image resolution and potentially causing misalignment errors between breath-hold and free-breathing data. Here, a 3D free-breathing PET-MR acquisition is presented, yielding T1 and T2-weighted MR images with an isotropic resolution of 1.5 mm3. In addition, nonrigid respiratory motion information and respiratory-resolved attenuation-correction maps are obtained without an increase in scan time. Motion information is used in motion-compensated image reconstructions to improve MR and PET image quality while shortening scan times. Results The proposed approach was evaluated in 11 oncology patients and provided respiratory motion information with an accuracy of 1.3 ± 0.1 mm. Sharpness of anatomical features was increased by 19 ± 13% compared with the uncorrected MR images in a 54 ± 26% shorter scan time than a gated MR acquisition. The MR-based motion information improved uptake values (75 ± 94%) and resolution (16 ± 27%) of simultaneously acquired PET images. Conclusions The proposed method provides motion-compensated 3D high-quality MR and PET images in a comprehensive and highly efficient examination. Magn Reson Med, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Purpose 4D flow MRI permits to quantify non-invasively time-dependent velocity vector fields, but it demands long acquisition times. 2D-selective excitation allows to accelerate the acquisition by reducing the FOV in both phase encoding directions. In this study, we investigate 2D-selective excitation with reduced FOV imaging for fast 4D flow imaging while obtaining correct velocity quantification. Methods Two different 2D-selective excitation pulses were designed using spiral k-space trajectories. Further, their isophase time point was analyzed using simulations that considered both stationary and moving spins. On this basis, the 2D-selective RF pulses were implemented into a 4D flow sequence. A flow phantom study and seven 4D flow in vivo measurements were performed to assess the accuracy of velocity quantification by comparing the proposed technique to non-selective and conventional 1D slab-selective excitation. Results The isophase time point for spiral 2D-selective RF pulses was found to be located at the end of excitation for both stationary and moving spins. Based on that, 2D-selective excitation with reduced FOV allowed us to successfully quantify velocities both in a flow phantom and in vivo. In a flow phantom, the velocity difference Δ v ¯ = 0.8 ± 5.3 cm/s between the smaller reduced FOV and the reference scan was similar to the inter-scan variability of Δ v ¯ = - 1.0 ± 2.3 cm/s . In vivo, the differences in flow (P = 0.995) and flow volume (P = 0.469) between the larger reduced FOV and the reference scan were non-significant. By reducing the FOV by two-thirds, acquisition time was halved. Conclusion A reduced field-of-excitation allows to limit the FOV and therefore shorten 4D flow acquisition times while preserving successful velocity quantification.

PURPOSE Myocardial blood-flow simulation based on laws of fluid mechanics is a valuable tool for understanding tissue behavior. Our aim is to evaluate the ability of a porous-media flow model approach to reflect disturbed washout of contrast media (CM) from the myocardium as observed by cardiovascular MR. METHODS A coupled advection-diffusion model is used to describe the CM flow in the vascular and extravascular space as separate compartments. Their exchange of CM is controlled by the exchange rate ExR , which in turn determines the washout behavior. We fitted simulations to CM concentration measurements, derived from T1 maps of the midventricular slice. The CM concentration was extracted from 18 patients with myocarditis in the acute phase and during follow-up after 6 months. The results were compared with 18 sex- and age-matched controls. For each subject, the measurements were acquired before and during the first 10 minutes at 5 time points after CM administration, representing CM washout. Image registration was applied to compensate for motion between different time points. RESULTS Eight matched data sets had to be excluded due to low registration quality. Processing was successful in n = 10 matched data sets of acute and healed myocarditis as well as controls. Significant differences in ExR were observed when comparing patients with acute myocarditis to controls (P < .001), to their follow-up (P < .05), and the follow-up to controls (P < .05). CONCLUSION Our study suggests the feasibility of using the proposed porous-medium flow framework for the simulation of pathologic myocardial tissue.

Purpose To improve velocity vector field reconstruction from undersampled four-dimensional (4D) flow MRI by penalizing divergence of the measured flow field. Theory and Methods Iterative image reconstruction in which magnitude and phase are regularized separately in alternating iterations was implemented. The approach allows incorporating prior knowledge of the flow field being imaged. In the present work, velocity data were regularized to reduce divergence, using either divergence-free wavelets (DFW) or a finite difference (FD) method using the l1-norm of divergence and curl. The reconstruction methods were tested on a numerical phantom and in vivo data. Results of the DFW and FD approaches were compared with data obtained with standard compressed sensing (CS) reconstruction. Results Relative to standard CS, directional errors of vector fields and divergence were reduced by 55–60% and 38–48% for three- and six-fold undersampled data with the DFW and FD methods. Velocity vector displays of the numerical phantom and in vivo data were found to be improved upon DFW or FD reconstruction. Conclusion Regularization of vector field divergence in image reconstruction from undersampled 4D flow data is a valuable approach to improve reconstruction accuracy of velocity vector fields. Magn Reson Med 75:115–125, 2016. © 2015 Wiley Periodicals, Inc.

Purpose: To develop a flexible image navigator for 3D coronary MR angiography that allows respiratory motion of variable complexity to be compensated for on different temporal scales. Methods: A two-dimensional (2D) golden radial image navigator is proposed; translational motion is compensated for on a beat-to-beat basis, and residual affine motion is compensated for on a bin-to-bin basis in a two-step procedure. The method does not use a respiratory gating window and therefore achieves 100% scan efficiency and a predictable scan time. The proposed method was tested in 11 healthy volunteers and compared against a navigator-gated and tracked acquisition. Results: The proposed method achieved comparable quantitative and qualitative image quality to a 6-mm navigator-gated and tracked scan while reducing the scan time by a factor of approximately 2. Combined motion correction using image navigators improved visualization of the coronary arteries in four of 11 subjects in comparison with translational correction only, and image quality was maintained in the remaining cases. In one case, visualization of the left anterior descending coronary artery was degraded using combined correction compared with translation correction only. Conclusions: The feasibility of correcting for 2D translational

Purpose Echo planar-based diffusion-weighted MRI (DW-MRI) requires robust suppression of fat signal. Fat suppression techniques such as inversion recovery or spectrally selective excitation with subsequent gradient spoiling can extend scan time or perform suboptimally in the presence of strong main field inhomogeneities. Chemical shift-encoded water–fat separation using iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL) is robust in areas of main field inhomogeneity but requires accurate phase information, which can be distorted by patient motion during diffusion-weighting gradients. A method is proposed to overcome this with the use of image navigators. Theory and Methods A spin echo planar imaging (SE-EPI) diffusion-weighted sequence was modified to incorporate IDEAL acquisition in combination with an image navigator to correct for patient motion-induced phase effects. Images were acquired in phantoms and in healthy volunteers in brain, pelvic, and abdominal regions. Results Without navigator, diffusion-weighted IDEAL created artifacts in areas of motion. These were removed when the two-dimensional navigator was used to correct the phase, resulting in correct water–fat separation. Conclusion DW-EPI with IDEAL and an integrated image navigator allows for robust water and fat separation in different body areas and are a time-efficient alternative to standard fat-suppression techniques in DW-MRI. Magn Reson Med 73:964–972, 2015. © 2014 Wiley Periodicals, Inc.

PURPOSE: To extend SPIRiT to additionally exploit temporal correlations for highly accelerated generalized phase-contrast MRI and to compare the performance of the proposed radial k-t SPIRiT method relative to frame-by-frame SPIRiT and radial k-t GRAPPA reconstruction for velocity and turbulence mapping in the aortic arch. THEORY AND METHODS: Free-breathing navigator-gated two-dimensional radial cine imaging with three-directional multi-point velocity encoding was implemented and fully sampled data were obtained in the aortic arch of healthy volunteers. Velocities were encoded with three different first gradient moments per axis to permit quantification of mean velocity and turbulent kinetic energy. Velocity and turbulent kinetic energy maps from up to 14-fold undersampled data were compared for k-t SPIRiT, frame-by-frame SPIRiT, and k-t GRAPPA relative to the fully sampled reference. RESULTS: Using k-t SPIRiT, improvements in magnitude and velocity reconstruction accuracy were found. Temporally resolved magnitude profiles revealed a reduction in spatial blurring with k-t SPIRiT compared with frame-by-frame SPIRiT and k-t GRAPPA for all velocity encodings, leading to improved estimates of turbulent kinetic energy. CONCLUSION: k-t SPIRiT offers improved reconstruction accuracy at high radial undersampling factors and hence facilitates the use of generalized phase-contrast MRI for routine use. Magn Reson Med, 2013. © 2013 Wiley Periodicals, Inc.

Purpose To present and validate a manifold learning (ML)-based method that estimates the respiratory signal directly from undersampled k-space data and that can be applied for respiratory self-gated liver MRI. Methods ML methods embed high-dimensional space data in a low-dimensional space while preserving their characteristic properties. These methods have been used to estimate one-dimensional respiratory motion (low-dimensional manifold) from a set of high-dimensional free-breathing abdominal MR images. These approaches require MR images to be reconstructed first from the acquired undersampled data. Recently, the concept of compressive manifold learning (CML) has been introduced that combines compressed sensing with ML by learning low-dimensional manifolds directly from a partial set of compressed measurements, provided that the sampling satisfies the restricted isometry property. We propose to use the CML concept to extract the respiratory signal directly from undersampled k-space data. Results Simulation results from free-breathing abdominal MR data show that CML can accurately estimate respiratory motion from highly retrospectively undersampled k-space (up to 25-fold acceleration under ideal assumptions). Prospective free-breathing golden-angle radial two-dimensional (2D) acquisitions further demonstrate the feasibility of the CML method for respiratory self-gating acquisition, estimating the respiratory motion from up to 15-fold accelerated MR data. Conclusion The proposed method performs accurate respiratory signal estimation from highly undersampled k-space data and can be used for respiratory self-navigated 2D liver MRI. Magn Reson Med 72:1130–1140, 2014. © 2013 Wiley Periodicals, Inc.

Purpose Thoracic pulse wave velocity (PWV) variation due to modulated trans-mural pressure (TMP) may indicate mechanical properties of the aorta. Our aim was to measure beat-to-beat thoracic PWV and TMP to observe its normal variation during respiratory maneuvers. Methods We validated PWV measurements from a real-time velocity projection MRI scan in a pulsatile phantom. A volunteer study showed inter-scan repeatability of steady-state PWV, and observed PWV variation when performing Mueller and Valsalva maneuvers. Synchronized to the real-time projection velocity data, TMP was measured using a mouth piece and pressure sensor arrangement monitoring the intra-thoracic pressure and a single arterial pressure measurement. Results In the phantom, beat-to-beat PWV derived from real-time projection (5.33 ± 0.32 m s−1) agreed well with experimentally derived PWV using ultrasound probes (5.72 ± 0.50 m s−1). The within-subject PWV variation between scans was 0.28 m s−1. Volunteers' PWVs increased during Mueller maneuver (TMP increase of 14.67 ± 10.69 mmHg) by 32% (P

To provide 3D multicontrast anatomical MR with high isotropic resolution and metabolic positron emission tomography (PET) images using a respiratory motion-compensated simultaneous PET-MR examination with high scan efficiency.Standard abdominal PET-MR examinations combine MR data obtained during multiple breath-holds with free-breathing PET acquisitions, limiting the achievable image resolution and potentially causing misalignment errors between breath-hold and free-breathing data. Here, a 3D free-breathing PET-MR acquisition is presented, yielding TThe proposed approach was evaluated in 11 oncology patients and provided respiratory motion information with an accuracy of 1.3 ± 0.1 mm. Sharpness of anatomical features was increased by 19 ± 13% compared with the uncorrected MR images in a 54 ± 26% shorter scan time than a gated MR acquisition. The MR-based motion information improved uptake values (75 ± 94%) and resolution (16 ± 27%) of simultaneously acquired PET images.The proposed method provides motion-compensated 3D high-quality MR and PET images in a comprehensive and highly efficient examination. Magn Reson Med 79:900-911, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Compressed sensing (CS) has been demonstrated to accelerate MRI acquisitions by reconstructing sparse images of good quality from highly undersampled data. Motion during MR scans can cause inconsistencies in k-space data, resulting in strong motion artifacts in the reconstructed images. For CS to be useful in these applications, motion correction techniques need to be combined with the undersampled reconstruction. Recently, joint motion correction and CS approaches have been proposed to partially correct for effects of motion. However, the main limitation of these approaches is that they can only correct for affine deformations. In this work, we propose a novel motion corrected CS framework for free-breathing dynamic cardiac MRI that incorporates a general motion correction formulation directly into the CS reconstruction. This framework can correct for arbitrary affine or nonrigid motion in the CS reconstructed cardiac images, while simultaneously benefiting from highly accelerated MR acquisition. The application of this approach is demonstrated both in simulations and in vivo data for 2D respiratory self-gated free-breathing cardiac CINE MRI, using a golden angle radial acquisition. Results show that this approach allows for the reconstruction of respiratory motion corrected cardiac CINE images with similar quality to breath-held acquisitions.

Compressed sensing has been of great interest to speed up the acquisition of MR images. The k-t group sparse (k-t GS) method has recently been introduced for dynamic MR images to exploit not just the sparsity, as in compressed sensing, but also the spatial group structure in the sparse representation. k-t GS achieves higher acceleration factors compared to the conventional compressed sensing method. However, it assumes a spatial structure in the sparse representation and it requires a time consuming hard-thresholding reconstruction scheme. In this work, we propose to modify k-t GS by incorporating prior information about the sorted intensity of the signal in the sparse representation, for a more general and robust group assignment. This approach is referred to as group sparse reconstruction using intensity-based clustering. The feasibility of the proposed method is demonstrated for static 3D hyperpolarized lung images and applications with both dynamic and intensity changes, such as 2D cine and perfusion cardiac MRI, with retrospective undersampling. For all reported acceleration factors the proposed method outperforms the original compressed sensing method. Improved reconstruction over k-t GS method is demonstrated when k-t GS assumptions are not satisfied. The proposed method was also applied to cardiac cine images with a prospective sevenfold acceleration, outperforming the standard compressed sensing reconstruction.

The applicability of cine blood flow measurements in a clinical setting is often compromised by the long scan times associated with phase-contrast imaging. In this work, we propose an extension to the k-t principal component analysis method and demonstrate that by definition of spatial compartment-dependent temporal basis functions, significant improvements in reconstruction accuracy can be achieved relative to the original k-t principal component analysis and k-t SENSE formulations. Using this method, it is shown that prospective nominal undersampling of up to 16 corresponding to a net acceleration factor of 8 including training data acquisition can be realized while keeping the error in stroke volume below 5%. As a practical application, the acquisition of cine flow data in the aorta is demonstrated permitting assessment of two-dimensional velocity images and pulse wave velocities at 100 frames per second in a single breathhold per slice.

Cardiovascular diseases, including arrhythmias and heart failure, are commonly treated with percutaneous procedures guided by X-ray fluoroscopy. The visualization of the targeted structures can be enhanced using preacquired respiratory-resolved anatomic data (dynamic roadmap), which is displayed as an overlay onto X-ray fluoroscopy images. This article demonstrates how dynamic roadmaps using an affine motion model can be obtained from one respiratory-resolved three-dimensional whole-heart acquisition using the previously introduced Radial Phase Encoding-Phase Ordering with Automatic Window Selection method. Respiratory motion in different regions of the heart was analyzed in 10 volunteers, and it was shown that the use of dynamic roadmaps can reduce misalignment errors from more than 10 down to less than 1.5 mm. Furthermore, the results suggest that reliable motion information can be obtained from highly undersampled images due to the advantageous undersampling properties of the radial phase encoding trajectory. Finally, results of a three-dimensional dynamic roadmap obtained from a patient before catheter ablation for atrial fibrillation treatment are presented.

The value of phase-contrast magnetic resonance imaging for quantifying tissue motion and blood flow has been long recognized. However, the sensitivity of the method to system imperfections can lead to inaccuracies limiting its clinical acceptance. A key source of error relates to eddy current-induced phase fluctuations, which can offset the measured object velocity significantly. A higher-order dynamic field camera was used to study the spatiotemporal evolution of background phases in cine phase-contrast measurements. It is demonstrated that eddy current-induced offsets in phase-difference data are present up to the second spatial order. Oscillatory temporal behaviors of offsets in the kHz range suggest mechanical resonances of the MR system to be non-negligible in phase-contrast imaging. By careful selection of the echo time, their impact can be significantly reduced. When applying field monitoring data for correcting eddy current and mechanically induced velocity offsets, errors decrease to less than 0.5% of the maximum velocity for various sequence settings proving the robustness of the correction approach. In vivo feasibility is demonstrated for aortic and pulmonary flow measurements in five healthy subjects. Using field monitoring data, mean error in stroke volume was reduced from 10% to below 3%. Magn Reson Med, 2012. © 2011 Wiley Periodicals, Inc.

Three dimensional (3D) whole-heart magnetic resonance imaging (MRI) has become an important imaging modality to assess cardiovascular diseases. The main challenges for 3D whole-heart MRI are long acquisition times, required to achieve high spatial resolution, and image artefacts due to physiological motion. Here we propose to overcome these problems by the combination of an interleaved Radial Phase Encoding trajectory and the Phase Ordering with Automatic Window Selection method. This Radial Phase Encoding-Phase Ordering with Automatic Window Selection approach yields fast 3D whole-heart imaging with a high isotropic resolution and high navigator efficiency even for extremely irregular breathing. Numerical simulations were performed and Radial Phase Encoding-Phase Ordering with Automatic Window Selection was implemented on a clinical scanner. A comparison between the proposed method and a respiratory gated 3D Cartesian approach was carried out. Radial Phase Encoding-Phase Ordering with Automatic Window Selection leads to a better depiction of coronary arteries and an increase in navigator efficiency. In addition to a high resolution image, this method also provides dynamic respiratory information without an increase in scan time. Magn Reson Med, 2011. © 2011 Wiley-Liss, Inc.

Compressed sensing (CS) is a data-reduction technique that has been applied to speed up the acquisition in MRI. However, the use of this technique in dynamic MR applications has been limited in terms of the maximum achievable reduction factor. In general, noise-like artefacts and bad temporal fidelity are visible in standard CS MRI reconstructions when high reduction factors are used. To increase the maximum achievable reduction factor, additional or prior information can be incorporated in the CS reconstruction. Here, a novel CS reconstruction method is proposed that exploits the structure within the sparse representation of a signal by enforcing the support components to be in the form of groups. These groups act like a constraint in the reconstruction. The information about the support region can be easily obtained from training data in dynamic MRI acquisitions. The proposed approach was tested in two-dimensional cardiac cine MRI with both downsampled and undersampled data. Results show that higher acceleration factors (up to 9-fold), with improved spatial and temporal quality, can be obtained with the proposed approach in comparison to the standard CS reconstructions. Magn Reson Med, 2011. © 2011 Wiley-Liss, Inc.

Fluorine MRI offers broad potential for specific detection and quantification of molecularly targeted agents in diagnosis and therapy planning or monitoring. Because non-proton MRI applications lack morphological information, accompanying proton images are needed to elucidate the spatial tissue context. Furthermore, low concentrations typical of targeted molecular imaging agents require long examinations for signal averaging during which physiological motion may lead to blurring, underestimation in signal quantification, and erroneous localization of the agent distribution. Novel methods for truly simultaneous acquisition of dual-nuclei MR data are presented that offer efficient and precise anatomical localization of fluorine signals using accurate motion correction based on contemporaneous proton signals. The feasibility of simultaneous dual-nuclei MRI motion correction and corresponding dual-resolution reconstruction, providing nuclei-specific spatial resolution to retrospectively optimize the balance between signal-to-noise ratio and resolution, is shown on a clinical 3 T MR system.

In magnetic resonance imaging, implantable devices are usually visualized with a negative contrast. Recently, positive contrast techniques have been proposed, such as susceptibility gradient mapping (SGM). However, SGM reduces the spatial resolution making positive visualization of small structures difficult. Here, a development of SGM using the original resolution (SUMO) is presented. For this, a filter is applied in k-space and the signal amplitude is analyzed in the image domain to determine quantitatively the susceptibility gradient for each pixel. It is shown in simulations and experiments that SUMO results in a better visualization of small structures in comparison to SGM. SUMO is applied to patient datasets for visualization of stent and prostate brachytherapy seeds. In addition, SUMO also provides quantitative information about the number of prostate brachytherapy seeds. The method might be extended to application for visualization of other interventional devices, and, like SGM, it might also be used to visualize magnetically labelled cells.

An MR-electrophysiology (EP) catheter is presented that provides full diagnostic EP functionality and a high level of radiofrequency safety achieved by custom-designed transmission lines. Highly resistive wires transmit intracardiac electrograms and currents for intracardiac pacing. A transformer cable transmits the localization signal of a tip coil. Specific absorption rate simulations and temperature measurements at 1.5 T demonstrate that a wire resistance > 3 kΩ/m limits dielectric heating to a physiologically irrelevant level. Additional wires do not increase tip specific absorption rate significantly, which is important because some clinical catheters require up to 20 electrodes. It is further demonstrated that radiofrequency-induced and pacing-induced resistive heating of the wires is negligible under clinical conditions. The MR-EP catheters provided uncompromised recording of electrograms and cardiac pacing in combination with a standard EP recorder in MR-guided in vivo EP studies, and the tip coil enabled fast and robust catheter localization. In vivo temperature measurements during such a study did not detect any device-related heating, which confirms the high level of safety of the catheter, whereas unacceptable heating was found with a standard EP catheter. The presented concept for the first time enables catheters with full diagnostic EP functionality and active tracking and at the same time a sufficient level of radiofrequency safety for MRI without specific absorption rate-related limitations.

Whole-heart isotropic nonangulated cardiac magnetic resonance (CMR) is becoming an important protocol in simplifying MRI, since it reduces the need of cumbersome planning of angulations. However the acquisition times of whole-heart MRI are prohibitive due to the large fields of view (FOVs) and the high spatial resolution required for depicting small structures and vessels. To address this problem, we propose a three-dimensional (3D) acquisition scheme that combines Cartesian sampling in the readout direction with an undersampled radial scheme in the phase-encoding plane. Different undersampling patterns were investigated in combination with an iterative sensitivity encoding (SENSE) reconstruction and a 32-channel cardiac coil. Noise amplification maps were calculated to compare the performance of the different patterns using iterative SENSE reconstruction. The radial phase-encoding (RPE) scheme was implemented on a clinical MR scanner and tested on phantoms and healthy volunteers. The proposed method exhibits better image quality even for high acceleration factors (up to 12) in comparison to Cartesian acquisitions.

Four-dimensional (4D) flow imaging has been used to study flow patterns and pathophysiology, usually focused on specific thoracic vessels and cardiac chambers. Whole-heart 4D flow at high measurement accuracy covering the entire thoracic cardiovascular system would be desirable to simplify and improve hemodynamic assessment. This has been a challenge because compensation of respiratory motion is difficult to achieve, but it is paramount to limit artifacts and improve accuracy. In this work we propose a self-gating technique for respiratory motion-compensation integrated into a whole-heart 4D flow acquisition that overcomes these challenges. Flow components are measured in all three directions for each pixel over the complete cardiac cycle, and 1D volume projections are obtained at certain time intervals for respiratory gating in real time during the acquisition. The technique was tested in 15 volunteers, in which stroke volumes (SVs) in the great arteries showed excellent agreement with standard 2D phase-contrast (PC) scans. In contrast, nongated 4D flow with two averages had substantial disagreement with 2D flow. Applied to patients with congenital cardiac left-to-right shunting, the precision of flux data was highly beneficial. The methodology presented here has the potential to allow a complete study of flow pathophysiology of the thoracic cardiovascular system from a single free-breathing scan.

A crucial requirement in MR-guided interventions is the visualization of catheter devices in real time. However, true 3D visualization of the full length of catheters has hitherto been impossible given scan time constraints. Compressed sensing (CS) has recently been proposed as a method to accelerate MR imaging of sparse objects. Images acquired with active interventional devices exhibit a high CNR and are inherently sparse, therefore rendering CS ideally suited for accelerating data acquisition. A framework for true visualization of active catheters in 3D is proposed employing CS to gain high undersampling factors making real-time applications feasible. Constraints are introduced taking into account prior knowledge of catheter geometry and catheter motion over time to improve and accelerate image reconstruction. The potential of the method is demonstrated using computer simulations and phantom experiments and in vivo feasibility is demonstrated in a pig experiment. Magn Reson Med, 2009. © 2009 Wiley-Liss, Inc.

Delayed contrast-enhanced magnetic resonance imaging (DCE-MRI) provides prognostic information by delineating regions of myocardial scar. The mechanism of this delayed enhancement in myocardial infarctions (MIs) is hypothesized to result from altered kinetics and changes in the volumes of distribution in the myocardium. Pharmacokinetic models with two and three compartments were fitted to the concentration-time curves of dynamic contrast-enhanced MRI data obtained from five patients with known MI. Furthermore, the parameter stability was investigated in simulations for the two different models. The transfer constants and volumes of distribution showed a good correlation with imaging findings on early and delayed contrast-enhanced MRI. The two compartment model showed higher parameter stability. The three compartment model allows a more in-depth quantification of myocardial scarring. These models have the potential to improve the diagnosis of myocardial pathologies involving scar, with differing kinetics and volumes of distribution such as infarction or cardiomyopathy.

Local susceptibility gradients result in a dephasing of the precessing magnetic moments and thus in a fast decay of the NMR signals. In particular, cells labeled with superparamagnetic iron oxide particles (SPIOs) induce hypointensities, making the in vivo detection of labeled cells from such a negative image contrast difficult. In this work, a new method is proposed to selectively turn this negative contrast into a positive contrast. The proposed method calculates the susceptibility gradient and visualizes it in a parametric map directly from a regular gradient-echo image dataset. The susceptibility gradient map is determined in a postprocessing step, requiring no dedicated pulse sequences or adaptation of the sequence before and during image acquisition. Phantom experiments demonstrated that local susceptibility differences can be quantified. In vivo experiments showed the feasibility of the method for tracking of SPIO-labeled cells. The method bears the potential also for usage in other applications, including the detection of contrast agents and interventional devices as well as metal implants.

T(2)* relaxometry for quantitative MR imaging is strongly hampered by large-scale field inhomogeneities, which lead to signal losses and an overestimation of the relaxation rate R(2)*. This is of particular importance for the sensitive detection of iron oxide contrast agent distributions. To derive an accurate measurement of T(2)*, a main field inhomogeneity correction is applied: the main field inhomogeneity is derived from multislice T(2)* relaxometry data and used as an initial value for an iterative optimization, by which the relaxation signal is corrected for each voxel. These corrected T(2)* maps show reduced influence of the local field variation and contain information about the local SPIO concentration. The method was tested on phantoms and the limit of detection of SPIO labeled cells using T(2)* relaxometry was estimated in volunteers to be 120 x 10(3) cells/mL (2.4 microg Fe/mL) in the brain and 385 x 10(3) cells/mL (8 microg Fe/mL) in the liver.

This work presents an investigation into catheter visualization and localization using 19F nuclear magnetic resonance (NMR) in conjunction with proton imaging. For this purpose, the imaging capabilities of a standard system were extended to allow for 19F excitation and signal detection. Two modes of operation were implemented: 1) a real-time tracking mode that provides tip tracking and automatic slice position updates interleaved with real-time, interactive proton imaging; and 2) a non-real-time catheter length visualization mode in which the entire length of a catheter can be assessed. Initial phantom experiments were conducted with the use of an angiographic balloon catheter filled with the blood substitute perfluorooctylbromide (PFOB). Using limited bandwidth excitation centered at the resonances of the CF2 groups of PFOB, we found that sufficient signal could be received to facilitate tip tracking during catheter motion and length visualization for various catheter configurations. The present approach is considered a promising alternative to existing methods, which either are associated with safety concerns (if active markers are employed) or suffer from insufficient, direction-dependent contrast (if passive visualization is used). Furthermore, our approach enables visualization of the entire length of the catheter. The proposed method provides a safe technique that, unlike electrical or optical devices, does not require modification of commercially available catheters.

The purpose of this study was to test the in vivo feasibility of safe automatic catheter tracking based on an optically detunable resonant marker installed on the catheter tip, and also to test the compatibility of this approach with guidewire materials. The design of the resonant marker and the integration into the real-time MR environment is described. The catheter was used for real-time MR-guided catheterization of the aorta, left ventricle, and carotid in two swine. For in-plane visualization, the marker was repeatedly detuned. For automatic slice tracking, a projection difference measurement including detuning was interleaved with the imaging sequence. In vitro experiments were conducted to investigate the RF-safety of the marker and the effect of the guidewires on the signal intensity. For all orientations the marker provided excellent in vivo contrast using a radial steady-state free-precession sequence. Flashing of the marker by repetitive tuning/detuning further improved the in-plane visualization. Automatic slice tracking during real-time imaging was successfully performed. The plastic guidewires did not interfere with the marker, and detuning by guidewires containing nitinol could be compensated. In conclusion, automatic slice tracking as well as excellent in-plane visualization can be achieved with this approach and it is safe with respect to RF transmission.

A new method for the selective spin labeling of left- or right-sided supplying arteries of the brain without the need for additional RF coils is demonstrated. A clinical 1.5 T scanner was used. The spatial selectivity of the labeling process is based on the limited coverage of the excitation field of a standard send/receive head coil together with an oblique positioning of the labeling plane. A computer simulation was used to optimize key labeling parameters under the condition of laminar flow. The validity of the computer model results was confirmed by MRI measurements with a flow model. For human studies, a double-inversion continuous arterial spin labeling (CASL) sequence was modified to allow for arbitrary positioning of the labeling plane. The obtained perfusion-weighted images showed a clear delineation of the perfusion territories of the selected arteries in the anterior circulation of the brain and good gray/white matter contrast. Magn Reson Med 52:1443–1447, 2004. © 2004 Wiley-Liss, Inc.

To improve velocity vector field reconstruction from undersampled four-dimensional (4D) flow MRI by penalizing divergence of the measured flow field.Iterative image reconstruction in which magnitude and phase are regularized separately in alternating iterations was implemented. The approach allows incorporating prior knowledge of the flow field being imaged. In the present work, velocity data were regularized to reduce divergence, using either divergence-free wavelets (DFW) or a finite difference (FD) method using the ℓ1-norm of divergence and curl. The reconstruction methods were tested on a numerical phantom and in vivo data. Results of the DFW and FD approaches were compared with data obtained with standard compressed sensing (CS) reconstruction.Relative to standard CS, directional errors of vector fields and divergence were reduced by 55-60% and 38-48% for three- and six-fold undersampled data with the DFW and FD methods. Velocity vector displays of the numerical phantom and in vivo data were found to be improved upon DFW or FD reconstruction.Regularization of vector field divergence in image reconstruction from undersampled 4D flow data is a valuable approach to improve reconstruction accuracy of velocity vector fields.

To develop a flexible image navigator for 3D coronary MR angiography that allows respiratory motion of variable complexity to be compensated for on different temporal scales.A two-dimensional (2D) golden radial image navigator is proposed; translational motion is compensated for on a beat-to-beat basis, and residual affine motion is compensated for on a bin-to-bin basis in a two-step procedure. The method does not use a respiratory gating window and therefore achieves 100% scan efficiency and a predictable scan time. The proposed method was tested in 11 healthy volunteers and compared against a navigator-gated and tracked acquisition.The proposed method achieved comparable quantitative and qualitative image quality to a 6-mm navigator-gated and tracked scan while reducing the scan time by a factor of approximately 2. Combined motion correction using image navigators improved visualization of the coronary arteries in four of 11 subjects in comparison with translational correction only, and image quality was maintained in the remaining cases. In one case, visualization of the left anterior descending coronary artery was degraded using combined correction compared with translation correction only.The feasibility of correcting for 2D translational motion on a beat-to-beat basis as well as affine motion on a bin-to-bin basis has been demonstrated.

A balanced fast field echo (FFE) sequence (also referred to as true fast imaging with steady precession (true FISP)), based on projection reconstruction (PR) is evaluated in combination with real-time reconstruction and interactive scanning capabilities for cardiac function studies. Cardiac image sequences obtained with the balanced PR-FFE method are compared with images obtained with a spin-warp (2D Fourier transform (2DFT)) technique. In particular, the representation of motion artifacts in both techniques is investigated. Balanced PR-FFE provides a similar contrast to spin-warp-related techniques, but is less sensitive to motion artifacts. The use of angular undersampling within balanced PR-FFE is examined as a means to increase temporal resolution while causing only minor artifacts. Furthermore, a modification of the profile order allows the reconstruction of PR images at different spatial and temporal resolution levels from the same data. This study shows that balanced PR-FFE is a robust tool for cardiac function studies.

Echo planar-based diffusion-weighted MRI (DW-MRI) requires robust suppression of fat signal. Fat suppression techniques such as inversion recovery or spectrally selective excitation with subsequent gradient spoiling can extend scan time or perform suboptimally in the presence of strong main field inhomogeneities. Chemical shift-encoded water-fat separation using iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL) is robust in areas of main field inhomogeneity but requires accurate phase information, which can be distorted by patient motion during diffusion-weighting gradients. A method is proposed to overcome this with the use of image navigators.A spin echo planar imaging (SE-EPI) diffusion-weighted sequence was modified to incorporate IDEAL acquisition in combination with an image navigator to correct for patient motion-induced phase effects. Images were acquired in phantoms and in healthy volunteers in brain, pelvic, and abdominal regions.Without navigator, diffusion-weighted IDEAL created artifacts in areas of motion. These were removed when the two-dimensional navigator was used to correct the phase, resulting in correct water-fat separation.DW-EPI with IDEAL and an integrated image navigator allows for robust water and fat separation in different body areas and are a time-efficient alternative to standard fat-suppression techniques in DW-MRI.

To extend SPIRiT to additionally exploit temporal correlations for highly accelerated generalized phase-contrast MRI and to compare the performance of the proposed radial k-t SPIRiT method relative to frame-by-frame SPIRiT and radial k-t GRAPPA reconstruction for velocity and turbulence mapping in the aortic arch.Free-breathing navigator-gated two-dimensional radial cine imaging with three-directional multi-point velocity encoding was implemented and fully sampled data were obtained in the aortic arch of healthy volunteers. Velocities were encoded with three different first gradient moments per axis to permit quantification of mean velocity and turbulent kinetic energy. Velocity and turbulent kinetic energy maps from up to 14-fold undersampled data were compared for k-t SPIRiT, frame-by-frame SPIRiT, and k-t GRAPPA relative to the fully sampled reference.Using k-t SPIRiT, improvements in magnitude and velocity reconstruction accuracy were found. Temporally resolved magnitude profiles revealed a reduction in spatial blurring with k-t SPIRiT compared with frame-by-frame SPIRiT and k-t GRAPPA for all velocity encodings, leading to improved estimates of turbulent kinetic energy.k-t SPIRiT offers improved reconstruction accuracy at high radial undersampling factors and hence facilitates the use of generalized phase-contrast MRI for routine use.

To present and validate a manifold learning (ML)-based method that estimates the respiratory signal directly from undersampled k-space data and that can be applied for respiratory self-gated liver MRI.ML methods embed high-dimensional space data in a low-dimensional space while preserving their characteristic properties. These methods have been used to estimate one-dimensional respiratory motion (low-dimensional manifold) from a set of high-dimensional free-breathing abdominal MR images. These approaches require MR images to be reconstructed first from the acquired undersampled data. Recently, the concept of compressive manifold learning (CML) has been introduced that combines compressed sensing with ML by learning low-dimensional manifolds directly from a partial set of compressed measurements, provided that the sampling satisfies the restricted isometry property. We propose to use the CML concept to extract the respiratory signal directly from undersampled k-space data.Simulation results from free-breathing abdominal MR data show that CML can accurately estimate respiratory motion from highly retrospectively undersampled k-space (up to 25-fold acceleration under ideal assumptions). Prospective free-breathing golden-angle radial two-dimensional (2D) acquisitions further demonstrate the feasibility of the CML method for respiratory self-gating acquisition, estimating the respiratory motion from up to 15-fold accelerated MR data.The proposed method performs accurate respiratory signal estimation from highly undersampled k-space data and can be used for respiratory self-navigated 2D liver MRI.

Purpose Electrocardiogram (ECG)-gated cine MRI provides highly accurate functional assessment of the heart. Nevertheless, reliable ECG signals are not always available due to patient's electrophysiology or due to high MR field strengths. Here, a novel framework for cardiac functional assessment using physiological information is presented, which is obtained from MR image data. Methods Multiple long-axis slices rotated around the center axis of the left ventricle are acquired using a 2D Golden Radial acquisition scheme. This sampling approach allows for both real-time data and retrospectively reordered cine images with different temporal resolutions. Functional information from the left ventricle is used for retrospective reordering of the data to reconstruct cine images without an external ECG signal. Afterward, individual 2D cine slices are synchronized using physiological information on the mitral valve closure. The proposed approach was assessed in 15 volunteers and applied in four patients for feasibility. Results Physiological gating signals obtained with our approach show great correlation with an ECG reference signal. Functional parameters determined with the presented method show a relative difference of less than 1.3% when compared with an ECG-gated approach. Conclusion It is successfully demonstrated that functional assessment of the heart is possible using physiological information obtained directly from MR image data. Magn Reson Med 71:942–954, 2014. © 2013 Wiley Periodicals, Inc.

Thoracic pulse wave velocity (PWV) variation due to modulated trans-mural pressure (TMP) may indicate mechanical properties of the aorta. Our aim was to measure beat-to-beat thoracic PWV and TMP to observe its normal variation during respiratory maneuvers.We validated PWV measurements from a real-time velocity projection MRI scan in a pulsatile phantom. A volunteer study showed inter-scan repeatability of steady-state PWV, and observed PWV variation when performing Mueller and Valsalva maneuvers. Synchronized to the real-time projection velocity data, TMP was measured using a mouth piece and pressure sensor arrangement monitoring the intra-thoracic pressure and a single arterial pressure measurement.In the phantom, beat-to-beat PWV derived from real-time projection (5.33 ± 0.32 m s(-1) ) agreed well with experimentally derived PWV using ultrasound probes (5.72 ± 0.50 m s(-1) ). The within-subject PWV variation between scans was 0.28 m s(-1) . Volunteers' PWVs increased during Mueller maneuver (TMP increase of 14.67 ± 10.69 mmHg) by 32% (P 0.001), and during Valsalva maneuver (TMP decrease of TMP = 17.01 ± 12.91 mmHg), PWV response were inconsistent with an average increase of 14% (P 0.05).Gating TMP to beat-to-beat PWV allows insight into how aortic stiffness varies with strain. However, quantifying nonlinear arterial stiffness requires real-time arterial pressure measurement.

Compressed sensing (CS) has been demonstrated to accelerate MRI acquisitions by reconstructing sparse images of good quality from highly undersampled data. Motion during MR scans can cause inconsistencies in k-space data, resulting in strong motion artifacts in the reconstructed images. For CS to be useful in these applications, motion correction techniques need to be combined with the undersampled reconstruction. Recently, joint motion correction and CS approaches have been proposed to partially correct for effects of motion. However, the main limitation of these approaches is that they can only correct for affine deformations. In this work, we propose a novel motion corrected CS framework for free-breathing dynamic cardiac MRI that incorporates a general motion correction formulation directly into the CS reconstruction. This framework can correct for arbitrary affine or nonrigid motion in the CS reconstructed cardiac images, while simultaneously benefiting from highly accelerated MR acquisition. The application of this approach is demonstrated both in simulations and in vivo data for 2D respiratory self-gated free-breathing cardiac CINE MRI, using a golden angle radial acquisition. Results show that this approach allows for the reconstruction of respiratory motion corrected cardiac CINE images with similar quality to breath-held acquisitions.

The value of phase-contrast magnetic resonance imaging for quantifying tissue motion and blood flow has been long recognized. However, the sensitivity of the method to system imperfections can lead to inaccuracies limiting its clinical acceptance. A key source of error relates to eddy current-induced phase fluctuations, which can offset the measured object velocity significantly. A higher-order dynamic field camera was used to study the spatiotemporal evolution of background phases in cine phase-contrast measurements. It is demonstrated that eddy current-induced offsets in phase-difference data are present up to the second spatial order. Oscillatory temporal behaviors of offsets in the kHz range suggest mechanical resonances of the MR system to be non-negligible in phase-contrast imaging. By careful selection of the echo time, their impact can be significantly reduced. When applying field monitoring data for correcting eddy current and mechanically induced velocity offsets, errors decrease to less than 0.5% of the maximum velocity for various sequence settings proving the robustness of the correction approach. In vivo feasibility is demonstrated for aortic and pulmonary flow measurements in five healthy subjects. Using field monitoring data, mean error in stroke volume was reduced from 10% to below 3%.

Fluorine MRI offers broad potential for specific detection and quantification of molecularly targeted agents in diagnosis and therapy planning or monitoring. Because non-proton MRI applications lack morphological information, accompanying proton images are needed to elucidate the spatial tissue context. Furthermore, low concentrations typical of targeted molecular imaging agents require long examinations for signal averaging during which physiological motion may lead to blurring, underestimation in signal quantification, and erroneous localization of the agent distribution. Novel methods for truly simultaneous acquisition of dual-nuclei MR data are presented that offer efficient and precise anatomical localization of fluorine signals using accurate motion correction based on contemporaneous proton signals. The feasibility of simultaneous dual-nuclei MRI motion correction and corresponding dual-resolution reconstruction, providing nuclei-specific spatial resolution to retrospectively optimize the balance between signal-to-noise ratio and resolution, is shown on a clinical 3 T MR system.

Three dimensional (3D) whole-heart magnetic resonance imaging (MRI) has become an important imaging modality to assess cardiovascular diseases. The main challenges for 3D whole-heart MRI are long acquisition times, required to achieve high spatial resolution, and image artefacts due to physiological motion. Here we propose to overcome these problems by the combination of an interleaved Radial Phase Encoding trajectory and the Phase Ordering with Automatic Window Selection method. This Radial Phase Encoding-Phase Ordering with Automatic Window Selection approach yields fast 3D whole-heart imaging with a high isotropic resolution and high navigator efficiency even for extremely irregular breathing. Numerical simulations were performed and Radial Phase Encoding-Phase Ordering with Automatic Window Selection was implemented on a clinical scanner. A comparison between the proposed method and a respiratory gated 3D Cartesian approach was carried out. Radial Phase Encoding-Phase Ordering with Automatic Window Selection leads to a better depiction of coronary arteries and an increase in navigator efficiency. In addition to a high resolution image, this method also provides dynamic respiratory information without an increase in scan time.

One of the current limitations of dynamic contrast-enhanced MR angiography is the requirement of both high spatial and high temporal resolution. Several undersampling techniques have been proposed to overcome this problem. However, in most of these methods the tradeoff between spatial and temporal resolution is constant for all the time frames and needs to be specified prior to data collection. This is not optimal for dynamic contrast-enhanced MR angiography where the dynamics of the process are difficult to predict and the image quality requirements are changing during the bolus passage. Here, we propose a new highly undersampled approach that allows the retrospective adaptation of the spatial and temporal resolution. The method combines a three-dimensional radial phase encoding trajectory with the golden angle profile order and non-Cartesian Sensitivity Encoding (SENSE) reconstruction. Different regularization images, obtained from the same acquired data, are used to stabilize the non-Cartesian SENSE reconstruction for the different phases of the bolus passage. The feasibility of the proposed method was demonstrated on a numerical phantom and in three-dimensional intracranial dynamic contrast-enhanced MR angiography of healthy volunteers. The acquired data were reconstructed retrospectively with temporal resolutions from 1.2 sec to 8.1 sec, providing a good depiction of small vessels, as well as distinction of different temporal phases. Magn Reson Med 64:514-526, 2010. (C) 2010 Wiley-Liss, Inc.

An MR-electrophysiology (EP) catheter is presented that provides full diagnostic EP functionality and a high level of radiofrequency safety achieved by custom-designed transmission lines. Highly resistive wires transmit intracardiac electrograms and currents for intracardiac pacing. A transformer cable transmits the localization signal of a tip coil. Specific absorption rate simulations and temperature measurements at 1.5 T demonstrate that a wire resistance > 3 kΩ/m limits dielectric heating to a physiologically irrelevant level. Additional wires do not increase tip specific absorption rate significantly, which is important because some clinical catheters require up to 20 electrodes. It is further demonstrated that radiofrequency-induced and pacing-induced resistive heating of the wires is negligible under clinical conditions. The MR-EP catheters provided uncompromised recording of electrograms and cardiac pacing in combination with a standard EP recorder in MR-guided in vivo EP studies, and the tip coil enabled fast and robust catheter localization. In vivo temperature measurements during such a study did not detect any device-related heating, which confirms the high level of safety of the catheter, whereas unacceptable heating was found with a standard EP catheter. The presented concept for the first time enables catheters with full diagnostic EP functionality and active tracking and at the same time a sufficient level of radiofrequency safety for MRI without specific absorption rate-related limitations.

In magnetic resonance imaging, implantable devices are usually visualized with a negative contrast. Recently, positive contrast techniques have been proposed, such as susceptibility gradient mapping (SGM). However, SGM reduces the spatial resolution making positive visualization of small structures difficult. Here, a development of SGM using the original resolution (SUMO) is presented. For this, a filter is applied in k-space and the signal amplitude is analyzed in the image domain to determine quantitatively the susceptibility gradient for each pixel. It is shown in simulations and experiments that SUMO results in a better visualization of small structures in comparison to SGM. SUMO is applied to patient datasets for visualization of stent and prostate brachytherapy seeds. In addition, SUMO also provides quantitative information about the number of prostate brachytherapy seeds. The method might be extended to application for visualization of other interventional devices, and, like SGM, it might also be used to visualize magnetically labelled cells.