<div dir="ltr"><div>Segue algumas alternativas:</div><div><br></div><div>library(repmis) # reads text data directly from dropbox - no need to download any file</div><div># read data</div><div>url <- '<a href="https://dl.dropboxusercontent.com/u/27700634/precip.txt">https://dl.dropboxusercontent.com/u/27700634/precip.txt</a>'</div><div>tmp <- repmis::source_data(url, sep = '', header = TRUE)</div><div><br></div><div># convert julian day to date</div><div>date <- as.Date(tmp$julian, origin='2011-12-31')</div><div>data <- cbind(date, tmp)</div><div>head(data)</div><div># define vector with depth of soil layers</div><div>depths <- c(0.05,0.10,0.20,0.30,</div><div> 0.40,0.60,0.80,1.00,</div><div> 1.50,2.00,2.50)</div><div><br></div><div># teste com lattice</div><div>cols <- ncol(data):4</div><div>levelplot(as.matrix(data[, cols]), </div><div> aspect = "fill", xlab = "date", ylab = "soil Moist.",</div><div> scales = list(x = list(at = which(date %in% pretty(date)), </div><div> labels = pretty(date)),</div><div> y = list(at = length(cols):1, labels = depths)),</div><div> contour = T)</div><div><br></div><div># Plot soil moisture profile</div><div><br></div><div># como são profundidades (abaixo da superfície) use a convenção de z negativo </div><div># acrescentei a chuva, precisa ajustar ylabels, mtext() resolve</div><div><br></div><div>x11()</div><div>op <- par(mar = c(3,4,12,3))</div><div>plot <- filled.contour(date, sort(-depths), as.matrix(data[, cols]), ylab = expression(theta(m^3~~m^-3)))</div><div>par(new=T,mar=c(23,4,1,8))</div><div> with(data,</div><div> plot(date, precip, type = "h", xaxt = "n", xpd = T, xaxs="i", ylab = "Precip (mm)", xlab = ""))</div><div>par(op)</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-09-10 16:15 GMT-03: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>:<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>Olá pessoal,</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"><span>Estou tentando fazer um gráfico de contorno no R mas ainda não encontrei a função certa para isso. </span><span style="background-color:transparent">Basicamente, eu tenho uma simulação de umidade em 11 camadas do solo
para cada dia do ano, e gostaria de mostrar isso na forma de um perfil de solo.</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 style="background-color:transparent"><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"><span style="background-color:transparent">Por favor leia os comentários ao reproduzir os meus dados com o código abaixo:</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 style="background-color:transparent"><br></span></div><div style="background-color:transparent">library(repmis) # reads text data directly from dropbox - no need to download any file</div><div style="background-color:transparent"><br></div><div style="background-color:transparent"># read data</div><div style="background-color:transparent">url <- '<a href="https://dl.dropboxusercontent.com/u/27700634/precip.txt" target="_blank">https://dl.dropboxusercontent.com/u/27700634/precip.txt</a>'</div><div style="background-color:transparent">tmp <- repmis::source_data(url, sep = '', header = TRUE)</div><div style="background-color:transparent"><br></div><div style="background-color:transparent"># convert julian day to date</div><div style="background-color:transparent">date <- as.Date(tmp$julian, origin='2011-12-31')</div><div style="background-color:transparent">data <-
cbind(date, tmp)</div><div style="background-color:transparent">head(data)</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"><br></div><div style="background-color:transparent"># now, convert soil layers to matrix and transpose it</div><div style="background-color:transparent">mat <- t(as.matrix(data[, 4:14]))</div><div style="background-color:transparent"># essa é a matriz que eu quero plotar. Reparem que ela já está organizada como um "perfil",</div><div style="background-color:transparent"><span style="background-color:transparent"></span></div><div style="background-color:transparent"># com as linhas representando as camadas de solo e as colunas representando os dias do ano</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 style="background-color:transparent"><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"><span style="background-color:transparent">A minha primeira tentativa foi usar a função </span><span style="font-size:12pt">filled.contour</span><span style="font-size:12pt"> do pacote graphics. Eu defino um vetor com o nome das camadas do solo e então eu tento plotar, mas a função diz que a matriz está com as dimensoes invertidas:</span><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"><span style="font-size:12pt"><br></span></div><div style="background-color:transparent"># define vector with depth of soil layers</div><div style="background-color:transparent">depths <- c(0.05,0.10,0.20,0.30,</div><div style="background-color:transparent"> 0.40,0.60,0.80,1.00,</div><div style="background-color:transparent"> 1.50,2.00,2.50)</div><div style="background-color:transparent"><br></div><div style="background-color:transparent"># Plot soil moisture profile</div><div style="background-color:transparent"><span></span></div><div style="background-color:transparent">plot <-
filled.contour(data$julian, depths, mat)</div><div><div><br></div><div>#Error in .filled.contour(x, y, z, levels, col) : dimension mismatch</div><div><br></div></div><div>Obviamente eu consigo transpor a matrix para forçar o gráfico, mas a figura final não é o que eu preciso - o perfil fica de "cabeça para baixo".</div><div><div style="margin-top:0.1em;margin-bottom:0.1em;background-color:transparent">plot <- filled.contour(data$julian, depths, t(mat))</div><div style="margin-top:0.1em;margin-bottom:0.1em;background-color:transparent"><br></div><div style="margin-top:0.1em;margin-bottom:0.1em;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">Como essa função só aceita os eixos em ordem crescente, eu não consigo inverter o eixo y para mostrar a primeira camada do solo no topo do gráfico.</div><div></div></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"><span>Poranto, eu aprecio qualquer sugestão de funções e/ou pacotes para construir o gráfico que eu preciso. </span></div><div></div><div> </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: <a href="tel:%28612%29%20323%209898" value="+16123239898" target="_blank">(612) 323 9898</a></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><br clear="all"><div><br></div>-- <br><div dir="ltr"><div><font style="font-family:arial">###############################################################</font><br><span style="font-family:arial;font-size:small">## </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Jônatan Dupont Tatsch</span></div><div><span style="font-family:arial;font-size:small">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Professor do Departamento de Física</span><font face="arial" style="font-size:small"><br></font></div><div><span style="font-size:small;font-family:arial">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Centro de Ciências Exatas e Naturais (CCNE)</span><br></div><div><span style="font-size:small;font-family:arial">## </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Universidade Federal de Santa Maria</span></div><div><span style="font-family:arial;font-size:small">## </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Santa Maria, RS, Brasil - 97105-900</span><br></div><div style="font-family:arial;font-size:small">## <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Telefone: +55(55)33012083</span><br></div><div style="font-family:arial;font-size:small">## <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px"><a href="http://www.ufsm.br/meteorologia" target="_blank">www.ufsm.br/meteorologia</a></span></div><div style="font-family:arial;font-size:small">###############################################################</div></div>
</div>