pcolormesh с отсутствующими значениями?

У меня есть 3 1-D ndarrays: x, y, z

и следующий код:

import numpy as np
import matplotlib.pyplot as plt
import scipy.interpolate as spinterp

## define data
npoints = 50
xreg = np.linspace(x.min(),x.max(),npoints)
yreg = np.linspace(y.min(),y.max(),npoints)
X,Y = np.meshgrid(xreg,yreg)
Z = spinterp.griddata(np.vstack((x,y)).T,z,(X,Y),
                      method='linear').reshape(X.shape)

## plot
plt.close()
ax = plt.axes()
col = ax.pcolormesh(X,Y,Z.T)
plt.draw()

Мой график выходит пустым, и я подозреваю, что это потому, что интерполяция method = 'linear' выходит с nans. I ' Я пробовал преобразовать в замаскированный массив, но безрезультатно - график все еще пустой. Можете ли вы сказать мне, что я делаю не так? Спасибо.