Study on automatic measurement method of elastic modulus for optical glass
-
摘要: 针对现有光学玻璃弹性模量测量周期长、测量过程复杂等问题，提出了一种利用计算机图像处理技术快速测量小样品光学平板玻璃弹性模量的方法。基于接触力学理论和牛顿环干涉原理推导出测量光学玻璃弹性模量的解析式，利用图像处理技术的优势，通过摄像头拍摄牛顿环干涉图像，设计MFC监控与单张图片拍摄界面程序，自动测量牛顿环干涉图像中心黑斑半径，最后根据光学玻璃弹性模量与牛顿环干涉图像中心黑斑半径及牛顿环中心应力之间的关联关系，实现小样品光学玻璃弹性模量的快速测量。实验结果表明：在应力31.17 N～55.11 N范围内，光学玻璃弹性模量的测量相对误差不超过±8.8%；在应力55.11 N～71.07 N范围内，测量相对误差不超过±16%。Abstract: In order to solve the problems of long measuring period and complex measurement process of elastic modulus for optical glass, a method of quickly measuring the elastic modulus of small sample optical flat glass by computer image processing technology was proposed. Firstly, based on contact mechanics theory and Newton ring interference principle, the analytical formula for measuring the elastic modulus of optical glass was deduced. Secondly, taking advantage of the image processing technology, the Newton ring interference image was taken by camera, and the interface program of MFC monitoring and single image shooting was designed to automatically measure the center black spot radius of Newton ring interference image. Finally, according to the relationship between the elastic modulus of optical glass, the center black spot radius of Newton ring interference image, and the center stress of Newton ring, the fast measurement of the elastic modulus of small sample optical glass was realized. The experimental results show that the relative error of the elastic modulus for optical glass is no more than ± 8.8% in the range of stress 31.17 N~55.11 N, and ± 16% in the range of stress 55.11 N~71.07 N.
-
Keywords:
- applied optics /
- elastic modulus /
- image processing /
- contact mechanics /
- Newton ring
-
引言
高光谱遥感器是指利用分光技术获取连续窄带光谱的遥感器，其通道光谱带宽Δλ与中心波长λ的关系为Δλ/λ≈0.01[1]。在高光谱遥感器使用前，需要进行光谱定标和辐射定标。光谱定标的目的是建立高光谱遥感器光谱通道与不同入瞳波长的对应关系，辐射定标时建立高光谱遥感器的输出值与相应已知的、用国际单位(SI)表示的辐射量基准之间定量关系[2-3]。光谱定标和辐射定标是保证高光谱遥感器获取数据的科学可信度和最大限度的定量化应用基础。
高光谱遥感器的光谱定标和辐射定标通常在实验室内进行，其中光谱定标通常利用基于元素特征谱的光谱定标参考光源[4](低压汞灯和高压钠灯等特征谱线灯安装在积分球里面)或经过波长校准的单色仪或可调谐激光器[5]等精细扫描的方法实现，辐射定标通常采用光谱平坦性较好的标准灯板系统或传统卤素灯积分球光源完成。由于高光谱遥感器工作环境复杂，通常需要在野外进行观测，尤其是长期观测，受本身老化和外界干扰等因素的影响，其光谱中心波长漂移会发生漂移，其辐射响应特性会不断衰减[6]。前期高光谱遥感器应用研究表明，中心波长漂移与通道光谱带宽比值在超出1%情况下，在光谱形状变化剧烈的谱段，辐射测量的偏差有可能超过5%，甚至达到20%[7]。因此有必要对高光谱遥感器光谱辐射定标精度在现场环境下进行精度验证，实现光谱定标精度优于±0.5 nm，确保辐射定标不确定小于5%。
目前在现场环境下，可利用大气廓线中的太阳夫琅禾费线以及氧气、二氧化碳和水汽等气体吸收光谱特征对高光谱遥感器进行光谱定标精度验证，例如Gao等人采用光谱匹配技术实现了机载高光谱成像光谱仪PHILLS的光谱定标[8-9]。通过选择大气干洁、地表开阔的场地，利用太阳光经过不同反射率漫反射参考板或靶标的反射信号对高光谱遥感器进行辐射定标精度验证[10-11]。这种定标方法较精确，但是在实际操作中，因为很难能找到稳定、晴好、无云的天气条件，进而导致定标数据质量下降。
现阶段尚未有高精度的光谱定标和辐射定标实现现场环境下光谱辐射定标验证方法。本文设计了一种高光谱遥感器现场光谱辐射定标精度的验证方法，以SVC光谱辐亮度计为例，利用新型具有高稳定性的光谱可调积分球参考光源STIS(a spectrum-tunable reference light source) [12]，在待验证波段内输出光谱辐亮度形状单调上升或单调下降状态，通过比对多波段辐亮度标准传递探测器MRSTD(a muti-band radiance standard transfer detector)[11]和SVC光谱辐亮度计辐亮度测量值的相对偏差，实现SVC光谱辐亮度计光谱辐射定标的精度验证和校正，其中光谱定标验证精度小于±0.2 nm, 辐射定标验证精度优于5%。
1 精度验证的方法
本文使用STIS和MRSTD验证高光谱光学遥感器实验室光谱定标精度。STIS采用多种波长的LED密集阵列发光单元和卤钨灯发光单元，在400 nm~800 nm具有分立波段输出功能，可输出半高带宽约20 nm~40 nm近似高斯形状的光谱辐亮度状态。STIS采用循环水冷温控和高精度可编程直流电源驱动技术，保障了光源的稳定性。STIS角度均匀性在±60°以内优于99.4%，面均匀性优于99.5%，2 h内稳定性优于99.8%。STIS积分球内径250 mm，出光口直径为60 mm，如图 1所示。
MRSTD是基于低温辐射计传递的多波段辐亮度标准传递探测器[13]，用于绝对辐射定标以及监测参考目标的辐射特性变化。MRSTD的工作波段有412 nm、443 nm、490 nm、520 nm、565 nm、670 nm、750 nm和865 nm，其通道带宽约为(10±2)nm。MRSTD采用激光导入积分球的单色面光源进行光谱辐射定标，其中心波长定标准确度优于±0.1 nm，其绝对不确定度 < 1.26%。
设计的高光谱遥感器光谱定标的精度验证方法如下：在MRSTD各通道的带内区域，利用STIS输出光谱形状单调上升和单调下降的光谱辐亮度状态。使用MRSTD和SVC光谱辐亮度计交替测量STIS输出的光谱辐亮度，如图 2所示。通过比较MRSTD和SVC光谱辐亮度计辐亮度测量结果的相对偏差，分析SVC光谱辐亮度计的光谱辐射定标精度。
利用标准传递探测器第i通道响应DNi和其带内响应度为Ri，测量光谱可调积分球参考光源输出辐亮度Lc(λi)作为参考值，如(1)式和(2)式所示：
$$ {{L}_{c}}\left( {{\lambda }_{i}} \right)={D{{N}_{i}}}/{{{R}_{i}}}\; $$ (1) $$ {{R}_{i}}=\int_{{{\lambda }_{l}}}^{{{\lambda }_{u}}}{R\left( {{\lambda }_{i}} \right)}\text{d}\lambda $$ (2) SVC光谱辐亮度计在350 nm~1 000 nm范围内的光谱分辨率，即半高带宽优于3.5 nm，小于MRSTD通道带宽。根据前期研究[14]，SVC光谱辐亮度计通道光谱响应度一般接近于sinc函数，与MRSTD通道的类似高斯光谱响应度曲线形状不同，因此两种仪器测量的辐亮度无法直接进行比对。
解决方法：首先对SVC光谱辐亮度计测量的光谱辐亮度曲线进行线性插值，插值波长间隔为0.2 nm；其次卷积SVC光谱辐亮度计测量的光谱辐亮度Lsvc(λi)和MRSTD第i通道的相对光谱响应度R(λi)，得到光谱辐亮度计对应MRSTD第i通道辐亮度测量值LBSR(λi)：
$$ {{L}_{BSR}}\left( {{\lambda }_{i}} \right)=\int_{{{\lambda }_{l}}}^{{{\lambda }_{u}}}{{{L}_{SVC}}\left( {{\lambda }_{i}} \right)R\left( {{\lambda }_{i}} \right)}{\text{d}\lambda }/{{{R}_{i}}}\; $$ (3) 计算MRSTD的参考辐亮度Lc(λi)和SVC光谱辐亮度计辐亮度LBSR(λi)测量值之间的相对偏差Δi：
$$ {{\Delta }_{i}}=\frac{{{L}_{BSR}}\left( {{\lambda }_{i}} \right)-{{L}_{c}}\left( {{\lambda }_{i}} \right)}{{{L}_{c}}\left( {{\lambda }_{i}} \right)}\times 100% $$ (4) 当SVC光谱辐亮度计存在波长偏移时，参考辐亮度Lc(λi)和辐亮度LBSR (λi)测量值之间的相对偏差Δi差异明显。在MRSTD通道对应(10±2)nm光谱范围，SVC光谱辐亮度计的光谱带宽可默认为不发生变化。此时参考前期研究方法[8-9]，利用光谱匹配技术，即对SVC光谱辐亮度计测量的光谱辐亮度LSVC(λi)进行波长平移，判断Δi的变化情况。
根据测量不确定度合格评定通用比对标准[15]，当SVC光谱辐亮度计不存在波长偏移时，参考辐亮度Lc(λi)和辐亮度LBSR (λi)测量值的相对偏差Δi应满足下式：
$$ {\Delta _i} \le \sqrt {{u^2}\left( {{\rm{SVC}}} \right) + {u^2}\left( {{\rm{MRSTD}}} \right)}<5\% $$ (5) 式中：u(SVC)为SVC光谱辐亮度计的辐射定标不确定度；u(MRSTD)为MRSTD辐射定标不确定度。设计最大相对偏差Δi＜5%为判据，认为此时平移后的波长为SVC光谱辐亮度计真实的波长，进而验证SVC光谱辐亮度计的光谱辐射定标精度。
2 实验及数据处理
利用STIS在MRSTD第2通道的带内区域输出4种光谱辐亮度度状态，对应的峰值中心波长分别为474.5 nm、475.9 nm、496.8 nm和499.7 nm。MRSTD测量的辐射亮度值分别为4.409μW·cm-2·sr-1、6.0377μW·cm-2·sr-1、13.950 7 μW·cm-2·sr-1和13.967 7 μW·cm-2·sr-1。图 3(a)为SVC光谱辐亮度计测量的光谱辐亮度曲线。图 3(b)为MRSTD第2通道490 nm波段的相对光谱响应度曲线，其辐亮度测量值为参考值，如表 1的第一行所示。
![]()
图 3 SVC光谱辐亮度计测量的光谱辐亮度曲线和MRSTD 490 nm波段的相对光谱响应度Figure 3. Measured spectral radiance curve by SVC and relative spectral responsivity of MRSTD 490 nm band表 1 490 nm波段各波长相对偏移量(%) Table 1. Relative deviation of each wavelength in 490 nm band波长平移量/nm 相对偏差Δi 474.5 nm 475.9 nm 496.8 nm 499.7 nm 1.0 2.68 3.13 -4.65 -5.80 0.8 1.10 1.62 -3.66 -4.43 0.6 -0.47 0.14 -2.69 -3.08 0.4 -2.01 -1.32 -1.76 -1.74 0.2 -3.53 -2.76 -0.86 -0.43 0 -5.04 -4.18 0 0.85 -0.2 -6.53 -5.58 0.83 2.12 SVC光谱辐亮度计的波长不作平移时，计算MRSTD的参考辐亮度Lc(λi)和SVC光谱辐亮度计辐亮度LBSR (λi)测量值之间的相对偏差Δi最大为5.04%。将SVC波长向长波方向平移0.4 nm后，针对4种光谱辐亮度状态，计算的相对偏差Δi均小于2.01%。SVC光谱辐亮度计波长平移(长波方向为正)与相对偏差Δi的变化情况如图 4所示，相对偏差Δi变化的具体数值见表 1。
![]()
图 4 490nm波段MRSTD与SVC光谱辐亮度计波长平移的相对偏差Δi变化情况Figure 4. Changes of relative deviation of MRSTD and SVC wavelength translation in 490 nm band从表 1可以看出，针对STIS输出的474.5 nm和475.9 nm两种光谱辐亮度状态，SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律趋向单调上升，对于STIS输出的496.8 nm和499.7 nm2种光谱辐亮度状态，SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律呈单调下降趋势。通过两种STIS分别输出光谱形状单调上升和单调下降的两种光谱辐亮度状态，结合SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律和相对偏差Δi优于5%的判据，得到SVC光谱辐亮度计的波长偏移了0.4 nm。
SVC光谱辐亮度计760 nm波段的光谱和辐射定标精度验证情况如图 5所示，SVC的波长偏移了1.8 nm。SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律明显。760 nm波段测量与490 nm波段测量结果相比，SVC光谱辐亮度计波长平移与相对偏差Δi的变化规律呈单调上升和下降趋势较为平坦，主要是由于与峰值中心波长474.5 nm、475.9 nm、496.8 nm和499.7 nm的光谱辐亮度状态相比，STIS输出峰值中心波长748.3 nm和759.9 nm状态的光谱形状单调上升和单调下降的变化较为平缓。
![]()
图 5 760 nm波段MRSTD与SVC光谱辐亮度计波长平移的相对偏差Δi变化情况Figure 5. Changes of relative deviation of MRSTD and SVC wavelength translation in 760 nm band同样我们验证了SVC光谱辐亮度计在350 nm~800 nm光谱范围内MRSTD其余波段的光谱辐射定标精度。受制于STIS内部LED的种类，在800 nm以后无法输出半高带宽约20 nm~40 nm近似高斯形状的光谱辐亮度状态。
3 结论
本文以SVC野外光谱辐亮度计模拟高光谱遥感器，通过STIS和MRSTD设计了一种高光谱光学遥感器光谱辐射定标的精度验证方法。借助2种STIS在待验证的光谱通道内分别输出光谱形状单调上升和单调下降的两种光谱辐亮度状态，利用SVC光谱辐亮度计波长平移与MRSTD的参考辐亮度和SVC辐亮度测量值相对偏差Δi的变化规律，可有效验证SVC光谱辐亮度计光谱和辐射定标精度，其中光谱定标验证精度优于0.2 nm，辐射定标精度优于5%。这种光谱辐射定标精度验证方法可用于高光谱遥感器350 nm~800 nm光谱范围内实验室定标和现场定标，同时也为光谱连续型光学遥感器与窄带滤光片型光学遥感器之间的辐亮度测量提供一种参考方法。
-
![]()
图 2 牛顿环应力测量及图像采集装置结构示意图
Figure 2. Structure diagram of Newton ring stress measurement and image acquisition device
![]()
图 4 去除中心黑色小联通区域后的二值图像
Figure 4. Binary image after removing central black small connection area
![]()
图 7 光学平板玻璃弹性模量实验值与标准值对比
Figure 7. Comparison of experimental and standard values of elastic modulus for optical flat glass
![]()
图 8 光学平板玻璃弹性模量相对误差绝对值和应力的关系
Figure 8. Relationship between absolute value of relative error and stress of elastic modulus for optical flat glass
![]()
图 9 牛顿环中心黑斑半径和应力的关系
Figure 9. Relationship between radius of black spot and stress in center of Newton ring
表 1 光学平板玻璃弹性模量测量结果
Table 1 Measurement results of elastic modulus for optical flat glass
序号 应力/N 黑斑半径/mm 弹性模量实验值/${\rm{N}} \cdot {\rm{m} }{ {\rm{m} }^{\rm{-2} } }$ 弹性模量标准值/${\rm{N}} \cdot {\rm{m} }{ {\rm{m} }^{\rm{-2} } }$ 绝对误差/${\rm{N}} \cdot {\rm{m} }{ {\rm{m} }^{\rm{-2} } }$ 相对误差/% 1 4.57 0.631 67 637.537 55 000 12 637.54 23 2 7.23 0.857 47 310.654 55 000 −7 689.35 −14 3 9.89 0.969 42 029.226 55 000 −12 970.77 −23 4 12.55 1.064 42 413.707 55 000 −12 586.29 −23 5 15.21 1.128 44 421.487 55 000 −10 578.51 −19 6 17.87 1.184 43 510.926 55 000 −11 489.07 −21 7 20.53 1.204 48 689.329 55 000 −6 310.67 −12 8 23.19 1.264 48 374.456 55 000 −6 625.54 −12 9 25.85 1.323 45 863.630 55 000 −9 136.37 −17 10 28.51 1.352 46 408.676 55 000 −8 591.32 −16 11 31.17 1.361 50 609.136 55 000 −4 390.86 −8.0 12 33.83 1.393 51 906.154 55 000 −3 093.85 −5.6 13 36.49 1.414 54 146.236 55 000 −853.76 −1.6 14 39.15 1.493 51 137.903 55 000 −3 862.10 −7.0 15 41.81 1.499 53 009.162 55 000 −1 990.84 −3.6 16 44.47 1.516 54 885.973 55 000 −114.03 −0.2 17 47.13 1.532 53 328.514 55 000 −1 671.486 −3.0 18 49.79 1.579 52 898.537 55 000 −2 101.46 −3.8 19 52.45 1.596 54 391.243 55 000 −608.76 −1.1 20 55.11 1.566 59 828.659 55 000 4 828.66 8.8 21 57.77 1.567 63 030.544 55 000 8 030.54 15 22 60.43 1.606 61 624.111 55 000 6 624.11 12 23 63.09 1.608 63 437.057 55 000 8 437.06 15 24 65.75 1.633 63 426.907 55 000 8 426.91 15 25 68.41 1.647 63 748.102 55 000 8 748.10 16 26 71.07 1.689 61 791.874 55 000 6 791.87 12 
下载: 导出CSV
-
[1] 干福熹, 林凤英, 高文燕. 无机玻璃弹性的测试方法和计算方法[J]. 硅酸盐学报,1978(4):256-265. GAN Fuxi, LIN Fengying, GAO Wenyan. Test method and calculation method of inorganic glass elasticity[J]. Journal of the Chinese Ceramic Society,1978(4):256-265.
[2] 赵宏伟, 董晓龙, 张霖, 等. 块体材料弹性模量的四点弯曲自动测试[J]. 吉林大学学报(工学版),2016,46(1):140-145. ZHAO Hongwei, DONG Xiaolong, ZHANG Lin, et al. Automatic four-point bending test of elastic modulus of block materials[J]. Journal of Jilin University (Engineering and Technology edition),2016,46(1):140-145.
[3] 申志荣, 刘建科. 测定玻璃体弹性模量的一种新方法[J]. 陕西科技大学学报(自然科学版),1999,17(1):116-118. SHEN Zhirong, LIU Jianke. A new method of measuringbulk modulus of glass[J]. Journal of Shaanxi University of Science & Technology,1999,17(1):116-118.
[4] 万德田, 包亦望, 刘小根, 等. 玻璃材料的弹性模量评价技术和影响因素分析[J]. 门窗,2012(8):31-35. WAN Detian, BAO Yiwang, LIU Xiaogen, et al. Evaluation of elastic modulus of glass materials and analysis of influencing factors[J]. Doors & Windows,2012(8):31-35.
[5] 侯俊杰. 深入浅出MFC[M]. 武汉: 华中科技大学出版社, 2002. HOU Junjie. Deep penetration of MFC[M]. Wuhan: Huazhong University of Science and Technology Press, 2002.
[6] 毛星云., 冷雪飞, 王碧辉. OpenCV3编程入门[M]. 北京: 电子工业出版社, 2015: 1-437. MAO Xingyun. Introduction to OpenCV3 programming[M]. Beijing: Electronic Industry Press, 2015: 1-437.
[7] 何东健. 数字图像处理[M]. 西安: 西安电子科技大学出版社, 2008: 1-329. HE Dongjian. Digital image processing[M]. Xi'an:Xi 'an University of Electronic Science and Technology Press, 2008: 1-329.
[8] 杨易, 郭长立, 郭朝霞, 等. 基于图像处理的牛顿环应力测量方法[J]. 实验室研究与探索,2017,36(8):71-76. YANG Yi, GUO Changli, GUO Chaoxia, et al. Newton ring stress measurement method based on image processing[J]. Laboratory Research and Exploration,2017,36(8):71-76.
[9] 郭长立. 大学物理教程[M]. 西安: 西北工业大学出版社, 2009. GUO Changli. University physics course[M]. Xi 'an: Northwest Polytechnical University Press, 2009.
[10] 郭长立, 杨曼, 郭朝霞, 等. 基于光干涉法的光学玻璃应力测量方法研究[J]. 应用光学,2015,36(6):923-930. doi: 10.5768/JAO201536.0603005 GUO Changli, YANG Man, GUO Zhaoxia, et al. Study on optical glass stress measurement method based on optical interference method[J]. Journal of Applied Optics,2015,36(6):923-930. doi: 10.5768/JAO201536.0603005
[11] 王必利. 利用牛顿环实验测定玻璃的弹性模量[J]. 大学物理实验,2012,25(2):44-46. doi: 10.3969/j.issn.1007-2934.2012.02.016 WANG Bili. Determination of elastic modulus of glass by Newton's ring experiment[J]. University Physics Experiment,2012,25(2):44-46. doi: 10.3969/j.issn.1007-2934.2012.02.016
[12] 郭长立, 冯小强. 光学玻璃挠曲刚度的光学测量方法研究[J]. 光学技术,2018,44(2):231-236. GUO Changli, FENG Xiaoqiang. Study on optic measurement method of optical glass flexural rigidity[J]. Optical Technique,2018,44(2):231-236.
[13] 宋淑珍, 张永利, 王教方. 牛顿环中心暗斑大小对测量结果影响的研究[J]. 大学物理实验,2006,19(3):33-35. doi: 10.3969/j.issn.1007-2934.2006.03.010 SONG Shuzhen, ZHANG Yongli, WANG Jiaofang. A study of the effects of the size of the Newton rings center dark spot on the measuring results[J]. Physical Experiment of College,2006,19(3):33-35. doi: 10.3969/j.issn.1007-2934.2006.03.010
[14] 熊有德. 材料强度力学[M]. 北京: 科学出版社, 2009: 218-260. XIONG Youde. Strength mechanics of materials[M]. Beijing: Science Press, 2009: 218-260.
[15] 郭长立, 王守全, 冯小强, 等. 一种牛顿环应力测量装置: 中国, ZL105865686A[P]. 2016-08-17 GUO Changli, WANG Shouquan, FENG Xiaoqiang, et al. A Newton ring stress measuring device: China, ZL105865686A[P]. 2016-08-17.
-
期刊类型引用(5)
1. 方知. 高光谱图像降维压缩比自动检测数学模型仿真. 计算机仿真. 2021(04): 487-491 . 
百度学术
2. 陈敏，仲兆祥. 电感耦合等离子光谱仪的实验室定标方法. 实验室研究与探索. 2020(05): 12-16 . 百度学术
3. 张朝忙，叶远智，邓轶，王建邦. 我国自然资源监测技术装备发展综述. 国土资源遥感. 2020(03): 8-14 . 百度学术
4. 冯姗，曾祥忠. 四通道可见光光谱相机的设计. 应用光学. 2019(03): 393-398 . 本站查看
5. 谢臣瑜，翟文超，李健军，高放，李越，吴浩宇，郑小兵. 基于超连续激光-单色仪的绝对辐亮度响应度定标. 光学学报. 2019(07): 187-193 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 923
- HTML全文浏览量: 234
- PDF下载量: 25
- 被引次数: 7
陕公网安备 61011302001501号 