Metpy: анализ снегопада строится, но проблема с цветовой картой ⇐ Python

[Anonymous]

График, который я построил для анализа снегопада, выглядит хорошо, но когда я рисую цветовую полосу, по какой-то причине он пропускает предпоследний цвет, и с помощью функции расширения = 'max' он отображает только последний цвет на палитре для последние два интервала на цветовой панели.
Я пробовал использовать cmap.set_over, но получаю те же результаты.
Набор данных: набор данных Snowfall< /p>

Код: Выделить всё

from datetime import datetime

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
import xarray as xr
import numpy as np
import metpy.calc as mpcalc
import netCDF4
from netCDF4 import Dataset
from metpy.units import units
from scipy.ndimage import gaussian_filter
import scipy.ndimage as ndimage
import matplotlib.colors as mcolors
from metpy.plots import USCOUNTIES
from herbie import Herbie, paint
from herbie.toolbox import EasyMap, pc

snowdata = xr.open_dataset('sfav2_CONUS_48h_2025010612.nc')

lat = snowdata['lat'][:].squeeze()
lon = snowdata['lon'][:].squeeze()
snow = snowdata['Data'][:].squeeze()
lon_2d, lat_2d = np.meshgrid(lon, lat)

snowin = snow*39.37

nws_snowfall_colors = [
"#ffffff",# 0 inches (white)
"#daeef5", # 0.1 inches
"#bdd7e7", # 1 inch
"#6baed6", # 2 inches
"#08519c", # 3 inches
"#082694", # 4 inches
"#ffff96", # 6 inches
"#ffc400", # 8 inches
"#ff8700", # 12 inches
"#db1400", # 18 inches
"#9e0000", # 24 inches
"#690000", # 30 inches
"#360000", # 36 inches
]
# Define the corresponding snowfall values (in inches)
nws_snowfall_cmap = mcolors.ListedColormap(nws_snowfall_colors)

clev_snowfall =  np.concatenate((np.arange(0, .1, .1), np.arange(.1, 1, .9), np.arange(1, 4, 1), np.arange(4, 8, 2), np.arange(8,12,4), np.arange(12,18,6), np.arange(18,24,6), np.arange(24,30,6), np.arange(30,36,6), np.arange(36,42,6)))
norm = mcolors.BoundaryNorm(clev_snowfall, 13)
#nws_snowfall_cmap.set_under("#ffffff")
#nws_snowfall_cmap.set_over("#360000")

datacrs = ccrs.PlateCarree()
plotcrs = ccrs.LambertConformal(central_latitude=35, central_longitude=-100,standard_parallels=(30, 60))
bounds = ([-110, -85, 33, 43])
fig = plt.figure(figsize=(14,12))
ax = fig.add_subplot(1,1,1, projection=plotcrs)
ax.set_extent(bounds, crs=ccrs.PlateCarree())
ax.add_feature(cfeature.COASTLINE.with_scale('50m'), linewidth = 0.75)
ax.add_feature(cfeature.STATES, linewidth = 2)
ax.add_feature(USCOUNTIES, edgecolor='black', linewidth = 1)
cf = ax.contourf(lon_2d, lat_2d, snowin, clev_snowfall, cmap = nws_snowfall_cmap, norm=norm, extend='max', transform=datacrs)
cb = plt.colorbar(cf, orientation='horizontal', ticks = [0, .1, 1, 2, 3, 4, 6, 8, 12, 18, 24, 30, 36])

colorbar: Colorbar
plot: график снегопада

Подробнее здесь: https://stackoverflow.com/questions/793 ... h-colormap