高精度配准与高精度鲁棒的相位解缠算法与软件——结题报告.ppt

2. 基于最大堆排序的质量图导引快速相位解缠 建立初始最大堆示意图 最大堆删除根结点示意图 在最大堆的最后插入新结点的二叉树 将插入新结点的二叉树调整为最大堆 质量指导的路径跟踪示意图 Sigma自适应滤波 Sigma自适应综合滤波 非线性ECS自配准成像算法 + 基于最大相关函数的分块多级配准处理进行配准； 对干涉复数据首先进行5×5快速均值滤波与7×7快速准中值滤波，进而对滤波后的实部和虚部进行质量图指导的Sigma自适应滤波获取高精度干涉相位； 对干涉相位利用基于最大堆的质量图导引法进行相位展开； 对展开后的干涉相位进行5×5快速均值滤波和7×7中值滤波； 利用初始基线参数和展开相位等进行干涉定标，获取干涉定标参数：水平基线长度 、垂直基线长度及基准相位； 利用定标后的参数进行DEM反演获取场景DEM； 对DEM进行5×5均值滤波和7×7中值滤波，获取最终DEM。 干涉处理算法： 大小：12000×7000，场景大小为：715m×938m。 具体位置：方位向1：12000，距离向：1001：8000。 数据块一： 定标参数为水平基线0.09023m，垂直基线0.29502m，基准相位12.55918rad。定标后反演的DEM精度为0.3114m，再对DEM进行5×5均值滤波及7×7中值滤波后DEM精度为0.2963m。若进行5×5均值和9×9中值滤波后DEM精度为0.2885m，由处理结果分析，不会损坏DEM的细节信息。 高程伪彩色图 灰度三维透视图 大小：16384×7000，场景大小为：976m×938m 具体位置：方位向1:16384，距离向：1:7000 处理数据块二： 定标水平基线0.09039m，垂直基线0.29544m，基准相位6.27465rad。定标后反演得到的DEM精度为0.3228m，再对DEM进行5×5均值滤波后高程精度为0.3070m，再进行9×9中值滤波后DEM精度为0.3024m。若进行7×7均值和9×9中值滤波后DEM精度为0.2896m。 灰度三维透视图 三基线高程重建算法 三基线高程重建算法 基于闭合形式鲁棒性中国余数定理的多通道干涉SAR 高程重建方法：也是先求解模糊数 ，然后再估计出同余方程组的解。还将该定理扩展到了实数范围（即 和 可以为实数），并且通过最优余数的选择和同余方程的差分操作使其鲁棒性和计算效率均有了更进一步的提高。 首先获取多个基线对应的滤波后的干涉相位； 利用参考平地获取去除平地后的干涉相位； 计算每个基线的高度模糊数： 求每个点的模糊数： 定义 可得同余方程组 令 ，则同余方程组化为 按照闭和形式鲁棒性中国余数定理的步骤 对目标像素的高程进行求解了 (9) 再对高程进行聚类分析, 去除噪声模糊类，进一步高程重建 仿真数据处理 高度变化剧烈的城市建筑物。最小高程为16米，最大高程为144米。场景既有非常平坦的区域，也有边沿十分陡峭的区域，相邻两点的最大高程差达128米。仿真了三个干涉通道的干涉条纹图，其对应的高度模糊数分别为21.4米、32.1米和53.5米，均小于最大高程差，由此导致相邻两点的最大干涉相位差分别达到11.96π、7.98π和4.79π。不符合相位连续性的假设，因此不能采用传统的基于梯度差的单基线相位解缠的方法实现高程重建。三幅干涉条纹图的相关系数均为0.95，添加的相位噪声仍然按照相干系数对干涉相位概率密度函数关系生成。 实测数据处理 天线12，13和23所形成干涉系统的模糊高度分别为57.469m 、41.0223m和143.3422m，为满足互素条件，选择为：58.8, 42, 142.8。 问题探讨： （1）多个基线最优化布置的问题； （2）闭合式鲁棒性CRT算法的进一步完善； （3）单基线与多基线联合高程重建算法的研究。 （4）多基线干涉SAR层析成像技术：通过在垂直于视线 方向依次增加多个基线，获得高度维的分辨能力， 使多基线SAR系统实现3D成像。 （差分层析成像技术具有基线与时间域进行联合数据 处理能力，可为未来高度-速度应用领域的需求提供 整体解决方案）。 （5）多基线重建相对高程精度的研究 对于叠掩区域，可以利用多基线InSAR 技术来检测分辨多个散射源并进行高程估计，可以考虑两种方法：一是利用信号谱估计的方法估计出叠掩区域多个散射源的干涉相位频率，进而恢复各自的绝对相位和高程（包括Beamforming，Capon 等非参量估计方法以及 Music，Relax 等参量估计方法）；二是利用统计方法对多基线InSAR 信号进行建模，利用最大似然估计、最大后验概率估计