发明名称 一种多视角拍摄进行孔径合成成像方法
摘要 本发明公开了一种多视角拍摄进行孔径合成成像方法,属于图像处理技术领域,包括以下几个步骤:步骤一:采用平面射影变换,对原始图像序列进行坐标变换,重采样合成对焦到非感兴趣平面的图像序列;步骤二:根据方差阈值法,对对焦到非感兴趣平面的像点进行识别,并在原始图像序列里去除非感兴趣目标平面∏的像点;步骤三:根据平面射影变换,对去除非感兴趣目标的原始图像序列进行坐标变换;步骤四:综合所有经变换后的图像坐标,重采样合成感兴趣目标平面,得的目标图像;本发明通过多次平面射影变换的方法,准确地对感兴趣目标平面进行透视成像,校正了平面相机阵列的安装误差,提高了感兴趣目标的对比度和清晰度。
申请公布号 CN102968784B 申请公布日期 2015.06.17
申请号 CN201210393904.6 申请日期 2012.10.17
申请人 北京航空航天大学 发明人 袁艳;颜玉慧;苏丽娟
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文利
主权项 一种多视角拍摄进行孔径合成成像方法,其特征在于,包括以下几个步骤:步骤一:采用平面射影变换,对原始图像序列进行坐标变换,重采样合成对焦到非感兴趣平面的图像序列;对不同视角的原始图像序列进行平面射影变换,将所有图像对焦到非感兴趣目标平面Π上,具体为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>&Pi;i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>&Pi;i</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><msub><mi>H</mi><mi>&Pi;i</mi></msub><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000664456740000011.GIF" wi="380" he="233" /></maths>其中:[x<sub>i</sub> y<sub>i</sub> 1]<sup>T</sup>是相机阵列中第i个相机所获得的包含感兴趣目标场景的原始图像齐次坐标,其坐标单位为像素,H<sub>Πi</sub>是原始图像对焦到非感兴趣平面Π的平面射影矩阵,[x<sub>Πi</sub> y<sub>Πi</sub> 1]<sup>T</sup>是原始图像坐标,[x<sub>i</sub> y<sub>i</sub> 1]<sup>T</sup>为经平面射影变换H<sub>Πi</sub>对应的齐次坐标;取第i张原始图像,将平面射影变换后的坐标[x<sub>Πi</sub> y<sub>Πi</sub> 1]<sup>T</sup>进行双线性重采样,得到原始图像对焦到非感兴趣目标平面Π的图像I<sub>Πi</sub>;同理,对其它原始图像序列进行上述操作,得到对焦到非感兴趣目标平面的图像序列;所述的步骤一中,H<sub>Πi</sub>由第i个相机的视角和对焦到非感兴趣目标平面的视角的相对位置共同决定,计算H<sub>Πi</sub>采用计算机视觉中的直接线性变换;步骤二:根据方差阈值法,对对焦到非感兴趣平面的像点进行识别,并在原始图像序列里去除非感兴趣目标平面Π的像点;采用方差阈值法判断非感兴趣目标平面Π的坐标在对焦到非感兴趣目标平面Π的图像序列I<sub>pi</sub>的坐标[x<sub>oc</sub> y<sub>oc</sub> 1]<sup>T</sup>;则非感兴趣目标平面Π的坐标在第i张原始图像的坐标[x<sub>oci</sub> y<sub>oci</sub> 1]<sup>T</sup>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>oci</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>oci</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><msubsup><mi>H</mi><mi>&Pi;i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>oc</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>oc</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000664456740000012.GIF" wi="404" he="236" /></maths>其中:在原图I<sub>i</sub>中将坐标[x<sub>oci</sub> y<sub>oci</sub> 1]<sup>T</sup>的邻域整像素全部剔除,即灰度设置为0,生成在原图基础上去除非感兴趣目标平面Π的图像I′<sub>i</sub>;当非感兴趣目标为覆盖了一定的深度范围时,判断所有深度非感兴趣平面是否全部去除,如果是,则进入步骤三,否则,返回步骤一,在合成目标图像之前,反复使用步骤一、二,逐次消除各深度平面内的非感兴趣目标;步骤三:根据平面射影变换,对去除非感兴趣目标的原始图像序列进行坐标变换;对去除非感兴趣目标的图像I′<sub>i</sub>进行平面射影变换,对焦到目标平面上;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>obji</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>obji</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><msub><mi>H</mi><mi>obji</mi></msub><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>x</mi><mi>i</mi><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000664456740000021.GIF" wi="413" he="245" /></maths>其中:H<sub>obji</sub>是图像I′<sub>i</sub>对焦到目标的平面射影矩阵,[x<sub>obji</sub> y<sub>obji</sub> 1]<sup>T</sup>是图像I′<sub>i</sub>坐标[x<sub>i</sub> y<sub>i</sub> 1]<sup>T</sup>经目标平面射影变换对应的齐次坐标;步骤四:综合所有经变换后的图像坐标,重采样合成感兴趣目标平面,得的目标图像;通过上述步骤得到N组去除非感兴趣目标且对焦到目标平面的坐标[x<sub>obj1</sub> y<sub>obj1</sub> 1]<sup>T</sup>,…[x<sub>obji</sub> y<sub>obji</sub> 1]<sup>T</sup>,…x<sub>objN</sub> y<sub>objN</sub> 1]<sup>T</sup>,对这N组坐标数据进行双线性重采样,合成对焦到目标平面的图像,最后,得到目标图像;当目标物体有一定的深度范围时,判断感兴趣目标物体各个深度平面的目标图像是否全部获取,如果是,得到目标物体在各深度处的合成孔径切片图像,对各切片图像进行图像融合,得到被遮拦目标物体的完整合成大孔径图像;否则,返回步骤三。
地址 100191 北京市海淀区学院路37号