[每日一码] 求Line和Su**ce的交点

admin

下面是示例代码，求端点是(0 0 0) (0 0 1)的LINE和Su**ce的交点。

注意：include和lib文件在  ObjectARX SDK\\utils\brep 目录下。

1. static AcGePoint3dArray Intersect(AcDbSu**ce* pSu**ce,AcGeLine3d line)
2. {
3.     AcGePoint3dArray returnPtArray;
4.   AcDbBody* pBody = new AcDbBody();
6.   // 2013
7.   // Acad::ErrorStatus es = pBody->setASMBody(pSu**ce->ASMBodyCopy());
8.   //before 20123
9.    Acad::ErrorStatus es = pBody->setBody(pSu**ce->body());
11.    //build AcBrBrep
12.    AcBrBrep* pBrep = new AcBrBrep();
13.    //
14.   if(AcBr::eOk == pBrep->set(*pBody))
15.   {
16.     AcBrBrepFaceTraverser* pFaceTrav = new AcBrBrepFaceTraverser;
17.     if(AcBr::eOk == pFaceTrav->setBrep(*pBrep))
18.    {
19.        for(pFaceTrav->restart();!pFaceTrav->done();pFaceTrav->next())
20.        {
21.         AcB**ce face;
23.         if(AcBr::eOk == pFaceTrav->getFace(face))
24.         {
25.           double area = 0.0f;
26.          face.getSu**ceArea(area);
28.          acutPrintf(L"\nSu**ce Area: %f", area);
30.         //*****whole su**ce of the Brep face******        
31.         //AcGeNurbSu**ce nurbSu**ce;
32.         //face.getSu**ceAsNurb(nurbSu**ce);
33.         //AcGeCurveSurfInt curveSI;
34.         ////input the curve and line
35.         //curveSI.set(line,nurbSu**ce);
36.         ////get the count of intersect points
37.         //int count = curveSI.numIntPoints(err采用1);
38.         //    if(err采用1 == AcGe::kXXOk && count >0 )
39.         //    {      
40.         //        AcGeIntersectError err采用2;
41.         //        for(int index = 0 ;index < count; index ++)
42.         //        {
43.         //           AcGePoint3d pt =
44.         //             curveSI.intPoint(index,err采用2);      
45.         //            returnPtArray.append(pt);
46.         //        }      
47.         //      
48.         //    }
49.         //**********
51.         //****real su**ce of the orignal AcDbSu**ce      
52.           AcGeExternalBoundedSu**ce** nurbs = NULL;
53.           Adesk::UInt32 numNurbs = 0;
54.           face.getSu**ceAsTrimmedNurbs(numNurbs,nurbs);
55.           //*****
57.         for (Adesk::UInt32 i = 0; i < numNurbs; i++)
58.         {
59.             AcGeCurveSurfInt curveSI;
60.             AcGeIntersectError  err采用1 = AcGe::kXXOk;         
61.              //input the curve and line
62.             curveSI.set(line,*nurbs<i>);
64.             //get the count of intersect points
65.             int count = curveSI.numIntPoints(err采用1);
66.             if(err采用1 == AcGe::kXXOk && count >0 )
67.             {
68.                 AcGeIntersectError err采用2;
69.                 for(int index = 0 ;index < count; index ++)
70.                 {
71.                     AcGePoint3d pt =
72.                         curveSI.intPoint(index,err采用2);      
73.                     returnPtArray.append(pt);
74.                 }      
75.             }
76.             delete nurbs<i>;
77.         }
78.         // your responsibility to delete the
79.        // array of su**ces
80.         delete[] nurbs;
81.       }
82.     }
83.   }
84.   delete pFaceTrav;
85.   }
86.   delete pBrep;
89.   return returnPtArray;
91. }

复制代码