search for: precipitation_mapname_in_grass

...ASK map of a specific area. # 1) a landuse map landuse<-readRAST6("landuse_mapname_in_grass") # 2) a catchment map which divide the area in several catchements catchment<-readRAST6("catchement_mapname_in_grass") # 3) a precipitation map precipitation<-readRAST6("precipitation_mapname_in_grass") # then I would like to sum the precipitation spatialy over each catchment and landuse. So, first I cbind all maps with cbind MAP<-cbind(landuse,catchment,precipitation) # then I use the aggregate function SUM<-aggregate(MAP[3],by=list(MAP[1],MAP[2]),sum,na.rm=TRUE) # here is the pr...