曲线按特定方式拟合

风树

如何实现用多段线拟合样条曲线和圆的lisp
希望能实现的是：将一条曲线用一条多段线近似表示，多段线由与x轴夹角为0,x,90度（x大于0,小于90）的线段（共四种形态）构成
参数输入为：
1：k,，k单位为1个图形单位，不小于0,代表构成多段线的子线段的最小长度，
子线段等长，非等长均可。
2:x

拟合精度控制不考虑数学指标，通过调整参数可以基本在视觉上过得去即可。

自贡黄明儒

这个早就人写过，ET工具也有，晓东上的newer版主写了很多曲线相互转化

风树

自贡黄明儒 发表于 2014-7-12 21:16
这个早就人写过，ET工具也有，晓东上的newer版主写了很多曲线相互转化

固定角度的

wzg356

很抱歉，记不清原帖子的位置了
适用范围比较大
曲线、椭圆、圆、圆弧等转换的数学指标按相似度调整，如果是直线、未拟合的多线段不受相似度影响

2. ;;将各种线段转为多线段，将多线段上的弧线段也转为多线段
   
3. ;;(CurveToPoline en);;参数：图元名;;示例(CurveToPoline (car(entsel "\n选择对象: ")))
   
4. ;;适用"line,lwpolyline,polyline,spline,circle,arc,ellipse"
   
5. (defun CurveToPoline ( en / c2 C ob a n pt D l m  q p b o r F )
   
6.         (vl-load-com)                        
   
7.         (setq C2 0.65)        ;;可理解为相似度，不要小于0.65，不能为1
   
8.         (setq C(/(- 1.0 C2)2))        
   
9.     (setq ob (vlax-ename->vla-object en)
   
10.               a 0
    
11.               n 1
    
12.               pt nil
    
13.               D (vla-get-objectname ob)
    
14.               l (vlax-curve-getdistatparam ob(vlax-curve-getendparam ob))
    
15.              m (cond((= D"AcDb3dPolyline")3)((= D"AcDb2dPolyline")2)((= D"AcDbPolyline")2))
    
16.         )
    
17.         (if(or(= D"AcDb3dPolyline")(= D"AcDbPolyline"))
    
18.           (if(vl-catch-all-error-p(vl-catch-all-apply 'vla-get-Coordinates(list ob)))(setq Pt(StarEnd ob))
    
19.                   (progn
    
20.                           (setq q(vlax-safearray->list(vlax-variant-value(vla-get-Coordinates ob))))
    
21.                           (repeat(/(length q)m)
    
22.                                     (cond
    
23.                                           ((= m 2)(setq p(list(nth a q)(nth(+ a 1)q))))
    
24.                                           ((= m 3)(setq p(list(nth a q)(nth(+ a 1)q)(nth(+ a 2)q))))
    
25.                                         )
    
26.                                         (setq a(+ a m)
    
27.                                                   Pt(if(or(null Pt)(>(distance p(cdr(car(reverse Pt))))0.00005))
    
28.                                                       (append Pt(list(append(list(vlax-curve-getparamatpoint ob p))p)))pt)
    
29.                                         )
    
30.                                 )
    
31.                         )
    
32.                 )
    
33.                 (setq Pt(StarEnd ob))
    
34.         )
    
35.         (while(< n(length Pt))
    
36.                 (setq m(1- n)
    
37.                           a(vlax-curve-getdistAtparam ob(car(nth m Pt)))
    
38.                           b(-(vlax-curve-getdistAtparam ob(car(nth n Pt)))a)
    
39.                 )
    
40.                 (if(<= b 0)(setq b(- l a)))
    
41.                 (setq o(vlax-curve-getparamatdist ob(+(/ b 2)a))
    
42.                           r(append(list o)(vlax-curve-getpointatparam ob o))
    
43.                           F(count-f ob (* b c2)(cdr(nth m Pt))(cdr r)(cdr(nth n Pt))c)
    
44.                 )
    
45.                 (cond
    
46.                         ((= F 1)(setq n(1+ n)))
    
47.                         ((= F 2)(setq Pt(insertlist Pt r n)))
    
48.                 )
    
49.         )
    
50.         (setq F(if(vlax-curve-isClosed ob)1 0)
    
51.                   b(entget(vlax-vla-object->ename ob)'("*"))
    
52.                   Pt(if(and(= F 1)(NULL(vl-string-search "Polyline" D)))(reverse(cdr(reverse Pt)))Pt)
    
53.         )
    
54.         (setq a (list(cons 0 "LWPOLYLINE")
    
55.                        (cons 8(cdr(assoc 8 b)))
    
56.                        (cons 6(if(assoc 6 b)(cdr(assoc 6 b))"Bylayer"))
    
57.                        (cons 62(if(assoc 62 b)(cdr(assoc 62 b))256))
    
58.                        (cons 370(if(assoc 370 b)(cdr(assoc 370 b))0))
    
59.                        (cons 48(if(assoc 48 b)(cdr(assoc 48 b))1))
    
60.                        (cons 100 "AcDbEntity")
    
61.                        (cons 100 "AcDbPolyline")
    
62.                        (cons 90(-(length Pt)F))
    
63.                        (cons 70 (+ 128 F))
    
64.                        (cons 43(if(assoc 43 b)(cdr(assoc 43 b))0))
    
65.                        (cons 38(caddr(vlax-curve-getstartpoint ob)))
    
66.                        (cons 39(if(assoc 39 b)(cdr(assoc 39 b))0))
    
67.                     )
    
68.                   a (append a(foreach b Pt(setq Pt(subst(subst 10(car b)b)b Pt))))
    
69.         )
    
70.     (entmake(if(assoc -3 b)(append a(list(assoc -3 b)))a))
    
71.     (prin1)
    
72. )
    
73. ;;以下为3个子函数
    
74. (defun StarEnd(ob)
    
75.    (list(append(list(vlax-curve-getstartparam ob))(vlax-curve-getstartpoint ob))
    
76.    (append(list(vlax-curve-getendparam ob))(vlax-curve-getendpoint ob)))
    
77. )
    
78. 
    
79. (defun count-f (OBJ l p q o C / F)
    
80.   (setq F (if(and (>(+(distance p q)(distance q o))l)
    
81.                 (<(abs(*(distance p q)(sin(-(angle p q)(angle p o)))))C)
    
82.                 )
    
83.                 1
    
84.                 2
    
85.               )
    
86.         )
    
87. )
    
88. 
    
89. (defun insertlist(Lst en n / l m s)
    
90.         (cond
    
91.                 ((= n 0)(setq s(append(list en)LST)))
    
92.                 ((= n(setq l(length LST)))(setq s(append LST(list en))))
    
93.                 (T(setq m 1 s(list(car LST)))
    
94.                         (while(< m n)(setq s(append s(list(nth m LST)))m(1+ m)))
    
95.                         (setq s(append s (list en)))
    
96.                         (while(< m l)(setq s(append s(list(nth m LST)))m(1+ m)))
    
97.                 )
    
98.         )
    
99.         (setq Lst s)
    
100. )

;;将各种线段转为多线段，将多线段上的弧线段也转为多线段<br /> ;;(CurveToPoline en);;参数：图元名;;示例(CurveToPoline (car(entsel "\n选择对象: ")))<br /> ;;适用"line,lwpolyline,polyline,spline,circle,arc,ellipse"<br /> (defun CurveToPoline ( en / c2 C ob a n pt D l m  q p b o r F )<br />         (vl-load-com)                        <br />         (setq C2 0.65)        ;;可理解为相似度，不要小于0.65，不能为1<br />         (setq C(/(- 1.0 C2)2))        <br />     (setq ob (vlax-ename->vla-object en)<br />               a 0 <br />               n 1 <br />               pt nil <br />               D (vla-get-objectname ob)<br />               l (vlax-curve-getdistatparam ob(vlax-curve-getendparam ob))<br />              m (cond((= D"AcDb3dPolyline")3)((= D"AcDb2dPolyline")2)((= D"AcDbPolyline")2))<br />         )<br />         (if(or(= D"AcDb3dPolyline")(= D"AcDbPolyline"))<br />           (if(vl-catch-all-error-p(vl-catch-all-apply 'vla-get-Coordinates(list ob)))(setq Pt(StarEnd ob))<br />                   (progn<br />                           (setq q(vlax-safearray->list(vlax-variant-value(vla-get-Coordinates ob))))<br />                           (repeat(/(length q)m)<br />                                     (cond<br />                                           ((= m 2)(setq p(list(nth a q)(nth(+ a 1)q))))<br />                                           ((= m 3)(setq p(list(nth a q)(nth(+ a 1)q)(nth(+ a 2)q))))<br />                                         )<br />                                         (setq a(+ a m)<br />                                                   Pt(if(or(null Pt)(>(distance p(cdr(car(reverse Pt))))0.00005))<br />                                                       (append Pt(list(append(list(vlax-curve-getparamatpoint ob p))p)))pt)<br />                                         )<br />                                 )<br />                         )<br />                 )<br />                 (setq Pt(StarEnd ob))<br />         )<br />         (while(< n(length Pt))<br />                 (setq m(1- n)<br />                           a(vlax-curve-getdistAtparam ob(car(nth m Pt)))<br />                           b(-(vlax-curve-getdistAtparam ob(car(nth n Pt)))a)<br />                 )<br />                 (if(<= b 0)(setq b(- l a)))<br />                 (setq o(vlax-curve-getparamatdist ob(+(/ b 2)a))<br />                           r(append(list o)(vlax-curve-getpointatparam ob o))<br />                           F(count-f ob (* b c2)(cdr(nth m Pt))(cdr r)(cdr(nth n Pt))c)<br />                 )<br />                 (cond<br />                         ((= F 1)(setq n(1+ n)))<br />                         ((= F 2)(setq Pt(insertlist Pt r n)))<br />                 )<br />         )<br />         (setq F(if(vlax-curve-isClosed ob)1 0)<br />                   b(entget(vlax-vla-object->ename ob)'("*"))<br />                   Pt(if(and(= F 1)(NULL(vl-string-search "Polyline" D)))(reverse(cdr(reverse Pt)))Pt)<br />         )<br />         (setq a (list(cons 0 "LWPOLYLINE")<br />                        (cons 8(cdr(assoc 8 b)))<br />                        (cons 6(if(assoc 6 b)(cdr(assoc 6 b))"Bylayer"))<br />                        (cons 62(if(assoc 62 b)(cdr(assoc 62 b))256))<br />                        (cons 370(if(assoc 370 b)(cdr(assoc 370 b))0))<br />                        (cons 48(if(assoc 48 b)(cdr(assoc 48 b))1))<br />                        (cons 100 "AcDbEntity")<br />                        (cons 100 "AcDbPolyline")<br />                        (cons 90(-(length Pt)F))<br />                        (cons 70 (+ 128 F))<br />                        (cons 43(if(assoc 43 b)(cdr(assoc 43 b))0))<br />                        (cons 38(caddr(vlax-curve-getstartpoint ob)))<br />                        (cons 39(if(assoc 39 b)(cdr(assoc 39 b))0))<br />                     )<br />                   a (append a(foreach b Pt(setq Pt(subst(subst 10(car b)b)b Pt))))<br />         )<br />     (entmake(if(assoc -3 b)(append a(list(assoc -3 b)))a))<br />     (prin1)<br /> )<br /> ;;以下为3个子函数<br /> (defun StarEnd(ob)<br />    (list(append(list(vlax-curve-getstartparam ob))(vlax-curve-getstartpoint ob))<br />    (append(list(vlax-curve-getendparam ob))(vlax-curve-getendpoint ob)))<br /> )<br /> <br /> (defun count-f (OBJ l p q o C / F)<br />   (setq F (if(and (>(+(distance p q)(distance q o))l)<br />                 (<(abs(*(distance p q)(sin(-(angle p q)(angle p o)))))C)<br />                 )<br />                 1<br />                 2<br />               )<br />         )<br /> )<br /> <br /> (defun insertlist(Lst en n / l m s)<br />         (cond<br />                 ((= n 0)(setq s(append(list en)LST)))<br />                 ((= n(setq l(length LST)))(setq s(append LST(list en))))<br />                 (T(setq m 1 s(list(car LST)))<br />                         (while(< m n)(setq s(append s(list(nth m LST)))m(1+ m)))<br />                         (setq s(append s (list en)))<br />                         (while(< m l)(setq s(append s(list(nth m LST)))m(1+ m)))<br />                 )<br />         )<br />         (setq Lst s)<br /> )

qyming

曲线拟和，整理

ivde

楼主想输出pat还是shp?

风树

ivde 发表于 2014-7-13 07:38
楼主想输出pat还是shp?

想大致试试制作填充的一种方向

xyp1964

填充图案制作 http://bbs.mjtd.com/thread-110273-1-1.html

1mi8

wzg356 发表于 2014-7-13 01:22
很抱歉，记不清原帖子的位置了
适用范围比较大
曲线、椭圆、圆、圆弧等转换的数学指标按相似度调整，如果 ...

这个程序的命令不是“CurveToPoline”吗，为什么我显示未知命令啊。

luchangliang

同楼上，显示未知命令