主权项 |
利用非线性矢量曲面构建InSAR相位图像模型的方法,其特征在于,包括以下步骤:步骤一:将相位干涉图以复数的浮点形式读取到Matlab中,标记矩阵行、列大小为M、N,单位为像素;将干涉相位图划分大小为k×l的初始估计窗口,并计算值得信任的第一个估计窗口,将窗口内相位点(m,n)单频信号形式转换为三角函数弦分量的复数相位形式:φ(m,n)=angle(exp[j(f<sub>x</sub>m+f<sub>y</sub>n)])=cos(f<sub>x</sub>m+f<sub>y</sub>n)+jsin(f<sub>x</sub>m+f<sub>y</sub>n),其中,f<sub>x</sub>和f<sub>y</sub>分别为相位点行向量和列向量的频率;步骤二:在估计窗口内,利用离散小波变换对干涉图像正弦分量和余弦分量进行多层分解,将相位表示为低频部分和高频部分和的形式;步骤三:在估计窗口内,利用相位点高频分量的最大值确定主矢量的方向,也即是下一步搜索窗口的方向;步骤四:在估计窗口内,计算相位主矢量,利用高频部分主矢量平滑窗口内的相位频率;步骤五:在估计窗口内,根据各层分解分量的高频矢量和构建顾及相位细节变化的非线性相位曲面模型,将相位表示成非线性相位和线性相位和的形式Φ=φ<sub>linear</sub>+φ<sub>non‑linear</sub>,对非线性部分利用<img file="FDA0001001442270000011.GIF" wi="963" he="135" />对j层分解的三个高频方向求矢量和,其中j=1,2…J,<img file="FDA0001001442270000012.GIF" wi="301" he="71" />分别表示对j层高频的水平、垂直和视线方向分量求矢量和,矢量求和保证了非线性曲面的特性;步骤六:对小波分解相位进行逆变换得到相位估计模型,在估计窗口k×l内,对相位点ψ<sub>k,l</sub>(m,n),利用<img file="FDA0001001442270000013.GIF" wi="1102" he="83" />可得到估计窗口内非线性相位估计模型,其中,WT<sup>‑1</sup>(·)为小波逆变换,φ<sub>Lp</sub>为小波分解低频部分,<img file="FDA0001001442270000014.GIF" wi="109" he="63" />为小波分解高频部分和。 |