Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging.

PURPOSE: To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques. METHODS: Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method. RESULTS: The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m. CONCLUSION: A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

Dynamic mode decomposition of dynamic MRI for assessment of pulmonary ventilation and perfusion.

PURPOSE: To introduce dynamic mode decomposition (DMD) as a robust alternative for the assessment of pulmonary functional information from dynamic non-contrast-enhanced acquisitions. METHODS: Pulmonary fractional ventilation and normalized perfusion maps were obtained using DMD from simulated phantoms as well as in vivo dynamic acquisitions of healthy volunteers at 1.5T. The performance of DMD was compared with conventional Fourier decomposition (FD) and matrix pencil (MP) methods in estimating functional map values. The proposed method was evaluated based on estimated signal amplitude in functional maps across varying number of measurements. RESULTS: Quantitative assessments performed on phantoms and in vivo measurements indicate that DMD is capable of successfully obtaining pulmonary functional maps. Specifically, compared to FD and MP methods, DMD is able to reduce variations in estimated amplitudes across different number of measurements. This improvement is evident in the fractional ventilation and normalized perfusion maps obtain from phantom simulations with frequency variations and noise, as well as in the maps obtained from in vivo measurements. CONCLUSIONS: A robust method for accurately estimating pulmonary ventilation and perfusion related signal changes in dynamic acquisitions is presented. The proposed method uses DMD to obtain functional maps reliably, while reducing amplitude variations caused by differences in number of measurements.

Functional lung imaging of 2-year-old children after congenital diaphragmatic hernia repair using dynamic mode decomposition MRI.

OBJECTIVES: To investigate the feasibility of non-contrast-enhanced functional lung imaging in 2-year-old children after congenital diaphragmatic hernia (CDH) repair. METHODS: Fifteen patients after CDH repair were examined using non-contrast-enhanced dynamic magnetic resonance imaging (MRI). For imaging two protocols were used during free-breathing: Protocol A with high temporal resolution and Protocol B with high spatial resolution. The dynamic images were then analysed through a recently developed post-processing method called dynamic mode decomposition (DMD) to obtain ventilation and perfusion maps. The ventilation ratios (VRatio) and perfusion ratios (QRatio) of ipsilateral to contralateral lung were compared to evaluate functional differences. Lastly, DMD MRI-based perfusion results were compared with perfusion parameters obtained using dynamic contrast-enhanced (DCE) MRI to assess agreement between methods. RESULTS: Both imaging protocols successfully generated pulmonary ventilation (V) and perfusion (Q) maps in all patients. Overall, the VRatio and QRatio values were 0.84 ± 0.19 and 0.70 ± 0.24 for Protocol A, and 0.88 ± 0.18 and 0.72 ± 0.23 for Protocol B, indicating reduced ventilation ([Formula: see text]) and perfusion ([Formula: see text]) on the ipsilateral side. Moreover, there is a very strong positive correlation ([Formula: see text]) and close agreement between DMD MRI-based perfusion values and DCE MRI-based perfusion parameters. CONCLUSIONS: DMD MRI can obtain pulmonary functional information in 2-year-old CDH patients. The results obtained with DMD MRI correlate with DCE MRI, without the need for ionising radiation or exposure to contrast agents. While further studies with larger cohorts are warranted, DMD MRI is a promising option for functional lung imaging in CDH patients. CLINICAL RELEVANCE STATEMENT: We demonstrate that pulmonary ventilation and perfusion information can be obtained in 2-year-old patients after CDH repair, without the need for ionising radiation or contrast agents by utilising non-contrast-enhanced MRI acquisitions together with dynamic mode decomposition analysis. KEY POINTS: • Non-contrast-enhanced functional MR imaging is a promising option for functional lung imaging in 2-year-old children after congenital diaphragmatic hernia. • DMD MRI can generate pulmonary ventilation and perfusion maps from free-breathing dynamic acquisitions without the need for ionising radiation or contrast agents. • Lung perfusion parameters obtained with DMD MRI correlate with perfusion parameters obtained using dynamic contrast-enhanced MRI.

Phase-cycled balanced SSFP imaging for non-contrast-enhanced functional lung imaging.

PURPOSE: To introduce phase-cycled balanced SSFP (bSSFP) acquisition as an alternative in Fourier decomposition MRI for improved robustness against field inhomogeneities. METHODS: Series 2D dynamic lung images were acquired in 5 healthy volunteers at 1.5 T and 3 T using bSSFP sequence with multiple RF phase increments and compared with conventional single RF phase increment acquisitions. The approach was evaluated based on functional map homogeneity analysis, while ensuring image and functional map quality by means of SNR and contrast-to-noise ratio analyses. RESULTS: At both field strengths, functional maps obtained with phase-cycled acquisitions displayed improved robustness against local signal losses compared with single-phase acquisitions. The coefficient of variation (mean ± SD, across volunteers) measured in the ventilation maps resulted in 29.7 ± 2.6 at 1.5 T and 37.5 ± 3.1 at 3 T for phase-cycled acquisitions, compared with 39.9 ± 5.2 at 1.5 T and 49.5 ± 3.7 at 3 T for single-phase acquisitions, indicating a significant improvement ( p<0.05$$ p<0.05 $$ ) in ventilation map homogeneity. CONCLUSIONS: Phase-cycled bSSFP acquisitions improve robustness against field inhomogeneity artifacts and significantly improve ventilation map homogeneity at both field strengths. As such, phase-cycled bSSFP may serve as a robust alternative in lung function assessments.

The number of glomeruli and pyruvate metabolism is not strongly coupled in the healthy rat kidney.

PURPOSE: The number of glomeruli is different in men and women, as they also present different prevalence and progression of chronic kidney disease. A recent study has demonstrated a potential difference in renal metabolism between sexes, and a potential explanation could be the differences in glomeruli number. This study investigates the potential correlation between glomerular number and pyruvate metabolism in healthy kidneys. METHODS: This study is an experimental study with rats (N = 12). We used cationized-ferritin MRI to visualize and count glomeruli and hyperpolarized [1-13 C]pyruvate to map the metabolism. Dynamic contrast-enhanced MRI was used to analyze kidney hemodynamics using gadolinium tracer. RESULTS: Data showed no or subtle correlation between the number of glomeruli and the pyruvate metabolism. Minor differences were observed in the number of glomeruli (female = 24,509 vs. male = 26 350; p = .16), renal plasma flow (female = 606.6 vs. male= 455.7 ml/min/100 g; p = .18), and volume of distribution (female = 87.44 vs. male = 76.61 ml/100 ml; p = .54) between sexes. Mean transit time was significantly prolonged in males compared with females (female = 8.868 s vs. male = 10.63 s; p = .04). CONCLUSION: No strong statistically significant correlation between the number of glomeruli and the pyruvate metabolism was found in healthy rat kidneys.

Reconstruction by calibration over tensors for multi-coil multi-acquisition balanced SSFP imaging.

PURPOSE: To develop a rapid imaging framework for balanced steady-state free precession (bSSFP) that jointly reconstructs undersampled data (by a factor of R) across multiple coils (D) and multiple acquisitions (N). To devise a multi-acquisition coil compression technique for improved computational efficiency. METHODS: The bSSFP image for a given coil and acquisition is modeled to be modulated by a coil sensitivity and a bSSFP profile. The proposed reconstruction by calibration over tensors (ReCat) recovers missing data by tensor interpolation over the coil and acquisition dimensions. Coil compression is achieved using a new method based on multilinear singular value decomposition (MLCC). ReCat is compared with iterative self-consistent parallel imaging (SPIRiT) and profile encoding (PE-SSFP) reconstructions. RESULTS: Compared to parallel imaging or profile-encoding methods, ReCat attains sensitive depiction of high-spatial-frequency information even at higher R. In the brain, ReCat improves peak SNR (PSNR) by 1.1 ± 1.0 dB over SPIRiT and by 0.9 ± 0.3 dB over PE-SSFP (mean ± SD across subjects; average for N = 2-8, R = 8-16). Furthermore, reconstructions based on MLCC achieve 0.8 ± 0.6 dB higher PSNR compared to those based on geometric coil compression (GCC) (average for N = 2-8, R = 4-16). CONCLUSION: ReCat is a promising acceleration framework for banding-artifact-free bSSFP imaging with high image quality; and MLCC offers improved computational efficiency for tensor-based reconstructions. Magn Reson Med 79:2542-2554, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Profile-encoding reconstruction for multiple-acquisition balanced steady-state free precession imaging.

PURPOSE: The scan-efficiency in multiple-acquisition balanced steady-state free precession imaging can be maintained by accelerating and reconstructing each phase-cycled acquisition individually, but this strategy ignores correlated structural information among acquisitions. Here, an improved acceleration framework is proposed that jointly processes undersampled data across N phase cycles. METHODS: Phase-cycled imaging is cast as a profile-encoding problem, modeling each image as an artifact-free image multiplied with a distinct balanced steady-state free precession profile. A profile-encoding reconstruction (PE-SSFP) is employed to recover missing data by enforcing joint sparsity and total-variation penalties across phase cycles. PE-SSFP is compared with individual compressed-sensing and parallel-imaging (ESPIRiT) reconstructions. RESULTS: In the brain and the knee, PE-SSFP yields improved image quality compared to individual compressed-sensing and other tested methods particularly for higher N values. On average, PE-SSFP improves peak SNR by 3.8 ± 3.0 dB (mean ± s.e. across N = 2-8) and structural similarity by 1.4 ± 1.2% over individual compressed-sensing, and peak SNR by 5.6 ± 0.7 dB and structural similarity by 7.1 ± 0.5% over ESPIRiT. CONCLUSION: PE-SSFP attains improved image quality and preservation of high-spatial-frequency information at high acceleration factors, compared to conventional reconstructions. PE-SSFP is a promising technique for scan-efficient balanced steady-state free precession imaging with improved reliability against field inhomogeneity. Magn Reson Med 78:1316-1329, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Targeted vessel reconstruction in non-contrast-enhanced steady-state free precession angiography.

Image quality in non-contrast-enhanced (NCE) angiograms is often limited by scan time constraints. An effective solution is to undersample angiographic acquisitions and to recover vessel images with penalized reconstructions. However, conventional methods leverage penalty terms with uniform spatial weighting, which typically yield insufficient suppression of aliasing interference and suboptimal blood/background contrast. Here we propose a two-stage strategy where a tractographic segmentation is employed to auto-extract vasculature maps from undersampled data. These maps are then used to incur spatially adaptive sparsity penalties on vascular and background regions. In vivo steady-state free precession angiograms were acquired in the hand, lower leg and foot. Compared with regular non-adaptive compressed sensing (CS) reconstructions (CSlow ), the proposed strategy improves blood/background contrast by 71.3 ± 28.9% in the hand (mean ± s.d. across acceleration factors 1-8), 30.6 ± 11.3% in the lower leg and 28.1 ± 7.0% in the foot (signed-rank test, P < 0.05 at each acceleration). The proposed targeted reconstruction can relax trade-offs between image contrast, resolution and scan efficiency without compromising vessel depiction.

Automated Parameter Selection for Accelerated MRI Reconstruction via Low-Rank Modeling of Local k-Space Neighborhoods.

PURPOSE: Image quality in accelerated MRI rests on careful selection of various reconstruction parameters. A common yet tedious and error-prone practice is to hand-tune each parameter to attain visually appealing reconstructions. Here, we propose a parameter tuning strategy to automate hybrid parallel imaging (PI) - compressed sensing (CS) reconstructions via low-rank modeling of local k-space neighborhoods (LORAKS) supplemented with sparsity regularization in wavelet and total variation (TV) domains. METHODS: For low-rank regularization, we leverage a soft-thresholding operation based on singular values for matrix rank selection in LORAKS. For sparsity regularization, we employ Stein's unbiased risk estimate criterion to select the wavelet regularization parameter and local standard deviation of reconstructions to select the TV regularization parameter. Comprehensive demonstrations are presented on a numerical brain phantom and in vivo brain and knee acquisitions. Quantitative assessments are performed via PSNR, SSIM and NMSE metrics. RESULTS: The proposed hybrid PI-CS method improves reconstruction quality compared to PI-only techniques, and it achieves on par image quality to reconstructions with brute-force optimization of reconstruction parameters. These results are prominent across several different datasets and the range of examined acceleration rates. CONCLUSION: A data-driven parameter tuning strategy to automate hybrid PI-CS reconstructions is presented. The proposed method achieves reliable reconstructions of accelerated multi-coil MRI datasets without the need for exhaustive hand-tuning of reconstruction parameters.

Projection onto Epigraph Sets for Rapid Self-Tuning Compressed Sensing MRI.

The compressed sensing (CS) framework leverages the sparsity of MR images to reconstruct from the undersampled acquisitions. CS reconstructions involve one or more regularization parameters that weigh sparsity in transform domains against fidelity to acquired data. While parameter selection is critical for reconstruction quality, the optimal parameters are subject and dataset specific. Thus, commonly practiced heuristic parameter selection generalizes poorly to independent datasets. Recent studies have proposed to tune parameters by estimating the risk of removing significant image coefficients. Line searches are performed across the parameter space to identify the parameter value that minimizes this risk. Although effective, these line searches yield prolonged reconstruction times. Here, we propose a new self-tuning CS method that uses computationally efficient projections onto epigraph sets of the l1 and total-variation norms to simultaneously achieve parameter selection and regularization. In vivo demonstrations are provided for balanced steady-state free precession, time-of-flight, and T1-weighted imaging. The proposed method achieves an order of magnitude improvement in computational efficiency over line-search methods while maintaining near-optimal parameter selection.

Learning-Based Compressive MRI.

In the area of magnetic resonance imaging (MRI), an extensive range of non-linear reconstruction algorithms has been proposed which can be used with general Fourier subsampling patterns. However, the design of these subsampling patterns has typically been considered in isolation from the reconstruction rule and the anatomy under consideration. In this paper, we propose a learning-based framework for optimizing MRI subsampling patterns for a specific reconstruction rule and anatomy, considering both the noiseless and noisy settings. Our learning algorithm has access to a representative set of training signals, and searches for a sampling pattern that performs well on average for the signals in this set. We present a novel parameter-free greedy mask selection method and show it to be effective for a variety of reconstruction rules and performance metrics. Moreover, we also support our numerical findings by providing a rigorous justification of our framework via statistical learning theory.