其实很简单
#先导入会用到的库
import geopandas as gpd
import matplotlib.pyplot as plt
import cartopy.io.shapereader as shpreader
# 整两个字典用于控制颜色图例之类的东西
color_dic = {}
zone_new_name_dic = {}
#建个画布
fig, ax = plt.subplots(1,1, figsize=(8, 6))
#然后就画吧
#这是好多个shape文件的路径
zone_folder = os.path.join("I:\IPCC6\Analysis2\\0.6 zone_ssp_trend\\a zone_com_extent")
#循环所有的shp文件
for zone in [x for x in os.listdir(zone_folder) if x.endswith(".shp")]:
# 把每个shp都画出来
zone_data = gpd.read_file(os.path.join(zone_folder, zone), encoding='utf-8').to_crs(4326)
print(color_dic[zone[1:2]])
ax = zone_data.to_crs(4326).plot(column=var,
linewidth=0,
legend=False,
ax=ax,
facecolor = "none", #这个就是设置pylogon是透明的
edgecolor = "k", #这个是设置边为黑的
lw = 1,
zorder = 0)
#最后show一下()
#也可以设置pylogon有颜色
#然后就画吧
#这是好多个shape文件的路径
zone_folder = os.path.join("I:\IPCC6\Analysis2\\0.6 zone_ssp_trend\\a zone_com_extent")
#循环所有的shp文件
for zone in [x for x in os.listdir(zone_folder) if x.endswith(".shp")]:
# 把每个shp都画出来
zone_data = gpd.read_file(os.path.join(zone_folder, zone), encoding='utf-8').to_crs(4326)
print(color_dic[zone[1:2]])
ax = zone_data.to_crs(4326).plot(column=var,
linewidth=0,
legend=False,
ax=ax,
color = color_dic[str(zone[1:2])], #这就有颜色了
edgecolor = "k", #这个是设置边为黑的
lw = 1,
zorder = 0)
#最后show一下()