CN103345099B - 一种基于偏振光消光的成像方法 - Google Patents

一种基于偏振光消光的成像方法 Download PDF

Info

Publication number
CN103345099B
CN103345099B CN201310312020.8A CN201310312020A CN103345099B CN 103345099 B CN103345099 B CN 103345099B CN 201310312020 A CN201310312020 A CN 201310312020A CN 103345099 B CN103345099 B CN 103345099B
Authority
CN
China
Prior art keywords
pixel
window
gray
theta
image
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
Application number
CN201310312020.8A
Other languages
English (en)
Other versions
CN103345099A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201310312020.8A priority Critical patent/CN103345099B/zh
Publication of CN103345099A publication Critical patent/CN103345099A/zh
Application granted granted Critical
Publication of CN103345099B publication Critical patent/CN103345099B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Polarising Elements (AREA)

Abstract

本发明公开一种基于偏振光消光的成像方法,它解决了现有的在饱和光下对目标物体拍照成像清晰度低、不具有实时性的问题。其方法:控制拟被拍照物体到光源和相机连线的夹角为两倍的布儒斯特角,在相机镜头的前方设置一个垂直相机光轴的偏振片,通过旋转偏振片,分时拍摄多幅具有不同偏振角度的偏振图像,对偏振图像划定像素点窗口,将窗口内像素点灰度值与阀值进行比较,选取灰度值大于阀值的像素点个数最小的窗口内像素点融合成像。<!--1-->

Description

一种基于偏振光消光的成像方法
技术领域
本发明涉及一种基于偏振光消光的成像方法。
背景技术
由于光强太强或者物体镜面反射导致成像饱和,拍摄到的图像会成为刺眼的一片白色,不能分辨目标物体上的细节。现有方法一般是在相机镜头与物体之间增加一个衰减片以降低光强。但是,这样获得的图像整体光强变小,对比度也减弱,目标物体上的细节不能被完全分辨。
发明内容
本发明所要解决的技术问题是提供一种在饱和光下对目标物体通过偏振光消光的方法成像,从而提高目标物体成像的清晰度。
为了解决上述技术问题,本发明提供一种基于偏振光消光的成像方法,其步骤包括:
步骤一:控制被拍照物体(1)到光源(2)连线与被拍照物体(1)到相机镜头(3)连线所成的夹角为两倍的布儒斯特角;
步骤二:在相机镜头(3)的前方布置一个垂直相机光轴的偏振片(4),转动偏振片(4),选取m个不同偏振角度,拍摄m幅具有不同偏振角度的偏振图像,m≥2;
步骤三:用具有同一大小的像素点窗口将步骤二中所述的每幅偏振图像分别划分为k个像素点窗口,k≥2,在m幅偏振图像中形成k组、每组m个具有同一位置的对应像素点窗口;
步骤四:在步骤三中所述的每组具有同一位置的m个对应像素点窗口中,选定一个像素点窗口作为该组所在同一位置的图像融合窗口,所述选定方法是,该组中每个像素点窗口的全部像素点灰度值与事先设定的灰度阈值进行比较,计算该组中每个像素点窗口的像素点灰度值大于灰度阈值的像素点个数,选定像素点灰度值大于灰度阈值的像素点个数最小的窗口作为该组所在同一位置的图像融合窗口;
步骤五:将步骤四选定的k个图像融合窗口中的全部像素点融合生成完整图像。
本发明的原理是:自然光在电介质界面上反射和折射时,一般情况下反射光和折射光都是部分偏振光,只有当入射角为某特定角时反射光才是线偏振光,此规律被称为布儒斯特定律,此特定角称为布儒斯特角,布儒斯特角的大小由物体和物体所在介质的折射率决定。在日常应用中,需要使反射的线偏振光通过转动的偏振片获得具有特定偏振角度的线偏振光,这一过程就是偏振消光的过程。通过偏振消光获得多幅具有不同偏振角度的偏振图像,通过对不同偏振角度的偏振图像素点灰度值的比较,选取每幅图像最清晰的部分融合成最终图像。
本发明与现有技术相比,其显著优点在于通过偏振消光获得不同偏振角度的偏振图像进行图像融合,能够在不减弱反射光强的情况下获得清晰度高的图像。
附图说明
图1是偏振光消光原理图,其中标记θb是布儒斯特角,标记a是处于线偏振光状态的反射光,标记b是通过偏振片的具有特定偏振方向的线偏振光。
图2是本发明图像融合方法的流程示意图。
具体实施方式
实施例一、一种基于偏振光消光的成像方法,包括以下步骤:
步骤一:控制被拍照物体1到光源2和相机镜头3连线的夹角为两倍的布儒斯特角,所述布儒斯特角的大小由物体和物体所在介质的折射率决定;
步骤二:在相机镜头3的前方布置一个垂直相机光轴的偏振片4,转动偏振片4,选取m个不同偏振角度,拍摄m幅具有不同偏振角度的偏振图像,m≥2;
步骤三:用具有同一大小的像素点窗口将步骤二中所述的每幅偏振图像分别划分为k个像素点窗口,k≥2,在每幅偏振图像的相同区域形成k组,每组m个,具有同一位置的对应像素点窗口;
步骤四:在步骤三中所述的每组具有同一位置的m个对应像素点窗口中,选定一个像素点窗口作为该组所在同一位置的图像融合窗口,所述选定方法是,该组中每个像素点窗口的全部像素点灰度值与事先设定的灰度阀值进行比较,计算该组中每个像素点窗口的像素点灰度值大于阀值的像素点个数,选定像素点灰度值大于阀值的像素点个数最小的窗口作为该组所在同一位置的图像融合窗口。
步骤五:将步骤四选定的k个图像融合窗口中的全部像素融合生成完整图像。
实施例二、步骤一、步骤二与实施例一相同,步骤三、步骤四和步骤五的图像融合公式为:
D ( i , j ) = D &theta; 1 ( i , j ) i f ( s u m ( I &theta; 1 ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) D &theta; 2 ( i , j ) i f ( s u m ( I &theta; 2 ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) . . . D &theta; n ( i , j ) i f ( s u m ( I &theta; n ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) - - - ( 1 )
其中,
sum(Iθ1(i,j))表示偏振角度为θ1时对应图片第i行第j列像素点周围n×n窗口内灰度值大于阈值的像素点的个数,其中n必须为大于1的奇数,以此类推,
sum(Iθn(i,j))表示偏振角度为θn时对应图片第i行第j列像素点周围n×n窗口内灰度值大于阈值的像素点的个数,
阈值为250至255之间的任意一个实数,
Dθ1(i,j)表示偏振角度为θ1时对应图片第i行j列的灰度值,以此类推,
Dθn(i,j)表示偏振角度为θn时对应图片第i行j列的灰度值,
D(i,j)表示融合后的图像第i行j列的灰度值。
实施例三、步骤一至步骤五与实施例一相同,所述图像为彩色图像,对彩色图像按照公式(1)进行图像融合时,需要把每个像素分成R、G、B三通道像素,分别比较对应像素R、G、B的灰度值,然后把通过融合公式获得的R、G、B合成彩色像素点。
实施例四、步骤一至步骤五与实施例一相同,步骤四所述事先设定的灰度阈值为250至255之间任意实数,最优值为255。
实施例五、步骤一与实施例一相同;
步骤二:在相机镜头3的前方布置一垂直相机光轴的偏振片4,转动偏振片4,选取0°、60°和120°三个不同偏振角度,拍摄3幅偏振图像;
步骤三:用大小为3×3像素点窗口将步骤二中所述的每幅偏振图像分别划分为k个像素点窗口,k的数值由偏振图像大小和像素点窗口大小决定,在每幅偏振图像的相同区域形成k组,每组3个,具有同一位置的对应像素点窗口;
步骤四:在步骤三中所述的每组具有同一位置的3个对应像素点窗口中,选定一个像数点窗口作为该组所在相同区域的图像融合窗口,所述选定方法是,该组中每个窗口的全部像素点灰度值与事先设定的灰度阀值进行比较,计算改组中每个窗口的像素点灰度值大于阀值的像素点个数,选定像素点灰度值大于阀值的像素点个数最小的窗口作为该组所在相同区域的图像融合窗口。
步骤五:将步骤四选定的k个图像融合窗口中的全部像素融合生成完整图像。
步骤三至步骤五用融合公式表达为公式(2):
D ( i ) = I 0 ( i ) , i f ( s u m ( I 0 ( i ) ) < = s u m ( I 60 ( i ) ) , s u m ( I 0 ( i ) ) < = s u m ( I 120 ( i ) ) ) I 60 ( i ) , i f ( s u m ( I 60 ( i ) ) < = s u m ( I 0 ( i ) ) , s u m ( I 60 ( i ) ) < = s u m ( I 120 ( i ) ) ) I 120 ( i ) , i f ( s u m ( I 120 ( i ) ) < = s u m ( I 0 ( i ) ) , s u m ( I 120 ( i ) ) < = s u m ( I 60 ( i ) ) ) - - - ( 2 )
其中i是图像中某个像素点,I0(i)为偏振角为0°的图像在像素点i处的灰度值,I60(i)为偏振角为60°的图像在像素点i处的灰度值,I120(i)为偏振角120°的图像在像素点i处的灰度值,sum(I0(i))表示偏振角度为0°时对应图片第i行第j列像素点周围3×3的灰度值大于阈值的像素点的个数,(sum(I60(i))、sum(I120(i))依次类推,D(i)是用于最后合成图像的像素,将通过公式(2)计算获得的像素融合成最终图像。
实施例六、步骤一与实施例一相同;
步骤二:在相机镜头3的前方布置一垂直相机光轴的偏振片4,转动偏振片4,选取0°、45°、90°和135°四个不同偏振角度,拍摄4幅偏振图像;
步骤三:用大小为5×5像素点窗口将步骤二中所述的每幅偏振图像分别划分为k个像素点窗口,k的数值由偏振图像大小和像素点窗口大小决定,在每幅偏振图像的相同区域形成k组,每组4个,具有同一位置的对应像素点窗口;
步骤四:在步骤三中所述的每组具有同一位置的4个对应像素点窗口中,选定一个像数点窗口作为该组所在相同区域的图像融合窗口,所述选定方法是,该组中每个窗口的全部像素点灰度值与事先设定的灰度阀值进行比较,计算改组中每个窗口的像素点灰度值大于阀值的像素点个数,选定像素点灰度值大于阀值的像素点个数最小的窗口作为该组所在相同区域的图像融合窗口。
步骤五:将步骤四选定的k个图像融合窗口中的全部像素融合生成完整图像。
步骤三至步骤五用融合公式表达为公式(3):
D ( i ) = { I 0 ( i ) , i f ( s u m ( I 0 ( i ) ) < = s u m ( I 45 ( i ) ) , s u m ( I 0 ( i ) ) < = s u m ( I 90 ( i ) ) , s u m ( I 0 ( i ) ) < = s u m ( I 135 ( i ) ) ) I 45 ( i ) , i f ( s u m ( I 45 ( i ) ) < = s u m ( I 0 ( i ) ) , s u m ( I 45 ( i ) ) < = s u m ( I 90 ( i ) ) , s u m ( I 45 ( i ) ) < = s u m ( I 135 ( i ) ) ) I 90 ( i ) , i f ( s u m ( I 90 ( i ) ) < = s u m ( I 0 ( i ) ) , s u m ( I 90 ( i ) ) < = s u m ( I 45 ( i ) ) , s u m ( I 90 ( i ) ) < = s u m ( I 135 ( i ) ) ) I 135 ( i ) , i f ( s u m ( I 135 ( i ) ) < = s u m ( I 0 ( i ) ) , s u m ( I 135 ( i ) ) < = s u m ( I 45 ( i ) ) , s u m ( I 135 ( i ) ) < = s u m ( I 90 ( i ) ) ) - - - ( 3 )
偏振是各种矢量波共有的一种性质,通常选用电场强度E来定义偏振态,而描述偏振态经常使用斯托克斯矢量(Stokes)描述法。Stokes描述法是由I、Q、U、V四个参数来描述。Stokes可采用两种方法表示如公式(4)和公式(5)所示;因为自然中圆偏振很少,所以假设V=0。
I=I0+I90
Q=I0-I90
U=I45-I135(4)
由公式(4)和(5)可知采集偏振角为0°、60°和120°三幅图获得的效果和采集偏振角为0°、45°、90°和135°四幅图所获得的效果是等价的。

Claims (7)

1.一种基于偏振光消光的成像方法,其特征在于:
步骤一:控制被拍照物体(1)到光源(2)连线与被拍照物体(1)到相机镜头(3)连线所成的夹角为两倍的布儒斯特角;
步骤二:在相机镜头(3)的前方布置一个垂直相机光轴的偏振片(4),转动偏振片(4),选取m个不同偏振角度,拍摄m幅具有不同偏振角度的偏振图像,m≥2;
步骤三:用具有同一大小的像素点窗口将步骤二中所述的每幅偏振图像分别划分为k个像素点窗口,k≥2,在m幅偏振图像中形成k组、每组m个具有同一位置的对应像素点窗口;
步骤四:在步骤三中所述的每组具有同一位置的m个对应像素点窗口中,选定一个像素点窗口作为该组所在同一位置的图像融合窗口,所述选定方法是,该组中每个像素点窗口的全部像素点灰度值与事先设定的灰度阈值进行比较,计算该组中每个像素点窗口的像素点灰度值大于灰度阈值的像素点个数,选定像素点灰度值大于灰度阈值的像素点个数最小的窗口作为该组所在同一位置的图像融合窗口;
步骤五:将步骤四选定的k个图像融合窗口中的全部像素点融合生成完整图像。
2.根据权利要求1所述的基于偏振光消光的成像方法,其特征在于所述偏振图像为RGB色彩模式,所述像素点灰度值是红、绿、蓝三个颜色通道对应像素点的灰度值。
3.根据权利要求1或2所述的基于偏振光消光的成像方法,其特征在于步骤三所述像素点窗口大小为3×3。
4.根据权利要求1或2所述的基于偏振光消光的成像方法,其特征在于步骤四所述事先设定的灰度阈值为250至255之间的任意一个实数。
5.根据权利要求1或2所述的基于偏振光消光的成像方法,其特征在于步骤二所述m=3,偏振角度分别为0°、60°和120°。
6.根据权利要求1或2所述的基于偏振光消光的成像方法,其特征在于步骤二所述m=4,偏振角度分别为0°、45°、90°和135°。
7.根据权利要求1所述的基于偏振光消光的成像方法,其特征在于步骤四和步骤五的图像融合公式为:
D ( i , j ) = D &theta; 1 ( i , j ) i f ( s u m ( I &theta; 1 ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) D &theta; 2 ( i , j ) i f ( s u m ( I &theta; 2 ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) &CenterDot; &CenterDot; &CenterDot; D &theta; n ( i , j ) i f ( s u m ( I &theta; n ( i , j ) ) = min ( s u m ( I &theta; 1 ( i , j ) ) , s u m ( I &theta; 2 ( i , j ) ) ... s u m ( I &theta; n ( i , j ) ) ) ) - - - ( 1 )
其中,
sum(Iθ1(i,j))表示偏振角度为θ1时对应偏振图像第i行第j列像素点周围n×n窗口内灰度值大于灰度阈值的像素点的个数,其中n必须为大于1的奇数,以此类推,
sum(Iθn(i,j))表示偏振角度为θn时对应偏振图像第i行第j列像素点周围n×n窗口内灰度值大于灰度阈值的像素点的个数,
阈值为250至255之间的任意一个实数,
Dθ1(i,j)表示偏振角度为θ1时对应偏振图像第i行j列的灰度值,以此类推,
Dθn(i,j)表示偏振角度为θn时对应偏振图像第i行j列的灰度值,
D(i,j)表示融合后的偏振图像第i行j列的灰度值。
CN201310312020.8A 2013-07-24 2013-07-24 一种基于偏振光消光的成像方法 Expired - Fee Related CN103345099B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310312020.8A CN103345099B (zh) 2013-07-24 2013-07-24 一种基于偏振光消光的成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310312020.8A CN103345099B (zh) 2013-07-24 2013-07-24 一种基于偏振光消光的成像方法

Publications (2)

Publication Number Publication Date
CN103345099A CN103345099A (zh) 2013-10-09
CN103345099B true CN103345099B (zh) 2016-03-23

Family

ID=49279909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310312020.8A Expired - Fee Related CN103345099B (zh) 2013-07-24 2013-07-24 一种基于偏振光消光的成像方法

Country Status (1)

Country Link
CN (1) CN103345099B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104777624B (zh) * 2015-03-25 2017-07-04 清华大学深圳研究生院 一种偏振成像设备及方法
CN104902195A (zh) * 2015-04-02 2015-09-09 中国人民解放军63655部队 基于空中运动平台的红外偏振成像方法
CN106303396B (zh) * 2015-06-05 2019-03-05 南京理工大学 一种实时双目偏振消光探测方法
CN106292128A (zh) * 2015-06-05 2017-01-04 南京理工大学 一种基于偏振消光的成像方法
CN106303395B (zh) * 2015-06-05 2023-04-25 南京理工大学 一种实时双目偏振消光探测装置
CN107643292B (zh) * 2016-07-20 2021-05-04 南京造币有限公司 一种基于液晶玻璃的表面检测装置及其方法
CN106791339B (zh) * 2017-03-22 2019-12-13 北京小米移动软件有限公司 成像系统及其控制方法
CN107101948B (zh) * 2017-06-06 2019-12-27 南京理工大学 一种基于偏振吸收特性的海面溢油检测方法
CN107592465B (zh) * 2017-10-10 2020-05-26 联想(北京)有限公司 一种成像系统以及成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202024877U (zh) * 2011-03-16 2011-11-02 中国科学院上海技术物理研究所 一种标定偏振片透光轴的装置
CN102279052A (zh) * 2011-06-21 2011-12-14 中国科学院上海技术物理研究所 一种实时测量偏振光特性的方法
CN102401916A (zh) * 2011-11-30 2012-04-04 福建福晶科技股份有限公司 一种高消光比的偏振分光棱镜

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010121935A (ja) * 2007-11-05 2010-06-03 Nippon Sheet Glass Co Ltd 偏光画像撮像装置および画像処理装置
WO2008102562A1 (ja) * 2007-02-23 2008-08-28 Nippon Sheet Glass Company, Limited 偏光画像撮像装置、画像処理装置、偏光画像撮像方法および画像処理方法
US8072599B2 (en) * 2008-03-14 2011-12-06 Teledyne Scientific & Imaging, Llc Real-time, hybrid amplitude-time division polarimetric imaging camera

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202024877U (zh) * 2011-03-16 2011-11-02 中国科学院上海技术物理研究所 一种标定偏振片透光轴的装置
CN102279052A (zh) * 2011-06-21 2011-12-14 中国科学院上海技术物理研究所 一种实时测量偏振光特性的方法
CN102401916A (zh) * 2011-11-30 2012-04-04 福建福晶科技股份有限公司 一种高消光比的偏振分光棱镜

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
可见光偏振图像的自适应融合算法;曾海芳等;《光子学报》;20110131;第40卷(第1期);第132-135页 *
基于小波变换的红外偏振图像融合算法;虞文俊等;《激光技术》;20130531;第37卷(第3期);第289-292页 *
应用磁光调制器的高分辨率偏振消光测量系统;张建华;《光电子激光》;20011031;第12卷(第10期);第1041-1042页 *

Also Published As

Publication number Publication date
CN103345099A (zh) 2013-10-09

Similar Documents

Publication Publication Date Title
CN103345099B (zh) 一种基于偏振光消光的成像方法
CN101394573B (zh) 一种基于特征匹配的全景图生成方法及系统
CN106060493A (zh) 多源投影无缝边缘拼接方法及系统
CN103327220B (zh) 低照度Bayer图像上以绿色通道为导向的去噪方法
CN103384848A (zh) 用于确定与使用者佩戴的眼镜架内的眼镜片的安装相关联的至少一个几何/相貌参数的方法
CN108805839A (zh) 基于卷积神经网络的联合估计图像去雾方法
CN103763479A (zh) 实时高速高清全景视频的拼接装置及其方法
CN105245786A (zh) 一种基于智能测光的自拍方法、自拍系统及拍摄终端
CN109255758A (zh) 基于全1*1卷积神经网络的图像增强方法
CN103065334A (zh) 一种基于hsv颜色空间的偏色检测、校正方法及装置
CN114697623B (zh) 投影面选取和投影图像校正方法、装置、投影仪及介质
US11763432B2 (en) Multi-exposure image fusion method based on feature distribution weight of multi-exposure image
CN103886565A (zh) 基于目的性优化及直方图均衡化的夜间彩色图像增强方法
CN105222725B (zh) 一种基于光谱分析的高清图像动态采集方法
CN102867295A (zh) 一种彩色图像颜色校正方法
CN110223251A (zh) 适用于人工与自然光源的卷积神经网络水下图像复原方法
CN106331662A (zh) 一种图像撷取装置及图像撷取方法
CN104950476A (zh) 用于摄像装置的电控调光单元及摄像装置
CN103839245B (zh) 基于统计规律的Retinex夜间彩色图像增强方法
CN105678240B (zh) 一种用于道路去反光的图像处理方法
Wu et al. Underwater image restoration using color correction and non-local prior
CN104143203A (zh) 一种图像编辑传播方法
CN104683780A (zh) 一种视频监控摄像机的自动白平衡方法和装置
CN109035331B (zh) 一种信号灯组的位置矫正方法及设备
CN105427299A (zh) 一种基于畸变校正的像机焦距求解方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160323

Termination date: 20170724