#!/bin/bash

# This file executes the post-processing steps for a task after a dataset 
# has been processed by OpenDroneMap. It generates secondary outputs.
# 
# As a general rule, post-processing commands should never fail the task 
#(for example, if a point cloud could not be generated, the PotreeConverter 
# step will fail, but the script should still continue processing the rest and
# return a 0 code). The idea is to post-process as much as possible, knowing 
# that some parts might fail and that partial results should be returned in such cases.

if [ -z "$1" ]; then
	echo "Usage: $0 <projectPath>"
	exit 0
fi

# Switch to project path folder (data/<uuid>/)
script_path=$(realpath $(dirname "$0"))
cd "$script_path/../$1"
echo "Postprocessing: $(pwd)"

# Generate colored shaded relief for DTM/DSM files if available
dem_products=()
if [ -e "odm_dem/dsm.tif" ]; then dem_products=(${dem_products[@]} dsm); fi
if [ -e "odm_dem/dtm.tif" ]; then dem_products=(${dem_products[@]} dtm); fi

if hash gdaldem 2>/dev/null; then
	for dem_product in ${dem_products[@]}; do
		dem_path="odm_dem/""$dem_product"".tif"

		gdaldem color-relief $dem_path $script_path/color_relief.txt "odm_dem/""$dem_product""_colored.tif" -alpha -co ALPHA=YES
		gdaldem hillshade $dem_path "odm_dem/""$dem_product""_hillshade.tif" -z 1.0 -s 1.0 -az 315.0 -alt 45.0
		python "$script_path/hsv_merge.py" "odm_dem/""$dem_product""_colored.tif" "odm_dem/""$dem_product""_hillshade.tif" "odm_dem/""$dem_product""_colored_hillshade.tif"
	done
else
	echo "gdaldem is not installed, will skip colored hillshade generation"
fi

# Generate Tiles
g2t_options="--processes $(nproc) -z 5-21 -n -w none"
orthophoto_path="odm_orthophoto/odm_orthophoto.tif"

if [ -e "$orthophoto_path" ]; then
	python "$script_path/gdal2tiles.py" $g2t_options $orthophoto_path orthophoto_tiles

    # Check for DEM tiles also
    for dem_product in ${dem_products[@]}; do
        colored_dem_path="odm_dem/""$dem_product""_colored_hillshade.tif"
        if [ -e "$colored_dem_path" ]; then
            python "$script_path/gdal2tiles.py" $g2t_options $colored_dem_path "$dem_product""_tiles"
        else
            echo "No $dem_product found at $colored_dem_path: will skip tiling"
        fi
    done
else
	echo "No orthophoto found at $orthophoto_path: will skip tiling"
fi

# Generate point cloud (if entwine or potreeconverter is available)
pointcloud_input_path=""
for path in "odm_georeferencing/odm_georeferenced_model.laz" \
            "odm_georeferencing/odm_georeferenced_model.las" \
            "odm_filterpoints/point_cloud.ply" \
            "opensfm/depthmaps/merged.ply" \
            "smvs/smvs_dense_point_cloud.ply" \
            "mve/mve_dense_point_cloud.ply" \
            "pmvs/recon0/models/option-0000.ply"; do
    if [ -e $path ]; then
        echo "Found point cloud: $path"
        pointcloud_input_path=$path
        break
    fi
done

if [ ! -z "$pointcloud_input_path" ]; then
    # Convert the failsafe PLY point cloud to laz in odm_georeferencing 
    # if necessary, otherwise it will not get zipped
    if [ "$pointcloud_input_path" == "odm_filterpoints/point_cloud.ply" ] || [ "$pointcloud_input_path" == "opensfm/depthmaps/merged.ply" ] || [ "$pointcloud_input_path" == "pmvs/recon0/models/option-0000.ply" ]; then
        echo "Converting $pointcloud_input_path to odm_georeferencing/odm_georeferenced_model.laz, even though it's not georeferenced..."
        if hash pdal 2>/dev/null; then
            pdal translate -i "$pointcloud_input_path" -o "odm_georeferencing/odm_georeferenced_model.laz"
            if [ -e "odm_georeferencing/odm_georeferenced_model.laz" ]; then
                pointcloud_input_path="odm_georeferencing/odm_georeferenced_model.laz"
            fi
        else
            echo "Cannot find pdal, skipping conversion"
        fi
    fi
    
    if hash entwine 2>/dev/null; then
        if [ ! -e "entwine_pointcloud" ]; then
            entwine build --threads $(nproc) --tmp "entwine_pointcloud-tmp" -i "$pointcloud_input_path" -o entwine_pointcloud
        else
            echo "Entwine point cloud is already built."
        fi
        
        # Cleanup
        if [ -e "entwine_pointcloud-tmp" ]; then
            rm -fr "entwine_pointcloud-tmp"
        fi
    else
        echo "Entwine is not installed, checking if PotreeConverter is available instead..."
        if hash PotreeConverter 2>/dev/null; then
            PotreeConverter "$pointcloud_input_path" -o potree_pointcloud --overwrite -a RGB CLASSIFICATION
        else
            echo "PotreeConverter is also not installed, will skip generation of Potree point cloud"
        fi
    fi
else
    echo "Point cloud tiles will not be generated"
fi


echo "Postprocessing: done (•̀ᴗ•́)و!"
exit 0