当循环多边形时，Python Shapely无法在多边形中找到点

Question

我有一个带有lat和lon的商店位置csv文件。我还有一个geojson文件，其中包含美国人口普查区的多边形特征。我想使用Python查看每个位置存在哪些多边形。

我使用Shapely Python库的包含（）来循环通过存储位置的csv文件，获取lat和lon坐标，然后检查这些坐标是否存在于Geojson文件中的一个多边形中。

如果我先遍历每个位置/坐标，然后遍历每个多边形，使用contains（）或infrance（）检查多边形是否包含该点，那么我的代码将成功运行。

但是这太慢了，所以我想反转这个过程，首先遍历polygons/geojson文件，然后根据数千个坐标检查数十万个多边形，而不是相反。

但是，当我简单地切换循环以对照坐标检查多边形时，contains（）不会找到任何匹配项，即使当我首先循环坐标并对照多边形检查它们时，contains（）会找到正确的匹配项。代码是一样的，只是颠倒了。

代码（contains显示false）：

with open(mapFile) as f:
    data = json.load(f)

    #Loop through features in Geojson file
    for feature in data['features']:

        polygon = shape(feature["geometry"])

        for row in locationsFile:
            #get point coordinates from file
            pointLat = float(row[13])
            pointLon = float(row[14])

            point = Point(pointLon, pointLat)

            print(polygon.contains(point))   
            if polygon.contains(point):

                #Grab data
                newx = feature["properties"]["x"]
                newy = feature["properties"]["y"]
                newz = feature["properties"]["z"]
                #append data
                row.append(newx)
                row.append(newy)
                row.append(newz)
                #update file
                newFile.writerow(row)
                break

此代码生成准确的结果：

with open(mapFile) as f:
    data = json.load(f)

    #Loop through coordinates in CSV file
    for row in locationsFile:
        #get point coordinates from file
        pointLat = float(row[13])
        pointLon = float(row[14])

        point = Point(pointLon, pointLat)

        #Loop through features in Geojson file
        for feature in data['features']:

            polygon = shape(feature["geometry"])

            if polygon.contains(point):

                #Grab data
                newx = feature["properties"]["x"]
                newy = feature["properties"]["y"]
                newz = feature["properties"]["z"]
                #append data
                row.append(newx)
                row.append(newy)
                row.append(newz)
                #update file
                newFile.writerow(row)
                break

Answer 1

break在这两种情况下的行为都不同。它只会中断最里面的循环，在“正确”的情况下，是通过数据['features']的循环，但在不正确的情况下，是通过位置文件的循环