The RGB2GRAY conversion model is the most popular and classical tool for image decolorization. A recent study showed that adapting the three weighting parameters in this first-order linear model with a discrete searching solver has a great potential in its conversion ability. In this paper, we present a two-step strategy to efficiently extend the parameter searching solver to a two-order multivariance polynomial model, as a sum of three subspaces. We show that the first subspace in the two-order model is the most important and the second one can be seen as a refinement. In the first stage of our model, the gradient correlation similarity (Gcs) measure is used on the first subspace to obtain an immediate grayed image. Then, Gcs is applied again to select the optimal result from the immediate grayed image plus the second subspace-induced candidate images. Experimental results show the advantages of the proposed approach in terms of quantitative evaluation, qualitative evaluation, and algorithm complexity.