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