Matplotlib 参数曲面图出乎意料的结果，为什么？

Question

我正在尝试绘制一个椭球体，所以我想我会修改 matplotlib 3D 绘图页面中的球体示例代码。

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
from matplotlib import cm
import numpy as np

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

# Ellipsoid
u = np.linspace(-np.pi/2.0,np.pi/2.0,100)
v = np.linspace(-np.pi,np.pi,100)
x = 10 * np.outer(np.cos(u), np.cos(v))
y = 10 * np.outer(np.cos(u), np.sin(v))
z = 10 * np.outer(np.ones(np.size(u)), np.sin(v))


# Sphere
#u = np.linspace(0, 2 * np.pi, 100)
#v = np.linspace(0, np.pi, 100)
#x = 10 * np.outer(np.cos(u), np.sin(v))
#y = 10 * np.outer(np.sin(u), np.sin(v))
#z = 10 * np.outer(np.ones(np.size(u)), np.cos(v))

ax.plot_surface(x, y, z,  rstride=4, cstride=4, cmap = cm.copper)
ax.set_xlabel('x-axis')
ax.set_ylabel('y-axis')
ax.set_zlabel('z-axis')
plt.show()

如果你运行代码你会看到情节 returns 一个美观的半里外外的船状表面但遗憾的是不是椭圆体。

已包含 sphere 代码（已注释掉）以供比较。

这里有什么我明显遗漏的地方吗？

Answer 1

为什么要更改参数化？以球体为例，只需改变半轴长度即可：

# Ellipsoid
u = np.linspace(0, 2.*np.pi, 100)
v = np.linspace(0, np.pi, 100)
x = 60 * np.outer(np.cos(u), np.sin(v))
y = 20 * np.outer(np.sin(u), np.sin(v))
z = 10 * np.outer(np.ones(np.size(u)), np.cos(v))

Matplotlib 参数曲面图出乎意料的结果，为什么？

Matplotlib parametric surface plot unexpected results, why?

python

3d

matplotlib

surface

parametric-equations