我有大量多边形(~100000)，并尝试找到一种聪明的方法来计算它们与规则网格单元的相交面积。目前，我正在使用shapely(基于它们的角坐标)创建多边形和网格单元。然后，我使用一个简单的for循环遍历每个多边形并将其与附近的网格单元格进行比较。只是一个小例子来说明多边形/网格单元。fromshapely.geometryimportbox,Polygon#Examplepolygonxy=[[130.21001,27.200001],[129.52,27.34],[129.45,27.1],[130.13,26.950001]]polygon_shape=Polygon(xy)#Ex