为Cartopy地图中的缺失数据指定非透明颜色

Question

我在将Cartopy与matplotlib的pcolormesh一起使用以制作地图时，为指定缺失值的颜色时遇到麻烦。我想将所有“ nan”值设为灰色，而不是白色/透明。为此，我尝试使用matplotlib颜色图的set_bad方法，但这似乎破坏了Cartopy。该问题触发了Cartopy的UserWarning，但我不知道解决方法是什么。

这是一个甚至没有nan值但显示问题的示例：

import xarray as xr
import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import matplotlib as mpl
import cartopy
import sys
print(f"PYTHON: {sys.version}")  # PYTHON: 3.8.1 | packaged by conda-forge | (default, Jan 29 2020, 15:06:10) [Clang 9.0.1 ]
print(f" xarray {xr.__version__}")  # xarray 0.14.1
print(f" numpy {np.__version__}")  # numpy 1.17.3
print(f" matplotlib {mpl.__version__}")  # matplotlib 3.1.2
print(f" cartopy {cartopy.__version__}")  # cartopy 0.17.0

# construct fake data:
sig = np.sin(np.radians(np.linspace(-90, 90, 192)))
sig = np.tile(sig, [288,1]).T
darr = xr.DataArray(sig, dims=("lat", "lon"), coords={"lat":np.linspace(-90, 90, 192), "lon":np.linspace(0,360,288)})
# darr.plot.pcolormesh()  # Looks fine.

# fig, ax = plt.subplots(subplot_kw={"projection":ccrs.PlateCarree()})
# lons, lats = np.meshgrid(darr['lon'], darr['lat'])
# ax.pcolormesh(lons, lats, darr, transform=ccrs.PlateCarree())  # Looks fine.

cm = plt.get_cmap("viridis")
cm.set_bad("gray")

fig, ax = plt.subplots(subplot_kw={"projection":ccrs.PlateCarree()})
lons, lats = np.meshgrid(darr['lon'], darr['lat'])
ax.pcolormesh(lons, lats, darr, transform=ccrs.PlateCarree())  # Right side of map is gray / userwarning triggered

用户警告是：“用户警告：已经设置了颜色图的'坏'，但是为了将pcolormesh包裹在地图上，它必须是完全透明的。”

因此警告很明确，结果是一张地图，其中一半的数据被遮盖。我的问题是如何解决Cartopy的这种局限性？

Answer 1

尝试遇到问题的代码行

cm.set_bad("gray")

与

cm.set_bad( "gray", alpha=0 )

代表“完全透明”的alpha参数应该可以解决您的问题。