R语言中的空间插值

邻近多边形插值

这个方法可用于类别变量的插值。实际上是利用点生产沃诺尼多边形，每个点所在的沃诺尼多边形的值就等于点值，这样就实现了由点到面的变化，完成了插值。代码里的p是一个由很多点组成的空间对象，boundry即是我们的研究区域。p里的所有点应该在boundry的范围之内。在制作栅格、插值的过程中，都是基于p的外廓边界（bbox，这个中文是我随便起的，不知道对不对），即由所有坐标的极值围合形成矩形，所以为获得我们需要的研究边界范围内的插值数据还需要用boundry来进行裁剪。
1. sp对象版本

library(dismo)
pv <- voronoi(p)
ba <- aggregate(boundry)#boundry如果内部包含了许多更小的多边形，需要先进行合并
pv_b <- intersect(v, ba)

2.sf对象版本

sf对象完成同样的事情似乎要复杂的多，首先是对sfc对象的类型有着严格的要求，必须是MULTIPOINT类型,在st_voronoi函数里的对象也必须是sfc对象，返回的对象类型是GOMETRYCOLLECTION，完全不方便，所以其实更简单的方法应该是用sp的方法最后结果转成sf类型。纯sf类型的方法就比较繁琐了，以下是我摸索出来的。