Thiago,
No meu caso eu extrai qual é o pixel pela função cellFromXY() ,
como o exemplo do Luciano, mas alem do número do pixel gostaria que
fosse retornada as coordenadas do pixel, teria alguma sugestão,
segue CRM,
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))#cbind-junta duas colunas em
uma, no caso as coordenadas
cellFromXY(r,pontos)
plot(r)
points(pontos, col="red")
Obrigado,
Alexandre
Em 16/01/2014 16:00, Luciano F. Sgarbi
escreveu:
Muito
obrigado Thiago.
Era
exatamente isso.
Abraço
_______________________________________________
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.
--
======================================================================
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
======================================================================