Python中底图上的3D散点图

时间:2014-12-05 11:08:36

标签: python matplotlib-basemap

我正在尝试在python中绘制一个烟雾羽流的3D散点图,下面是北美地图。 我尝试了各种各样的方法,但似乎都有垮台。 我的代码目前是:

fig = plt.figure(figsize=(22,4))
ax = Axes3D(fig)

M = Basemap(llcrnrlon=-140,llcrnrlat=10,urcrnrlon=-40,urcrnrlat=70,
        projection='cyl', resolution ='c',area_thresh=100.)
ax.add_collection3d(M.drawcoastlines(linewidth=0.5))
ax.add_collection3d(M.drawcountries(linewidth=0.5))
ax.add_collection3d(M.drawparallels(np.arange(35,65,5),labels=[1,0,0,0]))

date = save_dir[-10:]
tau = tau_mod.int2tau(date)
dt = tau_mod.tau2time(tau)

q = ax.plot(max_x, max_y, max_z, c = 'black', lw = 3)

y,x=M(y,x)

p = ax.scatter(x,y,z,lw = 0,c = conc ,alpha = 0.75, s = 40,
                norm=matplotlib.colors.LogNorm(),vmin = 20,
               vmax = 4000, marker = ',',zorder=2)


ax.view_init(azim = a ,elev = e)


cb = fig.colorbar(p)
cb.set_label('CO [ppb]')
cmin , cmax = cb.get_clim()
ticks = np.linspace(cmin,cmax,10)
for t in range(len(ticks)):
    ticks[t] = int("%i" % ticks[t])
cb.set_ticks(ticks)
cb.set_ticklabels(ticks)

ax.set_xlabel('Longitude')
ax.set_ylabel('Latitude')
ax.set_zlabel('Altitude [km]')
plt.show()

我得到的两个主要问题是: 散点图不会出现在地图上 - 您可以在轴边界外的距离看到它。 drawparallels不起作用 - 它出现在这条消息中: 'T = ax.text(xlab,YY [n]的,latlab,的Horizo​​ntalAlignment =' 右”,verticalalignment = '中心',** kwargs) TypeError:text()至少需要5个参数(给定4个)'(我在2D绘图中使用此命令时从未遇到此问题)。

我还尝试使用ax.plot_surface绘制地图的png图像,但这会产生图像顺序的问题 - 无论我使用什么zorders,它似乎都会在散点图上绘制表面图像。

非常感谢任何帮助!

1 个答案:

答案 0 :(得分:0)

这看起来(至少部分)在分散函数中缺少latlon = True。

相关问题