pcolormesh（）和contourf（）不起作用

Question

尊敬的专家，我回想起两个月前我提出的一个问题，我一直在努力，因为没有成功。这涉及到底图上轮廓的叠加。我已经看过很多这方面的例子，例如：这里的示例：http://nbviewer.ipython.org/github/Unidata/tds-python-workshop/blob/master/matplotlib.ipynb 数据样本位于我以前的一篇帖子中：Contours with map overlay on irregular grid in python。 准备好数据后，这里有绘图方法：

# Setting the plot size and text
fig = plt.figure(figsize=(10,8))

lev = [15, 20, 25, 30, 35, 40,45]
norm1 = colors.BoundaryNorm(lev, 256)

# Draw filled contours

# 1. pcolor does not show the filled contours 
#cs = plt.pcolor(x,y,zi, cmap = cm.jet, norm = norm1)

# 2. pcolormesh does not show the filled contours
#cs = plt.pcolormesh(x,y,zi, shading = "flat", cmap=cmap)

# 3. contourf does not show the filled contours
#cs = plt.contourf(xi, yi, zi) #, levels=np.linspace(zi.min(),zi.max(),5))
cs = plt.contourf(xi, yi, zi, cmap = cm.jet, levels = lev, norm = norm1)

# 4. Draw line contours with contour()
#cs = m.contour(x,y,zi,linewidths=1.2)              # This works

plt.scatter(data.Lon, data.Lat, c=data.Z, s=100,
            vmin=zi.min(), vmax=zi.max())           # Does not work at all

# Color bar
#cbar = m.colorbar(fig,location='right',pad="10%")
fig.colorbar(cs)

# Plot a title
plt.figtext(.5,.05,'Figure 1. Mean Rainfall Onset Dates',fontsize=12,ha='center')

plt.show()

抱歉，我无法发布情节示例，但是：

上面的
• pcolor，pcolormesh和contourf都提供了没有任何填充轮廓但带有颜色条的地图
• 上面没有地图对象的图给出了填充轮廓，包括散点图（没有地图背景）
• 轮廓给出了叠加轮廓线的地图：

我感到困惑，因为这是从上面引用的链接中的示例中复制粘贴的示例。

任何有关问题可能原因的提示都将受到赞赏 Zilore Mumba

Answer 1

你需要使用底图来绘制轮廓与使用matplotlib.pyplot。看看我的一些代码示例。

#Set basemap and grid
px,py=n.meshgrid(x,y)

m=Basemap(projection='merc',llcrnrlat=20,urcrnrlat=55,
   llcrnrlon=230,urcrnrlon=305,resolution='l')

X,Y=m(px,py)


#Draw Latitude Lines
#labels[left,right,top,bottom]  1=True 0=False
parallels = n.arange(0.,90,10.)
m.drawparallels(parallels,labels=[1,0,0,0],fontsize=10,linewidth=0.)

# Draw Longitude Lines
#labels[left,right,top,bottom]  1=True 0=False
meridians = n.arange(180.,360.,10.)
m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10,linewidth=0)

#Draw Map
m.drawcoastlines()
m.drawcountries()
m.drawstates()
m.fillcontinents(color='grey',alpha=0.1,lake_color='aqua')

#Plot Contour lines and fill
levels=[5.0,5.1,5.2,5.3,5.4,5.6,5.7,5.8,5.9,6.0]
cs=m.contourf(px,py,thickness,levels,cmap=p.cm.RdBu,latlon=True,extend='both')
cs2=m.contour(px,py,thickness,levels,latlon=True,colors='k')

#Plot Streamlines
m.streamplot(px,py,U,V,latlon=True,color='k')

#Add Colorbar
cbar = p.colorbar(cs)
cbar.add_lines(cs2)
cbar.ax.set_ylabel('1000 hPa - 500 hPa Thickness (km)')

#Title
p.title('Geostrophic Winds with Geopotential Thickness')


p.show()

Answer 2

如果不知道您的数据是什么样子，那么回答您的问题有点困难，但无论如何我都会尝试。您可能希望对数据进行网格化，例如，使用直方图，然后对结果进行轮廓分析。

例如，如果您对绘制具有坐标（x，y）和第三个属性（z）的点的2D轮廓感兴趣，则需要使用对于颜色，你可以尝试一下

from numpy import *

H=histogram2d(x,y,weights=z)

contourf(H[0].T,origin='lower')

但是，就像我说的那样，如果您没有提供有关数据的详细信息，很难理解您正在寻找的内容。有关更多示例http://matplotlib.org/examples/pylab_examples/contourf_demo.html

，请查看matplotlib指南