Как импортировать файл .csv со столбцом координат для полигонов в геоданный? - PullRequest
Новая
       11

Как импортировать файл .csv со столбцом координат для полигонов в геоданный?

0 голосов
/ 13 декабря 2018

У меня есть CSV-файл, который включает столбец с координатами для ряда полигонов.

Например: 

  patch_ID  geometry                                                 flow   hab_area    Total   Portion of Flow     flow_area   Dispersal   mean_flow
1   156     POLYGON ((369000.0000019073 548000.0000019073,...   1.620000e-12    3874.00     3.740000e-11    4.328570    4.180000e-16    d4  3.081707
3   65      POLYGON ((362000.0000019073 560000.0000019073,...   1.010000e-12    404.25      3.740000e-11    2.709681    2.510000e-15    d4  1.447226
4   170     POLYGON ((383000.0000019073 520000.0000019073,...   8.770000e-13    5397.25     3.740000e-11    2.342323    1.620000e-16    d4  3.270493
6   179     POLYGON ((390000.0000019073 484000.0000019073,...   7.350000e-13    12341.50    3.740000e-11    1.963216    5.960000e-17    d4  3.796241
10  162     POLYGON ((385000.0000019073 541000.0000019073,...   3.530000e-13    1546.75     3.740000e-11    0.943354    2.280000e-16    d4  1.139087

Я импортировал его как фрейм данных pandas с помощью pd.read_csv (), как мне теперь преобразовать его в геоданные с такими геопандами, чтобы столбец Geometry былгеометрия многоугольников?Когда я пытаюсь: 

df1 = gdp.GeoDataFrame(df, geometry = 'geometry')

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

Затем, когда я пытаюсь проверить достоверность геометрии с помощью некоторого кода, который я нашел здесь 

for index, row in blbog_unpro.iterrows():
geom = row['geometry']
if len(geom.coords) <= 2:
      print("This row has an invalid polygon geometry")

Я получаю еще одну ошибку: у объекта 'str' нет атрибута 'координаты'

Я в тупике.

1 Ответ

0 голосов
/ 13 декабря 2018

Как уже упоминал DatHydroGuy, хитрость заключалась в использовании формата Well Know Text (WKT).Следующий код решил проблему: 

from shapely import wkt
df['geometry'] = df['geometry'].apply(wkt.loads)
df1 = gpd.GeoDataFrame(df, geometry = 'geometry')

Спасибо DatHydroGuy

...