Correcting image blur in spiral, retraced in/out (RIO) acquisitions using a maximized energy objective.

PURPOSE: Images acquired with spiral k-space trajectories can suffer from off-resonance image blur. Previous work showed that averaging 2 images acquired with a retraced, in/out (RIO) trajectory self-corrects image blur so long as off-resonant spins accrue less than 1 half-cycle of relative phase over the readout. Practical scenarios frequently exceed this threshold. Here, we derive and characterize a more-robust off-resonance image blur correction method for RIO acquisitions. METHODS: Phantom and human volunteer data were acquired using a RIO trajectory with readout durations ranging from 4 to 60 ms. The resulting images were deblurred using 3 candidate methods: conventional linear correction of the component images; semiautomatic deblurring of the component images using an established minimized phase objective function; and semiautomatic deblurring of the average of the component images using a maximized energy objective function, derived below. Deblurring errors were estimated relative to images acquired with 4 ms readouts. RESULTS: All 3 methods converged to similar solutions in cases where less than 2 and 4 cycles of phase accrued over the readout in in vivo and phantom images, respectively (<13 ms readout at 3T). Above this threshold, the linear and minimized phase methods introduced several errors. The maximized energy function provided accurate deblurring so long as less than 6 and 10 cycles of phase accrued over the readout in in vivo and phantom images, respectively (<34 ms readout at 3T). CONCLUSION: The maximized energy objective function can accurately deblur RIO acquisitions over a wide spectrum of off resonance frequencies.

A spiral-based volumetric acquisition for MR temperature imaging.

PURPOSE: To develop a rapid pulse sequence for volumetric MR thermometry. METHODS: Simulations were carried out to assess temperature deviation, focal spot distortion/blurring, and focal spot shift across a range of readout durations and maximum temperatures for Cartesian, spiral-out, and retraced spiral-in/out (RIO) trajectories. The RIO trajectory was applied for stack-of-spirals 3D imaging on a real-time imaging platform and preliminary evaluation was carried out compared to a standard 2D sequence in vivo using a swine brain model, comparing maximum and mean temperatures measured between the two methods, as well as the temporal standard deviation measured by the two methods. RESULTS: In simulations, low-bandwidth Cartesian trajectories showed substantial shift of the focal spot, whereas both spiral trajectories showed no shift while maintaining focal spot geometry. In vivo, the 3D sequence achieved real-time 4D monitoring of thermometry, with an update time of 2.9-3.3 s. CONCLUSION: Spiral imaging, and RIO imaging in particular, is an effective way to speed up volumetric MR thermometry. Magn Reson Med 79:3122-3127, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Automated Segmentation of Head Computed Tomography Images Using FSL.

OBJECTIVE: The aim of this study was to investigate the use of one magnetic resonance image-processing tool, FSL, in its ability to perform automated segmentation of computed tomographic images of the brain. METHODS: Head computed tomography (CT) images were brain extracted and segmented using the FSL tools BET and FAST, respectively. The products of segmentation were analyzed by histogram. The impact of image intensity inhomogeneity correction was investigated using simulated bias fields, 14 routine head CT scans, and selected illustrative clinical cases. RESULTS: FSL FAST performs direct segmentation of head CT images, permitting quantitation of gray and white matter densities and volumes, achieving a more complete segmentation than masking methods. "Bias field correction" reduced the covariance of image signal intensities of the total brain and gray matter images (P < 0.01). Correction is larger when the effects of beam hardening and radiation scatter are larger, resulting in improved segmentation. CONCLUSIONS: FSL FAST enables direct segmentation of head CT images.

Rapid 3D dynamic arterial spin labeling with a sparse model-based image reconstruction.

Dynamic arterial spin labeling (ASL) MRI measures the perfusion bolus at multiple observation times and yields accurate estimates of cerebral blood flow in the presence of variations in arterial transit time. ASL has intrinsically low signal-to-noise ratio (SNR) and is sensitive to motion, so that extensive signal averaging is typically required, leading to long scan times for dynamic ASL. The goal of this study was to develop an accelerated dynamic ASL method with improved SNR and robustness to motion using a model-based image reconstruction that exploits the inherent sparsity of dynamic ASL data. The first component of this method is a single-shot 3D turbo spin echo spiral pulse sequence accelerated using a combination of parallel imaging and compressed sensing. This pulse sequence was then incorporated into a dynamic pseudo continuous ASL acquisition acquired at multiple observation times, and the resulting images were jointly reconstructed enforcing a model of potential perfusion time courses. Performance of the technique was verified using a numerical phantom and it was validated on normal volunteers on a 3-Tesla scanner. In simulation, a spatial sparsity constraint improved SNR and reduced estimation errors. Combined with a model-based sparsity constraint, the proposed method further improved SNR, reduced estimation error and suppressed motion artifacts. Experimentally, the proposed method resulted in significant improvements, with scan times as short as 20s per time point. These results suggest that the model-based image reconstruction enables rapid dynamic ASL with improved accuracy and robustness.

A simple acquisition strategy to avoid off-resonance blurring in spiral imaging with redundant spiral-in/out k-space trajectories.

PURPOSE: The major hurdle to widespread adoption of spiral trajectories has been their poor off-resonance performance. Here we present a self-correcting spiral k-space trajectory that avoids much of the well-known spiral blurring during data acquisition. THEORY AND METHODS: In comparison with a traditional spiral-out trajectory, the spiral-in/out trajectory has improved off-resonance performance. By combining two spiral-in/out acquisitions, one rotated 180° in k-space compared with the other, multishot spiral-in/out artifacts are eliminated. A phantom was scanned with the center frequency manually tuned 20, 40, 80, and 160 Hz off-resonance with both a spiral-out gradient echo sequence and the redundant spiral-in/out sequence. The phantom was also imaged in an oblique orientation in order to demonstrate improved concomitant gradient field performance of the sequence. Additionally, the trajectory was incorporated into a spiral turbo spin echo sequence for brain imaging. RESULTS: Phantom studies with manually tuned off-resonance agree well with theoretical calculations, showing that moderate off-resonance is well-corrected by this acquisition scheme. Blur due to concomitant fields is reduced, and good results are obtained in vivo. CONCLUSION: The redundant spiral-in/out trajectory results in less image blur for a given readout length than a traditional spiral-out scan, reducing the need for complex off-resonance correction algorithms.

Noncontrast peripheral MRA with spiral echo train imaging.

PURPOSE: To develop a spin echo train sequence with spiral readout gradients with improved artery-vein contrast for noncontrast angiography. THEORY: Venous T2 becomes shorter as the echo spacing is increased in echo train sequences, improving contrast. Spiral acquisitions, due to their data collection efficiency, facilitate long echo spacings without increasing scan times. METHODS: Bloch equation simulations were performed to determine optimal sequence parameters, and the sequence was applied in five volunteers. In two volunteers, the sequence was performed with a range of echo times and echo spacings to compare with the theoretical contrast behavior. A Cartesian version of the sequence was used to compare contrast appearance with the spiral sequence. Additionally, spiral parallel imaging was optionally used to improve image resolution. RESULTS: In vivo, artery-vein contrast properties followed the general shape predicted by simulations, and good results were obtained in all stations. Compared with a Cartesian implementation, the spiral sequence had superior artery-vein contrast, better spatial resolution (1.2 mm(2) versus 1.5 mm(2) ), and was acquired in less time (1.4 min versus 7.5 min). CONCLUSION: The spiral spin echo train sequence can be used for flow-independent angiography to generate three-dimensional angiograms of the periphery quickly and without the use of contrast agents.

Refocused turbo spin-echo for noncontrast peripheral MR angiography.

PURPOSE: To develop and assess a three-dimensional refocused turbo spin-echo (rTSE) sequence for generating peripheral angiograms. This sequence combines the rapid T2 -weighting of TSE and the better flow performance of the fully-refocused gradients of balanced steady state free precession (bSSFP), along with bSSFP-style phase alternation of refocusing radiofrequency (RF) pulses. MATERIALS AND METHODS: The signal behavior generated by such a sequence was explored through Bloch equation simulations. The rTSE and TSE sequences were both used to generate peripheral angiograms in nine normal volunteers. The signal to noise ratio, contrast resolution, and vessel sharpness of the resulting images were used as bases for comparison. Additionally, the rTSE sequence was applied in four patients with peripheral artery disease to preliminarily assess its efficacy in a clinical setting through quality scoring by two experienced radiologists. RESULTS: The rTSE's RF phase alternation approach out-performs a simple balanced-gradient CPMG (Carr-Purcell-Meiboom-Gill) -style TSE sequence in the presence of B0 and B1 inhomogeneities. In volunteers, the rTSE sequence yielded better arterial-venous contrast (0.378 ± 0.145 versus 0.155 ± 0.202; P < 0.01) and increased vessel sharpness (0.340 ± 0.034 versus 0.263 ± 0.034; P < 0.005) over TSE images. Stenoses visible in conventional angiographic images in patients were successfully imaged with the rTSE sequence; however, image quality scores in patients were lower than in volunteers (1.2 ± 0.38 versus 3.0 ± 1.0; P < 0.05). CONCLUSION: The rTSE sequence generates nonsubtractive, flow-independent, peripheral MR angiograms with better arterial-venous contrast and vessel sharpness in normal volunteers than a conventional TSE sequence.

AAPM Task Group Report 325: MRI static magnetic field homogeneity measurement and evaluation procedures - Guidance and resources.

Measurement of static magnetic field (B0) homogeneity is an essential component of routine MRI system evaluation. This report summarizes the work of AAPM Task Group (TG) 325 on vendor-specific methods of B0 homogeneity measurement and evaluation. TG 325 was charged with producing a set of detailed, step-by-step instructions to implement B0 homogeneity measurement methods discussed in the American College of Radiology (ACR) MRI Quality Control Manual using specific makes and models of MRI scanners. The TG produced such instructions for as many approaches as was relevant and practical on six currently available vendor platforms including details of software/tools, settings, phantoms, and other experimental details needed for a reproducible protocol. Because edits to these instructions may need to be made as vendors enter and exit the market and change available tools, interfaces, and access levels over time, the step-by-step instructions are published as a living document on the AAPM website. This summary document provides an introduction to B0 homogeneity testing in MRI and several of the common methods for its measurement and evaluation. A living document on the AAPM website provides vendor-specific step-by-step instructions for performing these tests to facilitate accurate and reproducible B0 homogeneity evaluation on a routine basis.

The effects of the skull on CT imaging of the brain: a skull and brain phantom study.

OBJECTIVE: To investigate the effects of beam hardening by the skull on the measured radiodensity of the brain. To test a hypothesis that these effects of beam hardening are decreased using a monochromatic energy source. METHODS: Selected clinical cases were reviewed in illustration. An anthropomorphic skull and brain phantom was created and scanned in a clinical CT scanner with skull, without skull, and with hemicraniectomy. The effects of beam hardening were illustrated by scanning the phantom with mono- and poly-chromatic X-ray sources. RESULTS: In clinical cases, the HU values of the brain were consistently lower when the X-ray beam traversed the skull than when it did not. An anthropomorphic skull-and-brain phantom further demonstrated these effects, which were evident with a polychromatic energy source and absent with a virtual monochromatic energy source. CONCLUSIONS: Beam hardening by the skull lowers the measured HU values of the brain. The effects, which can impact quantitative imaging, may be mitigated by a virtual monochromatic energy source. ADVANCES IN KNOWLEDGE: Beam hardening by the skull lowers the measured radiodensity of the brain. The effects may be mitigated by a virtual monochromatic energy source.

3D-Encoded DENSE MRI with Zonal Excitation for Quantifying Biventricular Myocardial Strain During a Breath-Hold.

PURPOSE: Right ventricular (RV) function is increasingly recognized for its prognostic value in many disease states. As with the left ventricle (LV), strain-based measurements may have better prognostic value than typical chamber volumes or ejection fraction. Complete functional characterization of the RV requires high-resolution, 3D displacement tracking methods, which have been prohibitively challenging to implement. Zonal excitation during Displacement ENcoding with Stimulated Echoes (DENSE) magnetic resonance imaging (MRI) has helped reduce scan time for 2D LV strain quantification. We hypothesized that zonal excitation could alternatively be used to reproducibly acquire higher resolution, 3D-encoded DENSE images for quantification of bi-ventricular strain within a single breath-hold. METHODS: We modified sequence parameters for a 3D zonal excitation DENSE sequence to achieve in-plane resolution < 2 mm and acquired two sets of images in eight healthy adult male volunteers with median (IQR) age 32.5 (32.0-33.8) years. We assessed the inter-test reproducibility of this technique, and compared computed strains and torsion with previously published data. RESULTS: Data for one subject was excluded based on image artifacts. Reproducibility for LV (CoV: 6.1-9.0%) and RV normal strains (CoV: 6.3-8.2%) and LV torsion (CoV = 7.1%) were all very good. Reproducibility of RV torsion was lower (CoV = 16.7%), but still within acceptable limits. Computed global strains and torsion were within reasonable agreement with published data, but further studies in larger cohorts are needed to confirm. CONCLUSION: Reproducible acquisition of 3D-encoded biventricular myocardial strain data in a breath-hold is feasible using DENSE with zonal excitation.

Genomic Screening for Pathogenic Transthyretin Variants Finds Evidence of Underdiagnosed Amyloid Cardiomyopathy From Health Records.

BACKGROUND: New treatments for transthyretin amyloidosis improve survival, but diagnosis remains challenging. Pathogenic or likely pathogenic (P/LP) variants in the transthyretin (TTR) gene are one cause of transthyretin amyloidosis, and genomic screening has been proposed to identify at-risk individuals. However, data on disease features and penetrance are lacking to inform the utility of such population-based genomic screening for TTR. OBJECTIVES: This study characterized the prevalence of P/LP variants in TTR identified through exome sequencing and the burden of associated disease from electronic health records for individuals with these variants from a large (N = 134,753), primarily European-ancestry cohort. METHODS: We compared frequencies of common disease features and cardiac imaging findings between individuals with and without P/LP TTR variants. RESULTS: We identified 157 of 134,753 (0.12%) individuals with P/LP TTR variants (43% male, median age 52 [Q1-Q3: 37-61] years). Seven P/LP variants accounted for all observations, the majority being V122I (p.V142I; 113, 0.08%). Approximately 60% (n = 91) of individuals with P/LP TTR variants (all V122I) had African ancestry. Diagnoses of amyloidosis were limited (2 of 157 patients), although related heart disease diagnoses, including cardiomyopathy and heart failure, were significantly increased in individuals with P/LP TTR variants who were aged >60 years. Fourteen percent (7 of 49) of individuals aged ≥60 or older with a P/LP TTR variant had heart disease and ventricular septal thickness >1.2 cm, only one of whom was diagnosed with amyloidosis. CONCLUSIONS: Individuals with P/LP TTR variants identified by genomic screening have increased odds of heart disease after age 60 years, although amyloidosis is likely underdiagnosed without knowledge of the genetic variant.

Automated Estimation of Acute Infarct Volume from Noncontrast Head CT Using Image Intensity Inhomogeneity Correction.

Identification of early ischemic changes (EIC) on noncontrast head CT scans performed within the first few hours of stroke onset may have important implications for subsequent treatment, though early stroke is poorly delimited on these studies. Lack of sharp lesion boundary delineation in early infarcts precludes manual volume measures, as well as measures using edge-detection or region-filling algorithms. We wished to test a hypothesis that image intensity inhomogeneity correction may provide a sensitive method for identifying the subtle regional hypodensity which is characteristic of early ischemic infarcts. A digital image analysis algorithm was developed using image intensity inhomogeneity correction (IIC) and intensity thresholding. Two different IIC algorithms (FSL and ITK) were compared. The method was evaluated using simulated infarcts and clinical cases. For synthetic infarcts, measured infarct volumes demonstrated strong correlation to the true lesion volume (for 20% decreased density "infarcts," Pearson r = 0.998 for both algorithms); both algorithms demonstrated improved accuracy with increasing lesion size and decreasing lesion density. In clinical cases (41 acute infarcts in 30 patients), calculated infarct volumes using FSL IIC correlated with the ASPECTS scores (Pearson r = 0.680) and the admission NIHSS (Pearson r = 0.544). Calculated infarct volumes were highly correlated with the clinical decision to treat with IV-tPA. Image intensity inhomogeneity correction, when applied to noncontrast head CT, provides a tool for image analysis to aid in detection of EIC, as well as to evaluate and guide improvements in scan quality for optimal detection of EIC.

A Radiodensity Histogram Study of the Brain in Multiple Sclerosis.

Multiple sclerosis (MS) is a progressive neurodegenerative disease, affecting 1 million Americans and 2.5 million people globally. Although the diagnosis is made clinically, imaging plays a major role in diagnosing and monitoring disease progression and treatment response. Magnetic resonance imaging (MRI) has proven sensitive in imaging MS lesions, but the characterization offered by routine clinical MRI remains qualitative and with discrepancies between imaging and clinical findings. We investigated the ability of digital analysis of noncontrast head computed tomography (CT) images to detect global brain changes of MS. All routine diagnostic head CTs obtained on patients with known MS obtained from 1 of 2 scan platforms from 6/1/2011 to 6/1/2015 were reviewed. Head CT images from 54 patients with MS met inclusion criteria. Head CT images were processed and histogram metrics were compared to age- and gender- matched control subjects from the same CT scanners during the same time interval. Histogram metrics were correlated with plaque burden as seen on MRI studies. Compared with control subjects, patients had increased total brain radiodensity (P < .0001), further characterized as an increased histogram modal radiodensity (P < .0001) with decrease in histogram skewness (P < .0001). Radiodensity decreased with increasing plaque burden. Similar findings were seen in the patients with only mild plaque burden sub- group. Radiodensity is a unique tissue metric that is not measured by other imaging techniques. Our study finds that brain radiodensity histogram metrics highly correlate with MS, even in cases with minimal plaque burden.

Typical readout durations in spiral cine DENSE yield blurred images and underestimate cardiac strains at both 3.0 T and 1.5 T.

INTRODUCTION: Displacement encoding with stimulated echoes (DENSE) is a phase contrast technique that encodes tissue displacement into phase images, which are typically processed into measures of cardiac function such as strains. For improved signal to noise ratio and spatiotemporal resolution, DENSE is often acquired with a spiral readout using an 11.1 ms readout duration. However, long spiral readout durations are prone to blurring due to common phenomena such as off-resonance and T2* decay, which may alter the resulting quantifications of strain. We hypothesized that longer readout durations would reduce image quality and underestimate cardiac strains at both 3.0 T and 1.5 T and that using short readout durations could overcome these limitations. MATERIAL AND METHODS: Computational simulations were performed to investigate the relationship between off-resonance and T2* decay, the spiral cine DENSE readout duration, and measured radial and circumferential strain. Five healthy participants subsequently underwent 2D spiral cine DENSE at both 3.0 T and 1.5 T with several different readout durations 11.1 ms and shorter. Pearson correlations were used to assess the relationship between cardiac strains and the spiral readout duration. RESULTS: Simulations demonstrated that long readout durations combined with off-resonance and T2* decay yield blurred images and underestimate strains. With the typical 11.1 ms DENSE readout, blurring was present in the anterior and lateral left ventricular segments of participants and was markedly improved with shorter readout durations. Radial and circumferential strains from those segments were significantly correlated with the readout duration. Compared to the 1.9 ms readout, the 11.1 ms readout underestimated radial and circumferential strains in those segments at both field strengths by up to 19.6% and 1.5% (absolute), or 42% and 7% (relative), respectively. CONCLUSIONS: Blurring is present in spiral cine DENSE images acquired at both 3.0 T and 1.5 T using the typical 11.1 ms readout duration, which yielded substantially reduced radial strains and mildly reduced circumferential strains. Clinical studies using spiral cine DENSE should consider these limitations, while future technical advances may need to leverage accelerated techniques to improve the robustness and accuracy of the DENSE acquisition rather than focusing solely on reduced acquisition time.

Spatially-segmented undersampled MRI temperature reconstruction for transcranial MR-guided focused ultrasound.

BACKGROUND: Volumetric thermometry with fine spatiotemporal resolution is desirable to monitor MR-guided focused ultrasound (MRgFUS) procedures in the brain, but requires some form of accelerated imaging. Accelerated MR temperature imaging methods have been developed that undersample k-space and leverage signal correlations over time to suppress the resulting undersampling artifacts. However, in transcranial MRgFUS treatments, the water bath surrounding the skull creates signal variations that do not follow those correlations, leading to temperature errors in the brain due to signal aliasing. METHODS: To eliminate temperature errors due to the water bath, a spatially-segmented iterative reconstruction method was developed. The method fits a k-space hybrid signal model to reconstruct temperature changes in the brain, and a conventional MR signal model in the water bath. It was evaluated using single-channel 2DFT Cartesian, golden angle radial, and spiral data from gel phantom heating, and in vivo 8-channel 2DFT data from a FUS thalamotomy. Water bath signal intensity in phantom heating images was scaled between 0-100% to investigate its effect on temperature error. Temperature reconstructions of retrospectively undersampled data were performed using the spatially-segmented method, and compared to conventional whole-image k-space hybrid (phantom) and SENSE (in vivo) reconstructions. RESULTS: At 100% water bath signal intensity, 3 ×-undersampled spatially-segmented temperature reconstruction error was nearly 5-fold lower than the whole-image k-space hybrid method. Temperature root-mean square error in the hot spot was reduced on average by 27 × (2DFT), 5 × (radial), and 12 × (spiral) using the proposed method. It reduced in vivo error 2 × in the brain for all acceleration factors, and between 2 × and 3 × in the temperature hot spot for 2-4 × undersampling compared to SENSE. CONCLUSIONS: Separate reconstruction of brain and water bath signals enables accelerated MR temperature imaging during MRgFUS procedures with low errors due to undersampling using Cartesian and non-Cartesian trajectories. The spatially-segmented method benefits from multiple coils, and reconstructs temperature with lower error compared to measurements from SENSE-reconstructed images. The acceleration can be applied to increase volumetric coverage and spatiotemporal resolution.

A new method for the determination of aortic pulse wave velocity using cross-correlation on 2D PCMR velocity data.

PURPOSE: To evaluate the reproducibility of a new multisite axial pulse wave velocity (PWV) measurement technique that makes use of 2D PCMR data and cross-correlation analysis. MATERIALS AND METHODS: PWV was estimated with MRI in 13 healthy volunteers by a transit-time technique (TT), a multisite technique utilizing 1D PCMR data in the descending aorta (FOOT), and a new multisite axial technique that uses 2D PCMR data over the ascending, transverse, and descending sections of the aorta (2D-XC). RESULTS: No significant difference was observed between PWV measurements values measured by the three techniques. However, 2D-XC displayed significantly better intertest reproducibility than either the TT or FOOT methodologies. Average percent difference between scans: TT: 15.8% +/- 13.4%, FOOT: 21.3% +/- 16.9%, 2D-XC: 7.72% +/- 4.73%. P = 0.02 for both 2D-XC/TT comparison and 2D-XC/FOOT comparison. CONCLUSION: 2D-XC is a more reproducible method than either the established TT or FOOT methods to estimate the aortic PWV.