Readout Direction Research Articles

Artifact characterization of Nitinol needles in magnetic resonance imaging-guided musculoskeletal interventions at 3.0 tesla: a phantom study.

To characterize the artifacts of an 18-gauge coaxial nickel-titanium needle using a balanced steady-state free precession sequence in magnetic resonance imaging-guided interventions at 3.0 tesla. The influence of flip angle (FA), bandwidth, matrix, slice thickness (ST), and read-out direction on needle artifact behavior was investigated for different intervention angles (IA). Artifact diameters were rated at predefined positions. Subgroup differences were assessed using Bonferroni-corrected non-parametric tests and correlations between continuous variables were expressed using the Bravais-Pearson coefficient. Interrater reliability was quantified using intraclass correlation coefficients (ICCs), and a contrast-enhanced target lesion to non-enhanced muscle tissue contrast ratio was quantified. The artifact diameters decreased with an increase in FA for all IAs (P < 0.001) and with an increase in ST for IAs of 45°-90° (all P < 0.05). Tip artifacts occurred at low IAs (0°-45°) and gradually increased in size with a decrease in IA (P = 0.022). The interrater reliability was high (ICC: 0.994-0.999). The contrast-enhanced target lesion to non-enhanced muscle tissue contrast ratio presented positive correlations with increasing FAs and matrices (P < 0.001; P = 0.003) and a negative correlation with increasing STs (P = 0.007). To minimize needle artifacts, it is recommended to use FAs of 40°-60°, a ST of >7 mm, and, if possible, an IA of 45°-60°. The visibility of the target lesion and the needle's artifact behavior must be weighed up against each other when choosing the ST, while higher FAs (40°-60°) and matrices (224 × 224/256 × 256) are associated with low artifacts and sufficient lesion visibility.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL).

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

Improving spreading projection algorithm for rapid k-space sampling trajectories through minimized off-resonance effects and gridding of low frequencies.

Non-Cartesian MRI with long arbitrary readout directions are susceptible to off-resonance artifacts due to patient induced inhomogeneities. This results in degraded image quality with strong signal losses and blurring. Current solutions to address this issue involve correcting the off-resonance artifacts during image reconstruction or reducing inhomogeneities through improved shimming. The recently developed SPARKLING algorithm is extended to drastically diminish off-resonance artifacts by generating temporally smooth k-space sampling patterns. For doing so, the cost function which is optimized in SPARKLING is modified using a temporal weighting factor. Additionally, oversampling of the center of k-space beyond the Nyquist criteria is prevented through the use of gridded sampling in the region, enforced with affine constraints. Prospective k-space data was acquired at 3 T on new trajectories, and we show robustness to inhomogeneities through in silico experiments by adding through artificial degradation of system shimming. Later on, in vivo experiments were carried out to optimize parameters of the new improvements and benchmark the gain in performance. The improved trajectories allowed for the recovery of signal dropouts observed on original SPARKLING acquisitions at larger field inhomogeneities. Furthermore, imposing gridded sampling at the center of k-space provided improved reconstructed image quality with limited artifacts. These advancements allowed us for nearly shorter scan time compared to GRAPPA-p4x1, allowing us to reach 600 µm isotropic resolution in 3D -w imaging in just 3.3 min at 3 T with negligible degradation in image quality.

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PurposeTo develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution. MethodsWe first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI. ResultsBoth simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches. ConclusionsThe proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

Evaluation of simultaneous multi-slice readout-segmented diffusion-weighted MRI acquisition in human optic nerve measurements

Diffusion-weighted magnetic resonance imaging (dMRI) is the only available method to measure the tissue properties of white matter tracts in living human brains and has opened avenues for neuroscientific and clinical studies on human white matter. However, dMRI using conventional simultaneous multi-slice (SMS) single-shot echo planar imaging (ssEPI) still presents challenges in the analyses of some specific white matter tracts, such as the optic nerve, which are heavily affected by susceptibility-induced artifacts. In this study, we evaluated dMRI data acquired by using SMS readout-segmented EPI (rsEPI), which aims to reduce susceptibility-induced artifacts by dividing the acquisition space into multiple segments along the readout direction to reduce echo spacing. To this end, we acquired dMRI data from 11 healthy volunteers by using SMS ssEPI and SMS rsEPI, and then compared the dMRI data of the human optic nerve between the SMS ssEPI and SMS rsEPI datasets by visual inspection of the datasets and statistical comparisons of fractional anisotropy (FA) values. In comparison with the SMS ssEPI data, the SMS rsEPI data showed smaller susceptibility-induced distortion and exhibited a significantly higher FA along the optic nerve. In summary, this study demonstrates that despite its prolonged acquisition time, SMS rsEPI is a promising approach for measuring the tissue properties of the optic nerve in living humans and will be useful for future neuroscientific and clinical investigations of this pathway.

Development of an ultrafast brain MR neuronavigation protocol for ventricular shunt placement.

Advancements in MRI technology have provided improved ways to acquire imaging data and to more seamlessly incorporate MRI into modern pediatric surgical practice. One such situation is image-guided navigation for pediatric neurosurgical procedures, including intracranial catheter placement. Image-guided surgery (IGS) requires acquisition of CT or MR images, but the former carries the risk of ionizing radiation and the latter is associated with long scan times and often requires pediatric patients to be sedated. The objective of this project was to circumvent the use of CT and standard-sequence MRI in ventricular neuronavigation by investigating the use of fast MR sequences on the basis of 3 criteria: scan duration comparable to that of CT acquisition, visualization of ventricular morphology, and image registration with surface renderings comparable to standard of care. The aim of this work was to report image development, implementation, and results of registration accuracy testing in healthy subjects. The authors formulated 11 candidate MR sequences on the basis of the standard IGS protocol, and various scan parameters were modified, such as k-space readout direction, partial k-space acquisition, sparse sampling of k-space (i.e., compressed sensing), in-plane spatial resolution, and slice thickness. To evaluate registration accuracy, the authors calculated target registration error (TRE). A candidate sequence was selected for further evaluation in 10 healthy subjects. The authors identified a candidate imaging protocol, termed presurgical imaging with compressed sensing for time optimization (PICO). Acquisition of the PICO protocol takes 25 seconds. The authors demonstrated noninferior TRE for PICO (3.00 ± 0.19 mm) in comparison with the default MRI neuronavigation protocol (3.35 ± 0.20 mm, p = 0.20). The developed and tested sequence of this work allowed accurate intraoperative image registration and provided sufficient parenchymal contrast for visualization of ventricular anatomy. Further investigations will evaluate use of the PICO protocol as a substitute for CT and conventional MRI protocols in ventricular neuronavigation.

Artifact reduction of coaxial needles in magnetic resonance imaging-guided abdominal interventions at 1.5 T: a phantom study

Needle artifacts pose a major limitation for MRI-guided interventions, as they impact the visually perceived needle size and needle-to-target-distance. The objective of this agar liver phantom study was to establish an experimental basis to understand and reduce needle artifact formation during MRI-guided abdominal interventions. Using a vendor-specific prototype fluoroscopic T1-weighted gradient echo sequence with real-time multiplanar acquisition at 1.5 T, the influence of 6 parameters (flip angle, bandwidth, matrix, slice thickness, read-out direction, intervention angle relative to B0) on artifact formation of 4 different coaxial MR-compatible coaxial needles (Nitinol, 16G–22G) was investigated. As one parameter was modified, the others remained constant. For each individual parameter variation, 2 independent and blinded readers rated artifact diameters at 2 predefined positions (15 mm distance from the perceived needle tip and at 50% of the needle length). Differences between the experimental subgroups were assessed by Bonferroni-corrected non-parametric tests. Correlations between continuous variables were expressed by the Bravais–Pearson coefficient and interrater reliability was quantified using the intraclass classification coefficient. Needle artifact size increased gradually with increasing flip angles (p = 0.002) as well as increasing intervention angles (p < 0.001). Artifact diameters differed significantly between the chosen matrix sizes (p = 0.002) while modifying bandwidth, readout direction, and slice thickness showed no significant differences. Interrater reliability was high (intraclass correlation coefficient 0.776–0.910). To minimize needle artifacts in MRI-guided abdominal interventions while maintaining optimal visibility of the coaxial needle, we suggest medium-range flip angles and low intervention angles relative to B0.

Quantitative effects of off-resonance related distortion on brain mechanical property estimation with magnetic resonance elastography.

Off-resonance related geometric distortion can impact quantitative MRI techniques, such as magnetic resonance elastography (MRE), and result in errors to these otherwise sensitive metrics of brain health. MRE is a phase contrast technique to determine the mechanical properties of tissue by imaging shear wave displacements and estimating tissue stiffness through inverse solution of Navier's equation. In this study, we systematically examined the quantitative effects of distortion and corresponding correction approaches on MRE measurements through a series of simulations, phantom models, and in vivo brain experiments. We studied two different k-space trajectories, echo-planar imaging and spiral, and we determined that readout time, off-resonance gradient strength, and the combination of readout direction and off-resonance gradient direction, impact the estimated mechanical properties. Images were also processed through traditional distortion correction pipelines, and we found that each of the correction mechanisms works well for reducing stiffness errors, but are limited in cases of very large distortion. The ability of MRE to detect subtle changes to neural tissue health relies on accurate, artifact-free imaging, and thus off-resonance related geometric distortion must be considered when designing sequences and protocols by limiting readout time and applying correction where appropriate.

Effect of beam quality and readout direction in the edge profile on the modulation transfer function of photostimulable phosphor systems via the edge method.

Purpose: The objective of this study was to investigate the effects of (i) the difference in the beam qualities on the presampled modulation transfer function (MTF) using the edge method and (ii) the readout direction of the edge profile for the photostimulable phosphor (PSP) system. Approach: The International Electrotechnical Commission (IEC) defined a technique using the "radiation qualities based on a phantom made up of an aluminum added filter" (RQA). A general radiographic system with a tube voltage of 50kV and 9.7mm of additional aluminum filtration was used to conform the x-ray to the IEC-specified beam quality definition RQA3. Additionally, we employed two different beams with tube voltages of 60 and 70kV using a dental x-ray unit. The MTF was measured in the readout direction from low-to-high exposure regions and vice versa with respect to the scanning and subscanning directions of the PSP system. Results: The difference in the 50%MTF value for all directions between 60 and 70kV averaged less than . The 50%MTF of RQA3 was on average 0.2 lower than the value for 60kV for each direction. For all beam qualities, no difference was observed between the MTFs measured in the readout direction from the low-to-high exposure regions and vice versa. Conclusions: The MTFs, measured using the dental x-ray unit, were unaffected by the tube voltage, and they were slightly higher than those measured using the RQA3. Furthermore, the MTF was unaffected by the differences in the readout directions of the edge profile.

Effects of B0 eddy currents on imaging isocenter shifts in 0.35‐T MRI‐guided radiotherapy (MR‐IGRT) system

PurposeThe purpose of this study was to measure gantry angle‐related eddy currents in a 0.35‐T MRI‐Linac and determine if B0 (zeroth order) eddy currents are the primary cause of gantry angle‐dependent imaging isocenter shifts vs other potential causes like B0 inhomogeneities and gradient (first order) eddy currents. For conventional Cartesian acquisitions, B0 eddy currents can cause imaging isocenter shifts along both phase encode and readout directions. Gradient eddy currents can cause spatial distortion along both the phase encode and readout directions. Center frequency offsets can cause imaging isocenter shifts along the readout direction that vary with readout gradient polarity.MethodsMRI‐related eddy currents and imaging isocenter shifts were measured on a 0.35‐T MRI‐Linac at gantry angles from 0° to 330° in increments of 30°. All measurements were made after gradient shimming and center frequency tuning at each planned gantry angle. Eddy current and field homogeneity measurements were conducted using a 24‐cm diameter spherical phantom. Gradient and B0 eddy currents were calculated from the free induction decays (FIDs) resulting from selective excitation of slices located ±5 cm from isocenter. B0 eddy currents were also calculated from FIDs acquired with nonselective excitation and compared with B0 eddy current values derived using selective excitation. B0 inhomogeneities and center frequency offsets were measured by acquiring FIDs with nonselective excitation. Imaging isocenter shifts were measured using a 33x33x10.5 cm3 uniformity linearity (grid) phantom and a 3D true fast imaging with steady‐state precession (TrueFISP) sequence used in MRI‐guided radiation therapy. Eddy currents were compared to vendor specifications and correlated with the imaging isocenter shifts. Measurements were conducted before and after the MRI‐Linac’s waveguide was replaced with an updated design to reduce eddy currents.ResultsB0 eddy currents were highly correlated (r = 0.986, P &lt;&lt; 0.001) for measurements made with vs without selective excitation. Transverse (X and Y) axis B0 eddy currents before and after the waveguide upgrade were out of specification (specification: ≤0.1 μT m/mT for delays &lt; 10 ms) for most of the measured gantry angles. Gradient eddy currents before and after the upgrade were within specifications for the measured gantry angles (≤0.1% for delays &lt; 10 ms). B0 eddy currents and imaging isocenter shifts were highly correlated (r = 0.965, P &lt;&lt; 0.001). After the Linac waveguide upgrade, root mean square (RMS) peak B0 and gradient eddy currents dropped 45% and 11%, respectively, for delays &lt;10 ms, while imaging isocenter shifts dropped 53%. Isocenter shifts were observed in both phase encode and readout directions. Center frequency offsets were &lt;26 Hz while B0 inhomogeneities were &lt;33 Hz full width at half maximum (FWHM).ConclusionsImaging isocenter shifts measured in a 0.35‐T MRI‐Linac were highly correlated with B0 eddy currents. The eddy currents and imaging isocenter shifts decreased after the MRI‐Linac’s waveguide was replaced.

A novel quality assurance procedure for trajectory log validation using phantom-less real-time latency corrected EPID images.

The use of trajectory log files for routine patient quality assurance is gaining acceptance. Such use requires the validation of the trajectory log itself. However, the accurate localization of a multileaf collimator (MLC) leaf while it is in motion remains a challenging task. We propose an efficient phantom‐less technique using the EPID to verify the dynamic MLC positions with high accuracy. Measurements were made on four Varian TrueBeams equipped with M120 MLCs. Two machines were equipped with the S1000 EPID; two were equipped with the S1200 EPID. All EPIDs were geometrically corrected prior to measurements. Dosimetry mode EPID measurements were captured by a frame grabber card directly linked to the linac. All leaf position measurements were corrected both temporally and geometrically. The readout latency of each panel, as a function of pixel row, was determined using a 40 × 1.0 cm2 sliding window (SW) field moving at 2.5 cm/s orthogonal to the row readout direction. The latency of each panel type was determined by averaging the results of two panels of the same type. Geometric correction was achieved by computing leaf positions with respect to the projected isocenter position as a function of gantry angle. This was determined by averaging the central axis position of fields at two collimator positions of 90° and 270°. The radiological to physical leaf end position was determined by comparison of the measured gap with that determined using a feeler gauge. The radiological to physical leaf position difference was found to be 0.1 mm. With geometric and latency correction, the proposed method was found to be improve the ability to detect dynamic MLC positions from 1.0 to 0.2 mm for all leaves. Latency and panel residual geometric error correction improve EPID‐based MLC position measurement. These improvements provide for the first time a trajectory log QA procedure.

Minimizing echo and repetition times in magnetic resonance imaging using a double half-echo k-space acquisition and low-rank reconstruction.

Sampling k-space asymmetrically (ie, partial Fourier sampling) in the readout direction is a common way to reduce the echo time (TE) during magnetic resonance image acquisitions. This technique requires overlap around the center of k-space to provide a calibration region for reconstruction, which limits the minimum fractional echo to ~60% before artifacts are observed. The present study describes a method for reconstructing images from exact half echoes using two separate acquisitions with reversed readout polarity, effectively providing a full line of k-space without additional data around central k-space. This approach can benefit sequences or applications that prioritize short TE, short inter-echo spacing or short repetition time. An example of the latter is demonstrated to reduce banding artifacts in balanced steady-state free precession.

Measurements of venous oxygen saturation in the superior sagittal sinus using conventional 3D multiple gradient-echo MRI: Effects of flow velocity and acceleration.

This work investigates the effects of flow acceleration in the superior sagittal sinus on slice-dependent variations in venous oxygen saturation (SvO2 ) estimations using susceptibility-based MR oximetry. Three-dimensional multiple gradient-echo images, with first-order flow compensation along the anterior-posterior readout direction for the first echo, were acquired twice from 15 healthy volunteers. For all slices, phases within the superior sagittal sinus were fitted using linear regression across four TEs to obtain the Pearson's correlation coefficients (PCCs), the largest of which corresponded to minimum acceleration influence. SvO2 derived from odd echoes on this slice was used to assess interscan difference, and compared with the central 15th slice for slice-dependent difference, both using Bland-Altman analysis. Within-scan interslice SvO2 consistency was examined versus PCC. Multislice-averaged SvO2 values were then computed from slices with PCCs above a certain threshold. Slice-dependent difference in SvO2 varied from -16.2% to 21.5% at two SDs, in agreement with a recent report, and about twice larger than interscan differences for the automatically selected slice (-7.5% to 10.3%) and for the central 15th slice (-8.0% to 8.8%). For slices with PCCs higher than -0.98, interslice SvO2 deviations were all found to be less than 5.0%. Multislice-averaged SvO2 with PCCs higher than -0.98 further reduced interscan difference to -4.7% to 8.2%. Slice-dependent variations in SvO2 may partly be explained by the effects of flow acceleration. Our method may enable conventional 3D multiple gradient echo to be used for SvO2 estimations in the presence of pulsatile flow.

Applicator visualization using ultrashort echo time MRI for high-dose-rate endorectal brachytherapy

PurposeThe individual channels in an endorectal applicator for high-dose-rate endorectal brachytherapy are not visible on standard MRI sequences. The aim of this study was to test whether an ultrashort echo time (UTE) MRI sequence could be used to visualize the individual channels to enable MR-only treatment planning for rectal cancer. Methods and MaterialsWe used a radial three-dimensional (3D) UTE pulse sequence and acquired images of phantoms and two patients with rectal cancer. We rigidly registered a UTE image and CT scan of an applicator phantom, based on the outline of the applicator. One observer compared channel positions on the UTE image and CT scan in five slices spaced 25 mm apart. To quantify geometric distortions, we scanned a commercial 3D geometric quality assurance phantom and calculated the difference between detected marker positions on the UTE image and corresponding marker positions on two 3D T1-weighted images with opposing readout directions. ResultsOn the UTE images, there is sufficient contrast to discern the individual channels. The difference in channel positions on the UTE image compared with the CT was on average −0.1 ± 0.1 mm (left-right) and 0.1 ± 0.3 mm (anteroposterior). After rigid registration to the 3D T1-weighted sequences, the residual 95th percentile of the geometric distortion inside a 550-mm-diameter sphere was 1.0 mm (left-right), 0.9 mm (anteroposterior), and 0.9 mm (craniocaudal). ConclusionsWith a UTE sequence, the endorectal applicator and individual channels can be adequately visualized in both phantom and patients. The geometrical fidelity is within an acceptable range.

Appearance Learning for Image-Based Motion Estimation in Tomography.

In tomographic imaging, anatomical structures are reconstructed by applying a pseudo-inverse forward model to acquired signals. Geometric information within this process is usually depending on the system setting only, i.e., the scanner position or readout direction. Patient motion therefore corrupts the geometry alignment in the reconstruction process resulting in motion artifacts. We propose an appearance learning approach recognizing the structures of rigid motion independently from the scanned object. To this end, we train a siamese triplet network to predict the reprojection error (RPE) for the complete acquisition as well as an approximate distribution of the RPE along the single views from the reconstructed volume in a multi-task learning approach. The RPE measures the motion-induced geometric deviations independent of the object based on virtual marker positions, which are available during training. We train our network using 27 patients and deploy a 21-4-2 split for training, validation and testing. In average, we achieve a residual mean RPE of 0.013mm with an inter-patient standard deviation of 0.022mm. This is twice the accuracy compared to previously published results. In a motion estimation benchmark the proposed approach achieves superior results in comparison with two state-of-the-art measures in nine out of twelve experiments. The clinical applicability of the proposed method is demonstrated on a motion-affected clinical dataset.

A locally segmented reconstruction method for parallel imaging.

A locally segmented parallel imaging reconstruction method is proposed that efficiently utilizes sensitivity distribution of multichannel receiver coil. A method of locally segmenting a MR signal is introduced to maximize the differences in sensitivity between receiver channels. A 1D Fourier transformation of the undersampled k-space data is performed along the readout direction, which generates a hybrid 2D space. The hybrid space is partitioned into localized segments along the readout direction. In every localized segment, kernels representing relation between adjacent signals are estimated from autocalibration signals, and data at unsampled points are estimated using the kernels. Then, the images are reconstructed from full k-space data that consists of the sampled data and the estimated data at unsampled points. In a computer simulation and in vivo experiments, the locally segmented reconstruction method produced fewer residual artifacts compared to the conventional parallel imaging reconstruction methods with the same kernel geometry. The performance gain of the proposed method comes from maximizing encoding capability of receiver channels, thus resulting in the accurately estimated kernel weights that reflect the relation between adjacent signals. The proposed spatial segmentation method maximally utilizes differences in the sensitivity of receiver channels to reconstruct images with reduced artifacts.

Application of compressed sensing using chirp encoded 3D GRE and MPRAGE sequences

AbstractAn implementation of Non‐Fourier chirp‐encoding in 3D Gradient Recalled Echo (GRE), susceptibility‐weighted imaging (SWI) and Magnetization Prepared Rapid Gradient Echo (MPRAGE) sequences is presented with compressive sensing reconstruction. 3D GRE and MPRAGE sequences were designed, in which the phase encoding (PE) direction was encoded with spatially selective chirp encoding Radio Frequency (RF) pulses, while the slice and the readout directions were Fourier encoded using gradients. During each excitation along the PE direction, a different spatially‐selective RF excitation pulse was used to encode the PE direction with a complete set of unitary chirp encoding basis. Multichannel compressive sensing reconstruction on the undersampled in vivo data demonstrated that images reconstructed from chirp encoded data were able to preserve the spatial resolution better than the Fourier encoding. The mean Structural Similarity (SSIM) across five subjects at the acceleration factor of 6, for chirp encoded MPRAGE was 0.934 compared to 0.912 for Fourier encoded MPRAGE. The implementation of prospective undersampling demonstrated the feasibility of using chirp encoding in clinical practice for accelerated imaging. The minimum intensity projection of the compressive sensing (CS) reconstructed susceptibility weighted images revealed that chirp encoding is able to delineate small vessels better than the Fourier encoding with the SSIM of 0.960 for chirp encoding compared to the SSIM of 0.949 for the Fourier encoding. Improved performance of chirp encoding for CS reconstruction and SWI, along with the feasibility of implementation makes them a practical candidate for clinical MRI scans.

A new approach to accelerate readout segmented EPI with compressed sensing.

High resolution diffusion-weighted imaging is limited by susceptibility-induced distortions and relaxation-induced blurring. Segmented acquisition techniques can address these limitations at the expense of a prolonged scan time. If segmentation is performed along the readout direction, e.g., in RESOLVE (readout segmentation of long and variable echo-trains), scan time can be reduced by readout (RO) partial Fourier methods, or simultaneous multi-slice (SMS) methods. In this paper, we present a new approach to additionally accelerate the image acquisition called variable segment (VASE) RESOLVE. To avoid discontinuities at the boundaries of the segments, the phase evolution and therefore the effective echo-spacing needs to be adjusted. To achieve this, we use higher undersampling factors in the outer parts of k-space. Simultaneously we increase the width of the outer segments resulting in an increase of the echo-spacing. Because of this variation, we introduce a kind of randomization to the sampling scheme. This enables the use of compressed sensing reconstruction techniques, which results in improved image quality compared to standard parallel imaging methods. The RMS errors for the VASE RESOLVE acquisitions were lower compared to the standard reconstructions. The VASE RESOLVE in vivo images show a higher apparent signal to noise ratio. VASE RESOLVE is a new approach to further decrease the acquisition time of RO segmented acquisitions. Compared to RESOLVE with SMS, VASE RESOLVE additionally reduces the acquisition time by a factor of 2.

Velocity encoding and velocity compensation for multi-spoke RF excitation

PurposeTo investigate velocity encoded and velocity compensated variants of multi-spoke RF pulses that can be used for flip-angle homogenization at ultra-high fields (UHF). Attention is paid to the velocity encoding for each individual spoke pulse and to displacement artifacts that arise in Fourier transform imaging in the presence of flow. Theory and methodsA gradient waveform design for multi-spoke excitation providing an algorithm for minimal TE was proposed that allows two different encodings. Such schemes were compared to an encoding approach that applies an established scheme to multi-spoke excitations. The impact on image quality and quantitative velocity maps was evaluated in phantoms using single- and two-spoke excitations. Additional validation measurements were obtained in-vivo at 7 T. ResultsPhantom experiments showed that keeping the first gradient moment constant for all k-space lines eliminates any displacements in phase-encoding and slice-selection direction for all spoke pulses but leads to artifacts for non-zero velocity components along readout direction. Introducing variable but well-defined first gradient moments in the phase-encoding direction creates displacements along the velocity vector and thus minimizes velocity-induced geometrical distortions. Phase-resolved mean volume flow in the ascending and descending aorta obtained from two-spoke excitation showed excellent agreement with single-spoke excitation over the cardiac cycle (mean difference 0.8 ± 16.2 ml/s). ConclusionsThe use of single- and multi-spoke RF pulses for flow quantification at 7 T with controlled displacement artifacts has been successfully demonstrated. The presented techniques form the basis for correct velocity quantification and compensation not only for conventional but also for multi-spoke RF pulses allowing in-plane B1+ homogenization using parallel transmission at UHF.

Lower Oxygen Saturation in the Internal Cerebral Vein of Patients with Sickle Cell Disease Revealed By Qsm-MRI

Introduction: Imaging of brain oxygenation can be a powerful tool to assess the risk of stroke in patients with sickle cell disease (SCD), and may potentially aid in stroke prevention. Currently, there is disagreement in oxygen extraction fraction estimated by techniques that measure tissue oxygenation from deep brain structures, and techniques that measure blood oxygenation in the sagittal sinus. Quantitative susceptibility mapping (QSM) is an MRI technique that can noninvasively measure venous blood oxygen saturation (SvO2) in deep cerebral structures based on the magnetic susceptibility shift of venous blood caused by deoxyhemoglobin. In this study, our goal was to compare venous oxygen saturation measured by QSM in the internal cerebral vein (ICV) to venous oxygenation in the sagittal sinus (SS) in SCD patients, anemia subjects with normal hemoglobin and healthy controls. Methods: 16 SCD patients, 7 non-sickle anemia patients and 11 healthy controls participated in the study(Table 1). All subjects provided written consent to a protocol approved by the Committee on Clinical Investigation (CCI#11-00083). Images were acquired on a clinical 3 T Philips system with a 32-channel RF coil. The 3D gradient echo sequence had parameters: TR = 30 ms, α = 25°, 2 echoes: TE1 = 4.94 ms, ΔTE = 5.2 ms, FOV = 210 x 190 x 120 mm, spatial resolution: 0.6 x 0.6 x 1.3 mm, SENSE acceleration rate = 2 in the phase-encoding direction, BW = 289 Hz/pix and total acquisition time = 6 mins 50 seconds. Flow-compensation was added in the readout direction only, which was the anterior-posterior (AP) direction. For each subject, phase images were fitted to generate a B0 field map. Brain extraction and phase unwrapping was performed using FSL. Background field was removed using PDF. Unreliable phase voxels were identified and removed from the brain mask for subsequent processing. L1-regularized field-to-susceptibility inversion was performed to derive the susceptibility map (lambda = 4*10-4). Venous oxygen saturation was computed based on: χ = (1 - SvO2)χd-oHct + χo-wHct [1] where χ is the susceptibility measurement of the ICV, χd-o is the susceptibility shift per unit hematocrit between fully oxygenated and fully deoxygenated erythrocytes, and χo-w is the susceptibility shift between oxygenated blood cells and water. Values of 0.27 ppm (c.g.s unit) and -0.03 ppm were used for χd-o and χo-w. The ROI mask of the ICV was manually selected based on the susceptibility map that was threshold at 0.1 ppm to avoid partial-volume effect. Only the segment that had an angle below 30 degree with respect to the AP axis was included. The purpose was to exclude regions that were susceptible to flow artifact. Intra-ROI average venous oxygen saturation were computed. Group difference of venous oxygen saturation was tested using Student t-test. Results: Figure 2a shows ICV magnetic susceptibility increases linearly with hematocrit as predicted by Eq [1]. Mean ICV venous oxygen saturation in the three groups was 73.9% in CTL, 73.6% in ACTL, compared with SCD 69.3%, p&amp;lt;0.05 by analysis of variance. Figure 3 compares the ICV and SS oxygen saturation as a function of hemoglobin concentration. At normal hemoglobin values, SS saturation was lower than in the ICV, but this difference disappeared with increasing anemia severity. After controlling for hemoglobin value, there were no patient group differences in oxygen saturation in either location. Discussion: In non-anemic subjects, SS oxygen saturation was lower than in the ICV, potentially reflecting increased oxygen demand of the cerebral cortex compared with deep brain structures. ICV is one of the major deep cortical veins, and oxygen saturation of this vein can indicate the oxygenation of deep structures of the cerebral hemispheres, including the basal ganglia, corpus callosum and thalamus. With progressive anemia, brain blood flow increases to preserve cerebral oxygen delivery (not shown), but favors grey matter perfusion at the expense of deep white matter structures. This could explain the decline in venous saturation with hemoglobin seen in the ICV compared with the rise in venous saturation observed in the SS. Alternatively, global increases in flow may homogenize and shorten microvascular transit times, equalizing the oxygen extraction across brain vascular territories. Disclosures Coates: celgene: Consultancy, Honoraria, Other: steering committee of clinical study; agios pharma: Consultancy, Honoraria; vifor: Consultancy, Honoraria; apo pharma: Consultancy, Honoraria, Speakers Bureau. Wood:Apopharma: Consultancy; BiomedInformatics: Consultancy; National Institutes of Health: Research Funding; Imago Biosciences: Consultancy; WorldcareClinical: Consultancy; Philips Healthcare: Research Funding; BluebirdBio: Consultancy; Celgene: Consultancy.