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