ArcEngine 反距离权重插值法计算简单地面高程

手上有一张测绘的高程点的矢量图，需要做一个简单的高程插值算法，在鼠标移动是根据其周边的点的高程计算出当前的高程值，

不想做的太复杂，基本思路是按一定的半径搜索高程点，根据每一个高程点的距离远近进行加权平均，计算出高程值，

精确一点，就是使用反距离权重插值法，距离越远的高程点，影响度越低，如下图：

 

在ArcEngine中实现，以鼠标点为中心搜索 指定半径的图元，使用如下方法：

/// <summary> /// 搜寻高程点 /// </summary> /// <param name="elevationLayer">高程图层</param> /// <param name="currentPosition">当前鼠标位置</param> /// <param name="searchRadius">搜寻半径</param> /// <returns></returns> private IFeatureCache SearchPoint(ILayer elevationLayer, IPoint currentPosition, int searchRadius) { IFeatureLayer ftrLyr = elevationLayer as IFeatureLayer; IFeatureCache m_featureCache = new ESRI.ArcGIS.Carto.FeatureCacheClass(); m_featureCache.Initialize(currentPosition, searchRadius); m_featureCache.AddFeatures(ftrLyr.FeatureClass); return m_featureCache; }

 

然后遍历图元，使用IFeature的get_value方法获取图元的高程值字段内容，按算法计算一下，就出来了，凑合可以用了，至少比直接平均要来的好一点，使用的搜索半径默认为90M，这个是参考了国际DEM的精度，以及地图高程点的疏密程度来设置的，可根据具体情况设置。

 

参考资料：《地理信息系统算法基础》科学出版社 2006  p166