[R-br] Extrair coordenadas de um raster dada uma condição

ASANTOS alexandresantosbr em yahoo.com.br
Domingo Fevereiro 1 23:48:11 BRST 2015


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:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================




---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com



Mais detalhes sobre a lista de discussão R-br