根据南方cass高程点 内插高程点

yanshengjiang

luzhmu 发表于 2025-9-29 08:29
三角网内插么？

根据greed鼠标位置创建矩形高程点选择集、得到三角网点表、判断在哪个三角形内、重心坐标插值法求点在面的投影高程
大概这几个步骤

yanshengjiang

测试过程中发现被广泛引用的zInsert函数有个bug，就是当三角形有一边在竖直方向上的时候，会出现除数为0的错误。需要先判断除数为0的时候 直接给ZL赋值

或者用下面这个deep seek代码代替这个功能，我简单测试这个函数是没有问题的，不知道到底对不对。

2. ;重心坐标插值法   代替双线性内插  避免垂足三角网发生除数为0的错误
   
3. (defun 重心坐标插值法求点在面的投影高程 (lst / p a b c denom lambda1 lambda2 lambda3 h)
   
4.   (setq p(car lst)
   
5.         a(car(cadr lst))
   
6.         b(cadr(cadr lst))
   
7.         c(caddr(cadr lst)))
   
8.   (setq denom (+ (* (- (cadr B) (cadr C)) (- (car A) (car C)))
   
9.                  (* (- (car C) (car B)) (- (cadr A) (cadr C)))))
   
10.   
    
11.   (setq lambda1 (/ (+ (* (- (cadr B) (cadr C)) (- (car P) (car C)))
    
12.                       (* (- (car C) (car B)) (- (cadr P) (cadr C))))
    
13.                  denom))
    
14.   
    
15.   (setq lambda2 (/ (+ (* (- (cadr C) (cadr A)) (- (car P) (car C)))
    
16.                       (* (- (car A) (car C)) (- (cadr P) (cadr C))))
    
17.                  denom))
    
18.   
    
19.   (setq lambda3 (- 1 lambda1 lambda2))
    
20.   
    
21.   ; 计算点P的高程
    
22.   (setq h(+ (* lambda1 (caddr A))
    
23.      (* lambda2 (caddr B))
    
24.      (* lambda3 (caddr C))))
    
25.   (list(car(car lst))(cadr(car lst))h)
    
26. )

;重心坐标插值法   代替双线性内插  避免垂足三角网发生除数为0的错误<br /> (defun 重心坐标插值法求点在面的投影高程 (lst / p a b c denom lambda1 lambda2 lambda3 h)<br />   (setq p(car lst)<br />         a(car(cadr lst))<br />         b(cadr(cadr lst))<br />         c(caddr(cadr lst)))<br />   (setq denom (+ (* (- (cadr B) (cadr C)) (- (car A) (car C))) <br />                  (* (- (car C) (car B)) (- (cadr A) (cadr C)))))<br />   <br />   (setq lambda1 (/ (+ (* (- (cadr B) (cadr C)) (- (car P) (car C))) <br />                       (* (- (car C) (car B)) (- (cadr P) (cadr C)))) <br />                  denom))<br />   <br />   (setq lambda2 (/ (+ (* (- (cadr C) (cadr A)) (- (car P) (car C))) <br />                       (* (- (car A) (car C)) (- (cadr P) (cadr C)))) <br />                  denom))<br />   <br />   (setq lambda3 (- 1 lambda1 lambda2))<br />   <br />   ; 计算点P的高程<br />   (setq h(+ (* lambda1 (caddr A)) <br />      (* lambda2 (caddr B)) <br />      (* lambda3 (caddr C))))<br />   (list(car(car lst))(cadr(car lst))h)<br /> )