Next Article in Journal
A Decade of Ground Deformation in Kunming (China) Revealed by Multi-Temporal Synthetic Aperture Radar Interferometry (InSAR) Technique
Next Article in Special Issue
A Semi-Automated Object-Based Gully Networks Detection Using Different Machine Learning Models: A Case Study of Bowen Catchment, Queensland, Australia
Previous Article in Journal
Sensitivity Improvement to Active Piezoresistive AFM Probes Using Focused Ion Beam Processing
Previous Article in Special Issue
Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery

1
School of Geomatics and Urban Spatial Informatics, Beijing University of Civil Engineering and Architecture, Beijing 100044, China
2
Beijing Advanced Innovation Center for Future Urban Design, Beijing University of Civil Engineering and Architecture, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(20), 4430; https://doi.org/10.3390/s19204430
Submission received: 3 September 2019 / Revised: 27 September 2019 / Accepted: 9 October 2019 / Published: 12 October 2019

Abstract

:
Sample size estimation is a key issue for validating land cover products derived from satellite images. Based on the fact that present sample size estimation methods account for the characteristics of the Earth’s subsurface, this study developed a model for estimating sample size by considering the scale effect and surface heterogeneity. First, we introduced a watershed with different areas to indicate the scale effect on the sample size. Then, by employing an all-subsets regression feature selection method, three landscape indicators describing the aggregation and diversity of the land cover patches were selected (from 14 indicators) as the main factors for indicating the surface heterogeneity. Finally, we developed a multi-level linear model for sample size estimation using explanatory variables, including the estimated sample size (n) calculated from the traditional statistical model, size of the test region, and three landscape indicators. As reference data for developing this model, we employed a case study in the Jiangxi Province using a 30 m spatial resolution global land cover product (Globeland30) from 2010 as a classified map, and national 30 m land use/cover change (LUCC) data from 2010 in China. The results showed that the adjusted square coefficient of R2 is 0.79, indicating that the joint explanatory ability of all predictive variables in the model to the sample size is 79%. This means that the predictability of this model is at a good level. By comparing the sample size NS obtained by the developed multi-level linear model and n as calculated from the statistics model, we find that NS is much smaller than n, which mainly contributes to the concerns regarding surface heterogeneity in this study. The validity of the established model is tested and is proven as effective in the Anhui Province. This indicates that the estimated sample size from considering the scale effect and spatial heterogeneity in this study achieved the same accuracy as that calculated from a probability statistical model, while simultaneously saving more time, labour, and money in the accuracy assessment of a land cover dataset.

1. Introduction

Land cover provides basic geospatial information for applications in the fields of global environmental change, natural resources management, carbon and nitrogen cycle, and ecological monitoring [1,2,3]. Because of the continuous earth surface scanning and the correspondingly long-term data archives, satellite remote sensing is proven as an effective way in mapping global land cover and measuring land cover dynamic change [4]. Currently, national and international agencies have successfully created no less than ten global scale land cover datasets with spatial resolutions of 1 km, 500 m, 300 m, 30 m, and 12 m. These existing land cover datasets provide basic geographic information for studying climate, hydrology, environment, ecology, and urban regions [5,6,7,8]. Their accuracies are undoubtedly one of the most concerning issues for the potential users. Although data accuracy or uncertainty information is given at the same time as the land cover product releases [9,10], users in different application fields generally need to verify the accuracy before making a decision in using the land cover product [11,12].
Sampling inspection is a commonly used method for verifying the accuracy of land cover products. It provides reliable information on product quality and uncertainty [13,14]. Determining the sample size, as a key procedure in the sampling scheme, lays the foundation for the later stages of sample displacement and verification. A reasonably-estimated sample size is an effective way to avoid the phenomena of over-sampling or under-sampling [15,16]. In addition, the sample size also affects the number of investigators, and the cost and time of the field investigation. Therefore, the estimation of sample size is not only important theoretically, but also plays a guiding role in scientific research and field work investigation.
Currently, there are three major sample size estimation methods: empirical values, fixed-grid sampling, and calculation from a statistical model. The empirical values are the sample numbers required by researchers to test the accuracy of each category of classified satellite images. For example, Hay provided at least 50 empirical values for samples in each category initially [17], although the sample size could be enlarged with an increase in spatial regions and/or the amount of categories involved in the image classification. Congalton provided 75–100 sample data for each category of classified image, which are common empirical values for testing classified images [18]. Empirical sample sizes are simple and can achieve the goal of sample representation through spatial optimization in the process of sample placement. It is usually used to determine the sample size in remote sensing products, especially in an accuracy assessment of classified images. However, with the sprawl of spatial regions, the spatial heterogeneity can vary to a large extent around the entire region(s).
A fixed grid is another commonly used method for determining a sample size for performing an accuracy assessment. It divides the study area into regular grids with a certain size, for instance, 1 km or 10 km, and a sample from each grid is required [19,20]. Ridder used a 10 km × 10 km grid and randomly selected 9000 samples to assess the accuracy of a global forest dataset [21]. Stehman designed a 5 km × 5 km grid, in which 500 grids were randomly selected by synthesizing information on global climate zones and population density [22,23]. Through designing a sample encryption algorithm, a dataset for the verification of global land cover products was generated. This dataset can be used to validate other 1 km or 500 m resolution land cover datasets. The fixed grid sampling method is easy to implement, but the determination of grid specifications relies more on expert knowledge.
A statistical model is widely used to calculate a sample size. This method is based on theories of traditional probability and statistics [24,25]. Tong calculates sample sizes from two scales, one on the level of map divisions, and the other on the level of map elements [26]. By establishing probability distribution functions in each scale, the calculated sample size overcomes the problem of ‘too strict in map divisions and too wide in map elements’, which commonly exists in classical sampling schemes. Based on probability and statistics theory, Olofsson calculated a sample size for simple random sampling and stratified sampling methods respectively, using a statistical model [27]. As compared with the empirical or fixed-grid method, the statistical model is based on theories of statistics and probability. However, the sample size obtained by this method relies on the expected classification accuracy and sampling error information of the products. One of the main problems is that the same remote sensing product in different test sites with different areas will obtain the same numbers of samples, indicating that the surface spatial characteristics, such as patch numbers, aggregation, and diversity, from different study sites, are not involved in the calculation of sample size.
Until now, the determination of a sample size has primarily relied on expert knowledge or conditional assumptions. This often makes it difficult to ensure the rationality of the sample size. Therefore, this study addresses how to determine a number of samples while considering surface spatial heterogeneity.
Statistical theory is the foundation for determining the sample size. This study derived a sample size estimation using a stratified sampling approach. Then, a multi-level linear sample size estimation model was developed by considering the scale effect and surface spatial heterogeneity, with emphasis on two aspects of these issues. First, a watershed unit with ecological and geographical significance was introduced in this study as the basic spatial unit for performing the accuracy assessment, avoiding the objectivity issues existing in current spatial units of pixels or polygons [28]. Second, landscape indicators were employed to describe the surface heterogeneity and complexity. As the characteristics of the spatial heterogeneity would inevitably affect the sample size used for validating the land cover dataset, this study computed several major landscape indicators and assessed their impacts on the surface heterogeneity in watershed units, thereby reaching the goal of this study (to develop a reasonable model to estimate the sample size).
The remainder of the paper is organised as follows. Section 2 introduces the study area, data sources, and data pre-processing methods. Section 3 describes the commonly used statistical model of sample size estimation, the selection of scale factor and landscape indicators, and the development of the multi-level linear regression model. Section 4 presents results and an analysis of the developed model, and Section 5 provides preliminary conclusions.

2. Study Area and Data

2.1. Study Area

The Jiangxi Province is located in south-eastern China, with a total area of 169,900 square kilometres. It belongs to a subtropical humid climate with abundant rainfall. The landform is dominated by mountains (36%) and hills (42%). The main land use types are forest lands and crop lands, with the forest coverage rate reaching 60%, ranking as first in China. The typical geomorphological and climatic features cause Jiangxi to be covered by various types of land covers. Therefore, we selected the Jiangxi Province as the study site to develop the sample size estimation model [29,30]. The Anhui province which is adjacent to the Jiangxi Province is selected to testify the developed model (Figure 1).

2.2. Data Sources

In this study, the global land cover dataset ‘Globeland30’ from 2010 was selected as the test dataset (http://www.globeland30.org). Globeland30 is a high-resolution land cover mapping product developed by the National Geomatics Centre of China, with a spatial resolution of 30 m. It has attracted the attention of researchers at home and abroad. The Globeland30 dataset includes 10 first-level classes. An approach based on the integration of pixel- and object-based methods with knowledge (POK) was used to extract land categories, effectively improving the classification accuracy [31].
The reference data is the national land use/cover change (LUCC) dataset from 2010 with a scale of 1:100,000, which was provided by the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences (CAS). The data was produced mainly by visual interpretation from remote sensing experts in China and was updated every five years. The LUCC data has six first-level classes and 25 second-level classes. Field investigation has shown that the LUCC classification accuracy is very high, with more than 90% overall accuracy (OA) for the second-level class. Both of these two datasets use the same major data sources of Landsat imagery [32,33]. In addition, the digital elevation model (DEM) dataset of the ‘Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model’ (ASTGTM) (http://datamirror.csdb.cn) in the study area was downloaded to perform the division of the watershed unit.

2.3. Data Processing

2.3.1. Classification System Transformation

The Globeland30 and LUCC datasets must be pre-processed before accuracy assessment, owing to their differences in projection and classification systems. For avoiding area deformation, the Globeland30 dataset was projected to the Albert equal area projection (European Petroleum Survey Group (EPSG) code: 3857), which is the projection system employed for the LUCC dataset. Then, these two products were unified in the classification system. Based on the LUCC definition for each product, the classification transformation is presented in Table 1 [34,35]. Finally, we obtained the test and reference data, as shown in Figure 2.

2.3.2. Digital Elevation Model (DEM) Processing

The DEM was projected into the Albert equal area projection system. A hydrological analysis was implemented by using ArcGIS 10.2 to calculate the watershed units. The watershed unit number was decided by the threshold value of the catchment surface [36,37,38]. After repeated experiments and visualization on the division of the watershed unit, we used a threshold value of flow accumulation of 4,000,000 for the catchment surface to obtain 48 basin units. Some units with small areas emerged adjacent to the unit that has the largest area. Finally, we collected 30 basin units as shown in Figure 3b, to perform the following study.

3. Methodology

3.1. Sample Size Determination from Probability Statistical Model

The estimation of the sample size is dependent on the sampling design method. For a land cover dataset with multiple categories, a stratified random sampling design has often been recommended in related research [39,40,41]. Cochran provided a set of rigorous calculation equations to obtain a sample size from the perspective of statistical theory, which is the foundation of latter research on sampling techniques [27,42,43]. According to the good practices for estimating accuracy recommended by Olofsson et al. [27], the commonly used equation for calculating sample size is:
n = ( W h S h ) 2 V ( y ¯ s t ) + ( 1 / N ) W h S h 2
where, V( y ¯ s t ) is the standard error of the estimated OA, and is generally designated as 0.01, W h is the stratum weight (proportion of area of class i in the map), Sh = p i ( 1 p i )   [43], S h represents the level i standard deviation, and p i is the user accuracy, which can be obtained through experiments [42,44].

3.2. Determination of Variables in a Multi-Level Linear Model

On the basis of the estimated sample size from the statistical model, this study attempted to develop a sample size estimation model by considering the scale effects of spatial units and the spatial heterogeneity reflected from the land cover product. As at least three factors, i.e., sample size from the statistical model, spatial effect, and heterogeneity characteristics, impacted the number of samples, this study employed a multi-level linear regression method to develop the sample size estimation, expressed as follows:
N S = A 0 + n A 1 + C A 2 + L A 3 = A 0 + ( W h S h ) 2 V ( y ¯ s t ) + ( 1 / N ) W h S h 2 A 1 + C A 2 + L A 3
In the above, NS is the estimated sample size, n is the sample size calculated by the probability sampling method (Equation (1)), and C represents the sampling constraint on each watershed unit. Theoretically, the larger the test region is, the more samples are needed, and vice versa. In that regard, the ‘best’ sample constraints for regions with different areas should be discussed in this study. L is a set of spatial heterogeneity factors, including a number of landscape indicators that can indicate information on the fragmentation, diversity, and stability of the watershed units. A0, A1, A2, and A3 are regression coefficients, and will be discussed in Section 4.
With the help of watershed units, we obtained Globeland30 and LUCC data as test and reference data respectively, for each single watershed. This following part mainly describes the calculation and determination of n, NS, C, L, and the other variables needed in Equation (2). The main procedures are shown in Figure 4.

3.2.1. Determining NS, n, and C

Equation (1) in Section 3.1 is used to calculate the sample size n. We can see that all parameters in Equation (1) are available except for S h . To ensure the optimization and specificity of the n value in each watershed unit, S h is determined by stratified random sampling for the land cover products of Globeland30. The main steps are as follows (Figure 5): First, each basin unit in the study area is sampled by the stratified random sampling method. Sample NS values are assigned by 50, 100, 150, …, 1000, i.e., with an interval of 50. Second, a user’s accuracy pi from different samples NS in each basin unit is calculated from the constructed error matrix, while using LUCC as a reference data. Third, S h is obtained for each single basin unit according to the evaluated results.
As an example, Table 2 describes the S h results for the Number 2 basin unit. The value of h ranges from 1 to 8, representing eight different strata or land types of Globeland30. W h is the area weight information of each stratum in the basin unit. The S h for each basin unit can be obtained by the above-mentioned stratified random sampling.
After obtaining the parameter S h , the sample size n of each NS can be calculated from Equation (1). As there are 10 different values of NS in each basin unit, there are 10 different corresponding values of n. The determination of n used in the simulation of the multi-level linear sample size model for a single basin unit should be based on the determination of NS.
Figure 6 shows broken-line maps between the NS values obtained from the successive stratified random sampling of 30 basin units and the OA information. The blue line is the result from the stratified random sampling method, the green line represents the result of the accuracy assessed by the whole sample, the black dotted line gives the allowable precision range under the premise of an absolute error of 0.01, and the red circled part is the selected region where the OA varies with NS in the acceptable precision range. The first NS value in the circled region is regarded as the estimated sample size [45], aiming to obtain a reasonable result in accuracy assessment by using as few samples as possible. The sample size n corresponds to N S . Once NS is determined, n can also be obtained. Finally, we obtained 30 records of NS and n from the 30 basin units.
Sample constraint C describes how large of a region can allocate a sample, which is a factor affecting the number of samples. In this study, the sample constraints of each basin unit were obtained by the following expression:
C = NS/S
Here, S is the basin unit area (km2).

3.2.2. Selection of Landscape Indicators

The landscape index, an index for quantitative analysis of landscape patterns, can measure the type, quantity, shape, spatial distribution, and complexity of the analysis units [46,47]. In recent years, an increasing number of studies have used the landscape index to describe spatial heterogeneity information, although their focus is not on the estimation of sample size, but on the layout of sample points [45,48] or land cover extraction [49,50]. According to the target of estimating the sample size of the surface coverage data, 14 landscape indicators were selected to describe the spatial heterogeneity information of the landscape levels in the watershed units from seven categories: area metrics, contrast metrics, edge metrics, shape metrics, proximity metrics, aggregation metrics, and diversity indexes. The ecological significance and descriptions of these indicators are presented in Table 3.
We need to identify the most representative indicators from the 14 landscape indexes in the 7 categories. All-subsets regression, a commonly used method for feature selection, was employed to select the indicators. By adjusting the values of R2, the ‘best’ model was selected to determine the variables of the fitting model. As shown in Figure 7, we found that the adjusted R2 value of the ‘best’ model was 0.78, and the corresponding landscape indices were the landscape shape index (LSI), contagion index (CONTAG), Shannon’s evenness index (SHEI), area-weighted mean shape index (AWMSI), area-weighted mean patch fractal dimension (AWMPFD), and patch richness density (PRD). In addition, the R2 value of the sample size from the probability sampling theory and sample constraints in the fitting model is greater than 0.7, indicating that both of them are independent in the model for sample size estimation. Therefore, the 6 indicators of LSI, CONTAG, SHEI, AWMSI, AWMPFD, and PRD were selected for the following multi-level regression analysis.

3.2.3. Regression Analysis

Equation (2) indicates that the sample size estimation is supposed to be a multi-level linear function of sample size from probability sampling theory, sample constraints, and spatial heterogeneity. Among them, coefficients A0, A1, A2, and A3 are determined by regression analysis. In this study, an ordinary least squares (OLS) regression model was used to determine these regression coefficients [51,52]. OLS is one of the most commonly used core methods in multi-level linear regression models. Its form is as follows:
Y ^ i = B ^ 0 + B ^ 1 X 1 i + + B ^ k X k i
where, i = 1, 2, …, n, Y ^ i is the predicted value of the dependent variable corresponding to the ith observation, Xki denotes the value of the jth predictive variable corresponding to the ith observation, and B ^ 0 denotes the intercept term, i.e., the predicted value of Y when all of the predicted variables are zero. B ^ k is the regression coefficient of the predictive variable j, i.e., the change of Y caused by Xj changing a unit.
The OLS model obtains regression coefficients by reducing the difference between the real values of response variables and the predicted values, i.e., by minimizing the sum of squares of residual errors:
i = 1 n ( Y i Y ^ i ) 2 = i = 1 n ( Y i ( B ^ 0 + B ^ 1 X 1 i + + B ^ k X k i ) ) 2 = i = 1 n ε i 2

4. Result and Analysis

4.1. Multi-Level Regression

As shown in Figure 3b, there are 30 basin units in the study area. Therefore, we obtained a total of 30 records with optimised sample size NS. To show the correlations between sample size NS and the explanatory variables, a scatterplot matrix was obtained, and is presented in Figure 8. The diagonal area is the density map of the variables, whereas the blue line in the non-diagonal area represents the linear and smooth-fitting curves. From the matrix, every predictive variable has a tilt trend to some extent. The sample size NS decreases with the increase of sample constraints (C), SHEI, and AWMPFD, whereas it increases with an increase of CONTAG and PRD. This means that the relationship between the sample size NS and the independent variables is not a phenomenon of simple positive or negative correspondence.
After analysing the scatterplot matrix, the multi-level linear regression analysis was implemented, and the coefficients are shown in Table 4.
The column for Pr(>t) shows the significance of the regression coefficients of the independent variables, and the column of ‘Significance codes’ represents the degree of significance. The regression coefficients of AWMSI, AWMPFD, and PRD are not significant enough to pass the t test, and the multiple R2 of the model is different from the adjusted R2. This reflects a problem of instability. Therefore, the model needs to be improved. From additional experiments, we found that the regression coefficients of the independent variables are more significant when AWMSI, AWMPFD, and PRD are removed. Although the R2 value decreases, the stability performance improves (Table 5). As a result, only 3 indicators were involved in the final estimation of sample size.
Therefore, the multi-level linear regression equation can be expressed as follows:
N S = 7159 0.255 n 9.261 C + 4.21 L S I + 77.9 C O N T A G + 5085 S H E I
By substituting Equation (1) into Equation (6), we finally obtain the estimation of the sample size by considering the scale effect and spatial heterogeneity characteristics, and it is expressed as follows:
N S = 7159 0.255 ( W h S h ) 2 V ( y ¯ s t ) + ( 1 / N ) W h S h 2 9.261 C + 4.21 L S I + 77.9 C O N T A G + 5085 S H E I
This relationship shows that LSI, CONTAG, and SHEI contribute positively to the sample size NS, whereas the sample size n of probability sampling theory and the sample constraint C contribute negatively to the sample size NS. CONTAG and SHEI contribute the most and play a vital role in the change of sample size. All of the predictive variables explain the variance of 79% of the sample size NS.

4.2. Model Verification

This study used a cross-validation method to test the generalization ability of the established OLS regression model. Cross-validation, as a commonly used model validation technology, has the advantage of high prediction accuracy [53]. Considering the small amount of data used to fit the model, three-fold cross-validation was used. Three-fold cross-validation was used to divide the original sample into three equal-sized sub-samples, one of which was retained as test data for model validation, while the other two sub-samples were used as training data. This process was repeated three times, and each of the three sub-samples was used as the validation data. The average R2 of the three cross-validation results was taken as the final estimation solution, as shown in Table 6. The results show that there are some differences between the original R2 and the three-fold cross-validation R2. The R2 value changes greatly after cross-validation, indicating that the stability between the variables and the generalization ability of the models is less than satisfactory.
In addition to the cross-validation, we applied the method of sample size estimation to the Anhui Province to test the developed model. Using the same method as mentioned above, we divided the Anhui province into 32 watersheds and selected 14 units to perform the model validation (Figure 9). The same datasets as those used in the Jiangxi Province were used to test this model. First, we calculated the sample size using a statistical model, i.e., Equation (1), and the developed model in this study, i.e., Equation (7), respectively. Then, we assessed the accuracy of Globeland30 by using the two above-calculated sample sizes. Finally, we compared the accuracies and the results are presented in Table 7. We can see that the difference of OA obtained from two different approaches is very small in most watersheds, the maximum and minimum difference is 4% and 0.1% respectively, with an average of 1.21%. In contrast to the similar accuracy, the sample sizes calculated by the developed model are smaller than that from the statistical model. For the employed 14 watersheds in the Anhui Province, there are 24,425 sample points calculated by the statistical model, while it is 12,399 computed by the developed model in this work. Compared with the statistical model, our developed model decreases the sample size by 49%. This indicates that a smaller sample size can achieve the same performance as the statistical model by considering the scale effect and surface heterogeneity.

4.3. Relative Importance of Predictor Variables

The OLS regression method is used to analyse the influence of the sample size n of the probability sampling theory, the sample constraints C, and spatial heterogeneity factors on the estimated sample size NS, and the relationship between them is obtained in Equation (7). The standardised regression coefficient method and the relative weight method are employed to evaluate the relative importance of the predictor variables in the multi-level regression analysis. The results are shown in Table 8 and Figure 10.
The standardised regression coefficient is the simplest method for predicting the relative importance of variables. It represents the expected change of response variables caused by the change of one standard deviation of a predictor variable when other predictor variables remain unchanged. Table 8 shows that when the other variables remain unchanged, SHEI changes by a standard deviation, and that the sample size will increase by 2.13 standard deviations, i.e., the most important relative to the sample size NS. In contrast, the sample size n of the probability sampling theory has the least relative importance to the sample size.
As compared with the standardized regression coefficient, the relative weight is a more promising method for predicting the relative importance of factors [54]. It ranks variables according to their contribution to R2. Figure 10 shows the relative importance of each factor. The results show that the sample constraint C accounts for 44% of R2, which is of the greatest relative importance to the sample size NS. It shows that the regional scale cannot be neglected in the sample size determination. SHEI and CONTAG, as spatial heterogeneity factors describing the aggregation and diversity of watershed units, explain 25% and 15% of R2, respectively. The remaining factors are LSI and n. Therefore, in terms of relative importance, the regional scale effect and spatial heterogeneity have an influence on the determination of the sample size in accuracy assessment of remotely sensed land cover products.

4.4. Analysis

Figure 11 is a histogram of the sample size NS obtained from the multi-level linear model and sample size n and based on probability and statistics theory. Figure 12 and Figure 13 are polygonal maps of the landscape index and sample constraint factor C, respectively. LSI reflects the complexity of the overall landscape shape of the watershed units. The larger the value, the simpler the overall landscape shape. CONTAG describes the degree of agglomeration of patch types in the watershed units. The higher the value, the better the degree of agglomeration. SHEI describes the diversity of the patch distribution in the watershed units. The smaller the value, the simpler the patch type, and the smaller the diversity.
NS and n are obtained under the premise that the OA standard deviation V( y ¯ s t ) equals 0.01. According to Figure 11, the value of n is often more than 1000, and most of them are distributed near 1750, whereas the values of NS are much smaller. This shows that the sample size NS, as obtained by the multi-level linear model considering spatial heterogeneity, can save more manpower and material resources in the accuracy assessment of a land cover dataset.
The green squares in Figure 11 indicate the relative size of the basin unit area, i.e., the scale difference. According to Figure 3b, there are evident scale differences among the basin units. In theory, the larger the scale of the watershed unit, the more samples are extracted, and the n values obtained based on the probability and statistics theory should show similar laws. However, the larger the scale of the unit, the smaller the NS values needed, such as the basin unit number (No.) 11, No. 20, and No. 35, and the smaller the unit, the larger the NS values, e.g., basin unit No. 14 and No. 17. According to Figure 12 and Figure 13, we can explain why NS changes with areas of the spatial unit. The values of the LSI, CONTAG, and SHEI in units No. 11, No. 20, and No. 35 are higher, indicating that these few units have a high patch aggregation, uniform patch distribution, and simple landscape shape, and thus weak spatial heterogeneity. Although the scale is large, the required sample size NS is small. However, the values of LSI and CONTAG in units No. 14 and No. 17 are low, and the SHEI values are high. They need a large sample size of NS. Therefore, by analysing the relationship between the sample size NS of the watershed units No. 11, No. 14, No. 17, No. 20, No. 35 and the landscape index, it is further demonstrated that the values calculated from the multi-level sample size estimation model (considering the spatial heterogeneity of land cover) are more reasonable than those calculated from the probability statistical model.
As noted above, C indicates how large of a sample is taken. If the land cover in different areas is homogeneous and the patch types and distribution are the same, the C values are similar. However, in Figure 12, the trend of C variation is similar to that of LSI, because the spatial distribution of the patches in different basin units is different. In an area with simple spatial distribution and weak spatial heterogeneity, NS is small, and the C value is high, i.e., C has a negative correlation to NS, whereas when NS is large and the C value is low, C has a positive correlation to NS.
Where the C values are the same, e.g., in basin units No. 16 and No. 17, because the LSI value of NO. 16 is high, the SHEI value is low, and the spatial heterogeneity of the land cover is weaker than that of No. 17. Because the two scales are similar, the sample size of No. 17 is larger. Although the C values in basin units No. 38 and No. 39 are the same, the LSI values of No. 38 are large, but the CONTAG and SHEI values are less than those of No. 39. Therefore, the difference in spatial heterogeneity between No. 38 and No. 39 is smaller than that between No. 16 and No. 17, resulting in a small sample size gap between No. 38 and No. 39. As a result, C, as a factor affecting the sample size of NS, is also influenced by the surface spatial heterogeneity.

5. Conclusions

The accuracy of a dataset is often the first problem to be considered in scientific research and field applications. Sample size calculation is the first step in performing an accuracy assessment of land cover products from satellite imagery. On the basis of a statistical model for the estimation of sample size, this study establishes a multi-level linear model for estimation of sample size by considering the scale effect and spatial heterogeneity. A watershed unit was introduced to obtain a spatial analysis unit, to avoid subjectivity in the selection of assessment units. Landscape indices were selected to indicate the spatial heterogeneity of the region.
The multi-level linear sample size estimation model shows that:
(1): All predictive variables can explain 79% of the variance of the sample size NS. The coefficients of the predictive variables of the model are significant, indicating that there is a strong relationship between the sample size NS and the independent variables. By comparing the sample size NS from the multi-level linear model with a sample size n based on probability and statistics theory, we see that the sample size of NS is much smaller than that of n. The smaller sample size can achieve the same performance as the statistical model and it contributes to the consideration of surface heterogeneity. The relative importance of the predicted variables in the model is calculated by using standardised regression coefficients and relative weights. The results show that the CONTAG and SHEI indicators (describing the diversity and dispersion of basin units, respectively) are relatively important, followed by LSI, sample constraint C, and sample size n, as calculated from probability sampling theory. According to the validation of the developed model, we can conclude that the smaller sample size from the developed estimation model can achieve the same performance as the statistical model while saving more time, cost, and energy in the accuracy assessment of land cover products.
(2): After performing three-fold cross-validation, the R2 value changes from 0.79 to 0.63. This means that the generalization of the sample size estimation model is still a problem, although the test of the model in the Anhui Province proved the validity of this estimation of sample size. Therefore, we need more work on the improvement and testing of the developed model for sample size estimation.
For a specific work on accuracy assessment, although the model established in this study cannot be directly applied, it is expected to provide an approach for the determination of sample size, by considering the study areas and the characteristics of the surface heterogeneity. In the future, we need to improve the developed model by applying this surface heterogeneity-concerned sample size estimation model to other study sites, aiming to assess the accuracy of a land cover dataset with as few samples as possible.

Author Contributions

H.R., and G.C. designed the research method, processed the land cover products, performed the model development and wrote the manuscript, M.D. proposed the requirement and revised the manuscript.

Funding

This work is financially supported by National Key Research and Development Program of China (No. 2016YFB0501404) and the National Natural Science Foundation (NSFC) (Key Project # 41930650).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maria, B.; Monia, M.; Eman, H.; Jun, C.; Ran, L. The First Comprehensive Accuracy Assessment of GlobeLand30 at a National Level: Methodology and Results. Remote Sens. 2015, 7, 4191–4212. [Google Scholar] [Green Version]
  2. Gallego Pinilla, F.J. Comparing CORINE Land Cover with a more detailed database in Arezzo (Italy). Towards Agric.-Environ. Indic. 2001, 1, 118–125. [Google Scholar]
  3. Ren, H.; Cai, G.; Zhao, G.; Li, Z. Accuracy assessment of the globeland30 dataset in jiangxi province. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, 3, 1481–1487. [Google Scholar] [CrossRef]
  4. Chen, J.; Jin, C.; Liao, A.; Xin, C.; Chen, L.; Chen, X.; He, C.; Gang, H.; Shu, P.; Miao, L. Global land cover mapping at 30 m resolution: A POK-based operational approach. ISPRS J. Photogramm. Remote Sens. 2015, 103, 7–27. [Google Scholar] [CrossRef]
  5. Global Urban Footprint. Available online: http://www.dlr.de/eoc/en/desktopdefault.aspx/tabid-9628/16557_read-40454/ (accessed on 5 August 2018).
  6. Esch, T.; Heldens, W.; Hirne, A.; Keil, M.; Marconcini, M.; Roth, A.; Zeidler, J.; Dech, S.; Strano, E. Breaking new ground in mapping human settlements from space—The Global Urban Footprint. ISPRS J. Photogramm. Remote Sens. 2017, 134, 30–42. [Google Scholar] [CrossRef]
  7. Ban, Y.; Gong, P.; Giri, C. Global land cover mapping using Earth observation satellite data: Recent progresses and challenges. ISPRS J. Photogramm. Remote Sens. 2015, 103, 1–6. [Google Scholar] [CrossRef] [Green Version]
  8. Yang, Y.; Xiao, P.; Feng, X.; Li, H. Accuracy assessment of seven global land cover datasets over China. ISPRS J. Photogramm. Remote Sens. 2017, 125, 156–173. [Google Scholar] [CrossRef]
  9. Friedl, M.A.; Mciver, D.K.; Hodges, J.C.F.; Zhang, X.Y.; Muchoney, D.; Strahler, A.H.; Woodcock, C.E.; Gopal, S.; Schneider, A.; Cooper, A. Global land cover mapping from MODIS: Algorithms and early results. Remote Sens. Environ. 2002, 83, 287–302. [Google Scholar] [CrossRef]
  10. Mayaux, P.; Eva, H.; Gallego, J.; Strahler, A.H.; Herold, M.; Agrawal, S.; Naumov, S.; De Miranda, E.E.; Di Bella, C.M.; Ordoyne, C. Validation of the global land cover 2000 map. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1728–1739. [Google Scholar] [CrossRef] [Green Version]
  11. Tsendbazar, N.E.; Bruin, S.D.; Herold, M. Assessing global land cover reference datasets for different user communities. ISPRS J. Photogramm. Remote Sens. 2014, 103, 93–114. [Google Scholar] [CrossRef]
  12. Foody, G. Harshness in image classification accuracy assessment. Int. J. Remote Sens. 2008, 29, 3137–3158. [Google Scholar] [CrossRef] [Green Version]
  13. Herold, M.; Mayaux, P.; Woodcock, C.E.; Baccini, A.; Schmullius, C. Some challenges in global land cover mapping: An assessment of agreement and accuracy in existing 1 km datasets. Remote Sens. Environ. 2008, 112, 2538–2556. [Google Scholar] [CrossRef]
  14. Curran, P.J.; Williamson, H.D. Sample size for ground and remotely sensed data. Remote Sens. Environ. 1986, 20, 31–41. [Google Scholar] [CrossRef]
  15. Foody, G.M. Sample size determination for image classification accuracy assessment and comparison. Int. J. Remote Sens. 2009, 30, 5273–5291. [Google Scholar] [CrossRef]
  16. Martínez-Abraín, A. Are there any differences? A non-sensical question in ecology. Acta Oecol. 2007, 32, 203–206. [Google Scholar] [CrossRef]
  17. Hay, A.M. Sampling Designs to Test Land-use Map Accuracy. Photogramm. Eng. Remote Sens. 1979, 45, 529–533. [Google Scholar]
  18. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 1991, 37, 270–279. [Google Scholar] [CrossRef]
  19. Liu, H. Sampling Method with Remote Sensing for Monitoring of Cultivated Land Changes on Large Scale. Trans. Chin. Soc. Agric. Eng. 2001, 17, 168–171. [Google Scholar]
  20. Stehman, S.; Selkowitz, D. A spatially stratified, multi-stage cluster sampling design for assessing accuracy of the Alaska (USA) National Land Cover Database (NLCD). Int. J. Remote Sens. 2010, 31, 1877–1896. [Google Scholar] [CrossRef]
  21. Ridder, R.M. Options and Recommendations for a Global Remote Sensing Survey of Forests; Forest Resources Assessment Programme Working Paper 141; FAO: Rome, Italy, 2007. [Google Scholar]
  22. Stehman, S.V.; Olofsson, P.; Woodcock, C.E.; Herold, M.; Friedl, M.A. A global land-cover validation data set, II: Augmenting a stratified sampling design to estimate accuracy by region and land-cover class. Int. J. Remote Sens. 2012, 33, 6975–6993. [Google Scholar] [CrossRef]
  23. Stehman, S.V. Impact of sample size allocation when using stratified random sampling to estimate accuracy and area of land-cover change. Remote Sens. Lett. 2012, 3, 111–120. [Google Scholar] [CrossRef]
  24. Fitzpatrick-Lins, K. Comparison of sampling procedures and data analysis for land-use and land-cover map. Photogramm. Eng. Remote Sens. 1981, 47, 343–351. [Google Scholar]
  25. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data—Principles and Practice, 2nd ed.; Lewis Publishers: New York, NY, USA, 1999. [Google Scholar]
  26. Tong, X.; Wang, Z.; Xie, H.; Dan, L.; Jiang, Z.; Li, J.; Li, J. Designing a two-rank acceptance sampling plan for quality inspection of geospatial data products. Comput. Geosci. 2011, 37, 1570–1583. [Google Scholar] [CrossRef]
  27. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  28. Stehman, S.V.; Wickham, J.D. Pixels, blocks of pixels, and polygons: Choosing a spatial unit for thematic accuracy assessment. Remote Sens. Environ. 2011, 115, 3044–3055. [Google Scholar] [CrossRef]
  29. Yin, S.; Chen, X.; Chuanqing, W.U.; Yao, Y.; Wang, X. Spatial-temporal analysis on the variations of the vegetation in Jiangxi Province based on NDVI time series. J. Huazhong Norm. Univ. 2013, 47, 129–135. [Google Scholar]
  30. Zhao, G.; Cai, G.; Mingyi, D.U. Accuracy Assessment for Land Cover Remote Sensing Mapping Product Based on Landscape Shape Index. Beijing Surv. Mapp. 2017, 1, 271–277. [Google Scholar]
  31. Chen, J.; Chen, J.; Gong, P.; Liao, A.P.; Chao-Ying, H.E. Higher Resolution Global Land Cover Mapping. Geomat. World 2011, 2, 12–14. [Google Scholar]
  32. Wang, S.; Liu, J.; Zhang, Z.; Zhou, Q.; Wang, C. Spatial Pattern Change of Land Use in China in Recent 10 Years. Acta Geogr. Sin. 2002, 57, 523–530. [Google Scholar]
  33. Wang, S.; Liu, J.; Zhang, Z.; Zhou, Q.; Wang, C. Study on spatial pattern and change of land use in recent ten years, China. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002; pp. 2369–2371. [Google Scholar]
  34. Ma, J.; Sun, Q.; Xiao, P.; Wen, B. Accuracy Assessment and Comparative Analysis of GlobeLand30 Dataset in Henan Province. J. Geogr.-Inf. Sci. 2016, 18, 1563–1572. [Google Scholar]
  35. Huang, Y.; Liao, S. Regional accuracy assessments of the first global land cover dataset at 30-meter resolution: A case study of Henan province. Geogr. Res. 2016, 35, 1433–1446. [Google Scholar]
  36. Mao, N.; Liu, W.; Wang, H.; Dai, H. Arcgis 10 Tutorial: From Beginner to Master; Surveying and Mapping Publisher: Beijing, China, 2012. [Google Scholar]
  37. Yi, W.; Yang, P. Determination of drainage area threshold for extraction of DEM-based digital drainage network. Jiangxi Hydraul. Sci. Technol. 2008, 34, 259–262. [Google Scholar]
  38. Zongmei, L.I.; Zhouqin, M.A.; Nie, Q.; Man, W.; Huang, Y. Method of Ecological Watershed Partitioning. J. China Hydrol. 2017, 37, 27–30. [Google Scholar]
  39. Stehman, S.V.; Zhang, J.; Foody, G.M. Sampling designs for accuracy assessment of land cover. Int. J. Remote Sens. 2009, 30, 5243–5272. [Google Scholar] [CrossRef]
  40. Cakir, H.I.; Khorram, S.; Nelson, S.A.C. Correspondence analysis for detecting land cover change. Remote Sens. Environ. 2006, 102, 306–317. [Google Scholar] [CrossRef]
  41. Olofsson, P.; Kuemmerle, T.; Griffiths, P.; Knorn, J.; Baccini, A.; Gancz, V.; Blujdea, V.; Houghton, R.A.; Abrudan, I.V.; Woodcock, C.E. Carbon implications of forest restitution in post-socialist Romania. Environ. Res. Lett. 2011, 6, 45–202. [Google Scholar] [CrossRef]
  42. Badjana, H.M.; Olofsson, P.; Woodcock, C.E.; Helmschrot, J.; Wala, K.; Akpagana, K. Mapping and estimating land change between 2001 and 2013 in a heterogeneous landscape in West Africa: Loss of forestlands and capacity building opportunities. Int. J. Appl. Earth Obs. Geoinf. 2017, 63, 15–23. [Google Scholar] [CrossRef]
  43. Cochran, W.G. Sampling Techniques; John Wiley & Sons.: New York, NY, USA, 1977. [Google Scholar]
  44. Olofsson, P.; Foody, G.M.; Stehman, S.V.; Woodcock, C.E. Making better use of accuracy data in land change studies: Estimating accuracy and area and quantifying uncertainty using stratified estimation. Remote Sens. Environ. 2013, 129, 122–131. [Google Scholar] [CrossRef]
  45. Chen, J.; Chen, J.; Hao, W.; Hou, D.Y.; Zhang, W.W.; Zhang, J.; Zhou, X.G.; Chen, L.J. A landscape shape index-based sampling approach for land cover accuracy assessment. Sci. China Earth Sci. 2016, 59, 1–12. [Google Scholar] [CrossRef]
  46. Crews-Meyer, K.; Hudson, P.; Colditz, R.R. Landscape Complexity and Remote Classification in Eastern Coastal Mexico: Applications of Landsatâ-7 ETM+ Data. Geocarto Int. 2004, 19, 45–56. [Google Scholar] [CrossRef]
  47. Neel, M.C.; Mcgarigal, K.; Cushman, S.A. Behavior of class-level landscape metrics across gradients of class aggregation and area. Landsc. Ecol. 2003, 19, 435–455. [Google Scholar] [CrossRef]
  48. Huang, D.; Ke, C.; Wang, Z.; Shuang, L.; Information, S.O.; University, S.O. Accuracy assessment method for remote sensing image classification results based on spatial sampling theory. Comput. Appl. Softw. 2016, 33, 190–194. [Google Scholar]
  49. Sun, W.; Du, B.; Xiong, S. Quantifying Sub-Pixel Surface Water Coverage in Urban Environments Using Low-Albedo Fraction from Landsat Imagery. Remote Sens. 2017, 9, 428. [Google Scholar] [CrossRef]
  50. Sun, W.; Halevy, A.; Benedetto, J.J.; Czaja, W.; Li, W.; Liu, C.; Shi, B.; Wang, R. Nonlinear Dimensionality Reduction via the ENH-LTSA Method for Hyperspectral Image Classification. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 375–388. [Google Scholar] [CrossRef]
  51. Cai, G.; Du, M.; Gao, Y. City block-based assessment of land cover components’ impacts on the urban thermal environment. Remote Sens. Appl. Soc. Environ. 2019, 13, 85–96. [Google Scholar] [CrossRef]
  52. Kabacoff, R. R in Action: Data Analysis and Graphics with R, 2nd ed.; Manning Publications: Shelter Island, NY, USA, 2011. [Google Scholar]
  53. Kohavi, R. A Study of Cross-Validation and Bootstrap for Accuracy Estimation and Model Selection. Int. Jt Conf. Artif. Intell. 1995, 2, 1137–1143. [Google Scholar]
  54. Johnson, J.W. Factors Affecting Relative Weights: The Influence of Sampling and Measurement Error. Organ. Res. Methods 2004, 7, 283–299. [Google Scholar] [CrossRef]
Figure 1. The study site (a) locations in China of (b) the Jiangxi Province and the Anhui Province.
Figure 1. The study site (a) locations in China of (b) the Jiangxi Province and the Anhui Province.
Sensors 19 04430 g001
Figure 2. Globeland30 dataset and the reclassified land cover from land use/cover change (LUCC).
Figure 2. Globeland30 dataset and the reclassified land cover from land use/cover change (LUCC).
Sensors 19 04430 g002
Figure 3. (a) 30 m ‘ASTER Global Digital Elevation Model’ (ASTGTM) digital elevation model (DEM) data and (b) the watershed units (calculated from a).
Figure 3. (a) 30 m ‘ASTER Global Digital Elevation Model’ (ASTGTM) digital elevation model (DEM) data and (b) the watershed units (calculated from a).
Sensors 19 04430 g003
Figure 4. Procedures of parameter determination used in a multi-level linear regression method.
Figure 4. Procedures of parameter determination used in a multi-level linear regression method.
Sensors 19 04430 g004
Figure 5. Flowchart of determining NS, n, and C.
Figure 5. Flowchart of determining NS, n, and C.
Sensors 19 04430 g005
Figure 6. Determination of sample size Ns.
Figure 6. Determination of sample size Ns.
Sensors 19 04430 g006
Figure 7. Result of all-subsets regression method.
Figure 7. Result of all-subsets regression method.
Sensors 19 04430 g007
Figure 8. Scatterplot matrix between sample size NS and the predictive variables.
Figure 8. Scatterplot matrix between sample size NS and the predictive variables.
Sensors 19 04430 g008
Figure 9. The divided watershed units in the Anhui Province.
Figure 9. The divided watershed units in the Anhui Province.
Sensors 19 04430 g009
Figure 10. Relative importance of the predictor variables used in this study.
Figure 10. Relative importance of the predictor variables used in this study.
Sensors 19 04430 g010
Figure 11. Histogram of sample sizes n and NS.
Figure 11. Histogram of sample sizes n and NS.
Sensors 19 04430 g011
Figure 12. Diagrams of landscape shape index (LSI), contagion index (CONTAG), and C.
Figure 12. Diagrams of landscape shape index (LSI), contagion index (CONTAG), and C.
Sensors 19 04430 g012
Figure 13. Diagram of Shannon’s evenness index (SHEI).
Figure 13. Diagram of Shannon’s evenness index (SHEI).
Sensors 19 04430 g013
Table 1. Land cover dataset reclassification and corresponding relation.
Table 1. Land cover dataset reclassification and corresponding relation.
Globeland30Land Use/Cover Change (LUCC)
IndexClass NameIndexClass Name
10Cultivated Land11Paddy Land
12Dry Land
20Forest21Forest
23Woods
24Others
30Grassland31Dense Grass
32Moderate Grass
33Sparse Grass
40Shrubland22Shrub
50Wetland46Bottom Land
64Swampland
60Water Bodies41Stream and Rivers
42Lakes
43Reservoir and Ponds
80Artificial Surfaces51Urban Built-up
52Rural Settlements
53Others
90Bare Land65Bare soil
67Others
Table 2. Parameter values from watershed unit Number 2.
Table 2. Parameter values from watershed unit Number 2.
WhNS
501001502002503003504004505005506006507007508008509009501000
S10.230.470.410.420.490.490.410.450.410.390.460.450.470.430.440.390.460.470.400.450.44
S20.670.410.290.320.330.250.280.340.300.350.310.290.320.300.340.320.320.310.290.310.33
S30.05-0.000.330.330.270.000.000.250.000.350.280.200.260.400.290.350.310.220.330.25
S40.00--------------------
S50.020.470.500.470.490.000.350.500.330.330.400.490.500.480.430.490.480.460.480.490.43
S60.02-0.000.470.500.350.500.500.450.470.460.500.490.480.500.390.500.450.500.500.47
S70.020.000.000.490.470.500.490.500.470.480.490.490.470.490.480.500.360.390.450.460.43
S80.00---------0.00-0.00-0.00-0.00--0.00-
Note: S 1 , …, S 8 corresponds to 10, 20, 30, 40, 50, 60, 80 and 90 land cover categories in Table 1 respectively, - represents that the stratum is not sampled in Globeland30.
Table 3. Landscape indicators used in this study.
Table 3. Landscape indicators used in this study.
ClassNameUnitRange
Area MetricsLargest Patch Index (LPI)%(0,1)
Contrast MetricsMean Patch Size (MPS)hm2>0
Edge MetricsEdge Density (ED)m/hm2--
Patch density (PD)-->0
Shape MetricsLandscape shape index (LSI)--≥1
Area-weighted mean shape index (AWMSI)--(1,2)
Area-weighted Mean patch fractal dimension (AWMPFD)--[1,2]
Proximity MetricsMean proximity index (MPI)-->0
Diversity MetricsShannon’s diversity index (SHDI)--≥0
Patch richness density (PRD)-->0
Shannon’s evenness index (SHEI)--(0,1)
Aggregation MetricsInterspersion and Juxtaposition index (IJI)%(0,100)
Contagion index (CONTAG)%(0,100)
Table 4. Multi-level regression coefficients of sample size NS.
Table 4. Multi-level regression coefficients of sample size NS.
EstimateStandard Errort ValuePr (>|t|)Significance Codes
(Intercept)−17,8805846−3.0585.97 × 10−30.01
n−0.2780.076−3.6581.47 × 10−30.01
C−9.3921.344−6.9876.72 × 10−70.001
LSI5.4791.0265.342.70 × 10−50.001
CONTAG109.90023.3804.6991.22 × 10−40.001
SHEI716314654.897.78 × 10−50.001
AWMSI−7.4244.932−1.5051.47 × 10−11
AWMPFD580330921.8777.45 × 10−20.1
PRD48,49020,1602.4052.55 × 10−20.05
Residual standard error95.28 on 21 degrees of freedom
Multiple R2 squared0.8384
Adjusted R2 squared0.7768
F-statistic13.62 on 8 and 21 DF
p-value<9.80 × 10−7
Table 5. Modified multi-level regression coefficient of sample size NS.
Table 5. Modified multi-level regression coefficient of sample size NS.
EstimateStandard Errort ValuePr(>|t|)Significance Codes
(Intercept)−71591632−4.3861.98 × 10−40.001
n−0.2550.070−3.6271.34 × 10−30.01
C−9.2611.366−6.7815.16 × 10−70.001
LSI4.2100.7625.5231.11 × 10−50.001
CONTAG77.90016.9604.5921.17 × 10−40.001
SHEI5085960.5005.2941.98 × 10−50.001
Residual standard error100.7 on 24 degrees of freedom
Multiple R2 squared0.7935
Adjusted R2 squared0.7505
F-statistic18.45 on 5 and 24 DF
p-value1.58 × 10−7
Table 6. Cross-validation output results.
Table 6. Cross-validation output results.
Original R-square0.7935414
Three-fold cross-validated R-square0.634445
Change0.1590964
Table 7. Comparison of sample size and accuracy calculated from statistical model and this study (No. = Number).
Table 7. Comparison of sample size and accuracy calculated from statistical model and this study (No. = Number).
MethodMulti-LEVEL ModelStatistical ModelAbsolute Value of OA
Difference (%)
Sample Size Difference
(NS-n)
Units Code NSOA (%)nOA (%)
No.1_Ah47579.62151179.560.061035
No.2_Ah67579.56138879.540.02714
No.3_Ah128856.06156156.630.57274
No.4_Ah111068.38162068.400.02510
No.5_Ah48577.11175977.370.261274
No.6_Ah95675.10155676.561.45601
No.7_Ah92273.64159169.943.71669
No.8_Ah87966.78204567.630.851166
No.9_Ah97869.19197967.661.531002
No.10_Ah110659.19190259.940.75796
No.11_Ah69652.65218648.863.801490
No.12_Ah112864.69138466.161.47256
No.13_Ah91261.95200663.241.291093
No.14_Ah79063.29193564.501.201145
Sum = 12,399 Sum= 24,425 Average = 1.21Sum = 12,026
Table 8. Standardised regression coefficients result.
Table 8. Standardised regression coefficients result.
Independent VariablenCLSICONTAGSHEI
Standardised regression coefficients−3.60 × 10−1−8.13 × 10−19.90 × 10−11.972.31
Relative importance level54321

Share and Cite

MDPI and ACS Style

Ren, H.; Cai, G.; Du, M. Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery. Sensors 2019, 19, 4430. https://doi.org/10.3390/s19204430

AMA Style

Ren H, Cai G, Du M. Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery. Sensors. 2019; 19(20):4430. https://doi.org/10.3390/s19204430

Chicago/Turabian Style

Ren, Huiqun, Guoyin Cai, and Mingyi Du. 2019. "Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery" Sensors 19, no. 20: 4430. https://doi.org/10.3390/s19204430

APA Style

Ren, H., Cai, G., & Du, M. (2019). Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery. Sensors, 19(20), 4430. https://doi.org/10.3390/s19204430

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop