[R-br] Extração de valores de pixels de interesse em uma área circular usando extract()

ASANTOS alexandre.santos em cas.ifmt.edu.br
Quinta Fevereiro 9 11:57:18 BRST 2017


Prezados membros,

             Gostaria de extrair a média de valores dos pixels de 
interesse em um raster, dadas as coordenadas centrais de áreas 
circulares e a determinação do raio, que é dada pelo buffer da função 
extract() do pacote raster, porém da forma como estou fazendo não esta 
funcionando, em meu exemplo:

require(raster)
require(sp)
require(plotrix)# Somente para representar a área alvo

## Crio um raster
r <- raster(nc=30, nr=30)
r <- setValues(r, round(runif(ncell(r))* 255))
plot(r)


#Seleciono o centro das áreas circulares
xd     <- c(-24.99270,45.12069,99.40321,73.64419)
yd  <- c(-45.435267,-88.369745,-7.086949,44.174530)
pts <- data.frame(xd,yd)
pts_s<- SpatialPoints(pts)

#Representação da área alvo
for (i in 1:length(xd)) {
draw.circle(pts$xd[i],pts$yd[i],20,border="blue",col="grey")
}


#Visualização das coordenadas do centro

points(pts_s, col="red", pch=16) #Visualização das coordenadas do centro


#Extração da média dos valores de pixels em uma parcela de raio = 20
r_sample <- extract(x = r, y = pts_s, buffer=20, fun=mean, df=TRUE)
r_sample

Não funciona quando realizado dessa maneira. Alguém poderia me dar um help?

Obrigado

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
OrcID: orcid.org/0000-0001-8232-6722   -   ResearcherID: A-5790-2016
Researchgate: www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
======================================================================



Mais detalhes sobre a lista de discussão R-br