过2点的构造线是否切割封闭多段线

/db_自贡黄明儒_ · 发表于 2021-2-7 08:51:10

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

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

×

本帖最后由 /db_自贡黄明儒_ 于 2021-2-9 19:38 编辑

已知：两点p1、P2和封闭多段线(不带弧)。判断:过p1 p2的构造线(不存在）是否切割封闭多段线。

Lisphk · 发表于 2021-2-7 08:59:27

判断过p1,p2的直线是否和多段线大于等于两个交点不就行了吗？

/db_自贡黄明儒_ · 发表于 2021-2-7 09:12:29

Lisphk 发表于 2021-2-7 08:59
判断过p1,p2的直线是否和多段线大于等于两个交点不就行了吗？

您的意思是：先产生一条直线，然后决断直线与多段线有无交点？

Lisphk · 发表于 2021-2-7 09:17:58

本帖最后由 Lisphk 于 2021-2-7 09:19 编辑

/db_自贡黄明儒_ 发表于 2021-2-7 09:12
您的意思是：先产生一条直线，然后决断直线与多段线有无交点？

可以啊，用完删除嘛
也可以用数学方法，就不用建立直线

lijiao · 发表于 2021-2-7 09:20:02

判断多段线的每一个端点相对于p1、p2是左转还是右转，如果所有的转向一致，则不相切割，如果转向不一致，则相割。

Lisphk · 发表于 2021-2-7 09:31:36

lijiao 发表于 2021-2-7 09:20
判断多段线的每一个端点相对于p1、p2是左转还是右转，如果所有的转向一致，则不相切割，如果转向不一致，则 ...

条件不够吧

/db_自贡黄明儒_ · 发表于 2021-2-7 10:00:50

本帖最后由 /db_自贡黄明儒_ 于 2021-2-7 19:55 编辑

lijiao 发表于 2021-2-7 09:20
判断多段线的每一个端点相对于p1、p2是左转还是右转，如果所有的转向一致，则不相切割，如果转向不一致，则 ...

版主就是版主，佩服！！我来一个罗索版的,似乎没有版主的严谨

;;-----------------------------------------------------------;;
;; 旋转一个向量或者点90度 By highflybird ;;
;; 输入: 一个向量 ;;
;; 输出: 被旋转90度后的向量 ;;
;;-----------------------------------------------------------;;
(defun MAT:Rot90 (vec)
(vl-list* (- (cadr vec)) (car vec) (cddr vec))
)
;;164.3 [功能] 多段线端点列表
;;示例(HH:PtLists (car (entsel)))
(defun HH:PtLists (en)
(mapcar 'cdr
(vl-remove-if-not '(lambda (x) (= (car x) 10)) (entget en))
)
)
;;判断过2的构造线是否切割封闭多段线
(defun 2P-LWPOLYLINE-P1 (pts p1 p2)
(apply 'or
(mapcar '(lambda (x y) (apply 'inters (list p1 p2 x y))) pts (cdr pts))
)
)
(defun 2P-LWPOLYLINE-P (pts p1 p2 / FLAG P3 P4 P5 P6)
(setq pts (cons (last pts) pts))
(setq Flag (2P-LWPOLYLINE-P1 pts p1 p2))
(if (not Flag)
(progn
(setq p3 (mapcar '+ p1 (MAT:Rot90 (mapcar '- p1 p2))))
(setq p4 (mapcar '+ p2 (MAT:Rot90 (mapcar '- p1 p2))))
(setq Flag (2P-LWPOLYLINE-P1 pts p3 p4))
(if (not Flag)
(progn
(setq p5 (mapcar '+ p1 (MAT:Rot90 (mapcar '- p2 p1))))
(setq p6 (mapcar '+ p2 (MAT:Rot90 (mapcar '- p2 p1))))
(setq Flag (2P-LWPOLYLINE-P1 pts p5 p6))
)
)
)
)
Flag
)
(defun C:w11 ()
(setq e (car (entsel)))
(setq pts (HH:PtLists e))
(setq p1 (getpoint))
(setq p2 (getpoint))
(princ (2P-LWPOLYLINE-P pts p1 p2))
(princ)
)

;;-----------------------------------------------------------;; 
;; 旋转一个向量或者点90度     By highflybird                 ;; 
;; 输入: 一个向量                                            ;; 
;; 输出: 被旋转90度后的向量                                  ;; 
;;-----------------------------------------------------------;; 
(defun MAT:Rot90 (vec) 
  (vl-list* (- (cadr vec)) (car vec) (cddr vec)) 
) 
;;164.3 [功能] 多段线端点列表 
;;示例(HH:PtLists (car (entsel))) 
(defun HH:PtLists (en) 
  (mapcar 'cdr 
    (vl-remove-if-not '(lambda (x) (= (car x) 10)) (entget en)) 
  ) 
) 
 
;;判断过2的构造线是否切割封闭多段线 
(defun 2P-LWPOLYLINE-P1  (pts p1 p2) 
  (apply 'or 
   (mapcar '(lambda (x y) (apply 'inters (list p1 p2 x y))) pts (cdr pts)) 
  ) 
) 
(defun 2P-LWPOLYLINE-P (pts p1 p2 / FLAG P3 P4 P5 P6) 
  (setq pts (cons (last pts) pts)) 
  (setq Flag (2P-LWPOLYLINE-P1 pts p1 p2)) 
  (if (not Flag) 
    (progn 
      (setq p3 (mapcar '+ p1 (MAT:Rot90 (mapcar '- p1 p2)))) 
      (setq p4 (mapcar '+ p2 (MAT:Rot90 (mapcar '- p1 p2)))) 
      (setq Flag (2P-LWPOLYLINE-P1 pts p3 p4)) 
      (if (not Flag) 
  (progn 
    (setq p5 (mapcar '+ p1 (MAT:Rot90 (mapcar '- p2 p1)))) 
    (setq p6 (mapcar '+ p2 (MAT:Rot90 (mapcar '- p2 p1)))) 
    (setq Flag (2P-LWPOLYLINE-P1 pts p5 p6)) 
  ) 
      ) 
    ) 
  ) 
  Flag 
) 
 
(defun C:w11 () 
  (setq e (car (entsel))) 
  (setq pts (HH:PtLists e)) 
  (setq p1 (getpoint)) 
  (setq p2 (getpoint)) 
  (princ (2P-LWPOLYLINE-P pts p1 p2)) 
  (princ) 
)

st788796 · 发表于 2021-2-7 15:42:04

本帖最后由 st788796 于 2021-2-7 15:43 编辑

(defun iscross (pline p1 p2 / n i pts tf)
(setq n (fix (vlax-curve-getendparam pline))
i 0
)
(repeat n
(setq pts (cons (vlax-curve-getpointatparam pline i) pts)
i (1+ i)
)
)
(if (or (vlax-curve-isclosed pline)
(equal (vlax-curve-getendpoint pline)
(vlax-curve-getstartpoint pline)
0.00001
)
)
(setq pts (cons (vlax-curve-getstartpoint pline) pts))
)
(setq tf t)
(mapcar '(lambda (x y)
(if (and tf
(inters x y p1 p2)
)
(setq tf nil)
)
)
pts
(cdr pts)
)
(not tf)
)

/db_自贡黄明儒_ · 发表于 2021-2-7 22:28:35

st788796 发表于 2021-2-7 15:42

怎么感觉您这代码同我上面的一样不严谨呢？是不是我的这题太简单了。

newer · 发表于 2021-2-7 23:33:55

/db_自贡黄明儒_ 发表于 2021-2-7 22:28
怎么感觉您这代码同我上面的一样不严谨呢？是不是我的这题太简单了。

这个或许能帮助

用WINDING数判断点是否在多边形内部的算法及C++/LISP实现代码

1 个回复 - 279 次查看

Determining the inclusion of a point P in a 2D planar polygon is a geometric problem that results in interesting algorithms. Two commonly used methods are:

The Crossing Number (cn) method - which ...

2017-7-14 17:57 - newer - 算法与数据结构

/db_自贡黄明儒_ · 发表于 2021-2-8 12:37:20

newer 发表于 2021-2-7 23:33
这个或许能帮助

2个点有一个在内部，肯定是切割多段线了。主要问题两个都在外时的判断。

;;;;;164.3 [功能] 多段线端点列表
;;;;;示例(HH:PtLists (car (entsel)))
;;;(defun HH:PtLists (en)
;;; (mapcar 'cdr
;;; (vl-remove-if-not '(lambda (x) (= (car x) 10)) (entget en))
;;; )
;;;)
;;;(defun det (p1 p2 p3 / x2 y2)
;;; (setq x2 (car p2)
;;; y2 (cadr p2)
;;; )
;;; (- (* (- x2 (car p3)) (- y2 (cadr p1)))
;;; (* (- x2 (car p1)) (- y2 (cadr p3)))
;;; )
;;;)
;;过2点的构造线与封闭多段线相交，返回T
(defun 2P-LWPOLYLINE-P (pts p1 p2 / L)
(setq L (mapcar '(lambda (x) (MINUSP (det p1 p2 x))) pts))
(setq L (cons (last L) L))
(not (VL-EVERY '= L (cdr L)))
)

st788796 · 发表于 2021-2-8 14:54:32

/db_自贡黄明儒_ 发表于 2021-2-8 12:37
2个点有一个在内部，肯定是切割多段线了。主要问题两个都在外时的判断。

和任意一个段有交点就是分割，如用xline，inters + nil 参数

lijiao · 发表于 2021-2-9 10:02:13

以P1为基准点,将多边形旋转一个角度,角度为(* (angle p1 p2) -1.0),然后求多边形的最小包围框,包围框的纵坐标y1<y2,如果P1点的纵坐标在y1与y2之间,则为相割

/db_自贡黄明儒_ · 发表于 2021-2-9 11:52:13

lijiao 发表于 2021-2-9 10:02
以P1为基准点,将多边形旋转一个角度,角度为(* (angle p1 p2) -1.0),然后求多边形的最小包围框,包围框的纵坐 ...

这个方法就不如你上次说的方法了。

还是直接计算，可能效率更高

newer · 发表于 2021-2-9 12:41:34

/db_自贡黄明儒_ 发表于 2021-2-9 11:52
这个方法就不如你上次说的方法了。

还是直接计算，可能效率更高

这要看你定义的两点是直线，还是线段了

如果是线段，那么就要判断的多一些

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

[研讨] 过2点的构造线是否切割封闭多段线

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

评分

评分

点评

点评