CN102663768B - 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 - Google Patents
偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 Download PDFInfo
- Publication number
- CN102663768B CN102663768B CN201210144015.6A CN201210144015A CN102663768B CN 102663768 B CN102663768 B CN 102663768B CN 201210144015 A CN201210144015 A CN 201210144015A CN 102663768 B CN102663768 B CN 102663768B
- Authority
- CN
- China
- Prior art keywords
- image
- polarization
- polarization image
- thoroughly
- fourier transform
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法,通过任意三幅已知透振方向的偏振图像,合成任意透振方向的偏振图像,通过对这些偏振图像进行清晰度的评价,从而为偏振视觉系统确定最佳的透振方向。本发明提供了一种廉价、客观、方便的方法,解决了现有偏振视觉系统中偏振镜方向难以确定的问题;而且确定了最佳透振方向后,可以最大程度消除杂散光影响,大大地提高成像质量。
Description
技术领域
本发明涉及偏振视觉系统领域,具体涉及一种偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法。
背景技术
偏振视觉系统的应用越来越广泛,如偏振识别、偏振导航,这类偏振视觉系统通过多幅已知透振方向的偏振图像,求取偏振光的stokes矢量或E矢量,作为特征进行目标识别或导航。这类视觉系统对偏振镜的透振方向没有特别的要求。偏振视觉另外一个重要的用途是减弱或消除杂散光对成像的影响。
很多被检测物体的表面会发出杂乱的眩光或环境中存在强的逆光,如透过水面或玻璃拍摄物体、对金属表面进行视觉检测、车载视觉或交通监控系统对面的车打强光灯、雨天的视觉监控系统等等,这对视觉检测会产生非常严重的负面影响。为了减弱或者消除杂散光、眩光、逆光等干扰,目前很多应用场合采用在镜头前面搭配偏振镜来解决问题。
为了达到最好的使用效果,需要旋转偏振镜片,选择最好的偏振透振方向。而对于不同的场合透振方向角是不同的,即便对同一场合透振方向角也是随着目标或光源的不同而实时发生改变的。而对于动态场景,是不允许花过多的时间调整偏振镜片的透振方向。
为了削弱杂散光的影响,现有的偏振视觉系统是通过如下2种方式确定偏振镜片的透振方向:(1)旋转偏振镜片,观察成像效果,从而确定最好的透振方向。(2)根据先验知识,确定偏振镜的透振方向,如司机佩戴的偏振太阳镜。
发明内容
为了克服现有偏振视觉系统无法快速客观地确定偏振镜透振方向的缺陷,本发明提供了一种能够提高成像质量,使得目前的偏振系统使用更加方便,应用更加广泛的偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法。
本发明解决其技术问题所采用的技术方案:
偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法,其特征在于,包括以下步骤:
1)利用任意三幅偏振镜的已知透振方向的偏振图像,合成任意透振方向的偏振图像;
当偏振图像的透振方向与参考方向(选定的0°方向)夹角为θ时,任意透振方向的偏振图像可表示为:
I、Q、U为偏振图像的Stokes矢量,I(θ)代表偏振图像在θ透振方向时出射光的光强,I表示总光强,Q表示X轴方向与Y轴方向偏振图像分量之差,U表示±450方向偏振图像分量之差;
(1)利用偏振装置采集某场景三幅不同透振方向的偏振图像,假设三幅已知透振方向的偏振图像的透振方向与参考方向的夹角分别为θ1、θ2、θ3,相应的偏振图像分别为I(θ1)、I(θ2)、I(θ3);
(2)对偏振图像进行配准
由于,同一场景多幅偏振图像的采集不是同一时刻的,场景中有运动目标或采集设备的震动会使得视场会轻微变化,因此采集的三幅图像并不严格一致,如果三幅图像不对准,将会给后面的计算stokes矢量带来误差,因此,需要先对偏振图像进行配准(如果成像系统非常稳定,此步骤可以省略);
假设I(θ1)为参考图像,I(θ2)和I(θ3)为待配准图像,相对于I(θ1),假设I(θ2)被平移(x0,y0),被旋转θ0角度;
①平移配准
参考图像I(θ1)记为I1(x,y),待配准图像I(θ2)记为I2(x,y)即:
I2(x,y)=I1(x-x0,y-y0) (2)
则2幅图像对应的傅里叶变换有如下关系:
由傅里叶变换的平移性质知,信号在空间域的平移对应于其在频域的相位移,这两幅图像的相位差就等于它们互功率谱的相位,归一化互功率谱计算如下:
其中,FI1是I1的傅里叶变换,FI2是I2傅里叶变换的复共轭,通过对互功率谱进行傅里叶逆变换得到一个脉冲函数δ(x-x0,y-y0),该函数的极值点位置就是对应的平移量(x0,y0);
②旋转配准
I1(x,y)存在旋转,则要先计算出旋转角度并且得到补偿之后才能进行平移计算;
假设I2(x,y)是I1(x,y)旋转θ0后得到的图像,根据傅里叶变换的旋转性质,旋转一幅图像,在频域相当于对其傅里叶变换作相同角度的旋转,M1、M2为FI1、FI2功率谱,满足如下关系:
M2(x,y)=M1(xcosθ0+ysinθ0,-xsinθ0+ycosθ0) (5)
为了突出图像中的高频信息,对傅里叶变换的功率谱M1、M2取对数运算,通过双线性插值将其转换到极坐标下:
M2(ρ,θ)=M1(ρ,θ-θ0) (6)
这样把直角坐标下求旋转量转换为极坐标下求平移量,采用平移配准介绍的方法即可求得θ0;
(3)只要求解出偏振图像的Stokes矢量I、Q、U,将求解Stokes矢量带入式(1),就可求得任意透振方向的偏振图像I(θ);
对偏振图像的Stokes矢量I、Q、U的主要求解过程为:
采用分水岭算法对偏振图像进行过分割,假设在每个分割块内的所有点的Stokes矢量是相同的,利用最小二乘法,使误差平方和最小,从而求解Stokes矢量;
2)采用图像清晰度评价函数对合成的偏振图像进行评价;
用图像对比度来衡量图像是否清晰,统计各个透振方向的偏振图像的对比度,采用统计图像水平、垂直和对角方向梯度表示图像对比度,图像对比度公式如下:
3)确定合成的偏振图像的最佳透振方向;
根据图像清晰度评价函数,即公式(7)评价偏振图像的清晰度,选择杂散光滤除最多的,目标最清晰的偏振图像,被选中的最清晰的图像为:
该偏振图像对应的透振方向即为最佳透振方向。
本发明的有益效果在于:
本发明仅需要任意三幅已知透振方向的偏振图像,就可编码合成各个透振方向的偏振图像,对图像清晰度进行客观评价,从而确定偏振镜的最佳透振方向。本发明并不需要实际采集各个偏振方向的图像,大大减少了时间的消耗。因此,本发明提供了一种廉价、客观、方便的方法,解决了现有偏振视觉系统中偏振镜方向难以确定的问题;而且确定了最佳透振方向后,可以最大程度消除杂散光影响,大大地提高成像质量。
附图说明
图1为本发明的平移配准算法流程图。
图2为本发明的旋转配准算法流程图。
图3为本发明的搜索策略流程图。
图4为本发明的整体流程图。
具体实施方式
如图4所示,偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法,包括以下步骤:
1)利用任意三幅偏振镜的已知透振方向的偏振图像,合成任意透振方向的偏振图像;
当偏振图像的透振方向与参考方向(选定的0°方向)夹角为θ时,任意透振方向的偏振图像可表示为:
I、Q、U为偏振图像的Stokes矢量,I(θ)代表偏振图像在θ透振方向时出射光的光强,I表示总光强,Q表示X轴方向与Y轴方向偏振图像分量之差,U表示±450方向偏振图像分量之差;
(1)利用偏振装置采集某场景三幅不同透振方向的偏振图像,假设三幅已知透振方向的偏振图像的透振方向与参考方向的夹角分别为θ1、θ2、θ3,相应的偏振图像分别为I(θ1)、I(θ2)、I(θ3);
(2)对偏振图像进行配准
由于,同一场景多幅偏振图像的采集不是同一时刻的,场景中有运动目标或采集设备的震动会使得视场会轻微变化,因此采集的三幅图像并不严格一致,如果三幅图像不对准,将会给后面的计算stokes矢量带来误差,因此,需要先对偏振图像进行配准(如果成像系统非常稳定,此步骤可以省略);
假设I(θ1)为参考图像,I(θ2)和I(θ3)为待配准图像,相对于I(θ1),假设I(θ2)被平移(x0,y0),被旋转θ0角度;
①平移配准
平移算法流程图见图1所示,参考图像I(θ1)记为I1(x,y),待配准图像I(θ2)记为I2(x,y)即:
I2(x,y)=I1(x-x0,y-y0) (2)
则2幅图像对应的傅里叶变换有如下关系:
由傅里叶变换的平移性质知,信号在空间域的平移对应于其在频域的相位移,这两幅图像的相位差就等于它们互功率谱的相位,归一化互功率谱计算如下:
其中,FI1是I1的傅里叶变换,FI2是I2傅里叶变换的复共轭,通过对互功率谱进行傅里叶逆变换得到一个脉冲函数δ(x-x0,y-y0),该函数的极值点位置就是对应的平移量(x0,y0);
②旋转配准
旋转算法流程图见图2所示,I1(x,y)存在旋转,则要先计算出旋转角度并且得到补偿之后才能进行平移计算;
假设I2(x,y)是I1(x,y)旋转θ0后得到的图像,根据傅里叶变换的旋转性质,旋转一幅图像,在频域相当于对其傅里叶变换作相同角度的旋转,M1、M2为FI1、FI2功率谱,满足如下关系:
M2(x,y)=M1(xcosθ0+ysinθ0,-xsinθ0+ycosθ0) (5)
为了突出图像中的高频信息,对傅里叶变换的功率谱M1、M2取对数运算,通过双线性插值将其转换到极坐标下:
M2(ρ,θ)=M1(ρ,θ-θ0) (6)
这样把直角坐标下求旋转量转换为极坐标下求平移量,采用平移配准介绍的方法即可求得θ0;
(3)只要求解出偏振图像的Stokes矢量I、Q、U,将求解Stokes矢量带入式(1),就可求得任意透振方向的偏振图像I(θ);
对偏振图像的Stokes矢量I、Q、U的主要求解过程为:
采用分水岭算法对偏振图像进行过分割,假设在每个分割块内的所有点的Stokes矢量是相同的,利用最小二乘法,使误差平方和最小,从而求解Stokes矢量;
2)采用图像清晰度评价函数对合成的偏振图像进行评价;
用图像对比度来衡量图像是否清晰,统计各个透振方向的偏振图像的对比度,采用统计图像水平、垂直和对角方向梯度表示图像对比度,图像对比度公式如下:
3)确定合成的偏振图像的最佳透振方向;
根据图像清晰度评价函数,即公式(7)评价偏振图像的清晰度,选择杂散光滤除最多的,目标最清晰的偏振图像,被选中的最清晰的图像为:
该偏振图像对应的透振方向即为最佳透振方向。
在寻找最佳透振方向时,并不需要遍历所有的透振方向的偏振图像;可采取如图3所示的策略进行快速搜索。
Claims (1)
1.一种偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法,其特征在于,包括以下步骤:
1)利用任意三幅偏振镜的已知透振方向的偏振图像,合成任意透振方向的偏振图像;
当偏振图像的透振方向与参考方向夹角为 时,任意透振方向的偏振图像可表示为:
(1)
I、Q、U为偏振图像的Stokes矢量,代表偏振图像在透振方向时出射光的光强,I表示总光强,Q表示X轴方向与Y轴方向偏振图像分量之差,U表示±450方向偏振图像分量之差;
(1)利用偏振装置采集某场景三幅不同透振方向的偏振图像,假设三幅已知透振方向的偏振图像的透振方向与参考方向的夹角分别为,相应的偏振图像分别为;
(2)对偏振图像进行配准
由于,同一场景多幅偏振图像的采集不是同一时刻的,场景中有运动目标或采集设备的震动会使得视场会轻微变化,因此采集的三幅图像并不严格一致,如果三幅图像不对准,将会给后面的计算stokes矢量带来误差,因此,需要先对偏振图像进行配准,如果成像系统非常稳定,此步骤省略;
假设为参考图像,为待配准图像,相对于,假设被平移,被旋转角度;
①平移配准
参考图像记为,待配准图像记为即:
(2)
则2幅图像对应的傅里叶变换有如下关系:
(3)
由傅里叶变换的平移性质知,信号在空间域的平移对应于其在频域的相位移,这两幅图像的相位差就等于它们互功率谱的相位,归一化互功率谱计算如下:
(4)
其中,是的傅里叶变换,是傅里叶变换的复共轭,通过对互功率谱进行傅里叶逆变换得到一个脉冲函数,该函数的极值点位置就是对应的平移量;
②旋转配准
存在旋转,则要先计算出旋转角度并且得到补偿之后才能进行平移计算;
假设是旋转后得到的图像,根据傅里叶变换的旋转性质,旋转一幅图像,在频域相当于对其傅里叶变换作相同角度的旋转,、为、功率谱,满足如下关系:
(5)
为了突出图像中的高频信息,对傅里叶变换的功率谱、取对数运算,通过双线性插值将其转换到极坐标下:
(6)
这样把直角坐标下求旋转量转换为极坐标下求平移量,采用平移配准介绍的方法即可求得;
(3)只要求解出偏振图像的Stokes矢量I、Q、U,将求解Stokes矢量带入式(1),就可求得任意透振方向的偏振图像;
对偏振图像的Stokes矢量I、Q、U的主要求解过程为:
采用分水岭算法对偏振图像进行过分割,假设在每个分割块内的所有点的Stokes矢量是相同的,利用最小二乘法,使误差平方和最小,从而求解Stokes矢量;
2)采用图像清晰度评价函数对合成的偏振图像进行评价;
用图像对比度来衡量图像是否清晰,统计各个透振方向的偏振图像的对比度,采用统计图像水平、垂直和对角方向梯度表示图像对比度,图像对比度公式如下:
(7)
3)确定合成的偏振图像的最佳透振方向;
根据图像清晰度评价函数,即公式(7)评价偏振图像的清晰度,选择杂散光滤除最多的,目标最清晰的偏振图像,被选中的最清晰的图像为:
(8)
该偏振图像对应的透振方向即为最佳透振方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210144015.6A CN102663768B (zh) | 2012-05-10 | 2012-05-10 | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210144015.6A CN102663768B (zh) | 2012-05-10 | 2012-05-10 | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102663768A CN102663768A (zh) | 2012-09-12 |
CN102663768B true CN102663768B (zh) | 2014-11-05 |
Family
ID=46773246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210144015.6A Expired - Fee Related CN102663768B (zh) | 2012-05-10 | 2012-05-10 | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102663768B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105530419B (zh) * | 2016-02-05 | 2019-12-31 | 京东方科技集团股份有限公司 | 图像采集系统、图像采集处理系统和图像采集处理方法 |
CN110332994B (zh) * | 2019-07-22 | 2021-09-14 | 天津大学 | 基于偏振相机和优化波片的全Stokes矢量测量方法 |
CN110488349B (zh) * | 2019-08-20 | 2021-11-02 | 福建省建筑设计研究院有限公司 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
CN111028192B (zh) * | 2019-12-18 | 2023-08-08 | 维沃移动通信(杭州)有限公司 | 一种图像合成方法及电子设备 |
CN112672008B (zh) * | 2020-12-14 | 2022-09-06 | 杭州联吉技术有限公司 | 一种镜片的调整方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276465A (zh) * | 2008-04-17 | 2008-10-01 | 上海交通大学 | 广角图像自动拼接方法 |
CN101739674A (zh) * | 2008-11-19 | 2010-06-16 | 深圳迈瑞生物医疗电子股份有限公司 | 影像序列自动拼接方法与装置以及拼接系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4324238B2 (ja) * | 2007-02-13 | 2009-09-02 | パナソニック株式会社 | 画像処理システム、方法及び装置 |
-
2012
- 2012-05-10 CN CN201210144015.6A patent/CN102663768B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276465A (zh) * | 2008-04-17 | 2008-10-01 | 上海交通大学 | 广角图像自动拼接方法 |
CN101739674A (zh) * | 2008-11-19 | 2010-06-16 | 深圳迈瑞生物医疗电子股份有限公司 | 影像序列自动拼接方法与装置以及拼接系统 |
Non-Patent Citations (4)
Title |
---|
Automated registration of polarimetric imagery using Fourier transform techniques;Christopher M. Persons et al;《Proceedings of SPIE》;20021231;第4819卷;107-117页 * |
Christopher M. Persons et al.Automated registration of polarimetric imagery using Fourier transform techniques.《Proceedings of SPIE》.2002,第4819卷107-117页. * |
田光 等.相位相关法及其在偏振图像配准中的应用研究.《大气与环境光学学报》.2006,第1卷(第3期),193-197页. * |
相位相关法及其在偏振图像配准中的应用研究;田光 等;《大气与环境光学学报》;20061130;第1卷(第3期);193-197页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102663768A (zh) | 2012-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10630962B2 (en) | Systems and methods for object location | |
CN102663768B (zh) | 偏振视觉系统中自适应确定偏振镜的最佳透振方向的方法 | |
US10529076B2 (en) | Image processing apparatus and image processing method | |
CN106204595B (zh) | 一种基于双目摄像机的机场场面三维全景监视方法 | |
CN102313536B (zh) | 基于机载双目视觉的障碍物感知方法 | |
Shimizu et al. | Precise subpixel estimation on area‐based matching | |
CN104931966B (zh) | 一种基于dcs算法的星载视频sar成像处理方法 | |
CN108805910A (zh) | 多目车载记录仪、目标检测方法、智能驾驶系统及汽车 | |
CN102901489B (zh) | 路面积水积冰检测方法及装置 | |
WO2018128667A1 (en) | Systems and methods for lane-marker detection | |
CN103813099A (zh) | 一种基于特征点匹配的视频防抖方法 | |
CN102694963B (zh) | 一种获取无阴影目标图像的方法 | |
CN104535188A (zh) | 一种空间频率调制的静态全偏振成像探测系统及方法 | |
Morita et al. | Measurement in three dimensions by motion stereo and spherical mapping | |
CN103083089B (zh) | 数字立体显微系统的虚拟标尺方法和系统 | |
CN104167001A (zh) | 基于正交补偿的大视场摄像机标定方法 | |
Jin et al. | Practical speed measurement for an intelligent vehicle based on double radon transform in urban traffic scenarios | |
CN111243021A (zh) | 基于多组合相机的车载视觉定位方法、系统及存储介质 | |
Bergamasco et al. | A robust multi-camera 3D ellipse fitting for contactless measurements | |
CN112446915A (zh) | 一种基于图像组的建图方法及装置 | |
Lelégard et al. | Multiscale Haar transform for blur estimation from a set of images | |
Burschka et al. | Direct pose estimation with a monocular camera | |
Mares et al. | Vehicle self-localization in GPS-denied zones by multi-band imaging and analysis of prominent scene features | |
CN111489398B (zh) | 一种成像设备标定方法以及装置 | |
Kuhl et al. | Monocular 3D scene reconstruction at absolute scales by combination of geometric and real-aperture methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20180510 |
|
CF01 | Termination of patent right due to non-payment of annual fee |