Recent Deforestation Pattern Changes (2000–2017) in the Central Carpathians: A Gray-Level Co-Occurrence Matrix and Fractal Analysis Approach

The paper explores the distribution of tree cover and deforested areas in the Central Carpathians in the central-east part of Romania, in the context of the anthropogenic forest disturbances and sustainable forest management. The study aims to evaluate the spatiotemporal changes in deforested areas due to human pressure in the Carpathian Mountains, a sensitive biodiverse European ecosystem. We used an analysis of satellite imagery with Landsat-7 Enhanced Thematic Mapper Plus (Landsat-7 ETM+) from the University of Maryland (UMD) Global Forest Change (GFC) dataset. The workflow started with the determination of tree cover and deforested areas from 2000–2017, with an overall accuracy of 97%. For the monitoring of forest dynamics, a Gray-Level Co-occurrence Matrix analysis (Entropy) and fractal analysis (Fractal Fragmentation-Compaction Index and Tug-of-War Lacunarity) were utilized. The increased fragmentation of tree cover (annually 2000–2017) was demonstrated by the highest values of the Fractal Fragmentation-Compaction Index, a measure of the degree of disorder (Entropy) and heterogeneity (Lacunarity). The principal outcome of the research reveals the dynamics of disturbance of tree cover and deforested areas expressed by the textural and fractal analysis. The results obtained can be used in the future development and adaptation of forestry management policies to ensure sustainable management of exploited forest areas.


Introduction
Forest ecosystems around the world have experienced a great deal of change [1]. The development of alternative forest management solutions was able to incorporate multiple stakeholder preferences [2]; the importance of ecological forest values primarily influenced the forest management attitude among private forest owners in Sweden [3]. On the other hand, in Indonesia, land cover changes are completed by participatory mapping of stakeholders [4]. The forest dynamics, caused by deforestation, determine a severe geomorphological impact, related to soil loss and water erosion, due to long-term logging activities and more clear-cut areas [5].
After the collapse of socialism, changes in forest dynamics were recorded, caused by increases in the human population, requiring more arable land for agriculture [6,7], illegal logging or legal wood harvesting [8,9], and even fluctuant changes in forestry legislation [10].
Deforestation causes significant perturbations in the integrity of forest ecosystems. The fragmentation of forests due to deforestation causes flash floods [11,12], landslides [13][14][15], soil degradation [5,16], loss of biodiversity [17,18] and increases in CO 2 emissions [19][20][21]. Deforestation may be the result of logging and harvesting practices or can be a permanent conversion to other land uses (i.e., fires, agriculture, and development of human settlements) [22]. According to FAO [23], deforestation represents the loss of forest cover and implies transformation into another land use caused by human or natural perturbations. Deforestation reduces the essential values of a forest's role as a climatic regulator [24,25], hydrological function [26,27], social and aesthetic values [28,29], and ecological values [3]. Deforestation has become a significant perturbation in forest stability for the Carpathian Mountains, one of the mountains that preserve virgin forests [30], protected areas [8], and endemic species of flora and fauna [17,31]. The Carpathian Mountains represent the largest mountain range in Central Europe and the second largest in Europe with a length of about 910 km within Romanian borders [32,33].
The changes in forest policy created the context for the tremendous ownership changes from state to private-owned in the mountain areas, as in Romanian Carpathians [8,34]. This situation is recorded in many ex-communist countries as Poland [35] or Slovakia [36,37].
The effects of deforestation and forest fragmentation were the subject of numerous research papers for the Carpathian Mountains. Most of the studies mapped the regional forest dynamics but with less focus on particular changes in the European Carpathian environment [38], quantifying spatial metrics such box-counting fractal dimensions from Amazonian forests [39], or added distinct indices in the fractal analysis (i.e., Fractal Fragmentation-Compaction Index-FFI) [40]. These studies focus on important methodological approaches, but there are needed improvements at the local scale using some fractal indices: Gray-Level Co-occurrence Matrix (GLCM), Fractal Fragmentation-Compaction Index (FFI) or Tug-of-War Lacunarity (Λ T−o−W ).
Most recent changes in forest fragmentation have not been studied using fractal analysis at a local scale or using GLCM and fractal analysis. The purpose of this paper is to quantify the amount of tree cover and deforested areas and its degree of fragmentation in the Central Carpathians in Romania, with an innovative methodological use of the spatial approaches that provide complementary information to traditional per-pixel deforestation mapping.

Study Area
The Central Carpathians represent the study area of the paper; a subunit of the Eastern Carpathians, bounded by the North Carpathian Mountains in the north and by the Curvature Carpathian Mountains in the south (Figure 1). With a diversity of landscapes due to the geologic mosaic, this mountainous area covers three sub-regions: the volcanic rocks in the western part (Călimani, Gurghiu and Harghita Mountains), the crystalline schist rocks (Bistrit , ei, Ceahlău and Tarcău Mountains), (Giurgeu, Ciuc and Nemira Mountains) and the flysch rocks (Stânis , oarei, Gos , manu and Berzunt , i Mountains). These sub-regions are separated mainly by three lower areas: Giurgeu, Ciuc and Comănes , ti. These mountains are covered by tree species as spruce (Picea abies), fir (Abies alba), oak (Quercus robur), and ash (Fraxinus excelsior) [30] and overlap with the four most deforested counties of Romania (Suceava, Neamt , , Bacău and Harghita).

Forest Imagery and Pre-Processing
The analysis of deforestation patterns in Central Carpathians followed three main steps as shown in the flowchart ( Figure 2). The first is related to the Landsat-7 Enhanced Thematic Mapper Plus (Landsat-7 ETM+) image data pre-processing and Geographic Information System Analysis (GIS analysis) for mapping the changes in deforestation patterns (Table 1), while the second part is the

Forest Imagery and Pre-Processing
The analysis of deforestation patterns in Central Carpathians followed three main steps as shown in the flowchart (Figure 2). The first is related to the Landsat-7 Enhanced Thematic Mapper Plus (Landsat-7 ETM+) image data pre-processing and Geographic Information System Analysis (GIS analysis) for mapping the changes in deforestation patterns (Table 1), while the second part is the statistical validation (classification accuracy assessment) of them. The last part was based on the GLCM and fractal analysis of forest area dynamics.  The process starts with the collection of satellite images from 2000-2017 (a spatial resolution of 30 m) from Global Forest Change (GFC) dataset from the Department of Geographical Sciences, University of Maryland (UMD) [41]. These images were used to extract the deforested areas (annually 2001 to 2017) and, tree cover areas. The primary data source for our study is the GFC, and the process for obtaining the numerical information followed an entire algorithm.
The first step consisted of downloading the Raster data. Having a different projection (World Geodetic System (WGS) 84-European Petroleum Standards Group (EPSG): 4326) than the other spatial data (Dealul Piscului 1970 (Stereo 70)-European Petroleum Standards Group (EPSG): 31700), used in the study, a projection transformation step followed.
Next, a mask was applied to delineate the study area. To extract the numerical information, the raster image was converted from raster to vector, each pixel being transformed into a point by classes, but keeping all their properties. So, 18 classes were obtained for the deforestation image (class 0was ignored because it corresponded to the entire surface of the study area-no interesting information, class 1 for deforestation of the year 2001, class 2 for deforestation of the year 2002, till class 17 for deforestation of 2017). From the tree cover image, of the same data source, 100 classes were obtained, each class corresponding to the degree of forest cover in each pixel (class 0 for pixels  [41]. These images were used to extract the deforested areas (annually 2001 to 2017) and, tree cover areas. The primary data source for our study is the GFC, and the process for obtaining the numerical information followed an entire algorithm.
The first step consisted of downloading the Raster data. Having a different projection (World Geodetic System (WGS) 84-European Petroleum Standards Group (EPSG): 4326) than the other spatial data (Dealul Piscului 1970 (Stereo 70)-European Petroleum Standards Group (EPSG): 31700), used in the study, a projection transformation step followed.
Next, a mask was applied to delineate the study area. To extract the numerical information, the raster image was converted from raster to vector, each pixel being transformed into a point by classes, but keeping all their properties. So, 18 classes were obtained for the deforestation image (class 0-was ignored because it corresponded to the entire surface of the study area-no interesting information, class 1 for deforestation of the year 2001, class 2 for deforestation of the year 2002, till class 17 for deforestation of 2017). From the tree cover image, of the same data source, 100 classes were obtained, each class corresponding to the degree of forest cover in each pixel (class 0 for pixels with 0% tree cover, class 100 to 100% tree cover). In all cases, the spatial dimension (resolution) of each pixel was taken by each point. For the numerical data extracted from the deforestation image, the number of the points for each class was multiplied with the surface of the point/spatial resolution of pixel (24.514 × 24.514 = 600.97 m 2 ), but for the tree cover, the entire surface (600.97 m 2 ) was giving only for pixels/points in class 100 (with 100% of tree cover). For the other classes, the surface decreasing proportionally, until class 0.
In order to have accurate information, and to compare with other economic indicators at a territorial administrative unit level, a spatial join was done. Then, all this information was compiled into Excel tables, which were used to continue the analysis.

Classification Accuracy Assessment
The GFC dataset for Central Carpathians results for 2001-2017 has been compared to a Google Earth image (February 2017) classified into deforested and non-deforested for validation purposes. Two hundred random ground reference test points have been selected and analysed, and classification error matrices have been calculated to estimate the accuracy assessment of the obtained GFC dataset (Figure 3a-c). In the selection of ground reference test points, some bias in the classification result may arise. It turns out that the agreement between the GFC dataset and the Google Earth Image mappings is very high (see Table 2) with an overall classification accuracy of 97%. In the GFC dataset, 96% of non-deforested pixels where correctly classified (Producer's Accuracy and User's Accuracy = 98% and commission error = 2.5% and omission error = 1.54%), whereas for the deforested pixels, the Producer's Accuracy = 60% and User's Accuracy = 50%. The confusion between deforested and non-deforested pixels is easily explained due to ground control errors or preprocessing from the GFC dataset (i.e., some pixels are classified as deforested ones even though there are observed compact forests; also, the pastures with clumps of forest are classified as deforested areas). with 0% tree cover, class 100 to 100% tree cover). In all cases, the spatial dimension (resolution) of each pixel was taken by each point. For the numerical data extracted from the deforestation image, the number of the points for each class was multiplied with the surface of the point/spatial resolution of pixel (24.514 × 24.514 = 600.97 m²), but for the tree cover, the entire surface (600.97 m²) was giving only for pixels/points in class 100 (with 100% of tree cover). For the other classes, the surface decreasing proportionally, until class 0.
In order to have accurate information, and to compare with other economic indicators at a territorial administrative unit level, a spatial join was done. Then, all this information was compiled into Excel tables, which were used to continue the analysis.

Classification Accuracy Assessment
The GFC dataset for Central Carpathians results for 2001-2017 has been compared to a Google Earth image (February 2017) classified into deforested and non-deforested for validation purposes. Two hundred random ground reference test points have been selected and analysed, and classification error matrices have been calculated to estimate the accuracy assessment of the obtained GFC dataset (Figure 3a-c). In the selection of ground reference test points, some bias in the classification result may arise. It turns out that the agreement between the GFC dataset and the Google Earth Image mappings is very high (see Table 2) with an overall classification accuracy of 97%. In the GFC dataset, 96% of non-deforested pixels where correctly classified (Producer's Accuracy and User's Accuracy = 98% and commission error = 2.5% and omission error = 1.54%), whereas for the deforested pixels, the Producer's Accuracy = 60% and User's Accuracy = 50%. The confusion between deforested and non-deforested pixels is easily explained due to ground control errors or preprocessing from the GFC dataset (i.e., some pixels are classified as deforested ones even though there are observed compact forests; also, the pastures with clumps of forest are classified as

Preprocessing of GFC for GLCM and Fractal Analysis
To perform the GLCM and fractal analysis, it was necessary to go through several processes for making the maps for analysis. The first step was to estimate the tree cover and loss areas for Central Carpathians in grayscale. The second requirement was to export images to the same proportional scale in TIFF format to preserve their properties.
Subsequently, images were transformed into binary images using ImageJ 1.52 [43] and analysed using a macro for calculating the area and the percentage of forest pixels. The binarization was made by converting foreground pixels (which in turn corresponds to tree cover and loss) of grey and white tones, with intensities between 1 and 255, only in white pixels with the intensity of 255. The black background pixels (non-tree cover and non-loss) with 0 intensity remained the same.
The resolution of the analyzed images was 2716 × 2716 pixels. Pixels representing the tree cover and loss areas are automatically extracted from the GFC dataset on a grey scale. These images were used for Entropy GLCM analysis, but for the Λ T−o−W and FFI, the same images were binarized, all pixels indicating tree cover and loss becoming white.
GFC dataset provides information about the 2000 tree cover and loss areas for 2001-2017. Using the function Image Calculator operator (from ImageJ) the cumulative loss and tree cover data were obtained. Cumulative loss (the summed values for each year) was obtained using the Add function in Image Calculator. i.e., Cumulative loss for 2017 was obtained by adding the loss of each year from 2001 to 2017. The tree cover for the years 2001-2017 was obtained by using the Difference function in the Image Calculator. i.e., Tree cover 2017 was obtained by using difference between Tree cover 2000 and Cumulative loss 2017.

GLCM and Fractal Analysis
GLCM analysis, also known as the grey-level spatial dependence matrix, is a statistical method of examining texture that considers the spatial relationship of pixels. The GLCM functions characterise the composition of an image by calculating how often pairs of the pixel of specific values and in a specified spatial relationship occur in an image, creating a GLCM, and then extracting statistical measures from this matrix. They provide information about shape, i.e., the spatial relationships of pixels in an image. The GLCM measures how often a pixel of grey-level 8-bits images (grayscale intensity or tone) value i occurs either horizontally, vertically, or diagonally to adjacent pixels with the value j.
To determine the Entropy, a significant method for non-texturally uniform images and for small values of elements, a GLCM analysis was used.
Complex textures tend to have high entropy values. The entropy is measured according to Equation (1) [44].
where p(i,j) are co-occurrence probabilities and i and j are coordinates of the co-occurrence matrix. Entropy measures the degree of disorder inside patches of tree cover and loss based on the relationship between pixels with different degrees of forest coverage. Complex texture tends to have high entropy values. The entropy expresses the degree of disturbance of the tree cover canopy of the forest, and it is in strong relation with the degree of fragmentation expressed by FFI. GLCM Entropy is complemented by FFI and Λ T−o−W , two binary analyzes, where all foreground pixels become white, measure the degree of fragmentation of the patches, some relative to each other.
For the analysis of fractal indices as FFI and Λ T−o−W , were used FFI plugin [45] and Frac2D plugin [46] in the IQM 3.5 software [47].
These indices are relevant because they quantify the spatial pattern of tree cover and loss areas by analysing the degree of fragmentation (FFI), and heterogeneity (Λ T−o−W ). It is expected that as the loss increases, fragmentation and heterogeneity of forests will increase and at the same time a process of homogenization and compacting of cumulative losses will occur.
Fractal fragmentation is very useful in estimating fragmentation or compaction of fractal and natural objects that do not follow the classical geometry. FFI is calculated according to the Equation (2) and can be interpreted as a compaction index [40]: where D B−C A is the box-counting fractal dimension of the summed-up areas; D B−C P is the box-counting fractal dimension of the summed-up perimeters; ε is the side length of the box; N(ε) is the number of contiguous and non-overlapping boxes required to cover the area of the object; and N'(ε) is the number of contiguous and non-overlapping boxes required to cover just the perimeter of the object [48,49]. The tree cover patches and loss are very small and highly fragmented and appear as point-like objects in the image; the degree of fragmentation is maximum, FFI = 0, according to Equation (2), in the situation where D B−C A = D B−C P. As the analyzed areas are more compact, the FFI value will increase to 1, and as they are more fragmented, irregular, the FFI will be closer to 0. The maximum compaction, when the areas are perfectly geometric, have an FFI = 1. However, self-similar objects, such as forests, with an identical fractal dimension may differ significantly in their textural appearance [50,51]. Therefore, the use of fractal fragmentation only is not useful for discriminating against objects, while the fractal fragmentation dimension quantifies the way space is occupied, and the lacunarity completes the fractal dimension with its ability to quantify how space is filled. Moreover, lacunarity discriminates the spatial distribution of gaps in texture at multiple scales, and is not sensitive to the edges of the images. In this paper, the Tug-of-War algorithm lacunarity [46] was calculated based on the equation: with N(r) is the number of boxes, Z 2 is the second moment for each width as the median of s 2 values, each is the mean of s 1 squares of the counter values. s 1 and s 2 are two random variables for each width that indicate the accuracy and confidence.

Validation of GLCM and Fractal Analysis Indices
To check the normality of the fractal and non-fractal indicators used in the present study, a Shapiro-Wilk statistical test has been applied, instead of the graphical one, due to the small amount of data (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017). In such cases, this statistical test is more easily understood instead of graphical ones. The data displayed in Table 3 indicate a normal distribution (W values-over 0.75). This hypothesis is also sustained by the p-values, which are above 0.05 in most cases. The p-values close to 0.05 are in direct accordance with the degree of compactness of tree cover or deforested areas.

The Analysis of Deforested Areas from Central Carpathians
In the Central Carpathians, the general decrease of the forest areas during 2000-2017 was observed, however with considerable inter-annual variability. In total, 100.613 ha was deforested during the entire period of analysis, corresponding to a reduction in forest area by 11%. (Figure 4).  In Figure 5a, the central aspect relates to the fact that the tree cover areas are diminishing, and the deforested areas are increasing with some annual variabilities. Also, the highest values of deforested areas were observed in 2007 (14.003 ha.) and 2012 (9.325 ha.).
An important aspect is related to the Annual Average Growth Deforestation Rate (AAGDR) (Figure 5b)    In Figure 5a, the central aspect relates to the fact that the tree cover areas are diminishing, and the deforested areas are increasing with some annual variabilities. Also, the highest values of deforested areas were observed in 2007 (14.003 ha.) and 2012 (9.325 ha.). The AAGDR values present the high rate of growth of deforestation because the deforested areas are reduced from year to year, mainly in 2007 (229.14%) and in 2016 (112.25%) Until 2004, the deforestation is shown to be dispersed in small patches, but from 2005 and onwards clustering of deforestation can be observed. This aspect coincides with an overall diminishing of clusters for the tree cover areas.

Correlation between Deforested Areas and Logging Activities
An important disturbance agent of forest dynamics is represented by the logging activities that are in strong relation with the high deforestation rates. To have a more quantitative causal link between these two components, we have correlated the indicators of economic activities from the Classification of the Activities from the National Economy (NACE code 0220-logging activities) with the deforestation rates, registered at each Territorial Administrative Unit, from Central Carpathians. The AAGDR values present the high rate of growth of deforestation because the deforested areas are reduced from year to year, mainly in 2007 (229.14%) and in 2016 (112.25%) Until 2004, the deforestation is shown to be dispersed in small patches, but from 2005 and onwards clustering of deforestation can be observed. This aspect coincides with an overall diminishing of clusters for the tree cover areas.

Correlation between Deforested Areas and Logging Activities
An important disturbance agent of forest dynamics is represented by the logging activities that are in strong relation with the high deforestation rates. To have a more quantitative causal link between these two components, we have correlated the indicators of economic activities from the Classification of the Activities from the National Economy (NACE code 0220-logging activities) with the deforestation rates, registered at each Territorial Administrative Unit, from Central Carpathians. Figure 6 (a-d) show that for all the four economic indicators, the deforestation rates have no relation with the dimension of logging activities, being distributed on both the two axes. Many territorial administrative units have low concentrations of companies and employees, and low values of profit and turnover with small logging activities. Only some separate territorial administrative-units appear with high values for the four economic indicators. It is those ones where the palette of the economic activities is reduced to forestry activities (the monospecialized ones, from the economic point of view). are in strong relation with the high deforestation rates. To have a more quantitative causal link between these two components, we have correlated the indicators of economic activities from the Classification of the Activities from the National Economy (NACE code 0220-logging activities) with the deforestation rates, registered at each Territorial Administrative Unit, from Central Carpathians. Figure 6 (a-d) show that for all the four economic indicators, the deforestation rates have no relation with the dimension of logging activities, being distributed on both the two axes. Many territorial

The Analysis of GLCM Entropy
The Entropy of tree cover areas was relatively high and steadily increasing from 2.

The Analysis of GLCM Entropy
The Entropy of tree cover areas was relatively high and steadily increasing from 2.

The Analysis of GLCM Entropy
The Entropy of tree cover areas was relatively high and steadily increasing from 2.40 in 2000 to 2.51 in 2017 (change of 0.11). The main reason is a very different type of deforestation practices. However, the entropy values of the deforested areas were relatively small and constant over the

The Analysed Fractal Indices
The fractal indicators examined were FFI and Λ T−o−W . The FFI values were reduced from 0. 22 (2000) to 0.21 (2017) indicating a continuous growth of fragmentation of tree cover areas, with a reduction or even elimination of clusters of tree cover areas. The decrease of FFI, which indicates the increase in fractal fragmentation, is generated by the fragmentation of the patches, and the increase in their number by the detachment of the compact surfaces following the deforestation. The

Deforestation changes in Central Carpathians
Results indicated an increase in deforested areas Central Carpathians from Romania during recent decades. By using satellite images from GFC, the information describing the situation of tree cover and deforested areas from a subdivision of Romanian Carpathians during the period 2000-2017 is provided. The satellite analysis of tree cover and deforested areas was complemented by GLCM and fractal analysis describing differentiations of forest disturbance, providing a quantification of changes in the heterogeneity and fragmentation of forest areas.
A general increment trend characterizes deforested areas in Romanian Carpathians due to high rates of illegal logging [8]. The significant changes in forestry legislation reveal that this growth in deforested areas resulted from the restitution of forest to former owners the Law 247/2005 [52] and privation of forests, especially in the mountain areas [8,10,34,53]. Among various natural and anthropic perturbations in forest areas from Central Carpathians, the logging activities represent an important factor of disturbance, determining the growth of logging activities [54].

Use of GLCM and Fractal Analysis for Quantification the Deforestation Changes
We have evaluated the forest and deforested areas [41] using the GLCM and fractal analysis, determining the entropy with FFI and Λ T-o-W. This research demonstrates the usefulness of GLCM and fractal analysis for the measurement of deforestation patterns [55]. The findings of this study confirmed the hypothesis that fractal analysis adds interesting ways of measuring the degree of fragmentation of tree cover areas caused by the deforestation. Here, GLCM and fractal analysis are used to study the characteristic of deforested and tree cover areas in Central Carpathians.
In the GLCM analysis of tree cover and deforested areas, the entropy analysis is used to describe the spatial relationship between pixels [56]. For the fractal analysis, the FFI and Λ T-o-W methods were performed to define the degree of fragmentation of tree cover areas determined by forest disturbances. In previous studies, fractal analysis of tree cover and deforested areas were calculated by using the FFI [40,57,58], Fixed Grid 2D Lacunarity (FG2DL) [59]. Until now, fractal analysis was used in urban agglomeration, urban growth [60][61][62], or green infrastructure models of cities [63].
The main findings are discussed in the following aspects: 1.
Fractal methods may provide valuable complementary information to currently available methodologies in the field of forestry research; 2.
A continuous decrease of tree cover areas has occurred during the period of analysis, however with considerable inter-annual variability; 3.
Results indicate that, as the loss areas increased, forest fragmentation and heterogeneity also increased. The process of homogenization and compaction of cumulative loss has also been confirmed; 4.
Differences between the fractal and GLCM indices arise from the type of image analyzed and from the information extracted. For the Entropy, the images are grey-scale 8-bits, and we obtain the clutter of spatial pixel distribution in grayscale within patches of tree cover and loss. Thus, the analysis is at the level of forest cover as the forest looks dense. For Λ T−o−W and FFI, the images are binary, and we extract information about how these patches are spatially distributed, regardless of how the forest looks or how loss has occurred within these patches, resulting in anti-parallel developments.
Spatial patterns changes can be distinguished between the analysis of deforested and tree cover areas. The period with more intense deforestation is characterized by a more compact and homogenous way, but disordered at the level of patch forest cover, and vice versa. As the deforestation increases, the cumulative deforested areas increase gradually, with a decrease in the heterogeneity of spatial distribution of the deforested patches and an increase in tree cover disorder around them. The newly deforested areas are made as new patches, alternating with agglutinations around patches previously deforested. Regarding the spatial patterns of tree cover areas, as the deforested areas increase, there is a tendency to increase fractal fragmentation of forests; this is associated with an increase in patch-level forest clutter. Λ T−o−W of tree cover shows tight yearly oscillations related to the degree of heterogeneity of deforested areas.

Prospects About Sustainable Forest Management
The current paper explains the importance of using the entropy, FFI and Λ T−o−W for the better understanding of tree cover areas and its immediate effects as deforestation. The method, based on GLCM and fractal analysis, exemplifies the advantage of making known the effects of forest fragmentation. The fragmentation affects the stability of forest ecosystems and determines the decrease of biodiversity and rare species of fauna and flora [17,18]. The analysis of deforestation by fractal analysis provides complementary information to the already known methods of deforestation mapping.

Conclusions
Here, an analysis of Landsat-7 ETM+ satellite images at 30 m resolution in Central Carpathians is provided, offering quantitative information describing the spatial and temporal dynamics of deforestation and tree cover areas during the period of analysis (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017). Based on the fractal analysis that can expertly analyse irregular spatial structures, new information about deforestation-as described by textural uniformity, compactness and chaotic distribution of the forest-disturbance processes were obtained. The analysis of fractal indicators (entropy, FFI and Λ T−o−W ) describes and quantifies the textural uniformity, compactness and chaotic distribution of the disturbance processes of deforestation. Overall, the Central Carpathians is a region with a highly deforested mountainous area. It is concluded that fractal analysis of deforested areas is an effective tool in quantifying the degree of homogeneity or spatial heterogeneity of deforested areas. Analysis of the degree of textural disorder of deforestation in tree cover areas provides essential information to identify efficient methods for managing legal as well as illegal logging. Such quantitative metrics can also assist in identifying the impact of forest disturbances on biodiversity and local economies and may act as guidance for policy makers from deforested territorial administrative units from Central Carpathians.