function color = ColorDistPlot(cb, cr)
% plot the chromatic color distribution of the image

chroma = zeros(256);
cb = round(cb);
cr = round(cr);
for i = 1:length(cb)
   chroma(cb(i), cr(i)) = chroma(cb(i), cr(i)) + 1;
end
surf(chroma)