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@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@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.