Directional and inter-acquisition variability in diffusion-weighted imaging and editing for restricted diffusion.

PURPOSE: To evaluate and quantify inter-directional and inter-acquisition variation in diffusion-weighted imaging (DWI) and emphasize signals that report restricted diffusion to enhance cancer conspicuity, while reducing the effects of local microscopic motion and magnetic field fluctuations. METHODS: Ten patients with biopsy-proven prostate cancer were studied under an Institutional Review Board-approved protocol. Individual acquisitions of DWI signal intensities were reconstructed to calculate inter-acquisition distributions and their statistics, which were compared for healthy versus cancer tissue. A method was proposed to detect and filter the acquisitions affected by motion-induced signal loss. First, signals that reflect restricted diffusion were separated from the acquisitions that suffer from signal loss, likely due to microscopic motion, by imposing a cutoff value. Furthermore, corrected apparent diffusion coefficient maps were calculated by employing a weighted sum of the multiple acquisitions, instead of conventional averaging. These weights were calculated by applying a soft-max function to the set of acquisitions per-voxel, making the analysis immune to acquisitions with significant signal loss, even if the number of such acquisitions is high. RESULTS: Inter-acquisition variation is much larger than the Rician noise variance, local spatial variations, and the estimates of diffusion anisotropy based on the current data, as well as the published values of anisotropy. The proposed method increases the contrast for cancers and yields a sensitivity of 98 . 8 % $$ 98.8\% $$ with a false positive rate of 3 . 9 % $$ 3.9\% $$ . CONCLUSION: Motion-induced signal loss makes conventional signal-averaging suboptimal and can obscure signals from areas with restricted diffusion. Filtering or weighting individual acquisitions prior to image analysis can overcome this problem.

Physically implausible signals as a quantitative quality assessment metric in prostate diffusion-weighted MR imaging.

PURPOSE: To provide a quantitative assessment of diffusion-weighted MR images of the prostate through identification of PIDS which clearly represents artifacts in the data. We calculated the percentage and distribution of PIDS in prostate DWI and compare the amount of PIDS between mpMRI images obtained with and without an endorectal coil. METHODS: This IRB approved retrospective study (from 03/03/2014 to 03/10/2020), included 40 patients scanned with endorectal coil (ERC) and 40 without ER coil (NERC). PIDS contains any voxel where: (1) the diffusion signal increases despite an increase in b-value; and/or (2) apparent diffusion coefficient (ADC) is more than 3.0 µm2/ms (the ADC of pure water at 37 °C and it is physically implausible for any material to have a higher ADC). PIDS for transition zone (TZ) and peripheral zone (PZ) was calculated using an in-house MATLAB program. DWI images were quantitatively inspected for noise, motion, and distortion. T-test was used to compare the difference between PIDS levels in ERC versus NERC and ANOVA to compare the PIDS levels in the anatomic zones. The images were evaluated by a fellowship-trained radiologist in Abdominal Imaging with more than 10 years of experience in reading prostate MRI. This was tested only in prostate in this study. RESULTS: 80 patients (58 ± 8 years old, 80 men) were evaluated. The percentage of voxels exhibiting PIDS was 17.1 ± 8.1% for the ERC cohort and 22.2 ± 15.5% for the NERC cohort. PIDS for NERC versus ERC were not significantly different (p = 0.14). The apex and base showed similar percentages of PIDS in ERC (p = 0.30) and NERC (p = 0.86). The mid (13.8 ± 8.6%) in ERC showed lower values (p = 0.02) of PIDS compared to apex (19.9 ± 11.1%) and base (17.5 ± 8.3%). CONCLUSION: PIDS maps provide a spatially resolved quantitative quality assessment for prostate DWI. Average PIDS over the entire prostate were similar for the ERC and NERC cohorts, and did not differ significantly across prostate zones. However, for many of the patients, PIDS was focally much higher in specific prostate zones. PIDS assessment can guide Radiologist's evaluation of images and the development of improved DWI sequences.

Long-Term Stability of Gradient Characteristics Warrants Model-Based Correction of Diffusion Weighting Bias.

The study aims to test the long-term stability of gradient characteristics for model-based correction of diffusion weighting (DW) bias in an apparent diffusion coefficient (ADC) for multisite imaging trials. Single spin echo (SSE) DWI of a long-tube ice-water phantom was acquired quarterly on six MR scanners over two years for individual diffusion gradient channels, along with B0 mapping, as a function of right-left (RL) and superior-inferior (SI) offsets from the isocenter. Additional double spin-echo (DSE) DWI was performed on two systems. The offset dependences of derived ADC were fit to 4th-order polynomials. Chronic shim gradients were measured from spatial derivatives of B0 maps along the tube direction. Gradient nonlinearity (GNL) was modeled using vendor-provided gradient field descriptions. Deviations were quantified by root-mean-square differences (RMSD), normalized to reference ice-water ADC, between the model and reference (RMSDREF), measurement and model (RMSDEXP), and temporal measurement variations (RMSDTMP). Average RMSDREF was 4.9 ± 3.2 (%RL) and -14.8 ± 3.8 (%SI), and threefold larger than RMSDEXP. RMSDTMP was close to measurement errors (~3%). GNL-induced bias across gradient systems varied up to 20%, while deviation from the model accounted at most for 6.5%, and temporal variation for less than 3% of ADC reproducibility error. Higher SSE RMSDEXP = 7.5-11% was reduced to 2.5-4.8% by DSE, consistent with the eddy current origin. Measured chronic shim gradients below 0.1 mT/m had a minor contribution to ADC bias. The demonstrated long-term stability of spatial ADC profiles and consistency with system GNL models justifies retrospective and prospective DW bias correction based on system gradient design models. Residual errors due to eddy currents and shim gradients should be corrected independent of GNL.

High spectral and spatial resolution MRI of prostate cancer: a pilot study.

PURPOSE: High spectral and spatial resolution (HiSS) MRI is a spectroscopic imaging method focusing on water and fat resonances that has good diagnostic utility in breast imaging. The purpose of this work was to assess the feasibility and potential utility of HiSS MRI for the diagnosis of prostate cancer. METHODS: HiSS MRI was acquired at 3 T from six patients who underwent prostatectomy, yielding a train of 127 phase-coherent gradient echo (GRE) images. In the temporal domain, changes in voxel intensity were analyzed and linear (R) and quadratic (R1, R2) quantifiers of signal logarithm decay were calculated. In the spectral domain, three signal scaling-independent parameters were calculated: water resonance peak width (PW), relative peak asymmetry (PRA), and relative peak distortion from ideal Lorentzian shape (PRD). Seven cancer and five normal tissue regions of interest were identified in correlation with pathology and compared. RESULTS: HiSS-derived quantifiers, except R2, showed high reproducibility (coefficients of variation, 5%-14%). Spectral domain quantifiers performed better than temporal domain quantifiers, with receiver operator characteristic areas under the curve ranging from of 0.83 to 0.91. For temporal domain parameters, the range was 0.74 to 0.91. Low absolute values of the coefficients of correlation between monoexponential decay markers (R, PW) and resonance shape markers (PRA, PRD) were observed (range, 0.23-0.38). CONCLUSION: The feasibility and potential diagnostic utility of HiSS MRI in the prostate at 3 T without an endorectal coil was confirmed. Weak correlation between well-performing markers indicates that complementary information could be leveraged to further improve diagnostic accuracy.

Transforming UTE-mDixon MR Abdomen-Pelvis Images Into CT by Jointly Leveraging Prior Knowledge and Partial Supervision.

Computed tomography (CT) provides information for diagnosis, PET attenuation correction (AC), and radiation treatment planning (RTP). Disadvantages of CT include poor soft tissue contrast and exposure to ionizing radiation. While MRI can overcome these disadvantages, it lacks the photon absorption information needed for PET AC and RTP. Thus, an intelligent transformation from MR to CT, i.e., the MR-based synthetic CT generation, is of great interest as it would support PET/MR AC and MR-only RTP. Using an MR pulse sequence that combines ultra-short echo time (UTE) and modified Dixon (mDixon), we propose a novel method for synthetic CT generation jointly leveraging prior knowledge as well as partial supervision (SCT-PK-PS for short) on large-field-of-view images that span abdomen and pelvis. Two key machine learning techniques, i.e., the knowledge-leveraged transfer fuzzy c-means (KL-TFCM) and the Laplacian support vector machine (LapSVM), are used in SCT-PK-PS. The significance of our effort is threefold: 1) Using the prior knowledge-referenced KL-TFCM clustering, SCT-PK-PS is able to group the feature data of MR images into five initial clusters of fat, soft tissue, air, bone, and bone marrow. Via these initial partitions, clusters needing to be refined are observed and for each of them a few additionally labeled examples are given as the partial supervision for the subsequent semi-supervised classification using LapSVM; 2) Partial supervision is usually insufficient for conventional algorithms to learn the insightful classifier. Instead, exploiting not only the given supervision but also the manifold structure embedded primarily in numerous unlabeled data, LapSVM is capable of training multiple desired tissue-recognizers; 3) Benefiting from the joint use of KL-TFCM and LapSVM, and assisted by the edge detector filter based feature extraction, the proposed SCT-PK-PS method features good recognition accuracy of tissue types, which ultimately facilitates the good transformation from MR images to CT images of the abdomen-pelvis. Applying the method on twenty subjects' feature data of UTE-mDixon MR images, the average score of the mean absolute prediction deviation (MAPD) of all subjects is 140.72 ± 30.60 HU which is statistically significantly better than the 241.36 ± 21.79 HU obtained using the all-water method, the 262.77 ± 42.22 HU obtained using the four-cluster-partitioning (FCP, i.e., external-air, internal-air, fat, and soft tissue) method, and the 197.05 ± 76.53 HU obtained via the conventional SVM method. These results demonstrate the effectiveness of our method for the intelligent transformation from MR to CT on the body section of abdomen-pelvis.

Retrospective Correction of ADC for Gradient Nonlinearity Errors in Multicenter Breast DWI Trials: ACRIN6698 Multiplatform Feasibility Study.

The presented analysis of multisite, multiplatform clinical oncology trial data sought to enhance quantitative utility of the apparent diffusion coefficient (ADC) metric, derived from diffusion-weighted magnetic resonance imaging, by reducing technical interplatform variability owing to systematic gradient nonlinearity (GNL). This study tested the feasibility and effectiveness of a retrospective GNL correction (GNC) implementation for quantitative quality control phantom data, as well as in a representative subset of 60 subjects from the ACRIN 6698 breast cancer therapy response trial who were scanned on 6 different gradient systems. The GNL ADC correction based on a previously developed formalism was applied to trace-DWI using system-specific gradient-channel fields derived from vendor-provided spherical harmonic tables. For quantitative DWI phantom images acquired in typical breast imaging positions, the GNC improved interplatform accuracy from a median of 6% down to 0.5% and reproducibility of 11% down to 2.5%. Across studied trial subjects, GNC increased low ADC (<1 µm2/ms) tumor volume by 16% and histogram percentiles by 5%-8%, uniformly shifting percentile-dependent ADC thresholds by ∼0.06 µm2/ms. This feasibility study lays the grounds for retrospective GNC implementation in multiplatform clinical imaging trials to improve accuracy and reproducibility of ADC metrics used for breast cancer treatment response prediction.

UTE-mDixon-based thorax synthetic CT generation.

PURPOSE: Accurate photon attenuation assessment from MR data remains an unmet challenge in the thorax due to tissue heterogeneity and the difficulty of MR lung imaging. As thoracic tissues encompass the whole physiologic range of photon absorption, large errors can occur when using, for example, a uniform, water-equivalent or a soft-tissue-only approximation. The purpose of this study was to introduce a method for voxel-wise thoracic synthetic CT (sCT) generation from MR data attenuation correction (AC) for PET/MR or for MR-only radiation treatment planning (RTP). METHODS: Acquisition: A radial stack-of-stars combining ultra-short-echo time (UTE) and modified Dixon (mDixon) sequence was optimized for thoracic imaging. The UTE-mDixon pulse sequence collects MR signals at three TE times denoted as UTE, Echo1, and Echo2. Three-point mDixon processing was used to reconstruct water and fat images. Bias field correction was applied in order to avoid artifacts caused by inhomogeneity of the MR magnetic field. ANALYSIS: Water fraction and R2* maps were estimated using the UTE-mDixon data to produce a total of seven MR features, that is UTE, Echo1, Echo2, Dixon water, Dixon fat, Water fraction, and R2*. A feature selection process was performed to determine the optimal feature combination for the proposed automatic, 6-tissue classification for sCT generation. Fuzzy c-means was used for the automatic classification which was followed by voxel-wise attenuation coefficient assignment as a weighted sum of those of the component tissues. Performance evaluation: MR data collected using the proposed pulse sequence were compared to those using a traditional two-point Dixon approach. Image quality measures, including image resolution and uniformity, were evaluated using an MR ACR phantom. Data collected from 25 normal volunteers were used to evaluate the accuracy of the proposed method compared to the template-based approach. Notably, the template approach is applicable here, that is normal volunteers, but may not be robust enough for patients with pathologies. RESULTS: The free breathing UTE-mDixon pulse sequence yielded images with quality comparable to those using the traditional breath holding mDixon sequence. Furthermore, by capturing the signal before T2* decay, the UTE-mDixon image provided lung and bone information which the mDixon image did not. The combination of Dixon water, Dixon fat, and the Water fraction was the most robust for tissue clustering and supported the classification of six tissues, that is, air, lung, fat, soft tissue, low-density bone, and dense bone, used to generate the sCT. The thoracic sCT had a mean absolute difference from the template-based (reference) CT of less than 50 HU and which was better agreement with the reference CT than the results produced using the traditional Dixon-based data. CONCLUSION: MR thoracic acquisition and analyses have been established to automatically provide six distinguishable tissue types to generate sCT for MR-based AC of PET/MR and for MR-only RTP.

Performance of T2 Maps in the Detection of Prostate Cancer.

RATIONALE AND OBJECTIVES: This study compares the performance of T2 maps in the detection of prostate cancer (PCa) in comparison to T2-weighted (T2W) magnetic resonance images. MATERIALS AND METHODS: The prospective study was institutional review board approved. Consenting patients (n = 45) with histologic confirmed PCa underwent preoperative 3-T magnetic resonance imaging with or without endorectal coil. Two radiologists, working independently, marked regions of interests (ROIs) on PCa lesions separately on T2W images and T2 maps. Each ROI was assigned a score of 1-5 based on the confidence in accurately detecting cancer, with 5 being the highest confidence. Subsequently, the histologically confirmed PCa lesions (n = 112) on whole-mount sections were matched with ROIs to calculate sensitivity, positive predictive value (PPV), and radiologist confidence score. Quantitative T2 values of PCa and benign tissue ROIs were measured. RESULTS: Sensitivity and confidence score for PCa detection were similar for T2W images (51%, 4.5 ± 0.8) and T2 maps (52%, 4.5 ± 0.6). However, PPV was significantly higher (P = .001) for T2 maps (88%) compared to T2W (72%) images. The use of endorectal coils nominally improved sensitivity (T2W: 55 vs 47%, T2 map: 54% vs 48%) compared to the use of no endorectal coils, but not the PPV and the confidence score. Quantitative T2 values for PCa (105 ± 28 milliseconds) were significantly (P = 9.3 × 10-14) lower than benign peripheral zone tissue (211 ± 71 milliseconds), with moderate significant correlation with Gleason score (ρ = -0.284). CONCLUSIONS: Our study shows that review of T2 maps by radiologists has similar sensitivity but higher PPV compared to T2W images. Additional quantitative information obtained from T2 maps is helpful in differentiating cancer from normal prostate tissue and determining its aggressiveness.

Diagnosis of Prostate Cancer with Noninvasive Estimation of Prostate Tissue Composition by Using Hybrid Multidimensional MR Imaging: A Feasibility Study.

Purpose To evaluate whether compartmental analysis by using hybrid multidimensional magnetic resonance (MR) imaging can be used to diagnose prostate cancer and determine its aggressiveness. Materials and Methods Twenty-two patients with prostate cancer underwent preoperative 3.0-T MR imaging. Axial images were obtained with hybrid multidimensional MR imaging by using all combinations of echo times (47, 75, 100 msec) and b values of 0, 750, 1500 sec/mm2, resulting in a 3 × 3 array of data associated with each voxel. Volumes of the tissue components stroma, epithelium, and lumen were calculated by fitting the hybrid data to a three-compartment signal model, with distinct, paired apparent diffusion coefficient (ADC) and T2 values associated with each compartment. Volume fractions and conventional ADC and T2 were measured for regions of interest in sites of prostatectomy-verified malignancy (n = 28) and normal tissue (n = 71). Receiver operating characteristic (ROC) analysis was used to evaluate the performance of various parameters in differentiating prostate cancer from benign tissue. Results Compared with normal tissue, prostate cancer showed significantly increased fractional volumes of epithelium (23.2% ± 7.1 vs 48.8% ± 9.2, respectively) and reduced fractional volumes of lumen (26.4% ± 14.1 vs 14.0% ± 5.2) and stroma (50.5% ± 15.7 vs 37.2% ± 9.1) by using hybrid multidimensional MR imaging. The fractional volumes of tissue components show a significantly higher Spearman correlation coefficient with Gleason score (epithelium: ρ = 0.652, P = .0001; stroma: ρ = -0.439, P = .020; lumen: ρ = -0.390, P = .040) compared with traditional T2 values (ρ = -0.292, P = .132) and ADCs (ρ = -0.315, P = .102). The area under the ROC curve for differentiation of cancer from normal prostate was highest for fractional volume of epithelium (0.991), followed by fractional volumes of lumen (0.800) and stroma (0.789). Conclusion Fractional volumes of prostatic lumen, stroma, and epithelium change significantly when cancer is present. These parameters can be measured noninvasively by using hybrid multidimensional MR imaging and have the potential to improve the diagnosis of prostate cancer and determine its aggressiveness. © RSNA, 2018 Online supplemental material is available for this article.

A new design and rationale for 3D orthogonally oversampled k-space trajectories.

A novel center-out 3D trajectory for sampling magnetic resonance data is presented. The trajectory set is based on a single Fermat spiral waveform, which is substantially undersampled in the center of k-space. Multiple trajectories are combined in a "stacked cone" configuration to give very uniform sampling throughout a "hub," which is very efficient in terms of gradient performance and uniform trajectory spacing. The fermat looped, orthogonally encoded trajectories (FLORET) design produces less gradient-efficient trajectories near the poles, so multiple orthogonal hub designs are shown. These multihub designs oversample k-space twice with orthogonal trajectories, which gives unique properties but also doubles the minimum scan time for critical sampling of k-space. The trajectory is shown to be much more efficient than the conventional stack of cones trajectory, and has nearly the same signal-to-noise ratio efficiency (but twice the minimum scan time) as a stack of spirals trajectory. As a center-out trajectory, it provides a shorter minimum echo time than stack of spirals, and its spherical k-space coverage can dramatically reduce Gibbs ringing.

Fast, simple gradient delay estimation for spiral MRI.

Timing delays between data acquisition and gradient transmission result in image degradation. This is especially true in spiral MRI, where delays can alter data in a nonuniform manner, generating significant artifact in the reconstructed data. The many methods that exist to mitigate these delays or measure the k-space coordinates require long measurement times, complicated analysis, specialized phantoms or hardware, or significant changes to the sequence of interest. A fast and simple method is proposed to measure delays on each gradient channel. It requires only minimal modification to an existing spiral sequence and can be used to measure independent delays on three gradient channels and any scan subject within six sequence repetition times. The effectiveness and accuracy of this method are analyzed.

Elliptical field-of-view PROPELLER imaging.

Traditionally two-dimensional scans are designed to support an isotropic field-of-view (iFOV). When imaging elongated objects, significant savings in scan time can potentially be achieved by supporting an elliptical field-of-view (eFOV). This work presents an empirical closed-form solution to adapt the PROPELLER trajectory for an eFOV. The proposed solution is built on the geometry of the PROPELLER trajectory permitting the scan prescription and data reconstruction to remain largely similar to standard PROPELLER. The achieved FOV is experimentally validated by the point spread function (PSF) of a phantom scan. The details of potential savings in scan time and the signal-to-noise ratio (SNR) performance in comparison to iFOV scans for both phantom and in-vivo images are also described.

Functional near-infrared neuroimaging.

Functional near-infrared spectroscopy (fNIR) is a neroimaging modality that enables continuous, noninvasive, and portable monitoring of changes in blood oxygenation and blood volume related to human brain function. Over the last decade, studies in the laboratory have established that fNIR spectroscopy provides a veridical measure of oxygenation and blood flow in the brain. Our recent findings indicate that fNIR can effectively monitor cognitive tasks such as attention, working memory, target categorization, and problem solving. These experimental outcomes compare favorably with functional magnetic resonance imaging (fMRI) studies, and in particular, with the blood oxygenation level dependent signal. Since fNIR can be implemented in the form of a wearable and minimally intrusive device, it has the capacity to monitor brain activity under real life conditions and in everyday environments. Moreover, the fNIR system is amenable to integration with other established physiological and neurobehavioral measures, including electroencephalogram, eye tracking, pupil reflex, heart rate variability, respiration, and electrodermal activity.

Motion artifact cancellation in NIR spectroscopy using Wiener filtering.

We present a Wiener filtering based algorithm for the elimination of motion artifacts present in Near Infrared (NIR) spectroscopy measurements. Until now, adaptive filtering was the only technique used in the noise cancellation in NIR studies. The results in this preliminary study revealed that the proposed method gives better estimates than the classical adaptive filtering approach without the need for additional sensor measurements. Moreover, this novel technique has the potential to filter out motion artifacts in functional near infrared (fNIR) signals, too.