17.30. 插补和等高线绘制

Paolo Cavallini贡献的模块- Faunalia

备注

本章介绍如何使用不同的后端来计算不同的插补。

17.30.1. 内插

该项目显示了降雨量的梯度，从南到北。让我们使用不同的方法进行内插，所有方法都基于向量 points.shp ，参数 RAIN ：

警告

将单元格大小设置为 500 用于所有分析。

• GRASS ► v.surf.rst

• SAGA ► Multilevel B-Spline Interpolation

• SAGA ► Inverse Distance Weighted [Inverse distance to a power; Power: 4; Search radius: Global; Search range: all points]

• GDAL ► Grid (Inverse Distance to a power) [Power:4]

• GDAL ► Grid (Moving average) [Radius1&2: 50000]

然后测量不同方法之间的差异，并将其与到点的距离相关联：

• GRASS ► r.series [Unselect Propagate NULLs, Aggregate operation: stddev]

• GRASS ► v.to.rast.value on points.shp

• GDAL ► Proximity

• GRASS ► r.covar 显示关联矩阵；检查关联的重要性，例如用http://vassarstats.net/rsig.html.

因此，远离点的区域的插补精度会较低。

17.30.2. 等高线

绘制轮廓线的各种方法 [always step= 10] 在 stddev 栅格：

• GRASS ► r.contour.step

• GDAL ► Contour

• SAGA ► Contour lines from grid [NB: in some older SAGA versions, output shp is not valid, known bug]