CODEWITHSUNDEEP
matlabcomputer-visionmatlab-cvstcamera-calibration
Pedro77
Pedro77

Reputation: 5294

Understanding extrinsic stereoParameters (rotation) on Matlab stereoCameraCalibrator

Consider the fallowing stereo camera system calibration parameters using matlab stereoCameraCalibrator app.

R1 = stereoParams.CameraParameters1.RotationMatrices(:,:,N);
R2 = stereoParams.CameraParameters2.RotationMatrices(:,:,N);
R12 = stereoParams.RotationOfCamera2;

Were:
R1: rotation from world coordinates (for image N) to camera 1.
R2: rotation from world coordinates (for image N) to camera 2.
R12: rotation from camera 1 coordinates to camera 2. As described on a related SO question

If that is correct, shouldn't R12*R1 == R2 ?

But I'm getting different values, so, what I'm missing here?

Edit

Well, it seams all matrices are transposed. So: R12'*R1' == R2' !
Why they are transposed?

Upvotes: 3

Views: 807

Answers (1)

rayryeng
rayryeng

Reputation: 104545

The reason why they are transposed is due to the fact that when performing geometric transformations between coordinates, MATLAB uses row vectors to perform the transformation whereas column vectors are traditionally used in practice.

In other words, to transform a coordinate from one point to another, you typically perform:

x' = A*x

A would be the transformation matrix and x is a column vector of coordinates. The output x' would be another column vector of coordinates. MATLAB in fact uses a row vector and so if you want to achieve the same effect in multiplication, you must transpose the matrix A (i.e. A^{T}) and pre-multiply by A instead of post-multiplying:

x' = x*A^{T}

Here x would be a row vector and to ensure that weighted combination of rows and columns is correctly accumulated, you must transpose A to maintain the same calculations. However, the shape of the output x' would be a row vector instead.

This can also be verified by transposing the product of two matrices. Specifically, if x' = A*x, then in order to transform the output into a row vector x'^{T}, we must transpose the matrix-vector product:

x'^{T} = (A*x)^{T} = x^{T}*A^{T}

The last statement is a natural property of transposing the product of two matrices. See point 3 at the Transpose Wikipedia article for more details: https://en.wikipedia.org/wiki/Transpose#Properties

The reason why the transpose is performed ultimately stems back to the way MATLAB handles how numbers are aligned in memory. MATLAB is a column-major based language which means that numbers are populated in a matrix column-wise. Therefore, if you were to populate a matrix one element at a time, this would be done in a column-wise fashion and so the coefficients are populated per column instead of per row as we are normally used to, ultimately leading to what we've concluded above.

Therefore, when you have transposed both R12 and R1, this brings back the representation into a row major setting where these matrices were originally column major for ease of MATLAB use. The row major setting thus allows you to use coordinates that are column vectors to facilitate the transformation. This column vector setting is what we are used to. Therefore, multiplying R12 and R1 after you transpose them both brings you to the correct transformation matrix R2 in the standard row major representation.

Upvotes: 3

Related Questions