This paper investigated minimum-norm least-squares solutions to the quaternion matrix equation $KXL+MYN=O$, where the unknown matrices $X$ and $Y$ were subject to bi-Hermitian or skew bi-Hermitian constraints. Building on the derived theoretical results, an explicit numerical algorithm was proposed. Moreover, several numerical examples were presented to validate the accuracy of these results.