使用descartes PolygonPatch与shapely绘图触发IndexError问题求助

核心原因 #

Descartes是已停止维护的老旧库，完全不兼容Shapely 2.0及以上版本。Shapely 2.0重构了几何对象的内部实现，原Descartes代码中直接将LinearRing（圆形的外环）转换为numpy数组的逻辑失效，导致索引越界错误。

替代解决方案 #

方案1：使用Shapely原生Matplotlib支持（推荐） #

Shapely 2.0+已内置Matplotlib集成，无需依赖Descartes：

import shapely.geometry as sg
import matplotlib.pyplot as plt

circle = sg.Point((0, 0)).buffer(1)

fig, ax = plt.subplots()
# 直接调用plot方法绘制几何对象
ax.plot(circle, facecolor='lightblue', edgecolor='darkblue')
ax.set_aspect('equal')  # 保证圆形不变形
plt.show()

方案2：手动构建Matplotlib Patch #

如果需要更精细控制Patch样式，可手动提取坐标构建PathPatch：

import shapely.geometry as sg
import matplotlib.pyplot as plt
from matplotlib.patches import PathPatch
from matplotlib.path import Path

circle = sg.Point((0, 0)).buffer(1)

# 提取外环的所有坐标点
exterior_coords = list(circle.exterior.coords)
# 创建Matplotlib Path对象
path = Path(exterior_coords)
# 创建Patch并设置样式
patch = PathPatch(path, facecolor='pink', alpha=0.6, edgecolor='red')

fig, ax = plt.subplots()
ax.add_patch(patch)
ax.set_aspect('equal')
ax.set_xlim(-2, 2)
ax.set_ylim(-2, 2)
plt.show()

方案3：降级Shapely到1.x版本（不推荐） #

如果必须使用Descartes，可降级Shapely到1.8.x版本：

pip install shapely==1.8.5

错误细节解释 #

原代码中descartes.PolygonPatch调用PolygonPath时，尝试执行asarray(t.exterior)[:, :2]，但在Shapely 2.0+中，t.exterior返回的是LinearRing对象，直接转numpy数组会得到0维数组，无法用二维索引[:, :2]访问，因此触发IndexError。

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