应用于阵列的双线性插值 (2D)

为了在 Python 中自己实现双线性插值，我读了很多书。目前，我有一个基于我在这里找到的函数 (How to perform bilinear interpolation in Python)。

主要问题是我必须将它应用于两个数组。例如，我有一个数组 X (len(X)=100) 和一个数组 Y (len(Y)=100)，我想要这些点中的值，所以我也有一个大小为 100 的输出 Z。

为此，我必须找到最近的 4 个点才能进行插值，这是我已经完成的。

因此，最明显的方法是使用循环并遍历应用该函数的 X 和 Y 值。

我的问题是是否有办法加快实现避免循环。所以我想要 Z (len(Z)=100) 直接没有循环，就像 Scipy 函数一样（我不能使用 Scipy）。

谢谢，如果需要，我会给出更多解释