规则格网记录点号，反距离加权插值求出高程

馨馨

离散点格网化，记录点号在哪个格网，用反距离加权插值求出高程
把每个格网点当做原点，在四个象限内寻找4-8个离散点参与插值运算
设格网点与离散点距离为d，若d<20,范围内有4-8个离散点，则用这些离散点参与插值运算
若离散点小于4个或者大于8个，动态调整距离，是d增加或者减小，  代码是记录点号，问题是不能找出每个格网点参与运算的离散点，求大家帮忙

2. Sub RecordPointPosition()
   
3. Dim NX As Integer, NY As Integer '离散点所在网格的列'离散点所在的网格的行
   
4. Dim i As Integer, j As Integer
   
5. Dim Xmax As Double, Xmin As Double, Ymax As Double, Ymin As Double
   
6. 
   
7. dbld = 20
   
8. Open "E:\zuizhi.txt" For Input As #1 '打开输入文件
   
9. Input #1, Xmax, Ymax, Xmin, Ymin '读取文件数据
   
10. Close #1
    
11. 
    
12. L = 0
    
13. Open "E:\demdata.txt" For Input As #1 '打开原始文件
    
14. Do While Not EOF(1)
    
15. Input #1, H(L), X(L), Y(L), Z(L)
    
16. L = L + 1
    
17. Loop
    
18. Close #1
    
19. 
    
20. M = Fix((Xmax - Xmin) / dbld) + 1 '列数
    
21. N = Fix((Ymax - Ymin) / dbld) + 1 '行数
    
22. 
    
23. ReDim IP(L) As Long, IREG(N, M) As Double
    
24. 
    
25. 
    
26. '循环判断
    
27. For i = 1 To L - 1
    
28. '计算样点所在网格
    
29. NX = Fix((X(i) - Xmin) / dbld)
    
30. NY = Fix((Y(i) - Ymin) / dbld)
    
31. If IREG(NY, NX) = 0 Then
    
32.   IREG(NY, NX) = H(i) '记录点号
    
33. Else
    
34.   j = IREG(NY, NX)
    
35. '寻找链表的末端
    
36. Do
    
37.   If IP(j) = 0 Then
    
38.     IP(j) = H(i) '记录点号
    
39. Exit Do
    
40.   Else
    
41.     j = IP(j)
    
42. End If
    
43. Loop
    
44. 
    
45. End If
    
46. Next i
    
47. end sub

Sub RecordPointPosition()<br /> Dim NX As Integer, NY As Integer '离散点所在网格的列'离散点所在的网格的行<br /> Dim i As Integer, j As Integer<br /> Dim Xmax As Double, Xmin As Double, Ymax As Double, Ymin As Double<br /> <br /> dbld = 20<br /> Open "E:\zuizhi.txt" For Input As #1 '打开输入文件<br /> Input #1, Xmax, Ymax, Xmin, Ymin '读取文件数据<br /> Close #1<br /> <br /> L = 0<br /> Open "E:\demdata.txt" For Input As #1 '打开原始文件<br /> Do While Not EOF(1)<br /> Input #1, H(L), X(L), Y(L), Z(L)<br /> L = L + 1<br /> Loop<br /> Close #1<br /> <br /> M = Fix((Xmax - Xmin) / dbld) + 1 '列数<br /> N = Fix((Ymax - Ymin) / dbld) + 1 '行数<br /> <br /> ReDim IP(L) As Long, IREG(N, M) As Double<br /> <br /> <br /> '循环判断<br /> For i = 1 To L - 1<br /> '计算样点所在网格<br /> NX = Fix((X(i) - Xmin) / dbld)<br /> NY = Fix((Y(i) - Ymin) / dbld)<br /> If IREG(NY, NX) = 0 Then<br />   IREG(NY, NX) = H(i) '记录点号<br /> Else<br />   j = IREG(NY, NX)<br /> '寻找链表的末端<br /> Do<br />   If IP(j) = 0 Then<br />     IP(j) = H(i) '记录点号<br /> Exit Do<br />   Else<br />     j = IP(j)<br /> End If<br /> Loop<br /> <br /> End If<br /> Next i<br /> end sub