[Qgis_fr] Extraire et utiliser des flux wms ppnc du serveur RW puis les utiliser

Message posté par : Pierre Coll (pierrecolla(a)gmail.com) ---------------------------------------- Bonjour, j'utilise qfield pour mon boulot de conservation de la nature. en son temps, j'ai extrait les images aériennes des serveurs de la région wallone puis je les ai traitées pour les alleger et les utiliser sur le terrain hors connection. j'ai fais cela en 2021 cela a bien marché mais c etait tres chronophage. Aujourd hui, les nouvelles photos aériennes sont sorties et j'aimerais bien recharger ces nouvelles données. je vous mets ci dessous ma procédure. La question est: est-il possible d'automatiser cela avec un script? D'avance merci pour vos réponses. Procédure: Extraire et utiliser des flux wms ppnc du serveur RW puis les utiliser sur qfield : 1 :extraire les données wms en TIFF Impossible d’automatiser et de faire un traitement par lot. Charger le wms sur le projet. Clic droit et exporter ==> enregistrer sous=> geotiff==>decocher « cree un VRT »==>enregistrer sous ==>zoomer sur un site N2K et cliquer sur etendue du canevas de la carte. Resolution H et V 0.5 (ou 0.25 mais plus lourd) Option de creation : profil : compression elevée (plutôt compression jpeg avec l’expérience) ==> Executer le script. Pour reduire fortement la taille du fichier il y a moyen de decouper le raster obtenu selon une couche masque : -Boite a outils de traitement==> GDAL==>decouper raster selon une couche masque. -choisire la couche source (le carré ortho du site N2K), la couche masque ( le périmetre N2k +20m de buffer), cocher créer une bande transparente, cocjer conserver la resolution du raster d’entree, choisir profil : compression elevée (plutôt compression jpeg avec l’expérience), choisir un dossier d’enregistrement. 2: Compresser les Tiff en Tiff conpressé jpeg Gdal==>conversio raster==> convertir Clique sur la roue dentee des paratmetres pour afficher le champs « option de création supplementaire (commande « COMPRESS=JPEG|JPEG_QUALITY=50 » dans option de création supplementaire en exécution de traitement de lot) 3: crée un virtual raster 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 31370 Cocher add aplpha mask Dans le champs no data ajouter : « 0 » 4)Modifier les bandes du VRT : Gdal==> conversion raster ==> reordonner les bandes==> ne garder que les bandes 1 (R) 2(G) 3(B) 4 (A) 5 : Exporter le VRT en gpkg Mettre tt les tiff et le VRT danc un dossier C:/traitement Ouvrire osgeo4 shell Taper : « cd traitement » (cd= change directory = changer de dossier Puis taper e code suivant : « gdal_translate --config OGR_SQLITE_SYNCHRONOUS OFF -co APPEND_SUBDATASET=YES -co TILE_FORMAT=WEBP -a_srs EPSG:31370 -of GPKG lenomdemonvrt.vrt lenomdemongpkensortie.gpkg gdaladdo --config OGR_SQLITE_SYNCHRONOUS OFF -r AVERAGE lenomdemongpkensortie.gpkg 2 4 8 16 32 64 128 256 » !!!!!!!!!!!!!!!!!!!!!!!!!Attention : il faut charger le gpkg par glisser du dossier vers deposer dans le projet !!!!!!!!!!!! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=373621#p373621 Pour y répondre : qgis_fr(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