我试图找到两组多边形的交点，其中一组是从人口普查数据产品（poly1）中选择的一些人口普查区域，另一组由我自己定义的多边形（poly2）。我有兴趣查看哪些人口普查区域落入我的每个多边形（交叉点）中，并从那里计算相交的人口普查区域的面积等等。
我使用了以下代码：
import geopandas as gp
gp.overlay(censustracts, mypolygons, how = 'intersection')
这适用于涉及的人口普查区域不多的情况。我精心挑选了一些案例（5000~40,000 个人口普查区）并在 ArcMap 中进行了交集计算并交叉检查了结果。似乎 geopandas.overlay 给出的结果不太准确：有时给定的“相交”人口普查区不应该算作相交；有时缺少应该相交的人口普查区域。
此外，如果我只从我的自定义多边形 (poly2) 中选择一个多边形并选择它应该与之相交的正确人口普查区域，geopandas.overlay 可以正常工作。
我不知道为什么会这样。似乎只有关于 geopandas 覆盖速度的讨论（https://github.com/geopandas/geopandas/issues/400），但我还没有发现任何关于准确性的信息。我错过了什么吗？