Cartopy Contourf和Coastlines，但海岸线显示出来

Question

我正在尝试使用Cartopy的轮廓和海岸线绘制从GESC Earthdata档案中下载的区域的降水数据，但是由于某种原因，只有数据的CMAP才显示出来，并且数据被数据掩盖了海岸线。这是我的绘图代码。

plt.clf() #clear figure before
fig=plt.figure()

ax = plt.axes(projection=ccrs.PlateCarree(central_longitude=100.0, globe=None))
ax.set_extent([95,106,0,9]) # lon_left, lon_right, lat_below, lat_upper
ax.coastlines("10m", alpha=0.1) # avail:110m, 50m, 10m..... '10m' is better resolution than default
ax.gridlines(linewidths=0.01, draw_labels=True, alpha= 0.3)
ax.xlocator = mticker.FixedLocator(np.arange(95.,106.,0.5))
ax.ylocator = mticker.FixedLocator(np.arange(0.,9.,0.5))

x,y=np.meshgrid(lat,lon) 
clevs = np.arange(5,800,50)
cs = ax.contourf(y, x, finalmean,  levels=clevs, extend="max", cmap='viridis_r', transform=ccrs.PlateCarree())
# draw legend -------------
#cax = fig.add_axes([0.1,0.018,0.8,0.05], projection=ccrs.PlateCarree()) #plot cb at the bottom [left, bottom, width, height] 
#aa=fig.colorbar(cs,cax=cax,orientation='horizontal')
aa.set_label('Rainfall (mm)')
#plt.savefig("trmm_clim_mo_1.jpg", bbox_inches='tight')
plt.show()

这是我的输出

输出

我还指定了plt.contourf中的“变换”参数为受到遇到同样问题但无济于事的人的建议。如果有人可以指出我缺少的/错的代码，那将非常感谢。

Answer 1

尝试使用matplotlib kwarg zorder 彼此之间的层面。在您的示例中：

ax.coastlines（“ 10m”，alpha = 0.1，zorder = 3）

应该做。具有较低Zorder的对象具有较高的Zorder图。

或者，添加您的海岸线 呼叫轮廓。默认情况下，层添加的层将添加到前任之上。