Mauro, falta o carregar o pacote GDAL. Por favor, insira o comando no script require(rgdal). Att, Alex Santos Por que não rodaram os últimos comandos?

require(raster) Loading required package: raster Loading required package: sp

r <- raster(nc=20,nr=10) r[] <- rep(1:5,each=40) plot(r) res <- data.frame(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) names(res) <- c("row","col","cell","v") ## criando stack com 4 layers s <- stack(r,r,r,r) ## atribuindo valores para cada layer baseado no res s[] <- as.matrix(res) plot(s) #p<-rasterToPoints(s, spatial=T) pp<- rasterToPolygons(s,dissolve=F) #spplot(pp) # writeOGR(p,getwd(),"p2",driver="ESRI Shapefile") writeOGR(pp,getwd(),"matriz_pl",driver="ESRI Shapefile") Error: could not find function "writeOGR"

pescrito <- readOGR(getwd(), "matriz_pl") Error: could not find function "readOGR" pescrito@data Error: object 'pescrito' not found