I want to get some statistics that r.stats can't directly come up
with, and I'm wondering if there is another function to do the trick.
I have two raster maps: one has polygons defining different "landscape
units". The second is a crop yield "surface" from the same area.
I want to know, for each landscape unit, not simply the total area,
but also the mean yield, and standard deviation, etc. So I need a
routine that uses the landscape units map as a "mask", and calculates
descriptive statistics from the raster values of the yield map found
within the area of each polygon.
Any suggestions?
Tim.
-----------------------------------------------------------------
Tim Martin *
Spatial Information Systems * If you want peace,
University of Alberta * prepare for justice.
martin@ulysses.sis.ualberta.ca *
-----------------------------------------------------------------