Insar处理过程及原理详解.pptVIP

Insar处理过程及原理详解 基于轨道信息的粗配准 流程图 计算偏移量多项式系数的精配准 干涉图和强度图的生成 从光盘读取SLC数据 选取主影像和从影像 ……… 读取数据及选取主从影像 原始光盘数据 解压缩并获取所需头文件及影像文件 例：SO_000012648_0002_1.tar.gz 例：头文件-- TSX1_SAR__SSC______SM_S_SRA221744221752.xml 影像文件-- IMAGE_HH_SRA_strip_012.cos 选取成像时间居中的影像为主影像，其余为从影像 例：2009.06.23的影像为主影像，其余为从影像，目的是为了降低影像配准时时间上可能产生的影响 读取数据 例slcslc.par 影像数据参数说明 29795*0.00026849613s=7.999842s ≈ 8s 80272-80264 = 8s dt dR dD R0 影像数据参数说明 41° Rnear = 652144.4173m Rcenter Rfar = 672850.1688m · (672850.1688-652144.4173)/sin41° ≈31km 41° 29795* 1.896414 ≈ 56km 给定行列号（r, c），则可计算出该像元的斜距R与成像时间t: R=Rnear+c×dR t=t0+r×dt 影像数据参数说明 轨道数据 轨道点个数 影像数据参数说明 12 11 10 9 8 7 6 5 4 3 2 1 12 11 10 9 8 7 6 5 4 3 2 1 (Xp, Yp, Zp) (Xm, Ym, Zm), (Vxm,Vym, Vzm) (Xs, Ys, Zs), (Vxs,Vys, Vzs) 6 7 Track 1 Track 2 基于轨道信息粗配准原理 X=a0+a1t+a2t2+a3t3 Y=b0+b1t+b2t2+b3t3 Z=c0+c1t+c2t2+c3t3 Vx=d0+d1t+d2t2+d3t3 Vy=e0+e1t+e2t2+e3t3 Vz=f0+f1t+f2t2+f3t3 12 11 10 9 8 7 6 5 4 3 2 1 12 11 10 9 8 7 6 5 4 3 2 1 (Xp, Yp, Zp) (Xm, Ym, Zm), (Vxm,Vym, Vzm) (Xs, Ys, Zs), (Vxs,Vys, Vzs) Track 1 Track 2 h 12 11 10 9 8 7 6 5 4 3 2 1 12 11 10 9 8 7 6 5 4 3 2 1 (Xp, Yp, Zp) Track 1 Track 2 基于轨道信息粗配准原理 (Xm, Ym, Zm), (Vxm,Vym, Vzm) (Xs, Ys, Zs), (Vxs,Vys, Vzs) R min R S P 基于轨道信息的粗配准 Gamma 执行操作语句： init_offset_orbit 例：输入数据slc.parslc.par 输出数据20090930.off 偏移量估计 Gamma 执行操作语句： offset_pwr 例： 输出数据20090930.offs,20090930.off.snr,20090930.offsets 20090930.offsets 偏移量估计 △R = a0+a1R2+a2C2+a3R+a4C+a5RC △C= b0+b1R2+b2C2+b3R+b4C+b5RC S m (Rm, Cm) (Rs, Cs) 偏移量多项式系数估计 Gamma 执行操作语句： offset_fit 例： 输出数据20090930.coffsets,20090930.coffs 20090930.coffsets 偏移量多项式系数估计 △R 1= a0+a1R12+a2C12+a3R1+a4C1+a5R1C1 △C1= b0+b1R12+b2C12+b3R1+b4C1+b5R1C1 △R 2= a0+a1R22+a2C22+a3R2+a4C2+a5R2C2 △C2= b0+b1R22+b2C22+b3R2+b4C2+b5R2C2 。。。 最小二乘法 a0 a1 a2 a3 a4 a5 b0 b1 b2 b3 b4 b5 干涉图和强度图的生成 Gamma 执行操作语句： interf_SLC 例： 输出数据20090930.int,20090930.p