2026《InSAR基本原理及误差来源分析》5600字.docxVIP

InSAR基本原理及误差来源分析

目录

TOC\o1-3\h\u708InSAR基本原理及误差来源分析 1

145561.1InSAR地表形变监测基本原理 1

133031.2InSAR复影像配准精度分析 2

186081.1.1InSAR系统模型 2

76391.1.2干涉相位误差统计特性 3

133451.1.3配准误差 5

2441.3干涉基线距误差分析 5

144871.3.1基线误差分析 5

40431.3.2基线误差对干涉相位的影响 6

180531.3.3基线校正 7

23741.4时间失相关 7

308261.5空间失相关 8

521.6大气延迟误差 9

140481.6.1大气延迟误差成因 9

326891.6.2大气延迟效应对干涉相位的影响 9

66401.6.3大气延迟改正 10

1.1InSAR地表形变监测基本原理

如图2-1所示为重复轨道干涉InSAR系统的几何示意图，S1和S2分别代表雷达两次经过相同区域成像时的天线的相位中心，S1和S2之间的距离是基线B，它们与地面目标P的距离分别为R1和R2，α是S1和S2之间的空间距离与水平方向的夹角，θ是S1处雷达侧视角，h为测量目标点P的高度，H是天线到参考椭球的垂直高度。在三角形S1S2P中，采用余弦定理可得：

R22=

其中R1=R2

δR=?B2

在雷达干涉测量系统中，可以认为基线B?R2，

δR=Bsinθ?α

相位差为：

δ?=?4πλ

高程为:

h=H?R1

图2-1InSAR几何示意图

Figure2-1InSARgeometrydiagram

理论上来说，InSAR干涉相位?由多个相位的分量组成，具体如下：

?=?ref

其中?ref为参考椭球相位，?top为地形相位，?def为形变相位，?atm为大气相位，?noi为噪声相位。参考椭球相位也就是平地相位。由平地相位引起，在雷达观测成像时因卫星传感器对观测点的入射角和斜距的不同而导致相位发生变化。加入精密的卫星轨道数据即可去除。地形相位由地形的变化引起，使用精密轨道数据和DEM数据可以去除部分的地形相位而得到残余地形相位。形变相位是由卫星两次观测地表时，地表发生了形变引起的相位变化，主要体现在沿雷达视线方向的变化量。而噪声相位以高频分量的形式表现在相位模型中，通过低通相位滤波来抑制。大气相位在后面小节有具体的描述。

1.2InSAR复影像配准精度分析

1.1.1InSAR系统模型

假设SAR获取到的由物体反射的回波信号为复平稳高斯过程，则天线获取的信号可表示为

P=fx,y

式中，f(x,y)为观测区域中散射单元；W(x,y)=sinc(πx/Rx)·sinc(πy/Ry)为SAR系统的脉冲响应函数，Rx为雷达系统方位向的分辨率，Ry为雷达系统距离向的分辨率；r为卫星到地面目标的距离；n为复高斯噪声。令雷达影像干涉对在方位向和距离向的配准误差为Δx和Δy,假设Δx=0，Δy≠0，且基线的水平分量Bx=B，垂直分量By=0，当观测目标由独立的均匀散射中心组成时，通过积分得到干扰回波信号的互相关函数为

P1?

式中，σ0为平均雷达散射截面，而相关系数γ可表示为

γ=γbaseline

式中，γbaseline为基线去相关；γreg为由配准误差导致的去相关。

1.1.2干涉相位误差统计特性

由式（2-8）和式（2-9）可知，复影像精度配准是直接影响相位精度的首个环节。因此一定要分析配准误差在干涉相位中造成的影响程度。干涉相位与相干系数的精度的相关性非常高。假设散射体为高斯分布，则干涉相位的概率密度函数为

pdf?=

式中，Ф0为真实相位。如果将相位限制在Ф0中的土π范围，有

E?

σ?2=

式中，Li2为二次对数函数，Li2γ2=

图2-2干涉系数估计随独立像素样本个数变化示意图

Figure2-2Schematicdiagramofinterferencecoefficientestimationwiththenumberofindependentpixelsamplesinthewindow

（a）干涉相位随相干系数的概率密度函数曲线

（a）Probabilitydensityfunctioncurveofinterferencephaseversuscoherencecoefficient

（b）干涉相位随相干系数和估计窗大小的关系

（b）Therelationshipbetweentheinterferencephaseandthecoher