#!/bin/csh -f

set master = S1_20230129_033427_F3.PRM
set SNR = 30.0
set nx = 512
set ny = 3072
echo $nx
set PRF = `grep PRF $master |awk -F"=" '{print $2}'`
set SC_vel = `grep SC_vel $master|awk -F"=" '{print $2}'`
set earth_radius = `grep earth_radius $master|awk -F"=" '{print $2}'`
set SC_height = `grep SC_height $master|awk -F"=" '{print $2}'`
set rng_samp_rate = `grep rng_samp_rate $master | head -1 | awk '{print $3}'`
echo "ground velocity: " $ground_vel
set rng_size = `echo $rng_samp_rate | awk '{printf "%10.3f", 299792458.0/$1/2}' `
#########################################################################3
echo "azi pixel size", $azi_size
echo "rng pixel size", $rng_size
awk '{if ($5 > '$SNR' ) print $1,$3,$4,$5}' freq_xcorr.dat > azi1.dat
awk '{if ($3>-1 && $3<1) print $0}' azi1.dat > azi.dat
awk '{ if ($5 > '$SNR' ) print $1,$3,$2,$5}' freq_xcorr.dat > rng1.dat
awk '{if ($3>-3 && $3<3) print $0}' rng1.dat > rng.dat
###################################################################### azi
#####################################range##################################################333
set xmin = `gmt gmtinfo rng.dat -C |awk '{print $1}'`
set xmax = `gmt gmtinfo rng.dat -C |awk '{print $2}'`
set ymin = `gmt gmtinfo rng.dat -C |awk '{print $3}'`
set ymax = `gmt gmtinfo rng.dat -C |awk '{print $4}'`
set xinc = `echo $xmax $xmin $nx |awk '{printf "%d", ($1-$2)/($3-1)}'`
set yinc = `echo $ymax $ymin $ny |awk '{printf "%d", ($1-$2)/($3-1)}'`
echo "Range: xinc:"$xinc"  yinc:"$yinc
gmt blockmedian rng.dat -R$xmin/$xmax/$ymin/$ymax -I$xinc/$yinc  -Wi | awk '{print $1, $2, $3}' > rng_b.dat
proj_ra2ll_ascii.csh trans.dat rng_b.dat roff.llo
set xmin2 = `gmt gmtinfo roff.llo -C |awk '{print $1}'`
set xmax2 = `gmt gmtinfo roff.llo -C |awk '{print $2}'`
set ymin2 = `gmt gmtinfo roff.llo -C |awk '{print $3}'`
set ymax2 = `gmt gmtinfo roff.llo -C |awk '{print $4}'`
set xinc2 = `echo $xmax2 $xmin2 $nx |awk '{printf "%12.5f", ($1-$2)/($3-1)}'`
set yinc2 = `echo $ymax2 $ymin2 $ny |awk '{printf "%12.5f", ($1-$2)/($3-1)}'`
set xinc2 = `echo $xinc2 | awk '{print $1}'`
set yinc2 = `echo $yinc2 | awk '{print $1}'`
echo "Range: xinc2:"$xinc2"  yinc2:"$yinc2
gmt xyz2grd roff.llo -R$xmin2/$xmax2/$ymin2/$ymax2  -I$xinc2/$yinc2  -Groff_ll.grd
gmt grdmath roff_ll.grd $azi_size MUL = rg_offset_ll.grd
rm -f tmp.grd aoff_ll.grd aoff.grd aoff.llo azi.dat dem_grd.grd grey_tmp.cpt ps2rast* raln* ralt* s_dem.grd temp.dat topo.cpt