CN109471131A - It is lain fallow the method and apparatus of situation by remote sensing satellite photo statistical monitoring crop rotation - Google Patents
It is lain fallow the method and apparatus of situation by remote sensing satellite photo statistical monitoring crop rotation Download PDFInfo
- Publication number
- CN109471131A CN109471131A CN201811363119.XA CN201811363119A CN109471131A CN 109471131 A CN109471131 A CN 109471131A CN 201811363119 A CN201811363119 A CN 201811363119A CN 109471131 A CN109471131 A CN 109471131A
- Authority
- CN
- China
- Prior art keywords
- photo
- ndvi
- days
- condition
- remote sensing
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- 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
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Abstract
It is lain fallow the method and apparatus of situation this application involves a kind of by remote sensing satellite photo statistical monitoring crop rotation, by obtaining annual satellite remote sensing photo of the areal current year by the end of June to second year by the end of June, pass through the extraction to NDVI value, Green value and Blue value is carried out in satellite remote sensing photo, and judged according to 6 conditions, finally can determine that whether corresponding region is to lie fallow to sun the earth which has been ploughed up to lie fallow with crop rotation, have the advantages that result is accurate, realizes regional arable land crop rotation and lie fallow using plot as the accurate monitoring and surveying of the remote sensing of unit.
Description
Technical field
The application belongs to remote sensing monitoring technical field, passes through remote sensing satellite photo statistical monitoring wheel more particularly, to one kind
The method and apparatus for situation of lying fallow.
Background technique
Currently, China's agricultural is faced with the problem of pattern of farming is unreasonable, and resource environment overloads.Carry out arable land crop rotation to stop
Ploughing is to promote agriculture Green Development, realizes " hiding grain is in ground, hiding grain in skill " target, ensures the crucial behave of permanent grain security.
Crop rotation (Crop rotation) refers to that on same arable land, some cycles lubrication groove turns growing different crops or combined plantation side
Formula (Zhao Qiguo, 2017)；Lie fallow (Land fallow) refer to arable land stopped in Growing Season of Crops to plough mode (Liu's xun, one of the Eight Diagrams for not planting
It is great, 1996).Crop rotation, which is lain fallow, can be divided into crop rotation, lie fallow and sun the earth which has been ploughed up and lie fallow the Three models of fertilizing.As crop rotation is lain fallow work
The development of work, remote sensing technology are lain fallow benefit compensation and effect because it obtains that information is fast, obtains the characteristics of abundant information and have become crop rotation
The important means of fruit assessment.It is expanded using the variation of remote sensing technology tracking and monitoring arable land long-term cropping and is widely ground both at home and abroad
Study carefully application.
Currently, various regions arable land situation is different, crop rotation is lain fallow situation difference, in the remote sensing monitoring for being directly facing crop rotation and lying fallow
Still lack perfect method system.
Summary of the invention
The technical problem to be solved by the present invention is to solve deficiency in the prior art, so that providing one kind passes through remote sensing
Satellite photo statistical monitoring crop rotation is lain fallow the method and apparatus of situation.
The technical solution adopted by the present invention to solve the technical problems is:
A method of it is lain fallow situation by remote sensing satellite photo statistical monitoring crop rotation, comprising the following steps:
S1: areal current year 1 year satellite remote sensing photo to second year by the end of June by the end of June is obtained, to satellite remote sensing
Arable land plot in photo carries out piecemeal or without piecemeal；
S2: the NDVI value of each piecemeal in each photo, the average value of Green value and Blue value, NDVI=are obtained
(ρNIR-ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρR+ρG+ρB), ρNIR、ρR、ρG、ρBIt respectively indicates
Near infrared band reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；When without piecemeal, obtain every
NDVI value, Green value and the Blue value of each pixel in a photo；
S3: carrying out the judgement of NDVI value, Green value and Blue value to each piecemeal in satellite remote sensing photo or pixel,
And piecemeal or pixel are demarcated, judgment basis is based on the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
It ploughs rape and green manuring is lain fallow for crop rotation；
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI (August 13-
Photo on the 15th)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24,
NDVI (photo on 10-12 days November) < 0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)
>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI (May 13-15
The photo of day) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI (May 13-15
Day photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (May
3-5 days photos) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March) < 0.40, NDVI (photograph on 8-10 days April
Piece)>0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI (May 13-15
The photo of day) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32；
S4: summarize piecemeal or pixel as a result, and exporting image.
Preferably, of the invention to be lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, in S3 step
In, eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow.
Preferably, of the invention to be lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, to satellite remote sensing
When photo is without piecemeal, in S4 step, using n × n as pixel window, moving step length 1, threshold value is that m divides adjoining
Block aggregates into bulk, will be filled to be formed inside regular edges continuously by the fritter for being unsatisfactory for condition surrounded in bulk
Piecemeal, and export image.
Preferably, of the invention to be lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, the n=3 or
4 or 5, m=4 or 5 or 6.
Preferably, of the invention to be lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, to satellite remote sensing
Arable land plot in photo carries out piecemeal and identifies to obtain using multi-scale segmentation method by the ridge road to arable land.
The application also provide it is a kind of lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, including, image
Acquisition module, image processing module, classification demarcating module and result output module；
Described image acquisition module is used to obtain annual satellite of the areal current year by the end of October to second year by the end of October
Remote sensing photo；
Described image processing module is used to obtain NDVI value, Green value and the Blue of each pixel in each photo
Value, or piecemeal first is carried out to the arable land plot in satellite remote sensing photo and obtains the NDVI of each piecemeal in each photo
The average value of value, Green value and Blue value,
Wherein, NDVI=(ρNIR-ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρR+ρG+ρB),
ρNIR、ρR、ρG、ρBRespectively indicate near infrared band reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；
The classification demarcating module is for demarcating the pixel in satellite remote sensing photo, the judgment basis base of calibration
In the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
It ploughs rape and green manuring is lain fallow for crop rotation；
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI (August 13-
Photo on the 15th)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24,
NDVI (photo on 10-12 days November) < 0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)
>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI (May 13-15
The photo of day) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI (May 13-15
Day photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (May
3-5 days photos) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March) < 0.40, NDVI (photograph on 8-10 days April
Piece)>0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI (May 13-15
The photo of day) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32；
The result output module is used to export the image demarcated according to the result of the classification demarcating module.
Preferably, of the invention to be lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, the contingency table
In cover half block, eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow.
Preferably, of the invention to be lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, at described image
When reason module is used to obtain the NDVI value, Green value and Blue value of each pixel in each photo, the result exports mould
In block, using n × n as pixel window, adjoining piecemeal is aggregated into bulk for m by moving step length 1, threshold value, will be by bulk
The fritter for being unsatisfactory for condition surrounded is filled to form continuous piecemeal inside regular edges, and exports image.
Preferably, of the invention to be lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, the n=3 or
4 or 5, m=4 or 5 or 6.
Preferably, of the invention to be lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, to satellite remote sensing
Arable land plot in photo carries out piecemeal and identifies to obtain using multi-scale segmentation method by the ridge road to arable land.
The beneficial effects of the present invention are:
The application is led to by obtaining annual satellite remote sensing photo of the areal current year by the end of June to second year by the end of June
Cross the extraction to NDVI value, Green value and Blue value is carried out in satellite remote sensing photo, and judged according to 6 conditions,
Finally can determine that whether corresponding region is to lie fallow to sun the earth which has been ploughed up to lie fallow with crop rotation, has the advantages that result is accurate, realizes region
Property arable land crop rotation lie fallow using plot as the accurate monitoring and surveying of the remote sensing of unit.
Detailed description of the invention
The technical solution of the application is further illustrated with reference to the accompanying drawings and examples.
Fig. 1 is the condition criterion rule of the embodiment of the present application；
Fig. 2 is Suzhou districts under city administration elevation distribution remote sensing image data using European Space Agency (ESA) sentry 2 (Sentinel-2)
Data；
Fig. 3 is arable land plot piecemeal result figure；
Fig. 4 is that the Suzhou 2017-2018 districts under city administration wheel for ploughing makees the mode lain fallow distribution；
Fig. 5 is that the Suzhou 2017-2018 districts under city administration wheel for ploughing is lain fallow stubble plantation distribution after rice.
Specific embodiment
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase
Mutually combination.
It is described in detail the technical solution of the application below with reference to the accompanying drawings and in conjunction with the embodiments.
Embodiment 1
It is lain fallow the method for situation the present embodiment provides a kind of by remote sensing satellite photo statistical monitoring crop rotation, including following
Step:
S1: areal current year 1 year satellite remote sensing photo to second year by the end of June by the end of June is obtained, to satellite remote sensing
Arable land plot in photo carries out piecemeal；(for example utilize the 0.8m resolution image of GF-2 based on the convolution algorithm of machine learning
Ground block message in arable land is identified, it, can also be further combined with to the ridge road in arable land using multi-scale segmentation method identification
High-resolution Google Earth image carries out artificial quality inspection and amendment to recognition result)；
S2: the NDVI value of each piecemeal in each photo, the average value of Green value and Blue value, NDVI=are obtained
(ρNIR-ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρR+ρG+ρB), ρNIR、ρR、ρG、ρBIt respectively indicates
Near infrared band reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；
S3: the judgement of NDVI value, Green value and Blue value is carried out to piecemeal each in satellite remote sensing photo, and to piecemeal
It is demarcated, judgment basis is based on the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
Eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow；
It ploughs rape and green manuring is lain fallow for crop rotation；
Different vegetation have different spectral signatures, and NDVI (normalized differential vegetation index) can intuitively reflect vegetation not
Same coverage condition.The feature restriction of multidate NDVI is one of the effective ways that different vegetation are extracted.And in rape full blossom
Its bluish-green color accounting of phase has significant difference with wheat.
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI (August 13-
Photo on the 15th)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24,
NDVI (photo on 10-12 days November) < 0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)
>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI (May 13-15
The photo of day) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI (May 13-15
Day photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (May
3-5 days photos) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March) < 0.40, NDVI (photograph on 8-10 days April
Piece)>0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI (May 13-15
The photo of day) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32.
S4: demarcate piecemeal as a result, and export image, can be obtained the area of various situations according to image.
Embodiment 2
It is lain fallow the method for situation the present embodiment provides a kind of by remote sensing satellite photo statistical monitoring crop rotation, including following
Step:
S1: areal current year 1 year satellite remote sensing photo to second year by the end of June by the end of June is obtained；
S2: the NDVI value, Green value and Blue value of each pixel in each photo, NDVI=(ρ are obtainedNIR-ρR)/
(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρ R+ρG+ρB), ρNIR、ρR、ρG、ρBRespectively indicate near infrared band
Reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；
S3: the judgement of NDVI value, Green value and Blue value, judgment basis are carried out to pixel each in satellite remote sensing photo
Based on the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
Eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow；
It ploughs rape and green manuring is lain fallow for crop rotation；
Different vegetation have different spectral signatures, and NDVI (normalized differential vegetation index) can intuitively reflect vegetation not
Same coverage condition.The feature restriction of multidate NDVI is one of the effective ways that different vegetation are extracted.And in rape full blossom
Its bluish-green color accounting of phase has significant difference with wheat.
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI (August 13-
Photo on the 15th)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24,
NDVI (photo on 10-12 days November) < 0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)
>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI (May 13-15
The photo of day) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI (May 13-15
Day photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (May
3-5 days photos) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March) < 0.40, NDVI (photograph on 8-10 days April
Piece)>0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI (May 13-15
The photo of day) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32.
S4: using n × n as pixel window, moving step length 1, threshold value is that adjoining piecemeal is aggregated into bulk by m, will be by
The fritter for being unsatisfactory for condition surrounded in bulk is filled to form continuous piecemeal inside regular edges, and exports image,
Middle n=3-5, m=4-6 can be obtained the area of various situations according to image.
Embodiment 3
It is lain fallow the device of situation the present embodiment provides a kind of by remote sensing satellite photo statistical monitoring crop rotation, comprising:
Including image capture module, image processing module, classification demarcating module and result output module；
It is distant to the annual satellite of second year by the end of June by the end of June that described image acquisition module is used to obtain areal current year
Feel photo；
Described image processing module is used to obtain NDVI value, Green value and the Blue of each pixel in each photo
Value, or piecemeal first is carried out to the arable land plot in satellite remote sensing photo and obtains the NDVI of each piecemeal in each photo
The average value of value, Green value and Blue value,
Wherein, NDVI=(ρNIR-ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB),
Blue=ρG/(ρR+ρG+ρB), ρNIR、ρR、ρG、ρBRespectively indicate near infrared band reflectivity, red wave band reflectivity,
Green wave band reflectivity, blue wave band reflectivity；
The classification demarcating module is for demarcating the pixel in satellite remote sensing photo, the judgment basis base of calibration
In the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
It ploughs rape and green manuring is lain fallow for crop rotation；
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI (August 13-
Photo on the 15th)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24,
NDVI (photo on 10-12 days November) < 0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)
>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI (May 13-15
The photo of day) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photograph on 8-10 days April
Piece)>0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI (May 13-15
Day photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (May
3-5 days photos) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March) < 0.40, NDVI (photograph on 8-10 days April
Piece)>0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI (May 13-15
The photo of day) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32；
The result output module is used to export the image demarcated according to the result of the classification demarcating module.
Preferably, in the classification demarcating module, eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition
Plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow.
Preferably, described image processing module is used to obtain NDVI value, the Green value of each pixel in each photo
When with Blue value, in the result output module, using n × n as pixel window, moving step length 1, threshold value is that m will be adjoining
Piecemeal aggregates into bulk, will be filled to be formed inside regular edges continuously by the fritter for being unsatisfactory for condition surrounded in bulk
Piecemeal, and export image.
Preferably, the n=3 or 4 or 5, m=4 or 5 or 6.
Preferably, piecemeal is carried out by the ridge road to arable land using more to the arable land plot in satellite remote sensing photo
Multi-scale segmentation method identifies to obtain.
It is illustrated by taking the remote sensing image in Suzhou districts under city administration concerning farmers as an example below, including Wujiang area, Wuzhong District, Huqiu District
(Fig. 1, elevation use SRTM 90m resolution ratio DEM data, and data source is believed in Chinese Academy of Sciences's computer network with Xiangcheng District
Breath center international scientific data image website, http://www.gscloud.cn).Research area is located in South of Jiangsu Province, Taihu Lake edge
Bank belongs to subtropical zone monsoon marine climate, and topography is low flat, mean annual precipitation 1099.6mm, abundant rainfall, rice cultivation,
Wheat, rape and green manure crop produce silkworm and mulberry, woods fruit and tealeaves；The domestic network of waterways is intensive, causes agriculture plantation plot broken.Suzhou
Crop rotation is implemented to wheat-rice rotation area arable land in 2017 and is lain fallow in city.
Remote sensing image data uses European Space Agency (ESA) sentry 2 (Sentinel-2) data.Sentinel-2 satellite includes
Sentinel-2A satellite and Sentinel-2B satellite, multispectral data include 13 wave bands, and providing spatial resolution is respectively
The data of 10m, 20m, 60m, wherein 10m resolution ratio wave band include blue (B2:458~523mm), green (B3:543~578mm),
Red (B4:650~680mm), 4 wave bands of near-infrared (B8:785~900mm).Sentinel-2A satellite emitted in 2015,
Sentinel-2B satellite emitted in 2017, and double star networking can revisit for 5 days.From European Space Agency's dissemination system (https: //
Scihub.copernicus.eu/dhus/#/home the 10 scape images in covering research area, in June, 2017-2018 time) are obtained
It is exported and is closed using the wave band that the Radiance Sentinel-2LIC module of ENVI5.3 carries out Sentinel-2 data June in year
At and radiation calibration, atmospheric correction carried out under FLAASH module.
For phenomena such as research area plot is crushed, pattern of farming is complicated, to guarantee to extract the integrality of ground block message and dividing
Class precision, using GF-2 0.8m resolution image based on the convolution algorithm of machine learning to arable land ground block message identify,
The ridge road in arable land is identified using multi-scale segmentation method, and combines high-resolution Google Earth image to recognition result
Quality inspection and amendment are carried out, 44082, plot is drawn in production altogether.
Above-mentioned pixel result is mapped on the arable land plot produced using the method for embodiment 1.It is superimposed pixel grid
Trrellis diagram layer and arable land plot VectorLayer, using Spatial Join and the Intersect tool of ArcGIS, automatically by overlapping
Arable land plot marks, again assignment objective attribute target attribute.Using verification mode on area and spatial position: i.e. based on statistical data and
Precision test based on spatial distribution.The verifying of area is lain fallow each mode area and statistical report with the crop rotation of remote sensing monitoring
Data subregion carries out evaluation verifying, and generally lie fallow error of sunning the earth which has been ploughed up of remote sensing monitoring is 0.22%, remote sensing monitoring crop rotation error
It is -0.19%；The verifying of spatial distribution is layouted on the spot with the higher Google image of resolution ratio and field operation, chooses rice cultivation
Plot, stubble plantation ploughs rape, does not plough the plot of rape, wheat, green manure after rice, and lie fallow sun the earth which has been ploughed up plot and other
The plot of situation, whether block type and verifying sampling point are consistent with verifying remote sensing monitoring, 90% or more overall matching precision.
It is enlightenment with the above-mentioned desirable embodiment according to the application, through the above description, relevant staff is complete
Full various changes and amendments can be carried out in the range of without departing from this item application technical idea.The technology of this item application
Property range is not limited to the contents of the specification, it is necessary to which the technical scope thereof is determined according to the scope of the claim.
Claims (10)
1. a kind of lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, which comprises the following steps:
S1: areal current year 1 year satellite remote sensing photo to second year by the end of June by the end of June is obtained, to satellite remote sensing photo
In arable land plot carry out piecemeal or without piecemeal；
S2: the NDVI value of each piecemeal in each photo, the average value of Green value and Blue value, NDVI=(ρ are obtainedNIR-
ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρR+ρG+ρB), ρNIR、ρR、ρG、ρBRespectively indicate near-infrared
Wave band reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；When without piecemeal, obtain in each photo
Each pixel NDVI value, Green value and Blue value；
S3: the judgement of NDVI value, Green value and Blue value is carried out to each piecemeal in satellite remote sensing photo or pixel, and right
Piecemeal or pixel are demarcated, and judgment basis is based on the following conditions:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
It ploughs rape and green manuring is lain fallow for crop rotation；
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI be (August 13-15 days
Photo)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24, NDVI (November
10-12 days photos)<0.31, NDVI (29-31 days photos of August)-NDVI (photo on 10-12 days November)>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photo on 8-10 days April) >
0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI be (13-15 days May
Photo) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photo on 8-10 days April) >
0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI be (13-15 days May
Photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (3-5 days May
Photo) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March)<0.40, NDVI (photo on 8-10 days April)>
0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI be (13-15 days May
Photo) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32；
S4: summarize piecemeal or pixel as a result, and exporting image.
2. according to claim 1 lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, feature exists
In, in the S3 step, eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow.
3. according to claim 1 or 2 lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, feature
It is, when to satellite remote sensing photo without piecemeal, in S4 step, using n × n as pixel window, moving step length 1, threshold value
Adjoining piecemeal is aggregated into bulk for m, will be filled to form edge by the fritter for being unsatisfactory for condition surrounded in bulk
The continuous piecemeal of regular interior, and export image.
4. according to claim 3 lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation, feature exists
In the n=3 or 4 or 5, m=4 or 5 or 6.
5. according to claim 1-4 lain fallow the method for situation by remote sensing satellite photo statistical monitoring crop rotation,
It is characterized in that, carrying out piecemeal to the arable land plot in satellite remote sensing photo uses multiple dimensioned point by the ridge road to arable land
Segmentation method identifies to obtain.
6. a kind of lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, which is characterized in that including Image Acquisition
Module, image processing module, classification demarcating module and result output module；
Described image acquisition module is used to obtain annual satellite remote sensing photograph of the areal current year by the end of June to second year by the end of June
Piece；
Described image processing module is used to obtain NDVI value, Green value and the Blue value of each pixel in each photo, or
Piecemeal first is carried out to the arable land plot in satellite remote sensing photo and obtains the NDVI value of each piecemeal in each photo, Green
The average value of value and Blue value,
Wherein, NDVI=(ρNIR-ρR)/(ρNIR+ρR), Green=ρG/(ρR+ρG+ρB), Blue=ρG/(ρR+ρG+ρB), ρNIR、ρR、
ρG、ρBRespectively indicate near infrared band reflectivity, red wave band reflectivity, green wave band reflectivity, blue wave band reflectivity；
For demarcating to the pixel in satellite remote sensing photo, the judgment basis of calibration is based on following the classification demarcating module
Condition:
Eligible 1 and condition 5 and it is ineligible 2 when sun the earth which has been ploughed up to lie fallow；
Eligible 1, condition 2 and condition 4 and it is ineligible 3 when to plough rape；
Eligible 1, condition 6 and be green manuring when ineligible 2 and 5；
It ploughs rape and green manuring is lain fallow for crop rotation；
Condition 1, while meeting the following conditions: 0.55 < NDVI (photo on 27-30 days July) < 0.75, NDVI be (August 13-15 days
Photo)>0.36, NDVI (29-31 days photos of August)>0.42, NDVI (29-31 days photos of September)<0.24, NDVI (November
10-12 days photos)<0.31, NDVI (29-31 days photos of August)-NDVI (photo in 10-12 days 11 months November)>0.30；
Condition 2, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.28, NDVI (photo on 8-10 days April) >
0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.72, NDVI be (13-15 days May
Photo) < 0.45, NDVI (photo on 12-14 days June) < 0.23；
Condition 3, while meeting the following conditions: Green (photo on 8-10 days April) > Blue (photo on 8-10 days April)；
Condition 4, while meeting the following conditions: NDVI (photo on 9-11 days March) > 0.37, NDVI (photo on 18-20 days April) >
0.47, NDVI (photo on 18-20 days April)>0.24, NDVI (photo on 3-5 days May)<0.16, NDVI be (13-15 days May
Photo) < 0.14, NDVI (photo on 12-14 days June) < 0.22, NDVI (photo on 13-15 days May)-NDVI (3-5 days May
Photo) < 0.2；
Condition 5, while meeting the following conditions: NDVI (photo on 9-11 days March)<0.40, NDVI (photo on 8-10 days April)>
0.05, NDVI (photo on 18-20 days April)<0.40, NDVI (photo on 3-5 days May)>0.05, NDVI be (13-15 days May
Photo) < 0.40, NDVI (photo on 12-14 days June) < 0.26；
Condition 6, ρNIR(photo on 8-10 days April) > 0.32；
The result output module is used to export the image demarcated according to the result of the classification demarcating module.
7. according to claim 6 lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, feature exists
In, in the classification demarcating module, eligible 1 and condition 2 and ineligible 3 and 4 Shi Weiwei of condition plough rape；
Eligible 1, condition 2 and condition 3 are to have planted the plot of wheat；
Ineligible 1 is non-aqueous rice field block；
Not ploughing rape or having planted wheat is that non-crop rotation is lain fallow.
8. according to claim 6 or 7 lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, feature
It is, when described image processing module is used to obtain the NDVI value, Green value and Blue value of each pixel in each photo,
In the result output module, using n × n as pixel window, moving step length 1, threshold value is that m aggregates into adjoining piecemeal greatly
Block will be filled to form continuous piecemeal inside regular edges by the fritter for being unsatisfactory for condition surrounded in bulk, and export
Image.
9. according to claim 8 lain fallow the device of situation by remote sensing satellite photo statistical monitoring crop rotation, feature exists
In the n=3 or 4 or 5, m=4 or 5 or 6.
10. it is lain fallow the device of situation according to claim 6-9 is described in any item by remote sensing satellite photo statistical monitoring crop rotation,
It is characterized in that, carrying out piecemeal to the arable land plot in satellite remote sensing photo uses multiple dimensioned point by the ridge road to arable land
Segmentation method identifies to obtain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811363119.XA CN109471131B (en) | 2018-11-15 | 2018-11-15 | Method and device for counting and monitoring crop rotation fallow situation through remote sensing satellite photo |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811363119.XA CN109471131B (en) | 2018-11-15 | 2018-11-15 | Method and device for counting and monitoring crop rotation fallow situation through remote sensing satellite photo |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109471131A true CN109471131A (en) | 2019-03-15 |
| CN109471131B CN109471131B (en) | 2022-05-24 |
Family
ID=65673774
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811363119.XA Active CN109471131B (en) | 2018-11-15 | 2018-11-15 | Method and device for counting and monitoring crop rotation fallow situation through remote sensing satellite photo |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109471131B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110909679A (en) * | 2019-11-22 | 2020-03-24 | 中国气象科学研究院 | Remote sensing identification method and system for fallow crop rotation information of winter wheat historical planting area |
| CN111007013A (en) * | 2019-11-01 | 2020-04-14 | 中科禾信遥感科技（苏州）有限公司 | Crop rotation fallow remote sensing monitoring method and device for northeast cold region |
| CN111444783A (en) * | 2020-03-11 | 2020-07-24 | 中科禾信遥感科技（苏州）有限公司 | Crop planting land parcel identification method and device based on pixel statistics |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106845806A (en) * | 2017-01-06 | 2017-06-13 | 中国科学院遥感与数字地球研究所 | The remote-sensing monitoring method and system of farmland planting state |
| CN106918816A (en) * | 2017-03-21 | 2017-07-04 | 中国科学院遥感与数字地球研究所 | crop growth monitoring system and its data processing method and module |
| CN107064932A (en) * | 2017-02-28 | 2017-08-18 | 华南农业大学 | A kind of construction land change detection method based on time series SAR remote sensing images |
| WO2018081043A1 (en) * | 2016-10-24 | 2018-05-03 | Board Of Trustees Of Michigan State University | Methods for mapping temporal and spatial stability and sustainability of a cropping system |
-
2018
- 2018-11-15 CN CN201811363119.XA patent/CN109471131B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018081043A1 (en) * | 2016-10-24 | 2018-05-03 | Board Of Trustees Of Michigan State University | Methods for mapping temporal and spatial stability and sustainability of a cropping system |
| CN106845806A (en) * | 2017-01-06 | 2017-06-13 | 中国科学院遥感与数字地球研究所 | The remote-sensing monitoring method and system of farmland planting state |
| CN107064932A (en) * | 2017-02-28 | 2017-08-18 | 华南农业大学 | A kind of construction land change detection method based on time series SAR remote sensing images |
| CN106918816A (en) * | 2017-03-21 | 2017-07-04 | 中国科学院遥感与数字地球研究所 | crop growth monitoring system and its data processing method and module |
Non-Patent Citations (1)
| Title |
|---|
| ANDREW D 等: "Tracking vegetation phenology across diverse North American biomes using PhenoCam imagery", 《SCIENTIFIC DATE》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111007013A (en) * | 2019-11-01 | 2020-04-14 | 中科禾信遥感科技（苏州）有限公司 | Crop rotation fallow remote sensing monitoring method and device for northeast cold region |
| CN111007013B (en) * | 2019-11-01 | 2022-10-14 | 中科禾信遥感科技（苏州）有限公司 | Crop rotation fallow remote sensing monitoring method and device for northeast cold region |
| CN110909679A (en) * | 2019-11-22 | 2020-03-24 | 中国气象科学研究院 | Remote sensing identification method and system for fallow crop rotation information of winter wheat historical planting area |
| CN110909679B (en) * | 2019-11-22 | 2020-09-18 | 中国气象科学研究院 | Remote sensing identification method and system for fallow crop rotation information of winter wheat historical planting area |
| CN111444783A (en) * | 2020-03-11 | 2020-07-24 | 中科禾信遥感科技（苏州）有限公司 | Crop planting land parcel identification method and device based on pixel statistics |
| CN111444783B (en) * | 2020-03-11 | 2023-11-28 | 中科禾信遥感科技（苏州）有限公司 | Crop planting land block identification method and device based on pixel statistics |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109471131B (en) | 2022-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhou et al. | Long-short-term-memory-based crop classification using high-resolution optical images and multi-temporal SAR data | |
| Azar et al. | Assessing in-season crop classification performance using satellite data: A test case in Northern Italy | |
| US8712148B2 (en) | Generating agricultural information products using remote sensing | |
| Long et al. | Object-oriented crop classification using multitemporal ETM+ SLC-off imagery and random forest | |
| Esquerdo et al. | Use of NDVI/AVHRR time-series profiles for soybean crop monitoring in Brazil | |
| CN104915674B (en) | The method that Landsat8 and MODIS constructs high-spatial and temporal resolution data identification autumn grain crop | |
| CN106918816B (en) | Crop growth monitoring system and its data processing method and module | |
| Vaudour et al. | Early-season mapping of crops and cultural operations using very high spatial resolution Pléiades images | |
| CN109471131A (en) | It is lain fallow the method and apparatus of situation by remote sensing satellite photo statistical monitoring crop rotation | |
| CN107527364A (en) | A kind of seaweed growing area monitoring method based on remote sensing images and lace curtaining information | |
| CN113033279A (en) | Crop fine classification method and system based on multi-source remote sensing image | |
| CN104951754A (en) | Sophisticated crop classifying method based on combination of object oriented technology and NDVI (normalized difference vegetation index) time series | |
| Hebbar et al. | Object oriented classification of high resolution data for inventory of horticultural crops | |
| Xie et al. | Suitable remote sensing method and data for mapping and measuring active crop fields | |
| Song et al. | Recognition of sunflower growth period based on deep learning from UAV remote sensing images | |
| Demir et al. | Determination of opium poppy (Papaver somniferum) parcels using high-resolution satellite imagery | |
| Li et al. | Potential Evaluation of High Spatial Resolution Multi-Spectral Images Based on Unmanned Aerial Vehicle in Accurate Recognition of Crop Types | |
| Zou et al. | The fusion of satellite and unmanned aerial vehicle (UAV) imagery for improving classification performance | |
| Cui et al. | Optimal spatial resolution of remote-sensing imagery for monitoring cantaloupe greenhouses | |
| Szatmári et al. | The Role of Field Survey in the Identification of Farmland Abandonment in Slovakia Using Sentinel-2 Data | |
| Yang et al. | Extraction of multiple cropping information at the Sub-pixel scale based on phenology and MODIS NDVI time-series: a case study in Henan Province, China | |
| Shan et al. | Study on extraction methods of paddy rice area based on GF-6 satellite image | |
| CN109214287A (en) | Crops decomposition method and system based on RapidEye satellite remote-sensing image | |
| Sun et al. | A spatiotemporal collaborative approach for precise crop planting structure mapping based on multi-source remote-sensing data | |
| Gorokhova et al. | Application of geoinformational technologies and aerospace photography materials for monitoring the land irrigated by the Svetlyi Yar irrigation system (Volgograd oblast) |
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 |