Altimetry Quality of SRTM and ASTER GDEM Products for Areas with Different Reliefs

Hydroelectric power plants cause impacts that are usually estimated by an indicator, the hydroelectric power generation per hectare of flooded area. And, although entrepreneurs use quality-declared cartographic bases to comply with the standards, at the project stages of a hydroelectric plant it is not a priori determined whether the altimetry tolerance of such bases is sufficient to ensure that the impacts will not increase, which in the project consolida-tion becomes a huge problem. This work aims to define the altimetry quality of SRTM Digital Elevation Models (DEM) and ASTER GDEM and whether they are a priori sufficient to simulate the flood level of hydroelectric power plants in different reliefs. To accomplish this objective, a morphological method of assessment of the DEM quality was developed, through Geographic Information Systems (GIS), so that the altimetry information generated by the models and field surveys, when compared, would show their actual differences also in relation to their areas. To this end, two study areas were used: one with a slightly undulated relief and the other with undulated relief. To validate the models, quality assessments were carried out: based on points, according to the Brazilian Map Accuracy Standard (MAS) and STANAG 2215; and based on surfaces, according to Article 500 of the Brazilian Civil Code and the morphological method. At the end, practical applications relating to the M&P indicator and hydroelectric power plants projects were also carried out. The results presented demonstrate that the quality of an SRTM DEM when used in undulating or gently undulating reliefs can be used up to 1:80,000 scale. Already for DEM ASTER under the same conditions,


Theoretical Framework
The concepts described below serve as the basis for the understanding of the research. Several researches have been developed by scientific community to evaluate the altimetry accuracy of the free distribution DEM, as well as for identification of the possible factors that can influence this altimetry data quality. In the search of the [6] was presented results above expectations, in the search of the [7] was evaluated SRTM for a smooth relief area reporting the smallest errors (Root Mean Square Error (RMSE) = 9.85 m) and, also, for an area with mountainous relief, where the biggest errors were identified (RMSE = 22.96). In the search of the [8], when evaluating an DEM SRTM, verified a quality compatible with the scale 1:250,000.
About the DEM ASTER, some evaluations regarding the PEC performed by [9] showed that the average of the altimetric differences between the DEM coordinates and the reference coordinates was smaller than the USGS value of 20 m and in the [10] was evaluated ASTER GDEM data for areas with different types of relief and soil cover, concluding that the typical urban areas presented the highest altimetry errors (RMSE 21.21 m). In the paper of the [11] was found RMSE of 12.12 m for DEM ASTER and in the [12] was found PEC, Class B, associated to mapping scale 1:100.000 for DEM ASTER.
The consistency of the altimetry errors observed in evaluated DEM for any application is directly related to the objective of the DEM use associated to the mapping scale applied [13]. DEM SRTM and DEM ASTER were evaluated were evaluated by the authors that realized that the values of RMSE increased with increasing slope of the land. In the paper of the [6] and [14] was indicate a Class A Cartographic Accuracy Standard (PEC), associated to 1:100,000 mapping scale, for both DEM SRTM and DEM ASTER, that converged with results by [15].
This research aims to meet a demand related to the application of DEM SRTM and DEM ASTER in hydroelectric plants, because very little on the subject is found in the available literature. In the search of the [16] was verified that DEM SRTM was not adequate to identify places with potential for hydroelectric power generation for the studied region. The main factor that stands out and becomes relevant in this research is the fact that quality assessment was performed through topographic surveys via geometric accuracy leveling, being 5 cm the quality of both planimetric and altimetric reference bases.

Studies and Projects Stages for Implementation of a Hydroelectric Enterprise
According to manuals [17] [18] and [19], in short, the studies and projects stages for implementation of a hydroelectric project are divided into various phases, as follows: Journal of Geographic Information System

Geometrical Elements of Relief
The geometrical elements of a relief, according to the [20], can be classified geometrically according to its inclinations and slopes. Inclination expresses the mean angle formed by the intersection of the vertical slope line and the horizontal axis usually measured from its base. The slope represents the angle inclination expressed in percent ratio between the vertical slope line and the length of the hillslope horizontal line, according to the following equations: (1) 100 S H L = × (2) where: • I is inclination. • S is slope.
• H is the vertical slope line.
• L is the length of the horizontal line. Table 1 shows the ratio values of inclination and slope.
Reference [21] also defined classifications of relief slopes, but as a function of the terrain geomorphology, as follows: Table 1. Inclination to slope ratio [20]. • Flat: slopes lower than 3%.

Cartographic Quality
The product quality must be assessed whenever a document that is being used does not specify its geometrical quality, which is the case of orbital images [22].
The process quality must be assessed throughout the work stages by means of criteria that ensure data accuracy.  [25], and describes the methodology for obtaining planialtimetric data from a sample of 167 points. STANAG considers absolute vertical accuracy as a linear error along Z-axis. The products ratings for the necessary altimetry accuracies are listed in Table 3.  In paragraph 1, it is specified that the difference must exceed 1/20 of the area in order that the difference in dimensions can be formalized. One twentieth of the area corresponds to 5% of the area.
The definition of criteria for choosing the number of points for quality assessment may be performed using a statistical method, specifically the statistical function of Pereira [27]. : where: n = control points. N = base points (vectors). ε = control point error (m). εr = relative control point error.

Spatial Interpolation and the Topo to Raster Method
There are various interpolation tools for diverse purposes, but [28] observed that the ANUDEM method is indicated to generate accurate elevation models, since the erosive effect of water alters considerably the relief contours. The method was specifically designed to create a ground surface that best represents natural drainage and best preserves the ridgelines. The current version of ANUDEM used in ArcGIS it the Topo to Raster.

Radio Detection and Ranging (RADAR) Space Shuttle Topography Mission (SRTM)
The Space Shuttle Radar Topography Mission (SRTM) of the National Aeronautics and Space Administration (NASA) was launched in February 11, 2000 in a period of 11 days. SAR radars in C and X bands were used to obtain topographic data on more than 80 percent of the Earth's emerged area, between 60 degrees North and 56 degrees South [29]. The DEM derived from SRTM of a great part of the planet are now available with a spatial resolution of 30 m. Data can be obtained on http://www-radar.jpl.nasa.gov/srtm [30].
According to [31] and [32], C-band radar data provided by the SRTM mission, with spatial resolution of 90 m, has an absolute vertical and horizontal accuracy of 16 and 20 m, respectively, with 90% of confidence.
Diverse works have been developed to assess the quality of the DEM obtained via SRTM, and in some of them the results were higher than the expected ones, e.g., [6].
Others, such as in [8], the DEM was suitable for utilization only on scale 1:250,000 or lower, and in [16], the DEM from SRTM was not adequate to identify places with potential for hydroelectric power generation for the region studied.

Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM)
ASTER is a cooperation effort between NASA and the Ministry of International Trade and Industry in Japan (MITI), the Earth Resources Data Analysis Center (ERSDAC) and the United States Geological Survey (USGS) to build a global, free access DEM. Since 2009, the DEM built from stereoscopic pairs of images provided by the ASTER instrument have been available at no cost and without restrictions [33].
To obtain surface elevation data, the sensor operates at the stereoscopic mode with 3N (nadir) bands and 3B (backward) bands corresponding to the near infrared range. This product has a 30-m spatial resolution and covers the entire Journal of Geographic Information System globe. It is the most complete land map of Earth, covering 99% of the its surface. The GDEM covers the terrestrial surface between latitudes of 83 degrees North and 83 degrees South, with images of 1 degree per 1 degree [9].
On the USGS website, it can be found the accuracies of ASTER GDEM, the vertical accuracy of the DEM ASTER GDEM being of 20 m for 95% of the points sampled. In the same document, the stated horizontal accuracy of this product is 30 m (1 pixel) also for 95% of the sample [34].
The evaluations related to the MAS that reference [9] conducted show that the mean altimetric differences between the images tested and the reference coordinates were lower than the value shown by USGS.

Materials and Methods
This section presents the characteristics of the study areas, a description of the materials (cartographic products, software and hardware units) used and the method developed in the research.

Stydy Areas
Two Eletrosul-owned HPPs were defined as the study areas. All information relating to the altimeter bases was kindly provided by this company.
The Passo São João HPP (PSJHPP), built in the Ijuí River, between the municipalities of Roque Gonzales, Dezesseis de Novembro, São Luiz Gonzaga, São Pedro do Butiá and Rolador (Northwest of Rio Grande do Sul), is a run-of-theriver power plant with an installed capacity of 77 MW [34]. The plant has a total flooded area of 25.24 km 2 , perimeter of de 181 km and reached 518 urban and rural properties during its implementation. It started operations in 2012. The Flood Level (FL) of the plant is 128.25 m. See a panoramic picture of the PSJHPP in Figure 1.
The São Domingos HPP (SDHPP) built in the Verde River, between the municipalities of Ribas do Rio Pardo and Água Clara (Mato Grosso do Sul), is also a run-of-the-river power plant with an installed capacity of 48 MW [35]. The plant has a total flooded area of 18.64 km 2 , perimeter of de 76.50 km and covered 25 rural properties during implementation. It started operations in 2013. The FL of the plant is 345 m. See a simulated picture of SDHPP ( Figure 2).

Materials
The materials (cartographic products) and equipment units (hardware and software) used are listed below:

Method
In this section, the research phases are described. First, an evaluation of data quality per control point was conducted, using as reference the MAS and the STANAG, and then existing trends were checked by hypothesis testing. Subsequently, the terrains morphology in the study areas was evaluated, and, finally, practical applications were made to the M&P indicator and the hydroelectric project stages.

1) Quality Assessment per Control Point
It was carried out by control points only to verify the altimetry quality of the products tested without considering the relief of the studied regions.
First, the products' declared accuracy was checked to ensure, in the evaluation, control points with accuracy three times better than the one declared by the bases tested. As declared accuracies we considered the definitions contained in item 2.3-MAS planialtimetry Class A standard and STANAG 2215 Class A planialtimetry and Class 0 altimetry applied to the products scales. The plane coordinates derived from reference bases surveys were transformed into coordinates matrices by the DXF2XYZ software, where the number of coordinates of each base was quantified. Each planimetric coordinate from the bases reference received a number in ascending order.
To define the number of control points, the equations 3, 4 and 5 [27] were used to define the number of points to be used as a sample. Approximately 32 (thirty two) control points were indicated for each base. However, as a guarantee, 40 (forty) were used for each study area.
Having considered as homogeneous the errors distribution in all bases tested, a distribution method of control points also homogeneous was adopted. To this end, a simple, unbiased random draw method was utilized, by simple manual draw without replacement, to indicate the 40 vectors (control points) of the reference bases that should be entered into the tested bases.
The 40 reference planimetric coordinates drawn were included into each base to be tested, and manual interpolations were performed using the AutodeskMap software to calculate the new coordinate altitudes. At the points not fitting in between contours of different altitudes, the value of the curve involving the point was defined with no need for interpolation.
To validate the bases tested and verify whether they had the same quality as stated, which was verified a priori, the evaluation method followed the standards defined in item 2.3.1. For this, the tested bases were classified according to the following criteria: • The tested control points utilized must represent 100% of the universe of points; • The coordinates of the control points tested cannot exceed the maximum acceptable error in relation to MAS and the SE identified in the calculation of residuals in number higher than 10% of total.
• The coordinates of the control points tested cannot exceed the maximum acceptable error in relation to the STANAG identified in the calculation of residuals in number higher than 10% of total; • The products that have a number of points with error above the acceptable limit in quantity equal to or below that specified by MAS must be classified into A, B or C cartographic classes; • The products with a number of points with error above the acceptable limit, in an amount above that specified by MAS, cannot have a defined class; • The products with a number of points with error above the acceptable limit, in an amount equal to or below that specified by the STANAG 2215 altimetry must be classified into 0, 1, 2, 3 or 4 standards. Each study area showed 40 planimetric coordinates with two altitudes (one of reference and one of test) for each product tested. The differences of each one were calculated and then compared with MAS and STANAG 2215.
For the difference values, the means and sample standard deviations were calculated, and hypothesis tests were applied first for the mean and T-Student Distribution and, then, by application of Normal Distribution (with attribution Journal of Geographic Information System of the declared SE of the products as population standard deviation). In the bases where systematic trend was found, extraction of the trend was performed by basic arithmetic (sum or subtraction) through maps algebra, and again the quality standards were applied to the bases and new comparisons were performed.

2) Morphological Evaluation
At this stage, the altimetry quality of the products was verified, considering the dimensions formed by the surface covered by the simulated curves in relation to the reference curve for each product and study area tested. For this purpose, several GIS tools were used.
First, the study areas had their reliefs classified in relation to inclination/slope. The calculation and respective classification were performed according to the definitions described in section 2.2. by using the formulations of Equations (1) and (2).
From each DEM, the contour lines that simulated the FL of each HPP were extracted. Extraction of the contours was done automatically through the Contour function of the Spatial Analyst tool of the ArcGIS software. For the bases from which systematic trends were removed, the contours were also extracted, and this stage was named "Extraction of re-simulation curves".
Cartographic generalizations were conducted to eliminate islands and also the contours that exceeded the area demarcated by the FL. The limits defined were the dams of each HPP as well as the end of the flooded area upstream each river.
To check for the quality of the products tested, that is, for the homogeneous distribution of the area formed by the simulated/re-simulated curves in relation to the reference curves, the method of evaluation followed the standards defined in 2.3., "Article 500 of the Brazilian Civil Code". The products were then rated according to the following criteria: • It was considered as given dimension (extension measure usually recorded on the property deed), that one formed by the curves extracted directly from the products tested (DEM) in their respective FL; • It was considered as determined dimension the extension measured on field by geometric leveling and planialtimetric survey of the Reference Bases; • The products that exhibited differences between the dimensions formed by the simulated/re-simulated curves and the reference curve of less than one twentieth of the total stated area, i.e., less than 5%, were approved; • The products that exhibited differences between the dimensions formed by the simulated/re-simulated curves and the reference curve of more than one twentieth of the total stated area, i.e., over 5%, were disapproved; • It was defined as an absolute evaluation method the one which considered only the absolute dimension differences (area and perimeter) of the curves generated by the products tested in the measurement, without considering the location/spatialization of the same in relation to the curve of reference, as if the simulated/re-simulated curves had been produced by any expedient Journal of Geographic Information System method, without considering neighboring boundary lines or any kind of linkage or georeferencing; • It was defined as a relative evaluation method the one which considered, in addition to dimensions (area and perimeter) of the curves generated by the products tested in the measurement, the relative location/spatialization (georeferencing), i.e. the total area that is formed by joining the areas and, in the common area, by the areas intersection, between the simulated/re-simulated curves and the reference ones.
The differences between the simulated/re-simulated areas and the reference ones in the absolute form were calculated using the Autodesk software in the area function. To generate the polygons of the total area and common area of the simulated/re-simulated curves and the reference ones, the geoprocessing tool of the ArcGIS software was used, in the merge and intersect functions, respectively.
To this end, the following relationships were used: where: ∪ = set of union of all spatially-related polygons; ∩ = intersection set of spatially-related common polygons; After that, a morphological evaluation was carried out, which is related to the evaluation of the geometric shape of the polygons formed by the simulated/resimulated curves (contours) and their respective quantitative differences in relation to the polygons formed by the reference curves. To this end, the dimensions formed between the reference curves and the simulated/re-simulated ones, except for the common areas, were calculated by simple logic using the merge, intersect and clip functions of ArcGIS software. The relations of Equations (6) and (7), as well as the following ones, were used: where: IA = inside area formed by the polygons positioned between the reference curve and the simulated/re-simulated curve; OA = outside area formed by the polygons positioned between the simulated/re-simulated curve and the reference one.
V. da Silva Celestino Reginato Journal of Geographic Information System

3) Applications
The results from the simulated/re-simulated curves of the products tested enabled a practical application to two situations usually used in engineering, as defined below.
First, the products tested were classified according to the indicator proposed by [3], herein called M&P indicator, which considered the output, or power generation, of the future HPP and probable flooding area. The values that varied up to more or less than 20% in relation to the calculated reference value were used as parameters.
The second application consisted of fitting the bases into the following projects stages (the "Executive Project" stage, because it is the consolidated stage, will not be part of the assessment because it requires, for implementation, large topographic scales, thus not fitting, a priori into any product tested): • Office stages: Estimation of the Hydroelectric Potential and Hydroelectric Inventory; • Stages that require complementation with field works: Feasibility and Basic Design.
In this work, the optimal curves to be applied in hydroelectric projects were the ones that had larger surfaces in common with those formed by the reference curves, followed by those that, which had no common surfaces exhibited surfaces distribution either inside or outside the reference altitude at a rate of 50% each, to compensate for losses and gains.
In the case of FL simulation/re-simulation by different bases, the areas formed either outside or inside the graphic environment in relation to the reference FL (the altitude effectively defined on field and its consequent filling), must be analyzed carefully because, depending on what it is considered as the best or worst case, from the entrepreneur's point of view or that affected by damming (the owner of the land to be inundated), it may interfere with the results and analyses.
In large plants, if the areas estimated by simulation/re-simulation in the phases of study are larger than the effective one found in the consolidated phases of the project, it could affect the other phases, such as, for instance, the inventory and feasibility study. In these cases, it could and should be recommended to carry out more studies to obtain higher quality information before considering the project unviable.
On the other hand, if the areas estimated by simulation/re-simulation in the phases of study are smaller than the effective one found in the consolidated phases of the project, it could be a problem because a larger flooded area may result in a greater number of affected properties, with larger areas of suppressed vegetation, among others, which will ultimately cause losses to entrepreneurs and an adverse impact on the affected community.

Results
This section contains the results obtained from the evaluation stages per control point and from the morphological evaluation as well as the analyses of the products quality, their matching points and applications.

Quality Evaluation per Control Point
As the SRTM DEM has a reported altimeter accuracy (a priori) of 16 m (according to item 2.5), a value that can be considered as one-half equidistance on contour line, for calculation of the expected accuracy and SE, its compatible scale will be taken into consideration, i.e. 1:80,000. The a priori expected accuracy and SE for the tested SRTM bases are shown in Table 4.
As the ASTER DEM has a reported altimeter accuracy (a priori) of 20 m (according to item 2.6), a value that can be considered as one-half equidistance on the contour, for calculation of expected accuracy and SE, its compatible scale will be taken into consideration, i.e. 1:100,000. The a priori expected accuracy and SE for the tested ASTER bases are shown in Table 5.
Each study area has a specific flood altitude called Flood Level (FL). For the PSJHPP, the FL considered was 128.5 m, and for the SDHPP the FL was 345 m.
The information provided by the field topography generated altimetric bases which, due to their high accuracy, were utilized as reference for the extraction of the control points used in this study.   Figure 4 and Figure 5.
The location of the 40 points for each reference base was analyzed in relation to the contours of each based tested and followed by interpolation. The mean differences of the 40 points of SRTM and ASTER bases for each studied area can be seen in Table 6.
Hypothesis testing was applied to check for systematic trends for the dataset with 39 degrees of freedom and confidence intervals of 90%, 95% and 98%, both to '"t" and "z" statistics, in this case using the values of MAS, Class A Standard, for the products' declared SE as population standard deviation. The calculated values for both statistics for statistics are shown in Table 7.
According to what is shown in Table 7, the hypothesis that the bases could be considered as free from systematic trends was rejected. The 90%, 95% and 97.5% The new difference means for the 40 points of the SRTM and ASTER bases, after removal of systematic trend for each study area, can be seen in Table 8.

1) Evaluation according to MAS
The four highest altitude differences were not considered for calculation of MAS (four points out of 40). The fifth highest value in descending order was considered threshold, i.e., the accuracy value. Based on these criteria, the MAS values obtained are described in Table 9, for the SRTM bases, and in Table 10 for the ASTER bases.     By observing the values described in Table 9, it can be concluded that the SRTM Bases PSJHPP and SDHPP met MAS-Class A standard for altimetry with regard to accuracy and SE. By observing the values described in Table 10, it can be concluded that both ASTER bases met MAS-Class A standard for altimetry with regard to accuracy and SE.
Taking into account the systematic trend shown in Table 6, after the trend removal, the MAS values are obtained, which are shown in Table 11 for the SRTM bases, and in Table 12 for the ASTER bases.
By observing the values of Table 11, it can be concluded that both SRTM bases had an increase in their scale and could then meet MAS-Class A standard for altimetry with better accuracy values when compared with the values shown in Table 9.
By observing the values of Table 12, it can be concluded that the ASTER Base of SDHPP has an increase in its scale and could then meet MAS-Class A Standard for altimetry with better accuracy values when compared with the values shown in Table 10. The ASTER Base of PSJHPP continued to meet the standard with the same accuracy as before.

2) Evaluation according to STANAG 2215
At this stage, the evaluation was also performed by extracting the planialtimetric points in the terrain and the respective inclusion of the horizontal coordinates of these points into the generated products. To fit the product into STANAG, thus, in a given altimetry class, it was verified whether 90% (ninety percent) of the points were within the contour intervals described in Table 3.
In the same way as of item 4.1.a), for calculation of the expected accuracy and SE with regard to the STANAG of SRTM, the scale 1:80,000 and 1:100,000 for ASTER was considered. The a priori expected accuracy for the SRTM and ASTER DEM bases tested are described in Table 13 and Table 14.
Based specifically on the five greatest altitude differences, where the values with results above 10% maximum (four points out of 40) were disregarded, the fifth highest value in descending order is considered threshold, i.e., the accuracy value. Based on these values, the resulting STANAG values are presented in Table 15, for the SRTM bases, and in Table 16 for the ASTER bases. Journal of Geographic Information System     Considering the systematic trend shown in Table 7, after removal of the trend, the STANAG values were obtained, which are described in Table 17 for the SRTM bases and in Table 18 for the ASTER bases.
By observing the values described in

Morphological Evaluation
The relief inclination/slope of each studied area was characterized based on the contour lines of each base tested. To this end, the Equations (1) and (2) were applied, resulting in the slope and inclination values described in Table 19.
According to the classifications developed by [21] and the values shown in Table 19, the relief of the PSJHPP was classified as undulated, with slopes ranging from 8% to 20%. The SDHPP relief, in turn, was classified as gently undulated, with slopes between 3% and 8%.  Table 20.
The simulated dimensions correspond to the curves generated from the tested bases at the FL altitudes.  for SDHPP were also extracted. Table 22 contains the dimensions formed by the re-simulated curves for each area and DEM.
To calculate the absolute dimensions, basic arithmetic was used by first subtracting the values for the simulated bases, of FL and SI, as shown in Table 20 and Table 21. The calculation resulted in the dimensions shown in Table 23.
Afterwards, the same procedure was conducted between the re-simulated values of Table 20 and Table 22. The calculation resulted in the dimensions shown in Table 24.
To calculate the relative dimensions, the Equations (3) and (4) Table 25, for the simulated values, with their respective percentages shown in Table 27.
Afterwards, the same calculations were done for the re-simulated values, which are shown in Table 26, with their respective percentages presented in     According to the dimensions shown in the column FL-SI (m 2 ) Area% of Table 23 and in column TA -CA (m 2 )% of Table 27, relative to the areas, for the bases tested there was a significant difference, much above the acceptable level (5%), and, therefore, the Simulated Curves by SRTM and ASTER, in absolute and relative terms, are not in conformity with Article 500 of the Civil Code. According to the dimensions presented in column FL -SI (m 2 ) Area% of Table 24 and in column TA -CA (m 2 )% of Table 28, relative to the areas, for the re-simulated bases tested, there was a significant difference, much above the acceptable level (5%); therefore, the Re-Simulated Curves by SRTM and ASTER, in absolute and relative terms, are not in conformity with Article 500 of the Civil Code.
When the values of the simulated curves, as presented in Table 27, are compared with the values of re-simulated curves in Table 28, it can be seen that there were significant differences. In most of the re-simulated bases, the differences between the FL and SI diminished, except for the ASTER Base 1, in which the values increased in absolute terms.
By observing the results of Table 27, specifically the values in column CA (m 2 )%, it can be seen that the areas formed in common in the bases tested is relatively small. Only SRTM Base 1 and SRTM Base 2 exhibited common areas slightly above 50% of their total.
By observing the results of Table 28, specifically the values in column CA (m 2 )%, it can be seen that the area formed in all re-simulated bases tested enlarged, when compared with the values presented in the same column of Table  27, exhibiting common areas above 50% in most of them.
At the quantitative morphological evaluation stage, the dimensions formed between the reference and simulated/re-simulated curves were calculated, except for the common areas formed. The results can be seen in Table 29, where in column OA% is shown the percent value relative to the area formed between the simulated/re-simulated curve and the reference curve, considering the total value TA = CA + OA + IA. In column OA/TA% is shown the percent value relative to the area formed between the simulated/re-simulated curve, considering only the areas that are not in common with the simulated/re-simulated curves and reference curves (TA -CA). Column IA% contains the percent value relative to the area formed between the reference curve and the simulated/re-simulated curve, considering the total value TA = CA + OA + IA. In column IA/TA% is shown the percent value relative to the area formed between the reference curve and the simulated/re-simulated curve, considering only the areas that are not in common with the simulated/re-simulated curves and reference curves (TA -CA). The re-simulated bases are indicated by (R) in the "Bases" column.
By observing the results shown in Table 29, it can be seen, especially in columns OA/TA% and IA/TA%, that the majority of the non-common areas of the simulated bases formed were inside the reference curve, with nearly 100% of the areas formed inwardly, as were the cases of the SRTM 2 and ASTER 2 bases. This result shows a diminishing/narrowing tendency of the curve that simulates the inundation and consequent formation of islands.
Therefore, by observing the results relative to the re-simulated bases, a more homogeneous distribution of non-common areas between outside and inside of V. da Silva Celestino Reginato Journal of Geographic Information System Details of each study area and simulated/re-simulated bases can be seen in Also, with regard to the study area 2, but relative to the DEM of ASTER Base 2, the results are justified because the DEM seems not to match perfectly the        For both study areas and tested bases, the inconsistencies in absolute and relative forms in relation to the Civil Code, evaluations already described in Section 4.2, seem compatible with the results obtained in the quality morphological evaluation.

Applications
First, the results were applied to the M&P indicator. The reference values are described in Table 30.
For each tested base, an indicator was calculated based on the extracted curve that simulates/re-simulates the flood. The calculated indicators are described in Table 31, with its respective percentage of 20%.
According to the values presented in Table 31, there was a significant difference in the indicators calculated for both simulated bases when compared with the reference values shown in Table 30, and the simulated curves from SRTM and ASTER were not within the established limit.
However, according to the values presented in the same table for the re-simulated bases (R), there was no significant difference for the calculated indicators The second application was related to the stages of hydroelectric projects, as follows: • As the SRTM Base 1 curve for the simulated height of 130 m virtually does not follow the reference curve and indicates great downward discrepancies and the formation of islands, with only 55.15% of common areas in relation to the total areas, a non-homogeneous distribution of non-common areas (13.23% OA and 86.77% IA), with most of the non-common areas inside the reference curve, as shown in Table 27 and Table 29, it can be concluded that the simulated curve cannot be used to simulate a flood area at any stage of a HPP project; • As the curve of the SRTM Base 1 for the re-simuated height of 135 m, in turn, follows most of the reference curve line and indicates formation of some islands, with 65.05% of common areas in relation to the total and a slightly more homgeneous distribtuion of non-common areas (69.42% OA and 30.58% IA) with most of non-common areas outside the curve, as shown in Table 27 and Table 29, it can be concluded that the re-simulated curve can be used from the Estimation of Hydroeletric Output stage to the Hydroelectric Inventory stage; • As the ASTER Base 1 curve for the simulated altitude of 130 m virtually does not follows the reference curve and indicates great downward discrepancies and formation of inslands, with common areas for only 49.87% of the total, a non-homogeneous distribution of non-common areas (30.13% OA and 69.87% IA), with most of non-common areas within the reference curve, as shown in Table 27 and Table 29, it can be concluded that the simulated curve cannot be used to simulate a flood area at any stage of a HPP project; • As the ASTER Base 1 curve for the re-simulated height of 135 m, in turn, follows most of the reference curve, but indicates formation of some islands, with 53.26% of common areas in relation to the total, and a non-homogeneous distribution of non-common areas (71.25% OA and 28.75% IA) and the majority of non-common areas outside the reference curve, as shown in Table 27 and Table 29, it can be concluded that the re-simulated curve can be used from the Estimation of Hydroelectric Output stage to the first stage, Study Planning, of the Hydroelectric Inventory stage; • As the SRTM Base 2 curve for the simulated altidude of 345 m virtually does not follow the reference curve and indicates great downward discrepancies and islands formation, with only 55.26% of common areas in relation to the Journal of Geographic Information System total, with non-homogeneous distribution of non-common areas (0.80% OA and 99.20% IA) with most of the non-common areas within the reference curve, as shown in Table 27 and Table 29, it can be concluded that the simulated curve cannot be used to simulate a flood area at any stage of a HPP project; • As the SRTM Base 2 curve for the re-simulated height of 350 m, in turn, follows most of the reference curve, both in an out the reference curve, with 79.54% of common areas in relation to the total, with no indication of islands formation, and a slightly more homogeneous distribution of non-common areas (34.40% OA and 65.60% IA), even with most of non-common areas within the reference curve, as shown in Table 27 and Table 29, it can be concluded that the re-simulated curve can be used from the Estimation of Hydroelectric Output stage to the Feasibility stage; • As the ASTER Base 2 curve for the simulated height of 350 m virtually does not follow the reference curve and indicates great downard discrepancies and islands formation, with only 28.93% of common areas in relation to the total, a non-homogeneous distribution of non-common areas (0.06% OA and 99.94% IA) with most of the non-common areas within the reference curve, as shown in Table 27 and Table 29, it can be concluded that the simulated curve cannot be used to simulate the flood area at any stage of the project of a hydroelectric power plant; • As the ASTER Base 2 curve for the re-simulated height of 352 m, in turn, follows partially the reference curve and indicates formation of islands, with 68.78% of common areas in relation to the total and a non-homogeneous distribution of non-common areas (19.69% OA and 80.31% IA), even with most of the non-common areas within the curve, as shown in Table 27 and   Table 29, it can be concluded that the re-simulated curve can be used from the Estimation of Hydroelectric Output stage to the Hydroelectric Inventory stage.

Conclusions and Recommendations
It can be concluded that the main objective was achieved considering that a method based on points was proposed to perform the evaluation of the altimetry quality of altimetric bases and to extract contour lines that simulate the FL, indicating their real differences in areas. The GIS-developed method was applied in two study areas with different reliefs. In bases where altitudes are extracted from SRTM models, after removal of the inherent systematic error found in the product (around −7.3), in slightly undulated and undulated reliefs, quality can be assured for levels between 65% and 79% of the common areas in the flood height simulation.
In bases where the altitudes are extracted from ASTER, after removal of the inherent systematic error found in the product (around −6.2), in slightly undulated and undulated reliefs, quality can be assured for levels between 53% and 68% of common areas in the flood height simulation.
It is a limitation of this study not to have performed the quality assessment by the same point and morphological method study areas that had flat and/or strongly wavy reliefs.
It is recommended to conduct studies using the method described in this paper to assess areas with plane relief, highly undulated and steep reliefs, as well as the use of other products for the same areas tested, e.g., laser scanning, stereoscopic models derived from remote sensing images or radar missions or provided by remotely-piloted air-crafts.
For other applications, such as landslides or inundations, more in-depth studies are recommended.