CN107024445A - 蔬菜中硝酸盐的快速检测的建模方法和检测方法 - Google Patents
蔬菜中硝酸盐的快速检测的建模方法和检测方法 Download PDFInfo
- Publication number
- CN107024445A CN107024445A CN201710249286.0A CN201710249286A CN107024445A CN 107024445 A CN107024445 A CN 107024445A CN 201710249286 A CN201710249286 A CN 201710249286A CN 107024445 A CN107024445 A CN 107024445A
- Authority
- CN
- China
- Prior art keywords
- nitrate
- spectrum
- modeling
- vegetables
- sample
- 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
- 229910002651 NO3 Inorganic materials 0.000 title claims abstract description 98
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 title claims abstract description 98
- 235000013311 vegetables Nutrition 0.000 title claims abstract description 85
- 238000001514 detection method Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000001228 spectrum Methods 0.000 claims abstract description 74
- 239000007788 liquid Substances 0.000 claims abstract description 25
- 239000000126 substance Substances 0.000 claims abstract description 19
- 238000012216 screening Methods 0.000 claims abstract description 15
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 12
- 230000003595 spectral effect Effects 0.000 claims abstract description 11
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 150000002823 nitrates Chemical class 0.000 claims abstract description 4
- 235000011389 fruit/vegetable juice Nutrition 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 238000011156 evaluation Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- BFJSDDYECHPBIO-UHFFFAOYSA-N 2-hydroxybenzoic acid;sulfuric acid Chemical compound OS(O)(=O)=O.OC(=O)C1=CC=CC=C1O BFJSDDYECHPBIO-UHFFFAOYSA-N 0.000 claims description 5
- 238000004611 spectroscopical analysis Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000000205 computational method Methods 0.000 claims description 4
- 235000021384 green leafy vegetables Nutrition 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 235000014121 butter Nutrition 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000009499 grossing Methods 0.000 claims description 2
- 239000003153 chemical reaction reagent Substances 0.000 abstract description 6
- 244000221633 Brassica rapa subsp chinensis Species 0.000 description 12
- 235000010149 Brassica rapa subsp chinensis Nutrition 0.000 description 12
- 235000013305 food Nutrition 0.000 description 9
- 238000010521 absorption reaction Methods 0.000 description 5
- IOVCWXUNBOPUCH-UHFFFAOYSA-M Nitrite anion Chemical compound [O-]N=O IOVCWXUNBOPUCH-UHFFFAOYSA-M 0.000 description 4
- 241000219315 Spinacia Species 0.000 description 4
- 235000009337 Spinacia oleracea Nutrition 0.000 description 4
- 238000002790 cross-validation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 150000002500 ions Chemical class 0.000 description 4
- HEMHJVSKTPXQMS-UHFFFAOYSA-M Sodium hydroxide Chemical compound [OH-].[Na+] HEMHJVSKTPXQMS-UHFFFAOYSA-M 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 230000008520 organization Effects 0.000 description 3
- 238000010802 RNA extraction kit Methods 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 238000004737 colorimetric analysis Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000003651 drinking water Substances 0.000 description 2
- 235000020188 drinking water Nutrition 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000002329 infrared spectrum Methods 0.000 description 2
- 238000001834 photoacoustic spectrum Methods 0.000 description 2
- 230000037039 plant physiology Effects 0.000 description 2
- FGIUAXJPYTZDNR-UHFFFAOYSA-N potassium nitrate Chemical compound [K+].[O-][N+]([O-])=O FGIUAXJPYTZDNR-UHFFFAOYSA-N 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 235000000536 Brassica rapa subsp pekinensis Nutrition 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 240000008436 Ipomoea aquatica Species 0.000 description 1
- 235000019004 Ipomoea aquatica Nutrition 0.000 description 1
- GRYLNZFGIOXLOG-UHFFFAOYSA-N Nitric acid Chemical compound O[N+]([O-])=O GRYLNZFGIOXLOG-UHFFFAOYSA-N 0.000 description 1
- IOVCWXUNBOPUCH-UHFFFAOYSA-N Nitrous acid Chemical class ON=O IOVCWXUNBOPUCH-UHFFFAOYSA-N 0.000 description 1
- 229910019142 PO4 Inorganic materials 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 230000003698 anagen phase Effects 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 230000000711 cancerogenic effect Effects 0.000 description 1
- 231100000357 carcinogen Toxicity 0.000 description 1
- 239000003183 carcinogenic agent Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 235000005911 diet Nutrition 0.000 description 1
- 230000000378 dietary effect Effects 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 239000002778 food additive Substances 0.000 description 1
- 235000013373 food additive Nutrition 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- DOTMOQHOJINYBL-UHFFFAOYSA-N molecular nitrogen;molecular oxygen Chemical compound N#N.O=O DOTMOQHOJINYBL-UHFFFAOYSA-N 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 229910017604 nitric acid Inorganic materials 0.000 description 1
- 239000000618 nitrogen fertilizer Substances 0.000 description 1
- GQPLMRYTRLFLPF-UHFFFAOYSA-N nitrous oxide Inorganic materials [O-][N+]#N GQPLMRYTRLFLPF-UHFFFAOYSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- NBIIXXVUZAFLBC-UHFFFAOYSA-K phosphate Chemical compound [O-]P([O-])([O-])=O NBIIXXVUZAFLBC-UHFFFAOYSA-K 0.000 description 1
- 239000010452 phosphate Substances 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000002798 spectrophotometry method Methods 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000000870 ultraviolet spectroscopy Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
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/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法:蔬菜分为两份,一份化学法测得实验对照数据;另一份取6.5‑8 μm波数特征吸收光谱图为建模样本;配制一系列硝酸盐溶液标准溶液并采集光谱信息作建模样本;将二曲线预处理扣除干扰;将蔬菜硝酸盐光谱曲线及标准溶液光谱曲线与化学参考值建模确立相关性;获得待测蔬菜的特征吸收光谱图并将光谱预处理扣除干扰;将硝酸盐光谱曲线与模型智能鉴定分类,选择最佳匹配建模样本及建模样本数,通过智能算法获取最具相似特性的样本,预测其硝酸盐含量进行比对获得蔬菜的食用信息。本发明安全、不使用化学试剂,适合现场快速获取检测效率高。
Description
技术领域
本发明涉及食品中硝酸盐的检测的建模方法和检测方法,具体涉及一种基于红外光谱的硝酸盐的快速定量测定方法,属农业食品检测领域。
背景技术
随着居民生活水平的提高和科技的发展,食品健康越来越受到关注。研究表明,硝酸盐在人体经转化可被还原为亚硝酸盐,亚硝酸盐是强致癌物,持续摄入0.3-0.5g亚硝酸盐可引起中毒,摄入3.0g亚硝酸盐可致人死亡;硝酸盐进入人体内被肠胃细菌还原为亚硝酸盐,亚硝酸盐过量可诱发人体消化系统的癌变等。增加人体摄入硝酸盐的潜在危害的因素一方面是我国水体中硝酸盐污染日益严重,部分省份地下水监测调研结果显示硝酸盐含量平均超标率达12.06-31.3%,硝酸盐最高含量达169mg/L,超出了世界卫生组织或欧盟应用水硝酸盐卫生标准。另一方面是我国居民消费较大的蔬菜均存在不同程度的硝酸盐污染现象。有研究表明,过量使用氮肥导致蔬菜中硝酸盐含量严重超标。人体摄入的硝酸盐中,72%-94%来自蔬菜,而叶菜是极易富集硝酸盐的蔬菜。有调查表明,青菜,空心菜和菠菜等绿叶蔬菜菜叶中硝酸盐含量普遍偏高。蔬菜硝酸盐含量的对身体健康的影响引发居民的关注意识,及时了解食品及饮用水中硝酸盐含量,对居民选择饮食种类具有一定的指导作用。
目前在食品中测定硝酸盐的方法有比色法、紫外分光光度法、离子色谱法、和离子专用电极法,还有定性的液体试剂法、固体试剂法和试剂盒法等。其中比色法、紫外分光光度法、离子色谱法、和离子专用电极法等需要特殊的仪器设备与专业检测人员,样品前处理步骤繁琐;试剂盒法需要使用化学试剂,通过颜色比对比得出硝酸盐在某一范围内的结果,以上方法均在检测普及上缺乏可操作性,制约了蔬菜和饮用水中硝酸盐在市场快速检测需求方面及普及性的应用。目前市场缺少不使用化学试剂、操作简单快速、现场安全、快速检测蔬菜中硝酸盐的检测设备和方法。
发明内容
为了克服现有技术的上述不足,本发明的目的是提供一种基于红外光谱的蔬菜和水体中硝酸盐快速检测的建模方法和检测方法。本发明提供一种利用基于硝酸盐中的氮氧键在中红外波段光谱中具有振动特征,在特定红外光源照射样品时,硝酸根分子中的N-O基团发生震动吸收,根据其特殊的频率6.5-8μm的位置产生特征吸收,获得硝酸盐的信息,再通过光谱吸收信息,通过自适应模型,实现对蔬菜和水体中硝酸盐的快速检测。测定方法快速、安全无污染、易于维护,有效提高蔬菜和水体中硝酸盐的检测效率,在农业和食品检测领域显示很好的应用潜力。
完成本申请发明任务的技术方案如下:
一种基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,包括如下步骤:
步骤一:购买市售不同种类绿色菜叶蔬菜,每种蔬菜分为两份,一份采用水杨酸-硫酸法测定,得到不同样品的硝酸盐的化学参考值,作为实验对照。另一份挤压取汁,将汁液置于样品采集池中,获得该类蔬菜在6.5-8μm波数的特征吸收光谱图,作为建模样本;
步骤二:配制一系列硝酸盐溶液不同浓度的标准溶液,并采集标准溶液的光谱信息,作为建模样本;
步骤三:将步骤一与步骤二所得曲线统一进进行光谱预处理,扣除干扰;
步骤四:将步骤三中处理好的市售蔬菜菜叶中硝酸盐光谱曲线及标准溶液光谱曲线与步骤二中通过化学分析法测得硝酸盐的化学参考值建模确立相关性;
步骤五:将待测蔬菜菜叶挤压取汁,将待测液置于样品采集池中,获得待测液的6.5-8μm波数的特征吸收光谱图;并将待测光谱进行预处理,扣除干扰;
步骤六:将待测蔬菜菜叶中硝酸盐光谱曲线与步骤四中所得模型智能鉴定分类,选择最佳匹配建模样本及建模样本数,通过智能算法和相关参数,获取建模样本中最具相似特性的样本,预测待测蔬菜菜叶中硝酸盐含量。
步骤七:将预测出的硝酸盐含量值与相关食用参考进行比对,获得蔬菜的食用信息。
更具体地说,以上步骤五至步骤七是一种基于筛选范围的红外全反射光谱的蔬菜和水体中硝酸盐现场安全、快速的检测方法,其特征在于,采用中红外量子级激光源,激发硝酸盐在6.5-8μm内产生振动,获得硝酸根的中红外特征吸收光谱。
其中,光谱采集过程中,连续扫描32次,扫描分辨率为4cm-1,取平均光谱,得到标准硝酸盐溶液的6.5-8μm波数的特征吸收光谱图,作为建模样本。
选择最佳匹配建模样本及建模样本数,通过欧氏距离计算方法获得,并通过模型参数如相关系数R2,均方根误差RMSE以及模型预测性能综合性评价指标RPD,建模样本数,通过最佳计算模型,预测待测蔬菜菜叶中硝酸盐含量。
以上方案中,步骤一所述的水杨酸-硫酸法测定,请参考:王学奎,《植物生理生化实验原理和技术》,高等教育出版社,北京,2015;
步骤三所述的光谱预处理,扣除干扰,请参考:杜昌文,《土壤红外光声光谱原理及应用》,科学出版社,北京,2012。
步骤五所述的光谱预处理方法同步骤三。
优选的,步骤一和步骤五中,由于硝酸盐富集在蔬菜菜叶部分,故采集待测蔬菜菜叶,并挤压取汁,新鲜蔬菜汁现场检测;
优选的,步骤一中,采集本方法的中红外衰减全反射光源为特定中红外量子级激光源,能够激发硝酸盐在6.5-8μm内产生振动,频率范围窄,获得光谱信噪比高,且重现性好;
优选的,步骤二中,配制一系列硝酸盐溶液不同浓度的标准溶液,推荐采用以下不同浓度的标准溶液:0mg/kg、500mg/kg、1000mg/kg、1500mg/kg、2000mg/kg、2500mg/kg、3000mg/kg、3500mg/kg、4000mg/kg、4500mg/kg的标准溶液。
优选的,步骤一、步骤二和步骤五中,扫描方法为连续扫描32次,扫描分辨率为4cm-1,取平均光谱;
优选的,计算过程通过Matlab软件进行代码编写,
原始光谱采用Savitzky-Golay(S-G)卷积平滑法,其语法结构为:
Spectrum%原始光谱;
[sm_data]=smooth(Spectrum,window,deg,type);
平滑预处理,其语法结构为:
[b,a]=butter(n,wn,’low’)%a,b为巴特沃兹滤波器参数;
ATRspectra=filtfilt(b,a,Spectrum);
扣除水分干扰,其语法结构为:
ref_water_ATR%纯水(硝酸盐含量为0时)的光谱曲线;
ATR(:,k)%不同种类蔬菜光谱曲线;
ATRnew=ATR(:,k)-factor*ref_water_ATR;
P=polyfit(wavenumbrt,ATRnew,1)。
优选的,当建模样本数据足够多及模型精确度达标后,可将建模样本数据预存于算法中。在预测待测蔬菜硝酸盐含量时可直接从步骤五开始进行。
优选的,步骤六中,将待测蔬菜菜叶中硝酸盐光谱曲线与步骤四中所得模型智能鉴定分类,通过欧氏距离计算方法,选择最佳匹配建模样本:
两样本之间的欧氏距离D计算公式如下:
其中Xis为建模样本的光谱数据,Xim为待测样本的光谱数据;将训练样本特征光谱和待测样本的光谱数据矩阵计算距离D,并根据光谱相似性按照从近到远排序,重新组合与待预测光谱最相近的训练样本矩阵。其语法结构为:
D=pdist(ATRnew,'euclidean');
优选的,步骤六中,模型参数如相关系数R2,均方根误差RMSE以及模型预测性能综合性评价指标RPD的相关参数如下:相关系数R2接近于1,均方根误差RMSE接近于0,模型预测性能综合性评价指标RPD大于2,通过最佳计算模型,预测待测蔬菜菜叶中硝酸盐含量。具体参数计算方程如下:
其中y和y'为个样本的化学参考值和偏最小二乘模型预测值,为样品化学参考值的均值,n为样品容量,SD为样本化学参考值的标准差。过校正集样本计算得到的均方根误差称为均方根校正误差(RMSEC);通过交叉验证计算得到的均方根误差称为交叉验证均方根误差(RMSECV);通过预测集样本计算得到的均方根误差称为均方根预测误差RMSEP)。在评价模型时,RMSE值越小而R2或R值越大,模型的预测能力越好;RPD被看作是评价模型预测能力的综合指标,其值越大说明模型越好。根据RDP值的范围,模型的评价结果如下
表1依据RPD值模型适用情况
本发明的创新点在于:
一、中红外量子级激光源,激发硝酸盐在6.5-8μm内产生振动,由于减少了光学元件,没有光栅或棱镜分光器,降低光损耗,使到达检测器的辐射强度大,信噪比高,且重现性好;
二、硝酸盐的光谱范围在6.5-8μm,频率范围窄,采集光谱数据只需数秒,大大提高的检测效率;
三、采集市售蔬菜获得硝酸盐含量化学参考值,并与光谱信息数字信号进行模型构建,内置智能模型,测定时无需再计算,可直接定量检测待测蔬菜叶中硝酸盐含量,结果更加接近市售蔬菜中硝酸盐含量的真实情况,结果更加可信。
本发明的有益效果在于:
一、本发明利用特定中红外激发光源获取在硝酸盐特征吸收波段的光谱吸收,简化光谱处理程序,所需样本量少,无需繁琐的预处理步骤,无需专业技术人员操作;
二、检测速度快、不使用化学试剂,整个过程仅需数分钟,将待测蔬菜挤压出少量汁液即可检测,适合现场、安全、快速获取,易于推广;
三、在农业生产和食品检测中,能够实现检测蔬菜中硝酸盐含量、检测水体中硝酸盐含量、以及蔬菜不同生长阶段硝酸盐含量变化,既能满足市民对蔬菜食品质量安全性的要求,也可根据硝酸盐诊断结果指导蔬菜种植过程中肥料的施加,符合市场需求,具有较大的社会效益。
具体实施方式
下面对于本发明所提出的一种基于筛选范围的红外全反射光谱的蔬菜和水体中硝酸
盐现场安全、快速的检测方法,结合具体实施例,作进一步详细阐述。以下实施例用于说明本发明,但并不用来限制本发明的范围。
若未特别指明,实施例中所有的技术手段为本领域技术人员所熟知的常规手段。
实施例1
一种基于筛选范围的红外全反射光谱的蔬菜和水体中硝酸盐现场安全、快速的检测方法,具体方式如下:
步骤一:购买市售不同种类绿色叶菜蔬菜(青菜、菠菜和小白菜)各30份,每种蔬菜分别用两种方法测定。一种采用水杨酸-硫酸法测定,参考(王学奎,《植物生理生化实验原理和技术》,高等教育出版社,北京,2015):选取1-2株新鲜蔬菜,擦净叶片表面污物,将叶片剪碎混匀,称取1-2g放入刻度试管中,加入10mL去离子水并用封口膜封口,置于沸水浴中提取30min后取出,冷取,将提取液过滤到25mL容量瓶中,反复冲洗残渣,最后定容至刻度。吸取样品液0.1mL分别于刻度试管中,然后加入5%水杨酸-硫酸溶液0.4mL,混匀后置室温下20min,再慢慢加入9.5mL 8%NaOH溶液,待冷却至室温后,以空白做参比,在410nm波长下测其吸光度。通过标准曲线的标准线性方程计算硝酸盐含量,得到硝酸盐的化学参考值。分别测得青菜中硝酸盐含量的平均值为3030mg/kg;菠菜中硝酸盐含量的平均值为2480mg/kg;小白菜中硝酸盐含量的平均值为900mg/kg;作为实验对照。另一份挤压取汁,将汁液置于样品采集池中,连续扫描32次,扫描分辨率为4cm-1,取平均光谱,分别得到蔬菜在红外光谱范围为6.5-8μm波数的特征吸收光谱图,作为建模样本;
步骤二:配制一系列硝酸盐溶液不同浓度的标准溶液,(0mg/kg、500mg/kg、1000mg/kg、1500mg/kg、2000mg/kg、2500mg/kg、3000mg/kg、3500mg/kg、4000mg/kg、4500mg/kg的标准溶液),并采集标准溶液的光谱信息,连续扫描32次,扫描分辨率为4cm-1,取平均光谱,得到标准硝酸盐溶液的6.5-8μm波数的特征吸收光谱图,作为建模样本;
步骤三:将步骤一与步骤二所得曲线统一进进行光谱预处理,扣除干扰(杜昌文,《土壤红外光声光谱原理及应用》,科学出版社,北京,2012);
通过Matlab软件进行代码编写,其语法结构为:
消噪、平滑预处理:
Spectrum%原始光谱
[b,a]=butter(2,0.01,’low’)%a,b为2阶截止频率为0.01的巴特沃兹滤波器
ATRspectra=filtfilt(b,a,Spectrum)
扣除水分干扰:
ref_water_ATR%纯水(硝酸盐含量为0时)的光谱曲线
ATR(:,k)%不同种类蔬菜光谱曲线
ATRnew=ATR(:,k)-factor*ref_water_ATR;
P=polyfit(wavenumbrt,ATRnew,1);
步骤四:将步骤三中处理好的市售蔬菜菜叶中硝酸盐光谱曲线及标准溶液光谱曲线与步骤二中通过化学分析法测得硝酸盐的化学参考值建模确立相关性;利用光谱数值与步骤一化学分析法检测的硝酸盐含量,采用偏最小二乘法建立样品光谱信息模型;
步骤五:将购买小白菜中另取出10份作为待测蔬菜,将待测小白菜菜叶挤压取汁,将待测液置于样品采集池中,连续扫描32次,扫描分辨率为4cm-1,取平均光谱,得到标准硝酸盐溶液的6.5-8μm波数的特征吸收光谱图;并将待测光谱进行预处理,扣除干扰,预处理方法同步骤三;
步骤六:将待测蔬菜菜叶中硝酸盐光谱曲线与步骤四中所得模型智能鉴定分类,通过欧氏距离计算方法,选出与待测小白菜中最相关的蔬菜菜叶的光谱图(步骤一中采集的光谱图),选择最佳匹配建模样本,该智能模型可优先选择出步骤一中小白菜的光谱曲线,基本不选择步骤一中菠菜与青菜的光谱曲线,并通过模型参数如相关系数R2,均方根误差RMSE以及模型预测性能综合性评价指标RPD的计算方程:
其中y和y'为个样本的化学参考值和偏最小二乘模型预测值,为样品化学参考值的均值,n为样品容量,SD为样本化学参考值的标准差。过校正集样本计算得到的均方根误差称为均方根校正误差(RMSEC);通过交叉验证计算得到的均方根误差称为交叉验证均方根误差(RMSECV);通过预测集样本计算得到的均方根误差称为均方根预测误差RMSEP)。
可得出相关结果:通过智能算法,获取建模样本中最具相似特性的建模样本数目为20个,相关系数R2为0.9724,均方根误差RMSE为20mg/kg,模型预测性能综合性评价指标为2.21,通过最佳计算模型,预测待测蔬菜菜叶中小白菜的硝酸盐含量912mg/kg。按照2002年世界卫生组织和联合国粮农组织(WHO)食品添加剂联合专家委员会JECFA第59次会议建议,参照居民日常小白菜的消费量,折合小白菜中成人每日硝酸盐容许摄入量约为456mg/人。
步骤七:将预测出的硝酸盐含量值与相关食用参考进行比对,最后根据小白菜硝酸盐含量判断食用信息。相关参考信息如下:
表2蔬菜食用参考
可得出结论,居民每人每日经蔬菜摄入的硝酸盐含量介于432-785mg/人.日的范围内,该类小白菜不可生食。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进接润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (9)
1.一种基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,步骤如下:
步骤一:购买市售不同种类绿色菜叶蔬菜,每种蔬菜分为两份,一份采用水杨酸-硫酸法测定,得到不同样品的硝酸盐的化学参考值,作为实验对照;另一份挤压取汁,将汁液置于样品采集池中,获得该类蔬菜在6.5-8μm波数的特征吸收光谱图,作为建模样本;
步骤二:配制一系列硝酸盐溶液不同浓度的标准溶液,并采集标准溶液的光谱信息,作为建模样本;
步骤三:将步骤一与步骤二所得曲线统一进进行光谱预处理,扣除干扰;
步骤四:将步骤三中处理好的市售蔬菜菜叶中硝酸盐光谱曲线及标准溶液光谱曲线与步骤二中通过化学分析法测得硝酸盐的化学参考值建模确立相关性;
步骤五:将待测蔬菜菜叶挤压取汁,将待测液置于样品采集池中,获得待测液的6.5-8μm波数的特征吸收光谱图;并将待测光谱进行预处理,扣除干扰;
步骤六:将待测蔬菜菜叶中硝酸盐光谱曲线与步骤四中所得模型智能鉴定分类,选择最佳匹配建模样本及建模样本数,通过智能算法和相关参数,获取建模样本中最具相似特性的样本,预测待测蔬菜菜叶中硝酸盐含量;
步骤七:将预测出的硝酸盐含量值与相关食用参考进行比对,获得蔬菜的食用信息。
2.根据权利要求1所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,所述步骤五至步骤七是一种基于筛选范围的红外全反射光谱的蔬菜和水体中硝酸盐现场安全、快速的检测方法:采用中红外量子级激光源,激发硝酸盐在6.5-8μm内产生振动,获得硝酸根的中红外特征吸收光谱。
3.根据权利要求1所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,其中,光谱采集过程中,连续扫描32次,扫描分辨率为4cm-1,取平均光谱,得到标准硝酸盐溶液的6.5-8μm波数的特征吸收光谱图,作为建模样本。
4.根据权利要求1所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,选择最佳匹配建模样本及建模样本数,通过欧氏距离计算方法获得,并通过模型参数如相关系数R2,均方根误差RMSE以及模型预测性能综合性评价指标RPD,建模样本数,通过最佳计算模型,预测待测蔬菜菜叶中硝酸盐含量。
5.根据权利要求1所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,步骤一和步骤五中,由于硝酸盐富集在蔬菜菜叶部分,故采集待测蔬菜菜叶,并挤压取汁,新鲜蔬菜汁现场检测。
6.根据权利要求1所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,步骤二中,配制一系列硝酸盐溶液不同浓度的标准溶液,推荐采用以下不同浓度的标准溶液:0mg/kg、500mg/kg、1000mg/kg、1500mg/kg、2000mg/kg、2500mg/kg、3000mg/kg、3500mg/kg、4000mg/kg、4500mg/kg的标准溶液。
7.根据权利要求1-6之一所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,计算过程通过Matlab软件进行代码编写,
原始光谱采用Savitzky-Golay卷积平滑法,其语法结构为:
Spectrum%原始光谱;
[sm_data]=smooth(Spectrum,window,deg,type);
平滑预处理,其语法结构为:
[b,a]=butter(n,wn,’low’)%a,b为巴特沃兹滤波器参数;
ATRspectra=filtfilt(b,a,Spectrum);
扣除水分干扰,其语法结构为:
ref_water_ATR%纯水的光谱曲线;
ATR(:,k)%不同种类蔬菜光谱曲线;
ATRnew=ATR(:,k)-factor*ref_water_ATR;
P=polyfit(wavenumbrt,ATRnew,1)。
8.根据权利要求7所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,将待测蔬菜菜叶中硝酸盐光谱曲线与步骤四中所得模型智能鉴定分类,通过欧氏距离计算方法,选择最佳匹配建模样本:
两样本之间的欧氏距离D计算公式如下:
其中Xis为建模样本的光谱数据,Xim为待测样本的光谱数据;将训练样本特征光谱和待测样本的光谱数据矩阵计算距离D,并根据光谱相似性按照从近到远排序,重新组合与待预测光谱最相近的训练样本矩阵;其语法结构为:
D=pdist(ATRnew,'euclidean')。
9.根据权利要求8所述的基于筛选范围的红外全反射光谱的蔬菜硝酸盐现场安全、快速检测的建模方法和检测方法,其特征在于,步骤六中,模型参数如相关系数R2,均方根误差RMSE以及模型预测性能综合性评价指标RPD的相关参数如下:相关系数R2接近于1,均方根误差RMSE接近于0,模型预测性能综合性评价指标RPD大于2,通过最佳计算模型获得智能模型,预测待测蔬菜菜叶中硝酸盐含量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710249286.0A CN107024445A (zh) | 2017-04-17 | 2017-04-17 | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710249286.0A CN107024445A (zh) | 2017-04-17 | 2017-04-17 | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107024445A true CN107024445A (zh) | 2017-08-08 |
Family
ID=59526398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710249286.0A Pending CN107024445A (zh) | 2017-04-17 | 2017-04-17 | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107024445A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505283A (zh) * | 2017-09-04 | 2017-12-22 | 上海交通大学 | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 |
CN110887792A (zh) * | 2019-12-17 | 2020-03-17 | 深圳慧格科技服务咨询有限公司 | 一种基于无人机高光谱影像的河流水质快速监测系统 |
CN111239071A (zh) * | 2020-02-19 | 2020-06-05 | 中国科学院烟台海岸带研究所 | 一种光谱法海水硝酸盐浓度检测方法 |
CN112750507A (zh) * | 2021-01-15 | 2021-05-04 | 中南大学 | 基于混合机器学习模型的同时检测水中硝酸盐和亚硝酸盐含量的方法 |
CN113092407A (zh) * | 2021-04-20 | 2021-07-09 | 深圳市国科光谱技术有限公司 | 基于近红外光谱的五常稻花香水稻原产地品种识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710071A (zh) * | 2009-12-09 | 2010-05-19 | 川渝中烟工业公司 | 建立烟草近红外模型的选样方法 |
CN105445217A (zh) * | 2015-07-27 | 2016-03-30 | 南京财经大学 | 基于衰减全反射傅里叶变换红外光谱技术的糙米中黄曲霉毒素含量的快速检测方法 |
CN105445218A (zh) * | 2015-11-25 | 2016-03-30 | 中国科学院南京土壤研究所 | 中红外光谱油菜籽蛋白质含量检测自适应模型的建立方法 |
CN105486655A (zh) * | 2015-11-23 | 2016-04-13 | 中国科学院南京土壤研究所 | 基于红外光谱智能鉴定模型的土壤有机质快速检测方法 |
CN105866056A (zh) * | 2015-03-25 | 2016-08-17 | 山东翰能高科科技有限公司 | 一种基于近红外光谱的杂交种纯度鉴别方法 |
CN106525761A (zh) * | 2016-11-08 | 2017-03-22 | 浙江大学 | 基于太赫兹光谱扫描的亚硝酸盐检测方法 |
-
2017
- 2017-04-17 CN CN201710249286.0A patent/CN107024445A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710071A (zh) * | 2009-12-09 | 2010-05-19 | 川渝中烟工业公司 | 建立烟草近红外模型的选样方法 |
CN105866056A (zh) * | 2015-03-25 | 2016-08-17 | 山东翰能高科科技有限公司 | 一种基于近红外光谱的杂交种纯度鉴别方法 |
CN105445217A (zh) * | 2015-07-27 | 2016-03-30 | 南京财经大学 | 基于衰减全反射傅里叶变换红外光谱技术的糙米中黄曲霉毒素含量的快速检测方法 |
CN105486655A (zh) * | 2015-11-23 | 2016-04-13 | 中国科学院南京土壤研究所 | 基于红外光谱智能鉴定模型的土壤有机质快速检测方法 |
CN105445218A (zh) * | 2015-11-25 | 2016-03-30 | 中国科学院南京土壤研究所 | 中红外光谱油菜籽蛋白质含量检测自适应模型的建立方法 |
CN106525761A (zh) * | 2016-11-08 | 2017-03-22 | 浙江大学 | 基于太赫兹光谱扫描的亚硝酸盐检测方法 |
Non-Patent Citations (1)
Title |
---|
杨家宝 等: ""基于傅里叶变换中红外光谱的小白菜硝酸盐含量的快速测定"", 《分析化学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107505283A (zh) * | 2017-09-04 | 2017-12-22 | 上海交通大学 | 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法 |
CN110887792A (zh) * | 2019-12-17 | 2020-03-17 | 深圳慧格科技服务咨询有限公司 | 一种基于无人机高光谱影像的河流水质快速监测系统 |
CN111239071A (zh) * | 2020-02-19 | 2020-06-05 | 中国科学院烟台海岸带研究所 | 一种光谱法海水硝酸盐浓度检测方法 |
CN111239071B (zh) * | 2020-02-19 | 2023-03-31 | 中国科学院烟台海岸带研究所 | 一种光谱法海水硝酸盐浓度检测方法 |
CN112750507A (zh) * | 2021-01-15 | 2021-05-04 | 中南大学 | 基于混合机器学习模型的同时检测水中硝酸盐和亚硝酸盐含量的方法 |
CN112750507B (zh) * | 2021-01-15 | 2023-12-22 | 中南大学 | 基于混合机器学习模型的同时检测水中硝酸盐和亚硝酸盐含量的方法 |
CN113092407A (zh) * | 2021-04-20 | 2021-07-09 | 深圳市国科光谱技术有限公司 | 基于近红外光谱的五常稻花香水稻原产地品种识别方法 |
CN113092407B (zh) * | 2021-04-20 | 2024-01-30 | 深圳华创智信科技有限公司 | 基于近红外光谱的五常稻花香水稻原产地品种识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107024445A (zh) | 蔬菜中硝酸盐的快速检测的建模方法和检测方法 | |
CN101915744B (zh) | 物质成分含量的近红外光谱无损检测方法及装置 | |
Xu et al. | An overview on nondestructive spectroscopic techniques for lipid and lipid oxidation analysis in fish and fish products | |
CN108875913B (zh) | 一种基于卷积神经网络的松茸快速无损检测系统及方法 | |
Cozzolino et al. | Quantitative analysis of minerals and electric conductivity of red grape homogenates by near infrared reflectance spectroscopy | |
CN104048941B (zh) | 采用近红外光谱快速测定麦冬中多指标性成分含量的方法 | |
CN102788752A (zh) | 基于光谱技术的农作物内部信息无损检测装置及方法 | |
CN107014769B (zh) | 一种基于THz光谱分析的鲜肉K值快速无损检测模型 | |
Varrà et al. | Multi-element signature of cuttlefish and its potential for the discrimination of different geographical provenances and traceability | |
CN103954734B (zh) | 一种适宜果汁加工的荔枝品质的测定和评价方法 | |
CN108181262A (zh) | 一种利用近红外光谱技术快速测定铜藻纤维素含量的方法 | |
Song et al. | Spatial variability analysis of within-field winter wheat nitrogen and grain quality using canopy fluorescence sensor measurements | |
CN108956584B (zh) | 一种桑椹中重金属元素铬的快速准确检测方法 | |
CN105136737A (zh) | 基于近红外光谱的馒头中马铃薯全粉含量快速测定的方法 | |
CN108519339A (zh) | 一种基于WT-LSSVR的叶片镉含量Vis-NIR光谱特征建模方法 | |
CN107505179A (zh) | 一种土壤预处理与养分近红外光谱检测方法 | |
Yang et al. | Detection of total viable count in spiced beef using hyperspectral imaging combined with wavelet transform and multiway partial least squares algorithm | |
CN104655585B (zh) | 一种基于近红外光谱的pse肉筛选方法 | |
CN109060771A (zh) | 一种基于光谱不同特征集的共识模型构建方法 | |
Chen et al. | Advanced nondestructive detection technologies in food | |
Zhang et al. | Evaluating maize grain quality by continuous wavelet analysis under normal and lodging circumstances | |
CN208476785U (zh) | 一种适用于水产品重金属检测的激光诱导击穿光谱检测装置 | |
Shi et al. | Nondestructive diagnostics of magnesium deficiency based on distribution features of chlorophyll concentrations map on cucumber leaf | |
CN106960091A (zh) | 一种鲜肉新鲜度K值的THz光谱快速无损检测方法 | |
CN108593596A (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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Du Changwen Inventor after: Ma Fei Inventor before: Du Changwen Inventor before: Ma Fei |
|
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170808 |