使用GDAL在Python中按功能分割shapefile？
是否有可能在python中按功能分割一个shapefile？（最好是将临时生成的矢量对象临时保存到内存而不是磁盘的解决方案）。 原因：我想将gdal rasterizeLayer函数与shapefile的几个不同子集一起使用。该函数需要一个osgeo.ogr.Layer对象。 mkay，我尝试了一下，它可能如下工作。您可以按以下方式获取每个要素的gdal图层对象的几何形状。 # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() 现在，我只需要知道如何基于此几何形状创建osgeo.ogr.layer对象。 为了澄清。我需要使用纯Ogr / gdal代码编写的函数！这似乎也引起了其他人的兴趣，我仍然希望没有任何辅助模块的解决方案（尽管来自此处的任何解决方案都将在免费的可用qgis插件中使用）。