Message posté par : Jérémie ---------------------------------------- Bonjour, Il y a plusieurs sujets dans votre message : - gagner du temps à l'ouverture de rasters - ouvrir les raster correspondant à des emprises vecteur sélectionnées à la main - fusionner les raster ou vrt pour stocker les raster dans une base postgres Pour ce qui est de gagner du temps, j'ai testé récemment la méthode de Klas Karlsson expliquée dans la vidéo suivante pour les tuiles d'un département (132 tuiles SCAN25 pour ~5,8Go de base) : https://www.youtube.com/watch?v=cRgkJzbx8gc L'idée est de stocker et tuiler + construire des pyramides pour les rasters, dans un gpkg. L'affichage est ensuite très fluide et comme on compresse au passage, on réduit la taille nécessaire, au dépit d'une perte en qualité qu'il faut essayer de limité par rapport à son besoin. Résumé des actions : 1- Ouvrir tous les fichiers raster dans QGIS 2- Créer un raster virtuel .vrt avec le menu de QGIS Raster > Divers > Contruire un raster virtuel (qui se base sur la commande gdalbuildvrt) Ne PAS cocher "Place each input file into a separate band" Préciser le SCR correspondant aux tuiles dans "Override projection", dans mon cas 2154 Lambert 93 3- Exporter le vrt en gpkg avec compression Clic droit sur la couche du vrt > Exporter > Enregistrer sous Choisir Geopackage et préciser l'emplacement et le nom du fichier souhaité Cocher Pyramides Pour Format de pyramide, choisir Interne (si possible) Cocher Niveaux personnalisés et indiquer dans le champ texte : 2 4 8 16 32 64 128 256 Dans Profil, choisir compression JPEG Modifier dans le champ texte en-dessous la valeur de la compression en 50 au lieu de 75 par défaut "JPEG_QUALITY_OVERVIEW=50" La création va prendre assez longtemps, surtout si comme pour le SCAN25 il y a beaucoup de rasters à la base. Concernant l'ouverture de raster par rapport aux emprises sélectionnées, je ne sais pas comment faire, peut-être via une action ou un bout de code python ? Concernant la fusion des raster, vrt en vue du stockage dans une base postgres, j'avoue ne pas trop comprendre votre question ? ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=336503#p336503 Pour y répondre : donnees(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net

Message posté par : EmilieDB ---------------------------------------- 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 : ----------------- Code : 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 ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=336508#p336508 Pour y répondre : donnees(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net