[R-br] Remoção de coordinates() em SpatialGridDataFrame de diferentes dimensões

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Setembro 30 19:05:59 BRT 2015


Caros Listeiros,

             Eu criei dos raster r1 e r2, respectivamente com dimensões 
10x10 e 15x15, porém eu converti os rasters em SpatialGridDataFrame e 
gostaria de remover as coordenadas e respectivos valores do objeto 15x15 
para que ele fique com o mesmo tamanho e nas mesmas coordenadas que o 
objeto 10x10, porém não estou conseguindo determinar quais pixels estão 
contidos em no objeto r1.grid que não existe em r2.grid, segue CRM 
abaixo, alguém poderia me ajudar?


### <code r>
require(raster); require(rgdal); require(dplyr)

## Criando dois raster GTiff
#10x10
r1 <- raster(nc=10, nr=10)
r1 <- setValues(r1, round(runif(ncell(r))* 255))
writeRaster(r1, filename="r1.tif", format="GTiff", overwrite=TRUE)

#15x15
r2 <- raster(nc=15, nr=15)
r2 <- setValues(r2, round(runif(ncell(r))* 255))
writeRaster(r2, filename="r2.tif", format="GTiff", overwrite=TRUE)

#Converter para SpatialGridDataFrame
r1.gtiff<-readGDAL('r1.tif')
r1.grid<- as(r1.gtiff, 'SpatialGridDataFrame')
r2.gtiff<-readGDAL('r2.tif')
r2.grid<- as(r2.gtiff, 'SpatialGridDataFrame')

#Detecção de quais pixels estão contidos em r1.grid e r2.grid
df.pixels<-anti_join(coordinates(r1.grid),coordinates(r2.grid), 
copy=TRUE, replace = TRUE)

#Remoção das 125 coordenadas e do valor dos pixels de r2.grid para 
deixar os dois SpatialGridDataFrame com mesmo tamanho
new.r2.grid<-r2.grid[which(df.pixels), ]


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 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
======================================================================



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