Повышение производительности наложения геоданных (пересечение)

Question

У меня есть два набора шейп-файлов с полигонами. Один набор шейп-файлов - это только те округа США, которые меня интересуют, и это зависит от разных фирм и лет. Другой набор шейп-файлов - это бизнес-направление фирм, и, конечно же, это варьируется в зависимости от фирм и лет. Мне нужно каждый раз пересекать эти два слоя для каждой фирмы. До сих пор функция overlay (df1, df2, how = 'intersection') достигла моей цели. Но для каждого года фирмы требуется около 300 долларов. Учитывая, что у меня длинный список фирм и много лет, мне потребовались бы несколько дней, чтобы закончить. Есть ли способ повысить производительность?

Я замечаю, что если я сделаю то же самое в ArcGIS, 300 секунд сработают до нескольких секунд. Но я новый пользователь ArcGIS, еще не знакомый с питоном.

Answer 1

Надеюсь, вы уже поняли это, но решение заключается в использовании пространственного индекса R-дерева Geopanda. Вы можете добиться улучшения порядка, реализуя его надлежащим образом.

Goeff Boeing написал excellent tutorial.

http://geoffboeing.com/2016/10/r-tree-spatial-index-python/