step_dem

StaMps最新测试版

#!/bin/tcsh -f

#set WORKDIR = `cat .workdir`
set WORKDIR = ".."
set first = `grep 'First_pixel ' $WORKDIR/master.res | gawk 'END {print $4}'`
set last = `grep 'Last_pixel ' $WORKDIR/master.res | gawk 'END {print $4}'`
@ width = ($last - $first) + 1
echo $width > ../width.txt
                                                                                         
set first_l = `grep 'First_line ' $WORKDIR/master.res | gawk 'END {print $4}'`
set last_l = `grep 'Last_line ' $WORKDIR/master.res | gawk 'END {print $4}'`
@ numlines = ($last_l - $first_l) + 1
@ numlines_ml= $numlines / 5

if (! -e $WORKDIR/dem.dorisin) then
    cp $DORIS_SCR/dem.dorisin $WORKDIR
endif                                                                                         
echo "*****************************************************" > fake_int.out
echo "*_Start_interfero:" >> fake_int.out
echo "*****************************************************" >> fake_int.out
echo "Data_output_file:                       cint.raw" >> fake_int.out
echo "Data_output_format:                     complex_real4" >> fake_int.out
echo "Flatearth correction subtracted:        no" >> fake_int.out
echo "First_line (w.r.t. original_master):   " $first_l >> fake_int.out
echo "Last_line (w.r.t. original_master):    " $last_l >> fake_int.out
echo "First_pixel (w.r.t. original_master):  " $first >> fake_int.out
echo "Last_pixel (w.r.t. original_master):   " $last >> fake_int.out
echo "Multilookfactor_azimuth_direction:      5" >> fake_int.out
echo "Multilookfactor_range_direction:        1" >> fake_int.out
echo "Number of lines (multilooked):         " $numlines_ml >> fake_int.out
echo "Number of pixels (multilooked):        " $width >> fake_int.out
echo "*****************************************************" >> fake_int.out
echo "* End_interfero:_NORMAL" >> fake_int.out
echo "*****************************************************" >> fake_int.out

get_sar_parms
cat coreg.out fake_int.out > dem.out
rm -f fake_int.out
doris $WORKDIR/dem.dorisin > step_dem.log
echo "interp_dem $width" >matbgparms.txt      
matlab -nojvm -nosplash -nodisplay < $STAMPS/matlab/batchjob.m > interp_dem.log