How do i generate a rotation matrix iteratively.

Question

Stewart Tan 2019-9-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479834-how-do-i-generate-a-rotation-matrix-iteratively

评论： Walter Roberson 2019-9-11

I have two matrices containing coordinates:

example1 = [1 3; 
            4 5; 
            2 3; 
            6 17];
        
example2 = [1 4; 
            6 2; 
            8 9; 
            10 11];

and I'm implementing the RANSAC algorithm to remove outlier coordinates. To do that, I will need to include a factor of whether a not one coordinate in example1 is either a rotation, scaling or translation which result in the coordinate in example2. (1 3 -> 1 4, 4 5 -> 6 2 etc..) Hence, I will need to have a rotation matrix, scaling and translation matrix to be built iteratively within the RANSAC algorithm.

For now, I'll just keep it simple by just using a rotation matrix.

How do i generate a rotation matrix [cos θ -sin θ; sin θ cos θ] if i do not have anything else other than the two coordinate matrices above?

Edit: sorry for the confusion but the goal is to "estimate" the transformation matrix in every N iteration. The following shows:

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2019-9-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479834-how-do-i-generate-a-rotation-matrix-iteratively#answer_391324

在 MATLAB Online 中打开

example1 = [1 3; 
            4 5; 
            2 3; 
            6 17];
        
example2 = [1 4; 
            6 2; 
            8 9; 
            10 11];
        
syms theta
RM = [cos(theta) -sin(theta); sin(theta) cos(theta)];
residue = simplify( sum(sum((example2 - (RM * example1.').').^2)) );
best_theta = vpasolve(diff(residue,theta));

It is also possible to code it as a numeric minimization of residue over a bounded range, without using the symbolic toolbox.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2019-9-11

Your original goal was restricted to rotation matrix, and my code finds the rotation matrix that produces the least squared error of rotation between the given sets of coordinates.