CN110399860A - A kind of corn damage caused by waterlogging monitoring method and system - Google Patents
A kind of corn damage caused by waterlogging monitoring method and system Download PDFInfo
- Publication number
- CN110399860A CN110399860A CN201910710761.9A CN201910710761A CN110399860A CN 110399860 A CN110399860 A CN 110399860A CN 201910710761 A CN201910710761 A CN 201910710761A CN 110399860 A CN110399860 A CN 110399860A
- Authority
- CN
- China
- Prior art keywords
- corn
- rainfall
- water body
- ndwi
- ndvi
- 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.)
- Granted
Links
- 240000008042 Zea mays Species 0.000 title claims abstract description 159
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 title claims abstract description 157
- 235000002017 Zea mays subsp mays Nutrition 0.000 title claims abstract description 157
- 235000005822 corn Nutrition 0.000 title claims abstract description 157
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000012544 monitoring process Methods 0.000 title claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 128
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 239000000284 extract Substances 0.000 claims abstract description 8
- 239000002689 soil Substances 0.000 claims description 32
- 238000010606 normalization Methods 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 241000209094 Oryza Species 0.000 description 4
- 235000007164 Oryza sativa Nutrition 0.000 description 4
- 235000009566 rice Nutrition 0.000 description 4
- 239000010977 jade Substances 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Animal Husbandry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mining & Mineral Resources (AREA)
- General Business, Economics & Management (AREA)
- Marine Sciences & Fisheries (AREA)
- Human Resources & Organizations (AREA)
- Agronomy & Crop Science (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention discloses a kind of corn damage caused by waterlogging monitoring method and system.This method comprises: obtain the image in different time periods by flooded corn, the different time sections include: rainfall before, in rainfall, after rainfall stopping and rainfall；The image in different time periods is pre-processed, pretreated image is obtained；To the method that the pretreated image uses supervised classification, the planting area of corn crop is obtained；Water body index calculation method is normalized using multidate according to the planting area of the corn crop, extracts corn by flooded region；Corn is divided by flooded region using multidate normalized differential vegetation index calculation method by flooded region according to the corn, obtains the flooded region of corn weight and corn gently flooded region.The flooded region of corn weight and corn gently flooded region can accurately be marked off using method or system of the invention, to be to determine the disaster-stricken range of corn crop, and provide reference frame for assessment corn crop damage situations.
Description
Technical field
The present invention relates to crops damages caused by waterlogging to monitor field, more particularly to a kind of corn damage caused by waterlogging monitoring method and system.
Background technique
The method of previous identification damage caused by waterlogging has Meteorological Index and remote sensing；Meteorological Index is the reality based on meteorological site data
It surveys what rainfall product data was calculated, detection promptly and accurately, Er Qieduo cannot be carried out over time and space to damage caused by waterlogging process
Number Meteorological Index such as SPI needs the station data of many years, and station data distribution is sparse, and using effect is poor.
In terms of remote sensing, the method for being commonly used to extract water body has coloration diagnostic method, single band threshold method, band combination
Method, multispectral Mixed method, Classification of TM Imagery, using relationship between spectrum establish model, DEM takes sounding method.Single band
Middle near infrared band is obvious for the response characteristic of water body, is commonly used to identification, monitoring water body, but identifies in crops damage caused by waterlogging
On, since there are more scramble datas, it cannot reflect the moisture of soil, so that the flooded region of weight and light flood cannot be marked off effectively
Region.
Summary of the invention
The object of the present invention is to provide a kind of corn damage caused by waterlogging monitoring method and systems, can accurately mark off corn weight flood
The light flooded region of region and corn thus for the disaster-stricken range for determining corn crop, and mentions for assessment corn crop damage situations
Foundation for reference.
To achieve the above object, the present invention provides following schemes:
A kind of corn damage caused by waterlogging monitoring method, comprising:
Obtain the image in different time periods by flooded corn, the different time sections include: rainfall before, in rainfall, rainfall
After stopping and rainfall；
The image in different time periods is pre-processed, pretreated image is obtained；
To the method that the pretreated image uses supervised classification, the planting area of corn crop is obtained；
According to the planting area of the corn crop using multidate normalize water body index calculation method, extract corn by
Flooded region；
Corn is divided by flooded region using multidate normalized differential vegetation index calculation method by flooded region according to the corn,
Obtain the flooded region of corn weight and corn gently flooded region.
Optionally, the image in different time periods obtained by flooded corn, specifically includes:
The image in different time periods by flooded corn is obtained by using the WFV sensor of high score No.1 satellite.
Optionally, described that the image in different time periods is pre-processed, pretreated image is obtained, it is specific to wrap
It includes:
Radiation calibration, atmospheric correction and ortho-rectification processing are carried out to the image in different time periods respectively, obtained pre-
Treated image.
Optionally, the method for using supervised classification to the pretreated image, obtains the plantation of corn crop
Region specifically includes:
To the pretreated image using the minimum distance method in supervised classification method, the corn of survey region is extracted
The range of crop；
To the method that the range of the corn crop of the survey region uses blotch removal, the growing area of corn crop is obtained
Domain.
Optionally, the planting area according to corn crop normalizes water body index calculation method using multidate, mentions
It takes corn by flooded region, specifically includes:
Formula is used according to the planting area of the corn cropObtain multidate normalization
Water body index NDWI, the multidate normalization water body index include the water body index NDWI before rainfallBefore, the water body in rainfall refers to
Number NDWIIn, rainfall stop water body index NDWIStopWith the water body index NDWI after rainfallAfterwards, wherein NDWI is water body index,
Green is the green band value of image picture element, and Nir is the near infrared band value of image picture element；
The water body index NDWI stopped according to the rainfallStopWith the water body index NDWI before the rainfallBefore, using formula Δ
NDWI=NDWIStop-NDWIBefore, obtain Soil water diffevence value Δ NDWI；
Obtain the water body index threshold value in rainfall, the water body index threshold value and Soil water diffevence threshold value that rainfall stops；
By the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the soil moisture
Difference value Δ NDWI respectively in the corresponding rainfall water body index threshold value, the rainfall stop water body index threshold value and
The Soil water diffevence threshold value compares, and obtains comparison result；
According to the comparison result, determine corn by flooded region.
Optionally, the water body index NDWI by the rainfallIn, the water body index NDWI that stops of the rainfallStopWith
The water that the Soil water diffevence value Δ NDWI stops with the water body index threshold value in the corresponding rainfall, the rainfall respectively
Body index threshold and the Soil water diffevence threshold value compare, and obtain comparison result, specifically include:
If the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the soil moisture
Difference value Δ NDWI be respectively less than water body index threshold value in the corresponding rainfall, the water body index threshold value that the rainfall stops and
The Soil water diffevence threshold value, it is determined that the planting area of the corn crop is not by flooded region；
If the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the soil moisture
Any value of difference value Δ NDWI stops more than or equal to the water body index threshold value in the corresponding rainfall, the rainfall
Water body index threshold value and the Soil water diffevence threshold value, it is determined that the planting area of the corn crop is by flooded region.
Optionally, described that jade is divided using multidate normalized differential vegetation index calculation method by flooded region according to the corn
Rice is obtained the flooded region of corn weight and corn gently flooded region, specifically included by flooded region:
Formula is used by flooded region according to the cornNormalized differential vegetation index NDVI is obtained,
The normalized differential vegetation index includes the vegetation index NDVI before rainfallBefore, rainfall stop vegetation index NDVIStopAfter rainfall
Vegetation index NDVIAfterwards, wherein Red is the red band value of pixel where corn, and Nir is the near infrared band of pixel where corn
Value；
According to the vegetation index NDVI before the rainfallBeforeWith the vegetation index NDVI after the rainfallAfterwardsUsing formula Δ
NDVI1=NDVIAfterwards-NDVIBefore, obtain the vegetation index difference DELTA NDVI in rainfall later period1；
According to the vegetation index NDVI after the rainfallAfterwardsThe vegetation index NDVI stopped with the rainfallStopUsing formula Δ
NDVI2=NDVIAfterwards-NDVIStop, obtain the vegetation index difference DELTA NDVI of rainfall early period2；
To the vegetation index difference DELTA NDVI in the rainfall later period1With the vegetation index difference DELTA NDVI of the rainfall early period2Point
It Cai Yong not formulaWith
It is standardized, the vegetation index difference ST Δ NDVI in the rainfall later period after being standardized1With the rainfall after standardization
The vegetation index difference ST Δ NDVI of early period2；
According to the vegetation index difference ST Δ NDVI in the rainfall later period after the standardization1With the rainfall after the standardization
The vegetation index difference ST Δ NDVI of early period2Using formulaObtain combined influence
Difference | Δ D |；
According to the combined influence difference | Δ D | use formulaObtain standard summary
Influence difference ST | Δ D |；
When the standard summary influence difference ST | Δ D | be less than given threshold when, determine the corn by flooded region for jade
The flooded region of rice weight；
When the standard summary influences difference ST | Δ D | when being greater than or equal to given threshold, determine the corn by flooded area
Domain is corn gently flooded region.
A kind of corn damage caused by waterlogging monitoring system, comprising:
Module is obtained, for obtaining the image in different time periods by flooded corn, the different time sections include: rainfall
Before, in rainfall, after rainfall stopping and rainfall；
Preprocessing module obtains pretreated image for pre-processing to the image in different time periods；
Categorization module, the method for using supervised classification to the pretreated image, obtains the kind of corn crop
Growing area domain；
Water body index computing module, for being referred to according to the planting area of the corn crop using multidate normalization water body
Number calculating method extracts corn by flooded region；
Vegetation index computing module, for being calculated by flooded region using multidate normalized differential vegetation index according to the corn
Method divides corn by flooded region, obtains the flooded region of corn weight and corn gently flooded region.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
The present invention provides a kind of corn damage caused by waterlogging monitoring method, comprising: obtains the image in different time periods by flooded corn, no
Before including: rainfall with the period, in rainfall, after rainfall stopping and rainfall；Image in different time periods is pre-processed, is obtained
Pretreated image；To the method that pretreated image uses supervised classification, the planting area of corn crop is obtained；According to
The planting area of corn crop normalizes water body index calculation method using multidate, extracts corn by flooded region；According to corn
Corn is divided by flooded region using multidate normalized differential vegetation index calculation method by flooded region, obtains the flooded region of corn weight and jade
The light flooded region of rice.The flooded region of corn weight and corn gently flooded region can accurately be marked off using method of the invention, to be
It determines the disaster-stricken range of corn crop, and provides reference frame for assessment corn crop damage situations.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is corn damage caused by waterlogging monitoring method flow chart of the present invention；
Fig. 2 is that corn damage caused by waterlogging of the present invention monitors system construction drawing.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of corn damage caused by waterlogging monitoring method and systems, can accurately mark off corn weight flood
The light flooded region of region and corn thus for the disaster-stricken range for determining corn crop, and mentions for assessment corn crop damage situations
Foundation for reference.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is corn damage caused by waterlogging monitoring method flow chart of the present invention.The corn damage caused by waterlogging monitoring method includes:
Step 101: obtain the image in different time periods by flooded corn, different time sections include: rainfall before, in rainfall,
After rainfall stopping and rainfall.
By using high score No.1 satellite WFV sensor obtain before by flooded corn rainfall, in rainfall, rainfall stop and drop
Image after rain；The image data spatial resolution is 16m, including 4 wave bands, respectively blue wave band (0.45~0.25 μm),
Green wave band (0.52~0.59 μm), red wave band (0.63~0.69 μm), nearly red wave band (0.77~0.89 μm).
Step 102: image in different time periods is pre-processed, pretreated image is obtained, specifically includes:
It carries out radiation calibration, atmospheric correction and ortho-rectification processing respectively to image in different time periods, is pre-processed
Image afterwards.
Step 103: the method for using supervised classification to pretreated image obtains the planting area of corn crop, has
Body includes:
To pretreated image using the minimum distance method in supervised classification method, the corn crop of survey region is extracted
Range.
To the method that the range of the corn crop of survey region is removed using blotch removal or manually, corn crop is obtained
Planting area.
Step 104: water body index calculation method being normalized using multidate according to the planting area of corn crop, is extracted beautiful
Rice is specifically included by flooded region:
Formula is used according to the planting area of corn cropObtain multidate normalization water body
Index NDWI, multidate normalization water body index NDWI include the water body index NDWI before rainfallBefore, water body index in rainfall
NDWIIn, rainfall stop water body index NDWIStopWith the water body index NDWI after rainfallAfterwards, wherein NDWI is multidate normalization
Water body index, Green are the green band value of image picture element, and Nir is the near infrared band value of image picture element；Previous research table
Bright, near-infrared and visible light can protrude water body for the different response contrasts of water body well, and NDWI uses near infrared band
With the prominent water body of green band combination, vegetation information can be inhibited well.
The water body index NDWI stopped according to rainfallStopWith the water body index NDWI before rainfallBefore, using formula Δ NDWI=
NDWIStop-NDWIBefore, obtain Soil water diffevence value Δ NDWI；Δ NDWI considers the difference of soil moisture before rainfall, and reflects
The retention degree of the extra water body of earth's surface after rainfall.
Obtain the water body index threshold value in rainfall, the water body index threshold value and Soil water diffevence threshold value that rainfall stops.
By the water body index NDWI in rainfallIn, rainfall stop water body index NDWIStopWith Soil water diffevence value Δ NDWI
The water body index threshold value and Soil water diffevence threshold value stopped respectively with the water body index threshold value in corresponding rainfall, rainfall does ratio
Compared with obtaining comparison result.
According to comparison result, determine corn by flooded region.
Wherein, by the water body index NDWI in rainfallIn, rainfall stop water body index NDWIStopWith Soil water diffevence value
The water body index threshold value and Soil water diffevence threshold that Δ NDWI stops with the water body index threshold value in corresponding rainfall, rainfall respectively
Value compares, and obtains comparison result, specifically includes:
If the water body index NDWI in rainfallIn, rainfall stop water body index NDWIStopWith Soil water diffevence value Δ NDWI
The water body index threshold value and Soil water diffevence threshold value that water body index threshold value, rainfall in respectively less than corresponding rainfall stop, then
The planting area for determining corn crop is not by flooded region.
If the water body index NDWI in rainfallIn, rainfall stop water body index NDWIStopWith Soil water diffevence value Δ NDWI
Any value is greater than or equal to water body index threshold value, the water body index threshold value and the soil water of rainfall stopping in corresponding rainfall
Divide discrepancy threshold, it is determined that the planting area of corn crop is by flooded region.
Step 105: corn is divided by flood using multidate normalized differential vegetation index calculation method by flooded region according to corn
Region obtains the flooded region of corn weight and corn gently flooded region, specifically includes:
Formula is used by flooded region according to cornObtain normalized differential vegetation index NDVI, normalizing
Changing vegetation index NDVI includes the vegetation index NDVI before rainfallBefore, rainfall stop vegetation index NDVIStopWith the vegetation after rainfall
Index NDVIAfterwards, wherein Red is the red band value of pixel where corn, and Nir is the near infrared band value of pixel where corn.
According to the vegetation index NDVI before rainfallBeforeWith the vegetation index NDVI after rainfallAfterwardsUsing formula Δ NDVI1=
NDVIAfterwards-NDVIBefore, obtain the vegetation index difference DELTA NDVI in rainfall later period1。
According to the vegetation index NDVI after rainfallAfterwardsThe vegetation index NDVI stopped with rainfallStopUsing formula Δ NDVI2=
NDVIAfterwards-NDVIStop, obtain the vegetation index difference DELTA NDVI of rainfall early period2。
To the vegetation index difference DELTA NDVI in rainfall later period1With the vegetation index difference DELTA NDVI of rainfall early period2It is respectively adopted
FormulaWithIt is marked
Quasi-ization processing, the vegetation index difference ST Δ NDVI in the rainfall later period after being standardized1With rainfall early period after standardization
Vegetation index difference ST Δ NDVI2。
According to the vegetation index difference ST Δ NDVI in the rainfall later period after standardization1With the plant of rainfall early period after standardization
Formula is used by index difference value ST Δ NDVI2Obtain combined influence difference | Δ D
|, | Δ D | indicate comprehensive ST △ NDVI1With the influence of ST △ two differences of NDVI2.
According to the combined influence difference | Δ D | use formulaObtain standard summary
Influence difference ST Δ D |.Single NDVI just has difference before moisture is affected, and uses the difference of the NDVI of two phases can be with
Eliminate the influence of corn normal growth, therefore the difference of comprehensive two NDVI is as a result, introduce ST Δ D |, ST Δ D | value 0~1
Between.It can be made to calculate duration all between 0 to 1 in any region ∣ △ D ∣ standardization, and consider itself and territory of use
The difference of maximum value and minimum value uses it with more popularity in other regions.
When the standard summary influence difference ST | Δ D | be less than given threshold when, determine corn by flooded region be corn weight
Flooded region.
When the standard summary influences difference ST | Δ D | when being greater than or equal to given threshold, determine that corn is by flooded region
The light flooded region of corn.
Given threshold may be provided in 0.4 to 0.6 range according to different regions, if there is sample point data can accomplish
More accurate setting.
Due to the otherness of Regional Rainfall, the time of rainfall is largely concentrated between different regions and stops the time of rainfall
Difference, for concentrating rainfall will be earlier than concentration rainfall close to rear half of rainfall close to the regional response time of preceding half of rainfall cycle
The area in period, the regional influence in order to consider this factor long for rainfall cycle, selects two NDVI difference results
Reflect the influence of rainfall early period and later period.STΔNDVI1、STΔNDVI2, ST | Δ D | these three indexes can be respectively intended to draw
Divide threshold value, but combines ST Δ NDVI1、STΔNDVI2And the ST calculated | Δ D | it is more preferable in precision.
Corn in various degree can have some impact on yield by flood, by extracting corn by flooded range and dividing by flood
Degree can settle a claim for agricultural insurance and provide foundation.
Fig. 2 is that corn damage caused by waterlogging of the present invention monitors system construction drawing.The corn damage caused by waterlogging monitors system
Obtain module 201, for obtaining the image in different time periods by flooded corn, different time sections include: rainfall before,
In rainfall, after rainfall stopping and rainfall；
Preprocessing module 202 obtains pretreated image for pre-processing to image in different time periods；
Categorization module 203, the method for using supervised classification to pretreated image, obtains the plantation of corn crop
Region；
Water body index computing module 204, for being referred to according to the planting area of corn crop using multidate normalization water body
Number calculating method extracts corn by flooded region；
Vegetation index computing module 205, for being calculated by flooded region using multidate normalized differential vegetation index according to corn
Method divides corn by flooded region, obtains the flooded region of corn weight and corn gently flooded region.
The above method not only monitors corn damage caused by waterlogging and is applicable in, while being also suitable to the monitoring of the damage caused by waterlogging of other crops.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other
The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment
For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part
It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention；At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.
Claims (8)
1. a kind of corn damage caused by waterlogging monitoring method characterized by comprising
Obtain the image in different time periods by flooded corn, the different time sections include: rainfall before, in rainfall, rainfall stops
After rainfall；
The image in different time periods is pre-processed, pretreated image is obtained；
To the method that the pretreated image uses supervised classification, the planting area of corn crop is obtained；
Water body index calculation method is normalized using multidate according to the planting area of the corn crop, extracts corn by flooded area
Domain；
Corn is divided by flooded region using multidate normalized differential vegetation index calculation method by flooded region according to the corn, is obtained
Corn weighs flooded region and corn gently flooded region.
2. corn damage caused by waterlogging monitoring method according to claim 1, which is characterized in that when the difference obtained by flooded corn
Between section image, specifically include:
The image in different time periods by flooded corn is obtained by using the WFV sensor of high score No.1 satellite.
3. corn damage caused by waterlogging monitoring method according to claim 1, which is characterized in that described to the shadow in different time periods
Picture is pre-processed, and is obtained pretreated image, is specifically included:
Radiation calibration, atmospheric correction and ortho-rectification processing are carried out to the image in different time periods respectively, pre-processed
Image afterwards.
4. corn damage caused by waterlogging monitoring method according to claim 1, which is characterized in that described to the pretreated image
Using the method for supervised classification, the planting area of corn crop is obtained, is specifically included:
To the pretreated image using the minimum distance method in supervised classification method, the corn crop of survey region is extracted
Range；
To the method that the range of the corn crop of the survey region uses blotch removal, the planting area of corn crop is obtained.
5. corn damage caused by waterlogging monitoring method according to claim 1, which is characterized in that the growing area according to corn crop
Domain normalizes water body index calculation method using multidate, extracts corn by flooded region, specifically includes:
Formula is used according to the planting area of the corn cropObtain multidate normalization water body
Index NDWI, the multidate normalization water body index NDWI include the water body index NDWI before rainfallBefore, the water body in rainfall refers to
Number NDWIIn, rainfall stop water body index NDWIStopWith the water body index NDWI after rainfallAfterwards, wherein NDWI is multidate normalizing
Change water body index, Green is the green band value of image picture element, and Nir is the near infrared band value of image picture element；
The water body index NDWI stopped according to the rainfallStopWith the water body index NDWI before the rainfallBefore, using formula Δ NDWI
=NDWIStop-NDWIBefore, obtain Soil water diffevence value Δ NDWI；
Obtain the water body index threshold value in rainfall, the water body index threshold value and Soil water diffevence threshold value that rainfall stops；
By the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the Soil water diffevence
It is worth the water body index threshold value and described that Δ NDWI stops with the water body index threshold value in the corresponding rainfall, the rainfall respectively
Soil water diffevence threshold value compares, and obtains comparison result；
According to the comparison result, determine corn by flooded region.
6. corn damage caused by waterlogging monitoring method according to claim 5, which is characterized in that the water body by the rainfall refers to
Number NDWIIn, the water body index NDWI that stops of the rainfallStopWith the Soil water diffevence value Δ NDWI respectively with it is corresponding described
The water body index threshold value and the Soil water diffevence threshold value that water body index threshold value, the rainfall in rainfall stop compare,
Comparison result is obtained, is specifically included:
If the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the Soil water diffevence
Value Δ NDWI is respectively less than water body index threshold value in the corresponding rainfall, the water body index threshold value that the rainfall stops and described
Soil water diffevence threshold value, it is determined that the planting area of the corn crop is not by flooded region；
If the water body index NDWI in the rainfallIn, the water body index NDWI that stops of the rainfallStopWith the Soil water diffevence
It is worth any value of Δ NDWI more than or equal to the water body that the water body index threshold value in the corresponding rainfall, the rainfall stop
Index threshold and the Soil water diffevence threshold value, it is determined that the planting area of the corn crop is by flooded region.
7. corn damage caused by waterlogging monitoring method according to claim 1, which is characterized in that it is described according to the corn by flooded region
Corn is divided by flooded region using multidate normalized differential vegetation index calculation method, obtains the flooded region of corn weight and the area corn Qing Lao
Domain specifically includes:
Formula is used by flooded region according to the cornNormalized differential vegetation index NDVI is obtained, it is described
Normalized differential vegetation index NDVI includes the vegetation index NDVI before rainfallBefore, rainfall stop vegetation index NDVIStopAfter rainfall
Vegetation index NDVIAfterwards, wherein Red is the red band value of pixel where corn, and Nir is the near infrared band of pixel where corn
Value；
According to the vegetation index NDVI before the rainfallBeforeWith the vegetation index NDVI after the rainfallAfterwardsUsing formula Δ NDVI1=
NDVIAfterwards-NDVIBefore, obtain the vegetation index difference DELTA NDVI in rainfall later period1；
According to the vegetation index NDVI after the rainfallAfterwardsThe vegetation index NDVI stopped with the rainfallStopUsing formula Δ NDVI2
=NDVIAfterwards-NDVIStop, obtain the vegetation index difference DELTA NDVI of rainfall early period2；
To the vegetation index difference DELTA NDVI in the rainfall later period1Distinguish with the vegetation index difference DELTA NDVI2 of the rainfall early period
Using formulaWithInto
Row standardization, the vegetation index difference ST Δ NDVI in the rainfall later period after being standardized1Before the rainfall after standardization
The vegetation index difference ST Δ NDVI of phase2；
According to the vegetation index difference ST Δ NDVI in the rainfall later period after the standardization1With rainfall early period after the standardization
Vegetation index difference ST Δ NDVI2Using formulaObtain combined influence difference
|ΔD|；
According to combined influence difference | Δ D | use formulaObtaining standard summary influences difference
ST|ΔD|；
When the standard summary influence difference ST | Δ D | be less than given threshold when, determine the corn by flooded region be corn weight
Flooded region；
When the standard summary influences difference ST | Δ D | when being greater than or equal to given threshold, determine that the corn is by flooded region
The light flooded region of corn.
8. a kind of corn damage caused by waterlogging monitors system characterized by comprising
Obtain module, for obtaining the image in different time periods by flooded corn, the different time sections include: rainfall before, drop
In the rain, after rainfall stopping and rainfall；
Preprocessing module obtains pretreated image for pre-processing to the image in different time periods；
Categorization module, the method for using supervised classification to the pretreated image, obtains the growing area of corn crop
Domain；
Water body index computing module, by being normalized based on water body index according to the planting area of the corn crop using multidate
Calculation method extracts corn by flooded region；
Vegetation index computing module, for being used multidate normalized differential vegetation index calculation method by flooded region according to the corn
Corn is divided by flooded region, obtains the flooded region of corn weight and corn gently flooded region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710761.9A CN110399860B (en) | 2019-08-02 | 2019-08-02 | Corn flood monitoring method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910710761.9A CN110399860B (en) | 2019-08-02 | 2019-08-02 | Corn flood monitoring method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110399860A true CN110399860A (en) | 2019-11-01 |
| CN110399860B CN110399860B (en) | 2021-04-09 |
Family
ID=68327283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910710761.9A Active CN110399860B (en) | 2019-08-02 | 2019-08-02 | Corn flood monitoring method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110399860B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113205006A (en) * | 2021-04-12 | 2021-08-03 | 武汉大学 | Multi-temporal remote sensing image rice extraction method based on rice indexes |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070003107A1 (en) * | 2005-07-01 | 2007-01-04 | Deere & Company, A Delaware Corporation | Method and system for vehicular guidance using a crop image |
| CN105389559A (en) * | 2015-11-12 | 2016-03-09 | 中国科学院遥感与数字地球研究所 | System and method for identifying agricultural disaster scope based on high-resolution remote sensing image |
| US9292747B2 (en) * | 2013-03-15 | 2016-03-22 | The Boeing Company | Methods and systems for automatic and semi-automatic geometric and geographic feature extraction |
| CN105809140A (en) * | 2016-03-18 | 2016-07-27 | 华南农业大学 | Method and device for extracting surface water body information based on remote sensing model |
| CN106022341A (en) * | 2016-05-03 | 2016-10-12 | 浙江海洋大学 | High resolution optical remote sensing image post-disaster water body information extracting method and system |
| CN107067334A (en) * | 2017-04-21 | 2017-08-18 | 中国科学院遥感与数字地球研究所 | Farmland fire and crop straw burning monitoring method and system based on multi- source Remote Sensing Data data |
| CN107966116A (en) * | 2017-11-20 | 2018-04-27 | 苏州市农业科学院 | The remote-sensing monitoring method and system of a kind of Monitoring of Paddy Rice Plant Area |
| CN108229356A (en) * | 2017-12-23 | 2018-06-29 | 航天恒星科技有限公司 | Dynamic integrity natural calamity remote-sensing monitoring method |
| CN108256534A (en) * | 2018-01-29 | 2018-07-06 | 中国科学院地理科学与资源研究所 | A kind of raft formula marine cultivation region extracting method based on high score remote sensing image |
| CN108985260A (en) * | 2018-08-06 | 2018-12-11 | 航天恒星科技有限公司 | A kind of remote sensing and meteorological integrated rice yield estimation method |
| CN109614891A (en) * | 2018-11-27 | 2019-04-12 | 北京师范大学 | Crops recognition methods based on phenology and remote sensing |
| CN110059890A (en) * | 2019-04-26 | 2019-07-26 | 北京观微科技有限公司 | County Scale agricultural flood monitor method and system |
-
2019
- 2019-08-02 CN CN201910710761.9A patent/CN110399860B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070003107A1 (en) * | 2005-07-01 | 2007-01-04 | Deere & Company, A Delaware Corporation | Method and system for vehicular guidance using a crop image |
| US9292747B2 (en) * | 2013-03-15 | 2016-03-22 | The Boeing Company | Methods and systems for automatic and semi-automatic geometric and geographic feature extraction |
| CN105389559A (en) * | 2015-11-12 | 2016-03-09 | 中国科学院遥感与数字地球研究所 | System and method for identifying agricultural disaster scope based on high-resolution remote sensing image |
| CN105809140A (en) * | 2016-03-18 | 2016-07-27 | 华南农业大学 | Method and device for extracting surface water body information based on remote sensing model |
| CN106022341A (en) * | 2016-05-03 | 2016-10-12 | 浙江海洋大学 | High resolution optical remote sensing image post-disaster water body information extracting method and system |
| CN107067334A (en) * | 2017-04-21 | 2017-08-18 | 中国科学院遥感与数字地球研究所 | Farmland fire and crop straw burning monitoring method and system based on multi- source Remote Sensing Data data |
| CN107966116A (en) * | 2017-11-20 | 2018-04-27 | 苏州市农业科学院 | The remote-sensing monitoring method and system of a kind of Monitoring of Paddy Rice Plant Area |
| CN108229356A (en) * | 2017-12-23 | 2018-06-29 | 航天恒星科技有限公司 | Dynamic integrity natural calamity remote-sensing monitoring method |
| CN108256534A (en) * | 2018-01-29 | 2018-07-06 | 中国科学院地理科学与资源研究所 | A kind of raft formula marine cultivation region extracting method based on high score remote sensing image |
| CN108985260A (en) * | 2018-08-06 | 2018-12-11 | 航天恒星科技有限公司 | A kind of remote sensing and meteorological integrated rice yield estimation method |
| CN109614891A (en) * | 2018-11-27 | 2019-04-12 | 北京师范大学 | Crops recognition methods based on phenology and remote sensing |
| CN110059890A (en) * | 2019-04-26 | 2019-07-26 | 北京观微科技有限公司 | County Scale agricultural flood monitor method and system |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113205006A (en) * | 2021-04-12 | 2021-08-03 | 武汉大学 | Multi-temporal remote sensing image rice extraction method based on rice indexes |
| CN113205006B (en) * | 2021-04-12 | 2022-07-19 | 武汉大学 | Multi-temporal remote sensing image rice extraction method based on rice indexes |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110399860B (en) | 2021-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hufkens et al. | Monitoring crop phenology using a smartphone based near-surface remote sensing approach | |
| Qiao et al. | UAV-based chlorophyll content estimation by evaluating vegetation index responses under different crop coverages | |
| Zeng et al. | A hybrid approach for detecting corn and soybean phenology with time-series MODIS data | |
| US9734400B2 (en) | System and method for field variance determination | |
| Gnyp et al. | Hyperspectral canopy sensing of paddy rice aboveground biomass at different growth stages | |
| CN110796001B (en) | Satellite image film-covering farmland identification and extraction method and system | |
| CN109389049A (en) | Crop Classification in Remote Sensing Image method based on multidate SAR data and multispectral data | |
| CN113850139B (en) | Multi-source remote sensing-based forest annual phenological monitoring method | |
| CN106780091A (en) | Agricultural disaster information remote sensing extracting method based on vegetation index time space statistical nature | |
| CN106897707A (en) | Characteristic image time series synthetic method and device based in multi-source points | |
| CN116188465B (en) | Crop growth state detection method based on image processing technology | |
| CN112085781B (en) | Method for extracting winter wheat planting area based on spectrum reconstruction technology | |
| Wang et al. | Mapping paddy rice with the random forest algorithm using MODIS and SMAP time series | |
| CN103969632A (en) | Device and method of using radar remote sensing data for monitoring wheat lodging | |
| CN109211791A (en) | Crop condition monitoring method and system | |
| Li et al. | Rapid diagnosis of agricultural soil health: A novel soil health index based on natural soil productivity and human management | |
| Jiang et al. | Mapping interannual variability of maize cover in a large irrigation district using a vegetation index–phenological index classifier | |
| CN107273797B (en) | Rice sub-pixel identification method based on water body index variation coefficient | |
| Wang et al. | An automated extraction of small-and middle-sized rice fields under complex terrain based on SAR time series: A case study of Chongqing | |
| Beneduzzi et al. | Temporal variability in active reflectance sensor-measured NDVI in soybean and wheat crops | |
| Kumari et al. | Soybean cropland mapping using multi-temporal sentinel-1 data | |
| Yang et al. | Fraction vegetation cover extraction of winter wheat based on RGB image obtained by UAV | |
| Pauly | Applying conventional vegetation vigor indices to UAS-derived orthomosaics: issues and considerations | |
| CN110399860A (en) | A kind of corn damage caused by waterlogging monitoring method and system | |
| CN117274798B (en) | Remote sensing rice identification method based on regularized time sequence variation model |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |