[分享]：通过算法精确测试点与多边形的位置关系

wkai · 发表于 2003-12-10 16:24:29

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

200312101551:
以测试点为起点作超过多边形外框的线段A。
除了测线段A与多边形的各边的交点数（M）,
还需要测落在在线段A上的多边形的端点数（N）,
以及重合在线段A的多边形边线的条数（O）
如果（REM （- （+ M N） O）2）为1则在多边形内部，否则就在外部。

200312101730:
应 zhynt 的要求改为纯 LISP 版，请大家指正。

[FONT=courier new]
;_精确测试点与多边形的位置关系__纯lisp版__BY WKAI __晓东CAD论坛
;_晓东CAD论坛
;_200312102049(LM)
;_200312101730
;_200312101551
(defun c:tt (/ plendpt x_lst e op
n fst_pt int_num online pc
pn min_num edge_int_num pt_online_num
)
(setq min_num 0.0001)
(princ "\n选择多义线:")
(while (not e) (setq e (car (entsel))))
(if (= "LWPOLYLINE" (cdr (assoc 0 (entget e))))
(progn
(redraw e 3)
(while (not op) (setq op (getpoint "\n测试点:")))
(foreach n (entget e)
(if (= 10 (car n))
(if plendpt
(setq plendpt (append plendpt (list (cdr n))))
(setq plendpt (list (cdr n)))
)
)
)
(setq X_lst (mapcar '(lambda (x) (car x)) plendpt))
(setq delta_X (abs (- (apply 'max x_lst) (car op))))
(setq plendpt (append plendpt (list (nth 0 plendpt)))
edge_int_num 0
pt_online_num 0
)
(while (> (length plendpt) 1)
(setq pc (nth 0 plendpt)
pn (nth 1 plendpt)
)
(if (inters op
(polar op 0 (* delta_x 2))
pc
pn
)
(setq edge_int_num (+ 1 edge_int_num))
)
(if (equal (angle op pc) 0 min_num)
(setq pt_online_num (+ 1 pt_online_num))
)
(if (and (equal (angle op pc) 0 min_num)
(equal (angle op pn) 0 min_num)
)
(setq pt_online_num (- pt_online_num 1))
)
(if (equal (+ (distance op pc) (distance op pn))
(distance pc pn)
min_num
)
(setq online "t")
)
(setq plendpt (cdr plendpt))
)
(redraw e 4)
(if online
(princ "\n在多边形的边界上.")
(if (= (rem (+ pt_online_num edge_int_num) 2) 1)
(princ "\n内部")
(princ "\n外部")
)
)
)
)
(princ)
)
[/FONT]

[FONT=courier new] 
 ;_精确测试点与多边形的位置关系__纯lisp版__BY WKAI __晓东CAD论坛 
 ;_晓东CAD论坛 
 ;_200312102049(LM) 
 ;_200312101730 
 ;_200312101551 
(defun c:tt (/                plendpt           x_lst   e                 op 
             n                fst_pt           int_num    online         pc 
             pn                min_num           edge_int_num                 pt_online_num 
            ) 
  (setq min_num 0.0001) 
  (princ "\n选择多义线:") 
  (while (not e) (setq e (car (entsel)))) 
  (if (= "LWPOLYLINE" (cdr (assoc 0 (entget e)))) 
    (progn 
      (redraw e 3) 
      (while (not op) (setq op (getpoint "\n测试点:"))) 
      (foreach n (entget e) 
        (if (= 10 (car n)) 
          (if plendpt 
            (setq plendpt (append plendpt (list (cdr n)))) 
            (setq plendpt (list (cdr n))) 
          ) 
        ) 
      ) 
      (setq X_lst (mapcar '(lambda (x) (car x)) plendpt)) 
      (setq delta_X (abs (- (apply 'max x_lst) (car op)))) 
      (setq plendpt          (append plendpt (list (nth 0 plendpt))) 
            edge_int_num  0 
            pt_online_num 0 
      ) 
      (while (> (length plendpt) 1) 
        (setq pc (nth 0 plendpt) 
              pn (nth 1 plendpt) 
        ) 
        (if (inters op 
                    (polar op 0 (* delta_x 2)) 
                    pc 
                    pn 
            ) 
 
          (setq edge_int_num (+ 1 edge_int_num)) 
        ) 
        (if (equal (angle op pc) 0 min_num) 
          (setq pt_online_num (+ 1 pt_online_num)) 
        ) 
        (if (and (equal (angle op pc) 0 min_num) 
                 (equal (angle op pn) 0 min_num) 
            ) 
          (setq pt_online_num (- pt_online_num 1)) 
        ) 
        (if (equal (+ (distance op pc) (distance op pn)) 
                   (distance pc pn) 
                   min_num 
            ) 
          (setq online "t") 
        ) 
        (setq plendpt (cdr plendpt)) 
      ) 
      (redraw e 4) 
      (if online 
        (princ "\n在多边形的边界上.") 
        (if (= (rem (+ pt_online_num edge_int_num) 2) 1) 
          (princ "\n内部") 
          (princ "\n外部") 
        ) 
      ) 
    ) 
  ) 
  (princ) 
) 
[/FONT]

wkai · 发表于 2003-12-10 16:26:15

这个命题很有用，在这以前我一直没有很好的解决，大家帮我测试一下吧。
我做试验用的dwg文件。

zhynt · 发表于 2003-12-10 16:48:10

你里面使用了vl-sort 函数，如果不用呢？用纯lisp

wkai · 发表于 2003-12-10 17:06:49

调用你的排序函数就成了纯LISP呀.
：）

aeo · 发表于 2003-12-10 22:30:59

只能为多边形,不要有曲线段.

写的还能简单一点.(好理解一点)

陌生人 · 发表于 2003-12-14 05:55:26

好啊,谁写个支持曲线段的? :)

lsjjm · 发表于 2003-12-14 06:17:44

最初由 wkai 发布
[B]200312101551: 通过算法精确测试点与多边形的位置关系...[/B]

用(bpoly point) 来判别岂不是更简单?

wkai · 发表于 2003-12-14 08:40:17

最初由 lsjjm 发布
[B]

用(bpoly point) 来判别岂不是更简单? [/B]

愿闻其详
请介绍一下好么?

陌生人 · 发表于 2003-12-14 16:53:40

用bpoly可以选择实体范围么?

(续)今天用bply试了一下,觉得用这个方法还是有些漏洞.以下分析,望能赐教
1.不知bploy有没有选择物体的参数,因为用(bpoly pt)的话,就默认以全图实体为搜索边界. 当要测试的实体与其他实体相交(有交集)时候,生成的新曲线面积无法利用,如果用是否生成新物体来判断,显然又是不行,因为点在外面也可能会生成新的曲线
2.当(bpoly pt)生成实体的时候,如果碰到ellipse,spline等,不会生成新的闭合lwpolyline,而是跳出提示,问是否要生成regon. 遇到这种情况如何在程序中判断并处理?我还没想好.

除了以上用数学方法分析外(当然这样写程序会长些),我倾向于用offset一段很小的距离的方法来判断,这样写程序会短些,但是有个担忧,就是offset虽然很小,但是还是有可能(虽然可能性不大)不能生成新实体

另外想问:谁知道xdapi里面的判断是否支持曲线段?
(待续)

做了个测试:

;|方法2: offset后判断面积
|;
(defun c:ptin2 ()
;(setvar "OFFSETDIST" 0.01)
(setq plsel (entsel)
pl (car plsel)
plobj (vlax-ename->vla-object pl)
area1 (vla-get-area plobj))
(vl-cmdf ".offset" "1" plsel (getpoint "\n取点:") "");;1可按需要调整
(setq plobj2 (vlax-ename->vla-object (entlast))
area2 (vla-get-area plobj2))
(entdel (entlast))
(if (> area1 area2)
T ;点在曲线内
nil ;点在曲线外
)
)

以上程序一般情况下可用,并且支持有曲线段的封闭曲线.
但是:
1.当实体是regon时候不适用
2.当封闭曲线有自相交的时候结果不正常

lsjjm · 发表于 2003-12-15 05:22:20

最初由陌生人发布
[B]...bply试了一下,觉得用这个方法还是有些漏洞... [/B]

1. BPOLY对1楼所说"多边形" zu gou了
2. 其他曲线可使用BHATCH命令.

陌生人 · 发表于 2003-12-15 14:15:54

不够,因为用bpoly无法判断(如上所述),另外补充,即便用bpoly,bhatch,boundary 在实际使用的时候常常在封闭区域内无法得出正确的封闭曲线.当图面复杂,曲线实体多的时候尤其如此.
另外探讨一下,用bhatch,或者boundary 如果在程序中调用,使用命令行方式, 输入 /a/b/n...选择物体后无法用"" 终结选取进程,试了n次也许是bug.有谁试过是可以的么?
长老你的关于(bpoly)方面的使用是从哪里得来的?我找了下cad帮助, 以及上 point a 都没找到相关信息.

lsjjm · 发表于 2003-12-15 16:20:44

最初由陌生人发布
[B]不够,因为用bpoly无法判断...[/B]

能否举个具体的例子来说明?

陌生人 · 发表于 2003-12-15 23:32:09

实例1

lsjjm · 发表于 2003-12-16 03:18:06

我试了, 没有问题啊!

陌生人 · 发表于 2003-12-16 23:50:30

你怎么试的?能把程序贴出来么?
我不是说在这两个点无法生成封闭曲线,而是说用bpoly方法无法判断这两个点是否在那个近似星行的闭合曲线内

账号		自动登录	找回密码
密码			立即注册

[分享]：通过算法精确测试点与多边形的位置关系

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

相关帖子

Re: [分享]：通过算法精确测试点与多边形的位置关系

Re: Re: [分享]：通过算法精确测试点与多边形的位置关系