如何使地图线在 cartopy 中的数据上可见？

Question

我试图使地图线始终位于数据之上，但在 cartopy 中找不到合适的 command/options 来这样做。在下面的图中，我想要粗蓝线 "under" 黑色状态线但在米色状态之上。

代码：

import matplotlib.pyplot as plt
import numpy as np
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import os,sys

conus_proj = ccrs.LambertConformal(central_longitude=-96,central_latitude=39.0)
fig = plt.figure(figsize=(10,8))
ax = fig.add_subplot(1,1,1,projection=conus_proj)
ax.set_extent([-120,-70,22,50])
#ax.add_feature(cfeature.BORDERS)
ax.add_feature(cfeature.COASTLINE)
ax.add_feature(cfeature.OCEAN, facecolor='#CCFEFF')
ax.add_feature(cfeature.LAKES, facecolor='#CCFEFF')
ax.add_feature(cfeature.RIVERS, facecolor='#CCFEFF')
ax.add_feature(cfeature.LAND, facecolor='#FFE9B5')
state_borders = cfeature.NaturalEarthFeature(category='cultural', name='admin_1_states_provinces_lakes', scale='50m', facecolor='#FFE9B5')
ax.add_feature(state_borders, edgecolor='black')

plt.plot([-120,-70],[35,45],linewidth=8, transform=ccrs.PlateCarree())
plt.show()

我曾尝试更改 ax.add_feature(state_borders...) 和 plt.plot(...) 的 zorder，但收到了奇怪的结果。默认情况下，州边界在 pcolormesh 的顶部，但在 plt.plot

时似乎并非如此

以上代码的输出如下：

Answer 1

pcolormesh 与 plot 的结果不同的原因是它们具有不同的默认 zorders。如果我将州边界的 zorder（注意下面我使用 Cartopy 对州的内置支持）设置为 10，我会让它们出现在图的顶部：

import matplotlib.pyplot as plt
import numpy as np
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import os,sys

conus_proj = ccrs.LambertConformal(central_longitude=-96,central_latitude=39.0)
fig = plt.figure(figsize=(10,8))
ax = fig.add_subplot(1,1,1,projection=conus_proj)
ax.set_extent([-120,-70,22,50])

ax.add_feature(cfeature.COASTLINE)
ax.add_feature(cfeature.OCEAN, facecolor='#CCFEFF')
ax.add_feature(cfeature.LAKES, facecolor='#CCFEFF')
ax.add_feature(cfeature.RIVERS, edgecolor='#CCFEFF')
ax.add_feature(cfeature.LAND, facecolor='#FFE9B5')
ax.add_feature(cfeature.STATES, edgecolor='black', zorder=10)

plt.plot([-120,-70],[35,45],linewidth=8, transform=ccrs.PlateCarree())
plt.show()

我还必须删除某些特征的面部颜色。

如何使地图线在 cartopy 中的数据上可见？

how to make map lines visible over data in cartopy?

python

cartopy