如何计算出一条直线与一个元素(体、曲面等)的所有交点？[MicroStation V8i C++]

Question

函数mdlIntersect_allBetweenElms能够求得某些情况下两个元素之间的交点，但该函数局限性较大，例如只能得到从某个角度看过去的表象交点。同时该函数也要求操作元素必须是已经add到dgn中，而仅存在内存中的元素则无法操作。 
 现在将条件放宽，将其中一个元素确定为直线(且垂直于Z=0平面)，那么是否有比较稳定的办法可以求出该直线与某元素之间的交点呢？

Yongan.Fu · Accepted Answer

经测试，这个mdlKISolid_rayTestBody确实有问题。根据其帮助文档找到了一个替代函数mdlKISolid_getEntities能正确返回交点坐标。测试代码如下： 
 void lineSrufaceIntersect()
 {
 MSElementDescrP edP = NULL;
 mdlAssoc_getElementDescr(&edP, NULL, 621, ACTIVEMODEL, FALSE);

DPoint3d intersectPt, rayOrg = { 1929, 11371, 25396 }, rayDir = { 0, 0, -1 };
 KIBODY* pBody = NULL;
 double* params = NULL;
 KIENTITY_LIST* hits = NULL;
 Transform trans;
 int status = mdlKISolid_elementToBody2(&pBody, &trans, edP, MASTERFILE, 1, FALSE);
 Transform invBodyTransform;
 mdlTMatrix_getInverse(&invBodyTransform, &trans);
 mdlKISolid_beginCurrTrans(MASTERFILE);
 mdlCurrTrans_invtransPointArray(&rayOrg, &rayOrg, 1);
 mdlTMatrix_transformPoint(&rayOrg, &invBodyTransform);
 mdlKISolid_endCurrTrans();
 mdlTMatrix_transformVector((DVec3d*)&rayDir, &invBodyTransform);
 
 mdlKISolid_listCreate(&hits);
 if (SUCCESS == mdlKISolid_getEntities(hits, &params, &rayOrg, &rayDir, 1.0E-8,
 mdlKISolid_faceType(), pBody, ACTIVEMODEL))
 {
 int hitCnt;
 mdlKISolid_listCount(&hitCnt, hits);
 for (int i = 0; i<hitCnt; i++)
 {
 mdlVec_projectPoint(&intersectPt, &rayOrg, (DVec3d*)&rayDir, params[i]);
 mdlKISolid_beginCurrTrans(MASTERFILE);
 mdlTMatrix_transformPoint(&intersectPt, &trans);
 mdlCurrTrans_transformPointArray(&intersectPt, &intersectPt, 1);
 mdlKISolid_endCurrTrans();
 displayMsg("x=%.2f, y=%.2f, z=%.2f", intersectPt.x, intersectPt.y, intersectPt.z);
 }
 dlmSystem_mdlFree(params);
 mdlKISolid_listDelete(&hits);
 }
 else mdlDialog_dmsgsPrint("ERROR in mdlKISolid_getEntities");
 mdlKISolid_freeBody(pBody);
 mdlElmdscr_freeAll(&edP);
 }