[R-br] excluir mar no ETOPO

[Thiago Veloso]

Carlos,

Eu tentaria algo assim:

library (raster) #carrega pacote exigido
r <- raster('ETOPO1_Ice_g_geotiff.tif') #carrega arquivo de trabalho
r [ r<0 ] <- NA #atribui NA a todos os valores negativos no raster.

E caso você queira atribuir NA a um escopo de valores, pode ser mais fácil criar uma função:


vals <- c(0,24,32,56) # estes são os valores que você quer eliminar da tabela de dados do raster
fun <- function (x) {
  x <- ifelse (x %in% vals, 1, NA) # aqui você substitui os valores acima por "1" e atribui NA para todo o resto. Você pode substituir o "1" por vals
}

E então aplicá-la sobre o raster:
r.new <- calc (r, fun)


Uma vantagem de usar funções é pode realizar uma tarefa para diversos arquivos sem precisar repetir comandos para cada um deles.

Espero que ajude,
Thiago.


________________________________
 From: D'Apolito <carlosdapolito em gmail.com>
To: r-br em listas.c3sl.ufpr.br 
Sent: Monday, November 19, 2012 8:15 PM
Subject: [R-br] excluir mar no ETOPO
 

Prezados,

Estou graficando mapas com os dados do ETOPO01 de elevação e gostaria de excluir o mar, para que os dados iniciem em zero ao invés de números negativos (profundidades dos oceanos).
Qualquer ajuda é bem vinda, desde já agradeço

Att
Carlos

#código:

library(raster) 
library(rgdal)

etopo<-"ETOPO1_Ice_g_geotiff.tif"
r.etopo<-raster(etopo)

ext <- extent(-95,-35,-31,31)

zoom(r.etopo, ext=ext,col=rev(gray.colors(12)))





_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/675bda4e/attachment.html>