Center Of K-space Research Articles

Using High-Pass Filter to Enhance Scan Specific Learning for MRI Reconstruction without Any Extra Training Data.

In accelerated MRI, the robust artificial-neural-network for k-space interpolation (RAKI) method is an attractive learning-based reconstruction that does not require additional training data. This study was focused on obtaining high quality MR images from regular under-sampled multi-coil k-space data using a high-pass filtered RAKI (HP-RAKI) reconstruction without any extra training data. MRI scan from human subjects was under-sampled with a regular pattern using skipped phase encoding and a fully sampled k-space center. A high-pass (HP) filter was applied in k-space to reduce image support to facilitate linear prediction. The HP filtered k-space center was used to train the RAKI network without any extra training data. The unacquired k-space data can be predicted from a trained RAKI network with optimized parameters. Final reconstruction was obtained after performing an inverse HP filtering for the predicted k-space data. This HP-RAKI method can be extended to corresponding residual structure (HP-rRAKI). HP-RAKI was compared with GRAPPA, HP-GRAPPA, RAKI and MW-RAKI algorithms, and HP-rRAKI was compared with corresponding residual extensions, including rRAKI and MW-rRAKI, all qualitatively and quantitatively using visual inspection and such metrics as SSIM and PSNR. HP-RAKI and HP-rRAKI were found to be effective in reconstructing MR images even at high acceleration factors. HP-RAKI and HP-rRAKI compared favorably with other algorithms. Using high-pass filtered central k-space data for training, HP-RAKI offers higher reconstruction quality for regularly under-sampled multi-coil k-space data without any extra training data. It has shown promising capabilities for fast MRI applications, especially those lacking fully sampled training data.

A hybrid PCA acceleration method for rapid real-time 2D MRI.

Objective.To develop a 2D MR acceleration method utilizing principal component analysis (PCA) in a hybrid fashion for rapid real-time applications.Approach.Retrospective testing was performed on 10 lung, 10 liver and 10 prostate 3T MRI data sets for image quality and target contourability. Sampling of k-space is performed by acquiring central (low-frequency) data in every frame while the high-frequency data is incoherently undersampled such that all of k-space is acquired in a pre-determined number of frames. Firstly, principal components (PCs) representative of intra-frame correlations between central and outer k-space data are used to estimate unsampled data in the frame of interest. Then to add further stability, PCs representative of time-domain fluctuations within a reconstruction window of the most recent frames are fit to outer k-space data (including above estimations) to obtain final estimates in the frame of interest. Accelerated reconstructions between 3x and 8x were tested for image quality and contourability along with the optimal number of PCs for fitting.Main results.It was found that at higher acceleration rates, image quality did not deteriorate significantly. Similarly, it was found that the images were of sufficient quality to contour a target using auto-contouring software at all tested acceleration rates and sites. SSIM values were found to be ⩾0.91 at all accelerations tested. Similarly dice coefficients at the different sites were found to be ⩾0.89 even at 8x accelerations which is on par with or better than intra-observer variation.Significance.This method appears to produce improved image quality and contourability compared to previous PCA methods while also allowing a greater number of PCs to be used in reconstruction. The method can be run using a simple single-channel coil and does not require significant computing power to meet real-time interventional standards (reconstruction times ∼60 ms/frame on Intel i5 CPU).

Fast, automated, real-time 3D passive balloon catheter tracking during MRI-guided cardiac catheterization using orthogonal projection imaging and real-time image-based catheter detection.

MRI-guidance of cardiac catheterization is currently performed using one or multiple 2D imaging planes, which may be suboptimal for catheter navigation, especially in patients with complex anatomies. The purpose of the work was to develop a robust real-time 3D catheter tracking method and 3D visualization strategy for improved MRI-guidance of cardiac catheterization procedures. A fast 3D tracking technique was developed using continuous acquisition of two orthogonal 2D-projection images. Each projection corresponds to a gradient echo stack of slices with only the central k-space lines being collected for each slice. To enhance catheter contrast, a saturation pulse is added ahead of the projection pair. An offline image processing algorithm was developed to identify the 2D coordinates of the balloon in each projection image and to estimate its corresponding 3D coordinates. Post-processing includes background signal suppression using an atlas of background 2D-projection images. 3D visualization of the catheter and anatomy is proposed using three live sagittal, coronal, and axial (MPR) views and 3D rendering. The technique was tested in a subset of a catheterization step in three patients undergoing MRI-guided cardiac catheterization using a passive balloon catheter. The extraction of the catheter balloon 3D coordinates was successful in all patients and for the majority of time-points (accuracy >96%). This tracking method enabled a novel 3D visualization strategy for passive balloon catheter, providing enhanced anatomical context during catheter navigation. The proposed tracking strategy shows promise for robust tracking of passive balloon catheter and may enable enhanced visualization during MRI-guided cardiac catheterization.

Dixon-based B0 self-navigation in radial stack-of-stars multi-echo gradient echo imaging.

To develop a Dixon-based self-navigation approach to estimate and correct temporal variations in radial stack-of-stars gradient echo imaging for quantitative body MRI. The proposed method estimates temporal variations using a self-navigator estimated by a graph-cut-based water-fat separation algorithm on the oversampled k-space center. The self-navigator was employed to correct for phase differences between radial spokes (one-dimensional [1D] correction) and to perform a motion-resolved reconstruction to correct spatiotemporal pseudo-periodic variations (three-dimensional [3D] correction). Numerical simulations, phantom experiments and in vivo neck scans were performed to evaluate the effects of temporal variations on the field-map, proton density fat fraction (PDFF) and map, and to validate the proposed method. Temporal variations were found to cause signal loss and phase shifts on the multi-echo images that lead to an underestimation of , while PDFF mapping was less affected. The self-navigator captured slowly varying temporal drifts and temporal variations caused by respiratory motion. While the 1D correction effectively corrected drifts in phantom studies, it was insufficient in vivo due to 3D spatially varying temporal variations with amplitudes of up to 25 Hz at 3 T near the lungs. The proposed 3D correction locally improved the correction of field-map and and reduced image artifacts. Temporal variations particularly affect mapping in radial stack-of-stars imaging. The self-navigation approach can be applied without modifying the MR acquisition to correct for drift and physiological motion-induced variations, especially in the presence of fat.

Fast and motion-robust saturation transfer MRI with inherent B0 correction using rosette trajectories and compressed sensing.

To implement rosette readout trajectories with compressed sensing reconstruction for fast and motion-robust CEST and magnetization transfer contrast imaging with inherent correction of B0 inhomogeneity. A pulse sequence was developed for fast saturation transfer imaging using a stack of rosette trajectories with a higher sampling density near the k-space center. Each rosette lobe was segmented into two halves to generate dual-echo images. B0 inhomogeneities were estimated using the phase difference between the images and corrected subsequently. The rosette-based imaging was evaluated in comparison to a fully sampled Cartesian trajectory and demonstrated on CEST phantoms (creatine solutions and egg white) and healthy volunteers at 3 T. Compared with the conventional Cartesian acquisition, compressed sensing reconstructed rosette images provided image quality with overall higher contrast-to-noise ratio and significantly faster readout time. Accurate B0 map estimation was achieved from the rosette acquisition with a negligible bias of 0.01 Hz between the rosette and dual-echo Cartesian gradient echo B0 maps, using the latter as ground truth. The water-saturation spectra (Z-spectra) and amide proton transfer weighted signals obtained from the rosette-based sequence were well preserved compared with the fully sampled data, both in the phantom and human studies. Fast, motion-robust, and inherent B0-corrected CEST and magnetization transfer contrast imaging using rosette trajectories could improve subject comfort and compliance, contrast-to-noise ratio, and provide inherent B0 homogeneity information. This work is expected to significantly accelerate the translation of CEST-MRI into a robust, clinically viable approach.

Fast 3D fMRI acquisition with high spatial resolutions over a reduced FOV.

To develop and demonstrate a fast 3D fMRI acquisition technique with high spatial resolution over a reduced FOV, named k-t 3D reduced FOV imaging (3D-rFOVI). Based on 3D gradient-echo EPI, k-t 3D-rFOVI used a 2D RF pulse to reduce the FOV in the in-plane phase-encoding direction, boosting spatial resolution without increasing echo train length. For image acceleration, full sampling was applied in the central k-space region along the through-slab direction (kz) for all time frames, while randomized undersampling was used in outer kz regions at different time frames. Images were acquired at 3T and reconstructed using a method based on partial separability. fMRI detection sensitivity of k-t 3D-rFOVI was quantitively analyzed with simulation data. Human visual fMRI experiments were performed to evaluate k-t 3D-rFOVI and compare it with a commercial multiband EPI sequence. The simulation data showed that k-t 3D-rFOVI can detect 100% of fMRI activations with an acceleration factor (R) of 2 and ˜80% with R = 6. In the human fMRI data acquired with 1.5-mm spatial resolution and 800-ms volume TR (TRvol), k-t 3D-rFOVI with R = 4 detected 46% more activated voxels in the visual cortex than the multiband EPI. Additional fMRI experiments showed that k-t 3D-rFOVI can achieve TRvol of 480 ms with R = 6, while reliably detecting visual activation. k-t 3D-rFOVI can simultaneously achieve a high spatial resolution (1.5-mm isotropically) and short TRvol (480-ms) at 3T. It offers a robust acquisition technique for fast fMRI studies over a focused brain volume.

Accelerated multi-b-value multi-shot diffusion-weighted imaging based on EPI with keyhole and a low-rank tensor constraint

PurposeMulti-Shot (MS) Echo-Planar Imaging (EPI) may improve the in-plane resolution of multi-b-value DWI, yet it also considerably increases the scan time. Here we explored the combination of EPI with Keyhole (EPIK) and a calibrationless reconstruction algorithm for acceleration of multi-b-value MS-DWI. MethodsWe firstly analyzed the impact of nonuniform phase accrual in EPIK on the reconstructed image. Based on insights gained from the analysis, we developed a calibrationless reconstruction algorithm based on a Space-Contrast-Coil Locally Low-Rank Tensor (SCC-LLRT) constraint for reconstruction of EPIK-acquired data. We compared the algorithm with a modified SPatial-Angular Locally Low-Rank (SPA-LLR) algorithm through simulations, phantoms, and in vivo study. We then compared EPIK with uniformly undersampled EPI for accelerating multi-b-value DWI in 6 healthy subjects. ResultsThrough theoretical derivations, we found that the reconstruction of EPIK with a SENSE-encoding-based algorithm, such as SPA-LLR, may cause additional aliasing artifacts due to the frequency-dependent distortion of the coil sensitivity. Results from simulations, phantoms, and in vivo study verified the theoretical finding by showing that the calibrationless SCC-LLRT algorithm reduced aliasing artifacts compared with SPA-LLR. Finally, EPIK with SCC-LLRT substantially reduced the ghosting artifacts compared with uniform undersampled multi-b-value DWI, decreasing the fitting errors in ADC (0.05 ± 0.01 vs 0.10 ± 0.01, P < 0.001) and IVIM mapping (0.026 ± 0.004 vs 0.06 ± 0.006, P < 0.001). ConclusionThe SCC-LLRT algorithm reduced the aliasing artifacts of EPIK by using a calibrationless modeling of the multi-coil data. The dense sampling of k-space center offers EPIK a potential to improve image quality for acceleration of multi-b-value MS-DWI.

A characterization of cardiac-induced noise in R2 * maps ofthe brain.

Cardiac pulsation increases the noise level in brain maps of the transverse relaxation rate R2 *. Cardiac-induced noise is challenging to mitigate during the acquisition of R2 * mapping data because its characteristics are unknown. In this work, we aim to characterize cardiac-induced noise in brain maps of the MRI parameter R2 *. We designed a sampling strategy to acquire multi-echo 3D data in 12 intervals of the cardiac cycle, monitored with a fingertip pulse-oximeter. We measured the amplitude of cardiac-induced noise in this data and assessed the effect of cardiac pulsation on R2 * maps computed across echoes. The area of k-space that contains most of the cardiac-induced noise in R2 * maps was then identified. Based on these characteristics, we introduced a tentative sampling strategy that aims to mitigate cardiac-induced noise in R2 * maps of the brain. In inferior brain regions, cardiac pulsation accounts for R2 * variations of up to 3 s-1 across the cardiac cycle (i.e., ∼35% of the overall variability). Cardiac-induced fluctuations occur throughout the cardiac cycle, with a reduced intensity during the first quarter of the cycle. A total of 50% to 60% of the overall cardiac-induced noise is localized near the k-space center (k < 0.074 mm-1 ). The tentative cardiac noise mitigation strategy reduced the variability of R2 * maps across repetitions by 11% in the brainstem and 6% across the whole brain. We provide a characterization of cardiac-induced noise in brain R2 * maps that can be used as a basis for the design of mitigation strategies during data acquisition.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction. Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF), , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and between and QSM in in vivo study. Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction ( = 0.31 0.05, = 0.999) and Cartesian without motion ( = 0.32 0.04, = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit ( = 0.00261 0.524, = 0.977) compared to motion-averaged reconstruction ( = 0.0021 0.572, = 0.723) in ROI-based linear regression analysis between and QSM. Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Retrospective frequency drift correction of rosette MRSI data using spectral registration.

Frequency drift correction is an important postprocessing step in MRS that yields improvements in spectral quality and metabolite quantification. Although routinely applied in single-voxel MRS, drift correction is much more challenging in MRSI due to the presence of phase-encoding gradients. Thus, separately acquired navigator scans are normally required for drift estimation. In this work, we demonstrate the use of self-navigating rosette MRSI trajectories combined with time-domain spectral registration to enable retrospective frequency drift corrections without the need for separately acquired navigator echoes. A rosette MRSI sequence was implemented to acquire data from the brains of 5 healthy volunteers. FIDs from the center of k-space ( FIDs) were isolated from each shot of the rosette acquisition, and time-domain spectral registration was used to estimate the frequency offset of each FID relative to a reference scan (the first FID in the series). The estimated frequency offsets were then used to apply corrections throughout -space. Improvements in spectral quality were assessed before and after drift correction. Spectral registration resulted in significant improvements in signal-to-noise ratio(12.9%) and spectral linewidths (18.5%). Metabolite quantification was performed using LCModel, and the average Cramer-Rao lower bounds uncertainty estimates were reduced by 5.0% for all metabolites, following field drift correction. This study demonstrated the use of self-navigating rosette MRSI trajectories to retrospectively correct frequency drift errors in in vivo MRSI data. This correction yields meaningful improvements in spectral quality.

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.

Compartment-based reconstruction of 3D acquisition-weighted 31 P cardiac magnetic resonance spectroscopic imaging at 7 T: A reproducibility study.

Even at 7 T, cardiac 31 P magnetic resonance spectroscopic imaging (MRSI) is fundamentally limited by low signal-to-noise ratio (SNR), leading to long scan times and poor temporal and spatial resolutions. Compartment-based reconstruction algorithms such as magnetic resonance spectroscopy with linear algebraic modeling (SLAM) and spectral localization by imaging (SLIM) may improve SNR or reduce scan time without changes to acquisition. Here, we compare the repeatability and SNR performance of these compartment-based methods, applied to three different acquisition schemes at 7 T. Twelve healthy volunteers were scanned twice. Each scan session consisted of a 6.5-min 3D acquisition-weighted (AW) cardiac 31 P phase encode-based MRSI acquisition and two 6.5-min truncated k-space acquisitions with increased averaging (4 × 4 × 4 central k-space phase encodes and fractional SLAM [fSLAM] optimized k-space phase encodes). Spectra were reconstructed using (i) AW Fourier reconstruction; (ii) AW SLAM; (iii) AW SLIM; (iv) 4 × 4 × 4 SLAM; (v) 4 × 4 × 4 SLIM; and (vi) fSLAM acquisition-reconstruction combinations. The phosphocreatine-to-adenosine triphosphate (PCr/ATP) ratio, the PCr SNR, and spatial response functions were computed, in addition to coefficients of reproducibility and variability. Using the compartment-based reconstruction algorithms with the AW 31 P acquisition resulted in a significant increase in SNR compared with previously published Fourier-based MRSI reconstruction methods while maintaining the measured PCr/ATP ratio and improving interscan reproducibility. The alternative acquisition strategies with truncated k-space performed no better than the common AW approach. Compartment-based spectroscopy approaches provide an attractive reconstruction method for cardiac 31 P spectroscopy at 7 T, improving reproducibility and SNR without the need for a dedicated k-space sampling strategy.

Joint Cardiac T1 Mapping and Cardiac Cine Using Manifold Modeling

The main focus of this work is to introduce a single free-breathing and ungated imaging protocol to jointly estimate cardiac function and myocardial T1 maps. We reconstruct a time series of images corresponding to k-space data from a free-breathing and ungated inversion recovery gradient echo sequence using a manifold algorithm. We model each image in the time series as a non-linear function of three variables: cardiac and respiratory phases and inversion time. The non-linear function is realized using a convolutional neural networks (CNN) generator, while the CNN parameters, as well as the phase information, are estimated from the measured k-t space data. We use a dense conditional auto-encoder to estimate the cardiac and respiratory phases from the central multi-channel k-space samples acquired at each frame. The latent vectors of the auto-encoder are constrained to be bandlimited functions with appropriate frequency bands, which enables the disentanglement of the latent vectors into cardiac and respiratory phases, even when the data are acquired with intermittent inversion pulses. Once the phases are estimated, we pose the image recovery as the learning of the parameters of the CNN generator from the measured k-t space data. The learned CNN generator is used to generate synthetic data on demand by feeding it with appropriate latent vectors. The proposed approach capitalizes on the synergies between cine MRI and T1 mapping to reduce the scan time and improve patient comfort. The framework also enables the generation of synthetic breath-held cine movies with different inversion contrasts, which improves the visualization of the myocardium. In addition, the approach also enables the estimation of the T1 maps with specific phases, which is challenging with breath-held approaches.

Calibrationless reconstruction of uniformly-undersampled multi-channel MR data with deep learning estimated ESPIRiT maps.

To develop a truly calibrationless reconstruction method that derives An Eigenvalue Approach to Autocalibrating Parallel MRI (ESPIRiT) maps from uniformly-undersampled multi-channel MR data by deep learning. ESPIRiT, one commonly used parallel imaging reconstruction technique, forms the images from undersampled MR k-space data using ESPIRiT maps that effectively represents coil sensitivity information. Accurate ESPIRiT map estimation requires quality coil sensitivity calibration or autocalibration data. We present a U-Net based deep learning model to estimate the multi-channel ESPIRiT maps directly from uniformly-undersampled multi-channel multi-slice MR data. The model is trained using fully-sampled multi-slice axial brain datasets from the same MR receiving coil system. To utilize subject-coil geometric parameters available for each dataset, the training imposes a hybrid loss on ESPIRiT maps at the original locations as well as their corresponding locations within the standard reference multi-slice axial stack. The performance of the approach was evaluated using publicly available T1-weighed brain and cardiac data. The proposed model robustly predicted multi-channel ESPIRiT maps from uniformly-undersampled k-space data. They were highly comparable to the reference ESPIRiT maps directly computed from 24 consecutive central k-space lines. Further, they led to excellent ESPIRiT reconstruction performance even at high acceleration, exhibiting a similar level of errors and artifacts to that by using reference ESPIRiT maps. A new deep learning approach is developed to estimate ESPIRiT maps directly from uniformly-undersampled MR data. It presents a general strategy for calibrationless parallel imaging reconstruction through learning from the coil and protocol-specific data.

Cartesian vs radial MR-STAT: An efficiency and robustness study

MR Spin TomogrAphy in Time-domain (“MR-STAT”) is quantitative MR technique in which multiple quantitative parameters are estimated from a single short scan by solving a large-scale non-linear optimization problem. In this work we extended the MR-STAT framework to non-Cartesian gradient trajectories. Cartesian MR-STAT and radial MR-STAT were compared in terms of time-efficiency and robustness in simulations, gel phantom measurements and in vivo measurements.In simulations, we observed that both Cartesian and radial MR-STAT are highly robust against undersampling. Radial MR-STAT does have a lower spatial encoding power because the outer corners of k-space are never sampled. However, especially in T2, this is compensated by a higher dynamic encoding power that comes from sampling the k-space center with each readout. In gel phantom measurements, Cartesian MR-STAT was observed to be robust against overfitting whereas radial MR-STAT suffered from high-frequency artefacts in the parameter maps at later iterations. These artefacts are hypothesized to be related to hardware imperfections and were (partially) suppressed with image filters. The time-efficiencies were higher for Cartesian MR-STAT in all vials. In-vivo, the radial reconstruction again suffered from overfitting artefacts. The robustness of Cartesian MR-STAT over the entire range of experiments may make it preferable in a clinical setting, despite radial MR-STAT resulting in a higher T1 time-efficiency in white matter.

Improved Parallel Magnetic Resonance Imaging Reconstruction with Sampling based on Normal Distribution

Generalized autocalibrating partially parallel acquisitions (GRAPPA) is a k-space-based magnetic resonance imaging (MRI) reconstruction algorithm that obtains coefficients by acquiring k-space center data with the Nyquist frequency requirement, fitting these data as auto-calibration signal (ACS). The under-sampled k-space data are then filled using the linear correlation of the neighboring points of the k-space. At high reduction factors, images reconstructed by using GRAPPA have high levels of noise even when the number of ACS lines is significantly increased. NL-GRAPPA (Non-linear GRAPPA) can improve GRAPPA for a comparatively impressive image quality. However, it is extremely difficult to reduce the sampling time significantly because of the high number of ACS lines that must be sampled. In this paper, we proposed an undersampling parallel MRI (pMRI) reconstruction technique based on the ideology of normal distribution to improve the traditional GRAPPA and NL-GRAPPA methods. The proposed method mainly takes advantage of the fact that the data in the center of the k-space had a greater impact on the imaging quality than the data at the edges, and then samples the data based on the ideology of normal distribution, i.e., discarding part of the data at the edges of the k-space. The experimental results show that our proposed method can significantly reduce the number of sampled lines and have higher imaging quality with a lower number of ACS lines.

Free-breathing myocardial T1 mapping using inversion-recovery radial FLASH and motion-resolved model-based reconstruction.

To develop a free-breathing myocardial mapping technique using inversion-recovery (IR) radial fast low-angle shot (FLASH) and calibrationless motion-resolved model-based reconstruction. Free-running (free-breathing, retrospective cardiac gating) IR radial FLASH is used for data acquisition at 3T. First, to reduce the waiting time between inversions, an analytical formula is derived that takes the incomplete recovery into account for an accurate calculation. Second, the respiratory motion signal is estimated from the k-space center of the contrast varying acquisition using an adapted singular spectrum analysis (SSA-FARY) technique. Third, a motion-resolved model-based reconstruction is used to estimate both parameter and coil sensitivity maps directly from the sorted k-space data. Thus, spatiotemporal total variation, in addition to the spatial sparsity constraints, can be directly applied to the parameter maps. Validations are performed on an experimental phantom, 11 human subjects, and a young landrace pig with myocardial infarction. In comparison to an IR spin-echo reference, phantom results confirm good accuracy, when reducing the waiting time from 5 s to 1 s using the new correction. The motion-resolved model-based reconstruction further improves precision compared to the spatial regularization-only reconstruction. Aside from showing that a reliable respiratory motion signal can be estimated using modified SSA-FARY, in vivo studies demonstrate that dynamic myocardial maps can be obtained within 2 min with good precision and repeatability. Motion-resolved myocardial mapping during free-breathing with good accuracy, precision and repeatability can be achieved by combining inversion-recovery radial FLASH, self-gating and a calibrationless motion-resolved model-based reconstruction.

Single projection driven real-time multi-contrast (SPIDERM) MR imaging using pre-learned spatial subspace and linear transformation

Objective. To develop and test the feasibility of a novel Single ProjectIon DrivEn Real-time Multi-contrast (SPIDERM) MR imaging technique that can generate real-time 3D images on-the-fly with flexible contrast weightings and a low latency. Approach. In SPIDERM, a ‘prep’ scan is first performed, with sparse k-space sampling periodically interleaved with the central k-space line (navigator data), to learn a subject-specific model, incorporating a spatial subspace and a linear transformation between navigator data and subspace coordinates. A ‘live’ scan is then performed by repeatedly acquiring the central k-space line only to dynamically determine subspace coordinates. With the ‘prep’-learned subspace and ‘live’ coordinates, real-time 3D images are generated on-the-fly with computationally efficient matrix multiplication. When implemented based on a multi-contrast pulse sequence, SPIDERM further allows for data-driven image contrast regeneration to convert real-time contrast-varying images into contrast-frozen images at user’s discretion while maintaining motion states. Both digital phantom and in-vivo experiments were performed to evaluate the technical feasibility of SPIDERM. Main results. The elapsed time from the input of the central k-space line to the generation of real-time contrast-frozen 3D images was approximately 45 ms, permitting a latency of 55 ms or less. Motion displacement measured from SPIDERM and reference images showed excellent correlation (). Geometric variation from the ground truth in the digital phantom was acceptable as demonstrated by pancreas contour analysis (Dice ≥ 0.84, mean surface distance ≤ 0.95 mm). Quantitative image quality metrics showed good consistency between reference images and contrast-varying SPIDREM images in in-vivo studies (mean ). Significance. SPIDERM is capable of generating real-time multi-contrast 3D images with a low latency. An imaging framework based on SPIDERM has the potential to serve as a standalone package for MR-guided radiation therapy by offering adaptive simulation through a ‘prep’ scan and real-time image guidance through a ‘live’ scan.

Improved signal-to-noise performance of MultiNet GRAPPA 1 H FID MRSI reconstruction with semi-synthetic calibration data.

To further develop MultiNet GRAPPA, a neural-network-based reconstruction, for lower SNR proton MRSI (1 H MRSI) data using adapted undersampling schemes and improved training sets. 1 H FID-MRSI data and an anatomical image for GRAPPA reconstruction were acquired in two slices in the human brain (n = 6) at 7T. MRSI data were retrospectively undersampled for a 4×, 6×, and 7× acceleration rate. Signal-to-noise, relative error (RE) between accelerated and fully sampled metabolic maps, RMS of the lipid artifacts, and fitting reliability were compared across acceleration rates, to the fully sampled data, and with different kinds and amounts of training images. Training with semi-synthetic images resulted in higher SNR and lower lipid RMS relative to training with acquired images from one or several subjects. SNR increased with the number of semi-synthetic training images and the 4× accelerated data retains ∼30% more SNR than other accelerated data. Spectra reconstructed with 20 semi-synthetic averages retained ∼100% more SNR and had ∼5% lower lipid RMS than those reconstructed with the center k-space points of one image as was originally proposed for very high SNR MRSI data and had higher fitting reliability. The metabolite RE was lowest when training with 20-semi-synthetic training images and highest when training with the center k-space points of one image. MultiNet GRAPPA is feasible with lower SNR 1 H MRSI data if 20-semi-synthetic training images are used at a 4× acceleration rate. This acceleration rate provided the best trade-off between scan time and spectral SNR.

Development of a novel 10-echo multi-contrast sequence based on EPIK to deliver simultaneous quantification of T2 and T2* with application to oxygen extraction fraction.

The simultaneous quantification of T2 and T2 * maps based on fast sequences for combined GE and SE acquisition has rekindled research and clinical interest by offering a wide range of attractive applications, e.g., dynamic tracking of oxygen extraction fraction (OEF). However, previously published methods based on EPI-readouts have been hindered by resolution and the number of acquired echoes. This work presents a novel 10-echo GE-SE EPIK (EPI with keyhole) sequence for the rapid quantification of T2 '. T2 /T2 * maps from the GE-SE EPIK sequence were validated using three phantoms and 15 volunteers at 3T. The incorporation of a sliding window approach, combined with the full sampling of the k-space center inherent to EPIK, enables a high effective temporal resolution. That is, for an eight-slice breath-hold experiment, a temporal sampling rate of eight reconstructed slices per 1.1s. In comparison with repeated single-echo SE, multi-echo GE, and spectroscopy methods, the GE-SE EPIK sequence shows good agreement in quantifying T2 /T2 * values, while the gray matter/white matter separation yielded the expected contrast differentiation. The OEF was calculated with a view to an initial application with clinical relevance, producing results comparable to those in the literature and with good sensitivity in breath-hold experiments. GE-SE EPIK provides increased resolution and more echoes, including two SEs, than comparable sequences. Moreover, GE-SE EPIK achieves this within an acquisition time of 57 s for 20 slices (matrix size=128×128; FOV=24 cm) and with a reasonably short TE for the final echo (114 ms). The sequence can dynamically track OEF changes in a breath-hold experiment.