Geopandas叠加交叉口

Question

我试图找到两组多边形的交集，其中一组是从人口普查数据产品（poly1）中选择的一些普查区域，另一组是我自己定义的多边形集（poly2）。我感兴趣的是，哪些普查区域属于我的每个多边形（交叉点），然后从那里计算出相交的普查区域的面积，等等。在

我使用了以下代码：

import geopandas as gp
gp.overlay(censustracts, mypolygons, how = 'intersection')

这对于没有太多涉及普查区的情况很有效。在ArcMap中手工选取了一些案例（5000~40000个普查小区），并在ArcMap中进行了交叉口计算，并对结果进行了交叉检验。看来geopandas.覆盖给出不太准确的结果：有时给定的“交叉”普查区域不应算作交叉点；有时，应相交的普查区域会丢失。在

此外，如果我只从我的自定义多边形（poly2）中选择一个多边形，并选择它应该与之相交的正确普查区域，geopandas.覆盖工作正常。在

我不知道为什么会这样。似乎只有关于geopandas叠加速度的讨论（https://github.com/geopandas/geopandas/issues/400），但我没有发现任何关于精度的东西。我错过什么了吗？在