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

Alex alxcart em gmail.com
Quarta Julho 3 16:24:25 BRT 2013


Primeiramente obrigado a todos!
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 num shapefile (tipo polígono ou ponto).
matriz
         row col
    [1,]   1   1     1     0
    [2,]   1   2     2     0
    [3,]   1   3     3     0
(....)

      Por favor, em relação a conversão para shapefile, qual seria a melhor
função no R para criar o shapefile com os atributos da matriz acima.

      Abaixo está o código elaborado por Jônatam, com alterações apenas na
quantidade de linhas e colunas, agora resta apenas converter o resultado
"matriz" para shapefile.

## criando um raster
r <- raster(ncol=319,nrow=319)
## atribuindo valores
r[] <- (1:ncell(r))
## linha, coluna, celula, valor
cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))

# resultado matriz
matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
matriz


Mais uma vez agradeço a lista.
Att,
Alex Santos
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130703/3b11b95c/attachment.html>


Mais detalhes sobre a lista de discussão R-br