Caros Listeiros,
Fiz a seleção da coordenada dos pixels (p.cen) em um raster
que continham pontos de interesse (pontos) no CRM abaixo, gostaria de
saber como impor uma condição para selecionar as coordenadas dos pixels
que não contém nenhum ponto. Tentei criar um novo objeto com todas as
coordenadas no raster e remover as coordenadas p.cen sem sucesso, alguém
poderia me ajudar? Obrigado,
### <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,]
p.cen <- xyFromCell(r, cells); p.cen
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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com
_______________________________________________
R-br mailing list
[email protected]
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.
