如何基于光栅遮罩删除边？

Question

我是Python地理空间工具箱的新手，所以我的问题听起来可能很基本，与OSMNX包无关。 我正在尝试识别并移除被洪水破坏的路段（即图中的边缘）。之后，我可以做我的网络分析

我使用OSMNX下载了我感兴趣的区域的道路图和geopanda数据框，我有一个光栅文件，显示了同一区域的水深。我知道什么样的水深使道路无法使用（比如说超过50厘米）。我想知道如何使用光栅文件将其与道路图相交，并识别必须删除的道路段？任何建议都会很有帮助，因为我目前一直在坚持这一点

Answer 1

我还不能留下评论，所以我加上它作为回答。除了@gboeing的答案外，您还应该能够根据所需水位对光栅进行多边形化。（使用光栅功能模块）

您可以使用多边形（来自光栅）和边（来自图形）来排除使用geopandas的sjoin或overlay功能超出标准的边

最后，通过OSMnx的graph\u from_gdfs再次构建该图

Answer 2

这是一个有点模糊的问题，所以我恐怕我能做的最好的回答就是有点模糊。我建议你研究一下：

1. rasterio用于python光栅处理
2. OSMnx的graph_to_gdfs函数，用于将图形转换为（向量）地理数据帧
3. 在python（tutorial）中使用光栅+矢量数据