《GAMMA软件的InSAR处理流程》由会员分享，可在线阅读，更多相关《GAMMA软件的InSAR处理流程(27页珍藏版)》请在人人文库网上搜索。

1、基于GAMMA软件的InSAR处理流程,大纲,GAMMA简介 InSAR干涉处理流程 -步骤 -示例 差分干涉测量及地理编码流程 -步骤 -示例,GAMMA简介,Source: Swiss corporation (Aktiengesellschaft - AG) founded by Dr. Charles Werner and Dr. Urs Wegmuller Environment: Unix, Linux or Windows Packages: MSP,ISP,DIFF&GEO,LAT,IPTA,DISP Object of dealing with: SAR data from 。

2、ENVISAT, ERS1/2, Radarsat . THE WEBSITE OF GAMMA: http:/www.gamma-rs.ch,SAR干涉处理流程 (Interferometric SAR Processor-ISP),(1)多视处理及显示 multi_look 960523.slc 960523.par 960523.mli 960523.mli.par 1 5 multi_look 960523.slc 960523.par 960523.mli 960523.mli.par 1 5 raspwr 960523.mli 2500 raspwr 960524.mli 2500。

3、 dis2pwr 960523.mli 960524.mli 2500 2500 (2)SLC影像偏移量估计：从影像相对于主影像 create_offset 960524.par 960523.par off 1创建偏移量文件 利用强度相关特性对偏移量初始估计： init_offset 960524.slc 960523.slc 960524.par 960523.par off 1 5 配准多项式的估计： offset_pwr 960524.slc 960523.slc 960524.par 960523.par off offs snr offset_fit offs snr off co。

4、ffs coffsets 7.0 3精度达到要求，继续下面操作,ISP处理步骤,(3)干涉纹图的生成 interf_SLC 960524.slc 960523.slc 960524.par 960523.par off pwr1 pwr2 int 1 5 rasmph_pwr int pwr1 2500 (4)基线估算基线估算方法的选择！ base_init 960524.par 960523.par off int base 2计算初始基线 (5)平地效应的去除 ph_slope_base int 960524.par off base flt 1 0 rasmph_pwr flt pwr1。

5、 2500,(6)自适应滤波降低噪声 adf flt flt_sm flt.smcc 2500 0.5 32 7 4 0 0 0.25 rasmph_pwr flt_sm pwr1 2500 (7)相位解缠 相位解缠有两种方法：枝切法和最小费用流法(mcf) 枝切法： UNWRAP_PAR off flt flt.smcc pwr1 fltsm.unw flag - - 1271 1081 最小费用流法： rascc_mask smcc pwr1 2500 1 1 0 1 1 0.7 mcf flt_sm flt.smcc mask.ras fltsm.unw 2500 1 - - - - 1。

6、 1 1271 1081 1,(8)基线的精密估算利用已知高程控制点数据 GCPs可以用命令gcp_ras从地图中提取；也可以从SAR坐标系下的DEM提取，extract_gcp。本实例中给定的gcp文件是用第二种方法得到的，extract_gcp 960523.dem off gcp 24 24 提取控制点处解缠相位： gcp_phase fltsm.unw off 960524_960523.gcp gcp_ph 基线的最小二乘估计： base_ls 960524.par off gcp_ph base 1 1 1 1 1 (9)将解缠相位转换为高程并生成一个地距文件 hgt_map fl。

7、tsm.unw 960524.par off base hgt grd 1 rashgt hgt pwr1 2500 1 1 0 1 1 160. 1. .35 -1 我们可以发现，flt_sm.ras与hgt.ras有地形上的差异，这是由于基线估算方法的不同。 disras hgt.ras 注明：此时的高程指的是斜距向上的的高程,(10)将斜距向高程转换为正射高程 res_map hgt grd pwr1 960524.par off rhgt rpwr1 7 7 20(读取rhgt的width为2638) 注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。 rashgt rhgt。

8、 rpwr1 2638,示例(上海TanDEM数据),init_offset_orbit offset_pwr offset_fit interf_SLC raspwr rasmph,base_orbit ph_slope_base rasmph cc_wave adf rasmph,UNWRAP_PAR bridge,差分干涉测量及地理编码流程 (Differential Interferometry & Geogoding -DIFF&GEO),DIFF&GEO处理步骤,(1)生成转换表 : gc_map 960523.mli.par LV.dem_par LV.dem dem_seg_p。

9、ar dem_seg geocode eqa_rdc sim_sar 3160 nsim_sar 2500 1800 1 0 create_diff_par 960523.mli.par 960523.diff_par 1 init_offsetm 960523.mli nsim_sar 960523.diff_par offset_pwrm 960523.mli nsim_sar 960523.diff_par offs snr 512 512 offsets 1 8 8 7.0 offset_fitm offs snr 960523.diff_par coffs coffsets 7.0 。

10、3 gc_map_fine eqa_rdc 3160 960523.diff_par neqa_rdc 1 (LUT的优化),(2)将DEM从WGS84坐标系转换到SAR坐标系 : geocode neqa_rdc dem_seg 3160 960523.dem 2500 1800 1 0 (3)模拟解缠的地形相位 : phase_sim 960523.par off base 960523.dem sim_unw 0 0 (4)平地及地形相位的剔除 : sub_phase int sim_unw diff_par diff 1 0,(5)相位解缠 为了减少数据量，提高解缠效率，我们把质量差的。

11、点掩模掉： rascc_mask smcc pwr1 2500 1 1 0 1 1 0.7 mcf diffsm smcc smcc_mask.ras diffsm.unw 2500 1 - - - - 1 1 - 1014 915 0 (相位参考点为(1014，915)) rasrmg diffsm.unw pwr1 2500显示解缠相位 (6)解缠相位转换到地表形变 dispmap diffsm.unw 960523.dem 960523.par off disp 1 (7)地理编码 geocode_back 960523.mli 2500 neqa_rdc 960523.mli_eqa 。

12、3160 2028 0 2 geocode_back disp 2500 neqa_rdc disp_eqa 3160 2028 0 2,示例(台湾集集地震),19990715-19991028,SRTM data SRTM fill data http:/www2.jpl.nasa.gov/srtm/,create_dem_par gc_map geocode create_diff_par init_offsetm offset_pwrm gc_map_fine geocode phase_sim,diff_ls_fit diff_ls_unw dispmap geocode_back,9.21地震的同震地表形变位移场,The End。