CN106404675A - 一种测量粗糙表面面外偏振双向反射函数的方法 - Google Patents

一种测量粗糙表面面外偏振双向反射函数的方法 Download PDF

Info

Publication number
CN106404675A
CN106404675A CN201510478598.XA CN201510478598A CN106404675A CN 106404675 A CN106404675 A CN 106404675A CN 201510478598 A CN201510478598 A CN 201510478598A CN 106404675 A CN106404675 A CN 106404675A
Authority
CN
China
Prior art keywords
light
polarization
angle
detector
detected materials
Prior art date
Application number
CN201510478598.XA
Other languages
English (en)
Inventor
潘佳惠
陈钱
吕芳
俞晓东
Original Assignee
南京理工大学
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 南京理工大学 filed Critical 南京理工大学
Priority to CN201510478598.XA priority Critical patent/CN106404675A/zh
Publication of CN106404675A publication Critical patent/CN106404675A/zh

Links

Abstract

本发明提供一种测量粗糙表面面外偏振双向反射函数的方法。使光源产生的光线经偏振态发生器到达待测材料,待测材料反射的光线经偏振态分析器到达探测器;将待测材料围绕一定光路旋转一定角度,且保证光路在待测材料上的光斑位置不发生变化;旋转偏振态发生器或偏振态分析器,使偏振态发生器穆勒矩阵或偏振态分析器的穆勒矩阵发生变化;探测器采集一次出射光数据并存储,直至探测器至少采集十六次出射光数据;根据探测器采集并存储的出射光数据获得面外偏振双向反射函数。本发明方法原理简单,易实施,只需改变辅助角即可改变待测材料相对于光源或者探测器的天顶角或方位角,同时能有效测量方位角改变对偏振双向反射函数的影响。

Description

一种测量粗糙表面面外偏振双向反射函数的方法
技术领域
[0001] 本发明属于光学技术领域,具体涉及一种测量粗糙表面面外偏振双向反射函数的 方法。
背景技术
[0002] 入射光与界面发生相互作用产生的偏振效应是独特且有辨识性的,它能有效改善 人造物体的检测效果,同时对用于探测和识别的光谱算法的改进也有所帮助。双向反射分 布函数能有效的描述材料表面的散射光特性,而偏振双向反射分布函数是标量双向反射分 布函数的一般化形式,能够模拟表面散射光学中的偏振效应。
发明内容
[0003] 本发明的目的在于克服现有技术的不足,提出一种测量粗糙表面面外偏振双向反 射函数的方法,该方法原理简单,易实施,只需改变辅助角即可改变待测材料相对于光源或 者探测器的天顶角或方位角,同时能有效测量方位角改变对偏振双向反射函数的影响,利 于后续的科研。
[0004] 为了解决上述技术问题,本发明提供一种测量粗糙表面面外偏振双向反射函数的 方法,包括以下步骤:
[0005] 步骤1,构造光路,使光源产生的光线经偏振态发生器到达待测材料,待测材料反 射的光线经偏振态分析器到达探测器;
[0006] 步骤2,将待测材料围绕一定光路旋转一定角度,且保证在旋转过程中,光路在待 测材料上的光斑位置不发生变化;
[0007] 步骤3,旋转偏振态发生器或偏振态分析器,使偏振态发生器穆勒矩阵MP或偏振态 分析器的穆勒矩阵M A发生变化;
[0008] 步骤4,探测器采集一次出射光数据并存储,然后重复步骤3,直至探测器至少采 集十六次出射光数据;
[0009] 步骤5,根据探测器采集并存储的出射光数据获得面外偏振双向反射函数。
[0010] 进一步,在步骤2中,所述的待测材料围绕入射光光路旋转,此时旋转角为微面元 法线与宏观面法线在入射光方向的夹角^,相对于待测材料宏观面,入射光的天顶角 不变,出射光的天顶角Θ/变化。
[0011] 进一步,在步骤2中,所述的待测材料围绕出射光光路旋转,此时旋转角为微面元 法线与宏观面法线在入射光方向的夹角rU,相对于待测材料宏观面,出射光的天顶角 不变,入射光的天顶角变化。
[0012] 进一步,在步骤5中,根据光路传输方程S°ut= M AFsMFStn,直接联立出射光的16个 方程组解得面外偏振双向反射函数F s。
[0013] 进一步,在步骤5中,根据如下所示近似傅里叶级数的光强I (t)关系式,求解傅里 叶系数α2η和β 2n,获得用傅里叶系数表示的面外偏振双向反射函数Fs,
[0014] ~ ^')||
[0015] 式中,n是探测器采集出射光数据的次数,I。为辐照度对时间的平均值,a 2n和β 2n 为傅里叶系数,C为偏振态发生器或偏振态分析器中两块波片共有的单位角度,Φ2η为相位 修正角度。
[0016] 上本发明与现有技术相比,其显著优点在于,(1)原理简单,只需改变辅助角即可 改变待测材料相对于光源或者探测器的天顶角或方位角,(2)容易实施并且实施条件简单, 所需测量设备不再局限于圆弧导轨。
附图说明
[0017] 图1是待测材料粗糙表面的宏观示意图。
[0018] 图2是待测材料粗糙表面的微面元示意图。
具体实施方式
[0019] 容易理解,依据本发明的技术方案,在不变更本发明的实质精神的情况下,本领域 的一般技术人员可以想象出本发明测量粗糙表面面外偏振双向反射函数的方法的多种实 施方式。因此,以下具体实施方式和附图仅是对本发明的技术方案的示例性说明,而不应当 视为本发明的全部或者视为对本发明技术方案的限制或限定。
[0020] 基于微面元模型的偏振双向反射函数在学术界得到广泛的认可,其假设一般的粗 糙材料表面是由一系列的微面元构成,而每一个微面元认为是遵循菲尼尔方程的镜面反射 面,此外微面元的整体分布遵循一定的统计分布函数。4*4的偏振双向反射函数F s具体定 义式如下:
[0021]
Figure CN106404675AD00041
[0025] 其中,如图1所示,喝分别代表入射光的天顶角和方位角,L分别代表反射 光的天顶角和方位角,方位角~%如图2所示,Θ是微面元法向量和宏观面法向量 之间的夹角,β是微面元的天顶角。就本质上而言,穆勒矩阵Μ是考虑偏振和微面元效应 后的反射率。Ρ(σ,Θ)是高斯分布函数,σ作为可变的宽度分布因子,实质的物理意义与 表面粗糙度相关。前述方程(1)中的穆勒矩阵Μ-般由琼斯矩阵Τ推导而得,其简单推导 如下所述:
Figure CN106404675AD00051
[0029] 其中,參_分别代表入射光和反射光中的s光的电场分量,而考分别代表 入射光和反射光中的P光的电场分量。S光和p光的菲涅尔反射系数rs、rp将入射电场与 反射电场两者关联。辅助角1、ru分别是宏观面法线和微面元法线之间沿入射方向、出射 方向的夹角。一般而言,1、L两者不等,但是两者都是由决定的唯一值。最后 通过穆勒矩阵Μ与琼斯矩阵Τ的关系式即可得到粗糙表面的偏振双向反射函数。
[0030] 为测得面外偏振双向反射函数,传统的测量方法是将样本固定于图1的χ-y平面, 调节光源和探测器的位置实现方位角的变化。这样的方法易于理解,实际操作却比较复杂。 为了保证在光源和探测器的位置变化时探测器采集到的待测材料无位移,就需要价格昂贵 的圆弧导轨,并且待测材料必须保证固定于圆弧导轨的圆心。
[0031] 本发明所述测量粗糙表面面外偏振双向反射函数的方法,具体步骤如下:
[0032] 步骤1,构造光路,使光源产生的光线经偏振态发生器到达待测材料,待测材料反 射的光线经偏振态分析器到达探测器;
[0033] 步骤2,将待测材料围绕入射光路旋转且在旋转过程中入射光路在待测材料上的 光斑位置不发生变化,此时旋转角为微面元法线与宏观面法线在入射光路方向的夹角η ρ 相对于待测材料宏观面而言,入射光路的天顶角不变,出射光路的天顶角θ /变化。在 前述方程(6)中,入射光路的天顶角和微面元的天顶角β不变,已知辅助角η i即可推 导出相应的出射光路的天顶角θρ再代入前述方程(3)求得方位角'这样就间接改变入 射光路方位角:¾与出射光路方位角象<的夹角&即_ ~爲丨。
[0034] 在上述步骤2中,所述的待测材料还可以围绕出射光路旋转且在旋转过程中出射 光路在待测材料上的光斑位置不发生变化,此时旋转角为微面元法线与宏观面法线在入射 光路方向的夹角相对于待测材料宏观面而言,出射光路的天顶角θι^不变,入射光路的 天顶角变化。在前述方程(7)中,出射光路的天顶角Θ 微面元的天顶角β不变,已 知辅助角L即可推导出相应的入射光路的天顶角Θ i,再代入前述方程(3)求得方位角 这样就间接改变入射光路方位角崎与出射光路方位角瓢的夹角秦即渗《料-N。
[0035] 步骤3,旋转偏振态发生器或偏振态分析器,使偏振态发生器穆勒矩阵MP或偏振态 分析器的穆勒矩阵M A发生变化;偏振态发生器和偏振态分析器都是由一块偏振片和一块波 片构成,对于入射光先经过偏振态发生器的偏振片再经过波片,而出射光先经过偏振态分 析器的波片再经过偏振片到达探测器。
[0036] 步骤4,探测器采集一次出射光数据并存储,重复步骤3。由于偏振双向反射函数& 是4*4的矩阵,至少共旋转偏振态发生器或偏振态分析器十六次,探测器至少再次十六次 下分别采集十六次出射光数据才能求得该函数;
[0037] 步骤5,将探测器采集并存储的所有出射光数据根据光路传输方程S°ut = MAFsMFStn,直接联立16个出射光方程组解得面外偏振双向反射函数F s。重复步骤2,进行下 一方位角的数据测量。
[0038] 在步骤5中,所述的数据处理还可以是傅里叶分析。由于探测器获得的数据是出 射光的光强信息,所以将光路传输方程:S° ut= MAFsMFStn按第一行展开,得到近似傅里叶级 数的光强关系式:
[0039]
Figure CN106404675AD00061
[0040] 式中,η是探测器采集出射光数据的次数,I。为辐照度对时间的平均值,α 2n和β 2n 为傅里叶系数,C为偏振态发生器或偏振态分析器中两块波片共有的单位角度,Φ2η为相位 修正角度。

Claims (5)

1. 一种测量粗糙表面面外偏振双向反射函数的方法,其特征在于,包括以下步骤: 步骤1,构造光路,使光源产生的光线经偏振态发生器到达待测材料,待测材料反射的 光线经偏振态分析器到达探测器; 步骤2,将待测材料围绕一定光路旋转一定角度,且保证在旋转过程中,光路在待测材 料上的光斑位置不发生变化; 步骤3,旋转偏振态发生器或偏振态分析器,使偏振态发生器穆勒矩阵MP或偏振态分析 器的穆勒矩阵MA发生变化; 步骤4,探测器采集一次出射光数据并存储,然后重复步骤3,直至探测器至少采集 十六次出射光数据; 步骤5,根据探测器采集并存储的出射光数据获得面外偏振双向反射函数。
2. 根据权利1所述面外偏振双向反射函数的方法,其特征在于,步骤2中,所述的待 测材料围绕入射光光路旋转,此时旋转角为微面元法线与宏观面法线在入射光方向的夹角 n i,相对于待测材料宏观面,入射光的天顶角不变,出射光的天顶角Θ /变化。
3. 根据权利1所述面外偏振双向反射函数的方法,其特征在于,步骤2中,所述的待 测材料围绕出射光光路旋转,此时旋转角为微面元法线与宏观面法线在入射光方向的夹角 ru,相对于待测材料宏观面,出射光的天顶角θι^不变,入射光的天顶角Θ i变化。
4. 根据权利1至3所述面外偏振双向反射函数的方法,其特征在于,步骤5中,根据光 路传输方程S°ut= MAFsMPStn,直接联立出射光的16个方程组解得面外偏振双向反射函数F s。
5. 根据权利1至3所述面外偏振双向反射函数的方法,其特征在于,步骤5中,根据如 下所示近似傅里叶级数的光强I(t)关系式,求解傅里叶系数α 2η和β 2n,获得用傅里叶系 数表示的面外偏振双向反射函数Fs,
Figure CN106404675AC00021
式中,η是探测器采集出射光数据的次数,I。为辐照度对时间的平均值,α 2n和β 2n为 傅里叶系数,C为偏振态发生器或偏振态分析器中两块波片共有的单位角度,Φ2η为相位修 正角度。
CN201510478598.XA 2015-08-03 2015-08-03 一种测量粗糙表面面外偏振双向反射函数的方法 CN106404675A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510478598.XA CN106404675A (zh) 2015-08-03 2015-08-03 一种测量粗糙表面面外偏振双向反射函数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510478598.XA CN106404675A (zh) 2015-08-03 2015-08-03 一种测量粗糙表面面外偏振双向反射函数的方法

Publications (1)

Publication Number Publication Date
CN106404675A true CN106404675A (zh) 2017-02-15

Family

ID=58007534

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510478598.XA CN106404675A (zh) 2015-08-03 2015-08-03 一种测量粗糙表面面外偏振双向反射函数的方法

Country Status (1)

Country Link
CN (1) CN106404675A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6618145B1 (en) * 2000-01-19 2003-09-09 The United States Of America As Represented By The Secretary Of The Air Force Spectroplarimetric reflectometer
CN102175650A (zh) * 2011-01-30 2011-09-07 哈尔滨工业大学 连续光谱双向反射分布函数测量装置
CN102590150A (zh) * 2012-03-01 2012-07-18 浙江大学 室内高光谱brdf测定系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6618145B1 (en) * 2000-01-19 2003-09-09 The United States Of America As Represented By The Secretary Of The Air Force Spectroplarimetric reflectometer
CN102175650A (zh) * 2011-01-30 2011-09-07 哈尔滨工业大学 连续光谱双向反射分布函数测量装置
CN102590150A (zh) * 2012-03-01 2012-07-18 浙江大学 室内高光谱brdf测定系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JIAHUI PAN ET AL.: "Results of a new polarimetric BRDF simulation of metallic surfaces", 《INFRARED PHYSIC & TECHNOLOGY》 *
M.W.HYDE IV ET AL.: "A geometrical optics polarimetric bidirectional reflectance distribution function for dielectric and metallic surfaces", 《OPTICS EXPRESS》 *
刘卿 等: "粗糙表面偏振二向反射特性建模方法研究", 《飞行器测控学报》 *
王霞 等: "粗糙表面反射辐射偏振特性研究", 《北京理工大学学报》 *

Similar Documents

Publication Publication Date Title
Garcia-Caurel et al. Application of spectroscopic ellipsometry and Mueller ellipsometry to optical characterization
Hipps et al. Applications of the photoelastic modulator to polarization spectroscopy
Naftaly Terahertz metrology
US6384916B1 (en) Parallel detecting, spectroscopic ellipsometers/polarimeters
Hug et al. Raman circular intensity differential spectroscopy. Spectra of (-)-. alpha.-pinene and (+)-. alpha.-phenylethylamine
Mehta et al. Polarized light imaging of birefringence and diattenuation at high resolution and high sensitivity
Chiu et al. Complex refractive-index measurement based on Fresnel’s equations and the uses of heterodyne interferometry
US6256097B1 (en) Ellipsometer and ellipsometry method
CN105829866B (zh) 用于使用太赫兹辐射确定片状电介质样本的至少一种性质的系统
Patat et al. Error Analysis for Dual‐Beam Optical Linear Polarimetry
Reid et al. Terahertz birefringence and attenuation properties of wood and paper
CN103335819B (zh) 一种用于高精度角锥棱镜光学检测的装置与方法
CN101498654B (zh) 使用太赫兹波的检查装置和检查方法
TWI486550B (zh) 厚度線上即時檢測之光學干涉裝置及其方法
CN102175636B (zh) 一种中草药检测鉴别的系统和方法
Wang et al. Terahertz polarization real-time imaging based on balanced electro-optic detection
US20010035954A1 (en) Method and apparatus for measuring particle size distributions using light scattering
EP3093636B1 (en) Optical polarisation modulation and detection apparatus and detection method
Snik et al. Spectral modulation for full linear polarimetry
CN204535899U (zh) 玻璃表面应力检测装置
TWI331672B (en) Focused-beam ellipsometer
US3653767A (en) Particle size distribution measurement using polarized light of a plurality of wavelengths
Sourisseau Polarization measurements in macro-and micro-Raman spectroscopies: molecular orientations in thin films and azo-dye containing polymer systems
CN106415249B (zh) 在晶片检验期间确定定位于收集孔隙中的光学元件的配置
CN100380193C (zh) 液晶面板的评价方法和评价装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170215