Effect of beam attenuation on muscle ultrasound echogenicity measurement in muscular dystrophies.

PURPOSE: In advanced muscular dystrophies (AMD), quantification of muscle echo-intensity (EI) may be influenced by ultrasound beam attenuation, due to fibrosis and fatty infiltration of muscle tissue. Objective of the study was to compare EI measurements using grayscale analysis between a superficial and whole-muscle region of interest (ROI) in subjects with advanced and mild-to-moderate muscular dystrophy (MMD). METHODS: Thirty-two adult subjects diagnosed with a muscular dystrophy and twenty-five matched healthy controls underwent ultrasound assessment of the biceps brachii (BB), rectus femoris (RF) and tibialis anterior (TA) muscles. Based on Heckmatt grading scale of muscles, two disease groups, an AMD (Heckmatt grades 3 or 4) and a MMD (Heckmatt grade 2), were analyzed. Superficial ROI was set as one-fourth of the whole-muscle area, located immediately below the superficial fascia and always inside muscle boundaries. RESULTS: Muscle EI was significantly higher in the superficial compared to whole-muscle ROI, in all evaluated muscle groups of AMD subjects (BB, p = 0.004/RF, p = 0.027/TA, p = 0.002). EI values in superficial ROIs, for individual muscle analysis using z-scores, were more representative in assessments of muscle abnormality in advanced stages of the disease course (Heckmatt grades 3 and 4). In MMD and healthy muscles, no statistical difference was found in EI measurements between the two ROI types. CONCLUSIONS: In AMD, selection of superficial ROI is better representative of changes in muscle texture, although caution should be exercised when comparing ROIs of different sizes.

Incorporating radiomics into clinical trials: expert consensus endorsed by the European Society of Radiology on considerations for data-driven compared to biologically driven quantitative biomarkers.

Existing quantitative imaging biomarkers (QIBs) are associated with known biological tissue characteristics and follow a well-understood path of technical, biological and clinical validation before incorporation into clinical trials. In radiomics, novel data-driven processes extract numerous visually imperceptible statistical features from the imaging data with no a priori assumptions on their correlation with biological processes. The selection of relevant features (radiomic signature) and incorporation into clinical trials therefore requires additional considerations to ensure meaningful imaging endpoints. Also, the number of radiomic features tested means that power calculations would result in sample sizes impossible to achieve within clinical trials. This article examines how the process of standardising and validating data-driven imaging biomarkers differs from those based on biological associations. Radiomic signatures are best developed initially on datasets that represent diversity of acquisition protocols as well as diversity of disease and of normal findings, rather than within clinical trials with standardised and optimised protocols as this would risk the selection of radiomic features being linked to the imaging process rather than the pathology. Normalisation through discretisation and feature harmonisation are essential pre-processing steps. Biological correlation may be performed after the technical and clinical validity of a radiomic signature is established, but is not mandatory. Feature selection may be part of discovery within a radiomics-specific trial or represent exploratory endpoints within an established trial; a previously validated radiomic signature may even be used as a primary/secondary endpoint, particularly if associations are demonstrated with specific biological processes and pathways being targeted within clinical trials. KEY POINTS: • Data-driven processes like radiomics risk false discoveries due to high-dimensionality of the dataset compared to sample size, making adequate diversity of the data, cross-validation and external validation essential to mitigate the risks of spurious associations and overfitting. • Use of radiomic signatures within clinical trials requires multistep standardisation of image acquisition, image analysis and data mining processes. • Biological correlation may be established after clinical validation but is not mandatory.

Transauricular balloon angioplasty in rabbit thoracic aorta: a novel model of experimental restenosis.

BACKGROUND: The aim of this study was to demonstrate a percutaneous transauricular method of balloon angioplasty in high-cholesterol fed rabbits, as an innovative atherosclerosis model. METHODS: Twenty male New Zealand rabbits were randomly divided into two groups of ten animals, as follows: atherogenic diet plus balloon angioplasty (group A) and atherogenic diet alone (group B). Balloon angioplasty was performed in the descending thoracic aorta through percutaneous catheterization of the auricular artery. Eight additional animals fed regular diet were served as long term control. At the end of 9 week period, rabbits were euthanized and thoracic aortas were isolated for histological, immunohistochemical and biochemical analysis. RESULTS: Atherogenic diet induced severe hypercholesterolemia in both group A and B (2802 ± 188.59 and 4423 ± 493.39 mg/dl respectively) compared to the control animals (55.5 ± 11.82 mg/dl; P < 0.001). Group A atherosclerotic lesions appeared to be more advanced histologically (20% type IV and 80% type V) compared to group B lesions (50% type III and 50% type IV). Group A compared to group B atherosclerotic lesions demonstrated similar percentage of macrophages (79.5 ± 9.56% versus 84 ± 12.2%; P = 0.869), more smooth muscle cells (61 ± 14.10% versus 40.5 ± 17.07; P = 0.027), increased intima/media ratio (1.20 ± 0.50 versus 0.62 ± 0.13; P = 0.015) despite the similar degree of intimal hyperplasia (9768 ± 1826.79 µm² versus 12205 ± 8789.23 µm²; P = 0.796), and further significant lumen deterioration (23722 ± 4508.11 versus 41967 ± 20344.61 µm²; P = 0.05) and total vessel area reduction (42350 ± 5819.70 versus 73190 ± 38902.79 µm²; P = 0.022). Group A and B animals revealed similar nitrated protein percentage (P = NS), but significantly higher protein nitration compared to control group (P < 0.01; P < 0.01, respectively). No deaths or systemic complications were reported. CONCLUSION: Transauricular balloon angioplasty constitutes a safe, minimally invasive and highly successful model of induced atherosclerosis in hyperlipidaemic rabbits.

Automated 3D ιnterstitial lung disease εxtent quantification: performance evaluation and correlation to PFTs.

In this study, the performance of a recently proposed computer-aided diagnosis (CAD) scheme in detection and 3D quantification of reticular and ground glass pattern extent in chest computed tomography of interstitial lung disease (ILD) patients is evaluated. CAD scheme performance was evaluated on a dataset of 37 volumetric chest scans, considering five representative axial anatomical levels per scan. CAD scheme reliability analysis was performed by estimating agreement (intraclass correlation coefficient, ICC) of automatically derived ILD pattern extent to semi-quantitative disease extent assessment in terms of 29-point rating scale provided by two expert radiologists. Receiver operating characteristic (ROC) analysis was employed to assess CAD scheme accuracy in ILD pattern detection in terms of area under ROC curve (A z ). Correlation of reticular and ground glass volumetric pattern extent to pulmonary function tests (PFTs) was also investigated. CAD scheme reliability was substantial for ILD extent (ICC = 0.809) and distinct reticular pattern extent (0.806) and moderate for distinct ground glass pattern extent (0.543), performing within inter-observer agreement. CAD scheme demonstrated high accuracy in detecting total ILD (A z = 0.950 ± 0.018), while accuracy in detecting distinct reticular and ground glass patterns was 0.920 ± 0.023 and 0.883 ± 0.024, respectively. Moderate and statistically significant negative correlation was found between reticular volumetric pattern extent and diffusing capacity, forced expiratory volume in 1 s, forced vital capacity, and total lung capacity (R = -0.581, -0.513, -0.494, and -0.446, respectively), similar to correlations found between radiologists' semi-quantitative ratings with PFTs. CAD-based quantification of disease extent is in agreement with radiologists' semi-quantitative assessment and correlates to specific PFTs, suggesting a potential imaging biomarker for ILD staging and management.

The impact of Platelet Rich Plasma (PRP) in osseointegration of oral implants in dental panoramic radiography: texture based evaluation.

PURPOSE: In this study the temporal texture differentiation associated with the bone formation properties, around loaded oral implants after Platelet Rich Plasma (PRP) employment, was investigated in Panoramic Radiographs. MATERIALS AND METHODS: Thirty eligible patients are randomly assigned to two groups. The test group received PRP application around new implants, while in the control group no PRP treatment was made. The bone-to-implant contact region was analyzed in a clinical sample of 60 Digitized Panoramic Radiographs, 30 corresponding to immediate implant loading (Class-I) and 30 after an 8 month follow-up period (Class-II). This region was sampled by 1146 circular Regions-of-Interest (ROIs), resulting from a specifically designed segmentation scheme based on Markov-Random-Fields (MRF). From each ROI, 41 textural features were extracted, then reduced to a subset of 4 features due to redundancy and employed as input to Receiver-Operating-Characteristic (ROC) analysis, to assess the textural differentiation between two classes. RESULTS: The selected subset, achieved Area-Under-Curve (AUC) values ranging from 0.77-0.81 in the PRP group, indicating the significant temporal textural differentiation has been made. In the control group, the AUC values ranged from 0.56-0.68 demonstrating lesser osseo integration activity. CONCLUSION: This study provides evidences that PRP application may favor bone formation around loaded dental implants that could modify the dental treatment planning.

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.

Evaluation of the spline reconstruction technique for preclinical PET imaging.

BACKGROUND AND OBJECTIVE: The Spline Reconstruction Technique (SRT) is a fast algorithm based on a novel numerical implementation of an analytic representation of the inverse Radon transform. The purpose of this study is to provide a comparison between SRT, Filtered Back-Projection (FBP), Ordered Subset Expectation Maximization 2D (2D-OSEM), and the Tera-Tomo 3D algorithm, using phantom data at various acquisition durations as well as small-animal data obtained from the Mediso nanoScan® PET/CT scanner. METHODS: For this purpose, the "NEMA NU 4-2008 standards" protocol was employed at five different realizations and acquisition durations. In addition to the image quality metrics described by the NEMA protocol, Cold Region Contrast was also considered as a figure-of-merit. Furthermore, Cold Region Contrast was measured in the myocardial infarction region of six male Wistar rats. The volumetric defect quantification was assessed with dedicated computer software. Lastly, plots of Recovery Coefficient and Spill-Over Ratio as a function of the Percentage Standard Deviation were generated, after smoothing the phantom reconstructions with four different Gaussian filters. Statistical significance was determined by employing the Kruskal-Wallis test or One-way Analysis of Variance depending on the normality of the variable's distribution. RESULTS: The present study revealed that, at the expense of slightly increased noise in the reconstructed images, SRT resulted in higher Recovery Coefficient values for small hot regions of interest, when compared with FBP and 2D-OSEM at all acquisition durations. Furthermore, SRT reconstructed images exhibit higher Recovery Coefficient values, for all hot regions of interest, when compared to the other 2D algorithms at short acquisition durations. In both phantom and animal studies, SRT achieved a significant improvement over 2D-OSEM for the Spill-Over Ratio and the Cold Region Contrast. These advantages were maintained even after comparing the algorithms at equal noise levels. The Tera-Tomo 3D algorithm (4 subsets, iterations≥ 13) performed significantly better compared to the other algorithms for all figures-of-merit. No statistically significant differences regarding the myocardial defect size were observed between the algorithms investigated. CONCLUSIONS: Overall, SRT appears that could be useful for the quantification of small hot regions of interest, cold regions of interest, as well as in low-count imaging applications.

Standardised lesion segmentation for imaging biomarker quantitation: a consensus recommendation from ESR and EORTC.

BACKGROUND: Lesion/tissue segmentation on digital medical images enables biomarker extraction, image-guided therapy delivery, treatment response measurement, and training/validation for developing artificial intelligence algorithms and workflows. To ensure data reproducibility, criteria for standardised segmentation are critical but currently unavailable. METHODS: A modified Delphi process initiated by the European Imaging Biomarker Alliance (EIBALL) of the European Society of Radiology (ESR) and the European Organisation for Research and Treatment of Cancer (EORTC) Imaging Group was undertaken. Three multidisciplinary task forces addressed modality and image acquisition, segmentation methodology itself, and standards and logistics. Devised survey questions were fed via a facilitator to expert participants. The 58 respondents to Round 1 were invited to participate in Rounds 2-4. Subsequent rounds were informed by responses of previous rounds. RESULTS/CONCLUSIONS: Items with ≥ 75% consensus are considered a recommendation. These include system performance certification, thresholds for image signal-to-noise, contrast-to-noise and tumour-to-background ratios, spatial resolution, and artefact levels. Direct, iterative, and machine or deep learning reconstruction methods, use of a mixture of CE marked and verified research tools were agreed and use of specified reference standards and validation processes considered essential. Operator training and refreshment were considered mandatory for clinical trials and clinical research. Items with a 60-74% agreement require reporting (site-specific accreditation for clinical research, minimal pixel number within lesion segmented, use of post-reconstruction algorithms, operator training refreshment for clinical practice). Items with ≤ 60% agreement are outside current recommendations for segmentation (frequency of system performance tests, use of only CE-marked tools, board certification of operators, frequency of operator refresher training). Recommendations by anatomical area are also specified.

Texture-based quantification of lumbar intervertebral disc degeneration from conventional T2-weighted MRI.

BACKGROUND: Disc degeneration quantification is important for monitoring the effects of new therapeutic methods, such as cell and growth factor therapy. Magnetic resonance (MR) image texture reflects biochemical and structural tissue properties and has been used for differentiating between normal and pathological status in a variety of medical applications. PURPOSE: To investigate the suitability of textural descriptors for the quantification of intervertebral disc degeneration using conventional T2-weighted magnetic resonance images of the lumbar spine. MATERIAL AND METHODS: A 3 Tesla scanner was used, and conventional T2- weighted MR images were obtained, and a total of 255 lumbar discs were analyzed. An atlas-based method was used for segmenting the disc regions from the images. A set of first and second order statistics describing texture of each region were calculated. The validity and reliability of these descriptors for disc degeneration severity quantification was tested through their correlation with patient age and qualitative clinical grading of degeneration severity. Texture quantification results were compared to a widely accepted method for disc degeneration quantification based on the measurement of disc's mean signal intensity. RESULTS: Out of the set of texture descriptors tested, two descriptors quantifying image intensity inhomogeneity, i.e. the grey level standard deviation and co-occurrence derived sum of squares displayed the strongest association to patient age and clinical grading of disc degeneration severity (P<0.001). This is attributed to these inhomogeneity descriptors' capability to capture the progressive loss of nucleus-annulus distinction in the degenerative progress. Statistical analysis indicates that these descriptors can effectively separate between early stages of degeneration. Quantitative measurements are highly repeatable (intraclass correlation >0.98). CONCLUSION: Inhomogeneity descriptors could be a valuable tool for tracking the evolution of disc degeneration and monitoring the response to treatment in a simple, precise and repeatable manner.

Integrating segmentation information into CNN for breast cancer diagnosis of mammographic masses.

Background and ObjectivesSegmentation of mammographic lesions has been proven to be a valuable source of information, as it can assist in both extracting shape-related features and providing accurate localization of the lesion. In this work, a methodology is proposed for integrating mammographic mass segmentation information into a convolutional neural network (CNN), aiming to improve the diagnosis of breast cancer in mammograms. MethodsThe proposed methodology involves modification of each convolutional layer of a CNN, so that information of not only the input image but also the corresponding segmentation map is considered. Furthermore, a new loss function is introduced, which adds an extra term to the standard cross-entropy, aiming to steer the attention of the network to the mass region, penalizing strong feature activations based on their location. The segmentation maps are acquired either from the provided ground-truth or from an automatic segmentation stage. ResultsPerformance evaluation in diagnosis is conducted on two mammographic mass datasets, namely DDSM-400 and CBIS-DDSM, with differences in quality of the corresponding ground-truth segmentation maps. The proposed method achieves diagnosis performance of 0.898 and 0.862 in terms AUC when using ground-truth segmentation maps and a maximum of 0.880 and 0.860 when a U-Net-based automatic segmentation stage is employed, for DDSM-400 and CBIS-DDSM, respectively. ConclusionsThe experimental results demonstrate that integrating segmentation information into a CNN leads to improved performance in breast cancer diagnosis of mammographic masses.

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.

Deep Learning for Breast Cancer Diagnosis from Mammograms-A Comparative Study.

Deep convolutional neural networks (CNNs) are investigated in the context of computer-aided diagnosis (CADx) of breast cancer. State-of-the-art CNNs are trained and evaluated on two mammographic datasets, consisting of ROIs depicting benign or malignant mass lesions. The performance evaluation of each examined network is addressed in two training scenarios: the first involves initializing the network with pre-trained weights, while for the second the networks are initialized in a random fashion. Extensive experimental results show the superior performance achieved in the case of fine-tuning a pretrained network compared to training from scratch.

B-spline active rays segmentation of microcalcifications in mammography.

Accurate segmentation of microcalcifications in mammography is crucial for the quantification of morphologic properties by features incorporated in computer-aided diagnosis schemes. A novel segmentation method is proposed implementing active rays (polar-transformed active contours) on B-spline wavelet representation to identify microcalcification contour point estimates in a coarse-to-fine strategy at two levels of analysis. An iterative region growing method is used to delineate the final microcalcification contour curve, with pixel aggregation constrained by the microcalcification contour point estimates. A radial gradient-based method was also implemented for comparative purposes. The methods were tested on a dataset consisting of 149 mainly pleomorphic microcalcification clusters originating from 130 mammograms of the DDSM database. Segmentation accuracy of both methods was evaluated by three radiologists, based on a five-point rating scale. The radiologists' average accuracy ratings were 3.96 +/- 0.77, 3.97 +/- 0.80, and 3.83 +/- 0.89 for the proposed method, and 2.91 +/- 0.86, 2.10 +/- 0.94, and 2.56 +/- 0.76 for the radial gradient-based method, respectively, while the differences in accuracy ratings between the two segmentation methods were statistically significant (Wilcoxon signed-ranks test, p < 0.05). The effect of the two segmentation methods in the classification of benign from malignant microcalcification clusters was also investigated. A least square minimum distance classifier was employed based on cluster features reflecting three morphological properties of individual microcalcifications (area, length, and relative contrast). Classification performance was evaluated by means of the area under ROC curve (Az). The area and length morphologic features demonstrated a statistically significant (Mann-Whitney U-test, p < 0.05) higher patient-based classification performance when extracted from microcalcifications segmented by the proposed method (0.82 +/- 0.06 and 0.86 +/- .05, respectively), as compared to segmentation by the radial gradient-based method (0.71 +/- 0.08 and 0.75 +/- 0.08). The proposed method demonstrates improved segmentation accuracy, fulfilling human visual criteria, and enhances the ability of morphologic features to characterize microcalcification clusters.

Texture classification-based segmentation of lung affected by interstitial pneumonia in high-resolution CT.

Accurate and automated lung field (LF) segmentation in high-resolution computed tomography (HRCT) is highly challenged by the presence of pathologies affecting lung borders, also affecting the performance of computer-aided diagnosis (CAD) schemes. In this work, a two-dimensional LF segmentation algorithm adapted to interstitial pneumonia (IP) patterns is presented. The algorithm employs k-means clustering followed by a filling operation to obtain an initial LF order estimate. The final LF border is obtained by an iterative support vector machine neighborhood labeling of border pixels based on gray level and wavelet coefficient statistics features. A second feature set based on gray level averaging and gradient features was also investigated to evaluate its effect on segmentation performance of the proposed method. The proposed method is evaluated on a dataset of 22 HRCT cases spanning a range of IP patterns such as ground glass, reticular, and honeycombing. The accuracy of the method is assessed using area overlap and shape differentiation metrics (d(mean), d(rms), and d(max)), by comparing automatically derived lung borders to manually traced ones, and further compared to a gray level thresholding-based (GLT-based) method. Accuracy of the methods evaluated is also compared to interobserver variability. The proposed method incorporating gray level and wavelet coefficient statistics demonstrated the highest segmentation accuracy, averaged over left and right LFs (overlap=0.954, d(mean)=1.080 mm, d(rms)=1.407 mm, and d(max)=4.944 mm), which is statistically significant (two-tailed student's t test for paired data, p<0.0083) with respect to all metrics considered as compared to the proposed method incorporating gray level averaging and gradient features (overlap=0.918, d(mean)=2.354 mm, d(rms)=3.711 mm, and d(max)=14.412 mm) and the GLT-based method (overlap=0.897, d(mean)=3.618 mm, d(rms)=5.007 mm, and d(max)=16.893 mm). The performance of the three segmentation methods, although decreased as IP pattern severity level (mild, moderate, and severe) was increased, did not demonstrate statistically significant difference (two-tailed student's t test for unpaired data, p>0.0167 for all metrics considered). Finally, the accuracy of the proposed method, based on gray level and wavelet coefficient statistics ranges within interobserver variability. The proposed segmentation method could be used as an initial stage of a CAD scheme for IP patterns.

Computer-aided grading and quantification of hip osteoarthritis severity employing shape descriptors of radiographic hip joint space.

A computer-based system was designed for the grading and quantification of hip osteoarthritis (OA) severity. Employing an active-contours segmentation model, 64 hip joint space (HJS) images (18 normal, 46 osteoarthritic) were obtained from the digitized radiographs of 32 unilateral and bilateral OA-patients. Shape features, generated from the HJS-images, and a hierarchical decision tree structure was used for the grading of OA. A shape features based regression model quantified the OA-severity. The system accomplished high accuracies in characterizing hips as "Normal" (100%), of "mild/moderate"-OA (93.8%) or "severe"-OA (96.7%). OA-severity values, as expressed by HJS-narrowing, correlated highly (r=0.9,p<0.001) with the values predicted by the regression model. The system may contribute to OA-patient management.

Selecting registration schemes in case of interstitial lung disease follow-up in CT.

PURPOSE: Primary goal of this study is to select optimal registration schemes in the framework of interstitial lung disease (ILD) follow-up analysis in CT. METHODS: A set of 128 multiresolution schemes composed of multiresolution nonrigid and combinations of rigid and nonrigid registration schemes are evaluated, utilizing ten artificially warped ILD follow-up volumes, originating from ten clinical volumetric CT scans of ILD affected patients, to select candidate optimal schemes. Specifically, all combinations of four transformation models (three rigid: rigid, similarity, affine and one nonrigid: third order B-spline), four cost functions (sum-of-square distances, normalized correlation coefficient, mutual information, and normalized mutual information), four gradient descent optimizers (standard, regular step, adaptive stochastic, and finite difference), and two types of pyramids (recursive and Gaussian-smoothing) were considered. The selection process involves two stages. The first stage involves identification of schemes with deformation field singularities, according to the determinant of the Jacobian matrix. In the second stage, evaluation methodology is based on distance between corresponding landmark points in both normal lung parenchyma (NLP) and ILD affected regions. Statistical analysis was performed in order to select near optimal registration schemes per evaluation metric. Performance of the candidate registration schemes was verified on a case sample of ten clinical follow-up CT scans to obtain the selected registration schemes. RESULTS: By considering near optimal schemes common to all ranking lists, 16 out of 128 registration schemes were initially selected. These schemes obtained submillimeter registration accuracies in terms of average distance errors 0.18 ± 0.01 mm for NLP and 0.20 ± 0.01 mm for ILD, in case of artificially generated follow-up data. Registration accuracy in terms of average distance error in clinical follow-up data was in the range of 1.985-2.156 mm and 1.966-2.234 mm, for NLP and ILD affected regions, respectively, excluding schemes with statistically significant lower performance (Wilcoxon signed-ranks test, p < 0.05), resulting in 13 finally selected registration schemes. CONCLUSIONS: Selected registration schemes in case of ILD CT follow-up analysis indicate the significance of adaptive stochastic gradient descent optimizer, as well as the importance of combined rigid and nonrigid schemes providing high accuracy and time efficiency. The selected optimal deformable registration schemes are equivalent in terms of their accuracy and thus compatible in terms of their clinical outcome.

A digital equalisation technique improving visualisation of dense mammary gland and breast periphery in mammography.

INTRODUCTION: In mammographic imaging, use of high contrast screen-film combinations results in under-exposed and over-exposed film areas corresponding to dense mammary gland and breast periphery (BP), respectively, characterised by degraded image contrast. A digital equalisation technique was designed and developed in order to deal with the problem of poor visualisation of these regions. METHODS AND MATERIAL: The technique is based on the film-digitiser characteristic curve and a layer model of the breast region, as depicted on a mammogram. It remaps each layer grey level (GL) values by a correction factor that accounts for thickness variation in BP and the presence of dense fibroglandular tissues at the mammary gland. The major steps of the technique are segmentation, to isolate the breast region from mammogram background, and adaptive layer GL remapping. RESULTS: The performance of the technique was initially evaluated on a sample of 60 mammograms. Comparative evaluation between the initial and processed images was performed on the basis of nine anatomical features situated at dense mammary gland and BP. The mammographic images resulting from application of the proposed technique are GL equalised and the visualisation improvement of all anatomical features was found to be statistically significant (P<0.05) or highly significant (P<0.0001). The proposed technique was also compared with contrast-limited adaptive histogram equalisation (CLAHE) and found to be more effective in the visualisation of all anatomical features examined, for both dense breast (DB) and BP. DISCUSSION AND CONCLUSION: Application of the proposed technique results in improved visualisation of both dense mammary gland and BP regions. The proposed technique is independent of breast size, breast symmetry and mammographic view. The technique contributes to breast dose minimisation by eliminating the need for a second acquisition.

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.