高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于Tchebichef图像矩的氧化锌太赫兹光谱定量研究

殷贤华 郭超 奉慕霖 贺微

引用本文:
Citation:

基于Tchebichef图像矩的氧化锌太赫兹光谱定量研究

    作者简介: 殷贤华(1974-), 男, 副教授, 主要研究领域为太赫兹检测技术。E-mail:cpxu_ck@163.com.
  • 基金项目:

    广西壮族自治区自然科学基金资助项目 2018GXNSFAA281341

    桂林电子科技大学研究生教育创新计划资助项目 2017YJCX97

    国家自然科学基金资助项目 61841502

    国家自然科学基金资助项目 11574059

  • 中图分类号: O433.4

Quantitative study on terahertz spectra of zinc oxide based on Tchebichef image moments

  • CLC number: O433.4

  • 摘要: 为了提高多组分混合物太赫兹光谱定量分析的准确性和稳定性, 将Tchebichef图像矩结合偏最小二乘回归(PLS)的建模方法引入到太赫兹光谱定量分析中, 以橡胶添加剂氧化锌、白炭黑与丁腈橡胶的混合物为实验样本, 对混合物中的氧化锌进行定量分析。获取样本太赫兹吸光度数据后, 用3维光谱自构建方法构建样本的3维光谱, 利用Tchebichef图像矩提取样本3维光谱灰度图的特征信息, 然后用PLS建立定量模型对氧化锌进行定量分析。结果表明, 由本方法得到的样本预测集相关系数和均方根误差分别为0.9993和0.0351, 与用PLS和支持向量回归直接对吸光度数据进行建模的方法相比, 结果更加准确和稳定。本方法能够对混合物中氧化锌进行准确、稳定的定量分析, 为多组分混合物太赫兹光谱定量分析提供了新思路。
  • Figure 1.  Schematic of terahertz time-domain spectroscopy system

    Figure 2.  a—absorption spectra of three pure samples b—absorption spectra of the mixture samples with different proportions of zinc oxide content

    Figure 3.  a—THz absorbance spectra of samples (one sample) b—3-D spectra of samples c—gray image

    Figure 4.  Sample prediction set results from three methods

    a—PLS b—SVR c—Tchebichef image moments

    Table 1.  Proportional information of components in experimental samples

    number NBR ZnO silica
    1 0.50 0.10 0.40
    2 0.50 0.14 0.36
    3 0.50 0.18 0.32
    4 0.50 0.22 0.28
    5 0.50 0.26 0.24
    6 0.50 0.30 0.20
    7 0.50 0.34 0.16
    8 0.50 0.38 0.12
    9 0.50 0.42 0.08
    10 0.50 0.46 0.04
    下载: 导出CSV

    Table 2.  Results of regression analysis

    method calibration prediction
    Rc ERMS, c Rp ERMS, p
    PLS 0.9854 0.0798 0.9695 0.0921
    SVR 0.9964 0.0384 0.9949 0.0438
    Tchebichef 0.9996 0.0126 0.9993 0.0351
    下载: 导出CSV
  • [1]

    CHEN H. Application of organic zinc instead of indirect zinc oxide in tire tread compound[J].China Rubber, 2016, 32(10):45-48(in Chinese). 
    [2]

    CHEN Zh H. Development of green tire in china and demand for rubber additives[J].Rubber Science and Technology, 2013, 11(08):5-12(in Chinese). 
    [3]

    ZHANG Zh, MENG K, ZHU L G, et al. Study on terahertz absorption characteristics of ischemic rat brain tissue[J].Laser Technology, 2016, 40(3):372-376(in Chinese). 
    [4]

    ZHANG H, LI Zh, HU F, et al. Sensitive distinction between herbs by terahertz spectroscopy and a metamaterial resonator[J].Spectroscopy Letters, 2018, 51(4):174-178. doi: 10.1080/00387010.2018.1442357
    [5]

    HUANG X B. The study on dection methods of solid drug constituents based on terahertz time-domain spectroscopy[D].Hangzhou: Zhejiang University, 2017: 17-19(in Chinese).
    [6]

    MUKUNDAN R, ONG S H, LEE P A. Image analysis by Tchebichef moments [J]. IEEE Transactions on Image Processing, 2001, 10(9):1357-1364. doi: 10.1109/83.941859
    [7]

    LIU J J, LI B Q, WANG X, et al. Applying Tchebichef image moments to the simultaneous quantitative analysis of the four components in corn based on raw NIR spectra[J]. Chemometrics and Intelligent Laboratory Systems, 2018(173):14-20. 
    [8]

    PAN Y, CUI Y Y, WU X J, et al. Krawtchouk moment method for the quantitative analysis of polycyclic aromatic hydrocarbons based on fluorescence three-dimensional spectra[J]. Spectroscopy and Spectral Analysis, 2018, 38(12):3785-3789(in Chinese). 
    [9]

    XUE W, BAO Q, LI H, et al. An efficient approach to the quantitative analysis of humic acid in water [J]. Food Chemistry, 2016, 190:1033-1039. doi: 10.1016/j.foodchem.2015.06.074
    [10]

    YIN X H, JIANG Y, LV B Ch, et al. Quantitative analysis of 2-mercaptobenzothiazole based on terahertz time-domain spectroscopy [J]. Laser Technology, 2019, 43(1):83-87(in Chinese). 
    [11]

    DUVILLARET L, GREET F, COUTAZ J L. A reliable method for extraction of material parameters in terahertz time-domain spectroscopy [J]. IEEE Journal of Sekected Topics in Quantum Electronics, 1996, 2(3):739-746. doi: 10.1109/2944.571775
    [12]

    DUVILLARET L, GREET F, COUTAZ J L. Highly precise determination of optical constants and sample thickness in terahertz time-domain spectroscopy [J]. Applied Optics, 1999, 38(2):409-415. doi: 10.1364/AO.38.000409
    [13]

    QIN B Y, LI Z, CHEN T, et al. Identification of genetically modified cotton seeds by terahertz spectroscopy with MPGA-SVM [J]. Optik, 2017, 142: 576-582. doi: 10.1016/j.ijleo.2017.06.030
    [14]

    MUKUNDAN R, ONG S H, LEE P A. Image analysis by Tchebichef moments[J].IEEE Transactions on Image Processing, 2001, 10(9):1357-1364. doi: 10.1109/83.941859
    [15]

    WANG X. Research on the application of several image moment characteristics in analysis [D]. Lanzhou: Lanzhou University, 2017: 22-23(in Chinese).
    [16]

    YAP P T, PARAMESRAN R, ONG S H. Image analysis by Krawtchouk moments[J].IEEE Transactions on Image Processing, 2003, 12(11):1367-1377. doi: 10.1109/TIP.2003.818019
    [17]

    SEE K W, LOKE K S, LEE P A, et al. Image reconstruction using various discrete orthogonal polynomials in comparison with DCT [J]. Applied Mathematics and Computation, 2007, 193 (2): 346-359. doi: 10.1016/j.amc.2007.03.080
    [18]

    XU M L. The new applications of Tchebichef image moments in analysis. [D]. Lanzhou: Lanzhou University, 2017: 34-35(in Chinese).
    [19]

    LI B, LONG Y, LIU H Sh, et al. The determination of glucose based on terahertz spectroscopy [J]. Spectroscopy and Spectral Analysis, 2017, 37(7): 2165-2170(in Chinese). 
    [20]

    LIU J J. The applications of Tchebichef image moments in analysis of near infrared spectra[D]. Lanzhou: Lanzhou University, 2018: 12-13(in Chinese).
  • [1] 郝晓剑任龙杨彦伟孙永凯 . 煤中激光诱导击穿光谱的碳元素定量分析. 激光技术, 2020, 44(2): 232-236. doi: 10.7510/jgjs.issn.1001-3806.2020.02.017
    [2] 王绍龙王阳恩陈奇陈善俊 . 激光诱导击穿光谱技术定量分析原油金属元素. 激光技术, 2015, 39(1): 104-108. doi: 10.7510/jgjs.issn.1001-3806.2015.01.021
    [3] 殷贤华姜燕吕斌川陈德勇陈涛 . 2-巯基苯并噻唑的太赫兹时域光谱定量研究. 激光技术, 2019, 43(1): 83-87. doi: 10.7510/jgjs.issn.1001-3806.2019.01.017
    [4] 欧阳爱国蔡会周李斌刘刚 . 玉米粉中苯甲酸的太赫兹光谱定量检测研究. 激光技术, 2020, 44(4): 478-484. doi: 10.7510/jgjs.issn.1001-3806.2020.04.014
    [5] 陆俊高淑梅熊婕杨幼益陈国庆 . 女性尿液荧光光谱学特性及机理分析. 激光技术, 2010, 34(1): 45-47,84. doi: 10.3969/j.issn.1001-3806.2010.01.013
    [6] 罗宇恒万恩来刘玉柱 . 利用LIBS技术对电烙铁的烟雾进行在线分析. 激光技术, 2022, 46(5): 663-667. doi: 10.7510/jgjs.issn.1001-3806.2022.05.014
    [7] 李江涛鲁翠萍沙文 . 复合肥中磷元素的激光诱导击穿光谱定量分析. 激光技术, 2019, 43(5): 601-607. doi: 10.7510/jgjs.issn.1001-3806.2019.05.003
    [8] 崔祖文向昱霖张宇璇吴玉兰何宇刘作业孙少华 . 铀酰溶液中铀含量的激光诱导击穿光谱测量. 激光技术, 2021, 45(3): 331-335. doi: 10.7510/jgjs.issn.1001-3806.2021.03.012
    [9] 付奎娄本浊孙彦清龙姝明黄朝军 . Zn0.95-xBe0.05MnxSe稀磁半导体的光谱特性分析. 激光技术, 2015, 39(1): 135-139. doi: 10.7510/jgjs.issn.1001-3806.2015.01.027
    [10] 王龙沈学举张维安董红军 . 高斯光束的光谱传输特性分析. 激光技术, 2012, 36(5): 700-703. doi: 10.3969/j.issn.1001-3806.2012.05.032
    [11] 余洋赵南京王寅方丽孟德硕胡丽马明俊刘建国 . 含铅污泥中Pb的激光诱导击穿光谱定量反演研究. 激光技术, 2015, 39(4): 537-540. doi: 10.7510/jgjs.issn.1001-3806.2015.04.024
    [12] 王召兵刘涛郝殿中彭捍东张霞吴闻迪 . 有色方解石晶体的光谱分析. 激光技术, 2008, 32(6): 596-597,604.
    [13] 徐远泽郭建强高晓蓉王黎王泽男 . 温度对CO光谱线吸收的影响分析. 激光技术, 2010, 34(6): 778-780,784. doi: 10.3969/j.issn.1001-3806.2010.06.016
    [14] 武传龙冯国英韩旭姜海涛欧群飞王建军李密 . 微型光纤光谱仪的波长定标分析. 激光技术, 2012, 36(5): 682-685. doi: 10.3969/j.issn.1001-3806.2012.05.027
    [15] 关世豪杨桄李豪付严宇 . 3维卷积递归神经网络的高光谱图像分类方法. 激光技术, 2020, 44(4): 485-491. doi: 10.7510/jgjs.issn.1001-3806.2020.04.015
    [16] 杨桄田张男李豪关世豪 . 基于线性解混的高光谱图像目标检测研究. 激光技术, 2020, 44(2): 143-147. doi: 10.7510/jgjs.issn.1001-3806.2020.02.001
    [17] 朱潘雨黄敏赵鑫 . 基于SMOTE-UVE-SVM的小麦种子纯度高光谱图像检测. 激光技术, 2024, 48(2): 281-287. doi: 10.7510/jgjs.issn.1001-3806.2024.02.021
    [18] 赵佳乐王广龙周冰应家驹王强辉李秉璇 . 基于边缘剔除的陆基高光谱图像噪声评估方法. 激光技术, 2023, 47(1): 121-126. doi: 10.7510/jgjs.issn.1001-3806.2023.01.019
    [19] 王超王兴田进寿卢裕曹希斌王俊锋徐向晏温文龙 . 高能阿秒脉冲聚焦及光谱分析复合系统设计. 激光技术, 2014, 38(2): 251-254. doi: 10.7510/jgjs.issn.1001-3806.2014.02.022
    [20] 谷玲黄奇峰马超群史院平陈国庆高淑梅 . 液态纯牛奶的3维荧光光谱特性与成分分析. 激光技术, 2011, 35(6): 742-744,769. doi: 10.3969/j.issn.1001-3806.2011.06.005
  • 加载中
图(4) / 表(2)
计量
  • 文章访问数:  5927
  • HTML全文浏览量:  4275
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-21
  • 录用日期:  2019-02-27
  • 刊出日期:  2019-11-25

基于Tchebichef图像矩的氧化锌太赫兹光谱定量研究

    作者简介: 殷贤华(1974-), 男, 副教授, 主要研究领域为太赫兹检测技术。E-mail:cpxu_ck@163.com
  • 1. 桂林电子科技大学 电子工程与自动化学院, 桂林 541004
  • 2. 广西自动检测技术与仪器重点实验室, 桂林 541004
基金项目:  广西壮族自治区自然科学基金资助项目 2018GXNSFAA281341桂林电子科技大学研究生教育创新计划资助项目 2017YJCX97国家自然科学基金资助项目 61841502国家自然科学基金资助项目 11574059

摘要: 为了提高多组分混合物太赫兹光谱定量分析的准确性和稳定性, 将Tchebichef图像矩结合偏最小二乘回归(PLS)的建模方法引入到太赫兹光谱定量分析中, 以橡胶添加剂氧化锌、白炭黑与丁腈橡胶的混合物为实验样本, 对混合物中的氧化锌进行定量分析。获取样本太赫兹吸光度数据后, 用3维光谱自构建方法构建样本的3维光谱, 利用Tchebichef图像矩提取样本3维光谱灰度图的特征信息, 然后用PLS建立定量模型对氧化锌进行定量分析。结果表明, 由本方法得到的样本预测集相关系数和均方根误差分别为0.9993和0.0351, 与用PLS和支持向量回归直接对吸光度数据进行建模的方法相比, 结果更加准确和稳定。本方法能够对混合物中氧化锌进行准确、稳定的定量分析, 为多组分混合物太赫兹光谱定量分析提供了新思路。

English Abstract

    • 氧化锌(ZnO)是轮胎制造过程中重要的添加剂,氧化锌过量释放对人体和环境,尤其是对锌较为敏感的水生动植物会造成不良影响。美国加利福尼亚州提出过相应的法案,建议在轮胎中限制使用氧化锌,同时可能禁售含有一定锌的轮胎产品[1]。轮胎中氧化锌含量的多少将影响轮胎质量是否达到标准[2]。因此,对橡胶添加剂氧化锌含量的检测具有重要意义。

      太赫兹(teranerta, THz)波是频率在0.1THz~10THz之间的电磁辐射[3],太赫兹时域光谱(terahertz time-domain spectroscopy,THz-TDS)技术作为一种新型的光谱分析技术,因其具有强大的穿透能力、安全性和“指纹”特性[4]等独特优点,已经在物质的定量分析中得到应用。在利用太赫兹光谱技术对多组分混合物进行定量分析时,由于混合物中不同物质的波长相近或相互干扰等原因会导致光谱重叠和失真[5],此时采用常规的建模方法(如偏最小二乘回归(partial least squares,PLS)和支持向量回归(support vector regression,SVR))往往不能得到满意的分析结果。

      图像矩作为图像处理的方法已在计算机视觉等领域得到广泛应用。近年来,因其具有多分辨率及固有不变性等优点[6],图像矩已经应用于利用光谱技术对物质进行定量或定性分析中。LIU等人[7]将Tchebichef图像矩应用于利用红外光谱对混合物的定量分析中,解决了红外光谱由于光谱重叠和移位等原因造成的定量分析结果不准确的问题。PAN等人[8]提出了一种将Krawtchouk图像矩和荧光光谱、广义回归神经网络相结合的定量分析方法,并对多环芳烃进行了定量分析,得到了精确的分析结果。XUE等人[9]在用3维荧光光谱对水中腐殖酸含量进行定量分析时,利用Zernike矩对3维荧光光谱灰度图提取特征值,建立了腐殖酸的定量模型,得到了比多维偏最小二乘和交替三线性分解方法更加可靠和准确的结果。

      针对多组分混合物太赫兹光谱定量分析中由于光谱重叠和失真等原因而容易造成分析结果不准确的问题,本文中拟在太赫兹光谱定量分析中引入Tchebichef图像矩和PLS建模方法,对橡胶添加剂氧化锌、白炭黑和丁腈橡胶(nitrile butadiene rubber, NBR)组成的混合物中氧化锌的含量进行预测,以提高氧化锌定量分析的准确性和稳定性。

    • 丁腈橡胶在汽车行业中应用广泛,白炭黑是制造“绿色轮胎”的重要添加剂,将氧化锌与上述两种材料按照不同比例进行混合,制作实验样本。实验中所用丁腈橡胶的质量分数为0.9375,氧化锌和白炭黑的质量分数均大于0.995。实验样本采用压片法进行制备,首先按照设计好的样本配置比例进行称量,然后将样本在研钵中进行充分的研磨混合,最后利用压片机在666Pa的压力下压制成片。压制成的每个样片厚度1mm左右,直径12mm左右,重量为180mg。在实际的同种规格橡胶轮胎中,橡胶所占的比例往往是固定的,为了使实验样本比较接近实际,把样片中丁腈橡胶所占质量分数固定为0.50,剩余0.50则由氧化锌和白炭黑按比例构成,基于实际问题和本机构实验仪器精度,将混合物实验样本中氧化锌最低质量分数设为0.1,样片中氧化锌质量分数以0.04梯度增加,共设置10种不同质量分数,每个样片分别压制6个,共60个样本。样本中各成分的质量分数详见表 1

      Table 1.  Proportional information of components in experimental samples

      number NBR ZnO silica
      1 0.50 0.10 0.40
      2 0.50 0.14 0.36
      3 0.50 0.18 0.32
      4 0.50 0.22 0.28
      5 0.50 0.26 0.24
      6 0.50 0.30 0.20
      7 0.50 0.34 0.16
      8 0.50 0.38 0.12
      9 0.50 0.42 0.08
      10 0.50 0.46 0.04
    • 实验仪器采用美国Zomega公司研制的Z-3太赫兹时域光谱系统,系统频谱分辨率为5GHz,信噪比大于70dB,激光器采用波长为780nm、脉宽小于100fs的超快飞秒激光器,太赫兹时域光谱系统原理如图 1所示。由于水分会对THz波造成吸收[10],为减少水分对实验的影响,实验时使用除湿机对实验环境进行除湿干燥,使实验环境中H2O的质量分数控制在0.01以下。实验样本在测量之前需要在真空恒温箱干燥2h左右,以减少实验样本中水分对实验结果的影响,实验样本干燥完成后,用太赫兹时域光谱系统对氧化锌、白炭黑和丁腈橡胶的3种纯样本以及制备的60个混合物样本分别进行测量,为减少实验误差,对每个样本测量3次并取平均以获得样本的最终光谱信号数据,实验中还需对无样本时的环境进行测量以作为参考信号。

      Figure 1.  Schematic of terahertz time-domain spectroscopy system

    • 在得到样本的光谱数据后,需对数据进行处理和提取光学参量,本研究中采用DUVILLARET等人提出的模型[11-12]作为数据处理和光学参量的提取方法,对由实验中获得的参考信号和样本信号进行快速傅里叶变换,得到频域信号,如下式所示:

      $ \frac{E_{\mathrm{s}}(\omega)}{E_{\mathrm{r}}(\omega)}=A(\omega) \exp [\mathrm{i} \varphi(\omega)] $

      (1)

      式中,ω为太赫兹波震动的角频率,Es(ω)为样本信号,Er(ω)为参考信号,A(ω)为振幅,φ(ω)为相位差。

      由于对实验样本厚度的测量会存在一定的误差,为避免样本厚度测量误差对后续定量分析结果的影响,本研究中采用无量纲的吸光度B作为定量分析模型建模参量。吸光度表示光波被材料吸收的程度,其计算公式如下[13]

      $ B=-\lg \left[\frac{E_{{\rm s}}(\omega)^{2}}{E_{\mathrm{r}}(\omega)^{2}}\right] $

      (2)
    • 2001年,MUKUNDAN提出基于离散Tchebichef多项式的离散正交矩函数[14],除了具有传统图像矩的平移、旋转、缩放不变性和正交矩冗余量小等优势外,由于Tchebichef图像矩是建立在图像坐标空间的离散域上的,不包括任何数值上的近似,也不需要坐标空间的转化,因而性能更优[15]。其具体的定义及相关计算如下[14, 16]

      对于一个给定的N×M大小的图像,其在阶数为(m+n)时的Tchebichef图像矩定义为:

      $ \begin{array}{c} {T_{n, m}} = \frac{1}{{\rho (n, N)\rho (m, M)}}\sum\limits_{x = 0}^{N - 1} {\sum\limits_{y = 0}^{M - 1} {{t_n}} } (x){t_m}(y)f(x, y)\\ (n = 0, 1, 2, \cdots , N - 1;m = 0, 1, 2, \cdots , M - 1) \end{array} $

      (3)

      式中, tn(y)和tm(y)是阶数为n阶和阶数为m阶对应的标准正交多项式,ρ(n, N)和ρ(m, M)是标准范数的平方,f(x, y)是图像中的坐标为(x, y)处的图像强度。

      为了能够方便计算Tchebichef多项式,使用如下的递推关系式[17]

      $ \begin{array}{c} t_{n}(x) =\frac{(2 n-1)(2 x-N+1) t_{n-1}(x)}{n \times N}-\\ \frac{(n-1)\left[N^{2}-(n-1)^{2}\right] t_{n-2}(x)}{n \times N^{2}} \end{array} $

      (4)

      $ t_{0}(x)=1 $

      (5)

      $ t_{1}(x)=\frac{2 x+1-N}{N} $

      (6)

      $ \rho(n)=\left(\frac{2 n-1}{2 n+1}\right)\left(1-\frac{n^{2}}{N^{2}}\right) \rho(n-1) $

      (7)

      $ \rho(0)=N $

      (8)

      Tchebichef图像矩的阶数nm决定着提取的图像特征精细程度。低阶矩时,矩代表的是趋势信息,而高阶矩时,矩代表的是细节信息[18]

    • 研究中所建立的定量分析模型以相关系数R和均方根(root mean square, RMS)误差ERMS作为评价指标。相关系数表示样本的预测值和参考值之间的线性相关程度,均方根误差反映样本的预测值和参考值的离散程度,相关系数的值越接近1,均方根误差的值越小,说明所建立的模型性能越好[19]。其计算公式如下所示:

      $ R=\sqrt{1-\frac{\sum\limits_{i=1}^{n}\left(y_{i}-\hat{y}_{i}\right)^{2}}{\sum\limits_{i=1}^{n}\left(y_{i}-\bar{y}\right)^{2}} \times 100 \%} $

      (9)

      $ {E_{{\rm{RMS}}}} = \sqrt {\frac{1}{n}\sum\limits_{i = 1}^n {{{\left( {{y_i} - {{\hat y}_i}} \right)}^2}} } $

      (10)

      式中,n为样本总数,yi为第i个样本的参考值,$ \hat{y}_{i}$为第i个样本的预测值,yn个参考值的平均值。

    • 本研究利用太赫兹时域光谱系统分别对ZnO、白炭黑(NBR)和丁腈橡胶(silica)的纯样本和3种物质的混合物样本进行了测量,通过计算得到了3种纯样本的吸光度谱, 如图 2a所示。为了方便观察混合物样本质量分数和吸光度之间的关系,得到10种氧化锌质量分数所占比例不同的混合物样本的吸光度谱, 如图 2b所示。

      Figure 2.  a—absorption spectra of three pure samples b—absorption spectra of the mixture samples with different proportions of zinc oxide content

      由于不同分子具有其独特的震动模式,特定分子在太赫兹波段内具有其特有的吸收峰。从图 2a中可以看出,氧化锌在1.38THz处有明显的吸收峰,白炭黑在1.33THz和1.53THz处有明显吸收峰,丁腈橡胶在1.51THz处有明显吸收峰;从图 2b中可以看出,由氧化锌、白炭黑和丁腈橡胶按不同比例混合而成的混合物在1.26THz和1.48THz处有明显吸收峰,并且随着氧化锌含量的增加,吸光度也随之变化,说明混合物中氧化锌的含量和吸光度存在相互关系,初步说明了利用太赫兹光谱对混合物中氧化锌含量进行定量分析的可行性。同样也可以看出,不同样本的光谱之间存在重叠。

    • 由实验得到的样本的太赫兹吸光度谱是2-D光谱图(见图 3a),为了能够用Tchebichef图像矩提取其特征信息,需要构建样本的3-D光谱图,本文中采用3-D光谱的自构建方法[20]构建太赫兹吸光度图谱的3-D光谱。3-D光谱的自构建方法是对2-D光谱进行复制来构建样本的3-D光谱,该方法得到的3-D光谱只包含样本信息,并没有增加其它信息,从而为Tchebichef图像矩方法提取光谱的特征信息提供了方便[20]。在本文中,对每一个样本的吸光度谱复制两次,共得到3条吸光度谱,以此来构建样本的3-D光谱图,如图 3b所示。得到的3-D光谱图是一个大小为3×98(太赫兹频率点)的矩阵,此矩阵可以表示为3-D光谱的灰度图,如图 3c所示。此时,便可以利用Tchebichef图像矩提取其特征信息。

      Figure 3.  a—THz absorbance spectra of samples (one sample) b—3-D spectra of samples c—gray image

    • 为了能够准确地对混合物中橡胶添加剂氧化锌的含量进行预测,在得到样本3-D光谱图的灰度图后,利用Tchebichef图像矩对其提取特征信息,将提取的特征信息采用PLS建立定量分析模型。在10种质量分数的混合物样本中每种质量分数选取4个样片,共40个样本作为校正集用来建立定量模型,剩余的20个样本作为预测集。为了比较本方法的优越性,作者还采用PLS和SVR直接对得到的吸光度数据建立定量模型,对氧化锌进行定量分析。3种定量模型得到的样本预测集的结果如图 4所示,其相关系数和均方根误差如表 2所示。

      Figure 4.  Sample prediction set results from three methods

      Table 2.  Results of regression analysis

      method calibration prediction
      Rc ERMS, c Rp ERMS, p
      PLS 0.9854 0.0798 0.9695 0.0921
      SVR 0.9964 0.0384 0.9949 0.0438
      Tchebichef 0.9996 0.0126 0.9993 0.0351

      图 4中的斜线表示零误差线,由定量模型得到的预测值越接近零误差线, 表明预测结果越接近真实值。从图 4中可以看出,用Tchebichef图像矩结合PLS方法所建立的定量模型得到的氧化锌的含量预测值比利用SVR和PLS所建立的定量模型得到的预测值更接近零误差线,且波动相对较小,说明预测结果更加准确和稳定。

      表 2中可以看出,由PLS所建立的定量模型的校正集的相关系数Rc和均方根误差ERMS, c分别为0.9854和0.0798,预测集的相关系数Rp和均方根误差ERMS, p分别为0.9695和0.0921;SVR所建立的定量模型的校正集的相关系数Rc和其均方根误差ERMS, c分别为0.9964和0.0384,预测集的相关系数Rp和其均方根误差ERMS, p分别为0.9949和0.0438;而由Tchebichef图像矩和PLS所建立的定量模型的校正集和预测集的相关系数RcRp均大于0.9993,校正集和预测集的均方根误差ERMS, cERMS, p均小于0.0351,对比可见,本模型性能更优。所以,本方法对混合物中橡胶添加剂氧化锌的定量分析更有效。

    • 本文中以橡胶添加剂氧化锌、白炭黑和丁腈橡胶的混合物为实验对象,利用3-D光谱的自构建方法构建样本的3-D光谱,将Tchebichef图像矩与PLS建模方法相结合建立了定量分析模型,实现了对多组分混合物中橡胶添加剂氧化锌含量的预测。由定量模型得到的样本预测集相关系数和均方根误差分别为0.9993和0.0351,均优于SVR和PLS对样本吸光度数据所建立的定量模型。以上研究为太赫兹光谱技术用于多组分混合物的定量分析提供了新的参考方法。

参考文献 (20)

目录

    /

    返回文章
    返回