快速遍历2D数组的方法

Question

我正在使用噪声库根据数组的索引生成2D Perlin噪声。

当前，我正在遍历x，y索引并为每个索引集调用函数。有没有一种更快的方法可以在numpy中做到这一点。

import numpy as np
import noise

def func(x, y, a, b, T):
    return noise.pnoise3(x/a, y/b, T)

我想做的就是将此函数的结果分配给2D数组的x，y索引。我目前正在做的是：

n = np.empty(a, b)

for ix, iy in np.ndindex((a, b)):
    n[ix, iy] = func(ix, iy, a, b, T)

预期结果是一个[a x b]数组，并在这些点上评估了噪声值。

在仍然使用噪声库的同时，有没有一种更快的方法？