Ola Eder, 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. att Carlos Em 20 de novembro de 2012 11:30, Eder Comunello <ecomunel@gmail.com>escreveu:

Carlos e colegas, bom dia!

Vou aproveitar um script anterior, fazendo algumas alterações. Como da outra vez, creio que dá pra adaptar facilmente pro seu caso.

Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no objeto origem ao final...

##### <begin> ##### require(sp); require(rgdal)

### Download do exemplo (399 Kb) download.file('https://dl.dropbox.com/u/117618178/dourados/DDOS.tif', file.path(getwd(),'DDOS.tif'), mode='wb')

alt = readGDAL("DDOS.tif"); gridded(alt) ### 'importa' o geotiff image(alt, col=terrain.colors(21)) ## visualização slotNames(alt) ### nome dos slots (S4) proj4string(alt) ### projeçãoo definida no geotiff head(alt@data) ### observe os valores dos dados head(alt@data$band1) ### observe os valores dos dados

### Opção para fatiar (0 ou 1)

alt2 <- alt ### para manipular uma cópia do objeto z <- alt2@data$band1 ### para manipular apenas o vetor de interesse n <- length(z); n

for(i in 1:n){if (z[i]<400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1 alt2@data$band1 <- z ### substitui o vetor manipulado no objeto image(alt2) ## visualização

### Opção para mascarar valores abaixo de critério (com 0 ou NA)

alt3 <- alt ### para manipular outra cópia do objeto z <- alt3@data$band1 ### para manipular apenas o vetor de interesse n <- length(z); n

for(i in 1:n){if (z[i]<400){z[i]=0}}; head(z) ### 0 ou 1 #for(i in 1:n){if (z[i]<400){z[i]=NA}}; head(z) ### 0 ou 1 alt3@data$band1 <- z ### substitui o vetor manipulado no objeto image(alt3) ## visualização ##### <end> #####

================================================ Éder Comunello

Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

_______________________________________________ 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.

-- _____________________________________ *Carlos D'Apolito - PhD candidate - ID:1239480* Earth Sciences School of Geography & Environmental Sciences Aston Webb Building University of Birmingham Birmingham B15 2TT UK phone: 44 - 07 919 673 173 CXD280@bham.ac.uk CV: http://lattes.cnpq.br/4557754746424026