[R-sig-Geo] getting mean values for each raster in a stack

[Antonio Silva]

Hi,

After many trials and some reading I got what I need

apply(vals[[1]],2,mean)

Anyway I would welcome comments to improve this solution.

Best regards

Antonio Olinto


2016-05-16 15:32 GMT-03:00 Antonio Silva <aolinto.lst at gmail.com>:

> Dear R users
>
> I have a raster stack and a polygon.
>
> I want the mean value inside this polygon for each raster in the stack.
>
> I'm doing the following:
>
> vals <- extract(stack,polygon)
> colMeans(r.vals[[1]]) # only for mean values
> # or
> for (i in 1:nlayers(stack)) {print(mean(vals[[1]][,i]))} # mean or other
> parameter
>
> I also tried to use lapply, but I got only the overall mean with
> lapply(vals, FUN=mean)
>
> I wonder if there is a way to get the mean (and other parameters) for each
> layer using lapply function.
>
> Thanks in advance.
>
> Antonio Olinto
>
>


-- 
Antônio Olinto Ávila da Silva
Biólogo / Oceanógrafo
Instituto de Pesca (Fisheries Institute)
São Paulo, Brasil

	[[alternative HTML version deleted]]