Insar处理过程及原理详解

1、Insar处理过程及原理详解1基于轨道信息的粗配准流程图计算偏移量多项式系数的精配准干涉图和强度图的生成从光盘读取SLC数据选取主影像和从影像2读取数据及选取主从影像原始光盘数据解压缩并获取所需头文件及影像文件例：SO_000012648_0002_1.tar.gz例：头文件- TSX1_SAR_SSC_SM_S_SRA_20090930T221744_20090930T221752.xml影像文件- IMAGE_HH_SRA_strip_012.cos选取成像时间居中的影像为主影像，其余为从影像例：2009.06.23的影像为主影像，其余为从影像，目的是为了降低影像配准时时间上可能产生的影响

2、读取数据例：20090930.slc及20090930.slc.par3影像数据参数说明29795*0.s=7.999842s 8s80272-80264 = 8sdtdRdDR04影像数据参数说明41Rnear =652144.4173mRcenterRfar =672850.1688m(672850.1688-652144.4173)/sin41 31km4129795* 1.896414 56km给定行列号（r, c），则可计算出该像元的斜距R与成像时间t:R=Rnear+cdR t=t0+rdt5影像数据参数说明轨道数据轨道点个数6影像数据参数说明1211109876543211211

3、10987654321(Xp, Yp, Zp)(Xm, Ym, Zm), (Vxm,Vym, Vzm)(Xs, Ys, Zs), (Vxs,Vys, Vzs)67Track 1Track 27基于轨道信息粗配准原理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)

4、, (Vxs,Vys, Vzs)Track 1Track 28h121110987654321121110987654321(Xp, Yp, Zp)Track 1Track 2基于轨道信息粗配准原理(Xm, Ym, Zm), (Vxm,Vym, Vzm)(Xs, Ys, Zs), (Vxs,Vys, Vzs)RminRSP9基于轨道信息的粗配准Gamma 执行操作语句： init_offset_orbit例：输入数据20090623.slc.par及20090930.slc.par输出数据20090.off10偏移量估计Gamma 执行操作语句： offset_pwr例：输出数据20090.o

5、ffs, 20090.off.snr, 20090.offsets 20090.offsets11偏移量估计R = a0+a1R2+a2C2+a3R+a4C+a5RCC= b0+b1R2+b2C2+b3R+b4C+b5RCSm(Rm, Cm)(Rs, Cs)12偏移量多项式系数估计Gamma 执行操作语句： offset_fit例：输出数据20090.coffsets, 20090.coffs 20090.coffsets13偏移量多项式系数估计R 1= a0+a1R12+a2C12+a3R1+a4C1+a5R1C1C1= b0+b1R12+b2C12+b3R1+b4C1+b5R1C1R 2= a0+a1R22+a2C22+a3R2+a4C2+a5R2C2C2= b0+b1R22+b2C22+b3R2+b4C2+b5R2C2。最小二乘法a0a1a2a3a4a5b0b1b2b3b4b