Correcting time-intensity curves in dynamic contrast-enhanced breast MRI for inhomogeneous excitation fields at 7T.

PURPOSE: Inhomogeneous excitation at ultrahigh field strengths (7T and above) compromises the reliability of quantified dynamic contrast-enhanced breast MRI. This can hamper the introduction of ultrahigh field MRI into the clinic. Compensation for this non-uniformity effect can consist of both hardware improvements and post-acquisition corrections. This paper investigated the correctable radiofrequency transmit ( B1+ ) range post-acquisition in both simulations and patient data for 7T MRI. METHODS: Simulations were conducted to determine the minimum B1+ level at which corrections were still beneficial because of noise amplification. Two correction strategies leading to differences in noise amplification were tested. The effect of the corrections on a 7T patient data set (N = 38) with a wide range of B1+ levels was investigated in terms of time-intensity curve types as well as washin, washout and peak enhancement values. RESULTS: In simulations assuming a common amount of T1 saturation, the lowest B1+ level at which the SNR of the corrected images was at least that of the original precontrast image was 43% of the nominal angle. After correction, time-intensity curve types changed in 24% of included patients, and the distribution of curve types corresponded better to the distribution found in literature. Additionally, the overlap between the distributions of washin, washout, and peak enhancement values for grade 1 and grade 2 tumors was slightly reduced. CONCLUSION: Although the correctable range varies with the amount of T1 saturation, post-acquisition correction for inhomogeneous excitation was feasible down to B1+ levels of 43% of the nominal angle in vivo.

Untangling the diffusion signal using the phasor transform.

Separating the decay signal from diffusion-weighted scans into two or more components can be challenging. The phasor technique is well established in the field of optical microscopy for visualization and separation of fluorescent dyes with different lifetimes. The use of the phasor technique for separation of diffusion-weighted decay signals was recently proposed. In this study, we investigate the added value of this technique for fitting decay models and visualization of decay rates. Phasor visualization was performed in five glioblastoma patients. Using simulations, the influence of incorrect diffusivity values and of the number of b-values on fitting a three-component model with fixed diffusivities (dubbed "unmixing") was investigated for both a phasor-based fit and a linear least squares (LLS) fit. Phasor-based intravoxel incoherent motion (IVIM) fitting was compared with nonlinear least squares (NLLS) and segmented fitting (SF) methods in terms of accuracy and precision. The distributions of the parameter estimates of simulated data were compared with those obtained in a healthy volunteer. In the phasor visualizations of two glioblastoma patients, a cluster of points was observed that was not seen in healthy volunteers. The identified cluster roughly corresponded to the enhanced edge region of the tumor of two glioblastoma patients visible on fluid-attenuated inversion recovery (FLAIR) images. For fitting decay models the usefulness of the phasor transform is less pronounced, but the additional knowledge gained from the geometrical configuration of phasor space can aid fitting routines. This has led to slightly improved fitting results for the IVIM model: phasor-based fitting yielded parameter maps with higher precision than the NLLS and SF methods for parameters f and D (interquartile range [IQR] for f: NLLS 27, SF 12, phasor 5.7%; IQR for D: NLLS 0.28, SF 0.18, phasor 0.10 µm2 /s). For unmixing, LLS fitting slightly but consistently outperformed phasor-based fitting in all of the tested scenarios.

Harmonization of Quantitative Parenchymal Enhancement in T1 -Weighted Breast MRI.

BACKGROUND: Differences in imaging parameters influence computer-extracted parenchymal enhancement measures from breast MRI. PURPOSE: To investigate the effect of differences in dynamic contrast-enhanced MRI acquisition parameter settings on quantitative parenchymal enhancement of the breast, and to evaluate harmonization of contrast-enhancement values with respect to flip angle and repetition time. STUDY TYPE: Retrospective. PHANTOM/POPULATIONS: We modeled parenchymal enhancement using simulations, a phantom, and two cohorts (N = 398 and N = 302) from independent cancer centers. SEQUENCE FIELD/STRENGTH: 1.5T dynamic contrast-enhanced T1 -weighted spoiled gradient echo MRI. Vendors: Philips, Siemens, General Electric Medical Systems. ASSESSMENT: We assessed harmonization of parenchymal enhancement in simulations and phantom by varying the MR parameters that influence the amount of T1 -weighting: flip angle (8°-25°) and repetition time (4-12 msec). We calculated the median and interquartile range (IQR) of the enhancement values before and after harmonization. In vivo, we assessed overlap of quantitative parenchymal enhancement in the cohorts before and after harmonization using kernel density estimations. Cohort 1 was scanned with flip angle 20° and repetition time 8 msec; cohort 2 with flip angle 10° and repetition time 6 msec. STATISTICAL TESTS: Paired Wilcoxon signed-rank-test of bootstrapped kernel density estimations. RESULTS: Before harmonization, simulated enhancement values had a median (IQR) of 0.46 (0.34-0.49). After harmonization, the IQR was reduced: median (IQR): 0.44 (0.44-0.45). In the phantom, the IQR also decreased, median (IQR): 0.96 (0.59-1.22) before harmonization, 0.96 (0.91-1.02) after harmonization. Harmonization yielded significantly (P < 0.001) better overlap in parenchymal enhancement between the cohorts: median (IQR) was 0.46 (0.37-0.58) for cohort 1 vs. 0.37 (0.30-0.44) for cohort 2 before harmonization (57% overlap); and 0.35 (0.28-0.43) vs. .0.37 (0.30-0.44) after harmonization (85% overlap). DATA CONCLUSION: The proposed practical harmonization method enables an accurate comparison between patients scanned with differences in imaging parameters. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY STAGE: 4.

Reducing distortions in echo-planar breast imaging at ultrahigh field with high-resolution off-resonance maps.

PURPOSE: DWI is a promising modality in breast MRI, but its clinical acceptance is slow. Analysis of DWI is hampered by geometric distortion artifacts, which are caused by off-resonant spins in combination with the low phase-encoding bandwidth of the EPI sequence used. Existing correction methods assume smooth off-resonance fields, which we show to be invalid in the human breast, where high discontinuities arise at tissue interfaces. METHODS: We developed a distortion correction method that incorporates high-resolution off-resonance maps to better solve for severe distortions at tissue interfaces. The method was evaluated quantitatively both ex vivo in a porcine tissue phantom and in vivo in 5 healthy volunteers. The added value of high-resolution off-resonance maps was tested using a Wilcoxon signed rank test comparing the quantitative results obtained with a low-resolution off-resonance map with those obtained with a high-resolution map. RESULTS: Distortion correction using low-resolution off-resonance maps corrected most of the distortions, as expected. Still, all quantitative comparison metrics showed increased conformity between the corrected EPI images and a high-bandwidth reference scan for both the ex vivo and in vivo experiments. All metrics showed a significant improvement when a high-resolution off-resonance map was used (P < 0.05), in particular at tissue boundaries. CONCLUSION: The use of off-resonance maps of a resolution higher than EPI scans significantly improves upon existing distortion correction techniques, specifically by superior correction at glandular tissue boundaries.

Estimating B1+ in the breast at 7 T using a generic template.

Dynamic contrast-enhanced MRI is the workhorse of breast MRI, where the diagnosis of lesions is largely based on the enhancement curve shape. However, this curve shape is biased by RF transmit (B1+ ) field inhomogeneities. B1+ field information is required in order to correct these. The use of a generic, coil-specific B1+ template is proposed and tested. Finite-difference time-domain simulations for B1+ were performed for healthy female volunteers with a wide range of breast anatomies. A generic B1+ template was constructed by averaging simulations based on four volunteers. Three-dimensional B1+ maps were acquired in 15 other volunteers. Root mean square error (RMSE) metrics were calculated between individual simulations and the template, and between individual measurements and the template. The agreement between the proposed template approach and a B1+ mapping method was compared against the agreement between acquisition and reacquisition using the same mapping protocol. RMSE values (% of nominal flip angle) comparing individual simulations with the template were in the range 2.00-4.01%, with mean 2.68%. RMSE values comparing individual measurements with the template were in the range8.1-16%, with mean 11.7%. The agreement between the proposed template approach and a B1+ mapping method was only slightly worse than the agreement between two consecutive acquisitions using the same mapping protocol in one volunteer: the range of agreement increased from ±16% of the nominal angle for repeated measurement to ±22% for the B1+ template. With local RF transmit coils, intersubject differences in B1+ fields of the breast are comparable to the accuracy of B1+ mapping methods, even at 7 T. Consequently, a single generic B1+ template suits subjects over a wide range of breast anatomies, eliminating the need for a time-consuming B1+ mapping protocol.