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)