Python 来自scipy.interpolate的bisplrep和bisplev问题

我从

scipy.interpolate

模块收到了一些关于

bisplrep

和

bisplev

的非常奇怪的行为。我尝试从scipy主页（）复制该示例。我有一个函数

f（x，y）

，它计算出一些z值，然后用

bisplrep

插值。如果我使用

bisplev

重新计算数据并绘制数据，结果值将旋转约90°。即使我交换x和y值，绘图仍会旋转。有人能告诉我我是不是做错了什么吗？以下代码应足以再现错误。我使用的是最新版本的scipy，错误发生在Jupyter、Spyder和IDLE中

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

def f(x, y):
    return x**2 + y**2

x, y = np.linspace(0, 5, 15), np.linspace(-2*np.pi, 2*np.pi, 15)
xx, yy = np.meshgrid(x, y)
zz = f(xx, yy)

tck = interpolate.bisplrep(xx, yy, zz)

plt.pcolor(xx, yy, zz)

x_new, y_new = np.linspace(0, 5, 100), np.linspace(-2*np.pi, 2*np.pi, 100)

z_new = interpolate.bisplev(x_new, y_new, tck)

plt.figure()    
plt.pcolor(x_new, y_new, z_new)

plt.figure()    
plt.pcolor(y_new, x_new, z_new)

plt.show()

---

我对网格使用了另一种定义。现在，它应该起作用了：

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

def f(x, y):
    return x**2 + y**2

#x, y = np.linspace(0, 5, 15), np.linspace(-2*np.pi, 2*np.pi, 15)
xx, yy = np.mgrid[0:5:15j, -2*np.pi:2*np.pi:15j]
zz = f(xx, yy)

tck = interpolate.bisplrep(xx, yy, zz)

plt.pcolor(xx, yy, zz)

#x_new, y_new = np.linspace(0, 5, 100), np.linspace(-2*np.pi, 2*np.pi, 100)
xx_new, yy_new = np.mgrid[0:5:100j, -2*np.pi:2*np.pi:100j]

zz_new = interpolate.bisplev(xx_new[:,0], yy_new[0,:], tck)

plt.figure()    
plt.pcolor(xx_new, yy_new, zz_new)

plt.figure()    
plt.pcolor(yy_new, xx_new, zz_new)

plt.show()

---

谢谢你的回答。我还认识到，当我转置

z_new

数组时（例如，如果我写入

z_new.T

），它会起作用。