<div dir="ltr">Ana,<br>Que tal você usar o pacote netcdf? Ele é muito mais adequado para isso que você deseja fazer (além de ser mais simples).<br>Segue o manual (é bem simples de fazer a extraçao de um ponto)<br><a href="https://cran.r-project.org/web/packages/RNetCDF/RNetCDF.pdf">https://cran.r-project.org/web/packages/RNetCDF/RNetCDF.pdf</a><br><br><br><br><div><br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature">Tito Conte<br><br></div></div>
<br><div class="gmail_quote">Em 9 de outubro de 2016 20:09, Mauro Sznelwar via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:medium;color:#0000ff">Poderia enviar o arquivo correspondente para poder rodar?</span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"> </span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"> </span></div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Prezados colegas, bom dia!</span></div>
<div id="m_-8654081205235771673divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif">
<div style="color:#000000">
<div>
<div id="m_-8654081205235771673divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<div style="color:#000000">
<div id="m_-8654081205235771673divRplyFwdMsg" dir="ltr">
<div> </div>
</div>
<div>
<div style="color:#000;background-color:#fff;font-family:times new roman,new york,times,serif;font-size:16px">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3550" class="m_-8654081205235771673yahoo_quoted" style="display:block">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3549" style="font-family:times new roman,new york,times,serif;font-size:16px">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3548" style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3547" dir="ltr"><span id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3551" style="font-family:Arial;font-size:small"><br> </span></div>
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3552" class="m_-8654081205235771673y_msg_container">
<div id="m_-8654081205235771673yiv0170119630">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3554">
<div id="m_-8654081205235771673yui_3_16_0_ym19_1_1474633639940_3553" style="color:#000;background-color:#fff;font-family:times new roman,new york,times,serif;font-size:16px">
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2915"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2915" dir="ltr">Então, é o seguinte, no site <a id="m_-8654081205235771673LPlnk55456" class="m_-8654081205235771673yiv0170119630" href="http://berkeleyearth.org/" rel="nofollow" target="_blank">http://berkeleyearth.org/</a><wbr> 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. </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_3476"> </div>
<div>Segue o script utilizado!</div>
<div> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2917">library(raster)</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2918"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2920">b <- brick("D:/Documentos/Desktop/<wbr>analises_cap2/dados_<wbr>climaticos/<a href="http://precipt.nc" target="_blank">precipt.nc</a>", varname <- "precip")</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2921">print(b)</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2922">b</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2923"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2925">idx <- getZ(b)</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2928">coords <- matrix(c(-55, -10), ncol = 2) </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2931">vals <- extract(b, coords, df = T)</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2934">df <- data.frame(idx, t(vals)[-1,])</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2935"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2937">rownames(df) <- NULL</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2938">names(df) <- c('date','value')</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2941">head(df)</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2942"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2944" dir="ltr">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. </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2944" dir="ltr"> </div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2944" dir="ltr">Desde já, muito obrigada.</div>
<div id="m_-8654081205235771673yiv0170119630yui_3_16_0_ym19_1_1474594594781_2944" dir="ltr"> </div>
<div dir="ltr">Ana Paula</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div><br>______________________________<wbr>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>