将径向数据转换为笛卡尔网格以绘制曲面图

Question

我有一个数据列表，其分量对应于二维网格上某个径向距离 r 处的电位。数据对应极坐标中的数据点，在theta分量上是对称的

import numpy as np
import matplotlib.pyplot as plt

r = np.linspace(1.0, 5.0, 99)
#Data looks like:
V = np.array([9.0,...,0.0])

x = np.linspace(-5.0,5.0,99)
y = np.linspace(-5.0,5.0,99)

xx,yy = np.meshgrid(x,y)

我想在 (x,y) 空间中创建数据的表面图，但是要使用 matplotlib，您需要一个与每个 (x,y) 位置的电位相对应的数据点网格。鉴于我有一组在 (r,theta) 空间中测量的数据，我该如何创建曲面图？

Answer 1

您需要先将数据转换为笛卡尔坐标，然后通过为每个 theta 重复值 V 创建一个与网格网格形状相同的值数组。

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

fig = plt.figure()

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

r = np.linspace(1.0, 5.0, 99)
V = np.sqrt(np.sinc(r-0.5)**2) # your data here (same length as r)


theta = np.linspace(0,2*np.pi,50)
R, Theta = np.meshgrid(r,theta)
X = R * np.cos(Theta)
Y = R * np.sin(Theta)
Z = np.tile(V,(len(theta),1))

norm = plt.Normalize(vmin=Z.min(), vmax=Z.max())
ax.plot_surface(X,Y,Z, facecolors=plt.cm.RdYlGn(norm(Z)))

plt.show()