优化多段线的倒角

雪山飞狐_lzh

在QQ群为风吹桥动改写优化多段线的倒角程序，
原代码使用DB库的曲线并用纯数学计算方式

3.         /// <summary>
4.         /// 通过一个顶点将多义线倒圆角,本程序不支持多义线起点和终点处倒圆角
5.         /// </summary>
6.         /// <param name="polyline">多义线对象</param>
7.         /// <param name="vertexIndex">顶点索引号</param>
8.         /// <param name="radius">圆弧半径</param>
9.         public static void FilletatVertex(Polyline polyline, int vertexIndex, double radius)
10.         {
11.             Database db = HostApplicationServices.WorkingDatabase;//获取数据库
12.             Editor ed = Application.DocumentManager.MdiActiveDocument.Editor;//获取editor对象
13.             using (Transaction trans = db.TransactionManager.StartTransaction())//开始事务     
14.             {
15.                 try
16.                 {
17.                     //判断顶点索引号是否属于多义线；
18.                     if (vertexIndex <= 0 || vertexIndex >= polyline.NumberOfVertices)
19.                     {
20.                         ed.WriteMessage("\n顶点无效。");
21.                         return;
22.                     }
23.                     else if ((polyline.GetBulgeAt(vertexIndex - 1) != 0) && (polyline.GetBulgeAt(vertexIndex) != 0))
24.                     {
25.                         ed.WriteMessage("\n顶点两侧不是直线。");
26.                         return;
27.                     }
28.                     else
29.                     {
30.                         Point2d pnt = polyline.GetPoint2dAt(vertexIndex);
31.                         Vector2d vec1 = polyline.GetPoint2dAt(vertexIndex - 1) - pnt;//由顶点与前一个顶点组成的向量
32.                         Vector2d vec2 = polyline.GetPoint2dAt(vertexIndex + 1) - pnt;//由顶点与后一个顶点组成的向量
33.                         double angle = vec1.GetAngleTo(vec2);//获得角点两侧线段间的角度变化
34.                         ed.WriteMessage(angle.ToString());
35.                         ed.WriteMessage("\n" + vec1.Length.ToString());
36.                         ed.WriteMessage("\n" + (radius / Math.Tan(angle / 2)).ToString());
37.                         ed.WriteMessage("\n" + vec2.Length.ToString());
38.                         ed.WriteMessage("\n" + (radius / Math.Tan(angle / 2)).ToString());
39.                         if ((vec1.Length < radius / Math.Tan(angle / 2)) || (vec2.Length < radius / Math.Tan(angle / 2)))
40.                         {
41.                             ed.WriteMessage("\n半径太大。");
42.                             return;
43.                         }
44.                         else
45.                         {
46.                             double lentemp = radius / Math.Tan(angle / 2);
47.                             Point2d addedVertex1 = pnt + vec1.GetNormal() * lentemp;
48.                             Point2d addedVertex2 = pnt + vec2.GetNormal() * lentemp;
49.                             polyline.RemoveVertexAt(vertexIndex);
50.                             polyline.AddVertexAt(vertexIndex, addedVertex1, 0, 0, 0);
51.                             polyline.AddVertexAt(vertexIndex + 1, addedVertex2, 0, 0, 0);
52.                             //polyline.AddVertexAt(vertexIndex, addedVertex1.Convert2d(new Plane()), 0, 0, 0);
53.                             //polyline.AddVertexAt(vertexIndex + 1, addedVertex2.Convert2d(new Plane()), 0, 0, 0);
55.                             //根据前后三顶点判断该顶点处多义线的前进方向是顺时针还是逆时针
56.                             CircularArc2d cir = new CircularArc2d(polyline.GetPoint2dAt(vertexIndex - 1), pnt, polyline.GetPoint2dAt(vertexIndex + 1));
57.                             double ang = Math.PI - angle;
58.                             double bulge = Math.Abs(Math.Tan(ang / 4));
59.                             if (cir.IsClockWise)
60.                             {
61.                                 polyline.SetBulgeAt(vertexIndex, -1 * bulge);//若为顺时针，则凸度为负；否则为正
62.                             }
63.                             else
64.                             {
65.                                 polyline.SetBulgeAt(vertexIndex, bulge);
66.                             }
67.                         }
68.                     }
69.                 }
70.                 catch (Autodesk.AutoCAD.Runtime.Exception e)
71.                 {
72.                     ed.WriteMessage(e.ToString());
73.                     trans.Abort();
74.                 }
75.                 finally
76.                 {
77.                     trans.Commit();
78.                 }
79.             }
80.         }

        /// <summary>         /// 通过一个顶点将多义线倒圆角,本程序不支持多义线起点和终点处倒圆角         /// </summary>         /// <param name="polyline">多义线对象</param>         /// <param name="vertexIndex">顶点索引号</param>         /// <param name="radius">圆弧半径</param>         public static void FilletatVertex(Polyline polyline, int vertexIndex, double radius)         {             Database db = HostApplicationServices.WorkingDatabase;//获取数据库             Editor ed = Application.DocumentManager.MdiActiveDocument.Editor;//获取editor对象             using (Transaction trans = db.TransactionManager.StartTransaction())//开始事务                  {                 try                 {                     //判断顶点索引号是否属于多义线；                     if (vertexIndex <= 0 || vertexIndex >= polyline.NumberOfVertices)                     {                         ed.WriteMessage("\n顶点无效。");                         return;                     }                     else if ((polyline.GetBulgeAt(vertexIndex - 1) != 0) && (polyline.GetBulgeAt(vertexIndex) != 0))                     {                         ed.WriteMessage("\n顶点两侧不是直线。");                         return;                     }                     else                     {                         Point2d pnt = polyline.GetPoint2dAt(vertexIndex);                         Vector2d vec1 = polyline.GetPoint2dAt(vertexIndex - 1) - pnt;//由顶点与前一个顶点组成的向量                         Vector2d vec2 = polyline.GetPoint2dAt(vertexIndex + 1) - pnt;//由顶点与后一个顶点组成的向量                         double angle = vec1.GetAngleTo(vec2);//获得角点两侧线段间的角度变化                         ed.WriteMessage(angle.ToString());                         ed.WriteMessage("\n" + vec1.Length.ToString());                         ed.WriteMessage("\n" + (radius / Math.Tan(angle / 2)).ToString());                         ed.WriteMessage("\n" + vec2.Length.ToString());                         ed.WriteMessage("\n" + (radius / Math.Tan(angle / 2)).ToString());                         if ((vec1.Length < radius / Math.Tan(angle / 2)) || (vec2.Length < radius / Math.Tan(angle / 2)))                         {                             ed.WriteMessage("\n半径太大。");                             return;                         }                         else                         {                             double lentemp = radius / Math.Tan(angle / 2);                             Point2d addedVertex1 = pnt + vec1.GetNormal() * lentemp;                             Point2d addedVertex2 = pnt + vec2.GetNormal() * lentemp;                             polyline.RemoveVertexAt(vertexIndex);                             polyline.AddVertexAt(vertexIndex, addedVertex1, 0, 0, 0);                             polyline.AddVertexAt(vertexIndex + 1, addedVertex2, 0, 0, 0);                             //polyline.AddVertexAt(vertexIndex, addedVertex1.Convert2d(new Plane()), 0, 0, 0);                             //polyline.AddVertexAt(vertexIndex + 1, addedVertex2.Convert2d(new Plane()), 0, 0, 0);                             //根据前后三顶点判断该顶点处多义线的前进方向是顺时针还是逆时针                             CircularArc2d cir = new CircularArc2d(polyline.GetPoint2dAt(vertexIndex - 1), pnt, polyline.GetPoint2dAt(vertexIndex + 1));                             double ang = Math.PI - angle;                             double bulge = Math.Abs(Math.Tan(ang / 4));                             if (cir.IsClockWise)                             {                                 polyline.SetBulgeAt(vertexIndex, -1 * bulge);//若为顺时针，则凸度为负；否则为正                             }                             else                             {                                 polyline.SetBulgeAt(vertexIndex, bulge);                             }                         }                     }                 }                 catch (Autodesk.<a href="http://bbs.mjtd.com/forum-41-1.html" target="_blank" class="relatedlink">AutoCAD</a>.Runtime.Exception e)                 {                     ed.WriteMessage(e.ToString());                     trans.Abort();                 }                 finally                 {                     trans.Commit();                 }             }         }

joyhai

你好，ToCurve() 是自定义的函数吗？ 我这也提示缺这个函数，请指教，谢谢。

雪山飞狐_lzh

感觉直接用Ge库的曲线应该也可以做，改写代码如下

3.         /// <summary>
4.         /// 为优化多段线倒角
5.         /// </summary>
6.         /// <param name="polyline">优化多段线</param>
7.         /// <param name="index">顶点索引号</param>
8.         /// <param name="radius">倒角半径</param>
9.         /// <param name="isFillet">倒角类型</param>
10.         public static void ChamferAt(this Polyline polyline, int index, double radius, bool isFillet)
11.         {
12.             if (index < 1 || index > polyline.NumberOfVertices - 2)
13.                 throw new System.Exception("错误的索引号");
14.             if (polyline.GetSegmentType(index - 1) != SegmentType.Line || polyline.GetSegmentType(index) != SegmentType.Line)
15.                 throw new System.Exception("非直线段不能倒角");
16.             //获取当前索引号的前后两段直线,并组合为Ge复合曲线
17.             Curve3d[] c3ds =
18.                 new Curve3d[]
19.                     {
20.                         polyline.GetLineSegmentAt(index - 1),
21.                         polyline.GetLineSegmentAt(index)
22.                     };
23.             var cc3d = new CompositeCurve3d(c3ds);
24.             //试倒直角
25.             //子曲线的个数有三种情况:
26.             //1、=3时倒角方向正确
27.             //2、=2时倒角方向相反
28.             //3、=0或为直线时失败
29.             c3ds =
30.                 cc3d.GetTrimmedOffset
31.                 (
32.                     radius,
33.                     Vector3d.ZAxis,
34.                     OffsetCurveExtensionType.Chamfer
35.                 );
36.             if (c3ds.Length > 0 && c3ds[0] is CompositeCurve3d)
37.             {
38.                     var newcc3d = c3ds[0] as CompositeCurve3d;
39.                     c3ds = newcc3d.GetCurves();
40.                     if (c3ds.Length == 3)
41.                     {
42.                         c3ds =
43.                             cc3d.GetTrimmedOffset
44.                             (
45.                                 -radius,
46.                                 Vector3d.ZAxis,
47.                                 OffsetCurveExtensionType.Chamfer
48.                             );
49.                         if (c3ds.Length == 0 || c3ds[0] is LineSegment3d)
50.                         {
51.                             throw new System.Exception("倒角半径过大");
52.                         }
53.                     }
54.                     else if (c3ds.Length == 2)
55.                     {
56.                         radius = -radius;
57.                     }
58.             }
59.             else
60.             {
61.                 throw new System.Exception("倒角半径过大");
62.             }
63.             //GetTrimmedOffset会生成倒角+偏移，故先反方向倒角,再倒回
64.             c3ds =
65.                 cc3d.GetTrimmedOffset
66.                 (
67.                     -radius,
68.                     Vector3d.ZAxis,
69.                     OffsetCurveExtensionType.Extend
70.                 );
71.             OffsetCurveExtensionType type =
72.                 isFillet ?
73.                 OffsetCurveExtensionType.Fillet : OffsetCurveExtensionType.Chamfer;
74.             c3ds =
75.                 c3ds[0].GetTrimmedOffset
76.                 (
77.                     radius,
78.                     Vector3d.ZAxis,
79.                     type
80.                 );
81.             //将结果Ge曲线转为Db曲线,并将相关的数值反映到原曲线
82.             Polyline plTemp = c3ds[0].ToCurve() as Polyline;
83.             polyline.RemoveVertexAt(index);
84.             polyline.AddVertexAt(index, plTemp.GetPoint2dAt(1), plTemp.GetBulgeAt(1), 0, 0);
85.             polyline.AddVertexAt(index + 1, plTemp.GetPoint2dAt(2), 0, 0, 0);
86.         }
90.         [CommandMethod("CTest")]
91.         public void test()
92.         {
93.             var db = HostApplicationServices.WorkingDatabase;
94.             var doc = Application.DocumentManager.GetDocument(db);
95.             var ed = doc.Editor;
96.             var resEnt = ed.GetEntity("\n请选择优化多段线:");
97.             var resInt = ed.GetInteger("\n请输入索引号:");
98.             bool isFillet = false;
99.             PromptDoubleOptions opt = new PromptDoubleOptions("\n请输入倒角半径:");
100.             opt.Keywords.Add("F");
101.             var resDbl = ed.GetDouble(opt);
102.             if (resDbl.Status == PromptStatus.Keyword)
103.             {
104.                 isFillet = true;
105.                 resDbl = ed.GetDouble(opt);
106.             }
107.             using (Transaction tr = db.TransactionManager.StartTransaction())
108.             {
109.                 Polyline pl = resEnt.ObjectId.GetObject(OpenMode.ForWrite) as Polyline;
110.                 ChamferAt(pl, resInt.Value, resDbl.Value, isFillet);
111.                 tr.Commit();
112.             }
113.         }

        /// <summary>         /// 为优化多段线倒角         /// </summary>         /// <param name="polyline">优化多段线</param>         /// <param name="index">顶点索引号</param>         /// <param name="radius">倒角半径</param>         /// <param name="isFillet">倒角类型</param>         public static void ChamferAt(this Polyline polyline, int index, double radius, bool isFillet)         {             if (index < 1 || index > polyline.NumberOfVertices - 2)                 throw new System.Exception("错误的索引号");             if (polyline.GetSegmentType(index - 1) != SegmentType.Line || polyline.GetSegmentType(index) != SegmentType.Line)                 throw new System.Exception("非直线段不能倒角");             //获取当前索引号的前后两段直线,并组合为Ge复合曲线             Curve3d[] c3ds =                 new Curve3d[]                     {                         polyline.GetLineSegmentAt(index - 1),                         polyline.GetLineSegmentAt(index)                     };             var cc3d = new CompositeCurve3d(c3ds);             //试倒直角             //子曲线的个数有三种情况:             //1、=3时倒角方向正确             //2、=2时倒角方向相反             //3、=0或为直线时失败             c3ds =                 cc3d.GetTrimmedOffset                 (                     radius,                     Vector3d.ZAxis,                     OffsetCurveExtensionType.Chamfer                 );             if (c3ds.Length > 0 && c3ds[0] is CompositeCurve3d)             {                     var newcc3d = c3ds[0] as CompositeCurve3d;                     c3ds = newcc3d.GetCurves();                     if (c3ds.Length == 3)                     {                         c3ds =                             cc3d.GetTrimmedOffset                             (                                 -radius,                                 Vector3d.ZAxis,                                 OffsetCurveExtensionType.Chamfer                             );                         if (c3ds.Length == 0 || c3ds[0] is LineSegment3d)                         {                             throw new System.Exception("倒角半径过大");                         }                     }                     else if (c3ds.Length == 2)                     {                         radius = -radius;                     }             }             else             {                 throw new System.Exception("倒角半径过大");             }             //GetTrimmedOffset会生成倒角+偏移，故先反方向倒角,再倒回             c3ds =                 cc3d.GetTrimmedOffset                 (                     -radius,                     Vector3d.ZAxis,                     OffsetCurveExtensionType.Extend                 );             OffsetCurveExtensionType type =                 isFillet ?                 OffsetCurveExtensionType.Fillet : OffsetCurveExtensionType.Chamfer;             c3ds =                 c3ds[0].GetTrimmedOffset                 (                     radius,                     Vector3d.ZAxis,                     type                 );             //将结果Ge曲线转为Db曲线,并将相关的数值反映到原曲线             Polyline plTemp = c3ds[0].ToCurve() as Polyline;             polyline.RemoveVertexAt(index);             polyline.AddVertexAt(index, plTemp.GetPoint2dAt(1), plTemp.GetBulgeAt(1), 0, 0);             polyline.AddVertexAt(index + 1, plTemp.GetPoint2dAt(2), 0, 0, 0);         }         [CommandMethod("CTest")]         public void test()         {             var db = HostApplicationServices.WorkingDatabase;             var doc = Application.DocumentManager.GetDocument(db);             var ed = doc.Editor;             var resEnt = ed.GetEntity("\n请选择优化多段线:");             var resInt = ed.GetInteger("\n请输入索引号:");             bool isFillet = false;             PromptDoubleOptions opt = new PromptDoubleOptions("\n请输入倒角半径:");             opt.Keywords.Add("F");             var resDbl = ed.GetDouble(opt);             if (resDbl.Status == PromptStatus.Keyword)             {                 isFillet = true;                 resDbl = ed.GetDouble(opt);             }             using (Transaction tr = db.TransactionManager.StartTransaction())             {                 Polyline pl = resEnt.ObjectId.GetObject(OpenMode.ForWrite) as Polyline;                 ChamferAt(pl, resInt.Value, resDbl.Value, isFillet);                 tr.Commit();             }         }

命令行的输入：
命令: ctest
请选择优化多段线:
请输入索引号: 3
请输入倒角半径 [F]: f
请输入倒角半径 [F]: 5
命令: ctest
请选择优化多段线:
请输入索引号: 2
请输入倒角半径 [F]: 5

倒角效果：

single-yu

顶飞狐版主呀，学习了！

hmxmylove

lincoln_ma

楼主强大，顶

雪山飞狐_lzh

修正试倒角时的Bug，2楼代码已更改

wenxi0504

太好的东西，正是需要的资料。

sicky111

要是有LISP代码就好了。

sdaulj

路过并表示强烈支持！