Multi-parametric MRI lesion heterogeneity biomarkers for breast cancer diagnosis.

PURPOSE: To identify intra-lesion imaging heterogeneity biomarkers in multi-parametric Magnetic Resonance Imaging (mpMRI) for breast lesion diagnosis. METHODS: Dynamic Contrast Enhanced (DCE) and Diffusion Weighted Imaging (DWI) of 73 female patients, with 85 histologically verified breast lesions were acquired. Non-rigid multi-resolution registration was utilized to spatially align sequences. Four (4) DCE (2nd post-contrast frame, Initial-Enhancement, Post-Initial-Enhancement and Signal-Enhancement-Ratio) and one (1) DWI (Apparent-Diffusion-Coefficient) representations were analyzed, considering a representative lesion slice. 11 1st-order-statistics and 16 texture features (Gray-Level-Co-occurrence-Matrix (GLCM) and Gray-Level-Run-Length-Matrix (GLRLM) based) were derived from lesion segments, provided by Fuzzy C-Means segmentation, across the 5 representations, resulting in 135 features. Least-Absolute-Shrinkage and Selection-Operator (LASSO) regression was utilized to select optimal feature subsets, subsequently fed into 3 classification schemes: Logistic-Regression (LR), Random-Forest (RF), Support-Vector-Machine-Sequential-Minimal-Optimization (SVM-SMO), assessed with Receiver-Operating-Characteristic (ROC) analysis. RESULTS: LASSO regression resulted in 7, 6 and 7 features subsets from DCE, DWI and mpMRI, respectively. Best classification performance was obtained by the RF multi-parametric scheme (Area-Under-ROC-Curve, (AUC) ± Standard-Error (SE), AUC ± SE = 0.984 ± 0.025), as compared to DCE (AUC ± SE = 0.961 ± 0.030) and DWI (AUC ± SE = 0.938 ± 0.032) and statistically significantly higher as compared to DWI. The selected mpMRI feature subset highlights the significance of entropy (1st-order-statistics and 2nd-order-statistics (GLCM)) and percentile features extracted from 2nd post-contrast frame, PIE, SER maps and ADC map. CONCLUSION: Capturing breast intra-lesion heterogeneity, across mpMRI lesion segments with 1st-order-statistics and texture features (GLCM and GLRLM based), offers a valuable diagnostic tool for breast cancer.

Quantitative assessment of microcalcification cluster image quality in digital breast tomosynthesis, 2-dimensional and synthetic mammography.

Quantitative assessment of microcalcification (MC) cluster image quality is presented, in terms of cluster signal-difference-to-noise ratio (SDNR) intercomparison among digital breast tomosynthesis (DBT) and 2-dimensional (2D) and synthetic-2-dimensional (s2D) mammography. A phantom that provides realistic appearance of MC clusters located in uniform and nonuniform background was imaged in 2D and DBT, considering various scattering conditions. MC cluster SDNR differentiation is investigated with respect to MC particle size (uniform background) and surrounding parenchyma density (nonuniform background). An accurate MC cluster segmentation method was used to delineate individual MC particles and estimate MC cluster SDNR. Analysis of the uniform part of the phantom indicated higher performance of DBT and 2D over s2D for the smallest cluster size (106-177 µm), no difference among mammographic modes for the largest MC cluster (224-354 µm), and enhanced role of 2D for decreasing cluster size and increasing scattering. Analysis of the nonuniform part of the phantom indicated DBT performed better than 2D and s2D in case of dense parenchyma pattern, while 2D and s2D did not differ across parenchyma density patterns and scattering conditions. The presented MC cluster SDNR analysis was capable of revealing subtle differences among mammographic modes and suggests a methodology for clinical image quality assessment. Graphical abstract.

Quantitative visually lossless compression ratio determination of JPEG2000 in digitized mammograms.

The current study presents a quantitative approach towards visually lossless compression ratio (CR) threshold determination of JPEG2000 in digitized mammograms. This is achieved by identifying quantitative image quality metrics that reflect radiologists' visual perception in distinguishing between original and wavelet-compressed mammographic regions of interest containing microcalcification clusters (MCs) and normal parenchyma, originating from 68 images from the Digital Database for Screening Mammography. Specifically, image quality of wavelet-compressed mammograms (CRs, 10:1, 25:1, 40:1, 70:1, 100:1) is evaluated quantitatively by means of eight image quality metrics of different computational principles and qualitatively by three radiologists employing a five-point rating scale. The accuracy of the objective metrics is investigated in terms of (1) their correlation (r) with qualitative assessment and (2) ROC analysis (A z index), employing pooled radiologists' rating scores as ground truth. The quantitative metrics mean square error, mean absolute error, peak signal-to-noise ratio, and structural similarity demonstrated strong correlation with pooled radiologists' ratings (r, 0.825, 0.823, -0.825, and -0.826, respectively) and the highest area under ROC curve (A z , 0.922, 0.920, 0.922, and 0.922, respectively). For each quantitative metric, the highest accuracy values of corresponding ROC curves were used to define metric cut-off values. The metrics cut-off values were subsequently used to suggest a visually lossless CR threshold, estimated to be between 25:1 and 40:1 for the dataset analyzed. Results indicate the potential of the quantitative metrics approach in predicting visually lossless CRs in case of MCs in mammography.

Vessel tree segmentation in presence of interstitial lung disease in MDCT.

The automated segmentation of vessel tree structures is a crucial preprocessing stage in computer aided diagnosis (CAD) schemes of interstitial lung disease (ILD) patterns in multidetector computed tomography (MDCT). The accuracy of such preprocessing stages is expected to influence the accuracy of lung CAD schemes. Although algorithms aimed at improving the accuracy of lung fields segmentation in presence of ILD have been reported, the corresponding vessel tree segmentation stage is under-researched. Furthermore, previously reported vessel tree segmentation methods have only dealt with normal lung parenchyma. In this paper, an automated vessel tree segmentation scheme is proposed, adapted to the presence of pathologies affecting lung parenchyma. The first stage of the method accounts for a recently proposed method utilizing a 3-D multiscale vessel enhancement filter based on eigenvalue analysis of the Hessian matrix and on unsupervised segmentation. The second stage of the method is a texture-based voxel classification refinement to correct possible over-segmentation. The performance of the proposed scheme, and of the previously reported technique, in vessel tree segmentation was evaluated by means of area overlap (previously reported: 0.715 ± 0.082, proposed: 0.931 ± 0.027), true positive fraction (previously reported: 0.968 ± 0.019, proposed: 0.935 ± 0.036) and false positive fraction (previously reported: 0.400 ± 0.181, proposed: 0.074 ± 0.031) on a dataset of 210 axial slices originating from seven ILD affected patient scans (used for performance evaluation out of 15). The proposed method demonstrated a statistically significantly ( p < 0.05) higher performance as compared to the previously reported vessel tree segmentation technique. The impact of suboptimal vessel tree segmentation in a reticular pattern quantification system is also demonstrated.

Texture-based identification and characterization of interstitial pneumonia patterns in lung multidetector CT.

Identification and characterization of diffuse parenchyma lung disease (DPLD) patterns challenges computer-aided schemes in computed tomography (CT) lung analysis. In this study, an automated scheme for volumetric quantification of interstitial pneumonia (IP) patterns, a subset of DPLD, is presented, utilizing a multidetector CT (MDCT) dataset. Initially, lung-field segmentation is achieved by 3-D automated gray-level thresholding combined with an edge-highlighting wavelet preprocessing step, followed by a texture-based border refinement step. The vessel tree volume is identified and removed from lung field, resulting in lung parenchyma (LP) volume. Following, identification and characterization of IP patterns is formulated as a three-class pattern classification of LP into normal, ground glass, and reticular patterns, by means of k-nearest neighbor voxel classification, exploiting 3-D cooccurrence features. Performance of the proposed scheme in indentifying and characterizing ground glass and reticular patterns was evaluated by means of volume overlap (ground glass: 0.734 +/- 0.057, reticular: 0.815 +/- 0.037), true-positive fraction (ground glass: 0.638 +/- 0.055, reticular: 0.942 +/- 0.023) and false-positive fraction (ground glass: 0.361 +/- 0.027, reticular: 0.147 +/- 0.032) on five MDCT scans.

Breast cancer diagnosis: analyzing texture of tissue surrounding microcalcifications.

The current study investigates texture properties of the tissue surrounding microcalcification (MC) clusters on mammograms for breast cancer diagnosis. The case sample analyzed consists of 85 dense mammographic images, originating from the Digital Database for Screening Mammography. Mammograms analyzed contain 100 subtle MC clusters (46 benign and 54 malignant). The tissue surrounding MCs is defined on original and wavelet decomposed images, based on a redundant discrete wavelet transform. Gray-level texture and wavelet coefficient texture features at three decomposition levels are extracted from surrounding tissue regions of interest (ST-ROIs). Specifically, gray-level first-order statistics, gray-level cooccurrence matrices features, and Laws' texture energy measures are extracted from original image ST-ROIs. Wavelet coefficient first-order statistics and wavelet coefficient cooccurrence matrices features are extracted from subimages ST-ROIs. The ability of each feature set in differentiating malignant from benign tissue is investigated using a probabilistic neural network. Classification outputs of most discriminating feature sets are combined using a majority voting rule. The proposed combined scheme achieved an area under receiver operating characteristic curve ( A(z)) of 0.989. Results suggest that MCs' ST texture analysis can contribute to computer-aided diagnosis of breast cancer.