Consistency Analysis of Remote Sensing Land Cover Products in the Tropical Rainforest Climate Region: A Case Study of Indonesia

: Land cover changes in tropical rainforest climate zones play an important role in global climate change and the functioning of the Earth’s natural system. Existing research on the consistency of di ﬀ erent land cover products has mainly focused on administrative divisions (continental or national scales). However, the ongoing production of large regional or global land cover products with higher resolutions requires us to have a better grasp of confusing land types and their geographical locations for di ﬀ erent zoning (e.g., geographical zoning) in order to guide the optimization of strategies such as zoning and sample selection in automated land cover classiﬁcation. Therefore, we selected the GlobeLand30-2010, GLC_FCS30-2015, and FROM_GLC2015 global land cover products with a 30-m resolution covering Indonesia, which has a tropical rainforest climate, as a case study, and then analyzed these products in terms of areal consistency, spatial consistency, and accuracy evaluation. The results revealed that (a) all three land cover products revealed that forest is the main land cover type in Indonesia. The area correlation coe ﬃ cient of any two products is better than 0.89; (b) the areas that are completely consistent among the three products account for 58% of the total area of Indonesia, mainly distributed in the central and northern parts of Kalimantan and Papua, which are dominated by forest land types. The spatial consistency of the three products is low, however, due to the complex surface types and staggered distributions of grassland, shrub, cultivated land, artiﬁcial surface, and other land cover types in Java, eastern Sumatra, and the eastern, southern, and northwestern sections of Kalimantan, where the elevation is less than 200 m. Given these results, land cover producers should take heed of the classiﬁcation accuracy of these areas; (c) the absolute accuracy evaluation demonstrated that the GLC_FCS30-2015 product has the highest overall accuracy (65.59%), followed by the overall accuracy of the GlobeLand30-2010 product (61.65%), while the FROM_GLC2015 exhibits the lowest overall accuracy (57.71%). The mapping accuracy of the three products is higher for forests and artiﬁcial surfaces. The cropland mapping accuracy of the GLC_FCS30-2015 product is higher than those of the other two products. The mapping accuracy of all products is low for grassland, shrubland, bareland, and wetland. The classiﬁcation accuracy of these land cover types requires further improvement and cannot be used directly by land cover users when conducting relevant research in tropical rainforest climate zones, since the utilization of these products could lead to serious errors.

However, most of these studies focused on administrative divisions (continental or national scales). In fact, there are some types of surface cover under different zoning (e.g., geographical or ecological zoning) with poor accuracy due to the cognitive standards and complexity of the spatial patterns of ground objects. The current production of land cover products with a high resolution of 30 m in large regions or globally requires a better grasp of the confusing land cover types and their geographical locations in different zones in order to guide the optimization of strategies such as zoning and sample selection during automated land cover classification. Hua et al. [35] performed consistency studies on different land cover data at global and continental scales through climate and altitude zoning. Their results revealed that the spatial consistency of Europe is high, at 66.57%, and the overall consistency of the frost climate is as high as 95%. Research on the consistency of different land cover products in tropical rainforest climate zones is lacking, however, even though land cover changes in tropical rainforest climate zones will impact the greenhouse effect, energy balance, and water transport, thereby affecting climate change on a regional or even global scale. Therefore, land cover and its changes in tropical rainforest climate zones have important implications for global climate change. In addition, existing data sources used in the consistency study of different land cover products focus mainly on products with lower spatial resolution. However, low-resolution products, such as MODIS, CCI-LC2000, and GLC2000 have some limitations when carrying out in-depth research on natural geography, the ecological environment, and global change. Hence, there is an urgent need to evaluate and analyze the consistency of current higher-resolution land cover products. Indonesia is located around the equator and is the largest archipelagic country in the world. Its territory spans Asia and Oceania, and it plays an important role in global strategy. Due to its typical tropical rainforest climate, Indonesia is rich in resources, and its land cover types are complex and diverse. The distributions of population, agriculture, forestry, mining, and other resources on each of its islands are unique.
Therefore, in order to make up for the shortcomings of the existing research, this study used Indonesia, with its typical tropical rainforest climate characteristics, as the research area, and the global land cover products GlobeLand30-2010, FROM-GLC2015, and GLC_FCS30-2015 as the data sources to carry out the analysis of areal consistency, spatial consistency, relative accuracy, and absolute accuracy evaluation methods. The study results can provide valuable reference material for research in the field of global ecological environment change and climate change.

Study Area
Indonesia ( Figure 1) is located in Southeast Asia and is connected to countries such as Papua New Guinea, Timor-Leste, and Malaysia. It consists of 17,508 islands and is the largest archipelagic country in the world. The larger islands include Kalimantan, Sumatra, Sulawesi, and so on. The interior of the islands is mostly rugged mountains and hills, with only narrow plains along the coast. Kalimantan has mountains extending from its center to the west, with vast coastal plains and swampy southern regions. On Sumatra, the mountain range is oblique, extending from northwest to southeast, with hills and wider coastal alluvial plain in the northeast of mountain range. Sulawesi is mostly mountainous, with only narrow plains along the coast. Indonesia mainly has a tropical rainforest climate, and seasonal changes in precipitation due to monsoon. The annual rainfall in the plain areas ranges from 1780 to 3175 mm, while yearly totals in the mountain areas can reach 6100 mm.

Data and Preprocessing
We selected three datasets of higher-resolution global land cover products that are currently available free of charge for consistent research, i.e., the GlobeLand30-2010 (http://www.globallandcover. com/), FROM-GLC2015 (http://data.ess.tsinghua.edu.cn/), and GLC_FCS30-2015 (http://data.casearth. cn/sdo/detail/5d904b7a0887164a5c7fbfa0). There are differences among the data sources, classification systems, and classification methods used in the production of these products. The information regarding the main parameters of the three selected land cover products is listed in Table 1. These products require preprocessing prior to accuracy evaluation and consistency analysis, including data clipping, projection conversion, and classification system merging. In this study, the Geographic Information System (GIS) software ArcGIS developed by the Environmental Systems Research Institute (ESRI) was used to preprocess the original data. We used Indonesia's vector border cropping to obtain three land cover products encompassing the study area. The three original land cover products selected in this study are all in the same reference systems, and their geographic coordinate systems are all WGS84. In order to carry out subsequent accuracy evaluation and consistency analysis of different products, the coordinate system of three products is projected as World Mercator. Since the classification system types as well as the formulations of related standards were inconsistent among the three products (Table 2), this study selected the GlobeLand30-2010 classification system as a reference  to standardize the code of the FROM_GLC2015 and GLC_FCS30-2015 land cover types, ultimately  forming a new classification system (Table 3). The spatial distributions of the three land cover products in Indonesia after preprocessing are shown in Figure 2.

Area Composition Similarity
The correlation coefficient of the land cover type area can quantitatively evaluate the similarity degree of the same land cover type between different data [36]. The formula is: where R i is the area correlation coefficient of two land cover data, i is the land cover type, X i is the total area of type i in land cover dataset X, Y i is the total area of type i in land cover dataset Y, X is the average of the total area of all types in land cover dataset X, Y is the average of the total area of all types in land cover dataset Y, and n is the total number of land cover types.

Spatial Pattern Consistency
In order to intuitively express the spatial pattern consistency of different land cover data, the spatial superposition method was employed to obtain the spatial correspondence of different land cover data pixel-by-pixel. Based on the number of matching coverage types of different land cover data determined pixel-by-pixel, the degree of consistency was classified into three levels, in descending order: (1) High consistency: the coverage types of the three land cover products are exactly the same at the same pixel; (2) Moderate consistency: any two products have the same coverage type at the same pixel; (3) Low consistency: the three land cover products all have different coverage types at the same pixel.

Consistency Distribution of Topographic Features
Topographic and geomorphic characteristics influence the precision of land cover products [37]. In this study, elevation data were introduced in order to analyze the consistency distribution of multi-source remote sensing land cover products at different elevations. The digital elevation model (DEM) data used in the study are Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model Version 2 (ASTER GDEM V2). The product is based on the calculation of ASTER data, and is the only high-resolution elevation image data covering the global land surface. We downloaded the 30 m resolution ASTER GDEM V2 data covering the study area from the Geospatial Data Cloud (http://www.gscloud.cn/). Based on the topographic and geomorphic characteristics of the study area, in order to provide the land cover types with obvious gradient characteristics in the elevation distribution, we divided the elevation into five grades (Table 4).

Relative Precision Evaluation
The precision evaluation of land cover data based on the confusion matrix is the most common and important method used [38][39][40][41]. The precision evaluation indices obtained by the confusion matrix are producer accuracy (PA), user accuracy (UA), overall accuracy (OA), and kappa coefficient. The calculation formulas of each index are as follows [42]: where n is the total number of pixels, r is the number of types, x ii is the number of pixels correctly classified, x +i is the number of pixels of a particular type in the reference data, and x i+ is the number of pixels of a particular type in the data to be evaluated. We arbitrarily selected two of the three land cover products to calculate their confusion matrix and determine the relative accuracy evaluation results of the different products.

Absolute Precision Evaluation
In order to evaluate the absolute accuracy of the three land cover products, this study took the field survey samples from Indonesia in 2018 ( Figure 3) as the reference data, and combined them using manual encryption in order to obtain the final accuracy evaluation verification samples. The field survey samples referenced are collected by GPS. The main steps include five parts: sample point selection, coordinate positioning, taking photos, information recording, and sample point correction. Sample point selection: the basic principles to be followed are that the collected sample points are typical and representative; the collected sample points should choose the larger type patches as far as possible; and the collected sample points try not to be close to the distribution boundary of the type to avoid the influence of the mixed type. Coordinate positioning: use GPS to locate geographical coordinates at selected points. Taking photos: take actual photos at each selected point and save these photos as original materials. Information recording: record the coordinates (latitude and longitude), type, and other related information of each sample point. Sample point correction: refer to high-resolution images (such as Google Earth) and other information to correct the coordinates of the sample points indoors. When obtaining verification samples, there is a certain time interval between the base year of the land cover products to be evaluated and the investigation time. In addition, the Google Earth high-resolution image repository was one of the main data sources for accuracy evaluation given its advantages of accurate positioning, abundant time phase, and wide coverage [43]. Therefore, we transform sample points into the Keyhole Markup Language format (KML) and imported into Google Earth, and then the sample points are interpreted online with Google Earth high-resolution images in 2010 and 2015. If there is no Google Earth historical image at a sample point during interpretation, the sample is discarded. In order to reduce the negative impact of positioning error and interpretation error on sample quality, the following principles were followed during sample selection and interpretation: (1) Since the positioning error of the high-resolution image samples of Google Earth is approximately 15 m, while the spatial resolution of the land cover products to be evaluated is 30 m, the center point of a 210 × 210 m homogenized area was selected as the sample in order to reduce the impact of positioning error [44]; (2) for samples that were difficult to interpret, reference was made to other auxiliary information, such as Geo-Wiki [45]; (3) the interpretation adopted the multi-person independent interpretation method, and samples were abandoned when the interpretation results could not be unified after negotiation.    Figure 5 presents the area statistics of the three land cover products in Indonesia. This figure reveals that all of these products generated basically the same type of Indonesia land cover pattern, i.e., the distribution of forest land cover types in the study area predominates, while the cropland, grassland, shrubland, and other types have smaller distributions. The GlobeLand30-2010 and FROM_GLC2015 products produced highly consistent forest area coverage, with area proportions of 77.39% and 79.07%, respectively. The water and artificial surface types of the three products were highly consistent, with water area proportions of 1.31% (GlobeLand30-2010), 1.26% (FROM_GLC2015), and 0.71% (GLC_FCS30-2015), and the proportions of artificial surface area were 0.85% (GlobeLand30-2010), 1.13% (FROM_GLC2015), and 1.88% (GLC_FCS30-2015). The consistencies of cropland, grassland, shrubland, and other land cover types among the three products were relatively low. By calculating the correlation coefficient of area composition for any two of the three land cover products (Table 5), the area composition correlation coefficient of the GlobeLand30-2010 and FROM_GLC2015 products was determined to be the highest, 0.99, while that of the GLC_FCS30-2015 and FROM_GLC2015 products was found to be the lowest, 0.89.

Consistency Analysis of Spatial Patterns
As the largest country in Southeast Asia, Indonesia is rich in natural resources due to its location in the contact zone between the Eurasian and Pacific plates, a region featuring frequent crustal movements and numerous volcanoes [46]. Figures 6 and 7 illustrate the spatial consistency distributions and area statistics of several major cover types from the three land cover products in the study area. For cropland types (Figure 6a), the highly consistent regions of the three products were mainly distributed on northwestern Java, eastern Sumatra, and southwestern Sulawesi. The percentage of highly consistent area to the total area of the study area was 3.15% (Figure 7a). On eastern, western, and northwestern Kalimantan, central Sumatra, Java, and Nusa Tenggara, the consistency of cropland for the three products was low. For forest types (Figure 6b), all of the products exhibited high consistency on Papua, Maluku, central Sulawesi, and central and northern Kalimantan. From Figure 7b, the percentage of highly consistent area of forest to the total area of the study area was 54.21%. For grassland types (Figure 6c), the three land cover products differed greatly, with the percentage of highly consistent area to the total area of the study area only 0.004% (Figure 7c). The low consistency areas were distributed in the border areas of forest and cropland types on Java and other coastal areas. For the artificial surface type (Figure 6d), the consistency of the three products was found to be relatively low. The highly consistent area was mainly distributed in the Jakarta area of Java, and the consistency area accounted for 0.29% of the study area (Figure 7d).   Based on the consistency analysis of each land cover type, the consistency of all types for the three products in Indonesia was further analyzed (Figures 8 and 9). The results revealed that the area ratios of the three land cover products for high-consistency areas, moderate-consistency areas, and low-consistency areas were 58.14%, 33.25%, and 8.61%, respectively. The high-consistency areas were mainly distributed on Papua Island, with forest land as the main cover type, and in the central and northern regions of Kalimantan island, with the proportions of the high-consistency areas on each island to the total study area totaling 18.05% and 17.05%, respectively. Meanwhile, low consistency of the three land cover products was found in the eastern, southern, and northwestern regions of Kalimantan, and the alluvial plains of eastern Sumatra and Java. The main cover types in these low-consistency areas were cropland, artificial surface, and grassland.

Spatial Consistency at Different Elevations
In order to analyze the influence of topographic and geomorphologic features on the precision of land cover products, the spatial consistency distribution results of the three land cover products obtained in Section 4.2 were superimposed with different elevations to obtain the spatial consistency distribution law under different elevations. Figure 10 shows the distribution law of low consistency for several major land cover types at different elevations. It was discovered that the low consistency of cropland, forest, grassland, and artificial surface types among the three products was mainly concentrated in the elevation range of <200 m. The proportions of the areas of low-consistency regions of each type to the total area of the corresponding low-consistency regions were 79.30%, 83.64%, 72.69%, and 85.25%, respectively. With reference to Figure 12, the consistency of each type of the three land cover products is low in the area with the elevation is <200 m which accounts for 61.19% of the total area of the study area. When the elevation is >200 m, however, the low-consistency area of each type decreases significantly, with the area ratios of the low-consistency areas <≈10%.   Figure 11 shows the distribution law of the overall low consistency area of the three products at different elevations. For all of Indonesia, the low-consistency area was mainly distributed at elevations <200 m, accounting for 80.22% of the total low-consistency area. These areas were mainly located on Sumatra, Java, and Kalimantan. Therefore, as with the distribution law of low consistency of different types at different elevation levels, the overall consistency of the three land cover products was the lowest in the area with the elevation of <200 m which accounts for 61.19% (Figure 12) of the total area of the study area, and the land cover types varied greatly among the different products. The low consistency of the three products across the entire study area was significantly less when the elevation was >200 m, particularly when the elevation was >2000 m, where the low-consistency area only accounted for 1.50% of the total low-consistency area, and forest was the main land cover type.
From the above analysis, the low consistency of each type or total of the three land cover products decreases with the decrease of the percentage of elevation grade area. That is, the low consistency is highest in the range of elevation <200 m with the largest proportion of area, and decreases with the decrease of the percentage of elevation grade area.

Relative Precision
Tables 6-8 are the confusion matrices between combinations of two of the three products, showing that the overall accuracy of any two of the three products was <73%, and the kappa coefficient was less than 0.30. The overall accuracy and kappa coefficient between the FROM_GLC2015 and GlobeLand30-2010 products were the highest, with respective values of 72.24% and 0.22, which is consistent with the conclusion regarding the area correlation coefficients of the different product types obtained in Section 4.1. The analysis of the various land cover types revealed that the consistency of the forest land cover type between the FROM_GLC2015 and GLC_FCS30-2015 products was the highest, and the mapping accuracy was 95.81%, while the mapping accuracy of the forest between the GLC_FCS30-2015 and GlobeLand30-2010 was the lowest, with a value of 70.03%. The consistency levels of the grassland and water types between the FROM_GLC2015 and GLC_FCS30-2015 products were high, with mapping accuracies of 73.33% and 83.31%, respectively. For cropland, the consistency between the GLC_FCS30-2015 and GlobeLand30-2010 products was high, and the mapping accuracy was 68.36%, but the user accuracy was only 29.35%.

Absolute Precision
The absolute accuracy of the three land cover products was evaluated using the sample points obtained by field investigation and artificial encryption. The results revealed that the overall accuracy and kappa coefficient of the GLC_FCS30-2015 product were the highest, with values of 65.59% and 0.55, respectively (Table 9), followed by the GlobeLand30-2010 product, with overall accuracy and kappa coefficient values of 61.65% and 0.49, respectively (Table 10). The overall accuracy and kappa coefficient values of the FROM_GLC2015 product were the lowest, with respective values of 57.71% and 0.46 (Table 11). For the various land cover types, the mapping accuracies of the cropland, forest, and artificial surface types of the GLC_FCS30-2015 product were high, with levels ranging from 83.21% to 99.07%, while the mapping accuracies of the grassland, shrubland, and bareland types were <15%. For the FROM_GLC2015 and GlobeLand30-2010 products, the mapping accuracies of the forest and artificial surface types ranged from 78.50% to 96.25%, while the grassland and shrubland types were lower. In addition, the accuracy of the bareland for GlobeLand30 and the wetland for FROM_GLC is 0, indicating that these two types have not been correctly classified.

Discussion
Given that the tropical rainforest climate zone is the world's largest biological gene bank, changes in its land cover not only reflect its surface natural environment, but also directly affect the global environment, especially human living conditions. Remote sensing mapping technology can now provide data support for related research. In this study, we performed a comprehensive evaluation and analysis of three freely available global land cover products, in general, there is low consistency among these products, possibly as a result of the following: (1) Data sources. Differences in remote sensing images will affect the consistency between different products. The GlobeLand30-2010 product is classified using single Landsat TM/ETM remote sensing images combined with HJ-1A/B remote sensing images, the FROM_GLC2015 product is classified using only single Landsat TM/ETM/OL remote sensing images, while the GLC_FCS30-2015 product is classified using time series Landsat OLI remote sensing images. Our absolute accuracy evaluation revealed that the mapping accuracy of the land cover products obtained with time series Landsat imagery (Table 9) was higher than that obtained with single-period Landsat imagery (Tables 10 and 11). At the same time, the fact that the complex cloudy and rainy weather conditions in the tropical rainforest climate area seriously impact the acquisition of high-quality remote sensing images cannot be ignored, which directly affects the accuracy of subsequent classification. Moreover, we found that the area composition correlation coefficient between the GlobeLand30-2010 and FROM_GLC2015 products (Table 5) as well as the relative accuracy were highest (Table 6). It shows that the inconsistency caused by the dynamic change of land cover is much smaller than that caused by different institutions, different data sources, and different classification methods.
(2) Classification systems. The global land cover classification system was established for global classification by taking full account of global land cover characteristics, which inevitably leads to the limitations of its application to unique geographical locations such as tropical rainforest climate zones [47]. The GlobeLand30-2010 product classification system only includes 10 first-class categories, while the classification systems of the FROM_GLC2015 and GLC_FCS30-2015 products are more refined. There are differences in the definitions of some land cover types among the classification systems of these products. For example, for the shrubs of various vegetation type in the tropical rainforest climate zone, the GlobeLand30-2010 product clearly defines land with shrub coverage > 30% and desert shrub coverage > 10% in desert areas as shrubs, while the GLC_FCS30-2015 and FROM_GLC2015 products only define evergreen shrubland and deciduous shrubland, and do not explicitly provide the vegetation coverage values. Differences in the definitions of these vegetation types resulted in low consistencies of the grassland, shrubland, and bareland types among the three products.
(3) Classification strategies and methods. There are some differences in the classification strategies and methods used in the three global land cover products, which will have an impact on consistency. The GlobeLand30-2010 product is constructed using a single-type classification, followed by an integrated classification strategy. The classification method based on "pixel-object-knowledge" is used to classify each type one-by-one, including pixel-based classification, object-based filtering, and human-computer interaction verification to fully utilize the advantages of various classification algorithms and make full use of knowledge and human experience in order to improve classification quality. In addition, this method effectively reduces the classification errors caused by the same objects having different spectra and different objects having the same spectrum [18]. However, the classification strategies and methods adopted by the GlobeLand30-2010 product not only create a large workload, but also require a great deal of human input. The FROM_GLC2015 and GLC_FCS30-2015 products are generated using random forest classification algorithms that are rated as the most robust in global land cover mapping [48,49]. The FROM_GLC2015 product first divides the Earth into 16 regions, and then DEM data, Slope data, Landsat images, and the first global multi-season sample set [50] are used for training classification. Moreover, nighttime light data are introduced to improve the accuracy of impervious types. The FROM_GLC2015 product, however, is extracted scene-by-scene and stitched scene-by-scene during the entire extraction process. The final classification results exhibit severe banding, which seriously affects the overall accuracy. The GLC_FCS30-2015 product proposes a global spatial-temporal spectral library (GSPECLib) using temporal MCD43A4 [20] reflectance products and CCI_LC2015 [51] surface cover products. When constructing the GSPECLib, CCI_LC coverage products are utilized to provide category information for the spectral features satisfying the conditions. Then, based on the geographic location of the GSPECLib, along with the time series Landsat data and the auxiliary terrain data integrated by the Google Earth Engine (GEE) cloud platform, the random forest multi-temporal classification model is trained region-by-region, and the global 30-m high-resolution surface cover classification results are obtained. The classification technology of the GLC_FCS30-2015 product obtains a continuous spatial distribution pattern of ground objects, effectively eliminating the banding problem in the single-period classification results and achieving higher overall accuracy than the other two products.
(4) Other factors. The consistency of different products is high for land cover types with obvious spectral and texture features (such as water and artificial surfaces). For the various vegetation types in the tropical rainforest climate zone, such as forest, shrubs, and grassland, due to the small differences in spectral and textural characteristics and the similar life forms, different types often display the same spectrum, making it more difficult for optical remote sensing to yield accurate classifications. Therefore, the consistency of these confusing types is low.
Analysis revealed that these global land cover products exhibit better mapping accuracy for the dominant forest types in tropical rainforest climate regions such as Indonesia. Users can select one of the products as their data source for relevant research. However, the special climatic and geographical features in the tropical rainforest climate region, such as cloudy and rainy weather, as well as the common occurrence of the same object having different spectra and different objects having the same spectrum, increase the difficulty of remote sensing interpretation. Therefore, for remote sensing interpretation of the tropical rainforest climate region, only relying on a single optical remote sensing image for classification leads to a very large error. It is therefore necessary to integrate multiple data sources (such as SAR and LiDAR data) in order to avoid interference from clouds and rain. Furthermore, the introduction of auxiliary data (such as DEM data) to prior regional divisions of tropical rainforest climate zones may help to improve classification accuracy.

Conclusions
In order to provide a reference for the selection of suitable land cover data required by many studies in tropical rainforest climate regions, this study performed precision evaluation and consistency analysis of three global land cover products, reaching the following conclusions: (1) the performances of the three products on Indonesia's overall land cover type are basically the same. There is a strong area correlation between different products, with a correlation coefficient better than 0.89. (2) The overall spatial consistency of the three products is low, with the highly consistent area accounting for 58% of the total area of Indonesia. This high-consistency area is mainly distributed on Papua Island, as well as the central and northern parts of Kalimantan, with forest as the main land cover type. The low-consistency area is mainly distributed in the more complex areas with elevations < 200 m, such as eastern Sumatra, Java, and the eastern, southern, and northwestern parts of Kalimantan. The grassland, shrubland, forest, cropland, artificial surface, and other land cover types of these areas are staggered and affected by human activities. Therefore, attempts to improve classification accuracy should focus on these areas. (3) The absolute accuracy evaluation experiment of three products shows that the GLC_FCS30-2015 product displays the highest overall accuracy in Indonesia (65.59%), the GlobeLand30-2010 product exhibits the second-highest overall accuracy (61.65%), and the FROM_GLC2015 product has the lowest overall accuracy (57.71%). For forests and artificial surfaces, the classification accuracy of the three products is higher and land cover users can choose any of these products when conducting related research on these land cover types. The GLC_FCS30-2015 product exhibits better classification accuracy for cropland types, with its mapping accuracy in Indonesia reaching 85.28%. The mapping accuracies of the three products are low for grassland, shrubland, bareland, and wetland types, however, indicating that land cover mappers should focus on ways to improve the accuracy of these types. In general, the accuracy levels of these three global land cover products in Indonesia are not ideal, especially for grassland, shrubland, bareland, and wetland land cover types. Hence, these products are not suitable for land cover change studies in Indonesia.
With the development of remote sensing technology, different production agencies and organizations are constantly introducing new and higher-resolution land cover products. In view of the differences between these multi-source remote sensing land cover products, the use of data fusion technology to improve the accuracy of these products will be the main development trend of global land cover mapping in the future.