CN112945902A - 一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺 - Google Patents
一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺 Download PDFInfo
- Publication number
- CN112945902A CN112945902A CN202110308931.8A CN202110308931A CN112945902A CN 112945902 A CN112945902 A CN 112945902A CN 202110308931 A CN202110308931 A CN 202110308931A CN 112945902 A CN112945902 A CN 112945902A
- Authority
- CN
- China
- Prior art keywords
- azodicarbonamide
- flour
- spectrum
- pure
- pixels
- 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.)
- Withdrawn
Links
- 239000004156 Azodicarbonamide Substances 0.000 title claims abstract description 138
- XOZUGNYVDXMRKW-AATRIKPKSA-N azodicarbonamide Chemical compound NC(=O)\N=N\C(N)=O XOZUGNYVDXMRKW-AATRIKPKSA-N 0.000 title claims abstract description 138
- 235000019399 azodicarbonamide Nutrition 0.000 title claims abstract description 138
- 235000013312 flour Nutrition 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005516 engineering process Methods 0.000 title claims abstract description 28
- 238000000701 chemical imaging Methods 0.000 title claims abstract description 25
- 238000001228 spectrum Methods 0.000 claims abstract description 122
- 230000003595 spectral effect Effects 0.000 claims abstract description 52
- 238000004458 analytical method Methods 0.000 claims abstract description 36
- 238000010521 absorption reaction Methods 0.000 claims abstract description 28
- 239000000203 mixture Substances 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 238000002310 reflectometry Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 238000002835 absorbance Methods 0.000 claims 1
- 230000002708 enhancing effect Effects 0.000 abstract description 3
- 239000002245 particle Substances 0.000 description 16
- 239000000843 powder Substances 0.000 description 6
- 238000000985 reflectance spectrum Methods 0.000 description 6
- 230000001965 increasing effect Effects 0.000 description 5
- 235000013305 food Nutrition 0.000 description 4
- 150000001408 amides Chemical class 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 235000017060 Arachis glabrata Nutrition 0.000 description 2
- 244000105624 Arachis hypogaea Species 0.000 description 2
- 235000010777 Arachis hypogaea Nutrition 0.000 description 2
- 235000018262 Arachis monticola Nutrition 0.000 description 2
- 240000004160 Capsicum annuum Species 0.000 description 2
- 229920000877 Melamine resin Polymers 0.000 description 2
- 241000209140 Triticum Species 0.000 description 2
- 235000021307 Triticum Nutrition 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 229910052736 halogen Inorganic materials 0.000 description 2
- 150000002367 halogens Chemical class 0.000 description 2
- 238000004128 high performance liquid chromatography Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- JDSHMPZPIAZGSV-UHFFFAOYSA-N melamine Chemical compound NC1=NC(N)=NC(N)=N1 JDSHMPZPIAZGSV-UHFFFAOYSA-N 0.000 description 2
- 235000013336 milk Nutrition 0.000 description 2
- 239000008267 milk Substances 0.000 description 2
- 210000004080 milk Anatomy 0.000 description 2
- 235000020232 peanut Nutrition 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- DUIOPKIIICUYRZ-UHFFFAOYSA-N semicarbazide Chemical compound NNC(N)=O DUIOPKIIICUYRZ-UHFFFAOYSA-N 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000001237 Raman spectrum Methods 0.000 description 1
- 208000031320 Teratogenesis Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000007844 bleaching agent Substances 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000010411 cooking Methods 0.000 description 1
- XUJNEKJLAYXESH-UHFFFAOYSA-N cysteine Natural products SCC(N)C(O)=O XUJNEKJLAYXESH-UHFFFAOYSA-N 0.000 description 1
- 235000018417 cysteine Nutrition 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002019 doping agent Substances 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000013373 food additive Nutrition 0.000 description 1
- 239000002778 food additive Substances 0.000 description 1
- 239000003205 fragrance Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000000589 high-performance liquid chromatography-mass spectrometry Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000012880 independent component analysis Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002703 mutagenesis Methods 0.000 description 1
- 231100000350 mutagenesis Toxicity 0.000 description 1
- IAIWVQXQOWNYOU-FPYGCLRLSA-N nitrofural Chemical class NC(=O)N\N=C\C1=CC=C([N+]([O-])=O)O1 IAIWVQXQOWNYOU-FPYGCLRLSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012628 principal component regression Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000000273 veterinary drug Substances 0.000 description 1
- 230000000007 visual 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/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/38—Diluting, dispersing or mixing samples
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (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
本发明公开了一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，包括如下步骤：S1：取纯偶氮甲酰胺样品用高光谱成像仪获得纯偶氮甲酰胺的平均漫反射光谱,并从纯偶氮甲酰胺的平均漫反射光谱中提取出明显的吸收峰；本发明探讨应用近红外高光谱技术检测面粉中偶氮甲酰胺。首先，采集纯偶氮甲酰胺、纯面粉和含不同浓度偶氮甲酰胺面粉样品的高光谱图像，找到偶氮甲酰胺区别于面粉的特征吸收波段；其次，通过二阶导数预处理增强光谱中吸收峰的分辨，依据纯偶氮甲酰胺和纯面粉样品确定SAM、SCA、SCM的阈值；最后，分析三种光谱相似性分析方法对10种不同浓度混合样品的分类结果，并对其分类正确性进行验证。
Description
技术领域
本发明涉及监测技术领域，更具体地涉及应用近红外高光谱成像技术检测面粉中偶氮甲酰胺。
背景技术
高光谱成像技术不仅可以获得样品的图像信息，每个像素点还包含了光谱信息，因其信息量丰富的特点已在食品安全检测方面得到了应用。本发明应用近红外高光谱成像技术检测面粉中偶氮甲酰胺。分别采集纯偶氮甲酰胺、纯面粉和面粉中10种不同浓度偶氮甲酰胺混合样品的高光谱图像。通过比较纯偶氮甲酰胺和纯面粉的平均漫发射光谱，找到两者区分度较大的4处吸收波段：1574.38nm、2038.55nm、2166.88nm和2269.91nm。采用二阶导数对样品图像中的像素点光谱进行预处理，通过光谱角制图、光谱相关角和光谱相关性度量三种光谱相似性分析方法对混合样品中的偶氮甲酰胺像素和面粉像素进行检测。结果表明，预处理后的平均光谱不能有效检测面粉中偶氮甲酰胺；单像素点光谱结合光谱相似性分析实现了混合样品中偶氮甲酰胺像素和面粉像素的分类；分类结果的验证显示了偶氮甲酰胺像素和面粉像素的正确分类。发明结果为利用高光谱技术检测面粉中添加剂提供了方法支持，为食品中掺杂物的检测提供参考。
面粉作为世界上主要食品原料之一，其制品以色、香、味俱全著称。面粉的质量直接关乎人们的健康，食品添加剂的使用备受关注。偶氮甲酰胺作为一种漂白剂和面粉改良剂用于面粉的熟化处理，它能够氧化小麦粉中的半胱氨酸，从而使面粉筋度增加，提高面团气体保留量，增加烘焙制品的弹性和韧性，改善面团的可操作性和调理性。但偶氮甲酰胺极其不稳定，在加热过程中它会分解从而形成微量的氨基脲，氨基脲是呋喃西林兽药的代谢物，具有很强的致癌、致畸、致突变等副作用。
目前，对面粉中偶氮甲酰胺的传统检测方法主要有高效液相色谱法和高效液相色谱-质谱联用法，这两种方法具有较高的检测精度，但前处理复杂、检测时间长、需要熟练的操作人员，一般只限于实验室应用。随着光谱技术的发展，一些光谱方法在面粉添加剂检测中展示了良好的应用，如近红外光谱、太赫兹光谱和拉曼光谱等。这些方法均是对试样的单点检测，覆盖的空间范围不能满足对整个样品的检测要求。
高光谱成像技术是是光谱分析和机器视觉技术的高精度融合技术，兼有两种技术的优势，采集到的高光谱图像中既包含了待测物内部品质的光谱信息，又包含了待测物外部品质的图像信息。根据光源和光谱相机位置的不同，高光谱图像的获取方式可分为反射、透射和漫反射三种，其中漫反射成像在检测食品污染物方面展示了良好的应用。Mishra等应用独立成分分析中的特征矩阵联合近似对角化算法处理近红外高光谱图像以检测面粉中痕量花生粉，实现了花生粉痕量分布的可视化识别。Lim等使用近红外高光谱成像技术结合偏最小二乘模型的回归系数检测奶粉中的三聚氰胺颗粒，实现了奶粉中三聚氰胺的有效检测，最低检测浓度为0.02％。周瑶等利用高光谱检测辣椒粉中的苏丹红一号，通过比较偏最小二乘回归法、多元线性回归法和主成分回归法建立预测模型的结果，选择多元线性回归法实现辣椒粉中的苏丹红一号的无损检测。
发明内容
本发明为了克服现有技术的不足，提供一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺
为了实现上述目的，本发明采用以下技术方案：
一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，包括如下步骤：
S1：取纯偶氮甲酰胺样品用高光谱成像仪获得纯偶氮甲酰胺的平均漫反射光谱,并从纯偶氮甲酰胺的平均漫反射光谱中提取出明显的吸收峰；
S2：配置偶氮甲酰胺和面粉不同比例的混合物，并获取纯面粉和上述混合物的平均漫反射光谱，和各个像素点的漫反射光谱；
S3：对混合物各个像素点光谱进行二阶导数预处理，以获得混合物的各个像素点的二阶导数光谱，二阶导数预处理纯面粉混合物的平均漫反射光谱，以获得纯面粉的二阶导数光谱，将混合物各个像素点的二阶导数光谱分别减去纯面粉的二阶导数光谱后得到的光谱作用于SAM、SCA以及SCM三种光谱相似性分析算法中，以分析出混合物的各个像素点中的面粉像素和偶氮甲酰胺像素；
S4：将所述S3中的面粉像素和偶氮甲胺像素的平均二阶导数光谱和所述S1中纯偶氮甲酰胺的吸收峰进行对比，以验证偶氮甲酰胺像素和面粉像素的分类结果。
本发明探讨应用近红外高光谱技术检测面粉中偶氮甲酰胺。首先，采集纯偶氮甲酰胺、纯面粉和含不同浓度偶氮甲酰胺面粉样品的高光谱图像，找到偶氮甲酰胺区别于面粉的特征吸收波段；其次，通过二阶导数预处理增强光谱中吸收峰的分辨，依据纯偶氮甲酰胺和纯面粉样品确定SAM、SCA、SCM的阈值；最后，分析三种光谱相似性分析方法对10种不同浓度混合样品的分类结果，并对其分类正确性进行验证。
可选的，所述S3中，对混合样品进行相似性分析的具体方式如下，首先计算纯偶氮甲酰胺样品中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，然后计算纯面粉中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，阈值即为纯偶氮甲酰胺结果的最大值和纯面粉结果最小值的中值。
可选的，SAM、SCA和SCM三种光谱相似性分析算法的阈值分别为1.4109、0.9525和0.1590，因此，具有<1.4109的角度值，<0.9525的角度值和>0.1590的相关值被分类为偶氮甲酰胺像素。
可选的，所述S4中，由SCM光谱相似性分析算法获得的偶氮甲酰胺像素的光谱于226.91nm处显示出强吸收峰，且SCM光谱相似性分析算法获得的面粉像素的光谱未在226.91nm处显示出吸收峰。
可选的，所述S1中，需要采集白色和黑色标定参考图像，以校正原始图像中CCD检测器暗电流和光强度分布不均匀的影响，校正公式如下：
其中，Rraw为高光谱采集的原始的图像；Rdark为关闭光源，盖子覆盖镜头时采集的黑色参考图像；Rwhite为采集的99％反射率标准白板的白色参考图像；Rcal为校正后获得的相对光谱图像。
可选的，所述S1中，偶氮甲酰胺的吸收峰为1574.38nm、2038.55nm、2166.88nm和2269.91nm。
可选的，所述S1和所述S2中，纯偶氮甲酰胺、纯面粉以及上述混合物的平均漫反射的光谱获取方式如下：将需要采集平均漫反射光谱的样品放置到所述高光谱成像仪相机镜头的正下方，调整焦距使镜头到样品表面的距离为40cm；光谱采集范围为：1398.11～2502.89nm，光谱分辨率为8nm，曝光时间为2ms，所述高光谱成像仪的定位平台的移动速度为43mm·s-1。
综上所述，本发明探讨应用近红外高光谱技术检测面粉中偶氮甲酰胺。首先，采集纯偶氮甲酰胺、纯面粉和含不同浓度偶氮甲酰胺面粉样品的高光谱图像，找到偶氮甲酰胺区别于面粉的特征吸收波段；其次，通过二阶导数预处理增强光谱中吸收峰的分辨，依据纯偶氮甲酰胺和纯面粉样品确定SAM、SCA、SCM的阈值；最后，分析三种光谱相似性分析方法对10种不同浓度混合样品的分类结果，并对其分类正确性进行验证。
附图说明
图1近红外高光谱成像系统示意图
图2纯偶氮甲酰胺、纯面粉和偶氮甲酰胺-面粉混合样品的平均漫反射光谱
图4纯偶氮甲酰胺、纯面粉和偶氮甲酰胺-面粉混合样品在特征吸收波段处的灰度图像
图3面粉中不同浓度偶氮甲酰胺混合样品的平均二阶导数光谱
图5 1％浓度混合样品感兴趣区域内各像素点的原始光谱(A)和二阶导数光谱(B)
图6经SAM，SCA和SCM方法分类为偶氮甲酰胺像素和面粉像素的图像
图7组2中不同浓度混合样品经SAM分类为偶氮甲酰胺像素和面粉像素的平均二阶导数光谱
具体实施方式
为了使本技术领域的人员更好的理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述。
1实验部分
1.1仪器与试剂
高光谱成像系统如图1所示，该系统由焦平面阵列相机(Xeva-FPA-2.5-320(100Hz)，Xenics Ltd.，Leuven，Belgium)、成像光谱仪3(ImSpector N25E，SpectralImaging Ltd.，Oulu,Finland)、变焦相机镜头4(OLES30f-2.0/30mm，Spectral ImagingLtd.，Oulu，Finland)、两个150W卤素灯5(AT-500WB，Antefore International CO.，LTD.Taiwan)、步进电机8、单轴运动平台7和计算机1组成。其工作原理为：卤素灯4发出的光线投射到样品表面，光线和样品6相互作用后发生漫反射，该信号由相机2和光谱仪3采集，数据经USB电缆传输至电脑。样品放置于单轴运动平台7上，移动范围为0～40cm，通过与计算机连接的步进电机进行控制。整个系统除计算机1外，放置在一个黑色的箱子中，以避免采集过程中环境光的影响。系统的控制和数据的采集通过Spectral Image软件(五铃光学股份有限公司开发)完成。
偶氮甲酰胺，分析纯，购于山东西亚化学工业有限公司；家用小麦粉，购于当地超市，经高效液相色谱法检测不含有偶氮甲酰胺；培养皿(外径39mm，内径37mm，深5mm)，购于科晶生物科技有限公司；涡旋混合器，Vortex-Genie2型，美国Scientific Industries公司。
1.2实验方法
1.2.1样品制备
称取偶氮甲酰胺样品10g，然后将其转入培养皿中，使样品填满整个培养皿并使得样品表面与培养皿的上边缘平齐。培养皿的样品容量约为3g，为充分利用制备的10g样品，每个样品分为3组进行采样。以同样的方式制备面粉样品。
将一定量的偶氮甲酰胺和面粉样品装入50ml离心管中，通过涡旋混合器混合均匀，制得10种不同浓度(0.02％、0.04％、0.06％、0.08％、0.1％、0.2％、0.4％、0.6％、0.8％和1.0％，W/W)的偶氮甲酰胺-面粉混合样品。每种混合样品的质量为10g，例如，最高混合浓度(1％)通过0.1g偶氮甲酰胺和9.9g面粉混合制得；最低混合浓度(0.01％)通过0.001g偶氮甲酰胺和9.999g面粉混合制得。将混合样品转入培养皿中，使混合样品填满整个培养皿并使得样品表面与培养皿的上边缘平齐。培养皿的样品容量约为3g，为充分利用制备的10g样品，每个混合样品分为3组进行采样。
1.2.2数据采集
将制备好的样品放置于单轴机动定位平台上相机镜头的正下方位置，调整焦距使镜头到样品表面的距离为40cm。光谱采集范围为：1398.11～2502.89nm，光谱分辨率为8nm，曝光时间为2ms，定位平台的移动速度为43mm·s-1。
1.2.3数据处理
(1)图像校正
在与采集样本图像相同的实验条件下，采集白色和黑色标定参考图像，以校正原始图像中CCD检测器暗电流和光强度分布不均的影响。校正方法如式(1)
其中，Rraw为高光谱采集的原始的图像；Rdark为关闭光源，盖子覆盖镜头时采集的黑色参考图像；Rwhite为采集的99％反射率标准白板的白色参考图像；Rcal为校正后获得的相对光谱图像。
(2)光谱数据处理
在每个直径为37mm的圆形样品区域内选择直径约为33mm(图像中心周围)的感兴趣区域(Region of Interest，ROI)，以确保只包含面粉区域，排除背景区域的影响。将三维高光谱数据转换为单个像素点的二维漫反射光谱，对每条光谱进行二阶导数(SecondDerivatives，2nd)预处理以消除基线漂移的影响。从偶氮甲酰胺-面粉混合样品的单个像素光谱中减去纯面粉的平均光谱以去除混合样品中的面粉背景，通过光谱相似性分析将得到的差谱与偶氮甲酰胺的平均光谱进行比较。所有数据分析基于MATAB 7.11(Mathworks,USA)软件完成。
1.2.4光谱相似性分析
本发明采用三种光谱相似性分析算法，包括光谱角制图(Spectral AngleMapper，SAM)、光谱相关角(Spectral Correlation Angle，SCA)和光谱相关性度量(Spectral Correlation Measure，SCM)。
(1)光谱角制图(SAM)
SAM是由Kruse,F.A.等提出的一种光谱匹配技术，其原理是根据目标光谱和参考光谱之间的光谱角度来判断两者的相似性。SAM的计算方法如公式(1)：
其中，t为目标光谱的矢量，r为参考光谱的矢量，n为波段数。α为两条光谱之间的夹角，其取值范围位于[0，π/2]之间，其值越小，表示两条光谱的相似度越高。
(2)光谱相关角(SCA)
SCA是目标光谱和参考光谱间Pearson相关系数的余弦角，其计算公式为：
SCA值越小，表示目标光谱和参考光谱的相似度越高。SC表示两光谱的相关系数，其值介于-1和1之间，计算方法如下：
(3)光谱相关性度量(SCM)
Meer等将n定义为重叠位置，引入了用于测量目标光谱t和参考光谱r的相似性方法，计算方法为：
SCA值越大，表示目标光谱和参考光谱的相似度越高。
根据三种光谱相似性算法的结果，为每种算法设置阈值，通过阈值分割实现ROI内各像素点光谱的分类，从而创建面粉中偶氮甲酰胺颗粒的分布图像。
2结果与讨论
2.1样品的光谱和图像
图2所示为纯偶氮甲酰胺、纯面粉和混合样品的平均漫反射光谱。由图可看出，偶氮甲酰胺的漫反射光谱同面粉具有显著差异，偶氮甲酰胺在1574.38nm、2038.55nm、2166.88nm和2269.91nm处具有明显的吸收峰。偶氮甲酰胺分子式为C2H4N4O2，NH是主要的含氢基团，1574.38nm为NH对称振动和NH反对称振动的合频，2038.55nm为NH对称振动和酰胺II变形振动的合频，2166.88nm为酰胺I和酰胺III变形振动的合频，2269.91nm为NH对称振动和NH2摇摆振动的合频。这4处吸收峰反映了偶氮甲酰胺分子结构内化学键的振动，也反映出与面粉组分的区别。
在图2中，面粉中不同浓度偶氮甲酰胺混合样品的平均漫反射光谱之间无显著差异，未表现出偶氮甲酰胺的光谱特征，且类似于纯面粉的漫反射光谱。对混合样品的像素点光谱进行二阶导数预处理，以消除基线漂移的影响，更好的分辨光谱中吸收峰。图3所示为面粉中不同浓度偶氮甲酰胺混合样品的平均二阶导数光谱。由图可看出，在特征波段2269.91nm处，当偶氮甲酰胺浓度从1％下降到0.2％时，反射率不断降低；当偶氮甲酰胺浓度小于0.2％时，无明显反射率变化，平均二阶导数光谱相似。这表明，基于整个区域的平均光谱不能有效检测面粉中偶氮甲酰胺含量，因此，本发明采用对ROI内单个像素点的光谱进行评估，以更好的检测面粉中偶氮甲酰胺颗粒。
图4所示为纯偶氮甲酰胺、纯面粉和混合样品在偶氮甲酰胺特征吸收波段(1574.38nm、2038.55nm、2166.88nm和2269.91nm)的灰度图像。由图可看出，不同波段对应的图像的亮度不同，这是因为同一样品在不同波段下样品的漫反射率不同，以偶氮甲酰胺样品为例，其图像由明亮到灰暗的变化顺序为：1574.38nm、2166.88nm、2038.55nm、2269.91nm，这和图2中偶氮甲酰胺漫反射率由强到弱的变化顺序相一致。在同一波段下，纯偶氮甲酰胺、纯面粉和混合样品的图像在亮度存在差异，纯偶氮甲酰胺图像均比纯面粉和混合样品的图像更暗，然而在混合样品的图像中仍无法区分面粉和偶氮甲酰胺像素。因此，采用光谱相似性分析方法增加面粉和偶氮甲酰胺的对比度，以识别出偶氮甲酰胺像素。
2.2阈值的选取
图5所示为1％浓度混合样品ROI内各像素点的原始光谱(A)和二阶导数光谱(B)。对比图5A和5B可看出，经二阶导数预处理后，消除了基线漂移的影响，增强了光谱的分辨，一些偶氮甲酰胺的吸收峰可观察到。对不同浓度混合样品ROI像素点的光谱均进行二阶导数预处理，从预处理光谱中减去纯面粉预处理后的平均光谱，以减弱面粉背景信号的影响。将相减后的光谱用于SAM、SCA和SCM三种光谱相似性分析算法中，以识别混合样品中的偶氮甲酰胺像素。
在对不同浓度的混合样品进行光谱相似性分析前需确定阈值，以实现面粉像素和偶氮甲酰胺像素的分类。阈值采用中值法确定，首先计算纯偶氮甲酰胺样品中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，然后计算纯面粉中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，阈值即为纯偶氮甲酰胺结果的最大值和纯面粉结果最小值的中值。通过计算，得到SAM、SCA和SCM三种光谱相似性分析算法的阈值分别为1.4109、0.9525和0.1590，因此，具有<1.4109的角度值，<0.9525的角度值和>0.1590的相关值被分类为偶氮甲酰胺像素。基于三种光谱相似性分析的阈值，纯偶氮甲酰胺和纯面粉ROI中的所有像素被正确地分类为偶氮甲酰胺像素和面粉像素。
2.3光谱相似性分析结果与验证
图6所示为三组混合样品经SAM，SCA和SCM方法分类后的图像。其中，深色像素表示偶氮甲酰胺像素，浅色表示面粉像素；各样品的排列顺序同图4中样品排列顺序一致。在各组样品中，三种光谱相似性算法均较好的实现了不同浓度混合样品中偶氮甲酰胺颗粒的分类。对于同一组样品，SAM，SCA和SCM方法的分类结果基本一致。在0.2％～1％浓度的混合样品中，图像中显示出较多的偶氮甲酰胺像素分布且递增趋势明显；在0.02％～0.2％浓度的混合样品中，图像中仅有较少的偶氮甲酰胺像素分布，对各混合样品中分类为偶氮甲酰胺的像素进行统计，其占ROI区域内总像素点(5016个)的比例如表1所示。从表中可看出，同一浓度的三组样品中，被分类为偶氮甲酰胺颗粒的数量存在较大差别，这是因为近红外高光谱系统不仅从样品的表面采集漫发射信号，而且还可以采集样品内部的漫反射信号，但在样品内部，不同层所含有的偶氮甲酰胺颗粒的数量不同，与其他层中偶氮甲酰胺颗粒重叠的数量也不同。不同浓度的样品分类为偶氮甲酰胺颗粒的比例大于面粉中偶氮甲酰胺的实际浓度，在浓度大于0.1％的混合样品中更为明显，这是因为高浓度的混合样品中含有较多的偶氮甲酰胺颗粒，不同层内均可采集到偶氮甲酰胺颗粒的漫反射信号，图像中偶氮甲酰胺颗粒的分布来自于多层的映射。非最低浓度样品中存在没有分类为偶氮甲酰胺颗粒的现象(如组1中SAM方法分类的0.04％浓度混合样品)，出现这种现象的原因有两种，一种是混合样品虽然经过混匀，但浓度较低的样品转入培养皿后表面不含有偶氮甲酰胺颗粒，由于面粉层过厚未采集到位于底层的偶氮甲酰胺颗粒的漫发射信号；一种是出现了假阴性的可能。因此，为避免高浓度样品过高估计样品及样品过厚导致底部信号无法采集，下一步可对样品的穿透深度进行发明，以采用较薄的样品厚度进行高光谱成像。
表1 SAM，SCA和SCM分类为偶氮甲酰胺像素占总像素的百分比
Table 1 Percent of total pixels classified as azodicarbonamide pixelsby SAM,SCA,and SCM.
为进一步评估光谱相似性分析的分类结果，对分类为偶氮甲酰胺像素和面粉像素的光谱进行验证。在三组样品中，三种光谱相似性分析算法的分类结果基本一致，组1和组3的低浓度混合样品中存在偶氮甲酰胺颗粒分类为零的情况，因此选择SAM算法对组2中不同浓度混合的分类结果进行光谱分析。图7所示为组2中不同浓度混合样品经SAM分类为偶氮甲酰胺像素和面粉像素的平均二阶导数光谱。图中仅显示了偶氮甲酰胺最强吸收峰226.91nm附近的光谱，其中分类为偶氮甲酰胺像素的平均光谱以实线表示，分类为面粉像素的平均光谱以虚线表示。由图可看出，分类为偶氮甲酰胺像素的平均光谱显示出在2269.91nm的吸收峰，而面粉像素的平均光谱中未出现该吸收峰，说明利用光谱相似性分析方法可以正确分类偶氮甲酰胺像素和面粉像素。
3结论
本文通过近红外高光谱技术实现了面粉中偶氮甲酰胺的检测。主要结论有：
(1)通过分析纯偶氮甲酰胺的平均漫反射光谱，找到了区别于面粉漫反射光谱的4处特征吸收波段：1574.38nm、2038.55nm、2166.88nm和2269.91nm。
(2)面粉中含不同浓度偶氮甲酰胺的混合样品的平均光谱未能显示出偶氮甲酰胺吸收特征，特征波段对应灰度图像无法区分面粉和偶氮甲酰胺像素，拟采用光谱相似性分析对单像素点光谱进行计算，从而实现不同浓度偶氮甲酰胺-面粉混合样品的分类。
(3)三种光谱相似性分析算法(SAM，SCA和SCM)实现了面粉中不同浓度偶氮甲酰胺的分类，但对偶氮甲酰胺的分类结果中存在数量差别大、浓度过高估计和分类为零的问题，后续的发明中，采用较薄的样品厚度进行高光谱成像可实现更好的分类效果。
(4)对分类为偶氮甲酰胺像素和面粉像素所对应光谱的验证结果显示，光谱相似性分析方法可以正确分类混合样品中偶氮甲酰胺像素和面粉像素。
显然，所描述的实施例仅仅本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
Claims (8)
1.一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，包括如下步骤：
S1：取纯偶氮甲酰胺样品用高光谱成像仪获得纯偶氮甲酰胺的平均漫反射光谱,并从纯偶氮甲酰胺的平均漫反射光谱中提取出明显的吸收峰；
S2：配置偶氮甲酰胺和面粉不同比例的混合物，并获取纯面粉和上述混合物的平均漫反射光谱，和各个像素点的漫反射光谱；
S3：对混合物各个像素点光谱进行二阶导数预处理，以获得混合物的各个像素点的二阶导数光谱，二阶导数预处理纯面粉混合物的平均漫反射光谱，以获得纯面粉的二阶导数光谱，将混合物各个像素点的二阶导数光谱分别减去纯面粉的二阶导数光谱后得到的光谱作用于SAM、SCA以及SCM三种光谱相似性分析算法中，以分析出混合物的各个像素点中的面粉像素和偶氮甲酰胺像素；
S4：将所述S3中的面粉像素和偶氮甲胺像素的平均二阶导数光谱和所述S1中纯偶氮甲酰胺的吸收峰进行对比，以验证偶氮甲酰胺像素和面粉像素的分类结果。
2.根据权利要求1所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S3中，所述SAM的计算公式如下，
其中，t为目标光谱的矢量，r为参考光谱的矢量，n为波段数。α为两条光谱之间的夹角，其取值范围位于[0，π/2]之间；
光谱相关角(SCA)的计算公式如下，
SCA是目标光谱和参考光谱间Pearson相关系数的余弦角，其计算公式为：
SCA值越小，表示目标光谱和参考光谱的相似度越高。SC表示两光谱的相关系数，其值介于-1和1之间，计算方法如下：
光谱相关性度量(SCM)的计算公式如下，
Meer等将n定义为重叠位置，引入了用于测量目标光谱t和参考光谱r的相似性方法，计算方法为：
3.根据权利要求2所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S3中，对混合样品进行相似性分析的具体方式如下，首先计算纯偶氮甲酰胺样品中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，然后计算纯面粉中各像素点和偶氮甲酰胺平均光谱的光谱相似性分析的最大值和最小值，阈值即为纯偶氮甲酰胺结果的最大值和纯面粉结果最小值的中值。
4.根据权利要求3所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，SAM、SCA和SCM三种光谱相似性分析算法的阈值分别为1.4109、0.9525和0.1590，因此，具有<1.4109的角度值，<0.9525的角度值和>0.1590的相关值被分类为偶氮甲酰胺像素。
5.根据权利要求4所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S4中，由SCM光谱相似性分析算法获得的偶氮甲酰胺像素的光谱于226.91nm处显示出强吸收峰，且SCM光谱相似性分析算法获得的面粉像素的光谱未在226.91nm处显示出吸收峰。
6.根据权利要求1所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S1中，需要采集白色和黑色标定参考图像，以校正原始图像中CCD检测器暗电流和光强度分布不均匀的影响，校正公式如下：
其中，Rraw为高光谱采集的原始的图像；Rdark为关闭光源，盖子覆盖镜头时采集的黑色参考图像；Rwhite为采集的99％反射率标准白板的白色参考图像；Rcal为校正后获得的相对光谱图像。
7.根据权利要求1所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S1中，偶氮甲酰胺的吸收峰为1574.38nm、2038.55nm、2166.88nm和2269.91nm。
8.根据权利要求1所述的一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺，其特征在于，所述S1和所述S2中，纯偶氮甲酰胺、纯面粉以及上述混合物的平均漫反射的光谱获取方式如下：将需要采集平均漫反射光谱的样品放置到所述高光谱成像仪相机镜头的正下方，调整焦距使镜头到样品表面的距离为40cm；光谱采集范围为：1398.11～2502.89nm，光谱分辨率为8nm，曝光时间为2ms，所述高光谱成像仪的定位平台的移动速度为43mm·s-1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110308931.8A CN112945902A (zh) | 2021-03-23 | 2021-03-23 | 一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110308931.8A CN112945902A (zh) | 2021-03-23 | 2021-03-23 | 一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112945902A true CN112945902A (zh) | 2021-06-11 |
Family
ID=76228472
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110308931.8A Withdrawn CN112945902A (zh) | 2021-03-23 | 2021-03-23 | 一种应用近红外高光谱成像技术检测面粉中偶氮甲酰胺 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112945902A (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104749132A (zh) * | 2015-04-08 | 2015-07-01 | 中国计量学院 | 一种测定面粉中偶氮甲酰胺的含量的方法 |
| CN106353274A (zh) * | 2016-11-17 | 2017-01-25 | 福州大学 | 一种面粉中偶氮甲酰胺的比色检测方法 |
-
2021
- 2021-03-23 CN CN202110308931.8A patent/CN112945902A/zh not_active Withdrawn
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104749132A (zh) * | 2015-04-08 | 2015-07-01 | 中国计量学院 | 一种测定面粉中偶氮甲酰胺的含量的方法 |
| CN106353274A (zh) * | 2016-11-17 | 2017-01-25 | 福州大学 | 一种面粉中偶氮甲酰胺的比色检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王晓彬 等: "应用近红外高光谱成像技术检测面粉中偶氮甲酰胺", 《光谱学与光谱分析》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zheng et al. | Detection of adulteration with duck meat in minced lamb meat by using visible near-infrared hyperspectral imaging | |
| Zhang et al. | Nondestructive measurement of soluble solids content in apple using near infrared hyperspectral imaging coupled with wavelength selection algorithm | |
| ElMasry et al. | Chemical-free assessment and mapping of major constituents in beef using hyperspectral imaging | |
| Wu et al. | Application of time series hyperspectral imaging (TS-HSI) for determining water distribution within beef and spectral kinetic analysis during dehydration | |
| ElMasry et al. | Principles of hyperspectral imaging technology | |
| Wu et al. | Advanced applications of hyperspectral imaging technology for food quality and safety analysis and assessment: A review—Part I: Fundamentals | |
| Wu et al. | Rapid prediction of moisture content of dehydrated prawns using online hyperspectral imaging system | |
| Burger et al. | Hyperspectral NIR imaging for calibration and prediction: a comparison between image and spectrometer data for studying organic and biological samples | |
| Park et al. | Performance of hyperspectral imaging system for poultry surface fecal contaminant detection | |
| CN105181606B (zh) | 基于高光谱成像技术检测花生中蔗糖含量分布的方法 | |
| CN112014331A (zh) | 一种水体污染的检测方法、装置、设备以及存储介质 | |
| Kim et al. | Line-scan hyperspectral imaging platform for agro-food safety and quality evaluation: System enhancement and characterization | |
| Wang et al. | Smartphone imaging spectrometer for egg/meat freshness monitoring | |
| CN113607676B (zh) | 基于异形比色皿和近红外图像的牛奶成分快速检测方法 | |
| CN1553165A (zh) | 一种多光谱细胞涂片自动分析仪及其用于宫颈细胞的分析方法 | |
| CN111537469A (zh) | 一种基于近红外技术的苹果品质快速无损检测方法 | |
| Yu et al. | Rapid and visual measurement of fat content in peanuts by using the hyperspectral imaging technique with chemometrics | |
| Duan et al. | Sensitive variables extraction, non-destructive detection and visualization of total viable count (TVC) and pH in vacuum packaged lamb using hyperspectral imaging | |
| Baek et al. | A novel hyperspectral line-scan imaging method for whole surfaces of round shaped agricultural products | |
| Mäkelä et al. | Hyperspectral near infrared image calibration and regression | |
| CN109991174A (zh) | 一种高光谱相机与普通光谱仪之间的模型转移方法 | |
| Schults et al. | Digital colorimetry: analytical possibilities and prospects of use | |
| Wang et al. | Near‐infrared hyperspectral imaging for detection and quantification of azodicarbonamide in flour | |
| CN111521575B (zh) | 质控物质选择方法及装置 | |
| KR102099164B1 (ko) | 초분광 nir 영상을 이용한 인삼의 내부 내백 등급 판정시스템 및 판정방법 |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210611 |
|
| WW01 | Invention patent application withdrawn after publication |