函数分享：AcDbPolyline 节点除重

Love-Lisp

和大家分享一下我写的AcDbPolyline 节点除重函数！

1. //多段线除重
   
2. //输入 AcDbPolyline 、AcGeTol、angPrec = 角度偏差值//功能：AcDbPolyline 消除重复节点和同直线段或同圆弧上的节点
   
3. Acad::ErrorStatus CCurve::RemoveDumpPoly(AcDbPolyline *&pPolyine,AcGeTol Tol,double angPrec)
   
4. {
   
5.   if (!pPolyine->isKindOf(AcDbPolyline::desc()))
   
6.   {
   
7.     return Acad::eNotImplementedYet;
   
8.   }
   
9.   bool isclosed = pPolyine->isClosed();
   
10.   AcGePoint3d p0 ,p1,p2;
    
11.   AcGeVector3d v1,v2;
    
12.   int k =0;
    
13.   if (isclosed)
    
14.   {
    
15.     pPolyine->getPointAt(pPolyine->numVerts()-1,p0);
    
16.   }
    
17.   pPolyine->getPointAt(k,p1);
    
18.   pPolyine->getPointAt(k+1,p2);
    
19.   while (k < pPolyine->numVerts())
    
20.   {
    
21.     int n = k+1;
    
22.     if (n == pPolyine->numVerts())
    
23.     {
    
24.       n = 0;
    
25.     }
    
26.     int m = k - 1;
    
27.     if (m < 0)
    
28.     {
    
29.       m = pPolyine->numVerts()-1;
    
30.     }
    
31.     pPolyine->getPointAt(m,p0);
    
32.     pPolyine->getPointAt(k,p1);
    
33.     pPolyine->getPointAt(n,p2);
    
34.     if (p1.isEqualTo(p2,Tol))
    
35.     {
    
36.       pPolyine->removeVertexAt(k);
    
37.       if (n == 0)
    
38.       {
    
39.         pPolyine->setClosed(true);
    
40.       }
    
41.     }
    
42.     else
    
43.     {
    
44.       v1 =p1 - p0 ;
    
45.       v2 = p2 - p1;
    
46.       if ((v1.angleTo(v2)<= angPrec && k != 0) || (v1.angleTo(v2)<= angPrec && k == 0 && isclosed))
    
47.       {
    
48.         pPolyine->removeVertexAt(k);
    
49.       }
    
50.       else
    
51.       k++;
    
52.     }
    
53. 
    
54.   }
    
55.   //删除同弧段
    
56.     k = 0;
    
57.   double bulge1=0,bulge2=0;
    
58.   AcGePoint3d cp1,cp2;
    
59.   while (k < pPolyine->numVerts() -1)
    
60.   {
    
61.     pPolyine->getBulgeAt(k,bulge2);
    
62.     pPolyine->getPointAt(k,p1);
    
63.     if (!g_Equal(0,bulge2,0.00001))
    
64.     {
    
65.       if (g_Equal(0,bulge1,0.00001))
    
66.       {
    
67.         pPolyine->getPointAt(k,p0);
    
68.         pPolyine->getPointAt(k+1,p1);
    
69.         bulge1 = bulge2;
    
70.         //计算弧段的圆心
    
71.         double ang,rad;
    
72.         ang = 2 * atan(bulge1);
    
73.         rad = p0.distanceTo(p1) / (2 * sin(ang));
    
74.         ads_point cp;
    
75.         ang = (Pi2 - ang) + acutAngle(asDblArray(p0),asDblArray(p1));
    
76.         acutPolar(asDblArray(p0),ang,rad,cp);
    
77.         cp1 = asPnt3d(cp);
    
78.         k++;
    
79. 
    
80.       }
    
81.       else
    
82.       {
    
83.         pPolyine->getPointAt(k,p1);
    
84.         pPolyine->getPointAt(k+1,p2);
    
85.         //计算弧段的圆心
    
86.         double ang,rad,ang1,ang2 ;
    
87.         ang = 2 * atan(bulge2);
    
88.         rad = p1.distanceTo(p2) / (2 * sin(ang));
    
89.         ads_point cp;
    
90.         ang = (Pi2 - ang) + acutAngle(asDblArray(p1),asDblArray(p2));
    
91.         acutPolar(asDblArray(p1),ang,rad,cp);
    
92.         cp2 = asPnt3d(cp);
    
93.         if (cp1.isEqualTo(cp2,Tol))
    
94.         {
    
95.           pPolyine->removeVertexAt(k);
    
96.           //重新计算弓弦比
    
97.           ang1 = acutAngle(asDblArray(cp1),asDblArray(p0));
    
98.           ang2 = acutAngle(asDblArray(cp1),asDblArray(p2));
    
99.           if (bulge1 > 0)
    
100.           {
     
101.             ang = ang2 -ang1;
     
102.             if (ang > _2Pi)
     
103.               ang = ang - _2Pi;
     
104.             else if (ang < 0)
     
105.             ang = ang + _2Pi;
     
106.           }
     
107.           else
     
108.           {
     
109.             ang = ang1 - ang2;
     
110.             if (ang > _2Pi)
     
111.               ang = ang - _2Pi;
     
112.             else if (ang < 0)
     
113.               ang = ang + _2Pi;
     
114.           }
     
115.           ang = ang / 4; //包角
     
116.           if (bulge1 > 0)
     
117.           {
     
118.             bulge1 = tan(ang);
     
119.           }
     
120.           else
     
121.           {
     
122.             bulge1 = -1 * tan(ang);
     
123.           }
     
124.           pPolyine->setBulgeAt(k-1,bulge1);
     
125. 
     
126.            
     
127. 
     
128.         }
     
129.         else
     
130.         {
     
131.           p0 = p1;
     
132.           cp1 = cp2;
     
133.           bulge1 = bulge2;
     
134.           k++;
     
135.         }
     
136. 
     
137.       }
     
138.     }
     
139.     else
     
140.     {
     
141.       p0 = p1;
     
142.       bulge1 = 0;
     
143.       k++;
     
144.     }
     
145.   }
     
146.   if (pPolyine->isClosed() && pPolyine->numVerts() > 2)
     
147.   {
     
148.     int numVerts = pPolyine->numVerts();
     
149.      pPolyine->getBulgeAt(numVerts-1,bulge1);
     
150.      pPolyine->getBulgeAt(0,bulge2);
     
151.     if (!g_Equal(bulge1,0,0.00001) && !g_Equal(bulge2,0,0.00001))
     
152.     {
     
153.       pPolyine->getPointAt(numVerts-1,p0); /*最后一点*/
     
154.       pPolyine->getPointAt(0,p1); /*起始点*/
     
155.       pPolyine->getPointAt(1,p2);/* 第二点*/
     
156.       //计算p0p1弧段的圆心
     
157.       double ang,rad,ang1,ang2;
     
158.       ang = 2 * atan(bulge1);
     
159.       rad = p0.distanceTo(p1) / (2 * sin(ang));
     
160.       ads_point cp;
     
161.       ang = (Pi2 - ang) + acutAngle(asDblArray(p0),asDblArray(p1));
     
162.       acutPolar(asDblArray(p0),ang,rad,cp);
     
163.       cp1 = asPnt3d(cp);
     
164.       //计算p1p2弧段的圆心
     
165.       ang = 2 * atan(bulge2);
     
166.       rad = p1.distanceTo(p2) / (2 * sin(ang));
     
167.       ang = (Pi2 - ang) + acutAngle(asDblArray(p1),asDblArray(p2));
     
168.       acutPolar(asDblArray(p1),ang,rad,cp);
     
169.       cp2 = asPnt3d(cp);
     
170.       if (cp1.isEqualTo(cp2,Tol))
     
171.       {
     
172.         pPolyine->removeVertexAt(0);
     
173.         //重新计算弓弦比
     
174.         ang1 = acutAngle(asDblArray(cp1),asDblArray(p0));
     
175.         ang2 = acutAngle(asDblArray(cp1),asDblArray(p2));
     
176.         if (bulge1 > 0)
     
177.         {
     
178.           ang = ang2 -ang1;
     
179.           if (ang > _2Pi)
     
180.             ang = ang - _2Pi;
     
181.           else if (ang < 0)
     
182.             ang = ang + _2Pi;
     
183.         }
     
184.         else
     
185.         {
     
186.           ang = ang1 - ang2;
     
187.           if (ang > _2Pi)
     
188.             ang = ang - _2Pi;
     
189.           else if (ang < 0)
     
190.             ang = ang + _2Pi;
     
191.         }
     
192.         ang = ang / 4; //包角
     
193.         if (bulge1 > 0)
     
194.         {
     
195.           bulge1 = tan(ang);
     
196.         }
     
197.         else
     
198.         {
     
199.           bulge1 = -1 * tan(ang);
     
200.         }
     
201.         pPolyine->setBulgeAt(numVerts-2,bulge1);
     
202. 
     
203.       }
     
204.     }
     
205. 
     
206. 
     
207.   }
     
208.   return Acad::eOk;
     
209. }
     

//多段线除重<br /> //输入 AcDbPolyline 、AcGeTol、angPrec = 角度偏差值//功能：AcDbPolyline 消除重复节点和同直线段或同圆弧上的节点<br /> Acad::ErrorStatus CCurve::RemoveDumpPoly(AcDbPolyline *&pPolyine,AcGeTol Tol,double angPrec)<br /> {<br />   if (!pPolyine->isKindOf(AcDbPolyline::desc()))<br />   {<br />     return Acad::eNotImplementedYet;<br />   }<br />   bool isclosed = pPolyine->isClosed();<br />   AcGePoint3d p0 ,p1,p2;<br />   AcGeVector3d v1,v2;<br />   int k =0;<br />   if (isclosed)<br />   {<br />     pPolyine->getPointAt(pPolyine->numVerts()-1,p0);<br />   }<br />   pPolyine->getPointAt(k,p1);<br />   pPolyine->getPointAt(k+1,p2);<br />   while (k < pPolyine->numVerts())<br />   {<br />     int n = k+1;<br />     if (n == pPolyine->numVerts())<br />     {<br />       n = 0;<br />     }<br />     int m = k - 1;<br />     if (m < 0)<br />     {<br />       m = pPolyine->numVerts()-1;<br />     }<br />     pPolyine->getPointAt(m,p0);<br />     pPolyine->getPointAt(k,p1);<br />     pPolyine->getPointAt(n,p2);<br />     if (p1.isEqualTo(p2,Tol))<br />     {<br />       pPolyine->removeVertexAt(k);<br />       if (n == 0)<br />       {<br />         pPolyine->setClosed(true);<br />       }<br />     }<br />     else<br />     {<br />       v1 =p1 - p0 ;<br />       v2 = p2 - p1;<br />       if ((v1.angleTo(v2)<= angPrec && k != 0) || (v1.angleTo(v2)<= angPrec && k == 0 && isclosed))<br />       {<br />         pPolyine->removeVertexAt(k);<br />       }<br />       else<br />       k++;<br />     }<br /> <br />   }<br />   //删除同弧段<br />     k = 0;<br />   double bulge1=0,bulge2=0;<br />   AcGePoint3d cp1,cp2;<br />   while (k < pPolyine->numVerts() -1)<br />   {<br />     pPolyine->getBulgeAt(k,bulge2);<br />     pPolyine->getPointAt(k,p1);<br />     if (!g_Equal(0,bulge2,0.00001))<br />     {<br />       if (g_Equal(0,bulge1,0.00001))<br />       {<br />         pPolyine->getPointAt(k,p0);<br />         pPolyine->getPointAt(k+1,p1);<br />         bulge1 = bulge2;<br />         //计算弧段的圆心<br />         double ang,rad;<br />         ang = 2 * atan(bulge1);<br />         rad = p0.distanceTo(p1) / (2 * sin(ang));<br />         ads_point cp;<br />         ang = (Pi2 - ang) + acutAngle(asDblArray(p0),asDblArray(p1));<br />         acutPolar(asDblArray(p0),ang,rad,cp);<br />         cp1 = asPnt3d(cp);<br />         k++;<br /> <br />       }<br />       else<br />       {<br />         pPolyine->getPointAt(k,p1);<br />         pPolyine->getPointAt(k+1,p2);<br />         //计算弧段的圆心<br />         double ang,rad,ang1,ang2 ;<br />         ang = 2 * atan(bulge2);<br />         rad = p1.distanceTo(p2) / (2 * sin(ang));<br />         ads_point cp; <br />         ang = (Pi2 - ang) + acutAngle(asDblArray(p1),asDblArray(p2));<br />         acutPolar(asDblArray(p1),ang,rad,cp);<br />         cp2 = asPnt3d(cp);<br />         if (cp1.isEqualTo(cp2,Tol))<br />         {<br />           pPolyine->removeVertexAt(k);<br />           //重新计算弓弦比<br />           ang1 = acutAngle(asDblArray(cp1),asDblArray(p0));<br />           ang2 = acutAngle(asDblArray(cp1),asDblArray(p2));<br />           if (bulge1 > 0)<br />           {<br />             ang = ang2 -ang1;<br />             if (ang > _2Pi)<br />               ang = ang - _2Pi;<br />             else if (ang < 0)<br />             ang = ang + _2Pi;<br />           } <br />           else<br />           {<br />             ang = ang1 - ang2;<br />             if (ang > _2Pi)<br />               ang = ang - _2Pi;<br />             else if (ang < 0)<br />               ang = ang + _2Pi;<br />           }<br />           ang = ang / 4; //包角<br />           if (bulge1 > 0)<br />           {<br />             bulge1 = tan(ang);<br />           } <br />           else<br />           {<br />             bulge1 = -1 * tan(ang);<br />           }<br />           pPolyine->setBulgeAt(k-1,bulge1);<br /> <br />            <br /> <br />         }<br />         else<br />         {<br />           p0 = p1;<br />           cp1 = cp2;<br />           bulge1 = bulge2;<br />           k++;<br />         }<br /> <br />       }<br />     }<br />     else<br />     {<br />       p0 = p1;<br />       bulge1 = 0;<br />       k++;<br />     }<br />   }<br />   if (pPolyine->isClosed() && pPolyine->numVerts() > 2)<br />   {<br />     int numVerts = pPolyine->numVerts();<br />      pPolyine->getBulgeAt(numVerts-1,bulge1);<br />      pPolyine->getBulgeAt(0,bulge2);<br />     if (!g_Equal(bulge1,0,0.00001) && !g_Equal(bulge2,0,0.00001))<br />     {<br />       pPolyine->getPointAt(numVerts-1,p0); /*最后一点*/<br />       pPolyine->getPointAt(0,p1); /*起始点*/<br />       pPolyine->getPointAt(1,p2);/* 第二点*/<br />       //计算p0p1弧段的圆心<br />       double ang,rad,ang1,ang2;<br />       ang = 2 * atan(bulge1);<br />       rad = p0.distanceTo(p1) / (2 * sin(ang));<br />       ads_point cp;<br />       ang = (Pi2 - ang) + acutAngle(asDblArray(p0),asDblArray(p1));<br />       acutPolar(asDblArray(p0),ang,rad,cp);<br />       cp1 = asPnt3d(cp);<br />       //计算p1p2弧段的圆心<br />       ang = 2 * atan(bulge2);<br />       rad = p1.distanceTo(p2) / (2 * sin(ang));<br />       ang = (Pi2 - ang) + acutAngle(asDblArray(p1),asDblArray(p2));<br />       acutPolar(asDblArray(p1),ang,rad,cp);<br />       cp2 = asPnt3d(cp);<br />       if (cp1.isEqualTo(cp2,Tol))<br />       {<br />         pPolyine->removeVertexAt(0);<br />         //重新计算弓弦比<br />         ang1 = acutAngle(asDblArray(cp1),asDblArray(p0));<br />         ang2 = acutAngle(asDblArray(cp1),asDblArray(p2));<br />         if (bulge1 > 0)<br />         {<br />           ang = ang2 -ang1;<br />           if (ang > _2Pi)<br />             ang = ang - _2Pi;<br />           else if (ang < 0)<br />             ang = ang + _2Pi;<br />         } <br />         else<br />         {<br />           ang = ang1 - ang2;<br />           if (ang > _2Pi)<br />             ang = ang - _2Pi;<br />           else if (ang < 0)<br />             ang = ang + _2Pi;<br />         }<br />         ang = ang / 4; //包角<br />         if (bulge1 > 0)<br />         {<br />           bulge1 = tan(ang);<br />         } <br />         else<br />         {<br />           bulge1 = -1 * tan(ang);<br />         }<br />         pPolyine->setBulgeAt(numVerts-2,bulge1);<br /> <br />       }<br />     }<br /> <br /> <br />   }<br />   return Acad::eOk;<br /> }<br />

XDSoft

你算是真上路了，而且跑的还挺快{:soso_e113:}

Free-Lancer

以后就是 Love ARX 不是 Love Lisp 了:lol

Love-Lisp

XDSoft 发表于 2013-9-11 22:22
你算是真上路了，而且跑的还挺快

呵呵，谢谢夸奖，刚开始入门，发现要学的东西太多了！{:soso_e127:}
目前学着写了大概有二三十个函数了吧，主要是曲线方面的函数：
        {_T("XLRX-CURVEINTERS"), XL_CurveInters},/*计算交点 (XLRX-CURVEINTERS e1 e2) (XLRX-CURVEINTERS ss) */
        {_T("XLRX-SortListUp"), XL_QSortLispUp},/*表升序排序 （XLRX-SortListUp list)*/
        {_T("XLRX-SortPtsOnCurve"), XL_SortPtsOnCurve}, /*点沿曲线排序 (XLRX-SortPtsOnCurve curve PtList)*/
        {_T("XLRX-SplitCurves"), XL_BreakCurves}, /*按点表打断曲线 (XLRX-BreakCurves curve PtList)*/
        {_T("XLRX-SSBreak"), XL_SsIntersectBreak}, /*交点打断，参数(XLRX-SSBreak SS [tol])*/
        {_T("XLRX-SETPREC"), XL_SetGolblePrec}, /*设置点、向量的全局精度 (XLRX-SETPREC prec) */
        {_T("XLRX-GETSSBOX"), XL_GetSsBox}, /*返回选择集对象包围框 (XLRX-GETSSBOX ss) */
        {_T("XLRX-GETSAMPLEPT"), XL_Curvesampt},/*返回多段线节点 （XLRX-GETSAMPLEPT curve [容差值]) 返回值不含闭合点*/
        {_T("XLRX-Region2Poly"), XL_RegionToPoly},/*Region转LWPOLYLINE (XLRX-Region2Poly ename/ss [容差值])*/
        {_T("XLRX-DumpPoly"), XL_RemoveDumpPolyPoint},/*多段线除重 参数 (XLRX-DumpPoly ename [点的容差值] [角度容差值])*/
        {_T("XLRX-Curve2Region"), XL_CurveToRegion},//曲线转面域 (XLRX-Curve2Region ss [容差值])
        {_T("XLRX-SYS-ZoomWinDow"), XL_ZoomWinDow},//窗口显示（XLRX-SYS-ZoomWinDow p1 p2 [扩展比率 = 1.0])
        {_T("XLRX-SYS-GETVAR"), XL_GetVar},/*(xlrx-getvar sysname)*/
        {_T("xlrx-Curve-Centroid"), XL_GetCentroid}, /*计算封闭曲线的形心坐标 (xlrx-Curve-Centroid curve)*/
        {_T("xlrx-Pts2Region"), XL_MakeRegionFromPoints},/*点表绘制region*/
        {_T("xlrx-DumpLine"), XL_DumpLineArcCircle},/*直线消重 (xlrx-DumpLine ss [Prec] [blayer] [bLType]) */
        {_T("xlrx-PtinPoly"), XL_PtIsinPoly}, /*判断点是否在封闭曲线内，模拟点法 (xlrx-ptInPoly  curve pt [prec = 模拟点精度]) (xlrx-ptInPoly pts pt)*/
        {_T("xlrx-PtinPoly1"), XL_PtIsinPoly1}, /*判断点是否在封闭曲线内，切线法！ (xlrx-ptInPoly pt curve) (xlrx-ptInPoly pt pts)*/
        {_T("xlrx-ClockWiseP"), XL_ClockWiseP}, /*判断点表或曲线或3点是否顺时针 (xlrx-clockwisep pts) (xlrx-clockwisep curve ) (xlrx-clockwisep p1 p2 p3 )*/
        {_T("xlrx-Curve2Spline"), XL_CurveToSpline}, /*曲线转Spline (xlrx-Curve2Spline curve[ss])*/
        {_T("xlrx-PolyRelation"), XL_PolyRelation}, /*判断两封闭曲线之间的相互关系 (xlrx-PolyRelation curve1 curve2 [prec])*/
        {_T("xlrx-OnorInPoly"), XL_OnorInPoly}, /*判断封闭曲线1是否在封闭曲线2内 (xlrx-OnorInPoly poly1 poly2 [prec]) */
        {_T("xlrx-GetHull"), XL_GetHull}, //根据点表计算凸包 (xlrx-GetHull ptlist [DrawBox Flag])
        {_T("xlrx-Spl_Ell2LWPoly"), XL_Spl_Ell2LWPoly}, //LWPOLYLINE模拟SPline (xlrx-Spl_Ell2LWPoly spline [ArcFlag t= 圆弧模拟 nil = 直线模拟] [弦长距离 = 0.1])

Love-Lisp

Free-Lancer 发表于 2013-9-11 22:30
以后就是 Love ARX 不是 Love Lisp 了

lovearx 名字已经有老大用了！{:soso_e113:}

XDSoft

Love-Lisp 发表于 2013-9-11 22:39
呵呵，谢谢夸奖，刚开始入门，发现要学的东西太多了！
目前学着写了大概有二三十个函数了 ...

wow!!

可否给我寄来，有些我想完善到API里面{:soso_e163:}

Love-Lisp

XDSoft 发表于 2013-9-11 22:49
wow!!

可否给我寄来，有些我想完善到API里面

没有问题！
能够入老大法眼，无尚荣幸！{:soso_e113:}

XDSoft

Love-Lisp 发表于 2013-9-11 22:58
没有问题！
能够入老大法眼，无尚荣幸！

{:soso_e163:}收到。

newer

楼主，可否把你的API PLINE模拟的SPLINE,分别用折线模拟和弧线模拟的，带原来的SPLINE，传个DWG图上来吗？ 看看弧线模拟和PLINE模拟的区别。

Free-Lancer

newer 发表于 2013-9-12 16:30
楼主，可否把你的API PLINE模拟的SPLINE,分别用折线模拟和弧线模拟的，带原来的SPLINE，传个DWG图上来吗？  ...

最简单的实体应该是 Line 吧！ 模拟点好像更好的就是股票的曲线，偏差

Love-Lisp

newer 发表于 2013-9-12 16:30
楼主，可否把你的API PLINE模拟的SPLINE,分别用折线模拟和弧线模拟的，带原来的SPLINE，传个DWG图上来吗？  ...

直线模拟和弧线模拟的节点数是相同的，弧线模拟是将直线段变成圆弧，弧线的中点落在曲线上，这样看上去模拟线和原曲线更接近！