step_ifg

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
                                                                                           
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/interferogram.dorisin) then
    cp $DORIS_SCR/interferogram.dorisin $WORKDIR
endif
                                                                                           
echo "*****************************************************" > fake_dem.out
echo "*_Start_comp_refdem:" >> fake_dem.out
echo "*****************************************************" >> fake_dem.out
echo "Method:                                 dem gtopo30 (phase w.r.t. ellipsoid.)" >> fake_dem.out
echo "DEM source file:                        dummy" >> fake_dem.out
echo "Data_output_file:                       refdem_1li.raw" >> fake_dem.out
echo "Data_output_format:                     real4" >> fake_dem.out
echo "First_line (w.r.t. original_master):   " $first_l >> fake_dem.out
echo "Last_line (w.r.t. original_master):    " $last_l >> fake_dem.out
echo "First_pixel (w.r.t. original_master):  " $first >> fake_dem.out
echo "Last_pixel (w.r.t. original_master):   " $last >> fake_dem.out
echo "Multilookfactor_azimuth_direction:      1" >> fake_dem.out
echo "Multilookfactor_range_direction:        1" >> fake_dem.out
echo "Number of lines (multilooked):         " $numlines >> fake_dem.out
echo "Number of pixels (multilooked):        " $width >> fake_dem.out
echo "*****************************************************" >> fake_dem.out
echo "* End_comp_refdem:_NORMAL" >> fake_dem.out
echo "*****************************************************" >> fake_dem.out
cat coreg.out fake_dem.out > interferogram.out
rm -f fake_dem.out
doris $WORKDIR/interferogram.dorisin >> step_ifg.log

cpxfiddle -w$width -fcr4 -qmixed -M4/20 -osunraster -e0.2 -s1.8 cint.minrefpha.raw > cint.minrefpha_4l.ras

\rm -f cint.minrefpha.raw

cpxfiddle -w$width -fcr4 -qmixed -M4/20 -osunraster -e0.2 -s1.8 cint.minrefdem.raw > cint.minrefdem_4l.ras