[R-br] Extração de dados climáticos

[Tito Conte]

Ana,
Que tal você usar o pacote netcdf? Ele é muito mais adequado para isso que
você deseja fazer (além de ser mais simples).
Segue o manual (é bem simples de fazer a extraçao de um ponto)
https://cran.r-project.org/web/packages/RNetCDF/RNetCDF.pdf






Tito Conte


Em 9 de outubro de 2016 20:09, Mauro Sznelwar via R-br <
r-br em listas.c3sl.ufpr.br> escreveu:

> Poderia enviar o arquivo correspondente para poder rodar?
>
>
> Prezados colegas, bom dia!
>
>
>
> Então, é o seguinte, no site http://berkeleyearth.org/ eu tenho dados de
> temperatura média de toda a terra. Eu preciso extrair de acordo com os meus
> pontos geográficos, a temperatura média anual para cada ponto. Só que eu
> tenho os específicos que eu preciso. Comecei o script de forma que eu
> pudesse extrair para um ponto específico e para todos os anos, para depois
> arrumar o script para fazer para todos os meus pontos e anos de interesse.
>
> Segue o script utilizado!
>
> library(raster)
>
> b <- brick("D:/Documentos/Desktop/analises_cap2/dados_climaticos/
> precipt.nc", varname <- "precip")
> print(b)
> b
>
> idx <- getZ(b)
> coords <- matrix(c(-55, -10), ncol = 2)
> vals <- extract(b, coords, df = T)
> df <- data.frame(idx, t(vals)[-1,])
>
> rownames(df) <- NULL
> names(df) <- c('date','value')
> head(df)
>
> Seguindo o script no local onde deveria apresentar as variáveis
> extraídas , aparece NA. Alguém já fez este tipo de extração? Peço ajuda.
>
> Desde já, muito obrigada.
>
> Ana Paula
>
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161010/0fa4e3ae/attachment.html>