JP2022136057A - Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program - Google Patents

Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program Download PDF

Info

Publication number
JP2022136057A
JP2022136057A JP2022033768A JP2022033768A JP2022136057A JP 2022136057 A JP2022136057 A JP 2022136057A JP 2022033768 A JP2022033768 A JP 2022033768A JP 2022033768 A JP2022033768 A JP 2022033768A JP 2022136057 A JP2022136057 A JP 2022136057A
Authority
JP
Japan
Prior art keywords
year
prediction model
data
prediction
production performance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2022033768A
Other languages
Japanese (ja)
Inventor
涼子 森岡
Ryoko Morioka
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.)
National Agriculture and Food Research Organization
Original Assignee
National Agriculture and Food Research Organization
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 National Agriculture and Food Research Organization filed Critical National Agriculture and Food Research Organization
Publication of JP2022136057A publication Critical patent/JP2022136057A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

To implement a technology for predicting production performance of crops.SOLUTION: An apparatus (20) for generating a prediction model for predicting production performance of crops includes a prediction model generation unit (23) which generates a prediction model (24) which predicts production performance of the year from at least a part of actual measurement values of production performance of previous years, the prediction model (24) referring to actual measurement value data of past annual production performance and representing interannual change of production performance.SELECTED DRAWING: Figure 1

Description

新規性喪失の例外適用申請有り There is an application for exception to loss of novelty

本発明は、農作物の生産成績を予測する予測モデルの生成方法、生成装置、及び生成プログラムに関する。 The present invention relates to a prediction model generation method, generation device, and generation program for predicting the production performance of agricultural products.

様々な農作物の栽培現場において、栽培管理、出荷計画、販売交渉等に活用するために、気象データを用いて生育期や収穫期の農作物の成分含有量を予測することが試みられている。 At the cultivation sites of various crops, attempts have been made to predict the ingredient content of crops during the growing and harvesting periods using weather data for use in cultivation management, shipping plans, sales negotiations, and the like.

特許文献1には、作物の糖度を目標糖度にするための日射情報及び気温情報を予測し、予測した日射情報及び気温情報に基づき温室内の環境を制御する技術が記載されている。特許文献1に記載された技術では、作物固有の基準糖度と目標糖度との差を、収穫前4週間の日射及び気温の平均値と基準値との差に関連させて、糖度予測に利用している。 Patent Literature 1 describes a technique for predicting solar radiation information and temperature information for setting the sugar content of crops to a target sugar content, and controlling the environment in a greenhouse based on the predicted solar radiation information and temperature information. In the technique described in Patent Document 1, the difference between the standard sugar content specific to crops and the target sugar content is used for sugar content prediction in association with the difference between the average values of solar radiation and temperature for four weeks before harvest and the standard values. ing.

非特許文献1には、過去数年間の気象変動の時系列データ及びミカンの品質の時系列データを用いて学習したモデルを構築したことが記載されている。非特許文献1に記載された技術では、構築したモデルを用いて、当年の収穫時期までの気象及びミカンの品質の時系列データから当年の収穫時期におけるミカンの品質を予測している。 Non-Patent Document 1 describes that a trained model was constructed using time-series data of climate change and time-series data of mandarin orange quality over the past several years. In the technique described in Non-Patent Document 1, the constructed model is used to predict the quality of mandarin oranges at the harvest time of the current year from time-series data of the weather and the quality of mandarin oranges until the harvest time of the current year.

特開2020-48551号公報JP 2020-48551 A

森本ら、植物環境工学17(2):90-98,2005Morimoto et al., Plant Environmental Engineering 17(2):90-98, 2005

農作物の成分含有量の予測における問題点として、以下の3つが挙げられる。1点目は、データ測定点が少ないことである。生育過程の農作物の成分含有量を高頻度に測定してデータを蓄積する例は稀であり、多くは収穫期又は出荷時にのみ成分含有量が測定されている。2点目は、成分含有量が気象以外の影響も受けることである。例えば、多年生の作物では、経年のゆらぎが成分含有量に影響する場合があり、また、単年性の作物では消費者のし好を反映した栽培管理の経年傾向が成分含有量に影響する場合がある。3点目は、成分含有量に影響する気象データの変数の種類や時期が不明な場合が多いことである。 There are three problems in predicting the component content of agricultural products. The first point is that the number of data measurement points is small. It is rare to frequently measure the ingredient content of agricultural products during the growth process and accumulate data, and in many cases, the ingredient content is measured only at the time of harvest or shipping. The second point is that the component content is also affected by factors other than weather. For example, in perennial crops, fluctuations over time may affect the component content, and in single-year crops, the aging trend of cultivation management that reflects consumer preferences may affect the component content. There is The third point is that in many cases the type and timing of weather data variables that affect component content are unknown.

したがって、過去の成分含有量の測定値と気象データとを用いて、作物生産環境の経年変化を反映した予測モデルが実現できれば、栽培現場に大きく貢献する。 Therefore, if a prediction model that reflects changes in the crop production environment over time can be realized using past measured values of component contents and meteorological data, it will greatly contribute to cultivation sites.

特許文献1に記載された技術では、作物固有の基準糖度が不明な場合には糖度予測に適用することが困難である。また、特許文献1に記載された技術では、収穫日が未定の場合には、予測に利用する日射及び気温のデータの起点が定まらないため、糖度予測に適用することが困難である。 It is difficult to apply the technique described in Patent Document 1 to sugar content prediction when the standard sugar content specific to crops is unknown. In addition, with the technique described in Patent Document 1, when the harvest date is undecided, the starting point of the solar radiation and temperature data used for prediction is not determined, so it is difficult to apply it to sugar content prediction.

非特許文献1に記載された技術では、みかんの品質の時系列データを作成するために、生育過程で複数回成分含有量の測定を行う必要がある。また、非特許文献1に記載されたモデルは、栽培暦上の次の時期の品質予測には利用できるが、より早期に品質を予測することは困難である。 In the technique described in Non-Patent Document 1, it is necessary to measure the component content multiple times during the growth process in order to create time-series data on the quality of mandarin oranges. In addition, the model described in Non-Patent Document 1 can be used for quality prediction for the next season on the cultivation calendar, but it is difficult to predict quality earlier.

本発明の一態様は、上述した問題点を解決するためになされたものであり、その目的は、農作物の生産成績を予測する技術を実現することを目的とする。 One aspect of the present invention has been made to solve the above-described problems, and an object thereof is to realize a technique for predicting the production results of agricultural crops.

本発明の一態様に係る予測モデルの生成方法は、農作物の生産成績を予測する予測モデルの生成方法であって、過去の年毎の前記生産成績の実測値データを参照して、前記生産成績の経年変化を表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部から当年の前記生産成績を予測する予測モデルを生成する工程を包含する。 A method of generating a prediction model according to one aspect of the present invention is a method of generating a prediction model for predicting the production performance of agricultural products, wherein the actual measurement data of the production performance for each past year is referenced, and the production performance and generating a prediction model that predicts the production performance of the current year from at least a portion of the measured production performance of previous years.

本発明の一態様に係る予測モデルの生成装置は、農作物の生産成績を予測する予測モデルの生成装置であって、過去の年毎の前記生産成績の実測値データを参照して、前記生産成績の経年変化を表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部から当年の前記生産成績を予測する予測モデルを生成する生成部を備えている。 A prediction model generation device according to an aspect of the present invention is a prediction model generation device for predicting the production results of agricultural products, wherein the actual measurement data of the production results for each past year is referred to, and The forecasting model representing the secular change of is provided with a generation unit that generates a forecasting model for predicting the production performance of the current year from at least a part of the measured values of the production performance up to the previous year.

本発明の各態様に係る予測モデルの生成装置は、コンピュータによって実現してもよく、この場合には、コンピュータを前記生成装置が備える各部(ソフトウェア要素)として動作させることにより前記生成装置をコンピュータにて実現させる生成装置の制御プログラム、およびそれを記録したコンピュータ読み取り可能な記録媒体も、本発明の範疇に入る。 The predictive model generation device according to each aspect of the present invention may be implemented by a computer. A control program for a generating device realized by a computer and a computer-readable recording medium recording it are also included in the scope of the present invention.

本発明の一態様によれば、農作物の生産成績を予測する技術を実現することができる。 ADVANTAGE OF THE INVENTION According to one aspect of the present invention, it is possible to realize a technique for predicting the production results of agricultural products.

本発明の一態様に係る予測システムの要部構成の一例を示すブロック図である。It is a block diagram showing an example of the important section composition of the prediction system concerning one mode of the present invention. 本発明の一態様に係る予測装置で利用する前年までの年毎の生産成績の実測値データの一例を示す図である。FIG. 5 is a diagram showing an example of measured value data of production results for each year up to the previous year, which is used by the prediction device according to one aspect of the present invention; 農作物の生産成績の経年変化及び気象が生産成績に及ぼす影響について説明する図である。It is a figure explaining the secular change of the production result of agricultural products, and the influence which a weather has on a production result. 本発明の一態様に係る気象データの予測モデル生成装置において利用する、気象データの変数が生産成績に及ぼす影響を説明する図である。FIG. 5 is a diagram for explaining the effects of weather data variables on production results, which are used in the weather data prediction model generation device according to one aspect of the present invention. 本発明の一態様に係る予測システムが実行する、予測モデルを用いた農作物の生産成績の予測の概念を説明する図である。FIG. 2 is a diagram illustrating the concept of predicting the production results of agricultural products using a prediction model, which is executed by the prediction system according to one aspect of the present invention; 本発明の一態様に係る予測装置が実行する予測処理の一例を示すフローチャートである。4 is a flowchart showing an example of prediction processing executed by a prediction device according to one aspect of the present invention; 本発明の一態様に係る予測モデルの生成装置が実行する予測モデル生成処理の一例を示すフローチャートである。4 is a flowchart showing an example of a prediction model generation process executed by a prediction model generation device according to an aspect of the present invention; 本発明の一態様に係る気象データの予測モデル生成装置が実行する予測モデルの生成処理の一例を示すフローチャートである。4 is a flowchart showing an example of a prediction model generation process executed by a weather data prediction model generation device according to an aspect of the present invention. 実施例の糖度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the sugar content prediction accuracy of an Example. 実施例の糖度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the sugar content prediction accuracy of an Example. 実施例の酸度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the acidity prediction accuracy of an Example. 実施例の酸度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the acidity prediction accuracy of an Example. 実施例の糖度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the sugar content prediction accuracy of an Example. 実施例の酸度予測精度の検証結果を示す図である。It is a figure which shows the verification result of the acidity prediction accuracy of an Example. 実施例における、降水量、平均気温、及び日照時間の糖度への影響を示す図である。FIG. 4 is a diagram showing the effects of rainfall, average temperature, and hours of sunshine on sugar content in Examples. 土壌水分値が生産成績に及ぼす影響を説明する図である。It is a figure explaining the influence which a soil water content value has on production results. 本発明の他の態様に係る予測システムの要部構成の一例を示すブロック図である。FIG. 11 is a block diagram showing an example of a main configuration of a prediction system according to another aspect of the present invention;

本発明の一態様は、農作物の生産成績を予測するための技術を実現するものである。本発明の一態様により予測する農作物の生産成績には、農作物の収穫量;糖度、酸度、タンパク含有量、質量分析等で計測される分子成分、色素量等の成分含有量;色、形、大きさ、長さ等の形態;倒伏度、葉・花・茎・果実の成長指数(茎の長さや葉の出開き度、花や実の色づき具合)、着花数、着果数、葉の数が含まれるがこれに限定されない。すなわち、本発明の一態様によれば、生産された農作物を評価する評価項目の少なくとも1つを予測することができる。本発明の一態様により生産成績を予測する農作物は特に限定されず、ミカン、リンゴ、モモ、ブドウ、びわ、梅、コーヒー、緑茶、ビーツ、トウキ、カンゾウ、シャクヤク、ダイオウ等の多年生植物;イネ、麦、大豆、落花生、きび、トマト、ホウレンソウ、ナス、カボチャ、イチゴ、シソ、センキュウ、サイコ、ジャガイモ、サツマイモ、観賞用花き等の一年生植物等が例として挙げられる。 One aspect of the present invention provides a technique for predicting the production performance of agricultural crops. The production results of agricultural crops predicted by one embodiment of the present invention include yield of agricultural crops; sugar content, acidity, protein content, content of components such as molecular components and pigment amounts measured by mass spectrometry, etc.; color, shape, Morphology such as size and length; lodging rate, leaf/flower/stem/fruit growth index (stem length, leaf opening, flower and fruit coloring), number of flowers, number of fruits, leaves including but not limited to the number of That is, according to one aspect of the present invention, it is possible to predict at least one evaluation item for evaluating produced crops. Agricultural crops whose production results are predicted according to one embodiment of the present invention are not particularly limited, and perennial plants such as mandarin oranges, apples, peaches, grapes, loquats, plums, coffee, green tea, beets, angelica, licorice, peony, and rhubarb; Examples include annual plants such as barley, soybean, peanut, millet, tomato, spinach, eggplant, pumpkin, strawberry, perilla, cnidium, rhizome, potato, sweet potato, and ornamental flowers.

〔予測システム100〕
図1に基づいて、農作物の生産成績を予測するために用いる予測システム100について説明する。予測システム100は、予測装置10、予測モデル生成装置20、及び気象データの予測モデル生成装置30を備えている。また、予測システム100は、さらに、入力装置40、記憶装置50、及び出力装置60を備えている。予測システム100は、予測装置10、予測モデル生成装置20、及び予測モデル生成装置30をそれぞれ独立した装置として備えていてもよいし、一の装置内に一体として備えていてもよい。
[Prediction system 100]
Based on FIG. 1, a prediction system 100 used for predicting the production performance of agricultural products will be described. The prediction system 100 includes a prediction device 10 , a prediction model generation device 20 , and a weather data prediction model generation device 30 . Moreover, the prediction system 100 further includes an input device 40 , a storage device 50 and an output device 60 . The prediction system 100 may include the prediction device 10, the prediction model generation device 20, and the prediction model generation device 30 as independent devices, or may be integrated in one device.

入力装置40は、ユーザによる予測システム100に対する入力操作を受け付ける。入力装置40は、一例として、予測装置10において農作物の生産成績を予測するために用いるデータの入力を受け付ける。また、入力装置40は、予測モデル生成装置20において予測モデルを生成するために用いるデータの入力を受け付ける。さらに、入力装置40は、予測モデル生成装置30において気象データの予測モデルを生成するために用いるデータの入力を受け付ける。入力装置40は、例えば、キーボード、マウス、タッチセンサなどであってもよい。 The input device 40 receives an input operation to the prediction system 100 by the user. As an example, the input device 40 receives input of data used by the prediction device 10 to predict the production results of agricultural products. The input device 40 also receives input of data used for generating a prediction model in the prediction model generation device 20 . Furthermore, the input device 40 accepts input of data used by the prediction model generation device 30 to generate a prediction model of weather data. The input device 40 may be, for example, a keyboard, mouse, touch sensor, or the like.

記憶装置50は、予測システム100にて使用されるプログラム及びデータを記憶する。記憶装置50は、一例として、入力装置40を介して入力された各種データを記憶している。また、記憶装置50は、一例として、予測装置10において農作物の生産成績を予測するために用いる予測モデル、入力情報、及び出力情報を記憶している。さらに、記憶装置50は、一例として、予測モデル生成装置20において、予測モデルの生成に使用する学習データ及び生成した予測モデルを記憶している。また、記憶装置50は、一例として、予測モデル生成装置30において、予測モデルの生成に使用する学習データ及び生成した予測モデルを記憶している。記憶装置50は、各種データを記憶するデータベースをクラウド又はサーバ上に有していてもよい。 Storage device 50 stores programs and data used in prediction system 100 . The storage device 50 stores, for example, various data input via the input device 40 . The storage device 50 also stores, as an example, a prediction model, input information, and output information used by the prediction device 10 to predict the production results of agricultural products. Further, the storage device 50 stores, as an example, the learning data used to generate the prediction model and the generated prediction model in the prediction model generation device 20 . In addition, the storage device 50 stores, as an example, the learning data used to generate the prediction model and the generated prediction model in the prediction model generation device 30 . The storage device 50 may have a database on a cloud or a server that stores various data.

出力装置60は、予測装置10が予測した結果を出力する。また、出力装置60は、予測装置10が予測した結果に基づく栽培管理指導情報を出力してもよい。栽培管理指導情報には、例えば、マルチング実施時期、排水路の整備に関する情報、摘果数や摘果時期、施肥に関する情報、断根の実施に関する情報等が含まれる。 The output device 60 outputs the result predicted by the prediction device 10 . Moreover, the output device 60 may output cultivation management guidance information based on the result predicted by the prediction device 10 . The cultivation management guidance information includes, for example, mulching implementation time, information on the maintenance of drainage channels, information on the number and time of fruit thinning, information on fertilization, information on implementation of pruning, and the like.

出力装置60による出力の態様は特に限定されない。出力装置60は、例えば、当該情報を画像として表示する表示装置、当該情報を印刷する印刷装置、又は、当該情報を音声として出力する警報装置であってもよい。また、出力装置60は、予測装置10が予測した結果や、結果に基づく栽培管理指導情報等を表示する、スマートフォンのようなモバイルデバイスのディスプレイであってもよい。 The mode of output by the output device 60 is not particularly limited. The output device 60 may be, for example, a display device that displays the information as an image, a printing device that prints the information, or an alarm device that outputs the information as sound. Also, the output device 60 may be a display of a mobile device such as a smart phone that displays the results predicted by the prediction device 10, cultivation management guidance information based on the results, and the like.

(予測装置10)
予測装置10は、農作物の生産成績を予測する予測装置である。予測装置10は、過去の年毎の生産成績の実測値データを参照して生成された、生産成績の経年変化を表す予測モデルを用いて、前年までの生産成績の実測値の少なくとも一部から当年の生産成績を予測する。
(Prediction device 10)
The prediction device 10 is a prediction device that predicts the production results of agricultural products. The prediction device 10 uses a prediction model representing changes in production performance over time, which is generated with reference to actual measurement data of production performance for each past year. Predict production performance for the current year.

予測装置10は、制御部11を備えている。制御部11は、予測装置10の各部を統括して制御するものであり、一例として、プロセッサ及びメモリにより実現される。この例において、プロセッサはストレージ(不図示)にアクセスし、ストレージに格納されているプログラム(不図示)をメモリにロードし、当該プログラムに含まれる一連の命令を実行する。これにより、制御部11の各部が構成される。当該各部として、制御部11は、入力データ取得部12及び予測部13を備えている。 The prediction device 10 includes a control section 11 . The control unit 11 comprehensively controls each unit of the prediction device 10, and is realized by, for example, a processor and a memory. In this example, a processor accesses storage (not shown), loads a program (not shown) stored in the storage into memory, and executes a sequence of instructions contained in the program. Each unit of the control unit 11 is thus configured. The control unit 11 includes an input data acquisition unit 12 and a prediction unit 13 as the respective units.

入力データ取得部12は、生産成績を予測する予測モデルの入力情報である入力データを取得する。入力データ取得部12は、入力装置40からの予測の開始指示を表す入力信号に基づき、前年までの生産成績の実測値の少なくとも一部を含む実測値データを記憶装置50から読み出す。また、入力データ取得部12は、入力装置40を介して入力された、実測値データを取得してもよい。入力データ取得部12は、取得した入力データを予測部13へ出力する。 The input data acquisition unit 12 acquires input data, which is input information for a prediction model for predicting production results. The input data acquiring unit 12 reads from the storage device 50 the measured value data including at least a part of the measured values of the production results up to the previous year based on the input signal indicating the prediction start instruction from the input device 40 . Also, the input data acquisition unit 12 may acquire measured value data input via the input device 40 . The input data acquisition unit 12 outputs the acquired input data to the prediction unit 13 .

ここで、前年までの生産成績の実測値とは、生産成績を予測する対象となる年である当年よりも前に測定された生産成績の実測値が意図される。生産成績の実測値は、各年において、少なくとも生育期や収穫期又は出荷時の1回測定されていればよい。実測値データは、少なくとも前年の生産成績の実測値のみを含んでいればよいが、前年までの複数年の生産成績の実測値を含んでいてもよい。 Here, the measured values of the production results up to the previous year are intended to mean the measured values of the production results measured before the current year, which is the target year for predicting the production results. Actual production results may be measured at least once each year during the growing season, harvesting season, or shipment. The measured value data may include at least the measured values of the production results of the previous year only, but may also include the measured values of the production results of several years up to the previous year.

予測部13は、過去の年毎の生産成績の実測値データを参照して生成された、生産成績の経年変化を表す予測モデル24を用いて、前年までの生産成績の実測値の少なくとも一部から当年の生産成績を予測する。予測部13は、記憶装置50から予測モデル24を読み出し、入力データ取得部12から送られた実測値データを予測モデル24に入力して、予測モデル24から出力された当年の生産成績を取得する。予測部13は、取得した当年の生産成績を表すデータを、予測結果として、出力装置60へ出力する。 The prediction unit 13 uses a prediction model 24 representing changes in production results over time, which is generated with reference to actual measurement data of production results for each past year, to obtain at least part of the actual measurement values of production results up to the previous year. Predict the production performance for the current year from The prediction unit 13 reads the prediction model 24 from the storage device 50, inputs the measured value data sent from the input data acquisition unit 12 to the prediction model 24, and acquires the current year's production results output from the prediction model 24. . The prediction unit 13 outputs the acquired data representing the production results for the current year to the output device 60 as a prediction result.

予測部13が生産成績を予測するために用いる予測モデル24は、後述する予測モデル生成装置20により生成された予測モデル24であり得る。予測モデル24は、一例として、実測値データを学習データとして機械学習を行うことによって生成されており、前年までの生産成績の実測値の少なくとも一部が入力情報であり、当年の生産成績が出力情報である予測モデルである。すなわち、予測モデル24を用いることにより、農作物の年毎の生産成績の変動(生産成績のゆらぎ)を経年変化として考慮して、生産成績を予測することができる。 The prediction model 24 used by the prediction unit 13 to predict production results may be a prediction model 24 generated by a prediction model generation device 20 described later. As an example, the prediction model 24 is generated by performing machine learning using measured value data as learning data, at least a part of the measured values of the production results up to the previous year is input information, and the production results of the current year are output. It is a predictive model that is information. That is, by using the predictive model 24, it is possible to predict the production results by taking into consideration the annual fluctuations in the production results of agricultural products (fluctuations in the production results) as changes over time.

一例として、ミカンのような多年生植物は、生産成績の良かった年の翌年は、木の貯蔵養分が減少するなどの樹勢の低下により生産成績が低下する傾向にあることが知られている。そして、多年生植物は、表年(生産成績の良い年)と裏年(生産成績の悪い年)とを交互に繰り返す。このように、農作物には年単位での生産成績の変動が生じるため、前年までの生産成績の実測値から当年の生産成績を予測する予測モデル24を用いることで、年毎のゆらぎを考慮した生産成績の予測が可能である。 As an example, it is known that perennial plants such as mandarin oranges tend to have a lower production performance in the year following a year of good production performance due to a decrease in tree vigor, such as a decrease in the stored nutrients of the tree. Perennial plants alternate between front years (years with good production results) and back years (years with poor production results). In this way, since the production results of agricultural products fluctuate on an annual basis, annual fluctuations are taken into consideration by using the prediction model 24 that predicts the production results of the current year from the measured values of the production results up to the previous year. Prediction of production performance is possible.

また、他の例として、単年生植物では、栽培管理方法に年単位での変動が生じ得、この栽培管理方法の変動が、年毎の生産成績のゆらぎに繋がる。したがって、単年生植物においても、前年までの生産成績の実測値から当年の生産成績を予測する予測モデル24を用いることで、栽培管理方法の変動に起因する年毎のゆらぎを考慮した生産成績の予測が可能である。 As another example, in the case of monoennial plants, the method of cultivation management may vary from year to year, and this variation in the method of cultivation management may lead to fluctuations in production results from year to year. Therefore, even for single-year plants, by using the prediction model 24 that predicts the production results of the current year from the measured values of the production results up to the previous year, it is possible to predict the production results that take into account the annual fluctuations caused by changes in cultivation management methods. Predictable.

ここで、生産成績の予測に利用する実測値データについて、図2を参照して説明する。図2は、予測装置10で利用する前年までの年毎の生産成績の実測値データの一例を示す図である。図2の左側のグラフにおいて、X軸は暦、Y軸は糖度、Z軸は年次を示しており、当該グラフは各年の年間の糖度変化を表している。このような各年の年間の糖度変化のデータから、特定の日における年毎の糖度変化のデータを取り出すことで、図2の右側のグラフのように、特定の日における年毎の糖度変化を表すグラフが得られる。なお、図2の左側のグラフは、各年において複数回糖度を測定した結果から得られたものである。図2は、本発明の概念を説明するために示すものであり、本発明を図2に示す態様に限定するものではない。 Here, the measured value data used for predicting production results will be described with reference to FIG. FIG. 2 is a diagram showing an example of actual measurement value data of production results for each year up to the previous year, which is used by the prediction device 10. As shown in FIG. In the graph on the left side of FIG. 2, the X axis indicates the calendar, the Y axis indicates the sugar content, and the Z axis indicates the year. By extracting data on annual changes in sugar content on specific days from such data on annual changes in sugar content in each year, we can obtain yearly changes in sugar content on specific days, as shown in the graph on the right side of Fig. 2. A graph representing The graph on the left side of FIG. 2 was obtained from the results of measuring the sugar content multiple times in each year. FIG. 2 is shown for explaining the concept of the present invention and does not limit the present invention to the embodiment shown in FIG.

一例として、図2においては、過去15年間の各年の年間の糖度変化を表すグラフ(左側のグラフ)と、10月30日の過去15年間の糖度変化を表すグラフ(右側のグラフ)とを示している。予測装置10は、このように、特定の日における年毎の糖度変化を表すデータを用いることで、年毎のゆらぎを考慮して生産成績を予測することができる。 As an example, in FIG. 2, a graph (left graph) representing annual changes in sugar content for the past 15 years and a graph (right graph) representing changes in sugar content for the past 15 years on October 30 are shown. showing. In this way, the prediction device 10 can predict the production results in consideration of annual fluctuations by using the data representing the annual sugar content change on a specific day.

生産成績の年間の変化を実測する場合、基準となる植物個体(決められた1本の基準木等)において年に複数回測定を繰り返すことになる。つまり、サンプリングのために、基準となる植物個体から葉や実を年間に複数回間引くことになる。サンプリングにより葉や実がなくなったことで、植物個体内の光合成産物の転流が変わったり、サンプリングによる傷が傷害ストレスとなったりして、植物個体の状態が変化し得る。そのため、通常の生育とは異なった状態で生産成績を測定することとなり、実際の生育過程を表していないデータになる可能性がある。このようにして得られたデータに基づき農作物の生産成績を予測すれば、予測精度が低下する虞がある。予測装置10では、年間の生産成績の変化ではなく、特定の日における年毎の生産成績の変化を考慮して生産成績を予測するので、精度よく予測することが可能である。 When actually measuring annual changes in production results, the measurement is repeated multiple times a year using a reference plant individual (one predetermined reference tree or the like). In other words, for sampling, leaves and fruits are thinned out several times a year from a reference plant individual. Loss of leaves and fruits due to sampling can change the translocation of photosynthetic products in individual plants, and injury stress caused by sampling can change the state of individual plants. Therefore, the production results are measured in a state different from normal growth, and the data may not represent the actual growth process. If the production results of agricultural products are predicted based on the data obtained in this way, there is a risk that the prediction accuracy will deteriorate. The prediction device 10 predicts the production result by taking into consideration the yearly change in the production result on a specific day rather than the yearly change in the production result, so it is possible to make a highly accurate prediction.

また、予測部13は、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを参照して生成された、生産成績の経年変化と気象の影響とを表す予測モデル24を用いて、当年の生産成績を予測し得る。予測部13は、記憶装置50から予測モデル24を読み出し、入力データ取得部12から送られた実測値データと当年の気象データの予報値とを予測モデル24に入力して、予測モデル24から出力された当年の生産成績を取得する。気象データの予報値は、一例として、気象庁等の気象情報を提供する機関が公開しているデータベースから取得することができる。予測装置10は、取得した気象データの予報値を記憶装置50に格納し、予測部13において生産成績を予測する際に、記憶装置50から読み出してもよい。 In addition, the prediction unit 13 is generated by referring to measurement data including measured value data of production results for each past year and observation value data of meteorological data for those years, and predicts changes in production results over time and meteorological data. A predictive model 24 representing the impact can be used to predict the current year's production performance. The prediction unit 13 reads out the prediction model 24 from the storage device 50, inputs the actually measured value data sent from the input data acquisition unit 12 and the forecast value of the weather data for the current year into the prediction model 24, and outputs from the prediction model 24. Get the production results for the current year. Forecast values of weather data can be obtained, for example, from a database published by an organization that provides weather information, such as the Meteorological Agency. The prediction device 10 may store the forecast value of the obtained weather data in the storage device 50 and read it from the storage device 50 when the prediction unit 13 predicts the production results.

予測モデル24は、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを学習データとして機械学習を行うことによって生成される。そして、このような予測モデル24は、前年までの生産成績の実測値の少なくとも一部と当年の気象データの予報値とが入力情報であり、当年の生産成績が出力情報である予測モデルである。このような予測モデル24を用いることにより、農作物の年毎の生産成績のゆらぎと共に、生産成績に影響を及ぼす気象データを考慮して、生産成績を予測することができる。農業に密接に関連する気象データをさらに考慮することによって、年毎の生産成績のゆらぎのみを考慮して生産成績を予測する場合よりも、より精度よく予測することができる。 The predictive model 24 is generated by performing machine learning using measured data including actual measured value data of production results for each past year and observed value data of weather data for those years as learning data. Such a prediction model 24 is a prediction model whose input information is at least a part of the measured values of the production results up to the previous year and the forecast values of the current year's weather data, and whose output information is the production results of the current year. . By using such a prediction model 24, it is possible to predict the production results in consideration of the fluctuations in the production results of agricultural products from year to year and the weather data that affects the production results. By additionally considering meteorological data closely related to agriculture, it is possible to make more accurate predictions than when predicting production results by considering only fluctuations in production results from year to year.

経年変化及び気象が生産成績に及ぼす影響について、図3を参照して説明する。図3は、農作物の生産成績の経年変化及び気象が生産成績に及ぼす影響について説明する図である。図3のグラフは、X軸が年次、Y軸が糖度を示しており、当該グラフは特定の日における年毎の糖度変化を表すグラフである。この年毎の糖度変化には、生産成績の経年変化だけでなく、気象が影響を及ぼしていると考えられる。そこで、図3の右側のグラフのように、年毎の糖度変化(実線のデータ)を、経年変化の影響による糖度変化(破線のデータ)と、破線のデータと実線のデータとの差により表される気象の影響による糖度変化(矢印の範囲)とに分解する。 The effects of aging and weather on production performance are described with reference to FIG. FIG. 3 is a diagram for explaining changes over time in the production results of agricultural products and the effects of weather on the production results. In the graph of FIG. 3, the X-axis indicates the year and the Y-axis indicates the sugar content, and the graph is a graph showing changes in the sugar content for each year on a specific day. It is thought that not only changes in production results over time but also weather influences the change in sugar content from year to year. Therefore, as shown in the graph on the right side of Fig. 3, the change in sugar content (solid line data) for each year is expressed by the sugar content change (dashed line data) due to aging and the difference between the broken line data and the solid line data. It decomposes into sugar content changes (range of arrows) due to the influence of weather.

図3の破線で表される経年変化の影響を表す糖度は、上述した年毎の生産成績のゆらぎを考慮して予測される糖度であり、これを基準糖度(基準生産成績)とする。基準糖度は、年毎のゆらぎのみを考慮して予測されるので、図3の実線で表される実際の糖度とは異なり得る。この実際の糖度と年毎のゆらぎを考慮して予測される基準糖度との差は、気象が糖度に及ぼす影響分であると考えられる。 The sugar content indicated by the dashed line in FIG. 3, which represents the effect of aging, is a sugar content that is predicted in consideration of the above-mentioned fluctuations in production results from year to year, and is defined as a standard sugar content (standard production performance). Since the standard sugar content is predicted by considering only annual fluctuations, it may differ from the actual sugar content represented by the solid line in FIG. The difference between the actual sugar content and the reference sugar content predicted by considering annual fluctuations is considered to be the influence of the weather on the sugar content.

予測モデル24は、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データから当年の生産成績を予測する予測モデルであるため、年毎のゆらぎのみならず、気象の影響も考慮した生産成績の予測が可能である。したがって、図3の右側のグラフにおいて、破線により表される基準生産成績に、気象の影響を表す矢印の範囲の補正が加えられて、実線により表される実際の生産成績により近い、精度のよい生産成績の予測が可能である。 Since the prediction model 24 is a prediction model that predicts the current year's production results from the measured data including the measured value data of the production results for each past year and the observed value data of the meteorological data of those years, the annual fluctuation In addition, it is possible to predict production results taking into account the influence of weather. Therefore, in the graph on the right side of FIG. 3, the reference production performance represented by the dashed line is corrected in the range of the arrows representing the influence of weather, resulting in a more accurate production performance closer to the actual production performance represented by the solid line. Prediction of production performance is possible.

また、予測モデル24は、予測対象地点毎に生成した予測モデルであり、当該予測対象地点における前年までの生産成績の実測値の少なくとも一部と当該予測対象地点における当年の気象データの予報値とから、当該予測対象地点における当年の生産成績を予測するモデルであり得る。 In addition, the forecast model 24 is a forecast model generated for each forecast target point, and includes at least a part of the actual measurement value of the production performance up to the previous year at the forecast target point and the forecast value of the current year's weather data at the forecast target point. Therefore, it can be a model that predicts the production performance of the current year at the prediction target location.

生産成績に及ぼす経年変化や気象の影響は、農作物を生育する地域毎に異なり得る。したがって、生産成績を予測する農作物を生育する地域である予測対象地点毎に作成された予測モデル24を用いて生産成績を予測することで、より精度よく生産成績を予測することができる。 The effects of aging and weather on production performance can vary from region to region where crops are grown. Therefore, by predicting the production results using the prediction model 24 created for each prediction target point, which is the region where the crops whose production results are to be predicted, are grown, it is possible to predict the production results with higher accuracy.

予測装置10によれば、前年までの生産成績の実測値の少なくとも一部により、当年の生産成績を予測することができるので、栽培期間の初期であっても、その年の収穫期又は出荷日の生産成績を予測することができる。これにより、栽培期間の初期にその年の生産成績を予測し、予測結果に基づいて、施肥、灌水、摘果等の栽培管理を変更することで、その年の生産成績を向上させることができる。また、予測した生産成績に基づいて、最適な作業時期や収穫時期を決定することもできる。 According to the prediction device 10, it is possible to predict the production results of the current year based on at least part of the measured values of the production results up to the previous year. It is possible to predict the production performance of As a result, the year's production results can be predicted at the beginning of the cultivation period, and cultivation management such as fertilization, irrigation, and fruit thinning can be changed based on the prediction results, thereby improving the year's production results. In addition, it is possible to determine the optimum working time and harvesting time based on the predicted production results.

(予測モデル生成装置20)
予測モデル生成装置20は、農作物の生産成績を予測する予測モデルの生成装置である。予測モデル生成装置20は、過去の年毎の生産成績の実測値データを参照して、生産成績の経年変化を表す予測モデルであって、前年までの生産成績の実測値の少なくとも一部から当年の生産成績を予測する予測モデルを生成する。
(Prediction model generation device 20)
The predictive model generating device 20 is a predictive model generating device that predicts the production results of agricultural products. The predictive model generation device 20 refers to actual measured value data of production results for each past year, and is a predictive model representing changes in production results over time. Generate a predictive model that predicts the production performance of

予測モデル生成装置20は、制御部21を備えている。制御部21は、予測モデル生成装置20の各部を統括して制御するものであり、一例として、プロセッサ及びメモリにより実現される。この例において、プロセッサはストレージ(不図示)にアクセスし、ストレージに格納されているプログラム(不図示)をメモリにロードし、当該プログラムに含まれる一連の命令を実行する。これにより、制御部21の各部が構成される。当該各部として、制御部21は、学習データ取得部22及び予測モデル生成部(生成部)23を備えている。 The predictive model generation device 20 includes a control section 21 . The control unit 21 controls each unit of the prediction model generation device 20 in an integrated manner, and is realized by, for example, a processor and a memory. In this example, a processor accesses storage (not shown), loads a program (not shown) stored in the storage into memory, and executes a sequence of instructions contained in the program. Thereby, each part of the control part 21 is comprised. The control unit 21 includes a learning data acquisition unit 22 and a prediction model generation unit (generation unit) 23 as the respective units.

学習データ取得部22は、予測モデル24を学習するための学習データを取得する。学習データ取得部22は、入力装置40からの学習の開始指示を表す入力信号に基づき、学習データを記憶装置50から読み出す。また、学習データ取得部22は、入力装置40を介して入力された学習データを取得してもよい。学習データ取得部22は、取得した学習データを予測モデル生成部23へ出力する。 The learning data acquisition unit 22 acquires learning data for learning the prediction model 24 . The learning data acquisition unit 22 reads out learning data from the storage device 50 based on an input signal representing an instruction to start learning from the input device 40 . Also, the learning data acquisition unit 22 may acquire learning data input via the input device 40 . The learning data acquisition unit 22 outputs the acquired learning data to the prediction model generation unit 23 .

予測モデル生成部23は、過去の年毎の生産成績の実測値データを参照して、生産成績の経年変化を表す予測モデルであって、前年までの生産成績の実測値の少なくとも一部から当年の生産成績を予測する予測モデル24を生成する。 The prediction model generation unit 23 refers to the measured value data of the production results for each past year to create a prediction model representing changes in the production results over time. generate a predictive model 24 that predicts the production performance of

予測モデル生成部23は、学習データ取得部22から取得した学習データを用いて機械学習を行うことにより、予測モデル24を生成する。ここで学習データは、実測値データである。予測モデル生成部23は、一例としては、ニューラルネットワーク、決定木、ランダムフォレスト、サポートベクトルマシン等の既知の機械学習方法を用いて、予測モデル24を生成する。予測モデル生成装置20により生成された予測モデル24は、前年までの生産成績の実測値の少なくとも一部を入力情報として、当年の生産成績を出力情報とする予測モデルである。 The prediction model generation unit 23 generates the prediction model 24 by performing machine learning using the learning data acquired from the learning data acquisition unit 22 . Here, the learning data are measured value data. The prediction model generation unit 23 generates the prediction model 24 using, for example, known machine learning methods such as neural networks, decision trees, random forests, and support vector machines. The prediction model 24 generated by the prediction model generation device 20 is a prediction model that uses as input information at least a part of the measured values of the production results up to the previous year and outputs the current year's production results.

また、予測モデル生成部23は、一例として、ある年の前年までの年毎の生産成績とその年の生産成績との相関を表す統計モデルを生成してもよい。すなわち、予測モデル24は、線形回帰モデル、スパースモデリング、一般化線形モデル、状態空間モデル、階層ベイズモデル、時系列モデル、クラスタリング等の統計モデルであってもよい。 In addition, the predictive model generation unit 23 may generate, as an example, a statistical model representing the correlation between the production results for each year up to the previous year and the production results for that year. That is, the predictive models 24 may be statistical models such as linear regression models, sparse modeling, generalized linear models, state space models, hierarchical Bayesian models, time series models, clustering, and the like.

予測モデル24は、過去の年毎の生産成績の実測値データを参照して、生産成績の経年変化を表すモデルであるため、農作物の年毎の生産成績の変動を経年変化として考慮した、生産成績の予測に利用できる。予測モデル24は、上述した予測装置10において生産成績を予測するために用いられる予測モデルである。予測モデル生成部23は、生成した予測モデル24を記憶装置50に格納してもよい。 Since the prediction model 24 is a model that represents changes in production results over time by referring to actual measurement data of production results for each year in the past, production It can be used to predict grades. The prediction model 24 is a prediction model used for predicting production results in the prediction device 10 described above. The prediction model generator 23 may store the generated prediction model 24 in the storage device 50 .

また、予測モデル生成部23は、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを参照して、生産成績の経年変化と気象の影響とを表す予測モデル24を生成し得る。このような予測モデル24は、前年までの生産成績の実測値の少なくとも一部と当年の気象データの予報値とから、当年の生産成績を予測する予測モデルである。気象データの観測値データは、一例として、気象庁等の気象情報を提供する機関が公開しているデータベースから取得することができる。予測モデル生成装置20は、取得した気象データの観測値を記憶装置50に格納し、予測モデル生成部23において予測モデル24を生成する際に、記憶装置50から読み出してもよい。 In addition, the prediction model generating unit 23 refers to measurement data including measured value data of production results for each past year and observation value data of weather data for those years, and determines changes in production results over time and effects of weather. can generate a predictive model 24 that represents Such a prediction model 24 is a prediction model that predicts the current year's production performance based on at least part of the actual production performance values up to the previous year and the current year's weather data forecast values. Observed value data of meteorological data can be obtained, for example, from a database published by an organization that provides meteorological information, such as the Meteorological Agency. The predictive model generation device 20 may store the obtained observation values of the weather data in the storage device 50 and read them from the storage device 50 when the prediction model generation unit 23 generates the prediction model 24 .

予測モデル生成部23は、測定データを学習データとして機械学習を行うことによって、前年までの生産成績の実測値の少なくとも一部と当年の気象データの予報値とが入力情報であり、当年の生産成績が出力情報である予測モデルを生成する。また、予測モデル生成部23は、一例として、ある年の前年までの年毎の生産成績及びその年の気象データの観測値と、その年の生産成績との相関を表す統計モデルを予測モデル24として生成してもよい。 By performing machine learning using measured data as learning data, the prediction model generation unit 23 uses as input information at least a part of the actual measurement values of the production results up to the previous year and the forecast values of the weather data for the current year. Generate a predictive model whose output is grades. In addition, the prediction model generation unit 23 generates, as an example, a statistical model representing the correlation between the production results for each year up to the previous year and the observation values of the weather data for that year, and the production results for that year. may be generated as

予測モデル24は、農作物の年毎の生産成績のゆらぎと共に、生産成績に影響を及ぼす気象データを考慮した、生産成績の予測に利用できる。農業に密接に関連する気象データを考慮することによって、年毎の生産成績のゆらぎのみを考慮して生産成績を予測する場合よりも、より精度よく予測することに利用できる。 The prediction model 24 can be used for forecasting production results, taking into consideration weather data that affects production results, as well as annual fluctuations in the production results of crops. By considering meteorological data closely related to agriculture, it can be used for more accurate prediction than when predicting production results by considering only fluctuations in production results from year to year.

予測モデル生成部23は、予測対象地点毎に予測モデル24を生成することが好ましい。予測モデル生成部23は、当該予測対象地点における過去の年毎の生産成績の実測値と当該予測対象地点におけるそれらの年の気象データの観測値とを含む測定データを参照して、予測対象地点毎の予測モデル24を生成する。 The prediction model generator 23 preferably generates the prediction model 24 for each prediction target point. The prediction model generation unit 23 refers to the measurement data including the actual measurement values of the production results of each past year at the prediction target point and the observed values of the weather data of those years at the prediction target point, and the prediction target point Generate a prediction model 24 for each.

生産成績に及ぼす経年変化や気象の影響は、農作物を生育する地域毎に異なり得る。したがって、生産成績を予測する農作物を生育する地域である予測対象地点毎に予測モデル24を生成することで、より精度の高い生産成績の予測に利用できる。 The effects of aging and weather on production performance can vary from region to region where crops are grown. Therefore, by generating the prediction model 24 for each prediction target point, which is an area where the agricultural products whose production performance is to be predicted, is grown, it can be used for more accurate production performance prediction.

予測モデル生成部23は、予測精度を向上させるために、予測結果に基づいて予測モデル24を再生成してもよい。予測モデル24が機械学習により生成された学習モデルである場合、予測モデル24を再学習してもよい。一例として、予測モデル生成部23は、過去のある年の生産成績の実測値と、予測モデル24を用いて予測した、その年の生産成績の予測値との差が最小になるように、予測モデル24を再生成する。 The prediction model generator 23 may regenerate the prediction model 24 based on the prediction result in order to improve the prediction accuracy. If the prediction model 24 is a learning model generated by machine learning, the prediction model 24 may be re-learned. As an example, the prediction model generation unit 23 makes a prediction so that the difference between the measured value of the production performance of a certain past year and the predicted value of the production performance of that year predicted using the prediction model 24 is minimized. Regenerate model 24 .

予測モデル生成部23は、一例として、前年までの各年の生産成績の実測値と、予測モデル24を用いて予測した、それらの年の生産成績の予測値との差を年毎に算出し、全ての年においてその差が最小になるように、予測モデル24を再生成する。また、予測モデル生成部23は、過去のある年の生産成績の実測値と、予測モデル24を用いて予測した、その年の生産成績の予測値との差の二乗和が最小になるように予測モデル24を再生成してもよい。 As an example, the prediction model generation unit 23 calculates, for each year, the difference between the measured production results for each year up to the previous year and the predicted production results for those years predicted using the prediction model 24. , regenerate the forecast model 24 such that the difference is minimized in all years. In addition, the prediction model generation unit 23 generates a value that minimizes the sum of the squares of the differences between the measured production results of a certain past year and the predicted production results of that year predicted using the prediction model 24. Predictive model 24 may be regenerated.

さらに、予測モデル生成部23は、予測モデル24を再生成する場合に、過去のある年の生産成績の実測値と、予測モデル24により予測した、その年の生産成績の予測値との差を、その年の気象データの観測値に基づいて回帰する回帰分析を用いて、予測モデルを再生成してもよい。 Furthermore, when regenerating the prediction model 24, the prediction model generation unit 23 calculates the difference between the measured value of the production performance of a certain past year and the predicted value of the production performance of that year predicted by the prediction model 24. , the forecast model may be regenerated using regression analysis that recurs on observations of weather data for the year.

予測装置10に関して説明したように、図3を参照すれば、年毎の生産成績の変化を、経年変化の影響による生産成績の変化と、気象の影響による生産成績の変化とに分解することができる。予測モデル生成部23は、経年変化の影響による生産成績の変化と気象の影響による生産成績の変化とのバランスを考慮して、年毎の生産成績の変化を、経年変化の影響による生産成績の変化と、気象の影響による生産成績の変化とに分解する。 As described with respect to the prediction device 10, with reference to FIG. 3, it is possible to decompose changes in production results from year to year into changes in production results due to aging and changes in production results due to weather effects. can. The predictive model generating unit 23 considers the balance between changes in production results due to the effects of aging and changes in production results due to the effects of weather, and calculates changes in production results from year to year as changes in production results due to the effects of changes over time. and changes in production performance due to weather effects.

予測モデル生成部23は、一例として、正則化を用いて、経年変化を考慮した生産成績の予測値である基準生産成績を算出する。まず、正則化の強さを様々に変化させて、正則化係数毎に基準生産成績を算出する。そして、生産成績の測定値と正則化係数毎の基準生産成績との差を求め、この差を生産成績に及ぼす気象の影響分としてそれぞれ算出する。そして、気象の影響分をその年の気象データの観測値を用いて回帰する。予測モデル生成部23は、このような処理を繰り返して、生産成績の予測値と実測値との差が最小になる正則化係数、基準生産成績、気象データの回帰係数を決定し、予測モデル24を再生成する。 As an example, the predictive model generation unit 23 uses regularization to calculate a standard production result, which is a predicted value of production results in consideration of aging. First, the standard production performance is calculated for each regularization coefficient by varying the strength of regularization. Then, the difference between the measured value of the production performance and the standard production performance for each regularization coefficient is obtained, and this difference is calculated as the influence of the weather on the production performance. Then, the meteorological influence is regressed using the observed meteorological data for that year. The prediction model generation unit 23 repeats such processing to determine the regularization coefficient that minimizes the difference between the predicted value of the production performance and the measured value, the standard production performance, and the regression coefficient of the weather data. to regenerate.

なお、予測モデル生成部23が予測モデル24を再生成する場合に行う、気象データを用いた生産成績に及ぼす気象の影響分の回帰分析の一例は、後述する予測モデル生成装置30において予測モデル生成部33が行う回帰分析であり得る。 An example of regression analysis of the influence of weather on production performance using weather data, which is performed when the prediction model generation unit 23 regenerates the prediction model 24, is generated by the prediction model generation device 30 described later. It can be a regression analysis performed by the part 33 .

このように、予測モデル24を再生成することによって、より精度よく生産成績を予測することが可能な予測モデル24を生成することができる。 By regenerating the prediction model 24 in this way, it is possible to generate the prediction model 24 that is capable of predicting production results with higher accuracy.

予測モデル生成装置20によれば、前年までの生産成績の実測値の少なくとも一部により、当年の生産成績を予測する予測モデルを生成することができるので、栽培期間の初期にその年の収穫期又は出荷日の生産成績を予測することにも利用することができる。これにより、栽培期間の初期にその年の生産成績を予測し、予測結果に基づいて、施肥、灌水、摘果等の栽培管理を変更することで、その年の生産成績を向上させることに利用することができる。また、予測した生産成績に基づいて、最適な作業時期や収穫時期を決定することにも利用できる。 According to the prediction model generating device 20, it is possible to generate a prediction model for predicting the production results of the current year based on at least part of the measured values of the production results up to the previous year. Alternatively, it can be used to predict the production performance on the shipping date. As a result, it is possible to predict the production results of the year at the beginning of the cultivation period, and to improve the production results of the year by changing cultivation management such as fertilization, irrigation, and fruit thinning based on the prediction results. be able to. It can also be used to determine the optimum working and harvesting times based on predicted production results.

(予測モデル生成装置30)
気象データの予測モデル生成装置30は、農作物の生産成績を予測する予測モデルの生成装置である。予測モデル生成装置30は、予測対象日以前の気象データの観測値データを参照して、生産成績と気象の影響とを表す予測モデルであって、予測対象日の気象データの予報値から生産成績を予測する予測モデル34を生成する。予測モデル生成装置30は、過去のある年の生産成績の実測値と、予測モデル34を用いて予測したその年の生産成績の予測値との差が最小となるように、気象データの変数の組み合わせを選択して、予測モデル34を生成する。
(Prediction model generation device 30)
The prediction model generation device 30 for weather data is a prediction model generation device for predicting the production results of agricultural products. The prediction model generation device 30 is a prediction model representing the production performance and the influence of weather by referring to the observation value data of the weather data before the forecast target date, and is based on the forecast values of the weather data for the forecast target date. generates a prediction model 34 that predicts The predictive model generation device 30 adjusts variables of the weather data so that the difference between the measured value of the production performance of a certain past year and the predicted value of the production performance of that year predicted using the forecast model 34 is minimized. A combination is selected to generate a predictive model 34 .

予測モデル生成装置30は、制御部31を備えている。制御部31は、予測モデル生成装置30の各部を統括して制御するものであり、一例として、プロセッサ及びメモリにより実現される。この例において、プロセッサはストレージ(不図示)にアクセスし、ストレージに格納されているプログラム(不図示)をメモリにロードし、当該プログラムに含まれる一連の命令を実行する。これにより、制御部31の各部が構成される。当該各部として、制御部31は、学習データ取得部32及び予測モデル生成部33を備えている。 The predictive model generation device 30 includes a control section 31 . The control unit 31 controls each unit of the prediction model generation device 30 in an integrated manner, and is realized by, for example, a processor and a memory. In this example, a processor accesses storage (not shown), loads a program (not shown) stored in the storage into memory, and executes a sequence of instructions contained in the program. Each unit of the control unit 31 is thus configured. The control unit 31 includes a learning data acquisition unit 32 and a prediction model generation unit 33 as the respective units.

学習データ取得部32は、予測モデル34を学習するための学習データを取得する。学習データ取得部32は、入力装置40からの学習の開始指示を表す入力信号に基づき、学習データを記憶装置50から読み出す。また、学習データ取得部32は、入力装置40を介して入力された学習データを取得してもよい。学習データ取得部32は、取得した学習データを予測モデル生成部33へ出力する。 The learning data acquisition unit 32 acquires learning data for learning the prediction model 34 . The learning data acquisition unit 32 reads out learning data from the storage device 50 based on an input signal representing a learning start instruction from the input device 40 . Also, the learning data acquisition unit 32 may acquire learning data input via the input device 40 . The learning data acquisition unit 32 outputs the acquired learning data to the prediction model generation unit 33 .

予測モデル生成部33は、予測対象日以前の気象データの観測値を参照して、生産成績と気象の影響とを表す予測モデルであって、予測対象日の気象データの予報値から生産成績を予測する予測モデル34を生成する。 The prediction model generation unit 33 is a prediction model representing the production results and the effects of weather by referring to the observed values of the weather data before the prediction target day, and calculates the production results from the forecast values of the weather data for the prediction target day. Generate a prediction model 34 to predict.

予測モデル生成部33は、日毎の気象データの観測値とそれらの日の生産成績の実測値とを含む学習データを参照して、生産成績と気象の影響とを表す予測モデル34を生成する。予測モデル生成部33は、学習データ取得部22から取得した学習データを用いて機械学習を行うことにより、予測モデル34を生成する。予測モデル生成部33は、一例としては、ニューラルネットワーク、決定木、ランダムフォレスト、サポートベクトルマシン等の既知の機械学習方法を用いて、予測モデル34を生成する。 The prediction model generating unit 33 generates a prediction model 34 that represents the production results and the effects of weather by referring to the learning data including the daily observed weather data and the actual measurement values of the production results for those days. The prediction model generation unit 33 generates the prediction model 34 by performing machine learning using the learning data acquired from the learning data acquisition unit 22 . The prediction model generation unit 33 generates the prediction model 34 using, for example, known machine learning methods such as neural networks, decision trees, random forests, and support vector machines.

また、予測モデル生成部33は、一例として、日毎の気象データとそれらの日の生産成績との相関を表す統計モデルを生成してもよい。すなわち、予測モデル34は、線形回帰モデル、スパースモデリング、一般化線形モデル、状態空間モデル、階層ベイズモデル、時系列モデル、クラスタリング等の統計モデルであってもよい。 In addition, the predictive model generation unit 33 may generate, as an example, a statistical model representing the correlation between the daily weather data and the production results for those days. That is, the predictive model 34 may be a statistical model such as a linear regression model, sparse modeling, generalized linear model, state-space model, hierarchical Bayesian model, time series model, clustering, or the like.

予測モデル生成装置30は、過去のある年の生産成績の実測値と、予測モデル34を用いて予測したその年の生産成績の予測値との差が最小となるように、気象データの変数の組み合わせを選択して、予測モデル34を生成する。農作物の生産成績に影響することが予想される気象データには、降水量、平均気温、最高気温、最低気温、全天日射量、日照時間等の多くの変数が含まれる。気象データの変数には、相互依存する変数が存在し、一例として、雨が降る日は気温が低め(降水量と気温との関連性)のように、相互依存性の高い変数の組み合わせが存在する。したがって、生産成績の予測に最適な気象データの変数の組み合わせを選択することで、精度よく生産成績を予測可能な予測モデルを生成することができる。 The predictive model generation device 30 adjusts variables of the weather data so that the difference between the measured value of the production performance of a certain past year and the predicted value of the production performance of that year predicted using the forecast model 34 is minimized. A combination is selected to generate a predictive model 34 . Meteorological data expected to affect crop production performance includes many variables such as precipitation, average temperature, maximum temperature, minimum temperature, global solar radiation, and sunshine hours. Variables in meteorological data include interdependent variables. For example, there are combinations of highly interdependent variables, such as lower temperatures on rainy days (relationship between precipitation and temperature). do. Therefore, by selecting the optimal combination of weather data variables for predicting production results, it is possible to generate a prediction model that can accurately predict production results.

気象データの変数の組み合わせが生産成績に及ぼす影響について、図4を参照して説明する。図4は、気象データの変数が生産成績に及ぼす影響を説明する図である。図4において、上段左側のグラフは降水量の糖度への影響を表しており、上段中央のグラフは平均気温の糖度への影響を表しており、上段右側のグラフは日照時間の糖度への影響を表している。これらのグラフにおいて、X軸は日付、Y軸は糖度の増減値、Z軸は測定対象地点を示している。すなわち、これらのグラフは、測定対象地点毎に、気象データの変数が生産成績に及ぼす影響を示している。 The effects of combinations of weather data variables on production performance will now be described with reference to FIG. FIG. 4 is a diagram for explaining the effects of weather data variables on production results. In Fig. 4, the upper left graph shows the effect of precipitation on sugar content, the upper middle graph shows the effect of average temperature on sugar content, and the upper right graph shows the effect of sunshine hours on sugar content. represents. In these graphs, the X axis indicates dates, the Y axis indicates increases and decreases in sugar content, and the Z axis indicates points to be measured. That is, these graphs show the effects of weather data variables on production performance for each measurement point.

図4において、下段のグラフは、上段の3つのグラフにおいて、特定の測定対象地点のグラフのみを取出して重ねたグラフである。すなわち、このグラフは、気象データのそれぞれの変数が糖度に及ぼす影響をまとめて示すものである。なお、下段のグラフにおいては、実線は降水量の影響、点線は平均気温の影響、○印プロットは日照時間の影響を示している。 In FIG. 4, the lower graph is a graph obtained by superimposing only the graphs of a specific measurement target point in the upper three graphs. That is, this graph collectively shows the influence of each variable of the weather data on the sugar content. In the lower graph, the solid line indicates the effect of precipitation, the dotted line indicates the effect of average temperature, and the circle plot indicates the effect of sunshine hours.

図4の上段のグラフに示すように、場所毎及び日毎に生産成績が変数の影響を受ける程度が異なっている。図4の下段のグラフに示すように、各変数のグラフを重ねることで、同じ場所及び同じ日であっても、影響を受ける変数の種類やその程度が異なっていることが明白である。 As shown in the upper graph of FIG. 4, the degree to which the production results are affected by variables varies from place to place and day to day. As shown in the lower graph of FIG. 4, by superimposing the graphs of each variable, it is clear that the types and degrees of affected variables are different even at the same place and on the same day.

このように、気象データの変数が生産成績に及ぼす影響の程度に基づいて生産成績の予測に最適な変数の組み合わせを選択すれば、より精度よく生産成績を予測可能な予測モデルを生成することができる。また、最適な変数の組み合わせを、日毎又は場所毎に選択することで、より精度よく生産成績を予測可能な予測モデルを生成する事ができる。 In this way, by selecting the optimal combination of variables for predicting production results based on the degree of influence of weather data variables on production results, it is possible to generate a prediction model that can predict production results with higher accuracy. can. In addition, by selecting the optimum combination of variables for each day or each place, it is possible to generate a prediction model that can predict production results with higher accuracy.

予測モデル生成部33は、予測対象日以前の気象データの各変数について、それぞれの積算値又は平均値に基づく回帰分析により、気象データの変数の組み合わせを選択する。予測モデル生成部33は、各変数について、それぞれを所定期間積算した積算値、又は、それぞれの所定期間内の平均値を算出し、算出した積算値又は平均値に基づく回帰分析を行って、生産成績の実測値と予測値との差が最小となる変数の組み合わせを選択する。予測モデル生成部33は、一例として、相互依存性の高い変数を解析可能な機械学習による回帰分析を行うことで、最適な変数の組み合わせを選択する。 The prediction model generation unit 33 selects a combination of weather data variables by regression analysis based on each integrated value or average value for each variable of the weather data before the prediction target date. The prediction model generation unit 33 calculates an integrated value obtained by integrating each variable for a predetermined period or an average value within a predetermined period, performs regression analysis based on the calculated integrated value or average value, and performs production Select the combination of variables that minimizes the difference between the actual and predicted performance. As an example, the prediction model generation unit 33 selects an optimum combination of variables by performing regression analysis by machine learning that can analyze highly interdependent variables.

予測モデル生成部33は、各変数についての積算値又は平均値を算出する起算日及び終算日、並びに算出日数に基づく回帰分析により、気象データの変数の組み合わせを選択する。生産成績に及ぼす気象の影響は、年単位だけでなく、開花期や着果期のような短い期間単位で変化し得る。したがって、各変数について積算を開始する起算日及び積算を終了する終算日、又は、各変数について平均値の算出を開始する起算日及び算出を終了する終算日に基づく回帰分析を行うことで、期間毎に最適な変数の組み合わせを選択することができる。同様に、各変数について積算する日数又は平均する日数(算出日数)に基づく回帰分析を行うことで、より短い期間単位で最適な変数の組み合わせを選択することができる。 The prediction model generation unit 33 selects a combination of weather data variables by regression analysis based on the start date and end date for calculating the integrated value or average value for each variable, and the number of calculation days. The effects of weather on production performance can vary not only over years, but over short periods such as flowering and fruiting. Therefore, by performing a regression analysis based on the start date and end date of calculation for each variable, or the start date and end date of calculation for each variable, , the optimal combination of variables can be selected for each period. Similarly, by performing regression analysis based on the number of days to be integrated or the number of days to average for each variable (calculation number of days), it is possible to select an optimal combination of variables for a shorter period.

予測モデル生成部33は、気象データの各変数について、起算日から終算日までの各日における算出日数の積算値又は平均値に基づく回帰分析により、前記気象データの変数の組み合わせを選択する。このように、起算日から終算日までの算出期間内の各日において算出日数分の各変数の積算値又は平均値を算出することで、算出日数分のウインドウを一日ずつずらし、算出期間内の毎日の各変数が生産成績に及ぼす影響を表すベクトルを作成することができる。そして、作成した毎日のベクトルを用いて回帰分析を行い、生産成績の実測値と予測値との差が最小となる変数の組み合わせを選択することで、生産成績の予測に最適な変数の組み合わせを選択することができる。また、毎日のベクトルを用いた回帰分析を行うことで、連日続く猛暑の影響のように、日毎の変数の連続性が生産成績に及ぼす影響についても考慮して、変数の組み合わせを選択することができる。 The prediction model generation unit 33 selects a combination of the weather data variables by regression analysis based on the integrated value or average value of the calculated days for each day from the start date to the end date for each variable of the weather data. In this way, by calculating the integrated value or average value of each variable for the number of calculation days on each day in the calculation period from the start date to the end date, the window for the number of calculation days is shifted by one day, and the calculation period A vector can be created representing the impact of each daily variable in the production performance. Then, by performing regression analysis using the created daily vector and selecting the combination of variables that minimizes the difference between the actual measurement value and the predicted value of production performance, the optimal combination of variables for forecasting production performance is determined. can be selected. In addition, by performing regression analysis using daily vectors, it is possible to select a combination of variables by considering the effects of the continuity of daily variables on production performance, such as the effects of intense heat that continues for consecutive days. can.

予測モデル生成部33は、過去の年毎の前記生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを参照して生成した、生産成績の経年変化と気象の影響とを表す予測モデルであって、前年までの生産成績の実測値の少なくとも一部と予測対象日の気象データの予報値とから、生産成績を予測する予測モデル34を生成してもよい。また、このような予測モデル34は、測定データを学習データとして機械学習を行うことによって生成し、これにより、前年までの前記生産成績の実測値の少なくとも一部と前記予測対象日の気象データの予報値とが入力情報であり、当年の前記生産成績が出力情報である予測モデル34を生成してもよい。すなわち、予測モデル生成部33は、上述した予測モデル生成装置20の予測モデル生成部23において、年毎の生産成績の変化を経年変化の影響分と気象の影響分とに分解したときの、気象の影響分を回帰して、生産成績の経年変化と気象の影響を表す予測モデル34を生成し得る。 The prediction model generation unit 33 generates changes in production performance over time and weather forecasts by referring to measurement data including measured value data of production performance for each past year and observation value data of weather data for those years. A prediction model 34 for predicting production results may be generated from at least a part of actual measurement values of production results up to the previous year and forecast values of weather data for the day to be predicted. In addition, such a prediction model 34 is generated by performing machine learning using measurement data as learning data, and as a result, at least a part of the actual measurement value of the production performance up to the previous year and the weather data of the prediction target day. A prediction model 34 may be generated in which forecast values are input information and the production results for the current year are output information. That is, the predictive model generation unit 33 determines the weather conditions when the yearly changes in production performance are decomposed into the effects of aging and the effects of weather in the prediction model generation unit 23 of the prediction model generation device 20 described above. can be regressed to generate a predictive model 34 that represents changes in production performance over time and the effects of weather.

また、予測モデル生成部33は、予測対象地点毎に予測モデル34を生成することが好ましい。予測モデル生成部33は、当該予測対象地点における過去の年毎の生産成績の実測値と当該予測対象地点におけるそれらの年の気象データの観測値とを含む測定データを参照して、予測対象地点毎の予測モデル24を生成する。生産成績を予測する農作物を生育する地域である予測対象地点毎に予測モデル34を生成することで、より精度の高い生産成績の予測に利用できる。 Moreover, it is preferable that the prediction model generation unit 33 generates the prediction model 34 for each prediction target point. The prediction model generation unit 33 refers to the measurement data including the actual measured values of the production results for each past year at the prediction target point and the observed values of the weather data for those years at the prediction target point, and the prediction target point Generate a prediction model 24 for each. By generating the prediction model 34 for each prediction target point, which is an area where the crops whose production performance is to be predicted, is grown, it can be used for more accurate prediction of production performance.

予測モデルの生成装置30によれば、予測対象日の気象データの予報値から、当年の生産成績を予測することができるので、栽培期間の初期にその年の収穫期又は出荷日の生産成績を予測することにも利用することができる。これにより、栽培期間の初期にその年の生産成績を予測し、予測結果に基づいて、施肥、灌水、摘果等の栽培管理を変更することで、その年の生産成績を向上させることに利用することができる。また、予測した生産成績に基づいて、最適な作業時期や収穫時期を決定することにも利用できる。 According to the prediction model generation device 30, it is possible to predict the production results of the current year from the forecast values of the weather data for the target day of prediction. It can also be used for prediction. As a result, it is possible to predict the production results of the year at the beginning of the cultivation period, and to improve the production results of the year by changing cultivation management such as fertilization, irrigation, and fruit thinning based on the prediction results. be able to. It can also be used to determine the optimum working and harvesting times based on predicted production results.

(予測システム100における生産成績の予測)
予測システム100における予測モデルを用いた生産成績の予測処理の流れの概念について、図5を参照して説明する。図5は、本発明の一態様に係る予測システムが実行する、予測モデルを用いた農作物の生産成績の予測の概念を説明する図である。なお、図5に示す気象データの変数は一例であり、気象データの変数はこれに限定されない。
(Prediction of production results in prediction system 100)
The concept of the flow of production result prediction processing using the prediction model in the prediction system 100 will be described with reference to FIG. FIG. 5 is a diagram illustrating the concept of predicting the production performance of agricultural products using a prediction model, which is executed by the prediction system according to one aspect of the present invention. Note that the weather data variables shown in FIG. 5 are merely an example, and the weather data variables are not limited to these.

図5に示すように、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを用いて学習済みの予測モデルを用いる。当該予測モデルに、前年までの(過去の)生産成績の実測値と当年の気象データの予報値とを入力する。予測モデルは、生産成績の予測値を、当年の基準生産成績と当年の気象の影響分とに分割する。予測モデルは、気象データの変数の組み合わせにより、当年の気象の影響分を回帰し、最適な変数の組み合わせを選択する。予測モデルは、選択された変数の組み合わせに基づいて当年の気象の影響分を回帰する。予測モデルは、当年の基準生産成績と当年の気象の影響分とから、当年の生産成績の予測値を出力する。 As shown in FIG. 5, a prediction model that has been trained using measured data including actual measured value data of production results for each past year and observed value data of weather data for those years is used. Into the prediction model, actual measurement values of (past) production results up to the previous year and forecast values of weather data for the current year are input. The forecast model divides the predicted value of production performance into the current year's baseline production performance and the current year's weather impact. The prediction model regresses the influence of the current year's weather on the combination of weather data variables and selects the optimum combination of variables. The forecast model regresses the current year's weather impact on the selected combination of variables. The prediction model outputs a predicted value of the current year's production performance based on the current year's standard production performance and the influence of the current year's weather.

(予測処理の流れ)
予測装置10による予測処理(予測方法)の流れについて、図6を参照して説明する。図6は、本発明の一態様に係る予測装置が実行する予測処理の一例を示すフローチャートである。図6に示すように、まず、入力データ取得部12は、前年までの(過去の)生産成績の実測値を取得する(ステップS1)。次に、入力データ取得部12は、当年の気象データの予報値を取得する(ステップS2)。そして、予測部13は、取得した前年までの生産成績の実測値と当年の気象データの予報値とを予測モデル24に入力し、出力された当年の生産成績の予測値を取得する(ステップS3、予測する工程)。予測部は、取得した予測値を予測結果として出力装置60に出力し(ステップS4)、処理を終了する。
(Flow of prediction processing)
A flow of prediction processing (prediction method) by the prediction device 10 will be described with reference to FIG. FIG. 6 is a flowchart showing an example of prediction processing executed by the prediction device according to one aspect of the present invention. As shown in FIG. 6, first, the input data acquiring unit 12 acquires the measured values of the (past) production results up to the previous year (step S1). Next, the input data acquisition unit 12 acquires forecast values of weather data for the current year (step S2). Then, the prediction unit 13 inputs the acquired measured values of the production results up to the previous year and the forecast values of the weather data for the current year to the prediction model 24, and acquires the output predicted values of the production results for the current year (step S3 , the process of predicting). The prediction unit outputs the obtained prediction value to the output device 60 as a prediction result (step S4), and ends the process.

(予測モデル生成処理の流れ)
予測モデル生成装置20による予測モデル24の生成処理(予測モデルの生成方法)の流れについて、図7を参照して説明する。図7は、本発明の一態様に係る予測モデルの生成装置が実行する予測モデル生成処理の一例を示すフローチャートである。図7に示すように、まず、学習データ取得部22は、過去の年毎の生産成績の実測値を取得する(ステップS11)。次に、学習データ取得部22は、それらの年の気象データの観測値を取得する(ステップS12)。そして、学習データ取得部22は、過去の年毎の生産成績に実測値とそれらの年の気象データの観測値とを対応付け、学習データ(測定データ)を生成する(ステップS13)。予測モデル生成部23は、生成した学習データを用いて、前年までの(過去の)生産成績の実測値と当年の気象データの予報値とが入力情報であり、当年の生産成績が出力情報である予測モデル24を生成する(ステップS14、生成する工程)。予測モデル生成部23は、生成した予測モデル24を記憶装置50へ格納する。
(Flow of prediction model generation processing)
The flow of the prediction model generation process (prediction model generation method) by the prediction model generation device 20 will be described with reference to FIG. 7 . FIG. 7 is a flowchart showing an example of a prediction model generation process executed by the prediction model generation device according to one aspect of the present invention. As shown in FIG. 7, first, the learning data acquisition unit 22 acquires actual measurement values of production results for each past year (step S11). Next, the learning data acquisition unit 22 acquires observation values of weather data for those years (step S12). Then, the learning data acquisition unit 22 generates learning data (measurement data) by associating actual measurement values with past yearly production results and observation values of weather data for those years (step S13). Using the generated learning data, the prediction model generation unit 23 uses the actual measurement values of production results up to the previous year (past) and forecast values of weather data for the current year as input information, and the production results for the current year as output information. A certain prediction model 24 is generated (step S14, step of generating). The prediction model generator 23 stores the generated prediction model 24 in the storage device 50 .

(気象データの予測モデル生成処理の流れ)
予測モデル生成装置30による予測モデル34の生成処理(予測モデルの生成方法)の流れについて、図8を参照して説明する。図8は、本発明の一態様に係る気象データの予測モデル生成装置が実行する予測モデルの生成処理の一例を示すフローチャートである。図8に示すように、まず、予測モデル生成部33は、気象データの各変数の積算値又は平均値を算出する起算日及び終算日、並びに、積算する又は平均する日数(算出日数)を設定する(ステップS21)。次に、予測モデル生成部33は、各変数について、起算日から終算日までの各日における算出日数の積算値又は平均値を算出する(ステップS22)。そして、予測モデル生成部33は、算出した値を用いて予測モデル34を回帰し、予測値と実測値との誤差を求める(ステップS23)。予測モデル生成部33は、誤差が最小となる変数の組み合わせ、積算値又は平均値の起算日及び終算日、並びに、算出日数を求める(ステップS24)。予測モデル生成部33は、得られた変数の組み合わせ、積算値又は平均値の起算日及び終算日、並びに、算出日数に基づいて、予測モデル34を生成する(ステップS25、生成する工程)。予測モデル生成部33は、生成した予測モデル34を記憶装置50へ格納する(ステップS26)。
(Flow of weather data prediction model generation processing)
The flow of the prediction model 34 generation process (prediction model generation method) by the prediction model generation device 30 will be described with reference to FIG. 8 . FIG. 8 is a flowchart showing an example of a prediction model generation process executed by the weather data prediction model generation device according to one aspect of the present invention. As shown in FIG. 8 , first, the prediction model generation unit 33 sets the starting date and ending date for calculating the integrated value or average value of each variable of the weather data, and the number of days to be integrated or averaged (calculated number of days). Set (step S21). Next, the predictive model generation unit 33 calculates the integrated value or average value of the number of calculated days in each day from the start date to the end date for each variable (step S22). Then, the prediction model generator 33 regresses the prediction model 34 using the calculated values, and obtains the error between the predicted value and the measured value (step S23). The prediction model generation unit 33 obtains the combination of variables that minimizes the error, the start date and end date of the integrated value or the average value, and the number of calculation days (step S24). The prediction model generation unit 33 generates the prediction model 34 based on the obtained combination of variables, the start date and end date of the integrated value or the average value, and the number of calculation days (step S25, step of generating). The prediction model generator 33 stores the generated prediction model 34 in the storage device 50 (step S26).

〔ソフトウェアによる実現例〕
予測装置10、予測モデル生成装置20、及び予測モデル生成装置30(以下、「装置」と呼ぶ)の機能は、これらの装置としてコンピュータを機能させるためのプログラムであって、これらの装置の各制御ブロック(特に制御部11、制御部21、及び制御部31に含まれる各部)としてコンピュータを機能させるためのプログラムにより実現することができる。
[Example of realization by software]
The functions of the prediction device 10, the prediction model generation device 20, and the prediction model generation device 30 (hereinafter referred to as "devices") are programs for causing a computer to function as these devices, and each control of these devices It can be realized by a program for causing a computer to function as a block (especially each unit included in the control unit 11, the control unit 21, and the control unit 31).

この場合、これらの装置は、上記プログラムを実行するためのハードウェアとして、少なくとも1つの制御装置(例えばプロセッサ)と少なくとも1つの記憶装置(例えばメモリ)を有するコンピュータを備えている。この制御装置と記憶装置により上記プログラムを実行することにより、上記実施形態で説明した各機能が実現される。 In this case, these devices comprise a computer having at least one control device (eg processor) and at least one storage device (eg memory) as hardware for executing the program. Each function described in the above embodiment is realized by executing the above program using the control device and the storage device.

上記プログラムは、一時的ではなく、コンピュータ読み取り可能な、1または複数の記録媒体に記録されていてもよい。この記録媒体は、これらの装置が備えていてもよいし、備えていなくてもよい。後者の場合、上記プログラムは、有線または無線の任意の伝送媒体を介してこれらの装置に供給されてもよい。 The program may be recorded on one or more computer-readable recording media, not temporary. This recording medium may or may not be included in these devices. In the latter case, the program may be supplied to these devices via any transmission medium, wired or wireless.

また、上記各制御ブロックの機能の一部または全部は、論理回路により実現することも可能である。例えば、上記各制御ブロックとして機能する論理回路が形成された集積回路も本発明の範疇に含まれる。この他にも、例えば量子コンピュータにより上記各制御ブロックの機能を実現することも可能である。 Also, part or all of the functions of the above control blocks can be realized by logic circuits. For example, integrated circuits in which logic circuits functioning as the control blocks described above are formed are also included in the scope of the present invention. In addition, it is also possible to implement the functions of the control blocks described above by, for example, a quantum computer.

また、上記実施形態で説明した各処理は、AI(Artificial Intelligence:人工知能)に実行させてもよい。この場合、AIは上記制御装置で動作するものであってもよいし、他の装置(例えばエッジコンピュータまたはクラウドサーバ等)で動作するものであってもよい。 Further, each process described in the above embodiment may be executed by AI (Artificial Intelligence). In this case, the AI may operate on the control device, or may operate on another device (for example, an edge computer or a cloud server).

〔実施例〕
ミカンにおいて、複数の予測対象地点における過去の糖度及び酸度データを用いて生成した予測モデルを用いて糖度及び酸度を予測し、その精度を検証した。糖度予測及び酸度予測の精度は、令和2年度の実測値と予測モデルを用いた予測値とを比較して検証した。平成30年度及び令和1年度の気象観測値、並びに、令和2年度の3月22日時点の気象予報値を用いた。
〔Example〕
For mandarin oranges, we predicted sugar content and acidity using a prediction model generated using past sugar content and acidity data at multiple prediction target points, and verified the accuracy. The accuracy of sugar content prediction and acidity prediction was verified by comparing the measured values in 2020 and the predicted values using the prediction model. The weather observation values in 2018 and 2019, and the weather forecast values as of March 22 in 2020 were used.

2か年の全予測対象地点及び全予測対象日において、予測値と実測値とのRMSE(二乗平均平方根誤差)を求めた。糖度予測において、気象観測値を使用した場合、RMSEは0.6561であった(図9)。糖度予測において、気象予報値を使用した場合、RMSEは1.5267であった(図10)。次に、酸度予測において、気象観測値を使用の場合、RMSEは0.2517であった(図11)。酸度予測において、気象予報値を使用した場合、RMSEは0.2636であった(図12)。 The RMSE (root mean square error) between the predicted values and the measured values was obtained at all the prediction target points and all the prediction target days for two years. When meteorological observations were used in sugar content prediction, the RMSE was 0.6561 (Fig. 9). The RMSE was 1.5267 when the weather forecast value was used in the sugar content prediction (Fig. 10). Next, in the acidity prediction, the RMSE was 0.2517 when meteorological observations were used (Fig. 11). In the acidity prediction, the RMSE was 0.2636 when the weather forecast values were used (Fig. 12).

また、令和2年3月22日時点の気象予報値に基づき、令和2年度の糖度及び酸度を予測し、令和2年11月9日時点までの測定値を用いて精度を検証した。令和2年11月9日時点までの精度の検証では、気象予報値を使用した。糖度予測におけるRMSEは0.7355であった(図13)。酸度予測におけるRMSEは0.1615(図14)であった。 In addition, based on the weather forecast values as of March 22, 2020, the sugar content and acidity in 2020 were predicted, and the accuracy was verified using the measured values up to November 9, 2020. . Weather forecast values were used to verify accuracy up to November 9, 2020. The RMSE in sugar content prediction was 0.7355 (Fig. 13). The RMSE for acidity prediction was 0.1615 (Figure 14).

このように、予測モデルを用いることにより、3月時点で翌年度の生産成績の傾向を把握するために十分な精度で予測できることが示された。 Thus, it was shown that by using the prediction model, it is possible to make predictions with sufficient accuracy as of March to grasp trends in production results for the next fiscal year.

図15は、ある予測対象地点における令和2年度高糖度系統選果予測における気象3変数の時期別の影響を示す。図15において、実線は降水量の影響、点線は平均気温の影響、○印プロットは日照時間の影響を示している。なお、令和2年3月22日時点の気象予報値を使用しているため、令和2年度の実際の気象の影響とは異なる。平年値となる4月下旬以降が平年並みの気象であったとき、降水量、平均気温、及び日照時間の中では、降水量の糖度への影響が最も大きい。図15に示すように、令和2年6月末頃までの降水量は、選果時糖度を日毎で最大0.01上方修正する効果が見られた。一方で、令和2年7月及び8月の降水量は、選果時糖度を日毎で0.01から最大0.025程度下方修正する効果が見られた。この時期は、マルチシート被覆を行っているが、7月及び8月の降水量の影響はゼロではないことが示唆された。 FIG. 15 shows the seasonal effects of the three weather variables in 2020 high sugar content system fruit sorting prediction at a certain prediction target point. In FIG. 15, the solid line indicates the effect of precipitation, the dotted line indicates the effect of average temperature, and the circle plot indicates the effect of sunshine hours. In addition, since the weather forecast value as of March 22, 2020 is used, it differs from the actual weather effects in 2020. When the weather is similar to that of a normal year after late April, which is the normal value, rainfall has the greatest effect on sugar content among rainfall, average temperature, and hours of sunshine. As shown in FIG. 15, the amount of rainfall until the end of June 2020 had the effect of adjusting the sugar content at the time of fruit sorting upward by up to 0.01 on a daily basis. On the other hand, the rainfall in July and August 2020 had the effect of lowering the sugar content at the time of fruit sorting from 0.01 to a maximum of 0.025 on a daily basis. Multi-sheet covering was applied at this time, but it was suggested that the influence of precipitation in July and August was not zero.

〔他の実施形態〕
(生産成績の平均値を用いる態様)
予測装置10は、過去の年毎の生産成績の実測値データを参照して生成された、生産成績の平均値を表す予測モデルを用いて、前年までの生産成績の実測値の少なくとも一部から当年の生産成績を予測してもよい。一例として、過去の年毎の生産成績の蓄積年数が2年から5年程度と少ない場合には、生産成績の経年変化を表す予測モデルではなく、生産成績の平均値を表す予測モデルを用いてもよい。ここで、生産成績の平均値は、所定の年数における各年の生産成績を平均した値であり得る。
また、予測モデル生成装置20が生成する予測モデル24及び予測モデル生成装置30が生成する予測モデル34は、過去の年毎の生産成績の実測値データとそれらの年の気象データの観測値データとを含む測定データを参照して生成された、生産成績の平均値と気象の影響とを表す予測モデルであってもよい。このような予測モデル24及び予測モデル34は、農作物の年毎の変動を平均値として考慮して、生産成績を予測することができる。
[Other embodiments]
(Aspect using the average value of production results)
The prediction device 10 uses a prediction model representing the average value of the production results generated with reference to the actual measurement data of the production results for each past year, based on at least a part of the actual measurement values of the production results up to the previous year. The production performance for the current year may be predicted. As an example, if the number of years of accumulated production results for each past year is as small as 2 to 5 years, a prediction model that expresses the average value of production results is used instead of a prediction model that expresses changes in production results over time. good too. Here, the average value of production results may be a value obtained by averaging production results for each year in a predetermined number of years.
In addition, the prediction model 24 generated by the prediction model generation device 20 and the prediction model 34 generated by the prediction model generation device 30 are based on actual measured value data of production results for each past year and observed value data of weather data for those years. It may be a prediction model representing the average value of production performance and the influence of weather, generated with reference to measurement data including: Such prediction model 24 and prediction model 34 are capable of predicting production results by taking into account annual variations in crops as average values.

すなわち、上述した各実施形態において、生産成績の経年変化に関する説明の全てを、生産成績の平均値に読み替えた実施形態についても、本発明の技術的範囲に含まれる。 That is, in each of the above-described embodiments, the technical scope of the present invention also includes an embodiment in which all the explanations regarding the aging of the production results are replaced with average values of the production results.

(土壌成分の測定値データを測定データとして用いる態様)
予測装置10は、過去の年毎の生産成績の実測値データと、それらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値データの少なくとも一方とを含む測定データを参照して生成された、生産成績の経年変化又は平均値と気象の影響とを表す予測モデルを用いて、当年の生産成績を予測してもよい。ここで、土壌成分の測定値データは、一例として、農作物を生育させる圃場に設けられた、土壌水分計のような土壌センサにより測定された土壌成分値(土壌水分値、土壌pH値等)であり得る。
(Aspect using measured value data of soil components as measured data)
The prediction device 10 refers to measured data including actual measured value data of production results for each past year and at least one of observed value data of meteorological data of those years and measured value data of soil components of those years. A forecasting model of production performance over time or averages and the effects of weather may be used to forecast production performance for the current year. Here, the measured value data of soil components are, for example, soil component values (soil moisture value, soil pH value, etc.) measured by a soil sensor such as a soil moisture meter provided in a field where crops are grown. could be.

予測モデル生成装置20が生成する予測モデル24は、過去の年毎の生産成績の実測値データと、それらの年の気象データの観測値データ及び土壌成分の測定データの少なくとも一方とを含む測定データを参照して生成された、生産成績の経年変化又は平均値と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルであってもよい。このような予測モデル24は、測定データを学習データとして機械学習を行うことによって生成されており、前年までの前記生産成績の実測値の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とが入力情報であり、当年の前記生産成績が出力値である。予測モデル24によれば、土壌成分を考慮して、生産成績を予測することができる。 The prediction model 24 generated by the prediction model generation device 20 is measured data including measured value data of production results for each past year, and at least one of observed value data of meteorological data and measured data of soil components for those years. may be a predictive model representing changes in production performance over time or averages and effects of weather and/or soil composition, generated with reference to. Such a prediction model 24 is generated by performing machine learning using measurement data as learning data, and includes at least a part of the actual measurement values of the production results up to the previous year, the forecast values of the current year's weather data, and the current year's soil At least one of the component measurement values is the input information, and the production result for the current year is the output value. According to the prediction model 24, it is possible to predict the production performance in consideration of soil components.

予測モデル24は、予測対象地点毎に生成した予測モデルであり、当該予測対象地点における前年までの前記生産成績の実測値の少なくとも一部と当該予測対象地点における当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とから、当該予測対象地点における当年の生産成績を予測してもよい。
予測モデル24は、過去のある年の生産成績の実測値と、予測モデル24により予測した、その年の生産成績の予測値との差を、その年の気象データの観測値及びその年の土壌成分の測定値の少なくとも一方に基づいて回帰する回帰分析を用いて、再生成してもよい。
The forecast model 24 is a forecast model generated for each forecast target point, and includes at least a part of the actual measurement values of the production performance up to the previous year at the forecast target point, the forecast value of the current year's weather data at the forecast target point, and the current year At least one of the measured values of the soil components in the prediction target location may be used to predict the production results for the current year.
The prediction model 24 calculates the difference between the actual measured value of the production performance of a certain past year and the predicted value of the production performance of that year predicted by the prediction model 24 by calculating the difference between the observed value of the weather data of that year and the soil of that year. It may be regenerated using a regression analysis that recurs on at least one of the component measurements.

予測モデル生成装置30が生成する予測モデル34は、予測対象日以前の気象データの観測値データ及び土壌成分の測定値データの少なくとも一方を参照して生成された、生産成績と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルであって、予測対象日の気象データの予報値及び土壌成分の測定値の少なくとも一方から生産成績を予測する予測モデルであり得る。 The prediction model 34 generated by the prediction model generation device 30 is generated by referring to at least one of the observed value data of the weather data before the target date of prediction and the measured value data of the soil components, and the production results and the weather and soil components. The prediction model may be a prediction model that expresses at least one of the effects and predicts the production performance from at least one of the predicted value of the weather data for the prediction target day and the measured value of the soil component.

このような予測モデル34は、過去のある年の生産成績の実測値と、予測モデル34を用いて予測したその年の生産成績の予測値との差が最小となるように、気象データの変数の組み合わせ及び土壌成分の変数の組み合わせを選択して生成されている。これらの変数の組み合わせは、予測対象日以前の気象データ及び土壌成分の各変数について、それぞれの変数の積算値又は平均値に基づく回帰分析により選択されてもよい。また、これらの変数の組み合わせは、上記積算値又は平均値を算出する起算日及び終算日、並びに算出日数に基づく回帰分析により選択されてもよい。さらに、これらの変数の組み合わせは、気象データ及び土壌成分の各変数について、起算日から終算日までの各日における算出日数分の積算値又は平均値に基づく回帰分析により選択されてもよい。 Such a prediction model 34 uses variables of meteorological data such that the difference between the actual production performance of a certain past year and the predicted production performance of that year predicted using the prediction model 34 is minimized. and a combination of soil component variables. A combination of these variables may be selected by regression analysis based on the integrated value or average value of each variable for each variable of weather data and soil components before the prediction target date. Also, the combination of these variables may be selected by regression analysis based on the start date and end date for calculating the integrated value or average value, and the number of days for calculation. Furthermore, the combination of these variables may be selected by regression analysis based on the integrated values or average values for the number of days calculated for each day from the start date to the end date for each variable of weather data and soil components.

また、予測モデル34は、過去の年毎の前記生産成績の実測値データとそれらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値データの少なくとも一方とを含む測定データを参照して生成した、生産成績の経年変化又は平均値と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルであって、前年までの生産成績の実測値の少なくとも一部と予測対象日の気象データの予報値及び土壌成分の測定値の少なくとも一方とから、生産成績を予測する予測モデルである。このような予測モデル34は、測定データを学習データとして機械学習を行うことによって生成することができる。このような予測モデル34は、前年までの生産成績の実測値の少なくとも一部と予測対象日の気象データの予報値及び土壌成分の測定値の少なくとも一方とが入力情報であり、当年の生産成績が出力情報である予測モデルである。 In addition, the prediction model 34 stores measurement data including at least one of actual measurement data of production results for each past year, observation data of meteorological data of those years, and measurement data of soil components of those years. A prediction model generated by reference, representing changes in production results over time or average values and the influence of at least one of weather and soil components, at least part of the actual production results up to the previous year and the target date for prediction A prediction model for predicting production results from at least one of forecast values of weather data and measured values of soil components. Such a prediction model 34 can be generated by performing machine learning using measured data as learning data. Such a prediction model 34 uses as input information at least a part of the actual measured values of the production results up to the previous year and at least one of the predicted values of the weather data for the prediction target day and the measured values of the soil components, and the production results of the current year. is the prediction model whose output information is

土壌成分が生産成績に及ぼす影響について、図16を参照して説明する。図16は、土壌水分値が生産成績に及ぼす影響を説明する図である。図16に示すグラフにおけるX軸は日付、Y軸は土壌水分値の糖度換算影響値を示している。図16に示すグラフは、土壌水分値が収穫時糖度に及ぼす影響を示している。 The effects of soil components on production performance will be described with reference to FIG. FIG. 16 is a diagram for explaining the influence of soil moisture value on production results. In the graph shown in FIG. 16, the X-axis indicates the date, and the Y-axis indicates the effect value of the soil water content converted into sugar content. The graph shown in FIG. 16 shows the influence of the soil moisture value on the sugar content at harvest.

図16に示すように、日毎に土壌水分値が収穫時糖度に影響する程度が異なっている。一例として、図16によれば、7月~9月の土壌水分値が収穫時糖度に大きな影響を与えていることが分かる。この知見は、7月及び8月の乾燥ストレスにより糖度が上昇し、それ以降はその効果が飽和するという、果樹分野における既存の知見に一致する。 As shown in FIG. 16, the extent to which the soil moisture value affects the sugar content at harvest differs from day to day. As an example, according to FIG. 16, it can be seen that the soil moisture value from July to September has a great effect on the sugar content at harvest. This finding agrees with the existing knowledge in the field of fruit trees that the sugar content increases due to drought stress in July and August, and the effect saturates thereafter.

このように、土壌成分を考慮して生産成績を予測することにより、より精度よく生産成績を予測することができる。また、気象と土壌成分との両方を考慮して生産生成器を予測することにより、予測精度のさらなる向上が期待できる。 In this way, by predicting the production results in consideration of the soil components, it is possible to predict the production results with higher accuracy. In addition, by predicting production generators in consideration of both weather and soil components, further improvement in prediction accuracy can be expected.

すなわち、本明細書に記載された各実施形態において、気象データに関する説明の全てを、(i)土壌成分の測定値データ、又は、(ii)気象データ及び土壌成分の測定値データに読み替えた実施形態についても、本発明の技術的範囲に含まれる。 That is, in each embodiment described in this specification, all the explanations about the weather data are read as (i) the measured value data of the soil components, or (ii) the measured value data of the weather data and the soil components. The form is also included in the technical scope of the present invention.

(情報量基準に基づき変数を選択する態様)
予測モデル生成装置30が生成する予測モデル34における気象データの変数の組み合わせ及び土壌成分の変数の組み合わせの選択は、情報量基準に基づいて実行してもよい。すなわち、予測モデルの生成方法において、情報量基準に基づき、気象データの変数の組み合わせ及び土壌成分の変数の組み合わせを選択して、予測モデル34を生成する態様についても、本発明の技術的範囲に含まれる。ここで、情報量基準としては、従来公知の情報量基準を用いることが可能であり、一例として、カルバック-ライブラー情報量基準、赤池情報量基準、ベイズ情報量基準等が挙げられる。
(Mode for selecting variables based on information content criteria)
The combination of weather data variables and the combination of soil component variables in the prediction model 34 generated by the prediction model generation device 30 may be selected based on information amount criteria. That is, in the prediction model generation method, a combination of weather data variables and a combination of soil component variables are selected based on the information amount criteria to generate the prediction model 34, which is also within the technical scope of the present invention. included. Here, as the information criterion, a conventionally known information criterion can be used, and examples thereof include the Kullback-Leibler information criterion, the Akaike information criterion, the Bayesian information criterion, and the like.

一例として、気象データの変数及び土壌成分の変数の組み合わせにおいて生成した各予測モデルの推定アルゴリズム内にモデル選択基準が内包されていない場合には、情報量基準に基づいて気象データの変数及び土壌成分の変数の組み合わせを選択してもよい。これにより、変数の選択処理に操作者の知識や熟練度等は影響せず、変数の選択処理を画一的に実行することができる。 As an example, if the model selection criteria are not included in the estimation algorithm of each prediction model generated in combination of weather data variables and soil component variables, meteorological data variables and soil components based on information criteria A combination of variables may be selected. As a result, the operator's knowledge, skill level, etc. do not affect the variable selection process, and the variable selection process can be uniformly executed.

(生育ステージを予測する態様)
予測装置10は、予測された生産成績と農作物の生育ステージ毎の成績を表す基準値とを比較して、農作物の生育ステージを予測してもよい。ここで、農作物の生育ステージとは、農作物の開花期、収穫期、移植期等が意図される。予測装置10における生育ステージの予測は、予測部13において実行し得る。一例として、予測装置10の予測部13は、予測モデル24により得られた糖度や酸度のような生産成績の予測値が、予め設定された基準値(しきい値)を超えた場合に、収穫基準を満たしたとして、当該生産成績が予測された時期を収穫時期であると予測する。
(Aspect of Predicting Growth Stage)
The prediction device 10 may predict the growth stage of the crop by comparing the predicted production performance with a reference value representing the performance of each growth stage of the crop. Here, the growing stage of crops means the blooming period, the harvesting period, the transplanting period, and the like of the crops. Prediction of the growth stage in prediction device 10 can be performed in prediction section 13 . As an example, the prediction unit 13 of the prediction device 10 predicts the harvest when the predicted value of the production performance such as the sugar content and acidity obtained by the prediction model 24 exceeds a preset reference value (threshold value). Assuming that the criteria are met, the time when the relevant production performance is predicted is predicted as the harvest time.

また、予測装置10は、予測モデルにおける、気象及び土壌成分の少なくとも一方の影響量を算出し、算出した影響量を参照して、農作物の生育ステージを予測してもよい。 The prediction device 10 may also calculate the influence amount of at least one of weather and soil components in the prediction model, and refer to the calculated influence amount to predict the growth stage of the crop.

農作物には、室内で苗を育ててから野外に移植される農作物や、収穫後に屋根のある日陰で貯蔵される農作物等がある。このような農作物に関して、予測モデル生成部23及び予測モデル生成部33において、生産成績に及ぼす気象又は土壌成分の影響量(影響分)を算出すれば、影響量が離散的に切り替わる日付を検出することができる。例えば、収穫後に屋根のある日陰で貯蔵されるような農作物の場合、日照時間の影響量が急激にほぼゼロになる日付が、収穫後に屋根のある貯蔵環境に移送された日に対応する。したがって、予測装置10は、このような影響量を参照することで、農作物の生育ステージを予測することができる。 Agricultural crops include crops that are grown indoors and then transplanted to the outdoors, crops that are stored in the shade with a roof after harvesting, and the like. With respect to such crops, the predictive model generation unit 23 and the predictive model generation unit 33 calculate the amount of influence (influence) of weather or soil components on production performance, and detect the date when the amount of influence discretely changes. be able to. For example, for crops that are stored in covered shade after harvest, the date on which the amount of sunshine duration abruptly drops to near zero corresponds to the date of transfer to the covered storage environment after harvest. Therefore, the prediction device 10 can predict the growth stage of crops by referring to such an influence amount.

(生育ステージを予測する予測装置)
なお、過去の年毎の農作物の実際の生育ステージの推移を表すデータと、それらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値データの少なくとも一方とを含む測定データとを参照して生成された、生育ステージの推移と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルを用いて、前年までの実際の生育ステージの推移の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とから、当年の生育ステージの推移を予測する、生育ステージの予測装置についても、本発明の技術的範囲に含まれ得る。ここで、生育ステージの推移とは、一例として、開花日、収穫日、移植日等の農作物の各生育ステージの切り替わりの日付に関する情報が意図される。
(Prediction device for predicting growth stage)
In addition, data representing the transition of the actual growth stage of agricultural products for each past year, and measurement data including at least one of observation value data of weather data of those years and measured value data of soil components of those years At least part of the actual growth stage transition up to the previous year and the current year's weather data The technical scope of the present invention also includes a growth stage prediction device that predicts the transition of the growth stage of the current year from at least one of the predicted values and the measured values of the soil components of the current year. Here, the growth stage transition is intended to be, for example, information regarding the date of switching of each growth stage of crops, such as flowering date, harvest date, transplant date, and the like.

すなわち、このような生育ステージの予測装置が用いる予測モデルは、農作物の生育ステージの推移の経年変化又は平均値と気象及び土壌成分の少なくとも一方の影響とを表す生育ステージの予測モデルである。生育ステージの予測装置は、前年までの実際の生育ステージの推移の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とを、生育ステージの予測モデルに入力して、出力される生育ステージの推移を取得する。一例として、生育ステージの予測装置は、生育ステージの予測モデルを用いることで、開花日の日付、所定の日付からの開花日までの日数、播種日から開花日までの日数等の予測結果を取得することができる。 That is, the prediction model used by such a growth stage prediction device is a growth stage prediction model that expresses changes over time or the average value of changes in the growth stages of agricultural crops and the influence of at least one of weather and soil components. The growth stage prediction device inputs at least a part of the actual growth stage transition up to the previous year and at least one of the current year's weather data forecast values and current year's measured soil component values into the growth stage prediction model. to acquire the transition of the output growth stage. As an example, the growth stage prediction device uses a growth stage prediction model to obtain prediction results such as the date of flowering, the number of days from a predetermined date to the flowering date, and the number of days from the sowing date to the flowering date. can do.

なお、本明細書に記載された実施形態において、生産成績の予測に関する説明の全てを生育ステージの推移の予測に読み替えた実施形態についても、本発明の技術的範囲に含まれる。 In addition, in the embodiment described in this specification, the technical scope of the present invention also includes an embodiment in which all the explanations related to the prediction of the production performance are replaced with the prediction of the transition of the growth stage.

(生産成績のポテンシャル値を算出する態様)
予測装置の他の形態について、図17を参照して説明する。なお、説明の便宜上、上述した予測システム100及び予測装置10にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を繰り返さない。図17は、本発明の他の形態に係る予測システム200の要部構成を示すブロック図である。予測システム200は、予測装置210の制御部211に、入力値生成部212、出力値取得部213、設定部214、及び算出部215を備える点において、図1に示す予測システム100と異なっている。
(Method of calculating potential value of production results)
Another form of prediction device will be described with reference to FIG. For convenience of explanation, members having the same functions as the members explained in the prediction system 100 and the prediction device 10 described above are denoted by the same reference numerals, and the explanation thereof will not be repeated. FIG. 17 is a block diagram showing the main configuration of a prediction system 200 according to another embodiment of the invention. The prediction system 200 differs from the prediction system 100 shown in FIG. 1 in that the control unit 211 of the prediction device 210 includes an input value generation unit 212, an output value acquisition unit 213, a setting unit 214, and a calculation unit 215. .

入力値生成部212は、過去の気象データの観測値及び土壌成分の測定値の少なくとも一方に基づき、気象データの観測値及び土壌成分の測定値の少なくとも一方のシミュレーション入力値を生成する(入力値を生成する工程)。入力値生成部212は、入力データ取得部12が取得した気象データの観測値及び土壌成分の測定値の少なくとも一方において、一日単位でデータを入れ替え、シミュレーション入力値を生成する。データの入れ替えは、同一場所における過去の別の年の同じ日のデータと入れ替えるように実行する。すなわち、入力値生成部212が生成するシミュレーション入力値は、複数年分のデータを一日単位でつぎはぎした気象データの観測値及び土壌成分の測定値の少なくとも一方のデータであり得る。入力値生成部212は、生成したシミュレーション入力値を、出力値取得部213へ出力する。 The input value generation unit 212 generates a simulation input value for at least one of the observed meteorological data and the measured value of the soil component based on at least one of the observed value of the past meteorological data and the measured value of the soil component (input value step of generating). The input value generation unit 212 generates a simulation input value by replacing data in at least one of the observed values of the weather data and the measured values of the soil components acquired by the input data acquisition unit 12 on a daily basis. The replacement of the data is performed so as to replace the data on the same day of another past year in the same location. In other words, the simulation input values generated by the input value generation unit 212 can be at least one of observed values of weather data and measured values of soil components obtained by patching data for a plurality of years on a daily basis. The input value generation unit 212 outputs the generated simulation input values to the output value acquisition unit 213 .

出力値取得部213は、予測モデル24又は予測モデル34を用いて、シミュレーション入力値から予測された、生産成績のシミュレーション出力値を取得する(出力値を取得する工程)。出力値取得部213は、予測モデル24又は予測モデル34に、シミュレーション入力値を入力し、出力される生産成績の予測値をシミュレーション出力値として取得する。出力値取得部213は、取得したシミュレーション出力値を設定部214へ出力する。 The output value acquiring unit 213 acquires the simulation output value of the production performance predicted from the simulation input value using the prediction model 24 or the prediction model 34 (step of acquiring the output value). The output value acquiring unit 213 inputs a simulation input value to the prediction model 24 or the prediction model 34, and acquires the output predicted value of production performance as a simulation output value. The output value acquisition unit 213 outputs the acquired simulation output value to the setting unit 214 .

設定部214は、得られたシミュレーション出力値が最も高くなるときのシミュレーション入力値をポテンシャル値に設定する(設定する工程)。設定部214は、さまざまなシミュレーション入力値に対して得られたシミュレーション出力値を比較して、シミュレーション出力値が最も高くなる時のシミュレーション入力値をポテンシャル値に設定する。すなわち、設定部214は、つぎはぎの態様が異なる複数の気象データ又は土壌成分のデータを用いて予測された生産成績を比較して、ポテンシャル値を設定する。設定部214は設定したポテンシャル値を算出部215へ出力する。 The setting unit 214 sets the simulation input value when the obtained simulation output value is the highest as the potential value (setting step). The setting unit 214 compares the simulation output values obtained for various simulation input values, and sets the simulation input value when the simulation output value is the highest as the potential value. That is, the setting unit 214 sets the potential value by comparing production performance predicted using a plurality of pieces of meteorological data or soil component data with different patches. The setting unit 214 outputs the set potential value to the calculation unit 215 .

算出部215は、ポテンシャル値が得られたときのシミュレーション入力値を最適環境値として算出する(算出する工程)。すなわち、算出部215は、最もよい生産成績が得られたときの気象データ又は土壌成分のデータを、当該生産成績を得るための気象条件及び土壌成分の条件の少なくとも一方を、最適環境値として算出する。 The calculation unit 215 calculates the simulation input value when the potential value is obtained as the optimum environment value (step of calculating). That is, the calculation unit 215 calculates the weather data or the soil component data when the best production results are obtained, and calculates at least one of the weather conditions and the soil component conditions for obtaining the best production results as the optimum environmental value. do.

予測装置210においては、気象データの観測値及び土壌成分の測定値の少なくとも一方について、同一場所における過去の別の年の同じ日のデータと入れ替えたシミュレーション入力値を仮想環境条件とし、仮想環境条件により予測される生産成績をポテンシャル値とする。このようなシミュレーション入力値は、例えば、「2010年3月3日の最高気温、日照時間及び降水量」を「2011年3月3日の最高気温、日照時間及び降水量」と入れ替えるというような処理を生育期間のすべての日付において実行することで、生成する。このように、仮想環境条件ごとに生産成績を比較することで、最も理想的な環境条件を予測することができる。 In the prediction device 210, for at least one of the observed meteorological data and the measured value of the soil component, the simulation input value that is replaced with data on the same day in another past year at the same place is set as the virtual environment condition. The potential value is the production performance predicted by Such simulation input values are, for example, replacing "maximum temperature, sunshine hours and precipitation on March 3, 2010" with "maximum temperature, sunshine hours and precipitation on March 3, 2011". Generated by running the treatment on all dates of the growing season. In this way, by comparing the production results for each virtual environment condition, it is possible to predict the most ideal environment condition.

なお、入力値生成部212は、予測モデル生成部23及び予測モデル生成部33から得られる生産成績に及ぼす気象又は土壌成分の影響分に関する情報から、どの時期のどの気象条件又は土壌成分の条件を入れ替えるかを決定してもよい。また、入力値生成部212は、生産成績に及ぼす気象又は土壌成分の影響分に関する情報に基づき、入れ替えるデータを決定するようなアルゴリズムを用いて、シミュレーション入力値を生成してもよい。 Note that the input value generating unit 212 determines which weather conditions or soil component conditions at what time based on the information on the influence of the weather or soil components on the production performance obtained from the prediction model generating unit 23 and the predictive model generating unit 33. You can decide to replace it. The input value generation unit 212 may also generate simulation input values using an algorithm that determines data to be replaced based on information about the influence of weather or soil components on production performance.

本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。 The present invention is not limited to the above-described embodiments, but can be modified in various ways within the scope of the claims, and can be obtained by appropriately combining technical means disclosed in different embodiments. is also included in the technical scope of the present invention.

10、210 予測装置
13 予測部
20 予測モデル生成装置
23 予測モデル生成部(生成部)
30 予測モデル生成装置
33 予測モデル生成部
100、200 予測システム
10, 210 prediction device 13 prediction unit 20 prediction model generation device 23 prediction model generation unit (generation unit)
30 prediction model generation device 33 prediction model generation unit 100, 200 prediction system

Claims (10)

農作物の生産成績を予測する予測モデルの生成方法であって、
過去の年毎の前記生産成績の実測値データを参照して、前記生産成績の経年変化又は平均値を表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部から当年の前記生産成績を予測する予測モデルを生成する工程を包含する、生成方法。
A method for generating a prediction model for predicting the production performance of agricultural products, comprising:
A prediction model that represents changes over time or an average value of the production performance by referring to the actual measurement data of the production performance for each past year, wherein at least a part of the production performance of the previous year is converted to the current year A method of generating, comprising generating a predictive model that predicts the production performance.
前記生成する工程においては、前記実測値データを学習データとして機械学習を行うことによって、前年までの前記生産成績の実測値の少なくとも一部が入力情報であり、当年の前記生産成績が出力情報である予測モデルを生成する、請求項1に記載の生成方法。 In the generating step, by performing machine learning using the measured value data as learning data, at least part of the measured values of the production results up to the previous year are input information, and the production results of the current year are output information. 2. The generating method of claim 1, wherein a predictive model is generated. 前記生成する工程においては、過去の年毎の前記生産成績の実測値データとそれらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値データの少なくとも一方とを含む測定データを参照して、前記生産成績の経年変化又は平均値と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とから、当年の前記生産成績を予測する予測モデルを生成する、請求項1に記載の生成方法。 In the generating step, measurement data including at least one of actual measurement data of production results for each past year, observation data of weather data for those years, and measurement data of soil components for those years. Reference, a prediction model representing changes over time or an average value of the production performance and the influence of at least one of weather and soil components, wherein at least a portion of the production performance measured up to the previous year and weather data for the current year 2. The generation method according to claim 1, wherein the prediction model for predicting the production performance for the current year is generated from at least one of the predicted value of , and the measured value of the soil component for the current year. 前記生成する工程においては、前記測定データを学習データとして機械学習を行うことによって、前年までの前記生産成績の実測値の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とが入力情報であり、当年の前記生産成績が出力情報である予測モデルを生成する、請求項3に記載の生成方法。 In the generating step, by performing machine learning using the measured data as learning data, at least a part of the measured production results up to the previous year, forecast values of weather data for the current year, and measured values of soil components for the current year 4. The generating method according to claim 3, wherein at least one of is input information, and the production result for the current year is output information. 前記生成する工程においては、予測対象地点における過去の年毎の前記生産成績の実測値データと当該予測対象地点におけるそれらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値の少なくとも一方とを含む測定データを参照して、予測対象地点毎の予測モデルを生成する、請求項3又は4に記載の生成方法。 In the generating step, actual measurement data of the production results for each past year at the target point for prediction, observed value data of meteorological data for those years at the target point for prediction, and measured values of soil components for those years 5. The generation method according to claim 3 or 4, wherein a prediction model for each prediction target point is generated with reference to measurement data including at least one of 過去のある年の前記生産成績の実測値と、前記予測モデルを用いて予測した、その年の前記生産成績の予測値との差が最小になるように、前記予測モデルを再生成する工程をさらに包含する、請求項3から5のいずれか1項に記載の生成方法。 a step of regenerating the prediction model so that the difference between the measured value of the production performance for a certain past year and the predicted value of the production performance for that year predicted using the prediction model is minimized; 6. The production method of any one of claims 3-5, further comprising. 前記再生成する工程においては、過去のある年の前記生産成績の実測値と、前記予測モデルにより予測した、その年の前記生産成績の予測値との差を、その年の気象データの観測値及びその年の土壌成分の測定値の少なくとも一方に基づいて回帰する回帰分析を用いて、前記予測モデルを再生成する、請求項6に記載の生成方法。 In the regenerating step, the difference between the actual production performance of a certain past year and the predicted production performance of that year predicted by the prediction model is calculated as the observed meteorological data of that year. 7. The generation method of claim 6, wherein the predictive model is regenerated using a regression analysis that regresses based on at least one of measurements of soil constituents for that year. 農作物の生産成績を予測する予測モデルの生成装置であって、
過去の年毎の前記生産成績の実測値データを参照して、前記生産成績の経年変化又は平均値を表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部から当年の前記生産成績を予測する予測モデルを生成する生成部を備えた、生成装置。
A prediction model generation device for predicting the production performance of agricultural products,
A prediction model that represents changes over time or an average value of the production performance by referring to the actual measurement data of the production performance for each past year, wherein at least a part of the production performance of the previous year is converted to the current year A generation device comprising a generation unit that generates a prediction model for predicting the production performance.
農作物の生産成績を予測する予測モデルの生成プログラムであって、
過去の年毎の前記生産成績の実測値データを参照して、前記生産成績の経年変化又は平均値を表す予測モデルであって、前年までの前記生産成績の実測値の少なくとも一部から当年の前記生産成績を予測する予測モデルを生成するステップを包含する、生成プログラム。
A prediction model generation program for predicting the production performance of agricultural products,
A prediction model that represents changes over time or an average value of the production performance by referring to the actual measurement data of the production performance for each past year, wherein at least a part of the production performance of the previous year is converted to the current year A generating program comprising generating a predictive model that predicts the production performance.
農作物の生育ステージを予測する予測モデルの生成方法であって、
過去の年毎の農作物の実際の生育ステージの推移を表すデータと、それらの年の気象データの観測値データ及びそれらの年の土壌成分の測定値データの少なくとも一方とを含む測定データを参照して、前記生育ステージの推移と気象及び土壌成分の少なくとも一方の影響とを表す予測モデルであって、前年までの実際の前記生育ステージの推移の少なくとも一部と当年の気象データの予報値及び当年の土壌成分の測定値の少なくとも一方とから、当年の前記生育ステージの推移を予測する予測モデルを生成する工程を包含する、生成方法。
A method for generating a prediction model for predicting the growth stage of agricultural crops, comprising:
Refers to measurement data that includes data representing changes in the actual growth stage of crops for each year in the past, and at least one of observed meteorological data for those years and measured value data for soil components for those years. A prediction model representing the transition of the growth stage and the influence of at least one of weather and soil components, wherein at least a part of the actual transition of the growth stage up to the previous year, the forecast value of the current year's weather data, and the current year generating a predictive model that predicts the transition of said growth stage for the current year from at least one of the measured values of the soil components of .
JP2022033768A 2021-03-05 2022-03-04 Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program Pending JP2022136057A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2021035833 2021-03-05
JP2021035833 2021-03-05

Publications (1)

Publication Number Publication Date
JP2022136057A true JP2022136057A (en) 2022-09-15

Family

ID=83232405

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022033768A Pending JP2022136057A (en) 2021-03-05 2022-03-04 Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program

Country Status (1)

Country Link
JP (1) JP2022136057A (en)

Similar Documents

Publication Publication Date Title
ur Rahman et al. Multi-model projections of future climate and climate change impacts uncertainty assessment for cotton production in Pakistan
US11617313B2 (en) Controlling agricultural production areas
Dias et al. Evaluation of three sugarcane simulation models and their ensemble for yield estimation in commercially managed fields
JP2022136058A (en) Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program
Vezy et al. DynACof: A process-based model to study growth, yield and ecosystem services of coffee agroforestry systems
KR20180086738A (en) Apparatus for providing guide information for growing crops, and method thereof
Bannayan et al. A stochastic modelling approach for real-time forecasting of winter wheat yield
Höglind et al. Process-based simulation of growth and overwintering of grassland using the BASGRA model
Seyoum et al. Application of the APSIM model to exploit G× E× M interactions for maize improvement in Ethiopia
JP2021056573A (en) Crop growth prediction program, crop growth prediction method, and crop growth prediction device
JP2019030253A (en) Culture support device, culture support method, computer program and recording medium
JP2019170359A (en) Plant cultivation result prediction system
KR20160040672A (en) Cultivation assistance device, cultivation assistance method, and recording medium for storing program
AU2023203307A1 (en) Controlling Agricultural Production Areas
Chopard et al. Development of a decision support system to evaluate crop performance under dynamic solar panels
Dinh et al. Statistical analysis of the weather impact on Robusta Coffee Yield in Vietnam
Hopf et al. Development and improvement of the CROPGRO-Strawberry model
JP2021179983A (en) Method and program for predicting growth, and information processing device
JP7208503B2 (en) Machine learning program, machine learning method and machine learning apparatus
CN113010529A (en) Crop management method and device based on knowledge graph
JP2022136057A (en) Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program
Bhengra et al. Calibration and validation study of sugarcane (DSSAT-CANEGRO V4. 6.1) model over North Indian region
JP2022136056A (en) Prediction method for predicting crop production performance, prediction apparatus, and prediction program
Susilo et al. Decision Tree-Based Bok Choy Growth Prediction Model for Smart Farm
Spiegel et al. Step-by-step development of a model simulating returns on farm from investments: the example of hazelnut plantation in Italy

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230622

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20230726

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231024

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231219

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240319

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240515