1. A method of optimizing gray component replacement in a color image, the method: comprising: receiving a color image from an input imaging device; determining a color gamut of an output imaging device; generating a CIELab to CMY numerical model to convert the color image to a color space within the color gamut of the output imaging device; generating a CMYK to CIELab numerical conversion model; removing a portion of CMY color values generated from the CIELab to CMY numerical conversion model; replacing the removed portion of the CMY color values with black color values; and optimizing CMYK color values via a multidimensional optimization using the CMYK to CIELab numerical model.