Carregar MDT no formato *.asc no R
Olá Pessoal, Preciso carregar um Modelo Digital do Terreno (raster) no formato *.asc no R gostaria de saber qual pacote é necessário e qual o procedimento. O arquivo quando aberto possui o cabeçalho abaixo: ncols 141 nrows 144 xllcorner 190611.999818 yllcorner 7536478.000059 cellsize 100.000000 nodata_value -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 Att, Hélder Gramacho *agrohelder@gmail.com <agrohelder@hotmail.com>*
eu uso o pacote "raster" entre as diversas opções disponíveis. 2013/11/14 Helder Gramacho <agrohelder@gmail.com>
Olá Pessoal,
Preciso carregar um Modelo Digital do Terreno (raster) no formato *.asc no R gostaria de saber qual pacote é necessário e qual o procedimento. O arquivo quando aberto possui o cabeçalho abaixo:
ncols 141 nrows 144 xllcorner 190611.999818 yllcorner 7536478.000059 cellsize 100.000000 nodata_value -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999
Att, Hélder Gramacho * agrohelder@gmail.com <agrohelder@hotmail.com>*
_______________________________________________ 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)33012083 ## www.ufsm.br/meteorologia ###############################################################
Boa noite, Outras opções de pacotes/comandos para leitura de grids (e rasters)... require(sp) ?read.asciigrid require(rgdal) ?readGDAL ### readGDAL() reconhece o formato automaticamente! Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 14 de novembro de 2013 17:17, Helder Gramacho <agrohelder@gmail.com>escreveu:
Olá Pessoal,
Preciso carregar um Modelo Digital do Terreno (raster) no formato *.asc no R gostaria de saber qual pacote é necessário e qual o procedimento. O arquivo quando aberto possui o cabeçalho abaixo:
ncols 141 nrows 144 xllcorner 190611.999818 yllcorner 7536478.000059 cellsize 100.000000 nodata_value -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999
Att, Hélder Gramacho * agrohelder@gmail.com <agrohelder@hotmail.com>*
_______________________________________________ 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.
Olá Éder, O readGDAL carregou e já transformou em *spatialgridDataframe* que era exatamente o que eu precisava, obrigado. Estou tentando realizar a krigagem interpolando a variavel* "ALF"* utilizando como grid o *"MDT"* da dúvida anterior, entretanto, está ocorrendo o seguinte erro:
OK = krige(Z ~ 1, ALF, MDT, model = modelo_ajustado, nmax=40) [using ordinary kriging]
"chfactor.c", *line 131: singular matrix in function LDLfactor()* Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, : LDLfactor Sabe dizer qual pode ser o problema e como pode ser contornado? Hélder Gramacho *agrohelder@gmail.com <agrohelder@hotmail.com>* Em 14 de novembro de 2013 20:34, Éder Comunello <comunello.eder@gmail.com>escreveu:
Boa noite,
Outras opções de pacotes/comandos para leitura de grids (e rasters)...
require(sp) ?read.asciigrid
require(rgdal) ?readGDAL ### readGDAL() reconhece o formato automaticamente!
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Em 14 de novembro de 2013 17:17, Helder Gramacho <agrohelder@gmail.com>escreveu:
Olá Pessoal,
Preciso carregar um Modelo Digital do Terreno (raster) no formato *.asc no R gostaria de saber qual pacote é necessário e qual o procedimento. O arquivo quando aberto possui o cabeçalho abaixo:
ncols 141 nrows 144 xllcorner 190611.999818 yllcorner 7536478.000059 cellsize 100.000000 nodata_value -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999
Att, Hélder Gramacho * agrohelder@gmail.com <agrohelder@hotmail.com>*
_______________________________________________ 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.
_______________________________________________ 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.
Hélder, Posso dar uma olhada, mas vou precisar de um código mínimo (CMR) e dos dados pra tentar reproduzir o erro aqui. Pra disponibilizar os dados, sugiro o dropbox. As cotas são mais generosas quando você usa um email de universidade. Éder -- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
participantes (3)
-
Helder Gramacho -
Jônatan -
Éder Comunello