[R-br] excluir mar no ETOPO

[D'Apolito]

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 em 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 em data) ### observe os valores dos dados
> head(alt em 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 em 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 em 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 em 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 em 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 em 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 em bham.ac.uk
CV: http://lattes.cnpq.br/4557754746424026
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/615af15d/attachment.html>