geopandas栅格化shpefile

问题描述 投票:0回答:1

[我正在寻找在geopandas中光栅化shp文件的最简单方法-等效于arcpy PolygonToRaster_conversion()在一行中完成工作。

我发现了一些相对复杂的方法,例如https://snorfalorpagus.net/blog/2014/11/09/masking-rasterio-layers-with-vector-features/

这复杂吗?还是有一个单行选项,例如arcpy的PolygonToRaster_conversion()

我正在寻找最简单的起点来理解这一点

我一直在探索rasterio来做到这一点,但也许还有其他方法

我只是刚刚开始使用Geopandas,并且希望您能使用任何指针

geopandas rasterizing rasterio
1个回答
0
投票

您是否要一步一步将一组具有唯一值的多边形栅格化?如果是这样,您想rasterize为每个多边形使用该唯一值,但要注意栅格化到给定像素的最后一个多边形将“声明”它(即,多个多边形可能会碰到一个像素,但是列表中的最后一个多边形的特征将在那里栅格化的值)。

或者您是否要独立栅格化每个多边形(或同时对所有多边形进行栅格化,就好像它们是单个多边形一样,以便可以从栅格中提取统计信息?为此,Mask可能会在每个功能上循环播放。

最可能接近单行操作的地方是使用rasterio's rio maskrio rasterize操作。链接到示例的原因更多,是因为您需要做一些额外的事情来提取原始栅格的子集。 rasterio中现在有一些额外的方法可以使操作起来更加容易(docs)。

[从几何熊猫开始,您的几何图形位于GeoSeries中。我尚未直接对此进行测试,但是您可能需要调用该系列的__geo_interface__来获取rasterio期望作为输入的类似于GeoJSON的形状。

© www.soinside.com 2019 - 2024. All rights reserved.