GAMMA软件的InSAR处理流程课件

基于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——创建偏移量文件利用强度相关特性对偏移量初始估计：init_offset960524.slc960523.slc960524.par960523.paroff15

配准多项式的估计：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度达到要求，继续下面操作ISP处理步骤（1）多视处理及显示ISP处理步骤5（3）干涉纹图的生成interf_SLC960524.slc960523.slc960524.par960523.paroffpwr1pwr2int15rasmph_pwrintpwr12500（4）基线估算————基线估算方法的选择！！base_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）基线的精密估算——利用已知高程控制点数据

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）基线的精密估算——利用已知高程控制点数据8（10）将斜距向高程转换为正射高程

res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（读取rhgt的width为2638）

注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。rashgtrhgtrpwr12638（10）将斜距向高程转换为正射高程9示例（上海TanDEM数据）示例（上海TanDEM数据）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的优化)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）相位解缠为了减少数据量，提高解缠效率，我们把质量差的点掩模掉：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——创建偏移量文件利用强度相关特性对偏移量初始估计：init_offset960524.slc960523.slc960524.par960523.paroff15

配准多项式的估计：offset_pwr960524.slc960523.slc960524.par960523.paroffoffssnroffset_fitoffssnroffcoffscoffsets7.03——精度达到要求，继续下面操作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）基线的精密估算——利用已知高程控制点数据

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）基线的精密估算——利用已知高程控制点数据35（10）将斜距向高程转换为正射高程

res_maphgtgrdpwr1960524.paroffrhgtrpwr17720（读取rhgt的width为2638）

注：pwr1是斜距向上的强度文件，rpwr1是正射强度文件。rashgtrhgtrpwr12638（10）将斜距向高程转换为正射高程36示例（上海TanDEM数据）示例（上海TanDEM数据）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的优化)DIFF&GEO处理步骤（1）生成转换表:46（2）将DEM从WGS84坐标系转换到SAR坐标系:geocodeneqa_rdcdem_seg3160960523.dem2500180010（3）模拟解缠的地形相位:phase_sim960523.paroffbase