step_geo

来自「StaMps最新测试版」· 代码 · 共 41 行

#!/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
                                                                                           
echo "*****************************************************" > fake_geo.out
echo "*_Start_slant2h:" >> fake_geo.out
echo "*****************************************************" >> fake_geo.out
echo "Method:                                 schwabisch" >> fake_geo.out
echo "Data_output_file:                       dem_radar_i.raw" >> fake_geo.out
echo "Data_output_format:                     real4" >> fake_geo.out
echo "First_line (w.r.t. original_master):   " $first_l >> fake_geo.out
echo "Last_line (w.r.t. original_master):    " $last_l >> fake_geo.out
echo "First_pixel (w.r.t. original_master):  " $first >> fake_geo.out
echo "Last_pixel (w.r.t. original_master):   " $last >> fake_geo.out
echo "Multilookfactor_azimuth_direction:      1" >> fake_geo.out
echo "Multilookfactor_range_direction:        1" >> fake_geo.out
echo "Ellipsoid (name,a,b):                   WGS84 6.37814e+06 6.35675e+06" >> fake_geo.out
echo "*****************************************************" >> fake_geo.out
echo "* End_slant2h:_NORMAL" >> fake_geo.out
echo "*****************************************************" >> fake_geo.out

echo "interp_dem" $width "dem_radar.raw" "dem_radar_i.raw"  >matbgparms.txt
matlab -nosplash -nodisplay < $STAMPS/matlab/batchjob.m > step_geo.log
cat interferogram.out fake_geo.out > geocode.out
rm -f fake_geo.out
doris $WORKDIR/geocode.dorisin >> step_geo.log
mv lon.raw $WORKDIR
mv lat.raw $WORKDIR
mv dem_radar_i.raw $WORKDIR