Le portail francophone de la géomatique

C'est possible en utilisant une action python. Il faut que la table SHP des orthos ait un champ "chemin" avec le chemin et le nom de chaque fichier.
C'est un peu mariol mais il existe pas mal de tuto pour utiliser les actions.
Il est alors possible de charger chaque ortho à la demande.

Bonjour,
Il y a bien ce vieil article du portailSIG pour QGIS 2.x

Euh... console Python? Je n'en trouve pas référence dans le lien envoyé. L'article passe par le menu Raster. Sans la phase 1, je ne pense pas que la 2 soit réellement possible/utile.
Vous êtes bien en version 2.x? Auquel cas avez-vous (aussi) essayé l'extension Image Boundary?

J'ai également tenté avec Raster/Divers/ et ça bloque.

En ligne de commande ça fonctionne: (Tu fermes Qgis)

1 - Tu listes tes images dans un fichier texte avec cette commande dos. (Tu dois être dans le répertoire où sont les images.)

dir *.ext c:\vrt> listeraster.txt

remplaces ext par l'extension de tes images

2 - Tu l'ouvres avec un éditeur de texte type notepad++ et tu gardes uniquement les chemins de fichiers.

3 - Depuis le répertoire où se situe gdalbuildvrt.exe C:\Program Files\QGIS 3.4\bin normalement, dans la fenêtre dos tu tapes cette ligne de commande.

gdalbuildvrt -resolution highest -r nearest -input_file_list c:\vrt\listeraster.txt c:\vrt\listeraster.vrt

Après tu ajoutes le fichier listeraster.vrt comme une couche dans Qgis.
Avant de l’insérer, tu zoomes de façon à ne pas charger toutes les images.

Dernière modification par michel56 (Mon 14 January 2019 16:59)

Pour la console dos, tu vas dans démarrer et tu tapes cmd

l'extension c'est jp2. (L'ign fournie ce format maintenant.)

dir *.jp2 c:\vrt> listeraster.txt

pour tu déplacer dans un répertoire tu tapes: cd et le nom du répertoire.

Dernière modification par michel56 (Mon 14 January 2019 17:23)

par contre il y a un truc qui peut peut etre faire bloquer, c'est que les orthos sont dans un dossier sur un serveur externe (et pas sur mon bureau/pc = pas dans c)

du coup le chemin d'accès via C me semble faux


edit : je n'ai pas fais de vrt (je le dis vu qu'il y a vrt dans la ligne de commande).

Je crains que ca ne soit trop complexe pour moi cette affaire à moins d'être patient ^^

Dernière modification par bill08 (Mon 14 January 2019 17:33)

je n'arrive pas à créer de vrt, chaque fois j'ai droit à un message d'erreur : 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.


je crois que c'est le fait qu'il y ai trop de fichiers (1500 dalles) . j'ai essayé avec 2 rasters ca a marché



je vais essayer en 2/3 fois ?

Dernière modification par bill08 (Mon 14 January 2019 17:41)

Pour faire suite à ce qu'à expliqué E.Gonthier, vu que j'ai déjà testé la méthode des différents tutos il y a pas trop longtemps voilà ce que tu peux faire :

1-tu fais une copie de ton shp table d'assemblage un niveau au dessus de ton répertoire de dalles raster si ce n'est pas déjà le cas.

2-ouvrir cette couche sous qgis--> là tu as déjà normalement un champ avec le nom des dalles.jp2 ou autre on appellera "dalle" ; soit tu crées un deuxième champs "repertoire" et tu le remplis avec le chemin du répertoire contenant, en absolu, ou relatif par rapport à ta copie de table d'assemblage avec de souvenir des / au lieux de .

3-puis dans les propriétés de ta couche (clic droit) dans rubrique Action tu définis type : Python ; Description (nom de l'action): Ouvre dalle par ex

Texte de l'action (copié collé si tu as bien fait avec deux champs repertoire+dalle) :

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

Voilà, paré pour ouvrir en un clic tes ortho via ton vecteur d'assemblage par une action python ; Remarque tu peux shorter l'étape 2 et mettre directement à la place de [%"repertoire"%] -->[%'/Monrepertoiredortho/'%] comme çà tu ne modifie même pas la table d'assemblage fournie avec tes orthos si les raster sont tous au même endroit.

Je te laisse imaginer pour faire une action pour refermer les dalles.

Vraiment pratique sous QGIS par contre sur des tablettes Android avec QField je ne pense pas, mais ce n'est peut être pas ton cas.
Cdlt

Tu regardes avec ce code, le fichier listeraster.txt va être créé dans le même répertoire que tes orthos.

dir *.jp2 > listeraster.txt

Dernière modification par michel56 (Tue 15 January 2019 10:23)

Qu'entendez vous par 2 - Tu l'ouvres avec un éditeur de texte type notepad++ et tu gardes uniquement les chemins de fichiers.

ouvrir avec un editeur de texte (ici j'ai wordpad). J'ouvre ce doc, j'ai donc en première ligne, l'intégralité du chemin d'accès à mes orthos, ensuite 3 colonnes Mode / LastWriteTime / Length Name puis les infos pour chaque orthos.

Que dois-je supprimer ? j'ai essayé différents trucs, ca a pas marché ensuite avec la commande.
Si je dois supprimer toute la colonne lastwrite time (par exemple), comment faire ? car je ne peux pas sélectionner la colonne pour la supprimer (ca sélectionne chaque ligne).
et en mode wordpad, ca met une ligne avec le nom de l'ortho et en dessous des dates et heures...donc supprimer 1 par 1 j'en ai pour 107ans

merci par avance

Dernière modification par bill08 (Tue 15 January 2019 10:48)

la j'ai ca :
il faut donc que je supprime les données horaires/dates, je peux pas automatiser leur suppression ? car faire ca 1500 fois



    Répertoire: \\adresse IP\documents\SIG \FONDS DE PLAN\Données 2015\Orthophotos
    2013\ortho\2013\25cm\Nouveau dossier


Mode                LastWriteTime         Length Name                                                                 
----                -------------         ------ ----                                                                 
-a----       22/02/2014     19:55       11517300 6922-64692.jp2                                                       
-a----       22/02/2014     19:55       10095921 6922-64707.jp2

ok j'ai saisi le truc, faut que je me replonge dans la manip mais j'ai eu fait.

par contre, sur ton exemple, tu as le chemin d'accès pour chaque orthos, moi dans mon truc je l'ai juste une première fois, est-ce que je dois avoir comme toi ou quoi ? (c'est à dire un chemin d'accès complet/ligne)

bon j'ai bien modifié le chemin d'accès pour l'adapter à mon pc ^^
je me met dans le fichier avec mes orthos, pas de message d'erreur mais pas de vrt non plus
ca me met juste >>

et ca me remet meme pas le chemin qui se met automatiquement à l'ouverture de la commande

voici ce que je tape : C:\Program Files\QGIS 2.18\bin\gdalbuildvrt.exe" fichier.vrt *.jp2

Dernière modification par bill08 (Tue 15 January 2019 11:24)

effectivement, le double tiret manquait, ca donne un message d'erreur :

Au caractère Ligne:1 : 51
+ "C:\Program Files\QGIS 2.18\bin\gdalbuildvrt.exe" fichier.vrt *.jp2
+                                                   ~~~~~~~~~~~
Jeton inattendu « fichier.vrt » dans l’expression ou l’instruction.
    + CategoryInfo          : ParserError: ( [], ParentContainsErrorRecordException
    + FullyQualifiedErrorId : UnexpectedToken