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