C++ 如何到达交叉点？射线三角形交点C++；

当光线穿过三角形时，我想得到它的交点。我按照在线教程制作了一个函数，但我无法获得交点的正确坐标

例如，如果使用光线原点（0，2，-1）、光线方向（0，1，0）、三角形顶点p0（0，3，0）、p1（-0.5，3，-1）、p2（0.5，3，-1），则应获取交点（0，3，-1）。但是，我得到了交点（0，7，-1），这是不正确的

感谢您的时间和帮助。：-）

float kEpsilon=0.000001；
V3f交叉积（V3f点1、V3f点2）{
V3f矢量；
vector.x=point1.y*point2.z-point2.y*point1.z；
vector.y=point2.x*point1.z-point1.x*point2.z；
vector.z=point1.x*point2.y-point1.y*point2.x；
返回向量；
}
浮动点积（V3f点1、V3f点2）{
浮点点=dot1.x*dot2.x+dot1.y*dot2.y+dot1.z*dot2.z；
返回点；
}
//原点：光线原点，方向：光线方向，三角形顶点：p0，p1，p2。
bool-raytrialleintersect（V3f-orig，V3f-dir，V3f-p0，V3f-p1，V3f-p2）{
//计算平面法线
V3f p0p1、p0p2；
p0p1.x=p1.x-p0.x；
p0p1.y=p1.y-p0.y；
p0p1.z=p1.z-p0.z；
p0p2.x=p2.x-p0.x；
p0p2.y=p2.y-p0.y；
p0p2.z=p2.z-p0.z；
//不需要正常化
V3f N=叉积（p0p1，p0p2）；//N
//步骤1：找到P
//检查光线和平面是否平行？
float NdotRayDirection=dotProduct（N，dir）；//如果结果为0，函数将返回值false（无交点）。
if（fabs（NdotRayDirection）

---

---

谢谢你的帮助。

你得到了t=5和（0,7，-1），但是t的正确数字是1，原因是你得到了

float t=（dotProduct（N，orig）+d）/NdotRayDirection而正确的代码是
float t=-（（点积（N，orig）-d）/NdotRayDirection）这将有望解决您的问题。
您得到了t=5和（0,7，-1），但t的正确数字是1，原因是您有
浮点t=（dotProduct（N，orig）+d）/NdotRayDirection而正确的代码是
float t=-（（点积（N，orig）-d）/NdotRayDirection）这将有望解决您的问题。
您得到了t=5和（0,7，-1），但t的正确数字是1，原因是您有
浮点t=（dotProduct（N，orig）+d）/NdotRayDirection而正确的代码是
float t=-（（点积（N，orig）-d）/NdotRayDirection）这将有望解决您的问题。
您得到了t=5和（0,7，-1），但t的正确数字是1，原因是您有
浮点t=（dotProduct（N，orig）+d）/NdotRayDirection而正确的代码是
float t=-（（点积（N，orig）-d）/NdotRayDirection）这将有望解决您的问题。
如果您不想重新发明轮子，请查看。您使用了哪一个在线教程？嗨，Angew和tobi303，谢谢您的回复。Angew，我对计算机图形学很陌生，我在那个网站上找不到相关的代码。tobi303，这是链接（）。谢谢您的帮助。请尝试将
N
和
dir
正常化，看看是否有帮助。我怀疑
t
缩放错误。网站导航：浏览源代码>数学>交点>线性平面>
IntrRay3Triangle3.h
如果你不想重新发明轮子，请查看。你使用了哪一个在线教程？嗨，Angew和tobi303，谢谢你的回复。Angew，我对计算机图形学很陌生，我在那个网站上找不到相关的代码。tobi303，这是链接（）。谢谢您的帮助。请尝试将
N
和
dir
正常化，看看是否有帮助。我怀疑
t
缩放错误。网站导航：浏览源代码>数学>交点>线性平面>
IntrRay3Triangle3.h
如果你不想重新发明轮子，请查看。你使用了哪一个在线教程？嗨，Angew和tobi303，谢谢你的回复。Angew，我对计算机图形学很陌生，我在那个网站上找不到相关的代码。tobi303，这是链接（）。谢谢您的帮助。请尝试将
N
和
dir
正常化，看看是否有帮助。我怀疑
t
缩放错误
float kEpsilon = 0.000001;

V3f crossProduct(V3f point1, V3f point2){

  V3f vector; 

  vector.x = point1.y * point2.z - point2.y * point1.z; 
  vector.y = point2.x * point1.z - point1.x * point2.z; 
  vector.z = point1.x * point2.y - point1.y * point2.x; 

  return vector;
}

float dotProduct(V3f dot1, V3f dot2){

  float dot = dot1.x * dot2.x + dot1.y * dot2.y + dot1.z * dot2.z; 

  return dot;
}

//orig: ray origin, dir: ray direction, Triangle vertices: p0, p1, p2.  
bool rayTriangleIntersect(V3f orig, V3f dir, V3f p0, V3f p1, V3f p2){ 

// compute plane's normal

  V3f p0p1, p0p2;

  p0p1.x = p1.x - p0.x; 
  p0p1.y = p1.y - p0.y; 
  p0p1.z = p1.z - p0.z; 

  p0p2.x = p2.x - p0.x;
  p0p2.y = p2.y - p0.y; 
  p0p2.z = p2.z - p0.z;

  // no need to normalize
  V3f N = crossProduct(p0p1, p0p2); // N 

  // Step 1: finding P

  // check if ray and plane are parallel ?
  float NdotRayDirection = dotProduct(N, dir); // if the result is 0, the function will return the value false (no intersection).

  if (fabs(NdotRayDirection) < kEpsilon){ // almost 0 

      return false; // they are parallel so they don't intersect ! 
  }

  // compute d parameter using equation 2
  float d = dotProduct(N, p0); 

  // compute t (equation P=O+tR P intersection point ray origin O and its direction R)

  float t = (dotProduct(N, orig) + d) / NdotRayDirection; 

  // check if the triangle is in behind the ray
  //if (t < 0){ return false; } // the triangle is behind 

  // compute the intersection point using equation
  V3f P; 

  //this part should do the work, but it does not work.
  P.x = orig.x + t * dir.x; 
  P.y = orig.y + t * dir.y; 
  P.z = orig.z + t * dir.z; 


  // Step 2: inside-outside test
  V3f C; // vector perpendicular to triangle's plane 

  // edge 0
  V3f edge0; 

  edge0.x = p1.x - p0.x;
  edge0.y = p1.y - p0.y;
  edge0.z = p1.z - p0.z;

  V3f vp0; 

  vp0.x = P.x - p0.x;
  vp0.y = P.y - p0.y; 
  vp0.z = P.z - p0.z; 

  C = crossProduct(edge0, vp0); 

  if (dotProduct(N, C) < 0) { return false; }// P is on the right side 

  // edge 1
  V3f edge1;

  edge1.x = p2.x - p1.x;
  edge1.y = p2.y - p1.y;
  edge1.z = p2.z - p1.z;

  V3f vp1; 

  vp1.x = P.x - p1.x; 
  vp1.y = P.y - p1.y; 
  vp1.z = P.z - p1.z; 

  C = crossProduct(edge1, vp1); 

  if (dotProduct(N, C) < 0) { return false; } // P is on the right side 

  // edge 2
  V3f edge2;

  edge2.x = p0.x - p2.x;    
  edge2.y = p0.y - p2.y;
  edge2.z = p0.z - p2.z;

  V3f vp2; 

  vp2.x = P.x - p2.x;
  vp2.y = P.y - p2.y;
  vp2.z = P.z - p2.z;

  C = crossProduct(edge2, vp2);

  if (dotProduct(N, C) < 0) { return false; } // P is on the right side; 

  return true; // this ray hits the triangle 
}