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@yahoo.com.br alexandre.santos@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/ ======================================================================