if [ -z "$1" ]; then

echo "Usage: $0 <projectPath>"

exit 0

script_path=$(realpath $(dirname "$0"))

cd "$script_path/../$1"

echo "Postprocessing: $(pwd)"

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

echo "gdaldem is not installed, will skip colored hillshade generation"

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

echo "No orthophoto found at $orthophoto_path: will skip tiling"

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

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

echo "Point cloud tiles will not be generated"

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

exit 0