Revision d171ca2ae8f1280e83f5bf85aceb2fb92972bbf6 authored by Zhang Yunjun on 06 June 2019, 23:40:29 UTC, committed by GitHub on 06 June 2019, 23:40:29 UTC
Move the files related with github pages from root to docs directory, for more simple structure. Add github user name (insarlab) to the page title.
1 parent 3438022
plot_smallbaselineApp.sh
#! /bin/sh
###############################################################
# Plot Results from Routine Workflow with smallbaselineApp.py
# Author: Zhang Yunjun, 2017-07-23
# Latest update: 2019-05-18
###############################################################
## 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'
dem_file='./inputs/geometryRadar.h5'
if [ ! -f $dem_file ]; then
dem_file='./inputs/geometryGeo.h5'
fi
# list of GAM models
tropo_models=(
'ECMWF'
'MERRA'
'NARR'
'ERA5'
)
## Log File
log_file='plot_smallbaselineApp.log'
touch $log_file
echo "\n\n\n\n\n" >> $log_file
echo "######################## ./plot_smallbaselineApp.sh ########################" >> $log_file
date >> $log_file
echo "##############################################################################" >> $log_file
#use "echo 'yoyoyo' | tee -a log" to output message to both screen and file.
#use "echo 'yoyoyo' >> log" to output message to file only.
## Create pic folder
if [ ! -d "pic" ]; then
echo 'Create ./pic folder'
mkdir pic
fi
## common view.py option for all files
view='view.py --nodisplay --dpi 150 --update '
## Plot Key files
opt=' --dem '$dem_file' --mask '$mask_file' -u cm '
#opt=' --dem '$dem_file' --mask '$mask_file' -u cm --vlim -2 2'
if [ $plot_key_files -eq 1 ]; then
file=velocity.h5; test -f $file && $view $file $opt >> $log_file
file=temporalCoherence.h5; test -f $file && $view $file -c gray --vlim 0 1 >> $log_file
file=maskTempCoh.h5; test -f $file && $view $file -c gray --vlim 0 1 >> $log_file
file=inputs/geometryRadar.h5; test -f $file && $view $file >> $log_file
file=inputs/geometryGeo.h5; test -f $file && $view $file >> $log_file
fi
## Loaded Dataset
if [ $plot_loaded_data -eq 1 ]; then
file=inputs/ifgramStack.h5
test -f $file && $view $file unwrapPhase- --zero-mask --wrap >> $log_file
test -f $file && $view $file unwrapPhase- --zero-mask >> $log_file
test -f $file && $view $file coherence- --mask no >> $log_file
fi
## Auxliary Files from loaded dataset
if [ $plot_loaded_data_aux -eq 1 ]; then
file=avgPhaseVelocity.h5; test -f $file && $view $file >> $log_file
file=avgSpatialCoh.h5; test -f $file && $view $file -c gray --vlim 0 1 >> $log_file
file=maskConnComp.h5; test -f $file && $view $file -c gray --vlim 0 1 >> $log_file
fi
## Time-series files
opt='--mask '$mask_file' --noaxis -u cm '
#opt='--mask '$mask_file' --noaxis -u cm --vlim -10 10 '
if [ $plot_timeseries -eq 1 ]; then
file=timeseries.h5; test -f $file && $view $file $opt >> $log_file
#LOD for Envisat
file=timeseries_LODcor.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_LODcor_ECMWF.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_LODcor_ECMWF_demErr.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_LODcor_ECMWF_ramp.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_LODcor_ECMWF_ramp_demErr.h5; test -f $file && $view $file $opt >> $log_file
#w trop delay corrections
for trop in "${tropo_models[@]}" "tropHgt"; do
file=timeseries_${trop}.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_${trop}_demErr.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_${trop}_ramp.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_${trop}_ramp_demErr.h5; test -f $file && $view $file $opt >> $log_file
done
#w/o trop delay correction
file=timeseries_ramp.h5; test -f $file && $view $file $opt >> $log_file
file=timeseries_demErr_ramp.h5; test -f $file && $view $file $opt >> $log_file
fi
## Geo coordinates for UNAVCO Time-series InSAR Archive Product
if [ $plot_geocoded_data -eq 1 ]; then
file=./geo/geo_maskTempCoh.h5; test -f $file && $view $file -c gray >> $log_file
file=./geo/geo_temporalCoherence.h5; test -f $file && $view $file -c gray >> $log_file
file=./geo/geo_velocity.h5; test -f $file && $view $file velocity >> $log_file
file=./geo/geo_timeseries_ECMWF_demErr_ramp.h5; test -f $file && $view $file --noaxis >> $log_file
file=./geo/geo_timeseries_ECMWF_demErr.h5; test -f $file && $view $file --noaxis >> $log_file
file=./geo/geo_timeseries_demErr_ramp.h5; test -f $file && $view $file --noaxis >> $log_file
file=./geo/geo_timeseries_demErr.h5; test -f $file && $view $file --noaxis >> $log_file
fi
if [ $plot_the_rest -eq 1 ]; then
for trop in "${tropo_models[@]}"; do
file=velocity${trop}.h5; test -f $file && $view $file --mask no >> $log_file
done
file=numInvIfgram.h5; test -f $file && $view $file --mask no >> $log_file
fi
## Move picture files to pic folder
echo "Move *.png/pdf/kmz files into ./pic folder."
mv *.png *.pdf *.kmz ./geo/*.kmz pic/
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...