Motion-compensated image reconstruction for improved kidney function assessment using dynamic contrast-enhanced MRI.

Accurately measuring renal function is crucial for pediatric patients with kidney conditions. Traditional methods have limitations, but dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a safe and efficient approach for detailed anatomical evaluation and renal function assessment. However, motion artifacts during DCE-MRI can degrade image quality and introduce misalignments, leading to unreliable results. This study introduces a motion-compensated reconstruction technique for DCE-MRI data acquired using golden-angle radial sampling. Our proposed method achieves three key objectives: (1) identifying and removing corrupted data (outliers) using a Gaussian process model fitting with a k -space center navigator, (2) efficiently clustering the data into motion phases and performing interphase registration, and (3) utilizing a novel formulation of motion-compensated radial reconstruction. We applied the proposed motion correction (MoCo) method to DCE-MRI data affected by varying degrees of motion, including both respiratory and bulk motion. We compared the outcomes with those obtained from the conventional radial reconstruction. Our evaluation encompassed assessing the quality of images, concentration curves, and tracer kinetic model fitting, and estimating renal function. The proposed MoCo reconstruction improved the temporal signal-to-noise ratio for all subjects, with a 21.8% increase on average, while total variation values of the aorta, right, and left kidney concentration were improved for each subject, with 32.5%, 41.3%, and 42.9% increases on average, respectively. Furthermore, evaluation of tracer kinetic model fitting indicated that the median standard deviation of the estimated filtration rate ( σ F T ), mean normalized root-mean-squared error (nRMSE), and chi-square goodness-of-fit of tracer kinetic model fit were decreased from 0.10 to 0.04, 0.27 to 0.24, and, 0.43 to 0.27, respectively. The proposed MoCo technique enabled more reliable renal function assessment and improved image quality for detailed anatomical evaluation in the case of bulk and respiratory motion during the acquisition of DCE-MRI.

Prospective motion correction in kidney MRI using FID navigators.

PURPOSE: Abdominal MRI scans may require breath-holding to prevent image quality degradation, which can be challenging for patients, especially children. In this study, we evaluate whether FID navigators can be used to measure and correct for motion prospectively, in real-time. METHODS: FID navigators were inserted into a 3D radial sequence with stack-of-stars sampling. MRI experiments were conducted on 6 healthy volunteers. A calibration scan was first acquired to create a linear motion model that estimates the kidney displacement due to respiration from the FID navigator signal. This model was then applied to predict and prospectively correct for motion in real time during deep and continuous deep breathing scans. Resultant images acquired with the proposed technique were compared with those acquired without motion correction. Dice scores were calculated between inhale/exhale motion states. Furthermore, images acquired using the proposed technique were compared with images from extra-dimensional golden-angle radial sparse parallel, a retrospective motion state binning technique. RESULTS: Images reconstructed for each motion state show that the kidneys' position could be accurately tracked and corrected with the proposed method. The mean of Dice scores computed between the motion states were improved from 0.93 to 0.96 using the proposed technique. Depiction of the kidneys was improved in the combined images of all motion states. Comparing results of the proposed technique and extra-dimensional golden-angle radial sparse parallel, high-quality images can be reconstructed from a fraction of spokes using the proposed method. CONCLUSION: The proposed technique reduces blurriness and motion artifacts in kidney imaging by prospectively correcting their position both in-plane and through-slice.

Learning the Regularization in DCE-MR Image Reconstruction for Functional Imaging of Kidneys.

Kidney DCE-MRI aims at both qualitative assessment of kidney anatomy and quantitative assessment of kidney function by estimating the tracer kinetic (TK) model parameters. Accurate estimation of TK model parameters requires an accurate measurement of the arterial input function (AIF) with high temporal resolution. Accelerated imaging is used to achieve high temporal resolution, which yields under-sampling artifacts in the reconstructed images. Compressed sensing (CS) methods offer a variety of reconstruction options. Most commonly, sparsity of temporal differences is encouraged for regularization to reduce artifacts. Increasing regularization in CS methods removes the ambient artifacts but also over-smooths the signal temporally which reduces the parameter estimation accuracy. In this work, we propose a single image trained deep neural network to reduce MRI under-sampling artifacts without reducing the accuracy of functional imaging markers. Instead of regularizing with a penalty term in optimization, we promote regularization by generating images from a lower dimensional representation. In this manuscript we motivate and explain the lower dimensional input design. We compare our approach to CS reconstructions with multiple regularization weights. Proposed approach results in kidney biomarkers that are highly correlated with the ground truth markers estimated using the CS reconstruction which was optimized for functional analysis. At the same time, the proposed approach reduces the artifacts in the reconstructed images.

SNR Weighting for Shear Wave Speed Reconstruction in Tomoelastography.

In tomoelastography, to achieve a final wave speed map by combining reconstructions obtained from all spatial directions and excitation frequencies, the use of weights is inevitable. Here, a new weighting scheme, which maximizes the signal-to-noise ratio (SNR) of the final wave speed map, has been proposed. To maximize the SNR of the final wave speed map, the use of squares of estimated SNR values of reconstructed individual maps has been proposed. Therefore, derivations of the SNR of the reconstructed wave speed maps have become necessary. Considering the noise on the complex MRI signal, the SNR of the reconstructed wave speed map was formulated by an analytical approach assuming a high SNR, and the results were verified using Monte Carlo simulations (MCSs). It has been assumed that the noise remains approximately Gaussian when the image SNR is high enough, despite the nonlinear operations in tomoelastography inversion. Hence, the SNR threshold was determined by comparing the SNR computed by MCSs and analytical approximations. The weighting scheme was evaluated for accuracy, spatial resolution and SNR performances on simulated phantoms. MR elastography (MRE) experiments on two different phantoms were conducted. Wave speed maps were generated for simulated 3D human abdomen MRE data and experimental human abdomen MRE data. The simulation results demonstrated that the SNR-weighted inversion improved the SNR performance of the wave speed map by a factor of two compared to the performance of the original (i.e., amplitude-weighted) reconstruction. In the case of a low SNR, no bias occurred in the wave speed map when SNR weighting was used, whereas 10% bias occurred when the original weighting (i.e., amplitude weighting) was used. Thus, while not altering the accuracy or spatial resolution of the wave speed map with the proposed weighting method, the SNR of the wave speed map has been significantly improved.

Fast Robust Dejitter and Interslice Discontinuity Removal in MRI Phase Acquisitions: Application to Magnetic Resonance Elastography.

MRI phase contrast imaging methods that assemble slice-wise acquisitions into volumes can contain interslice phase discontinuities (IPDs) over the course of the scan from sources, including unavoidable physiological activity. In magnetic resonance elastography (MRE), this can alter wavelength and tissue stiffness estimates, invalidating the analysis. We first model this behavior as jitter along the z-axis of the phase of 3D complex-valued wave volumes. A two-step image processing pipeline is then proposed that removes IPDs. First, constant slicewise phase shift is removed with a novel, non-convex dejittering algorithm. Then, regional physiological noise artifacts are removed with novel filtering of 3D wavelet coefficients. Calibration of two pipeline coefficients, the dejitter parameter α and the wavelet band high-pass coefficient ωc , was first performed on a finite-element method brain phantom. A comparative investigation was then performed, on a cohort of 48 brain acquisitions, of four approaches to IPDs: 1) the proposed method; 2) a "control" condition of neglect of IPDs; 3) an anisotropic wavelet-based method; and 4) a method of in-plane (2D) processing. The present method showed medians of [Formula: see text] Pa for a multifrequency wave inversion centered at 40 Hz which was within 6% of methods 3) and 4), while neglect produced [Formula: see text] estimates a mean of 17% lower. The proposed method reduced the value range of the cohort against methods 3) and 4) by 29% and 31%, respectively. Such reduction in variance enhances the ability of brain MRE to predict subtler physiological changes. Our theoretical approach further enables more powerful applications of fundamental findings in noise and denoising to MRE.

Heterogeneous Multifrequency Direct Inversion (HMDI) for magnetic resonance elastography with application to a clinical brain exam.

A new viscoelastic wave inversion method for MRE, called Heterogeneous Multifrequency Direct Inversion (HMDI), was developed which accommodates heterogeneous elasticity within a direct inversion (DI) by incorporating first-order gradients and combining results from a narrow band of multiple frequencies. The method is compared with a Helmholtz-type DI, Multifrequency Dual Elasto-Visco inversion (MDEV), both on ground-truth Finite Element Method simulations at varied noise levels and a prospective in vivo brain cohort of 48 subjects ages 18-65. In simulated data, MDEV recovered background material within 5% and HMDI within 1% of prescribed up to SNR of 20 dB. In vivo HMDI and MDEV were then combined with segmentation from SPM to create a fully automated "brain palpation" exam for both whole brain (WB), and brain white matter (WM), measuring two parameters, the complex modulus magnitude |G*| , which measures tissue "stiffness", and the slope of |G*| values across frequencies, a measure of viscous dispersion. |G*| values for MDEV and HMDI were comparable to the literature (for a 3-frequency set centered at 50 Hz, WB means were 2.17 and 2.15 kPa respectively, and WM means were 2.47 and 2.49 kPa respectively). Both methods showed moderate correlation to age in both WB and WM, for both |G*| and |G*| slope, with Pearson's r ≥ 0.4 in the most sensitive frequency sets. In comparison to MDEV, HMDI showed better preservation of recovered target shapes, more noise-robustness, and stabler recovery values in regions with rapid property change, however summary statistics for both methods were quite similar. By eliminating homogeneity assumptions within a fast, fully automatic, regularization-free direct inversion, HMDI appears to be a worthwhile addition to the MRE image reconstruction repertoire. In addition to supporting the literature showing decrease in brain viscoelasticity with age, our work supports a wide range of inter-individual variation in brain MRE results.