将自定义Shapefile导入Basemap后如何设置填充颜色？

嘿，这个问题我之前也踩过坑！map.readshapefile() 方法本身确实不直接支持 fill_color 参数——它的核心设计是用来绘制矢量边界线条的，不像内置的drawcountries()、drawstates()这类方法已经封装了填充逻辑。不过咱们完全可以通过手动创建面要素集合的方式实现填充，下面给你两种可行的方案：

方案一：使用Matplotlib的PatchCollection（无需额外库） #

这是最直接的方法，步骤是读取Shapefile的坐标后，手动生成多边形对象，再打包成集合添加到地图上：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from matplotlib.patches import Polygon
from matplotlib.collections import PatchCollection

# 初始化你的Basemap对象（根据自己的需求设置投影、范围等）
m = Basemap(projection='merc', llcrnrlat=-30, urcrnrlat=30,
            llcrnrlon=-120, urcrnrlon=-60, resolution='l')

# 读取Shapefile，这里把drawbounds设为False，避免重复画边界
m.readshapefile('square', 'dir', drawbounds=False)

# 遍历每个形状的坐标，创建Polygon对象
patches = []
for coords in m.dir:
    # 把Shapefile的坐标转化为Matplotlib可用的Polygon
    poly = Polygon(coords, closed=True)
    patches.append(poly)

# 创建PatchCollection，设置填充色、边界色、线宽等样式
patch_collection = PatchCollection(patches, 
                                   facecolor='#87CEEB',  # 填充色
                                   edgecolor='black',    # 边界色
                                   linewidth=1.5,        # 边界线宽
                                   alpha=0.7)            # 透明度

# 将集合添加到当前地图的坐标轴上
ax = plt.gca()
ax.add_collection(patch_collection)

# 如果需要单独控制边界样式，也可以再调用一次readshapefile画边界
# m.readshapefile('square', 'dir', drawbounds=True, linewidth=1.5, color='k')

plt.title("自定义Shapefile填充示例")
plt.show()

方案二：结合Shapely和Descartes（适合复杂要素） #

如果你的Shapefile包含MultiPolygon这类复杂要素，用Shapely处理会更稳妥，不过需要额外安装shapely和descartes库：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from shapely.geometry import shape
from descartes import PolygonPatch
import shapefile

# 初始化地图
m = Basemap(projection='merc', llcrnrlat=-30, urcrnrlat=30,
            llcrnrlon=-120, urcrnrlon=-60, resolution='l')

# 用shapefile库读取要素（也可以用Basemap的readshapefile，不过这里直接读更灵活）
sf = shapefile.Reader('square')
for sr in sf.shapeRecords():
    shp_geom = shape(sr.shape)
    # 转换坐标到Basemap的投影下
    x, y = zip(*shp_geom.exterior.coords)
    proj_x, proj_y = m(x, y)
    proj_geom = shape({'type': 'Polygon', 'coordinates': [list(zip(proj_x, proj_y))]})
    # 创建填充Patch
    patch = PolygonPatch(proj_geom, facecolor='green', edgecolor='black', linewidth=1.5)
    plt.gca().add_patch(patch)

plt.show()

为什么内置形状能用fill_color？ #

内置的绘图方法（比如drawcountries()）是Basemap已经帮我们封装了从边界到面要素的转换逻辑，而readshapefile()只是把原始的坐标数据读取出来并绘制线条，所以需要我们自己手动处理填充的部分。

内容的提问来源于stack exchange，提问作者MattM