发明名称 基于双相机阵列的合成孔径目标多视角成像方法
摘要 本发明公开了一种基于双相机阵列的合成孔径目标多视角成像方法,用于解决现有目标成像方法成像精度低的技术问题。技术方案是采用前后对称的双面棋盘标定板,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug布局对双相机阵列进行标定,依据各相机投影在聚焦平面的方差对各个深度合成孔径图像进行统计,最终获得目标所在最优聚焦平面的成像结果。在充分利用多视角信息的基础上,不依赖于目标状态信息,解决多个相机阵列联动时合成孔径目标成像问题,即在一定的范围内,用双相机阵列的互信息共享,实现目标多视角成像,提高了目标间深度差距较小时成像的精度。
申请公布号 CN105589072A 申请公布日期 2016.05.18
申请号 CN201610044284.3 申请日期 2016.01.22
申请人 陕西师范大学 发明人 裴炤;张艳宁;孙莉;马苗;汪西莉;郭敏
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于双相机阵列的合成孔径目标多视角成像方法,其特征在于包括以下步骤:步骤一、采用前后对称的双面棋盘标定版,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug配置对双相机阵列进行标定;平面π<sub>r</sub>为选定的参考平面,k为两相机阵列平面的交线;利用双面黑白格对称标定板,分别在距平面π<sub>r</sub>,不同角度且均相交于k的若干平面上对双相机阵列进行标定,获得相机间距离参数ΔX,应用公式(1)得到视差Δp<maths num="0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>p</mi><mo>=</mo><mi>&Delta;</mi><mi>X</mi><mo>&CenterDot;</mo><mfrac><mrow><mi>&theta;</mi><mo>-</mo><mi>r</mi></mrow><mi>&theta;</mi></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000011.GIF" wi="1221" he="113" /></maths>式中,θ为任意平面与某一相机阵列的夹角大小,r为参考平面与该相机阵列的夹角;基于视差,将各相机阵列聚焦平面,投影至相交于直线k的平面簇上,应用公式(2)获得在穿过k任意倾斜角度平面<img file="FDA00009123604200000110.GIF" wi="51" he="79" />的合成孔径成像结果<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>i</mi><mi>i</mi><mi>&theta;</mi></msubsup><mo>=</mo><msubsup><mi>I</mi><mi>i</mi><mi>r</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>&Delta;</mi><mi>p</mi><mo>)</mo></mrow><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000012.GIF" wi="1246" he="85" /></maths>式中,<img file="FDA0000912360420000013.GIF" wi="145" he="78" />为相机阵列中各视角投影在参考平面π<sub>r</sub>的投影图像<img file="FDA0000912360420000014.GIF" wi="567" he="79" />N为单相机阵列中的相机数目,i为相机序号;步骤二、采用基于方差的合成孔径成像算法对各个深度合成孔径图像进行统计,应用公式(3)计算各视角各对应像素m的方差<maths num="0003"><math><![CDATA[<mrow><mi>var</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msup><mrow><mo>(</mo><msubsup><mi>I</mi><mi>i</mi><mi>&theta;</mi></msubsup><mo>(</mo><mi>m</mi><mo>)</mo><mo>-</mo><msubsup><mi>u</mi><mi>i</mi><mi>&theta;</mi></msubsup><mo>(</mo><mi>m</mi><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000015.GIF" wi="1472" he="115" /></maths>其中<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>i</mi><mi>&theta;</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>I</mi><mi>i</mi><mi>&theta;</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000016.GIF" wi="1312" he="113" /></maths>式中,<img file="FDA0000912360420000017.GIF" wi="159" he="85" />为聚焦角度θ上对应m点的均值;设一个角度范围为D;在统计的基础上,当焦平面接近最优聚焦角度时,目标上对应像素点m的方差将会越接近于H,H为一接近于0的值,应用公式(5)得到m在角度θ的标记值<img file="FDA0000912360420000018.GIF" wi="71" he="87" /><maths num="0005"><math><![CDATA[<mrow><msubsup><mi>f</mi><mi>m</mi><mi>&theta;</mi></msubsup><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>var</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>&lt;</mo><mi>H</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></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><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000019.GIF" wi="1333" he="135" /></maths>在当前存在目标区域范围内,搜索角度范围中最优角度是D中具有最小方差的角度,即d=arg min<sub>D</sub>{var(I)}.  (6)式中,d为某一相机阵列在探测区域D中的最优角度;步骤三、从每个相机的参考视角中的位置确定合成孔径成像中的各像素值;<maths num="0006"><math><![CDATA[<mrow><msubsup><mi>I</mi><mi>i</mi><mi>r</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>I</mi><mi>i</mi><mi>r</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mi>&Delta;p</mi><mo>)</mo></mrow><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000912360420000021.GIF" wi="1285" he="79" /></maths>通过被标记的合成孔径成像得到每个相机视角在不同深度的被标记的参考平面图像;再应用公式(8)将参考平面上每个相机视角的所有深度范围内的标记单独融合;<img file="FDA0000912360420000022.GIF" wi="1256" he="90" />式中,&amp;&amp;表示在同一位置不同深度的每个标记的与操作;f表示某像素点的标记,若θ从0到π,<img file="FDA0000912360420000023.GIF" wi="149" he="87" />中有一个标记为0,则<img file="FDA0000912360420000024.GIF" wi="135" he="86" />为0,否则为1;将被标记的参考图像投影至第二相机阵列目标平面;计算每个相机视角的图像的平均值并去除被标记为0的像素;通过应用公式(9)去除遮挡后得到最终的透视结果;<img file="FDA0000912360420000025.GIF" wi="1604" he="213" />式中,<img file="FDA0000912360420000026.GIF" wi="160" he="77" />为最终的透视结果。
地址 710119 陕西省西安市长安区西长安街620号