# two overlapping horizons
z <- data.frame(
id = 'SPC',
top = c(0, 25, 25, 50, 75, 100, 100),
bottom = c(25, 50, 50, 75, 100, 125, 125)
)
# init SPC
depths(z) <- id ~ top + bottom
# flag perfectly overlapping horizons
z$.overlapFlag <- flagOverlappingHz(z)
# thematic sketches
plotSPC(z, color = '.overlapFlag', hz.depths = TRUE,
depth.axis = FALSE, cex.names = 0.85)
Run the code above in your browser using DataLab