特許 7694953 タマネギの収穫情報予測方法及びタマネギの収穫情報予測プログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2025-06-10

(45)【発行日】2025-06-18

(54)【発明の名称】タマネギの収穫情報予測方法及びタマネギの収穫情報予測プログラム

(51)【国際特許分類】

   G06Q 50/02 20240101AFI20250611BHJP

   G06Q 10/04 20230101ALI20250611BHJP

【ＦＩ】

G06Q50/02

G06Q10/04

【請求項の数】 11

(21)【出願番号】P 2021191111

(22)【出願日】2021-11-25

(65)【公開番号】P2022122822

(43)【公開日】2022-08-23

【審査請求日】2024-05-08

(31)【優先権主張番号】P 2021020072

(32)【優先日】2021-02-10

(33)【優先権主張国・地域又は機関】JP

(73)【特許権者】

【識別番号】501203344

【氏名又は名称】国立研究開発法人農業・食品産業技術総合研究機構

(74)【代理人】

【識別番号】110004370

【氏名又は名称】弁理士法人片山特許事務所

(72)【発明者】

【氏名】山本 岳彦

(72)【発明者】

【氏名】室 崇人

(72)【発明者】

【氏名】木下 貴文

(72)【発明者】

【氏名】山内 大輔

【審査官】山崎 誠也

(56)【参考文献】

【文献】特開２０２０－０２４７０３（ＪＰ，Ａ）

【文献】特開２０１９－１０３５２８（ＪＰ，Ａ）

【文献】特開２０２０－０５４２８９（ＪＰ，Ａ）

【文献】中ネギの９～10月まき作型における第１回目収穫日の推定法，[online]，2017年，[令和７年２月２１日検索], インターネット<URL:https://www.pref.kagoshima.jp/ag11/pop-tech/nenndo/documents/documents/105866_20230418141122-1.pdf>

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｑ １０／００－９９／００

(57)【特許請求の範囲】

【請求項1】

タマネギの葉鞘基部径又はりん茎の球径の実測値と実測時期の情報の入力を受け付け、
前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、
前記実測値と前記実測値に対応する積算気温とを用いて、りん茎が肥大を開始した後の積算気温と前記りん茎の球径との関係を示す第１式を特定し、
前記第１式を用いて、予め定められた収穫時の積算気温に到達するときのりん茎の球径を特定する、
処理をコンピュータが実行することを特徴とするタマネギの収穫情報予測方法。

【請求項2】

前記予め定められた収穫時の積算気温に到達するときのりん茎の球径から、収穫時のりん茎の重量を特定する、処理を前記コンピュータが更に実行することを特徴とする請求項１に記載のタマネギの収穫情報予測方法。

【請求項3】

前記第１式における、前記積算気温の変化に対する前記りん茎の球径の変化を示す値が、予め定められた値である、ことを特徴とする請求項１又は２に記載のタマネギの収穫情報予測方法。

【請求項4】

前記実測時期が前記りん茎が肥大開始するまでの時期であった場合、
前記実測値と前記実測値に対応する積算気温とに基づいて、りん茎が肥大開始するまでの積算気温と葉鞘基部径との関係を示す第２式を特定するとともに、前記第２式を用いて、前記葉鞘基部径が前記りん茎の肥大開始を示す第１の値になるときの積算気温を特定し、
前記第１の値と、前記葉鞘基部径が前記第１の値になるときの積算気温とを用いて、前記第１式を特定する、ことを特徴とする請求項１～３のいずれか一項に記載のタマネギの収穫情報予測方法。

【請求項5】

前記第２式における、前記積算気温の変化に対する前記葉鞘基部径の変化を示す値が、予め定められた値である、ことを特徴とする請求項４に記載のタマネギの収穫情報予測方法。

【請求項6】

前記タマネギの栽培地点における気温の実測データ及び予測データと、前記予め定められた収穫時の積算気温と、に基づいて収穫に適した日を予測する、処理を前記コンピュータが更に実行することを特徴とする請求項１～５のいずれか一項に記載のタマネギの収穫情報予測方法。

【請求項7】

タマネギの葉鞘基部径又はりん茎の球径の実測値と実測時期の情報の入力を受け付け、
前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、
前記実測値と前記実測値に対応する積算気温とを用いて、りん茎が肥大を開始した後の積算気温と前記りん茎の球径との関係を示す第１式を特定し、
前記第１式を用いて、予め定められた収穫時の積算気温に到達するときのりん茎の球径を特定する、
処理をコンピュータに実行させることを特徴とするタマネギの収穫情報予測プログラム。

【請求項8】

タマネギの葉鞘基部径又はりん茎の球径の実測値と実測時期の情報の入力を受け付け、
前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、
前記実測値と前記実測値に対応する積算気温とを用いて、定植時以降の積算気温と前記りん茎の球径との関係を示す指数関数式を特定し、
前記指数関数式を用いて、予め定められた収穫時の積算気温に到達するときのりん茎の球径を特定する、
処理をコンピュータが実行することを特徴とするタマネギの収穫情報予測方法。

【請求項9】

タマネギに関する生育指標の実測値と実測時期の情報の入力を受け付け、
前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、
前記実測値と前記実測値に対応する積算気温とを用いて、りん茎が肥大を開始した後の積算気温と前記生育指標との関係を示す第１式を特定し、
前記第１式を用いて、予め定められた収穫時の積算気温に到達するときの前記生育指標を特定する、
処理をコンピュータが実行し、
前記実測時期が前記りん茎が肥大開始するまでの時期であった場合、
前記実測値と前記実測値に対応する積算気温とに基づいて、りん茎が肥大開始するまでの積算気温と生育指標との関係を示す第２式を特定するとともに、前記第２式を用いて、前記生育指標が前記りん茎の肥大開始を示す第１の値になるときの積算気温を特定し、
前記第１の値と、前記生育指標が前記第１の値になるときの積算気温とを用いて、前記第１式を特定する、ことを特徴とするタマネギの収穫情報予測方法。

【請求項10】

タマネギに関する乾物重、新鮮重、草丈、葉面積のいずれかである生育指標の実測値と実測時期の情報の入力を受け付け、
前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、
前記実測値と前記実測値に対応する積算気温とを用いて、定植時以降の積算気温と前記生育指標との関係を示す指数関数式を特定し、
前記指数関数式を用いて、予め定められた収穫時の積算気温に到達するときの前記生育指標を特定する、
処理をコンピュータが実行することを特徴とするタマネギの収穫情報予測方法。

【請求項11】

前記生育指標は、乾物重、新鮮重、草丈、葉面積のいずれかであることを特徴とする請求項９に記載のタマネギの収穫情報予測方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、タマネギの収穫情報予測方法及びタマネギの収穫情報予測プログラムに関する。

【背景技術】

【0002】

近年、東北地域においてタマネギ栽培が広まりつつあるが、東北地域は広いため多様な気象条件があり、また、新興産地であるために生産者が経験不足であるなどの理由から、安定生産が難しい。このため、栽培前に予め産地ごとの収穫時期や収穫量を把握して作付計画を策定できることが好ましい。

【0003】

従来、タマネギの収穫日や収穫量は、タマネギ産地近傍の公設試験研究機関やＪＡ等で実施した栽培試験の結果（定植時期、品種、栽植密度、平均収量等）に基づく栽培暦に、産地の気象条件を勘案することで予測していた。

【0004】

一方、最近では、一部の野菜品目ではあるが、過去の気象・生育データ等から生育モデルを用いて収量を予測する取り組みが進められている（例えば、特許文献１等参照）。

【先行技術文献】

【特許文献】

【0005】

【文献】特開２０２０－０５３０３９号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

上述した過去の栽培試験の結果に基づく予測方法では、近年の異常気象によるずれや新興産地、新作型への対応が難しく、それぞれについて栽培試験を行い、栽培暦を作成する必要がある。また、栽培試験を行った地点を含む地域における収穫日や収穫量の推定精度は高いものの、土壌や気象条件が異なると推定精度が低くなる。

【0007】

また、生育モデルを用いて収量を予測する方法は、膨大なデータを入力する必要があり、煩雑である。

【0008】

そこで、本発明は、簡易に、タマネギの収穫情報を予測することが可能なタマネギの収穫情報予測方法及びタマネギの収穫情報予測プログラムを提供することを目的とする。

【課題を解決するための手段】

【0009】

本発明のタマネギの収穫情報予測方法は、タマネギの葉鞘基部径又はりん茎の球径の実測値と実測時期の情報の入力を受け付け、前記実測時期に対応する定植時からの積算気温を、前記実測値に対応する積算気温として特定し、前記実測値と前記実測値に対応する積算気温とを用いて、りん茎が肥大を開始した後の積算気温と前記りん茎の球径との関係を示す第１式を特定し、前記第１式を用いて、予め定められた収穫時の積算気温に到達するときのりん茎の球径を特定する、処理をコンピュータが実行するタマネギの収穫情報予測方法である。

【発明の効果】

【0010】

本発明のタマネギの収穫情報予測方法及びタマネギの収穫情報予測プログラムは、簡易に、タマネギの収穫情報を予測することができるという効果を奏する。

【図面の簡単な説明】

【0011】

【図1】第１の実施形態に係る生育状態予測システムの構成が概略的に示す図である。

【図2】図２（ａ）は、サーバのハードウェア構成を示す図であり、図２（ｂ）は、利用者端末のハードウェア構成を示す図である。

【図3】第１の実施形態に係るサーバの機能ブロック図である。

【図4】ある品種のタマネギを所定条件で栽培したときの、定植時からの積算気温と葉鞘基部径・球径との関係を示すグラフである。

【図5】第１の実施形態におけるサーバの処理を示すフローチャートである。

【図6】処理に用いる定数を示す図である。

【図7】図７（ａ）、図７（ｂ）は、図５のステップＳ２０，Ｓ２２の処理を説明するための図である。

【図8】図５のステップＳ２４の処理を説明するための図である。

【図9】図９（ａ）、図９（ｂ）は、図５のステップＳ１６の判断が否定されたときの処理を説明するための図である。

【図10】図１０（ａ）～図１０（ｃ）は、第１の実施形態の変形例について説明するための図である。

【図11】第２の実施形態に係るサーバの機能ブロック図である。

【図12】図１２（ａ）～図１２（ｃ）は、第２の実施形態について説明するための図である。

【図13】第２の実施形態におけるサーバの処理を示すフローチャートである。

【図14】第２の実施形態の変形例を説明するための図である。

【発明を実施するための形態】

【0012】

《第１の実施形態》
以下、第１の実施形態に係るタマネギの収穫情報予測方法について、詳細に説明する。

【0013】

図１には、第１の実施形態に係る収穫情報予測システムの構成が概略的に示されている。本実施形態の収穫情報予測システム１００は、タマネギを栽培する生産者がタマネギの作付け計画や出荷計画を策定する際に、ある栽培地点において、いつ定植すると、いつ収穫日（収穫適期）が到来し、そのときのタマネギのりん茎の球径がどのくらいの大きさで、収量がどの程度になるのかを確認するためのシステムである。

【0014】

収穫情報予測システム１００は、図１に示すように、サーバ１０と、利用者端末７０と、を備える。サーバ１０及び利用者端末７０は、インターネットなどのネットワーク８０に接続されている。

【0015】

サーバ１０は、データセンタ等に設置される情報処理装置であり、タマネギを栽培する地点（栽培地点）の環境情報（平年値や、過去データ、予測データなど）や、利用者端末７０から入力される情報を用い、栽培地点におけるタマネギの収穫日や、収穫日におけるタマネギのりん茎の大きさ（球径）、収量を推定し、推定結果を利用者端末７０に出力する。

【0016】

図２（ａ）には、サーバ１０のハードウェア構成が示されている。図２（ａ）に示すように、サーバ１０は、コンピュータとしてのＣＰＵ（Central Processing Unit）９０、ＲＯＭ（Read Only Memory）９２、ＲＡＭ（Random Access Memory）９４、記憶部（ＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等）９６、ネットワークインタフェース９７、及び可搬型記憶媒体用ドライブ９９等を備えている。これらサーバ１０の構成各部は、バス９８に接続されている。サーバ１０では、ＲＯＭ９２あるいは記憶部９６に格納されているプログラム（収穫情報予測プログラムを含む）、或いは可搬型記憶媒体用ドライブ９９が可搬型記憶媒体９１から読み取ったプログラム（収穫情報予測プログラムを含む）をＣＰＵ９０が実行することにより、図３に示す各部の機能が実現される。なお、図３の各部の機能は、例えば、ＡＳＩＣ（Application Specific Integrated Circuit）やＦＰＧＡ（Field Programmable Gate Array）等の集積回路により実現されてもよい。なお、図３の各部の詳細については、後述する。

【0017】

利用者端末７０は、生産者等が利用可能なＰＣ（Personal Computer）等の端末であり、生産者等による作物の生体情報の入力を受け付けてサーバ１０に対して送信したり、サーバ１０の推定結果を受信して、表示したりする。利用者端末７０は、図２（ｂ）に示すようなハードウェア構成を有する。具体的には、利用者端末７０は、図２（ｂ）に示すように、ＣＰＵ１９０、ＲＯＭ１９２、ＲＡＭ１９４、記憶部（ＨＤＤ等）１９６、ネットワークインタフェース１９７、表示部１９３、入力部１９５、及び可搬型記憶媒体１９１の読み取りが可能な可搬型記憶媒体用ドライブ１９９等を備えている。これら利用者端末７０の構成各部は、バス１９８に接続されている。表示部１９３は、液晶ディスプレイ等を含み、入力部１９５は、キーボード、マウス、タッチパネル等を含む。

【0018】

図３には、サーバ１０の機能ブロック図が示されている。図３に示すように、サーバ１０では、ＣＰＵ９０がプログラムを実行することにより、入力受付部２０、実測値取得時期判定部２２、肥大開始前関係式特定部２４、肥大開始後関係式特定部２６、収穫時球径特定部２８、収穫時りん茎重特定部３０、収量推定部３２、収穫日特定部３４、出力部３６としての機能が実現されている。なお、図３には、サーバ１０の記憶部９６等に格納されている定数テーブル４２や気温ＤＢ４０も図示されている。

【0019】

入力受付部２０は、利用者端末７０を介して生産者等が入力する情報を取得する。生産者等が入力する情報には、圃場の位置情報、栽培しているタマネギの品種情報、定植日の情報、栽植密度の情報などが含まれる。また、生産者等が入力する情報には、生産者等が葉鞘基部径又はりん茎の球径を計測した結果（実測値）と、測定日の情報が含まれる。

【0020】

ここで、定数テーブル４２には、圃場の位置やタマネギの品種、定植時期（春まき作か、秋まき作か）に紐づけて各種定数（例えば、図６に示す定数）が格納されている。入力受付部２０は、生産者等が入力した情報に対応する定数を定数テーブル４２から取得し、実測値取得時期判定部２２、肥大開始前関係式特定部２４、肥大開始後関係式特定部２６、収穫時りん茎重特定部３０、収穫日特定部３４に受け渡す。

【0021】

実測値取得時期判定部２２は、生産者等が入力した実測値がりん茎肥大の開始前の実測値であるのか、りん茎肥大の開始後の実測値であるのかを判定する。

【0022】

図４は、ある品種のタマネギを所定条件で栽培したときの、定植時からの積算気温（日平均気温の積算値）と葉鞘基部径及びりん茎の球径の大きさ（葉鞘基部径・球径）との関係を示すグラフである。図４において２直線が交わる点がりん茎の肥大開始時を示している。この肥大開始時の葉鞘基部径（Ｄｂ）は、品種や栽培時期、土壌などにより定まる定数であると考えられる。

【0023】

したがって、実測値取得時期判定部２２は、実測値がＤｂ以下であれば、測定日がりん茎肥大の開始前であると判定し、葉鞘基部径が所定値よりも大きければ、測定日がりん茎肥大の開始後であると判定する。

【0024】

肥大開始前関係式特定部２４は、実測値の測定日がりん茎肥大の開始前である場合に、実測値と測定日における定植日からの積算気温とに基づいて、りん茎肥大が起こる前の時期における積算気温と葉鞘基部径の関係を示す第２式（肥大開始前関係式と呼ぶ）を特定する。

【0025】

本発明者らは、定植時からの積算気温の変化に対するタマネギのりん茎径の変化を示す傾きは、肥大開始前と後では大きく異なるものの、それぞれの傾きは、品種や定植場所が同じであれば、定植日が多少ずれてもほとんど変化しないことを見出した。すなわち、図４において、肥大開始前の期間の直線の傾き（積算気温の変化に対する葉鞘基部径の大きさの変化の割合）ａ１は、品種や栽培時期、土壌などにより定まるが、切片は定植日によって異なることを意味する。

【0026】

したがって、肥大開始前関係式特定部２４は、式ｙ＝ａ１・ｘ＋ｂ１の、ｘに実測値、ｙに測定日における定植時からの積算気温を代入することでｂ１を求め、肥大開始前関係式（ｙ＝ａ１・ｘ＋ｂ１）を特定する。

【0027】

肥大開始後関係式特定部２６は、りん茎肥大の開始後の時期における、積算気温とりん茎の球径との関係を示す第１式（肥大開始後関係式と呼ぶ）を特定する。

【0028】

収穫時球径特定部２８は、定植時からの積算気温が所定の値（後述するＴｈｖ）に到達するときを収穫時（収穫適期）とみなし、肥大開始後関係式特定部２６が特定した肥大開始後関係式を用いて、収穫時におけるりん茎の球径を特定する。なお、値Ｔｈｖは、品種や栽培時期、土壌などに固有の値である。

【0029】

収穫時りん茎重特定部３０は、収穫時球径特定部２８が特定した収穫時におけるりん茎の球径を用いて、収穫時におけるりん茎重を次式（１）から特定する。
りん茎重＝４／３×π×（りん茎の球径／２）³×ｋ …（１）

【0030】

なお、上式（１）の係数ｋは、事前の栽培試験等において定まる値であるものとする。

【0031】

収量推定部３２は、収穫時におけるりん茎重と、生産者等が入力した栽植密度とに基づいて、収量を次式（２）から推定する。
収量＝りん茎重×栽植密度×商品球率 …（２）

【0032】

なお、商品球率（収穫したタマネギのうち商品として販売可能なタマネギの割合）は、事前の調査により定められた値であるとする。

【0033】

収穫日特定部３４は、気温ＤＢ４０を参照して、定植時からの積算気温が所定値（Ｔｈｖ）になる日を収穫日として特定する。収穫日特定部３４は、定植日からの積算気温が所定値（Ｔｈｖ）になる日を、気温ＤＢ４０に格納されている気温の過去データや将来の予測データに基づいて特定する。

【0034】

出力部３６は、収量推定部３２が推定した収量や、収穫日特定部３４が特定した収穫日の情報を、利用者端末７０に対して出力する。これにより、利用者端末７０の表示部１９３上に、収量や収穫日の推定結果の情報が表示されるようになっている。

【0035】

（サーバ１０の処理について）
以下、サーバ１０の処理について、図５のフローチャートに沿って、その他図面を適宜参照しつつ詳細に説明する。

【0036】

図５の処理においては、まずステップＳ１０において、入力受付部２０が、利用者端末７０から入力があるまで待機する。利用者端末７０において、生産者等が情報を入力すると、ステップＳ１２に移行する。なお、生産者等が入力する情報には、前述のように、圃場の位置情報、栽培しているタマネギの品種の情報、定植日の情報、栽植密度の情報、葉鞘基部径又はりん茎の球径の実測値、測定日の情報が含まれる。

【0037】

ステップＳ１２に移行すると、入力受付部２０は、利用者端末７０から入力された情報を取得し、取得した情報に対応する定数を取得する。例えば、入力受付部２０は、定数テーブル４２を参照して、取得した圃場の位置情報や、品種の情報、定植日の情報に対応する定数として、図６に示すような定数ａ１，ａ２，Ｄｂ，Ｔｈｖを取得したものとする。なお、入力受付部２０は、取得した定数を実測値取得時期判定部２２、肥大開始前関係式特定部２４、肥大開始後関係式特定部２６、収穫時りん茎重特定部３０、収穫日特定部３４に受け渡す。

【0038】

次いで、ステップＳ１４では、実測値取得時期判定部２２が、実測値を取得した時期が肥大開始前か後かを判定する。ここでは、実測値取得時期判定部２２は、ステップＳ１２で取得したＤｂよりも取得した実測値が大きいか小さいかにより、実測値を取得した時期が肥大開始前か後かを判定する。

【0039】

次いで、ステップＳ１６では、実測値取得時期判定部２２が、肥大開始前か否かを判断する。例えば、実測値がＤｂ（＝１７ｍｍ）よりも小さい場合には、ステップＳ１６の判断は肯定され、ステップＳ１８に移行する。

【0040】

ステップＳ１８に移行すると、肥大開始前関係式特定部２４は、実測値から肥大開始前関係式を特定する。具体的には、肥大開始前関係式特定部２４は、気温ＤＢ４０の過去データに基づいて、測定日における定植日からの積算気温を算出する。そして、肥大開始前関係式特定部２４は、図７（ａ）に示すように、実測値と算出した積算気温とを、式ｙ＝ａ１・ｘ＋ｂ１に代入することで、切片ｂ１を求め、肥大開始前関係式を特定する。例えば、肥大開始前関係式として、次式（３）が特定されたとする。
ｙ＝０．０１４×ｘ＋４．９５ …（３）

【0041】

次いで、ステップＳ２０では、肥大開始後関係式特定部２６が、肥大開始前関係式（上式（３））から肥大開始時の積算気温を求め、肥大開始後関係式を特定する。具体的には、肥大開始後関係式特定部２６は、上式（３）のｙにＤｂ＝１７を代入することで、肥大開始時の積算気温ｘ≒８６０を求める。そして、肥大開始後関係式特定部２６は、図７（ｂ）に示すように、（ｘ，ｙ）＝（８６０，１７）を肥大開始後関係式ｙ＝ａ２・ｘ＋ｂ２に代入することで、切片ｂ２を求め、肥大開始後関係式を特定する。例えば、肥大開始後関係式として、次式（４）が特定されたとする。
ｙ＝０．０８５×ｘ－５３．６１ …（４）

【0042】

次いで、ステップＳ２４では、収穫時球径特定部２８が、肥大開始後関係式（上式（４））を用いて、収穫時球径を特定する。この場合、収穫時球径特定部２８は、定植時からの積算気温が図６のＴｈｖに到達するときを収穫時（収穫適期）とみなし、肥大開始後関係式（上式（４））にｘ＝Ｔｈｖを代入することで、収穫時におけるりん茎の球径を特定する。図８に示すように、Ｔｈｖ＝１７００℃の場合には、この値を上式（４）のｘに代入することで、収穫時におけるりん茎の球径（ｙ≒９１ｍｍ）を得る。

【0043】

次いで、ステップＳ２６では、収穫時りん茎重特定部３０が、収穫時球径を用いて、収穫時りん茎重を特定する。この場合、収穫時りん茎重特定部３０は、上式（１）に、収穫時におけるりん茎の球径（ｙ≒９１ｍｍ）を代入することで、収穫時りん茎重を特定する。

【0044】

次いで、ステップＳ２８では、収量推定部３２が、ステップＳ２６で特定された収穫時りん茎重を用いて、上式（２）より、収量を推定する。

【0045】

次いで、ステップＳ３０では、収穫日特定部３４が、収穫日を特定する。収穫日特定部３４は、気温ＤＢ４０の過去データや将来データを参照して、定植日からの積算気温がＴｈｖ＝１７００℃になる日が２０１９年７月２４であった場合には、その日を収穫適期（収穫日）と特定する。

【0046】

次いで、ステップＳ３２では、出力部３６が、利用者端末７０に対して、ステップＳ２８で推定された収量の情報やステップＳ３０で特定された収穫日の情報を表示する画面を出力する。これにより、生産者等は、いつ、どのくらいの収量が得られるかを確認することができる。以上により、図５の処理が終了する。

【0047】

一方、ステップＳ１６の判断が否定された場合、すなわち、実測値の測定日が肥大開始後であった場合には、ステップＳ２２に移行する。

【0048】

この場合、ステップＳ２２に移行すると、肥大開始後関係式特定部２６は、実測値から肥大開始後関係式を特定する。具体的には、肥大開始後関係式特定部２６は、気温ＤＢ４０の過去データに基づいて、測定日における定植日からの積算気温を算出する。そして、肥大開始後関係式特定部２６は、算出した積算気温と実測値とを、肥大開始後関係式ｙ＝ａ２・ｘ＋ｂ２に代入することで、図９（ａ）に示すように、肥大開始後関係式を特定する（この場合、肥大開始前関係式は特定しない）。ここでは、肥大開始後関係式として、次式（５）が特定できたとする。
ｙ＝０．０８５×ｘ－５７．６４ …（５）

【0049】

次いで、ステップＳ２４では、収穫時球径特定部２８が、肥大開始後関係式（上式（５））を用いて、収穫時球径を特定する。この場合、収穫時球径特定部２８は、定植時からの積算気温が図６のＴｈｖに到達するときを収穫時（収穫適期）とみなし、肥大開始後関係式（上式（５））にｘ＝Ｔｈｖを代入することで、収穫時におけるりん茎の球径を特定する。図９（ｂ）に示すように、Ｔｈｖ＝１７００℃の場合には、この値を上式（５）のｘに代入することで、収穫時におけるりん茎の球径（ｙ≒８７ｍｍ）を得る。

【0050】

その後の、ステップＳ２６、Ｓ２８，Ｓ３０，Ｓ３２の処理については、上述したのと同様の処理となっている。

【0051】

以上詳細に説明したように、本第１の実施形態によると、入力受付部２０が、タマネギの葉鞘基部径又はりん茎の球径の実測値と測定日（実測時期）の情報を取得し（Ｓ１２）、肥大開始前関係式特定部２４又は肥大開始後関係式特定部２６が、測定日の定植時からの積算気温（実測値に対応する積算気温）を特定する（Ｓ１８、Ｓ２２）。そして、肥大開始後関係式特定部２６は、実測値と実測値に対応する積算気温とを用いて、肥大開始後関係式を特定し（Ｓ１８及びＳ２０、又はＳ２２）、収穫時球径特定部２８は、肥大開始後関係式を用いて、予め定められた収穫時の積算気温（Ｔｈｖ）になるときのりん茎の球径を特定する（Ｓ２４）。このように、本第１の実施形態では、少なくとも１回の葉鞘基部径又はりん茎の球径の実測値を用いることで、簡易に、収穫時りん茎重を精度よく推定することができる。また、図６に示すような定数を、様々な品種のタマネギを各種条件下で栽培したときに求め、定数テーブル４２に格納しておくことで、様々な品種のタマネギを様々な条件下で栽培した場合でも、収穫時の球径や重量、収量や収穫日を簡易に算出することができる。

【0052】

また、本第１の実施形態では、収穫時りん茎重特定部３０は、定植時からの積算気温が予め定めた値（Ｔｈｖ）になるときのりん茎の球径から、上式（１）に基づいて、収穫時のりん茎の重量を特定する。これにより、収穫時のりん茎の重量を精度よく特定することができる。

【0053】

また、本第１の実施形態では、肥大開始後関係式の傾きａ２が、品種や栽培時期、土壌などに応じた定数とされている。このように、品種や栽培時期、土壌などの組み合わせごとに傾きを事前に用意しておくことで、肥大開始後に複数回実測しなくても肥大開始後関係式を特定することが可能となる。

【0054】

また、本第１の実施形態では、入力された実測値がりん茎の肥大開始前に得られた値であった場合には、肥大開始前関係式特定部２４は、実測値から肥大開始前関係式を特定する（Ｓ１８）。そして、肥大開始後関係式特定部２６が、肥大開始前関係式から肥大開始時の積算気温を求め、求めた肥大開始時の積算気温を用いて肥大開始後関係式を特定する（Ｓ２０）。これにより、りん茎が肥大する前の葉鞘基部径の実測値が入力された場合でも、肥大開始後関係式を特定することができる。

【0055】

また、本第１の実施形態では、肥大開始前関係式の傾きａ１が、品種や栽培時期、土壌などに応じた定数とされている。このように、品種や栽培時期、土壌などの組み合わせごとに傾きａ１を事前に用意しておくことで、肥大開始前に複数回実測しなくても肥大開始前関係式を特定することが可能となる。

【0056】

また、本第１の実施形態では、生産者等が入力する実測値は、りん茎の肥大開始前の実測値であっても肥大開始後の実測値であってもよい。このため、生産者等が葉鞘基部径又はりん茎の球径を実測するタイミングに自由度を持たせることができる。

【0057】

（変形例）
なお、上記第１の実施形態では、サーバ１０は、タマネギの葉鞘基部径又はりん茎の球径の実測値と、実測値の測定日における定植時からの積算気温の情報と、に基づいて、収穫時のりん茎の球径を特定（推定）する場合について説明した。しかしながら、これに限られるものではなく、サーバ１０は、タマネギに関する生育指標（例えば、乾物重、新鮮重、草丈、葉面積など）の実測値と、実測値の測定日における定植時からの積算気温の情報と、に基づいて、収穫時の生育指標（例えば、乾物重、新鮮重、草丈、葉面積など）を特定することも可能である。

【0058】

以下、タマネギに関する生育指標が、乾物重である場合について説明する。

【0059】

本発明者らは、定植時からの積算気温の変化に対するタマネギの乾物重の変化を示す傾きは、りん茎の肥大開始前と後では大きく異なるものの、それぞれの傾きは、品種や定植場所が同じであれば、定植日が多少ずれてもほとんど変化しないことを見出した。すなわち、図１０（ａ）において、りん茎の肥大開始前の期間の直線の傾き（積算気温の変化に対する乾物重の変化の割合）α１は、品種や栽培時期、土壌などにより定まるが、切片β１は定植日によって異なることを意味する。したがって、式ｙ＝α１・ｘ＋β１の、ｘに乾物重の実測値、ｙに実測値の測定日における定植時からの積算気温を代入することでβ１を求めれば、肥大開始前関係式（ｙ＝α１・ｘ＋β１）を特定することができる。また、図１０（ａ）において、肥大開始後の期間の直線の傾き（積算気温の変化に対する乾物重の変化の割合）α２は、品種や栽培時期、土壌などにより定まるが、切片は定植日によって異なる。したがって、式ｙ＝α２・ｘ＋β２のｘ，ｙに値を代入することで、切片β２を求めれば、肥大開始後関係式（ｙ＝α２・ｘ＋β２）を特定することができる。なお、肥大開始後関係式のｘ，ｙに代入する値は、生産者等が肥大開始後の実測値を入力した場合には、測定日の定植日からの積算気温と乾物重の実測値であり、生産者等が肥大開始前の実測値を入力した場合には、肥大開始前関係式から得られる肥大開始時の積算気温とδｂ（肥大開始時の乾物重）の値である。

【0060】

例えば、取得した圃場の位置情報や、品種の情報、定植日の情報に対応する定数として、図１０（ｂ）示すような定数α１，α２，δｂ，Ｔｈｖ（Ｔｈｖは、予め定められた収穫時の積算気温）が定数テーブル４２に格納されていたとする。そして、生産者等が入力した情報に基づいて、図１０（ｃ）に示すような乾物重と定植時からの積算気温との関係式が得られたとする。この場合、肥大開始後関係式にｘ＝Ｔｈｖ＝１７００を代入することで、収穫時における乾物重ｙをｙ＝０．０３４２×１７００－２４．７４６≒３３．４（ｇ）と特定（推定）することができる。

【0061】

以上のように、本変形例では、少なくとも１回の乾物重の実測値を用いることで、簡易に、収穫時の乾物重を精度よく推定することができる。また、図１０（ｂ）に示すような定数を、様々な品種のタマネギを各種条件下で栽培したときに求めておくことで、様々な品種のタマネギを様々な条件下で栽培した場合でも、収穫時の乾物重を簡易に推定することができる。

【0062】

なお、タマネギに関するその他の生育指標（新鮮重、草丈、葉面積など）についても、乾物重の場合と同様である。各生育指標の積算気温に対する変化は、図１０（ａ）と同様に変化するため、定数α１，α２，δｂ，Ｔｈｖを定めておくことで、収穫時の生育指標を簡易に精度よく推定することができる。

【0063】

《第２の実施形態》
次に、第２の実施形態について説明する。上記第１の実施形態においては、肥大開始前関係式（一次式）と肥大開始後関係式（一次式）の２つの式を用いて、収穫時のりん茎の球径を特定（推定）することとしたが、本第２の実施形態では、１つの指数関数式を用いて収穫時のりん茎の球径を特定する点が第１の実施形態と異なる。なお、本第２の実施形態のサーバ１０は、図１１に示すように、図３の実測値取得時期判定部２２、肥大開始前関係式特定部２４及び肥大開始後関係式特定部２６に代えて、肥大関係式特定部２５を有している点以外は、第１の実施形態と同様である。

【0064】

本発明者らは、図１２（ａ）に示すように、タマネギの葉鞘基部径及びりん茎の球径の変化（ｙ）を、横軸（ｘ）を積算気温とする片対数グラフにプロットした。その結果、直線状に近似できること、すなわち、次式（６）のような指数関数式で表現できることを見出した。
ｙ＝ｂ１・ｅｘｐ（ｃ１・ｘ） …（６）

【0065】

また、本発明者らは、それぞれの傾きは、品種や定植場所が同じであれば、定植日が多少ずれてもほとんど変化しないことを見出した。すなわち、図１２（ｂ）において、ｃ１は、品種や栽培時期、土壌などにより定まるが、切片ｂ１は実測値によって異なることとなる。

【0066】

以下、本第２の実施形態においてサーバ１０が実行する処理について、図１３のフローチャートに沿って説明する。なお、図１３の太実線で記載したステップＳ１３、Ｓ２５以外は、第１の実施形態と同様であるので、説明を省略又は簡略化するものとする。

【0067】

図１３の処理においては、まずステップＳ１０において、入力受付部２０が、利用者端末７０から入力があるまで待機する。利用者端末７０から入力があり、ステップＳ１２に移行すると、入力受付部２０は、利用者端末７０から入力された情報を取得し、取得した情報に対応する定数を取得する。例えば、入力受付部２０は、図１２（ｃ）の定数テーブル４２を参照して、取得した圃場の位置情報や、品種の情報、定植日の情報に対応する定数ｃ１、Ｔｈｖを取得する。図１２（ｃ）の定数ｃ１は、上式（６）の傾きｃ１であり、Ｔｈｖは、予め定められた収穫時の積算気温である。

【0068】

次いで、ステップＳ１３では、肥大関係式特定部２５が、上式（６）のｘに測定日における定植時からの積算気温、ｙに実測値（葉鞘基部径又はりん茎の球径）を代入することで切片ｂ１を求め、肥大関係式を特定する。

【0069】

例えば、ある測定日における実測値ａｌｙがａ１ｙ＝３１．７６であり、気温ＤＢ４０の過去データに基づいて算出される測定日における定植日からの積算気温ａｌｘが、ａｌｙ＝１０４９．７であったとする。この場合、肥大関係式特定部２５は、上式（６）のｙにａｌｙ、ｘにａｌｘを代入し（1049.7＝ｂ１・ｅｘｐ（0.0016398251×31.76））、切片ｂ１を求め（ｂ１＝5.68021）、肥大関係式として次式（６）’を特定する。
ｙ＝5.68021・ｅｘｐ(0.001639825x) …（６）’

【0070】

次いで、ステップＳ２５では、収穫時球径特定部２８が、肥大関係式（上式（６）’）を用いて収穫時のりん茎の球径を特定（推定）する。この場合、収穫時球径特定部２８は、定植時からの積算気温が図１１（ｃ）のＴｈｖ＝１７００に到達するときを収穫時（収穫適期）とみなし、肥大関係式（上式（６）’）にｘ＝Ｔｈｖ＝１７００を代入することで、収穫時におけるりん茎の球径（ｙ≒９２ｍｍ）を特定（推定）する。

【0071】

その後の処理（Ｓ２６、Ｓ２８，Ｓ３０，Ｓ３２）は、上記第１の実施形態と同様である。

【0072】

以上説明したように、本第２の実施形態によると、入力受付部２０が、タマネギの葉鞘基部径又はりん茎の球径の実測値と測定日（実測時期）の情報を取得し（Ｓ１２）、肥大関係式特定部２５が、測定日の定植時からの積算気温（実測値に対応する積算気温）を特定し、実測値と実測値に対応する積算気温とを用いて、肥大関係式（指数関数式）を特定する（Ｓ１３）。そして、収穫時球径特定部２８は、肥大関係式を用いて、予め定められた収穫時の積算気温（Ｔｈｖ）になるときのりん茎の球径を特定する（Ｓ２５）。このように、本第２の実施形態では、少なくとも１回の葉鞘基部径又はりん茎の球径の実測値を用いることで、簡易に、収穫時のりん茎重を精度よく推定することができる。また、図１２（ｃ）に示すような定数を、様々な品種のタマネギを各種条件下で栽培したときに求め、定数テーブル４２に格納しておくことで、様々な品種のタマネギを様々な条件下で栽培した場合でも、収穫時の球径や重量、収量や収穫日を簡易に算出することができる。

【0073】

なお、サーバ１０が、上記第１の実施形態の処理と第２の実施形態の処理のいずれを実行するかは、生産者等が選択できるようにしてもよい。

【0074】

（変形例）
なお、上記第２の実施形態では、サーバ１０は、タマネギの葉鞘基部径又はりん茎の球径の実測値と、測定日における定植時からの積算気温の情報と、に基づいて、収穫時のりん茎の球径を特定（推定）する場合について説明した。しかしながら、これに限られるものではなく、第１の実施形態の変形例と同様、収穫時におけるタマネギに関する生育指標（例えば、乾物重、新鮮重、草丈、葉面積など）を特定（推定）することも可能である。

【0075】

例えば、図１４に示すように、タマネギの乾物重の変化（ｙ）は、横軸（ｘ）を積算気温とする片対数グラフにプロットすると、直線状に近似することができる。

【0076】

したがって、収穫時の乾物重についても、上記第２の実施形態と同様に、乾物重の実測値と、実測値の測定日における定植時からの積算気温とに基づいて、肥大関係式を求め、求めた肥大関係式から収穫時（積算気温がＴｈｖのとき）の乾物重を特定（推定）することができる。これにより、収穫時の乾物重を簡易かつ精度よく推定することができる。

【0077】

なお、タマネギに関するその他の生育指標（新鮮重、草丈、葉面積など）についても、乾物重の場合と同様である。各生育指標の積算気温に対する変化は、図１４と同様に変化するため、定数ｃ１、Ｔｈｖを定めておくことで、収穫時の生育指標を簡易かつ精度よく推定することができる。

【0078】

なお、上記実施形態では、サーバ１０が図３の各部の機能を有する場合について説明したが、これに限られるものではない。例えば、利用者端末７０が図３の各部の機能を有していてもよい。この場合、利用者端末７０は、ネットワーク８０に接続されていなくてもよい。また、利用者端末７０は、ネットワーク８０を介してサーバ１０に格納されている各種データを取得して、図５の処理を実行するようにしてもよい。

【0079】

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。

【符号の説明】

【0080】

１０ サーバ
２０ 入力受付部
２２ 実測値取得時期判定部
２４ 肥大開始前関係式特定部
２５ 肥大関係式特定部
２６ 肥大開始後関係式特定部
２８ 収穫時球径特定部
３０ 収穫時りん茎重特定部
３２ 収量推定部
３４ 収穫日特定部
３６ 出力部
９０ ＣＰＵ（コンピュータ）
１００ 収穫情報予測システム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】