Python plot 使用面

Question

我有這個文件： https://pastebin.com/WK6K97jv （列分別是 x,y,z=f(x,y)）。

我想 plot 使用 matplotlib (pyplot) 不幸的是我永遠不必處理 3d plot 我試過這樣做：

import matplotlib.pyplot as plt
import numpy as np
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter

fig = plt.figure()
ax = fig.gca(projection='3d')

X=np.genfromtxt('./jpdf1/jpdf1.000533048.dat',usecols=(0,))
Y=np.genfromtxt('./jpdf1/jpdf1.000533048.dat',usecols=(1,))
X, Y = np.meshgrid(X, Y)
Z=np.genfromtxt('./jpdf1/jpdf1.000533048.dat',usecols=(2,))

surf = ax.plot_surface(X, Y, Z, cmap=cm.coolwarm,
                       linewidth=0, antialiased=False)

但是口譯員支持我：

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/opt/homebrew/lib/python3.9/site-packages/matplotlib/_api/deprecation.py", line 431, in wrapper
    return func(*inner_args, **inner_kwargs)
  File "/opt/homebrew/lib/python3.9/site-packages/mpl_toolkits/mplot3d/axes3d.py", line 1658, in plot_surface
    raise ValueError("Argument Z must be 2-dimensional.")
ValueError: Argument Z must be 2-dimensional.

Answer 1

檢查Z.shape是否等於X.shape或Y.shape 。

Z.shape == X.shape

然后舉個例子。

# import package
import matplotlib.pyplot as plt
import numpy as np
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter

fig = plt.figure()
ax = fig.gca(projection='3d')


# generate x and y data
X = np.arange(20) # shape: 20
Y = np.arange(30) # shape:30

X, Y = np.meshgrid(X, Y) # convert x y data shape to new shape

print(X.shape) # (20,30)
print(Y.shape) # (20, 30)

Z= np.cos(X**2 + Y**2) # shape =(20,30)
printt(Z.shape) # (20, 30)

surf = ax.plot_surface(X, Y, Z, cmap=cm.coolwarm,
                       linewidth=0, antialiased=False)