Analysis of asymmetry in mammograms via directional filtering with Gabor wavelets

Automated Diagnosis of Breast Cancer from Mammogram Using Wavelet, Curvelet Features, and Convolutional Neural Network

Breast cancer is the utmost generally occurring cancer in women and the second most communal cancer. The ground truth standard used in real-time clinical application for the diagnosis is a mammogram. A novel approach is projected in this paper for the automated diagnosis of breast cancer from mammogram images composed from the MIAS data set using curvelet/wavelet transform-based features and a convolutional neural network. The following sequences of operations are involved, namely pre-processing, application of curvelet/wavelet transform, statistical and gray level co-occurrence matrix-based features extracted from curvelet/wavelet coefficients followed by a selection of highly discriminative features by statistical p-test. Initially, pre-trained models VGG16 and VGG19 are used for classification, and Deep convolutional neural network architecture is constructed for which feature matrix is given as input. Pretrained models are used for classification using the concept of transfer learning. The constructed architecture hyperparameters are adjusted, and the highest classification precision of 93% is achieved. The obtained results outperform the state of art methods available in the state of art.

Breast Cancer Detection in Thermal Infrared Images Using Representation Learning and Texture Analysis Methods

Nowadays, breast cancer is one of the most common cancers diagnosed in women. Mammography is the standard screening imaging technique for the early detection of breast cancer. However, thermal infrared images (thermographies) can be used to reveal lesions in dense breasts. In these images, the temperature of the regions that contain tumors is warmer than the normal tissue. To detect that difference in temperature between normal and cancerous regions, a dynamic thermography procedure uses thermal infrared cameras to generate infrared images at fixed time steps, obtaining a sequence of infrared images. In this paper, we propose a novel method to model the changes on temperatures in normal and abnormal breasts using a representation learning technique called learning-to-rank and texture analysis methods. The proposed method generates a compact representation for the infrared images of each sequence, which is then exploited to differentiate between normal and cancerous cases. Our method produced competitive (AUC = 0.989) results when compared to other studies in the literature.

Classification of breast lesions in ultrasonography using sparse logistic regression and morphology-based texture features

PURPOSE

This work proposes a new reliable computer-aided diagnostic (CAD) system for the diagnosis of breast cancer from breast ultrasound (BUS) images. The system can be useful to reduce the number of biopsies and pathological tests, which are invasive, costly, and often unnecessary.

METHODS

The proposed CAD system classifies breast tumors into benign and malignant classes using morphological and textural features extracted from breast ultrasound (BUS) images. The images are first preprocessed to enhance the edges and filter the speckles. The tumor is then segmented semiautomatically using the watershed method. Having the tumor contour, a set of 855 features including 21 shape-based, 810 contour-based, and 24 textural features are extracted from each tumor. Then, a Bayesian Automatic Relevance Detection (ARD) mechanism is used for computing the discrimination power of different features and dimensionality reduction. Finally, a logistic regression classifier computed the posterior probabilities of malignant vs benign tumors using the reduced set of features.

RESULTS

A dataset of 104 BUS images of breast tumors, including 72 benign and 32 malignant tumors, was used for evaluation using an eightfold cross-validation. The algorithm outperformed six state-of-the-art methods for BUS image classification with large margins by achieving 97.12% accuracy, 93.75% sensitivity, and 98.61% specificity rates.

CONCLUSIONS

Using ARD, the proposed CAD system selects five new features for breast tumor classification and outperforms state-of-the-art, making a reliable and complementary tool to help clinicians diagnose breast cancer.

Multidetector Computed Tomography Imaging: Effect of Sparse Sampling and Iterative Reconstruction on Trabecular Bone Microstructure

Multidetector computed tomography-based trabecular bone microstructure analysis ensures promising results in fracture risk prediction caused by osteoporosis. Because multidetector computed tomography is associated with high radiation exposure, its clinical routine use is limited. Hence, in this study, we investigated in 11 thoracic midvertebral specimens whether trabecular texture parameters are comparable derived from (1) images reconstructed using statistical iterative reconstruction (SIR) and filtered back projection as criterion standard at different exposures (80, 150, 220, and 500 mAs) and (2) from SIR-based sparse sampling projections (12.5%, 25%, 50%, and 100%) and equivalent exposures as criterion standard. Twenty-four texture features were computed, and those that showed similar values between (1) filtered back projection and SIR at the different exposure levels and (2) sparse sampling and equivalent exposures and reconstructed with SIR were identified. These parameters can be of equal value in determining trabecular bone microstructure with lower radiation exposure using sparse sampling and SIR.

Effect of radiation dose reduction on texture measures of trabecular bone microstructure: an in vitro study

Osteoporosis is characterized by bone loss and degradation of bone microstructure leading to fracture particularly in elderly people. Osteoporotic bone degeneration and fracture risk can be assessed by bone mineral density and trabecular bone score from 2D projection dual-energy X-ray absorptiometry images. However, multidetector computed tomography image based quantification of trabecular bone microstructure showed significant improvement in prediction of fracture risk beyond that from bone mineral density and trabecular bone score; however, high radiation exposure limits its use in routine clinical in vivo examinations. Hence, this study investigated reduction of radiation dose and its effects on image quality of thoracic midvertebral specimens. Twenty-four texture features were extracted to quantify the image quality from multidetector computed tomography images of 11 thoracic midvertebral specimens, by means of statistical moments, the gray-level co-occurrence matrix, and the gray-level run-length matrix, and were analyzed by an independent sample t-test to observe differences in image texture with respect to radiation doses of 80, 150, 220, and 500 mAs. The results showed that three features-namely, global variance, energy, and run percentage, were not statistically significant ([Formula: see text]) for low doses with respect to 500 mAs. Hence, it is evident that these three dose-independent features can be used for disease monitoring with a low-dose imaging protocol.

Feasibility of opportunistic osteoporosis screening in routine contrast-enhanced multi detector computed tomography (MDCT) using texture analysis

This study investigated the feasibility of opportunistic osteoporosis screening in routine contrast-enhanced MDCT exams using texture analysis. The results showed an acceptable reproducibility of texture features, and these features could discriminate healthy/osteoporotic fracture cohort with an accuracy of 83%.

INTRODUCTION

This aim of this study is to investigate the feasibility of opportunistic osteoporosis screening in routine contrast-enhanced MDCT exams using texture analysis.

METHODS

We performed texture analysis at the spine in routine MDCT exams and investigated the effect of intravenous contrast medium (IVCM) (n = 7), slice thickness (n = 7), the long-term reproducibility (n = 9), and the ability to differentiate healthy/osteoporotic fracture cohort (n = 9 age and gender matched pairs). Eight texture features were extracted using gray level co-occurrence matrix (GLCM). The independent sample t test was used to rank the features of healthy/fracture cohort and classification was performed using support vector machine (SVM).

RESULTS

The results revealed significant correlations between texture parameters derived from MDCT scans with and without IVCM (r up to 0.91) slice thickness of 1 mm versus 2 and 3 mm (r up to 0.96) and scan-rescan (r up to 0.59). The performance of the SVM classifier was evaluated using 10-fold cross-validation and revealed an average classification accuracy of 83%.

CONCLUSIONS

Opportunistic osteoporosis screening at the spine using specific texture parameters (energy, entropy, and homogeneity) and SVM can be performed in routine contrast-enhanced MDCT exams.

Assessment of global and local region-based bilateral mammographic feature asymmetry to predict short-term breast cancer risk

This study aims to develop and test a new imaging marker-based short-term breast cancer risk prediction model. An age-matched dataset of 566 screening mammography cases was used. All 'prior' images acquired in the two screening series were negative, while in the 'current' screening images, 283 cases were positive for cancer and 283 cases remained negative. For each case, two bilateral cranio-caudal view mammograms acquired from the 'prior' negative screenings were selected and processed by a computer-aided image processing scheme, which segmented the entire breast area into nine strip-based local regions, extracted the element regions using difference of Gaussian filters, and computed both global- and local-based bilateral asymmetrical image features. An initial feature pool included 190 features related to the spatial distribution and structural similarity of grayscale values, as well as of the magnitude and phase responses of multidirectional Gabor filters. Next, a short-term breast cancer risk prediction model based on a generalized linear model was built using an embedded stepwise regression analysis method to select features and a leave-one-case-out cross-validation method to predict the likelihood of each woman having image-detectable cancer in the next sequential mammography screening. The area under the receiver operating characteristic curve (AUC) values significantly increased from 0.5863  ±  0.0237 to 0.6870  ±  0.0220 when the model trained by the image features extracted from the global regions and by the features extracted from both the global and the matched local regions (p  =  0.0001). The odds ratio values monotonically increased from 1.00-8.11 with a significantly increasing trend in slope (p  =  0.0028) as the model-generated risk score increased. In addition, the AUC values were 0.6555  ±  0.0437, 0.6958  ±  0.0290, and 0.7054  ±  0.0529 for the three age groups of 37-49, 50-65, and 66-87 years old, respectively. AUC values of 0.6529  ±  0.1100, 0.6820  ±  0.0353, 0.6836  ±  0.0302 and 0.8043  ±  0.1067 were yielded for the four mammography density sub-groups (BIRADS from 1-4), respectively. This study demonstrated that bilateral asymmetry features extracted from local regions combined with the global region in bilateral negative mammograms could be used as a new imaging marker to assist in the prediction of short-term breast cancer risk.

Classifying symmetrical differences and temporal change for the detection of malignant masses in mammography using deep neural networks

We investigate the addition of symmetry and temporal context information to a deep convolutional neural network (CNN) with the purpose of detecting malignant soft tissue lesions in mammography. We employ a simple linear mapping that takes the location of a mass candidate and maps it to either the contralateral or prior mammogram, and regions of interest (ROIs) are extracted around each location. Two different architectures are subsequently explored: (1) a fusion model employing two datastreams where both ROIs are fed to the network during training and testing and (2) a stagewise approach where a single ROI CNN is trained on the primary image and subsequently used as a feature extractor for both primary and contralateral or prior ROIs. A "shallow" gradient boosted tree classifier is then trained on the concatenation of these features and used to classify the joint representation. The baseline yielded an AUC of 0.87 with confidence interval [0.853, 0.893]. For the analysis of symmetrical differences, the first architecture where both primary and contralateral patches are presented during training obtained an AUC of 0.895 with confidence interval [0.877, 0.913], and the second architecture where a new classifier is retrained on the concatenation an AUC of 0.88 with confidence interval [0.859, 0.9]. We found a significant difference between the first architecture and the baseline at high specificity with [Formula: see text]. When using the same architectures to analyze temporal change, we yielded an AUC of 0.884 with confidence interval [0.865, 0.902] for the first architecture and an AUC of 0.879 with confidence interval [0.858, 0.898] in the second setting. Although improvements for temporal analysis were consistent, they were not found to be significant. The results show our proposed method is promising and we suspect performance can greatly be improved when more temporal data become available.

Towards localization of malignant sites of asymmetry across bilateral mammograms

BACKGROUND AND OBJECTIVES

The analysis of patterns of asymmetry between the left and right mammograms of a patient can provide meaningful insights into the presence of an underlying tumor in its early stage. However, the identification of breast cancer by investigating bilateral asymmetry is difficult to perform due to the indistinct and borderline nature of the asymmetric signs as they appear on mammograms.

METHODS

In this study, to increase the positive-predictive value of asymmetry in mammographic screening, a novel computerized approach for the automatic localization of malignant sites of asymmetry in mammograms is proposed. The sites of anatomical correspondence between the right and left regions of each radiographic projection were extracted by means of two bilateral masking procedures, inspired by radiologists' criteria in interpreting mammograms and based on the use of detected landmarking structures. Relative variations of spatial patterns of intensity values and of orientations of directional components within each site were quantified by combining multidirectional Gabor filters and indices of structural similarity. The localization of the sites of malignant asymmetry was performed by coupling two quadratic discriminant analysis classifiers, one for each masking procedure, that assigned the likelihood of malignancy to each site of correspondence.

RESULTS

The performance of the proposed method was assessed on 94 mammographic images from two publicly available databases and containing at least one asymmetric site. Sensitivity, specificity and balanced accuracy levels of 0.83 (0.09), 0.75 (0.06), and 0.79 (0.04), respectively were obtained in the classification of malignant asymmetric sites vs benign/normal sites using cross-validation. In addition, a further blind test on a dataset of Full Field Digital Mammograms achieved levels of sensitivity, specificity, and balanced accuracy of 0.86, 0.65, and 0.75, respectively.

CONCLUSIONS

The achieved performance indicates that the proposed system is effective in localizing sites of malignant asymmetry and it is expected to improve computer-aided diagnosis of breast cancer.

A new computer-aided detection scheme based on assessment of local bilateral mammographic feature asymmetry - a preliminary evaluation

Accurate segmentation of breast lesions depicting on two-dimensional projection mammograms has been proven very difficult and unreliable. In this study we investigated a new approach of a computer-aided detection (CAD) scheme of mammograms without lesion segmentation. Our scheme was developed based on the detection and analysis of region-of-interest (ROI)-based bilateral mammographic tissue or feature asymmetry. A bilateral image registration, image feature selection process, and naïve Bayes linear classifier were implemented in CAD scheme. CAD performance predicting the likelihood of either an ROI or a subject (case) being abnormal was evaluated using 161 subjects from the mini-MIAS database and a leave-one-out testing method. The results showed that areas under receiver operating characteristic (ROC) curves were 0.87 and 0.72 on the ROI-based and case-based evaluation, respectively. The study demonstrated that using ROI-based bilateral mammographic tissue asymmetry can provide supplementary information with high discriminatory power in order to improve CAD performance.

AUTOMATED CHARACTERIZATION AND DETECTION OF DIABETIC RETINOPATHY USING TEXTURE MEASURES

The chronic and uncontrolled diabetes mellitus (DM) damages the retinal blood vessels leading to diabetic retinopathy (DR). The advanced stage of DR leads to loss of vision and subsequently blindness. The morphological changes during the progression of DR can be diagnosed using digital fundus images. The pathological changes in the retina influence the variations in pixel patterns which can be quantified using texture measures. In this paper, we have explored different texture measures namely statistical moments, gray level co-occurrence matrix (GLCM), gray level run length matrix (GLRLM), local binary pattern (LBP), laws mask energy (LME), fractal dimension (FD), fourier spectrum (FS) and Gabor wavelet to characterize and classify the normal and DR classes. We have tabulated 109 texture parameters for the normal and DR classes. Further, these features were subjected to empirical receiver operating characteristic (ROC) based ranking to select optimal feature set. The ranked nested features were fed to the support vector machine (SVM) classifier with different kernel functions to evaluate the highest performance measure using the least number of features to discriminate normal and DR classes. Our proposed system was evaluated using two different databases Kasturba Medical College Hospital (KMCH) and Tan Tock Seng Hospital (TTSH), each with 340 images (170 normal and 170 DR). We have also formulated an integrated index called as diabetic retinopathy risk index (DRRI) using selected texture features to discriminate normal and DR classes using single number. The proposed frame work can be used to help the clinicians and also for mass DR screening programs.

Analysis of structural similarity in mammograms for detection of bilateral asymmetry

We hypothesize that quantification of structural similarity or dissimilarity between paired mammographic regions can be effective in detecting asymmetric signs of breast cancer. Bilateral masking procedures are applied for this purpose by using automatically detected anatomical landmarks. Changes in structural information of the extracted regions are investigated using spherical semivariogram descriptors and correlation-based structural similarity indices in the spatial and complex wavelet domains. The spatial distribution of grayscale values as well as of the magnitude and phase responses of multidirectional Gabor filters are used to represent the structure of mammographic density and of the directional components of breast tissue patterns, respectively. A total of 188 mammograms from the DDSM and mini-MIAS databases, consisting of 47 asymmetric cases and 47 normal cases, were analyzed. For the combined dataset of mammograms, areas under the receiver operating characteristic curves of 0.83, 0.77, and 0.87 were obtained, respectively, with linear discriminant analysis, the Bayesian classifier, and an artificial neural network with radial basis functions, using the features selected by stepwise logistic regression and leave-one-patient-out cross-validation. Two-view analysis provided accuracy up to 0.94, with sensitivity and specificity of 1 and 0.88, respectively.

Preparation of 2D sequences of corneal images for 3D model building

A confocal microscope provides a sequence of images, at incremental depths, of the various corneal layers and structures. From these, medical practioners can extract clinical information on the state of health of the patient's cornea. In this work we are addressing problems associated with capturing and processing these images including blurring, non-uniform illumination and noise, as well as the displacement of images laterally and in the anterior-posterior direction caused by subject movement. The latter may cause some of the captured images to be out of sequence in terms of depth. In this paper we introduce automated algorithms for classification, reordering, registration and segmentation to solve these problems. The successful implementation of these algorithms could open the door for another interesting development, which is the 3D modelling of these sequences.

Detection of architectural distortion in prior mammograms via analysis of oriented patterns

We demonstrate methods for the detection of architectural distortion in prior mammograms of interval-cancer cases based on analysis of the orientation of breast tissue patterns in mammograms. We hypothesize that architectural distortion modifies the normal orientation of breast tissue patterns in mammographic images before the formation of masses or tumors. In the initial steps of our methods, the oriented structures in a given mammogram are analyzed using Gabor filters and phase portraits to detect node-like sites of radiating or intersecting tissue patterns. Each detected site is then characterized using the node value, fractal dimension, and a measure of angular dispersion specifically designed to represent spiculating patterns associated with architectural distortion. Our methods were tested with a database of 106 prior mammograms of 56 interval-cancer cases and 52 mammograms of 13 normal cases using the features developed for the characterization of architectural distortion, pattern classification via quadratic discriminant analysis, and validation with the leave-one-patient out procedure. According to the results of free-response receiver operating characteristic analysis, our methods have demonstrated the capability to detect architectural distortion in prior mammograms, taken 15 months (on the average) before clinical diagnosis of breast cancer, with a sensitivity of 80% at about five false positives per patient.

Hybrid Discrete Wavelet Transform and Gabor Filter Banks Processing for Features Extraction from Biomedical Images

A new methodology for automatic feature extraction from biomedical images and subsequent classification is presented. The approach exploits the spatial orientation of high-frequency textural features of the processed image as determined by a two-step process. First, the two-dimensional discrete wavelet transform (DWT) is applied to obtain the HH high-frequency subband image. Then, a Gabor filter bank is applied to the latter at different frequencies and spatial orientations to obtain new Gabor-filtered image whose entropy and uniformity are computed. Finally, the obtained statistics are fed to a support vector machine (SVM) binary classifier. The approach was validated on mammograms, retina, and brain magnetic resonance (MR) images. The obtained classification accuracies show better performance in comparison to common approaches that use only the DWT or Gabor filter banks for feature extraction.

Pectoral muscle segmentation: a review

Mammograms are X-ray images of breasts which are used to detect breast cancer. The pectoral muscle is a mass of tissue on which the breast rests. During routine mammographic screenings, in medio-lateral oblique (MLO) views, the pectoral muscle turns up in the mammograms along with the breast tissues. The pectoral muscle has to be segmented from the mammogram for an effective automated computer aided diagnosis (CAD). This is due to the fact that pectoral muscles have pixel intensities and texture similar to that of breast tissues which can result in awry CAD results. As a result, a lot of effort has been put into the segmentation of pectoral muscles and finding its contour with the breast tissues. To the best of our knowledge, currently there is no definitive literature available which provides a comprehensive review about the current state of research in this area of pectoral muscle segmentation. We try to address this shortcoming by providing a comprehensive review of research papers in this area. A conscious effort has been made to avoid deviating into the area of automated breast cancer detection.

Automated oral cancer identification using histopathological images: a hybrid feature extraction paradigm

Oral cancer (OC) is the sixth most common cancer in the world. In India it is the most common malignant neoplasm. Histopathological images have widely been used in the differential diagnosis of normal, oral precancerous (oral sub-mucous fibrosis (OSF)) and cancer lesions. However, this technique is limited by subjective interpretations and less accurate diagnosis. The objective of this work is to improve the classification accuracy based on textural features in the development of a computer assisted screening of OSF. The approach introduced here is to grade the histopathological tissue sections into normal, OSF without Dysplasia (OSFWD) and OSF with Dysplasia (OSFD), which would help the oral onco-pathologists to screen the subjects rapidly. The biopsy sections are stained with H&E. The optical density of the pixels in the light microscopic images is recorded and represented as matrix quantized as integers from 0 to 255 for each fundamental color (Red, Green, Blue), resulting in a M×N×3 matrix of integers. Depending on either normal or OSF condition, the image has various granular structures which are self similar patterns at different scales termed "texture". We have extracted these textural changes using Higher Order Spectra (HOS), Local Binary Pattern (LBP), and Laws Texture Energy (LTE) from the histopathological images (normal, OSFWD and OSFD). These feature vectors were fed to five different classifiers: Decision Tree (DT), Sugeno Fuzzy, Gaussian Mixture Model (GMM), K-Nearest Neighbor (K-NN), Radial Basis Probabilistic Neural Network (RBPNN) to select the best classifier. Our results show that combination of texture and HOS features coupled with Fuzzy classifier resulted in 95.7% accuracy, sensitivity and specificity of 94.5% and 98.8% respectively. Finally, we have proposed a novel integrated index called Oral Malignancy Index (OMI) using the HOS, LBP, LTE features, to diagnose benign or malignant tissues using just one number. We hope that this OMI can help the clinicians in making a faster and more objective detection of benign/malignant oral lesions.

Automated Diagnosis of Oral Cancer Using Higher Order Spectra Features and Local Binary Pattern: A Comparative Study

In the field of quantitative microscopy, textural information plays a significant role very often in tissue characterization and diagnosis, in addition to morphology and intensity. The objective of this work is to improve the classification accuracy based on textural features for the development of a computer assisted screening of oral sub-mucous fibrosis (OSF). In fact, the approach introduced is used to grade the histopathological tissue sections into normal, OSF without dysplasia (OSFWD) and OSF with dysplasia (OSFD), which would help the oral onco-pathologists to screen the subjects rapidly. The main objective of this work is to evaluate the use of Higher Order Spectra (HOS) features and Local Binary Pattern (LBP) features extracted from the epithelial layer in classifying normal, OSFWD and OSFD. For this purpose, we extracted twenty three HOS features and nine LBP features and fed them to a Support Vector Machine (SVM) for automated diagnosis. One hundred and fifty eight images (90 normal, 42 OSFWD and 26 OSFD images) were used for analysis. LBP features provide a good sensitivity of 82.85% and specificity of 87.84%, and the HOS features provide higher values of sensitivity (94.07%) and specificity (93.33%) using SVM classifier. The proposed system, can be used as an adjunct tool by the onco-pathologists to cross-check their diagnosis.

A fully automated scheme for mammographic segmentation and classification based on breast density and asymmetry

This paper presents a fully automated segmentation and classification scheme for mammograms, based on breast density estimation and detection of asymmetry. First, image preprocessing and segmentation techniques are applied, including a breast boundary extraction algorithm and an improved version of a pectoral muscle segmentation scheme. Features for breast density categorization are extracted, including a new fractal dimension-related feature, and support vector machines (SVMs) are employed for classification, achieving accuracy of up to 85.7%. Most of these properties are used to extract a new set of statistical features for each breast; the differences among these feature values from the two images of each pair of mammograms are used to detect breast asymmetry, using an one-class SVM classifier, which resulted in a success rate of 84.47%. This composite methodology has been applied to the miniMIAS database, consisting of 322 (MLO) mammograms -including 15 asymmetric pairs of images-, obtained via a (noisy) digitization procedure. The results were evaluated by expert radiologists and are very promising, showing equal or higher success rates compared to other related works, despite the fact that some of them used only selected portions of this specific mammographic database. In contrast, our methodology is applied to the complete miniMIAS database and it exhibits the reliability that is normally required for clinical use in CAD systems.

Computer-aided detection of architectural distortion in prior mammograms of interval cancer

Architectural distortion is an important sign of breast cancer, but because of its subtlety, it is a common cause of false-negative findings on screening mammograms. This paper presents methods for the detection of architectural distortion in mammograms of interval cancer cases taken prior to the detection of breast cancer using Gabor filters, phase portrait analysis, fractal analysis, and texture analysis. The methods were used to detect initial candidates for sites of architectural distortion in prior mammograms of interval cancer and also normal control cases. A total of 4,224 regions of interest (ROIs) were automatically obtained from 106 prior mammograms of 56 interval cancer cases, including 301 ROIs related to architectural distortion, and from 52 prior mammograms of 13 normal cases. For each ROI, the fractal dimension and Haralick's texture features were computed. Feature selection was performed separately using stepwise logistic regression and stepwise regression. The best results achieved, in terms of the area under the receiver operating characteristics curve, with the features selected by stepwise logistic regression are 0.76 with the Bayesian classifier, 0.73 with Fisher linear discriminant analysis, 0.77 with an artificial neural network based on radial basis functions, and 0.77 with a support vector machine. Analysis of the performance of the methods with free-response receiver operating characteristics indicated a sensitivity of 0.80 at 7.6 false positives per image. The methods have good potential in detecting architectural distortion in mammograms of interval cancer cases.

Automatic classification for pathological prostate images based on fractal analysis

Accurate grading for prostatic carcinoma in pathological images is important to prognosis and treatment planning. Since human grading is always time-consuming and subjective, this paper presents a computer-aided system to automatically grade pathological images according to Gleason grading system which is the most widespread method for histological grading of prostate tissues. We proposed two feature extraction methods based on fractal dimension to analyze variations of intensity and texture complexity in regions of interest. Each image can be classified into an appropriate grade by using Bayesian, k-NN, and support vector machine (SVM) classifiers, respectively. Leave-one-out and k-fold cross-validation procedures were used to estimate the correct classification rates (CCR). Experimental results show that 91.2%, 93.7%, and 93.7% CCR can be achieved by Bayesian, k-NN, and SVM classifiers, respectively, for a set of 205 pathological prostate images. If our fractal-based feature set is optimized by the sequential floating forward selection method, the CCR can be promoted up to 94.6%, 94.2%, and 94.6%, respectively, using each of the above three classifiers. Experimental results also show that our feature set is better than the feature sets extracted from multiwavelets, Gabor filters, and gray-level co-occurrence matrix methods because it has a much smaller size and still keeps the most powerful discriminating capability in grading prostate images.

Automatic detection of the breast border and nipple position on digital mammograms using genetic algorithm for asymmetry approach to detection of microcalcifications

The presence of microcalcifications in breast tissue is one of the most incident signs considered by radiologist for an early diagnosis of breast cancer, which is one of the most common forms of cancer among women. In this paper, the Genetic Algorithm (GA) is proposed for automatic look at commonly prone area the breast border and nipple position to discover the suspicious regions on digital mammograms based on asymmetries between left and right breast image. The basic idea of the asymmetry approach is to scan left and right images are subtracted to extract the suspicious region. The proposed system consists of two steps: First, the mammogram images are enhanced using median filter, normalize the image, at the pectoral muscle region is excluding the border of the mammogram and comparing for both left and right images from the binary image. Further GA is applied to magnify the detected border. The figure of merit is calculated to evaluate whether the detected border is exact or not. And the nipple position is identified using GA. The some comparisons method is adopted for detection of suspected area. Second, using the border points and nipple position as the reference the mammogram images are aligned and subtracted to extract the suspicious region. The algorithms are tested on 114 abnormal digitized mammograms from Mammogram Image Analysis Society database.

Radon-domain detection of the nipple and the pectoral muscle in mammograms

In this paper, methods are presented for automatic detection of the nipple and the pectoral muscle edge in mammograms via image processing in the Radon domain. Radon-domain information was used for the detection of straight-line candidates with high gradient. The longest straight-line candidate was used to identify the pectoral muscle edge. The nipple was detected as the convergence point of breast tissue components, indicated by the largest response in the Radon domain. Percentages of false-positive (FP) and false-negative (FN) areas were determined by comparing the areas of the pectoral muscle regions delimited manually by a radiologist and by the proposed method applied to 540 mediolateral-oblique (MLO) mammographic images. The average FP and FN were 8.99% and 9.13%, respectively. In the detection of the nipple, an average error of 7.4 mm was obtained with reference to the nipple as identified by a radiologist on 1,080 mammographic images (540 MLO and 540 craniocaudal views).

Gabor filters and phase portraits for the detection of architectural distortion in mammograms

Architectural distortion is a subtle abnormality in mammograms, and a source of overlooking errors by radiologists. Computer-aided diagnosis (CAD) techniques can improve the performance of radiologists in detecting masses and calcifications; however, most CAD systems have not been designed to detect architectural distortion. We present a new method to detect and localise architectural distortion by analysing the oriented texture in mammograms. A bank of Gabor filters is used to obtain the orientation field of the given mammogram. The curvilinear structures (CLS) of interest (spicules and fibrous tissue) are separated from confounding structures (pectoral muscle edge, parenchymal tissue edges, breast boundary, and noise). The selected core CLS pixels and the orientation field are filtered and downsampled, to reduce noise and also to reduce the computational effort required by the subsequent methods. The downsampled orientation field is analysed to produce three phase portrait maps: node, saddle, and spiral. The node map is further analysed in order to detect the sites of architectural distortion. The method was tested with 19 mammograms containing architectural distortion. In a preliminary experiment, a sensitivity of 84% was obtained at 7.8 false positives per image.

Linear structures in mammographic images: detection and classification

We describe methods for detecting linear structures in mammograms, and for classifying them into anatomical types (vessels, spicules, ducts, etc). Several different detection methods are compared, using realistic synthetic images and receiver operating characteristic (ROC) analysis. There are significant differences (p < 0.001) between the methods, with the best giving an Az value for pixel-level detection of 0.943. We also investigate methods for classifying the detected linear structures into anatomical types, using their cross-sectional profiles, with particular emphasis on recognising the "spicules" and "ducts" associated with some of the more subtle abnormalities. Automatic classification results are compared with expert annotations using ROC analysis, demonstrating useful discrimination between anatomical classes (Az = 0.746). Some of this discrimination relies on simple attributes such as profile width and contrast, but important information is also carried by the shape of the profile (Az = 0.653). The methods presented have potentially wide application in improving the specificity of abnormality detection by exploiting additional anatomical information.

Segmentation of the fibro-glandular disc in mammograms using Gaussian mixture modelling

The paper presents a technique for the segmentation of the fibro-glandular disc in mammograms based upon a statistical model of breast density. The density function of the model was represented by a mixture of up to four weighted Gaussians, each one corresponding to a specific density class in the breast. The parameters of the model and the number of tissue classes in the breast were determined using the expectation-maximisation algorithm and the minimum description length method. Grey-level statistics of the pectoral muscle were used to determine the tissue categories that are likely to represent the fibro-glandular disc. The method was applied to 84 medio-lateral oblique mammograms from the Mini-MIAS database. The results of the segmented fibro-glandular disc were assessed by a radiologist using the original and the segmented images, with reference to a ranking table categorising the results of segmentation as: 1: excellent; 2: good; 3: average; 4: poor; and 5: complete failure. Of the 84 cases analysed, 64.3% were rated as excellent, 16.7% were rated as good, 10.7% were rated as average, and 4.7% were rated as poor; only 3.6% of the cases were rated as a complete failure with regard to segmentation of the fibro-glandular disc.

Estimation of the tissue composition of the tumour mass in neuroblastoma using segmented CT images

Neuroblastoma is the most common extra-cranial, solid, malignant tumour in children. Advances in radiology have made possible the detection and staging of the disease. Nevertheless, there is no method available at present that can go beyond detection and qualitative analysis, towards quantitative assessment of the tissue composition of the primary tumour mass in neuroblastoma. Such quantitative analysis could provide important information and serve as a decision-support tool to the radiologist and the oncologist, result in better treatment and follow-up and even lead to the avoidance of delayed surgery. The problem investigated was the improvement of the analysis of the primary tumour mass, in patients with neuroblastoma, using X-ray computed tomography (CT) images. A methodology was proposed for the estimation of the tissue content of the mass: it comprised a Gaussian mixture model for estimation, from segmented CT images, of the tissue composition of the primary tumour. To demonstrate the potential of the method, the results are presented of its application to ten CT examinations of four patients. The method provides quantitative information, and it was observed that the tumour in one of the patients reduced from 523 cm3 to 81 cm3 in volume, with an increase in calcification from about 20% to about 88% of the tumour volume, in response to chemotherapy over a period of five months. Results indicate that the proposed technique may be of considerable value in assessing the response to therapy of patients with neuroblastoma.

Identification of the breast boundary in mammograms using active contour models

A method for the identification of the breast boundary in mammograms is presented. The method can be used in the preprocessing stage of a system for computer-aided diagnosis (CAD) of breast cancer and also in the reduction of image file size in picture archiving and communication system applications. The method started with modification of the contrast of the original image. A binarisation procedure was then applied to the image, and the chain-code algorithm was used to find an approximate breast contour. Finally, the identification of the true breast boundary was performed by using the approximate contour as the input to an active contour model algorithm specially tailored for this purpose. After demarcation of the breast boundary, all artifacts outside the breast region were eliminated. The method was applied to 84 medio-lateral oblique mammograms from the Mini-MIAS database. Evaluation of the detected breast boundary was performed based upon the percentage of false-positive and false-negative pixels determined by a quantitative comparison between the contours identified by a radiologist and those identified by the proposed method. The average false positive and false negative rates were 0.41% and 0.58%, respectively. The two radiologists who evaluated the results considered the segmentation results to be acceptable for CAD purposes.

Automatic identification of the pectoral muscle in mammograms

The pectoral muscle represents a predominant density region in most medio-lateral oblique (MLO) views of mammograms; its inclusion can affect the results of intensity-based image processing methods or bias procedures in the detection of breast cancer. Local analysis of the pectoral muscle may be used to identify the presence of abnormal axillary lymph nodes, which may be the only manifestation of occult breast carcinoma. We propose a new method for the identification of the pectoral muscle in MLO mammograms based upon a multiresolution technique using Gabor wavelets. This new method overcomes the limitation of the straight-line representation considered in our initial investigation using the Hough transform. The method starts by convolving a group of Gabor filters, specially designed for enhancing the pectoral muscle edge, with the region of interest containing the pectoral muscle. After computing the magnitude and phase images using a vector-summation procedure, the magnitude value of each pixel is propagated in the direction of the phase. The resulting image is then used to detect the relevant edges. Finally, a post-processing stage is used to find the true pectoral muscle edge. The method was applied to 84 MLO mammograms from the Mini-MIAS (Mammographic Image Analysis Society, London, U.K.) database. Evaluation of the pectoral muscle edge detected in the mammograms was performed based upon the percentage of false-positive (FP) and false-negative (FN) pixels determined by comparison between the numbers of pixels enclosed in the regions delimited by the edges identified by a radiologist and by the proposed method. The average FP and FN rates were, respectively, 0.58% and 5.77%. Furthermore, the results of the Gabor-filter-based method indicated low Hausdorff distances with respect to the hand-drawn pectoral muscle edges, with the mean and standard deviation being 3.84 +/- 1.73 mm over 84 images.