<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">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>