Ola Eder,<div><br></div><div>O arquivo de origem dos dados (ETOPO1) `e muito grande, de maneira que tenho obtido erro constantemente, Tentei aumentar a capacidade da memoria com memory.limit() mas tambem nao funcionou.</div>
<div><br></div><div><br></div><div>att<br>Carlos </div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 20 de novembro de 2012 11:30, Eder Comunello <span dir="ltr"><<a href="mailto:ecomunel@gmail.com" target="_blank">ecomunel@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Carlos e colegas, bom dia!</div><div><br></div><div>Vou aproveitar um script anterior, fazendo algumas alterações. Como da outra vez, creio que dá pra adaptar facilmente pro seu caso.</div>
<div><br></div><div>Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no objeto origem ao final...</div>
<div><br></div><div>##### <begin> #####<br></div><div>require(sp); require(rgdal)</div><div><br></div><div>### Download do exemplo (399 Kb)</div><div>download.file('<a href="https://dl.dropbox.com/u/117618178/dourados/DDOS.tif" target="_blank">https://dl.dropbox.com/u/117618178/dourados/DDOS.tif</a>', </div>
<div> file.path(getwd(),'DDOS.tif'), mode='wb')</div><div><br></div><div>alt = readGDAL("DDOS.tif"); gridded(alt) ### 'importa' o geotiff</div><div>image(alt, col=terrain.colors(21)) ## visualização</div>
<div>slotNames(alt) ### nome dos slots (S4)</div><div>proj4string(alt) ### projeçãoo definida no geotiff</div><div>head(alt@data) ### observe os valores dos dados</div><div>head(alt@data$band1) ### observe os valores dos dados</div>
<div><br></div><div>### Opção para fatiar (0 ou 1)</div><div><br></div><div>alt2 <- alt ### para manipular uma cópia do objeto</div><div>z <- alt2@data$band1 ### para manipular apenas o vetor de interesse</div><div>
n <- length(z); n</div><div><br></div><div>for(i in 1:n){if (z[i]<400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1</div><div>alt2@data$band1 <- z ### substitui o vetor manipulado no objeto</div><div>image(alt2) ## visualização</div>
<div><br></div><div>### Opção para mascarar valores abaixo de critério (com 0 ou NA)</div><div><br></div><div>alt3 <- alt ### para manipular outra cópia do objeto</div><div>z <- alt3@data$band1 ### para manipular apenas o vetor de interesse</div>
<div>n <- length(z); n</div><div><br></div><div>for(i in 1:n){if (z[i]<400){z[i]=0}}; head(z) ### 0 ou 1</div><div>#for(i in 1:n){if (z[i]<400){z[i]=NA}}; head(z) ### 0 ou 1</div><div>alt3@data$band1 <- z ### substitui o vetor manipulado no objeto</div>
<div>image(alt3) ## visualização</div><div>##### <end> #####<br></div><div class="gmail_extra"><br clear="all"><font face="arial, helvetica, sans-serif"><br><span style="font-size:small">================================================</span><br>
Éder Comunello<br><br>Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)<br>Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]<br><br>Researcher at Embrapa Western Region Agriculture</font><div><font face="arial, helvetica, sans-serif">Dourados, MS, Brazil [</font>22 16.5'S, 54 49.0'W<span style="font-family:arial,helvetica,sans-serif">]</span></div>
<div><font face="arial, helvetica, sans-serif"><span style="font-size:small">================================================</span><br>UTC-03:00</font></div><br></div>
<br>_______________________________________________<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" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<font face="Tahoma">_____________________________________</font><br style="text-indent:0px!important"><b style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">Carlos D'Apolito - <span style="text-indent:0px!important;font-weight:normal">PhD candidate - ID:1239480</span></font></b><div>
<span style="font-family:Arial;font-size:x-small">Earth Sciences</span></div><div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">School of Geography & Environmental Sciences</font></div>
<div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">Aston Webb Building</font></div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important">
<font size="1" face="Arial" style="text-indent:0px!important">University of Birmingham</font></div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">Birmingham</font></div>
<div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">B15 2TT</font></div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important">
<font size="1" face="Arial" style="text-indent:0px!important">UK</font></div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important"><br></font></div>
<div style="font-family:Tahoma;font-size:13px;text-indent:0px!important"><font size="1" face="Arial" style="text-indent:0px!important">phone: 4</font><font size="1" face="Arial">4 - </font><font face="Arial" size="1">07 919 673 173</font></div>
<div style="font-family:Tahoma;text-indent:0px!important"><font size="1"><a href="mailto:CXD280@bham.ac.uk" target="_blank">CXD280@bham.ac.uk</a></font></div><div style="font-family:Tahoma;font-size:13px;text-indent:0px!important">
<span style="font-size:8pt">CV:</span> <span style="font-family:Arial,Helvetica,sans-serif;font-size:11px;color:rgb(102,102,102)"><a href="http://lattes.cnpq.br/4557754746424026" target="_blank">http://lattes.cnpq.br/4557754746424026</a></span> </div>
<div style="text-indent:0px!important"><br></div></div><br>
</div>