[分享]：spline->pl的高效代码

cy956

是否还要加入曲率变化判断？
或者判断p1p2曲线长度与直线长度的关系？
这个问题再发展下去，是不是也要来个“优选法”模型？

aeo

最初由 cy956 发布
[B]是否还要加入曲率变化判断？
或者判断p1p2曲线长度与直线长度的关系？
这个问题再发展下去，是不是也要来个“优选法”模型？ [/B]

曲率和半径是没关系的，所以大圆和小圆等分下去是一样的(份数)
这就失去了意义.

wkai

这段代码的特点是能够根据曲线的曲率分配顶点数，曲率变化大的地方顶点的密度大。

2. 
   
3.   [FONT=courier new]
   
4. ;__2004.01.12.09.43_____________________________________
   
5. ;__样条曲线--〉复义线___BY__WKAI__晓东CAD论坛___________
   
6. (defun c:tt (/ pts min_distacne vlo st_par ed_par ent n max_angle_rad max_an)
   
7.   (setq        vlo              (vlax-ename->vla-object (setq ent (car (entsel "\n选择spline："))))
   
8.         st_par              (vlax-curve-getstartparam vlo)
   
9.         ed_par              (vlax-curve-getendparam vlo)
   
10.         pts              (list (vlax-curve-getstartpoint vlo))
    
11.         max_angle     (if (setq        max_an (getreal        (strcat        "\n许可角度误差<"
    
12.                                                         (if max_angle
    
13.                                                           (rtos max_angle 2 2)
    
14.                                                           (rtos (setq max_angle 5) 2 2)
    
15.                                                         )
    
16.                                                         ">:"
    
17.                                                 )
    
18.                                        )
    
19.                           )
    
20.                         max_an
    
21.                         max_angle
    
22.                       )
    
23.         max_angle_rad (* pi (/ max_angle 180.0)) ;_最大许可误差角度弧度
    
24.   )
    
25.   (setq zhf_time_dot nil)
    
26.   (z_timer)
    
27.   (t1 vlo st_par ed_par)
    
28.   (command "pline")
    
29.   (foreach n pts (command "non" n))
    
30.   (command)
    
31.   (princ (strcat "\n创建一条顶点数为" (itoa (length pts)) "的复义线。"))
    
32.   (princ (z_timer))
    
33.   (princ)
    
34. )
    
35. ;________________________________________________________
    
36. ;________________________________________________________
    
37. (defun t1 (vl_obj start_par end_par / dis_m_m ed_p ed_pt mid_p mid_pt mid_pt_vlo st_p st_pt)
    
38.   (setq        st_p           start_par
    
39.         ed_p           end_par
    
40.         st_pt           (vlax-curve-getpointatparam vl_obj st_p)
    
41.         ed_pt           (if (vlax-curve-getpointatparam vl_obj ed_p)
    
42.                      (vlax-curve-getpointatparam vl_obj ed_p)
    
43.                      (vlax-curve-getendpoint vl_obj)
    
44.                    )
    
45.         mid_p           (/ (+ st_p ed_p) 2.0)
    
46.         mid_pt           (mapcar '(lambda (x y) (/ (+ x y) 2.0)) st_pt ed_pt)
    
47.         mid_pt_vlo (vlax-curve-getpointatparam vl_obj mid_p)
    
48.         dis_m_m           (distance mid_pt_vlo mid_pt)
    
49.   )
    
50.   (if (or (> (t3 ed_pt st_pt mid_pt_vlo t) max_angle_rad)
    
51.           (and (< (t3 ed_pt st_pt mid_pt_vlo t) max_angle_rad)
    
52.                (t2 vl_obj st_p ed_p st_pt ed_pt)
    
53.           )
    
54.       )
    
55.     (progn
    
56.       (t1 vl_obj st_p mid_p)
    
57.       (t1 vl_obj mid_p ed_p)
    
58.     )
    
59.     (setq pts (append pts (list mid_pt_vlo) (list ed_pt)))
    
60.   )
    
61. )
    
62. ;_____________________________________________________________
    
63. (defun t2 (vl st_par ed_par s_pt e_pt / run delta_p n lst pt)
    
64.   (setq        run        t
    
65.         delta_p        (- ed_par st_par)
    
66.         lst        '(0.0001 0.01 0.25) ;_测试采样表
    
67.         n        -1
    
68.   )
    
69.   (while (and run (< (setq n (1+ n)) (length lst)))
    
70.     (setq pt (vlax-curve-getpointatparam vl (+ st_par (* delta_p (nth n lst)))))
    
71.     (if        (> (t3 e_pt s_pt pt t) max_angle_rad)
    
72.       (setq run nil)
    
73.     )
    
74.     (if        (> (t3 s_pt e_pt pt t) max_angle_rad)
    
75.       (setq run nil)
    
76.     )
    
77.     (setq pt (vlax-curve-getpointatparam vl (- ed_par (* delta_p (nth n lst)))))
    
78.     (if        (> (t3 e_pt s_pt pt t) max_angle_rad)
    
79.       (setq run nil)
    
80.     )
    
81.     (if        (> (t3 s_pt e_pt pt t) max_angle_rad)
    
82.       (setq run nil)
    
83.     )
    
84.   )
    
85.   (not run)
    
86. )
    
87. ;_____________________________________________________________
    
88. (defun t3 (p1 p2 p3 rad)
    
89.   (if rad
    
90.     (abs (- (angle p2 p3) (angle p2 p1)))
    
91.     (* 180 (/ (abs (- (angle p2 p3) (angle p2 p1))) pi))
    
92.   )
    
93. )
    
94. ;________________________________________________________
    
95. ;________________________________________________________
    
96. (princ "\n样条曲线--〉复义线___BY__WKAI__晓东CAD论坛")  [/FONT]
    

<br />   [FONT=courier new]<br /> ;__2004.01.12.09.43_____________________________________<br /> ;__样条曲线--〉复义线___BY__WKAI__晓东CAD论坛___________<br /> (defun c:tt (/ pts min_distacne vlo st_par ed_par ent n max_angle_rad max_an)<br />   (setq        vlo              (vlax-ename->vla-object (setq ent (car (entsel "\n选择spline："))))<br />         st_par              (vlax-curve-getstartparam vlo)<br />         ed_par              (vlax-curve-getendparam vlo)<br />         pts              (list (vlax-curve-getstartpoint vlo))<br />         max_angle     (if (setq        max_an (getreal        (strcat        "\n许可角度误差<"<br />                                                         (if max_angle<br />                                                           (rtos max_angle 2 2)<br />                                                           (rtos (setq max_angle 5) 2 2)<br />                                                         )<br />                                                         ">:"<br />                                                 )<br />                                        )<br />                           )<br />                         max_an<br />                         max_angle<br />                       )<br />         max_angle_rad (* pi (/ max_angle 180.0)) ;_最大许可误差角度弧度<br />   )<br />   (setq zhf_time_dot nil)<br />   (z_timer)<br />   (t1 vlo st_par ed_par)<br />   (command "pline")<br />   (foreach n pts (command "non" n))<br />   (command)<br />   (princ (strcat "\n创建一条顶点数为" (itoa (length pts)) "的复义线。"))<br />   (princ (z_timer))<br />   (princ)<br /> )<br /> ;________________________________________________________<br /> ;________________________________________________________<br /> (defun t1 (vl_obj start_par end_par / dis_m_m ed_p ed_pt mid_p mid_pt mid_pt_vlo st_p st_pt)<br />   (setq        st_p           start_par<br />         ed_p           end_par<br />         st_pt           (vlax-curve-getpointatparam vl_obj st_p)<br />         ed_pt           (if (vlax-curve-getpointatparam vl_obj ed_p)<br />                      (vlax-curve-getpointatparam vl_obj ed_p)<br />                      (vlax-curve-getendpoint vl_obj)<br />                    )<br />         mid_p           (/ (+ st_p ed_p) 2.0)<br />         mid_pt           (mapcar '(lambda (x y) (/ (+ x y) 2.0)) st_pt ed_pt)<br />         mid_pt_vlo (vlax-curve-getpointatparam vl_obj mid_p)<br />         dis_m_m           (distance mid_pt_vlo mid_pt)<br />   )<br />   (if (or (> (t3 ed_pt st_pt mid_pt_vlo t) max_angle_rad)<br />           (and (< (t3 ed_pt st_pt mid_pt_vlo t) max_angle_rad)<br />                (t2 vl_obj st_p ed_p st_pt ed_pt)<br />           )<br />       )<br />     (progn<br />       (t1 vl_obj st_p mid_p)<br />       (t1 vl_obj mid_p ed_p)<br />     )<br />     (setq pts (append pts (list mid_pt_vlo) (list ed_pt)))<br />   )<br /> )<br /> ;_____________________________________________________________<br /> (defun t2 (vl st_par ed_par s_pt e_pt / run delta_p n lst pt)<br />   (setq        run        t<br />         delta_p        (- ed_par st_par)<br />         lst        '(0.0001 0.01 0.25) ;_测试采样表<br />         n        -1<br />   )<br />   (while (and run (< (setq n (1+ n)) (length lst)))<br />     (setq pt (vlax-curve-getpointatparam vl (+ st_par (* delta_p (nth n lst)))))<br />     (if        (> (t3 e_pt s_pt pt t) max_angle_rad)<br />       (setq run nil)<br />     )<br />     (if        (> (t3 s_pt e_pt pt t) max_angle_rad)<br />       (setq run nil)<br />     )<br />     (setq pt (vlax-curve-getpointatparam vl (- ed_par (* delta_p (nth n lst)))))<br />     (if        (> (t3 e_pt s_pt pt t) max_angle_rad)<br />       (setq run nil)<br />     )<br />     (if        (> (t3 s_pt e_pt pt t) max_angle_rad)<br />       (setq run nil)<br />     )<br />   )<br />   (not run)<br /> )<br /> ;_____________________________________________________________<br /> (defun t3 (p1 p2 p3 rad)<br />   (if rad<br />     (abs (- (angle p2 p3) (angle p2 p1)))<br />     (* 180 (/ (abs (- (angle p2 p3) (angle p2 p1))) pi))<br />   )<br /> )<br /> ;________________________________________________________<br /> ;________________________________________________________<br /> (princ "\n样条曲线--〉复义线___BY__WKAI__晓东CAD论坛")  [/FONT]<br />

陌生人

想问问
1.转换的意义？什么情况下需要转换
2。转换优化的方法，距离逼近好，还是角度逼近好，曲率控制？其它的方法？
3.角度控制和曲率控制等同么？

另，没仔细看完程序，问wkai，'(0.0001 0.01 0.25) ;_测试采样表
是采样时候的分级取点用的么？也就是说，你是用逼近角度的方法？

aeo

最初由 wkai 发布
[B]这段代码的特点是能够根据曲线的曲率分配顶点数，曲率变化大的地方顶点的密度大。
[CODE]
  
;__2004.01.12.09.43_____________________________________
;__样条曲线--〉复义线___BY__WKAI... [/B]

spline谈不上曲率,因为你测试的点天知道在哪一个弯上面，也许你的两点已经经过了好几个弯了.

我上面说了:曲率,对大的弧合小的弧是一样的(和半径无关)
所以长1000的线和1的线,等分下去是一样的.

其实就是个模拟，大差不差就行了(不可能讲"正确")
也就是要用尽量少的点，来模拟出适合你要求精度的线来.
精度可以设，关键是取点.

其实规定：一天要三个苹果的营养，你要是吃三个苹果吧，可能营养不够，你要吃到那么多营养吧，又不知道是几个.因为规定就是个模糊概念.
这里的模拟本身就是个模糊的事情,挑3个中等大小的吃了就可以了,当然吃3个大的更好（有时就浪费),当你在考虑苹果的产地时，就复杂了.



我的方法是:吃一个，吃到小的再吃一个。（ 好象精度不够:)  )
你的方法：拨一半，看多了，再拨...
当然还可以:切碎，用天平:) 我看没必要了吧.

njsabc

对于以上同志提出的解答，本人深表谢意。该程序正式我想要的代码。谢谢

陌生人

为什么没一个帖子模拟的时候用bulge？不好用么？还是有什么缺陷？