Research Article Statistical Analysis of Haralick Texture Features to Discriminate Lung Abnormalities

Similar documents
Study and Comparison of Various Techniques of Image Edge Detection

Journal of Engineering Science and Technology Review 11 (2) (2018) Research Article

Copy Number Variation Methods and Data

*VALLIAPPAN Raman 1, PUTRA Sumari 2 and MANDAVA Rajeswari 3. George town, Penang 11800, Malaysia. George town, Penang 11800, Malaysia

Improvement of Automatic Hemorrhages Detection Methods using Brightness Correction on Fundus Images

FAST DETECTION OF MASSES IN MAMMOGRAMS WITH DIFFICULT CASE EXCLUSION

Using the Perpendicular Distance to the Nearest Fracture as a Proxy for Conventional Fracture Spacing Measures

A Computer-aided System for Discriminating Normal from Cancerous Regions in IHC Liver Cancer Tissue Images Using K-means Clustering*

Gene Selection Based on Mutual Information for the Classification of Multi-class Cancer

A New Machine Learning Algorithm for Breast and Pectoral Muscle Segmentation

Modeling Multi Layer Feed-forward Neural. Network Model on the Influence of Hypertension. and Diabetes Mellitus on Family History of

A Geometric Approach To Fully Automatic Chromosome Segmentation

AUTOMATED DETECTION OF HARD EXUDATES IN FUNDUS IMAGES USING IMPROVED OTSU THRESHOLDING AND SVM

Semantics and image content integration for pulmonary nodule interpretation in thoracic computed tomography

A New Diagnosis Loseless Compression Method for Digital Mammography Based on Multiple Arbitrary Shape ROIs Coding Framework

Detection of Lung Cancer at Early Stage using Neural Network Techniques for Preventing Health Care

Research Article Statistical Segmentation of Regions of Interest on a Mammographic Image

Optimal Planning of Charging Station for Phased Electric Vehicle *

Available online at ScienceDirect. Procedia Computer Science 46 (2015 )

Automatic Labelling and BI-RADS Characterisation of Mammogram Densities

EXAMINATION OF THE DENSITY OF SEMEN AND ANALYSIS OF SPERM CELL MOVEMENT. 1. INTRODUCTION

Proceedings of the 6th WSEAS Int. Conf. on EVOLUTIONARY COMPUTING, Lisbon, Portugal, June 16-18, 2005 (pp )

Survival Rate of Patients of Ovarian Cancer: Rough Set Approach

International Journal of Emerging Technologies in Computational and Applied Sciences (IJETCAS)

Biomarker Selection from Gene Expression Data for Tumour Categorization Using Bat Algorithm

Shape-based Retrieval of Heart Sounds for Disease Similarity Detection Tanveer Syeda-Mahmood, Fei Wang

DETECTION AND CLASSIFICATION OF BRAIN TUMOR USING ML

INITIAL ANALYSIS OF AWS-OBSERVED TEMPERATURE

310 Int'l Conf. Par. and Dist. Proc. Tech. and Appl. PDPTA'16

Prediction of Total Pressure Drop in Stenotic Coronary Arteries with Their Geometric Parameters

Comparative Analysis of Feature Extraction Methods for Optic Disc Detection

Myocardial Motion Analysis of Echocardiography Images using Optical Flow Radial Direction Distribution

Lymphoma Cancer Classification Using Genetic Programming with SNR Features

Using Past Queries for Resource Selection in Distributed Information Retrieval

Investigation of zinc oxide thin film by spectroscopic ellipsometry

Fast Algorithm for Vectorcardiogram and Interbeat Intervals Analysis: Application for Premature Ventricular Contractions Classification

Modeling the Survival of Retrospective Clinical Data from Prostate Cancer Patients in Komfo Anokye Teaching Hospital, Ghana

Parameter Estimates of a Random Regression Test Day Model for First Three Lactation Somatic Cell Scores

An Approach to Discover Dependencies between Service Operations*

Joint Modelling Approaches in diabetes research. Francisco Gude Clinical Epidemiology Unit, Hospital Clínico Universitario de Santiago

Evaluation of the generalized gamma as a tool for treatment planning optimization

Perceptual image quality: Effects of tone characteristics

Algorithms 2009, 2, ; doi: /a OPEN ACCESS

Jurnal Teknologi USING ASSOCIATION RULES TO STUDY PATTERNS OF MEDICINE USE IN THAI ADULT DEPRESSED PATIENTS. Full Paper

Classification of Breast Tumor in Mammogram Images Using Unsupervised Feature Learning

An Improved Time Domain Pitch Detection Algorithm for Pathological Voice

Dr.S.Sumathi 1, Mrs.V.Agalya 2 Mahendra Engineering College, Mahendhirapuri, Mallasamudram

Towards Automated Pose Invariant 3D Dental Biometrics

Recognition of ASL for Human-robot Interaction

Automated and ERP-Based Diagnosis of Attention-Deficit Hyperactivity Disorder in Children

AUTOMATED CHARACTERIZATION OF ESOPHAGEAL AND SEVERELY INJURED VOICES BY MEANS OF ACOUSTIC PARAMETERS

Computing and Using Reputations for Internet Ratings

Physical Model for the Evolution of the Genetic Code

Experimental Study of Dielectric Properties of Human Lung Tissue in Vitro

ARTICLE IN PRESS Biomedical Signal Processing and Control xxx (2011) xxx xxx

ARTICLE IN PRESS. computer methods and programs in biomedicine xxx (2007) xxx xxx. journal homepage:

Association Analysis and Distribution of Chronic Gastritis Syndromes Based on Associated Density

A MIXTURE OF EXPERTS FOR CATARACT DIAGNOSIS IN HOSPITAL SCREENING DATA

Arrhythmia Detection based on Morphological and Time-frequency Features of T-wave in Electrocardiogram ABSTRACT

Research Article Segmentation of Bone with Region Based Active Contour Model in PD Weighted MR Images of Shoulder

Using a Wavelet Representation for Classification of Movement in Bed

IMPROVING THE EFFICIENCY OF BIOMARKER IDENTIFICATION USING BIOLOGICAL KNOWLEDGE

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

Sparse Representation of HCP Grayordinate Data Reveals. Novel Functional Architecture of Cerebral Cortex

JOINT SUB-CLASSIFIERS ONE CLASS CLASSIFICATION MODEL FOR AVIAN INFLUENZA OUTBREAK DETECTION

Insights in Genetics and Genomics

THE NORMAL DISTRIBUTION AND Z-SCORES COMMON CORE ALGEBRA II

Balanced Query Methods for Improving OCR-Based Retrieval

Detection of Cancer Metastasis Using a Novel Macroscopic Hyperspectral Method

A SIMULATION STUDY OF MECHANISM OF POSTFLIGHT ORTHOSTATIC INTOLERANCE

Resampling Methods for the Area Under the ROC Curve

Functional and Molecular Imaging for Radiation Therapy Guidance

Heart Rate Variability Analysis Diagnosing Atrial Fibrillation

A GEOGRAPHICAL AND STATISTICAL ANALYSIS OF LEUKEMIA DEATHS RELATING TO NUCLEAR POWER PLANTS. Whitney Thompson, Sarah McGinnis, Darius McDaniel,

PERFORMANCE EVALUATION OF DIVERSIFIED SVM KERNEL FUNCTIONS FOR BREAST TUMOR EARLY PROGNOSIS

An Introduction to Modern Measurement Theory

Subject-Adaptive Real-Time Sleep Stage Classification Based on Conditional Random Field

Appendix for. Institutions and Behavior: Experimental Evidence on the Effects of Democracy

CLUSTERING is always popular in modern technology

The Limits of Individual Identification from Sample Allele Frequencies: Theory and Statistical Analysis

EVALUATION OF BULK MODULUS AND RING DIAMETER OF SOME TELLURITE GLASS SYSTEMS

Prognosis and Diagnosis of Breast Cancer Using Interactive Dashboard Through Big Data Analytics

ARTICLE IN PRESS Neuropsychologia xxx (2010) xxx xxx

Project title: Mathematical Models of Fish Populations in Marine Reserves

Evaluation of Literature-based Discovery Systems

AN ENHANCED GAGS BASED MTSVSL LEARNING TECHNIQUE FOR CANCER MOLECULAR PATTERN PREDICTION OF CANCER CLASSIFICATION

Reconstruction of gene regulatory network of colon cancer using information theoretic approach

Journal of Engineering Science and Technology Review 11 (2) (2018) Research Article

RENAL FUNCTION AND ACE INHIBITORS IN RENAL ARTERY STENOSISA/adbon et al. 651

Machine Understanding - a new area of research aimed at building thinking/understanding machines

FEATURE EXTRACTION AND SELECTION FOR AUTOMATIC SLEEP STAGING USING EEG

Richard Williams Notre Dame Sociology Meetings of the European Survey Research Association Ljubljana,

A Classification Model for Imbalanced Medical Data based on PCA and Farther Distance based Synthetic Minority Oversampling Technique

Introduction ORIGINAL RESEARCH

Statistical Analysis on Infectious Diseases in Dubai, UAE

Decreased Nailfold Capillary Density in Limited Scleroderma with Pulmonary Hypertension. and a longer disease duration. 3,4

Does reporting heterogeneity bias the measurement of health disparities?

Efficiency Considerations for the Purely Tapered Interference Fit (TIF) Abutments Used in Dental Implants

Submitted for Presentation 94th Annual Meeting of the Transportation Research Board January 11-15, 2015, Washington, D.C.

A-UNIFAC Modeling of Binary and Multicomponent Phase Equilibria of Fatty Esters+Water+Methanol+Glycerol

Transcription:

Hndaw Publshng Corporaton Internatonal Journal of Bomedcal Imagng Volume 2015, Artcle ID 267807, 7 pages http://dx.do.org/10.1155/2015/267807 Research Artcle Statstcal Analyss of Haralck Texture Features to Dscrmnate Lung Abnormaltes Nourhan Zayed and Heba A. Elnemr Computer & Systems Department, Electroncs Research Insttute, Caro 12611, Egypt Correspondence should be addressed to Nourhan Zayed; nourhan@er.sc.eg Receved 23 May 2015; Revsed 10 September 2015; Accepted 15 September 2015 Academc Edtor: Tange Zhuang Copyrght 2015 N. Zayed and H. A. Elnemr. Ths s an open access artcle dstrbuted under the Creatve Commons Attrbuton Lcense, whch permts unrestrcted use, dstrbuton, and reproducton n any medum, provded the orgnal work s properly cted. The Haralck texture features are a well-known mathematcal method to detect the lung abnormaltes and gve the opportunty to the physcan to localze the abnormalty tssue type, ether lung tumor or pulmonary edema. In ths paper, statstcal evaluaton of the dfferent features wll represent the reported performance of the proposed method. Thrty-seven patents CT datasets wth ether lung tumor or pulmonary edema were ncluded n ths study. The CT mages are frst preprocessed for nose reducton and mage enhancement, followed by segmentaton technques to segment the lungs, and fnally Haralck texture features to detect the type of the abnormalty wthn the lungs. In spte of the presence of low contrast and hgh nose n mages, the proposed algorthms ntroduce promsng results n detectng the abnormalty of lungs n most of the patents n comparson wth the normal and suggest that some of the features are sgnfcantly recommended than others. 1. Introducton The lung s an organ that performs a multtude of vtal functons every second of our lves. Ths fact leads to consderng lung abnormaltes, lfe-sustaned dseases that have hgh prorty n detecton, dagnoss, and treatment f possble. Our focus n ths paper wll be on two popular abnormaltes wthn the lung, whch are pulmonary edema and lung tumor. Pulmonary edema (water n the lungs) s caused by flud buldng up n the ar sacs of the lungs [1, 2]. On the other hand, lung cancer/tumor s a dsease where uncontrolled cell growth n tssues of the lung occurred [3]. Computer-aded dagnoss (CAD) schemes for thoracc computed tomography (CT) are wdely used to characterze, quantfy, and detect numerous lung abnormaltes, such as pulmonary edema and lung cancer [4, 5]. An accurate lung segmentaton method s always a crtcal frst step n these CAD schemes and can sgnfcantly mprove the performance level of these schemes. Although manual or semautomatc lung segmentaton methods for CT mages were used n some early CAD schemes [6 10], they are mpractcal for current CAD schemes because multdetector CT (MDCT) scanners can generate hundreds of CT slces for a patent. An automated method for lung segmentaton s needed for MDCT. In addton, the eye dentfcaton/detecton of the abnormalty type (pulmonary edema or tumors) n computed tomography (CT) mages s very dffcult even for the experenced clncans because of ts varable shape along wth low contrast and hgh nose assocated wth t. As the fnal stage of treatng the lung cancer s surgcal removal of the dseased lung, hence t s necessary to dentfy the cancer locaton, whch can be useful before they plan for the surgery. The am of our work s to develop an automated novel textureanalyssbasedmethodforthesegmentatonofthelungs and the detecton of the abnormaltes, whether pulmonary edema or lung tumor. Haralck s features based on the gray level cooccurrence matrx (GLCM) are appled to capture textural patterns n lung mages. The obectve of ths work s the selecton of the most dscrmnatng and fndng out the sgnfcant texture features that can dfferentate between thesetwotypesofabnormaltes,ncomparsontonormal. Haralck features are statstcal features that are computed over the entre mage. These measurements are utlzed to descrbe the overall texture of the mage usng measures

2 Internatonal Journal of Bomedcal Imagng (a) (b) (c) Fgure 1: (a) The lung CT mage; (b) the hstogram equalzed mage; (c) the Wener fltered output mage. such as entropy and sum of varance. Chaddad et al. propose an approach, based on Haralck s features, to detect and classfy colon cancer cells. Ths work amed to select the most dscrmnatng parameters for cancer cells [11]. A study to nvestgate the feasblty of usng Haralck features to dscrmnate between cancer and normal submages wthn a patent s llustrated n [12]. In ths paper, CT mages are frst preprocessed for nose reducton and mage enhancement, followed by segmentaton technques, as the tools to segment the lungs, and fnally Haralck texture features [13 15] are calculated. Statstcal analyss s done to detect the most sgnfcant Haralck features that wll characterze the type of the abnormalty wthn the lungs. Despte the low contrast and hgh nose exstence n the mages, the proposed algorthms ntroduce promsng results n detectng the abnormalty of lungs n most of the patents n comparson wth the normal. 2. Materals and Methods Ths paper presents a new automatc lung cancer detecton system based on Haralck texture features extracted from the slce of DICOM Lung CT mages. The proposed system s accomplshed n four stages: mage preprocessng, lung mage segmentaton, feature extracton, and classfcaton. Statstcalanalysssusedtoobtanthebestfeaturesfor classfcaton to dfferentate between lung cancer patents, ordered edema patents, and control subects. The followng sectons wll descrbe n detal these stages. All mage analyses were acheved wthout any knowledge of patent clncal characterstcs or status. 2.1. Dataset. Patents wth ether a lung cancer tumor or pulmonary edema were encompassed n the study. Ths study ncluded two datasets, the frst dataset referred to the Radology Department at New Elkasr ElAny teachng hosptal,unverstyofcaro.theotherdatasetwasobtaned from The Cancer Imagng Archve (TCIA) sponsored by the SPIE, NCI/NIH, AAPM, and the Unversty of Chcago [16]. The two datasets of 532 CT mages from 37 dfferent patents were ncluded. The mages are 512 512storednDICOM format. For each lung CT mage, we separate the left lung from the rght lung automatcally, and each separated lung s labeled as normal or edema/cancer based on the dataset nformaton. 2.2. Preprocessng. The man goal of preprocessng s to mprove the qualty of an mage as well as make t n a form suted for further processng by human or machne [17]. Ths s accomplshed by enhancng the vsual appearance of an mage besdes removng the rrelevant nose and unwanted parts n the background. The proposed enhancement process, whch s based on combnng flters and nose reducton technques for pre- and postprocessng as well, s carred out applyng hstogram equalzaton (HE) [18 20] followed by Wener flterng [21, 22]. Fgure 1 presents the enhancement n the lung mage contrast attaned by applyng the hstogram equalzaton. However, the obtaned gray scale mage contans noses such as whte nose and salt and pepper nose. Thus, Wener flter s utlzed to remove these noses from the enhanced lung mage. Fgure 1(c) shows the effect of applyng Wener flter on the contrast enhanced lung mage. 2.3. Lung Segmentaton. Lung segmentaton step ams to bascally extract the voxels correspondng to the lung cavty n the axal CT scan slces from the surroundng lung anatomy. The segmentaton technque proposed n [23] s utlzed. Ths technque s based on the fact that there s a large densty dfference between ar-flled lung tssues and surroundng tssues. Furthermore, both lungs are almost lookng lke mrror mages of themselves n a human body. The segmentaton of lung regons s acheved through the followng steps. In thefrststep,thepreprocessedctmagesconvertedntoa bnary mage; a threshold of 128 was selected. Values greater than the threshold are mapped to whte, whle others less than that are marked as black. Consequently, the two lungs are marked and the area around them s cropped out. Second, an eroson morphologcal operaton s employed n order to elmnate any whte pxels wthn the two lungs. Afterward,

Internatonal Journal of Bomedcal Imagng 3 (a) (b) (c) (d) (e) Fgure 2: (a) The threshold mage; (b) the eroded mage; (c) the lung mask mrror; (d) the mask proecton of the correspondng lungs mages; (e) the extracted lungs. the eroded and the orgnal mages are both dvded nto two equal regons. Black pxels for each regon n the eroded mage are counted; the regon wth the largest black area wll be deemed as a lung mask. The attaned lung mask s reflected n the opposte drecton. As a result, rght and left lungmasksareobtaned.thesemasksaremultpledwth the correspondng orgnal mage regons; ths wll proect the lung masks on the orgnal two lungs mages. Fnally, update each black pxel n the obtaned mages by ts orgnal value; other pxels are set to 255. Fgure 2 llustrates the lung extracton process. 2.4. Feature Extracton. Feature extracton s the process of obtanng hgher-level nformaton of an mage such as color, shape, and texture. Texture s a key component of human vsual percepton. Statstcal texture methods analyze the spatal dstrbuton of gray values, by computng local features at each pont n the mage and nferrng a set of statstcs from the dstrbutons of the local features. Haralck et al. ntroduced Gray Level Cooccurrence Matrx (GLCM) and texturefeaturesbackn1973[13].thstechnquehasbeen wdely used n mage analyss applcatons, especally n the bomedcal feld. It conssts of two steps for feature extracton. The GLCM s computed n the frst step, whle the texture featuresbasedontheglcmarecalculatednthesecondstep. GLCM shows how often each gray level occurs at a pxel located at a fxed geometrc poston relatve to each other pxel, as a functon of the gray level [13]. The horzontal drecton 00 wth a range of 1 (nearest neghbor) was used n ths work. The 9 texture descrptons used are presented n (4) to (13), where N g sthenumberofgraylevels,p d s the normalzed symmetrc GLCM of dmenson N g N g, and p d (, ) s the (, )th element of the normalzed GLCM [13]. Contrast (Moment 2 or standard devaton) s a measure of ntensty or gray level varatons between the reference pxel and ts neghbor. Large contrast reflects large ntensty dfferences n GLCM: Contrast = ( ) 2 p d (, ). Homogenety measures how close the dstrbuton of elementsntheglcmstothedagonalofglcm.ashomogenety ncreases, the contrast, typcally, decreases: 1 Homogenety = 1+( ) 2 p d (, ). (2) Entropy s the randomness or the degree of dsorder present nthemage.thevalueofentropysthelargestwhenall elementsofthecooccurrencematrxarethesameandsmall when elements are unequal: Entropy = p d (, ) ln p d (, ). Energy s derved from the Angular Second Moment (ASM). The ASM measures the local unformty of the gray levels. (1) (3)

4 Internatonal Journal of Bomedcal Imagng When pxels are very smlar, the ASM value wll be large. Consder Energy = ASM ASM = p 2 d (, ). (4) Correlaton feature shows the lnear dependency of gray level values n the cooccurrence matrx: Correlaton = p d (, ) ( μ x)( μ y ), (5) σ x σ y where μ x ; μ y and σ x ; σ y are the means and standard devatons and are expressed as μ x = μ y = σ x = σ y = p d (, ) p d (, ) ( μ x ) 2 p d (, ) ( μ y ) 2 p d (, ). The moments are the statstcal expectaton of certan power functons of a random varable and are characterzed as follows. Moment 1 (m 1 )sthemeanwhchstheaverageofpxel values n an mage and t s represented as m 1 = ( ) p d (, ). Moment 2 (m 2 ) s the standard devaton that can be denoted as m 2 = ( ) 2 p d (, ). Moment 3 (m 3 ) measures the degree of asymmetry n the dstrbuton and t s defned as m 3 = ( ) 3 p d (, ). And fnally Moment 4 (m 4 ) measures the relatve peak or flatness of a dstrbuton and s also known as kurtoss: m 4 = ( ) 4 p d (, ). (6) (7) (8) (9) (10) Furthermore, dfference statstcs that are a subset of the cooccurrencematrxarealsoused.thesefeaturesarebased on the dstrbuton of probablty P x y (k) whch s defned as follows: P x y (k) = C d (,), k=0,1,...,n g 1, (11) where C d (, ) s the (, )th element of the GLCM. The most basc dfference statstc texture descrptons are the ASM, mean, and entropy: ASM = (P x y (k)) 2. (12) k When the P x y (k) values are very smlar or close, ASM s small. ASM s large when certan values are hgh and others are low: Mean = kp x y (k). (13) k When P x y (k) values are concentrated near the orgn, mean s small and mean s large when they are far from the orgn: Entropy = P x y (k) log P x y (k). (14) k Entropy s smallest when P x y (k) values are unequal and largest when P x y (k) values are equal. The calculaton of the Haralck texture features usng the prevous equatons for the CT mages volume sequences foreverysegmentedlung(rghtandleft)separatelywas performed. For each partcpant the gray level cooccurrence texture features: contrast, homogenety, entropy, energy, correlaton, and m 1, m 2, m 3,andm 4 accompaned by the dfference statstcal features: ASM, contrast, mean, and entropy were obtaned for each segmented lung (rght and left). 2.5. Statstcal Analyss. For the purpose of random lung assgnment n healthy volunteers, the left lung represented the dseased lung n the same percentage of cases as the patent populaton. For the acute data, two sngle factor analyses of varance (ANOVA) tests were conducted for each Haralck texture feature measurement between affected (ether left or rght) and fellow lung (ether left or rght) for both categores cancer and edema patents. A sngle factor analyss of varance (ANOVA) was conducted as well between patents and controls. Other between-subect sngle factor analyses were conducted to fnd out the sgnfcant Haralck features that could dfferentate cancer from edema. 3. Expermental Results Two datasets of 532 CT mages were ncluded. For each lung CT preprocessed mage, we separate the left lung from the rght lung automatcally as dscussed before n Secton 2.3, and each separated lung s labeled as normal or edema/cancer based on the dataset nformaton. The Haralck texture features measurements for each lung separately are calculated (the gray level cooccurrence texture features: contrast, homogenety, entropy, energy, correlaton, and moments along wth the dfference statstcal features: ASM, mean, and entropy). The mean and the standard devaton of the Haralck texture features measurements calculated as well as the ANOVA results are gven for tumor patents affected lung versus fellow lung n Table 1 and for pulmonary edema patents n Table 2. The ANOVA summary of statstcs for ether pulmonary

Internatonal Journal of Bomedcal Imagng 5 Table 1: ANOVA (1 wthn-subect factor) results for cancer patents Haralck texture features (comparson between AL and FL). AL: affected lung; FL: fellow lung. Feature name AL (average ± SEM) FL (average ± SEM) AL versus FL Homogenety 0.511 ± 0.01 0.517 ± 0.01 F(1, 426) = 22.0 p< 0.000004 Energy 0.372 ± 0.01 0.374 ± 0.01 F(1, 426) = 15.1 p< 0.0001 Correlaton 0.964 ± 0.001 0.965 ± 0.001 F(1, 426) = 6.15 p < 0.013 Contrast 231.98 ± 4.54 231.76 ± 4.54 F(1, 426) = 0.012 p < 0.911 Entropy 8.0 ± 0.19 7.94 ± 0.19 F(1, 426) = 11.8 p< 0.0007 m 1 0.003 ± 0.02 0.007 ± 0.02 F(1, 426) = 0.029 p < 0.88 m 2 231.13 ± 4.54 231.75 ± 4.01 F(1, 426) = 0.012 p < 0.911 m 3 164 ± 190.79 683.99 ± 155.33 F(1, 426) = 2.65 p < 0.09 m 4 1784467 ± 83311 1654941 ± 56455 F(1, 426) = 6.25 p < 0.19 Dff ASM 0.226963389 ± 0.006 0.229353096 ± 0.005 F(1, 426) = 3.18 p < 0.06 Dff Mean 6.195 ± 0.08 6.28 ± 0.09 F(1, 426) = 2.16 p < 0.12 Dff Entropy 3.159 ± 0.03 3.55 ± 0.03 F(1, 426) = 2.32 p < 0.10 Table 2: ANOVA (1 wthn-subect factor) results for edema patents Haralck texture features (comparson between AL and FL). AL: affected lung; FL: fellow lung. Feature name AL (average ± SEM) FL (average ± SEM) AL versus FL Homogenety 0.64 ± 0.013 0.60 ± 0.020 F(1, 105) = 2.16 p < 0.15 Energy 0.428 ± 0.01 0.429.01 ± 0.01 F(1, 105) = 0.029 p < 0.87 Correlaton 0.006 ± 0.001 0.008 ± 0.001 F(1, 105) = 2.32 p < 0.141 Contrast 177.07 ± 5.89 188.58 ± 4.26 F(1, 105) = 15.1 p< 0.0002 Entropy 2.10 ± 0.04 2.19 ± 0.067 F(1, 105) = 1.28 p < 0.269 m 1 0.52 ± 0.03 0.47 ± 0.02 F(1, 105) =41.8 p< 0.000001 m 2 199.975 ± 9.658 218.583 ± 10.085 F(1, 105) = 2.20 p < 0.152 m 3 5219 ± 1436 7539 ± 885 F(1, 105) =41.8 p< 0.000001 m 4 2854294 ± 208886 2382237 ± 263250 F(1, 105) = 2.12 p < 0.158 Dff ASM 0.377 ± 0.01 0.288 ± 0.08 F(1, 105) =4.56 p < 0.043 Dff Mean 4.07 ± 0.4379 4.89 ± 0.48478 F(1, 105) = 7.87 p < 0.01 Dff Entropy 2.96 ± 0.05 3.29 ± 0.05 F(1, 105) = 4.73 p < 0.039 edema or tumor patents versus normal s gven n Table 3. The sgnfcant Haralck texture features that can dfferentate between pulmonary edema and tumor are found n Table 4. From Table 1, we can conclude that Haralck texture features measurements (homogenety, energy, correlaton, and entropy) of the affected cancer lung were sgnfcantly dfferent than that of the fellow lung. The homogenety, energy, and correlaton were sgnfcantly less than those of thenormalfellowlung.whleentropyofthecancerouslung s approachng beng sgnfcantly more than that of the fellow lung, Moment 3 and the dfference statstcal feature ASM (dff ASM) texture feature measurement of the cancerous lung s approachng beng sgnfcantly less than that of the normal lung. Table 2 showed that Haralck texture features measurements (homogenety, entropy, and moments calculated from the cooccurrence matrx as well as mean and ASM computed from the dfference statstcs) of the pulmonary edema affected lung were also sgnfcantly dfferent than those of the control subect lung; moreover contrast and entropy computed from the dfference statstcs were sgnfcantly more than those of the fellow lung. Consderng Tables 1 and 2, we can conclude from Table 3 that the homogenety, energy, entropy, m 3, m 4,dffASM, dff mean, and dff entropy are good bomarkers to sgnfcantly dfferentate between dseased and normal lungs wthout any dsease specfcaton. On the other hand, the results llustrated n Tables 1, 2, and 4 show that entropy and the entropy calculated from the dfference statstcs would be a good canddate to sgnfcantly dfferentate between pulmonary edema and cancer. 4. Concluson and Dscusson The texture features analyses are well known approaches to quantfy and express the heterogenety that may not be apprecated by clncal naked eyes, and t was presented before as good magng bomarkers to dfferentate between dseases. InthspaperanevaluatonoftheHaralcktexturefeatures s done n order to dentfy the most sgnfcant features that can be used n order to detect and dfferentate abnormaltes wthn the lungs for cancer and edema versus normal. Our results ndcate that entropy determned by gray level cooccurrence matrx and ASM s sgnfcantly dfferent n

6 Internatonal Journal of Bomedcal Imagng Table 3: ANOVA (1 wthn-subect factor) results summary of statstcs p value for patents (ether edema or cancer) Haralck texture features versus normal controls. Feature name Dseased versus normal controls (p value) Feature name Dseased versus normal controls (p value) Homogenety p < 0.00002 m 2 p < 0.229 Energy p < 0.0006 m 3 p < 0.0002 Correlaton p < 0.485 m 4 p < 0.002 Contrast p < 0.229 Dff ASM p < 0.000001 Entropy p < 0.0004 Dff Mean p < 0.000005 m 1 p < 0.0007 Dff Entropy p < 0.0004 Table 4: ANOVA (1 between-subect factor) results summary of statstcs p value for patents Haralck texture features cancer versus edema patents. Feature name Cancer versus edema patents (p value) Feature name Cancer versus edema patents (p value) Homogenety p < 0.0002 m 2 p < 0.69 Energy p < 0.065 m 3 p < 0.01 Correlaton p < 0.179 m 4 p < 0.89 Contrast p < 0.699 ASM p < 0.73 Entropy p < 0.017 Mean p < 0.032 m 1 p < 0.0004 Entropy p < 0.007 edema patents versus normal whle t s not n cancer patents versus normal. Snce the entropy s the degree of randomness or the degree of dsorder n the mage, and the angular second moment represents the unformty n the mage, ths may be nterpreted as the cancer dsease causng a localzed heterogenety n the dseased specfed area n the lung whle the edema causes heterogeneous dsorder n the whole lung mage. Hgh entropy values calculated mples that the elevated level of dsorder and dsorganzaton occurred due to the edema dseased lung versus the cancer dseased lung. The energy feature that s derved from the angular second moment measures and representng the local unformty of the gray levels s a good bomarker to dfferentate between cancer and edema dseases. From Table 2, contrast s a good bomarker for the pulmonary edema dsease and ths agrees wth the texture feature meanng whch means hgh contrast values for heavy texture changes. Gray level cooccurrence matrx textural propertes such as homogenety, correlaton, mean, and moments are good sgnfcant bomarkers for dseased lung versus normal ones n general wthout any specfcaton for the dsease type. Our results agree wth other artcles ndcatng that textural analyss has the potental to develop nto a valuable clncal tool that mproves the dagnoss, tumor stagng, and therapy assessment. Whle our results are promsng, there s stll further work that can be done n the detectng of the abnormalty wthn the lungs to detect the type of that abnormalty whether t wll be a lung cancer or edema. A prelmnary nvestgaton has been done usng statstcal analyss to dentfy the most useful texturefeaturesthatcanbefedtoanyclassfcatontechnque later. Ths statstcal analyss s done usng ANOVA. After selectng these features we can feed them for better localzaton and classfcaton as further work. Conflct of Interests The authors declare that there s no conflct of nterests regardng the publcaton of ths paper. Acknowledgments The authors acknowledge the SPIE, the NCI, the AAPM, and The Unversty of Chcago for provdng publc access to the lung cancer dataset. References [1] L.B.WareandM.A.Matthay, Acutepulmonaryedema, The New England Journal of Medcne, vol.353,no.26,pp.2788 2796, 2005. [2] J. Šedý, J. Zcha, J. Kuneš, P. Jendelová, and E. Syková, Mechansms of neurogenc pulmonary edema development, Physologcal Research,vol.57,no.4,pp.499 506,2008. [3] R. S. Herbst, J. V. Heymach, and S. M. Lppman, Lung cancer, The New England Journal of Medcne,vol.359,no.13,pp.1367 1380, 2008. [4] N. Hollngs and P. Shaw, Dagnostc magng of lung cancer, European Respratory Journal,vol.19,no.4,pp.722 742,2002. [5] T. Mankandan and N. Bharath, Lobar fssure extracton n sotropc CT lung mages an applcaton to cancer dentfcaton, Internatonal Journal of Computer Applcatons,vol.33,no. 6, pp. 17 21, 2011. [6] J. Wang, F. L, and Q. L, Automated segmentaton of lungs wth severe ntersttal lung dsease n CT, Medcal Physcs, vol. 36, no.10,pp.4592 4599,2009. [7] U. Baĝc,M.Bray,J.Caban,J.Yao,andD.J.Mollura, Computerasssted detecton of nfectous lung dseases: a revew, Computerzed Medcal Imagng and Graphcs, vol.36,no.1,pp.72 84, 2012.

Internatonal Journal of Bomedcal Imagng 7 [8] D. Sharma and G. Jndal, Computer aded dagnoss system for detecton of lung cancer n CT scan mages, Internatonal Journal of Computer and Electrcal Engneerng,vol.3,no.5,pp. 714 718, 2011. [9] J. Wang, F. L, K. Do, and Q. L, Computerzed detecton of dffuse lung dsease n MDCT: the usefulness of statstcal texture features, Physcs n Medcne and Bology, vol.54,no. 22, pp. 6881 6899, 2009. [10] J. Wang, F. L, K. Do, and Q. L, A novel scheme for detecton of dffuse lung dsease n MDCT by use of statstcal texture features, n Medcal Imagng 2009: Computer-Aded Dagnoss, 726039, vol. 7260 of Proceedngs of SPIE, Lake Buena Vsta, Fla, USA, February 2009. [11] A. Chaddad, C. Tanougast, A. Dandache, and A. Bourdane, Extracton of haralck features from segmented texture multspectral bo-mages for detecton of colon cancer cells, n Proceedngs of the 1st Internatonal Conference on Informatcs and Computatonal Intellgence (ICI 11), pp.55 59,Bandung, Indonesa, December 2011. [12] B.D.Fleet,J.Yan,D.B.Knoester,M.Yao,J.R.DellerJr.,andE. D. Goodman, Breast cancer detecton usng haralck features of mages reconstructed from ultra wdeband mcrowave scans, n Clncal Image-Based Procedures. Translatonal Research n Medcal Imagng,vol.8680ofLecture Notes n Computer Scence, pp. 9 16, Sprnger Internatonal, Cham, Swtzerland, 2014. [13] R. M. Haralck, K. Shanmugam, and I. Dnsten, Textural features for mage classfcaton, IEEE Transactons on Systems, Man and Cybernetcs,vol.3,no.6,pp.610 621,1973. [14] R. M. Haralck, Statstcal and structural approaches to texture, Proceedngs of the IEEE, vol. 67, no. 5, pp. 786 804, 1979. [15] J. M. H. du Buf, M. Kardan, and M. Spann, Texture feature performance for mage segmentaton, Pattern Recognton,vol. 23, no. 3-4, pp. 291 309, 1990. [16] S. G. Armato III, SPIE-AAPM-NCI Lung Nodule Classfcaton Challenge Dataset, The Cancer Imagng Archve, 2015. [17] R. C. Gonzalez and R. E. Woods, Dgtal Image Processng, Addson-Wesley, 1992. [18] C. Wang and Z. Ye, Brghtness preservng hstogram equalzaton wth maxmum entropy: a varatonal perspectve, IEEE Transactons on Consumer Electroncs, vol.51,no.4,pp.1326 1334, 2005. [19] I. Jafar and H. Yng, Image contrast enhancement by constraned varatonal hstogram equalzaton, n Proceedngs of the IEEE Internatonal Conference on Electro/Informaton Technology (EIT 07), pp. 120 125, IEEE, Chcago, Ill, USA, May 2007. [20] S.-D. Chen and A. R. Raml, Mnmum mean brghtness error b-hstogram equalzaton n contrast enhancement, IEEE Transactons on Consumer Electroncs, vol.49,no.4,pp.1310 1319, 2003. [21] D. Sharma and G. Jndal, Computer aded dagnoss system for detecton of lung cancer n CT scan mages, Internatonal Journal of Computer and Electrcal Engneerng,vol.3,no.5,2011. [22] S. V. Vasegh, Advanced Dgtal Sgnal Processng and Nose Reducton, John Wley & Sons, New York, NY, USA, 2nd edton, 2000. [23] H. A. Elnemr, Statstcal analyss of law s mask texture features for cancer and water lung detecton, Internatonal Journal of Computer Scence Issues,vol.10,no.6,2013.