GAMMA軟件的InSAR處理流程

基于GAMMA軟件的InSAR處理流程1整理ppt大綱GAMMA簡(jiǎn)介InSAR干預(yù)處理流程-步驟-例如差分干預(yù)測(cè)量及地理編碼流程-步驟-例如2整理pptGAMMA簡(jiǎn)介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干預(yù)處理流程

〔InterferometricSARProcessor-ISP〕

4整理ppt〔1〕多視處理及顯示multi_look960523.slc960523.par960523.mli960523.mli.par15multi_look960523.slc960523.par960523.mli960523.mli.par15raspwr960523.mli2500raspwr960524.mli2500dis2pwr960523.mli960524.mli25002500〔2〕SLC影像偏移量估計(jì)：從影像相對(duì)于主影像create_offset960524.par960523.paroff1——?jiǎng)?chuàng)立偏移量文件利用強(qiáng)度相關(guān)特性對(duì)偏移量初始估計(jì)：init_offset960524.slc960523.slc960524.par960523.paroff15

配準(zhǔn)多項(xiàng)式的估計(jì)：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度到達(dá)要求，繼續(xù)下面操作ISP處理步驟5整理ppt〔3〕干預(yù)紋圖的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500〔4〕基線估算————基線估算方法的選擇！！base_init960524.par960523.paroffintbase2——計(jì)算初始基線〔5〕平地效應(yīng)的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr125006整理ppt〔6〕自適應(yīng)濾波——降低噪聲adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500〔7〕相位解纏相位解纏有兩種方法：枝切法和最小費(fèi)用流法〔mcf〕枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag--12711081最小費(fèi)用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001----11–1271108117整理ppt〔8〕基線的精密估算——利用高程控制點(diǎn)數(shù)據(jù)GCPs可以用命令gcp_ras從地圖中提??；也可以從SAR坐標(biāo)系下的DEM提取，extract_gcp。本實(shí)例中給定的gcp文件是用第二種方法得到的，extract_gcp960523.demoffgcp2424提取控制點(diǎn)處解纏相位：gcp_phasefltsm.unwoff960524_960523.gcpgcp_ph基線的最小二乘估計(jì)：base_ls960524.paroffgcp_phbase11111〔9〕將解纏相位轉(zhuǎn)換為高程并生成一個(gè)地距文件hgt_mapfltsm.unw960524.paroffbasehgtgrd1rashgthgtpwr1250011011160.1..35-1我們可以發(fā)現(xiàn)，flt_sm.ras與hgt.ras有地形上的差異，這是由于基線估算方法的不同。disrashgt.ras注明：此時(shí)的高程指的是斜距向上的的高程8整理ppt〔10〕將斜距向高程轉(zhuǎn)換為正射高程res_maphgtgrdpwr1960524.paroffrhgtrpwr17720〔讀取rhgt的width為2638〕注：pwr1是斜距向上的強(qiáng)度文件，rpwr1是正射強(qiáng)度文件。rashgtrhgtrpwr126389整理ppt例如〔上海TanDEM數(shù)據(jù)〕10整理ppt11整理ppt12整理ppt13整理ppt14整理pptinit_offset_orbitoffset_pwroffset_fitinterf_SLCraspwrrasmph15整理pptbase_orbitph_slope_baserasmphcc_waveadfrasmph16整理pptUNWRAP_PARbridge17整理ppt差分干預(yù)測(cè)量及地理編碼流程

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

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform18整理pptDIFF&GEO處理步驟〔1〕生成轉(zhuǎn)換表: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的優(yōu)化)19整理ppt〔2〕將DEM從WGS84坐標(biāo)系轉(zhuǎn)換到SAR坐標(biāo)系:geocodeneqa_rdcdem_seg3160960523.dem2500180010〔3〕模擬解纏的地形相位:phase_sim960523.paroffbase960523.demsim_unw00〔4〕平地及地形相位的剔除:sub_phaseintsim_unwdiff_pardiff1020整理ppt〔5〕相位解纏為了減少數(shù)據(jù)量，提高解纏效率，我們把質(zhì)量差的點(diǎn)掩模掉：rascc_masksmccpwr12500110110.7mcfdiffsmsmccsmcc_mask.rasdiffsm.unw25001----11-10149150(相位參考點(diǎn)為〔1014，915〕)rasrmgdiffsm.unwpwr12500——顯示解纏相位〔6〕解纏相位轉(zhuǎn)換到地表形變dispmapdiffsm.unw960523.dem960523.paroffdisp1〔7〕地理編碼geocode_back960523.mli2500neqa_rdc960523.mli_eqa3160202802geocode_backdisp2500neqa_rdcdisp_eqa316020280221整理ppt例如〔臺(tái)灣集集地震〕19990715-1999102822整理pptSRTMdataSRTMfilldata:///srtm/23整理pp