用于快速BRDF测量的子孔径扫描傅里叶变换系统 下载: 856次
1 引言
近年来,材料表面散射特性的模型与测量在许多研究领域都得到了积极的应用,如材料物理学[1-2]、工业工程[3-4],以及图像渲染技术[5-6]。科技的发展对高精度、快速的散射测量技术提出了更高的要求。
以球面度的倒数为单位的双向反射分布函数(BRDF)可以描述一般的光学性质,这一术语是由美国国家标准局(NIST)的研究人员在1977年定义的[7],它取决于入射光的性质(波长、偏振态)以及入射光和散射光的方向。对BRDF的研究主要集中在散射通量随散射角度的变化上。BRDF的测量需要精密的仪器,这些仪器主要是由计量实验室或仪器实验室开发的,它们被称为角度细分光谱辐射计。这些装置面临的困难是以一定方向照射样品时如何测量散射在各个方向的光。光谱辐射计需要测量无数次才能够保证测量模型更加精确,且由于所测量的信号总是很小,因此这些装置面临的更大挑战在于如何使探测器足够灵敏,以达到足够大的信噪比。自2000年以来,这些方面都有了改进,随着精密机械结构的发展,测量所需时间也逐渐减少。欧洲计量研究项目已经运行,以改善参考和测量技术[8]。为了使设备更紧凑,French Pierre、French Pierre Boher等开发了使用傅里叶变换光学方法的BRDF测量设备,实现了测角角度细分光谱辐射计测量难以实现的测角分辨率[9](视场角高达±88°,是一个超广角傅里叶光学系统),将BRDF测量技术又提升了一个量级。Kim首次提出采用子孔径拼接方法测量大口径平面镜,并将被测面划分为多个小孔径进行测量,获得的子图像在扫描探测器的整个范围内创建了一个合成孔径,从而得到高分辨率图像[10-11]。
本文提出一种基于子孔径扫描傅里叶变换系统的 BRDF 快速测量方法。针对超广角系统中较大像差难以校正的情况,提出使用一个普通单透镜作为傅里叶光学元件的新思路。由于单透镜在入射角40°以上时工作效率较低,故采用子孔径扫描方法来克服这一角度范围的限制。本文叙述了傅里叶变换方法测量角散射光谱的方法与子孔径扫描扩展探测角度的方法,并通过实验对BRDF数据采集过程进行分析。为了验证该方法的有效性,本文测量了三个不同的样品,并将其数据与角度细分光谱辐射计BRDF测量得到的数据进行了比较。结果表明,该技术可以为今后的快速、简便的BRDF测量提供更好的思路。
2 用于BRDF测量的傅里叶变换法
2.1 BRDF定义
双向反射分布函数(BRDF)是表征材料表面的散射特性的一种最常用的方法。BRDF用特定波长下的反射辐亮度Lr(θi,φi,θr,φr,λ)与入射辐照度Ei(θi,φi,λ)的比值表示,定义为
该式描述了在入射方向(θi,φi)和给定波长λ条件下由表面材料反射,在反射方向(θr,φr)的辐射强度的变化。
传统的BRDF测量是在特定的照明条件下,测量不同观测角度下的辐射量,但是使用这样的设备通过密集测量来获取更精确的BRDF是非常耗时的。如
图 1. BRDF测量原理及BiRT-2.0 BRDF测量系统
Fig. 1. Principle of BRDF measurement and the BiRT-2.0 BRDF measurement system
2.2 傅里叶变换原理
与直接测量技术相比,只需简单地应用二维傅里叶变换,就可以数值计算出样品的远场光散射模式[12]。傅里叶变换具有独特的优点:1)可以在一次测量中获得大角度范围的散射图案;2)由于充分利用了阵列探测器,信噪比极高;3)可以获得微米量级纹理的光散射图案。
本文的目的是使用单个傅里叶变换透镜,在合理的测量时间测量光谱BRDF。其原理基于光学傅里叶变换,使用普通的单透镜,将角度场转换为傅里叶平面图[12]。傅里叶平面图被一个阵列探测器接收,能够实现视场内的所有散射光的同时测量。
图 2. 空间和傅里叶域中的光场分布坐标。(a)空间坐标;(b)傅里叶域坐标
Fig. 2. Coordinates of field distributions in the spatial and Fourier domains. (a) Space coordinates; (b) Fourier domain coordinates
傅里叶域中每个单元的位置由平行投影的对应方向余弦决定。对应投影方向的单元中心空间频率为
式中:Ez(x,y)为初始场E(x0,y0)的远场分布;λ为波长。电场分布Ez(x,y)与初始电场的二维傅里叶变换(FT)成正比,
式中:u和v分别表示x和y方向上的空间频率;θ和φ分别表示散射矢量的天顶角和方位角。严格地说,(3) 式仅在近轴近似中有效,采用广角傅里叶变换方法时会产生大的测量误差。更一般地说,空间频率和散射角之间是一一对应的。因此,在傅里叶平面上测量得到的简单傅里叶变换电场E(x,y)给出了样品的角散射光谱。
2.3 子孔径扫描傅里叶变换系统
傅里叶光学器件是一个透镜(或一组透镜),它可以收集待测微面元散射的光,然后将光线重新聚焦在一个平面(称为傅里叶平面)上。
傅里叶平面中的单元取决于散射光相对该平面的入射角和方位角。由半径r和方位角φ定义的傅里叶平面上的基本表面照明度,即
式中:L(r,φ)是待测微面元的亮度;S是待测微面元的表面积。
如
其中S0为光源垂直入射时的光斑面积,则辐射亮度可以表示为
cos4(θ)随入射角的变化如
因此提出了子孔径扫描傅里叶变换系统的方法,该方法不用考虑超广角光学系统中的大像差问题,仅使用单个普通透镜,通过球面子孔径扫描来扩展探测角度范围。散射矢量的最大可测角度θmax由透镜的数值孔径(NA)决定,数值孔径的表达式为
式中:D为透镜口径。
为了克服这一角度限制,本文采用子孔径扫描方法。在本应用中,傅里叶变换透镜和探测器采取同轴安装,如
线阵探测器已经足以满足本方案的验证需求,但是本文选择面阵探测器来清晰显示散射光谱。这是便于在以后的研究中,将这项工作与球面子孔径图像拼接相结合,以进一步简化数据采集的扫描过程。
如
3 实验验证
3.1 实验原理
本文选择3种不同样品,使用角度细分光谱辐射计和子孔径扫描傅里叶变换系统两种方法进行BRDF测量,样品分别为实验室标准漫反射白板、镜面和划痕镜面样品,如
实验验证过程如
图 6. 待测样品以及显微镜测量的划痕尺寸
Fig. 6. Samples to be tested and the scratch measured under microscope
测量机械系统主要由1台ZX110-300、1台ZX110-100和1台RS01-DZ的精密电控转台、探测器与光源旋臂、样品托等组成,可以测试材料不同光照与探测视点的变化。照明旋臂与探测器悬臂的角度分辨率为36″,重复定位精度为18″,精密电控转台的角度分辨率为6',重复定位精度为36″。为了实现自动测量,本课题借助北京联英精机科技公司编制的运动卡应用软件将光谱仪测量软件、探测器图像采集软件与运动控制软件整合在一起,实现了运动路线的主动控制与测量数据的自动保存。
3.2 BRDF数据获取
根据亮度守恒定律,可以得到反射辐射亮度与探测器入瞳面处光谱辐射照度之间的关系为
式中:Er为探测器入瞳处的光谱辐照度;ΔApupol为探测器入瞳的面积;ΔAs为样品表面照明面源面积;ΔΩ是探测立体角。
根据定义,可得
已知实际测量过程中样品的探测距离和照明面积,可得BRDF测量方程为
式中:Ei是材料表面在光照下的辐照度,是通过测量垂直入射辐照度并乘以入射角的余弦值获得的;R为探测距离。
3.3 数据结果分析
在角度细分光谱辐射计对比实验中,测量整个半球空间的BRDF数据量很大。为了有效对比方法的有效性,本文只选取一个入射角度、一个反射方位角,来测量反射方位角改变的散射情况。测试的角度选取如下:光源天顶角为45°,探测器方位角为180°,探测器天顶角为0°~90°,角度细分光谱辐射计的探测步长为1°(每1°采集一次辐射数据),子孔径扫描傅里叶变换系统步长为10°。测量波长为654.988 nm。通过调整样品和检测装置的相对位置,基于3.2节介绍的数据处理方法,可以得到标准漫反射白板、镜面和划痕镜面样品的角散射光谱。结果如
为验证本文方法的有效性,定义了两个比较数据结果相对误差的评价函数,如
式中:n为实验数据组数;fGS(θi,θr,φr)和fFT(θi,θr,φr)表示Gonio角度细分光谱辐射计和子孔径扫描傅里叶变换方法测量的BRDF值;ERR表示两种方法获得实验结果的相对误差,ERR越小,两种方法的结果越接近。
据,此角度分辨率可以达到0.005°。由于子孔径扫描傅里叶变换系统的步长为角度细分光谱辐射计系统步长的10倍,测量时间是角度细分光谱辐射计的1/10,故所需的扫描时间更短,效率更高。结果表明,本文建立的子孔径扫描傅里叶变换系统对测量BRDF具有良好的效果。
4 结论
基于傅里叶变换原理,通过单次测量即可获得多角度的角散射信息。使用单个普通透镜进行测量,不仅扩展了测量视场,还具有工作距离远、价格低廉、光学系统简单等优点。使用子孔径扫描和傅里叶空间数值拼接克服了单透镜傅里叶变换的角度限制。在分析了实验方案和BRDF数据采集过程后,将角度细分光谱辐射计中的光纤角度细分光谱替换为薄透镜与阵列探测器组合。通过对两种系统得到的标准漫反射白板、镜面和划痕镜面样品的BRDF数据进行对比,验证了该方法的有效性。实验结果表明,傅里叶变换法与原模拟分析和角度细分光谱辐射计结果吻合良好,两种系统的相对误差均在1%以下。此外,该方法可以获得更小的角度分辨率,检测用时更短,效率更高,由此验证了子孔径扫描傅里叶变换方法在BRDF测量中以及在角散射测量中的可行性。
[1] 黄聪, 游兴海, 张彬. 光学元件表面洁净度对其表面散射特性的影响[J]. 红外与激光工程, 2019, 48(1): 0120002.
[2] 史卫朝, 郑建明, 李言, 等. 加工表面双向反射分布函数的测量与建模[J]. 光学学报, 2018, 38(10): 1029001.
[3] 陈川, 易维宁, 崔文煜, 等. 大气散射对地表双向反射分布函数反演的影响与修正[J]. 光学学报, 2018, 38(7): 0729001.
[4] 汪夏, 张雅声, 徐灿, 等. 空间目标褶皱材质BRDF建模方法[J]. 光学学报, 2019, 39(10): 1024001.
[5] Chermain X, Claux F, Mérillou S. Glint rendering based on a multiple-scattering patch BRDF[J]. Computer Graphics Forum, 2019, 38(4): 27-37.
[6] Raymond B, Guennebaud G, Barla P. Multi-scale rendering of scratched materials using a structured SV-BRDF model[J]. ACM Transactions on Graphics, 2016, 35(4): 1-11.
[7] Nicodemus FE, Richmond JC, Hsia JJ, et al. Geometrical considerations and nomenclature for reflectance[M] ∥Radiometry. Sudbury, MA: Jones and Bartlett Publishers, Inc., 1992: 94- 145.
[9] Simonot L, Obein G, Bringier B, et al. Modeling, measuring, and using BRDFs: significant French contributions[J]. Journal of the Optical Society of America A, 2019, 36(11): C40.
[10] 于瀛洁, 齐特, 武欣. 大尺寸光学元件在位动态干涉拼接测量系统[J]. 光学精密工程, 2017, 25(7): 1764-1770.
[11] 刘昂, 何宇航, 李强, 等. 大口径光学元件功率谱密度的拼接干涉检测[J]. 中国激光, 2019, 46(2): 0204008.
[12] Ding H F, Wang Z, Nguyen F, et al. Fourier transform light scattering of inhomogeneous and dynamic structures[J]. Physical Review Letters, 2008, 101(23): 238102.
[13] 冯绚, 韩昌佩, 邹曜璞, 等. 红外傅里叶光谱仪在轨光谱定标算法研究[J]. 光学学报, 2019, 39(6): 0630002.
Article Outline
陆敏, 王治乐, 高萍萍, 张树青, 郭继锴. 用于快速BRDF测量的子孔径扫描傅里叶变换系统[J]. 光学学报, 2020, 40(13): 1329001. Min Lu, Zhile Wang, Pingping Gao, Shuqing Zhang, Jikai Guo. A Sub-Aperture Scanning Fourier Transform System for Fast BRDF Measurements[J]. Acta Optica Sinica, 2020, 40(13): 1329001.