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