CN112001809A - Method for acquiring farmland returning information of agriculture and forestry area - Google Patents
Method for acquiring farmland returning information of agriculture and forestry area Download PDFInfo
- Publication number
- CN112001809A CN112001809A CN202010761233.9A CN202010761233A CN112001809A CN 112001809 A CN112001809 A CN 112001809A CN 202010761233 A CN202010761233 A CN 202010761233A CN 112001809 A CN112001809 A CN 112001809A
- Authority
- CN
- China
- Prior art keywords
- year
- ndvi
- information
- value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000009336 multiple cropping Methods 0.000 claims abstract description 48
- 238000012544 monitoring process Methods 0.000 claims abstract description 24
- 238000001914 filtration Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 8
- 238000010899 nucleation Methods 0.000 claims description 7
- 238000003971 tillage Methods 0.000 abstract description 6
- 238000000605 extraction Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 10
- 244000025254 Cannabis sativa Species 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 241000498849 Chlamydiales Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 241000209504 Poaceae Species 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Marine Sciences & Fisheries (AREA)
- Human Resources & Organizations (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Multimedia (AREA)
- Economics (AREA)
- Agronomy & Crop Science (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Image Processing (AREA)
Abstract
The invention provides a method for acquiring farmland returning information in an agricultural and forestry area, which comprises the following steps: obtaining an NDVI sequence value within one year according to the remote sensing image with low space and high time resolution; calculating a CV value of any year and extracting crop phenological information according to the NDVI sequence value, and further determining planting information of any pixel area; determining a multiple cropping index of any pixel area within the monitoring years according to the planting information and the crop phenological information; based on the multiple cropping index, the remote sensing image sequence with medium and high spatial resolution is combined to accurately determine the tillage withdrawal time, the tillage withdrawal area and the tillage withdrawal position. The invention utilizes the agricultural and forestry area remote sensing images with low space and high time resolution to obtain the NDVI value of a time sequence, calculates the land retreating position, the land retreating time and the land retreating area, then applies the remote sensing images with medium and high resolution and combines different multiple cropping index remote sensing models to obtain the land retreating information, thereby improving the accuracy of information extraction.
Description
Technical Field
The invention relates to the technical field of agriculture and forestry protection, in particular to a method for acquiring farmland returning information in an agriculture and forestry area.
Background
For the requirements of protecting some ecological areas with fragile ecological environment or important ecological functional areas, some farmlands need to be used as return lands for returning forest and grass to realize the maximization of the value of the ecological system.
However, from the management perspective, specific spatial positions and areas of returning ploughs at different times need to be counted in time, certain economic compensation is given to places or farmers according to the time and the areas of returning ploughs, and meanwhile, important support is provided for monitoring the effect of returning the returning farmlands to the grasses.
Due to the characteristics of multiple times, multiple spaces, multiple spectrums and the like, the remote sensing technology can effectively acquire relevant information of land surface coverage and change, and therefore, the remote sensing technology has irreplaceable potential for identification of land withdrawal and time withdrawal in agricultural and forestry grass areas.
However, according to the results of the related remote sensing research and application at present, the remote sensing technology is mostly used for extracting information such as land withdrawal and the like by classifying land utilization through spectrum information or combining texture information and the like, and then obtaining land information. The problem of doing so is that the types of crops are different, the planting time difference may be large, and the classification method is mostly obtained by adopting one time phase, so that the possibility that some farmland information cannot be effectively extracted may exist; secondly, land is frequently abandoned, and the characteristics that the spatial distribution is dispersed, the land size is small, and the mixed situation of farmland, woodland, grassland, bare soil, water body and the like is possible, so that the extracted result has great uncertainty; thirdly, the land use classification method has large calculation amount and high calculation resource consumption, and besides, due to the selection of training samples and other reasons, the land use classification results in each time have certain difference.
Based on the problems, the accuracy of the traditional method for extracting the fallow land by remote sensing has certain uncertainty.
Disclosure of Invention
The embodiment of the invention provides a method for acquiring land returning information of an agricultural and forestry area, which improves the accuracy of information extraction of land returning information.
The embodiment of the invention provides a method for acquiring farmland returning information in an agricultural and forestry area, which comprises the following steps:
taking a year as a unit, extracting an NDVI value of each pixel region in a remote sensing image sequence according to the remote sensing image sequence with low space and high time resolution acquired in one year to obtain an NDVI sequence value corresponding to the remote sensing image sequence;
calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year;
determining planting information of any one of the image element areas according to the CV value of any one year and the NDVI value of any one of the image element areas of any one year;
determining a multiple cropping index of any pixel area in the monitored years according to the planting information of any pixel area in each year and the crop phenological information of any year in the monitored years;
and determining the back-plowing time, the back-plowing ground area and the back-plowing area position based on the multiple cropping index of each pixel region in the monitoring years and by combining the remote sensing image sequence with the medium and high spatial resolution acquired in each year.
On the basis of the above technical solutions, the embodiments of the present invention may be further improved as follows.
Optionally, taking the year as a unit, extracting an NDVI value of each pixel region in the remote sensing image sequence according to the remote sensing image sequence with low space and high time resolution acquired within one year, and obtaining an NDVI sequence value corresponding to the remote sensing image sequence includes:
extracting the NDVI value of each pixel region according to the remote sensing image sequence with low space and high time resolution in one year to obtain the NDVI sequence value corresponding to the remote sensing image sequence;
carrying out filtering processing on the NDVI sequence value to obtain a filtered NDVI sequence value;
correspondingly, calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year comprises the following steps:
and calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value after filtering processing of any year.
Optionally, the calculating the CV value of any year according to the NDVI sequence value after filtering processing of any year includes:
wherein NDVIiIs the NDVI value, CV, of the i-th pixel area of the low-resolution remote sensing image sequence of the j-th yearjIs the coefficient of variation, σ, of the NDVI value of the pixel region of the j-th low-resolution remote-sensing image sequencejIs the variance, mu, of the NDVI value of the pixel region of the j-th remote sensing image sequencejThe average value of the NDVI values of the pixel areas of the j-th remote sensing image sequence is shown.
Optionally, the determining the planting information of any one of the image element areas according to the CV value of any one year and the NDVI value of any one of the image element areas of any one year includes:
if the ith pixel area in the jth year meets the NDVIiNot less than b and CVjAnd a is more than or equal to a, the crop planting information of the ith pixel area is planting, wherein a and b are preset threshold values.
Optionally, the multiple cropping index of any pixel area in the monitored years is determined according to the planting information of any pixel area in each year and the crop phenological information of any year in the monitored years:
preliminarily determining a multiple cropping index of any pixel area in the monitored years based on the planting information of any pixel area in the remote sensing image sequence in each year;
and accurately determining the multiple cropping index of any pixel area in the monitoring years based on the crop phenological information of any year and the preliminarily determined multiple cropping index of any pixel area in the monitoring years.
Optionally, the preliminarily determining the multiple cropping index of any pixel area within the monitored years based on the planting information of any pixel area in the remote sensing image sequence in each year includes:
if the planting information of any pixel area in each year in the monitored years is planting, the multiple planting index of any pixel area is more than 1;
if the planting information of one part of the year is planting and the planting information of the other part of the year is not planting in any pixel area within the monitoring years, the multiple planting index of any pixel area is between 0 and 1.
Optionally, the accurately determining the multiple cropping index of any one of the pixel areas in the monitored years based on the crop phenological information of any year and the preliminarily determined multiple cropping index of any one of the pixel areas in the monitored years includes:
based on the remote sensing image sequence with low space and high time resolution, extracting an NDVI sequence value corresponding to the remote sensing image sequence of each year;
generating a crop phenological curve of each year according to the NDVI sequence value;
determining a plurality of NDVI values of a plurality of growth starting points of the crops and a plurality of NDVI values of a high-value area in a growth vigorous period according to the crop climatic curve;
and determining the multiple cropping index of each remote sensing pixel area within the monitoring years to be 2 or 1 or between 0 and 1 according to the NDVI values of the multiple growth starting points of the crops and the NDVI values of the high-value areas in the vigorous growth period.
Optionally, the determining of the plowing-down time, the plowing-down ground area and the position of the plowing-down position based on the multiple cropping index of each pixel region in the number of monitored years by combining the remote sensing image sequence with the medium and high spatial resolution acquired every year comprises:
for the pixel area with multiple seeding index of 1 or 2, if any pixel area is the t th pixel areajPlanting information of year is planting from the tj+1If the planting information of the beginning of the year is no planting, the pixel area is determined as withdrawn farmland, and the time of withdrawing the farmland is tthjYear;
for the pixel area with multiple seeding index of 0-1, if any pixel area is the t th pixel areajPlanting information of year is planting, and within d years later, if planting information is not planting, the time of back plowing is tthjYear;
calculating the land retreating area by the following formula:
S1＝e.p2；
S1the area of the abandoned land is shown as e, the number of pixel areas of the abandoned land is shown as p, and the pixel spatial resolution of the remote sensing image is shown as p.
The method for acquiring the land retreating information of the agricultural and forestry area, the electronic device and the storage medium provided by the embodiment of the invention are combined with a remote sensing identification technology, utilize the agricultural and forestry area remote sensing image with low space and high time resolution to extract the corresponding NDVI value, automatically calculate the position of the land retreated, the time of the land retreated and the area of the land retreated to acquire, and utilize the automatic remote sensing identification technology to improve the accuracy of the information extraction of the land retreated.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for acquiring information of returning plowing in an agricultural and forestry area according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, a method for acquiring land withdrawal information of an agricultural and forestry area is provided, which includes:
s1, taking a year as a unit, extracting the NDVI value of each pixel region in the remote sensing image sequence according to the remote sensing image sequence with low space and high time resolution acquired in one year to obtain the NDVI sequence value corresponding to the remote sensing image sequence;
s2, calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year;
s3, determining planting information of any one of the pixel areas according to the CV value of any one year and the NDVI value of any one of the pixel areas of any one year;
s4, determining a multiple cropping index of any pixel area in the monitored years according to the planting information of any pixel area in each year and the crop phenological information of any year in the monitored years;
and S5, determining the back-plowing time, the back-plowing ground area and the back-plowing area position based on the multiple cropping index of each pixel region in the number of monitored years and by combining the remote sensing image sequence with the medium and high spatial resolution acquired in each year.
It can be understood that the remote sensing image data of the agriculture and forestry area is shot and collected by using a remote sensing technology, and when shooting is carried out, the remote sensing image data of the agriculture and forestry area can be shot by adopting low space high time resolution and medium high space resolution, so that the remote sensing image data of the low space high time resolution and the remote sensing image data of the medium and high space resolution are respectively obtained. The density of the remote sensing image data with low space and high time resolution is relatively high, and the density of the remote sensing image data with medium and high space resolution is relatively low.
When the remote sensing technique is used for shooting, regular shooting may be performed at the same time every year, for example, every 8 days every year, and the time point of each shooting is the same for different years, that is, the same time point of each year.
And acquiring a remote sensing image sequence for each year, wherein the remote sensing image sequence consists of a series of remote sensing pixels. According to the embodiment of the invention, based on the remote sensing image sequence data with low space and high time resolution, for each remote sensing pixel in the remote sensing image sequence data, a corresponding NDVI (Normalized Difference Vegetation index) value is extracted, and the NDVI sequence value corresponding to the remote sensing image sequence of each year can be obtained.
For any year, calculating a CV value of the year according to the NDVI sequence value of the year, wherein the CV value is the coefficient of variation of the NDVI sequence value corresponding to the remote sensing image sequence of any year; and extracting the phenological information of the crop in any year according to the NDVI sequence value of the year. And determining planting information of each pixel area, namely whether each pixel area is in a planting state or a non-planting state according to the CV value of the year and the NDVI value of each pixel area of the year. Determining a multiple cropping index of any pixel area in the monitored years according to the planting information and the crop phenological information of any pixel area in each year in the monitored years; and determining the back-plowing time, the back-plowing ground area and the back-plowing area position according to the multiple cropping index of each pixel area within the monitoring years and the remote sensing image sequence with medium and high spatial resolution.
The embodiment of the invention combines a remote sensing identification technology, utilizes the remote sensing image of the agriculture and forestry area with low space and high time resolution to extract the corresponding NDVI value, automatically calculates the position of the abandoned land, the abandoned land time and the area of the abandoned land, and utilizes the automatic remote sensing identification technology to improve the extraction precision of the abandoned land information.
As a possible implementation manner, taking a year as a unit, extracting an NDVI value of each pixel region in a remote sensing image sequence according to a remote sensing image sequence with low space and high time resolution acquired within one year, and obtaining an NDVI sequence value corresponding to the remote sensing image sequence includes:
extracting the NDVI value of each pixel region according to the remote sensing image sequence with low space and high time resolution in one year to obtain the NDVI sequence value corresponding to the remote sensing image sequence;
carrying out filtering processing on the NDVI sequence value to obtain a filtered NDVI sequence value;
correspondingly, calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year comprises the following steps:
and calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value after filtering processing of any year.
It can be understood that remote sensing data with high time resolution such as MODIS are selected, time-series NDVI is extracted, and a time-series crop growth period growth curve is obtained through a filtering method. The filtering selects a Savitzky-Golay filtering model. Savitzky-Golay filter primitive:
where Y is the original value of NDVI, Y is the fitted value of NDVI, CiIs the coefficient when filtering the ith NDVI value, N is the number of convolutions and is also equal to the width of the sliding array (2m + 1). The coefficient j refers to the coefficient of the original NDVI array. The sliding array contains (2m +1) points. The method is essentially a kind of smooth filtering, so two parameters control the filtering effect, one is m, i.e. the size of the filtering window; the second is the degree of the smoothing polynomial.
The NDVI sequence values of each year are filtered by adopting the filtering method to obtain the filtered NDVI sequence values of each year, the CV value corresponding to each year is calculated based on the filtered NDVI sequence values of each year, and the phenological information of the crops is extracted based on the filtered NDVI sequence values.
As a possible implementation, calculating the CV value of any year according to the filtered NDVI sequence value of any year includes:
wherein NDVIiIs the NDVI value, CV, of the i-th pixel area of the low-resolution remote sensing image sequence of the j-th yearjIs the coefficient of variation, σ, of the NDVI value of the pixel region of the j-th low-resolution remote-sensing image sequencejIs the variance, mu, of the NDVI value of the pixel region of the j-th remote sensing image sequencejThe average value of the NDVI values of the pixel areas of the j-th remote sensing image sequence is shown.
It can be understood that, for the NDVI sequence values corresponding to the remote sensing image sequence with low resolution and high time resolution, the CV value of each year is calculated according to the above formula.
And determining planting information of each pixel area in the remote sensing image sequence according to the NDVI sequence value of the remote sensing image sequence of each year.
Wherein, NDVI for the ith pixel area of the j yearsiNot less than b and CVjAnd a is larger than or equal to a, the planting information of the ith pixel area is planting, wherein a and b are preset threshold values. That is, taking the jth year as an example, wherein the NDVI value of the ith pixel area is greater than or equal to b, and the CV value of the jth year is greater than or equal to a, then the planting information of the ith pixel area is planting, that is, the pixel area is a farmland and crops are planted; otherwiseThe image element area is other vegetation.
As a possible embodiment, the multiple cropping index of any one pixel area in the monitored years is determined according to the planting information of any one pixel area in each year and the crop phenological information of any one year in the monitored years:
preliminarily determining a multiple cropping index of any pixel area in the monitored years based on the planting information of any pixel area in the remote sensing image sequence in each year;
and accurately determining the multiple cropping index of any pixel area in the monitoring years based on the crop phenological information of any year and the preliminarily determined multiple cropping index of any pixel area in the monitoring years.
It can be understood that planting information of each pixel area can be determined based on the low-space high-time-resolution remote sensing image sequence of each year. And preliminarily determining the multiple cropping index of each pixel area within the monitoring years based on the planting information of each pixel area determined by the remote sensing image sequence with low space and high time resolution.
And then accurately determining the multiple cropping index of each pixel area in the monitoring years based on the preliminarily determined multiple cropping index of each pixel area in the monitoring years and the crop phenological information.
As a possible implementation, the preliminary determination of the cropping index of any one of the image element regions in the monitored years based on the planting information of any one of the image element regions in the low-space high-time-resolution remote sensing image sequence in each year includes:
if the planting information of any remote sensing pixel area in each year in the monitored years is planting, the multiple planting index of any remote sensing pixel area is more than 1;
if the planting information of one part of the year is planting and the planting information of the other part of the year is not planting in any remote sensing pixel area within the monitored year, the multiple planting index of any remote sensing pixel area is between 0 and 1.
It can be understood that the multiple information of each pixel region is preliminarily determined based on the remote sensing image sequence with low resolution and high time resolution as follows: in the monitoring years, crops are planted in the pixel area at the same position in a set number of years (monitoring years), namely planting information of the pixel area at the same position in each year is planted, and then the multiple cropping index of the pixel area is determined to be greater than 1.
In the monitoring years, the pixel area at the same position is in a set year (monitoring years), planting is carried out in one part of the year, and no planting is carried out in the other part of the year, so that the multiple cropping index of the pixel area is determined to be between 0 and 1.
Vectorizing the image element region with the multiple seeding index larger than 1 and between 0 and 1.
As a possible embodiment, the accurately determining the multiple cropping index of any one picture element region in the monitored years based on the crop phenological information in any year and the preliminarily determined multiple cropping index of any one picture element region in the monitored years comprises:
based on the remote sensing image sequence with low space and high time resolution, extracting an NDVI sequence value corresponding to the remote sensing image sequence of each year;
generating a crop phenological curve of each year according to the NDVI sequence value;
determining a plurality of NDVI values of a plurality of growth starting points of the crops and a plurality of NDVI values of a high-value area in a growth vigorous period according to the crop climatic curve;
and determining the multiple cropping index of each remote sensing pixel area within the monitoring years to be 2 or 1 or between 0 and 1 according to the NDVI values of the multiple growth starting points of the crops and the NDVI values of the high-value areas in the vigorous growth period.
It can be understood that, when the multiple indexes of each pixel region are accurately determined, for the remote sensing image sequence with low space and high time resolution, the NDVI sequence value corresponding to the remote sensing image sequence of each year is extracted.
Connecting each NDVI value in the NDVI sequence values to generate a crop phenological curve, analyzing the phenological curves in different growth periods by combining the phenological information of the crops recorded in the local for many years, and extracting the characteristics of time points of the curves in different stages of crop planting, ridge sealing, jointing and the like.
And determining the NDVI value corresponding to the growth starting point and the NDVI value corresponding to the growth vigorous point of the crops for different crops according to the crop climatic curve.
Determining the growth start point NDVI value Vg of different kinds of crop phenological curves assuming that the area is planted with at most three crops (such as winter wheat and corn)m、Vgn、VgpAnd the NDVI value Vf of the high-value region with vigorous growthm、Vfn、VfpWherein, Vgm、VgnAnd VgpExpressing the NDVI value, Vf value of the starting point of growth of each crop in one yearm、VfnAnd VfpNDVI values for high value areas where the plants grow vigorously per season.
According to the NDVI values of the crops, the multiple cropping index of each image element area can be further determined, and according to the result of the preliminary determination, whether the multiple cropping index of each image element area is 2 or 1 or is between 0 and 1 is further determined.
The specific judgment process is that the following conditions are set firstly:
condition 1: simultaneously, the following conditions are met:
Vfm-Vgm≥a1，Vfm-Vgn≥a2，2Vfm-Vgm-Vgn≥a3；
condition 2: simultaneously, the following conditions are met:
Vfn-Vgn≥b1；Vfp-Vgp≥b2；2Vfn-Vgp-Vgn≥b3；
condition 3: simultaneously, the following conditions are met:
Vfp-Vgp≥c1；Vfn-Vgp≥c2；2Vfp-Vgp-Vgn≥c3。
wherein, a1、a2、a3,b1、b2、b3,c1、c2、c3Is a preset threshold.
If the pixel area of a certain spatial position in each year in a set period (within the monitoring years) simultaneously meets the conditions 1,2 and 3 or simultaneously meets the conditions 1 and 2, judging that the multiple cropping index of the pixel area is 2; if only one of the conditions 1,2 or 3 is satisfied, the complex seed index of the chlamydial region is judged to be 1; and if the number of years only meets 1 or 2 or 3 and the number of years does not meet one year, judging that the multiple cropping index of the pixel region is between 0 and 1.
As a possible implementation, the determining of the plowing time, the plowing ground area and the position of the plowing location based on the multiple cropping index of each pixel area in the number of monitored years and the combination of the remote sensing image sequence with the medium and high spatial resolution acquired in each year comprises:
for the pixel area with multiple seeding index of 1 or 2, if any pixel area is the t th pixel areajPlanting information of year is planting from the tj+1If the planting information of the beginning of the year is no planting, the pixel area is determined as withdrawn farmland, and the time of withdrawing the farmland is tthjYear;
for the pixel area with multiple seeding index of 0-1, if any pixel area is the t th pixel areaiPlanting is carried out according to annual planting information, planting information is not planted within d years later, the pixel area is determined as withdrawn farmland, and the time of withdrawing the farmland is tthjAnd (5) year.
It can be understood that, according to the above further determined multiple cropping index of each pixel region, combining with the remote sensing image sequence data of the agriculture and forestry regions with medium and high spatial resolution, for example, extracting remote sensing images with different multiple cropping indexes by using medium and high spatial resolution, and determining the tillage withdrawal time and the tillage withdrawal area according to the remote sensing images, wherein, in any pixel region with multiple cropping index of 1 or 2, if the tth pixel region of any pixel region is found, the tillage withdrawal time is determined byjAnnual planting, tj+1Starting to stop planting, determining any pixel area as a withdrawn farmland, wherein the withdrawal time is tthjAnd (5) year. For any image element area with multiple cropping index of 0-1, if any, thenNow the t-th image element region of any one of the image elementsjPlanting in the year, stopping planting in the following d years, and judging that the back-plowing time is the tthjAnd (5) year.
Calculating the farmland returning area according to the pixel area judged as farmland returning, and calculating the farmland returning area by the following formula:
S1＝e.p2；
S1the area of the abandoned land is shown as e, the number of pixel areas of the abandoned land is shown as p, and the pixel spatial resolution of the remote sensing image is shown as p.
The method for returning the agricultural and forestry areas to the agricultural and forestry areas provided by the embodiment of the invention is combined with a remote sensing identification technology, adopts a time sequence agricultural and forestry area remote sensing image, utilizes a corresponding NDVI value to automatically calculate the positions of the returned cultivated lands, the returning cultivation time and the area of the returned cultivated lands, and utilizes an automatic remote sensing identification technology to improve the accuracy of information extraction of the returned cultivated lands.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (8)
1. A method for acquiring land withdrawal information of an agricultural and forestry area is characterized by comprising the following steps:
taking a year as a unit, extracting an NDVI value of each pixel region in a remote sensing image sequence according to the remote sensing image sequence with low space and high time resolution acquired in one year to obtain an NDVI sequence value corresponding to the remote sensing image sequence;
calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year;
determining planting information of any one of the image element areas according to the CV value of any one year and the NDVI value of any one of the image element areas of any one year;
determining a multiple cropping index of any pixel area in the monitored years according to the planting information of any pixel area in each year and the crop phenological information of any year in the monitored years;
and determining the back-plowing time, the back-plowing ground area and the back-plowing area position based on the multiple cropping index of each pixel region in the monitoring years and by combining the remote sensing image sequence with the medium and high spatial resolution acquired in each year.
2. The method for acquiring land evacuation information in an agricultural and forestry area according to claim 1, wherein the extracting an NDVI value of each pixel area in a remote sensing image sequence according to a low-space high-time resolution remote sensing image sequence acquired within one year by taking the year as a unit to obtain the NDVI sequence value corresponding to the remote sensing image sequence comprises:
extracting the NDVI value of each pixel region according to the remote sensing image sequence with low space and high time resolution in one year to obtain the NDVI sequence value corresponding to the remote sensing image sequence;
carrying out filtering processing on the NDVI sequence value to obtain a filtered NDVI sequence value;
correspondingly, calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value of any year comprises the following steps:
and calculating the CV value of any year and extracting the crop phenological information of any year according to the NDVI sequence value after filtering processing of any year.
3. The method for acquiring abandoned land information in an agricultural and forestry area according to claim 2, wherein the calculating the CV value of any year according to the NDVI sequence value after filtering processing in any year comprises:
wherein NDVIiIs the NDVI value, CV, of the i-th pixel area of the low-resolution remote sensing image sequence of the j-th yearjIs the coefficient of variation, σ, of the NDVI value of the pixel region of the j-th low-resolution remote-sensing image sequencejIs the variance, mu, of the NDVI value of the pixel region of the j-th remote sensing image sequencejThe average value of the NDVI values of the pixel areas of the j-th remote sensing image sequence is shown.
4. The method for acquiring the farmland returning information of the agricultural and forestry regions according to claim 2 or 3, wherein the determining the planting information of any one of the pixel areas according to the CV value of any one year and the NDVI value of any one of the pixel areas of any one year comprises:
if the ith pixel area in the jth year meets the NDVIiNot less than b and CVjAnd a is more than or equal to a, the crop planting information of the ith pixel area is planting, wherein a and b are preset threshold values.
5. The method for acquiring the land evacuation information of the agricultural and forestry area according to claim 4, wherein the multiple cropping index of any one pixel area in the monitored years is determined according to the planting information of any one pixel area in each year in the monitored years and the crop phenological information of any one year:
preliminarily determining a multiple cropping index of any pixel area in the monitored years based on the planting information of any pixel area in the remote sensing image sequence in each year;
and accurately determining the multiple cropping index of any pixel area in the monitoring years based on the crop phenological information of any year and the preliminarily determined multiple cropping index of any pixel area in the monitoring years.
6. The method for acquiring land evacuation information in an agricultural and forestry area according to claim 5, wherein the preliminary determination of the multiple cropping index of any pixel area within the monitored years based on the planting information of any pixel area in the remote sensing image sequence in each year comprises:
if the planting information of any pixel area in each year in the monitored years is planting, the multiple planting index of any pixel area is more than 1;
if the planting information of one part of the year is planting and the planting information of the other part of the year is not planting in any pixel area within the monitoring years, the multiple planting index of any pixel area is between 0 and 1.
7. The method for acquiring farmland returning information in an agricultural and forestry area according to claim 6, wherein the accurately determining the cropping index of any one of the pixel areas in the monitored years based on the crop phenological information of any one year and the preliminarily determined cropping index of any one of the pixel areas in the monitored years comprises:
based on the remote sensing image sequence with low space and high time resolution, extracting an NDVI sequence value corresponding to the remote sensing image sequence of each year;
generating a crop phenological curve of each year according to the NDVI sequence value;
determining a plurality of NDVI values of a plurality of growth starting points of the crops and a plurality of NDVI values of a high-value area in a growth vigorous period according to the crop climatic curve;
and determining the multiple cropping index of each remote sensing pixel area within the monitoring years to be 2 or 1 or between 0 and 1 according to the NDVI values of the multiple growth starting points of the crops and the NDVI values of the high-value areas in the vigorous growth period.
8. The method for acquiring land evacuation information of an agricultural and forestry area according to claim 7, wherein the determining of the time, the area and the position of the land evacuation includes, based on the multiple cropping index of each pixel area within the number of monitored years, in combination with a remote sensing image sequence of medium and high spatial resolution acquired each year:
for the pixel area with multiple seeding index of 1 or 2, if any pixel area is the t th pixel areajPlanting information of year is planting from the tj+1If the planting information of the beginning of the year is no planting, the pixel area is determined as withdrawn farmland, and the time of withdrawing the farmland is tthjYear;
for the pixel area with multiple seeding index of 0-1, if any pixel area is the t th pixel areajPlanting information of year is planting, and within d years later, if planting information is not planting, the time of back plowing is tthjYear;
calculating the land retreating area by the following formula:
S1＝e.p2；
S1for back-tillingThe area of the land, e is the number of pixel areas of the cultivated land, and p is the pixel spatial resolution of the remote sensing image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010761233.9A CN112001809A (en) | 2020-07-31 | 2020-07-31 | Method for acquiring farmland returning information of agriculture and forestry area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010761233.9A CN112001809A (en) | 2020-07-31 | 2020-07-31 | Method for acquiring farmland returning information of agriculture and forestry area |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112001809A true CN112001809A (en) | 2020-11-27 |
Family
ID=73463614
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010761233.9A Pending CN112001809A (en) | 2020-07-31 | 2020-07-31 | Method for acquiring farmland returning information of agriculture and forestry area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112001809A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113538388A (en) * | 2021-07-23 | 2021-10-22 | 中国电子科技集团公司第五十四研究所 | Arable land loss assessment method based on MODIS NDVI time sequence data |
| CN113609899A (en) * | 2021-06-24 | 2021-11-05 | 河海大学 | Backing-out land information positioning display system based on remote sensing time sequence analysis |
| CN114782816A (en) * | 2022-04-28 | 2022-07-22 | 天津大学 | Remote sensing extraction method for crop multiple cropping index |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102592181A (en) * | 2011-12-30 | 2012-07-18 | 中国农业科学院农业资源与农业区划研究所 | Method for optimizing spatial distribution of statistical data about crop planting area |
| CN104850694A (en) * | 2015-05-13 | 2015-08-19 | 福州大学 | Winter wheat remote sensing monitoring method based on vegetation index increment in growing period |
| US20180189564A1 (en) * | 2016-12-30 | 2018-07-05 | International Business Machines Corporation | Method and system for crop type identification using satellite observation and weather data |
| CN108345992A (en) * | 2018-01-31 | 2018-07-31 | 北京师范大学 | A kind of multiple crop index extracting method and device |
| CN109829234A (en) * | 2019-01-30 | 2019-05-31 | 北京师范大学 | A kind of across scale Dynamic High-accuracy crop condition monitoring and yield estimation method based on high-definition remote sensing data and crop modeling |
| 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 |
-
2020
- 2020-07-31 CN CN202010761233.9A patent/CN112001809A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102592181A (en) * | 2011-12-30 | 2012-07-18 | 中国农业科学院农业资源与农业区划研究所 | Method for optimizing spatial distribution of statistical data about crop planting area |
| CN104850694A (en) * | 2015-05-13 | 2015-08-19 | 福州大学 | Winter wheat remote sensing monitoring method based on vegetation index increment in growing period |
| US20180189564A1 (en) * | 2016-12-30 | 2018-07-05 | International Business Machines Corporation | Method and system for crop type identification using satellite observation and weather data |
| CN108345992A (en) * | 2018-01-31 | 2018-07-31 | 北京师范大学 | A kind of multiple crop index extracting method and device |
| CN109829234A (en) * | 2019-01-30 | 2019-05-31 | 北京师范大学 | A kind of across scale Dynamic High-accuracy crop condition monitoring and yield estimation method based on high-definition remote sensing data and crop modeling |
| 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 |
Non-Patent Citations (2)
| Title |
|---|
| 王璇: "基于MODIS时间序列的河南主要农作物种植信息提取", 《中国优秀硕士学位论文全文数据库农业科技辑》 * |
| 许玉凤 等: "贵州省黔南州植被覆盖时空变化分析", 《地理空间信息》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113609899A (en) * | 2021-06-24 | 2021-11-05 | 河海大学 | Backing-out land information positioning display system based on remote sensing time sequence analysis |
| CN113609899B (en) * | 2021-06-24 | 2024-05-14 | 河海大学 | Remote sensing time sequence analysis-based tilling land information positioning display system |
| CN113538388A (en) * | 2021-07-23 | 2021-10-22 | 中国电子科技集团公司第五十四研究所 | Arable land loss assessment method based on MODIS NDVI time sequence data |
| CN114782816A (en) * | 2022-04-28 | 2022-07-22 | 天津大学 | Remote sensing extraction method for crop multiple cropping index |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tarolli et al. | High-resolution morphologic characterization of conservation agriculture | |
| Tellaeche et al. | A new vision-based approach to differential spraying in precision agriculture | |
| CN112001809A (en) | Method for acquiring farmland returning information of agriculture and forestry area | |
| CN110909679B (en) | Remote sensing identification method and system for fallow crop rotation information of winter wheat historical planting area | |
| CN109492619B (en) | Variable pesticide application method and system integrating remote sensing, model and algorithm | |
| Vélez et al. | Mapping the spatial variability of Botrytis bunch rot risk in vineyards using UAV multispectral imagery | |
| CN113222991A (en) | Deep learning network-based field ear counting and wheat yield prediction | |
| CN113221806B (en) | Cloud platform fusion multi-source satellite image and tea tree phenological period based automatic tea garden identification method | |
| Ma et al. | Automatic detection of crop root rows in paddy fields based on straight-line clustering algorithm and supervised learning method | |
| CN110728446A (en) | County scale crop yield estimation method based on CNN-LSTM | |
| CN108345992A (en) | A kind of multiple crop index extracting method and device | |
| CN114818909A (en) | Weed detection method and device based on crop growth characteristics | |
| CN112232230B (en) | Remote sensing identification method and device for land withdrawal information | |
| CN112418050B (en) | Remote sensing identification method and device for land withdrawal information | |
| JP2022082636A (en) | Information processing device | |
| Narra et al. | A data driven approach to decision support in farming | |
| Salma et al. | Identifying rice crop flooding patterns using Sentinel-1 SAR data | |
| Zhang et al. | Automatic counting of lettuce using an improved YOLOv5s with multiple lightweight strategies | |
| CN116385885A (en) | Method and system for identifying plant diseases and insect pests | |
| CN115527111A (en) | Remote sensing extraction method, device and equipment for rice planting information and multiple cropping index | |
| Mangla et al. | Statistical growth prediction analysis of rice crop with pixel-based mapping technique | |
| Bulanon et al. | Predicting fruit yield using shallow neural networks | |
| CN113984212A (en) | Agricultural irrigated area extraction method and system | |
| Lou et al. | Combining shape and crop models to detect soybean growth stages | |
| Poleshchenko et al. | Development of a System for Automated Control of Planting Density, Leaf Area Index and Crop Development Phases by UAV Photos |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201127 |