发明名称 一种基于颜色内相关和自适应支撑权重的立体匹配方法
摘要 本发明提出一种基于颜色内相关和自适应支撑权重的立体匹配方法。本发明方法提出了像素点颜色内相关函数,并结合颜色相似性、欧式距离相似性和颜色内相关相似性来确定匹配窗口内像素点的权重大小,从而获得立体匹配的初始视差值,然后对初始视差图进行三步优化,分别为左右一致性检查、中值滤波、直方图统计。本发明平均准确率高,而且具有很强的鲁棒性。
申请公布号 CN103985128B 申请公布日期 2017.03.15
申请号 CN201410223156.6 申请日期 2014.05.23
申请人 南京理工大学 发明人 龚文彪;任建乐;陆恺立;刘琳;顾国华;钱惟贤;路东明;任侃;于雪莲;吕芳
分类号 G06T7/90(2017.01)I;G06T7/593(2017.01)I 主分类号 G06T7/90(2017.01)I
代理机构 南京理工大学专利中心 32203 代理人 唐代盛;孟睿
主权项 一种基于颜色内相关和自适应支撑权重的立体匹配方法,其特征在于,步骤如下:步骤一、定义左图像和右图像的像素点颜色内相关函数,获得颜色内相关分量rin,左图像与右图像的颜色内相关分量rin表达式相同,均如公式(1)所示,rin=[rce1,rce2,rce3]      (1)公式(1)中,rce1、rce2、rce3依次为左图像或者右图像的三个颜色内相关函数,具体如公式(2)、(3)、(4)所示,rce1=I<sub>R</sub>(x,y)‑I<sub>G</sub>(x,y)      (2)rce2=I<sub>G</sub>(x,y)‑I<sub>B</sub>(x,y)      (3)rce3=I<sub>B</sub>(x,y)‑I<sub>R</sub>(x,y)      (4)公式(2)、(3)、(4)中,I<sub>R</sub>(x,y)、I<sub>G</sub>(x,y)、I<sub>B</sub>(x,y)依次为左图像或者右图像像素点的R通道分量、G通道分量、B通道分量的像素值,x、y表示当前像素点在行方向和列方向上的坐标值;步骤二、建立如公式(5)所示的匹配窗口权重函数,<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>w</mi><mrow><mo>(</mo><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><msub><mi>&Delta;c</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>c</mi></msub></mfrac></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><msub><mi>&Delta;d</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>d</mi></msub></mfrac></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><msub><mi>&Delta;r</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>r</mi></msub></mfrac></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>&Delta;c</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>c</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>&Delta;d</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>d</mi></msub></mfrac><mo>+</mo><mfrac><mrow><msub><mi>&Delta;r</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub></mrow><msub><mi>&tau;</mi><mi>r</mi></msub></mfrac></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194419380000011.GIF" wi="1454" he="295" /></maths>公式(5)中,w(p,q)为匹配窗口权重函数,exp为指数函数符号,p为左图像匹配窗口的中心像素点,q为左图像匹配窗口中除中心像素点之外的像素点,Δc<sub>pq</sub>为像素点p和q的RGB颜色相似性,Δd<sub>pq</sub>为像素点p和q的欧式距离相似性,Δr<sub>pq</sub>为像素点p和q的颜色内相关分量rin的相似性,τ<sub>c</sub>为颜色相似性所占匹配窗口权重函数的比例,τ<sub>d</sub>为欧式距离相似性所占匹配窗口权重函数的比例,τ<sub>r</sub>为颜色内相关相似性所占匹配窗口权重函数的比例;公式(5)中,颜色相似性Δc<sub>pq</sub>、欧式距离相似性Δd<sub>pq</sub>、颜色内相关相似性Δr<sub>pq</sub>的计算方法依次如公式(6)、(7)、(8)所示:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;c</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mo>|</mo><mo>|</mo><mi>p</mi><mo>-</mo><mi>q</mi><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><mrow><msub><mi>p</mi><mi>R</mi></msub><mo>-</mo><msub><mi>q</mi><mi>R</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>p</mi><mi>G</mi></msub><mo>-</mo><msub><mi>q</mi><mi>G</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><msub><mi>p</mi><mi>B</mi></msub><mo>-</mo><msub><mi>q</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194419380000012.GIF" wi="1470" he="167" /></maths><maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;d</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mo>|</mo><mo>|</mo><msub><mi>d</mi><mi>p</mi></msub><mo>-</mo><msub><mi>d</mi><mi>q</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>-</mo><msub><mi>x</mi><mi>q</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>p</mi></msub><mo>-</mo><msub><mi>y</mi><mi>q</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194419380000021.GIF" wi="1486" he="191" /></maths><maths num="0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;r</mi><mrow><mi>p</mi><mi>q</mi></mrow></msub><mo>=</mo><mo>|</mo><mo>|</mo><msub><mi>rin</mi><mi>p</mi></msub><mo>-</mo><msub><mi>rin</mi><mi>q</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>1</mn><mi>p</mi></msub><mo>-</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>1</mn><mi>q</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>2</mn><mi>p</mi></msub><mo>-</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>2</mn><mi>q</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>3</mn><mi>p</mi></msub><mo>-</mo><mi>r</mi><mi>c</mi><mi>e</mi><msub><mn>3</mn><mi>q</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001194419380000022.GIF" wi="1483" he="191" /></maths>公式(6)、(7)、(8)中,p<sub>R</sub>、p<sub>G</sub>、p<sub>B</sub>依次为像素点p点的R、G、B颜色像素值,q<sub>R</sub>、q<sub>G</sub>、q<sub>B</sub>依次为像素点q点的R、G、B颜色像素值,(x<sub>p</sub>,y<sub>p</sub>)为像素点p的坐标值,(x<sub>q</sub>,y<sub>q</sub>)表示像素点q的坐标值,rin<sub>p</sub>为像素点p处的颜色内相关分量值,rin<sub>q</sub>为像素点q处的颜色内相关分量值;步骤三、根据匹配窗口权重函数w(p,q)计算左右图像立体匹配的初始视差值;步骤四、对初始视差值依次进行左右一致性检查、中值滤波和直方图统计,获得最终精确的视差结果。
地址 210094 江苏省南京市孝陵卫200号