CN104833425B - 一种偏振成像系统中偏振片角度的组合方法 - Google Patents
一种偏振成像系统中偏振片角度的组合方法 Download PDFInfo
- Publication number
- CN104833425B CN104833425B CN201510201708.8A CN201510201708A CN104833425B CN 104833425 B CN104833425 B CN 104833425B CN 201510201708 A CN201510201708 A CN 201510201708A CN 104833425 B CN104833425 B CN 104833425B
- Authority
- CN
- China
- Prior art keywords
- angle
- polaroid
- analyzing
- transmission direction
- time
- 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
Landscapes
- Studio Devices (AREA)
Abstract
本发明公开一种偏振成像系统中偏振片角度的组合方法,该方法包括:根据偏振成像系统的类型,从第一映射关系中选择线性偏振度误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;根据偏振成像系统的类型,从第二映射关系中选择偏振角误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;其中,第一映射关系为线性偏振度误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系;第二映射关系为偏振角误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系。本发明的偏振成像系统中偏振片角度的组合方法与现有技术相比能够明显减少线性偏振度近似值与偏振角近似值的噪声,提高探测偏振辐射时的信噪比。
Description
技术领域
本发明涉及偏振探测技术领域,具体涉及一种偏振成像系统中偏振片角度的组合方法。
背景技术
大多数情况下,尤其是在被动成像中,测量室外场景的发射辐射和反射辐射时,表示辐射的斯托克斯矢量中代表椭圆偏振成分的第四个分量可以忽略(参见文献[1]和[2])。由其他三个分量可以计算出辐射的线偏振度与偏振角,并用于3D重建与目标探测当中(参见文献[3]和[4])。
由三幅或以上的偏振图像可以计算得到斯托克斯矢量的前三个分量(参见文献[1]和[2])。文献[5]中基于主元素分析法,介绍了N通道偏振成像系统的最优线性组合法。另有多位相关研究人员提出了测量完全斯托克斯矢量及穆勒矩阵的优化方法(参见文献[5]至[7])。
常用线偏振成像系统按结构主要分为两种:一种是采用多个偏振方向固定的偏振片构成的多通道系统,可以一次性获得不同的偏振图像;另外一种是只有一片偏振片的单通道系统,通过旋转偏振片的偏振方向到不同的角度,获得相应的偏振图像,见附图1。对于采集三个不同角度的偏振图像的系统,称为三次检偏线偏振成像系统,采集四个不同角度图像的称为四次检偏线偏振成像系统。关于偏振片角度的优化组合方法中主要考虑线性偏振度(DoLP)和偏振角(AoP)的误差。文献[8]至[10]基于误差传播理论(参见文献[11])提出了表示斯托克斯矢量、线性偏振度以及偏振角噪声特性的方法。文献[9]基于误差传播理论提出了测量线性偏振度和偏振角时的最优通道数(Tyo于文献[12]中提及)。这些学者都指出了利用这种理论的局限性,但为了简化计算,他们并没有使用更合理的方法,而是通过限制计算范围以确保结果有效。
背景技术中涉及的参考文献如下:
[1]J.S.Tyo,D.L.Goldstein,D.B.Chenault,and J.A.Shaw,“Review of passiveimaging polarimetry for remote sensing applications,”Applied Optics 45(22),5453-5469(2006).
[2]J.R.Schott,“Fundamentals of polarimetric remote sensing,”(SPIEPress,2009).
[3]F.Sadjadi,and C.Chun,“Automatic detection of small objects fromtheir infrared state-of-polarization vectors,”Optics Letters 28(7),531-533(2003).
[4]D.Miyazaki,M.Kagesawa,and K.Ikeuchi,“Transparent surface modelingfrom a pair of polarization images,”IEEE Trans.Pattern Anal.MachineIntell.26,73–82(2004).
[5]J.S.Tyo,“Design of optimal polarimeters:maximization of signal-to-noise ratio and minimization of systematic error,”Applied Optics 41(4),619-630(2002)
[6]A.Ambirajan and D.C.Look,“Optimum angles for a polarimeter:part1and 2,”Opt.Eng.34(6),1651-1658(1995).
[7]D.S.Sabatke,M.R.Descour,“Optimization of retardance for a completeStokes polarimeter,”Optics Letters 25(11),802-804(2000).
[8]P.Robert,and V.Gruev.“Signal-to-noise analysis of Stokesparameters in division of focal plane polarimeters,”Optics express 18(25),25815-25824(2010).
[9]G.and A.Bénière,“Estimation precision of the degree oflinear polarization and of the angle of polarization in the presence ofdifferent sources of noise,”Applied optics 49(4),683-693(2010).
[10]W Feng,Y.Ji,and L.Chen."The impact of signal–noise ratio ondegree of linear polarization measurement,"Optik 124(3),192-194(2013).
[11]P.R.Bevington,and D.K.Robinson,“Data reduction and error analysisfor the physical sciences,”(McGraw-Hill,1969).
[12]J.S.Tyo,“Optimum linear combination strategy for an N-channelpolarization-sensitive imaging or vision system,”JOSA 15(2),359-366(1998).
发明内容
本发明针对现有线偏振成像系统测量线性偏振度和偏振角时优化方法的不足,提出一种偏振片透射方向的最优组合方法,在检偏次数一定的情况下,得到线偏振度和偏振角的误差具有最小值。
为此目的,本发明提出一种偏振成像系统中偏振片角度的组合方法,包括:
根据偏振成像系统的类型,从第一映射关系中选择线性偏振度误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
根据偏振成像系统的类型,从第二映射关系中选择偏振角误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
其中,所述偏振成像系统的类型包括:采用加性噪声相机的三次检偏线偏振成像系统、采用散粒噪声相机的三次检偏线偏振成像系统、采用加性噪声相机的四次检偏线偏振成像系统以及采用散粒噪声相机的四次检偏线偏振成像系统;
其中,第一映射关系为线性偏振度误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系;
第二映射关系为偏振角误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系。
可选的,所述第一映射关系具体包括:
采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtda,第三次检偏的偏振片的透射方向的角度为α+90°-βtda,其中α和βtda为预设常数;
采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtds,第三次检偏的偏振片的透射方向的角度为α+90°-βtds,其中α和βtds预设常数;
采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfda1,第二次检偏的偏振片的透射方向的角度为α-βfda1,第三次检偏的偏振片的透射方向的角度为α+90°+βfda2,第四次检偏的偏振片的透射方向的角度为α+90°-βfda2,其中α、βfda1和βfda2预设常数;
采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfda3,第四次检偏的偏振片的透射方向的角度为α+90°-βfda3,其中α和βfda3预设常数;
采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfds1,第二次检偏的偏振片的透射方向的角度为α-βfds1,第三次检偏的偏振片的透射方向的角度为α+90°+βfds2,第四次检偏的偏振片的透射方向的角度为α+90°-βfds2,其中α、βfds1和βfds2预设常数;
采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfds3,第四次检偏的偏振片的透射方向的角度为α+90°-βfds3,其中α和βfds3预设常数。
可选的,所述第二映射关系具体包括:
采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-45°,第二次检偏的偏振片的透射方向的角度为α+45°+βtaa,第三次检偏的偏振片的透射方向的角度为α+45°-βtaa;或第一次检偏的偏振片的透射方向的角度为α+45°,第二次检偏的偏振片的透射方向的角度为α-45°+βtaa,第三次检偏的偏振片的透射方向的角度为α-45°-βtaa,其中α和βtaa预设常数;
采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-γ1,第二次检偏的偏振片的透射方向的角度为α+γ1+βtas,第三次检偏的偏振片的透射方向的角度为α+γ1-βtas;或第一次检偏的偏振片的透射方向的角度为α+γ1,第二次检偏的偏振片的透射方向的角度为α-γ1+βtas,第三次检偏的偏振片的透射方向的角度为α-γ1-βtas,其中,γ1根据线性偏振度P确定,α和βtas预设常数;
采用加性噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+45°+βfaa,第二次检偏的偏振片的透射方向的角度为α+45°-βfaa,第三次检偏的偏振片的透射方向的角度为α-45°+βfaa,第四次检偏的偏振片的透射方向的角度为α-45°-βfaa,其中α和βfaa预设常数;
采用散粒噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+γ2+βfas,第二次检偏的偏振片的透射方向的角度为α+γ2-βfas,第三次检偏的偏振片的透射方向的角度为α-γ2+βfas,第四次检偏的偏振片的透射方向的角度为α-γ2-βfas,其中,γ2根据线性偏振度P确定,α和βfas预设常数。
相比于现有技术,本发明的偏振成像系统中偏振片角度的组合方法根据偏振成像系统的类型,在确定类型的基础上,即在检偏次数一定的情况下,针对线性偏振度误差的最小值和偏振角误差的最小值,通过变换不同偏振片的透射方向,得到不同偏振成像系统与偏振片角度的组合的对应关系。与传统的Pickering法、Fessenkov法以及改进的Pickering法(见参考文献[2])相比能够明显减少线性偏振度近似值与偏振角近似值的噪声,提高探测偏振辐射时的信噪比。
附图说明
图1为背景技术中只有一片偏振片的单通道偏振成像系统结构示意图;
图2为本发明实施例公开的一种偏振成像系统中偏振片角度的组合方法流程图;
图3a和3b为本发明实施例公开的使用加性噪声相机的偏振成像系统在不同θ1的情况下,变换θ2、θ3时线性偏振度近似值的标准差的最小值示意图;
图4a和4b为本发明实施例公开的使用加性噪声相机的偏振成像系统在不同θ1的情况下,变换θ2、θ3时偏振角近似值的标准差的最小值示意图;
图5为本发明实施例公开的使用加性噪声相机的偏振成像系统在线性偏振度增加部分βtda的值的示意图;
图6为本发明实施例公开的一种偏振成像系统中各参数计算流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图2所示,本实施例公开一种偏振成像系统中偏振片角度的组合方法,该方法可包括步骤201和202:
201、根据偏振成像系统的类型,从第一映射关系中选择线性偏振度误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
202、根据偏振成像系统的类型,从第二映射关系中选择偏振角误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
其中,所述偏振成像系统的类型包括:采用加性噪声相机的三次检偏线偏振成像系统、采用散粒噪声相机的三次检偏线偏振成像系统、采用加性噪声相机的四次检偏线偏振成像系统以及采用散粒噪声相机的四次检偏线偏振成像系统;
其中,第一映射关系为线性偏振度误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系;
其中,第二映射关系为偏振角误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系。
本实施例中,所述第一映射关系具体包括A1、B1、C1、D1、E1和F1:
A1:采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtda,第三次检偏的偏振片的透射方向的角度为α+90°-βtda,其中α和βtda为预设常数。
其中,第一次检偏的偏振片为检偏线偏振成像系统第一次检偏时对应的偏振片,第二次检偏的偏振片为检偏线偏振成像系统第二次检偏时对应的偏振片,第三次检偏的偏振片为检偏线偏振成像系统第三次检偏时对应的偏振片。
在具体应用中,α由场景辐射确定,即α为场景辐射的偏振角,本文中所有α的含义均相同。
B1:采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtds,第三次检偏的偏振片的透射方向的角度为α+90°-βtds,其中α和βtds预设常数。
C1:采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfda1,第二次检偏的偏振片的透射方向的角度为α-βfda1,第三次检偏的偏振片的透射方向的角度为α+90°+βfda2,第四次检偏的偏振片的透射方向的角度为α+90°-βfda2,其中α、βfda1和βfda2预设常数。
其中,第一次检偏的偏振片为检偏线偏振成像系统第一次检偏时对应的偏振片,第二次检偏的偏振片为检偏线偏振成像系统第二次检偏时对应的偏振片,第三次检偏的偏振片为检偏线偏振成像系统第三次检偏时对应的偏振片,第四次检偏的偏振片为检偏线偏振成像系统第四次检偏时对应的偏振片。
D1:采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfda3,第四次检偏的偏振片的透射方向的角度为α+90°-βfda3,其中α和βfda3预设常数。
E1:采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfds1,第二次检偏的偏振片的透射方向的角度为α-βfds1,第三次检偏的偏振片的透射方向的角度为α+90°+βfds2,第四次检偏的偏振片的透射方向的角度为α+90°-βfds2,其中α、βfds1和βfds2预设常数。
F1:采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfds3,第四次检偏的偏振片的透射方向的角度为α+90°-βfds3,其中α和βfds3预设常数。
本实施例中,所述第二映射关系具体包括A2、B2、C2和D2:
A2:采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-45°,第二次检偏的偏振片的透射方向的角度为α+45°+βtaa,第三次检偏的偏振片的透射方向的角度为α+45°-βtaa;或第一次检偏的偏振片的透射方向的角度为α+45°,第二次检偏的偏振片的透射方向的角度为α-45°+βtaa,第三次检偏的偏振片的透射方向的角度为α-45°-βtaa,其中α和βtaa预设常数;
B2:采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-γ1,第二次检偏的偏振片的透射方向的角度为α+γ1+βtas,第三次检偏的偏振片的透射方向的角度为α+γ1-βtas;或第一次检偏的偏振片的透射方向的角度为α+γ1,第二次检偏的偏振片的透射方向的角度为α-γ1+βtas,第三次检偏的偏振片的透射方向的角度为α-γ1-βtas,其中,γ1根据线性偏振度P确定,α和βtas预设常数;
C2:采用加性噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+45°+βfaa,第二次检偏的偏振片的透射方向的角度为α+45°-βfaa,第三次检偏的偏振片的透射方向的角度为α-45°+βfaa,第四次检偏的偏振片的透射方向的角度为α-45°-βfaa,其中α和βfaa预设常数;
D2:采用散粒噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+γ2+βfas,第二次检偏的偏振片的透射方向的角度为α+γ2-βfas,第三次检偏的偏振片的透射方向的角度为α-γ2+βfas,第四次检偏的偏振片的透射方向的角度为α-γ2-βfas,其中,γ2根据线性偏振度P确定,α和βfas预设常数。
在一个具体的例子中,使用加性噪声相机的线偏振成像系统,βtda,βtaa,βfda1,βfda2,βfda3,βfaa只与所测场景的线性偏振度以及E(S0)/σa相关。
使用散粒噪声相机的线偏振成像系统,βtds,βtas,βfds1,βfds2,βfds3,βfas只与所测场景的线偏振度以及E(S0)/(E(S0)/g)0.5相关。
S0为线偏振成像系统所测场景辐射的斯托克斯矢量的第一个分量,E(S0)为测得的S0的灰度值的期望,σa、g为偏振成像系统参数。
用SNR1表示E(S0)/σa或者SNR2表示E(S0)/(E(S0)/g)0.5。
与线性偏振度测量相关的βtda、βtds、βfda1、βfda2、βfda3、βfds1、βfds2和βfds3表示为:
β*d*=(SNR1)-1/3f*d*(DoLP)
与偏振角测量相关的βtaa、βtas、βfaa和βfas表示为:
β*a*=(SNR2)-1/2f*a*(DoLP)
其中,f*d*(DoLP)与f*a*(DoLP)皆为线性偏振度的函数。
对于使用加性噪声相机的线偏振成像系统,与线性偏振度测量相关的βtda、βfda1、βfda2、βfda3表示为:
β*da=(SNR1)-1/3f*da(DoLP)
对于使用加性噪声相机的线偏振成像系统,与偏振角测量相关的βtaa、βfaa表示为:
β*aa=(SNR1)-1/2f*aa(DoLP)
对于使用散粒噪声相机的线偏振成像系统,与线性偏振度测量相关的βtds、βfds1、βfds2、βfds3表示为:
β*ds=(SNR2)-1/3f*ds(DoLP)
对于使用散粒噪声相机的线偏振成像系统,与偏振角测量相关的βtas、βfas表示为:
β*as=(SNR2)-1/2f*as(DoLP)
其中,f***(DoLP)为与线性偏振度相关的变量,为预设参数。
γ1与γ2大致相等,实际环境中,场景辐射的线偏振度一般小于0.8,此时γ1或γ2可由下述经验公式进行计算:
γ=45+17.5*(DoLP)
上述偏振成像系统中偏振片角度的组合方法,主要考虑了线性偏振度和偏振角的测量误差。上述方法中第一映射关系和第二映射关系通过下述方式得到:由不同透射方向偏振片后被偏振成像系统探测到的辐射强度,计算偏振成像系统接收的场景辐射的斯托克斯矢量的前三个线性分量,并得到线性偏振度近似值与偏振角近似值测得的辐射强度的误差因此会传递到和中去,使用离散概率密度代替连续概率密度,考虑与待测信号强度无关的加性噪声与散粒噪声,在检偏次数一定的情况下,针对不同透射方向的偏振片组合,求取线性偏振度误差的最小值和偏振角误差的最小值,并将线性偏振度误差最小值和偏振角误差最小值对应的偏振片角度的组合分别作为第一映射关系和第二映射关系,即得到了不同偏振成像系统对应的偏振片角度的最优组合。
举例来说,用斯托克斯矢量的前三个分量S0、S1和S2构成的向量S=(S0,S1,S2)T表示线性偏振辐射。偏振成像系统中的探测器的响应在线性区域时,假设第n次检偏时偏振片滤光后所成图像的灰度值为In,其期望<In>用式(1)表示:
式中,θn是第n次检偏时偏振片透射方向的角度。偏振成像系统中的探测器的噪声接近于高斯分布时,经N次检偏滤光后所成图像的灰度值可以表示为:
I~N(WS,σ) (2)
式中,I=[I1,…,IN],σ是由偏振成像系统的噪声特性和入射辐射的强度的所决定的参数。W是一个N×3的矩阵:
在理想情况下,线性斯托克斯矢量可以由式(1)计算,不包括探测器引入的误差。当考虑类似于式(2)的扰动时,如果W的逆存在,S的近似值可以表示为:
如果W的逆不存在,需用广义逆矩阵代替W-1,从而式(4)可以表示为:
得到线性斯托克斯矢量后,DoLP与AoP的近似值可以表示为:
式中,arctan是四象限的反正切函数,函数值处于闭区间[-π,π]中。
假设X为随机向量,p(X1,…,XN)为X1,…,XN的联合概率密度函数,XN为X的第N个分量。y=f(X)是X向量的函数,此函数仅说明y和X有一定函数关系,对于具体的函数形式不作讨论,仅关注期望和方差,因此y的期望和方差可以表示为:
假设所有分量相互独立,联合概率密度函数可以表示为:
p(X1,…,XN)=p(X1)…p(XN) (10)
本实施例中X为图像的灰度值,各分量都处于零到最大灰度值gmax之间。图像的灰度值是离散的,因此必须对式(8)和(9)离散化,认为各通道采集图像的过程相互独立,式(8)和(9)可以改写为:
已知偏振成像系统的相机的噪声概率密度函数的情况下,将式(6)、(7)分别代入到式(11)、(12),可以得到线性偏振度近似值与偏振角近似值的期望和方差。
本实施例中考虑两类主要的噪声:加性噪声和泊松散粒噪声。
通常认为加性噪声呈高斯分布,即期望为0,方差为σa 2,本实施例中把连续区间(In–0.5,In+0.5)都认为是灰度值In,因此测得的第n次检偏时偏振片滤光后所成图像的灰度值In(期望为<In>)的概率密度函数可以写做:
散粒噪声是由信号决定的噪声,通常认为其服从泊松分布。当探测器处于线性响应范围时,图像的灰度值I与响应生成的电子数nphoton-electron成正比:
在只考虑泊松散粒噪声的情况下,测得In(期望为<In>)的概率密度函数可以写做:
式中erf为误差函数,erfc为余误差函数。
对于三次检偏线偏振成像系统,获得偏振片最佳组合的步骤为:首先,将θ1,θ2,θ3在[0°,180°]的范围内均分成48个不同的值,从而获得48×48×48组不同的角度组合,本实施例中分割成48个角度值只是为了使计算机仿真具有较快的速度与较高的精度,但并不只限于这种分割的方法;然后,利用式(11)和式(12)对不同组合计算出线性偏振度近似值与偏振角近似值的方差;最后重复计算出所有的结果比较得到偏振片角度的最佳组合。
以使用加性噪声相机为例。如图3a和3b所示,为不同θ1的情况下,变换θ2、θ3时线性偏振度近似值的标准差的最小值,图3a和3b中的辐射强度相同而偏振特性不同。三角箭头所指即为达到最小方差时的θ1、θ2、θ3最优组合,记为(θ1,θ2,θ3)1:
(θ1,θ2,θ3)1=(α,(α+90°)+βtda,(α+90°)-βtda) (16)
如图4a和4b所示,为不同θ1的情况下,变换θ2、θ3时偏振角近似值的标准差的最小值,此时有两种最优组合,分别用三角箭头和四角箭头标识,分别记为(θ1,θ2,θ3)2和(θ1,θ2,θ3)3。
当线性偏振度增加时βtda和βtaa都减小,部分βtda的值如图5所示。
使用散粒噪声相机时,能够获得具有最小方差的线性偏振度近似值与偏振角近似值的最佳组合模式表示为:
(θ1,θ2,θ3)4=(α,(α+90°)+βtds,(α+90°)-βtds) (18)
对于四次检偏线偏振成像系统,获得偏振片最佳组合的步骤与三次检偏的系统一样,只是由于有四次检偏,因此有48×48×48×48组不同的角度组合,本实施例中同样分割成48个角度值也只是为了使计算机仿真具有较快的速度与较高的精度,但并不只限于这种分割的方法。
使用加性噪声相机时,能够获得具有最小方差的线性偏振度近似值与偏振角近似值的最佳组合模式表示为:
(θ1,θ2,θ3,θ4)8=((α-45°)+βfaa,(α-45°)-βfaa,(α+45°)+βfaa,(α+45°)-βfaa) (21)
使用散粒噪声相机时,能够获得具有最小方差的线性偏振度近似值与偏振角近似值的最佳组合模式表示为:
(θ1,θ2,θ3,θ4)10=((α-γ2)+βfas,(α-γ2)-βfas,(α+γ2)+βfas,(α+γ2)-βfas) (23)
当使用一种确定的线偏振成像系统时(比如说基于加性噪声相机的三次检偏线偏振成像系统,σa=4),首先用传统的角度组合方法测量估计出场景辐射的线性斯托克斯矢量与场景辐射强度期望值E(S0),也就能算出与SNR的值,再按照前述方法变换偏振片的角度组合,计算出当线偏振度或偏振角具有最小方差时对应的偏振片的角度值,相应的β或γ的值也就由此确定,这就得到了在该场景条件下使用该系统的偏振片角度的最优组合。通过这种方法,可以事先计算出场景线偏振度不同且辐射强度不同的情况下,对应的计算β的经验公式中的可变系数f*d*(DoLP)或f*a*(DoLP),这样在实际应用中,首先用传统方法快速估测出场景的辐射强度、线偏振度及偏振角信息,就能插值得到f*d*(DoLP)或f*a*(DoLP)的值,并由前面的经验公式计算出β和γ的值,从而确定相应的偏振片角度的最优组合,相关流程图见图6。
如表1所示,为本发明的最优组合方法与使用传统方法时与的标准差的比值的示例,以Fessenkov法、改进的Pickering法为例,关于的标准差的比值分别表示为 关于的标准差的比值分别表示为可以看到比值皆小于0.8,因此运用最优组合法能够明显减少与的噪声。
表1使用最优组合方法与传统方法计算时与的标准差的比值
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方法,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种偏振成像系统中偏振片角度的组合方法,其特征在于,包括:
根据偏振成像系统的类型,从第一映射关系中选择线性偏振度误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
根据偏振成像系统的类型,从第二映射关系中选择偏振角误差最小值对应的所述偏振成像系统中偏振片角度的组合类型;
其中,所述偏振成像系统的类型包括:采用加性噪声相机的三次检偏线偏振成像系统、采用散粒噪声相机的三次检偏线偏振成像系统、采用加性噪声相机的四次检偏线偏振成像系统以及采用散粒噪声相机的四次检偏线偏振成像系统;
其中,第一映射关系为线性偏振度误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系;
第二映射关系为偏振角误差最小值时,所述偏振成像系统的类型与偏振片角度的组合的对应关系。
2.根据权利要求1所述的方法,其特征在于,所述第一映射关系具体包括:
采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtda,第三次检偏的偏振片的透射方向的角度为α+90°-βtda,其中βtda为预设常数,α为场景辐射的偏振角;
采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°+βtds,第三次检偏的偏振片的透射方向的角度为α+90°-βtds,其中βtds预设常数,α为场景辐射的偏振角;
采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfda1,第二次检偏的偏振片的透射方向的角度为α-βfda1,第三次检偏的偏振片的透射方向的角度为α+90°+βfda2,第四次检偏的偏振片的透射方向的角度为α+90°-βfda2,其中βfda1和βfda2预设常数,α为场景辐射的偏振角;
采用加性噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.27时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfda3,第四次检偏的偏振片的透射方向的角度为α+90°-βfda3,其中βfda3预设常数,α为场景辐射的偏振角;
采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P小于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+βfds1,第二次检偏的偏振片的透射方向的角度为α-βfds1,第三次检偏的偏振片的透射方向的角度为α+90°+βfds2,第四次检偏的偏振片的透射方向的角度为α+90°-βfds2,其中βfds1和βfds2预设常数,α为场景辐射的偏振角;
采用散粒噪声相机的四次检偏线偏振成像系统且线性偏振度P大于0.5时,对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α,第二次检偏的偏振片的透射方向的角度为α+90°,第三次检偏的偏振片的透射方向的角度为α+90°+βfds3,第四次检偏的偏振片的透射方向的角度为α+90°-βfds3,其中βfds3预设常数,α为场景辐射的偏振角。
3.根据权利要求1所述的方法,其特征在于,所述第二映射关系具体包括:
采用加性噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-45°,第二次检偏的偏振片的透射方向的角度为α+45°+βtaa,第三次检偏的偏振片的透射方向的角度为α+45°-βtaa;或第一次检偏的偏振片的透射方向的角度为α+45°,第二次检偏的偏振片的透射方向的角度为α-45°+βtaa,第三次检偏的偏振片的透射方向的角度为α-45°-βtaa,其中βtaa预设常数,α为场景辐射的偏振角;
采用散粒噪声相机的三次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α-γ1,第二次检偏的偏振片的透射方向的角度为α+γ1+βtas,第三次检偏的偏振片的透射方向的角度为α+γ1-βtas;或第一次检偏的偏振片的透射方向的角度为α+γ1,第二次检偏的偏振片的透射方向的角度为α-γ1+βtas,第三次检偏的偏振片的透射方向的角度为α-γ1-βtas,其中,γ1根据线性偏振度P确定,βtas预设常数,α为场景辐射的偏振角;
采用加性噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+45°+βfaa,第二次检偏的偏振片的透射方向的角度为α+45°-βfaa,第三次检偏的偏振片的透射方向的角度为α-45°+βfaa,第四次检偏的偏振片的透射方向的角度为α-45°-βfaa,其中βfaa预设常数,α为场景辐射的偏振角;
采用散粒噪声相机的四次检偏线偏振成像系统对应的偏振片角度的组合为:第一次检偏的偏振片的透射方向的角度为α+γ2+βfas,第二次检偏的偏振片的透射方向的角度为α+γ2-βfas,第三次检偏的偏振片的透射方向的角度为α-γ2+βfas,第四次检偏的偏振片的透射方向的角度为α-γ2-βfas,其中,γ2根据线性偏振度P确定,βfas预设常数,α为场景辐射的偏振角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510201708.8A CN104833425B (zh) | 2015-04-24 | 2015-04-24 | 一种偏振成像系统中偏振片角度的组合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510201708.8A CN104833425B (zh) | 2015-04-24 | 2015-04-24 | 一种偏振成像系统中偏振片角度的组合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104833425A CN104833425A (zh) | 2015-08-12 |
CN104833425B true CN104833425B (zh) | 2017-01-04 |
Family
ID=53811452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510201708.8A Active CN104833425B (zh) | 2015-04-24 | 2015-04-24 | 一种偏振成像系统中偏振片角度的组合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104833425B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297910B (zh) * | 2018-09-25 | 2021-05-18 | 中国人民解放军海军航空大学 | 一种分时成像偏振图像的偏振角度误差校正方法 |
CN112672008B (zh) * | 2020-12-14 | 2022-09-06 | 杭州联吉技术有限公司 | 一种镜片的调整方法及装置 |
CN113055566A (zh) * | 2021-02-23 | 2021-06-29 | 中国科学院合肥物质科学研究院 | 一种基于分时法的快速偏振成像装置及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441169B (zh) * | 2008-12-31 | 2010-09-01 | 合肥工业大学 | 平面四通道大气偏振信息检测传感器 |
-
2015
- 2015-04-24 CN CN201510201708.8A patent/CN104833425B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104833425A (zh) | 2015-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3371548B1 (en) | 3-d polarimetric imaging using a microfacet scattering model to compensate for structured scene reflections | |
CN106896069B (zh) | 一种基于彩色数码相机单幅rgb图像的光谱重建方法 | |
Hoekstra et al. | Masses of galaxy clusters from gravitational lensing | |
Huff et al. | Seeing in the dark–II. Cosmic shear in the Sloan Digital Sky Survey | |
CN111383285B (zh) | 一种基于毫米波雷达与摄像机传感器融合标定方法及系统 | |
CN102183214B (zh) | 一种大口径非球面镜结构光检测方法 | |
CN104793343A (zh) | 三通道单渥拉斯顿棱镜偏振成像装置及偏振信息探测方法 | |
AU2011242408B2 (en) | Shape and photometric invariants recovery from polarisation images | |
CN104833977A (zh) | 一种基于微波片阵列的瞬时遥感偏振成像装置及其实现方法 | |
CN104833425B (zh) | 一种偏振成像系统中偏振片角度的组合方法 | |
CN102509304A (zh) | 基于智能优化的摄像机标定方法 | |
CN102072706A (zh) | 一种多相机定位与跟踪方法及系统 | |
CN104363369A (zh) | 一种光场相机的图像恢复方法及装置 | |
CN106680812A (zh) | 一种基于解析面元的微波关联成像仿真方法 | |
CN103900696A (zh) | 一种仿螳螂虾视觉偏振拮抗感知的水下偏振成像方法 | |
Hayasaka et al. | Mobile visualization of density fields using smartphone background-oriented schlieren | |
CN116519257A (zh) | 基于单光场相机双视角背景纹影的三维流场测试方法及系统 | |
CN108318458B (zh) | 一种适用于不同天气条件下的室外典型地物pBRDF的测量方法 | |
Dong et al. | Non-iterative spot center location algorithm based on Gaussian for fish-eye imaging laser warning system | |
CN104813217B (zh) | 用于设计能够估计景深的无源单通道成像器的方法 | |
CN103873773B (zh) | 基于主辅协同双光路设计的全向成像方法 | |
CN104111080A (zh) | 一种凝视卫星面阵ccd相机的mtf在轨检测方法 | |
Lelégard et al. | Multiscale Haar transform for blur estimation from a set of images | |
CN109813429B (zh) | 一种偏振成像系统逐点标定方法 | |
CN106530351B (zh) | 一种基于图像传感器像素内量子效率获取的质心定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |