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.