CN105241450A - 基于四象限偏振片的天空偏振模式探测方法与系统 - Google Patents

基于四象限偏振片的天空偏振模式探测方法与系统 Download PDF

Info

Publication number
CN105241450A
CN105241450A CN201510639950.3A CN201510639950A CN105241450A CN 105241450 A CN105241450 A CN 105241450A CN 201510639950 A CN201510639950 A CN 201510639950A CN 105241450 A CN105241450 A CN 105241450A
Authority
CN
China
Prior art keywords
sky
pixel
image
polarization
polaroid
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
CN201510639950.3A
Other languages
English (en)
Other versions
CN105241450B (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510639950.3A priority Critical patent/CN105241450B/zh
Publication of CN105241450A publication Critical patent/CN105241450A/zh
Application granted granted Critical
Publication of CN105241450B publication Critical patent/CN105241450B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Polarising Elements (AREA)
  • Studio Devices (AREA)

Abstract

本发明提供一种天空偏振模式探测系统和探测方法。技术方案一是:天空偏振模式探测系统包括镜头,四象限偏振片,手持式光场相机,四象限偏振片位于镜头光圈所在平面,由四个材料相同的偏振片组成。其中,相邻两个偏振片的偏振方向之间的夹角为45°。技术方案二是:天空偏振模式探测方法,包括下述步骤:第一步,利用天空偏振模式探测系统拍摄一幅天空图像,并将天空图像的灰度值转化为目标辐射的光强值。第二步,获得天空偏振度和偏振化方向角。本发明解决了现有方法应用于天空偏振模式探测时实时性差和系统结构复杂的问题。

Description

基于四象限偏振片的天空偏振模式探测方法与系统
技术领域
本发明的应用属于天文导航技术领域,涉及一种利用四象限偏振片和手持式光场相机探测天空偏振模式的方法与系统。
背景技术
利用天空光的偏振分布信息可以进行导航,它具有抗干扰能力强、成本较低等特点。为实现精确的天空偏振光导航,需要准确探测天空在较大视场范围内的偏振度和偏振化方向角信息。通常把偏振度和偏振化方向角信息一起称为偏振模式。目前用于大视场天空偏振模式探测的方法主要有两种,分别是基于光电二极管的方法和相机成像的方法。基于光电二极管的方法需要很多个光电二极管指向不同的天区测量,因而系统结构复杂,实现也较为困难。相机成像的方法借助于广角镜头可对大视场的天空区域直接成像,但是为了测量偏振模式需要获得多幅不同偏振方向下的天空图像,这往往需要通过一个相机多次拍摄或者多个相机同时拍摄两种途径来实现,前一实现途径实时性不强,后一实现途径会导致系统结构复杂。
总之,在实现大视场天空偏振模式探测时,现有的方法要么实时性差,要么系统结构复杂。手持式光场相机的出现为实现经过一次成像即可测量大视场天空偏振模式提供了可能。手持式光场相机是一种在普通相机的光敏感元件前安装有微透镜阵列的新型成像器件,它由华裔科学家RenNg于2005年发明。由于微透镜阵列的存在,手持式光场相机在对天空拍摄时每个微透镜都将某个特定视点下的天区成像为弥散光斑。将该弥散光斑按照不同区域分割利用,实现天空同一视点不同偏振方向下图像的获取,则可达到一次成像即完成大视场天空偏振模式探测的目的。
发明内容
发明要解决的技术问题是:本发明基于一种天空偏振模式探测系统,提供一种天空偏振模式探测方法。本发明通过单个手持式光场相机对大视场天空的一次成像,即实现了天空偏振模式的探测,用于解决现有方法在天空偏振模式探测时要么因需要多次测量而使得系统实时性差,要么因需要多套测量单元而导致系统结构复杂的问题。
本发明的技术方案一是:
一种天空偏振模式探测系统,包括镜头,四象限偏振片,手持式光场相机。其特征在于,四象限偏振片位于镜头光圈所在平面,由四个材料相同的偏振片组成。其中,相邻两个偏振片的偏振方向之间的夹角为45°。
本发明的技术方案二是:
一种天空偏振模式探测方法,利用技术方案一提供的天空偏振模式探测系统,具体包括下述步骤:
第一步,利用天空偏振模式探测系统拍摄一幅天空图像,并将天空图像的灰度值转化为目标辐射的光强值。
第二步,获得天空偏振度和偏振化方向角。
设天空偏振模式探测系统中手持式光场相机的微透镜阵列包括M行、N列微透镜;任意微透镜Lenm,n对应的成像圆斑中心的像素在图像坐标系下的坐标为(Xm,n,Ym,n),其中第一个坐标元素表示像素在图像矩阵中的行序号,第二个坐标元素表示像素在图像矩阵中的列序号。图像坐标系约定如下:图像左上角为坐标系的原点,图像中竖直向下的方向为第一个坐标元素的正向,图像中水平向右的方向为第二个坐标元素的正向。成像圆斑半径为d个像素,其中d≥5,1≤m≤M,1≤n≤N。本发明中所有微透镜的焦距相同,尺寸相同,成像圆斑直径也相同。
利用下式计算微透镜Em,n对应的天区偏振度Pm,n和偏振化方向角θm,n
(公式一)
公式一中,I表示天空图像中所有与像素点(Xm,n-D,Ym,n-D)的距离小于等于k(1≤k≤D)个像素的像素点的光强平均值;其中, D = [ 2 4 d ] , []表示四舍五入取整。
公式一中,I45°表示天空图像中所有与像素点(Xm,n-D,Ym,n+D)的距离小于等于k(1≤k≤D)个像素的像素点的光强平均值。
公式一中,I90°表示天空图像中所有与像素点(Xm,n+D,Ym,n+D)的距离小于等于k(1≤k≤D)个像素的像素点的光强平均值。
公式一中,I135°表示天空图像中所有与像素点(Xm,n+D,Ym,n-D)的距离小于等于k(1≤k≤D)个像素的像素点的光强平均值。
对所有的微透镜Len1,1,...Lenm,n...,LenM,N按照上述方法进行计算,从而得到所有微透镜对应视角的天区偏振度矩阵P(P1,1,...Pm,n,...PM,N)和偏振化方向角矩阵θ(θ1,1,...θm,n,...θM,N)。偏振度矩阵P和偏振化方向角矩阵θ即为需要的天空偏振模式测量结果。
本发明的有益效果是:利用手持式光场相机和在光圈位置加入四象限偏振片的镜头,通过对大视场天空的一次拍摄即得到一幅在任意视角下均包含四个偏振方向的天空图像,进而通过对获取图像的处理计算出天空任意视角的偏振度和偏振化方向角信息。因此,本发明通过单个相机对大视场天空的一次成像即实现了天空偏振模式的探测,解决了现有方法应用于天空偏振模式探测时实时性差和系统结构复杂的问题。
附图说明
图1为本发明提供的天空偏振模式探测系统原理示意图;
图2为具体实施流程图;
图3为四象限偏振片实物图;
图4为在镜头光圈处插入四象限偏振片的鱼眼镜头实物图;
图5为天空偏振模式探测系统实物图;
图6为室内实验场景图;
图7为天空偏振模式探测系统在对液晶显示器一次拍摄时获取的图像;
图8为图7局部放大的结果;
图9为液晶显示器的偏振化方向角分布的测量结果;
图10为液晶显示器的偏振度分布的测量结果;
具体实施方式
以下将结合说明书附图对本发明做进一步详细说明。
图1为本发明提供的天空偏振模式探测系统(以下简称系统)原理示意图。如图所示,系统包括镜头1,四象限偏振片3,手持式光场相机11。为了尽可能的获得较大的对天空成像范围,镜头1在本实施实例中建议选用视场角达到180°的鱼眼镜头。四象限偏振片3位于镜头1的光圈2所在位置。四象限偏振片3包括第一偏振片4,第二偏振片5,第三偏振片6和第四偏振片7。四象限偏振片3为圆形,第一偏振片4、第二偏振片5、第三偏振片6和第四偏振片7是四个位于同一平面的形状为扇形且尺寸相同、材料相同的偏振片,每个扇形的张角均为90°,其中双向箭头指示了它们各自的偏振方向,相邻象限偏振片的偏振方向之间的夹角为45°。需要说明的是,各偏振方向的值并不需要像图1所示的均为45°的整数倍,只需要满足相邻夹角为45°这一个条件即可。微透镜阵列8位于镜头1和相机的感光面10之间,微透镜阵列8由许多焦距相同,尺寸相同的微透镜排列而成,见微透镜阵列细节放大图9。微透镜阵列细节放大图9是微透镜按照正方形排布的结果,实际应用时也可以设计为按照六边形蜂窝状排布以更高效的利用相机的感光面10的空间。微透镜阵列8和相机感光面10组合在一起构成了手持式光场相机11。关于手持式光场相机结构的详细介绍参见RenNg的博士学位论文《Digitallightfieldphotography》。为了方便地使用系统,将系统获得的天空图像进行后续处理,可以给系统再连接计算机12。计算机12负责相机控制、图像采集存储、图像处理和数据解算的工作。计算机12可以给手持式光场相机11发送控制信号13,控制信号13的箭头指示信号流向,手持式光场相机11也可以打包向计算机12发送图像数据包14,图像数据包14的箭头指示数据流向。
图2给出了本发明方法的工作流程图,其中第一步中将天空图像的灰度值转化为辐射光强值这一过程又称为相机响应曲线标定,使用Debevec和Malik所提出的方法(1997年ACMSIGGRAPH97会议论文,Recoveringhighdynamicrangeradiancemapsfromphotographs),可得到天空偏振模式探测系统拍摄得到的图像灰度值与拍摄目标辐射光强之间满足的光强响应曲线。天空偏振模式探测系统搭建完毕后,只需要在第一次使用的时候实施一次相机响应曲线标定即可,以后使用时可直接省略该步骤,无需重复标定。在本步骤中手持式光场相机11拍摄一副天空图像时,应注意控制相机的曝光量,以尽可能缩小图像中过曝光或欠曝光区域大小,因为曝光不良的图像会显著降低本发明获取的天空偏振模式的精度。第二步获得天空偏振度和偏振化方向角,该步骤中图像坐标系的定义与《数字图像处理》一书中的定义相同(作者:冈萨雷斯,电子工业出版社2005年出版,第8页)。此外,为了第二步中描述的方便和对应公式的简洁,假定四象限偏振片3固定在相机1的光圈2中时,其中有两个象限之间的连接缝与相机成像传感器10的任一边缘平行或者垂直。
图3为我们自行加工的四象限偏振片3的实物图,四象限偏振片3的各个象限是从同一片完整的偏振片上按照不同的方向裁剪然后拼接在一起得到的。
图4为在光圈2处插入四象限偏振片3的镜头1实物。使用了双面胶将四象限偏振片3的边缘与镜头1中光圈2的边缘相互粘结固定,以避免实际使用中四象限偏振片3在光圈2中发生相对位置的移动。
图5为天空偏振模式探测系统实物图。
利用上述的方法和系统开展了室内的验证性实验。我们将液晶显示器作为天空偏振模式探测系统的测量对象。由于液晶固有的特征,液晶显示器发出的光可以认为是完全偏振光。图6为室内实验场景图。
图7给出了天空偏振模式探测系统在对液晶显示器一次拍摄时获取的图像。图8为图7局部放大的结果。图8中可以看到,每个微透镜所成的子像中有四个不同的亮度区域,这分别对应了四象限偏振片的四个偏振方向。
图9给出了天空偏振模式探测系统对液晶显示器的偏振化方向角分布的测量结果。图像中的灰度值表示了偏振化方向角的数值(单位是度),可见对于液晶显示器的偏振化方向角的测量结果比较均匀,与真实情况吻合,其数值约为0°,与实际测量得到的液晶显示器真实的偏振化方向角数值吻合。约定沿着显示器竖直边缘向上的方向为0°参考方向。图中箭头所指示为液晶显示器的偏振化方向角测量结果数值对应在右侧灰度指示条上的位置。
图10给出了天空偏振模式探测系统对液晶显示器的偏振度分布的测量结果。图像中的灰度值表示了偏振度的数值(偏振度没有单位),可见对于液晶显示器的偏振度的测量结果同样比较均匀,其数值约为0.9,这与真实情况是较为吻合的。图中箭头所指示为液晶显示器的偏振化方向角测量结果数值对应在右侧灰度指示条上的位置。

Claims (2)

1.一种天空偏振模式探测系统,包括镜头、四象限偏振片、手持式光场相机,其特征在于,四象限偏振片位于镜头光圈所在平面,由四个材料相同的偏振片组成;其中,相邻两个偏振片的偏振方向之间的夹角为45°。
2.一种天空偏振模式探测方法,利用权利要求1提供的天空偏振模式探测系统,具体包括下述步骤:
第一步,利用天空偏振模式探测系统拍摄一幅天空图像,并将天空图像的灰度值转化为目标辐射的光强值;
第二步,获得天空偏振度和偏振化方向角:
设天空偏振模式探测系统中手持式光场相机的微透镜阵列包括M行、N列微透镜,所有微透镜均相同;任意微透镜Lenm,n对应的成像圆斑中心的像素在图像坐标系下的坐标为(Xm,n,Ym,n),其中第一个坐标元素表示像素在图像矩阵中的行序号,第二个坐标元素表示像素在图像矩阵中的列序号;图像坐标系定义如下:图像左上角为坐标系的原点,图像中竖直向下的方向为第一个坐标元素的正向,图像中水平向右的方向为第二个坐标元素的正向;成像圆斑半径为d个像素,其中d≥5,1≤m≤M,1≤n≤N;
利用下式计算微透镜Em,n对应的天区偏振度Pm,n和偏振化方向角θm,n
公式一中,I表示天空图像中所有与像素点(Xm,n-D,Ym,n-D)的距离小于等于k个像素的像素点的光强平均值,1≤k≤D;其中,[]表示四舍五入取整;
公式一中,I45°表示天空图像中所有与像素点(Xm,n-D,Ym,n+D)的距离小于等于k个像素的像素点的光强平均值;
公式一中,I90°表示天空图像中所有与像素点(Xm,n+D,Ym,n+D)的距离小于等于k个像素的像素点的光强平均值;
公式一中,I135°表示天空图像中所有与像素点(Xm,n+D,Ym,n-D)的距离小于等于k个像素的像素点的光强平均值;
对所有的微透镜Len1,1,...Lenm,n...,LenM,N按照上述方法进行计算,从而得到所有微透镜对应视角的天区偏振度矩阵P=(P1,1,...Pm,n,...PM,N)和偏振化方向角矩阵θ=(θ1,1,...θm,n,...θM,N),偏振度矩阵P和偏振化方向角矩阵θ即为需要的天空偏振模式测量结果。
CN201510639950.3A 2015-09-30 2015-09-30 基于四象限偏振片的天空偏振模式探测方法与系统 Active CN105241450B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510639950.3A CN105241450B (zh) 2015-09-30 2015-09-30 基于四象限偏振片的天空偏振模式探测方法与系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510639950.3A CN105241450B (zh) 2015-09-30 2015-09-30 基于四象限偏振片的天空偏振模式探测方法与系统

Publications (2)

Publication Number Publication Date
CN105241450A true CN105241450A (zh) 2016-01-13
CN105241450B CN105241450B (zh) 2019-02-05

Family

ID=55039180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510639950.3A Active CN105241450B (zh) 2015-09-30 2015-09-30 基于四象限偏振片的天空偏振模式探测方法与系统

Country Status (1)

Country Link
CN (1) CN105241450B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918845A (zh) * 2017-05-04 2017-07-04 哈尔滨工业大学 基于偏振探测的远距离激光半主动方位测量装置
CN109059898A (zh) * 2018-07-18 2018-12-21 大连理工大学 一种偏振光导航传感器及天空偏振光导航方法
CN111856770A (zh) * 2019-04-26 2020-10-30 曹毓 一种偏振成像设备
CN113008423A (zh) * 2021-02-24 2021-06-22 北京航空航天大学 一种玻璃幕墙应力检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1971246A (zh) * 2006-05-09 2007-05-30 中国科学院合肥物质科学研究院 大气偏振模式检测装置及其检测方法
CN102967311A (zh) * 2012-11-30 2013-03-13 中国科学院合肥物质科学研究院 基于天空偏振分布模型匹配的导航定位方法
US8530817B1 (en) * 2008-10-27 2013-09-10 Lockheed Martin Corporation Field of view limit detection enhancement for quadrant-array detector systems
CN103335714A (zh) * 2013-07-09 2013-10-02 中国科学院合肥物质科学研究院 图像式天空偏振光分布模式实时同步获取装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1971246A (zh) * 2006-05-09 2007-05-30 中国科学院合肥物质科学研究院 大气偏振模式检测装置及其检测方法
US8530817B1 (en) * 2008-10-27 2013-09-10 Lockheed Martin Corporation Field of view limit detection enhancement for quadrant-array detector systems
CN102967311A (zh) * 2012-11-30 2013-03-13 中国科学院合肥物质科学研究院 基于天空偏振分布模型匹配的导航定位方法
CN103335714A (zh) * 2013-07-09 2013-10-02 中国科学院合肥物质科学研究院 图像式天空偏振光分布模式实时同步获取装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIJUAN SU ET AL.: "A Light Field Imaing Polarimeter Using a Division of Focal-Plane Device", 《COMPUTATIONAL OPTICAL SENSING AND IMAGING》 *
关桂霞 等: "天空偏振光分布的实验研究", 《兵工学报》 *
赵开春 等: "天空光偏振模式自动探测装置", 《光学精密工程》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918845A (zh) * 2017-05-04 2017-07-04 哈尔滨工业大学 基于偏振探测的远距离激光半主动方位测量装置
CN109059898A (zh) * 2018-07-18 2018-12-21 大连理工大学 一种偏振光导航传感器及天空偏振光导航方法
CN109059898B (zh) * 2018-07-18 2020-05-22 大连理工大学 一种偏振光导航传感器及天空偏振光导航方法
CN111856770A (zh) * 2019-04-26 2020-10-30 曹毓 一种偏振成像设备
CN113008423A (zh) * 2021-02-24 2021-06-22 北京航空航天大学 一种玻璃幕墙应力检测方法

Also Published As

Publication number Publication date
CN105241450B (zh) 2019-02-05

Similar Documents

Publication Publication Date Title
CN106989824B (zh) 一种红外测温成像装置及检测方法
CN105488810B (zh) 一种聚焦光场相机内外参数标定方法
CN105678742B (zh) 一种水下相机标定方法
CN110211043B (zh) 一种用于全景图像拼接的基于网格优化的配准方法
CN102692214B (zh) 一种狭窄空间双目视觉测量定位装置及方法
CN101630406B (zh) 摄像机的标定方法及摄像机标定装置
CN105654502A (zh) 一种基于多镜头多传感器的全景相机标定装置和方法
CN105716542B (zh) 一种基于柔性特征点的三维数据拼接方法
CN104463949A (zh) 一种基于光场数字重聚焦的快速三维重建方法及其系统
US10021340B2 (en) Method and an apparatus for generating data representative of a light field
CN106886979A (zh) 一种图像拼接装置及图像拼接方法
CN101271187B (zh) 无死角的双目立体全方位视觉传感装置
CN105241450A (zh) 基于四象限偏振片的天空偏振模式探测方法与系统
CN202172446U (zh) 一种广角拍照装置
CN109325981B (zh) 基于聚焦像点的微透镜阵列型光场相机几何参数标定方法
CN107481288A (zh) 双目摄像头的内外参确定方法和装置
KR20180053668A (ko) 화소 빔을 나타내는 데이터를 생성하기 위한 방법 및 장치
CN107991665A (zh) 一种基于定焦相机的对标靶三维坐标连续测量方法
CN107610182A (zh) 一种光场相机微透镜阵列中心的标定方法
Anderson et al. Omnidirectional real time imaging using digital restoration
CN102628693A (zh) 一种用于摄像机主轴与激光束进行平行配准的方法
CN105203102B (zh) 基于s-波片的天空偏振模式探测方法与系统
CN107527323B (zh) 镜头畸变的标定方法及装置
CN105300377A (zh) 基于三象限偏振片的天空偏振模式探测方法与系统
CN105352499A (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
GR01 Patent grant
GR01 Patent grant