GAMMA软件的InSAR处理流程

基于GAMMA软件的InSAR处理流程1整理ppt大纲GAMMA简介InSAR干预处理流程-步骤-例如差分干预测量及地理编码流程-步骤-例如2整理pptGAMMA简介Source:Swisscorporation(Aktiengesellschaft-AG)foundedbyDr.CharlesWernerandDr.UrsWegmullerEnvironment:Unix,LinuxorWindowsPackages:MSP,ISP,DIFF&GEO,LAT,IPTA,DISPObjectofdealingwith:SARdatafromENVISAT,ERS1/2,Radarsat….THEWEBSITEOFGAMMA:://gamma-rs.ch3整理pptSAR干预处理流程

〔InterferometricSARProcessor-ISP〕

4整理ppt〔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整理ppt〔3〕干预纹图的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500〔4〕基线估算————基线估算方法的选择！！base_init960524.par960523.paroffintbase2——计算初始基线〔5〕平地效应的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr125006整理ppt〔6〕自适应滤波——降低噪声adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500〔7〕相位解缠相位解缠有两种方法：枝切法和最小费用流法〔mcf〕枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag--12711081最小费用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001----11–1271108117整理ppt〔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整理ppt〔10〕将斜距向高程转换为正射高程res_maphgtgrdpwr1960524.paroffrhgtrpwr17720〔读取rhgt的width为2638〕注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。rashgtrhgtrpwr126389整理ppt例如〔上海TanDEM数据〕10整理ppt11整理ppt12整理ppt13整理ppt14整理pptinit_offset_orbitoffset_pwroffset_fitinterf_SLCraspwrrasmph15整理pptbase_orbitph_slope_baserasmphcc_waveadfrasmph16整理pptUNWRAP_PARbridge17整理ppt差分干预测量及地理编码流程

〔DifferentialInterferometry&Geogoding-DIFF&GEO〕InSAR-derivedDEMOrbitaldataorExternalDEMDEMpairInterferogram(topography+deformation)TopographicfringesRangechangemap

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform18整理pptDIFF&GEO处理步骤〔1〕生成转换表:gc_map960523.mli.par–LV.dem_parLV.demdem_seg_pardem_seggeocodeeqa_rdcsim_sar3160nsim_sar2500180010create_diff_par960523.mli.par–960523.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整理ppt〔2〕将DEM从WGS84坐标系转换到SAR坐标系:geocodeneqa_rdcdem_seg3160960523.dem2500180010〔3〕模拟解缠的地形相位:phase_sim960523.paroffbase960523.demsim_unw00〔4〕平地及地形相位的剔除:sub_phaseintsim_unwdiff_pardiff1020整理ppt〔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_eqa316020280221整理ppt例如〔台湾集集地震〕19990715-1999102822整理pptSRTMdataSRTMfilldata:///srtm/23整理pp