利用GZXL程序：根据图面三角网查询任意点高程

yjycad

你好，提示（内插高程点; 错误: 除数为零），怎么处理啊？谢谢！！！

skg123

今天一朋友说有这种方法加密，算法不同，各位按需研究

反距离加权插值方法研究   
https://wenku.baidu.com/view/ae6 ... l?qq-pf-to=pcqq.c2c

树櫴希德

73哥  另类(作图法)三角网内插高程点

1. (defun PoInPl(pt lst / i p1 p2 an anl ret)
   
2.     (setq i -1 p1 (last lst))
   
3.     (while(and(not ret)(setq p2(nth(setq i(1+ i))lst)))
   
4.       (cond((equal p2 pt 1e-6)(setq ret t))
   
5.      (t(setq an(-(angle pt p1)(angle pt p2)))
   
6.       (if(equal pi(abs an) 1e-6)
   
7.         (setq ret t)
   
8.         (setq anl(cons(rem an PI)anl)))))
   
9.       (setq p1 p2))
   
10.     (cond(ret 0);线上;
    
11.    (t(if(equal PI(abs(apply'+ anl))1e-6)1 -1))))
    
12. 
    
13. (defun xyofen(e / q p i)
    
14.   (if(member(cdr(assoc 0(setq e(entget e))))'("LINE""TEXT""MTEXT""INSERT""SOLID"))
    
15.     (vl-remove'nil(mapcar'(lambda(x)(cdr(assoc x e)))'(10 11 12 13 14)))
    
16.     (progn(setq i -1 q(if(=(logand(cdr(assoc 70 e))1)1)0 1)e(cdr(assoc -1 e)))
    
17.       (repeat(fix(+(vlax-curve-getEndParam e)q))
    
18. (setq i(1+ i)q(vlax-curve-getPointAtParam e i)
    
19.       p(if(or(equal(last p)q 1e-5)(equal(car p)q 1e-5))p(cons q p))))
    
20.       (reverse p)
    
21.     )))
    
22. 
    
23. (defun c:tt(/ p e q q1 a);;另类(作图法)三角网内插高程点
    
24.   (while(setq p(getpoint))
    
25.     (if(and(setq e(ssget"F"(list'(-1e18 0)p)'((0 . "polyline")(100 . "AcDb3dPolyline"))))
    
26.    (setq e(car(vl-remove-if'(lambda(x)(=(poinpl p(setq pl(xyofen(cadr x))))-1))(ssnamex e)))));;p点所在三角网图元e
    
27.       (setq q(osnap(inters p(polar p(angle(mapcar'+(car pl)'(0 0))(cadr pl))(distance(mapcar'+(car pl)'(0 0))(cadr pl)))
    
28.    (mapcar'+(cadr pl)'(0 0))(mapcar'+(last pl)'(0 0))nil)"nea");;过p作1,2点的平行线相交于2，3所在边q
    
29.     q1(osnap(polar(cadr pl)(angle(cadr pl)(mapcar'+(car pl)'(0 0)))(distance(mapcar'+'(0 0)p)q))"nea");;在2 1边上截取pq长度于q1
    
30.     a(entmakex(list'(0 . "text")(cons 10 p)(cons 1(rtos(+(last q)(last q1)(-(last(cadr pl))));q点高程加上2，q1高差
    
31.        2 3))'(40 . 2)'(7 . "standard")))))))

(defun PoInPl(pt lst / i p1 p2 an anl ret)<br />     (setq i -1 p1 (last lst))<br />     (while(and(not ret)(setq p2(nth(setq i(1+ i))lst)))<br />       (cond((equal p2 pt 1e-6)(setq ret t))<br />      (t(setq an(-(angle pt p1)(angle pt p2)))<br />       (if(equal pi(abs an) 1e-6)<br />         (setq ret t)<br />         (setq anl(cons(rem an PI)anl)))))<br />       (setq p1 p2))<br />     (cond(ret 0);线上;<br />    (t(if(equal PI(abs(apply'+ anl))1e-6)1 -1))))<br /> <br /> (defun xyofen(e / q p i)<br />   (if(member(cdr(assoc 0(setq e(entget e))))'("LINE""TEXT""MTEXT""INSERT""SOLID"))<br />     (vl-remove'nil(mapcar'(lambda(x)(cdr(assoc x e)))'(10 11 12 13 14)))<br />     (progn(setq i -1 q(if(=(logand(cdr(assoc 70 e))1)1)0 1)e(cdr(assoc -1 e)))<br />       (repeat(fix(+(vlax-curve-getEndParam e)q))<br /> (setq i(1+ i)q(vlax-curve-getPointAtParam e i)<br />       p(if(or(equal(last p)q 1e-5)(equal(car p)q 1e-5))p(cons q p))))<br />       (reverse p)<br />     )))<br /> <br /> (defun c:tt(/ p e q q1 a);;另类(作图法)三角网内插高程点<br />   (while(setq p(getpoint))<br />     (if(and(setq e(ssget"F"(list'(-1e18 0)p)'((0 . "polyline")(100 . "AcDb3dPolyline"))))<br />    (setq e(car(vl-remove-if'(lambda(x)(=(poinpl p(setq pl(xyofen(cadr x))))-1))(ssnamex e)))));;p点所在三角网图元e<br />       (setq q(osnap(inters p(polar p(angle(mapcar'+(car pl)'(0 0))(cadr pl))(distance(mapcar'+(car pl)'(0 0))(cadr pl)))<br />    (mapcar'+(cadr pl)'(0 0))(mapcar'+(last pl)'(0 0))nil)"nea");;过p作1,2点的平行线相交于2，3所在边q<br />     q1(osnap(polar(cadr pl)(angle(cadr pl)(mapcar'+(car pl)'(0 0)))(distance(mapcar'+'(0 0)p)q))"nea");;在2 1边上截取pq长度于q1<br />     a(entmakex(list'(0 . "text")(cons 10 p)(cons 1(rtos(+(last q)(last q1)(-(last(cadr pl))));q点高程加上2，q1高差<br />        2 3))'(40 . 2)'(7 . "standard")))))))

zxvy

树櫴希德 发表于 2023-5-26 15:53
73哥  另类(作图法)三角网内插高程点

高程展点，对齐到中线，要是考虑到三角网，有何更好的思路？