GAMMA软件的InSAR处理流程

基于GAMMA软件的InSAR处理流程大纲GAMMA简介InSAR干涉处理流程-步骤-示例差分干涉测量及地理编码流程-步骤-示例GAMMA简介Source:Swisscorporation(Aktiengesellschaft-AG)foundedbyDr.CharlesWernerandDr.UrsWegmullerEnvironment:Unix,LinuxorWindowsPackages:MSP,ISP,DIFF&GEO,LAT,IPTA,DISPObjectofdealingwith:SARdatafromENVISAT,ERS1/2,Radarsat….THEWEBSITEOFGAMMA:干涉处理流程（InterferometricSARProcessor-ISP）（1）多视处理及显示multi_look960523.slc960523.par960523.mli960523.mli.par15multi_look960523.slc960523.par960523.mli960523.mli.par15raspwr960523.mli2500raspwr960524.mli2500dis2pwr960523.mli960524.mli25002500（2）SLC影像偏移量估计：从影像相对于主影像create_offset960524.par960523.paroff1——创建偏移量文件利用强度相关特性对偏移量初始估计：init_offset960524.slc960523.slc960524.par960523.paroff15配准多项式的估计：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度达到要求，继续下面操作ISP处理步骤（3）干涉纹图的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基线估算————基线估算方法的选择！！base_init960524.par960523.paroffintbase2——计算初始基线（5）平地效应的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500（6）自适应滤波——降低噪声adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500（7）相位解缠相位解缠有两种方法：枝切法和最小费用流法（mcf）枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag--12711081最小费用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001----11–127110811（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.paroffbasehgtgrd1rashgthgtpwr1250011011160.1..35-1我们可以发现，flt_sm.ras与hgt.ras有地形上的差异，这是由于基线估算方法的不同。disrashgt.ras注明：此时的高程指的是斜距向上的的高程（10）将斜距向高程转换为正射高程res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（读取rhgt的width为2638）注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。rashgtrhgtrpwr12638示例（上海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差分干涉测量及地理编码流程（DifferentialInterferometry&Geogoding-DIFF&GEO）InSAR-derivedDEMOrbitaldataorExternalDEMDEMpairInterferogram(topography+deformation)TopographicfringesRangechangemapDifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmapGeometrictransformDIFF&GEO处理步骤（1）生成转换表:gc_map960523.mli.par–LV.dem_parLV.demdem_seg_pardem_seggeocodeeqa_rdcsim_sar3160nsim_sar2500180010create_diff_par960523.mli.par–960523.diff_par1init_offsetm960523.mlinsim_sar960523.diff_paroffset_pwrm960523.mlinsim_sar960523.diff_paroffssnr512512offsets1887.0offset_fitmoffssnr960523.diff_parcoffscoffsets7.03gc_map_fineeqa_rdc3160960523.diff_parneqa_rdc1(LUT的优化)（2）将DEM从WGS84坐标系转换到SAR坐标系:geocodeneqa_rdcdem_seg3160960523.dem2500180010（3）模拟解缠的地形相位:phase_sim960523.paroffbase960523.demsim_unw00（4）平地及地形相位的剔除:sub_phaseintsim_unwdiff_pardiff10（5）相位解缠为了减少数据量，提高解缠效率，我们把质量差的点掩模掉：rascc_masksmccpwr12500110110.7mcfdiffsmsmccsmcc_mask.rasdiffsm.unw25001----11-10149150(相位参考点为（1014，915）)rasrmgdiffsm.unwpwr12500——显示解缠相位（6）解缠相位转换到地表形变dispmapdiffsm.unw960523.dem960523.paroffdisp1（7）地理编码geocode_back960523.mli2500neqa_rdc960523.mli_eqa3160202802geocode_backdisp2500neqa_rdcdisp_eqa3160202802示例（台湾集集地震）19990715-19991028SRTMdataSRTMfilldata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_back9.21地震的同震地表形变位移场