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

[Elias Teixeira Krainski]

Ver também a função over(), do pacote sp.

On 06/07/13 17:46, Alex wrote:
> 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)
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130708/78c71e3d/attachment.html>