1. A method, comprising: for each of a plurality of initial colors, estimating color space coordinates, including a lightness value (L*), a red-green value (a*), and a yellow-blue value (b*), that correspond to that initial color; for each of the estimated color space coordinates, converting the red-green value (a*) and the yellow-blue value (b*) into a chroma value and a hue angle, wherein converting the red-green value (a*) and the yellow-blue value (b*) into a chroma value comprises calculating the chroma value as the square root of the sum of the square of the red-green value (a*) and the square of the yellow-blue value (b*); applying a desaturation factor to each of the chroma values to generate corresponding desaturated chroma values; for each of the desaturated chroma values, converting that desaturated chroma value and the associated hue angle into a desaturated red-green value (a*) and a desaturated yellow-blue value (b*); converting each of the desaturated red-green values (a*), the desaturated yellow-blue values (b*), and the associated lightness values (L*) into a respective one of a plurality of desaturated colors; generating a color correction matrix relating the initial colors to the desaturated colors; and with image processing circuitry, applying the color correction matrix to at least one image to reduce saturation levels of the image without substantially altering hues in the image.