我正在执行ray和矩形相交测试。为此,我首先测试射线是否与平面相交,然后查看射线是否位于矩形范围内。
以下是代码:
float intersectQuad(Ray r, float3 p1, float3 p2, float3 p3, float3 p4, float3* normal)
{
float3 x1 = p2 - p1;
float3 x2 = p4 - p1;
float t;
float3 n = normalize(cross(x2, x1));
float denom = dot(n, r.dir);
if (denom > 0.00000000001)
{
float3 p0l0 = normalize(p1 - r.origin);
t = dot(p0l0, n) / denom;
printf(" %f ", t);
if( t > 0.000000000001f )
{
float3 hitPoint = r.origin + r.dir * t;
float3 V1 = normalize(p2 - p1);
float3 V2 = normalize(p3 - p2);
float3 V3 = normalize(p4 - p3);
float3 V4 = normalize(p1 - p4);
float3 V5 = normalize(hitPoint - p1);
float3 V6 = normalize(hitPoint - p2);
float3 V7 = normalize(hitPoint - p3);
float3 V8 = normalize(hitPoint - p4);
if (dot(V1,V5) < 0.0f) return 0.0f;
if (dot(V2,V6) < 0.0f) return 0.0f;
if (dot(V3,V7) < 0.0f) return 0.0f;
if (dot(V4,V8) < 0.0f) return 0.0f;
*normal = n;
return t;
}
}
return 0.0f;
}
我的理解是t
的值应为0.0f < t < 1.0f
但是当我打印值t
时,有时会看到值> 1.0f
。
我的代码有什么问题吗?
答案 0 :(得分:1)
首先,如果将射线方向矢量针对r.origin + r.direction*t
进行了归一化,则t
是从原点到hitPoint
的距离。该距离可以有任何值。
第二,我将返回值类型更改为bool
,目的是避免浮点运算,以便在函数的调用范围内再次将结果与0.0f进行比较。.
第三,更多地出于个人喜好,请使用早期退出技术以下列方式重新格式化代码,并在可能的地方标记const
值。
bool intersectQuad(Ray r, float3 p1, float3 p2, float3 p3, float3 p4, float3* outNormal, float* outT)
{
const float3 x1 = p2 - p1;
const float3 x2 = p4 - p1;
const float3 n = normalize(cross(x2, x1));
const float denom = dot(n, r.dir);
if (denom < 0.00000000001) return false;
const float3 p0l0 = normalize(p1 - r.origin);
const float t = dot(p0l0, n) / denom;
printf(" %f ", t);
if( t < 0.000000000001f ) return false;
const float3 hitPoint = r.origin + r.dir * t;
const float3 V1 = normalize(p2 - p1);
const float3 V2 = normalize(p3 - p2);
const float3 V3 = normalize(p4 - p3);
const float3 V4 = normalize(p1 - p4);
const float3 V5 = normalize(hitPoint - p1);
const float3 V6 = normalize(hitPoint - p2);
const float3 V7 = normalize(hitPoint - p3);
const float3 V8 = normalize(hitPoint - p4);
if (dot(V1,V5) < 0.0f) return false;
if (dot(V2,V6) < 0.0f) return false;
if (dot(V3,V7) < 0.0f) return false;
if (dot(V4,V8) < 0.0f) return false;
*outNormal = n;
*outT = t;
return true;
}