Wavefront coding refers to the use of a phase modulating element in conjunction with deconvolution to extend the depth of focus of an imaging system. The coding element is an asymmetrical phase plate shape, for most applications in the form of a trefoil or a cubic polynomial. Phase plates with trefoil shape generate not only the desired amount of trefoil aberration but also spherical aberration. It has been recently shown that a wavefront coding based optical system shows high tolerance to spherical aberration for monochromatic images; however, the depth of focus is considerably shortened for color images. In this work, we will show how to modify the shape of a phase plate in order to optimize its performance for color imaging. The design parameters of the phase plate are obtained by minimizing a merit function by means of genetic algorithms developed for this purpose. The evaluation of the optical characteristics of the phase plates for a feedback with the optimization algorithm is obtained by Zemax. Results will be illustrated by numerical simulations of color images.