发明名称 UWB InSAR干涉相位真值计算方法
摘要 本发明提供一种适用于低频UWB InSAR的干涉相位真值计算方法。技术方案是:首先对InSAR数据(两幅SAR图像)进行配准,从而生成干涉相位观测值(也称原始干涉相位)。再利用配准结果生成一个参考相位,将该参考相位从原始干涉相位中减去,得到的相位称之为误差相位。接着,对误差相位进行噪声抑制与相位解缠,并估计残余相位模糊,得到误差相位真值。最后将误差相位真值与参考相位相加得到干涉相位真值。本发明能够大大降低低频UWBInSAR干涉相位真值计算的复杂度,且普适性强,适用于具有很大的相对带宽的情况。
申请公布号 CN104459634A 申请公布日期 2015.03.25
申请号 CN201410742321.9 申请日期 2014.12.08
申请人 中国人民解放军国防科学技术大学 发明人 安道祥;许军毅;黄晓涛;黎向阳;李悦丽;周智敏
分类号 G01S7/02(2006.01)I;G01S13/90(2006.01)I 主分类号 G01S7/02(2006.01)I
代理机构 国防科技大学专利服务中心 43202 代理人 王文惠
主权项 一种UWB SAR干涉相位真值计算方法,其中UWB是指超宽带,SAR是指合成孔径雷达,其特征在于,包括下述步骤:第一步,图像配准与原始干涉相位生成:对两幅SAR图像,任选一幅SAR图像为主图像,另外一幅为辅图像进行配准,设经过配准后,主图像中的任意一点(x,y)对应于辅图像中的点(x′,y′);利用配准结果将辅图像重采样得到配准图像;将主图像与配准图像进行共轭相乘,得到原始干涉相位<img file="FDA0000627360610000011.GIF" wi="210" he="83" />利用主图像与配准图像计算相关系数图γ(x,y);第二步,计算误差相位:利用下式计算参考相位<img file="FDA0000627360610000012.GIF" wi="210" he="85" /><img file="FDA0000627360610000013.GIF" wi="793" he="148" />其中<img file="FDA0000627360610000014.GIF" wi="176" he="83" />指(x,y)点的参考相位,R<sub>m</sub>(x,y)指主图像中的点(x,y)的斜距,R<sub>s</sub>(x′,y′)指辅图像中点(x′,y′)的斜距,λ是SAR信号的波长;利用下式计算(x,y)点的误差相位<img file="FDA0000627360610000015.GIF" wi="211" he="84" /><img file="FDA0000627360610000016.GIF" wi="695" he="92" />其中W(·)指对相位进行缠绕;第三步,对误差相位进行噪声抑制与解缠:首先对误差相位进行噪声抑制;然后,利用最小不连续相位解缠方法对误差相位进行解缠,得到任意点(x,y)的解缠相位<img file="FDA0000627360610000017.GIF" wi="248" he="87" />第四步:估计解缠相位的残余相位模糊:利用下式估计残余相位模糊<img file="FDA0000627360610000018.GIF" wi="97" he="61" /><img file="FDA0000627360610000019.GIF" wi="723" he="112" />其中,<img file="FDA00006273606100000110.GIF" wi="224" he="101" />指对所有相关系数大于0.5的点(x,y)计算<img file="FDA00006273606100000111.GIF" wi="470" he="91" />的众数值;利用下式计算误差相位真值:<img file="FDA0000627360610000021.GIF" wi="566" he="84" />第五步,计算干涉相位真值利用下式计算(x,y)点的干涉相位真值<img file="FDA0000627360610000022.GIF" wi="194" he="83" /><img file="FDA0000627360610000023.GIF" wi="1220" he="168" />
地址 410073 湖南省长沙市开福区德雅路109号