Convertir des images TIFF en MBtiles dans QGIS

 Convertir images TIFF en MBtiles dans QGIS

Ceci peut intéresser les utilisateurs d'OruxMaps (application mobile de géolocalisation).
Cette application permet d'importer de multiples fonds de cartes, soit en-ligne (WMS/WMTS par exemple), soit hors-ligne. Pour Oruxmaps, le format privilégié pour les cartes hors-ligne est le  format Oruxmaps SQL ou Oruxmaps SQL lite. Le logiciel MOBAC permet de convertir de nombreuses sources de cartes en-ligne au format Orux SQL lite et donc d'en disposer hors-ligne sur le terrain.

Malheureusement, OruxMaps est incapable de lire directement des tuiles au format TIFF. Mais il est capable d'ouvrir des tuiles au format MBtiles. Une conversion TIFF --> MBtiles est donc nécessaire.

Pour cela, on pourra utiliser QGIS (version 3.2.4).

1) Glisser le fichier image TIFF géolocalisée du disque dur vers la fenêtre "couches". Vérifier le code EPSG.

2) Sélectionner cette image

3) Dans le menu : Traitement > Boîte à outils > Outils Raster > Génère des tuiles XYZ (MBtiles)

4) Dans la boîte de dialogue qui s'ouvre, introduire les informations :
    
    - Extent : cliquer sur la flèche "Set to current map canvas extent"
    - Zoom minimal : 12
    - Zoom maximal : 19
    - DPI : 100 (test provisoire)
    - Format de tuile : PNG
    - Fichier en sortie (pour MBtiles) : Choisir "enregistrer vers un fichier" et fixer son chemin.

5) Exécuter.

    

    

Des courbes de niveaux équidistantes de 1 mètre dans QGIS

1.  - Accéder au site geoservices.ign.fr

- dans le menu : cliquer sur "catalogue"
- choisir RGE ALTI

La précision du MNT varie en fait selon les zones et les enjeux : elle est portée à 20 cm EMQ dans les zones inondables ou littorales pour répondre aux besoins liés à la mise en oeuvre de la directive européenne inondation. Le RGE ALTI® est mis à jour à partir des levés obtenus par LIDAR aéroporté ou par corrélation d'images aériennes.

a)  Dans la section "Documentation" (bouton, colonne de droite)

     voir Tableau d'assemblage et récupérer le tableau d'assemblage à 1 m.
     https://geoservices.ign.fr/sites/default/files/2021-07/IGNF_RGEALTIr_2-0.xml

b) Récupérer ensuite les tuiles par département. (par exemple le 76).
      sur  RGE ALTI® 1 m
    Il s'agit de liens ftp
      exemple :
       ftp://RGE_ALTI_ext:Thae5eerohsei8ve@ftp3.ign.fr/RGEALTI_2-0_1M_ASC_LAMB93-                      IGN69_D076_2020-10-20.7z
       
        Pour télécharger ces tuiles : utiliser pour cela Filezilla (client FTP)

       Copier le lien précédent dans la boîte "hôte" de Filezilla.
       Le téléchargement est long (plusieurs heures si la connexion internet est mauvaise).
       Le fichier est zippé et fait quelques Go.
       Dézipper ce fichier.
       Le dossier résultant est très lourd (>50 Go)

    2. - Ouvrir QGIS

a) Ajouter un des fonds de carte suivants :

a1) Ajouter la couche OpenStreetMap dans XYZ Tiles
https://tile.openstreetmap.org/{z}/{x}/{y}.png
Assez long à afficher
Décocher le rendu en zoomant Accéder au menu SCR du projet et mettre EPSG 2154 (= Lambert 93)
a2) Ajouter la couche Google Satellite dans XYZ Tiles
http://mt0.google.com/vt/lyrs=s&hl=en&x={x}&y={y}&z={z}
a3) Ajouter la couche Google Maps dans XYZ Tiles
https://mt1.google.com/vt/lyrs=r&x={x}&y={y}&z={z}
a4) Ajouter la couche Bing dans XYZ Tiles
http://ecn.t3.tiles.virtualearth.net/tiles/a{q}.jpeg?g=1
a5) Ajouter Google Satellite Hybride dans XYZ Tiles
https://mt1.google.com/vt/lyrs=y&x={x}&y={y}&z={z}
a6)Ajouter Google roads dans XYZ Tiles
https://mt1.google.com/vt/lyrs=h&x={x}&y={y}&z={z}
a7 ) Ajouter Plan IGN v2 ou orthophotos dans flux WMS/WMTS
https://wxs.ign.fr/decouverte/geoportail/wmts?SERVICE=WMTS&REQUEST=GetCapabilities


b) Ajouter le GRAPHE MOSAIC (tableau d'assemblage) aux couches précédentes
Donner de la transparence à la mosaïque (propriétés de la couche)

c) Ajouter l'extension DB Manager à QGIS - installer le plugin French Collector
Repérer l'ensemble des dalles : identifier les numéros des dalles
cliquer sur l'icone "identifier des entités"

d) Dans le menu QGIS :
Raster > Extraction > Contour
Régler intervalle des courbes de niveau à 1 m
Régler Enregistrer dans un fichier > mettre le nom du fichier
Cocher "Ouvrir le fichier .. après exécution de l'algorithme"
Exécuter.

Pour afficher les valeurs des courbes de niveau
Sélectionner la couche "Courbes de niveau"
Ouvrir les propriétés
Choisir Etiquettes > Etiquettes simples
Valeur > Choisir Elévation
Position > Choisir "incurvé"

Voir également cette vidéo Youtube

ou cette autre vidéo pour afficher les valeurs des courbes