Revision f24076bbf4a4a4af8e90100f344e10fb237e9b3b authored by Zhang Yunjun on 08 November 2022, 07:17:03 UTC, committed by GitHub on 08 November 2022, 07:17:03 UTC
+ .circleci/config.yml update: - use $BASH_ENV to set and share environment variable (PATH) among multiple run steps - rename to workflow/job to "unit-n-workflow-tests" + cli/save_gmt.py: fix a typo in the module import + cli/load_data: check `-t smallbaselineApp.cfg` existence and print out error msg + tsview: show the reference index/date info on the slider as the title
1 parent 64ffdef
plot_smallbaselineApp.sh
#! /bin/sh
###############################################################
# Plot Results from Routine Workflow with smallbaselineApp.py
# Author: Zhang Yunjun, 2017-07-23
# Latest update: 2021-03-08
###############################################################
# Update the date above to enable auto copyover/overwrite
## 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
## Log File
log_file='plot_smallbaselineApp.log'
touch $log_file
printf "\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 && h5ls $file/unwrapPhase && $view $file unwrapPhase- --zero-mask --wrap -c cmy >> $log_file
test -f $file && h5ls $file/unwrapPhase && $view $file unwrapPhase- --zero-mask >> $log_file
test -f $file && h5ls $file/coherence && $view $file coherence- --mask no -v 0 1 >> $log_file
test -f $file && h5ls $file/connectComponent && $view $file connectComponent- --mask no >> $log_file
# phase-unwrapping error correction
for dset in 'unwrapPhase_bridging' 'unwrapPhase_phaseClosure' 'unwrapPhase_bridging_phaseClosure'; do
test -f $file && h5ls $file/$dset && $view $file $dset- --zero-mask >> $log_file
done
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 --wrap --wrap-range -10 10 '
if [ $plot_timeseries -eq 1 ]; then
file='timeseries.h5'; test -f $file && $view $file $opt >> $log_file
find . -name 'timeseries_*.h5' -exec $view {} $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
find . -name './geo/geo_timeseries_*.h5' -exec $view {} $opt >> $log_file \;
fi
if [ $plot_the_rest -eq 1 ]; then
for tropo in 'ERA5' 'ERAI' 'ECMWF' 'MERRA' 'NARR'; do
file='velocity'${tropo}'.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/copy picture files to pic folder
echo "Copy *.txt files into ./pic folder."
cp *.txt pic/
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 ...