Gray-level Co-occurrence Matrix Analysis Research Articles

Age-related changes in ultrasound-assessed muscle composition and postural stability

While the simultaneous degradation of muscle composition and postural stability in aging are independently highly investigated due to their association with fall risk, the interplay between the two has received little attention. Thus, the purpose of this study is to explore how age-related changes in muscle composition relate to postural stability. To that aim, we collected posturography measures and ultrasound images of the dominant Vastus Lateralis and Biceps Brachii from 32 young (18–35 year old) and 34 older (65–85 year old) participants. Muscle properties were quantified with echo-intensity and texture-based metrics derived from gray-level co-occurrence matrix analysis, and postural stability with the variability of the center of pressure during bipedal stance tasks. Ultrasound parameters revealed that young muscle possessed lower echo-intensity and higher homogeneity compared to the elderly. Echo-intensity and muscle thickness, and several texture-based parameters possessed outstanding young versus older classification performance. A canonical correlation analysis demonstrated a significant relationship between ultrasound and postural measures only within the young group (r = 0.53, p < 0.002), where those with ‘better’ muscle composition displayed larger postural sways. Our results indicate that, in older individuals, postural stability and muscle composition, two common fall risk factors, are unrelated. In view of this decoupling, both may contribute independently to fall risk. Furthermore, our data support the view that texture-based parameters provide a robust alternative to echo-intensity in providing markers of muscle composition.

Full-Automatic High-Efficiency Mueller Matrix Microscopy Imaging for Tissue Microarray Inspection.

This paper proposes a full-automatic high-efficiency Mueller matrix microscopic imaging (MMMI) system based on the tissue microarray (TMA) for cancer inspection for the first time. By performing a polar decomposition on the sample's Mueller matrix (MM) obtained by a transmissive MMMI system we established, the linear phase retardance equivalent waveplate fast-axis azimuth and the linear phase retardance are obtained for distinguishing the cancerous tissues from the normal ones based on the differences in their polarization characteristics, where three analyses methods including statistical analysis, the gray-level co-occurrence matrix analysis (GLCM) and the Tamura image processing method (TIPM) are used. Previous MMMI medical diagnostics typically utilized discrete slices for inspection under a high-magnification objective (20×-50×) with a small field of view, while we use the TMA under a low-magnification objective (5×) with a large field of view. Experimental results indicate that MMMI based on TMA can effectively analyze the pathological variations in biological tissues, inspect cancerous cervical tissues, and thus contribute to the diagnosis of postoperative cancer biopsies. Such an inspection method, using a large number of samples within a TMA, is beneficial for obtaining consistent findings and good reproducibility.

Revealing GLCM Metric Variations across a Plant Disease Dataset: A Comprehensive Examination and Future Prospects for Enhanced Deep Learning Applications

This study highlights the intricate relationship between Gray-Level Co-occurrence Matrix (GLCM) metrics and machine learning model performance in the context of plant disease identification. It emphasizes the importance of rigorous dataset evaluation and selection protocols to ensure reliable and generalizable classification outcomes. Through a comprehensive examination of publicly available plant disease datasets, focusing on their performance as measured by GLCM metrics, this research identified dataset_2 (D2), a database of leaf images, as the top performer across all GLCM analyses. These datasets were then utilized to train the DarkNet19 deep learning model, with D2 exhibiting superior performance in both GLCM analysis and DarkNet19 training (achieving about 91% testing accuracy) according to performance metrics such as accuracy, precision, recall, and F1-score. The datasets other than dataset_1 and 2 exhibited significantly low classification performance, particularly in supporting GLCM analysis. The findings underscore the need for transparency and rigor in dataset selection, particularly given the abundance of similar datasets in the literature and the growing trend of utilizing deep learning methods in future scientific research.

Anthelmintic efficacy of Holarrhena pubescens against Raillietina spp. of domestic fowl through ultrastructural, histochemical, biochemical and GLCM analysis.

Globally, traditional knowledge systems are a powerhouse of information which can revolutionise the world, if decoded accurately and logically. Plant-based ethno-traditional and folklore curatives/medicines has a firm basis in the psyche of the common masses of West Bengal and Holarrhena pubescens is a representative example of it. This article communication on depicting the anthelmintic efficacy of ethanolic extract and Ethyl acetate fraction of the stem bark of Holarrhena pubescens against the cestode Raillietina spp. through efficacy studies, ultra-structural observations, histochemical and biochemical analysis on some tegumental enzymes i.e., Acid Phosphatase (AcPase), Alkaline Phosphatase (AlkPase), Adenosine Triphosphatase (ATPase) and 5'-Nucleotidase (5'-Nu) along with Gray Level Co-occurrence Matrix (GLCM) analysis of histochemical study. Praziquantel was used as the reference drug. Investigations revealed 10mg/ml dosage of crude extract was the most efficacious dose and amongst the fractions the ethyl acetate fraction showed the most anthelmintic property. Ultrastructural studies through Scanning Electron Microscope (SEM) and Transmission Electron Microscope (TEM) clearly depicted the damage in head, sucker, proglottids, proximal and distal cytoplasm (DC), microtriches (MT), basal lamina (BL), nuclear membrane (NM), and, nucleolus (NL) in the treated worms. Histochemical studies revealed decrease in staining intensity for all the tegumental enzymes in the treated worms compared to control. The GLCM analysis strongly supported the result of histochemical studies. Biochemical studies revealed marked reduction in enzyme activity in the treated worms with maximum reduction in the activity of 5'- Nu (77.8%) followed by ATPase (63.17%).

Identification of EIMD Level Differences Between Long- and Short Head of Biceps Brachii Using Echo Intensity and GLCM Texture Features

ABSTRACT Purpose: This study aimed to compare the time-course changes of exercise-induced muscle damage (EIMD) levels in the long head of biceps brachii (LHB) and short head of the biceps brachii (SHB) using echo intensity (EI) and to determine the efficiency of the gray level co-occurrence matrix (GLCM) texture parameters. Methods: The participants performed 30 maximal eccentric contractions of the elbow flexor. Along with muscle damage indicators, including circumference, range of motion, muscle soreness, and maximal voluntary isometric contraction (MVIC), the EI and GLCM texture features of the LHB and SHB was also assessed using B-mode ultrasonography. All measurements were assessed pre- and immediately post-exercise and after 24, 48, 72, and 96 h. Results: The muscle damage indicators indicated significant changes after the eccentric contractions (p < 0.01 for circumference, range of motion, muscle soreness, and MVIC). The EI of LHB significantly increased following the contractions (p < 0.01), but that of SHB did not (p > 0.05). In contrast, for the GLCM texture parameters, there were significant changes in the SHB (p < 0.01 for homogeneity, energy, and entropy). Conclusion: Thus, this study demonstrated that EIMD severity is different between LHB and SHB even within the same muscle. In the GLCM features, the time course of SHB after eccentric contraction revealed different patterns compared with those of LHB. Therefore, even if there are no changes in EI within a target muscle following muscle contractions, new information on muscle quality can be obtained through GLCM analysis.

Echo intensity and gray-level co-occurrence matrix analysis of soft tissue grafting biomaterials and dental implants: an in vitro ultrasonographic pilot study.

To characterize different allogeneic and xenogeneic soft tissue graft substitutes and to assess their echo intensity and grayscale texture-related outcomes by using high-frequency ultrasonography (HFUS). Ten samples from each of the following biomaterials were scanned using HFUS: bilayered collagen matrix (CM), cross-linked collagen matrix (CCM), multilayered cross-linked collagen matrix (MCCM), human-derived acellular dermal matrix (HADM), porcine-derived acellular dermal matrix (PADM), collagen tape dressing (C) and dental implants (IMPs). The obtained images were then imported in a commercially available software for grayscale analysis. First-order grayscale outcomes included mean echo intensity (EI), standard deviation, skewness, and kurtosis, while second-order grayscale outcomes comprised entropy, contrast, correlation, energy and homogeneity derive from the gray-level co-occurrence matrix analysis. Descriptive statistics were performed for visualization of results, and one-way analysis of variance with Bonferroni post-hoc tests were performed to relative assessments of the biomaterials. The statistical analysis revealed a statistically significant difference among the groups for EI (p < .001), with the group C showing the lowest EI, and the IMP group presenting with the greatest EI values. All groups showed significantly higher EI when compared with C (p < .001). No significant differences were observed for energy, and correlation, while a statistically significant difference among the groups was found in terms of entropy (p < 0.01), contrast (p < .001) and homogeneity (p < .001). IMP exhibited the highest contrast, that was significantly higher than C, HADM, PADM, CCM and CM. HFUS grayscale analysis can be applied to characterize the structure of different biomaterials and holds potential for translation to in-vivo assessment following soft tissue grafting-related procedures.

Diagnosis of Laryngopharyngeal Reflux Disease Based on Gray and Texture Changes of Laryngoscopic Images

ObjectiveThis study aimed to compare the changing trends of gray and texture values of laryngoscopic images in patients with laryngopharyngeal reflux (LPR) and non-LPR. MethodsA total of 3428 laryngoscopic images were selected and divided into two groups, non-LPR and LPR groups based on the reflux symptom index. Gray histogram and gray-level co-occurrence matrix (GLCM) were used to quantify gray and texture features, and the model was trained based on these features. The total laryngoscopic images dataset was proportionally split into two parts including the training set and the test set according to the ratio of 7:3. Four different machine learning algorithms, including decision tree, naive Bayes, linear regression, and K-nearest neighbors, were applied to classify non-LPR or LPR laryngoscopic images. ResultsThe results showed that different classification algorithms are used to classify laryngoscopic image dataset and promising classification accuracy are obtained. Specifically, the accuracy of K-nearest neighbors was 83.38% for the gray histogram-only classification, that of linear regression was 88.63% for the GLCM-only classification, and that of the decision tree was 98.01% for the combined gray histogram and GLCM analysis. ConclusionGray histogram and GLCM analysis of the laryngoscopic images may be used as auxiliary tools to detect laryngopharyngeal mucosal damage in patients with LPR. Measurement of gray and texture feature values is an objective and convenient method, which may serve as a reference baseline for clinicians and have potential clinical usefulness.

Gray-Level Co-occurrence Matrix Analysis of Nuclear Textural Patterns in Laryngeal Squamous Cell Carcinoma: Focus on Artificial Intelligence Methods.

Gray-level co-occurrence matrix (GLCM) and discrete wavelet transform (DWT) analyses are two contemporary computational methods that can identify discrete changes in cell and tissue textural features. Previous research has indicated that these methods may be applicable in the pathology for identification and classification of various types of cancers. In this study, we present findings that squamous epithelial cells in laryngeal carcinoma, which appear morphologically intact during conventional pathohistological evaluation, have distinct nuclear GLCM and DWT features. The average values of nuclear GLCM indicators of these cells, such as angular second moment, inverse difference moment, and textural contrast, substantially differ when compared to those in noncancerous tissue. In this work, we also propose machine learning models based on random forests and support vector machine that can be successfully trained to separate the cells using GLCM and DWT quantifiers as input data. We show that, based on a limited cell sample, these models have relatively good classification accuracy and discriminatory power, which makes them suitable candidates for future development of AI-based sensors potentially applicable in laryngeal carcinoma diagnostic protocols.

Quantitative Evaluation of the Effectiveness of Chemical Peelings in Reducing Acne Lesions Based on Gray-Level Co-Occurrence Matrix (GLCM).

PurposeAcne vulgaris is a chronic, inflammatory disease accompanied by lesions affecting the structure of the skin. Chemical peels are one of the methods of reducing acne vulgaris. There is still a lack of quantitative methods of assessing impact of cosmetic procedure on the skin. Skin condition depends on skin texture characterization; therefore, the analysis that provides data about the textures can be helpful in assessing the effectiveness of cosmetic treatments.Patients and MethodsThe study involved 24 volunteers with acne lesions. Each participant underwent 4 treatments using chemical peels at two-week intervals. Before, during and after procedure clinical photography were made. To assess effectiveness of chemical peeling in acne lesion reduction, we were used gray-level co-occurrence matrix (GLCM) analysis. Qualitative assessment of acne severity was made by 12 experts in dermatology.ResultsAfter a series of treatments, the GLCM contrast value decreased in each area of the face, and the GLCM homogeneity value increased, which means that the number of acne lesions was reduced. Expert assessment according to the IGA scale confirms the effectiveness of therapy with both salicylic and glycolic acid and pyruvic acid.ConclusionThe results of this study prove that GLCM analysis is a useful tool for assessing the effectiveness of chemical peel treatments. It can also be used for quantitative assessment of skin texture.

Gray-Level Co-occurrence Matrix Analysis for the Detection of Discrete, Ethanol-Induced, Structural Changes in Cell Nuclei: An Artificial Intelligence Approach.

Gray-level co-occurrence matrix (GLCM) analysis is a contemporary and innovative computational method for the assessment of textural patterns, applicable in almost any area of microscopy. The aim of our research was to perform the GLCM analysis of cell nuclei in Saccharomyces cerevisiae yeast cells after the induction of sublethal cell damage with ethyl alcohol, and to evaluate the performance of various machine learning (ML) models regarding their ability to separate damaged from intact cells. For each cell nucleus, five GLCM parameters were calculated: angular second moment, inverse difference moment, GLCM contrast, GLCM correlation, and textural variance. Based on the obtained GLCM data, we applied three ML approaches: neural network, random trees, and binomial logistic regression. Statistically significant differences in GLCM features were observed between treated and untreated cells. The multilayer perceptron neural network had the highest classification accuracy. The model also showed a relatively high level of sensitivity and specificity, as well as an excellent discriminatory power in the separation of treated from untreated cells. To the best of our knowledge, this is the first study to demonstrate that it is possible to create a relatively sensitive GLCM-based ML model for the detection of alcohol-induced damage in Saccharomyces cerevisiae cell nuclei.

Using Google Earth Engine to classify unique forest and agroforest classes using a mix of Sentinel 2a spectral data and topographical features: a Sri Lanka case study

Global land cover classifications may lead to the loss of important local and national nuances such as forest and agroforestry classes. These classes are important to local contexts because they contribute to sustainable land management systems. This paper demonstrates the application of Sentinel-2A satellite images, elevation data, and the Google Earth Engine platform to generate more detailed, specialist land cover classification for forestry classes important in Sri Lanka deriving ten spectral, 16 textural, and three topographical features from the input datasets. The random forest classification model discriminates vegetation types as forest, forest plantations, shrub, grassland, home garden, and cultivation with an overall accuracy of 94% and kappa value of 0.91. Results indicate the elevation feature contributes the most to discriminate forest and agroforestry classes, and red band (664.6 nm) textural metrics derived from grey-level co-occurrence matrix analysis are more useful for separating the home garden from other land cover classes.

Association between Chromatin Structural Organization of Peripheral Blood Neutrophils and Self-Perceived Mental Stress: Gray-Level Co-occurrence Matrix Analysis.

Methods based on the evaluation of textural patterns in microscopy, such as the “gray-level co-occurrence matrix” (GLCM) analysis are modern and innovative computer and mathematical algorithms that can be used to quantify subtle structural changes in cells and their organelles. Potential application of GLCM method in the fields of psychophysiology and psychiatry to this date has not been systematically investigated. The main objective of our study was to test the existence and strength of the association between chromatin structural organization of peripheral blood neutrophils and levels of self-perceived mental stress. The research was done on a sample of 100 healthy student athletes, and the Depression, Anxiety, and Stress Scales (DASS-21) were used for the estimation of psychological distress. Chromatin textural homogeneity and uniformity were negatively correlated (p < 0.01) with mental distress and had relatively good discriminatory power in differentiating participants with normal and elevated stress levels. As an addition, we propose the creation of a machine learning model based on binomial logistic regression that uses these and other GLCM features to predict stress elevation. To the best of our knowledge, these results are one of the first to establish the link between neutrophil chromatin structural organization quantified by the GLCM method and indicators of normal psychological functioning.

Differentiation of Amyloid Plaques Between Alzheimer's Disease and Non-Alzheimer's Disease Individuals Based on Gray-Level Co-occurrence Matrix Texture Analysis.

Amyloid plaques, one of the main hallmarks of Alzheimer's disease (AD), are classified into diffuse (associated with cognitive impairment) and dense-core types (a common finding in brains of people without Alzheimer's disease (non-AD) and without impaired cognitive function) based on their morphology. We tried to determine the usability of gray-level co-occurrence matrix (GLCM) texture parameters of homogeneity and heterogeneity for the differentiation of amyloid plaque images obtained from AD and non-AD individuals. Images of amyloid-β (Aβ) immunostained brain tissue samples were obtained from the Aging, Dementia and Traumatic Brain Injury Project. A total of 1,039 plaques were isolated from different brain regions of 69 AD and non-AD individuals and used for further GLCM analysis. Images of Aβ stained plaques show higher values of heterogeneity parameters and lower values of homogeneity parameters in AD patients, and vice versa in non-AD patients. Additionally, GLCM analysis shows differences in Aβ plaque texture between different brain regions in non-AD patients and correlates with variables that characterize patient's dementia status. The present study shows that GLCM texture analysis is an efficient method to discriminate between different types of amyloid plaques based on their morphology and thus can prove as a valuable tool in the neuropathological investigation of dementia.

Effects of Iron Oxide Nanoparticles on Structural Organization of Hepatocyte Chromatin: Gray Level Co-occurrence Matrix Analysis.

Gray level co-occurrence matrix (GLCM) analysis is a contemporary and innovative computer-based algorithm that can be used for the quantification of subtle changes in a cellular structure. In this work, we use this method for the detection of discrete alterations in hepatocyte chromatin distribution after in vivo exposure to iron oxide nanoparticles (IONPs). The study was performed on 40 male, healthy C57BL/6 mice divided into four groups: three experimental groups that received different doses of IONPs and 1 control group. We describe the dose-dependent reduction of chromatin textural uniformity measured as GLCM angular second moment. Similar changes were detected for chromatin textural uniformity expressed as the value of inverse difference moment. To the best of our knowledge, this is the first study to investigate the impact of iron-based nanomaterials on hepatocyte GLCM parameters. Also, this is the first study to apply discrete wavelet transform analysis, as a supplementary method to GLCM, for the assessment of hepatocyte chromatin structure in these conditions. The results may present the useful basis for future research on the application of GLCM and DWT methods in pathology and other medical research areas.

Optimization of knee cartilage texture analysis of quantitative MRI T2 maps

Optimization of knee cartilage texture analysis of quantitative MRI T2 maps

A novel approach to thermographic images analysis of equine thoracolumbar region: the effect of effort and rider\u2019s body weight on structural image complexity

BackgroundThe horses’ backs are particularly exposed to overload and injuries due to direct contact with the saddle and the influence of e.g. the rider’s body weight. The maximal load for a horse’s back during riding has been suggested not to exceed 20% of the horses’ body weight. The common prevalence of back problems in riding horses prompted the popularization of thermography of the thoracolumbar region. However, the analysis methods of thermographic images used so far do not distinguish loaded horses with body weight varying between 10 and 20%.ResultsThe superficial body temperature (SBT) of the thoracolumbar region of the horse’s back was imaged using a non-contact thermographic camera before and after riding under riders with LBW (low body weight, 10%) and HBW (high body weight, 15%). Images were analyzed using six methods: five recent SBT analyses and the novel approach based on Gray Level Co-Occurrence Matrix (GLCM) and Gray Level Run Length Matrix (GLRLM). Temperatures of the horse’s thoracolumbar region were higher (p < 0.0001) after then before the training, and did not differ depending on the rider’s body weight (p > 0.05), regardless of used SBT analysis method. Effort-dependent differences (p < 0.05) were noted for six features of GLCM and GLRLM analysis. The values of selected GLCM and GLRLM features also differed (p < 0.05) between the LBW and HBW groups.ConclusionThe GLCM and GLRLM analyses allowed the differentiation of horses subjected to a load of 10 and 15% of their body weights while horseback riding in contrast to the previously used SBT analysis methods. Both types of analyzing methods allow to differentiation thermal images obtained before and after riding. The textural analysis, including selected features of GLCM or GLRLM, seems to be promising tools in considering the quantitative assessment of thermographic images of horses’ thoracolumbar region.

Second-order grey-scale texture analysis of pleural ultrasound images to differentiate acute respiratory distress syndrome and cardiogenic pulmonary edema

Discriminating acute respiratory distress syndrome (ARDS) from acute cardiogenic pulmonary edema (CPE) may be challenging in critically ill patients. Aim of this study was to investigate if gray-level co-occurrence matrix (GLCM) analysis of lung ultrasound (LUS) images can differentiate ARDS from CPE. The study population consisted of critically ill patients admitted to intensive care unit (ICU) with acute respiratory failure and submitted to LUS and extravascular lung water monitoring, and of a healthy control group (HCG). A digital analysis of pleural line and subpleural space, based on the GLCM with second order statistical texture analysis, was tested. We prospectively evaluated 47 subjects: 16 with a clinical diagnosis of CPE, 8 of ARDS, and 23 healthy subjects. By comparing ARDS and CPE patients’ subgroups with HCG, the one-way ANOVA models found a statistical significance in 9 out of 11 GLCM textural features. Post-hoc pairwise comparisons found statistical significance within each matrix feature for ARDS vs. CPE and CPE vs. HCG (P ≤ 0.001 for all). For ARDS vs. HCG a statistical significance occurred only in two matrix features (correlation: P = 0.005; homogeneity: P = 0.048). The quantitative method proposed has shown high diagnostic accuracy in differentiating normal lung from ARDS or CPE, and good diagnostic accuracy in differentiating CPE and ARDS. Gray-level co-occurrence matrix analysis of LUS images has the potential to aid pulmonary edemas differential diagnosis.

The Fractal and GLCM Textural Parameters of Chromatin May Be Potential Biomarkers of Papillary Thyroid Carcinoma in Hashimoto's Thyroiditis Specimens.

Occasionally, Hashimoto's thyroiditis (HT) and papillary thyroid carcinoma (PTC) share similar nuclear features. The current study aims to quantify the differences between the investigated specimens of HT-associated PTC versus the HT alone, to reduce the subjective experience of an observer, by the use of fractal parameters as well as gray-level co-occurrence matrix (GLCM) textural parameters. We have analyzed 250 segmented nuclei per group (nn = 25 per patient and np = 10 patients per group) using the ImageJ software (NIH, Bethesda, MD, USA) as well as an in-house written code for the GLCM analysis. The mean values of parameters were calculated for each patient. The results demonstrated that the malignant cells from the HT-associated PTC specimens showed lower chromatin fractal dimension (p = 0.0321) and higher lacunarity (p = 0.0038) compared with the corresponding cells from the HT specimens. Also, there was a statistically significant difference between the investigated specimens, in the contrast, correlation, angular second moment, and homogeneity, of the GLCM corresponding to the visual texture of follicular cell chromatin. The differences in chromatin fractal and GLCM parameters could be integrated with other diagnostic methods for the improved evaluation of distinctive features of the HT-associated PTC versus the HT in cytology and surgical pathology specimens.

Gray-Level Co-Occurrence Matrix Analysis of Granule Neurons of the Hippocampal Dentate Gyrus Following Cortical Injury.

Traumatic brain injury (TBI) is a main cause of death and disabilities in young adults. Although learning and memory impairments are a major clinical manifestation of TBI, the consequences of TBI on the hippocampus are still not well understood. In particular, how lesions to the sensorimotor cortex damage the hippocampus, to which it is not directly connected, is still elusive. Here, we study the effects of sensorimotor cortex ablation (SCA) on the hippocampal dentate gyrus, by applying a highly sensitive gray-level co-occurrence matrix (GLCM) analysis. Using GLCM analysis of granule neurons, we discovered, in our TBI paradigm, subtle changes in granule cell (GC) morphology, including textual uniformity, contrast, and variance, which is not detected by conventional microscopy. We conclude that sensorimotor cortex trauma leads to specific changes in the hippocampus that advance our understanding of the cellular underpinnings of cognitive impairments in TBI. Moreover, we identified GLCM analysis as a highly sensitive method to detect subtle changes in the GC layers that is expected to significantly improve further studies investigating the impact of TBI on hippocampal neuropathology.

Application of Fractal and Gray-Level Co-Occurrence Matrix Indices to Assess the Forest Dynamics in the Curvature Carpathians—Romania

The mountain ecosystems face significant damage from deforestation and environmental forest changes. We investigated the evolution of tree types of cover areas, deforested areas and total deforested areas from Curvature Carpathians using Gray-Level Co-occurrence Matrix and fractal analysis. The forest dynamics mapping was one of the main objectives of this study and it was carried out using multiple fractal and GLCM indices. We approached the analysis of satellite forest images by calculation of four fractal indices such as Pyramid dimension, Cube Counting Dimension, Fractal Fragmentation-Compaction Index and Tug-of-War lacunarity. We also calculated fractal dimension because it is an index of complexity comparing how the detail in a pattern changes with the scale at which it is measured. Fractal dimension is useful for estimation of irregularity or roughness of fractal and natural objects that do not conform to Euclidian geometry. While the fractal dimension quantifies how much space is occupied, the Tug-of-War lacunarity complements fractal dimension with its ability to quantify how space is occupied. Analysis was further supplemented by the Gray-Level Co-occurrence Matrix analysis because it quantifies spatial probability distributions of gray level values between pixel pairs within an image. The calculated Gray-Level Co-occurrence Matrix features included Angular Second Moment, Contrast, Correlation, Inverse Difference Moment and Entropy. Such comprehensive analysis has the advantage of combining fractal analysis that extracts quantitative information about the morphological complexity of the image with the spatial distribution of the gray pixel intensities as calculated by the co-occurrence features provided by Gray-Level Co-occurrence Matrix. Evolution of deforested areas, expansion of agricultural land and the increased demand for quality timber have affected the forests ecosystems and, the regional sustainable development of local communities.