请问怎样将shp图置于顶层。雷达反射率置于低层呢？方便查看强回波具体位置

lixiliang2008

from pycwr.io import read_auto
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from pycwr.draw.RadarPlot import GraphMap
import cartopy.io.shapereader as shpreader
import cartopy.feature as cfeature  

plt.rcParams['font.sans-serif']=['SimHei']#用来正常显示中文
plt.rcParams['axes.unicode_minus']=False#用来正常显示负号

file = r"E:/Z_RADR_I_Z9311_20220727150001_O_DOR_SAD_CAP.bin.bz2"
NRadar = read_auto(file)
PyartRadar = NRadar.ToPyartRadar()


graph = GraphMap(NRadar, ccrs.PlateCarree())
ax = plt.axes(projection=ccrs.PlateCarree(), zorder = 0)
graph.plot_ppi_map(ax, 0, "dBZ", cmap="pyart_NWSRef")
ax.set_title("example of PPI with map", fontsize=16)
shpname = r'E:/shp/dijishi_2004.shp'
reader = shpreader.Reader(shpname)
states_provinces = cfeature.ShapelyFeature(reader.geometries(), crs = ccrs.PlateCarree(), edgecolor = 'k', facecolor = 'None')
ax.add_feature(states_provinces, zorder = 1)
plt.show()

lixiliang2008

已解决，将zorder设置的大一点。因为雷达图有好几层仰角！

任性小王子

图的色调不错哟

枝桠桠

请问方便分享下脚本文件吗？

muggle

请教楼主，graph.plot_ppi_map(ax, 0, "dBZ", cmap="pyart_NWSRef")这一行实现的是“第0个”仰角的反射率对吗？组合反射率怎么画呢？