CN113822360B - Method and device for determining phenological period of crops - Google Patents

Method and device for determining phenological period of crops Download PDF

Info

Publication number
CN113822360B
CN113822360B CN202111123393.1A CN202111123393A CN113822360B CN 113822360 B CN113822360 B CN 113822360B CN 202111123393 A CN202111123393 A CN 202111123393A CN 113822360 B CN113822360 B CN 113822360B
Authority
CN
China
Prior art keywords
curve
phenological
target crop
crop
period
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.)
Active
Application number
CN202111123393.1A
Other languages
Chinese (zh)
Other versions
CN113822360A (en
Inventor
李福根
王宏斌
丁传艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinochem Agriculture Holdings
Original Assignee
Sinochem Agriculture Holdings
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinochem Agriculture Holdings filed Critical Sinochem Agriculture Holdings
Priority to CN202111123393.1A priority Critical patent/CN113822360B/en
Publication of CN113822360A publication Critical patent/CN113822360A/en
Application granted granted Critical
Publication of CN113822360B publication Critical patent/CN113822360B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/22—Matching criteria, e.g. proximity measures
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Animal Husbandry (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Agronomy & Crop Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method and a device for determining the phenological period of crops. The method comprises the steps of obtaining a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve; according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve; determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve; and determining the phenological period of the target crop according to the time alignment relation information. The method determines the phenological period of the target crop by using the target crop phenological curve of the target crop and the reference crop phenological curve corresponding to the target crop phenological curve, is more accurate and has more reference value, avoids the problems of subjectivity, sidedness and limitation caused by manual evaluation according to experience, and improves the determining efficiency of the phenological period of the target crop.

Description

Method and device for determining phenological period of crops
Technical Field
The invention relates to the technical field of agriculture, in particular to a method and a device for determining the phenological period of crops.
Background
Accurate determination of crop phenological information is critical to field management and decision making. The phenological period condition of crops is an important link between the crop ecosystem and yield estimation, so that accurate determination of crop phenological information is of great importance to field management and decision making. The traditional methods for extracting the crop phenological information mainly adopt manual recording, meteorological station observation and the like, and the methods are often accompanied by the defects of high cost, long time consumption and the like, so that the efficiency of extracting the crop phenological information is low, and the traditional methods are not suitable for acquiring the phenological information of crops in a large area. Meanwhile, due to the influence of latitude and climate factors, it is difficult to directly extend the phenological period of a single-point crop into a large-area. Therefore, a new method for determining the phenological period of crops is needed.
Disclosure of Invention
In view of this, embodiments of the present invention provide a method, an apparatus, a computer device, and a computer readable storage medium for determining a phenological period of a crop, so as to solve the problem in the prior art that the efficiency of extracting phenological information of the crop is low.
In a first aspect of the embodiments of the present invention, there is provided a method for determining a phenological period of a crop, the method including:
acquiring a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve;
according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve;
determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve;
and determining the phenological period of the target crop according to the time alignment relation information.
Optionally, the obtaining a target crop phenological curve of the target crop and a reference crop phenological curve corresponding to the target crop phenological curve includes:
acquiring a plurality of target crop phenological curves of the target crops;
generating an average crop phenological curve of the target crop according to a plurality of target crop phenological curves of the target crop;
and taking the average crop phenological curve of the target crop as a reference crop phenological curve of the target crop.
Optionally, the performing waveform adjustment on the target crop phenological curve according to the reference crop phenological curve to obtain an adjusted target crop phenological curve includes:
determining a first minimum value and a second minimum value of the reference crop phenological curve according to the maximum value of the reference crop phenological curve;
determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve;
and adjusting the amplitude of the waveform of the target crop phenological curve according to the maximum value, the first minimum value and the second minimum value of the reference crop phenological curve and the maximum value, the first minimum value and the second minimum value of the target crop phenological curve to obtain an adjusted target crop phenological curve.
Optionally, the determining the time alignment relationship information between the adjusted target crop phenological curve and the reference crop phenological curve includes:
determining a plurality of distance elements according to the adjusted target crop phenological curve, the reference crop phenological curve and a preset time weight factor;
determining a plurality of distance matrixes according to the distance elements;
and determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve according to the distance matrixes.
Optionally, the determining the phenological period of the target crop according to the time alignment relationship information includes:
determining the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve;
determining the green turning period node information, heading period node information and maturity period node information of the target crop phenological curve according to the time alignment relation information, the green turning period node information, heading period node information and maturity period node information of the reference crop phenological curve;
and determining the phenological period of the target crop according to the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve.
Optionally, determining the green-turning period node information, heading period node information and maturity period node information of the target crop phenological curve according to the time alignment relationship information, the green-turning period node information, heading period node information and maturity period node information of the target crop phenological curve includes:
according to the time alignment relation information and the green return period node information of the target crop phenological curve, determining alignment information of the green return period node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining the node information of the green turning period of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the heading stage node information of the target crop phenological curve, determining alignment information of the heading stage node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining heading stage node information of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the maturity node information of the target crop phenological curve, determining the alignment information of the maturity node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; and determining the maturity node information of the objective crop phenological curve according to the alignment information.
Optionally, the determining the phenological period of the target crop according to the node information of the green-turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green-turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve includes:
for each target crop phenological curve, determining difference information between the target crop phenological curve and the reference crop phenological curve according to the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the target crop phenological curve, and the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the reference crop phenological curve;
and determining the phenological period of the target crop according to the difference information corresponding to each target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve.
In a second aspect of embodiments of the present invention, there is provided an apparatus for determining a phenological period of a crop, the apparatus comprising:
the curve acquisition unit is used for acquiring a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve;
the curve adjusting unit is used for carrying out waveform adjustment on the target crop phenological curve according to the reference crop phenological curve to obtain an adjusted target crop phenological curve;
the information determining unit is used for determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve;
and the phenological period determining unit is used for determining the phenological period of the target crop according to the time alignment relation information.
Optionally, the curve obtaining unit is configured to:
acquiring a plurality of target crop phenological curves of the target crops;
generating an average crop phenological curve of the target crop according to a plurality of target crop phenological curves of the target crop;
and taking the average crop phenological curve of the target crop as a reference crop phenological curve of the target crop.
Optionally, the curve adjusting unit is configured to:
determining a first minimum value and a second minimum value of the reference crop phenological curve according to the maximum value of the reference crop phenological curve;
determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve;
and adjusting the amplitude of the waveform of the target crop phenological curve according to the maximum value, the first minimum value and the second minimum value of the reference crop phenological curve and the maximum value, the first minimum value and the second minimum value of the target crop phenological curve to obtain an adjusted target crop phenological curve.
Optionally, the information determining unit is configured to:
determining a plurality of distance elements according to the adjusted target crop phenological curve, the reference crop phenological curve and a preset time weight factor;
determining a plurality of distance matrixes according to the distance elements;
and determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve according to the distance matrixes.
Optionally, the phenological period determining unit is configured to:
determining the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve;
determining the green turning period node information, heading period node information and maturity period node information of the target crop phenological curve according to the time alignment relation information, the green turning period node information, heading period node information and maturity period node information of the reference crop phenological curve;
and determining the phenological period of the target crop according to the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve.
Optionally, the phenological period determining unit is configured to:
according to the time alignment relation information and the green return period node information of the target crop phenological curve, determining alignment information of the green return period node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining the node information of the green turning period of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the heading stage node information of the target crop phenological curve, determining alignment information of the heading stage node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining heading stage node information of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the maturity node information of the target crop phenological curve, determining the alignment information of the maturity node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; and determining the maturity node information of the objective crop phenological curve according to the alignment information.
Optionally, if the target crop corresponds to a plurality of target crop phenological curves; the phenological period determination unit is configured to:
for each target crop phenological curve, determining difference information between the target crop phenological curve and the reference crop phenological curve according to the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the target crop phenological curve, and the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the reference crop phenological curve;
and determining the phenological period of the target crop according to the difference information corresponding to each target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve.
In a third aspect of the embodiments of the present invention, there is provided a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the above method when executing the computer program.
In a fourth aspect of the embodiments of the present invention, a computer-readable storage medium is provided, which stores a computer program, and the computer program realizes the steps of the above method when being executed by a processor.
Compared with the prior art, the embodiment of the invention has the following beneficial effects: in the embodiment, a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve are obtained; according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve; determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve; and determining the phenological period of the target crop according to the time alignment relation information. Therefore, in the embodiment, the objective crop phenological curve of the objective crop and the reference crop phenological curve corresponding to the objective crop phenological curve can be utilized, so that the phenological period of the objective crop can be automatically extracted, the problems of subjectivity, sidedness and limitation caused by manual experience evaluation are avoided, and the phenological period efficiency and accuracy of the objective crop are improved; in this way, the phenological period of the target crop determined by using the target crop phenological curve of the target crop and the reference crop phenological curve corresponding to the target crop phenological curve is more accurate and has a higher reference value, and the determining efficiency of the phenological period of the target crop is also improved, so that the user experience is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed for the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a flow chart of a method for determining the phenological period of a crop according to an embodiment of the present invention;
FIG. 2 is a scenario diagram of an application scenario of an embodiment of the present invention;
FIG. 3 is a first schematic diagram of a crop phenology curve provided by an embodiment of the present invention;
FIG. 4 is a second schematic diagram of a crop phenology curve provided by an embodiment of the present invention;
FIG. 5 is a third schematic view of a crop phenology curve provided by an embodiment of the present invention;
FIG. 6 is a block diagram of a phenological period determination apparatus for crops according to an embodiment of the present invention;
fig. 7 is a schematic diagram of a computer device provided by an embodiment of the invention.
Reference numerals:
101, a terminal device; 301, target crop phenological curve; 302, referring to a crop phenological curve; 501, referring to a crop phenological curve; 502, target crop phenological curve.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
Hereinafter, a method and an apparatus for determining a phenological period of a crop according to an embodiment of the present invention will be described in detail with reference to the accompanying drawings.
In the prior art, the traditional methods for extracting the crop phenological information mainly adopt manual recording, meteorological station observation and the like, and the methods are often accompanied with the defects of high cost, long time consumption and the like, so that the efficiency of extracting the crop phenological information is low, and the methods are not suitable for acquiring the phenological information of crops in a large area. Meanwhile, due to the influence of latitude and climate factors, it is difficult to directly extend the phenological period of a single-point crop into a large-area. Therefore, a new method for determining the phenological period of crops is needed.
To solve the above problems. The invention provides a method for determining the phenological period of crops, which comprises the steps of obtaining a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve; according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve; determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve; and determining the phenological period of the target crop according to the time alignment relation information. Therefore, in the embodiment, the objective crop phenological curve of the objective crop and the reference crop phenological curve corresponding to the objective crop phenological curve can be utilized, so that the phenological period of the objective crop can be automatically extracted, the problems of subjectivity, sidedness and limitation caused by manual experience evaluation are avoided, and the phenological period efficiency and accuracy of the objective crop are improved; in this way, the phenological period of the target crop determined by using the target crop phenological curve of the target crop and the reference crop phenological curve corresponding to the target crop phenological curve is more accurate and has a higher reference value, and the determining efficiency of the phenological period of the target crop is also improved, so that the user experience is improved.
For example, the embodiment of the present invention may be applied to an application scenario as shown in fig. 2. In this scenario, terminal device 101 may be included.
The terminal device 1 may be hardware or software. When the terminal device 1 is hardware, it may be various electronic devices having a display screen and supporting data processing, including but not limited to a smart phone, a tablet computer, a laptop portable computer, a desktop computer, and the like; when the terminal device 1 is software, it may be installed in an electronic device as described above. The terminal device 1 may be implemented as a plurality of software or software modules, or may be implemented as a single software or software module, which is not limited in this embodiment of the present invention. Of course, the embodiment of the present invention may also be applied to a server.
Specifically, in the application scenario shown in fig. 2, the user may input the target crop through the terminal device 1. After the terminal device 1 acquires the target crop, the terminal device 1 may first acquire a target crop phenological curve of the target crop and a reference crop phenological curve corresponding to the target crop phenological curve; then, according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve; then, the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve can be determined; finally, the phenological period of the target crop can be determined according to the time alignment relation information. The method can realize automatic extraction of the phenological period of the target crop, and avoid problems of subjectivity, sidedness and limitation caused by manual evaluation according to experience, thereby improving the efficiency and accuracy of the phenological period of the target crop. Therefore, the phenological period of the target crop determined by the target crop phenological curve of the target crop and the reference crop phenological curve corresponding to the target crop phenological curve can be more accurate and have more reference value, and the determining efficiency of the phenological period of the target crop is improved, so that the user experience is improved.
It should be noted that the specific type, number, and combination of the terminal devices 1 may be adjusted according to the actual requirements of the application scenarios, which is not limited in this embodiment of the present invention.
It should be noted that the above application scenarios are only presented to facilitate understanding of the present invention, and the embodiments of the present invention are not limited in any way in this respect. Rather, embodiments of the present invention may be applied to any scenario where applicable.
Fig. 1 is a flowchart of a method for determining a phenological period of a crop according to an embodiment of the present invention. The method for determining the phenological period of the crop in fig. 1 may be executed by the terminal device (or server) in fig. 2. As shown in fig. 1, the method for determining the phenological period of the crop includes:
s201: and acquiring a target crop phenological curve of the target crop and a reference crop phenological curve corresponding to the target crop phenological curve.
In this embodiment, the crop whose phenological period needs to be determined may be referred to as a target crop, and may be, for example, winter wheat, corn, and the like. In this embodiment, the crop phenology curve may be a curve capable of reflecting phenology information of the crop, that is, a curve capable of reflecting the response of the growth, development, activity rule and abiotic change of the crop to the joint, for example, the crop phenology curve may be an NDVI (Normalized Difference Vegetation Index) curve, where the abscissa of the curve is days and the ordinate is NDVI value. For convenience of description, a curve formed according to the NDVI values of the collected target crop over a period of time (e.g., one year, 190 days) may be referred to as a target crop phenological curve; for example, because the NDVI value of the target crop in a certain area over a period of time can be obtained by using satellite data, and a target crop phenological curve is generated by using the NDVI value, a daily vegetation index time series can be reconstructed by polynomial fitting; wherein the satellite data has a certain time interval, and therefore an n +1 order polynomial is required to fit the vegetation index time series, where n is the satellite data time interval, and if n >10, the monitoring significance of the satellite in the phenology is greatly reduced, so in one implementation of the present application, the target crop phenology curve does not include a crop phenology curve generated according to satellite data with a revisitation period greater than 10 days (i.e., does not include a revisitation period greater than the NDVI value of the collected target crop over a period of time (e.g., one year, 190 days)). In one implementation, the reference crop phenology curve to which the target crop phenology curve corresponds may be an average crop phenology curve of all target crop phenology curves.
As an example, after determining the target crop and the region of the target crop, a plurality of target crop phenological curves of the target crop may be obtained, for example, a plurality of target crop phenological curves of the target crop may be generated according to phenological parameters of the target crop in the region. Then, an average crop phenological curve of the target crop may be generated according to the multiple target crop phenological curves of the target crop, for example, phenological parameters of the target crop on each day may be calculated, and then an average crop phenological curve of the target crop may be generated according to the phenological parameters of the target crop on each day. Also, the average crop phenology curve of the target crop may be used as a reference crop phenology curve of the target crop.
S202: and according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve.
Because factors such as temperature, longitude and latitude affect the vegetation index value of crops, even in the same phenological period, the vegetation index value may be different in different regions, for example, in some regions, when winter wheat enters the heading period, the NDVI reaches 0.8, and in other regions, the NDVI only reaches 0.7. Therefore, it is necessary to adjust the vegetation index time series curve to reduce the effect of different amplitudes on the two climatic curves. Specifically, the waveform of each target crop phenological curve may be adjusted according to the reference crop phenological curve to obtain an adjusted target crop phenological curve, for example, by stretching the amplitude of the target crop phenological curve, the maximum value (i.e., maximum value) and the minimum value (i.e., minimum value) of the target crop phenological curve may be the same as the maximum value and the minimum value on the reference crop phenological curve.
As an example, a first minimum and a second minimum of the reference crop phenology curve are determined from a maximum of the reference crop phenology curve. For example, a maximum value (i.e., a point at which the NDVI value is the greatest) may be determined in the reference crop phenology curve, and the reference crop phenology curve may be divided into two parts according to the point at which the maximum value is determined, wherein the minimum value of NDVI in the front part of the curve may be used as the first minimum value, and the minimum value of NDVI in the rear part of the curve may be used as the second minimum value.
And determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve. For example, a maximum value (i.e., a point with the largest NDVI value) may be determined in the target crop phenology curve, and the target crop phenology curve may be divided into two parts according to the maximum value, wherein the minimum value of NDVI in the front part of the curve may be used as the first minimum value, and the minimum value of NDVI in the rear part of the curve may be used as the second minimum value.
Then, the amplitude of the waveform of the target crop phenological curve may be adjusted according to the maximum value, the first minimum value, and the second minimum value of the reference crop phenological curve, and the maximum value, the first minimum value, and the second minimum value of the target crop phenological curve, so as to obtain an adjusted target crop phenological curve. In one implementation, the adjusted target crop phenology curve may be determined using the following formula,
Figure BDA0003277996870000121
wherein the content of the first and second substances,
Figure BDA0003277996870000122
the adjusted target crop phenological curve is obtained; sB1(t) is the front part of the phenological curve of the target crop, SB2(t) is the rear part of the objective crop phenological curve, BmaxIs the maximum value of the objective crop phenological curve, Bmin1Is the first minimum value of the objective crop phenological curve, Bmin2A second minimum value of the objective crop phenological curve; a. themaxIs the maximum of the phenological curve of the reference crop, Amin1Is the first minimum value of the reference crop phenological curve, Amin2Is the second minimum of the reference crop phenological curve.
S203: and determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve.
The time alignment relationship information of the adjusted target crop phenological curve and the reference crop phenological curve can be understood as the corresponding relationship of each coordinate point in the adjusted target crop phenological curve and each coordinate point in the reference crop phenological curve in terms of time. In this embodiment, the time alignment relationship information between the adjusted target crop phenological curve and the reference crop phenological curve may be determined according to the NDVI values of each point in the adjusted target crop phenological curve and the reference crop phenological curve.
As an example, a plurality of distance elements may be determined according to the adjusted target crop phenology curve, the reference crop phenology curve and a preset time weighting factor. When calculating the similarity between the sequences of the adjusted target crop phenological curve and the reference crop phenological curve, it is important to consider not only the numerical values of the matching points, but also the period between the matching points, so as to avoid serious mismatch between the sequences. Assuming adjusted target crop phenological curveThe time sequence is X ═ X (X)1,X2,X3,...,Xm) The time sequence of the reference crop phenological curve is Y ═ (Y)1,Y2,Y3,...,Yn) (ii) a First, a distance element matrix D of two time series is constructedm×nEach distance element in the distance element matrix is Di，j＝wi,j×|Xi-YjL where wi，jIs a temporal weighting factor that conforms to a logistic function.
Wherein the distance element matrix D of two time seriesm×nAs can be shown in the following formula,
Figure BDA0003277996870000131
temporal weighting factor wi,jMay be calculated in a manner of
Figure BDA0003277996870000132
Wherein α represents XiTime and YjThe steepness of the time series of the instants, beta representing the median of the time series, tiRepresents XiTime of day tjRepresents YjThe time of day.
Then, several distance matrices may be determined from the several distance elements. For example, formula C may be utilizedi,j＝Di,j+min[Ci-1,j-1,Ci-1,j,Ci,j-1]And respectively determining a distance matrix corresponding to each distance element.
Then, the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve can be determined according to the distance matrixes. In one implementation, the peak point in the adjusted target crop phenological curve may be aligned with the peak point in the reference crop phenological curve according to the principle of minimum distance, that is, for each peak point in the adjusted target crop phenological curve, the peak point in the reference crop phenological curve with the minimum distance from the peak point is used as the aligned peak point corresponding to the peak point; the valley point in the adjusted target crop phenological curve and the valley point in the reference crop phenological curve may be aligned according to the principle of minimum distance, that is, for each valley point in the adjusted target crop phenological curve, the valley point in the reference crop phenological curve having the minimum distance to the valley point is taken as the aligned valley point corresponding to the valley point; for a point between the adjacent peak point and valley point, the points in the bands corresponding to the peak point and the valley point in the two curves can be aligned according to the minimum distance. For example, as shown in fig. 3, the curve 301 located at the upper side is a target crop climate curve, the curve 302 located at the lower side is a reference crop climate curve, the first valley point 1 of the curve 301 may be aligned with the first valley point 1 with the smallest distance in the curve 302, the first peak point 6 of the curve 301 may be aligned with the first peak point 6 with the smallest distance in the curve 302, there are 1 point between the first valley point 1 and the first peak point 6 of the curve 301, and there are 2 points between the first valley point 1 and the first peak point 6 of the curve 302; the minimum distance between 1 point between the first valley point 1 and the first peak point 6 in the curve 301 and 2 points between the first valley point 1 and the first peak point 6 in the curve 302 is defined as the minimum distance, for example, the point 3 between the first valley point 1 and the first peak point 6 in the curve 301 is closest to the point 3 between the first valley point 1 and the first peak point 6 in the curve 302, therefore, a point 3 between the first valley point 1 and the first peak point 6 in the curve 301 is aligned with a point 3 between the first valley point 1 and the first peak point 6 in the curve 302, and a second point 1 between the first valley point 1 and the first peak point 6 in the curve 302 is aligned with the point 1 with the minimum distance between the first valley point 1 and the first peak point 6 in the curve 301; the other bands establish alignment relations in the same way.
S204: and determining the phenological period of the target crop according to the time alignment relation information.
After determining the time alignment relationship information of the adjusted target crop phenological curve and the reference crop phenological curve, determining the phenological period of the target crop according to the time alignment relationship information.
As an example, the green stage node information, heading stage node information, and maturity stage node information of the reference crop phenological curve may be determined first. It should be noted that the determination rules of the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve may be preset. For example, as shown in fig. 4, the NDVI curve of winter wheat has a maximum value, and the whole curve is divided into two parts, namely a curve rising period and a curve falling period, and the green-turning period (i.e., the green-turning period node information) is defined as the day when the NDVI first reaches 10% of the average NDVI curve rising period, i.e., the green-turning period is 10% · (NDVI)max-NDVImin 1)，NDVImaxAt the maximum of the curve, NDVImin1Is the first minimum (i.e., the NDVI value of the first valley point) of the curve, NDVImin2The second minimum of the curve (i.e., the NDVI value at the second valley point). The heading period (i.e., heading period node information) is defined as the day on which NDVI reaches a maximum value, i.e., heading period — NDVImax(ii) a The maturity period (i.e., maturity node information) is defined as the period after the heading, when NDVI falls to 10% of the curve descent period, such as can be obtained by subtracting NVDImin2 from NDVImax, and the maturity period is 10% · (NDVI)max-NDVImin 2)。
Then, the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the target crop phenological curve can be determined according to the time alignment relationship information, the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the reference crop phenological curve.
Specifically, the alignment information between the node information of the green turning period of the target crop phenological curve and the coordinate point of the target crop phenological curve may be determined according to the time alignment relationship information and the node information of the green turning period of the reference crop phenological curve; and determining the node information (such as the node of the green turning period) of the objective crop phenological curve according to the alignment information.It is understood that if it can be determined that the node of the green-turning period of the reference crop phenological curve is aligned with only one point of the target crop phenological curve according to the alignment information, the point in the target crop phenological curve can be used as the node of the green-turning period of the target crop phenological curve; if a plurality of points of the reference crop phenological curve including the nodes of the green-turning period can be determined to be aligned with one point of the target crop phenological curve according to the alignment information, the point in the target crop phenological curve can be used as the node of the green-turning period in the target crop phenological curve; if it can be determined that the nodes of the green-turning period of the reference crop phenological curve are aligned with the points in the target crop phenological curve according to the alignment information, the center point of the points in the target crop phenological curve can be used as the nodes of the green-turning period in the target crop phenological curve. For example, as shown in fig. 5, the upper curve 501 is a reference crop phenology curve, the lower curve 502 is a target crop phenology curve, and a plurality of points on the reference crop phenology curve are matched with one point on the target crop phenology curve; node GUD of green turning period on reference crop phenological curveRAnd the position of the GUD on the target crop phenological curve (GUD)T) This is the case; GUDRAnd three points on the left can match the GUDTPoint, from an alignment perspective, though GUDTAligned with four points, but only GUDRCan be combined with GUDTAnd matching, so the GUDT can be regarded as the green turning stage node GUD of the objective crop phenological curve.
According to the time alignment relation information and the heading stage node information of the target crop phenological curve, determining alignment information of the heading stage node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; and determining heading stage node information (such as heading stage nodes) of the target crop phenological curve according to the alignment information. It is understood that if it can be determined from the alignment information that the heading stage node of the reference crop phenological curve is aligned with only one point in the target crop phenological curve, that point in the target crop phenological curve can be taken as the heading stage node in the target crop phenological curve; if according to the alignment messageThe information can determine that a plurality of points of the reference crop phenological curve including the heading stage node are aligned with one point in the target crop phenological curve, and then the point in the target crop phenological curve can be used as the heading stage node in the target crop phenological curve; if it can be determined that the heading stage node of the reference crop phenological curve is aligned with the plurality of points in the target crop phenological curve according to the alignment information, the center point of the plurality of points in the target crop phenological curve can be used as the heading stage node in the target crop phenological curve. For example, as shown in fig. 5, the upper curve 501 is a reference crop phenological curve, the lower curve 502 is a target crop phenological curve, and one point on the reference crop phenological curve is matched with one point on the target crop phenological curve; heading date node HD of reference crop phenological curveRAnd the location of the heading stage node HD on the target crop phenological curve (HD)T) This is the case; that is, only one line segment connects HDRAnd HDT(ii) a This indicates that, from the viewpoint of waveform similarity, HDTAnd HDRHave similar characteristics, therefore, can convert HD intoTAnd (4) considering the heading stage node HD of the objective crop phenological curve.
According to the time alignment relation information and the maturity node information of the target crop phenological curve, determining the alignment information of the maturity node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; and determining maturity node information (namely maturity nodes) of the objective crop phenological curve according to the alignment information. It is understood that if it can be determined from the alignment information that the maturity node of the reference crop phenology curve is aligned with only one point of the target crop phenology curve, that point in the target crop phenology curve can be taken as the maturity node in the target crop phenology curve; if a plurality of points of the reference crop phenological curve including the maturity node can be determined to be aligned with one point of the target crop phenological curve according to the alignment information, the point of the target crop phenological curve can be used as the maturity node of the target crop phenological curve; if the reference crop can be determined according to the alignment informationThe maturity node of the phenological curve is aligned with the plurality of points in the target crop phenological curve, and then the central point of the plurality of points in the target crop phenological curve can be used as the maturity node in the target crop phenological curve. For example, as shown in fig. 5, the upper curve 501 is a reference crop phenological curve, the lower curve 502 is a target crop phenological curve, and one point on the reference crop phenological curve is matched with a plurality of points on the target crop phenological curve; maturity node MD of reference crop phenological curveRAnd the position of MD on the target crop phenology curve is the case: MDRAligning with three points on a target crop phenological curve; in this case, the average of the three points is selected as the maturity node MD of the objective crop phenological curveT。
Then, the phenological period of the target crop can be determined according to the node information of the green-turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green-turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve.
As an example, if the target crop corresponds to a plurality of target crop phenological curves, for each target crop phenological curve, difference information between the target crop phenological curve and the reference crop phenological curve may be determined according to the node information of the green-turning period, the node information of the heading period, and the node information of the maturity period of the target crop phenological curve, and the node information of the green-turning period, the node information of the heading period, and the node information of the maturity period of the reference crop phenological curve. For example, the difference information between the target crop phenology curve and the reference crop phenology curve may be calculated by the following formula, where n represents the number of points between the corresponding phenology periods of the reference crop phenology curve and the target crop phenology curve, and T represents the time resolution of the NDVI phenology curve.
Then, the phenological period of the target crop may be determined according to difference information corresponding to each target crop phenological curve (for example, the difference information may include difference information corresponding to each of the node information of the green-turning period, the node information of the heading period, and the node information of the maturity period), and the node information of the green-turning period, the node information of the heading period, and the node information of the maturity period of the reference crop phenological curve. For example, the phenological period of the target crop may be obtained by first calculating an average value of difference information corresponding to each target crop phenological curve, and then adding the average value of the difference information to the node information of the green-turning period, the node information of the heading period, and the node information of the maturity period of the reference crop phenological curve.
Compared with the prior art, the embodiment of the invention has the following beneficial effects: in the embodiment, a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve are obtained; according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve; determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve; and determining the phenological period of the target crop according to the time alignment relation information. Therefore, in the embodiment, the objective crop phenological curve of the objective crop and the reference crop phenological curve corresponding to the objective crop phenological curve can be utilized, so that the phenological period of the objective crop can be automatically extracted, the problems of subjectivity, sidedness and limitation caused by manual experience evaluation are avoided, and the phenological period efficiency and accuracy of the objective crop are improved; in this way, the phenological period of the target crop determined by using the target crop phenological curve of the target crop and the reference crop phenological curve corresponding to the target crop phenological curve is more accurate and has a higher reference value, and the determining efficiency of the phenological period of the target crop is also improved, so that the user experience is improved.
All the above-mentioned optional technical solutions can be combined arbitrarily to form the optional embodiments of the present invention, and are not described herein again.
The following are embodiments of the apparatus of the present invention that may be used to perform embodiments of the method of the present invention. For details which are not disclosed in the embodiments of the apparatus of the present invention, reference is made to the embodiments of the method of the present invention.
Fig. 6 is a schematic view of a phenological period determination apparatus for crops according to an embodiment of the present invention. As shown in fig. 6, the phenological period determination apparatus for crops includes:
a curve acquiring unit 601, configured to acquire a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve;
a curve adjusting unit 602, configured to perform waveform adjustment on the target crop phenological curve according to the reference crop phenological curve, so as to obtain an adjusted target crop phenological curve;
an information determining unit 603, configured to determine time alignment relationship information between the adjusted target crop phenological curve and the reference crop phenological curve;
and a phenological period determining unit 604, configured to determine a phenological period of the target crop according to the time alignment relationship information.
Optionally, the curve obtaining unit 601 is configured to:
acquiring a plurality of target crop phenological curves of the target crops;
generating an average crop phenological curve of the target crop according to a plurality of target crop phenological curves of the target crop;
and taking the average crop phenological curve of the target crop as a reference crop phenological curve of the target crop.
Optionally, the curve adjusting unit 602 is configured to:
determining a first minimum value and a second minimum value of the reference crop phenological curve according to the maximum value of the reference crop phenological curve;
determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve;
and adjusting the amplitude of the waveform of the target crop phenological curve according to the maximum value, the first minimum value and the second minimum value of the reference crop phenological curve and the maximum value, the first minimum value and the second minimum value of the target crop phenological curve to obtain an adjusted target crop phenological curve.
Optionally, the information determining unit 603 is configured to:
determining a plurality of distance elements according to the adjusted target crop phenological curve, the reference crop phenological curve and a preset time weight factor;
determining a plurality of distance matrixes according to the distance elements;
and determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve according to the distance matrixes.
Optionally, the waiting period determining unit 604 is configured to:
determining the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve;
determining the green turning period node information, heading period node information and maturity period node information of the target crop phenological curve according to the time alignment relation information, the green turning period node information, heading period node information and maturity period node information of the reference crop phenological curve;
and determining the phenological period of the target crop according to the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve.
Optionally, the waiting period determining unit 604 is configured to:
according to the time alignment relation information and the green return period node information of the target crop phenological curve, determining alignment information of the green return period node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining the node information of the green turning period of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the heading stage node information of the target crop phenological curve, determining alignment information of the heading stage node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining heading stage node information of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the maturity node information of the target crop phenological curve, determining the alignment information of the maturity node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; and determining the maturity node information of the objective crop phenological curve according to the alignment information.
Optionally, the waiting period determining unit 604 is configured to:
for each target crop phenological curve, determining difference information between the target crop phenological curve and the reference crop phenological curve according to the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the target crop phenological curve, and the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the reference crop phenological curve;
and determining the phenological period of the target crop according to the difference information corresponding to each target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve.
It should be understood that, the sequence numbers of the steps in the foregoing embodiments do not imply an execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
Fig. 7 is a schematic diagram of a computer device 7 according to an embodiment of the present invention. As shown in fig. 7, the computer device 7 of this embodiment includes: a processor 701, a memory 702, and a computer program 703 stored in the memory 702 and executable on the processor 701. The steps in the various method embodiments described above are implemented when the computer program 703 is executed by the processor 701. Alternatively, the processor 701 implements the functions of each module/module in each device embodiment described above when executing the computer program 703.
Illustratively, the computer program 703 may be partitioned into one or more modules/modules, which are stored in the memory 702 and executed by the processor 701 to implement the present invention. One or more modules/modules may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution of the computer program 703 in the computer device 7.
The computer device 7 may be a desktop computer, a notebook, a palm computer, a cloud server, or other computer devices. The computer device 7 may include, but is not limited to, a processor 701 and a memory 702. Those skilled in the art will appreciate that fig. 7 is merely an example of a computer device 7 and does not constitute a limitation of the computer device 7 and may include more or fewer components than shown, or some components may be combined, or different components, e.g., the computer device may also include input output devices, network access devices, buses, etc.
The Processor 701 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The storage 702 may be an internal storage module of the computer device 7, for example, a hard disk or a memory of the computer device 7. The memory 702 may also be an external storage device of the computer device 7, such as a plug-in hard disk provided on the computer device 7, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like. Further, the memory 702 may also include both internal and external memory modules of the computer device 7. The memory 702 is used to store computer programs and other programs and data required by the computer device. The memory 702 may also be used to temporarily store data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-mentioned functional modules and modules are illustrated as examples, and in practical applications, the above-mentioned functional allocation may be performed by different functional modules and modules according to requirements, that is, the internal structure of the apparatus is divided into different functional modules or modules to perform all or part of the above-mentioned functions. In the embodiments, each functional module and each module may be integrated into one processing module, or each module may exist alone physically, or two or more modules are integrated into one module, and the integrated modules may be implemented in a form of hardware or a form of software functional modules. In addition, specific names of the functional modules and modules are only used for distinguishing the functional modules and the modules from each other, and are not used for limiting the protection scope of the invention. The modules and the specific working processes of the modules in the system may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or illustrated in a certain embodiment.
Those of ordinary skill in the art will appreciate that the various illustrative modules and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided by the present invention, it should be understood that the disclosed apparatus/computer device and method may be implemented in other ways. For example, the above-described apparatus/computer device embodiments are merely illustrative, e.g., a division of modules or modules into only one logical division, another division may be present in an actual implementation, multiple modules or components may be combined or integrated into another system, or some features may be omitted, or not implemented. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or modules, and may be in an electrical, mechanical or other form.
Modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical modules, may be located in one place, or may be distributed on a plurality of network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
In addition, functional modules in the embodiments of the present invention may be integrated into one processing module, or each of the modules may exist alone physically, or two or more modules are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode.
The integrated modules/modules, if implemented in the form of software functional modules and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow in the method of the embodiments described above can be realized by the present invention, and the computer program can be stored in a computer readable storage medium to instruct related hardware, and when the computer program is executed by a processor, the steps of the method embodiments described above can be realized. The computer program may comprise computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer readable medium may include: any entity or device capable of carrying computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain suitable additions or additions that may be required in accordance with legislative and patent practices within the jurisdiction, for example, in some jurisdictions, computer readable media may not include electrical carrier signals or telecommunications signals in accordance with legislative and patent practices.
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; such modifications and substitutions do not substantially depart from the spirit and scope of the embodiments of the present invention, and are intended to be included within the scope of the present invention.

Claims (9)

1. A method of determining the phenological period of a crop, the method comprising:
acquiring a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve;
according to the reference crop phenological curve, performing waveform adjustment on the target crop phenological curve to obtain an adjusted target crop phenological curve;
determining time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve;
determining the phenological period of the target crop according to the time alignment relation information;
the waveform adjustment of the target crop phenological curve according to the reference crop phenological curve to obtain an adjusted target crop phenological curve comprises:
determining a first minimum value and a second minimum value of the reference crop phenological curve according to the maximum value of the reference crop phenological curve;
determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve;
and adjusting the amplitude of the waveform of the target crop phenological curve according to the maximum value, the first minimum value and the second minimum value of the reference crop phenological curve and the maximum value, the first minimum value and the second minimum value of the target crop phenological curve to obtain an adjusted target crop phenological curve.
2. The method of claim 1, wherein said obtaining a target crop phenology curve of a target crop corresponding to a reference crop phenology curve of the target crop phenology curve comprises:
acquiring a plurality of target crop phenological curves of the target crops;
generating an average crop phenological curve of the target crop according to a plurality of target crop phenological curves of the target crop;
and taking the average crop phenological curve of the target crop as a reference crop phenological curve of the target crop.
3. The method of claim 1, wherein determining the time alignment relationship information of the adjusted target crop climate curve and the reference crop climate curve comprises:
determining a plurality of distance elements according to the adjusted target crop phenological curve, the reference crop phenological curve and a preset time weight factor;
determining a plurality of distance matrixes according to the distance elements;
and determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve according to the distance matrixes.
4. The method of claim 3, wherein determining the phenological period of the target crop based on the time alignment relationship information comprises:
determining the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve;
determining the green turning period node information, heading period node information and maturity period node information of the target crop phenological curve according to the time alignment relation information, the green turning period node information, heading period node information and maturity period node information of the reference crop phenological curve;
and determining the phenological period of the target crop according to the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve, and the node information of the green turning period, the node information of the heading period and the node information of the maturity period of the reference crop phenological curve.
5. The method of claim 4, wherein determining the green-turning stage node information, heading stage node information, and maturity stage node information of the target crop phenology curve from the time-alignment relationship information, the green-turning stage node information, heading stage node information, and maturity stage node information of the reference crop phenology curve comprises:
according to the time alignment relation information and the green return period node information of the reference crop phenological curve, determining alignment information of the green return period node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; determining the node information of the green turning period of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the heading stage node information of the reference crop phenological curve, determining alignment information of the heading stage node information of the target crop phenological curve and a coordinate point of the target crop phenological curve; determining heading stage node information of the target crop phenological curve according to the alignment information;
according to the time alignment relation information and the maturity node information of the reference crop phenological curve, determining alignment information of the maturity node information of the target crop phenological curve and the coordinate point of the target crop phenological curve; and determining the maturity node information of the objective crop phenological curve according to the alignment information.
6. The method of claim 4, wherein determining the phenological period of the target crop based on the node information of the green-turning period, the node information of the heading period and the node information of the maturity period of the target crop phenological curve and the node information of the green-turning period, the node information of heading period and the node information of maturity period of the reference crop phenological curve comprises:
for each target crop phenological curve, determining difference information between the target crop phenological curve and the reference crop phenological curve according to the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the target crop phenological curve, and the node information of the green stage, the node information of the heading stage and the node information of the maturity stage of the reference crop phenological curve;
and determining the phenological period of the target crop according to the difference information corresponding to each target crop phenological curve and the node information of the green turning period, the node information of the heading period and the node information of the mature period of the reference crop phenological curve.
7. An apparatus for determining the phenological period of a crop, characterized in that it comprises:
the curve acquisition unit is used for acquiring a target crop phenological curve of a target crop and a reference crop phenological curve corresponding to the target crop phenological curve;
the curve adjusting unit is used for carrying out waveform adjustment on the target crop phenological curve according to the reference crop phenological curve to obtain an adjusted target crop phenological curve; the waveform adjustment of the target crop phenological curve according to the reference crop phenological curve to obtain an adjusted target crop phenological curve comprises:
determining a first minimum value and a second minimum value of the reference crop phenological curve according to the maximum value of the reference crop phenological curve;
determining a first minimum value and a second minimum value of the target crop phenological curve according to the maximum value of the target crop phenological curve;
adjusting the amplitude of the waveform of the target crop phenological curve according to the maximum value, the first minimum value and the second minimum value of the reference crop phenological curve and the maximum value, the first minimum value and the second minimum value of the target crop phenological curve to obtain an adjusted target crop phenological curve;
the information determining unit is used for determining the time alignment relation information of the adjusted target crop phenological curve and the reference crop phenological curve;
and the phenological period determining unit is used for determining the phenological period of the target crop according to the time alignment relation information.
8. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any one of claims 1 to 6 when executing the computer program.
9. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 6.
CN202111123393.1A 2021-09-24 2021-09-24 Method and device for determining phenological period of crops Active CN113822360B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111123393.1A CN113822360B (en) 2021-09-24 2021-09-24 Method and device for determining phenological period of crops

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111123393.1A CN113822360B (en) 2021-09-24 2021-09-24 Method and device for determining phenological period of crops

Publications (2)

Publication Number Publication Date
CN113822360A CN113822360A (en) 2021-12-21
CN113822360B true CN113822360B (en) 2022-04-22

Family

ID=78915498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111123393.1A Active CN113822360B (en) 2021-09-24 2021-09-24 Method and device for determining phenological period of crops

Country Status (1)

Country Link
CN (1) CN113822360B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106908415A (en) * 2017-02-24 2017-06-30 郑州大学 A kind of big region crops time of infertility Soil Moisture Monitoring method based on amendment NDVI time serieses
CN108805079A (en) * 2018-06-12 2018-11-13 中国科学院地理科学与资源研究所 The recognition methods of winter wheat and device
CN109115770A (en) * 2018-06-14 2019-01-01 中科禾信遥感科技（苏州）有限公司 A kind of a wide range of crops remote-sensing monitoring method and device
CN109211791A (en) * 2018-08-10 2019-01-15 北京观微科技有限公司 Crop condition monitoring method and system
CN111104858A (en) * 2019-11-19 2020-05-05 浙江大学 Large-scale crop phenological extraction method based on morphological model method
CN112304902A (en) * 2020-11-02 2021-02-02 航天宏图信息技术股份有限公司 Real-time monitoring method and device for crop phenology

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2013101105A4 (en) * 2013-08-20 2013-11-07 Beijing Normal University Method for detecting phenological change from time series of vegetation index
CN106845806A (en) * 2017-01-06 2017-06-13 中国科学院遥感与数字地球研究所 The remote-sensing monitoring method and system of farmland planting state
US10779476B2 (en) * 2018-09-11 2020-09-22 Pollen Systems Corporation Crop management method and apparatus with autonomous vehicles
CN112667708A (en) * 2020-12-24 2021-04-16 河南中原光电测控技术有限公司 Winter wheat phenological period monitoring method and device based on NDVI time sequence
CN113050158B (en) * 2021-03-19 2022-05-31 中国科学院武汉岩土力学研究所 Analysis method, device and equipment for near-field microseismic signal waveform and storage medium

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106908415A (en) * 2017-02-24 2017-06-30 郑州大学 A kind of big region crops time of infertility Soil Moisture Monitoring method based on amendment NDVI time serieses
CN108805079A (en) * 2018-06-12 2018-11-13 中国科学院地理科学与资源研究所 The recognition methods of winter wheat and device
CN109115770A (en) * 2018-06-14 2019-01-01 中科禾信遥感科技（苏州）有限公司 A kind of a wide range of crops remote-sensing monitoring method and device
CN109211791A (en) * 2018-08-10 2019-01-15 北京观微科技有限公司 Crop condition monitoring method and system
CN111104858A (en) * 2019-11-19 2020-05-05 浙江大学 Large-scale crop phenological extraction method based on morphological model method
CN112304902A (en) * 2020-11-02 2021-02-02 航天宏图信息技术股份有限公司 Real-time monitoring method and device for crop phenology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Monitoring Crop Phenology Using NDVI Time Series from Sentinel 2 Satellite Data;Mukesh Singh Boori et al;《2019 5th International Conference on Frontiers of Signal Processing (ICFSP)》;20191223;第62-66页 *
基于时序 Landsat遥感数据的新疆开孔河流域农作物类型识别;汪小钦等;《农业工程学报》;20190831;第35卷(第16期);正文第2.1节 *

Also Published As

Publication number Publication date
CN113822360A (en) 2021-12-21

Similar Documents

Publication Publication Date Title
Chong et al. Monthly composites from Sentinel-1 and Sentinel-2 images for regional major crop mapping with Google Earth Engine
Jiang et al. GPhenoVision: a ground mobile system with multi-modal imaging for field-based high throughput phenotyping of cotton
CN112101256B (en) Garlic crop identification method based on coupling active and passive remote sensing images of cloud platform
CN111126662B (en) Irrigation decision making method, device, server and medium based on big data
CN116152673B (en) Method and device for identifying ratoons by utilizing multi-time-phase high-resolution satellite images
CN110647932B (en) Planting crop structure remote sensing image classification method and device
CN114005039B (en) Crop growth assessment method and device based on remote sensing image and electronic equipment
CN112488297B (en) Neural network pruning method, model generation method and device
CN113378793A (en) Crop yield prediction method and device
CN113743819B (en) Crop yield estimation method, device, electronic equipment and storage medium
CN113822360B (en) Method and device for determining phenological period of crops
Guo et al. Panicle Ratio Network: streamlining rice panicle measurement by deep learning with ultra-high-definition aerial images in the field
CN117456367A (en) Crop planting area identification method, device, equipment and storage medium
CN116862106A (en) Urban construction land updating potential analysis method, device, terminal and medium
CN116523145A (en) Photovoltaic power ultra-short-term prediction method and device, computer equipment and storage medium
CN108171417B (en) Planting task adjusting method, electronic device and storage medium
CN113052433B (en) Crop unit yield estimation method based on key time phase and farmland landscape characteristic parameters
CN113240340B (en) Soybean planting area analysis method, device, equipment and medium based on fuzzy classification
CN111199345B (en) Method for measuring and calculating design rainfall of sponge city and terminal equipment
CN115311336A (en) Image registration method, device and equipment of multiple cameras and storage medium
CN113139563B (en) Optimization method and device for image classification model
Hosokawa et al. Contrasting area and yield responses to extreme climate contributes to climate-resilient rice production in Asia
CN114708495B (en) Multi-source irrigation information fusion decision method and system
CN115690585B (en) Method and system for extracting wheat tillering number based on digital photo
CN116936103B (en) User health prediction management method and system based on homodromous network

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