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

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.

Alex, A princípio a função rasterToPolygons faria o que você precisa. Existe algum motivo especial para manter os índices das células no shapefile? Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/inde... Phone: (612) 323 9898 ________________________________ From: Alex <alxcart@gmail.com> To: r-br@listas.c3sl.ufpr.br Sent: Wednesday, June 19, 2013 7:49 AM Subject: [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 _______________________________________________ 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.
participantes (3)
-
Alex
-
Elias Krainski
-
Thiago V. dos Santos