[求助]：这样判断两个线段是否共线为什么会经常漏掉一些？

liuyj · 发表于 2004-6-30 14:00:50

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

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

×

设p1，p2为线段一的两端点；p3，p4为线段二的两端点，agl是两线的角度，gcl是精度。
用如下的代码判断两线段是否共线有时会失败，为什么？

(setq p00 (polar p1 (+ (* pi 0.5) agl) 2000)
p_1 (xdrx_interst p1 p2 p1 p00 nil)
p_2 (xdrx_interst p3 p4 p1 p00 nil)
)
(equal p_1 p_2 gcl)

XDSoft · 发表于 2004-6-30 20:58:30

最初由 liuyj 发布
[B]设p1，p2为线段一的两端点；p3，p4为线段二的两端点，agl是两线的角度，gcl是精度。
用如下的代码判断两线段是否共线有时会失败，为什么？
(setq p00 (polar p1 (+ (* pi 0.5) agl) 2000)
p_1 (xd... [/B]

和大家先说说，你遇到的，什么情况下的测试会失败？最好附个图上来。

其实判断是否共线的最好方法是判断一个直线的两个端点到另一条直线的距离都等于0。

Free-Lancer · 发表于 2004-6-30 22:57:14

还可以这样

(defun lnonln (ln1 ln2 / sp1 ep1 sp2 ep2 v1 v2)
(setq sp1 (vlax-curve-getstartpoint ln1)
ep1 (vlax-curve-getendpoint ln1)
v1 (mapcar '- ep1 sp1)
sp2 (vlax-curve-getstartpoint ln2)
ep2 (vlax-curve-getendpoint ln2)
v2 (mapcar '- ep2 sp2)
)
(equal (- (* (car v1) (cadr v2))
(* (car v2) (cadr v1))
)
0.
0.00001
)
)

lijiao · 发表于 2004-7-1 08:10:26

最初由 Free-Lancer 发布
[B]还可以这样
(defun lnonln (ln1 ln2 / sp1 ep1 sp2 ep2 v1 v2)
(setq sp1 (vlax-curve-getstartpoint ln1)
ep1 (vlax-curve-getendpoint ln1)
v1 (mapcar '- ep1 sp1)
sp2 (vlax-curve-getstartp... [/B]

好像有点问题，你的方法只能判断出两线是否平行。

******************************************************************

我是用三角形法则，及两边之和大于第三边。

(defun lnonln (ln1 ln2 / sp1 ep1 sp2 ep2 v1 v2)
(setq sp1 (vlax-curve-getstartpoint ln1)
ep1 (vlax-curve-getendpoint ln1)
dis1 (distance sp1 ep1)
sp2 (vlax-curve-getstartpoint ln2)
ep2 (vlax-curve-getendpoint ln2)
dis11 (distance sp1 sp2)
dis12 (distance sp1 ep2)
dis21 (distance ep1 sp2)
dis22 (distance ep1 ep2)
)
(and (equal (* 2.0 (max dis1 dis11 dis21))
(+ dis1 dis11 dis21)
0.00001
)
(equal (* 2.0 (max dis1 dis12 dis22))
(+ dis1 dis12 dis22)
0.00001
)
)
)

lsjjm · 发表于 2004-7-1 10:26:32

3楼的程序不完全对.
请参考:
[php]
(defun c:test ()
  (setq ln1 (car (entsel)))
  (setq ln2 (car (entsel)))
  (setq d1 (dispt2line (vlax-curve-getstartpoint ln1) ln2))
  (setq d2 (dispt2line (vlax-curve-getendpoint ln1) ln2))
  (if (and (equal d1 0 0.001)(equal d2 0 0.001))
T
nil
  )
)
(defun dispt2line (p0 lline)
  (distance p0 (vlax-curve-getclosestpointto lline (trans p0 1 0) acExtendBoth))
)
[/php]

雪山飞狐(lzh) · 发表于 2004-7-1 10:44:05

偶的办法是做一条辅助线，与两条线交与同一点，且两条线的角度相等

lijiao · 发表于 2004-7-1 11:04:26

最初由 lzh741206 发布
[B]偶的办法是做一条辅助线，与两条线交与同一点，且两条线的角度相等 [/B]

这倒是个不错的方法，不过最好不要用角度来判断，用两条辅助线应该更好。

lsjjm · 发表于 2004-7-1 11:29:48

无需辅助线, 角度法:
[php]
(defun c:test ()
  (setq l1 (car (entsel)))
  (setq l2 (car (entsel)))
  (setq p1 (cdr (assoc 10 (entget l1))))
  (setq p2 (cdr (assoc 11 (entget l1))))
  (setq p3 (cdr (assoc 10 (entget l2))))
  (setq p4 (cdr (assoc 11 (entget l2))))
  (and (or (equal (angle p1 p2)(angle p3 p4) 0.001)
            (equal (angle p1 p2)(angle p4 p3) 0.001))
      (equal (angle p1 p2)(angle p1 p3) 0.001))
)
[/php]

雪山飞狐(lzh) · 发表于 2004-7-1 11:36:42

And的第二个条件(equal (angle p1 p2)(angle p1 p3) 0.001))有点问题，不过思路很好

Free-Lancer · 发表于 2004-7-1 12:16:59

看来大家很热心，俺就在说几句，直线的方程为 y=kx+b ，对两个线段，只要判断 k 斜率、b 截距相同就是共线的，已知两点求 k b 是很容易的

liuyj · 发表于 2004-7-1 12:45:53

谢谢大家！其实我是用了天正的一个子函数，用xdapi替换后就是：
;;检验两线段是否在一直线上;name1,name2--实体名;

(defun l_line1 (name1 name2 gcl / agl p1 p2 p3 p4 p00 p_1 p_2)
(xdrx_setenttodb name1)
(setq p1 (xdrx_getentdxf 10)
p2 (xdrx_getentdxf 11)
agl (angle p1 p2)
)
(xdrx_setenttodb name2)
(setq p3 (xdrx_getentdxf 10)
p4 (xdrx_getentdxf 11)
p00 (polar p1 (+ (* pi 0.5) agl) 2000)
p_1 (xdrx_interst p1 p2 p1 p00 nil)
p_2 (xdrx_interst p3 p4 p1 p00 nil)
)
(equal p_1 p_2 gcl)
)

有时明明两条线段是共线的，可是他就判断不出，从思路上讲程序是可行的，我找不到原因。
从明经上看到过lucas的一段程序：

(setq VLINE1_10 (cdr (assoc 10 (entget VLINE1))))
(setq VLINE1_11 (cdr (assoc 11 (entget VLINE1))))
(setq VLINE2_10 (cdr (assoc 10 (entget VLINE2))))
(setq VLINE2_11 (cdr (assoc 11 (entget VLINE2))))
(setq ANG1 (angle VLINE1_10 VLINE1_11))
(setq ANG2 (angle VLINE1_10 VLINE2_10))
(setq ANG3 (+ pi ANG1))
(if (or (equal ANG1 ANG2 0.000000001)
(equal ANG1 (+ ANG2 pi) 0.000000001)
(equal ANG3 ANG2 0.000000001)
(equal ANG3 (+ ANG2 pi) 0.000000001)
)......)

和lsjjm长老的思路大体相同，不过我觉得效率低了点，如果框选很多条平行的中间断开的线段，
明显会慢，前些日子ea斑竹做的东东速度很快，不知是用什么算法？

Free-Lancer · 发表于 2004-7-1 13:11:34

最初由 liuyj 发布
[B]谢谢大家！其实我是用了天正的一个子函数，用xdapi替换后就是：
;;检验两线段是否在一直线上;name1,name2--实体名;
(defun l_line1 (name1 name2 gcl / agl p1 p2 p3 p4 p00 p_1 p_2)
(xdrx_setenttodb name1)
... [/B]

唉，你怎么不试试10楼的方法？

一般说来，速度的快慢和函数中的求值的步骤多少有关，比如你写的 l_line1 没必要 setenttodb ，直接用 curve-getstartpoint 即可，你的写法反倒浪费了一点点时间。程序的运行效率都是这样一点点积累起来的。

其实3楼的方法只要最后的abd判断中再加一句也可以

(or
(vlax-curve-ison ln1 sp2)
(vlax-curve-ison ln1 ep2)
(vlax-curve-ison ln2 sp1)
(vlax-curve-ison ln2 ep1)
)

复制代码

Free-Lancer · 发表于 2004-7-1 13:33:11

最初由陌生人发布
[B]用xd说的方法,在简化一下，不求距离，直接判断点是否相等：
[code]
(defun atline (e1 e2 / p2s p2e)
(setq p2s (vlax-curve-getstartpoint e2)
p2e (vlax-curve-getendpoint e2))
(and (equal p2s (vlax-... [/B]

你的方法不对，在一条线上时 closestpointto 返回的不是垂点，而是另外一条线的起始点和终点，api 中有点到直线的垂点函数吗？

陌生人 · 发表于 2004-7-1 13:36:48

哈哈，你说话的语气让我想起某个人。
不过，在你说不对之后。请你好好看看vlax-curve-getclosestpointto 最后一个参数是 T 的时候，帮助文件里面是怎么说的，然后再试试上面的这段程序

Free-Lancer · 发表于 2004-7-1 14:49:01

最初由陌生人发布
[B]哈哈，你说话的语气让我想起某个人。
不过，在你说不对之后。请你好好看看vlax-curve-getclosestpointto 最后一个参数是 T 的时候，帮助文件里面是怎么说的，然后再试试上面的这段程序 [/B]

呵呵，你对了。建议加精华并奖励

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

[求助] [求助]：这样判断两个线段是否共线为什么会经常漏掉一些？

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

相关帖子

Re: [求助]：这样判断两个线段是否共线为什么会经常漏掉一些？