В моем понимании поиск отличительных цветов связан с эффективным поиском из единичного куба, где 3 измерения куба - это три вектора по красной, зеленой и синей осям. Это можно упростить для поиска в цилиндре (аналогия HSV), где вы фиксируете насыщенность (S) и значение (V) и находите случайные значения оттенка. Это работает во многих случаях, и посмотрите это здесь:
https://martin.ankerl.com/2009/12/09/how-to-create-random-colors-programmatically/
В R
get_distinct_hues <- function(ncolor,s=0.5,v=0.95,seed=40) {
golden_ratio_conjugate <- 0.618033988749895
set.seed(seed)
h <- runif(1)
H <- vector("numeric",ncolor)
for(i in seq_len(ncolor)) {
h <- (h + golden_ratio_conjugate) %% 1
H[i] <- h
}
hsv(H,s=s,v=v)
}
Альтернативный способ - использовать пакет R «единообразно»
https://cran.r-project.org/web/packages/uniformly/index.html
и эта простая функция может генерировать отличительные цвета:
get_random_distinct_colors <- function(ncolor,seed = 100) {
require(uniformly)
set.seed(seed)
rgb_mat <- runif_in_cube(n=ncolor,d=3,O=rep(0.5,3),r=0.5)
rgb(r=rgb_mat[,1],g=rgb_mat[,2],b=rgb_mat[,3])
}
Можно придумать немного более сложную функцию с помощью поиска по сетке:
get_random_grid_colors <- function(ncolor,seed = 100) {
require(uniformly)
set.seed(seed)
ngrid <- ceiling(ncolor^(1/3))
x <- seq(0,1,length=ngrid+1)[1:ngrid]
dx <- (x[2] - x[1])/2
x <- x + dx
origins <- expand.grid(x,x,x)
nbox <- nrow(origins)
RGB <- vector("numeric",nbox)
for(i in seq_len(nbox)) {
rgb <- runif_in_cube(n=1,d=3,O=as.numeric(origins[i,]),r=dx)
RGB[i] <- rgb(rgb[1,1],rgb[1,2],rgb[1,3])
}
index <- sample(seq(1,nbox),ncolor)
RGB[index]
}
проверьте эти функции:
ncolor <- 20
barplot(rep(1,ncolor),col=get_distinct_hues(ncolor)) # approach 1
barplot(rep(1,ncolor),col=get_random_distinct_colors(ncolor)) # approach 2
barplot(rep(1,ncolor),col=get_random_grid_colors(ncolor)) # approach 3
Однако обратите внимание, что определить отдельную палитру цветов, воспринимаемых человеком, непросто. Какой из вышеперечисленных подходов генерирует разнообразный набор цветов, еще предстоит проверить.
