The content of the invention
The technical problem to be solved in the present invention, is to provide a kind of forestry pests & diseases based on multi-source image information intelligently to know
Other system, relative satellite remote sensing, airborne and spaceborne RS cost are lower, and accuracy is high, easily controllable.
What the present invention was realized in:A kind of forestry pests & diseases intelligent identifying system based on multi-source image information, it is described
System includes unmanned plane and earth station, the equipment of UAV flight have humiture module, communication module, control module and with control
The connected d GPS locating module of molding block, sensing module and image capture module, the humiture module are connected with control module;
The communication module, realizes the communication of the control module and earth station of unmanned plane;By satellite remote-sensing image, the remote sensing of unmanned plane
Image and image recognition technology to forestry pests & diseases information extraction, using spectral information and phase in time series image data
Information, image automatic identification acquisition pest and disease damage is carried out to the image that image capture module is gathered, and a situation arises.
Further, the sensing module includes many ripples of the pulse radar of monitoring forestry pests & diseases and the monitoring forest reserves
Section spectrum scanner, infrared spectrometer.
Further, it is described sick to forestry by satellite remote-sensing image, the remote sensing image of unmanned plane and image recognition technology
The mode that insect pest information is extracted specifically includes following steps:
Step 1, download satellite remote-sensing image, and to obtain satellite remote-sensing image pre-process;
Step 2, the coverage for extracting forestry pests & diseases;
The light of Mono temporal and multi-temporal vegetation indices in step 3, the time series image data of extraction satellite remote-sensing image
Spectrum signature;
Step 4, disease monitor the phase, synchronized in image capturing unmanned plane low latitude investigation;
Forest reserves image and state that step 5, unmanned plane collection draw the line interior;
Step 6, the image of control module analysis image capture module collection and remote sensing images and light from sensing module
Modal data, the pulse radar of sensing module whether there is insect, the infrared spectrometer of sensing module according to the electromagnetic wave detection being reflected back
With the spectral signature of the different-waveband of the multiband spectrum scanner collection forest reserves, and collection of illustrative plates is generated;
Step 7, pass the data in step 6 back earth station and further analyze;With reference to unmanned plane low latitude survey data, utilize
Image recognition technology, screens the spectral signature of disease monitoring;
It is step 9, different according to forest reserves growth period and whether there is pest and disease damage to be respectively established;
Step 10, a situation arises for the pest and disease damage based on model, spectral information divergence and atlas analysis forestry, the spectrum letter
Breath includes:The spectral signature of the spectral signature of vegetation index, the spectral signature of different-waveband and disease monitoring.
Further, the step 3 is specially:Employ alternative features of 13 spectral signatures as disease monitoring, 13
Individual spectral signature includes the blue RB compatible with the multispectral satellite image of most middle high-resolutions, green RG, red RR, near-infrared RNIR
Passage primary reflection rate, and normalization difference vegetation index NDVI, ratio vegetation index SR, green normalization difference vegetation
Index GNDVI, adjustment soil lightness vegetation index SAVI, triangle vegetation index TVI, improve red side ratio vegetation index MSR, non-
Linear vegetation index NLI, renormalization vegetation index RDVI, soil regulation vegetation index nine broadband vegetation indexs of OSAVI;
The vegetation index that Mono temporal and multidate two kinds of version is respectively adopted is analyzed, wherein, Mono temporal vegetation index by certain for the moment
Phase image wave band reflectivity is calculated, for reflecting Physiology and biochemistry state of the vegetation on certain time point;Multidate vegetation
Index is normalized according to the Mono temporal vegetation index of two phases and is calculated, and becomes for reflecting that pest and disease damage develops in woodland
The characteristics of change.
Further, the scope of the image wave band includes visible ray and near infrared band.
Further, the satellite remote-sensing image of described pair of acquisition carries out pretreatment includes radiation calibration, atmospheric correction, geometry
Correction and cloud removal.
Further, the coverage of the forestry pests & diseases is according to forest land classification vector figure or multi_temporal images enter
Row classification is obtained.
Further, the multi_temporal images need to combine in assorting process forest land use pattern data, terrain data and
Phenology experience, takes decision tree, maximum likelihood classification or neutral net to carry out planting range extraction.
Further, the Spectra feature extraction mode of the disease monitoring is:According to the emphasis area that unmanned plane low latitude is investigated
A situation arises for domain disease, and key area is divided into normal sample and sample two parts of catching an illness;
The Mono temporal and multidate for extracting two class sample point multi-form spectral signatures from image identification system respectively are special
Value indicative;
Mono temporal or multidate version to every kind of spectral signature, using independent samples t test is relatively normal and sample of catching an illness
This difference degree;
The difference degree of a certain feature is characterized using the preset value p of t inspections, and accordingly generates an all kinds of multi-forms
Spectral signature different phases and when it is combined in p value statistical table, wherein, p value is smaller, normal and sample of catching an illness difference
Bigger, response of the feature to defect information is stronger.
The invention has the advantages that:Satellite image of the present invention based on multidate, the remote sensing image of unmanned plane, image are known
The forestry pests & diseases intelligent identification technology of other technology, make full use of the spectral information in time series image data and when believe
Breath, spectral information divergence analysis is introduced forestry by combining geographic information system GIS, global position system GPS, remote sensing technology
Pest and disease damage generation area, the system that proposition is monitored on a large scale using the star in certain area-ground synchrodata to disease, has
Effect reduces the cost of the field work of disease monitoring, and the extension by point and face has been carried out to traditional disease monitoring mode, is easy to
Government department and Management offorestry department in time, accurately grasp and understand the important informations such as the generation of region disease and the order of severity.
Specific embodiment
Refer to shown in Fig. 1 to Fig. 3, a kind of forestry pests & diseases intelligent identifying system based on multi-source image information, it is special
Levy and be:The system includes unmanned plane and earth station, and the equipment of UAV flight has humiture module, communication module, control
Module and the d GPS locating module, sensing module and the image capture module that are connected with control module, the humiture module and control
Molding block is connected;The communication module, realizes the communication of the control module and earth station of unmanned plane;By satellite remote-sensing image,
The remote sensing image and image recognition technology of unmanned plane to forestry pests & diseases information extraction, using the light in time series image data
Spectrum information and when phase information, to image capture module collection image carry out image automatic identification obtain pest and disease damage a situation arises.
The sensing module includes the multiband spectrum scanner, infrared of the pulse radar of monitoring forestry pests & diseases and the monitoring forest reserves
Spectrometer.Humiture module is used to gather the humiture situation of forestry, and d GPS locating module is used to position unmanned plane in forestry
Position.
Wherein, it is described by satellite remote-sensing image, the remote sensing image of unmanned plane and image recognition technology to forestry pests & diseases
The mode of information extraction specifically includes following steps:
Step 1, order download satellite remote-sensing image, and to obtain satellite remote-sensing image pre-process;
Step 2, the coverage for extracting forestry pests & diseases;
The light of Mono temporal and multi-temporal vegetation indices in step 3, the time series image data of extraction satellite remote-sensing image
Spectrum signature;Step 3 extracting mode:On the basis of comprehensive survey is suitable to the documents and materials of crop disease monitoring both at home and abroad, respectively
Alternative features of 13 spectral signatures as disease monitoring are employed, including it is simultaneous with the multispectral satellite image of most middle high-resolutions
Blue (RB), green (RG), red (RR), the passage primary reflection rate of near-infrared (RNIR) held, and (difference vegetation refers to for normalization
Number) NDVI, (ratio vegetation index) SR, (green vegetation index) GNDVI, (adjustment soil lightness vegetation index)
SAVI, (triangle vegetation index) TVI, (improve red side ratio vegetation index) MSR, (non-linear vegetation index) NLI, (return one
Change vegetation index) RDVI, (soil regulation vegetation index) nine broadband vegetation indexs of OSAVI;(each index fixed is shown in Table 1).
The present invention is respectively adopted Mono temporal and the vegetation index of multidate two kinds of version is analyzed.Wherein, Mono temporal vegetation index by
A certain phase image wave band reflectivity is calculated, for reflecting Physiology and biochemistry state of the vegetation on certain time point;When many
Phase vegetation index is normalized according to the Mono temporal vegetation index of two phases and is calculated, for reflecting pest and disease damage in woodland
The characteristics of development and change.Assuming that there is M and N to carry out two phases for before and after, for a certain vegetation index form V1, multidate reflectivity
Computing formula is (VIN-VIM)/(VIN+VIM).Wherein, the period of VIN vegetation indexs ordinal number N, VIM vegetation index ordinal numbers M's
Phase, for n phase, is carried out combination of two by the period, and C2n multidate spectral signature is can obtain altogether.It is above-mentioned in this method
Mono temporal feature and multidate feature are as the alternative features of disease monitoring.
The vegetation characteristics of table 1 are defined
Step 4, disease monitor the phase, synchronized in image capturing unmanned plane low latitude investigation;
Forest reserves image and state that step 5, unmanned plane collection draw the line interior;
Step 6, the image of control module analysis image capture module collection and remote sensing images and light from sensing module
Modal data, pulse radar whether there is insect according to the electromagnetic wave detection being reflected back, and infrared spectrum and multiband spectrum scanner are gathered
The spectral signature of the different-waveband of the forest reserves, and generate collection of illustrative plates;
Step 7, pass the data in step 6 back earth station and further analyze.
Step 8, with reference to unmanned plane low latitude survey data, using image recognition technology, screen the spectral signature of disease monitoring;
It is step 9, different according to forest reserves growth period and whether there is pest and disease damage to be respectively established;
Step 10, a situation arises for the pest and disease damage based on model, spectral information divergence and atlas analysis forestry, the spectrum letter
Breath includes:The spectral signature of the spectral signature of vegetation index, the spectral signature of different-waveband and disease monitoring.
Wherein, the scope of the image wave band includes visible ray and near infrared band.
Wherein, the satellite remote-sensing image of described pair of acquisition carries out pretreatment includes radiation calibration, atmospheric correction, geometric correction
And cloud removal.
Wherein, the coverage of the forestry pests & diseases is according to forest land classification vector figure or multi_temporal images are divided
Class is obtained.The multi_temporal images need to combine forest land use pattern data, terrain data and phenology experience in assorting process, adopt
Taking decision tree, maximum likelihood classification or neutral net carries out planting range extraction.
Wherein, the Spectra feature extraction mode of the disease monitoring is:According to the key area disease that unmanned plane low latitude is investigated
Harmful a situation arises, and key area is divided into normal sample and sample two parts of catching an illness;
The Mono temporal and multidate for extracting two class sample point multi-form spectral signatures from image identification system respectively are special
Value indicative;
Mono temporal or multidate version to every kind of spectral signature, using independent samples t test is relatively normal and sample of catching an illness
This difference degree;
The difference degree of a certain feature is characterized using the preset value p of t inspections, and accordingly generates an all kinds of multi-forms
Spectral signature different phases and when it is combined in p value statistical table, wherein, p value is smaller, normal and sample of catching an illness difference
Bigger, response of the feature to defect information is stronger.
For example:
Combined ground points for investigation data screening disease monitoring spectral signature.According to investigation sampling point disease, a situation arises, will
Sample point is divided into normal sample and sample two parts of catching an illness.According to third step method, two class sample point each spectrum are extracted respectively
The Mono temporal of feature and the characteristic value of multidate version, and calibration is compared using independent samples t test (Independentt_test)
Often with the difference degree of sample of catching an illness, the p value (p-value) checked using t characterizes the difference degree of a certain feature, and can be accordingly
The p that one all kinds of multi-form spectral signature of generation are rented in closing (multidate feature) in different phases (Mono temporal feature) and phase
Data-Statistics form.Wherein, p value is smaller, and the difference of sample of normally catching an illness is bigger, and response of the feature to defect information is stronger.
Therefore, the method for being set by threshold value is screened to all kinds of spectral signatures of Mono temporal and multidate, is retained to defect information
The strong spectral signature form of response.Under normal circumstances, p value can be set the threshold to less than 0.05,0.01 or 0.001 three
Kind.Defect information in image is weaker, and the threshold value of p value sets bigger.
Generally, postponing with phase, each spectral signature is improved constantly in the significance of difference of normal and disease sampling point,
The difference of wherein the 4th phase each feature is the most notable, and now powdery mildew is by the lower shape for having caused wheat strain value from upper infecting
State has significantly change, the even scab up to boot leaf in heavier field is infected in part, this and woodland factual survey
Situation is consistent.By setting a threshold value of P values < 0.001 in this example, it is retained in and normally reaches pole with disease training sample
15 Mono temporals and multidate feature of significant difference level, 15 Mono temporals and multidate feature include respectively:Blue wave band
4th test phase RB (RB-T4), second test phase RB-T4T2 of blue wave band, the 3rd test phase of green wave band
RG-T3, the 4th test phase RG-T4 of green wave band, the 3rd of red wave band tests phase RR-T3, the 4th of red wave band
Test phase RR-T4, the 4th test phase SR-T4 of ratio vegetation index, during the 4th test of vegetation index of normalization difference
Phase NDVI-T4, second and the 4th test phase vegetation index NDVI-T4T2 of normalization difference vegetation index, green wave band is returned
One changes the 4th test phase GNDVI-T4 of difference vegetation index, green wave band normalization difference vegetation index the 4th and second
Test phase GNDVI-T4T2, the 4th test phase OSAVI-T4 of soil regulation vegetation index, improves red side ratio vegetation and refers to
The 3rd test phase NLI-T3 of the 4th test phase MSR-T4 of number and non-linear vegetation index.
The process of the spectral information divergence analysis is:Degree of correlation between by judging two pixels, by picture to be sorted
Degree of correlation highest classification is included into unit.
Below in conjunction with the accompanying drawings and specific embodiment, the present invention is described in further detail.Following examples are used to illustrate
The present invention, is not only not limited to the scope of the present invention.
Realization of the invention is adopted the following technical scheme that:
The first step:Multidate satellite remote-sensing image is ordered and downloads and pre-process.There is process according to most forestry pests & diseases
Fast feature and it is currently available that satellite remote sensing date source, it is proposed that using the middle high-resolution satellite image of revisiting period high.Wave band
Scope needs covering visible light and near infrared band.The most suitable forecasting stage of local forestry pests & diseases is determined first.Obtain forestry disease pest
Evil occurs to the satellite image data of multiple phases in this period.The preprocessing process of image includes radiation calibration, air school
Just, geometric correction and cloud are removed.Follow-up defect information extracts the reflectivity image data based on the multidate obtained after pretreatment
Carry out.
Second step:Application region forest reserves coverage is extracted.Forest land classification vector figure can be combined or according to many
Phase image carries out classification acquisition.Assorting process needs forest land use pattern data in connected applications region, terrain data and thing
Time experience etc., crop-planting scope extraction is carried out using supervised classification methods such as decision tree, maximum likelihood classification or neutral nets.
Follow-up forestry pests & diseases information extraction is obtained being carried out in forest reserves coverage in classification, and it is come to reduce
The interference of his atural object or agrotype.
3rd step:The Mono temporal and multidate Spectra feature extraction of disease monitoring.The method is fitted in domestic object for appreciation that comprehensively light a cigarette
On the basis of the documents and materials of crop disease monitoring, alternative features of 13 spectral signatures as disease monitoring are respectively adopted, wrap
Include the indigo plant (RB) compatible with the multispectral satellite influence of most middle high-resolutions, green (RG), red (RR), near-infrared (RNIR) passage original
Beginning reflectivity, and normalization difference vegetation index (NDVI), ratio vegetation index (SR), green wave band normalization difference vegetation refer to
Number (GNDVI), triangle vegetation index (TVI) improves red side ratio vegetation index (MSR), non-linear vegetation index (NLI), soil
Earth regulation vegetation index (OSAVI) nine broadband vegetation indexs (each index fixed is shown in Table 1).This method is respectively adopted Mono temporal
The Physiology and biochemistry state navigated at certain time point with the vegetation index of multidate breeding version;Multi-temporal vegetation indices are according to two
The Mono temporal vegetation index of individual phase is normalized and is calculated, for reflecting pest and disease damage the characteristics of development and change.Assuming that
There is M and N to carry out two phases for before and after, for a certain vegetation index form V1 (V1 is the name of vegetation index form), multidate
Reflectivity computing formula is (VIN-VIM)/(VIN+VIM).For n phase, phase is carried out into combination of two, it is available altogether
C2n multidate spectral signature.In this method, above-mentioned Mono temporal feature and multidate feature are as the alternative spy of disease monitoring
Levy.
4th step:There is critical period development in disease to be investigated with the unmanned plane low latitude of image capturing time synchronized.Nobody
The investigation of machine low latitude is separated by no more than 3 days with corresponding period satellite shooting date.According to the area of application region, in sampling point setting
The density of 1 sampling point/10 sq-km should be not less than.Meanwhile, total investigation sampling point number should be no less than 30.The field of investigation includes institute
Selectable sampling point is the continuous overlay area of the forest reserves of the diameter more than 30m, and investigation content is gloomy in survey area
The onset grade of woods resource.For ease of woodland investigation and pathology management on a large scale, the forest zone block that will catch an illness be divided into gently, weigh two
Rank.For the different diseases of Different Forest resource, specific disease level delimitation canonical reference disease observes and predicts national standard to be carried out.
5th step:According to the multispectral principle of satellite remote-sensing image, suspicious forestry pests & diseases generation area scope is found out,
In conjunction with GIS forestry topographic maps, the forestry topographic map of suspicious region is marked, and show its state;
6th step:The present invention includes control module and the d GPS locating module, sensing module and the image that are connected with control module
Acquisition module, sensing module includes the multiband spectrum scanning of the pulse radar of monitoring forestry pests & diseases and monitoring forestry pests & diseases
Instrument, infrared spectrometer.Also include humiture module, be connected with control module;Communication module, for connecting main control module and ground
Stand.
Remote sensing technique predicts that pest and disease damage in addition to from spectrum angle analysis, can also recognize disease pest using remote sensing image
The harmful forest reserves.
7th step:By step 6 is by earth station and is mounted in the multispectral image data of monitoring unmanned camera and passes
Earth station is gone back to further to analyze.
8th step:Forest reserves pest and disease damage does not recognize the pest and disease damage initial stage come also in naked eyes, and it is in the outer shadow of rainbow
Dark red tone is showed on picture, and normal crop is then red.Tone difference according to produced by image contrast can be judged
Damaged forest resource.
Because spectral reflectivity of the crop in near-infrared wavelength domain is often very high, when the forest reserves receive pest and disease damage, blade
Water content reduction, cell is collapsed therewith, and chlorophyll is reduced, therefore, their spectral reflectance meeting natures in infrared wavelength domain are reduced,
There is reflection infrared spectrum decay, cause outside rainbow in shooting, the tone of the tone than normal crop of damaged forest resource
Dark, when the normal forest reserves are for red, damaged forest resource is then rendered as kermesinus.When pest and disease damage further develops, make
The chlorophyll of forest reserves blade disappear totally, blade construction is when suffering thoroughly to destroy, image tone can then become darker so that
Adjusted in cyan is presented.According to the difference of damaged forest resource on the outer image of rainbow and normal crop tone, we quickly just can be by
The damaged forest resource judgment that naked eyes are not yet aware of is out, as shown in table 2 below.
The image feature of the healthy forest reserves of table 2 and damaged forest resource
Tenth step:It is all gloomy in large geographical area with reference to GIS technology according to forest reserves growth period different characteristics
The relevant pest and disease damage of woods resource crop rotation refers on same soil, with the time to the forest reserves cycle, in large geographical area
For being diagnosed to be geographical position and the order of severity that all pest and disease damages relevant with forest reserves crop rotation occur, the forest reserves are set up
Different growth times whether the model of pest and disease damage, so that relevant departments are consulted;
11st step:With reference to forestry pests & diseases data screening disease monitoring spectral signature.According to the hair of investigation sampling point disease
Raw situation, normal sample is divided into and sample two parts of catching an illness by sample point.According to the three, the 8th one step process, two classes are extracted respectively
The Mono temporal of sample point each spectral signature and the characteristic value of multidate version, and use independent samples t test
(Independent t-test) relatively more normal and sample of catching an illness difference degree.The p value (p-value) checked using t is represented
The difference degree of a certain feature, and an all kinds of multi-form spectral signatures can be accordingly generated in different phases (Mono temporal feature)
With when combined (multidate feature) in p value statistical table.Wherein, p value is smaller, and normal and sample of catching an illness difference is bigger,
Response month of the feature to defect information is strong.Therefore, by the method for threshold setting to Mono temporal and all kinds of spectrum of multidate
Feature is screened, and is retained and is responded strong spectral signature form to defect information.Under normal circumstances, p can be set the threshold to
Value is less than 0.05,0.01 or 0.001 3 kinds.Defect information in image is weaker, and the threshold value of p value sets bigger.
As shown in figure 3, the A samples of on May 20th, 1,2016 ground unmanned plane multispectral image, extracts normalization difference vegetation index
(NDVI)>0.7 forest reserves region.
2nd, in 1 forest reserves regional extent extracted, then by near infrared spectrum feature extraction (Nir)<0.44 it is gloomy
Woods overlay area.
3rd, in 2 forest reserves regional extents extracted, characters of ground object value (DEM) is extracted by characters of ground object<100m's
Forest reserves region.
4th, the satellite remote-sensing image of on May 1st, 2016 and the unmanned aerial vehicle remote sensing A samples of May 20 ground multispectral image are folded
Plus, carry out sorting out the aggrieved insect pest forest reserves and the healthy forest reserves by Maximum likelihood classification (MLC).
Although the foregoing describing specific embodiment of the invention, those familiar with the art should manage
Solution, the specific embodiment described by us is merely exemplary, and rather than for the restriction to the scope of the present invention, is familiar with this
The technical staff in field should be covered of the invention in the equivalent modification and change made according to spirit of the invention
In scope of the claimed protection.