CN105589072A - 基于双相机阵列的合成孔径目标多视角成像方法 - Google Patents

基于双相机阵列的合成孔径目标多视角成像方法 Download PDF

Info

Publication number
CN105589072A
CN105589072A CN201610044284.3A CN201610044284A CN105589072A CN 105589072 A CN105589072 A CN 105589072A CN 201610044284 A CN201610044284 A CN 201610044284A CN 105589072 A CN105589072 A CN 105589072A
Authority
CN
China
Prior art keywords
camera
angle
formula
synthetic aperture
plane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610044284.3A
Other languages
English (en)
Other versions
CN105589072B (zh
Inventor
裴炤
张艳宁
孙莉
马苗
汪西莉
郭敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Wheatstone Intelligent Instrument Co.,Ltd.
Original Assignee
Shaanxi Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaanxi Normal University filed Critical Shaanxi Normal University
Priority to CN201610044284.3A priority Critical patent/CN105589072B/zh
Publication of CN105589072A publication Critical patent/CN105589072A/zh
Application granted granted Critical
Publication of CN105589072B publication Critical patent/CN105589072B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9058Bistatic or multistatic SAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9029SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Studio Devices (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于双相机阵列的合成孔径目标多视角成像方法,用于解决现有目标成像方法成像精度低的技术问题。技术方案是采用前后对称的双面棋盘标定板,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug布局对双相机阵列进行标定,依据各相机投影在聚焦平面的方差对各个深度合成孔径图像进行统计,最终获得目标所在最优聚焦平面的成像结果。在充分利用多视角信息的基础上,不依赖于目标状态信息,解决多个相机阵列联动时合成孔径目标成像问题,即在一定的范围内,用双相机阵列的互信息共享,实现目标多视角成像,提高了目标间深度差距较小时成像的精度。

Description

基于双相机阵列的合成孔径目标多视角成像方法
技术领域
本发明涉及一种目标成像方法,特别是涉及一种基于双相机阵列的合成孔径目标多视角成像方法。
背景技术
文献“Anovelmulti-objectdetectionmethodincomplexsceneusingsyntheticapertureimaging,PatternRecognition,2012,p1637–1658,April2012”公开了一种基于单相机阵列的动目标成像方法。该方法首先选定某一个相机视角作为参考视角,通过基于混合高斯模型的背景差对各个视角获得的图像数据进行前景提取,然后将获得的各个前景概率图像通过平面加视差的方法投影到不同深度,并加和求均值。采用一个给定的阈值对合成后的前景概率图进行统计即可获得最终目标所在的深度,从而较好的解决了相机阵列合成孔径成像中动目标成像问题。然而,该方法仅适用于动态变化的目标成像,且当目标间深度差距较小时,这种方法将会失效。
发明内容
为了克服现有目标成像方法成像精度低的不足,本发明提供一种基于双相机阵列的合成孔径目标多视角成像方法。该方法采用前后对称的双面棋盘标定板,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug布局对双相机阵列进行标定,依据各相机投影在聚焦平面的方差对各个深度合成孔径图像进行统计,最终获得目标所在最优聚焦平面的成像结果。在充分利用多视角信息的基础上,不依赖于目标状态信息,解决多个相机阵列联动时合成孔径目标成像问题,即在一定的范围内,用双相机阵列的互信息共享,实现目标多视角成像,可以提高目标间深度差距较小时成像的精度。
本发明解决其技术问题所采用的技术方案是:一种基于双相机阵列的合成孔径目标多视角成像方法,其特点是包括以下步骤:
步骤一、采用前后对称的双面棋盘标定版,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug配置对双相机阵列进行标定。平面πr为选定的参考平面,k为两相机阵列平面的交线。利用双面黑白格对称标定板,分别在距平面πr不同角度且均相交于k的若干平面上对双相机阵列进行标定,获得相机间距离参数ΔX,应用公式(1)得到视差Δp
Δ p = Δ X · θ - r θ , - - - ( 1 )
式中,θ为任意平面与某一相机阵列的夹角大小,r为参考平面与该相机阵列的夹角。
基于视差,将各相机阵列聚焦平面,投影至相交于直线k的平面簇上,应用公式(2)获得在穿过k任意倾斜角度平面的合成孔径成像结果
I i θ = I i r ( m + Δ p ) . - - - ( 2 )
式中,为相机阵列中各视角投影在参考平面πr的投影图像N为单相机阵列中的相机数目,i为相机序号。
步骤二、采用基于方差的合成孔径成像算法对各个深度合成孔径图像进行统计,应用公式(3)计算各视角各对应像素m的方差
var ( m ) = 1 N Σ i = 1 N ( I i θ ( m ) - u i θ ( m ) ) 2 , - - - ( 3 )
其中
u i θ ( m ) = 1 N Σ i = 1 N I i θ ( m ) . - - - ( 4 )
式中,为聚焦角度θ上对应m点的均值。设一个角度范围为D。
在统计的基础上,当焦平面接近最优聚焦角度时,目标上对应像素点m的方差将会越接近于H,H为一接近于0的值,应用公式(5)得到m在角度θ的标记值
f m &theta; = 1 var ( m ) < H 0 o t h e r w i s e , - - - ( 5 )
在当前存在目标区域范围内,搜索角度范围中最优角度是D中具有最小方差的角度,即
d=argminD{var(I)}.(6)
式中,d为某一相机阵列在探测区域D中的最优角度。
步骤三、从每个相机的参考视角中的位置确定合成孔径成像中的各像素值。
I i r ( m ) = I i r ( m - &Delta; p ) . - - - ( 7 )
通过被标记的合成孔径成像得到每个相机视角在不同深度的被标记的参考平面图像。再应用公式(8)将参考平面上每个相机视角的所有深度范围内的标记单独融合。
式中,&&表示在同一位置不同深度的每个标记的与操作。f表示某像素点的标记,若θ从0到π,中有一个标记为0,则为0,否则为1。将被标记的参考图像投影至第二相机阵列目标平面。计算每个相机视角的图像的平均值并去除被标记为0的像素。通过应用公式(9)去除遮挡后得到最终的透视结果。
式中,为最终的透视结果。
本发明的有益效果是:该方法采用前后对称的双面棋盘标定板,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug布局对双相机阵列进行标定,依据各相机投影在聚焦平面的方差对各个深度合成孔径图像进行统计,最终获得目标所在最优聚焦平面的成像结果。在充分利用多视角信息的基础上,不依赖于目标状态信息,解决多个相机阵列联动时合成孔径目标成像问题,即在一定的范围内,用双相机阵列的互信息共享,实现目标多视角成像。本发明采用双相机阵列布局模式,充分利用多相机视角信息,依据各相机投影在聚焦平面的方差对各个深度合成孔径图像进行统计,对目标最优深度进行估计;不依赖于目标状态信息,解决多个相机阵列联动时合成孔径成像中目标深度近似导致的成像失败问题,有效提高了目标间深度差距较小时多相机阵列视角成像的精度。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
本发明基于双相机阵列的合成孔径目标多视角成像方法具体步骤如下:
1、双相机阵列布局及参数优化。
首先对双相机阵列的布局及参数进行优化。假定从每一个相机中捕获图像序列Ci,i=1…N,两相机平面分别为πa和πb,将双相机阵列相机平面互成90度夹角布局,同时假设参考平面为πr其角度为r,任意焦平面为πθ,a其与相机阵列平面πa平面夹角为θ,在该布局中参考平面πr及任意焦平面πθ均相交于双相机平面交线k。在该布局中聚焦平面投影变换被简化为从πr的坐标系到指定的焦平面πθ,a的平移变换。利用双面黑白格对称标定板,分别在距平面πr不同角度且均相交于k的若干平面上对双相机阵列进行标定,获得相机间距离参数ΔX,然后应用公式(1)得到视差Δp。
&Delta; p = &Delta; X &CenterDot; &theta; - r &theta; , - - - ( 1 )
基于视差,可将各相机阵列聚焦平面,投影至相交于直线k的平面簇上,应用公式(2)
I i &theta; = I i r ( m + &Delta; p ) . - - - ( 2 )
获得在穿过k任意倾斜角度平面的合成孔径成像结果,式中,为相机阵列中各视角投影在参考平面πr的投影图像N为单相机阵列中的相机数目,i为相机序号。
2、双相机阵列合成孔径目标最优深度估计。
采用合成孔径成像可以透视遮挡物体,然而当目标与遮挡间深度接近时,传统的合成孔径成像算法将会失效。采用双相机阵列,同时对探测范围进行计算。应算公式(3)计算各视角各对应像素m的方差
var ( m ) = 1 N &Sigma; i = 1 N ( I i &theta; ( m ) - u i &theta; ( m ) ) 2 , - - - ( 3 )
其中
u i &theta; ( m ) = 1 N &Sigma; i = 1 N I i &theta; ( m ) . - - - ( 4 )
式中,为聚焦角度θ上对应m点的均值。
设一个角度范围为D,在统计的基础上,当焦平面接近最优聚焦角度时,目标上对应像素点m的方差将会越接近于H,H为一接近于0的值,应用公式(5)得到m在角度θ的标记值
f m &theta; = 1 var ( m ) < H 0 o t h e r w i s e , - - - ( 5 )
在当前存在目标区域范围内,搜索角度范围中最优角度是D中具有最小方差的角度,即
d=argminD{var(I)}.(6)
式中,d为某一相机阵列在探测区域D中的最优角度。
3、双相机阵列合成孔径成像。
通过被标记的合成孔径成像得到每个相机视角在不同深度的被标记的参考平面图像。应用公式(7),从每个相机视角的参考视角中的位置确定合成孔径成像中的各像素值。
I i r ( m ) = I i r ( m - &Delta; p ) . - - - ( 7 )
然后应用公式(8)将参考平面上每个相机视角的所有深度范围内的标记单独融合。
式中,&&表示在同一位置不同深度的每个标记的与操作。f表示某像素点的标记,若θ从0到π,中有一个标记为0,则为0,否则为1。将被标记的参考图像投影至第二相机阵列目标平面。计算每个相机视角的图像的平均值并去除被标记为0的像素。因此通过去除遮挡应用公式(9)可得最终的透视结果。
式中,为最终的透视结果。
因采用斜平面簇加视差法,且使用双相机阵列,解决了目标深度近似问题,因此优于单相机阵列探测方法。

Claims (1)

1.一种基于双相机阵列的合成孔径目标多视角成像方法,其特征在于包括以下步骤:
步骤一、采用前后对称的双面棋盘标定版,选定两相机阵列平面的交线作为所有聚焦平面簇的交线,利用Scheimpflug配置对双相机阵列进行标定;平面πr为选定的参考平面,k为两相机阵列平面的交线;利用双面黑白格对称标定板,分别在距平面πr,不同角度且均相交于k的若干平面上对双相机阵列进行标定,获得相机间距离参数ΔX,应用公式(1)得到视差Δp
&Delta; p = &Delta; X &CenterDot; &theta; - r &theta; , - - - ( 1 )
式中,θ为任意平面与某一相机阵列的夹角大小,r为参考平面与该相机阵列的夹角;
基于视差,将各相机阵列聚焦平面,投影至相交于直线k的平面簇上,应用公式(2)获得在穿过k任意倾斜角度平面的合成孔径成像结果
i i &theta; = I i r ( m + &Delta; p ) . - - - ( 2 )
式中,为相机阵列中各视角投影在参考平面πr的投影图像N为单相机阵列中的相机数目,i为相机序号;
步骤二、采用基于方差的合成孔径成像算法对各个深度合成孔径图像进行统计,应用公式(3)计算各视角各对应像素m的方差
var ( m ) = 1 N &Sigma; i = 1 N ( I i &theta; ( m ) - u i &theta; ( m ) ) 2 , - - - ( 3 )
其中
u i &theta; ( m ) = 1 N &Sigma; i = 1 N I i &theta; ( m ) . - - - ( 4 )
式中,为聚焦角度θ上对应m点的均值;设一个角度范围为D;
在统计的基础上,当焦平面接近最优聚焦角度时,目标上对应像素点m的方差将会越接近于H,H为一接近于0的值,应用公式(5)得到m在角度θ的标记值
f m &theta; = 1 var ( m ) < H 0 o t h e r w i s e , - - - ( 5 )
在当前存在目标区域范围内,搜索角度范围中最优角度是D中具有最小方差的角度,即
d=argminD{var(I)}.(6)
式中,d为某一相机阵列在探测区域D中的最优角度;
步骤三、从每个相机的参考视角中的位置确定合成孔径成像中的各像素值;
I i r ( m ) = I i r ( m - &Delta;p ) . - - - ( 7 )
通过被标记的合成孔径成像得到每个相机视角在不同深度的被标记的参考平面图像;再应用公式(8)将参考平面上每个相机视角的所有深度范围内的标记单独融合;
式中,&&表示在同一位置不同深度的每个标记的与操作;f表示某像素点的标记,若θ从0到π,中有一个标记为0,则为0,否则为1;将被标记的参考图像投影至第二相机阵列目标平面;计算每个相机视角的图像的平均值并去除被标记为0的像素;通过应用公式(9)去除遮挡后得到最终的透视结果;
式中,为最终的透视结果。
CN201610044284.3A 2016-01-22 2016-01-22 基于双相机阵列的合成孔径目标多视角成像方法 Active CN105589072B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610044284.3A CN105589072B (zh) 2016-01-22 2016-01-22 基于双相机阵列的合成孔径目标多视角成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610044284.3A CN105589072B (zh) 2016-01-22 2016-01-22 基于双相机阵列的合成孔径目标多视角成像方法

Publications (2)

Publication Number Publication Date
CN105589072A true CN105589072A (zh) 2016-05-18
CN105589072B CN105589072B (zh) 2017-11-17

Family

ID=55928798

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610044284.3A Active CN105589072B (zh) 2016-01-22 2016-01-22 基于双相机阵列的合成孔径目标多视角成像方法

Country Status (1)

Country Link
CN (1) CN105589072B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033166A (zh) * 2012-12-13 2013-04-10 南京航空航天大学 一种基于合成孔径聚焦图像的目标测距方法
GB2517659A (en) * 1993-06-23 2015-03-04 Dassault Electronique Radar for helicopters

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2517659A (en) * 1993-06-23 2015-03-04 Dassault Electronique Radar for helicopters
CN103033166A (zh) * 2012-12-13 2013-04-10 南京航空航天大学 一种基于合成孔径聚焦图像的目标测距方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SI-SENG WONG ET AL.: "Parasitic Calibration by Two-Step Ratio Approaching Technique for Split Capacitor Array SAR ADCs", 《IEEE》 *
ZHAO PEI ET AL.: "A novel multi-object detection method in complex scene using synthetic aperture imaging", 《PATTERN RECOGNITION》 *
裴炤 等: "多相机实时拼接视频影像技术", 《电子产品世界》 *

Also Published As

Publication number Publication date
CN105589072B (zh) 2017-11-17

Similar Documents

Publication Publication Date Title
US9451236B2 (en) Apparatus for synthesizing three-dimensional images to visualize surroundings of vehicle and method thereof
Saxena et al. Learning 3-d scene structure from a single still image
CN104616284A (zh) 彩色深度相机的彩色图像到深度图像的像素级对准算法
CN104050662B (zh) 一种用光场相机一次成像直接获取深度图的方法
CN103414872B (zh) 一种目标位置驱动ptz摄像机的方法
CN104156972A (zh) 基于激光扫描测距仪与多相机融合的透视成像方法
CN102609941A (zh) 基于ToF深度相机的三维注册方法
CN105258673B (zh) 一种基于双目合成孔径聚焦图像的目标测距方法、装置
CN105513063A (zh) Veronese映射和棋盘格标定拋物折反射摄像机
CN104598744A (zh) 一种基于光场的深度估计方法
CN107370994A (zh) 海域全景监控方法、装置、服务器及系统
CN109461185A (zh) 一种适用于复杂场景的机器人目标主动避障方法
CN104008542A (zh) 一种针对特定平面图形的快速角点匹配方法
CN104318576A (zh) 一种超像素级别的图像全局匹配方法
Ramisa et al. Combining invariant features and the ALV homing method for autonomous robot navigation based on panoramas
CN103593641A (zh) 基于立体摄像机的物体检测方法和装置
CN109961092A (zh) 一种基于视差锚点的双目视觉立体匹配方法及系统
Chan et al. An improved method for fisheye camera calibration and distortion correction
CN110148173B (zh) 车路协同中目标定位的方法和装置、电子设备、介质
CN104331890A (zh) 一种全局视差估计方法和系统
Román et al. Automatic Multiperspective Images.
CN103426170A (zh) 基于非结构光场合成孔径成像的隐藏目标成像方法
CN105589072A (zh) 基于双相机阵列的合成孔径目标多视角成像方法
CN103413304A (zh) 基于彩色深度融合的虚拟阵列合成孔径透视成像方法
Park et al. 48.2: Light field rendering of multi‐view contents for high density light field 3D display

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201224

Address after: No. 403, zone a, Shaanxi science and Technology Resources Coordination Center, No. 10, zhangbawu Road, high tech Zone, Xi'an City, Shaanxi Province, 710000

Patentee after: Xi'an Wheatstone Intelligent Instrument Co.,Ltd.

Address before: 710119 No. 620 West Chang'an Street, Chang'an District, Xi'an City, Shaanxi Province

Patentee before: Shaanxi Normal University

TR01 Transfer of patent right