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