Rotate the coordinate system to align an existing plane with Y'Z' plane

Question

PASUNURU SAI VINEETH 2022-5-15

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1719515-rotate-the-coordinate-system-to-align-an-existing-plane-with-y-z-plane

评论： PASUNURU SAI VINEETH 2022-5-17

I have a plane z = (-47.407313)*x + (0.322175)*y + (-3.333979) and data points of a ball flight trajectory on it. Even though the best fit is known to be a parabola, since the plane in discussion doesn't align with any of the XY, YZ and ZX planes, it is infeasible to define an explicit relationship z = F(x,y) for obtaining such a curve using curve fitting. Now, my only idea is to rotate the coordinate system so that the plane is parallel to Y'Z' where I can apply (z')=a(y')^2+b(y')+c. After finding the best fit and sampling the points on the fit, I need to retract the coordinate system to get the actual coordinates. Even though the idea seems clear, the application in MATLAB doesn't and I'm hoping someone can guide me with the rotation.

I would be grateful if anyone can suggest an alternate, easier approach to get the 3D best fit equation without this hassle.

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

PASUNURU SAI VINEETH 2022-5-16

编辑：PASUNURU SAI VINEETH 2022-5-16

One of the answers helped me. After obtaining the curve fit and sampling, I transformed them back to the original coordinate system where I could apply equations of motion.

Btw, do you know of a way I can directly fit a 3D curve (not the surface or a spline along the surface that isn't smooth)? The process is straightforward for 2D and I couldn't do the 3D fit even after projecting all the points to the best fit plane (angled to the XY, YZ and ZX). I am looking to eliminate this hassle of plane projection and double transformation.

Matt J 2022-5-16

You would need a parametric model for the 3D trajectory (x(t), y(t),z(t))

请先登录，再进行评论。

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

Answer 1

Torsten 2022-5-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1719515-rotate-the-coordinate-system-to-align-an-existing-plane-with-y-z-plane#answer_964520

编辑：Torsten 2022-5-15

https://stackoverflow.com/questions/9423621/3d-rotations-of-a-plane

But are you sure that the parabola in the YZ plane is not also rotated ?

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

PASUNURU SAI VINEETH 2022-5-16

Okay, I understand your point now. Yeah the parabola need not be aligned with the Z-axis. It somehow worked for this case but do you suggest a general formula ay^2+bz^2+cyz+dy+ez+f=0 when the direction of directrix is arbritary?

Torsten 2022-5-16

在 MATLAB Online 中打开

If you have a 3d-model of the parabola of the form

(a*t^2+b*t+c;d*t^2+e*t+f;g*t^2+h*t+i)

you can directly fit the 9 parameters in question.

请先登录，再进行评论。

Answer 2

Bruno Luong 2022-5-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1719515-rotate-the-coordinate-system-to-align-an-existing-plane-with-y-z-plane#answer_964655

Take a look at this thread and see if you can use it

2 个评论
显示无隐藏无

PASUNURU SAI VINEETH 2022-5-16

Appreciate your response. The thread Torsten shared helped me.

Btw, do you know of a way I can directly fit a 3D curve (not the surface or a spline along the surface that isn't smooth)? The process is straightforward for 2D and I couldn't do the 3D fit even after projecting all the points to the best fit plane (angled to the XY, YZ and ZX). I am looking to eliminate this hassle of plane projection and double transformation.

Bruno Luong 2022-5-16

The answer I gave in the above thread shows to use SVD to fit a point clouds, the input point cloud can be a non degenerated curve (i.e. not a line) and not necessary a plane.

请先登录，再进行评论。

Answer 3

Matt J 2022-5-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1719515-rotate-the-coordinate-system-to-align-an-existing-plane-with-y-z-plane#answer_964720

编辑：Matt J 2022-5-16

在 MATLAB Online 中打开

If you didn't obtain your plane fit with planarFit() from,

https://www.mathworks.com/matlabcentral/fileexchange/87584-object-oriented-tools-for-fitting-conics-and-quadrics

then I suggest you do so. If you did use planarFit(), then the rotation is easy:

pfit=planarFit(xyz);
yz_prime = pfit.R(:,2:3)'*xyz;

However, once you've rotated the data this way, you still can't be sure in which direction gravity pulls (see also my comment above).

One solution would be to post-optimize the in-plane rotation of yz

theta=fminbnd(@(theta) postfun(theta,yz_prime),  0,pi);
[~,abc]=postfun(theta,yz_prime)

using the 1D objective function,

function [fval,abc]=postfun(theta,yz_prime)
 R=makehgtform('zrotate',theta);
 
 yz=R(1:2,1:2)*yz_prime;
 
 [abc,S]=polyfit(yz(:,1),yz(:,2),2)
 fval=S.normr;
 
end