[GRASSLIST:5810] create new map of TS with min value month

Martin Wegmann wegmann at biozentrum.uni-wuerzburg.de
Wed Feb 16 18:12:45 EST 2005


Hello, 

I have a time series and I want to create a map containing the information 
which decade has the annual minimum value. 

e.g.:

values for one pixel

jan 1 = 10
jan 2 = 11
jan 3 = 18
feb 1 = 9
feb 2 = 3
feb 3 = 5
.....

output =   22 or 5 (either 2= feb ; 2=2. decade OR 5th decade)

The ouput map would hold the information for each pixel which decade has the 
annual minimum value.


I tried a few approaches with 
r.mapcalc (if(map_x==map_annual_min,x,null()) )
r.mapcalc (if(map_y==map_annual_min,y,null()) )

and r.patch for all raster maps

however before I start writing a script with a if condition for each decade 
map I would like to know if somebody knows a more straight forward method. 

cheers Martin

P.S.: would be a nice extra functionality for r.series ,-)








-- 
Dipl.-Biol. Martin Wegmann

DLR - German Aerospace Center
German Remote Sensing Data Center
@
Dept.of Geography
Remote Sensing and Biodiversity Unit
University of Wuerzburg
Am Hubland
97074 Würzburg

phone: +49-(0)931 - 888 4797
fax:   +49-(0)931 - 888 4961
http://www.biota-africa.org
http://www.biogis.de




More information about the grass-user mailing list