A method for forming a color transform for a color printing device having a plurality of colorant channels that includes applying a tone scale adjustment function and a total colorant load limit. An optimization process is used to optimize a merit function including a tone scale reproduction accuracy term and a gamut shape term. The method involves determining an initial merit function value responsive to an initial total colorant load limit and an initial set of tone scale adjustment parameters for a tone scale adjustment function. Modified tone scale adjustment parameters and a modified total colorant load limit are determined that provide an improved merit function value subject to a constraint that an average total colorant load is no more than a predefined threshold for a predefined distribution of input colors.