Le portail francophone de la géomatique

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 !!!!!!!!!!!!