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

lsjjm

最初由 陌生人 发布
[B]你怎么试的?能把程序贴出来么?...[/B]

很简单的:
(defun ptinside (pt)
  (if (bpoly p)
    (progn (vl-cmdf "erase" (entlast) "") T)
    nil
  )
)

wkai

最初由 lsjjm 发布
[B]
很简单的:
(defun ptinside (pt)
  (if (bpoly p)
    (progn (vl-cmdf "erase" (entlast) "") T)
    nil
  )
) [/B]

bpoly这个函数没有见过呀.
请问，lsjjm 大侠 AUTOLISP 中到底还有多少这种没有介绍的函数呢？

lsjjm

BPOLY有两种形式: 命令BPOLY和函数(BPOLY PT).
这是R12以后增加的用于补充和加强HATCH的功能.新版中其功能为BOUNADRY所代替, 旧的命令及函数仍保留, 但在帮助文件里是看不到了.

陌生人

最初由 lsjjm 发布
[B]
很简单的:
(defun ptinside (pt)
  (if (bpoly p)
    (progn (vl-cmdf "erase" (entlast) "") T)
    nil
  )
) [/B]

对啊,这正是我要说的,这样的测试无法确定点是否在某个封闭曲线内部.因为如果你点在这个曲线的外面,也有可能获得(bpoly pt)生成的实体,(当这个点的位置被其他实体"封闭"在某个封闭空间内的时候)
所以,bpoly 是无法用来判定的.
如图例中的两个点,都能正确生成新的曲线,但是相对于星形的封闭曲线,他们一个在内,一个在外.
看来你没明白我所提出的问题.


来了3次,我还是最后一贴,只好用编辑:(
今天做了个测试点在曲线内外的程序,是用构造xline,求交点排序位置的方法.

测试:
(c:ptin)
出错! 未选中曲线
_$ (c:ptin)
点在曲线外
_$ (c:ptin)
点在曲线内
_$ (c:ptin)
点在曲线上:
_$

附件为.fas文件. 拖入cad作图区即可加载. 调用用: ptin

lsjjm

最初由 陌生人 发布
[B]...他们一个在内,一个在外.看来你没明白我所提出的问题...[/B]

我怎么会不明白?
如此简单的程序只是显示一个基本的功能或者si 路.  还需要你发hui和完shan. 例如隐藏其他实体...?

陌生人

听你一说,我茅塞顿开.
这个方法我原来连想都没有去想,因为为了一个检测而对全图其他实体隐藏再显示,真是杀鸡用牛刀.如果图大,redraw一便也是会感觉明显的延迟,即便机子再快.
不过作为思维的严谨性,这个确实是一种思路,虽然打死我不会这样去作.

再:我的观点:用bpoly还是无法做出测试程序!
因为如果对 封闭的spline,ellipse测试, 会被迫中断退出.必须回答那个
讨厌的"Polyline boundary could not be derived. Create Region? <Y>"
我试了,这个回答(y)无法在程序的命令行语句中完成.
这样程序就无法正常进行下去,而且后续的恢复实体显示也只好在退出后用另外的程序来完成.
测试如下:
Command: el ELLIPSE
Specify axis endpoint of ellipse or [Arc/Center]:
Specify other endpoint of axis:
Specify distance to other axis or [Rotation]:
Command: (bpoly (getpoint))
<图元名: 400ad9f8>
Polyline boundary could not be derived. Create Region? <Y>

1 loop extracted.

lsjjm

1. 这是一种思路, 但更是一种解决问题的方法, 而且是不错的方法( 喜好繁琐, 有误的计算?)
  解决问题有许多方法, gan ma fei de zai yi ke shu shang diao si ?!
2. 我测试了ELLIPSE, SPLINE完全可以 !
3. 测试上面那个程序了吗?

陌生人

测试了.我用的是cad2002  ,用 (ptinside (getpoint)) 调用
对el,sp无论点在内,在外都返回nil, 你用cad2004?
其他人效果如何?不会是我的cad问题吧?呵呵,别的功能没什么异样啊.或者是有什么参数?
不管了,先睡觉.88

lsjjm

我用的是2000i和2004都没问题, 2002也不会有问题的! 检查检查你的CAD系统吧?

陌生人

还有没有其他兄弟?也来测试下?

小菜

楼主的程序对于附图的情况，直线的两个端点在矩形外面，可是测出来一个在内，一个在外

小菜

斗胆改了一下程序，不知道改得对不对，请大家指教:

2. 
   
3.   [FONT=courier new]
   
4. ;以测试点为起点作超过多边形外框的线段A。
   
5. ;除了测线段A与多边形的各边的交点数（M）,
   
6. ;还需要测落在在线段A上的多边形的端点数（N）,
   
7. ;以及重合在线段A的多边形边线的条数（O）
   
8. ;如果（REM （- （+ M N） O）2）为1则在多边形内部，否则就在外部。
   
9. 
   
10. ;_精确测试点与多边形的位置关系__纯lisp版__BY WKAI __晓东CAD论坛
    
11. ;_晓东CAD论坛
    
12. ;_200312102049(LM)
    
13. ;_200312101730
    
14. ;_200312101551
    
15. ;_20040301修改--by 小菜
    
16. (defun c:tt (/                plendpt           x_lst   e                 op
    
17.              n                fst_pt           int_num    online         pc
    
18.              pn                min_num           edge_int_num                 pt_online_num
    
19.              op1 )
    
20.   (setq min_num 0.0001)
    
21.   (princ "\n选择多义线:")
    
22.   (while (not e) (setq e (car (entsel))))
    
23.   (if (= "LWPOLYLINE" (cdr (assoc 0 (entget e))))
    
24.     (progn
    
25.       (redraw e 3)
    
26.       (while (not op) (setq op (getpoint "\n测试点:")))
    
27.       (foreach n (entget e)
    
28.         (if (= 10 (car n))
    
29.           (if plendpt
    
30.             (setq plendpt (append plendpt (list (cdr n))))
    
31.             (setq plendpt (list (cdr n)))
    
32.           )
    
33.         )
    
34.       )
    
35.       (setq X_lst (mapcar '(lambda (x) (car x)) plendpt))
    
36.       (setq delta_X (abs (- (apply 'max x_lst) (car op))))
    
37.       (setq plendpt          (append plendpt (list (nth 0 plendpt)))
    
38.             edge_int_num  0
    
39.             pt_online_num 0
    
40.       )
    
41.       (while (> (length plendpt) 1)
    
42.         (setq pc (nth 0 plendpt)
    
43.               pn (nth 1 plendpt)
    
44.         )
    
45.         (setq op1 (polar op 0 (* delta_x 2)))
    
46.         (if (inters op
    
47.                     op1
    
48.                     pc
    
49.                     pn
    
50.             )
    
51. 
    
52.           (setq edge_int_num (+ 1 edge_int_num))
    
53.         )
    
54.         (if (equal (+ (distance op pc) (distance pc op1))
    
55.                    (distance op op1)
    
56.                    min_num)
    
57.           (setq pt_online_num (+ 1 pt_online_num))
    
58.         )
    
59. ;        (if (and (equal (angle op pc) pi min_num)
    
60. ;                 (equal (angle op pn) pi min_num)
    
61. ;            )
    
62. ;          (setq pt_online_num (- pt_online_num 1))
    
63. ;        )
    
64.         (if (equal (+ (distance op pc) (distance op pn))
    
65.                    (distance pc pn)
    
66.                    min_num
    
67.             )
    
68.           (setq online "t")
    
69.         )
    
70.         (setq plendpt (cdr plendpt))
    
71.       )
    
72.       (redraw e 4)
    
73.       (if online
    
74.         (princ "\n在多边形的边界上.")
    
75.         (if (= (rem (+ pt_online_num edge_int_num) 2) 1)
    
76.           (princ "\n内部")
    
77.           (princ "\n外部")
    
78.         )
    
79.       )
    
80.     )
    
81.   )
    
82.   (princ)
    
83. )
    
84. 
    
85.   [/FONT]
    

<br />   [FONT=courier new]<br /> ;以测试点为起点作超过多边形外框的线段A。<br /> ;除了测线段A与多边形的各边的交点数（M）,<br /> ;还需要测落在在线段A上的多边形的端点数（N）,<br /> ;以及重合在线段A的多边形边线的条数（O）<br /> ;如果（REM （- （+ M N） O）2）为1则在多边形内部，否则就在外部。<br /> <br /> ;_精确测试点与多边形的位置关系__纯lisp版__BY WKAI __晓东CAD论坛<br /> ;_晓东CAD论坛<br /> ;_200312102049(LM)<br /> ;_200312101730<br /> ;_200312101551<br /> ;_20040301修改--by 小菜<br /> (defun c:tt (/                plendpt           x_lst   e                 op<br />              n                fst_pt           int_num    online         pc<br />              pn                min_num           edge_int_num                 pt_online_num<br />              op1 )<br />   (setq min_num 0.0001)<br />   (princ "\n选择多义线:")<br />   (while (not e) (setq e (car (entsel))))<br />   (if (= "LWPOLYLINE" (cdr (assoc 0 (entget e))))<br />     (progn<br />       (redraw e 3)<br />       (while (not op) (setq op (getpoint "\n测试点:")))<br />       (foreach n (entget e)<br />         (if (= 10 (car n))<br />           (if plendpt<br />             (setq plendpt (append plendpt (list (cdr n))))<br />             (setq plendpt (list (cdr n)))<br />           )<br />         )<br />       )<br />       (setq X_lst (mapcar '(lambda (x) (car x)) plendpt))<br />       (setq delta_X (abs (- (apply 'max x_lst) (car op))))<br />       (setq plendpt          (append plendpt (list (nth 0 plendpt)))<br />             edge_int_num  0<br />             pt_online_num 0<br />       )<br />       (while (> (length plendpt) 1)<br />         (setq pc (nth 0 plendpt)<br />               pn (nth 1 plendpt)<br />         )<br />         (setq op1 (polar op 0 (* delta_x 2)))<br />         (if (inters op<br />                     op1<br />                     pc<br />                     pn<br />             )<br /> <br />           (setq edge_int_num (+ 1 edge_int_num))<br />         )<br />         (if (equal (+ (distance op pc) (distance pc op1))<br />                    (distance op op1)<br />                    min_num)<br />           (setq pt_online_num (+ 1 pt_online_num))<br />         )<br /> ;        (if (and (equal (angle op pc) pi min_num)<br /> ;                 (equal (angle op pn) pi min_num)<br /> ;            )<br /> ;          (setq pt_online_num (- pt_online_num 1))<br /> ;        )<br />         (if (equal (+ (distance op pc) (distance op pn))<br />                    (distance pc pn)<br />                    min_num<br />             )<br />           (setq online "t")<br />         )<br />         (setq plendpt (cdr plendpt))<br />       )<br />       (redraw e 4)<br />       (if online<br />         (princ "\n在多边形的边界上.")<br />         (if (= (rem (+ pt_online_num edge_int_num) 2) 1)<br />           (princ "\n内部")<br />           (princ "\n外部")<br />         )<br />       )<br />     )<br />   )<br />   (princ)<br /> )<br /> <br />   [/FONT]<br />

lsjjm

最初由 小菜 发布
[B]楼主的程序对于附图的情况，直线的两个端点在矩形外面，可是测出来一个在内，一个在外 [/B]

试试这个:
(defun isinside (p ent)
  (command "_.area" "_O" ent)
  (setq earea (getvar "AREA"))
  (command "_.offset" 0.01 ent p "" "_.area" "_O"
    (setq nent (entlast))
  )
  (entdel nent)
  (> earea (getvar "AREA"))
)

小菜

楼上这个程序不对吧？如果点在另一个矩形内，而在我想测的这个矩形外面，结果就不对了

lsjjm

最初由 小菜 发布
[B]楼上这个程序不对吧？如果点在另一个矩形内，而在我想测的这个矩形外面，结果就不对了 [/B]

改了一下再试试.