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