make_small_baselines

StaMps最新测试版

#!/bin/tcsh -f
#
#
# Andy Hooper, Jul 2006
#

set WORKDIR=$cwd

if (! -e SMALL_BASELINES) then
    mkdir SMALL_BASELINES
endif

set SB_DIR = $cwd/SMALL_BASELINES

if (! -e $SB_DIR/filtrange.dorisin) then
    cp $DORIS_SCR/filtrange.dorisin $SB_DIR
endif
if (! -e $SB_DIR/filtazi.dorisin) then
    cp $DORIS_SCR/filtazi.dorisin $SB_DIR
endif

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
cd $WORKDIR
set masterdate = `\ls *_crop.slc | gawk '{print substr($0,1,8)}'` 

#cp -f $WORKDIR/small_baselines.list $SB_DIR
set i = 0

foreach slave(`cat $WORKDIR/small_baselines.list`)
   echo $slave
   if ($i == 1) then
      set IFG_DIR = $SB_DIR/$master"_"$slave
      if (! -e $IFG_DIR) then
          mkdir $IFG_DIR
      endif
      if ($master == $masterdate) then
#          cp $WORKDIR/$slave/cint.minrefdem.raw $IFG_DIR
          cp $WORKDIR/master.res  $IFG_DIR/master.res
      else
          sed -n '1,/_Start_crop:/p' $WORKDIR/$master/slave.res > $IFG_DIR/tmp_master.res
          sed -e '/resample/d' $IFG_DIR/tmp_master.res > $IFG_DIR/master.res
          rm -f $IFG_DIR/tmp_master.res
          echo "******************************************************************* " >> $IFG_DIR/master.res 
          echo "Data_output_file: " $WORKDIR"/"$master"/"slave_res.slc >> $IFG_DIR/master.res 
          echo "Data_output_format:         complex_real4 " >> $IFG_DIR/master.res 
          grep '(w.r.t. original_image):' $WORKDIR/master.res >> $IFG_DIR/master.res
          echo "******************************************************************* " >> $IFG_DIR/master.res 
          echo "* End_crop:_NORMAL" >> $IFG_DIR/master.res 
          echo "******************************************************************* " >> $IFG_DIR/master.res 
          sed -n '/_Start_precise_orbits/,/End_precise_orbits/p' $WORKDIR/$master/slave.res >> $IFG_DIR/master.res
      endif
      if ($slave == $masterdate) then
#          cpxsum $WORKDIR/$master/cint.minrefdem.raw $WORKDIR/$master/cint.minrefdem.raw tmp_zeros.raw $width cr4 -1 1
#          cpxsum tmp_zeros.raw $WORKDIR/$master/cint.minrefdem.raw $IFG_DIR/cint.minrefdem.raw $width cr4 -1 1 
#          rm -f tmp_zeros.raw
          cp $WORKDIR/master.res $IFG_DIR/master.res
          #echo "s#slave_res.slc#"$WORKDIR"/"$master"/slave_res.slc#" > $WORKDIR/tmp_sed.in
          #sed -f $WORKDIR/tmp_sed.in $WORKDIR/$master/slave.res > $IFG_DIR/slave.res
          set pslave = $master
      else
          set pslave = $slave
      endif
          sed -n '1,/_Start_crop:/p' $WORKDIR/$pslave/slave.res > $IFG_DIR/slave.res
          echo "******************************************************************* " >> $IFG_DIR/slave.res 
          echo "Data_output_file: " $WORKDIR"/"$pslave"/"slave_res.slc >> $IFG_DIR/slave.res 
          echo "Data_output_format:         complex_real4 " >> $IFG_DIR/slave.res 
          grep '(w.r.t. original_image):' $WORKDIR/master.res >> $IFG_DIR/slave.res
          echo "******************************************************************* " >> $IFG_DIR/slave.res 
          #echo "* End_crop:_NORMAL" >> $IFG_DIR/slave.res 
          #echo "******************************************************************* " >> $IFG_DIR/slave.res 
          #sed -n '/_Start_precise_orbits/,/End_resample/p' $WORKDIR/$pslave/slave.res > $IFG_DIR/tmp_slave.res
          sed -n '/End_crop:/,/End_resample/p' $WORKDIR/$pslave/slave.res > $IFG_DIR/tmp_slave.res
          echo "s#slave_res.slc#slave_azifilt.slc#" > $WORKDIR/tmp_sed.in
          sed -f $WORKDIR/tmp_sed.in $IFG_DIR/tmp_slave.res >> $IFG_DIR/slave.res 
          rm -f $IFG_DIR/tmp_slave.res
          cp $WORKDIR/$pslave/coreg.out $IFG_DIR/interferogram.out
          #echo "s#slave_res.slc#"$WORKDIR"/"$slave"/slave_res.slc#" > $WORKDIR/tmp_sed.in
          #sed -f $WORKDIR/tmp_sed.in $WORKDIR/$slave/slave.res > $IFG_DIR/slave.res
      cd $IFG_DIR 
          #sed -n '/Estimated_coefficients_flat/,/End/p' $WORKDIR/$master/interferogram.out > tmp_coeffs.in
          #tail +2 tmp_coeffs.in > tmp_coeffs2.in
          #head -9 tmp_coeffs2.in > master_coeffs.in
          #sed -n '/Estimated_coefficients_flat/,/End/p' $WORKDIR/$slave/interferogram.out > tmp_coeffs.in
          #tail +2 tmp_coeffs.in > tmp_coeffs2.in
          #head -9 tmp_coeffs2.in > slave_coeffs.in
          #matlab -nojvm -nodisplay -nosplash < $STAMPS/matlab/subtract_coeffs.m > subtract_coeffs.log     
     
          #echo "******************************************************************* " >> $IFG_DIR/interferogram.out 
          #echo "*_Start_comp_refphase:" >> $IFG_DIR/interferogram.out 
          #echo "******************************************************************* " >> $IFG_DIR/interferogram.out 
          #echo "Degree_flat:    3" >> $IFG_DIR/interferogram.out 
          #echo "Estimated_coefficients_flatearth:" >> $IFG_DIR/interferogram.out 
          #cat new_coeffs.out >> $IFG_DIR/interferogram.out 
          #echo "******************************************************************* " >> $IFG_DIR/interferogram.out 
          #echo "* End_comp_refphase:_NORMAL" >> $IFG_DIR/interferogram.out 
          #echo "******************************************************************* " >> $IFG_DIR/interferogram.out 


          doris $SB_DIR/filtazi.dorisin
          doris $SB_DIR/filtrange.dorisin
          

#cp interferogram.out subtrrefpha.out
          #sed -n '/*_Start_comp_refphase/,/End/p' $WORKDIR/$master/interferogram.out > tmp_coeffs
          #sed -e 's/-/+/' -e 's/ [^ ]\./-&/' -e 's/- /-/' tmp_coeffs >> subtrrefpha.out
          #doris $DORIS_SCR/subtrrefpha.dorisin
          #mv -f cint.minrefpha.raw cint.raw
          #cp interferogram.out subtrrefpha.out
          #sed -n '/*_Start_comp_refphase/,/End/p' $WORKDIR/$slave/interferogram.out >> $IFG_DIR/subtrrefpha.out 
          #doris $DORIS_SCR/subtrrefpha.dorisin
          if ($master == $masterdate) then
              cpxsum cint.raw $WORKDIR/$slave/cint.raw tmp_cint3.raw $width cr4 -1 1
              cpxsum tmp_cint3.raw $WORKDIR/$slave/cint.minrefdem.raw cint.minrefdem.raw $width cr4 1 1
          else if ($slave == $masterdate) then
              cpxsum $WORKDIR/$master/cint.raw cint.raw tmp_cint.raw $width cr4 -1 1
              cpxsum tmp_cint.raw $WORKDIR/$master/cint.minrefdem.raw cint.minrefdem.raw $width cr4 -1 1
          else
         
              cpxsum cint.raw $WORKDIR/$master/cint.raw tmp_cint.raw $width cr4 1 1
              cpxsum tmp_cint.raw $WORKDIR/$master/cint.minrefdem.raw tmp_cint2.raw $width cr4 -1 1
              cpxsum tmp_cint2.raw $WORKDIR/$slave/cint.raw tmp_cint3.raw $width cr4 -1 1
              cpxsum tmp_cint3.raw $WORKDIR/$slave/cint.minrefdem.raw cint.minrefdem.raw $width cr4 1 1
          endif
          rm -f tmp_cint*.raw
      rm -f cint.raw
      rm -f *azifilt.slc
      cpxfiddle -w$width -fcr4 -qmixed -M4/20 -osunraster -e0.2 -s1.8 cint.minrefdem.raw > cint.minrefdem_4l.ras
      cpxfiddle -w$width -fcr4 -qnormal -M4/20 -ofloat cint.minrefdem.raw > cint.minrefdem_4l.raw
       set i = 0
   else
       set i = 1
       set master = $slave
   endif
end
rm -f tmp_sed.in
cd $WORKDIR