CN111256827B - 一种基于矢量光场调制与图像处理的光波偏振态检测方法 - Google Patents

一种基于矢量光场调制与图像处理的光波偏振态检测方法 Download PDF

Info

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
Application number
CN202010033630.4A
Other languages
English (en)
Other versions
CN111256827A (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 CN202010033630.4A priority Critical patent/CN111256827B/zh
Publication of CN111256827A publication Critical patent/CN111256827A/zh
Application granted granted Critical
Publication of CN111256827B publication Critical patent/CN111256827B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J4/00Measuring 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.1采集并处理光强图像I1,得到光强图像I1的亮区方位角
Figure BDA0002365239900000023
与亮暗区域对比度C1,对比度C1的定义为:
Figure BDA0002365239900000021
3.2计算四分之一波片之后光波P、S分量的振幅比;
利用图像处理得到的亮区方位角
Figure BDA0002365239900000022
与亮暗区域对比度C1计算四分之一波片之后光波P分量和S分量的振幅比τ1
Figure BDA0002365239900000031
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.1采集并处理光强图像I2得到光强图像I2亮区方位角
Figure BDA0002365239900000032
与亮暗区域对比度C2
5.2计算四分之一波片之后光波P、S分量的相位差;
利用图像处理得到的亮区方位角
Figure BDA0002365239900000033
与亮暗区域对比度C2,计算四分之一波片之后光波P分量和S分量的相位差δ2
Figure BDA0002365239900000034
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ11)和(θ22)关系确定为:
Figure BDA0002365239900000035
Figure BDA0002365239900000036
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ11)和(θ22)计算出待测光波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.1采集并处理光强图像I1,得到光强图像I1的亮区方位角
Figure BDA0002365239900000041
与亮暗区域对比度C1,对比度C1的定义为:
Figure BDA0002365239900000051
3.2计算四分之一波片之后光波P、S分量的振幅比;
利用图像处理得到的亮区方位角
Figure BDA0002365239900000052
与亮暗区域对比度C1计算四分之一波片之后光波P分量和S分量的振幅比τ1
Figure BDA0002365239900000053
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.1采集并处理光强图像I2得到光强图像I2亮区方位角
Figure BDA0002365239900000054
与亮暗区域对比度C2
5.2计算四分之一波片之后光波P、S分量的相位差;
利用图像处理得到的亮区方位角
Figure BDA0002365239900000055
与亮暗区域对比度C2,计算四分之一波片之后光波P分量和S分量的相位差δ2
Figure BDA0002365239900000056
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ11)和(θ22)关系确定为:
Figure BDA0002365239900000057
Figure BDA0002365239900000058
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ11)和(θ22)计算出待测光波P、S分量的振幅比τ与相位差δ,进而得到待测光波的偏振态。
采用本发明检测方法进行了仿真验证实验,结果表明本发明提供的方法可有效检测光波的偏振态。
仿真验证实例具体实施中所采用的待测光波的P、S分量振幅比和相位差分别为:τ=3,δ=120°;仿真实验过程中,在光强图像中加入高斯噪声以模拟实际采集图像中的噪声,噪声参数设置为0.01。
首先旋转四分之一波片使得入射光经过检测系统后的光强图像如图3所示(亮区中心线在135度附近),记录此时四分之一波片的方位角θ1=27.0°。
对图3所示的光强图像采集并进行处理即可得到亮区方位角
Figure BDA0002365239900000061
如图4所示,同时得到光强图像的亮暗对比度C1
Figure BDA0002365239900000062
待测光波经过四分之一波片后P分量与S分量振幅比可由(式7)计算:
Figure BDA0002365239900000063
再次旋转四分之一波片使得入射光经过偏振态检测系统后,在相机上采集的图像如图5所示(亮区中心线在0度附近),记录此时四分之一波片的方位角θ2=51.3°。
对图5所示的光强图像进行处理即可得到亮区方位角
Figure BDA0002365239900000064
如图6所示,同时得到光强图像的亮暗对比度C2
Figure BDA0002365239900000065
待测光波经过四分之一波片后P分量与S分量相位差可由(式8)计算:
Figure BDA0002365239900000066
利用(式9)、(式10)与测量中得到的参量(θ11)和(θ22)即可计算出待测光波P、S分量的振幅比τ与相位差δ:
Figure BDA0002365239900000071
从图3与图5可以看出,所采集的图像中虽然包含了大量的噪声,但是利用本发明所提供的检测光波偏振态方法,检测得到的光波P、S分量的振幅比τ与相位差δ误差分别为0.0005与0.007°以内,表现出了极高的检测精度。图7为此待测光波偏振态的椭圆表示,图8为此待测光波偏振态的庞加莱球表示。
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (1)

1.一种基于矢量光场调制与图像处理的光波偏振态检测方法,其特征在于,从左至右依次并排设置待测光源、四分之一波片、零级涡旋半波片、检偏器、相机和计算机,待测光源发出偏振态未知的光波经四分之一波片后入射至零级涡旋半波片,零级涡旋半波片将入射光转化为矢量偏振光场,此矢量偏振光场经检偏器检偏后形成亮暗呈楔形分布的图像并被相机采集,相机的输出端连接计算机,计算机对相机采集的光强图像进行处理,具体包括如下步骤:
第一步,实时监测光强图像亮区方位角;
第二步,调整四分之一波片的角度,使得图像亮区方位角为45°或135°或两者附近,记录当前四分之一波片的角度θ1
第三步,采集光强图像I1并计算四分之一波片之后光波P、S分量的振幅比;
3.1采集并处理光强图像I1,得到光强图像I1的亮区方位角
Figure FDA0003332654840000011
与亮暗区域对比度C1,对比度C1的定义为:
Figure FDA0003332654840000012
3.2计算四分之一波片之后光波P、S分量的振幅比;
利用图像处理得到的亮区方位角
Figure FDA0003332654840000013
与亮暗区域对比度C1计算四分之一波片之后光波P分量和S分量的振幅比τ1
Figure FDA0003332654840000014
第四步,调整四分之一波片的角度,使得图像亮区方位角为0°或90°或两者附近,记录当前四分之一波片的角度θ2
第五步,采集光强I2并计算四分之一波片之后光波P、S分量的相位差;
5.1采集并处理光强图像I2得到光强图像I2亮区方位角
Figure FDA0003332654840000015
与亮暗区域对比度C2
5.2计算四分之一波片之后光波P、S分量的相位差;
利用图像处理得到的亮区方位角
Figure FDA0003332654840000016
与亮暗区域对比度C2,计算四分之一波片之后光波P分量和S分量的相位差δ2
Figure FDA0003332654840000021
第六步,计算待测光波的偏振参数,得到光波的偏振态
待测光波的P、S分量的振幅比为τ,相位差为δ,则(τ,δ)与测量中得到的参量(θ11)和(θ22)关系确定为:
Figure FDA0003332654840000022
Figure FDA0003332654840000023
其中,(式4)中|()|表示求括号内复数的模;(式5)中arg()表示求括号内复数的辐角;i表示虚数单位;利用(式4)、(式5)与测量中得到的参量(θ11)和(θ22)计算出待测光波P、S分量的振幅比τ与相位差δ,进而得到待测光波的偏振态。
CN202010033630.4A 2020-01-13 2020-01-13 一种基于矢量光场调制与图像处理的光波偏振态检测方法 Active CN111256827B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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