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

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@gmail.com> To: r-br@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@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.
participantes (2)
-
D'Apolito
-
Thiago Veloso