一种多视角卫星影像立体像对优化选择方法技术

本发明专利技术属于摄影测量与遥感领域，公开了一种多视角卫星影像立体像对优化选择方法，包括基于重叠度的像对初选、基于交会角和特征匹配数量的像对优选、基于测区分块的像对精选，构建了完善的两级三类优化选择策略，两级包括整景影像和测区分块两种尺度，三类策略包括影像重叠度、交会角和特征匹配数量。该发明专利技术能够综合考虑匹配效果和交会精度的要求，在整景影像优化选择的基础上，考虑测区分块的实际覆盖情况，在更小的测区影像块范围内实现了立体像对的二次精选排序，得到针对每个测区分块对应的若干优化像对组合，在保证重建精度和完整度的同时，显著缩短了三维重建时间。显著缩短了三维重建时间。显著缩短了三维重建时间。

【技术实现步骤摘要】
一种多视角卫星影像立体像对优化选择方法


[0001]本专利技术属于摄影测量与遥感领域，公开了一种多视角卫星影像立体像对优化选择方法。

技术介绍

[0002]多视角光学遥感卫星影像三维重建是一种根据多张(至少两张)不同视角的光学卫星影像，利用高精度定位和密集匹配技术，生成三维点云并通过纹理映射重建三维模型表面的技术。对多视角卫星影像三维重建，通常采用如下成像模式：相邻三条轨道三视立体(通常称为三像对)附加沿着轨道方向前后拍摄的两景单片影像，共11景影像。这11景影像两两之间均可以构成立体，如果没有任何约束，则有55种立体像对的组合。多视角光学遥感卫星影像三维重建的时间、精度和完整度与立体像对的数量紧密相关。如果立体像对数量过多，可能重建的完整度较好，但是三维重建会花费大量的时间，并且质量较差的立体像对会影响三维重建精度的效果。如果像对数量过少，可以减少计算量，但无法充分发挥多视角卫星影像的信息优势，无法确保重建的完整性和精细程度，影响三维重建的效果。
[0003]现有的立体像对优化选择主要是针对多视角航空影像，主要基于传统交会角的方式来进行立体像对的选择。多视角航空相机通常附带5个不同角度的镜头：一个垂直下视、四个前后左右各侧摆45度的镜头，在立体像对构成中，通常采用同一镜头的两张影像构成立体，它们摄影角度基本一致，辐射和几何特性基本一致，重叠度、交会角基本固定不变。因此单纯利用交会角的优化选择就能满足实际的应用需求。对多视角卫星影像立体像对的构成，可以是同轨的两张影像，也可能是异轨的两张影像，也有可能是两种不同类型卫星影像，它们之间差别比较大，包括重叠度、辐射和几何差异以及交会角度等因素。这些因素都会对后续的三维重建造成影响。因此必须针对多视角卫星影像，综合考虑多种因素，进行立体像对的优化选择才能满足后续三维重建的需求。

技术实现思路

[0004]本专利技术的目的在于提供一种多视角卫星影像立体像对优化选择方法，旨在通过最优立体像对自动选择技术，能够从多视角卫星影像中选择适合匹配和质量较好的立体像对，有效兼顾处理效率和三维重建效果，在保证重建精度和重建完整度的同时，显著缩短三维重建时间，进而满足三维重建对影像数据的实际需求。
[0005]为了实现上述技术效果，本专利技术采用的技术方案是：
[0006]一种多视角卫星影像立体像对优化选择方法，包括如下步骤：
[0007]S1、输入若干景卫星影像，计算任意两张影像之间的重叠度，根据重叠度的阈值，筛选出一系列重叠度大于阈值的初始立体像对；
[0008]S2、对初始立体像对进行基于交会角和特征匹配数量的像对优选；具体包括如下步骤：
[0009]a1、选择初始立体像对的左像中心像点坐标q，通过左像RPC参数确定出通过q和左
像投影中心的一条物方直线l，并利用平均高程求出该直线上一物方点，将该点投影到立体像对的右像上的像点为q
′
；然后通过右像RPC参数确定出通过q
′
和右像投影中心的一条物方直线l
′
；最后计算两条共面直线的夹角作为两景影像的交会角；
[0010]a2、对初始像对的左像和右像采用SIFT特征匹配获取大量的同名点，并统计同名点数量作为特征匹配数量；
[0011]a3、根据交会角和特征匹配数量进行综合评价，将匹配点数量与交会角分值通过数学模型进行计算，得到该像对的最终分值；评价采用数学模型为：
[0012][0013][0014]式中，i,j表示影像ID；num
match
表示特征匹配点数量；angle表示交会角；s
angle
表示跟交会角相关的分值；σ1和σ2表示与交会角相关的平滑项因子，其中σ1表示小于10度时的平滑因子；σ2表示大于45度时的平滑因子；
[0015]a4、针对于每一个初始立体像对，都执行步骤a1
‑
a3的处理，以获得每一个初始立体像对评价分值，按照最终分值由大到小进行排序，取前M项立体像对作为筛选后立体像对；
[0016]S3、基于测区分块的像对精选
[0017]①
对整个测区按照物方空间进行分块，根据分块找出经过S2筛选后立体像对对应的影像块，每一个物方分块得到若干立体像对影像块；
[0018]②
对每一个物方分块区域，采用与S2中a1
‑
a3同样方法，对若干立体像对影像块进行基于交会角和特征匹配数量的评价；
[0019]③
根据评价结果对每一个物方分块找出若干优化的立体像对，选择前N项的立体像对作为每个测区分块对应的最优立体像对组合，其中N≤M。
[0020]进一步地，步骤S1中初始立体像对筛选流程具体包括：
[0021]1)对输入的若干景卫星影像的每一种立体像对组合，利用其定向参数将其左右影像四个角点按照测区平均高度投影到地面，获得对应的物方重叠范围；
[0022]2)根据物方重叠范围再投影到左右影像，计算重叠区域与原始影像的区域像素个数比例，得到左右影像各自的重叠度，选择左右影像中重叠度的较小值作为立体像对的重叠度指标；
[0023]3)根据设置阈值进行筛选，如果重叠度指标大于阈值，则将对应的立体像对列入候选初始像对，否则排除该像对；
[0024]4)完成所有立体像的重叠度筛选，得到筛选后的一系列立体像对。
[0025]与现有技术相比，本专利技术所具备的有益效果是：
[0026]1.将重叠度、交会角、匹配特征点数量作为评价度量引入卫星像对优化选择中，构建了完善的两级三类优化选择策略；两级包括整景影像和测区分块两种尺度，三类策略包括影像重叠度、交会角和特征匹配数量的评价策略。
[0027]2.能够综合考虑匹配效果和交会精度的要求，为每块区域选出足够数量的最优立体像对覆盖，在更小的测区影像块范围内实现了立体像对的二次精选排序，得到针对每个测区分块对应的若干优化像对组合，在保证重建精度和重建完整度的同时，显著缩短三维重建时间。
附图说明
[0028]图1为实施例中一种多视角卫星影像立体像对优化选择方法流程图。
具体实施方式
[0029]下面结合实施例及附图对本专利技术作进一步的详细描述。但不应将此理解为本专利技术上述主题的范围仅限于以下的实施例，凡基于本
技术实现思路
所实现的技术均属于本专利技术的范围。
[0030]实施例
[0031]参见图1，一种多视角卫星影像立体像对优化选择方法，包括如下步骤：
[0032]S1、输入若干景卫星影像，计算任意两张影像之间的重叠度，根据重叠度的阈值，筛选出一系列重叠度大于阈值的初始立体像对。
[0033]作为可实施方式的举例，本步骤通过以下处理过程实现：
[0034]①
对输入的若干景卫星影像的每一种立体像对组合，利用其定向参数将其左右影像四个角点按照测区平均高度投影到地面，获得对应的物方重叠范围；
[0035]②
根据物方重叠范围再投本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多视角卫星影像立体像对优化选择方法，其特征在于，包括如下步骤：S1、输入若干景卫星影像，计算任意两张影像之间的重叠度，根据重叠度的阈值，筛选出一系列重叠度大于阈值的初始立体像对；S2、对初始立体像对进行基于交会角和特征匹配数量的像对优选；具体包括如下步骤：a1、选择初始立体像对的左像中心像点坐标q，通过左像RPC参数确定出通过q和左像投影中心的一条物方直线l，并利用平均高程求出该直线上一物方点，将该点投影到立体像对的右像上的像点为q
′
；然后通过右像RPC参数确定出通过q
′
和右像投影中心的一条物方直线l
′
；最后计算两条共面直线的夹角作为两景影像的交会角；a2、对初始像对的左像和右像采用SIFT特征匹配获取大量的同名点，并统计同名点数量作为特征匹配数量；a3、根据交会角和特征匹配数量进行综合评价，将匹配点数量与交会角分值通过数学模型进行计算，得到该像对的最终分值；评价采用数学模型为：模型进行计算，得到该像对的最终分值；评价采用数学模型为：式中，i,j表示影像ID；num
match
表示特征匹配点数量；angle表示交会角；s
angle
表示跟交会角相关的分值；σ1和σ2表示与交会角相关的平滑项因子，其中σ1表示小于10度时的平滑因子；σ2表示大于45度时的平滑因子；a4、针对于每一个初始立体像对，...

【专利技术属性】
技术研发人员：巩丹超，刘松林，张威，单非非，
申请(专利权)人：中国人民解放军六一五四零部队，
类型：发明
国别省市：