Oi Alex,

Veja este exemplo aplicado ao raster do exemplo da funcao citada.

  require(raster)
  example(rasterToPolygons)

  pol@data$i <- col(matrix(1L,36,18))[which(r@data@values>6)]
  pol@data$j <- row(matrix(1L,36,18))[which(r@data@values>6)]

  head(pol@data, 3)
  r[1, head(pol@data, 3)$j]

  tail(pol@data,3)
  r[18,   tail(pol@data,3)$j]


On 06/19/2013 09:49 AM, Alex wrote:
Aos membros da lista.
Por favor, eu gostaria de criar a partir de um dado raster um shapefile (polígono)  com os seguintes atributos:
- A posição (i, j) de cada célula e os respectivos valores de pixel para compor os atributos X, Y, Z com estas informações.
           No R existe a conversão de raster para polígono (rasterToPolygons), mas não encontrei uma função para extrair os índices (i,j) das células do raster e colocar como atributo no shapefile.

Nota: A princípio parece ser uma rotina simples, mas inexistente até mesmos nos melhores softwares de geoprocessamento, encontrei apenas no FME (comercial) mas gostaria de realizá-la em software livre. 

Desde já agradeço, 
Alex Santos



_______________________________________________
R-br mailing list
R-br@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.