LSGE LidarHD Ortho4XP Tuto,

[jasum]

Exemple d'utilisation de données LidarHD à l'approche de LSGE_Ecuvillens_Fribourg.
Personnelement je n'utilise pas de données LIDAR sur la totalité d'une tuile, mais j'utilise volontiers des LidarHD pour les approches de mes aérodrome préférés.
Contrairement à notre "IGN" qui actuelement nous propose en opendata des MNT à 1métre où seule l'hydrographie est réelement Lidar à 1mètre, le reste étant dans la plupart des zones qu'un rééchantillonnage de données classiques, Swisstopo nous offre en opendata de vrai LIDAR_HD à 0.5m et 2m sur tout son teritoire.

L'opération consiste à demander à Ortho4XP d'utiliser un fichier LSGE_WGS84.tif (haute définition) https://www.dropbox.com/s/o2d47uk75i..._LSGE.zip?dl=0 en complément pour la construction du mesh.

<<<<<<<<Non necessaire si vous utilsez le fichier en pièce jointe.>>>>>>>>>>>>>>>>>>>>>

J'ai obtenu ce fichier comme cela:
Avoir la librairie "gdal" installée et fonctionnelle (cela peut-être délicat sous Windows).
Télécharger dans un répertoire, les données LIDAR en geotif sur le site de Swisstopo: selectionner le fichier qui couvre l'aérodrome et ceux adjacents à l'Est et à l'Ouest cad 3 fichiers au total.
Depuis ce répertoire, construire un geotiff LSGE_WGS84 qui sera utilisé par Ortho4XP:
## on construit un fichier virtuel qui comprend les géotifs contenus dans ce répertoire.

gdalbuildvrt LSGE_swisstopo.vrt *.tif

## on convertit le fichier virtuel du systéme de référence suisse ESPG:2056 vers celui utilisé par Ortho4XP EPSG:4326 , on spécifie les zones de nodata à -9999 , on choisi le resampling en "cubic" (cela limite les artefacts, si non mentionné il sera en "near" plus proche voisin) , on obtient le fichier utilisable par Otho4XP:

gdalwarp -s_srs EPSG:2056 -t_srs EPSG:4326 -dstnodata -9999 -rc LSGE_swisstopo.vrt LSGE_WGS84.tif

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >> >>>>>>>>>>>>>>>>>>

Placer le fichier LSGE_WGS84.tif dans le répertoire "Elevation_data" d'Ortho4XP

Construire la tuile +46+007 avec les options suivantes:

>> Ortho4XP_Config (icône avec le tournevis) >> "Application" cleaning_level: 0

Configurer les autres options à votre convenance (ou laisser par default),
Enregister la configuration >> Write Tile Cfg (en bas à gauche)
Editer Tiles >> zOrtho4XP_+46+007 >> Ortho4XP_+46+007.cfg (sous Windows utiliser notepad ou notepad++), à la ligne "custom_dem=" indiquer les fichiers d'élévation que vous utiliserez avec un point virgule entre eux:

custon_dem= chemin_vers>>Elevation_data>>N46E007.hgt;chemin_ve rs_>>Elevation_data>>LSGE_WGS84.tif

(il y a ; entre les deux )

Enregistrer.
Vérifier la configuration de la tuile pour qu'Ortho4XP prenne en compte la config que l'on vient d'éditer:

>> Ortho4XP_Config (icône avec le tournevis) >> Load_Tile_cfg (en bas à gauche) >> Apply (en bas à droite) puis Exit pour fermer la fenêtre.

Faire étape 1 & 2, "Assemble Vector data" et "Triangulate 3D Mesh" comme d'habitude

Lorque l'étape "Triangulate 3D Mesh" est terminée (Ortho4XP a utisé N46E007.hgt)

Ouvrir >> Ortho4XP_Config (icône avec le tournevis) >> Mesh >> (mettre option iterate à 1) iterate:1 puis >> Apply >> exit

Renouveler l'étape2 "Triangulate 3D Mesh" (Ortho4XP va "affiner" avec LSGE_WGS84.tif) et continuer classiquement les autres étapes .

[hawkertempest]

merci... je me demande si je vais pas m'installer un ti linux virtuel rien pour ne plus galérer à faire tourner python, Gdal, et ortho4xp.py...
"Avoir la librairie "gdal" installée et fonctionnelle (cela peut-être délicat sous Windows)." je confirme, même sous python, il me signale "deprecated" partout ... etc etc...

[Glah_Salamanthe]

Exemple d'utilisation de données LidarHD à l'approche de LSGE_Ecuvillens_Fribourg.

jasum : Hello
Pourquoi ce soudain intérêt pour les aérodromes de Suisse ? Sont-ce de ceux que tu préfères ? ou est-ce la qualité du Lidar HD que tu ne trouverai pas ailleurs ?
LSGE, c'est là que j'ai effectué mon saut en chute libre (lâché sol 3000m , ouverture à 1000m !) bien sûr en biplace (génial)!

Je vais donc tenter (pas tout de suite) de suivre à la lettre toute ta démarche , fort bien décrite et peut-être améliorer mes terrains helvétiques, si cela ne me pose pas trop de soucis sur Mac .

Amitiés .

Glah

[jojo64]

j ou est-ce la qualité du Lidar HD que tu ne trouverai pas ailleurs ?

Si je peux te donner un avis là dessus, j'ai pas mal bricolé de données alti sur l'Europe voir même le monde, et la Suisse possède sans hésitation une des meilleures bd altimétrique ouverte, si ce n'est la meilleur. A la limite il n'y a que l'Italie qui arrive à faire aussi bien.

La seule chose que j'aurai à redire, c'est la forme de distribution des données.
Là où l'Italie ou la France propose de télécharger le pays entiers sur la base de quelques archives, 96 archives pour la France et même pas 12 pour l'Italie, la Suisse elle propose ça sous la forme de plus de 60 000 archives, à télécharger une par une à partir de 60 000 liens fournis sur un fichier .txt. Autant dire que ce n'est pas pratique du tout, il m'a fallut créer un programme pour le faire automatiquement (et ça m'a pris plus de 3h à tout télécharger/décompresser).
Mais c'est juste pour trouver un truc à redire, à part ca la Suisse est impressionnante quand à la qualité de ses données géographique! A la hauteur de la beauté de sa géographie!

[hawkertempest]

youhouhou, merci pour l'astuce, .. voici le code sous Windows : (je travaille dans la racine de S: mon DD de 16 To...)

import os
import csv
import requests

# Chemin du fichier CSV contenant les adresses internet des fichiers TIF à télécharger
csv_file = r'S:\suisse.csv'

# Boucle sur chaque ligne du fichier CSV
with open(csv_file, newline='') as csvfile:
    reader = csv.reader(csvfile)
    for row in reader:
        # Récupération de l'adresse internet de la ligne courante
        url = row[0]
        # Récupération du nom du fichier TIF dans l'adresse internet
        filename = os.path.basename(url).split('\\')[-1]
        # Téléchargement du fichier TIF à partir de l'adresse internet
        r = requests.get(url)
        with open(filename, 'wb') as f:
            f.write(r.content)​

[jasum]

Je ne comprends pas de quoi tu parles!
Code: sous Windows ?
Quel code pour faire quoi?

[hawkertempest]

en effet je n'ai pas désigné la personne à qui je répondais dans le message. En fait c'est ce que jojo64 a proposé: " la Suisse elle propose ça sous la forme de plus de 60 000 archives, à télécharger une par une à partir de 60 000 liens fournis sur un fichier .txt. Autant dire que ce n'est pas pratique du tout, il m'a fallut créer un programme pour le faire automatiquement (et ça m'a pris plus de 3h à tout télécharger/décompresser)."
ça m'a donné l'idée d'essayer et de me remettre au Python (c'est bon contre Alzheimer)
J'ai juste renommé le nom de la liste des URL des TIF du CSV (téléchargé ici : https://www.swisstopo.admin.ch/fr/ge...ht/alti3d.html) au préalable (trop long) en 'suisse.csv'

(au fait j'espère que c'est gratuit )

[jojo64]

Super, c'est à peu prêt ce que j'avais fait, mais en Java.
Comme je suis un bourin, j'avais tout codé en dure en désignant directement l'url de mon fichier dans le programme, du coup inutilisable pour quelqu'un d'autre. Je ne l'ai donc jamais partagé, et comme en plus j'ai l'impression que ça n’intéresse pas grand monde ce genre de truc.

Par contre, tu vas in fine te retrouver à devoir retaper dans du Gdal. Jasum est sympa, il t'as déjà expliqué comment faire!

[hawkertempest]

Ah oui, vous aurez tous reconnu que c'est du Python, les gamins de collège savent coder ça maintenant...
Assez bizarrement Gdal produit des erreurs à l'installation sous Python 3.11 sous Windows 10..... J'ai désinstallé Anaconda: trop usine à gaz! Python est si simple à partir du DOS dès les "path" sont paramétrés...
Mais QGIS fait le TAF pour fusionner les TIF, donc pas de souci..

[papet30]

Bonsoir en complément à la présentation complète de Jasum, vous pouvez visualiser le petit tuto, sur le sujet, réalisé par Oscar Pilote.

https://www.youtube.com/watch?v=o38f7n90Pmw&t=10s

PS : désolé avec youtube on échappe rarement à la pub, je n'y suis pour rien.

[jasum]

Salut hawkertempest

Puisque tu utilises QGIS, il y a une extension QGIS "swiss geo downloader" pour les données opendata de swisstopo.
Tu risques un banissement de ton adresse IP en voulant télécharger la totalité des 60.000 fichiers.
Il faudrait mieux diviser ta liste en morceaux raisonnables.

[hawkertempest]

ok merci de vos conseils papet30 t30 et jasum
Je vais essayer d'avoir de beaux "pics" pour survoler la Suisse (mais ça va ramer)