Insar处理过程及原理详解

Insar处理过程及原理详解基于轨道信息的粗配准流程图计算偏移量多项式系数的精配准干涉图和强度图的生成从光盘读取SLC数据选取主影像和从影像………读取数据及选取主从影像原始光盘数据解压缩并获取所需头文件及影像文件例：SO_000012648_0002_1.tar.gz例：头文件--TSX1_SAR__SSC______SM_S_SRA_20090930T221744_20090930T221752.xml影像文件--IMAGE_HH_SRA_strip_012.cos选取成像时间居中的影像为主影像，其余为从影像例：2009.06.23的影像为主影像，其余为从影像，目的是为了降低影像配准时时间上可能产生的影响读取数据例：20090930.slc及20090930.slc.par影像数据参数说明29795*0.00026849613s=7.999842s≈8s80272-80264=8sdtdRdDR0影像数据参数说明41°Rnear=652144.4173mRcenterRfar=672850.1688m·(672850.1688-652144.4173)/sin41°≈31km41°29795*1.896414≈56km给定行列号（r,c），则可计算出该像元的斜距R与成像时间t:R=Rnear+c×dRt=t0+r×dt影像数据参数说明轨道数据轨道点个数影像数据参数说明121110987654321121110987654321(Xp,Yp,Zp)(Xm,Ym,Zm),(Vxm,Vym,Vzm)(Xs,Ys,Zs),(Vxs,Vys,Vzs)67Track1Track2基于轨道信息粗配准原理X=a0+a1t+a2t2+a3t3Y=b0+b1t+b2t2+b3t3Z=c0+c1t+c2t2+c3t3Vx=d0+d1t+d2t2+d3t3Vy=e0+e1t+e2t2+e3t3Vz=f0+f1t+f2t2+f3t3121110987654321121110987654321(Xp,Yp,Zp)(Xm,Ym,Zm),(Vxm,Vym,Vzm)(Xs,Ys,Zs),(Vxs,Vys,Vzs)Track1Track2h121110987654321121110987654321(Xp,Yp,Zp)Track1Track2基于轨道信息粗配准原理(Xm,Ym,Zm),(Vxm,Vym,Vzm)(Xs,Ys,Zs),(Vxs,Vys,Vzs)RminR222))(())(())((tZZtYYtXXspspsp1)()(:Ellipsoid22222hbZhaYXPPP||:sphererangeSlantSPRSPVfDop/2:coneDopplerSP基于轨道信息的粗配准Gamma执行操作语句：init_offset_orbit例：输入数据—20090623.slc.par及20090930.slc.par输出数据—20090623-20090930.off偏移量估计Gamma执行操作语句：offset_pwr例：输出数据—20090623-20090930.offs,20090623-20090930.off.snr,20090623-20090930.offsets20090623-20090930.offsets偏移量估计△R=a0+a1R2+a2C2+a3R+a4C+a5RC△C=b0+b1R2+b2C2+b3R+b4C+b5RCSm(Rm,Cm)(Rs,Cs)偏移量多项式系数估计Gamma执行操作语句：offset_fit例：输出数据—20090623-20090930.coffsets,20090623-20090930.coffs20090623-20090930.coffsets偏移量多项式系数估计△R1=a0+a1R12+a2C12+a3R1+a4C1+a5R1C1△C1=b0+b1R12+b2C12+b3R1+b4C1+b5R1C1△R2=a0+a1R22+a2C22+a3R2+a4C2+a5R2C2△C2=b0+b1R22+b2C22+b3R2+b4C2+b5R2C2。。。最小二乘法a0a1a2a3a4a5b0b1b2b3b4b5干涉图和强度图的生成Gamma执行操作语句：interf_SLC例：输出数据—20090623-20090930.int,20090623-20090930.pwr1,20090623-20090930.pwr2强度图