#!/usr/bin/env bash 

###################################################################
#Author         : Joaquín M.C. Belart
#Email          : joaquin@lmi.is
###################################################################

source ~/shgen/lasfunciones.sh
holupre=holuhraun_pre_eruption_median_fx.tif
glo30=/u01/jmcb/1tdx/tdx30/GLO30_DGED_Iceland/tdx30_20131001_iceland_zmae_isn93.tif
livatna=/u01/jmcb/0adem8/dems/islandsdemv1/coarse_res_mosaics/IslandsDEMv1.0_20x20m_isn93_zmae.tif
plvatna=/u01/jmcb/0pl/vatna2021/vatna21_10x10m_isn93_align_lidar.tif

main()
{
    mkdir -p dems_isn93
    for dem in hofs/*_utm.nc
    do
        nam=$(basename $dem _latlon.nc | sed 's/raster//g')
        code=$(echo $nam | awk -F"_" '{print $2"_"$3"_"$4}' | sed -e 's/L//g' -e 's/R//g')
        date=$(grep "$code" ../tiles_cycles_dates.list | awk '{print $2}')
        gdalwarp -overwrite -r bilinear -tr 20 20 -t_srs EPSG:3057 -co compress=lzw -co predictor=2 $dem dems_isn93/vat_${date}${nam}_isn93.tif
    done
    
    #v20231209=dems_isn93/vat_20231209_007_449_270L_isn93.tif
    #v20231230=dems_isn93/vat_20231230_008_449_270L_isn93.tif
    #v20240120=dems_isn93/vat_20240120_009_449_270L_isn93.tif
    #v20240210=dems_isn93/vat_20240210_010_449_270L_isn93.tif
    #v20240301=dems_isn93/vat_20240301_011_449_270L_isn93.tif
    #v20240524=dems_isn93/vat_20240524_015_449_270L_isn93.tif
    #v20240614=dems_isn93/vat_20240614_016_449_270L_isn93.tif
    #v20240705=dems_isn93/vat_20240705_017_449_270L_isn93.tif
    #
    #v20240130=dems_isn93/vat_20240130_010_143_270L_isn93.tif
    #v20240219=dems_isn93/vat_20240219_011_143_270L_isn93.tif
    #v20240311=dems_isn93/vat_20240311_012_143_270L_isn93.tif
    #v20240401=dems_isn93/vat_20240401_013_143_270L_isn93.tif
    #v20240422=dems_isn93/vat_20240422_014_143_270L_isn93.tif
    #v20240513=dems_isn93/vat_20240513_015_143_270L_isn93.tif
    #v20240603=dems_isn93/vat_20240603_016_143_270L_isn93.tif
    #v20240624=dems_isn93/vat_20240624_017_143_270L_isn93.tif

    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20231209 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20231230 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240120 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240210 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240301 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240524 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240614 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240705 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240130 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240219 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240311 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240401 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240422 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240513 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240603 -max_iter 10 -outdir coreg_GLO30
    #dem_align.py -res max -mask_list "glaciers" -max_offset 300 -max_dz 20 -slope_lim 0.1 30 $glo30 $v20240624 -max_iter 10 -outdir coreg_GLO30

    #v20231209=coreg_GLO30/vat_20231209_007_449_270L_isn93coreg*align.tif
    #v20231230=coreg_GLO30/vat_20231230_008_449_270L_isn93coreg*align.tif
    #v20240120=coreg_GLO30/vat_20240120_009_449_270L_isn93coreg*align.tif
    #v20240210=coreg_GLO30/vat_20240210_010_449_270L_isn93coreg*align.tif
    #v20240301=coreg_GLO30/vat_20240301_011_449_270L_isn93coreg*align.tif
    #v20240524=coreg_GLO30/vat_20240524_015_449_270L_isn93coreg*align.tif
    #v20240614=coreg_GLO30/vat_20240614_016_449_270L_isn93coreg*align.tif
    #v20240705=coreg_GLO30/vat_20240705_017_449_270L_isn93coreg*align.tif
    #
    #v20240130=coreg_GLO30/vat_20240130_010_143_270L_isn93coreg*align.tif
    #v20240219=coreg_GLO30/vat_20240219_011_143_270L_isn93coreg*align.tif
    #v20240311=coreg_GLO30/vat_20240311_012_143_270L_isn93coreg*align.tif
    #v20240401=coreg_GLO30/vat_20240401_013_143_270L_isn93coreg*align.tif
    #v20240422=coreg_GLO30/vat_20240422_014_143_270L_isn93coreg*align.tif
    #v20240513=coreg_GLO30/vat_20240513_015_143_270L_isn93coreg*align.tif
    #v20240603=coreg_GLO30/vat_20240603_016_143_270L_isn93coreg*align.tif
    #v20240624=coreg_GLO30/vat_20240624_017_143_270L_isn93coreg*align.tif

    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240524 $v20231209 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20231230 $v20231209 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240120 $v20231230 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240210 $v20240120 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240301 $v20240210 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240524 $v20240301 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240614 $v20240524 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240705 $v20240614 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240130 $v20240603 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240130 $v20240219 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240219 $v20240311 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240311 $v20240401 -outdir swot_vs_swot
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240401 $v20240422 -outdir swot_vs_swot 
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240422 $v20240513 -outdir swot_vs_swot       
    #dem_align.py -mask_list "glaciers" -max_offset 100 -max_dz 10 -slope_lim 1 30 -max_iter 10 $v20240513 $v20240603 -outdir swot_vs_swot

}
main