如何使用scipy.interpolate中的griddata

Question

我有一個三列（x像素，y像素，z值）數據，有一百萬行。 數據來自圖像，並且存在重復的z值。 現在我需要制作表面圖。 這個圖像是一個很好的例子。 但現在輸出圖像為空。 有人可以查一下代碼嗎？

import numpy as np
from enthought.mayavi import mlab
from scipy.interpolate import griddata
x,y,z = np.loadtxt('test.csv',delimiter=',',usecols=(0,1,2),unpack=True) 
xi,yi = np.mgrid[0:3000:3000j, 0:3000:3000j]
zi = griddata((x, y), z, (xi, yi),method='linear')
mlab.surf(xi,yi,zi)
mlab.show()

Answer 1

我不能在沒有數據的情況下檢查代碼，但我懷疑問題是你使用默認的fill_value=nan作為griddata參數，所以如果你有網格點超出（x，y）的空間）點，網格中有NaNs，mlab可能無法處理（matplotlib不容易）。 嘗試設置fill_value=0或其他合適的實數。