GAMMA軟件的InSAR處理流程課件

基于GAMMA軟件的InSAR處理流程基于GAMMA軟件的InSAR處理流程1大綱GAMMA簡(jiǎn)介InSAR干涉處理流程

-步驟

-示例差分干涉測(cè)量及地理編碼流程

-步驟

-示例大綱GAMMA簡(jiǎn)介2GAMMA簡(jiǎn)介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.chGAMMA簡(jiǎn)介Source:3SAR干涉處理流程

（InterferometricSARProcessor-ISP）

SAR干涉處理流程

（InterferometricSAR4（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處理步驟（1）多視處理及顯示ISP處理步驟5（3）干涉紋圖的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基線估算————基線估算方法的選擇??！base_init960524.par960523.paroffintbase2——計(jì)算初始基線（5）平地效應(yīng)的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500（3）干涉紋圖的生成6（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–127110811（6）自適應(yīng)濾波——降低噪聲7（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）基線的精密估算——利用已知高程控制點(diǎn)數(shù)據(jù)8（10）將斜距向高程轉(zhuǎn)換為正射高程

res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（讀取rhgt的width為2638）

注：pwr1是斜距向上的強(qiáng)度文件，rpwr1是正射強(qiáng)度文件。rashgtrhgtrpwr12638（10）將斜距向高程轉(zhuǎn)換為正射高程9示例（上海TanDEM數(shù)據(jù)）示例（上海TanDEM數(shù)據(jù)）10GAMMA軟件的InSAR處理流程課件11GAMMA軟件的InSAR處理流程課件12GAMMA軟件的InSAR處理流程課件13GAMMA軟件的InSAR處理流程課件14init_offset_orbitoffset_pwroffset_fitinterf_SLCraspwrrasmphinit_offset_orbit15base_orbitph_slope_baserasmphcc_waveadfrasmphbase_orbit16UNWRAP_PARbridgeUNWRAP_PAR17差分干涉測(cè)量及地理編碼流程

（DifferentialInterferometry&Geogoding-DIFF&GEO）InSAR-derivedDEMOrbitaldataorExternalDEMDEMpairInterferogram(topography+deformation)TopographicfringesRangechangemap

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform差分干涉測(cè)量及地理編碼流程

（DifferentialI18DIFF&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)化)DIFF&GEO處理步驟（1）生成轉(zhuǎn)換表:19（2）將DEM從WGS84坐標(biāo)系轉(zhuǎn)換到SAR坐標(biāo)系:geocodeneqa_rdcdem_seg3160960523.dem2500180010（3）模擬解纏的地形相位:phase_sim960523.paroffbase960523.demsim_unw00（4）平地及地形相位的剔除:

sub_phaseintsim_unwdiff_pardiff10（2）將DEM從WGS84坐標(biāo)系轉(zhuǎn)換到SAR坐標(biāo)系:20（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_eqa3160202802（5）相位解纏21示例（臺(tái)灣集集地震）19990715-19991028示例（臺(tái)灣集集地震）19990715-1999102822SRTMdataSRTMfilldata/srtm/SRTMdata23create_dem_pargc_mapgeocodecreate_diff_parinit_offsetmoffset_pwrmgc_map_finegeocodephase_simcreate_dem_par24diff_ls_fitdiff_ls_unwdispmapgeocode_backdiff_ls_fit259.21地震的同震地表形變位移場(chǎng)9.21地震的同震地表形變位移場(chǎng)26TheEndTheEnd27基于GAMMA軟件的InSAR處理流程基于GAMMA軟件的InSAR處理流程28大綱GAMMA簡(jiǎn)介InSAR干涉處理流程

-步驟

-示例差分干涉測(cè)量及地理編碼流程

-步驟

-示例大綱GAMMA簡(jiǎn)介29GAMMA簡(jiǎn)介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.chGAMMA簡(jiǎn)介Source:30SAR干涉處理流程

（InterferometricSARProcessor-ISP）

SAR干涉處理流程

（InterferometricSAR31（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處理步驟（1）多視處理及顯示ISP處理步驟32（3）干涉紋圖的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基線估算————基線估算方法的選擇！！base_init960524.par960523.paroffintbase2——計(jì)算初始基線（5）平地效應(yīng)的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500（3）干涉紋圖的生成33（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–127110811（6）自適應(yīng)濾波——降低噪聲34（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）基線的精密估算——利用已知高程控制點(diǎn)數(shù)據(jù)35（10）將斜距向高程轉(zhuǎn)換為正射高程

res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（讀取rhgt的width為2638）

注：pwr1是斜距向上的強(qiáng)度文件，rpwr1是正射強(qiáng)度文件。rashgtrhgtrpwr12638（10）將斜距向高程轉(zhuǎn)換為正射高程36示例（上海TanDEM數(shù)據(jù)）示例（上海TanDEM數(shù)據(jù)）37GAMMA軟件的InSAR處理流程課件38GAMMA軟件的InSAR處理流程課件39GAMMA軟件的InSAR處理流程課件40GAMMA軟件的InSAR處理流程課件41init_offset_orbitoffset_pwroffset_fitinterf_SLCraspwrrasmphinit_offset_orbit42base_orbitph_slope_baserasmphcc_waveadfrasmphbase_orbit43UNWRAP_PARbridgeUNWRAP_PAR44差分干涉測(cè)量及地理編碼流程

（DifferentialInterferometry&Geogoding-DIFF&GEO）InSAR-derivedDEMOrbitaldataorExternalDEMDEMpairInterferogram(topography+deformation)TopographicfringesRangechangemap

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform差分干涉測(cè)量及地理編碼流程

（DifferentialI45DIFF&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)化)DIFF&GEO處理步驟（1）生成轉(zhuǎn)換表:46（2）將DEM從WGS84坐標(biāo)系轉(zhuǎn)換到SAR坐標(biāo)系:geocodeneqa_rdcdem_seg3160960523.dem2500180010（3）模擬解纏的地形相位:phase_sim960523.paroffbase