Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

I have a matrix equation as follows from a sensor calibration problem. At the heart of this problem is a number of 3x3 rotation matrices:

$ \mathbf{R}_{b1}^{b2} = \mathbf{R}_{v}^{b} \mathbf{R}_{v1}^{v2} (\mathbf{R}_{v}^{b})^T $

where the following is known: $ \mathbf{R}_{v1}^{v2} = \begin{pmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end{pmatrix} $

And the matrix $\mathbf{R}_{b1}^{b2}$ has already been estimated (and is necessarily symmetric)

What I would like to do is to estimate $\mathbf{R}_{b}^{v}$. I'm guessing that there will be no unique solution, but even if it can be constrained to a family of solutions.

Is there anyone that would be able to give a few hints as to how to go about solving this problem?

share|improve this question
add comment

3 Answers 3

up vote 2 down vote accepted

It's taken some time, but I have found a nice solution to the problem.

Utimately, the key is in the rotation matrix: $ \mathbf{R}_{v1}^{v2} = \begin{pmatrix} -1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end{pmatrix} $

As it turns out, the eigenvalues of this rotation matrix is -1, -1 and 1. We know that a similarity transform must preserve the eigenvalues, and therefore, the eigenvalues of $ \mathbf{R}_{b1}^{b2}$ must also be -1, -1 and 1.

Since $ \mathbf{R}_{b1}^{b2}$ is symmetric, we know that there must be an eigendecomposition as follows:

$ \mathbf{A} = \mathbf{Q} \mathbf{\Lambda} \mathbf{Q}^T$.

where $\mathbf{\Lambda}$ is a diagonal matrix of eigenvalues and $\mathbf{Q}$ is an orthogonal matrix constructed from the eigenvectors of $\mathbf{A}$.

In this case $\mathbf{\Lambda} = \mathbf{R}_{v1}^{v2}$ and therefore $\mathbf{Q} = \mathbf{R}_{v}^{b}$.

Hence, calculating the eigenvectors of matrix $\mathbf{R}_{b1}^{b2}$ forms a set of possible solutions. Note that $\mathbf{Q}$ is not necessarily special orthogonal - the eigenvectors may be multipled by -1 as required to get the +1 determinate.

Of course, there will be ambiguity in the solution, but choosing which rotation matrix is correct is problem-specific.

share|improve this answer
    
See also the solution from @AlbertH below, noting the eigenvector associated with the eigenvalue +1. –  Damien Jul 16 '11 at 4:07
add comment

$\mathbf R_{b1}^{b2}$ is a $\pi$ rotation. You can find its axis $a$, determining the eigenvector (with norm $1$) relative to the eigenvalue $1$.
$\mathbf R_v^b$ is the rotation that transforms the $z$ axis into the $a$ axis.

share|improve this answer
add comment

This is an instance of the equation $AX=XB$ (if you right-multiply by $X = {R}_{v}^{b}$), which arises in the so-called hand-eye-calibration problem. You find here a nice list of papers of how to solve this problem.

share|improve this answer
    
Actually, I really like the list of papers, but I have found a much easier way to solve this particular problem. –  Damien Jul 16 '11 at 3:50
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.