Salvando um índice espacial na API QGIS?

Estou tentando criar um índice espacial para alguns shapefiles usando a API Qgis. Eu segui os passos descritos no blog de Nathan Woodrow ( https://nathanw.net/2013/01/04/using-a-qgis-spatial-index-to-speed-up-your-code/ ) e funciona muito bem :

layer = QgsVectorLayer(path, name, 'ogr')
idx = QgsSpatialIndex()
all_features = layer.getFeatures()
map(idx.insertFeature, all_features)

Meu problema é que não tenho nenhum arquivo resultante (.qix? .Sbn? .Sbx?)

Como posso salvar esse índice espacial para que futuros usuários dos shapefiles possam se beneficiar?

Se você deseja criar um índice espacial para todos os recursos da sua camada, tente usar o seguinte:

layer.dataProvider().createSpatialIndex()

Isso criará o .qixarquivo de índice espacial.

Editar:

Conforme mencionado por Matthias Kuhn em seus comentários abaixo, a classe QgsSpatialIndex é usada apenas temporariamente, como visto no blog ao qual você vinculou; enquanto a função createSpatialIndex () cria um índice espacial na fonte de dados e resulta em um qixarquivo para shapefiles.