如何使用多色热图解决 MemoryError

我试图通过从文件中读取数据来绘制具有多种颜色的热图。我可以很好地生成 2D 和法线热图，但无法绘制如附图所示的图像。当使用随机数时，我可以绘制此图，但在从文件中读取数据时，它显示错误。

上面的热图是用随机数据生成的

Input: col[1] and col[2]是x and y坐标

00022d9064bc 819251 440006 1073260801 1073260803 2.0 
00022dba8f51 819251 440006 1073260801 1073260803 2.0 
00022de1c6c1 819251 440006 1073260801 1073260803 2.0 
003065f30f37 819251 440006 1073260801 1073260803 2.0 
00904b48a3b6 819251 440006 1073260801 1073260803 2.0 
00904b83a0ea 819213 439954 1073260803 1073260810 7.0 
00904b85d3cf 817526 439458 1073260803 1073261920 1117.0 
00904b14b494 817558 439525 1073260804 1073265410 4606.0 
00904b99499c 817558 439525 1073260804 1073262625 1821.0 
00904bb96e83 817558 439525 1073260804 1073265163 4359.0 
00904bf91b75 817558 439525 1073260804 1073263786 2982.0 
00022d36a6df 820428 438735 1073260807 1073260809 2.0 

Code:

from matplotlib import pyplot as plt 
from matplotlib import cm as CM
from matplotlib import mlab as ml
import numpy as np 

data = np.loadtxt('inputfile', unpack=True, dtype='str, int, int, int, int, float')

x  = data[1]
y  = data[2]

X, Y = np.meshgrid(x,y)

x = X.ravel()
y = Y.ravel()

gridsize = 30 
plt.subplot(111)

cb = plt.colorbar()
cb.set_label('density')
plt.show() 

Error:

Traceback (most recent call last):
  File "heat3.py", line 11, in <module>
    X, Y = np.meshgrid(x,y)
  File "/usr/lib/python2.7/dist-packages/numpy/lib/function_base.py", line 3106, in meshgrid
    mult_fact = np.ones(shape, dtype=int)
  File "/usr/lib/python2.7/dist-packages/numpy/core/numeric.py", line 178, in ones
    a = empty(shape, dtype, order)
MemoryError

任何有用的建议表示赞赏

您看到的错误来自meshgrid试图建立可能相当庞大的矩阵。如果您的数据包含N行，矩阵将为 N×N。根据您拥有的点数量以及它们的密集程度，您将需要热图中的两件事之一。任何一个

1. 您想要在相距较远的点之间进行插值以形成光滑的表面，或者
2. 您希望通过计算特定区域（2D 直方图）中有多少个点来聚合密集点。

我修改了下面的代码以适用于第二种情况（hexbin自动执行此操作），因为您似乎没有引用数据中的第三个值来进行插值。

from matplotlib import pyplot as plt 
from matplotlib import cm as CM
from matplotlib import mlab as ml
import numpy as np 

data = np.loadtxt('inputfile', unpack=True, dtype='str, int, int, int, int, float')

x = data[1]
y = data[2]
z = data[5]

# These lines are completely unnecessary and perhaps come
# from a different solution which was interpolating between points
#X, Y = np.meshgrid(x,y)
#x = X.ravel()
#y = Y.ravel()

gridsize = 30 
#plt.subplot(111)  # <- You don't need this as it is one plot anyway
plt.hexbin(x, y, C=z)   # <- You need to do the hexbin plot
cb = plt.colorbar()
cb.set_label('density')
plt.show() 

残余的meshgrid我上面注释掉的调用可能来自您在某处找到的一段代码，该代码执行第一个选项（在间隔点之间插值），也许我使用griddata http://docs.scipy.org/doc/scipy-0.15.1/reference/generated/scipy.interpolate.griddata.html。如果这确实是您想要的，您可以看看这个食谱条目 http://wiki.scipy.org/Cookbook/Matplotlib/Gridding_irregularly_spaced_data对不规则间隔的数据进行网格化