i = 0
pointsList=[]
x_list = []
y_list = []
z_list = []
for val in listSkv:
	pointsList.append([float(val.X/10000), float(val.Y/10000)])
	x_list.append(float(val.X))
	y_list.append(float(val.Y))
	z_list.append(float(listC01[i].h_oil/10))
	i += 1
	if i==201:
		break

x = np.asarray(x_list)
y = np.asarray(y_list)
z = np.asarray(z_list)
points = np.asanyarray(pointsList)

# new grid
x2, y2 = np.mgrid[0.1:0.9:201j, 0.1:0.9:201j]

# interpolate onto the new grid
z2 = griddata((x.ravel(), y.ravel()), z.ravel(), (x2, y2), method='cubic')
grid_z2 = griddata(points, z, (x2, y2), method='cubic')

print(z2)
print(grid_z2)

[[ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 ..., 
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]]
[[ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 ..., 
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]
 [ nan  nan  nan ...,  nan  nan  nan]]
[Finished in 7.5s]

i = 0
pointsList=[]
x_list = []
y_list = []
z_list = []
for val in listSkv:
	pointsList.append([float(val.X), float(val.Y)])
	x_list.append(float(val.X))
	y_list.append(float(val.Y))
	z_list.append(float(listC01[i].h_oil/10))
	i += 1
	if i==1000:
		break

x = np.asarray(x_list)
y = np.asarray(y_list)
values = np.asarray(z_list)

pointsList.sort(key=lambda x: x[0])
points = np.asanyarray(pointsList)

indent_x = (maxValX - minValX)*0.01
indent_y = (maxValY - minValY)*0.01
grid_x, grid_y = np.mgrid[minValX-indent_x:maxValX+indent_x:1000j, minValY-indent_y:maxValY+indent_y:1000j]

grid_z2 = griddata(points, values, (grid_x, grid_y), method='cubic')

plt.imshow(grid_z2.T, extent=(0,maxValX-minValX,0,maxValY-minValY), origin='lower')
plt.title('Cubic')

plt.show()