基于GAMMA軟件的InSAR處理流程基于GAMMA軟件的InSAR處理流程1大綱GAMMA簡介InSAR干涉處理流程

-步驟

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

-步驟

-示例大綱GAMMA簡介2GAMMA簡介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簡介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影像偏移量估計：從影像相對于主影像create_offset960524.par960523.paroff1——創(chuàng)建偏移量文件利用強度相關特性對偏移量初始估計：init_offset960524.slc960523.slc960524.par960523.paroff15

配準多項式的估計：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度達到要求，繼續(xù)下面操作ISP處理步驟（1）多視處理及顯示ISP處理步驟5（3）干涉紋圖的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基線估算————基線估算方法的選擇?。ase_init960524.par960523.paroffintbase2——計算初始基線（5）平地效應的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500（3）干涉紋圖的生成6（6）自適應濾波——降低噪聲adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500（7）相位解纏相位解纏有兩種方法：枝切法和最小費用流法（mcf）枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag--12711081最小費用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001----11–127110811（6）自適應濾波——降低噪聲7（8）基線的精密估算——利用已知高程控制點數(shù)據(jù)

GCPs可以用命令gcp_ras從地圖中提?。灰部梢詮腟AR坐標系下的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

我們可以發(fā)現(xiàn)，flt_sm.ras與hgt.ras有地形上的差異，這是由于基線估算方法的不同。

disrashgt.ras注明：此時的高程指的是斜距向上的的高程（8）基線的精密估算——利用已知高程控制點數(shù)據(jù)8（10）將斜距向高程轉換為正射高程

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

注：pwr1是斜距向上的強度文件，rpwr1是正射強度文件。rashgtrhgtrpwr12638（10）將斜距向高程轉換為正射高程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差分干涉測量及地理編碼流程

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

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform差分干涉測量及地理編碼流程

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

sub_phaseintsim_unwdiff_pardiff10（2）將DEM從WGS84坐標系轉換到SAR坐標系:20（5）相位解纏為了減少數(shù)據(jù)量，提高解纏效率，我們把質(zhì)量差的點掩模掉：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（5）相位解纏21示例（臺灣集集地震）19990715-19991028示例（臺灣集集地震）19990715-1999102822SRTMdataSRTMfilldata/srtm/SRTMdata23create_dem_pargc_mapgeocodecreate_diff_parinit_offsetmoffset_pwrmgc_map_finegeocodephase_simcreate_dem_par24diff_ls_fitdiff_ls_unwdispmapgeocode_backdiff_ls_fit259.21地震的同震地表形變位移場9.21地震的同震地表形變位移場26TheEndTheEnd27基于GAMMA軟件的InSAR處理流程基于GAMMA軟件的InSAR處理流程28大綱GAMMA簡介InSAR干涉處理流程

-步驟

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

-步驟

-示例大綱GAMMA簡介29GAMMA簡介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簡介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影像偏移量估計：從影像相對于主影像create_offset960524.par960523.paroff1——創(chuàng)建偏移量文件利用強度相關特性對偏移量初始估計：init_offset960524.slc960523.slc960524.par960523.paroff15

配準多項式的估計：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度達到要求，繼續(xù)下面操作ISP處理步驟（1）多視處理及顯示ISP處理步驟32（3）干涉紋圖的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基線估算————基線估算方法的選擇！！base_init960524.par960523.paroffintbase2——計算初始基線（5）平地效應的去除ph_slope_baseint960524.paroffbaseflt10rasmph_pwrfltpwr12500（3）干涉紋圖的生成33（6）自適應濾波——降低噪聲adffltflt_smflt.smcc25000.53274000.25rasmph_pwrflt_smpwr12500（7）相位解纏相位解纏有兩種方法：枝切法和最小費用流法（mcf）枝切法：UNWRAP_PARofffltflt.smccpwr1fltsm.unwflag--12711081最小費用流法：rascc_masksmccpwr12500110110.7mcfflt_smflt.smccmask.rasfltsm.unw25001----11–127110811（6）自適應濾波——降低噪聲34（8）基線的精密估算——利用已知高程控制點數(shù)據(jù)

GCPs可以用命令gcp_ras從地圖中提?。灰部梢詮腟AR坐標系下的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

我們可以發(fā)現(xiàn)，flt_sm.ras與hgt.ras有地形上的差異，這是由于基線估算方法的不同。

disrashgt.ras注明：此時的高程指的是斜距向上的的高程（8）基線的精密估算——利用已知高程控制點數(shù)據(jù)35（10）將斜距向高程轉換為正射高程

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

注：pwr1是斜距向上的強度文件，rpwr1是正射強度文件。rashgtrhgtrpwr12638（10）將斜距向高程轉換為正射高程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差分干涉測量及地理編碼流程

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

DifferentialinterferogramDeformationpairSimulatingOrbitaldataPhaseunwrappingSARimagepairFilteringDeformationmap

Geometrictransform差分干涉測量及地理編碼流程

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