Python Matplotlib曲面图，内存错误。三列数据

我有一个大数据文件，包含x、y、z的3列（1,6mx3数组，128MB）。我试图在Matplotlib中生成一个轮廓图，但是当转换X，Y网格的数据时，我的内存用完了，这当然很可能发生。 下面是一个示例脚本（我还尝试了其他用于插值的numpy/scipy方法，但使用了相同的MemoryError）：

因此，我尝试使用MATLAB进行相同的程序，最后，要么我必须使用一个非常小的矩阵进行插值，从而丢失我感兴趣的图片（数据显示非常快的振荡）。下面是一个MATLAB脚本（变量“a”也表示数据数组）：

显然，它看起来与Python中的相同：） 为了证明得到我想要的东西是可行的，一个来自原点的结果是

我在网上找到的解决方案都不能解决我的问题

import numpy as np
import matplotlib.pyplot as plt    
a=np.loadtxt("data.dat")
xi = np.linspace(min(a[:,0]), max(a[:,0]), 40000)
yi = np.linspace(min(a[:,1]), max(a[:,1]), 40000)

X, Y = np.meshgrid(xi, yi)
Z = griddata(a[:,0], a[:,1], a[:,2], xi, yi, interp='nearest')
plt.contourf(X,Y,Z)
plt.show()

x1=a(:,1);
y1=a(:,2);
z1=a(:,3);
xi=linspace(min(x1),max(x1),40000);
yi=linspace(min(y1),max(y1),40000);
[XI YI]=meshgrid(xi,yi);
ZI = griddata(x1,y1,z1,XI,YI,'nearest');
contourf(XI,YI,ZI)