Spatial Pattern Consistency among Different Remote-Sensing Land Cover Datasets: A Case Study in Northern Laos

: Comparisons of the accuracy and consistency of different remote-sensing land cover datasets are important for the rational application of multi-source land cover datasets to regional development, or to studies of global or local environmental change. Existing comparisons of accuracy or spatial consistency among land cover datasets primarily use confusion or transfer matrices and focus on the type and area consistency of land cover. However, less attention has been paid to the consistency of spatial patterns, and quantitative analyses of spatial pattern consistency are rare. However, when proportions of land cover types are similar, spatial patterns are essential for studies of the ecological functions of a landscape system. In this study, we used classical landscape indices that quantifies spatial patterns to analyze the spatial pattern consistency among different land cover datasets, and chose three datasets (GlobeLand30-2010, FROM-GLC2010, and SERVIR MEKONG2010) in northern Laos as a case study. We also analyzed spatial pattern consistency at different scales after comparing the landscape indices method with the confusion matrix method. We found that the degree of consistency between GlobeLand30-2010 and SERVIR MEKONG2010 was higher than that of GlobeLand30-2010 and FROM-GLC2010, FROM-GLC2010, and SERVIR MEKONG2010 based on the confusion matrix, mainly because of the best forest consistency and then water. However, the spatial consistency results of the landscape indices analysis show that the three datasets have large differences in the number of patches (NP), patch density (PD), and landscape shape index (LSI) at the original scale of 30 m, and decrease with the increase of the scale. Meanwhile, the aggregation index (AI) shows different changes, such as the changing trend of the forest aggregation index increasing with the scale. Our results suggested that, when using or producing land cover datasets, it is necessary not only to ensure the consistency of landscape types and areas, but also to ensure that differences among spatial patterns are minimized, especially those exacerbated by scale. Attention to these factors will avoid larger deviations and even erroneous conclusions from these data products. in land cover type), or completely consistent (all datasets agree on land cover type). These results were then aggregated to determine the proportion of the study region represented by these three classes of spatial consistency.

spatial patterns across different land cover datasets produced from different data sources or by different methods.
To compare spatial pattern consistency among different land cover datasets, we herein used classical landscape indices, which quantifies spatial distributions in landscape ecology, and then conduct quantitative contrast experiments. Using northern Laos as a case study area, we quantitatively compare spatial pattern consistency among three high-resolution land cover datasets (GlobeLand30-2010, FROM-GLC2010, and SERVIR MEKONG2010). Since spatial distribution is still heavily dependent on spatial scale, we also experimentally compare and analyze spatial pattern consistency at different scales.

Study Area
Laos, which is a landlocked country located in the northern part of the Indo-Chinese peninsula (between 100°05-107°38'E and 13°54'-22°30'), has an area of 236,800 km2. Laos is bordered by the western Yunnan Plateau (Yunnan, China) to the north; by Vietnam to the east; by Thailand to the west; by Myanmar to the northwest; and by Cambodia to the south. Northern Laos is connected by the Mekong River, Asia's most important transnational water system. The land cover research in Northern Laos has important implications for ecosystem productivity, biodiversity, and biogeochemical cycles in the Mekong River Basin. Our study area included nine first-level administrative districts, including the provinces Phôngsali, Luang Namtha, Oudômxai, and Vientiane ( Figure 1). The elevation of Laos is higher in the north and lower in the south; the country is inclined to the southeast in the northwest. Xiangkhoang is the highest province in Laos, with an average elevation of 1200 m. Laos has a tropical and subtropical monsoon climate, with a dry season (November-April) and a rainy season (May-October). During the rainy season, the average temperature is 24.2 °C, the average annual rainfall is about 1700 mm, and the rainfall on the plateau and mountains is about 1300 mm. During the dry season, the average temperature is 27.3 °C and there is almost no rainfall due to the dry, cool northeast wind. The little industry in the study area is primarily wood processing, rice milling, and tin mining. The main crops of Laos include rice, coffee, tobacco, and cotton.

Formatting of Mathematical Components
To compare spatial patterns in northern Laos, we used three global and regional land cover datasets: GlobeLand30-2010, produced by the national geomatics center of China (http://www.globallandcover.com/); FROM-GLC2010, produced by Tsinghua University, China (http://data.ess.tsinghua.edu.cn/); and the regional-scale MEKONG2010 dataset, produced by USAID, NASA, and other teams (https://rlcms-servir.adpc.net/en/landcover/); the main parameters of each dataset are summarized in Table 1. These datasets were selected because they all have 30 m resolution and were produced in the same year. Since spatial distribution is still heavily dependent on spatial scale, it is necessary to analyze them in a multiscale way (For more details, please refer to the Section 3.3).

Data Preprocessing
To compare the consistency of the three datasets, we preprocessed the original data by removing all of the extraneous data (selecting only the study area), converting the projections, and merging the classification systems. We used ArcGIS software to trim the datasets and to generate new datasets with consistent boundaries (i.e., those of the study area). We used the universal transverse mercator (UTM) projection and the world geodetic system-1984 (WGS84) coordinate system.
The use of a classification system is critical for global and regional land cover mapping. Various land cover classification systems have been proposed, based on the ability of remote sensing to acquire the attributes of surface features [30][31][32]. In order to compare different land cover datasets, we required a unified classification system. As the classification systems used by the three land cover data sets were different, the formulations of the relevant standards were not uniform, and errors or inaccuracies may occur if there is no general category correspondence rule for direct comparisons [33]. Therefore, we standardized the code used for each land cover type across the three datasets and merged the individual sets of land cover types to establish a unified classification system. The relationship between the original classification (Table 2) and the unified system is shown in Table 3. Importantly, the classification system used by the GlobeLand30-2010 dataset (i.e., only 10 types of first-class categories) was simpler than the other two classification systems. Therefore, we reduced the complexity of the other two systems to the GlobeLand30-2010 system. Other types of merges were performed according to the GlobeLand30-2010 type definition. In addition, there were cloud coverage areas in the SERVIR MEKONG2010 data, so we remove the area data before merging the types, and other products did not consider the grid, and the final results do not include the grid in the evaluation. The three preprocessed land cover products are shown in Figure 2.

Analysis Based on Confusion Matrix and Spatial Overlay
In order to describe the consistency or confusion degree of the same land cover type among different products, every two of three remote sensing land cover products were selected, and their error matrices were calculated to obtain the consistency and confusion degree [34]. The metrics derived from error matrix are: overall accuracy, producer accuracy, and user accuracy. The calculation formulas are: where is the pixel number of correctly classified of type i; n is the total number of pixels in the study area; + is the total number of pixels of type i in the data to be verified; + is the total number of pixels of type i in the reference data.
To determine the consistency in land cover types between the three datasets, the land cover assumed at each pixel was compared and classified as either completely inconsistent (no agreement in land cover type), basically consistent (some agreement in land cover type), or completely consistent (all datasets agree on land cover type). These results were then aggregated to determine the proportion of the study region represented by these three classes of spatial consistency.

Selection of Landscape Pattern Indices
A landscape pattern is the concrete manifestation of landscape heterogeneity, and is the result of various ecological processes acting on different scales [35,36]. The landscape pattern indices are vital for quantifying and analyzing landscape patterns, dynamic changes in landscape patterns, and the factors driving land cover composition [37][38][39][40]. At present, the number of landscape pattern indices, and the diversity of those indices, are increasing. As a single index cannot accurately reflect the real characteristics of a landscape, two or more representative indices should be used to analyze patterns and changes in the landscape. Due to the relative complexity of regional landscapes and the correlation between different landscape pattern indices, some redundancy or similarity appears in the selection of landscape pattern indices. Thus, these indices may not fully reflect the spatial heterogeneity and diversity of the landscape. Therefore, it is necessary to understand the ecological significance of each landscape pattern index, as well as the landscape structure prerequisites that are required for use. Based on previous studies [41,42] and the characteristics of the study area, we selected the number of patches (NP) as the representative index of landscape or habitat fragmentation; patch density (PD) as the representative index of landscape spatial heterogeneity; the landscape shape index (LSI) as the representative index of shape complexity; and the aggregation index (AI) as the representative index of patch connectivity. Each index was calculated using the software Fragstats version 4.2 [43] based on eight domain rules [44]. The definition of landscape indices and their ecological significance are given in Table 4.
where N is the total number of patches, and A is the total area of the landscape PD represents the number of patches per square kilometer 0  PD , No Upper Limit PD is often used to describe the degree of spatial heterogeneity and fragmentation of the landscape patches. Larger PD values imply spatial patch heterogeneity and patch fragmentation.
landscape shape index (LSI) Formula: where E is the total length of all patch boundaries in the landscape, and A is the total area of the landscape 1  LSI , No Upper Limit LSI reflects the complexity of landscape shape. As LSI approaches 1, the overall landscape shape becomes increasingly simple. When LSI = 1, there is only one patch of this type in the landscape, which is square or close to square. As LSI increases, the more complex the boundary shapes of the landscape patches are.
aggregation index (AI) Formula: Represents the number of pixel nodes of the patch type divided by the maximum number of nodes when the patch type is first aggregated 100 0   AI AI reflects the degree of non-randomness or aggregation of different patch types in the landscape, and shows the connectivity between patches of each landscape type. If a landscape is composed of many discrete patches, the aggregation index is low; when a few large patches dominate, or the same types of patches are highly linked, the aggregation index is high.

Multiscale Comparison of Spatial Pattern Consistence
The results of the spatial pattern analysis are heavily dependent on the spatial scale. This is a very well-known problem in spatial data analysis, called the modifiable area unit problem (MAUP), which leads to the inconsistent conclusions when the unit size and partition of the analysis are different [45]. In order to overcome the MAUP and comprehensively compare the consistency of the spatial patterns, this paper makes a comparison of the three land cover data on different spatial scales. Scaling in landscape ecology includes scaling up and scaling down. Scaling up means to move observations, tests, and simulation results from a fine scale to a larger scale [46]. Here, we focused on scaling up. We chose an appropriate scale based on previous studies [47,48]. We used ArcGIS software to resample the original land cover data using the mode sampling method, and obtained 20 different scales of land cover raster data (30,60,90,120,150,180,210,240,270, 300, 350, 400, 450, 500, 550, 600, 700, 800, 900, and 1000 m).

Consistency by Confusion Matrix and Overlap
We selected two of the three types of land cover products as the reference data and the other product as the data needing evaluation and calculated the confusion matrix. Considering that the error matrix values calculated by any two kinds of data reference data and the data exchange role to be evaluated are the same for consistency analysis, we only calculated the confusion matrix once for any two data. Tables 5, 6, and 7 are the confusion matrices calculated by treating FROM-GLC2010, SERVIR MEKONG2010, and Globeland30-2010 as reference data, and Globeland30-2010, FROM-GLC2010, and SERVIR MEKONG2010 as corresponding data to be evaluated. The experimental results show that Globeland30-2010 and SERVIR MEKONG2010 have the highest overall consistency among land cover data, while Globeland30-2010 and FROM-GLC2010 have the lowest overall consistency. From the comparative analysis of each type of consistency, the forest and water type consistency between the two types of land cover data is higher, and the wetland type consistency is lower. Bareland-type consistency is low between Globeland30-2010 and SERVIR MEKONG2010 land cover data. The consistency of shrub, grassland, and bareland types between Globeland30-2010 and FROM-GLC2010 data is low. The reason for the low consistency between shrubs, grasslands, wetlands, and bareland types between any two types of data is mainly because the classification systems and classification methods used in the three products are different, resulting in the three products being produced on the target recognition. The huge difference is followed by the ambiguity in the classification of shrubs, woodlands, grasslands, and bareland, which also leads to large differences in the recognition of different products in remote sensing. Based on the confusion matrix in Tables 5, 6, and 7, the area consistency of the three land cover data was compared and analyzed. The three land cover data show that the forest areas account for the largest proportion. For other types, the water area consistency of the three land cover data is the highest, accounting for 0.7% of the total area of the study area. The types with large differences are shrublands, grasslands, wetlands, and bareland. SERVIR MEKONG2010 data fail to identify shrublands and grasslands in the study area. The grassland area in Globeland30-2010 is about twice that of FROM-GLC2010, and this difference in shrubland is opposite. The percentage of bareland in Globeland30-2010 is only 0.0005%, while the bareland in FROM-GLC2010 accounts for up to 8.81%. For the wetland type, Globeland30-2010 and SERVIR MEKONG2010 have good data consistency, and the area ratio is 0.003% and 0.006%, respectively.
In order to facilitate the visual analysis of the spatially consistent distribution characteristics of multi-source land cover remote sensing data, the three products are spatially superimposed in Figure 3. The results show that the complete spatial consistency area accounts for 67.83% of the total area of the study area. The northern, central, southeastern, and northwestern areas of the study area are distributed, mainly because of the pattern of forest dominating, and water and cropland being relatively less in the study area. So, the land cover type in the completely consistent area is relatively simple, and the spectral features are simple or express the obvious phenological characteristics. The distribution of the basic consistent areas, accounting for 26.54% of the total area of the study area, conforms to the regional variable characteristics of the land cover category, which is mainly located around the completely consistent area. It is mainly distributed in the southern and western parts of the study area and few in the central and eastern parts. The land cover types in these areas are mainly cropland and forest. The area of the inconsistent area accounts for 5.63% of the total area of the study area, mainly distributed in the eastern, southern, and central areas of the study area. The land cover types in these areas are more complex, and the land surface types show obvious heterogeneity characteristics, and cropland, forest, shrubland, grassland, and bareland are staggered.

Spatial Pattern Consistency Analysis Based on Landscape Indices at the Original Scale
We determined the landscape indices values for each type of land cover in the study area based on each dataset (Table 6). There were some differences in the spatial patterns of croplands across the three datasets. The NP, PD, and LSI values calculated using the FROM-GLC2010 dataset were the highest across all three datasets, while the AI value was the lowest. This indicated that the FROM-GLC2010 data identified more cropland fragmentation and spatial patch heterogeneity. In the FROM-GLC2010 data, patch shapes were more complicated, and the degree of aggregation between the patches was low. In the Globeland30-2010 data, the cropland landscape was not highly fragmented, the patch shapes were simple, and the patches were highly connected. Therefore, the Globeland30-2010 and SERVIR MEKONG2010 datasets were similar with respect to cropland identification.
More than 70% of northern of Laos is forest land, and this was the main type of land cover in the study area. All three datasets showed a high degree of aggregation and patch self-similarity for forest land. However, the three datasets differed in the degree of fragmentation and shape complexity of the forest patches. The FROM-GLC2010 data suggested a high degree of fragmentation and complexity, while the SERVIR MEKONG2010 data suggested the opposite.
The FROM-GLC2010 data identified a high degree of shrubland and grassland fragmentation and spatial heterogeneity, with complex patch shapes and little aggregation between patches. This dataset suggested that shrubland in the study area were small and scattered.
Wetland landscapes accounted for less of the study area. Wetland patch fragmentation, spatial heterogeneity, and patch shape complexity were basically consistent between Globeland30-2010 and FROM-GLC2010 datasets, but patch connectivity was higher in the Globeland30-2010 datasets than in the other two data sets. The spatial patterns of water differed little among the three datasets, especially between FROM-GLC2010 and SERVIR MEKONG2010 data. The dispersion of the water landscape was low in the Globeland30-2010 dataset.
The patch fragmentation, patch shape complexity, and dispersion between patches were low for surface landscapes in the Globeland30-2010 dataset. Surface landscapes in the SERVIR MEKONG2010 dataset were the most fragmented and complex, as compared to the other two datasets, although the degree of surface landscape aggregation was higher than that of the FROM-GLC2010 dataset. The landscape indices for bareland differed greatly among datasets. The FROM-GLC2010 dataset suggested the highest patch aggregation, with high fragmentation and patch shape complexity, while the Globeland30-2010 data suggested the opposite. Note: symbol "-" indicates no data for the study area, G represents Globeland30-2010 data, F represents FROM-GLC2010 data, and S represents SERVIR MEKONG2010 data.

Comparison of Spatial Patterns at Different Scales
Changes in the landscape scale may alter patch boundaries, divide or fuse patches, and may thus alter landscape patterns, leading to corresponding changes in the indices of these patterns [49,50]. The landscape pattern indices of each type of land cover data changed differently as the scale increased.
As the spatial scale granularity of croplands increased, NP, PD, and LSI tended to converge among datasets, but AI values differed (Figure 4). At scales of 30-60 m, NP, PD, and LSI values based on FROM-GLC2010 and SERVIR MEKONG2010 data decreased rapidly; these values decreased more slowly at scales >60 m. Scale changes had little effect on the NP, PD, and LSI values of the Globeland30-2010 data. At 60 m, the AI value of the FROM-GLC2010 and SERVIR MEKONG2010 data suddenly increased, while the AI value of three datasets show a steady decreasing trend at other scales. The three datasets suggested that landscape pattern indices for forest lands have high consistency at increasing spatial scales (Figure 5). At a scale of 60 m, the NP, PD, and LSI values decreased sharply, while the AI value increased sharply. Similar to the cropland landscape, as the spatial scale increased, NP, PD, and LSI tended to converge across datasets, while the AI values differed. As scale increased, wetlands indices values differed among the three datasets ( Figure 6). When the scale increased gradually, the NP and PD values for wetlands based on the FROM-GLC2010 and Globeland30-2010 datasets showed high consistency trends across all scales, and scale had only a very slight influence on the indices values. The LSI value of the three datasets showed little difference in landscape index curve with scale increasing. The effects of scale on AI showed low consistency across the three datasets, and they have different fluctuation points in the whole scale range. The effects of scale on the water landscape pattern indices were consistent among the three datasets (Figure 7). When the scale was 60 m, the NP, PD, and LSI indices values dropped sharply, while other indices were gradually decreasing. The indices values became gradually more similar among datasets with scale increasing.  The increase in spatial scale affected the NP, PD, and LSI bareland indices similarly to the cropland indices across all three datasets (Figure 9). At a scale of 30-150 m, the AI values in Globe-land 30-2010 decrease rapidly and tend to be stable when the scale is larger than 150 m; The turning point of SERVIR MEKONG2010 data is 60 m, while the FROM-GLC2010 data decreases slowly in the whole scale.

Consistency Analysis at the Original Dataset Scale
Confusion matrix and spatial overlay visual analysis (after normalizing the original datasets) indicated serious inconsistencies among different landscape types (i.e., forest, shrublands, grasslands, croplands, and bareland). Shrublands and grasslands were not identified by the SERVIR MEKONG2010 dataset. Although the other two datasets identified shrublands and grasslands, the consistency of this type between them was low. The primary reason for these misclassifications is that there are no obvious differences in life form, spectral, and textural characteristics among shrublands, grasslands, and forest. During the imaging process, it is easy to capture different objects with the same spectrum that are difficult to distinguish. The accuracy of the classification system also affects land cover classification [51]. For example, the vegetation coverage and tree height thresholds are not clearly defined for forests, shrublands, and grasslands in the SERVIR MEKONG2010 dataset. Therefore, this classification system may require further revision to reduce the number of fuzzy concepts used. In addition, croplands are typically distributed as concentrated sheets, shrublands are often mixed with forests and croplands. Differences in vegetation thresholds of the different datasets further aggravates the incongruence. In the future, other auxiliary data should be considered to increase the accuracy of classification. Such problems are common in optical remote-sensing classification applications with multi-temporal spectral characteristics. Therefore, it is necessary to improve the accuracy of land cover classification recognition algorithms and by introducing multi-source data and knowledge.
The above consistency analysis was mainly based on dataset generation standards, confusion matrix, and spatial overlay visual analysis. Fine-scale visual comparisons of spatial consistency are subjective. By using the landscape indices to quantify spatial distributions, many subtle differences can be identified. For example, in the SERVIR MEKONG 2010 datasets, shrublands and grasslands were primarily identified as forest lands, leading to increases in forest land connectivity, patch closeness, and AI. Therefore, the loss of shrublands and grasslands impacts the AI of the forest landscape (as shown in Table 8), affecting ecosystem stability, biodiversity conservation, animal reproduction, and biological resource management. Data users must deliberately choose suitable land cover datasets based on research objectives and research area characteristics. The indices used here might be useful for assessing aspects of ecosystem function, including the determination of spatial distribution characteristics for various species and secondary species in the landscape, changes in the interactions among species, and the stability of synergistic symbiosis. Therefore, when selecting multi-source remote-sensing land cover datasets, researchers must not only consider the consistency of areas and qualitative spatial patterns, but also quantitative spatial pattern distribution, especially for studies of ecosystems and ecosystem functions.

Consistency of Spatial Patterns at Different Scales
Scale strongly influences spatial pattern analyses [52]. By comparing the spatial patterns of the three datasets at different scales, we found that the NP, PD, and LSI indices values differed greatly at scales between 30 m and 150 m. At scales >150 m, the difference gradually decreases and eventually tends to be consistent, but there is still some difference in AI value with scale increasing. For example, in Vientiane Province (Figure 10), indices values varied greatly among datasets at 30 m, due to the different classification methods and classification systems; as scale increased gradually, small patches were integrated into large patches, decreasing the number of patches of each type, and leading to decreases in NP and PD values. Through constant fragmentation and fusion, original patch shapes were destroyed, and new patch shapes tended to be regular. This substantially reduced LSI values. So, differences between the datasets gradually decreased. At scales greater than 150 m, the three datasets differed little in their assessments of ecosystem function, including landscape fragmentation, heterogeneity, and shape complexity. However, at scales <150 m, differences between datasets were large. Thus, the differences among datasets might strongly influence ecosystem-based landscape analyses. In addition, AI values remained incongruent among the three datasets, even at large scales, because the AI index is affected by the richness of the structural components and their spatial allocations [45]. As scale increased, the spatial allocation of new patches caused by grid cell aggregation might have led to the substantial differences in AI values observed among the three datasets. Therefore, studies of ecosystem connectivity should carefully consider the differences among datasets across all spatial scales.
Scale Globeland30-2010 FROM-GLC2010 SERVIR MEKONG2010 Figure 10. Spatial pattern distributions of three land cover datasets forest type at different scales.
In the comparative analysis of spatial patterns at different scales, we found that the values of landscape indices changed abruptly at certain spatial scales. The main turning point was 60 m. It was possible that, when the resolution was close to the intrinsic scale of small terrestrial patches, the number and connectivity of these patches changed substantially. As scale continued to increase, the intrinsic scale of the patch changes, and the indices gradually declines. Population dynamics, biodiversity, and ecosystems are inevitably affected by the spatial pattern of the landscape. Therefore, it is necessary to identify differences in the spatial pattern distributions captured by different datasets at different scales. Indeed, although the overall differences among landscape indices were low with scale increasing, future studies should be mindful of each of these change points.

Spatial Consistency Analysis and Absolute Accuracy Evaluation
Absolute accuracy evaluation results are of great significance for selecting and comparing different land cover products, especially when two or more products are inconsistent. It is a meaningful study to provide a standard for researchers to choose the right data. There has been a lot of research on this [25,53].
However, the study of land cover product consistency also has its scientific value. For example, when the spatial consistency of multiple sets of land cover products in a region is good, it indicates that the mapping accuracy of this area is relatively high, and these land use products can be used at a high degree; if the consistency is low, it indicates that the accuracy of mapping in this area is likely to be low. At this time, we may adopt the following suggestion of Giri et al. [21]: "Data producers may use the areas of spatial agreement for training area selection and pay special attention to areas of disagreement for further improvement in future land cover 60 m 150m 800 m characterization and mapping. Users can conveniently use the findings in the areas of agreement, whereas users might need to verify the information in the areas of disagreement with the help of secondary information." The paper aims to reveal the phenomenon or fact that product consistency comparison should not only consider type or area, but also their spatial pattern. This phenomenon has not been noticed in previous studies like Hua et al. [24]. We used the landscape indices, which can quantitatively express the spatial pattern and perform a quantitative case study on the inconsistency of the spatial pattern. We believe that it is useful to study this phenomenon, because: (1) it can remind producers not only to pursue the accuracy of type and area, but also the accuracy of spatial pattern when producing data products; (2) it also reminds users that when using different land use data, they need to consider the difference of spatial pattern relating to function analysis; (3) Landscape indices can measure the spatial pattern difference, thus can be referred by both producers and users.
Besides, the landscape indices analysis method used in this study can also be applied to the absolute accuracy evaluation of different data products. More specifically, researchers can compare the landscape indices of multiple land cover data with the landscape indices of the correct reference data. The closer the indices are, the better the consistency can reach. However, to achieve this method, we need to pay attention to the following two problems: (1) it is very difficult to obtain the correct reference data covering the whole region in practice, and we need to use a statistical inference method which requires a reasonable layout of sample areas. This is completely different from the traditional layout of sample points and requires a new spatial sampling method; (2) For different types of land cover and application purposes, the focus of the spatial pattern may be different. It is necessary for researchers to select landscape indices with corresponding meanings for specific analysis according to specific purpose. However, the solutions of the above two problems need more innovative ideas and study cases.

Conclusions
To test the spatial pattern consistency of different land cover datasets, we herein used the quantitative landscape indices to compare the spatial pattern consistency at different spatial scales. We used three land cover datasets covering northern Laos in the same year (2010) as a case study. Our results indicated that the relative areas identified for each landscape types were similar across the three datasets, but the recovered spatial patterns differed substantially. With the exception of AI, most of the differences among datasets decreased as spatial scale increased. However, some land types still exhibited a spatial pattern mutation. Our results indicated that, when studying ecosystem function based on landscape patterns, land cover datasets should be chosen according to the research aims. In addition, the impact of changes in spatial pattern, especially scale, requires special attention. In addition, dataset producers should evaluate the accuracy between quantity and spatial patterns. In this way, inconsistencies among datasets might be reduced, and comparability increased.