python使用双线性插值计算网格内数据

最近在处理再分析资料网格数据时，需要用到插值方法，双线性插值速度比较快，所以考虑双线性插值方法。

双线性插值，故名思义就是在x和y方向分别进行线性插值。双线性插值一般用于图像处理，对一个个像素块进行线性插值。

1.双线性插值

已知网格点Q12，Q22，Q11，Q21，但是要插值的点为P点，这就要用双线性插值了，首先在x轴方向上，对R1和R2两个点进行插值，这个很简单，然后根据R1和R2对P点进行插值，这就是所谓的双线性插值。在数学上，双线性插值是有两个变量的插值函数的线性插值扩展，其核心思想是在两个方向分别进行一次线性插值。

假如我们想得到未知函数 f 在点 P = (x, y) 的值，假设我们已知函数 f 在 Q11 = (x1, y1)、Q12 = (x1, y2), Q21 = (x2, y1) 以及 Q22 = (x2, y2) 四个点的值。

1） 在x方向上进行线性插值


2） 在y方向上进行线性插值

3） 综合起来，就是双线性插值的结果

其中(x2 - x1)， (y2 - y1) 为网格间距。这样只需要根据网格点上的值，网格间距，以及网格内需要插值的点的位置就可以计算插值结果。

2.python实现双线性插值

一般使用的再分析资料如NCEP，分辨率为1° x 1°的，网格间距即为1，则(x2 - x1)， (y2 - y1)都为1，

则上面的公式分母为1，结果为：
首