Быстрый расчет пересечений - PullRequest
Новая
       45

Быстрый расчет пересечений

0 голосов
/ 23 октября 2019

Я хотел бы рассчитать площадь пересечения между файлом формы и сеткой рыболовной сети.

Поэтому я использую этот подход (в среде Python 3.7): 

import geopandas as gpd

fnet = gpd.read_file("foo/bar/fnet.shp")
shapes = gpd.read_file("foo/bar/shapes.shp")

geom_fnet  = [x for x in fnet.geometry]
geom_shapes   = [x for x in shapes.geometry]

intersections = []
for i, geolist1 in enumerate(geom_fnet):
    for j, geolist2 in enumerate(geom_shapes):
        if geolist1.intersects(geolist2):
            intersections.append((i, (geolist1.intersection(geolist2).area/geolist1.area)))

# next steps of analysis

У рыболовной сети большое количество плиток, и я хотел бы использовать подход быстрее, чем цикл for,Есть предложения?

...