基于python和GDAL实现对shape文件要素进行查询

前言 最近使用GDAL对shape文件进行查询,走了很多弯路,查看网上的代码,发现很多代码都是错误的,或者不详细,浪费了很多时间。这里简单记录一下使用GDAL查询和读取要素的代码。 python版本3.5 GDAL版本2.1.0 代码实现 读取shape代码 from osgeo import ogr # 注册所有的驱动 ogr.RegisterAll() ds = ogr.Open(path, 0) if ds is None: # 打开失败 print "打开失败" return # 获取该数据源中的图层个数,一般shp数据图层只有一个,如果是mdb、dxf等图层就