Custom search

Samstag, 10. Mai 2014

Sun, terrain, open data and open source - some ingredients for fieldwork

This year growing season has already begun a few weeks ago. So it's time for field work.

A preliminary visual delineation of habitats based upon digital orthophotos is a common task for an ecologist before fieldwork.

But how valid are these visual delineations?

Ecological site conditions of habitats are often influenced amongst others by terrain forms, solar irradiance and irradiation.

Let's try some nice ingredients to facilitate and prepare fieldwork.

Open data

Open Source

and its nice tools and addons
  • r.sun.hourly - Runs r.sun in loop for given time range.
  • r.sun - Solar irradiance and irradiation model.
  • r.geomorphon - Calculate geomorphons (terrain forms)and associated geometry using machine vision approach.
  • Temporal data processing in GRASS GIS.
  • g.gui.animation - Tool for animating a series of raster and vector maps or a space time raster or vector dataset.
and for scripting

Step 1: Visually delineation of some habitats

Orthophotos of Tyrol as WMS-service.

Step 2: Load a DEM for further analysis

Digital terrain model of Tyrol (10m x 10m).

Step 3: v.habitat.dem - a tool for DEM derived characteristics of habitats

v.habitat.dem-GUI

v.habitat.dem-GUI

There are a terrain derived raster results, some of them are uploaded to the attribute table of the habitat vector:
Calculation of accumulation, drainage direction, topographic index by r.watershed ...
[...]
Calculation of slope and aspect by r.slope.aspect ...
[...]
Calculate habitat's areas and populate it to the attribute table ...
[...]
Habitat areas smaller than pixel area...
[...]
Upload DEM zonal statistics to the attribute table ...
[...]
Upload slope zonal statistics to the attribute table ...
[...]
Upload aspect zonal statistics to the attribute table ...
[...]
Do some simple checks regarding aspect range and populate it to the attribute table...
aspect range 100-200 *
aspect range 201-300 **
aspect range >= 300 ***
[...]
Do some simple checks regarding aspect range and slope median and populate it to the attribute table...
aspect range 100-200 and median slope < 5 *
aspect range 201-300 and median slope < 5 **
aspect range >= 300 and median slope < 5 ***
[...]
Calculate the most common geomorphons ...
[...]
Calculate mutual occurrences of habitats and geomorphons ...
[...]
Join geomorphon information to habitat attribute table ...
[...]
Run r.sun.hourly in binary mode for light/shadow for a certain day in the year ...
[...]
Registering created maps into temporal dataset...

 The most common terrain forms calculated by r.geomorphon are:
  • flat
  • summit
  • ridge
  • shoulder
  • spur
  •  slope
  •  hollow
  • footslope
  • valley
  • depression
for more details see the r.geomorphon-manual


r.geomorphon-results

The accumulation indicates the number of cells that drain through each cell and subsequently some information about moisture.

accumulation map

Step 4: Animation of direct sunlight / shadows by terrain for a certain day in the year


Temporal data processing was newly added in GRASS GIS. The results of the r.sun.hourly-analysis are autmatically registered into a temporal database. The space time raster dataset can be easily animated in the g.gui.animation-tool.


animated r.sun.hourly-results: 1st January 2014,
1 hour time step, start time 8:00 morning, end time 18:00 evening

Keep in mind the results are as good as the DEM is. But there are some nice additional informations for planning field work.

2 Kommentare:

  1. svn commit :-)

    http://grasswiki.osgeo.org/wiki/AddOns/GRASS7/vector#v.habitat.dem

    AntwortenLöschen
  2. a simple form of the LS factor (ratio of soil loss ) included now in calculations

    AntwortenLöschen