The function MAXDISTCOLOR generates an RGB colormap of maximally distinct colors. Its optional input arguments are explained in this illustrated document.

Although it is tempting to create colors which are extremely distinct, such sets of colors are garish and visually complex. Restricting the lightness and chroma ranges provides much more pleasant sets of distinct colors: see the end of this document for an example of this.

Also note different monitors and screens show colors very differently, so without a calibrated monitor it is worth comparing several monitors.


Getting Started

The simplest usage is to call MAXDISTCOLOR with the number of colors and the SRGB_TO_LAB colorspace conversion function (which is distributed with MAXDISTCOLOR and converts from sRGB to the non-uniform, inadequate, outdated L*a*b* colorspace):

N = 10;
rgb = maxdistcolor(N,@srgb_to_Lab)
X = linspace(0,10,1000);
Y = bsxfun(@(x,n)n*sin(x+2*n*pi/N), X(:), 1:N);
axes('ColorOrder',rgb, 'NextPlot','replacechildren', 'XTick',[],'YTick',[])
plot(X,Y, 'linewidth',4)
rgb =

         0         0    1.0000
         0    1.0000         0
    1.0000         0         0
    1.0000         0    0.6190
    1.0000    0.9606         0
         0    0.6220    1.0000
         0    0.6535    0.3810
    0.5556    0.3780    0.0159
         0         0    0.4127
         0    0.1102    0.1270

Recommended Uniform Colorspace: CAM02

The color distinctiveness depends significantly on the quality of the colorspace used: it is strongly recommended to use a uniform colorspace, such as CAM02-UCS or CAM02-LCD (which are derived from the color appearance model CIECAM02). You can download my implementation of CIECAM02 and use the function SRGB_TO_JAB with the recommended option 'LCD' (for very dissimilar colors) or the default 'UCS' (for closer colors):

fun = @(m)srgb_to_Jab(m,true,'LCD'); % recommended = Large Color Difference.
rgb = maxdistcolor(N,fun)
axes('ColorOrder',rgb, 'NextPlot','replacechildren', 'XTick',[],'YTick',[])
plot(X,Y, 'linewidth',4)
rgb =

    0.8095         0         0
    0.6190    0.0630    1.0000
         0    0.5512         0
    0.5873    1.0000         0
         0    0.6693    0.9206
    0.4286    0.3307    0.3333
    1.0000    0.4803    0.7619
         0    0.0157    0.6349
    1.0000    0.6299    0.0317
         0    0.0079         0

Specify Excluded Colors

Option exc excludes color/s from the output colormap. This is useful for specifying the background color of the plot (thus ensuring that the plotted colors are distinguishable from the background color), or for adding new distinct data to an already existing plot. By default the color white is excluded (i.e. exc=[1,1,1]), so that the output colors are visible on the white axes background or when printed on white paper.

rgb = maxdistcolor(N,fun, 'exc',[0,0,0]) % Exclude black (e.g. background).
axes('ColorOrder',rgb, 'NextPlot','replacechildren', 'XTick',[],'YTick',[], 'Color',[0,0,0])
plot(X,Y, 'linewidth',4)
rgb =

    1.0000    1.0000         0
    0.8571    0.3780    1.0000
    0.0317    0.7638    1.0000
    0.6032    0.3622         0
    0.3492         0    1.0000
         0    0.6772    0.0159
    0.4603         0    0.3651
    1.0000    0.8425    0.8571
         0    0.3780    0.3810
    1.0000         0    0.3651

Specify Included Colors

Option inc includes color/s in the output colormap: for example it can be used to find a set of distinct colors that includes a corporate or university colorscheme. By default no specific colors are included.

Note that both options inc and exc can be provided as:

rgb = maxdistcolor(N,fun, 'inc',[1,0,0]) % Include red.
axes('ColorOrder',rgb, 'NextPlot','replacechildren','XTick',[],'YTick',[])
plot(X,Y, 'linewidth',4)
rgb =

    1.0000         0         0
    0.5873    0.0472    1.0000
         0    0.4882         0
         0    0.6378    0.7937
         0    1.0000         0
    0.5079         0    0.2698
         0         0    0.5397
    1.0000    0.4803    0.8095
    0.8730    0.6299         0
    0.0952    0.1024         0

Limit Lightness Range

To create attractive sets of colors, or to match document requirements, it can be useful to limit the lightness range of the output colors. The output lightness range is controlled using options Lmin and Lmax. Note that the lightness limits are scaled so that 0==black and 1==white.

N = 9;
[rgb,ucs] = maxdistcolor(N,fun, 'Lmin',0.4, 'Lmax',0.6);
scatter3(ucs(:,3),ucs(:,2),ucs(:,1), 256, rgb, 'filled')
% Plot outline of RGB cube:
M = 23;
[X,Y,Z] = ndgrid(linspace(0,1,M),0:1,0:1);
mat = fun([X(:),Y(:),Z(:);Y(:),Z(:),X(:);Z(:),X(:),Y(:)]);
J = reshape(mat(:,1),M,[]);
a = reshape(mat(:,2),M,[]);
b = reshape(mat(:,3),M,[]);
% Add Labels:
title('Colors with Limited Lightness, Inside RGB Cube')
xlabel('Hue Dim.2 (b'')')
ylabel('Hue Dim.1 (a'')')
zlabel('Lightness (J'')')

Limit Chroma Range

To create attractive sets of colors, or to match document requirements, it can be useful to limit the chroma range of the output colors. The output chroma range is controlled using options Cmin and Cmax. Note that the chroma limits are scaled so that 1==max(Chroma).

[rgb,ucs] = maxdistcolor(N,fun, 'Cmin',0.5, 'Cmax',0.6);
scatter3(ucs(:,3),ucs(:,2),ucs(:,1), 256, rgb, 'filled')
% Plot outline of RGB cube:
% Add Labels:
title('Colors with Limited Chroma, Inside RGB Cube')
xlabel('Hue Dim.2 (b'')')
ylabel('Hue Dim.1 (a'')')
zlabel('Lightness (J'')')

RGB Bit Depth

By default MAXDISTCOLOR creates the RGB gamut using [6,7,6] bits for the red, green, and blue channels. The bit options can be used to:

Specifying 2 bits per channel defines an RGB gamut with 64 colors, so requesting 64 colors will return every color from that RGB gamut:

rgb = maxdistcolor(64,@srgb_to_Jab, 'bitR',2,'bitG',2,'bitB',2, 'exc',[]);
scatter3(rgb(:,3),rgb(:,2),rgb(:,1), 256, rgb, 'filled')

Sort the Colormap

The greedy algorithm returns an undefined, non-random color order. Some color orders that might be more attractive or useful can be selected using the sort option: see the m-file help for the complete list of sort orders. Note that some orders use an algorithm that checks all permutations: those orders will only work for nine or fewer colors.

The example here uses 'maxmin', which maximizes the minimum adjacent color difference, so that adjacent colors are maximally dissimilar:

rgb = maxdistcolor(N,fun, 'sort','maxmin', 'exc',[0,0,0;1,1,1]);
title('Sorted Colors')
ylabel('Colormap Index')

Bonus: Color Names

Sometimes it can be useful to give colors names, e.g. when providing a list of colors for a user to select from, or for defining HTML colors. One easy way to find a color name is to download my FEX submission COLORNAMES:

text(ones(1,N), 1:N, colornames('CSS',rgb),...
	'HorizontalAlignment','center', 'BackgroundColor','white')
title('Colors with Color Names')
ylabel('Colormap Index')

Sort Path Open or Closed

For sort types 'maxmin', 'minmax', 'longest', and 'shortest' the default sort algorithm treats the end colors as being disconnected from each other, forming an 'open' path. Selecting the path option as 'closed' treats the end colors as being adjacent to each other, forming one continuous loop (e.g. for data using modulo arithmetic):

N = 7;
txt = cellstr(num2str((1:N).'));
opt = {'Lmin',0.42, 'Lmax',0.57, 'bitR',7,'bitG',8,'bitB',7};
% Path Closed:
[~,ucs] = maxdistcolor(N,fun, 'sort','shortest', 'path','closed', opt{:});
idx = [1:N,1]; % closed loop
plot3(ucs(idx,3)-1,ucs(idx,2)-1,ucs(idx,1)-1,'-b', 'linewidth',2)
text(ucs(:,3),ucs(:,2),ucs(:,1)-9, txt, 'Color','b', 'HorizontalAlignment','center')
% Path Open:
[rgb,ucs] = maxdistcolor(N,fun, 'sort','shortest', 'path','open', opt{:});
idx = 1:N; % open loop
plot3(ucs(idx,3)+1,ucs(idx,2)+1,ucs(idx,1)+1,'-r', 'linewidth',2)
text(ucs(:,3),ucs(:,2),ucs(:,1)+9, txt, 'Color','r', 'HorizontalAlignment','center')
% Legend and Colors:
legend({'Closed','Open'}, 'Location','NorthWest')

Visualize Color Distances

It is easy to plot the distances between the returned color nodes, providing a visual confirmation that the colors are maximally distinct. The vertical gap between the mono-colored points (along the bottom) and the bi-colored points indicates the closest distance of the colors in the colorspace (this distance is maximized by the repeated greedy algorithm):

for k = 1:N
	dst = sqrt(sum(bsxfun(@minus,ucs,ucs(k,:)).^2,2));
	scatter(1:N, dst, 123, rgb,...
		'MarkerFaceColor',rgb(k,:), 'LineWidth',2.8, 'Marker','o')
	hold on
title('Euclidean Distances in CAM02-LCD Colorspace')
ylabel('Euclidean Distance')
xlabel('Colormap Index')

Bonus: Interactive Tool

An easy way to vary the optional arguments and see how they affect the output colormap is to use the included function MAXDISTCOLOR_VIEW. Simply adjust the sliders, menu options, and the included/excluded colors to create new sets of colors. For example, many attractive color sets can be found by limiting the lightness and chroma to quite tight ranges:

maxdistcolor_view(7,fun, 'Cmin',0.2,'Cmax',0.3, 'Lmin',0.7,'Lmax',0.8, 'sort','lightness')