主权项 |
一种基于双相机阵列的合成孔径目标多视角成像方法,其特征在于包括以下步骤:步骤一、采用前后对称的双面棋盘标定版,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug配置对双相机阵列进行标定;平面π<sub>r</sub>为选定的参考平面,k为两相机阵列平面的交线;利用双面黑白格对称标定板,分别在距平面π<sub>r</sub>,不同角度且均相交于k的若干平面上对双相机阵列进行标定,获得相机间距离参数ΔX,应用公式(1)得到视差Δp<maths num="0001"><math><![CDATA[<mrow><mi>Δ</mi><mi>p</mi><mo>=</mo><mi>Δ</mi><mi>X</mi><mo>·</mo><mfrac><mrow><mi>θ</mi><mo>-</mo><mi>r</mi></mrow><mi>θ</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>θ</mi></msubsup><mo>=</mo><msubsup><mi>I</mi><mi>i</mi><mi>r</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>Δ</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>Σ</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>θ</mi></msubsup><mo>(</mo><mi>m</mi><mo>)</mo><mo>-</mo><msubsup><mi>u</mi><mi>i</mi><mi>θ</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>θ</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msubsup><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>I</mi><mi>i</mi><mi>θ</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>θ</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><</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>Δ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" />式中,&&表示在同一位置不同深度的每个标记的与操作;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" />为最终的透视结果。 |