XD::Geom:QueryElev 查询三角网指定位置的高程值

newer

1. ;|
   
2. 查询三角网中指定位置的高程
   
3. 参数： lst -- 3D FACE实体 or 3点表
   
4. 返回值：实数
   
5. |;
   
6. (defun XD::Geom:QueryElev (lst p / ge int l)
   
7.   (and
   
8.     (or
   
9.       (and
   
10.         (= (type lst) 'ENAME)
    
11.         (setq ge (xdge::constructor lst))
    
12.       )
    
13.       (and
    
14.         (= (type lst) 'LIST)
    
15.         (setq ge (xdge::constructor "kBoundedPlane" lst))
    
16.       )
    
17.     )
    
18.     (setq l (xdge::constructor "kLine3d" p (mapcar
    
19.                                              '+
    
20.                                              p
    
21.                                              (getvar "viewdir")
    
22.                                            )
    
23.             )
    
24.     )
    
25.     (setq int (xdge::getpropertyvalue ge "Intersectwith" l))
    
26.     (xdge::free)
    
27.   )
    
28.   (last int)
    
29. )
    
30. 
    

;|<br /> 查询三角网中指定位置的高程<br /> 参数： lst -- 3D FACE实体 or 3点表<br /> 返回值：实数<br /> |;<br /> (defun XD::Geom:QueryElev (lst p / ge int l)<br />   (and<br />     (or<br />       (and<br />         (= (type lst) 'ENAME)<br />         (setq ge (xdge::constructor lst))<br />       )<br />       (and<br />         (= (type lst) 'LIST)<br />         (setq ge (xdge::constructor "kBoundedPlane" lst))<br />       )<br />     )<br />     (setq l (xdge::constructor "kLine3d" p (mapcar<br />                                              '+<br />                                              p<br />                                              (getvar "viewdir")<br />                                            )<br />             )<br />     )<br />     (setq int (xdge::getpropertyvalue ge "Intersectwith" l))<br />     (xdge::free)<br />   )<br />   (last int)<br /> )<br /> <br />

st788796

贴个以前写的采用用平面方程计算

1. ;|
   
2. 过三点的空间平面方程 p p1 p2 p3 WCS点
   
3. | x -x1 y -y1 z -z1 |
   
4. | x2-x1 y2-y1 z2-z1 | = 0
   
5. | x3-x1 y3-y1 z3-z1 |
   
6. |;
   
7. (defun Geom:QueryElev (p p1 p2 p3 / x1 x2 y1 y2 z1 z2 a b c d h1 h2 h3)
   
8.   (setq        h1 (caddr p1)
   
9.         h2 (caddr p2)
   
10.         h3 (caddr p3)
    
11.         x1 (- (car p2) (car p1))
    
12.         x2 (- (car p3) (car p1))
    
13.         y1 (- (cadr p2) (cadr p1))
    
14.         y2 (- (cadr p3) (cadr p1))
    
15.         z1 (- h2 h1)
    
16.         z2 (- h3 h1)
    
17.   )
    
18.   (setq        a (- (* y1 z2) (* y2 z1))
    
19.         b (- (* z1 x2) (* x1 z2))
    
20.         c (- (* y2 x1) (* y1 x2))
    
21.   )
    
22.   (if (< c 0.)
    
23.     (setq a (- a)
    
24.           b (- b)
    
25.           c (- c)
    
26.     )
    
27.   )
    
28.   (setq        d (- (+        (* (car p1) a)
    
29.                 (* (cadr p1) b)
    
30.                 (* h1 c)
    
31.              )
    
32.           )
    
33.   )
    
34.   (- (/ (+ (* a (car p)) (* b (cadr p)) d) c))
    
35. )
    

;|<br /> 过三点的空间平面方程 p p1 p2 p3 WCS点<br /> | x -x1 y -y1 z -z1 |<br /> | x2-x1 y2-y1 z2-z1 | = 0<br /> | x3-x1 y3-y1 z3-z1 |<br /> |;<br /> (defun Geom:QueryElev (p p1 p2 p3 / x1 x2 y1 y2 z1 z2 a b c d h1 h2 h3)<br />   (setq        h1 (caddr p1)<br />         h2 (caddr p2)<br />         h3 (caddr p3)<br />         x1 (- (car p2) (car p1))<br />         x2 (- (car p3) (car p1))<br />         y1 (- (cadr p2) (cadr p1))<br />         y2 (- (cadr p3) (cadr p1))<br />         z1 (- h2 h1)<br />         z2 (- h3 h1)<br />   )<br />   (setq        a (- (* y1 z2) (* y2 z1))<br />         b (- (* z1 x2) (* x1 z2))<br />         c (- (* y2 x1) (* y1 x2))<br />   )<br />   (if (< c 0.)<br />     (setq a (- a)<br />           b (- b)<br />           c (- c)<br />     )<br />   )<br />   (setq        d (- (+        (* (car p1) a)<br />                 (* (cadr p1) b)<br />                 (* h1 c)<br />              )<br />           )<br />   )<br />   (- (/ (+ (* a (car p)) (* b (cadr p)) d) c))<br /> )<br />