[讨论]：如何快速判断那些线相交

dubing

大家好！我有这样一个问题

1.我们首先在CAD中创建3个图层。Layer1,Layer2,Layer3
2. 任意在这三个图层中画 Line线(描画时故意让这些线彼此相交
但是，故意使一些交点不是所有线集合的首尾点

3.采用什么算法可以非常快速的判断出这些交点，并在交点处画一个圆

说明一下，如果您是从事GIS工作的人员，会理解这个检查程序目的是检查应该进行溶点的地方没有溶点的错误。

我有一个办法但太慢了！

我说说我的想法，遍历所有要处理的线实体，以首尾点做一个类似Fence的选择集（选择所有经过这两个点的实体）循环判断每一个返回的结果选择集中的线与当前线段求交点（使用CAD求交函数）如果交点不是结果选择集中线的首尾点，就报错，画圆。

但是如果我们需要处理的线是，没有规律的线，实体数量多达4万 - 5万 个Line实体， 如何才能更快的判断未溶断的点？



         希望大家踊跃发言   谢谢！

eachy

最初由 dubing 发布
[B]大家好！我有这样一个问题

1.我们首先在CAD中创建3个图层。Layer1,Layer2,Layer3
2. 任意在这三个图层中画 Line线(描画时故意让这些线彼此相交
但是，故意使一些交点不是所有线集合的首尾点

3.采用什么算法?.. [/B]

1 求选择集的交点集（去除重复）
2 求线的起始点集（去除重复）
3 两个点集运算

再快就只能求助ARX了。

dubing

请教 如何求选择集的交点集？

eachy

最初由 dubing 发布
[B]请教 如何求选择集的交点集？ [/B]

1 xdapi 中的 xdrx_getinterss 函数 使用方法 (xdrx_getinterss ss) ss---选择集

2 VLA 方法，论坛有现成的代码，搜索LUCAS，但是那个程序没有过滤重复点。

自由的鱼

想不到可以这样实现，又学了一招，谢谢楼上的几位！

dubing

不太理解 eachy 所说的 原理
1 xdapi 中的 xdrx_getinterss 函数 使用方法 (xdrx_getinterss ss) ss---选择集
2 VLA 方法，论坛有现成的代码，搜索LUCAS，但是那个程序没有过滤重复点。

如果使用选择集，我想基本的方法是这样：
  遍历，选择集中的实体 取出其中的一个Line线，接下来就是与谁求交点的问题
  .以这条Line线为范围选择所有与其相交的实体(参数CP)再求该线与所有结果集中的线的交点，这样做的结果是频繁调用选择集，如果当前数据库中的实体为万级，每做一次都需要时间，累级效应将是几个小时，而AutoCAD Map 的CleanUp功能中就有类似的功能（交点打断）它用几分钟甚至几秒钟就能够判断出来这是为什么？

   希望大家开动脑筋提出自己的好算法，算法是程序的灵魂。

                                                        沈阳 杜冰
                                                        dubing771216@hotmail.com

eachy

最初由 dubing 发布
[B]不太理解 eachy 所说的 原理
1 xdapi 中的 xdrx_getinterss 函数 使用方法 (xdrx_getinterss ss) ss---选择集
2 VLA 方法，论坛有现成的代码，搜索LUCAS，但是那个程序没有过滤重复点。

如果使用选择集，我想?.. [/B]

2. 
   
3. [PHP]
   
4. xdrx_getinters
   
5. 
   
6. 功    能：1.求两个AcDbCurve(曲线)实体的交点.
   
7.           2.求一个AcDbCurve(曲线)实体和一个选择集中所有AcDbCurve(曲线)的交点。
   
8.           3.求一个选择集中所有AcDbCurve(曲线)实体的交点.
   
9.           4.求一个选择集SS1中的所有AcDbCurve实体和另个选择集SS2所有AcDbCurve实体的交点。
   
10. 
    
11. 调用格式： 1. (xdrx_getinters <实体图元1> <实体图元2> [扩展方式])
    
12.            2. (xdrx_getinters <实体图元名> <选择集> [扩展方式])
    
13.            3. (xdrx_getinters <选择集> [扩展方式])
    
14.            4. (xdrx_getinters <选择集1> <选择集2> [扩展方式])
    
15.            5. (xdrx_getinters (<pt1>...<ptn>) <曲线实体名> [扩展方式])
    
16.            6. (xdrx_getinters (<pt1>...<ptn>) <选择集> [扩展方式])
    
17.            7. (xdrx_getinters (<pt1>...<ptn>) (<pt1>...<ptn>) [扩展方式])
    
18. 
    
19. 说    明：1.程序能求AcDbCurve(包括：AcDbLine、AcDbArc、AcDbCircle、AcDbEllipse、AcDbSpline、
    
20.             AcDbPolyline(lwpolyline)交点。
    
21.           2.参数[扩展方式]：
    
22.           对于1种情况：
    
23.                 0：不延伸实体
    
24.                 1：延伸第一个实体
    
25.                 2：延伸第二个实体
    
26.                 3: 两个都延伸
    
27.           对于第二、三、四种情况(选择集):
    
28.                 0: 不延伸实体
    
29.                 1: 两个都延伸
    
30. 
    
31. 示    例： 1. 求天正建筑软件轴线层上，所有轴线的交点
    
32.                 (if (setq ss (ssget "x" (list (cons 8 (xdrx_getlyrname "轴线")) '((0 . "line")))))
    
33.                     (progn
    
34.                         (setq intl (xdrx_getinters ss))
    
35.                     )
    
36.                 )
    
37. 
    
38.            2. 点取一根轴线，然后求出与这个轴线相交的所有轴线交点。
    
39.                 (setq lyr (xdrx_getlyrname "轴线"))
    
40.                 (if (setq en (xdrx_entsel "\n请点取一个轴线<退出>:" (list (cons 8 lyr) '(0 . "line"))))
    
41.                     (progn
    
42.                         (xdrx_setentodb)
    
43.                         (setq p1 (xdrx_getentdxf 10)
    
44.                               p2 (xdrx_getentdxf 11)
    
45.                               ss (ssget "f" p1 p2 (list (cons 8 lyr) '(0 . "line")))
    
46.                         )
    
47.                         (if (setq intl (xdrx_getinters (car en) ss))
    
48.                             (progn
    
49.                                .....;;intl中为交点表
    
50.                             )
    
51.                         )
    
52.                     )
    
53.                 )
    
54. [/PHP]

<br /> [PHP]<br /> xdrx_getinters<br /> <br /> 功    能：1.求两个AcDbCurve(曲线)实体的交点.<br />           2.求一个AcDbCurve(曲线)实体和一个选择集中所有AcDbCurve(曲线)的交点。<br />           3.求一个选择集中所有AcDbCurve(曲线)实体的交点.<br />           4.求一个选择集SS1中的所有AcDbCurve实体和另个选择集SS2所有AcDbCurve实体的交点。<br /> <br /> 调用格式： 1. (xdrx_getinters <实体图元1> <实体图元2> [扩展方式])<br />            2. (xdrx_getinters <实体图元名> <选择集> [扩展方式])<br />            3. (xdrx_getinters <选择集> [扩展方式])<br />            4. (xdrx_getinters <选择集1> <选择集2> [扩展方式])<br />            5. (xdrx_getinters (<pt1>...<ptn>) <曲线实体名> [扩展方式])<br />            6. (xdrx_getinters (<pt1>...<ptn>) <选择集> [扩展方式])<br />            7. (xdrx_getinters (<pt1>...<ptn>) (<pt1>...<ptn>) [扩展方式])<br /> <br /> 说    明：1.程序能求AcDbCurve(包括：AcDbLine、AcDbArc、AcDbCircle、AcDbEllipse、AcDbSpline、<br />             AcDbPolyline(lwpolyline)交点。<br />           2.参数[扩展方式]：<br />           对于1种情况：<br />                 0：不延伸实体<br />                 1：延伸第一个实体<br />                 2：延伸第二个实体<br />                 3: 两个都延伸<br />           对于第二、三、四种情况(选择集):<br />                 0: 不延伸实体<br />                 1: 两个都延伸<br /> <br /> 示    例： 1. 求天正建筑软件轴线层上，所有轴线的交点<br />                 (if (setq ss (ssget "x" (list (cons 8 (xdrx_getlyrname "轴线")) '((0 . "line")))))<br />                     (progn<br />                         (setq intl (xdrx_getinters ss))<br />                     )<br />                 )<br /> <br />            2. 点取一根轴线，然后求出与这个轴线相交的所有轴线交点。<br />                 (setq lyr (xdrx_getlyrname "轴线"))<br />                 (if (setq en (xdrx_entsel "\n请点取一个轴线<退出>:" (list (cons 8 lyr) '(0 . "line"))))<br />                     (progn<br />                         (xdrx_setentodb)<br />                         (setq p1 (xdrx_getentdxf 10)<br />                               p2 (xdrx_getentdxf 11)<br />                               ss (ssget "f" p1 p2 (list (cons 8 lyr) '(0 . "line")))<br />                         )<br />                         (if (setq intl (xdrx_getinters (car en) ss))<br />                             (progn<br />                                .....;;intl中为交点表<br />                             )<br />                         )<br />                     )<br />                 ) <br /> [/PHP]

2 LUCAS

2. 
   
3. [PHP]
   
4. ;;求所有线的交点
   
5. (defun C:TT (/ SS M E1 N E2 PTS PL PT1)
   
6.   (vl-load-com)
   
7.   (setq        SS (ssget '((0 . "*LINE,ARC,CIRCLE,ELLIPSE")))
   
8.         M  0
   
9.   )
   
10.   (while (< M (1- (sslength SS)))
    
11.     (setq E1 (vlax-ename->vla-object (ssname SS M))
    
12.           N  (1+ M)
    
13.     )
    
14.     (while (< N (sslength SS))
    
15.       (setq E2        (vlax-ename->vla-object (ssname SS N))
    
16.             PTS        (vlax-variant-value (vla-intersectwith E1 E2 0))
    
17.       )
    
18.       (if (> (vlax-safearray-get-u-bound PTS 1) 0)
    
19.         (progn
    
20.           (setq PTS (vlax-safearray->list PTS))
    
21.           (while (> (length PTS) 0)
    
22.             (setq PT1 (list (car PTS) (cadr PTS) (caddr PTS)))
    
23.             (if        (not (member PT1 PL))
    
24.               (setq PL (cons PT1 PL))
    
25.               (setq PTS (cdddr PTS))
    
26.             )
    
27.           )
    
28.         )
    
29.       )
    
30.       (setq N (1+ N))
    
31.     )
    
32.     (setq M (1+ M))
    
33.   )
    
34.   PL
    
35. )[/PHP]

<br /> [PHP]<br /> ;;求所有线的交点 <br /> (defun C:TT (/ SS M E1 N E2 PTS PL PT1)<br />   (vl-load-com)<br />   (setq        SS (ssget '((0 . "*LINE,ARC,CIRCLE,ELLIPSE")))<br />         M  0<br />   )<br />   (while (< M (1- (sslength SS)))<br />     (setq E1 (vlax-ename->vla-object (ssname SS M))<br />           N  (1+ M)<br />     )<br />     (while (< N (sslength SS))<br />       (setq E2        (vlax-ename->vla-object (ssname SS N))<br />             PTS        (vlax-variant-value (vla-intersectwith E1 E2 0))<br />       )<br />       (if (> (vlax-safearray-get-u-bound PTS 1) 0)<br />         (progn<br />           (setq PTS (vlax-safearray->list PTS))<br />           (while (> (length PTS) 0)<br />             (setq PT1 (list (car PTS) (cadr PTS) (caddr PTS)))<br />             (if        (not (member PT1 PL))<br />               (setq PL (cons PT1 PL))<br />               (setq PTS (cdddr PTS))<br />             )<br />           )<br />         )<br />       )<br />       (setq N (1+ N))<br />     )<br />     (setq M (1+ M))<br />   )<br />   PL<br /> )[/PHP]