(57)【要約】【課題】リモートセンシング画像、作物生理、環境データ及び気象を用いて、時系列分析で植物の常時観測、収穫管理などを実現する。【解決手段】プロセッサと、プロセッサがアクセスする記憶装置と、を有する作物管理システムであって、記憶装置は、農場の複数の区画を過去のそれぞれ異なる時期に撮影した複数の画像を含む画像データと、過去のそれぞれ異なる時期に計測された区画ごとの気象データと、を保持し、プロセッサは、気象データに基づいて、作物の生長状況を示す指標を推定するための生長モデルを生成し、画像データに基づいて、区画ごとに、生長モデルを補正する。【選択図】図１
2020140347.tif 000002
【特許請求の範囲】
【請求項１】
プロセッサと、前記プロセッサがアクセスする記憶装置と、を有する作物管理システムであって、
前記記憶装置は、
農場の複数の区画を過去のそれぞれ異なる時期に撮影した複数の画像を含む画像データと、
過去のそれぞれ異なる時期に計測された前記区画ごとの気象データと、を保持し、
前記プロセッサは、
前記気象データに基づいて、作物の生長状況を示す指標を推定するための生長モデルを生成し、
前記画像データに基づいて、前記区画ごとに、前記生長モデルを補正することを特徴とする作物管理システム。
【請求項２】
請求項１に記載の作物管理システムであって、
前記プロセッサは、
前記気象データに基づいて、時間ごとの前記作物の生長状況を示す指標を計算し、時間ごとの気温を取得し、
前記時間ごとの前記作物の生長状況を示す指標及び前記時間ごとの気温に基づいて、積算温度に対する前記作物の生長状況を示す指標を、前記生長モデルとして生成し、
前記積算温度に基づいて、前記作物の生長段階の区分を抽出することを特徴とする作物管理システム。
【請求項３】
請求項２に記載の作物管理システムであって、
前記プロセッサは、前記生長段階の区分ごとに、前記積算温度と前記作物の生長状況を示す指標との関係を示す関数を生成して前記関数のパラメータを前記生長モデルとして前記記憶装置に保持することを特徴とする作物管理システム。
【請求項４】
請求項２に記載の作物管理システムであって、
前記プロセッサは、
前記画像データに基づいて、前記区画ごとに、前記作物の量及び活力を示す植生指標を計算し、
前記植生指標に基づいて、前記作物の生長状況を示す指標を計算し、
前記植生指標に基づいて計算した前記作物の生長状況を示す指標を用いて、前記生長モデルを補正することを特徴とする作物管理システム。
【請求項５】
請求項４に記載の作物管理システムであって、
前記作物の生長状況を示す指標は、葉面積指数、バイオマス又は糖度のいずれかであり、
前記植生指標は、可視領域の赤の反射率と近赤外領域の反射率とに基づく正規化差植生指数であり、
前記積算温度は、前記作物の種類に応じた基準温度に基づく有効積算温度であることを特徴とする作物管理システム。
【請求項６】
請求項１に記載の作物管理システムであって、
前記プロセッサは、補正された前記生長モデルに基づいて、生長段階が統一された、前記複数の区画における前記作物の生長状況を示す画像を生成することを特徴とする作物管理システム。
【請求項７】
請求項６に記載の作物管理システムであって、
前記プロセッサは、
前記複数の区画のうち第１の区画を基準区画として選択し、
第１の時点の、前記複数の区画のうち前記第１の区画とは異なる第２の区画における積算温度に基づく生長段階と、前記第１の時点の、前記第１の区画における積算温度に基づく生長段階との相違を計算し、
前記第２の区画の生長モデルと、計算された前記生長段階の相違と、に基づいて、前記第２の区画における前記作物の生長状況を示す指標を、前記第２の区画における前記生長段階が前記第１の区画における前記生長段階と同一となる時点の、前記第２の区画における前記作物の生長状況を示す値に補正し、
前記第１の時点の、前記第１の区画における前記作物の生長状況を示す指標の値と、補正された前記第２の区画における前記作物の生長状況を示す指標の値と、に基づいて、前記複数の区画における前記作物の生長状況を示す画像を生成することを特徴とする作物管理システム。
【請求項８】
請求項７に記載の作物管理システムであって、
前記プロセッサは、
前記生長段階の区分ごとに、前記積算温度と前記作物の生長状況を示す指標との関係を示す関数を生成して前記関数のパラメータを前記生長モデルとして前記記憶装置に保持し、
前記第１の区画における前記生長段階を含む前記区分から前記第２の区画における前記生長段階を含む前記区分までの全ての前記区分の前記関数のパラメータを使用して、前記第２の区画における前記作物の生長状況を示す指標の補正を行うことを特徴とする作物管理システム。
【請求項９】
請求項７に記載の作物管理システムであって、
前記プロセッサは、他の区画の前記積算温度との関係が所定の条件を満たす積算温度を有する区画を前記基準区画として選択することを特徴とする作物管理システム。
【請求項１０】
請求項７に記載の作物管理システムであって、
前記プロセッサは、前記第２の区画に対応するポリゴン内の複数のピクセルに同一の補正を行うことによって、前記複数の区画における前記作物の生長状況を示す画像を生成することを特徴とする作物管理システム。
【請求項１１】
請求項１に記載の作物管理システムであって、
前記プロセッサは、補正された前記生長モデルに基づいて、前記複数の区画における前記作物の収穫順序を示す画像を生成することを特徴とする作物管理システム。
【請求項１２】
プロセッサと、前記プロセッサがアクセスする記憶装置と、を有する計算機システムが実行する作物管理方法であって、
前記記憶装置は、
農場の複数の区画を過去のそれぞれ異なる時期に撮影した複数の画像を含む画像データと、
過去のそれぞれ異なる時期に計測された前記区画ごとの気象データと、を保持し、
前記作物管理方法は、
前記プロセッサが、前記気象データに基づいて、作物の生長状況を示す指標を推定するための生長モデルを生成する手順と、
前記プロセッサが、前記画像データに基づいて、前記区画ごとに、前記生長モデルを補正する手順と、を含むことを特徴とする作物管理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、作物の生長状況及び収穫等を管理するためのシステム及び方法に関する。
【背景技術】
【０００２】
大規模な農場においては、現場で詳細を把握することは難しく、ＧＩＳ（Ｇｅｏｇｒａｐｈｉｃａｌ  Ｉｎｆｏｒｍａｔｉｏｎ  Ｓｙｓｔｅｍ）や衛星を利用した管理手法が求められている。常時監測、コスト削減及び早期収量予測を実現するため、衛星画像又はドローンが撮影した画像などを利用する分析手法が使われている。例えば、特開２０１６−１２３３６９号公報（特許文献１）には、「リモートセンシング画像に基づいて、植物の成長を分析する植物成長分析システムであって、植物の所定時期からの経過時間を示す時間情報に対応して、植物の成長の特徴量の変化が登録された植物成長モデルと、リモートセンシング画像の一部である植物の複数の生育場所に対応する画像に基づいて、特徴量を算出する特徴量算出部と、植物成長モデルを参照し、複数の生育場所の特徴量を基準時間情報に対応する特徴量に補正し、複数の生育場所の生育差を補正する生育差補正部と、リモートセンシング画像の一部である複数の生育場所に対応する画像が補正後の特徴量となるように、新たな画像を生成する画像生成部と、を備えることを特徴とする。」と記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−１２３３６９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＧＩＳデータに基づき、圃場を地図に表示し、圃場データをマッピングすることによって、農場の管理を行う場合、作物の植え付け日の差に起因して、生長段階にばらつきがあり、そのために分析及び予測などの観測管理が困難になっている。
【０００５】
一方、現状では、早期収量予測を行い、収穫プランニングを作成しているが、収穫時間については農家の経験に依存する。作物には、最適な収穫時点があり、その時点をこえると収量が下がるのが普通である。また、収穫順序のプランニングによって、収穫機器の待機消耗が大きいことがコストとなる。
【０００６】
上記の特許文献１には、生育段階にズレがある植物の成長の分析において、リモートセンシング画像から算出した特徴量を補正することが記載されているが、作物の生長状況を推定するための生長モデルを補正することは記載していない。
【課題を解決するための手段】
【０００７】
本願において開示される発明のうち、代表的なものの概要を簡単に説明すれば、下記のとおりである。すなわち、プロセッサと、前記プロセッサがアクセスする記憶装置と、を有する作物管理システムであって、前記記憶装置は、農場の複数の区画を過去のそれぞれ異なる時期に撮影した複数の画像を含む画像データと、過去のそれぞれ異なる時期に計測された前記区画ごとの気象データと、を保持し、前記プロセッサは、前記気象データに基づいて、作物の生長状況を示す指標を推定するための生長モデルを生成し、前記画像データに基づいて、前記区画ごとに、前記生長モデルを補正することを特徴とする。
【発明の効果】
【０００８】
本発明の一態様によれば、生長ステージ統一画像を生成し、作物の生長段階を補正でき、画像内の作物の特徴量のバラツキが生長段階ズレの影響を受けなくなり、解析の精度が向上でき、管理の精度も向上する。
【０００９】
上記した以外の課題、構成及び効果は、以下の実施例の説明によって明らかにされる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施例における作物管理収穫システムの構成を表すブロック図である。
【図２】本発明の実施例における作物管理収穫システムの全体的なハードウェア構成を示すブロック図である。
【図３】本発明の実施例における作物管理収穫システム全体の処理例を説明するフローチャートである。
【図４】本発明の実施例における生長モデル生成の処理を表す説明図である。
【図５】本発明の実施例における生長ステージ抽出の処理を表す説明図である。
【図６】本発明の実施例における生長ステージ解析の処理を表す説明図である。
【図７】本発明の実施例における生長モデル補正の処理を表す説明図である。
【図８】本発明の実施例における生長差計算の処理を表すフローチャートである。
【図９】本発明の実施例における生長差補正の処理の第１の例を表す説明図である。
【図１０】本発明の実施例における生長差補正の処理の第２の例を表す説明図である。
【図１１】本発明の実施例における生長差補正の処理の第３の例を表す説明図である。
【図１２】本発明の実施例における生長差補正の処理の第４の例を表す説明図である。
【図１３】本発明の実施例における生長差補正の処理の各ケースを表す説明図である。
【図１４Ａ】本発明の実施例におけるデータベースの例を表す説明図である。
【図１４Ｂ】本発明の実施例におけるデータベースに格納される時系列衛星画像データの例を表す説明図である。
【図１４Ｃ】本発明の実施例におけるデータベースに格納されるＧＩＳデータの例を表す説明図である。
【図１４Ｄ】本発明の実施例におけるデータベースに格納される植物生長データの例を表す説明図である。
【図１４Ｅ】本発明の実施例におけるデータベースに格納される気象データの例を表す説明図である。
【図１４Ｆ】本発明の実施例におけるデータベースに格納される土壌及び他の相関データの例を表す説明図である。
【図１４Ｇ】本発明の実施例におけるデータベースに格納される処理パラメータの例を表す説明図である。
【図１５】本発明の実施例における収穫マップの処理の各ケースを表す説明図である。
【発明を実施するための形態】
【００１１】
以下に本発明の作物管理システム及び方法を適用した実施例を説明する。
【００１２】
本実施例の作物管理システムは、作物の育成及び収穫等を管理するシステムであり、時系列リモートセンシング画像データの対象地域及び圃場の生長ステージを補正した新たな作物生理特徴量マップを出力する。図１３及び図１５は本実施例における入力と出力を表す模式図である。本実施例では、解析単位は圃場とするが、他の単位で解析を行ってもよい。
【００１３】
なお、圃場とは、農場の区画の一例である。例えば、各圃場は道路又は水路等によって区切られていててもよいが、圃場の大きさ、形状及び境界線の決め方は限定されない。以下の説明における圃場は、任意の区画に置き換えることができる。
【００１４】
図１は、本発明の実施例における作物管理システム１００の構成を表すブロック図である。
【００１５】
作物管理システム１００は、入力データ１０１、作物生理・環境データ取得部１０２、気象・ＧＩＳ（Ｇｅｏｇｒａｐｈｉｃａｌ  Ｉｎｆｏｒｍａｔｉｏｎ  Ｓｙｓｔｅｍ）データ取得部１０３、リモートセンシング画像データ取得部１０４、生長モデル生成部１０５、生長ステージ抽出部１０６、生長ステージ解析部１０７、生長モデル補正部１０８、生長差計算部１０９及び生長差補正部１１０を備える。入力データ１０１は、作物生理データ及び環境データ（以下、作物生理・環境データ１０１ａと記載する）、気象データ及びＧＩＳデータ（以下、気象・ＧＩＳデータ１０１ｂと記載する）、並びに、時系列リモートセンシング画像データ１０１ｃを含む。
【００１６】
作物管理システム１００は、作物生長時期の時系列リモートセンシングデータ、気象・ＧＩＳデータ、作物生理・環境データを取得し、データベース１１５に提供する。リモートセンシング画像データ１０１ｃ、気象・ＧＩＳデータ１０１ｂ及び作物生理・環境データ１０１ａは図１で説明する。また、データベース１１５は図１４で説明する。
【００１７】
データベース１１５が管理するデータは、生長ステージ統一画像生成部１１１及び作物収穫画像生成部１１２に提供される。生長ステージ統一画像生成部１１１及び作物収穫画像生成部１１２は、提供された情報を使用して、推計対象（圃場又は地域）の生長段階のズレを補正した新たな分析マップ（生長ズレ補正マップ）と、作物収穫画像とを出力する。
【００１８】
生長ステージ統一画像生成部１１１で生成した新たな生長ズレ補正マップは、生長分析システム使用部に提供され、作物の管理に使用される。作物収穫画像生成部１１２が生成した作物収穫画像は、物流配車システム使用部に提供され、作物の収穫プランに使用される。
【００１９】
図２は、本発明の実施例における作物管理システム１００の全体的なハードウェア構成を示すブロック図である。
【００２０】
リモートセンシング観測装置１２０から、時系列リモートセンシング画像を撮影する様子を示す。リモートセンシング観測装置１２０の例としては、観測衛星又は航空機などがあり、種類としては限定してない。本実施例は観測衛星が取得した衛星画像を利用する例について説明する。図２は、作物管理システム１００がコンピュータシステム１２１に実装されている例を示す。すなわち、図１の作物管理システム１００の主要部は、コンピュータシステム１２１で処理されている。このコンピュータシステムは、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１２２、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）１２３及び記憶部１２４などの一般的なハードウェア構成を有する。
【００２１】
図１に戻り説明を続ける。入力データ１０１（作物生理・環境データ１０１ａ、気象・ＧＩＳデータ１０１ｂ及び時系列リモートセンシング画像データ１０１ｃ）はリモートセンシングデータであり、ＲａｐｉｄＥｙｅ、ＭＯＤＩＳもしくはＬａｎｄＳａｔなどのような人工衛星からのデータ又は航空写真等である。なお、本実施例では時系列リモートセンシング画像を示したが、データの種類は限定せず、入力データ１０１が他の種類のデータであってもよい。また、同一のリモートセンシング画像における時系列画像を用いてもよい。なお、解像度の違いは相対比較である。
【００２２】
気象データとしては、アメダスのデータ、ＭＯＤＩＳ表面温度衛星からのデータ、又は気象局データなど様々ありうるが、種類としては限定してない。ＧＩＳデータとは、圃場ポリゴン又はＧＰＳデータなど、現地データと衛星データとを共用する際に不可欠なデータである。作物生理データとは、葉面積、茎重さ及びバイオマスなどの生理特性データであるが、種類としては限定してない。環境データとは、土壌データ及び生態データなどの作物生長の環境に関するデータであるが、種類としては限定してない。データは前処理を行った後にデータベース１１５に格納される。データベースの詳細は図１４で説明する。
【００２３】
作物生理・環境データ取得部１０２、気象・ＧＩＳデータ取得部１０３、リモートセンシング画像データ取得部１０４、生長モデル生成部１０５、生長ステージ抽出部１０６、生長ステージ解析部１０７、生長モデル補正部１０８、生長差計算部１０９、生長差補正部１１０、生長ステージ統一画像生成部１１１、作物収穫画像生成部１１２、生長分析システム使用部１１３及び物流配車システム使用部１１４は、各演算処理を行うため、役割毎に分けられた複数のＣＰＵ及びＲＡＭの組み合わせとして実装される。前記の各部は外部記憶装置としてのハードディスク、ＵＳＢメモリなどを採用する。
【００２４】
例えば、図１に示した作物生理・環境データ取得部１０２、気象・ＧＩＳデータ取得部１０３、リモートセンシング画像データ取得部１０４、生長モデル生成部１０５、生長ステージ抽出部１０６、生長ステージ解析部１０７、生長モデル補正部１０８、生長差計算部１０９、生長差補正部１１０、生長ステージ統一画像生成部１１１、作物収穫画像生成部１１２、生長分析システム使用部１１３（後述）及び物流配車システム使用部１１４（後述）は、それぞれ、ＲＡＭ１２３に格納されたプログラムによって実現されてもよい。例えば、それらのプログラムは、記憶部１２４（例えばハードディスク装置又はフラッシュメモリ等）に格納され、それらの少なくとも一部が必要に応じて読み出され、ＲＡＭ１２３に格納されてもよい。すなわち、以下の説明において上記の各部が実行する処理は、実際には、ＣＰＵ１２２（すなわちプロセッサ）がＲＡＭ１２３に格納されたプログラムに従って実行する。
【００２５】
同様に、データベース１１５は、記憶部１２４に格納され、その少なくとも一部が必要に応じて読み出され、ＲＡＭ１２３に格納されてもよい。
【００２６】
図３は、本発明の実施例における作物管理システム１００全体の処理例を説明するフローチャートである。
【００２７】
処理が開始されると（Ｓ３０１）、まず、Ｓ３０２において、時系列リモートセンシング画像データ１０１ｃがリモートセンシング画像データ取得部１０４に、気象・ＧＩＳデータ１０１ｂが気象・ＧＩＳデータ取得部１０３に、作物生理・環境データが作物生理・環境データ取得部１０２に、それぞれ入力される。本実施例では、気象データは圃場の表面温度のデータであるが、他の方法で計測した温度データでもよい。
【００２８】
時系列リモートセンシング画像データ、温度データ、作物生理データ及び環境データは、ＧＩＳデータと対応付けられて、データベース１１５に提供される。ＧＩＳデータには、位置の情報と形状の情報が含まれ、位置と形状によって、リモートセンシングデータに対して圃場のマスキングを行う。マスキング圃場が抽出され、その圃場のリモートセンシングデータと表面温度とがデータベース生成部から生成され、データベース１１５にデータが提供される。
【００２９】
Ｓ３０３では、生長モデル生成部１０５が、データベース１１５に格納されたＧＩＳデータの位置情報及び形状情報からマスキング処理を行い、その結果の圃場の作物生理データ、環境データ及び表面温度データから、各圃場の時間ごとの作物生理特徴量及び有効積算温度データを算出してデータベース１１５に格納する。
【００３０】
Ｓ３０４では、データベース１１５に格納された作物生理特徴量及び有効積算温度データを生長モデル生成部１０５に提供し、生長モデル生成部１０５がそれらを作物特徴量及び有効積算温度の座標系の曲線（例えば、横軸を有効積算温度、縦軸を作物特徴量とする曲線）に変換する。Ｓ３０４での作物特徴量及び有効積算温度データ算出とＳ３０５での生長曲線生成の詳細については、後に図４を参照して説明する。
【００３１】
Ｓ３０５では、生長モデル生成部１０５によって格納された作物特徴量及び有効積算温度の推移モデルを、生長ステージ抽出部１０６に提供し、生長ステージ抽出部１０６が作物の栽培特性から作物の生長ステージを抽出する。Ｓ３０５での作物生長ステージ抽出の詳細については、後に図５を参照して説明する。
【００３２】
Ｓ３０６では、生長ステージ抽出部１０６に格納された生長ステージ抽出済みの作物特徴量、有効積算温度の推移モデルを、生長ステージ解析部１０７に提供し、生長ステージ解析部１０７が作物の各ステージの生物特性に応じて、各ステージの生長モデルを解析する。Ｓ３０６での作物ステージの生物特性に応じた各ステージの生長モデルのパラメータ化の詳細については、後に図６を参照して説明する。
【００３３】
そして、Ｓ３０７では、生長モデル補正部１０８が、生長モデル生成部１０５によって生成された作物の生長曲線モデルに対して、生長ステージ抽出部１０６が抽出した生長ステージを生長ステージ解析部１０７が解析した結果に基づいて計算したパラメータを用いて、時系列リモートセンシング画像データから抽出したスペクトル特徴量／有効積算温度の分布を補正する。Ｓ３０７での生長モデル補正の詳細については、後に図７を参照して説明する。
【００３４】
Ｓ３０８では、生長差計算部１０９が、高解像度リモートセンシング画像データに対して、ＧＩＳ圃場マスキング及び作物生理特徴量を計算し、データベース１１５で格納した作物生理特徴量の各圃場の中に、ある圃場を基準とし、他の圃場と基準圃場の生長段階の差を計算する。Ｓ３０８での生長差計算の詳細については、後に図８参照して説明する。
【００３５】
Ｓ３０９では、生長差補正部１１０が、生長ステージ解析部１０７及び生長モデル補正部１０８によって計算された各圃場のパラメータと、生長差計算部１０９によって計算された各圃場の基準圃場に対する生長差とを用いて、作物生理特徴量を補正する。Ｓ３０９での生長差補正の詳細については、後に図９、図１０、図１１及び図１２を参照して説明する。
【００３６】
Ｓ３１０では、生長ステージ統一画像生成部１１１が、生長差補正部１１０によって計算された各圃場に対する補正した作物生理特徴量を有する圃場について、生長ズレ補正マップを生成する。Ｓ３１０での生長ステージ統一画像生成の詳細については、後に図１３を参照して説明する。
【００３７】
Ｓ３１１では、生長分析システム使用部１１３が、Ｓ３１０で生成された生長ステージ統一画像を用いて、作物生長管理を行う。
【００３８】
Ｓ３１２では、作物収穫画像生成部１１２が、生長モデル補正部１０８によって補正された作物生長モデルを用いて、作物収穫画像を生成する。Ｓ３１２での作物収穫画像生成の詳細については、後に図１５を参照して説明する。
【００３９】
Ｓ３１３では、物流配車システム使用部１１４が、Ｓ３１２で生成された作物収穫画像を用いて、物流配車最適化を行う。
【００４０】
すべてのデータはデータベース１１５に格納される。データベース１１５の一例は図１４に示す。
【００４１】
以下、上記の処理の詳細を説明する。
【００４２】
作物生理・環境データ取得部１０２は、作物生理・環境データ１０１ａを取得し、データベース１１５に提供する。
【００４３】
気象・ＧＩＳデータ取得部１０３は、気象・ＧＩＳデータ１０１ｂを取得し、データベース１１５に提供する。
【００４４】
リモートセンシング画像データ取得部１０４は、リモートセンシング画像データ１０１ｃを取得し、データベース１１５に提供する。
【００４５】
続いて、生長モデル生成部１０５は、データベース１１５に格納された圃場の作物生理特徴量データと表面温度データ（気象・ＧＩＳデータ１０１ｂに含まれる）とを用いて、各分析対象（圃場、地域など）の作物生理特徴量、有効積算温度生長曲線を作成する。
【００４６】
作物生理特徴量は、作物指標などの作物の生長状態を示すパラメータである。本実施例では、作物生理指標を用いて説明する。作物生理指標とは、農業分野で、地面センサーから取ったデータに基づいて、様々の作物生長状況を反映する指標である。最も使われている作物指標の一つは、ＬＡＩ（葉面積指数）、バイオマス、又は糖度などであり、本実施例では、ＬＡＩを例として説明する。ただし、他の作物指標を使用してもよい。生長モデル生成部１０５は、数式（１）を用いて有効積算温度（ＥＡＴ）を算出する。
【００４７】
【数１】
2020140347.tif 000003
【００４８】
ここで、ｔｄは一日の平均温度である。ｔ０は有効積算温度の基準温度であり、作物によって異なる。本実施例ではｔ０は１５度に設定される。ｉは日数である。
【００４９】
ＥＡＴは、ある期間の日々の平均気温のうち、一定の基準値を超えた分を取り出して合計したものである。ＥＡＴは植物の生長に必要な熱量の大小の目安に用いられ、基準値として、たとえば春コムギで３℃以上、トウモロコシで１３℃以上といった値が用いられる。
【００５０】
本実施例では、サトウキビを例として説明する。ほかの品種及び植物でもよい。
【００５１】
作物指標／有効積算温度生長モデルは、作物の特性として重要な時期の日付から、有効積算温度を積算し、その値を座標系の横軸に設定することで生成される。
【００５２】
そして、生長モデル生成部１０５は、各時期の作物指標をデータベース１１５から取得し、その値を座標系の縦軸に設定する。植物の植え付け日の日付がある場合は、植え付け日からの曲線が生成される。ただし、圃場地域によっては、植え付け日が分からない場合もあるので、開花期などのほか重要時期の日付から計算してもよい。データがある時期の点を座標系の中に入れ、全ての点に基づいて曲線を生成する。曲線の生成手法の例の一つは、二次近似であるが、スプライン補間など、他の手法でもよい。
【００５３】
図４は、本発明の実施例における生長モデル生成の処理を表す説明図である。
【００５４】
生長モデル生成部１０５は、時間ごとの生理データと、時間ごとの気温データと、に基づいて、積算温度に対する生理データを生成する（図４（ａ））。具体的には、生長モデル生成部１０５は、植物の重要時期の日からの積算温度を計算し、それを横軸に設定し、それに対応する作物指標ＬＡＩを縦軸に設定して、ＬＡＩ／積算温度分析曲線を生成する（図４（ｂ））。重要時期とは、作物の生長フェーズが変わる時期である。
【００５５】
本実施例のサトウキビの場合では、Ｒ１、Ｒ２、Ｒ３及びＲ４いう生長時期がある。Ｒ１とは、最初の花が出る時期である。Ｒ２とは、茎が生長する時期である。Ｒ３とは、糖分を蓄積する時期である。Ｒ４とは、葉が５０％以上黄色になり、糖分が減り、収穫される時期である。それぞれの重要時期では、この時期になるために必要なＥＡＴがあり、そのＥＡＴの時期のＬＡＩ値を算出できる。
【００５６】
生長モデルとするＬＡＩは、公知である植生生長モデルから作成できる。例えばＣＡＮＥＧＲＯモデルにさまざまな生理パラメータ及び環境パラメータを入力すれば、サトウキビの各生理特徴量に対して、植え付けから収穫までの一サイクルの生長曲線を推定できる。
【００５７】
続いて、生長ステージ抽出部１０６は、生長モデル生成部１０５で生成したある対応地域、圃場の生長モデルに基づき、植物の栽培生長特性を考量し、生長ステージを抽出する。
【００５８】
図５は、本発明の実施例における生長ステージ抽出の処理を表す説明図である。
【００５９】
生長ステージは、作物の栽培生育特性を考慮した上で抽出される。例えばサトウキビの場合は、Ｒ１、Ｒ２、Ｒ３及びＲ４という生長時期があり、Ｒ１−Ｒ２、Ｒ２−Ｒ３、及びＲ３−Ｒ４の三つの生長ステージに分けられる。Ｒ１−Ｒ２とは、花が出る時期から花が完全に成長し終わるまでの開花期である。Ｒ２−Ｒ３とは、サトウキビの茎が最大に生長するまでの時期である。Ｒ３−Ｒ４とは、体の生長が止まり、殆どの糖分が蓄積される時期である。Ｒ４以後とは、葉が５０％以上黄色になり、糖分が減り、収穫される時期である。
【００６０】
データは農家又は作物の基本生長データから入手し、データベース１１５に格納する。作物の基本生長データとしては、既に説明したように、各時期の日数及び必要な積算温度に関する公知のデータがある場合がある。
【００６１】
ほかの植物を対象とする場合は、作物の基本生長データが変わるが、抽出原理は変わらない。その作物の生長フェーズを調べた上で、各重要時期の必要積算温度を計算し、その積算温度のＬＡＩを算出することができる。
【００６２】
生長ステージ解析部１０７は、生長ステージ抽出部で抽出した作物の栽培生育時期に基づき、各ステージの生物性を分析し、各ステージの作物生長モデルを分解する。
【００６３】
図６は、本発明の実施例における生長ステージ解析の処理を表す説明図である。
【００６４】
例えばサトウキビの場合は、Ｒ１、Ｒ２、Ｒ３、Ｒ４という生長時期（言い換えると生長段階の区切り）があり、Ｒ１−Ｒ２、Ｒ２−Ｒ３、Ｒ３−Ｒ４の三つの生長ステージ（すなわち生長段階の区分）に分けられる。
【００６５】
生長モデル補正部１０８は、生長モデル生成部１０５から生成した作物の生長曲線モデルに対して、生長ステージ抽出部１０６で抽出した生長ステージを生長ステージ解析部１０７の解析に基づいて計算したパラメータを用いて、リモートセンシング画像データから抽出したＮＤＶＩ（正規化差植生指数）／有効積算温度の分布（図６（ａ））を補正して、ＬＡＩ＿ＮＤＶＩ／有効積算温度の分布（図６（ｂ））を生成する。ＮＤＶＩは、下記の数式（２）によって計算される。
【００６６】
【数２】
2020140347.tif 000004
【００６７】
ここで、Ｒは、可視領域の赤の反射率であり、ＩＲは、近赤外領域の反射率である。これらの領域ごとの反射率は衛星画像から抽出できる。
【００６８】
図７は、本発明の実施例における生長モデル補正の処理を表す説明図である。
【００６９】
生長モデル補正部１０８は、リモートセンシングのデータを取得し、作物生理データから得られた生長モデルを補正するため、スペクトル特徴量と作物生理特徴量との数学的関係性を推定する。そして、生長モデル補正部１０８は、観測されたリモートセンシングデータから得られたスペクトル特徴量曲線を、作物生理特徴量曲線に変換する。変換された作物生理特徴量曲線と元の作物生長モデルの曲線とをフィッティングし、元の作物生長モデルを補正する。
【００７０】
例えば、ある年（当該年度）の、作物が生長の途上にある時点に、その作物の生長状況を知るために生長モデルを利用しようとする場合に、生長モデル補正部１０８は、前年度までの気象データ等に基づいて生成された生長モデル（図４、図５）を、当該年度の現時点までのリモートセンシング画像データから得られたＬＡＩ＿ＮＤＶＩ（図６）によって補正してもよい。
【００７１】
補正の例として、サトウキビの場合では生長ステージがＲ１−Ｒ２、Ｒ２−Ｒ３、Ｒ３−Ｒ４の三つに分けられ、各ステージにおけるスペクトルモデルと作物生理モデルとの数学的関係式が異なる。生長モデル補正部１０８は、各ステージのＮＤＶＩ（図７（ａ））からＬＡＩ＿ＮＤＶＩ（図７（ｂ））を計算し、その結果をＬＡＩとフィッティングすることで、ＬＡＩ曲線のパラメータを最適化する（図７（ｃ））。ＬＡＩ曲線のパラメータ最適化によって、作物生理特徴量曲線は補正され、収量、糖度などの曲線も推定できるようになる。他の植物または他の作物特徴量の場合には、各ステージの生物性及び生物性から作物特徴量に反映する関係が異なるため、それらを用いることで、上記の解析手段をサトウキビに限らずに適用することができる。
【００７２】
生長差計算部１０９は、作物生理データ及び環境データに対して、ＧＩＳ圃場マスキング及び作物生理特徴量を計算し、データベース１１５に格納されたリモートセンシング画像データの各圃場の中で、ある圃場を基準とし、他の圃場と基準圃場の生長段階の差を計算する。
【００７３】
図８は、本発明の実施例における生長差計算の処理を表すフローチャートである。
【００７４】
計算の方法は次のとおりである。データベース１１５に格納された各圃場の植え付け時期又はある重要時期を用いて、基準圃場と他の圃場の植え付け時期又はある重要時期との積算温度の差を計算する。この差を各圃場の補正基準の値としてデータベース１１５に提供する。基準圃場の選び方として、二つの方法がある。リモートセンシング画像からある特定の時期の画像を生成したい場合には、この時期の圃場を基準圃場とする。そうでない場合には、各圃場の植え付け時期又はある重要時期との積算温度に基づき、積算温度の差の合計が最小になる圃場を基準としてもよい（例えば数式（３））。あるいはさらに他の方法で基準圃場を選択してもよい。
【００７５】
【数３】
2020140347.tif 000005
【００７６】
生長差補正処理Ｓ１０６０は、各圃場の生長差を計算する処理である。Ｓ１０６１において、生長差計算部１０９は、各圃場の積算温度値を計算する。Ｓ１０６２において、生長差補正部１１０は、基準圃場を前述の方法で選択する。Ｓ１０６３において、生長差補正部１１０は、各圃場と基準圃場の積算温度差を計算する。Ｓ１０６３において、生長差補正部１１０は、各圃場の積算温度と基準圃場の積算温度との差を生長差としてデータベース１１５に格納する。
【００７７】
生長差補正部１１０は、生長ステージ解析部１０７と生長モデル補正部１０８が計算した各圃場のパラメータと、生長差計算部１０９が計算した各圃場の基準圃場にたいする生長差とを用いて、各圃場の作物特徴量を補正する。
【００７８】
図９、図１０、図１１及び図１２は、それぞれ、本発明の実施例における生長差補正の処理の第１から第４の例を表す説明図である。
【００７９】
本実施例の補正では、サトウキビの生長ステージがＲ１−Ｒ２、Ｒ２−Ｒ３、及びＲ３−Ｒ４の三つに分けられ、各ステージに対して計算及び補正する生長曲線がある。図９から図１２に示す生長曲線はいずれも同じものであり、図７を参照して説明した方法で補正されたＬＡＩを示す。
【００８０】
図９から図１２のいずれの例においても、基準圃場はいずれもＲ２−Ｒ３の生長ステージにある（すなわち基準圃場の積算温度がＲ２からＲ３の範囲内にある）。一方、補正対象である他の圃場（以下、該当圃場とも記載する）は、図９の例においては、Ｒ２−Ｒ３の生長ステージであって、基準圃場の左側にあり、図１０の例においては、Ｒ２−Ｒ３の生長ステージであって、基準圃場の右側にあり、図１１の例においては、Ｒ３−Ｒ４のステージにあり、図１２の例においては、Ｒ１−Ｒ２のステージにある。図９及び図１２の例では、該当圃場の積算温度が基準圃場の積算温度より小さい（すなわち該当圃場の生長段階が基準圃場の生長段階より遅い）。一方、図１０及び図１１の例では、該当圃場の積算温度が基準圃場の積算温度より大きい（すなわち該当圃場の生長段階が基準圃場の生長段階より早い）。
【００８１】
図９の例では、生長差補正部１１０は、該当圃場の作物特徴量を、積算温度差とＲ２−Ｒ３ステージの生長モデル関数のパラメータとに応じて横軸のプラス方向（すなわち基準圃場の積算温度の方向）に補正する計算をする。
【００８２】
図９において、黒丸の点は、該当圃場の補正された生長モデルにおける、該当圃場の積算温度に対応するＬＡＩを示し、白丸は、基準圃場の積算温度に対応するＬＡＩを示す。黒丸の点のＬＡＩが、白丸の点のＬＡＩに補正される。両者の差が補正量である。図１０〜図１２でも同様である。
【００８３】
図１０の例では、生長差補正部１１０は、該当圃場の作物特徴量を、積算温度差とＲ２−Ｒ３ステージの生長モデル関数のパラメータとに応じて横軸のマイナス方向（すなわち基準圃場の積算温度の方向）に補正する計算をする。
【００８４】
図１１の例では、生長差補正部１１０は、該当圃場の積算温度と基準圃場の積算温度との差をＲ２−Ｒ３ステージ部分とＲ３−Ｒ４ステージ部分とに分ける。そして、生長差補正部１１０は、該当圃場の作物特徴量を、Ｒ３−Ｒ４ステージの生長モデル関数のパラメータを用いてＲ６時点に補正し、続いてＲ２−Ｒ３ステージの生長モデル関数のパラメータを用いて基準圃場の積算温度の方向に補正する計算をする。
【００８５】
図１２の例では、生長差補正部１１０は、該当圃場の積算温度と基準圃場の積算温度との差をＲ２−Ｒ３ステージ部分とＲ１−Ｒ２ステージ部分とに分ける。そして、生長差補正部１１０は、該当圃場の作物特徴量を、Ｒ１−Ｒ２ステージの生長モデル関数のパラメータを用いてＲ２時点に補正し、続いてＲ２−Ｒ３ステージの生長モデル関数のパラメータを用いて基準圃場の積算温度の方向に補正する計算をする。
【００８６】
基準圃場がＲ１−Ｒ２ステージ又はＲ３−Ｒ４ステージにある場合は、生長差補正部１１０は、上記の処理と同じように、基準圃場に向けて、補正対象の圃場のあるステージに応じて、補正値を計算し、補正圃場の作物特徴量を計算する。
【００８７】
図１１のケースを例として、ステージの生長モデル関数のパラメータと積算温度差とから補正する方法を数式（４）に示す。他の植物または他の作物特徴量の場合には、各ステージの生物性及び生物性から作物特徴量に反映する関係が異なるため、線形関数ではない場合があり、その関数パラメータに応じて補正を行う。
【００８８】
【数４】
2020140347.tif 000006
【００８９】
ここで、
（ｘ１，ｙ１）は補正前のデータ
（ｘ２，ｙ２）は補正後のデータ
ｘ２は基準圃場の積算温度値（既知）
ｙ２は算出対象
（ｘＲ３，ｙＲ３）はＲ３時期のデータ
ｘＲ３は既知
ｙＲ３は算出対象
である。
【００９０】
生長ステージ統一画像生成部１１１は、生長差補正部１１０が計算した各圃場の補正した作物生理特徴量を有する圃場によって、リモートセンシング画像の各対応圃場を置き替え、生長ズレ補正のマップを生成する。
【００９１】
例えば、生長ステージ統一画像生成部１１１は、基準圃場以外の各圃場の作物生理特徴量の補正量（具体的には、例えば、図９〜図１２の例における黒丸の点と白丸の点のＬＡＩの値の差）に基づいて、各圃場のポリゴン内の画素（ピクセル）の値の変更量を決定し、その変更量に従って各画素の値を変更してもよい。補正値は圃場ごとに計算されるため、圃場に対応するポリゴン内の画素に同一の変更量が適用される。これによって、圃場内の作物生理特徴量のばらつきに対応するポリゴン内の画素の値のばらつきは維持したまま、ポリゴン全体の画素の値のレベルを変更することができる。
【００９２】
図１３は、本発明の実施例における生長差補正の処理の各ケースを表す説明図である。
【００９３】
図１３（ａ）に示したように、元の画像は時間的に統一された画像である。一方、図１３（ｂ）に示した処理後の画像は、時間にズレがあるが、生長段階が統一された画像になる。このような画像は、例えば、各圃場の補正値を圃場内のピクセルにかけて、生長差を補正することによって生成してもよい。
【００９４】
生長ステージ統一画像生成部１１１は、元の各該当圃場の画像をそれぞれに対応する補正後の新画像に差し替えることで圃場マップを生成し、それをデータベース１１５に提供する。このようにしてデータベース１１５に格納された圃場マップは、その後、生長分析システム使用部１１３に提供され、農業支援解析に使用される。
【００９５】
作物収穫画像生成部１１２は、生長差補正部１１０が計算した各圃場の補正した作物生理特徴量を有する圃場によって、リモートセンシング画像の各対応圃場を置き替え、収穫時期及び収量のマップを生成する。
【００９６】
図１５は、本発明の実施例における収穫マップの処理の各ケースを表す説明図である。
【００９７】
図１５（ａ）に示した元の生理特性生長マップ、及び、図１５（ｂ）に示した収穫マップは、両方とも、時間的に統一された画像である。図１５（ｂ）に示した処理後の画像は圃場ごとに色がつけられ、その色が、収穫時期によって収穫順番及び収穫量を表すようになる。作物収穫画像生成部１１２は、元の各該当圃場の画像をそれぞれに対応する補正後の新画像に差し替えることで圃場マップを生成し、それをデータベース１１５に提供する。このようにしてデータベース１１５に格納された圃場マップは、その後、物流配車システム使用部１１４に提供され、収穫プラン作成に使用される。
【００９８】
図１４Ａは、本発明の実施例におけるデータベース１１５の例を表す説明図である。
【００９９】
データベース１１５には、時系列衛星画像データ１４０１、ＧＩＳデータ（位置及び形状データ）１４０２、植物生長データ（データベース生成部でＧＩＳデータ及び衛星、気象データから算出したデータ）１４０３、気象データ（本実施例では表面温度データ）１４０４、土壌及び他の相関データ１４０５、並びに、処理パラメータ１４０６を格納する。図示は省略するが、データベース１１５は上記以外のデータ（例えば災害データ等）を含んでもよい。
【０１００】
図１４Ｂは、本発明の実施例におけるデータベース１１５に格納される時系列衛星画像データ１４０１の例を表す説明図である。
【０１０１】
リモートセンシング画像データ取得部１０４が取得したリモートセンシング画像データ１０１ｃが、時系列衛星画像データ１４０１としてデータベース１１５に格納される。時系列衛星画像データ１４０１は、同じ場所（例えば作物管理システム１００による管理対象の圃場を含む農場）をそれぞれ異なる時期に撮影した複数の画像群であり、例えば、複数年の年間を通した複数リモートセンシング画像を含んでもよい。
【０１０２】
時系列衛星画像データ１４０１は、生長モデルの補正（Ｓ３０７）、及び、生長段階が統一された画像の生成（Ｓ３０８〜Ｓ３１０）に用いられる。
【０１０３】
図１４Ｃは、本発明の実施例におけるデータベース１１５に格納されるＧＩＳデータ１４０２の例を表す説明図である。
【０１０４】
気象・ＧＩＳデータ取得部１０３が取得した気象・ＧＩＳデータ１０１ｂのうち、ＧＩＳデータが、ＧＩＳデータ１４０２としてデータベース１１５に格納される。ＧＩＳデータ１４０２には、各圃場の位置データ及び形状データが登録される。ＧＩＳデータ１４０２は、圃場ＩＤ、ＧＰＳデータ、及び形状データを含む。圃場ＩＤには圃場の識別情報が登録される。ＧＰＳデータには、圃場の位置を示す緯度及び経度が登録される。形状データには圃場の輪郭を示すデータ（圃場ポリゴン）が登録される。
【０１０５】
ＧＩＳデータ１４０２を参照することによって、時系列衛星画像データ１４０１に含まれる画像から、圃場の画像の輪郭である圃場のポリゴンを抽出することができる。
【０１０６】
図１４Ｄは、本発明の実施例におけるデータベース１１５に格納される植物生長データ１４０３の例を表す説明図である。
【０１０７】
植物生長データ１４０３は、植生指標時系列データ及び作物生長基準データを含む。
【０１０８】
植生指標時系列データには、各圃場の植生の状況を示す特徴量（植生指標）が含まれる。本実施例では、特徴量としてＮＤＶＩが保持される。例えば、圃場１、圃場２といった圃場ごとの、時期ごとのＮＤＶＩが植生指標時系列データとして保持される。
【０１０９】
成長ステージ基準データ４６０は、例えば大豆の各重要な時期に達したと判断される積
算温度、及び各成長ステージの植物特性が登録される。
【０１１０】
作物生長基準データは、各生長ステージの積算温度及び植物特性を含む。例えば、大豆の生長ステージには、開花する開花期、莢を付ける着莢期、種が完全な状態になるまでの時期、及び葉の５０％以上黄色になり、種が完熟する完熟期がある。サトウキビなど、他の作物にも同様にそれぞれの生長ステージがある。作物生長基準データとして、対象の作物がそれぞれの生長ステージに達するまでの積算温度、及び、各生長ステージの植物特性が保持される。
【０１１１】
図１４Ｅは、本発明の実施例におけるデータベース１１５に格納される気象データ１４０４の例を表す説明図である。
【０１１２】
気象データ１４０４は、例えば、降水量データ及び気温データを含む。降水量データとして、ある地域（例えば管理対象の農場を含む地域）の所定期間の降水量が保持される。降水量データは、例えば、気象庁等の降水量のデータを利用することができる。
【０１１３】
また、気温データとして、ある地域の所定期間の最高気温及び最低気温が保持される。例えば、最高気温及び最低気温は、ある地域の圃場の表面温度であってもよい。気温データは、例えば、ＭＯＤＩＳ表面温度衛星によるデータを利用することができる。
【０１１４】
また、図１４Ｅでは省略しているが、気象データ１４０４は、さらに日射量のデータ等を含んでもよい。
【０１１５】
気温データは、積算温度の算出に利用することができる。また、何らかのモデルを使用してＬＡＩを算出する場合には、気象データ１４０４に含まれる気温データ、日射量データ及び降水量データ等を利用することができる。
【０１１６】
図１４Ｆは、本発明の実施例におけるデータベース１１５に格納される土壌及び他の相関データ１４０５の例を表す説明図である。
【０１１７】
図１４Ｆに示す土壌及び他の相関データ１４０５は、土壌データ及び作物栽培データを含む。
【０１１８】
土壌データには、ある地域の各圃場の土壌の成分のうち、植物の成長に必要な成分の植物の成長に寄与する係数が登録される。図１４Ｆに示す土壌データは、圃場ごとのリン酸係数を含む。
【０１１９】
作物栽培データには、各圃場の作物が所定の段階に達した時期（例えば日付）が登録される。図１４Ｆに示す作物栽培データは、圃場ごとの植付時期、開花時期、着莢時期及び収穫時期を含む。これらの時期は、各圃場の作物の種類（例えば大豆又はサトウキビ等）に応じたものが登録される。これらの時期は、例えば、農家から得られる情報に基づいて登録されてもよい。
【０１２０】
図１４Ｇは、本発明の実施例におけるデータベース１１５に格納される処理パラメータ１４０６の例を表す説明図である。
【０１２１】
処理パラメータ１４０６には、各圃場の生長モデルに関するパラメータ等が登録される。
【０１２２】
処理パラメータ１４０６は、圃場ごとの生長ステージパラメータ、生長差及び補正値を含む。生長ステージパラメータは、植物の生長ステージに対応するパラメータ（例えば、各ステージの生長モデルを表すためのパラメータ等）が登録される。生長差には、各圃場の積算温度と基準圃場の積算温度との差が生長差として登録される。補正値には、各圃場の特徴量（例えばＬＡＩ）を基準圃場の積算温度に対応させるための補正値が登録される。
【０１２３】
以上に説明した本発明の態様の代表的な例をまとめると次の通りとなる。すなわち、プロセッサ（例えばＣＰＵ１２２）と、プロセッサがアクセスする記憶装置（例えばＲＡＭ１２３又は記憶部１２４）と、を有する作物管理システムであって、記憶装置は、農場の複数の区画（例えば圃場）を過去のそれぞれ異なる時期に撮影した複数の画像を含む画像データ（例えばリモートセンシング画像データ取得部１０４が取得してデータベース１１５に格納した画像データ）と、過去のそれぞれ異なる時期に計測された区画ごとの気象データ（例えば気象・ＧＩＳデータ取得部１０３が取得してデータベース１１５に格納した気象データ）と、を保持し、プロセッサは、気象データに基づいて、作物の生長状況を示す指標を推定するための生長モデルを生成し（例えばＳ３０４）、画像データに基づいて、区画ごとに、生長モデルを補正（例えばＳ３０７）してもよい。
【０１２４】
これによって、実際の圃場の状況が反映された精度の高い生長モデルを得ることができる。
【０１２５】
このとき、プロセッサは、気象データに基づいて、時間ごとの作物の生長状況を示す指標を計算し、時間ごとの気温を取得し、時間ごとの作物の生長状況を示す指標及び時間ごとの気温に基づいて、積算温度に対する作物の生長状況を示す指標を、生長モデルとして生成し、積算温度に基づいて、作物の生長段階の区分（例えば図５のＲ１〜Ｒ４によって区切られた生長ステージ）を抽出してもよい。
【０１２６】
これによって、同種類の植物の植える時間などの差異によって生長状況が異なる問題が解決できる。
【０１２７】
また、プロセッサは、生長段階の区分ごとに、積算温度と作物の生長状況を示す指標との関係を示す関数を生成して、関数のパラメータを生長モデルとして記憶装置に保持してもよい。
【０１２８】
これによって、数学的に生長状況を表現することができ、表現された生長状況を、生育段階が統一された画像の生成のベースとすることができる。
【０１２９】
また、プロセッサは、画像データに基づいて、区画ごとに、作物の量及び活力を示す植生指標を計算し、植生指標に基づいて、作物の生長状況を示す指標を計算し、植生指標に基づいて計算した作物の生長状況を示す指標を用いて、生長モデルを補正してもよい（例えばＳ３０７、図７）。
【０１３０】
これによって、画像から得られる実際の植生の状況が反映され、より精度の高い生長状況の推定を行うことができる。
【０１３１】
また、作物の生長状況を示す指標は、葉面積指数、バイオマス又は糖度のいずれかであり、植生指標は、可視領域の赤の反射率と近赤外領域の反射率とに基づく正規化差植生指数であり、積算温度は、作物の種類に応じた基準温度に基づく有効積算温度であってもよい。
【０１３２】
これによって、作物の管理に有効な情報を得ることができる。
【０１３３】
また、プロセッサは、補正された生長モデルに基づいて、生長段階が統一された、複数の区画における作物の生長状況を示す画像を生成してもよい（例えばＳ３１０、図１３（ｂ））。
【０１３４】
これによって、植え付け時期の違いの影響が取り除かれた生長状況を視覚的に把握することが可能になる。
【０１３５】
また、プロセッサは、複数の区画のうち第１の区画を基準区画として選択し（例えばＳ１０６２）、第１の時点の、複数の区画のうち前記第１の区画とは異なる第２の区画における積算温度に基づく生長段階（例えば図９〜図１２の黒丸の点の生長段階）と、第１の時点の、第１の区画における積算温度に基づく生長段階（例えば図９〜図１２の白丸の点の生長段階）との相違を計算し（例えばＳ１０６３）、第２の区画の生長モデルと、計算された生長段階の相違と、に基づいて、第２の区画における作物の生長状況を示す指標を、第２の区画における生長段階が第１の区画における生長段階と同一となる時点の、第２の区画における作物の生長状況を示す値に補正し（例えばＳ１０６４、図９〜図１２）、第１の時点の、第１の区画における作物の生長状況を示す指標の値と、補正された第２の区画における作物の生長状況を示す指標の値と、に基づいて、複数の区画における作物の生長状況を示す画像を生成してもよい（例えばＳ３１０）。
【０１３６】
このとき、プロセッサは、生長段階の区分ごとに、積算温度と作物の生長状況を示す指標との関係を示す関数を生成して、関数のパラメータを生長モデルとして記憶装置に保持し、第１の区画における生長段階を含む区分から第２の区画における生長段階を含む区分までの全ての区分の関数のパラメータを使用して、第２の区画における作物の生長状況を示す指標の補正をしてもよい。これは、図１１の例においてＲ２−Ｒ３の生長ステージの関数パラメータ及びＲ３−Ｒ４の生長ステージの関数パラメータを使用してＬＡＩが補正され（式（４））、図１２の例においてＲ１−Ｒ２の生長ステージの関数パラメータ及びＲ２−Ｒ３の生長ステージの関数パラメータを使用してＬＡＩが補正されることに相当する。
【０１３７】
これによって、生長状況を示す指標の値が、基準圃場の生長段階に整合するように精度よく計算される。
【０１３８】
また、プロセッサは、他の区画の積算温度との関係が所定の条件を満たす積算温度を有する区画を基準区画として選択してもよい（例えばＳ１０６２、式（３））。
【０１３９】
これによって、適切な基準区画を選択することができる。
【０１４０】
また、プロセッサは、第２の区画に対応するポリゴン内の複数のピクセルに同一の補正を行うことによって、複数の区画における作物の生長状況を示す画像を生成してもよい。
【０１４１】
これによって、ポリゴン内（例えば同一圃場内）のピクセルの値のばらつきは維持したまま、全体のレベルを、植え付け日の相違の影響が取り除かれるように補正することができる。
【０１４２】
また、プロセッサは、補正された生長モデルに基づいて、複数の区画における作物の収穫順序を示す画像を生成してもよい（例えばＳ３１２、図１５）。
【０１４３】
これによって、収穫プランニングが容易になり、収穫機器の待機等によるコストの増加を抑えることができる。
【０１４４】
また、本発明の態様の代表的な例である作物管理システムは、リモートセンシング画像取得部、気象データ取得部、作物生長データ環境データで取得したリモートセンシング画像、気温データ、作物生育データ環境データを用いて参考になる作物データを保存し、作物生育、環境データと気象データを併用して作物生理特性の有効積算温度によりの変化推移を生成する生長モデル生成部と、この有効積算温度ごとの指標変化推移により作物生長の生長ステージを抽出する生長ステージ抽出部と、抽出した生長ステージに対して解析し各ステージの生長をパラメータ化する生長ステージ解析部と、リモートセンシング特徴量を用いた生長モデルを補正する生長ステージ補正部と、各圃場に対して生育差を計算する生育差計算部と、生育差を補正する生育差補正部と、生育補正に応じて生育ステージ統一画像を生成する生育ステージ統一画像生成部と、生育モデル補正に応じて作物収穫画像生成部と、生育ステージ統一画像を用いた生育分析システム使用部と、作物収穫画像を用いた物流配車システム使用部と、を有してもよい。
【０１４５】
リモートセンシング画像取得部は、時系列リモートセンシング画像を取得し、その画像を反映した作物特徴量を抽出してもよい。気象データ取得部としては、分析対象の地域または圃場の気温データを取得してもよい。作物生育、環境データ取得部としては、作物自身の生理特性パラメータ及び分析対象の土壌などの環境データを取得してもよい。そして、リモートセンシング観測装置から送信された、リモートセンシング画像と、作物生長、環境データを用いて、対象地域を分析する作物管理システムにおいて、該観測装置から送信された画像をデータベースに保存してもよい。
【０１４６】
また、観測装置から送信された、リモートセンシング画像と、作物生長、環境データを用いて、対象地域を分析する作物管理システムは、該観測装置から送信された画像を保存するデータベースと、該データベースに保存されたリモートセンシング画像のうち、時系列でとれる画像と作物生長、環境データを用いて、予め設定された画像取得期間内の時系列で演算処理して、気温データを合成し生長モデルを生成する生長モデル生成部を有してもよい。このシステムにおいて、生長モデル生成部は、通常の作物指標／時間の変化推移ではなく、時間データと気温データを併用し、作物の植え付け時間からの毎日有効温度の積算値を算出し、作物本質を示す有効積算温度分析モデルに変換し、この積算温度により変化を分析してもよい。有効積算温度は、同種類の作物に対して、生長段階を反映することができるので、本発明の有効積算温度の使用を通じて、同種類の植物の植える時間などの差異によって生長状況が異なる問題が解決できると考えられる。
【０１４７】
また、上記の作物管理システムは、生長モデル生成部が生成した作物指標／積算温度の変化推移曲線から生長ステージを抽出する生長ステージ抽出部を有してもよい。このシステムにおいて、生長ステージ抽出部が、すべての同種類の作物栽培特性のステージを明らかにすることによって、精密な分析が可能になる。
【０１４８】
また、上記の作物管理システムは、生長ステージ抽出部が抽出したステージを解析し、生物性と作物特徴量の関係を考量する上各ステージの生長モデルを関数パラメータ化する生長ステージ解析部を有してもよい。このシステムにおいて、生長ステージ解析部は、生長モデルを本質な作物生長特性と合わせ、数学的に生長状況を示すことができ、これを生育統一画像の生成のベースとすることができる。
【０１４９】
また、上記の作物管理システムは、生長モデル生成部が生成した作物指標／積算温度の変化推移曲線から生長ステージ抽出部が抽出したステージを解析し、生物性と作物特徴量の関係を考量する各ステージの生長モデルを関数パラメータ化し、本年度のデータから生長モデルを補正してもよい。このシステムにおいて、過去データと本年度データの融合にって、より精密な作物生長モデルを生成でき、生育補正の精度を高めることができる。
【０１５０】
また、作物管理システムは、上記のように生成した補正後の作物生長モデルを用いて、各解析単位に対して基準解析単位との生育差を計算し、基準解析単位に生育差を補正し、補正をかけた各解析単位を元の画像の対応解析単位を替えて、生育ステージ統一画像と生長モデル時間統一収穫マップを生成してもよい。このシステムにおいて、作物生長モデルの活用を通じて、生育ズレがある各解析単位の生育補正ができ、生育統一画像を生成でき、収量予測及び生長状況分析などによって高い精度と信頼度が得られる。
【０１５１】
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明のより良い理解のために詳細に説明したのであり、必ずしも説明の全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることが可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【０１５２】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によってハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによってソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、不揮発性半導体メモリ、ハードディスクドライブ、ＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）等の記憶デバイス、または、ＩＣカード、ＳＤカード、ＤＶＤ等の計算機読み取り可能な非一時的データ記憶媒体に格納することができる。
【０１５３】
また、制御線及び情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線及び情報線を示しているとは限らない。実際にはほとんど全ての構成が相互に接続されていると考えてもよい。
【符号の説明】
【０１５４】
１００：作物管理システム
１０１ａ：作物生長データ、環境データ
１０１ｂ：気象、ＧＩＳデータ
１０１ｃ：時系列リモートセンシング画像データ
１０２：作物生理、環境データ取得部
１０３：気象、ＧＩＳデータ取得部
１０４：リモートセンシング画像データ取得部
１０５：生長モデル生成部
１０６：生長ステージ抽出部
１０７：生長ステージ解析部
１０８：生長モデル補正部
１０９：生長差計算部
１１０：生長差補正部
１１１：生長ステージ統一画像生成部
１１２：作物収穫画像生成部
１１３：生長分析システム使用部
１１４：物流配車システム使用部
【図１】
2020140347.tif 000007
【図２】
2020140347.tif 000008
【図３】
2020140347.tif 000009
【図４】
2020140347.tif 000010
【図５】
2020140347.tif 000011
【図６】
2020140347.tif 000012
【図７】
2020140347.tif 000013
【図８】
2020140347.tif 000014
【図９】
2020140347.tif 000015
【図１０】
2020140347.tif 000016
【図１１】
2020140347.tif 000017
【図１２】
2020140347.tif 000018
【図１３】
2020140347.tif 000019
【図１４Ａ】
2020140347.tif 000020
【図１４Ｂ】
2020140347.tif 000021
【図１４Ｃ】
2020140347.tif 000022
【図１４Ｄ】
2020140347.tif 000023
【図１４Ｅ】
2020140347.tif 000024
【図１４Ｆ】
2020140347.tif 000025
【図１４Ｇ】
2020140347.tif 000026
【図１５】
2020140347.tif 000027