基于Kubelka-Munk理论的涂层表面多参量偏振双向反射分布函数模型 下载: 1377次
1 引言
随着现代科学技术的发展,涂层目标越来越受到人们的关注。由于涂层目标具有不同的折射率和表面粗糙度等特征,因此其表面散射特性很复杂,采用传统的光强探测方法会严重影响涂层目标的识别效果。对于人造目标的识别,偏振探测能够提供更多区分目标的特征信息,相比于传统的光强探测,偏振探测具有明显的优势,可以提高目标的探测精度和识别准确度。在目标偏振探测过程中,通常采用偏振度(DoP)来描述目标偏振态的变化情况,DoP被广泛用于目标探测、大气遥感以及生物医学[1]等领域。
目前,偏振度的获取方法有实验测量和理论建模2种。实验测量是指通过实验手段得到偏振度与观测几何、材料因素之间的变化关系,最常用的方法为斯托克斯(Stokes)矢量法[2];而理论建模是指通过基本理论推导得到偏振度与各变量之间的关系,可以计算任意条件下的偏振度。由于偏振度受观测几何的影响很大,需要对目标进行每个角度及方位的测量,这显然是不切实际的,因此利用有限的实验测量值建立相关的模型是必要的[3]。国内外学者对偏振度建模开展了相关研究。Thilak等[4]基于微面元理论建立了偏振双向反射分布函数(BRDF)模型,利用该模型推导出了入射光为自然光条件下散射光的偏振度表达式。微面元偏振BRDF模型忽略了体散射对偏振散射光整体能量的影响,理论计算的偏振度值与实测结果存在较大误差,且该模型局限于非镜面反射材料。Hyde等[5]在微面元偏振BRDF模型基础上引入方向半球反射率来表征体散射效应,该模型仅适用于理想的电导体表面,并且方向半球反射率是一个关于观测几何、表面粗糙度以及折射率的整个半球空间内的二重积分函数,形式较为复杂,计算量较大。巩蕾等[6]提出了一种粗糙基底上不同涂层的偏振BRDF模型,通过对不同涂层的光学特性进行研究来反演基底的光学信息。冯巍巍等[7]提出一种涂层表面的偏振BRDF模型,该模型着重考虑了涂层表面的体散射,并采用一个未知参数表示,但它不能描述体散射的产生以及入射光在涂层内部的吸收、散射过程。
针对以上问题,本文在微面元偏振BRDF模型的基础上引入体散射分量和镜向系数,并根据Kubelka-Munk(KM)理论[8-9]对体散射分量进行建模,建立了一种多参量偏振BRDF模型,推导出了散射偏振度关于表面粗糙度、复折射率实部和虚部、镜向系数、相对漫反射率系数这5个关键参数的表达式。针对典型的涂层样品,采用遗传算法从一部分训练集数据中反演出多参量偏振BRDF模型的关键参数,用另一部分测试集数据验证模型的准确性和可靠性,详细分析了镜向系数对散射偏振度的影响。
2 理论基础
2.1 BRDF
BRDF能够表征材料的散射特性和反射特性,它描述了某入射方向的光波经目标表面反射后,反射能量在上半球空间的分布情况[10],即
式中
2.2 基于KM理论的多参量偏振BRDF模型
偏振BRDF是标量BRDF更为一般的形式,可以用穆勒矩阵
式中
一般情况下,偏振BRDF包含偏振的表面散射分量
根据微面元理论,可以得到表面散射分量
式中
KM理论被广泛应用于模拟材料表面的体散射情况,le Hors等[9]在KM理论的基础上,引入一个反射率参数
式中
式中
结合(4)式和(7)式,多参量偏振BRDF模型可写为
根据(2)式和(8)式可以得到入射、散射Stokes矢量与多参量偏振BRDF之间的关系为
式中
假设入射到涂层表面的辐射是非偏的自然光,其Stokes矢量为
式中
由(10)式可知,涂层表面的散射偏振度与
3 实验
3.1 偏振探测系统
天顶角和方位角的精度均为0.2°。偏振光谱仪由便携式SVC-HR1024型地物光谱仪改造而成。由于在户外进行实验,因此为了避免自然环境变化对偏振测量结果的影响,通过设计实现自动测量,自动探测装置通过驱动步进电机带动偏振片旋转,实现目标样品在0°、60°、120°这3个偏振方向的起偏。主控系统作为偏振探测系统中的控制中心,主要实现半球扫描机构、偏振光谱仪之间的协调工作,测量数据的实时传输和存储,以及目标样品散射偏振度的全自动测量,最大限度地节省测量时间,减小太阳运动和天气变化带来的影响。
3.2 测量系统偏振定标
为了保证测量数据的精度和可靠性,需要对偏振光谱仪进行偏振定标。实验中采用高精度大动态范围可调偏振度参考光源(VPOLS-II)[16],其在0.46~2.00 μm波段内的线性偏振度调节范围可达0~0.72。偏振光谱仪定标流程图如
3.3 实验测量
由于户外实验条件受自然因素的影响很大,因此要求测量的时间越短越好,每种涂层样品的实际测量时间为17 min,测得绿漆和黑漆涂层目标的太阳入射天顶角分别为48.514°和41.476°,测量示意图如
4 模型参数反演与分析
4.1 模型参数反演结果
为了研究所提多参量偏振BRDF模型的准确性,分别选取绿漆和黑漆涂层表面的散射偏振度进行分析。在该模型中共有5个未知参数,最佳模型参数选择的标准是模型仿真值与实验测量值的标准均方差最小,因此,构建优化目标函数为
式中
表 1. 黑漆和绿漆涂层表面的参数反演结果
Table 1. Parameters inversion for black and green painted surfaces
|
利用
图 6. 不同相对方位角时DoLP仿真结果与实测数据的比较。(a)绿漆涂层, Δ?=120°;(b)黑漆涂层, Δ?=120°;(c)绿漆涂层, Δ?=180°;(d)黑漆涂层, Δ?=180°
Fig. 6. Comparison between simulated and measured degree of linear polarization at different relative azimuth angles. (a) Green painted surface, Δ?=120°; (b) black painted surface, Δ?=120°; (c) green painted surface, Δ?=180°; (d) black painted surface, Δ?=180°
4.2 镜向系数对仿真结果的影响
通过以上分析可知,不同涂层表面的镜向系数存在很大差异,且在基于KM理论的多参量偏振BRDF模型中,镜向系数对模型的仿真结果影响很大。
图 7. 含有镜向系数ks与不含镜向系数ks时DoLP仿真结果对比。(a)绿漆涂层, Δ?=120°;(b)黑漆涂层, Δ?=120°;(c)绿漆涂层, Δ?=180°;(d)黑漆涂层, Δ?=180°
Fig. 7. Comparison on simulated degree of linear polarization with and without mirror parameter ks. (a) Green painted surface, Δ?=120°; (b) black painted surface, Δ?=120°; (c) green painted surface, Δ?=180°; (d) black painted surface, Δ?=180°
为了更加直观地反映
表 2. 不同涂层表面的2种DoLP均方根误差
Table 2. Root-mean-square error of two degree of linear polarization for different painted surfaces
|
5 结论
提出了一种基于KM理论的多参量偏振BRDF模型,用于表征涂层表面的散射偏振特性。该模型利用KM理论来描述体散射分量,并引入镜向系数来反映表面散射分量的贡献程度,从而改进了传统的偏振BRDF模型,使其更加准确、合理。通过对比模型仿真和实测数据后发现,引入镜向系数使得模型的仿真值与实测值吻合得更好,提高了仿真精度,在相对方位角为120°的平面内,效果更明显。可见,基于KM理论的多参量偏振BRDF模型可以准确地描述涂层表面的散射偏振特征。需要指出的是,基于半经验的偏振BRDF模型可用于描述目标表面的散射偏振特性,若要获取更精确的偏振BRDF模型,需要考虑不同材料之间的模型参量存在的差异(可利用实际测量数据拟合得到)。另外,对于自然地物的偏振BRDF模型还需要考虑后向散射效应的影响,这是本课题组后续研究的方向。
[3] 陈超, 赵永强, 罗丽, 等. 基于遗传LM算法的涂层目标光谱偏振BRDF建模分析[J]. 光谱学与光谱分析, 2010, 30(3): 729-734.
陈超, 赵永强, 罗丽, 等. 基于遗传LM算法的涂层目标光谱偏振BRDF建模分析[J]. 光谱学与光谱分析, 2010, 30(3): 729-734.
陈超, 赵永强, 罗丽, 等. 基于遗传LM算法的涂层目标光谱偏振BRDF建模分析[J]. 光谱学与光谱分析, 2010, 30(3): 729-734.
[6] 巩蕾, 吴振森. 不同涂层对微粗糙基底的极化光散射影响[J]. 光学学报, 2011, 31(10): 1029001.
巩蕾, 吴振森. 不同涂层对微粗糙基底的极化光散射影响[J]. 光学学报, 2011, 31(10): 1029001.
巩蕾, 吴振森. 不同涂层对微粗糙基底的极化光散射影响[J]. 光学学报, 2011, 31(10): 1029001.
[7] 冯巍巍, 魏庆农, 汪世美, 等. 涂层表面偏振双向反射分布函数模的模型研究[J]. 光学学报, 2008, 28(2): 290-294.
冯巍巍, 魏庆农, 汪世美, 等. 涂层表面偏振双向反射分布函数模的模型研究[J]. 光学学报, 2008, 28(2): 290-294.
冯巍巍, 魏庆农, 汪世美, 等. 涂层表面偏振双向反射分布函数模的模型研究[J]. 光学学报, 2008, 28(2): 290-294.
[9] le HorsL, HartemannP, DolfiD, et al. Phenomenological model of paints for multispectral polarimetric imaging[C]. SPIE, 2001, 4370: 94- 105.
le HorsL, HartemannP, DolfiD, et al. Phenomenological model of paints for multispectral polarimetric imaging[C]. SPIE, 2001, 4370: 94- 105.
le HorsL, HartemannP, DolfiD, et al. Phenomenological model of paints for multispectral polarimetric imaging[C]. SPIE, 2001, 4370: 94- 105.
[10] Nicodemus FE, Richmond JC, Hsia JJ, et al. Geometrical considerations and nomenclature for reflectance[M]. Washington, D.C.: National Bureau of Standards, 1977: 1- 7.
Nicodemus FE, Richmond JC, Hsia JJ, et al. Geometrical considerations and nomenclature for reflectance[M]. Washington, D.C.: National Bureau of Standards, 1977: 1- 7.
Nicodemus FE, Richmond JC, Hsia JJ, et al. Geometrical considerations and nomenclature for reflectance[M]. Washington, D.C.: National Bureau of Standards, 1977: 1- 7.
[12] 张朝阳, 程海峰, 陈朝辉, 等. 伪装涂层材料的二向偏振散射研究[J]. 光学技术, 2009, 35(1): 117-119.
张朝阳, 程海峰, 陈朝辉, 等. 伪装涂层材料的二向偏振散射研究[J]. 光学技术, 2009, 35(1): 117-119.
张朝阳, 程海峰, 陈朝辉, 等. 伪装涂层材料的二向偏振散射研究[J]. 光学技术, 2009, 35(1): 117-119.
Zhang Z Y, Cheng H F, Chen Z H, et al. Study of bi-directional polarimetric scattering of camouflage coatings[J]. Optical Technique, 2009, 35(1): 117-119.
[13] Priest RG, Germer TA. Polarimetric BRDF in the microfacet model: theory and measurement[C]. Proceedings of the Meeting of the Military Sensing Symposia Specialty Group on Passive Sensors, 2000, 1: 169- 181.
Priest RG, Germer TA. Polarimetric BRDF in the microfacet model: theory and measurement[C]. Proceedings of the Meeting of the Military Sensing Symposia Specialty Group on Passive Sensors, 2000, 1: 169- 181.
Priest RG, Germer TA. Polarimetric BRDF in the microfacet model: theory and measurement[C]. Proceedings of the Meeting of the Military Sensing Symposia Specialty Group on Passive Sensors, 2000, 1: 169- 181.
[14] 王启超, 汪家春, 赵大鹏, 等. 光谱偏振探测对伪装网的识别研究[J]. 激光与红外, 2013, 43(11): 1260-1264.
王启超, 汪家春, 赵大鹏, 等. 光谱偏振探测对伪装网的识别研究[J]. 激光与红外, 2013, 43(11): 1260-1264.
王启超, 汪家春, 赵大鹏, 等. 光谱偏振探测对伪装网的识别研究[J]. 激光与红外, 2013, 43(11): 1260-1264.
Wang Q C, Wang J C, Zhao D P, et al. Research on the detection of camouflage net with spectral polarization detection system[J]. Laser & Infrared, 2013, 43(11): 1260-1264.
[15] 李新. 可见-短波红外双向反射分布函数测量系统的研制与应用[D]. 合肥: 中国科学院合肥物质科学研究院, 2008.
李新. 可见-短波红外双向反射分布函数测量系统的研制与应用[D]. 合肥: 中国科学院合肥物质科学研究院, 2008.
李新. 可见-短波红外双向反射分布函数测量系统的研制与应用[D]. 合肥: 中国科学院合肥物质科学研究院, 2008.
LiX. Development and application of VNIR-SWIR bidirectional reflectance distribution function (BRDF) measurement system[D]. Hefei: Hefei Institutes of Physical Science,Chinese Academy of Sciences, 2008.
LiX. Development and application of VNIR-SWIR bidirectional reflectance distribution function (BRDF) measurement system[D]. Hefei: Hefei Institutes of Physical Science,Chinese Academy of Sciences, 2008.
LiX. Development and application of VNIR-SWIR bidirectional reflectance distribution function (BRDF) measurement system[D]. Hefei: Hefei Institutes of Physical Science,Chinese Academy of Sciences, 2008.
[16] 康晴, 李健军, 陈立刚, 等. 大动态范围可调线性偏振度参考光源检测与不确定度分析[J]. 光学学报, 2015, 35(4): 0412003.
康晴, 李健军, 陈立刚, 等. 大动态范围可调线性偏振度参考光源检测与不确定度分析[J]. 光学学报, 2015, 35(4): 0412003.
康晴, 李健军, 陈立刚, 等. 大动态范围可调线性偏振度参考光源检测与不确定度分析[J]. 光学学报, 2015, 35(4): 0412003.
[17] 孙学金, 李浩, 唐丽萍. 可见光和红外波段大气体散射偏振度特性研究[J]. 光学学报, 2011, 31(5): 0501001.
孙学金, 李浩, 唐丽萍. 可见光和红外波段大气体散射偏振度特性研究[J]. 光学学报, 2011, 31(5): 0501001.
孙学金, 李浩, 唐丽萍. 可见光和红外波段大气体散射偏振度特性研究[J]. 光学学报, 2011, 31(5): 0501001.
Article Outline
杨敏, 方勇华, 吴军, 崔方晓, 李大成, 吴越, 王涛. 基于Kubelka-Munk理论的涂层表面多参量偏振双向反射分布函数模型[J]. 光学学报, 2018, 38(1): 0126002. Min Yang, Yonghua Fang, Jun Wu, Fangxiao Cui, Dacheng Li, Yue Wu, Tao Wang. Multiple-Component Polarized Bidirectional Reflectance Distribution Function Model for Painted Surfaces Based on Kubelka-Munk Theory[J]. Acta Optica Sinica, 2018, 38(1): 0126002.