Feasibility of high temporal resolution breast DCE-MRI using compressed sensing theory

Digital reference object toolkit of breast DCE MRI for quantitative evaluation of image reconstruction and analysis methods

PURPOSE

To develop a digital reference object (DRO) toolkit to generate realistic breast DCE-MRI data for quantitative assessment of image reconstruction and data analysis methods.

METHODS

A simulation framework in a form of DRO toolkit has been developed using the ultrafast and conventional breast DCE-MRI data of 53 women with malignant (n = 25) or benign (n = 28) lesions. We segmented five anatomical regions and performed pharmacokinetic analysis to determine the ranges of pharmacokinetic parameters for each segmented region. A database of the segmentations and their pharmacokinetic parameters is included in the DRO toolkit that can generate a large number of realistic breast DCE-MRI data. We provide two potential examples for our DRO toolkit: assessing the accuracy of an image reconstruction method using undersampled simulated radial k-space data and assessing the impact of theB 1 + $$ {\mathrm{B}}_1^{+} $$ field inhomogeneity on estimated parameters.

RESULTS

The estimated pharmacokinetic parameters for each region showed agreement with previously reported values. For the assessment of the reconstruction method, it was found that the temporal regularization resulted in significant underestimation of estimated parameters by up to 57% and 10% with the weighting factor λ = 0.1 and 0.01, respectively. We also demonstrated that spatial discrepancy ofv p $$ {v}_p $$ andPS $$ \mathrm{PS} $$ increase to about 33% and 51% without correction forB 1 + $$ {\mathrm{B}}_1^{+} $$ field.

CONCLUSION

We have developed a DRO toolkit that includes realistic morphology of tumor lesions along with the expected pharmacokinetic parameter ranges. This simulation framework can generate many images for quantitative assessment of DCE-MRI reconstruction and analysis methods.

Feasibility of spinal cord imaging at 7 T using rosette trajectory with magnetization transfer preparation and compressed sensing

MRI is a valuable diagnostic tool to investigate spinal cord (SC) pathology. SC MRI can benefit from the increased signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) at ultra-high fields such as 7 T. However, SC MRI acquisitions with routine Cartesian readouts are prone to image artifacts caused by physiological motion. MRI acquisition techniques with non-Cartesian readouts such as rosette can help reduce motion artifacts. The purpose of this study was to demonstrate the feasibility of high-resolution SC imaging using rosette trajectory with magnetization transfer preparation (MT-prep) and compressed sensing (CS) at 7 T. Five healthy volunteers participated in the study. Images acquired with rosette readouts demonstrated reduced motion artifacts compared to the standard Cartesian readouts. The combination of multi-echo rosette-readout images improved the CNR by approximately 50% between the gray matter (GM) and white matter (WM) compared to single-echo images. MT-prep images showed excellent contrast between the GM and WM with magnetization transfer ratio (MTR) and cerebrospinal fluid normalized MT signal (MTCSF) = 0.12 ± 0.017 and 0.74 ± 0.013, respectively, for the GM; and 0.18 ± 0.011 and 0.58 ± 0.009, respectively, for the WM. Under-sampled acquisition using rosette readout with CS reconstruction demonstrated up to 6 times faster scans with comparable image quality as the fully-sampled acquisition.

Survival rate prediction of nasopharyngeal carcinoma patients based on MRI and gene expression using a deep neural network

To achieve a better treatment regimen and follow-up assessment design for intensity-modulated radiotherapy (IMRT)-treated nasopharyngeal carcinoma (NPC) patients, an accurate progression-free survival (PFS) time prediction algorithm is needed. We propose developing a PFS prediction model of NPC patients after IMRT treatment using a deep learning method and comparing that with the traditional texture analysis method. One hundred and fifty-one NPC patients were included in this retrospective study. T1-weighted, proton density and dynamic contrast-enhanced magnetic resonance (MR) images were acquired. The expression level of five genes (HIF-1α, EGFR, PTEN, Ki-67, and VEGF) and infection of Epstein-Barr (EB) virus were tested. A residual network was trained to predict PFS from MR images. The output as well as patient characteristics were combined using a linear regression model to provide a final PFS prediction. The prediction accuracy was compared with that of the traditional texture analysis method. A regression model combining the deep learning output with HIF-1α expression and Epstein-Barr infection provides the best PFS prediction accuracy (Spearman correlation R²  = 0.53; Harrell's C-index = 0.82; receiver operative curve [ROC] analysis area under the curve [AUC] = 0.88; log-rank test hazard ratio [HR] = 8.45), higher than a regression model combining texture analysis with HIF-1α expression (Spearman correlation R²  = 0.14; Harrell's C-index =0.68; ROC analysis AUC = 0.76; log-rank test HR = 2.85). The deep learning method does not require a manually drawn tumor region of interest. MR image processing using deep learning combined with patient characteristics can provide accurate PFS prediction for nasopharyngeal carcinoma patients and does not rely on specific kernels or tumor regions of interest, which is needed for the texture analysis method.

Dynamic contrast-enhanced MRI parametric mapping using high spatiotemporal resolution Golden-angle RAdial Sparse Parallel MRI and iterative joint estimation of the arterial input function and pharmacokinetic parameters

The aim of this work is to develop a data-driven quantitative dynamic contrast-enhanced (DCE) MRI technique using Golden-angle RAdial Sparse Parallel (GRASP) MRI with high spatial resolution and high flexible temporal resolution and pharmacokinetic (PK) analysis with an arterial input function (AIF) estimated directly from the data obtained from each patient. DCE-MRI was performed on 13 patients with gynecological malignancy using a 3-T MRI scanner with a single continuous golden-angle stack-of-stars acquisition and image reconstruction with two temporal resolutions, by exploiting a unique feature in GRASP that reconstructs acquired data with user-defined temporal resolution. Joint estimation of the AIF (both AIF shape and delay) and PK parameters was performed with an iterative algorithm that alternates between AIF and PK estimation. Computer simulations were performed to determine the accuracy (expressed as percentage error [PE]) and precision of the estimated parameters. PK parameters (volume transfer constant [K^trans ], fractional volume of the extravascular extracellular space [v_e ], and blood plasma volume fraction [v_p ]) and normalized root-mean-square error [nRMSE] (%) of the fitting errors for the tumor contrast kinetic data were measured both with population-averaged and data-driven AIFs. On patient data, the Wilcoxon signed-rank test was performed to compare nRMSE. Simulations demonstrated that GRASP image reconstruction with a temporal resolution of 1 s/frame for AIF estimation and 5 s/frame for PK analysis resulted in an absolute PE of less than 5% in the estimation of K^trans and v_e , and less than 11% in the estimation of v_p . The nRMSE (mean ± SD) for the dual temporal resolution image reconstruction and data-driven AIF was 0.16 ± 0.04 compared with 0.27 ± 0.10 (p < 0.001) with 1 s/frame using population-averaged AIF, and 0.23 ± 0.07 with 5 s/frame using population-averaged AIF (p < 0.001). We conclude that DCE-MRI data acquired and reconstructed with the GRASP technique at dual temporal resolution can successfully be applied to jointly estimate the AIF and PK parameters from a single acquisition resulting in data-driven AIFs and voxelwise PK parametric maps.

Comparison of data-driven and general temporal constraints on compressed sensing for breast DCE MRI

PURPOSE

Current breast DCE-MRI strategies provide high sensitivity for cancer detection but are known to be insufficient in fully capturing rapidly changing contrast kinetics at high spatial resolution across both breasts. Advanced acquisition and reconstruction strategies aim to improve spatial and temporal resolution and increase specificity for disease characterization. In this work, we evaluate the spatial and temporal fidelity of a modified data-driven low-rank-based model (known as MOCCO, model consistency condition) compressed-sensing (CS) reconstruction compared to CS with temporal total variation with radial acquisition for high spatial-temporal breast DCE MRI.

METHODS

Reconstruction performance was characterized using numerical simulations of a golden-angle stack-of-stars breast DCE-MRI acquisition at 5-second temporal resolution. Specifically, MOCCO was compared with CS total variation and conventional SENSE reconstructions. The temporal model for MOCCO was prelearned over the source data, whereas CS total variation was performed using a first-order temporal gradient sparsity transform.

RESULTS

The MOCCO reconstruction was able to capture rapid lesion kinetics while providing high image quality across a range of optimal regularization values. It also recovered kinetics in small lesions (1.5 mm) in line-profile analysis and error images, whereas g-factor maps showed relatively low and constant values with no significant artifacts. The CS-TV method demonstrated either recovery of high spatial resolution with reduced temporal accuracy using large regularization values, or recovery of rapid lesion kinetics with reduced image quality using low regularization values.

CONCLUSION

Simulations demonstrated that MOCCO with radial acquisition provides a robust imaging technique for improving temporal fidelity, while maintaining high spatial resolution and image quality in the setting of bilateral breast DCE MRI.

Evaluating feasibility of high resolution T1-perfusion MRI with whole brain coverage using compressed SENSE: Application to glioma grading

PURPOSE

To evaluate the efficacy of optimized T1-Perfusion MRI protocol (protocol-2) with whole brain coverage and improved spatial resolution using Compressed-SENSE (CSENSE) to differentiate high-grade-glioma (HGG) and low-grade-glioma (LGG) and to compare it with the conventional protocol (protocol-1) with partial brain coverage used in our center.

METHODS

This study included MRI data from 5 healthy volunteers, a phantom and 126 brain tumor patients. Current study had two parts: To analyze the effect of CSENSE on 3D-T1-weighted (W) fast-field-echo (FFE) images, T1-W, dual-PDT2-W turbo-spin-echo images and T1 maps, and to evaluate the performance of high resolution T1-Perfusion MRI protocol with whole brain coverage optimized using CSENSE. Coefficient-of-Variation (COV), Relative-Percentage-Error (RPE), Normalized-Mean-Squared-Error (NMSE) and qualitative scoring were used for the former study. Tracer-kinetic (K^trans,v_e,v_p) and hemodynamic (rCBV,rCBF) parameters computed from both protocols were used to differentiate LGG and HGG.

RESULTS

The image quality of all structural images was found to be of diagnostic quality till R = 4. NMSE in healthy T1-W-FFE images and COV in phantom images increased with-respect-to R and images provided optimum quality till R = 4. Structural images and maps exhibited artefacts from R = 6. All parameters in tumor tissue and hemodynamic parameters in healthy gray matter tissue computed from both protocols were not significantly different. Parameters computed from protocol-2 performed better in terms of glioma grading. For both protocols, rCBF performed least (AUC = 0.759 and 0.851) and combination of all parameters performed best (AUC = 0.890 and 0.964).

CONCLUSION

CSENSE (R = 4) can be used to improve the resolution and brain coverage for T1-Perfusion analysis used to differentiate gliomas.

Time-domain principal component reconstruction (tPCR): A more efficient and stable iterative reconstruction framework for non-Cartesian functional MRI

PURPOSE

To improve the reconstruction efficiency (i.e., computational load) and stability of iterative reconstruction for non-Cartesian fMRI when using high undersampling rates and/or in the presence of strong off-resonance effects.

THEORY AND METHODS

The magnetic resonance encephalography (MREG) sequence with 3D non-Cartesian trajectory and 0.1s repetition time (TR) was applied to acquire fMRI datasets. Different from a conventional time-point-by-time-point sequential reconstruction (SR), the proposed time-domain principal component reconstruction (tPCR) performs three steps: (1) decomposing the k-t-space fMRI datasets into time-domain principal component space using singular value decomposition, (2) reconstructing each principal component with redistributed computation power according to their weights, and (3) combining the reconstructed principal components back to image-t-space. The comparison of reconstruction accuracy was performed by simulation experiments and then verified in real fMRI data.

RESULTS

The simulation experiments showed that the proposed tPCR was able to significantly reduce reconstruction errors, and subsequent functional activation errors, relative to SR at identical computational cost. Alternatively, at fixed reconstruction accuracy, computation time was greatly reduced. The improved performance was particularly obvious for L1-norm nonlinear reconstructions relative to L2-norm linear reconstructions and robust to different regularization strength, undersampling rates, and off-resonance effects intensity. By examining activation maps, tPCR was also found to give similar improvements in real fMRI experiments.

CONCLUSION

The proposed proof-of-concept tPCR framework could improve (1) the reconstruction efficiency of iterative reconstruction, and (2) the reconstruction stability especially for nonlinear reconstructions. As a practical consideration, the improved reconstruction speed promotes the application of highly undersampled non-Cartesian fast fMRI.

Impact of the Number of Iterations in Compressed Sensing Reconstruction on Ultrafast Dynamic Contrast-enhanced Breast MR Imaging

PURPOSE

To assess the impact of the number of iterations of compressed sensing (CS) reconstruction on the kinetic parameters and image quality in dynamic contrast-enhanced (DCE)-MRI of the breast, with prospectively undersampled CS-accelerated scans.

MATERIALS AND METHODS

Breast examinations including ultrafast DCE-MRI using CS were conducted for 21 patients. Images were reconstructed with different numbers of iterations. The peak enhancement ratio of the aorta and wash-in slope, initial area under the curve, and Ktrans of the breast lesions were measured. The root mean square error and structural similarity between the images using 50 iterations and images with a lower number of iterations were evaluated as criterion for quantitative image evaluation.

RESULTS

Using an insufficient number of iterations, the contrast-enhanced effect was highly underestimated. In all semi-quantitative parameters, the number of iterations that stabilized the parameters in malignant lesions was higher than that in benign lesions. At least 15 iterations were needed for semi-quantitative parameters. For Ktrans, there were no significant differences between 10 and 50 iterations in both malignant and benign lesions.

CONCLUSION

The kinetic parameters using ultrafast DCE-MRI with CS are affected by the number of iterations, especially in malignant lesions. However, if the images are reconstructed with an adequate number of iterations, ultrafast DCE-MRI with CS can be a powerful technique having high temporal and spatial resolution.

New parameters of ultrafast dynamic contrast-enhanced breast MRI using compressed sensing

BACKGROUND

Ultrafast dynamic contrast-enhanced (UF-DCE) breast MRI is considered a promising method of accelerated breast MRI. However, the value of new kinetic parameters derived from UF-DCE need clinical evaluation.

PURPOSE

To evaluate the diagnostic performance of the maximum slope (MS), time to enhancement (TTE), and time interval between arterial and venous visualization (AVI) derived from UF-DCE MRI using compressed sensing (CS).

STUDY TYPE

Retrospective.

POPULATION

Seventy-five patients with histologically proven breast lesions. The total number of analyzed lesions was 90 (61 malignant and 29 benign).

FIELD STRENGTH/SEQUENCE

3T MRI with UF-DCE MRI based on the 3D gradient-echo volumetric interpolated breath-hold examination (VIBE) sequence using incoherent k-space sampling combined with a CS reconstruction followed by conventional DCE MRI.

ASSESSMENT

The diagnostic performance of the MS, TTE, AVI, and conventional kinetic analysis was analyzed and compared with histology.

STATISTICAL TESTS

Wilcoxon rank sum test, receiver operating characteristic analysis.

RESULTS

The MS was larger and the TTE and AVI were smaller for malignant lesions compared with benign lesions: MS: 29.3%/s and 18.4%/s (P < 0.001), TTE: 7.0 and 12.0 seconds (P < 0.001), AVI: 2.7 and 4.4 frames (P = 0.006) for malignant and benign lesions. The discriminating power of the MS (area under the curve [AUC], 0.76) was slightly better than that of conventional kinetic analysis (AUC, 0.69) and comparable to that of the TTE and AVI (AUC, 0.78 and 0.76 for TTE and AVI, respectively). Invasive lobular carcinoma had smaller MS (21.8%/s) among malignant lesions (29.3%/s).

DATA CONCLUSION

The MS, TTE, and AVI can be used to evaluate breast lesions with clinical performance equivalent to that of conventional kinetic analysis. These parameters vary among histologies.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:164-174.

Magnetic Resonance Imaging of the Brain Using Compressed Sensing - Quality Assessment in Daily Clinical Routine

PURPOSE

To assess the effect of compressed sensing (CS) on image quality and acquisition speed in routine brain magnetic resonance imaging (MRI).

METHODS

During a 2-month implementation period of CS, two senior neuroradiologists, one MRI physicist and one application specialist optimized the CS acceleration factor to reduce scan time and improve spatial resolution, while maintaining image quality. Afterwards, two neuroradiologists independently scored image quality on a 5-point Likert scale in 3‑dimensional (3D) fluid attenuation inversion recovery (FLAIR), 3D double inversion recovery (DIR), 3D T2, 3D T1, 3D T1 + gadoteric acid, axial T2, axial FLAIR, axial T2*, and 3D arterial time-of-flight MR angiography (art. TOF) sequences acquired during 1 week before (CS-) and after (CS+) the implementation of CS. Time of acquisition was recorded for all sequences.

RESULTS

A total of 51 CS- and 48 CS+ patients were included. The median scan time reduction was 29.3% (range 0.0-58.4%), median voxel size reduction was 10.5% (0.0-33.3%). The CS+ image quality was rated superior for 3D FLAIR (p < 0.001), 3D T2 (p = 0.001), and axial T2* sequences (p = 0.024). For all other sequences, no statistical difference in image quality was observed. Interreader agreement regarding image quality was good for all sequences (weighted Cohen's κ > 0.5).

CONCLUSION

The use of CS saves considerable imaging time while allowing to increase spatial resolution in routine clinical brain MRI without loss in image quality.

Fast quantitative susceptibility reconstruction via total field inversion with improved weighted L0 norm approximation

Quantitative susceptibility mapping (QSM) is a meaningful MRI technique owing to its unique relation to actual physical tissue magnetic properties. The reconstruction of QSM is usually decomposed into three sub-problems, which are solved independently. However, this decomposition does not conform to the causes of the problems, and may cause discontinuity of parameters and error accumulation. In this paper, a fast reconstruction method named fast TFI based on total field inversion was proposed. It can accelerate the total field inversion by using a specially selected preconditioner and advanced solution of the weighted L₀ regularization. Due to the employment of an effective model, the proposed method can efficiently reconstruct the QSM of brains with lesions, where other methods may encounter problems. Experimental results from simulation and in vivo data verified that the new method has better reconstruction accuracy, faster convergence ability and excellent robustness, which may promote clinical application of QSM.

Abbreviated MRI Protocols in Breast Cancer Diagnostics

Oncologic imaging focused on the detection of breast cancer is of increasing importance, with over 1.7 million new cases detected each year worldwide. MRI of the breast has been described to be one of the most sensitive imaging modalities in breast cancer detection; however, clinical use is limited due to high costs. In the past, the objective and clinical routine of oncologic imaging was to provide one extended imaging protocol covering all potential needs and clinical implications regardless of the specific clinical indication or question. Future protocols might be more focused according to a "keep it short and simple" approach, with a reduction of patient magnet time and a limited number of images to review. Rather than replacing conventional full-diagnostic breast MRI protocols, these approaches aim at introducing a new thinking in oncologic imaging using a diversification of available imaging approaches targeted to the dedicated clinical needs of the individual patient. Here we review current approaches on using abbreviated protocols that aim to increase the clinical availability and use of breast MRI for improved early detection of breast cancer. Level of Evidence: 2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019;49:647-658.

Undersampling patterns in k-space for compressed sensing MRI using two-dimensional Cartesian sampling

In compressed sensing magnetic resonance imaging (CS-MRI), undersampling of k-space is performed to achieve faster imaging. For this process, it is important to acquire data randomly, and an optimal random undersampling pattern is required. However, random undersampling is difficult in two-dimensional (2D) Cartesian sampling. In this study, the effect of random undersampling patterns on image reconstruction was clarified using phantom and in vivo MRI, and a sampling pattern relevant for 2D Cartesian sampling in CS-MRI is suggested. The precision of image restoration was estimated with various acceleration factors and extents for the fully sampled central region of k-space. The root-mean-square error, structural similarity index, and modulation transfer function were measured, and visual assessments were also performed. The undersampling pattern was shown to influence the precision of image restoration, and an optimal undersampling pattern should be used to improve image quality; therefore, we suggest that the ideal undersampling pattern in CS-MRI for 2D Cartesian sampling is one with a high extent for the fully sampled central region of k-space.

Dynamic field-of-view imaging to increase temporal resolution in the early phase of contrast media uptake in breast DCE-MRI: A feasibility study

PURPOSE

To increase diagnostic accuracy of breast MRI by increasing temporal resolution and more accurately sampling the early kinetics of contrast media uptake. We tested the feasibility of accelerating bilateral breast DCE-MRI by reducing the FOV, allowing aliasing, and unfolding the resulting images.

METHODS

Previous experience with an "ultrafast" protocol for bilateral breast DCE-MRI (6-10 s temporal resolution) showed that the number of significantly enhancing voxels is very low in the first 30-45 s after contrast media injection. This suggests that overlap of enhancing voxels in aliased images will be very infrequent. Therefore, aliased images can be acquired during the first 30-45 s after contrast media injection and unfolded to produce full-FOV images with few errors. In a proof-of-principle test, aliased images were simulated from the first 30 s of full-FOV acquisitions. Cases with relatively dense early enhancement were selected to test this method in a worst-case scenario. In an initial test, an FOV of 60% the size of the full FOV was simulated. To reduce the probability of errors due to overlapping voxels in aliased images, we then tested a dynamic FOV approach. The FOV was progressively increased so that enhancing voxels could not overlap at multiple time-points, and areas where enhancing voxels overlapped at a given time-point could be unfolded by interpolating between the preceding and subsequent time-points (acquired with different FOVs). The simulated FOV sizes for each of the time-points were 31%, 44%, and 77% of the full FOV. Subtraction images (post- minus precontrast) were generated for aliased images and filtered to select significantly enhancing voxels. Comparison of early, highly aliased images, with later, less aliased images then helped to identify the true locations of enhancing voxels.

RESULTS

In the initial aliasing simulations, an average of 2.9% of the enhancing voxels above the chest wall overlapped in the aliased images (range 0.1%-6.7%). The similarity between simulated unfolded images and the correct full-FOV images, evaluated using CW-SSIM (complex wavelet similarity index), was 0.50 ± 0.26, 0.76 ± 0.09, and 0.80 ± 0.10 for the first, second, and third time-point, respectively (numbers closer to 1 indicate more similar images). For the dynamic FOV tests, an average of 11% of the enhancing voxels above the chest wall overlapped (range 0%-40%) due to greater aliasing at early time-points. Despite more voxels overlapping, the CW-SSIM values for the data acquired with dynamic FOVs were 0.64 ± 0.25, 0.93 ± 0.04, and 0.97 ± 0.02 for the first, second, and third time-points, respectively.

CONCLUSIONS

Dynamic FOV imaging allows accelerated bilateral breast DCE-MRI during the early contrast media uptake phase. This method relies on the sparsity of enhancement at the early phases of DCE-MRI of the breast. The results of simulations suggest that dynamic FOV imaging and unfolding produces images that are very close to fully sampled images, and allows temporal resolution as high as 2 s per image.

Feasibility of high spatiotemporal resolution for an abbreviated 3D radial breast MRI protocol

PURPOSE

To develop a volumetric imaging technique with 0.8-mm isotropic resolution and 10-s/volume rate to detect and analyze breast lesions in a bilateral, dynamic, contrast-enhanced MRI exam.

METHODS

A local low-rank temporal reconstruction approach that also uses parallel imaging and spatial compressed sensing was designed to create rapid volumetric frame rates during a contrast-enhanced breast exam (vastly undersampled isotropic projection [VIPR] spatial compressed sensing with temporal local low-rank [STELLR]). The dynamic-enhanced data are subtracted in k-space from static mask data to increase sparsity for the local low-rank approach to maximize temporal resolution. A T₁ -weighted 3D radial trajectory (VIPR iterative decomposition with echo asymmetry and least squares estimation [IDEAL]) was modified to meet the data acquisition requirements of the STELLR approach. Additionally, the unsubtracted enhanced data are reconstructed using compressed sensing and IDEAL to provide high-resolution fat/water separation. The feasibility of the approach and the dual reconstruction methodology is demonstrated using a 16-channel breast coil and a 3T MR scanner in 6 patients.

RESULTS

The STELLR temporal performance of subtracted data matched the expected temporal perfusion enhancement pattern in small and large vascular structures. Differential enhancement within heterogeneous lesions is demonstrated with corroboration from a basic reconstruction using a strict 10-second temporal footprint. Rapid acquisition, reliable fat suppression, and high spatiotemporal resolution are presented, despite significant data undersampling.

CONCLUSION

The STELLR reconstruction approach of 3D radial sampling with mask subtraction provides a high-performance imaging technique for characterizing enhancing structures within the breast. It is capable of maintaining temporal fidelity, while visualizing breast lesions with high detail over a large FOV to include both breasts.

High-resolution whole-brain DCE-MRI using constrained reconstruction: Prospective clinical evaluation in brain tumor patients

PURPOSE

To clinically evaluate a highly accelerated T1-weighted dynamic contrast-enhanced (DCE) MRI technique that provides high spatial resolution and whole-brain coverage via undersampling and constrained reconstruction with multiple sparsity constraints.

METHODS

Conventional (rate-2 SENSE) and experimental DCE-MRI (rate-30) scans were performed 20 minutes apart in 15 brain tumor patients. The conventional clinical DCE-MRI had voxel dimensions 0.9 × 1.3 × 7.0 mm(3), FOV 22 × 22 × 4.2 cm(3), and the experimental DCE-MRI had voxel dimensions 0.9 × 0.9 × 1.9 mm(3), and broader coverage 22 × 22 × 19 cm(3). Temporal resolution was 5 s for both protocols. Time-resolved images and blood-brain barrier permeability maps were qualitatively evaluated by two radiologists.

RESULTS

The experimental DCE-MRI scans showed no loss of qualitative information in any of the cases, while achieving substantially higher spatial resolution and whole-brain spatial coverage. Average qualitative scores (from 0 to 3) were 2.1 for the experimental scans and 1.1 for the conventional clinical scans.

CONCLUSIONS

The proposed DCE-MRI approach provides clinically superior image quality with higher spatial resolution and coverage than currently available approaches. These advantages may allow comprehensive permeability mapping in the brain, which is especially valuable in the setting of large lesions or multiple lesions spread throughout the brain.

Rapid reconstruction of quantitative susceptibility mapping via improved ℓ0 norm approximation

Quantitative susceptibility mapping (QSM) reconstruction is a well-known ill-posed problem. Various regularization techniques have been proposed for solving this problem. In this paper, a rapid method is proposed that uses ℓ₀ norm minimization in a gradient domain. Because ℓ₀ minimization is an NP-hard problem, a special alternating optimization strategy is employed to simplify the reconstruction algorithm. The proposed algorithm uses only simple point-wise multiplications and thresholding operations, and significantly speeds up the calculation. Both numerical simulations and in vivo experiments demonstrate that the proposed method can reconstruct susceptibility fast and accurately. Because morphology information weighted methods have achieved considerable success in QSM, we performed a quantitative comparison with some typical weighted methods, such as MEDI (morphology enabled dipole inversion), iLSQR (improved least squares algorithm), and wℓ₁ (weighted ℓ₁ norm minimization). The reconstructed results show that the proposed method can provide accurate results with a satisfactory speed.

Comparison of conventional DCE-MRI and a novel golden-angle radial multicoil compressed sensing method for the evaluation of breast lesion conspicuity

PURPOSE

To compare a novel multicoil compressed sensing technique with flexible temporal resolution, golden-angle radial sparse parallel (GRASP), to conventional fat-suppressed spoiled three-dimensional (3D) gradient-echo (volumetric interpolated breath-hold examination, VIBE) MRI in evaluating the conspicuity of benign and malignant breast lesions.

MATERIALS AND METHODS

Between March and August 2015, 121 women (24-84 years; mean, 49.7 years) with 180 biopsy-proven benign and malignant lesions were imaged consecutively at 3.0 Tesla in a dynamic contrast-enhanced (DCE) MRI exam using sagittal T1-weighted fat-suppressed 3D VIBE in this Health Insurance Portability and Accountability Act-compliant, retrospective study. Subjects underwent MRI-guided breast biopsy (mean, 13 days [1-95 days]) using GRASP DCE-MRI, a fat-suppressed radial "stack-of-stars" 3D FLASH sequence with golden-angle ordering. Three readers independently evaluated breast lesions on both sequences. Statistical analysis included mixed models with generalized estimating equations, kappa-weighted coefficients and Fisher's exact test.

RESULTS

All lesions demonstrated good conspicuity on VIBE and GRASP sequences (4.28 ± 0.81 versus 3.65 ± 1.22), with no significant difference in lesion detection (P = 0.248). VIBE had slightly higher lesion conspicuity than GRASP for all lesions, with VIBE 12.6% (0.63/5.0) more conspicuous (P < 0.001). Masses and nonmass enhancement (NME) were more conspicuous on VIBE (P < 0.001), with a larger difference for NME (14.2% versus 9.4% more conspicuous). Malignant lesions were more conspicuous than benign lesions (P < 0.001) on both sequences.

CONCLUSION

GRASP DCE-MRI, a multicoil compressed sensing technique with high spatial resolution and flexible temporal resolution, has near-comparable performance to conventional VIBE imaging for breast lesion evaluation.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2017;45:1746-1752.

Kinetic Analysis of Benign and Malignant Breast Lesions With Ultrafast Dynamic Contrast-Enhanced MRI: Comparison With Standard Kinetic Assessment

OBJECTIVE

The purposes of this study were to evaluate diagnostic parameters measured with ultrafast MRI acquisition and with standard acquisition and to compare diagnostic utility for differentiating benign from malignant lesions.

MATERIALS AND METHODS

Ultrafast acquisition is a high-temporal-resolution (7 seconds) imaging technique for obtaining 3D whole-breast images. The dynamic contrast-enhanced 3-T MRI protocol consists of an unenhanced standard and an ultrafast acquisition that includes eight contrast-enhanced ultrafast images and four standard images. Retrospective assessment was performed for 60 patients with 33 malignant and 29 benign lesions. A computer-aided detection system was used to obtain initial enhancement rate and signal enhancement ratio (SER) by means of identification of a voxel showing the highest signal intensity in the first phase of standard imaging. From the same voxel, the enhancement rate at each time point of the ultrafast acquisition and the AUC of the kinetic curve from zero to each time point of ultrafast imaging were obtained.

RESULTS

There was a statistically significant difference between benign and malignant lesions in enhancement rate and kinetic AUC for ultrafast imaging and also in initial enhancement rate and SER for standard imaging. ROC analysis showed no significant differences between enhancement rate in ultrafast imaging and SER or initial enhancement rate in standard imaging.

CONCLUSION

Ultrafast imaging is useful for discriminating benign from malignant lesions. The differential utility of ultrafast imaging is comparable to that of standard kinetic assessment in a shorter study time.

Assessment of Treatment Response With Diffusion-Weighted MRI and Dynamic Contrast-Enhanced MRI in Patients With Early-Stage Breast Cancer Treated With Single-Dose Preoperative Radiotherapy: Initial Results

Single-dose preoperative stereotactic body radiotherapy is a novel radiotherapy technique for the early-stage breast cancer, and the treatment response pattern of this technique needs to be investigated on a quantitative basis. In this work, dynamic contrast-enhanced magnetic resonance imaging and diffusion-weighted magnetic resonance imaging were used to study the treatment response pattern in a unique cohort of patients with early-stage breast cancer treated with preoperative radiation. Fifteen female qualified patients received single-dose preoperative radiotherapy with 1 of the 3 prescription doses: 15 Gy, 18 Gy, and 21 Gy. Magnetic resonance imaging scans including both diffusion-weighted magnetic resonance imaging and dynamic contrast-enhanced magnetic resonance imaging were acquired before radiotherapy for planning and after radiotherapy but before surgical resection. In diffusion-weighted magnetic resonance imaging, the regional averaged apparent diffusion coefficient was calculated. In dynamic contrast-enhanced magnetic resonance imaging, quantitative parameters K (trans) and v e were evaluated using the standard Tofts model based on the average contrast agent concentration within the region of interest, and the semiquantitative initial area under the concentration curve (iAUC6min) was also recorded. These parameters' relative changes after radiotherapy were calculated for gross tumor volume, clinical target volume, and planning target volume. The initial results showed that after radiotherapy, initial area under the concentration curve significantly increased in planning target volume (P < .006) and clinical target volume (P < .006), and v e significantly increased in planning target volume (P < .05) and clinical target volume (P < .05). Statistical studies suggested that linear correlations between treatment dose and the observed parameter changes exist in most examined tests, and among these tests, the change in gross tumor volume regional averaged apparent diffusion coefficient (P < .012) and between treatment dose and planning target volume K (trans) (P < .029) were found to be statistically significant. Although it is still preliminary, this pilot study may be useful to provide insights for future works.

Ultrafast Bilateral DCE-MRI of the Breast with Conventional Fourier Sampling: Preliminary Evaluation of Semi-quantitative Analysis

RATIONALE AND OBJECTIVES

The study aimed to evaluate the feasibility and advantages of a combined high temporal and high spatial resolution protocol for dynamic contrast-enhanced magnetic resonance imaging of the breast.

MATERIALS AND METHODS

Twenty-three patients with enhancing lesions were imaged at 3T. The acquisition protocol consisted of a series of bilateral, fat-suppressed "ultrafast" acquisitions, with 6.9- to 9.9-second temporal resolution for the first minute following contrast injection, followed by four high spatial resolution acquisitions with 60- to 79.5-second temporal resolution. All images were acquired with standard uniform Fourier sampling. A filtering method was developed to reduce noise and detect significant enhancement in the high temporal resolution images. Time of arrival (TOA) was defined as the time at which each voxel first satisfied all the filter conditions, relative to the time of initial arterial enhancement.

RESULTS

Ultrafast images improved visualization of the vasculature feeding and draining lesions. A small percentage of the entire field of view (<6%) enhanced significantly in the 30 seconds following contrast injection. Lesion conspicuity was highest in early ultrafast images, especially in cases with marked parenchymal enhancement. Although the sample size was relatively small, the average TOA for malignant lesions was significantly shorter than the TOA for benign lesions. Significant differences were also measured in other parameters descriptive of early contrast media uptake kinetics (P < 0.05).

CONCLUSIONS

Ultrafast imaging in the first minute of dynamic contrast-enhanced magnetic resonance imaging of the breast has the potential to add valuable information on early contrast dynamics. Ultrafast imaging could allow radiologists to confidently identify lesions in the presence of marked background parenchymal enhancement.

Accelerated Brain DCE-MRI Using Iterative Reconstruction With Total Generalized Variation Penalty for Quantitative Pharmacokinetic Analysis: A Feasibility Study

PURPOSE

To investigate the feasibility of using undersampled k-space data and an iterative image reconstruction method with total generalized variation penalty in the quantitative pharmacokinetic analysis for clinical brain dynamic contrast-enhanced magnetic resonance imaging.

METHODS

Eight brain dynamic contrast-enhanced magnetic resonance imaging scans were retrospectively studied. Two k-space sparse sampling strategies were designed to achieve a simulated image acquisition acceleration factor of 4. They are (1) a golden ratio-optimized 32-ray radial sampling profile and (2) a Cartesian-based random sampling profile with spatiotemporal-regularized sampling density constraints. The undersampled data were reconstructed to yield images using the investigated reconstruction technique. In quantitative pharmacokinetic analysis on a voxel-by-voxel basis, the rate constant K^trans in the extended Tofts model and blood flow F_B and blood volume V_B from the 2-compartment exchange model were analyzed. Finally, the quantitative pharmacokinetic parameters calculated from the undersampled data were compared with the corresponding calculated values from the fully sampled data. To quantify each parameter's accuracy calculated using the undersampled data, error in volume mean, total relative error, and cross-correlation were calculated.

RESULTS

The pharmacokinetic parameter maps generated from the undersampled data appeared comparable to the ones generated from the original full sampling data. Within the region of interest, most derived error in volume mean values in the region of interest was about 5% or lower, and the average error in volume mean of all parameter maps generated through either sampling strategy was about 3.54%. The average total relative error value of all parameter maps in region of interest was about 0.115, and the average cross-correlation of all parameter maps in region of interest was about 0.962. All investigated pharmacokinetic parameters had no significant differences between the result from original data and the reduced sampling data.

CONCLUSION

With sparsely sampled k-space data in simulation of accelerated acquisition by a factor of 4, the investigated dynamic contrast-enhanced magnetic resonance imaging pharmacokinetic parameters can accurately estimate the total generalized variation-based iterative image reconstruction method for reliable clinical application.

3D Cartesian MRI with compressed sensing and variable view sharing using complementary poisson-disc sampling

PURPOSE

To enable robust, high spatio-temporal-resolution three-dimensional Cartesian MRI using a scheme incorporating a novel variable density random k-space sampling trajectory allowing flexible and retrospective selection of the temporal footprint with compressed sensing (CS).

METHODS

A complementary Poisson-disc k-space sampling trajectory was designed to allow view sharing and varying combinations of reduced view sharing with CS from the same prospective acquisition. These schemes were used for two-point Dixon-based dynamic contrast-enhanced MRI (DCE-MRI) of the breast and abdomen. Results were validated in vivo with a novel approach using variable-flip-angle data, which was retrospectively accelerated using the same methods but offered a ground truth.

RESULTS

In breast DCE-MRI, the temporal footprint could be reduced 2.3-fold retrospectively without introducing noticeable artifacts, improving depiction of rapidly enhancing lesions. Further, experiments with variable-flip-angle data showed that reducing view sharing improved accuracy in reconstruction and T₁ mapping. In abdominal MRI, 2.3-fold and 3.6-fold reductions in temporal footprint allowed reduced motion artifacts.

CONCLUSION

The complementary-Poisson-disc k-space sampling trajectory allowed a retrospective spatiotemporal resolution tradeoff using CS and view sharing, imparting robustness to motion and contrast enhancement. The technique was also validated using a novel approach of fully acquired variable-flip-angle acquisition. Magn Reson Med 77:1774-1785, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Pattern identification of biomedical images with time series: Contrasting THz pulse imaging with DCE-MRIs

OBJECTIVE

We provide a survey of recent advances in biomedical image analysis and classification from emergent imaging modalities such as terahertz (THz) pulse imaging (TPI) and dynamic contrast-enhanced magnetic resonance images (DCE-MRIs) and identification of their underlining commonalities.

METHODS

Both time and frequency domain signal pre-processing techniques are considered: noise removal, spectral analysis, principal component analysis (PCA) and wavelet transforms. Feature extraction and classification methods based on feature vectors using the above processing techniques are reviewed. A tensorial signal processing de-noising framework suitable for spatiotemporal association between features in MRI is also discussed.

VALIDATION

Examples where the proposed methodologies have been successful in classifying TPIs and DCE-MRIs are discussed.

RESULTS

Identifying commonalities in the structure of such heterogeneous datasets potentially leads to a unified multi-channel signal processing framework for biomedical image analysis.

CONCLUSION

The proposed complex valued classification methodology enables fusion of entire datasets from a sequence of spatial images taken at different time stamps; this is of interest from the viewpoint of inferring disease proliferation. The approach is also of interest for other emergent multi-channel biomedical imaging modalities and of relevance across the biomedical signal processing community.

Compressed sensing MRI: a review of the clinical literature

MRI is one of the most dynamic and safe imaging techniques available in the clinic today. However, MRI acquisitions tend to be slow, limiting patient throughput and limiting potential indications for use while driving up costs. Compressed sensing (CS) is a method for accelerating MRI acquisition by acquiring less data through undersampling of k-space. This has the potential to mitigate the time-intensiveness of MRI. The limited body of research evaluating the effects of CS on MR images has been mostly positive with regards to its potential as a clinical tool. Studies have successfully accelerated MRI with this technology, with varying degrees of success. However, more must be performed before its diagnostic efficacy and benefits are clear. Studies involving a greater number radiologists and images must be completed, rating CS based on its diagnostic efficacy. Also, standardized methods for determining optimal imaging parameters must be developed.

High temporal resolution dynamic contrast-enhanced MRI using compressed sensing-combined sequence in quantitative renal perfusion measurement

PURPOSE

To demonstrate the feasibility of the improved temporal resolution by using compressed sensing (CS) combined imaging sequence in dynamic contrast-enhanced MRI (DCE-MRI) of kidney, and investigate its quantitative effects on renal perfusion measurements.

MATERIALS AND METHODS

Ten rabbits were included in the accelerated scans with a CS-combined 3D pulse sequence. To evaluate the image quality, the signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were compared between the proposed CS strategy and the conventional full sampling method. Moreover, renal perfusion was estimated by using the separable compartmental model in both CS simulation and realistic CS acquisitions.

RESULTS

The CS method showed DCE-MRI images with improved temporal resolution and acceptable image contrast, while presenting significantly higher SNR than the fully sampled images (p<.01) at 2-, 3- and 4-X acceleration. In quantitative measurements, renal perfusion results were in good agreement with the fully sampled one (concordance correlation coefficient=0.95, 0.91, 0.88) at 2-, 3- and 4-X acceleration in CS simulation. Moreover, in realistic acquisitions, the estimated perfusion by the separable compartmental model exhibited no significant differences (p>.05) between each CS-accelerated acquisition and the full sampling method.

CONCLUSION

The CS-combined 3D sequence could improve the temporal resolution for DCE-MRI in kidney while yielding diagnostically acceptable image quality, and it could provide effective measurements of renal perfusion.

An efficient calculation method for pharmacokinetic parameters in brain permeability study using dynamic contrast-enhanced MRI

PURPOSE

To develop an efficient method for calculating pharmacokinetic (PK) parameters in brain DCE-MRI permeability studies.

METHODS

A linear least-squares fitting algorithm based on a derivative expression of the two-compartment PK model was proposed to analytically solve for the PK parameters. Noise in the expression was minimized through low-pass filtering. Simulation studies were conducted in which the proposed method was compared with two existing methods in terms of accuracy and efficiency. Five in vivo brain studies were demonstrated for potential clinical application.

RESULTS

In the simulation studies using chosen parameter values, the calculated percent difference of K(trans) by the proposed method was <5.0% with a temporal resolution (Δt) < 5 s, and the accuracies of all parameter results were better or comparable to existing methods. When analyzed within certain parameter intensity ranges, the proposed method was more accurate than the existing methods and improved the efficiency by a factor of up to 458 for a Δt = 1 s and up to 38 for a Δt = 5 s. In the in vivo study, the calculated parameters using the proposed method were comparable to those using the existing methods with improved efficiencies.

CONCLUSIONS

An efficient method was developed for the accurate and efficient calculation of parameters in brain DCE-MRI permeability studies.

Review of treatment assessment using DCE-MRI in breast cancer radiation therapy

As a noninvasive functional imaging technique, dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is being used in oncology to measure properties of tumor microvascular structure and permeability. Studies have shown that parameters derived from certain pharmacokinetic models can be used as imaging biomarkers for tumor treatment response. The use of DCE-MRI for quantitative and objective assessment of radiation therapy has been explored in a variety of methods and tumor types. However, due to the complexity in imaging technology and divergent outcomes from different pharmacokinetic approaches, the method of using DCE-MRI in treatment assessment has yet to be standardized, especially for breast cancer. This article reviews the basic principles of breast DCE-MRI and recent studies using DCE-MRI in treatment assessment. Technical and clinical considerations are emphasized with specific attention to assessment of radiation treatment response.

Potential of compressed sensing in quantitative MR imaging of cancer

Classic signal processing theory dictates that, in order to faithfully reconstruct a band-limited signal (e.g., an image), the sampling rate must be at least twice the maximum frequency contained within the signal, i.e., the Nyquist frequency. Recent developments in applied mathematics, however, have shown that it is often possible to reconstruct signals sampled below the Nyquist rate. This new method of compressed sensing (CS) requires that the signal have a concise and extremely dense representation in some mathematical basis. Magnetic resonance imaging (MRI) is particularly well suited for CS approaches, owing to the flexibility of data collection in the spatial frequency (Fourier) domain available in most MRI protocols. With custom CS acquisition and reconstruction strategies, one can quickly obtain a small subset of the full data and then iteratively reconstruct images that are consistent with the acquired data and sparse by some measure. Successful use of CS results in a substantial decrease in the time required to collect an individual image. This extra time can then be harnessed to increase spatial resolution, temporal resolution, signal-to-noise, or any combination of the three. In this article, we first review the salient features of CS theory and then discuss the specific barriers confronting CS before it can be readily incorporated into clinical quantitative MRI studies of cancer. We finally illustrate applications of the technique by describing examples of CS in dynamic contrast-enhanced MRI and dynamic susceptibility contrast MRI.

Compressed sensing reconstruction for undersampled breath-hold radial cine imaging with auxiliary free-breathing data

PURPOSE

To improve compressed sensing (CS) reconstruction of accelerated breath-hold (BH) radial cine magnetic resonance imaging (MRI) by exploiting auxiliary data acquired between different BHs.

MATERIALS AND METHODS

Cardiac function is usually assessed using segmented cine acquisitions over multiple BHs to cover the entire left ventricle (LV). Subjects are given a resting period between adjacent BHs, when conventionally no data are acquired and subjects rest in the scanner. In this study the resting periods between BHs were used to acquire additional free-breathing (FB) data, which are subsequently used to generate a sparsity constraint for each cardiac phase. Images reconstructed using the proposed sparsity constraint were compared with conventional CS using a composite image generated by averaging different cardiac phases. The efficacy of the proposed reconstruction was compared using indices of LV function and blood-myocardium sharpness.

RESULTS

The proposed method provided accurate LV ejection fraction measurements for 33% and 20% sampled datasets compared with fully sampled reference images, and showed 14% and 11% higher blood-myocardium border sharpness scores compared to the conventional CS.

CONCLUSION

The FB data acquired during resting periods can be efficiently used to improve the image quality of the undersampled BH data without increasing the total scan time.

Differentiation of breast cancer from fibroadenoma with dual-echo dynamic contrast-enhanced MRI

Dynamic contrast-enhanced magnetic resonance imaging (DCE MRI) of the breast is a routinely used imaging method which is highly sensitive for detecting breast malignancy. Specificity, though, remains suboptimal. Dynamic susceptibility contrast magnetic resonance imaging (DSC MRI), an alternative dynamic contrast imaging technique, evaluates perfusion-related parameters unique from DCE MRI. Previous work has shown that the combination of DSC MRI with DCE MRI can improve diagnostic specificity, though an additional administration of intravenous contrast is required. Dual-echo MRI can measure both T₁W DCE MRI and T₂*W DSC MRI parameters with a single contrast bolus, but has not been previously implemented in breast imaging. We have developed a dual-echo gradient-echo sequence to perform such simultaneous measurements in the breast, and use it to calculate the semi-quantitative T₁W and T₂*W related parameters such as peak enhancement ratio, time of maximal enhancement, regional blood flow, and regional blood volume in 20 malignant lesions and 10 benign fibroadenomas in 38 patients. Imaging parameters were compared to surgical or biopsy obtained tissue samples. Receiver operating characteristic (ROC) curves and area under the ROC curves were calculated for each parameter and combination of parameters. The time of maximal enhancement derived from DCE MRI had a 90% sensitivity and 69% specificity for predicting malignancy. When combined with DSC MRI derived regional blood flow and volume parameters, sensitivity remained unchanged at 90% but specificity increased to 80%. In conclusion, we show that dual-echo MRI with a single administration of contrast agent can simultaneously measure both T₁W and T₂*W related perfusion and kinetic parameters in the breast and the combination of DCE MRI and DSC MRI parameters improves the diagnostic performance of breast MRI to differentiate breast cancer from benign fibroadenomas.

Under-sampling trajectory design for compressed sensing based DCE-MRI

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) needs high temporal and spatial resolution to accurately estimate quantitative parameters and characterize tumor vasculature. Compressed Sensing (CS) has the potential to accomplish this mutual importance. However, the randomness in CS under-sampling trajectory designed using the traditional variable density (VD) scheme may translate to uncertainty in kinetic parameter estimation when high reduction factors are used. Therefore, accurate parameter estimation using VD scheme usually needs multiple adjustments on parameters of Probability Density Function (PDF), and multiple reconstructions even with fixed PDF, which is inapplicable for DCE-MRI. In this paper, an under-sampling trajectory design which is robust to the change on PDF parameters and randomness with fixed PDF is studied. The strategy is to adaptively segment k-space into low-and high frequency domain, and only apply VD scheme in high-frequency domain. Simulation results demonstrate high accuracy and robustness comparing to VD design.

Location constrained approximate message passing for compressed sensing MRI

Iterative thresholding methods have been extensively studied as faster alternatives to convex optimization methods for solving large-sized problems in compressed sensing. A novel iterative thresholding method called LCAMP (Location Constrained Approximate Message Passing) is presented for reducing computational complexity and improving reconstruction accuracy when a nonzero location (or sparse support) constraint can be obtained from view shared images. LCAMP modifies the existing approximate message passing algorithm by replacing the thresholding stage with a location constraint, which avoids adjusting regularization parameters or thresholding levels. This work is first compared with other conventional reconstruction methods using random one-dimention signals and then applied to dynamic contrast-enhanced breast magnetic resonance imaging to demonstrate the excellent reconstruction accuracy (less than 2% absolute difference) and low computation time (5-10 s using Matlab) with highly undersampled three-dimentional data (244 × 128 × 48; overall reduction factor = 10).

A feasible high spatiotemporal resolution breast DCE-MRI protocol for clinical settings

Three dimensional bilateral imaging is the standard for most clinical breast dynamic contrast-enhanced (DCE) MRI protocols. Because of high spatial resolution (sRes) requirement, the typical 1-2 min temporal resolution (tRes) afforded by a conventional full-k-space-sampling gradient echo (GRE) sequence precludes meaningful and accurate pharmacokinetic analysis of DCE time-course data. The commercially available, GRE-based, k-space undersampling and data sharing TWIST (time-resolved angiography with stochastic trajectories) sequence was used in this study to perform DCE-MRI exams on thirty one patients (with 36 suspicious breast lesions) before their biopsies. The TWIST DCE-MRI was immediately followed by a single-frame conventional GRE acquisition. Blinded from each other, three radiologist readers assessed agreements in multiple lesion morphology categories between the last set of TWIST DCE images and the conventional GRE images. Fleiss' κ test was used to evaluate inter-reader agreement. The TWIST DCE time-course data were subjected to quantitative pharmacokinetic analyses. With a four-channel phased-array breast coil, the TWIST sequence produced DCE images with 20 s or less tRes and ~ 1.0×1.0×1.4 mm(3) sRes. There were no significant differences in signal-to-noise (P=.45) and contrast-to-noise (P=.51) ratios between the TWIST and conventional GRE images. The agreements in morphology evaluations between the two image sets were excellent with the intra-reader agreement ranging from 79% for mass margin to 100% for mammographic density and the inter-reader κ value ranging from 0.54 (P<.0001) for lesion size to 1.00 (P<.0001) for background parenchymal enhancement. Quantitative analyses of the DCE time-course data provided higher breast cancer diagnostic accuracy (91% specificity at 100% sensitivity) than the current clinical practice of morphology and qualitative kinetics assessments. The TWIST sequence may be used in clinical settings to acquire high spatiotemporal resolution breast DCE-MRI images for both precise lesion morphology characterization and accurate pharmacokinetic analysis.

Temporal/spatial resolution improvement of in vivo DCE-MRI with compressed sensing-optimized FLASH

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides critical information regarding tumor perfusion and permeability by injecting a T(1) contrast agent, such as Gd-DTPA, and making a time-resolved measurement of signal increase. Both temporal and spatial resolutions are required to be high to achieve an accurate and reproducible estimation of tumor perfusion. However, the dynamic nature of the DCE experiment limits simultaneous improvement of temporal and spatial resolution by conventional methods. Compressed sensing (CS) has become an important tool for the acceleration of imaging times in MRI, which is achieved by enabling the reconstruction of subsampled data. Similarly, CS algorithms can be utilized to improve the temporal/spatial resolution of DCE-MRI, and several works describing retrospective simulations have demonstrated the feasibility of such improvements. In this study, the fast low angle shot sequence was modified to implement a Cartesian, CS-optimized, sub-Nyquist phase encoding acquisition/reconstruction with multiple two-dimensional slice selections and was tested on water phantoms and animal tumor models. The mean voxel-level concordance correlation coefficient for Ak(ep) values obtained from ×4 and ×8 accelerated and the fully sampled data was 0.87±0.11 and 0.83±0.11, respectively (n=6), with optimized CS parameters. In this case, the reduction of phase encoding steps made possible by CS reconstruction improved effectively the temporal/spatial resolution of DCE-MRI data using an in vivo animal tumor model (n=6) and may be useful for the investigation of accelerated acquisitions in preclinical and clinical DCE-MRI trials.

Fast 3D contrast enhanced MRI of the liver using temporal resolution acceleration with constrained evolution reconstruction

Time-resolved imaging is crucial for the accurate diagnosis of liver lesions. Current contrast enhanced liver magnetic resonance imaging acquires a few phases in sequential breath-holds. The image quality is susceptible to bolus timing errors, which could result in missing the critical arterial phase. This impairs the detection of malignant tumors that are supplied primarily by the hepatic artery. In addition, the temporal resolution may be too low to reliably separate the arterial phase from the portal venous phase. In this study, a method called temporal resolution acceleration with constrained evolution reconstruction was developed with three-dimensional volume coverage and high-temporal frame rate. Data is acquired using a stack of spirals sampling trajectory combined with a golden ratio view order using an eight-channel coil array. Temporal frames are reconstructed from vastly undersampled data sets using a nonlinear inverse algorithm assuming that the temporal changes are small at short time intervals. Numerical and phantom experimental validation is presented. Preliminary in vivo results demonstrated high spatial resolution dynamic three-dimensional images of the whole liver with high frame rates, from which numerous subarterial phases could be easily identified retrospectively.

Quantitative effects of using compressed sensing in dynamic contrast enhanced MRI

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) involves the acquisition of images before, during and after the injection of a contrast agent. In order to perform quantitative modeling on the resulting signal intensity time course, data must be acquired rapidly, which compromises spatial resolution, signal to noise and/or field of view. One approach that may allow for gains in temporal or spatial resolution or signal to noise of an individual image is to use compressed sensing (CS) MRI. In this study, we demonstrate the accuracy of extracted pharmacokinetic parameters from DCE-MRI data obtained as part of pre-clinical and clinical studies in which fully sampled acquisitions have been retrospectively undersampled by factors of 2, 3 and 4 in Fourier space and then reconstructed with CS. The mean voxel-level concordance correlation coefficient for K(trans) (i.e. the volume transfer constant) obtained from the 2× accelerated and the fully sampled data is 0.92 and 0.90 for mouse and human data, respectively; for 3×, the results are 0.79 and 0.79, respectively; for 4×, the results are 0.64 and 0.70, respectively. The mean error in the tumor mean K(trans) for the mouse and human data at 2× acceleration is 1.8% and -4.2%, respectively; at 3×, 3.6% and -10%, respectively; at 4×, 7.8% and -12%, respectively. These results suggest that CS combined with appropriate reduced acquisitions may be an effective approach to improving image quality in DCE-MRI.

Nonconvex prior image constrained compressed sensing (NCPICCS): theory and simulations on perfusion CT

PURPOSE

To present and evaluate a new image reconstruction method for dynamic CT based on a nonconvex prior image constrained compressed sensing (NCPICCS) algorithm. The authors systematically compared the undersampling potential, functional information recovery, and solution convergence speed of four compressed sensing (CS) based image reconstruction methods using perfusion CT data: Standard l1-based CS, nonconvex CS (NCCS), and l1-based and nonconvex CS, including an additional constraint based on a prior image (PICCS and NCPICCS, respectively).

METHODS

The Shepp-Logan phantom was modified such that its uppermost ellipses changed attenuation through time, simulating both an arterial input function (AIF) and a homogeneous tissue perfusion region. Data were simulated with and without Poisson noise added to the projection data and subsequently reconstructed with all four CS-based methods at four levels of undersampling: 20, 12, 6, and 4 projections. Root mean squared (RMS) error of reconstructed images and recovered time attenuation curves (TACs) were assessed as well as convergence speed. The performance of both PICCS and NCPICCS methods were also evaluated using a kidney perfusion animal experiment data set.

RESULTS

All four CS-based methods were able to reconstruct the phantoms with 20 projections, with similar results on the RMS error of the recovered TACs. NCCS allowed accurate reconstructions with as few as 12 projections, PICCS with as few as six projections, and NCPICCS with as few as four projections. These results were consistent for noise-free and noisy data. NCPICCS required the fewest iterations to converge across all simulation conditions, followed by PICCS, NCCS, and then CS. On animal data, at the lowest level of undersampling tested (16 projections), the image quality of NCPICCS was better than PICCS with fewer streaking artifacts, while the TAC accuracy on the selected region of interest was comparable.

CONCLUSIONS

The authors have presented a novel method for image reconstruction using highly undersampled dynamic CT data. The NCPICCS method takes advantage of the information provided by a prior image, as in PICCS, but employs a more general nonconvex sparsity measure [such as the l(p)-norm (0 < p < or = 1)] rather than the conventional convex l1-norm. Despite the lack of guarantees of a globally optimal solution, the proposed nonconvex extension of PICCS consistently allowed for image reconstruction from fewer samples than the analogous l1-based PICCS method. Both nonconvex sparsity measures as well as prior image information (when available) significantly reduced the number of iterations required for convergence, potentially providing computational advantages for practical implementation of CS-based image reconstruction techniques.