VBA写程序

馨馨

各位大神有会用反距离加权插值法建立规则格网的吗，我已经将样点格网化处理了，不知道怎么将样点记录在每个格网中，或者用链数组，不知道有会的没


该贴已经同步到 馨馨的微博

gzxl

没看明白，是要用距离加权平均法计算网格点的属性值吗？

馨馨

gzxl 发表于 2016-5-6 12:13
没看明白，是要用距离加权平均法计算网格点的属性值吗？

计算格网点的高程值，现在需要记录每个样点在哪个格网，不知道用什么方法，书上说设链数组，但是我看懂了不会写代码，感觉大神好像会，能不能给点思路提醒，那种方法容易些，谢谢啦

馨馨

有没有人啊，求指点

gzxl

馨馨 发表于 2016-5-6 22:27
有没有人啊，求指点

用离网格点最近的几个点进行计算
下面是c++的 用距离加权平均法计算网格点的属性值的代码

1. // 用距离加权平均法计算网格点的属性值
   
2. AcGePoint3dArray CCFGWDLG::CalculateValueByJL(const AcGePoint3dArray &gridPts, const AcGePoint3dArray &gcdPts)
   
3. {
   
4.     AcGePoint3dArray newgridPts;
   
5.     int kk = gcdPts.length();
   
6.     if (kk >= 10)
   
7.     {
   
8.         kk = 10;
   
9.     }
   
10.     for (int i = 0; i < gridPts.length(); i++)
    
11.     {
    
12.         AcGePoint2dArray distanceArray;
    
13.         AcGePoint3d pt1 = gridPts;
    
14.         for (int j = 0; j < gcdPts.length(); j++)
    
15.         {
    
16.             AcGePoint3d pt2 = gcdPts[j];
    
17.             double dist = GetDistance(pt1.x, pt1.y, pt2.x, pt2.y);
    
18.             AcGePoint2d temp((double)j, dist);
    
19.             distanceArray.append(temp);
    
20.         }
    
21. 
    
22.         SortPointsOnY(distanceArray, 0, distanceArray.length() - 1);
    
23.         double up = 0, down = 0, dz;
    
24. 
    
25.         for (int k = 0; k < kk; k++)
    
26.         {
    
27.             AcGePoint2d temp = distanceArray[k];
    
28.             int index = (int)temp.x;
    
29.             double ptz = gcdPts[index].z;
    
30.             double dis = distanceArray[index].y;
    
31.             up += ptz / dis;
    
32.             down += 1 / dis;
    
33.         }
    
34.         dz = up / down;
    
35.         AcGePoint3d pt3d(pt1.x, pt1.y, dz);
    
36.         newgridPts.append(pt3d);
    
37.     }
    
38.     return newgridPts;
    
39. }

// 用距离加权平均法计算网格点的属性值<br /> AcGePoint3dArray CCFGWDLG::CalculateValueByJL(const AcGePoint3dArray &gridPts, const AcGePoint3dArray &gcdPts)<br /> {<br />     AcGePoint3dArray newgridPts;<br />     int kk = gcdPts.length();<br />     if (kk >= 10)<br />     {<br />         kk = 10;<br />     }<br />     for (int i = 0; i < gridPts.length(); i++)<br />     {<br />         AcGePoint2dArray distanceArray;<br />         AcGePoint3d pt1 = gridPts<i>;<br />         for (int j = 0; j < gcdPts.length(); j++)<br />         {<br />             AcGePoint3d pt2 = gcdPts[j];<br />             double dist = GetDistance(pt1.x, pt1.y, pt2.x, pt2.y);<br />             AcGePoint2d temp((double)j, dist);<br />             distanceArray.append(temp);<br />         }<br /> <br />         SortPointsOnY(distanceArray, 0, distanceArray.length() - 1);<br />         double up = 0, down = 0, dz;<br /> <br />         for (int k = 0; k < kk; k++)<br />         {<br />             AcGePoint2d temp = distanceArray[k];<br />             int index = (int)temp.x;<br />             double ptz = gcdPts[index].z;<br />             double dis = distanceArray[index].y;<br />             up += ptz / dis;<br />             down += 1 / dis;<br />         }<br />         dz = up / down;<br />         AcGePoint3d pt3d(pt1.x, pt1.y, dz);<br />         newgridPts.append(pt3d);<br />     }<br />     return newgridPts;<br /> }

馨馨

gzxl 发表于 2016-5-6 23:40
用离网格点最近的几个点进行计算
下面是c++的 用距离加权平均法计算网格点的属性值的代码

谢谢你

馨馨

gzxl 发表于 2016-5-6 23:40
用离网格点最近的几个点进行计算
下面是c++的 用距离加权平均法计算网格点的属性值的代码

谢谢你

hounengwei

非常感谢楼主提供分享！