3D Rotation Matrix

Hi all, I'm new here and am trying to obtain a 3D rotation matrix.

Currently, I have a set of n points: (x1,y1,z1), (x2,y2,z2)...(xn,yn,zn)

from which I have designated one particular point as the origin, and obtained an x'- and y'- vector, based on some features. It is unimportant, for this question, how I chose my origin and got the x'- and y' vector. From the x'- and y' vector, I can also obtain the z'-vector by taking cross product of x' and y'.

Now, I would like to rotate all these points about the origin I have defined, such that the x'-, y'- and z'- vectors are in line with the x-, y- and z- axes respectively.

Hence, I am looking for a 3x3 rotation matrix, R, that can be applied to all points such that:

Transformed (x,y,z) = R * Original (x,y,z)

An example of the vectors to be mapped would be say:

x': [-0.2831, -0.9246, 0.2548]' to be mapped to x: [1 0 0]

y': [0.9242, -0.1919, 0.3303]' to be mapped to y: [0 1 0]

z' will be mapped automatically if the above 2 are satisfied.

Thank you for your help!

[Walter Roberson]

Duplicates https://www.mathworks.com/matlabcentral/answers/349968-3d-rotation-matrix but this one appears to have been updated so I have closed that.

Answer 1

[Honglei Chen]

I'm not sure if I fully follow your question, but here is my understanding. The first step is to find the coordinates of all points in the new coordinate system. For that you'll have to first subtract origin from all points and then convert all points to the local coordinate system. That conversion matrix should just be R1 = [x' y' z']'. If you have access to Phased Array System Toolbox, you can consider using the following function

https://www.mathworks.com/help/phased/ref/global2localcoord.html

You can then rotate the points in the new coordinate system using Euler matrices. If you want again you can use the following functions:

https://www.mathworks.com/help/phased/ref/rotx.html

https://www.mathworks.com/help/phased/ref/roty.html

https://www.mathworks.com/help/phased/ref/rotz.html

HTH

Answer 2

[Herve Hugonnet]

The simplest if you do not have the phase array system toolbox is to define anonymous functions :

rotx = @(t) [1 0 0; 0 cos(t) -sin(t) ; 0 sin(t) cos(t)] ;
roty = @(t) [cos(t) 0 sin(t) ; 0 1 0 ; -sin(t) 0  cos(t)] ;
rotz = @(t) [cos(t) -sin(t) 0 ; sin(t) cos(t) 0 ; 0 0 1] ;

[Alec]

Thanks for this