CN109270009A - A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm - Google Patents
A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm Download PDFInfo
- Publication number
- CN109270009A CN109270009A CN201811360899.2A CN201811360899A CN109270009A CN 109270009 A CN109270009 A CN 109270009A CN 201811360899 A CN201811360899 A CN 201811360899A CN 109270009 A CN109270009 A CN 109270009A
- Authority
- CN
- China
- Prior art keywords
- dust
- spectral reflectivity
- plant
- laying
- wave band
- 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/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- 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/55—Specular reflectivity
-
- 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/84—Systems specially adapted for particular applications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/84—Systems specially adapted for particular applications
- G01N2021/8466—Investigation of vegetal material, e.g. leaves, plants, fruits
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm that the present invention relates to a kind of, comprising: obtain plant in the single order inverse of the spectral reflectivity of each wave band and the amount of the laying the dust data of corresponding characterization plant regulation of absorbing dust capability;The single order inverse of the spectral reflectivity of each wave band and the relative coefficient of the corresponding amount of laying the dust data are calculated, and screens acquisition wavelength ranges;Using the corresponding spectral reflectivity single order inverse of wavelength ranges as input, using the corresponding amount of the laying the dust data of wavelength ranges as anticipated output, the inverse model that lays the dust is constructed using grid-search algorithm and Taxonomy and distribution algorithm;Single order inverse to measuring plants in the spectral reflectivity of wavelength ranges is input in the inverse model that lays the dust, estimation obtains lay the dust amount data of the characterization to measuring plants regulation of absorbing dust capability.It is quickly screened by relative coefficient and obtains wavelength ranges, accurately constructed the inverse model that lays the dust in conjunction with grid-search algorithm and Taxonomy and distribution algorithm, realize quick and precisely determining to plant regulation of absorbing dust capability.
Description
Technical field
The present invention relates to Environmental Monitoring and Assessment technical field, more particularly to a kind of based on Taxonomy and distribution algorithm
Plant regulation of absorbing dust capability detection method.
Background technique
With the rapid development of industrialization and urbanization, urban atmospheric pollution is on the rise.Dust is urban atmospheric pollution
In major pollutants, also containing carcinogen and bacterial virus etc., human health is caused in dust in addition to containing heavy metal
Greatly threaten.
Urban greenery plants can effectively in blocks air dust, improve the eco-environmental quality in city.Currently, planting
The dust retention of object has become an important indicator of screening urban greenery plants.Therefore, the regulation of absorbing dust capability of plant how is assessed
Have become the emphasis of current research.
Summary of the invention
Based on this, the plant regulation of absorbing dust capability inspection based on Taxonomy and distribution algorithm that the object of the present invention is to provide a kind of
Survey method has the advantages that quick and precisely obtain the amount of the laying the dust data to measuring plants.
A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm, includes the following steps:
Plant is obtained the single order of the spectral reflectivity of each wave band is reciprocal and corresponding characterization plant regulation of absorbing dust capability
The amount of laying the dust data;
Calculate the single order inverse of the spectral reflectivity of each wave band and the relative coefficient of the corresponding amount of laying the dust data, and root
It is screened according to the relative coefficient and obtains wavelength ranges;Wherein, the absolute value of the relative coefficient is bigger, the corresponding wave
The single order inverse of the spectral reflectivity of section can more characterize the plant amount of laying the dust data;
Using the corresponding spectral reflectivity single order inverse of wavelength ranges as input, with the corresponding amount of the laying the dust data of wavelength ranges
As anticipated output, the inverse model that lays the dust is constructed using grid-search algorithm and Taxonomy and distribution algorithm;
The single order inverse to measuring plants in the spectral reflectivity of wavelength ranges is obtained, and is input in the inverse model that lays the dust,
The acquisition amount of laying the dust data.
The present invention is related to the corresponding amount of laying the dust data by the single order inverse for calculating the spectral reflectivity of each wave band
Property coefficient realizes that quickly screening obtains wavelength ranges, in conjunction with grid-search algorithm and the accurate structure of Taxonomy and distribution algorithm
The inverse model that lays the dust is built, and then realizes quick and precisely determining for the regulation of absorbing dust capability for treating measuring plants.
In one embodiment, the step of the single order inverse of spectral reflectivity of the acquisition plant in each wave band, packet
It includes:
The spectral reflectivity of plant is obtained, and spectral reflectivity is pre-processed;
Spectrum transform is carried out to pretreated spectral reflectivity, the single order for obtaining the spectral reflectivity of each wave band falls
Number.
It is in one embodiment, described that pretreated step is carried out to spectral reflectivity, comprising:
The spectral reflectivity of plant is obtained, and rejects the spectral reflectivity that error is more than the first given threshold, will be remained later
Remaining spectral reflectance data is averaging, using average value as the spectral reflectivity of plant;
Actual spectrum of the nearby principle divided by corresponding blank value, as plant is used to the spectral reflectivity of the plant
Reflectivity;
According to the actual spectrum reflectivity of plant, water vapor absorption wave band is obtained, and rejects the spectrum under water vapor absorption wave band
Reflectivity, the curve of spectrum fitting completion further along spectral reflectivity obtain pretreated spectral reflectivity.
The error of DATA REASONING can be reduced by way of rejecting error information and averaging again;By using nearby principle
Divided by the mode of corresponding blank value, the data obtained under different time sections, different experimental conditions can be made to be comparable, and same
When eliminate experimental situation background brought by error;Data processing by rejecting water vapor absorption wave band, after can facilitating.
In one embodiment, the single order inverse and the corresponding amount of laying the dust number of the spectral reflectivity for calculating each wave band
According to relative coefficient mode are as follows:
Wherein, rijFor the spectral reflectivity first derivative of j-th of wave band in i-th of sample;For the spectrum of j-th of wave band
Reflection
The average value of rate first derivative;xijFor the amount of the laying the dust data of j-th of wave band in i-th of sample;For j-th of wave
The average value of the section amount of laying the dust data;N is the quantity of sample;Each sample includes that the single order of the spectral reflectivity of each wave band is reciprocal
And the corresponding amount of laying the dust data.
It is in one embodiment, described that the step of obtaining wavelength ranges is screened according to relative coefficient, comprising:
The absolute value of the relative coefficient of each wave band is ranked up in the way of from large to small;
The absolute value for screening relative coefficient meets the wave band of the first given threshold as wavelength ranges.
The plant regulation of absorbing dust capability detection device based on Taxonomy and distribution algorithm that the present invention also provides a kind of, comprising:
Data acquisition module, for obtaining plant in the single order inverse of the spectral reflectivity of each wave band and corresponding table
Levy the amount of the laying the dust data of plant regulation of absorbing dust capability;
Wavelength ranges obtain module, for calculate the spectral reflectivity of each wave band single order inverse and the corresponding amount of laying the dust
The relative coefficient of data, and screened according to the relative coefficient and obtain wavelength ranges;Wherein, the relative coefficient is exhausted
Bigger to being worth, the single order inverse of the spectral reflectivity of the corresponding wave band can more characterize the plant amount of laying the dust data;
The inverse model that lays the dust constructs module, is used for using the corresponding spectral reflectivity single order inverse of wavelength ranges as input,
Using the corresponding amount of the laying the dust data of wavelength ranges as anticipated output, using grid-search algorithm and Taxonomy and distribution algorithm
Construct the inverse model that lays the dust;
The amount of laying the dust data acquisition module, it is reciprocal in the single order of the spectral reflectivity of wavelength ranges to measuring plants for obtaining,
And be input in the inverse model that lays the dust, obtain the amount of laying the dust data.
The present invention is related to the corresponding amount of laying the dust data by the single order inverse for calculating the spectral reflectivity of each wave band
Property coefficient realizes that quickly screening obtains wavelength ranges, in conjunction with grid-search algorithm and the accurate structure of Taxonomy and distribution algorithm
The inverse model that lays the dust is built, and then realizes quick and precisely determining for the regulation of absorbing dust capability for treating measuring plants.
The present invention also provides a kind of computer readable storage mediums, store computer program thereon, the computer program
The plant regulation of absorbing dust capability detection based on Taxonomy and distribution algorithm as described in above-mentioned any one is realized when being executed by processor
The step of method.
The present invention also provides a kind of computer equipments comprising reservoir, processor and is stored in the reservoir
And the computer program that can be executed by the processor, the processor are realized when executing the computer program as above-mentioned any
The step of plant regulation of absorbing dust capability detection method described in one based on Taxonomy and distribution algorithm.
In order to better understand and implement, the invention will now be described in detail with reference to the accompanying drawings.
Detailed description of the invention
Fig. 1 is that the present invention is based on the flow charts of the plant regulation of absorbing dust capability detection method of Taxonomy and distribution algorithm;
Fig. 2 is that the present invention obtains plant in the flow chart of the single order inverse of the spectral reflectivity of each wave band;
Fig. 3 is that the present invention carries out pretreated flow chart to spectral reflectivity;
Fig. 4 is to compare figure before and after the present invention rejects water vapor absorption wave band;
Fig. 5 is the flow chart that the present invention obtains wavelength ranges.
Specific embodiment
Referring to Fig. 1, it is the present invention is based on the processes of the plant regulation of absorbing dust capability detection method of Taxonomy and distribution algorithm
Figure.The plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm, includes the following steps:
Step S1: it obtains plant and lays the dust in the single order inverse of the spectral reflectivity of each wave band and corresponding characterization plant
The amount of the laying the dust data of ability.
Wherein, the spectral reflectivity of the plant is plant in the luminous flux of the reflection of each wave band and is incident on object
The ratio between luminous flux.The spectral reflectivity of the specific plant described in the present embodiment be the plant in 350nm-2500nm wave band and
Between be divided into the ratio between the reflected flux of 1nm wave band and the incident flux of the wave band, that is, include reflected flux of the plant in 350nm
It is described with the ratio between the incident flux of the wave band, the plant in the ratio between the reflected flux of 351nm and the incident flux of the wave band ...
Plant is in the ratio between the reflected flux of 2500nm and the incident flux of the wave band.
Step S2: the single order inverse of the spectral reflectivity of each wave band and the correlation system of the corresponding amount of laying the dust data are calculated
Number, and screened according to the relative coefficient and obtain wavelength ranges;Wherein, the absolute value of the relative coefficient is bigger, corresponding
The single order inverse of spectral reflectivity of the wave band more can characterize the plant amount of laying the dust data;
Step S3: corresponding stagnant with wavelength ranges using the corresponding spectral reflectivity single order inverse of wavelength ranges as input
Dust quantity data construct the inverse model that lays the dust as anticipated output, using grid-search algorithm and Taxonomy and distribution algorithm;
Step S5: the single order inverse to measuring plants in the spectral reflectivity of wavelength ranges is obtained, and is input to the inverting that lays the dust
In model, the amount of laying the dust data are obtained.
The present invention is related to the corresponding amount of laying the dust data by the single order inverse for calculating the spectral reflectivity of each wave band
Property coefficient realizes that quickly screening obtains wavelength ranges, in conjunction with grid-search algorithm and the accurate structure of Taxonomy and distribution algorithm
The inverse model that lays the dust is built, and then realizes quick and precisely determining for the regulation of absorbing dust capability for treating measuring plants.
Referring to Fig. 2, it obtains plant in the flow chart of the single order inverse of the spectral reflectivity of each wave band for the present invention.
In one embodiment, in step S1, the single order of spectral reflectivity of the acquisition plant in each wave band is reciprocal
The step of, comprising:
Step S11: the spectral reflectivity of plant is obtained, and spectral reflectivity is pre-processed.
Step S12: spectrum transform is carried out to pretreated spectral reflectivity, obtains the spectral reflectivity of each wave band
Single order is reciprocal.
Wherein, by carrying out spectrum transform to pretreated spectral reflectivity, compressible ambient noise is to target information
Influence, improve follow-up data processing accuracy.
Referring to Fig. 3, it carries out pretreated flow chart to spectral reflectivity for the present invention.
In one embodiment, pretreated step is carried out to spectral reflectivity described in step S11, comprising:
Step S111: obtaining the spectral reflectivity of plant, and rejects the spectral reflectivity that error is more than the first given threshold,
Remaining spectral reflectance data is averaging later, using average value as the spectral reflectivity of plant.
In one embodiment, the blade in the multiple orientation of single plant is chosen for a kind of plant, every leaf is chosen solid
Fixed certain point is as measurement point, 5 groups of spectral reflectivities of each measurement point duplicate measurements, then rejects in this 5 groups of spectral reflectivities
It is averaging after the value in the first given threshold, then by remaining spectral reflectivity, the light of this piece leaf is characterized with average value
Reflectivity is composed, and then reduces the error of DATA REASONING, improves the accuracy of DATA REASONING.Multiple orientation therein can plant to be single
The orientation such as the upper and lower, left and right of strain.First given threshold is 5 groups of spectral reflectivities according to measurement, selected by experience
For rejects with other data with large error data threshold value.
In one embodiment, the spectral reflectivity of plant is obtained using the measurement of 3 spectrometer of ASD Field-Spec, by
It is .asp format in the data format of measurement, the data of this format read and handle, therefore, measuring with being unable to direct convenience
Before carrying out data processing after the spectral reflectivity of acquisition plant, View Spectral Pro software spectral reflectivity is also utilized
Data Format Transform be .txt format.
Step S112: use nearby principle divided by corresponding blank value the spectral reflectivity of the plant, as plant
Actual spectrum reflectivity.
To be comparable the data obtained under different time sections, different experimental conditions, and experimental situation is eliminated simultaneously
Error brought by background, section obtains different blank values to the application in different times, hence in some period
The spectral reflectivity of the plant of acquisition then needs for the spectral reflectivity of the plant to be in the blank value of the period, as place
The spectral reflectivity of plant after reason.
Step S113: according to the actual spectrum reflectivity of plant, water vapor absorption wave band is obtained, and rejects water vapor absorption wave band
Under spectral reflectivity, the curve of spectrum fitting completion further along spectral reflectivity obtains pretreated spectral reflectivity.
Due to determining instrument of spectral reflectivity itself, the spectral reflectivity determined can be due to the absorption of steam
It is affected to the spectral reflectivity of plant, and the wavelength band studying also without too big meaning, therefore to vegetation spectrum
Data processing by rejecting water vapor absorption wave band, after can facilitating.Referring to Fig. 4, it rejects water vapor absorption wave for the present invention
Compare figure in the front and back of section.Specifically, the actual spectrum reflectivity data of plant is fitted to spectral reflectivity curve, pass through observation
Are there is abnormal wave band as water vapor absorption wave band by spectral reflectivity curve in data in spectral reflectivity curve;Again by the water
Spectral reflectivity in vapour absorption bands is rejected, and the curve of spectrum fitting completion acquisition further along spectral reflectivity is pretreated
Spectral reflectivity.
In one embodiment, spectrum transform is carried out to pretreated spectral reflectivity described in step S22, obtained each
The calculation of the single order inverse of the spectral reflectivity of a wave band are as follows:
Wherein, λi+1, λi, λi-1For adjacent wavelength, dR (λi) it is wavelength XiFirst derivative spectrum, R (λi+1), R (λi), R
(λi-1) be respectively wavelength be λi+1, λi, λi-1The reflectivity at place.
In one embodiment, the calculation of the amount of the laying the dust data of the plant regulation of absorbing dust capability of acquisition characterization described in step S1
Are as follows:
D=Δ W/S
In above-mentioned formula, D is the amount of the laying the dust data for characterizing plant regulation of absorbing dust capability;Δ W is a period of time implants blade
The amount of laying the dust, S are blade area.
Wherein, the amount of the laying the dust data of each blade of plant are consistent, i.e., the spectrum of each wave band is anti-in each blade
The single order corresponding amount of the laying the dust data reciprocal for penetrating rate are identical, it can think the corresponding sample of each blade, each sample
The single order inverse of the spectral reflectivity of each wave band including the blade and the corresponding amount of laying the dust data.
In one embodiment, the single order inverse of the spectral reflectivity of each wave band of calculating described in step 2 and corresponding
The mode of the relative coefficient of the amount of laying the dust data are as follows:
Wherein, rijFor the spectral reflectivity first derivative of j-th of wave band in i-th of sample;For the spectrum of j-th of wave band
The average value of reflectivity first derivative;xijFor the amount of the laying the dust data of j-th of wave band in i-th of sample;It is stagnant for j-th of wave band
The average value of dust quantity data;N is the quantity of sample;Each sample include the single order of the spectral reflectivity of each wave band it is reciprocal and
The corresponding amount of laying the dust data.
Above-mentioned formula can measure the spectral reflectivity of some wave band single order inverse and the corresponding amount of laying the dust data the two
The linear index of variable, value range are [- 1,1].And the absolute value of the calculated relative coefficient of above-mentioned formula is bigger,
The correlation of the single order inverse and the amount of laying the dust data that indicate the spectral reflectivity of the wave band is stronger, i.e., the wave band can more characterize plant
The amount of laying the dust data, the correlation intensity of two variables can be judged by following value range, referring specifically to table 1- correlation
Coefficient table corresponding with correlation intensity.
Table 2- relative coefficient table corresponding with correlation intensity
Relative coefficient (absolute value) | Strength of correlation |
0.8-1.0 | Extremely strong correlation |
0.6-0.8 | Strong correlation |
0.4-0.6 | Moderate correlation |
0.2-0.4 | Weak correlation |
0.0-0.2 | It is extremely weak related or without correlation |
Referring to Fig. 5, its flow chart for obtaining wavelength ranges for the present invention.
In one embodiment, the step of obtaining wavelength ranges, is screened according to relative coefficient described in step 2, comprising:
Step 21: the absolute value of the relative coefficient of each wave band is ranked up in the way of from large to small;
Step 22: the absolute value for screening relative coefficient meets the wave band of the first given threshold as wavelength ranges.
First given threshold can be the specific value about relative coefficient, and such as the first given threshold is correlation system
Several absolute values is greater than 0.8, i.e. wave band of the absolute value of relative coefficient greater than 0.8 is as wavelength ranges;First setting
Threshold value may be ranking value, and such as the first given threshold is that the absolute value of relative coefficient sorts within preceding 10, i.e. correlation
The sequence of the absolute value of coefficient from large to small is correlation using corresponding wave band as wavelength ranges or the first given threshold 10
The absolute value sequence quantity of property coefficient accounts for the 10% of all relative coefficient quantity using corresponding wave band as wavelength ranges etc..
In one embodiment, Taxonomy and distribution (Classification and regression tree, CART)
Using two points of recursive subdivision technologies, the dividing method of the gini index estimation function based on minimum range, by current sample set
It is divided into 2 sub- sample sets, the non-leaf nodes of each generation there are 2 branches.The decision tree that CART algorithm generates is a letter
Clean binary tree is classification tree (classification tree), if target is continuous if target is discrete variable
Variable is then regression tree (regression tree).
Wherein, the setting of the parameter of Taxonomy and distribution algorithm is most important to the training precision of algorithm, based on classification with
The main parameter setting range of parameter for returning tree algorithm, is found most using grid-search algorithms (grid-search algorithm)
Excellent parameter.Grid-search algorithm uses the means of exhaustive search, i.e., in the parameter selection of all candidates, by following
Ring traversal, attempts each possibility, the parameter to behave oneself best is exactly final as a result, its principle is like to look for most in array
Big value.
The thought of Taxonomy and distribution algorithm is so that the variance in group is reached minimum, and corresponding between-group variance reaches most
Greatly, even if the left and right branch differentiation of tree division is maximum.It is described to be constructed using grid-search algorithm and Taxonomy and distribution algorithm
Lay the dust inverse model the step of, comprising: passing through impurity level index, to calculate the corresponding spectral reflectivity single order of wavelength ranges reciprocal
Various cutting/combined situations find out optimum value combination/cut-off of the corresponding spectral reflectivity single order inverse of wavelength ranges;So
Optimum value combination/cut-off of the corresponding spectral reflectivity single order inverse of more each wavelength ranges afterwards;It finally finds out best anti-
Drill the optimum value group of the corresponding spectral reflectivity single order inverse of wave band spectral reflectivity single order inverse corresponding with the wavelength ranges
Conjunction/cut-off.
The growth entirely set in Taxonomy and distribution algorithm is a recursive procedure, until termination condition;The termination item
Part includes: that node is pure node, i.e., the target variable value of all records is identical;The depth of tree has reached preassigned maximum
Value;The maximum drop-out value of hybrid UV curing is less than preassigned value;The recorded amounts of node is recorded less than preassigned minimum node
Amount;The predictive variable value of all records is identical in one node.
It is constructed in the inverse model that lays the dust using grid-search algorithm and Taxonomy and distribution algorithm, each individual node
Minimum samples be sequentially 1,2,3,4,5,6,7,8,9,10.When determining parameter by grid-search algorithm, specifically
Step are as follows: import sample, sample is divided into training set and test set;It arranges in pairs or groups alternative parameter, and certain candidate is set separately
Value;The methods of marking of defined parameters;GridSearchCV is called, by SVC (), tuned_parameters, cv=5, also
Scoring transmitting is entered, and with training set training this learner clf, recalling clf.best_params_ can be directly obtained
Best parameter collocation result.
CART algorithm computation complexity is not high, and regular gradual change is it can be readily appreciated that very steady to intermediate missing values.
The result shows that constructing quantitative estimation model precision with higher, models fitting using Taxonomy and distribution algorithm
It mutually shuts, the R of obtained model2Greater than 0.9.
For the precision for verifying inverse model, the actual measured value of the amount of laying the dust data is compared with predicted value.Model
Precision test is measured using quantitative accuracy evaluation index, specially determines coefficient (the coefficient of
Determination, R2), root-mean-square error (the root-mean-square error, RMSE), mean absolute error
(the mean absolute error, MAE), average deviation (Bias), calculation formula is as follows.
Wherein, xiFor the actual measurement amount of laying the dust of i-th of sample, xi' for i-th of sample model prediction the amount of laying the dust,For institute
There is the average value of the actual measurement amount of laying the dust,For the average value of all model prediction amounts of laying the dust, N is the quantity of sample.
R2Indicate the linearly related degree between two data, MAE and RMSE are used to the integral level of evaluated error, and Bias is anti-
Estimation data and the departure degree of measured data numerically are reflected.
The result shows that the verifying precision RMSE of quantitative estimation model is respectively less than 0.36g/m2, efficiency of inverse process is good.
In one embodiment, after the step of amount of laying the dust data are obtained described in step S4, further includes: obtain multiple groups and wait for
The amount of the laying the dust data of measuring plants, and corresponding plant is planted in setting regions according to the amount of the laying the dust data.Specifically, according to stagnant
Dust quantity data are bigger, illustrate that the regulation of absorbing dust capability to measuring plants is stronger, and such as there is the area of sandstorm can plant according to setting regions
Regulation of absorbing dust capability is stronger to measuring plants.
The plant regulation of absorbing dust capability detection device based on Taxonomy and distribution algorithm that the present invention also provides a kind of, comprising:
Data acquisition module, for obtaining plant in the single order inverse of the spectral reflectivity of each wave band and corresponding table
Levy the amount of the laying the dust data of plant regulation of absorbing dust capability;
Wavelength ranges obtain module, for calculate the spectral reflectivity of each wave band single order inverse and the corresponding amount of laying the dust
The relative coefficient of data, and screened according to the relative coefficient and obtain wavelength ranges;Wherein, the relative coefficient is exhausted
Bigger to being worth, the single order inverse of the spectral reflectivity of the corresponding wave band can more characterize the plant amount of laying the dust data;
The inverse model that lays the dust constructs module, is used for using the corresponding spectral reflectivity single order inverse of wavelength ranges as input,
Using the corresponding amount of the laying the dust data of wavelength ranges as anticipated output, using grid-search algorithm and Taxonomy and distribution algorithm
Construct the inverse model that lays the dust;
The amount of laying the dust data acquisition module, it is reciprocal in the single order of the spectral reflectivity of wavelength ranges to measuring plants for obtaining,
And be input in the inverse model that lays the dust, obtain the amount of laying the dust data.
The present invention is related to the corresponding amount of laying the dust data by the single order inverse for calculating the spectral reflectivity of each wave band
Property coefficient realizes that quickly screening obtains wavelength ranges, in conjunction with grid-search algorithm and the accurate structure of Taxonomy and distribution algorithm
The inverse model that lays the dust is built, and then realizes quick and precisely determining for the regulation of absorbing dust capability for treating measuring plants.
In one embodiment, the data acquisition module includes:
Preprocessing module is pre-processed for obtaining the spectral reflectivity of plant, and to spectral reflectivity.
Spectrum transform module obtains the light of each wave band for carrying out spectrum transform to pretreated spectral reflectivity
The single order for composing reflectivity is reciprocal.
Wherein, the preprocessing module includes:
Mean value calculation module, for obtaining the spectral reflectivity of plant, and rejecting error is more than the first given threshold
Remaining spectral reflectance data is averaging by spectral reflectivity later, using average value as the spectral reflectivity of plant.
Actual spectrum reflectivity computing module, for the spectral reflectivity to the plant using nearby principle divided by correspondence
Blank value, the actual spectrum reflectivity as plant.
Water vapor absorption wave band rejects module, for the actual spectrum reflectivity according to plant, obtains water vapor absorption wave band, and
Reject water vapor absorption wave band under spectral reflectivity, further along spectral reflectivity the curve of spectrum fitting completion pre-processed after
Spectral reflectivity.
In one embodiment, described that spectrum transform is carried out to pretreated spectral reflectivity, obtain each wave band
The calculation of the single order inverse of spectral reflectivity are as follows:
Wherein, λi+1, λi, λi-1For adjacent wavelength, dR (λi) it is wavelength XiFirst derivative spectrum, R (λi+1), R (λi), R
(λi-1) be respectively wavelength be λi+1, λi, λi-1The reflectivity at place.
In one embodiment, the calculation of the amount of the laying the dust data for obtaining characterization plant regulation of absorbing dust capability are as follows:
D=Δ W/S
In above-mentioned formula, D is the amount of the laying the dust data for characterizing plant regulation of absorbing dust capability;Δ W is a period of time implants blade
The amount of laying the dust, S are blade area.
Wherein, the amount of the laying the dust data of each blade of plant are consistent, i.e., the spectrum of each wave band is anti-in each blade
The single order corresponding amount of the laying the dust data reciprocal for penetrating rate are identical, it can think the corresponding sample of each blade, each sample
The single order inverse of the spectral reflectivity of each wave band including the blade and the corresponding amount of laying the dust data.
In one embodiment, the single order inverse and the corresponding amount of laying the dust number of the spectral reflectivity for calculating each wave band
According to relative coefficient mode are as follows:
Wherein, rijFor the spectral reflectivity first derivative of j-th of wave band in i-th of sample;For the spectrum of j-th of wave band
The average value of reflectivity first derivative;xijFor the amount of the laying the dust data of j-th of wave band in i-th of sample;It is stagnant for j-th of wave band
The average value of dust quantity data;N is the quantity of sample.
In one embodiment, the wavelength ranges obtain the step of module, comprising:
Sorting module, for arranging the absolute value of the relative coefficient of each wave band in the way of from large to small
Sequence;
Screening module, the absolute value for screening relative coefficient meet the wave band of the first given threshold as inverting wave
Section.
First given threshold can be the specific value about relative coefficient, and such as the first given threshold is correlation system
Several absolute values is greater than 0.8, i.e. wave band of the absolute value of relative coefficient greater than 0.8 is as wavelength ranges;First setting
Threshold value may be ranking value, and such as the first given threshold is that the absolute value of relative coefficient sorts within preceding 10, i.e. correlation
The sequence of the absolute value of coefficient from large to small is correlation using corresponding wave band as wavelength ranges or the first given threshold 10
The absolute value sequence quantity of property coefficient accounts for the 10% of all relative coefficient quantity using corresponding wave band as wavelength ranges etc..
The present invention also provides a kind of computer readable storage mediums, store computer program thereon, the computer program
The plant regulation of absorbing dust capability detection based on Taxonomy and distribution algorithm as described in above-mentioned any one is realized when being executed by processor
The step of method.
It wherein includes storage medium (the including but not limited to disk of program code that the present invention, which can be used in one or more,
Memory, CD-ROM, optical memory etc.) on the form of computer program product implemented.Computer-readable storage media packet
Permanent and non-permanent, removable and non-removable media is included, can be accomplished by any method or technique information storage.Letter
Breath can be computer readable instructions, data structure, the module of program or other data.The example packet of the storage medium of computer
Include but be not limited to: phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM),
Other kinds of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory
(EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM), digital versatile disc
(DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices or any other non-biography
Defeated medium, can be used for storage can be accessed by a computing device information.
The present invention also provides a kind of computer equipment, including reservoir, processor and it is stored in the reservoir simultaneously
The computer program that can be executed by the processor, the processor are realized when executing the computer program as above-mentioned any one
The step of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm described in item.
The embodiments described above only express several embodiments of the present invention, and the description thereof is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
It says, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to protection of the invention
Range.
Claims (10)
1. a kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm, which comprises the steps of:
Plant is obtained in the single order inverse of the spectral reflectivity of each wave band and laying the dust for corresponding characterization plant regulation of absorbing dust capability
Measure data;
The single order inverse of the spectral reflectivity of each wave band and the relative coefficient of the corresponding amount of laying the dust data are calculated, and according to institute
It states relative coefficient screening and obtains wavelength ranges;Wherein, the absolute value of the relative coefficient is bigger, the corresponding wave band
The single order inverse of spectral reflectivity can more characterize the plant amount of laying the dust data;
Using the corresponding spectral reflectivity single order inverse of wavelength ranges as input, using the corresponding amount of the laying the dust data of wavelength ranges as
Anticipated output constructs the inverse model that lays the dust using grid-search algorithm and Taxonomy and distribution algorithm;
The single order inverse to measuring plants in the spectral reflectivity of wavelength ranges is obtained, and is input in the inverse model that lays the dust, is obtained
The amount of laying the dust data.
2. the plant regulation of absorbing dust capability detection method according to claim 1 based on Taxonomy and distribution algorithm, feature exist
In the step of the single order inverse of spectral reflectivity of the acquisition plant in each wave band, comprising:
The spectral reflectivity of plant is obtained, and spectral reflectivity is pre-processed;
Spectrum transform is carried out to pretreated spectral reflectivity, the single order for obtaining the spectral reflectivity of each wave band is reciprocal.
3. the plant regulation of absorbing dust capability detection method according to claim 2 based on Taxonomy and distribution algorithm, feature exist
In described to carry out pretreated step to spectral reflectivity, comprising:
The spectral reflectivity of plant is obtained, and rejects the spectral reflectivity that error is more than the first given threshold, it later will be remaining
Spectral reflectance data is averaging, using average value as the spectral reflectivity of plant;
Use nearby principle divided by corresponding blank value the spectral reflectivity of the plant, the actual spectrum as plant reflects
Rate;
According to the actual spectrum reflectivity of plant, water vapor absorption wave band is determined, and it is anti-to reject the corresponding spectrum of water vapor absorption wave band
Rate is penetrated, the curve of spectrum fitting completion further along spectral reflectivity obtains pretreated spectral reflectivity.
4. the plant regulation of absorbing dust capability detection method according to claim 1 based on Taxonomy and distribution algorithm, feature exist
In the mode of the single order inverse of the spectral reflectivity for obtaining each wave band are as follows:
Wherein, λi+1, λi, λi-1For adjacent wavelength, dR (λi) it is wavelength XiFirst derivative spectrum, R (λi+1), R (λi), R (λi-1)
Be respectively wavelength be λi+1, λi, λi-1The reflectivity at place.
The calculation of the amount of the laying the dust data for obtaining characterization plant regulation of absorbing dust capability are as follows:
D=Δ W/S
Wherein, D is the amount of the laying the dust data for characterizing plant regulation of absorbing dust capability;Δ W is the amount of laying the dust of a period of time implants blade;S is
Blade area.
5. the plant regulation of absorbing dust capability detection method according to claim 1 based on Taxonomy and distribution algorithm, feature exist
In the mode of the relative coefficient of the single order inverse and corresponding amount of laying the dust data of the spectral reflectivity for calculating each wave band
Are as follows:
Wherein, rijFor the spectral reflectivity first derivative of j-th of wave band in i-th of sample;For the spectral reflectance of j-th of wave band
The average value of rate first derivative;xijFor the amount of the laying the dust data of j-th of wave band in i-th of sample;For j-th of wave band amount of laying the dust
The average value of data;N is the quantity of sample;Each sample includes that the single order of the spectral reflectivity of each wave band is reciprocal and corresponding
The amount of laying the dust data.
6. the plant regulation of absorbing dust capability detection method according to claim 1 based on Taxonomy and distribution algorithm, feature exist
In described to screen the step of obtaining wavelength ranges according to relative coefficient, comprising:
The absolute value of the relative coefficient of each wave band is ranked up in the way of from large to small;
The absolute value for screening relative coefficient meets the wave band of the first given threshold as wavelength ranges.
7. the plant regulation of absorbing dust capability detection method according to claim 1 based on Taxonomy and distribution algorithm, feature exist
After in, the acquisition amount of the laying the dust data the step of, further includes: obtain the amount of the laying the dust data that multiple groups wait for measuring plants, and according to institute
It states the amount of laying the dust data and plants corresponding plant in setting regions.
8. a kind of plant regulation of absorbing dust capability detection device based on Taxonomy and distribution algorithm characterized by comprising
Data acquisition module is planted for obtaining plant in the single order inverse of the spectral reflectivity of each wave band and corresponding characterization
The amount of the laying the dust data of object regulation of absorbing dust capability;
Wavelength ranges obtain module, for calculate the spectral reflectivity of each wave band single order inverse and the corresponding amount of laying the dust data
Relative coefficient, and according to the relative coefficient screen obtain wavelength ranges;Wherein, the absolute value of the relative coefficient
Bigger, the single order inverse of the spectral reflectivity of the corresponding wave band can more characterize the plant amount of laying the dust data;
The inverse model that lays the dust constructs module, is used for using the corresponding spectral reflectivity single order inverse of wavelength ranges as input, with anti-
The corresponding amount of the laying the dust data of wave band are drilled as anticipated output, are constructed using grid-search algorithm and Taxonomy and distribution algorithm
Lay the dust inverse model;
The amount of laying the dust data acquisition module, it is reciprocal and defeated in the single order of the spectral reflectivity of wavelength ranges to measuring plants for obtaining
Enter into the inverse model that lays the dust, obtains the amount of laying the dust data.
9. a kind of computer readable storage medium, stores computer program thereon, which is characterized in that the computer program is located
Reason device realizes that plant as claimed in any of claims 1 to 7 in one of claims based on Taxonomy and distribution algorithm lays the dust energy when executing
The step of power detection method.
10. a kind of computer equipment, which is characterized in that including reservoir, processor and be stored in the reservoir and can
The computer program executed by the processor, the processor realize such as claim 1 to 7 when executing the computer program
Any one of described in the plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm the step of.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811360899.2A CN109270009A (en) | 2018-11-15 | 2018-11-15 | A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811360899.2A CN109270009A (en) | 2018-11-15 | 2018-11-15 | A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109270009A true CN109270009A (en) | 2019-01-25 |
Family
ID=65190226
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811360899.2A Pending CN109270009A (en) | 2018-11-15 | 2018-11-15 | A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109270009A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110849334A (en) * | 2019-09-30 | 2020-02-28 | 广州地理研究所 | Island reef shallow sea water depth prediction method based on classification and regression tree |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697964A (en) * | 2015-03-10 | 2015-06-10 | 西北大学 | Method for quantitative analysis of sulfur and phosphorus in steel and iron by combining random forest algorithm with laser induced breakdown spectroscopy |
CN104897592A (en) * | 2015-06-11 | 2015-09-09 | 石河子大学 | Monitoring method of salt ion content in saline soil based on hyperspectral technology |
-
2018
- 2018-11-15 CN CN201811360899.2A patent/CN109270009A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697964A (en) * | 2015-03-10 | 2015-06-10 | 西北大学 | Method for quantitative analysis of sulfur and phosphorus in steel and iron by combining random forest algorithm with laser induced breakdown spectroscopy |
CN104897592A (en) * | 2015-06-11 | 2015-09-09 | 石河子大学 | Monitoring method of salt ion content in saline soil based on hyperspectral technology |
Non-Patent Citations (3)
Title |
---|
CHONGYANG WANG等: "ESTIMATION MODEL FOR DUST-RETENTION CONTENT OF MAIN GREEN PLANTS IN SOUTH CHINA BASED ON THE RED EDGE OF REFLECTANCE", 《IEEE XPLORE》 * |
WENLONG JING等: "Machine Learning for Estimating Leaf Dust Retention Based on Hyperspectral Measurements", 《HINDAWI JOURNAL OF SENSORS》 * |
肖慧玲等: "园林植物滞尘能力分析及其叶面滞尘高光谱定量遥感模型", 《长江流域资源与环境》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110849334A (en) * | 2019-09-30 | 2020-02-28 | 广州地理研究所 | Island reef shallow sea water depth prediction method based on classification and regression tree |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jumin et al. | Solar radiation prediction using boosted decision tree regression model: A case study in Malaysia | |
CN106815643B (en) | Infrared spectroscopy Model Transfer method based on random forest transfer learning | |
Calders et al. | Nondestructive estimates of above‐ground biomass using terrestrial laser scanning | |
CN112464920B (en) | FY-3D infrared hyperspectral cloud detection method based on extreme random tree | |
Dye et al. | Examining the utility of random forest and AISA Eagle hyperspectral image data to predict Pinus patula age in KwaZulu-Natal, South Africa | |
CN116205310B (en) | Soil water content influence factor sensitive interval judging method based on interpretable integrated learning model | |
CN110162872B (en) | Biomass estimation method fusing sample plot data and forest resource clearing data | |
CN110569605B (en) | NSGA 2-ELM-based non-glutinous rice leaf nitrogen content inversion model method | |
CN106952170A (en) | A kind of county domain Natural quality of cultivated land evaluation method based on random forest | |
Richter et al. | Goodness-of-fit measures: what do they tell about vegetation variable retrieval performance from Earth observation data | |
Zervas et al. | HeLM: A macrophyte-based method for monitoring and assessment of Greek lakes | |
CN103278467A (en) | Rapid nondestructive high-accuracy method with for identifying abundance degree of nitrogen element in plant leaf | |
CN109521437A (en) | Multiplex physicotherapy laser radar wavelength selection method towards the detection of vegetation biochemical parameter | |
CN109270012A (en) | A kind of plant regulation of absorbing dust capability detection method based on relative coefficient | |
Evans et al. | A comparison between data requirements and availability for calibrating predictive ecological models for lowland UK woodlands: learning new tricks from old trees | |
Gontijo et al. | Electricity price forecasting on electricity spot market: A case study based on the Brazilian Difference Settlement Price | |
CN109270009A (en) | A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm | |
CN113139717A (en) | Crop seedling condition grading remote sensing monitoring method and device | |
Perraut et al. | Benchmarking the fundamental parameters of Ap stars with optical long-baseline interferometric measurements | |
CN109270010A (en) | A kind of plant regulation of absorbing dust capability detection method based on random forests algorithm | |
Li et al. | Use of terrestrial laser scanning to obtain the stem diameters of Larix olgensis and construct compatible taper-volume equations | |
CN109270011A (en) | A kind of plant regulation of absorbing dust capability detection method based on machine learning algorithm | |
CN107977727A (en) | A kind of method that probability is blocked based on social development and climatic factor prediction cable network | |
Hachicha et al. | Prediction of plant growth based on statistical methods and remote sensing data | |
Beyene | Estimation of forest variable and aboveground biomass using terrestrial laser scanning in the tropical rainforest |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190125 |