Blog

Ebergötzen is 10x10 km study area in the vicinity of the city of Göttingen in Central Germany (51°30’03.16’’–51°35’31.45’‘N; 10°00’28.67’’–10°09’15.21’‘E). This area has been extensively surveyed over the years, mainly for the purposes of developing operational digital soil mapping techniques. The dataset has also been frequently used by the SAGA development team and the SciLands GmbH in many of their demonstrations and documents.

Courtesy of Gehrt Ernst, the State Authority for Mining, Energy and Geology, Hannover, Germany.

The final Google Earth layout showing predicted soil texture fractions in topsoil.

Available layers:

- points.dbf - the point dataset consists of lab measurements four variables are available: SAND, SILT and CLAY (all expressed as % of mass measured for the 0-30 cm layer of soil) and SOILTYPE (type of soil based on the German classification system).
- DEM25.asc - 25 m DEM derived from the topo-maps;
- DEM100.asc - 100 m SRTMDEM;
- landimg.lan - LANDSAT image bands obtained from the http://image2000.jrc.it Corine Land Cover 2000 Project. The image consists of seven bands.
- ZONES.asc - 1:50.000 geological map of Germany.

Grid definition:

ncols: 400
nrows: 400
xllcorner: 3570000
yllcorner: 5708000
cellsize: 25 mproj4:+init=epsg:31467

Lineage: All input raster maps are in ArcInfo *.asc format, and the point data (tables) are in a *.dbf format. All coordinates are in the official German coordinate system, zone 3 (germany3): Transverse Mercator Projection, central meridian is 9°, false easting 3500000, Bessel 1841 ellipsoid with Potsdam datum. The bounding coordinates of the study area are: XMIN=3570000, YMIN=5708000, XMAX=3580000, YMAX=5718000. The input raster maps are available in two grid resolutions: 25 m (fine) and 100 m (coarse). The sand, silt and clay values have been determined by using the so-called_texture by hand_method: a surveyor distinguishes to which of the 32 texture classes a soil samples belongs to, and then estimates the content of fractions. E.g. texture classSt2has 10% clay, 25% silt and 65% sand.

Data owner: State Authority for Mining, Energy and Geology, Hannover, Germany
Reference: Gehrt, E., Buhner, J., (2001) Vom punkt zur flache — probleme des ‘upscaling’ in der bodenkartierung. In: Diskussionsforum Bodenwissenschaften: Vom Bohrstock zum Bildschirm. FH, Osnabruck, pp. 17-34.

Location:Ebergotzen, Germany
51° 34’ 19.2324” N,10° 6’ 28.7964” E
See map:Google Maps

Attachment:

ebergotzen.zip

ebergotzen_input.zip

Short title:BARSOIL

Inputs:baranja.txt - 59 field profile observations; 15 raster maps (predictors);
Outputs:Predictions and simulations of SOLUM and GLEY_P.

Purpose and use:

Script to interpolate soil thickness (SOLUM) and occurrence of gley horizon (GLEY_P) using land-surface parameters and mapping units. Prepared for the needs of book ‘Geomorphometry: concepts, software, applications’.

Programming environment:R / S language
Status of work:Public Domain

Reference: Geomorphometry: Concepts, Software, Applications
Data set name: Baranja hill

Attachment:

BARSOIL.zip

Start: Aug 31 2009 - 09:15End: Sep 2 2009 17:30
Timezone: Europe/Amsterdam

Programme

Monday 31.08.2009

   Tuesday 1.09.2009

   Wednesday 2.09.2009

The case study area ‘Zlatibor’ is located in the South-western part of Serbia (centred at 43°43’44.6’‘N and 19°42’37.8’‘E). The area is mainly hilly plateau, with the exception of the north-east part where the slopes are much steeper. Elevations range from 850 m to a maximum of 1174 m; the total size of the area is 13.5 square kilometers.

Study area Zlatibor: (a) perspective view on the area (1:25,000 topo-map) and location of 1020 error assessment points, (b) a preview of the auxiliary predictors used for the geostatistical modeling

Available layers:

- elevations.txt - a set of 2051 height measurements used for generation of DEMs;
- control.txt - a set of 1020 very precise spot heights used for error assessment;
- dem30.asc - the original topo-map DEM at 30 m resolution;
- SRTMDEM.asc - 3 arcsec (90 m) SRTM DEM;

Grid definition:

ncols: 150
nrows: 100
xllcorner: 7394249
yllcorner: 4841999
cellsize: 30 m

proj4:+proj=tmerc +lat_0=0 +lon_0=21 +k=0.9999 +x_0=7500000 +y_0=0 +ellps=bessel +towgs84=574.027,170.175,401.545,4.88786,-0.66524,-13.24673,0.99999311067 +units=m

Lineage:

The original topo-map DEM was produced by digitizing contour layers from two adjacent sheets of the 1:5000 topographic maps with contour interval of 5 m. Two sheets were scanned by ANATech Evolution scanner with 400 DPI resolution, then georeferenced to the Gauss-Krüger coordinate system (7th zone) and converted to a point map using a semi-automated digitalization of contour lines (I/GEOVEC, Intergraph program module). The faults obtained during automated digitalization were removed by 3D editing of contour lines. The final DEM was produced using the ArcGIS 3D analyst: first a TIN was produced, which was then converted to a regular grid of 30~m resolution. This will be referred to as the topo-DEM in further text. The original points extracted from the topo-map (51,847 points) were sub-sampled (for computational efficiency) to 2051 points.

A set of 1020 photogrammetric control points was provided with the help of the Geodetic governmental authority of Serbia. These were obtained throughout the orthophoto map production of the scale 1:1000 for local municipality. Aerial images were obtained using the RMK 21/23 analog camera with calibrated focal length f=207.96 mm. The average flying altitude was 1040 m. A stereoscopic model was produced using the WILD A10 analog stereo-restitution instrument and MapSoft2000 program package. The land-surface points were measured manually from the stereoscopic model with average lag of 15 m. This gave a total number of 46,021 points that were sub-sampled to 1020 points for faster processing. The estimated height accuracy of control points is 15 cm which allows us to use it as ground truth for the topo-DEM.

Data owner: Geodetic governmental authority of Serbia
Reference: Geostatistical modelling of topography using auxiliary maps
Location: Zlatibor, Serbia
43° 43’ 0.0012” N,19° 41’ 60” E See map: Google Maps

Attachment:

Download

The case study “fishcamp” is of size 1x2 km; located at 37.46353 N; 119.6119 W. The coordinate system used is the UTM NAD83 zone 11 North i.e. EPSG:26911. The complete dataset was obtained from the USGS National Map seamless server.

Fishcamp data layers.

• 2.5 m LiDAR DEM

• 1 m color orthoimage

• Topo-map 1:25k (contour lines)

• 10 and 25 m DEMs

Available layers:

- DEM2m.asc - 2.5 m DEM derived from the LiDAR survey (ground reflectance);
- DEMNED03.asc - 10 m National Elevation Dataset;
- DEMSRTM1.asc - 1 arcsec SRTM DEM (finilized);
- lidar.shp - subsampled LiDAR point measurements (the original dataset consist of over 5 milion of points);
- orthoimg.lan - National Agriculture Imagery Program (NAIP) Orthoimagery for Zone 11;
- topo24k.lan - 1:24k topo-map “White Chief Mountain, CA (37119-D5-TF-024);
- contours.shp - contours digitized from the topo24k map; - tstreams.shp - contour lines digitized from the topo map;
- soilmu.asc - 5 m gridded soil map for the area (1=”HOLLAND FAMILY, 35 TO 65 PERCENT SLOPES”, 2=”CHAIX-CHAWANAKEE FAMILIES-ROCK OUTCROP COMPLEX”, 3=”CHAIX FAMILY, DEEP, 5 TO 25 PERCENT SLOPES”, 4=”CHAIX FAMILY, DEEP, 15 TO 45 PERCENT SLOPES”, 5=”HOLLAND FAMILY, 5 TO 65 PERCENT SLOPES (VALLEYS)”, 6=”CHAIX-CHAWANAKEE FAMILIES-ROCK OUTCROP COMPLEX (HILLTOPS)”)

Grid definition:
ncols: 400
nrows: 200
xllcorner: 267999
yllcorner: 4148999
cellsize: 5 m

proj4:+init=epsg:26911 +proj=utm +zone=11 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0

Lineage:
The data set was obtained from the USGS National Map seamless server (http://seamless.usgs.gov). The map of soil mapping units was obtained from the Natural Resources Conservation Service (NRCS) Soil Data Mart (http://soildatamart.nrcs.usda.gov). The scripts used to predict soil mapping units and extract landforms are available via the authors website. The elevations range from 1400 to 1800 meters. There are six soil mapping units: (1) Holland family, 35 to 65% slopes; (2) Chaix-chawanakee family-rock outcrop complex; (3) Chaix family, deep, 5 to 25% slopes; (4) Chaix family, deep, 15 to 45% slopes, (5) Holland family, 5 to 65% slopes, valleys; (6) Chaix-chawanakee families-rock outcrop complex, hilltops.

Data owner: USGS
Reference: A Practical Guide to Geostatistical Mapping of Environmental Variables

Location:Fishcamp, United States
37° 27’ 48.708” N,119° 36’ 42.84” W
See map: Google Maps

Attachment:

fishcamp.zip

fishcamp_orthoimg.zip