Comment effectuer des tâches répétitives dans QGIS?

J'essaie de faire un traitement sur de nombreux fichiers de fonctionnalités, donc je voudrais l'automatiser.

En fait, j'ai un fichier de formes avec la distribution spatiale de certaines espèces et un avec le type de végétation.

Je voudrais sélectionner (par attribut) une espèce dans le fichier de formes des espèces, puis sélectionner (par localité) toutes les zones de végétation qui se croisent avec leur aire de répartition. Enfin, j'aimerais avoir un fichier de formes avec, comme nom, le nom de l'espèce, et avec les attributs et les formes des types de végétation fréquentés. Et je voudrais répéter cela pour toutes les espèces (plus de 100), et si possible, faites-le d'une manière facile (afin que cela puisse être fait par une autre personne).

J'ai déjà essayé cette tâche en utilisant le plugin Sextante, mais je ne peux pas avoir le nom de l'espèce comme nom de fichier de formes à la fin.

Quelqu'un peut-il suggérer une méthode pour cela?

Cette entrée de blog peut aider à comprendre comment le faire dans SEXTANTE:

http://qgissextante.blogspot.fr/2013/01/using-selection-algorithms.html

J'espère que ça aide

Cela nécessite un petit script. Pour le rendre reproductible je voudrais essayer de l' accomplir dans R . Cela devrait également être possible avec QGis et Sextante en utilisant l'exécution par lots (fonction clic droit sur la fonction) dans un modèle Sextante. Ici, vous pouvez d'abord utiliser l'outil d'intersection vectorielle et ensuite une sorte de jointure spatiale.

Dans R, je l'essayerais comme ça. Vous devrez peut-être modifier le code car je ne connais pas votre structure de données et vos variables.

library(raster);library(rgdal);library(sp)         # Load in required packages

vegetation <- readOGR("H:/Examplefolder",          # To load a vegetation polygon shape here 
                      "vegi")                      # called vegi.shp    

setwd(harddriveD)                                  # Now, switch to the directory containing your species shapes
                                                   # and use the following code 
species_files <- grep( ".shp",                     # to extract all shape names
                       dir(),
                       ignore.case=T,
                       value=T)

                                                   # Now read in your speciesfiles in a loop 
for(name in species_files){                        # and do a  vegetation data
                                                   # overlay with your basename
    spec_name <- strsplit(name,split=".shp")[[1]]  # to get only the load in
                                                   # your species name shape. 

    spec_shp <- readOGR(".",spec_name)             # I assume that you have point data. Otherwise change the code.
    ov <- over(spec_shp,vegetation)                # Make a overlay with your vegetation data, 
                                                   # returns a dataframe of the overlaying vegetation shape attributes, 
                                                   # which are within your species shape. 
                                                   # This dataframe has the same number of rows as your input species shape. 
   cd <- coordinates(spec_shp);                    # Therefore you could just append the needed information to your species shape.
   proj <- proj4string(spec_shp)                   # save coordinates and proj.

                                                   # Append your vegetation data to your species shape
   spec_shp$Vegetation <- ov$ShrubSpecies          # Or whatever you want to get. 

   spp <- SpatialPointsDataFrame(                  # In the end export your shape again. 
                    coords=cd,                     # I would advise you to use a different folder. 
                    data=as.data.frame(spec_shp),  # In case you have polygons instead of Point data
                    proj4string=CRS(proj) )        # use the SpatialPolygonDataFrame function. -> help ?SpatialPolygonDataFrame
  writeOGR(spp,                                    #Should result in a new shape 
           "foldername",                           # which has your species name.
           spec_name,
           driver="ESRI Shapefile")                      

}

J'ai émis de nombreuses hypothèses sur votre objectif et la structure de votre ensemble de données. Vous devez probablement corriger le code selon vos besoins avant de l'essayer.