CN112161707A - 一种微偏振阵列设计模型及其排布模式获取方法 - Google Patents
一种微偏振阵列设计模型及其排布模式获取方法 Download PDFInfo
- Publication number
- CN112161707A CN112161707A CN202011019076.0A CN202011019076A CN112161707A CN 112161707 A CN112161707 A CN 112161707A CN 202011019076 A CN202011019076 A CN 202011019076A CN 112161707 A CN112161707 A CN 112161707A
- Authority
- CN
- China
- Prior art keywords
- micro
- polarization
- array
- design
- arrangement mode
- 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
Links
- 238000013461 design Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000010287 polarization Effects 0.000 claims abstract description 46
- 238000003384 imaging method Methods 0.000 claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims description 26
- 238000001914 filtration Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001276 controlling effect Effects 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000003491 array Methods 0.000 abstract description 13
- 238000001514 detection method Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 8
- 230000003287 optical effect Effects 0.000 abstract description 8
- 238000011161 development Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/447—Polarisation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/28—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for polarising
- G02B27/286—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for polarising for controlling or changing the state of polarisation, e.g. transforming one polarisation state into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
- G01J4/04—Polarimeters using electric detection means
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种微偏振阵列设计模型及基于该模型的阵列排布模式获取方法,属于偏振光学器件领域,主要涉及偏振光学技术、分焦平面型偏振成像探测技术、傅里叶频域光学分析技术等。该模型,可针对2×N系列和N×N系列的微偏振阵列进行结构设计,通过改变设计模型的参数a,b,c,优化二维频域中各个组成成分之间的距离和位置,增加组成成分间的不受串扰影响的带宽范围,以此来提升斯托克斯参数图像及DoLP图像的重构质量。该设计模型从二维频率域分析的角度出发,简化了微偏振阵列排布模式的设计及获取步骤,拓展了阵列排布模式设计范围,为高性能的微偏振阵列的研制提供了理论依据,实现了高质量斯托克斯参数图像重构,促进了分焦平面型偏振成像系统的实际应用。
Description
所属技术领域
本发明属于偏振光学器件领域,主要涉及偏振光学技术、分焦平面型偏振成像探测技术、傅里叶频域光学分析技术等。
背景技术
偏振成像技术与强度成像、光谱成像等技术相比,具有独特的目标识别优势:除了获取传统目标成像信息外,还能够额外获取目标及背景的偏振度(DoLP)、偏振角(AoP)等有效信息,可以凸显目标物体的轮廓边缘与细节信息,提高目标识别精度与概率,相对于其它探测手段在目标探测识别方面具有无法替代的优越性。现有的各类偏振成像设备中,分焦平面型偏振成像系统实现了微偏振阵列芯片与探测器的像素级集成,可以对不同偏振状态的入射光进行采集,具有结构紧凑、集成度高等优点,适用于静态和动态场景目标探测成像应用。该成像系统的组成参阅图1,主要包括目标1、前置透镜组2、微偏振阵列结构3和探测器阵列4。作为分焦平面型偏振成像系统的核心光学元件,微偏振阵列的排布模式对偏振成像性能的影响成为当前研究的热点之一。
1994年C.S.LChun等人(Proc.SPIE,1994,2234:275.)研制了第一个基于金属微纳米光栅的2×2排布模式的微偏振阵列,实现了场景目标的线性偏振度测量,自此分焦平面型偏振成像系统均采用这种典型的2×2阵列排布模式。然而,这种典型的阵列结构存在三个不足:第一个参阅图2所示,通过2×2阵列得到的图像在频域5内有三个频率分量6,在一个周期的频域内各频率分量会与中心位置的频率分量产生严重的混叠,混叠区域为7,进而影响斯托克斯参数图像(s0、s1、s2)、及DoLP图像的重构质量;第二个是不同方向角度的微偏振片会“遮蔽”大部分的入射光,存在空间分辨率损失的问题,限制了其在高分辨成像和低光照成像中的应用;第三是基于2×2排布模式的微偏振阵列在重构全分辨率斯托克斯参数图像时存在瞬时视场(IFOV)误差,而去除瞬时视场误差不仅与重构算法有关还与微偏振阵列的排布模式有关。为了提升分焦平面型偏振成像系统中斯托克斯参数图像的重构质量,Le Master等人(Opt.Lett.,2014,39(7):1811.)设计出了2×4排布模式的新型微偏振阵列,对应的频谱图参阅图3所示,通过优化频域内各个频率成分之间的距离来增加斯托克斯参数图像不受串扰影响的带宽范围,实验结果表明:利用2×4阵列重构的斯托克斯参数图像质量有了较大提升。2015年赵永强等人(Chin.Control Conf.IEEE,2015,pp.4574.)在频域内设计彩色滤波阵列方法的启发下,提出了频域内微偏振阵列结构的设计方法,并设计了5×5排布模式的微偏振阵列,通过仿真实验研究分析不同的排布模式对偏振信息获取性能的影响。2017年Alenin等人(Opt.Lett.,2017,42(3):458.)在传统2×2和2×4排布模式的基础上建立了2×N排布系列的微偏振阵列设计模型,通过傅里叶频域分析扩大频域内的带宽范围来优化设计了能产生最佳带宽的新型2×N微偏振阵列结构,实验结果表明:设计的新型2×N阵列能减少不同频率成分间的混叠,其重构的斯托克斯参数图像质量较传统的2×2微偏振阵列和2×4微偏振阵列有了进一步提升。
上述的科研人员都指出了传统的2×2微偏振阵列排布模式存在的局限性,并且都基于傅里叶频域设计了新型的阵列结构,但尚未对更丰富的微偏振阵列排布模式进行优化设计和性能比较。作为分焦平面型偏振成像系统的核心光学器件,新型微偏振阵列结构的探索和研究是十分必要的。
发明内容
发明目的
本发明目的在于针对典型的2×2微偏振阵列排布模式存在频域混叠误差及图像重构质量低等缺点,以及现有的设计模型仅能设计2×N系列的微偏振阵列结构的局限性,提出一种可涵盖2×N系列和N×N系列微偏振阵列排布模式的设计模型,通过优化二维频域中各个组成成分的位置和带宽范围,设计新型微偏振阵列结构,旨在提升斯托克斯参数图像及DoLP图像重构质量,满足偏振成像领域对低瞬时视场误差和高偏振图像重构质量的要求。
技术方案
为了对本发明专利提出的微偏振阵列结构设计模型及设计方法进行说明,需要对微偏振阵列的成像方式进行合理的表示。
斯托克斯矢量S=[s0、s1、s2、s3]可以表示入射光的偏振状态,对于分焦平面型偏振成像系统,偏振成像探测器是通过微偏振器阵列来间接获取斯托克斯矢量,可由如下公式计算求解:
其中,I0°,I45°,I90°,和I135°表示通过偏振角度分别为0°,45°,90°,和135°方向上的测量强度,表示椭圆偏振成分的分量s3在线偏振成像探测系统中可以忽略,入射光的DoLP和AoP可以由其他三个分量s0,s1和s2计算求解。为了获取场景目标的偏振态信息,采用多个方向角度固定的微偏振片构成的微偏振阵列结构,可以对未知偏振态的入射光进行间接测量,一次性获得不同角度的偏振图像,公式如下:
I(x,y)=SA(x,y)·Si(x,y)+N(x,y)
其中,I(x,y)表示经过微偏振阵列的光强,Si(x,y)表示入射光的偏振状态,i=0,1,2,SA(x,y)为微偏振阵列排布模式的数学表达式,即结构设计模型,N(x,y)表示探测系统的噪声。同时,经过角度为θ的偏振片的光强度可以表示为:
Iθ=(s0+s1 cos2θ+s2 sin 2θ)/2
在分焦平面型偏振成像系统中,Si(x,y)随着探测的场景目标变化而变化,但通过微偏振阵列得到图像的二维频谱所包含的频率成分以及这些成分所在的二维位置和相对距离是不变的,而这三个方面均由微偏振阵列的排布模式所决定。即,SA(x,y)完全由微偏振阵列的排布模式决定。
下面结合说明书附图来说明本发明的具体内容。
本发明专利提出的微偏振阵列结构设计模型,不仅可以设计2×N排布模式的微偏振阵列结构,而且可以设计N×N对角排布模式的微偏振阵列,其具体的SA(x,y)表示式如下所示:
其中,x,y分别表示微偏振阵列结构中的水平和垂直轴的坐标,a,b,c为微偏振阵列排布模式的设计参数,表示通过微偏振阵列得到的图像对应的二维频谱图中不同频率分量的位置信息,通过调控这三个参数可以改变频率分量的位置及距离,进而实现不同阵列结构的设计。现有的2×N和新型的N×N微偏振阵列的具体设计参数及对应的排布模式如下表所示:
所述阵列结构类型参阅图4(a)为2×3阵列结构,图4(b)为2×4阵列结构,图4(c)为2×7阵列结构,图4(d)为2×2×2阵列结构,图4(e)为3×3阵列结构,图4(f)为5×5阵列结构,图4(g)为7×7阵列结构。
需要注意的是,2×2×2微偏振阵列不适应于快照式偏振成像探测系统。
进一步的,当设计2×N排布模式的微偏振阵列时,所述设计模型中的参数a、b满足:a=N,b=N/2,参数c可为c=1,也可为c=2、3…(N-1)。
进一步的,当设计N×N对角排布模式的微偏振阵列时,所述设计模型中的参数a应满足:a=N,参数b、c可为b=c=1,也可为b=c=2、3…(N-1)。
本发明提出的微偏振阵列结构优化设计基本原理是:对所提出的微偏振阵列排布模式设计模型SA(x,y)进行傅里叶变换,可以得到如下公式:
其中,ξ和η分别表示水平和竖直方向的频域坐标。
从上式可以看出,通过微偏振阵列得到的图像对应的二维频谱图,在一个频谱周期中包含0.5S0频率分量,(S1-iS2)频率分量和(S1+iS2)频率分量,且该频谱在二维空间上周期重复;在二维频谱的一个周期中,通过改变模型参数a,b,c,优化二维频域中各个组成成分之间的距离和位置,增加组成成分间的不受串扰影响的带宽范围,以此来提升斯托克斯参数图像重构质量。该设计模型从二维频率域分析的角度出发,简化了微偏振阵列排布模式的设计步骤,拓展了阵列排布模式设计范围,具有理论指导意义。
本发明专利提出一种基于设计模型进行2×N排布系列和N×N对角排布系列的微偏振阵列结构排布模式获取方法,其基本设计流程参阅图5,具体包括如下基本步骤:
步骤1:选定频谱图中频率分量的空间位置;
微偏振阵列对应的二维频谱周期为(1,1),在该二维频谱的一个周期中包含不同的频率分量,其中频率分量0.5S0的位置坐标为(0,0),设定频率分量(S1-iS2)的位置坐标为(x0,y0),频率分量(S1+iS2)的位置坐标为(x1,y1)。根据分析可知频率分量0.5S0在水平和竖直方向上能量较高,为了减小频率分量间的混叠程度,设计性能优于典型2×2阵列的微偏振阵列结构,其余两个频率分量的中心位置应当避免位于水平轴和垂直轴上,并且位置坐标分别满足:
步骤2:选定微偏振阵列结构的排布模式;
根据所选定的频谱结构形式,选定微偏振阵列结构的排布模式。其中,除去中心位置处的频率分量,当|x0|=|x1|,且|y0|=|y1|=1/2时,其余频率分量呈水平轴和垂直轴对称分布,所设计的阵列结构为2×N排布模式,其中N=1/|x0|=1/|x1|;当|x0|=|x1|=|y0|=|y1|≠1/2时,其余频率分量呈中心对称分布,所设计的阵列结构为N×N对角排布模式,且N=1/|x0|=1/|y0|=1/|x1|=1/|y1|。
步骤3:确定微偏振阵列结构的设计参数;
根据提出的设计模型及选定的阵列结构排布模式可知,当设计最小周期为2×N排布系列的微偏振阵列时,阵列结构的设计模型参数为:a=N,b=N/2,c=1。
当设计最小周期为N×N对角排布系列的微偏振阵列时,阵列结构的设计模型参数为:a=N,b=c=N*|x0|=N*|y0|。
步骤4:确定微偏振阵列结构最小周期内不同位置的偏振片角度;
根据微偏振阵列结构的设计参数,结合本发明提出的设计模型表达公式SA(x,y),可以得到所设计的微偏振阵列结构表达公式为:
其中,
根据公式:
Iθ=(s0+s1 cos2θ+s2 sin2θ)/2
即可求解出微偏振阵列结构最小周期内的不同位置处的偏振片角度θ。
步骤5:获取微偏振阵列的排布模式;
将所述得到的最小周期微偏振阵列结构周期延拓,即可实现所设计的微偏振阵列结构排布模式。
基于获得的微偏振阵列排布模式对应的频谱结构,还能够得到一种在频率域上重构斯托克斯参数图像及DoLP图像的方法。例如,根据传统的2×2微偏振阵列的频谱结构可知,进行目标物体偏振信息图像重构的最直接方法,就是先通过设计合适带宽的滤波器进行频域滤波得到各频率成分的频率信息,然后借助傅里叶反变换以及变换矩阵求解即可:
基于微偏振阵列的频谱结构,本发明专利给出一种在频率域上重构斯托克斯参数图像及DoLP图像的方法,此方法对本发明提出的模型所设计的任意微偏振阵列适用,具体流程参阅图6,包括如下基本步骤:
步骤1:根据设计参数,确定微偏振阵列结构的排布模式,得到经所述微偏振阵列偏振成像后的马赛克式图像;
步骤2:对所述微偏振阵列偏振成像后的马赛克式图像进行离散傅里叶变换,得到各频率分量的系数及各频率分量所处的频率调制位置;
步骤3:设计合适的频率滤波器在所述得到的各频率分量所处的频率调制位置进行带通滤波,得到所述的各频率分量的频谱信息;
步骤4:利用变换矩阵,对得到的各频率分量进行求解,得到所述图像对应的斯托克斯参数图像的频谱;
步骤5:对所述得到的斯托克斯参数图像的频谱进行傅里叶反变换,即可得到重构的斯托克斯参数s0、s1、s2图像,进一步可获得DoLP图像。
有益效果
相比于现有的微偏振阵列设计模型和方法,本发明专利提出微偏振阵列设计模型,可针对2×N系列和N×N系列的微偏振阵列进行结构设计,通过改变设计模型的参数a,b,c,优化二维频域中各个组成成分之间的距离和位置,增加组成成分间的不受串扰影响的带宽范围,以此来提升斯托克斯参数图像及DoLP图像的重构质量。该设计模型从二维频率域分析的角度出发,简化了微偏振阵列排布模式的设计及获取步骤,拓展了阵列排布模式设计范围,为高性能的微偏振阵列的研制提供了理论依据,实现了高质量斯托克斯参数图像重构,促进了分焦平面型偏振成像系统的实际应用。其具体的有益效果总结如下:
(1)本发明突破了现有设计模型仅能设计2×N系列的微偏振阵列结构的局限性,提出一种可涵盖2×N系列和N×N系列微偏振阵列排布模式的设计模型;
(2)本发明的设计模型是基于频率域进行多种类型微偏振阵列结构设计,简化了阵列排布模式的设计及获取步骤;
(3)利用本发明设计的新型3×3、5×5和7×7对角排布微偏振阵列重构的斯托克斯参数s0、s1、s2图像和DoLP图像质量均优于典型的2×2对角排布微偏振阵列;对于快照式偏振成像系统,新型3×3对角排布微偏振阵列重构的斯托克斯参量s1、s2和DoLP图像质量优于现有的2×N排布系列的微偏振阵列;新型7×7对角排布和2×3微偏振阵列重构的斯托克斯参量s0图像质量最优,具体参阅图8(a)。
附图说明
图1为分焦平面型偏振成像系统示意图
图2为经过典型2×2微偏振阵列的图像对应的二维频谱结构示意图
图3为经过2×4微偏振阵列的图像对应的二维频谱结构示意图
图4为利用本发明提出的设计模型获取的不同排布模式的微偏振阵列结构图
图5为利用本发明提出的设计模型进行微偏振阵列结构设计流程图
图6为基于频域滤波的斯托克斯参数图像重构流程图
图7为本发明实施例提出的3×3微偏振阵列结构设计方法
图8为本发明设计模型获取的不同类型微偏振阵列在不同信噪比噪声下的图像重构性能对比图
其中1.场景目标;2.前置透镜组;3.微偏振阵列结构;4.探测器阵列;5.傅里叶频域;6.频率分量;7.混叠区域。
具体实施方式
本发明实施例给出一个利用设计模型进行3×3对角排布的微偏振阵列排布模式获取的例子。该阵列结构由0°、60°和120°方向角度的微偏振片组成,其中SA如下所示:
对所提出的3×3对角排布的微偏振阵列结构进行傅里叶变换,可以得到如下公式:
由上式可以看出,该3×3阵列对应的二维频谱图中避免了两个频率分量位于水平和竖直方向上,增大了各频率分量之间的距离,减小了混叠误差,可以提升斯托克斯参数图像及DoLP图像的重构质量。
本发明实施的一种3×3对角排布的微偏振阵列排布模式获取,其基本设计流程参阅图5,包括如下基本步骤:
步骤1:选定频谱图中频率分量的空间位置;
参阅图7(a),微偏振阵列对应的二维频谱周期为(1,1),在该二维频谱的一个周期中包含三个频率分量,其中频率分量0.5S0的中心坐标为(0,0),为了减小频率分量间的混叠程度,选取频率分量(S1-iS2)的位置坐标为(-1/3,-1/3),频率分量(S1+iS2)的位置坐标为(1/3,1/3)。根据分析可知频率分量0.5S0在水平和竖直方向上能量较高,为了减小频率分量间的混叠程度,设计性能优于典型2×2阵列的微偏振阵列结构,选定的其余两个频率分量的中心位置避免了位于水平轴和垂直轴上,并且位置坐标分别满足:
步骤2:选定微偏振阵列结构的排布模式;
根据所选定的频谱结构形式,选定微偏振阵列结构的排布模式。根据二维频谱图中频率分量的分布,当除去中心位置处的频率成分时,其余频率成分呈中心对称分布,则选定该微偏振阵列结构为3×3对角排布模式。
步骤3:确定微偏振阵列结构的设计参数;
根据提出的设计模型及选定的阵列结构排布模式可知,当设计最小周期为3×3对角排布系列的微偏振阵列时,设计模型参数满足:a=3,b=c=1。
步骤4:确定微偏振阵列结构最小周期内不同位置的偏振片角度;
参阅图7(b),根据微偏振阵列结构的设计参数,结合本发明提出的设计模型表达公式SA(x,y),可以得到3×3对角排布的微偏振阵列结构表达式为:
那么,
根据公式:
Iθ=(s0+s1 cos2θ+s2 sin2θ)/2
可以求解出不同位置处的偏振片角度为:
θ(0,0)=0° | θ(0,1)=60° | θ(0,2)=120° |
θ(1,0)=60° | θ(1,1)=120° | θ(1,2)=0° |
θ(2,0)=120° | θ(2,1)=0° | θ(2,2)=60° |
步骤5:获取微偏振阵列的排布模式;
参阅图7(c),将所述得到的3×3最小周期排布结构进行周期延拓,即可获取3×3对角微偏振阵列的排布模式。
为了获取所设计的微偏振阵列结构对应的斯托克斯参数及DoLP图像,本发明实施例基于频率域滤波给出了一种偏振图像重构方法,对经过3×3对角排布的微偏振阵列得到的图像做离散傅里叶变换,得到所述图像对应的二维频谱图,该频谱在二维空间上周期重复,且在任意一个周期中,包含可以求解出斯托克斯参数的频率分量;随后利用带通滤波器进行频域滤波,重构斯托克斯参数图像及DoLP图像,其基本设计流程参阅图6,具体包括以下步骤:
步骤1:获取经过3×3对角排布的微偏振阵列成像后的马赛克式图像;
步骤2:对马赛克式图像进行离散傅里叶变换,得到3×3阵列对应的频域内各频率分量的系数及其所处的频率调制位置,各频率分量的系数及位置坐标如下所示:
步骤3:利用频域滤波方法,设计合适的频率滤波器在所述得到的各组成成分所处的频率调制位置进行带通滤波,得到所述各频率分量的频谱信息;
步骤4:利用如下变换矩阵,对得到的各频率分量进行求解,得到所述图像对应的斯托克斯参数s0、s1、s2的频谱;
步骤5:对所述得到的斯托克斯参数图像的频谱进行傅里叶反变换,即可得到重构的斯托克斯参数s0、s1、s2图像,进一步可获得DoLP图像。
为了进一步验证利用本发明专利提出的模型所设计的微偏振阵列的实用性、有效性和先进性,通过对原始图像添加不同程度的高斯噪声,进行偏振图像重构性能比较,其中添加的噪声信噪比(SNR)以2dB间隔从0dB增加至60dB。偏振图像重构质量用峰值信噪比(PSNR)来表示,其中PSNR值越高,表示重构图像质量越好。参阅图8所示,噪声SNR值较低时,不同类型的微偏振阵列抗噪性能相近;随着SNR值的增加,不同类型的阵列重构的偏振图像的PSNR随着SNR的增加而增加;当噪声SNR值较高时,PSNR随着SNR的增加而近乎不变。此时,通过比较PSNR值可知,2×3和7×7阵列重构的s0图像质量最高,2×2×2阵列重构的s1、s2和DoLP图像质量最高,而对于快照式偏振成像系统,3×3对角排布的微偏振阵列重构的s1、s2和DoLP图像优于现有的2×N系列阵列。因此,本发明专利提出的模型所设计的新型3×3阵列具有良好的抗噪性能,能有效保证重构图像的清晰度和细节信息,综合表现最好。
Claims (3)
1.一种微偏振阵列设计模型,其特征在于,所述微偏振阵列的SA(x,y)表示式如下所示:
其中,x,y分别表示微偏振阵列结构中的水平和垂直轴的坐标,a,b,c为微偏振阵列排布模式的设计参数,表示通过微偏振阵列得到的图像对应的二维频谱图中不同频率分量的位置信息,通过调控这三个参数可以改变频率分量的位置及距离,进而实现不同阵列结构的设计,则2×N和N×N微偏振阵列的具体设计参数及对应的排布模式如下表所示:
当设计类型为2×N排布模式的微偏振阵列时,所述设计模型中的参数a、b满足:a=N,b=N/2,c=1或c=2、3…(N-1);
当设计类型为N×N对角排布模式的微偏振阵列时,所述设计模型中的参数a应满足:a=N,参数b=c1或者2、3…(N-1)。
2.一种基于如权利要求1所述的模型进行微偏振阵列结构排布模式获取方法,其特征在于,具体包括如下基本步骤:
步骤1:选定频谱图中频率分量的空间位置;
微偏振阵列对应的二维频谱周期为(1,1),在该二维频谱的一个周期中包含不同的频率分量,其中频率分量0.5S0的位置坐标为(0,0),设定频率分量(S1-iS2)的位置坐标为(x0,y0),频率分量(S1+iS2)的位置坐标为(x1,y1)。根据分析可知频率分量0.5S0在水平和竖直方向上能量较高,为了减小频率分量间的混叠程度,设计性能优于典型2×2阵列的微偏振阵列结构,其余两个频率分量的中心位置应当避免位于水平轴和垂直轴上,并且位置坐标分别满足:
步骤2:选定微偏振阵列结构的排布模式;
根据所选定的频谱结构形式,选定微偏振阵列结构的排布模式。其中,除去中心位置处的频率分量,当|x0|=|x1|,且|y0|=|y1|=1/2时,其余频率分量呈水平轴和垂直轴对称分布,所设计的阵列结构为2×N排布模式,其中N=1/|x0|=1/|x1|;当|x0|=|x1|=|y0|=|y1|≠1/2时,其余频率分量呈中心对称分布,所设计的阵列结构为N×N对角排布模式,且N=1/|x0|=1/|y0|=1/|x1|=1/|y1|;
步骤3:确定微偏振阵列结构的设计参数;
根据提出的设计模型及选定的阵列结构排布模式可知,当设计最小周期为2×N排布系列的微偏振阵列时,阵列结构的设计模型参数为:a=N,b=N/2,c=1。
当设计最小周期为N×N对角排布系列的微偏振阵列时,阵列结构的设计模型参数为:a=N,b=c=N*|x0|=N*|y0|。
步骤4:确定微偏振阵列结构最小周期内不同位置的偏振片角度;
根据微偏振阵列结构的设计参数,结合本发明提出的设计模型表达公式SA(x,y),可以得到所设计的微偏振阵列结构表达公式为:
其中,
根据公式:
Iθ=(s0+s1cos2θ+s2sin2θ)/2
即可求解出微偏振阵列结构最小周期内的不同位置处的偏振片角度θ。
步骤5:获取微偏振阵列的排布模式;
将所述得到的最小周期微偏振阵列结构周期延拓,即可实现所设计的微偏振阵列结构排布模式。
3.一种在频率域上重构斯托克斯参数图像及DoLP图像的方法,此方法对本发明提出的模型所设计的任意微偏振阵列适用,具体流程参阅图6,包括如下基本步骤:
步骤1:根据设计参数,确定微偏振阵列结构的排布模式,得到经所述微偏振阵列偏振成像后的马赛克式图像;
步骤2:对所述微偏振阵列偏振成像后的马赛克式图像进行离散傅里叶变换,得到各频率分量的系数及各频率分量所处的频率调制位置;
步骤3:设计合适的频率滤波器在所述得到的各频率分量所处的频率调制位置进行带通滤波,得到所述的各频率分量的频谱信息;
步骤4:利用变换矩阵,对得到的各频率分量进行求解,得到所述图像对应的斯托克斯参数图像的频谱;
步骤5:对所述得到的斯托克斯参数图像的频谱进行傅里叶反变换,即可得到重构的斯托克斯参数s0、s1、s2图像,进一步可获得DoLP图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011019076.0A CN112161707B (zh) | 2020-09-25 | 2020-09-25 | 一种微偏振阵列及其排布模式获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011019076.0A CN112161707B (zh) | 2020-09-25 | 2020-09-25 | 一种微偏振阵列及其排布模式获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112161707A true CN112161707A (zh) | 2021-01-01 |
CN112161707B CN112161707B (zh) | 2023-04-14 |
Family
ID=73862849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011019076.0A Active CN112161707B (zh) | 2020-09-25 | 2020-09-25 | 一种微偏振阵列及其排布模式获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112161707B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106802184A (zh) * | 2017-02-22 | 2017-06-06 | 西安交通大学 | 快照式紧凑噪声免疫型光场成像全偏振光谱探测装置及方法 |
CN107014490A (zh) * | 2017-04-17 | 2017-08-04 | 大连理工大学 | 一种分焦面型实时偏振成像系统 |
CN108007575A (zh) * | 2017-11-29 | 2018-05-08 | 中国科学院光电技术研究所 | 一种基于二元数字编码双折射晶体的小型化全斯托克斯矢量偏振成像装置 |
CN111562223A (zh) * | 2019-03-25 | 2020-08-21 | 上海昊量光电设备有限公司 | 一种基于微偏振片阵列的偏振成像装置及方法 |
-
2020
- 2020-09-25 CN CN202011019076.0A patent/CN112161707B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106802184A (zh) * | 2017-02-22 | 2017-06-06 | 西安交通大学 | 快照式紧凑噪声免疫型光场成像全偏振光谱探测装置及方法 |
CN107014490A (zh) * | 2017-04-17 | 2017-08-04 | 大连理工大学 | 一种分焦面型实时偏振成像系统 |
CN108007575A (zh) * | 2017-11-29 | 2018-05-08 | 中国科学院光电技术研究所 | 一种基于二元数字编码双折射晶体的小型化全斯托克斯矢量偏振成像装置 |
CN111562223A (zh) * | 2019-03-25 | 2020-08-21 | 上海昊量光电设备有限公司 | 一种基于微偏振片阵列的偏振成像装置及方法 |
Non-Patent Citations (1)
Title |
---|
ANDREY S. ALENIN等: "《Optimal bandwidth micropolarizer arrays》", 《OPTICS LETTERS》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112161707B (zh) | 2023-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112180329B (zh) | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 | |
US20120076362A1 (en) | Coded aperture camera with adaptive image processing | |
CN106803892A (zh) | 一种基于光场测量的光场高清晰成像方法 | |
US20220270362A1 (en) | Method for processing hyperspectral images | |
CN111380612A (zh) | 一种高光谱成像系统 | |
CN107765221A (zh) | 适用于相干和非相干声源的反卷积声源成像算法 | |
CN111031264A (zh) | 一种基于透射式红外孔径编码成像系统及其超分辨方法 | |
EP3700197A1 (en) | Imaging device and method, and image processing device and method | |
CN110333189B (zh) | 基于压缩感知原理的光子集成干涉成像高分辨重构方法 | |
CN112161707B (zh) | 一种微偏振阵列及其排布模式获取方法 | |
CN110837090A (zh) | 一种星载高光谱红外图像干涉纹波的校正方法 | |
CN115685518A (zh) | 基于远心光学系统的反射式傅里叶叠层成像装置及方法 | |
CN109443250A (zh) | 一种基于s变换的结构光三维面形垂直测量方法 | |
CN110307781B (zh) | 干涉成像系统 | |
CN112097904A (zh) | 基于衍射透镜/变焦透镜阵列的光谱成像系统及方法 | |
CN111476738A (zh) | 一种高光谱图像Etalon效应纹波的分离方法 | |
CN1837763A (zh) | 采用菲涅尔双面镜的全反射式傅立叶变换成像光谱仪 | |
CN107664533B (zh) | 一种三狭缝高光谱运动目标探测方法及装置 | |
Bousi et al. | Lateral resolution improvement in Optical Coherence Tomography (OCT) images | |
CN113259558B (zh) | 一种光场相机的无损全聚焦方法及装置 | |
CN110543028B (zh) | 一种满足微光成像的实时偏振成像阵列结构 | |
CN111854956A (zh) | 一种基于微透镜阵列的多光谱成像系统及图像重建方法 | |
CN113900262A (zh) | 一种基于广义涡旋光束的超颖表面设计方法及制备方法 | |
CN108665422B (zh) | 傅里叶域内逆向感知的单帧红外非均匀性检测方法 | |
JP7373399B2 (ja) | 高精細度で電磁波を解析するための方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |