多段线优化去除非关键点，求优化

springwillow

编了一个多段线去除非关键点的小东西，但是感觉算法冗余，求优化。

2. ;;;-----------------------------------------------------------;;
   
3. ;;; 两向量相减  subtraction                                   ;;
   
4. ;;; Input: v1,v2 -vectors in R^n                              ;;
   
5. ;;; OutPut: A vector                                          ;;
   
6. ;;;-----------------------------------------------------------;;
   
7. (defun MAT:v-v (v1 v2)
   
8.   (mapcar '- v1 v2)
   
9. )
   
10. ;;;-----------------------------------------------------------;;
    
11. ;;; 两个2d向量的叉积的数值                                    ;;
    
12. ;;; 输入: 两个点（或者两个向量）                              ;;
    
13. ;;; 输出: 一个数值.如果为正则是逆时针,两向量形成的平面法线向量;;
    
14. ;;;       向上，为负则是顺时针，为零则两向量共线或平行。      ;;
    
15. ;;;       这个数值也为原点,P1,P2三点面积的两倍。              ;;
    
16. ;;;-----------------------------------------------------------;;
    
17. (defun MAT:Det2V (v1 v2)
    
18.   (- (* (car v1) (cadr v2)) (* (car v2) (cadr v1)))
    
19. )
    
20. 
    
21. ;判断多段线折向
    
22. (defun plzhexiang (p0 p1 p2)
    
23.   (MAT:Det2V(MAT:v-v p2 p0)(MAT:v-v p1 p0))
    
24.   )
    
25. (defun c:jfpl (/ ss pts sslen ename ent ent1 pt1 pt2 pt3 i)
    
26.   (if (setq ss (ssget '((0 . "LWPOLYLINE"))))
    
27.     (progn
    
28.       (repeat (setq sslen (sslength ss))
    
29.   (setq pts '())
    
30.   (setq ename (ssname ss (setq sslen (1- sslen))))
    
31.   (setq ent (entget ename))
    
32.   (mapcar '(lambda (x)(if (= (car x) 10)(setq pts (cons (cdr x) pts)))) ent)
    
33.   (cond ((and (= (cdr (assoc 70 ent)) 0)(equal (car pts)(last pts)))
    
34.       (entmod (append (vl-remove (cons 10 (car pts))(subst (cons 70 1) (assoc 70 ent) ent)) (list (cons 10 (car pts)))))
    
35.       (entupd ename)
    
36.       (setq pts (append pts (list (cadr pts))))
    
37.       )
    
38.         ((= (cdr (assoc 70 ent)) 1) (setq pts (append pts (list (car pts)(cadr pts)))))
    
39.     )
    
40.   (setq  ent1 (entget ename))
    
41.   (setq i 0)
    
42.   (if (> (length pts) 2)
    
43.     (progn
    
44.       (repeat (-(length pts)2)
    
45.         (setq pt1 (nth i pts) pt2 (nth (1+ i) pts) pt3 (nth (+ i 2) pts))
    
46.         (if (equal (plzhexiang pt1 pt2 pt3) 0 1e-6)
    
47.     (setq ent1 (vl-remove (cons 10 pt2) ent1)))
    
48.         (setq i (1+ i))
    
49.         )
    
50.       (entmod ent1)
    
51.       )
    
52.     )
    
53.   )
    
54.       )
    
55.     )
    
56.   )
    
57.       

;;;-----------------------------------------------------------;;<br /> ;;; 两向量相减  subtraction                                   ;;<br /> ;;; Input: v1,v2 -vectors in R^n                              ;;<br /> ;;; OutPut: A vector                                          ;;<br /> ;;;-----------------------------------------------------------;;<br /> (defun MAT:v-v (v1 v2)<br />   (mapcar '- v1 v2)<br /> )<br /> ;;;-----------------------------------------------------------;;<br /> ;;; 两个2d向量的叉积的数值                                    ;;<br /> ;;; 输入: 两个点（或者两个向量）                              ;;<br /> ;;; 输出: 一个数值.如果为正则是逆时针,两向量形成的平面法线向量;;<br /> ;;;       向上，为负则是顺时针，为零则两向量共线或平行。      ;;<br /> ;;;       这个数值也为原点,P1,P2三点面积的两倍。              ;;<br /> ;;;-----------------------------------------------------------;;<br /> (defun MAT:Det2V (v1 v2)<br />   (- (* (car v1) (cadr v2)) (* (car v2) (cadr v1)))<br /> )<br /> <br /> ;判断多段线折向<br /> (defun plzhexiang (p0 p1 p2)<br />   (MAT:Det2V(MAT:v-v p2 p0)(MAT:v-v p1 p0))<br />   )<br /> (defun c:jfpl (/ ss pts sslen ename ent ent1 pt1 pt2 pt3 i)<br />   (if (setq ss (ssget '((0 . "LWPOLYLINE"))))<br />     (progn<br />       (repeat (setq sslen (sslength ss))<br />   (setq pts '())<br />   (setq ename (ssname ss (setq sslen (1- sslen))))<br />   (setq ent (entget ename))<br />   (mapcar '(lambda (x)(if (= (car x) 10)(setq pts (cons (cdr x) pts)))) ent)<br />   (cond ((and (= (cdr (assoc 70 ent)) 0)(equal (car pts)(last pts)))<br />       (entmod (append (vl-remove (cons 10 (car pts))(subst (cons 70 1) (assoc 70 ent) ent)) (list (cons 10 (car pts)))))<br />       (entupd ename)<br />       (setq pts (append pts (list (cadr pts))))<br />       )<br />         ((= (cdr (assoc 70 ent)) 1) (setq pts (append pts (list (car pts)(cadr pts)))))<br />     )<br />   (setq  ent1 (entget ename))<br />   (setq i 0)<br />   (if (> (length pts) 2)<br />     (progn<br />       (repeat (-(length pts)2)<br />         (setq pt1 (nth i pts) pt2 (nth (1+ i) pts) pt3 (nth (+ i 2) pts))<br />         (if (equal (plzhexiang pt1 pt2 pt3) 0 1e-6)<br />     (setq ent1 (vl-remove (cons 10 pt2) ent1)))<br />         (setq i (1+ i))<br />         )<br />       (entmod ent1)<br />       )<br />     )<br />   )<br />       )<br />     )<br />   )<br />       

springwillow

fan_zh 发表于 2015-9-1 09:32
我发过一个多段线的优化程序，可以看我的贴

你的不是源码，看不到算法

sunjun2009

有没有把曲线变圆滑一些的啊

myjping

overkill不是有定功能吗

edata

非关键点是什么意思，同一直线上的点都去掉，貌似你这个没考虑凸度吧，如果不考虑凸度的话，我觉得求交点就可以了。

springwillow

edata 发表于 2015-9-1 00:51
非关键点是什么意思，同一直线上的点都去掉，貌似你这个没考虑凸度吧，如果不考虑凸度的话，我觉得求交点就 ...

基本是这个意思，非关键点就是指非顶点上的点，我这个凸度没影响吧，只有同直线的才有影响，对同直线的线宽会有影响。

fan_zh

我发过一个多段线的优化程序，可以看我的贴

llsheng_73

1去除直线段上中间的点：依次取1，2，3点，如果把2投影到13，距离为0，并且它在13之内，去掉2
2.折返点同上，如果2在13之外，去掉2
3.整合凸度，连续两点凸度值不为0且相等，只要第二点不是最后一点将其去掉？

最后，1，2两步是否能否综合。。。

自贡黄明儒

我 写了个多段线去重，可能跟你的差不多，出没有考虑线宽，凸度的的问题。其实要考虑也不复杂。见你用点积，应该是不错的程序了。

springwillow

llsheng_73 发表于 2015-9-1 10:33
1去除直线段上中间的点：依次取1，2，3点，如果把2投影到13，距离为0，并且它在13之内，去掉2
2.折返点同上 ...

感谢回贴，凸度和线宽确实没有考虑，只要是在一条直线上的多余点全部去除了，不在同一直线上的凸度和线宽无影响。

springwillow

自贡黄明儒 发表于 2015-9-1 10:52
我 写了个多段线去重，可能跟你的差不多，出没有考虑线宽，凸度的的问题。其实要考虑也不复杂。见你用点积， ...

黄大师，我觉得主要是在更新对象表的时候复杂了，大师有没有好的方法?

llsheng_73

springwillow 发表于 2015-9-1 12:21
黄大师，我觉得主要是在更新对象表的时候复杂了，大师有没有好的方法?

去除一个顶点的时候，需要连续去掉4个组。。。。相信这个不难
至于总点数，好象可以不用管它，基本上如果把数据分成3部分，10组之前的，第一个10组到最后一个42组之间的，最后的
再把中间部分4个一组，这样随便想去除哪个点都特别方便了，最后append起来后进行entmod就了事了
可以参考下
http://bbs.mjtd.com/forum.php?mo ... &fromuid=202795

另外凸度也是必须要考虑的，不然可能把连续几个曲线段弄成一个直线段导致彻底变形。。。。