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

[Thiago Veloso]

  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 em gmail.com>
To: Thiago Veloso <thi_veloso em 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 em 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 em gmail.com>
>To: r-br em 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 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.
>
>


-- 
_____________________________________
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 em bham.ac.uk
CV: http://lattes.cnpq.br/4557754746424026