用XDGE几何库算多段线自交点参数值

newer

http://bbs.xdcad.net/thread-711477-1-1.html
上面 Lispboy 的帖子链接介绍了API函数 xdrx_polyline_isSelfIntersect 判断多段线是否自交，并返回自交参数表

下面介绍下，使用XDGE几何库，用晓东通用LISP函数模拟ARX解算的过程。

1. (defun c:tt ()
   
2.   (if
   
3.     (and (setq
   
4.            e (car
   
5.                (xdrx_entsel "\n拾取多段线<退出>:" '((0 . "LWPOLYLINE")))
   
6.              )
   
7.          )
   
8.          (setq ints (xdrx_entity_intersectwith e e))
   
9.          ;;自身求交点（包括自交点和顶点）
   
10.          (setq verts (xdrx_getpropertyvalue e "vertices"))
    
11.          ;;获得多段线顶点表
    
12.          (setq g (xdge::constructor e))
    
13.          ;;构建多段线几何实体模型
    
14.          (setq gl (xdge::getpropertyvalue g "getcurvelist"))
    
15.          ;;获得几何多段线各段的几何实体（包括LINE,ARC）
    
16.     )
    
17.      (progn
    
18.        (setq params nil)
    
19.        (setq
    
20.          ints (vl-remove-if
    
21.                 '(lambda (x) (XD::List:ItemAllPosition verts x 1e-3))
    
22.                 ints
    
23.               )
    
24.        )
    
25.        ;;删除掉顶点，剩下自交点
    
26.        (foreach        n ints
    
27.          (setq inx 0)
    
28.          (foreach m gl
    
29.            ;;循环处理每段，求交点对应的两个参数
    
30.            (if (setq pa (xdge::getpropertyvalue m "ison" n))
    
31.              ;;判断自交点是否在当前段上，返回参数值（ON）
    
32.              (progn
    
33.                (setq dis (xdge::getpropertyvalue
    
34.                            g
    
35.                            "localToGlobalParam"
    
36.                            pa
    
37.                            inx
    
38.                          )
    
39.                )
    
40.                ;;返回当前段参数对应全局多段线成参数值
    
41.                (setq pa (xdrx_curve_getparamatdist e dis))
    
42.                ;;dis全局参数（几何实体）对应数据库实体的长度值，由长度算出数据库实体参数值
    
43.                (setq params (cons pa params))
    
44.              )
    
45.            )
    
46.            (setq inx (1+ inx))
    
47.          )
    
48.        )
    
49.        (setq params (reverse params))
    
50.        (princ "\n计算参数值表:")
    
51.        (princ params)
    
52.        (princ "\nAPI 参数值表:")
    
53.        (princ (xdrx_polyline_isSelfIntersect e 2))
    
54.        (xdge::free g gl) ;;释放几何对象内存
    
55.      )
    
56.   )
    
57.   (princ)
    
58. )
    

(defun c:tt ()<br />   (if<br />     (and (setq<br />            e (car<br />                (xdrx_entsel "\n拾取多段线<退出>:" '((0 . "LWPOLYLINE")))<br />              )<br />          )<br />          (setq ints (xdrx_entity_intersectwith e e))<br />          ;;自身求交点（包括自交点和顶点）<br />          (setq verts (xdrx_getpropertyvalue e "vertices"))<br />          ;;获得多段线顶点表<br />          (setq g (xdge::constructor e))<br />          ;;构建多段线几何实体模型<br />          (setq gl (xdge::getpropertyvalue g "getcurvelist"))<br />          ;;获得几何多段线各段的几何实体（包括LINE,ARC）<br />     )<br />      (progn<br />        (setq params nil)<br />        (setq<br />          ints (vl-remove-if<br />                 '(lambda (x) (XD::List:ItemAllPosition verts x 1e-3))<br />                 ints<br />               )<br />        )<br />        ;;删除掉顶点，剩下自交点<br />        (foreach        n ints<br />          (setq inx 0)<br />          (foreach m gl<br />            ;;循环处理每段，求交点对应的两个参数<br />            (if (setq pa (xdge::getpropertyvalue m "ison" n))<br />              ;;判断自交点是否在当前段上，返回参数值（ON）<br />              (progn<br />                (setq dis (xdge::getpropertyvalue<br />                            g<br />                            "localToGlobalParam"<br />                            pa<br />                            inx<br />                          )<br />                )<br />                ;;返回当前段参数对应全局多段线成参数值<br />                (setq pa (xdrx_curve_getparamatdist e dis))<br />                ;;dis全局参数（几何实体）对应数据库实体的长度值，由长度算出数据库实体参数值<br />                (setq params (cons pa params))<br />              )<br />            )<br />            (setq inx (1+ inx))<br />          )<br />        )<br />        (setq params (reverse params))<br />        (princ "\n计算参数值表:")<br />        (princ params)<br />        (princ "\nAPI 参数值表:")<br />        (princ (xdrx_polyline_isSelfIntersect e 2))<br />        (xdge::free g gl) ;;释放几何对象内存<br />      )<br />   )<br />   (princ)<br /> )<br />

执行结果：
命令: TT
拾取多段线<退出>:
计算参数值表:(1.17218 3.62837 1.81808 5.18041 3.48957 6.13156)
API 参数值表:(1.17218 3.62837 1.81808 5.18041 3.48957 6.13156)


得到的结果完全相同。