Re: [R-br] Como extrair índices linha e coluna de um raster
Aos membros do R-br por favor, voltando ao meu propósito inicial eu preciso apenas de uma orientação em relação a uma função que retorne o índice (i,j) e o valor do pixel de todo raster no R, nada mais do que isto. O R possui bastante recursos para executar esta tarefa, mas infelizmente não encontrei solução. Depois disto exportar/converter para shapefile será mais simples. Grato, 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]
require(raster) ## criando um raster r <- raster(ncol=20,nrow=20) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) 2013/6/20 Elias Krainski <eliaskrainski@yahoo.com.br>
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]
______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012084 ## www.ufsm.br/meteorologia ###############################################################
good! On 06/20/2013 03:21 PM, Jônatan wrote:
require(raster) ## criando um raster r <- raster(ncol=20,nrow=20) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
2013/6/20 Elias Krainski <eliaskrainski@yahoo.com.br <mailto:eliaskrainski@yahoo.com.br>>
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]
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012084 ## www.ufsm.br/meteorologia <http://www.ufsm.br/meteorologia> ###############################################################
_______________________________________________ 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 -
Jônatan