Re: [R-br] Fwd: como extrair dados de arquivos GEOtiff

Olá Carlos, Estou respondendo também para a lista para com vistas a guardar a discussão que no futuro pode ajudar alguém. Fica difícil analisar o problema sem ter uma amostra real dos dados, mas pelo que você relata eu tentaria mudar a projeção da imagem. A mais comum, e que eu tentaria logo de cara, é a lat/lon WGS84. Você pode usar a função "projection" do próprio pacote raster, em um comando parecido com: projection(r) <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 " Saudações, Thiago. ________________________________ From: D'Apolito <carlosdapolito@gmail.com> To: Thiago Veloso <thi_veloso@yahoo.com.br> Sent: Tuesday, November 13, 2012 5:17 PM Subject: Re: [R-br] Fwd: como extrair dados de arquivos GEOtiff Olá Thiago, Realmente, dividindo-se tudo por 10 os valores condizem perfeitamente com o esperado. Eles devem ter subido os dados assim para pesar menos, ou sei lá por que. Agora meu problema é o seguinte, as coordenadas geográficas aparecem de 0-3600 pontos correspondendo a 30 graus = 120/grau. Mas quando uso plot(r) a escala de Lat e Lon aparecem normais, haverá uma maneira fácil de transformar esses dados, do contrário as tabelas em excel ficam demasiado grandes. Muito obrigado Carlos Em 13 de novembro de 2012 14:05, Thiago Veloso <thi_veloso@yahoo.com.br> escreveu: Carlos,
O procedimento parece estar correto. Se você está recebendo valores que não esperava, pode ser algum problema com o arquivo tif de entrada. Você tem certeza de que ele já expressa temperatura? Dependendo da fonte desse arquivo, você pode ter que fazer uma calibração dos números digitais para obter valores de temperatura fisicamente consistentes.
Como se trata de uma imagem, a primeira coisa que eu faria seria plotar o objeto raster (imagem) para fazer uma análise exploratória:
plot(r)
Thiago.
________________________________ From: D'Apolito <carlosdapolito@gmail.com> To: r-br@listas.c3sl.ufpr.br Sent: Tuesday, November 13, 2012 12:21 PM Subject: [R-br] Fwd: como extrair dados de arquivos GEOtiff
Prezados,
Tenho arquivos no formato GEOtiff (ex. anexo) e quero extrair dados de temperatura média para cada coordenada geográfica.
Estou usando isso até o momento:
library(raster) f <- "tmean1_23.tif" r <- raster(f) v = getValues(r, format='matrix')
Mas o que obtenho com getValues() é apenas uma matriz enorme com números que não tem nada a ver com temperatura. Meu nível de R não é muito avançado, por isso não estou conseguindo sair desse passo, agradeço qualquer ajuda.
Att Carlos
_______________________________________________ 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.
-- _____________________________________ Carlos D'Apolito - PhD candidate - ID:1239480 Earth Sciences School of Geography & Environmental Sciences Aston Webb Building University of Birmingham Birmingham B15 2TT UK phone: 44 - 07 919 673 173 CXD280@bham.ac.uk CV: http://lattes.cnpq.br/4557754746424026
participantes (1)
-
Thiago Veloso