Идея состоит в том, чтобы построить слои, имеющие базовую карту (r5_base), а затем наложить на них gpd 'regionbybranch'.
Как видно из изображений ниже, были созданы две фигуры, но я ожидаю только одну,Желаемая карта была построена правильно со вторым рисунком, но она не соответствовала желаемому «фигурному размеру»
# Isolate Region on map
r5_base = bars[bars.Reg_Name=='REGION V (BICOL REGION)']
# Isolate region from OP data
region=mpabars_gdf[mpabars_gdf.RegionName=='Bicol Region']
# Plot area.
regionbybranch = region.dissolve(by='BranchName', aggfunc='sum',as_index=False)
fig, ax = plt.subplots(1,1,sharex=True, sharey= True, figsize=(30,30))
base = r5_base.plot(color='gray')
regionbybranch.plot(ax=base, linewidth=0.5, edgecolor='black', legend= False, column='Penetration',cmap='viridis')
# plt.title('Market Penetration on Bicol Region');
# fig.savefig(“map_export.png”, dpi=300)