发明名称 一种基于改进贝叶斯模型的立体图像生成方法
摘要 本发明属于图像处理领域,具体涉及一种基于改进贝叶斯模型的立体图像生成方法,可以更加高效、直观的生成立体图像。首先,由给定的视差计算出合成视角的逐像素视点位置;然后,通过改进的贝叶斯模型生成合成图像;最后,利用参考图像和合成图像生成红蓝立体图。与现有方法相比,本发明提出的基于改进贝叶斯模型的立体图像生成方法在时间性能、生成图像质量以及通用性等方便具有明显优势;该方法还可以直观逐像素指定视差,便于立体图像的后期处理,具有广阔的应用前景。
申请公布号 CN104902253A 申请公布日期 2015.09.09
申请号 CN201510064469.6 申请日期 2015.02.09
申请人 北京理工大学;北京理工大学深圳研究院 发明人 张磊;张驭航;黄华
分类号 H04N13/00(2006.01)I 主分类号 H04N13/00(2006.01)I
代理机构 代理人
主权项 一种基于改进贝叶斯模型的立体图像生成方法,其特征在于:其包括以下步骤:步骤一、通过参考图像得到合成图像的逐像素视点位置;所述参考图像是一幅指定的真实图像,参考图像所有像素视点位置相同;所述合成图像是通过一组输入图像得到的一幅合成图像,合成图像逐像素指定视点位置;步骤二、通过改进的贝叶斯模型生成合成图像;步骤2.1:针对合成图像中每个像素的视点位置s′(x),确定一个新成像平面Γ<sub>j</sub>;步骤2.2:给定场景的几何估计及一组输入图像v<sub>i</sub>,计算出每个输入图像的成像平面Ω<sub>i</sub>中每个像素点到新成像平面Γ<sub>j</sub>的映射关系φ<sub>ij</sub>以及输入图像的成像平面Ω<sub>i</sub>中每个像素点映射到新成像平面Γ<sub>j</sub>后是否可见的可见性算子m<sub>ij</sub>;所述输入图像v<sub>i</sub>中包含的像素数与参考图像v<sub>0</sub>中像素数相同;其中,i∈[1,n],n≥1;j∈[1,N];步骤2.3:通过公式(4)和公式(5)计算出输入图像成像平面Ω<sub>i</sub>中每个像素点到合成图像的成像平面Γ<sub>0</sub>的映射关系τ<sub>i</sub>及其逆映射关系β<sub>i</sub>;由于输入图像的像素数和合成图像的像素数均与参考图像的像素数相同,因此x也表示输入图像和合成图像中像素点的索引值;τ<sub>i</sub>(x)=φ<sub>ix</sub>(x)    (4)其中,x表示输入图像和合成图像中的像素点的索引值;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&tau;</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000113918010000011.GIF" wi="1861" he="58" /></maths>步骤2.4:根据公式(6)计算合成图像u的估计<img file="FSA0000113918010000012.GIF" wi="53" he="44" /><img file="FSA0000113918010000013.GIF" wi="1979" he="107" />其中,n为输入图像的数量;о为函数复合运算符;步骤2.5:根据公式(7)计算每幅输入图像v<sub>i</sub>相对于合成图像u的误差因子κ<sub>i</sub>;<img file="FSA0000113918010000014.GIF" wi="1862" he="138" />其中,<img file="FSA0000113918010000015.GIF" wi="40" he="53" />为梯度运算符;τ<sub>i</sub>是由τ<sub>i</sub>(x)构成的向量,x∈[1,N];z<sub>i</sub>为输入图像v<sub>i</sub>中像素点的深度值;将误差因子κ<sub>j</sub>由小到大进行排序,选出误差因子最小的M幅图像作为最终输入图像,并记录下图像编号l<sub>k</sub>,如公式(8)所示;l<sub>k</sub>=i    (8)其中,k∈[1,M],l<sub>k</sub>∈[1,n];步骤2.6:构造如公式(9)所示的能量项E(u);E(u)=E<sub>data</sub>(u)+λE<sub>prior</sub>(u)    (9)其中,E<sub>data</sub>(u)可通过公式(10)得到;其由贝叶斯图像生成模型推理得出;λ为调配系数,λ>0;E<sub>prior</sub>(u)是平滑项,用来填充合成图像u中缺少输入图像信息的部分,E<sub>prior</sub>(u)可通过公式(12)得到;<img file="FSA0000113918010000021.GIF" wi="1862" he="133" />其中,<img file="FSA0000113918010000022.GIF" wi="440" he="95" /><img file="FSA0000113918010000023.GIF" wi="51" he="57" />为高斯传感噪声方差,为常量;<img file="FSA0000113918010000024.GIF" wi="66" he="53" />为输入图像l<sub>k</sub>的图像误差,其值可通过公式(11)计算得到;<img file="FSA0000113918010000025.GIF" wi="54" he="46" />为可见性算子;<img file="FSA0000113918010000026.GIF" wi="1863" he="151" />其中,<img file="FSA0000113918010000027.GIF" wi="61" he="52" />为几何误差<img file="FSA0000113918010000028.GIF" wi="94" he="53" />为已知量;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>prior</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><msub><mo>&Integral;</mo><msub><mi>&Gamma;</mi><mn>0</mn></msub></msub><mo>|</mo><mo>&dtri;</mo><mi>u</mi><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000113918010000029.GIF" wi="1863" he="112" /></maths>步骤2.7:对公式(9)构造的能量项E(u)进行最小化,求解合成图像u;具体为:采用步骤2.4得到的合成图像的估计<img file="FSA00001139180100000210.GIF" wi="27" he="40" />作为迭代初始值,利用快速迭代收缩阈值FISTA算法进行迭代求解;步骤三、利用参考图像和合成图像生成满足图像视差D(x)的红蓝立体图。
地址 100081 北京市海淀区中关村南大街5号
您可能感兴趣的专利