Obrigado Éder, Queria também extrair o valor dos pixels e adicionei a rotina, vou deixar o scrip para depósito na lista, ### <code r> require(raster); require(sp) ## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255)) ##Pontos selecionados x <- c(44.25,33.06,63.22,71.35) y <- c(22.00,81.90,15.18,71.03) pontos <- SpatialPoints(cbind(x,y)) cells <- cellFromXY(r, pontos); cells coordinates(r)[cells,] ## alternativa 1 p.cen <- xyFromCell(r, cells); p.cen ## alternativa2 plot(r) points(pontos, col="red") points(p.cen) values <- extract(r,p.cen) ##Extrai as coordendas dados<-cbind(p.cen,cells,values) dados ### </code> -- ====================================================================== 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== Em 16/01/2014 21:45, Éder Comunello escreveu:
Alexandre,
Pelo que entendi você quer a coordenada do centro do pixel/célula identificado. Se for isso mesmo posso sugerir duas opções, conforme código abaixo:
### <code r> require(raster); require(sp)
## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255))
##Pontos selecionados x <- c(44.25,33.06,63.22,71.35) y <- c(22.00,81.90,15.18,71.03) pontos <- SpatialPoints(cbind(x,y))
cells <- cellFromXY(r, pontos); cells coordinates(r)[cells,] ## alternativa 1 p.cen <- xyFromCell(r, cells); p.cen ## alternativa2
plot(r) points(pontos, col="red") points(p.cen) ### </code>
Éder Comunello <c <mailto:comunello.eder@gmail.com>omunello.eder@gmail.com <mailto:omunello.eder@gmail.com>> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ 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 forneça código mínimo reproduzível.