Python处理空间数据：GeoPandas读取SHP与地图绘制

本文深入解析了使用GeoPandas处理Shapefile空间数据时最常遭遇的四大痛点：路径含中文或空格引发的“DriverError”、坐标系未统一导致的地图严重歪斜、plot绘图卡顿与图例异常，以及叠加底图时因CRS不匹配而出现的白屏或错位问题；通过实操级建议——如用pathlib安全构造路径、强制检查并转换至EPSG:3857、明确指定column/cmap/scheme参数、简化几何精度及切换底图源——手把手帮你绕过环境陷阱与认知盲区，让空间数据读得进、转得准、画得快、叠得对。

GeoPandas读SHP文件报错“DriverError: unable to open”

常见于路径含中文、空格或特殊字符，或GDAL底层驱动没加载成功。GeoPandas依赖fiona读取SHP，而fiona对路径敏感程度远超普通Python库。

• 用os.path.abspath()转成绝对路径，再传给gpd.read_file()，避免相对路径跳转出错
• 路径里有中文？先用pathlib.Path构造路径对象，再转字符串——比直接拼接字符串更稳
• 检查fiona.list_drivers()是否包含"ESRI Shapefile"，不包含说明GDAL环境异常，不是GeoPandas本身的问题
• 别用Windows资源管理器复制的路径（带隐藏的\r\n或全角空格），手动重敲一遍最可靠

读进来的GeoDataFrame坐标系是乱的，plot出来地图歪斜

Shapefile自带.prj文件定义坐标系，但GeoPandas默认不强制校验或转换。很多公开数据用WGS84（EPSG:4326），但绘图时误当Web Mercator（EPSG:3857）处理，就会拉伸变形。

• 读取后立刻检查gdf.crs，不是None才说明有坐标系信息
• 若gdf.crs为空，用gdf.set_crs(epsg=4326, inplace=True)手动设——但得确认原始数据确实是WGS84经纬度
• 绘图前统一转到Web Mercator：gdf.to_crs(epsg=3857)，否则matplotlib或contextily底图会错位
• 别在plot()里加transform=...参数试图“动态投影”，GeoPandas的plot不支持运行时重投影

用gdf.plot()画图卡顿、颜色不对、图例不显示

默认行为太“智能”，反而容易踩坑：自动按几何类型选色、忽略NaN、不处理分类字段的离散映射。

• 颜色混乱？加column="your_column"明确指定字段，并用cmap="viridis"或cmap="tab10"控制色阶类型
• 想画分类图但颜色混成渐变？必须加scheme="quantiles"或scheme="natural_breaks"，且legend=True才生效
• 图例文字太小？用legend_kwds={"loc": "lower left", "fmt": "{:.0f}"}微调，别指望默认适配所有屏幕
• 大文件（>10MB）直接plot()会卡死——先gdf.simplify(tolerance=0.001)降点精度，尤其处理高精度测绘数据时

叠加底图时出现“CRS mismatch”或白底无图

底图（比如contextily）严格要求Web Mercator（EPSG:3857），而你的GeoDataFrame还在WGS84里，坐标根本对不上。

• 底图加载失败，先确认gdf.crs.to_epsg() == 3857，不是就gdf = gdf.to_crs(epsg=3857)
• ctx.add_basemap(ax)必须在gdf.plot(ax=ax)之后调用，顺序反了底图会盖住矢量图
• 国内用户常遇到contextily连不上OpenStreetMap——换源：ctx.add_basemap(ax, source=ctx.providers.Esri.WorldImagery)
• 别用plt.show()前反复调用add_basemap，它只渲染一次，重复调用不报错但也没用

坐标系和路径是两个最常被跳过的环节，一错就整条链路失效。不是代码写得不够多，而是这两步没盯住。

今天关于《Python处理空间数据：GeoPandas读取SHP与地图绘制》的内容介绍就到此结束，如果有什么疑问或者建议，可以在golang学习网公众号下多多回复交流；文中若有不正之处，也希望回复留言以告知！