JP7044285B2 - 農作物関連値導出装置および農作物関連値導出方法 - Google Patents

Description

本発明は、農作物に関連する農作物関連値を導出する農作物関連値導出装置および農作物関連値導出方法に関する。

従来、圃場で栽培される農作物について、圃場に対する遠隔からの観測により得られた所定の指標の指標値と、追肥時の適切な施肥量との間に相関関係があることが分かっており、圃場への追肥にあたって、圃場に対する遠隔からの観測により所定の指標の指標値の実際の値（実際値）を取得し、実際値に基づいて追肥時の施肥量を求めることが行われている。一例として、ドローンを用いた上空からのセンシングによりＮＤＶＩに係る実際値を取得し、これを用いて追肥時の施肥量を求めることが行われている。ここで特許文献１には、圃場に対する遠隔からの観測により得られる実際値に関して、圃場全体が撮影可能な程度に圃場から遠隔に設けられたカメラにより圃場全体を撮影し、撮影結果に基づいて圃場の作物情報（実際値）を求める技術が記載されている。

特開２００８－１３６４１１号公報

追肥時の施肥量は、追肥後の農作物の成長や、最終的な収穫物の質や量等に強い影響を与える。従って、遠隔からの観測により得られる実際値に基づいて施肥量を導出する装置については、遠隔からの観測により実際値を得るという特徴を踏まえた処理を行って、高い精度で施肥量を導出することが求められる。

以上の課題は、農作物を栽培する農作物栽培領域を構成する領域を遠隔から観測して実際値を取得し、この実際値に基づいて農作物に関連する農作物関連値を導出する技術において、広く発生し得る。農作物関連値は、収穫物の収量、農作物の現在或いは将来の品質、農作物が所定の状態となると予測される時期、農作物についての被害の状況、および、農作物に影響を与える環境の状態等である。

本発明は、このような問題を解決するために成されたものであり、遠隔からの観測により実際値を得るという特徴を踏まえた処理を行って、高い精度で農作物関連値を導出できるようにすることを目的とする。

上記した課題を解決するために、本発明では、農作物栽培領域を構成する領域のうち、Ｎ（ただし、Ｎ≧２）個の領域のそれぞれについて、Ｎ個の領域に対して行われた、高度を変更可能な飛行物からの所定高度よりも高い第１高度からの観測により得られた第１指標の指標値の実際の値である第１実際値を取得し、Ｎ個の領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の領域のそれぞれについて、Ｍ個の領域に対して個別に行われた所定高度からの観測により得られた第１指標の指標値の実際の値である第２実際値を取得する。そしてＮ個の領域のうち、第２実際値が取得されていない一の領域について、農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する際に、当該一の領域について取得した第１実際値を、Ｍ個の領域についての第１実際値と第２実際値との関係に基づいて所定高度で観測が行われたとした場合の値となるように補正して補正値を導出し、第１指標の指標値を入力とし農作物関連値を出力とする計算式であって、所定高度で観測が行われた場合の第１指標の指標値を入力したときに適切な値の農作物関連値を出力するように調整された計算式である第１計算式に、当該一の領域について導出した補正値を入力することによって、農作物関連値を導出するようにしている。

農作物栽培領域を構成する領域に対する遠隔からの観測により実際値を得る場合において、観測が行われたときの状況に応じて実際値の値が変化する場合があり、観測が行われたときの状況による実際値の変化を考慮せずに農作物関連値を導出すると、導出された農作物関連値の精度が低下することがある。一例として同じ圃場を同じスペックの観測機で同じ位置から観測した場合であっても、日射量が変化すると実際値の値が変化し、このような実際値の変化を考慮せずに農作物関連値を導出すると、農作物関連値の精度が低下してしまう。また、領域に対する遠隔からの観測により実際値を得る場合において、観測が行われたときの状況に応じて実際値自体の精度の高低が変化する場合があり、観測が行われたときの状況による実際値の精度の高低の変化を考慮せずに農作物関連値を導出すると、導出された農作物関連値の精度が低下することがある。一例として、同じ圃場を同じスペックの観測機器で同じ自然環境で観測した場合であっても、観測位置の高度が高いほど、実際値の精度が低くなり、このような実際値の精度の高低の変化を考慮せずに農作物関連値を導出すると、農作物関連値の精度が低下してしまう。

以上を踏まえ、上記のように構成した本発明によれば、農作物関連値を導出するにあたって、領域に対する遠隔からの観測により得られた第１実際値がそのまま第１計算式に入力されるのではなく、所定の状況で観測が行われた場合の値となるように補正された上で入力される。このため、観測が行われたときの状況に応じて実際値の値が変化する場合に、観測時の実際の状況にかかわらず、所定の状況で観測が行われた場合の値となるように実際値を補正し、補正後の実際値を第１計算式に入力することができ、これにより観測時の状況に応じた実際値の変化に起因する精度の低下を抑制した状態で農作物関連値を導出することができる。

また、観測が行われたときの状況に応じて実際値の精度の高低が変化する場合に、得られる実際値の精度が低い状況で行われた観測に基づいて導出された実際値を、得られる実際値の精度が高い状況で観測が行われた場合の値となるように補正し、補正後の実際値を第１計算式に入力することができ、これにより観測時の状況に応じた実際値の精度の高低の変化に起因する精度の低下を抑制した状態で農作物関連値を導出することができる。すなわち、本発明によれば、遠隔からの観測により実際値が得るという特徴を踏まえた処理を行って、高い精度で農作物関連値を導出できる。

本発明の第１実施形態に係る農作物関連値導出装置の機能構成例を示すブロック図である。 縦軸が補正ＮＤＶＩ値、横軸が低レベル時ＮＤＶＩ値の２次元空間で低レベル時ＮＤＶＩ値変換式を示す図である。 測定タイミングのＮＤＶＩと追肥時の施肥量との関係を示す図である。 本発明の第１実施形態に係る農作物関連値導出装置の動作例を示すフローチャートである。 本発明の第２実施形態に係る農作物関連値導出装置の機能構成例を示すブロック図である。 対象圃場を識別情報と共に示す図である。 補正用圃場の低空ＮＤＶＩ値と上空ＮＤＶＩ値との関係を示す図である。 補正用計算式の一例によって求められる補正ＮＤＶＩ値を示す図である。 本発明の第２実施形態に係る農作物関連値導出装置の動作例を示すフローチャートである。

＜第１実施形態＞
以下、本発明の第１実施形態を図面に基づいて説明する。図１は、本実施形態に係る農作物関連値導出装置１の構成例を示すブロック図である。農作物関連値導出装置１は、稲作に関して、追肥（追肥についての詳細は後述）時に施す肥料の量（以下「追肥時施肥量」という）を導出し、追肥時施肥量を示す情報をユーザに提供する装置である。追肥時施肥量は、特許請求の範囲の「農作物関連値」に相当する。特に本実施形態に係る農作物関連値導出装置１は、複数の圃場（領域）を対象として、複数の圃場のそれぞれについて追肥時施肥量を導出する。ユーザは、追肥にあたって農作物関連値導出装置１により提供された追肥時施肥量を示す情報を参考にして、実際に施す肥料の量を適切に決定することができる。なお、本実施形態で「施肥量」という場合、特に説明がない限り、単位面積あたりの肥料の量を指す。

以下の説明では、対象とする稲および稲作について、稲の品種や、稲作が行われる環境（地域（寒冷地域や、温暖地域、乾燥地域等）、標高、圃場の規模）、稲作の方法、肥料の種類、その他の稲の生育に影響を与える外部的な要素は共通しているものとする。従って、農作物関連値導出装置１により導出された追肥時施肥量は、特定の環境で特定の方法により栽培される特定の品種の稲について、特定の種類の肥料を施す場合に適切な施肥量ということである。

本実施形態において、「追肥」とは以下を意味する。すなわち一般に、ある１つの期において、稲は、田植えによって苗代から圃場に移植された後、分げつ期→幼穂形成期→出穂期→開花・受粉期→穂揃期→登熟期→成熟期という段階を経て成長する。幼穂形成期前後から出穂期前後に至るまでの期間（以下「追肥対象期間」という）では、通常、穂に着生する籾の個数の増大や、籾に詰まるデンプンの量の増大等を目的として、所定のタイミングで穂肥と呼ばれる肥料の追加が行われる。そして本実施形態では、「追肥」とは、追肥対象期間における所定のタイミングで行われる肥料の追加を意味する。

図１で示すように、農作物関連値導出装置１には、液晶ディスプレイや、有機ＥＬディスプレイ等の表示装置２と、マウスやキーボード等の入力装置３とが接続されている。また図１で示すように、農作物関連値導出装置１は、機能構成として、第１実際値取得部１０、補正部１１および農作物関連値導出部１２を備えている。上記各機能ブロック１０～１２は、ハードウェア、ＤＳＰ（Digital Signal Processor）、ソフトウェアの何れによっても構成することが可能である。例えばソフトウェアによって構成する場合、上記各機能ブロック１０～１２は、実際にはコンピュータのＣＰＵ、ＲＡＭ、ＲＯＭなどを備えて構成され、ＲＡＭやＲＯＭ、ハードディスクまたは半導体メモリ等の記録媒体に記憶されたプログラムが動作することによって実現される。また、農作物関連値導出装置１は、記憶手段として記憶部１３を備えている。記憶部１３に記憶されるデータについては後述する。

本実施形態では、農作物関連値導出装置１は、ユーザから入力装置３に対して追肥時施肥量を示す情報の提供の指示（以下「情報提供指示」という）があったことをトリガとして、圃場毎の追肥時施肥量を導出し、当該情報を提供する。以下、ユーザから入力装置３に対して情報提供指示があったときの農作物関連値導出装置１の動作について詳述する。以下の説明では、説明の便宜のため、農作物関連値導出装置１は、１０個の圃場を対象として追肥時施肥量を導出するものとする。また追肥時施肥量を導出する対象の１０個の圃場を「対象圃場」という。なお追肥時施肥量の導出に関する処理を開始するトリガは本実施形態で例示するものに限られず、例えば、予め定められた日時が到来したことをトリガとしてもよい。

第１実際値取得部１０は、入力装置３に対して情報提供指示があったことを検出すると、以下の処理を実行する。すなわち、まず第１実際値取得部１０は、記憶部１３に記憶された上空データを取得する。上空データとは、追肥を行うタイミングの前の所定のタイミング（以下「測定タイミング」という）で、ドローンに搭載されたマルチスペクトルセンサにより、１０個の圃場の稲の群落が上空（遠隔）から一括して撮影されることによって生成された画像データである。測定タイミングは、追肥時施肥量の導出に利用する上空データを得るための撮影を行うタイミングとして予め定められたタイミングであり、例えば、追肥を行うタイミングの直前や、出穂の３週間前、出穂の４週間前等とされる。なお、本実施形態では、上空データの生成に際し、対象となる全ての圃場が一括して撮影されたこととしているが、圃場の群が上空から複数回にわたって撮影され、その複数回の撮影結果に基づいて上空データが生成されてもよい。この場合に、上空データは、単一のデータである必要はなく、複数のデータの集合であってもよい。

ただし本実施形態では、ドローンに搭載されたマルチスペクトルセンサによる撮影は、撮影タイミング（特定の時期）において特定のスペックのマルチスペクトルセンサを利用して、特定の時間帯で特定の位置（高度を含む３次元的な位置）から行われる。例えば、ドローンが無線操縦機による遠隔操縦によって飛行するものである場合、ドローンの３次元的な位置を示す情報が常時、操縦者に対して提供され、操縦者は、提供された情報に基づいて、ドローンの３次元的な位置を把握しつつ、ドローンを特定の位置に至らせ、その後、マルチスペクトルセンサによる撮影を実行させる。また、ドローンが自律飛行可能な場合は、特定の位置が目標飛行位置として設定され、ドローンがこの目標飛行位置に至った後、撮影が実行される。

第１実際値取得部１０は、取得した上空データを分析し、周知の方法で対象圃場毎のＮＤＶＩの指標値の実際の値（実際値）を導出し、取得する。なお、上空データにおける個々の対象圃場の領域を示す情報等の、上空データに基づいて対象圃場毎のＮＤＶＩに係る実際値を導出するために必要な情報は事前に登録されている。以下、上空データに基づいて第１実際値取得部１０により取得された対象圃場毎のＮＤＶＩに係る実際値を特に「上空ＮＤＶＩ値」という。ＮＤＶＩ（指標としてのＮＤＶＩ）は、特許請求の範囲の「第１指標」に相当し、特に「ドローンを用いた上空からのセンシングにより得られるＮＤＶＩ」に相当する。また上空ＮＤＶＩ値は、特許請求の範囲の「第１実際値」に相当し、

ここで本実施形態では、上空ＮＤＶＩ値の導出にあたって使用される上空データは、ドローンに搭載されたマルチスペクトルセンサの撮影結果に基づいて生成されるものである。このように、上空データの生成にあたってドローンを使用することにより、簡易な手段で複数の圃場にわたる広い領域を撮影することが可能であり、作業にかかる時間を短くし、また作業にかかる労力を低減することができる。なお、観測（センシング）に用いる装置は、マルチスペクトルセンサに限られず、一例として、赤外線サーモグラフィであってもよい。また、指標（第１指標）は、ＮＤＶＩに限られない。また、観測は、必ずしもドローンを用いて行われる必要はなく、ドローン以外の飛行物に搭載された装置や、塔の上部に設置された装置、専用の設置台に設置された装置等を用いて行われたりしてもよい。また、上空ＮＤＶＩ値（またはこれに相当する第１実際値）の導出に、衛星写真を用いるようにしてもよい。

上述の通りドローンを利用して圃場を撮影し、その撮影結果に基づいて上空ＮＤＶＩ値を導出する場合、作業にかかる時間の短縮、労力の低減を図ることができる。その一方で、撮影タイミングにおいて同じ圃場を、同じスペックのマルチスペクトルセンサにより、同じ時間帯で同じ位置（上述したように高度を含む３次元的な位置）から撮影し、その撮影結果に基づいて対象圃場毎の上空ＮＤＶＩ値を導出したとしても、撮影が行われたときの日射量（観測が行われたときの状況）によって、導出される上空ＮＤＶＩ値の値が変化する。これは、上空ＮＤＶＩ値は、基本的には太陽を光源とする反射光の観測により値が求められるため、撮影時の日射量が少ないときには値が小さくなり、逆に日射量が多いときには値が大きくなるからである。つまり日射量は、対象となる領域（圃場）、観測機器（マルチスペクトルセンサ）のスペック、観測の時期（測定タイミング）および時間帯、および、観測の位置が同じ場合に、観測される第１指標（ＮＤＶＩ）の指標値に変化を与える自然環境の要素である。

なお本実施形態では第１実際値取得部１０が上空ＮＤＶＩ値を導出し取得する構成であった。この点に関し、外部の装置により事前に対象圃場毎の上空ＮＤＶＩ値が導出され、第１実際値取得部１０がそれを取得する構成としてもよい。ただし、この場合においても、上空ＮＤＶＩ値の生成のための対象圃場の撮影は、特定のスペックの撮影装置により、特定の時間帯に特定の位置から行われる。

さて第１実際値取得部１０は、対象圃場のそれぞれの上空ＮＤＶＩ値を取得した後、取得した上空ＮＤＶＩ値を各圃場の識別情報と対応付けて補正部１１に出力する。なお、１０個の対象圃場のそれぞれには識別情報が事前に付与されており、個々の対象圃場は識別情報によって管理される。

補正部１１は、１０個の対象圃場のそれぞれについて、識別情報と共に上空ＮＤＶＩ値を入力すると、以下の処理を実行する。すなわち補正部１１は、記憶部１３にアクセスし、記憶部１３に記憶された撮影時日射量情報を取得する。撮影時日射量情報は、上空ＮＤＶＩ値の導出の元となった上空データについての撮影が行われたときの全天日射量（以下、単に「日射量」という）を示す情報である。撮影が行われたときの日射量は、例えば、ドローンに設けられた日射センサにより計測される。ただし、日射量の計測方法は、どのような方法でもよい。日射量は、撮影が行われたときの天候の外、撮影が行われたときの時期、撮影が行われたときの時間帯等によっても変化し得る。

次いで、補正部１１は、撮影時日射量情報が示す日射量（以下「撮影時日射量」という）の日射量レベル（特許請求の範囲の「日光の状況」に相当）を判定する。本実施形態では、日射量レベルとして、低レベル、中レベルおよび高レベル（ただし、低レベル＜中レベル＜高レベル）の３つのレベルが設定されており、低レベルと中レベルとの境の閾値、および、中レベルと高レベルとの境の閾値とが事前に定められている。一例として、「日射量≦２．０ＭＪ／ｍ^２／ｈ」のときは低レベル、「２．０ＭＪ／ｍ^２／ｈ≦日射量＜３．０ＭＪ／ｍ^２／ｈ」のときは中レベル、「３．０ＭＪ／ｍ^２／ｈ≦日射量」のときは高レベルと定められている場合において、撮影日射量が「１．５ＭＪ／ｍ^２／ｈ」の場合、補正部１１は、撮影時日射量の日射量レベルを低レベルと判定する。なお、本実施形態では、「中レベル」が特許請求の範囲の「所定日射量レベル」に相当する。

日射量レベルの判定後、補正部１１は、日射量レベルによって異なる処理を実行する。以下、補正部１１の処理を日射量レベル毎に説明する。以下では、日射量が低レベル、中レベルおよび高レベルのときの撮影結果に基づいて導出された上空ＮＤＶＩ値をそれぞれ、「低レベル時ＮＤＶＩ値」、「中レベル時ＮＤＶＩ値」および「高レベル時ＮＤＶＩ値」という。また各日射量レベルは、特許請求の範囲の「なり得る日光の状況」に相当する。

＜低レベルの場合＞
撮影時日射量の日射量レベルが低レベルの場合、補正部１１は、対象圃場のそれぞれについて、低レベル時ＮＤＶＩ値変換式（特許請求の範囲の「第２計算式」に相当）を用いて、低レベル時ＮＤＶＩ値から補正ＮＤＶＩ値（特許請求の範囲の「補正値」に相当）を導出する。詳述すると、低レベル時ＮＤＶＩ値変換式は、低レベル時ＮＤＶＩ値を、日射量が中レベルのときに撮影が行われていたとしたときに導出される値（厳密には、導出される値にできるだけ近づけた値）に変換する式である。

従って、あるタイミングのある対象圃場について、仮に日射量の日射量レベルが低レベルであった場合に導出されるＮＤＶＩがα値であり、仮に日射量レベルが中レベルであった場合に導出されるＮＤＶＩがβ値であったとする。この場合において、α値を低レベル時ＮＤＶＩ値変換式に入力すると、当該計算式からβ値或いはβ値に近い値が出力される。低レベル時ＮＤＶＩ値変換式により変換された値（つまり、低レベル時変換式の出力）が「補正ＮＤＶＩ値」である。なお、本実施形態（後述する第２実施形態も同様）において、計算式は、単純な数式だけでなく、プログラム上で定義された関数（モデル）や、プログラムから呼び出し可能な関数（モデル）を概念的に表している。

低レベル時ＮＤＶＩ値変換式は、以下の方法で事前に導出される。すなわち、随時行われる研究、実験により、「生育の状態が同じ（或いは近似する）圃場に対して、日射量が異なる状況で撮影が行われて導出された上空ＮＤＶＩ値の群」というデータが相当数、蓄積される。生育の状態が同じ圃場は、例えば、生育の状態を示す草丈、茎数、葉色等が実地で観測され、観測結果に基づいて判定される。また例えば、生育の状態が同じ圃場は、日射量に影響を受けない（或いは影響が限定的な）指標値が共通する圃場とされる。日射量に影響を受けない（或いは影響が限定的な）指標値は、例えば、光源を有する携帯型ＮＤＶＩ測定機により測定されたＮＤＶＩの指標値である。そして、対象となる圃場の生育の状態が共通し、かつ、ＮＤＶＩの指標値の導出に影響を与える環境のうち、日射量のみが変動し得るという状況における「日射量が低レベルのときに行われた撮影に基づくＮＤＶＩの指標値」と、「日射量が中レベルのときに行われた撮影に基づくＮＤＶＩの指標値」との相関関係が分析される。そして、分析された相関関係に基づいて、日射量レベルが低レベルであったときの撮影結果に基づいて導出されるＮＤＶＩの指標値を、中レベルで撮影が行われた場合の値（＝補正ＮＤＶＩ値）に変換する低レベル時ＮＤＶＩ値変換式が導出される。

図２は、横軸が低レベル時ＮＤＶＩ値であり、縦軸が補正ＮＤＶＩ値である二次元空間に、低レベル時ＮＤＶＩ値変換式を示す線分（低レベル時ＮＤＶＩ値と補正ＮＤＶＩ値との対応関係を示す点の集合よりなる線分）を描画した図である。図２の例では、低レベル時ＮＤＶＩ値変換式は「補正ＮＤＶＩ値＝低レベル時ＮＤＶＩ値+０．１」であり、例えば、低レベル時ＮＤＶＩ値が「０．４」の場合、補正ＮＤＶＩ値として「０．５」が導出される。なお、図２で示す低レベル時ＮＤＶＩ値変換式は、単純化した一例であり、二次以上の関数、その他の形式の関数として表される式であってもよく、また、関数ではなく、低レベル時ＮＤＶＩ値がとり得る値毎に、補正ＮＤＶＩ値を対応付けた表によって、低レベル時ＮＤＶＩ値と補正ＮＤＶＩ値との対応関係を定義付けたものであってもよい。

撮影時日射量の日射量レベルが低レベルの場合、補正部１１は、対象圃場のそれぞれについて、上述した低レベル時ＮＤＶＩ値変換式に上空ＮＤＶＩ値（低レベル時ＮＤＶＩ値）を入力することによって、補正ＮＤＶＩ値を導出する。

＜中レベルの場合＞
撮影時日射量の日射量レベルが中レベルの場合、補正部１１は、対象圃場のそれぞれについて、上空ＮＤＶＩ値（中レベル時ＮＤＶＩ値）をそのまま、補正ＮＤＶＩ値とする。なお、値の補正がなされないものについて「補正ＮＤＶＩ値」とするのは説明の便宜上のことである。

＜高レベルの場合＞
撮影時日射量の日射量レベルが高レベルの場合、補正部１１は、対象圃場のそれぞれについて、高レベル時ＮＤＶＩ値変換式（特許請求の範囲の「第２計算式」に相当）を用いて、上空ＮＤＶＩ値に基づいて補正ＮＤＶＩ値を導出する。高レベル時ＮＤＶＩ値変換式は、高レベル時ＮＤＶＩ値を、日射量が中レベルのときに撮影が行われていたとした場合に導出される値（或いは導出される値に近い値）に変換する式である。高レベル時ＮＤＶＩ値変換式は、低レベル時ＮＤＶＩ値変換式と同様の方法で事前に導出される。

以上のように、補正部１１は、撮影時日射量の日射量レベルに応じた方法で対象圃場毎の補正ＮＤＶＩ値を導出する。補正部１１は、導出した補正ＮＤＶＩ値を各圃場の識別情報と対応付けて農作物関連値導出部１２に出力する。

農作物関連値導出部１２は、補正部１１から、対象圃場毎に、補正ＮＤＶＩ値と識別情報との組み合わせを入力する。次いで、農作物関連値導出部１２は、対象圃場のそれぞれについて、施肥量導出計算式（特許請求の範囲の「第１計算式」に相当）に補正ＮＤＶＩ値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出する。

施肥量導出計算式は、ＮＤＶＩの指標値を入力し、追肥時施肥量を出力する計算式である。施肥量導出計算式に一のＮＤＶＩの指標値を入力したときに出力される施肥量は、追肥時のＮＤＶＩの指標値が入力された値であった場合に、所定のカテゴリについての農作物（農作物の収穫物も含む）の状態を目標状態とするために必要かつ適切な追肥時施肥量である。特に本実施形態に係る施肥量導出計算式は、日射量レベルが中レベルのときに行われた撮影結果に基づいて導出されるＮＤＶＩの指標値を入力し（所定の状況で観測が行われた場合の第１指標の指標値を入力し）、追肥時施肥量を出力する計算式である。

所定のカテゴリについての農作物の状態は例えば、農作物の外観的あるいは食味的な品質や、収穫物の収量、籾数、追肥後の生育過程の指標等である。追肥後の生育過程の指標は、例えば、追肥した後の穂揃期の葉緑素計値（ＳＰＡＤ値）である。なお、穂揃期の葉緑素計値（ＳＰＡＤ値）が一定値以上の場合は、品質が良いということができる。

例えば施肥量導出計算式の出力は、収穫物の収量（単位当たりの収量）を○○ｋｇとするために必要かつ適切な追肥時施肥量である。この場合、収穫物の収量が「所定のカテゴリについての農作物の状態」に相当し、収量が○○ｋｇであるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって、収穫物の収量が○○ｋｇとなることが期待される。

また例えば、施肥量導出計算式の出力は、収穫物全体に対する低外観品質の収穫物の割合（以下「低品質割合」という）を○○％とするために必要かつ適切な追肥時施肥量である。なお、低外観品質の収穫物は、例えば、白未熟粒（腹白粒や、背白粒、基部未熟粒、心白粒、乳白粒等）や、その他の未熟粒、被害粒、死米、着色粒等とである。この場合、低品質割合が「所定のカテゴリについての農作物の状態」に相当し、低品質割合が○○％であるという状態が「目標状態」に相当する。そしてこの場合、施肥量導出計算式の出力値分の肥料を追肥時に施すことによって、低品質割合が○○％となることが期待される。なお、低品質割合に代えて、収穫物の「整粒」の割合としてもよいことは勿論である。

図３は、測定タイミングにおけるＮＤＶＩの指標値を縦軸とし、追肥時の施肥量を横軸とする２次元空間に描画された１次式によって、上述した低品質割合（所定のカテゴリについての農作物の状態）を５％（目標状態）とするときの、測定タイミングにおけるＮＤＶＩの指標値と追肥時の施肥量との関係を単純化して示している。図３で示すように、低品質割合を５％とするときの追肥時の施肥量と、測定タイミングにおけるＮＤＶＩの指標値との間には、ＮＤＶＩの指標値が大きくなるほど施肥量が小さくなる関係がある。このように所定のカテゴリについての農作物の状態を目標状態とするときの、測定タイミングにおけるＮＤＶＩの指標値と追肥時の施肥量とには強い相関関係がある。

本実施形態では、施肥量導出計算式は、蓄積された過去の実際の追肥時の施肥量および測定タイミングにおけるＮＤＶＩの指標値の組み合わせと、所定のカテゴリについての農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングにおけるＮＤＶＩの指標値を説明変数とする重回帰分析を行って求められた重回帰式に基づいて求められる。ただし、本実施形態に係る施肥量導出計算式は、日射量レベルが中レベルのときに行われた撮影結果に基づいて導出されるＮＤＶＩの指標値が入力されることを前提とする計算式とされ、このようなＮＤＶＩの指標値が入力されたときに、適切な施肥量を出力するよう調整される。

詳述すると、事前に、実験や既存の圃場に対する観察により、追肥時の施肥量および測定タイミングのＮＤＶＩの指標値の組み合わせと、所定のカテゴリについての農作物の状態との関係を示すデータ、つまり、追肥時の施肥量が○○であり、測定タイミングのＮＤＶＩの指標値が○○であった場合に、所定のカテゴリについての農作物の状態は○○であったというデータが蓄積される。このときの追肥時のＮＤＶＩの指標値は、上述した低レベル時ＮＤＶＩ値変換式や高レベル時ＮＤＶＩ値変換式を利用して、撮影時の日射量レベルが中レベルであったとした場合のＮＤＶＩの指標値に適切に変換される。所定のカテゴリについての農作物の状態は、様々な観点から定義可能である。

そして蓄積されたデータについて、所定のカテゴリについての農作物の状態を目的変数とし、追肥時の施肥量および測定タイミングのＮＤＶＩの指標値を説明変数とする重回帰分析が行われ、農作物の状態を目標状態とするための追肥時の施肥量と測定タイミングのＮＤＶＩの指標値との関係が定義された重回帰式が生成される。例えば農作物の状態が上述した低品質割合である場合には、低品質割合を目標状態（一例として５％）とするための、追肥時の施肥量と測定タイミングのＮＤＶＩの指標値との関係が定義された重回帰式が生成される。そして重回帰式に基づいて、ＮＤＶＩの指標値を入力とし、重回帰式により表されるＮＤＶＩの指標値と施肥量との関係から決定される施肥量を出力とする施肥量導出計算式が生成される。このようにして生成される施肥量導出計算式の出力は、測定タイミングのＮＤＶＩの指標値が入力された値である場合に農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量に相当する。なお重回帰式を、稲の生育時期毎、例えば出穂前３０日、２５日、２０日、１５日、１０日、５日といったように作成するようにしてもよい。この場合において、生育診断する時期が出穂前何日なのかを判定するには、作付け地点、品種、移植日、苗の葉齢及び苗姿、気温や、日長時間等から推定したり、幼穂の長さから推定したりすることができる。

なお本実施形態の施肥量導出計算式の内容は、あくまで単純化した一例である。施肥量導出計算式は、測定タイミングのＮＤＶＩの指標値を入力の１つとし、農作物の状態を目標状態とするために必要かつ適切な追肥時の施肥量を出力する計算式であれば、どのようなものであってもよい。一例として測定タイミングのＮＤＶＩ以外に、農作物の状態に影響を与える要素の要素値を入力とし、当該要素を加味して追肥時の施肥量を出力する式であってもよい。当該要素は、一例として、特定の期間の気温（到来していない期間（例えば登熟期）については予想気温）や、特定の期間における雨量等である。また、施肥量導出計算式は重回帰分析により生成されるものである必要はなく、施肥量導出計算式の生成には様々な機械学習（一例としてニューラルネットワークを用いた機械学習）を応用可能である。

さて１０個の対象圃場のそれぞれについて最終的な追肥時施肥量を導出した後、農作物関連値導出部１２は、表示装置２を制御して、１０個の対象圃場のそれぞれについて識別情報と対応付けて、導出した追肥時施肥量を示す情報を表示装置２の表示領域に表示する。ユーザは、表示装置２を参照することにより、１０個の対象圃場のそれぞれについて農作物関連値導出装置１が導出した追肥時施肥量を認識でき、追肥の際に実際に供給する施肥量を決定するときの参考とすることができる。

以上の通り、本実施形態に係る農作物関連値導出装置１は、対象圃場の上空ＮＤＶＩ値（第１実際値）を取得すると共に、上空ＮＤＶＩ値を、対象圃場に対するドローンに搭載されたセンサからの撮影（遠隔からの観測）が行われたときの日射量（状況、日光の状況）に基づいて、日射量レベルが中レベル（所定の状況、所定の日光の状況）で撮影が行われた場合の値となるように補正して補正ＮＤＶＩ値（補正値）を導出する。そして、日射量レベルが中レベルのときに撮影が行われた場合（所定の状況で観測が行われた場合）のＮＤＶＩの指標値を入力とし追肥時の施肥量（農作物関連値）を出力とする施肥量導出計算式（第１計算式）に補正ＮＤＶＩ値を入力することによって、対象圃場への追肥時の施肥量を導出する。

この構成によれば、追肥時の施肥量を導出するにあたって、圃場に対する遠隔からの観測により得られた上空ＮＤＶＩ値がそのまま施肥量導出計算式に入力されるのではなく、日射量レベルが中レベルのときに撮影が行われた場合の値となるように補正された上で入力される。このため、撮影が行われたときの日射量に応じて実際値の値が変化するという中で、日射量レベルが中レベルという統一された状況で観測が行われた場合の値となるように上空ＮＤＶＩ値を補正し、補正によって生成された補正ＮＤＶＩ値を施肥量導出計算式に入力することができ、これにより撮影時の日射量に応じた上空ＮＤＶＩ値の変化に起因する精度の低下を抑制した状態で施肥量を導出することができる。

次に本実施形態に係る農作物関連値導出装置１の動作（特に、情報提供指示があった場合の動作）について図４のフローチャートを用いて説明する。以下の説明では、農作物関連値導出装置１は、複数個の圃場を対象として追肥時施肥量を決定するものとする。

図４で示すように、第１実際値取得部１０は、入力装置３に対して情報提供指示があったことを検出すると、上空データに基づいて、対象圃場のそれぞれの上空ＮＤＶＩ値を取得する（ステップＳＡ１）。対象圃場のそれぞれの上空ＮＤＶＩ値は、第１実際値取得部１０から補正部１１に対して出力される。補正部１１は、記憶部１３に記憶された撮影時日射量情報を取得し、当該情報に基づいて撮影時日射量の日射量レベルを判定する（ステップＳＡ２）。

次いで、補正部１１は、日射量レベルに応じた方法で対象圃場毎の補正ＮＤＶＩ値を導出する（ステップＳＡ３）。ステップＳＡ３において、補正部１１は、日射量レベルが低レベルの場合には低レベル時ＮＤＶＩ値変換式を用いて補正ＮＤＶＩ値を導出し、日射量レベルが高レベルの場合には高レベル時ＮＤＶＩ値変換式を用いて補正ＮＤＶＩ値を導出し、日射量レベルが中レベルの場合には上空ＮＤＶＩ値をそのまま補正ＮＤＶＩ値とする。対象圃場のそれぞれの補正ＮＤＶＩ値は、補正部１１から農作物関連値導出部１２に対して出力される。農作物関連値導出部１２は、対象圃場のそれぞれについて、補正ＮＤＶＩ値に基づいて、施肥量導出計算式を用いて、追肥時施肥量を導出する（ステップＳＡ４）。農作物関連値導出部１２は、ステップＳＡ４で導出した追肥時施肥量を示す情報を表示する（ステップＳＡ５）。

＜第１実施形態の変形例＞
次に、第１実施形態の変形例について説明する。第１実施形態では、日射量レベルを、低、中、高の３つのレベルに分ける構成としたが、これは説明の簡単のためであり、日射量レベルを２つに分ける構成としてもよく、４つ以上のレベルに分ける構成としてもよい。この場合、何れか１つのレベル（「基準レベル」とする）を「所定日射量レベル」に相当させ、他の日射量レベルに対応する上空ＮＤＶＩ値を、基準レベルに対応する補正ＮＤＶＩ値に変換する必要がある。また、施肥量導出計算式は、基準レベルに対応するＮＤＶＩの指標値を入力したときに、適切な追肥時施肥量を導出する計算式とされる。

また日射量をレベル分けせずに、補正部１１が、上空ＮＤＶＩ値を、所定日射量（レベルのように幅を持った概念ではなく、ピンポイントの値）で撮影が行われた場合の値となるように補正する構成としてもよい。この場合、例えば、撮影時日射量および上空ＮＤＶＩ値を入力とし、補正ＮＤＶＩ値を出力とする計算式（特許請求の範囲の「第２計算式」に相当）が研究やシミュレーションにより事前に導出され、補正部１１は、この計算式を用いて上空ＮＤＶＩ値から補正ＮＤＶＩ値を導出する。また、この場合、施肥量導出計算式は、所定日射量で行われた撮影結果に基づくＮＤＶＩの指標値を入力とし、施肥量を出力とする計算式とされ、農作物関連値導出部１２は、このような施肥量導出計算式を用いて、補正ＮＤＶＩ値から追肥時施肥量を導出する。

また、上記第１実施形態では、「領域に対する遠隔からの観測が行われたときの状況」は、日射量であり、日射量が、特許請求の範囲の「対象となる領域、観測機器のスペック、観測の時期および時間帯、および、観測の位置が同じ場合に、観測される第１指標の指標値に変化を与える自然環境の要素の状況」、「日光の状況」に相当した。しかしながら、自然環境の要素の状況は、日射量に限らず、例えば、気温や湿度、天候の状態が第１指標値に変化を与える場合に、これらを自然環境の要素の状況に含めてもよい。また、自然環境の要素の状況は１つに限らず、複数の組み合わせ（例えば、日射量と温度との組み合わせ）であってもよい。

特に「日光の状況」については、上記第１実施形態の日射量（全天日射量）のほか、全天日射量、散乱光量或いはその割合、直達光量或いはその割合、または、太陽高度の何れか１つまたは２つ以上の組み合わせであってもよい。

また「領域に対する遠隔からの観測が行われたときの状況」は、領域（第１実施形態では圃場）に対する遠隔からの観測が行われたときの時間的な状況であってもよい。この場合、補正部１１は、第１実際値取得部１０により取得された上空ＮＤＶＩ値（第１実際値）を、圃場（領域）に対する遠隔からの観測が行われたときの時間的な状況に基づいて、時間的な状況が所定の時間的な状況で観測が行われた場合の値となるように補正して補正値を導出する。農作物関連値導出部１２は、所定の時間的な状況で観測が行われた場合のＮＤＶＩ（第１指標）の指標値を入力とし追肥時の施肥量（農作物関連値）を出力とする第１計算式に補正部１１により導出された補正値を入力することによって、追肥時の施肥量（農作物関連値）を導出する。また、なり得る時間的な状況ごとに、各時間的な状況で観測が行われた場合のＮＤＶＩの指標値（第１指標の指標値）を、所定の時間的な状況で観測が行われた場合のＮＤＶＩの指標値（第１指標の指標値）に変換する第２計算式が事前に導出され、補正部１１は、第２計算式を用いて補正値を導出する。

時間的な状況は例えば、観測が行われた時期、または、観測が行われた時間帯である。例えば時間的な状況が「観測が行われた時期」である場合、マルチスペクトルセンサによる撮影が行われる時期が固定されていない状況下において（その他の状況は原則的に一致し或いは近似した状態とされる）、補正部１１は、事前に用意された第２計算式を用いて第１実際値取得部１０により取得された上空ＮＤＶＩ値を、予め定められた特定の時期で撮影されたとした場合の値となるように補正して補正値を導出する。また例えば時間的な状況が「観測が行われた時間帯」である場合、マルチスペクトルセンサによる撮影が行われる時間帯が固定されていない状況下において（その他の状況は原則的に一致し或いは近似した状態とされる）、補正部１１は、事前に用意された第２計算式を用いて第１実際値取得部１０により取得された上空ＮＤＶＩ値を、予め定められた特定の時間帯で撮影されたとした場合の値となるように補正して補正値を導出する。

また「領域に対する遠隔からの観測が行われたときの状況」は、領域（第１実施形態では圃場）に対する遠隔からの観測が行われたときの場所的な状況であってもよい。この場合、補正部１１は、第１実際値取得部１０により取得された上空ＮＤＶＩ値（第１実際値）を、圃場（領域）に対する遠隔からの観測が行われたときの場所的な状況に基づいて、場所的な状況が所定の場所的な状況で観測が行われた場合の値となるように補正して補正値を導出する。農作物関連値導出部１２は、所定の場所的な状況で観測が行われた場合のＮＤＶＩ（第１指標）の指標値を入力とし追肥時の施肥量（農作物関連値）を出力とする第１計算式に補正部１１により導出された補正値を入力することによって、追肥時の施肥量（農作物関連値）を導出する。また、なり得る場所的な状況ごとに、各場所的な状況で観測が行われた場合のＮＤＶＩの指標値（第１指標の指標値）を、所定の場所的な状況で観測が行われた場合のＮＤＶＩの指標値（第１指標の指標値）に変換する第２計算式が事前に導出され、補正部１１は、第２計算式を用いて補正値を導出する。

場所的な状況は例えば、観測が行われた場所の緯度である。例えば場所的な状況が「観測が行われた場所の緯度」である場合、マルチスペクトルセンサによる撮影が行われる場所の緯度が固定されていない状況下において（その他の状況は原則的に一致し或いは近似した状態とされる）、補正部１１は、事前に用意された第２計算式を用いて第１実際値取得部１０により取得された上空ＮＤＶＩ値を、予め定められた特定の緯度で撮影されたとした場合の値となるように補正して補正値を導出する。

＜第２実施形態＞
次に第２実施形態について説明する。図５は、本実施形態に係る農作物関連値導出装置１Ａの機能構成例を示すブロック図である。図５で示すように、本実施形態に係る農作物関連値導出装置１Ａは、機能構成として、第１実際値取得部１０Ａと第２実際値取得部２０と補正部１１Ａと農作物関連値導出部１２Ａとを備えている。また、農作物関連値導出装置１Ａは、記憶手段として、記憶部１３Ａを備えている。以下、農作物関連値導出装置１Ａは、第１実施形態に係る農作物関連値導出装置１と同様、情報提供指示があった場合に、対象となる圃場毎の追肥時施肥量を導出し、追肥時施肥量を示す情報を提供する。以下、追肥時施肥量を決定する対象となる圃場が１０個であるものとして、情報提供指示があった後の農作物関連値導出装置１Ａの動作について説明する

第１実際値取得部１０Ａは、入力装置３に対して情報提供指示があったことを検出すると、上空データを取得する。上述の通り、上空データは、追肥を行うタイミングに近いタイミングで、ドローンに搭載されたマルチスペクトルセンサにより、１０個の圃場の稲の群落が上空から一括して撮影されることによって生成された画像データである（ただし、第１実施形態で説明した通り、上空データは、必ずしも「一括して行われた撮影」に基づいて生成される必要はなく、複数回の撮影に基づいて生成されてもよい。）。ただし、第１実施形態では、上空データのための撮影は、必ずしもドローンを用いて行われる必要はないとしたが、本実施形態では、当該撮影は、ドローン或いはドローンに準じた装置のように、高度を変更可能な飛行物に搭載されたセンサにより行われている。上空データのための撮影は、１０個の圃場を一括して撮影可能な程度に十分に高い高度から行われることになる。本実施形態では、上空データのための撮影は、高度「１００メートル」において行われたものとする。上空データのための撮影が行われたときのドローンの高度は、特許請求の範囲の「第１高度」に相当する。上空データを取得した後、第１実際値取得部１０Ａは、上空データに基づいて、第１実施形態に係る第１実際値取得部１０と同様の方法で、対象圃場毎の上空ＮＤＶＩ値を導出し、取得する。

上空ＮＤＶＩ値の導出にあたって使用される上空データは、ドローンに搭載されたマルチスペクトルセンサの撮影結果に基づいて生成されるものである。上述の通り、ドローンを使用し、複数の圃場を一括して撮影することによって上空データを生成することにより、作業にかかる時間を短くすることができ、作業にかかる労力を低減することができる。しかしながら、複数の圃場を一括して撮影できる程度にドローンの高度を高くすることに起因して、導出される上空ＮＤＶＩ値の精度が、高度が低い状態での撮影結果に基づくＮＤＶＩの精度と比較して低下してしまう。同じ自然環境（日射量を含む）の下、同じ圃場を同じスペックのマルチスペクトルセンサによって撮影した場合であっても、撮影が行われたときのドローンの高度が高いほど、上空ＮＤＶＩ値の精度が低下するからである。これは、ドローンの高度が高いほど、観測対象（被写体）たる圃場とマルチスペクトルセンサとの離間距離が大きくなり、当該離間距離が大きくなると、解像度が低下し（＝１つの圃場に割り当てられる画素数が低下し）、また、センサが検出する反射光が弱まることに起因する。なお、本実施形態では、上空データの生成にあたり、１０個の圃場を高度「１００メートル」の位置にて撮影したこととしているが、これはあくまで説明の単純化を考慮した一例である。撮影位置の高度や、１回の撮影により撮影する範囲、その他の撮影の具体的な態様は、状況に応じて適切に定められるべきものである。また、広い領域（一例として１００ｈａを超える領域）を対象とする際に、十分に高い高度（第１高度）から、複数回にわたって領域を分割して撮影し、撮影結果に基づく各画像データをオルソ補正により結合しこれを上空ＮＤＶＩ値とするようにしてもよい。

さて、第１実際値取得部１０Ａは、対象圃場毎の上空ＮＤＶＩ値を導出し、取得した後、取得した上空ＮＤＶＩ値を各圃場の識別情報と対応付けて補正部１１Ａに出力する。

第２実際値取得部２０は、入力装置３に対して情報提供指示があったことを検出すると、記憶部１３Ａに記憶された３個の対象圃場のそれぞれについての低空データのそれぞれを取得する。ある１つの対象圃場についての低空データとは、ドローンがその対象圃場に十分に近い特定の高度に位置した状態でマルチスペクトルセンサにより実行された撮影により生成された画像データである。低空データに係る高度（特定の高度）は、１０個の対象圃場全体を撮影可能な程度に高い必要はなく、撮影対象の１つの対象圃場の一部が撮影できる高さであればよいため、上空データに係る高度よりも十分に低い。本実施形態では、低空データに係る高度は、「３０メートル」であるものとする。低空データは、３個の対象圃場のそれぞれが低空で撮影されることによって、３個の対象圃場についてのみ生成され、それ以外の７個の対象圃場についての低空データは生成されない。

ここで、３個の対象圃場は、１０個の対象圃場から稲の生育量の相違に着目して選択された３個の対象圃場であり、本実施形態では、生育量が最も多い圃場（以下「最大圃場」という）と、生育量が最も少ない圃場（以下「最小圃場」という）と、生育量が中間程度の圃場（以下「中間圃場」という）との３個の圃場が選択される。中間圃場は例えば、１０個の対象圃場を生育量の多さで昇順に並べたときの５番目または６番目に位置する対象圃場とされる。また、中間圃場は例えば、最大圃場の生育量と最小圃場の生育量との中間の生育量に、生育量が最も近い圃場とされる。なお、各対象圃場の生育量は、例えば既存の技術による観測によって把握され、また例えば対象圃場毎に生育量が計画的に調整されている場合には、生育量の計画に基づいて把握される。以下、最大圃場、最小圃場および中間圃場を総称して「補正用圃場」という。

なお、補正用圃場の個数は、補正用圃場を含む対象圃場全体の個数よりも少ないものとされる。本実施形態では、補正用圃場の個数は３個となるため、対象圃場は４個以上とされる。ただし、後に明らかとなる通り、補正用圃場の個数に比べて対象圃場の個数が十分に多い場合に、ユーザは本実施形態に係る農作物関連値導出装置１Ａによるメリットを効果的に得ることができる。このため、対象圃場の個数は、補正用圃場の個数よりも十分に多い方が好ましい。

記憶部１３Ａに記憶された３つの補正用圃場のそれぞれの低空データのそれぞれを取得した後、第２実際値取得部２０は、取得した低空データのそれぞれを分析し、周知の方法で補正用圃場毎のＮＤＶＩ（第１指標）の指標値の実際の値を導出し、取得する。以下、低空データに基づいて導出された補正用圃場のＮＤＶＩの指標値の実際の値を「低空ＮＤＶＩ値」という。低空ＮＤＶＩ値は、特許請求の範囲の「第２実際値」に相当する。上述した通り、低空ＮＤＶＩ値を生成する基礎となる低空データのための撮影は、上空データに係る撮影位置の高度（本実施形態では「１００メートル」）と比較して十分に低い高度（本実施形態では「３０メートル」）で行われたものである。従って、低空データに基づいて導出される低空ＮＤＶＩ値は、上空ＮＤＶＩ値と比較して、センサが検出する圃場の反射光の弱まりや解像度の低下に起因したＮＤＶＩの精度の低下が抑制されており、上空ＮＤＶＩ値よりも精度が高いと言える。

第２実際値取得部２０は、３つの補正用圃場の低空ＮＤＶＩ値を導出し、取得した後、取得した低空ＮＤＶＩ値を各圃場の識別情報と対応付けて補正部１１Ａに出力する。

補正部１１Ａは、第１実際値取得部１０Ａから１０個の対象圃場のそれぞれの上空ＮＤＶＩ値を入力し、第２実際値取得部２０から３個の補正用圃場のそれぞれの低空ＮＤＶＩ値を入力すると、補正用計算式生成処理を実行する。以下、補正用計算式生成処理について詳述する。

補正用計算式生成処理において、補正部１１Ａは、３個の補正用圃場のそれぞれについて、低空ＮＤＶＩ値および上空ＮＤＶＩ値を取得する。次いで、補正部１１Ａは、横軸に上空ＮＤＶＩ値を取り縦軸に低空ＮＤＶＩ値を取った２次元空間に、３個の補正用圃場のそれぞれについての上空ＮＤＶＩ値と低空ＮＤＶＩ値との組み合わせを示す点をプロットする。次いで、補正部１１Ａは、プロットした各点の近似曲線を求める。補正部１１Ａは、近似曲線を表す式を、上空ＮＤＶＩ値を入力値として入力して、低空ＮＤＶＩ値を出力値として出力する計算式へと変換し、これを補正用計算式（特許請求の範囲の「第３計算式」に相当）とする。

例えば、図６で示すように、１０個の圃場のそれぞれの識別情報がＨ０１、Ｈ０２、Ｈ０３・・・Ｈ１０であり、そのうち識別番号がＨ０１、Ｈ０２およびＨ０９の３個の圃場が補正用圃場であったとする。また、識別番号がＨ０１～Ｈ１０の圃場についての上空ＮＤＶＩ値が図に示す通りであり、識別番号がＨ０１、Ｈ０２、Ｈ０９の圃場について低空ＮＤＶＩ値が図に示す通りであったとする。この場合において、横軸に上空ＮＤＶＩ値を取り、縦軸に低空ＮＤＶＩ値を取った２次元空間に、識別番号がＨ０１、Ｈ０２およびＨ０９の３個の補正用圃場のそれぞれについて、低空ＮＤＶＩ値と上空ＮＤＶＩ値との組み合わせを示す点をプロットし、またその近似曲線を求めると図７で示す通りになる。そして近似曲線に基づいて求められる補正用計算式は、「出力値（低空ＮＤＶＩ値）＝入力値（上空ＮＤＶＩ値）＋０．１０」となる。

以下、高度１００メートルからの撮影結果により得られたＮＤＶＩの指標値を「高高度指標値」といい、高度３０メートルからの撮影結果により得られたＮＤＶＩの指標値を「低高度指標値」という。このような方法で生成される補正用計算式に一の高高度指標値を入力することによって得られる出力値は、補正用圃場を対象とした上空ＮＤＶＩ値と低空ＮＤＶＩ値との相関関係に準じて当該一の高高度指標値を仮想的に低高度指標値へと変換したものということができる。そして、上述の通り低高度指標値は、高高度指標値よりも精度が高く、補正用圃場を対象とした上空ＮＤＶＩ値と低空ＮＤＶＩ値との相関関係に準じて補正された出力値は、入力値よりも精度が高いＮＤＶＩの指標値ということができる。

ここで、上述したように本実施形態では、補正用圃場として、生育量が最も多い圃場と生育量が最も少ない圃場と生育量が中間程度の圃場との３個の圃場が選択される。これは以下の理由による。すなわち、周知の通り、圃場における生育量は、圃場についてのＮＤＶＩの指標値に影響を与え、生育量とＮＤＶＩの指標値とは正の相関関係がある。そして、上記の方針で補正用圃場を選択することにより、横軸に上空ＮＤＶＩ値を取り、縦軸に低空ＮＤＶＩ値を取った２次元空間で適度に離間した３つの点をプロットすることが可能となり、各点の位置に基づいて導出される近似曲線の妥当性を高めることができるからである。

以上の処理が補正用計算式生成処理である。補正用計算式生成処理により補正用計算式を生成した後、補正部１１Ａは、以下の処理を実行する。すなわち、補正部１１Ａは、補正用圃場については、低空ＮＤＶＩ値をそのまま、補正ＮＤＶＩ値（特許請求の範囲の「補正値」に相当）とする。更に補正部１１Ａは、補正用圃場ではない７個の圃場（以下「被補正圃場」という）のそれぞれについては、補正用計算式に上空ＮＤＶＩ値を入力し、その出力値を補正ＮＤＶＩ値とする。

例えば、図６を参照し、７個の被補正圃場について、上空ＮＤＶＩ値が図６の通りとし、補正用計算式が「出力値＝入力値＋０．１０」であるとすると、これら７個の被補正圃場については、補正ＮＤＶＩ値は図８に示す通りとなる。識別番号がＨ０３の被補正圃場に着目すると、当該被補正圃場の上空ＮＤＶＩ値は「０．４０」であり、これを入力値として補正用計算式に入力することによって出力される出力値は、「０．４０＋０．１０＝０．５０」となる。

対象圃場（３個の補正用圃場＋７個の被補正圃場）のそれぞれについて補正ＮＤＶＩ値を導出した後、補正部１１Ａは、導出した補正ＮＤＶＩ値を各圃場の識別情報と対応付けて農作物関連値導出部１２Ａに出力する。

農作物関連値導出部１２Ａは、補正部１１Ａから対象圃場のそれぞれについての補正ＮＤＶＩ値を入力すると、対象圃場のそれぞれについて、施肥量導出計算式に補正ＮＤＶＩ値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出する。施肥量導出計算式は、第１実施形態で説明したように、ＮＤＶＩの指標値を入力し、測定タイミングのＮＤＶＩの指標値が入力された値であった場合に、農作物の状態を目標状態とするために必要かつ適切な追肥時施肥量を出力する計算式である。ただし、本実施形態に係る施肥量導出計算式は、ドローンの高度が「３０メートル」のときに行われた撮影結果に基づいて導出される低高度指標値が入力されることを前提とする計算式であり、このような低高度指標値が入力されたときに、適切な施肥量を出力するように調整されている。

１０個の対象圃場のそれぞれについて最終的な追肥時施肥量を決定した後、農作物関連値導出部１２Ａは、表示装置２を制御して、１０個の対象圃場のそれぞれについて識別情報と対応付けて、導出した追肥時施肥量を示す情報を表示装置２の表示領域に表示する。ユーザは、表示装置２を参照することにより、１０個の対象圃場のそれぞれについて農作物関連値導出装置１Ａが導出した追肥時施肥量を認識でき、追肥の際に実際に供給する施肥量を決定するときの参考とすることができる。

以上の通り、本実施形態に係る農作物関連値導出装置１Ａは、対象圃場の上空ＮＤＶＩ値（圃場に対する遠隔からの観測により得られた第１実際値）を取得すると共に、上空ＮＤＶＩ値を、対象圃場に対するドローン（高度を変更可能な飛行物）からの撮影（遠隔からの観測）が行われたときの高度（状況）に基づいて、高度が３０メートルという状況（所定の状況）で撮影が行われた場合の値となるように補正して補正ＮＤＶＩ値（補正値）を導出する。そして、農作物関連値導出装置１Ａは、高度が３０メートルのときに撮影が行われた場合（所定の状況で観測が行われた場合）のＮＤＶＩの指標値を入力とし追肥時の施肥量を出力とする施肥量導出計算式に補正ＮＤＶＩ値を入力することによって、対象圃場への追肥時の施肥量を導出する。

この構成によれば、追肥時の施肥量を導出するにあたって、上空ＮＤＶＩ値がそのまま施肥量導出計算式に入力されるのではなく、高度が３０メートルのときに撮影が行われた場合の値となるように補正された補正ＮＤＶＩ値が入力される。このため、観測が行われたときのドローンの高度に応じて上空ＮＤＶＩ値の精度の高低が変化するという中で、高い高度（得られるＮＤＶＩの精度が低い状況）で行われた撮影に基づいて導出された実際値を、低い高度（得られるＮＤＶＩの精度が高い状況）で観測が行われた場合の値となるように補正し、補正によって生成された補正ＮＤＶＩ値を施肥量導出計算式に入力することができ、これにより撮影時の状況に応じた実際値の精度の高低の変化に起因する精度の低下を抑制した状態で施肥量を導出することができる。

ここで、全ての圃場について、低空で撮影を行って低空データを生成し、低空データに基づいて低空ＮＤＶＩ値を導出するようにすれば、撮影時の高度が高いことに起因するＮＤＶＩの指標値の精度の低下による悪影響は抑制できる。しかしながら、この場合、全ての圃場のそれぞれ（上記実施形態では、１０個としたが、当然、１０個よりも多い場合もある）について、低空の所定の位置で撮影を行う必要があり、作業が煩雑化し、作業に要する時間が長期化してしまう。一方で、本実施形態によれば、限定された個数（本実施形態では３個）の圃場についてのみ低空の撮影を行えばよく、作業の煩雑化および作業に要する時間の長期化を抑制しつつ、撮影時の高度が高いことに起因するＮＤＶＩの指標値の精度の低下による悪影響を抑制することが可能である。

次に本実施形態に係る農作物関連値導出装置１Ａの動作（特に、情報提供指示があった場合の動作）について図９のフローチャートを用いて説明する。以下の説明では、農作物関連値導出装置１Ａは、Ｎ（ただし、Ｎ≧２）個の対象圃場を対象として追肥時施肥量を決定するものとし、Ｎ個の対象圃場のうち、Ｍ（ただし、Ｍ＜Ｎ）個が補正用圃場であるものとする。

図９で示すように、ステップＳＢ１において、第１実際値取得部１０Ａは、入力装置３に対して情報提供指示があったことを検出すると、上空データを取得し、Ｎ個の対象圃場毎の上空ＮＤＶＩ値を取得する。またステップＳＢ１において、第２実際値取得部２０は、入力装置３に対して情報提供指示があったことを検出すると、Ｍ個の補正用圃場毎の低空ＮＤＶＩ値を取得する。ステップＳＢ１において第１実際値取得部１０Ａにより取得された上空ＮＤＶＩ値のそれぞれ、および、第２実際値取得部２０により取得された低空ＮＤＶＩ値のそれぞれは、補正部１１Ａに出力される。

補正部１１Ａは、第１実際値取得部１０ＡからＮ個の対象圃場のそれぞれの上空ＮＤＶＩ値を入力し、第２実際値取得部２０からＭ個の補正用圃場のそれぞれの低空ＮＤＶＩ値を入力すると、補正用計算式生成処理を実行し、補正用計算式を生成する（ステップＳＢ２）。次いで、補正部１１Ａは、補正用計算式を用いて、対象圃場のそれぞれについて、補正ＮＤＶＩ値を導出する（ステップＳＢ３）。ステップＳＢ３で導出された対象圃場毎の補正ＮＤＶＩ値は、農作物関連値導出部１２Ａに出力される。農作物関連値導出部１２Ａは、補正部１１Ａから対象圃場のそれぞれについての補正ＮＤＶＩ値を入力すると、対象圃場のそれぞれについて、施肥量導出計算式に補正ＮＤＶＩ値を入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出する（ステップＳＢ４）。ついて、農作物関連値導出部１２Ａは、ステップＳＢ４で導出した追肥時施肥量を示す情報を表示する（ステップＳＢ５）。

＜第２実施形態の変形例＞
次に、第２実施形態の変形例について説明する。上記第２実施形態では、撮影時のドローンの高度の影響を加味して追肥時の施肥量を導出する構成であったが、当然、ドローンの高度に加え、日射量を加味して施肥量を導出する構成でもよい。この場合、農作物関連値導出装置１Ａは、以下の処理を実行する。

すなわち、補正部１１Ａは、第１実際値取得部１０Ａにより取得された上空ＮＤＶＩ値に基づいて、第２実施形態で説明した方法で、対象圃場毎の「高度３０メートル（所定高度）で撮影が行われた場合のＮＤＶＩの指標値」（第２実施形態の補正ＮＤＶＩ値に相当）を導出する。説明の便宜上、ここで導出されたＮＤＶＩを中間ＮＤＶＩ値とする。更に、補正部１１Ａは、撮影が行われたときの日射量レベルに基づいて、所定日射量レベル（第１実施形態の中レベルに相当するレベル）或いは所定日射量で観測が行われた場合の値となるように中間ＮＤＶＩ値を補正して補正ＮＤＶＩ値を導出する。施肥量導出計算式は、高度３０メートル（所定高度）において所定日射量レベル或いは所定日射量のときに観測が行われた場合のＮＤＶＩの指標値を入力とし追肥時の施肥量を出力とする計算式とされる。そして、農作物関連値導出部１２Ａは、施肥量導出計算式に補正部１１Ａにより導出された補正ＮＤＶＩ値を入力することによって、圃場への追肥時の施肥量を導出する。なお、この構成において、自然環境の要素を日射量以外の要素（或いは複数の要素の組み合わせ）としてもよいことは第１実施形態の変形例で述べた通りである。

また、第２実施形態において、対象圃場のうち、補正用圃場を選択するときのルールは、第２実施形態で例示したものに限らない。例えば、生育量以外の要素に基づいて補正用圃場を選択するようにしてもよい。ただし補正用圃場の個数が少ないほど、作業に要する時間が短くなるものの、補正用計算式を生成するときのサンプル数が少なくなるため補正用計算式の精度が低くなり、逆に、補正用圃場の個数が多いほど補正用計算式の精度が高くなる一方、作業に要する時間が長くなる。補正用圃場の個数は、このことを考慮して定められる必要がある。

また、第２実施形態において、上空ＮＤＶＩ値を低空ＮＤＶＩ値によって補正する方法は、上記実施形態で例示した方法に限られない。すなわち、上空ＮＤＶＩ値の値に低空ＮＤＶＩ値を反映させることにより、上空ＮＤＶＩ値の精度が向上する方法であれば、どのような方法であってもよい。一例として、農作物関連値導出部１２Ａが、補正用圃場の上空ＮＤＶＩ値と低空ＮＤＶＩ値との相関関係を分析することによって、上空ＮＤＶＩ値を低空ＮＤＶＩ値に近づけるための加算値または減算値を固定値として導出し、被補正圃場の上空ＮＤＶＩ値に対して加算値を加算しまたは減算値を減算して補正ＮＤＶＩ値を求めるようにしてもよい。

また、第２実施形態では、農作物関連値導出装置１Ａは、３個（Ｍ個）の補正用圃場については、３個の補正用圃場についての低空ＮＤＶＩ値をそのまま補正ＮＤＶＩ値として施肥量導出計算式に入力し、出力値として得られる施肥量を最終的な追肥時施肥量として導出した。この点に関し、３個の補正用圃場についても、上空ＮＤＶＩ値から補正用計算式を用いて補正ＮＤＶＩ値を取得し、補正ＮＤＶＩ値を施肥量導出計算式に入力し、当該計算式の出力値を最終的な追肥時施肥量として導出する構成でもよい。

以上、本発明の一実施形態（変形例を含む）を説明したが、上記各実施形態は、本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその要旨、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。

例えば上記各実施形態では、圃場毎に、追肥時の施肥量（農作物関連値）を導出した。この点について、圃場毎に農作物関連値を導出する構成以外の構成でもよいことは勿論である。例えば、１つの圃場を複数の領域に区分けし、領域毎に追肥時の施肥量（農作物関連値）を導出する構成でもよい。

また例えば、第１実施形態について以下の構成としてもよい。すなわち補正部１１は、第１実際値取得部１０により取得された上空ＮＤＶＩ値（第１実際値）を、観測が行われたときの日光の状況、時間的な状況、場所的な状況および飛行物の高度の２つ以上の組み合わせに基づいて、所定の日光の状況、所定の時間的な状況、所定の場所的な状況および所定高度のうち、観測が行われたときの状況についての組み合わせに対応する２つ以上の組み合わせで観測が行われた場合の値となるように補正して補正ＮＤＶＩ値（補正値）を導出する。農作物関連値導出部１２は、所定の日光の状況、所定の時間的な状況、所定の場所的な状況および所定高度の２つ以上の組み合わせのときに観測が行われた場合のＮＤＶＩ（第１指標）の指標値を入力とし追肥時の施肥量（農作物関連値）を出力とする計算式に補正部１１により導出された補正値を入力することによって、追肥時の施肥量を導出する。

また第１実施形態では、農作物関連値導出装置１は、農作物関連値として稲（農作物）についての追肥時の施肥量（所定のタイミングにおける施肥量）を導出した。しかしながら、農作物は稲に限られるものではなく、また、農作物関連値は「追肥時の施肥量」に限られるものではない。まず、農作物関連値として追肥時の施肥量を導出する構成の場合、本発明は、栽培に際して追肥が行われる農作物に広く適用可能である。一例として小麦、大麦、大豆その他の穀物や、野菜、果物等に適用可能である。第２実施形態についても同様である。

農作物関連値は、農作物栽培領域で栽培される農作物に関連する値であって、第１指標の指標値（ＮＤＶＩに限られない。また、複数種類の指標の指標値であってもよい）から導出可能なものであればよい。

例えば、農作物関連値に、所定のタイミングにおいて供給する堆肥、土壌改良資材、植物生育調整剤或いは薬剤（各種の農薬を含む）の量が含まれていてもよい。所定のタイミングにおいて供給する堆肥の量とは、所定のカテゴリについての農作物（農作物の収穫物も含む）の状態を目標状態とするために、所定のタイミングにおいて施すべき堆肥の量のことである。所定のタイミングにおいて供給する土壌改良資材、植物生育調整剤或いは薬剤の量についても同様である。

この場合において例えば、農作物が稲であり、農作物関連値が所定のタイミングにおける堆肥の量である場合、施肥量導出計算式は、予め定められた測定タイミングにおけるＮＤＶＩの指標値を入力し、所定のタイミングにおける堆肥の量を出力する式（以下「堆肥量導出計算式」という）とされる。つまり堆肥量導出計算式は、測定タイミングにおけるＮＤＶＩの指標値が○○であったら、所定のタイミングにおける適切な堆肥の量は○○である、ということを計算する式である。測定タイミングにおけるＮＤＶＩの指標値と所定のタイミングにおける堆肥の量とには相関関係があることが分かっており（土壌改良資材、植物生育調整剤および薬剤の量も同様）、この堆肥量導出計算式は、測定タイミングにおけるＮＤＶＩの指標値と、所定のタイミングにおける堆肥の量との相関関係に基づいて事前に導出される。

また例えば農作物関連値を、収穫物の収量（単位面積当たりの予測収量を意味する）とすることができる。この場合、施肥量導出計算式は、出穂前の所定のタイミングにおけるＮＤＶＩの指標値を入力し、稲の収量を出力する式（以下「予測収量導出計算式」という）とされる。つまり予測収量導出計算式は、出穂前の所定のタイミングにおけるＮＤＶＩの指標値が○○であったら、稲の収量は○○と予測される、ということを計算する式である。発明者らの研究により、（出穂後ではなく）出穂前の所定のタイミングにおけるＮＤＶＩの指標値と稲の収量とには特に強い相関関係があることが分かっており、この予測収量導出計算式は、出穂前の所定のタイミングにおけるＮＤＶＩの指標値と稲の収量との相関関係に基づいて事前に導出される。

また、農作物関連値を、農作物の現在或いは将来の品質（外観品質のみならず、所定の栄養素の状態（例えばタンパク質の含有率）や、含有水分の状況などの成分的、内部的な品質も含む概念）とすることができる。農作物の現在或いは将来の品質は例えば、農作物が稲である場合において、現在或いは将来の籾水分や、現在或いは将来の玄米水分、現在或いは将来の玄米のタンパク質含有率である。例えば農作物の現在或いは将来の品質は、農作物が稲である場合において、収穫時における低外観品質の収穫物の割合（或いは整粒の割合）である。

また、農作物関連値を、「農作物が所定の状態となると予測される時期」とすることができる。「農作物が所定の状態となると予測される時期」は例えば、農作物が収穫可能な状態となる時期や、農作物が稲である場合において、出穂する時期である。例えば、農作物が稲であり、農作物関連値が「収穫時期」（稲を収穫する適切な時期）である場合、施肥量導出計算式は、測定タイミングのＮＤＶＩの指標値を入力し、収穫時期を出力する式とされる。ＮＤＶＩの指標値と収穫時期とには相関関係があることが分かっており、当該所定の計算式は、ＮＤＶＩの指標値と収穫時期との相関関係に基づいて事前に生成される。

また、農作物関連値を、「農作物についての被害の状況」とすることができる。「農作物についての被害の状況」は例えば、倒伏の状況（程度）や、枯れ被害の状況（程度）、病害の状況(程度)、害虫による被害の状況（程度）等である。

また、農作物関連値を「農作物に影響を与える環境の状態」とすることができる。「農作物に影響を与える環境の状態」は例えば、圃場における病害虫の発生の程度や、圃場における雑草の発生の程度等である。

また、記憶部１３、１３Ａは、農作物関連値導出装置１とは異なる外部記憶装置が有する構成でもよい。また、農作物関連値導出装置１、１Ａを、インターネット上のサーバとして構成し、インターネットを介してクライアントから追肥時施肥量を示す情報の提供の要請があったときに、上述した手法で追肥時施肥量を導出し、クライアントに提供する構成としてもよい。また、農作物関連値導出装置１、１Ａは、単体のコンピュータである必要はなく、複数のコンピュータにより構成されていてもよい。例えば、インターネットを介して接続された端末とクラウドサーバとが協働して農作物関連値導出装置１、１Ａとして機能し、端末が適宜クラウドサーバと協働して処理を実行する構成でもよい。

また第１実施形態では、ドローンを用いた上空からのセンシングにより得られるＮＤＶＩの値が第１指標の指標値であった。この点に関し、第１指標の指標値は、ＮＤＶＩの値に限らず、ＮＤＶＩ以外の他の指標の指標値であってもよい。また、第１指標の指標値は、衛星からの観測による衛星データに基づいて導出された値、或いは、ドローンを用いた上空からのセンシングにより得られるデータと前記衛星データとを少なくとも利用して導出された値であってもよい。

１、１Ａ 農作物関連値導出装置
１０、１０Ａ 第１実際値取得部
１１、１１Ａ 補正部
１２、１２Ａ 農作物関連値導出部
２０ 第２実際値取得部

Claims (6)

1. 農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する農作物関連値導出装置であって、
   前記農作物栽培領域を構成する領域に対する、高度を変更可能な飛行物からの観測により得られた第１指標の指標値の実際の値である第１実際値を取得する第１実際値取得部と、
   前記第１実際値取得部により取得された第１実際値を、観測が行われたときの前記飛行物の高度に基づいて、所定高度で観測が行われた場合の値となるように補正して補正値を導出する補正部と、
   第１指標の指標値を入力とし前記農作物関連値を出力とする計算式であって、前記所定高度で行われた観測によって得られた第１指標の指標値を入力したときに適切な値の前記農作物関連値を出力するように調整された計算式である第１計算式により前記農作物関連値を導出する農作物関連値導出部とを備え、
   前記第１実際値取得部は、Ｎ（ただし、Ｎ≧２）個の前記領域のそれぞれについて、Ｎ個の前記領域に対して行われた、前記所定高度よりも高い第１高度からの観測により得られた第１実際値を取得し、
   Ｎ個の前記領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の前記領域のそれぞれについて、Ｍ個の前記領域に対して個別に行われた前記所定高度からの観測により得られた第１指標の指標値の実際の値である第２実際値を取得する第２実際値取得部を更に備え、
   Ｎ個の前記領域のうち、第２実際値が取得されていない一の前記領域の前記農作物関連値を導出する際、前記補正部は、前記第１実際値取得部により取得された当該一の前記領域についての第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて前記所定高度で観測が行われた場合の値となるように補正して補正値を導出し、前記農作物関連値導出部は、前記補正部により導出された補正値を前記第１計算式に入力することによって前記農作物関連値を導出する
   ことを特徴とする農作物関連値導出装置。
2. 前記補正部は、
   Ｍ個の前記領域のそれぞれについて前記第１実際値取得部により取得された第１実際値と、Ｍ個の前記領域のそれぞれについて前記第２実際値取得部により取得された第２実際値との相関関係を分析し、分析結果に基づいて、前記第１高度からの観測により得られた第１指標の指標値を前記所定高度から観測した場合に得られる値に変換する第３計算式を生成し、
   Ｎ個の前記領域のうち、第２実際値が取得されていない一の前記領域についての補正値を導出する際、前記第３計算式に当該一の前記領域についての第１実際値を入力し、出力された値を補正値とする
   ことを特徴とする請求項１に記載の農作物関連値導出装置。
3. 前記農作物関連値は、追肥時の施肥量であり、
   前記農作物関連値導出部は、蓄積された過去の実際の追肥時の施肥量および追肥時において前記所定高度で観測が行われたときの第１指標の指標値の組み合わせと、農作物の状態との関係について、農作物の状態を目的変数とし、追肥時の施肥量および第１指標の指標値を説明変数とする重回帰分析を行って求められた重回帰式に基づく計算式を前記第１計算式として利用して追肥時の施肥量を導出する
   ことを特徴とする請求項１または２に記載の農作物関連値導出装置。
4. 前記農作物関連値導出部は、前記農作物関連値として、所定のタイミングにおいて供給する肥料、堆肥、土壌改良資材、植物生育調整剤或いは薬剤の量、収穫物の収量、農作物の現在或いは将来の品質、農作物が所定の状態となると予測される時期、農作物についての被害の状況、および、農作物に影響を与える環境の状態の何れか１つを少なくとも導出する
   ことを特徴とする請求項１または２に記載の農作物関連値導出装置。
5. 第１指標の指標値は、前記飛行物を用いた上空からのセンシングにより得られるＮＤＶＩの値である
   ことを特徴とする請求項１から４の何れか１項に記載の農作物関連値導出装置。
6. 農作物栽培領域を構成する領域のうち、Ｎ（ただし、Ｎ≧２）個の前記領域のそれぞれについて、Ｎ個の前記領域に対して行われた、高度を変更可能な飛行物からの所定高度よりも高い第１高度からの観測により得られた第１指標の指標値の実際の値である第１実際値を取得するステップと、
   Ｎ個の前記領域のうち、Ｍ（ただし、Ｍ＜Ｎ）個の前記領域のそれぞれについて、Ｍ個の前記領域に対して個別に行われた前記所定高度からの観測により得られた第１指標の指標値の実際の値である第２実際値を取得するステップと、
   Ｎ個の前記領域のうち、第２実際値が取得されていない一の前記領域について、前記農作物栽培領域で栽培される農作物に関連する農作物関連値を導出する際に、当該一の前記領域について取得した第１実際値を、Ｍ個の領域についての前記第１実際値と前記第２実際値との関係に基づいて前記所定高度で観測が行われたとした場合の値となるように補正して補正値を導出するステップと、
   当該一の前記領域の前記農作物関連値を導出する際に、第１指標の指標値を入力とし前記農作物関連値を出力とする計算式であって、前記所定高度で観測が行われた場合の第１指標の指標値を入力したときに適切な値の前記農作物関連値を出力するように調整された計算式である第１計算式に、当該一の前記領域について導出した補正値を入力することによって、前記農作物関連値を導出するステップとを含む
   ことを特徴とする農作物関連値導出方法。
   

Images