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.