将 geopandas 几何体转换为 shapleys 多边形

问题描述

总的来说，我对 python 非常陌生，但我有不错的 Fortran 和 Matlab 编程背景。我想弄清楚什么似乎是一个简单的任务：

我正在使用 GeoPandas read_file 命令导入 Shapefile。它没有问题，并且包括几何列，该列包含由所有顶点表示的多边形，正如您所期望的那样。数据是单行。

我需要的是一个 Shapely 的多边形来传递到 py3dep，这样我就可以收集我需要的空间栅格，但我不知道如何将几何条目转换为多边形。

解决方案必须很简单，但我已经搜索了几个小时，尝试了一千种不同的方法，但都无济于事。如果有人可以帮助进行简单的演示，我将永远感激不尽。

这是我正在尝试做的一个片段，当我将元组传递给它时，我已经验证代码有效，如果我尝试将“shape”传递给“geom”则失败，因为 shape 不是 shapelys 多边形. （我知道这是一次蹩脚的尝试，只是显示为“我正在尝试做的事情”）

根据@PaulH 提供的解决方案更新以下代码

import py3dep
import geopandas as gpd


shape = gpd.read_file("inputshape.shp")
print(shape.type)

geom = shape.loc[0,'geometry']
dem = py3dep.get_map("DEM",geom,resolution=10,geo_crs="epsg:4326",crs="epsg:3857",output_dir="raster")

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

geopandas gis shapely