- 1、本文档共27页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
GAMMA软件的InSAR处理流程概要1
基于GAMMA软件的InSAR处理流程;大纲;GAMMA简介;SAR干涉处理流程(Interferometric SAR Processor-ISP);(1)多视处理及显示
multi_look 960523.slc 960523.par 960523.mli 960523.mli.par 1 5
multi_look 960523.slc 960523.par 960523.mli 960523.mli.par 1 5
raspwr 960523.mli 2500
raspwr 960524.mli 2500
dis2pwr 960523.mli 960524.mli 2500 2500
(2)SLC影像偏移量估计:从影像相对于主影像
create_offset 960524.par 960523.par off 1——创建偏移量文件
利用强度相关特性对偏移量初始估计:
init_offset 960524.slc 960523.slc 960524.par 960523.par off 1 5?
配准多项式的估计:
offset_pwr 960524.slc 960523.slc 960524.par 960523.par off offs snr
offset_fit offs snr off coffs coffsets 7.0 3——精度达到要求,继续下面操作
;(3)干涉纹图的生成
interf_SLC 960524.slc 960523.slc 960524.par 960523.par off pwr1 pwr2 int 1 5
rasmph_pwr int pwr1 2500
(4)基线估算————基线估算方法的选择!!
base_init 960524.par 960523.par off int base 2——计算初始基线
(5)平地效应的去除
ph_slope_base int 960524.par off base flt 1 0
rasmph_pwr flt pwr1 2500
;(6)自适应滤波——降低噪声
adf flt flt_sm flt.smcc 2500 0.5 32 7 4 0 0 0.25
rasmph_pwr flt_sm pwr1 2500
(7)相位解缠
相位解缠有两种方法:枝切法和最小费用流法(mcf)
枝切法:
UNWRAP_PAR off flt flt.smcc pwr1 fltsm.unw flag - - 1271 1081
最小费用流法:
rascc_mask smcc pwr1 2500 1 1 0 1 1 0.7
mcf flt_sm flt.smcc mask.ras fltsm.unw 2500 1 - - - - 1 1 – 1271 1081 1
;(8)基线的精密估算——利用已知高程控制点数据
GCPs可以用命令gcp_ras从地图中提取;也可以从SAR坐标系下的DEM提取,extract_gcp。本实例中给定的gcp文件是用第二种方法得到的,extract_gcp 960523.dem off gcp 24 24
提取控制点处解缠相位:
gcp_phase fltsm.unw off 960524_960523.gcp gcp_ph
基线的最小二乘估计:
base_ls 960524.par off gcp_ph base 1 1 1 1 1
(9)将解缠相位转换为高程并生成一个地距文件
hgt_map fltsm.unw 960524.par off base hgt grd 1
rashgt hgt pwr1 2500 1 1 0 1 1 160. 1. .35 -1
我们可以发现,flt_sm.ras与hgt.ras有地形上的差异,这是由于基线估算方法的不同。
disras hgt.ras
注明:此时的高程指的是斜距向上的的高程;(10)将斜距向高程转换为正射高程
res_map hgt grd pwr1 960524.par off rhgt rpwr1 7 7 20(读取rhgt的width为2638)
注:pwr1是斜距向上的强度文件,rpwr1是正射强度文件。
rashgt rhgt rpwr1 2638
;示例(上海TanDEM数据);;;;;;;;差分干涉测量及地理编码流程 (Differential Interferometry Geogoding -DIFFGEO);DIFFGEO处理步骤;;;示例(台湾集集地震); SRTM data SRTM fill data
/srtm/;
文档评论(0)