基于相位信息另权的合成多重网格干涉SAR相位展开算法.pdfVIP

基于相位信息加权的合成多重网格干涉SAR相位展开算法 李笑郁毛士艺 (北京航空航天大学电子工程系203教研室，100083) 冀要：本文基于残差理论及HelⅢ_llol乜分解分析了干涉台成孔径雷达二维相位展开的理论 基础．进而论证了其最小均方意义下的解实际上是Poiss咖方程的解，使得相位展开的问题 转换成为Poisson方程的求解问题。为充分结合加权最小均方相位展开算法的优势，本文采 用了收敛速度快的多重网格方法，并在此基础上提出了一种基于相位信息加权的合成多重 网格算法。该算法结合了干涉sAR相位展开中常用的分支阻断方法简单快速并有一致解的 特点，使它为多重网格最小均方算法提供初始值，加速了多重阿格算法的收敛；并采用基 于相位导数方差的二值质量图替代相关系数作为权值输入，增强了算法的适应能力及可靠 性。实际x．sAR数据的相位展开处理结果验证了本文的分析和结论。 关链词：干涉sAR二维相位展开，相位导数方差二值质量图，合成加权多重网格算法 引言 干涉合成孔径雷达(IFsARl最主要的应用是利用同一地区的至少两幅在不同轨道位置或不 同时间获取的复数sAR图像形成干涉图像，重建地形拓扑图．今天它己被公认为是全天候 条件下获得大测绘带宽地球表面高度信息的强大手段【9l。由于系统的原因，干涉相位被限 制(缠绕)在(一石，石1之间。如何恢复模糊掉的相位周期，得到与数字高度图∞EM)相关的绝 对解缠相位．即相位展开。叶me unwrappin曲处理，成为干涉sAR处理的最关键步骤之～。 的相位不～致现象，即残差现象，使得对测量相位梯度进行简单的积分会得到与路径相关 的相位展开结果，导致误差的传递与积累。在两类最主要的干涉sAR相位展开方法中，积 分路径方法(Pa血．f01lo_Ivi唯Ⅲmod)建立在残差定理的基础上．试图实现与路径无关的1D 路径积分相位展开．基本思路是根据残差点或相位质量评价图设置分支或屏障于符号相反 的残差点之间，使积分路径绕过而不是穿过他们，从而得到与积分路径无关的一致解。其 代表算法GoldsIein的分支阻断(BrallchCul)法【‘J具有运算速度快，存储空间小的优点，缺点 是缺少基于先验知识的相位质量评价图对分支选择进行有效指导，从而导致展开错误及孤 立区域．另一类主要的相位展开方法是基于最小均方意义的．其思想是最小化测量相位梯 度与解平面梯度的均方差．许多相位展开的算法是在此基础上发展起来的，如mⅫero， Fr¨cesch甜i及Lanari的格林函数法m，Ghi画lia，Rornero咖的先验条件共轭梯度法俨cG)， 及PrInⅢ的加权多重网格算法。无加权最小均方方法可利用快速变换例．运算速度较快， 但不足之处在于相位展开时并未避开相位不一致的残差点，而是穿过这些区域．造成误差 在展开相位平面中的传播。GtIigiliaRomero【21等人引入了加权最小均方方法．以分离残差 该项目受到国家自然科学基金重点项目支持 点。加权最小均方算法的数值稳定性及避免误差传捶的性质，引起人们的极大关注·但加 权最小均方方法不能利用直接快速算法，只能利用Ga∞§-sc“ael迭代等经典迭代方法，而 此类迭代方法在相位存在不连续时收敛速度较慢．于是Prnt跚于1996年提出以多重网格算 法加速收敛过程的思想。该思想成功地解决了迭代收敛慢的缺点．获得了较理想的计算结 果．本文将加权最小均方相位展开算法具有数值稳定解的优势与分支阻断方法计算速度快 且通常有可靠解的优势相结合，使后者为加权多重网格最小均方方法提供初始值，从而得 到新的合成加权多重网格算法，新的合成算法不但继承了原来最小均方方法的稳定性及多 重网格的快速收敛的特性，而且由于利用分支截断方法提供初值，进一步提高了迭代的收 敛速度。另外，考虑到Prnt的多重网格方法未用任何先验信息．并且加权处理中的权系数 是根据复数sAR图像之间的相关系数得到的m，增加了计算量且对干涉SAR图象相关系数 的估计提出了要求，本文采用了基于相位数据的稳定的二值质量评价图取代相关系数权值， 增强了新算法的适应性和稳定性．实际典型山区x．SAR干涉数据的处理结果显示了本文提 出的算法的稳定性和有效性． 1．干涉相位残差分析与HeImhoI乜分解 在干涉处理过程中，我们得到的干涉相位矿与真实绝对相位≯之间的关系可由下式给 出： y=}y{≯)=mod{声+牙，2石)一石