射线与球的相交
字数1183 阅读142 评论0 喜欢1
今天来说说射线和球的相交检测。
从图形来说
![射线和圆相交, origin 是射线起点， dir 是射线的方向向量。p0，p1是两个交点，center 为圆心，半径为R ，d 为圆心到射线的距离][ray-sphere] 我们先以2D 切面图来说明，当射线和圆相交的时候，可以看到，球心 center 到射线 ray 的距离 d
1. 设圆心在射线上的投影为c' ，则 origin ，center, c' 形成了一个直角三角形。
2. 获得射线起点到圆心的向量 Voc = Vcenter - Vorigin
3. 在射线方向上的投影为: Poc= Voc·(Voc·dir)
4. 勾股定理：d·d = Voc·Voc - Poc·Poc
可以求出d 的数值，
∙
∙
∙ d R，射线在圆外，没有交点。
接下来求P0，P1：
1. c' ，center ，P0 or P1点构成直角三角形。
2. P0 or P1到c' 的距离 tca·tca = R·R - d·d;
3. 有如下式子 P0 = dir·( |Poc| - tca );
P1 = dir·( |Poc| + tca );
要注意，没有交点的时候， tca·tca
推导三维情况可以照上面的去做，dot 能保证投影点在同一个平面上的。 附代码
bool Intersect (const Ray& ray, const Sphere& sphere, float & t0, float & t1) {
Vector3 oc = sphere.GetCenter() - ray.GetOrigin();
float projoc = dot(ray.GetDirection(), oc);
if (projoc
return false;
float oc2 = dot(oc, oc);
float distance2 = oc2 - projoc * projoc;
if (distance2 > sphere.GetRadiusSquare())
return false;
float discriminant = sphere.GetRadiusSquare() - distance2;
if (discriminant
t0 = t1 = projoc;
else
{
discriminant = sqrt (discriminant);
t0 = projoc - discriminant;
t1 = projoc + discriminant;
if (t0
t0 = t1;
}
return true;
}
从方程角度来看
射线方程：ray : P(t) = O + D·t ( t >= 0 )
球的方程：sphere : sqr( P-C ) = R·R (sqr(x) = x^2 = x·x)
O=origin, D=direction, C=center, R=radius
射线方程表明的是如下一个点的集合P ，当t 从零增大时, D·t 会沿着D 向量的方向从零逐步变长，t 取值无限表示了射线单方向。从O 点开始在D 方向上无限个点构成了一条射线。
球的方程表明了任何点P ，只要到C 点的距离等于半径R ，则表明点在球面上，这么一个球面上的点的集合。
因此当射线与球相交的时候，这个点既在射线上，又在球面上。等式射线的P(t) = 球的P 成立。
联立两个方程，试着求解 t 有：
s qr( O + D·t - C ) = R·R
设 O-C=OC，有：
sqr( OC+D·t ) - R·R = 0
//展开得到如下式子
=>D ·D ·t ·t + 2·OC·D ·t + OC·OC - R·R = 0
=> (D ·D )·t ·t + 2·(OC·D )·t + OC·OC - R·R = 0
因为 D 是单位向量有D·D = dot(D, D) = 1最后方程为：
t·t + 2·(OC·D)·t + OC·OC - R·R = 0;
这是一个关于 t 的二次方程at^2 + bt + c = 0那么解就已经出来了： ∙
∙ t0 = -(b + √Δ) / 2a t1 = -(b - √Δ) / 2a
o
o
o
o a = D·D = dot(D, D) = 1; b = 2·OC·D = 2·dot(OC, D); c = OC·OC - R·R = dot(OC, OC) - R·R; 判别式 Δ = sqr(b) - 4ac
= 4·sqr( OC·D ) - 4·( OC·OC - R·R )
= 4·( sqr( OC·D ) - OC·OC + R·R );
如果判别式 Δ > 0，则表明球与射线相交。 根据以上方程，我们其中试着展开 t 的式子
t0 = -(b + √Δ) / 2a = -(b + √Δ) / 2·1
= -b/2 - √(Δ/4)
= -dot(OC, D) - √( sqr( dot(OC, D) ) - dot(OC, OC) + R·R ) 求出 t 后可以根据 P(t) = O + D * t 得到交点。
附代码
bool Intersect (const Ray& ray, const Sphere& sphere, float & t0, float & t1) {
Vector3 oc = ray.GetOrigin() - sphere.GetCenter();
float dotOCD = dot(ray.GetDirection(), oc);
if (dotOCD >0)
return false;
float dotOC = dot(oc, oc);
float discriminant = dotOCD*dotOCD - dotOC + sphere.GetRadiusSquare();
if (discriminant
return false;
elseif (discriminant
t0 = t1 = -dotOCD;
else
{
discriminant = sqrt (discriminant);
t0 = -dotOCD - discriminant;
t1 = -dotOCD + discriminant;
if (t0
t0 = t1;
}
return true;
}
补充一些内容：
交点的法线
因为交点在球面上，球面法线反向是从球心指向球面的点。 设交点为 IntersecionP ，只需要简单的计算：
normal = ( IntersectionP - sphere.GetCenter() ).Normalize ();
射线与球的相交
字数1183 阅读142 评论0 喜欢1
今天来说说射线和球的相交检测。
从图形来说
![射线和圆相交, origin 是射线起点， dir 是射线的方向向量。p0，p1是两个交点，center 为圆心，半径为R ，d 为圆心到射线的距离][ray-sphere] 我们先以2D 切面图来说明，当射线和圆相交的时候，可以看到，球心 center 到射线 ray 的距离 d
1. 设圆心在射线上的投影为c' ，则 origin ，center, c' 形成了一个直角三角形。
2. 获得射线起点到圆心的向量 Voc = Vcenter - Vorigin
3. 在射线方向上的投影为: Poc= Voc·(Voc·dir)
4. 勾股定理：d·d = Voc·Voc - Poc·Poc
可以求出d 的数值，
∙
∙
∙ d R，射线在圆外，没有交点。
接下来求P0，P1：
1. c' ，center ，P0 or P1点构成直角三角形。
2. P0 or P1到c' 的距离 tca·tca = R·R - d·d;
3. 有如下式子 P0 = dir·( |Poc| - tca );
P1 = dir·( |Poc| + tca );
要注意，没有交点的时候， tca·tca
推导三维情况可以照上面的去做，dot 能保证投影点在同一个平面上的。 附代码
bool Intersect (const Ray& ray, const Sphere& sphere, float & t0, float & t1) {
Vector3 oc = sphere.GetCenter() - ray.GetOrigin();
float projoc = dot(ray.GetDirection(), oc);
if (projoc
return false;
float oc2 = dot(oc, oc);
float distance2 = oc2 - projoc * projoc;
if (distance2 > sphere.GetRadiusSquare())
return false;
float discriminant = sphere.GetRadiusSquare() - distance2;
if (discriminant
t0 = t1 = projoc;
else
{
discriminant = sqrt (discriminant);
t0 = projoc - discriminant;
t1 = projoc + discriminant;
if (t0
t0 = t1;
}
return true;
}
从方程角度来看
射线方程：ray : P(t) = O + D·t ( t >= 0 )
球的方程：sphere : sqr( P-C ) = R·R (sqr(x) = x^2 = x·x)
O=origin, D=direction, C=center, R=radius
射线方程表明的是如下一个点的集合P ，当t 从零增大时, D·t 会沿着D 向量的方向从零逐步变长，t 取值无限表示了射线单方向。从O 点开始在D 方向上无限个点构成了一条射线。
球的方程表明了任何点P ，只要到C 点的距离等于半径R ，则表明点在球面上，这么一个球面上的点的集合。
因此当射线与球相交的时候，这个点既在射线上，又在球面上。等式射线的P(t) = 球的P 成立。
联立两个方程，试着求解 t 有：
s qr( O + D·t - C ) = R·R
设 O-C=OC，有：
sqr( OC+D·t ) - R·R = 0
//展开得到如下式子
=>D ·D ·t ·t + 2·OC·D ·t + OC·OC - R·R = 0
=> (D ·D )·t ·t + 2·(OC·D )·t + OC·OC - R·R = 0
因为 D 是单位向量有D·D = dot(D, D) = 1最后方程为：
t·t + 2·(OC·D)·t + OC·OC - R·R = 0;
这是一个关于 t 的二次方程at^2 + bt + c = 0那么解就已经出来了： ∙
∙ t0 = -(b + √Δ) / 2a t1 = -(b - √Δ) / 2a
o
o
o
o a = D·D = dot(D, D) = 1; b = 2·OC·D = 2·dot(OC, D); c = OC·OC - R·R = dot(OC, OC) - R·R; 判别式 Δ = sqr(b) - 4ac
= 4·sqr( OC·D ) - 4·( OC·OC - R·R )
= 4·( sqr( OC·D ) - OC·OC + R·R );
如果判别式 Δ > 0，则表明球与射线相交。 根据以上方程，我们其中试着展开 t 的式子
t0 = -(b + √Δ) / 2a = -(b + √Δ) / 2·1
= -b/2 - √(Δ/4)
= -dot(OC, D) - √( sqr( dot(OC, D) ) - dot(OC, OC) + R·R ) 求出 t 后可以根据 P(t) = O + D * t 得到交点。
附代码
bool Intersect (const Ray& ray, const Sphere& sphere, float & t0, float & t1) {
Vector3 oc = ray.GetOrigin() - sphere.GetCenter();
float dotOCD = dot(ray.GetDirection(), oc);
if (dotOCD >0)
return false;
float dotOC = dot(oc, oc);
float discriminant = dotOCD*dotOCD - dotOC + sphere.GetRadiusSquare();
if (discriminant
return false;
elseif (discriminant
t0 = t1 = -dotOCD;
else
{
discriminant = sqrt (discriminant);
t0 = -dotOCD - discriminant;
t1 = -dotOCD + discriminant;
if (t0
t0 = t1;
}
return true;
}
补充一些内容：
交点的法线
因为交点在球面上，球面法线反向是从球心指向球面的点。 设交点为 IntersecionP ，只需要简单的计算：
normal = ( IntersectionP - sphere.GetCenter() ).Normalize ();