前言 最近使用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 "打开失败" r