(57)【特許請求の範囲】
【請求項１】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域についての前記補正値を導出する際、
Ｍ個の領域について前記第１指標値取得部により取得された前記第１実際値と、前記第２指標値取得部により取得された前記第２実際値との相関関係を分析し、分析結果に基づいて、前記第１指標の指標値を入力とし前記第２指標の指標値を出力とする第２計算式を生成し、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値から前記第２計算式により前記補正値を導出する
ことを特徴とする農作物関連値導出装置。
【請求項２】
Ｍ個の領域は、Ｎ個の領域のうち、農作物の生育量の相違に着目して選択された複数の領域であることを特徴とする請求項１に記載の農作物関連値導出装置。
【請求項３】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第２指標値取得部は、葉緑素計を用いた計測により得られたＳＰＡＤの指標値の実際の値、地上における近傍からのセンシングにより得られるＮＤＶＩの指標値の実際の値、または、葉色板を用いた地上での計測により得られる葉色板値の実際の値を前記第２指標の指標値として取得することを特徴とする農作物関連値導出装置。
【請求項４】
前記農作物関連値導出部は、前記農作物関連値として、所定のタイミングにおいて供給する肥料、堆肥、土壌改良資材、植物生育調整剤或いは薬剤の量、収穫物の収量、農作物の現在或いは将来の品質、農作物が所定の状態となると予測される時期、農作物についての被害の状況、および、農作物に影響を与える環境の状態の何れか１つを少なくとも導出する
ことを特徴とする請求項１から３の何れか１項に記載の農作物関連値導出装置。
【請求項５】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、追肥時の施肥量であり、
前記農作物関連値導出部は、蓄積された過去の実際の追肥時の施肥量および前記第２指標の指標値の組み合わせと、農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および前記第２指標の指標値を説明変数とする重回帰分析を行って求められた重回帰式に基づく計算式を前記第１計算式として用いて追肥時の施肥量を導出する
ことを特徴とする農作物関連値導出装置。
【請求項６】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての収穫物の収量であり、
前記第１指標値取得部は、Ｎ個の領域について、出穂前の所定のタイミングの前記第１実際値を取得し、
前記第２指標値取得部は、Ｍ個の領域について、出穂前の所定のタイミングの前記第２実際値を取得する
ことを特徴とする農作物関連値導出装置。
【請求項７】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての現在或いは将来のタンパク質の状態であり、
前記第２指標値取得部は、Ｍ個の領域について、農作物の葉について測定された葉緑素に関する指標の指標値の実際の値を前記第２実際値として取得し、
前記農作物関連値導出部は、Ｎ個の領域について、ＳＰＡＤの指標値から、稲についての現在或いは将来のタンパク質の状態を前記農作物関連値として導出する前記第１計算式を用いて、稲についての現在或いは将来のタンパク質の状態を導出することを特徴とする農作物関連値導出装置。
【請求項８】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての将来の籾水分或いは玄米水分の状況であり、
前記第２指標値取得部は、Ｍ個の領域について、籾水分或いは玄米水分に係る指標の指標値の実際の値を前記第２実際値として取得し、
前記農作物関連値導出部は、Ｎ個の領域について、籾水分或いは玄米水分に係る指標の指標値から、稲についての将来の籾水分或いは玄米水分の状況を前記農作物関連値として導出する前記第１計算式を用いて、稲についての将来の籾水分或いは玄米水分の状況を導出する
ことを特徴とする農作物関連値導出装置。
【請求項９】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、将来における農作物の被害の状況であり、
前記第２指標値取得部は、Ｍ個の領域について、農作物の被害の状況に係る指標の指標値の実際の値を前記第２実際値として取得し、
前記農作物関連値導出部は、Ｎ個の領域について、農作物の被害の状況に係る指標の指標値から、将来における農作物の被害の状況を前記農作物関連値として導出する前記第１計算式を用いて、将来における農作物の被害の状況を導出する
ことを特徴とする農作物関連値導出装置。
【請求項１０】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値導出部は、前記農作物関連値として、農作物の現在或いは将来の品質に係る前記農作物関連値を導出し、
農作物の現在或いは将来の品質に係る前記農作物関連値には、現在或いは将来における農作物の栄養素の状態、および、現在或いは将来における農作物の含有水分の状況の何れかが含まれる
ことを特徴とする農作物関連値導出装置。
【請求項１１】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値導出部は、前記農作物関連値として、農作物が所定の状態となると予測される時期に係る前記農作物関連値を導出し、
農作物が所定の状態となると予測される時期に係る前記農作物関連値には、農作物が収穫可能な状態となる時期が含まれる
ことを特徴とする農作物関連値導出装置。
【請求項１２】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値導出部は、前記農作物関連値として、農作物についての被害の状況に係る前記農作物関連値を導出し、
農作物についての被害の状況に係る前記農作物関連値には、現在或いは将来における倒伏の状況、現在或いは将来における枯れ被害の状況、現在或いは将来における病害の状況、および、現在或いは将来における害虫による被害の状況の何れかが含まれる
ことを特徴とする農作物関連値導出装置。
【請求項１３】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得する第１指標値取得部と、
Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第２指標値取得部と、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する農作物関連値導出部とを備え、
前記農作物関連値導出部は、
Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、
前記第１指標値取得部により取得された当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値導出部は、前記農作物関連値として、農作物に影響を与える環境の状態に係る前記農作物関連値を導出し、
農作物に影響を与える環境の状態に係る前記農作物関連値には、農作物が栽培される領域における病害虫の発生の程度、および、圃場における雑草の発生の程度の何れかが含まれる
ことを特徴とする農作物関連値導出装置。
【請求項１４】
前記第１指標値取得部は、ドローンを用いた上空からのセンシングにより得られるＮＤＶＩの指標値の実際の値、または、衛星写真に基づくＮＤＶＩの指標値の実際の値を前記第１実際値として取得することを特徴とする請求項１から１３の何れか１項に記載の農作物関連値導出装置。
【請求項１５】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域についての前記補正値を導出する際、Ｍ個の領域についての前記第１実際値と前記第２実際値との相関関係を分析し、分析結果に基づいて、前記第１指標の指標値を入力とし前記第２指標の指標値を出力とする第２計算式を生成し、当該一の領域についての前記第１実際値から前記第２計算式により前記補正値を導出する
ことを特徴とする農作物関連値導出方法。
【請求項１６】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第１ステップにおいて、葉緑素計を用いた計測により得られたＳＰＡＤの指標値の実際の値、地上における近傍からのセンシングにより得られるＮＤＶＩの指標値の実際の値、または、葉色板を用いた地上での計測により得られる葉色板値の実際の値を前記第２指標の指標値として取得する
ことを特徴とする農作物関連値導出方法。
【請求項１７】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、追肥時の施肥量であり、
前記第２ステップにおいて、蓄積された過去の実際の追肥時の施肥量および前記第２指標の指標値の組み合わせと、農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および前記第２指標の指標値を説明変数とする重回帰分析を行って求められた重回帰式に基づく計算式を前記第１計算式として用いて追肥時の施肥量を導出する
ことを特徴とする農作物関連値導出方法。
【請求項１８】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての収穫物の収量であり、
前記第１ステップにおいて、Ｎ個の領域について、出穂前の所定のタイミングの前記第１実際値を取得すると共に、Ｍ個の領域について、出穂前の所定のタイミングの前記第２実際値を取得する
ことを特徴とする農作物関連値導出方法。
【請求項１９】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての現在或いは将来のタンパク質の状態であり、
前記第１ステップにおいて、Ｍ個の領域について、農作物の葉について測定された葉緑素に関する指標の指標値の実際の値を前記第２実際値として取得し、
前記第２ステップにおいて、Ｎ個の領域について、ＳＰＡＤの指標値から、稲についての現在或いは将来のタンパク質の状態を前記農作物関連値として導出する前記第１計算式を用いて、稲についての現在或いは将来のタンパク質の状態を導出することを特徴とする農作物関連値導出方法。
【請求項２０】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、稲についての将来の籾水分或いは玄米水分の状況であり、
前記第１ステップにおいて、Ｍ個の領域について、籾水分或いは玄米水分に係る指標の指標値の実際の値を前記第２実際値として取得し、
前記第２ステップにおいて、Ｎ個の領域について、籾水分或いは玄米水分に係る指標の指標値から、稲についての将来の籾水分或いは玄米水分の状況を前記農作物関連値として導出する前記第１計算式を用いて、稲についての将来の籾水分或いは玄米水分の状況を導出する
ことを特徴とする農作物関連値導出方法。
【請求項２１】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記農作物関連値は、将来における農作物の被害の状況であり、
前記第１ステップにおいて、Ｍ個の領域について、農作物の被害の状況に係る指標の指標値の実際の値を前記第２実際値として取得し、
前記第２ステップにおいて、Ｎ個の領域について、農作物の被害の状況に係る指標の指標値から、将来における農作物の被害の状況を前記農作物関連値として導出する前記第１計算式を用いて、将来における農作物の被害の状況を導出する
ことを特徴とする農作物関連値導出方法。
【請求項２２】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第２ステップにおいて、前記農作物関連値として、農作物の現在或いは将来の品質に係る前記農作物関連値を導出し、
農作物の現在或いは将来の品質に係る前記農作物関連値には、現在或いは将来における農作物の栄養素の状態、および、現在或いは将来における農作物の含有水分の状況の何れかが含まれる
ことを特徴とする農作物関連値導出方法。
【請求項２３】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第２ステップにおいて、前記農作物関連値として、農作物が所定の状態となると予測される時期に係る前記農作物関連値を導出し、
農作物が所定の状態となると予測される時期に係る前記農作物関連値には、農作物が収穫可能な状態となる時期が含まれる
ことを特徴とする農作物関連値導出方法。
【請求項２４】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第２ステップにおいて、前記農作物関連値として、農作物についての被害の状況に係る前記農作物関連値を導出し、
農作物についての被害の状況に係る前記農作物関連値には、現在或いは将来における倒伏の状況、現在或いは将来における枯れ被害の状況、現在或いは将来における病害の状況、および、現在或いは将来における害虫による被害の状況の何れかが含まれる
ことを特徴とする農作物関連値導出方法。
【請求項２５】
農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出方法であって、
前記農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値の実際の値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、前記領域の広域に比して限定された狭い範囲である領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値の実際の値である第２実際値を取得する第１ステップと、
Ｎ個の領域について、前記第２指標の指標値から前記農作物関連値を導出する計算式であって、前記農作物関連値と前記第２指標の指標値との関係に基づいて導出された計算式である第１計算式を用いて前記農作物関連値を導出する第２ステップとを含み、
前記第２ステップにおいて、Ｎ個の領域のうち前記第２実際値が取得されていない一の領域の前記農作物関連値を導出する際、当該一の領域についての前記第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した前記補正値から前記第１計算式により前記農作物関連値を導出し、
前記第２ステップにおいて、前記農作物関連値として、農作物に影響を与える環境の状態に係る前記農作物関連値を導出し、
農作物に影響を与える環境の状態に係る前記農作物関連値には、農作物が栽培される領域における病害虫の発生の程度、および、圃場における雑草の発生の程度の何れかが含まれる
ことを特徴とする農作物関連値導出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、農作物に関連する農作物関連値を導出する農作物関連値導出装置および農作物関連値導出方法に関する。
【背景技術】
【０００２】
従来、圃場で栽培される農作物について、圃場に対する遠方からの観測により得られた所定の指標の指標値と、追肥時の適切な施肥量との間に相関関係があることが分かっており、圃場に対して追肥するにあたって、圃場に対する遠方からの観測により所定の指標の指標値の実際の値（実際値）を取得し、当該相関関係を反映した計算式を用いて、取得した実際値から適切な施肥量を求めることが行われている。例えば、ドローンや衛星を用いた上空からのセンシングによりＮＤＶＩに係る実際値を取得し、これを用いて施肥量を求めることが行われている。また、特許文献１には、圃場に対する遠方からの観測により得られる実際値に関して、圃場全体が撮影可能な程度に圃場から遠方に設けられたカメラにより圃場全体を撮影し、撮影結果に基づいて圃場の作物情報（実際値）を求める技術が記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−１３６４１１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
実際値を利用して追肥時の適切な施肥量を高い精度で求めるためには、実際値の精度が重要になる。圃場に対する遠方からの観測によって実際値を得る場合、短時間の作業で多くの圃場についての実際値を得ることができる。しかしながら、天候、空気の状態等の環境による影響を受けやすく、得られた実際値の精度について改善の余地がある。これを改善するために、葉緑素計を用いた農作物の計測等、農作物の近傍で農作物に直接アクセスして農作物を観測し、実際値を得るようにした場合、実際値の精度は向上するものの、多くの圃場についての実際値を得るためには、作業に多大な時間を要する結果となってしまう。
【０００５】
以上の問題は、農作物を栽培する農作物栽培領域を遠方から観測して実際値を取得し、この実際値に基づいて農作物に関連する農作物関連値を導出する技術において、広く発生し得る。農作物関連値は、収穫物の収量、農作物の現在或いは将来の品質、農作物が所定の状態となると予測される時期、農作物についての被害の状況、および、農作物に影響を与える環境の状態等である。
【０００６】
本発明は、このような問題を解決するために成されたものであり、農作物関連値を求めるために必要な実際値を取得するための作業にかかる時間を短くしつつ、高い精度で農作物関連値を求めることができるようにすることを目的とする。
【課題を解決するための手段】
【０００７】
上記した課題を解決するために、本発明では、農作物栽培領域を構成するＮ個の領域について、領域の広域を遠方から観測することにより得られた第１指標の指標値である第１実際値を取得すると共に、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域について、農作物を計測器で計測することにより、領域の狭域を近傍から観測することにより、または、農作物を現場で観察することにより得られた第２指標の指標値である第２実際値を取得し、Ｎ個の領域について、第２指標の指標値から農作物関連値を導出する第１計算式を用いて農作物関連値を導出する。そして第２実際値が取得されていない一の領域の農作物関連値の導出に際し、当該一の領域についての第１実際値を、Ｍ個の領域についての第１実際値と第２実際値との関係に基づいて、第１指標の指標値から第２指標の指標値へと変換されるように補正することによって補正値を導出し、導出した補正値から第１計算式により農作物関連値を導出するようにしている。
【発明の効果】
【０００８】
上記のように構成した本発明によれば、第１指標の指標値である第１実際値のみが取得された一の領域について農作物関連値を導出するにあたって、第１実際値がそのまま用いられて農作物関連値が導出されるのではなく、以下の処理が行われる。すなわち、まず、第１実際値および第２実際値が取得された他の領域についての第１実際値と第２実際値との関係に基づいて、当該一の領域に係る第１実際値が第１指標の指標値から第２指標の指標値へと変換されるように補正されて補正値が導出される。このようにして導出された補正値は、第２指標の指標値に相当し、第１指標の指標値たる第１実際値よりも精度が高い。そしてこの補正値が用いられて農作物関連値が導出される。このため、高い精度で農作物関連値を導出することができる。更に計測器を用いた農作物の計測または領域の狭域に対する近傍からの観測は、Ｎ個の全ての圃場について行われるのではなく、Ｎ個よりも少ないＭ個の圃場についてのみ行われるため、実際値を取得するための作業にかかる時間を短くすることができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の一実施形態に係る農作物関連値導出装置の機能構成例を示すブロック図である。
【図２】対象圃場を識別情報と共に示す図である。
【図３】補正用圃場の地上ＮＤＶＩと上空ＮＤＶＩとの関係を示す図である。
【図４】補正用計算式の一例によって求められる補正値を示す図である。
【図５】測定タイミングのＮＤＶＩと追肥時の施肥量との関係を示す図である。
【図６】本発明の一実施形態に係る農作物関連値導出装置の動作例を示すフローチャートである。
【発明を実施するための形態】
【００１０】
以下、本発明の一実施形態を図面に基づいて説明する。図１は、本実施形態に係る農作物関連値導出装置１の構成例を示すブロック図である。農作物関連値導出装置１は、稲作に関して、追肥（追肥についての詳細は後述）時に施す肥料の量（以下「追肥時施肥量」という）を導出し、追肥時施肥量を示す情報をユーザに提供する装置である。追肥時施肥量は、特許請求の範囲の「農作物関連値」に相当する。特に本実施形態に係る農作物関連値導出装置１は、複数の圃場を対象として複数の圃場のそれぞれについて、追肥時施肥量を導出する。ユーザは、追肥にあたって農作物関連値導出装置１により提供された追肥時施肥量を示す情報を参考にして、実際に施す肥料の量を適切に決定することができる。なお、本実施形態で「施肥量」という場合、特に説明がない限り、単位面積あたりの肥料の量を指す。
【００１１】
以下の説明では、対象とする稲、稲作について、稲の品種や、稲作が行われる環境（地域（寒冷地域や、温暖地域、乾燥地域等）、標高、圃場の規模）、稲作の方法、肥料の種類、その他の稲の生育に影響を与える外部的な要素は共通しているものとする。従って、農作物関連値導出装置１により導出された追肥時施肥量は、特定の環境で特定の方法により栽培される特定の品種の稲について、特定の種類の肥料を施す場合に適切な施肥量ということになる。
【００１２】
本実施形態において、「追肥」とは以下を意味する。すなわち一般に、ある１つの期において、稲は、田植えによって苗代から圃場に移植された後、分げつ期→幼穂形成期→出穂期→開花・受粉期→穂揃期→登熟期→成熟期という段階を経て成長する。幼穂形成期前後から出穂期前後に至るまでの期間（以下「追肥対象期間」という。）では、穂に着生する籾の個数の増大や、籾に詰まるデンプンの量の増大等を目的として、所定のタイミングで穂肥と呼ばれる肥料の追加が行われる場合がある。そして本実施形態では、「追肥」とは追肥対象期間における所定のタイミングで行われる肥料の追加を意味する。
【００１３】
図１で示すように、農作物関連値導出装置１には、液晶ディスプレイや、有機ＥＬディスプレイ等の表示装置２と、マウスやキーボード等の入力装置３とが接続されている。また図１で示すように、農作物関連値導出装置１は、機能構成として、農作物関連値導出部１０と、第１指標値取得部１１と、第２指標値取得部１２とを備えている。上記各機能ブロック１０〜１２は、ハードウェア、ＤＳＰ（Digital Signal Processor）、ソフトウェアの何れによっても構成することが可能である。例えばソフトウェアによって構成する場合、上記各機能ブロック１０〜１２は、実際にはコンピュータのＣＰＵ、ＲＡＭ、ＲＯＭなどを備えて構成され、ＲＡＭやＲＯＭ、ハードディスクまたは半導体メモリ等の記録媒体に記憶されたプログラムが動作することによって実現される。また、農作物関連値導出装置１は、記憶手段として記憶部１３を備えている。記憶部１３に記憶されるデータについては後述する。
【００１４】
本実施形態では、農作物関連値導出装置１は、ユーザから入力装置３に対して追肥時施肥量を示す情報の提供の指示（以下「情報提供指示」という）があったことをトリガとして、圃場毎の追肥時施肥量を導出し、当該情報を提供する。以下、ユーザから入力装置３に対して情報提供指示があったときの農作物関連値導出装置１の動作について詳述する。ただし以下の説明では、説明の便宜のため、農作物関連値導出装置１は、１０個の圃場を対象として追肥時施肥量を導出するものとする。また、追肥時施肥量を導出する対象の１０個の圃場を「対象圃場」という。なお、追肥時施肥量の導出に関する処理を開始するトリガは本実施形態で例示するものに限られず、例えば、予め定められた日時が到来したことをトリガとしてもよい。本実施形態では、１０個の対象圃場を含む領域が特許請求の範囲の「農作物栽培領域」に相当する。また１０個の対象圃場のそれぞれが、特許請求の範囲における「Ｎ個の領域のそれぞれ」に相当する。
【００１５】
農作物関連値導出部１０は、入力装置３に対して情報提供指示があったことを検出すると、第１指標値取得部１１に上空ＮＤＶＩ値（後述）の取得を要求する。
【００１６】
第１指標値取得部１１は、農作物関連値導出部１０から上空ＮＤＶＩ値の取得の要求があった場合、以下の処理を実行する。すなわち、まず、第１指標値取得部１１は、記憶部１３に記憶された１０個の対象圃場のそれぞれの上空データを取得する。１個の対象圃場についての上空データとは、追肥を行うタイミングの前の所定のタイミング（以下「測定タイミング」という）で、ドローンに搭載されたマルチスペクトルセンサにより圃場の稲の群落を上空（遠方）から撮影することによって生成された画像データである。なお、測定タイミングは、追肥時施肥量の導出に利用する上空データを得るための撮影を行うタイミングとして予め定められたタイミングであり、例えば、追肥を行うタイミングの直前や、出穂の３週間前、出穂の４週間前等とされる。なお、１０個の対象圃場のそれぞれには識別情報が事前に付与されており、個々の対象圃場は識別情報によって管理され、また、上空データや、後述する地上データには、対応する圃場の識別情報がメタ情報として記録されている。
【００１７】
第１指標値取得部１１は、取得した上空データのそれぞれを分析し、周知の方法で対象圃場毎のＮＤＶＩの指標値の実際の値（実際値）を導出し、取得する。以下、第１指標値取得部１１により取得された、対象圃場毎のＮＤＶＩの指標値の実際の値を「上空ＮＤＶＩ値」という。ここで、ドローンによる上空からの撮影結果に基づくＮＤＶＩ（指標としてのＮＤＶＩ）を説明の便宜上、「上空ＮＤＶＩ」と表現する。上空ＮＤＶＩは、後述する地上ＮＤＶＩとは異なる指標を意味している。上空ＮＤＶＩは、特許請求の範囲の「第１指標」に相当する。また上空ＮＤＶＩ（指標）の指標値の実際の値として、上空データに基づいて導出される上空ＮＤＶＩ値は、特許請求の範囲の「第１実際値」および「ドローンを用いた上空からのセンシングにより得られるＮＤＶＩの指標値の実際の値」に相当する。
【００１８】
ここで本実施形態では、上空ＮＤＶＩ値の導出にあたって使用される上空データは、ドローンに搭載されたマルチスペクトルセンサの撮影結果に基づいて生成されるものである。ドローンを使用することにより、１つの圃場全体にわたる広い領域を上空から撮影することが可能であり、多くの圃場のそれぞれについてスピーディに上空データを生成することができ、上空ＮＤＶＩ値を取得するための作業にかかる時間を短くすることができる。
【００１９】
このような効果は、例えば葉緑素計を用いて葉を一枚ずつ測定した結果に基づいて得られる実際値を第１指標値取得部１１が取得する値とするのではなく、圃場（領域）の広域を遠方から観測することにより得られる実際値を第１指標値取得部１１が取得する値とすることにより得ることが可能である。従って、観測（センシング）に用いる装置は、マルチスペクトルセンサに限られず、例えば、赤外線サーモグラフィであってもよい。また、第１指標値取得部１１に係る指標は、上空ＮＤＶＩに限られない。また、観測は、必ずしもドローンを用いて行われる必要はなく、塔の上部に設置された装置や、専用の設置台に設置された装置を用いて行われたりしてもよい。また、上空ＮＤＶＩ値（またはこれに代替する実際値）の導出に、衛星写真を用いるようにしてもよい。また第１指標値取得部１１が上空ＮＤＶＩ値を導出するのではなく、外部の装置により上空ＮＤＶＩ値が導出され、第１指標値取得部１１が当該外部の装置からそれを取得する構成としてもよい。また、ある１つの圃場について１つの上空データが用意されるのではなく、その１つの圃場について複数の上空データが用意され、第１指標値取得部１１が、その１つの圃場に関する複数の上空データに基づいて、その１つの圃場の上空ＮＤＶＩ値を導出する構成でもよい。また、第１指標値取得部１１が、Ｘ個の上空データから、Ｙ（Ｙ＞Ｘ）個の圃場の上空ＮＤＶＩ値を導出する構成でもよい。
【００２０】
一方、ドローンを利用して圃場を撮影し、その撮影結果に基づいて上空ＮＤＶＩ値を導出する場合や、これに準じて圃場を遠方から観測し、観測結果に基づいて上空ＮＤＶＩ値に準じた実際値を導出する場合、観測結果が天候や、空気の状態等の環境に影響を受け、導出される上空ＮＤＶＩ値（またはこれに準じた実際値）の精度が低くなる可能性がある。
【００２１】
第１指標値取得部１１は、対象圃場のそれぞれの上空ＮＤＶＩ値を取得した後、取得した上空ＮＤＶＩ値を各圃場の識別情報と対応付けて農作物関連値導出部１０に応答する。
【００２２】
また、農作物関連値導出部１０は、入力装置３に対して情報提供指示があったことを検出すると、第２指標値取得部１２に地上ＮＤＶＩ値（後述）の取得を要求する。
【００２３】
第２指標値取得部１２は、農作物関連値導出部１０から地上ＮＤＶＩ値の取得の要求があった場合、以下の処理を実行する。すなわち、第２指標値取得部１２は、記憶部１３に記憶された３個の対象圃場のそれぞれの地上データを取得する。３個の対象圃場は、１０個の対象圃場から稲の生育量の相違に着目して選択された３個の対象圃場であり、本実施形態では、生育量が最も多い圃場（以下「最大圃場」という）と、生育量が最も少ない圃場（以下「最小圃場」という）と、生育量が中間程度の圃場（以下「中間圃場」という）との３個の圃場が選択される。中間圃場は例えば、１０個の対象圃場を生育量の多さで昇順に並べたときの５番目または６番目に位置する対象圃場とされる。また、中間圃場は例えば、最大圃場の生育量と最小圃場の生育量との中間の生育量に、生育量が最も近い圃場とされる。なお、各対象圃場の生育量は、例えば既存の技術による観測によって把握され、また例えば対象圃場毎に生育量が計画的に調整されている場合には、生育量の計画に基づいて把握される。以下、最大圃場、最小圃場および中間圃場を総称して「補正用圃場」という。
【００２４】
ここで、補正用圃場の個数は、補正用圃場を含む対象圃場全体の個数よりも少ないものとされる。本実施形態では、補正用圃場の個数は３個であるため、対象圃場は４個以上とされる。ただし、後に明らかとなる通り、補正用圃場の個数に比べて対象圃場の個数が十分に多い場合に、ユーザは本実施形態に係る農作物関連値導出装置１によるメリットを効果的に得ることができる。このため、本実施形態のように、対象圃場の個数は、補正用圃場の個数よりも十分に多い方が好ましい。
【００２５】
地上データは、地上において、圃場における稲の群生の近傍で、当該群生に対する観測機器を用いた観測の観測結果に基づくデータであって、ＮＤＶＩの導出に利用可能なデータである。例えば、地上データは、人間である作業者が、マルチスペクトルセンサが搭載された可搬型の撮像装置（当然、他のセンサであってもよい）を対象となる圃場の近傍まで運び、群生を自ら撮影することによって得られたデータである。また例えば、地上データは、作業者により圃場の中央部に撮像装置を支持するための支柱が立設され、その支柱の頭頂部に取り付けた撮像装置によって稲の群生が撮影され、その撮影の撮影結果に基づいて生成されたデータである。作業者は、１０個の対象圃場のうち、どの圃場が補正用圃場であるかを認識した上で、補正用圃場のそれぞれについて上述したような作業を行って地上データを生成し、記憶部１３に記憶する。逆に言えば、作業者は、補正用圃場以外の圃場については、地上データを作成するための作業を行わない。なお上述した上空データは、圃場の広い範囲（広域）を遠方から撮影（観測）した撮影結果により得られるものであるが、地上データは、上空データに係る撮影範囲に比して限定された狭い範囲（狭域）を近傍から撮影（観測）した撮影結果により得られるものである。
【００２６】
記憶部１３に記憶された補正用圃場のそれぞれの地上データを取得した後、第２指標値取得部１２は、取得した地上データのそれぞれを分析し、周知の方法で補正用圃場毎のＮＤＶＩの指標値の実際の値を導出し、取得する。以下、第２指標値取得部１２により取得された、補正用圃場毎のＮＤＶＩの指標値の実際の値を「地上ＮＤＶＩ値」という。ここで、地上データに基づくＮＤＶＩ（指標としてのＮＤＶＩ）を説明の便宜上、「地上ＮＤＶＩ」と表現する。地上ＮＤＶＩは、上述した上空ＮＤＶＩとは明確に異なる指標を意味している。地上ＮＤＶＩは、特許請求の範囲の「第２指標」に相当する。また地上ＮＤＶＩ（指標）の指標値の実際の値として、地上データに基づいて導出される地上ＮＤＶＩ値は、特許請求の範囲の「第２実際値」、特に「領域の狭域を近傍から観測することにより得られた第２指標の指標値の実際の値」および「地上における近傍からのセンシングにより得られるＮＤＶＩの指標値の実際の値」に相当する。
【００２７】
上述した通り、地上ＮＤＶＩ値を生成する基礎となる地上データは、圃場における稲の群生の近傍で、当該群生に対する観測機器を用いた観測の観測結果に基づくデータである。従って、ドローンを用いた撮影により生成される上空データと比較して、天候、空気の状態等の影響を受けにくい。このため、地上データに基づいて導出される地上ＮＤＶＩ値は、上空データに基づいて導出される上空ＮＤＶＩ値と比較して精度が高い。
【００２８】
第２指標値取得部１２は、補正用圃場のそれぞれの地上ＮＤＶＩ値を取得した後、取得した地上ＮＤＶＩ値を各圃場の識別情報と対応付けて農作物関連値導出部１０に応答する。
【００２９】
さて、第１指標値取得部１１から１０個の対象圃場のそれぞれの上空ＮＤＶＩ値の応答を受け、第２指標値取得部１２から３個の補正用圃場のそれぞれの地上ＮＤＶＩ値の応答を受けると、農作物関連値導出部１０は、補正用計算式生成処理を実行する。以下、補正用計算式生成処理について詳述する。
【００３０】
補正用計算式生成処理において、農作物関連値導出部１０は、３個の補正用圃場のそれぞれについて、地上ＮＤＶＩ値に加え、上空ＮＤＶＩ値を取得する。次いで、農作物関連値導出部１０は、横軸に地上ＮＤＶＩの指標値を取り、縦軸に上空ＮＤＶＩの指標値を取った２次元空間に、３個の補正用圃場のそれぞれについて、地上ＮＤＶＩ値と上空ＮＤＶＩ値との組み合わせを示す点をプロットする。次いで、農作物関連値導出部１０は、プロットした各点の近似曲線を求める。農作物関連値導出部１０は、近似曲線を表す式を、上空ＮＤＶＩの指標値を入力値として入力して、地上ＮＤＶＩの指標値を出力値として出力する計算式へと変換し、これを補正用計算式（特許請求の範囲の「第２計算式」に相当）とする。なお、補正用計算式は、単純な数式だけでなく、プログラム上で定義された関数（モデル）や、プログラムから呼び出し可能な関数（モデル）を概念的に表している。このことは他の計算式についても同様である。
【００３１】
例えば、図２で示すように、１０個の圃場のそれぞれの識別情報がＨ０１、Ｈ０２、Ｈ０３・・・Ｈ１０であり、そのうち識別番号がＨ０１、Ｈ０２およびＨ０９の３個の圃場が補正用圃場であったとする。また、識別番号がＨ０１〜Ｈ１０の圃場についての上空ＮＤＶＩ値が図に示す通りであり、識別番号がＨ０１、Ｈ０２、Ｈ０９の圃場について地上ＮＤＶＩ値が図に示す通りであったとする。
【００３２】
この場合において、横軸に地上ＮＤＶＩの指標値を取り、縦軸に上空ＮＤＶＩの指標値を取った２次元空間に、識別番号がＨ０１、Ｈ０２およびＨ０９の３個の補正用圃場のそれぞれについて、地上ＮＤＶＩ値と上空ＮＤＶＩ値との組み合わせを示す点をプロットし、またその近似曲線を求めると図３で示す通りになる。そして近似曲線に基づいて求められる補正用計算式は、「出力値（地上ＮＤＶＩの指標値）＝入力値（上空ＮＤＶＩの指標値）−０．１０」となる。
【００３３】
このような方法で生成される補正用計算式に一の上空ＮＤＶＩ値を入力することによって得られる出力値は、当該一の上空ＮＤＶＩ値を、補正用圃場についての上空ＮＤＶＩ値と地上ＮＤＶＩ値との相関関係に準じて、上空ＮＤＶＩの指標値から地上ＮＤＶＩの指標値へと変換されるように補正した値ということができる。そして、上述したように地上ＮＤＶＩの指標値は、上空ＮＤＶＩの指標値よりも精度が高く、補正用圃場についての上空ＮＤＶＩと地上ＮＤＶＩとの相関関係に準じて補正された出力値は、入力値よりも精度が高い値ということができる。以上の通り、農作物関連値導出部１０は、補正用圃場（Ｍ個の領域）についての上空ＮＤＶＩ値と地上ＮＤＶＩ値との相関関係を分析し、分析結果に基づいて、上空ＮＤＶＩ（第１指標）の指標値を入力とし地上ＮＤＶＩ（第２指標）の指標値を出力とする補正用計算式（第２計算式）を生成している。
【００３４】
ここで、上述したように本実施形態では、補正用圃場として、生育量が最も多い圃場と生育量が最も少ない圃場と生育量が中間程度の圃場との３個の圃場が選択される。これは以下の理由による。すなわち、周知の通り、圃場における生育量は、圃場についてのＮＤＶＩの指標値に影響を与え、生育量とＮＤＶＩの指標値とは正の相関関係がある。そして、上記の方針で補正用圃場を選択することにより、横軸に地上ＮＤＶＩの指標値を取り、縦軸に上空ＮＤＶＩの指標値を取った２次元空間で適度に離間した３つの点をプロットすることが可能となり、各点の位置に基づいて導出される近似曲線の妥当性を高めることができるからである。
【００３５】
なお、本実施形態では、補正用圃場を３つとしているが、これは説明の単純化のためであり、当然、４つ以上の圃場を補正用圃場としてもよい。すなわち、補正用圃場の個数は、高い精度で近似曲線を導出できるのに十分な個数とされる。補正用圃場の個数が４つ以上の場合も、生育量が分散した圃場を補正用圃場として選択することにより、上述した効果を奏することができる。
【００３６】
以上の処理が補正用計算式生成処理である。補正用計算式生成処理により補正用計算式を生成した後、農作物関連値導出部１０は、以下の処理を実行する。すなわち、補正用圃場ではない７個の圃場（以下「被補正圃場」という）のそれぞれについて、補正用計算式に上空ＮＤＶＩ値を入力し、その出力値を補正値として取得する。
【００３７】
例えば、図２を参照し、７個の被補正圃場について、上空ＮＤＶＩ値が図２の通りであるものとし、補正用計算式が「出力値＝入力値−０．１０」であるとする。この場合、これら７個の被補正圃場については、補正値は図４に示す通りとなる。識別番号がＨ０３の被補正圃場に着目すると、当該被補正圃場の上空ＮＤＶＩ値は「０．５０」であり、これを入力値として補正用計算式に入力することによって出力される出力値は、「０．５０−０．１０＝０．４０」となる。
【００３８】
以上のようにして、農作物関連値導出部１０は、３個の補正用圃場について地上ＮＤＶＩ値を得て、更に７個の被補正圃場について補正値を得る。そして、農作物関連値導出部１０は、３個の補正用圃場については、施肥量導出計算式（内容については後述）に地上ＮＤＶＩ値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出する。一方、農作物関連値導出部１０は、７個の被補正圃場については、施肥量導出計算式に補正値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出する。
【００３９】
施肥量導出計算式は、地上ＮＤＶＩの指標値を入力し、施肥量を出力する計算式である。施肥量導出計算式に一の地上ＮＤＶＩの指標値を入力したときに出力される施肥量は、測定タイミングの地上ＮＤＶＩの指標値の実際の値が、当該計算式に入力された地上ＮＤＶＩの指標値（と同じ値）であった場合に、所定のカテゴリについての農作物（農作物の収穫物も含む）の状態を目標状態とするために必要かつ適切な追肥時施肥量である。本実施形態において、所定のカテゴリについての農作物の状態は例えば、農作物の外観的あるいは食味的な品質や、収穫物の収量、籾数、追肥後の生育過程の指標等である。追肥後の生育過程の指標は、例えば、追肥した後の穂揃期の葉緑素計値（ＳＰＡＤ値）である。なお、穂揃期の葉緑素計値（ＳＰＡＤ値）が一定値以上の場合は、品質が良いということができる。
【００４０】
例えば、施肥量導出計算式の出力は、収穫物の収量（単位当たりの収量）を○○ｋｇとするために必要かつ適切な追肥時施肥量である。この場合、収穫物の収量が「所定のカテゴリについての農作物の状態」に相当し、収量が○○ｋｇであるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって、収穫物の収量が○○ｋｇとなることが期待される。また例えば、施肥量導出計算式の出力は、収穫物全体に対する低外観品質の収穫物の割合（以下「低品質割合」という）を○○％とするために必要かつ適切な追肥時施肥量である。なお、低外観品質の収穫物は、例えば、白未熟粒（腹白粒や、背白粒、基部未熟粒、心白粒、乳白粒等）や、その他の未熟粒、被害粒、死米、着色粒等とである。この場合、低品質割合が「所定のカテゴリについての農作物の状態」に相当し、低品質割合が○○％であるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって、低品質割合が○○％となることが期待される。なお、低品質割合に代えて、収穫物の「整粒」の割合としてもよいことは勿論である。
【００４１】
図５は、測定タイミングにおける地上ＮＤＶＩの指標値を縦軸とし、追肥時の施肥量を横軸とする２次元空間に描画された１次式によって、上述した低品質割合（所定のカテゴリについての農作物の状態）を５％（目標状態）とするときの、測定タイミングにおける地上ＮＤＶＩの指標値と追肥時の施肥量との関係を単純化して示している。図５で示すように、低品質割合を５％とするときの追肥時の施肥量と、測定タイミングにおけるＮＤＶＩとの間には、地上ＮＤＶＩ値が大きくなるほど施肥量が小さくなる関係がある。このように、所定のカテゴリについての農作物の状態を目標状態とするときの、測定タイミングにおける地上ＮＤＶＩと追肥時の施肥量とには強い相関関係がある。
【００４２】
本実施形態では、施肥量導出計算式は、蓄積された過去の実際の追肥時の施肥量および測定タイミングにおける地上ＮＤＶＩの指標値の組み合わせと、所定のカテゴリについての農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングにおける地上ＮＤＶＩの指標値を説明変数とする重回帰分析を行って求められた重回帰式に基づいて求められる。
【００４３】
詳述すると、事前に、実験や既存の圃場に対する観察により、追肥時の施肥量および測定タイミングの地上ＮＤＶＩの指標値の実際の値の組み合わせと、所定のカテゴリについての農作物の状態との関係を示すデータが蓄積される。つまり、追肥時の施肥量が○○であり、測定タイミングの地上ＮＤＶＩの指標値の実際の値が○○であった場合に、所定のカテゴリについての農作物の状態は○○となったというデータが蓄積される。上述した通り、所定のカテゴリについての農作物の状態は、様々な観点から定義可能である。
【００４４】
そして、蓄積されたデータについて、所定のカテゴリについての農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングの地上ＮＤＶＩの指標値を説明変数とする重回帰分析が行われ、農作物の状態を目標状態とするための追肥時の施肥量と測定タイミングの地上ＮＤＶＩの指標値との関係が定義された重回帰式が生成される。例えば、農作物の状態が上述した低品質割合である場合には、低品質割合を目標状態（一例として５％）とするための、追肥時の施肥量と測定タイミングの地上ＮＤＶＩの指標値との関係が定義された重回帰式が生成される。そして重回帰式に基づいて、地上ＮＤＶＩの指標値を入力とし、重回帰式により表される地上ＮＤＶＩの指標値と追肥時の施肥量との関係から導かれる施肥量（（本例では、低品質割合を目標状態とするための施肥量）を出力とする施肥量導出計算式が生成される。このようにして生成される施肥量導出計算式の出力は、測定タイミングの地上ＮＤＶＩの指標値が入力された値である場合に農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量に相当する。なお重回帰式を、稲の生育時期毎、例えば出穂前３０日、２５日、２０日、１５日、１０日、５日といったように作成するようにしてもよい。この場合において、生育診断する時期が出穂前何日なのかを判定するには、作付け地点、品種、移植日、苗の葉齢及び苗姿、気温や、日長時間等から推定したり、幼穂の長さから推定したりすることができる。
【００４５】
なお、本実施形態の施肥量導出計算式の内容は、あくまで単純化した一例である。施肥量導出計算式は、測定タイミングの地上ＮＤＶＩの指標値を入力の１つとし、農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量を出力する計算式であれば、どのようなものであってもよい。一例として、測定タイミングの地上ＮＤＶＩの指標値以外に、農作物の状態に影響を与える要素の要素値を入力とし、当該要素を加味して追肥時の施肥量を出力する式であってもよい。当該要素は、一例として、特定の期間の気温（到来していない期間（例えば登熟期）については予想気温）や、特定の期間における雨量等である。また、施肥量導出計算式は重回帰分析により生成されるものである必要はなく、施肥量導出計算式の生成には様々な機械学習（一例としてニューラルネットワークを用いた機械学習）を応用可能である。
【００４６】
ここで施肥量導出計算式を利用して追肥時施肥量を高い精度で求めるためには、入力値の精度が重要になる。上述したように上空ＮＤＶＩ値を導出するために用いる上空データは、ドローンを利用した圃場に対する遠方からの観測によって得ることが可能であるため、短時間の作業で多くの圃場について得ることができる。しかしながら、撮影結果が天候、空気の状態等の環境による影響を受けやすく、得られた上空ＮＤＶＩ値の精度について改善の余地がある。
【００４７】
これを改善するためには、全ての圃場について地上ＮＤＶＩ値を導出するようにすればよいが、上述したように、地上ＮＤＶＩ値の導出に用いられる地上データを準備するためには、作業者が機材をもって圃場の近傍に赴き、撮影装置その他のセンサを使用して圃場の群生を観測するといった作業や、圃場の中央部に撮像装置が取り付けられた支柱を立設する等の作業を１つ１つの圃場について行う必要があり、作業に非常に時間がかかる。
【００４８】
一方で本実施形態では、被補正圃場については、被補正圃場に係る上空ＮＤＶＩ値が単純に施肥量導出計算式に入力されるのではなく（この場合には、当然、上空ＮＤＶＩの指標値を入力し、追肥時の施肥量を導出する施肥量導出式が用意されることになる）、以下の処理が行われる。すなわち、被補正圃場に係る上空ＮＤＶＩ値が３個の補正用圃場に係る上空ＮＤＶＩ値と地上ＮＤＶＩ値との関係に基づいて上述した方法で補正され、その精度が高められた上で施肥量導出計算式に入力される。このため、施肥量導出計算式を利用して導出される追肥時施肥量補正量の精度を高めることができる。更に、地上での観測が行われる圃場を１０個ではなく３個に限定できるため、作業に要する時間を短縮することができる。
【００４９】
さて１０個の対象圃場のそれぞれについて最終的な追肥時施肥量を導出した後、農作物関連値導出部１０は、表示装置２を制御して、１０個の対象圃場のそれぞれについて識別情報と対応付けて、導出した追肥時施肥量を示す情報を表示装置２の表示領域に表示する。ユーザは、表示装置２を参照することにより、１０個の対象圃場のそれぞれについて農作物関連値導出装置１が決定した追肥時施肥量を認識でき、追肥の際に実際に供給する施肥量を決定するときの参考とすることができる。
【００５０】
次に本実施形態に係る農作物関連値導出装置１の動作について図６のフローチャートを用いて説明する。図６（Ａ）は農作物関連値導出部１０の動作を示し、図６（Ｂ）は第１指標値取得部１１の動作を示し、図６（Ｃ）は第２指標値取得部１２の動作を示している。以下の説明では、農作物関連値導出装置１は、Ｎ（ただし、Ｎ≧２）個の圃場を対象として追肥時施肥量を決定するものとし、Ｎ個の対象圃場のうち、Ｍ（ただし、Ｍ＜Ｎ）個の圃場が補正用圃場であるものとする。
【００５１】
図６（Ａ）で示すように、農作物関連値導出部１０は、情報提供指示があったか否かを監視する（ステップＳＡ１）。情報提供指示があった場合（ステップＳＡ１：ＹＥＳ）、農作物関連値導出部１０は、第１指標値取得部１１に上空ＮＤＶＩ値の取得を要求する（ステップＳＡ２）。更に農作物関連値導出部１０は、第２指標値取得部１２に地上ＮＤＶＩ値の取得を要求する（ステップＳＡ３）。ステップＳＡ２とステップＳＡ３の処理の順番は順不同であり、また、並行して実行されてもよい。
【００５２】
図６（Ｂ）で示すように、第１指標値取得部１１は、農作物関連値導出部１０から上空ＮＤＶＩ値の取得の要求があったか否かを監視する（ステップＳＢ１）。当該要求があった場合（ステップＳＢ１：ＹＥＳ）、第１指標値取得部１１は、記憶部１３に記憶されたＮ個の対象圃場の上空データを取得する（ステップＳＢ２）。次いで、第１指標値取得部１１は、ステップＳＢ２で取得した上空データのそれぞれを分析し、対象圃場毎の上空ＮＤＶＩ値を取得する（ステップＳＢ３）。次いで、第１指標値取得部１１は、ステップＳＢ３で取得した対象圃場毎の上空ＮＤＶＩ値のそれぞれを圃場の識別情報と対応付けて、農作物関連値導出部１０に応答する（ステップＳＢ４）。
【００５３】
図６（Ｃ）で示すように、第２指標値取得部１２は、農作物関連値導出部１０から地上ＮＤＶＩ値の取得の要求があったか否かを監視する（ステップＳＣ１）。当該要求があった場合（ステップＳＣ１：ＹＥＳ）、第２指標値取得部１２は、記憶部１３に記憶されたＭ個の補正用圃場のそれぞれの地上データを取得する（ステップＳＣ２）。次いで、第２指標値取得部１２は、ステップＳＣ２で取得した地上データのそれぞれを分析し、Ｍ個の補正用圃場毎の地上ＮＤＶＩ値を取得する（ステップＳＣ３）。次いで、第２指標値取得部１２は、ステップＳＣ３で取得した補正用圃場の地上ＮＤＶＩ値のそれぞれを圃場の識別情報と対応付けて、農作物関連値導出部１０に応答する（ステップＳＣ４）。
【００５４】
図６（Ａ）で示すように、ステップＳＡ３の処理後、農作物関連値導出部１０は、第１指標値取得部１１からの上空ＮＤＶＩ値の応答、および、第２指標値取得部１２からの地上ＮＤＶＩ値の応答の双方があったか否かを監視する（ステップＳＡ４）。双方の応答があった場合（ステップＳＡ４：ＹＥＳ）、農作物関連値導出部１０は、補正用計算式生成処理を実行し、補正用計算式を生成する（ステップＳＡ５）。上述したように補正用計算式生成処理では、補正用圃場についての上空ＮＤＶＩ値と地上ＮＤＶＩ値との相関関係が分析され、上空ＮＤＶＩの指標値を入力とし、相関関係を反映して補正された補正値を出力とする補正用計算式が生成される。
【００５５】
次いで、農作物関連値導出部１０は、「Ｎ−Ｍ」個の被補正圃場のそれぞれについて、ステップＳＡ５で生成した補正用計算式に上空ＮＤＶＩ値を入力することによって、補正値を得る（ステップＳＡ６）。次いで、農作物関連値導出部１０は、施肥量導出計算式を用いてＮ個の対象圃場のそれぞれについて最終的な追肥時施肥量を導出する（ステップＳＡ７）。詳述すると農作物関連値導出部１０は、補正用圃場については、施肥量導出計算式に地上ＮＤＶＩ値を入力して得られる出力値を最終的な追肥時施肥量として導出し、一方、被補正圃場については、施肥量導出計算式に補正値を入力して得られる出力値を最終的な追肥時施肥量として導出する。次いで、農作物関連値導出部１０は、導出した追肥時施肥量を示す情報を表示装置２に表示する（ステップＳＡ８）。
【００５６】
以上詳しく説明したように、農作物関連値導出装置１は、Ｎ（ただし、Ｎ≧２）個の対象圃場のそれぞれについて、遠方からの観測により得られた上空ＮＤＶＩ値を取得すると共に、Ｎ個の圃場のうち、Ｍ（ただし、Ｍ＜Ｎ）個の補正用圃場のそれぞれについて、地上ＮＤＶＩ値を取得する。そして農作物関連値導出装置１は、被補正圃場の追肥時施肥量を導出する際、被補正圃場についての上空ＮＤＶＩ値を、３個の補正用圃場の上空ＮＤＶＩ値と地上ＮＤＶＩ値との関係に基づいて補正することによって補正値を導出し、導出した補正値から施肥量導出計算式により追肥時の施肥量を導出するようにしている。
【００５７】
この構成によれば、被補正用圃場のそれぞれの追肥時の施肥量を導出するにあたって、遠方からの観測により得られた上空ＮＤＶＩ値がそのまま用いられて追肥時の施肥量が導出されるのではなく、地上ＮＤＶＩの指標値という精度の高い値に変換されてから施肥量導出計算式に入力され、導出されるため、高い精度で施肥量を導出できる。また、地上ＮＤＶＩ値の導出は、Ｎ個の全ての圃場について行われるのではなく、Ｎ個よりも少ないＭ個の圃場についてのみ行われるため、地上ＮＤＶＩ値を取得するための作業にかかる時間を短くすることができる。
【００５８】
＜変形例＞
次に上記実施形態の変形例について図１を援用し、上記実施形態と相違する部分を中心として説明する。本変形例に係る農作物関連値導出装置１は、Ｎ（ただし、Ｎ≧２）個の圃場を対象として追肥時施肥量を導出するものとし、Ｎ個の対象圃場のうち、Ｍ（ただし、Ｍ＜Ｎ）個の圃場が補正用圃場であるものとする。
【００５９】
上記実施形態では、第２指標値取得部１２が取得する値は、補正用圃場のそれぞれの地上ＮＤＶＩ値であった。一方、本変形例に係る第２指標値取得部１２は、補正用圃場のそれぞれのＳＰＡＤ（指標）の指標値の実際の値（以下「ＳＰＡＤ値」という）を取得する。詳述すると、記憶部１３には、補正用圃場のそれぞれについて、ＳＰＡＤ値が記録されたＳＰＡＤ値データが事前に記憶される。第２指標値取得部１２は、農作物関連値導出部１０からの要求に応じて、記憶部１３に記憶されたＳＰＡＤ値データに基づいて、補正用圃場毎のＳＰＡＤ値を取得する。
【００６０】
各補正用圃場についてのＳＰＡＤ値は、測定タイミングにおいて作業者が各補正用圃場に実際に赴き、葉緑素計（計測器）を用いて測定される。本変形例において、（指標としての）ＳＰＡＤは、特許請求の範囲の「第２指標」に相当する。また、ＳＰＡＤ値は、「領域の農作物を計測器で計測することにより得られた第２指標の指標値の実際の」および「葉緑素計を用いた計測により得られたＳＰＡＤの指標値の実際の値」に相当する。
【００６１】
また、本変形例に係る農作物関連値導出部１０は、以下の処理を実行する。すなわち、第１指標値取得部１１からＮ個の対象圃場についての上空ＮＤＶＩ値の応答を受け、第２指標値取得部１２からＭ個の補正用圃場についてのＳＰＡＤ値の応答を受けると、農作物関連値導出部１０は、まず、補正用計算式生成処理を実行する。
【００６２】
本変形例に係る補正用計算式生成処理において、本変形例に係る農作物関連値導出部１０は、補正用圃場のそれぞれの上空ＮＤＶＩ値とＳＰＡＤ値とに基づいて、上記実施形態と同様の方法で、上空ＮＤＶＩ値とＳＰＡＤ値との相関関係を分析する。本変形例に係る農作物関連値導出部１０は、分析結果に基づいて、上空ＮＤＶＩの指標値を入力値として入力し、補正用圃場を対象とした上空ＮＤＶＩ値とＳＰＡＤ値との相関関係に準じて、上空ＮＤＶＩ値を補正値へと変換し、この補正値を出力する補正用計算式を生成する。この補正値は、ＳＰＡＤ値と同じ単位で表される値である。
【００６３】
ここでＳＰＡＤの指標値は、葉緑素計を用いた計測に基づいて導出されるものであり、計測時に天候や気温、気象の状態によって影響を受けにくいという点において上空ＮＤＶＩの指標値よりも精度が高く、補正用圃場を対象とした上空ＮＤＶＩ値とＳＰＡＤ値との相関関係に準じて補正され、出力された補正値は、入力値よりも精度が高いということができる。
【００６４】
補正用計算式を生成した後、本変形例に係る農作物関連値導出部１０は、「Ｎ−Ｍ」個の被補正圃場のそれぞれについて、補正用計算式に上空ＮＤＶＩ値を入力し、その出力値を補正値として取得する。次いで、本変形例に係る農作物関連値導出部１０は、施肥量導出計算式を利用して、対象圃場の追肥時施肥量を導出する。本変形例では、施肥量導出計算式は、測定タイミングのＳＰＡＤ値を入力とし、農作物の状態を目標状態とするために必要かつ適切な追肥時施肥量を出力する計算式である。施肥量導出計算式は、上記実施形態で説明した方法と同様の方法で事前に生成される。
【００６５】
追肥時施肥量を導出する方法は上記実施形態と同様であり、本変形例に係る農作物関連値導出部１０は、補正用圃場については、施肥量導出計算式にＳＰＡＤ値を入力し、その出力を追肥時施肥量として導出し、また、被補正圃場については、施肥量導出計算式に補正値を入力し、その出力を追肥時施肥量として導出する。本変形例に係る農作物関連値導出装置１が奏する効果は、上記実施形態と同様である。なお、第２指標の指標値の指標値として、変形例のＳＰＡＤの指標値に代えて、他の指標値を用いてもよいことは勿論である。一例として、葉色板を用いた地上での計測により得られる葉色板値を第２指標の指標値の指標値として用いることができる。
【００６６】
以上、本発明の一実施形態（変形例を含む）を説明したが、上記実施形態は、本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその要旨、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００６７】
例えば上記実施形態では、農作物関連値導出装置１は、農作物関連値として稲（農作物）についての追肥時の施肥量（所定のタイミングにおける施肥量）を導出した。しかしながら、農作物は稲に限られるものではなく、また、農作物関連値は「追肥時の施肥量」に限られるものではない。まず、農作物関連値として追肥時の施肥量を導出する構成の場合、本発明は、栽培に際して追肥が行われる農作物に広く適用可能である。一例として小麦、大麦、大豆その他の穀物や、野菜、果物等に適用可能である。
【００６８】
農作物関連値は、所定のタイミングで測定された第２指標の指標値（上記実施形態では地上ＮＤＶＩの指標値、変形例ではＳＰＡＤの指標値）に基づいて導出可能な、農作物に関連する何らかの値であればよい。
【００６９】
例えば、農作物関連値に、所定のタイミングにおいて供給する堆肥、土壌改良資材、植物生育調整剤或いは薬剤（各種の農薬を含む）の量が含まれていてもよい。所定のタイミングにおいて供給する堆肥の量とは、所定のカテゴリについての農作物（農作物の収穫物も含む）の状態を目標状態とするために、所定のタイミングにおいて施すべき堆肥の量のことである。所定のタイミングにおいて供給する土壌改良資材、植物生育調整剤或いは薬剤の量についても同様である。以下、農作物が稲であり、農作物関連値が所定のタイミングにおける堆肥の量である場合の農作物関連値導出装置１の動作の一例について、上記実施形態を援用して説明する。
【００７０】
この例において、上記実施形態の施肥量導出計算式は、予め定められた測定タイミングにおける地上ＮＤＶＩの指標値を入力し、所定のタイミングにおける堆肥の量を出力する式（以下「堆肥量導出計算式」という）とされる。つまり堆肥量導出計算式は、測定タイミングにおける地上ＮＤＶＩの指標値が○○であったら、所定のタイミングにおける適切な堆肥の量は○○である、ということを計算する式である。発明者らの研究により、測定タイミングにおける地上ＮＤＶＩの指標値と所定のタイミングにおける堆肥の量とには相関関係があることが分かっており（土壌改良資材、植物生育調整剤および薬剤の量も同様）、この堆肥量導出計算式は、測定タイミングにおける地上ＮＤＶＩの指標値と、所定のタイミングにおける堆肥の量との相関関係に基づいて事前に導出される。
【００７１】
本例において農作物関連値導出部１０は、被補正圃場について上記実施形態と同様の方法で、上空ＮＤＶＩ値から補正値を導出する。そして、農作物関連値導出部１０は、導出した補正値を堆肥量導出計算式に入力することによって、所定のタイミングにおける堆肥の量を導出する。以上、農作物関連値が「稲についての所定のタイミングにおける堆肥の量」である場合の農作物関連値導出装置１の動作について説明した。ただし、対象となる農作物は稲でなくてもよく、また、導出される農作物関連値は、堆肥の量ではなく、土壌改良資材、植物生育調整剤或いは薬剤の量であってもよい。また堆肥量導出計算式（或いはこれに相当する計算式）に入力する補正値は、「地上ＮＤＶＩの指標値に対応する値」以外の値（一例として、「ＳＰＡＤの指標値に対応する値」）であってもよい。また、農作物関連値（上記例では所定のタイミングにおける堆肥の量）の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００７２】
また例えば農作物関連値を、収穫物の収量（単位面積当たりの予測収量を意味する）とすることができる。以下、農作物が稲であり、農作物関連値が稲の収量である場合の農作物関連値導出装置１の動作の一例について、上記実施形態を援用して説明する。この例において、第１指標値取得部１１は、対象圃場について、出穂前の所定のタイミング（例えば、出穂前の直前のタイミングや、出穂から１０日前のタイミング、出穂から２０日前のタイミング等。測定タイミングに相当）の上空ＮＤＶＩ値を取得する。また、第２指標値取得部１２は、補正用圃場について、出穂前の所定のタイミングの地上ＮＤＶＩ値を取得する。また本例において、上記実施形態の施肥量導出計算式は、出穂前の所定のタイミングにおける地上ＮＤＶＩの指標値を入力し、稲の収量を出力する式（以下「予測収量導出計算式」という）とされる。つまり予測収量導出計算式は、出穂前の所定のタイミングにおける地上ＮＤＶＩの指標値が○○であったら、稲の収量は○○と予測される、ということを計算する式である。発明者らの研究により、（出穂後ではなく）出穂前の所定のタイミングにおける地上ＮＤＶＩの指標値と稲の収量とには特に強い相関関係があることが分かっており、この予測収量導出計算式は、出穂前の所定のタイミングにおける地上ＮＤＶＩの指標値と稲の収量との相関関係に基づいて事前に導出される。
【００７３】
本例において農作物関連値導出部１０は、被補正圃場について上記実施形態と同様の方法で、上空ＮＤＶＩ値から補正値を導出する。そして、農作物関連値導出部１０は、導出した補正値を予測収量導出計算式に入力することによって、稲の収量を導出する。なお本例では、予測収量導出計算式は、出穂前の所定のタイミングの地上ＮＤＶＩの指標値を入力とするものであったが、これに代えて、出穂前の所定のタイミングのＳＰＡＤの指標値を入力とするようにしてもよい。発明者らの研究により、出穂前の所定のタイミングのＳＰＡＤの指標値と稲についての収穫物の収量とには強い相関関係が分かっている。また施肥量導出計算式は、対象となる稲の品種毎に生成されることが望ましい。この他、収穫物の収量と、施肥窒素や吸収窒素との関係式を予め作成し、これを基に翌年の施肥の量や時期を決定するようにしてもよい。この場合において、この関係式は、品種毎に関係式を作成されることが望ましい。
【００７４】
以上、農作物関連値が収穫物の収量である場合について、収穫物の収量が稲の収量である場合を例として説明したが、農作物関連値としての収穫物の収量は、第２指標の指標値との間での所定の関係性があり、第２指標の指標値との関係性に基づいて導出可能な収量であればよく、当然、稲以外の農作物についての収穫物の収量であってもよい。ただし上記実施形態でも説明した通り、農作物関連値の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００７５】
また、農作物関連値を、農作物の現在或いは将来の品質（外観品質のみならず、所定の栄養素の状態（例えばタンパク質の含有率）や、含有水分の状況などの成分的、内部的な品質も含む概念）とすることができる。農作物の現在或いは将来の品質は例えば、農作物が稲である場合において、現在或いは将来の籾水分や、現在或いは将来の玄米水分、現在或いは将来の玄米のタンパク質含有率である。また例えば農作物の現在或いは将来の品質は、農作物が稲である場合において、収穫時における低外観品質の収穫物の割合（或いは整粒の割合）である。農作物関連値としての「農作物の現在或いは将来の品質」は、第２指標の指標値との間で所定の関係性があり、第２指標の指標値との関係性に基づいて導出可能な品質であればよく、当然、稲以外の農作物についての品質であってもよい。ただし上記実施形態でも説明した通り、農作物関連値の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００７６】
なお農作物が稲であり、農作物関連値がタンパク質の状態（例えば、タンパク質の含有率）の場合は、測定タイミングは出穂前３０日から成熟期に至るまでの期間の所定のタイミングとされる。また第２指標は「ＳＰＡＤ」とされる。つまり第２指標値取得部１２は、補正用圃場について、ＳＰＡＤの指標値の実際の値を第２実際値として取得する。またＳＰＡＤの指標値を入力し、現在或いは将来（例えば収穫時）のタンパク質の状態を出力する計算式が事前に生成される。第２指標がＳＰＡＤとされる理由は、発明者の研究により、測定タイミングにおける「農作物の葉について測定された葉緑素に関する指標」（当然、ＳＰＡＤが含まれる）の値と、タンパク質の状態とには特に強い相関関係があることが分かっており、第２指標を「農作物の葉について測定された葉緑素に関する指標」とし、計算式として、「農作物の葉について測定された葉緑素に関する指標」の指標値を入力し、現在或いは将来（例えば収穫時）のタンパク質の状態を出力する式を用意することによって、より高い精度でタンパク質の状態を出力できるからである。計算式として、異なる複数の測定タイミングについての計算式が用意されてもよく、また、品種毎に異なる計算式が用意されてもよい。なお、タンパク質と施肥窒素や吸収窒素との関係式を予め作成し、これを基に翌年の施肥の量や時期を決定してもよい。この場合、品種毎に関係式を作成することが望ましい。このほか、施肥の量に代えて或いは施肥の量と併せて、堆肥や土壌改良資材の量等を決定するようにしてもよい。
【００７７】
ユーザは、導出されたタンパク質の状態に基づいて、例えば、タンパク質の高い圃場（領域）、中程度の圃場、低い圃場毎に分けて収穫・乾燥・調整し、出荷するといった取り組みを実施することができ、こういった取り組みにより食味の優れる米を提供することが可能となる。
【００７８】
次に農作物関連値が「農作物の現在或いは将来の品質」の場合における農作物関連値導出装置１の動作について、農作物が稲であり、農作物関連値が「収穫時（将来）に予想される籾水分の値」である場合を例にして、上記実施形態を援用して説明する。なお、以下の例において、ある圃場についての籾水分とは、その圃場の相当数（例えば３００〜５００粒）の籾の籾水分の平均値に相当する値を表しているものとする。
【００７９】
本例において、測定タイミングは、出穂期から成熟期に至るまでの期間の所定のタイミングとされる。本例に係る第１指標値取得部１１は、上記実施形態と同様、１０個の対象圃場のそれぞれの上空ＮＤＶＩ値を取得し、取得した上空ＮＤＶＩ値を各圃場の識別情報と対応付けて本例に係る農作物関連値導出部１０に出力する。また本例に係る第２指標値取得部１２は、３つの補正用圃場について、籾水分実際値を取得し、取得した籾水分実際値を各圃場の識別情報と対応付けて本例に係る農作物関連値導出部１０に出力する。なお籾水分実際値とは、「指標としての籾水分」（以下、便宜的に「籾水分指標」という）の指標値の実際の値を意味する。つまり、ある補正用圃場についての籾水分実際値とは、その補正用圃場についての籾水分の実際の値を意味する。各補正用圃場についての籾水分実際値は、圃場の籾が採取された上で、水分測定器を利用して導出される。各補正用圃場についての籾水分実際値を示す情報が記憶部１３に記憶されている。籾水分指標は特許請求の範囲の「第２指標」に相当し、籾水分実際値は「第２実際値」に相当する。また籾水分実際値は、「農作物を計測器で計測することにより得られた第２指標の指標値の実際の値」に相当する。
【００８０】
本例に係る農作物関連値導出部１０は、本例に係る第１指標値取得部１１から取得した１０個の対象圃場についての上空ＮＤＶＩ値、および、本例に係る第２指標値取得部１２から取得した３つの補正用圃場についての籾水分実際値に基づいて補正用計算式を生成する。補正用計算式は、圃場の上空ＮＤＶＩの指標値を入力値とし、補正値を出力する計算式である。補正用計算式の出力値（補正値）は、上空ＮＤＶＩ値を、補正用圃場についての上空ＮＤＶＩ値と籾水分実際値との相関関係に準じて、上空ＮＤＶＩの指標値から籾水分指標の指標値へと変換されるように補正した値である。本例に係る農作物関連値導出部１０は、上記実施形態と同様の方法で、補正用計算式を生成する。
【００８１】
補正用計算式を生成した後、本例に係る農作物関連値導出部１０は、補正用圃場ではない７個の被補正圃場のそれぞれについて、補正用計算式に上空ＮＤＶＩ値を入力し、その出力値を補正値として取得する。次いで本例に係る農作物関連値導出部１０は、３個の補正用圃場については、籾水分導出計算式（内容については後述）に籾水分実際値を入力し、「収穫時に予想される籾水分の値」（以下「籾水分予測値」という）を導出する。一方、本例に係る農作物関連値導出部１０は、７個の被補正圃場については、籾水分導出計算式に補正値を入力し、籾水分予測値を導出する。籾水分導出計算式とは、測定タイミングにおける籾水分指標の指標値を入力し、収穫時に予想される籾水分の値を出力する計算式である。つまり籾水分導出計算式は、測定タイミングにおける籾水分指標の指標値（籾水分の値）が○○であったら、収穫時の籾水分の値は○○と予測される、ということを計算する式である。発明者らの研究により、測定タイミングにおける籾水分指標の指標値（籾水分の値）と、収穫時の籾水分の値とには相関関係があることが分かっており、この籾水分導出計算式は、測定タイミングにおける籾水分指標の指標値と収穫時の籾水分の値との相関関係に基づいて事前に生成される。
【００８２】
本例に係る農作物関連値導出部１０により導出された籾水分予測値を示す情報がユーザに提供されることにより、ユーザは、各圃場の籾水分予測値を考慮して、各圃場の収穫時期を決定することができ、また、籾水分予測値が低い圃場から収穫を行っていく等の対応が可能となる。なお本例では、農作物関連値が「収穫時に予想される籾水分の値」であったが、これを「収穫時に予想される玄米水分の値」とし、籾水分予測値と同様の方法で「収穫時に予想される玄米水分の値」を求める構成でもよい。
【００８３】
また、農作物関連値を、「農作物が所定の状態となると予測される時期」とすることができる。「農作物が所定の状態となると予測される時期」は例えば、農作物が収穫可能な状態となる時期や、農作物が稲である場合において、出穂する時期である。農作物関連値としての「農作物が所定の状態となると予測される時期」は、第２指標の指標値との間で所定の関係性があり、第２指標の指標値との関係性に基づいて導出可能な時期であればよく、当然、稲以外の農作物についての時期であってもよい。ただし上記実施形態でも説明した通り、農作物関連値の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００８４】
本例における農作物関連値導出装置１の動作について、農作物が稲であり、農作物関連値が「収穫時期」（稲を収穫する適切な時期）である場合を例にして簡単に説明する。農作物関連値導出部１０は、被補正圃場について上記実施形態と同様の方法で、上空ＮＤＶＩ値から補正値を導出する。農作物関連値導出部１０は、補正値を所定の計算式に入力することによって、被補正圃場についての収穫時期を導出する。当該所定の計算式は、事前に定められた測定タイミングの地上ＮＤＶＩの指標値を入力し、収穫時期を出力する式である。地上ＮＤＶＩの指標値と収穫時期とには相関関係があることが分かっており、当該所定の計算式は、地上ＮＤＶＩの指標値と収穫時期との相関関係に基づいて事前に生成される。以上の処理により、被補正圃場についての収穫時期が導出される。
【００８５】
また、農作物関連値を、「農作物についての被害の状況」とすることができる。「農作物についての被害の状況」は例えば、倒伏の状況（程度）や、枯れ被害の状況（程度）、病害の状況(程度)、害虫による被害の状況（程度）等である。農作物関連値としての「農作物についての被害の状況」は、第２指標の指標値との間で所定の関係性があり、第２指標の指標値との関係性に基づいて導出可能なものであればよく、当然、稲以外の農作物についての被害の状況であってもよい。ただし上記実施形態でも説明した通り、農作物関連値の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００８６】
以下、農作物関連値が「農作物についての被害の状況」の場合における農作物関連値導出装置１の動作について、農作物が稲であり、農作物関連値が「収穫時に予想される倒伏の状況」である場合を例にして、上記実施形態を援用して説明する。なお本例において、ある１つの圃場についての倒伏の状況とは、以下を意味する。すなわち、倒伏の程度を示す複数のレベル（以下「倒伏レベル」という）が事前に設定されている。そして、ある１つの圃場についての倒伏の状況は、倒伏の程度に対応するレベルを意味する。
【００８７】
本例において、測定タイミングは、出穂期から成熟期に至るまでの期間の所定のタイミングとされる。本例に係る第１指標値取得部１１は、上記実施形態と同様、１０個の対象圃場のそれぞれの上空ＮＤＶＩ値を取得し、取得した上空ＮＤＶＩ値を各圃場の識別情報と対応付けて本例に係る農作物関連値導出部１０に出力する。また本例に係る第２指標値取得部１２は、３つの補正用圃場について、倒伏レベル値を取得し、取得した倒伏レベル値を各圃場の識別情報と対応付けて本例に係る農作物関連値導出部１０に出力する。なお倒伏レベル値とは、倒伏レベルの具体的な値であり、例えば、倒伏レベルとして倒伏の程度が低い順にレベル１〜６が設定されている場合におけるレベル３というものである。各補正用圃場についての倒伏レベル値は、担当者が現場で倒伏の状況を目視にて確認した上で、担当者により設定される。以下、指標としての倒伏レベルを便宜的に「倒伏指標」という。各補正用圃場についての倒伏レベル値を示す情報が記憶部１３に記憶されている。倒伏指標は特許請求の範囲の「第２指標」（農作物の被害の状況に係る指標）に相当し、倒伏レベル値は「第２実際値」に相当する。また倒伏レベル値は、「農作物を現場で観察することにより得られた第２指標の指標値の実際の値」に相当する。
【００８８】
本例に係る農作物関連値導出部１０は、本例に係る第１指標値取得部１１から取得した１０個の対象圃場についての上空ＮＤＶＩ値、および、本例に係る第２指標値取得部１２から取得した３つの補正用圃場についての倒伏レベル値に基づいて補正用計算式を生成する。補正用計算式は、圃場の上空ＮＤＶＩの指標値を入力値とし、補正値を出力する計算式である。補正用計算式の出力値（補正値）は、上空ＮＤＶＩ値を、補正用圃場についての上空ＮＤＶＩ値と倒伏レベル値との相関関係に準じて、上空ＮＤＶＩの指標値から倒伏指標の指標値へと変換されるように補正した値である。本例に係る農作物関連値導出部１０は、上記実施形態と同様の方法で、補正用計算式を生成する。
【００８９】
補正用計算式を生成した後、本例に係る農作物関連値導出部１０は、補正用圃場ではない７個の被補正圃場のそれぞれについて、補正用計算式に上空ＮＤＶＩ値を入力し、その出力値を補正値として取得する。次いで本例に係る農作物関連値導出部１０は、３個の補正用圃場については、倒伏レベル導出計算式（内容については後述）に倒伏レベル値を入力し、「収穫時に予想される倒伏レベルの値」（以下「倒伏レベル予測値」という）を導出する。一方、本例に係る農作物関連値導出部１０は、７個の被補正圃場については、倒伏レベル導出計算式に補正値を入力し、倒伏レベル予測値を導出する。倒伏レベル導出計算式とは、測定タイミングにおける倒伏指標の指標値を入力し、収穫時に予想される倒伏レベルの値（倒伏の状況）を出力する計算式である。つまり倒伏レベル導出計算式は、測定タイミングにおける倒伏指標の指標値（倒伏レベルの値）が○○であったら、収穫時の倒伏のレベルは○○と予測される、ということを計算する式である。発明者らの研究により、測定タイミングにおける倒伏指標の指標値（倒伏レベルの値）と、収穫時の倒伏レベルの値とには相関関係があることが分かっており、この倒伏レベル導出計算式は、測定タイミングにおける倒伏指標の指標値と収穫時の倒伏レベルの値との相関関係に基づいて事前に生成される。ただし、倒伏レベル予測値の導出にあたって、他の要素値を利用してもよいことは勿論である。
【００９０】
本例に係る農作物関連値導出部１０により導出された倒伏レベル予測値を示す情報がユーザに提供されることにより、ユーザは、各圃場の倒伏レベル予測値を考慮して、各圃場の収穫時期を決定することができ、また、倒伏程度の低い圃場から収穫を行っていく等の対応が可能となる。以上、農作物関連値が「収穫時に予想される倒伏の状況」の場合における農作物関連値導出装置１の動作を説明したが、枯れ被害の状況（程度）、病害の状況(程度)、害虫による被害の状況（程度）等についても同様の処理を行って、導出することが可能である。
【００９１】
また、農作物関連値を「農作物に影響を与える環境の状態」とすることができる。「農作物に影響を与える環境の状態」は例えば、圃場における病害虫の発生の程度や、圃場における雑草の発生の程度等である。農作物関連値としての「農作物に影響を与える環境の状態」は、第２指標の指標値との間で所定の関係性があり、第２指標の指標値との関係性に基づいて導出可能なものであればよく、当然、稲以外の農作物に関する「農作物に影響を与える環境の状態」であってもよい。ただし上記実施形態でも説明した通り、農作物関連値の導出にあたって、第２指標の指標値に加え、第２指標の指標値以外の要素値を利用してもよい。
【００９２】
本例における農作物関連値導出装置１の動作について、農作物が稲であり、農作物関連値が「圃場における病害虫の発生の程度」である場合を例にして簡単に説明する。なお病害虫の発生の程度は、例えば、予め範囲が定められたレベルによって表される。農作物関連値導出部１０は、被補正圃場について上記実施形態と同様の方法で、上空ＮＤＶＩ値から補正値を導出する。農作物関連値導出部１０は、補正値を所定の計算式に入力することによって、被補正圃場についての病害虫の発生の程度を導出する。当該所定の計算式は、所定のタイミングの地上ＮＤＶＩの指標値を入力し、そのタイミングにおける「圃場における病害虫の発生の程度」を出力する式である。地上ＮＤＶＩの指標値と病害虫の発生の程度とには相関関係があることが分かっており、当該所定の計算式は、地上ＮＤＶＩの指標値と病害虫の発生の程度との相関関係に基づいて事前に生成される。以上の処理により、被補正圃場についての病害虫の発生の程度が導出される。
【００９３】
また上記実施形態では、農作物栽培領域は、１０個の圃場によって構成されており、農作物関連値導出部１０は、圃場毎に、追肥時の施肥量（農作物関連値）を導出した。この点について、圃場毎に農作物関連値を導出する構成以外の構成でもよいことは勿論である。例えば、１つの圃場をＮ個の領域に区分けし、このうちのＭ個の領域を上記実施形態における補正用圃場に相当する領域とし、他の（Ｎ−Ｍ）個の領域を上記実施形態における被補正圃場に相当する領域としてもよい。
【００９４】
また、対象圃場のうち、補正用圃場を選択するときのルールは、上記実施形態で例示したものに限らない。例えば、生育量以外の要素に基づいて補正用圃場を選択するようにしてもよい。ただし補正用圃場の個数が少ないほど、作業に要する時間が短くなるものの、補正用計算式を生成するときのサンプル数が少なくなるため補正用計算式の精度が低くなり、逆に、補正用圃場の個数が多いほど補正用計算式の精度が高くなる一方、作業に要する時間が長くなる。補正用圃場の個数は、このことを考慮して定められる必要がある。
【００９５】
また、上空ＮＤＶＩ値を、補正用圃場の上空ＮＤＶＩ値と地上ＮＤＶＩとの関係によって補正する方法は、上記実施形態で例示した方法に限られない。すなわち、補正用圃場の上空ＮＤＶＩ値と地上ＮＤＶＩとの関係に基づいて上空ＮＤＶＩ値が補正されることにより、上空ＮＤＶＩ値の精度が向上する方法であれば、どのような方法であってもよい。一例として、農作物関連値導出部１０が、補正用圃場の上空ＮＤＶＩ値と地上ＮＤＶＩ値との相関関係を分析することによって、上空ＮＤＶＩの指標値を地上ＮＤＶＩの指標値に近づけるための加算値または減算値を固定値として導出し、被補正圃場の上空ＮＤＶＩ値に対して加算値を加算しまたは減算値を減算して補正値を求めるようにしてもよい。
【００９６】
また、記憶部１３は、農作物関連値導出装置１とは異なる外部記憶装置が有する構成でもよい。
【００９７】
また、農作物関連値導出装置１を、インターネット上のサーバとして構成し、インターネットを介してクライアントから追肥時施肥量を示す情報の提供の要請があったときに、上述した手法で追肥時施肥量を導出し、クライアントに提供する構成としてもよい。
【００９８】
また、農作物関連値導出装置１は、単体のコンピュータである必要はなく、複数のコンピュータにより構成されていてもよい。例えば、インターネットを介して接続された端末とクラウドサーバとが協働して農作物関連値導出装置１として機能し、端末が適宜クラウドサーバと協働して処理を実行する構成でもよい。
【００９９】
また、上記実施形態では、農作物関連値導出装置１は、３個（Ｍ個）の補正用圃場については、施肥量導出計算式に３個の補正用圃場についての地上ＮＤＶＩ値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出決定した。この点に関し、３個の補正用圃場についても、上空ＮＤＶＩ値から補正用計算式を用いて補正値を取得し、補正値を施肥量導出計算式に入力し、当該計算式の出力値を最終的な追肥時施肥量として導出する構成でもよい。
【符号の説明】
【０１００】
１  農作物関連値導出装置
１０  農作物関連値導出部
１１  第１指標値取得部
１２  第２指標値取得部
【図１】
0006960698.tif 000002
【図２】
0006960698.tif 000003
【図３】
0006960698.tif 000004
【図４】
0006960698.tif 000005
【図５】
0006960698.tif 000006
【図６】
0006960698.tif 000007