Direct3D中实现图元的鼠标拾取

在观察坐标系中，Dir是一条从观察坐标原点出发的射线，所以我们只需要再确定一个该射线经过的点，就可以得到它在观察坐标系中的表示。假设我们要求的射线上的另外一点为该射线与透视投影平截头体近剪切面的交点，针对最普遍的透视投影而言，透视投影平截头体经投影变换后，变成一个1/2立方体（请允许我这么叫^_^，因为它的大小为一个正方体的一半，x,y方向边长为2，z方向为1）如图：

投影坐标系以近剪切面中心为坐标原点，该立方体从z轴负向看过去与图形程序视区相对应，最终近剪切面（前剪切面）上一点与屏幕坐标之间的对应关系如下图所示：

根据比例关系，screenPt与投影空间上的点projPt之间的关系为

假设图形程序窗口的宽为screenWidth,高为screenHeight,

projPt.x = (screenPt.x-screenWidth/2)/screenWidth*2; （公式1）

projPt.y = (screenPt.y-screenHeight/2)/screenHeight*2; （公式2）

projPt.z =0;（实际该值可任意取，不影响最终结果。为了处理简单，我们取改值为0，表示该点取在近剪切面上）

得到projPt后，我们需要做的是把该点坐标从投影空间转换到观察空间(view space),

根据透视投影的定义，可假设点(projPt.x，projPt.y，projPt.z)

对应的其次坐标为

(projPt.x*projPt.w，projPt.y*projPt.w，projPt.z*projPt.w，projPt.w)

我们可以通过GetTransform(D3DTS_PROJECTION,&ProjMatrix)函数获得投影矩阵ProjMatrix,则根据观察空间到投影空间的变换关系则

(projPt.x*projPt.w，projPt.y*projPt.w，projPt.z*projPt.w，projPt.w)

= (viewPt.x，viewPt.y，viewPt.z, 1)*pProjMatrx;

根据定义和图形学原理

ProjMatrix = <shape id="_x0000_s1029" style="WIDTH: 101.25pt; HEIGHT: 81pt" coordsize="21600,21600" type="#_x0000_t75"><imagedata o:title="" src="pick_2004_529.files/image007.png"></imagedata></shape>=<shape id="_x0000_s1028" style="WIDTH: 103.5pt; HEIGHT: 85.5pt" coordsize="21600,21600" type="#_x0000_t75"><imagedata o:title="" src="pick_2004_529.files/image009.png"></imagedata></shape>

所以,

(projPt.x*projPt.w，projPt.y*projPt.w，projPt.z*projPt.w，projPt.w)

= ( viewPt.x*ProjMatrix._m11,

viewPt.y*ProjMatrix._m22,

viewPt.z*Q-QZn,

viewPt.z)

所以

projPt.x*projPt.w = viewPt.x*ProjMatrix._m11

projPt.y*projPt.w = viewPt.y*ProjMatrix._m22

projPt.z*projPt.w = viewPt.z*Q-QZn （注意projPt.z = 0）

projPt.w = viewPt.z;

解得

viewPt.x = projPt.x*Zn/ ProjMatrix._m11;

viewPt.y = projPt.y*Zn/ ProjMatrix._m22;

viewPt.z = Zn;

好了，到这里为止我们终于求出了射线与近剪切面交点在观察坐标系中的坐标，现在我们拥有了射线的出发点(0,0,0)和射线方向上另外一点(viewPt.x,viewPt.y,viewPt.z),则该射线的方向矢量在观察空间中的表示可确定为（viewPt.x-0,viewPt.y-0,viewPt.z-0）,化简一下三个分量同除近剪切面z坐标Zn，该方向矢量可写作

DIRview = (projPt.x/projMatrix._m11,projPt.y/projMatrix._m22,1)

代入公式1，公式2

DIRview.x = (2*screenPt.x/screenWidth-1)/projMatrix._m11;

DIRview.y = (2*screenPt.y/screenHeight-1)/projMatrix._m22;

DIRview.z = 1;

其中screenWidth和screenHeight可以通过图像显示的backBuffer的目标表面（D3DSURFACE_DESC）来获得，该表面在程序初始化时由用户创建。

2.1.3 转换Dir到世界坐标空间，并得到观察点在世界坐标系中的坐标

由于最终的运算要在世界坐标空间中进行，所以我们还需要把矢量DIRview从观察空间转换为世界坐标空间中的矢量DIRworld。

因为

DIRview = DIRworld*ViewMatrix;

其中ViewMatrix为观察矩阵，在D3D中可以用函数GetTransform( D3DTS_VIEW, &ViewMatrix )得到。

所以DIRworld = DIRview * inverse_ViewMatrix,其中inverse_ViewMatrix为

ViewMatrix的逆矩阵。

观察点在观察坐标系中坐标为OriginView（0，0，0，1），所以其在世界坐标系中的坐标同样可以利用ViewMatrix矩阵，反变换至世界坐标系中，事实上我们可以很简单的判断出,其在世界坐标系中的表示为:

OriginWorld = (inverse_ViewMatrix._41,

inverse_ViewMatrix._42,

inverse_ViewMatrix._43,

1);

到这里为止，判断射线与三角面是否相交的条件就完全具备了。

2.2 使用射线矢量对场景中的所有三角形图元求交，获得三角形索引值和重心坐标。

这一步骤地实现由两种途径:

第一种方法非常简单，利用D3D提供的扩展函数D3DXIntersect可以轻松搞定一切。见2.1

第二种方法就是我们根据空间解析几何的知识，自己来完成射线三角形的求交算法。一般来讲，应用上用第一种方法就足够了，但是我们如果要深入的话，必须理解相交检测的数学算法，这样才能自由的扩展，面对不同的需求，内容见2.2

下面分别讲解两种实现途径：

2.2.1 D3D扩展函数实现求交

这种方法很简单也很好用，对于应用来说应尽力是用这种方式来实现，毕竟效率比自己写得要高得多。

实际上其实没什么好讲的，大概讲一下函数D3DXIntersect吧

D3D SDK该函数声明如下

HRESULTD3DXIntersect(

LPD3DXBASEMESHpMesh,

CONST D3DXVECTOR3*pRayPos,

CONST D3DXVECTOR3*pRayDir,

BOOL*pHit,

DWORD*pFaceIndex,

FLOAT*pU,

FLOAT*pV,

FLOAT*pDist,

LPD3DXBUFFER*ppAllHits,

DWORD*pCountOfHits

);

l pMesh指向一个ID3DXBaseMesh的对象，最简单的方式是从.x文件获得，描述了要进行相交检测的三角面元集合的信息，具体规范参阅direct9 SDK

l pRayPos 指向射线发出点

l pRayDir 指向前面我们辛辛苦苦求出的射线方向的向量

l pHit 当检测到相交图元时，指向一个true,不与任何图元相交则为假

l pU 用于返回重心坐标U分量

l pV返回重心坐标V分量

l pDist 返回射线发出点到相交点的长度

注意：以上红色字体部分均指最近的一个返回结果（即*pDist最小）

l ppAllHits用于如果存在多个相交三角面返回相交的所有结果

l pCountOfHits 返回共有多少个三角形与该射线相交

补充：重心坐标的概念

其中pU和pV用到了重心坐标的概念，下面稍作描述

一个三角形有三个顶点，在迪卡尔坐标系中假设表示为V1(x1,y1,z1),V2(x2,y2,z2),V3(x3,y3,z3),则三角形内任意一点的坐标可以表示为 pV = V1 + U(V2-V1) + V(V3-V1),所以已知三个顶点坐标的情况下，任意一点可用坐标(U,V)来表示，其中参数U控制V2在结果中占多大的权值，参数V控制V3占多大权值，最终1－U－V控制V1占多大权值，这种坐标定义方式就叫重心坐标。

2..2.2射线三角面相交的数学算法

使用d3d扩展函数，毕竟有时不能满足具体需求，掌握了该方法，我们才能够获得最大的控制自由度，任意修改算法。

已知条件: 射线源点orginPoint,三角形三个顶点 v1,v2,v3,射线方向 Dir

（均以三维坐标向量形式表示）

算法目的: 判断射线与三角形是否相交，如果相交求出交点的重心坐标(U,V)和射线原点到交点的距离T。

我们可先假设射线与三角形相交则交点(注以下均为向量运算，*数乘，dot(X,Y) X，Y 点乘，cross（X，Y）X，Y叉乘；U，V，T为标量)

则：

IntersectPoint = V1 + U*(V2-V1) + V*(V3-V1) ;

IntersectPoint = originPoint + T*Dir；

所以

orginPoint + T*Dir = V1 + U*(V2-V1) + V*(V3-V1);

整理得：

这是一个简单的线性方程组，若有解则行列式<shape id="_x0000_s1026" style="WIDTH: 115.5pt; HEIGHT: 15.75pt" coordsize="21600,21600" type="#_x0000_t75"><imagedata o:title="" src="pick_2004_529.files/image013.png" grayscale="t"></imagedata></shape>不为0。

根据T,U,V的含义当T>0, 0<U<1,0<V<1,0<U+V<1时该交点在三角形内部，

解此方程组即可获得我们关心的值,具体解法不再赘述，克莱姆法则就够了（详细见线性代数）:射线原点到相交点的距离T,和交点的中心坐标(U,V)。

下面给出Direct 9 SDK示例程序中的实现代码

IntersectTriangle( const D3DXVECTOR3& orig,

const D3DXVECTOR3& dir, D3DXVECTOR3& v0,

D3DXVECTOR3& v1, D3DXVECTOR3& v2,

FLOAT* t, FLOAT* u, FLOAT* v )

{

// 算出两个边的向量

D3DXVECTOR3 edge1 = v1 - v0;

D3DXVECTOR3 edge2 = v2 - v0;

D3DXVECTOR3 pvec;

D3DXVec3Cross( &pvec, &dir, &edge2 );

// 如果det为0，或接近于零则射线与三角面共面或平行，不相交

//此处det就相当于上面的<shape id="_x0000_s1025" style="WIDTH: 115.5pt; HEIGHT: 15.75pt" coordsize="21600,21600" type="#_x0000_t75"><imagedata o:title="" src="pick_2004_529.files/image013.png" grayscale="t"></imagedata></shape>，

FLOAT det = D3DXVec3Dot( &edge1, &pvec );

D3DXVECTOR3 tvec;

if( det > 0 )

{

tvec = orig - v0;

}

else

{

tvec = v0 - orig;

det = -det;

}

if( det < 0.0001f )

return FALSE;

// 计算u并测试是否合法（在三角形内）

*u = D3DXVec3Dot( &tvec, &pvec );

if( *u < 0.0f || *u > det )

return FALSE;

// Prepare to test V parameter

D3DXVECTOR3 qvec;

D3DXVec3Cross( &qvec, &tvec, &edge1 );

//计算u并测试是否合法（在三角形内）

*v = D3DXVec3Dot( &dir, &qvec );

if( *v < 0.0f || *u + *v > det )

return FALSE;

/*计算t,并把t,u,v放缩为合法值（注意前面的t,v,u不同于算法描述中的相应量，乘了一个系数det）,注意：由于该步运算需要使用除法，所以放到最后来进行，避免不必要的运算，提高算法效率*/

*t = D3DXVec3Dot( &edge2, &qvec );

FLOAT fInvDet = 1.0f / det;

*t *= fInvDet;

*u *= fInvDet;

*v *= fInvDet;

return TRUE;

}

2.2.3 拾取完成根据获得的中心坐标计算我们关心的常见量，。

根据重心坐标（U,V）,我们可以很容易的算出各种相关量比如纹理坐标和交点的差值颜色，假设以纹理坐标为例设V1,V2,V3的纹理坐标分别为T1(tu1,tv1),T2(tu2,tv2),T3(tu3,tv3)则交点的坐标为

IntersectPointTexture = T1 + U(T2-T1) + V(T3-T1)

3、结束及声明

Ok, 到这里为止关于拾取的相关知识就介绍完了，小弟第一次写这种文章，不知道有没有把问题说清楚，希望对大家有所帮助，有任何问题可以给我发email: jzhang1@mail.xidian.edu.cn

或者到我的网站留言： www.heavysword.com

声明：

本文写作的目的是为了广大D3D学习者方便学习服务，文中算法为作者参考相关文献总结，作者无意把这些据为自己的成果，所有权原算法提出者所有（参阅参考文献），文中代码为D3d SDK的示例内容，由笔者进行了必要的解释，代码版权归microsoft所有。

4、参考文献

【1】Microsoft DirectX 9.0 SDK,microsoft

【2】fast,Minimun Storage Ray/Triangle Intersection,Tomas Moler,Ben Trumbore

(转载请注明出处）

分享到：

VC中Unicode字符串的处理 | 二维离散付立叶变换及其性质

2007-04-23 23:30
浏览 549
评论(0)
查看更多

发表评论

您还没有登录,请您登录后再发表评论

最近访客更多访客>>

博主相关

文章分类

社区版块

存档分类

最新评论