一种光学卫星遥感影像数字表面模型提取方法及系统与流程

1.本发明属于遥感数据处理技术领域，尤其涉及一种光学卫星遥感影像数字表面模型提取方法及系统。


背景技术：

2.密集匹配(dense matching)在光学卫星遥感影像数字表面模型(digital surface model,dsm)提取的核心环节，匹配结果与匹配效率直接决定最终的作业成果质量和作业效率。密集匹配通常可分为代价计算、代价聚合、视差计算和视差优化四部分。按照代价聚合的策略不同，可分为局部方法和(半)全局方法。相较于局部方法，(半)全局方法计算复杂，可靠性更高。2008年经典的半全局匹配方法(semi-global matching,sgm)提出以来，其利用多路径一维代价聚合近似实现二维代价聚合，兼顾匹配精度与效率，在dsm提取等领域得到了广泛应用，目前也已被geomatica、inpho、erdas imaging、rpc stereo processor(rsp)、micmac及sure等业内领先的遥感影像处理软件所采用。
3.sgm中代价聚合根据邻域点匹配代价和平滑项更新当前点匹配代价，通过多个一维方向代价聚合，实现二维平滑约束。该策略需要对所有像素的所有候选视差值对应的匹配代价进行聚合处理，以单路径代价聚合为例，其计算量为n
×
d，n为单路径像素数目，d为视差范围，sgm中的代价聚合策略将单个像素的视差范围内的所有视差值同等对待，存在大量的计算冗余，特别是当视差范围较大时，严重影响匹配效率。


技术实现要素：

4.本发明解决的技术问题是：克服现有技术的不足，提供了一种光学卫星遥感影像数字表面模型提取方法及系统，解决了半全局密集匹配存在大量计算冗余，导致密集匹配耗时较长，dsm提取作业效率低的问题。
5.本发明目的通过以下技术方案予以实现：一种光学卫星遥感影像数字表面模型提取方法，包括：对每个像点采用census测度方法得到每个像点的匹配代价信息；根据每个像点的匹配代价信息基于最小代价路径的匹配代价聚合得到每个像点的候选视差值集合；从每个像点的候选视差值集合中选取最大支持度的候选视差值为每个像点的最终视差值，根据每个像点的最终视差值得到视差图；对视差图进行精化处理得到精化后的视差图；基于精化后的视差图和预设影像有理多项式参数，通过空间前方交会得到物方点云，对物方点云通过格网化处理得到光学卫星遥感影像数字表面模型。
6.上述光学卫星遥感影像数字表面模型提取方法中，所述census测度方法为：以两模板内对应的像点与各自中心像点间灰度值的大小关系的二值化信息构建特征空间，从而得到两个census序列，而后以census序列间的汉明距离为相似测度。
7.上述光学卫星遥感影像数字表面模型提取方法中，若某像点的灰度值大于等于中心像点灰度值，则该像点在特征空间的取值为1，否则为0。
8.上述光学卫星遥感影像数字表面模型提取方法中，根据每个像点的匹配代价信息
基于最小代价路径的匹配代价聚合得到每个像点的候选视差值集合包括：对于每个像点，通过确定每个像点的8邻域方向上的最小代价路径，获取每个像点的8领域方向上的候选视差值，形成每个像点的候选视差值集合。
9.上述光学卫星遥感影像数字表面模型提取方法中，最小代价路径通过如下公式得到：
[0010][0011]
其中，l为最小代价路径，i为路径上一点，n为路径上点数，j为相邻两像点之间的连线，ci为像点i选定视差值的匹配代价，cj为相邻两像点选定视差值时的平滑值。
[0012]
上述光学卫星遥感影像数字表面模型提取方法中，从每个像点的候选视差值集合中选取最大支持度的候选视差值为每个像点的最终视差值包括：先统计每个像点候选视差值集合中每个视差候选值的支持度，而后采用wta策略，选择支持度最大的候选视差值作为每个像点的最终视差值。
[0013]
上述光学卫星遥感影像数字表面模型提取方法中，每个视差候选值的支持度为：每个候选视差值与其他候选视差值的差值绝对值之和的倒数作为每个候选视差值的支持度。
[0014]
上述光学卫星遥感影像数字表面模型提取方法中，所述精化处理包括：中值滤波、左右一致性检测、孤点检测与内插。
[0015]
上述光学卫星遥感影像数字表面模型提取方法中，空间前方交会的数学模型通过如下公式得到：
[0016][0017][0018]
其中，(ln，sn)为像点的归一化像方坐标，(uv，w)为物方点的归一化物方坐标；
[0019]
num
l
(u，v，w)＝a1+a2v+a3u+a4w+a5vu+a6vw+a7uw+a
sv2
+a9u2+a
10
w2+a
11
vuw+a
12v3
+a
13
vu2+a
14
vw2+a
15v2
u+a
16
u3+a
17
uw2+a
18v2
w+a
19
u2w+a
20
w3；
[0020]
den
l
(u，v，w)＝b1+b2v+b3u+b4w+b5vu+b6vw+b7uw+b
8v2
+b9u2+b
10
w2+b
11
vuw+b
12v3
+b
13
vu2+b
14
vw2+b
15v2
u+b
16
u3+b
17
uw2+b
18v2
w+b
19
u2w+b
20
w3；
[0021]
nums(u，v，w)＝c1+c2v+c3u+c4w+c5vu+c6vw+c7uw+c
sv2
+c9u2+c
10
w2+c
11
vuw+c
12v3
+c
13
vu2+c
14
vw2+c
15v2
u+c
16
u3+c
17
uw2+c
18v2
w+c
19
u2w+c
20
w3；
[0022]
dens(u，v，w)＝d1+d2v+d3u+d4w+d5vu+d6vw+d7uw+d
8v2
+d9u2+d
10
w2+d
11
vuw+d
12v3
+d
13
vu2+d
14
vw2+d
15v2
u+d
16
u3+d
17
uw2+d
18v2
w+d
19
u2w+d
20
w3；
[0023]
其中，num
l
(u，v，w)、den
l
(u，v，w)、nums(u，v，w)和dens(u，v，w)均为物方点的归一化物方坐标的多项式，aq、bq、cq和dq均为影像有理多项式参数，其中，q＝1、2、
…
、20，q为下标。
[0024]
一种光学卫星遥感影像数字表面模型提取系统，包括：第一模块，用于对每个像点采用census测度方法得到每个像点的匹配代价信息；第二模块，用于根据每个像点的匹配代价信息基于最小代价路径的匹配代价聚合得到每个像点的候选视差值集合；第三模块，
用于每个像点的候选视差值集合中选取最大支持度的候选视差值为每个像点的最终视差值，根据每个像点的最终视差值得到视差图；第四模块，用于对视差图进行精化处理得到精化后的视差图；第五模块，用于基于精化后的视差图和预设影像rpc参数，通过空间前方交会得到物方点云，对物方点云通过格网化处理得到光学卫星遥感影像数字表面模型。
[0025]
本发明与现有技术相比具有如下有益效果：
[0026]
本发明使视差图计算复杂度不再与视差范围相关，大幅降低了计算冗余，提高了处理效率。
附图说明
[0027]
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
[0028]
图1是本发明实施例提供的光学卫星遥感影像数字表面模型提取方法的流程示意图；
[0029]
图2是本发明实施例提供的census匹配代价示意图；
[0030]
图3是本发明实施例提供的最小代价路径示意图；
[0031]
图4是本发明实施例提供的点云格网化示意图。
具体实施方式
[0032]
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
[0033]
图1是本发明实施例提供的光学卫星遥感影像数字表面模型提取方法的流程示意图。如图1所示，该方法包括如下步骤：
[0034]
步骤1：对每个像点采用census测度方法得到每个像点的匹配代价信息。
[0035]
采用census测度逐像点计算其在视差范围内每一个视差值对应的匹配代价，获取所有像点在所有视差范围内的匹配代价。
[0036]
步骤2：根据每个像点的匹配代价信息基于最小代价路径的匹配代价聚合得到每个像点的候选视差值集合。
[0037]
对于单个像素，通过分别确定其8邻域方向上的最小代价路径，获取其8领域方向上的候选视差值，形成该像素的候选视差值集合，逐像素处理完毕后得到所有像素的候选视差值集合。
[0038]
步骤3：从每个像点的候选视差值集合中选取最大支持度的候选视差值为每个像点的最终视差值，根据每个像点的最终视差值得到视差图。
[0039]
对于单个像素在代价聚合后得到的候选视差值集合，以单个视差值与其余视差值之间的差值的绝对值累加值，作为该视差值的支持度，之后采用“赢者通吃”(winner take all,wta)的策略，选择支持度最大的候选视差值作为该像素的最终视差值，逐像素处理完
毕后得到视差图。
[0040]
步骤4：对视差图进行精化处理得到精化后的视差图。
[0041]
采用中值滤波等策略，对视差图中的噪点进行精化处理，得到精化后的视差图。
[0042]
步骤5：基于精化后的视差图和预设影像有理多项式参数(rational polynomial coefficients，rpc)，通过空间前方交会得到物方点云，对物方点云通过格网化处理得到光学卫星遥感影像数字表面模型。
[0043]
具体的，匹配代价计算：
[0044]
匹配代价用于衡量两点在相应特征空间的相似性，采用census测度计算。如图2所示，census以两模板内对应点与各自中心点间灰度值的大小关系的二值化信息构建特征空间，从而得到两个census序列，而后以census序列间的汉明距离为相似测度。若某点的灰度值大于等于中心点灰度值，则该点在特征空间的取值为1，否则为0。
[0045]
基于最小代价路径的代价聚合：
[0046]
最小代价路径指单条聚合路径上点与点之间连线的代价值之和最小，如下式所示：
[0047][0048]
其中，中，l为最小代价路径，i为路径上一点，n为路径上点数，j为相邻两点之间的连线，ci为点i选定视差值的匹配代价，cj为相邻两点选定视差值时的平滑值。
[0049]
如图3所示，与sgm代价聚合策略相同，mcp单路径起始点p视差值采用其最小匹配代价对应的视差值d0。但与sgm对当前路径此后每个点的所有视差候选值对应的匹配代价均进行聚合不同，mcp仅从后点q的d0,d0+1,d
0-1及dc四个视差候选值中选择最小代价路径，其中dc为点q最小匹配代价对应的视差值。路径代价为q点对应视差值匹配代价及其与p点确定视差值d0之间路径的平滑值之和。以q点d0视差值为例，其中平滑项与sgm一致。依次向后重复处理，直至当前路径所有点处理完毕。可以看出在单路径代价聚合中，与sgm单个像点的所有视差候选值均进行处理不同，mcp单个像点仅有四个视差候选值参与处理，其计算量为n
×
4，n为单路径像素数目，将有效减少计算冗余，提高处理效率，特别是当视差范围d取值较大时。
[0050]
基于最大支持度的视差确定：
[0051]
与sgm代价聚合得到新的匹配代价不同，mcp得到n个视差候选值，n为路径数目，即每条路径提供一个关于像点最终视差值的候选值。此时不能直接采用wta确定视差，首先统计每个视差候选值的支持度，采用单个视差候选值与其他视差候选值的差值绝对值之和的倒数作为其支持度，如下式所示。而后采用wta策略，选择支持度最大的视差候选值作为最终视差值。
[0052]
d＝argmax 1/(∑|d
i-dj|)
i，j∈d
；
[0053]
其中，d视差候选值集合，d为最终视差值。
[0054]
视差精化：
[0055]
采用中值滤波、左右一致性检测、孤点检测与内插等策略对视差图进行精化。
[0056]
dsm生成：
[0057]
根据密集匹配得到的视差图和影像rpc参数，通过空间前方交会和格网化生成
dsm。空间前方交会是指在相对定向后，根据影像内外方位元素、同名点坐标，利用同名光线前交物方点坐标的过程。其数学模型为有理函数模型。
[0058]
num
l
(u，v，w)＝a1+a2v+a3u+a4w+a5vu+a6vw+a7uw+a
8v2
+a9u2+a
10
w2+a
11
vuw+a
12v3
+a
13
vu2+a
14
vw2+a
15v2
u+a
16
u3+a
17
uw2+a
18v2
w+a
19
u2w+a
20
w3[0059]
den
l
(u，v，w)＝b1+b2v+b3u+b4w+b5vu+b6vw+b7uw+b
8v2
+b9u2+b
10
w2+b
11
vuw+b
12v3
+b
13
vu2+b
14
vw2+b
15v2
u+b
16
u3+b
17
uw2+b
18v2
w+b
19
u2w+b
20
w3nums(u，v，w)＝c1+c2v+c3u+c4w+c5vu+c6vw+c7uw+c
8v2
+c9u2+c
10
w2+c
11
vuw+c
12v3
+c
13
vu2+c
14
vw2+c
15v2
u+c
16
u3+c
17
uw2+c
18v2
w+c
19
u2w+c
20
w3[0060]
dens(u，v，w)＝d1+d2v+d3u+d4w+d5vu+d6vw+d7uw+d
8v2
+d9u2+d
10
w2+d
11
vuw+d
12v3
+d
13
vu2+d
14
vw2+d
15v2
u+d
16
u3+d
17
uw2+d
18v2
w+d
19
u2w+d
20
w3[0061]
其中，aq、bq、cq和dq均为影像有理多项式参数，其中，q＝1、2、
…
、20，q为下标，一般情况下，b1,d1均取值为1。
[0062]
前方交会得到物方点云，但由于物方点云并非规则排列，而dsm本质上是一系列规则格网对物方空间的离散表达，每个格网内采用固定的高程值，因此格网化是由物方点云生成dsm的必要处理。如图4所示，在物方局部范围内，采用3*3的格网r进行离散表达，每个格网单元的高程值则由落入其中的物方点的高程值确定。常用的高程值确定策略包括不规则三角网内插法、中值法等。
[0063]
本实施例还提供了一种光学卫星遥感影像数字表面模型提取系统，包括：第一模块，用于对每个像点采用census测度方法得到每个像点的匹配代价信息；第二模块，用于根据每个像点的匹配代价信息基于最小代价路径的匹配代价聚合得到每个像点的候选视差值集合；第三模块，用于每个像点的候选视差值集合中选取最大支持度的候选视差值为每个像点的最终视差值，根据每个像点的最终视差值得到视差图；第四模块，用于对视差图进行精化处理得到精化后的视差图；第五模块，用于基于精化后的视差图和预设影像rpc参数，通过空间前方交会得到物方点云，对物方点云通过格网化处理得到光学卫星遥感影像数字表面模型。
[0064]
本发明首先通过基于最小代价路径的匹配代价聚合，将密集匹配中单个像素在其视差范围内的聚合代价转变为候选视差值集合；而后通过统计单个像素候选视差值的支持度，采用“赢者通吃”(winner take all,wta)策略，选择支持度最大的候选视差值作为该像素的最终视差值，逐像素处理完毕后得到视差图；使视差图计算复杂度不再与视差范围相关，大幅降低了计算冗余，提高了处理效率。
[0065]
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。