Readout Gradient Research Articles

Seed point discontinuity-based segmentation method for the substantia nigra and the red nucleus in quantitative susceptibility maps.

The automatic segmentation of cerebral nuclei in the quantitative susceptibility mapping (QSM) images can provide assistance for surgical treatment and pathological mechanism studies. However, as the most frequently used segmentation method, the atlas method provides unsatisfactory results when segmenting the substantia nigra (SN) and the red nucleus (RN). To propose and evaluate an improved automatic method based on seed points-discontinuity for segmentations of the SN and the RN in QSM images. Prospective. In all, 22 subjects, 11 patients with Parkinson's disease (PD), and 11 healthy subjects (mean age of 68.0 ± 6.9 years) underwent MR scans. 3T system and a 3D multiecho gradient echo sequence with monopolar readout gradient. Manual segmentations by two radiologists (both with over 10 years of experience in neuroimaging) were used to establish a baseline for assessment. The Dice coefficient and the center-of-gravity distance was employed to evaluate the segmentation accuracy. The mean value and standard deviation of the Dice coefficient and center-of-gravity distance were calculated separately to compare segmentation results from the proposed method, the level set method, the atlas method (including the single-atlas method and the multi-atlas majority voting method). The statistical results of Dice coefficient of the SN and the RN between the ground truth and the segmentation were 0.79 ± 0.14 and 0.77 ± 0.06 for the proposed method, 0.40 ± 0.10 and 0.65 ± 0.09 for the level set method, 0.68 ± 0.09 and 0.64 ± 0.07 for the single-atlas method, 0.70 ± 0.06 and 0.68 ± 0.05 for the multi-atlas majority voting method, respectively. The proposed method also provides the lowest center-of-gravity distance value (1.05 ± 0.71 for the SN and 0.74 ± 0.35 for the RN). The segmentation results of the proposed method performed well on the in vivo data and were closer to the manual segmentation than the atlas method. 1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;48:1112-1119.

Increasing robustness of radial GRASE acquisition for SAR-reduced brain imaging

PurposeTo improve a radial multi-slice 2D gradient- and spin-echo (GRASE) sequence and provide an appropriate image reconstruction technique for SAR-reduced high-resolution neuroimaging. MethodsAdditional readout gradients per radio-frequency (RF) refocusing allow for a reduced number of RF pulses. In this way, a specific absorption rate (SAR) reduction is achieved and the application at high-field systems becomes more feasible. A phase insensitive image reconstruction is proposed to reduce signal dropout artifacts originating from opposite readout polarities. In addition, the image reconstruction allows for the calculation of images with varying contrast from one measurement. ResultsResults obtained at 3T and 7T demonstrate a SAR-reduction of at least 66% for a single-slice experiment with radial GRASE. The reduced SAR is used for an increased spatial coverage without increasing the measurement time. Experiments at 3T and 7T showed that the visual image quality is comparable to standard TSE and GRASE sequences with the same measurement parameters. Using higher EPI factors and the presented image reconstruction, artifact-free images with a significant SAR-reduction can be achieved. ConclusionRadial GRASE enables SAR-reduced acquisitions of high-resolution brain images with different contrasts from one measurement and is a promising sequence for high-field neuroimaging.

Phosphorus magnetic resonance spectroscopic imaging using flyback echo planar readout trajectories.

To present and evaluate a fast phosphorus magnetic resonance spectroscopic imaging (MRSI) sequence using echo planar spectroscopic imaging with flyback readout gradient trajectories. Waveforms were designed and implemented using a 3 Tesla MRI system. 31P spectra were acquired with 2×2cm2 and 3×3cm2 resolution over a 20- and 21-cm field of view and spectral bandwidths up to 1923Hz. The sequence was first tested using a 20-cm-diameter phosphate phantom, and subsequent in vivo tests were performed on healthy human calf muscles and brains from five volunteers. Flyback EPSI achieved 10× and 7× reductions in acquisition time, with 68.0±1.2 and 69.8±2.2% signal-to-noise ratio (SNR) per unit of time efficiency (theoretical SNR efficiency was 74.5 and 76.4%) for the in vivo experiments, compared to conventional phase-encoded MRSI for the 2×2cm2 and 3×3cm2 resolution waveforms, respectively. Statistical analysis showed no difference in the quantification of most metabolites. Time savings and SNR comparisons were consistent across phantom, leg and brain experiments. EPSI using flyback readout trajectories was found to be a reliable alternative for acquiring 31P-MRSI data in a shorter acquisition time.

Improved three-dimensional multi-echo gradient echo based myelin water fraction mapping with phase related artifact correction

Myelin water fraction (MWF) mapping with a multi-echo gradient echo (mGRE) sequence using complex-value based model fitting approach was recently described, in which the phase of the data plays an important role in the accuracy of the fitting results. Thus, influences coming from the undesirable phase components need to be reduced. Targeted for improved MWF mapping, methods to combat these phase related issues which include offset correction, main magnetic field (B0) inhomogeneity related correction, and flow compensation have been developed. For offset correction, a coil combine method with bipolar readout gradients has been devised. For B0 related artifact, a voxel spread function correction approach along with a navigator echo acquisition was utilized. Finally, flow compensation gradients were used to reduce the effects of flow from regions including the veins. The correction methods led to reduced residual fitting error and improved quality of the resultant MWF maps. Quantitative analysis shows increased reliability when the corrections were used. Data from normal volunteers showed improved MWF mapping with the integrated method developed in this work.

Inner‐volume echo volumar imaging (IVEVI) for robust fetal brain imaging

Fetal functional MRI studies using conventional 2-dimensional single-shot echo-planar imaging sequences may require discarding a large data fraction as a result of fetal and maternal motion. Increasing the temporal resolution using echo volumar imaging (EVI) could provide an effective alternative strategy. Echo volumar imaging was combined with inner volume (IV) imaging (IVEVI) to locally excite the fetal brain and acquire full 3-dimensional images, fast enough to freeze most fetal head motion. IVEVI was implemented by modifying a standard multi-echo echo-planar imaging sequence. A spin echo with orthogonal excitation and refocusing ensured localized excitation. To introduce T2* weighting and to save time, the k-space center was shifted relative to the spin echo. Both single and multi-shot variants were tested. Acoustic noise was controlled by adjusting the amplitude and switching frequency of the readout gradient. Image-based shimming was used to minimize B0 inhomogeneities within the fetal brain. The sequence was first validated in an adult. Eight fetuses were scanned using single-shot IVEVI at a 3.5 × 3.5 × 5.0 mm3 resolution with a readout duration of 383 ms. Multishot IVEVI showed reduced geometric distortions along the second phase-encode direction. Fetal EVI remains challenging. Although effective echo times comparable to the T2* values of fetal cortical gray matter at 3 T could be achieved, controlling acoustic noise required longer readouts, leading to substantial distortions in single-shot images. Although multishot variants enabled us to reduce susceptibility-induced geometric distortions, sensitivity to motion was increased. Future studies should therefore focus on improvements to multishot variants. Magn Reson Med 80:279-285, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Multi-Rate Acquisition for Dead Time Reduction in Magnetic Resonance Receivers: Application to Imaging With Zero Echo Time.

For magnetic resonance imaging of tissues with very short transverse relaxation times, radio-frequency excitation must be immediately followed by data acquisition with fast spatial encoding. In zero-echo-time (ZTE) imaging, excitation is performed while the readout gradient is already on, causing data loss due to an initial dead time. One major dead time contribution is the settling time of the filters involved in signal down-conversion. In this paper, a multi-rate acquisition scheme is proposed to minimize dead time due to filtering. Short filters and high output bandwidth are used initially to minimize settling time. With increasing time since the signal onset, longer filters with better frequency selectivity enable stronger signal decimation. In this way, significant dead time reduction is accomplished at only a slight increase in the overall amount of output data. Multi-rate acquisition was implemented with a two-stage filter cascade in a digital receiver based on a field-programmable gate array. In ZTE imaging in a phantom and in vivo, dead time reduction by multi-rate acquisition is shown to improve image quality and expand the feasible bandwidth while increasing the amount of data collected by only a few percent.

Echo-planar spectroscopic imaging with dual-readout alternated gradients (DRAG-EPSI) at 7 T: Application for 2-hydroxyglutarate imaging in glioma patients.

To develop echo-planar spectroscopic imaging (EPSI) with large spectral width and accomplish high-resolution imaging of 2-hydroxyglutarate (2HG) at 7 T. We designed a new EPSI readout scheme at 7 T. Data were recorded with dual-readout alternated gradients and combined according to the gradient polarity. Following validation of its performance in phantoms, the new readout scheme, together with previously reported 2HG-optimized magnetic resonance spectroscopy (point-resolved spectroscopy echo time of 78 ms), was used for time-efficient and high-resolution imaging of 2HG and other metabolites in five glioma patients before treatment. Unsuppressed water, acquired with EPSI, was used as reference for multichannel combination, eddy-current compensation, and metabolite quantification. Spectral fitting was conducted with the LCModel using in-house calculated basis sets. Using a readout gradient strength of 9.5 mT/m and slew rate of 90 mT/m/ms, dual-readout alternated gradients EPSI permitted 1638-Hz spectral width with 6 × 6 mm2 in-plane resolution at 7 T. Phantom data indicated that dual-readout alternated gradients EPSI provides proper metabolite signals and induces much less frequency drifts than conventional EPSI. For a spatial resolution of 0.5 mL, 2HG was detected in tumors with precision (Cramer-Rao lower bound < 10%). The 2HG was estimated to be 2.3 to 3.3 mM in tumors of three patients with biopsy-proven isocitrate dehydrogenase (IDH) mutant gliomas. The 2HG was undetectable in an IDH wild-type glioblastoma. For a radiographically suggested glioma, the estimated 2HG of 2.3 ± 0.2 mM (Cramer-Rao lower bound < 10%) indicated that the lesion may be an IDH mutant glioma. The data indicated that the dual-readout alternated gradients EPSI can provide reliable high-resolution imaging of 2HG in glioma patients at 7 T in vivo. Magn Reson Med 79:1851-1861, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Silent navigator-triggered silent MRI of the abdomen.

To develop and demonstrate the feasibility of a silent respiratory navigator technique for prospective triggering, which was incorporated into a three-dimensional radial zero-echo-time sequence for respiratory navigated silent abdominal imaging. A nonselective hard excitation radiofrequency pulse was used for the navigator sequence with a derated readout gradient, to avoid generation of high levels of acoustic noise. The acquired navigator signals were processed in real time and used for prospective triggering of the zero-echo-time sequence. Ten healthy volunteers were scanned using the proposed and conventional techniques at 1.5 T. An acoustic noise measurement with A-weighted continuous equivalent sound pressure level was also performed. The sound pressure-level values of the background noise, zero-echo-time imaging, conventional, and silent navigators were 68.3, 68.4, 102.5, and 69.4 dBA, respectively. Excellent correlation with correlation coefficients greater than 0.9 was observed between the bellows signals and displacement values calculated from the navigators. Sharpness of the portal vein of both conventional and silent navigator-triggered images was significantly higher than those of nontriggered images. The silent navigator-triggered zero-echo-time technique is feasible and might improve image quality and workflow of abdominal MRI of patients who are prone to acoustic noise. Magn Reson Med 79:2170-2175, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Echo-Planar J-resolved Spectroscopic Imaging using Dual Read-outs: Implementation and Quantitation of Human Brain Metabolites

Attempts have been made to reduce the total scan time in multi-dimensional J-resolved spectroscopic imaging (JRESI) using an echo-planar (EP) readout gradient, but acquisition duration remains a limitation for routine clinical use in the brain. We present here a significant acceleration achieved with a 4D EP-JRESI sequence that collects dual phase encoded lines within a single repetition time (TR) using two bipolar read-out trains. The performance and reliability of this novel 4D sequence, called Multi-Echo based Echo-Planar J-resolved Spectroscopic Imaging (ME-EP-JRESI), was evaluated in 10 healthy controls and a brain phantom using a 3 T MRI/MRS scanner. The prior knowledge fitting (ProFit) algorithm, with a new simulated basis set consisting of macromolecules and lipids apart from metabolites of interest, was used for quantitation. Both phantom and in-vivo data demonstrated that localization and spatial/spectral profiles of metabolites from the ME-EP-JRESI sequence were in good agreement with that of the EP-JRESI sequence. Both in the occipital and temporal lobe, metabolites with higher physiological concentrations including Glx (Glu+Gln), tNAA (NAA+NAAG), mI all had coefficient of variations between 9–25%. In summary, we have implemented, validated and tested the ME-EP-JRESI sequence, demonstrating that multi-echo acquisition can successfully reduce the total scan duration for EP-JRESI sequences.

GABA editing with macromolecule suppression using an improved MEGA-SPECIAL sequence.

The most common γ-aminobutyric-acid (GABA) editing approach, MEGA-PRESS, uses J-editing to measure GABA distinct from larger overlapping metabolites, but suffers contamination from coedited macromolecules (MMs) comprising 40 to 60% of the observed signal. MEGA-SPECIAL is an alternative method with better MM suppression, but is not widely used primarily because of its relatively poor spatial localization. Our goal was to develop an improved MM-suppressed GABA editing sequence at 3 Tesla. We modified a single-voxel MEGA-SPECIAL sequence with an oscillating readout gradient for improved spatial localization, and used very selective 30-ms editing pulses for improved suppression of coedited MMs. Simulation and in vivo experiments confirmed excellent MM suppression, insensitive to the range of B0 frequency drifts typically encountered in vivo. Both intersubject and intrasubject studies showed that MMs, when suppressed by the improved MEGA-SPECIAL method, contributed approximately 40% to the corresponding MEGA-PRESS measurements. From the intersubject study, the coefficient of variation for GABA+/Cre (MEGA-PRESS) was 11.2% versus 7% for GABA/Cre (improved MEGA-SPECIAL), demonstrating significantly reduced variance (P = 0.005), likely coming from coedited MMs. This improved MEGA-SPECIAL sequence provides unbiased GABA measurements with reduced variance as compared with conventional MEGA-PRESS. This approach is also relatively insensitive to the range of B0 drifts typically observed in in vivo human studies. Magn Reson Med 79:41-47, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Accelerated susceptibility-based positive contrast imaging of MR compatible metallic devices based on modified fast spin echo sequences

This study aims to develop an accelerated susceptibility-based positive contrast MR imaging method for visualizing MR compatible metallic devices. A modified fast spin echo sequence is used to accelerate data acquisition. Each readout gradient in the modified fast spin echo is slightly shifted by a short distance Tshift. Phase changes accumulated within Tshift are then used to calculate the susceptibility map by using a kernel deconvolution algorithm with a regularized ℓ1 minimization. To evaluate the proposed fast spin echo method, three phantom experiments were conducted and compared to a spin echo based technique and the gold standard CT for visualizing biopsy needles and brachytherapy seeds. Compared to the spin echo based technique, the data sampling speed of the proposed method was faster by 2–4 times while still being able to accurately visualize and identify the location of the biopsy needle and brachytherapy seeds. These results were confirmed by CT images of the same devices. Results also demonstrated that the proposed fast spin echo method can achieve good visualization of the brachytherapy seeds in positive contrast and in different orientations. It is also capable of correctly differentiating brachytherapy seeds from other similar structures on conventional magnitude images.

Rapid chemical shift encoding with single-acquisition single-slab 3D GRASE.

To investigate the feasibility of chemical shift encoded, single-slab 3D GRASE for rapid fat-water separation within a single acquisition. The proposed method incorporates signal-to-noise-ratio-optimal chemical shift encoding into single-slab 3D GRASE with variable flip angles. Chemical shift induced phase information was encoded in succession to different positions in k-space by inserting phase encoding blips between adjacent lobes of the oscillating readout gradients. To enhance imaging efficiency, signal prescription-based variable flip angles were used in the long refocusing pulse train. After echo-independent phase correction, missing signals in k-echo space were interpolated using convolution kernels that span over all echoes. Fat-water separation in a single acquisition was performed using both multi-echo fast spin echo and GRASE as compared to conventional multiacquisition fast spin echo with echo shifts. The proposed single-slab 3D GRASE shows superior performance in accurately delineating cartilage structures compared to its counterpart, multi-echo 3D fast spin echo. Compared with multiacquisition fast spin echo with three echo shifts (63 min), the proposed method substantially speeds up imaging time (7 min), and achieves 0.6 mm isotropic resolution in knee imaging with reduced artifacts and noise. We successfully demonstrated the feasibility of rapid chemical shift encoding and separation using the proposed, single-acquisition single-slab 3D GRASE for high resolution isotropic imaging within clinically acceptable time. Magn Reson Med 78:1852-1861, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Image-based gradient non-linearity characterization to determine higher-order spherical harmonic coefficients for improved spatial position accuracy in magnetic resonance imaging

PurposeSpatial position accuracy in magnetic resonance imaging (MRI) is an important concern for a variety of applications, including radiation therapy planning, surgical planning, and longitudinal studies of morphologic changes to study neurodegenerative diseases. Spatial accuracy is strongly influenced by gradient linearity. This work presents a method for characterizing the gradient non-linearity fields on a per-system basis, and using this information to provide improved and higher-order (9th vs. 5th) spherical harmonic coefficients for better spatial accuracy in MRI. MethodsA large fiducial phantom containing 5229 water-filled spheres in a grid pattern is scanned with the MR system, and the positions all the fiducials are measured and compared to the corresponding ground truth fiducial positions as reported from a computed tomography (CT) scan of the object. Systematic errors from off-resonance (i.e., B0) effects are minimized with the use of increased receiver bandwidth (±125kHz) and two acquisitions with reversed readout gradient polarity. The spherical harmonic coefficients are estimated using an iterative process, and can be subsequently used to correct for gradient non-linearity. Test-retest stability was assessed with five repeated measurements on a single scanner, and cross-scanner variation on four different, identically-configured 3T wide-bore systems. ResultsA decrease in the root-mean-square error (RMSE) over a 50cm diameter spherical volume from 1.80mm to 0.77mm is reported here in the case of replacing the vendor's standard 5th order spherical harmonic coefficients with custom fitted 9th order coefficients, and from 1.5mm to 1mm by extending custom fitted 5th order correction to the 9th order. Minimum RMSE varied between scanners, but was stable with repeated measurements in the same scanner. ConclusionsThe results suggest that the proposed methods may be used on a per-system basis to more accurately calibrate MR gradient non-linearity coefficients when compared to vendor standard corrections.

Current-induced alternating reversed dual-echo-steady-state for joint estimation of tissue relaxation and electrical properties.

To develop a current-induced, alternating reversed dual-echo-steady-state-based magnetic resonance electrical impedance tomography for joint estimation of tissue relaxation and electrical properties. The proposed method reverses the readout gradient configuration of conventional, in which steady-state-free-precession (SSFP)-ECHO is produced earlier than SSFP-free-induction-decay (FID) while alternating current pulses are applied in between the two SSFPs to secure high sensitivity of SSFP-FID to injection current. Additionally, alternating reversed dual-echo-steady-state signals are modulated by employing variable flip angles over two orthogonal injections of current pulses. Ratiometric signal models are analytically constructed, from which T1 , T2 , and current-induced Bz are jointly estimated by solving a nonlinear inverse problem for conductivity reconstruction. Numerical simulations and experimental studies are performed to investigate the feasibility of the proposed method in estimating relaxation parameters and conductivity. The proposed method, if compared with conventional magnetic resonance electrical impedance tomography, enables rapid data acquisition and simultaneous estimation of T1 , T2 , and current-induced Bz , yielding a comparable level of signal-to-noise ratio in the parameter estimates while retaining a relative conductivity contrast. We successfully demonstrated the feasibility of the proposed method in jointly estimating tissue relaxation parameters as well as conductivity distributions. It can be a promising, rapid imaging strategy for quantitative conductivity estimation. Magn Reson Med 78:107-120, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Thrombus-mimicking artifacts in two-point Dixon MRI: Prevalence, appearance, and severity.

To evaluate the incidence and severity of potentially thrombus mimicking, flow-induced misallocation artifacts in a clinical setting. Two-point "Dixon" fat-water separation methods, with bipolar readout gradients, may suffer from flow-induced fat-water misallocation artifacts. If these artifacts occur within blood vessels, they may mimic thrombus. Two-point Dixon coronal and axial images acquired in 102 consecutive patients were retrospectively evaluated for the presence of flow-induced artifacts in arteries and veins. Artifacts were graded on a 3-point scale (none, mild, severe) by two independent readers. Interreader agreement was evaluated with kappa statistics. Reader 1 reported 63 artifacts in 46 (45%) of the cases (severe in 19 cases, 18.6%). Reader 2 reported 51 artifacts in 43 (42.2%) of the cases (severe in 18 cases, 17.6%). Misallocation of fat and water was apparent in all datasets with severe artifacts, whereas variable signal intensity changes in water and fat images were observed in mild artifacts. Interreader agreement was good for artifacts appearing in coronal images (κ = 0.7) and fair for artifact appearance in axial images (κ = 0.24). Our study shows a high incidence of flow-induced mild and severe artifacts in a two-point Dixon method with bipolar readout gradients. This artifact should not be misinterpreted as intravascular thrombus. 3 J. Magn. Reson. Imaging 2017;45:229-236.

A new technique for MR elastography of the supraspinatus muscle: A gradient-echo type multi-echo sequence

Magnetic resonance elastography (MRE) can measure tissue stiffness quantitatively and noninvasively. Supraspinatus muscle injury is a significant problem among throwing athletes. The purpose of this study was to develop an MRE technique for application to the supraspinatus muscle by using a conventional magnetic resonance imaging (MRI). MRE acquisitions were performed with a gradient-echo type multi-echo MR sequence at 100Hz pneumatic vibration. A custom-designed vibration pad was used as a pneumatic transducer in order to adapt to individual shoulder shapes. In a gradient-echo type multi-echo MR sequence, without motion encoding gradient (MEG) that synchronizes with vibrations, bipolar readout gradient lobes achieved a similar function to MEG (MEG-like effect). In other words, a dedicated MRE sequence (built-in MEG) is not always necessary for MRE. In this study, 7 healthy volunteers underwent MRE. We investigated the effects of direction of the MEG-like effect and selected imaging planes on the patterns of wave propagation (wave image). The results indicated that wave images showed clear wave propagation on a condition that the direction of the MEG-like effect was nearly perpendicular to the long axis of the supraspinatus muscle, and that the imaging plane was superior to the proximal supraspinatus muscle. This limited condition might be ascribed to specific features of fibers in the supraspinatus muscle and wave reflection from the boundaries of the supraspinous fossa. The mean stiffness of the supraspinatus muscle was 10.6±3.17kPa. Our results demonstrated that using MRE, our method can be applied to the supraspinatus muscle by using conventional MRI.

Probing lactate secretion in tumours with hyperpolarised NMR.

Most tumours exhibit a high rate of glycolysis and predominantly produce energy by lactic acid fermentation. To maintain energy production and prevent toxicity, the lactate generated needs to be rapidly transported out of the cell. This is achieved by monocarboxylate transporters (MCTs), which therefore play an essential role in cancer metabolism and development. In vivo experiments were performed on eight male Fisher F344 rats bearing a subcutaneous mammary carcinoma after injection of hyperpolarised [1-(13) C]pyruvate. A Gd(III)DO3A complex that binds to pyruvate and its metabolites was used to efficiently destroy the extracellular magnetisation after hyperpolarised lactate had been formed. Moreover, a pulse sequence including a frequency-selective saturation pulse was designed so that the pyruvate magnetisation could be destroyed to exclude effects arising from further conversion. Given this preparation, metabolite transport out of the cell manifested as additional decay and apparent cell membrane transporter rates could thus be obtained using a reference measurement without a relaxation agent. In addition to slice-selective spectra, spatially resolved maps of apparent membrane transporter activity were acquired using a single-shot spiral gradient readout. A considerable increase in decay rate was detected for lactate, indicating rapid transport out of the cell. The alanine signal was unaltered, which corresponds to a slower efflux rate. This technique could allow for better understanding of tumour metabolism and progression, and enable treatment response measurements for MCT-targeted cancer therapies. Moreover, it provides vital insights into the signal kinetics of hyperpolarised [1-(13) C]pyruvate examinations. Copyright © 2016 John Wiley & Sons, Ltd.

Joint correction of Nyquist artifact and minuscule motion-induced aliasing artifact in interleaved diffusion weighted EPI data using a composite two-dimensional phase correction procedure

Diffusion-weighted imaging (DWI) obtained with interleaved echo-planar imaging (EPI) pulse sequence has great potential of characterizing brain tissue properties at high spatial-resolution. However, interleaved EPI based DWI data may be corrupted by various types of aliasing artifacts. First, inconsistencies in k-space data obtained with opposite readout gradient polarities result in Nyquist artifact, which is usually reduced with 1D phase correction in post-processing. When there exist eddy current cross terms (e.g., in oblique-plane EPI), 2D phase correction is needed to effectively reduce Nyquist artifact. Second, minuscule motion induced phase inconsistencies in interleaved DWI scans result in image-domain aliasing artifact, which can be removed with reconstruction procedures that take shot-to-shot phase variations into consideration. In existing interleaved DWI reconstruction procedures, Nyquist artifact and minuscule motion-induced aliasing artifact are typically removed subsequently in two stages. Although the two-stage phase correction generally performs well for non-oblique plane EPI data obtained from well-calibrated system, the residual artifacts may still be pronounced in oblique-plane EPI data or when there exist eddy current cross terms. To address this challenge, here we report a new composite 2D phase correction procedure, which effective removes Nyquist artifact and minuscule motion induced aliasing artifact jointly in a single step. Our experimental results demonstrate that the new 2D phase correction method can much more effectively reduce artifacts in interleaved EPI based DWI data as compared with the existing two-stage artifact correction procedures. The new method robustly enables high-resolution DWI, and should prove highly valuable for clinical uses and research studies of DWI.

Investigation of the 4D composite MR image distortion field associated with tumor motion for MR-guided radiotherapy.

Magnetic resonance (MR) images are affected by geometric distortions due to the specifics of the MR scanner and patient anatomy. Quantifying the distortions associated with mobile tumors is particularly challenging due to real anatomical changes in the tumor's volume, shape, and relative location within the MR imaging volume. In this study, the authors investigate the 4D composite distortion field, which combines the effects of the susceptibility-induced and system-related distortion fields, experienced by mobile lung tumors. The susceptibility (χ) effects were numerically simulated for two specific scenarios: (a) a full motion cycle of a lung tumor due to breathing as depicted on ten phases of a 4D CBCT data set and (b) varying the tumor size and location in lung tissue via a synthetically generated sphere with variable diameter (4-80 mm). The χ simulation procedure relied on the segmentation and generation of 3D susceptibility (χ) masks and computation of the magnetic field by means of finite difference methods. A system-related distortion field, determined with a phantom and image processing algorithm, was used as a reference. The 4D composite distortion field was generated as the vector summation of the χ-induced and system-related fields. The analysis was performed for two orientations of the main magnetic field (B0), which correspond to several MRIgRT system configurations. Specifically, B0 was set along the z-axis as in the case of a cylindrical-bore scanner and in the (x,y)-plane as for a biplanar MR. Computations were also performed for a full revolution at 15° increments in the case of a rotating biplanar magnet. Histograms and metrics such as maximum, mean, and range were used to evaluate the characteristics of the 4D distortion field. The χ-induced field depends on the change in volume and shape of the moving tumor as well as the local surrounding anatomy. In the case of system-related distortions, the tumor experiences increased field perturbations as it moves further away from the MR isocenter. For a mobile lung tumor, the 4D composite field, corresponding to a 1.5 T field and a readout gradient of 5 mT/m, amounts to 3.0 and 2.8 mm for the MRIgRT system designs featuring B0 oriented along the z-axis (cylindrical-bore scanner) and in the (x,y)-plane (biplanar scanner), respectively. For a rotating biplanar scanner, the composite distortion field varied nonlinearly with the rotation angle. Overall, the dominant contribution to the composite field was from the system-related distortion field. The tumor centroid experienced a systematic shift of 2 mm and showed a negligible perturbation for different B0 values. The dependency on the tumor size was also investigated, namely the max values varied from 1.2 to 2.5 mm for spherical volumes with a diameter between 4 and 80 mm. The composite distortion field requires adequate quantification for lung radiation therapy applications such as treatment planning, pretreatment patient setup verification, and real-time treatment delivery guidance. For certain scenarios such as small tumor volumes, the spatial distortions may be corrected by applying systematic shifts derived from a single tumor motion phase. In the case of high readout gradients common to fast imaging applications, the χ distortions were found to be less than 1 mm irrespective of scanner configuration.

Enhanced magnetic flux density mapping using coherent steady state equilibrium signal in MREIT

Measuring the z-component of magnetic flux density B = (Bx, By, Bz) induced by transversally injected current, magnetic resonance electrical impedance tomography (MREIT) aims to visualize electrical property (current density and/or conductivity distribution) in a three-dimensional imaging object. For practical implementations of MREIT technique, it is critical to reduce injection of current pulse within safety requirements. With the goal of minimizing the noise level in measured Bz data, we propose a new method to enhance the measure Bz data using steady-state coherent gradient multi-echo (SSC-GME) MR pulse sequence combining with injection current nonlinear encoding (ICNE) method in MREIT, where the ICNE technique injects current during a readout gradient to maximize the signal intensity of phase signal including Bz. The total phase offset in SSC-GME includes additional magnetic flux density due to the injected current, which is different from the phase signal for the conventional spoiled MR pulse sequence. We decompose the magnetization precession phase from the total phase offset including Bz and optimize Bz data using the steady-state equilibrium signal. Results from a real phantom experiment including different kinds of anomalies demonstrated that the proposed method enhanced Bz comparing to a conventional spoiled pulse sequence.