GAMMA软件的InSAR处理流程课件.ppt

大纲?GAMMA简介?InSAR干涉处理流程-步骤-示例?差分干涉测量及地理编码流程-步骤-示例GAMMA简介?Source: Swiss corporation (Aktiengesellschaft - AG) founded by Dr. Charles Werner and Dr. Urs Wegmuller?Environment: Unix, Linuxor Windows?Packages:MSP,ISP,DIFF&GEO,LAT,IPTA,DISP?Object of dealing with:SAR data from 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?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 coffs 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 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 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 fltsm.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 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）InSAR-derived DEMOrbital dataorExternal DEMDEM pairInterferogram (topography+deformation)TopographicfringesRange changemapDifferential interferogramDeformation pairSimulatingOrbital dataPhase unwrappingSAR imagepairFilteringDeformationmapGeometrictransformDIFF&GEO处理步骤（1）生成转换表 : ?gc_map 960523.mli.par –LV.dem_par LV.dem dem_seg_par 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 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）相位解缠?为了减少数据量，提高解缠效率，我们把质量差的点掩模掉：?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 3160 2028 0 2?geocode_back disp 2500 neqa_rdc disp_eqa 3160 2028 0 2示例（台湾集集地震）19990715-19991028SRTM data SRTM fill datahttp://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地震的同震地表形变位移场。