#! /bin/sh
###############################################################
# Plot pysarApp routine processing results.
# Author: Zhang Yunjun, 2017-07-23
###############################################################


## Change to 0 if you do not want to re-plot loaded dataset again
plot_key_files=1
plot_loaded_data=1
plot_loaded_data_aux=1
plot_timeseries=1
plot_geocoded_data=1
plot_the_rest=1


# Default file name
mask_file='maskTempCoh.h5'
geo_mask_file='./GEOCODE/geo_maskTempCoh.h5'

## Log File
log_file='plot_pysarApp.log'
echo "touch log file: "$log_file
touch $log_file

## Create PIC folder
if [ ! -d "PIC" ]; then
    echo 'Create PIC folder'
    mkdir PIC
fi


## Plot Key files
opt=' --dem INPUTS/geometryRadar.h5 --mask '$mask_file' -u cm '
#opt=' --dem INPUTS/geometryRadar.h5 --mask '$mask_file' -u cm --vlim -2 2'
if [ $plot_key_files -eq 1 ]; then
    view.py --nodisplay velocity.h5           $opt               | tee -a $log_file
    view.py --nodisplay temporalCoherence.h5  -c gray --vlim 0 1 | tee -a $log_file
    view.py --nodisplay maskTempCoh.h5        -c gray --vlim 0 1 | tee -a $log_file
    view.py --nodisplay INPUTS/geometryGeo.h5                    | tee -a $log_file
    view.py --nodisplay INPUTS/geometryRadar.h5                  | tee -a $log_file
fi


## Loaded Dataset
if [ $plot_loaded_data -eq 1 ]; then
    view.py --nodisplay INPUTS/ifgramStack.h5  unwrapPhase-        --zero-mask | tee -a $log_file
    view.py --nodisplay INPUTS/ifgramStack.h5  unwrapPhase_unwCor  --zero-mask | tee -a $log_file
    view.py --nodisplay INPUTS/ifgramStack.h5  coherence    --mask no   | tee -a $log_file
fi


## Auxliary Files from loaded dataset
if [ $plot_loaded_data_aux -eq 1 ]; then
    file=avgPhaseVelocity.h5; test -f $file && view.py --nodisplay  $file | tee -a $log_file
    view.py --nodisplay avgSpatialCoherence.h5 -c gray --vlim 0 1 | tee -a $log_file
    view.py --nodisplay mask.h5                -c gray --vlim 0 1 | tee -a $log_file
fi


## Time-series files
view='view.py --nodisplay --mask '$mask_file' --noaxis -u cm '
#view='view.py --nodisplay --mask '$mask_file' --noaxis -u cm --vlim -10 10 '
if [ $plot_timeseries -eq 1 ]; then
    file=timeseries.h5;                                       test -f $file && $view $file | tee -a $log_file

    file=timeseries_LODcor_ECMWF.h5;                          test -f $file && $view $file | tee -a $log_file
    file=timeseries_LODcor_ECMWF_demErr.h5;                   test -f $file && $view $file | tee -a $log_file
    file=timeseries_LODcor_ECMWF_demErr_refDate.h5;           test -f $file && $view $file | tee -a $log_file
    file=timeseries_LODcor_ECMWF_demErr_refDate_plane.h5;     test -f $file && $view $file | tee -a $log_file
    file=timeseries_LODcor_ECMWF_demErr_refDate_quadratic.h5; test -f $file && $view $file | tee -a $log_file

    file=timeseries_ECMWF.h5;                                 test -f $file && $view $file | tee -a $log_file
    file=timeseries_ECMWF_demErr.h5;                          test -f $file && $view $file | tee -a $log_file
    file=timeseries_ECMWF_demErr_refDate.h5;                  test -f $file && $view $file | tee -a $log_file
    file=timeseries_ECMWF_demErr_refDate_plane.h5;            test -f $file && $view $file | tee -a $log_file
    file=timeseries_ECMWF_demErr_refDate_quadratic.h5;        test -f $file && $view $file | tee -a $log_file

    file=timeseries_demErr.h5;                                test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_refDate.h5;                        test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_refDate_plane.h5;                  test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_refDate_quadratic.h5;              test -f $file && $view $file | tee -a $log_file

    file=timeseries_demErr.h5;                                test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_tropHgt.h5;                        test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_tropHgt_refDate.h5;                test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_tropHgt_refDate_plane.h5;          test -f $file && $view $file | tee -a $log_file
    file=timeseries_demErr_tropHgt_refDate_quadratic.h5;      test -f $file && $view $file | tee -a $log_file
fi


## Geo coordinates for UNAVCO Time-series InSAR Archive Product
if [ $plot_geocoded_data -eq 1 ]; then
    view.py --nodisplay --lalo-label ./GEOCODE/geo_maskTempCoh.h5          -c gray --vlim 0 1           | tee -a $log_file
    view.py --nodisplay --lalo-label ./GEOCODE/geo_temporalCoherence.h5    -c gray --vlim 0 1           | tee -a $log_file
    view.py --nodisplay --lalo-label ./GEOCODE/geo_velocity.h5              --mask $geo_mask_file -u cm | tee -a $log_file
    view.py --nodisplay --lalo-label ./GEOCODE/geo_timeseries_*.h5 --noaxis --mask $geo_mask_file -u cm | tee -a $log_file
fi


if [ $plot_the_rest -eq 1 ]; then
    file=velocityEcmwf.h5;    test -f $file && $view $file | tee -a $log_file
fi


## Move picture files to PIC folder
echo "Move *.png *.pdf into PIC folder"
mv *.png PIC/
mv *.pdf PIC/
mv *.kmz PIC/
mv $log_file PIC/