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

[Jônatan]

#Tu olhou  ?rasterToshapefile ?
# com essa funcao tu não precisa gerar a matriz, faz direto a conversão
raster p/ pol.
#Exemplo:

require(raster)
require(rgdal)
r[] <- rep(1:5,each=40)
plot(r)
as.matrix(r)
## covertendo raster para poligonos
 p<-rasterToPolygons(r,dissolve=T)
class(p)
plot(r)
plot(p, add=T,lwd=2)
 writeOGR(p,getwd(),"p",driver="ESRI Shapefile")
pescrito <- readOGR(getwd(), "p")
pescrito em data
str(pescrito)





2013/7/3 Alex <alxcart em gmail.com>

> Jônatan,
> respondendo as suas perguntas:
> ## estás concluindo isso da saída da linha de comando que há um limite de
> linhas para print.
> Sim, realmente estou concluindo em função do 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)
> Abaixo está os resultados do código que você solicitou.
> [24999,]  84  99 24999  9895
>  [ atingido getOption("max.print") -- últimas 65001 linhas omitidas ]]
>
> > matriz<-cbind(rowColFromCell(r,1:ncell(r)), 1:ncell(r),getValues(r))
>
> > dim(matriz)
> [1] 90000     4
>
> > tail(matriz)
>          row col
> [89995,] 300 295 89995 0
> [89996,] 300 296 89996 0
> [89997,] 300 297 89997 0
> [89998,] 300 298 89998 0
> [89999,] 300 299 89999 0
> [90000,] 300 300 90000 0
>
> Em relação a dica do shapefile, tentarei utilizar a função writeOGR()
> sugerida pelo autor, apenas não sei como criar os atributos. Conforme você
> sugeriu o link:
> http://listas.inf.ufpr.br/pipermail/r-br/2013-April/010226.html
>
> A questão agora é transferir este resultado "matriz" para um shapefile.
> Caso eu tenha êxito ao gerar o shapefile, a coluna referente aos valores
> do pixel do meus raster contem valores que eu posso utilizar como join.
>
> Mais uma vez obrigado!
> Alex Santos
>
> _______________________________________________
> R-br mailing list
> R-br em 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
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130706/cc49d5fe/attachment-0001.html>