求Line和Su**ce的交点

LoveArx

下面是示例代码，求端点是(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();
   
5. 
   
6.   // 2013
   
7.   // Acad::ErrorStatus es = pBody->setASMBody(pSu**ce->ASMBodyCopy());
   
8.   //before 20123
   
9.    Acad::ErrorStatus es = pBody->setBody(pSu**ce->body());
   
10. 
    
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;
    
22. 
    
23.         if(AcBr::eOk == pFaceTrav->getFace(face))
    
24.         {
    
25.           double area = 0.0f;
    
26.          face.getSu**ceArea(area);
    
27. 
    
28.          acutPrintf(L"\nSu**ce Area: %f", area);
    
29. 
    
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.         //**********
    
50. 
    
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.           //*****
    
56. 
    
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);
    
63. 
    
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;
    
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;
    
87. 
    
88. 
    
89.   return returnPtArray;
    
90. 
    
91. }
    

static AcGePoint3dArray Intersect(AcDbSu**ce* pSu**ce,AcGeLine3d line)<br /> { <br />     AcGePoint3dArray returnPtArray;<br />   AcDbBody* pBody = new AcDbBody();<br /> <br />   // 2013<br />   // Acad::ErrorStatus es = pBody->setASMBody(pSu**ce->ASMBodyCopy());<br />   //before 20123<br />    Acad::ErrorStatus es = pBody->setBody(pSu**ce->body());<br /> <br />    //build AcBrBrep<br />    AcBrBrep* pBrep = new AcBrBrep();<br />    //<br />   if(AcBr::eOk == pBrep->set(*pBody))<br />   {<br />     AcBrBrepFaceTraverser* pFaceTrav = new AcBrBrepFaceTraverser;<br />     if(AcBr::eOk == pFaceTrav->setBrep(*pBrep))<br />    {<br />        for(pFaceTrav->restart();!pFaceTrav->done();pFaceTrav->next())<br />        {<br />         AcB**ce face;<br /> <br />         if(AcBr::eOk == pFaceTrav->getFace(face))<br />         {<br />           double area = 0.0f;<br />          face.getSu**ceArea(area);<br /> <br />          acutPrintf(L"\nSu**ce Area: %f", area);<br /> <br />         //*****whole su**ce of the Brep face******        <br />         //AcGeNurbSu**ce nurbSu**ce;<br />         //face.getSu**ceAsNurb(nurbSu**ce);<br />         //AcGeCurveSurfInt curveSI;<br />         ////input the curve and line<br />         //curveSI.set(line,nurbSu**ce);<br />         ////get the count of intersect points<br />         //int count = curveSI.numIntPoints(err_1);<br />         //    if(err_1 == AcGe::kXXOk && count >0 )<br />         //    {       <br />         //        AcGeIntersectError err_2;<br />         //        for(int index = 0 ;index < count; index ++)<br />         //        {<br />         //           AcGePoint3d pt =<br />         //             curveSI.intPoint(index,err_2);      <br />         //            returnPtArray.append(pt);<br />         //        }      <br />         //      <br />         //    }<br />         //**********<br /> <br />         //****real su**ce of the orignal AcDbSu**ce       <br />           AcGeExternalBoundedSu**ce** nurbs = NULL;<br />           Adesk::UInt32 numNurbs = 0;<br />           face.getSu**ceAsTrimmedNurbs(numNurbs,nurbs);<br />           //*****<br /> <br />         for (Adesk::UInt32 i = 0; i < numNurbs; i++)<br />         {<br />             AcGeCurveSurfInt curveSI;<br />             AcGeIntersectError  err_1 = AcGe::kXXOk;          <br />              //input the curve and line<br />             curveSI.set(line,*nurbs<i>);<br /> <br />             //get the count of intersect points<br />             int count = curveSI.numIntPoints(err_1);<br />             if(err_1 == AcGe::kXXOk && count >0 )<br />             {<br />                 AcGeIntersectError err_2;<br />                 for(int index = 0 ;index < count; index ++)<br />                 {<br />                     AcGePoint3d pt =<br />                         curveSI.intPoint(index,err_2);      <br />                     returnPtArray.append(pt);<br />                 }      <br />             }<br />             delete nurbs<i>;<br />         }<br />         // your responsibility to delete the<br />        // array of su**ces<br />         delete[] nurbs;<br />       }<br />     }<br />   }<br />   delete pFaceTrav;<br />   }<br />   delete pBrep;<br /> <br /> <br />   return returnPtArray;<br /> <br /> }<br />

游客，如果您要查看本帖隐藏内容请回复

sdlsxygh

学习学习学习