CN113447126B - 一种空间调制偏振检测系统及设计方法 - Google Patents
一种空间调制偏振检测系统及设计方法 Download PDFInfo
- Publication number
- CN113447126B CN113447126B CN202110790721.7A CN202110790721A CN113447126B CN 113447126 B CN113447126 B CN 113447126B CN 202110790721 A CN202110790721 A CN 202110790721A CN 113447126 B CN113447126 B CN 113447126B
- Authority
- CN
- China
- Prior art keywords
- polarization
- polarization state
- wave plate
- quarter
- fast axis
- 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 169
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001514 detection method Methods 0.000 title claims abstract description 35
- 238000013461 design Methods 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 11
- 230000003287 optical effect Effects 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 abstract description 13
- 239000011159 matrix material Substances 0.000 description 13
- 238000005259 measurement Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011478 gradient descent method Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000028161 membrane depolarization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000017105 transposition 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
技术领域
本发明属于偏振测量的技术领域,具体涉及一种空间调制偏振检测系统及设计方法。
背景技术
偏振是光的重要属性,偏振光对光学元件、材料、生物组织等样品中微结构特性十分敏感。光与样品相互作用会发生折射、反射、散射等过程进而改变入射光的偏振态。样品对光的偏振态的改变能力用穆勒Mueller矩阵表示,Mueller矩阵包含了样品的全部偏振信息,可进一步分解为退偏、相位延迟、二向色性、快轴方向角、旋光等与样品微结构密切相关的、具有实际物理意义的、可量化的偏振参数,可用于获取样品的偏振特性和结构参数。偏振测量作为光和样品偏振特性分析的重要工具,已经被广泛的应用于生物医学、量子通信、激光雷达等领域。
如图1所示,目前现有的空间调制偏振检测系统包括光源101、偏振态产生器102、空间光滤波器103、准直透镜104、相位延迟器105、偏振态分析器106、聚焦透镜107、探测器108及数据处理模块109;待测光束依次透过偏振态产生器102、空间光滤波器103、准直透镜104、空间变化的相位延迟器105及偏振态分析器106,之后光束被聚焦透镜107成像到探测器108;探测器108接收一幅待测光的强度分布,将获取的强度图发送到数据处理模块109进行数据处理,进一步给出待测光的偏振信息。在该系统中,偏振态分析器快轴方位角设置位置不相同时,所检测的结果也存在不同,因此偏振态分析器快轴方位角的朝向,决定了系统测量的精准度。
专利文件CN110806266A公开了一种偏振态检测系统中偏振态分析器的选择方法,此方法通过计算偏振几率密度的偏导数,判断偏导数是否有唯一零点来确定合适的偏振态分析器,但是该方法存在如下问题:
第一、当多次旋转偏振态分析器的波片得到多个波片快轴方位角,由此得到的多个偏振几率密度函数的偏导数同时有唯一零点时,此时无法判断出哪个唯一零点所对应波片快轴方位角最优,因此上述方法在给定几个偏振态分析器时,仅能选出合适的偏振态分析器,不能够遍历所有偏振态分析模式找出最优的偏振态分析器。
第二、设定约束条件为:选择偏振态分析器的波片可变快轴方位角,其偏振片透光轴方向固定为水平方向;因此该技术中仅考虑偏振态分析器波片的不同快轴方位角对偏振检测的影响,没有考虑偏振片不同旋转角度的影响,通过此方法得到的偏振态分析模式仅得到偏振态分析器中合适的四分之一波片快轴方位角,不能够获得合适的偏振片角度。
论文文件“Spatially modulated polarimetry based on a vortex retarderand Fourier analysis”公开了一种使用涡旋半波片作为空间变化调制器的空间调制偏振计,该方法光路结构简单通过单次拍摄获取待测光偏振态,系统稳定性、鲁棒性高。然而,该方法仅使用涡旋半波片不能够实现对圆偏振分量的空间调制,因此该方法不能直接检测斯托克斯的s3参量,仅能够实现完全偏振光的偏振态检测不能够实现部分偏振光的偏振态检测。因此,需要对空间调制器涡旋半波片进行设计优化,以实现准确测量全部归一化的斯托克斯参量。
发明内容
鉴于此,本发明提供一种空间调制偏振检测系统及设计方法,利用该系统能够实现对待测样品的偏振特性进行精准的测量。
实现本发明的技术方案如下:
一种空间调制偏振检测系统,用于对待测样品的偏振特性进行检测,包括光源、偏振态产生器、空间光滤波器、准直透镜、相位延迟器、偏振态分析器、聚焦透镜、探测器及数据处理模块;该系统还包括四分之一波片,所述四分之一波片设置于相位延迟器与准直透镜之间;所述四分之一波片的快轴方位角为所述由所述偏振态系数对比度最大来确定。
进一步地,本发明所述偏振态分析器中偏振片透光轴和偏振波片的快轴方位角为水平方向。
一种空间调制偏振检测系统的设计方法,用于对待测样品的偏振特性进行检测,具体过程为:
在空间调制偏振检测系统中相位延迟器与准直透镜之间设置四分之一波片;
多次调整所述四分之一波片的快轴方位角θ,计算每一次调整所对应的偏振态的偏振系数ν1、ν2、ν3;
进一步地,本发明所述相位延迟器采用涡旋半波片。
进一步地,本发明将所述偏振态分析器中偏振片透光轴和偏振波片的快轴方位角固定为水平方向。
进一步地,本发明所述偏振态系数ν1、ν2、ν3的表示形式如下:
所述最优偏振态分析模型为:
其中,νimax表示模型参数νi(i=1,2,3)中的较大值,mean(νimax)表示νimax的平均值,νimin表示模型参数νi中较小的值,mean(νimin)表示νimax平均值,两平均值的差值表示模型参数νi的对比度。
进一步地,本发明所述模型参数νi中的较大值是指模型参数νi所包含的n1×n2个元素,元素值较大的(n1×n2)/2个元素,模型参数νi中的较小值是指模型参数νi所包含的n1×n2个元素,元素值较大的(n1×n2)/2个元素。
设定θ的初始值;
当所计算的损失函数大于零时,增大所述θ,反之,减小所述θ;
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为目前现有的空间调制偏振检测系统的示意图;
图2为本发明设计的空间调制偏振检测系统的示意图;
其中,201-光源,202-偏振态产生器,203-空间光滤波器,204-准直透镜,205-四分之一波片,206-涡旋半波片,207-偏振态分析器,208-聚焦透镜,209-探测器,210-数据处理模块。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合;并且,基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本公开保护的范围。
需要说明的是,下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见,本文中所描述的方面可体现于广泛多种形式中,且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本公开,所属领域的技术人员应了解,本文中所描述的一个方面可与任何其它方面独立地实施,且可以各种方式组合这些方面中的两者或两者以上。举例来说,可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外,可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
本发明的设计思想为:当待测光透过系统后在探测面强度分布差异越大,对比度就越大,在利用拟合算法解算偏振态收敛性越好,解得的偏振态更为准确,因此本申请在系统中设置四分之一波片,待测光束依次透过四分之一波片、涡旋半波片、聚焦透镜、偏振态分析器,最终成像到探测器。探测器接收一幅待测光的强度分布,将获取的强度图发送到数据处理模块进行数据处理,进一步解算出待测光的偏振信息。探测面强度分布同四分之一波片、涡旋半波片及偏振态分析器的穆勒矩阵、入射光偏振态之间存在关系,涡旋半波片及偏振态分析器相对固定,多次调整四分之一波片快轴方位角,当在某一方位角下对比度达到最大,表明在该方位角下解得的偏振态的准确度越高。
基于上述设计思想,本申请实施例一种空间调制偏振检测系统,用于对待测样品的偏振特性进行检测,包括光源、偏振态产生器、空间光滤波器、准直透镜、相位延迟器、偏振态分析器、聚焦透镜、探测器及数据处理模块;该系统还包括四分之一波片,所述四分之一波片设置于相位延迟器与准直透镜之间;所述四分之一波片的快轴方位角为所述由所述偏振态系数对比度最大来确定。
本发明一实施例中,所述偏振态分析器中偏振片透光轴和偏振波片的快轴方位角为水平方向。
基于上述设计思想,本申请实施例一种空间调制偏振检测系统的设计方法,用于对待测样品的偏振特性进行检测,具体过程为:
如图2所示,在空间调制偏振检测系统中相位延迟器与准直透镜之间设置四分之一波片;
将所述偏振态分析器中偏振片透光轴和偏振波片的快轴方位角固定为水平方向;
多次调整所述四分之一波片的快轴方位角θ,计算每一次调整所对应的偏振态的偏振系数ν1、ν2、ν3;
本实施例在系统中设置四分之一波片,由于不同快轴方位角的四分之一波片,会产生不同模型参数,模型参数空间变化特性直接影响待测偏振态的检测结果,若模型参数无空间变化特性或是空间变化特性很小(对比度较小),那不同偏振态对应的强度分布差异不明显,在使用拟合算法解算待测偏振态的过程中,会导致偏振态解算不收敛,因此偏振态解算误差较大。本申请实施例中选取所述四分之一波片的最优快轴方位角是使强度分布模型中的模型参数(即偏振态系数)有最大的对比度,以使得拟合解算偏振态过程中有全局最优解,来提高偏振态的解算精度。利用本申请实施例的设计方法所设计的系统,能够精确测量偏振光的偏振态。
本申请一实施例中,相位延迟器采用涡旋半波片。
本实施例中将四分之一波片同涡旋半波片相结合,通过使用模型参数对比度最大化的优化算法,选出最优的四分之一波片的快轴方位角。该发明解决了原空间调制偏振计不能获取空间调制圆偏振分量的问题,实现了对全偏振分量的空间调制,能够直接解算出归一化的斯托克参量。对于待测光偏振态的检测,本方法不需要旋转复杂光学器件,能够实现高精度、实时动态测量。
本申请一实施例中,在本系统中,所述偏振态系数ν1、ν2、ν3的表示形式如下:
本实施例中,偏振态系数ν1、ν2、ν3表达式的推导过程为:
步骤一,获取空间变化涡旋半波片的穆勒矩阵值M;
偏振分析器的偏振片透光轴方向和波片的快轴方位角固定为水平方向,偏振片和波片的穆勒矩阵可表示为:
步骤二,建立探测面强度分布模型,利用穆勒-斯托克斯理论关系,将上述穆勒矩阵与待测光斯托克斯矢量相乘,即能得到偏振态分析器出射光的偏振态,取出射光偏振态的第一个斯托克斯分量来表示所述探测面的强度值I表示为:
其中,所指强度分布是相对强度,使用探测器获取图像的灰度值表示。考虑到器件加工和对准误差以及系统微弱变动,为弥补系统误差,上述空间变化涡旋半波片的穆勒矩阵值可以实验测得,使用测得的实验值代替理论值带入到上述强度分布模型,以实现器件加工和对准误差的标定。
步骤二中的公式可如下表示:
其中,ν1,ν2,ν3可表示为:
模型参数ν1、ν2、ν3是偏振态(s1,s2,s3)的系数,其与所选空间变化调制器的穆勒矩阵和偏振态分析器的检偏方式相关的物理量,通过上式可以将入射光偏振态映射到一幅强度分布,不同偏振态对应的强度分布伴随不同的空间变化特性。获取探测面强度图后,进一步通过拟合算法可以解算出入射光的偏振态,偏振态的解算精度取决于理论模型的模型参数ν1、ν2、ν3是否具有强的空间变化属性,以实现对全偏振分量(s1,s2,s3)的空间调制。ν1、ν2、ν3是具有强的空间变化说明模型参数的对比度较大。
所述最优偏振态分析模型为:
其中,νimax表示模型参数νi(i=1,2,3)中的较大值,mean(νimax)表示νimax的平均值,νimin表示模型参数νi中较小的值,mean(νimin)表示νimax平均值,两平均值的差值表示模型参数νi的对比度。Argmax表示取最大值函数,表示当四分之一波片的快轴方位角θ变化时取对比度的最大值。通过上式可以选择出最优的θ用表示,也即选出最优的偏振态分析模式。
本实施例中,将所述偏振系数对比度确定为:偏振系数矩阵中的较大值的均值减较小值的均值,两者之差对应模型参数对比度,由于所述偏振系数有3个,即ν1、ν2、ν3,因此针对每一偏振系数计算差值,三个偏振系数差值之和越大,则认为对比度就越大,通过该方式所计算出的对比度,可以很好的反应探测面光强差异,对比度越大,差异越大,偏振态的收敛性就越好,检测结果就越准确。
本是申请又一实施例中,模型参数νi中的较大值是指模型参数νi所包含的n1×n2个元素,元素值较大的(n1×n2)/2个元素,模型参数νi中的较小值是指模型参数νi所包含的n1×n2个元素,元素值较大的(n1×n2)/2个元素。
例如,上模型参数较大和较小灰度值确定方式为:模型参数νi,其本质是一幅强度图(假设该强度图的像素是40*40),每一个像素有其对应的灰度值,模型参数较大灰度值是1600个像素点中较大的800个像素值,模型参数较小灰度值是1600个像素点中较小的800个像素值,较大、较小像素灰度值的划分可以通过MATLAB中排序算法来实现。
设定θ的初始值;
当所计算的损失函数大于零时,增大所述θ,反之,减小所述θ;
具体实施时:使用梯度下降法快速搜索最优的梯度下降法能够沿着模型参数对比度梯度上升的方向快速寻找出模型参数对比度最大时对应的首先给θ赋值初始值,本发明给定的初始值是θ=0;进一步的,设定损失函数的改变量,上式所得结果大于零,模型参数对比度与θ成正比关系,欲最大化模型参数对比度,在迭代过中应该增大θ;反之,应该减小θ;进一步每次更新迭代的θ,找到模型参数的最大值,更新迭代如下式所示:
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
2.根据权利要求1所述空间调制偏振检测系统,其特征在于,所述偏振态分析器中偏振片透光轴和偏振波片的快轴方位角为水平方向。
4.根据权利要求3所述空间调制偏振检测系统的设计方法,其特征在于,所述相位延迟器采用涡旋半波片。
7.根据权利要求6所述空间调制偏振检测系统的设计方法,其特征在于,所述模型参数νi中的较大值是指模型参数νi所包含的n1×n2个元素,元素值较大的(n1×n2)/2个元素,模型参数νi中的较小值是指模型参数νi所包含的n1×n2个元素,元素值较小的(n1×n2)/2个元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110790721.7A CN113447126B (zh) | 2021-07-13 | 2021-07-13 | 一种空间调制偏振检测系统及设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110790721.7A CN113447126B (zh) | 2021-07-13 | 2021-07-13 | 一种空间调制偏振检测系统及设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113447126A CN113447126A (zh) | 2021-09-28 |
CN113447126B true CN113447126B (zh) | 2022-02-22 |
Family
ID=77816064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110790721.7A Active CN113447126B (zh) | 2021-07-13 | 2021-07-13 | 一种空间调制偏振检测系统及设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113447126B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114152578B (zh) * | 2021-11-26 | 2023-10-20 | 北京理工大学 | 基于涡旋波片的空间调制偏振检测方法 |
CN114200686B (zh) * | 2021-12-23 | 2023-05-30 | 西安交通大学 | 一种偏振光散射光谱成像内窥装置与探测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006153503A (ja) * | 2004-11-25 | 2006-06-15 | Osaka Industrial Promotion Organization | 複屈折測定装置および応力分布測定装置 |
CN105203209A (zh) * | 2015-08-31 | 2015-12-30 | 天津大学 | 基于光强测量积分时间优化的Stokes矢量测量系统及方法 |
CN110806266A (zh) * | 2019-11-11 | 2020-02-18 | 北京理工大学 | 一种偏振态检测系统中偏振态分析器的选择方法 |
CN111256827A (zh) * | 2020-01-13 | 2020-06-09 | 中国人民解放军国防科技大学 | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 |
-
2021
- 2021-07-13 CN CN202110790721.7A patent/CN113447126B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006153503A (ja) * | 2004-11-25 | 2006-06-15 | Osaka Industrial Promotion Organization | 複屈折測定装置および応力分布測定装置 |
CN105203209A (zh) * | 2015-08-31 | 2015-12-30 | 天津大学 | 基于光强测量积分时间优化的Stokes矢量测量系统及方法 |
CN110806266A (zh) * | 2019-11-11 | 2020-02-18 | 北京理工大学 | 一种偏振态检测系统中偏振态分析器的选择方法 |
CN111256827A (zh) * | 2020-01-13 | 2020-06-09 | 中国人民解放军国防科技大学 | 一种基于矢量光场调制与图像处理的光波偏振态检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113447126A (zh) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113447126B (zh) | 一种空间调制偏振检测系统及设计方法 | |
US10365163B2 (en) | Optical critical dimension metrology | |
US8009292B2 (en) | Single polarizer focused-beam ellipsometer | |
CN108918425A (zh) | 一种穆勒矩阵测量系统及方法 | |
CN110806266A (zh) | 一种偏振态检测系统中偏振态分析器的选择方法 | |
CN108871579A (zh) | 一种偏振成像系统的标定方法 | |
JPS6134442A (ja) | 試料表面ないしは試料の表面膜層の物理的特性を検査するためのエリプソメトリ測定法とその装置 | |
CN105628343A (zh) | 一种波片检测装置及方法 | |
KR20190118603A (ko) | 높은 공간 해상도의 일립소메트리에서 사용하기 위한 시스템 및 방법 | |
CN113538381B (zh) | 一种弱光场样品Mueller矩阵快速检测方法及检测系统 | |
CN107490547B (zh) | 一种Mueller矩阵型椭偏仪椭偏参数测量的优化方法 | |
CN109781317B (zh) | 光学玻璃应力检测系统及检测方法 | |
CN101520323B (zh) | 傅立叶光谱仪中平面动镜倾斜角的大范围测角方法 | |
CN105758625B (zh) | 一种测量遥感仪器的线偏振灵敏度的装置及方法 | |
CN109580182B (zh) | 基于布儒斯特定律的曲面光学元件折射率测量方法和装置 | |
CN114136894B (zh) | 一种基于涡旋波片的偏振检测系统的误差校准方法及装置 | |
CN111457911A (zh) | 一种基于偏振二维残差信息的仿生偏振罗盘标定方法 | |
CN116802465A (zh) | 偏振测量法 | |
CN114264632A (zh) | 一种角分辨式散射仪中物镜偏振效应的原位校准方法 | |
KR101036455B1 (ko) | 하프 미러를 이용한 타원계측기 | |
CN111307722B (zh) | 一种多功能Stokes-Mueller矩阵成像偏振仪的标定方法 | |
CN114152578B (zh) | 基于涡旋波片的空间调制偏振检测方法 | |
CN116754492A (zh) | 一种成像椭偏仪测量系统及系统参数校准方法 | |
KR102176199B1 (ko) | 타원 계측기 | |
WO2020200330A2 (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 |