First of all, by characterizing solutions of the quaternion equation a x = x ? b , this paper studies consimilarity of quaternions and some related consequences. For the important role of coneigenvalues in consimilarity tranformations of quaternion matrices, this paper further derives the relations between principle right coneigenvalues of a quaternion matrix and eigenvalues of the corresponding real representation matrix. Then, based on the real representation matrix, an effective algorithm is presented to calculate all coneigenvalues and the associated coneigenvectors of a quaternion matrix. Finally, two numerical examples are given to verify the effectiveness of the proposed algorithm.