Le portail francophone de la géomatique

Merci Patrik,

   En fait je pense que l'altitude n'est pas modifiée avec EPSG:4326

oui, c'est exactement ca - si tu ne modifies pas les Z tu auras une décalage avec les données "globale" de Cesium ION (le terrain mondiale, les bâtis de Google ou OSM etc) qui sont rélatif à l'ellipsoide. Mais ca dépend aussi de quoi tu mets sur l'import de tes données dans Cesium ION. A l'époque (il y a env 4 ans, quand je l'utiliserai beaucoup en Suède, Cesium ION) on avait l'option de mettre les z des données source relatif à l'ellipsoide ou "sea level". Le bon choix là serait "sea level" si tu fais pas de translation Z entre la géoide et l'ellipsoide. ION le fera pendant l'import

Encore, à l'époque, Cesium Ion utiliserait une rélation entre la géoide et l'ellipsoide qui n'était pas très précis, donc j'avais décidé de faire le translation Z moi-même avant l'import, en utilisant des données depuis le "IGN" suédois

Mais cela dit, si tu as un décalage avec tes propres données "francaises" avec un z déjà rélatif à la géoide, et les données sont tous importés dans ION dans la même façon, le souci c'est qqc autre. Mais bon, je ne crois pas qu'on a vraiment besoin de reprojeter 2154->4326 avant de faire l'import dans tous les cas? Si je ne me trompe pas ION le fera aussi

En plus, si le souci altimétrique et vraiment grand, possible que ION crois que les z sont des dégrés, pas des mètres

Mais là je commence de kill-gissa (expression suédois - un mec qui sait rien mais propose des solutions quand même)

Désolé si je viens tout confondre

Bonne journée!!

Patrik

Si on peut mettre une valeur constante, pourquoi utiliser des scans laser de qualité. la BD Alti sera meilleure.
[
Même en ne regardant que les zones côtières, il y a pas loin de 7 mètres d'écart entre le Pas de Calais et le Finistère. Et ne parlons pas des Alpes Maritimes, où les variations à courte distance sont loin d'être négligeables. Il y a tout ce qu'il faut sur le site de la géodésie IGN, ça s'appelle RAF20 ...https://geodesie.ign.fr/index.php?page=grilles

Dernière modification par Yves Egels (Mon 17 June 2024 16:40)

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