GDAL/OGR 创建空间索引

GDAL/OGR 创建空间索引

命令行的方式

$ ogrinfo -sql "CREATE SPATIAL INDEX ON property_view" "property_view.shp"

Python的方式

from osgeo import ogr
drv = ogr.GetDriverByName('ESRI Shapefile')
ds = drv.Open(shp_file, 1)
lyr = ds.GetLayer()
lyr_name = lyr.GetName()
sql = 'CREATE SPATIAL INDEX ON ' + lyr_name 
ds.ExecuteSQL(sql)
ds = None