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

时间:2021-04-03 00:32:58

标签: gis geopandas shapely

总的来说,我对 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")

0 个答案:

没有答案