The use of SRTM DEM data to delineate the SOTER Terrain Unit system for the European states Endre Dobos & Joël Daroussin.
Download ReportTranscript The use of SRTM DEM data to delineate the SOTER Terrain Unit system for the European states Endre Dobos & Joël Daroussin.
The use of SRTM DEM data to delineate the SOTER Terrain Unit system for the European states Endre Dobos & Joël Daroussin SOTER (SOil and TERrain Digital Database) To replace the only existing global soil database at the scale of 1:1 Million Based on existing soil databases Methodology developed between 1988 and 1995 Status of SOTER in 2004 Examples of uses: assessing land degradation processes PROBLEMS…. Inconsistency of the data Interpretation of the map content Outdated variables Edgematching Lack of data Unknown accuracy of the source data SOTER Procedure Modification workshop 20-21October, 2004, Ispra, Italy Agreement on the need to update and revise the procedure published in 1995 Create a guide line (a cooking book) to harmonize the use of DEM to delinate physiogarphic units, which maintain the original meaning and idea of the SOTER methodology. Test the procedure on Europe and derive a SOTER database for Europe Material SRTM Digital Elevation model - 90 m Pilot area Carpathian-basin Pilot area Carpathian-basin SOTER Mapping Units are defined by physiography and lithology Physisography is characterized by 4 differentiating features: 1. Slope intensity 2. Relief intensity 3. Hypsometry (altitude in consideration of the relief intensity) 4. Dissection (PDD: potential drainage density) Objectives To build a quantitative, DEM based procedure for the delineation of SOTER physiographic units at the 1:1 and 1:5 million scales, following and translating the original criteria defined in the SOTER Manual of Procedure. SRTM DEM (90 m resolution) Continuous slope layer Continuous RI layer Continuous hypsometry layer Continuous PDD layer Block majority Reclassified (block size: layers 990 m,(90 grid m)resolution 90 m) BlockResampling majority (block to the size: block990 sizem,(resolution grid resolution 990 m)90 m) Resampling to with the block (resolution 990 m) Fitering (focal majority 4 andsize 6 cells radius circles ) (990 m) Fitering (focal majority with 4 and 6 cells radius circles) (990 m) The SRTM DEM based procedure for deriviring physiographic units Combination of the four layers (990 m) Fitering (focal majority with 3 cells radius circles) Raster based generalization • Elimination of the polygons under the size threshold using the minimum Euclidean distance procedure • Line simplification Vector based generalization Final physiographic units polygon system Regional slope classification Original SOTER 0-2 % Quantitative procedure 0-2 % Gently undulating 2-5 % 2-5 % Undulating 5-8 % 5-8 % Rolling 8-15 % 8-15 % Moderately steep 15-30 % 15-30 % Steep 30-60 % 30-60 % Flat Slopes Resampled slopes Classified slopes Filtered slopes Relief Intensity …median difference between the highest and lowest point within the terrain per specified distance. Units are m/km, m/slope unit, m/2 km Changes of the approach: interpret relief intensity on an aerial basis. 0-50 m/area of a 1 km diameter circle 50-100 m/area of a 1 km diameter circle 100-300 m/area of a 1 km diameter circle 300- m/area of a 1 km diameter circle Continuous relief intensity Classified relief intensity Classified and resampled relief intensity map Hypsometry Level and sloping lands with RI<50m/slope unit Sloping lands with RI>50 m/slope unit Steep and sloping lands with RI>600m/2 km < 300 m 300 - 600 m 600 -1500 m < 200 m 200 - 400 m 400 < …m 600 -1500 m 1500 - 3000 m 3000 - 5000 m 1500 - 3000 m 3000 < ….m Hypsometry class Elevation range (meters above sea level) 1 Up to 10 2 > 10 - 50 3 > 50 - 100 4 > 100 - 200 5 > 200 - 300 6 > 300 - 600 7 > 600 - 1500 8 > 1500 - 3000 9 > 3000 - 5000 10 Above 5000 Hypsometry classes Dissection 0-10, 10-25, and over 25 km/km2 Potential drainage density (PDD) – 0-7 and 8-49 – I. Deriving surface drainage system based on digital elevation data. 1. flowdirection, 2. flowaccumulation, 3. recoding cells having flowaccumulation value higher than a certain - resolution dependent - value are assigned a value of “1” (stream cells) while all others are assigned value “0” (background). – II. A specified sized moving window is sent through the image for counting the stream cells within a given area around each cell and the output value of this function is assigned to the corresponding cell. Potential drainage network Potential drainage network density (PDD) Potential drainage network and classified PDD The PDD class ranges Class 1: less dissected areas, convex surfaces 2: more dissected areas or depressions, concave surfaces PDD value range 0-90 Above 90 The resulting PDD layer Combining the 4 layers Slope Hypsometry Relief intensity to produce… PDD … the raw SOTER physiographic units Vector based generalization procedure Select polygons under the minimum size limit Apply the minimum Euclidean distance method: – Calculate the 4 mean terrain variables for each polygon – Calculate the Euclidean distance for each polygon pairs – Merge each selected small polygon into its neighbour having the smallest Euclidean distance Eliminating small polygons Simplifying lines 0.2 mm separability distance between features on the printout Displacement of the vertices with maximum 200 and 1000 m ground units respectively for the 1:1 and 1:5 million scales The resulting 1:5 million SOTER physiographic units The resulting SOTER physiographic units at 1:1 million scale The resulting SOTER physiographic units at 1:5 million scale Comparing levels of generalization raw 1/1M 1/5M Conclusions There is a strong need for standardized, quantitative data based methodologies to derive soil-landscape information at specified levels of resolution. The 90 m resolution SRTM DEM is a promising data source for 1/1 down to 1/5 million scale soilscape delineation. The procedure can be translated to a purely quantitative manner, maintaining all the advantages of the quantitaive approaches. However, some criteria had to be modified to meet the needs of a DEM-derived approach, like Relief and PDD Thank You for Your attention!