Extração de valores de pixels de interesse em uma área circular usando extract()
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/ ======================================================================
Leia o help da extract, em partcular sobre o arg buffer; veja a resolução de seu raster do exemplo. A extract está funcionando como esperado para o valor de buffer fornecido. Att. On Thursday, February 9, 2017, ASANTOS via R-br <r-br@listas.c3sl.ufpr.br> wrote:
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/ ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo mnimo reproduzvel.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria - UFSM ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012083 ## www.ufsm.br/meteorologia ###############################################################
participantes (2)
-
ASANTOS -
Jônatan