sharpshootR: A Soil Survey Toolkit

library(aqp)
library(venn)
library(sharpshootR)

venn(7, box = FALSE)

x <- as.data.frame(matrix(sample(0:1, 150, replace = TRUE), ncol = 5))
names(x) <- c('DB', 'DSM', 'DSS', 'Initial', 'KSSL')
venn(x)


x <- as.data.frame(matrix(sample(0:1, 1500, replace = TRUE), ncol = 5))
names(x) <- c('Life', 'Work', 'Hobbies', 'Family', 'Sleep')

venn(x, ellipse = TRUE, zcolor = 'style', box = FALSE)
venn(x, ellipse = FALSE, zcolor = 'style')

venn(5, ilabels = TRUE, zcolor = "style", ellipse = TRUE)


# chips <- c('10YR 8/2', '10R 4/8', '5Y 8/8', '5B 4/8', '5GY 4/8', '5Y 3/1', '7.5P 4/8')
# chips <- c('10P 4/8', '10R 4/8', '5Y 8/8', '5B 4/8', '5GY 4/8', '5Y 3/1')
# chips <- c('10YR 8/2', '10R 4/8', '5Y 8/8', '5B 4/8', '5GY 4/8')
chips <- c('10YR 8/2', '10R 4/8', '5Y 8/8', '5B 4/8')
# chips <- c('2.5R 4/6', '2.5B 4/6', '2.5Y 8/2')
# chips <- c('10R 6/8', '5GY 4/4')

# # reference diagram
# venn(4, zcolor = 'bw', box = FALSE, ellipse = TRUE, ilabels = TRUE)

colorMixtureVenn(chips)
colorMixtureVenn(chips, mixingMethod = 'exact')
colorMixtureVenn(chips, labels = FALSE)

colorMixtureVenn(chips, ellipse = TRUE)

colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'exact')
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'reference')
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'estimate')


par(mfcol = c(1,2))
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'reference')
title('Reference Spectra', line = -3)
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'exact')
title('Exact', line = -3)

mixMunsell(c('10R 4/8', '5B 4/8'), mixingMethod = 'reference')
mixMunsell(c('10R 4/8', '5B 4/8'), mixingMethod = 'exact')

par(bg = 'black', fg = 'white')
colorMixtureVenn(chips, ellipse = TRUE)


par(mfcol = c(1,2))
chips <- c('10YR 8/1', '2.5YR 3/6', '10YR 2/2')

colorMixtureVenn(chips, mixingMethod = 'reference')
title('Reference', line = -3)

colorMixtureVenn(chips, mixingMethod = 'exact')
title('Exact', line = -3)


colorMixtureVenn(chips)
title('Simulation of Subtractive Mixture', line = -4)

colorMixtureVenn(chips, mixingMethod = 'estimate')
title('Weighted Mean CIELAB', line = -4)



chips <- c('10YR 10/1', '2.5YR 3/6', '10YR 1/1')
colorMixtureVenn(chips)
title('Simulation of Subtractive Mixture', line = -4)

colorMixtureVenn(chips, mixingMethod = 'estimate')
title('Weighted Mean CIELAB', line = -4)



par(mfcol = c(1,2))
chips <- c('10YR 8/1', '5Y 8/8', '5B 4/8')

colorMixtureVenn(chips, mixingMethod = 'reference')
title('Simulation of Subtractive Mixture\nReference', line = -3)

colorMixtureVenn(chips, mixingMethod = 'exact')
title('Simulation of Subtractive Mixture\nExact', line = -3)

colorMixtureVenn(chips, mixingMethod = 'exact')
title('Simulation of Subtractive Mixture\nExact', line = -3)

colorMixtureVenn(chips, mixingMethod = 'estimate')
title('Weighted Mean CIELAB', line = -3)




chips <- c('10YR 8/1', '2.5YR 3/6', '10YR 8/10', '5G 4/8')
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'exact')

chips <- c('2.5YR 8/2', '2.5B 3/6', '2.5G 4/4', '2.5R 4/8', '2.5Y 6/10')
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'exact')

chips <- c('10YR 8/2', '5G 4/8', '2.5Y 8/2', '5R 4/8')
colorMixtureVenn(chips, ellipse = TRUE, mixingMethod = 'exact')

library(soilDB)
x <- fetchOSD('musick')

chips <- unique(sprintf("%s %s/%s", x$hue, x$value, x$chroma))

colorMixtureVenn(chips[1:5], mixingMethod = 'exact', ellipse = TRUE, labels = FALSE)
# colorMixtureVenn(chips[1:5], mixingMethod = 'reference')