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

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