CN111256827B - 一种基于矢量光场调制与图像处理的光波偏振态检测方法 - Google Patents
一种基于矢量光场调制与图像处理的光波偏振态检测方法 Download PDFInfo
- Publication number
- CN111256827B CN111256827B CN202010033630.4A CN202010033630A CN111256827B CN 111256827 B CN111256827 B CN 111256827B CN 202010033630 A CN202010033630 A CN 202010033630A CN 111256827 B CN111256827 B CN 111256827B
- Authority
- CN
- China
- Prior art keywords
- light
- wave
- wave plate
- quarter
- 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.)
- Active
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 66
- 238000001514 detection method Methods 0.000 title claims abstract description 43
- 238000012545 processing Methods 0.000 title claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000000034 method Methods 0.000 claims description 20
- 238000012544 monitoring process Methods 0.000 claims description 4
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 238000010191 image analysis Methods 0.000 abstract 1
- 230000003287 optical effect Effects 0.000 description 16
- 239000013307 optical fiber Substances 0.000 description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 208000000453 Skin Neoplasms Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002875 fluorescence polarization Methods 0.000 description 1
- 238000003018 immunoassay Methods 0.000 description 1
- GQYHUHYESMUTHG-UHFFFAOYSA-N lithium niobate Chemical compound [Li+].[O-][Nb](=O)=O GQYHUHYESMUTHG-UHFFFAOYSA-N 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 201000000849 skin cancer Diseases 0.000 description 1
- 230000003746 surface roughness 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
- G01J4/00—Measuring polarisation of light
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于偏振检测技术领域,公开了一种基于矢量光场调制与图像处理的光波偏振态检测方法。待测光源发出的偏振态未知的光波经四分之一波片后入射至零级涡旋半波片,零级涡旋半波片将入射光场转化为矢量偏振光场,此矢量偏振光场经检偏器作用后形成亮暗呈楔形分布的光强图像。在检测时,通过前后两次调整四分之一波片的角度,使光强图像亮区分别到达两种特定的方位并采集光强图像,将两次采集的光强图像进行图像分析处理,并结合两次测量中四分之一波片的角度即可解算出待测光波的偏振态。本发明操作便捷、检测精度高,检测结果对光源的功率和波长变化不敏感。
Description
技术领域
本发明涉及偏振检测技术领域,进一步是涉及一种利用矢量偏振光场调制与数字图像处理技术实现光波偏振态检测的方法。
背景技术
偏振是光波最基本的特征之一,偏振态的检测在很多领域中都有重要的应用,如在遥感探测中,可以从目标反射光的偏振状态中得知目标的方向朝向、形状、表面粗糙度等信息;在大容量、高速率、远距离光纤通信系统中需要对光纤中的偏振效应(如偏振模色散,偏振相关调制及损耗等)做出及时的检测、抑制与补偿,以提高信息的传输质量与传输距离;在生物医学领域,偏振检测可以用于皮肤癌的诊断、荧光偏振免疫分析等;在光学检测领域,可以利用椭偏仪检测出光波经薄膜样品透射或反射前后偏振态的变化,解算出薄膜样品的厚度与折射率等光学特性。因此实现对光波偏振态快速准确的检测具有重要的应用价值与实际意义。
目前实现偏振检测的方法主要有:调制型、分光型和干涉型三类。调制型检测系统通过旋转光路中的光学元件(如偏振片,波片等)或在光路中引入调制器件(如电光调制,磁光调制等)对光强进行调制,测量调制后的光强而得到斯托克斯参量。这类检测系统结构简单,但其一般适用于测量偏振态变化较慢的光波,且光学元件的机械旋转和光源功率的波动会引入测量误差。分光型检测方法进一步可分为分振幅和分波前两种方法,分振幅法利用部分反射镜或其它分光器件将待测光分为设置有不同检偏器件的多路,利用多个光电探测器同时测量各路功率,进而得到某一时刻光波的斯托克斯参量;分波前法通常是用同尺寸的多个光阑将光波的波前按功率等分为多路,然后对每路光信号进行相互关联的调制,再同时对各路光强进行探测并解算得到斯托克斯参量。分光型检测系统采用多通道同时测量的方法,有效地提高了测量速度,适用于实时检测场合,但降低了能量利用率并牺牲部分空间分辨率,且对各通道的空间配准要求较高。干涉型检测系统通过构造干涉光路,将多路调制有偏振信息的光束汇聚在探测面使其发生干涉,通过分析干涉图像获得光波的偏振状态,该方法仅适用于单波长或窄线宽光波的偏振检测,且从干涉图像中提取偏振信息的运算量较大。
近年来广泛研究的偏振态检测方法,基本都可归结为以上三类。调制型如对比文件1中公布的(中国发明专利,201310176136.3)“光波偏振态高速静态测量装置及测量方法”,对铌酸锂晶体进行横向电光调制来实现光波偏振态的静态测量,测量过程中没有任何机械旋转,系统稳定性好,测量速度快,但其缺点是晶体自然双折射的存在会导致调制光发生畸变,因此测量过程需要采用其它方法补偿自然双折射的影响。分光型如对比文件2中公布的(中国发明专利,201410779529.8)“一种快速偏振检测仪和检测方法”,通过在检偏模块设置至少3个通道,同时对同一待测偏振光束进行分束检测,来获得待测光束的偏振状态,该方法可同时测量出待测光的全部斯托克斯参量,测量速度快,但其对各通道的空间配准要求较高。干涉型如对比文件3中公布的(中国发明专利,201110226648.7)“监测光纤偏振变化的方法与光路系统”,采用光纤耦合器构成干涉光路,通过对干涉光强进行检测来监测光纤偏振状态的变化,但该方案仅适用于单波长或窄线宽光波的偏振检测。
发明内容
本发明的目的在于针对以上不足,提供一种基于矢量光场调制与图像处理的光波偏振态检测方法。本发明操作简单便捷、检测精度高,检测结果对光源的功率和波长变化不敏感。
本发明提供的一种基于矢量光场调制与图像处理的光波偏振态检测方法,技术方案是:待测光源发出的偏振态未知的光波经四分之一波片后入射至零级涡旋半波片,零级涡旋半波片将入射光转化为矢量偏振光场,此矢量偏振光场经检偏器检偏后形成亮暗呈楔形分布的图像并被相机采集,相机的输出端连接计算机,计算机对相机采集的光强图像进行处理。
本发明的详细技术方案为:
一种基于矢量光场调制与图像处理的光波偏振态检测方法,包括如下步骤:
第一步,实时监测光强图像亮区方位角;
第二步,调整四分之一波片的角度,使得图像亮区方位角为45°或135°或两者附近,记录当前四分之一波片的角度θ1;
第三步,采集光强图像I1并计算四分之一波片之后光波P、S分量的振幅比;
3.2计算四分之一波片之后光波P、S分量的振幅比;
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2;
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.2计算四分之一波片之后光波P、S分量的相位差;
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ1,τ1)和(θ2,δ2)关系确定为:
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ1,τ1)和(θ2,δ2)计算出待测光波P、S分量的振幅比τ与相位差δ,进而得到待测光波的偏振态。
本发明具有以下效益:
(1)本发明的检测结果对光源功率和波长变化不敏感,避免了光源功率和波长波动带来的检测误差;
(2)本发明提供的检测方法,其光路紧凑、操作简便、光能利用率高;
(3)通过两次对特征位置的测量,显著降低了图像对比度测量误差对检测结果的影响,提高了检测精度。
附图说明
图1为本发明提供的一种基于矢量光场调制与图像处理的光波偏振态检测方法的光路示意图;
图2为本发明提供的一种基于矢量光场调制与图像处理的光波偏振态检测方法的具体实施流程图;
图3为调节四分之一波片使亮区中心线在135°附近时的光强图像;
图4为对图3中图像处理得到亮区中心线的方位角;
图5为调节四分之一波片使亮区中心线在0°附近时的光强图像;
图6为对图5中图像处理得到亮区中心线的方位角;
图7为实施例待测光波偏振态的椭圆表示方式;
图8为实施例待测光波偏振态的庞加莱球表示方式;
其中:1—待测光源,2—四分之一波片,3—零级涡旋半波片,4—检偏器,5—相机,6—计算机。
具体实施方式
以下将结合说明书附图对本发明做进一步说明:
图1为本发明提供的一种基于矢量光场调制与图像处理的光波偏振态检测方法的光路示意图,从左至右依次并排设置待测光源1、四分之一波片2、零级涡旋半波片3、检偏器4、相机5和计算机6。待测光源1发出偏振态未知的光波经四分之一波片2后入射至零级涡旋半波片3,零度快轴设置在0°方向的零级涡旋半波片3将入射光场转化为矢量偏振光场,此矢量偏振光场经水平方向透光的检偏器4检偏后形成亮暗呈楔形分布的图像并被相机5采集,相机的输出端连接计算机6,计算机对相机采集的光强图像进行处理。
图2为本发明提供的一种基于矢量光场调制与图像处理的光波偏振态检测方法的具体实施流程图。本发明提供的检测方法具体步骤包括:
第一步,实时监测光强图像亮区方位角;
第二步,调整四分之一波片的角度,使得图像亮区方位角为45°或135°或两者附近,记录当前四分之一波片的角度θ1;
第三步,采集光强图像I1并计算四分之一波片之后光波P、S分量的振幅比;
3.2计算四分之一波片之后光波P、S分量的振幅比;
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2;
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.2计算四分之一波片之后光波P、S分量的相位差;
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ1,τ1)和(θ2,δ2)关系确定为:
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ1,τ1)和(θ2,δ2)计算出待测光波P、S分量的振幅比τ与相位差δ,进而得到待测光波的偏振态。
采用本发明检测方法进行了仿真验证实验,结果表明本发明提供的方法可有效检测光波的偏振态。
仿真验证实例具体实施中所采用的待测光波的P、S分量振幅比和相位差分别为:τ=3,δ=120°;仿真实验过程中,在光强图像中加入高斯噪声以模拟实际采集图像中的噪声,噪声参数设置为0.01。
首先旋转四分之一波片使得入射光经过检测系统后的光强图像如图3所示(亮区中心线在135度附近),记录此时四分之一波片的方位角θ1=27.0°。
待测光波经过四分之一波片后P分量与S分量振幅比可由(式7)计算:
再次旋转四分之一波片使得入射光经过偏振态检测系统后,在相机上采集的图像如图5所示(亮区中心线在0度附近),记录此时四分之一波片的方位角θ2=51.3°。
待测光波经过四分之一波片后P分量与S分量相位差可由(式8)计算:
利用(式9)、(式10)与测量中得到的参量(θ1,τ1)和(θ2,δ2)即可计算出待测光波P、S分量的振幅比τ与相位差δ:
从图3与图5可以看出,所采集的图像中虽然包含了大量的噪声,但是利用本发明所提供的检测光波偏振态方法,检测得到的光波P、S分量的振幅比τ与相位差δ误差分别为0.0005与0.007°以内,表现出了极高的检测精度。图7为此待测光波偏振态的椭圆表示,图8为此待测光波偏振态的庞加莱球表示。
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种基于矢量光场调制与图像处理的光波偏振态检测方法,其特征在于,从左至右依次并排设置待测光源、四分之一波片、零级涡旋半波片、检偏器、相机和计算机,待测光源发出偏振态未知的光波经四分之一波片后入射至零级涡旋半波片,零级涡旋半波片将入射光转化为矢量偏振光场,此矢量偏振光场经检偏器检偏后形成亮暗呈楔形分布的图像并被相机采集,相机的输出端连接计算机,计算机对相机采集的光强图像进行处理,具体包括如下步骤:
第一步,实时监测光强图像亮区方位角;
第二步,调整四分之一波片的角度,使得图像亮区方位角为45°或135°或两者附近,记录当前四分之一波片的角度θ1;
第三步,采集光强图像I1并计算四分之一波片之后光波P、S分量的振幅比;
3.2计算四分之一波片之后光波P、S分量的振幅比;
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2;
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.2计算四分之一波片之后光波P、S分量的相位差;
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ1,τ1)和(θ2,δ2)关系确定为:
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ1,τ1)和(θ2,δ2)计算出待测光波P、S分量的振幅比τ与相位差δ,进而得到待测光波的偏振态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010033630.4A CN111256827B (zh) | 2020-01-13 | 2020-01-13 | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010033630.4A CN111256827B (zh) | 2020-01-13 | 2020-01-13 | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111256827A CN111256827A (zh) | 2020-06-09 |
CN111256827B true CN111256827B (zh) | 2022-03-11 |
Family
ID=70945271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010033630.4A Active CN111256827B (zh) | 2020-01-13 | 2020-01-13 | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111256827B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326201B (zh) * | 2020-10-13 | 2023-02-14 | 中国科学院上海光学精密机械研究所 | 四分之一波片快轴方位角和相位延迟量分布的测量装置和方法 |
CN112629662A (zh) * | 2020-12-24 | 2021-04-09 | 华中科技大学 | 一种可重构分时偏振分析系统及探测方法 |
CN113447126B (zh) * | 2021-07-13 | 2022-02-22 | 北京理工大学 | 一种空间调制偏振检测系统及设计方法 |
CN113959957A (zh) * | 2021-10-20 | 2022-01-21 | 中国人民解放军国防科技大学 | 基于角向偏振片的旋光色散曲线测量系统及方法 |
CN114660822B (zh) * | 2022-02-28 | 2024-03-15 | 中国人民解放军战略支援部队航天工程大学 | 一种基于涡旋半波片制备多种混合庞加莱球光束的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106225928A (zh) * | 2016-07-08 | 2016-12-14 | 南京航空航天大学 | 一种矢量光束偏振分布检测装置和方法 |
CN108801464A (zh) * | 2018-04-13 | 2018-11-13 | 中国人民解放军国防科技大学 | 一种计算光波偏振方向与偏振度的图像处理方法 |
CN109990736A (zh) * | 2019-03-28 | 2019-07-09 | 华中科技大学 | 一种基于斯托克斯矢量的滚转角测量方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8004675B2 (en) * | 2007-09-20 | 2011-08-23 | Boss Nova Technologies, LLC | Method and system for stokes polarization imaging |
-
2020
- 2020-01-13 CN CN202010033630.4A patent/CN111256827B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106225928A (zh) * | 2016-07-08 | 2016-12-14 | 南京航空航天大学 | 一种矢量光束偏振分布检测装置和方法 |
CN108801464A (zh) * | 2018-04-13 | 2018-11-13 | 中国人民解放军国防科技大学 | 一种计算光波偏振方向与偏振度的图像处理方法 |
CN109990736A (zh) * | 2019-03-28 | 2019-07-09 | 华中科技大学 | 一种基于斯托克斯矢量的滚转角测量方法及装置 |
Non-Patent Citations (1)
Title |
---|
Novel ellipacmetry based on the manipulation of vectorial optical field and digital image processing;Gao CHao 等;《PROCEEDINGS OF SPIE》;20191218;第1133828-1页至第1133828-10页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111256827A (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111256827B (zh) | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 | |
CN100395538C (zh) | 一种新型快速椭圆偏振光测量系统 | |
CN103776537B (zh) | 一种偏振光斯托克斯参量的测量装置及其优化方法 | |
CN103389284B (zh) | 表面等离子体共振系统和其检测方法 | |
CN113777049B (zh) | 一种角分辨快照椭偏仪及其测量系统与方法 | |
CN1144906A (zh) | 薄膜厚度和折射率的成像检测法及其设备 | |
CN106093599B (zh) | 一种光学探头与电磁场测量设备及它们的测量方法 | |
CN108548658A (zh) | 一种单层膜光学元件应力和光学损耗同时测量的方法 | |
CN101609250A (zh) | 相机摆镜摆角扫描特性测试装置 | |
CN109115690A (zh) | 实时偏振敏感的太赫兹时域椭偏仪及光学常数测量方法 | |
CN103389285B (zh) | 表面等离子体共振系统及其检测方法 | |
CN105758625B (zh) | 一种测量遥感仪器的线偏振灵敏度的装置及方法 | |
CN105928605B (zh) | 一种检测水中声场信息的方法、装置及水下声传感器 | |
Gao et al. | Spatially modulated polarimetry based on a vortex retarder and Fourier analysis | |
CN114136894B (zh) | 一种基于涡旋波片的偏振检测系统的误差校准方法及装置 | |
CN111579075B (zh) | 一种基于傅里叶分析的光波偏振态快速检测方法 | |
CN101995292B (zh) | 反射法测量有机聚合物薄膜材料的电光系数的方法及装置 | |
CN111964873B (zh) | 一种用于保偏光纤的高精度分布式消光比测量方法 | |
CN106908002B (zh) | 一种基于光谱干涉装置的测量方法 | |
CN115325933B (zh) | 一种连续扫描式聚焦激光差分干涉仪及测量方法 | |
CN104215432B (zh) | 光源偏振状态动态反馈的相位延迟器特性检测装置及方法 | |
CN116718566A (zh) | 基于量子弱测量技术的平板玻璃折射率梯度测量装置 | |
CN116735156A (zh) | 一种多自由度望远镜测试系统 | |
CN113820051B (zh) | 材料的互补干涉应力测量装置 | |
CN109781317A (zh) | 光学玻璃应力检测系统及检测方法 |
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 |