Textural and shape analysis is gaining considerable interest in medical imaging, particularly to identify parameters characterizing tumor heterogeneity and to feed radiomic models. Here, we present a free, multiplatform, and easy-to-use freeware called LIFEx, which enables the calculation of conventional, histogram-based, textural, and shape features from PET, SPECT, MR, CT, and US images, or from any combination of imaging modalities. The application does not require any programming skills and was developed for medical imaging professionals. The goal is that independent and multicenter evidence of the usefulness and limitations of radiomic features for characterization of tumor heterogeneity and subsequent patient management can be gathered. Many options are offered for interactive textural index calculation and for increasing the reproducibility among centers. The software already benefits from a large user community (more than 800 registered users), and interactions within that community are part of the development strategy. Significance: This study presents a user-friendly, multi-platform freeware to extract radiomic features from PET, SPECT, MR, CT, and US images, or any combination of imaging modalities. Cancer Res; 78(16); 4786–9. ©2018 AACR.

Background: Radiomics extracts features from medical images more precisely and more accurately than visual assessment. However, radiomics features are affected by CT scanner parameters such as reconstruction kernel or section thickness, thus obscuring underlying biologically important texture features. Purpose: To investigate whether a compensation method could correct for the variations of radiomic feature values caused by using different CT protocols. Materials and Methods: Phantom data involving 10 texture patterns and 74 patients in cohorts 1 (19 men; 42 patients; mean age, 60.4 years; September-October 2013) and 2 (16 men; 32 patients; mean age, 62.1 years; January-September 2007) scanned by using different CT protocols were retrospectively included. For any radiomic feature, the compensation approach identified a protocol-specific transformation to express all data in a common space that were devoid of protocol effects. The differences in statistical distributions between protocols were assessed by using Friedman tests before and after compensation. Principal component analyses were performed on the phantom data to evaluate the ability to distinguish between texture patterns after compensation. Results: In the phantom data, the statistical distributions of features were different between protocols for all radiomic features and texture patterns (P , .05). After compensation, the protocol effect was no longer detectable (P. .05). Principal component analysis demonstrated that each texture pattern was no longer displayed as different clusters corresponding to different imaging protocols, unlike what was observed before compensation. The correction for scanner effect was confirmed in patient data with 100% (10 of 10 features for cohort 1) and 98% (87 of 89 features for cohort 2) of P values less than .05 before compensation, compared with 30% (three of 10) and 15% (13 of 89) after compensation. Conclusion: Image compensation successfully realigned feature distributions computed from different CT imaging protocols and should facilitate multicenter radiomic studies.

No abstract

Elasticity imaging is based on the measurements of local tissue deformation. The approach to ultrasound elasticity imaging presented in this paper relies on the estimation of dense displacement fields by a coarse-to-fine minimization of an energy function that combines constraints of conservation of echo amplitude and displacement field continuity. The multiscale optimization scheme presents several characteristics aimed at improving and accelerating the convergence of the minimization process. This includes the nonregularized initialization at the coarsest resolution and the use of adaptive configuration spaces. Parameters of the energy model and optimization were adjusted using data obtained from a tissue-like phantom material. Elasticity images from normal in vivo breast tissue were subsequently obtained with these parameters. Introducing a smoothness constraint into motion field estimation helped solve ambiguities due to incoherent motion, leading to elastograms less degraded by decorrelation noise than the ones obtained from correlation-based techniques.

Objective: Test a practical realignment approach to compensate the technical variability of MR radiomic features.Methods: T1 phantom images acquired on 2 scanners, FLAIR and contrast enhanced T1-weighted (CE-T1w) images of 18 brain tumor patients scanned on both 1.5T and 3T scanners, and 36 T2weighted (T2w) images of prostate cancer patients scanned in one of two centers were investigated. The ComBat procedure was used for harmonizing radiomic features. Differences in statistical distributions in feature values between 1.5 and 3T images were tested before and after harmonization. The prostate studies were used to determine the impact of harmonization to distinguish between Gleason grades (GG).Results: In the phantom data, 40 out of 42 radiomic feature values were significantly different between the 2 scanners before harmonization and none after. In white matter regions, the statistical distributions of features were significantly different (P<0.05) between the 1.5 and 3T images for 37 out of 42 features in both FLAIR and CE-T1w images. After harmonization, no statistically significant differences were observed. In brain tumors, 41 (FLAIR) or 36 (CE-T1w) out of 42 features were significantly different between the 1.5 and 3T images without harmonization, against 1 (FLAIR) or none (CE-T1w) with harmonization. In prostate studies, 636 radiomic features were significantly different between GG after harmonization against 461 before. The ability to distinguish between GG using radiomic features was increased after harmonization. Conclusion:ComBat harmonization efficiently removes inter-center technical inconsistencies in radiomic feature values and increases the sensitivity of studies using data from several scanners.

Purpose: To assess if segmentation of the aorta can be accurately achieved using the modulus image of phase contrast (PC) magnetic resonance (MR) acquisitions.Materials and Methods: PC image sequences containing both the ascending and descending aorta of 52 subjects were acquired using three different MR scanners. An automated segmentation technique, based on a 2Dþt deformable surface that takes into account the features of PC aortic images, such as flow-related effects, was developed. The study was designed to: 1) assess the variability of our approach and its robustness to the type of MR scanner, and 2) determine its sensitivity to aortic dilation and its accuracy against an expert manual tracing.Results: Interobserver variability in the lumen area was 0.59 6 0.92% for the automated approach versus 10.09 6 8.29% for manual segmentation. The mean Dice overlap measure was 0.945 6 0.014. The method was robust to the aortic size and highly correlated (r ¼ 0.99) with the manual tracing in terms of aortic area and diameter.Conclusion: A fast and robust automated segmentation of the aortic lumen was developed and successfully tested on images provided by various MR scanners and acquired on healthy volunteers as well as on patients with a dilated aorta.

BackgroundArterial stiffness is considered as an independent predictor of cardiovascular mortality, and is increasingly used in clinical practice. This study aimed at evaluating the consistency of the automated estimation of regional and local aortic stiffness indices from cardiovascular magnetic resonance (CMR) data.ResultsForty-six healthy subjects underwent carotid-femoral pulse wave velocity measurements (CF_PWV) by applanation tonometry and CMR with steady-state free-precession and phase contrast acquisitions at the level of the aortic arch. These data were used for the automated evaluation of the aortic arch pulse wave velocity (Arch_PWV), and the ascending aorta distensibility (AA_Distc, AA_Distb), which were estimated from ascending aorta strain (AA_Strain) combined with either carotid or brachial pulse pressure. The local ascending aorta pulse wave velocity AA_PWVc and AA_PWVb were estimated respectively from these carotid and brachial derived distensibility indices according to the Bramwell-Hill theoretical model, and were compared with the Arch_PWV. In addition, a reproducibility analysis of AA_PWV measurement and its comparison with the standard CF_PWV was performed. Characterization according to the Bramwell-Hill equation resulted in good correlations between Arch_PWV and both local distensibility indices AA_Distc (r = 0.71, p < 0.001) and AA_Distb (r = 0.60, p < 0.001); and between Arch_PWV and both theoretical local indices AA_PWVc (r = 0.78, p < 0.001) and AA_PWVb (r = 0.78, p < 0.001). Furthermore, the Arch_PWV was well related to CF_PWV (r = 0.69, p < 0.001) and its estimation was highly reproducible (inter-operator variability: 7.1%).ConclusionsThe present work confirmed the consistency and robustness of the regional index Arch_PWV and the local indices AA_Distc and AA_Distb according to the theoretical model, as well as to the well established measurement of CF_PWV, demonstrating the relevance of the regional and local CMR indices.

No abstract