Я делаю мелкомасштабный рабочий процесс в блокноте ipython с геопандами и фигурными данными, извлекая кучу геопространственных данных, иногда из шейп-файлов, иногда из postgis (где выполняется более дорогая обработка).
Теперь я перетаскиваю таблицы postgis в python с помощью sqlalchemy
преобразования геометрии в WKT, получая что-то вроде:
sql = """
SELECT ST_AsText(ST_Transform(the_geom,4326)) as newgeom,*
FROM public.parcels2010_small limit 5;
"""
parcels = pd.read_sql(sql, engine)
parcels
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| id | newgeom | the_geom | parcel_id | osm_node_id |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| 0 | MULTIPOLYGON(((-122.991093691444 38.4878691106... | 01060000209C0E00000100000001030000000100000097... | 1805792 | 66237 |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| 1 | MULTIPOLYGON(((-122.444576448624 37.7346386006... | 01060000209C0E0000010000000103000000010000008A... | 10435 | 123826 |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| 2 | MULTIPOLYGON(((-122.796785208193 38.5427593334... | 01060000209C0E0000010000000103000000010000007D... | 1817842 | 313047 |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| 3 | MULTIPOLYGON(((-122.695538506163 38.3618570798... | 01060000209C0E0000010000000103000000010000009B... | 1934612 | 63776 |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
| 4 | MULTIPOLYGON(((-122.223424422869 37.8416019090... | 01060000209C0E00000100000001030000000100000072... | 861785 | 26369 |
+----+---------------------------------------------------+---------------------------------------------------+-----------+-------------+
Это выглядит очень похоже на геометрию при загрузке из шейп-файла, но это не приведено как геометрическая фигура. Я не мог найти канонический способ сделать это, используя только форму, или, возможно, с Descartes
.