Cry*_*tie 5 python intersection geopandas
我试图找到两组多边形的交集，其中一组是从人口普查数据产品（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），但我没有发现任何关于准确性的信息。我错过了什么吗？