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