Le portail francophone de la géomatique

Bonjour à tous,

Je partage un outil libre que je développe pour mes randos, qui pourrait intéresser ceux qui manipulent le LiDAR HD et les MNT IGN : GPXSolar.

[img]https://raw.githubusercontent.com/nico579/gpxsolar/main/screenshots/georezo_montage.jpg[/img]

Ce que ça fait — à partir d'une trace GPX + une date et une heure, il calcule pour chaque point du parcours s'il sera au soleil ou à l'ombre : un rayon est lancé vers le soleil (position via pysolar) et testé contre le relief et la végétation, avec correction de la courbure terrestre. L'animation ci-dessus, c'est la même rando à 9h / 13h / 18h — on voit l'ombre balayer le versant.

Sources d'altitude (au choix) :
- IGN LiDAR HD 0,5 m (MNT/MNS/MNH), RGE ALTI 5 m, BD ALTI 25 m — France, licence Etalab 2.0
- SRTM / Copernicus DEM — monde
- Végétation : ESA WorldCover (hauteur de canopée estimée)

Sorties : GeoTIFF de la carte d'ombre (EPSG:2154), KMZ pour Google Earth, overlay MBTiles + trace KML pour Locus Map / OsmAnd, et un CSV récap (% ensoleillé / ombragé par type d'ombre).

Quelques points techniques susceptibles de vous intéresser :
- ray-tracing vectorisé NumPy + Numba ;
- tuilage MBTiles aligné sur la grille de tuiles (façon gdal2tiles), resampling nearest car le raster est catégoriel (soleil / relief / végétation / nuit) ;
- pour le KMZ, placement de l'overlay par gx:LatLonQuad plutôt que LatLonBox : sinon la convergence des méridiens (~2° à l'est de 3°E) décale l'image de ~100 m aux coins.

Technique : Python (numpy / numba / rasterio / pyproj / shapely / pysolar), GPL v3, interface graphique (pywebview) ou ligne de commande, Windows / Linux / macOS.

→ https://github.com/nico579/gpxsolar

Projet jumeau, qui parlera peut-être plus aux archéo / prospecteurs LiDAR : lidar2map (dalles LiDAR HD IGN → ombrages multi-directionnels / SVF / LRM / RRIM → MBTiles / RMAP / Mapsforge pour OsmAnd, Locus, QGIS) → https://github.com/nico579/lidar2map

C'est du développement de loisir : retours et critiques bienvenus, en particulier sur la chaîne de reprojection et le tuilage.

Merci !