function color = ColorDistPlot(cb, cr)% plot the chromatic color distribution of the imagechroma = zeros(256);cb = round(cb);cr = round(cr);for i = 1:length(cb)   chroma(cb(i), cr(i)) = chroma(cb(i), cr(i)) + 1;endsurf(chroma)