带有彩色海洋的极地立体投影中的 Geopandas 世界地图

Question

对此添加进一步的要求，我还需要将海洋设置为蓝色（或任何其他颜色）。

对于 'PlateCarree' 投影，我可以简单地这样做

crs = ccrs.PlateCarree()
crs_proj4 = crs.proj4_init
world = gpd.read_file(gpd.datasets.get_path("naturalearth_lowres"))
w = world.to_crs(crs_proj4)
g = w.plot(facecolor='sandybrown', edgecolor='black')

现在添加海洋颜色

g.set_facecolor('#A8C5DD')

如果我现在想使用极立体投影

ccrs.NorthPolarStereo()

或

ccrs.SouthPolarStereo()

投影不起作用。当应用这个的答案时，我无法让海洋着色

Answer 1

您需要在 Cartopy geoaxes 上绘制地图几何图形，并使用 cartopy.feature.OCEAN 绘制海洋。这是您可以尝试的工作代码。阅读代码中的注释以进行说明。

import geopandas as gpd
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import cartopy

facecolor = 'sandybrown'
edgecolor = 'black'
ocean_color = '#A8C5DD'

#crs1 = ccrs.SouthPolarStereo()
crs1 = ccrs.NorthPolarStereo()

world = gpd.read_file(gpd.datasets.get_path("naturalearth_lowres"))
w1 = world.to_crs(crs1.proj4_init)

fig1, ax1 = plt.subplots(figsize=(7,7), subplot_kw={'projection': crs1})

# useful code to set map extent,
# --- if you want maximum extent, comment out the next line of code ---
ax1.set_extent([-60.14, 130.4, -13.12, -24.59], crs=ccrs.PlateCarree())

# at maximum extent, the circular bound trims map features nicely
ax1.add_geometries(w1['geometry'], crs=crs1, \
                facecolor=facecolor, \
                edgecolor=edgecolor, \
                linewidth=0.5)

# this adds the ocean coloring
ax1.add_feature(cartopy.feature.OCEAN, facecolor=ocean_color, edgecolor='none')

plt.show()

输出图将是：

带有彩色海洋的极地立体投影中的 Geopandas 世界地图

Geopandas world map in Polar Stereographic projection with coloured oceans

python

cartopy

geopandas