Unbalanced SSFP for super-resolution in MRI.

PURPOSE: To achieve rapid, low specific absorption rate (SAR) super-resolution imaging by exploiting the characteristic magnetization off-resonance profile in SSFP. THEORY AND METHODS: In the presented technique, low flip angle unbalanced SSFP imaging is used to acquire a series of images at a low nominal resolution that are then combined in a super-resolution strategy analogous to non-linear structured illumination microscopy. This is demonstrated in principle via Bloch simulations and synthetic phantoms, and the performance is quantified in terms of point-spread function (PSF) and SNR for gray and white matter from field strengths of 0.35T to 9.4T. A k-space reconstruction approach is proposed to account for B0 effects. This was applied to reconstruct super-resolution images from a test object at 9.4T. RESULTS: Artifact-free super-resolution images were produced after incorporating sufficient preparation time for the magnetization to approach the steady state. High-resolution images of a test object were obtained at 9.4T, in the presence of considerable B0 inhomogeneity. For gray matter, the highest achievable resolution ranges from 3% of the acquired voxel dimension at 0.35T, to 9% at 9.4T. For white matter, this corresponds to 3% and 10%, respectively. Compared to an equivalent segmented gradient echo acquisition at the optimal flip angle, with a fixed TR of 8 ms, gray matter has up to 34% of the SNR at 9.4T while using a ×10 smaller flip angle. For white matter, this corresponds to 29% with a ×11 smaller flip angle. CONCLUSION: This approach achieves high degrees of super-resolution enhancement with minimal RF power requirements.

Magnetic resonance imaging-based measurement of internal deformation of vibrating vocal fold models.

A method is presented for tracking the internal deformation of self-oscillating vocal fold models using magnetic resonance imaging (MRI). Silicone models scaled to four times life-size to lower the flow-induced vibration frequency were embedded with fiducial markers in a coronal plane. Candidate marker materials were tested using static specimens, and two materials, cupric sulfate and glass, were chosen for testing in the vibrating vocal fold models. The vibrating models were imaged using a gated MRI protocol wherein MRI acquisition was triggered using the subglottal pressure signal. Two-dimensional image slices at different phases during self-oscillation were captured, and in each phase the fiducial markers were clearly visible. The process was also demonstrated using a three-dimensional scan at two phases. The benefit of averaging to increase signal-to-noise ratio was explored. The results demonstrate the ability to use MRI to acquire quantitative deformation data that could be used, for example, to validate computational models of flow-induced vocal fold vibration and quantify deformation fields encountered by cells in bioreactor studies.

Feasibility of Diffusion Tensor and Morphologic Imaging of Peripheral Nerves at Ultra-High Field Strength.

OBJECTIVES: The aim of this study was to describe the development of morphologic and diffusion tensor imaging sequences of peripheral nerves at 7 T, using carpal tunnel syndrome (CTS) as a model system of focal nerve injury. MATERIALS AND METHODS: Morphologic images were acquired at 7 T using a balanced steady-state free precession sequence. Diffusion tensor imaging was performed using single-shot echo-planar imaging and readout-segmented echo-planar imaging sequences. Different acquisition and postprocessing methods were compared to describe the optimal analysis pipeline. Magnetic resonance imaging parameters including cross-sectional areas, signal intensity, fractional anisotropy (FA), as well as mean, axial, and radial diffusivity were compared between patients with CTS (n = 8) and healthy controls (n = 6) using analyses of covariance corrected for age (significance set at P < 0.05). Pearson correlations with Bonferroni correction were used to determine association of magnetic resonance imaging parameters with clinical measures (significance set at P < 0.01). RESULTS: The 7 T acquisitions with high in-plane resolution (0.2 × 0.2mm) afforded detailed morphologic resolution of peripheral nerve fascicles. For diffusion tensor imaging, single-shot echo-planar imaging was more efficient than readout-segmented echo-planar imaging in terms of signal-to-noise ratio per unit scan time. Distortion artifacts were pronounced, but could be corrected during postprocessing. Registration of FA maps to the morphologic images was successful. The developed imaging and analysis pipeline identified lower median nerve FA (pisiform bone, 0.37 [SD 0.10]) and higher radial diffusivity (1.08 [0.20]) in patients with CTS compared with healthy controls (0.53 [0.06] and 0.78 [0.11], respectively, P < 0.047). Fractional anisotropy and radial diffusivity strongly correlated with patients' symptoms (r = -0.866 and 0.866, respectively, P = 0.005). CONCLUSIONS: Our data demonstrate the feasibility of morphologic and diffusion peripheral nerve imaging at 7 T. Fractional anisotropy and radial diffusivity were found to be correlates of symptom severity.

Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank.

UK Biobank is a large-scale prospective epidemiological study with all data accessible to researchers worldwide. It is currently in the process of bringing back 100,000 of the original participants for brain, heart and body MRI, carotid ultrasound and low-dose bone/fat x-ray. The brain imaging component covers 6 modalities (T1, T2 FLAIR, susceptibility weighted MRI, Resting fMRI, Task fMRI and Diffusion MRI). Raw and processed data from the first 10,000 imaged subjects has recently been released for general research access. To help convert this data into useful summary information we have developed an automated processing and QC (Quality Control) pipeline that is available for use by other researchers. In this paper we describe the pipeline in detail, following a brief overview of UK Biobank brain imaging and the acquisition protocol. We also describe several quantitative investigations carried out as part of the development of both the imaging protocol and the processing pipeline.

Effect of slice excitation profile on ungated steady state cardiac perfusion imaging.

In cardiac perfusion imaging, choice of flip angle is an important factor for steady state acquisition. This work focuses on presenting an analytical framework for understanding how non-ideal slice excitation profiles affect contrast in ungated 2D steady state cardiac perfusion studies, and to study a technique for estimating flip angle that maximizes enhanced/unenhanced myocardial contrast-to-noise ratio (CNR) in single slice and multi-slice acquisitions. A numerical simulation of ungated 2D golden ratio radial spoiled gradient echo (SPGR) was created that takes into consideration the actual (Bloch simulated) slice excitation profile. The effect of slice excitation profile on myocardial CNR as a function of flip angle was assessed in phantoms and in-vivo. For fast RF pulses, the flip angle that yields maximum CNR (considering the actual slice excitation profile) was considerably higher than expected, assuming an ideal excitation. The simulation framework presented accurately predicts the flip angle yielding maximum CNR when the actual slice excitation profile is taken into consideration. The prescribed flip angle for optimal contrast in ungated 2D steady-state SPGR cardiac perfusion studies can vary significantly from that calculated when an ideal slice excitation profile is assumed. Consideration of the actual slice excitation can yield a more optimal flip angle estimate in both the single slice and multi-slice cases.

Multimodal population brain imaging in the UK Biobank prospective epidemiological study.

Medical imaging has enormous potential for early disease prediction, but is impeded by the difficulty and expense of acquiring data sets before symptom onset. UK Biobank aims to address this problem directly by acquiring high-quality, consistently acquired imaging data from 100,000 predominantly healthy participants, with health outcomes being tracked over the coming decades. The brain imaging includes structural, diffusion and functional modalities. Along with body and cardiac imaging, genetics, lifestyle measures, biological phenotyping and health records, this imaging is expected to enable discovery of imaging markers of a broad range of diseases at their earliest stages, as well as provide unique insight into disease mechanisms. We describe UK Biobank brain imaging and present results derived from the first 5,000 participants' data release. Although this covers just 5% of the ultimate cohort, it has already yielded a rich range of associations between brain imaging and other measures collected by UK Biobank.

Radial simultaneous multi-slice CAIPI for ungated myocardial perfusion.

OBJECTIVE: Simultaneous multi-slice (SMS) imaging is a slice acceleration technique that acquires multiple slices in the same time as a single slice. Radial controlled aliasing in parallel imaging results in higher acceleration (radial CAIPIRINHA or CAIPI) is a promising SMS method with less severe slice aliasing artifacts as compared to its Cartesian counterpart. Here we use radial CAIPI with data undersampling and constrained reconstruction to improve the utility of ungated cardiac perfusion acquisitions. We test the proposed framework with a traditional saturation recovery fast low-angle shot (turboFLASH) sequence and also without saturation recovery as a steady-state spoiled gradient echo (SPGR) sequence on animal and human studies. METHODS: Simulations and phantom studies were performed for both the turboFLASH and the SPGR radial CAIPI methods. Ungated undersampled golden ratio radial CAIPI data with saturation recovery were acquired in 8 dogs and 2 human subjects. The CAIPI data without saturation pulses were acquired in 4 human subjects. For both methods, slice acceleration factors of two and three were used. A new spatio-temporal reconstruction using total variation and patch-based low rank constraints was used to jointly reconstruct the multi-slice multi-coil images. RESULTS: Phantom scans and computer simulations showed that ungated SPGR generally provides better contrast to noise ratio (CNR) than the saturation recovery sequence if the saturation recovery time is less than 100ms. Both of the ungated radial CAIPI methods demonstrated promising image quality in terms of preserving dynamics of the contrast agent and maintaining anatomical structures, even with three slices acquired simultaneously. CONCLUSION: Ungated simultaneous multi-slice acquisitions with either a saturation recovery turboFLASH sequence or a steady-state gradient echo SPGR sequence are feasible and provide increased slice coverage without loss of temporal resolution. Compared with a sensitivity encoding (SENSE) SMS reconstruction, the constrained reconstruction method provides better image quality for undersampled radial CAIPI data.

Quantitative sodium magnetic resonance imaging of cartilage, muscle, and tendon.

Sodium magnetic resonance imaging (MRI), or imaging of the 23Na nucleus, has been under exploration for several decades, and holds promise for potentially revealing additional biochemical information about the health of tissues that cannot currently be obtained from conventional hydrogen (or proton) MRI. This additional information could serve as an important complement to conventional MRI for many applications. However, despite these exciting possibilities, sodium MRI is not yet used routinely in clinical practice, and will likely remain strictly in the domain of exploratory research for the coming decade. This paper begins with a technical overview of sodium MRI, including the nuclear magnetic resonance (NMR) signal characteristics of the sodium nucleus, the challenges associated with sodium MRI, and the specialized pulse sequences, hardware, and reconstruction techniques required. Various applications of sodium MRI for quantitative analysis of the musculoskeletal system are then reviewed, including the non-invasive assessment of cartilage degeneration in vivo, imaging of tendinopathy, applications in the assessment of various muscular pathologies, and assessment of muscle response to exercise.

Quantitative techniques for musculoskeletal MRI at 7 Tesla.

Whole-body 7 Tesla MRI scanners have been approved solely for research since they appeared on the market over 10 years ago, but may soon be approved for selected clinical neurological and musculoskeletal applications in both the EU and the United States. There has been considerable research work on musculoskeletal applications at 7 Tesla over the past decade, including techniques for ultra-high resolution morphological imaging, 3D T2 and T2* mapping, ultra-short TE applications, diffusion tensor imaging of cartilage, and several techniques for assessing proteoglycan content in cartilage. Most of this work has been done in the knee or other extremities, due to technical difficulties associated with scanning areas such as the hip and torso at 7 Tesla. In this manuscript, we first provide some technical context for 7 Tesla imaging, including challenges and potential advantages. We then review the major quantitative MRI techniques being applied to musculoskeletal applications on 7 Tesla whole-body systems.

Sodium MRI radiofrequency coils for body imaging.

The proliferation of high-field whole-body systems, advances in gradient performance and refinement of signal-to-noise ratio (SNR)-efficient short-TE sequences suitable for sodium imaging have led to a resurgence of interest in sodium imaging for body applications. With this renewed interest has come increased demand for SNR-efficient sodium coils. Efficient coils can significantly increase SNR in sodium imaging, allowing higher resolutions and/or shorter scan times. In this work, we focus on body imaging applications of sodium MRI, and review developments in MRI radiofrequency (RF) coil topologies for sodium imaging. We first provide a brief discussion of RF coil design considerations in sodium imaging. This is followed by an overview of common coil topologies, their advantages and disadvantages, and examples of each.

Methodology for computing white matter nerve fiber orientation in human histological slices.

BACKGROUND: The gold standard for mapping nerve fiber orientation in white matter of the human brain is histological analysis through biopsy. Such mappings are a crucial step in validating non-invasive techniques for assessing nerve fiber orientation in the human brain by using diffusion MRI. However, the manual extraction of nerve fiber directions of histological slices is tedious, time consuming, and prone to human error. NEW METHOD: The presented semi-automated algorithm first creates a binary-segmented mask of the nerve fibers in the histological image, and then extracts an estimate of average directionality of nerve fibers through a Fourier-domain analysis of the masked image. It also generates an uncertainty level for its estimate of average directionality. RESULTS AND COMPARISON WITH EXISTING METHODS: The average orientations of the semi-automatic method were first compared to a qualitative expert opinion based on visual inspection of nerve fibers. A weighted RMS difference between the expert estimate and the algorithmically determined angle (weighted by expert's confidence in his estimate) was 15.4°, dropping to 9.9° when only cases with an expert confidence level of greater than 50% were included. The algorithmically determined angles were then compared with angles extracted using a manual segmentation technique, yielding an RMS difference of 11.2°. CONCLUSION: The presented semi-automated method is in good agreement with both qualitative and quantitative manual expert-based approaches for estimating directionality of nerve fibers in white matter from images of stained histological slices of the human brain.

Synchronous Radial 1H and 23Na Dual-Nuclear MRI on a Clinical MRI System, Equipped With a Broadband Transmit Channel.

The purpose of this work was to synchronously acquire proton (1H) and sodium (23Na) image data on a 3T clinical MRI system within the same sequence, without internal modification of the clinical hardware, and to demonstrate synchronous acquisition with 1H/23Na-GRE imaging with Cartesian and radial k-space sampling. Synchronous dual-nuclear imaging was implemented by: mixing down the 1H signal so that both the 23Na and 1H signal were acquired at 23Na frequency by the conventional MRI system; interleaving 1H/23Na transmit pulses in both Cartesian and radial sequences; and using phase stabilization on the 1H signal to remove mixing effects. The synchronous 1H/23Na setup obtained images in half the time necessary to sequentially acquire the same 1H and 23Na images with the given setup and parameters. Dual-nuclear hardware and sequence modifications were used to acquire 23Na images within the same sequence as 1H images, without increases to the 1H acquisition time. This work demonstrates a viable technique to acquire 23Na image data without increasing 1H acquisition time using minor additional custom hardware, without requiring modification of a commercial scanner with multinuclear capability.

Accuracy of 3D dual echo steady state (DESS) MR arthrography to quantify acetabular cartilage thickness.

PURPOSE: To deploy and quantify the accuracy of 3D dual echo steady state (DESS) MR arthrography with hip traction to image acetabular cartilage. Clinical magnetic resonance imaging (MRI) sequences used to image hip cartilage often have reduced out-of-plane resolution and may lack adequate signal-to-noise to image cartilage. MATERIALS AND METHODS: Saline was injected into four cadaver hips placed under traction. 3D DESS MRI scans were obtained before and after cores of cartilage were harvested from the acetabulum; the two MRIs were spatially aligned to reference core positions. The thickness of cartilage cores was measured under microscopy to serve as the reference standard. 3D reconstructions of cartilage and subchondral bone were generated using automatic and semiautomatic image segmentation. Cartilage thickness estimated from the 3D reconstructions was compared to physical measurements using Bland-Altman plots. RESULTS: As revealed by the automatic segmentation mask, saline imbibed the joint space throughout the articulating surface, with the exception of the posteroinferior region in two hips. Locations where air bubbles were introduced and regions of suspected low density bone disrupted an otherwise smooth automatic segmentation mask. Automatic and semiautomatic segmentation yielded a bias ± repeatability coefficient (95% limits of agreement) of 0.10 ± 0.51 mm (-0.41 to 0.61 mm) and 0.06 ± 0.43 mm (-0.37 to 0.49 mm), respectively. CONCLUSION: Cartilage thickness can be estimated to within ∼0.5 mm of the physical value with 95% confidence using 3D reconstructions of 3D DESS MR arthrography images. Manual correction of the automatic segmentation mask may improve reconstruction accuracy.

Comparison of centric and reverse-centric trajectories for highly accelerated three-dimensional saturation recovery cardiac perfusion imaging.

PURPOSE: Highly undersampled three-dimensional (3D) saturation-recovery sequences are affected by k-space trajectory since the magnetization does not reach steady state during the acquisition and the slab excitation profile yields different flip angles in different slices. This study compares centric and reverse-centric 3D cardiac perfusion imaging. METHODS: An undersampled (98 phase encodes) 3D ECG-gated saturation-recovery sequence that alternates centric and reverse-centric acquisitions each time frame was used to image phantoms and in vivo subjects. Flip angle variation across the slices was measured, and contrast with each trajectory was analyzed via Bloch simulation. RESULTS: Significant variations in flip angle were observed across slices, leading to larger signal variation across slices for the centric acquisition. In simulation, severe transient artifacts were observed when using the centric trajectory with higher flip angles, placing practical limits on the maximum flip angle used. The reverse-centric trajectory provided less contrast, but was more robust to flip angle variations. CONCLUSION: Both of the k-space trajectories can provide reasonable image quality. The centric trajectory can have higher CNR, but is more sensitive to flip angle variation. The reverse-centric trajectory is more robust to flip angle variation.

A 3 T sodium and proton composite array breast coil.

PURPOSE: The objective of this study was to determine whether a sodium phased array would improve sodium breast MRI at 3 T. The secondary objective was to create acceptable proton images with the sodium phased array in place. METHODS: A novel composite array for combined proton/sodium 3 T breast MRI is compared with a coil with a single proton and sodium channel. The composite array consists of a 7-channel sodium receive array, a larger sodium transmit coil, and a 4-channel proton transceive array. The new composite array design utilizes smaller sodium receive loops than typically used in sodium imaging, uses novel decoupling methods between the receive loops and transmit loops, and uses a novel multichannel proton transceive coil. The proton transceive coil reduces coupling between proton and sodium elements by intersecting the constituent loops to reduce their mutual inductance. The coil used for comparison consists of a concentric sodium and proton loop with passive decoupling traps. RESULTS: The composite array coil demonstrates a 2-5× improvement in signal-to-noise ratio for sodium imaging and similar signal-to-noise ratio for proton imaging when compared with a simple single-loop dual resonant design. CONCLUSION: The improved signal-to-noise ratio of the composite array gives breast sodium images of unprecedented quality in reasonable scan times.

A statistical analysis of the Bloch-Siegert B1 mapping technique.

A number of B1 mapping methods have been introduced. A model to facilitate choice among these methods is valuable, as the performance of each technique is affected by a variety of factors, including acquisition signal-to-noise ratio (SNR). The Bloch-Siegert shift B1 mapping method has recently garnered significant interest. In this paper, we present a statistical model suitable for analysis of the Bloch-Siegert shift method. Unlike previously presented models, the analysis is valid in both low SNR and high SNR regimes. We present a detailed analysis of the performance of the Bloch-Siegert shift B1 mapping method across a broad range of acquisition scenarios, and compare it to two other B1 mapping techniques (the dual angle method and the phase sensitive method). Further validation of the model is presented through both Monte Carlo simulations and experimental results. The simulations and experimental results match the model well, lending confidence to its accuracy. Each technique is found to perform well with high acquisition SNR. However, our results suggest that the dual angle method is not reliable in low SNR environments. Furthermore, the phase sensitive method appears to outperform the Bloch-Siegert shift method in these low-SNR cases, although variations of the Bloch-Siegert method may be possible that improve its performance at low SNR.

Musculoskeletal MRI at 3.0 T and 7.0 T: a comparison of relaxation times and image contrast.

OBJECTIVE: The purpose of this study was to measure and compare the relaxation times of musculoskeletal tissues at 3.0 T and 7.0 T, and to use these measurements to select appropriate parameters for musculoskeletal protocols at 7.0 T. MATERIALS AND METHODS: We measured the T1 and T2 relaxation times of cartilage, muscle, synovial fluid, bone marrow and subcutaneous fat at both 3.0 T and 7.0 T in the knees of five healthy volunteers. The T1 relaxation times were measured using a spin-echo inversion recovery sequence with six inversion times. The T2 relaxation times were measured using a spin-echo sequence with seven echo times. The accuracy of both the T1 and T2 measurement techniques was verified in phantoms at both magnetic field strengths. We used the measured relaxation times to help design 7.0 T musculoskeletal protocols that preserve the favorable contrast characteristics of our 3.0 T protocols, while achieving significantly higher resolution at higher SNR efficiency. RESULTS: The T1 relaxation times in all tissues at 7.0 T were consistently higher than those measured at 3.0 T, while the T2 relaxation times at 7.0 T were consistently lower than those measured at 3.0 T. The measured relaxation times were used to help develop high resolution 7.0 T protocols that had similar fluid-to-cartilage contrast to that of the standard clinical 3.0 T protocols for the following sequences: proton-density-weighted fast spin-echo (FSE), T2-weighted FSE, and 3D-FSE-Cube. CONCLUSION: The T1 and T2 changes were within the expected ranges. Parameters for musculoskeletal protocols at 7.0 T can be optimized based on these values, yielding improved resolution in musculoskeletal imaging with similar contrast to that of standard 3.0 T clinical protocols.

Simultaneous fat suppression and band reduction with large-angle multiple-acquisition balanced steady-state free precession.

Balanced steady-state free precession (bSSFP) MRI is a rapid and signal-to-noise ratio-efficient imaging method, but suffers from characteristic bands of signal loss in regions of large field inhomogeneity. Several methods have been developed to reduce the severity of these banding artifacts, typically involving the acquisition of multiple bSSFP datasets (and the accompanying increase in scan time). Fat suppression with bSSFP is also challenging; most existing methods require an additional increase in scan time, and some are incompatible with bSSFP band-reduction techniques. This work was motivated by the need for both robust fat suppression and band reduction in the presence of field inhomogeneity when using bSSFP for flow-independent peripheral angiography. The large flip angles used in this application to improve vessel conspicuity and contrast lead to specific absorption rate considerations, longer repetition times, and increased severity of banding artifacts. In this work, a novel method that simultaneously suppresses fat and reduces bSSFP banding artifact with the acquisition of only two phase-cycled bSSFP datasets is presented. A weighted sum of the two bSSFP acquisitions is taken on a voxel-by-voxel basis, effectively synthesizing an off-resonance profile at each voxel that puts fat in the stop band while keeping water in the pass band. The technique exploits the near-sinusoidal shape of the bSSFP off-resonance spectrum for many tissues at large (>50°) flip angles.

Three-dimensional fluid-suppressed T2-prep flow-independent peripheral angiography using balanced SSFP.

Accurate depiction of the vessels of the lower leg, foot or hand benefits from suppression of bright MR signal from lipid (such as bone marrow) and long-T1 fluid (such as synovial fluid and edema). Signal independence of blood flow velocities, good arterial/muscle contrast and arterial/venous separation are also desirable. The high SNR, short scan times and flow properties of balanced steady-state free precession (SSFP) make it an excellent candidate for flow-independent angiography. In this work, a new magnetization-prepared 3D SSFP sequence for flow-independent peripheral angiography is presented. The technique combines a number of component techniques (phase-sensitive fat detection, inversion recovery, T2-preparation and square-spiral phase-encode ordering) to achieve high-contrast peripheral angiograms at only a modest scan time penalty over simple 3D SSFP. The technique is described in detail, a parameter optimization performed and preliminary results presented achieving high contrast and 1-mm isotropic resolution in a normal foot.

Phase-sensitive sodium B1 mapping.

Quantitative sodium MRI requires accurate knowledge of factors affecting the sodium signal. One important determinant of sodium signal level is the transmit B(1) field strength. However, the low signal-to-noise ratio typical of sodium MRI makes accurate B(1) mapping in reasonable scan times challenging. A new phase-sensitive B(1) mapping technique has recently been shown to work better than the widely used dual-angle method in low-signal-to-noise ratio situations and over a broader range of flip angles. In this work, the phase-sensitive B(1) mapping technique is applied to sodium, and its performance compared to the dual-angle method through both simulation and phantom studies. The phase-sensitive method is shown to yield higher quality B(1) maps at low signal-to-noise ratio and greater consistency of measurement than the dual-angle method. An in vivo sodium B(1) map of the human breast is also shown, demonstrating the phase-sensitive method's feasibility for human studies.