JP7385931B2 - Information processing device, information processing method, and program - Google Patents

Information processing device, information processing method, and program Download PDF

Info

Publication number
JP7385931B2
JP7385931B2 JP2021039461A JP2021039461A JP7385931B2 JP 7385931 B2 JP7385931 B2 JP 7385931B2 JP 2021039461 A JP2021039461 A JP 2021039461A JP 2021039461 A JP2021039461 A JP 2021039461A JP 7385931 B2 JP7385931 B2 JP 7385931B2
Authority
JP
Japan
Prior art keywords
growth
crop
information processing
trait index
processing device
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
JP2021039461A
Other languages
Japanese (ja)
Other versions
JP2022139192A (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.)
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
Priority to JP2021039461A priority Critical patent/JP7385931B2/en
Publication of JP2022139192A publication Critical patent/JP2022139192A/en
Application granted granted Critical
Publication of JP7385931B2 publication Critical patent/JP7385931B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、情報処理装置、情報処理方法、およびプログラムに関する。 The present invention relates to an information processing device, an information processing method, and a program.

従来、栽培されている農作物に関する情報を取得する情報検出手段と、農作物の生育状況を予測するための農作物登録情報を記憶する記憶手段と、農作物に関する情報に基づき農作物登録情報を参照して、農作物の生育状況を予測する生育予測手段と、生育予測手段で予測された生育状況を表示する表示手段と、を有する圃場管理システムが知られている(特許文献1)。 Conventionally, an information detecting means for acquiring information on cultivated crops, a storage means for storing crop registration information for predicting the growth status of the crops, and a storage means for storing crop registration information based on the information on the crops, A field management system is known that includes a growth prediction means for predicting the growth situation of a field, and a display means for displaying the growth situation predicted by the growth prediction means (Patent Document 1).

特開2016-49102号公報Japanese Patent Application Publication No. 2016-49102

従来の技術では、様々な条件や外乱要素などの影響要因を考慮して個体生育モデルを作成しておらず、そのため作物の状態を精度よく予測することができない場合があった。 In conventional techniques, individual growth models are not created in consideration of influencing factors such as various conditions and disturbance elements, and as a result, it is sometimes not possible to accurately predict the condition of crops.

本発明は、このような事情を考慮してなされたものであり、作物の状態を精度よく予測することができる情報処理装置、情報処理方法、およびプログラムを提供することを目的の一つとする。 The present invention has been made in consideration of such circumstances, and one of its objects is to provide an information processing device, an information processing method, and a program that can accurately predict the condition of crops.

本発明の第1の態様である情報処理装置は、作物に与えられる条件および外乱要素を含む影響要因の情報を取得する取得部と、前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出する進行速度導出部と、基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出する状態値導出部と、を備えるものである。 The information processing device according to the first aspect of the present invention includes an acquisition unit that acquires information on influencing factors including conditions given to crops and disturbance elements, and a step to the next growth stage based on the influencing factors and the number of growing days. A progress speed derivation unit that derives the progress speed toward the transition for each growth trait index of the crop at least, and a vector in which the integral values of the progress speed from the reference time to the predicted time are arranged, and a state value deriving unit that derives a predicted state vector including a state value for each growth trait index of the crop by multiplying by a connection matrix obtained by subtracting the input arc from the output arc in a transition event.

本発明の第2の態様である情報処理装置は、第1の態様において、前記進行速度導出部は、同じ生育形質指標に関する前記進行速度を、生育段階ごとに異なる傾向で導出し、前記接続マトリクスは、生育段階を問わず前記積分値の用法を定義するものである。 In the information processing device according to a second aspect of the present invention, in the first aspect, the progression rate derivation unit derives the progression rate regarding the same growth trait index with a different tendency for each growth stage, and defines the usage of the integral value regardless of the growth stage.

本発明の第3の態様である情報処理装置は、第1または第2の態様において、前記取得部は、環境測定値に基づいて前記影響要因の情報を生成する処理を行うものである。 In the information processing device according to a third aspect of the present invention, in the first or second aspect, the acquisition unit performs a process of generating information on the influencing factor based on an environmental measurement value.

本発明の第4の態様である情報処理装置は、第1から第3の態様において、前記取得部は、環境測定値に基づいて前記影響要因の情報を生成する処理の手法を変えながら予測状態ベクトルを導出し、前記予測状態ベクトルの実際の作物への当てはまり度合いに基づいて前記影響要因の情報を生成する処理の手法を評価する評価部を更に備えるものである。 In the information processing device according to a fourth aspect of the present invention, in the first to third aspects, the acquisition unit changes the predicted state while changing a processing method for generating information on the influencing factor based on the environmental measurement value. The present invention further includes an evaluation unit that derives a vector and evaluates a processing method for generating information on the influencing factor based on the degree of applicability of the predicted state vector to the actual crop.

本発明の第5の態様である情報処理装置は、少なくとも前記影響要因、前記遷移イベント、および前記生育形質指標の関係をグラフ化した画像を表示装置に表示させる表示制御部を更に備えるものである。 The information processing device according to a fifth aspect of the present invention further includes a display control unit that causes a display device to display an image that is a graph of the relationship among at least the influencing factor, the transition event, and the growth trait index. .

本発明の第6の態様である情報処理装置は、生育形質指標の変化をグラフ化した画像を表示装置に表示させる表示制御部を更に備えるものである。 The information processing device according to the sixth aspect of the present invention further includes a display control unit that causes a display device to display an image in which changes in the growth trait index are graphed.

本発明の第7の態様である情報処理方法は、コンピュータが、作物に与えられる条件および外乱要素を含む影響要因の情報を取得し、前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出し、基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出するものである。 In the information processing method according to the seventh aspect of the present invention, a computer acquires information on influencing factors including conditions given to crops and disturbance elements, and proceeds to the next growth stage based on the influencing factors and the number of growing days. The rate of progress toward the transition of the crop is derived for at least each growth trait index of the crop, and the output arc at the transition event of the growth stage is created in a vector in which the integral values of the rate of progress from the reference time to the predicted time are arranged. By multiplying by a connection matrix obtained by subtracting the input arc from , a predicted state vector including state values for each growth trait index of the crop is derived.

本発明の第8の態様であるプログラムは、コンピュータに、作物に与えられる条件および外乱要素を含む影響要因の情報を取得させ、前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出させ、基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出させるものである。 A program according to an eighth aspect of the present invention causes a computer to acquire information on influencing factors including conditions given to crops and disturbance elements, and determines the transition to the next growth stage based on the influencing factors and the number of growing days. The rate of progress towards the crop is derived at least for each growth trait index of the crop, and inputted from the output arc at the transition event of the growth stage into a vector in which integral values of the rate of progress from the reference time to the predicted time are arranged. By multiplying the connection matrix by subtracting the arc, a predicted state vector including the state value for each growth trait index of the crop is derived.

上記各態様によれば、作物の状態を精度よく予測することができる。 According to each of the above aspects, the condition of crops can be predicted with high accuracy.

実施例に係る情報処理装置100の構成と処理対象の一例を示す図である。FIG. 1 is a diagram illustrating an example of a configuration of an information processing device 100 and processing targets according to an embodiment. 取得部110の機能について説明するための図である。FIG. 2 is a diagram for explaining the functions of an acquisition unit 110. FIG. 小麦の生育プロセスを示す図である。It is a diagram showing the growth process of wheat. 進行速度導出部120および状態値導出部130の機能について説明するための図である。FIG. 3 is a diagram for explaining the functions of a traveling speed deriving section 120 and a state value deriving section 130. t=0すなわち生育開始日の状態をモデルによって計算した結果を示す図である。It is a figure which shows the result of calculating the state of t=0, ie, the growth start day, by a model. t=5.5のときの状態をモデルによって計算した結果を示す図である。FIG. 7 is a diagram showing the results of calculating the state when t=5.5 using a model. t=6のときの状態をモデルによって計算した結果を示す図である。FIG. 6 is a diagram showing the results of calculating the state when t=6 using a model. t=30のときの状態をモデルによって計算した結果を示す図である。FIG. 7 is a diagram showing the result of calculating the state when t=30 using a model. t=30のときの状態をモデルによって計算した結果を全体的に示す図である。FIG. 3 is a diagram showing the overall result of calculating the state when t=30 using a model. t=33のときに鳥獣による被害が発生した状態をモデルによって計算した結果を全体的に示す図である。FIG. 7 is a diagram showing the overall results of calculation using a model of a state in which damage caused by birds and beasts occurs at t=33. 鳥獣による被害が発生したことを反映させて、作物の成熟(t=166)まで計算を行った結果をグラフで表した図である。It is a diagram showing in a graph the results of calculations performed until the maturity of crops (t=166), reflecting the occurrence of damage caused by birds and animals.

以下、図面を参照し、本発明の情報処理装置、情報処理方法、およびプログラムの実施形態について説明する。 DESCRIPTION OF THE PREFERRED EMBODIMENTS Embodiments of an information processing apparatus, an information processing method, and a program according to the present invention will be described below with reference to the drawings.

図1は、実施例に係る情報処理装置100の構成と処理対象の一例を示す図である。情報処理装置100は、例えば、自然環境の耕作地における対象圃場の個体株ごとに、作物の生育形質指標ごとの状態値を導出する装置である。情報処理装置100は、例えば、取得部110と、進行速度導出部120と、状態値導出部130と、評価部140と、表示制御部150とを備える。これらの構成要素は、例えば、CPU(Central Processing Unit)などのハードウェアプロセッサがプログラム(ソフトウェア)を実行することにより実現される。これらの構成要素のうち一部または全部は、LSI(Large Scale Integration)やASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)、GPU(Graphics Processing Unit)などのハードウェア(回路部;circuitryを含む)によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。プログラムは、予めHDD(Hard Disk Drive)やフラッシュメモリなどの記憶装置(非一過性の記憶媒体を備える記憶装置)に格納されていてもよいし、DVDやCD-ROMなどの着脱可能な記憶媒体(非一過性の記憶媒体)に格納されており、記憶媒体がドライブ装置に装着されることで記憶装置にインストールされてもよい。 FIG. 1 is a diagram illustrating an example of the configuration and processing targets of an information processing apparatus 100 according to an embodiment. The information processing device 100 is, for example, a device that derives state values for each crop growth trait index for each individual plant in a target field on cultivated land in a natural environment. The information processing device 100 includes, for example, an acquisition section 110, a progress speed derivation section 120, a state value derivation section 130, an evaluation section 140, and a display control section 150. These components are realized by, for example, a hardware processor such as a CPU (Central Processing Unit) executing a program (software). Some or all of these components are hardware (circuit parts) such as LSI (Large Scale Integration), ASIC (Application Specific Integrated Circuit), FPGA (Field-Programmable Gate Array), and GPU (Graphics Processing Unit). (including circuitry), or may be realized by collaboration between software and hardware. The program may be stored in advance in a storage device (a storage device with a non-transitory storage medium) such as an HDD (Hard Disk Drive) or flash memory, or may be stored in a removable storage device such as a DVD or CD-ROM. It may be stored in a medium (non-transitory storage medium), and installed in the storage device by loading the storage medium into a drive device.

以下、情報処理装置100の各部の機能について説明する。図2は、取得部110の機能について説明するための図である。 The functions of each part of the information processing device 100 will be explained below. FIG. 2 is a diagram for explaining the functions of the acquisition unit 110.

取得部110は、作物に与えられる影響要因の情報を取得する。影響要因は、後述する進行速度の計算過程に与えられる入力情報である。影響要因には、環境測定値に基づいて取得部110が生成するものと、イベントの発生状況を示すものとがある。前者について取得部110は、環境測定値から環境要因の情報を生成する。環境測定値とは、例えば、時系列の気温データから求められる積算温度量、時系列の日照データから求められる日照量、時系列の給水量データから求められる水量、時系列の酸素供給量データから求められる酸素量、時系列の施肥データから求められる積算温度量である。 The acquisition unit 110 acquires information on influencing factors given to crops. The influencing factor is input information given to the process of calculating the traveling speed, which will be described later. The influencing factors include those generated by the acquisition unit 110 based on environmental measurement values and those that indicate the occurrence of an event. Regarding the former, the acquisition unit 110 generates information on environmental factors from the environmental measurement values. Environmental measurements include, for example, the cumulative temperature obtained from time-series temperature data, the amount of sunlight obtained from time-series sunlight data, the amount of water obtained from time-series water supply data, and the amount of oxygen supply obtained from time-series data. This is the amount of oxygen required and the cumulative temperature amount determined from time-series fertilization data.

取得部110は、積算温度量に基づいて影響要因の一つである温度条件を生成する。取得部110は、日照量に基づいて影響要因の一つである日照条件を生成する。取得部110は、水量に基づいて影響要因の一つである水条件を生成する。取得部110は、酸素量に基づいて影響要因の一つである酸素条件を生成する。取得部110は、肥料量に基づいて影響要因の一つである肥料条件を生成する。温度条件、日照条件、水条件、酸素条件、肥料条件のそれぞれは、進行速度の計算過程に与えることが可能な任意の形態の情報である。取得部110は、これらの影響要因の生成手法を、例えばパラメータの変更等によって変更することが可能である。 The acquisition unit 110 generates a temperature condition, which is one of the influencing factors, based on the cumulative temperature amount. The acquisition unit 110 generates a sunshine condition, which is one of the influencing factors, based on the amount of sunlight. The acquisition unit 110 generates water conditions, which are one of the influencing factors, based on the amount of water. The acquisition unit 110 generates an oxygen condition, which is one of the influencing factors, based on the amount of oxygen. The acquisition unit 110 generates fertilizer conditions, which are one of the influencing factors, based on the amount of fertilizer. Each of the temperature conditions, sunlight conditions, water conditions, oxygen conditions, and fertilizer conditions is any form of information that can be provided to the progress rate calculation process. The acquisition unit 110 can change the generation method of these influencing factors by, for example, changing parameters.

また、イベントの発生状況を示す影響要因には、病虫害発生状況、鳥獣による被害状況、人的関与(防除)条件、雑草発生状況、その他がある。これらの影響要因は、例えば発生有りを示す1、発生無しを示すゼロなどで表される情報である。 Further, influencing factors indicating the occurrence status of an event include the occurrence status of pests and diseases, the status of damage caused by birds and animals, human involvement (control) conditions, the status of weed occurrence, and others. These influencing factors are information represented by, for example, 1 indicating occurrence, zero indicating non-occurrence, etc.

このようにして取得された環境要因は、進行速度導出部120の処理に使用される。なお図中、丸で表されているものは作物自らの生育条件や外的要因による条件または状態を示し、丸の中の黒点は一つの条件または状態を表す。環境要因は、遷移イベントに影響するものである。遷移イベントとは、例えば、発芽、分げつ、出穂、登熟といったものを含む。遷移イベント間の状態を、生育段階と称する。図2では、生育段階の例として発芽後、分げつ後(出穂前)を記載している。遷移イベントが到来すると、生育段階が次の生育段階に遷移する。 The environmental factors obtained in this way are used for processing by the progress speed deriving unit 120. In the figure, circles represent conditions or states caused by the crop's own growth conditions or external factors, and black dots within the circles represent one condition or state. Environmental factors are those that influence transition events. Transition events include, for example, germination, tillering, heading, and ripening. The states between transition events are referred to as growth stages. In FIG. 2, post-germination and post-tillering (before heading) are shown as examples of growth stages. When a transition event arrives, a growth stage transitions to the next growth stage.

以下、作物の一例として小麦を例にとって説明する。図3は、小麦の生育プロセスを示す図である。小麦は、発芽、分げつ、出穂、登熟という生育段階を経て成熟に至る。各生育段階において、最高分げつ、穂ばらみ、開花、黄熱といった状態が出現する。図4以降では、小麦に対してモデルが適用されることを前提に、各部の機能について説明する。なお本発明は小麦に限らず、あらゆる作物あるいは植物に適用可能である。 Below, wheat will be explained as an example of a crop. FIG. 3 is a diagram showing the growth process of wheat. Wheat reaches maturity through the growth stages of germination, tillering, heading, and ripening. At each growth stage, conditions such as peak tillering, ear opening, flowering, and yellow fever occur. From FIG. 4 onwards, the functions of each part will be explained on the premise that the model is applied to wheat. Note that the present invention is applicable not only to wheat but also to any crops or plants.

図4は、進行速度導出部120および状態値導出部130の機能について説明するための図である。進行速度導出部120は、影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも作物の生育形質指標ごとに導出する。生育形質指標とは、例えば、草丈、葉数、茎数、穂数などである。生育形質指標は、図4において二重丸で表されている。生育形質指標には、穂数のように特定の生育段階に至った後に出現するものも存在する。 FIG. 4 is a diagram for explaining the functions of the advancing speed deriving section 120 and the state value deriving section 130. The progress rate deriving unit 120 derives, based on the influencing factors and the number of growing days, the rate of progress toward transition to the next growth stage for at least each growth trait index of the crop. Growth trait indicators include, for example, plant height, number of leaves, number of stems, and number of panicles. Growth trait indicators are represented by double circles in FIG. Some growth trait indicators, such as the number of panicles, appear after reaching a specific growth stage.

状態値導出部130は、基準時点(例えば種まき時点)から予測時点までの進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスAを乗算することで、作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出する。 The state value derivation unit 130 multiplies a vector in which the integral values of the progress speed from the reference time point (for example, the time of sowing seeds) to the predicted time point are arranged by a connection matrix A obtained by subtracting the input arc from the output arc in the transition event of the growth stage. By doing so, a predicted state vector including state values for each growth trait index of the crop is derived.

以下、進行速度導出部120および状態値導出部130の機能について、より簡易な例と計算結果を示しながら説明する。以下の説明において、eは影響因子、tは生育日数、jは生育段階、kは生育形質指標の識別子を表す。図5~8は、進行速度導出部120および状態値導出部130の処理について例示した図である。これらの図の例では、T1とT2の二つの遷移イベントがある。また、説明を簡略化するために生育段階を分げつ後(出穂前)まで、生育形質指標を草丈に絞って説明するが、他の部分についても同様に展開できる。図5は、t=0すなわち生育開始日の状態をモデルによって計算した結果を示す図である。 Hereinafter, the functions of the advancing speed deriving section 120 and the state value deriving section 130 will be explained while showing a simpler example and calculation results. In the following description, e represents an influencing factor, t represents the number of growing days, j represents a growth stage, and k represents an identifier of a growth trait index. 5 to 8 are diagrams illustrating the processing of the advancing speed deriving section 120 and the state value deriving section 130. In the example of these figures, there are two transition events, T1 and T2. Furthermore, in order to simplify the explanation, the growth stage will be explained up to after tillering (before heading), and the growth trait index will be focused on plant height, but the same can be applied to other parts. FIG. 5 is a diagram showing the results of calculation using a model for the state of t=0, that is, the growth start date.

遷移イベントの切り替わり日数は予め定められており、T1が6日後、T2が120日後に発生すると定義されている。切り替わり日数は可変であってもよいし、当該作物に関して経験から求められた適合値であってもよい。この場合、進行速度導出部120が導出する進行速度は、以下の式で表される。式中、E(e,t,j,k)は、影響因子e、生育日数t、生育段階j、生育形質指標kに基づいて進行速度を導出する関数またはアルゴリズムである。E(e,t,j,k)は、例えば、当該作物に関して経験から求められたものである。 The switching days of the transition event are predetermined, and it is defined that T1 will occur after 6 days and T2 will occur after 120 days. The number of switching days may be variable or may be an appropriate value determined from experience for the crop in question. In this case, the traveling speed derived by the traveling speed deriving unit 120 is expressed by the following formula. In the formula, E (e, t, j, k) is a function or algorithm that derives the progress rate based on the influencing factor e, the number of growing days t, the growing stage j, and the growth trait index k. E(e,t,j,k) is, for example, determined from experience regarding the crop in question.

Figure 0007385931000001
Figure 0007385931000001

図5~8に示すモデルに関して状態値導出部130が用いる接続マトリクス(接続行列)Aは、以下の式で表される。接続マトリクスAは、生育段階の遷移イベントにおける出力アークOから入力アークIを差し引いたものである。入力アークとは、遷移イベントよりも前の状態であって、当該遷移イベントまで作用する状態について1、そうでない状態についてゼロが定義されるベクトルを、遷移イベントごとに並べたマトリクスである。出力アークとは、遷移イベントよりも後の状態であって、当該遷移イベントの直後から作用する状態について1、そうでない状態についてゼロが定義されるベクトルを、遷移イベントごとに並べたマトリクスである。 A connection matrix (connection matrix) A used by the state value deriving unit 130 regarding the models shown in FIGS. 5 to 8 is expressed by the following equation. The connectivity matrix A is the output arc O minus the input arc I at the growth stage transition event. The input arc is a matrix in which vectors are arranged for each transition event, with vectors defined as 1 for states that are prior to a transition event and that affect the transition event, and zero for states that do not. An output arc is a matrix in which vectors are arranged for each transition event, with vectors defined as 1 for states that act immediately after the transition event and zero for other states.

Figure 0007385931000002
Figure 0007385931000002

進行速度導出部120および状態値導出部130が行う処理を一般式で表すと、以下の式のようになる。式中、Sは初期状態の予測状態ベクトルであり、Sは生育日数tにおける予測状態ベクトルであり、uは時間増分である。SおよびSの要素には生育形質指標が少なくとも含まれ、例えば、{P(t),P(t),P(t),P(t),P2-1(t),P3-1(t)}で表される。このように、進行速度導出部120は、同じ生育形質指標に関する進行速度を、生育段階ごとに異なる傾向で導出し、接続マトリクスAは、生育段階を問わず積分値の用法(取捨選択、およびプラスで作用させるかマイナスで作用させるか)を定義するものである。 If the processing performed by the advancing speed deriving unit 120 and the state value deriving unit 130 is expressed in a general formula, it will be as shown in the following formula. where S 0 is the predicted state vector at the initial state, S t is the predicted state vector at growing days t, and u is the time increment. The elements of S 0 and S t include at least a growth trait index, for example, {P 1 (t), P 2 (t), P 3 (t), P e (t), P 2-1 (t ),P 3-1 (t)} Represented by T. In this way, the progress rate derivation unit 120 derives the progress rate regarding the same growth trait index with different tendencies for each growth stage, and the connection matrix A is configured to calculate the usage of integral values (selecting and adding This defines whether the effect is to be applied in a positive or negative manner.

Figure 0007385931000003
Figure 0007385931000003

図6は、t=5.5のときの状態をモデルによって計算した結果を示す図である。t=5.5のとき、P2-1すなわち生育段階が「発芽後」である作物の生育形質指標「草丈」として1.0が算出されている。この場合の計算式は以下の式で表される。 FIG. 6 is a diagram showing the results of calculating the state when t=5.5 using a model. When t=5.5, 1.0 is calculated as P2-1 , that is, the growth trait index "plant height" of the crop whose growth stage is "post-germination." The calculation formula in this case is expressed by the following formula.

Figure 0007385931000004
Figure 0007385931000004

図7は、t=6のときの状態をモデルによって計算した結果を示す図である。t=6のとき、P2-1すなわち生育段階が「発芽後」である作物の生育形質指標「草丈」として1.2が算出されている。この場合の計算式は以下の式で表される。 FIG. 7 is a diagram showing the result of calculating the state when t=6 using a model. When t=6, P2-1 , that is, 1.2 is calculated as the growth trait index "plant height" of the crop whose growth stage is "post-germination." The calculation formula in this case is expressed by the following formula.

Figure 0007385931000005
Figure 0007385931000005

図8は、t=30のときの状態をモデルによって計算した結果を示す図である。t=30のとき、P3-1すなわち生育段階が「分げつ後(出穂前)」である作物の生育形質指標「草丈」として11.5が算出されている。この場合の計算式は以下の式で表される。t=6で遷移の条件を満たしたことで遷移イベント(分げつ)が発生したため、t=6以降の分の生育形質指標は、v(p3-1)に基づいて決定される。 FIG. 8 is a diagram showing the result of calculating the state when t=30 using a model. When t=30, 11.5 is calculated as the growth trait index "plant height" of P3-1 , that is, the crop whose growth stage is "post-tillering (before heading)". The calculation formula in this case is expressed by the following formula. Since the transition event (tillering) occurred when the transition conditions were satisfied at t=6, the growth trait index for t=6 and thereafter is determined based on v 2 (p 3-1 ).

Figure 0007385931000006
Figure 0007385931000006

上記のように計算を行うことで、任意の生育日数について各生育形質指標を求めることができる。 By performing calculations as described above, each growth trait index can be determined for any number of growth days.

再度、生育段階の全体に関して説明する。図9は、t=30のときの状態をモデルによって計算した結果を全体的に示す図である。図示するように、草丈だけでなく葉数や茎数も同様の計算で算出されている。 The entire growth stage will be explained again. FIG. 9 is a diagram showing the overall result of calculating the state when t=30 using a model. As shown in the figure, not only the plant height but also the number of leaves and stems are calculated using the same method.

ここで、病虫害発生状況、鳥獣による被害状況、人的関与(防除)条件、雑草発生状況などの影響要因を反映させる手法について説明する。例えば、t=33のときに鳥獣による被害が発生したと仮定すると、取得部110は、例えば、所定の生育形質指標に関して、実際の計測値や予め定められた被害後の生育形質値または被害度合いから計算するという形で影響要因を定義する。図10は、t=33のときに鳥獣による被害が発生した状態をモデルによって計算した結果を全体的に示す図である。図示するように、草丈と葉数に関してt=30のときの状態から値が減少しているのが分かる。 Here, we will explain a method that reflects influencing factors such as the occurrence of pests and diseases, damage caused by birds and animals, human involvement (control) conditions, and weed occurrence. For example, assuming that damage caused by birds and animals occurs at time t=33, the acquisition unit 110 acquires, for example, actual measured values, predetermined post-damage growth trait values, or damage levels with respect to a predetermined growth trait index. Influence factors are defined by calculating them from FIG. 10 is a diagram generally showing the results of calculation using a model of the state in which damage caused by birds and beasts occurs at t=33. As shown in the figure, it can be seen that the values of plant height and number of leaves have decreased from the state when t=30.

図11は、鳥獣による被害が発生したことを反映させて、作物の成熟(t=166)まで計算を行った結果をグラフで表した図である。図示するように、t=33における被害によって減少した草丈と葉数がその後に回復し、作物が成長しているのが分かる。 FIG. 11 is a graph showing the results of calculations performed until crop maturity (t=166), reflecting the occurrence of damage caused by birds and animals. As shown in the figure, it can be seen that the plant height and number of leaves that were reduced due to damage at t=33 have since recovered, and the crop is growing.

図1に戻り、取得部110は、環境測定値に基づいて影響要因の情報を生成する処理の手法を変えながら予測状態ベクトルを導出する。評価部140は、予測状態ベクトルの実際の作物への当てはまり度合いに基づいて、影響要因の情報を生成する処理の手法を評価する。そして、評価部140は、対象圃場において当てはまりの良い(相似性の高い)所定数のモデル(影響要因の仕方)を抽出して群落の将来状態を予測する。当てはまりの良い所定数のモデルとは、最大評価値から一定の決定枠範囲内(例えば最大評価値の90%~100%の範囲内)である。評価部140は、例えば、遺伝的アルゴリズム、焼き鈍し法のような発見的手法やモンテカルロ法などを用いて評価を行ってよい。また、評価部140は、状態値導出部130の導出結果に基づいて、作物の生育量や活性化指数などによる作物の収穫量を増加するための施肥量を定めても良いし、病虫害を防ぐための予察に利用しても良いし、播種時期毎の最適な植え付け本数を定めるのに利用しても良いし、収穫物の品質・等級を優先するための肥培管理を求めても良い。 Returning to FIG. 1, the acquisition unit 110 derives a predicted state vector while changing the processing method for generating information on influencing factors based on environmental measurement values. The evaluation unit 140 evaluates the processing method for generating influencing factor information based on the degree of applicability of the predicted state vector to the actual crop. Then, the evaluation unit 140 extracts a predetermined number of models (methods of influencing factors) that are well-applicable (highly similar) to the target field and predicts the future state of the community. The predetermined number of models with good fit are those within a certain decision frame range from the maximum evaluation value (for example, within the range of 90% to 100% of the maximum evaluation value). The evaluation unit 140 may perform the evaluation using, for example, a genetic algorithm, a heuristic method such as an annealing method, a Monte Carlo method, or the like. Furthermore, the evaluation unit 140 may determine the amount of fertilizer applied to increase the crop yield based on the crop growth amount or activation index based on the derivation result of the state value derivation unit 130, or may determine the amount of fertilizer applied to increase the crop yield based on the crop growth amount or activation index, or prevent pest damage. It can be used for forecasting, determining the optimal number of plants to be planted at each sowing period, and determining fertilization management to prioritize the quality and grade of the harvest.

表示制御部150は、例えば、図4~10で示したように、少なくとも影響要因、遷移イベント、および生育形質指標の関係をグラフ化した画像を表示装置200に表示させる。また、表示制御部150は、図11で示したように、上記説明した各部の機能によって生成される、生育形質指標の変化をグラフ化した画像を表示装置200に表示させてもよい。表示装置200は、情報処理装置100に付設されるものであってもよいし、インターネットやLAN(Local Area Network)などのネットワークを介して情報処理装置100から画像が提供されるもの(例えば、パーソナルコンピュータ、タブレット端末、スマートフォンなど)であってもよい。 For example, as shown in FIGS. 4 to 10, the display control unit 150 causes the display device 200 to display an image in which the relationships among at least influencing factors, transition events, and growth trait indicators are graphed. Further, as shown in FIG. 11, the display control unit 150 may cause the display device 200 to display an image that is a graph of changes in the growth trait index, which is generated by the functions of each unit described above. The display device 200 may be attached to the information processing device 100, or may be a device that provides images from the information processing device 100 via a network such as the Internet or a LAN (Local Area Network) (for example, a personal computer). (computers, tablet terminals, smartphones, etc.).

以上説明した実施形態によれば、作物に与えられる条件および外乱要素を含む影響要因の情報を取得する取得部110と、影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも作物の生育形質指標ごとに導出する進行速度導出部120と、基準時点から予測時点までの進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出する状態値導出部130と、を備えることにより、作物の状態を精度よく予測することができる。 According to the embodiment described above, the acquisition unit 110 acquires information on influencing factors including conditions given to crops and disturbance elements, and the acquisition unit 110 acquires information on influencing factors including conditions given to crops and disturbance elements, and information on transition to the next growth stage based on the influencing factors and the number of growing days. A progress speed derivation unit 120 derives the progress speed for each growth trait index of the crop at least, and a vector in which integral values of the progress speed from the reference time point to the predicted time point are arranged from the output arc at the transition event of the growth stage. A state value deriving unit 130 that derives a predicted state vector including a state value for each growth trait index of the crop by multiplying by a connection matrix obtained by subtracting the input arc, thereby accurately predicting the state of the crop. be able to.

以上、本発明を実施するための形態について実施形態を用いて説明したが、本発明はこうした実施形態に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々の変形及び置換を加えることができる。 Although the mode for implementing the present invention has been described above using embodiments, the present invention is not limited to these embodiments in any way, and various modifications and substitutions can be made without departing from the gist of the present invention. can be added.

100 情報処理装置
110 取得部
120 進行速度導出部
130 状態値導出部
140 評価部
150 表示制御部
200 表示装置
100 Information processing device 110 Acquisition unit 120 Progressing speed derivation unit 130 State value derivation unit 140 Evaluation unit 150 Display control unit 200 Display device

Claims (8)

作物に与えられる条件および外乱要素を含む影響要因の情報を取得する取得部と、
前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出する進行速度導出部と、
基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出する状態値導出部と、
を備える情報処理装置。
an acquisition unit that acquires information on influencing factors including conditions and disturbance factors given to crops;
a progress rate derivation unit that derives a progress rate toward transition to the next growth stage for at least each growth trait index of the crop, based on the influencing factors and the number of growing days;
The state of each growth trait index of the crop can be determined by multiplying a vector listing the integral values of the progress speed from the reference time to the predicted time by a connection matrix obtained by subtracting the input arc from the output arc at the transition event of the growth stage. a state value derivation unit that derives a predicted state vector including the value;
An information processing device comprising:
前記進行速度導出部は、同じ生育形質指標に関する前記進行速度を、生育段階ごとに異なる傾向で導出し、
前記接続マトリクスは、生育段階を問わず前記積分値の用法を定義するものである、
請求項1記載の情報処理装置。
The progression rate derivation unit derives the progression rate regarding the same growth trait index with different trends for each growth stage,
The connection matrix defines the usage of the integral value regardless of the growth stage,
The information processing device according to claim 1.
前記取得部は、環境測定値に基づいて前記影響要因の情報を生成する処理を行う、
請求項1または2記載の情報処理装置。
The acquisition unit performs a process of generating information on the influencing factor based on environmental measurement values.
The information processing device according to claim 1 or 2.
前記取得部は、環境測定値に基づいて前記影響要因の情報を生成する処理の手法を変えながら予測状態ベクトルを導出し、
前記予測状態ベクトルの実際の作物への当てはまり度合いに基づいて前記影響要因の情報を生成する処理の手法を評価する評価部を更に備える、
請求項1から3のうちいずれか1項記載の情報処理装置。
The acquisition unit derives a predicted state vector while changing a processing method for generating information on the influencing factor based on the environmental measurement value,
further comprising an evaluation unit that evaluates a processing method for generating information on the influencing factors based on the degree of applicability of the predicted state vector to the actual crop;
An information processing device according to any one of claims 1 to 3.
少なくとも前記影響要因、前記遷移イベント、および前記生育形質指標の関係をグラフ化した画像を表示装置に表示させる表示制御部を更に備える、
請求項1から4のうちいずれか1項記載の情報処理装置。
further comprising a display control unit that causes a display device to display an image that is a graph of the relationship between at least the influencing factor, the transition event, and the growth trait index;
An information processing device according to any one of claims 1 to 4.
生育形質指標の変化をグラフ化した画像を表示装置に表示させる表示制御部を更に備える、
請求項1から5のうちいずれか1項記載の情報処理装置。
further comprising a display control unit that causes the display device to display an image in which changes in the growth trait index are graphed;
An information processing device according to any one of claims 1 to 5.
コンピュータが、
作物に与えられる条件および外乱要素を含む影響要因の情報を取得し、
前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出し、
基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出する、
情報処理方法。
The computer is
Obtain information on influencing factors including conditions and disturbance factors given to crops,
Based on the influencing factors and the number of growing days, the rate of progress toward transition to the next growth stage is derived for at least each growth trait index of the crop,
The state of each growth trait index of the crop can be determined by multiplying a vector listing the integral values of the progress speed from the reference time to the predicted time by a connection matrix obtained by subtracting the input arc from the output arc at the transition event of the growth stage. Derive a predicted state vector containing the values,
Information processing method.
コンピュータに、
作物に与えられる条件および外乱要素を含む影響要因の情報を取得させ、
前記影響要因および生育日数に基づいて、次の生育段階への遷移に向けて進行する進行速度を、少なくとも前記作物の生育形質指標ごとに導出させ、
基準時点から予測時点までの前記進行速度の積分値を並べたベクトルに、生育段階の遷移イベントにおける出力アークから入力アークを差し引いた接続マトリクスを乗算することで、前記作物の生育形質指標ごとの状態値を含む予測状態ベクトルを導出させる、
プログラム。
to the computer,
Obtain information on influencing factors including conditions and disturbance factors given to crops,
Based on the influencing factors and the number of growing days, the rate of progress toward transition to the next growth stage is derived for at least each growth trait index of the crop;
The state of each growth trait index of the crop can be determined by multiplying a vector listing the integral values of the progress speed from the reference time to the predicted time by a connection matrix obtained by subtracting the input arc from the output arc at the transition event of the growth stage. Derive a predicted state vector containing the values,
program.
JP2021039461A 2021-03-11 2021-03-11 Information processing device, information processing method, and program Active JP7385931B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021039461A JP7385931B2 (en) 2021-03-11 2021-03-11 Information processing device, information processing method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021039461A JP7385931B2 (en) 2021-03-11 2021-03-11 Information processing device, information processing method, and program

Publications (2)

Publication Number Publication Date
JP2022139192A JP2022139192A (en) 2022-09-26
JP7385931B2 true JP7385931B2 (en) 2023-11-24

Family

ID=83399096

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021039461A Active JP7385931B2 (en) 2021-03-11 2021-03-11 Information processing device, information processing method, and program

Country Status (1)

Country Link
JP (1) JP7385931B2 (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020218157A1 (en) 2019-04-25 2020-10-29 国立大学法人静岡大学 Prediction system, prediction method, and prediction program

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020218157A1 (en) 2019-04-25 2020-10-29 国立大学法人静岡大学 Prediction system, prediction method, and prediction program

Also Published As

Publication number Publication date
JP2022139192A (en) 2022-09-26

Similar Documents

Publication Publication Date Title
EP3550498A1 (en) Information processing device and information processing system
Nischitha et al. Crop prediction using machine learning approaches
Johnson et al. Can autocorrelated recruitment be estimated using integrated assessment models and how does it affect population forecasts?
US20190095556A1 (en) Information processing device, simulation method, and non-transitory recording medium storing simulation program
US20180173823A1 (en) Simulation device, simulation method, and recording medium for storing program
JP2017163956A (en) Prediction method, prediction program, and prediction apparatus
US20200250597A1 (en) Information processing apparatus, information processing method, and non-transitory recording medium
JP2019128904A (en) Prediction system, simulation system, method and program
Polacheck et al. An initial evaluation of management strategies for the southern bluefin tuna fishery
Lodovichi et al. Operational planning of herbicide-based weed management
Thierry et al. Simulating spatially-explicit crop dynamics of agricultural landscapes: The ATLAS simulator
CN111095314A (en) Yield estimation for crop plant planting
JP7385931B2 (en) Information processing device, information processing method, and program
JP2002136223A (en) Method for predicting crop yield of field crop, apparatus for predicting crop yield of field crop and recording medium
JP7059789B2 (en) Sequential control program, sequential control method and sequential control device
JP6908111B2 (en) Predictor, Prediction Method, Prediction Program, and Gene Estimator
JP6920151B2 (en) Forecasting device, forecasting method and forecasting program
De Oliveira et al. Potential improvements in the management of Bay of Biscay anchovy by incorporating environmental indices as recruitment predictors
Mertens et al. Data-based design of an intelligent control chart for the daily monitoring of the average egg weight
JP2022039406A (en) Information processing device, information processing method, and program
Leroux et al. Optimal investment in ecological rehabilitation under climate change
CN116579521A (en) Yield prediction time window determining method, device, equipment and readable storage medium
JP6570929B2 (en) Characteristic estimation model generation apparatus and method, analysis target characteristic estimation apparatus and method
JP7024733B2 (en) Information processing equipment, information processing method, and information processing program
Takimoto et al. Impact assessment and adaptation simulation for chalky rice grain in the cultivar ‘Koshihikari’in Japan using large ensemble climate projection data sets

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230626

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20230803

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: 20231024

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231106

R150 Certificate of patent or registration of utility model

Ref document number: 7385931

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150