Le portail francophone de la géomatique

Merci pour votre réponse, mais je n'arrive justement pas a crée le shape avec la colonne location. Je possède néanmoins un shape livrée par IGN qui comporte le nom des rasters.
Donc je pense qu'il existe 2 solutions, soit réussir à faire ce fameux shape avec les chemins d'accès ( ce que je préférerais).
Soi utilisé des chemins relatifs et mettre le shape dans le même dossier que les rasters. Dans quel cas, seul nom des fichiers suffirait ?

Bonjour,

C'est justement cet outil : raster > divers > index des tuiles, que je n'arrive pas a faire fonctionner. En effet ce dernier me renvoie toujours l'erreur suivante ( dans les deux versions de Qgis) :
"Le processus n'a pu être lancé. Soit le programme requis est manquant, soit vous n'avez pas les droits nécessaires pour l'utiliser."

Non, malheureusement ce serrait trop simple

Bonjour,

Merci beaucoup, en effet cela fonctionne. J'avais déjà tenté la chose et le même message d'erreur (cité au dessus. Signe qu'il ne peux pas gérer plus de X fichiers ?) m’était apparût. En le fessant en 2 fois pas de message d'erreur... Ce serait plus pratique avec 1 seul fichier mais bon ... C'est mieux que rien.

Autre choses, a une grande échelle, les temps de chargement sont très long. J'ai vu sur votre lien qu'il parlaient d'un fichier unique qui serrait plus optimisé au codage même du Jpg2. Le fait d'avoir un fichier unique serait-il donc utile dans on cas ? De plus cela me permettrait de créer le VRT en 1 fois et donc d'avoir un seul fichier.

Et dernier point, le fichier VRT crée des "dalles noires" autours des Rasters, est il possible de ne pas les créer ou les supprimer ?

Bonjour,

En fin de message, vous trouverez un code action fonctionnel pour QGIS3. Il permet d'ouvrir ET de fermer au clic, une dalle raster à partir d'un quadrillage vecteur. Le code est à insérer dans "Propriétés de la couche" > "Actions" (type Python).

Typiquement, les couches "dalles.shp" des dossiers par départements (ex : "OHR_RVB_0M20_JP2-E080_LAMB93_D75-2014") de "3_SUPPLEMENTS_LIVRAISON" de la BDORTHO de l'IGN peuvent faire l'objet de ce code (ajouté manuellement ou via l'application d'un style) pour ouvrir/fermer les jp2 livrés dans les dossiers par départements homonymes de "1_DONNEES_LIVRAISON".

A noter que "DOSSIER" et "NOM2" sont des variables à adapter.
Dans mon cas, ce sont des champs virtuels générés par la calculatrice de champs, pour une utilisation en relatif :
- DOSSIER : chemin du dossier où sont situées les dalles raster. Exemple indiquant le chemin du dossier à partir de celui dans lequel est contenu "dalles.shp" :

replace((layer_property(@layer_id,'source')), array('3_SUPPLEMENTS_LIVRAISON','dalles.shp'), array('1_DONNEES_LIVRAISON',''))

- NOM2 : nom de la dalle raster à afficher. Exemple basé sur le champ "NOM" déjà inclus dans la table attributaire de "dalles.shp" : 

replace ( "NOM" , './', '' )

L'avantage vis à vis d'un VRT ou d'un projet incluant déjà les jp2 est l'ouverture bien plus rapide (seule les couches "dalle.shp" sont chargées). Les jp2 sont appelés au cas par cas selon les besoin, mais cela nécessite évidement une activation par l'utilisateur (le clic pour activer l'action).

En espérant que cela réponde à votre besoin !

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="[% "DOSSIER" %]".replace("\\","/")+ "[% "NOM2" %]"
instRegistry = QgsProject.instance()
macouche = instRegistry.mapLayersByName( "[% "NOM2" %]")
courant=qgis.utils.iface.mapCanvas().currentLayer()
raster_ouvert=getVectorLayerByName("[% "NOM2" %]")
if raster_ouvert is not None:
    QgsProject.instance().removeMapLayer(raster_ouvert.id())
    qgis.utils.iface.mapCanvas().refresh()
else:
    qgis.utils.iface.addRasterLayer(mypath,"[% "NOM2" %]")
    qgis.utils.iface.setActiveLayer(courant)

Dernière modification par lylambda (Fri 17 May 2019 00:22)

Un bémol pour cette méthode, les pixels noirs de l'orthophoto deviendront aussi transparents.

Bonjour,
j'essaye d'exploiter le code fourni par lylambda, mais sans résultats.
Ma dalle shp contient les champs : Img_2010, ex 18-2010-0630-6670-LA93.ecw et
Path_2010, ex \\cen-linda\gb_ref\BDORTHO\N_ORTHO_COUL\N_ORTHO_COUL_DALLES\N_ORTHO_COUL_018_2010_DALLES

J'ai personnalisé le code :
# le chemin du fichier à afficher:
fichier=r"[% c %]\[% Img_2010 %]"
# Nom à afficher dans la legende de Qgis:
titre="[% Img_2010 %]"
# nom du groupe où seront placées les dalles:
nomDuGroupe=u"Dalles"

lay=qgis.utils.iface.activeLayer()
titre=titre.decode('utf8')
fichier=fichier.decode('utf8')
if not os.path.exists(fichier):
  from PyQt4.QtGui import QMessageBox
  QMessageBox.information(qgis.utils.iface.mainWindow(), u"Attention la dalle ci-dessous est introuvable", fichier)
else:
  raster=qgis.utils.iface.addRasterLayer( fichier, titre )
  # Si le groupe voulu n'existe pas, il faut le créer puis déplacer la nouvelle dalle dans ce groupe :
  li = qgis.utils.iface.legendInterface()
  i = 0
  index = -1
  import qgis.core
  if QGis.QGIS_VERSION_INT>10800:
    for grp in li.groups():
      if grp ==nomDuGroupe: index = i
      i = i + 1
  else:
    for item in li.groupLayerRelationship():
      if item[0] == nomDuGroupe: index = i
      i = i + 1
  if index == -1: index = li.addGroup( nomDuGroupe )
  li.moveLayer(raster, index)
  # Redonner le statut 'Active' à la couche (pour continuer à utiliser cette Action) :
  qgis.utils.iface.setActiveLayer(lay)


MAIS QGIS 3.4.5 affiche une fenêtre Erreur Python (voir PJ)

D'avance merci

Cordialement

gerard