点到直线的距离，求大神优化，第一次写的

13648893846

点到线的距离，求大神优化，第一次写的

xyp1964

;; 点到两点连线的距离
(defun abc (p0 p1 p2 / p3)
  (setq p3 (polar p0 (+ (angle p1 p2) (* pi 0.5)) 1))
  (distance (inters p1 p2 p0 p3 nil) p0)
)

vectra

1. (defun c:tt ()
   
2.   (setq p1 (getpoint "\n请选择直线起点："))
   
3.   (setq p2 (getpoint "\n请选择直线第二点："))
   
4.   (setq p3 (getpoint "\n请输入待求点："))
   
5. 
   
6.   (abs (car (trans (mapcar '- p3 p1) 0 (mapcar '- p2 p1))))
   
7. )

(defun c:tt ()<br />   (setq p1 (getpoint "\n请选择直线起点："))<br />   (setq p2 (getpoint "\n请选择直线第二点："))<br />   (setq p3 (getpoint "\n请输入待求点："))<br /> <br />   (abs (car (trans (mapcar '- p3 p1) 0 (mapcar '- p2 p1))))<br /> )

http://bbs.mjtd.com/forum.php?mo ... 954&highlight=trans

yuanziyou

用内置的曲线函数：

1. (defun c:fgg ();点到直线的距离
   
2.   (vl-load-com)
   
3.   (setq objen (car (entsel "请选择直线")))
   
4.   (setq pt2 (getpoint "\n指定点:"))
   
5.   (setq pt (vlax-curve-getClosestPointTo objen pt2 T))
   
6.   (setq d1 (rtos (distance pt pt2) 2 3))
   
7.   (setq d2 (rtos (distance pt (vlax-curve-getStartPoint objen)) 2 3))
   
8.   (princ "\n垂足到起点的距离D2为=") (princ d2) (princ)  
   
9.   (princ "\n待求点到直线垂足的距离D1为=") (princ d1) (princ)
   
10. )

(defun c:fgg ();点到直线的距离<br />   (vl-load-com)<br />   (setq objen (car (entsel "请选择直线")))<br />   (setq pt2 (getpoint "\n指定点:"))<br />   (setq pt (vlax-curve-getClosestPointTo objen pt2 T))<br />   (setq d1 (rtos (distance pt pt2) 2 3))<br />   (setq d2 (rtos (distance pt (vlax-curve-getStartPoint objen)) 2 3))<br />   (princ "\n垂足到起点的距离D2为=") (princ d2) (princ)  <br />   (princ "\n待求点到直线垂足的距离D1为=") (princ d1) (princ)<br /> )

因为我最近也在写关于点到多段线的,也可以用矢量法：

1. (defun c:p2l ();点到直线的距离  (setq pt1 (getpoint "\n请选择直线起点："))
   
2.   (setq pt2 (getpoint "\n请选择直线第二点："))  
   
3.   (setq pt3 (getpoint "\n请输入待求点："))
   
4.   (setq x1 (car pt1));获取点列表的第一个元素
   
5.   (setq y1 (cadr pt1));获取点列表的第二个元素
   
6.   (setq x2 (car pt2));获取点列表的第一个元素
   
7.   (setq y2 (cadr pt2));获取点列表的第二个元素
   
8.   (setq x3 (car pt3));获取点列表的第一个元素
   
9.   (setq y3 (cadr pt3));获取点列表的第二个元素
   
10.   (setq d1 (P2LDist x1, y1, x2, y2, x3, y3)
    
11.                 d3 (distance pt1 pt3)
    
12.   )
    
13.   (if (< (- d3 (abs d1)) 0.00001)
    
14.         (setq d2 0)
    
15.           (setq d2 (sqrt (- (* d3 d3) (* d1 d1))))
    
16.    )
    
17.   (princ "\n垂足到起点的距离D2为=") (princ (rtos d2 2 3)) (princ)  
    
18.   (princ "\n待求点到直线垂足的距离D1为=") (princ (rtos (abs d1) 2 3)) (princ)
    
19. )
    
20. (defun P2LDist(x1, y1, x2, y2, x3, y3);点到直线的距离左+右-
    
21.         (setq r0(Sqrt (+ (* ( - y2 y1) ( - y2 y1)) (* (- x1 x2) (- x1 x2)))));直线起点到终点的距离
    
22.         (setq a1(-(* x1 y2) (* x3 (- y2 y1)) (* y3 (- x1 x2)) (* x2 y1)))
    
23.         (setq distL (/ a1 r0))
    
24. )

(defun c:p2l ();点到直线的距离  (setq pt1 (getpoint "\n请选择直线起点："))<br />   (setq pt2 (getpoint "\n请选择直线第二点："))  <br />   (setq pt3 (getpoint "\n请输入待求点："))<br />   (setq x1 (car pt1));获取点列表的第一个元素<br />   (setq y1 (cadr pt1));获取点列表的第二个元素<br />   (setq x2 (car pt2));获取点列表的第一个元素<br />   (setq y2 (cadr pt2));获取点列表的第二个元素<br />   (setq x3 (car pt3));获取点列表的第一个元素<br />   (setq y3 (cadr pt3));获取点列表的第二个元素<br />   (setq d1 (P2LDist x1, y1, x2, y2, x3, y3)<br />                 d3 (distance pt1 pt3)<br />   )<br />   (if (< (- d3 (abs d1)) 0.00001) <br />         (setq d2 0)<br />           (setq d2 (sqrt (- (* d3 d3) (* d1 d1))))<br />    )<br />   (princ "\n垂足到起点的距离D2为=") (princ (rtos d2 2 3)) (princ)  <br />   (princ "\n待求点到直线垂足的距离D1为=") (princ (rtos (abs d1) 2 3)) (princ)<br /> )<br /> (defun P2LDist(x1, y1, x2, y2, x3, y3);点到直线的距离左+右-<br />         (setq r0(Sqrt (+ (* ( - y2 y1) ( - y2 y1)) (* (- x1 x2) (- x1 x2)))));直线起点到终点的距离<br />         (setq a1(-(* x1 y2) (* x3 (- y2 y1)) (* y3 (- x1 x2)) (* x2 y1)))<br />         (setq distL (/ a1 r0))<br /> )

asen

太复杂, 给你个简单思路

假设A为点, B和C为两点直线

求BC角度, 然后加90度求过A点的直线上的一点(polar), 然后求四点交点(inter)即为垂足点, 然后求A到垂足点距离(distance)

asen

另一种方法求垂足点

(defun ::perPoint(p0 p1 p2 / a)
    (setq a(mapcar '- p1 p2))
    (inters(mapcar '+ (list(-(cadr a))(car a)(caddr a)) p0) p0 p1 p2 nil)
)

13648893846

asen 发表于 2017-10-6 15:50
太复杂, 给你个简单思路

假设A为点, B和C为两点直线

那个已经能算，因为是第一次写，只是要优化而已

13648893846

asen 发表于 2017-10-6 17:35
另一种方法求垂足点

(defun ::perPoint(p0 p1 p2 / a)

谢谢，又学得一个大招

13648893846

yuanziyou 发表于 2017-10-6 16:38
用内置的曲线函数：

你看看你的能计算这个位置吗

yuanziyou

13648893846 发表于 2017-10-6 19:01
你看看你的能计算这个位置吗

重新修改了下,你再试试

13648893846

yuanziyou 发表于 2017-10-6 19:49
重新修改了下,你再试试

谢谢你大神

13648893846

yuanziyou 发表于 2017-10-6 19:49
重新修改了下,你再试试

还有另一个BUG,当待求点垂直于起点时还是错误