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