[原创] 判断点在封闭曲线的位置

caoyin · 发表于 2009-9-16 21:35:00

本帖最后由作者于 2009-9-16 21:55:26 编辑

回不死猫：拟合曲线还是存在精度问题，但速度应该很快

前几天的讨论http://bbs.mjtd.com/forum.php?mod=viewthread&tid=77921，都是求点集，直接求图形对象，特别好的算法可能不是很多，等待高手...

自贡黄明儒 · 发表于 2011-5-10 18:13:57

太好了，用此程序可以完成我的"局部放大"了.谢谢caoyin

狂刀lxx · 发表于 2011-5-10 22:05:15

我是射线党。

luoguoqianga · 发表于 2011-5-13 08:28:37

选择封闭曲线:选择一个圆，然后点圆心处，发现：
拾取任意点: ......在曲线外。
这显然错误。

chlh_jd · 发表于 2011-8-5 20:53:39

本帖最后由 chlh_jd 于 2011-8-5 20:54 编辑

楼主的判别函数判别效果非常精确，似乎漏了回收2个变量LST1 LST2；
如果需要大量运行，效率可能会略低些，有没有办法调用ACAD自身的呢？

chlh_jd · 发表于 2011-8-5 21:01:23

本帖最后由 chlh_jd 于 2011-8-5 21:03 编辑

就多边形而言，我一直用狂刀前辈的代码；期待能早日见到支持各类型曲线的纯LISP版或调用ACAD自身的版本

;;;判断点是否在封闭多边形内 .
;;;Function : judge a point location with polygon
;;;Arg : pt -- a point
;;; pts -- points of polygon
;;; eps -- allowance
;;;return :
;;; -1 -- out of polygon , 0 -- at , 1 -- in
(defun pipl?(pt pts eps / is at a)
;; by 狂刀
;; Edit by GSLS(SS) 2011.03.28
;; Solved the problem : if a point at the given polygon , it perhap return T or NIL .
(setq pt(trans pt 0 0)pts(mapcar(function(lambda(x)(trans x 0 0)))pts)Eps(abs Eps))
(if(equalmember pt pts eps)0
(progn(setq is(equal PI (abs(apply(function +)(mapcar(function(lambda (x y / a)
(setq a (rem (- (angle pt x) (angle pt y)) PI))
(if(equal (+(distance pt x)(distance pt y))(distance x y)Eps)(setq at T))a));_(check-pt (list x y))
(cons (last pts) pts) pts))) eps))(cond (at 0)(is 1)(T -1)))))
(defun equalmember (e lst tor /)
(apply (function or)
(mapcar (function (lambda (x) (equal x e tor))) lst)
)
)

;;;判断点是否在封闭多边形内 . 
;;;Function : judge a point location with polygon 
;;;Arg : pt -- a point 
;;;      pts -- points of polygon 
;;;      eps -- allowance 
;;;return : 
;;;     -1 -- out of polygon , 0 -- at , 1 -- in 
(defun pipl?(pt pts eps / is at a) 
  ;; by 狂刀   
  ;; Edit by GSLS(SS) 2011.03.28 
  ;; Solved the problem : if a point at the given polygon , it perhap return T or NIL . 
  (setq pt(trans pt 0 0)pts(mapcar(function(lambda(x)(trans x 0 0)))pts)Eps(abs Eps)) 
  (if(equalmember pt pts eps)0 
    (progn(setq is(equal PI (abs(apply(function +)(mapcar(function(lambda (x y / a) 
        (setq a (rem (- (angle pt x) (angle pt y)) PI)) 
                (if(equal (+(distance pt x)(distance pt y))(distance x y)Eps)(setq at T))a));_(check-pt (list x y)) 
                   (cons (last pts) pts) pts))) eps))(cond (at 0)(is 1)(T -1))))) 
(defun equalmember (e lst tor /) 
  (apply (function or) 
         (mapcar (function (lambda (x) (equal x e tor))) lst) 
  ) 
)

yanglin112 · 发表于 2011-10-27 11:33:02

好帖，学习一下！

434939575 · 发表于 2012-3-23 23:43:22

这个什么用啊？

nzl1116 · 发表于 2012-3-24 09:02:33

最简单的方法用(bpoly pt)

wjl1014 · 发表于 2013-3-11 16:47:44

可否用多边形选择集判断？选中在内，没选中在外？

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