曲线选择集交点打断的源码

Love-Lisp

XDSoft 发表于 2013-8-24 22:02
看了下API的代码，里面多了次重复的选择集遍历到AcDbObjectIdArray数组，API代码为了考虑编程的方便和通 ...

深有同感！代码越通用，必然降低运行效率！对于非大量计算的代码，基本无所谓，无需考虑效率问题！对于需要大量计算的程序来说，还真不能为了考虑代码的简洁性，而损失了效率！

XDSoft

Love-Lisp 发表于 2013-8-24 21:52
原先计算交点是nXn次计算，
for (int i=0,i < len,i++) //循环n次
{

另外，API考虑的参数多，还有一种自身打断，你的 i=0 第二次i=i+1, 就会落下自身对自身的交点的处理。就是你省下的那些次数。代码可以在综合考虑完善下算法。

Love-Lisp

XDSoft 发表于 2013-8-24 22:21
另外，API考虑的参数多，还有一种自身打断，你的 i=0 第二次i=i+1, 就会落下自身对自身的交点的处理。就 ...

哦，XDRX_CURVE_INTERSECTBREAK函数对自己也进行了打断！但是发现XDRX_CURVE_INTERSECTBREAK 函数对于随便画的自交曲线，打断的结果毫无规律可循，在交点处的两条线只打断一条，多段线自身打断，有的每个节点间都打断，但是有些节点间却不能打断!我觉得自身打断应该是在自交点处将曲线打断才对！

XDSoft

Love-Lisp 发表于 2013-8-24 23:03
哦，XDRX_CURVE_INTERSECTBREAK函数对自己也进行了打断！但是发现XDRX_CURVE_INTERSECTBREAK 函数对于随 ...

自身打断处理的不是自交，API处理的是比如多段线，在每个顶点打断，但实体还是多段线。不改变实体类型（区别于EXPLODE）

API函数已经优化了，下个版本提供xdrx_curve_intersectbreak测试， 打开新图出现对话框提示的也解决了，谢谢你的测试和修改建议。

XDSoft

修改了代码，速度大大提高了。

liuyj

先放个改进版吧{:soso_e100:}

eachy

仅Line构成的Net，没有调用 Command break 的spline，看来 Lisp 效率还可以接受{:soso_e100:}
命令: TT
选择对象: 指定对角点: 找到 1236 个
选择对象:
14.473
命令:
命令: U INTELLIZOOM '.ZOOM
命令: SSBRK
选择Line,ellipse,arc,circle,spline,lwpolyline,polyline.....
选择对象: 指定对角点: 找到 1236 个
选择对象:
共 1236 个实体排序中，请稍候.....
Ok, 共耗时35.526秒.
命令:

漂浮不定

[ 本帖最后由 LoveArx 于 2016-11-27 12:54 编辑 ]\n\n在楼主的代码基础上优化了一下，当1000*1000网格，即2002条图元相交，速度提升了差不多一半。代码如下：

1. struct IntersectOnCurve
   
2. {
   
3.         AcGeDoubleArray IntersectParamat;
   
4. };
   
5. 
   
6. int doublecmp( const void *a , const void *b )
   
7. {
   
8. return *(double *)a > *(double *)b ? 1 : -1;
   
9. }
   
10. 
    
11. 
    
12. void WTUModifyEnt::IntersectBreakSS (ads_name ss,AcDbVoidPtrArray& curveSegments,ads_real rTol)
    
13. {
    
14.         Acad::ErrorStatus es;
    
15.         ads_name e1;
    
16.         AcGePoint3dArray nPoints;
    
17.         AcGeTol tol;
    
18.         tol.setEqualPoint(rTol);
    
19.         tol.setEqualVector(rTol);
    
20.         if (!(nPoints.isEmpty()))
    
21.         {
    
22.                 nPoints.setPhysicalLength(0);
    
23.                 nPoints.setLogicalLength(0);
    
24.         }
    
25.         long len;
    
26.         acedSSLength(ss,&len);
    
27.         //实体个数小于两个，直接返回值
    
28.         if (len < 2)
    
29.         {
    
30.                 if (len = 1)
    
31.                 {
    
32.                         AcDbEntity *pEnt;
    
33.                         AcDbObjectId Id;
    
34.                         acedSSName(ss,0,e1);
    
35.                         acdbGetObjectId(Id,e1);
    
36.                         es = acdbOpenObject(pEnt,Id,AcDb::kForRead);
    
37.                         if (es != Acad::eOk)
    
38.                         {
    
39.                                 pEnt->close();
    
40. 
    
41.                         }
    
42.                         curveSegments.append(pEnt);
    
43.                         pEnt->close();
    
44.                 }
    
45.                 return;
    
46.         }
    
47.         //将选择集实体转换到对象数组IdArray
    
48.         AcDbObjectIdArray IdArray;
    
49.         AcDbEntity *pEnt,*pEnt1;
    
50.         AcDbObjectId Id,Id1;
    
51.         AcDbCurve *pCurve,*pCurve1;
    
52.         for (long i =0;i < len;i++)
    
53.         {
    
54.                 acedSSName(ss,i,e1);
    
55.                 acdbGetObjectId(Id,e1);
    
56.                 es = acdbOpenObject(pEnt,Id,AcDb::kForRead);
    
57.                 if (es != Acad::eOk)
    
58.                 {
    
59.                         pEnt->close();
    
60. 
    
61.                 }
    
62.                 pCurve = AcDbCurve::cast(pEnt);
    
63.                 if (pCurve == NULL)
    
64.                 {
    
65.                         pEnt->close();
    
66.                         continue;
    
67.                 }
    
68.                 IdArray.append(Id);
    
69.                 pEnt->close();
    
70.         }
    
71.         //曲线总个数
    
72.         len = IdArray.length();
    
73. 
    
74.         //创建曲线的IntersectOnCurve结构数组,每个元素保存对应曲线上的交点Param,以AcGeDoubleArray保存
    
75.         IntersectOnCurve * intersectpt = new IntersectOnCurve[len];
    
76. 
    
77.         //初始化intersectpt
    
78.         for (long i = 0;i < len;i++)
    
79.         {
    
80.                 if (!(intersectpt.IntersectParamat.isEmpty()))
    
81.                 {
    
82.                         intersectpt.IntersectParamat.setPhysicalLength(0);
    
83.                         intersectpt.IntersectParamat.setLogicalLength(0);
    
84.                 }
    
85.         }
    
86. 
    
87.         for (long i = 0;i < len-1;i++)
    
88.         {
    
89.                 Id = IdArray.at(i);
    
90.                 acdbOpenObject(pEnt,Id,AcDb::kForRead);
    
91.                 pCurve = AcDbCurve::cast(pEnt);
    
92.                 AcGePoint3d pstart0,pend0;
    
93.                 pCurve->getStartPoint(pstart0);
    
94.                 pCurve->getEndPoint(pend0);
    
95.                 for (long j = i+1;j < len;j++)
    
96.                 {
    
97.                         Id1 = IdArray.at(j);
    
98.                         acdbOpenObject(pEnt1,Id1,AcDb::kForRead);
    
99.                         pCurve1 = AcDbCurve::cast(pEnt1);
    
100.                         AcGePoint3dArray nPts;
     
101.                         //计算交点
     
102.                         pEnt->intersectWith(pEnt1,AcDb::kOnBothOperands,nPts);
     
103.                         AcGePoint3d pstart1,pend1,pt;
     
104.                         pCurve1->getStartPoint(pstart1);
     
105.                         pCurve1->getEndPoint(pend1);      
     
106.                         for (int m = 0;m < nPts.length();m++)
     
107.                         {
     
108.                                 pt = nPts[m];
     
109.                                 if (!(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))
     
110.                                 {
     
111.                                         //交点添加到第i条曲线上
     
112.                                         double paramat;
     
113.                                         pCurve->getParamAtPoint(pt,paramat);
     
114.                                         intersectpt.IntersectParamat.append(paramat);
     
115.                                 }        
     
116.                                 if (!(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))
     
117.                                 {
     
118.                                         //交点添加到第j条曲线上
     
119.                                         //Pt3DArrays[j].append(pt);
     
120.                                         double paramat;
     
121.                                         pCurve1->getParamAtPoint(pt,paramat);
     
122.                                         intersectpt[j].IntersectParamat.append(paramat);
     
123.                                 }      
     
124.                         }
     
125.                         //判断pCurve1的端点是否在pCurve容差范围之内，添加到交点表第i条曲线上
     
126.                         pCurve->getClosestPointTo(pstart1,pt);
     
127.                         if (pstart1.isEqualTo(pt,tol) && !(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))
     
128.                         {
     
129.                                 double paramat;
     
130.                                 pCurve->getParamAtPoint(pt,paramat);
     
131.                                 intersectpt.IntersectParamat.append(paramat);
     
132.                         }
     
133.                         pCurve->getClosestPointTo(pend1,pt);
     
134.                         if (pend1.isEqualTo(pt,tol) && !(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))
     
135.                         {
     
136.                                 double paramat;
     
137.                                 pCurve->getParamAtPoint(pt,paramat);
     
138.                                 intersectpt.IntersectParamat.append(paramat);
     
139.                         }  
     
140.                         //判断pCurve的端点是否在pCurve1容差范围之内，添加到交点表第j条曲线上
     
141.                         pCurve1->getClosestPointTo(pstart0,pt);
     
142.                         if (pstart0.isEqualTo(pt,tol) && !(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))
     
143.                         {
     
144.                                 double paramat;
     
145.                                 pCurve1->getParamAtPoint(pt,paramat);
     
146.                                 intersectpt[j].IntersectParamat.append(paramat);
     
147.                         }
     
148.                         pCurve1->getClosestPointTo(pend0,pt);
     
149.                         if (pend0.isEqualTo(pt,tol) && !(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))
     
150.                         {
     
151.                                 double paramat;
     
152.                                 pCurve1->getParamAtPoint(pt,paramat);
     
153.                                 intersectpt[j].IntersectParamat.append(paramat);
     
154.                         }  
     
155.                         pEnt1->close();
     
156.                 }
     
157.                 pEnt->close();
     
158.         }
     
159.         for (long i = 0;i < len;i++)
     
160.         {
     
161.                 //取出第i条曲线pCurve
     
162.                 Id = IdArray.at(i);
     
163.                 acdbOpenObject(pEnt,Id,AcDb::kForRead);
     
164.                 pCurve = AcDbCurve::cast(pEnt);
     
165. 
     
166.                 //根据点intersectpt打断第i条曲线
     
167.                 if (!(intersectpt.IntersectParamat.isEmpty()))
     
168.                 {
     
169.                         AcGeDoubleArray nPts;
     
170.                         nPts.setPhysicalLength(0);
     
171.                         nPts.setLogicalLength(0);
     
172.                         long inptlen=intersectpt.IntersectParamat.length();
     
173.                         double *intdouble=new double[inptlen];
     
174.                         for (long j = 0;j < inptlen;j++)
     
175.                         {
     
176.                                 intdouble[j]=intersectpt.IntersectParamat.at(j);
     
177.                         }
     
178.                         qsort(intdouble,inptlen,sizeof(intdouble[0]),doublecmp);
     
179.                         for (long i = 0;i < inptlen;i++)
     
180.                         {
     
181.                                 nPts.append(intdouble);
     
182.                         }
     
183.                         delete []intdouble;
     
184.                         //len1打断曲线前curveSegments的长度
     
185.                         long len1 = curveSegments.length();
     
186.                         pCurve->upgradeOpen();
     
187.                         //按参数表params打断曲线pCurve
     
188.                         es = pCurve->getSplitCurves(nPts,curveSegments);
     
189.                         //len2打断曲线后curveSegments的长度
     
190.                         long len2 = curveSegments.length();
     
191. 
     
192.                         if ((len2 - len1) > 1)/*生成了新的曲线*/
     
193.                         {
     
194.                                 //acutPrintf(_T("\n共生成直线 %d  条！"),(len2 - len1));
     
195.                                 pCurve->erase();
     
196.                                 pEnt->close();
     
197.                                 for (int ii = len1; ii <len2; ii++)
     
198.                                 {
     
199.                                         AcDbEntity *pEntity = NULL;
     
200.                                         pEntity = (AcDbEntity *)curveSegments[ii];
     
201.                                         //将生成的曲线添加到数据库
     
202.                                         WTUcreateEnt::PostToModelSpace(pEntity);
     
203.                                 }
     
204. 
     
205.                         }
     
206.                         else /*没有生成新曲线*/
     
207.                         {
     
208.                                 pEnt->close();
     
209.                                 AcDbEntity *pEnt1;
     
210.                                 pEnt1  = (AcDbEntity *)curveSegments[(len2 - 1)];
     
211.                                 delete(pEnt1);
     
212.                                 curveSegments.setAt(len2 - 1,pCurve);
     
213.                         }
     
214.                 }
     
215.                 else /*曲线没有交点，无需打断*/
     
216.                 {      
     
217.                         pEnt->close();
     
218.                         curveSegments.append(pCurve);
     
219.                 }
     
220. 
     
221.         }
     
222.         delete []intersectpt;
     
223. }

struct IntersectOnCurve<br /> { <br />         AcGeDoubleArray IntersectParamat; <br /> };<br /> <br /> int doublecmp( const void *a , const void *b ) <br /> { <br /> return *(double *)a > *(double *)b ? 1 : -1;<br /> } <br /> <br /> <br /> void WTUModifyEnt::IntersectBreakSS (ads_name ss,AcDbVoidPtrArray& curveSegments,ads_real rTol)<br /> {<br />         Acad::ErrorStatus es;<br />         ads_name e1;<br />         AcGePoint3dArray nPoints;<br />         AcGeTol tol;<br />         tol.setEqualPoint(rTol);<br />         tol.setEqualVector(rTol);<br />         if (!(nPoints.isEmpty()))<br />         {<br />                 nPoints.setPhysicalLength(0);<br />                 nPoints.setLogicalLength(0);<br />         }<br />         long len;<br />         acedSSLength(ss,&len);<br />         //实体个数小于两个，直接返回值<br />         if (len < 2)<br />         {<br />                 if (len = 1)<br />                 {<br />                         AcDbEntity *pEnt;<br />                         AcDbObjectId Id;<br />                         acedSSName(ss,0,e1);<br />                         acdbGetObjectId(Id,e1);<br />                         es = acdbOpenObject(pEnt,Id,AcDb::kForRead);<br />                         if (es != Acad::eOk)<br />                         {<br />                                 pEnt->close();<br /> <br />                         }<br />                         curveSegments.append(pEnt);<br />                         pEnt->close();<br />                 }<br />                 return;<br />         }<br />         //将选择集实体转换到对象数组IdArray<br />         AcDbObjectIdArray IdArray;<br />         AcDbEntity *pEnt,*pEnt1;<br />         AcDbObjectId Id,Id1;<br />         AcDbCurve *pCurve,*pCurve1;<br />         for (long i =0;i < len;i++)<br />         {<br />                 acedSSName(ss,i,e1);<br />                 acdbGetObjectId(Id,e1);<br />                 es = acdbOpenObject(pEnt,Id,AcDb::kForRead);<br />                 if (es != Acad::eOk)<br />                 {<br />                         pEnt->close();<br /> <br />                 }<br />                 pCurve = AcDbCurve::cast(pEnt);<br />                 if (pCurve == NULL)<br />                 {<br />                         pEnt->close();<br />                         continue;<br />                 }<br />                 IdArray.append(Id);<br />                 pEnt->close();<br />         }<br />         //曲线总个数<br />         len = IdArray.length();<br /> <br />         //创建曲线的IntersectOnCurve结构数组,每个元素保存对应曲线上的交点Param,以AcGeDoubleArray保存<br />         IntersectOnCurve * intersectpt = new IntersectOnCurve[len];<br /> <br />         //初始化intersectpt<br />         for (long i = 0;i < len;i++)<br />         {<br />                 if (!(intersectpt.IntersectParamat.isEmpty()))<br />                 {<br />                         intersectpt.IntersectParamat.setPhysicalLength(0);<br />                         intersectpt.IntersectParamat.setLogicalLength(0);<br />                 }<br />         }<br /> <br />         for (long i = 0;i < len-1;i++)<br />         {<br />                 Id = IdArray.at(i);<br />                 acdbOpenObject(pEnt,Id,AcDb::kForRead);<br />                 pCurve = AcDbCurve::cast(pEnt);<br />                 AcGePoint3d pstart0,pend0;<br />                 pCurve->getStartPoint(pstart0);<br />                 pCurve->getEndPoint(pend0);<br />                 for (long j = i+1;j < len;j++)<br />                 {<br />                         Id1 = IdArray.at(j);<br />                         acdbOpenObject(pEnt1,Id1,AcDb::kForRead);<br />                         pCurve1 = AcDbCurve::cast(pEnt1);<br />                         AcGePoint3dArray nPts;<br />                         //计算交点<br />                         pEnt->intersectWith(pEnt1,AcDb::kOnBothOperands,nPts);<br />                         AcGePoint3d pstart1,pend1,pt;<br />                         pCurve1->getStartPoint(pstart1);<br />                         pCurve1->getEndPoint(pend1);      <br />                         for (int m = 0;m < nPts.length();m++)<br />                         {<br />                                 pt = nPts[m];<br />                                 if (!(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))<br />                                 {<br />                                         //交点添加到第i条曲线上<br />                                         double paramat;<br />                                         pCurve->getParamAtPoint(pt,paramat);<br />                                         intersectpt.IntersectParamat.append(paramat);<br />                                 }        <br />                                 if (!(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))<br />                                 {<br />                                         //交点添加到第j条曲线上<br />                                         //Pt3DArrays[j].append(pt);<br />                                         double paramat;<br />                                         pCurve1->getParamAtPoint(pt,paramat);<br />                                         intersectpt[j].IntersectParamat.append(paramat);<br />                                 }      <br />                         }<br />                         //判断pCurve1的端点是否在pCurve容差范围之内，添加到交点表第i条曲线上<br />                         pCurve->getClosestPointTo(pstart1,pt);<br />                         if (pstart1.isEqualTo(pt,tol) && !(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))<br />                         {<br />                                 double paramat;<br />                                 pCurve->getParamAtPoint(pt,paramat);<br />                                 intersectpt.IntersectParamat.append(paramat);<br />                         }<br />                         pCurve->getClosestPointTo(pend1,pt);<br />                         if (pend1.isEqualTo(pt,tol) && !(pstart0.isEqualTo(pt,tol) || pend0.isEqualTo(pt,tol)))<br />                         {<br />                                 double paramat;<br />                                 pCurve->getParamAtPoint(pt,paramat);<br />                                 intersectpt.IntersectParamat.append(paramat);<br />                         }  <br />                         //判断pCurve的端点是否在pCurve1容差范围之内，添加到交点表第j条曲线上<br />                         pCurve1->getClosestPointTo(pstart0,pt);<br />                         if (pstart0.isEqualTo(pt,tol) && !(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))<br />                         {<br />                                 double paramat;<br />                                 pCurve1->getParamAtPoint(pt,paramat);<br />                                 intersectpt[j].IntersectParamat.append(paramat);<br />                         }<br />                         pCurve1->getClosestPointTo(pend0,pt);<br />                         if (pend0.isEqualTo(pt,tol) && !(pstart1.isEqualTo(pt,tol) || pend1.isEqualTo(pt,tol)))<br />                         {<br />                                 double paramat;<br />                                 pCurve1->getParamAtPoint(pt,paramat);<br />                                 intersectpt[j].IntersectParamat.append(paramat);<br />                         }  <br />                         pEnt1->close();<br />                 }<br />                 pEnt->close();<br />         }<br />         for (long i = 0;i < len;i++)<br />         {<br />                 //取出第i条曲线pCurve<br />                 Id = IdArray.at(i);<br />                 acdbOpenObject(pEnt,Id,AcDb::kForRead);<br />                 pCurve = AcDbCurve::cast(pEnt); <br /> <br />                 //根据点intersectpt打断第i条曲线<br />                 if (!(intersectpt.IntersectParamat.isEmpty()))<br />                 {<br />                         AcGeDoubleArray nPts;<br />                         nPts.setPhysicalLength(0);<br />                         nPts.setLogicalLength(0);<br />                         long inptlen=intersectpt.IntersectParamat.length();<br />                         double *intdouble=new double[inptlen];<br />                         for (long j = 0;j < inptlen;j++)<br />                         {<br />                                 intdouble[j]=intersectpt.IntersectParamat.at(j);<br />                         }<br />                         qsort(intdouble,inptlen,sizeof(intdouble[0]),doublecmp);<br />                         for (long i = 0;i < inptlen;i++)<br />                         {<br />                                 nPts.append(intdouble);<br />                         }<br />                         delete []intdouble;<br />                         //len1打断曲线前curveSegments的长度<br />                         long len1 = curveSegments.length();<br />                         pCurve->upgradeOpen();<br />                         //按参数表params打断曲线pCurve<br />                         es = pCurve->getSplitCurves(nPts,curveSegments);<br />                         //len2打断曲线后curveSegments的长度<br />                         long len2 = curveSegments.length();<br /> <br />                         if ((len2 - len1) > 1)/*生成了新的曲线*/<br />                         {<br />                                 //acutPrintf(_T("\n共生成直线 %d  条！"),(len2 - len1));<br />                                 pCurve->erase();<br />                                 pEnt->close();<br />                                 for (int ii = len1; ii <len2; ii++)<br />                                 {<br />                                         AcDbEntity *pEntity = NULL;<br />                                         pEntity = (AcDbEntity *)curveSegments[ii];<br />                                         //将生成的曲线添加到数据库<br />                                         WTUcreateEnt::PostToModelSpace(pEntity);<br />                                 }<br /> <br />                         }<br />                         else /*没有生成新曲线*/<br />                         {<br />                                 pEnt->close();<br />                                 AcDbEntity *pEnt1;<br />                                 pEnt1  = (AcDbEntity *)curveSegments[(len2 - 1)];<br />                                 delete(pEnt1);<br />                                 curveSegments.setAt(len2 - 1,pCurve);<br />                         }<br />                 }<br />                 else /*曲线没有交点，无需打断*/<br />                 {      <br />                         pEnt->close();<br />                         curveSegments.append(pCurve);<br />                 }<br /> <br />         }<br />         delete []intersectpt;<br /> }

fyl891

没胡阅读权···············

chaolong

可惜级别不够，无法浏览

wxplodder

现在还看不到，过段时间才能看到，努力中

hwwrrbb

要求太高，看不了！

mh0801

我现在真的很需要啊  谢谢  

不懂暧

感谢分享，楼主真的太好啦啊

shaodeng000

不错，顶一个