Remote Sensing External Validation of the Aster Gdem2, Gmted2010 and Cgiar-csi-srtm V4.1 Free Access Digital Elevation Models (dems) in Tunisia and Algeria Abstract: Digital Elevation Models (dems) including Advanced Spaceborne Thermal Emission and Reflection Radiometer-global Digital Elevation Mode

Elevation Data 2010 (GMTED2010) are freely available for nearly the entire earth's surface. DEMs that are usually subject to errors need to be evaluated using reference elevation data of higher accuracy. This work was performed to assess the vertical accuracy of the ASTER GDEM version 2, (ASTER GDEM2), the Consultative Group on International Agriculture Research-Consortium for Spatial Information (CGIAR-CSI) SRTM version 4.1 (SRTM v4.1) and the systematic subsample GMTED2010, at their original spatial resolution, using Global Navigation Satellite Systems (GNSS) validation points. Two test sites, the Anaguid Saharan platform in southern Tunisia and the Tebessa basin in north eastern Algeria, were chosen for accuracy assessment of the above mentioned DEMs, based on geostatistical and statistical measurements. Within the geostatistical approach, empirical variograms of each DEM were compared with those of the GPS validation points. Statistical measures were computed from the elevation differences between the DEM pixel value and the corresponding GPS point. For each DEM, a Root Mean Square Error (RMSE) was determined for model validation. In addition, statistical tools such as frequency histograms and Q-Q plots were used to OPEN ACCESS 4601 evaluate error distributions in each DEM. The results indicate that the vertical accuracy of SRTM model is much higher than ASTER GDEM2 and GMTED2010 for both sites. In Anaguid test site, the vertical accuracy of SRTM is estimated 3.6 m (in terms of RMSE) 5.3 m and 4.5 m for the ASTERGDEM2 and GMTED2010 DEMs, respectively. In Tebessa test site, the overall vertical accuracy shows a RMSE of 9.8 m, 8.3 m and 9.6 m for ASTER GDEM 2, SRTM and GMTED2010 DEM, respectively. This work is the first study to report the lower accuracy of ASTER GDEM2 compared to the GMTED2010 data.


Introduction
Digital Elevation Model (DEM) and its derivative attributes (slope, curvature, roughness, local relief, etc.) constitute important parameters for an assessment of any process using digital terrain analysis.These morphometric parameters are commonly used in various applications including flood simulation [1], dune volume calculation [2], volcanic hazards mapping [3], seismic wave propagation [4], soil erosion mapping [5] and relative tectonic activity modeling [6].Summaries of DEM applications in hydrological, geomorphological, and biological applications can be found in Moore et al. [7].
DEMs can be generated using different techniques such as air-borne and satellite-borne stereoscopic photogrammetry, RADAR/SAR interferometry, Light Detection and Ranging (LIDAR), and conventional surveying techniques (e.g., GPS, levelling).These techniques can be compared considering four aspects (i.e., price, accuracy, sampling density, pre-processing requirements).Each technique has its exclusive advantages but also some disadvantages; for a comprehensive review, see [8].However, four main steps are encountered during the generation process of each DEM, regardless of which technology is used [9]: (1) data acquisition (source of elevation data); (2) resampling to required grid spacing (i.e., regular spaced grid); (3) interpolation to extract height of required point (i.e., in between two grid cell centers) and (4) DEM representation, editing and accuracy assessment.All of these steps mentioned above can introduce errors to the final DEM.
DEMs, like other spatial data sets, are subject to errors.Fisher and Tate [10] have investigated errors on gridded data sets and classified them into three main classes: (1) gross errors or blunders; (2) systematic errors and (3) random errors.
Causes and consequences of error in DEM have been widely explored [10].A DEM quality depends on several factors including the acquisition system; methodology and algorithms; complexity of the terrain; grid spacing and data characteristics [11,12].
A variety of DEMs including Advanced Spaceborne Thermal Emission and Reflection Radiometer-Global Digital Elevation Model (ASTER GDEM), Shuttle Radar Topography Mission (SRTM), and Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010) are freely available for the scientific community worldwide.Moreover, several projects are currently being undertaken to generate more precise DEMs.For instance, TanDEM-X mission was an efficient example that would cover the earth at 12 m grid spacing at an accuracy ranging between 4 and 10 m [13].Scientific user can submit, via German Aerospace Center (DLR) home page, a scientific proposal in order to get cost-free DEM extracts.These DEM products may constitute an alternative for researchers in many developing countries with poor geospatial infrastructure [14], including Tunisia and Algeria.However, it is advisable, before attempting to extract different information from these models, to validate their accuracies and to understand the potential and limitations of using these datasets for a specific area.
Techniques for DEM validation have been widely investigated.Gonga-Saholiariliva et al. [15] gave an overview and mentioned various papers related to this topic.One approach of investigation uses the terms of internal and external validation depending on whether or not independent reference data are included in the assessment procedure [16].Another way, often applied, is to group methods of DEM accuracy assessment into quantitative, based on statistics and accuracy measures, and into qualitative based on visual analysis.

Test Sites
Two test sites, representing different elevation ranges and land cover types, were selected in this study: (1) the Anaguid Saharan platform in southern Tunisia and (2) the Tebessa basin in north eastern Algeria (Figure 1).
The Anaguid test site was chosen to represent area in a smooth landscape with elevation differences of approximately 150 m.It is confined within a bounding area rectangle formed by the geographic coordinates 31°57ʹN to 31°44ʹN and 9°50ʹE to 9°33ʹE, covering an area of 280 km 2 .Its topography is flat and denuded; the elevation ranges between 239 m and 383 m with an average of 288 m and standard deviation of 23 m.The average slope, computed from SRTM DEM, varies between values of 4° and 6° (11° at maximum).Hard carbonaceous rocks and conglomerates deposits of the upper Cretaceous to the Neogene outcrop in the Anaguid site [35] The Tebessa test site represents a high Atlas environment with rugged terrain, steep slopes and high elevation differences.It is located in north eastern Algeria within a bounding rectangle defined by the geographic coordinates 35°29ʹ30ʹʹN to 35°16ʹ30ʹʹN and 8°0ʹ0ʹʹE to 8°15ʹ0ʹʹE, covering an area of 784 km 2 .The elevation ranges from approximately 762-1712 m.The average slope, computed from SRTM DEM, is 18°, showing a threefold higher than the Anaguid Saharan platform site.The Tebessa test site includes several different geological formations deposited from the Triassic to the Quaternary, which support different types of landforms such as low rolling hills, and steep mountain slopes [36].Bare ground, grasslands and urban areas are the predominant landcover with forests occupying mountain slopes at approximately 1250 m a.s.l.(Figure 2).The Advanced Spaceborne Thermal Emission and Reflection Radiometer-Global Digital Elevation Model (ASTER GDEM) was developed jointly by the METI (Ministry of Economy, Trade, and Industry) of Japan and the NASA (National Aeronautics and Space Administration).It was released to the public and made freely available on 29 June 2009 [36].The ASTER GDEM product is generated from automatic processing of 1.5 million stereo pairs by applying the stereo correlation methodology.The ASTER GDEM is provided at a one arc-second resolution (approximately 30 m) and referenced to the World Geodetic System 1984 (WGS84).Elevations are computed with respect to the Earth Geopotential Model 1996 (EGM96) geoid.The absolute vertical accuracy of ASTER GDEM is 20 m at 95% confidence level.In 2011, a second data product (ASTER GDEM2) was introduced in which the original data were: (1) augmented with an additional 260,000 overlapping images, and (2) refined in order to decrease the incidence of data artifacts, improve the spatial resolution, and increase the accuracy of water body coverage [37].The resulting DEM is available online through the Data Pool at the NASA Land Processes Distributed Active Archive Center (LP DAAC) [37], from which tiles No. 31E009 and 34E008 covering the Anaguid and Tebessa test sites respectively, were downloaded.
Shuttle Radar Topography Mission (SRTM) DEM represents a collaborative effort between the National Imagery and Mapping Agency (NIMA) and NASA of U.S. Data were collected over an 11-day mission in February 2000 and are available for approximately 80% of the globe (up to 60° north and 56° south).Elevations were measured via radar interferometry, using X-band and C-band Interferometric Synthetic Aperture Radar (InSAR) sensor (5.6 cm and 5.3 cm wavelength, respectively).SRTM C-band data are available at 1 arc-second for USA and 3 arc-seconds outside of the USA.The absolute vertical height accuracy is 16 m and the absolute horizontal accuracy is 20 m.Since its original release, SRTM has been updated in order to remove erroneous pixel values, better delineate water body boundaries, fill data voids, etc. [38].The most recent SRTM product has been treated by Consultative Group for International Agriculture Research Consortium for Spatial Information (CGIAR-CSI) and is available online through the CGIAR-CSI download database [38], from which tiles No. 3609 and 3605 were downloaded.

Reference Elevation Data
The accuracy assessment of the DEMs requires a large number of checkpoints with high accuracy to obtain reliable measures.According to Maune [40], the accuracy of the checkpoints should be at least three times more accurate than the DEM elevations.The American Society of Photogrammetry and Remote Sensing (ASPRS) recommended a minimum of 20 checkpoints in each major land cover category.In the case of three landcover classes (e.g., open terrain, forested areas, and urban areas), a minimum of 60 checkpoints are required [41].For the purpose of this study, the check points were collected using Global Navigation Satellite Systems (GNSS) survey techniques.For Anaguid test site (Figure 3b), GNSS data were selected from an intensive 3D seismic survey conducted on February 2011.The acquisition technique is the Real Time Kinematic Differential Global Positioning System (RTK-DGPS), consisting of 328 GNSS points measured by Leica GRX 1200 series receivers.This moving mode provides a precision and accuracy of 20 mm + 1 ppm vertically and 10 mm + 1 ppm horizontally [42].As for Tebessa test site (Figure 4b), the GNSS validation points are surveyed using Trimble GeoXH 6000 series receivers.They represented 60 ground control points (GCPs), measured using rapid static GNSS surveying techniques at decimeter-centimeter level accuracy [43].However, the accuracy of GNSS data itself depended on many factors including the positional dilution of precision (PDOP) and the number of satellites in view (SV).Lower PDOP value indicated precise GPS data [44].Those results were confirmed by Kaya and Saritas [45] who classified PDOP as (1) excellent if PDOP ranges between 1 and 2 and (2) good if PDOP varies between 2 and 5. Our results indicated a PDOP average value of 2.14 and 1.95 for Anaguid and Tebessa test sites, respectively.The average SV value was 7.82 and 9.25 for Anaguid and Tebessa, respectively (Figure 5).The Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010) was produced by the U.S. Geological Survey (USGS) and the National Geospatial-Intelligence Agency (NGA) and became available in 2010 [39].The mission's objective was the development of a global elevation model that replaces GTOPO 30 arc-seconds data (GTOPO30).GMTED2010 was derived from 11 raster based elevation sources using DEMs fusion technique.The GMTED2010 products were generated using the following aggregation methods: minimum elevation, maximum elevation, mean elevation, median elevation, standard deviation of elevation, systematic subsample, and breakline emphasis.Each product was generated at three separate resolutions of 30 arc-seconds, 15 arc-seconds, and 7.5 arc-seconds (approximately 1 km, 500 m and 250 m, respectively).All products are in a geographic coordinate system referenced to the WGS84 horizontal datum, with the horizontal coordinates expressed in decimal degrees.The vertical units for the elevation values are integer meters, referenced to the EGM96 geoid as the vertical datum.In the present study, the systematic subsample product was selected because of its lower RMSE (ranging between 26 and 29 m) in comparison with other GMTED2010 products [39].

Data Preparation and Processing
After downloading the tiles of the ASTER GDEM2, SRTM v4.1 and GMTED2010, the three DEMs are projected to Universal Transverse Mercator (UTM) zone 32, and then clipped to the extent of each test site.As mentioned earlier, the three DEMs were given in terms of orthometric heights with respect to the geoid model computed from the Earth Geopotential Model 1996 (EGM96).However, the GNSS used Word Geodetic System 1984 (WGS84) as default vertical datum to compute the relative heights [45].Therefore, GNSS geometric height had to be converted into orthometric height by subtracting the height of the geoid at each GNSS point location.For this reason, MS-DOS open source software developed by NASA Goddard Space Flight Center, the NIMA, and Ohio State University is used.This software computed geoid heights over land area with accuracy better than 1 m, except in areas lacking accurate surface gravity data [46].
GNSS data are loaded into the ArcGIS software as vector points [47].The -Extract Values to Points tool‖ in ArcGIS is used to extract the elevation value from the DEM at each control point location.The open source software R [48] is then used for model validation.

Validation Methods
In this study, the validation process was based on geostatistical and statistical measures.From a geostatistical point of view, the vertical accuracy of the three DEMs was evaluated by comparing the variogram of each DEM with that of GPS validation points.This step involved the computation of the semivariance between two sample points with respect to the distance that separates them.The gstat package for the statistics software R was used to produce the empirical variograms [49].Mathematically, empirical variogram can be computed using the following equation ( [50], Equation ( 1) In Equation ( 1) N(h) is the set of all pairwise Euclidean distances i − j = h, |N(h)|is the number of distinct pairs in N(h), and z i and z j are data values at spatial locations i and j, respectively.
From a statistical point of view, the vertical accuracy of the three DEM was assessed by computing the differences between the DEM pixel value and the corresponding GPS point.For each point, an elevation error was computed as the difference between explored and reference data (Equation ( 2)): In Equation ( 2), Z dif is the elevation error, Z ex is the elevation of the explored DEM, and Z ref is the elevation of the GPS points.Positive differences represent locations where the DEM elevation exceeded the GPS point elevation; and, conversely, negative errors occur at locations where the DEM elevation was below the GPS elevation.After that, the mean error (ME), standard deviation Error (STD), root mean square error (RMSE) and maximum (Max) and minimum (Min) error values were calculated as follows: STD and RMSE are measures of surface quality and provide insight into the distribution of deviations on either side of the mean value.The level of agreement between ASTER GDEM2, SRTM and GMTED2010 derived elevation values and GPS data is also evaluated in terms of linear regression and correlation.
In order to describe and compare the error distributions in each DEM, a normality test is made.Quantile-quantile plots (Q-Q plots) are generated for visual inspection of normality.The Q-Q plot is a scatter plot with the quantiles of the observed values on the horizontal axis and the expected normal value on the vertical axis.A dataset with a best-fit linear relationship indicated that the observed values were normally distributed [51].It is also possible to use statistical tests to investigate whether data originate from a normal distribution, but these tests are often rather sensitive in case of large data sets or outliers [52].Therefore, visual methods are preferred.More details about normality test can be found in [53], who recommended visual methods as a component of good data analysis.

Results and Discussions
Figure 6 shows the empirical variograms of the elevation computed for each test site and for each elevation data source.It can easily be observed that all variograms showed similar trend, but with different sills.The SRTM variogram is much closer to the GNSS variogram than those of ASTER and GMTED2010 for both sites.The ASTER model presents a much lower sill variogram than other elevation data regardless of the relief zone.Therefore, the plotted variograms indicate the superior vertical accuracy of SRTM v4.1 followed by GMTED2010 and ASTER GDEM2 compared with the GNSS elevations data.The correlation plots between GNSS data and each of the three DEMs obtained for both sites are shown in Figures 7 and 8a-c.As mentioned earlier, these plots are based on a selection of 328 and 60 GNSS points for Anaguid and Tebessa test sites, respectively.Results from site 1 (Figure 7a-c) indicate that SRTM v4.1 and GMTED2010 elevation values are slightly better correlated to the reference than ASTER.The correlation coefficients are 0.96, 0.99 and 0.97 for ASTER, SRTM v4.1 and GMTED2010, respectively.The plot for site 2 (Figure 8a-c) indicates, however, that the three DEMs have the same correlation coefficient with the reference elevation data.All stations extend the line of perfect fit, showing excellent coefficient of correlation (r = 0.99).This could be due to the number and distribution of selected points.The histograms of the elevation differences between GNSS data and the corresponding DEM are presented in Figures 3 and 4d-f.The Table 2 lists the minimum, maximum and mean error together with the RMSE and standard deviation values.The basic statistics of errors in ASTER show a negative mean error of −4.07 m and −1.8 m (i.e., biased negatively) for Anaguid and Tebessa test sites, respectively.The histograms of elevation differences (Figures 7d and 8d) present a slightly negative skew, indicating that the ASTER model underestimates the spatial distribution of terrain elevation, more pronounced on flat and less complex terrains (site 1) than in hilly and complex terrains (site 2).This underestimation of ASTER was noted in previous studies.Hirt et al. [23] reported a clear negative bias for ASTER with respect to GCPs.
SRTM errors show positive mean values of 2.9 m for the Anaguid and 0.48 m for the Tebessa test sites.The histograms of the differences between SRTM DEMs and the GNSS reference (Figures 7e and 8e) closely follow a normal distribution.However, there is a slight bias towards positive values.The elevation value in the SRTM DEMs is greater than the GNSS value.This can be explained by the difference between the surface model provided by the SRTM measurements and the ground elevations measured by the GNSS survey.These observations were further confirmed in many previous studies [27,33], but debated by Li et al. [28] and Zhao et al. [25] who did not find a positive bias for SRTM.
The statistical values of errors in GMTED2010 model present also a positive mean error of 3.4 m and 1.12 m for Anaguid and Tebessa test sites, respectively.The histograms (Figures 7f and 8f) present frequencies of positive errors greater than those of the negative errors.A large positive bias is clearly observed in both histograms, indicating that the GMTED2010 overestimated the terrain elevation.A general overestimation of terrain elevation by the systematic subsample GMTED2010 with respect to 1.6 million control points was reported by Danielson and Gesch [39].
Figure 9 shows the Q-Q plots of elevation errors in three datasets for each test site.These graphics were generated using -qqplot‖ commend in the optional -car‖ R add-on package [52].The default options of the qqplot procedure automatically generated a reference line, along with upper and lower 95% confidence intervals for the plot.The Q-Q plots for the Tebessa test site (Figure 9d-f) show best-fit linear relationship with the exception of few points on both ends of the distribution.However, the Q-Q plots for site 1 (Figure 9a-c) expose a sigmoid-type function with a considerable deviation from the straight line, indicating that the data are not normally distributed.Many observations show a strong deviation from the 95% confidence limits, suggesting presence of outliers and gross errors affecting the data.A simple approach based on Li et al. [28] was applied in this case study.It consists to use an 3σ elimination empirical rule [28]: (1) DEM was corrected by shifting the bias between DEM and GPS (i.e., (1) DEMcorrected = DEM − Mean (DEM-GPS); (2) the standard deviation between the corrected DEM and GPS was calculated; (3) all pixels with differences between the corrected DEM and GPS more than three times the standard deviation are considered to be gross errors, therefore removed.
After eliminating outliers, the comparison of the accuracy of all DEMs for both sites was recomputed.The 3σ elimination empirical approach greatly increased the accuracy of all DEMs, as indicated by all validation statistics reported in Table 2.For Anaguid test site, the RMSE value is improved by 3.2 m for ASTER, 1.9 m for SRTM and 1.7 m for GMTED2010.The SRTM model continues to yield a better accuracy (RMSE = 3.6 m) than GMTED2010 (RMSE = 4.5 m) and ASTER (RMSE = 5.3 m).Both SRTM and GMTED2010 models present a mean error near to zero and an RMSE approximately equal to the standard deviation, indicating that the estimator for these datasets is unbiased.For Tebessa test site, the RMSE measured for elevation differences between SRTM-GNSS is 8.3 m, while this value for ASTER-GNSS is 9.8 m and 9.6 m for GMTED2010-GNSS.Consequently, the SRTM DEM is more accurate than GMTED2010 and ASTER.
The vertical accuracy assessment of the three DEMs reveals that the SRTM data shows a better vertical accuracy than both ASTER GDEM2 and GMTED2010 for both sites.Surprisingly, the vertical accuracy of ASTER GDEM2 is less than the accuracy of GMTED2010, as it achieved a much higher RMSE compared to GMTED2010.The better accuracy of SRTM over ASTER has been noted in previous studies.Those results were contradicted by Mukherjee et al. [30] who mentioned that ASTER gives better accuracy (Table 3).However, the large variation in average global height precision found in the examined DEMs literature seems to indicate that the vertical precision of the DEMs depends considerably on several parameters (i.e., location, reference data errors, terrain characteristics and surface feature properties).Therefore, it is recommended to perform an in depth investigations about those factors.In addition, a quick look at Table 3 reveals that the RMSE values of 5.3-3.6 m for SRTM v4.1 and ASTER for Anaguid site are very similar to those reported by Rexer and Hirt [33] and Suwandana et al. [27], but much lower than those of Li et al. [28].However, this study is the first investigation to report the lower accuracy of ASTER compared to GMTED2010 data.An explanation for this finding can be attributed to the accuracy of the original data source used to produce GMTED2010.According to Danielson and Gesch [36], all of the 7.5-arc-second GMTED2010 products derived from the 1-arc-second SRTM DTED2 void-filled data, like our study areas, have an RMSE around of 9.7 m that was stated for the 1-arc-second SRTM mission specification.However, many authors have attributed the cause of the relatively low accuracy of ASTER to the varied stacking numbers from place to place.The stacking number (NUM files) is included in the quality assessment (QA) file and corresponds to the number of stacked ASTER DEM scenes used to compute the final ASTER elevation value for the corresponding pixel in the DEM file [21].Theoretically, the ASTER RMSE should decrease with increasing stack number.Miliaresis and Paraschou [25] have found that the vertical accuracy of ASTER GDEM2, relatively, is better with increasing stack number and increasing number of grid points used per NUM class.According to ASTER GDEM2 validation team [34], an increase in the number of scenes reduces error significantly between 1 and 10 scenes, but there is little improvement after about 15 scenes.

Conclusions
This study investigated the quality of three new digital elevation models ASTER GDEM2, CGIAR-CSI SRTM v4.1 and GDEMTED2010 over two different geographical localities in Tunisia and Algeria, all of which are available free of charge.First, the basic characteristics of the models were described.Then, comparisons among the three models were presented and their respective vertical accuracy was estimated by means of comparisons against GNSS reference data.Finally, model differences were discussed from statistical viewpoint.For the assessment of the accuracy, statistical methods based on frequency histograms and Q-Q plots were presented and the error characteristics in three datasets for each test site were investigated.After outliers' removal using a 3σ elimination rule, the validation statistics were recomputed for each DEM and for both sites.The results for site 1 reveal a RMSE value of about 5.3 m of ASTER GDEM2 vs. GNSS data.The RMSE value is higher than that computed for SRTM (3.6 m) and GMTED2010 (4.5 m).The results for site 2 show RMSE of terrain elevation ranging from 8.3 m for the SRTM v4.1 to 9.8 m for ASTER GDEM2.A RMSE of 9.6 m was calculated for the recently available GMTED2010.For both sites, the results suggest that ASTER GDEM2 slightly underestimates altitudes, while this trend was not found in SRTM or GMTED2010 DEM.The significant differences between the two test sites can probably be attributed to vegetation coverage.Tebessa is more vegetated (brushes and trees) than the Anaguid area.Both sensors, ASTER and SRTM are sensitive to the top of the canopy.That is why the negative bias of ASTER is not as severe (only −1.8 m) and the RMSE values of SRTM, GMTED2010 and ASTER in Tebessa are higher, in general.In conclusion, this study has demonstrated the importance of computing validation statistics for DEM before and after removing outliers.Further investigation can be made to evaluate the outliers and reasons for their occurrence.

Figure 1 .
Figure 1.Location of the test sites plotted on a shaded relief map obtained from SRTM DEMs.

Figure 2 .
Figure 2. Landscape characteristics of the (a) Anaguid Saharan platform and (b) Tebessa Basin.Note the GPS base station photo in (a).

Figure 5 .
Figure 5. Scatter plots of satellites in view values vs. PDOP for (a) Anaguid and (b) Tebessa test sites.

Figure 6 .
Figure 6.Empirical variograms of the three DEMs compared with that of GNSS elevations points for each test site: (a) Anaguid and (b) Tebessa.The red, black, blue and green variograms represent GNSS, SRTM v4.1, GMTED2010 and ASTER GDEM2, respectively.

Figure 7 .
Figure 7. Summary statistics for Anaguid test site.Scatter plots of GNSS elevations vs. (a) ASTER GDEM2, (b) SRTM v4.1 and (c) GMTED2010.The fine red line stands for the line of perfect fit.Histograms of elevation errors and relevant descriptive statistics.(d) ASTER GDEM2 minus GNSS elevations, (e) SRTM v4.1 minus GNSS elevations and (f) GMTED2010 minus GNSS elevations.The solid red line represents the fitted density curve.

Figure 8 .
Figure 8. Summary statistics for Tebessa test site.Scatter plots of GNSS elevations vs. (a) ASTER GDEM2, (b), SRTM v4.1 and (c) GMTED2010.The red line stands for the line of perfect fit.Histograms of elevation errors and relevant descriptive statistics: (d) ASTER GDEM2 minus GNSS elevations, (e) SRTM v4.1 minus GNSS elevations and (f) GMTED2010 minus GNSS elevations.The solid red line represents the fitted density curve.

Figure 9 .
Figure 9. Q-Q plots for the Anaguid test site (a, b and c) showing the error distribution for (a) ASTER GDEM2, (b) SRTM v4.1 and (c) GMTED2010.Q-Q plots for the Tebessa test site showing the error distribution for (d) ASTER GDEM2, (e) SRTM v4.1 and (f) GMTED2010.The solid and dashed red lines represent theoretical normal distribution and 95% confidence intervals, respectively.

Table 1 .
Input source data characteristics.

Table 2 .
Difference statistics of the study sites (Units in meters).

Table 3 .
Varying reported height accuracies for the ASTER GDEM2 and SRTM v4.1 DEMs (unit in meters).