CN110567883B - 实时测量Mueller矩阵光谱的方法 - Google Patents
实时测量Mueller矩阵光谱的方法 Download PDFInfo
- Publication number
- CN110567883B CN110567883B CN201910899485.5A CN201910899485A CN110567883B CN 110567883 B CN110567883 B CN 110567883B CN 201910899485 A CN201910899485 A CN 201910899485A CN 110567883 B CN110567883 B CN 110567883B
- Authority
- CN
- China
- Prior art keywords
- fourier transform
- mueller matrix
- spectrometer
- retarder
- inverse fourier
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/21—Polarisation-affecting properties
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Software Systems (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了实时测量Mueller矩阵光谱的系统,包括沿入射光线的入射方向依次设有光源、光纤、准直系统、起偏器、第一延迟器、第二延迟器、样品台、第三延迟器、Wollaston棱镜、光谱仪组。本发明还公开了实时测量Mueller矩阵光谱的方法,对光谱仪获取的光谱强度进行傅里叶变换、通道滤波与傅里叶逆变换可复原被测样品表面出射的Stokes矢量S,对S中的四个Stokes参数进行傅里叶变换、通道滤波与傅里叶逆变换可复原出16个Mueller矩阵元素的光谱。本发明避免了通道串扰,测量准确性高,可操作性强。
Description
技术领域
本发明属于偏振光谱测量技术领域,涉及实时测量Mueller矩阵光谱的系统,还涉及实时测量Mueller矩阵光谱的方法。
背景技术
Mueller矩阵可完全描述物体对光的偏振响应属性,通过获取被测样品的Mueller矩阵光谱可得被测样品的诸多光学特征,如退偏效应、表面刻蚀的各向异性以及断面的非对称性。Mueller矩阵光谱测量技术作为一种新型的光学探测方法在集成电路制造、平板显示、生物医学诊断等领域具有广泛的应用前景。目前国际上只有少数科研机构开展相关研究工作。尤其是实时获取Mueller矩阵光谱的仪器,目前仍处于原理探索与实验验证阶段,其技术手段尚未成熟。
目前国际上常规的Mueller矩阵光谱测量方法主要有两类:一类是是按时间顺序进行的,起偏臂与检偏臂的补偿器按一定速率比5ω:mω(ω为同轴电机的基频,m为非零整数且不等于5或10)旋转产生不同的时间调制频率,通过傅里叶变换解调获取样品全部16个穆勒矩阵元素需要探测器读数 N=4m+21次以保证在一个光学周期T=π/ω内完成完整的波形分析,耗费的时间为t=N×t0(t0为探测器的读出时间),测量时间较长不适用于被测穆勒矩阵随时间快速变化的情况,并且系统中两个运动部件增加了系统误差产生的概率。另一类是则通过频率调制将被测穆勒矩阵全部16个元素被调制到37 个不同的频率通道,采用通道滤波与傅里叶变换解调可获取全部16个穆勒矩阵元素的光谱,由于系统中不含运动部件,可实现静态实时测量(测量时间由探测器读出时间决定),但是当被测穆勒矩阵具有锐利特征峰时会出现严重的通道串扰从而影响测量准确度与精度,并且复原穆勒矩阵的光谱分辨率远小于光谱仪的光谱分辨率,因此其仅适用于被测穆勒矩阵随波长缓慢变化的情况。
发明内容
本发明的目的是提供实时测量Mueller矩阵光谱的系统,可实时获取被测Mueller矩阵光谱,同时可抑制通道串扰的产生。
本发明的另一目的是提供实时测量Mueller矩阵光谱的方法。
本发明所采用的第一种技术方案是,实时测量Mueller矩阵光谱的系统,包括沿入射光线的入射方向依次设有光源、光纤、准直系统、起偏器、第一延迟器、第二延迟器、样品台、第三延迟器、Wollaston棱镜、光谱仪组。
本发明第一种技术方案的特点还在于,
准直系统包括第一透镜和第二透镜,第一透镜和第二透镜分别靠近光纤和起偏器,光谱仪组包括第一光谱仪和第二光谱仪,第一光谱仪和第二光谱仪均位于Wollaston棱镜的出射侧以测量不同出射光的光谱强度。
其构建满足右手定则的xyz坐标系,以光线的主光轴方向为z轴正向,以起偏器的透光轴方向为x轴正向,第一延迟器的快轴方向与x轴正向的夹角为45°,第二延迟器的快轴方向与x轴正向的夹角为0°,第三延迟器的快轴方向与x轴正向的夹角为22.5°,Wollaston棱镜的两光轴分别位于yz与xz 平面且均与z轴垂直。
第一延迟器与第二延迟器相同。
本发明第二种技术方案是,实时测量Mueller矩阵光谱的方法,应用如本发明第一种方案的实时测量Mueller矩阵光谱的系统进行测量,具体按照以下步骤实施:
步骤1、测量从Wollaston棱镜(9)出射两束光的光谱强度为IWP1和IWP2,并根据四个Stokes参数建立IWP1和IWP2的表达式;
步骤2、对IWP1和IWP2的表达式分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到Stokes参数;
步骤3、对四个Stokes参数分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到十六个Mueller矩阵元素的光谱,Mueller矩阵光谱的测量完成。
本发明第二种技术方案的特点还在于,
步骤1具体为,打开光源,从光纤出射的非偏振光经准直系统准直后变为平行光,平行光依次通过起偏器、第一延迟器、第二延迟器、样品台、第三延迟器,由Wollaston棱镜射出两束光,采用第一光谱仪、第二光谱仪分别同时测量从Wollaston棱镜出射两束光的光谱强度为IWP1和IWP2,其表达式分别如下,
上式中,X0、X1、X2和X3为Stokes参数Xa,a=0,1,2,3,
步骤2具体为,
由式(15)~(17)可得Stokes参数:
式(18)~(21)中,real和imag分别表示实部和虚部。
步骤3具体为,
对式(18)~(21)式做傅里叶变换并进行通道滤波截取C0、C1、C2,并对C0、C1、C2做傅里叶逆变换可得被测的十六个穆勒矩阵元素:
此时,Mueller矩阵光谱的测量完成。
本发明的有益效果是:
本发明实时测量Mueller矩阵光谱的系统与现有时间调制型Mueller矩阵光谱探测系统相比可实时获取被测Mueller矩阵光谱,且系统中不含运动部件;与现有频率调制型Mueller矩阵光谱探测系统与方法相比,光谱通道由37个降低为15个,提高了复原光谱分辨率的同时可抑制通道串扰产生的概率。本发明实时测量Mueller矩阵光谱的方法通过起偏器、第一延迟器和第二延迟器将从被测样品表面出射的Stokes矢量S变为被测样品的全部16 个Mueller矩阵元素光谱的函数,四个Stokes参数均含有五个不同的光谱通道;再通过第三延迟器和Wollaston棱镜将这四个Stokes参数加载至三个不同的光谱通道;对光谱仪获取的光谱强度进行傅里叶变换、通道滤波与傅里叶逆变换可复原S,对S中的四个Stokes参数进行傅里叶变换、通道滤波与傅里叶逆变换可复原出16个Mueller矩阵元素的光谱;测量准确性高,且操作简单清楚。
附图说明
图1为本发明实时测量Mueller矩阵光谱的系统的结构示意图;
其中,1.光源,2.光纤,3.准直系统,301.第一透镜,302.第二透镜,4. 起偏器,5.第一延迟器,6.第二延迟器,7.样品台,8.第三延迟器,9.Wollaston 棱镜,10.第一光谱仪,11.第二光谱仪。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
本发明实时测量Mueller矩阵光谱的系统,其构建满足右手定则的xyz 坐标系,如图1所示,包括沿入射光线的入射方向依次设有光源1、光纤2、准直系统3、起偏器4、第一延迟器5、第二延迟器6、样品台7、第三延迟器8、Wollaston棱镜9、光谱仪组。以光线的主光轴方向为z轴正向,以起偏器4的透光轴方向为x轴正向,第一延迟器5的快轴方向与x轴正向的夹角为45°,第二延迟器6的快轴方向与x轴正向的夹角为0°,第三延迟器8 的快轴方向与x轴正向的夹角为22.5°,Wollaston棱镜9的两光轴分别位于 yz与xz平面且均与z轴垂直。
准直系统包括第一透镜301和第二透镜302,第一透镜301和第二透镜 302分别靠近光纤2和起偏器4,光谱仪组包括第一光谱仪10和第二光谱仪11,第一光谱仪10和第二光谱仪11均位于Wollaston棱镜9的出射侧以测量不同出射光的光谱强度。第一延迟器5与第二延迟器6相同。
本发明实时测量Mueller矩阵光谱的方法,应用本发明的实时测量 Mueller矩阵光谱的系统进行测量,具体按照以下步骤实施:
步骤1、测量从Wollaston棱镜(9)出射两束光的光谱强度为IWP1和IWP2,并根据四个Stokes参数建立IWP1和IWP2的表达式:
打开光源1,从光纤2出射的非偏振光经准直系统3准直后变为平行光,平行光依次通过起偏器4、第一延迟器5、第二延迟器6、样品台7、第三延迟器8,由Wollaston棱镜9射出两束光,采用第一光谱仪10、第二光谱仪 11分别同时测量从Wollaston棱镜9出射两束光的光谱强度为IWP1和IWP2,其表达式分别如下,
上式中,X0、X1、X2和X3为Stokes参数Xa,a=0,1,2,3,
步骤2、对IWP1和IWP2的表达式分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到Stokes参数:
由式(15)~(17)可得Stokes参数:
式(18)~(21)中,real和imag分别表示实部和虚部。
步骤3、对四个Stokes参数分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到十六个Mueller矩阵元素的光谱,Mueller矩阵光谱的测量完成:
对式(18)~(21)式做傅里叶变换并进行通道滤波截取C0、C1、C2,并对C0、C1、C2做傅里叶逆变换可得被测的十六个穆勒矩阵元素:
此时,Mueller矩阵光谱的测量完成。
通过上述方式,本发明实时测量Mueller矩阵光谱的方法,对光谱仪获取的光谱强度进行傅里叶变换、通道滤波与傅里叶逆变换可复原被测样品表面出射的Stokes矢量S,对S中的四个Stokes参数进行傅里叶变换、通道滤波与傅里叶逆变换可复原出16个Mueller矩阵元素的光谱。测量准确性高,可操作性强。
Claims (4)
1.实时测量Mueller矩阵光谱的方法,应用实时测量Mueller矩阵光谱的系统进行测量,其特征在于,实时测量Mueller矩阵光谱的系统,包括沿入射光线的入射方向依次设有光源(1)、光纤(2)、准直系统(3)、起偏器(4)、第一延迟器(5)、第二延迟器(6)、样品台(7)、第三延迟器(8)、Wollaston棱镜(9)、光谱仪组;
所述准直系统包括第一透镜(301)和第二透镜(302),所述第一透镜(301)和第二透镜(302)分别靠近光纤(2)和起偏器(4),所述光谱仪组包括第一光谱仪(10)和第二光谱仪(11),所述第一光谱仪(10)和第二光谱仪(11)均位于Wollaston棱镜(9)的出射侧以测量不同出射光的光谱强度;
实时测量Mueller矩阵光谱的系统构建满足右手定则的xyz坐标系,以光线的主光轴方向为z轴正向,以所述起偏器(4)的透光轴方向为x轴正向,第一延迟器(5)的快轴方向与x轴正向的夹角为45°,第二延迟器(6)的快轴方向与x轴正向的夹角为0°,所述第一延迟器(5)与第二延迟器(6)相同,第三延迟器(8)的快轴方向与x轴正向的夹角为22.5°,Wollaston棱镜(9)的两光轴分别位于yz与xz平面且均与z轴垂直;
具体按照以下步骤实施:
步骤1、测量从Wollaston棱镜(9)出射两束光的光谱强度为IWP1和IWP2,并根据四个Stokes参数建立IWP1和IWP2的表达式;
步骤2、对IWP1和IWP2的表达式分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到Stokes参数;
步骤3、对四个Stokes参数分别做傅里叶变换、通道滤波,对滤波结果做傅里叶逆变换,得到十六个Mueller矩阵元素的光谱,Mueller矩阵光谱的测量完成。
2.根据权利要求1所述的实时测量Mueller矩阵光谱的方法,其特征在于,所述步骤1具体为,打开光源(1),从光纤(2)出射的非偏振光经准直系统(3)准直后变为平行光,平行光依次通过起偏器(4)、第一延迟器(5)、第二延迟器(6)、样品台(7)、第三延迟器(8),由Wollaston棱镜(9)射出两束光,采用第一光谱仪(10)、第二光谱仪(11)分别同时测量从Wollaston棱镜(9)出射两束光的光谱强度为IWP1和IWP2,其表达式分别如下,
上式中,X0、X1、X2和X3为Stokes参数Xa,a=0,1,2,3,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910899485.5A CN110567883B (zh) | 2019-09-23 | 2019-09-23 | 实时测量Mueller矩阵光谱的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910899485.5A CN110567883B (zh) | 2019-09-23 | 2019-09-23 | 实时测量Mueller矩阵光谱的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110567883A CN110567883A (zh) | 2019-12-13 |
CN110567883B true CN110567883B (zh) | 2022-07-01 |
Family
ID=68781933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910899485.5A Active CN110567883B (zh) | 2019-09-23 | 2019-09-23 | 实时测量Mueller矩阵光谱的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110567883B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110261317B (zh) * | 2019-06-17 | 2021-11-16 | 西安理工大学 | 一种Mueller矩阵光谱的测量系统及方法 |
CN111413282A (zh) * | 2020-04-11 | 2020-07-14 | 华中科技大学 | 一种光弹型高速穆勒矩阵椭偏仪及其原位校准与测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107036713A (zh) * | 2017-05-17 | 2017-08-11 | 西安交通大学 | 图像、高分辨率强度光谱与线偏振光谱探测系统与方法 |
CN108918425A (zh) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | 一种穆勒矩阵测量系统及方法 |
CN110261319A (zh) * | 2019-06-24 | 2019-09-20 | 西安理工大学 | 基于四次测量Mueller矩阵光谱的装置及测量方法 |
CN110261317A (zh) * | 2019-06-17 | 2019-09-20 | 西安理工大学 | 一种Mueller矩阵光谱的测量系统及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7277172B2 (en) * | 2005-06-06 | 2007-10-02 | Kla-Tencor Technologies, Corporation | Measuring overlay and profile asymmetry using symmetric and anti-symmetric scatterometry signals |
US7298480B2 (en) * | 2005-12-23 | 2007-11-20 | Ecole Polytechnique | Broadband ellipsometer / polarimeter system |
-
2019
- 2019-09-23 CN CN201910899485.5A patent/CN110567883B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107036713A (zh) * | 2017-05-17 | 2017-08-11 | 西安交通大学 | 图像、高分辨率强度光谱与线偏振光谱探测系统与方法 |
CN108918425A (zh) * | 2018-06-01 | 2018-11-30 | 中国科学院西安光学精密机械研究所 | 一种穆勒矩阵测量系统及方法 |
CN110261317A (zh) * | 2019-06-17 | 2019-09-20 | 西安理工大学 | 一种Mueller矩阵光谱的测量系统及方法 |
CN110261319A (zh) * | 2019-06-24 | 2019-09-20 | 西安理工大学 | 基于四次测量Mueller矩阵光谱的装置及测量方法 |
Non-Patent Citations (2)
Title |
---|
Channeled spectropolarimetry with increased bandwidth and aliasing reduction;Naicheng Quan et al.;《Optik》;20190330;第181卷;第562-570页 * |
基于孔径分割与视场分割的通道型成像光谱偏振技术;权乃承等;《物理学报》;20160331;第64卷(第14期);第080703-1页-080703-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110567883A (zh) | 2019-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110261317B (zh) | 一种Mueller矩阵光谱的测量系统及方法 | |
Garcia-Caurel et al. | Advanced Mueller ellipsometry instrumentation and data analysis | |
CN1811357B (zh) | 分光偏振测定法 | |
CN107367329B (zh) | 一种图像、光谱、偏振态一体化获取装置及探测方法 | |
CN102332956B (zh) | 一种宽带光源的色散补偿方法 | |
CN103954360B (zh) | 一种基于偏振阵列的光谱偏振装置及探测方法 | |
CN110579443B (zh) | 一种Mueller矩阵光谱的测量系统及其测量方法 | |
CN102589850B (zh) | 一种波片相位延迟的精密测量系统及其实现方法 | |
CN111413282A (zh) | 一种光弹型高速穆勒矩阵椭偏仪及其原位校准与测量方法 | |
CN113777049B (zh) | 一种角分辨快照椭偏仪及其测量系统与方法 | |
JPS6134442A (ja) | 試料表面ないしは試料の表面膜層の物理的特性を検査するためのエリプソメトリ測定法とその装置 | |
CN102539119A (zh) | 基于可旋转波片的Mueller矩阵测试装置和方法 | |
CN110567883B (zh) | 实时测量Mueller矩阵光谱的方法 | |
CN109990736B (zh) | 一种基于斯托克斯矢量的滚转角测量方法及装置 | |
CN110261319B (zh) | 基于四次测量Mueller矩阵光谱的装置及测量方法 | |
CN102426058A (zh) | 一种静态干涉成像偏振仪及获得目标的偏振信息的方法 | |
CN105136299A (zh) | 一种基于pem的新型红外光谱反演方法及其装置 | |
Burtscher et al. | Observing faint targets with MIDI at the VLTI: the MIDI AGN large programme experience | |
CN109612585B (zh) | 一种基于四分区相位延迟阵列的全光信息获取装置及方法 | |
CN114152578B (zh) | 基于涡旋波片的空间调制偏振检测方法 | |
CN106442336A (zh) | 一种基于光谱分析的偏振特性测试方法 | |
Zhou et al. | Adaptive linear reconstruction with regularizer for channeled spectropolarimeter | |
CN103575662B (zh) | 光学测量系统 | |
CN108956480A (zh) | 一种实时获取光谱椭偏参数的装置及方法 | |
CN112710615B (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 |