如何将 Cartopy 等值线图更改为箱线图

Question

我绘制了一个 CartoPy 等高线图，如下所示：

使用以下脚本：

precip_full1 = xr.open_dataset('era_yr1979.nc')
precip_full2 = xr.open_dataset('era_yr1980.nc')
precip_full3 = xr.open_dataset('era_yr1981.nc')
precip_full4 = xr.open_dataset('era_yr1982.nc')
precip_full5 = xr.open_dataset('era_yr1983.nc')
precip_full6 = xr.open_dataset('era_yr1984.nc')

precip_full = xr.concat([precip_full1,precip_full2,precip_full3,precip_full4,precip_full5,precip_full6],dim = 'time')


output = []

for x in np.arange(6.5,10.25,0.25):
    for y in np.arange(-15,-9.75,0.25):
        precip = precip_full.where((precip_full.latitude==x)&(precip_full.longitude==y),drop=True)
        roll = precip.rolling(time=6,center=False).sum()
    

        annual = roll.groupby('time.year').max()

        tab = annual.to_dataframe().rename(columns={'tp':6})

    


output = pd.concat(output,1)


mean = output.mean()


data_mean = pd.DataFrame(mean, columns=['mean'])

df = data_mean.to_numpy()

new = [df[i:i+21] for i in range(0,len(df),21)]

new = np.reshape(new, [-1, 21])
df = pd.DataFrame(data=new, dtype=object)


lon2d, lat2d = np.meshgrid(lon, lat)


plt.figure(figsize=(6,5))
ax = plt.axes(projection=ccrs.PlateCarree())
ax.set_extent([-15,-10,6.5,10])
ax.coastlines()
ax.add_feature(cfeature.LAND)
ax.add_feature(cfeature.LAKES)
ax.add_feature(cfeature.RIVERS)
ax.add_feature(cfeature.BORDERS) 


gl = ax.gridlines(draw_labels=True, xlocs=np.arange(-180,180,0.25), ylocs=np.arange(-90,90,0.25),linewidth=0.4)
gl.top_labels   = False
gl.right_labels = False
plot = plt.contourf(lon2d, lat2d, df, cmap = 'jet', transform=ccrs.PlateCarree())

我现在意识到我更喜欢箱形图，每个网格框中都有一种纯色，我不再想要网格之间的插值点。

我发现我可以使用 pcolormesh 而不是轮廓来做到这一点。但是，当我更改最后一行代码时：

plot = plt.pcolormesh(lon2d, lat2d, df, cmap = 'jet', transform=ccrs.PlateCarree())

我收到以下错误：

TypeError: Dimensions of C (15, 21) are incompatible with X (15) and/or Y (15)

我看不到此错误意味着什么，也不知道如何修复它。有人做过类似的事情吗？

Answer 1

好的，既然您已经问过...这是一个快速示例，如何使用 eomaps 。 。

​
也可以通过m.new_layer_from_file.netcdf（...）

from eomaps import Maps
import numpy as np

# create some data in a regular lon/lat grid (=epsg 4326)
x, dx = np.linspace(-45, 45, 55, retstep=True)
y, dy = np.linspace(-20, 30, 25, retstep=True)
vals = np.random.randint(0,100, (x.size, y.size))

# plot the data as lon/lat rectangles on a map displayed in Orthographic projection.
m = Maps(Maps.CRS.Orthographic())
m.add_feature.preset.coastline()
m.set_data(vals, x, y, crs=4326)
m.set_shape.rectangles(radius=(dx/2, dy/2), radius_crs=4326)
m.plot_map()