CN112148785B - Crop distribution charting method - Google Patents
Crop distribution charting method Download PDFInfo
- Publication number
- CN112148785B CN112148785B CN202011016570.1A CN202011016570A CN112148785B CN 112148785 B CN112148785 B CN 112148785B CN 202011016570 A CN202011016570 A CN 202011016570A CN 112148785 B CN112148785 B CN 112148785B
- Authority
- CN
- China
- Prior art keywords
- crop
- mode
- area
- crops
- planting
- 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.)
- Expired - Fee Related
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims description 12
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 238000009336 multiple cropping Methods 0.000 claims description 22
- 238000009355 double cropping Methods 0.000 claims description 16
- 238000009356 triple cropping Methods 0.000 claims description 8
- 238000010899 nucleation Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000007670 refining Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000003306 harvesting Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000001932 seasonal effect Effects 0.000 claims description 2
- 241000405911 Rehmannia glutinosa Species 0.000 claims 1
- 241000209140 Triticum Species 0.000 description 11
- 235000021307 Triticum Nutrition 0.000 description 11
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000035800 maturation Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005070 ripening Effects 0.000 description 2
- 238000009331 sowing Methods 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/26—Visual data mining; Browsing structured data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a crop distribution map, which is characterized in that the area of various crops on each grid is calculated according to grid scale crop planting mode information extracted from remote sensing data and the area proportion of various crops on the planting mode of a corresponding administrative unit, so that a crop space distribution map is manufactured. The invention combines the crop phenology characteristics identified by remote sensing with the actual phenology of the crops to determine the spatial distribution condition of the crops, and can realize large-range and multi-scale crop spatial distribution mapping with higher spatial heterogeneity through the characteristics of macro, quick acquisition, multi-scale and the like of remote sensing.
Description
Technical Field
The invention relates to the field of data analysis, in particular to a crop distribution chart.
Background
The crop space-time distribution reflects the situation that human beings use land to engage in agricultural production in space at different periods, and is important basic data for developing researches on the patterns and functions of farmland ecosystems, the circulation of the land ecosystems, global changes and the like. The understanding of the crop space-time distribution has important significance for guaranteeing national food safety and sustainable development of resource environment.
At present, the crop space-time distribution information acquisition approach mainly comprises administrative statistics and remote sensing monitoring. The statistical method can only reflect the quantity change of the crops in a certain administrative unit, and the spatial difference of the crops in the statistical unit is difficult to reflect; moreover, the statistical data acquisition consumes a large amount of manpower, material resources and financial resources, and is also interfered by human factors. The simple remote sensing technology has the influence of factors such as atmospheric interference, scale conversion, mixed pixels and the like, and is difficult to acquire the spatial distribution information of the large-area-scale full crop series. Therefore, the statistical data and the remote sensing data are deeply fused, and the development of the spatial distribution mapping of multiple crops and even the whole crop type is particularly necessary.
Disclosure of Invention
The invention aims to provide a crop distribution mapping method,
in order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the invention comprises the following steps
A, acquiring remote sensing data information, constructing a vegetation index time sequence curve of a pixel scale, extracting multiple species information through the number of wave crests of the time sequence curve of the pixel scale, and extracting, classifying and refining planting season information reflected by the pixel scale according to the occurrence time of the wave crests to obtain planting mode information of each raster crop;
b, acquiring the seeding area data of the administrative unit crops, merging the crops in different seasons of the administrative unit, and obtaining the seeding area information of the crops in different seasons of various planting modes;
and C, calculating the areas of various crops on each grid according to the grid scale crop planting mode extracted by remote sensing, the corresponding cultivated area distribution of the grid scale crop planting mode and the area information of the seasonal crops on different planting modes, and manufacturing a crop space distribution map.
Further, the characteristic parameters comprise the number of time sequence curve peaks and time points of the time sequence curve peaks, and the grid planting mode information is obtained by presumption according to the characteristic parameters:
specifically, the preprocessing includes denoising the data by using a time series harmonic analysis method to smooth a time series curve.
Further, the specific flow of pixel scale extraction is as follows:
let each time phase pixel value of pixel i be a1,a2…a23,
Comparing EVI magnitude relation between two adjacent pixels, and using delta aiThe sequence data represents:
② calculating adjacent delta aiThe difference between:
ΔΔai=Δai+1-Δai (3-8)
wherein, if Δ Δ ai<0, then i-1 is the time point of the local peak, ai-1For local peak EVIA value; if Δ Δ ai>0, then i-1 is the time point of the local trough, ai-1Is the local trough EVI value.
And thirdly, because of the existence of the winter peak of the crop and the peak information of some interference, the real peak point of the growing period of the crop needs to be further extracted. If the peak point satisfies that the EVI value of two continuous time points increases before and the EVI value of two continuous time points decreases after, the point belongs to the EVI peak point of the growing period of the crops. The sum of the effective peak points is the multimaturity information of the pixel, namely the peak number is 1, and the pixel is mature for one year; the number of wave peaks is 2, and the product is twice cooked in one year; the number of wave peaks is more than 2, the temperature is more mature in one year, the temperature of the secondary maturation zone is required to be more than 10 ℃ and more than 3500 ℃, and the temperature of the secondary maturation zone is required to be more than 10 ℃ and more than 5000 ℃.
Specifically, the administrative unit divides crops into winter crops, spring crops and summer crops according to the lunar calendar information of the crop types, and the crops are corresponding to various planting mode refinements existing in the administrative unit, and the planting mode refinement method comprises the following steps
a, for the area with triple cropping potential, if the extracted multiple cropping result is triple cropping, dividing the plot into three-season crop multiple cropping in winter-spring-summer; if the extracted multiple cropping result is double cropping, extracting the time point of occurrence of the peak value of the extracted time series curve, determining whether the peak value of crop growth reflected by the curve occurs in spring, summer or autumn, and further dividing the plot into a winter-summer multiple cropping mode or a spring-summer multiple cropping mode; and if the extracted multiple cropping result is one-time cropping, extracting a peak occurrence point of the time series curve, and determining that the planting mode is a winter mode, a spring mode or a summer mode.
b, for the area with double cropping potential, if the double cropping result extracted in the step 2 is double cropping, dividing the plot into a winter-summer double cropping mode; and if the extracted multiple cropping result is one-time cropping, extracting a peak occurrence point of the time series curve, and determining that the planting mode is a winter mode, a spring mode or a summer mode.
And c, directly extracting peak occurrence points for the region which can only realize one-time harvest, and determining whether the peak occurrence points are in a spring mode or a summer mode.
Further, the crop space distribution map is based on the following formula, the grid scale crop area is calculated, and spatialization is carried out to form the crop space distribution map.
Ail=CAjl×HIiik
Wherein A isijRepresents the area of crop i in grid l; CAikRepresenting the area of the planting pattern j in the grid l; HI (high-intensity)iikRepresents the area proportion of the plot of the crop i in the planting pattern j in the grid l, wherein the grid 1 is extracted by the administrative unit.
In particular, the area ratio of various crops in a unit area of cultivation on a plot with different planting patterns (HI)ijk) The formula is as follows:
HIijk=HAikj/CAjk
wherein, HAilIndicating the area of sowing of the crop i in the administration unit k, CAklThe area of the manifold pattern j in the administration unit k is shown.
Further, in the step A, curve fitting is carried out according to different multiple cropping potential areas by adopting corresponding frequency numbers, a second-maturing area adopts 3, and a third-maturing area adopts 4.
Compared with the prior art, the invention has the following beneficial effects:
the invention combines the crop phenology characteristics identified by remote sensing with the actual phenology of the crops to determine the spatial distribution condition of the crops, realizes the crop spatial distribution mapping with large range, multiple scales and higher spatial heterogeneity by the characteristics of macro, quick acquisition, multiple scales and the like of remote sensing, and can acquire the spatial distribution information of the crops which have both spatial heterogeneity information and full crop information by the combination of deep mining of remote sensing data and statistical data.
Drawings
FIG. 1 is a schematic diagram of a crop distribution charting process provided by the present invention;
Detailed Description
The present invention will be further illustrated by the following examples and the accompanying drawings, which include, but are not limited to, the following examples.
Referring to FIG. 1, in the present embodiment, the winter wheat is used as an example to perform a crop distribution chart
Step 1: data are denoised by a time series harmonic analysis method (harmonic analysis of time series, HANTS) to obtain a smooth time series curve. And (3) performing curve fitting in different multi-cropping potential areas by adopting different frequencies (Numberoff sequences), wherein the second-maturing area adopts 3, and the third-maturing area adopts 4.
Step 2: after curve denoising, extracting multiple species information in pixel scale, wherein the specific flow is as follows:
let each time phase pixel value of pixel i be a1,a2…a23,
Comparing EVI magnitude relation between two adjacent pixels, and using delta aiThe sequence data represents:
② calculating adjacent delta aiThe difference between:
ΔΔai=Δai+1-Δai (3-8)
wherein, if Δ Δ ai<0, then i-1 is the time point of the local peak, ai-1Is the local peak EVI value; if Δ Δ ai>0, then i-1 is the time point of the local trough, ai-1Is the local trough EVI value.
And thirdly, because the information of the winter peak of winter wheat and some interfering wave peaks exists, the real peak point of the crop in the growing period needs to be further extracted. If the peak point satisfies that the EVI value of two continuous time points increases before and the EVI value of two continuous time points decreases after, the point belongs to the EVI peak point of the growing period of the crops. The sum of the effective peak points is the multimaturity information of the pixel, namely the peak number is 1, and the pixel is mature for one year; the number of wave peaks is 2, and the product is twice cooked in one year; if the number of peaks is greater than 2, the product is more mature than one year.
The temperature of the secondary ripening area is required to be higher than 3500 ℃ when the temperature is more than or equal to 10 ℃, and the temperature of the secondary ripening area is required to be higher than 5000 ℃ when the temperature is more than or equal to 10 ℃.
And step 3: on the basis of extracting the multiple cropping information, the planting season information reflected by the pixel is extracted, and then the planting mode is extracted.
1 for the area with triple cropping potential, if the triple cropping result extracted in the step 2 is triple cropping, dividing the plot into winter-spring-summer triple crop multiple cropping; if the extracted multiple cropping result is double cropping, extracting the time point of occurrence of the peak value of the time series curve extracted in the step 2, determining whether the peak value of the crop growth reflected by the curve occurs in spring, summer or autumn, and further dividing the plot into a winter-summer multiple cropping mode or a spring-summer multiple cropping mode; and if the extracted multiple cropping result is one-time cropping, extracting a peak occurrence point of the time series curve, and determining that the planting mode is a winter mode, a spring mode or a summer mode.
2 for the area with double cropping potential, if the double cropping result extracted in the step 2 is double cropping, dividing the plot into a winter-summer double cropping mode; and if the extracted multiple cropping result is one-time cropping, extracting a peak occurrence point of the time series curve, and determining that the planting mode is a winter mode, a spring mode or a summer mode.
And 3, directly extracting peak occurrence points for the region which can only realize the first-maturing, and determining the mode to be a spring mode or a summer mode.
And 4, step 4: and calculating the cultivated land area with the winter crop planting mode in each administrative unit by taking the administrative unit with the obtained statistical data of the winter wheat seeding area as a boundary.
And 5: according to the farming season information of the crop types contained in each administrative unit, the crops are divided into winter crops, spring crops and summer crops, and the winter crops correspond to the planting modes of the winter crops existing in the administrative unit.
Step 6: and calculating the area proportion of the winter wheat in each administrative statistical unit in the land plots of various planting modes. According to the planting mode corresponding to the winter crops in each administrative unit determined in the step 5, combining the area of the planting mode land parcel in the administrative unit obtained in the step 4, calculating the unit cultivated land area on the land parcels with different planting modesWheat area ratio (HI)ijk) The formula is as follows:
HIjk=HAkj/CAjk (1)
wherein, HAilIndicates the sowing area of winter wheat in administrative Unit k, CAklRepresents the area of the planting pattern j containing winter crops in the administration k.
And 7: and calculating the areas of the winter wheat on the grids in different planting modes according to the grid planting mode information and the corresponding farmland distribution and combining the areas of the winter wheat crops on the planting modes of the administrative units, thereby manufacturing a winter wheat space distribution map. Specifically, the area of the winter wheat is calculated in a grid scale based on the following formula, and spatialization is performed to form a crop space distribution map.
Al=CAjl×HIjl
Wherein A islRepresents the area of winter wheat in grid l; CAljRepresenting the area of the planting pattern j in the grid l; HI (high-intensity)ljAnd (3) representing the area proportion of the winter wheat in the plot of the planting mode j in the grid l, wherein the grid l belongs to the administrative unit k and is obtained in the step 6.
It will be clear to a person skilled in the art that the scope of the present invention is not limited to the examples discussed in the foregoing, but that several amendments and modifications thereof are possible without deviating from the scope of the present invention as defined in the attached claims. While the invention has been illustrated and described in detail in the drawings and the description, such illustration and description are to be considered illustrative or exemplary and not restrictive. The invention is not limited to the disclosed embodiments.
Variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the term "comprising" does not exclude other steps or elements. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage. Any reference signs in the claims shall not be construed as limiting the scope of the invention.
The above-mentioned embodiment is only one of the preferred embodiments of the present invention, and should not be used to limit the scope of the present invention, but all the insubstantial modifications or changes made within the spirit and scope of the main design of the present invention, which still solve the technical problems consistent with the present invention, should be included in the scope of the present invention.
Claims (6)
1. The crop distribution charting method is characterized by comprising the following steps
A, acquiring remote sensing data information, constructing a vegetation index time series curve of a pixel scale, extracting, classifying and refining planting season information reflected by the pixel scale according to peak occurrence time through extracting multiple planting information of the number of peaks of the pixel scale time series curve to obtain planting mode information of each raster crop, wherein the specific flow of extracting the multiple planting information by the pixel scale is as follows:
let each time phase pixel value of pixel i be a1,a2… a23,
Comparing the EVI magnitude relation between two adjacent pixels and usingThe sequence data represents:
wherein, if<0, theni+1 is the point in time where the local peak is located,a i+1is the local peak EVI value; if it is>0, then i+1 is the point in time at which the local trough is located, a i+1is the local trough EVI value;
thirdly, because of the existence of the winter peak of the crop and some interfered peak information, further extracting the real peak point of the growing period of the crop, if the peak point meets the condition that the EVI value of two continuous time points is increased before and the EVI value of two continuous time points is reduced later, the peak point belongs to the EVI peak point of the growing period of the crop, and the sum of the effective peak points is the multimaturity information of the pixel, namely the number of the peak is 1, and the maturity is one year; the number of wave peaks is 2, and the product is twice cooked in one year; if the number of wave peaks is more than 2, the product is more mature in one year;
b, acquiring seeding area data of crops in an administrative unit, merging the crops in different seasons in the administrative unit, obtaining seeding area information of the crops in different seasons in various planting modes, dividing the crops into winter crops, spring crops and summer crops by the administrative unit according to the information of the agricultural time of the types of the crops, and refining the crops to correspond to various planting modes existing in the administrative unit, wherein the method for refining the planting modes comprises the following steps:
a, for the area with triple cropping potential, if the extracted multiple cropping result is a triple cropping plot, dividing the triple cropping plot into winter-spring-summer three-season crop multiple cropping; if the extracted multiple cropping result is a double cropping plot, extracting the time point of occurrence of the time series curve peak value of the double cropping plot, determining whether the crop growth peak value reflected by the curve occurs in spring, summer or autumn, and further dividing the two rehmannia glutinosa plots into a winter-summer multiple cropping mode or a spring-summer multiple cropping mode; if the extracted multiple cropping result is a first-maturing plot, extracting a peak occurrence point of a time series curve of the first-maturing plot, and determining that the planting mode is a winter mode, a spring mode or a summer mode;
b, for the area with double cropping potential, if the extracted double cropping result is a double cropping plot, dividing the plot into a winter-summer double cropping mode; if the extracted multiple cropping result is a first-maturing plot, extracting a peak occurrence point of the time series curve, and determining that the planting mode is a winter mode, a spring mode or a summer mode;
c, directly extracting peak occurrence points for the region which can only realize one-time harvest, and determining the mode is a spring mode or a summer mode;
and C, calculating the areas of various crops on each grid according to the grid scale crop planting mode extracted by remote sensing, the corresponding cultivated area distribution of the grid scale crop planting mode and the area information of the seasonal crops on different planting modes, and manufacturing a crop space distribution map.
2. The method according to claim 1, wherein the characteristic parameters include the number of peaks of the time series curve and the time points of the peaks of the time series curve, and the information of the grid planting pattern is obtained by estimation according to the characteristic parameters.
3. The method of claim 1, wherein the preprocessing comprises denoising the data by time series harmonic analysis to smooth the time series curve.
4. The method for mapping the crop distribution according to claim 1, wherein the crop space distribution map is formed by calculating the area of the crop with the grid dimension based on the following formula and spatializing the area of the crop with the grid dimension:
A il = CA jl ×HI ijk
wherein A is il Presentation gridlChinese cropiThe area of (d); CA jl Presentation gridlMiddle planting modejThe area of (d); HI (high-intensity) ijk Presentation gridlChinese cropiIn the planting modejArea ratio of land mass, gridl Extracted by administrative unit k.
5. The method for crop distribution mapping according to claim 1, wherein the area ratio of each crop in unit area of cultivated land on the land of different planting modes is as follows:
HI ijk = HA ik /CA jk
wherein, HA ik Representing cropsiIn administrative unitkSeeding area of, CA jk Indicating planting patternsjIn administrative unitkThe area of (a).
6. The method for mapping crop distribution according to claim 1, wherein in step a, curve fitting is performed according to different multi-cropping potential areas by using corresponding frequency numbers, the second cropping area is 3, and the third cropping area is 4.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011016570.1A CN112148785B (en) | 2020-09-24 | 2020-09-24 | Crop distribution charting method |
AU2020103047A AU2020103047A4 (en) | 2020-09-24 | 2020-10-27 | Crop Distribution Mapping |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011016570.1A CN112148785B (en) | 2020-09-24 | 2020-09-24 | Crop distribution charting method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112148785A CN112148785A (en) | 2020-12-29 |
CN112148785B true CN112148785B (en) | 2021-09-28 |
Family
ID=73838804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011016570.1A Expired - Fee Related CN112148785B (en) | 2020-09-24 | 2020-09-24 | Crop distribution charting method |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112148785B (en) |
AU (1) | AU2020103047A4 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114092804A (en) * | 2021-11-06 | 2022-02-25 | 支付宝(杭州)信息技术有限公司 | Method and device for identifying remote sensing image |
CN114067158B (en) * | 2021-11-17 | 2022-09-23 | 江苏天汇空间信息研究院有限公司 | Farmland use state monitoring system and method applying multi-source remote sensing data |
CN114739943A (en) * | 2022-05-10 | 2022-07-12 | 福州大学 | Accurate extraction method of cooking information based on full life cycle of crop development and maturity |
CN114926565B (en) * | 2022-05-24 | 2024-04-12 | 中国科学院东北地理与农业生态研究所 | Soybean drafting method based on green degree and water comprehensive index |
CN114782838B (en) * | 2022-06-17 | 2022-10-18 | 中化现代农业有限公司 | Rice identification method and device, electronic equipment and storage medium |
Citations (7)
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 |
CN102855351A (en) * | 2012-08-09 | 2013-01-02 | 中国科学院地理科学与资源研究所 | Crop straw resource spatialization method based on statistical data and remotely-sensed data |
CN108345992A (en) * | 2018-01-31 | 2018-07-31 | 北京师范大学 | A kind of multiple crop index extracting method and device |
CN109115770A (en) * | 2018-06-14 | 2019-01-01 | 中科禾信遥感科技(苏州)有限公司 | A kind of a wide range of crops remote-sensing monitoring method and device |
CN109360117A (en) * | 2018-10-08 | 2019-02-19 | 西充恒河农牧业开发有限公司 | A kind of crop growing mode recognition methods |
US10303677B2 (en) * | 2015-10-14 | 2019-05-28 | The Climate Corporation | Computer-generated accurate yield map data using expert filters and spatial outlier detection |
CN110175931A (en) * | 2019-05-10 | 2019-08-27 | 北京师范大学 | A kind of method of a wide range of rapidly extracting crop acreage and phenology information |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103500421B (en) * | 2013-10-09 | 2017-01-11 | 福州大学 | Frequency characteristic-based farmland cropping index extraction method |
CN108764255B (en) * | 2018-05-21 | 2021-06-01 | 二十一世纪空间技术应用股份有限公司 | Method for extracting winter wheat planting information |
CN110443420B (en) * | 2019-08-05 | 2023-05-09 | 山东农业大学 | Crop yield prediction method based on machine learning |
CN111598019B (en) * | 2020-05-19 | 2023-05-26 | 华中农业大学 | Crop type and planting mode identification method based on multi-source remote sensing data |
-
2020
- 2020-09-24 CN CN202011016570.1A patent/CN112148785B/en not_active Expired - Fee Related
- 2020-10-27 AU AU2020103047A patent/AU2020103047A4/en not_active Ceased
Patent Citations (7)
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 |
CN102855351A (en) * | 2012-08-09 | 2013-01-02 | 中国科学院地理科学与资源研究所 | Crop straw resource spatialization method based on statistical data and remotely-sensed data |
US10303677B2 (en) * | 2015-10-14 | 2019-05-28 | The Climate Corporation | Computer-generated accurate yield map data using expert filters and spatial outlier detection |
CN108345992A (en) * | 2018-01-31 | 2018-07-31 | 北京师范大学 | A kind of multiple crop index extracting method and device |
CN109115770A (en) * | 2018-06-14 | 2019-01-01 | 中科禾信遥感科技(苏州)有限公司 | A kind of a wide range of crops remote-sensing monitoring method and device |
CN109360117A (en) * | 2018-10-08 | 2019-02-19 | 西充恒河农牧业开发有限公司 | A kind of crop growing mode recognition methods |
CN110175931A (en) * | 2019-05-10 | 2019-08-27 | 北京师范大学 | A kind of method of a wide range of rapidly extracting crop acreage and phenology information |
Non-Patent Citations (2)
Title |
---|
农作物空间分布遥感制图发展方向探讨;宋茜等;《中国农业资源与区划》;20200625;第41卷(第6期);第57-65页 * |
基于MODIS/EVI 的中国北方耕地复种指数提取;左丽君等;《农业工程学报》;20090831(第8期);第141-146页 * |
Also Published As
Publication number | Publication date |
---|---|
AU2020103047A4 (en) | 2020-12-24 |
CN112148785A (en) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112148785B (en) | Crop distribution charting method | |
CN105760978B (en) | One kind being based on the agricultural drought disaster grade monitoring method of temperature vegetation drought index (TVDI) | |
CN111368736B (en) | Rice refined estimation method based on SAR and optical remote sensing data | |
Tao et al. | Mapping winter wheat using phenological feature of peak before winter on the North China Plain based on time-series MODIS data | |
CN106529451B (en) | Winter wheat-summer corn planting mode remote sensing identification method | |
CN114387516B (en) | Single-season rice SAR (synthetic aperture radar) identification method for small and medium-sized fields in complex terrain environment | |
CN105718751B (en) | Deep scarification operation area calculating method based on area grid point coverage | |
CN111091079B (en) | TLS-based method for measuring vegetation advantage single plant structural parameters in friable region | |
CN109063553A (en) | Field-crop growth defect area's remote sensing fast diagnosis method after a kind of land control | |
Zhao et al. | Using satellite remote sensing to understand maize yield gaps in the North China Plain | |
Kotsuki et al. | SACRA–a method for the estimation of global high-resolution crop calendars from a satellite-sensed NDVI | |
Wang et al. | Spatial-temporal evolution of vegetation evapotranspiration in Hebei Province, China | |
Dan et al. | Delineating the rice crop activities in Northeast China through regional parametric synthesis using satellite remote sensing time-series data from 2000 to 2015 | |
Wang et al. | Rice yield estimation based on an NPP model with a changing harvest index | |
CN102156886A (en) | Statistical data and remote sensing image data based regional fertilizer application spatialization method | |
Sarkar et al. | Water quality impacts of converting intensively-managed agricultural lands to switchgrass | |
Anderson et al. | Evaluating global land degradation using ground-based measurements and remote sensing | |
CN112801924A (en) | Method for constructing surface pure microwave radiation background field | |
Bolun et al. | Estimating rice paddy areas in China using multi-temporal cloud-free normalized difference vegetation index (NDVI) imagery based on change detection | |
Dalposso et al. | Spatial autocorrelation of NDVI and GVI indices derived from Landsat/TM images for soybean crops in the western of the state of Paraná in 2004/2005 crop season | |
CN110929222A (en) | Irrigation farmland identification method based on remote sensing vegetation canopy moisture index | |
CN117611993B (en) | Method for estimating vegetation classification based on remote sensing actual evapotranspiration | |
Chen et al. | Tracking the spatio-temporal change of the main food crop planting structure in the Yellow River Basin over 2001–2020 | |
CN110598514A (en) | Method for monitoring plot scale crop seeding area of land reclamation project area | |
CN113570273A (en) | Spatial method and system for irrigation farmland statistical data |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210928 |