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.