Insar處理過程及原理詳解

1、Insar處理過程及原理詳解1基于軌道信息的粗配準(zhǔn)流程圖計算偏移量多項(xiàng)式系數(shù)的精配準(zhǔn)干涉圖和強(qiáng)度圖的生成從光盤讀取SLC數(shù)據(jù)選取主影像和從影像2讀取數(shù)據(jù)及選取主從影像原始光盤數(shù)據(jù)解壓縮并獲取所需頭文件及影像文件例：SO_000012648_0002_1.tar.gz例：頭文件- TSX1_SAR_SSC_SM_S_SRA_20090930T221744_20090930T221752.xml影像文件- IMAGE_HH_SRA_strip_012.cos選取成像時間居中的影像為主影像，其余為從影像例：2009.06.23的影像為主影像，其余為從影像，目的是為了降低影像配準(zhǔn)時時間上可能產(chǎn)生的影響

2、讀取數(shù)據(jù)例：20090930.slc及20090930.slc.par3影像數(shù)據(jù)參數(shù)說明29795*0.s=7.999842s 8s80272-80264 = 8sdtdRdDR04影像數(shù)據(jù)參數(shù)說明41Rnear =652144.4173mRcenterRfar =672850.1688m(672850.1688-652144.4173)/sin41 31km4129795* 1.896414 56km給定行列號（r, c），則可計算出該像元的斜距R與成像時間t:R=Rnear+cdR t=t0+rdt5影像數(shù)據(jù)參數(shù)說明軌道數(shù)據(jù)軌道點(diǎn)個數(shù)6影像數(shù)據(jù)參數(shù)說明1211109876543211211

3、10987654321(Xp, Yp, Zp)(Xm, Ym, Zm), (Vxm,Vym, Vzm)(Xs, Ys, Zs), (Vxs,Vys, Vzs)67Track 1Track 27基于軌道信息粗配準(zhǔn)原理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基于軌道信息粗配準(zhǔn)原理(Xm, Ym, Zm), (Vxm,Vym, Vzm)(Xs, Ys, Zs), (Vxs,Vys, Vzs)RminRSP9基于軌道信息的粗配準(zhǔn)Gamma 執(zhí)行操作語句： init_offset_orbit例：輸入數(shù)據(jù)20090623.slc.par及20090930.slc.par輸出數(shù)據(jù)20090.off10偏移量估計Gamma 執(zhí)行操作語句： offset_pwr例：輸出數(shù)據(jù)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偏移量多項(xiàng)式系數(shù)估計Gamma 執(zhí)行操作語句： offset_fit例：輸出數(shù)據(jù)20090.coffsets, 20090.coffs 20090.coffsets13偏移量多項(xiàng)式系數(shù)估計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