Re: [R-br] Como extrair índices linha e coluna de um raster
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
#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@data str(pescrito) 2013/7/3 Alex <alxcart@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@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