python - 大型不规则网格到规则网格的二维插值

Question

我有 2048x2048 的不规则数据网格zi = f(xi, yi)，它们本质上是三组独立的 2048 个实数值。我需要将其平滑地插值(可能是双三次样条)到 wi = f(ui, vi) 的常规网格中，其中 ui 和 vi 是从 0 到 2047 的整数值。

我已经尝试过 griddata，它似乎在小于 1000x1000 的图像上工作良好，但当你达到 1500x1500 时就会爆炸(显然是 Delaunay 网格的内存 qhull 错误)。我看过一些 ndimage 函数，即 geometric_transform、RectBivariateSpline 和 map_coordinates，但它们似乎都需要正则化数据作为输入。我也可能遗漏了一些东西并且只是实现错误!

我正在尝试使用 Python/SciPy 来执行我一直在使用 tformarray 和 makeresampler 执行的这个 Matlab 脚本。关于我可以使用什么功能来处理这个大数据集有什么建议吗？谢谢!

Answer 1

我尝试重现您的错误，但没有成功。你是32位系统吗？我在使用 scipy/numpy 和大型数组时遇到问题，所以切换到 64 位，从那以后就没有问题了。

这是我用来尝试重现您的错误的代码(它不会生成任何有用的东西，但至少应该遇到相同的错误):

y,x=indices([2048,2048],dtype='float64')
z = randn(2048,2048)
yr = y + randn(2048,2048)
xr = x + randn(2048,2048)
zn = griddata(xr.ravel(),yr.ravel(),z.ravel(),x,y)
zl = griddata(xr.ravel(),yr.ravel(),z.ravel(),x,y,interp='linear')

这适用于我的机器。

如果您不能运行 64 位版本的 python(这可能会很困难，具体取决于您使用的操作系统)，您能否将 2048x2048 网格分成 4 个 1024x1024 网格？