我正在寻找在geopandas中栅格化shpfile的最简单的方法-等同于在一行中执行操作的arcpy PolygonToRaster_conversion()。
我发现了一些相对复杂的方法,如https://snorfalorpagus.net/blog/2014/11/09/masking-rasterio-layers-with-vector-features/
有这么复杂吗?还是有像arcpy的PolygonToRaster_conversion()这样的一行选项?
我正在寻找最简单的起点来理解我的想法
我一直在探索rasterio来做这件事,但也许还有其他方法
我才刚刚开始使用Geopandas,如果有任何建议,我将不胜感激
发布于 2020-06-10 22:32:16
您是否正在尝试在一步中栅格化一组具有唯一值的多边形?如果是这样的话,你想对每个多边形使用唯一的值进行rasterize,但要注意最后一个被栅格化到给定像素的多边形将“声明”它(即,多个多边形可能会接触到一个像素,但您的特征列表中的最后一个将是那里的光栅化值)。
或者想要单独栅格化每个多边形(或同时栅格化所有多边形,就像它们是一个单独的多边形一样),以便可以从栅格中提取统计信息?Mask可能会在每个要素上循环执行此操作。
最接近单行操作的方法是使用rasterio's rio mask或rio rasterize操作。链接到的示例更复杂的原因是,您需要做一些额外的事情来提取原始栅格的子集。现在在rasterio中有一些额外的方法可以让这变得更容易(docs)。
在geopandas中,您的几何体位于GeoSeries中。我还没有直接测试过,但是您可能需要调用该系列的__geo_interface__来获取rasterio期望作为输入的类似GeoJSON的形状。
https://stackoverflow.com/questions/62242328
复制相似问题