Получить значения только в пределах геопанд - PullRequest
Новая
       73

Получить значения только в пределах геопанд

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

У меня большой набор данных с геоданными. Я хочу обнародовать данные, относящиеся только к части города. Поэтому мне нужно создать форму и проверить, есть ли моя точка.

Я пытался работать с этим ответом :

ТОЧКИ

1) 

turin.head() = 

       latitude     longitude   
    0   44.9125     7.7432  
    21  45.0764     7.5249  
    22  45.0764     7.5249  
    23  45.0755     7.5248 
    24  45.0718     7.5236

2) 

geometry = [Point(xy) for xy in zip(turin.longitude,turin.latitude)]
turin_point = gpd.GeoDataFrame(turin,crs=crs,geometry=geometry)
turin_point.head()


    latitude    longitude   geometry
0   44.9125     7.7432  POINT (7.74320 44.91250)
21  45.0764     7.5249  POINT (7.52490 45.07640)
22  45.0764     7.5249  POINT (7.52490 45.07640)
23  45.0755     7.5248  POINT (7.52480 45.07550)
24  45.0718     7.5236  POINT (7.52360 45.07180)

ГРАНИЦЫ

1) 

border.head()

    longitude   latitude    
0   7.577835    45.041828   
1   7.579849    45.039877   
2   7.580106    45.039628   
3   7.580852    45.038576   
4   7.580866    45.038556

2) 

geometry2 = [Point(xy) for xy in zip(border.longitude,border.latitude)]
border_point = gpd.GeoDataFrame(border,crs=crs,geometry=geometry2)
border_point.head() = 

        longitude   latitude    geometry
    0   7.577835    45.041828   POINT (7.57783 45.04183)
    1   7.579849    45.039877   POINT (7.57985 45.03988)
    2   7.580106    45.039628   POINT (7.58011 45.03963)
    3   7.580852    45.038576   POINT (7.58085 45.03858)
    4   7.580866    45.038556   POINT (7.58087 45.03856)

Тогда в соответствии с ответить : 

turin_final= border_point.geometry.unary_union
within_turin = turin_point[turin_point.geometry.within(turin_final)]

IndexError: слишком много индексов для массива

1 Ответ

1 голос
/ 09 октября 2019

Если вы хотите найти точки в пределах вашей границы, сама граница должна быть полигоном, а не другим набором точек. Если координаты граничных точек расположены в правильном порядке, вы можете попробовать заменить последние две строки на эти: 

from shapely.geometry import Polygon

turin_final = Polygon([[p.x, p.y] for p in border_point.geometry])
within_turin = turin_point[turin_point.geometry.within(turin_final)]

Часто это не так, но, как я понял из ваших комментариев, это решило вашу проблему.

...