The best approximation in the optimal solution set of the Hankel-norm model reduction problem is studied in this paper since the optimal solutions are not unique for linear multi-input-multi-output systems (matrix-value transfer functions). This kind of model reduction problems is defined properly and intuitively. The appropriate Schmidt pairs will characterize the sub-layers of the optimal model errors. The optimal solution set is also parametrized in the suitable domain in order to keep the reduced model with the constant order after the optimizations. The results from this proposed approach show that they are better than those from the other optimal approximate models in sense of the Hankel operator singular values.