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 PDF

Info

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
Application number
CN201811360899.2A
Other languages
Chinese (zh)
Inventor
荆文龙
周成虎
王重洋
韩留生
张晨
周霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou Institute of Geography of GDAS
Original Assignee
Guangzhou Institute of Geography of GDAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangzhou Institute of Geography of GDAS filed Critical Guangzhou Institute of Geography of GDAS
Priority to CN201811360899.2A priority Critical patent/CN109270009A/en
Publication of CN109270009A publication Critical patent/CN109270009A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/55Specular reflectivity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N2021/8466Investigation 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

A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm
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.
CN201811360899.2A 2018-11-15 2018-11-15 A kind of plant regulation of absorbing dust capability detection method based on Taxonomy and distribution algorithm Pending CN109270009A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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