多段线自交判断

sen.sam

在写一个小程序需要先判断多段线是否自交，找了一下论坛，很少相关的资料。按照某网友的提示凑了一个，但是只对闭合的完全自交的多段线（图1）有效，像图2有一个顶点落在边上的情况就不行了，另外对于虽然自交但不闭合（始点和终点重合），的情况也不支持。不知道各位有没更好的办法去判断这种情况？

2. (defun c:test(/ ptlist1 ent1 ent2 ent0 ps entlist listlen1 listlen2)
   
3. (setq ent1 (car (entsel "选择实体:")));
   
4. (command "_.copy" ent1 "" '(0 0) "")
   
5. (setq ent2 (entlast))
   
6. (setq ps (obj_int ent1 ent2))
   
7. (setq ent0  (entget ent1))
   
8. (mapcar '(lambda(X)(if (= 10 (car x))(setq ptlist1 (cons (cdr x) ptlist1))))
   
9. 
   
10. ent0)
    
11. (setq entlist(reverse ptlist1))
    
12. (setq listlen1 (length ptlist1))
    
13. (setq listlen2 (length ps))
    
14. (if (/= listlen2 listlen1 )
    
15.        (alert "相交")
    
16.       (alert "不相交")
    
17.      )  
    
18. (command "erase" ent2 "")
    
19. (princ)
    
20.   )
    
21. 
    
22. 
    
23. ;功能:返回两个对象的所有交点
    
24. ;参数: ent1、ent2 均为ename对象
    
25. (defun obj_int (ent1 ent2 / ax_ent_1 ax_ent_2 intpoints points i)
    
26.    (setq ax_ent_1 (vlax-ename->vla-object ent1)
    
27.          ax_ent_2 (vlax-ename->vla-object ent2)
    
28.    )
    
29.    (setq intpoints (vla-intersectwith ax_ent_1 ax_ent_2 acextendnone))
    
30.    (setq intpoints (vlax-variant-value intpoints))
    
31.    (setq i 0)
    
32.    (if (> (vlax-safearray-get-u-bound intpoints 1) 0)
    
33.      (repeat (/ (+ 1
    
34.                (- (vlax-safearray-get-u-bound intpoints 1)
    
35.                   (vlax-safearray-get-l-bound intpoints 1)
    
36.                )
    
37.             )
    
38.             3
    
39.          )
    
40.        (setq points (append points (list (list
    
41.          (vlax-safearray-get-element intpoints i)
    
42.          (vlax-safearray-get-element intpoints (+ i 1))
    
43.          (vlax-safearray-get-element intpoints (+ i 2))
    
44.        )))
    
45.        )
    
46.        (setq i (+ 3 i))
    
47.      )
    
48.    )
    
49.    points
    
50. )
    

(defun c:test(/ ptlist1 ent1 ent2 ent0 ps entlist listlen1 listlen2)<br /> (setq ent1 (car (entsel "选择实体:")));<br /> (command "_.copy" ent1 "" '(0 0) "")<br /> (setq ent2 (entlast))<br /> (setq ps (obj_int ent1 ent2))<br /> (setq ent0  (entget ent1))<br /> (mapcar '(lambda(X)(if (= 10 (car x))(setq ptlist1 (cons (cdr x) ptlist1))))<br /> <br /> ent0)<br /> (setq entlist(reverse ptlist1))<br /> (setq listlen1 (length ptlist1))<br /> (setq listlen2 (length ps))<br /> (if (/= listlen2 listlen1 )<br />        (alert "相交")<br />       (alert "不相交")<br />      )  <br /> (command "erase" ent2 "")<br /> (princ)<br />   )<br /> <br /> <br /> ;功能:返回两个对象的所有交点<br /> ;参数: ent1、ent2 均为ename对象<br /> (defun obj_int (ent1 ent2 / ax_ent_1 ax_ent_2 intpoints points i)<br />    (setq ax_ent_1 (vlax-ename->vla-object ent1)<br />          ax_ent_2 (vlax-ename->vla-object ent2)<br />    )<br />    (setq intpoints (vla-intersectwith ax_ent_1 ax_ent_2 acextendnone))<br />    (setq intpoints (vlax-variant-value intpoints))<br />    (setq i 0)<br />    (if (> (vlax-safearray-get-u-bound intpoints 1) 0)<br />      (repeat (/ (+ 1<br />                (- (vlax-safearray-get-u-bound intpoints 1)<br />                   (vlax-safearray-get-l-bound intpoints 1)<br />                )<br />             )<br />             3<br />          )<br />        (setq points (append points (list (list<br />          (vlax-safearray-get-element intpoints i)<br />          (vlax-safearray-get-element intpoints (+ i 1))<br />          (vlax-safearray-get-element intpoints (+ i 2))<br />        )))<br />        )<br />        (setq i (+ 3 i))<br />      )<br />    )<br />    points<br /> )<br />

jicqj

好难  先留个脚印 下次来看

vormittag

楼主的思路是比较(obj_int en en)得到的交点数目与en的顶点数目是否一致，但是由于en对应的多段线图元型式的不同，(obj_int en en)的结果并不稳定。
经测试：
1. 当en不闭合，(obj_int en en) 得到的结果不包含en的起终点，所以对于不闭合的多义线，只要自交点不是两个，都会判断为自交，自交点若是两个则反而判断不自交，有趣。
2. 当en闭合，不论是图面闭合（仅终点重合）还是标志闭合（图元闭合标志为1）或者二者兼而有之，(obj_int en en) 的结果都包含en的起终点。

因为求交点利用了vla-intersectwith（用(vla-invoke ... 'intersectwith ... ...)也一样），我想可能还是要研究一下函数返回值吧。

btw:

2. (mapcar '(lambda(X)(if (= 10 (car x))(setq ptlist1 (cons (cdr x) ptlist1)))) ent0)
   
3. (setq entlist(reverse ptlist1))
   
4. (setq listlen1 (length ptlist1))
   

(mapcar '(lambda(X)(if (= 10 (car x))(setq ptlist1 (cons (cdr x) ptlist1)))) ent0)<br /> (setq entlist(reverse ptlist1))<br /> (setq listlen1 (length ptlist1))<br />

是不是太麻烦了？

2. (setq listlen1 (cdr (assoc 90 ent0)))

(setq listlen1 (cdr (assoc 90 ent0)))

不是更简洁一些么？

sen.sam

有个想法，把多段线打散，逐条比较是否有交点。或者在多段线每两个顶点画一条直线，有凸度的地方计算出圆心半径画圆弧，再逐条判断每条线或者圆弧是否有交点。只是这样算起来太费劲了。不知道有没更好的方法？

sen.sam

自己顶一下！

caoyin

多段线还是用几何方法计算比较稳妥，效率也最高

vormittag

无意中在Express Tools的extrim.lsp中发现带有函数能够判断多段线自交-- acet-geom-self-intersect 函数。
不过这个函数的源代码及参数含义没有找到。

miraclelll

楼主, 我采用如下方法,可以解决这个问题, 且测试下来, 效率也是可以的~

解决方法:  对于多边形自交的定义是: 不相邻的线段有交点, 即为自交. 故此, 我采用对多边形的组成线段之间求交点的方式来判断.

编程思路:  
1 将多边形en 的节点组成一个表 ptlist, 首尾点重合(按先后顺序), 如: (pt1,pt2,pt3,...... ptn)
2 每两相邻点设定为 此多边形的 组成线段.
3 对ptlist中的点进行循环,从 第一个点pt1开始, 至pt(n-1) 为止
  取出其中的一点pt(i) 做为 线段L 的起点 pt1, pt(i+1)为终点pt2
   再取 pt(i+2) 为 对比判断线段M 的起点 pt3, pt(i+3) 为终点pt4
   然后 求L与M的交点, 如果有交点,则可判定 此多边形为自交
4 其中, 线段M 需要循环取, 但是注意不要和线段L 相邻

程序代码如下:

2. (defun selfint (en /
   
3.                        ptlist x n i j pt1 pt2 pt3 pt4 rtn pt_inter
   
4.                        )
   
5.   (setq rtn nil pt_inter nil)
   
6.   ;;; 点坐标串ptlist;
   
7.   (setq ptlist (mapcar 'cdr
   
8.                        (vl-remove-if-not
   
9.                          '(lambda (x) (= 10 (car x)))
   
10.                          (entget en))
    
11.                        )
    
12.              )
    
13.   ;;  首尾坐标不重合的,需要增加;;;
    
14.   (if (not (equal
    
15.         (nth 0 ptlist)
    
16.         (nth (1- (length ptlist)) ptlist)
    
17.         ))
    
18.     (setq ptlist (append ptlist (list (nth 0 ptlist))))
    
19.     )
    
20.   (setq n (length ptlist))
    
21.   (setq i 0
    
22.         j (+ i 2)
    
23.         )
    
24.   ;;; 开始判断;
    
25.   (while (and (<= i (- n 4)) (= rtn nil))
    
26.     ;;; 线段L 的起始点;;
    
27.     (setq pt1 (nth i ptlist)
    
28.           pt2 (nth (1+ i) ptlist)
    
29.           )
    
30.     ;;; 下面是循环取 线段M的起始点;
    
31.     (while (and (<= j (- n 2)) (= rtn nil))
    
32.       (setq pt3 (nth j ptlist)
    
33.             pt4 (nth (1+ j) ptlist)
    
34.             )
    
35.       (if (not (and (= i 0) (= j (- n 2)))) ;; 对于第一个线段L,与最后一个线段属于相邻,故需要排除;
    
36.         (if (setq pt_inter (inters pt1 pt2 pt3 pt4));;; 判断L与M 是否相交;
    
37.           (setq rtn T)
    
38.           )
    
39.         )
    
40.       (setq j (1+ j))
    
41.       )
    
42.     (setq i (1+ i)
    
43.           j (+ i 2)
    
44.           )
    
45.     )
    
46.   (list rtn pt_inter)
    
47.   )

(defun selfint (en /<br />                        ptlist x n i j pt1 pt2 pt3 pt4 rtn pt_inter<br />                        )<br />   (setq rtn nil pt_inter nil)<br />   ;;; 点坐标串ptlist;<br />   (setq ptlist (mapcar 'cdr<br />                        (vl-remove-if-not<br />                          '(lambda (x) (= 10 (car x)))<br />                          (entget en))<br />                        )<br />              )<br />   ;;  首尾坐标不重合的,需要增加;;;<br />   (if (not (equal<br />         (nth 0 ptlist)<br />         (nth (1- (length ptlist)) ptlist)<br />         ))<br />     (setq ptlist (append ptlist (list (nth 0 ptlist))))<br />     )<br />   (setq n (length ptlist))<br />   (setq i 0<br />         j (+ i 2)<br />         )<br />   ;;; 开始判断;<br />   (while (and (<= i (- n 4)) (= rtn nil))<br />     ;;; 线段L 的起始点;;<br />     (setq pt1 (nth i ptlist)<br />           pt2 (nth (1+ i) ptlist)<br />           )<br />     ;;; 下面是循环取 线段M的起始点;<br />     (while (and (<= j (- n 2)) (= rtn nil))<br />       (setq pt3 (nth j ptlist)<br />             pt4 (nth (1+ j) ptlist)<br />             )<br />       (if (not (and (= i 0) (= j (- n 2)))) ;; 对于第一个线段L,与最后一个线段属于相邻,故需要排除;<br />         (if (setq pt_inter (inters pt1 pt2 pt3 pt4));;; 判断L与M 是否相交;<br />           (setq rtn T)<br />           )<br />         )<br />       (setq j (1+ j))<br />       )<br />     (setq i (1+ i)<br />           j (+ i 2)<br />           )<br />     )<br />   (list rtn pt_inter)<br />   )

llsheng_73

如果多线段的任一子线段与多线段的交点数大于2那么它一定是自相交