Raster aggregate values

From GRASS-Wiki
Revision as of 08:48, 26 August 2010 by Neteler (talk | contribs) (new, based on Hamish's notes)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Q: How can I calculate the average map values by spatial unit (grid)? For example, I would like to calculate the average slope by square kilometer for a given area.

A: Procedure outline:

  • Make a slope map from the DEM with r.slope.aspect at the same resolution and bounds as the DEM. (run 'g.region rast=dem' first)
  • Make the grid with v.mkgrid, with box=1000,1000 (assuming the location is meters-based)
  • Draw a 1 km x 1 km grid (d.grid), and overlay the v.mkgrid vector map (area fill = none, or only display boundary features, colored red or so) just to verify where it will be and what it will look like.
  • Run v.rast.stats to calculate stats for each grid box, and upload to that area's attribute table.
Retrieved from "https://grasswiki.osgeo.org/w/index.php?title=Raster_aggregate_values&oldid=11338"