Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7747324B2 - Agricultural support program, agricultural support method, and agricultural support device - Google Patents
[go: Go Back, main page]

JP7747324B2 - Agricultural support program, agricultural support method, and agricultural support device - Google Patents

Agricultural support program, agricultural support method, and agricultural support device

Info

Publication number
JP7747324B2
JP7747324B2 JP2021185380A JP2021185380A JP7747324B2 JP 7747324 B2 JP7747324 B2 JP 7747324B2 JP 2021185380 A JP2021185380 A JP 2021185380A JP 2021185380 A JP2021185380 A JP 2021185380A JP 7747324 B2 JP7747324 B2 JP 7747324B2
Authority
JP
Japan
Prior art keywords
cultivation
information
growth
parameters
leaf
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
JP2021185380A
Other languages
Japanese (ja)
Other versions
JP2022078987A (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
Publication of JP2022078987A publication Critical patent/JP2022078987A/en
Priority to JP2024194384A priority Critical patent/JP2025010441A/en
Application granted granted Critical
Publication of JP7747324B2 publication Critical patent/JP7747324B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
    • Y02A40/25Greenhouse technology, e.g. cooling systems therefor

Landscapes

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

Description

本発明は、農業支援プログラム、農業支援方法及び農業支援装置に関する。 The present invention relates to an agricultural support program, an agricultural support method, and an agricultural support device.

従来、農作物の情報が掲載されたカタログ等においては、各品種の特徴や特性を言葉や画像、平均値などで表現するのが一般的である。したがって、農作業従事者は、栽培する品種を選定する際には、言葉や画像、平均値などで表現された各品種の特徴や特性を参考にしている。 Traditionally, in catalogs and other publications containing information about agricultural crops, the characteristics and properties of each variety are typically expressed using words, images, average values, etc. Therefore, when selecting varieties to cultivate, farmers refer to the characteristics and properties of each variety expressed using words, images, average values, etc.

また、近年においては、ユーザが栽培地、栽培時期、及び品目を指定すると、指定された品目のうち指定された栽培地、時期において栽培できる品種を一覧表示する技術が知られている(例えば、特許文献1等参照)。 In recent years, a technology has become known that, when a user specifies a cultivation area, cultivation period, and item, displays a list of varieties of the specified item that can be cultivated in the specified cultivation area and period (see, for example, Patent Document 1).

特開2002-203002号公報Japanese Patent Application Laid-Open No. 2002-203002

しかしながら、カタログには、品種の一般的な特徴や特性しか掲載されていないため、栽培しようとしている場所や環境条件に合っているかが正確に把握できない。このため、実際に栽培してみると、想定していたように農作物が育たないおそれがある。また、特許文献1のように品種が一覧表示されたとしても、実際にどの品種を栽培すべきかを判断することは難しい。 However, because the catalog only lists the general characteristics and properties of each variety, it is not possible to accurately determine whether the variety is suitable for the intended cultivation location and environmental conditions. As a result, when actually cultivating crops, there is a risk that they will not grow as expected. Furthermore, even if a list of varieties is displayed, as in Patent Document 1, it is difficult to determine which variety should actually be cultivated.

本発明は、品種選定又は栽培管理を支援するための適切な情報を出力することが可能な農業支援プログラム、農業支援方法及び農業支援装置を提供することを目的とする。 The present invention aims to provide an agricultural support program, agricultural support method, and agricultural support device that can output appropriate information to support variety selection or cultivation management.

一つの態様では、農業支援プログラムは、ユーザが選択した農作物の品種と栽培地点との複数の組み合わせの情報を受け付け、前記栽培地点それぞれにおける農作物の栽培環境に関する情報を取得し、ユーザが入力した栽培方法に関する情報を取得し、前記品種それぞれの特徴を示すパラメータを記憶部から読み出し、前記記憶部から読み出した前記パラメータと、取得した前記栽培環境に関する情報と、取得した前記栽培方法に関する情報と、に基づいて、農作物の各器官の日々の生育を予測する生育モデル作成し、前生育モデルから、農作物の各器官の日々の生育を前記品種、栽培地点、栽培方法に関する情報の組み合わせである栽培条件ごとに予測し、前記予測により得られた前記農作物の各器官の日々の生育を示す情報を比較可能に出力し、ユーザによって入力された農作物の栽培結果を示す情報のうち、所定以上の信頼度を有する予め定められたユーザによって入力された情報に基づいて前記パラメータを更新し、前記記憶部に記憶する、処理をコンピュータに実行させるためのである。 In one aspect, the agricultural support program receives information on multiple combinations of crop varieties and cultivation locations selected by a user, acquires information on the cultivation environment of the crop at each of the cultivation locations , acquires information on the cultivation method entered by the user, reads parameters indicating the characteristics of each of the varieties from a memory unit, creates a growth model that predicts the daily growth of each organ of the crop based on the parameters read from the memory unit, the acquired information on the cultivation environment, and the acquired information on the cultivation method, predicts the daily growth of each organ of the crop from the growth model for each cultivation condition that is a combination of information on the variety, cultivation location, and cultivation method , outputs information indicating the daily growth of each organ of the crop obtained by the prediction in a comparable manner, updates the parameters based on information entered by a predetermined user that has a reliability of a predetermined level or higher from the information indicating the cultivation results of the crop entered by users, and stores the updated parameters in the memory unit.

本発明の農業支援プログラム、農業支援方法及び農業支援装置は、品種選定又は栽培管理を支援する適切な情報を出力することができるという効果を奏する。 The agricultural support program, agricultural support method, and agricultural support device of the present invention have the effect of being able to output appropriate information to support variety selection or cultivation management.

一実施形態に係る農業システムの構成を示す図である。FIG. 1 is a diagram illustrating a configuration of an agricultural system according to an embodiment. 図2(a)は、図1の利用者端末のハードウェア構成を示す図であり、図2(b)は、図1のサーバのハードウェア構成を示す図である。2A is a diagram showing the hardware configuration of the user terminal of FIG. 1, and FIG. 2B is a diagram showing the hardware configuration of the server of FIG. サーバの機能ブロック図である。FIG. 2 is a functional block diagram of a server. サーバの処理を示すフローチャートである。10 is a flowchart showing a process of the server. 生育モデルの概要を示す図である。FIG. 1 is a diagram showing an overview of a growth model. シミュレーションにおける作物サイズの変化を概念的に示す図である。FIG. 10 is a diagram conceptually illustrating changes in crop size in a simulation. 葉や茎について、生育段階と大きさの関係を相対値で示したものである。This shows the relationship between growth stage and size of leaves and stems as relative values. シミュレーションに使用する気象データの一例を示す図である。FIG. 10 is a diagram illustrating an example of meteorological data used in a simulation. シミュレーションで用いるパラメータの一例を示す図である。FIG. 10 is a diagram illustrating an example of parameters used in a simulation. 図10(a)、図10(b)は、本シミュレーションにおける葉の成長過程と果実の成長過程を概念的に示す図である。10(a) and 10(b) are diagrams conceptually showing the leaf growth process and the fruit growth process in this simulation. シミュレーション結果を表示する画面を示す図(その1)である。FIG. 10 is a diagram (part 1) showing a screen displaying a simulation result. シミュレーション結果を表示する画面を示す図(その2)である。FIG. 10 is a diagram (part 2) showing a screen displaying a simulation result. 図13(a)、図13(b)は、シミュレーション結果を表示する画面を示す図(その3)である。13(a) and 13(b) are diagrams (part 3) showing a screen displaying the simulation results. シミュレーション結果を表示する画面を示す図(その4)である。FIG. 4 is a diagram showing a screen displaying a simulation result. 更新部の処理の一例を示すフローチャートである。10 is a flowchart illustrating an example of processing by an update unit.

以下、農業システムの一実施形態について、図1~図15に基づいて詳細に説明する。図1には、一実施形態に係る農業システム100の構成が概略的に示されている。本実施形態の農業システム100は、農家等(以下「生産者」と呼ぶ)がトマト、イチゴ、キュウリ、パプリカなどの果菜類を栽培する際に、生産者に対して栽培品種の選択や栽培管理を支援する情報を提供するためのシステムである。なお、本実施形態では、イチゴを栽培する生産者の場合について説明する。 One embodiment of an agricultural system will be described in detail below with reference to Figures 1 to 15. Figure 1 shows a schematic diagram of the configuration of an agricultural system 100 according to one embodiment. The agricultural system 100 of this embodiment is a system that provides farmers and other such individuals (hereinafter referred to as "producers") with information to assist them in selecting cultivars and managing cultivation when cultivating fruit vegetables such as tomatoes, strawberries, cucumbers, and peppers. This embodiment will be described in the case of a producer cultivating strawberries.

農業システム100は、図1に示すように、農業支援装置としてのサーバ10と、利用者端末70と、を備える。利用者端末70は、生産者が利用するPC(Personal Computer)やタブレット型端末、スマートフォン等の端末である。サーバ10と利用者端末70は、インターネットなどのネットワーク80に接続されており、各装置間において情報のやり取りが可能となっている。 As shown in FIG. 1, the agricultural system 100 comprises a server 10 as an agricultural support device and a user terminal 70. The user terminal 70 is a terminal used by producers, such as a personal computer (PC), tablet terminal, or smartphone. The server 10 and user terminal 70 are connected to a network 80 such as the Internet, enabling the exchange of information between the devices.

利用者端末70は、生産者が入力した情報をサーバ10に対して送信する。図2(a)には、利用者端末70のハードウェア構成が示されている。図2(a)に示すように、利用者端末70は、CPU(Central Processing Unit)190、ROM(Read Only Memory)192、RAM(Random Access Memory)194、記憶部(ここではSSD(Solid State Drive)やHDD(Hard Disk Drive))196、ネットワークインタフェース197、表示部193、入力部195、及び可搬型記憶媒体191の読み取りが可能な可搬型記憶媒体用ドライブ199等を備えている。これら利用者端末70の構成各部は、バス198に接続されている。表示部193は、液晶ディスプレイ等を含み、入力部195は、キーボード、マウス、タッチパネル等を含む。 The user terminal 70 transmits information entered by the producer to the server 10. Figure 2(a) shows the hardware configuration of the user terminal 70. As shown in Figure 2(a), the user terminal 70 includes a CPU (Central Processing Unit) 190, ROM (Read Only Memory) 192, RAM (Random Access Memory) 194, a storage unit (here, an SSD (Solid State Drive) or HDD (Hard Disk Drive)) 196, a network interface 197, a display unit 193, an input unit 195, and a portable storage medium drive 199 capable of reading the portable storage medium 191. These components of the user terminal 70 are connected to a bus 198. The display unit 193 includes an LCD display or the like, and the input unit 195 includes a keyboard, mouse, touch panel, etc.

サーバ10は、利用者端末70から情報を取得し、取得した情報に基づいて、イチゴの品種選定や栽培管理を支援する情報を生成し、当該情報を表示する画面を、生産者が利用する利用者端末70に対して出力する装置である。 The server 10 is a device that acquires information from the user terminal 70, generates information to support strawberry variety selection and cultivation management based on the acquired information, and outputs a screen displaying this information to the user terminal 70 used by the producer.

図2(b)には、サーバ10のハードウェア構成が示されている。図2(b)に示すように、サーバ10は、コンピュータとしてのCPU90、ROM92、RAM94、記憶部(ここではSSDやHDD)96、ネットワークインタフェース97、及び可搬型記憶媒体用ドライブ99等を備えている。これらサーバ10の構成各部は、バス98に接続されている。サーバ10では、ROM92あるいはHDD96に格納されているプログラム(栽培補助プログラムを含む)、或いは可搬型記憶媒体用ドライブ99が可搬型記憶媒体91から読み取ったプログラム(栽培補助プログラムを含む)をCPU90が実行することにより、図4に示す各部の機能が実現される。なお、図4の各部の機能は、例えば、ASIC(Application Specific Integrated Circuit)やFPGA(Field Programmable Gate Array)等の集積回路により実現されてもよい。 Figure 2(b) shows the hardware configuration of the server 10. As shown in Figure 2(b), the server 10 includes a computer-like component such as a CPU 90, ROM 92, RAM 94, a storage unit (here, an SSD or HDD) 96, a network interface 97, and a portable storage medium drive 99. These components of the server 10 are connected to a bus 98. In the server 10, the CPU 90 executes programs (including the cultivation assistance program) stored in the ROM 92 or HDD 96, or programs (including the cultivation assistance program) read from the portable storage medium 91 by the portable storage medium drive 99, thereby realizing the functions of the components shown in Figure 4. Note that the functions of the components shown in Figure 4 may also be realized by integrated circuits such as an ASIC (Application Specific Integrated Circuit) or FPGA (Field Programmable Gate Array).

図3には、サーバ10の機能ブロック図が示されている。サーバ10においては、CPU90がプログラムを実行することにより、図3に示すように、選択受付部41、取得部としての環境データ取得部42、栽培情報取得部43、シミュレーション部44、出力部45、更新部46、として機能する。なお、図3には、サーバ10の記憶部96等に格納されている記憶部としてのパラメータDB50についても図示されている。 Figure 3 shows a functional block diagram of the server 10. In the server 10, the CPU 90 executes a program to function as a selection receiving unit 41, an environmental data acquisition unit 42 as an acquisition unit, a cultivation information acquisition unit 43, a simulation unit 44, an output unit 45, and an update unit 46, as shown in Figure 3. Note that Figure 3 also shows a parameter DB 50 as a memory unit stored in the memory unit 96 of the server 10.

選択受付部41は、生産者が選択した品種(例えば、とちおとめ、こいみのりなど)と栽培地点(例えば、つくば、盛岡、久留米など)の組み合わせの情報を受け付ける。なお、生産者は、品種と地点の組み合わせを1つ選択してもよいし、複数選択してもよい。選択受付部41は、受け付けた情報をシミュレーション部44に送信する。 The selection receiving unit 41 receives information on the combination of variety (e.g., Tochiotome, Koiminori, etc.) and cultivation location (e.g., Tsukuba, Morioka, Kurume, etc.) selected by the producer. Note that the producer may select one or more combinations of variety and location. The selection receiving unit 41 transmits the received information to the simulation unit 44.

環境データ取得部42は、選択受付部41が受け付けた栽培地点の情報に基づいて、栽培地点における環境データを取得する。なお、環境データ取得部42は、サーバ10内又はサーバ10以外の装置において管理されている過去の気象データや、将来の気象データ(予測データ)の中から、選択された栽培地点に対応するデータを取得する。気象データには、屋外の日射及び気温、ビニールハウス内の日射及び気温、湿度、CO2濃度、地温、土壌水分などのデータが含まれる。環境データ取得部42は、取得したデータをシミュレーション部44に送信する。 The environmental data acquisition unit 42 acquires environmental data for the cultivation location based on the information about the cultivation location accepted by the selection acceptance unit 41. The environmental data acquisition unit 42 acquires data corresponding to the selected cultivation location from past weather data and future weather data (forecast data) managed in the server 10 or in a device other than the server 10. The weather data includes data such as outdoor solar radiation and air temperature, solar radiation and air temperature inside the greenhouse, humidity, CO2 concentration, soil temperature, and soil moisture. The environmental data acquisition unit 42 transmits the acquired data to the simulation unit 44.

栽培情報取得部43は、生産者が入力した栽培情報(例えば定植日、栽培密度、土耕/溶液栽培、定植時葉数、施肥量、培養液濃度等)を取得する。栽培情報取得部43は、取得した情報をシミュレーション部44に送信する。 The cultivation information acquisition unit 43 acquires cultivation information entered by the producer (e.g., planting date, cultivation density, soil/solution cultivation, number of leaves at planting, amount of fertilizer applied, concentration of nutrient solution, etc.). The cultivation information acquisition unit 43 transmits the acquired information to the simulation unit 44.

シミュレーション部44は、選択受付部41が受け付けた品種に対応するパラメータをパラメータDB50から読み出し、読み出したパラメータと、環境データ取得部42及び栽培情報取得部43から送信されてくるデータと、を用いて、生育モデルを作成する。ここで、パラメータDB50に格納されているパラメータは、例えば、図9に示すようなパラメータであり、品種ごとに定められているものとする。そして、シミュレーション部44は、作成した生育モデルを用いて、生産者が選択した品種を、生産者が選択した栽培地点、栽培方法で栽培したときの農作物の生育に関するシミュレーションを実行する。シミュレーション部44のシミュレーションの結果としては、葉面積、花房別の開花日、花房別の収量、果実乾物分配率、光合成量、生育量(茎葉)、生育量(果実)、養分吸収量(施肥量)等が得られる。シミュレーション部44は、シミュレーション結果を出力部45に送信する。 The simulation unit 44 reads parameters corresponding to the variety accepted by the selection acceptance unit 41 from the parameter DB 50, and creates a growth model using the read parameters and data transmitted from the environmental data acquisition unit 42 and the cultivation information acquisition unit 43. Here, the parameters stored in the parameter DB 50 are, for example, parameters such as those shown in FIG. 9, and are assumed to be determined for each variety. The simulation unit 44 then uses the created growth model to perform a simulation of the growth of the crop when the variety selected by the producer is cultivated at the cultivation location and using the cultivation method selected by the producer. The simulation results obtained by the simulation unit 44 include leaf area, flowering date per inflorescence, yield per inflorescence, fruit dry matter distribution rate, photosynthesis rate, growth amount (stems and leaves), growth amount (fruit), and nutrient absorption amount (fertilizer amount). The simulation unit 44 transmits the simulation results to the output unit 45.

出力部45は、シミュレーション部44から受信したシミュレーション結果を表示する画面を生成し、生産者が利用する利用者端末70に対して送信する。このとき、出力部45は、栽培条件(品種・栽培地点・栽培情報の組み合わせ)ごとにシミュレーション結果を纏め、異なる栽培条件のシミュレーション結果を比較可能な状態で表示する画面を生成する。 The output unit 45 generates a screen displaying the simulation results received from the simulation unit 44 and transmits them to the user terminal 70 used by the producer. At this time, the output unit 45 compiles the simulation results for each cultivation condition (combination of variety, cultivation location, and cultivation information) and generates a screen displaying the simulation results for different cultivation conditions in a manner that allows them to be compared.

更新部46は、生産者から、実際に栽培した品種の情報や、栽培結果(例えば、葉面積、花房別の開花日、花房別の収量など)を取得する。そして、更新部46は、実際の栽培結果と、シミュレーション結果とが近づくように、パラメータDB50に格納されているパラメータを更新する。なお、更新部46は、実際の栽培結果を入力した生産者に対して予め定められている信頼度に基づいて、当該実際の栽培結果を用いたパラメータの更新を行うか否かを判断することとしてもよい。 The update unit 46 acquires information on the varieties actually cultivated and cultivation results (e.g., leaf area, flowering date per inflorescence, yield per inflorescence, etc.) from the producer. The update unit 46 then updates the parameters stored in the parameter DB 50 so that the actual cultivation results and the simulation results are closer. Note that the update unit 46 may determine whether or not to update the parameters using the actual cultivation results based on a predetermined reliability for the producer who input the actual cultivation results.

(サーバ10の処理について)
次に、サーバ10の処理の詳細について説明する。
(Regarding the processing of the server 10)
Next, the processing of the server 10 will be described in detail.

図4には、サーバ10の処理がフローチャートにて示されている。図4の処理では、まず、ステップS10において、選択受付部41は、生産者が利用者端末70上で品種と栽培地点の組み合わせを選択するまで待機する。ここでは、例えば、生産者により、「とちおとめ・つくば」、「とちおとめ・盛岡」、「こいみのり・つくば」、「おいCベリー・つくば」、「さちのか・久留米」の各組み合わせが選択されたものとする。 Figure 4 shows the processing of the server 10 in a flowchart. In the processing of Figure 4, first, in step S10, the selection receiving unit 41 waits until the producer selects a combination of variety and cultivation location on the user terminal 70. Here, for example, it is assumed that the producer has selected the combinations of "Tochiotome - Tsukuba," "Tochiotome - Morioka," "Koiminori - Tsukuba," "Oi C Berry - Tsukuba," and "Sachinoka - Kurume."

ステップS10の判断が肯定されると、次のステップS12では、選択受付部41は、選択された品種と栽培地点の組み合わせを取得し、環境データ取得部42とシミュレーション部44に送信する。 If the determination in step S10 is positive, in the next step S12, the selection reception unit 41 acquires the selected combination of variety and cultivation location and transmits it to the environmental data acquisition unit 42 and the simulation unit 44.

次いで、ステップS14では、環境データ取得部42は、選択受付部41で受け付けられた、生産者によって選択された栽培地点に基づいて、当該栽培地点における気象データ(過去データや予測データ)を取得する。 Next, in step S14, the environmental data acquisition unit 42 acquires weather data (past data and forecast data) for the cultivation location based on the cultivation location selected by the producer and accepted by the selection acceptance unit 41.

次いで、ステップS16では、栽培情報取得部43は、生産者によって栽培情報(例えば、定植時期、定植密度、土耕/養分栽培、定植時の葉の枚数)が入力されるまで待機する。生産者が栽培情報を入力すると、ステップS18に移行し、栽培情報取得部43は、入力された栽培情報を取得し、シミュレーション部44に送信する。 Next, in step S16, the cultivation information acquisition unit 43 waits until the producer inputs cultivation information (e.g., planting time, planting density, soil/nutrient cultivation, number of leaves at planting). Once the producer inputs the cultivation information, the process proceeds to step S18, where the cultivation information acquisition unit 43 acquires the input cultivation information and sends it to the simulation unit 44.

次いで、ステップS20では、シミュレーション部44は、パラメータDB50に格納されている、生産者によって選択された各品種のパラメータを読み出し、当該パラメータと環境データ取得部42が取得した気象データ、栽培情報取得部43が取得した栽培情報、に基づいて、選択された品種と栽培地点の組み合わせに対応する生育モデルを作成する。 Next, in step S20, the simulation unit 44 reads the parameters for each variety selected by the producer that are stored in the parameter DB 50, and creates a growth model corresponding to the combination of the selected variety and cultivation location based on the parameters, the weather data acquired by the environmental data acquisition unit 42, and the cultivation information acquired by the cultivation information acquisition unit 43.

そして、シミュレーション部44は、作成した生育モデルを用いてシミュレーションを実行し、各品種が各栽培地点でどのように生育するかを示すシミュレーション結果を得る。シミュレーション部44が作成する生育モデルやシミュレーション結果の詳細については後述する。シミュレーション部44は、シミュレーション結果を出力部45に送信する。 The simulation unit 44 then runs a simulation using the created growth model and obtains simulation results that indicate how each variety will grow at each cultivation location. Details of the growth model created by the simulation unit 44 and the simulation results will be described later. The simulation unit 44 transmits the simulation results to the output unit 45.

次いで、ステップS22では、出力部45が、シミュレーション結果を表示する画面を生成し、利用者端末70に送信する。シミュレーション結果を表示する画面は、例えば、図11~図13に示すような画面である。 Next, in step S22, the output unit 45 generates a screen displaying the simulation results and transmits it to the user terminal 70. The screen displaying the simulation results is, for example, a screen such as those shown in Figures 11 to 13.

(シミュレーションについて)
以下、シミュレーション部44が実行するシミュレーションについて、詳細に説明する。
(About the simulation)
The simulation performed by the simulation unit 44 will be described in detail below.

(1)シミュレーションの基本的な考え方
本実施形態においては、品種間の重要な違いを説明することが可能な生育モデルを利用する。具体的には、図5に示すような、光合成に関する特性、光合成産物の分配、葉の発生と伸長、花の発生と肥大、養分吸収、形態的特性(葉の形状、受光態勢)等を説明することが可能な生育モデルを利用する。なお、本実施形態では、イチゴを対象としたシミュレーションについて説明するが、他の作物への転用できるようにするため、パラメータの意味が明確になるように設計している。
(1) Basic Concept of Simulation In this embodiment, a growth model capable of explaining important differences between varieties is used. Specifically, as shown in Fig. 5, a growth model capable of explaining photosynthetic characteristics, distribution of photosynthetic products, leaf development and elongation, flower development and enlargement, nutrient absorption, morphological characteristics (leaf shape, light reception posture), etc. is used. Note that although this embodiment describes a simulation targeting strawberry, the model is designed to clearly indicate the meaning of parameters so that it can be applied to other crops.

(a)成長パターンについて
シミュレーションにおいては、生育モデルに対して日単位で気象データ(日平均気温、および日屋外全天日射量)を与え、作物の器官(根、葉、クラウン、果実)ごとの生育状況を計算する。ここで、作物は実際には、常に生育状態にあり、作物サイズ(根、葉、クラウンのサイズ、果実数など)が常に変化するが、シミュレーションにおいては、0時で作物サイズを固定し、その状態で0~24時までの気温、日射を受けると仮定している。この仮定による誤差を減らしたい場合には、シミュレーションの計算間隔(本シミュレーションの計算間隔は1日)を短くするなどすればよい。
(a) Growth Pattern In the simulation, weather data (daily average temperature and daily outdoor global solar radiation) is provided to the growth model on a daily basis, and the growth status of each crop organ (roots, leaves, crown, fruit) is calculated. In reality, crops are always growing, and crop size (size of roots, leaves, crown, number of fruits, etc.) is always changing, but in the simulation, it is assumed that the crop size is fixed at midnight, and that the temperature and solar radiation are received in that state from midnight to midnight. To reduce errors due to this assumption, the calculation interval of the simulation (the calculation interval in this simulation is one day) can be shortened, for example.

本シミュレーションでは、N日目の計算を行う場合、N日0時にその日の作物サイズを決定する。すなわち、本シミュレーションでは、作物がN日目0~24時に気温、日射量を受け取るための状態(相対受光率:葉面積指数から求められる値)、および分配率(当日生成される光合成産物を葉、クラウン、花房、根の各器官に分配する割合)が固定されていることを意味している。そして、本シミュレーションでは、N日目24時の時点で、当日(N日目)の気温、日射量、CO2濃度により、光合成量が計算され、その値から作物の成長計算により作物サイズが変化する。図6には、本シミュレーションにおける作物サイズの変化が概念的に示されている。なお、計算の1日目については、0時における作物サイズ(0日目作物サイズ)が計算できないため、生産者等が値を設定する。 In this simulation, when calculating the Nth day, the crop size for that day is determined at midnight on the Nth day. This means that the temperature, solar radiation conditions (relative light interception rate: a value calculated from the leaf area index), and distribution rate (the proportion of photosynthetic products generated that day distributed to each organ: leaves, crown, inflorescence, and roots) for the crop are fixed from midnight to midnight on the Nth day. In this simulation, the amount of photosynthesis is calculated based on the temperature, solar radiation, and CO2 concentration for that day (day N), and the crop size changes based on this value through crop growth calculations. Figure 6 conceptually illustrates the change in crop size in this simulation. Note that for the first day of the calculation, the crop size at midnight (day 0 crop size) cannot be calculated, so the producer or other person sets the value.

(b)光合成産物量とサイズの関係
本シミュレーションにおいて、光合成産物を「ソース」と呼び、ソースを生産する器官を「ソース器官」と呼ぶ。ソース器官は、主に葉である。一方、ソースの貯蔵、および消費を行う器官を、「シンク器官」と呼ぶ。シンク器官は、根、クラウン、葉、果実であるが、最も消費量の多いシンク器官は果実である。
(b) Relationship between photosynthetic product amount and size In this simulation, photosynthetic products are called "sources" and organs that produce the sources are called "source organs." Source organs are mainly leaves. On the other hand, organs that store and consume the sources are called "sink organs." Sink organs are roots, crowns, leaves, and fruits, but the sink organ with the highest consumption is the fruit.

本シミュレーションでは、作物の成長は主として気温によって決まるとしている。例えば、一枚の葉の成長過程は積算気温を用いて次のように記述する。成長が停止するときの積算気温を1.0、成長を開始時の積算気温を0として、葉の生育段階を積算気温の相対値として示す。これを[Leaf Growth Index]と定義する。 In this simulation, crop growth is assumed to be primarily determined by temperature. For example, the growth process of a single leaf is described as follows using accumulated temperature. The accumulated temperature when growth stops is set to 1.0, and the accumulated temperature when growth begins is set to 0, and the leaf growth stage is shown as a relative value of the accumulated temperature. This is defined as the [Leaf Growth Index].

ここで、[Leaf Growth Index]に対応する相対的な葉の大きさ(サイズ)を [Leaf Size Index]と定義する(図7)。[Leaf Size Index]は、完全に成長しきった状態の大きさを1.0、発育開始前の大きさを0とする。 Here, the relative leaf size (size) corresponding to the [Leaf Growth Index] is defined as the [Leaf Size Index] (Figure 7). The [Leaf Size Index] is set to 1.0 for the size of a fully grown leaf, and 0 for the size before growth begins.

本シミュレーションでは、葉や果実など生育をシグモイド曲線で近似する。[Leaf Growth Index]と[Leaf Size Index]の関係を式(1)で表す(図7)。xは[Leaf Growth Index]で、S(x)は[Leaf Size Index]である。 In this simulation, the growth of leaves, fruits, etc. is approximated by a sigmoid curve. The relationship between [Leaf Growth Index] and [Leaf Size Index] is expressed by equation (1) (Figure 7). x is [Leaf Growth Index] and S(x) is [Leaf Size Index].

N日目の0時の時点と24時の時点におけるS(x)の値の差(図7参照)を、一日当たりの葉の成長増加量の指標値として[Leaf Increase Index]と定義する。花房についても同様に花房の相対成長量[Fruits Increase Index]を定義する。なお、本シミュレーション用のシグモイド曲線の値は、x=0でS(x)=0.007、x=1でS(x)=0.993となるため、x=0以下でS(x)=0、x=1以上でS(x)=1と変更して使用する。 The difference between the S(x) values at midnight and midnight on the Nth day (see Figure 7) is defined as the [Leaf Increase Index], an index of the amount of leaf growth increase per day. Similarly, the relative growth of the inflorescence is defined as the [Fruits Increase Index]. Note that the values of the sigmoid curve used in this simulation are S(x) = 0.007 at x = 0 and S(x) = 0.993 at x = 1, so we will use the values changed to S(x) = 0 below x = 0 and S(x) = 1 above x = 1.

(c)シミュレーションモデルの概要
シミュレーションは、気象データ及びパラメータ(定数)に基づいて実行する。シミュレーションに使用する気象データは、図8に示すように、日平均気温(℃)、日屋外全天日射量(MJ/m2)、CO2濃度(ppm)などとする。生産者等のユーザが設定するパラメータの例を図9に示す。本シミュレーションは、日単位で計算が実行され、葉一枚ごと、一花房ごとに計算される値、および作物全体に対して計算される値がある。生育は、作物の器官(葉、クラウン、果実、根、花房)ごとの状態が計算される。
(c) Overview of the Simulation Model The simulation is performed based on meteorological data and parameters (constants). The meteorological data used in the simulation includes the average daily temperature (°C), the amount of outdoor solar radiation (MJ/ m2 ), and the CO2 concentration (ppm), as shown in Figure 8. An example of parameters set by users such as producers is shown in Figure 9. This simulation is performed on a daily basis, with values calculated for each leaf and each inflorescence, as well as values calculated for the entire crop. Growth is calculated for the state of each organ of the crop (leaves, crown, fruit, roots, and inflorescence).

(2)シミュレーション方法について
(a)シミュレーションの概要
本シミュレーションは、大きく分けて、「光合成量計算」、「作物成長計算」、「果実収穫量計算」および「養分吸収量」のプロセスから構成される。これらの計算は、各日単位で行う。光合成量計算に使用される変数(相対受光量、分配率)は、計算日0時(計算日前日)における値である。相対受光量は、日射量から作物が光合成に使用する受光量を求めるための変数である。分配率は、その日に作物が合成するソースが各器官(根、葉、クラウン、果実)へ分配される割合である。分配率、およびこれらの値、および当日の気温、日射量により、計算日24時における光合成量が計算される。そして、光合成量をもとに計算日24時における作物の成長度が計算される。作物の成長量計算は、(a)分配量計算、(b)葉面積計算、(c)着果数計算、(d)果実乾物重量計算、(e)分配率計算に分かれる。
(2) Simulation Method: (a) Simulation Overview. This simulation is broadly divided into the following processes: photosynthesis calculation, crop growth calculation, fruit yield calculation, and nutrient absorption calculation. These calculations are performed on a daily basis. The variables used to calculate photosynthesis (relative light reception, distribution rate) are values at midnight on the calculation day (the day before the calculation day). Relative light reception is a variable used to calculate the amount of light received by the crop for photosynthesis from the amount of solar radiation. The distribution rate is the proportion of the source synthesized by the crop on that day distributed to each organ (root, leaf, crown, fruit). The photosynthesis rate at 24:00 on the calculation day is calculated using the distribution rate, these values, and the temperature and solar radiation on that day. The crop growth rate at 24:00 on the calculation day is then calculated based on the photosynthesis rate. The crop growth rate calculation is divided into (a) distribution rate calculation, (b) leaf area calculation, (c) fruit number calculation, (d) fruit dry weight calculation, and (e) distribution rate calculation.

なお、本明細書においては、変数名に付した下付き添字Nは、N日24時の結果を示すものとする。また、葉一枚ごとに計算する変数には下付き添字M(葉の順位)、花房ごとに計算する変数には下付き添字F(花房順位)を付すものとする。また、式で使用する変数、定数は、[]を付して示すものとする。なお、本シミュレーションでは、計算1日目0時の値を「初期値」とする。 In this specification, the subscript N attached to a variable name indicates the result at 24:00 on the Nth day. Variables calculated for each leaf are marked with the subscript M (leaf rank), and variables calculated for each inflorescence are marked with the subscript F (inflorescence rank). Variables and constants used in formulas are indicated with square brackets ([ ]). In this simulation, the value at 0:00 on the first day of calculation is considered the "initial value."

(b)作物成長過程
本シミュレーションでは、作物の成長は、前の葉、前の花房順位(葉や花房の発生順序。発生が早いほど順位は小さい)との関係、および積算温度により決定する。図10(a)、図10(b)は、本シミュレーションにおける葉と花房の成長過程を概念的に示す図である(以下に記載した数値は、架空の品種を想定した架空の値であり、実際の値ではない)。図10(a)に示すように、葉は、前の葉の成長開始時を起点として、積算温度が150℃になってから、次の葉の成長が開始する。この成長開始時を起点として、積算温度が450℃になった時点で葉の成長が終了する。一方、花房は、図10(b)に示すように、前順位花房の花芽分化開始時を起点として、積算温度が600℃になってから、次花房の花芽分化が開始する。この花芽分化を起点として、積算温度が600℃になった時点で開花期となり、花房を構成する果実の肥大が開始する。更に、開花時を起点として、積算温度が600℃になった時点で花房の成長が終了する。また、図10(b)に示すように、花芽分化開始からの積算温度150℃までの期間が花数決定期間であり、この期間の条件で着果数が決定する。
(b) Crop Growth Process In this simulation, crop growth is determined by the relationship between the previous leaf and the order of the previous inflorescence (the order in which leaves and inflorescences develop; the earlier they develop, the lower the order), and the accumulated temperature. Figures 10(a) and 10(b) are conceptual diagrams illustrating the leaf and inflorescence growth process in this simulation (the values listed below are fictitious values assuming a fictitious variety and are not actual values). As shown in Figure 10(a), the growth of the next leaf begins when the accumulated temperature reaches 150°C, starting from the start of growth of the previous leaf. Leaf growth ends when the accumulated temperature reaches 450°C, starting from the start of flower bud differentiation of the previous inflorescence. Meanwhile, as shown in Figure 10(b), the flower bud differentiation of the next inflorescence begins when the accumulated temperature reaches 600°C, starting from this flower bud differentiation. Starting from this flower bud differentiation, the flowering period begins when the accumulated temperature reaches 600°C, and the fruit that makes up the inflorescence begins to enlarge. Furthermore, starting from flowering, the growth of the inflorescence ends when the accumulated temperature reaches 600° C. As shown in FIG. 10( b), the period from the start of flower bud differentiation to an accumulated temperature of 150° C. is the flower number determination period, and the number of fruits is determined by the conditions during this period.

(c)光合成量計算
光利用効率LUE(gDW/MJ)は、受光量あたりの作物全体の乾物生産量を示す値である。気温またはCO2濃度の関数として与えられる。
(c) Calculation of photosynthesis amount Light use efficiency (LUE) (gDW/MJ) is a value that indicates the dry matter production of the entire crop per unit of received light. It is given as a function of temperature or CO2 concentration.

またハウス内の日射量[ハウス内全天日射量](MJ/m2)は、[屋外全天日射量](MJ/m2)にそのハウスの[日射透過率](定数)を乗じて求める(式(2))。
[ハウス内全天日射量]N=[日射透過率]・[屋外全天日射量]N …(2)
The amount of solar radiation inside the greenhouse [global solar radiation inside the greenhouse] (MJ/m 2 ) is calculated by multiplying the [global solar radiation outside] (MJ/m 2 ) by the [solar transmittance] (constant) of the greenhouse (equation (2)).
[Indoor global solar radiation] N = [Solar transmittance] × [Outdoor global solar radiation] N … (2)

また、作物が光合成に使用する[受光量](MJ/m2)は、次式(3)で求める。なお、相対受光率は、前日の葉面積から求められる値が用いられる。
[受光量]N=[相対受光量](N-1)・[ハウス内全天日射量]N …(3)
The amount of light received (MJ/m 2 ) used by crops for photosynthesis is calculated using the following formula (3): The relative light receiving rate is calculated from the leaf area on the previous day.
[Amount of received light] N = [Relative amount of received light] (N-1) × [Total amount of solar radiation in the greenhouse] N … (3)

更に、単位面積当たりの作物全体の光合成量である[乾物生産量](gDW/m2)は、次式(4)で求める。
[乾物生産量]N=[受光量]N・[LUE]N …(4)
Furthermore, the amount of photosynthesis of the entire crop per unit area, [dry matter production] (gDW/m 2 ), is calculated using the following formula (4).
[Dry matter production] N = [Amount of light received] N · [LUE] N ... (4)

また、作物1株当たりの光合成量である[光合成量](gDW/株)は、単位面積当たりの作物本数(定数)である[Plant Density](株/m2)を用い、次式(5)で求める。
[光合成量]N=[乾物生産量]N/[Plant Density] …(5)
The amount of photosynthesis per plant, [Photosynthesis Amount] (gDW/plant), is calculated using the number of plants per unit area (constant), [Plant Density] (plant/m 2 ), using the following formula (5).
[Photosynthesis rate] N = [Dry matter production rate] N / [Plant density] ... (5)

また、後で行う葉面積計算では、葉の乾物重量当たりの葉面積である比葉面積(Specific Leaf Area)である[SLA](m2/gDW)を用いる。[SLA]Nは気温の関数として表す。 In addition, in the leaf area calculations that will be performed later, we will use [SLA] (m 2 /gDW), which is the specific leaf area, which is the leaf area per dry weight of the leaf. [SLA] N is expressed as a function of temperature.

(3)作物成長計算について
(a)分配量計算
分配量計算では、生産された光合成量の各器官への分配量を求める。分配量は、作物全体の値として計算され、基本的に光合成量と分配率により求められるが、花房が受け取れるソースの最大値(分配量最大値)として[花房ポテンシャル成長量](gDW/株)を設定し補正を行う。[花房ポテンシャル成長量]は、次式(6)で求めることができる。なお、式(6)の[花房成長量係数](gDW/(株・℃))は、気温1℃に対し一花房が受けとれるソースの最大値を示す係数であり、[Total Fruits Increase Index]は、全花房の相対成長量の合計値である。
[花房ポテンシャル成長量]N
=[花房成長量係数]・[気温]N・[Total Fruits Increase Index]N-1 …(6)
(3) Crop growth calculations (a) Allocation calculation Allocation calculations determine the amount of photosynthesis produced distributed to each organ. Allocation is calculated as a value for the entire crop, and is basically determined from the amount of photosynthesis and the distribution rate. However, a correction is made by setting [inflorescence potential growth] (gDW/plant) as the maximum amount of source that can be received by an inflorescence (maximum amount of allocation). [Inflorescence potential growth] can be calculated using the following formula (6). Note that the [inflorescence growth coefficient] (gDW/(plant·°C)) in formula (6) is a coefficient that indicates the maximum amount of source that one inflorescence can receive per 1°C of temperature, and the [Total Fruits Increase Index] is the sum of the relative growth amounts of all inflorescences.
[Flower cluster potential growth] N
= [Flower cluster growth coefficient] × [Temperature] N × [Total Fruits Increase Index] N-1 … (6)

分配量は、3段階で計算する。このため、最初の2段階の変数には(1)、(2)を付し、最終結果と区別して表記するものとする。光合成量、分配量の単位は(gDW/株)である。[葉分配量(1)]、[クラウン分配量(1)]、[花房分配量(1)]、[根分配量(1)](gDW/株)は、前日(当日0時)の分配率に従い、当日0~24時の光合成量が各器官に割り振られるとして、次式(7A)~(7D)にて求める。
[葉分配量(1)]N=[葉分配率]N-1・[光合成量]N …(7A)
[クラウン分配量(1)]N=[クラウン分配率]N-1・[光合成量]N …(7B)
[花房分配量(1)]N=[花房分配率]N-1・[光合成量]N …(7C)
[根分配量(1)]N=[根分配率]N-1・[光合成量]N …(7D)
Allocation is calculated in three stages. For this reason, the variables for the first two stages are labeled (1) and (2) to distinguish them from the final results. The unit for photosynthesis and allocation is (gDW/plant). [Leaf allocation (1)], [Crown allocation (1)], [Inflorescence allocation (1)], [Root allocation (1)] (gDW/plant) are calculated using the following formulas (7A) to (7D), assuming that the amount of photosynthesis from 0:00 to 24:00 on the day is allocated to each organ according to the allocation rate from the previous day (0:00 on the day).
[Leaf distribution amount (1)] N = [Leaf distribution rate] N-1・[Photosynthesis amount] N … (7A)
[Crown distribution amount (1)] N = [Crown distribution rate] N-1 × [Photosynthetic amount] N … (7B)
[Inflorescence distribution amount (1)] N = [Inflorescence distribution rate] N-1・[Photosynthesis amount] N … (7C)
[Root distribution amount (1)] N = [Root distribution rate] N-1・[Photosynthesis amount] N … (7D)

また、[葉分配量(2)]、[クラウン分配量(2)]、[根分配量(2)](gDW/株)は、花房分配量が花房ポテンシャル成長量を超えないように以下のように補正する。
[葉分配量(2)]N=[葉分配量(1)]N …(7A)’
[クラウン分配量(2)]N=[クラウン分配量(1)]N …(7B)’
[根分配量(2)]N=[根分配量(1)]N …(7D)’
In addition, [Leaf distribution amount (2)], [Crown distribution amount (2)], and [Root distribution amount (2)] (gDW/plant) are corrected as follows so that the inflorescence distribution amount does not exceed the inflorescence potential growth amount.
[Leaf distribution amount (2)] N = [Leaf distribution amount (1)] N … (7A)'
[Crown distribution amount (2)] N = [Crown distribution amount (1)] N ... (7B)'
[Root distribution amount (2)] N = [Root distribution amount (1)] N … (7D)'

なお、[花房分配量(2)]については、[花房分配量(1)]N<[花房ポテンシャル成長量]Nであれば、
[花房分配量(2)N=[花房分配量(1)]N …(7C1)’
とし、その他の場合には、
[花房分配量(2)]N=[花房ポテンシャル成長量]N …(7C2)’
とする。
Regarding [Flower cluster allocation amount (2)], if [Flower cluster allocation amount (1)] N < [Flower cluster potential growth amount] N ,
[Inflorescence distribution amount (2) N = [Inflorescence distribution amount (1)] N … (7C1)'
In other cases,
[Flower cluster distribution amount (2)] N = [Flower cluster potential growth amount] N ... (7C2)'
Let's say.

ここで、花房ポテンシャル成長量を超えたため、花房に分配されなかった光合成量である[余剰乾物量](gDW/株)は、次式(8)にて表される。
[余剰乾物量]N=([葉分配量(1)]N+[クラウン分配量(1)]N+[花房分配量(1)]N+[根分配量(1)]N)-([葉分配量(2)]N+[クラウン分配量(2)]N+[花房分配量(2)]N+[根分配量(2)]N) …(8)
Here, the amount of photosynthesis that was not distributed to the inflorescence because it exceeded the inflorescence potential growth amount, [excess dry matter] (gDW/plant), is expressed by the following formula (8).
[Excess dry matter] N = ([Leaf allocation (1)] N + [Crown allocation (1)] N + [Inflorescence allocation (1)] N + [Root allocation (1)] N ) - ([Leaf allocation (2)] N + [Crown allocation (2)] N + [Inflorescence allocation (2)] N + [Root allocation (2)] N ) ... (8)

余剰乾物量を花房以外に再配分し、それぞれの分配量を次式(9A)~(9D)により求める。
[葉分配量]N=[葉分配量(2)]N+[葉余剰分配率]・[余剰乾物量]N …(9A)
[クラウン分配量]N
=[クラウン分配量(2)]N+[クラウン余剰分配率]・[余剰乾物量]N …(9B)
[花房分配量]N=[花房分配量(2)]N …(9C)
[根分配量]N=[根分配量(2)]N+[根余剰分配率]・[余剰乾物量]N …(9D)
The excess dry matter is redistributed to areas other than the inflorescence, and the respective distribution amounts are calculated using the following formulas (9A) to (9D).
[Leaf distribution amount] N = [Leaf distribution amount (2)] N + [Leaf surplus distribution rate] × [Excess dry matter amount] N … (9A)
[Crown distribution amount] N
= [crown distribution amount (2)] N + [crown surplus distribution rate] · [surplus dry matter amount] N ... (9B)
[Inflorescence distribution amount] N = [Inflorescence distribution amount (2)] N … (9C)
[Root distribution amount] N = [Root distribution amount (2)] N + [Root surplus distribution rate] · [Excess dry matter amount] N ... (9D)

(b)葉面積計算
葉面積計算では葉一枚ごとに面積、および葉、クラウンの乾物重量(DW)を計算する。まず、葉の生育段階を積算気温の相対値として示す[Leaf Growth Index](無次元)を次式(10A)で求める。なお、Mは葉の順位である。
[Leaf Growth Index]M,N
={[積算温度]N-(M-1)・150}/450 …(10A)
(b) Leaf Area Calculation The leaf area calculation involves calculating the area of each leaf and the dry weight (DW) of the leaf and crown. First, the Leaf Growth Index (dimensionless), which indicates the leaf growth stage as a relative value of the accumulated temperature, is calculated using the following equation (10A). M is the leaf rank.
[Leaf Growth Index] M,N
= {[Accumulated temperature] N - (M - 1) 150} / 450 ... (10 A)

葉のサイズ指数である[Leaf Size Index](無次元)を次式(10B)で計算する。なお、[Leaf Size Index]の初期値は、ゼロとする。 The leaf size index, [Leaf Size Index] (dimensionless), is calculated using the following formula (10B). Note that the initial value of [Leaf Size Index] is set to zero.

また、前日24時(当日0時)から当日24時における葉のサイズ指数変化から、葉の成長量指数である[Leaf Increase Index](無次元)を次式(11)で求める。
[Leaf Increase Index]M,N
=[Leaf Size Index]M,N-[Leaf Size Index]M,N-1 …(11)
In addition, the leaf growth index [Leaf Increase Index] (dimensionless) is calculated from the change in leaf size index from midnight of the previous day (midnight of the current day) to midnight of the current day using the following equation (11).
[Leaf Increase Index] M,N
=[Leaf Size Index] M,N -[Leaf Size Index] M,N-1 ...(11)

また、全ての葉の[Leaf Increase Index]の総計を[Total Leaf Increase Index]とし、[Total Leaf Increase Index]を次式(12)から求める。なお、NLは、葉数を意味する。 The sum of the [Leaf Increase Index] of all leaves is taken as the [Total Leaf Increase Index], which is calculated using the following formula (12). NL stands for the number of leaves.

そして、これらの値を用い、葉の乾物重量である[Leaf DW](gDW/株)を次のように計算する。なお、[Leaf DW]の初期値は、第1葉(ここでは定植後に新たに展開してきた葉を第1葉とする)のみ生産者等が設定し(例えば10.0)、それ以外はゼロとする。
[Leaf DW]M,N=[Leaf DW]M,N-1+[Leaf ΔDW]M,N …(13A)
なお、[Leaf ΔDW]M,Nは、次式(13B)で表される。
[Leaf ΔDW]M,N
=[葉分配量]M,N・[Leaf Increase Index]M,N/[Total Leaf Increase Index]N
…(13B)
These values are then used to calculate the leaf dry weight, [Leaf DW] (gDW/plant), as follows: The initial value of [Leaf DW] is set by the producer (e.g., 10.0) only for the first leaf (here, the first leaf is defined as the leaf that newly unfolds after planting), and all other values are set to zero.
[Leaf DW] M,N = [Leaf DW] M,N-1 + [Leaf ΔDW] M,N … (13A)
[Leaf ΔDW] M,N is expressed by the following formula (13B).
[Leaf ΔDW] M,N
= [Leaf Allocation Amount] M,N * [Leaf Increase Index] M,N / [Total Leaf Increase Index] N
…(13B)

また、クラウンは、葉1枚増加と同時にクラウンが一定割合増加するとみなし、乾物重量である[Stem DW](gDW/株)を[Leaf DW]と同様、以下のように求める。なお、[Stem DW]の初期値は、第1節(ここでは定植後に新たに発生した節を第1節とする)のみ生産者等が設定し(例えば5.0)、それ以外はゼロとする。
[Stem DW]M,N=[Stem DW]M,N-1+[Stem ΔDW]M,N …(14A)
なお、[Stem ΔDW]M,Nは、次式(14B)で表される。
[Stem ΔDW]M,N
=[クラウン分配量]M,N・[Leaf Increase Index]M,N/[Total Leaf Increase Index]N
…(14B)
In addition, the crown is assumed to increase by a fixed percentage with each additional leaf, and the dry weight [Stem DW] (gDW/plant) is calculated in the same way as [Leaf DW], as follows: The initial value of [Stem DW] is set by the producer (e.g., 5.0) only for the first node (here, the node newly generated after planting is defined as the first node), and all other nodes are set to zero.
[Stem DW] M,N = [Stem DW] M,N-1 + [Stem ΔDW] M,N … (14A)
[Stem ΔDW] M,N is expressed by the following formula (14B).
[Stem ΔDW] M,N
= [Crown Distribution Amount] M,N * [Leaf Increase Index] M,N / [Total Leaf Increase Index] N
…(14B)

葉面積[Leaf Area](m2/株)は、比葉面積[SLA](m2/ gDW)を用い、次式(15)にて求める。なお、[Leaf Area]は、第1葉のみ生産者等が設定し、それ以外の葉は0とする。
[Leaf Area]M,N=[Leaf Area]M,N-1+[SLA]N・[葉分配量]M,N・[Leaf Increase Index]M,N/[Total Leaf Increase Index]N …(15)
Leaf area ( /plant) is calculated using specific leaf area (SLA) ( /gDW) using the following formula (15): Leaf area is set by the producer only for the first leaf, and other leaves are set to 0.
[Leaf Area] M,N = [Leaf Area] M,N-1 + [SLA] N・[Leaf distribution amount] M,N・[Leaf Increase Index] M,N /[Total Leaf Increase Index] N … (15)

また、作物に着生している葉の面積[Leaf Area]の総計を[Total Leaf Area(1)](m2/株)とすると、[Total Leaf Area(1)]は、次式(16)から求められる。 Furthermore, if the total area of leaves attached to a crop [Leaf Area] is [Total Leaf Area(1)] (m 2 /plant), [Total Leaf Area(1)] can be calculated using the following formula (16).

更に、光合成量計算に用いる葉のパラメータを計算する。 Furthermore, leaf parameters used to calculate photosynthesis are calculated.

まず、単位土地面積当たりの葉の面積である葉面積指数[LAI](無次元)については、土地面積当たりの作物株数である[Plant Density](株/m2)を用いて、次式(17)から求める。
[LAI]N=[Total Leaf Area]N・[Plant Density] …(17)
First, the leaf area index [LAI] (dimensionless), which is the leaf area per unit land area, is calculated from the following equation (17) using the plant density (plants/m 2 ), which is the number of crop plants per land area.
[LAI] N = [Total Leaf Area] N・[Plant Density] …(17)

また、土地面積当たりの日射量に対し、作物が受光する日射量の比を表す[相対受光量](無次元)は、吸光係数である[K](無次元)、および前日の葉面積指数を用いて、次式(18)から求める。なお、吸光係数は、群落内部への光の届きやすさを表わす係数である。
[相対受光量]N=1-exp(-[K]・[LAI]N) …(18)
The relative amount of received light (dimensionless), which represents the ratio of the amount of solar radiation received by the crop to the amount of solar radiation per unit land area, is calculated using the light absorption coefficient [K] (dimensionless) and the leaf area index of the previous day using the following equation (18). The light absorption coefficient is a coefficient that indicates the ease with which light reaches the interior of the canopy.
[Relative light receiving amount] N = 1 - exp(-[K]·[LAI] N ) ... (18)

(c)着果数計算
着果数は、作物の花芽分化が開始した翌日から開始する特定期間(花数決定期間)の光合成量(期間乾物生産量)、及び期間乾物生産量と着果数の関係式から、花房ごとに求められる。花芽分化は、花房順位の順番で発生し、1つ少ない花房順位(前順位花房)が分化してから、一定量の気温が蓄積しないと、次の花房順位が分化しないものとする。
(c) Calculation of fruit number The number of fruits is calculated for each inflorescence from the amount of photosynthesis (period dry matter production) during a specific period (flower number determination period) starting the day after flower bud differentiation of the crop begins, and the relational equation between period dry matter production and fruit number. Flower bud differentiation occurs in the order of inflorescence rank, and after the differentiation of the next inflorescence rank (previous inflorescence), a certain amount of temperature must accumulate before the next inflorescence rank will differentiate.

分化が開始したことを表す[分化開始Index]は、以下で説明する[分化条件判定A](℃)、および[分化条件判定B]により決定する。ここで、[花芽分化期積算温度](℃)は、[分化開始Index]が最初に1になった日の積算温度とし、それ以前は-1.0とする。 The [Differentiation Start Index], which indicates the start of differentiation, is determined by [Differentiation Condition Judgment A] (°C) and [Differentiation Condition Judgment B], which are explained below. Here, the [Flower Bud Differentiation Stage Accumulated Temperature] (°C) is the accumulated temperature on the day when the [Differentiation Start Index] first becomes 1, and before that it is -1.0.

本シミュレーションでは、上記[花芽分化期積算温度]の値をもとに、以下のようにして、前順位花房の花芽分化が開始した翌日からの積算温度が一定値(例えば600℃)を超えたかどうかを判定する指標である[分化条件判定A]を求める。なお、600℃という値は、花房と花房の間には当該品種では通常4枚の葉が発生するため、葉が1枚発生するのに必要な積算温度を150℃として設定している。式中のFは花房順位を示す。
[分化開始Index]F-1,N=1の場合、
[分化条件判定A]F,N
=[積算温度]N-([花芽分化期積算温度]F-1,N+600) …(19A)
その他の場合、
[分化条件判定A]F,N=-1.0 …(19B)
と求める。
In this simulation, based on the value of the [Flower Bud Differentiation Period Accumulated Temperature], [Differentiation Condition Judgment A] is calculated as follows, which is an index for determining whether the accumulated temperature from the day after flower bud differentiation of the pre-order inflorescence begins exceeds a certain value (for example, 600°C). Note that the value of 600°C is set at 150°C, the accumulated temperature required for one leaf to emerge, since in this variety, four leaves normally emerge between inflorescences. F in the formula indicates the inflorescence order.
[Differentiation Start Index] F-1,N = 1
[Differentiation condition judgment A] F,N
= [accumulated temperature] N - ([accumulated temperature at flower bud differentiation stage] F-1,N + 600) ... (19A)
In other cases,
[Differentiation condition judgment A] F,N = -1.0...(19B)
He demands.

また、光合成量が花芽分化可能な状態にあるかの判定指標である[分化条件判定B](gDW/株)を次式(20)から求める。
[分化条件判定B]N=([光合成量]のN-6~N日平均値) …(20)
In addition, [Differentiation Condition Judgment B] (gDW/strain), which is an index for judging whether the amount of photosynthesis is in a state where flower bud differentiation is possible, is calculated using the following formula (20).
[Differentiation condition determination B] N = (N-6 to N-day average value of [photosynthetic amount]) ... (20)

なお、花芽分化が開始したことを示す[分化開始Index]は、一旦1になった場合、それ以降は1になるようにする。
すなわち、[分化開始Index]F,N-1=1であれば、
[分化開始Index]F,N=1 …(21A)
0.0≦[分化条件判定A]F,N、かつ1.0<[分化条件判定B]Nであれば、
[分化開始Index]F,N=1 …(21B)
その他の場合には、
[分化開始Index]F,N=0 …(21C)
とする。
In addition, once the [Differentiation Start Index], which indicates the start of flower bud differentiation, becomes 1, it will remain 1 thereafter.
That is, if [differentiation start index] F,N-1 = 1,
[Differentiation start index] F,N = 1...(21A)
If 0.0≦[Differentiation Condition Judgment A] F,N and 1.0<[Differentiation Condition Judgment B] N , then
[Differentiation start index] F,N = 1...(21B)
In other cases,
[Differentiation start index] F,N = 0...(21C)
Let's say.

また、上式(19A)の[花芽分化期積算温度]は、[分化開始Index]の前日値が0、かつ当日値が1となった日の積算温度である。なお、[花芽分化期積算温度]の初期値は-1.0とし、一旦[花芽分化期積算温度]が設定された後は、同じ値とする。初期値を負値に設定するのは、花芽分化期開始以降の[花芽分化期積算温度]を正値、それより前を負値として区分するためである。この正負値による区分は、開花期積算温度でも同様に行う。
すなわち、[分化開始Index]F,N-1=0、かつ[分化開始Index]F,N=1であれば、
[花芽分化期積算温度]F,N=[積算温度]N …(22A)
その他の場合には、
[花芽分化期積算温度]F,N=[花芽分化期積算温度]F,N-1 …(22B)
である。
Furthermore, the [Flower bud differentiation stage accumulated temperature] in the above formula (19A) is the accumulated temperature on the day when the value of the [Differentiation start index] the previous day was 0 and the value of the day was 1. The initial value of the [Flower bud differentiation stage accumulated temperature] is -1.0, and once the [Flower bud differentiation stage accumulated temperature] is set, it remains the same value. The reason for setting the initial value to a negative value is to classify the [Flower bud differentiation stage accumulated temperature] after the start of the flower bud differentiation stage as a positive value and those before that as negative values. This classification by positive and negative values is also performed on the flowering stage accumulated temperature.
That is, if [differentiation start index] F,N-1 = 0 and [differentiation start index] F,N = 1,
[Accumulated temperature during flower bud differentiation] F,N = [Accumulated temperature] N ... (22A)
In other cases,
[Cumulative temperature during flower bud differentiation stage] F,N = [cumulative temperature during flower bud differentiation stage] F,N-1 ...(22B)
is.

更に、本シミュレーションでは、花芽分化開始後の積算温度が0~150℃の範囲にある期間を花数決定期間とし、花数決定期間であるかを示す指標となる[花数決定期間Index(1)](無次元)を以下のようにして求める。なお、花数決定期間は、葉が1枚展開する期間と定義し、1枚当たりの積算温度を150℃と仮定して、花数決定期間の積算温度を150℃としている。
すなわち、[分化開始Index]F,N=0であれば、
[花数決定期間Index(1)]F,N=-1.0 …(23A)
その他の場合には、
[花数決定期間Index(1)]F,N
=([積算温度]N-[花芽分化期積算温度]F,N)/150 …(23B)
Furthermore, in this simulation, the period after the start of flower bud differentiation in which the accumulated temperature is in the range of 0 to 150°C is defined as the flower number determination period, and the [flower number determination period Index(1)] (dimensionless), which serves as an index indicating whether or not the period is in the flower number determination period, is calculated as follows: The flower number determination period is defined as the period in which one leaf unfolds, and assuming that the accumulated temperature per leaf is 150°C, the accumulated temperature during the flower number determination period is set to 150°C.
That is, if [differentiation start index] F,N = 0,
[Flower number determination period Index (1)] F,N = -1.0...(23A)
In other cases,
[Flower number determination period Index(1)] F,N
= ([accumulated temperature] N - [accumulated temperature at flower bud differentiation stage] F,N ) / 150 ... (23B)

そして、[花数決定期間Index](無次元)は、花芽分化開始からの積算気温が0~150℃の範囲にある場合に1となるように、以下のようにして求められる。
0.0≦[花数決定期間Index(1)]F,N≦1.0であれば、
[花数決定期間Index]F,N=1 …(24A)
その他の場合には、
[花数決定期間Index]F,N=0 …(24B)
The [Flower Number Determination Period Index] (dimensionless) is calculated as follows so that it is 1 when the accumulated temperature from the start of flower bud differentiation is in the range of 0 to 150°C.
If 0.0≦[Flower number determination period Index(1)] F,N ≦1.0,
[Flower number determination period Index] F,N = 1...(24A)
In other cases,
[Flower number determination period Index] F,N = 0...(24B)

ところで、花数決定期間における光合成量の積算値である[決定期間光合成量](gDW/株)は、次式(25)にて表される。
[決定期間光合成量]F,N=[決定期間光合成量]F,N-1+[花数決定期間Index]F,N・[光合成量]F,N …(25)
The integrated value of the amount of photosynthesis during the flower number determination period, [determination period photosynthesis amount] (gDW/plant), is expressed by the following formula (25).
[Photosynthesis amount during determination period] F,N = [Photosynthesis amount during determination period] F,N-1 + [Flower number determination period Index] F,N・[Photosynthesis amount] F,N … (25)

そして、[着果数(1)](個/株)は、上式(25)の[決定期間光合成量]、[決定期間乾物生産量当たりの着果数](m2・個/gDW・株)、及び[Plant Density](株/m2)を用いて、次式(26)のように求めることができる。
[着果数(1)]F,N=[決定期間光合成量]M,N・[決定期間乾物生産量当たりの着果数]・[Plant Density] …(26)
Then, the number of fruits (1) (pieces/plant) can be calculated using the [Photosynthetic amount during the determination period], [Number of fruits per dry matter production during the determination period] ( m2 ·pieces/gDW·plant), and [Plant Density] (plant/m2) from the above formula ( 25 ) as shown in the following formula (26).
[Number of fruits (1)] F,N = [Photosynthetic amount in the critical period] M,N × [Number of fruits per dry matter production in the critical period] × [Plant Density] ... (26)

なお、[着果数](個/株)は、以下のように、1.0~10.0の範囲となるように補正する。
すなわち、[着果数(1)]F,N<1.0であれば、
[着果数]F,N=1.0 …(27A)
1.0≦[着果数(1)]F,N≦10.0であれば、
[着果数]F,N=[着果数(1)]F,N …(27B)
その他の場合には、
[着果数]F,N=10.0 …(27C)
とする。
The number of fruits (pieces/plant) is corrected to fall within the range of 1.0 to 10.0 as follows:
That is, if [fruit number (1)] F,N < 1.0,
[Number of fruit set] F,N = 1.0…(27A)
If 1.0≦[Number of fruits (1)] F,N ≦10.0,
[Number of fruit set] F,N= [Number of fruit set (1)] F,N … (27B)
In other cases,
[Number of fruit set] F,N = 10.0…(27C)
Let's say.

なお、開花期については、花芽分化期からの積算温度が600℃に達した時点とし、その時点での積算温度を[開花期積算温度](℃)としている。なお、600℃という値は、花芽が分化してから葉が4枚で開花し、1枚当たり150℃の積算温度が必要であると仮定して、設定している。なお、開花期積算温度は、花芽分化期積算温度が決定された後にのみ有効になるようにし、初期値を-1.0とする。また、開花期積算温度が一旦設定された後は、同じ値であるものとする。
すなわち、[花芽分化期積算温度]F,N<0であれば、
[開花期積算温度]F,N=-1.0 …(28A)
その他の場合には、
[開花期積算温度]F,N=[花芽分化期積算温度]F,N+600 …(28B)
となる。
The flowering stage is defined as the point at which the accumulated temperature from the flower bud differentiation stage reaches 600°C, and the accumulated temperature at that point is defined as the [flowering stage accumulated temperature] (°C). The value of 600°C is set on the assumption that four leaves bloom after flower bud differentiation, and that an accumulated temperature of 150°C is required per leaf. The flowering stage accumulated temperature is set to be valid only after the flower bud differentiation stage accumulated temperature has been determined, and its initial value is set to -1.0. Once the flowering stage accumulated temperature is set, it remains the same value.
That is, if [flower bud differentiation stage accumulated temperature] F,N < 0,
[Flowering stage accumulated temperature] F,N = -1.0 ... (28A)
In other cases,
[Accumulated temperature at flowering stage] F,N = [Accumulated temperature at flower bud differentiation stage] F,N + 600 ... (28B)
This becomes:

(d)花房乾物重量計算
花房の相対的生育段階である[Fruits Growth Index](無次元)を以下のようにして求める。開花期積算温度が無効な場合(0未満の場合)は、Fruits Growth Indexの値は-1.0とする。また、Fruits Growth Indexの初期値も-1.0とする。
すなわち、[開花期積算温度]F,N<0であれば、
[Fruits Growth Index]F,N=-1.0 …(29A)
その他の場合には、
[Fruits Growth Index]F,N
=([積算温度]N-[開花期積算温度]F,N)/600 …(29B)
とする。
(d) Calculation of inflorescence dry weight The Fruits Growth Index (dimensionless), which is the relative growth stage of the inflorescence, is calculated as follows. If the flowering stage accumulated temperature is invalid (less than 0), the value of the Fruits Growth Index is set to -1.0. The initial value of the Fruits Growth Index is also set to -1.0.
That is, if [Flowering stage accumulated temperature] F,N < 0,
[Fruits Growth Index] F,N =-1.0...(29A)
In other cases,
[Fruits Growth Index] F,N
= ([accumulated temperature] N - [accumulated temperature at flowering stage] F,N ) / 600 ... (29B)
Let's say.

また、 [Fruits Size Index](無次元)を次式(30)により計算する。Fruits Size Indexの初期値は、ゼロとする。 Fruits Size Index (dimensionless) is calculated using the following equation (30). The initial value of Fruits Size Index is set to zero.

また、前日24時(当日0時)から当日24時における花房の相対的な変化の大きさから、花房の相対的な成長増加量として[Fruits Increase Index(1)](無次元)を次式(31)から求める。
[Fruits Increase Index(1)]F,N
=[Fruits Size Index]F,N-[Fruits Size Index]F,N-1 …(31)
In addition, the relative growth increase of the inflorescence is calculated from the magnitude of the relative change in the inflorescence from midnight the previous day (midnight on the current day) to midnight on the current day, and the [Fruits Increase Index(1)] (dimensionless) is calculated using the following equation (31).
[Fruits Increase Index(1)] F,N
=[Fruits Size Index] F,N -[Fruits Size Index] F,N-1 …(31)

ただし、花房は、着果数が10であることを前提にしているため、次式(32)のように補正し、[Fruits Increase Index](無次元)とする。
[Fruits Increase Index]F,N
=([Fruits Increase Index(1)]F,N・[着果数]F,N)/10 …(32)
However, since the inflorescence is based on the assumption that the number of fruits is 10, the index is corrected as shown in the following equation (32) to obtain the [Fruits Increase Index] (dimensionless).
[Fruits Increase Index] F,N
=([Fruits Increase Index(1)] F,N・[Number of fruit set] F,N )/10…(32)

また、全花房の[Fruits Increase Index]の総計を[Total Fruits Increase Index](無次元)とする(次式(33)参照)。なお、NFは、花房数を意味する。 The sum of the [Fruits Increase Index] of all inflorescences is defined as the [Total Fruits Increase Index] (dimensionless) (see the following equation (33)). NF stands for the number of inflorescences.

本シミュレーションでは、これらの値を用いて花房乾物重量である[Fruits DW](gDW/株)を求める。なお、基本的には、[Fruits DW]は、[Fruits Growth Index]が1より大きくなった後は、増加しない。また、[Fruits DW]の初期値は、ゼロとする。
[Fruits DW]F,N=[Fruits DW]F,N-1+[Fruits ΔDW]F,N …(34A)
なお、[Fruits ΔDW]F,Nは、次式(34B)で表される。
[Fruits ΔDW]F,N
=[花房分配量]N・[Fruits Increase Index]F,N/[Total Fruits Increase Index]N
…(34B)
In this simulation, these values are used to calculate the inflorescence dry weight, [Fruits DW] (gDW/plant). Basically, [Fruits DW] does not increase once [Fruits Growth Index] exceeds 1. The initial value of [Fruits DW] is set to zero.
[Fruits DW] F,N = [Fruits DW] F,N-1 + [Fruits ΔDW] F,N … (34A)
[Fruits ΔDW] F,N is expressed by the following equation (34B).
[Fruits ΔDW] F,N
=[Flower cluster distribution amount] N・[Fruits Increase Index] F,N /[Total Fruits Increase Index] N
…(34B)

(e)分配率計算
作物の[分配率](無次元)は、作物成長量をもとに作物全体の値として計算される。分配率は、花房におけるソースの受け取りやすさである[花房シンク強度](無次元)、および全器官のシンク強度の合計である[シンク強度合計]から求められる。
(e) Calculation of the distribution ratio The distribution ratio (dimensionless) of a crop is calculated as a value for the whole crop based on the amount of crop growth. The distribution ratio is calculated from the inflorescence sink strength (dimensionless), which is the ease with which the inflorescence receives the source, and the total sink strength, which is the sum of the sink strengths of all organs.

解体調査時などの実験データから[Total Fruits Increase Index]、[Total Leaf Increase Index]、[Total Stem Increase Index]の重みづけを行う。[果実分配調整係数]、[葉分配調整係数]、[クラウン分配調整係数]を[α]、[β]、[γ]とした。これらから、花房シンク強度、葉シンク強度、クラウンシンク強度、シンク強度合計を次式から求める。
[花房シンク強度]N=[Total Fruits Increase Index]N・[α] …(35A)
[葉シンク強度]N=[Total Leaf Increase Index]N・[β] …(35B)
[クラウンシンク強度]N=[Total Stem Increase Index]N・[γ] …(35C)
[シンク強度合計]N=[花房シンク強度]N+[葉シンク強度]N+[クラウンシンク強度]N +[根シンク強度] …(35D)
The [Total Fruits Increase Index], [Total Leaf Increase Index], and [Total Stem Increase Index] are weighted based on experimental data from dissection surveys, etc. The [Fruit Distribution Adjustment Coefficient], [Leaf Distribution Adjustment Coefficient], and [Crown Distribution Adjustment Coefficient] are designated as [α], [β], and [γ]. From these, the inflorescence sink strength, leaf sink strength, crown sink strength, and total sink strength are calculated using the following formula.
[Flower cluster sink strength] N = [Total Fruits Increase Index] N · [α] ... (35A)
[Leaf sink strength] N = [Total Leaf Increase Index] N · [β] ... (35B)
[Crown sink strength] N = [Total Stem Increase Index] N · [γ] ... (35C)
[Total sink strength] N = [Flower cluster sink strength] N + [Leaf sink strength] N + [Crown sink strength] N + [Root sink strength] ... (35D)

[花房分配率]、[根分配率]、[葉分配率]、[クラウン分配率](無次元)を次式(36A)~(36D)から求める。
[花房分配率]N=[花房シンク強度]N/[シンク強度合計]N …(36A)
[根分配率]N=0.05 …(36B)
[葉分配率]N=[葉シンク強度]N/[シンク強度合計]N…(36C)
[クラウン分配率]N=[クラウンシンク強度]N/[シンク強度合計]N …(36D)
The [inflorescence distribution rate], [root distribution rate], [leaf distribution rate], and [crown distribution rate] (dimensionless) are calculated from the following equations (36A) to (36D).
[Flower cluster distribution rate] N = [Flower cluster sink strength] N / [Total sink strength] N ... (36A)
[Root distribution ratio] N = 0.05...(36B)
[Leaf distribution ratio] N = [Leaf sink strength] N / [Total sink strength] N ... (36C)
[Crown distribution ratio] N = [Crown sink strength] N / [Total sink strength] N ... (36D)

(4)果実収穫量計算について
本シミュレーションでは、対象の農作物は、花房内において1日で収穫せず、徐々に収穫していくものであることを想定している。収穫が可能であることを示す[収穫開始Index](無次元)は、[Fruits Increase Index]が初めて1.0以上になった翌日以降に1とし、それ以前はゼロとする。初めて[収穫開始Index]が1となった日を収穫開始日とする。[収穫開始Index]は、以下のように表され、初期値はゼロである。
すなわち、[収穫開始Index]F,N-1=1であれば、
[収穫開始Index]F,N=1 …(37A)
[Fruits Source Index]F,N-1≧1.0であれば、
[収穫開始Index]F,N=1 …(37B)
その他の場合には、
[収穫開始Index]F,N=0 …(37C)
となる。
(4) Fruit Yield Calculation This simulation assumes that the target crops are not harvested in one day within the inflorescence, but are harvested gradually. The [Harvest Start Index] (dimensionless), which indicates that harvesting is possible, is set to 1 from the day after the [Fruits Increase Index] first reaches 1.0 or higher, and is set to zero before that. The first day the [Harvest Start Index] reaches 1 is considered the harvest start date. The [Harvest Start Index] is expressed as follows, and its initial value is zero.
That is, if [Harvesting Start Index] F,N-1 = 1,
[Harvest start index] F,N = 1...(37A)
[Fruits Source Index] If F,N-1 ≧ 1.0,
[Harvest start index] F,N = 1 ... (37B)
In other cases,
[Harvest start index] F,N = 0…(37C)
This becomes:

ここで、果実の花房ごとの総収穫量は、収穫開始日における[Fruits DW] (gDW/株)である。なお、基本的には、[Fruits DW]は、[Fruits Growth Index]が1より大きくなった後は、増加しないので、収穫開始日以降であれば、同じである。果実1個に対し、積算温度60℃になる期間をかけて収穫するものとし、以下のようにして、気温1℃当たりに収穫する果実の乾物重量である[収穫果実気温DW](gDW/株・℃)を計算する。なお、[収穫果実気温DW]の初期値は-1.0とし、有効な値が一旦設定された後は、同じ値とする。
すなわち、[収穫果実気温DW]F,N-1≧0.0であれば、
[収穫果実気温DW]F,N=[収穫果実気温DW]F,N-1 …(38A)
[収穫開始Index]F,N=1、かつ[着果数]F,N-1>0.0であれば、
[収穫果実気温DW]F,N=[Fruits DW]F,N-1/([着果数]F,N-1・60) …(38B)
その他の場合には、
[収穫果実気温DW]F,N=-1.0 …(38C)
と計算する。
Here, the total harvest volume per fruit inflorescence is [Fruits DW] (gDW/plant) on the start date of harvest. Basically, [Fruits DW] does not increase once the [Fruits Growth Index] exceeds 1, so it remains the same from the start date of harvest. Each fruit is harvested over a period of time until the accumulated temperature reaches 60°C. The dry weight of fruit harvested per 1°C of air temperature, [Harvested Fruit Temperature DW] (gDW/plant·°C), is calculated as follows. The initial value of [Harvested Fruit Temperature DW] is -1.0, and once a valid value is set, it remains the same.
That is, if [Harvested fruit temperature DW] F,N-1 ≧ 0.0,
[Harvested fruit temperature DW] F,N = [Harvested fruit temperature DW] F,N-1 ... (38A)
If [Harvesting Start Index] F,N = 1 and [Number of Fruits] F,N-1 > 0.0,
[Harvested Fruit Temperature DW] F,N = [Fruits DW] F,N-1 / ([Number of Fruits] F,N-1 x 60) ... (38B)
In other cases,
[Harvested fruit temperature DW] F,N = -1.0 ... (38C)
Calculate as follows.

そして、計算日における収穫作業前の収穫可能な果実乾物重量を[Fruits unHarvest](gDW/株)とし、[Fruits unHarvest]を以下のようにして求める。なお、[Fruits unHarvest]は、収穫開始日に前日の[Fruits DW]が与えられ、それ以降は収穫された分が減少していく。また、[Fruits unHarvest]の初期値は-1.0とし、収穫期間前後の値も-1.0である。
すなわち、[Fruits unHarvest]F,N-1≧0.0であれば、
[Fruits unHarvest]F,N
=[Fruits unHarvest]F,N-1-[Fruits HARVEST]F,N-1 …(39A)
[収穫開始Index]F,N-1=0、かつ[収穫開始Index]F,N=1であれば、
[Fruits unHarvest]F,N=[Fruits DW]F,N-1 …(39B)
その他の場合には、
[Fruits unHarvest]F,N=-1.0 …(39C)
となる。
The dry weight of harvestable fruit before harvesting on the calculation date is defined as [Fruits unHarvest] (gDW/plant), and [Fruits unHarvest] is calculated as follows: [Fruits unHarvest] is calculated by taking the [Fruits DW] of the previous day on the day harvesting begins, and then decreasing by the amount harvested thereafter. The initial value of [Fruits unHarvest] is -1.0, and the values before and after the harvest period are also -1.0.
That is, if [Fruits unHarvest] F,N-1 ≧0.0, then
[Fruits unHarvest] F,N
=[Fruits unHarvest] F,N-1 -[Fruits HARVEST] F,N-1 …(39A)
If [Harvesting Start Index] F,N-1 = 0 and [Harvesting Start Index] F,N = 1,
[Fruits unHarvest] F,N = [Fruits DW] F,N-1 …(39B)
In other cases,
[Fruits unHarvest] F,N =-1.0...(39C)
This becomes:

なお、収穫期間開始後、かつ収穫可能な果実がある場合には、以下のようにして、収穫された果実乾物重量である[Fruits HARVEST](gDW/株)が計算される。
すなわち、[収穫開始Index]F,N=1、かつ[Fruits unHarvest]F,N>0.0の場合、
[Fruits HARVEST]F,N
=min([収穫果実気温DW]F・[気温]N,[Fruits unHarvest]F,N) …(40A)
その他の場合には、
[Fruits HARVEST]F,N=0.0 …(40B)
After the start of the harvest period, if there are fruits available for harvest, the dry weight of the harvested fruits [Fruits HARVEST] (gDW/plant) is calculated as follows.
That is, if [Harvest Start Index] F,N = 1 and [Fruits unHarvest] F,N > 0.0,
[Fruits HARVEST] F,N
= min([Fruits unHarvest] F , N ) ... (40A)
In other cases,
[Fruits HARVEST] F,N =0.0…(40B)

また、全花房の[Fruits HARVEST]総計を[果実収穫乾物重量](gDW/株)とすると、[果実収穫乾物重量]は、次式(41)にて表すことができる。 Furthermore, if the total [Fruits Harvest] of all inflorescences is [Fruits Harvest Dry Weight] (gDW/plant), [Fruits Harvest Dry Weight] can be expressed by the following equation (41):

この[果実収穫乾物重量]に[果実生体重](g/gDW)を乗じることで、次式(42)のように、[果実収穫生体重量](g/株)が求められる。
[果実収穫生体重量]N=[果実収穫乾物重量]N・[果実生体重] …(42)
By multiplying this [fruit harvest dry weight] by the [fruit fresh weight] (g/gDW), the [fruit harvest fresh weight] (g/plant) can be calculated as shown in the following equation (42).
[Fruit harvest fresh weight] N = [Fruit harvest dry weight] N · [Fruit fresh weight] ... (42)

そして、[果実収穫生体重量]に[圃場株数](株)を乗じることで、次式(43)のように、圃場全体の収量である[圃場果実収穫生体重量](kg)が求められる。
[圃場果実収穫生体重量]N
=0.001・[果実収穫生体重量]N・[圃場株数] …(43)
Then, by multiplying the [live weight of fruit harvested] by the [number of plants in the field] (plants), the [live weight of fruit harvested in the field] (kg), which is the yield of the entire field, can be calculated as shown in the following equation (43).
[Field fruit harvested live weight] N
= 0.001 · [fruit harvested live weight] N · [number of plants in the field] ... (43)

シミュレーション部44は、環境データやパラメータを用いて上述したような生育モデルを作成し、当該生育モデルを用いて、生産者が選択した品種、栽培地点、栽培情報の組み合わせごとに生育に関するシミュレーションを実行する。そして、各組み合わせのシミュレーション結果(シミュレーションにより得られる値)を出力部45に送信する。なお、シミュレーション部44は、品目や品種に応じたパラメータを用いることで、イチゴ以外、例えばトマト、イチゴ、キュウリ、パプリカなどの果菜類のシミュレーションを行うことも可能である。 The simulation unit 44 creates a growth model as described above using environmental data and parameters, and uses the growth model to perform a growth simulation for each combination of variety, cultivation location, and cultivation information selected by the producer. The simulation results for each combination (values obtained by the simulation) are then sent to the output unit 45. Note that the simulation unit 44 can also simulate fruit and vegetable crops other than strawberries, such as tomatoes, strawberries, cucumbers, and peppers, by using parameters according to the item and variety.

(5)養分吸収量計算
葉、クラウン、果実、根の乾物増加量に養分含有率を乗じて、養分吸収量を求める。ここではN(窒素)を例として記載するが、他の元素でも同様な計算で求めることができる。
[Leaf ΔN]M,N=[Leaf ΔDW]M,N×[Leaf N%]M,N …(44A)
[Stem ΔN]M,N=[Stem ΔDW]M,N×[Stem N%]M,N …(44B)
[Fruits ΔN]F,N=[Fruits DW]F,N×[Fruits N%]F,N …(44C)
[Root ΔN]M,N=[Root DW]M,N×[Root N%]M,N …(44D)
(5) Calculation of nutrient absorption The amount of nutrient absorption is calculated by multiplying the increase in dry matter of leaves, crowns, fruits, and roots by the nutrient content. Here, N (nitrogen) is used as an example, but similar calculations can be used for other elements.
[Leaf ΔN] M,N = [Leaf ΔDW] M,N × [Leaf N%] M,N … (44A)
[Stem ΔN] M,N = [Stem ΔDW] M,N × [Stem N%] M,N … (44B)
[Fruits ΔN] F,N = [Fruits DW] F,N × [Fruits N%] F,N … (44C)
[Root ΔN] M,N = [Root DW] M,N × [Root N%] M,N … (44D)

なお、上式(44A)の[Leaf ΔN]M,Nが葉の養分吸収量であり、[Leaf ΔDW]M,Nが葉の乾物増加量であり、[Leaf N%]M,Nが葉の養分含有率である。その他の式(44B)~(44D)も同様である。 In the above formula (44A), [Leaf ΔN] M,N is the amount of nutrient absorbed by the leaf, [Leaf ΔDW] M,N is the increase in dry matter in the leaf, and [Leaf N%] M,N is the nutrient content in the leaf. The same applies to the other formulas (44B) to (44D).

ここで、養分含有率([Leaf N%]M,N等)はその器官の生育段階(Growth index)の関数として、例えば以下の式で表される。なお、β、γは定数である。
[Leaf N%]M,N=βL×ln([Leaf Growth Index]M,N)+γL …(45A)
[Stem N%]M,N=βS×ln([Stem Growth Index]M,N)+γS …(45B)
[Fruits N%]F,N=βF×ln([Fruits Growth Index]F,N)+γF …(45C)
[Root N%]M,N=γR …(45D)
Here, the nutrient content ( M, N , etc. [Leaf N%]) is expressed as a function of the growth index of the organ, for example, by the following formula: where β and γ are constants.
[Leaf N%] M,N = β L × ln ([Leaf Growth Index] M,N ) + γ L … (45A)
[Stem N%] M,N = β S × ln ([Stem Growth Index] M,N ) + γ S … (45B)
[Fruits N%] F,N = β F × ln ([Fruits Growth Index] F,N ) + γ F … (45C)
[Root N%] M,N = γ R … (45D)

また、葉、クラウン、花房、根の各養分吸収量は、以下の式で表される。 The amount of nutrients absorbed by leaves, crowns, inflorescences, and roots can be expressed by the following formula:

更に、株あたりの養分吸収量[Total ΔN]Nは以下の式(47)にて表される。
[Total ΔN]N
=[Total Leaf ΔN]N+[Total Stem ΔN]N+[Total Fruits ΔN]N+[Total Root ΔN]N
…(47)
Furthermore, the nutrient absorption amount per plant [Total ΔN] N is expressed by the following equation (47).
[Total ΔN] N
=[Total Leaf ΔN] N +[Total Stem ΔN] N +[Total Fruits ΔN] N +[Total Root ΔN] N
…(47)

なお、[Total ΔN]Nは株あたりの乾物重の増加量から推定した「養分吸収量」である。これに栽植密度を乗じて「面積当たりの養分吸収量」とすることができる。更に、「養分吸収量」を[肥料利用効率]で除することによって、「施肥量」とすることもできる。 [Total ΔN] N is the "nutrient absorption amount" estimated from the increase in dry weight per plant. This can be multiplied by the planting density to obtain the "nutrient absorption amount per area." Furthermore, the "nutrient absorption amount" can be calculated by dividing it by the [fertilizer use efficiency] to obtain the "fertilizer amount."

(出力部45の処理について)
出力部45は、シミュレーション部44から送信されてくるシミュレーション結果を受け付けると、シミュレーション結果を表示する画面を生成する。例えば、出力部45は、シミュレーション結果として、図11に示すような情報を表示する画面を生成する。この場合、利用者端末70の表示部193上には、図11の画面が表示されるので、生産者は、どの品種をどの栽培地点でどのように栽培すると、どのように生育するのか、どの程度の収量が得られるのかを確認することができる。また、複数の品目を異なる栽培条件で栽培したときのシミュレーション結果を並べて表示するため、生産者は、品種や栽培条件を異ならせた場合にどのような栽培結果が得られるのかを比較することができる。
(Regarding processing by the output unit 45)
When the output unit 45 receives the simulation results transmitted from the simulation unit 44, it generates a screen displaying the simulation results. For example, the output unit 45 generates a screen displaying information such as that shown in FIG. 11 as the simulation results. In this case, the screen shown in FIG. 11 is displayed on the display unit 193 of the user terminal 70, allowing the producer to confirm which variety, at which cultivation location, and how it will grow, as well as the yield that can be obtained. In addition, because the simulation results for cultivating multiple items under different cultivation conditions are displayed side by side, the producer can compare the cultivation results that can be obtained when the variety and cultivation conditions are changed.

また、出力部45は、例えば、生産者が選択した品種を選択した栽培地点で栽培したときの花房別出荷量を示す、図12に示すような画面を生成し、利用者端末70に出力することもできる。生産者は、図12の画面を参照することで、頂花房、第2花房、第3花房…における収穫量の推移を確認することができる。なお、図12においては、各花房の収穫量を別々のグラフに表示しているが、各グラフを1つのグラフにまとめて表示することとしてもよい。この場合、各花房の収穫量を示す棒グラフを色分け表示等して、わかりやすく表示してもよい。更に、出力部45は、例えば図13(a)に示すように、総収量の推移をグラフ表示してもよいし、図13(b)に示すように、LAIの推移をグラフ表示してもよい。また、その他の生育に関する情報(例えば、葉面積、花房別の開花日、光合成量、生育量(クラウン、葉、果実、着果負担)、養分吸収量(施肥量)など)を数値表示したり、グラフ表示してもよい。 The output unit 45 can also generate a screen such as that shown in FIG. 12, which shows the shipping volume by inflorescence when the variety selected by the producer is cultivated at the selected cultivation site, and output the screen to the user terminal 70. By referring to the screen in FIG. 12, the producer can check the trends in the yield for the terminal inflorescence, the second inflorescence, the third inflorescence, etc. While FIG. 12 shows the yield for each inflorescence in separate graphs, these graphs may be displayed together in a single graph. In this case, the bar graphs showing the yield for each inflorescence may be displayed in different colors for easier understanding. Furthermore, the output unit 45 may graphically display the trends in total yield, as shown in FIG. 13(a), or the trends in LAI, as shown in FIG. 13(b). Other growth-related information (e.g., leaf area, flowering date by inflorescence, photosynthesis rate, growth rate (crown, leaf, fruit, fruit load), nutrient absorption rate (fertilizer amount), etc.) may also be displayed numerically or graphically.

また、出力部45は、例えば図14に示すように、ある栽培地点で品種A、B、Cを栽培したときの、葉、クラウン、各花房の乾物重の推移のシミュレーション結果を表示することとしてもよい。 The output unit 45 may also display the simulation results of the changes in dry weight of leaves, crowns, and inflorescences when cultivars A, B, and C are cultivated at a certain cultivation location, as shown in Figure 14, for example.

以上のように、出力部45は、図11に示すように、品種、栽培地点、栽培情報の組み合わせごとにシミュレーション結果を比較可能に表示してもよいし、図12~図13(b)に示すように、ある品種をある栽培地点で栽培したときのシミュレーション結果を表示してもよい。また、図14に示すように、同一の栽培地点で異なる品種を栽培した場合のシミュレーション結果を比較可能に表示してもよい。更に、異なる栽培地点である品種を栽培した場合のシミュレーション結果を比較可能に表示してもよい。いずれにしても、従来のように品種の特徴を言葉や画像、平均値などで表現したカタログと比べ、品種の特徴を適切に表現することが可能となっている。したがって、生産者は、シミュレーション結果を確認することで、品種や栽培地点の選定を適切に行うことができる。また、将来において農作物がどのように生育するのかを把握することができるため、栽培管理(作業者の確保、資材等の調達など)を適切に行うことができる。 As described above, the output unit 45 may display simulation results for each combination of variety, cultivation location, and cultivation information in a comparable manner, as shown in FIG. 11, or may display simulation results for cultivating a certain variety at a certain cultivation location, as shown in FIGS. 12 to 13(b). Furthermore, as shown in FIG. 14, simulation results for cultivating different varieties at the same cultivation location may be displayed in a comparable manner. Furthermore, simulation results for cultivating a certain variety at different cultivation locations may be displayed in a comparable manner. In either case, it is possible to appropriately express the characteristics of a variety compared to conventional catalogs that express the characteristics of a variety using words, images, average values, etc. Therefore, by checking the simulation results, producers can appropriately select varieties and cultivation locations. Furthermore, by being able to understand how crops will grow in the future, they can appropriately manage cultivation (securing workers, procuring materials, etc.).

(更新部46の処理について)
更新部46は、生産者から、収量の推移やLAIの推移などの実測値(実際の栽培結果)が入力された場合に、当該実測値と、対応するシミュレーション結果とを比較し、シミュレーション結果が実測値に近づくように、各種パラメータを調整し、更新する。このようにすることで、実測値に基づいてパラメータが適切に更新されるため、次回以降のシミュレーションにおいて、適切なシミュレーション結果を得ることができるようになる。
(Regarding the processing of the update unit 46)
When the producer inputs actual measurements (actual cultivation results) such as changes in yield or LAI, the update unit 46 compares the actual measurements with the corresponding simulation results, and adjusts and updates various parameters so that the simulation results approach the actual measurements. In this way, the parameters are appropriately updated based on the actual measurements, making it possible to obtain appropriate simulation results in the next and subsequent simulations.

なお、全ての生産者から入力される情報に基づいて、パラメータを更新すると、パラメータが適切に更新されないおそれがある。このため、更新部46は、予め定めた生産者(信頼度の高い生産者)から情報が入力されたときにのみ、入力された情報に基づいてパラメータを更新するようにしてもよい。この場合、更新部46は、図15に示すようなフローチャートに沿った処理を実行する。 Note that if the parameters are updated based on information input from all producers, there is a risk that the parameters may not be updated appropriately. For this reason, the update unit 46 may update the parameters based on the input information only when information is input from a predetermined producer (a producer with a high level of reliability). In this case, the update unit 46 executes processing in accordance with the flowchart shown in FIG. 15.

図15の処理では、まず、ステップS30において、更新部46は、生産者が利用者端末70を介して、実際の栽培結果(実測値)を入力するまで待機する。そして、実際の栽培結果が入力されると、ステップS32に移行し、更新部46は、予め定められた生産者(信頼度の高い生産者)からの入力であったか否かを判断する。このステップS32の判断が否定された場合には、そのまま図15の全処理を終了するが、肯定された場合には、ステップS34に移行する。ステップS34に移行した場合、更新部46は、シミュレーション結果が、実際の栽培結果に近づくようにパラメータを更新し、パラメータDB50に格納し、図15の全処理を終了する。 In the processing of FIG. 15, first, in step S30, the update unit 46 waits until the producer inputs actual cultivation results (actual measurements) via the user terminal 70. Then, when the actual cultivation results are input, the processing proceeds to step S32, where the update unit 46 determines whether the input was from a predetermined producer (a producer with a high level of reliability). If the determination in step S32 is negative, the entire processing of FIG. 15 ends, but if the determination is positive, the processing proceeds to step S34. If the processing proceeds to step S34, the update unit 46 updates the parameters so that the simulation results approach the actual cultivation results, stores the updated parameters in the parameter DB 50, and then the entire processing of FIG. 15 ends.

なお、更新部46は、ある生産者から入力された情報を用いてパラメータを更新した場合には、更新後のパラメータを当該生産者専用のパラメータとして管理するようにしてもよい。これにより、各生産者が、実際の栽培結果に合わせて、パラメータをカスタマイズすることができるようになる。なお、ある生産者が別の生産者のパラメータを呼び出して利用できるようにしてもよい。 When the update unit 46 updates parameters using information input by a producer, the updated parameters may be managed as parameters exclusive to that producer. This allows each producer to customize parameters to suit their actual cultivation results. It may also be possible for one producer to call up and use the parameters of another producer.

なお、生産者は、パラメータを直接修正してもよい。生産者がパラメータを修正した場合には、更新部46は、パラメータDB50を修正内容に応じて更新するようにすればよい。 The producer may also modify the parameters directly. If the producer modifies the parameters, the update unit 46 may update the parameter DB 50 in accordance with the modifications.

これまでの説明から分かるように、本実施形態のシミュレーション部44は、農作物の品種ごとの特徴を示すパラメータを記憶部(パラメータDB50)から読み出す読み出し部として機能する。また、本実施形態のシミュレーション部44は、読み出したパラメータと、栽培環境に関する情報(環境データ)と、に基づいて、品種ごとの生育モデルを作成し、品種ごとの生育モデルから、品種ごとの生育に関する予測を実行する予測部として機能する。 As can be seen from the explanation so far, the simulation unit 44 of this embodiment functions as a readout unit that reads out parameters indicating the characteristics of each crop variety from the storage unit (parameter DB 50). The simulation unit 44 of this embodiment also functions as a prediction unit that creates a growth model for each variety based on the readout parameters and information about the cultivation environment (environmental data), and performs predictions regarding the growth of each variety from the growth model for each variety.

以上、詳細に説明したように、本実施形態によると、シミュレーション部44は、農作物の品種ごとの特徴を示すパラメータをパラメータDB50から読み出し、環境データ取得部42から、生産者が選択した栽培地点に対応する環境データを取得する。また、シミュレーション部44は、取得した環境データと、読み出したパラメータと、に基づいて、品種ごとの生育モデルを作成し、作成した生育モデルを用いてシミュレーションを実行する。そして、出力部45は、シミュレーション結果を品種選定や栽培管理を支援する情報として出力する。このように、本実施形態では、品種ごと及び栽培地点ごとに作成した生育モデルを用いて農作物の生育に関するシミュレーションを行った結果を品種選定や栽培管理を支援する情報として表示するため、従来のように品種の特徴を言葉や画像、平均値などで表現したカタログと比べ、品種の特徴を適切に表現することが可能である。したがって、生産者は、品種や栽培地点の選定を適切に行うことができるとともに、将来においてどのように生育するのかを把握することができるため、栽培管理を適切に行うことができる。 As described above in detail, in this embodiment, the simulation unit 44 reads parameters indicating the characteristics of each crop variety from the parameter DB 50 and obtains environmental data corresponding to the cultivation location selected by the producer from the environmental data acquisition unit 42. The simulation unit 44 also creates a growth model for each variety based on the acquired environmental data and the read parameters, and performs a simulation using the created growth model. The output unit 45 then outputs the simulation results as information to support variety selection and cultivation management. In this way, in this embodiment, the results of a simulation of crop growth using a growth model created for each variety and cultivation location are displayed as information to support variety selection and cultivation management. This makes it possible to appropriately express the characteristics of a variety compared to conventional catalogs that express variety characteristics using words, images, average values, etc. This allows producers to appropriately select varieties and cultivation locations, and to understand how the crops will grow in the future, thereby enabling appropriate cultivation management.

また、本実施形態では、更新部46は、生産者が入力した実際の栽培結果(収量等)に基づいて、パラメータを更新して、パラメータDB50を更新する。これにより、栽培実績に基づいてパラメータを適切な値に更新することができる。この場合、更新部46は、シミュレーション結果が栽培実績に近づくように、パラメータを更新することで、シミュレーションの精度を向上することができる。 In addition, in this embodiment, the update unit 46 updates the parameters based on the actual cultivation results (yield, etc.) entered by the producer, and updates the parameter DB 50. This allows the parameters to be updated to appropriate values based on the cultivation results. In this case, the update unit 46 can improve the accuracy of the simulation by updating the parameters so that the simulation results approach the cultivation results.

また、本実施形態では、更新部46は、生産者が予め定められた生産者(信頼度が所定以上の生産者)である場合(S32:肯定)に、当該生産者が入力した栽培実績に基づいてパラメータを更新する(S34)。これにより、パラメータの更新が適切に行われる可能性を高めることができる。 In addition, in this embodiment, if the producer is a predetermined producer (a producer with a predetermined level of reliability or higher) (S32: Yes), the update unit 46 updates the parameters based on the cultivation history entered by the producer (S34). This increases the likelihood that the parameters will be updated appropriately.

また、本実施形態では、出力部45は、生産者が選択した品種、栽培地点、栽培情報の組み合わせごとにシミュレーションを行い、シミュレーション結果を比較可能に出力する。これにより、生産者は、複数の品種、栽培地点、栽培情報の組み合わせごとのシミュレーション結果を比較して、どの品種をどのように栽培すべきかを判断することができる。 In addition, in this embodiment, the output unit 45 performs a simulation for each combination of variety, cultivation location, and cultivation information selected by the producer, and outputs the simulation results in a comparable manner. This allows the producer to compare the simulation results for multiple combinations of variety, cultivation location, and cultivation information, and determine which variety should be cultivated and how.

なお、上記実施形態では、サーバ10が更新部46を備える場合について説明したが、これに限られるものではない。すなわち、サーバ10は、更新部46を備えていなくてもよい。 In the above embodiment, the server 10 is described as having an update unit 46, but this is not limited to this. In other words, the server 10 does not have to have an update unit 46.

なお、上記実施形態では、本発明の農業支援装置の機能をサーバ10が有する場合について説明したが、これに限らず、利用者端末70が農業支援装置の機能を有していてもよい。すなわち、スタンドアローンの利用者端末70が単独で動作することで、上記処理を実現することとしてもよい。 In the above embodiment, the server 10 has the functions of the agricultural support device of the present invention, but this is not limited to this. The user terminal 70 may also have the functions of the agricultural support device. In other words, the above processing may be realized by the standalone user terminal 70 operating independently.

なお、上記の処理機能は、コンピュータによって実現することができる。その場合、処理装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記憶媒体(ただし、搬送波は除く)に記録しておくことができる。 The above processing functions can be realized by a computer. In this case, a program is provided that describes the processing details of the functions that the processing device should have. By executing the program on a computer, the above processing functions are realized on the computer. The program that describes the processing details can be recorded on a computer-readable storage medium (excluding carrier waves, however).

プログラムを流通させる場合には、例えば、そのプログラムが記録されたDVD(Digital Versatile Disc)、CD-ROM(Compact Disc Read Only Memory)などの可搬型記憶媒体の形態で販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。 When distributing a program, it is sold, for example, in the form of a portable storage medium on which the program is recorded, such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory). The program can also be stored in the storage device of a server computer, and then transferred from the server computer to other computers via a network.

プログラムを実行するコンピュータは、例えば、可搬型記憶媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記憶媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。 A computer that executes a program stores, for example, a program recorded on a portable storage medium or a program transferred from a server computer in its own storage device. The computer then reads the program from its own storage device and executes processing in accordance with the program. Note that the computer can also read the program directly from a portable storage medium and execute processing in accordance with that program. The computer can also execute processing in accordance with the received program each time a program is transferred from the server computer.

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。 The above-described embodiment is a preferred example of the present invention. However, it is not limited to this, and various modifications are possible within the scope of the gist of the present invention.

10 サーバ(農業支援装置)
42 環境データ取得部(取得部)
44 シミュレーション部(読み出し部、予測部)
45 出力部
50 パラメータDB(記憶部)
100 農業システム
10 Server (agricultural support device)
42 Environmental data acquisition unit (acquisition unit)
44 Simulation unit (reading unit, prediction unit)
45 Output section 50 Parameter DB (storage section)
100 Agricultural Systems

Claims (4)

ユーザが選択した農作物の品種と栽培地点との複数の組み合わせの情報を受け付け、
前記栽培地点それぞれにおける農作物の栽培環境に関する情報を取得し、
ユーザが入力した栽培方法に関する情報を取得し、
前記品種それぞれの特徴を示すパラメータを記憶部から読み出し、前記記憶部から読み出した前記パラメータと、取得した前記栽培環境に関する情報と、取得した前記栽培方法に関する情報と、に基づいて、農作物の各器官の日々の生育を予測する生育モデル作成し、前生育モデルから、農作物の各器官の日々の生育を前記品種、栽培地点、栽培方法に関する情報の組み合わせである栽培条件ごとに予測し、
前記予測により得られた前記農作物の各器官の日々の生育を示す情報を比較可能に出力し、
ユーザによって入力された農作物の栽培結果を示す情報のうち、所定以上の信頼度を有する予め定められたユーザによって入力された情報に基づいて前記パラメータを更新し、前記記憶部に記憶する、
処理をコンピュータに実行させるための農業支援プログラム。
Accepting information on multiple combinations of crop varieties and cultivation locations selected by the user;
Acquire information about the cultivation environment of the crops at each of the cultivation locations ;
Obtain information about the cultivation method entered by the user,
parameters indicating the characteristics of each of the varieties are read from a storage unit, a growth model is created that predicts the daily growth of each organ of the agricultural crop based on the parameters read from the storage unit, the acquired information on the cultivation environment, and the acquired information on the cultivation method, and the daily growth of each organ of the agricultural crop is predicted from the growth model for each cultivation condition that is a combination of information on the variety, cultivation location, and cultivation method;
outputting information indicating the daily growth of each organ of the crop obtained by the prediction in a comparable manner;
updating the parameters based on information input by a predetermined user having a predetermined reliability or higher among information indicating the cultivation results of the agricultural crops input by the user, and storing the updated parameters in the storage unit;
An agricultural support program that allows computers to perform processing.
前記予測の結果が前記予め定められたユーザによって入力された情報に近づくように、前記パラメータを更新する、ことを特徴とする請求項1に記載の農業支援プログラム。 The agricultural support program described in claim 1, characterized in that the parameters are updated so that the prediction results approach the information input by the predetermined user. ユーザが選択した農作物の品種と栽培地点との複数の組み合わせの情報を受け付け、
前記栽培地点それぞれにおける農作物の栽培環境に関する情報を取得し、
ユーザが入力した栽培方法に関する情報を取得し、
前記品種それぞれの特徴を示すパラメータを記憶部から読み出し、前記記憶部から読み出した前記パラメータと、取得した前記栽培環境に関する情報と、取得した前記栽培方法に関する情報と、に基づいて、農作物の各器官の日々の生育を予測する生育モデル作成し、前生育モデルから、農作物の各器官の日々の生育を前記品種、栽培地点、栽培方法に関する情報の組み合わせである栽培条件ごとに予測し、
前記予測により得られた前記農作物の各器官の日々の生育を示す情報を比較可能に出力し、
ユーザによって入力された農作物の栽培結果を示す情報のうち、所定以上の信頼度を有する予め定められたユーザによって入力された情報に基づいて前記パラメータを更新し、前記記憶部に記憶する、
処理をコンピュータが実行することを特徴とする農業支援方法。
Accepting information on multiple combinations of crop varieties and cultivation locations selected by the user;
Acquire information about the cultivation environment of the crops at each of the cultivation locations ;
Obtain information about the cultivation method entered by the user,
parameters indicating the characteristics of each of the varieties are read from a storage unit, a growth model is created that predicts the daily growth of each organ of the agricultural crop based on the parameters read from the storage unit, the acquired information on the cultivation environment, and the acquired information on the cultivation method, and the daily growth of each organ of the agricultural crop is predicted from the growth model for each cultivation condition that is a combination of information on the variety, cultivation location, and cultivation method;
outputting information indicating the daily growth of each organ of the crop obtained by the prediction in a comparable manner;
updating the parameters based on information input by a predetermined user having a predetermined reliability or higher among information indicating the cultivation results of the agricultural crops input by the user, and storing the updated parameters in the storage unit;
An agricultural support method characterized in that processing is executed by a computer.
ユーザが選択した農作物の複数の品種と栽培地点との複数の組み合わせの情報を受け付ける選択受付部と
前記栽培地点それぞれにおける農作物の栽培環境に関する情報を取得する環境情報取得部と、
ユーザが入力した栽培方法に関する情報を取得する栽培情報取得部と、
前記品種それぞれの特徴を示すパラメータを記憶部から読み出し、前記記憶部から読み出した前記パラメータと、取得した前記栽培環境に関する情報と、取得した前記栽培方法に関する情報と、に基づいて、農作物の各器官の日々の生育を予測する生育モデル作成し、前生育モデルから、農作物の各器官の日々の生育を前記品種、栽培地点、栽培方法に関する情報の組み合わせである栽培条件ごとに予測する予測部と、
前記予測により得られた農作物の各器官の日々の生育を示す情報を比較可能に出力する出力部と、
ユーザによって入力された農作物の栽培結果を示す情報のうち、所定以上の信頼度を有する予め定められたユーザによって入力された情報に基づいて前記パラメータを更新し、前記記憶部に記憶する更新部と、
を備える農業支援装置。
a selection receiving unit that receives information on a plurality of combinations of a plurality of varieties of agricultural crops and cultivation locations selected by a user;
an environmental information acquisition unit that acquires information about the cultivation environment of the agricultural crops at each of the cultivation locations ;
a cultivation information acquisition unit that acquires information about the cultivation method input by the user;
a prediction unit that reads out parameters indicating characteristics of each of the varieties from a storage unit, creates a growth model that predicts daily growth of each organ of the agricultural crop based on the parameters read out from the storage unit, the acquired information on the cultivation environment, and the acquired information on the cultivation method, and predicts daily growth of each organ of the agricultural crop from the growth model for each cultivation condition that is a combination of information on the variety, cultivation location, and cultivation method ;
an output unit that outputs information indicating the daily growth of each organ of the crop obtained by the prediction in a comparable manner;
an updating unit that updates the parameters based on information input by a predetermined user having a predetermined reliability or higher among information indicating the cultivation results of the agricultural crops input by the user, and stores the updated parameters in the storage unit;
An agricultural support device comprising:
JP2021185380A 2020-11-13 2021-11-15 Agricultural support program, agricultural support method, and agricultural support device Active JP7747324B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2024194384A JP2025010441A (en) 2020-11-13 2024-11-06 Agricultural support methods

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2020189462 2020-11-13
JP2020189462 2020-11-13
JP2021184702 2021-11-12
JP2021184702 2021-11-12

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2024194384A Division JP2025010441A (en) 2020-11-13 2024-11-06 Agricultural support methods

Publications (2)

Publication Number Publication Date
JP2022078987A JP2022078987A (en) 2022-05-25
JP7747324B2 true JP7747324B2 (en) 2025-10-01

Family

ID=81707211

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2021185380A Active JP7747324B2 (en) 2020-11-13 2021-11-15 Agricultural support program, agricultural support method, and agricultural support device
JP2024194384A Pending JP2025010441A (en) 2020-11-13 2024-11-06 Agricultural support methods

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2024194384A Pending JP2025010441A (en) 2020-11-13 2024-11-06 Agricultural support methods

Country Status (1)

Country Link
JP (2) JP7747324B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7748102B2 (en) * 2022-10-20 2025-10-02 国立研究開発法人農業・食品産業技術総合研究機構 Prediction method and prediction program, and environmental adjustment method and environmental adjustment program
CN116107248B (en) * 2023-02-01 2025-08-19 青岛海尔空调器有限总公司 Crop storage monitoring method, apparatus and storage medium
JP2026038460A (en) * 2024-08-23 2026-03-06 株式会社クボタ Optical simulation method, optical simulation system, and computer program
CN119313497A (en) * 2024-12-13 2025-01-14 凉山彝族自治州农业科学研究院 System and method for guiding crop-economy integrated planting based on meteorological data

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018013348A (en) 2016-07-19 2018-01-25 日本電信電話株式会社 Striking inspection system and percussion inspection method
JP2019219917A (en) 2018-06-20 2019-12-26 株式会社東芝 Display system, wearable device and monitoring control device
JP2019219704A (en) 2018-06-15 2019-12-26 株式会社オーガニックnico Farm management support system
JP2020024703A (en) 2018-08-03 2020-02-13 三菱ケミカル株式会社 Production management system, management method, and program

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5589493B2 (en) * 2010-03-26 2014-09-17 株式会社大林組 Turfgrass community growth prediction method
JP6634728B2 (en) * 2015-08-10 2020-01-22 富士電機株式会社 Environmental control system
JP2019030253A (en) * 2017-08-08 2019-02-28 株式会社誠和 Cultivation support apparatus, cultivation support method, computer program, and recording medium
JP7557853B2 (en) * 2019-08-02 2024-09-30 国立研究開発法人農業・食品産業技術総合研究機構 Yield estimation program, yield estimation method, and yield estimation device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018013348A (en) 2016-07-19 2018-01-25 日本電信電話株式会社 Striking inspection system and percussion inspection method
JP2019219704A (en) 2018-06-15 2019-12-26 株式会社オーガニックnico Farm management support system
JP2019219917A (en) 2018-06-20 2019-12-26 株式会社東芝 Display system, wearable device and monitoring control device
JP2020024703A (en) 2018-08-03 2020-02-13 三菱ケミカル株式会社 Production management system, management method, and program

Also Published As

Publication number Publication date
JP2025010441A (en) 2025-01-20
JP2022078987A (en) 2022-05-25

Similar Documents

Publication Publication Date Title
JP7747324B2 (en) Agricultural support program, agricultural support method, and agricultural support device
Craparo et al. Coffea arabica yields decline in Tanzania due to climate change: Global implications
Sattar et al. Responses of phenology, yield attributes, and yield of wheat varieties under different sowing times in Indo-Gangetic Plains
JP7544392B2 (en) Method, device, and program for generating a prediction model for predicting agricultural crop production performance
Heinemann et al. Characterization of drought stress environments for upland rice and maize in central Brazil
JP7557853B2 (en) Yield estimation program, yield estimation method, and yield estimation device
JP6825571B2 (en) Information system, farming support method, information processing device, its control method, and control program
Kang et al. The derivation of sink functions of wheat organs using the GREENLAB model
Singels et al. Farm level decision support for sugarcane irrigation management during drought
JP7604001B2 (en) Method, device, and program for predicting agricultural crop production performance
JP7505144B2 (en) Information processing method, information processing program, and information processing device
Lindemann-Zutz et al. Prediction of time to harvest and its variability of broccoli (Brassica oleracea var. italica) part II. Growth model description, parameterisation and field evaluation
Oliveira et al. Yield and weight distribution of two potato cultivars grown from seed potatoes of different physiological ages
JP7386503B2 (en) Yield estimation program, yield estimation method, and yield estimation device
Cornet et al. Expected yield and economic improvements of a yam seed system in West Africa using agro‐physiological modelling
JP7523841B2 (en) Agricultural support program, agricultural support method, and agricultural support device
JP7773242B2 (en) Parameter calculation method, agricultural support information output method, parameter calculation program, and agricultural support information output program
JP2022153260A (en) Agriculture support program and agriculture support method
JP7370519B2 (en) Prediction program, prediction method and prediction device
Ji et al. Prediction of grain number and grain weight across different spikelet positions in winter wheat under short-term low-temperature stress
CA3245761A1 (en) Non-transitory computer-readable recodring medium, agricultural support method, and agricultural support device
Kurniasih et al. Technical efficiency in several levels of adoption of garlic farming standard operating procedures in production centre in Indonesia
JP2021087424A (en) Yield regulation method, yield regulation supporting method, and yield regulation supporting program, and information processing device
JP2022093884A (en) Growth state prediction method and growth state prediction program
Cadena Modeling the Dynamics of Starch Content in Cassava

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240517

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20240517

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20240517

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20240517

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240611

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240711

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20240827

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241106

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20241114

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20241122

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250623

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250910

R150 Certificate of patent or registration of utility model

Ref document number: 7747324

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150