<html><head></head><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px"><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><span id="yui_3_16_0_1_1455747280644_30057">Na parte de leitura do ASCII, você não queria um stack a partir do r_ascii? Me parece que as duas funções extract estão fazendo a mesma coisa - extraindo os valores do stack s (criado do tiff) usando as coordenadas especificadas.</span></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><span><br></span></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><span id="yui_3_16_0_1_1455747280644_32072">Não rodei o script, mas tenho uns palpites baseados em experiências passadas. Além do tempo de processamento, veja o tamanho - em disco e memória - do objeto ASCII versus objeto raster. O tamanho em disco do ascii deve ficar na casa das unidades ou dezenas de vezes maior, e na memória (workspace do R) deve ser incomparável - da ordem de kilobytes para a imagem e centenas de megabytes para o ascii.</span></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><span><br></span></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058">O gerenciamento de memória do R é muito criticado em comparação a outras linguagens porque ele tende a carregar tudo na memória. Se você trabalha com grande quantidade de dados ou com dados de grande tamanho, você precisa tomar cuidado com isso sob o risco de encher a memória do computador com muita facilidade.</div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><br></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058">É por isso que as funções do pacote raster não carregam a imagem na memória por padrão. Digite r_tiff e veja o data source do objeto. Provavelmente deve ser um mero ponteiro para o arquivo em disco, ao invés de estar na memória. Usando opções pouco documentadas, você pode mudar esse comportamento do raster, permitindo a leitura direto na memória - o que acaba sendo mais rápido, mas de novo você precisa conhecer bem o tamanho do seus dados e da memória do seu computador para não acabar deixando-o uma carroça por falta de memória.</div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058"><br></div><div dir="ltr" id="yui_3_16_0_1_1455747280644_30058">Da maneira que eu vejo, essas funções para ler e escrever dados espaciais em ascii são para dados mais esparsos, como medidas de gps de um experimento de campo etc. Para trabalhar com dados espaciais mais densos, e em áreas geográficas mais extensas, como imagens de satélite ou dados climáticos, eu dificilmente enxergo uma justificativa para conversão para ascii.</div><div></div><div id="yui_3_16_0_1_1455747280644_29987"> </div><div class="signature" id="yui_3_16_0_1_1455747280644_29988">Greetings,<div id="yui_3_16_0_1_1455747280644_29989"> -- Thiago V. dos Santos</div><div id="yui_3_16_0_1_1455747280644_29990"><br></div><div id="yui_3_16_0_1_1455747280644_29991">PhD student</div><div id="yui_3_16_0_1_1455747280644_29992">Land and Atmospheric Science</div><div id="yui_3_16_0_1_1455747280644_29993">University of Minnesota</div></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> On Wednesday, February 17, 2016 4:58 PM, ASANTOS <alexandresantosbr@yahoo.com.br> wrote:<br></font></div> <br><br> <div class="y_msg_container"><div id="yiv0956943036"><div>
Boa noite Éder,<br clear="none">
<br clear="none">
Montei uma rotina abaixo para iniciarmos um novo tópico e é
possível observar o que você havia comentado, sendo:<br clear="none">
<br clear="none">
### <code r><br clear="none">
require(raster)<br clear="none">
require(sp)<br clear="none">
require(rgdal)<br clear="none">
<br clear="none">
# RasterLayer inventado<br clear="none">
r <- raster(nrows=10000, ncols=10000)<br clear="none">
r <- setValues(r, 1:ncell(r))<br clear="none">
r_SGDF<- as(r, 'SpatialGridDataFrame')<br clear="none">
writeGDAL(r_SGDF,"r_test.tif")### Cria o Tiff<br clear="none">
write.asciigrid(r_SGDF,"r_test.asc")### Cria o ASCII<br clear="none">
#<br clear="none">
<br clear="none">
# Lendo o tiff<br clear="none">
r_tif<-raster(c("r_test.tif"))<br clear="none">
s <- stack(r_tif)<br clear="none">
POI <- SpatialPoints(cbind(lon=c(-47.5,-47.5,-44.5,-40),
lat=c(-18.5,-19.5,-20.5,-22)))<br clear="none">
proj4string(POI) <- proj4string(s)<br clear="none">
e <- extract(s, POI, method= "bilinear", df=T) #ou
method="simple"<br clear="none">
proc.time()<br clear="none">
#<br clear="none">
<br clear="none">
# Lendo o ASCII<br clear="none">
r_ascii<-read.asciigrid("r_test.asc")<br clear="none">
POI <- SpatialPoints(cbind(lon=c(-47.5,-47.5,-44.5,-40),
lat=c(-18.5,-19.5,-20.5,-22)))<br clear="none">
proj4string(POI) <- proj4string(s)<br clear="none">
e <- extract(s, POI, method= "bilinear", df=T) #ou
method="simple"<br clear="none">
proc.time()<br clear="none">
### </code><br clear="none">
<br clear="none">
<pre class="yiv0956943036moz-signature">--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<a rel="nofollow" shape="rect" class="yiv0956943036moz-txt-link-abbreviated" ymailto="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a rel="nofollow" shape="rect" class="yiv0956943036moz-txt-link-abbreviated" ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a rel="nofollow" shape="rect" class="yiv0956943036moz-txt-link-freetext" target="_blank" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
OrcID: orcid.org/0000-0001-8232-6722
Researchgate: <a rel="nofollow" shape="rect" class="yiv0956943036moz-txt-link-freetext" target="_blank" href="https://www.researchgate.net/profile/Alexandre_Santos10">https://www.researchgate.net/profile/Alexandre_Santos10</a>
LinkedIn: <a rel="nofollow" shape="rect" class="yiv0956943036moz-txt-link-freetext" target="_blank" href="https://br.linkedin.com/in/alexandre-dos-santos-87961635">https://br.linkedin.com/in/alexandre-dos-santos-87961635</a>
======================================================================
</pre>
<div class="yiv0956943036moz-cite-prefix">Em 17/02/2016 16:07, Éder Comunello
escreveu:<br clear="none">
</div>
<div class="yiv0956943036yqt4825870984" id="yiv0956943036yqt40202"><blockquote type="cite">
<div dir="ltr">
<div class="yiv0956943036gmail_default" style="font-family:verdana, sans-serif;">Alexandre, boa tarde!</div>
<div class="yiv0956943036gmail_default" style="font-family:verdana, sans-serif;"><br clear="none">
</div>
<div class="yiv0956943036gmail_default" style="font-family:verdana, sans-serif;">Sobre o formato dos
dados, fiz uns testes aqui e dá pra tirar umas conclusões, mas
acho melhor criar um tópico novo pra não embolar.
Resumidamente, ao acessar usando o pacote raster (antes que
rgdal) você mantém a informação no disco e economiza memória.
Você poderia usar um formato ascii, mas os binários otimizam a
velocidade de leitura e acesso.</div>
<div class="yiv0956943036gmail_extra"><br clear="all">
<div>
<div class="yiv0956943036gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr"><font style="font-size:small;" face="arial, helvetica, sans-serif">
</font><div style="font-family:'trebuchet ms', sans-serif;display:inline;"></div>
<font style="font-size:12.8px;" face="arial, helvetica, sans-serif"><span style="font-size:small;">================================================</span><br clear="none">
Éder Comunello</font>
<div style="font-size:12.8px;"><span style="font-family:arial, helvetica, sans-serif;">Agronomist
(UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br clear="none">
</font></div>
<div style="font-size:12.8px;">
<div style="font-size:small;"><span style="font-family:arial, helvetica, sans-serif;">DSc
in Agricultural Systems
Engineering (USP/Esalq)</span><br clear="none">
</div>
<div><span style="font-size:small;">Brazilian
Agricultural Research Corporation (</span><font style="font-size:small;" face="arial, helvetica, sans-serif">Embrapa)</font></div>
<div>
<div><font style="font-size:small;" face="arial, helvetica, sans-serif">Dourados,
MS, Brazil </font><font size="2" face="arial, helvetica, sans-serif">|<O>|</font></div>
<div style="font-size:small;"><span style="font-family:arial, helvetica, sans-serif;">================================================</span><br clear="none">
</div>
</div>
</div>
<div style="font-size:small;"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752,
-54.8182, 408m<span style="font-family:arial, helvetica, sans-serif;"><br clear="none">
</span></div>
<div style="font-size:12.8px;"><span style="font-size:small;font-family:arial, helvetica, sans-serif;">UTC-04:00
/ DST: UTC-03:00</span></div>
<div>
<div>
<div><br clear="none">
</div>
<div><br clear="none">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote></div>
<br clear="none">
</div></div><br><div class="yqt4825870984" id="yqt11936">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne� c�igo m�imo reproduz�el.</div><br><br></div> </div> </div> </div></div></body></html>