OGR Info - Vérifier l'index spatial?

Si nous voulons obtenir des informations sur une couche vectorielle (fichier de formes par exemple), nous pourrions utiliser ogrinfo:

ogrinfo -so myshapefile.shp mylayer

Cela renvoie beaucoup d'informations pertinentes sur le fichier de formes. Est-il possible de renvoyer s'il existe ou non un indice spatial? Si oui, pouvons-nous accéder à des informations sur l'indice spatial?

Remarque: Je me rends compte qu'il est trivial d'ouvrir le dossier contenant le fichier de formes et de rechercher un .sbn ou .sbx et donc la question est spécifique à ogrinfo.

En réponse au commentaire de @Richard Law, le nom de la table virtuelle rtree qui contient l'index spatial dans GeoPackage est toujours nommé selon le modèle rtree_[table_name]_[geometry_column_name]

Par conséquent, vous pouvez effectuer une requête SQL qui vérifie l'existence de la table rtree.

Par exemple, vous pouvez vérifier si la table "table1" a un index spatial avec ogrinfo

ogrinfo -sql "SELECT EXISTS(SELECT 1 FROM sqlite_master WHERE tbl_name like 'rtree_table1_%') as has_spatial_index" mygeopackage.gpkg


INFO: Open of `mygeopackage.gpkg'
      using driver `GPKG' successful.

Layer name: SELECT
Geometry: None
Feature Count: 1
Layer SRS WKT:
(unknown)
has_spatial_index: Integer (0.0)
OGRFeature(SELECT):0
  has_spatial_index (Integer) = 1

Dans ce cas, l'index spatial existe. Sinon, le aurait été "0".

Dans GeoPackage, une table ne peut avoir qu'un seul champ de géométrie et il suffit donc de faire une vérification en "rtree" et le nom de la table sans connaître le nom du champ de géométrie comme rtree_table1_%.

Avec GDAL 2.2.0 , l'une des modifications apportées au pilote GPKG est:

• ajouter la fonction SQL HasSpatialIndex (tblname, geomcolname)

Par exemple:

ogrinfo -sql "SELECT HasSpatialIndex('some_layer', 'geom')" some-file.gpkg

Montre:

INFO: Open of `some-file.gpkg'
      using driver `GPKG' successful.

Layer name: SELECT
Geometry: Unknown (any)
Feature Count: 1
Layer SRS WKT:
(unknown)
HasSpatialIndex: Integer (0.0)
OGRFeature(SELECT):0
  HasSpatialIndex (Integer) = 1

Ou mieux, sélectionnez de gpkg_geometry_columnspour afficher les informations de toutes les colonnes de géométrie du fichier:

ogrinfo -sql "SELECT table_name, column_name, HasSpatialIndex(table_name, column_name) FROM gpkg_geometry_columns" some-file.gpkg

Ogr n'utilise pas les fichiers d'index spatial ESRI .sbn, sbx. Il crée un fichier .qix qui peut contenir un index spatial et / ou un index d'attribut. Vous pouvez créer le .qix via ogrinfo.
Pilote de fichier de formes Ogr

J'ai également trouvé un ticket GDAL mentionnant que le fichier .sbn est maintenant lisible. Je ne sais pas ce que cela signifie.
billet # 4719