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

[Elias Krainski]

Oi Alex,

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

   require(raster)
   example(rasterToPolygons)

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

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

   tail(pol em data,3)
   r[18,   tail(pol em 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 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/20130619/5a8c4258/attachment.html>