<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 href="https://dl.dropboxusercontent.com/u/27700634/precip.txt">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 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>