CN110243406A - Crop Estimation Method, device, electronic equipment and storage medium - Google Patents

Crop Estimation Method, device, electronic equipment and storage medium Download PDF

Info

Publication number
CN110243406A
CN110243406A CN201910547885.XA CN201910547885A CN110243406A CN 110243406 A CN110243406 A CN 110243406A CN 201910547885 A CN201910547885 A CN 201910547885A CN 110243406 A CN110243406 A CN 110243406A
Authority
CN
China
Prior art keywords
phenological period
phenological
crops
period
temperature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910547885.XA
Other languages
Chinese (zh)
Other versions
CN110243406B (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.)
Wuhan Union Space Information Technology Co ltd
Wuhan University WHU
Original Assignee
Wuhan Thought Public Space Mdt Infotech Ltd
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 Wuhan Thought Public Space Mdt Infotech Ltd filed Critical Wuhan Thought Public Space Mdt Infotech Ltd
Priority to CN201910547885.XA priority Critical patent/CN110243406B/en
Publication of CN110243406A publication Critical patent/CN110243406A/en
Application granted granted Critical
Publication of CN110243406B publication Critical patent/CN110243406B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00—Measuring or testing not otherwise provided for
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/188—Vegetation

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application provides a kind of Crop Estimation Method, device, electronic equipment and storage medium, which comprises obtains crops to be assessed in the remote sensing image and meteorological data in each phenological period；The remote sensing image and the meteorological data based on each phenological period determine the first net productivity of the crops in the phenological period；The described first net productivity and default production estimation model based on each phenological period, estimate the yield of the crops.Since crops are different in the growing environment of different phenological, nutrition condition or growing way situation, therefore, crops are different in the remote sensing image spectral information of different phenological, simultaneously, since the meteorological data of different phenological is different to the growth effect of crops, and the meteorological data and remote sensing image data of different phenological and the estimated result of crop yield are closely related, therefore, can be improved the reliability of yield estimated result through the above way.

Description

Crop Estimation Method, device, electronic equipment and storage medium
Technical field
This application involves production technology of crops fields, in particular to a kind of Crop Estimation Method, device, electricity Sub- equipment and storage medium.
Background technique
So far from the 1980s, China has made significant headway to the yield by estimation of crops, unites from traditional agronomy Meter method is assessed to satellite remote sensing, and the yield by estimation from monocrop to various crop is assessed, and from zonule to big region Establish multiple Remote Sensing Yield Estimation systems.With the rapid development of remote sensing technology, Remote Sensing Yield Estimation technology has greatly innovated crops Yield investigation method.
However, usually being calculated using the satellite remote-sensing image number of crops in existing Crop Estimation Method Photosynthetic radiation effective ratio, then the yield based on photosynthetically active radiation than estimating crops with solar radiation quantity, yield are estimated The factor considered in the process is relatively simple, causes the reliability of yield estimated result low.
Apply for content
In consideration of it, the embodiment of the present application is designed to provide a kind of Crop Estimation Method, device, electronic equipment and deposits Storage media, to improve the reliability of crop yield estimated result.
In a first aspect, the embodiment of the present application provides a kind of Crop Estimation Method, which comprises obtain to be assessed Remote sensing image and meteorological data of the crops in each phenological period；The remote sensing image and the meteorology based on each phenological period Data determine the first net productivity of the crops in the phenological period；The described first net productivity based on each phenological period With default production estimation model, the yield of the crops is estimated.
During above-mentioned realization, since crops are in the growing environment in different phenological periods, nutrition condition or growing way feelings Condition is different, and therefore, crops are different in the remote sensing image spectral information of different phenological, simultaneously as the gas of different phenological Image data is different to the growth effect of crops, and the meteorological data of different phenological and remote sensing image data and crop yield Estimated result it is closely related, therefore, remote sensing image and the meteorological data in each phenological period based on each phenological period, determine should The net productivity of the first of the crops in phenological period, the first net productivity and default yield then based on each phenological period estimate mould Type estimates the crop yield, fully considered crops different phenological growing state and meteorologic factor to farming The influence of object growth and development improves the reliability of yield estimated result.
Based in a first aspect, the meteorological data includes: actual temperature, solar radiation quantity in a kind of possible design With the ideal temperature of predetermined suitable for crop growth, the remote sensing image and the meteorological number based on each phenological period According to determining the first net productivity of the crops in the phenological period, comprising: the actual temperature based on the phenological period and should The ideal temperature in phenological period determines the temperature stress coefficient in the phenological period；Determine the phenological period temperature stress coefficient and The product of the predetermined ideal efficiency of light energy utilization is first efficiency of light energy utilization in the phenological period；It is described distant based on the phenological period Feel image, determines the photosynthetically active radiation absorptance in the phenological period；It is the solar radiation quantity based on the phenological period, described photosynthetic Net long wave radiation absorptance and first efficiency of light energy utilization determine the first net productivity of the crops in the phenological period.
During above-mentioned realization, due to the yield and the actual temperature in each phenological period of crops, solar radiation quantity with And the ideal temperature of suitable for crop growth is related, therefore, actual temperature and ideal temperature based on the phenological period determine the phenology The temperature stress coefficient of phase, it is then determined the product of the temperature stress coefficient in the phenological period and the preset efficiency of light energy utilization For first efficiency of light energy utilization in the phenological period, finally, the solar radiation quantity, the photosynthetically active radiation based on the phenological period Absorptance and first efficiency of light energy utilization determine that the first net productivity of the crops in the phenological period, which are abundant The different situation of the temperature of different phenological, influence of the solar radiation quantity to crop growth is considered, yield is then made Estimated result is more objective, and precision is higher.
Based in a first aspect, in a kind of possible design, the actual temperature and the phenological period based on the phenological period The ideal temperature, determine the temperature stress coefficient in the phenological period, comprising: the ideal temperature based on the phenological period and One preset algorithm determines the first temperature stress coefficient in the phenological period；The actual temperature, the phenological period based on the phenological period The ideal temperature and the second preset algorithm, determine the second temperature stress coefficient in the phenological period, wherein described second is default Algorithm is different with first preset algorithm；Determine the first temperature stress coefficient and the second temperature side of body in the phenological period The product for compeling coefficient is the temperature stress coefficient in the phenological period.
During above-mentioned realization, since the ideal temperature in the phenological period can generate the calculated result of temperature stress coefficient It influencing, both the ideal temperature in the phenological period and actual temperature jointly can have an impact the calculated result of temperature stress coefficient, Therefore, the ideal temperature and the first preset algorithm based on the phenological period determine the first temperature stress coefficient in the phenological period； The ideal temperature and the second preset algorithm of the actual temperature, the phenological period based on the phenological period, determine the phenological period Second temperature coerce coefficient, wherein second preset algorithm and first preset algorithm are different；Determine the phenological period The product of the first temperature stress coefficient and second temperature stress coefficient is the temperature stress coefficient in the phenological period, the party Formula considers influence of the temperature to temperature stress coefficient from different aspect, then makes the calculated result of temperature stress coefficient more objective It sees.
Based in a first aspect, in a kind of possible design, the meteorological data includes: actual temperature, precipitation and too Positive amount of radiation, the remote sensing image and the meteorological data based on the phenological period, determines the crops in the phenological period First net productivity, which comprises the precipitation of the actual temperature and the phenological period based on the phenological period, really The water stress factor in the fixed phenological period；Determine the water stress factor and predetermined ideal luminous energy benefit in the phenological period It is second efficiency of light energy utilization in the phenological period with the product of rate；The remote sensing image based on the phenological period, determines the phenological period Photosynthetically active radiation absorptance；The solar radiation quantity, the photosynthetically active radiation absorptance and institute based on the phenological period Second efficiency of light energy utilization is stated, determines the first net productivity of the crops in the phenological period.
During above-mentioned realization, due to the yield and the actual temperature in each phenological period of crops, precipitation and too Positive amount of radiation is related, therefore, firstly, the precipitation of the actual temperature based on the phenological period and the phenological period, determines The water stress factor in the phenological period；It is then determined the water stress factor in the phenological period and predetermined desired light The product of energy utilization rate is second efficiency of light energy utilization in the phenological period；Then, the solar radiation quantity based on the phenological period, institute Photosynthetically active radiation absorptance and second efficiency of light energy utilization are stated, determines the first net production of the crops in the phenological period Power.Which has fully considered the influence of precipitation, solar radiation quantity and the temperature of different phenological to crop growth Different situations then makes yield estimated result more objective.
Based in a first aspect, in a kind of possible design, the remote sensing image based on the phenological period determines the phenology The photosynthetically active radiation absorptance of phase, comprising: image procossing is carried out to the remote sensing image in the phenological period, obtains the phenological period Near infrared band reflected value and red spectral band reflected value；The near infrared band reflected value based on the phenological period and described red Optical band reflected value determines the vegetation index in the phenological period；The vegetation index based on the phenological period, determines the phenological period Photosynthetically active radiation absorptance.
During above-mentioned realization, since the photosynthetically active radiation absorptance in each phenological period is close red with each phenological period Wave section reflected value is related with red spectral band reflected value, therefore, is carried out at image by the remote sensing image to the phenological period Reason, obtains the near infrared band reflected value and red spectral band reflected value in the phenological period, then described close red based on the phenological period Wave section reflected value and the red spectral band reflected value, determine the photosynthetically active radiation absorptance in the phenological period.Which is being counted When calculating the photosynthetically active radiation absorptance in each phenological period, the near infrared band reflected value in each phenological period and red has been fully considered The different situation of optical band reflected value, so that the calculated result of photosynthetically active radiation absorptance is more objective.
Based in a first aspect, the default production estimation model includes: preset is somebody's turn to do in a kind of possible design Ratio coefficient, the preset crops in the dry matter of crops between carbon element content and the amount of dry matter are in The ratio coefficient of biomass and total biomass on ground, the preset water in the grain of the storage period crops The ratio coefficient for dividing content and the grain yield of the crops, the described first net productivity and default production based on each phenological period Appraising model is measured, the yield of the crops is estimated, comprising: determine the described first net productivity in all phenological periods and be the Two net productivity；Based on the described second net productivity, crop specie and default production estimation model, the production of crop is determined Amount.
During above-mentioned realization, it is contemplated that the type of production estimation model and crops, the crops dry matter in The biomass and total biomass on ground of ratio coefficient, the crops between carbon element content and the amount of dry matter Ratio coefficient, storage the period crops grain in moisture content and the crops grain yield ratio coefficient Relationship, therefore, so that yield estimated result accuracy is higher.
Based in a first aspect, in a kind of possible design, the described first net productivity based on each phenological period and pre- If production estimation model, the yield of the crops is estimated, comprising: be based on harvesting mode and crop specie, determine the first receipts Obtain coefficient；Mould is estimated based on first coefficient of harvest, the described first net productivity in each phenological period and the default yield Type estimates the crop yield.
It is related with harvesting mode and crop specie due to considering production estimation model during above-mentioned realization, because This, determines the first coefficient of harvest by harvesting mode and crop specie, and based on first coefficient of harvest, the first net production Power and the default production estimation model, estimate crop yield, then make the accuracy of yield estimation result higher.
Second aspect, the embodiment of the present application provide a kind of Crop Estimation device, and described device includes: acquiring unit, use In obtaining crops to be assessed in the remote sensing image in each phenological period and the meteorological data in the phenological period；Determination unit is used for The remote sensing image and the meteorological data based on each phenological period determine that the first of the crops in the phenological period is net raw Force of labor；Estimate unit, for based on each phenological period the described first net productivity and default production estimation model, estimate described The yield of crops.
Based on second aspect, in a kind of possible design, the meteorological data includes: actual temperature, solar radiation quantity With the ideal temperature of predetermined suitable for crop growth, the determination unit is specifically used for the reality based on the phenological period The ideal temperature of temperature and the phenological period determines the temperature stress coefficient in the phenological period；For determining the temperature in the phenological period The product of degree stress coefficient and the predetermined ideal efficiency of light energy utilization is first efficiency of light energy utilization in the phenological period；For being based on The remote sensing image in the phenological period determines the photosynthetically active radiation absorptance in the phenological period；And for being based on the phenological period The solar radiation quantity, the photosynthetically active radiation absorptance and first efficiency of light energy utilization, determine the institute in the phenological period State the first net productivity of crops.
Based on second aspect, in a kind of possible design, the determination unit is also used to based on described in the phenological period Ideal temperature and the first preset algorithm determine the first temperature stress coefficient in the phenological period；And based on described in the phenological period Actual temperature, the ideal temperature in the phenological period and the second preset algorithm determine the second temperature stress coefficient in the phenological period, Wherein, second preset algorithm is different with first preset algorithm；And determine first temperature side of body in the phenological period The product for compeling coefficient and second temperature stress coefficient is the temperature stress coefficient in the phenological period.
Based on second aspect, in a kind of possible design, the meteorological data includes: actual temperature, precipitation and too Positive amount of radiation, the determination unit are also used to the actual temperature based on the phenological period and the precipitation in the phenological period, Determine the water stress factor in the phenological period；For determining the water stress factor and predetermined ideal in the phenological period The product of the efficiency of light energy utilization is second efficiency of light energy utilization in the phenological period；For the remote sensing image based on the phenological period, really The photosynthetically active radiation absorptance in the fixed phenological period；And for the solar radiation quantity, described photosynthetic based on the phenological period Net long wave radiation absorptance and second efficiency of light energy utilization determine the first net productivity of the crops in the phenological period.
Based on second aspect, in a kind of possible design, the determination unit is also used to described distant to the phenological period Feel image and carry out image procossing, obtains the near infrared band reflected value and red spectral band reflected value in the phenological period；And it is based on being somebody's turn to do The near infrared band reflected value and the red spectral band reflected value in phenological period, determine the vegetation index in the phenological period；And The vegetation index based on the phenological period determines the photosynthetically active radiation absorptance in the phenological period.
Based on second aspect, in a kind of possible design, the default production estimation model includes: preset is somebody's turn to do Ratio coefficient, the preset crops in the dry matter of crops between carbon element content and the amount of dry matter are in The ratio coefficient of biomass and total biomass on ground, the preset water in the grain of the storage period crops Divide the ratio coefficient of content and the grain yield of the crops, it is described to estimate unit, it is also used to determine the described of all phenological periods First net productivity and be the second net productivity；Estimated based on the described second net productivity, crop specie and default yield Model is calculated, determines the yield of crop.
It is described to estimate unit in a kind of possible design based on second aspect, it is also used to based on harvesting mode and farming Species determine the first coefficient of harvest；And the described first net productivity based on first coefficient of harvest, each phenological period With the default production estimation model, the lane crop yield is estimated.
The third aspect, the embodiment of the present application provides a kind of electronic equipment, including processor and is connected to the processor Memory, computer program is stored in the memory, when the computer program is executed by the processor, so that institute It states electronic equipment and executes method described in first aspect.
Fourth aspect, the embodiment of the present application provide a kind of storage medium, are stored with computer program in the storage medium, When the computer program is run on computers, so that the computer executes method described in first aspect.
Other feature and advantage of the application will be illustrated in subsequent specification, also, partly be become from specification It is clear that by implementing the embodiment of the present application understanding.The purpose of the application and other advantages can be by written theorys Specifically noted structure is achieved and obtained in bright book, claims and attached drawing.
Detailed description of the invention
Technical solution in ord to more clearly illustrate embodiments of the present application, below will be to needed in the embodiment attached Figure is briefly described, it should be understood that the following drawings illustrates only some embodiments of the application, therefore is not construed as pair The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this A little attached drawings obtain other relevant attached drawings.
Fig. 1 is the structural schematic diagram of electronic equipment provided by the embodiments of the present application；
Fig. 2 is the flow diagram of Crop Estimation Method provided by the embodiments of the present application；
Fig. 3 is the detailed process schematic diagram of Crop Estimation Method provided by the embodiments of the present application；
Fig. 4 is the structural schematic diagram of Crop Estimation device provided by the embodiments of the present application.
Specific embodiment
Below in conjunction with attached drawing in the embodiment of the present application, technical solutions in the embodiments of the present application is described.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing.Meanwhile the application's In description, term " first ", " second " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Fig. 1 is please referred to, the embodiment of the present application provides the structural schematic diagram of a kind of electronic equipment 100, the electronic equipment 100 It can be PC (personal computer, PC), tablet computer, smart phone, personal digital assistant (personal Digital assistant, PDA) etc..
Electronic equipment 100 may include: memory 102, processing 101, communication interface 103 and communication bus, communication bus For realizing the connection communication of these components.
Memory 102 is for storing crops to be assessed in the remote sensing image in each phenological period, the meteorology in the phenological period Data, default production estimation model and Crop Estimation Method provided by the embodiments of the present application and the corresponding calculation procedure of device The various data such as instruction, wherein memory 102 may be, but not limited to, random access memory (Random Access Memory, RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM), electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Processor 101 is for it is real to execute the application when reading and running the computer program instructions being stored in memory The step of Crop Estimation Method of example offer is provided, it is distant in each phenological period to obtain crops to be assessed from memory Feel image, the meteorological data in the phenological period and default production estimation model, and the remote sensing image based on each phenological period With the meteorological data, the first net productivity of the crops in the phenological period is determined, finally, the institute based on each phenological period The first net productivity and default production estimation model are stated, the yield of the crops is estimated.
Wherein, processor 101 may be a kind of IC chip, the processing capacity with signal.Above-mentioned processor 101 can be general processor, including central processing unit (Central Processing Unit, CPU), network processing unit (Network Processor, NP) etc.；It can also be digital signal processor (DSP), specific integrated circuit (ASIC), scene Programmable gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hardware group Part.It may be implemented or execute disclosed each method, step and the logic diagram in the embodiment of the present application.General processor can be with It is that microprocessor or the processor are also possible to any conventional processor etc..
The device of any transceiver one kind can be used in communication interface 103, by the yield estimation results of crops send to The user terminal that electronic equipment 100 communicates to connect is shown.
Referring to figure 2., Fig. 2 is the flow chart of Crop Estimation Method provided by the embodiments of the present application, the method application In electronic equipment 100 as shown in Figure 1, process shown in Fig. 2 will be described in detail below, which comprises
S100: crops to be assessed are obtained in the remote sensing image in each phenological period and the meteorological data in the phenological period.
S200: the remote sensing image and the meteorological data based on each phenological period determine the agriculture in the phenological period The net productivity of the first of crop.
S300: the described first net productivity and default production estimation model based on each phenological period estimate the farming The yield of object.
Wherein, the phenological period refers to that the rule such as growth, development, activity of animals and plants reacts section time with the variation of biology, The phenological period is cried when generating this reaction.
Since crops to be assessed may not in the remote sensing image in each phenological period and the meteorological data in the phenological period Together, therefore, in a kind of possible embodiment, S100 can be implemented as follows: electronic equipment 100 receives unmanned plane Remote sensing image of the crops to be assessed taken in each phenological period, the remote sensing image in the phenological period is stored to memory 102, electronic equipment 100 obtains crops to be assessed in the meteorological data in each phenological period, by the phenological period from weather bureau Remote sensing image is stored to memory 102, when the yield for needing to treat the yield by estimation crops is estimated, is obtained from memory 102 Take crops to be assessed in the remote sensing image in each phenological period and the meteorological data in the phenological period.Wherein, in the present embodiment, The crops to be assessed are rice, and in other embodiments, the crops to be assessed may be wheat etc..
In alternatively possible embodiment, the remote sensing image of crops to be assessed in each phenological period can be logical Cross the remote sensing image that satellite takes.
Wherein, the resolution ratio of the remote sensing image got by unmanned plane is high by the remote sensing image that satellite is got, So the embodiment of the present application obtains the image data of crops using unmanned plane.
In alternatively possible embodiment, when the yield for needing to treat the yield by estimation crops is estimated, then from gas As office obtains crops to be assessed in the meteorological data in each phenological period.
It should be noted that each phenological period in S100 may include all phenological periods of crops to be assessed, it can also To only include at least two crucial phenological periods.The crucial phenological period is related to crop specie, and is the yield to the crops The phenological period being affected.If the result assessed in this way can be more acurrate including whole phenological periods.If only including crucial phenology Phase, can obtain it is more more accurate than the prior art the yield by estimation as a result, and computation complexity it is also smaller.
Specifically, the calculating for reducing yield estimation method is multiple in order on the basis of guaranteeing yield estimated result accuracy Therefore miscellaneous degree in alternatively possible embodiment, estimates crops wait produce for every kind, can only be existed by the crops The remote sensing image in crucial phenological period and the meteorological data in the phenological period realize accurately estimating to the yield of the crops, such as Rice can only need to obtain data of the rice at tillering stage, incubation period and heading stage, and realization is estimated rice yield, is not necessarily to Rice is obtained again in Seedling Stage, milk stage, dough stage and the data in maturity period.Such as wheat again, can be only small by obtaining Data of the wheat at jointing stage, boot stage and heading stage, realization rice yield is estimated, no longer need to obtain wheat sowing time, Seeding stage, tillering stage, Wintering Period, period of seedling establishment, stand up the phase (biology jointing), florescence, pustulation period and maturity period data.
Rice is described as follows in each phenological period growing state:
Seedling Stage: go out a piece of leaf within every 3-4 days；Seminal root is replaced by secondary root.
Tillering stage: when tiller occurs in rice individual, into tillering stage, until starting Spike development.
Boot stage: the sword-like leave pulvinus of rice is exposed to the first dew point of rice, about 1.2 leaf ages (9 days or so).
Heading stage: the spike of rice of 50% plant exposes leaf sheath.
Milk stage: seed content is full of glume in the middle part of 50% or more spike of rice, when content is creamy.
Dough stage: the dense knot of seed content in the middle part of 50% or more spike of rice, do nothing creamy object when.
Maturity period: every fringe has 90% grain yellow maturity, and it is the maturity period that spike of rice base portion blueness grain is also hard.
Wheat is described as follows in each phenological period growing state:
Sowing time: the date of sowing.
Seeding stage: when 50% seed rough leaf of Quan Tian exposes plumule and grows 2 centimetres of ground, spring wheat, March -4 Month；Winter wheat, or so October.
Tillering stage: when leaf sheath 1.5-2cm is stretched out in 50% plant of Quan Tian, first tiller, spring wheat, about ten days in April；Winter wheat, Or so mid or late October.
Wintering Period: daily mean temperature drops to 2 DEG C or so, the wheat plant substantially dormant date, spring wheat, nothing；Winter Wheat, at the beginning of 1 month by the end of November.
Period of seedling establishment: in second year spring, with the rise of temperature, wheat starts to grow, the leaf that 50% plant Nian Houxin is grown Piece (mostly Winter-Spring handover leaf) stretch out leaf sheath 1-2cm, and crop field from it is dark green become dark green when, spring wheat, nothing；Winter wheat, 2 months Under-March on.
Stand up the phase (biology jointing): wheat seeding is grown up by original sprawl growth, the first elongate leaf of Nian Hou, leaf sheath Significant elongation, the auricle of the first elongation leaf and before year last a piece of leaf auricle away from up to 1.5cm, the first internode of base portion is stretched slightly It is long, spring wheat, about the first tenday period of a month in May；Winter wheat about early April.
Jointing stage (agronomy jointing): it is frangible to pinch wheat base portion with finger by the first internode of stem of wheat 1.5-2cm from the ground It makes sound, spring wheat, about mid or late May；Winter wheat, the middle ten days and the last ten days in April.
Boot stage: plant boot leaf (last a piece of leaf) fully extended (auricle is visible), spring wheat, about early June；Winter is small Wheat, the first tenday period of a month in May.
Heading stage: tassel top or side (not referring to awns), when boot leaf sheath stretches out the half of spike length degree, spring wheat, about 6 Early July the middle ten days and the last ten days-the moon；Winter wheat, early June mid or late May-.
Florescence: the flowers are in blossom puts for 50% plant of Quan Tianyou first, spring wheat, about early July；Winter wheat, early June.
Pustulation period: seed shape is basically completed, and length reaches 3/4ths of maximum value, and thickness increases little, spring wheat, About mid-July；Winter wheat starts to be in the milk mid-June.
Maturity period: 1. dough stage: close to normal, inside is in wax-like, seed aqueous 22%, stem leaf for seed size, color Basic to dry out, it is suitable harvest time that wax ripeness latter stage seed dry weight, which reaches maximum value,.2. full ripe stage: it is normal that seed has had kind Size and color, inside are hardened, and moisture content is down to 20% hereinafter, dry-matter accumulation stops, and spring wheat, about the first tenday period of a month in August ten days, winter are small Wheat, early July.
S200: the remote sensing image and the meteorological data based on each phenological period determine the agriculture in the phenological period The net productivity of the first of crop.
In one possible implementation, the meteorological data include: actual temperature, solar radiation quantity and in advance really When the ideal temperature of fixed suitable for crop growth, corresponding, S200 may include:
A1: the ideal temperature of the actual temperature and the phenological period based on the phenological period determines the phenological period Temperature stress coefficient.
In the present embodiment, the actual temperature in the phenological period can be by calculating the highest temperature daily in the phenological period Degree and the average value of minimum temperature acquire, and ideal temperature of the crops in each phenological period data can obtain based on practical experience It takes.
For example, first according to the daily highest temperature and the lowest temperature, calculating daily temperature when having 15 days in the phenological period Then average value, divided by 15, will obtain the actual temperature in the phenological period after temperature averages summation daily in 15 days.
For example, 15 days highest temperatures and the lowest temperature are summed first when having 15 days in the phenological period, then, summation knot Fruit obtains the actual temperature in the phenological period divided by 30.
As a kind of possible embodiment, the actual temperature in the phenological period can be every in the phenological period by calculating The average value of it multiple groups temperature data acquires, wherein the multiple groups temperature data can be in one day at interval of 1 hour Temperature data, the multiple groups temperature data can be the temperature data in one day at interval of 2 hours.
For example, having 24 groups of temperature datas daily in 15 days, first according to daily temperature when having 15 days in the phenological period Data calculate daily temperature averages, then, divided by 15, will obtain the object after temperature averages summation daily in 15 days The actual temperature of time phase.
For example, when having 15 days in the phenological period, 15 days all temperature datas are summed first, then, summed result divided by 360, obtain the actual temperature in the phenological period.
In a kind of possible embodiment, A1 can be implemented as follows: the reality based on the phenological period Temperature, the ideal temperature in the phenological period and preset algorithm determine the temperature stress coefficient in the phenological period, the pre- imputation Method are as follows: Tε2-Pi(x, t)=1.184/ { 1+exp [0.2 × (Topt-Pi(x,t)-10-Temp_Pi(x,t))]}/{1+exp[0.3 ×(-Topt-Pi(x,t)-10+Temp_Pi(x, t)] }, wherein Tε2-Pi(x, t) is the temperature stress coefficient in i-th of phenological period, Topt-Pi(x, t) is the ideal temperature in i-th of phenological period, Temp_Pi(x, t) is the actual temperature in i-th of phenological period, wherein t For marking year locating for the phenological period and the moon, x indicates crops to be estimated represented by a pixel in remote sensing image Real area, in the present embodiment, a pixel indicate that area is 1 square metre of crops.
As another embodiment, in Temp_Pi(x,t)-Topt-Pi(x, t) > 10C ° or Topt-Pi(x,t)-Temp_Pi When (x, t) > 13C °,Wherein, TT is and the actual temperature average value immediate phenological period in phenological period Actual temperature, for example, the actual temperature data in 3 phenological periods is shared, respectively 15 degree, 20 degree and 30 degree, then, the phenological period Actual temperature average value 21.66, then 20 degree closest with 21.66 degree, therefore, TT is 20 degree.
As alternatively possible embodiment, A1 includes:
The ideal temperature and the first preset algorithm based on the phenological period determine the first temperature stress system in the phenological period Number, wherein first preset algorithm are as follows: Tε1-Pi(x, t)=0.8+0.02 × Topt-Pi(x,t)-0.0005×[Topt-Pi(x, t)]2, wherein Tε1-Pi(x, t) is the first temperature stress coefficient in i-th of phenological period.
As an implementation, in Temp_PiAt (x, t)≤- 10C °, T is enabledε1-Pi(x, t)=0.
The ideal temperature and the second preset algorithm of the actual temperature, the phenological period based on the phenological period determine The second temperature in the phenological period coerces coefficient, wherein second preset algorithm is different with first preset algorithm, wherein Second preset algorithm is identical as the preset algorithm.
The product of the first temperature stress coefficient and second temperature stress coefficient that determine the phenological period is the object The temperature stress coefficient of time phase.
The first temperature stress coefficient in the phenological period and the second temperature are coerced into multiplication, and product conduct The temperature system in the object later period coerces coefficient.
A2: the product of the temperature stress coefficient and the predetermined ideal efficiency of light energy utilization that determine the phenological period is the phenology First efficiency of light energy utilization of phase.
In a kind of possible embodiment, A2 can be implemented as follows, by the temperature stress system in the phenological period Several and predetermined ideal efficiency of light energy utilization εmaxIt is multiplied, and using product as first efficiency of light energy utilization in the phenological period, In, under ideal conditions, the ε of rice and wheatmaxValue be 2.8g/MJ.
A3: the remote sensing image based on the phenological period determines the photosynthetically active radiation absorptance in the phenological period.
As an implementation, A3 includes:
Image procossing is carried out to the remote sensing image in the phenological period, obtain the phenological period near infrared band reflected value and Red spectral band reflected value.
Radiometric calibration, geometric calibration, image joint and Image registration processing are carried out to the remote sensing image in i-th of phenological period Afterwards, the near infrared band reflected value NIR in i-th of phenological period is obtainedi(x, t) and red spectral band reflected value Ri(x,t)。
The near infrared band reflected value and the red spectral band reflected value based on the phenological period, determine the phenological period Vegetation index, the vegetation index based on the phenological period determine the photosynthetically active radiation absorptance in the phenological period.
As an implementation, firstly, passing through the near infrared band reflected value in i-th of phenological period, the feux rouges Wave band reflected value andDetermine that the normalization in i-th of phenological period is planted By index NDVIi(x,t).Pass through NDVIi(x, t) andDetermine the light in i-th of phenological period Close Net long wave radiation absorptanceWherein, NDVIi,minAnd NDVIi,maxIt is related with vegetation pattern, for example, rice and wheat NDVIi,min=0.023, NDVIi,max=0.634, the NDVI of bushesi,min=0.023, NDVIi,max=0.636, it is evergreen wealthy The NDVI of Ye Lini,min=0.023, NDVIi,max=0.676, FPARi,maxAnd FPARi,minIt is unrelated with vegetation pattern, FPARi,min =0.001, FPARi,max=0.95.
As an implementation, anti-by the near infrared band reflected value in i-th of phenological period, the red spectral band Penetrate value andOrDetermine the ratio vegetation index RVI in i-th of phenological periodi (x,t).Pass through RVIi(x, t) andReally The photosynthetically active radiation absorptance in fixed i-th of phenological periodWherein, RVIi,minAnd RVIi,maxHave with vegetation pattern It closes, for example, the RVI of rice and wheati,min=1.05, RVIi,max=4.46, the RVI of bushesi,min=1.05, RVIi,max= 5.17, the RVI of evergreen broadleaf foresti,min=0.95, RVIi,max=5.17.
In order to improve the phenological period photosynthetically active radiation absorptance computational accuracy, therefore, as another embodiment, Pass through the photosynthetically active radiation absorptance in i-th of phenological period AndDetermine the photosynthetically active radiation absorptance in i-th of phenological period FPAR_Pi(x, t), wherein α is correction factor, 0 < α < 1, in the present embodiment, and α=0.7, in other embodiments, α=0.6.
A4: the solar radiation quantity, the photosynthetically active radiation absorptance and first luminous energy based on the phenological period Utilization rate determines the first net productivity of the crops in the phenological period.
Wherein, the solar radiation quantity in the phenological period can be by the total radiation of the moon where the phenological period multiplied by the phenological period Number of days and the ratio of the number of days of this month obtain, for example, the number of days i-th of phenological period is 15 days, the moon where i-th phenological period Number of days when being 30 days, the ratio of the number of days of the number of days and this month in i-th of phenological period is 1/2, the moon where i-th of phenological period When total radiation is SOL (x, t), the solar radiation quantity in i-th of phenological period is
As a kind of possible embodiment, the solar radiation quantity in the phenological period can be by calculating the sun in the phenological period Amount of radiation and obtain.
In a kind of possible embodiment, A4 can be implemented as follows, pass through the sun spoke in i-th of phenological period The amount of penetrating andDetermine the photosynthetically active radiation PAR_P in i-th of phenological periodi(x, t) is based on i-th The product of the photosynthetically active radiation absorptance in the photosynthetically active radiation in phenological period and i-th of phenological period determines i-th phenological period The photosynthetically active radiation APAR_P of actual absorptioni(x, t), the photosynthetically active radiation of the actual absorption based on i-th of phenological period APAR_PiThe product of first efficiency of light energy utilization in (x, t) and i-th of phenological period obtains the crops in i-th of phenological period First net productivity NPP_Pi(x,t)。
In a kind of possible embodiment, A4 can be implemented as follows, the sun spoke based on i-th of phenological period The product of the amount of penetrating and the photosynthetically active radiation absorptance in i-th of phenological period, determines that the photosynthetic of the actual absorption in i-th of phenological period has Effect radiation, the photosynthetically active radiation of the actual absorption based on i-th of phenological period and first efficiency of light energy utilization in i-th of phenological period Product determines the first net productivity of the crops in i-th of phenological period.
In alternatively possible implementation, the meteorological data include: actual temperature, solar radiation quantity and in advance When the ideal temperature of determining suitable for crop growth, S200 includes:
B1: the precipitation of the actual temperature and the phenological period based on the phenological period determines the water in the phenological period Divide stress coefficient.
In a kind of possible embodiment, B1 can be implemented in the following way, based on by the reality in i-th of phenological period Temperature and third preset algorithm obtain the local potential steaming amount E in i-th of phenological periodpo-Pi(x, t), wherein in 0C ° < Temp_ PiWhen (x, t) < 26.5C °, the third preset algorithm are as follows: Epo-Pi(x, t)=16 × [10 × Temp_Pi(x,t)/I]β/ 2, In, β=[0.675 × I3-77.1×I2+17920×I+492390]×10-6, wherein when the quantity in phenological period only has 4, I=1, when 2,3,4. quantity in the phenological period only has 3, i=1,2,3. in the present embodiment, the phenology as Crop Estimation Phase is 3,In Temp_PiWhen (x, t) >=26.5C °, local Penman-Monteith formula is only in temperature It rises and increase and, the third preset algorithm unrelated with I value are as follows:In Temp_PiWhen (x, t)≤0C °, Epo-Pi(x, t)=0. Precipitation Precip_P based on i-th of phenological periodiThe local potential steaming amount E in (x, t), i-th crucial phenological periodpo-Pi(x, T) with the 4th preset algorithm, the surface net radiation amount R in i-th of phenological period is obtainedn-Pi(x, t), wherein the 4th preset algorithm Are as follows: Rn-Pi(x, t)=[Epo-Pi(x,t)×Precip_Pi(x,t)]0.5×{0.369+0.598×[Epo-Pi(x,t)/ Precip_Pi(x,t)]0.5}.Precipitation E based on i-th of phenological periodpo-PiThe surface net radiation amount in (x, t), i-th phenological period Rn-Pi(x, t) and the 5th preset algorithm obtains the region actual evapotranspiration EET_P in i-th of phenological periodi(x, t), wherein described 5th preset algorithm are as follows:? Temp_PiWhen (x, t)≤0C °, EET_Pi(x, t)=0.Local potential steaming amount E based on i-th of phenological periodpo-Pi(x,t)、 The region actual evapotranspiration EET_P in i-th of phenological periodiIt is potential to obtain i-th of phenology term area for (x, t) and the 6th preset algorithm Evapotranspiration PET_Pi(x, t), wherein the 6th preset algorithm are as follows: PET_Pi(x, t)=[Epo-Pi(x,t)+EET_Pi(x, t)]/2.Based on i-th of phenological period Regional potential evapotranspiration amount PET_Pi(x, t), i-th of phenology term area actual evapotranspiration EET_Pi (x, t) and the 7th preset algorithm simultaneously calculate, and obtain the water stress factor W in i-th of phenological periodε-Pi(x,t)。Wε-Pi(x, t) is got over Greatly, indicate that ground is more wet.Wherein, the 7th preset algorithm are as follows: Wε-Pi(x, t)=0.5+0.5 × EET_Pi(x,t)/ PET_Pi(x,t)。
As an implementation, in the water stress factor W for calculating i-th of phenological periodε-PiWhen (x, t), in EET_Pi (x,t)≥PET_PiWhen (x, t), EET_P is determinedi(x, t)=PET_Pi(x, t), in EET_Pi(x,t)<PET_PiWhen (x, t), Determine EET_Pi(x, t)=EET_Pi(x,t).Work as Wε-PiWhen (x, t)=0.5, indicate extremely arid, Wε-PiWhen (x, t)=1, table Show very wet.
B2: the product of the water stress factor and the predetermined ideal efficiency of light energy utilization that determine the phenological period is should Second efficiency of light energy utilization in phenological period.
As an implementation, B2 can be implemented as follows, the water stress factor based on the phenological period and Predetermined ideal efficiency of light energy utilization εmaxProduct, determine second efficiency of light energy utilization in the phenological period, wherein in ideal item Under part, the ε of rice and wheatmaxValue be 2.8g/MJ.
Alternatively, B2 can be implemented as follows, the water stress factor based on the phenological period and temperature The product of degree stress coefficient, determines the stress coefficient in the phenological period, based on the stress coefficient and predetermined ideal luminous energy benefit With rate εmaxProduct, determine second efficiency of light energy utilization in the phenological period.
B3: the remote sensing image based on the phenological period determines the photosynthetically active radiation absorptance in the phenological period.Wherein, B3 is identical as A3, please refers to content described in the embodiment of A3, therefore repeat no more.
B4: the solar radiation quantity, the photosynthetically active radiation absorptance and second luminous energy based on the phenological period Utilization rate determines the first net productivity of the crops in the phenological period.
In a kind of possible embodiment, B4 can be implemented as follows, the sun spoke based on i-th of phenological period The amount of penetrating andObtain the photosynthetically active radiation PAR_P in i-th of phenological periodi(x, t) is based on i-th The product of the photosynthetically active radiation absorptance in the photosynthetically active radiation in phenological period and i-th of phenological period obtains i-th phenological period The photosynthetically active radiation APAR_P of actual absorptioni(x, t), finally, photosynthetic effective spoke of the actual absorption based on i-th of phenological period Penetrate APAR_PiThe product of second efficiency of light energy utilization in (x, t) and i-th of phenological period obtains the crops in i-th of phenological period The first net productivity NPP_Pi(x,t)。
In a kind of possible embodiment, B4 can be implemented as follows, the sun spoke based on i-th of phenological period The product of the amount of penetrating and the photosynthetically active radiation absorptance in i-th of phenological period, the photosynthetic of actual absorption for obtaining i-th of phenological period have Effect radiation, the photosynthetically active radiation of the actual absorption based on i-th of phenological period and second efficiency of light energy utilization in i-th of phenological period Product obtains the first net productivity of the crops in i-th of phenological period.
S300: the described first net productivity and default production estimation model based on each phenological period estimate the farming The yield of object.
In one possible implementation, the default production estimation model includes: the preset crops Ratio coefficient, the preset crops in dry matter between carbon element content and the amount of dry matter on ground Biomass and total biomass ratio coefficient, it is preset storage the period crops grain in moisture content with The ratio coefficient of the grain yield of the crops, S300 include:
Determine the described first net productivity in all phenological periods and be the second net productivity NPP (x, t).
For example, treat the yield by estimation crops estimated when, if only needing in the data in all phenological periods of the crops The data in 3 phenological periods when can be achieved with accurately estimating the yield of the crops, it is net raw by described the first of 3 phenological periods Force of labor be added, and will with as the described second net productivity
Treat the yield by estimation crops estimated when, it is assumed that only need the data in 4 phenological periods of the crops, therefore, By the way that described the first of 4 phenological periods the net productivity is added, and will with as the described second net productivity
Based on the described second net productivity, crop specie and default production estimation model, the yield of crop is determined.
Firstly, according to the type of crops, determine in the dry matter of the crops carbon element content and the amount of dry matter it Between ratio coefficient T, the crops the biomass and total biomass on the ground ratio coefficient ρ, in storage period The ratio coefficient ω of the grain yield of moisture content and the crops in the grain of the crops and the harvest of the crops Then coefficient HI based on above-mentioned coefficient and the second net productivity and the default production estimation model, determines the production of crop Measure Y, wherein the default production estimation model are as follows: Y=a × B × HI × 0.001+b, whereinA and b is regression coefficient, and HI is coefficient of harvest, and the coefficient of harvest of different crops is different, harvest The value range of coefficient is 0.3 to 0.8, for example, the coefficient of harvest of rice is 0.45；The coefficient of harvest of wheat is 0.37.A and b Two equations are constructed by historical test data come simultaneous solution, there was only a and b in the two equations is unknown number, different Crops T it is not necessarily identical, for example, the T of rice is 0.38, the T of wheat is 0.39, and different crops ρ is not necessarily identical, For example, the ρ of rice is 0.91, the ρ of wheat is 0.9, and different crops ω is not necessarily identical, for example, the ω of rice is 13.0%, the ω of wheat is 12.5%.
As alternatively possible embodiment, S300 includes:
Due to coefficient of harvest and gather in mode and crops related, the different harvesting modes of different types of crops, Coefficient of harvest is different, therefore, in order to improve the precision of yield estimation results, is based on harvesting mode and crop specie, determining should The first coefficient of harvest HII of crops.
Wherein, when the mode of harvesting is harvester, the first coefficient of harvest of rice is 0.48；First coefficient of harvest of wheat It is 0.4.
When harvesting mode is cut for hand, the first coefficient of harvest of rice is 0.42；First coefficient of harvest of wheat is 0.36.
It is pre- based on the coefficient of harvest, the described first net productivity in each phenological period and the default production estimation model Estimate the crop yield, wherein the default production estimation model are as follows: Y=c [B × HI × 0.001 × 60%+ (a1 × B × HI × 0.001+b1) × 39%]+d, wherein a1, b1, c and d are regression coefficient, and regression coefficient can be by historical data come structure It builds four equations and carrys out simultaneous solution, wherein there was only a1, b1, c and d in this four equations is unknown number, returns system getting After number, it is based on regression coefficient, determines Y=c [B × HI × 0.001 × 60%+ (a1 × B × HI × 0.001+b1) × 39%]+d, Based on the coefficient of harvest, the described first net productivity in each phenological period and the default production estimation model, estimate described Crop yield.Referring to figure 3., the detail flowchart of Crop Estimation Method provided by the embodiments of the present application, wherein this method Calculating process include:
First stage: crops are obtained in the remote sensing image and meteorological data in each phenological period.
Wherein, the remote sensing image in each phenological period includes: multispectral image and RGB image, the meteorological number in each phenological period It is grown according to the precipitation, the actual temperature in the phenological period, the suitable for crop in the preset phenological period that include: the phenological period The solar radiation in ideal temperature and the phenological period；Wherein, the precipitation in the phenological period be the phenological period every intra day ward and, The actual temperature in the phenological period be the phenological period every daily temperature maximum temperature and minimum temperature average value, the phenological period Solar radiation is the sum of the daily solar radiation quantity in the phenological period.
Second stage, remote sensing image and meteorological data based on each phenological period obtain crops in each phenological period Water stress factor, the first temperature stress coefficient, second temperature coerce the ideal light-use of coefficient, predetermined crop Rate, photosynthetically active radiation and photosynthetically active radiation absorptance.
Wherein, the calculating process of the water stress factor in each phenological period is as follows: firstly, the practical temperature based on the phenological period Degree, calculates the local Penman-Monteith formula in the phenological period, then, gross precipitation and local Penman-Monteith formula based on the phenological period, The surface net radiation amount in the phenological period is calculated, then, precipitation and surface net radiation amount based on the phenological period calculate the phenology The region actual evapotranspiration of phase, then, region actual evapotranspiration and local Penman-Monteith formula based on the phenological period calculate the object The Regional potential evapotranspiration amount of time phase, finally, Regional potential evapotranspiration amount and region actual evapotranspiration based on the phenological period, calculating should The water stress factor in phenological period.
Wherein, the calculating process of the first temperature stress coefficient in each phenological period and second temperature stress coefficient is as follows, base In the ideal temperature and the first preset algorithm of the suitable for crop growth in the phenological period, the first temperature stress system in the phenological period is calculated Number；The ideal temperature and the second preset algorithm of the suitable for crop growth of actual temperature, the phenological period based on the phenological period, meter Calculate the second temperature stress coefficient in the phenological period.
Wherein, the calculating process of the photosynthetically active radiation absorptance in each phenological period is as follows: first to each phenological period Multispectral image and RGB image in remote sensing image carry out radiometric calibration, geometric correction, image joint and audio-visual registration, obtain Take the ratio vegetation index RVI in the phenological periodi(x, t) and normalized differential vegetation index NDVIi(x, t) is then based on RVIi(x,t) And NDVIi(x, t) calculates the photosynthetically active radiation absorptance in the phenological period.
Wherein, the photosynthetically active radiation in each phenological period is obtained based on the solar radiation in the phenological period.
Phase III, water stress factor, the first temperature stress coefficient, second temperature stress system based on each phenological period The ideal efficiency of light energy utilization, photosynthetically active radiation and the photosynthetically active radiation absorptance of several, predetermined crop obtain farming First net productivity of the object in the phenological period.
Wherein, the calculating process of the first of each phenological period the net productivity is as follows: firstly, the moisture based on each phenological period The ideal efficiency of light energy utilization of coefficient, the first temperature stress coefficient, second temperature stress coefficient and predetermined crop is coerced, really Then the practical efficiency of light energy utilization in the fixed phenological period is based on photosynthetically active radiation and photosynthetically active radiation absorptance, determines the object The actual absorption photosynthetically active radiation of time phase, finally, the practical efficiency of light energy utilization based on the phenological period and actual absorption is photosynthetic has Effect radiation, determines the first net productivity in the phenological period.
Fourth stage, the described first net productivity and default production estimation model based on each phenological period are estimated described The yield of crops, wherein the default production estimation model includes: carbon in the dry matter of the preset crops The biomass of ratio coefficient, the preset crops between content and the amount of dry matter being on ground and total life The grain of the ratio coefficient of object amount, the preset moisture content in the grain of the storage period crops and the crops The coefficient of harvest of the ratio coefficient of yield and the preset crops.
Wherein, the calculating process of the yield of the crops is as follows, firstly, determining described the first of all phenological periods Net productivity and be the second net productivity；Then, estimated based on the described second net productivity, crop specie and default yield Model is calculated, determines the yield of crop.
Referring to figure 4., Fig. 4 is a kind of structural block diagram of Crop Estimation device 400 provided by the embodiments of the present application.Below Structural block diagram shown in Fig. 4 will be illustrated, shown device includes:
Acquiring unit 410, for obtaining the remote sensing image and the phenological period of crops to be assessed in each phenological period Meteorological data.
Determination unit 420, for based on each phenological period the remote sensing image and the meteorological data, determine the phenology The first net productivity of the crops of phase.
Estimate unit 430, for based on each phenological period the described first net productivity and default production estimation model, in advance Estimate the yield of the crops.
The meteorological data includes: the ideal temperature of actual temperature, solar radiation quantity and the growth of predetermined suitable for crop Degree, determination unit 420, for determining the temperature stress coefficient in the phenological period；For determining the temperature stress coefficient in the phenological period Product with the predetermined ideal efficiency of light energy utilization is first efficiency of light energy utilization in the phenological period；For based on the phenological period The remote sensing image determines the photosynthetically active radiation absorptance in the phenological period；And for the sun based on the phenological period Amount of radiation, the photosynthetically active radiation absorptance and first efficiency of light energy utilization, determine the crops in the phenological period First net productivity.
Determination unit 420 is also used to the ideal temperature and the first preset algorithm based on the phenological period, determines the phenology The first temperature stress coefficient of phase；And the actual temperature, the phenological period based on the phenological period the ideal temperature and Second preset algorithm determines the second temperature stress coefficient in the phenological period, wherein second preset algorithm and described first is in advance Imputation method is different；And determine the first temperature stress coefficient in the phenological period and the product of second temperature stress coefficient For the temperature stress coefficient in the phenological period.
The meteorological data includes: actual temperature, precipitation and solar radiation quantity, determination unit 420, for being based on the object The actual temperature of time phase and the precipitation in the phenological period, determine the water stress factor in the phenological period；For determining The product of the water stress factor in the phenological period and the predetermined ideal efficiency of light energy utilization is second light in the phenological period It can utilization rate；For the remote sensing image based on the phenological period, the photosynthetically active radiation absorptance in the phenological period is determined；And For the solar radiation quantity, the photosynthetically active radiation absorptance and second efficiency of light energy utilization based on the phenological period, Determine the first net productivity of the crops in the phenological period.
The determination unit 420 is also used to carry out image procossing to the remote sensing image in the phenological period, obtains the phenology The near infrared band reflected value and red spectral band reflected value of phase；And the near infrared band reflected value based on the phenological period and The red spectral band reflected value, determines the vegetation index in the phenological period；And the vegetation index based on the phenological period, it determines The photosynthetically active radiation absorptance in the phenological period.
The default production estimation model includes: carbon element content and this is dry in the dry matter of the preset crops The ratio of the biomass and total biomass on ground of ratio coefficient, the preset crops between amount of substance The ratio of the grain yield of coefficient, the preset moisture content in the grain of the storage period crops and the crops The coefficient of harvest of coefficient and the preset crops estimates unit 430, is also used to determine described the first of all phenological periods Net productivity and be the second net productivity；Mould is estimated based on the described second net productivity, crop specie and default yield Type determines the yield of crop.
Unit 430 is estimated, is also used to determine the first coefficient of harvest based on harvesting mode and crop specie；And it is based on First coefficient of harvest, the described first net productivity in each phenological period and the default production estimation model are estimated described Crop yield.
Each functional unit of the present embodiment pair realizes the process of respective function, refers to and retouches in above-mentioned embodiment illustrated in fig. 2 The content stated, details are not described herein again.
In addition, it is stored with computer program in the storage medium the embodiment of the present application also provides a kind of storage medium, When the computer program is run on computers, so that the computer executes any one of the application embodiment and is provided Method.
In conclusion Crop Estimation Method, device, electronic equipment and storage medium that each embodiment of the application proposes, The described method includes: obtaining crops to be assessed in the remote sensing image in each phenological period and the meteorological data in the phenological period；Base The remote sensing image and the meteorological data in each phenological period determine the first net production of the crops in the phenological period Power；The described first net productivity and default production estimation model based on each phenological period, estimate the yield of the crops.
Since crops are different in the image data of the growing environment in different phenological periods, nutrition condition or growing way situation, And the meteorological data in each phenological period is different, and therefore, remote sensing image and the meteorology in each phenological period based on each phenological period Data determine the first net productivity of the crops in the phenological period, the first net productivity then based on each phenological period and pre- If production estimation model, the crop yield is estimated, has fully considered that crops are gentle in the growing state of different phenological As factor is on the growth and development of crops influence, the reliability of raising yield estimated result.
In embodiment provided herein, it should be understood that disclosed device and method, it can also be by other Mode realize.The apparatus embodiments described above are merely exemplary, for example, the flow chart and block diagram in attached drawing are shown According to device, the architectural framework in the cards of method and computer program product, function of multiple embodiments of the application And operation.In this regard, each box in flowchart or block diagram can represent one of a module, section or code Point, a part of the module, section or code includes one or more for implementing the specified logical function executable Instruction.It should also be noted that function marked in the box can also be attached to be different from some implementations as replacement The sequence marked in figure occurs.For example, two continuous boxes can actually be basically executed in parallel, they sometimes may be used To execute in the opposite order, this depends on the function involved.It is also noted that each of block diagram and or flow chart The combination of box in box and block diagram and or flow chart can be based on the defined function of execution or the dedicated of movement The device of hardware is realized, or can be realized using a combination of dedicated hardware and computer instructions.
In addition, each functional module in each embodiment of the application can integrate one independent portion of formation together Point, it is also possible to modules individualism, an independent part can also be integrated to form with two or more modules.

Claims (10)

1. a kind of Crop Estimation Method, which is characterized in that the described method includes:
Crops to be assessed are obtained in the remote sensing image and meteorological data in each phenological period；
The remote sensing image and the meteorological data based on each phenological period, determine the first of the crops in the phenological period Net productivity；
The described first net productivity and default production estimation model based on each phenological period, estimate the yield of the crops.
2. the method according to claim 1, wherein the meteorological data includes: actual temperature, solar radiation quantity With the ideal temperature of predetermined suitable for crop growth, the remote sensing image and the meteorological number based on each phenological period According to determining the first net productivity of the crops in the phenological period, comprising:
The actual temperature and the ideal temperature based on the phenological period, determine the temperature stress coefficient in the phenological period；
The product of the temperature stress coefficient and the predetermined ideal efficiency of light energy utilization that determine the phenological period is the of the phenological period One efficiency of light energy utilization；
The remote sensing image based on each phenological period determines the photosynthetically active radiation absorptance in the phenological period；
The solar radiation quantity, the photosynthetically active radiation absorptance and first efficiency of light energy utilization based on the phenological period, Determine the described first net productivity in the phenological period.
3. according to the method described in claim 2, it is characterized in that, the actual temperature and the ideal based on the phenological period Temperature determines the temperature stress coefficient in the phenological period, comprising:
The ideal temperature and the first preset algorithm based on the phenological period determine the first temperature stress coefficient in the phenological period；
The ideal temperature and the second preset algorithm of the actual temperature, the phenological period based on the phenological period, determine the object The second temperature of time phase coerces coefficient, wherein second preset algorithm is different with first preset algorithm；
The product of the first temperature stress coefficient and second temperature stress coefficient that determine the phenological period is the phenological period The temperature stress coefficient.
4. the method according to claim 1, wherein the meteorological data includes: actual temperature, precipitation and too Positive amount of radiation, the remote sensing image and the meteorological data based on the phenological period, determines the crops in the phenological period First net productivity, which comprises
The precipitation of the actual temperature and the phenological period based on the phenological period, determines the water stress system in the phenological period Number；
The product of the water stress factor and the predetermined ideal efficiency of light energy utilization that determine the phenological period is the phenological period Second efficiency of light energy utilization；
The remote sensing image based on the phenological period determines the photosynthetically active radiation absorptance in the phenological period；
The solar radiation quantity, the photosynthetically active radiation absorptance and second efficiency of light energy utilization based on the phenological period, Determine the described first net productivity in the phenological period.
5. method according to claim 2 or 4, which is characterized in that determine the photosynthetically active radiation absorptance in the phenological period, Include:
Image procossing is carried out to the remote sensing image in the phenological period, obtains the near infrared band reflected value and feux rouges in the phenological period Wave band reflected value；
The near infrared band reflected value and the red spectral band reflected value based on the phenological period, determine the vegetation in the phenological period Index；
The vegetation index based on the phenological period determines the photosynthetically active radiation absorptance in the phenological period.
6. the method according to claim 1, wherein the default production estimation model includes: preset The place of ratio coefficient, the preset crops in the dry matter of the crops between carbon element content and the amount of dry matter The ratio coefficient of biomass and total biomass on ground, it is preset storage the period crops grain in The coefficient of harvest of the ratio coefficient and the preset crops of moisture content and the grain yield of the crops, based on every The described first net productivity in a phenological period and default production estimation model, estimate the yield of the crops, comprising:
Determine the described first net productivity in all phenological periods and be the second net productivity；
Based on the described second net productivity, crop specie and default production estimation model, the yield is determined.
7. the method according to claim 1, wherein the described first net productivity based on each phenological period and pre- If production estimation model, the yield of the crops is estimated, comprising:
Based on harvesting mode and crop specie, the first coefficient of harvest is determined；
Based on first coefficient of harvest, the described first net productivity in each phenological period and the default production estimation model, Estimate the yield.
8. a kind of Crop Estimation device, which is characterized in that described device includes:
Acquiring unit, for obtaining crops to be assessed in the remote sensing image in each phenological period and the meteorological number in the phenological period According to；
Determination unit, for based on each phenological period the remote sensing image and the meteorological data, determine the institute in the phenological period State the first net productivity of crops；
Estimate unit, for based on each phenological period the described first net productivity and default production estimation model, estimate described The yield of crops.
9. a kind of electronic equipment, which is characterized in that including memory and processor, computer journey is stored in the memory Sequence instruction, which is characterized in that when the computer program instructions are read and run by the processor, execute such as claim 1- The step of method described in any one of 7.
10. a kind of storage medium, which is characterized in that be stored with computer program instructions, the computer on the storage medium When program instruction is readable by a computer and runs, execute such as the step of method of any of claims 1-7.
CN201910547885.XA 2019-06-21 2019-06-21 Crop yield estimation method and device, electronic equipment and storage medium Active CN110243406B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910547885.XA CN110243406B (en) 2019-06-21 2019-06-21 Crop yield estimation method and device, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910547885.XA CN110243406B (en) 2019-06-21 2019-06-21 Crop yield estimation method and device, electronic equipment and storage medium

Publications (2)

Publication Number Publication Date
CN110243406A true CN110243406A (en) 2019-09-17
CN110243406B CN110243406B (en) 2021-06-15

Family

ID=67889014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910547885.XA Active CN110243406B (en) 2019-06-21 2019-06-21 Crop yield estimation method and device, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN110243406B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112131948A (en) * 2020-08-25 2020-12-25 航天信德智图（北京）科技有限公司 Actual measurement data-based Guangxi camellia oleifera yield estimation method
CN113743819A (en) * 2021-09-15 2021-12-03 二十一世纪空间技术应用股份有限公司 Method and device for crop yield estimation, electronic equipment and storage medium
CN114510528A (en) * 2022-02-15 2022-05-17 平安科技（深圳）有限公司 Crop yield display method, device electronic equipment and storage medium
CN115577866A (en) * 2022-12-09 2023-01-06 中化现代农业有限公司 Method and device for predicting waiting period, electronic equipment and storage medium
CN116579521A (en) * 2023-05-12 2023-08-11 中山大学 Yield prediction time window determining method, device, equipment and readable storage medium
US11763557B1 (en) * 2020-09-22 2023-09-19 Sentera, Inc. Permanent crop and permanent cropland analysis using aerial vehicles

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101858971A (en) * 2010-06-02 2010-10-13 浙江大学 Rice yield remote sensing estimation method based on MODIS data
CN102650587A (en) * 2012-05-11 2012-08-29 中国农业大学 Crop biomass inversion method based on SEBAL-HJ model
CN104732299A (en) * 2015-04-03 2015-06-24 中国农业科学院农业信息研究所 Maize yield combined prediction system and method
CN106295865A (en) * 2016-08-02 2017-01-04 沈阳农业大学 A kind of Forecasting Methodology of rice yield
CN108985588A (en) * 2018-06-28 2018-12-11 中国科学院遥感与数字地球研究所 A kind of crop yield prediction remote sensing estimation method, device and system
CN109063893A (en) * 2018-06-25 2018-12-21 浙江大学 A kind of rice yield estimating and measuring method of dynamic harvest index in conjunction with net primary productivity
CN109272161A (en) * 2018-09-18 2019-01-25 三亚中科遥感研究所 Rice yield estimation method based on dynamic HI
CN109614891A (en) * 2018-11-27 2019-04-12 北京师范大学 Crops recognition methods based on phenology and remote sensing

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101858971A (en) * 2010-06-02 2010-10-13 浙江大学 Rice yield remote sensing estimation method based on MODIS data
CN102650587A (en) * 2012-05-11 2012-08-29 中国农业大学 Crop biomass inversion method based on SEBAL-HJ model
CN104732299A (en) * 2015-04-03 2015-06-24 中国农业科学院农业信息研究所 Maize yield combined prediction system and method
CN106295865A (en) * 2016-08-02 2017-01-04 沈阳农业大学 A kind of Forecasting Methodology of rice yield
CN109063893A (en) * 2018-06-25 2018-12-21 浙江大学 A kind of rice yield estimating and measuring method of dynamic harvest index in conjunction with net primary productivity
CN108985588A (en) * 2018-06-28 2018-12-11 中国科学院遥感与数字地球研究所 A kind of crop yield prediction remote sensing estimation method, device and system
CN109272161A (en) * 2018-09-18 2019-01-25 三亚中科遥感研究所 Rice yield estimation method based on dynamic HI
CN109614891A (en) * 2018-11-27 2019-04-12 北京师范大学 Crops recognition methods based on phenology and remote sensing

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112131948A (en) * 2020-08-25 2020-12-25 航天信德智图（北京）科技有限公司 Actual measurement data-based Guangxi camellia oleifera yield estimation method
US11763557B1 (en) * 2020-09-22 2023-09-19 Sentera, Inc. Permanent crop and permanent cropland analysis using aerial vehicles
CN113743819A (en) * 2021-09-15 2021-12-03 二十一世纪空间技术应用股份有限公司 Method and device for crop yield estimation, electronic equipment and storage medium
CN113743819B (en) * 2021-09-15 2024-03-26 二十一世纪空间技术应用股份有限公司 Crop yield estimation method, device, electronic equipment and storage medium
CN114510528A (en) * 2022-02-15 2022-05-17 平安科技（深圳）有限公司 Crop yield display method, device electronic equipment and storage medium
CN114510528B (en) * 2022-02-15 2023-11-17 平安科技（深圳）有限公司 Crop yield display method, device electronic equipment and storage medium
CN115577866A (en) * 2022-12-09 2023-01-06 中化现代农业有限公司 Method and device for predicting waiting period, electronic equipment and storage medium
CN116579521A (en) * 2023-05-12 2023-08-11 中山大学 Yield prediction time window determining method, device, equipment and readable storage medium
CN116579521B (en) * 2023-05-12 2024-01-19 中山大学 Yield prediction time window determining method, device, equipment and readable storage medium

Also Published As

Publication number Publication date
CN110243406B (en) 2021-06-15

Similar Documents

Publication Publication Date Title
CN110243406A (en) Crop Estimation Method, device, electronic equipment and storage medium
CN109829234B (en) A kind of across scale Dynamic High-accuracy crop condition monitoring and yield estimation method based on high-definition remote sensing data and crop modeling
Zeng et al. A hybrid approach for detecting corn and soybean phenology with time-series MODIS data
CN109614891A (en) Crops recognition methods based on phenology and remote sensing
Kumar et al. Growth and yield response of soybean (Glycine max L.) in relation to temperature, photoperiod and sunshine duration at Anand, Gujarat, India
Zhou et al. Detection of phenology using an improved shape model on time-series vegetation index in wheat
JP7176980B2 (en) Crop management system and crop management method
Coelho et al. Vegetation indices in the prediction of biomass and grain yield of white oat under irrigation levels
CN101595812A (en) A kind of method for acquiring crop harvest index
CN114254964B (en) Rice regional climate quality assessment method and system
Bai et al. Growth simulation and yield prediction for perennial jujube fruit tree by integrating age into the WOFOST model
Diao et al. Hybrid phenology matching model for robust crop phenological retrieval
CN115730523A (en) Near-real-time prediction method for regional scale crop yield based on deep learning
Gertsis et al. Precision agriculture applications in a high density olive grove adapted for mechanical harvesting in Greece
Shen et al. Optimization of Sowing Date, Irrigation, and Nitrogen Management of Summer Maize Using the DSSAT-CERES-Maize Model in the Guanzhong Plain, China
Song et al. Evaluation of Landsat 8 time series image stacks for predicitng yield and yield components of winter wheat
CN114219227A (en) Method and system for precise fertilization decision and plot level display
Tran et al. Response of maize yield under changing climate and production conditions in Vietnam.
MAK et al. Assessing agro-climatological model of summer mungbean
丸山篤志 et al. Modeling of phenological development in winter wheat to estimate the timing of heading and maturity based on daily mean air temperature and photoperiod
Inman-Bamber Some physiological factors affecting the optimum age and season for harvesting sugarcane
Jeuffroy et al. A model to predict seed number per pod from early pod growth rate in pea (Pisum sativum L.)
CN111160151A (en) Sugarcane continuous cropping range extraction method, device and equipment based on radar time sequence image
Kaur et al. Factors affecting paddy pant at different growth stages
Singh et al. Wheat yield prediction using CERES-Wheat model for decision support in agro-advisory

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
CB03 Change of inventor or designer information

Inventor after: Du Zhiqiang

Inventor after: Liu Zijun

Inventor after: Luo Jun

Inventor before: Du Zhiqiang

Inventor before: Liu Zifa

Inventor before: Luo Jun

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20230301

Address after: 430071 Rooms 02 and 03, 6/F, Research and Experiment Complex Building, Hubei Institute of Automation, No. 34, Xiaohongshan Road, Wuchang District, Wuhan, Hubei Province

Patentee after: WUHAN UNION SPACE INFORMATION TECHNOLOGY Co.,Ltd.

Patentee after: WUHAN University

Address before: No. B03, 26 Floor, Building A, Huicheng, Longan Port, National Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province

Patentee before: WUHAN UNION SPACE INFORMATION TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right