阅读说明：本技术 一种多源光谱融合的滇龙胆种类识别方法及系统 (Multi-source spectrum fused gentiana rigescens species identification method and system ) 是由 孙明华 孔汶汶 孙永祺 于 2020-06-04 设计创作，主要内容包括：本发明涉及一种多源光谱融合的滇龙胆种类识别方法及系统,包括：获取待识别滇龙胆切片样本；利用LIBS光谱仪采集待识别滇龙胆切片样本不同位置的光谱数据,得待识别滇龙胆LIBS光谱数据；根据波长与元素之间的对应关系,对待识别滇龙胆LIBS光谱数据进行筛选,确定筛选待识别滇龙胆LIBS光谱数据；利用NIR光谱仪采集待识别滇龙胆切片样本的像素点光谱数据,得待识别滇龙胆近红外光谱数据；采用主成分分析方法对上述光谱数据进行降维和融合,得待识别融合光谱数据；根据待识别融合光谱数据,采用偏最小二乘判别分析方程确定待识别滇龙胆切片样本的种类。通过本发明的上述方法及系统,实现野生滇龙胆快速的实时识别,提高识别的精确性。(The invention relates to a method and a system for identifying the type of gentiana rigescens with multi-source spectrum fusion, which comprises the following steps: obtaining a gentiana rigescens section sample to be identified; collecting spectrum data of different positions of the gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectrum data of the gentiana rigescens to be identified; screening LIBS spectral data of the gentiana rigescens to be identified according to the corresponding relation between the wavelength and the elements, and determining to screen the LIBS spectral data of the gentiana rigescens to be identified; collecting pixel point spectrum data of the gentiana rigescens sliced sample to be identified by using an NIR spectrometer to obtain gentiana rigescens near infrared spectrum data to be identified; performing dimensionality reduction and fusion on the spectral data by adopting a principal component analysis method to obtain fused spectral data to be identified; and determining the type of the gentiana rigescens section sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectral data to be identified. By the method and the system, the wild gentiana rigescens can be rapidly identified in real time, and the identification accuracy is improved.)

1. A multisource spectrum fused gentiana rigescens species identification method is characterized by comprising the following steps:

obtaining a gentiana rigescens section sample to be identified;

collecting spectrum data of different positions of the gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectrum data of the gentiana rigescens to be identified; the LIBS spectrum data of the gentiana rigescens to be identified comprise a plurality of LIBS spectra of the gentiana rigescens to be identified;

according to the corresponding relation between the wavelength and the elements, the wavelength of each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified is screened, and the Gentiana rigescens LIBS spectrum data to be identified is determined to be screened;

collecting pixel point spectrum data of the gentiana rigescens slice sample to be identified by using an NIR spectrometer to obtain gentiana rigescens near infrared spectrum data to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified;

performing dimensionality reduction and fusion on the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain fused spectral data to be identified;

determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

2. The multisource spectrum fused gentiana rigescens species identification method according to claim 1, wherein the screening of the LIBS spectral data of gentiana rigescens to be identified and the near infrared spectral data of gentiana rigescens to be identified are subjected to dimensionality reduction and fusion by a principal component analysis method to obtain fused spectral data to be identified, and the method specifically comprises the following steps:

determining a Gentiana rigescens LIBS spectral matrix to be identified according to the screened Gentiana rigescens LIBS spectral data to be identified;

obtaining a first projection vector set by adopting a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value;

selecting 8 first projection vectors with the maximum contribution value in the first projection vector set as LIBS spectral projection vectors of the gentiana rigescens to be identified;

determining a Gentiana rigescens LIBS scoring matrix to be identified according to the Gentiana rigescens LIBS spectral matrix to be identified and the Gentiana rigescens LIBS spectral projection vector to be identified;

determining a gentiana rigescens near infrared spectrum matrix to be identified according to the gentiana rigescens near infrared spectrum data to be identified;

obtaining a second projection vector set by adopting a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value;

selecting 3 second projection vectors with the maximum contribution value in the second projection vector set as the near infrared spectrum projection vectors of the gentiana rigescens to be identified;

determining a near infrared score matrix of the gentiana rigescens to be identified according to the near infrared spectrum matrix of the gentiana rigescens to be identified and the near infrared spectrum projection vector of the gentiana rigescens to be identified;

and fusing the LIBS scoring matrix of the gentiana rigescens to be identified and the near-infrared scoring matrix of the gentiana rigescens to be identified to obtain fused spectral data to be identified.

3. The method for identifying the type of the gentiana rigescens slice sample through multi-source spectrum fusion according to claim 1, wherein the determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectrum data to be identified specifically comprises the following steps:

according to the formulaDetermining a discrimination coefficient; wherein Y is a discrimination coefficient, p₁,p₂,p₃,p₄,p₅,p₆,p₇,p₈,p₉,p₁₀And p₁₁All are fused spectral data to be identified;

when the absolute value Y-1 is less than 0.5, determining the gentiana rigescens slice sample to be identified as the domestic gentiana rigescens;

and when the absolute value Y-2 is less than 0.5, determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens.

4. The multisource spectrum fused gentiana rigescens species identification method according to claim 1, wherein the screening of the LIBS spectral data of gentiana rigescens to be identified and the near infrared spectral data of gentiana rigescens to be identified are subjected to dimensionality reduction and fusion by a principal component analysis method to obtain fused spectral data to be identified, and the method further comprises the following steps:

and respectively processing the LIBS spectral data of the gentiana rigescens to be screened and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be normalized and the near infrared spectral data of the gentiana rigescens to be normalized.

5. The multisource spectrum fused gentiana rigescens species identification method according to claim 1, wherein the NIR spectrometer is used for collecting pixel point spectrum data of the gentiana rigescens sliced sample to be identified to obtain gentiana rigescens near infrared spectrum data to be identified, and then the method further comprises the following steps:

and denoising each gentiana rigescens near infrared spectrum to be identified in the gentiana rigescens near infrared spectrum data to be identified by adopting a wavelet denoising method and a moving window method to obtain the gentiana rigescens near infrared spectrum data to be denoised.

6. The multisource spectrum fused gentiana rigescens species identification method according to claim 1, wherein the collection of spectral data of different positions of the gentiana rigescens sliced sample to be identified by using a LIBS spectrometer obtains the LIBS spectral data of the gentiana rigescens to be identified, and then the method further comprises the following steps:

and denoising each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified by adopting a wavelet denoising method to obtain the Gentiana rigescens LIBS spectrum data to be denoised.

7. The utility model provides a yunnan gentian kind identification system of multisource spectrum fusion which characterized in that includes:

the cutting sample acquisition module is used for acquiring a gentiana rigescens cutting sample to be identified;

the LIBS spectrum data acquisition module is used for acquiring spectrum data of different positions of the Gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectrum data of the Gentiana rigescens to be identified; the LIBS spectrum data of the gentiana rigescens to be identified comprise a plurality of LIBS spectra of the gentiana rigescens to be identified;

the LIBS spectrum data determination module is used for screening the wavelengths of each LIBS spectrum to be identified in the LIBS spectrum data of the gentiana rigescens to be identified according to the corresponding relation between the wavelengths and the elements, and determining to screen the LIBS spectrum data of the gentiana rigescens to be identified;

the method comprises the steps that a gentiana rigescens near infrared spectrum data acquisition module to be identified collects pixel point spectrum data of a gentiana rigescens slice sample to be identified by using an NIR spectrometer to obtain gentiana rigescens near infrared spectrum data to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified;

the system comprises a module for obtaining fused spectral data to be identified, and a module for obtaining fused spectral data to be identified, wherein the module is used for reducing the dimension and fusing the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain the fused spectral data to be identified;

the identification module is used for determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least square discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

8. The multi-source spectrum fused gentiana rigescens species identification system according to claim 7, wherein the fused spectrum data obtaining module to be identified specifically comprises:

the LIBS spectral matrix determining unit is used for determining the LIBS spectral matrix of the gentiana rigescens to be identified according to the LIBS spectral data of the gentiana rigescens to be identified;

a first projection vector set obtaining unit, configured to obtain a first projection vector set by using a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value;

the LIBS spectral projection vector selecting unit is used for selecting 8 first projection vectors with the largest contribution value in the first projection vector set as LIBS spectral projection vectors of the gentiana rigescens to be identified;

the LIBS score matrix determining unit is used for determining the LIBS score matrix of the gentiana rigescens to be identified according to the LIBS spectrum matrix of the gentiana rigescens to be identified and the LIBS spectrum projection vector of the gentiana rigescens to be identified;

the near infrared spectrum matrix determining unit is used for determining a near infrared spectrum matrix of the gentiana rigescens to be identified according to the near infrared spectrum data of the gentiana rigescens to be identified;

a second projection vector set obtaining unit, configured to obtain a second projection vector set by using a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value;

the near infrared spectrum projection vector selecting unit is used for selecting 3 second projection vectors with the largest contribution value in the second projection vector set as near infrared spectrum projection vectors of the gentiana rigescens to be identified;

the near-infrared score matrix determining unit is used for determining a near-infrared score matrix of the gentiana rigescens to be identified according to the near-infrared spectrum matrix of the gentiana rigescens to be identified and the near-infrared spectrum projection vector of the gentiana rigescens to be identified;

and the spectrum data fusion unit is used for fusing the LIBS score matrix of the gentiana rigescens to be identified and the near-infrared score matrix of the gentiana rigescens to be identified to obtain fusion spectrum data to be identified.

9. The multi-source spectrum fused gentiana rigescens species identification system according to claim 7, wherein the identification module specifically comprises:

a discrimination coefficient determination unit forAccording to the formula

the family gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the family gentiana rigescens when the absolute value Y-1 is less than 0.5;

and the wild gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens when the absolute value Y-2 is less than 0.5.

10. The multi-source spectrally fused gentiana rigescens species identification system of claim 7, wherein said multi-source spectrally fused gentiana rigescens species identification system further comprises:

and the area normalization processing module is used for respectively processing the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified.

Technical Field

The invention relates to the technical field of gentiana rigescens identification, in particular to a gentiana rigescens type identification method and system based on multi-source spectrum fusion.

Background

Root and rhizome of Gentiana rigescens (Gentiana rigescens Franch) which is a perennial herb is used as Chinese Gentiana which is a good medicine for purging liver and gallbladder fire. Modern phytochemistry and pharmacological research shows that gentiana rigescens has pharmacological actions such as liver protection, inflammation diminishing, pain relieving, virus resisting and the like, and iridoid glycoside components such as gentiopicroside, swertiamarin, loganin, swertiamarin and the like are the main material basis for exerting the pharmacological actions.

The gentian plants mostly grow on mountains or plateaus, and the conditions of the growing environment, the growing climate and the like are extreme and variable. At present, the main distribution areas of gentiana rigescens are Yunnan, Sichuan, Guizhou, Guangxi and the like. Due to the difference of factors such as temperature, precipitation, soil organic content and the like, the growth and metabolism of the gentiana rigescens are influenced.

Due to the difference of the growing environment and the growing age of the gentiana rigescens, the family and the wild have obvious difference in some aspects: 1) in appearance, wild lean, weak and compact, and domestic rich and loose bubbles; 2) the growth period of wild plants is generally long, the components are balanced and rich, the curative effect is good, and the problems of non-standard planting technology, hormone abuse and the like exist in the medicinal materials of families. Therefore, the wild identification of the Chinese medicinal material gentiana rigescens is particularly important.

The traditional identification method of wild gentiana rigescens is mainly realized by a chemical method, namely, the source of gentiana rigescens is determined by measuring the composition proportion and content value of various nutrient components in gentiana rigescens.

Disclosure of Invention

The invention aims to provide a method and a system for identifying the type of gentiana rigescens through multi-source spectrum fusion, which solve the problems of high labor cost and complex operation in the prior art, realize quick real-time identification of wild gentiana rigescens and improve the identification accuracy.

In order to achieve the purpose, the invention provides the following scheme:

a multisource spectrum fused gentiana rigescens species identification method comprises the following steps:

obtaining a gentiana rigescens section sample to be identified;

collecting spectrum data of different positions of the gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectrum data of the gentiana rigescens to be identified; the LIBS spectrum data of the gentiana rigescens to be identified comprise a plurality of LIBS spectra of the gentiana rigescens to be identified;

according to the corresponding relation between the wavelength and the elements, the wavelength of each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified is screened, and the Gentiana rigescens LIBS spectrum data to be identified is determined to be screened;

collecting pixel point spectrum data of the gentiana rigescens slice sample to be identified by using an NIR spectrometer to obtain gentiana rigescens near infrared spectrum data to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified;

performing dimensionality reduction and fusion on the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain fused spectral data to be identified;

determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

Optionally, the screening of the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified are subjected to dimensionality reduction and fusion by a principal component analysis method to obtain fused spectral data to be identified, and the method specifically includes the following steps:

determining a Gentiana rigescens LIBS spectral matrix to be identified according to the screened Gentiana rigescens LIBS spectral data to be identified;

obtaining a first projection vector set by adopting a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value;

selecting 8 first projection vectors with the maximum contribution value in the first projection vector set as LIBS spectral projection vectors of the gentiana rigescens to be identified;

determining a Gentiana rigescens LIBS scoring matrix to be identified according to the Gentiana rigescens LIBS spectral matrix to be identified and the Gentiana rigescens LIBS spectral projection vector to be identified;

determining a gentiana rigescens near infrared spectrum matrix to be identified according to the gentiana rigescens near infrared spectrum data to be identified;

obtaining a second projection vector set by adopting a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value;

selecting 3 second projection vectors with the maximum contribution value in the second projection vector set as the near infrared spectrum projection vectors of the gentiana rigescens to be identified;

determining a near infrared score matrix of the gentiana rigescens to be identified according to the near infrared spectrum matrix of the gentiana rigescens to be identified and the near infrared spectrum projection vector of the gentiana rigescens to be identified;

and fusing the LIBS scoring matrix of the gentiana rigescens to be identified and the near-infrared scoring matrix of the gentiana rigescens to be identified to obtain fused spectral data to be identified.

Optionally, the determining the type of the gentiana rigescens section sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectral data to be identified specifically includes:

according to the formulaDetermining a discrimination coefficient; wherein Y is a discrimination coefficient, p₁,p₂,p₃,p₄,p₅,p₆,p₇,p₈,p₉,p₁₀And p₁₁All are fused spectral data to be identified;

when the absolute value Y-1 is less than 0.5, determining the gentiana rigescens slice sample to be identified as the domestic gentiana rigescens;

and when the absolute value Y-2 is less than 0.5, determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens.

Optionally, the screening of the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified are subjected to dimensionality reduction and fusion by using a principal component analysis method to obtain fused spectral data to be identified, and the method further includes the following steps:

and respectively processing the LIBS spectral data of the gentiana rigescens to be screened and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be normalized and the near infrared spectral data of the gentiana rigescens to be normalized.

Optionally, the NIR spectrometer is used to collect the pixel point spectral data of the gentiana rigescens slice sample to be identified, so as to obtain the gentiana rigescens near infrared spectral data to be identified, and then the NIR spectrometer further comprises:

and denoising each gentiana rigescens near infrared spectrum to be identified in the gentiana rigescens near infrared spectrum data to be identified by adopting a wavelet denoising method and a moving window method to obtain the gentiana rigescens near infrared spectrum data to be denoised.

Optionally, the acquiring, by using a LIBS spectrometer, spectral data of different positions of the gentiana rigescens slice sample to be identified to obtain the LIBS spectral data of the gentiana rigescens to be identified, and then further including:

and denoising each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified by adopting a wavelet denoising method to obtain the Gentiana rigescens LIBS spectrum data to be denoised.

A kind of gentiana rigescens kind identification system of multisource spectrum fusion, includes:

the cutting sample acquisition module is used for acquiring a gentiana rigescens cutting sample to be identified;

the LIBS spectrum data acquisition module is used for acquiring spectrum data of different positions of the Gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectrum data of the Gentiana rigescens to be identified; the LIBS spectrum data of the gentiana rigescens to be identified comprise a plurality of LIBS spectra of the gentiana rigescens to be identified;

the LIBS spectrum data determination module is used for screening the wavelengths of each LIBS spectrum to be identified in the LIBS spectrum data of the gentiana rigescens to be identified according to the corresponding relation between the wavelengths and the elements, and determining to screen the LIBS spectrum data of the gentiana rigescens to be identified;

the method comprises the steps that a gentiana rigescens near infrared spectrum data acquisition module to be identified collects pixel point spectrum data of a gentiana rigescens slice sample to be identified by using an NIR spectrometer to obtain gentiana rigescens near infrared spectrum data to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified;

the system comprises a module for obtaining fused spectral data to be identified, and a module for obtaining fused spectral data to be identified, wherein the module is used for reducing the dimension and fusing the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain the fused spectral data to be identified;

the identification module is used for determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least square discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

Optionally, the fused spectral data to be identified obtaining module specifically includes:

the LIBS spectral matrix determining unit is used for determining the LIBS spectral matrix of the gentiana rigescens to be identified according to the LIBS spectral data of the gentiana rigescens to be identified;

a first projection vector set obtaining unit, configured to obtain a first projection vector set by using a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value;

the LIBS spectral projection vector selecting unit is used for selecting 8 first projection vectors with the largest contribution value in the first projection vector set as LIBS spectral projection vectors of the gentiana rigescens to be identified;

the LIBS score matrix determining unit is used for determining the LIBS score matrix of the gentiana rigescens to be identified according to the LIBS spectrum matrix of the gentiana rigescens to be identified and the LIBS spectrum projection vector of the gentiana rigescens to be identified;

the near infrared spectrum matrix determining unit is used for determining a near infrared spectrum matrix of the gentiana rigescens to be identified according to the near infrared spectrum data of the gentiana rigescens to be identified;

a second projection vector set obtaining unit, configured to obtain a second projection vector set by using a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value;

the near infrared spectrum projection vector selecting unit is used for selecting 3 second projection vectors with the largest contribution value in the second projection vector set as near infrared spectrum projection vectors of the gentiana rigescens to be identified;

the near-infrared score matrix determining unit is used for determining a near-infrared score matrix of the gentiana rigescens to be identified according to the near-infrared spectrum matrix of the gentiana rigescens to be identified and the near-infrared spectrum projection vector of the gentiana rigescens to be identified;

and the spectrum data fusion unit is used for fusing the LIBS score matrix of the gentiana rigescens to be identified and the near-infrared score matrix of the gentiana rigescens to be identified to obtain fusion spectrum data to be identified.

Optionally, the identification module specifically includes:

a discrimination coefficient determining unit for determining the discrimination coefficient based on the formulaDetermining a discrimination coefficient; wherein Y is a discrimination coefficient, p₁,p₂,p₃,p₄,p₅,p₆,p₇,p₈,p₉,p₁₀And p₁₁All are fused spectral data to be identified;

the family gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the family gentiana rigescens when the absolute value Y-1 is less than 0.5;

and the wild gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens when the absolute value Y-2 is less than 0.5.

Optionally, the multi-source spectrum fused gentiana rigescens species identification system further comprises:

and the area normalization processing module is used for respectively processing the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified.

According to the specific embodiment provided by the invention, the invention discloses the following technical effects:

the invention provides a method and a system for identifying the type of gentiana rigescens with multi-source spectrum fusion, which are characterized in that a LIBS spectrometer and an NIR spectrometer are used for respectively carrying out spectrum acquisition on sliced samples, the acquired spectrum data are processed, then the processed spectrum data are subjected to dimensionality reduction and fusion by adopting a principal component analysis method, and finally the fusion data are input into a partial least square discriminant analysis equation to determine the type of the gentiana rigescens sliced samples to be identified. The near infrared band represents molecular composition and structure information of a substance, and the LIBS spectrum reflects element composition and content information of the substance, so that the invention utilizes the PCA method to perform dimensionality reduction and fusion on spectral data of the two technologies, fully combines the advantages of the two technologies, reduces spectral variables through the PCA method, shortens calculation time, and greatly improves identification accuracy.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.

Fig. 1 is a flowchart of a method for identifying a kind of gentiana rigescens with multi-source spectrum fusion according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of the relationship between the spectral wavelength and the element provided by the embodiment of the present invention;

FIG. 3 is a three-dimensional scattergram of the LIBS score matrix according to the embodiment of the present invention;

FIG. 4 is a three-dimensional scattergram of a near-infrared scoring matrix according to an embodiment of the present invention;

FIG. 5 is a schematic diagram of the degree of discrimination provided by the embodiment of the present invention;

fig. 6 is a schematic structural diagram of a gentiana rigescens species identification system with multi-source spectrum fusion provided by an embodiment of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The invention aims to provide a method and a system for identifying the type of gentiana rigescens through multi-source spectrum fusion, which solve the problems of high labor cost and complex operation in the prior art, realize quick real-time identification of wild gentiana rigescens and improve the identification accuracy.

In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.

Fig. 1 is a flowchart of a method for identifying a kind of a gentiana rigescens based on multi-source spectrum fusion according to an embodiment of the present invention, and as shown in fig. 1, the method for identifying a kind of gentiana rigescens based on multi-source spectrum fusion according to the present invention includes:

s101, obtaining a gentiana rigescens section sample to be identified. Specifically, the gentiana rigescens plant to be identified is collected, and the gentiana rigescens plant to be identified can be one plant or multiple plants, and one plant is taken as an example in the embodiment of the invention. Cleaning the roots of the plants, drying and grinding, then taking 0.15g of each sample of the plants, pressing into a piece, and preparing the gentiana rigescens slice sample. When the slice is manufactured, the surface of the slice is relatively flat, and the thickness of the slice is constant. The surface of the sliced piece is relatively flat, so that the experimental error caused by uneven powder surface can be reduced when the near infrared spectrum data of the gentiana rigescens is acquired; the thickness of the tablet is fixed, and when LIBS spectral data of the gentiana rigescens are collected, errors caused by the fact that the laser focus point is not consistent with the surface height of a sample can be reduced.

S102, collecting spectral data of different positions of the gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the LIBS spectral data of the gentiana rigescens to be identified; the LIBS spectrum data of the Gentiana rigescens to be identified comprise a plurality of LIBS spectra of the Gentiana rigescens to be identified.

Specifically, for LIBS spectral data, a LIBS spectrometer is used to perform spectral acquisition on different positions of each slice sample, and in the embodiment of the present invention, 16 different positions of one slice sample are selected to perform spectral acquisition. The spectra at each point are averaged over 5 laser hits, i.e. 16 spectra were collected for one sample section. Each spectrum has 22015 wavelengths and the spectral range is 229.1-877.48nm, namely each spectrum has 22015 spectral variables. The LIBS spectrometer parameters for spectrum acquisition were: the laser energy is 60mJ, the laser wavelength is 532nm, the focal depth is 2mm, the delay time is 2 mus, and the integration time is 20 mus.

After S102, further comprising: and denoising each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified by adopting a wavelet denoising method to obtain the Gentiana rigescens LIBS spectrum data to be denoised. Specifically, for the LIBS spectrum, the wavelet parameters with the highest signal-to-noise ratio are searched in different wavelet transformation parameters with the wavelet basis functions db 3-db 10 and the number of layers of 3-10, that is, the optimal wavelet parameters are as follows: the basis function is db5, and the number of layers is 3. On the basis, each LIBS spectrum of the gentiana rigescens to be identified is subjected to wavelet desizing.

S103, according to the corresponding relation between the wavelength and the elements, the wavelength of each Gentiana rigescens LIBS spectrum to be identified in the Gentiana rigescens LIBS spectrum data to be identified is screened, and the Gentiana rigescens LIBS spectrum data to be identified is determined to be screened.

Referring to the spectral wavelengths corresponding to elements provided by the national institute of technology (NIST), there are 42 spectral wavelengths corresponding to the nutritive elements and trace elements of gentiana yunnanensis, such that each gentiana yunnanensis LIBS spectrum is reduced from 22015 spectral variables to the spectral variables corresponding to the wavelengths of 42 relevant elements, i.e., to 42 spectral variables. The correspondence between elements and wavelengths is shown in table 1.

TABLE 1 element to wavelength correspondence

Specifically, 22015 spectrum variables of each gentiana rigescens LIBS collected by the LIBS spectrometer represent element types and relative content information, but the element types contained in the gentiana rigescens herb are limited, so that a large amount of background information is contained in the 22015 spectrum variables, specifically, no obvious excitation peak signal exists, the relative intensity of the spectrum is extremely low, and the redundant information has a strong interference effect on the optimization of the classification model, so that the spectrum wavelength corresponding to the specific element is extracted by referring to the information provided by the national labeling and technical research institute of america, and the LIBS spectrum variables are reduced from 22015 spectrum variables to 42 spectrum variables. And screening out the spectrums corresponding to the 42 spectrum variables from the Gentiana rigescens LIBS spectrum variables to form new Gentiana rigescens LIBS spectrum variables, so that each new Gentiana rigescens LIBS spectrum has 42 wavelengths, namely 42 spectrum variables.

S104, collecting pixel point spectral data of the gentiana rigescens slice sample to be identified by using an NIR spectrometer to obtain near infrared spectral data of the gentiana rigescens to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified.

Specifically, for near infrared spectrum data, an NIR spectrometer is used to collect pixel point spectra of a sliced sample in a line scanning manner. The NIR spectroscopy system mainly comprises: the imaging spectrometer has a spectrum signal receiving range of 874-1734nm and 256 wave bands in the spectrum range, namely each collected spectrum has 256 spectrum variables and tens of thousands of pixel point spectra are collected. The system is provided with two 150w halogen tungsten lamp line light sources, an IRCP0076 type electric control displacement platform which can drive a conveyor belt carrying samples, a computer for controlling the operation of the system and hyperspectral imaging system acquisition software provided by Taiwan Wuling optical corporation. The spectral resolution of the system is 5nm, and the image resolution is 320x256 pixels. Therefore, a total of ten thousand spectra were collected in the spectral range 874-1734nm, each spectrum having 256 spectral variations.

After S104, further comprising: and denoising each gentiana rigescens near infrared spectrum to be identified in the gentiana rigescens near infrared spectrum data to be identified by adopting a wavelet denoising method and a moving window method to obtain the gentiana rigescens near infrared spectrum data to be denoised.

Specifically, for each gentiana rigescens near infrared spectrum, because the noise near the first segment and the tail end of the wavelength in the spectral range is large and the stability is low, the spectrum near the first segment and the tail end is removed, 256 bands are included in the 874-1734nm spectral range, and the spectral data of 192 bands among 31-222 bands are selected, namely the gentiana rigescens near infrared spectrum is changed into 192 spectral variables.

Because the number of pixel point spectra collected by the NIR spectrometer is far more than that collected by the LIBS spectrometer, the number of gentiana rigescens near infrared spectra to be identified needs to be reduced to 16 in order to ensure that the number of spectra collected by the NIR spectrometer corresponds to that collected by the LIBS spectrometer. The spectral data collected by the instrument is in the form of,

x_jrepresents the jth spectral vector, j is 1,2,. n, and is divided into 16 intervals from top to bottom, and the number of spectral strips in each interval is approximately

(if N can not be evenly divided, a plurality of spectra are put in the last interval, and 16 intervals are ensured in total), and the spectra in each interval are averaged to obtain 16 spectra. Then, performing wavelet transformation on the 16 spectra, and searching a wavelet parameter with the highest signal-to-noise ratio in different wavelet transformation parameters with a wavelet basis function of db 3-db 10 and 3-10 layers, wherein the optimal wavelet parameter is as follows: the basis function is db3, and the number of layers is 3. Based on the wavelet parameters, the wavelet denoising is carried out on the 16 near infrared spectrums, and in order to eliminate the burr noise on the spectrum image, the moving average processing with the moving window of 7 is carried out on each spectrum.

Before S105, the method further includes: and respectively processing the LIBS spectral data of the gentiana rigescens to be screened and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be normalized and the near infrared spectral data of the gentiana rigescens to be normalized.

And S105, performing dimensionality reduction and fusion on the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain fused spectral data to be identified.

S105 specifically comprises the following steps:

1. the acquisition process of the LIBS score matrix of the gentiana rigescens to be identified comprises the following steps: 1-1) determining a Gentiana rigescens LIBS spectral matrix to be identified according to the screened Gentiana rigescens LIBS spectral data to be identified; 1-2) obtaining a first projection vector set by adopting a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value; 1-3) selecting 8 first projection vectors with the maximum contribution value in the first projection vector set as LIBS spectral projection vectors of the gentiana rigescens to be identified; 1-4) determining a Gentiana rigescens LIBS score matrix to be identified according to the Gentiana rigescens LIBS spectral matrix to be identified and the Gentiana rigescens LIBS spectral projection vector to be identified.

2. The process for acquiring the near-infrared score matrix of the gentiana rigescens to be identified comprises the following steps: 2-1) determining a near infrared spectrum matrix of the gentiana rigescens to be identified according to the near infrared spectrum data of the gentiana rigescens to be identified; 2-2) obtaining a second projection vector set by adopting a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value; 2-3) selecting 3 second projection vectors with the maximum contribution rate value in the second projection vector set as the near infrared spectrum projection vectors of the gentiana rigescens to be identified; and 2-4) determining a near infrared score matrix of the gentiana rigescens to be identified according to the near infrared spectrum matrix of the gentiana rigescens to be identified and the near infrared spectrum projection vector of the gentiana rigescens to be identified.

3. And fusing the LIBS scoring matrix of the gentiana rigescens to be identified and the near-infrared scoring matrix of the gentiana rigescens to be identified to obtain fused spectral data to be identified.

Specifically, the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified are screened and subjected to dimensionality reduction by a Principal Component Analysis (PCA) method. Firstly, preprocessing each spectrum in LIBS spectrum data of the gentiana rigescens and near infrared spectrum data of the gentiana rigescens by adopting an area normalization method, and eliminating errors caused by experiment parameters and uneven surfaces of press-cut pieces. And then, carrying out dimensionality reduction and fusion processing on the normalized data by adopting a PCA (principal component analysis) method.

The area normalization process for each spectrum is as follows:

according to the formulaArea normalization was performed for each spectrum, where X_NEWiFor area normalized spectra, x_iRepresents the ith spectral variable in a spectrum, i is 1,2,.. n, if the near infrared spectral data of the gentiana rigescens, n is 192, and if the LIBS spectral data of the gentiana rigescens, n is 42.

The PCA dimensionality reduction and fusion process is as follows:

recording the wavelength of each spectrum in the LIBS spectrum data after area normalization as lambda₁,λ₂,λ₃...λ_p(p-42), one spectrum of LIBS spectral data has 42 wavelengths corresponding to a line intensity I λ₁,Iλ₂,……Iλ_pA spectral vector of LIBS spectral data can be represented as X_j，

j 1, 2.. 16, j represents a spectrum number, 16 spectra are in total, and a composition spectrum matrix can be represented as

Recording projection vector V ═ a₁,a₂……a_p]^TThen score the variable X_LIBS× V variance σ²Can be expressed as:

the objective of the PCA analysis is to make σ²The larger the variance is, the better the projection vectors V with different variances are obtained through a Lagrange solution, the larger the variance is, the larger the contribution rate is, and the 8 projection vectors V with the maximum total contribution rate value are found₁,V₂,V₃,V₄,V₅,V₆,V₇,V₈Score matrix S_LIBSCan be expressed as: s_LIBS＝X_LIBS[V₁…… V₈]，S_LIBSIs 16 rows and 8 columns.

Similarly, a near-infrared scoring matrix S is obtained_NIR＝X_NIR[V₁,V₂,V₃]，S_NIRIs 16 rows and 3 columns, and the total contribution value is the largest of the 3 projection vectors when the near infrared spectrum data is subjected to PCA analysis.

Two kinds of data S_LIBS,S_NIRFusing 11 lines of score variables to obtain fused spectrum data to be identified, wherein the fused spectrum data to be identified can comprehensively represent LIBS and near-infrared spectrum information of the gentiana rigescens.

PCA is an unsupervised method, linear combination is carried out on spectral variables to form new linearly independent score variables capable of effectively representing gentiana rigescens information, the linear combination mode is vector projection, the goal is to enable the variance of the score vectors after the spectral variables are projected to be maximum, according to the variance, the PCA can generate projection vectors with different contribution rates (the larger the variance is, the larger the contribution rate is), and the score variables corresponding to the projection vectors with the total contribution rate value larger than 99% are taken to represent the original spectral information.

S106, determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least squares discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

S106 specifically comprises:

according to the formula

Determining a discrimination coefficient; wherein Y is a discrimination coefficient, p₁,p₂,p₃,p₄,p₅,p₆,p₇,p₈,p₉,p₁₀And p₁₁All are fused spectral data to be identified;

when the absolute value Y-1 is less than 0.5, determining the gentiana rigescens slice sample to be identified as the domestic gentiana rigescens;

and when the absolute value Y-2 is less than 0.5, determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens.

The specific determination process of the partial least squares discriminant analysis equation is as follows:

s201, obtaining a gentiana rigescens section sample. Specifically, the gentiana rigescens plants of the Yunan province posterior rock mountain head are collected, 10 plants of each of the family species and the wild species are collected, and then 20 sliced samples are prepared, wherein the 20 sliced samples comprise 10 sliced samples of the family species gentiana rigescens and 10 sliced samples of the wild gentiana rigescens, and the preparation method is the same as that in S101.

S202, collecting spectrum data of different positions of each slice sample by using the LIBS spectrometer, selecting 16 different positions of one slice sample for spectrum collection in the embodiment of the invention, and collecting 320 LIBS spectra in total, wherein the collection system and the collection method are the same as S102.

S203, denoising 320 collected LIBS spectrums by adopting a wavelet denoising method, and reducing the LIBS spectrum variables from 22015 spectrum variables to 42 spectrum variables according to the corresponding relation between the wavelength and the elements, wherein the relation between the spectrum wavelength and the elements is shown in FIG. 2.

S204, placing 20 slices in an NIR spectrum system, collecting pixel point spectrum data by using an NIR spectrometer, wherein each collected spectrum has 256 spectrum variables, and collecting tens of thousands of pixel point spectra, and the collecting system and the collecting method are the same as S104.

And then denoising the collected tens of thousands of pixel point spectrums to obtain 160 family NIR spectrums and 160 wild NIR spectrums.

S205, processing 320 LIBS spectra and 320 NIR spectra by adopting the method in S105 to obtain LIBS scoresMatrix S_LIBS＝X_LIBS[V₁…… V₈]，S_LIBSIs 320 rows, 8 columns; obtaining a near infrared score matrix S_NIR＝X_NIR[V₁,V₂,V₃]，S_NIRIs 320 rows and 3 columns. The scored three-dimensional scatter plots for both data are shown in fig. 3 and fig. 4, where PC1, PC2, and PC3 represent the scoring vectors for the first three variances.

And S206, dividing the data after PCA fusion into a modeling set and a prediction set based on a K-Means algorithm (K-Means). K-Means is an unsupervised learning algorithm, samples are grouped into one class according to a standard measure function (sample mean square error), the score variables of the family samples are taken as an example, the K-Means groups the samples into one class according to the standard measure function (sample mean square error), the score variables of the family samples are divided into M classes, each class is composed of a plurality of score vectors, 2/3 vectors are randomly selected for each class of data and put into a modeling set, and the rest are put into a prediction set, so that the ratio of the modeling set to the prediction set is 2:1, and the uniformity of the modeling set to the prediction set is ensured.

S207, establishing a Partial Least squares discriminant Analysis (PLS-DA) equation, wherein Y is α₀+α₁p₁+α₂p₂+……+α₁₁p₁₁Determining coefficients in the equation according to the modeling set data and the partial least squares discriminant analysis equation, taking the obtained coefficients into the equation, and then verifying the coefficients by using the prediction set to obtain the partial least squares discriminant analysis equation:

in fig. 5, the first 52 spectra are spectra of the family gentiana rigescens, the second 52 spectra are spectra of the wild gentiana rigescens, and it can be seen from fig. 5 that the difference between the Y value of the discrimination degree and the absolute value of 1 is less than 0.5, and the family gentiana rigescens is identified, and the wild gentiana rigescens is identified, and otherwise, the family gentiana rigescens and the wild gentiana rigescens are not identified.

The near infrared band represents molecular composition and structure information of a substance, the LIBS spectrum reflects element composition and content information of the substance, the molecular composition and the content information can reflect difference between samples to a certain extent, but the represented information is limited, and the wild identification rate of the gentiana rigescens based on near infrared data and LIBS data is 85% and 81% respectively. Wherein, the accuracy of the modeling set and the prediction set is as follows:

data source	Modeling Collection accuracy (%)	Prediction set accuracy (%)
			LIBS data	93.06	81.73
Near infrared data	95.39	85.44
			Fusion data based on PCA	96.76	95.19

As can be seen from the table data, the accuracy of the prediction set after data fusion is greatly improved.

The invention also provides a multi-source spectrum fused gentiana rigescens species identification system, as shown in fig. 6, the system comprises:

and the section sample acquisition module 1 is used for acquiring a gentiana rigescens section sample to be identified.

The Gentiana rigescens LIBS spectral data acquisition module 2 is used for acquiring spectral data of different positions of the Gentiana rigescens slice sample to be identified by using a LIBS spectrometer to obtain the Gentiana rigescens LIBS spectral data to be identified; the LIBS spectrum data of the Gentiana rigescens to be identified comprise a plurality of LIBS spectra of the Gentiana rigescens to be identified.

And the screening LIBS spectral data determining module 3 is used for screening the wavelength of each LIBS spectrum to be identified in the LIBS spectral data according to the corresponding relation between the wavelength and the element, and determining to screen the LIBS spectral data to be identified.

The near infrared spectrum data acquisition module 4 for the gentiana rigescens to be identified acquires pixel point spectrum data of the gentiana rigescens sliced sample to be identified by using an NIR spectrometer to obtain near infrared spectrum data of the gentiana rigescens to be identified; the near infrared spectrum data of the gentiana rigescens to be identified comprises a plurality of near infrared spectrums of the gentiana rigescens to be identified.

And the fused spectral data acquisition module 5 is used for performing dimensionality reduction and fusion on the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting a principal component analysis method to obtain fused spectral data to be identified.

The identification module 6 is used for determining the type of the gentiana rigescens slice sample to be identified by adopting a partial least square discriminant analysis equation according to the fused spectral data to be identified; the species include the family species gentiana rigescens and the wild gentiana rigescens.

Preferably, the fused spectral data to be identified obtaining module 5 specifically includes:

and the LIBS spectrum matrix determining unit is used for determining the LIBS spectrum matrix of the gentiana rigescens to be identified according to the LIBS spectrum data of the gentiana rigescens to be identified.

A first projection vector set obtaining unit, configured to obtain a first projection vector set by using a principal component analysis method according to the LIBS spectral matrix of the gentiana rigescens to be identified; the first projection vector set comprises a plurality of first projection vectors, and each first projection vector corresponds to a different contribution rate value.

And the LIBS spectrum projection vector selecting unit is used for selecting the 8 first projection vectors with the maximum contribution value in the first projection vector set as LIBS spectrum projection vectors of the gentiana rigescens to be identified.

And the LIBS score matrix determining unit is used for determining the LIBS score matrix of the gentiana rigescens to be identified according to the LIBS spectrum matrix of the gentiana rigescens to be identified and the LIBS spectrum projection vector of the gentiana rigescens to be identified.

And the near infrared spectrum matrix determining unit is used for determining the near infrared spectrum matrix of the gentiana rigescens to be identified according to the near infrared spectrum data of the gentiana rigescens to be identified.

A second projection vector set obtaining unit, configured to obtain a second projection vector set by using a principal component analysis method according to the near infrared spectrum matrix of the gentiana rigescens to be identified; the second projection vector set comprises a plurality of second projection vectors, and each second projection vector corresponds to a different contribution rate value.

And the near infrared spectrum projection vector selecting unit is used for selecting the 3 second projection vectors with the maximum contribution values in the second projection vector set as the near infrared spectrum projection vectors of the gentiana rigescens to be identified.

And the near-infrared score matrix determining unit is used for determining the near-infrared score matrix of the gentiana rigescens to be identified according to the near-infrared spectrum matrix of the gentiana rigescens to be identified and the near-infrared spectrum projection vector of the gentiana rigescens to be identified.

And the spectrum data fusion unit is used for fusing the LIBS score matrix of the gentiana rigescens to be identified and the near-infrared score matrix of the gentiana rigescens to be identified to obtain fusion spectrum data to be identified.

Preferably, the identification module 6 specifically includes:

a discrimination coefficient determining unit for determining the discrimination coefficient based on the formulaDetermining a discrimination coefficient; wherein Y is a discrimination coefficient, p₁,p₂,p₃,p₄,p₅,p₆,p₇,p₈,p₉,p₁₀And p₁₁All are fused spectral data to be identified.

And the family gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the family gentiana rigescens when the absolute value Y-1 is less than 0.5.

And the wild gentiana rigescens determination unit is used for determining that the gentiana rigescens slice sample to be identified is the wild gentiana rigescens when the absolute value Y-2 is less than 0.5.

Preferably, the multisource spectrum fused gentiana rigescens species identification system further comprises:

and the area normalization processing module is used for respectively processing the LIBS spectral data of the screened gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified by adopting an area normalization method to obtain the LIBS spectral data of the gentiana rigescens to be identified and the near infrared spectral data of the gentiana rigescens to be identified.

The invention has the advantages that:

(1) laser Induced Breakdown Spectroscopy (LIBS) is a novel spectroscopic measurement technique that uses laser to excite a plasma and collect the emission spectra of atoms or ions in the cooling process of the plasma to obtain the elemental composition and content information of the measured object. The near infrared spectrum technology is used for indirectly obtaining the component information of a substance by collecting a section of spectrum containing the molecular composition and the structure of a sample. The two technologies have the advantages of being a rapid nondestructive testing technology, low in requirement on the form of a sample, capable of meeting the requirements of two experiments simultaneously in the tabletting form, simple, rapid and efficient.

(2) The near infrared band represents molecular composition and structure information of a substance, the LIBS spectrum reflects element composition and content information of the substance, the LIBS spectrum and the LIBS spectrum can reflect difference between samples to a certain extent, but the represented information is limited.

(3) By fusing the principal component score vectors of the two data, a Gentiana rigescens wild rapid discrimination model combining LIBS and NIR technology, namely a partial least square discrimination analysis equation, can be established, the prediction effect of the model can be greatly improved, the accuracy of a prediction set is as high as 95%, and the method has the characteristics of high detection speed, few required samples, no pollution, high detection precision, strong reliability and low artificial interference factors. The method provides a technical means for distinguishing the wild gentiana rigescens, and meets the requirements of rapid market detection, simple sample processing, high efficiency, accuracy and the like.

The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.

The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.