Next Article in Journal
Study on Local to Global Radiometric Balance for Remotely Sensed Imagery
Previous Article in Journal
Hyperspectral Identification of Chlorophyll Fluorescence Parameters of Suaeda salsa in Coastal Wetlands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach to the Quality Control of Slope and Aspect Classes Derived from Digital Elevation Models

1
Department Statistic and O.R., University of Jaén, 23071 Jaén, Spain
2
Department Cartography, Geodesy, and Photogrammetry Engineering, University of Jaén, 23071 Jaén, Spain
3
Department Statistic and O.R., University of Sevilla, 41012 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(11), 2069; https://doi.org/10.3390/rs13112069
Submission received: 5 March 2021 / Revised: 12 May 2021 / Accepted: 18 May 2021 / Published: 24 May 2021

Abstract

:
The usefulness of the parameters (e.g., slope, aspect) derived from a Digital Elevation Model (DEM) is limited by its accuracy. In this paper, a thematic-like quality control (class-based) of aspect and slope classes is proposed. A product can be compared against a reference dataset, which provides the quality requirements to be achieved, by comparing the product proportions of each class with those of the reference set. If a distance between the product proportions and the reference proportions is smaller than a small enough positive tolerance, which is fixed by the user, it will be considered that the degree of similarity between the product and the reference set is acceptable, and hence that its quality meets the requirements. A formal statistical procedure, based on a hypothesis test, is developed and its performance is analyzed using simulated data. It uses the Hellinger distance between the proportions. The application to the slope and aspect is illustrated using data derived from a 2 × 2 m DEM (reference) and 5 × 5 m DEM in Allo (province of Navarra, Spain).

Graphical Abstract

1. Introduction

A Digital Elevation Model (DEM) is a generic term for digital topographic and/or bathymetric data [1] and normally implies elevations of the terrain (bare earth) void of vegetation and man-made features. DEMs describe the morphological characteristics of the landforms of a terrain (topography) such as hills, ridges, cliffs, valleys, rivers, etc., and a large number of parameters can be derived (e.g., slope, aspect, curvature, etc.) [2]. DEM data can be generated through many diverse technologies (classical surveying and photogrammetry, Lidar and InSAR systems), and topographic surface modeling can be achieved using very different options such as contour lines, triangular irregular networks, gridded surfaces, parametric surfaces, etc. [3].
We are going to deal with the gridded-model case as it is currently the most widely used for DEM treatment [4]. Grid models are in conformity as a regular-rectangular array of cells, each of which stores an attribute value, so the area of each cell adopts the same elevation value, producing a sudden change in the limits of each cell (a raster model) [5]. The elevation is represented by patches, which are the cells of the grid. Thus, the size of the cell, or resolution, is a key parameter. This parameter is related to the traditional concept of scale and detail of cartographic products. Small cell sizes, giving high-spatial-resolution, are linked to large scales and great spatial detail. Large cell sizes, giving small resolution, are linked to small scales and less spatial detail.
DEMs are geospatial data products that are mostly produced by national mapping agencies (e.g., the Instituto Geográfico Nacional (IGN) in Spain, the Ordnance Survey in the U.K., etc.) and in most developed countries, the territory is covered with high-spatial-resolution DEMs. DEMs are data products at the service of entrepreneurship, public administrations, the private sector, researchers, universities, etc, to generate added value [6].
DEMs are data products with applications in an extensive set of disciplines such as civil engineering [7], hydrology [8], geomorphology [9], agriculture [10], forestry [11], etc. Some of their applications include modelling for the prevention of natural disaster (floods, fires, etc.), soil erosion studies, weather forecasting, climate change, etc. DEMs have been included as a data theme for Europe (Inspire) and as a global fundamental geospatial data theme by the United Nations Committee of Experts on Global Geospatial Information Management. Therefore, one of the fundamental premises is the need for DEMs to have sufficient quality to meet the requirements of those applications.
Two major uses of DEMs are the calculation of the slope and aspect [3,5,12]. Using relatively simple algorithms, slope and orientation models are derived from the DEMs. The calculation of these parameters is based on a local analysis of the elevation derivative. So, neither of them is directly measurable on a single point on the ground (a neighborhood is required). However, while it is not very usual, their direct measurement in the field is also possible [13,14]. Another peculiarity of the calculation of slope and orientation is that small variations in the elevations can generate notable variations in the final values of slope and orientation. The slope is relevant for engineering works and natural hazards (e.g., landslides) [15]. Orientation is related to insolation and therefore has great relevance for plants and the use of solar energy [16].
Given two digital slope or orientation models, two scenarios can be considered. In the first one, the two models to be compared have the same consideration and the interest is focused on the similarity analysis between them, while in contrast, the second scenario is related to the quality control process when one model is considered as a reference and the other must achieve the quality requirements that arise from the reference. In either of the two contexts, to overcome the possible drawbacks derived from the method used to calculate slope and orientation from elevations, it is more robust to work with a categorization of slope or orientation than with the original magnitudes which are continuous and circular (slope and orientation, respectively).
The present work is centred in the second scenario, and we propose to perform a thematic-like quality control (class based) of slope and aspect classes through a user-defined categorization process of these initial magnitudes. A statistical hypothesis test is proposed for controlling the thematic-like accuracy. For this task, it is necessary to assume that there is available a reference DEM (e.g., a sample of greater accuracy). This assumption is not as restrictive as it may seem at first sight because it does not imply a complete knowledge of the reference dataset, it is sufficient to know or to have an estimation of the proportions associated with the categorization used. Because we work with a class-based control, proportions and a multinomial statistical law are considered. In this way, the similarity between the reference and the product can be measured by comparing proportions. The novelty of our approach is that agreement will be understood as both proportions are equal except for irrelevant deviations, which means it is allowed a small enough discrepancy between the observed data and the reference ones. This is the underlying idea of an equivalence test that will be proposed in Section 2.
The rest of the paper is organized as follows. Section 2 introduces how a thematic-accuracy quality control of slope and aspect classes can be performed. To do this, an equivalence test for multinomial data based on the square Hellinger distance is developed. Because the critical region is built by using asymptotic results, that is, it is valid for infinite sample sizes, several simulation experiments are carried out to study the behaviour of the proposal for finite sample sizes. The result of this simulation analysis will help the practical use of the proposed method. As a practical example, Section 3 includes the application of the proposed statistical test to accept/reject the similarity between the observed classifications for two slope and aspect models for the area of Allo (province of Navarra, Spain). Section 4 discusses the main characteristics of the proposal, and finally, Section 5 concludes the paper.

2. Methods and Materials

2.1. Method

Let Y be a continuous random variable, as slope or orientation, and let Y 1 , Y 2 , , Y k 1 be a set of thresholds defining a partition of the random variable into k disjoint intervals. So, k categories can be defined in the following way: C 1 for those values of Y less than or equal to Y 1 , the categories C i , i = 2 , , k 1 are defined by the interval ( Y i 1 , Y i ] and the category C k for the rest of values of Y (greater than Y k ). In accordance with this categorization, given a random sample of size n from Y, the random vector whose components count the number of events falling into each category, say X = ( X 1 , , X k ) , follows a multinomial distribution with parameters n and P = ( p 1 , , p k ) Δ k = { ( p 1 , , p k ) : p i > 0 , 1 i k , i = 1 k p i = 1 } , where p i stands for the true occurrence probability of the ith category. In short, X M k ( n ; P ) .
Let P 0 = ( p 01 , , p 0 k ) be a fixed point of Δ k which represents the classification adopted on the reference dataset. The hypothetical probabilities p 0 i , i = 1 , , k represent the quality requirements to be achieved by the product, that is, the ideal proportions. Deciding whether or not the quality requirements are fulfilled by a product is tantamount to testing the hypotheses
H 0 c : P = P 0 , H 1 c : P P 0 .
Recall that the critical region of a test is built so that the probability of rejecting the null hypothesis when it is true is bounded by a quantity, usually denoted by α , called the significance level. In other words, when the null hypothesis is rejected either we are making the correct decision or we are wrong, and the probability of the latter is less than α , that is, the probability of a false rejection is controlled. By contrast, when the null hypothesis is accepted, either we are making the correct decision or we are wrong and, in most cases, the probability of being wrong is unknown, that is, the probability of a false acceptance is unknown. Many tests are available for the testing problem (1) (see, e.g., references [17,18]). Because of the above explanation, those tests are applied when H 0 c is expected to be rejected, since when H 0 c is not rejected it cannot be concluded that P = P 0 because the probability of making an incorrect decision is unknown and could be arbitrarily high.
When the interest is centred in proving that P = P 0 , the hypotheses in (1) must be adequately switched, and an equivalence test should be applied, see [19]. With this aim, the user must fix a distance or dissimilarity measure on Δ k , say H, and a positive quantity, say ε > 0 , and now the hypotheses to be tested become
H 0 e : H ( P , P 0 ) ε , H 1 e : H ( P , P 0 ) < ε .
Notice that with this approach, the equality of probability vectors has been replaced by a larger subset of Δ k . Specifically, the subset { P = P 0 } has been enlarged through adding an indifference zone or a neighborhood around P 0 in the parametric space Δ k , H ( P , P 0 ) < ε . In this setting, if H 0 e is rejected then it is concluded that P and P 0 are equal except for practically irrelevant deviations whenever ε > 0 is small enough.
Therefore, when one is interested in proving that a product meets the quality requirements stated by P 0 , a test of (2) must be applied. In [19] a test of (2) is proposed for H the Euclidean distance. Nevertheless, other choices for H are possible. Previous research related to tests for the similarity between spatial point patterns and tests for the thematic-accuracy quality control suggested that the Hellinger distance may be an appealing choice for this goal [20,21,22]. For two multinomial probability distributions, their Hellinger distance is equal to the Euclidean norm of the difference of the square root vectors. To avoid dealing with square roots, from now on we will work with the square of the Hellinger distance, defined as follows:
H ( P , P 0 ) = i = 1 k ( p i p 0 i ) 2 .
Notice that 0 H ( P , Q ) 2 and that H ( P , P 0 ) = 0 if and only if P = P 0 . So, for fixed small ε > 0 , P and P 0 will be considered equivalent according to (3) if H ( P , P 0 ) < ε , implying that both distributions are equal except for insignificant deviations.
To derive a critical region for (2), with H as in (3), we proceed as follows. Let P ^ = ( p ^ 1 , p ^ 2 , , p ^ k ) be the vector of relative frequencies,
p ^ i = X i n , 1 i k .
From Corollary 3.1 in [23], for each P H 0 e ( P P 0 ), one gets
n H ( P ^ , P 0 ) H ( P , P 0 ) L N ( 0 , σ 2 ( P ) ) ,
as n , where
σ 2 ( P ) = 1 i = 1 k p i p 0 i 2 ,
and L denotes convergence in distribution. For σ ^ 2 = σ 2 ( P ^ ) , it follows that
n H ( P ^ , P 0 ) H ( P , P 0 ) σ ^ L N ( 0 , 1 ) ,
as n , for each P H 0 e . So, for α ( 0 , 1 ) , the decision rule is:
r e j e c t   H 0 e i f T n = n ( H ( P ^ , P 0 ) ε ) σ ^ Z α ,
where Z α stands for the α -percentile of the standard normal distribution. The critical point Z α in (4) has been taken from the asymptotic distribution of T n , that is, the level of the test is (exactly) equal to α for infinite sample sizes. To study the behaviour of test (4) for finite sample sizes, we carried out a simulation study, which is described in the next section.
Note that the equivalence test proposed is developed in the context of the multinomial law underlying a categorization scheme. So, it can be also applied not only for classifications of slopes and orientations, but also for other geospatial thematic products (e.g., land covers maps, land uses, maps, etc.).

Simulation Results

For feasibility reasons (e.g., economic, temporal, etc.), quality controls are usually carried out based on sampling. When interest is centred on proving that a product meets the quality requirements stated by P 0 , test (4) can be applied. As observed before, that test is based on asymptotic results. To study numerically its behavior for small and moderate sample sizes, we carried out some simulation experiments.
In these experiments, we considered the testing problem H 0 e vs. H 1 e with P 0 = ( 1 / k , , 1 / k ) for k = 3 , 5 , 7 . This particular choice of P 0 is not casual. It has been considered in previous simulation studies involving inference research over the multinomial distribution (see [24,25,26], among others). The values of k were chosen because they conform to the most usual number of categories in the slope classifications from the literature revised [27,28,29,30,31,32]. Finally, several values for ε were taken, in particular, ε = 0.1 , 0.15 , 0.2 covering reasonable thresholds of closeness between P and P 0 .
The first part of the simulation is devoted to studying the actual level of the equivalence test (4). For this task, as usual, we have to simulate data under the null hypothesis H 0 e , that is, for configurations of P such that H ( P , P 0 ) ε . Here, we considered a set of configurations of P verifying that H ( P , P 0 ) = ε for ε as before. Table 1, Table 2 and Table 3 display the considered configurations. For each scenario, a random sample of size n = 50 was generated and the decision rule (4) was applied for the nominal level α = 0.05 . After 100,000 repetitions, the rejection percentage was collected. The whole experiment was repeated for n = 100 , 200 , 500 for k = 3 , 5 and n = 1000 , 2000 for k = 7 , respectively. Table 4 and Table 5 show the rejection percentages obtained which estimate the actual level of the test for the nominal value α = 0.05 .
Looking at Table 4 and Table 5, we observe that for k = 3 small sample sizes are needed for the actual level to match the nominal level for all considered values of ε . However, when the number of categories increases, larger sample sizes are required for the simulation results to reach the nominal level ( n > 200 for k = 5 and n > 2000 for k = 7 ).
To study the power of the proposal, that is, the probability of rejecting the null hypothesis when it is false, we generated a random sample of size n = 50 from a M ( n , P 0 ) for k = 3 and we applied the decision rule in (4) for α = 0.05 and ε = 0.1 . After 100,000 repetitions, we collected the proportion of rejections which are now the estimated power associated with the nominal level α . The experiment was repeated for ε = 0.15 , 0.2 , n = 100 , 200 and k = 5 , 7 . Table 6 shows the estimated power. As can be seen, in all the tried cases, the power is very high, that is, the procedure makes the correct decision with a probability close to 1. According to these findings, the usefulness of test (4) as a tool for quality control is evidenced.

2.2. Materials

To show an example of an application on real data, a study area around Allo (Navarra, Spain) was considered using the following digital models:
  • SLM02 & ASM02. Both slope (SLM02) and aspect (ASM02) models were derived by GIS operations from a DEM (the DEM02). The DEM02 is a DEM grid with 2 × 2 m resolution, generated in the year 2017 and its primary data source is a Lidar survey (second coverage of the PNOA-LiDAR project https://pnoa.ign.es/estado-del-proyecto-lidar/segunda-cobertura, accessed on 21 May 2021). Both datasets are considered as references in this example.
  • SLM05 & ASM05. Both slope (SLM05) and aspect (ASM05) models were derived by GIS operations from a DEM (the DEM05). The DEM05 is a DEM grid with a resolution of 5 × 5 m, generated in the year 2012 and its primary data source is a Lidar survey (first coverage of the PNOA-LiDAR project https://pnoa.ign.es/estado-del-proyecto-lidar/primera-cobertura, accessed on 21 May 2021).
Both DEM datasets come from the IGN (Spain) and are freely available. In both cases, GIS operations (slope and aspect) are those implemented by ArcGIS (TM) in its Spatial Analyst toolbox.
Concerning the study area, Figure 1 shows its corresponding DEM, and Figure 2 a zoom to a smaller area to have a detailed view to appreciate differences. The area is 504 km 2 , and it has a varied relief, but not abrupt, with valleys of different widths, and areas with different degrees of undulation. The elevation is in the interval 316–1046 m, mean value is 468 m and the standard deviation 92.8 m.
Figure 3 and Figure 4 show the observed values of slopes and aspects, respectively. In general, there are small differences between the curves corresponding to SLM02 vs. SLM05 and ASM02 vs. ASM05, respectively. The differences are minor because they are high-quality products, relatively close in time and for a territory in which there have not been major changes. The curves of the slopes (SLM02 in blue and SLM05 in red) show considerable similarity between the two datasets. Clearly, areas are observed where one curve occurs over the other, and these areas have a certain amplitude. The curves of the orientation (ASM02 in blue and ASM05 in red) show more tremor than the slope curves; this is logical since the orientation is more sensitive than the slope to small changes in elevation. On the other hand, it is observed that the curve corresponding to case ASM02 has more variability than that of case ASM05, which is also logical because it comes from a more detailed model. Furthermore, it is observed that the two curves intersect, although these crosses occur with a relatively wide length between them, which implies a certain bias. It should be remembered that the proposed analysis method is not based on the use of these frequency distributions, which have helped us to better understand the data, but on the comparison of the proportions between categories defined on them and that the differences that are worked with will be those that exist between the proportions corresponding to those classes. In this case, given the statistical control to be applied, the most important thing is to analyze whether, given a categorization of the two variables of interest (slope and orientation) applied under the two models, the hypothetical probability vector for the reference DEM02 ( 2 × 2 m) is close enough the true one of DEM05 ( 5 × 5 m). If this fact is confirmed via an equivalence test, it will mean that the quality conditions of the reference DEM02 are achieved by the product DEM05 to be controlled, except for irrelevant deviations.

3. Results

Given that in the proposed method the number k of categories to be considered is relevant to the sample size, examples with a different number of classes are presented. Thus, for the slope, two analyses with different numbers of categories (k = 3 and k = 5) will be considered, while for the orientation k = 8. In this way, with k = 3,5 and 8, the applicability of the method for the usual numbers of categories is better evidenced. In most practical situations, the probability vector P 0 is obtained from a previous sampling procedure or it can even be proposed by an expert.

3.1. Slope

Over this area, we considered two classification schemes for the slope in the reference data (SLM02). From the reference DEM, such classifications give the percentages of each categories as are summarized in Table 7 and Table 8. Figure 5 and Figure 6 show the application of the classification in Table 7 and Table 8 over SLM02. These values (percentages) are considered the “pattern” ( P 0 ) that the product to be controlled (SLM05) must follow.
Notice that the value of ε is fixed by the user. Here, applying the equivalence test for ε = 0.1 ,
H 0 e : H ( P , P 0 ) 0.1 , H 1 e : H ( P , P 0 ) < 0.1 ,
which led to quality control terms raised as follows:
Null hypothesis. The distance between the distribution of proportions between the classes in the product (SLM05) and the reference (SLM02) exceeds the threshold considered; that is, the product does not achieve the quality requirements stated by the reference data.
Alternative hypothesis. The distance between the distribution of proportions between the classes in the product (SLM05) and the reference (SLM02) does not exceed the threshold considered. In other words, the quality requirements stated by the reference agree with those observed in the product (except for irrelevant deviations).
The equivalence test is applied for random samples of size 100, 200 and 500 drawn from SLM05. Table 9 and Table 10 show the values of the test statistic T n and the result of the decision rules for α = 0.05 for both classifications.
As can be seen in Table 9 and Table 10, in all cases the null hypothesis is rejected at the significance level α = 0.05 , which in this case means that the square Hellinger distance between both sets of DEM data is less than 0.1, so that the observed classifications on the product SLM05 agree with those indicated by the reference SLM02, except for irrelevant deviations.

3.2. Aspect

In this analysis, we considered only one classification scheme for the aspect in the reference data (ASM02). We adopted a classical eight orientation classification based on octants. From the reference model, such classification gives the percentages of each category as are summarized in Table 11 (the category “Flat” does not appear in Table 11 because it is practically null). These values (percentages) are considered the “pattern” of the distribution of class proportions ( P 0 ) that the product to be controlled (ASM05) must follow. Figure 7 shows the categorization in Table 11 over SLM02.
Again, the hypotheses to be tested are
H 0 e : H ( P , P 0 ) 0.1 , H 1 e : H ( P , P 0 ) < 0.1 ,
which led to quality control terms raised as follows:
Null hypothesis. The distance between the distribution of proportions between the classes in the product (ASM05) and the reference (ASM02) exceeds the threshold considered.
Alternative hypothesis. The distance between the distribution of proportions between the classes in the product (ASM05) and the reference (ASM02) does not exceed the threshold considered. To assure the correct application of the equivalence test in this case, we got a random sample of size n = 500 and we obtained the observed proportions given by P ^ = ( 0.070 , 0.152 , 0.124 , 0.134 , 0.148 , 0.134 , 0.098 , 0.088 , 0.052 ) . Considering again ε = 0.1 , it is obtained that T n = 58.30 , which implies the decision to reject the corresponding H 0 e . In other words, the observed classification of orientation in the DEM05 agrees with the percentages of the classification in the reference data DEM02, so that the quality requirements are fulfilled for the orientation classification.

4. Discussion

The discussion will take place along three lines, the first on the method itself, the second on relevant aspects of the proposal as the number of categories and the relationship between the distance and the tolerance ε , and the third on the application to the thematic-type control of slope and aspect classes.
On the first hand, recall the final objective is to decide whether a set of quality requirements stated by the reference dataset is achieved by a product or not. To reach this goal by using a hypothesis test, two points of views can be considered, say, if our interest is centred in proving that the quality requirements are not fulfilled ( P P 0 ) or, on the contrary, if our focus is centred in confirming that the quality requirements are achieved (in the sense of P and P 0 are close enough through the distance H which is bounded by ε ). Since in any hypothesis test the critical region is built so that the probability of rejecting the null hypothesis when it is true is bounded by the significance level ( α ), only the probability of a false rejection is controlled. As a consequence, for the first approach, one must test H 0 c versus H 1 c , whereas for the second one, the hypotheses to be tested must be H 0 e versus H 1 e .
Nevertheless, the difference in the formulation (and hence in the meaning) of the null and alternative hypotheses does not imply any type of change concerning the interpretation of the level of significance and power associated with any hypothesis test.
Secondly, we highlight that the number of categories, the square Hellinger distance and the threshold ε are interrelated in the method and they affect the actual level of the test (2), as is shown in the simulation results. So, in the practical applications, users have to take into account that a larger number of categories requires larger sample sizes and it is suggested to take the smallest number of categories possible. In any case, we consider that the values of k used in the applications cover the most common values of the category numbers.
Finally, the hypothesis test is not carried out on a natural measure to the problem (e.g., degrees, radians, meters, etc.). The square Hellinger distance can be unknown for some users and, although its mathematical basis is understood, it is needed to manage and understand the role of ε . That is, is ε = 0.1 a lot or a little? The answer to this question combines two things: First, the definition of H ( P , P 0 ) itself, which is related in some way to the Euclidean distance on the square roots of the observed frequencies and pattern, and second, the interest that the user puts in what should be the “indifference zone” to decide when both multinomial distributions are close enough to be considered as equal except for “irrelevant” deviations.
In any case, the use of a particular distance measure depending on the type of statistical tools is very common in their applications to geospatial data (e.g., cluster analysis, etc.), where several distances appear (Minkowski, Chebyshev, Manhattan, Mahalanobis, etc.), which are also not “natural” distances to the problem either. This entails the need to achieve sufficient sensitivity in the management of the distance that is applied. In Section 2.1 the study of the level and the power developed can guide the cases that have been considered, but it can also serve as a tool for those interested in the method who want to make their adjustments between the number of categories and ε . Finally, about the execution of these samplings in the field, it should not be forgotten that the measurement of the slope and orientation in the field is usually carried out simultaneously or conditioned by other types of work (e.g., edaphological surveys), so their availability will be at most times conditioned by these types of actions.

5. Conclusions

A statistical procedure was proposed for the quality control of a DEM dataset to meet a reference DEM. The method is based on a categorization scheme over a continuous magnitude and on the application of an equivalence test for multinomial distributions. So, to proceed it is enough to know the values of the proportions for the reference dataset. Since there is a tradeoff between the number of classes to be considered and the size of the sample, a simulation study was carried out which allows guidance concerning the sample sizes. Those interested in applying this method could carry out their simulations to adjust the most suitable size to their case.
On the other hand, the practical example presented illustrates that the application process is quite simple. In addition, and given that statistical tests are used, one could think of carrying out a joint analysis on the slope and the orientation, for this it is enough to apply simple rules such as Holm–Bonferroni correction. The framework of the application of this method is the comparison between the proportions that correspond to the categories of a given classification system of a product and those known from a reference dataset. In this sense, the method was applied to classifications of slopes and orientations but could be applied to compare proportions of any geospatial thematic product (e.g., land covers maps, land uses maps, etc.).
Finally, it is very relevant to indicate that the formulation of the hypothesis test carried out is adequate when interest is centred on proving that a product matches the required quality requirements. The parameter ε that is used must be chosen by users, allowing it to be applied to multiple situations in which greater or lower flexibility is required.

Author Contributions

Conceptualization, M.V.A.-F. and F.J.A.-L.; methodology, M.V.A.-F. and M.D.J.-G.; software, M.V.A.-F. and M.D.J.-G.; formal analysis, M.V.A.-F. and F.J.A.-L.; investigation, M.V.A.-F., M.D.J.-G. and F.J.A.-L.; resources, F.J.A.-L.; data curation, F.J.A.-L.; writing—review and editing, M.V.A.-F.; project administration, F.J.A.-L.; funding acquisition, F.J.A.-L. All authors have read and agreed to the published version of the manuscript.

Funding

M.V.A.-F. and F.J.A.-L. acknowledge financial support from grant PID2019-106195RB-100 (Spanish Ministerio de Ciencia, Innovación y Universidades). M.D.J.-G. acknowledges financial support from grants MTM2017-89422-P (Spanish Ministerio de Economía, Industria y Competitividad, Agencia Estatal de Investigación and European Regional Development Fund), and P18-FR-2369 (Junta de Andalucía).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pike, R.J.; Evans, I.; Hengl, T. Geomorphometry: A Brief Guide; Developments in Soil Science: Amsterdam, The Netherlands, 2009; pp. 3–30. [Google Scholar]
  2. Maune, D.F.; Amar, N. Digital Elevation Model Technologies and Applications. The DEM Users Manual, 3rd ed.; Asprs Publications: Falls Church, VI, USA, 2019. [Google Scholar]
  3. Ariza-López, F.J.; Chicaiza-Mora, E.G.; Mesa-Mingorance, J.L.; Cai, J.; Reinoso-Gordo, J.F. DEMs: An Approach to Users and Uses from the Quality Perspective. Int. J. Spat. Data Infrastruct. Res. 2018, 13, 131–171. [Google Scholar]
  4. Maune, D.F. Digital Elevation Model. Technologies and Applications: The DEM Users Manual, 2nd ed.; American Society for Photogrammetry and Remote Sensing: Bethesda, MD, USA, 2007. [Google Scholar]
  5. Mesa-Mingorance, J.L.; Ariza-López, F.J. Accuracy Assessment of Digital Elevation Models (DEMs): A Critical Review of Practices of the Past Three Decades. Remote Sens. 2020, 12, 2630. [Google Scholar] [CrossRef]
  6. Mesa-Mingorance, J.L.; Chicaiza-Mora, E.G.; Bueaño, X.; Cai, J.; Rodríguez-Pascual, A.F.; Ariza-López, F.J. Analysis of Users and Uses of DEMs in Spain. Int. J. Geogr. Inf. 2017, 6, 406. [Google Scholar] [CrossRef]
  7. Stroeven, P.; Li, K.; Le, N.L.B.; He, H.; Stroeven, M. Capabilities for property assessment on diferent levels of the micro-structure of DEM-simulated cementitious materials. Constr. Build. Mater. 2015, 88, 105–117. [Google Scholar] [CrossRef]
  8. Saksena, S.; Merwade, V. Incorporating the effect of DEM resolution and accuracy for improved flood inundation mapping. J. Hydrol. 2015, 530, 180–194. [Google Scholar] [CrossRef] [Green Version]
  9. Gomez, C.; Hayakawa, Y.; Obanawa, H. A study of Japanese landscapes using structure from motion derived DSMs and DEMs based on historical aerial photographs: New opportunities for vegetation monitoring and diachronic geomorphology. Geomorphology 2015, 242, 11–20. [Google Scholar] [CrossRef] [Green Version]
  10. Rekha, P.N.; Gangadharan, R.; Ravich, R.P.; Mahalakshmi, P.; Panigrahi, A.; Pillai, S.M. Assessment of impact of shrimp farming on coastal groundwater using Geographical Information System based Analytical Hierarchy Process. Aquaculture 2015, 448, 491–506. [Google Scholar] [CrossRef]
  11. Juel, A.; Groom, G.B.; Svenning, J.C.; Ejrnaes, R. Spatial Application of Random Forest models for fine-scale coastal vegetation classification using object based analysis of aerial orthophoto and DEM data. Int. J. Appl. Earth Obs. 2015, 42, 106–114. [Google Scholar] [CrossRef]
  12. Malaperdas, G.D.; Panagiotidis, V.V. The aspects of Aspect: Understanding land exposure and its part in geographic information systems analysis. Energy Environ. 2018, 29, 1022–1037. [Google Scholar] [CrossRef]
  13. Warrena, S.D.; Hohmannb, M.G.; Auerswaldc, K.; Mitasova, H. An evaluation of methods to determine slope using digital elevation data. CATENA 2004, 58, 215–233. [Google Scholar] [CrossRef]
  14. Ryder, W.H.; Voyadgis, D.E. Measuring the Performance of Algorithms for Generating Ground Slope. In Proceedings of the ESRI 1996 International User Conference, Palm Springs, CA, USA, 20–24 May 1996. [Google Scholar]
  15. Van Westen, C.J. Remote sensing and GIS for natural hazards assessment and disaster risk management. In Treatise on Geomorphology; Shroder, J., Bishop, M.P., Eds.; Academic Press: San Diego, CA, USA, 2013; Volume 3, pp. 259–298. [Google Scholar]
  16. Kumar, L.; Skidmore, A.K.; Knowles, E. Modelling topographic variation in solar radiation in a GIS environment. Int. J. Geogr. Inf. Sci. 1997, 11, 475–497. [Google Scholar] [CrossRef]
  17. Pardo, L. Statistical Inference Based on Divergence Measures; Chapman & Hall/CRC: Boca Raton, FL, USA, 2006. [Google Scholar]
  18. Jiménez-Gamero, M.D.; Pino-Mejías, R.; Alba-Fernández, M.V.; Moreno-Rebollo, J.L. Minimum ϕ-divergence estimation in misspecified multinomial models. Comput. Stat. Data Anal. 2011, 55, 3365–3378. [Google Scholar] [CrossRef]
  19. Wellek, S. Testing Statistical Hypotheses of Equivalence and Noninferiority; Chapman & Hall/CRC: Boca Raton, FL, USA, 2010. [Google Scholar]
  20. Alba-Fernández, M.V.; Ariza-López, F.J. A homogeneity test for comparing gridded-spatial-point patterns of human caused fires. Forests 2018, 9, 454. [Google Scholar] [CrossRef] [Green Version]
  21. García-Balboa, J.L.; Alba-Fernández, M.V.; Ariza-López, F.J.; Rodríguez-Avi, J. Analysis of Thematic Similarity Using Confusion Matrices. ISPRS Int. J. Geoinf. 2018, 7, 233. [Google Scholar] [CrossRef] [Green Version]
  22. Alba-Fernández, M.V.; Ariza-López, F.J.; Jiménez-Gamero, M.D.; Rodríguez-Avi, J. On the similarity analysis of spatial patterns. Spat. Stat. 2016, 18, 352–362. [Google Scholar] [CrossRef]
  23. Zografos, K.; Ferentinos, A.; Papaioannou, T. Divergence statistics: Sampling properties and multinomial goodness of fit and divergence tests. Commun. Stat. Theory Method 1990, 19, 1785–1802. [Google Scholar] [CrossRef]
  24. Pardo, M.C. On Burbea-Rao divergence based Goodness-of-fit tests for multinomials models. J. Multivar. Anal. 1999, 69, 65–87. [Google Scholar] [CrossRef] [Green Version]
  25. Alba-Fernández, V.; Jiménez-Gamero, M.D. Bootstrapping divergence statistics for testing homogeneity in multinomial populations. Math. Comput. Simul. 2009, 79, 3375–3384. [Google Scholar] [CrossRef]
  26. Basu, A.; Sarkar, S. On disparity based goodness-of-fit tests for multinomial models. Statist. Probab. Lett. 1994, 19, 307–312. [Google Scholar] [CrossRef]
  27. Sheng, T.C. Watershed Management Field Manual: Watershed Survey and Planning; FAO Conservation Guide 13/6; Food and Agriculture Organization of the United Nations: Rome, Italy, 1990. [Google Scholar]
  28. Elewa, H.; Qaddah, A. Groundwater potentiality mapping in the Sinai Peninsula, Egypt, using remote sensing and GIS-watershed-based modeling. J. Hydrol. 2011, 19, 613–628. [Google Scholar] [CrossRef] [Green Version]
  29. Sajnyg, N.; Krupowicz, W.; Mikulska, K.S. Methodology of validation of agricultural real properties in Poland with the use of geographic information system tools. In Proceedings of the 14th SGEM GeoConference on Informatics, Albena, Bulgaria, 17–26 June 2014; Volume 2, pp. 345–356. [Google Scholar]
  30. Mogaji, K.A.; Lim, H.S.; Abdullah, K. Regional prediction of groundwater potential mapping in a multifaceted geology terrain using GIS-based DempsterShafer model. Arab. J. Geosci. 2015, 8, 3235–3258. [Google Scholar] [CrossRef]
  31. CEN prEN 1998-1—Eurocode 8: Design of structures for earthquake resistance. In Part 1: General Rules, Seismic Actions and Rules for Resistance. Part 1: General Rules, Seismic Actions and Rules for Buildings; Draft No. 6, Doc CEN/TC250/SC8/N335; European Commission: Brussels, Belgium, 2003.
  32. NTC Nuove Norme Tecniche per le Costruzioni. In Part 3: Categorie di Sottosuolo e Condizioni Topografiche; NTC: Rome, Italy, 2008.
Figure 1. DEM over the study area (Allo, Navarra, Spain). Elevation range: [315, 1046] m.
Figure 1. DEM over the study area (Allo, Navarra, Spain). Elevation range: [315, 1046] m.
Remotesensing 13 02069 g001
Figure 2. Zoom to a smaller area (red-edged rectangle marked in Figure 1 for more detailed display purposes) Elevation range: [450, 560] m.
Figure 2. Zoom to a smaller area (red-edged rectangle marked in Figure 1 for more detailed display purposes) Elevation range: [450, 560] m.
Remotesensing 13 02069 g002
Figure 3. Frequency distributions of slopes (SLM02 in blue and SLM05 in red).
Figure 3. Frequency distributions of slopes (SLM02 in blue and SLM05 in red).
Remotesensing 13 02069 g003
Figure 4. Frequency distribution of aspects (ASM02 in blue and ASM05 in red).
Figure 4. Frequency distribution of aspects (ASM02 in blue and ASM05 in red).
Remotesensing 13 02069 g004
Figure 5. Classification of slopes: Case #1 (example in the sub-area shown in Figure 2).
Figure 5. Classification of slopes: Case #1 (example in the sub-area shown in Figure 2).
Remotesensing 13 02069 g005
Figure 6. Classification of slopes: Case #2 (example in the sub-area shown in Figure 2).
Figure 6. Classification of slopes: Case #2 (example in the sub-area shown in Figure 2).
Remotesensing 13 02069 g006
Figure 7. Classification of aspect (example in the sub-area shown in Figure 2).
Figure 7. Classification of aspect (example in the sub-area shown in Figure 2).
Remotesensing 13 02069 g007
Table 1. Probability distributions included in the simulation study for k = 3 .
Table 1. Probability distributions included in the simulation study for k = 3 .
Selected Cases p 1 p 2 p 3
Case 1: ε = 0.100.17580.17580.6483
Case 2: ε = 0.150.14330.14330.7133
Case 3: ε = 0.200.11670.11670.7667
Table 2. Probability distributions included in the simulation study for k = 5 .
Table 2. Probability distributions included in the simulation study for k = 5 .
Selected Cases p 1 p 2 p 3 p 4 p 5
Case 1: ε = 0.100.13500.13500.13500.13500.4600
Case 2: ε = 0.150.14000.15000.14000.15000.4200
Case 3: ε = 0.200.15350.15350.15350.15350.3860
Table 3. Probability distributions included in the simulation study for k = 7 .
Table 3. Probability distributions included in the simulation study for k = 7 .
Selected Cases p 1 p 2 p 3 p 4 p 5 p 6 p 7
Case 1: ε = 0.100.13000.13000.13000.13000.13000.13000.2200
Case 2: ε = 0.150.12290.12290.12290.12290.12290.12290.2629
Case 3: ε = 0.200.11550.11550.11550.11550.11550.11550.3071
Table 4. Rejection probabilities at the nominal level α = 0.05 under the configurations in Table 1 and Table 2.
Table 4. Rejection probabilities at the nominal level α = 0.05 under the configurations in Table 1 and Table 2.
k = 3 k = 5
nCase 1Case 2Case 3Case 1Case 2Case 3
ε 0.100.150.200.100.150.20
500.0570.0470.0520.0350.0370.038
1000.0560.0560.0490.0400.0400.040
2000.0570.0530.0490.0440.0440.042
5000.0540.0520.0500.0460.0460.047
Table 5. Rejection probabilities at the nominal level α = 0.05 under the configurations in Table 3.
Table 5. Rejection probabilities at the nominal level α = 0.05 under the configurations in Table 3.
nCase 1Case 2Case 3
ε 0.100.150.20
500.0190.0220.023
1000.0270.0280.030
2000.0330.0340.033
5000.0390.0400.039
10000.0430.0420.043
20000.0460.0440.044
Table 6. Estimated power.
Table 6. Estimated power.
k = 3 k = 5 k = 7
n \ ε 0.10.150.20.10.150.20.10.150.2
500.9900.9990.9970.9420.9950.9990.8330.9770.995
1000.9980.99910.9980.99910.8990.9970.999
2001111110.99911
Table 7. Classification of slopes: Case #1.
Table 7. Classification of slopes: Case #1.
DegreesClassification P 0
0–1Flat0.07
1–7.5Moderate0.56
>7.5Steeper0.37
Table 8. Classification of slopes: Case #2.
Table 8. Classification of slopes: Case #2.
DegreesClassification P 0
0–2Flat0.17
2–8Undulaing0.48
8–15Rolling0.19
15–30Moderately steep0.13
>30Steep0.03
Table 9. Application of the equivalence test to classification of slopes: Case #1.
Table 9. Application of the equivalence test to classification of slopes: Case #1.
Sample Size P ^ (Observed Proportions) T n Decision
100 ( 0.050 , 0.560 , 0.390 ) −27.89Reject H 0 e
200 ( 0.085 , 0.600 , 0.315 ) −20.83Reject H 0 e
500 ( 0.076 , 0.554 , 0.370 ) −10.26Reject H 0 e
Table 10. Application of the equivalence test to classification of slopes: Case #2.
Table 10. Application of the equivalence test to classification of slopes: Case #2.
Sample Size P ^ (Observed Proportions) T n Decision
100 ( 0.140 , 0.500 , 0.200 , 0.150 , 0.010 ) −10.51Reject H 0
200 ( 0.220 , 0.485 , 0.180 , 0.085 , 0.030 ) −14.21Reject H 0
500 ( 0.150 , 0.500 , 0.200 , 0.128 , 0.022 ) −72.18Reject H 0
Table 11. Classification of aspect.
Table 11. Classification of aspect.
AspectDegrees P 0
North337–360, 0–22.50.11
Northeast22.5–67.50.14
East67.5–112.50.13
Southeast112.5–157.50.13
South157.5–202.50.15
Southwest202.5–247.50.15
West247.5–292.50.10
Northest292.5–337.50.09
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alba-Fernández, M.V.; Ariza-López, F.J.; Jiménez-Gamero, M.D. A New Approach to the Quality Control of Slope and Aspect Classes Derived from Digital Elevation Models. Remote Sens. 2021, 13, 2069. https://doi.org/10.3390/rs13112069

AMA Style

Alba-Fernández MV, Ariza-López FJ, Jiménez-Gamero MD. A New Approach to the Quality Control of Slope and Aspect Classes Derived from Digital Elevation Models. Remote Sensing. 2021; 13(11):2069. https://doi.org/10.3390/rs13112069

Chicago/Turabian Style

Alba-Fernández, M. V., F. J. Ariza-López, and M. D. Jiménez-Gamero. 2021. "A New Approach to the Quality Control of Slope and Aspect Classes Derived from Digital Elevation Models" Remote Sensing 13, no. 11: 2069. https://doi.org/10.3390/rs13112069

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