一种多源多轨长条带卫星遥感影像联合并行匹配方法，引入特征评价算子，统计分析影像局部纹理的空间结构与反差；采用多轨道卫星影像联合匹配策略，保证了轨道间连接点的最大重叠次数；在全球SRTM数据的辅助下，利用地形起伏变化建立近似核线几何约束，极大地减少了由于高程误差造成的核线预测误差；综合运用了小面元几何纠正法、基于控制网的匹配生长算法，并通过不断精化差异补偿模型，消除了多源影像间的几何与辐射畸变，最终获得稳健地、高精度匹配结果；将匹配流程与匹配策略进行了整合与改进，为多源多轨长条带卫星影像的联合处理提供了较为完善的匹配方案，为海量卫星遥感数据的集成应用提供了保证。
【技术实现步骤摘要】
本专利技术属于测绘科学与
，涉及一种，主要应用于空中三角测量、数字表面模型生产等领域。
技术介绍
近20年来，随着航天技术的飞速发展，在轨卫星数量的急剧增加，越来越多的卫星数据应用到摄影测量领域，利用卫星数据进行测绘产品的生产已成为地理空间信息获取的重要手段。由于卫星成像系统在高空飞行时，具有覆盖范围广、运行周期短、姿态稳定、多余观测多等特点，这对于影像匹配中“病态解”的消除、匹配精度与可靠性的提高具有较大意义。另外，卫星上搭载的恒星定位仪与星敏器可获得固定采样间隔下扫描行的姿轨数据，从而通过几何约束方法来缩小匹配搜索范围，获得相对可靠的未知参数的初始值。然而卫星成像系统一般采用推扫式CCD线阵成像方式，具有不完全中心投影成像几何特点，其成像的几何特征远比传统的中心透视投影复杂，并且由于多视传感器(如前/正/后视)交会角、焦距、地面分辨率等差异，导致影像间存在复杂的几何与辐射畸变。这对于卫星影像的匹配研究带来了挑战。现有基于卫星遥感影像的匹配方法大致可分为三类:(I)基于物方的匹配:由于特殊的立体观测结构和线阵传感器上集成GPS接收器与星敏感器，基于物方的匹配算法应用较为广泛。但该算法依赖于精确姿态轨道参数获得的大致准确的同名预测点，考虑到国产卫星由于硬件系统的原因，其无控对地定位精度一般较差，单纯采用该算法很难取得较好的匹配结果。(2)基于不变性算子的匹配算法:由于基于不变性算子的匹配算法在计算机视觉领域取得了重大突破，凭借其强大的匹配性能在目标识别、全景图拼接、从运动恢复结构、遥感影像配准等领域得到成功应用。但该算法存在内存消耗多、运算速度慢的问题，对于大范围的卫星遥感影像匹配适用性不强，可靠性和有效性有待进一步检验。(3)基于特征与灰度属性的匹配:该算法结合基于特征的匹配和基于区域的匹配算法特点，首先利用特征匹配算子获得种子点，以种子点为基础，利用区域增长，进行最小二乘匹配。该算法适用于平坦地形下的影像匹配，但对居民地、山地、丘陵等复杂地形影像匹配效果较差。总的来说，现有匹配方法主要针对同一传感器在同轨或异轨获取的同源分景立体影像进行处理，而无法直接对不同传感器在不同轨道下获取的长条带卫星影像进行联合处理。如果利用分景后的影像进行匹配及区域网平差，则在实际生产中，景与景之间平差精度很难保持一致，从而出现影像拼接问题。如果直接利用长条带影像进行整体匹配，可对影像内部进行整体纠正，既有利于消除影像拼接问题，又可为稀少控制或无控地区定位精度的提高提供了保证。另外，现有方法对于影像间存在的空间分辨率、入射角、时相、成像机理的差异，无法完全消除，导致匹配算法的失效；而对于纹理贫乏区域，匹配点的加密一直以来都是亟待解决的难点问题
技术实现思路
本专利技术的目的是提供一种，以实现在多机多核情况下，全自动地联合不同传感器在不同轨道下获取的多源、多轨海量影像数据，进行快速并行匹配处理。本专利技术的技术方案为一种，包括以下步骤:步骤一，在每一轨道内，选取接近于垂直摄影的长条带影像作为基准影像，采用物理分块得到相互独立的基准影像块，并记录每个基准影像块与原始的基准影像的对应关系；物理分块时，块与块之间在列与行方向上均具有一定重叠；步骤二，利用特征提取算子，对各基准影像块提取特征点，引入特征评价算子统计特征点的二维直方图，根据二维直方图将位于纹理贫乏区域的特征点剔除，保留局部纹理丰富的特征点作为待匹配点；步骤三，确定各轨道内长条带影像间的重叠关系，建立轨道内匹配任务列表；轨道内匹配任务列表记录与每一轨道内的各基准影像块在同一轨道中有重叠关系的影像，作为相应搜索影像，每个基准影像块与相应所有搜索影像构成一个匹配任务；步骤四，根据轨道内匹配任务列表，对所有匹配任务分别通过基于地形起伏的近似核线约束匹配，获得轨道内的同名点匹配结果；步骤五，确定轨道间的长条带影像重叠关系，建立轨道间匹配模型任务列表；轨道间匹配模型任务列表记录与每一轨道内的各基准影像块在所有轨道中有重叠关系的影像，作为相应搜索影像，每个基准影像块与相应各搜索影像分别构成一个匹配任务；步骤六，根据轨道间匹配模型任务列表，对所有匹配任务分别通过基于地形起伏的近似核线约束匹配，获得轨道间的同名点匹配结果；步骤七，根据轨道内 的同名点匹配结果和轨道间的同名点匹配结果，合并所有轨道内与轨道间匹配点，进行整体的误匹配检测与剔除，最终输出匹配结果。而且，步骤四中对任一匹配任务进行基于地形起伏的近似核线约束匹配，包括以下子步骤，(I)对当前处理的基准影像块，基于全球地形高程数据计算高程范围，预测各搜索影像中与基准影像块重叠的范围并读取；以原始的基准影像块为底层建立多层金字塔影像，层数记为N ;令i=l ；(2)对各搜索影像分别执行以下子步骤，在对所有搜索影像处理完底层的原始影像后进入步骤(3)，(2.1)对第i层金字塔建立近似核线方程，得到近似核线；(2.2)利用差异补偿模型改正局部畸变；(2.3)进行基于近似核线约束的二维相关匹配，得到初始的匹配点；(2.4)从初始的匹配点中剔除误匹配的匹配点，处理完第i层金字塔；(2.5)判断当前是否处理完底层的原始影像，即判断是否i=N，否则令i=l+l，返回步骤(2.1);返回步骤(2.1)时更新约束条件，实现方式如下，处理完第i层金字塔后，构建数字表面模型，用于返回步骤(2.1)处理下层金字塔时建立近似核线方程，此处i的取值为1，2，…N-1 ;并且根据处理完第i层金字塔所得匹配点，重新计算得到处理下层金字塔时的差异补偿模型的系数，此处i的取值为1，2，…N-1 ；(3)对匹配点加密与精化；(4)对加密与精化后的匹配点剔除误匹配。而且，步骤六中对任一匹配任务进行基于地形起伏的近似核线约束匹配，包括以下子步骤，(I)对当前处理的基准影像块，基于全球地形高程数据计算高程范围预测搜索影像中与基准影像块重叠的范围并读取；以原始的基准影像块为底层建立多层金字塔影像，层数记为N ;令i=l ；(2)对第i层金字塔建立近似核线方程，得到近似核线；(3)利用差异补偿模型改正局部畸变；(4)进行基于近似核线约束的二维相关匹配，得到初始的匹配点；(5)从初始的匹配点中剔除误匹配的匹配点，处理完第i层金字塔；(6)判断当前是否处理完底层的原始影像，即判断是否i=N，是则进入下一步骤(7)，否则令i=l+l，返回步骤(2);返回步骤(2)时更新约束条件，实现方式如下，处理完第i层金字塔后，构建数字表面模型，用于返回步骤(2)处理下层金字塔时建立近似核线方程，此处i的取值为1，2，…N-1 ;并且根据处理完第i层金字塔所得匹配点，重新计算得到处理下层金字塔时的差异补偿模型的系数，此处i的取值为1，2，…N-1 ；(7)对匹配点加密与精化；`(8)对加密与精化后的匹配点剔除误匹配。而且，步骤二中利用LBP/C算子统计特征点的局部LBP图与LC图的像素值分布，得到LBP/C 二维直方图。本专利技术的优点在于直接利用长条带卫星影像，联合不同传感器在不同轨道下获取的多源、多轨海量数据，在多机多核的硬件条件下，完全自动化地进行快速并行匹配，避免了通过影像分景处理带来的景与景之间的影像拼接问题，实现了不同分辨率、不同视角、不同时相的卫星影像联合匹配，解决了纹理贫乏区域本文档来自技高网...
【技术保护点】
一种多源多轨长条带卫星遥感影像联合并行匹配方法，包括以下步骤：步骤一，在每一轨道内，选取接近于垂直摄影的长条带影像作为基准影像，采用物理分块得到相互独立的基准影像块，并记录每个基准影像块与原始的基准影像的对应关系；物理分块时，块与块之间在列与行方向上均具有一定重叠；步骤二，利用特征提取算子，对各基准影像块提取特征点，引入特征评价算子统计特征点的二维直方图，根据二维直方图将位于纹理贫乏区域的特征点剔除，保留局部纹理丰富的特征点作为待匹配点；步骤三，确定各轨道内长条带影像间的重叠关系，建立轨道内匹配任务列表；轨道内匹配任务列表记录与每一轨道内的各基准影像块在同一轨道中有重叠关系的影像，作为相应搜索影像，每个基准影像块与相应所有搜索影像构成一个匹配任务；步骤四，根据轨道内匹配任务列表，对所有匹配任务分别通过基于地形起伏的近似核线约束匹配，获得轨道内的同名点匹配结果；步骤五，确定轨道间的长条带影像重叠关系，建立轨道间匹配模型任务列表；轨道间匹配模型任务列表记录与每一轨道内的各基准影像块在所有轨道中有重叠关系的影像，作为相应搜索影像，每个基准影像块与相应各搜索影像分别构成一个匹配任务；步骤六，根据轨道间匹配模型任务列表，对所有匹配任务分别通过基于地形起伏的近似核线约束匹配，获得轨道间的同名点匹配结果；步骤七，根据轨道内的同名点匹配结果和轨道间的同名点匹配结果，合并所有轨道内与轨道间匹配点，进行整体的误匹配检测与剔除，最终输出匹配结果。...
【技术特征摘要】
【专利技术属性】
技术研发人员：张永军，熊金鑫，段延松，
申请(专利权)人：武汉大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。