如何解决将 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 举报,一经查实,本站将立刻删除。