Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos espaços no solo
Caros colegas do R,
Um amigo me pediu ajuda para fazer um diagrama volumétrico de algumas propriedades do solo (agua retida a diferentes potenciais, volume de sólidos, etc).
Este é uma figura que ilustra todo a distribuição do espaço poroso de um perfil de solo. Mas não consigo montar este diagrama no R. Talvez alguém aqui da lista possa me dar uma ajuda ou já tenha visto algo semelhante. Ela era usualmente feita como figura no Corel.
A figura/diagrama que comento está neste endereço
https://www.dropbox.com/l/zPSzOqAKPVSEhHFopbu3Xc
No diagrama no eixo x fica o volume do solo (variando de 0,0 - 1,0 m3/m3. No eixo y fica os horizontes/profundidades em cm. E o volume fica dividido entre os valores de i) CAD1 capacidade de água disponível – [0,010 – 1,5 MPA], ii) CAD2 Capacidade de água disponível [ 0,033 – 1,5 MPA], iii) Macroporos; iv) Água fortemente retida; v) água efêmera (conteúdo de água retida entre 0,01 e 0,033 MPA), vi) Matéria sólida. Alguns deste parâmetros se sobrepõe.
Sei que a regra e tentar enviar um CMR. Mas não sei por onde começar.
Cordialmente,
Wenceslau Teixeira
Embrapa Solos
Horiz.
Profund.
CAD1 (0,01)
CAD2 (0,03)
Macro
Água fortemente retida (PM)
Água retida entre (0,01 e 0,03)
Matéria sólida
A
0-18
0,079
0,062
0,311
0,206
0,08
0,396
AB
18-37
0,085
0,057
0,252
0,260
0,016
0,387
Bw1
37-90
0,079
0,051
0,251
0,265
0,017
0,388
Bw2
90-150
0,089
0,054
0,215
0,263
0,019
0,414
-- Wenceslau Teixeira
wgt007@gmail.com skype: wgt007
-- Wenceslau Teixeira wgt007@gmail.com skype: wgt007
A partir dessa tabela que você forneceu, você pode criar as listas que demarcam os vértices dos polígonos que representam cada porção do solo. Algo como plot(x=NULL, y=NULL, xlim=c(0,1), ylim=c(1.2,0), xlab="Umidade (%)", ylab="Profundidade (m)") ## points(expand.grid(seq(0,1,0.1), seq(0,1.2,0.2))) ## l <- locator() ## lines(l) ## l <- lapply(l, round, digits=1) ## dput(l) l1 <- list(x = c(0, 0.2, 0.2, 0.4, 0.4, 0.5, 0.5, 0), y = c(1.2, 1.2, 0.8, 0.8, 0.4, 0.4, 0, 0)) l2 <- list(x = c(0.2, 0.7, 0.7, 0.6, 0.6, 0.6, 0.5, 0.5, 0.4, 0.4, 0.2, 0.2), y = c(1.2, 1.2, 0.8, 0.8, 0.6, 0, 0, 0.4, 0.4, 0.8, 0.8, 1.2)) l3 <- list(x = c(0.7, 1, 1, 0.6, 0.6, 0.7), y = c(1.2, 1.2, 0, 0, 0.8, 0.8)) polygon(l1$x, l1$y, col="gray90") polygon(l2$x, l2$y, density=3, angle=33) polygon(l3$x, l3$y, density=12, angle=55, col="blue") À disposição. Walmes.
Vale a pena dar uma olhada também no pacote plotrix (http://cran.r-project.org/web/packages/plotrix/index.html). Ele tem muitos gráficos aplicados a solos: física, química, etc Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
participantes (3)
-
Jose Claudio Faria -
walmes . -
Wenceslau Teixeira