Comparison of Landscape Metrics for Three Different Level Land Cover/Land Use Maps

This research aims to investigate how different landscape metrics are affected by the enhancement of the thematic classes in land cover/land use (LC/LU) maps. For this aim, three different LC/LU maps based on three different levels of CORINE (Coordination of Information on The Environment) nomenclature were created for the selected study area using GEOBIA (Geographic Object Based Image Analysis) techniques. First, second and third level LC/LU maps of the study area have five, thirteen and twenty-seven hierarchical thematic classes, respectively. High-resolution Spot 7 images with 1.5 m spatial resolution were used as the main Earth Observation data to create LC/LU maps. Additional geospatial data from open sources (OpenStreetMap and Wikimapia) were also integrated to the classification in order to identify some of the 2nd and 3rd level LC/LU classes. Classification procedure was initially conducted for Level 3 classes in which we developed decision trees to be used in object-based classification. Afterwards, Level 3 classes were merged to create Level 2 LC/LU map and then Level 2 classes were merged to create the Level 1 LC/LU map according to CORINE nomenclature. The accuracy of Level 1, Level 2, Level 3 maps are calculated as; 93.50%, 89.00%, 85.50% respectively. At the last stage, several landscape metrics such as Number of Patch (NP), Edge Density (ED), Largest Patch Index (LPI), Euclidean Nearest Neighbor Distance (ENN), Splitting Index (SPLIT) and Aggregation Index (AI) metrics and others were calculated for different level LC/LU maps and landscape metrics values were compared to analyze the impact of changing thematic details on landscape metrics. Our results show that, increasing the thematic detail allows landscape characteristics to be defined more precisely and ensure comprehensive assessment of cause and effect relationships between classes.


Introduction
The Earth' surface changes constantly, not only by its natural dynamics, but also due to human activities.These changes, whether they are minor or major, local or global, affect the ecology and the human beings living on that ecosystem.In order to better understand the cause and effect relationships of different changes on the Earth' surface, such as deforestation, urbanization and other processes, some of the most crucial information is up to date spatio-temporal distribution of land cover and land use.Land cover is defined as the physical cover of Earth and land use describes how that land is used by humans.Together, they provide a baseline for many ecological and societal studies such as environmental models, weather and climate studies, hydrological and landscape planning studies, etc. [1][2][3].Therefore, mapping land cover/land use (LC/LU) characteristics of the Earth's surface becomes an important topic.There are different standard LC/LU systems developed by different agencies for various scales and needs, such as "National Land Cover Database (NLCD) of United States", "United Nations Land cover maps (LCCS) of Food and Agricultural Organization (FAO)" and CORINE (Coordination of Information on The Environment) project of European Union (EU) [4].
CORINE is a monitoring program that aims to create standardized seamless LC/LU maps throughout the EU and its candidate countries in order to aid environmental policies [5,6].It is updated every six years and uses satellite imagery as its main data source.CORINE land cover nomenclature provides a conceptual framework for understanding the different types of LC/LU classes and is widely used in related studies [7,8].It is based on the visual interpretation and on screen digitization of the LC/LU classes according to the defined nomenclature.Five main classes defined for the first level of this nomenclature are; Artificial Surfaces (1), Agricultural Areas (2), Forest and Semi-Natural Areas (3), Wetlands (4) and Water Bodies (5).These five classes become thematically more detailed at the second level with 15 classes and the most detailed thematic classes are defined at the third level with a total number of 44 [9].
Producing the LC/LU maps with on screen digitization could be very time consuming when applied to large areas.Thus, researchers have been exploring ways to include some level of automation into this process by classifying the satellite images [10,11].However, since the detailed LC/LU systems such as third level of CORINE have complicated class definitions, automatic or semi-automatic classification of these classes becomes a challenging task.Spatial resolution of the satellite images affects the level of detail in classification and accuracy of the final LC/LU map specifically for complex landscapes [11,12].High and very high resolution satellite images provide significant advantages, especially in urban mapping and analyzing spatial/temporal changes of urban areas, due to their enhanced ability to represent heterogeneous areas and the level of information they provide [13,14].However, their classification becomes more intricate as their resolution improves.
Geographic Object Based Image Analysis (GEOBIA) is a relatively recent classification approach compared to traditional pixel-based techniques.Pixel-based techniques have some limitations since they only use spectral and textural properties of the pixels [14][15][16].Compared to pixel-based classification methods, GEOBIA usually produce higher accuracy values [17][18][19].In the first process step of GEOBIA method, segmentation is applied to produce the image objects that are formed by groups of pixels according to a certain criteria set and then the classification is performed on these objects [18,20].In order to obtain better thematic results from the classification, different features could be added to spectral information such as shape, pattern, tone, texture, context, size, connectivity, proximity of pixels and others.[21].Geometry of a specific area can be determined by using shape and size information and the pattern can be detected with tone and texture information to improve classification results [21,22].GEOBIA takes advantage of all these spectral, spatial and contextual information by incorporating them into segmentation and classification process [23].GEOBIA has proven to be quite successful in creating LC/LU maps of urban areas [13,16,24] and it is quite successful with high and very high resolution satellite images.It provides effective results by producing finer details required for the urban studies, such as creating parcel borders and analyzing adjacent parcels in urban areas [25].Moreover, GEOBIA is preferred for its ability to easily integrate satellite data and different geospatial information obtained from other geographic data sources such as OpenStreetMap, European Environmental Agency Geodata, thematic maps and Wikimapia.
The LC/LU maps are used in various research areas related to dynamic monitoring of the environment and human interactions with the environment.Since a significant portion of human population are in urban environments, mapping of urban areas and their surroundings with LC/LU systems is one of the most important research topics on this regard.LC/LU maps are also useful to analyze the patterns in a landscape.The landscape metrics, which are the indices developed in order to characterize geometric and spatial properties of map patterns, are used effectively for that purpose [26].These metrics enable quantitative and objective analysis of the different landscape types.Thus, they are used to test future scenarios of environmental policies and to monitor environmental goals determined by international conventions and agreements [27].
Landscape is defined by Forman and Gordon as "a heterogeneous land area composed of a cluster of interacting ecosystems that is repeated in similar form throughout" [28].There are several spatial processes that transform the land by affecting landscape in different ways and cause changes in spatial patterns and ecological processes [29].Landscape pattern analysis, which studies landscape components and their spatial patterns, is one of the most effective methods defining these spatial processes [30].In order to analyze the interaction between spatial patterns and ecological processes, first, spatial patterns need to be defined [28,31].Landscape metrics are useful tools to measure arrangement of landscape components both in time and space and they are used to characterize spatial patterns [32].More than one metric is needed to define the landscape pattern.The metric group should describe the pattern variety seen throughout the landscape, but should be minimized in use, especially in indexes that are highly related to each other [31].
Previous studies have revealed that selecting the correct approach or combination of different approaches is an important task in landscape investigations.In recent years, Gradient Surface Model (GSM) has arisen as an alternative to Patch Mosaic Model (PMM) and it is becoming widespread.However, this model also has its shortcomings, which are similar to the unaddressed issues of PMM, such as the ecological relevancy of metrics.Therefore, applications of the GSM are also limited [33].Integration of Gradient surfaces into the ecological analyses is difficult.Additionally, many of the surface metrics might not have patch analogs and might not be ecologically relevant [34].Considering the above mentioned restrictions of GSM, metrics to be used in this study are chosen according to Patch Mosaic Model approach.Besides, PMM also allows for comparison according to CORINE LU/LC levels.
Landscape metrics can be defined at three levels which are namely; patch level, class level and landscape level.Patch level metrics are calculated for individual patches which represent discrete areas of similar characteristics.Class level metrics are calculated from all patches of a particular type, which, in the context of this study, are LC/LU classes.Landscape level metrics are the combination of all patch and class types in a given area [26].There are many landscape metrics such as Largest Patch Index (LPI), Euclidean Nearest Neighbor Distance (ENN), Percentage of Landscape (PLAND), Splitting Index (SPLIT) [35][36][37][38] that can be used in patch, class and landscape levels.
The main purpose of this study is to investigate how different landscape metrics are affected by changing thematic details of LC/LU maps.This investigation provides clues that can be reflected in the landscape in general, by revealing the causal relationships between different levels.With the aid of CORINE levels, the detailing of classification categories can increase the level of significance of metric interpretations, and also provides a basis for more qualified assessment of the overall landscape.To fulfill this aim, a study area with diverse urban structure is selected in Izmir city of Turkey.Next, two SPOT 7 images of the study area with different dates (winter and spring) are obtained and pre-processed.Object based classification is employed using various features and indexes in order to create the CORINE Level 3 LC/LU map of the area.After that, classes of the Level 3 map were merged according to CORINE nomenclature, in order to create Level 2 and Level 1 maps.Area-based accuracy assessment is performed to determine accuracy of these 3 maps with 200 random areas.For each level LC/LU map, landscape metrics were evaluated and results were compared to objectively evaluate the thematic detail and metric relationship.

Study Area
The study area is located in Izmir metropolitan city, which is located in the Aegean region of Turkey (38 • 24 N and 27 • 10 E).Izmir is the third most populated city of Turkey with the population of 4,279,677 [39].The city attracts national and international visitors with its long coastline along Aegean Sea, warm climate and historical sites.Additionally, the Izmir Port is an important transportation hub which contributes to the city's economy.All of these attractions cause a rapid growth in Izmir in recent years and create an urban sprawl.
The study area covers approximately 286 square kilometers.This region has complex landscape characteristics with several LC/LU classes such as water bodies, agricultural areas, forests and urban areas (Figure 1).The study area covers approximately 286 square kilometers.This region has complex landscape characteristics with several LC/LU classes such as water bodies, agricultural areas, forests and urban areas (Figure 1).

Data
SPOT 6/7 satellites have multispectral sensors with 6 m spatial resolution in Red, Green, Blue and Near-Infrared (NIR) regions of electromagnetic spectrum in addition to 1.5 m resolution Panchromatic sensor [40].In this study, two SPOT 7 images acquired on 24 April 2017 and 8 December 2016 were used as primary data source for creating LC/LU maps.The main reason for using two images with different dates was to take the advantage of capturing the temporal patterns, in addition to spectral responses for the class definitions.Temporal characteristics provide valuable information to extract the LU classes such as agricultural lands and separating them from spectrally similar classes.
Additionally, Open Street Map (OSM) vector data was used as a thematic layer for road extraction.Wikimapia was also used as a vector thematic layer for some of the artificial classes, which are basically representing the land use.Moreover, online maps such as Google Earth, Google Street Viewer and the CORINE 2012 dataset were used for visual interpretation of the study area to form

Data
SPOT 6/7 satellites have multispectral sensors with 6 m spatial resolution in Red, Green, Blue and Near-Infrared (NIR) regions of electromagnetic spectrum in addition to 1.5 m resolution Panchromatic sensor [40].In this study, two SPOT 7 images acquired on 24 April 2017 and 8 December 2016 were used as primary data source for creating LC/LU maps.The main reason for using two images with different dates was to take the advantage of capturing the temporal patterns, in addition to spectral responses for the class definitions.Temporal characteristics provide valuable information to extract the LU classes such as agricultural lands and separating them from spectrally similar classes.
Additionally, Open Street Map (OSM) vector data was used as a thematic layer for road extraction.Wikimapia was also used as a vector thematic layer for some of the artificial classes, which are basically representing the land use.Moreover, online maps such as Google Earth, Google Street Viewer and the CORINE 2012 dataset were used for visual interpretation of the study area to form the decision trees.Lastly, parcel data obtained from Ministry of Food, Agriculture and Livestock of Turkey is used for accuracy assessment of the classification results.

Methods
Figure 2 illustrates the flowchart of the process chain.Firstly, SPOT 7 images (winter and spring) were pre-processed and classified by use of various features and indexes in order to create the CORINE based Level 3 LC/LU map.The reason for choosing winter and spring images was to differentiate between agricultural areas with seasonal crops.In spring time, when crops are in growing stage, vegetation indices such as the Normalized Difference Vegetation Index (NDVI) gives higher results than winter time (after harvest).After that, classes of the Level 3 map were merged according to CORINE nomenclature in order to create Level 2 and Level 1 maps.Area-based accuracy assessment was conducted to determine accuracy of these three maps.Finally, landscape metrics were calculated and evaluated for each LC/LU map.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 5 of 21 the decision trees.Lastly, parcel data obtained from Ministry of Food, Agriculture and Livestock of Turkey is used for accuracy assessment of the classification results.

Methods
Figure 2 illustrates the flowchart of the process chain.Firstly, SPOT 7 images (winter and spring) were pre-processed and classified by use of various features and indexes in order to create the CORINE based Level 3 LC/LU map.The reason for choosing winter and spring images was to differentiate between agricultural areas with seasonal crops.In spring time, when crops are in growing stage, vegetation indices such as the Normalized Difference Vegetation Index (NDVI) gives higher results than winter time (after harvest).After that, classes of the Level 3 map were merged according to CORINE nomenclature in order to create Level 2 and Level 1 maps.Area-based accuracy assessment was conducted to determine accuracy of these three maps.Finally, landscape metrics were calculated and evaluated for each LC/LU map.

Pre-Processing
SPOT 7 images are acquired as ortho-rectified and pan-sharpened products with a 1.5-m spatial resolution.As the first step, image digital number values were converted into top of atmosphere (ToA) reflectance values.Geometric correction was performed for the images with taking the pre-ortho-rectified Pleiades imagery with 1 m locational accuracy as the reference data.Each SPOT 7 image was geometrically corrected using the first order polynomial model and use of homogenously distributed 25 ground control points (GCPs), collected from reference Pleaides imagery.During geometric correction, all RMSE (Root Mean Square Error) values were better than 2 m.
As the last step, different dated SPOT images were layer stacked in order to analyze the impact of different seasons on LC/LU classification simultaneously.Table 1 presents the bands of layer-stacked image.

Classification
In this study, object-based classification was performed using eCognition Developer© v.9.3 software.Performing GEOBIA on the satellite imagery consists of two steps; segmentation and classification.The purpose of the segmentation is to create image objects by grouping similar pixels.The most popular and widely used segmentation technique is multi-resolution segmentation.In this algorithm, each pixel is considered as an image object at first and homogeneous image objects are created by combining them with their neighbors in an iterative approach [41].In this process, five parameters called scale, shape, color, compactness and smoothness are defined by the user, based on the spatial resolution and spectral characteristics of the data and desired object size.The most important parameter is scale, which directly affects the size of the object [42].Using larger scale parameters would create larger objects and vice versa.Shape/color, and compactness/smoothness parameter pairs complement each other to 1, which means decreasing the value of one will increase the value of complementing parameter.Shape parameter is directly related to the geometry; while color, compactness and smoothness parameters are affected by the spectral characteristics of an object [14,42].
In order to obtain the most suitable objects for LC/LU classes with diverse characteristics, several multi resolution segmentations were performed.Parameters used for the segmentation of different classes are presented in Table 2.The reason for using various segmentation scales was the difference between the size and spectral characteristics of classes.For example, the objects belonging to Water Bodies (5) and Sea and Ocean (523) classes covered the biggest area and the objects of Water Courses (511) class were extremely small, compared to 523.Using different shape and compactness parameters changed the composition of objects.During the segmentation and classification processes, 2 thematic layers were used as ancillary vector data in order to obtain better results for some of the LC/LU classes.Table 3 shows the additional vector layers and their purpose in this study.In the classification procedure, using vector layers with related attributes is important to obtain accurate land use (such as schools, hospitals . . . ) information [13,43].This information cannot be determined by using only the satellite images, since creating land use information is a much more complex task than creating land cover information.Specifically, for some classes, considering the fact that same representations in satellite image could represent different land use, ancillary data was needed.For example, some vegetated regions located in the Artificial Areas were classified as Green Urban Areas (141) because they had same spectral characteristics as other 141 objects.However, these areas had different purposes, such as being a sport facility, therefore, these objects would be classified as Sport and Leisure Facilities (142).These examples revealed that, in some cases the actual class can only be identified with some ancillary data, which includes land use and management information.As a result, obtaining land use information using only satellite image data was not as easy as determining the land cover.Therefore, additional open source geo-spatial data were used as thematic layers in the GEOBIA procedure to improve the identification of LU classes.However, boundaries of these objects were deducted from satellite images; therefore, compatibility of thematic layers with satellite image was carefully evaluated.

Thematic Layer Purpose
OpenStreetMap Road Data Segmentation/Classification of Road and rail networks and associated land Wikimapia Open Source Data Segmentation/Classification of all Agricultural Areas sub-classes (in Level-3) After image segmentation, image objects were classified with various features and indices, as explained in Table 4 [13,44,45], in order to create the Level 3 LC/LU map of the study area.
Figure 3 shows the general structure of the classification process and features/indices used for the definition of general classes.It is worth noting that, while general classification procedure follows that structure, there are some sub-classes that required refinement in the end, which will be explained in detail further in this section.Additionally, a 0.25 ha minimum mapping unit rule for sub-classes of Artificial Areas (1) and a 1 ha minimum mapping unit rule for all other classes were applied both in classification and segmentation procedures.The classification process was started with the classification of Level 3 classes of Water Bodies (5) which are Water Courses (511), Water Bodies (512), and Sea and Ocean (523).At first, multi-temporal Normalized Difference Water Index (NDWI) indices were used to determine appropriate thresholds.In addition to NDWI, "Ratio" of both Layer 4 and 8 (Near Infrared bands of 2016 and 2017), "Area" and, "Distance to" functions were used for accurate identification of all Water Bodies (5)."Distance to" function was applied to Sea and Ocean (523) class in order to classify Water Courses (511) near Sea and Ocean (523).Classification of these areas were relatively simple, when compared to other classes.
Secondly, Road and Rail Networks and Associated Land (122) class was created with the help of Open Street Map (OSM) vector data.OSM was very beneficial in the segmentation step since it enabled the delineation of highways and crossroads more accurately."Asymmetry", "Length/Width", "Density", "Border index" functions were used to identify threshold values for (122) due to distinguishing geometries of the members of this class."Related Border to" function applied to (122) class was useful to assure continuity of the roads.
Normalized Difference Vegetation Index (NDVI) was applied as an arithmetic function to distinguish the vegetated and non-vegetated areas.Then, classification of the remaining Artificial Surfaces (1) conducted with the help of Wikimapia and various features shown in Figure 3. "Coordinate" function applied to "X and Y center", "Shape Index", "Rectangular Fit", "Area" and "Brightness" functions were the best performing features for this process.
Next, classification of Agricultural Areas (2) was performed.The difference between Near Infrared bands of 2016 and 2017 images (Layer 4 − Layer 8) gave the best results to identify Non-Irrigated Arable Land (211) because wheat and other crops cultivated in the study area without irrigation were in the growing stage in April and they were harvested in December.This phenological difference is a key information for the identification of different crop types.NDVI 2016, NDVI 2017 and the difference of these NDVI layers were used to identify Permanent Crops (22), Pastures (231) and Heterogeneous Agricultural Areas (24) classes.NDWI 2016, NDWI 2017, "Maximum difference" and "Mean Layer In the next step, after applying a bigger scale segmentation to the unclassified areas, classification of the Broad-Leaved Forest (311) and Coniferous Forest (312) were conducted with functions provided in Figure 3.After that, a new segmentation with a finer scale was executed in order to classify sub-classes of Shrub and/or Herbaceous Vegetation Associations (32) and Open Spaces with Little or No Vegetation (33).NDVI data of 2016 and 2017 images, "Coordinate", "Area", "Maximum difference" functions, "Skewness" and "Standard deviation" functions of Layer 4, were used in this step.
After applying these steps, most of the study area was classified into appropriate classes.However, there were two classes that were classified inaccurately during this process, due to their similarities.Some very narrow Water Courses (511) members mixed with Artificial Surfaces because of their similar spectral characteristics.Also, Green Urban Areas (141) were mixed with other vegetated areas were corrected according to their positions.These classes were further analyzed and classified by using geometrical characteristics.
After the segmentation and classification processes, Level 3 LC/LU map was created and classes of Level 3 were combined according to CORINE nomenclature to obtain Level 2 and Level 1 LC/LU maps.Finally, an error matrix was created by using 200 randomly selected reference areas for each level map.
After applying these steps, most of the study area was classified into appropriate classes.However, there were two classes that were classified inaccurately during this process, due to their similarities.Some very narrow Water Courses (511) members mixed with Artificial Surfaces because of their similar spectral characteristics.Also, Green Urban Areas (141) were mixed with other vegetated areas were corrected according to their positions.These classes were further analyzed and classified by using geometrical characteristics.
After the segmentation and classification processes, Level 3 LC/LU map was created and classes of Level 3 were combined according to CORINE nomenclature to obtain Level 2 and Level 1 LC/LU maps.Finally, an error matrix was created by using 200 randomly selected reference areas for each level map.

Landscape Metrics
Turner et al. (1989) states that qualitative and quantitative changes in measurements across spatial scales will differ depending on how scale is defined [46].There is no single correct or optimal scale to describe spatial heterogeneity.However, since landscape pattern is dependent on the scale, using different scales in analysis could result in different outcomes, which could affect landscape

Landscape Metrics
Turner et al. (1989) states that qualitative and quantitative changes in measurements across spatial scales will differ depending on how scale is defined [46].There is no single correct or optimal scale to describe spatial heterogeneity.However, since landscape pattern is dependent on the scale, using different scales in analysis could result in different outcomes, which could affect landscape management options.If it is desired to make a comparison between the landscape pattern levels, the data used must have the same spatial resolution [47].For that purpose, a more qualified landscape assessment by changing the classification level was conducted without changing the spatial resolution.Therefore, universal and consistent class level metrics [48] were used to determine the spatial structure of the study area because information accumulated in class level from patches does not exist at patch level.Landscape level metrics give information about whole region without detailing class relations.Most metrics at the class level are derived from patch level attributes and integrated over all the patches of a particular class [26,49].The class-level landscape metrics are more effective in defining ecological processes [50].Additionally, metrics, used in some studies on aggregation and subdivision, were sensitive to the number of classes [51].In this study, metric values at different levels of classification were interpreted comparatively.
In order to quantify and monitor large databases of landscape characteristics, landscape complexity must be defined and measured [52].Components that define landscape characteristics have different level of significance depending to classification categories.Thus, metrics to be used in this study are chosen so they could define characteristics in different levels and also be comparatively analyzed in order to assess the complexity in a landscape.
Class metrics used in this study are Percentage of Landscape (PLAND), Number of Patches (NP), Patch Density (PD), Largest Patch Index (LPI), Total Edge (TE), Edge Density (ED), Landscape Shape Index (LSI), Shape Index Area-Weighted Mean (SHAPE_AM), Total Core Area (TCA), Euclidean nearest Neighbor Distance Area-Weighted Mean (ENN_AM), Splitting Index (SPLIT), Aggregation Index (AI).These metrics are described in Table 5 and metrics values were calculated using 8 × 8 m cell neighborhood rule in FRAGSTAT [26] software for each level CLC map.After that, results are evaluated.
Table 5. Class metrics and descriptions [35,38,53].Landscape Shape Index (LSI) A standardized measure of patch compactness that adjusts for the size of the patch Area-Weighted Mean Shape Index (SHAPE_AM)

Metric
Weighting patches according to their size, on contrary to LSI in which the total length of edge is compared to a landscape with a standard shape (square) of the same size and without any internal edge Total Core Area (TCA) The sum of the core areas of each patch (m 2 ) Euclidean Nearest Neighbor Distance Area-Weighted Mean (ENN_AM) Shortest straight-line distance (m) between a focal patch and its nearest neighbor of the same class Splitting Index (SPLIT) The number of patches obtained with subdividing the landscape into equal-sized patches based on the effective mesh size Aggregation Index (AI) The ratio of the observed number of like adjacencies to the maximum possible number of like adjacencies given the proportion of the landscape comprised of each patch type, given as a percentage

LC/LU Maps
In this study, 3 different levels of LC/LU maps were produced according to CORINE nomenclature by applying GEOBIA technique on multi-temporal SPOT 7 images.Figure 4 illustrates the classification results.In the map of Level 1, Artificial Surfaces (1) cover 167 km 2 , Agricultural Areas (2) cover 25 km 2 , Forest and Semi-Natural Areas (3) cover 73 km 2 , and Water Bodies (5) cover 18 km 2 area.According to classification results, Artificial Surfaces is the dominant land cover type within the study area.A variety of artificial classes could be deducted from Level 2 and Level 3 classifications with more thematic details.As an example for Artificial Surfaces, 4 and 9 different artificial sub-classes were obtained in Level 2 and Level 3, respectively for the study area.These thematic details could provide efficient and precise information on urban, transport, industrial, port areas that could be used as an input for variety of different applications ranging from landscape architecture to environmental studies.Semi-Natural Areas (3) also cover a large proportion when compared to other classes.Evaluation of the generated maps showed that class diversity is quite significant especially at the Level 3 LC/LU map.This suggests the difficulty of keeping the accuracy high during the classification process.In addition, the SPOT 7 image was not enough by itself to determine some LC/LU classes of Level 3 LC/LU map, such as Industrial or Commercial Units (121), and Olive groves (223).Inclusion of different thematic layers shown in Table 3 into the classification improved the results.Considering the different number of patches and classes in different levels, landscape metrics of three different levels provided various information.
Selection of the most appropriate LC/LU classification system is an important step and there are different standardized nomenclatures available for this purpose.For most of the LC/LU systems, hierarchical levels are available in which the 1st Level represents general land cover classes and thematic detail, and complexity in class definition increases in 2nd and 3rd Levels.We used 3 Level hierarchical CORINE-based classes in this research.Although generating 1st Level classes such as Artificial Surfaces, Agricultural Areas, Forest and Semi-natural, Water Bodies are easy and could be automatically done by object based classification approaches in conjunction with some spectral indices, classification of Level 2 and Level 3 thematic classes is a more challenging task to automatize.Specifically, for some Level 3 Artificial Surfaces classes, land use information is crucial, which could not be directly produced from satellite images.In such cases, integration of open source geo-information to segmentation and classification stage improved the accuracy of classifications and lead to identify thematically detailed LC/LU classes.It is important to use multi-temporal data from different seasons to accurately identify agriculture and vegetation related 2nd and 3rd level classes such as Broad-Leaved Forest, Coniferous Forest, Non-irrigated Arable Land and Permanent Crops.Texture based Haralick features are used as an important indicator for the identification of some permanent crops such as Vineyards.
significant especially at the Level 3 LC/LU map.This suggests the difficulty of keeping the accuracy high during the classification process.In addition, the SPOT 7 image was not enough by itself to determine some LC/LU classes of Level 3 LC/LU map, such as Industrial or Commercial Units (121), and Olive groves (223).Inclusion of different thematic layers shown in Table 3 into the classification improved the results.Considering the different number of patches and classes in different levels, landscape metrics of three different levels provided various information.

Accuracy Assessment
An area-based accuracy assessment technique was performed to evaluate thematic accuracies of three LC/LU maps.Very high resolution Pleiades image of the area was used as reference data for ground truth selection.First, whole area was divided grid by grid; each grid having an area of 1 ha (100 m × 100 m).Then, 200 grids with pure LC/LU classes within its borders were randomly selected.Finally, these areas were used for the accuracy assessment of three different level LC/LU maps by visual comparison of Pleiades image.Error matrix of each thematic map was created and accuracy results derived from error matrix are presented in Tables 6-8.Overall accuracy values of Level 1, Level 2 and Level 3 were found as 93.50%, 88.00% and 83.50%, respectively.Reduction of thematic details towards to Level 1 resulted in better accuracy values.Level 3 classes are more complicated and detailed; therefore, obtaining accuracies similar to Level 1 was not possible at this level.Water related classes are the most accurate classes in all three levels.Although Artificial Surfaces (1) could be successfully classified with more than 95% accuracy in Level 1, accuracy values decreased at Level 2 and Level 3 since these classes are more complex in terms of spectral response and geometrical characteristics.

Landscape Metrics
Landscape metrics such as Number of Patch (NP), Edge Density (ED), Largest Patch Index (LPI), Euclidean Nearest Neighbor Distance (ENN), Splitting Index (SPLIT) and Aggregation Index (AI) metrics are useful indicators to assess the landscape configuration of related regions by using satellite based LC/LU maps.Landscape metric values change with classification levels; therefore, it is important to select the most appropriate classification level for a given area.Most of the studies conducted in the literature are based on Level 1 LC/LU class metrics, which could only provide very general information about the related region.Calculation of landscape metrics for Level 2 and Level 3 classification scheme leads to better understanding of the configuration and composition of the landscape.Tables 9-11 provide the class-level metric calculation results for each of the three different levels LC/LU maps.
Calculated metric values according to Level 1 classification show that Artificial Surfaces (1) is the dominant class in the study area (PLAND = 56.45).Although the NP value of the Artificial Surfaces (1) class is lower than the NP value of the Forest and semi-natural areas (3) class, the LPI, ED, TCA values are higher than the other classes.These values show that compact and related patches are concentrated on this class.ENN_AM, weighted with respect to the patch size, refers to the neighborhood relations of the patches in terms of proximity.According to this metric value, patches that belong to Artificial Surfaces (1) are in a close relation to each other.Also, SPLIT shows that the scattering in Artificial Surfaces (1) is low.Water bodies (5) was identified as the most passive class in terms of PD, NP and PLAND.Patches are too small and too far from each other to develop an edge-core relationship.This relationship is important as it presents the ability of a patch to interpret stability or change in time and it increases the meaning of shape metrics.
Even though it covers half of the surface compared to Artificial Surfaces (1) class, NP and PD values of the Forest and Semi-Natural Areas class (3) are quite high.When considered with the low level of LPI (6.68), these results suggest that this class is undergoing fragmentation over time.It can be also asserted that the Forest and Semi-Natural Areas (3) class, which has the highest LSI value, has a complex structure and formed by natural geometries with relatively natural lines with natural geometries form it.However, in terms of ED and TCA values, it is behind the Artificial Surfaces (1) class.In particular, TCA value is significantly lower (4039.79)when compared to Artificial Surfaces (15,050.68).This shows that patches of Forest and Semi-Natural Areas class (3), have started to lose their core properties and edge relations have been weakened.
Agricultural Areas class (2) was identified as the class with the highest class fragmentation in Level 1.Although value of NP is much higher compared to class Water Bodies (5), which covers the smallest surface in the study area, Agricultural Areas class (2) has the lowest LPI value (1.39).It also has the lowest value for TCA.These values show that the patches in Agricultural Areas (2) are small, undeveloped and according to the values of SPLIT and ENN_AM they are scattered and not connected.
Water bodies class (5), which covers the smallest area in the study region and also is the weakest class in terms of NP and PD has the closest LSI value (3.27) to the geometric form.ED values of patches that belong to this class were found to be particularly low (1.89).Patches are relatively far apart (ENN_AM = 24.96)and scattering level is very high (SPLIT = 266.54)though not as much as the Agricultural Areas class (2).
All of these evaluations of Level 1 class metrics show that, the patches of the Artificial Surfaces (1) in the landscape are compact, big, interrelated and they have mature edge-core relationships.Results conclude that Artificial Surfaces (1) dominate the landscape.Agricultural Areas (2) are scattered in the region as disjointed small units.Forests and Semi-Natural Areas (3) are suppressed by Artificial Surfaces (1) both in area and patch relation.Water Bodies (5) is the weakest class in terms of area, patch density and patch number.Geometry of the water surfaces in the area exhibits unnatural, geometrical forms.These results indicate that water sources in the area are insufficient in terms of ecology.(12) class in terms of TCA, LSI and LPI metrics.In terms of patch and percentage of landscape Urban Fabric (11), in terms of shape and dominant patch class Commercial, and Transport Unit (12) classes stand out.It had been determined that, the dominant class among the Agricultural Areas classes (21, 22, 23 and 24) is Arable Land (21).The values of PLAND, PD, NP and LPI metrics of this class are relatively higher.Heterogeneous Agricultural Areas (24) and Permanent Crops (22) are the classes that enables the general Agricultural Areas of Level 1 to have a heterogeneous form because of their high LSI values.On the other hand, the most disjointed sub-category in terms of the distance between the patches is Permanent Crops (22), while SPLIT shows that the most scattered is Pastures (23).
It is interesting that the class which causes the patch shape to be more indented (LSI = 600.10) is essentially a linear class.This can be explained by the fact that the Urban Fabric (11) class has a much higher SPLIT ratio than Industrial, Commercial and Transport units (12) classes.A large number of patches of class Urban Fabric (11) are scattered throughout landscape, while class Industrial, Commercial and Transport class (12), which has the lowest SPLIT ratio (57.91), forming more complete and uniform patches (LPI = 13.14) which means that shape irregularities and patch shape has become more heterogeneous.On the other hand, Construction Sites class (131) induces the fragmentation (SPLIT) of class Mine, Dump and Construction Sites class (13) in Level 2. In this class, which stands out with a very low PLAND (0.54), a few patches (NP = 8) are located far away from each other and scattered within the landscape.
Discussion of spatial metrics and their possible contribution to landscape ecology should also address the issue of what these measures are to be compared against [54].The findings of this research showed that selected metric sets are successful for revealing landscape characteristics at different levels of classification.Using detailed class definitions while keeping the scale level constant, enabled a more comprehensive interpretation of landscape characteristics.Additionally, performing these observations based on quantitative data is a valuable contribution for landscape analysis and evaluation.

Conclusions
High-resolution urban LC/LU maps are very important data sources for landscape and city planning, management and decision making progresses, considering the population growth in metropolitan cities, economic and industrial developments and increasing urbanization.Therefore, it is important to create up-to-date, thematically rich LC/LU maps using high-resolution satellite images and geographic object based image analysis.
The bottom-to-top (from Level 3 to Level 1) approach is more feasible than top-to-bottom approach and accurate for the creation of three different level LC/LU maps.Multi-temporal data from different seasons are needed to accurately identify agriculture and vegetation related thematic classes.Texture based Haralick features are useful for the identification of some permanent crops such as Vineyards when these crops have a regular planting pattern such as linear, gridwise.Classification accuracy values of different classes were adequate for this research, however very high spatial resolution (1 m or better) is needed to accurately identify Fruit Trees and Berry Plantations considering their lower accuracy values among all thematic classes.Open source geo-information could be a good data source to enhance the results of GEOBIA; however, accuracy and compliancy of these data sets have to be verified before integrating them to classification procedure.
Evaluation of the same set of landscape metrics calculated for each level LC/LU maps proved that meaningful interpretations based on cause and effect relationships could be achieved with increased thematic detail.Calculation of these metrics (PD, NP, PLAND, LPI, TE, ED, LSI, . . . ) used in the study has been effective for interpreting the spatial process of each LC/LU class at three different levels.The LPI and PLAND metrics determine the effectiveness of the classes at different levels in the landscape.Fragmentation of the related class could be inferred from the evaluation of NP, PD, ENN, SPLIT and AI indices together.ED, TCA and LSI provided information about the qualifications of the units that make up that class.Interpretations made in this way, with the help of metric groups, provide information on the locations of the classes on the landscape, structure of their quality, and the levels of fragmentation at each level.This interpretation approach can provide a more intensive cause-effect relationship-oriented assessment as the class level increases.It also helps clarify the landscape character more precisely, and to evaluate the spatial processes in the landscape in more detail.The 1st Level classification explanations are usually more general.The possibility of such a detailed Patch Mosaic Model (PMM) based landscape evaluation by using 2nd and 3rd Level classification has become one of the most important outputs of this research.It is important to find methods that will help to reveal the landscape character in landscape analysis and evaluations.Moreover, landscape metrics based evaluations are quantitative analysis and complex structured landscapes could be easily understandable with the aid of these quantitative information.
ISPRS Int.J. Geo-Inf.2018, 7, x FOR PEER REVIEW 4 of 21 transportation hub which contributes to the city's economy.All of these attractions cause a rapid growth in Izmir in recent years and create an urban sprawl.

Figure 2 .
Figure 2. Flowchart of the study.

Figure 2 .
Figure 2. Flowchart of the study.
Value of Layer 4" supported the classification of cultivated areas."Texture after Haralick" functions gave best results to differentiate permanent crops from each other and other green areas which are Vineyards (221), Fruit Trees and Berry Plantations (222) and Olive Groves (223).Because, plantation pattern of these crops cause them to have unique texture properties in an image.Texture functions used for this purpose are; Entropy, Homogeneity, Dissimilarity and Contrast in all directions.

Figure 3 .
Figure 3. Flowchart of the classification process.

Figure 3 .
Figure 3. Flowchart of the classification process.

Description
Percentage of Landscape (PLAND)The percentage of the landscape comprised of a particular patch type Number of Patches (NP) Number of patches of corresponding patch type (class) Patch Density (PD) Number of patches of corresponding patch type (class) per unit area Largest Patch Index (LPI) The area (m 2 ) of the largest patch in the landscape divided by total landscape area (m 2 ) Total Edge (TE) The sum of the lengths (m) of all edge segments in the landscape Edge Density (ED) The sum of the lengths (m) of all edge segments in the landscape, divided by the total landscape area (m 2 )

Table 1 .
Layer combinations of stacked image.

Table 3 .
Thematic Layers used in segmentation and classification processes.

Table 4 .
Features and indices used in this study.
X-position and Y-position of the center of an image object.The calculation is based on the center of gravity (geometric center) of the image object in the internal map.Related border toDetermines the relative border length an object shares with neighbor objects of a certain classDistance toThe distance (in pixels) of the image object's center concerned to the closest image object's center assigned to a defined class Texture after Haralick Texture features are used to evaluate the texture of image objects.Texture after Haralick features are calculated from gray level co-occurrence matrix.

Table 6 .
Accuracy assessment results of Level 3 classes.

Table 7 .
Accuracy assessment results of Level 2 classes.

Table 8 .
Accuracy assessment results of Level 1 classes.

Table 11 .
(11)l 1 landscape metrics results.Analysis of Level 2 and Level 3 class metrics show that Urban Fabric(11)and Continuous Urban Fabric (111) classes are dominant classes at their respective levels.At Level 2, Urban Fabric(11)which is dominant in terms of NP, PD, PLAND and ED metrics, is behind Industrial, Commercial, and Transport Unit of Level 3 metric results shows that, classes that dominate Level 3 in terms of PLAND and LPI are forest classes (311, 312 . . ., 333).Forest classes are in large units (LPI = 4.79) with centralized area development (TCA = 2397.37).Coniferous Forest (312) leads the forest classes in terms of both PLAND and LPI, and it also has developed edge-core relations (TCA and ED).Transitional Woodland Shrub (324) has the highest patch values (PLAND, PD, NP and LPI).LSI metric of class Moors and Heathland (322) shows that it is the most natural class among others.The interpretation that applies to the whole forest classes at Level 3 is, they have the largest fragmentation (SPLIT = 157,027.38)and the longest average distance between patches (ENN_AM = 18.37).Another notable class at Level 3 is the Olive Grooves class (223), which increases overall LSI value of Permanent Crops.These results revealed that, this class is in a more heterogeneous form than the other sub classes of Permanent Crops.An overall evaluation of Level 2 and Level 3 metrics shows that, according to PLAND, NP and PD metrics, Level 3 class Continuous Urban Fabric (111) induce the Level 2 class Urban Fabric (11) which means class 111 increases the metric values of class 11.Similarly, Road and Rail Networks and Associated Land class (122) induce Shape Metrics and LPI of Industrial, Commercial and Transport class