GAMMA地理编码教程

#定义路径

'p='文件路径',d='输出路径''

#第一步是进行多视，这一步需要输入slc影像和它的参数文件，生成它的多视影像和参数文件，命令末尾输入多视系数

'multi_look &p/20180108.slc.par $d/20180108.mli$d/20180108.mli.par 5 1 '

#保存并查看多视后的影像,末尾输入多视影像的宽度

'raspwr $p/20180108.mli 4306 '

#创建空DEM参数文件

'create_dem_par $d/20180108.dem_par$d/20180108.mli.par'

#生成初始查找表和地图结构下的强度模拟图，参数输入重采样因数，范围较大的DEM放在前面

'gc_map $d/20180108.mli.par -$p/shanxi.dem_par $p/shanxi.dem$p/20180108.dem_par $p/20180108.dem$d/20180108.rough.utm_to_rdc 3 2 $d/20180108.utm.sim_sar'

**#把模拟的 SAR 影像从地图结构转换到雷达结构,我们必须知道模拟地图结构 SAR 强度影像的列数(4678),将要获得的雷达结构影像的维数(4306 4175)

'geocode $d/20180108.rough.utm_to_rdc$d/20180108.utm.sim_sar 4678 $d/20180108.sim_sar 4306 4175 1 0 '

#生成 DIFF/GEO 参数文件

'create_diff_par $p/20180108.mli.par -$d/20180108.diff_par 1 '

#计算初始偏移值

'init_offsetm $d/20180108.mli$d/20180108.sim_sar $d/20180108.diff_par'

#计算局部偏移值,调整窗口大小 'offset_pwrm $d/20180108.mli$d/20180108.sim_sar $d/20180108.diff_par$d/offs $d/ccp 256 256 offsets 1 8 8 0.15 '

'offset_fitm $d/offs$d/ccp $d/20180108.diff_par$d/coffs $d/coffsets 0.15 ' #观察offset_fitm命令生成的结果中，range和azimuth是否小于0.2，若大于0.2则要调整窗口的大小，重复这两个命令直到小于0.2

#精化查找表

'gc_map_fine $d/20180108.rough.utm_to_rdc 4678$d/20180108.diff_par $d/20180108.utm_to_rdc 0 '

#向后编码，检查地理编码后的影像是否正确

'geocode_back $d/20180108.mli 4306$d/20180108.utm_to_rdc $d/20180108.mli.utm 4678 3128 ' #用程序raspwr可以查看该影像