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

[Alex]

Jônatan,
a função rasterToshapefile é interessante, mas não atribuí o índices i e j
ao shapefile criado.
       Encontrei outra solução no script abaixo e espero auxilar outras
pessoas.
       Apenas uma ressalva, ao final gero um shapefile de pontos, com os
atributos i, j e o valor da célula, mas depois associo ao shapefile de
polígonos no Quantum GIS através da função Join de Tabelas.
        Outra solução possível é juntar a sua solução de rasterToPolygon e
depois realizar o Join de tabela com os shapefiles de pontos.
        Realmente seria interessante que a função rasterToPolygon levasse
as informações de linha e coluna para o shapefile criado.

Mais uma vez obrigado Jôntan e Elias, vocês me orientaram bastante na
solução deste problema.

require(raster)
require(rgdal)
require(shapefiles)
require(maptools)

### ENTRADA DE DADOS ###
## criando um raster
r <- raster(ncol=319,nrow=319)
atribuindo valores
r[] <- (1:ncell(r))

## linha, coluna, celula, valor
cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
matriz
dim(matriz)
tail(matriz)
summary (matriz)
names(matriz)  # a matriz não possui nome de campos

# Data Frame possui nome dos atributos
class(matriz)
require(sp)
     pontos<-data.frame(matriz)
     head(pontos)
     class(pontos)
     dim(pontos)
     names(pontos)
     coordinates(pontos) <- c("row","col") # c("Norte_Sul","Leste_Oeste")
     class(pontos)
     dim(pontos)
     getSlots("SpatialPointsDataFrame")
     slot(pontos, "bbox")

require(maptools)
writePointsShape(pontos,"matriz319pt", factor2char = TRUE, max_nchar=254)
readShapePoints("matriz319pt", proj4string = CRS(as.character(NA)),verbose =
FALSE,repair=FALSE)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130706/f99e466d/attachment.html>