matlab克里格插值(基于扇区邻域点选取)

前言

此次博客是前面mink函数的博客的进阶版本，往往我们在进行克里格插值时需要计算目标点周围邻域内多个点来进行插值操作，一般情况下邻域内的点数并不是越多越好，有时要根据邻域内点的分布分布状况以及点的分布个数来评估克里格插值的结果，但是有时邻域的点对数量选取是基于距离最短的原则来选取，这样就有时会造成邻域点选取的结果分布不均匀，为了使插值更加精确，arcgis中已经封装了设置领域点选择的一些条件，如下图所示，arcgis中还考虑了扇区的数目，每个扇区内点的数量限制，潘老师的理解是把扇区的数目就是把整个区域搜索的圆分成$n$等份，如下图会有更直观的理解：

从上图中定义的扇区大小为4，则将其分为4等分，每个扇区至少要有4个相邻要素，最多8个，因此权重那里最少是16个，那么问题来了，怎样在Matlab中编程将每个扇区内这样的点对提取出来呢。

1方法原理

1.1坐标系平移

对于坐标点集$(x_i,y_i)$,取出要搜索的目标点$(x_0,y_0)$,剩余点集为$(x_{0i},y_{0i})$,进行坐标平移，如下，新坐标系以$(x_0,y_0)$为原点，因此剩余点集坐标有：
$$\begin{gathered} x_{0i}=x_{0i}-x_0 \\ {y}_{0i}=y_{0i}-y_0 \end{gathered}$$
如下图所示：

1.2方位角确定

方位角原本是在测绘中的一个定义，这里我们求解的平移坐标系之后的剩余点集与原点连线与正北方向即$(y轴)$的夹角 α ∈ [ 0 , 2 π ] \al