[Donnees] Conversion Litto3D en WGS84

Message posté par : Serge Morvan ---------------------------------------- Bonjour à tous, Je reviens vers vous pour vous remercier de vos conseils, j'ai, je pense réglé mon problème. J'avais l'habitude de mapper avec WMS par exemple de la 3D, d'où l'utilisation de l'EPSG:4326. Mais ici le problème est entièrement en 3D car il s'agit de fusionner des MNT. J'ai tout d'abord injecter le CRS EPSG:5698 dans les fichiers .asc, qui n'en avait pas, et ensuite reprojeté les données Lambert en WGS84 avec l'EPSG:4979. Reste le problème des trous dans le MNT et celui des abords, Litto3D étant livré en dalles rectangulaires, alors que les données suivent les contours de la côte. Ce n'est pas vraiment rédhibitoire pour mon projet. J'ai testé aussi l'utilisation de la RAF20 de l'IGN, mais sans gain apparent, gdal fait le travail. Voici mon script actuel : ####################### #!/bin/bash # Fonction pour parcourir l'arborescence des fichiers traverse_directory() { local directory="$1" # Boucle à travers chaque fichier/répertoire dans le répertoire courant for entry in "$directory"/*; do if [ -d "$entry" ]; then # echo "Directory: $entry" traverse_directory "$entry" # Appel récursif pour parcourir les sous-répertoires elif [ -f "$entry" ] ; then if [[ "$entry" = *"MNT_"* ]]; then #echo "File: $entry" cp $entry ./tmp/ fi fi done } #EPSG:2154 LAMBERT #EPSG:4326 WGS84 2D #EPSG:5720 IGN69 #EPSG:5698 RGF93 v1 / Lambert-93 + NGF-IGN69 height #EPSG:4979 WGS84 3D lambert2WGS84Tif(){ for i in *.asc do j=`basename $i .asc`.tif echo '' echo '$i : ' $i echo '$j : ' $j gdal_translate -of GTiff -a_srs EPSG:5698 $i $j k=`basename $j .tif`.asc echo '' echo '$j : ' $j echo '$k : ' $k gdal_translate -of AAIGrid -a_srs EPSG:5698 $j $k l=`basename $k .asc`_wgs84.tif echo '' echo '$k : ' $k echo '$l : ' $l gdalwarp -s_srs EPSG:5698 -t_srs EPSG:4979 -r bilinear -dstnodata -9999 $k $l done } # Point d'entrée du script if [ -z "$1" ]; then echo "Usage: $0 <directory>" exit 1 fi # Appel de la fonction traverse_directory avec le répertoire spécifié mkdir tmp traverse_directory "$1"; cd tmp mkdir wgs84Tif # Reprojection en WGS84 depuis Lambert 93 lambert2WGS84Tif mv *_wgs84.tif wgs84Tif # Fusion des images .tif cd wgs84Tif gdal_merge.py -o "$1".tif ./*.tif mv "$1".tif ../../ cd ../../ # Decommenter la ligne suivante pour suppression répertoires dans tmp #rm -R tmp ######################## ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=368569#p368569 Pour y répondre : donnees(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