Le portail francophone de la géomatique

Bonjour,

pour charger plusieurs dalles raster dans un projet qgis, j'utilise une action Python sur le shape contenant l'emprise des dalles. La table attributaire du shape contient un champ (dans l'exemple ci-après, "NOM_FICHIE") contenant lui même le nom du fichier raster (par exemple RGEALTI_FXX_0665_6640_MNT_LAMB93_IGN69.asc).
Au click sur la dalle avec l'icône "exécuter l'action", les dalles cliquées se chargent automatiquement. Si je clique à l'intersection de plusieurs dalles, plusieurs dalles se chargent.
Voici le code python, il doit être possible de l'adapter pour que les dalles à charger correspondent à la selection et non au click :

chemin = os.path.dirname( unicode( qgis.utils.iface.activeLayer().dataProvider().dataSourceUri() ) )
fic = chemin   + "/" + "[% "NOM_FICHIE" %]"
layer = qgis.utils.iface.activeLayer()
qgis.utils.iface.addRasterLayer(fic,"[% "NOM_FICHIE" %]")
qgis.utils.iface.setActiveLayer(layer)

Dans ce bout de code, le shape est stocké dans le même répertoire que les dalles raster.

En espérant que mes explications soient claires,

Emilie

Merci pour ces réponses.

En effet, il y a plusieurs sujets dans mon message.

Le premier étant d'ouvrir les fichier raster selon mes dalles déjà préparées (suppression des dalles inutiles). A la suite de quoi, je dois exporter mon raster pour le partager.
J'ai testé les deux codes python en vain. Surement une erreur de ma part (première fois que j'utilise cette fonctionnalité). Je pense que l'erreur est la même pour le code d'Emilie.

Lylambda voici mon code :

import os
def getVectorLayerByName(NomCouche):
    layermap=QgsProject.instance().mapLayers()
    for name, layer in layermap.items():
        if layer.name()==NomCouche:
            if layer.isValid():
               return layer
            else:
               return None
mypath=V:/DALLES_REG1/DALLE_SELEC_R1.shp.replace([%NOM%])
instRegistry = QgsProject.instance()
macouche = instRegistry.mapLayersByName([%NOM%])
courant=qgis.utils.iface.mapCanvas().currentLayer()
raster_ouvert=getVectorLayerByName([%NOM%])
if raster_ouvert is not None:
    QgsProject.instance().removeMapLayer(raster_ouvert.id())
    qgis.utils.iface.mapCanvas().refresh()
else:
    qgis.utils.iface.addRasterLayer(mypath,[%NOM%])
    qgis.utils.iface.setActiveLayer(courant)

Pour éviter les erreurs de frappes, j'ai cliqué sur l’icône "Dossier" pour créer mon chemin d'accès et cliqué sur "Insérer" pour déterminer l'attribut "NOM" de ma couche dalle. Mes fichiers raster jp2 se trouvent dans le même dossier "DALLES_REG1"

Merci !

Dernière modification par Bamakus (Mon 26 October 2020 16:50)

Sans trop creuser vos modifs du code, la variable "mypath" me parait incorrect :

originel :

mypath="[% "DOSSIER" %]".replace("\\","/")+ "[% "NOM2" %]"

> appel de la variable DOSSIER, qui est un chemin vers le dossier contenant la tuile, avec les \ (doublé pour échappement) remplacer par des / auquel est ajouté (+) la variable NOM2, qui est le nom du ficher jp2 de la tuile correspondant à l'emprise de la dalle

modifié :

mypath=V:/DALLES_REG1/DALLE_SELEC_R1.shp.replace([%NOM%])

> V:/DALLES_REG1/DALLE_SELEC_R1.shp est un chemin de fichier shp, et non un chemin vers un dossier. Il est écrit en dur et pas encadré de " (pas certain que ce soit nécessaire mais à tester). La fonction de remplacement veut y échanger la variable NOM sans préciser quoi.
Ce qui génère un chemin qui n'a pas de sens pour appeler une tuile.

Pour vous proposer une correction, je vais supposer :
- que V:/DALLES_REG1 est le dossier ...
- ...contenant votre couche de dallage DALLE_SELEC_R1.shp...
- ... et les dalles raster
- que la variable NOM est un attribut de la couche de dallage renseignant le nom des dalles sous la forme /XXXX.jp2 et rien d'autre

L'action à activer au clic doit donc avoir la ligne suivante :

mypath="V:/DALLES_REG1"+ "[% "NOM" %]"

> Cela permettra de générer pour chaque emprise, le lien d'appel suivant : V:/DALLES_REG1/XXXX.jp2

A adapter selon votre champ NOM et la véritable localisation du shp vis à vis des tuiles jp2.

Dernière modification par lylambda (Tue 27 October 2020 19:11)