PolygonSoup碰撞檢測學習筆記

解決PolygonSoup碰撞檢測的成熟方案是使用層次BV樹,根據BV的不同衍生出:AABB,OBB,K-dop三種方案.
在內存使用上: AABB -> OBB -> K-dop遞增.在效率上則尚無定論,根據具體的應用環境不同而不同.

Opcode: 開源,完全免費,採用AABB樹構造,內存耗用很小,支持複雜模型,起初想使用它,但發現其,基本上沒有文檔,Sample,1.3版本接口變動巨大,且無文檔,自2003年來再無更新,最致命的是不支持Instance transform的Scale,AABBCollider沒有Transform支持,而且是在看代碼時發現的,吐血.主要設計目標是少的內存耗用,但速度也很快.爲ODE,Tokamak,Crystal Space所用.

RAPID: 開源,免費,採用OBB構造,內存耗用巨大,學院派設計,最倒黴的庫,因爲其他每個庫都拿它作比較,都說比它快.

SOLID: GPL開源,商業軟件,採用AABB樹構造,2005.4月出3.5版本,效率不錯,內存耗用中等,在後面4.0會支持AABB數據壓縮,有文檔,接口簡單,代碼清晰.支持複雜模型,支持Instance transform的Scale,甚至支持各向不等縮放,支持Penetration Depth計算,cool.作者Gino van den Bergen是3D碰撞研究專家,著有:Collision Detection in Interactive 3D Environments一書,有e版.邏輯順暢地集成入DM,效果不錯,滿足使用要求.不過其IndexBuffer是unsigned int,Matrix是DX風格的,但代碼很清晰,很容易定製.
SOLID有個嚴重的Bug,必須確保DT_GetCommonPoint,DT_GetPenDepth,DT_GetClosestPair在檢測Complex和Convex對象時,第一個參數必須是Complex對象,否則崩潰.因爲SOLID的double dispatch實現有問題,不完整,沒有Convex + Complex順序的對應函數.

QuickCD: 開源,免費,K-dop樹構造.

ColDet,支持複雜模型,也不支持Instance Transform的Scale.

Reference:
http://www.codercorner.com/Opcode.htm
http://www.dtecta.com/
http://www.cs.unc.edu/~geom/OBB/OBBT.html
http://www.ams.sunysb.edu/~jklosow/quickcd/QuickCD.html
http://photoneffect.com/coldet/
http://www.merl.com/projects/vclip/
http://www.q12.org/ode/ode.html

void DmModel::SetBoundPolygonSoup(uint32 nPolygonSoupVertices, vec3_type* pVertices, uint32 nTrianglePrims, uint16* pTrianglePrims)
{    
    ASSERT(m_hSolidShape == NULL);
    
    if(nTrianglePrims == 0) return;
    if(nPolygonSoupVertices == 0) return;
            
    // Save polygon soup data for collision test!
    m_pPolygonSoupVertices = new vec3_type[nPolygonSoupVertices];
    memcpy(m_pPolygonSoupVertices, pVertices, sizeof(vec3_type) * nPolygonSoupVertices);
    
    // Build 
    DT_VertexBaseHandle hVertexBase = DT_NewVertexBase(m_pPolygonSoupVertices, 0);
    m_hSolidShape = DT_NewComplexShape(hVertexBase);
    m_hSolidVertexBase = hVertexBase;
        
    // Define mesh triangles
    for(uint32 i = 0; i < nTrianglePrims; i += 3)
    {    
        DT_VertexIndices(3, pTrianglePrims + i);
    }
    
    // End it
    DT_EndComplexShape();
}

void DmModelInstanceNode::BuildPhysics()
{    
    // Setup Phisic stuffs
    DT_ShapeHandle hShape = (DT_ShapeHandle)m_pRefModelInstance->m_pRefModel->m_hSolidShape;
    if(hShape)
    {    
        DT_ObjectHandle hObject = DT_CreateObject(NULL, hShape);
        DM_SetMatrixf(hObject, m_matTransform);
        m_hSolidObject = hObject;
    }
}

void DmGameCreature::BuildPhysics()
{    
    ASSERT(m_hSolidObject == NULL);
    
    vec3_type vMin, vMax;
    getBoundBox(vMin, vMax);
    DT_ShapeHandle hShape = DT_NewBox(vMax.x - vMin.x, vMax.y - vMin.y, vMax.z - vMin.z);
    DT_ObjectHandle hObject = DT_CreateObject(NULL, hShape);
    
    m_hSolidShape = hShape;
    m_hSolidObject = hObject;
}

bool DmModelInstanceNode::TestCollision(HANDLE hSolidObject)
{    
    if(!m_hSolidObject) return false;
    
    DT_ObjectHandle hTestObject = (DT_ObjectHandle)hSolidObject;
    DT_ObjectHandle hThisObject = (DT_ObjectHandle)m_hSolidObject;
        
    DT_Vector3 v;
    return (DT_GetCommonPoint(hThisObject, hTestObject, v) != 0);
}