Fast multi-planar gradient echo MR imaging: impact of variation in pulse sequence parameters on image quality and artifacts

The purpose of this investigation was to quantitatively evaluate the practical impact of alteration of key imaging parameters on image quality and artifacts in fast multi-planar gradient echo (GRE) pulse sequences. These include multi-planar GRASS (MPGR) and fast multi-planar spoiled GRASS (FMPSPGR). We developed a composite phantom with different T(1) and T(2) values comprising the range of common biological tissues, which was also subjected to periodic motion in order to evaluate motion effects. Magnetic resonance imaging was performed on a GE Signa 1.5-T system. Experimental variations in key parameters included excitation flip angle (FL), echo time (TE), repetition time (TR), and receive bandwidth (BW). Quantitative analysis consisted of signal-to-noise-ratio (SNR) and contrast (CN), image nonuniformity (NU), full-width-at-half-maximum (FWHM) (i.e., blurring or geometric distortion), and ghosting ratio (GR). We found that flip angle, TE, and TR play particularly critical roles in determining image signal, homogeneity, and ghosting artifact with these sequences. Optimum clinical application of these pulse sequences requires careful attention to these imaging parameters and to their complex interactions.     

An optimized pulse sequence for isotropically weighted diffusion imaging.

Single-shot echo-planar imaging is becoming the most widely used technique for magnetic resonance diffusion imaging, since it enables measurement of diffusion coefficients in human brain without motion artifacts. However, its reliability is limited by geometrical distortions due to eddy currents. In this report, an isotropically weighted echo-planar pulse sequence, optimized to give the maximum signal-to-noise ratio in the computed trace image and designed to produce inherently low distortions, is presented. It is also shown how the residual translational distortion can be easily characterized and removed by postprocessing. A full characterization of the distortion artifact involves a few measurements on a phantom, in order to estimate the distortion as a function of slice orientation, which can then be used to correct any slice orientation. Results of applying the image translation correction to data collected from a patient are presented.     

On the use of water phantom images to calibrate and correct eddy current induced artefacts in MR diffusion tensor imaging

The accurate determination of absolute measures of diffusion anisotropy in vivo using single-shot, echo-planar imaging techniques requires the acquisition of a set of high signal-to-noise ratio, diffusion-weighted images that are free from eddy current induced image distortions. Such geometric distortions can be characterized and corrected in brain imaging data using magnification (M), translation (T), and shear (S) distortion parameters derived from separate water phantom calibration experiments. Here we examine the practicalities of using separate phantom calibration data to correct high b-value diffusion tensor imaging data by investigating the stability of these distortion parameters, and hence the eddy currents, with time. It is found that M, T, and S vary only slowly with time (i.e., on the order of weeks), so that calibration scans need not be performed after every patient examination. This not only minimises the scan time required to collect the calibration data, but also the computational time needed to characterize these eddy current induced distortions. Examples of how measurements of diffusion anisotropy are improved using this post-processing scheme are also presented.     

Quantitative MRI: Defining repeatability,reproducibility and accuracy for prostate cancer imaging biomarker development

IntroductionQuantitative MRI (qMRI) parameters have been increasingly used to develop predictive models to accurately monitor treatment response in prostate cancer after radiotherapy. To reliably detect changes in signal due to treatment response, predictive models require qMRI parameters with high repeatability and reproducibility. The purpose of this study was to measure qMRI parameter uncertainties in both commercial and in-house developed phantoms to guide the development of robust predictive models for monitoring treatment response.Materials and methodsADC, T1, and R2* values were acquired across three 3 T scanners with a prostate-specific qMRI protocol using the NIST/ISMRM system phantom, RSNA/NIST diffusion phantom, and an in-house phantom. A B1 field map was acquired to correct for flip angle inhomogeneity in T1 maps. All sequences were repeated in each scan to assess within-session repeatability. Weekly scans were acquired on one scanner for three months with the in-house phantom. Between-session repeatability was measured with test-retest scans 6-months apart on all scanners with all phantoms. Accuracy, defined as percentage deviation from reference value for ADC and T1, was evaluated using the system and diffusion phantoms. Repeatability and reproducibility coefficients of variation (%CV) were calculated for all qMRI parameters on all phantoms.ResultsOverall, repeatability CV of ADC was <2.40%, reproducibility CV was <3.98%, and accuracy ranged between −8.0% to 2.7% across all scanners. Applying B1 correction on T1 measurements significantly improved the repeatability and reproducibility (p<0.05) but increased error in accuracy (p<0.001). Repeatability and reproducibility of R2* was <4.5% and <7.3% respectively in the system phantom across all scanners.ConclusionRepeatability, reproducibility, and accuracy in qMRI parameters from a prostate-specific protocol was estimated using both commercial and in-house phantoms. Results from this work will be used to identify robust qMRI parameters for use in the development of predictive models to longitudinally monitor treatment response for prostate cancer in current and future clinical trials.     

Inter-vender and test-retest reliabilities of resting-state functional magnetic resonance imaging: Implications for multi-center imaging studies

This prospective multi-center study aimed to evaluate the inter-vendor and test-retest reliabilities of resting-state functional magnetic resonance imaging (RS-fMRI) by assessing the temporal signal-to-noise ratio (tSNR) and functional connectivity. Study included 10 healthy subjects and each subject was scanned using three 3 T MR scanners (GE Signa HDxt, Siemens Skyra, and Philips Achieva) in two sessions. The tSNR was calculated from the time course data. Inter-vendor and test-retest reliabilities were assessed with intra-class correlation coefficients (ICCs) derived from variant component analysis. Independent component analysis was performed to identify the connectivity of the default-mode network (DMN). In result, the tSNR for the DMN was not significantly different among the GE, Philips, and Siemens scanners (P = 0.638). In terms of vendor differences, the inter-vendor reliability was good (ICC = 0.774). Regarding the test-retest reliability, the GE scanner showed excellent correlation (ICC = 0.961), while the Philips (ICC = 0.671) and Siemens (ICC = 0.726) scanners showed relatively good correlation. The DMN pattern of the subjects between the two sessions for each scanner and between three scanners showed the identical patterns of functional connectivity. The inter-vendor and test-retest reliabilities of RS-fMRI using different 3 T MR scanners are good. Thus, we suggest that RS-fMRI could be used in multicenter imaging studies as a reliable imaging marker.     

Use of USPIO-induced magnetic susceptibility artifacts to identity sentinel lymph nodes and lymphatic drainage patterns. I. dependence of artifact size with subcutaneous Combidex® dose in rats

Subcutaneously administered Combidex^® contrast agent produced characteristic magnetic susceptibility artifacts in gradient-echo (GE) images of rat brachial and axillary lymph nodes. These artifacts were useful in the rapid location and identification of normal sentinel lymph nodes. A linear dose response was observed with maximum artifact size in transverse images and was used noninvasively to study lymphatic drainage patterns.     

Fast T2-weighted MR imaging: impact of variation in pulse sequence parameters on image quality and artifacts

The purpose of this study was to quantitatively evaluate in a phantom model the practical impact of alteration of key imaging parameters on image quality and artifacts for the most commonly used fast T(2)-weighted MR sequences. These include fast spin-echo (FSE), single shot fast spin-echo (SSFSE), and spin-echo echo-planar imaging (EPI) pulse sequences. We developed a composite phantom with different T1 and T2 values, which was evaluated while stationary as well as during periodic motion. Experiments involved controlled variations in key parameters including effective TE, TR, echo spacing (ESP), receive bandwidth (BW), echo train length (ETL), and shot number (SN). Quantitative analysis consisted of signal-to-noise ratio (SNR), image nonuniformity, full-width-at-half-maximum (i.e., blurring or geometric distortion) and ghosting ratio. Among the fast T(2)-weighted sequences, EPI was most sensitive to alterations in imaging parameters. Among imaging parameters that we tested, effective TE, ETL, and shot number most prominently affected image quality and artifacts. Short T(2) objects were more sensitive to alterations in imaging parameters in terms of image quality and artifacts. Optimal clinical application of these fast T(2)-weighted imaging pulse sequences requires careful attention to selection of imaging parameters.     

Considerations in applying 3D PRESS H-1 brain MRSI with an eight-channel phased-array coil at 3 T

The purpose of this study was to assess the benefits of a 3 T scanner and an eight-channel phased-array head coil for acquiring three-dimensional PRESS (Point REsolved Spectral Selection) proton (H-1) magnetic resonance spectroscopic imaging (MRSI) data from the brains of volunteers and patients with brain tumors relative to previous studies that used a 1.5 T scanner and a quadrature head coil. Issues that were of concern included differences in chemical shift artifacts, line broadening due to increased susceptibility at higher field strengths, changes in relaxation times and the increased complexity of the postprocessing software due to the need for combining signals from the multichannel data. Simulated and phantom spectra showed that very selective suppression pulses with a thickness of 40 mm and an overpress factor of at least 1.2 are needed to reduce chemical shift artifact and lipid contamination at higher field strengths. Spectral data from a phantom and those from six volunteers demonstrated that the signal-to-noise ratio (SNR) in the eight-channel coil was more than 50% higher than that in the quadrature head coil. For healthy volunteers and eight patients with brain tumors, the SNR at 3 T with the eight-channel coil was on average 1.5 times higher relative to the eight-channel coil at 1.5 T in voxels from normal-appearing brains. In combination with the effect of a higher field strength, the use of the eight-channel coil was able to provide an increase in the SNR of more than 2.33 times the corresponding acquisition at 1.5 T with a quadrature head coil. This is expected to be critical for clinical applications of MRSI in patients with brain tumors because it can be used to either decrease acquisition time or improve spatial resolution.     

Reduced aliasing artifacts using shaking projection k-space sampling trajectory

Radial imaging techniques, such as projection-reconstruction （PR）, are used in magnetic resonance imaging （MRI） for dynamic imaging, angiography, and short-T2 imaging. They are less sensitive to flow and motion artifacts, and support fast imaging with short echo times. However, aliasing and streaking artifacts are two main sources which degrade radial imaging quality. For a given fixed number of k-space projections, data distributions along radial and angular directions will influence the level of aliasing and streaking artifacts. Conventional radial k-space sampling trajectory introduces an aliasing artifact at the first principal ring of point spread function （PSF）. In this paper, a shaking projection （SP） k-space sampling trajectory was proposed to reduce aliasing artifacts in MR images. SP sampling trajectory shifts the projection alternately along the k-space center, which separates k-space data in the azimuthal direction. Simulations based on conventional and SP sampling trajectories were compared with the same number projections. A significant reduction of aliasing artifacts was observed using the SP sampling trajectory. These two trajectories were also compared with different sampling frequencies. ASP trajectory has the same aliasing character when using half sampling frequency （or half data） for reconstruction. SNR comparisons with different white noise levels show that these two trajectories have the same SNR character. In conclusion, the SP trajectory can reduce the aliasing artifact without decreasing SNR and also provide a way for undersampling recon- struction. Furthermore, this method can be applied to three-dimensional （3D） hybrid or spherical radial k-space sampling for a more efficient reduction of aliasing artifacts.     

Contrast-enhanced MR angiography with frequency-dependent mask subtraction

One of the main challenges for high-resolution contrast-enhanced magnetic resonance angiography (CE-MRA) is the limited signal-to-noise ratio (SNR). In conventional CE-MRA, spoiled gradient-recalled echo sequence is typically used to acquire raw data before and after contrast material injection, followed by mask subtraction to improve vessel contrast at the cost of increased image noise. We reported here on a frequency-dependent mask subtraction technique where only low spatial frequency data were subject to mask subtraction to improve image contrast and reduce motion artifact while high spatial frequency data were not subtracted to preserve image noise level. The feasibility of this technique was demonstrated through phantom, animal, and human volunteer studies. The lowest half, quarter, 1/8 and 1/16 of all spatial frequencies were subject to mask subtraction, respectively. These partial subtraction techniques were compared with conventional mask subtraction in terms of SNR and artifact power. An SNR gain up to 37% was achieved without significant artifact power increase for quarter-subtraction where only the central one quarter k-space data were subject to mask subtraction.     

Ficoll as testing material for diffusion weighted imaging-quality assurance phantoms

PurposeThe aim of this work is to test the use of aqueous solutions of Ficoll®**, a highly branched polymer displaying crowding properties, to build a phantom suitable for Diffusion Weighted Imaging (DWI) in Magnetic Resonance Imaging (MRI).MethodsWe developed a test object made of a cylindrical plastic container with a precise geometrical arrangement suitable for measuring several samples at the same time. The container was designed to host single vials with variable geometry and number, and to fit inside common commercial head coils for MRI scanners.In our experiments, vials were filled with 8 aqueous solutions of Ficoll 70 and Ficoll 400 spanning a range of polymer concentration from 5 to 30% by weight. Vials containing ultra-pure water were also used as reference. Experiments were performed on both 1.5 and 3 T clinical scanners (GE, Philips and Siemens), under the conditions of a standard clinical examination.ResultsThe geometry of the phantom provided reduced imaging artifacts, especially image distortions at magnetic interfaces. We found that the Apparent Diffusion Coefficient (ADC) varied in the range of 0.00125–0.00223 mm²/s and decreased with Ficoll concentration. ADC vs Ficoll concentration exhibited a linear trend. Results were consistent over time and among different MRI clinical scanners, showing an average variability of 3% at 1.5 T and of 7.5% at 3 T. Moreover, no substantial difference was found between Ficoll 70 and 400. By varying Ficoll concentration, ADC can be modulated to approach tissue-mimicking values. Preliminary results for relaxation measurements proved that both T₁ and T₂ decreased with Ficoll concentration in the ranges 1.3–2.4 s and 150–800 ms respectively.ConclusionsIn this work, we propose a 3D phantom design based on the widespread crowding agent Ficoll, which is suitable for DWI quality assurance purposes in MRI acquisitions. Aqueous Ficoll solutions provide good performance in terms of stability, ease of preparation, and safety.     

Reduction of Artifacts in Capacitive Electrocardiogram Signals of Driving Subjects

The development of smart cars with e-health services allows monitoring of the health condition of the driver. Driver comfort is preserved by the use of capacitive electrodes, but the recorded signal is characterized by large artifacts. This paper proposes a method for reducing artifacts from the ECG signal recorded by capacitive electrodes (cECG) in moving subjects. Two dominant artifact types are coarse and slow-changing artifacts. Slow-changing artifacts removal by classical filtering is not feasible as the spectral bands of artifacts and cECG overlap, mostly in the band from 0.5 to 15 Hz. We developed a method for artifact removal, based on estimating the fluctuation around linear trend, for both artifact types, including a condition for determining the presence of coarse artifacts. The method was validated on cECG recorded while driving, with the artifacts predominantly due to the movements, as well as on cECG recorded while lying, where the movements were performed according to a predefined protocol. The proposed method eliminates 96% to 100% of the coarse artifacts, while the slow-changing artifacts are completely reduced for the recorded cECG signals larger than 0.3 V. The obtained results are in accordance with the opinion of medical experts. The method is intended for reliable extraction of cardiovascular parameters to monitor driver fatigue status.     

Clinical comparison of single-shot EPI,readout-segmented EPI and TGSE-BLADE for diffusion-weighted imaging of cerebellopontine angle tumors on 3 tesla

ObjectiveThe complex anatomical structures of cerebellopontine angle (CPA) pose a unique challenge to diffusion weighted imaging (DWI). This study aimed to compare the clinical utility of the prototypic 2D turbo gradient- and spin echo-BLADE-DWI (TGSE-BLADE-DWI) with that of readout-segmented echo-planar DWI (RESOLVE-DWI) and single-shot echo-planar DWI (SS-EPI-DWI) to visualize CPA anatomic structures and identify CPA tumors.MethodsA total of 8 volunteers and 36 patients with pathological CPA tumors were enrolled to perform the three DWI sequences at 3 T. Scan time of TGSE-BLADE-DWI, RESOLVE-DWI and SS-EPI-DWI was 5 min 51 s, 5 min 15 s and 1 min 22 s, respectively. Subjective analysis, including visualization of anatomical structures, geometric distortion, ghosting artifacts, lesion conspicuity, diagnostic confidence, and overall image quality of the three DWI sequences were scored and assessed. Signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and apparent diffusion coefficient (ADC) of CPA tumors were measured and compared.ResultsA total of 39 lesions were identified, TGSE-BLADE-DWI detected all of them, RESOLVE-DWI 36 and SS-EPI-DWI 27. Significant differences were found in all the subjective parameters among the three DWI sequences (all p < 0.001). TGSE-BLADE-DWI was significantly better than RESOLVE-DWI in visualization of CPA anatomical structures, geometric distortion, ghosting artifacts, lesion conspicuity, diagnostic confidence, and overall image quality (all p < 0.01), and RESOLVE-DWI showed significantly superior performance than SS-EPI-DWI in all parameters (all p < 0.001). CNRs and ADCs were not significantly different among the three DWI sequences (p = 0.355, p = 0.590, respectively). No significant differences were detected between TGSE-BLADE-DWI SNR and RESOLVE-DWI SNR (p = 0.058), or TGSE-BLADE-DWI SNR and SS-EPI-DWI SNR (p = 0.155).ConclusionCompared with RESOLVE-DWI and SS-EPI-DWI, TGSE-BLADE-DWI minimized geometric distortions and ghosting artifacts and demonstrated an improved ability for depicting CPA tumors with better lesion conspicuity.SummaryGeometric distortions and ghosting artifacts are found at bone-air interfaces using conventional diffusion-weighted imaging (DWI), which is a challenge for imaging cerebellopontine angle (CPA) tumors. Our study validated that geometric distortions and ghosting artifacts were not present on 2D turbo gradient- and spin-echo-BLADE-DWI scans, making this technique useful for visualizing CPA anatomic structures and diagnosing CPA tumors.     

Optimization of spoiler gradients in flash MRI

The FLASH technique for fast magnetic resonance (MR) imaging often employs strong magnetic field gradients, called spoiler gradients, to dephase the transverse magnetization after it has been measured. Otherwise, image artifacts can develop. The effectiveness of spoiler gradients at suppressing these artifacts was evaluated experimentally on two-dimensional MR images of a uniform phantom and patients. It was informative to compare the magnetization immediately before the RF excitation in each phase encoding step. Only spoiler gradients in the slice selection direction were effective. Spoiler gradients that decreased steadily from a large amplitude in the first phase encoding step to zero in the last minimized the transverse magnetization and suppressed the image artifact, without changing the image contrast.     

DC artifact correction for arbitrary phase-cycling sequence

In magnetic resonance imaging (MRI), a non-zero offset in the receiver baseline signal during acquisition results in a bright spot or a line artifact in the center of the image known as a direct current (DC) artifact. Several methods have been suggested in the past for the removal or correction of DC artifacts in MR images, however, these methods cannot be applied directly when a specific phase-cycling technique is used in the imaging sequence. In this work, we proposed a new, simple technique that enables correction of DC artifacts for any arbitrary phase-cycling imaging sequences. The technique is composed of phase unification, DC offset estimation and correction, and phase restoration. The feasibility of the proposed method was demonstrated via phantom and in vivo experiments with a multiple phase-cycling balanced steady-state free precession (bSSFP) imaging sequence. Results showed successful removal of the DC artifacts in images acquired using bSSFP with phase-cycling angles of 0°, 90°, 180°, and 270°, indicating potential feasibility of the proposed method to any imaging sequence with arbitrary phase-cycling angles.     

Mapping eddy current induced fields for the correction of diffusion-weighted echo planar images

Diffusion-weighted echo-planar magnetic resonance imaging is potentially of great importance as a diagnostic imaging tool; however, the technique currently suffers a number of limitations, including the image distortion caused by the eddy current induced fields when the diffusion-weighting magnetic field gradient pulses are applied. The distortions cause mis-registration between images with different diffusion-weighting, that then results in artifacts in quantitative diffusion images. A method is presented to measure the magnetic fields generated from the eddy currents for each of three orthogonal gradient pulse vectors, and then to use these to ascertain the image distortion that occurs in subsequent diffusion-weighted images with arbitrary gradient pulse vector amplitude and direction, and image plane orientation. The image distortion can then be reversed. Both temporal and spatial dependence of the residual eddy current induced fields are included in the analysis. Image distortion was substantially reduced by the correction scheme, for arbitrary slice position and angulation. This method of correction is unaffected by the changes in image contrast that occur due to diffusion weighting, and does not need any additional scanning time during the patient scan. It is particularly suitable for use with single-shot echo planar imaging.     

Semiautomated quality assurance for quantitative magnetic resonance imaging.

It is now well established that MRI can be used for quantitative (as opposed to simply qualitative) measurements, and good accuracy and precision have been obtained in phantom experiments. To make routine quantitative measurements as part of a clinical scanning protocol, however, quality assurance (QA) methods particularly suited to quantification must be developed. We describe a set of QA tests using clinical protocols on test phantoms, with which we have assessed quantitative performance of our Picker 0.5-T scanner (Picker International, Cleveland, OH) over 2 years. We also describe the automated data processing methods we have developed to deal with the large amounts of data generated by these tests.     

Localized high-resolution DTI of the human midbrain using single-shot EPI,parallel imaging,and outer-volume suppression at 7 T

Localized high-resolution diffusion tensor images (DTI) from the midbrain were obtained using reduced field-of-view (rFOV) methods combined with SENSE parallel imaging and single-shot echo planar (EPI) acquisitions at 7 T. This combination aimed to diminish sensitivities of DTI to motion, susceptibility variations, and EPI artifacts at ultra-high field. Outer-volume suppression (OVS) was applied in DTI acquisitions at 2- and 1-mm² resolutions, b = 1000 s/mm², and six diffusion directions, resulting in scans of 7- and 14-min durations. Mean apparent diffusion coefficient (ADC) and fractional anisotropy (FA) values were measured in various fiber tract locations at the two resolutions and compared. Geometric distortion and signal-to-noise ratio (SNR) were additionally measured and compared for reduced-FOV and full-FOV DTI scans. Up to an eight-fold data reduction was achieved using DTI-OVS with SENSE at 1 mm², and geometric distortion was halved. The localization of fiber tracts was improved, enabling targeted FA and ADC measurements. Significant differences in diffusion properties were observed between resolutions for a number of regions suggesting that FA values are impacted by partial volume effects even at a 2-mm² resolution. The combined SENSE DTI-OVS approach allows large reductions in DTI data acquisition and provides improved quality for high-resolution diffusion studies of the human brain.     

Geometric distortion in clinical MRI systems Part I: evaluation using a 3D phantom

Recently, a 3D phantom that can provide a comprehensive and accurate measurement of the geometric distortion in MRI has been developed. Using this phantom, a full assessment of the geometric distortion in a number of clinical MRI systems (GE and Siemens) has been carried out and detailed results are presented in this paper. As expected, the main source of geometric distortion in modern superconducting MRI systems arises from the gradient field nonlinearity. Significantly large distortions with maximum absolute geometric errors ranged between 10 and 25 mm within a volume of 240 x 240 x 240 mm(3) were observed when imaging with the new generation of gradient systems that employs shorter coils. By comparison, the geometric distortion was much less in the older-generation gradient systems. With the vendor's correction method, the geometric distortion measured was significantly reduced but only within the plane in which these 2D correction methods were applied. Distortion along the axis normal to the plane was, as expected, virtually unchanged. Two-dimensional correction methods are a convenient approach and in principle they are the only methods that can be applied to correct geometric distortion in a single slice or in multiple noncontiguous slices. However, these methods only provide an incomplete solution to the problem and their value can be significantly reduced if the distortion along the normal of the correction plane is not small.     

Resolution recovery in Turbo Spin Echo using segmented Half Fourier acquisition

Turbo Spin Echo (TSE) is a sequence of choice for obtaining T(2)-weighted images. TSE reduces acquisition time by acquiring several echoes within each TR, at the cost of introducing an exponential weighting in the k-space that leads to a certain image blurring. This is particularly important for short-T(2) structures, which can even disappear if their size in the phase encoding direction is comparable to the degree of blurring. This article suggests the use of a combination of Half Fourier (HF) and segmented (multishot) TSE (sHF-TSE) to recover the original resolution of the SE images. The improved symmetry of the dataset achieved by HF reconstruction is used to increase the resolution of the TSE images. The proposed combination, available in most clinical scanners, reduces the blurring artifact inherent to the TSE sequence without increasing the scan time or the number of acquisitions, but at the cost of a slight reduction of the signal-to-noise ratios (SNR). Qualitative and quantitative results are presented using both numerical simulation and imaging. Significant edge enhancement has been achieved for structures with short T(2), (narrowing of the full width at half maximum [FWHM] up to 45%). The proposed sequence is more sensitive to movement artifacts but has proven to be superior to the conventional TSE for imaging static structures.