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

基于GAMMA软件的InSAR处理流程,胡俊2008-12-22,1,大纲,GAMMA简介InSAR干涉处理流程-步骤-示例差分干涉测量及地理编码流程-步骤-示例,2,GAMMA简介,Source:Swisscorporation(Aktiengesellschaft-AG)foundedbyDr.CharlesWernerandDr.UrsWegmullerEnvironment:Unix,LinuxorWindowsPackages:MSP,ISP,DIFF&GEO,LAT,IPTA,DISPObjectofdealingwith:SARdatafromENVISAT,ERS1/2,Radarsat.THEWEBSITEOFGAMMA:http:/www.gamma-rs.ch,3,SAR干涉处理流程（InterferometricSARProcessor-ISP）,4,（1）多视处理及显示multi_look960523.slc960523.par960523.mli960523.mli.par15multi_look960523.slc960523.par960523.mli960523.mli.par15raspwr960523.mli2500raspwr960524.mli2500dis2pwr960523.mli960524.mli25002500（2）SLC影像偏移量估计：从影像相对于主影像create_offset960524.par960523.paroff1创建偏移量文件利用强度相关特性对偏移量初始估计：init_offset960524.slc960523.slc960524.par960523.paroff15配准多项式的估计：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03精度达到要求，继续下面操作,ISP处理步骤,5,（3）干涉纹图的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基线估算基线估算方法的选择！base_init960524.par960523.paroffintbase2计算初始基线（5）平地效应的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500,6,（6）自适应滤波降低噪声adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500（7）相位解缠相位解缠有两种方法：枝切法和最小费用流法（mcf）枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag-12711081最小费用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001-11127110811,7,（8）基线的精密估算利用已知高程控制点数据GCPs可以用命令gcp_ras从地图中提取；也可以从SAR坐标系下的DEM提取，extract_gcp。本实例中给定的gcp文件是用第二种方法得到的，extract_gcp960523.demoffgcp2424提取控制点处解缠相位：gcp_phasefltsm.unwoff960524_960523.gcpgcp_ph基线的最小二乘估计：base_ls960524.paroffgcp_phbase11111（9）将解缠相位转换为高程并生成一个地距文件hgt_mapfltsm.unw960524.paroffbasehgtgrd1rashgthgtpwr1250011011160.1.35-1我们可以发现，flt_sm.ras与hgt.ras有地形上的差异，这是由于基线估算方法的不同。disrashgt.ras注明：此时的高程指的是斜距向上的的高程,8,（10）将斜距向高程转换为正射高程res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（读取rhgt的width为2638）注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。rashgtrhgtrpwr12638,9,示例（上海TanDEM数据）,10,11,12,13,14,init_offset_orbitoffset_pwroffset_fitinterf_SLCraspwrrasmph,15,base_orbitph_slope_baserasmphcc_waveadfrasmph,16,UNWRAP_PARbridge,17,差分干涉测量及地理编码流程（DifferentialInterferometry&Geogoding-DIFF&GEO）,18,DIFF&GEO处理步骤,（1）生成转换表:gc_map960523.mli.parLV.dem_parLV.demdem_seg_pardem_seggeocodeeqa_rdcsim_sar3160nsim_sar2500180010create_diff_par960523.mli.par960523.diff_par1init_offsetm960523.mlinsim_sar960523.diff_paroffset_pwrm960523.mlinsim_sar960523.diff_paroffssnr512512offsets1887.0offset_fitmoffssnr960523.diff_parcoffscoffsets7.03gc_map_fineeqa_rdc3160960523.diff_parneqa_rdc1(LUT的优化),19,（2）将DEM从WGS84坐标系转换到SAR坐标系:geocodeneqa_rdcdem_seg3160960523.dem2500180010（3）模拟解缠的地形相位:phase_sim960523.paroffbase960523.demsim_unw00（4）平地及地形相位的剔除:sub_phaseintsim_unwdiff_pardiff10,20,（5）相位解缠为了减少数据量，提高解缠效率，我们把质量差的点掩模掉：rascc_masksmccpwr12500110110.7mcfdiffsmsmccsmcc_mask.rasdiffsm.unw25001-11-10149150(相位参考点为（1014，915）)rasrmgdiffsm.unwpwr12500显示解缠相位（6）解缠相位转换到地表形变dispmapdiffsm.unw960523.dem960523.paroffdisp1（7）地理编码geocode_back960523.mli2500neqa_rdc960523.mli_eqa3160202802geocode_backdisp2500neqa_rdcdisp_eqa3160202802,21,示例（台湾集集地震）,19990715-19991028,22,SRTMdataSRTMfilldata/srtm/,23,create_d