CN114112979A - 一种基于太赫兹光谱的农产品中农药残留定量检测方法 - Google Patents
一种基于太赫兹光谱的农产品中农药残留定量检测方法 Download PDFInfo
- Publication number
- CN114112979A CN114112979A CN202111463994.7A CN202111463994A CN114112979A CN 114112979 A CN114112979 A CN 114112979A CN 202111463994 A CN202111463994 A CN 202111463994A CN 114112979 A CN114112979 A CN 114112979A
- Authority
- CN
- China
- Prior art keywords
- frequency band
- terahertz
- spectrum
- absorption spectrum
- agricultural product
- Prior art date
- Legal status (The legal status 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 status listed.)
- Pending
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 83
- 239000000447 pesticide residue Substances 0.000 title claims abstract description 71
- 238000001514 detection method Methods 0.000 title claims abstract description 55
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 34
- 239000000575 pesticide Substances 0.000 claims abstract description 29
- 238000012937 correction Methods 0.000 claims abstract description 15
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 238000013528 artificial neural network Methods 0.000 claims abstract description 7
- 238000005457 optimization Methods 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 12
- NWBJYWHLCVSVIJ-UHFFFAOYSA-N N-benzyladenine Chemical compound N=1C=NC=2NC=NC=2C=1NCC1=CC=CC=C1 NWBJYWHLCVSVIJ-UHFFFAOYSA-N 0.000 claims description 9
- 239000005906 Imidacloprid Substances 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- YWTYJOPNNQFBPC-UHFFFAOYSA-N imidacloprid Chemical compound [O-][N+](=O)\N=C1/NCCN1CC1=CC=C(Cl)N=C1 YWTYJOPNNQFBPC-UHFFFAOYSA-N 0.000 claims description 8
- 229940056881 imidacloprid Drugs 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 239000003905 agrochemical Substances 0.000 claims description 5
- 238000002156 mixing Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 238000010521 absorption reaction Methods 0.000 abstract description 28
- 230000003595 spectral effect Effects 0.000 abstract description 14
- 238000013507 mapping Methods 0.000 abstract description 2
- 239000002245 particle Substances 0.000 description 18
- 239000000203 mixture Substances 0.000 description 14
- 241000209140 Triticum Species 0.000 description 5
- 235000021307 Triticum Nutrition 0.000 description 5
- 235000013312 flour Nutrition 0.000 description 5
- KHZWIIFEFQBNKL-UHFFFAOYSA-N 2-(2,6-dichlorophenoxy)acetic acid Chemical compound OC(=O)COC1=C(Cl)C=CC=C1Cl KHZWIIFEFQBNKL-UHFFFAOYSA-N 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004445 quantitative analysis Methods 0.000 description 4
- YOYAIZYFCNQIRF-UHFFFAOYSA-N 2,6-dichlorobenzonitrile Chemical compound ClC1=CC=CC(Cl)=C1C#N YOYAIZYFCNQIRF-UHFFFAOYSA-N 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010238 partial least squares regression Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Toxicology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于太赫兹光谱的农产品中农药残留定量检测方法。本发明利用BP神经网络(BPNN)对含农药农产品的太赫兹光谱数据进行建模，建立农药残留含量与太赫兹吸收光谱间定量模型，即BPNN模型；并结合光谱去噪、基线校正等光谱预处理方法提升光谱质量；同时根据光谱吸收峰特性，进行光谱波长选取，选取高信噪比频谱区间；最后对BPNN典型参数进行优化，完成完整的网络映射，实现农产品中农药残留的高精度检测。从实验结果看，本发明切实可行且高效，可为农产品的安全检测奠定基础。
Description
技术领域
本发明属于食品安全检测技术领域，具体涉及一种基于太赫兹光谱的农产品中农药残留定量检测方法。
背景技术
近些年，我国农产品生产规模正在逐年扩大，产生了巨大的经济效益和社会效益，但农产品生产面临着很多问题，其中最重要的就是农产品中农药残留超标问题屡禁不止。
与传统辐射源相比，太赫兹辐射具有强穿透性、瞬态性、宽带性、相干性、安全性、低能性等诸多特性，太赫兹光谱技术在农产品检测领域展现出巨大的应用潜力，结合相关机器学习算法，可实现农产品中农药残留含量的准确鉴定。目前，已有众多学者开展了基于太赫兹光谱的农产品中农药残留的定量检测研究。但是，目前研究主要具有如下两个问题：
(1)、基于太赫兹光谱的农产品中多种农药残留定量检测中，由于残留成分复杂，吸收谱线重叠混乱，导致光谱数据与农药残留浓度线性不可分，针对此特点的数据无法获得有效的定量回归模型，实现基于太赫兹光谱的农药残留含量的准确检测；
(2)、已报道的基于太赫兹光谱的农产品中多种农药残留定量检测中，农产品中农药残留检测浓度偏高，缺乏实际应用价值。
上述问题限制了太赫兹光谱在农药残留检测中的研究与应用。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足提供一种基于太赫兹光谱的农产品中农药残留定量检测方法，基于太赫兹光谱的农产品中农药残留定量检测方法可用于实现农药残留含量快速、准确、无损检测，具有广阔的应用潜力。
为实现上述技术目的，本发明采取的技术方案为：
一种基于太赫兹光谱的农产品中农药残留定量检测方法，包括：
步骤1：将农产品与已知品种的农药混合，且获取混合后的农产品中农药的残留浓度实际值，制作多组不同残留浓度的样本，采用太赫兹光谱仪对所有样本进行检测，得到太赫兹吸收光谱数据集；
步骤2：对步骤1中检测得到的吸收光谱数据进行光谱预处理；
步骤3：基于预处理后的吸收光谱数据，进行BP神经网络非线性回归模型的构建，其中模型的输入为步骤2预处理后的吸收光谱数据，模型的输出为农药的残留浓度，构建完成后，得到BPNN模型；
步骤4：对步骤3构建的BPNN模型进行参数优化；
步骤5：采用PSO算法从步骤2预处理后的吸收光谱数据中选取多个频段，将每个频段的吸收光谱数据分别输入到参数优化后的BPNN模型，预测得到每个频段对应的农药残留浓度，选择其中预测精度最高的频段，该频段记为最佳频段；
步骤6：选取与步骤1的农产品中具有相同农药品种的待测农产品，采用太赫兹光谱仪对该待测农产品进行检测，得到太赫兹吸收光谱数据，对该太赫兹吸收光谱数据进行光谱预处理，得到预处理后的吸收光谱数据，从预处理后吸收光谱数据中选取与步骤5中最佳频段相同的频段区间，将该频段区间的吸收光谱数据输入到参数优化后的BPNN模型中进行检测，得到农产品的农药残留浓度检测结果。
作为本发明进一步改进的技术方案，所述的光谱预处理具体包括：
依次采用光谱去噪法和基线校正法对太赫兹吸收光谱数据进行光谱预处理，所述光谱去噪法为小波软阈值去噪法，所述基线校正法为二阶导基线校正法。
作为本发明进一步改进的技术方案，所述步骤3中BPNN模型的构建如下：
P＝WI；
即:
其中，W表示神经网络权重矩阵，I表示光谱输入矩阵，P表示预测浓度矩阵；wmv表示第m层网络第v个神经元的权重系数，inputn表示第n个样本的太赫兹吸收光谱数据。
作为本发明进一步改进的技术方案，所述步骤4具体为：
通过误差反向传递法进行参数优化，其中误差表达式为：
ek＝yk-Ok；
其中，yk表示样本的检测浓度，Ok表示样本的预测浓度，ek表示样本的预测误差。
作为本发明进一步改进的技术方案，所述步骤1中的农产品含有的农药品种有1种或多种，所述步骤1中某一农产品包含的农药品种与步骤6中待测农产品包含的农药品种相同。
作为本发明进一步改进的技术方案，所述农产品中含有的已知品种的农药包含6-苄氨基嘌呤、吡虫啉和6-二氯苯腈。
作为本发明进一步改进的技术方案，所述的步骤5具体为：
在0-0.50THz峰宽范围内，以0.02THz为峰宽递增间隔，并选用PSO算法计算不同峰宽下频段起始位置，继而确定不同峰宽对应的多个频段，将每个频段的吸收光谱数据分别输入到参数优化后的BPNN模型，预测得到每个频段的吸收光谱数据对应的农药残留浓度，选择其中预测精度最高的频段，该频段记为最佳频段。
作为本发明进一步改进的技术方案，若样本中农产品含有的农药品种有1种，计算每个频段的吸收光谱数据对应的农药残留浓度预测值的相关系数或者均方根误差，选择其中相关系数最大或者均方根误差最小的频段，该频段为最佳频段；
若样本中农产品含有的农药品种有多种，计算单个频段的吸收光谱数据对应的所有农药残留浓度预测值的相关系数的平均值或者均方根误差的平均值，选择其中相关系数的平均值最大或者均方根误差的平均值最小的频段，该频段为最佳频段。
本发明的有益效果为：
(1)、本发明实现了多种农药残留成分的低浓度定量检测任务，通过光谱去噪、基线校正，解决由于农药残留成分复杂，造成的吸收谱线重叠混乱、光谱质量不佳的问题。
(2)、本发明通过BPNN的构建，解决光谱数据与农药残留浓度线性不可分的问题，且通过网络的调参，进一步提升网络映射的完整性，提高农药残留含量定量检测的精度。
(3)、本发明通过太赫兹光谱的波长选取，遴选出高质量的最佳频段区间，选取对应区间内的吸收光谱数据以定量分析，可获得更理想的检测效果。
附图说明
图1为农产品小麦粉中三种农药(6-BA、Imidacloprid和2，6-D)的制样浓度示意图。
图2为使用太赫兹光谱仪进行样本检测，三种农药的THz光谱检测结果示意图。
图3为不同浓度农药残留混合物的太赫兹吸收系数谱图。
图4为光谱小波去噪后的农药残留混合物的太赫兹吸收系数谱图。
图5为光谱小波去噪及基线校正后的农药残留混合物的太赫兹吸收系数谱图。
图6中(a)为使用MW-BPNN得到0.35THz峰宽下的理想吸收频段及对应检测结果示意图。
图6中(b)为使用MW-BPNN得到0.30THz峰宽下的理想吸收频段及对应检测结果示意图。
图6中(c)为使用MW-BPNN得到0.25THz峰宽下的理想吸收频段及对应检测结果示意图。
图6中(d)为使用MW-BPNN得到0.20THz峰宽下的理想吸收频段及对应检测结果示意图。
图7中(a)为使用PSO对混合物光谱进行波长选取，并给出了不同峰宽对应的6-BA定量检测结果示意图。
图7中(b)为使用PSO对混合物光谱进行波长选取，并给出了不同峰宽的Imidacloprid定量检测结果示意图。
图7中(c)为使用PSO对混合物光谱进行波长选取，并给出了不同峰宽的2,6-D定量检测结果示意图。
图7中(d)为使用PSO对混合物光谱进行波长选取，并给出了不同峰宽的三种农药平均定量检测结果示意图。
具体实施方式
下面根据附图对本发明的具体实施方式作出进一步说明：
本实施例中提供一种基于太赫兹光谱的农产品中农药残留定量检测方法，具体提供了一种基于太赫兹光谱的农产品中三种农药残留物，6-苄氨基嘌呤(6-BA)、吡虫啉(Imidacloprid)和2，6-二氯苯腈(2,6-D)的定量检测方法，图1为农产品小麦粉中三种农药(6-BA、Imidacloprid和2，6-D)的制样浓度，共28组不同农药残留含量浓度样本，农药浓度间隔为0.5mg，残留含量共计3mg。检测主要包括以下步骤：
步骤(1)：采用太赫兹光谱仪进行样本检测，构建农产品与农药的混合物光谱数据集：
测试得到三种农药在0.5-2.5THz频段内的吸收系数谱，如图2所示。6-苄氨基嘌呤在2.08THz处有一明显吸收峰；吡虫啉存在五个吸收峰，分别位于1.20,1.45,1.89,2.05及2.41THz；敌草腈存在3个吸收峰，分别位于1.58,1.98及2.32THz处。
在此基础上，测试得到了小麦粉中不同农药残留含量混合物光谱，如图3所示，包含28组样本的吸收光谱数据，由于三种农药均在2.0THz左右处存在一个吸收峰，其中，6-苄氨基嘌呤在2.08THz处有吸收峰，吡虫啉在2.05THz处有吸收峰，敌草腈在1.98THz处有吸收峰，图中可以清楚地观察到混合物光谱在2.0THz处展现了一个明显的吸收峰，其他位置由于样品浓度过低，且浓度间隔较少，无法观察到明显的吸收峰，且谱线出现交叉重叠，不易区分，在后续定量分析中，光谱交叉重叠会严重影响定量分析的精度。
步骤(2)：根据光谱特征进行相关光谱预处理(光谱去噪、基线校正)：
将步骤(1)中检测得到的吸收光谱数据进行相关光谱预处理，改善光谱信号质量，光谱去噪选取小波软阈值去噪法：
f(t)＝s(t)+e(t)；
WTf(y,x)＝WTs(y,x)+WTe(y,x)；
其中，f(t)表示样本信号；s(t)表示有效信号，e(t)表示噪声信号；x，y分别为原始函数中自变量及因变量，对应实例中农药残留浓度及太赫兹光谱吸收系数；WTf(y,x)、WTs(y,x)和WTe(y,x)分别表示样本信号、有效信号和噪声信号的小波变换函数，并选取自适应阈值进行软阈值去噪。
本实例中具体表现为：对输入的含噪光谱信号进行多层小波分解，并采用小波软阈值对各层小波的高频系数进行收缩处理，最后对处理后的小波系数进行小波反变换，得到去噪后的光谱信号。如图4所示为小麦粉中不同农药残留含量混合物光谱去噪后数据，可以看出去噪后光谱质量得到明显改善，特别地，在吸收峰2.0THz处附近，交叉重叠现象较原始光谱有所优化。
基线校正选取二阶导基线校正法：
其中，yk,t1表示x＝k处光谱信号的一阶导数，yk+△，yk-△表示与k距离为△的信号强度；yk,t2表示x＝k处光谱信号的二阶导数，yk+△，yk-△表示与k距离为△的信号强度，yk表示k处信号强度。
本实例中具体表现为：在频段范围内对各离散光谱点进行二阶导数求解，将二阶导数值为零的点进行记录，并以此点为中心实现基线校正，得到基线校正后的光谱信号。如图5所示为小麦粉中不同农药残留含量混合物光谱去噪并基线校正后数据，与未校正前相比，在吸收峰处，光谱信号质量得到较大改善，不同浓度光谱在吸收峰2.0THz附近区分度更加明显。
步骤(3)：基于优化后的农产品与农药的混合物光谱数据，进行非线性回归模型BP神经网络(BPNN)的构建。
BPNN的构建如下：
P＝WI；
即选取步骤(2)预处理后光谱数据为模型输入，农药残留浓度为输出，进行模型构建。
其中，W表示神经网络权重矩阵，I表示光谱输入矩阵，P表示预测浓度矩阵；wmv表示第m层网络第v个神经元的权重系数，inputn表示第n个样本的太赫兹吸收光谱数据。
步骤(4)：对步骤(3)构建的BPNN模型进行参数优化。
通过不断优化网络误差，反向传递误差实现参数优化：
ek＝yk-Ok；
其中，yk表示样本的检测浓度，Ok表示样本的预测浓度，ek表示样本的预测误差。
步骤(5)：结合移动窗口BPNN(MW-BPNN)及粒子群算法(PSO)，对混合物光谱进行波长选取，遴选出高质量的最佳频段光谱区间。
所述步骤5中MW-BPNN模型中信号窗口宽度的构建如下：
其中，Inputmw表示选定峰宽的光谱输入信号矩阵，xn,k表示第n个样本的第k个光谱数据点，p表示信号峰宽。
本实例中具体表现为：在1.00-2.30THz频段内，分别以0.20THz、0.25THz、0.30THz及0.35THz为对应峰宽、以0.02THz为移动间隔(如以0.20THz为对应峰宽，依次取1.0-1.20THz、1.02-1.22THz、1.04-1.24THz…频段；如以0.25THz为对应峰宽，依次取1.0-1.25THz、1.02-1.27THz、1.04-1.29THz…频段；如以0.30THz为对应峰宽，依次取1.0-1.30THz、1.02-1.32THz、1.04-1.34THz…频段；如以0.35THz为对应峰宽，依次取1.0-1.35THz、1.02-1.37THz、1.04-1.39THz…频段)，选取光谱输入信号矩阵，建立MW-BPNN模型进行农药残留含量定量分析，即将选取的每个频段作为光谱输入信号矩阵输入至参数优化后的BPNN内，得到农药残留成分定量检测结果，选取每个峰宽对应的检测精度高的频段。如图6所示，即四种特定峰宽表征下，农产品中农药残留成分定量检测结果。其中图6中(a)、(b)、(c)、(d)中的灰色斜线填充部分分别为对应峰宽下的检测精度高的频段。结果表明，不同峰宽表征下，定量检测精度最优区间均为包含吸收峰2.0THz的区间，表明了吸收峰在物质吸收光谱中的特殊性。其中Rp-avg为灰色斜线填充部分频段对应的农药残留浓度预测值(检测值)的相关系数平均值，RMSEp-avg为灰色斜线填充部分频段对应的农药残留浓度预测值(检测值)的均方根误差平均值。
所述步骤5中PSO的构建如下：适应度函数公式如下：
其中，f(a)表示适应度函数，Dmmax表示输入样本与随机样本的最大聚集距离，Dmmean表示平均聚集距离。
粒子群中每个粒子具有两个属性：速度和位置，速度代表移动的快慢，位置代表移动的方向。每个粒子在搜索空间中单独的搜寻最优解，粒子的速度及位置更新如下：
vi＝wvi+c1×rand()×(pbesti-xi)；
xi＝xi+vi；
其中，i＝1,2,3…n，i为当前粒子位置，n为粒子总数，vi为粒子运动记忆项，表示当前粒子运动速度大小和方向，pbesti表示设定的跟踪极值，表征当前粒子指向粒子自身最优点的一个矢量，xi表示粒子当前的位置，两者做差求解粒子下一步运动的速度及方向，其中，rand()为0至1之间的随机数，c1为学习因子,w为惯性权重；
惯性权重影响粒子速度及搜索精度，为两者间达到平衡，定义动态惯性权重w如下：
其中，wmax为最大惯性权重，wmax＝0.8，wmin为最小惯性权重，wmin＝0.4，run为当前迭代次数，runmax为最大迭代次数；
设置最大迭代次数300及优化目标相关系数R>0.97，满足其一，则PSO运行结束，获取最佳惯性权重w，并确定频段选取起始位置，否则更新粒子速度及位置，循环计算。其中run每循环一次，更新一次最佳惯性权重w，w替代到上式粒子位置和速度公式中，更新一次粒子位置和粒子速度，对应一个频段选取的起始点(吸收峰宽是确定的)，进而确定频段区间。
本实例中具体表现为：在0-0.50THz峰宽范围内，以0.02THz为峰宽递增间隔(即在0-0.50THz峰宽范围内，分别选取0.02THz、0.04THz、0.06THz、0.08THz、0.10THz…为峰宽，以上峰宽均以0.02THz为移动间隔)，并选用PSO算法计算不同峰宽下多个频段起始位置，继而确定不同峰宽对应的多个频段区间，将每个频段区间内的吸收光谱数据分别输入到参数优化后的BPNN模型，预测得到每个频段区间内的吸收光谱数据对应的农药残留浓度，选择其中预测精度最高的频段，该频段记为最佳频段，从而实现不同吸收谱频段表征下农药残留含量检测精度的讨论。图7给出了不同峰宽对应的每种农药残留浓度预测值的相关系数以及三种农药残留浓度预测值的相关系数平均值，其中，图7中的(d)中圆点表示理想定量检测精度的对应峰宽，可以观察到，当吸收峰宽选取0.20-0.25THz区间时，定量检测精度较高，表明该峰宽可较完整地表征光谱，同时验证了PSO算法在峰宽选取中的适用性。通过PSO选取一个峰宽对应的多个频段，从多个频段对应的农药残留浓度预测值中选取其中最大相关系数的值作为该峰宽对应的最优相关系数，然后再通过峰宽和对应的最优相关系数绘制图7中(a)、(b)、(c)，图7中(d)的纵坐标为图7中(a)、(b)、(c)的相关系数的平均。图7中纵坐标的相关系数对应的频段区间均在吸收峰附近。
步骤(6)：选用参数优化后BPNN模型对预处理及波长选取后光谱进行检测，并对样本数据采用留一交叉验证法(LOOCV)进行划分，得到预测集的定量回归系数(R)及均方根误差(RMSE)。
相关系数R、均方根误差RMSE的计算如下：
其中，n表示检测样本数，yi表示第i个样本对应的农药残留浓度的实际值，
表示所有样本对应的农药残留浓度的平均值，y′i表示第i个样本对应的农药残留浓度的预测值。
下面结合实验数据，对本案的发明效果给出说明。
表1为仅进行参数优化后BPNN农药残留物的定量检测准确率。
表1、参数优化后BPNN农药残留物的定量检测准确率：
表2为进行参数优化及波长选取后农药残留物的定量检测准确率，同时给出了两种传统定量检测算法(偏最小二乘回归PLS、支持向量机回归SVR)的定量检测精度，并与本案进行对比。
表2、参数优化及波长选取后农药残留物的定量检测准确率：
从表1和表2可以看出，对光谱进行相关预处理，可有效提升光谱质量；同时，对定量模型进行合理调参、光谱波长选取，可进一步提升检测精度。此外，与PLS、SVR对比，本案中BPNN具有较强的非线性拟合能力，优化后的网络可以更完整地映射太赫兹光谱与残留浓度间的对应关系。
此外，对于未知农药残留浓度的待测农产品(训练集内最低已知残留浓度<农药残留浓度<训练集内最高已知残留浓度)，可参考实例实施以下步骤，实现农药残留含量检测：
选取与步骤(1)的农产品中具有相同农药品种的待测农产品，采用太赫兹光谱仪对该待测农产品进行检测，得到太赫兹吸收光谱数据，对该太赫兹吸收光谱数据选用小波阈值去噪及直接差分二阶求导进行光谱预处理，得到预处理后的吸收光谱数据，从预处理后吸收光谱数据中选取与步骤(5)中最佳频段相同的频段区间，将该频段区间的吸收光谱数据输入到参数优化后的BPNN模型中进行检测，得到农产品的农药残留浓度检测结果。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，包括：
步骤1：将农产品与已知品种的农药混合，且获取混合后的农产品中农药的残留浓度实际值，制作多组不同残留浓度的样本，采用太赫兹光谱仪对所有样本进行检测，得到太赫兹吸收光谱数据集；
步骤2：对步骤1中检测得到的吸收光谱数据进行光谱预处理；
步骤3：基于预处理后的吸收光谱数据，进行BP神经网络非线性回归模型的构建，其中模型的输入为步骤2预处理后的吸收光谱数据，模型的输出为农药的残留浓度，构建完成后，得到BPNN模型；
步骤4：对步骤3构建的BPNN模型进行参数优化；
步骤5：采用PSO算法从步骤2的吸收光谱数据中选取多个频段，将每个频段的吸收光谱数据分别输入到参数优化后的BPNN模型，预测得到每个频段对应的农药残留浓度，选择其中预测精度最高的频段，该频段记为最佳频段；
步骤6：选取与步骤1的农产品中具有相同农药品种的待测农产品，采用太赫兹光谱仪对该待测农产品进行检测，得到太赫兹吸收光谱数据，对该太赫兹吸收光谱数据进行光谱预处理，得到预处理后的吸收光谱数据，从预处理后吸收光谱数据中选取与步骤5中最佳频段相同的频段区间，将该频段区间的吸收光谱数据输入到参数优化后的BPNN模型中进行检测，得到农产品的农药残留浓度检测结果。
2.根据权利要求1所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，所述的光谱预处理具体包括：
依次采用光谱去噪法和基线校正法对太赫兹吸收光谱数据进行光谱预处理，所述光谱去噪法为小波软阈值去噪法，所述基线校正法为二阶导基线校正法。
3.根据权利要求2所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，
所述步骤3中BPNN模型的构建如下：
P＝WI；
即:
其中，W表示神经网络权重矩阵，I表示光谱输入矩阵，P表示预测浓度矩阵；wmv表示第m层网络第v个神经元的权重系数，inputn表示第n个样本的太赫兹吸收光谱数据。
4.根据权利要求2所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，所述步骤4具体为：
通过误差反向传递法进行参数优化，其中误差表达式为：
ek＝yk-Ok；
其中，yk表示样本的检测浓度，Ok表示样本的预测浓度，ek表示样本的预测误差。
5.根据权利要求2所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，所述步骤1中的农产品含有的农药品种有1种或多种，所述步骤1中某一农产品包含的农药品种与步骤6中待测农产品包含的农药品种相同。
6.根据权利要求2所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，所述农产品中含有的已知品种的农药包含6-苄氨基嘌呤、吡虫啉和6-二氯苯腈。
7.根据权利要求2所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，所述的步骤5具体为：
在0-0.50THz峰宽范围内，以0.02THz为峰宽递增间隔，并选用PSO算法计算不同峰宽下频段起始位置，继而确定不同峰宽对应的多个频段，将每个频段的吸收光谱数据分别输入到参数优化后的BPNN模型，预测得到每个频段的吸收光谱数据对应的农药残留浓度，选择其中预测精度最高的频段，该频段记为最佳频段。
8.根据权利要求7所述的基于太赫兹光谱的农产品中农药残留定量检测方法，其特征在于，
若样本中农产品含有的农药品种有1种，计算每个频段的吸收光谱数据对应的农药残留浓度预测值的相关系数或者均方根误差，选择其中相关系数最大或者均方根误差最小的频段，该频段为最佳频段；
若样本中农产品含有的农药品种有多种，计算单个频段的吸收光谱数据对应的所有农药残留浓度预测值的相关系数的平均值或者均方根误差的平均值，选择其中相关系数的平均值最大或者均方根误差的平均值最小的频段，该频段为最佳频段。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463994.7A CN114112979A (zh) | 2021-12-02 | 2021-12-02 | 一种基于太赫兹光谱的农产品中农药残留定量检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111463994.7A CN114112979A (zh) | 2021-12-02 | 2021-12-02 | 一种基于太赫兹光谱的农产品中农药残留定量检测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114112979A true CN114112979A (zh) | 2022-03-01 |
Family
ID=80365761
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111463994.7A Pending CN114112979A (zh) | 2021-12-02 | 2021-12-02 | 一种基于太赫兹光谱的农产品中农药残留定量检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114112979A (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969214A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 |
| CN103969212A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 |
| US20200217790A1 (en) * | 2019-01-09 | 2020-07-09 | University Of Shanghai For Science And Technology | Method for qualitative and quantitative determination of key substances in mixture based on terahertz spectrum |
| CN111751316A (zh) * | 2020-05-30 | 2020-10-09 | 浙江大学山东工业技术研究院 | 基于太赫兹光谱技术的农产品中2,4-d农残检测的方法 |
| CN111812058A (zh) * | 2020-05-30 | 2020-10-23 | 浙江大学山东工业技术研究院 | 基于太赫兹成像技术的香椿中农残的定性检测方法 |
| AU2020102406A4 (en) * | 2020-09-29 | 2020-12-24 | Suzhou Polytechnic Institute Of Agriculture | Method for detecting pesticide residues in agricultural products |
-
2021
- 2021-12-02 CN CN202111463994.7A patent/CN114112979A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103969214A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 |
| CN103969212A (zh) * | 2014-04-17 | 2014-08-06 | 首都师范大学 | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 |
| US20200217790A1 (en) * | 2019-01-09 | 2020-07-09 | University Of Shanghai For Science And Technology | Method for qualitative and quantitative determination of key substances in mixture based on terahertz spectrum |
| CN111751316A (zh) * | 2020-05-30 | 2020-10-09 | 浙江大学山东工业技术研究院 | 基于太赫兹光谱技术的农产品中2,4-d农残检测的方法 |
| CN111812058A (zh) * | 2020-05-30 | 2020-10-23 | 浙江大学山东工业技术研究院 | 基于太赫兹成像技术的香椿中农残的定性检测方法 |
| AU2020102406A4 (en) * | 2020-09-29 | 2020-12-24 | Suzhou Polytechnic Institute Of Agriculture | Method for detecting pesticide residues in agricultural products |
Non-Patent Citations (1)
| Title |
|---|
| 瞿芳芳: "太赫兹光谱在农产品农药残留检测中的关键技术研究", 《万方学位导航》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106124449B (zh) | 一种基于深度学习技术的土壤近红外光谱分析预测方法 | |
| CN110174359B (zh) | 一种基于高斯过程回归的航空高光谱影像土壤重金属浓度评估方法 | |
| US20100297291A1 (en) | Visible/near-infrared spectrum analyzing method and grape fermenting method | |
| CN104020135B (zh) | 基于近红外光谱的校正模型建模方法 | |
| CN101915744A (zh) | 物质成分含量的近红外光谱无损检测方法及装置 | |
| CN109324013B (zh) | 利用高斯过程回归模型构建原油性质近红外快速分析方法 | |
| CN110006844A (zh) | 基于函数性主元分析的近红外光谱特征提取方法和系统 | |
| CN110579186B (zh) | 基于逆高斯过程反演叶面积指数的农作物长势监测方法 | |
| CN108169165B (zh) | 基于太赫兹光谱和图像信息融合的麦芽糖混合物定量分析方法 | |
| CN115436407A (zh) | 一种随机森林回归结合主成分分析的元素含量定量分析方法 | |
| CN108956584B (zh) | 一种桑椹中重金属元素铬的快速准确检测方法 | |
| CN116559110A (zh) | 一种基于相关性和高斯曲线拟合的自适应近红外光谱变换方法 | |
| WO2020248961A1 (zh) | 一种无参考值的光谱波数选择方法 | |
| CN106126879B (zh) | 一种基于稀疏表示技术的土壤近红外光谱分析预测方法 | |
| CN111999258A (zh) | 一种面向光谱基线校正的加权建模局部优化方法 | |
| CN103955711B (zh) | 一种成像光谱目标识别分析中的模式识别方法 | |
| CN112651173B (zh) | 一种基于跨域光谱信息的农产品品质无损检测方法及可泛化系统 | |
| Du et al. | Quantitative detection of azodicarbonamide in wheat flour by near-infrared spectroscopy based on two-step feature selection | |
| CN111766210B (zh) | 一种近岸复杂海水硝酸盐氮多光谱测量方法 | |
| CN114112979A (zh) | 一种基于太赫兹光谱的农产品中农药残留定量检测方法 | |
| CN114624402B (zh) | 一种基于近红外光谱的螺蛳粉酸笋品质评价方法 | |
| CN114034684B (zh) | 一种基于一站式光谱模型的土壤属性快速测定方法 | |
| CN115392305A (zh) | 基于改进时间卷积网络的土壤有机质含量高光谱建模方法 | |
| CN113406037B (zh) | 一种基于序列前向选择的红外光谱在线快速识别分析方法 | |
| CN115728290A (zh) | 一种土壤中铬元素的检测方法、系统、设备及存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for 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: 20220301 |