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)