CN104134003A - Crop single yield prediction method based on knowledge and data common drive - Google Patents

Crop single yield prediction method based on knowledge and data common drive Download PDF

Info

Publication number
CN104134003A
CN104134003A CN201410371605.1A CN201410371605A CN104134003A CN 104134003 A CN104134003 A CN 104134003A CN 201410371605 A CN201410371605 A CN 201410371605A CN 104134003 A CN104134003 A CN 104134003A
Authority
CN
China
Prior art keywords
crop
data
knowledge
crop yield
model
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
CN201410371605.1A
Other languages
Chinese (zh)
Other versions
CN104134003B (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.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
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 Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201410371605.1A priority Critical patent/CN104134003B/en
Publication of CN104134003A publication Critical patent/CN104134003A/en
Application granted granted Critical
Publication of CN104134003B publication Critical patent/CN104134003B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Breeding Of Plants And Reproduction By Means Of Culturing (AREA)

Abstract

The invention discloses a crop single yield prediction method based on knowledge and data common drive. The method includes the following steps that a crop single yield production potential submodel based on data drive is used, and single yield production potential E of crops is obtained according to input environment data; by means of a crop model based on a process, a crop single yield submodel based on knowledge drive is built, and then the single yield y of the crops is obtained according to the single yield production potential E of the crops and internal mechanism parameters of the crops; by means of area historical environment data and the single yield y of the crops, model parameters theta d and theta k are obtained through identification so that the single yield of the crops in a certain area is predicted according to meteorological data of the area in the future. By means of the method, historical environment data and a crop knowledge model can be used to the largest extent, the model parameters are estimated systematically, potential analysis is more accurate, and the method is easy to use and reliable and assists in greenhouse crop production, environment conditioning and cultivation management.

Description

Based on knowledge and the common crop yield amount Forecasting Methodology driving of data
Technical field
The invention belongs to data processing method and general botany technical field, relate in particular to a kind of based on knowledge and the common crop yield amount Forecasting Methodology driving of data.
Background technology
Chamber crop production, environment conditioning and cultivation management need reliable crop yield amount Forecasting Methodology.As a rule, crop yield amount Forecasting Methodology, according to whether comprising the knowledge relevant to crop growth rule, can be divided into crop yield amount Forecasting Methodology and the crop yield amount Forecasting Methodology based on knowledge driving based on data-driven.
Crop yield amount Forecasting Methodology based on data-driven is considered as black box by the complicated growth behavior of crop, growth course is not done to consider, the structure of its method place one's entire reliance upon history environment data and crop yield amount data (not comprising the knowledge that any crop is relevant), as adopted the crop yield amount Forecasting Methodology based on artificial neural network.This method can be processed by data mode of learning the nonlinear relationship of mechanism the unknown, and can reach higher precision of prediction to the data under certain environment growth.But the method can not the existing crop knowledge of maximum using removes to improve the precision of prediction of crop yield amount, and the complexity of model promotes rapidly when increasing output variable.
Be different from the method based on data-driven, the crop yield amount Forecasting Methodology driving based on knowledge is the model of setting up based on crop growth rule, the dynamic process such as biomass generation, distribution, leaf area formation that comprises crop, as the plant function structural model (GreenLab) coming based on cooperation beween China and France research and development, this model is a kind of general plant growth model, be successfully applied to the industrial crops such as corn, wheat, cucumber, tomato, and can have assisted chamber crop production, environment conditioning and cultivation management.But this model can not maximum using environmental data removes to improve the precision of prediction of crop yield amount.
Be different from traditional crop yield amount Forecasting Methodology, for maximum using crop knowledge and environmental data improve crop yield amount precision of prediction, the present invention proposes a kind of based on knowledge and the common crop yield amount Forecasting Methodology driving of data, can assist chamber crop production, environment conditioning and cultivation management.
Summary of the invention
The object of the present invention is to provide a kind of reliable crop yield amount Forecasting Methodology, and make the method can the existing crop knowledge of maximum using and environmental data, improve crop yield amount precision of prediction, and can assist chamber crop production, environment conditioning and cultivation management.
For achieving the above object, the invention provides a kind of crop yield amount Forecasting Methodology based on knowledge and the common driving of data, the method comprises the following steps:
Step 1, utilizes the crop yield amount production potential submodel based on data-driven, obtains the single rate production potential E of crop according to the environmental data of input;
Step 2, utilizes the crop modeling based on process, builds the crop yield quantum model obtaining based on knowledge driving, then according to the mechanism parameter of the single rate production potential E of described crop and crop inherence, obtains the single rate y of crop;
Step 3, utilizes the single rate y of the crop that region history environment data and described step 2 obtain, and identification obtains the model parameter θ of the described crop yield amount production potential submodel based on data-driven dmodel parameter θ with the described crop yield quantum model driving based on knowledge k, the single rate of this area crops is predicted according to the weather data in following certain region.
What the present invention proposed takes full advantage of based on knowledge and the common crop yield amount Forecasting Methodology driving of data method and the advantage separately of the method based on data-driven driving based on knowledge, can maximum using crop knowledge and environmental data, improve the reliability of crop yield amount precision of prediction.Difference with the prior art of the present invention is mainly reflected in the present invention and adopts the method based on data-driven, the crop yield amount production potential submodel of structure based on data-driven, the crop modeling of employing based on process, the crop yield quantum model that structure drives based on knowledge, two submodels adopt the mode of series connection to be coupled.Therefore, the present invention more flexibly applicating history environmental data builds the crop yield amount production potential model based on data-driven, and in conjunction with crop growth rule (crop knowledge), from mechanism angle, crop yield amount is predicted.Experiment shows, the advantage of the two that the inventive method is integrated, and there is higher crop yield amount precision of prediction.
The present invention utilizes region history environment data and crop yield amount data, picks out model parameter, for the single rate prediction of following this area crops.The present invention can maximum using history environment data and crop knowledge model, making model parameter obtain systematization estimates, Potentials is more accurate, is easy to use, and this is also for chamber crop production, environment conditioning and cultivation management provide method more reliably.
Accompanying drawing explanation
Fig. 1 is the theory diagram that the present invention is based on knowledge and the common crop yield amount Forecasting Methodology driving of data.
Fig. 2 is the block diagram based on the common crop yield amount Forecasting Methodology driving of knowledge and data according to an embodiment of the invention.
Embodiment
For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.
Fig. 1 is the theory diagram that the present invention is based on knowledge and the common crop yield amount Forecasting Methodology driving of data, as shown in Figure 1, said method comprising the steps of:
Step 1, utilizes the crop yield amount production potential submodel based on data-driven, obtains the single rate production potential E of crop according to the environmental data of input;
This step has reflected the impact of climatic environment for crop yield amount.
Wherein, described environmental data includes but not limited to temperature, illumination, CO 2the weather data such as concentration, relative humidity.
Wherein, the described crop yield amount production potential submodel based on data-driven is that the method structure based on data-driven obtains, and it is expressed as:
E=f d(x,θ d),
Wherein, f d() represents the method based on data-driven, and x is the environmental data of input, and E is the production potential of crop yield amount under this environmental baseline, θ dparameter for submodel E.
Step 2, utilizes the crop modeling based on process, builds the crop yield quantum model obtaining based on knowledge driving, then according to the mechanism parameter of the single rate production potential E of described crop and crop inherence, obtains the single rate y of crop;
Wherein, the described crop yield quantum model driving based on knowledge is expressed as:
y=f k(E,θ k),
Wherein, f k() represents the crop modeling based on process, θ kfor the parameter of submodel y, i.e. the mechanism parameter of crop inherence.
Described step 1 and step 2 also can be thought the crop yield amount production potential submodel based on data-driven and the crop yield quantum model series coupled based on knowledge driving, obtain based on knowledge and the common crop yield amount forecast model driving of data, for this model, input environment data, can obtain the single rate of crop.
Step 3, utilizes the single rate y of the crop that region history environment data and described step 2 obtain, and identification obtains the model parameter θ of the described crop yield amount production potential submodel based on data-driven dmodel parameter θ with the described crop yield quantum model driving based on knowledge k, for the single rate of this area crops being predicted according to the weather data in following certain region.
In an embodiment of the present invention, utilize gradient descent method identification to obtain based on knowledge and the common crop yield amount prediction model parameters θ driving of data dand θ k.
Fig. 2 is the block diagram based on the common crop yield amount Forecasting Methodology driving of knowledge and data according to an embodiment of the invention, as shown in Figure 2, in an embodiment of the present invention, based on radial basis function neural network (RBFN), build the crop yield amount production potential submodel based on data-driven, based on plant function structural model (GreenLab) or garden crop universal model (HortiSim), build the crop yield quantum model driving based on knowledge.
Based on radial basis function neural network (RBFN), building the crop yield amount production potential submodel based on data-driven obtaining can be expressed as:
E=f d(x,θ d)=Φ(x)θ d
Wherein, θ dfor the weighting parameter of RBFN, Φ (x)=[φ 1(x), φ 2(x) ..., φ h(x)] be radial basis function, and φ l ( x ) = 1 + ( x - c l σ l ) 2 .
In plant function structural model (GreenLab), the biomass increment Q (i) of i growth cycle is expressed as:
Q ( i ) = E ( i ) · r · Sp { 1 - exp [ - S ( i ) Sp ] } ,
Wherein, E (i) is i growth cycle crop yield amount production potential under this cycle weather environment condition, r and Sp are inherent mechanism parameter of growing, and usually according to actual production data, estimate to set the total leaf area of crop that S (i) is i growth cycle.
For growth cycle i, the biomass table of Different Crop organ o accumulation is shown:
q o ( i ) = Σ j Σ k = 1 j P o f o ( k ) Q ( i - j + k - 1 ) Σ o Σ k = 1 i - j + k - 1 N o ( i - j + k , k ) P o f o ( i - j + k ) ,
Wherein, P ofor the strong parameter in relative storehouse of crop organ o, usually according to crop actual production data, estimate to set N ofor the number of crop organ o, f ofor crop organ o spread function, the biomass increment that Q () is different growth periods, N o() is the number of organ o, and j is the growth age of organ o, and k is subscript, and crop organ o generally includes blade b, internode e and fruit f.GreenLab model belongs to common plant function structural model, as space is limited, about the more details of GreenLab model, repeats no more.
It should be noted that, for different industrial crops, the calculating object of crop yield amount is slightly different.Such as for tomato crop, what we paid close attention to is the output of tamato fruit, and its single rate computing formula is:
y ( i ) = q f ( i ) = Σ j Σ k = 1 j P f f f ( k ) Q ( i - j + k - 1 ) Σ o Σ k = 1 i - j + k - 1 N o ( i - j + k , k ) P o f o ( i - j + k ) ,
For romaine lettuce, what we mainly paid close attention to is the blade of romaine lettuce, and its single rate computing formula is:
y ( i ) = q d ( i ) = Σ j Σ k = 1 j P b f b ( k ) Q ( i - j + k - 1 ) Σ o Σ k = 1 i - j + k - 1 N o ( i - j + k , k ) P o f o ( i - j + k ) .
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of making, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (10)

1. based on knowledge and the common crop yield amount Forecasting Methodology driving of data, it is characterized in that, the method comprises the following steps:
Step 1, utilizes the crop yield amount production potential submodel based on data-driven, obtains the single rate production potential E of crop according to the environmental data of input;
Step 2, utilizes the crop modeling based on process, builds the crop yield quantum model obtaining based on knowledge driving, then according to the mechanism parameter of the single rate production potential E of described crop and crop inherence, obtains the single rate y of crop;
Step 3, utilizes the single rate y of the crop that region history environment data and described step 2 obtain, and identification obtains the model parameter θ of the described crop yield amount production potential submodel based on data-driven dmodel parameter θ with the described crop yield quantum model driving based on knowledge k, the single rate of this area crops is predicted according to the weather data in following certain region.
2. method according to claim 1, is characterized in that, described environmental data includes but not limited to weather data.
3. method according to claim 1, is characterized in that, the described crop yield amount production potential submodel based on data-driven is that the method structure based on data-driven obtains, and it is expressed as:
E=f d(x,θ d),
Wherein, f d() represents the method based on data-driven, and x is the environmental data of input, and E is the production potential of crop yield amount under corresponding environmental baseline, θ dparameter for submodel E.
4. method according to claim 1, is characterized in that, the described crop yield quantum model driving based on knowledge is expressed as:
y=f k(E,θ k),
Wherein, f k() represents the crop modeling based on process, θ kfor the parameter of submodel y, i.e. the mechanism parameter of crop inherence.
5. method according to claim 1, is characterized in that, in described step 3, utilizes gradient descent method identification to obtain the model parameter θ of the crop yield amount production potential submodel based on data-driven dmodel parameter θ with the crop yield quantum model driving based on knowledge k.
6. method according to claim 1, is characterized in that, based on radial basis function neural network (RBFN), builds the crop yield amount production potential submodel based on data-driven.
7. method according to claim 6, is characterized in that, based on radial basis function neural network (RBFN), building the crop yield amount production potential submodel based on data-driven obtaining can be expressed as:
E=f d(x,θ d)=Φ(x)θ d
Wherein, θ dfor the weighting parameter of RBFN, Φ (x)=[φ 1(x), φ 2(x) ..., φ h(x)] be radial basis function, and φ l ( x ) = 1 + ( x - c l σ l ) 2 .
8. method according to claim 1, is characterized in that, based on plant function structural model (GreenLab) or garden crop universal model (HortiSim), builds the crop yield quantum model driving based on knowledge.
9. method according to claim 8, is characterized in that, in described plant function structural model, the biomass increment Q (i) of i growth cycle is expressed as:
Q ( i ) = E ( i ) · r · Sp { 1 - exp [ - S ( i ) Sp ] } ,
Wherein, E (i) is i growth cycle crop yield amount production potential under this cycle weather environment condition, and r and Sp are inherent mechanism parameter of growing, the total leaf area of crop that S (i) is i growth cycle;
For growth cycle i, the biomass table of Different Crop organ o accumulation is shown:
q o ( i ) = Σ j Σ k = 1 j P o f o ( k ) Q ( i - j + k - 1 ) Σ o Σ k = 1 i - j + k - 1 N o ( i - j + k , k ) P o f o ( i - j + k ) ,
Wherein, P ofor the strong parameter in relative storehouse of crop organ o, N ofor the number of crop organ o, f ofor crop organ o spread function, the biomass increment that Q () is different growth periods, N o() is the number of organ o, and j is the growth age of organ o, and k is subscript.
10. method according to claim 9, is characterized in that, described crop organ o includes but not limited to blade b, internode e and fruit f.
CN201410371605.1A 2014-07-30 2014-07-30 The crop yield amount Forecasting Methodology that knowledge based drives jointly with data Active CN104134003B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410371605.1A CN104134003B (en) 2014-07-30 2014-07-30 The crop yield amount Forecasting Methodology that knowledge based drives jointly with data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410371605.1A CN104134003B (en) 2014-07-30 2014-07-30 The crop yield amount Forecasting Methodology that knowledge based drives jointly with data

Publications (2)

Publication Number Publication Date
CN104134003A true CN104134003A (en) 2014-11-05
CN104134003B CN104134003B (en) 2018-01-30

Family

ID=51806678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410371605.1A Active CN104134003B (en) 2014-07-30 2014-07-30 The crop yield amount Forecasting Methodology that knowledge based drives jointly with data

Country Status (1)

Country Link
CN (1) CN104134003B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656451A (en) * 2015-01-21 2015-05-27 中国科学院自动化研究所 Optimal regulating and controlling method based on crop model for environmental factors in closed system
CN104732299A (en) * 2015-04-03 2015-06-24 中国农业科学院农业信息研究所 Maize yield combined prediction system and method
CN105513096A (en) * 2015-11-18 2016-04-20 青岛农业大学 Method for estimating biomass of winter wheat
CN105494033A (en) * 2015-10-30 2016-04-20 青岛智能产业技术研究院 Intelligent water-saving irrigating method based on crop demand
JP2016146046A (en) * 2015-02-06 2016-08-12 株式会社Jsol Predictor, prediction method and program
CN106570768A (en) * 2016-11-02 2017-04-19 江苏大学 Representation method for greenhouse tomato sink strength
CN107368687A (en) * 2017-07-25 2017-11-21 中国农业科学院农业信息研究所 A kind of method for optimizing and device of meteorological yield model
CN112308284A (en) * 2019-08-02 2021-02-02 丰田自动车株式会社 Crop characteristic prediction system and method, and computer readable storage medium
CN113011683A (en) * 2021-04-26 2021-06-22 中国科学院地理科学与资源研究所 Crop yield estimation method and system based on corrected crop model
CN113379155A (en) * 2021-06-29 2021-09-10 哈尔滨工业大学 Method for estimating development suitability of biomass energy based on village and town population prediction

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110495408B (en) * 2019-09-20 2021-08-17 重庆工商大学 Fish, shrimp and ginseng breeding decision system and device based on common driving of data and knowledge

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101480143A (en) * 2009-01-21 2009-07-15 华中科技大学 Method for predicating single yield of crops in irrigated area
US20110313666A1 (en) * 2008-08-08 2011-12-22 Raisio Nutrition Ltd Method for monitoring environmental impacts of crop cultivation
CN103518516A (en) * 2013-10-21 2014-01-22 贵州省辣椒研究所 Method for measuring field yield of chilies

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110313666A1 (en) * 2008-08-08 2011-12-22 Raisio Nutrition Ltd Method for monitoring environmental impacts of crop cultivation
CN101480143A (en) * 2009-01-21 2009-07-15 华中科技大学 Method for predicating single yield of crops in irrigated area
CN103518516A (en) * 2013-10-21 2014-01-22 贵州省辣椒研究所 Method for measuring field yield of chilies

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刁军 等: "结构-功能模型在林木竞争研究中的应用", 《世界林业研究》 *
卢文娟 等: "奇台县粮食生产潜力变化过程分析", 《干旱区资源与环境》 *
杨丽丽 等: "不同植物密度番茄生长行为的结构功能模型模拟", 《农业机械学报》 *
韩泽群 等: "加工番茄产量组合预测模型研究", 《中国农学通报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656451B (en) * 2015-01-21 2017-06-09 中国科学院自动化研究所 A kind of closed system envirment factor optimization regulating method based on crop modeling
CN104656451A (en) * 2015-01-21 2015-05-27 中国科学院自动化研究所 Optimal regulating and controlling method based on crop model for environmental factors in closed system
JP2016146046A (en) * 2015-02-06 2016-08-12 株式会社Jsol Predictor, prediction method and program
CN104732299A (en) * 2015-04-03 2015-06-24 中国农业科学院农业信息研究所 Maize yield combined prediction system and method
CN105494033A (en) * 2015-10-30 2016-04-20 青岛智能产业技术研究院 Intelligent water-saving irrigating method based on crop demand
CN105494033B (en) * 2015-10-30 2018-06-01 青岛智能产业技术研究院 A kind of intelligent water-saving irrigation method based on crop demand
CN105513096A (en) * 2015-11-18 2016-04-20 青岛农业大学 Method for estimating biomass of winter wheat
CN106570768A (en) * 2016-11-02 2017-04-19 江苏大学 Representation method for greenhouse tomato sink strength
CN106570768B (en) * 2016-11-02 2019-08-27 江苏大学 A kind of characterizing method that greenhouse tomato library is strong
CN107368687A (en) * 2017-07-25 2017-11-21 中国农业科学院农业信息研究所 A kind of method for optimizing and device of meteorological yield model
CN107368687B (en) * 2017-07-25 2020-01-21 中国农业科学院农业信息研究所 Optimal selection method and device for meteorological unit production model
CN112308284A (en) * 2019-08-02 2021-02-02 丰田自动车株式会社 Crop characteristic prediction system and method, and computer readable storage medium
CN113011683A (en) * 2021-04-26 2021-06-22 中国科学院地理科学与资源研究所 Crop yield estimation method and system based on corrected crop model
CN113379155A (en) * 2021-06-29 2021-09-10 哈尔滨工业大学 Method for estimating development suitability of biomass energy based on village and town population prediction

Also Published As

Publication number Publication date
CN104134003B (en) 2018-01-30

Similar Documents

Publication Publication Date Title
CN104134003A (en) Crop single yield prediction method based on knowledge and data common drive
US11741555B2 (en) Crop yield estimation method based on deep temporal and spatial feature combined learning
Ariffin et al. Automatic climate control for mushroom cultivation using IoT approach
US20170273258A1 (en) System and method for dynamic irrigation management
Jiayu et al. Application of intelligence information fusion technology in agriculture monitoring and early-warning research
Pandit et al. Accounting for environmental conditions in data-driven wind turbine power models
Thilagavathi Online farming based on embedded systems and wireless sensor networks
Çakır et al. Yield prediction of wheat in south-east region of Turkey by using artificial neural networks
Wen et al. Application of ARIMA and SVM mixed model in agricultural management under the background of intellectual agriculture
Massaoudi et al. A novel approach based deep RNN using hybrid NARX-LSTM model for solar power forecasting
Guillén-Navarro et al. A decision support system for water optimization in anti-frost techniques by sprinklers
Praveena et al. IoT Based Water Showering Mechanism in Agriculture
Vladykin et al. Mathematical model of temperature mode for protected ground
CN105159247A (en) Information transmission based phase affiliation judgment method for real time sampling points in intermittent process
Wanyama et al. A Systematic Review of Fourth Industrial Revolution Technologies in Smart Irrigation: Constraints, Opportunities, and Future Prospects for sub-Saharan Africa
CN103903065A (en) Penicillin production process product quality prediction method based on kernel learning technology
Kempelis et al. Implementation of machine learning based approach in IoT network prototype
Ghareeb et al. Wireless Sensor Network-Based Artificial Intelligent Irrigation System: Challenges and Limitations.
Singh Sustainable and Smart Agriculture: A Holistic Approach
Karna et al. Temperature Prediction using Regression Model
Bakare Machine learning-based smart irrigation system and soil nutrients analysis to increase productivity in agriculture field
Bertelè et al. Non-deterministic wind observation from wind turbine loads
Gujar et al. Intelligent Irrigation System for Agriculture using IoT and Machine Learning
Fontanet et al. Simulating Irrigation Scheduling Scenarios and the Effect of Soil Hydraulic Properties on Water Balance Components and Crop Yield
Cisty et al. Support vector machine based model for water content in soil interpolation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant