[原创]：测试点和闭合曲线的位置关系

wkai

学习ARX一周，今天终于写出了第一段代码，请大家批评指正.

射线法测试点和闭合曲线的位置关系。
考虑了 射线与曲线相切的特殊情况。

1. 
   
2. Adesk::Boolean PointRelationToCurve(ads_point pt,AcDbCurve *pPline)
   
3. {
   
4.                 //acutPrintf("x %d",pt[X]);
   
5.                 ads_point tmpP;
   
6.                 AcDbRay *pRay =new AcDbRay();
   
7.                 pRay->setBasePoint(asPnt3d(pt));
   
8.                 pRay->setUnitDir(AcGeVector3d(1.0,0.0,0.0));
   
9. 
   
10. 
    
11.                 AcGePoint3d pointOnCurve;       
    
12.                 pPline->getClosestPointTo(asPnt3d(pt),pointOnCurve);
    
13. 
    
14.                 ads_point pt1 = {pointOnCurve.x,pointOnCurve.y,pointOnCurve.z};
    
15.                 if(acutDistance(pt,pt1)==0)
    
16.                 {
    
17.                         acutPrintf("\nON");
    
18.                         return Adesk::kTrue;
    
19.                 }
    
20.                 AcGePoint3dArray points,points1;
    
21.                 pPline->intersectWith(pRay,AcDb::kOnBothOperands,points);
    
22.                 acutPrintf("\n交点数目：%d", points.length());
    
23. 
    
24.                 AcGeVector3d fd,fd1,fd2;
    
25.                 pRay->getFirstDeriv(asPnt3d(pt),fd);
    
26.                 pRay->close();
    
27.                 AcDbRay *pRay1;
    
28.                 AcGePoint3d tmpPoint;
    
29.                 ads_name tmpE;
    
30.                 double min=0.00002;
    
31.                 double paramBegin,paramEnd,param,param1,param2,angle1,angle2;
    
32.                 for(int i=0;i<points.length();i++)
    
33.                 {
    
34.                         //acutPrintf("\n交点%d：%d,%d,%d",i,points[i].x,points[i].y,points[i].z);
    
35.                         //pPline->getFirstDeriv(points[i],fd1);
    
36. 
    
37.                         pPline->getStartParam(paramBegin);
    
38.                         pPline->getEndParam(paramEnd);
    
39.                         pPline->getParamAtPoint(points[i],param);
    
40. 
    
41.                         /*addline(points[i]);*/
    
42. 
    
43.                         if (param== paramBegin)
    
44.                         {
    
45.                                 param1=paramEnd-min;
    
46.                                 param2=paramBegin+min;
    
47.                         }
    
48.                         else if (param==  paramEnd )
    
49.                         {
    
50.                                 param1=paramEnd-min;
    
51.                                 param2=paramBegin+min;
    
52.                         }
    
53.                         else
    
54.                         {
    
55.                                 param1=param-min;
    
56.                                 param2=param+min;
    
57.                         }       
    
58.                         pPline->getFirstDeriv(param1,fd1);
    
59.                         pPline->getFirstDeriv(param2,fd2);       
    
60. 
    
61.                         /*acutPrintf("\n1 %f:%f",param1,fd1.angleTo(fd));       
    
62.                         acutPrintf(" %f %f %f",fd1.x,fd1.y,fd1.z);
    
63.                         acutPrintf("\n2 %f:%f",param2,fd2.angleTo(fd));               
    
64.                         acutPrintf(" %f %f %f",fd2.x,fd2.y,fd2.z);
    
65.                         acutPrintf("\n");*/                       
    
66. 
    
67.                         int valid=1;
    
68.                         if ((fd1.y>0)&&(fd2.y<0)||(fd1.y<0)&&(fd2.y>0))
    
69.                         {
    
70.                                 valid=0;
    
71.                                 acutPrintf("\n！！！");
    
72.                         }
    
73. 
    
74.                         /*acdbEntLast( tmpE);
    
75.                         acedGetPoint(NULL,"",tmpP);
    
76.                         acdbEntDel(tmpE);*/
    
77. 
    
78.                         if (valid==1)
    
79.                         {
    
80.                                 points1.append(points[i]);
    
81.                         }
    
82. 
    
83.                 }
    
84.                 //acutPrintf("\n交点数目：%d", points1.length());
    
85.                 if((points1.length()%2)==1)
    
86.                 {
    
87.                         acutPrintf("\nIN");
    
88.                         return Adesk::kTrue;
    
89.                 }
    
90.                 else
    
91.                 {
    
92.                         acutPrintf("\nOUT");
    
93.                         return Adesk::kFalse;
    
94.                 }
    
95.                
    
96. 
    
97. }
    

复制代码