возврат списка точек в геопандах полигона - PullRequest
Новая
       68

возврат списка точек в геопандах полигона

1 голос
/ 09 февраля 2020

У меня есть шейп-файл Лондона с более чем 4000 уникальными полигонами, представляющими области изоа. Я также создал геоданные с точками, представляющими датчики в городе. Мне нужно выяснить, в какой области ISO (ie какой полигон) каждый датчик принадлежит / находится внутри. Так что в идеале у меня был бы список каждого sesnsor_id_location и соответствующего номера lsoa (LSOA11D C).

Мои очки выглядят так: 

    pod_id_location Borough Latitude    Longitude   geometry
0   1245    Barnet  51.604486   -0.206551   POINT (-0.20655 51.60449)
1   2245    Camden  51.521880   -0.120434   POINT (-0.12043 51.52188)
2   3245    Camden  51.555485   -0.152338   POINT (-0.15234 51.55548)
3   5245    Wandsworth  51.440399   -0.186775   POINT (-0.18677 51.44040)
4   6245    Hounslow    51.468625   -0.359770   POINT (-0.35977 51.46863)

и мой файл лондонской формы выглядит так: 

    LSOA11CD    LSOA11NM    geometry
0   E01000001   City of London 001A POLYGON ((-0.09729 51.52158, -0.09652 51.52027...
1   E01000002   City of London 001B POLYGON ((-0.08813 51.51941, -0.08929 51.51752...
2   E01000003   City of London 001C POLYGON ((-0.09679 51.52325, -0.09647 51.52282...
3   E01000005   City of London 001E POLYGON ((-0.07323 51.51000, -0.07553 51.50974...
4   E01000006   Barking and Dagenham 016A   POLYGON ((0.09115 51.53909, 0.09326 51.53787, ...

Визуально мои данные выглядят так:

enter image description here

Спасибо за помощь.

1 Ответ

0 голосов
/ 09 февраля 2020

Я нашел очень простое решение с использованием пространственного соединения следующим образом: 

merge = gpd.sjoin(map_df, sensors, how="right", op='contains')


    index_left  LSOA11CD    LSOA11NM    pod_id_location Borough Latitude    Longitude   geometry
index_right                             
0   199 E01000204   Barnet 025C 1245    Barnet  51.604486   -0.206551   POINT (-0.20655 51.60449)
1   897 E01000915   Camden 027A 2245    Camden  51.521880   -0.120434   POINT (-0.12043 51.52188)
...