Distortion correction of echo planar images applying the concept of finite rate of innovation to point spread function mapping (FRIP)

Objective

Point spread function (PSF) mapping enables estimating the displacement fields required for distortion correction of echo planar images. Recently, a highly accelerated approach was introduced for estimating displacements from the phase slope of under-sampled PSF mapping data. Sampling schemes with varying spacing were proposed requiring stepwise phase unwrapping. To avoid unwrapping errors, an alternative approach applying the concept of finite rate of innovation to PSF mapping (FRIP) is introduced, using a pattern search strategy to locate the PSF peak, and the two methods are compared.

Materials and methods

Fully sampled PSF data was acquired in six subjects at 3.0 T, and distortion maps were estimated after retrospective under-sampling. The two methods were compared for both previously published and newly optimized sampling patterns. Prospectively under-sampled data were also acquired. Shift maps were estimated and deviations relative to the fully sampled reference map were calculated.

Results

The best performance was achieved when using FRIP with a previously proposed sampling scheme. The two methods were comparable for the remaining schemes. The displacement field errors tended to be lower as the number of samples or their spacing increased.

Conclusion

A robust method for estimating the position of the PSF peak has been introduced.

     

An improved PSF mapping method for EPI distortion correction in human brain at ultra high field (7T)

Object  

To further improve the quality and robustness of the point-spread function (PSF) mapping method for fully automatic and accurate correction of geometric distortions in EPI at ultra high field such as 7 Tesla with high fidelity.     

Effects of systematic partial volume errors on the estimation of gray matter cerebral blood flow with arterial spin labeling MRI

Objective

Partial volume (PV) correction is an important step in arterial spin labeling (ASL) MRI that is used to separate perfusion from structural effects when computing the mean gray matter (GM) perfusion. There are three main methods for performing this correction: (1) GM-threshold, which includes only voxels with GM volume above a preset threshold; (2) GM-weighted, which uses voxel-wise GM contribution combined with thresholding; and (3) PVC, which applies a spatial linear regression algorithm to estimate the flow contribution of each tissue at a given voxel. In all cases, GM volume is obtained using PV maps extracted from the segmentation of the T1-weighted (T1w) image. As such, PV maps contain errors due to the difference in readout type and spatial resolution between ASL and T1w images. Here, we estimated these errors and evaluated their effect on the performance of each PV correction method in computing GM cerebral blood flow (CBF).

Materials and methods

Twenty-two volunteers underwent scanning using 2D echo planar imaging (EPI) and 3D spiral ASL. For each PV correction method, GM CBF was computed using PV maps simulated to contain estimated errors due to spatial resolution mismatch and geometric distortions which are caused by the mismatch in readout between ASL and T1w images. Results were analyzed to assess the effect of each error on the estimation of GM CBF from ASL data.

Results

Geometric distortion had the largest effect on the 2D EPI data, whereas the 3D spiral was most affected by the resolution mismatch. The PVC method outperformed the GM-threshold even in the presence of combined errors from resolution mismatch and geometric distortions. The quantitative advantage of PVC was 16% without and 10% with the combined errors for both 2D and 3D ASL. Consistent with theoretical expectations, for error-free PV maps, the PVC method extracted the true GM CBF. In contrast, GM-weighted overestimated GM CBF by 5%, while GM-threshold underestimated it by 16%. The presence of PV map errors decreased the calculated GM CBF for all methods.

Conclusion

The quality of PV maps presents no argument for the preferential use of the GM-threshold method over PVC in the clinical application of ASL.

     

Alleviating artifacts in 1H MRI thermometry by single scan spatiotemporal encoding

Objective

Recent years have seen an increased interest in combining MRI thermometry with devices capable of destroying malignancies by heat ablation. Expected from the MR protocols are accurate and fast thermal characterizations, providing real time feedback on restricted tissue volumes and/or rapidly moving organs like liver. This article explores the potential advantages of relying on spatiotemporally encoded (SPEN) sequences for retrieving real-time thermometric images based on the water’s proton resonance frequency (PRF) shifts.

Materials and methods

Hybrid spatiotemporal/k-space encoding single-scan MRI experiments were implemented on animal and human scanners, and their abilities to deliver single- and multi-slice real-time thermometric measurements based on PRF-derived phase maps in phantoms and in vivo, were compared against echo planar imaging (EPI) and gradient-echo counterparts.

Results

Under comparable acquisition conditions, SPEN exhibited advantages vis-à-vis EPI in terms of dealing with inhomogeneous magnetic field distortions, with shifts arising due to changes in the central frequency offsets, with PRF distributions, and for zooming into restricted fields-of-view without special pulse sequence provisions.

Conclusion

This work confirms the ability of SPEN sequences, particularly when implemented under fully-refocused conditions, to exploit their built-in robustness to shift- and field-derived inhomogeneities for monitoring thermal changes in real-time under in vitro and in vivo conditions.     

A correction method for streak artifacts in gradient-echo EPI using spin-echo EPI reference data

Objective

To analyze the streak artifacts in a gradient-echo echo planar imaging (GE-EPI) sequence and to propose a correction method for the Nyquist ghost artifacts that does not cause streak artifacts in the GE-EPI imaging.

Materials and methods

Several GE-EPI imaging experiments with various reference scans, using both GE-EPI and SE-EPI scan data, were performed to analyze the streak artifacts and to investigate the spin dephasing phenomena of the GE-EPI reference scan. In addition, the analysis based on the spin dephasing was undertaken in order to demonstrate that the SE-EPI reference data can be used for the correction of the GE-EPI main scan data.

Results

The experimental results confirmed that the improvement of the reference data using either signal averaging or a large flip angle cannot guarantee perfect correction of the streak artifact if the noise is not completely removed. Due to the main field inhomogeneity, the spins of the GE-EPI reference data were dephased in multiple echo signals. The proposed correction method, which uses a SE-EPI reference scan for the GE-EPI images, eliminates the N/2 ghost artifacts without producing streak artifacts.

Conclusion

It is believed that the proposed phase error correction scheme can improve the EPI performance in high field MRIs with higher magnetic field inhomogeneities.     

Correction of <Emphasis Type="Italic">B</Emphasis><Subscript>0</Subscript>-induced geometric distortion variations in prospective motion correction for 7T MRI

Objective

Prospective motion correction can effectively fix the imaging volume of interest. For large motion, this can lead to relative motion of coil sensitivities, distortions associated with imaging gradients and B ₀ field variations. This work accounts for the B ₀ field change due to subject movement, and proposes a method for correcting tissue magnetic susceptibility-related distortion in prospective motion correction.

Materials and methods

The B ₀ field shifts at the different head orientations were characterized. A volunteer performed large motion with prospective motion correction enabled. The acquired data were divided into multiple groups according to the object positions. The correction of B ₀-related distortion was applied to each group of data individually via augmented sensitivity encoding with additionally integrated gradient nonlinearity correction.

Results

The relative motion of the gradients, B ₀ field and coil sensitivities in prospective motion correction results in residual spatial distortion, blurring, and coil artifacts. These errors can be mitigated by the proposed method. Moreover, iterative conjugate gradient optimization with regularization provided superior results with smaller RMSE in comparison to standard conjugate gradient.

Conclusion

The combined correction of B ₀-related distortion and gradient nonlinearity leads to a reduction of residual motion artifacts in prospective motion correction data.

     

Characterization of the impact to PET quantification and image quality of an anterior array surface coil for PET/MR imaging

Object

The aim of this study was to determine the impact to PET quantification, image quality and possible diagnostic impact of an anterior surface array used in a combined PET/MR imaging system.

Materials and methods

An extended oval phantom and 15 whole-body FDG PET/CT subjects were re-imaged for one bed position following placement of an anterior array coil at a clinically realistic position. The CT scan, used for PET attenuation correction, did not include the coil. Comparison, including liver SUV_mean, was performed between the coil present and absent images using two methods of PET reconstruction. Due to the time delay between PET scans, a model was used to account for average physiologic time change of SUV.

Results

On phantom data, neglecting the coil caused a mean bias of ?8.2 % for non-TOF/PSF reconstruction, and ?7.3 % with TOF/PSF. On clinical data, the liver SUV neglecting the coil presence fell by ?6.1 % (±6.5 %) for non-TOF/PSF reconstruction; respectively ?5.2 % (±5.3 %) with TOF/PSF. All FDG-avid features seen with TOF/PSF were also seen with non-TOF/PSF reconstruction.

Conclusion

Neglecting coil attenuation for this anterior array coil results in a small but significant reduction in liver SUV_mean but was not found to change the clinical interpretation of the PET images.     

k-t BLAST and SENSE accelerated time-resolved three-dimensional phase contrast MRI in an intracranial aneurysm

Objective

The objective of this study was to investigate the performance of k-t BLAST (Broad-use Linear Acquisition Speed-up Technique) accelerated time-resolved 3D PC-MRI compared to SENSE (SENSitivity Encoding) acceleration in an in vitro and in vivo intracranial aneurysm.

Materials and methods

Non-accelerated, SENSE and k-t BLAST accelerated time-resolved 3D PC-MRI measurements were performed in vivo and in vitro. We analysed the consequences of various temporal resolutions in vitro.

Results

Both in vitro and in vivo measurements showed that the main effect of k-t BLAST was underestimation of velocity during systole. In the phantom, temporal blurring decreased with increasing temporal resolution. Quantification of the differences between the non-accelerated and accelerated measurements confirmed that in systole SENSE performed better than k-t BLAST in terms of mean velocity magnitude. In both in vitro and in vivo measurements, k-t BLAST had higher SNR compared to SENSE. Qualitative comparison between measurements showed good similarity.

Conclusion

Comparison with SENSE revealed temporal blurring effects in k-t BLAST accelerated measurements.     

Lateralization of amygdala activation in fMRI may depend on phase-encoding polarity

Object

Susceptibility artifacts along the phase-encoding (PE) direction impact the activation pattern in the amygdala and may lead to systematic asymmetries. We implemented a triple-echo echo-planar imaging (EPI) sequence, acquiring opposite PE polarities along left?Cright PE direction in a single shot, to investigate its effects on amygdala lateralization.

Materials and Methods

Twelve subjects viewed emotional faces to evoke amygdala activation.

Results and Conclusion

A region of interest analysis revealed that the lateralization of amygdala responses depended on the PE polarity thus representing a pure method artifact. Alternating PE with multi-echo EPI reduced the artifact. Lateralized fMRI activation in areas with magnetic field inhomogeneities need to be interpreted with caution.     

Two-dimensional accelerated MP-RAGE imaging with flexible linear reordering

Object

Implementation of an accelerated Magnetization Prepared RApid Gradient Echo (MP-RAGE) sequence for T1 weighted neuroimaging; exploiting modern MRI technologies to minimize scan time while preserving the image quality.

Materials and methods

A custom MP-RAGE sequence was implemented on a state-of-the-art 3T MR scanner equipped with a 32-channel receiver array head coil. The sequence utilized a shifted CAIPIRINHA k _y –k _z under-sampling pattern combined with elliptical scanning and a two-dimensional view ordering scheme to achieve high parallel imaging acceleration factors at maintained image contrast.

Results

It could be shown that MP-RAGE accelerated in two k-space directions outperforms single direction acceleration, which is the common practice with standard view ordering. Applying the CAIPIRINHA technique in conjunction with elliptical scanning further increased this benefit.

Conclusion

By combining MP-RAGE with CAIPIRINHA sampling and elliptical scanning, the scan time can be reduced from 4–5 min to 2–3 min with insignificant reduction in image quality.     

Application of the limited-memory quasi-Newton algorithm for multi-dimensional,large flip-angle RF pulses at 7T

Objective

Ultrahigh field MRI provides great opportunities for medical diagnostics and research. However, ultrahigh field MRI also brings challenges, such as larger magnetic susceptibility induced field changes. Parallel-transmit radio-frequency pulses can ameliorate these complications while performing advanced tasks in routine applications. To address one class of such pulses, we propose an optimal-control algorithm as a tool for designing advanced multi-dimensional, large flip-angle, radio-frequency pulses. We contrast initial conditions, constraints, and field correction abilities against increasing pulse trajectory acceleration factors.

Materials and methods

On an 8-channel 7T system, we demonstrate the quasi-Newton algorithm with pulse designs for reduced field-of-view imaging with an oil phantom and in vivo with scans of the human brain stem. We used echo-planar imaging with 2D spatial-selective pulses. Pulses are computed sufficiently rapid for routine applications.

Results

Our dataset was quantitatively analyzed with the conventional mean-square-error metric and the structural-similarity index from image processing. Analysis of both full and reduced field-of-view scans benefit from utilizing both complementary measures.

Conclusion

We obtained excellent outer-volume suppression with our proposed method, thus enabling reduced field-of-view imaging using pulse trajectory acceleration factors up to 4.

     

Diffusion-weighted echo planar MR imaging of the neck at 3 T using integrated shimming: comparison of MR sequence techniques for reducing artifacts caused by magnetic-field inhomogeneities

Objective

Our objective was to compare available techniques reducing artifacts in echo planar imaging (EPI)-based diffusion-weighed magnetic resonance imaging MRI (DWI) of the neck at 3 Tesla caused by B0-field inhomogeneities.

Materials and methods

A cylindrical fat–water phantom was equipped with a Maxwell coil allowing for additional linear B0-field variations in z-direction. The effect of increasing strength of this superimposed gradient on image quality was observed using a standard single-shot EPI-based DWI sequence (sEPI), a zoomed single-shot EPI sequence (zEPI), a readout-segmented EPI sequence (rsEPI), and an sEPI sequence with integrated dynamic shimming (intEPI) on a 3-Tesla system. Additionally, ten volunteers were examined over the neck region using these techniques. Image quality was assessed by two radiologists. Scan durations were recorded.

Results

With increasing strength of the external gradient, marked distortions, signal loss, and failure of fat suppression were observed using sEPI, zEPI, and rsEPI. These artifacts were markedly reduced using intEPI. Significantly better in vivo image quality was also observed using intEPI compared with the other techniques. Scan time of intEPI was similar to sEPI and zEPI and shorter than rsEPI.

Conclusion

The use of integrated 2D shim and frequency adjustment for EPI-based DWI results in a significant improvement in image quality of the head/neck region at 3 Tesla. Combining integrated shimming with rsEPI or zEPI can be expected to provide additional improvements.

     

Combined acquisition technique (CAT) for high-field neuroimaging with reduced RF power

Object

Clinical 3 T MRI systems are rapidly increasing and MRI systems with a static field of 7 T or even more have been installed. The RF power deposition is proportional to the square of the static magnetic field strength and is characterized by the specific absorption rate (SAR). Therefore, there exist defined safety limits to avoid heating of the patient. Here, we describe a hybrid method to significantly reduce the SAR compared to a turbo-spin-echo (TSE) sequence.

Materials and methods

We investigate the potential benefits of a combined acquisition technique (CAT) for high-field neuroimaging at 3 and 7 T. The TSE/EPI CAT experiments were performed on volunteers and patients and compared with standard TSE and GRASE protocols. Problems and solutions regarding T2 weighted CAT imaging are discussed.

Results

We present in vivo images with T2 and proton density contrast obtained on 3 and 7 T with significant SAR reduction (up to 60 %) compared with standard TSE. Image quality is comparable to TSE but CAT shows fewer artifacts than a GRASE sequence.

Conclusion

CAT is a promising candidate for neuroimaging at high fields up to 7 T. The SAR reduction allows one to shorten the waiting time between two excitations or to image more slices thereby reducing the overall measurement time.     

An embedded optical tracking system for motion-corrected magnetic resonance imaging at 7T

Object

Prospective motion correction using data from optical tracking systems has been previously shown to reduce motion artifacts in MR imaging of the head. We evaluate a novel optical embedded tracking system.

Materials and methods

The home-built optical embedded tracking system performs image processing within a 7T scanner bore, enabling high speed tracking. Corrected and uncorrected in vivo MR volumes are acquired interleaved using a modified 3D FLASH sequence, and their image quality is assessed and compared.

Results

The latency between motion and correction of the slice position was measured to be (19?±?5)?ms, and the tracking noise has a standard deviation no greater than 10???m/0.005° during conventional MR scanning. Prospective motion correction improved the edge strength by 16?% on average, even though the volunteers were asked to remain motionless during the acquisitions.

Conclusion

Using a novel method for validating the effectiveness of in vivo prospective motion correction, we have demonstrated that prospective motion correction using motion data from the embedded tracking system considerably improved image quality.     

Improving fMRI in signal drop-out regions at 7 T by using tailored radio-frequency pulses: application to the ventral occipito-temporal cortex

Objective

Signal drop-off occurs in echo-planar imaging in inferior brain areas due to field gradients from susceptibility differences between air and tissue. Tailored-RF pulses based on a hyperbolic secant (HS) have been shown to partially recover signal at 3 T, but have not been tested at higher fields.

Materials and methods

The aim of this study was to compare the performance of an optimized tailored-RF gradient-echo echo-planar imaging (TRF GRE-EPI) sequence with standard GRE-EPI at 7 T, in a passive viewing of faces or objects fMRI paradigm in healthy subjects.

Results

Increased temporal-SNR (tSNR) was observed in the middle and inferior temporal lobes and orbitofrontal cortex of all subjects scanned, but elsewhere tSNR decreased relative to the standard acquisition. In the TRF GRE-EPI, increased functional signal was observed in the fusiform, lateral occipital cortex, and occipital pole, regions known to be part of the visual pathway involved in face-object perception.

Conclusion

This work highlights the potential of TRF approaches at 7 T. Paired with a reversed-gradient distortion correction to compensate for in-plane susceptibility gradients, it provides an improved acquisition strategy for future neurocognitive studies at ultra-high field imaging in areas suffering from static magnetic field inhomogeneities.

     

High-resolution 3D whole-heart coronary MRA: a study on the combination of data acquisition in multiple breath-holds and 1D residual respiratory motion compensation

Object

To study a scan protocol for coronary magnetic resonance angiography based on multiple breath-holds featuring 1D motion compensation and to compare the resulting image quality to a navigator-gated free-breathing acquisition. Image reconstruction was performed using L1 regularized iterative SENSE.

Materials and methods

The effects of respiratory motion on the Cartesian sampling scheme were minimized by performing data acquisition in multiple breath-holds. During the scan, repetitive readouts through a k-space center were used to detect and correct the respiratory displacement of the heart by exploiting the self-navigation principle in image reconstruction. In vivo experiments were performed in nine healthy volunteers and the resulting image quality was compared to a navigator-gated reference in terms of vessel length and sharpness.

Results

Acquisition in breath-hold is an effective method to reduce the scan time by more than 30 % compared to the navigator-gated reference. Although an equivalent mean image quality with respect to the reference was achieved with the proposed method, the 1D motion compensation did not work equally well in all cases.

Conclusion

In general, the image quality scaled with the robustness of the motion compensation. Nevertheless, the featured setup provides a positive basis for future extension with more advanced motion compensation methods.     

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms

Objectives

To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging.

Materials and methods

Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared.

Results

Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: ?3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, ?1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6).

Conclusion

The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

     

3D single point imaging with compressed sensing provides high temporal resolution <Emphasis Type="Italic">R</Emphasis><Subscript>2</Subscript>* mapping for in vivo preclinical applications

Objective

Purely phase-encoded techniques such as single point imaging (SPI) are generally unsuitable for in vivo imaging due to lengthy acquisition times. Reconstruction of highly undersampled data using compressed sensing allows SPI data to be quickly obtained from animal models, enabling applications in preclinical cellular and molecular imaging.

Materials and methods

TurboSPI is a multi-echo single point technique that acquires hundreds of images with microsecond spacing, enabling high temporal resolution relaxometry of large-R ₂* systems such as iron-loaded cells. TurboSPI acquisitions can be pseudo-randomly undersampled in all three dimensions to increase artifact incoherence, and can provide prior information to improve reconstruction. We evaluated the performance of CS-TurboSPI in phantoms, a rat ex vivo, and a mouse in vivo.

Results

An algorithm for iterative reconstruction of TurboSPI relaxometry time courses does not affect image quality or R ₂* mapping in vitro at acceleration factors up to 10. Imaging ex vivo is possible at similar acceleration factors, and in vivo imaging is demonstrated at an acceleration factor of 8, such that acquisition time is under 1 h.

Conclusions

Accelerated TurboSPI enables preclinical R ₂* mapping without loss of data quality, and may show increased specificity to iron oxide compared to other sequences.

     

In vivo visualization of cells labeled with superparamagnetic iron oxides by a sub-millisecond gradient echo sequence

Object

In vivo magnetic resonance imaging (MRI) of iron-labeled pancreatic islets (PIs) transplanted into the liver is still challenging in humans. The aim of this study was to develop and evaluate a double contrast method for the detection of PIs labeled with superparamagnetic iron oxide (SPIO) nanoparticles.

Materials and methods

A double-echo three-dimensional (3D) spoiled gradient echo sequence was adapted to yield a sub-millisecond first echo time using variable echo times and highly asymmetric Cartesian readout. Positive contrast was achieved by conventional and relative image subtraction. Experiments for cell detection efficiency were performed in vitro on gelatin phantoms, in vivo on a Lewis rat and on a patient 6 months after PI transplantation.

Results

It was demonstrated that the proposed method can be used for the detection of transplanted PIs with positive contrast in vitro and in vivo. For all experiments, relative subtraction yielded comparable and in some cases better contrast than conventional subtraction. For the first time, positive contrast imaging of transplanted human PIs was performed in vivo in patients.

Conclusion

The proposed method allows 3D data acquisition within a single breath-hold and yields enhanced contrast-to-noise ratios of transplanted SPIO labeled pancreatic islets relative to negative contrast images, therefore providing improved identification.