QGIS-2.6-PyQGISDeveloperCookbook-fr

PyQGIS developer cookbook, Version 2.6
pour une position donnée, vous devez déterminer les 10 points les plus proches dans une couche de points, dans
l’objectif d’utiliser ces points pour calculer une valeur interpolée. Sans index spatial, la seule méthode pour QGIS
de trouver ces 10 points est de calculer la distance entre tous les points de la couche et l’endroit indiqué et
de comparer ces distances entre-elles. Cela peut prendre beaucoup de temps spécialement si vous devez répeter
l’opération sur plusieurs emplacements. Si index spatial existe pour la couche, l’opération est bien plus efficace.
Vous pouvez vous représenter une couche sans index spatial comme un annuaire dans lequel les numéros de
téléphone ne sont pas ordonnés ou indexés. Le seul moyen de trouver le numéro de téléphone d’une personne est
de lire l’annuaire en commençant du début jusqu’à ce que vous le trouviez.
Les index spatiaux ne sont pas créés par défaut pour une couche vecteur mais vous pouvez le faire facilement de
cette manière :
1. créez l’index spatial — le code qui suit créé un index vide
index = QgsSpatialIndex()
2. ajouter les entités à l’index – l’index utilise des objets QgsFeature et les ajoute dans sa structure de
données interne. Vous pouvez créer les objets manuellement ou utiliser ceux qui sont issus de la méthode
nextFeature() du fournisseur de données :
index.insertFeature(feat)
3. Une fois que l’index est rempli avec des valeurs, vous pouvez lancer vos requêtes :
# returns array of feature IDs of five nearest features
nearest = index.nearestNeighbor(QgsPoint(25.4, 12.7), 5)
# returns array of IDs of features which intersect the rectangle
intersect = index.intersects(QgsRectangle(22.5, 15.3, 23.1, 17.2))
4.5 Ecrire dans des couches vecteur
Vous pouvez générer des fichiers de couche vecteur en utilisant la classe QgsVectorFileWriter. Elle gère
tous les formats vecteurs gérés par QGIS (fichier Shape, GeoJSON, KML, etc.).
Il y a deux façons d’exporter une couche vectorielle :
– A partir d’une instance de la classe QgsVectorLayer
error = QgsVectorFileWriter.writeAsVectorFormat(layer, "my_shapes.shp", "CP1250", None, "ESRI Sh
if error == QgsVectorFileWriter.NoError:
print "success!"
error = QgsVectorFileWriter.writeAsVectorFormat(layer, "my_json.json", "utf-8", None, "GeoJSON")
if error == QgsVectorFileWriter.NoError:
print "success again!"
Le troisième paramètre indique l’encodage du texte en sortie. Seuls certains pilotes ont besoin de ce
paramètre pour fonctionner correctement, les fichiers Shape sont dans ce cas. Néanmoins, vous ne devriez
pas rencontrer de problèmes tant que vous n’utilisez pas un jeu de caractères international. Le quatrième
paramètre que nous avons laissé à None peut indiquer un SCR de destination, si une instance valide de
QgsCoordinateReferenceSystem est utilisée, la couche est transformée dans ce SCR.
Consultez les formats gérés par OGR pour trouver les noms de pilote valides. Vous ne devez indiquer cette valeur
dans la colonne “Code”. Vous pouvez indiquer optionnellement d’exporter uniquement les entités sélectionnées
ou utiliser des options de création spécifiques à chaque pilote ou encore indiquer au pilote de ne pas créer
d’attributs. Consultez la documentation pour la syntaxe complète.
– Directement depuis les entités
4.5. Ecrire dans des couches vecteur 17