Re: [R-br] Como extrair índices linha e coluna de um raster
Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler? A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster. Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....) Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima. Mais uma vez agradeço a lista. Att, Alex Santos Jônatam escreveu (alterei apenas a quantidade de linhas e colunas): require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
## estás concluindo isso da saída da linha de comando que há um limite de linhas para print. ## jogue em uma variavel o resultado do cbind e verifique res <- cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r)) dim(res) tail(res) ## de uma olhada também em ?rasterToPolygons A questão do shape file tem que testar, há diversas opções conforme mensagem anterior da lista. http://listas.inf.ufpr.br/pipermail/r-br/2013-April/010226.html 2013/7/3 Alex <alxcart@gmail.com>
Primeiramente obrigado a todos! Jônatan, por favor, a função cbind tem um limite para escrita dos dados de um raster, pois meu raster contém 319 linhas x 319 colunas, e a função interrompe no registro 24999, faltando 65001. Esta função é limitada ou gera todos os dados que ela se propõe a ler?
A princípio esta função (cbind) atende ao que desejo, pois gera as posições i e j, além de ler o valor do pixel correspondente no raster.
Abaixo está uma pequena parte da matriz gerada, após a leitura do meu raster. A coluna com os valores 0 (zero), junto com os valores de row e col são os que desejo no shapefile. matriz row col [1,] 1 1 1 0 [2,] 1 2 2 0 [3,] 1 3 3 0 (....)
Por favor, assim, em relação a conversão para shapefile, qual seria a melhor função para criar o shapefile (preferencialmente tipo polígono) com os atributos da matriz acima.
Mais uma vez agradeço a lista. Att, Alex Santos
Jônatam escreveu (alterei apenas a quantidade de linhas e colunas):
require(raster) ## criando um raster r <- raster(ncol=319,nrow=319) ## atribuindo valores r[] <- (1:ncell(r))*100 ## linha, coluna, celula, valor cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
_______________________________________________ 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.
-- ############################################################### ## 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 ###############################################################
participantes (2)
-
Alex -
Jônatan