objectarx点集沿曲线排序太慢了，有其他方法吗

zdqwy19

；点集vector <AcGePoint3d> pts沿曲线排序
AcDbCurve *pCurve;
bool sortcurve(AcGePoint3d &a, AcGePoint3d &b) //沿曲线排序比较函数
{
AcGePoint3d outa, outb;
double outpa, outpb;
pCurve->getClosestPointTo(a, outa);
pCurve->getClosestPointTo(b, outb);
pCurve->getParamAtPoint(outa, outpa);
pCurve->getParamAtPoint(outb, outpb);
return outpa < outpb;
}
// 提示用户选择样条参照曲线  
  ads_name ent; ads_point pt;
  if(acedEntSel(ACRX_T("\n选择参照曲线:"), ent, pt) == RTNORM)
  {
   AcDbObjectId objId;
   acdbGetObjectId(objId, ent);
   // 获得指向当前元素的指针;
   AcDbEntity *pEnt;
   acdbOpenAcDbEntity(pEnt, objId, AcDb::kForRead);

   if (pEnt->isKindOf(AcDbCurve::desc()))//曲线;
   {
    pCurve = AcDbCurve::cast(pEnt);
    sort(pts.begin(), pts.end(), sortcurve );
    pCurve->close();
   }
   pEnt->close();
  }
这样可以沿曲线同向排序，就是速度太慢了，几十个点CAD的界面都看到非常明显的停顿，有其他方法可以加速吗？

newer

代码不全？

pts 从哪里来的？ sort函数？

你排序用的什么算法？

用快速排序改你的算法。

jiangnanemail

分析一下效率差的问题出在哪里，如果是算法的问题，可以使用 快速排序改你的算法；如果是“bool sortcurve(AcGePoint3d &a, AcGePoint3d &b) //沿曲线排序比较函数”这个比较函数效率差的的话，应当修改你的算法。我估计还是第二个原因导致的效率差。

zdqwy19

是bool sortcurve(AcGePoint3d &a, AcGePoint3d &b)的原因，沿曲线排序objectarx中有别的函数可以用吗？

LoveArx

zdqwy19 发表于 2017-1-16 10:18

问题不在sortcurve,而是前面的sort, sortcurve仅仅是值的比较而已，自己写排序算法，用快速排序。建议用 排序 搜下论坛看看。

marting

AEO版主写过一个 AcArray的快速排序模板，你去找找看看。

Highflybird

你先把点的参数存起来，不要在排序的比较函数中用到曲线的函数，这样就会节省大量的时间。节省的时间简直不止一个数量级.

csharp

1. 
   
2. typedef struct pointInfo
   
3. {
   
4.         double param;
   
5.         AcGePoint3d point;
   
6.         bool operator < ( pointInfo &ptInfo)
   
7.         {
   
8.                 return param < ptInfo.param;
   
9. 
   
10.         }
    
11.         bool operator > ( pointInfo &ptInfo)  
    
12.         {
    
13.                 return param > ptInfo.param;
    
14.         }
    
15. 
    
16. }PointInfo;
    
17. 
    
18. static void asdkMyGroupsSortpoints()
    
19.         {
    
20.                 // 提示用户选择样条参照曲线  
    
21.                 ads_name ent;
    
22.                 ads_point pt;
    
23.                 if(acedEntSel(ACRX_T("\n选择参照曲线:"), ent, pt) == RTNORM)
    
24.                 {
    
25.                         AcDbObjectId objId;
    
26.                         acdbGetObjectId(objId, ent);
    
27.                         // 获得指向当前元素的指针;
    
28. 
    
29.                         AcDbObjectPointer <AcDbCurve> pCurve(objId, AcDb::kForRead);
    
30.                         if(pCurve.openStatus()!=Acad::eOk)
    
31.                                 return;
    
32.                         ads_name ss;
    
33.                         resbuf *filter=acutBuildList(RTDXF0,_T("POINT"),RTNONE);
    
34.                         if (acedSSGet(NULL,NULL,NULL,filter,ss) != RTNORM)
    
35.                         {
    
36.                                 acutRelRb(filter);
    
37.                                 return;
    
38.                         }
    
39.                         acutRelRb(filter);
    
40.                         AcGePoint3dArray pts;
    
41.                         long slen;
    
42.                         acedSSLength(ss,&slen);
    
43.                         ads_name pnt;               
    
44.                         double pam;
    
45.                         Acad::ErrorStatus es;
    
46.                         std::vector<struct pointInfo> myvec;
    
47.                         struct pointInfo ptInfo;
    
48.                         for (long i=0;i<slen;i++)
    
49.                         {
    
50.                                 acedSSName(ss,i,pnt);
    
51.                                 acdbGetObjectId(objId,pnt);
    
52.                                 AcDbObjectPointer<AcDbPoint> pPoint(objId,AcDb::kForRead);
    
53.                                 AcGePoint3d pt(pPoint->position()),ptOn;
    
54.                                 es=pCurve->getClosestPointTo(pt,ptOn);
    
55.                                 if (es!=Acad::eOk)
    
56.                                 {
    
57.                                         continue;
    
58.                                 }
    
59.                                 es=pCurve->getParamAtPoint(ptOn,pam);
    
60.                                 if (es!=Acad::eOk)
    
61.                                 {
    
62.                                         continue;
    
63.                                 }
    
64.                                
    
65.                                 ptInfo.param=pam;
    
66.                                 ptInfo.point=pt;
    
67. 
    
68.                                 myvec.push_back(ptInfo);
    
69.                         }
    
70.                         std::sort(myvec.begin(),myvec.end());
    
71.                         std::vector<struct pointInfo>::iterator itor=myvec.begin();
    
72.                         for (;itor!=myvec.end();itor++)
    
73.                         {
    
74.                                 acutPrintf(_T("\n%.3f,%.3f,%.3f"),((*itor).point).x,((*itor).point).y,((*itor).point).z);
    
75.                         }                       
    
76.                 }
    

复制代码

Highflybird

arx的曲线操作函数特别费时，尤其是对样条曲线。所以在涉及到这类函数时，要小心使用。

zdqwy19

谢谢Highflybird和csharp的解答；样条曲线我感觉模拟比较理想，不容易出现飞点现象。我刚才按照二位的提示修改了一下代码，500个点排序用时3秒可以接受；比我原先的代码快多了！谢谢！

iLisp

3s还是慢啊！

marting

zdqwy19 发表于 2017-1-16 19:36

几十万个点，3S还差不多，你还是去好好改改快速排序的算法。

aeo

把点的参数求出来
按参数排序
参数反求点
。。。。。。。。。。。。。
如果你的点求了一点都不能动（参数反求点有误差)
那你参数按位置排序 (相当于vl-sort-i)

你如果直接按点排序，每次都要求参数，那还了得，吓死人了