小结空间插值分析

Posted by Xiaoyan(Sharon) Mu on 2017-10-15
Estimated Reading Time 6 Minutes
Words 1.8k In Total
Viewed Times

空间插值分析的基本原则是“地理学第一定律”,即“Everything is related to everything else, but near things are more related than distant things”,这是Waldo Tobler提出的。 例如,如果你所在的位置正在下雨,那么很有可能离你100米的位置也在下雨,而20公里以外的地方下雨就不一定了,或者说可能性就低了很多。 大多数插值方法遵循了这样的准则。

插值方法大致可以分为,确定性插值方法/Deterministic interpolation methods和随机性插值方法: 确定性插值方法使用了确定性算法,即,给予特定的输入数据我们会得到相同的结果; 而随机性插值方法是基于统计学的,通常用于更高级的表面预测建模,这种建模方法不仅可以预测表面的数值,也可能会带来误差或者不确定性,一些常见的地统计插值方法就属于这个范畴。

从参与计算预测值的采样点的范围来看,确定性插值方法可以分为全局/Globe局部/Local两种。 全局插值方法在计算预测值时使用了整个数据集,而局部插值方法是通过特定的领域来计算预测值。

从预测表面是否穿过已知采样点来看,确定性插值还可以分为精确/exact非精确/inexact两种。 如果插值算法预测的采样点位置的数值与采样点自身的数据一致,那么就是精确插值;否则为非精确插值。

  • 精确插值: 反距离权重法/IDW,径向基函数法/radial basis functions,样条函数法/Spline
  • 非精确插值: 全局或局部多项式法/global and local polynomial

IDW 插值方法

IDW是一种非常实用又很好理解的插值方法,适用于大多数情况。但是有关这个算法还有些注意事项:IDW并没有考虑数据的趋势,例如,如果数据在不同的方向上有不同的变化,那么这会被IDW方法忽略,IDW得到的是各个方向均一的预测插值;IDW插值仅在采样点的最大值和最小值之间进行预测插值,如果你是高程表面建模,那么山峰和河谷很有可能被抹平了。

IDW插值时,可以设置power参数来控制周围采样点对插值影响的显著性。 power值越大,较远的点对预测插值的影响就会越弱,也就是说,远处的采样点对预测值的影响随着power的增大而削减越快。power可以是任意正数,但是通常取值范围是0.5-3,默认为2。

有关IDW还有个常见小问题。既然IDW是精确插值,那为什么有时候看采样点所在像元的预测值与采样点值还是稍有差异?这是因为采样点并没有恰好落在像元的正中心。

样条函数插值方法/Spline Method

样条函数相比IDW更灵活,我们可以想象一下通过一张橡皮膜穿过所有已知采样点来预测表面。 样条函数法不像IDW,它可以是预测值高于或低于已知采样点的最大值或最小值。 然而,如果相聚非常近的采样点的值变化非常剧烈,样条函数就会出现异常,因为样条函数的算法中是基于坡度(单位距离的值变化量)来计算预测值的。 也就是说,例如,如果你的数据中如果含有悬崖,断层线,最好还是用IDW吧。

样条函数有两种类型: regularized 和 tension。 使用相同的采样点,tension spline比 regularized spline得到的预测表面更扁平,也就是使得预测值更接近采样值。 或者说tension spline更僵硬,而regularized spline 更有弹性,这个描述看起来更生动。

另外一个可以进一步控制预测表面的形态的参数是权重/weight。对于regularized spline,权重越高表面越平滑,常用的典型参数有0, 0.001、 0.01、 0.1、 0.5,默认是0.1;对于tension spline,权重越高表面约粗糙,常用的典型参数有0、 1、 5、10。

地形转拓扑/Topo to Raster

这是一种水文校正的插值方法,可以综合点线面输入数据,充当不同的角色或类别进行表面建模。

从帮助中摘过来的描述:

共有九种数据类型可作为地形类的输入:

  • PointElevation 表示表面高程的点要素类,field 用于存储点的高程。
  • Contour 表示高程等值线的线要素类, field 用于存储等值线的高程。
  • Stream 河流位置的线要素类。所有弧线必须定向为指向下游。要素类中应该仅包含单条弧线组成的河流。
  • Sink 表示已知地形凹陷的点要素类。
  • Boundary 边界是包含表示输出栅格外边界的单个面的要素类。在输出栅格中,位于此边界以外的像元将为 NoData。此选项可用于在创建最终输出栅格之前沿海岸线裁剪出水域。
  • Lake 指定湖泊位置的面要素类。湖面内的所有输出栅格像元均将指定为使用沿湖岸线所有像元高程值中最小的那个高程值。
  • Cliff 悬崖的线要素类。必须对悬崖线要素进行定向以使线的左侧位于悬崖的低侧,线的右侧位于悬崖的高侧。
  • Exclusion 其中的输入数据应被忽略的区域的面要素类。这些面允许从插值过程中移除高程数据。通常将其用于移除与堤壁和桥相关联的高程数据。这样就可以内插带有连续地形结构的基础山谷。
  • Coast 包含沿海地区轮廓的面要素类。位于这些面之外的最终输出栅格中的像元会被设置为小于用户所指定的最小高度限制的值。
    PointElevation、Contour 和 Sink 要素输入类型均可指定包含 z 值的字段。但对于 Boundary、Lake、Cliff、Coast、Exclusion 或 Stream 输入类型则不存在 Field 选项。

容差 1 和容差 2 设置的典型值为:

  • 对于比例为 1:100,000 的点数据,请分别使用 5.0 和 200.0。
  • 对于比例高达 1:500,000 的较不密集的点数据,请分别使用 10.0 和 400.0。
  • 对于等值线间距为 10 的等值线数据,请分别使用 5.0 和 100.0。
    容差 2 应至少比容差 1 大 6 倍。

插值方法其实还有很多,今天就针对Esri Training的课程整理这些吧。


如果您喜欢此博客或发现它对您有用,则欢迎对此发表评论。 也欢迎您共享此博客,以便更多人可以参与。 如果博客中使用的图像侵犯了您的版权,请与作者联系以将其删除。 谢谢 !