CN113610017B - 一种基于中红外光谱和svm的鹿角帽种类识别方法 - Google Patents
一种基于中红外光谱和svm的鹿角帽种类识别方法 Download PDFInfo
- Publication number
- CN113610017B CN113610017B CN202110918614.8A CN202110918614A CN113610017B CN 113610017 B CN113610017 B CN 113610017B CN 202110918614 A CN202110918614 A CN 202110918614A CN 113610017 B CN113610017 B CN 113610017B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- caps
- deer antler
- svm
- cap
- 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.)
- Active
Links
- 241000282994 Cervidae Species 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000002329 infrared spectrum Methods 0.000 title claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims abstract description 95
- 210000003056 antler Anatomy 0.000 claims abstract description 68
- 238000012706 support-vector machine Methods 0.000 claims abstract description 39
- 241000283007 Cervus nippon Species 0.000 claims abstract description 35
- 238000007637 random forest analysis Methods 0.000 claims abstract description 25
- 238000000513 principal component analysis Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 11
- 241000282985 Cervus Species 0.000 claims abstract description 8
- 238000001035 drying Methods 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 238000013144 data compression Methods 0.000 claims abstract description 4
- 238000013075 data extraction Methods 0.000 claims abstract description 4
- 230000003287 optical effect Effects 0.000 claims abstract description 3
- 238000012360 testing method Methods 0.000 claims description 27
- 238000012549 training Methods 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 19
- IOLCXVTUBQKXJR-UHFFFAOYSA-M potassium bromide Chemical compound [K+].[Br-] IOLCXVTUBQKXJR-UHFFFAOYSA-M 0.000 claims description 12
- 239000000843 powder Substances 0.000 claims description 12
- 230000002068 genetic effect Effects 0.000 claims description 10
- 230000035772 mutation Effects 0.000 claims description 6
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000002790 cross-validation Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 3
- 238000013506 data mapping Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 238000000227 grinding Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000010187 selection method Methods 0.000 claims description 3
- 238000007873 sieving Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000010998 test method Methods 0.000 abstract 1
- 210000004271 bone marrow stromal cell Anatomy 0.000 description 23
- 230000000694 effects Effects 0.000 description 9
- 230000001186 cumulative effect Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 8
- 239000000463 material Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 101000812677 Homo sapiens Nucleotide pyrophosphatase Proteins 0.000 description 2
- 101100353526 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) pca-2 gene Proteins 0.000 description 2
- 102100039306 Nucleotide pyrophosphatase Human genes 0.000 description 2
- 241001416181 Axis axis Species 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 206010061245 Internal injury Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000004907 gland Anatomy 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 239000003053 toxin Substances 0.000 description 1
- 231100000765 toxin Toxicity 0.000 description 1
- 201000008827 tuberculosis Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/16—Classification; Matching by matching signal segments
-
- 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/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于中红外光谱和SVM的鹿角帽种类识别方法,涉及鹿角帽识别技术领域。本发明先对样品进行烘干、粉碎和制片处理,采集漫反射光谱,并对采集的原始光谱进行多重数据处理,利用K‑S检验法对样品进行分类;然后采用归一化与主成分分析降维的方法,对光谱数据进行高维数据压缩和主要特征分量的提取;最后分别将主成分分析后全段SMC光谱和选出有明显差异波段后主成分分析的MSC光谱作为输入变量,建立SVM、ELM和RF梅花鹿鹿角帽和马鹿鹿角帽的识别模型。本发明通过中红外光谱技术以及支持向量机SVM、随机森林算法RF和极限学习机ELM建立模型的方法,实现了梅花鹿鹿角帽和马鹿鹿角帽的高效、准确和无损的鹿角帽品种识别。
Description
技术领域
本发明属于鹿角帽识别技术领域,特别是涉及一种基于中红外光谱和SVM的鹿角帽种类识别方法。
背景技术
鹿角帽是梅花鹿或马鹿采茸后,在雄鹿头上留下的平台状角盘,逐渐骨化,换角时,角炳和角盘分离断裂,脱落下来的角,其有较高的药用价值,可以用来治疗疮痛肿毒,瘀血作痛,虚痨内伤等症,对各种腺体炎症效果明显。随着国内对鹿药材的开发,鹿角帽开始在市场火热起来。虽然梅花鹿鹿角帽和马鹿鹿角帽都是国家规定的鹿药材,但梅花鹿鹿角帽在价格和药用价值上都高于马鹿鹿角帽且二者外形和内部分子结构差不多,导致市场上鹿角帽种类不符成为了常见问题之一,如何找到一种高效,在线,低成本的鹿角帽种类识别方法是重中之重。
目前,国内外学者主要是对鹿茸进行检测识别研究,尽管很多学者对鹿茸的检测识别研究已经有了较好的成果,可对药用鹿角帽检测识别还非常少。红外光谱技术具体有高效,快速,低成本,无损等等点,在农林产品检测及识别方面有着良好的前景。现在对鹿角帽检测识别研究有较强的局限性,测试样本需和标准图谱一模一样才能被认定为正品鹿角帽,否则又要进行其它检测对比且不能对鹿角帽种类进行区分。为此,本发明基于中红外光谱技术的优点与广泛应用结合数学建模寻找到了一种高效、快速且准确的鹿角帽品种识别方法。
发明内容
本发明的目的在于提供一种基于中红外光谱和SVM的鹿角帽种类识别方法,通过中红外光谱技术以及支持向量机SVM、随机森林算法RF和极限学习机ELM建立模型的方法,解决了现有技术中梅花鹿鹿角帽和马鹿鹿角帽无法进行高效、准确和无损的识别,进而无法为解决鹿角帽种类和质量检测问题提供了新的思路和方法的问题。
为解决上述技术问题,本发明是通过以下技术方案实现的:
本发明为一种基于中红外光谱和SVM的鹿角帽种类识别方法,包括以下步骤:
步骤1:对待测的梅花鹿鹿角帽和马鹿鹿角帽分别进行破碎处理,经200目筛筛分后得到梅花鹿鹿角帽和马鹿鹿角帽的粉末,再将鹿角帽粉末和溴化钾放到60℃恒温干燥箱内烘8-12h;
步骤2:精密称取烘干后的梅花鹿鹿角帽粉末1.8mg和溴化钾190mg,混合研磨均匀,研磨后的粉末置于红外压片模具中压制成片,得到梅花鹿鹿角帽压片,同样方式得到马鹿鹿角帽压片,分别将梅花鹿鹿角帽压片和马鹿鹿角帽压片放置于中红外光谱仪上,分别采集对应压片的漫反射光谱;
步骤3:采用The Unscrambler X 10.4软件对采集的原始光谱进行多重数据处理,得到预处理后的光谱,再将预处理后的光谱与对应的原始光谱进行对比;所述多重数据处理包括有多元散射校正MSC;
采用全段光谱数据波段和只从光谱数据中取出差异较大的所有波段进行后面的数据分析对比,从而选出最优的波段;
步骤4:采用K-S检验法,以训练集和测试集样品数为5:2的比例,将总样品数划分为若干训练集和若干测试集,其中训练集和测试集中的梅花鹿鹿角帽和马鹿鹿角帽的样品比例均为1:1;
步骤5:采用归一化与主成分分析降维的方法,对光谱数据进行高维数据压缩和主要特征分量的提取,主要包括如下步骤:
步骤51:采用matlab2014b中的mapminmax函数把MSC光谱数据进行归一化处理,并将数据映射范围设置为0~1;
步骤52:使用Python 3.7软件对归一化后的MSC鹿角帽光谱数据进行主成分分析,并分别绘制全段MSC光谱和选出有明显差异波段MSC光谱中的前两个主成分的散点图;
步骤6:以支持向量机SVM、随机森林RF和极限学习机ELM三个方法建模,分别将主成分分析后全段SMC光谱和选出有明显差异波段后主成分分析的MSC光谱作为输入变量,建立SVM、ELM和RF梅花鹿鹿角帽和马鹿鹿角帽的识别模型。
进一步地,所述步骤2中的漫反射光谱采集软件的波数范围为4000-400cm-1,分辨率为4cm-1,扫描次数为16次,每个样本重复扫描3次,取平均光谱;
光谱采集过程中,室内温度设置为25度,湿度设置为35%。
进一步地,所述步骤3中的多重数据处理还包括有标准正态变量变换SNV、平滑SG、一阶导数和二阶导数。
进一步地,所述步骤5中主成分的个数选择方法为方法一和方法二的结合;
其中,所述方法一为主成分的个数累积贡献率至少要大于等于85%,所述方法二为主成分特征值大于等于1。
进一步地,所述步骤6中的支持向量机SVM的具体方法为:首先训练集采用K-CV交叉验证同时支持向量机SVM需要确定最佳惩罚因子c、核函数参数g及最优核函数;
采用网络搜索法,所述最佳惩罚因子c设置为2-15~215,所述核函数参数g的范围设置为2-15~215,步长为0.1,使用径向基核函数作为最优核函数。
进一步地,所述步骤6中的随机森林RF的具体方法为:采用遗传算法寻找最优棵树及其它影响参数,遗传算法中待优化的变量个数设置为2,个体数目设置为20,最大遗传代数设置为200,变量的二进制位数设置为10,代沟设置为0.95,交叉变异概率设置为0.7,变异概率设置为0.01。
进一步地,所述步骤6中的极限学习机ELM中的隐藏节点数量设置为40-100,比较得出最优的隐藏节点数量。
本发明具有以下有益效果:
1、本发明利用中红外光谱技术以及支持向量机SVM、随机森林算法RF和极限学习机ELM建立模型的方法,实现了梅花鹿鹿角帽和马鹿鹿角帽的高效、准确和无损的鹿角帽品种识别,为解决鹿角帽种类和质量检测问题提供了新的思路和方法。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为MSC鹿角帽的平均光谱图。
图2为归一化后的平均MSC光谱数据图。
表1为前10个主成分特征值和累积贡献率。
图3为全段光谱测试集前2个主成分得分散点图。
图4为差异光谱测试集前2个主成分得分散点图。
图5为全段光谱网格搜索参数寻优适应度曲线。
图6为全段光谱测试集拟合结果。
图7为差异光谱网格搜索参数寻优适应度曲线
图8为差异光谱测试集拟合结果。
图9为全段光谱迭代误差变化曲线。
图10为全段光谱训练集拟合结果。
图11为全段光谱测试集拟合结果。
图12为差异光谱迭代误差变化曲线。
图13为差异光谱训练集拟合结果。
图14为差异光谱测试集拟合结果。
表2为ELM算法预测结果对比。
图15为全段光谱训练集拟合结果。
图16为全段光谱测试集拟合结果。
图17为差异光谱训练集拟合结果。
图18为差异光谱测试集拟合结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1-18和表1-2,本发明为一种基于中红外光谱和SVM的鹿角帽种类识别方法,包括以下步骤:
步骤1:对待测的梅花鹿鹿角帽和马鹿鹿角帽分别进行破碎处理,经200目筛筛分后得到梅花鹿鹿角帽和马鹿鹿角帽的粉末,再将鹿角帽粉末和溴化钾放到60℃恒温干燥箱内烘8-12h;采集梅花鹿鹿角帽与马鹿鹿角帽各42个样品,共84个样品;
步骤2:精密称取烘干后的梅花鹿鹿角帽粉末1.8mg和溴化钾190mg,混合研磨均匀,研磨后的粉末置于红外压片模具中压制成片,得到梅花鹿鹿角帽压片,同样方式得到马鹿鹿角帽压片,分别将梅花鹿鹿角帽压片和马鹿鹿角帽压片放置于中红外光谱仪上,分别采集对应压片的漫反射光谱;
步骤3:光谱信息易受高频随机噪声、基线漂移、样本本身和光散射等影响,需对原始光谱进行预处理,减少这些因素干扰,本发明采用The Unscrambler X 10.4软件对采集的原始光谱进行多重数据处理,得到预处理后的光谱,再将预处理后的光谱与对应的原始光谱进行对比;所述多重数据处理包括有多元散射校正MSC、标准正态变量变换SNV、平滑SG、一阶导数和二阶导数;经过各种预处理后的光谱进行对比,可看出来经过多元散射校正MSC处理的光谱差异性更为明显,如图1所示。由图1可以看出梅花鹿鹿角帽和马鹿鹿角帽在波段740-840,1260-1360,1420-1540,1900-3320,3700-4000中有明显的不同;
在光谱数据中特征峰是判断光谱区别的主要因素,这里采用全段光谱数据波段和只从光谱数据中取出差异较大的所有波段进行后面的数据分析对比,从而选出最优的波段;
步骤4:K-S检验是一种基于累计分布函数,快速检测训练集与测试集划分的方法,采用K-S检验法,以训练集和测试集样品数为5:2的比例,将84份样品划分为60个训练集和24个测试集,其中训练集和测试集中的梅花鹿鹿角帽和马鹿鹿角帽的样品比例均为1:1;训练集中梅花鹿和马鹿鹿角帽各30份,测试集中梅花鹿和马鹿鹿角帽各12份;
步骤5:中红外光谱波段范围为4000-400cm-1,其有波段多、数据量大、冗余性强的特点,采用归一化与主成分分析降维的方法,对光谱数据进行高维数据压缩和主要特征分量的提取,主要包括如下步骤:
步骤51:采用matlab2014b中的mapminmax函数把MSC光谱数据进行归一化处理,并将数据映射范围设置为0~1,如图2所示;
步骤52:使用Python 3.7软件对归一化后的MSC鹿角帽光谱数据进行主成分分析,分别展示全段光谱与选出有明显差异波段光谱的前10个主成分特征值和累积贡献率如表1所示,在全段MSC光谱中PCA1的贡献率最大,为55.62909%,PCA2的贡献率为20.63217%,前3个PC的累积贡献率为84.15974%,直到前8个PC的累积贡献率为97.16765%,之后的各PC贡献率都小于1%且累积贡献率提高速度逐步变小。在选出有明显差异波段MSC光谱中PCA1的贡献率最大,为59.52195%,PCA2的贡献率为29.89027%,前3个PC的累积贡献率为92.68794%,直到前6个PC的累积贡献率为97.92108%,之后的各PC贡献率都小于1%且累积贡献率提高速度逐步变小;
并分别绘制全段MSC光谱和选出有明显差异波段MSC光谱中的前两个主成分的散点图,如图3-4所示;
步骤6:以支持向量机SVM、随机森林RF和极限学习机ELM三个方法建模,分别将主成分分析后全段SMC光谱和选出有明显差异波段后主成分分析的MSC光谱作为输入变量,建立SVM、ELM和RF梅花鹿鹿角帽和马鹿鹿角帽的识别模型。
优选地,所述步骤2中的漫反射光谱采集软件的波数范围为4000-400cm-1,分辨率为4cm-1,扫描次数为16次,每个样本重复扫描3次,取平均光谱;
光谱采集过程中,室内温度设置为25度,湿度设置为35%。
优选地,所述步骤5中主成分的个数选择方法为方法一和方法二的结合;
其中,所述方法一为主成分的个数累积贡献率至少要大于等于85%,所述方法二为主成分特征值大于等于1;在全段光谱数据上,选择前8个主成分构成主成分降维后的光谱数据,又在选出有明显差异波段的光谱数据上,选择前6个主成分构成主成分降维后的光谱数据。
优选地,支持向量机SVM是最好的有监督学习算法之一,该算法基于内核对线性和非线性问题的极限进行建模,可以解决支持分类和回归问题,它对于“过拟合”也非常可行,尤其是在小样本中。所述步骤6中的支持向量机SVM的具体方法为:首先训练集采用K-CV交叉验证同时支持向量机SVM需要确定最佳惩罚因子c、核函数参数g及最优核函数;
采用网络搜索法,所述最佳惩罚因子c设置为2-15~215,所述核函数参数g的范围设置为2-15~215,步长为0.1,使用径向基核函数作为最优核函数;
SVM建模对比,基于全段SMC光谱,和选出有明显差异波段后主成分分析的MSC光谱不同模型训练集,测试集识别效果与确定的c,g如图5-8所示。由图5-8知全段SMC光谱和选出有明显差异波段后主成分分析的MSC光谱所建模型训练集和预测集识别率均为100%,说明SVM对鹿角帽种类识别具有良好效果。
优选地,随机森林RF是一种很灵活实用的方法,具有极好的准确率,能够评估各个特征在分类问题上的重要性,对于缺省值问题也能够获得很好得结果,在BF模型中,设置树的棵树影响着最终结果的好坏。所述步骤6中的随机森林RF的具体方法为:采用遗传算法寻找最优棵树及其它影响参数,遗传算法中待优化的变量个数设置为2,个体数目设置为20,最大遗传代数设置为200,变量的二进制位数设置为10,代沟设置为0.95,交叉变异概率设置为0.7,变异概率设置为0.01;
建立随机森林模型,用遗传算法寻找树的目数,基于全段SMC光谱树的目数为500和选出有明显差异波段后主成分分析的MSC光谱树的目数为600,建模结果如图9-14所示。由图9-14可知,在全段SMC光谱建模中训练集识别率为100%,测试集识别率为95.8333%,有1份梅花鹿鹿角帽识别错误。在选出有明显差异波段后主成分分析的MSC光谱建模中,训练集识别率为100%,测试集识别率为87.5%,有3份梅花鹿鹿角帽识别错误,说明所建立模型出现了过拟合的情况,整体来说全段SMC光谱建模识别率高于选出有明显差异波段后主成分分析的MSC光谱建模识别率。
优选地,极限学习机ELM是一种新型的快速学习算法,学习可以不需要调整隐层节点,也就是说ELM网络隐藏层节点的权重随机生成或者人工定义,学习过程仅需计算输出权重,具有训练参数少、学习速度快、泛化能力强的优点。ELM模型中隐藏节点设置数量与训练集和测试集正确率有直接关系,同时隐藏节点数量的增加又会使算法耗时时间延长,所以根据相关文献对隐藏节点数量的选择为隐层节点数大于训练数据集数量后,正确率增加不明显,且有波动,所述步骤6中的极限学习机ELM中的隐藏节点数量设置为40-100,比较得出最优的隐藏节点数量;
在ELM模型中,选择sigmoidal函数作为激活函数,设置隐藏节点数量为40-100,步长为1进行比较,得表2(为避免数据过于密集,此处仅展示隐藏节点数量从40起以步长为10的建模结果)和图15-18。由表2知当全段SMC光谱建模隐藏节点数量为60,训练集识别率为100%,测试集识别率为95%,有一份梅花鹿鹿角帽识别错误;当选出有明显差异波段后主成分分析的MSC光谱建模隐藏节点数量为50,训练集识别率为100%,测试集识别率为95.833%,有一份梅花鹿鹿角帽识别错误。总体来说ELM算法在2种建模情况下都有稳定的识别效果。
从RF建模效果可以看出,光谱数据特征的提取上采用全段光谱建模的预测集识别率为100%明显高于差异光谱建模的测试集识别率87.5%。从ELM模型中也能看出全段光谱建模预测集识别率为100%略高于差异光谱建模的预测集识别率95.8333%。由此可得全段光谱建模效果总体优于差异光谱建模效果。从建模方法上对比,SVM模型平均识别率为100%,RF模型平均识别率为93.75%,ELM模型平均识别率为97.91665%。可以得出SVM不论是在全段光谱还是差异光谱都有很好的识别效果,非常适用于小样本中红外光谱数据的分析。中红外全段光谱数据结合SVM建模识别效果最佳,为解决鹿角帽种类和质量检测问题提供了新的思路和方法。
在本说明书的描述中,参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。
Claims (1)
1.一种基于中红外光谱和SVM的鹿角帽种类识别方法,其特征在于:包括以下步骤:
步骤1:对待测的梅花鹿鹿角帽和马鹿鹿角帽分别进行破碎处理,经200目筛筛分后得到梅花鹿鹿角帽和马鹿鹿角帽的粉末,再将鹿角帽粉末和溴化钾放到60℃恒温干燥箱内烘8-12h;
步骤2:精密称取烘干后的梅花鹿鹿角帽粉末1.8mg和溴化钾190mg,混合研磨均匀,研磨后的粉末置于红外压片模具中压制成片,得到梅花鹿鹿角帽压片,同样方式得到马鹿鹿角帽压片,分别将梅花鹿鹿角帽压片和马鹿鹿角帽压片放置于中红外光谱仪上,分别采集对应压片的漫反射光谱;
步骤3:采用The UnscramblerX 10.4软件对采集的原始光谱进行多重数据处理,得到预处理后的光谱,再将预处理后的光谱与对应的原始光谱进行对比;所述多重数据处理包括有多元散射校正MSC;
采用全段光谱数据波段和只从光谱数据中取出差异较大的所有波段进行后面的数据分析对比,从而选出最优的波段;
步骤4:采用K-S检验法,以训练集和测试集样品数为5:2的比例,将总样品数划分为若干训练集和若干测试集,其中训练集和测试集中的梅花鹿鹿角帽和马鹿鹿角帽的样品比例均为1:1;
步骤5:采用归一化与主成分分析降维的方法,对光谱数据进行高维数据压缩和主要特征分量的提取,主要包括如下步骤:
步骤51:采用matlab2014b中的mapminmax函数把MSC光谱数据进行归一化处理,并将数据映射范围设置为0~1;
步骤52:使用Python 3.7软件对归一化后的MSC鹿角帽光谱数据进行主成分分析,并分别绘制全段MSC光谱和选出有明显差异波段MSC光谱中的前两个主成分的散点图;
步骤6:以支持向量机SVM、随机森林RF和极限学习机ELM三个方法建模,分别将主成分分析后全段SMC光谱和选出有明显差异波段后主成分分析的MSC光谱作为输入变量,建立SVM、ELM和RF梅花鹿鹿角帽和马鹿鹿角帽的识别模型;
所述步骤2中的漫反射光谱采集软件的波数范围为4000-400cm-1,分辨率为4cm-1,扫描次数为16次,每个样本重复扫描3次,取平均光谱;
光谱采集过程中,室内温度设置为25度,湿度设置为35%;
所述步骤3中的多重数据处理还包括有标准正态变量变换SNV、平滑SG、一阶导数和二阶导数;
所述步骤5中主成分的个数选择方法为方法一和方法二的结合;
其中,所述方法一为主成分的个数累积贡献率至少要大于等于85%,所述方法二为主成分特征值大于等于1;
所述步骤6中的支持向量机SVM的具体方法为:首先训练集采用K-CV交叉验证同时支持向量机SVM需要确定最佳惩罚因子c、核函数参数g及最优核函数;
采用网络搜索法,所述最佳惩罚因子c设置为2-15~215,所述核函数参数g的范围设置为2-15~215,步长为0.1,使用径向基核函数作为最优核函数;
所述步骤6中的随机森林RF的具体方法为:采用遗传算法寻找最优棵树及其它影响参数,遗传算法中待优化的变量个数设置为2,个体数目设置为20,最大遗传代数设置为200,变量的二进制位数设置为10,代沟设置为0.95,交叉变异概率设置为0.7,变异概率设置为0.01;
所述步骤6中的极限学习机ELM中的隐藏节点数量设置为40-100,比较得出最优的隐藏节点数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110918614.8A CN113610017B (zh) | 2021-08-11 | 2021-08-11 | 一种基于中红外光谱和svm的鹿角帽种类识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110918614.8A CN113610017B (zh) | 2021-08-11 | 2021-08-11 | 一种基于中红外光谱和svm的鹿角帽种类识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113610017A CN113610017A (zh) | 2021-11-05 |
CN113610017B true CN113610017B (zh) | 2024-02-02 |
Family
ID=78340265
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110918614.8A Active CN113610017B (zh) | 2021-08-11 | 2021-08-11 | 一种基于中红外光谱和svm的鹿角帽种类识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113610017B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2478823A1 (en) * | 2002-03-11 | 2003-09-18 | Allan L. Schaefer | Method for the evaluation of velvet antler |
CN109034261A (zh) * | 2018-08-10 | 2018-12-18 | 武汉工程大学 | 一种基于支持向量机的近红外光谱数据分析方法 |
CN109374573A (zh) * | 2018-10-12 | 2019-02-22 | 乐山师范学院 | 基于近红外光谱分析的黄瓜表皮农药残留识别方法 |
CN110765962A (zh) * | 2019-10-29 | 2020-02-07 | 刘秀萍 | 一种基于三维点云轮廓分维值的植物识别分类方法 |
CN111024645A (zh) * | 2020-01-08 | 2020-04-17 | 山东金璋隆祥智能科技有限责任公司 | 一种东阿阿胶、福胶、鹿角胶定量研究的方法 |
CN111458306A (zh) * | 2019-01-22 | 2020-07-28 | 吉林农业大学 | 一种红外光谱结合pls-da识别不同花鹿茸的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NZ535178A (en) * | 2002-03-11 | 2007-05-31 | Allan L Schaefer | Method for the evaluation of velvet antler |
-
2021
- 2021-08-11 CN CN202110918614.8A patent/CN113610017B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2478823A1 (en) * | 2002-03-11 | 2003-09-18 | Allan L. Schaefer | Method for the evaluation of velvet antler |
CN109034261A (zh) * | 2018-08-10 | 2018-12-18 | 武汉工程大学 | 一种基于支持向量机的近红外光谱数据分析方法 |
CN109374573A (zh) * | 2018-10-12 | 2019-02-22 | 乐山师范学院 | 基于近红外光谱分析的黄瓜表皮农药残留识别方法 |
CN111458306A (zh) * | 2019-01-22 | 2020-07-28 | 吉林农业大学 | 一种红外光谱结合pls-da识别不同花鹿茸的方法 |
CN110765962A (zh) * | 2019-10-29 | 2020-02-07 | 刘秀萍 | 一种基于三维点云轮廓分维值的植物识别分类方法 |
CN111024645A (zh) * | 2020-01-08 | 2020-04-17 | 山东金璋隆祥智能科技有限责任公司 | 一种东阿阿胶、福胶、鹿角胶定量研究的方法 |
Non-Patent Citations (3)
Title |
---|
中红外光谱的进口木材树种识别方法;冯国红;朱玉杰;李耀翔;;光谱学与光谱分析(07);全文 * |
基于分段主成分分析和高光谱技术的大豆品种识别;刘瑶;谭克竹;陈月华;王志朋;谢红;王立国;;大豆科学(04);全文 * |
太赫兹光谱技术在生物活性肽检测中应用研究;王璞;何明霞;李萌;曲秋红;刘锐;陈永德;;光谱学与光谱分析(09);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113610017A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Feng et al. | Investigation on data fusion of multisource spectral data for rice leaf diseases identification using machine learning methods | |
Pang et al. | Rapid vitality estimation and prediction of corn seeds based on spectra and images using deep learning and hyperspectral imaging techniques | |
Dong et al. | Discrimination of “Hayward” kiwifruits treated with forchlorfenuron at different concentrations using hyperspectral imaging technology | |
CN107271372A (zh) | 一种苹果叶片叶绿素遥感估测方法 | |
Liu et al. | Discrimination of producing areas of Auricularia auricula using visible/near infrared spectroscopy | |
CN112098358A (zh) | 基于四元数卷积神经网络的近红外光谱并行融合定量建模方法 | |
Wu et al. | Deep convolution neural network with weighted loss to detect rice seeds vigor based on hyperspectral imaging under the sample-imbalanced condition | |
CN103743705A (zh) | 一种假高粱及其近似物种的快速检测方法 | |
CN114112983A (zh) | 一种基于Python数据融合的藏药全缘叶绿绒蒿产地判别方法 | |
CN113610017B (zh) | 一种基于中红外光谱和svm的鹿角帽种类识别方法 | |
Liu et al. | Research on the online rapid sensing method of moisture content in famous green tea spreading | |
CN112782148B (zh) | 阿拉比卡和罗伯斯塔两种咖啡豆的快速鉴别方法 | |
CN116735527B (zh) | 一种近红外光谱优化方法、装置、***以及存储介质 | |
Li et al. | HSI combined with CNN model detection of heavy metal Cu stress levels in apple rootstocks | |
CN112881333B (zh) | 一种基于改进免疫遗传算法的近红外光谱波长筛选方法 | |
Li et al. | Data fusion of multiple‐information strategy based on Fourier transform near infrared spectroscopy and Fourier‐transform mid infrared for geographical traceability of Wolfiporia cocos combined with chemometrics | |
CN113049526A (zh) | 一种基于太赫兹衰减全反射的玉米种子水分含量测定方法 | |
CN114062306A (zh) | 一种近红外光谱数据分段预处理方法 | |
CN111595802A (zh) | 一种基于nir光谱的忧遁草种源地分类模型的构建方法及应用 | |
CN112763448A (zh) | 一种基于atr-ftir技术的米糠中多糖含量的快速检测方法 | |
Chen et al. | Classification of wheat grain varieties using terahertz spectroscopy and convolutional neural network | |
CN105651727A (zh) | 基于jade和elm的近红外光谱分析鉴别苹果货架期的方法 | |
Zhang et al. | Apple identity recognition based on SVM model parameter optimization and near infrared hyperspectral. | |
He et al. | A rapid method for identification of Lanxangia tsaoko origin and fruit shape: FT‐NIR combined with chemometrics and image recognition | |
CN108982408A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |