[编程申请]：求一个点到另两个点组成的直线的垂足

小菜

三个点均为空间点，怎么用lisp求一个点到另两个点组成的直线的垂足

陌生人

[php]
;返回结果：（世界坐标系的垂足  当前用户坐标系的垂足  以p1p2p3构成的坐标系的垂足）.
;测试:(setvar "pdmode" 3)(vl-cmdf ".point" (cadr (setq ptlst (3dper (getpoint)(getpoint)(getpoint)))))
(defun 3dper(p1 p2 p3 /   p3w pt ptw ptu)
  (setq p3w (trans p3 1 0))
  (vl-cmdf "ucs" "n" "3" p1 p2 p3)
  (setq pt (list (car (trans p3w 0 1)) 0 0)
        ptw (trans pt 1 0))
  (vl-cmdf "ucs" "p")
  (setq ptu (trans ptw 0 1))
  (list ptw ptu pt)
)
[/php]

lsjjm

试试这个:
(defun getperpt (p1 p2 p3)
  (setq ang (angle p2 p3))
  (setq p4 (inters p2 p3 p1 (polar p1 (+ ang (/ pi 2)) 1) nil))
  p4
)

小菜

这个不行，我原来也是这么写的，对于Z坐标不为0的点就不行了，比如P1的Z坐标为0，P2,P3的Z坐标不为0，返回nil，不正确
我想画一条P2,P3的直线，然后通过vlax-curve-getclosestpointto取得P1到这条直线的最近点，但这样好象太哆嗦了，有没有简单一点的算法？

eachy

最初由 小菜 发布
[B]这个不行，我原来也是这么写的，对于Z坐标不为0的点就不行了，比如P1的Z坐标为0，P2,P3的Z坐标不为0，返回nil，不正确
我想画一条P2,P3的直线，然后通过vlax-curve-getclosestpointto取得P1到这条直线的最近点，但?.. [/B]

将点转换为二维点。

aeo

p1 p2 p3

aX2 + k = Y2
aX3 + k = Y3
求K 和a
aX1 + k  --> Y1 -->p1

小菜

我数学没学好，五楼的公式没看懂
转换为二维点计算出的结果没有Z坐标，而且不是真正的垂足
我是这么写的，结果正确，但画条线太麻烦：

2. 
   
3.   [FONT=courier new]
   
4. (defun wllp1( c01 c02 c03 / c04 c05 c06) ;计算 c01到c02,c03线段的垂线的交点
   
5.   (command "line" c02 c03 "")
   
6.   (setq c04 (entlast)
   
7.         c05 (vlax-ename->vla-object c04)
   
8.         c06 (vlax-curve-getclosestpointto c05 (trans c01 1 0) acExtendBoth)
   
9.   )
   
10.   (entdel c04)
    
11.   (if c06 (setq c06 (trans c06 0 1)))
    
12. c06         
    
13. )  ;end of defun
    
14.   [/FONT]
    

<br />   [FONT=courier new]<br /> (defun wllp1( c01 c02 c03 / c04 c05 c06) ;计算 c01到c02,c03线段的垂线的交点<br />   (command "line" c02 c03 "")<br />   (setq c04 (entlast)<br />         c05 (vlax-ename->vla-object c04)<br />         c06 (vlax-curve-getclosestpointto c05 (trans c01 1 0) acExtendBoth)<br />   )<br />   (entdel c04)<br />   (if c06 (setq c06 (trans c06 0 1)))<br /> c06          <br /> )  ;end of defun<br />   [/FONT]<br />

wkai

最初由 小菜 发布
[B]转换为二维点计算出的结果没有Z坐标，而且不是真正的垂足
[B]

求出交点的投影后，可以再计算其Z坐标~

小菜

这样算出来的点不是垂足了吧？

ll_j

试试这样：

2. (defun 3dper(p1 p2 p3)
   
3.   (command "ucs" "3P" p1 p2 p3)
   
4.   (setq pd2 (trans p2 0 1)
   
5.         pd3 (trans p3 0 1)
   
6.         pd  (inters '(0 0 0)
   
7.                     (append (polar '(0 0) (+ (angle pd2 pd3) (/ pi 2)) 1) (list 0))
   
8.                     pd2
   
9.                     pd3
   
10.                     nil
    
11.             )
    
12.         pd  (trans pd 1 0)
    
13.   )
    
14.   (command "ucs" "")
    
15.   pd
    
16. )
    

(defun 3dper(p1 p2 p3)<br />   (command "ucs" "3P" p1 p2 p3)<br />   (setq pd2 (trans p2 0 1)<br />         pd3 (trans p3 0 1)<br />         pd  (inters '(0 0 0)<br />                     (append (polar '(0 0) (+ (angle pd2 pd3) (/ pi 2)) 1) (list 0))<br />                     pd2<br />                     pd3<br />                     nil<br />             )<br />         pd  (trans pd 1 0)<br />   )<br />   (command "ucs" "")<br />   pd<br /> )<br />

e2002

使用vl-curve类函数可以获得你要求的点

小菜

转换UCS的方法 是一个很好的想法，很简洁
我本想写个不用command命令的函数，但因为不知道怎么计算，不会写，就先用楼上的了，谢谢!