<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div class="" style=""><span class="" style="">Caro Wenceslau,</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;" class=""><span class="" style=""><br class="" style=""></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;" class="">Também estou buscando mudar a escala de cores para algo mais didático, mas ainda não acertei totalmente as configurações. Por enquanto, o que estou fazendo é:</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica
Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;" class=""><br class="" style=""></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;" class=""># criar uma escala de azul usando o RColorBrewer.</div><div style="background-color: transparent;" class=""># digite display.brewer.all() para ver as paletas disponíveis</div><div style="background-color: transparent;" class="">cols <- colorRampPalette(brewer.pal(9,"Blues"))(200)</div><div class="" style=""><br class="" style=""></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;" class=""># plotar o perfil especificando a paleta no argumento "cols"</div><div
style="background-color: transparent;" class="">plot <- filled.contour(x= date, </div><div style="background-color: transparent;" class=""> y= sort(-(as.numeric(names(depths)))), </div><div style="background-color: transparent;" class=""> z= t(soil_data[11:1,]),</div><div style="background-color: transparent;" class=""> zlim= c(0,1),</div><div style="background-color: transparent;" class=""> ylab= expression(theta(m^3~~m^-3)),</div><div style="background-color: transparent;" class=""> col=cols,</div><div style="background-color:
transparent;" class=""> nlevels= 100)</div><div class="" style=""><br class="" style=""></div><div class="" style=""></div><div class="" style="">Ainda não é o ideal, mas pelo menos uma escala azul torna mais fácil a interpretação de diferentes quantidades de umidade no solo.</div><div class="" style=""><br></div><div class="" style="">Saudações,<br class="" style="">--<br class="" style="">Thiago V. dos Santos<br class="" style="">PhD student<br class="" style="">Land and Atmospheric Science<br class="" style="">University of Minnesota<br class="" style="">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm<br class="" style="">Phone: (612) 323 9898</div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande,
sans-serif; font-size: 12pt;" class=""> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;" class=""> <div dir="ltr" class="" style=""> <font size="2" face="Arial" class="" style=""> On Thursday, September 18, 2014 9:46 AM, Wenceslau Teixeira <wgt007@gmail.com> wrote:<br class="" style=""> </font> </div> <br class="" style=""><br class="" style=""> <div class="" style=""><div id="yiv3021326806" class="" style=""><div dir="ltr" class="" style=""><div class="" style="">Caros colegas do R,<br class="" style=""><br class="" style="">muito bacana o script disponibilizado para traçar os perfis de umidade do solo. <br class="" style="">Será útil para nós. Só gostaria de trocar as cores, com valores mais alto (umido) para azul e valores mais baixos (solo seco)<br class="" style="">para vermelho, assim o gráfico fica mais cognitivo. <br class="" style="">Alguém sabe como mudar as
cores ?<br class="" style=""><br class="" style=""></div>Cordialmente,<br class="" style=""><br class="" style="">Wenceslau Teixeira<br class="" style="">Embrapa Solos<br class="" style=""><br class="" style=""><div class="" style=""><div class="" style=""><br class="" style=""><br class="" style="">library(repmis) # reads text data directly from dropbox - no need to download any file<br class="" style=""># read data<br class="" style="">url <- '<a rel="nofollow" target="_blank" href="https://dl.dropboxusercontent.com/u/27700634/precip.txt" class="" style="">https://dl.dropboxusercontent.com/u/27700634/precip.txt</a>'<br class="" style="">tmp <- repmis::source_data(url, sep = '', header = TRUE)<br class="" style=""><br class="" style=""># convert julian day to date<br class="" style="">date <- as.Date(tmp$julian, origin='2011-12-31')<br class="" style="">data <- cbind(date, tmp)<br class="" style="">head(data)<br class="" style=""># define
vector with depth of soil layers<br class="" style="">depths <- c(0.05,0.10,0.20,0.30,<br class="" style=""> 0.40,0.60,0.80,1.00,<br class="" style=""> 1.50,2.00,2.50)<br class="" style=""><br class="" style=""># teste com lattice<br class="" style="">cols <- ncol(data):4<br class="" style="">levelplot(as.matrix(data[, cols]),<br class="" style=""> aspect = "fill", xlab = "date", ylab = "soil Moist.",<br class="" style=""> scales = list(x = list(at = which(date %in% pretty(date)),<br class=""
style=""> labels = pretty(date)),<br class="" style=""> y = list(at = length(cols):1, labels = depths)),<br class="" style=""> contour = T)<br class="" style=""><br class="" style=""># Plot soil moisture profile<br class="" style=""><br class="" style=""># como são profundidades (abaixo da superfície) use a convenção de z negativo<br class="" style=""># acrescentei a chuva, precisa ajustar ylabels, mtext() resolve<br class="" style=""><br class="" style="">x11()<br class="" style="">op <- par(mar = c(3,4,12,3))<br class="" style="">plot <- filled.contour(date,
sort(-depths), as.matrix(data[, cols]), ylab = expression(theta(m^3~~m^-3)))<br class="" style="">par(new=T,mar=c(23,4,1,8))<br class="" style="">with(data,<br class="" style=""> plot(date, precip, type = "h", xaxt = "n", xpd = T, xaxs="i", ylab = "Precip (mm)", xlab = ""))<br class="" style="">par(op)<br class="" style=""><br clear="all" class="" style=""><br class="" style="">-- <br class="" style=""><div dir="ltr" class="" style=""><br class="" style="">Wenceslau Teixeira<br class="" style=""><br class="" style=""><a rel="nofollow" ymailto="mailto:wgt007@gmail.com" target="_blank" href="mailto:wgt007@gmail.com" class="" style="">wgt007@gmail.com</a><br class="" style="">skype: wgt007<br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class=""
style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""><br class="" style=""></div>
</div></div></div></div><br class="" style="">_______________________________________________<br class="" style="">R-br mailing list<br class="" style=""><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br" class="" style="">R-br@listas.c3sl.ufpr.br</a><br class="" style=""><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank" class="" style="">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br class="" style="">Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank" class="" style="">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br class="" style=""><br class="" style=""></div> </div> </div> </div> </div></body></html>