GAMMA提DEM脚本（ERS）

##1、读取ERS1/2数据
 
par_ESA_ERS E:\11GAMMA_ERS_Demo\ERS\19960415\LEA_01.001 19960415.slc.par E:\11GAMMA_ERS_Demo\ERS\19960415\DAT_01.001 19960415.slc
par_ESA_ERS E:\11GAMMA_ERS_Demo\ERS\19960416\LEA_01.001 19960416.slc.par E:\11GAMMA_ERS_Demo\ERS\19960416\DAT_01.001 19960416.slc

##2、ERS轨道加密

##添加30 1后估计基线不对？？
##DELFT_vec2 19960415.slc.par E:\11GAMMA_ERS_Demo\ERS\ODR_VOR\ODR\ERS1 30 1
##DELFT_vec2 19960416.slc.par E:\11GAMMA_ERS_Demo\ERS\ODR_VOR\ODR\ERS2 30 1

DELFT_vec2 19960415.slc.par E:\11GAMMA_ERS_Demo\ERS\ODR_VOR\ODR\ERS1
DELFT_vec2 19960416.slc.par E:\11GAMMA_ERS_Demo\ERS\ODR_VOR\ODR\ERS2


##3、初始偏移量估计

create_offset 19960415.slc.par 19960416.slc.par 19960415_19960416.off 1
init_offset_orbit 19960415.slc.par 19960416.slc.par 19960415_19960416.off
     
init_offset 19960415.slc 19960416.slc 19960415.slc.par 19960416.slc.par 19960415_19960416.off 1 5
init_offset 19960415.slc 19960416.slc 19960415.slc.par 19960416.slc.par 19960415_19960416.off 1 1

##4、偏移多项式精估计

offset_pwr 19960415.slc 19960416.slc 19960415.slc.par 19960416.slc.par 19960415_19960416.off offs snr 128 128 offsets 2 32 32 7.0
offset_fit offs snr 19960415_19960416.off coffs coffsets 7.0 3 0

##5、计算干涉图

#重采样
SLC_interp 19960416.slc 19960415.slc.par 19960415.slc.par 19960415_19960416.off 19960416.rslc 19960416.rslc.par

#生成干涉图
SLC_intf 19960415.slc 19960416.rslc 19960415.slc.par 19960416.rslc.par 19960415_19960416.off 19960415_19960416.int 1 5 - - 1 1


##生成干涉缩略图
rasmph 19960415_19960416.int 4900 - - - - - - - 19960415_19960416.int.bmp


#多视
multi_look 19960415.slc 19960415.slc.par 19960415.mli 19960415.mli.par 1 5
multi_look 19960416.slc 19960416.slc.par 19960416.mli 19960416.mli.par 1 5
multi_look 19960416.rslc 19960416.rslc.par 19960416.rmli 19960416.rmli.par 1 5

##5、干涉基线估计
# remove flat phase

  #******  2  1024  1024  用轨道信息估计平行基线，条纹频率估计垂直基线，FFT方法的窗口大小为1024*1024   ******
#  base_init 19960415.slc.par 19960416.slc.par 19960415_19960416.off 19960415_19960416.int 19960415_19960416.base 2 1024 1024

#考虑到ERS所覆盖的地区为山区，不适合FFT条纹方法，加上ERS轨道信息已经被精化，所以应采用轨道信息来估计垂直基线与平行基线
base_init 19960415.slc.par 19960416.slc.par 19960415_19960416.off 19960415_19960416.int 19960415_19960416.base 0

base_perp 19960415_19960416.base 19960415.slc.par 19960415_19960416.off > 219960415_19960416.base.perp

##6、去除平地相位
ph_slope_base 19960415_19960416.int 19960415.slc.par 19960415_19960416.off 19960415_19960416.base 19960415_19960416.flt


##保存为.ras格式，再用ENVI打开即可 ####################################################
rasmph_pwr 19960415_19960416.flt 19960415.mli 4900 - - - 4 4 - - - 19960415_19960416.flt.ras

##保存为.bmp格式，直接打开 ####################################################
rasmph_pwr 19960415_19960416.flt 19960415.mli 4900 - - - 4 4 - - - 19960415_19960416.flt.bmp

##7、求相干图
cc_wave 19960415_19960416.flt 19960415.mli 19960416.rmli 19960415_19960416.cc 4900 5 5 1

## 生成ras格式图片###

rascc 19960415_19960416.cc - 4900 - - - 4 4 - - - - - 19960415_19960416.cc.ras
rascc 19960415_19960416.cc - 4900 - - - 4 4 - - - - - 19960415_19960416.cc.bmp

##8、干涉图滤波

adf 19960415_19960416.flt 19960415_19960416.flt_sm 19960415_19960416.smcc 4900 0.5 32 7 8 0 0 0.25

##生成Bmp
rasmph 19960415_19960416.flt_sm 4900 - - 4 4 - - - 19960415_19960416.flt_sm.bmp
rasmph_pwr 19960415_19960416.flt_sm 19960415.mli 4900 - - - 4 4 - - - 19960415_19960416.flt_sm.bmp

##9、相位解缠
#生成掩膜
rascc_mask 19960415_19960416.cc 19960415.mli 4900 1 1 0 1 1 0.3 - - - - - 19960415_19960416.mask.ras
#rascc_mask 19960415_19960416.cc 19960415.mli 4900 1 1 0 1 1 0.3 - - - - - 19960415_19960416.mask.bmp

rascc_mask_thinning 19960415_19960416.mask.ras 19960415_19960416.cc 4900 19960415_19960416.mask_thinned.ras 3 0.3 0.4 0.5
#rascc_mask_thinning 19960415_19960416.mask.ras 19960415_19960416.cc 4900 19960415_19960416.mask_thinned.bmp 3 0.3 0.4 0.5


##注意：块个数设置为4*4,不然内存溢出#####

#有掩膜
mcf 19960415_19960416.flt_sm 19960415_19960416.cc 19960415_19960416.mask.ras 19960415_19960416.flt_sm.mcf.unw0 4900 1 9 18 - - 4 4 - - -

#无掩膜
mcf 19960415_19960416.flt_sm 19960415_19960416.cc - 19960415_19960416.flt_sm.mcf.unw 4900 1 9 18 - - 4 4 - 1200 1500


#生成ras,bmp###
rasrmg 19960415_19960416.flt_sm.mcf.unw 19960415.mli 4900 - - - 4 4 - - - - - 19960415_19960416.flt_sm.mcf.unw.ras

rasrmg 19960415_19960416.flt_sm.mcf.unw0 19960415.mli 4900 - - - 4 4 - - - - - 19960415_19960416.flt_sm.mcf.unw.bmp


interp_ad 19960415_19960416.flt_sm.mcf.unw0 19960415_19960416.flt_sm.mcf.unw0_interp 4900 32 8 16 2

rasrmg 19960415_19960416.flt_sm.mcf.unw0_interp 19960415.mli 4900 - - - 4 4 - - - - - 19960415_19960416.flt_sm.mcf.unw0_interp.ras
rasrmg 19960415_19960416.flt_sm.mcf.unw0_interp 19960415.mli 4900 - - - 4 4 - - - - - 19960415_19960416.flt_sm.mcf.unw0_interp.bmp

unw_model 19960415_19960416.flt_sm 19960415_19960416.flt_sm.mcf.unw0_interp 19960415_19960416.flt_sm.mcf.unw_mod 4900 - - -

rasrmg 19960415_19960416.flt_sm.mcf.unw_mod 19960415.mli 4900 - - - 4 4 - - - - - 19960415_19960416.flt_sm.mcf.unw_mod.bmp


##10、干涉估计高程和地距
hgt_map 19960415_19960416.flt_sm.mcf.unw_mod 19960415.slc.par 19960415_19960416.off 19960415_19960416.base 19960415_19960416.hgt 19960415_19960416.grd 1 0 - 19960416.rslc.par


rashgt 19960415_19960416.hgt 19960415.mli 4900 - - - 4 4 - - - - 19960415_19960416.hgt.bmp

res_map 19960415_19960416.hgt 19960415_19960416.grd 19960415.mli 19960415.slc.par 19960415_19960416.off 19960415_19960416.rhgt 19960415.grd 7 7 20

rashgt 19960415_19960416.rhgt 19960415.grd 5042 - - - 4 4 - - - - 19960415_19960416.rhgt.ras
rashgt 19960415_19960416.rhgt 19960415.grd 5042 - - - 4 4 - - - - 19960415_19960416.rhgt.bmp


rasshd 19960415_19960416.rhgt 5042 20 20 - - 4 4 - - - 19960415_19960416.rhgt_shd.ras
rasshd 19960415_19960416.rhgt 5042 20 20 - - 4 4 - - - 19960415_19960416.rhgt_shd.bmp



##11、地理编码
(利用19960415主影像地理编码生成的查找表，对hgt进行地理编码。并赋予投影)   （5826 6215为地理编码后的行列号，也是查找表的行列号）

geocode_back 19960415_19960416.hgt 4900 Tibet_cover.rough.utm_to_rdc 19960415_19960416.hgt.utm 5826 6215 2 0