发明名称 一种基于改进的自适应加权和置信度传播的视差估计方法
摘要 一种基于改进的自适应加权和置信度传播的视差估计方法,属于计算机立体视觉技术领域。本发明为了解决现有的匹配算法计算的视差在遮挡区域、深度不连续区域匹配精度低;高精度匹配算法的计算量大,可实施性差;基于分割的细化方法视差平面评估不准确的问题。本发明的技术要点为:利用加权等级变换方法计算匹配像素间的相关值;通过改进的自适应加权方法对遮挡像素进行重新匹配;利用改进的置信度传播算法对视差图进行全局优化;利用细化模块对视差图进行精细化;再次利用改进的置信度传播算法和相关值和视差图进行全局优化。本发明可应用于立体图像分割、立体视频编码、机器人视觉,目标跟踪等领域。
申请公布号 CN104166987B 申请公布日期 2017.01.25
申请号 CN201410323405.9 申请日期 2014.07.08
申请人 哈尔滨师范大学 发明人 荣宪伟;薛远洋;于晓艳
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于改进的自适应加权和置信度传播的视差估计方法,其特征在于所述方法是按照以下步骤实现的:步骤一、利用加权等级变换方法计算匹配像素间的相关值C<sub>L</sub>;步骤二、利用左右一致性检测方法检测图像中的遮挡像素,并通过改进的自适应加权方法对遮挡像素进行重新匹配,生成初始视差图D<sub>1</sub>和初始相关值C<sub>1</sub>;其中,自适应加权方法的改进过程如下:假设f(x,y)代表参考图像中某一点,f(x+i,y+j)代表以像素f(x,y)为中心的匹配窗口内的像素,窗口内像素权重的计算如公式(1)所示:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>W</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>&Delta;C</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><msub><mi>&beta;</mi><mrow><mi>a</mi><mi>w</mi></mrow></msub></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>&Delta;G</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><msub><mi>&gamma;</mi><mrow><mi>a</mi><mi>w</mi></mrow></msub></mfrac></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001113111350000011.GIF" wi="1174" he="143" /></maths>其中,f代表参考图像的像素值,β<sub>aw</sub>和γ<sub>aw</sub>分别为颜色相似度和空间邻近度的计算系数,ΔC代表中心像素f(x,y)与周围像素f(x+i,y+j)在Lab颜色空间的欧几里得距离,通过公式(2)计算,其中c代表Lab颜色通道,ΔG代表f(x,y)与f(x+i,y+j)间的空间距离,通过公式(3)计算:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>&Delta;C</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><munder><mo>&Sigma;</mo><mrow><mi>c</mi><mo>&Element;</mo><mo>{</mo><mi>L</mi><mi>a</mi><mi>b</mi><mo>}</mo></mrow></munder><msup><mrow><mo>(</mo><msub><mi>f</mi><mi>c</mi></msub><mo>(</mo><mi>x</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>j</mi><mo>)</mo><mo>-</mo><msub><mi>f</mi><mi>c</mi></msub><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001113111350000012.GIF" wi="1134" he="127" /></maths><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&Delta;G</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msup><mi>i</mi><mn>2</mn></msup><mo>+</mo><msup><mi>j</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001113111350000013.GIF" wi="978" he="86" /></maths>计算目标窗口的像素权重;步骤三、利用改进的置信度传播算法对视差图D<sub>1</sub>进行全局优化,生成新的视差图D<sub>2</sub>;其中,置信度传播算法的改进过程如下:在第t次迭代时刻,节点p传递到周围节点q的信息<img file="FDA0001113111350000014.GIF" wi="180" he="68" />的计算如公式(10)所示:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>m</mi><mrow><mi>p</mi><mo>&RightArrow;</mo><mi>q</mi></mrow><mi>t</mi></msubsup><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>d</mi></munder><msub><mi>&phi;</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><msub><mi>&psi;</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>I</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><munder><mo>&Pi;</mo><mrow><mi>s</mi><mo>&Element;</mo><mi>N</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>/</mo><mi>q</mi></mrow></munder><msubsup><mi>m</mi><mrow><mi>s</mi><mo>&RightArrow;</mo><mi>p</mi></mrow><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001113111350000015.GIF" wi="1267" he="111" /></maths>其中,N(p)/q代表除q以外的p的周围节点,φ<sub>p</sub>(d)代表节点q的局部匹配代价的聚合,利用相关值C<sub>1</sub>放大α<sub>bp</sub>倍后计算e的指数得到,如公式(11)所示:φ<sub>p</sub>(d)=exp(α<sub>bp</sub>×C<sub>p</sub>(d))                 (11)其中,α<sub>bp</sub>为恒定的算法系数,通过经验获得;Ψ<sub>p,q</sub>(I,d)代表平滑函数,当周围临近的像素违背平滑约束的时候做出惩罚,改进后如公式(12):ψ<sub>p,q</sub>(I,d)=exp(‑h<sub>1</sub>(△I<sub>pq</sub>)×h<sub>2</sub>(△d<sub>pq</sub>))                 (12)其中I代表像素的亮度值,d代表像素的视差,p和q代表两个相邻的像素,当周围临近的像素违背平滑约束的时候做出惩罚,其中函数h<sub>1</sub>和h<sub>2</sub>分别通过公式(13)、(14)表示:<maths num="0005"><math><![CDATA[<mrow><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>&Delta;I</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&mu;</mi><mrow><mi>b</mi><mi>p</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>if&Delta;I</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>&gt;</mo><msub><mi>&beta;</mi><mrow><mi>b</mi><mi>p</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&lambda;</mi><mrow><mi>b</mi><mi>p</mi></mrow></msub><mo>&times;</mo><msub><mi>&mu;</mi><mrow><mi>b</mi><mi>p</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001113111350000021.GIF" wi="1270" he="151" /></maths>h<sub>2</sub>(△d<sub>pq</sub>)=min(γ<sub>bp</sub>,|△d<sub>pq</sub>|)+δ<sub>bp</sub>                  (14)公式中ΔI<sub>pq</sub>代表图像中p和q的像素亮度差,β<sub>bp</sub>代表像素亮度差的阈值,λ<sub>bp</sub>和μ<sub>bp</sub>代表平滑惩罚系数,通过经验获得,Δd<sub>pq</sub>代表p和q视差的差,γ<sub>bp</sub>和δ<sub>bp</sub>为恒定的系数,通过经验获得,γ<sub>bp</sub>代表截断阈值,用于限制平滑强度的,而δ<sub>bp</sub>的设置用于保证一定量的平滑强度;步骤四、利用细化模块对视差图D<sub>2</sub>进行精细化,并生成相关值C<sub>2</sub>,视差图D<sub>3</sub>;步骤五、再次利用改进的置信度传播算法和相关值C<sub>2</sub>和视差图D<sub>3</sub>进行全局优化,得到最终的视差图D<sub>f</sub>。
地址 150080 黑龙江省哈尔滨市南岗区和兴路50号