CN105241450B - 基于四象限偏振片的天空偏振模式探测方法与*** - Google Patents
基于四象限偏振片的天空偏振模式探测方法与*** Download PDFInfo
- Publication number
- CN105241450B CN105241450B CN201510639950.3A CN201510639950A CN105241450B CN 105241450 B CN105241450 B CN 105241450B CN 201510639950 A CN201510639950 A CN 201510639950A CN 105241450 B CN105241450 B CN 105241450B
- Authority
- CN
- China
- Prior art keywords
- sky
- image
- polarization
- equal
- pixels
- 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.)
- Active
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 97
- 238000001514 detection method Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 15
- 238000003384 imaging method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000005855 radiation Effects 0.000 claims description 4
- 239000004973 liquid crystal related substance Substances 0.000 description 14
- 238000010586 diagram Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 2
- 239000002390 adhesive tape Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; 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
技术领域
本发明的应用属于天文导航技术领域,涉及一种利用四象限偏振片和手持式光场相机探测天空偏振模式的方法与***。
背景技术
利用天空光的偏振分布信息可以进行导航,它具有抗干扰能力强、成本较低等特点。为实现精确的天空偏振光导航,需要准确探测天空在较大视场范围内的偏振度和偏振化方向角信息。通常把偏振度和偏振化方向角信息一起称为偏振模式。目前用于大视场天空偏振模式探测的方法主要有两种,分别是基于光电二极管的方法和相机成像的方法。基于光电二极管的方法需要很多个光电二极管指向不同的天区测量,因而***结构复杂,实现也较为困难。相机成像的方法借助于广角镜头可对大视场的天空区域直接成像,但是为了测量偏振模式需要获得多幅不同偏振方向下的天空图像,这往往需要通过一个相机多次拍摄或者多个相机同时拍摄两种途径来实现,前一实现途径实时性不强,后一实现途径会导致***结构复杂。
总之,在实现大视场天空偏振模式探测时,现有的方法要么实时性差,要么***结构复杂。手持式光场相机的出现为实现经过一次成像即可测量大视场天空偏振模式提供了可能。手持式光场相机是一种在普通相机的光敏感元件前安装有微透镜阵列的新型成像器件,它由华裔科学家Ren Ng于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:
(公式一)
公式一中,I0°表示天空图像中所有与像素点(Xm,n-D,Ym,n-D)的距离小于等于k(1≤
k≤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。关于手持式光场相机结构的详细介绍参见Ren Ng的博士学位论文《Digital light field photography》。为了方便地使用***,将***获得的天空图像进行后续处理,可以给***再连接计算机12。计算机12负责相机控制、图像采集存储、图像处理和数据解算的工作。计算机12可以给手持式光场相机11发送控制信号13,控制信号13的箭头指示信号流向,手持式光场相机11也可以打包向计算机12发送图像数据包14,图像数据包14的箭头指示数据流向。
图2给出了本发明方法的工作流程图,其中第一步中将天空图像的灰度值转化为辐射光强值这一过程又称为相机响应曲线标定,使用Debevec和Malik所提出的方法(1997年ACM SIGGRAPH97会议论文,Recovering high dynamic range radiance maps fromphotographs),可得到天空偏振模式探测***拍摄得到的图像灰度值与拍摄目标辐射光强之间满足的光强响应曲线。天空偏振模式探测***搭建完毕后,只需要在第一次使用的时候实施一次相机响应曲线标定即可,以后使用时可直接省略该步骤,无需重复标定。在本步骤中手持式光场相机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:
公式一中,I0表示天空图像中所有与像素点(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和偏振化方向角矩阵θ即为需要的天空偏振模式测量结果。
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 CN105241450A (zh) | 2016-01-13 |
CN105241450B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106918845B (zh) * | 2017-05-04 | 2018-10-30 | 哈尔滨工业大学 | 基于偏振探测的远距离激光半主动方位测量装置 |
CN109059898B (zh) * | 2018-07-18 | 2020-05-22 | 大连理工大学 | 一种偏振光导航传感器及天空偏振光导航方法 |
CN111856770B (zh) * | 2019-04-26 | 2022-07-19 | 曹毓 | 一种偏振成像设备 |
CN113706398A (zh) * | 2020-05-22 | 2021-11-26 | 西北工业大学 | 在运动场景下生成高动态图像的装置及方法 |
CN113008423B (zh) * | 2021-02-24 | 2022-05-06 | 北京航空航天大学 | 一种玻璃幕墙应力检测方法 |
Citations (4)
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 | 中国科学院合肥物质科学研究院 | 图像式天空偏振光分布模式实时同步获取装置 |
-
2015
- 2015-09-30 CN CN201510639950.3A patent/CN105241450B/zh active Active
Patent Citations (4)
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 (2)
Title |
---|
天空偏振光分布的实验研究;关桂霞 等;《兵工学报》;20110430;第32卷(第4期);第459-463页 |
天空光偏振模式自动探测装置;赵开春 等;《光学精密工程》;20130228;第21卷(第2期);论文正文第3节 |
Also Published As
Publication number | Publication date |
---|---|
CN105241450A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105241450B (zh) | 基于四象限偏振片的天空偏振模式探测方法与*** | |
US10043290B2 (en) | Image processing to enhance distance calculation accuracy | |
EP3364643B1 (en) | Method for obtaining combined depth image, and depth camera | |
CN105827922B (zh) | 一种摄像装置及其拍摄方法 | |
CN103472592B (zh) | 一种快照式高通量的偏振成像方法和偏振成像仪 | |
CN106303228B (zh) | 一种聚焦型光场相机的渲染方法和*** | |
US10897608B2 (en) | Capturing light-field images with uneven and/or incomplete angular sampling | |
US9977981B2 (en) | Method and apparatus for calibrating a three-dimensional (3D) image in a tiled display | |
US10887576B2 (en) | Light field data representation | |
KR20170005009A (ko) | 3d 라돈 이미지의 생성 및 사용 | |
CN105488810A (zh) | 一种聚焦光场相机内外参数标定方法 | |
WO2021103872A1 (zh) | 图像传感器、摄像装置、电子设备和成像方法 | |
CN103592808A (zh) | 宽视场体积紧凑型偏振成像相机 | |
JPH04507301A (ja) | 高解像度三次元テレビジョン | |
CN105203102B (zh) | 基于s-波片的天空偏振模式探测方法与*** | |
CN113963065A (zh) | 一种基于外参已知的镜头内参标定方法及装置、电子设备 | |
CN203587939U (zh) | 宽视场体积紧凑型偏振成像相机 | |
US20170347046A1 (en) | Image sensor, image acquiring device comprising same, and portable terminal comprising the device | |
CN105300377B (zh) | 基于三象限偏振片的天空偏振模式探测方法与*** | |
CN105698751A (zh) | 一种测距***、测距方法、测距装置、以及拍摄装置 | |
WO2020181869A1 (zh) | 成像装置和电子设备 | |
CN108010071A (zh) | 一种利用3d深度测量的亮度分布测量***及方法 | |
US11882259B2 (en) | Light field data representation | |
CN105352499B (zh) | 基于多象限偏振片的天空偏振模式探测方法与*** | |
WO2017120640A1 (en) | Image sensor |
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 |