#!/bin/bash # # # Author: Mahmut Arikan # # TUDelft 2007 # # Trick Doris to use radarcoded SRTM heights to gecode. iFile=$1 radarcodedDEM=$2 # Functions goodfunc(){ grep _Start_slant2h $iFile && echo slant2h info already there, pls check && exit 0; m_s=`awk '/INTERFEROGRAM RESULTFILE/{print substr($3,1,11)}' *_*.res` l0=`grep First_line *_*.res | awk 'END{print $4}'` LN=`grep Last_line *_*.res | awk 'END{print $4}'` p0=`grep First_pixel *_*.res | awk 'END{print $4}'` PN=`grep Last_pixel *_*.res | awk 'END{print $4}'` mlA=`grep Multilookfactor_azimuth_direction *_*.res | awk 'END{print $2}'` mlR=`grep Multilookfactor_range_direction *_*.res | awk 'END{print $2}'` cat << END >> $iFile ******************************************************************* *_Start_slant2h: ******************************************************************* Method: schwabisch Data_output_file: Outdata/${m_s}_S2H.float Data_output_format: real4 First_line (w.r.t. original_master): $l0 Last_line (w.r.t. original_master): $LN First_pixel (w.r.t. original_master): $p0 Last_pixel (w.r.t. original_master): $PN Multilookfactor_azimuth_direction: $mlA Multilookfactor_range_direction: $mlR Ellipsoid (name,a,b): WGS84 6.37814e+06 6.35675e+06 ******************************************************************* * End_slant2h:_NORMAL ******************************************************************* END echo $iFile is updated tail -n 35 $iFile # cd Outdata && ln -s ../refDemLP.raw ${m_s}_S2H.float && echo "File: Outdata/${m_s}_S2H.float --> refDemLP.raw" #cd Outdata && ln -s refdem_hei.raw ${m_s}_S2H.float && echo "File: Outdata/${m_s}_S2H.float --> refDemLP.raw" cd Outdata && ln -s ${radarcodedDEM} ${m_s}_S2H.float && echo "File: Outdata/${m_s}_S2H.float --> ${radarcodedDEM}" } case $# in 2) goodfunc ;; *) echo -e "Usage: ${0##*/} \n" ;; esac