CN108898509A - 一种大米溯源模型的建立方法及其应用 - Google Patents
一种大米溯源模型的建立方法及其应用 Download PDFInfo
- Publication number
- CN108898509A CN108898509A CN201810758482.5A CN201810758482A CN108898509A CN 108898509 A CN108898509 A CN 108898509A CN 201810758482 A CN201810758482 A CN 201810758482A CN 108898509 A CN108898509 A CN 108898509A
- Authority
- CN
- China
- Prior art keywords
- rice
- layer
- output
- building
- hidden layer
- 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
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 130
- 235000009566 rice Nutrition 0.000 title claims abstract description 130
- 238000000034 method Methods 0.000 title claims abstract description 38
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 129
- 238000004519 manufacturing process Methods 0.000 claims abstract description 30
- 239000002689 soil Substances 0.000 claims abstract description 22
- 230000001427 coherent effect Effects 0.000 claims abstract description 21
- 239000000470 constituent Substances 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 15
- 239000000463 material Substances 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 238000000120 microwave digestion Methods 0.000 claims description 21
- 229910052712 strontium Inorganic materials 0.000 claims description 15
- 229910052725 zinc Inorganic materials 0.000 claims description 15
- 229910052689 Holmium Inorganic materials 0.000 claims description 13
- 229910052787 antimony Inorganic materials 0.000 claims description 13
- 229910052785 arsenic Inorganic materials 0.000 claims description 13
- 229910052793 cadmium Inorganic materials 0.000 claims description 13
- 229910052791 calcium Inorganic materials 0.000 claims description 13
- 229910052804 chromium Inorganic materials 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 13
- 229910052746 lanthanum Inorganic materials 0.000 claims description 13
- 229910052745 lead Inorganic materials 0.000 claims description 13
- 229910052748 manganese Inorganic materials 0.000 claims description 13
- 229910052709 silver Inorganic materials 0.000 claims description 13
- 229910052688 Gadolinium Inorganic materials 0.000 claims description 12
- 229910052779 Neodymium Inorganic materials 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 12
- 229910052692 Dysprosium Inorganic materials 0.000 claims description 11
- 229910052772 Samarium Inorganic materials 0.000 claims description 11
- 229910052749 magnesium Inorganic materials 0.000 claims description 11
- 229910052691 Erbium Inorganic materials 0.000 claims description 10
- 229910052770 Uranium Inorganic materials 0.000 claims description 10
- 229910052769 Ytterbium Inorganic materials 0.000 claims description 10
- 229910052701 rubidium Inorganic materials 0.000 claims description 10
- 229910052802 copper Inorganic materials 0.000 claims description 9
- 229910052759 nickel Inorganic materials 0.000 claims description 9
- 229910052693 Europium Inorganic materials 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 8
- 229910052788 barium Inorganic materials 0.000 claims description 8
- 239000007788 liquid Substances 0.000 claims description 8
- 229910052777 Praseodymium Inorganic materials 0.000 claims description 7
- 235000013312 flour Nutrition 0.000 claims description 7
- 229910052741 iridium Inorganic materials 0.000 claims description 7
- 229910052697 platinum Inorganic materials 0.000 claims description 7
- 229910052703 rhodium Inorganic materials 0.000 claims description 7
- GRYLNZFGIOXLOG-UHFFFAOYSA-N Nitric acid Chemical compound O[N+]([O-])=O GRYLNZFGIOXLOG-UHFFFAOYSA-N 0.000 claims description 6
- 229910052771 Terbium Inorganic materials 0.000 claims description 6
- 229910052782 aluminium Inorganic materials 0.000 claims description 6
- 229910017604 nitric acid Inorganic materials 0.000 claims description 6
- 229910052700 potassium Inorganic materials 0.000 claims description 6
- 229910052708 sodium Inorganic materials 0.000 claims description 6
- 229910052684 Cerium Inorganic materials 0.000 claims description 5
- 229910052776 Thorium Inorganic materials 0.000 claims description 5
- 229910052775 Thulium Inorganic materials 0.000 claims description 5
- 238000001816 cooling Methods 0.000 claims description 5
- 239000007789 gas Substances 0.000 claims description 5
- 229910052742 iron Inorganic materials 0.000 claims description 5
- 229910052750 molybdenum Inorganic materials 0.000 claims description 5
- 229910052763 palladium Inorganic materials 0.000 claims description 5
- 238000000918 plasma mass spectrometry Methods 0.000 claims description 5
- 229910052707 ruthenium Inorganic materials 0.000 claims description 5
- 229910052706 scandium Inorganic materials 0.000 claims description 5
- 229910052714 tellurium Inorganic materials 0.000 claims description 5
- 229910052718 tin Inorganic materials 0.000 claims description 5
- 229910052720 vanadium Inorganic materials 0.000 claims description 5
- 229910052727 yttrium Inorganic materials 0.000 claims description 5
- 239000012159 carrier gas Substances 0.000 claims description 4
- 239000000498 cooling water Substances 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 239000011259 mixed solution Substances 0.000 claims description 2
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims 1
- 230000008878 coupling Effects 0.000 claims 1
- 238000010168 coupling process Methods 0.000 claims 1
- 238000005859 coupling reaction Methods 0.000 claims 1
- 229910001882 dioxygen Inorganic materials 0.000 claims 1
- 230000001939 inductive effect Effects 0.000 claims 1
- 150000002500 ions Chemical class 0.000 claims 1
- 238000001819 mass spectrum Methods 0.000 claims 1
- 238000010219 correlation analysis Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 38
- 238000013528 artificial neural network Methods 0.000 description 11
- MHAJPDPJQMAIIY-UHFFFAOYSA-N Hydrogen peroxide Chemical compound OO MHAJPDPJQMAIIY-UHFFFAOYSA-N 0.000 description 9
- 239000002366 mineral element Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 5
- 235000013339 cereals Nutrition 0.000 description 4
- 210000002569 neuron Anatomy 0.000 description 4
- 229910052711 selenium Inorganic materials 0.000 description 4
- 239000004677 Nylon Substances 0.000 description 3
- 239000002253 acid Substances 0.000 description 3
- 229920001778 nylon Polymers 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 229910052797 bismuth Inorganic materials 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 239000000428 dust Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000013305 food Nutrition 0.000 description 2
- 229910052732 germanium Inorganic materials 0.000 description 2
- 229910052738 indium Inorganic materials 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000009423 ventilation Methods 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 230000003698 anagen phase Effects 0.000 description 1
- 239000012496 blank sample Substances 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005034 decoration Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 239000010903 husk Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 235000012054 meals Nutrition 0.000 description 1
- 238000003801 milling Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 235000013348 organic food Nutrition 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000010298 pulverizing process Methods 0.000 description 1
- 238000003044 randomized block design Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 229910052716 thallium Inorganic materials 0.000 description 1
- 229910021642 ultra pure water Inorganic materials 0.000 description 1
- 239000012498 ultrapure water Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/66—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light electrically excited, e.g. electroluminescence
- G01N21/68—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light electrically excited, e.g. electroluminescence using high frequency electric fields
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Business, Economics & Management (AREA)
- Molecular Biology (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Analytical Chemistry (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Immunology (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种大米溯源模型的建立方法及其应用，属于大米溯源技术领域。本发明所述建立方法包括以下步骤：1)收集稻谷及稻谷所在地的土壤母质样品，得检测样本；2)测定检测样本中的元素含量，经方差变异贡献率分析和相关性分析，筛选出相关元素；3)利用筛选出的所述相关元素建立大米溯源的前馈神经网络模型。利用本发明的溯源模型进行大米产地识别，结果平均相对误差值为17.14％，产地的整体识别准确率为100％，可见本发明提供的模型具有自适应性和稳定快速收敛特性，可作为大米产地的智能识别模型。
Description
技术领域
本发明属于大米溯源技术领域，具体涉及一种大米溯源模型的建立方法及其应用。
背景技术
水稻是我国的重要粮食之一，稻谷播种面积逐年增加，2011年我国水稻播种面积已突破30000千公顷，水稻产量连续十连增，国家统计局数字表明2012年我国稻谷产量达20,423.59万吨。全国约有19个省市以大米为主食，生产的水稻中85％作为口粮消费，因此水稻在我国具有很重要的位置。黑龙江水稻目前进入一个优质、高效、专业的新时期。黑龙江省地处纬度最北边的稻作区，具有光照充足、昼夜温差大，土质肥沃，水质优，污染少等条件，有利于发展水稻生产。由于产量稳定，米质优良，商品率高，黑龙江大米已成为我国重要的优质粳米生产基地，产品远销全国各地。2015年，黑龙江省水稻种植面积占全国水稻种植面积的17.0％，占东北三省水稻种植面积的69.25％，是东北最大的种植省份。黑龙江水稻产量占全国水稻产量的16.3％，占东北三省水稻产量的67.6％，在大米市场上具有重要的地位和作用。
黑龙江经国国家质量监督检验检疫总局批准的地理标志大米产品有9种，分别是五常大米、方正大米、响水大米、珍宝岛大米、建三江大米、泰来大米、通河大米、庆安大米和富拉尔基温水大米等。截止2016年底，全国经国国家质量监督检验检疫总局批准的地理标志大米产品有44种，其中东北三省的地理标志大米占40％以上，黑龙江省地理标志大米是全国地理标志大米产品中最多的省份。
黑龙江地理标志大米如五常大米、查哈阳大米、响水大米、方正大米、建三江大米等等，因其品质好、口感佳，在国内和地区间流动性较大，且销售价格占有一定优势。随着市场需求量的逐年增加，大米行业市场出现鱼龙混杂，良莠不齐，普遍存在地理标志产品、绿色食品和有机食品错误标识，假冒黑龙江省产地品牌大米销售现象，使消费者权益受到严重损害。
目前，并没有关于如何更好的开展水稻产地溯源技术，建立水稻指纹数据库的记载。
发明内容
有鉴于此，本发明的目的在于提供一种大米溯源模型，解决其它影响因素对大米溯源的影响问题，具有自适应分辨性和良好的容错能力，鲁棒性高且精确达100％。
为了实现上述发明目的，本发明提供以下技术方案：
本发明提供了一种大米溯源模型的建立方法，包括以下步骤：
1)收集稻谷及稻谷所在地的土壤母质样品，得检测样本；
2)测定步骤1)得到的所述检测样本中的元素含量，经方差变异贡献率分析和相关性分析，筛选出相关元素；所述元素包括Na、Mg、Al、K、Ca、Sc、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Mo、Ru、Rh、Pd、Ag、Cd、Sn、Sb、Te、Cs、Ba、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu、Hf、Ir、Pt、Au、Tl、Pb、Th和U；
3)利用步骤2)筛选出的所述相关元素建立大米溯源的前馈神经网络模型，包括输入层、隐含层和输出层，层与层之间为互连方式，每层分别有nI，nh和no节点数，第一层为输入层，所述输入层的节点数为筛选出的相关元素的个数，设网络输入向量为式中T为向量维度或指标个数，只是传递信号到隐含层，令为输入层第i节点的输出：
第二层为第一隐含层，令输入层到第一隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第三层为第二隐含层，令输入层到第二隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第四层为输出层，令第二隐含层到输出层的连接权为输出层的激励函数p，阈值为输出向量为其中则为输出层的第k节点输出：
优选的，步骤1)所述样品稻谷的水分含量为14％以下。
优选的，步骤1)所述粉碎为利用旋风磨粉碎，并过80～120目筛。
优选的，步骤2)所述测定前还包括微波消解过程，所述微波消解包括：将所述米粉与微波消解液按照1g：(15～35)mL的质量体积比混合，置于微波消解仪中进行微波消解。
优选的，所述微波消解液为包括浓硝酸和双氧水的混合溶液。
优选的，所述微波消解的微波程序为8min内从0w增到1600w，温度升到120℃，保持2min；在5min内从温度120℃升到160℃，保持5min；在5min内再从160℃升到180℃，并在此温度下消解15min；冷却。
优选的，步骤2)所述元素含量利用电感耦合等离子体质谱进行测定。
优选的，所述电感耦合等离子体质谱的工作参数为射频功率1280w，雾化室温度2℃，冷却水流量1.47L/min，载气流量1.0L/min，补偿气体流量1.0L/min。
优选的，当鉴别大米是否来源于建三江、五常或查哈阳地区时，步骤2)所述相关元素包括Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Te、Ba、La、Nd、Sm、Gd、Dy、Ho、Er、Yb、Pb和U。
本发明还提供了一种利用上述建立方法建立的大米溯源模型在大米溯源识别中的应用，设共有M样本，第p个样本，对应的期望实际输出为dp，计算网络输出为yp，网络各个输出总误差为
若记wjk数，即连接权值和阈值(θ和)，则学习规则为
式中：η为学习速度；α为惯性系数；t为学习次数；
对得到的所述期望实际输出dp和所述计算网络输出yp进行比对，从而得到相应的稻谷产地。
本发明提供了一种大米溯源模型的建立方法，利用前馈神经网络采用四层网络训练，得到大米溯源模型，将筛选得到的元素个数(ni)作为输入层节点数，将筛选元素的含量(x)输入网络。利用本发明的溯源模型进行黑龙江建三江、五常和查哈阳等地区的大米产地识别，结果平均相对误差值为17.14％，产地的整体识别准确率为100％，可见前馈神经网络利用其具有自适应性和稳定快速收敛特性，可作为大米产地的智能识别模型。
附图说明
图1为大米溯源的前馈神经网络模型；
图2为大米产地分类网络收敛曲线。
具体实施方式
本发明提供了一种大米溯源模型的建立方法，包括以下步骤：
1)收集稻谷及稻谷所在地的土壤母质样品，得检测样本；
2)测定步骤1)得到的所述检测样本中的元素含量，经方差变异贡献率分析和相关性分析，筛选出相关元素；所述元素包括Na、Mg、Al、K、Ca、Sc、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Mo、Ru、Rh、Pd、Ag、Cd、Sn、Sb、Te、Cs、Ba、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu、Hf、Ir、Pt、Au、Tl、Pb、Th和U；
3)利用步骤2)筛选出的所述相关元素建立大米溯源的前馈神经网络模型，包括输入层、隐含层和输出层，层与层之间为互连方式，每层分别有nI，nh和no节点数，第一层为输入层，所述输入层的节点数为筛选出的相关元素的个数，设网络输入向量为x，式中只是传递信号到隐含层，T为向量维度或指标个数，令为输入层第i节点的输出：
第二层为第一隐含层，令输入层到第一隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第三层为第二隐含层，令输入层到第二隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第四层为输出层，令第二隐含层到输出层的连接权为输出层的激励函数p，阈值为输出向量为其中则为输出层的第k节点输出：
本发明所述大米溯源模型的建立方法，首先收集稻谷及稻谷所在地的土壤母质样品，得检测样本。本发明所述稻谷的收集，优选采用代表性采样原则，更优选采用棋盘式采样法，具体包括：依照不同保护范围各个区域，依照种植范围大小设置采样点，每块地块随机设置分点5个重复，每个采集点沿植株中部割取稻穗1m2左右，每个小区收集1～2kg稻穗。本发明实施例中，所述土壤母质样本优选收集自五常、建三江管理局和查哈阳农场，各选3块试验田块，设置3次重复，随机区组设计，小区面积不少于10m2，四周设保护行，保护行品种与各对应品种相同。按照当地大田水稻的管理方式，统一农事管理，本发明对所述管理方式并没有特殊限定。本发明所述五常、建三江管理局和查哈阳农场的地理信息和水稻生长期天气条件信息如表1所示：
表1.试验地地理信息和水稻生长期天气条件信息表
得到检测样本后，本发明测定所述检测样本中的元素含量，经方差变异贡献率分析和相关性分析后，筛选出相关元素；所述元素包括Na、Mg、Al、K、Ca、Sc、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Mo、Ru、Rh、Pd、Ag、Cd、Sn、Sb、Te、Cs、Ba、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu、Hf、Ir、Pt、Au、Tl、Pb、Th和U。本发明当样本为稻谷时，在所述检测前，优选包括脱粒后粉碎后得米粉，再微波消解的过程。本发明收集样品稻谷后，优选的还包括晾晒过程，使所述样品稻谷的水分含量在14％以下。本发明对所述晾晒的晾晒场地并没有特殊要求，优选无扬尘、整洁和避光的晾晒场地。本发明对晾晒完成的样品稻谷进行脱粒，优选的还包括砻谷和碾米过程，获得二级精米。本发明对所述二级精米进行粉碎，所述粉碎为利用旋风磨粉碎，并过80～120目筛。本发明对所述筛的类型并没有特殊限定，优选为尼龙筛，所述尼龙筛的筛孔优选为100目。本发明所述微波消解的步骤优选包括：将所述米粉与微波消解液按照1g：(15～35)mL的质量体积比混合，置于微波消解仪中进行微波消解。本发明所述米粉与微波消解液的质量体积比优选为1g：(20～34)mL，更优选为1g：(25～33)mL，最优选为1g：32mL。本发明所述所述微波消解液优选为包括浓硝酸和双氧水的混合溶液，所述浓硝酸的质量浓度优选为65～80％，更优选为68～75％，最优选为70％。本发明所述双氧水的质量浓度优选为20～45％，更优选为25～38％，最优选为30％。在本发明所述微波消解液中，所述浓硝酸和双氧水的体积比优选为(2～5)：1，更优选为(2.5～4)：1，最优选为3:1。本发明所述微波消解的微波程序优选为8min内从0w增到1600w，温度升到120℃，保持2min；在5min内从温度120℃升到160℃，保持5min；在5min内再从160℃升到180℃，并在此温度下消解15min；冷却。本发明所述冷却的时间优选为15～30min，更优选为18～25min，最优选为20min。本发明在所述冷却后，优选的还包括赶酸的过程，具体包括于通风橱内将微波消解后的溶液置于精确控温电热消解器中进行赶酸，并利用超纯水(>18.2MΩ·cm)洗涤样品，定容至100.00g进行元素含量的测定。本发明对所述元素含量的测定方法并没有特殊限定，优选利用电感耦合等离子体质谱(ICP-MS)进行测定。本发明所述电感耦合等离子体质谱的工作参数优选设置为射频功率1280w，雾化室温度2℃，冷却水流量1.47L/min，载气流量1.0L/min，补偿气体流量1.0L/min。本发明在测定元素含量时，优选利用Ge、In和Bi三种元素作为内标元素，保证仪器的稳定性。当内标元素的RSD＞5％时需重新测定，且要求米粉中元素的回收率均大于90％。本发明对所述土壤母质样本中的元素含量检测方法优选与稻谷样本的检测方法相同。
本发明优选利用SPSS软件计算各样品元素间的相关系数，选定相关元素，当鉴别大米是否来源于建三江、五常或查哈阳地区时，本发明所述相关元素优选为Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Te、Ba、La、Nd、Sm、Gd、Dy、Ho、Er、Yb、Pb和U。
得到建三江、五常或查哈阳地区的相关元素后，本发明利用相关元素建立大米溯源的前馈神经网络模型，包括输入层、隐含层和输出层，层与层之间为互连方式，每层分别有nI，nh和no节点数，第一层为输入层，所述输入层的节点数为筛选出的相关元素的个数，设网络输入向量为式中T为向量维度或指标个数，只是传递信号到隐含层，令为输入层第i节点的输出：
第二层为第一隐含层，令输入层到第一隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第三层为第二隐含层，令输入层到第二隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第四层为输出层，令第二隐含层到输出层的连接权为输出层的激励函数p，阈值为输出向量为其中则为输出层的第k节点输出：
在本发明实施例中，对三类大米产地样本中矿物元素数据采集，由于提取了23个相关元素，即网络输入层节点数为23；研究中将类型分类的编码为100、010和001，分别表示查哈阳大米、五常大米和建三江大米，因此，输出节点数为3；根据Kolmogorov定理和大量实际训练经验，确定第一隐含层神经元节点为36个，第二隐含层神经元节点为28个；所以基于神经网络的大米产地识别分类的神经网络的拓扑结构为23-36-28-3型，具体模型如图1所示，并根据输出层的实际输出值与理论输出值(类型分类的编码，如100、010和001)进行比对，从而确定各样品稻谷的来源。
本发明还提供了一种利用上述建立方法建立的大米溯源模型在大米溯源识别中的应用，设共有M样本，第p个样本，对应的期望实际输出为dp，计算网络输出为yp，网络各个输出总误差为
若记wjk数，即连接权值和阈值(θ和)，则学习规则为
式中：η为学习速度；α为惯性系数；t为学习次数；
对得到的所述期望实际输出dp和所述计算网络输出yp进行比对，从而得到相应的稻谷产地。
下面结合实施例对本发明提供的大米溯源模型进行详细的说明，但是不能把它们理解为对本发明保护范围的限定。
实施例1
分别测定建三江、五常和查哈阳地区的大米样本和土壤母质样本中的矿质元素含量。选择建三江、五常和查哈阳地区的稻米在通风处晾晒至水分14％以下，要求晾晒场地无扬尘、整洁、避光。对稻穗进行脱粒、砻谷、碾米获得二级精米，二级精米参照GB/T 5502-2008染色法判定大米加工精度。并对精米旋风磨粉碎处理，过100目尼龙筛，得到米粉样本。
准确称取0.25g的大米粉样品，置于消化管中加入6mL浓硝酸(质量浓度70％，BV3级)和2mL双氧水(质量浓度30％，BV3级)，放入MARS微波消解仪中进行消解。微波设置程序为8min内从0w增到1600w，温度升到120℃，保持2min；在5min内从温度120℃升到了160℃，保5min；在5min内再从160℃升到180℃，并在此温度下消解15min；然后冷却20min，将微波消化管取出，于通风橱内打开塞子将微波消化管置于精确控温电热消解器中进行赶酸。超纯水(>18.2MΩ·cm)洗涤样品，定容至100.00g，采用同样方法进行空白样品和大米标准物样品消解。
ICP-MS测定各元素含量，选用Ge、In和Bi三种元素作为内标元素，设置工作参数为射频功率1280w，雾化室温度2℃，冷却水流量1.47L/min，载气流量1.0L/min，补偿气体流量1.0L/min，仪器测定样品和大米标准物中Na、Mg、Al、K、Ca、Sc、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Mo、Ru、Rh、Pd、Ag、Cd、Sn、Sb、Te、Cs、Ba、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu、Hf、Ir、Pt、Au、Tl、Pb、Th、U52种元素，要求大米标准物中元素的回收率均大于90％。
利用所述ICP-MS仪器测定多种矿物元素的含量如表1所示：
表1.不同产地矿物元素的含量
注：*元素单位为mg/kg，其余均为ug/kg。
对于不同产地试验田来源大米样品中矿物元素含量进行方差分析，结果显示，Na、K、Mn、Ni、Cu、Zn、As、Se、Rb、Sr、Mo、Cd、Sb、Te、Ba、Pr、Nd、Eu、Gd、Pb和Rh元素的含量在不同产地之间有极显著差异(P<0.01)，Ca、Al、Co、Ag、La和Ho元素在不同产地之间有显著差异(P<0.05)，K受农业施肥管理措施影响，不予分析。
通过SPSS软件一般线性模型实现的多变量分析，即主效应和交互效应的方差分析以及分析产地对各元素含量变异的影响。结果如表2所示：
表2.产地因素对元素含量变异的方差贡献率(％)
结果显示，受产地影响较大的元素有Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Ba、La、Sm、Dy、Ho、Er、Pb、U含量变异贡献率最大。产地因素对元素Na、Ca、Mn、Ni、Cu、Zn、As、Se、Rb、Sr、Mo、Cd、Sb、Te、Ba、Pr、Nd、Eu、Gd、Pb和Rh含量有极显著影响(P<0.01)，对元素Al、Co、Ag、La、Er和Ho含量有显著影响(P<0.05)。
通过以上建立不同产地研究，初步筛选到受产地影响较大的Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Ba、La、Sm、Dy、Ho、Er、Pb、U19种元素。
利用相同方法测量三个地区的母质土壤中的元素含量，测得结果见表3：
表3.不同产地母质土壤中的矿物元素含量统计分析
结果显示：19种元素在不同地区间存在差异，Cr、Ni、Cu、Zn、Se、Rb、Sr、Mo、Ba、La和Pb在不同地区母质土壤中存在显著差异(P＜0.05)，Eu、Tb、Ag、Cd、Nd、Ir、Pt和Tl在不同地区母质土壤中存在极显著差异(P＜0.01)。其中，不同的地区母质土壤有不同的元素含量差异特点。Cr、Ni、Cu、Zn、Sr、Ag、Cd元素含量在不同产地中含量高低依次为查哈阳＞五常＞建三江。Se、Rb、Ba、La、Eu、Tl和Pb元素含量高低依次为查哈阳＞建三江＞五常。五常地区母质土壤中Nd含量最高，建三江地区母质土壤中Mo、Tb、Ir、Pt含量最高。
利用SPSS相关性分析，分析大米和不同产区土壤母质的相关性，结果如表4所示：
表4.母质土壤和大米之间的相关系数
结果显示元素Mg和Nd在大米和母质土壤中呈显著正相关(P＜0.05)，Sr在大米和母质土壤中呈极显著正相关(P＜0.05)，元素Te、Gd、Yb和U在大米和母质土壤中呈极显著负相关性(P＜0.01)。
综上，结合大米受产地影响显著的19种元素和大米与母质土壤直接相关元素，最终选定Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Te、Ba、La、Nd、Sm、Gd、Dy、Ho、Er、Yb、Pb和U为相关元素。
对三类大米产地样本中矿物元素数据采集，由于提取了23个有效特征空间的维数，即网络输入层节点数为23；输出层的节点数由分类类型空间的维数决定，研究中将类型分类的编码为100、010和001，分别表示查哈阳大米、五常大米和建三江大米，因此，输出节点数为3；根据Kolmogorov定理和大量实际训练经验，确定第一隐含层神经元节点为36个，第二隐含层神经元节点为28个；所以基于神经网络的大米产地识别分类的神经网络的拓扑结构为23-36-28-3型。
实施例2
随机采集哈阳、五常和建三江地区的274份稻谷样品，利用梯度下降算法求解最优解，定义训练目标迭代精度为0.01，学习速度0.8，惯性系数0.5，最大学习次数4000000，收敛情况如图2所示：前馈神经网络迭代次数为201412次，满足精度误差为0.01。在网络学习过程中，学习算法都是使网络达到设定精度，网络收敛比较稳定，并且符合误差限制要求。
将相关元素信号随机进行归一化处理，如表5所示：
表5.归一化数据
注：编码100、010和001，分别表示查哈阳大米、五常大米和建三江大米。
利用训练好的神经网络的大米产地自动分类方法，将274组待测数据的筛选23中矿物元素归一化数据随机输入溯源模型中，具体产地识别结果如表6所示：
表6.大米产地识别结果
由上述可知，各组数据的平均相对误差值为17.14％，网络的计算输出值和期望输出值的误差较小；三种类型的整体识别准确率为100％。网络相对检测误差分散性较小，具有较好的泛化性和稳定性，是一种有效的大米产地识别方法。
综上所述，本发明提供了一种大米溯源模型，基于人工神经网络可以对决策因素同时处理，有机结合大规模数据的并行处理与串行处理，有效处理各种变化的信息，实现对任意复杂函数的映射，从而适应环境的变化因素。尤其前馈神经网络能够学习和存储大量输入-输出模式映射关系，通过提供大量的不同产地和年份的大米样本，供网络进行学习训练，将提取的样本对中的非线性映射关系存储在权值矩阵中，保证了识别精度和稳定性。建立的前馈神经网络，具有自适应分辨性和良好的容错能力，鲁棒性高且精确达100％。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种大米溯源模型的建立方法，包括以下步骤：
1)收集稻谷及稻谷所在地的土壤母质样品，得检测样本；
2)测定步骤1)得到的所述检测样本中的元素含量，经方差变异贡献率分析和相关性分析，筛选出相关元素；所述元素包括Na、Mg、Al、K、Ca、Sc、V、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Rb、Sr、Y、Mo、Ru、Rh、Pd、Ag、Cd、Sn、Sb、Te、Cs、Ba、La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu、Hf、Ir、Pt、Au、Tl、Pb、Th和U；
3)利用步骤2)筛选出的所述相关元素建立大米溯源的前馈神经网络模型，包括输入层、隐含层和输出层，层与层之间为互连方式，每层分别有nI，nh和no节点数，第一层为输入层，所述输入层的节点数为筛选出的相关元素的个数，设网络输入向量为x，式中T为向量维度或指标个数，只是传递信号到隐含层，令为输入层第i节点的输出：
第二层为第一隐含层，令输入层到第一隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第三层为第二隐含层，令输入层到第二隐含层的连接权为隐含层激励函数为f，阈值为该层输出向量则为隐含层第j节点的输出：
第四层为输出层，令第二隐含层到输出层的连接权为输出层的激励函数p，阈值为输出向量为其中则为输出层的第k节点输出：
2.根据权利要求1所述的建立方法，其特征在于，步骤1)所述样品稻谷的水分含量为14％以下。
3.根据权利要求1所述的建立方法，其特征在于，步骤1)所述粉碎为利用旋风磨粉碎，并过80～120目筛。
4.根据权利要求1所述的建立方法，其特征在于，步骤2)所述测定前还包括微波消解过程，所述微波消解包括：将所述米粉与微波消解液按照1g：(15～35)mL的质量体积比混合，置于微波消解仪中进行微波消解。
5.根据权利要求4所述的建立方法，其特征在于，所述微波消解液为包括浓硝酸和双氧水的混合溶液。
6.根据权利要求4所述的建立方法，其特征在于，所述微波消解的微波程序为8min内从0w增到1600w，温度升到120℃，保持2min；在5min内从温度120℃升到160℃，保持5min；在5min内再从160℃升到180℃，并在此温度下消解15min；冷却。
7.根据权利要求1所述的建立方法，其特征在于，步骤2)所述元素含量利用电感耦合等离子体质谱进行测定。
8.根据权利要求7所述的建立方法，其特征在于，所述电感耦合等离子体质谱的工作参数为射频功率1280w，雾化室温度2℃，冷却水流量1.47L/min，载气流量1.0L/min，补偿气体流量1.0L/min。
9.根据权利要求1所述的建立方法，其特征在于，当鉴别大米是否来源于建三江、五常或查哈阳地区时，步骤2)所述相关元素包括Mg、Ca、Cr、Mn、Zn、As、Rb、Sr、Ag、Cd、Sb、Te、Ba、La、Nd、Sm、Gd、Dy、Ho、Er、Yb、Pb和U。
10.权利要求1～9任一项所述的建立方法建立的大米溯源模型在大米溯源识别中的应用，设共有M样本，第p个样本，对应的期望实际输出为dp，计算网络输出为yp，网络各个输出总误差为
若记wjk数，即连接权值和阈值(θ和)，则学习规则为
式中：η为学习速度；α为惯性系数；t为学习次数；
对得到的所述期望实际输出dp和所述计算网络输出yp进行比对，从而得到相应的稻谷产地。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810758482.5A CN108898509A (zh) | 2018-07-11 | 2018-07-11 | 一种大米溯源模型的建立方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810758482.5A CN108898509A (zh) | 2018-07-11 | 2018-07-11 | 一种大米溯源模型的建立方法及其应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108898509A true CN108898509A (zh) | 2018-11-27 |
Family
ID=64349242
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810758482.5A Pending CN108898509A (zh) | 2018-07-11 | 2018-07-11 | 一种大米溯源模型的建立方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108898509A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111047223A (zh) * | 2019-12-31 | 2020-04-21 | 黑龙江八一农垦大学 | 一种预测水稻中砷元素含量风险评估方法 |
| CN116976915A (zh) * | 2023-05-24 | 2023-10-31 | 南京审计大学 | 地理信息支持的种植类地理标志产品溯源方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102455320A (zh) * | 2010-11-01 | 2012-05-16 | 中国食品发酵工业研究院 | 一种产地葡萄酒识别技术方法 |
| CN104200282A (zh) * | 2014-08-26 | 2014-12-10 | 南京信息工程大学 | 一种设施农业环境因子快速预测方法 |
| CN106770604A (zh) * | 2016-11-11 | 2017-05-31 | 吉林省产品质量监督检验院 | 利用同位素指纹分析溯源大米原产地的方法 |
| KR101887503B1 (ko) * | 2017-06-13 | 2018-08-10 | (주)다온정보 | 인공지능기술을 이용한 시설하우스 제어장치 |
-
2018
- 2018-07-11 CN CN201810758482.5A patent/CN108898509A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102455320A (zh) * | 2010-11-01 | 2012-05-16 | 中国食品发酵工业研究院 | 一种产地葡萄酒识别技术方法 |
| CN104200282A (zh) * | 2014-08-26 | 2014-12-10 | 南京信息工程大学 | 一种设施农业环境因子快速预测方法 |
| CN106770604A (zh) * | 2016-11-11 | 2017-05-31 | 吉林省产品质量监督检验院 | 利用同位素指纹分析溯源大米原产地的方法 |
| KR101887503B1 (ko) * | 2017-06-13 | 2018-08-10 | (주)다온정보 | 인공지능기술을 이용한 시설하우스 제어장치 |
Non-Patent Citations (3)
| Title |
|---|
| 夏立娅: ""大米产地特征因子及溯源方法研究"", 《万方学位论文》 * |
| 张铭钧主编: "《智能控制技术》", 31 December 2006 * |
| 陈小娟著: "《高校本科专业设置预测模型构建》", 31 December 2015 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111047223A (zh) * | 2019-12-31 | 2020-04-21 | 黑龙江八一农垦大学 | 一种预测水稻中砷元素含量风险评估方法 |
| CN116976915A (zh) * | 2023-05-24 | 2023-10-31 | 南京审计大学 | 地理信息支持的种植类地理标志产品溯源方法 |
| CN116976915B (zh) * | 2023-05-24 | 2024-04-05 | 南京审计大学 | 地理信息支持的种植类地理标志产品溯源方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Senay et al. | Combining AVHRR-NDVI and landuse data to describe temporal and spatial dynamics of vegetation | |
| CN107180076B (zh) | 基于高分辨质谱+互联网+地理信息的农药残留可视方法 | |
| Hardisky et al. | Remote sensing of biomass and annual net aerial primary productivity of a salt marsh | |
| Mitchell et al. | Interspecific and environmentally induced variation in foliar dark respiration among eighteen southeastern deciduous tree species | |
| CN107274297A (zh) | 一种土地作物种植适宜性评估方法 | |
| CN106547948A (zh) | 用于引种的园林绿化小乔木生长适应性的综合评价方法 | |
| Klinka et al. | Site diagnosis, tree species selection, and slashburning guidelines for the Vancouver Forest Region | |
| CN104914156A (zh) | 一种基于矿物质分析技术的大米产地鉴定方法与应用 | |
| Li et al. | Collaborative inversion heavy metal stress in rice by using two-dimensional spectral feature space based on HJ-1 A HSI and radarsat-2 SAR remote sensing data | |
| CN108898509A (zh) | 一种大米溯源模型的建立方法及其应用 | |
| KR102684015B1 (ko) | 탄소 흡수 증진을 위한 붉가시나무 다층 식재 모델 | |
| Webb et al. | Studies in the Numerical Analysis of Complex Rain-Forest Communities: V. A Comparison of the Properties of Floristic and Physiognomic-Structural Data | |
| O’Grady et al. | Photosynthesis and respiration decline with light intensity in dominant and suppressed Eucalyptus globulus canopies | |
| CN112666246A (zh) | 一种仿野生岩缝附生栽培铁皮石斛的甄别方法 | |
| CN105184234B (zh) | 一种冬小麦秸秆焚烧污染物排放量的测算方法及装置 | |
| Hutto et al. | Preliminary ecological land classification of the Chauga Ridges Region of South Carolina | |
| Hachmi et al. | Investigation of Flora, Soil, and Biodiversity of Ecosystems in Arid Eastern Moroccan Rangelands | |
| Senay et al. | Capability of AVHRR data in discriminating rangeland cover mixtures | |
| Irwin et al. | Measuring and modeling urban sprawl: Data, scale and spatial dependencies | |
| Fu et al. | Functional relationship between leaf/stem pseudobulb size and photosynthetic pathway in the Orchidaceae | |
| Polley et al. | Impacts of rising CO~ 2 concentration on water use efficiency of woody grassland invaders | |
| CN109241866A (zh) | 一种基于历史作物分布图的农作物自动分类方法及装置 | |
| Bohman | Functional and morphological diversity of trees in different land use types along a rainforest margin in Sulawesi, Indonesia | |
| CN111079672B (zh) | 基于最大熵特征的葡萄分类提取方法 | |
| Darbandi et al. | Flora and weed distribution of pistachio (Pistacia vera L.) orchards in Kerman. |
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 |
Application publication date: 20181127 |
|
| RJ01 | Rejection of invention patent application after publication |