[日积月累]：计算地形图中某范围内的现状高程

eachy

2. 
   
3. 
   
4. 
   
5. 过三点的空间平面方程
   
6. 
   
7. | x -x1 y -y1 z -z1 |
   
8. | x2-x1 y2-y1 z2-z1 | = 0
   
9. | x3-x1 y3-y1 z3-z1 |
   
10. 
    
11. 根据拾取地形图中的高程数字计算三角形内点的高程;;
    
12. 
    
13. 声明:
    
14.      1 本程序需在WCS下执行;
    
15.      2 高程算法可能与其它软件计算的不同
    
16. 
    
17. |;
    
18. (defun c:test ( /    t1         t2   t3   elst1     elst2     elst3         h1
    
19.                h2   h3         p1   p3   p3        x1   x2          x3   y1   y2         y3
    
20.                P    a         b    c           d
    
21.               )
    
22.   (if (and (setq t1 (entsel "\nFirst text: "))
    
23.            (setq t2 (entsel "\nSencond text: "))
    
24.            (setq t3 (entsel "\nThird text: "))
    
25.            (setq p (getpoint "\nTest point: "))
    
26.       )
    
27.     (progn
    
28.       (setq elst1 (entget (car t1))
    
29.             elst2 (entget (car t2))
    
30.             elst3 (entget (car t3))
    
31.             h1          (atof (cdr (assoc 1 elst1)))
    
32.             h2          (atof (cdr (assoc 1 elst2)))
    
33.             h3          (atof (cdr (assoc 1 elst3)))
    
34.             p1          (trans (cdr (assoc 10 elst1)) 1 0)
    
35.             p2          (trans (cdr (assoc 10 elst2)) 1 0)
    
36.             p3          (trans (cdr (assoc 10 elst3)) 1 0)
    
37.             x1          (- (car p2) (car p1))       
    
38.             x2          (- (car p3) (car p1))       
    
39.             y1          (- (cadr p2) (cadr p1))
    
40.             y2          (- (cadr p3) (cadr p1))
    
41.             z1          (- h2 h1)               
    
42.             z2          (- h3 h1)               
    
43.             p (trans p 1 0)
    
44.       )
    
45.       (setq a (- (* y1 z2) (* y2 z1))
    
46.             b (- (* z1 x2) (* x1 z2))
    
47.             c (- (* y2 x1) (* y1 x2))
    
48.             d (- (+ (* (car p1) a)
    
49.                     (* (cadr p1) b)
    
50.                     (* (last p1) c)
    
51.                  )
    
52.               )
    
53.       )
    
54.     )
    
55.   )
    
56.   (- (/  (+ (* a (car p)) (* b (cadr p)) d) c))
    
57. )

<br /> <br /> <br /> 过三点的空间平面方程<br /> <br /> | x -x1 y -y1 z -z1 |<br /> | x2-x1 y2-y1 z2-z1 | = 0<br /> | x3-x1 y3-y1 z3-z1 |<br /> <br /> 根据拾取地形图中的高程数字计算三角形内点的高程;;<br /> <br /> 声明:<br />      1 本程序需在WCS下执行;<br />      2 高程算法可能与其它软件计算的不同<br /> <br /> |;<br /> (defun c:test ( /    t1         t2   t3   elst1     elst2     elst3         h1<br />                h2   h3         p1   p3   p3        x1   x2          x3   y1   y2         y3<br />                P    a         b    c           d<br />               )<br />   (if (and (setq t1 (entsel "\nFirst text: "))<br />            (setq t2 (entsel "\nSencond text: "))<br />            (setq t3 (entsel "\nThird text: "))<br />            (setq p (getpoint "\nTest point: "))<br />       )<br />     (progn<br />       (setq elst1 (entget (car t1))<br />             elst2 (entget (car t2))<br />             elst3 (entget (car t3))<br />             h1          (atof (cdr (assoc 1 elst1)))<br />             h2          (atof (cdr (assoc 1 elst2)))<br />             h3          (atof (cdr (assoc 1 elst3)))<br />             p1          (trans (cdr (assoc 10 elst1)) 1 0)<br />             p2          (trans (cdr (assoc 10 elst2)) 1 0)<br />             p3          (trans (cdr (assoc 10 elst3)) 1 0)<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 />             p (trans p 1 0)<br />       )<br />       (setq a (- (* y1 z2) (* y2 z1))<br />             b (- (* z1 x2) (* x1 z2))<br />             c (- (* y2 x1) (* y1 x2))<br />             d (- (+ (* (car p1) a)<br />                     (* (cadr p1) b)<br />                     (* (last p1) c)<br />                  )<br />               )<br />       )<br />     )<br />   )<br />   (- (/  (+ (* a (car p)) (* b (cadr p)) d) c))<br /> )