如何获得平面上三个不共线的点？ - C++

我正在尝试实现线平面相交算法。根据维基百科我需要平面上的三个非共线点才能做到这一点。

因此我尝试实施这个算法然而在 C++ 中。肯定有问题，因为我可以选择任何 x 和 y 坐标并且它们适合平面，这是没有意义的。如果平面垂直且沿 x 轴怎么办？ y=1 的点将不在平面内。

我意识到这个问题已经在 StackOverflow 上发布了很多，并且我看到了很多解决方案，其中平面是由 3 个点定义的。但我只有正常和地位。在我整理出我的非共线点查找器之前，我无法测试我的线平面相交算法。

现在的问题是，我要除以normal.z，而当normal.z为0时，这显然不起作用。

我正在用这个平面进行测试： Plane* p = new Plane(Color(), Vec3d(0.0,0.0,0.0), Vec3d(0.0,1.0,0.0)); // 第二个参数：位置，第三个参数：正常

当前的代码给出了这个错误的答案：

{0 , 0 , 0} // alright, this is the original
{12.8377 , 17.2728 , -inf} // obviously this is not a non-colinear point on the given plane

这是我的代码：

std::vector<Vec3d>* Plane::getThreeNonColinearPoints() {
    std::vector<Vec3d>* v = new std::vector<Vec3d>();

    v->push_back(Vec3d(position.x, position.y, position.z)); // original position can serve as one of the three non-colinear points.

    srandom(time(NULL));

    double rx, ry, rz, start;

    rx = Plane::fRand(10.0, 20.0);
    ry = Plane::fRand(10.0, 20.0);
    // Formula from here: http://en.wikipedia.org/wiki/Plane_(geometry)#Definition_with_a_point_and_a_normal_vector
    // nx(x-x0) + ny(y-y0) + nz(z-z0) = 0
    // |-----------------| <- this is "start"
    //I'll try to insert position as x0,y0,z0 and normal as nx,ny,nz, and solve the equation
    start = normal.x * (rx - position.x) + normal.y * (ry - position.y);
    // nz(z-z0) = -start
    start = -start;
    // (z-z0) = start/nz
    start /= normal.z; // division by zero
    // z = start+z0
    start += position.z;
    rz = start;

    v->push_back(Vec3d(rx, ry, rz));

    // TODO one more point

    return v;
}

我意识到我可能试图完全错误地解决这个问题。如果是这样，请链接一个具体的实现。当我看到这么多线平面相交的实现时，我确信它一定存在。

提前致谢。

平面可以用多种方式定义。通常使用平面上的点和法向量。从三个点获取法线向量 (P1, P2, P3) 取三角形边的叉积

P1 = {x1, y1, z1};
P2 = {x2, y2, z2};
P3 = {x3, y3, z3};

N = UNIT( CROSS( P2-P1, P3-P1 ) );
Plane P = { P1, N }

反过来，从一个点开始P1和正常的N到三点，你从任何方向开始G not沿着正常的N这样DOT(G,N)!=0。那么沿着平面的两个正交方向是

//try G={0,0,1} or {0,1,0} or {1,0,0}
G = {0,0,1};
if( MAG(CROSS(G,N))<TINY ) { G = {0,1,0}; }
if( MAG(CROSS(G,N))<TINY ) { G = {1,0,0}; }
U = UNIT( CROSS(N, G) );  
V = CROSS(U,N);
P2 = P1 + U;
P3 = P1 + V;

线由点和方向定义。通常有两点（Q1, Q2) 定义线

Q1 = {x1, y1, z1};
Q2 = {x2, y2, z2};
E = UNIT( Q2-Q1 );
Line L = { Q1, E }

直线和平面的交点由直线上的点定义r=Q1+t*E与平面相交，使得DOT(r-P1,N)=0。这已解决标量距离t沿着线作为

t = DOT(P1-Q1,N)/DOT(E,N);

和位置为

r = Q1+(t*E);

注意：DOT()返回两个向量的点积，CROSS()叉积，以及UNIT()单位向量（幅度=1）。

DOT(P,Q) = P[0]*Q[0]+P[1]*Q[1]+P[2]*Q[2];
CROSS(P,Q) = { P[1]*Q[2]-P[2]*Q[1], P[2]*Q[0]-P[0]*Q[2], P[0]*Q[1]-P[1]*Q[0] };
UNIT(P) = {P[0]/sqrt(DOT(P,P)), P[1]/sqrt(DOT(P,P)), P[2]/sqrt(DOT(P,P))};
t*P =  { t*P[0], t*P[1], t*P[2] };
MAG(P) = sqrt(P[0]*P[0]+P[1]*P[1]+P[2]*P[2]);