CN105548037A - Method for nondestructively detecting raw traditional Chinese medicinal materials - Google Patents
Method for nondestructively detecting raw traditional Chinese medicinal materials Download PDFInfo
- Publication number
- CN105548037A CN105548037A CN201610020565.5A CN201610020565A CN105548037A CN 105548037 A CN105548037 A CN 105548037A CN 201610020565 A CN201610020565 A CN 201610020565A CN 105548037 A CN105548037 A CN 105548037A
- Authority
- CN
- China
- Prior art keywords
- sample
- pixel
- crude drug
- chinese medicine
- pls
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
Abstract
The invention belongs to the field of medicinal material identification, and discloses a method for nondestructively detecting whether raw traditional Chinese medicinal materials are true or false. The problems that detecting modes in the traditional technology are complex in operation and high in detecting cost, and need to consume samples are solved. The method includes the steps that a, positive samples and negative samples are selected; b, hyperspectral images of the samples are collected; c, the hyperspectral images are processed, and relative characteristics are extracted with a principal component analysis method; d, a PLS-DA prediction model is built; e, to-be-detected objects are detected through the PLS-DA prediction model; f, whether the to-be-detected objects are true or false is judged according to the output of the model. The method is suitable for rapidly and accurately identifying the raw traditional Chinese medicinal materials.
Description
Technical field
The invention belongs to field of crude medicine identification, be specifically related to the method for the Non-Destructive Testing Chinese medicine crude drug true and false.
Background technology
Cordyceps sinensis Cordycepssinensis (Berkeley) Saccardo is that section ergot fungus cordyceps sinensis bacterium colonizes in the entomogenous fungi that Lepidoptera Hepialidae Genus Hepialus (Hepialus) larva is formed, for the compound of the stroma of Clavicipitaceae Cordyceps sinensis fungus Cordyceps sinensis and the larva corpse of host Lepidoptera Hepialidae insect bat moth thereof, property sweet flat, invigorate the lung and the kidney, hemostasis and phlegm, breathe heavily for chronic cough void, phthisical cough is spat blood, impotence and seminal emission, soreness of waist and knee joint is the traditional rare traditional Chinese medicine of China.
Cordyceps sinensis growing environment is special, is mainly distributed in mesophorbium, coryphile and the alpine scrub of China's Qinghai-Tibet Platean height above sea level 3000 ~ 5000m.Qinghai, Tibet, Sichuan, Gansu and Yunnan are the major production areas of Cordyceps sinensis, and the seed output and quality in Qinghai occupies first of each provinces and regions, and cajaput, Golog etc. are Cordyceps of Qinghai Province main product ground.Cordyceps sinensis still can not carry out artificial culture at present, and the distributed areas of wild cordyceps are narrow and small, and natural parasitic rate is low, harsh to requirement for environmental conditions, in addition ecological disruption in recent years and predation formula are excavated, and the output of Cordyceps sinensis is declined year by year, and price constantly rises.Therefore, market exists serious mix the fakement phenomena such as puppet, weightening finish.Except Cordyceps sinensis, other Cordyceps sinensis fungus colonize in the complex that insect bodies is formed and are also referred to as " Chinese caterpillar fungus ", and some of them are often used as the adulterant of Cordyceps sinensis, and proterties is similar and be difficult to distinguish.
Therefore set up Chinese medicine crude drug discrimination method fast and accurately, to standard market, ensure that quality of medicinal material has very important significance.Authentication method in conventional art is mainly tested according to the dependence experience such as source, kind form, proterties, microscopic features, physics and chemistry discriminating, inspection, assay of Chinese medicine or exact instrument, especially the method such as microscopic features, physics and chemistry discriminating, inspection, assay all needs consumption of raw material, is difficult to identify by the feature on medicinal material surface.
In recent years, the detection mode of DNA molecular marker is used to carry out the detection of medicinal material and differentiate also common reporter.But, this detection mode complex operation, reagent and testing cost higher, preparation of samples program is complicated, is difficult to promote.
The most important thing is, current various detection methods all need to consume sample, therefore can only sampling Detection, can not accomplish the detection of all over products, therefore be difficult to the needs meeting Chinese Medicine Industry.
Summary of the invention
The technical problem to be solved in the present invention is: a kind of method proposing Non-Destructive Testing Chinese medicine crude drug true and false, solves detection mode in conventional art and there is the problem that complex operation, testing cost are high, need to consume sample.
The present invention solves the problems of the technologies described above adopted technical scheme: the method for Non-Destructive Testing Chinese medicine crude drug, comprises the following steps:
Set up certain Chinese medicine crude drug PLS-DA forecast model, when carrying out this kind Chinese medicine crude drug and differentiating, detected sample is placed in high spectrum reflection image capturing system, described high spectrum reflection image capturing system is utilized to gather the hyperspectral information of described detected sample, and input this kind Chinese medicine crude drug PLS-DA forecast model set up and predict after image procossing is carried out to the hyperspectral information gathered, according to prediction output valve, this sample is differentiated.
Concrete, when carrying out this kind Chinese medicine crude drug and differentiating, the pixel spectral information that the hyperspectral information utilizing high spectrum reflection image capturing system to gather is detected sample.
Concrete, the described step setting up certain Chinese medicine crude drug PLS-DA forecast model comprises:
A1. this kind crude drug sample is chosen respectively after pretreatment as positive sample sets;
A2. the sample easily obscured with this kind crude drug is chosen respectively after pretreatment as bearing sample sets;
A3. respectively the positive sample in positive sample sets and the negative sample in negative sample sets are placed in high spectrum reflection image capturing system, utilize the hyperspectral information of described high spectrum reflection image capturing system collected specimens;
A4. spectral signature is extracted after image procossing being carried out to the hyperspectral information gathered, by the spectral signature input database extracted;
A5. steps A 3-A4 is repeated, until complete extraction and the typing of the spectral signature of all samples in positive and negative sample sets;
A6. the spectral signature of random selecting some samples sets up PLS-DA forecast model.
Concrete, when setting up certain Chinese medicine crude drug PLS-DA forecast model, in steps A 3 and A4, the sample hyperspectral information of described high spectrum reflection image capturing system collection is the averaged spectrum information of this sample.
Concrete, in steps A 1, choose this kind crude drug sample of national Different sources respectively after pretreatment as positive sample sets; In steps A 2, choose on the market the multiple sample easily obscured with this kind crude drug respectively after pretreatment as bearing sample sets.
As preferably, described Chinese medicine crude drug is the former grass of Cordyceps sinensis.
In steps A 1, described in choose national Different sources this kind crude drug sample specifically comprise: be selected from national main producing region Yushu district, Qinghai, Qinghai Golog, Hainan, Qinghai, Qinghai Province east, Sichuan, Tibet, Gansu, the regional Cordyceps sinensis sample in eight, Yunnan.
In steps A 2, the described sample easily obscured with this kind crude drug includes but are not limited to: numb back, Cordyceps militaris, sub-fragrant excellent, liangshan cordyceps herb and Xinjiang Chinese caterpillar fungus etc.
Concrete, in steps A 1 and A2, described pre-service refers to successively through overdrying brush, cleaning, 40 DEG C of low temperature dryings.
Concrete, in steps A 4, maximum according to positive and negative sample spectra variance after image procossing is carried out to the hyperspectral information gathered
Change principle and extract spectral signature.
Concrete, in steps A 6, the method for the spectral signature of random selecting some samples is:
The spectral signature of the positive and negative sample of random selecting equal proportion.
Concrete, describedly carry out differentiating that concrete grammar is to this sample according to prediction output valve: if certain pixel of this sample through its output valve of PLS-DA model prediction between-1.5 ~ 0.45, then judge that this pixel is as pseudo-medicinal material pixel; If certain pixel of this sample between 0.55 ~ 1.5, then judges that this pixel is as true medicinal material pixel through its output valve of PLS-DA model prediction; If certain pixel of this sample between 0.45 ~ 0.55, then judges that this pixel is as unknown classification pixel through its output valve of PLS-DA model prediction;
Finally also calculate the ratio in the total pixel number of this sample by the number of system difference programming count true medicinal material pixel, pseudo-medicinal material pixel, unknown classification pixel, differentiate the true and false of this sample according to ratio distribution situation.
The described true and false differentiating this sample according to ratio distribution situation, concrete grammar is:
If ratio >=60% of the number of true medicinal material pixel in the total pixel number of this sample, then judge that this sample is as true.
Concrete, described high spectrum reflection image capturing system adopts mercury cadmium telluride two-dimensional array detector, and light source is quartz halogen lamp; Spectra collection scope is short infrared wave band 940 – 2537nm, and pixel is 320 × 256, pixel size 150 μm × 150 μm, adopts visual field to be 50mm camera lens; Scan mode is high speed push-broom type high light spectrum image-forming, pushes away and sweeps speed 3mm/s, picking rate 100fps.
Concrete, utilize described high spectrum reflection image capturing system to gather the hyperspectral information of this sample, concrete steps are as follows:
1. gather and obtain the continuous spectrum curve of this sample (m × n) individual pixel under k wave band, the spectral signal response that each wave band is corresponding is I
k, k=1,2 ... K;
2. utilize the light intensity value of standard white plate uncalibrated image, calculate the relative light intensity value of every width image high spectrum reflection image under a kth wave band
wherein
for the relative high light value of Cordyceps sinensis high spectrum reflection image each under a kth wave band; I
kfor the light intensity value of Cordyceps sinensis high spectrum reflection image each under a kth wave band;
for the light intensity value of standard white plate high spectrum reflection image under a kth wave band; D
kfor the complete black uncalibrated image light intensity value gathered under a kth wave band;
3. the relative light intensity value calculated is changed through A/D, be converted to the curve of spectrum.
Concrete, extracting composition correlated characteristic after the described hyperspectral information to gathering carries out image procossing, specifically comprising:
After average centralization conversion is carried out respectively by spectrum dimension to each pel data of collection, carry out PCA conversion, retain destination object place pixel spectrum and positional information, deduct useless background cell;
With each sample for unit, calculate the averaged spectrum obtaining each sample, and the data set that averaged spectrum is formed is carried out successively to Savitsky – Golay is level and smooth, standard just too variate calibration, average centralization process and PCA conversion, when ensureing cumulative variance >=90%, get the characteristic information of top n major component.
As further optimization, N=3.
Concrete, the method for the background cell that described deduction is useless is:
Man-machine interactive selects the ROI pixel of sample in PCA score space, calculates the Euclidean distance between backdrop pels and sample pixel and is shown as represented as histograms, finds the threshold value that background significantly can be separated with sample pixel, deletes useless background cell.
The invention has the beneficial effects as follows: the present invention adopts the averaged spectrum of whole root herb to set up PLS-DA forecast model, prediction adopts each pixel spectrum of medicinal material as input, while completing prediction, complete and remain medicinal material profile, the spatial positional informations such as surface texture featur.Accurately, directly perceived, easy, rapidly medicinal material is differentiated.Have objective quantification, result accurate, easy and simple to handle, test the plurality of advantages such as rapid, with low cost.
Accompanying drawing explanation
Fig. 1 is the method flow diagram of the Non-Destructive Testing Chinese medicine crude drug true and false in embodiment;
Fig. 2 is the Golog Cordyceps sinensis PLS-DA prognostic chart in embodiment 1;
Fig. 3 is the sub-fragrant excellent PLS-DA prognostic chart of puppet grass in embodiment 2.
Embodiment
The present invention is intended to a kind of method proposing Non-Destructive Testing Chinese medicine crude drug true and false, solves detection mode in conventional art and there is the problem that complex operation, testing cost are high, need to consume sample.
As shown in Figure 1, it comprises the method for the Non-Destructive Testing Chinese medicine crude drug true and false in the present invention:
A. positive and negative sample is chosen;
B. the high spectrum image of collected specimens;
C. high spectrum image processed and utilize principal component analysis (PCA) to extract correlated characteristic;
D. PLS-DA forecast model is set up;
E. utilize PLS-DA forecast model to treat detected object to detect;
F. the true and false differentiating object to be detected is exported according to model.
In specific implementation, to detect the former grass of Cordyceps sinensis, the present invention program comprises following components:
One, the step of the former careless PLS-DA forecast model of Cordyceps sinensis is set up:
1. preparation of samples:
1.1 choose main producing region Yushu district, Qinghai, Qinghai Golog, Hainan, Qinghai, Qinghai Province east, Sichuan, Tibet, Gansu, the former careless sample of 80 parts, area, eight, Yunnan Cordyceps sinensis, as positive sample sets after dry brush, cleaning, 40 DEG C of low temperature dryings.
1.2 30 parts, the samples choosing the former grass of commercially available similar Cordyceps sinensis (easily with the former grass of Cordyceps sinensis obscure), after dry brush, cleaning, 40 DEG C of low temperature dryings, sample sets is born in conduct.
2. image acquisition:
High spectrum reflection image capturing system of the present invention is the SisuCHEMA laboratory EO-1 hyperion scanner that Finland SPECIM produces, and it adopts mercury cadmium telluride two-dimensional array detector, and light source is quartz halogen lamp.Spectra collection scope is short infrared wave band (SWIR, 940 – 2537nm), and pixel is 320 (space) × 256 (spectrum), pixel size 30 μm × 30 μm, adopts visual field to be 10mm camera lens.Scan mode is high speed push-broom type high light spectrum image-forming, pushes away and sweeps speed 3mm/s, picking rate 100fps.
Respectively the sample prepared is placed on sample stage successively, require should there be the gap being not less than 5mm between every root Cordyceps sinensis or pseudo-careless sample, overlap or adhesion can not be had; High spectrum reflection image capturing system enters from sample the high spectrum image information that visual field starts each pixel of collected specimens:
1. gather and obtain the continuous spectrum curve of this sample (m × n) individual pixel under k wave band, the spectral signal response that each wave band is corresponding is I
k, k=1,2 ... K;
2. utilize the light intensity value of standard white plate uncalibrated image, calculate the relative light intensity value of every width image high spectrum reflection image under a kth wave band
wherein
for the relative high light value of Cordyceps sinensis high spectrum reflection image each under a kth wave band; I
kfor the light intensity value of Cordyceps sinensis high spectrum reflection image each under a kth wave band;
for the light intensity value of standard white plate high spectrum reflection image under a kth wave band; D
kfor the complete black uncalibrated image light intensity value gathered under a kth wave band;
3. the relative light intensity value calculated is changed through A/D, be converted to the curve of spectrum.
3. image procossing:
For noise decrease interference, data under 940 – 1000nm containing more noise and 2469 – 2537nm wave bands can be deleted; Then, after average centralization conversion being carried out respectively by spectrum dimension to each pel data, carry out PCA conversion, retain destination object place pixel spectrum and positional information, deduct useless background cell, obtain effective pixel spectrum;
In order to eliminate the interference such as baseline wander, light scattering, noise, sample surface morphology difference, Savitsky – Golay level and smooth (window value 11, polynomial expression exponent number is 3), standard normal variable correction (SNV) and average centralization process need be carried out to effective pixel spectrum; Then PCA (principal component analysis (PCA)) conversion is carried out to pretreated data.
4. model is set up:
When ensureing cumulative variance >=90%, the characteristic information getting top n major component (generally, first three major component can meet cumulative variance >=90%) sets up PLS-DA (partial least squares discriminant analysis) forecast model; Constantly improve and the reliability of verification model to keep model, the sample composing training collection of 2/3 quantity can be extracted from the positive and negative sample sets prepared, remaining 1/3 as test set, utilize the continuous sophisticated model of the sample in training set, utilize the sample in test set to verify model prediction accuracy.
Two, when carrying out the former grass of Cordyceps sinensis and differentiating, detected sample is placed in high spectrum reflection image capturing system, described high spectrum reflection image capturing system is utilized to gather the hyperspectral information of described detected sample, and input to the former careless PLS-DA forecast model of the Cordyceps sinensis set up in steps A and predict after image procossing is carried out to the hyperspectral information gathered, according to prediction output valve, this sample is differentiated.
Concrete discrimination method is: if certain pixel of this sample through its output valve of PLS-DA model prediction between-1.5 ~ 0.45, then judge that this pixel is as the careless pixel of puppet; If certain pixel of this sample between 0.55 ~ 1.5, then judges that this pixel is as very careless pixel through its output valve of PLS-DA model prediction; If certain pixel of this sample between 0.45 ~ 0.55, then judges that this pixel is as unknown classification pixel through its output valve of PLS-DA model prediction;
Finally also calculate the ratio in the total pixel number of this sample by the number of the very careless pixel of system difference programming count, pseudo-careless pixel, unknown classification pixel, differentiate the true and false of this sample according to ratio distribution situation.
With several specific embodiment, the PLS-DA forecast model set up in the present invention is verified below:
Embodiment 1:
This example is differentiated for the former grass of Golog Cordyceps sinensis, comprises following performing step:
1. preparation of samples:
By the former grass of Golog Cordyceps sinensis that collects after overdrying brush, cleaning, 40 DEG C of low temperature dryings, it is neatly put in sample panel, requires should there be the gap being not less than 5mm between every root Cordyceps sinensis sample, to guarantee can not there be overlapping adhesion between every root grass.
2. sample high-spectral data gathers:
Under sample being sent to EO-1 hyperion camera lens, by setup parameter, (spectra collection scope is short infrared wave band (SWIR, 940 – 2537nm), pixel is 320 (space) × 256 (spectrum), pixel size 30 μm × 30 μm, adopts visual field to be 10mm camera lens.Scan mode is high speed push-broom type high light spectrum image-forming, pushes away and sweeps speed 3mm/s, picking rate 100fps.) scan, preserve high spectrum image and the spectral information of sample.Image procossing adopts the built-in monochrome scale material of camera to carry out automatic calibration to image by the Evince software of Umbio company of Sweden, then light intensity value is converted to spectral absorption curve through A/D.
3. image procossing:
First delete containing data under 940 – 1000nm of more noise and 2469 – 2537nm wave bands.After carrying out average centralization conversion to each pel data of collection respectively by spectrum dimension, carry out PCA conversion, man-machine interactive selects the Cordyceps sinensis ROI in PCA score space (RegionofInteresting area-of-interest) pixel; Euclidean distance between further calculating backdrop pels and Cordyceps sinensis pixel, and be shown as represented as histograms, find the threshold value that background significantly can be separated with Cordyceps sinensis pixel, delete useless backdrop pels.
4.PLS-DA predicts:
The each pixel spectrum 3rd step being obtained every root Cordyceps sinensis carry out successively Savitsky – Golay level and smooth (window value 11, polynomial expression exponent number is 3), standard just too variate calibration (SNV) and average centralization process to eliminate the interference such as baseline wander, light scattering, noise, sample surface morphology difference.
Each for pretreated Cordyceps sinensis pixel spectrum is brought into the PLS-DA model established and carry out prediction and calculation (should ensure that Pixel domain position does not change).When the PLS model predication value of certain pixel is-1.5 ~ 0.45, judges that this pixel is as the careless pixel of puppet, in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as blueness; When the predicted value of certain pixel is 0.55 ~ 1.5, judges that this pixel is as very careless pixel, in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as green; When the predicted value of certain pixel is 0.45 ~ 0.55, is judged to be unknown classification (being likely silt or other pollutants), in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as redness.
After completing the prediction of all pixels, add up the ratio of the sum of very careless pixel, the sum of pseudo-careless pixel, the sum of unknown classification and the corresponding shared total pixel number of this sample; Predict the outcome as shown in table 1:
Table 1: Golog Cordyceps sinensis EO-1 hyperion pixel PLS-DA predicts the outcome
As can be seen from table 1.1, in the prediction of Golog Cordyceps sinensis, the ratio being predicted as Cordyceps sinensis pixel accounts for 83.25%, and pseudo-careless pixel scale is 11.76%.But consider the globality (it is very careless for can not there is a part, and a part is pseudo-grass) of Cordyceps sinensis, this Cordyceps sinensis can be judged to be very grass (RED sector may be other materials such as earth).
In addition, also from the PLS-DA prognostic chart Fig. 2, intuitively can find out that this sample overwhelming majority predicted picture is for green, smaller portions predicted picture is blue, and very fraction predicted picture is red, equally based on the consideration of Cordyceps sinensis globality, can judge that this sample is as very careless.
As can be seen from above, this model has image, the feature such as directly perceived, simple, stable, has good fault-tolerant ability, can be applicable to the automatic discriminating of the Cordyceps sinensis raw material true and false in actual production, reduce the interference of human factor, have objective, authenticity.
Embodiment 2:
This example is differentiated for the sub-fragrant rod of pseudo-grass, comprises following performing step:
1. preparation of samples:
After overdrying brush, cleaning, 40 DEG C of low temperature dryings, it is neatly put excellent for the sub-perfume of puppet grass collected in sample panel, require should there be the gap being not less than 5mm, to guarantee having overlapping adhesion between the sub-fragrant excellent sample of the pseudo-grass of every root.
2. sample high-spectral data gathers:
Under sample being sent to EO-1 hyperion camera lens, scan by setup parameter (setting with in embodiment 1), preserve high spectrum image and the spectral information of sample.Image procossing adopts the built-in monochrome scale material of camera to carry out automatic calibration to image by the Evince software of Umbio company of Sweden, and then light intensity value is converted to spectral absorption curve through A/D.
3. image procossing:
First delete containing data under 940 – 1000nm of more noise and 2469 – 2537nm wave bands.After carrying out average centralization conversion to each pel data of collection respectively by spectrum dimension, carry out PCA conversion, man-machine interactive selects the ROI pixel of sub-fragrant rod in PCA score space; Euclidean distance between the fragrant excellent pixel of further calculating backdrop pels and Asia, and be shown as represented as histograms, find the threshold value that fragrant to background and Asia excellent pixel significantly can be separated, delete useless backdrop pels.
4.PLS-DA predicts:
The each pixel spectrum 3rd step being obtained the sub-fragrant rod grass of every root carry out successively Savitsky – Golay level and smooth (window value 11, polynomial expression exponent number is 3), standard just too variate calibration (SNV) and average centralization process to eliminate the interference such as baseline wander, light scattering, noise, sample surface morphology difference.
Bring the careless each pixel spectrum of pretreated sub-fragrant rod into establish PLS-DA model and carry out prediction and calculation (should ensure that Pixel domain position does not change).When the PLS model predication value of certain pixel is-1.5 ~ 0.45, judges that this pixel is as the careless pixel of puppet, in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as blueness; When the predicted value of certain pixel is 0.55 ~ 1.5, judges that this pixel is as very careless pixel, in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as green; When the predicted value of certain pixel is 0.45 ~ 0.55, is judged to be unknown classification (being likely silt or other pollutants), in order to more intuitively manifest, at this pixel correspondence position, this pixel can be labeled as redness.
After completing the prediction of all pixels, add up the ratio of the sum of very careless pixel, the sum of pseudo-careless pixel, the sum of unknown classification and the corresponding shared total pixel number of this sample; Predict the outcome as shown in table 2:
Table 2: Cordyceps hawkesii Gary EO-1 hyperion pixel PLS-DA predicts the outcome
As can be seen from Table 2, the ratio being predicted as sub-fragrant excellent pixel accounts for 68.54%, and Cordyceps sinensis pixel scale is 8.99%.But consider the globality (it is very careless for can not there is a part, and a part is pseudo-grass) of the pseudo-grass of the fragrant rod in Asia, this grass can be judged to be pseudo-former grass (RED sector may be other materials such as earth).
In addition, also from the PLS-DA prognostic chart Fig. 3, intuitively can find out that this sample overwhelming majority predicted picture is for blue, smaller portions predicted picture is red, and very fraction predicted picture is green, equally based on the consideration of globality, can judge that this sample is as pseudo-grass.
As can be seen from above, this model has image, the feature such as directly perceived, simple, stable, has good fault-tolerant ability, can be applicable to the automatic discriminating of the Cordyceps sinensis raw material true and false in actual production, reduce the interference of human factor, have objective, authenticity.
Claims (10)
1. the method for Non-Destructive Testing Chinese medicine crude drug, is characterized in that, comprising:
Set up certain Chinese medicine crude drug PLS-DA forecast model, when carrying out this kind Chinese medicine crude drug and differentiating, detected sample is placed in high spectrum reflection image capturing system, described high spectrum reflection image capturing system is utilized to gather the hyperspectral information of described detected sample, and input this kind Chinese medicine crude drug PLS-DA forecast model set up and predict after image procossing is carried out to the hyperspectral information gathered, according to prediction output valve, this sample is differentiated; Hyperspectral information is the pixel spectral information of detected sample.
2. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 1, is characterized in that,
The described step setting up certain Chinese medicine crude drug PLS-DA forecast model comprises:
A1. this kind crude drug sample is chosen respectively after pretreatment as positive sample sets;
A2. the sample easily obscured with this kind crude drug is chosen respectively after pretreatment as bearing sample sets;
A3. respectively the positive sample in positive sample sets and the negative sample in negative sample sets are placed in high spectrum reflection image capturing system, utilize the hyperspectral information of described high spectrum reflection image capturing system collected specimens;
A4. spectral signature is extracted after image procossing being carried out to the hyperspectral information gathered, by the spectral signature input database extracted;
A5. steps A 3-A4 is repeated, until complete extraction and the typing of the spectral signature of all samples in positive and negative sample sets;
A6. the spectral signature of the positive and negative sample of random selecting equal proportion sets up PLS-DA forecast model;
In steps A 3 and A4, the sample hyperspectral information of described high spectrum reflection image capturing system collection is the averaged spectrum information of this sample.
3. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 2, is characterized in that,
Described Chinese medicine crude drug is the former grass of Cordyceps sinensis.
4. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, is characterized in that,
In steps A 1, choose this kind crude drug sample of national Different sources respectively after pretreatment as positive sample sets; Described this kind crude drug sample choosing national Different sources specifically comprises: be selected from national main producing region Yushu district, Qinghai, Qinghai Golog, Hainan, Qinghai, Qinghai Province east, Sichuan, Tibet, Gansu, the regional Cordyceps sinensis sample in eight, Yunnan;
In steps A 2, choose the multiple sample easily obscured with this kind crude drug on the market and, respectively after pretreatment as negative sample sets, comprising: be numb back, Cordyceps militaris, sub-fragrant excellent, liangshan cordyceps herb and Xinjiang Chinese caterpillar fungus.
5. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, is characterized in that,
In steps A 1 and A2, described pre-service refers to successively through overdrying brush, cleaning, 40 DEG C of low temperature dryings.
6. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, is characterized in that,
In steps A 4, after image procossing is carried out to the hyperspectral information gathered, extract spectral signature according to positive and negative sample spectra maximum variance principle.
7. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, it is characterized in that, described according to prediction output valve to this sample carry out differentiate concrete grammar be: if certain pixel of this sample through its output valve of PLS-DA model prediction between-1.5 ~ 0.45, then judge that this pixel is as pseudo-medicinal material pixel; If certain pixel of this sample between 0.55 ~ 1.5, then judges that this pixel is as true medicinal material pixel through its output valve of PLS-DA model prediction; If certain pixel of this sample between 0.45 ~ 0.55, then judges that this pixel is as unknown classification pixel through its output valve of PLS-DA model prediction;
Finally also calculate the ratio in the total pixel number of this sample by the number of system difference programming count true medicinal material pixel, pseudo-medicinal material pixel, unknown classification pixel, the true and false of this sample is differentiated: if ratio >=60% of the number of true medicinal material pixel in the total pixel number of this sample, then judge that this sample is as true according to ratio distribution situation.
8. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, is characterized in that,
Described high spectrum reflection image capturing system adopts mercury cadmium telluride two-dimensional array detector, and light source is quartz halogen lamp;
The spectra collection scope of described high spectrum reflection image capturing system is short infrared wave band 940 – 2537nm;
The pixel of described high spectrum reflection image capturing system is 320 × 256, pixel size 150 μm × 150 μm, adopts visual field to be 50mm camera lens;
Scan mode is high speed push-broom type high light spectrum image-forming, pushes away and sweeps speed 3mm/s, picking rate 100fps.
9. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, it is characterized in that, utilize described high spectrum reflection image capturing system to gather the hyperspectral information of this sample, concrete steps are as follows:
1. gather and obtain the continuous spectrum curve of this sample (m × n) individual pixel under k wave band, the spectral signal response that each wave band is corresponding is I
k, k=1,2 ... K;
2. utilize the light intensity value of standard white plate uncalibrated image, calculate the relative light intensity value of every width image high spectrum reflection image under a kth wave band
wherein
for the relative high light value of Cordyceps sinensis high spectrum reflection image each under a kth wave band; I
kfor the light intensity value of Cordyceps sinensis high spectrum reflection image each under a kth wave band;
for the light intensity value of standard white plate high spectrum reflection image under a kth wave band; D
kfor the complete black uncalibrated image light intensity value gathered under a kth wave band;
3. the relative light intensity value calculated is changed through A/D, be converted to the curve of spectrum.
10. the method for Non-Destructive Testing Chinese medicine crude drug as claimed in claim 3, is characterized in that, extracts spectral signature, specifically comprise described in steps A 4 to the hyperspectral information gathered after carrying out image procossing:
After average centralization conversion is carried out respectively by spectrum dimension to each pel data of collection, carry out PCA conversion, retain destination object place pixel spectrum and positional information, deduct useless background cell;
With each sample for unit, calculate the averaged spectrum obtaining each sample, and the data set that averaged spectrum is formed is carried out successively to Savitsky – Golay is level and smooth, standard just too variate calibration, average centralization process and PCA conversion, when ensureing cumulative variance >=90%, get the spectral signature information of front 3 major components;
The method of the background cell that described deduction is useless is:
Man-machine interactive selects the ROI pixel of sample in PCA score space, calculates the Euclidean distance between backdrop pels and sample pixel and is shown as represented as histograms, finds the threshold value that background significantly can be separated with sample pixel, deletes useless background cell.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610020565.5A CN105548037A (en) | 2015-01-14 | 2016-01-13 | Method for nondestructively detecting raw traditional Chinese medicinal materials |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510018096 | 2015-01-14 | ||
| CN2015100180969 | 2015-01-14 | ||
| CN201610020565.5A CN105548037A (en) | 2015-01-14 | 2016-01-13 | Method for nondestructively detecting raw traditional Chinese medicinal materials |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105548037A true CN105548037A (en) | 2016-05-04 |
Family
ID=55827378
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610020565.5A Pending CN105548037A (en) | 2015-01-14 | 2016-01-13 | Method for nondestructively detecting raw traditional Chinese medicinal materials |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105548037A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106404689A (en) * | 2016-09-23 | 2017-02-15 | 广东食品药品职业学院 | Identification method of components of exocarpium |
| CN109164059A (en) * | 2018-10-12 | 2019-01-08 | 四川农业大学 | A method of using infrared spectrum technology Rapid identification conyza blinii and its easily obscuring Chinese herbal medicine |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1702684A (en) * | 2005-04-06 | 2005-11-30 | 北京航空航天大学 | Strong noise image characteristic points automatic extraction method |
| US20060290933A1 (en) * | 2005-06-27 | 2006-12-28 | Paige Holm | System and method for monitoring plant conditions |
| CN101881726A (en) * | 2010-06-18 | 2010-11-10 | 北京农业智能装备技术研究中心 | Nondestructive detection method for comprehensive character living bodies of plant seedlings |
| CN102081039A (en) * | 2010-08-17 | 2011-06-01 | 江苏大学 | Environment-controllable hyperspectral image detecting device for crop nutrition and moisture |
| CN102495005A (en) * | 2011-11-17 | 2012-06-13 | 江苏大学 | Method for diagnosing crop water deficit through hyperspectral image technology |
| CN103411973A (en) * | 2013-09-03 | 2013-11-27 | 西北农林科技大学 | Method for measuring anthocyanin content in wine grape pericarp based on hyperspectrum |
-
2016
- 2016-01-13 CN CN201610020565.5A patent/CN105548037A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1702684A (en) * | 2005-04-06 | 2005-11-30 | 北京航空航天大学 | Strong noise image characteristic points automatic extraction method |
| US20060290933A1 (en) * | 2005-06-27 | 2006-12-28 | Paige Holm | System and method for monitoring plant conditions |
| CN101881726A (en) * | 2010-06-18 | 2010-11-10 | 北京农业智能装备技术研究中心 | Nondestructive detection method for comprehensive character living bodies of plant seedlings |
| CN102081039A (en) * | 2010-08-17 | 2011-06-01 | 江苏大学 | Environment-controllable hyperspectral image detecting device for crop nutrition and moisture |
| CN102495005A (en) * | 2011-11-17 | 2012-06-13 | 江苏大学 | Method for diagnosing crop water deficit through hyperspectral image technology |
| CN103411973A (en) * | 2013-09-03 | 2013-11-27 | 西北农林科技大学 | Method for measuring anthocyanin content in wine grape pericarp based on hyperspectrum |
Non-Patent Citations (2)
| Title |
|---|
| 双利合谱: "基于高光谱成像技术的虫草粉末无损鉴定的实验探索", 《双利合谱》 * |
| 高海龙: "基于透射和反射高光谱成像技术的马铃薯缺陷检测方法研究", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106404689A (en) * | 2016-09-23 | 2017-02-15 | 广东食品药品职业学院 | Identification method of components of exocarpium |
| WO2018054091A1 (en) * | 2016-09-23 | 2018-03-29 | 广州道地南药技术研究有限公司 | Method for identifying components of exocarpium |
| CN109164059A (en) * | 2018-10-12 | 2019-01-08 | 四川农业大学 | A method of using infrared spectrum technology Rapid identification conyza blinii and its easily obscuring Chinese herbal medicine |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wei et al. | Cloud detection for Landsat imagery by combining the random forest and superpixels extracted via energy-driven sampling segmentation approaches | |
| CN102435713B (en) | Automatic detection system for quality of traditional Chinese medicine | |
| CN109765194A (en) | Fructus lycii Production area recognition method based on high light spectrum image-forming technology | |
| Yuan et al. | Selecting key wavelengths of hyperspectral imagine for nondestructive classification of moldy peanuts using ensemble classifier | |
| Wang et al. | Landscape-level vegetation classification and fractional woody and herbaceous vegetation cover estimation over the dryland ecosystems by unmanned aerial vehicle platform | |
| CN105699343A (en) | Method for nondestructively identifying authenticity of agilawood | |
| CN109827957B (en) | Rice leaf SPAD value estimation method and system based on computer vision | |
| Li et al. | GIS-based detection of grain boundaries | |
| CN105466921A (en) | Simultaneous detection method of many samples | |
| CN105628708A (en) | Quick nondestructive testing method for multi-parameter quality of south Xinjiang red dates | |
| Li et al. | SPM-IS: An auto-algorithm to acquire a mature soybean phenotype based on instance segmentation | |
| CN103528967A (en) | Hyperspectral image based overripe Lonicera edulis fruit identification method | |
| CN105115910A (en) | Method for detecting protein content distribution in peanuts based on hyperspectral imaging technology | |
| CN105527241A (en) | Non-destructive method for detecting authenticity of raw cordyceps sinensis | |
| Raut et al. | Assessment of fruit maturity using digital image processing | |
| CN106845366A (en) | Sugarcane coverage automatic testing method based on image | |
| CN105606544A (en) | Nondestructive detection method of insect bodies of Cordyceps sinensis | |
| Tian et al. | The detection system for greenhouse tomato disease degree based on android platform | |
| Magiera et al. | Relating canopy reflectance to the vegetation composition of mountainous grasslands in the Greater Caucasus | |
| CN105548037A (en) | Method for nondestructively detecting raw traditional Chinese medicinal materials | |
| CN105466865A (en) | Nondestructive ophiocordyceps sinensis stroma adulteration detection method | |
| CN116862456B (en) | Traditional Chinese medicine production monitoring control system and method based on image processing | |
| CN106338488A (en) | Method for fast undamaged determination of transgenic soybean milk powder | |
| CN109308464A (en) | A kind of rice tassel recognition methods based on digital camera images | |
| CN107255624A (en) | A kind of fast and convenient method for determining protein content in fan's power mushroom entity |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160504 |
|
| RJ01 | Rejection of invention patent application after publication |