JP6849958B2 - Prediction method, prediction program, and prediction device - Google Patents

Prediction method, prediction program, and prediction device Download PDF

Info

Publication number
JP6849958B2
JP6849958B2 JP2016055340A JP2016055340A JP6849958B2 JP 6849958 B2 JP6849958 B2 JP 6849958B2 JP 2016055340 A JP2016055340 A JP 2016055340A JP 2016055340 A JP2016055340 A JP 2016055340A JP 6849958 B2 JP6849958 B2 JP 6849958B2
Authority
JP
Japan
Prior art keywords
year
target
crop
prediction
date
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016055340A
Other languages
Japanese (ja)
Other versions
JP2017163956A (en
Inventor
栄輔 北
栄輔 北
望 久保
望 久保
詠梨 吉田
詠梨 吉田
佐藤 仁志
仁志 佐藤
公佑 神部
公佑 神部
英己 北野
英己 北野
俊策 西内
俊策 西内
笑 亀岡
笑 亀岡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tokai National Higher Education and Research System NUC
Original Assignee
Tokai National Higher Education and Research System NUC
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 Tokai National Higher Education and Research System NUC filed Critical Tokai National Higher Education and Research System NUC
Priority to JP2016055340A priority Critical patent/JP6849958B2/en
Publication of JP2017163956A publication Critical patent/JP2017163956A/en
Application granted granted Critical
Publication of JP6849958B2 publication Critical patent/JP6849958B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本開示は、作物の収穫状況を予測する技術に関する。 The present disclosure relates to techniques for predicting crop yields.

ニューラルネットワーク又はカオス分析を用いて、農作物の収穫量を予測する技術が知られている(特許文献1参照)。該技術では、農作物の栽培計画に基づき、第1予測収穫量が算出される。また、ニューラルネットワーク等を用いて、天候データに基づく補正値が算出される。そして、第1予測収穫量と補正値とに基づき農作物の収穫量が算出される。 A technique for predicting the yield of agricultural products using a neural network or chaos analysis is known (see Patent Document 1). In this technique, the first predicted yield is calculated based on the cultivation plan of the crop. In addition, a correction value based on the weather data is calculated using a neural network or the like. Then, the yield of the agricultural product is calculated based on the first predicted yield and the corrected value.

特開2002−136223号公報JP-A-2002-136223

しかしながら、気象の変動が大きい場合には、農作物の収穫量を精度良く予測できない恐れがある。
本開示の一側面は、作物の収穫状況を精度良く予測するのが好ましい。
However, if the weather fluctuates greatly, there is a risk that the yield of crops cannot be predicted accurately.
One aspect of the present disclosure is preferably to accurately predict the crop yield.

本開示の一側面に係る予測方法は、予め定められた地域における予め定められた作物の収穫状況を予測する。該予測方法では、作物が予め定められた生長段階に到達する到達時期を、該到達時期を含む年である対象年と、対象年よりも前の任意の年とのうちの双方又は一方における地域の気象情報に基づき予測する。さらに、該予測方法では、対象年における収穫状況を、対象年と任意の年とのうちの双方又は一方における到達時期を基準として定められた測定時期の気象情報に基づき予測する。 The prediction method according to one aspect of the present disclosure predicts a predetermined crop yield in a predetermined area. In the prediction method, the arrival time when the crop reaches a predetermined growth stage is determined by the target year, which is the year including the arrival time, and any year before the target year, or one of the regions. Predict based on the weather information of. Further, in the prediction method, the harvest situation in the target year is predicted based on the meteorological information of the measurement time determined based on the arrival time in both or one of the target year and an arbitrary year.

作物が予め定められた生長段階に達する時期、若しくは、該時期の前又は後の時期の気象が、作物の生長に相対的に大きな影響を与える場合がある。特に、気象の変動が激しい場合、このような時期の気象が例年と大きく異なる可能性があり、その結果、収穫状況が例年と比べて大きく変動する恐れがある。 The weather at the time when the crop reaches a predetermined growth stage, or before or after that time, may have a relatively large effect on the growth of the crop. In particular, when the weather fluctuates drastically, the weather at such times may differ significantly from the usual year, and as a result, the harvest situation may fluctuate significantly compared to the usual year.

これに対し、上記予測方法によれば、作物が予め定められた生長段階に到達する到達時期が予測される。そして、予測された到達時期を基準として定められた測定時期の気象情報に基づき、作物の収穫状況が予測される。このため、作物の生長に相対的に大きな影響を与える時期の気象に基づき、収穫状況を予測することが可能となる。したがって、作物の収穫状況を精度良く予測することが可能となる。 On the other hand, according to the above prediction method, the arrival time when the crop reaches a predetermined growth stage is predicted. Then, the harvest situation of the crop is predicted based on the meteorological information of the measurement time determined based on the predicted arrival time. Therefore, it is possible to predict the harvest situation based on the weather at the time when the crop growth is relatively significantly affected. Therefore, it is possible to accurately predict the harvest status of crops.

なお、上記構成において、対象年における作物を植えた日を、開始日とし、到達時期を、対象年と任意の年とのうちの双方又は一方における開始日を基準に定められた日の気象情報に基づき予測しても良い。また、測定時期とは、対象年と任意の年とのうちの双方又は一方における到達時期を基準として定められた日であっても良い。 In the above configuration, the day when the crop was planted in the target year is set as the start date, and the arrival time is the meteorological information of the day determined based on the start date in both or one of the target year and any year. It may be predicted based on. Further, the measurement time may be a day determined based on the arrival time in both or one of the target year and the arbitrary year.

作物の栽培期間における特定の日の気象が、作物の生長に相対的に大きな影響を与える場合がある。これに対し、上記構成によれば、開始日を基準に定められた日の気象情報に基づき、到達時期が予測される。そして、予測された到達時期を基準として定められた日である測定時期の気象情報に基づき、収穫状況が予測される。このため、作物の生長への影響が相対的に大きい特定の日の気象に基づき、到達時期及び収穫状況を予測することが可能となる。したがって、収穫状況を精度良く予測することが可能となる。 The weather on a particular day during the growing season of the crop can have a relatively large effect on the growth of the crop. On the other hand, according to the above configuration, the arrival time is predicted based on the meteorological information of the day determined based on the start date. Then, the harvest situation is predicted based on the meteorological information of the measurement time, which is the day determined based on the predicted arrival time. Therefore, it is possible to predict the arrival time and the harvest situation based on the weather on a specific day, which has a relatively large effect on the growth of crops. Therefore, it is possible to accurately predict the harvest situation.

また、上記構成において、到達時期とは、作物の穂又は花が予め定められた生長段階に到達した時期であっても良い。
作物の穂又は花が予め定められた生長段階に到達した時期、若しくは、該時期の前又は後の時期の気象が、作物の生長に相対的に大きな影響を与える場合がある。これに対し、上記構成によれば、作物がこのような生長段階に到達した時期が、到達時期として予測される。そして、予測された到達時期を基準に定められた測定時期の気象情報により、収穫状況が予測される。このため、作物の収穫状況を精度良く予測することが可能となる。
Further, in the above configuration, the arrival time may be the time when the spikes or flowers of the crop reach a predetermined growth stage.
The weather at the time when the spikes or flowers of the crop reach a predetermined growth stage, or before or after that time, may have a relatively large effect on the growth of the crop. On the other hand, according to the above configuration, the time when the crop reaches such a growth stage is predicted as the time of arrival. Then, the harvest situation is predicted by the weather information of the measurement time determined based on the predicted arrival time. Therefore, it is possible to accurately predict the harvest status of crops.

また、上記構成において、作物とは、稲であっても良い。
このような構成によれば、稲の収穫状況を精度良く予測することが可能となる。
また、本開示の一側面は、上述した予測方法をコンピュータに実行させる予測プログラムに関するものである。
Further, in the above configuration, the crop may be rice.
With such a configuration, it is possible to accurately predict the rice harvest situation.
Further, one aspect of the present disclosure relates to a prediction program that causes a computer to execute the above-mentioned prediction method.

また、本開示の一側面は、上述した予測方法を実行する予測装置に関するものである。
このような予測プログラム又は予測装置によれば、作物の収穫状況を精度良く予測することが可能となる。
Further, one aspect of the present disclosure relates to a prediction device that executes the above-mentioned prediction method.
According to such a prediction program or a prediction device, it is possible to accurately predict the harvest situation of the crop.

図1Aは、本実施形態の予測システムのブロック図である。図1Bは、本実施形態の予測処理のフローチャートである。FIG. 1A is a block diagram of the prediction system of the present embodiment. FIG. 1B is a flowchart of the prediction process of the present embodiment. 予測処理にて用いられるニューラルネットワークの説明図である。It is explanatory drawing of the neural network used in the prediction processing.

以下、図面を参照しながら、発明を実施するための形態を説明する。
[1.構成]
図1Aに示す予測システム1は、サーバ10と、端末20と、測定装置30とを有する。これらは、インターネット等の通信網に接続されている。
Hereinafter, modes for carrying out the invention will be described with reference to the drawings.
[1. Constitution]
The prediction system 1 shown in FIG. 1A includes a server 10, a terminal 20, and a measuring device 30. These are connected to a communication network such as the Internet.

サーバ10は、通信網を介して端末20及び測定装置30等と通信を行う。サーバは、記憶部11、制御部、表示部、及び操作部等を備える。なお、本実施形態では、一例として、PCがサーバ10として用いられる。しかし、これに限らず、他の電子装置をサーバ10として用いても良い。 The server 10 communicates with the terminal 20, the measuring device 30, and the like via the communication network. The server includes a storage unit 11, a control unit, a display unit, an operation unit, and the like. In this embodiment, as an example, a PC is used as the server 10. However, the present invention is not limited to this, and another electronic device may be used as the server 10.

表示部は、液晶ディスプレイ等を有し、各種情報を表示する。また、操作部は、キーボードやマウス等を有し、各種操作を受け付ける。
記憶部11は、フラッシュメモリやHDD等の書き換え可能な不揮発性の記憶装置を有する。記憶部11は、非遷移的実体的記録媒体に相当し、各種データが記憶されている。記憶部11には、サーバ10にインストールされたアプリケーションである予測プログラムが記憶されている。予測プログラムは、予め定められた地域(以後、対象地域)における予め定められた農地(以後、対象農地)での作物(以後、対象作物)の収穫状況を予測する。また、記憶部11には、気象履歴及び栽培履歴が保存されている。気象履歴は、対象地域にて測定された複数の気象情報を有する。また、栽培履歴は、過去における対象地域の農地での対象作物の栽培に関する情報を示す。なお、栽培履歴は、一例として、過去における対象農地での対象作物の栽培に関する情報を示しても良い。
The display unit has a liquid crystal display or the like and displays various information. In addition, the operation unit has a keyboard, a mouse, and the like, and accepts various operations.
The storage unit 11 has a rewritable non-volatile storage device such as a flash memory or an HDD. The storage unit 11 corresponds to a non-transitional substantive recording medium and stores various data. The storage unit 11 stores a prediction program, which is an application installed on the server 10. The prediction program predicts the harvest status of crops (hereinafter, target crops) in a predetermined farmland (hereinafter, target farmland) in a predetermined area (hereinafter, target area). In addition, the weather history and cultivation history are stored in the storage unit 11. The meteorological history has a plurality of meteorological information measured in the target area. In addition, the cultivation history shows information on the cultivation of the target crop on the agricultural land in the target area in the past. As an example, the cultivation history may show information on the cultivation of the target crop on the target agricultural land in the past.

制御部は、CPUと、ROM、RAM、及びフラッシュメモリ等の半導体メモリとを有する周知のコンピュータを備える。サーバ10の機能は、CPUが半導体メモリに記憶されたプログラムに従い動作することで実現される。また、プログラムの実行により、プログラムに対応する方法が実行される。半導体メモリには、予測プログラムがロードされる。なお、サーバ10の機能を実現する方法は、ソフトウェアに限らない。該機能の一部又は全部を、サーバ10の論理回路等(換言すれば、ハードウェア)により実現してもよい。 The control unit includes a well-known computer having a CPU and a semiconductor memory such as a ROM, RAM, and a flash memory. The function of the server 10 is realized by operating the CPU according to a program stored in the semiconductor memory. In addition, by executing the program, the method corresponding to the program is executed. The prediction program is loaded in the semiconductor memory. The method of realizing the function of the server 10 is not limited to software. A part or all of the functions may be realized by a logic circuit or the like (in other words, hardware) of the server 10.

また、サーバ10は、出穂日予測部12と、収穫量予測部13と、品質予測部14とを有する。これらは、制御部が予測プログラムに従い動作することで実現される機能に相当する。これらの部位の機能については、後述する。 Further, the server 10 has a heading date prediction unit 12, a yield prediction unit 13, and a quality prediction unit 14. These correspond to the functions realized by the control unit operating according to the prediction program. The functions of these parts will be described later.

また、測定装置30は、対象地域の気象を測定し、対象地域における日ごとの気象情報を生成する。なお、測定装置30は、対象地域における対象農地の気象を測定し、対象農地における日ごとの気象情報を生成しても良い。具体的には、気象情報として、例えば、一日における平均気温、最高気温、最低気温、日照時間、降水量、及び気温日格差等の種類が設けられていても良い。なお、気温日格差とは、その日における最高気温と最低気温の差である。測定装置30は、サーバ10に対し、随時気象情報を送信する。一方、サーバ10は、測定装置30から受信した気象情報を、気象履歴の一部として記憶部11に記憶する。 In addition, the measuring device 30 measures the weather in the target area and generates daily weather information in the target area. The measuring device 30 may measure the weather of the target farmland in the target area and generate daily weather information of the target farmland. Specifically, as the meteorological information, for example, types such as average temperature, maximum temperature, minimum temperature, sunshine duration, precipitation amount, and temperature-daily difference in one day may be provided. The daily temperature difference is the difference between the maximum temperature and the minimum temperature on that day. The measuring device 30 transmits weather information to the server 10 at any time. On the other hand, the server 10 stores the weather information received from the measuring device 30 in the storage unit 11 as a part of the weather history.

この他にも、測定装置30からサーバ10に対し、気温、降水量、及び日照状態等の測定結果をリアルタイムに送信しても良い。そして、サーバ10にて、測定装置30から受信した測定結果に基づき、上述した気象情報を生成しても良い。 In addition, the measuring device 30 may transmit the measurement results such as the temperature, precipitation, and sunshine condition to the server 10 in real time. Then, the server 10 may generate the above-mentioned meteorological information based on the measurement result received from the measuring device 30.

また、端末20は、一例として、PCやスマートフォン等の電子装置として構成されている。端末20は、ユーザからの指示に応じてサーバ10にアクセスし、対象作物の収穫状況等の予測結果を取得する。なお、端末20は、さらに、対象地域の気象履歴又は栽培履歴等といった、対象作物の栽培に関する情報をサーバ10から取得しても良い。そして、端末20は、ブラウザ等により、対象作物の収穫状況等の予測結果、又は、対象作物の栽培に関する情報等を表示する。 Further, the terminal 20 is configured as an electronic device such as a PC or a smartphone as an example. The terminal 20 accesses the server 10 in response to an instruction from the user, and acquires a prediction result such as a harvest status of the target crop. The terminal 20 may further acquire information on the cultivation of the target crop, such as the weather history or cultivation history of the target area, from the server 10. Then, the terminal 20 displays the prediction result such as the harvest status of the target crop or the information regarding the cultivation of the target crop by using a browser or the like.

[2.処理]
次に、予測プログラムによる処理について説明する。本実施形態では、予測プログラムは、一例として、予め定められた品種の稲を対象作物として、対象年における収穫状況を予測する。なお、稲とは、水稲であっても良いし、陸稲であっても良い。また、収穫状況とは、一例として、対象作物の収穫量、及び、対象作物から収穫された収穫物の品質であっても良い。また、収穫物の品質とは、対象作物が稲の場合であれば、例えば、収穫物である米の等級に関する情報であっても良い。米の等級は、米に整粒が含まれる割合、及び、米における水分の含有率等に基づき定められる。なお、整粒とは、損傷を受けた米粒(具体的には、被害粒、死米、及び着色米)以外の米粒である。収穫物である米に整粒が含まれる割合、又は、該米の水分の含有率等を、対象作物である稲の品質としても良い。この他にも、例えば、収穫物の味、又は、収穫物にどの程度の栄養が含まれるか等を、品質としても良い。また、例えば、対象作物の収穫時期等を、収穫状況として予測しても良い。
[2. processing]
Next, the processing by the prediction program will be described. In the present embodiment, as an example, the prediction program predicts the harvest situation in the target year by using a predetermined variety of rice as the target crop. The rice may be paddy rice or upland rice. Further, the harvest status may be, for example, the yield of the target crop and the quality of the harvested product harvested from the target crop. Further, the quality of the harvested product may be information on the grade of the harvested rice, for example, when the target crop is rice. The grade of rice is determined based on the ratio of sized rice contained in the rice and the water content in the rice. The sizing is rice grains other than damaged rice grains (specifically, damaged grains, dead rice, and colored rice). The ratio of sized rice contained in the harvested rice, the water content of the rice, and the like may be used as the quality of the rice as the target crop. In addition to this, for example, the taste of the harvested product, the amount of nutrition contained in the harvested product, and the like may be used as the quality. Further, for example, the harvest time of the target crop may be predicted as the harvest status.

以後、対象年の直前の年を前年と記載する。記憶部11に記憶されている気象履歴は、過去から対象年にかけての複数年(以後、記録期間)における複数の気象情報を有する。以後、前年の気象履歴を前年気象履歴と記載し、対象年の気象履歴を対象年気象履歴と記載する。また、記憶部11に記憶されている栽培履歴は、記録期間における対象地域の農地での対象作物である稲の栽培に関する情報を示す。栽培履歴は、記録期間における各年の田植日、出穂日、収穫日、収穫量、及び、品質等を示す。なお、栽培記歴は、操作部を介して入力された情報に基づき生成されても良い。気象履歴及び栽培履歴は、収穫状況の予測に用いられる。 Hereinafter, the year immediately before the target year will be described as the previous year. The meteorological history stored in the storage unit 11 has a plurality of meteorological information for a plurality of years (hereinafter, a recording period) from the past to the target year. Hereinafter, the weather history of the previous year will be described as the previous year's weather history, and the weather history of the target year will be described as the target year's weather history. In addition, the cultivation history stored in the storage unit 11 indicates information on the cultivation of rice, which is the target crop, on the agricultural land in the target area during the recording period. The cultivation history shows the rice planting date, heading date, harvest date, harvest amount, quality, etc. of each year during the recording period. The cultivation history may be generated based on the information input via the operation unit. The weather history and cultivation history are used to predict the harvest situation.

予測プログラムの出穂日予測部12は、出穂日を予測する。出穂日とは、農地で栽培されている稲における出穂した稲が占める割合が、予め定められた閾値に達した日である。なお、該閾値は、例えば、40〜50%の値であっても良い。出穂日の予測には、予め定められた田植基準日における予め定められた種類の気象情報が、複数用いられる。田植基準日とは、対象年において、田植日を基準に定められた日である。田植日とは、対象年において、対象農地にて対象作物である稲の田植が行われた日(換言すれば、対象作物である稲を植えた日)である。出穂日の予測に用いられる複数の気象情報は、対象年気象履歴又は前年気象履歴から抽出される。なお、出穂日は、前年気象履歴又は対象年気象履歴における1つの田植基準日における1つの種類の気象情報に基づき予測されても良い。 The heading date prediction unit 12 of the prediction program predicts the heading date. The heading date is the day when the ratio of the headed rice to the rice cultivated on the agricultural land reaches a predetermined threshold value. The threshold value may be, for example, a value of 40 to 50%. A plurality of predetermined types of meteorological information on a predetermined rice planting reference date are used for predicting the heading date. The rice planting reference date is the date set based on the rice planting date in the target year. The rice planting date is the day when the rice, which is the target crop, was planted on the target farmland in the target year (in other words, the day when the rice, which is the target crop, was planted). A plurality of meteorological information used for predicting the heading date is extracted from the target year meteorological history or the previous year meteorological history. The heading date may be predicted based on one type of meteorological information on one rice planting reference date in the previous year's meteorological history or the target year's meteorological history.

そして、収穫量予測部13は、出穂日に基づき稲の収穫量を予測し、品質予測部14は、出穂日に基づき稲の品質を予測する。これらの予測には、予め定められた出穂基準日における予め定められた種類の気象情報が、複数用いられる。出穂基準日とは、対象年において、予測された出穂日を基準に定められた日である。このような複数の気象情報は、対象年気象履歴又は前年気象履歴から抽出される。なお、稲の収穫量又は品質は、前年気象履歴又は対象年気象履歴における1つの出穂基準日における1つの種類の気象情報に基づき予測されても良い。 Then, the yield prediction unit 13 predicts the yield of rice based on the heading date, and the quality prediction unit 14 predicts the quality of rice based on the heading date. A plurality of predetermined types of meteorological information on a predetermined heading reference date are used for these predictions. The heading reference date is a date determined based on the predicted heading date in the target year. Such a plurality of meteorological information is extracted from the target year's meteorological history or the previous year's meteorological history. The yield or quality of rice may be predicted based on one type of meteorological information on one heading reference date in the previous year's meteorological history or the target year's meteorological history.

以下では、対象作物である稲の収穫量及び品質を予測する予測処理について説明する(図1B)。本処理は、予測プログラムに従い動作するサーバ10により実行される。なお、予測処理では、様々な地域を対象地域としたり、他の品種の稲を対象作物としたり、様々な作物(例えば、麦等の穂を有する作物、又は、トマト又は枝豆等の野菜)の様々な品種を対象作物としたりして収穫状況を予測することができる。穂を有さない野菜等の作物が対象作物となる場合、出穂日に替えて、到達時期を用いて収穫状況が予測される。到達時期とは、対象作物が出穂以外の予め定められた生長段階に到達した時期である。また、予測処理に用いられる気象情報又は数式等は、対象地域又は対象作物に応じて定められる。 In the following, a prediction process for predicting the yield and quality of rice, which is a target crop, will be described (FIG. 1B). This process is executed by the server 10 that operates according to the prediction program. In the prediction processing, various areas are targeted, rice of other varieties is targeted as crops, and various crops (for example, crops having spikes such as wheat, or vegetables such as tomatoes or green soybeans). It is possible to predict the harvest situation by targeting various varieties. When crops such as vegetables without spikes are the target crops, the harvest status is predicted using the arrival time instead of the heading date. The arrival time is the time when the target crop reaches a predetermined growth stage other than heading. In addition, the meteorological information or mathematical formulas used for the prediction processing are determined according to the target area or the target crop.

S100では、サーバ10は、出穂日の予測に用いる入力変数を算出する。サーバ10は、予め定められた1又は複数の田植基準日における、予め定められた1又は複数の種類の気象情報を、説明変数とする。入力変数は、複数の異なる説明変数を用いて生成される。なお、該複数の説明変数の数を、Nとする。また、複数の説明変数の各々を、E0[i]とも記載する。iは1〜Nのうちのいずれかの整数である。 In S100, the server 10 calculates an input variable used for predicting the heading date. The server 10 uses the predetermined one or more types of weather information on the predetermined one or more rice planting reference dates as explanatory variables. Input variables are generated with a number of different explanatory variables. The number of the plurality of explanatory variables is N. In addition, each of the plurality of explanatory variables is also described as E0 [i]. i is an integer of any of 1 to N.

稲の出穂に相対的に大きな影響を与える気象情報が、説明変数となる。一例として、説明変数とは、田植日から3日目の日の降水量、田植日から22日目の日の最低気温、田植日から26日目の日の平均気温、及び、田植日から26日目の日の最高気温等であっても良い。説明変数に係る田植基準日及び説明変数となる気象情報の種類は、記録期間における気象履歴及び栽培履歴に基づき定められる。説明変数の数、説明変数に係る田植基準日、又は、説明変数となる気象情報の種類等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 Meteorological information that has a relatively large effect on the heading of rice is an explanatory variable. As an example, the explanatory variables are the precipitation on the 3rd day from the rice planting date, the minimum temperature on the 22nd day from the rice planting date, the average temperature on the 26th day from the rice planting date, and 26 from the rice planting date. It may be the maximum temperature of the day. The rice planting reference date related to the explanatory variable and the type of meteorological information used as the explanatory variable are determined based on the meteorological history and cultivation history during the recording period. The number of explanatory variables, the rice planting reference date related to the explanatory variables, the type of weather information used as the explanatory variables, etc. differ depending on the type of the target crop, the variety of the target crop, the target area, and the like.

サーバ10は、記憶部11に記憶されている前年気象履歴及び対象年気象履歴のうちの双方又は一方から、説明変数を抽出する。この時、サーバ10は、予測を行っている日よりも前の田植基準日の気象情報である説明変数については、対象年気象履歴から抽出しても良い。一方、予測を行っている日以後の田植基準日の気象情報である説明変数については、前年気象履歴から抽出しても良い。つまり、前年における該田植基準日に相当する日における該説明変数と同一種類の気象情報を、前年気象履歴から抽出しても良い。なお、田植基準日が対象年のX月Y日の場合であれば、他の年における該田植基準日に相当する日は、該他の年のX月Y日となる。すなわち、サーバ10は、原則として対象年気象履歴から説明変数を抽出し、抽出できない場合には、前年気象履歴から説明変数を抽出しても良い。無論、サーバ10は、前年気象履歴から全ての説明変数を抽出しても良い。 The server 10 extracts explanatory variables from both or one of the previous year's meteorological history and the target year's meteorological history stored in the storage unit 11. At this time, the server 10 may extract the explanatory variables that are the meteorological information of the rice planting reference date before the forecasting date from the meteorological history of the target year. On the other hand, the explanatory variables that are the meteorological information of the rice planting reference date after the forecasting date may be extracted from the previous year's meteorological history. That is, the same type of meteorological information as the explanatory variable on the day corresponding to the rice planting reference date in the previous year may be extracted from the previous year's meteorological history. If the rice planting reference date is the X month and Y day of the target year, the day corresponding to the rice planting reference date in another year is the X month and Y day of the other year. That is, in principle, the server 10 may extract explanatory variables from the target year's meteorological history, and if it cannot be extracted, it may extract explanatory variables from the previous year's meteorological history. Of course, the server 10 may extract all the explanatory variables from the previous year's weather history.

次に、サーバ10は、式(1)により説明変数E0[i]を正規化し、scale(E0[i])を得る。 Next, the server 10 normalizes the explanatory variable E0 [i] by the equation (1) to obtain a scale (E0 [i]).

Figure 0006849958
Figure 0006849958

ここで、記録期間に含まれる各年におけるE0[i]に係る田植基準日に相当する日を、対応日とする。また、対応日におけるE0[i]と同一種類の気象情報を、対応説明変数とする。ave(E0[i])は、記録期間に含まれる各年の対応説明変数の平均値であり、σ(E0[i])は、各年の対応説明変数の標準偏差である。 Here, the date corresponding to the rice planting reference date pertaining to E0 [i] in each year included in the recording period is defined as the corresponding date. Further, the same type of weather information as E0 [i] on the corresponding date is used as the corresponding explanatory variable. ave (E0 [i]) is the average value of the corresponding explanatory variables for each year included in the recording period, and σ (E0 [i]) is the standard deviation of the corresponding explanatory variables for each year.

そして、サーバ10は、N個のscale(E0[i])の主成分分析を行う。主成分分析により生成された合成変数が、入力変数X0となる。具体的には、式(2)により、M個の入力変数X0[1]〜[M]を生成する。なお、Mは整数であり、1<M<Nとなる。 Then, the server 10 analyzes the principal components of N scales (E0 [i]). The synthetic variable generated by the principal component analysis is the input variable X0. Specifically, M input variables X0 [1] to [M] are generated by the equation (2). Note that M is an integer and 1 <M <N.

Figure 0006849958
Figure 0006849958

なお、式(2)における主成分負荷量R011〜R0NMは、記録期間に含まれる各年における各説明変数の対応説明変数に基づき算出される。また、入力変数の数、又は、主成分負荷量の値等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 The principal component load amounts R0 11 to R0 NM in the equation (2) are calculated based on the corresponding explanatory variables of each explanatory variable in each year included in the recording period. Further, the number of input variables, the value of the principal component load, and the like differ depending on the type of the target crop, the variety of the target crop, the target area, and the like.

S105では、サーバ10は、多層パーセプトロンであるニューラルネットワークにより、出穂日を予測する。なお、図2Aの200は、このようなニューラルネットワークの一例である。ニューラルネットワークの入力層は、M個の入力変数に対応している。ニューラルネットワークは、1層の中間層を有し、中間層にはM個のニューロンが設けられている。また、ニューラルネットワークの出力層は、1つの出力変数を出力する。また、ニューラルネットワークでは、式(3)に示すシグモイド関数が出力関数g(x)として用いられる。 In S105, the server 10 predicts the heading date by a neural network that is a multi-layer perceptron. Note that 200 in FIG. 2A is an example of such a neural network. The input layer of the neural network corresponds to M input variables. The neural network has one intermediate layer, and M neurons are provided in the intermediate layer. Further, the output layer of the neural network outputs one output variable. Further, in the neural network, the sigmoid function shown in the equation (3) is used as the output function g (x).

Figure 0006849958
Figure 0006849958

ニューラルネットワークには、M個の入力変数X0[1]〜[M]が入力される。以後、j番目の入力変数を、X0[j]とする。なお、jは1〜Mのうちのいずれかの整数である。 M input variables X0 [1] to [M] are input to the neural network. Hereinafter, the j-th input variable is set to X0 [j]. In addition, j is an integer of any one of 1 to M.

また、中間層におけるM個のニューロンのうちのk番目のニューロンからの出力を、Y0[k]とする。なお、kは1〜Mのうちのいずれかの整数である。また、中間層におけるM個のニューロンのうちのk番目のニューロンにおけるX0[j]に対応する重みパラメータを、W0kjとする。また、該k番目のニューロンにおけるバイアスパラメータを、θ0とする。Y0[k]は、式(4)により算出される。 Further, the output from the k-th neuron among the M neurons in the intermediate layer is defined as Y0 [k]. In addition, k is an integer of any one of 1 to M. Further, the weight parameter corresponding to X0 [j] in the k-th neuron among the M neurons in the intermediate layer is W0 kj . Further, the bias parameter in the k-th neuron is set to θ0 k . Y0 [k] is calculated by the formula (4).

Figure 0006849958
Figure 0006849958

また、ニューラルネットワークの出力層は、出穂日を示すZ0を出力する。出力層のニューロンにおける重みパラメータをW0とし、バイアスパラメータをΘ0とする。Z0は、式(5)により算出される。 Further, the output layer of the neural network outputs Z0 indicating the heading date. Let the weight parameter in the neurons of the output layer be W0 k and the bias parameter be Θ0. Z0 is calculated by the formula (5).

Figure 0006849958
Figure 0006849958

Z0は、対象年の出穂日の予測値を正規化した値となっている。このため、式(6)によりZ0の形式を変更する。これにより、出穂日の予測値であるZ0´が得られる。なお、出穂日は、例えば、田植日等からの経過日数により表されても良い。 Z0 is a value obtained by normalizing the predicted value of the heading date of the target year. Therefore, the format of Z0 is changed by the equation (6). As a result, Z0', which is a predicted value of the heading date, can be obtained. The heading date may be represented by, for example, the number of days elapsed from the rice planting date or the like.

Figure 0006849958
Figure 0006849958

なお、σ(z0)は、記録期間における各年の出穂日の標準偏差を示す。また、ave(z0)は、記録期間における各年の出穂日の平均値を示す。また、ニューラルネットワークにおける重みパラメータW0kj,W0及びバイアスパラメータθ0,Θ0は、記録期間における栽培履歴が示す過去の出穂日を教師データとする機械学習により定められる。重みパラメータW0kj,W0、又は、バイアスパラメータθ0,Θ0は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる値となる。また、対象種類、対象作物の品種、又は、対象地域等に応じて、ニューラルネットワークの構成を変えても良い。 Note that σ (z0) indicates the standard deviation of the heading date of each year during the recording period. Further, ave (z0) indicates the average value of the heading date of each year in the recording period. Further, the weight parameters W0 kj and W0 k and the bias parameters θ0 k and Θ0 in the neural network are determined by machine learning using the past heading date indicated by the cultivation history in the recording period as teacher data. The weight parameters W0 kj , W0 k , or the bias parameters θ0 k , Θ0 have different values depending on the type of the target crop, the variety of the target crop, the target area, and the like. Further, the configuration of the neural network may be changed according to the target type, the variety of the target crop, the target area, and the like.

続いて、S110では、サーバ10は、収穫量の予測に用いる入力変数を算出する。サーバ10は、予め定められた1又は複数の出穂基準日における予め定められた1又は複数の種類の気象情報を、説明変数とする。なお、出穂基準日は、S105にて予測された出穂日に基づき定められる。入力変数は、複数の異なる説明変数を用いて生成される。また、該複数の説明変数の数を、Oとする。また、複数の説明変数の各々を、E1[i]とも記載する。iは1〜Oのうちのいずれかの整数である。 Subsequently, in S110, the server 10 calculates an input variable used for predicting the yield. The server 10 uses a predetermined one or more types of weather information on a predetermined heading reference date as an explanatory variable. The heading reference date is determined based on the heading date predicted in S105. Input variables are generated with a number of different explanatory variables. Further, let O be the number of the plurality of explanatory variables. In addition, each of the plurality of explanatory variables is also described as E1 [i]. i is an integer of any of 1 to O.

稲の収穫量に相対的に大きな影響を与える気象情報が、説明変数となる。一例として、説明変数とは、出穂日の13日前の日の降水量、出穂日から10日目の日の気温日格差、及び、出穂日から10日目の日の日照時間等であっても良い。説明変数に係る出穂基準日及び説明変数となる気象情報の種類は、記録期間における気象履歴及び栽培履歴に基づき定められる。説明変数の数、説明変数に係る出穂基準日、又は、説明変数となる気象情報の種類等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 Meteorological information, which has a relatively large effect on rice yield, is an explanatory variable. As an example, the explanatory variables may be precipitation on the day 13 days before the heading date, temperature-day difference on the 10th day from the heading date, and sunshine hours on the 10th day from the heading date. good. The heading reference date related to the explanatory variable and the type of meteorological information used as the explanatory variable are determined based on the meteorological history and cultivation history during the recording period. The number of explanatory variables, the heading reference date related to the explanatory variables, the type of weather information used as the explanatory variables, etc. differ depending on the type of the target crop, the variety of the target crop, the target area, and the like.

サーバ10は、記憶部11に記憶されている前年気象履歴及び対象年気象履歴のうちの双方又は一方から、説明変数を抽出する。この時、サーバ10は、S100と同様にして、原則として対象年気象履歴から説明変数を抽出し、抽出できない場合には、前年気象履歴から説明変数を抽出しても良い。無論、サーバ10は、前年気象履歴から全ての説明変数を抽出しても良い。 The server 10 extracts explanatory variables from both or one of the previous year's meteorological history and the target year's meteorological history stored in the storage unit 11. At this time, the server 10 may extract the explanatory variables from the target year's meteorological history in principle in the same manner as in S100, and if it cannot be extracted, the explanatory variables may be extracted from the previous year's meteorological history. Of course, the server 10 may extract all the explanatory variables from the previous year's weather history.

次に、サーバ10は、式(7)により説明変数E1[i]を正規化し、scale(E1[i])を得る。 Next, the server 10 normalizes the explanatory variable E1 [i] by the equation (7) to obtain a scale (E1 [i]).

Figure 0006849958
Figure 0006849958

ここで、記録期間に含まれる各年におけるE1[i]に係る出穂基準日に相当する日を、対応日とする。また、対応日におけるE1[i]と同一種類の気象情報を、対応説明変数とする。ave(E1[i])は、記録期間に含まれる各年の対応説明変数の平均値であり、σ(E1[i])は、各年の対応説明変数の標準偏差である。 Here, the date corresponding to the heading reference date related to E1 [i] in each year included in the recording period is defined as the corresponding date. Further, the same type of weather information as E1 [i] on the corresponding date is used as the corresponding explanatory variable. ave (E1 [i]) is the average value of the corresponding explanatory variables for each year included in the recording period, and σ (E1 [i]) is the standard deviation of the corresponding explanatory variables for each year.

そして、サーバ10は、O個のscale(E1[i])の主成分分析を行う。主成分分析により生成された合成変数が、入力変数X1となる。具体的には、式(8)により、P個の入力変数X1[1]〜[P]を生成する。なお、Pは整数であり、1<P<Oとなる。 Then, the server 10 analyzes the principal components of O scales (E1 [i]). The synthetic variable generated by the principal component analysis is the input variable X1. Specifically, P input variables X1 [1] to [P] are generated by the equation (8). Note that P is an integer and 1 <P <O.

Figure 0006849958
Figure 0006849958

なお、式(8)における主成分負荷量R111〜R1OPは、記録期間に含まれる各年における各説明変数の対応説明変数に基づき算出される。また、入力変数の数、又は、主成分負荷量の値等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 The principal component load amounts R1 11 to R1 OP in the equation (8) are calculated based on the corresponding explanatory variables of each explanatory variable in each year included in the recording period. Further, the number of input variables, the value of the principal component load, and the like differ depending on the type of the target crop, the variety of the target crop, the target area, and the like.

S115では、サーバ10は、多層パーセプトロンであるニューラルネットワークにより、収穫量を予測する。該ニューラルネットワークは、出穂日の予測に用いられたニューラルネットワークと同様に構成されているが、入力層が対応する入力変数の数、及び、中間層に含まれるニューロンの数が相違する。すなわち、収穫量の予測に用いられるニューラルネットワークの入力層は、P個の入力変数に対応している。また、該ニューラルネットワークは、中間層にはP個のニューロンが設けられている。また、該ニューラルネットワークでは、式(3)に示すシグモイド関数が出力関数g(x)として用いられる。 In S115, the server 10 predicts the yield by means of a neural network that is a multi-layer perceptron. The neural network is configured in the same manner as the neural network used for predicting the heading date, but the number of input variables corresponding to the input layer and the number of neurons contained in the intermediate layer are different. That is, the input layer of the neural network used for predicting the yield corresponds to P input variables. Further, in the neural network, P neurons are provided in the intermediate layer. Further, in the neural network, the sigmoid function shown in the equation (3) is used as the output function g (x).

ニューラルネットワークには、P個の入力変数X1[1]〜[P]が入力される。以後、j番目の入力変数を、X1[j]とする。なお、jは1〜Pのうちのいずれかの整数である。 P input variables X1 [1] to [P] are input to the neural network. Hereinafter, the j-th input variable will be referred to as X1 [j]. In addition, j is an integer of any one of 1 to P.

また、中間層におけるP個のニューロンのうちのk番目のニューロンからの出力を、Y1[k]とする。なお、kは1〜Pのうちのいずれかの整数である。また、中間層におけるP個のニューロンのうちのk番目のニューロンにおけるX1[j]に対応する重みパラメータを、W1kjとする。また、該k番目のニューロンにおけるバイアスパラメータを、θ1とする。Y1[k]は、式(9)により算出される。 Further, the output from the k-th neuron among the P neurons in the intermediate layer is defined as Y1 [k]. In addition, k is an integer of any one of 1 to P. Further, the weight parameter corresponding to X1 [j] in the k-th neuron among the P neurons in the intermediate layer is W1 kj . Further, the bias parameter in the k-th neuron is set to θ1 k . Y1 [k] is calculated by the equation (9).

Figure 0006849958
Figure 0006849958

また、ニューラルネットワークの出力層は、収穫量を示すZ1を出力する。出力層のニューロンにおける重みパラメータをW1とし、バイアスパラメータをΘ1とする。Z1は、式(10)により算出される。 Further, the output layer of the neural network outputs Z1 indicating the yield. Let the weight parameter in the neurons of the output layer be W1 k and the bias parameter be Θ1. Z1 is calculated by the formula (10).

Figure 0006849958
Figure 0006849958

Z1は、対象年の収穫量の予測値を正規化した値となっている。このため、式(11)によりZ1の形式を変更する。これにより、収穫量の予測値であるZ1´が得られる。 Z1 is a value obtained by normalizing the predicted value of the yield in the target year. Therefore, the format of Z1 is changed by the equation (11). As a result, Z1', which is a predicted value of the yield, can be obtained.

Figure 0006849958
Figure 0006849958

なお、σ(z1)は、記録期間における各年の収穫量の標準偏差を示す。また、ave(z1)は、記録期間における各年の収穫量の平均値を示す。また、ニューラルネットワークにおける重みパラメータW1kj,W1及びバイアスパラメータθ1,Θ1は、記録期間における栽培履歴が示す過去の収穫量を教師データとする機械学習により定められる。重みパラメータW1kj,W1及びバイアスパラメータθ1,Θ1は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる値となる。また、対象作物の種類、対象作物の品種、及び、対象地域等に応じて、予測に用いるニューラルネットワークの構成を変えても良い。 Note that σ (z1) indicates the standard deviation of the yield for each year during the recording period. In addition, ave (z1) indicates the average value of the yield of each year in the recording period. Further, the weight parameters W1 kj and W1 k and the bias parameters θ1 k and Θ1 in the neural network are determined by machine learning using the past yield indicated by the cultivation history in the recording period as teacher data. The weight parameters W1 kj and W1 k and the bias parameters θ1 k and Θ1 have different values depending on the type of the target crop, the variety of the target crop, the target area, and the like. Further, the configuration of the neural network used for prediction may be changed according to the type of the target crop, the variety of the target crop, the target area, and the like.

続いて、S120では、サーバ10は、品質の予測に用いる入力変数を算出する。サーバ10は、S110と同様、予め定められた1又は複数の出穂基準日における予め定められた1又は複数の種類の気象情報を、説明変数とする。入力変数は、複数の異なる説明変数を用いて生成される。また、該複数の説明変数の数を、Qとする。また、複数の説明変数の各々を、E2[i]とも記載する。iは1〜Qのうちのいずれかの整数である。 Subsequently, in S120, the server 10 calculates an input variable used for quality prediction. Similar to S110, the server 10 uses a predetermined one or a plurality of types of weather information on a predetermined heading reference date as an explanatory variable. Input variables are generated with a number of different explanatory variables. Further, let Q be the number of the plurality of explanatory variables. In addition, each of the plurality of explanatory variables is also described as E2 [i]. i is an integer of any of 1 to Q.

稲の品質に相対的に大きな影響を与える気象情報が、説明変数となる。一例として、説明変数とは、出穂日の10日前の日の平均気温、出穂日から8日目の日の最高気温、及び、出穂日から19日目の日の日照時間等であっても良い。説明変数に係る出穂基準日及び説明変数となる気象情報の種類は、記録期間における気象履歴及び栽培履歴に基づき定められる。説明変数の数、説明変数に係る出穂基準日、又は、説明変数となる気象情報の種類等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 Meteorological information, which has a relatively large effect on rice quality, is an explanatory variable. As an example, the explanatory variables may be the average temperature on the day 10 days before the heading date, the maximum temperature on the 8th day from the heading date, the sunshine hours on the 19th day from the heading date, and the like. .. The heading reference date related to the explanatory variable and the type of meteorological information used as the explanatory variable are determined based on the meteorological history and cultivation history during the recording period. The number of explanatory variables, the heading reference date related to the explanatory variables, the type of weather information used as the explanatory variables, etc. differ depending on the type of the target crop, the variety of the target crop, the target area, and the like.

サーバ10は、記憶部11に記憶されている前年気象履歴及び対象年気象履歴のうちの双方又は一方から、説明変数を抽出する。この時、サーバ10は、S100と同様にして、原則として対象年気象履歴から説明変数を抽出し、抽出できない場合には、前年気象履歴から説明変数を抽出しても良い。無論、サーバ10は、前年気象履歴から全ての説明変数を抽出しても良い。 The server 10 extracts explanatory variables from both or one of the previous year's meteorological history and the target year's meteorological history stored in the storage unit 11. At this time, the server 10 may extract the explanatory variables from the target year's meteorological history in principle in the same manner as in S100, and if it cannot be extracted, the explanatory variables may be extracted from the previous year's meteorological history. Of course, the server 10 may extract all the explanatory variables from the previous year's weather history.

そして、サーバ10は、Q個のE2[i]の主成分分析を行う。主成分分析により生成された合成変数が、入力変数X2となる。具体的には、式(12)により、R個の入力変数X2[1]〜[R]を生成する。なお、Rは整数であり、1<R<Qとなる。 Then, the server 10 analyzes the principal components of Q E2 [i]. The synthetic variable generated by the principal component analysis is the input variable X2. Specifically, R input variables X2 [1] to [R] are generated by the equation (12). R is an integer, and 1 <R <Q.

Figure 0006849958
Figure 0006849958

ここで、記録期間に含まれる各年におけるE2[i]に係る出穂基準日に相当する日を、対応日とする。また、対応日におけるE2[i]と同一種類の気象情報を、対応説明変数とする。式(12)における主成分負荷量R211〜R2QRは、記録期間に含まれる各年における各説明変数の対応説明変数に基づき算出される。なお、入力変数の数、又は、主成分負荷量の値等は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる。 Here, the date corresponding to the heading reference date related to E2 [i] in each year included in the recording period is defined as the corresponding date. Further, the same type of weather information as E2 [i] on the corresponding date is used as the corresponding explanatory variable. The principal component load amounts R2 11 to R2 QR in the formula (12) are calculated based on the corresponding explanatory variables of each explanatory variable in each year included in the recording period. The number of input variables, the value of the principal component load, etc. differ depending on the type of target crop, the variety of target crop, the target area, and the like.

S125では、サーバ10は、X2[j]を説明変数とした重回帰分析により品質を予測する。なお、jは1〜Rのうちのいずれかの整数である。具体的には、式(13)により、品質の予測値であるY2が算出される。 In S125, the server 10 predicts the quality by multiple regression analysis using X2 [j] as an explanatory variable. In addition, j is an integer of any one of 1 to R. Specifically, Y2, which is a predicted quality value, is calculated by the equation (13).

Figure 0006849958
Figure 0006849958

そして、サーバ10は、本処理を終了する。
なお、式(13)における偏回帰係数a〜aは、以下のようにして算出される。すなわち、サーバ10は、記録期間に含まれる各年について、S120と同様にして入力変数X2を生成する。つまり、サーバ10は、S120と同様にして、各年の気象履歴から該年に対応する説明変数E2[i]を抽出する。この時、各年のE2[i]に係る出穂基準日は、該年の出穂日に基づき定められる。そして、サーバ10は、S120と同様、式(12)により、説明変数E2[i]に基づきX2[j]を生成する。さらに、サーバ10は、各年のX2[j]と各年の品質とに基づき、偏回帰係数a〜aを算出する。なお、偏回帰係数a〜aは、対象作物の種類、対象作物の品種、又は、対象地域等に応じて異なる値となる。
Then, the server 10 ends this process.
The partial regression coefficients a 1 to a R in the equation (13) are calculated as follows. That is, the server 10 generates the input variable X2 in the same manner as in S120 for each year included in the recording period. That is, the server 10 extracts the explanatory variable E2 [i] corresponding to the year from the weather history of each year in the same manner as S120. At this time, the heading reference date for E2 [i] of each year is determined based on the heading date of the year. Then, the server 10 generates X2 [j] based on the explanatory variable E2 [i] according to the equation (12) as in S120. Further, the server 10 calculates the partial regression coefficients a 0 to a R based on X2 [j] of each year and the quality of each year. The partial regression coefficients a 1 to a R have different values depending on the type of the target crop, the variety of the target crop, the target area, and the like.

[3.効果]
以上詳述した上記実施形態によれば、以下の効果が得られる。
(1)作物が予め定められた生長段階に達する時期、若しくは、該時期の前又は後の時期の気象が、作物の生長に相対的に大きな影響を与える場合がある。特に、気象の変動が激しい場合、このような時期の気象が例年と大きく異なる可能性があり、その結果、収穫状況が例年と比べて大きく変動する恐れがある。
[3. effect]
According to the above-described embodiment described in detail above, the following effects can be obtained.
(1) The weather at the time when the crop reaches a predetermined growth stage or before or after the stage may have a relatively large influence on the growth of the crop. In particular, when the weather fluctuates drastically, the weather at such times may differ significantly from the usual year, and as a result, the harvest situation may fluctuate significantly compared to the usual year.

これに対し、本実施形態の予測処理によれば、対象作物である稲の出穂日が予測される。そして、予測された出穂日を基準として定められた出穂基準日等の気象情報に基づき、稲の収穫状況が予測される。このため、稲の生長に相対的に大きな影響を与える時期の気象に基づき、収穫状況を予測することが可能となる。したがって、稲の収穫状況を精度良く予測することが可能となる。 On the other hand, according to the prediction process of the present embodiment, the heading date of the target crop, rice, is predicted. Then, the rice harvest situation is predicted based on the weather information such as the heading reference date determined based on the predicted heading date. Therefore, it is possible to predict the harvest situation based on the weather at the time when the rice growth is relatively significantly affected. Therefore, it is possible to accurately predict the rice harvest situation.

(2)また、作物の栽培期間における特定の日の気象が、作物の生長に相対的に大きな影響を与える場合がある。これに対し、本実施形態の予測処理によれば、田植日を基準に定められた日である田植基準日等の気象情報に基づき出穂日が予測される。そして、予測された出穂日を基準として定められた日である出穂基準日等の気象情報に基づき、収穫状況が予測される。このため、対象作物の生長への影響が相対的に大きい特定の日の気象に基づき、出穂日及び収穫状況を予測することが可能となる。したがって、稲の収穫状況を精度良く予測することが可能となる。 (2) In addition, the weather on a specific day during the cultivation period of the crop may have a relatively large effect on the growth of the crop. On the other hand, according to the prediction processing of the present embodiment, the heading date is predicted based on the weather information such as the rice planting reference date, which is the date determined based on the rice planting date. Then, the harvest situation is predicted based on the weather information such as the heading reference date, which is the date determined based on the predicted heading date. Therefore, it is possible to predict the heading date and the harvest situation based on the weather on a specific day, which has a relatively large effect on the growth of the target crop. Therefore, it is possible to accurately predict the rice harvest situation.

(3)また、作物の穂又は花が予め定められた生長段階に到達した時期、若しくは、該時期の前又は後の時期の気象が、作物の生長に相対的に大きな影響を与える場合がある。これに対し、本実施形態の予測処理によれば、予測された出穂日を基準に定められた時期の気象情報により、稲の収穫状況が予測される。このため、収穫状況を精度良く予測することが可能となる。 (3) In addition, the weather when the spikes or flowers of the crop reach a predetermined growth stage, or before or after that stage may have a relatively large effect on the growth of the crop. .. On the other hand, according to the prediction process of the present embodiment, the rice harvest situation is predicted based on the weather information at the time determined based on the predicted heading date. Therefore, it is possible to accurately predict the harvest situation.

[4.他の実施形態]
以上、本開示を実施するための形態について説明したが、本開示は上記実施形態に限定されることなく、種々変形して実施することができる。
[4. Other embodiments]
Although the embodiment for carrying out the present disclosure has been described above, the present disclosure is not limited to the above embodiment, and can be implemented in various modifications.

(1)上記実施形態の予測処理では、対象年又は前年の気象情報に基づき出穂日が予測される。そして、予測された出穂日に基づき定められた出穂基準日、又は、前年における出穂基準日に相当する日の気象情報に基づき、収穫状況が予測される。以後、予測に用いる気象情報が測定された時期を、測定時期と記載する。しかしながら、予測プログラムは、予測処理において、対象作物が出穂以外の他の生長段階に到達する到達時期を予測しても良い。そして、予測プログラムは、本実施形態と同様にして、予測した到達時期を基準として測定時期を定め、測定時期における予め定められた種類の気象情報に基づき収穫状況を予測しても良い。 (1) In the prediction process of the above embodiment, the heading date is predicted based on the weather information of the target year or the previous year. Then, the harvest situation is predicted based on the weather information of the heading reference date determined based on the predicted heading date or the day corresponding to the heading reference date in the previous year. Hereinafter, the time when the meteorological information used for the prediction is measured will be referred to as the measurement time. However, the prediction program may predict the arrival time when the target crop reaches a growth stage other than heading in the prediction process. Then, the prediction program may determine the measurement time based on the predicted arrival time and predict the harvest situation based on the predetermined type of meteorological information at the measurement time in the same manner as in the present embodiment.

具体的には、到達時期とは、例えば、対象作物の穂又は花が、予め定められた生長段階に到達した時期であっても良い。具体的には、対象作物が、幼穂分化、花芽分化、開花、又は、落花等の生長段階に達した時期であっても良い。なお、幼穂分化とは、例えば、対象作物に穂が形成され始めた生長段階であり、花芽分化とは、例えば、対象作物に花芽が形成され始めた生長段階であっても良い。この他にも、到達時期に係る生長段階とは、例えば、対象作物の受粉であっても良いし、対象作物の茎、枝、又は葉等の大きさ、長さ、又は数等が、予め定められた基準に達した段階であっても良い。 Specifically, the arrival time may be, for example, the time when the spikes or flowers of the target crop reach a predetermined growth stage. Specifically, the target crop may reach a growth stage such as ear differentiation, flower bud differentiation, flowering, or falling flowers. The ear differentiation may be, for example, a growth stage at which ears have begun to be formed on the target crop, and the flower bud differentiation may be, for example, a growth stage at which flower buds have begun to be formed on the target crop. In addition to this, the growth stage related to the arrival time may be, for example, pollination of the target crop, and the size, length, number, etc. of the stem, branch, leaf, etc. of the target crop may be determined in advance. It may be at the stage when the set standard is reached.

また、到達時期とは、対象農地で栽培されている対象作物における上記生長段階に到達した対象作物が占める割合が、予め定められた閾値(例えば、40〜50%)に達した日であっても良い。また、到達時期又は測定時期とは、いずれかの日に限らず、例えば、複数の日を跨いだ期間、又は、いずれかの日における時間帯等であっても良い。到達時期及び測定時期は、対象作物の種類、対象作物の品種、又は、対象地域等に応じて定めるのが好適である。 The arrival time is the day when the ratio of the target crops cultivated on the target farmland to the target crops that have reached the growth stage reaches a predetermined threshold value (for example, 40 to 50%). Is also good. Further, the arrival time or the measurement time is not limited to any day, and may be, for example, a period spanning a plurality of days, a time zone on any day, or the like. It is preferable to determine the arrival time and the measurement time according to the type of the target crop, the variety of the target crop, the target area, and the like.

(2)上記実施形態の予測処理では、田植日に基づき定められた田植基準日の気象情報に基づき出穂日が予測される。しかしながら、予測プログラムは、予測処理において、例えば、対象年又は前年における田植日に関わらず定められた月日の気象情報に基づき、出穂日を予測しても良い。 (2) In the prediction process of the above embodiment, the heading date is predicted based on the weather information of the rice planting reference date determined based on the rice planting date. However, the forecasting program may predict the heading date in the forecasting process, for example, based on the meteorological information of the specified date regardless of the rice planting date in the target year or the previous year.

また、例えば、対象作物の種を対象農地に植えて栽培する場合等においては、予測処理において、田植基準日に替えて開始基準日が用いられても良い。開始基準日とは、対象農地に対象作物の種を植えた日である開始日を基準として定められた日である。すなわち、予測プログラムは、予測処理において、田植基準日に替えて開始基準日を用い、本実施形態と同様にして、対象作物が予め定められた生長段階に到達する到達時期を予測しても良い。そして、予測プログラムは、到達時期を基準として定められた測定時期の気象情報に基づき、対象作物の収穫状況を予測しても良い。 Further, for example, when the seeds of the target crop are planted and cultivated on the target agricultural land, the start reference date may be used instead of the rice planting reference date in the prediction processing. The start reference date is a date set based on the start date, which is the date when the seeds of the target crop are planted on the target farmland. That is, the prediction program may use the start reference date instead of the rice planting reference date in the prediction process to predict the arrival time when the target crop reaches a predetermined growth stage in the same manner as in the present embodiment. .. Then, the prediction program may predict the harvest status of the target crop based on the meteorological information of the measurement time determined based on the arrival time.

(3)上記実施形態の予測処理では、対象作物である稲の出穂日及び収穫量の予測にはニューラルネットワークが用いられ、品質の予測には重回帰分析が用いられる。しかしながら、重回帰分析により出穂日及び収穫量が予測されても良いし、ニューラルネットワークにより品質が予測されても良い。また、他の方法により、出穂日、収穫量、及び、品質が予測されても良い。 (3) In the prediction process of the above embodiment, a neural network is used for predicting the heading date and yield of rice as a target crop, and multiple regression analysis is used for predicting quality. However, the heading date and yield may be predicted by multiple regression analysis, and the quality may be predicted by a neural network. In addition, the heading date, yield, and quality may be predicted by other methods.

(4)上記実施形態の予測処理では、対象年の前年の気象情報である前年気象履歴を用いて収穫状況が予測される。しかし、予測処理では、前年気象履歴に限らず、対象年よりも前の任意の1又は複数の年の気象履歴を用いて、収穫状況が予測されても良い。つまり、予測処理では、前年気象履歴に替えて、前年よりも前のいずれかの年の気象履歴から、同様にして説明変数が抽出されても良い。この他にも、対象年よりも前の任意の複数年の気象履歴を用いて、収穫状況が予測されても良い。つまり、該複数の年の気象履歴に基づき、平均気象履歴が生成されても良い。平均気象履歴は、これらの年の各日における各種類の気象情報の平均値を示す。そして、予測処理では、前年気象履歴に替えて、平均気象履歴の中から同様にして説明変数が抽出されても良い。 (4) In the prediction process of the above embodiment, the harvest situation is predicted using the previous year's weather history, which is the weather information of the previous year of the target year. However, in the prediction process, the harvest situation may be predicted using not only the previous year's meteorological history but also the meteorological history of any one or more years before the target year. That is, in the prediction process, the explanatory variables may be extracted in the same manner from the meteorological history of any year before the previous year instead of the previous year's meteorological history. In addition to this, the harvest situation may be predicted using the meteorological history of any plurality of years prior to the target year. That is, the average weather history may be generated based on the weather history of the plurality of years. The average weather history shows the average value of each type of weather information on each day of these years. Then, in the prediction process, the explanatory variables may be extracted from the average weather history in the same manner instead of the previous year's weather history.

(5)上記実施形態における1つの構成要素が有する複数の機能を、複数の構成要素によって実現したり、1つの構成要素が有する1つの機能を、複数の構成要素によって実現したりしてもよい。また、複数の構成要素が有する複数の機能を、1つの構成要素によって実現したり、複数の構成要素によって実現される1つの機能を、1つの構成要素によって実現したりしてもよい。また、上記実施形態の構成の一部を省略してもよい。また、上記実施形態の構成の少なくとも一部を、他の上記実施形態の構成に対して付加又は置換してもよい。なお、特許請求の範囲に記載した文言のみによって特定される技術思想に含まれるあらゆる態様が本開示の実施形態である。 (5) A plurality of functions possessed by one component in the above embodiment may be realized by a plurality of components, or one function possessed by one component may be realized by a plurality of components. .. Further, a plurality of functions possessed by the plurality of components may be realized by one component, or one function realized by the plurality of components may be realized by one component. Further, a part of the configuration of the above embodiment may be omitted. In addition, at least a part of the configuration of the above embodiment may be added or replaced with the configuration of the other above embodiment. It should be noted that all aspects included in the technical idea specified only by the wording described in the claims are embodiments of the present disclosure.

(6)上述した予測システム1の他、当該予測システム1に相当する装置、予測システム1における予測プログラムを記録した半導体メモリ等の非遷移的実態的記録媒体、予測処理に相当する方法など、種々の形態で本開示を実現することもできる。 (6) In addition to the prediction system 1 described above, there are various devices such as a device corresponding to the prediction system 1, a non-transitional actual recording medium such as a semiconductor memory in which the prediction program in the prediction system 1 is recorded, and a method corresponding to the prediction processing. The present disclosure can also be realized in the form of.

[特許請求の範囲との対応]
上記実施形態の説明で用いた用語と、特許請求の範囲の記載に用いた用語との対応を示す。
[Correspondence with claims]
The correspondence between the terms used in the description of the above-described embodiment and the terms used in describing the scope of claims is shown.

予測処理により実現される方法が予測方法の一例に、サーバ10が予測装置の一例に相当する。
また、対象作物である稲の出穂が生長段階の一例に、出穂日が到達時期の一例に相当する。また、出穂基準日、及び、前年における出穂基準日に相当する日が、測定時期の一例に相当する。また、田植日が開始日の一例に相当する。
The method realized by the prediction process corresponds to an example of the prediction method, and the server 10 corresponds to an example of the prediction device.
In addition, the heading of rice, which is the target crop, corresponds to an example of the growth stage, and the heading date corresponds to an example of the arrival time. In addition, the heading reference date and the day corresponding to the heading reference date in the previous year correspond to an example of the measurement time. In addition, the rice planting date corresponds to an example of the start date.

また、予測処理のS100,S105が到達時期予測部の一例に、S110〜S125が収穫状況予測部の一例に相当する。 Further, S100 and S105 of the prediction process correspond to an example of the arrival time prediction unit, and S110 to S125 correspond to an example of the harvest status prediction unit.

1…予測システム、10…サーバ、11…記憶部、12…出穂日予測部、13…収穫量予測部、14…品質予測部、20…端末、30…測定装置。 1 ... Prediction system, 10 ... server, 11 ... storage unit, 12 ... heading date prediction unit, 13 ... yield prediction unit, 14 ... quality prediction unit, 20 ... terminal, 30 ... measuring device.

Claims (8)

予め定められた地域における予め定められた作物の収穫状況を予測する予測方法であって、
前記作物が予め定められた生長段階に到達する到達時期を、該到達時期を含む年である対象年と、前記対象年よりも前の任意の年とのうちの双方又は一方における前記地域の気象情報に基づき予測し、
前記対象年における前記収穫状況を、多層パーセプトロンであるニューラルネットワークにより予測し、該ニューラルネットワークに入力される第1入力変数は、前記対象年と前記任意の年とのうちの双方又は一方における前記到達時期を基準として定められた特定の日である第1基準日における前記気象情報であって、前記収穫状況への影響が相対的に大きい前記気象情報に基づき生成されること、
を特徴とする予測方法。
It is a predictive method for predicting the harvest situation of a predetermined crop in a predetermined area.
The time of arrival of the crop at a predetermined growth stage is the weather in the area in both or one of the target year, which is the year including the arrival time, and any year prior to the target year. Informed predictions
The harvest situation in the target year is predicted by a neural network which is a multi-layer perceptron, and the first input variable input to the neural network is the arrival in both or one of the target year and the arbitrary year. The meteorological information on the first reference date, which is a specific day determined based on the time , and is generated based on the meteorological information having a relatively large influence on the harvest situation.
A prediction method characterized by.
請求項1に記載の予測方法において、In the prediction method according to claim 1,
前記第1入力変数は、前記第1基準日における前記気象情報に対し主成分分析を行うことで、生成されること、The first input variable is generated by performing principal component analysis on the meteorological information on the first reference date.
を特徴とする予測方法。A prediction method characterized by.
請求項1又は請求項2に記載の予測方法において、
前記対象年における前記作物を植えた日を、開始日とし、
前記到達時期を、多層パーセプトロンであるニューラルネットワークにより予測し、該ニューラルネットワークに入力される第2入力変数は、前記対象年と前記任意の年とのうちの双方又は一方における前記開始日を基準に定められた特定の日である第2基準日における前記気象情報であって、前記作物が前記生長段階に到達することへの影響が相対的に大きい前記気象情報に基づき生成されること、
を特徴とする予測方法。
In the prediction method according to claim 1 or 2.
The day when the crop was planted in the target year was set as the start date.
The arrival time is predicted by a neural network that is a multi-layer perceptron, and the second input variable input to the neural network is based on the start date in both or one of the target year and the arbitrary year. The meteorological information on the second reference date, which is a specified specific day, and is generated based on the meteorological information having a relatively large influence on the crop reaching the growth stage.
A prediction method characterized by.
請求項3に記載の予測方法において、In the prediction method according to claim 3,
前記第2入力変数は、前記第2基準日における前記気象情報に対し主成分分析を行うことで、生成されること、The second input variable is generated by performing principal component analysis on the meteorological information on the second reference date.
を特徴とする予測方法。A prediction method characterized by.
請求項1から請求項4のうちのいずれか1項に記載の予測方法において、
前記到達時期とは、前記作物の穂又は花が予め定められた生長段階に到達した時期であること、
を特徴とする予測方法。
In the prediction method according to any one of claims 1 to 4,
The arrival time is the time when the spikes or flowers of the crop reach a predetermined growth stage.
A prediction method characterized by.
請求項1から請求項のうちのいずれか1項に記載の予測方法において、
前記作物とは、稲であること、
を特徴とする予測方法。
In the prediction method according to any one of claims 1 to 5,
The crop is rice
A prediction method characterized by.
予め定められた地域における予め定められた作物の収穫状況を予測する予測プログラムであって、
前記作物が予め定められた生長段階に到達する到達時期を、該到達時期を含む年である対象年と、前記対象年よりも前の任意の年とのうちの双方又は一方における前記地域の気象情報に基づき予測する到達時期予測部と、
前記対象年における前記収穫状況を、多層パーセプトロンであるニューラルネットワークにより予測する部位であって、該ニューラルネットワークに入力される第1入力変数は、前記対象年と前記任意の年とのうちの双方又は一方における前記到達時期を基準として定められた特定の日である第1基準日における前記気象情報であって、前記収穫状況への影響が相対的に大きい前記気象情報に基づき生成される部位である収穫状況予測部として、
コンピュータを動作させることを特徴とする予測プログラム。
A forecasting program that predicts the harvest status of predetermined crops in a predetermined area.
The time of arrival of the crop at a predetermined growth stage is the weather in the area in both or one of the target year, which is the year including the arrival time, and any year prior to the target year. The arrival time prediction unit that predicts based on information,
The part that predicts the harvest situation in the target year by a neural network that is a multi-layer perceptron, and the first input variable input to the neural network is both the target year and the arbitrary year, or On the other hand, it is the meteorological information on the first reference date, which is a specific day determined based on the arrival time , and is a part generated based on the meteorological information having a relatively large influence on the harvest situation. As a harvest status prediction department
A prediction program characterized by running a computer.
予め定められた地域における予め定められた作物の収穫状況を予測する予測装置であって、
前記作物が予め定められた生長段階に到達する到達時期を、該到達時期を含む年である対象年と、前記対象年よりも前の任意の年とのうちの双方又は一方における前記地域の気象情報に基づき予測する到達時期予測部と、
前記対象年における前記収穫状況を、多層パーセプトロンであるニューラルネットワークにより予測する部位であって、該ニューラルネットワークに入力される第1入力変数は、前記対象年と前記任意の年とのうちの双方又は一方における前記到達時期を基準として定められた特定の日である第1基準日における前記気象情報であって、前記収穫状況への影響が相対的に大きい前記気象情報に基づき生成される部位である収穫状況予測部と、
を備える予測装置。
It is a predictor that predicts the harvest situation of a predetermined crop in a predetermined area.
The time of arrival of the crop at a predetermined growth stage is the weather in the area in both or one of the target year, which is the year including the arrival time, and any year prior to the target year. The arrival time prediction unit that predicts based on information,
The part that predicts the harvest situation in the target year by a neural network that is a multi-layer perceptron, and the first input variable input to the neural network is both the target year and the arbitrary year, or On the other hand, it is the meteorological information on the first reference date, which is a specific day determined based on the arrival time , and is a part generated based on the meteorological information having a relatively large influence on the harvest situation. Harvest status prediction department and
Predictor device.
JP2016055340A 2016-03-18 2016-03-18 Prediction method, prediction program, and prediction device Active JP6849958B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016055340A JP6849958B2 (en) 2016-03-18 2016-03-18 Prediction method, prediction program, and prediction device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016055340A JP6849958B2 (en) 2016-03-18 2016-03-18 Prediction method, prediction program, and prediction device

Publications (2)

Publication Number Publication Date
JP2017163956A JP2017163956A (en) 2017-09-21
JP6849958B2 true JP6849958B2 (en) 2021-03-31

Family

ID=59909403

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016055340A Active JP6849958B2 (en) 2016-03-18 2016-03-18 Prediction method, prediction program, and prediction device

Country Status (1)

Country Link
JP (1) JP6849958B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6850715B2 (en) * 2017-11-16 2021-03-31 ヤンマーパワーテクノロジー株式会社 Carry-in order determination device
JP2019170359A (en) * 2018-03-28 2019-10-10 国立大学法人名古屋大学 Plant cultivation result prediction system
KR102184813B1 (en) * 2018-12-13 2020-11-30 대한민국 Method for predicting of rice crops using similar weather information
JP7075127B2 (en) * 2019-02-22 2022-05-25 株式会社ナイルワークス Field analysis method, field analysis program, field analysis device, drone system and drone
JP2020171217A (en) * 2019-04-09 2020-10-22 株式会社セラク Method, apparatus and program for predicting yield
JP7229864B2 (en) * 2019-06-28 2023-02-28 株式会社日立製作所 REMOTE SENSING IMAGE ACQUISITION TIME DETERMINATION SYSTEM AND CROPT GROWTH ANALYSIS METHOD
JP7163881B2 (en) * 2019-08-02 2022-11-01 トヨタ自動車株式会社 Crop Characteristics Prediction System, Crop Characteristics Prediction Method, and Crop Characteristics Prediction Program

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006212008A (en) * 2005-02-01 2006-08-17 Nagano Prefecture Crop cultivation supporting device
JP2009106261A (en) * 2007-10-30 2009-05-21 Asuzac Inc Crop cultivation-supporting device
CN101480143B (en) * 2009-01-21 2010-11-03 华中科技大学 Method for predicating single yield of crops in irrigated area
JP6147579B2 (en) * 2013-06-18 2017-06-14 株式会社日立製作所 Yield prediction system and yield prediction device
JP6394445B2 (en) * 2014-06-16 2018-09-26 Tdk株式会社 Crop harvest time estimation device and crop harvest time estimation program

Also Published As

Publication number Publication date
JP2017163956A (en) 2017-09-21

Similar Documents

Publication Publication Date Title
JP6849958B2 (en) Prediction method, prediction program, and prediction device
US10942113B2 (en) Methods, systems, and components thereof relating to using multi-spectral imaging for improved cultivation of cannabis and other crops
KR101936317B1 (en) Method for smart farming
Sankar et al. Non-destructive estimation of leaf area of durian (Durio zibethinus)–An artificial neural network approach
Bushong et al. Development of an in-season estimate of yield potential utilizing optical crop sensors and soil moisture data for winter wheat
JP2000300077A (en) Method for determining fertilizing quantity for grain crop, method for estimating quality and yield of grain and apparatus for providing production information on grain
Sagar et al. Agriculture data analytics in crop yield estimation: a critical review
Çakır et al. Yield prediction of wheat in south-east region of Turkey by using artificial neural networks
Lambert et al. Behavior patterns related to the agricultural practices in the production of Persian lime (Citrus latifolia tanaka) in the seasonal orchard
Banerjee et al. Crop Status Index as an indicator of wheat crop growth condition under abiotic stress situations
Molinari et al. Towards an integrated weed management decision support system: A simulation model for weed-crop competition and control
Yadav et al. Calibration and validation of PNUTGRO (DSSAT v4. 5) model for yield and yield attributing characters of kharif groundnut cultivars in middle Gujarat region
Mishra et al. Modeling and forecasting of sugarcane production in South Asian countries
Saxena et al. Weather based model development for outbreak of mustard aphid (L. erysimi, kalt) using artificial neural network
Higgins et al. Forecasting maturity of green peas: An application of neural networks
Banakara et al. Evaluation of weather parameter-based pre-harvest yield forecast models for wheat crop: a case study in Saurashtra region of Gujarat
Garde et al. Pre harvest forecasting of kharif rice yield using weather parameters for strategic decision making in agriculture
Pasley et al. APSIM next generation mungbean model: A tool for advancing mungbean production
Banakara et al. Pre-harvest forecast of kharif rice yield using PCA and MLR technique in Navsari district of Gujarat
JP2022136058A (en) Method of generating prediction model for predicting crop production performance, generation apparatus, and generation program
Salam et al. Climate recommender system for wheat cultivation in North Egyptian Sinai Peninsula
JP7140086B2 (en) Fruit Vegetable Plant and Fruit Tree Cultivation Management Apparatus, Learning Device, Fruit Vegetable Plant and Fruit Tree Cultivation Management Method, Learning Model Generation Method, Fruit Vegetable Plant and Fruit Tree Cultivation Management Program, and Learning Model Generation Program
Agarawal Forecasting techniques in crops
Sapak et al. POMICS: A simulation disease model for timing fungicide applications in management of powdery mildew of cucurbits
Brotons-Martínez et al. The Use of Fuzzy Decoupled Net Present Value in pepper production

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190315

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20191218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200114

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200316

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200623

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210126

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210218

R150 Certificate of patent or registration of utility model

Ref document number: 6849958

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250