[R-br] Como extrair índices linha e coluna de um raster

[Alex]

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 em data
> Error: object 'pescrito' not found
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130707/6eac1f5d/attachment.html>