Итак, у меня есть ситуация, когда у меня есть тонна линейных разрывов разбитого маршрута, где мне нужно объединить их, используя линейку Shapely или Union или Postgis st_union. < /p>
< P> Сейчас моя идея состоит в том, чтобы использовать Shipely, чтобы импортировать линейные разрез в качестве типов геометрии. Соедините их или объедините их с помощью Shipely, а затем экспортируйте обратно в таблицу результатов в базе данных. Полем Как ... < /p>
01020000020e61000....
< /code>
Как я могу перевести это из базы данных в тип геометрии Python, используя Shapely, сделать некоторые манипуляции, а затем экспортировать его в базу данных? < /p>
В настоящее время это мой код, он просто импортирует эту строку geom object string из базы данных прямо сейчас и бросает ошибки, потому что это не геометрия. < /p>
def create_shortest_route_geom(shortest_routes):
conn = connect_to_database()
cur = conn.cursor()
shortest_route_geoms = []
for route in shortest_routes:
source = str(int(route[1]))
target = str(int(route[2]))
query = 'SELECT the_geom FROM public.ways WHERE target_osm = ' + target + ' AND source_osm = ' + source + ' OR target_osm = ' + source + ' AND source_osm = ' + target + ';'
cur.execute(query)
total_geom = cur.fetchone()
for index, node in enumerate(route):
try:
source = str(int(node))
target = str(int(route[index + 1]))
query = 'SELECT the_geom FROM public.ways WHERE target_osm = ' + target + ' AND source_osm = ' + source + ' OR target_osm = ' + source + ' AND source_osm = ' + target + ';'
cur.execute(query)
geom = cur.fetchone()
query = "SELECT ST_Union("+str(geom[0])+","+str(total_geom[0])+")"
cur.execute(query)
total_geom = cur.fetchone()
except IndexError:
print "Last element"
shortest_route_geoms.insert(total_geom)
return shortest_route_geoms
< /code>
Редактировать: я мог бы найти свой ответ здесь, подробнее рассмотреть его и обновить мой вопрос с ответом, если я выясню это. < /p>
Подробнее здесь: https://stackoverflow.com/questions/468 ... om-shapely
Импорт типа геометрии Postgis в Python в виде геометрии от Shapely? ⇐ Python
-
- Похожие темы
- Ответы
- Просмотры
- Последнее сообщение
-
-
Отношение «schema.table» не существует. Тестовый контейнер Hibernate Spatial Postgis.
Anonymous » » в форуме JAVA - 0 Ответы
- 30 Просмотры
-
Последнее сообщение Anonymous
-