calculate the normal for a plane passing through more than three points

Question

Alberto Acri 2023-9-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027374-calculate-the-normal-for-a-plane-passing-through-more-than-three-points

编辑： Matt J 2023-9-30

采纳的回答： Bruno Luong

V1_in.mat

在 MATLAB Online 中打开

Hi! I have a plane passing through three points (triangular plane) and from which I determine the normal.

P = [ 0.034488  0.036484  0.006912; ...
      0.019104  0.055041  0.047894; ...
     -0.008596  0.123650  0.033786];
N = cross(P(1,:) - P(2,:), P(3,:) - P(2,:));
N = N/norm(N);

Is it possible to do the same thing with a circular plane ('V1_in')?

If it is not possible to do this, how can I create the matrix equal to 'P' using the 'V1_in' plane?

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

请先登录，再进行评论。

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

Answer 1

Bruno Luong 2023-9-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027374-calculate-the-normal-for-a-plane-passing-through-more-than-three-points#answer_1322384

在 MATLAB Online 中打开

V1_in.mat

load('V1_in.mat')
[~,~,V]=svd( V_1-mean(V_1));
N=V(:,3)
N = 3×1
    0.2147
   -0.1801
    0.9599

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

Bruno Luong 2023-9-29

在 MATLAB Online 中打开

% Check againts other method for 3 points

P = [ 0.034488  0.036484  0.006912; ...
      0.019104  0.055041  0.047894; ...
     -0.008596  0.123650  0.033786];
N = cross(P(1,:) - P(2,:), P(3,:) - P(2,:));
N = N/norm(N)
N = 1×3
    0.9037    0.3976    0.1592
[~,~,V]=svd( P-mean(P));
N=V(:,3)% opposite sign, which is arbitrary choice for a 2D plane
N = 3×1
   -0.9037
   -0.3976
   -0.1592

请先登录，再进行评论。

Answer 2

Torsten 2023-9-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027374-calculate-the-normal-for-a-plane-passing-through-more-than-three-points#answer_1322389

编辑：Torsten 2023-9-29

If the points in V1_in.mat all lie in a common plane, you can arbitrarily pick three points and define these points as P.

If the points are only "approximately" in a common plane, you have to determine this plane via regression (like @Bruno Luong did).

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

请先登录，再进行评论。

Answer 3

Matthew Blomquist 2023-9-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027374-calculate-the-normal-for-a-plane-passing-through-more-than-three-points#answer_1322394

Hi Alberto,

If you think of the circular plane as a collection of points with coordinates (x, y, z), you can create a bunch of "triangular planes" by selecting any three coordinates to create a triangle. If all of the coordinates lie in the same plane, computing the normal of any triangular plane will be the same as computing the normal of the circular plane.

So, you can use the same code, but substitute three points from the circular plane (in V1_in) to your P matrix. Also, I often use the functions "patch" and "quiver3" to help visualize the normals.

Hope that helps!

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

请先登录，再进行评论。

Answer 4

Matt J 2023-9-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2027374-calculate-the-normal-for-a-plane-passing-through-more-than-three-points#answer_1322514

编辑：Matt J 2023-9-30

在 MATLAB Online 中打开

Use planarFit() from this FEX package,

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

fitObj=planarFit(V_1')

fitObj =

planarFit with properties:

normal: [0.2147 -0.1801 0.9599]

distance: 82.5216

[hL,hD]=plot(fitObj); %Visualize the fit
legend([hL,hD],'Plane fit', 'XYZ samples', 'Location','northoutside','FontSize',15);
axis padded; 
view(-70,15);drawnow