歌德的地形图形状和外部不透明度的应用

Question

我正在绘制一个城市边界(geopandas dataframe)，我使用上下文添加了一个基本地图。我想把不透明度应用于城市界限以外的地图区域。下面的例子显示了与预期效果相反的效果，因为除了城市界限之外，不透明度应该适用于任何地方。

import osmnx as ox
import geopandas as gpd
import contextily as cx

berlin = ox.geocode_to_gdf('Berlin,Germany')
fig, ax = plt.subplots(1, 1, figsize=(10,10))
_ = ax.axis('off')

berlin.plot(ax=ax,
            color='white',
            edgecolor='black',
            alpha=.7,
           )

# basemap 
cx.add_basemap(ax,crs=berlin.crs,)

plt.savefig('stackoverflow_question.png',
            dpi=100,
            bbox_inches='tight',
           )

显示与预期结果相反的地块：

​

Answer 1

您可以创建一个新的多边形，它是在几何图形减去几何图形的全部边界上的缓冲区。

import osmnx as ox
import geopandas as gpd
import contextily as cx
import matplotlib.pyplot as plt
from shapely.geometry import box


berlin = ox.geocode_to_gdf("Berlin,Germany")
notberlin = gpd.GeoSeries(
    [
        box(*box(*berlin.total_bounds).buffer(0.1).bounds).difference(
            berlin["geometry"].values[0]
        )
    ],
    crs=berlin.crs,
)


fig, ax = plt.subplots(1, 1, figsize=(10, 10))
_ = ax.axis("off")

notberlin.plot(
    ax=ax,
    color="white",
    edgecolor="black",
    alpha=0.7,
)

# basemap
cx.add_basemap(
    ax,
    crs=berlin.crs,
)

# plt.savefig('stackoverflow_question.png',
#             dpi=100,
#             bbox_inches='tight',
#            )

​

​