JP7394458B2 - 推定装置、推定方法、及びプログラム - Google Patents

Description

特許法第３０条第２項適用 発行日 令和元年６月１０日 刊行物 ＪＳＴｎｅｗｓ ２０１９年６月号（Ｊｕｎｅ ２０１９） Ｐ５－７（ウェブ版ＵＲＬ：ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔ．ｇｏ．ｊｐ／ｐｒ／ｊｓｔ－ｎｅｗｓ／ｂａｃｋｎｕｍｂｅｒ／２０１９／２０１９０６／ｉｎｄｅｘ．ｈｔｍｌ ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔ．ｇｏ．ｊｐ／ｐｒ／ｊｓｔ－ｎｅｗｓ／ｂａｃｋｎｕｍｂｅｒ／２０１９／２０１９０６／ｐｄｆ／２０１９＿０６．ｐｄｆ） 国立研究開発法人科学技術振興機構 発行

特許法第３０条第２項適用 ウェブサイト公開日 令和元年７月９日 ウェブサイトアドレス ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔａｇｅ．ｊｓｔ．ｇｏ．ｊｐ／ｂｒｏｗｓｅ／ｊｓｔｎｅｗｓ／－ｃｈａｒ／ｊａ ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔａｇｅ．ｊｓｔ．ｇｏ．ｊｐ／ａｒｔｉｃｌｅ／ｊｓｔｎｅｗｓ／２０１９／６／２０１９＿５／＿ａｒｔｉｃｌｅ／－ｃｈａｒ／ｊａ／

特許法第３０条第２項適用 発行日 令和元年８月６日 刊行物 国立研究開発法人科学技術振興機構 戦略的創造研究推進事業（さきがけ） 「情報科学との協働による革新的な農作物栽培手法を実現するための技術基盤の創出」 研究領域 領域活動・評価報告書－２０１８年度終了研究課題－国立研究開発法人科学技術振興機構 発行 （ウェブ版ＵＲＬ）ｈｔｔｐｓ：／／ｗｗｗ．ｊｓｔ．ｇｏ．ｊｐ／ｋｉｓｏｋｅｎ／ｐｒｅｓｔｏ／ｅｖａｌｕａｔｉｏｎ／ｐｏｓｔｅｒｉｏｒｉ／ｈ３０／ＪＳＴ＿１１１２０７５＿２０１８＿ＰＥＥ．ｐｄｆ

本発明は、作物における個体毎の生長領域を推定する推定装置、推定方法、及びプログラムに関する。

画像から作物の生育状態を観察する研究が作物を問わず古くから行われている。これらの研究は、特定の面積にある複数の個体をまとめて観察する方法、つまり、群落の観察によるものがほとんどである。しかし、究極の精密農業である個体管理を想定すると、画像によって作物個体の生育状態を観察する方が望ましい。画像から個々の作物における生長領域を定義することができれば、個体の生育状態を把握しやすくなり、有用な情報となり得る。

非特許文献１には、デジタルカメラで撮影した圃場作物の画像から作物個体の位置を検出する技術が開示されている。非特許文献１では、作物としてキャベツとハクサイを対象にして所定の処理アルゴリズムにより作物個体の位置を決定する。非特許文献２には、個体の位置からその個体の生長領域を決定する技術として、ボロノイ分割を使用した技術が開示されている。

北海道農業研究センター １９９５年の成果情報、ファジィ画像処理による作物配列の認識、［online］、農研機構、［令和２年４月２日検索］，インターネット＜https://www.naro.affrc.go.jp/project/results/laboratory/harc/1995/cryo95-016.html＞ Mithen et al. 1984. Growth and mortality of individual plants as a function of "available area". Oecologia vol.62, 57-60.

しかしながら、非特許文献１に記載の検出方法は、作物が格子状に配置されていることを前提とする。圃場に植えられる作物の全般においては、作物が格子状に配置されない場合も多く、圃場作物の全般に非特許文献１の技術を適用しようとすれば、適用することができる作物が限定されてしまう。非特許文献２に記載の生長領域の決定方法は、作物の個体の大きさを考慮しておらず、個体が利用できる生長領域の実態が反映されていないという問題があった。

本発明は、このような事情に鑑みてなされたもので、作物が格子状に配置されていない場合であっても、作物の個体の大きさを考慮した生長領域を推定することができる推定装置、推定方法、及びプログラムを提供する。

本発明の一態様によれば、推定装置は、対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定する個体間境界決定部と、前記個体間境界決定部によって決定された前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出する生長領域算出部と、を備える。

前記対象画像を取得する画像取得部と、前記画像取得部によって取得された前記対象画像に基づいて、前記対象場所において作物が植生している領域である植生部分を抽出する植生部分抽出部と、前記植生部分抽出部によって抽出された前記植生部分に含まれる作物個体の個体数を推定する個体数推定部と、前記個体数推定部によって推定された前記植生部分に含まれる作物個体の個体数に基づいて、前記植生部分に含まれる作物個体の位置である個体位置を決定する個体位置決定部と、を更に備え、前記個体間境界決定部は、前記個体位置決定部によって決定された前記個体位置を、前記対象場所において隣り合う作物のそれぞれの位置とする。

前記個体数推定部は、前記植生部分における面積、及び形状に関する特徴に基づいた分類を行うことにより、前記植生部分における個体数を推定するようにしてもよい。

前記個体数推定部は、前記植生部分における形状に関する特徴として、当該形状における円形度、及び最小外接矩形の長辺長と短辺長とを用いるようにしてもよい。

前記個体数推定部は、サポートベクタマシンに、前記植生部分における形状に関する特徴を入力して得られる分類結果を用いて、前記植生部分における個体数を推定するようにしてもよい。

前記個体位置決定部は、前記個体数推定部により前記植生部分に複数の個体が含まれると推定された場合、当該植生部分における主軸に垂直な方向に存在する植生部分の画素数の分布に基づいて前記個体位置を決定するようにしてもよい。

前記個体位置決定部は、前記植生部分に含まれる画素の分布を主成分分析して得られる第１主成分軸を、前記植生部分における主軸とするようにしてもよい。

前記個体位置決定部は、前記植生部分を貫通する直線のうち、前記植生部分を貫通する距離が最も長い直線を、前記植生部分における主軸とするようにしてもよい。

前記個体間境界決定部は、隣り合う２つの作物個体の個体位置を焦点とし、当該２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を頂点とする双曲線を、前記個体間境界とするようにしてもよい。

前記個体間境界決定部は、隣り合う２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を通る線を、前記個体間境界とするようにしてもよい。

本発明の一態様によれば、推定方法は、個体間境界決定部が、対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定し、生長領域算出部が、前記個体間境界決定部によって決定された前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出する。

本発明の一態様によれば、プログラムは、コンピュータに、対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定させ、前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出させる、ためのプログラムである。

本発明によれば、作物が格子状に配置されていない場合であっても、作物の個体の大きさを考慮した生長領域を推定することができる。

実施形態に係る推定装置１の構成の例を示すブロック図である。 実施形態に係る植生部分抽出部１２が行う処理を説明する図である。 実施形態に係る個体数推定部１３が行う処理を説明する図である。 実施形態に係る個体数推定部１３が行う処理を説明する図である。 実施形態に係る個体数推定部１３が行う処理を説明する図である。 実施形態に係る個体位置決定部１４が行う処理を説明する図である。 実施形態に係る個体位置決定部１４が行う処理を説明する図である。 実施形態に係る個体位置決定部１４が行う処理を説明する図である。 実施形態に係る個体間境界決定部１５が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。 実施形態に係る推定装置１が行う処理の流れを示すフロー図である。

以下、本発明の実施形態を説明する。図１は、実施形態に係る推定装置１の構成の例を示すブロック図である。推定装置１は、圃場に植生する個々の作物（作物個体）の生長領域を推定するコンピュータ装置である。作物の生長領域とは、作物の生長に利用することが可能な領域である。推定装置１は、圃場の作物が発芽、或いは萌芽してから、定期的（例えば、１週間毎）に作物個体の生長領域を推定する。これにより、作物個体それぞれの生育状況に応じた生長領域の実態を把握する。

推定装置１は、例えば、画像取得部１０と、画像合成部１１と、植生部分抽出部１２と、個体数推定部１３と、個体位置決定部１４と、個体間境界決定部１５と、個体生長領域算出部１６と、記憶部１７とを備える。

画像取得部１０は、空撮画像を取得する。空撮画像は上空からみた圃場に植生している作物の様子が撮像された画像であり、例えば、デジタルカメラなどの撮像装置が搭載されたドローンなどを用いて撮像される。この場合、ドローンは、圃場の上空を往復飛行しながら、定期的（例えば、１～２秒毎）に静止画像（空撮画像）を連続的に撮像し、複数に分けて圃場の全体が撮像されるようにする。画像取得部１０は、取得した空撮画像を画像合成部１１に出力する。

画像合成部１１は、画像取得部１０によって得られた空撮画像に基づいて、圃場全体の合成画像（以下、単に「合成画像」という）を生成する。画像合成部１１は、複数の空撮画像の中心から所定の範囲にある画像をつなぎ合わせたオルソ画像を生成し、生成したオルソ画像を合成画像とする。この場合において、画像合成部１１は複数の空撮画像から得られる土壌や作物個体の高さ情報等に基づいて画像に生じる歪を補正するようにしてもよい。

本実施形態において、合成画像は、後述する画像処理に要求される精度を満たす空間分解能であることが好ましい。後述する画像処理とは、画像合成部１１、植生部分抽出部１２、個体数推定部１３、個体位置決定部１４、及び個体生長領域算出部１６により行われる画像処理である。合成画像に求められる空間分解能は、例えば、１［ｃｍ／ｐｉｘｅｌ］、すなわち、画像の１［ピクセル］に、実際の空間における１［ｃｍ］程度の範囲が撮像されるような分解能である。この場合、合成画像に要求される空間分解能に応じて、例えば、空撮画像の空間分解能が決定され、決定された空間分解能により空撮画像が撮像される。

以下、合成画像を用いて、推定装置１が生長領域を推定する場合を例に説明する。しかしながらこれに限定されることはない。推定装置１は、少なくとも生長領域を推定する対象とする作物とその周囲の領域が撮像された画像を用いて生長領域を推定できればよい。例えば、推定装置１は、圃場の一部のみが撮像された画像を用いて、その画像に撮像されている作物個体の生長領域を推定するようにしてもよい。また、推定装置１は、空撮画像に生じ得る歪が低減される観点から、オルソ画像を用いて生長領域を推定することが好ましいが、これに限定されることはない。推定装置１は、合成前の空撮画像を用いて作物個体の生長領域を推定してもよい。

植生部分抽出部１２は、合成画像から植生部分を抽出する。例えば、植生部分抽出部１２は、合成画像における植生部分と土壌部分との色の相違に基づいて両者を分離する。この場合、植生部分抽出部１２は、植生部分と土壌部分の色の特徴に基づいて、画素ごとの色情報（例えば、ＲＧＢ値）を２値化する。

例えば、植生部分抽出部１２は、植生部分と想定される標準的な色（例えば、緑色：ＲＧＢ値（０、１２８、０））と、色空間上の距離が所定の範囲内にある色の画素を特定の第１色（例えば、白色）に変換する。植生部分抽出部１２は、植生部分における標準的な色と、色空間上の距離が所定の範囲外にある色の画素を特定の第２色（例えば、黒色）に変換する。植生部分抽出部１２は、合成画像において、第１色（例えば、白色）に変換された部分を、植生部分として抽出する。植生部分抽出部１２は、合成画像において植生部分を示す情報（例えば、画素の座標値）を、個体数推定部１３に出力する。

ここで、合成画像から抽出される植生部分としては、１つの作物個体からなる植生部分として抽出されるもの（図２の植生部分Ｒ２１参照）のみならず、２つ以上の作物個体がつながった状態で抽出されるもの（図２の植生部分Ｒ２２参照）があり得る。本実施形態の推定装置１では、作物個体の生長領域を推定する。このため、推定装置１では、合成画像から２つ以上の作物個体がつながった状態で植生部分が抽出された場合、抽出された植生部分を、１つの作物個体に相当する植生部分ごとに分離する必要がある。

この対策として、推定装置１では、個体数推定部１３により、植生部分に含まれる作物の個体数を推定する。そして、推定装置１では、作物個体の生長領域を推定する方法を、植生部分に含まれる作物の個体数に応じて、それぞれの個体数に適した方法を用いて作物個体の生長領域を推定する。これにより、推定装置１は、合成画像から２つ以上の作物個体がつながった状態で植生部分が抽出された場合であっても、１つの作物ごとに、その作物個体の生長領域を推定することが可能となる。

個体数推定部１３は、植生部分抽出部１２によって抽出された植生部分に含まれる作物の個体数を推定する。一般に、植生部分の面積が大きい程、その植生部分に含まれる作物の個体数が多いことが推定される。また、１つの作物個体からなる植生部分は、略円形の形状となり、複数の作物個体の植生部分がつながって１つの植生部分が構成される場合、個々の作物個体の略円形の形状における周縁部が重なった形状と面積が大きい程、その植生部分に含まれる作物の個体数が多いことが推定される。

個体数推定部１３は、この性質を利用し、植生部分のそれぞれの大きさや形状に基づいて、当該植生部分に含まれる作物の個体数を推定する。例えば、個体数推定部１３は、サポートベクタマシンによる分類結果を利用して、植生部分に含まれる作物の個体数を２つのクラスに分離する。この場合、サポートベクタマシンに入力する説明変数は、植生部分の面積、及び形状に関する特徴である。形状に関する特徴は、任意の特徴であってよいが、特に、植生部分に含まれる作物個体の個体数に応じた形状の特徴であることが好ましい。これにより、形状に関する特徴を入力することにより、サポートベクタマシンから得られる結果を、植生部分に含まれる作物個体の個体数に応じて分離された結果とすることが可能である。

個体数推定部１３は、例えば、植生部分の形状に関する特徴として、植生部分の周囲長、円形度、最小外接矩形の長辺長と短辺長、回帰円からのＲＭＳ（Root Mean Square）誤差、及び回帰円の半径などを用いる。周囲長は、植生部分の周囲の長さであり、例えば、植生部分の輪郭（エッジ）に相当する画素を抽出し、抽出した画素の数に応じて導出される。円形度は、植生部分の形状における円形らしさを示す指標であり。例えば、円形度は、植生部分の形状が円形であった場合の面積と周囲長の関係に対する、実際の植生部分の面積と周囲長の関係、に応じて導出される。最小外接矩形は、植生部分に外接する最小の矩形である。回帰円は、植生部分の形状に当てはめることができる円である。

まず、個体数推定部１３は、合成画像から抽出された植生部分の全てを対象として、植生部分の面積、及び形状に関する特徴をサポートベクタマシンに入力する。これにより、サポートベクタマシンは、入力された説明変数に基づいて、例えば、植生部分を、１つの作物個体からなる植生部分と、複数（２つ以上）の作物個体からなる植生部分とに分離する（図３参照）。

次に、個体数推定部１３は、サポートベクタマシンにより、複数（２つ以上）の作物個体からなる植生部分に分類された植生部分を対象として、植生部分の面積、及び形状に関する特徴をサポートベクタマシンに入力する。これにより、サポートベクタマシンは、入力された説明変数に基づいて、例えば、植生部分を、２つの作物個体からなる植生部分と、３つ以上の作物個体からなる植生部分とに分離する（図４参照）。

さらに、個体数推定部１３は、サポートベクタマシンにより、３つ以上の作物個体からなる植生部分に分類された植生部分を対象として、植生部分の面積、及び形状に関する特徴をサポートベクタマシンに入力する。これにより、サポートベクタマシンは、入力された説明変数に基づいて、例えば、植生部分を、３つの作物個体からなる植生部分と、それ以上（４つ以上）の作物個体からなる植生部分とに分離する（図５参照）。

このように、例えば、個体数推定部１３は、分離する対象を絞り込みながら、サポートベクタマシンから得られる結果を用いて、植生部分に含まれる作物個体の個体数を推定する。

なお、植生部分に含まれる作物の個体数が増加した場合、サポートベクタマシンを用いた分離の精度が劣化することが考えられる。例えば、４つの作物個体からなる植生部分と、それ以上（５つ以上）の作物個体からなる植生部分における植生部分の面積、及び形状に関する特徴が類似することが考えられる。植生部分を個体数ごとに分離し得る特徴が見いだされなかった場合、サポートベクタマシンを用いた分離の精度が劣化する可能性がある。

この対策として、個体数推定部１３は、サポートベクタマシンにより分離する個体数の上限を予め設定するようにしてもよい。例えば、個体数推定部１３は、サポートベクタマシンにより精度よく分離可能な範囲において、サポートベクタマシンを用いた個体数の推定を行う。例えば、個体数推定部１３は、個体数が既知である植生部分を、サポートベクタマシンに分離させた場合に、実際の個体数と整合するか否かに基づいて、サポートベクタマシンを用いた分離を行う個体数の上限を決定する。

なお、上記では、個体数推定部１３がサポートベクタマシンを用いた分類を行う場合を例に説明したが、これに限定されない。個体数推定部１３は、少なくとも、植生部分の面積、及び形状に基づいて分類した結果を用いて、個体数を推定すればよい。例えば、個体数推定部１３は、植生部分の面積及び形状の組合せに、その植生部分の個体数が対応づけられたテーブルを参照することにより、植生部分に含まれる作物個体の個体数を推定するようにしもよい。

個体位置決定部１４は、植生部分における作物個体の個体位置を決定する。個体位置は、生長領域を推定する際に用いられる作物個体の位置である。個体位置決定部１４は、植生部分に含まれる作物個体が、
（ケース１）１個体である場合
（ケース２）２個体又は３個体である場合
（ケース３）４個体以上である場合
のそれぞれに応じて、互いに異なる方法を用いて、個体位置を決定する。

まず、（ケース１）植生部分に含まれる作物個体が１個体である場合において、個体位置決定部１４が、植生部分における作物個体の個体位置を決定する方法について説明する。植生部分に含まれる作物個体が１個体である場合、個体位置決定部１４は、植生部分の重心位置を、その植生部分における作物個体の個体位置とする（図６参照）。

次に、（ケース２）植生部分に含まれる作物個体が２個体又は３個体である場合において、個体位置決定部１４が、植生部分における作物個体の個体位置を決定する方法について説明する。まず、植生部分に含まれる作物個体が２個体である場合について説明する。個体位置決定部１４は、植生部分における主軸を定義する。ここでの主軸は、植生部分に含まれる二つの作物個体の個体位置と想定される点を通る直線である（図７の主軸ＳＧ参照）。

植生部分は、作物個体の個体位置を中心として同心円状に広がる傾向にある。そして、植生部分に二つの作物個体が含まれている場合には、二つの作物個体それぞれに相当する二つの略円形が、互いが隣接する周縁部において重なり、楕円の中央部分が窪んだ、ひょうたんのような形状となる。この場合における二つの円の中心を通る直線が、本実施形態における主軸に相当する。

上記のように定義される主軸は、植生部分における画素の分布のばらつきを最大にする軸とみなすことができる。この性質を利用して、個体位置決定部１４は、例えば、主成分分析を用いて主軸を決定する。この場合、植生部分における画素の分布のばらつきを主成分分析して得られる第１主成分軸が、本実施形態における、植生部分の主軸に相当する。

具体的に、個体位置決定部１４は、まず、植生部分の輪郭を構成する画素の位置座標の平均座標を原点とする座標系（以下、ｘｙ座標系という）を定義する。個体位置決定部１４は、植生部分の輪郭を構成する画素のそれぞれの位置座標を、ｘｙ座標系を用いて表す。個体位置決定部１４は、ｘｙ座標系の原点を通る直線ＴＳ（ａｘ＋ｂｙ＝０）を定義する。ａ、及びｂは任意の実数である。個体位置決定部１４は、植生部分の輪郭を構成する画素のｘｙ座標（ｘ_ｉ、ｙ_ｉ）から直線ＴＳまでの最短距離を誤差ｅ_ｉとする。但し、iは、植生部分の輪郭を構成する画素の総数をＮとした場合における、１～Ｎまでの任意の自然数である。誤差ｅ_ｉは、下記（式１）で導出される。

ｅ_ｉ＝｜ａｘ_ｉ、ｂｙ_ｉ｜／√（ａ^２＋ｂ^２） …（式１）

個体位置決定部１４は、植生部分を構成する全ての画素における誤差ｅ_ｉの和が最小となる条件を満たすａ、及びｂの組合せを導出する。この条件を満たすａ、及びｂの組合せは、下記（式２）で表現される。（式２）におけるα、及びβは、誤差ｅ_ｉの和が最小となる条件を満たすａ、及びｂの組合せである。また、（式２）におけるΣで加算する範囲は、（ｉ＝１）から（ｉ＝Ｎ）まで、である。Ｎは、植生部分の輪郭を構成する画素の総数である。（式２）を変形すると、元のデータ（植生部分の輪郭を構成する画素）について分散共分散の固有ベクトル（α、β）を求める問題となる。

α、β＝ａｒｇｍｉｎ｛Σｅ_ｉ ^２｝ …（式２）

個体位置決定部１４は、（式２）を満たすα、βで定義される直線ＴＳ（αｘ＋βｙ＝０）を主軸にする。

或いは、個体位置決定部１４は、植生部分を貫通する直線のうち、植生部分を通る距離が長いものを主軸としてもよい。この場合、個体位置決定部１４は、植生部分の輪郭を構成する画素から抽出した任意の２点を結ぶ直線のうち、直線の距離が最大となるものを主軸とする。

或いは、個体位置決定部１４は、植生部分の輪郭を構成する画素の回帰直線を主軸とすることが考えられる。上述した主成分分析における主軸と、回帰直線とは、直線までの距離の和が最小となる直線ＴＳを求めるという点において概念が共通するためである。

しかしながら、植生部分の主軸を求める手法に回帰直線を採用することは困難である。主成分分析の主軸と回帰直線では、誤差ｅ_ｉの取り方が異なるためである。主成分分析では、画素から直線までの最短距離を誤差と定義する。一方、回帰直線では、画素から直線までのｙ軸方向の距離を誤差と定義している。このため、植生部分が、ｙ軸方向に縦長の形状である場合、誤差の和を最小とする直線の傾き（ａ／ｂの絶対値）が大きくなってしまう傾向になり、意図しない直線が主軸として導出されてしまう可能性がある。このため、本実施形態において、植生部分がｙ軸方向に縦長の形状である場合、植生部分における回帰直線を主軸として採用することは好ましくない。

個体位置決定部１４は、主軸に対して水平な方向（図７の方向ＳＨ参照）に、植生部分の画素数を投影させたヒストグラム（図７のヒストグラムＨＧ参照）を作成する。個体間境界決定部１５は、作成したヒストグラムにおいて、分布を２分割する判別分析を行う。例えば、個体位置決定部１４は、ヒストグラムのうち、画素の分布数が減少から増加に転じる位置に対応する箇所を、判別分析により分布を２分割する境界（図７の境界閾値ＫＳ参照）とする。

個体位置決定部１４は、植生部分の輪郭と主軸が交差する２点（図７の端点ＥＰ１、及びＥＰ２）を、植生部分における端点とする。個体位置決定部１４は、主軸上の端点ＥＰ１から、主軸と境界閾値ＫＳの交点を結んだ直線の中間点（図７における点ＧＰ１）を、植生部分に含まれる二つの作物個体のうち、一方の作物個体の個体位置とする。また、個体位置決定部１４は、主軸上の端点ＥＰ２から、主軸と境界閾値ＫＳの交点を結んだ直線の中間点（図７における点ＧＰ２）を、植生部分に含まれる二つの作物個体のうち、他方の作物個体の個体位置とする。

植生部分に含まれる作物個体が３個体である場合にも、植生部分に含まれる作物個体が２個体である場合と同様の手法を適用して、作物個体の個体位置を決定することができる。この場合、個体位置決定部１４は、作成したヒストグラムにおいて、分布を３分割する判別分析を行う。

次に、（ケース３）植生部分に含まれる作物個体が４個体以上である場合において、個体位置決定部１４が、植生部分における作物個体の個体位置を決定する方法について説明する。

植生部分に含まれる作物個体が４個体以上である場合、植生部分における長い楕円のような形状において、３つ以上の窪んだ部分を判別することが困難となる（図８における領域Ｒ５１参照）。このため、（ケース２）と同様の方法を適用して個体位置を決定することが難しい。

ここで、一般的に、作物の生育がある程度進んだ場合に、４個体以上の作物が互いにつながって見えるようになってくると考えられる。個体位置決定部１４は、この性質を利用し、前回以前の合成画像における作物個体の位置を用いて、今回の個体位置を決定する。

具体的に、個体位置決定部１４は、４個体以上の作物がつながる以前に撮像された空撮画像に基づく合成画像に基づいて、個体作物の個体位置を決定する。この場合、植生部分には、４個体未満の作物が含まれていることから、（ケース１）又は（ケース２）を用いて個体位置を決定することができる。個体位置決定部１４は、決定した個体位置を、例えば、合成画像に対応させ記憶させておく。

個体位置決定部１４は、今回の合成画像において抽出された植生部分に、４個体以上の作物が含まれていると推定された場合、その植生部分を用いた個体位置を決定しない。代わりに、個体位置決定部１４は、前回以前の合成画像において、今回個体位置を決定しなかった植生部分に相当する領域における個体位置（図８における領域Ｒ３１、或いはＲ４１）を、当該植生部分の個体位置とする。

個体間境界決定部１５は、個体間の境界（以下、個体間境界ともいう）を決定する。個体間境界とは、作物個体と、それに隣り合う作物個体との生長領域の境界である。一般に、植生部分の面積が大きい作物個体は、植生部分の面積が小さい作物個体と比較して、より広い生長領域を有すると考えられる。この性質が反映されるように、個体間境界決定部１５は、個体位置決定部１４により決定された作物個体の個体位置、及び作物個体の被覆面積（植生部分の面積）に基づいて、個体間境界を決定する。

個体間境界決定部１５は、個体間境界を決定する対象とする、二つの作物個体（植生部分ＳＢ３、及びＳＢ４）に着目する。個体間境界決定部１５は、二つの個体位置のそれぞれを結んだ線分の中間点を原点とし、二つの個体位置のそれぞれを結んだ直線をｘ軸とする直行座標系を定義する（図９参照）。

個体間境界決定部１５は、二つの個体位置を焦点とする双曲線（ｘ^２／Ａ^２－ｙ^２／Ｂ^２＝１）を、個体間境界（図９の境界ＲＫ３４）として、係数Ａ、Ｂを以下の方法で決定する。

まず、個体間境界決定部１５は、二つの作物個体の被覆面積（植生部分の面積）を算出し、それぞれ面積ｗ１、ｗ２とする。また、個体間境界決定部１５は、二つの作物個体の位置のｘｙ座標を、それぞれ（－Ｌ、０）、（Ｌ、０）とする。また、個体間境界決定部１５は、二つの個体位置の距離２Ｌを、面積ｗ１：ｗ２の割合で分割した距離をｐ、及びｑとする（図９参照）。この場合、ｐ、及びｑは以下の（式３）により表すことができる。

ｐ＝｛２Ｌ／（ｗ１＋ｗ２）｝ｗ１、
ｑ＝｛２Ｌ／（ｗ１＋ｗ２）｝ｗ２ …（式３）

個体間境界決定部１５は、（式３）で示されるｐ、及びｑを用いて、双曲線の頂点のｘｙ座標が（ｐ－Ｌ、０）となるように、以下の（式４）を示すようにＡ、及びＢを決定する。

Ａ＝（ｐ－ｑ）／２、
Ｂ＝√（ｐｑ）
但し、ｗ１＞ｗ２の場合、ｘ＞０
ｗ１＜ｗ２の場合、ｘ＜０
ｗ１＝ｗ２の場合、個体間境界はｙ軸に一致する …（式４）

個体生長領域算出部１６は、作物個体の生長領域を算出する。個体生長領域算出部１６は、生長領域を算出する対象とする作物個体に着目する（図１０の植生部分ＳＢ５）。個体生長領域算出部１６は、着目した作物個体と隣り合う作物個体（図１０の植生部分ＳＢ７～ＳＢ１０）それぞれとの個体間境界を取得する。ここでの個体間境界は、個体間境界決定部１５によって決定される。個体生長領域算出部１６は、取得した境界によって囲まれる領域（図１０の領域ＳＲ）を、その作物個体の生長領域とする。

上述した、推定装置１の機能部（画像取得部１０、画像合成部１１、植生部分抽出部１２、個体数推定部１３、個体位置決定部１４、個体間境界決定部１５、及び個体生長領域算出部１６）は、例えば、ＣＰＵ（Central Processing Unit）などのハードウェアプロセッサがプログラム（ソフトウェア）を実行することにより実現される。これらの構成要素のうち一部または全部は、ＬＳＩ（Large Scale Integration）やＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）、ＧＰＵ（Graphics Processing Unit）などのハードウェア（回路部；circuitryを含む）によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。プログラムは、予めＨＤＤ（Hard Disk Drive）やフラッシュメモリなどの記憶装置（非一過性の記憶媒体を備える記憶装置）に格納されていてもよいし、ＤＶＤやＣＤ－ＲＯＭなどの着脱可能な記憶媒体（非一過性の記憶媒体）に格納されており、記憶媒体がドライブ装置に装着されることでインストールされてもよい。

記憶部１７は、例えば、ＨＤＤ（Hard Disk Drive）やフラッシュメモリ、ＲＡＭ（Random Access Memory）などで構成される。記憶部１７は、例えば、画像情報を記憶する。例えば、記憶部１７は、空撮画像の画像情報を記憶する。この場合、画像取得部１０は、取得した空撮画像の情報を１７に記憶させる。記憶部１７は、合成画像の画像情報を記憶するようにしてもよい。この場合、画像合成部１１は、生成した合成画像の画像情報を記憶部１７に記憶させる。

記憶部１７は、画像合成部１１により抽出された植生部分が示された合成画像の画像情報を記憶するようにしてもよい。記憶部１７は、個体数推定部１３により推定された植生部分に含まれる個体数を示す情報を記憶するようにしてもよい。記憶部１７は、個体位置決定部１４により決定された個体位置を示す情報を記憶するようにしてもよい。記憶部１７は、個体間境界決定部１５により決定された個体間境界が示された合成画像の画像情報を記憶するようにしてもよい。記憶部１７は、個体生長領域算出部１６により算出された作物個体の生長領域が示された合成画像の画像情報を記憶するようにしてもよい。記憶部１７には、他の情報、例えば、推定装置１のＣＰＵに実行させるプログラムや、プログラムが実行される際に使用される各種パラメータなども記憶される。

図２は、実施形態に係る植生部分抽出部１２が行う処理を説明する図である。図２の左側には合成画像の一部を示す画像Ｇ１が示されている。図２の右側には画像Ｇ１から抽出した植生部分を白色、植生部分とは異なる部分を黒色で示した二値の画像Ｇ２が示されている。このように、植生部分抽出部１２は、画像Ｇ１における画素ごとの色を、植生部分の色に基づいて二値化することにより、植生部分を抽出する。なお、画像Ｇ２における領域Ｒ２０には、二つの個体作物が含まれる植生部分の例が示されている。また、画像Ｇ２における領域Ｒ２１には、１つの個体作物からなる植生部分の例が示されている。

図３～図５は、実施形態に係る個体数推定部１３が行う処理を説明する図である。図３から図５には、植生部分の形状などの特徴に基づいてサポートベクタマシンにより生成された多次元のベクトル空間が示されている。それぞれのベクトル空間には、合成画像から抽出された植生部分のそれぞれの特徴が、ベクトル空間における特徴軸に応じた位置座標で示されている。

図３には、植生部分に含まれる個体数が１つであるクラスＣ１と、二つ以上であるクラスＣ２とに分離された様子が示されている。サポートベクタマシンは、合成画像から抽出された植生部分の全てを対象として入力された説明変数に基づいて、例えば、図３に示すような、クラスＣ１とＣ２との分離曲面を提示する。

図４には、植生部分に含まれる個体数が二つ以上であるクラスＣ２が、個体数が二つであるクラスＣ２１と、３つ以上であるクラスＣ２２とに分離された様子が示されている。サポートベクタマシンは、クラスＣ２に分類された植生部分を対象として入力された説明変数に基づいて、例えば、図４に示すような、クラスＣ２１とＣ２２との分離曲面を提示する。

図５には、植生部分に含まれる個体数が３つ以上であるクラスＣ２２が、個体数が３つであるクラスＣ２２１と、４つ以上であるクラスＣ２２２とに分離された様子が示されている。サポートベクタマシンは、クラスＣ２２に分類された植生部分を対象として入力された説明変数に基づいて、例えば、図５に示すような、クラスＣ２２１とＣ２２２との分離曲面を提示する。

図６～図８は、実施形態に係る個体位置決定部１４が行う処理を説明する図である。図６には、植生部分に含まれる作物の個体数が１つである場合における個体位置、つまり（ケース１）の手法により決定された個体位置が示されている。図６における画像Ｇ２は図２における画像Ｇ２と同様に、植生部分が白色に変換された合成画像である。画像Ｇ２０は、画像Ｇ２における植生部分（１つの作物個体が含まれる植生部分）を拡大させたものである。この場合、個体位置決定部１４は、（ケース１）の手法に基づいて、植生部分ＳＢ１の重心ＧＰ１の位置を、個体位置とする。

図７には、植生部分に含まれる作物の個体数が２個体又は３個体である場合における個体位置、つまり（ケース２）の手法により決定された個体位置が示されている。図７における画像Ｇ２は図２における画像Ｇ２と同様に、植生部分が白色に変換された合成画像である。画像Ｇ２１は、画像Ｇ２における植生部分（２つつの作物個体が含まれる植生部分）を拡大させたものである。この場合、個体位置決定部１４は、（ケース２）の手法に基づいて、植生部分ＳＢ２の主軸ＳＧを決定する。個体位置決定部１４は、決定した主軸ＳＧに垂直な方向ＳＨに沿って、植生部分ＳＢを構成する画素のヒストグラムＨＧを作成する。個体位置決定部１４は、作成したヒストグラムに基づいて、境界閾値ＫＳの位置を決定する。個体位置決定部１４は、植生部分ＳＢ２の端点ＥＰ１から、主軸と境界閾値ＫＳの交点までの距離を等分する点を、一方の個体位置ＧＰ１とする。個体位置決定部１４は、植生部分ＳＢ２の端点ＥＰ２から、主軸と境界閾値ＫＳの交点までの距離を等分する点を、他方の個体位置ＧＰ２とする。

図８には、植生部分に含まれる作物の個体数が４個体以上である場合における個体位置、つまり（ケース３）の手法により決定された個体位置が示されている。図８には、撮像日が異なる３つの合成画像である画像Ｇ３～Ｇ５が示されている。画像Ｇ３～Ｇ５は、いずれも圃場における同じ領域が示されている。画像Ｇ３～Ｇ５のうち、画像Ｇ３が最も早い日付（例えば、６月７日）、画像Ｇ４が次点の日付（例えば、６月９日）、画像Ｇ５が最も遅い日付（例えば、６月１１日）に撮像されたものである。

図８の領域Ｒ５１のように、４つ以上の作物がつながった植生部分が観察された場合、個体位置決定部１４は、このままでは領域Ｒ５１の個体位置を決定することができない。そこで、個体位置決定部１４は、２日前に撮像された画像Ｇ４の領域Ｒ４１、或いは、４日前に撮像された画像Ｇ３の領域Ｒ３１で決定された個体位置を、画像Ｇ５における領域Ｒ５１の個体位置とする。

（ケース３）における個体位置の決定をより安定して行うためには、より頻繁に（例えば、毎日、或いは数日おき）空撮画像を撮像することが考えられる。作物は固定され、生育の過程で個体位置が動くことがない。このため、作物の生育期間中に少なくとも１回個体位置を決定すれば、その後日以降の合成画像から個体位置を決定する必要がないとも考えられる。

しかしながら、毎回、個体位置を決定することにより、個体位置の精度を向上させることが可能である。例えば、毎回決定した個体位置について、複数回（例えば、３回）の個体位置がともに所定の範囲の領域内であれば、その個体位置は正確な位置を示しているとみなすことが可能である。

これを利用して、例えば、個体位置決定部１４は、毎回決定した個体位置について、複数回（例えば、３回）における個体位置のそれぞれのばらつき度合いに基づき、個体位置の尤度（確からしさ）を判定するようにしてもよい。

図９は、実施形態に係る個体間境界決定部１５が行う処理を説明する図である。図９には、２つの植生部分ＳＢ３、ＳＢ４における個体間境界である境界ＲＫ３４が示されている。図９に示すように、個体間境界決定部１５は、２つの植生部分ＳＢ３、ＳＢ４の個体位置（－Ｌ、０）、（Ｌ、０）を焦点とし、植生部分ＳＢ３、ＳＢ４間の距離２Ｌを、植生部分ＳＢ３、ＳＢ４の面積ｗ１、ｗ２に応じた割合（ｗ１：ｗ２）で分割する点（ｐ－Ｌ、０）が頂点となる双曲線を、植生部分ＳＢ３とＳＢ４との個体間境界である境界ＲＫ３４とする。

図１０から図１６は、実施形態に係る個体生長領域算出部１６が行う処理を説明する図である。図１０には、作物個体の植生部分ＳＢ５の生長領域ＳＲが示されている。図１０に示すように、個体生長領域算出部１６は、植生部分ＳＢ５と、その隣り合う作物個体の植生部分ＳＢ７～ＳＢ１０のそれぞれとの間の境界を、個体間境界決定部１５から取得する。植生部分ＳＢ５とＳＢ６との個体間境界は境界ＲＫ５６である。植生部分ＳＢ５とＳＢ７との個体間境界は境界ＲＫ５７である。植生部分ＳＢ５とＳＢ８との個体間境界は境界ＲＫ５８である。植生部分ＳＢ５とＳＢ９との個体間境界は境界ＲＫ５９である。植生部分ＳＢ５とＳＢ１０との個体間境界は境界ＲＫ５１０である。個体生長領域算出部１６は、植生部分ＳＢ５と、その隣り合う作物個体の植生部分ＳＢ７～ＳＢ１０のそれぞれとの間の境界（境界ＲＫ５６、ＲＫ５７、ＲＫ５８、ＲＫ５９、及びＲＫ５１０）で囲まれる領域を、生長領域ＳＲとする。

図１１には、画像Ｇ１０において抽出された植生部分のそれぞれに含まれる作物個体ごとの生長領域が示されている。図１１に示すように、圃場の端に植えられた作物については、隣り合う作物個体が存在しない側の境界として、圃場と圃場外との境界を適用するようにしてもよい。

図１２には、作物の生育の過程により、その生長領域が変化していく様子が示されている。隣り合う作物の植生部分の面積が同程度である場合、境界は、個体位置を結んだ直線の中心を通り、個体位置を結んだ直線にほぼ垂直な線となる（撮影日１参照）。作物の生長が進み、隣り合う作物の植生部分の面積に差が出てくると、その個体間境界は、個体位置を結んだ直線を、面積に比率で分割した点（面積の小さい方の側によった位置）を頂点とする双曲線となる（撮影日２参照）。隣り合う作物の植生部分の面積の差が大きくなると、その個体間境界は、双曲線の曲率が更に大きくしたものとなる（撮影日３参照）。隣り合う作物の植生部分がつながった場合、植生部分がつながる前に決定した境界を、最終的な境界として採用する（撮影日４参照）。

図１３～図１６には、圃場における作物が発芽し始める段階から、植生部分がつながるまでの各段階において決定された個体間境界が示されている。作物が発芽し始める段階においては、空撮画像の撮像が行われる度に新しい個体位置が検出されるため、新たに検出された個体位置に基づいて新たな生長領域が決定される（図１３、及び図１４参照）。個体間境界は、植生部分の面積を重みとして決定されるため、作物が生長する段階では、撮像が行われる度に生長領域の形状が変化する（図１４、及び１５参照）。生長領域の形状は、隣り合う作物の植生部分がつながるまで更新される（図１５、及び１６参照）。

図１７は、実施形態に係る推定装置１が行う処理の流れを示すフロー図である。推定装置１は、ドローン等によって上空から圃場が撮像された空撮画像を取得する（ステップＳ１）。推定装置１は、取得した空撮画像をつなぎ合わせて合成し、圃場全体の合成画像（オルソ画像）を生成する（ステップＳ２）。推定装置１は、合成画像の色を二値化するなどして、合成画像における植生部分を抽出する（ステップＳ３）。推定装置１は、抽出した植生部分のそれぞれに含まれる作物個体の個体数を、サポートベクタマシンを利用するなどして推定する（ステップＳ４）。推定装置１は、推定した個体数に応じて、主成分分析を利用するなどして、作物個体の個体位置を決定する（ステップＳ５）。推定装置１は、決定した個体位置を用いて、隣り合う作物個体の被覆面積（植生部分の面積）を重みづけした境界（個体間境界）を決定する（ステップＳ６）。推定装置１は、作物個体に隣り合う作物との間で決定した境界により囲まれる領域を、その作物個体の生長領域として算出する（ステップＳ７）。

以上、説明した通り、実施形態の推定装置１は、個体間境界決定部１５と、個体生長領域算出部１６とを備える。個体間境界決定部１５は、圃場が撮像された空撮画像に基づいて、圃場において隣り合う作物の位置及び被覆面積に基づいて、隣り合う作物における生長領域の境界である個体間境界を決定する。圃場は「対象場所」の一例である。空撮画像は、「対象画像」の一例である。個体生長領域算出部１６は、個体間境界決定部１５によって判定された個体間境界を用いて、圃場に植生している作物のそれぞれの生長領域を算出する。

これにより、実施形態の推定装置１では、隣り合う作物個体のそれぞれの個体位置及び被覆面積に基づいて、個体間境界を決定する。このため、作物の個体の大きさを考慮した生長領域を推定することができる。したがって、圃場に植えられた作物が格子状に配置されていない場合であっても、作物の個体の大きさを考慮した生長領域を推定することが可能である。

また、実施形態の推定装置１では、画像取得部１０と、植生部分抽出部１２と、個体数推定部１３と、個体位置決定部１４とを更に備える。画像取得部１０は、圃場が空撮された空撮画像を取得する。植生部分抽出部１２は、画像取得部１０によって取得された空撮画像において作物が植生している領域である植生部分を抽出する。個体数推定部１３は、植生部分抽出部１２によって抽出された植生部分に含まれる作物の個体数を推定する。個体位置決定部１４は、個体数推定部１３によって推定された、植生部分に含まれる作物個体の個体数に基づいて、植生部分に含まれる個体作物の位置である個体位置を決定する。個体間境界決定部１５は、個体位置決定部１４によって決定された個体位置を、圃場において隣り合う作物の位置とする。

これにより、実施形態の推定装置１では、画像から抽出した植生部分に含まれる作物の個体数に応じて、個体位置を推定する。このため、作物が格子状に配置されていない場合であっても、個体位置を決定できる。このようにして決定された個体位置を用いて、個体間境界を決定することができるため、作物の個体の大きさを考慮した生長領域をより精度よく推定することができる。

また、実施形態の推定装置１では、個体数推定部１３は、植生部分における面積、及び形状に関する特徴に基づいた分類を行うことにより、植生部分における個体数を推定する。これにより、実施形態の推定装置１では、植生部分の形状がもつ特徴に基づいて個体数を推定することができる。一般に、植生部分は、植生部分に含まれる作物個体の個体数によってその形状が互いに異なる。したがって、植生部分の形状がもつ特徴に基づいて個体数を推定すれば、植生部分に含まれる作物個体の個体数を精度よく推定することが可能となる。

また、実施形態の推定装置１では、個体数推定部１３は、植生部分における形状に関する特徴として、当該形状における円形度、及び最小外接矩形の長辺長と短辺長とを用いる。これにより、実施形態の推定装置１では、植生部分の形状に係る円形度、及び最小外接矩形の長辺長と短辺長に基づいて個体数を推定することができる。一般に、作物は、個体位置を中心とする略円形に育成していく。このため、植生部分に含まれる個体数が１つであればその形状は略円形に近い形状となり、複数の作物が含まれる場合には楕円に近い形状となると考えられる。したがって、植生部分の形状の円形度、及び最小外接矩形の長辺長と短辺長の特徴に基づいて個体数を推定すれば、植生部分に含まれる作物個体の個体数を精度よく推定することが可能となる。

また、実施形態の推定装置１では、個体数推定部１３は、サポートベクタマシンに、植生部分における形状に関する特徴を入力して得られる分類結果を用いて、植生部分における個体数を推定する。これにより、実施形態の推定装置１では、サポートベクタマシンを用いて個体数を推定することができる。したがって、サポートベクタマシンに、形状に関する特徴を入力するという簡単な方法により、植生部分に含まれる作物個体の個体数を精度よく推定することが可能となる。

また、実施形態の推定装置１では、個体位置決定部１４は、個体数推定部１３により植生部分に複数の個体が含まれると推定された場合、当該植生部分における主軸に垂直な方向に存在する植生部分の画素数の分布に基づいて個体位置を決定する。これにより、実施形態の推定装置１では、主軸、及び主軸に垂直な方向に存在する画素数の分布に応じて個体位置を決定することができる。植生部分に複数の作物個体が含まれる場合、互いの作物がつながってできる形状が、楕円の中央部分が窪んだ、ひょうたんのような形状となる。したがって、主軸に垂直な方向に存在する画素数の分布に応じて、この窪んだ部分を検出することができ、検出した部分を境界とみなすことで、個体位置を精度よく推定することが可能となる。

また、実施形態の推定装置１では、個体位置決定部１４は、植生部分に含まれる画素の分布を主成分分析して得られる第１主成分軸を、植生部分における主軸とする。これにより、実施形態の推定装置１では、植生部分を構成する画素群の分布のうち、最もばらつきが大きい（分散が最大となる）軸を主軸とすることができる。植生部分に複数の作物個体が含まれる場合、植生部分の形状は、個体位置を結んだ直線に沿った方向に幅広となると考えられる。したがって、主軸を主成分分析の第１主成分軸とすることにより、個体位置を通る可能性が高い直線を検出でき、この直線を用いて個体位置を精度よく推定することが可能となる。

また、実施形態の推定装置１では、個体位置決定部１４は、植生部分を貫通する直線のうち、植生部分を貫通する距離が最も長い直線を、植生部分における主軸とする。これにより、実施形態の推定装置１では、上述した効果と同様の効果を奏することができる。

また、実施形態の推定装置１では、個体間境界決定部１５は、隣り合う２つの作物個体の個体位置を焦点とし、当該２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を頂点とする双曲線を、個体間境界とする。これにより、実施形態の推定装置１では、隣り合う２つの作物個体の面積に応じて、面積がより大きい作物個体がより広い生長領域となり、面積がより小さい作物個体がより狭い生長領域となるように、個体間境界を決定することができる。したがって、作物個体が利用できる生長領域の実態が反映された個体間境界を決定することが可能である。

なお、上述した実施形態では、個体間境界を双曲線とする場合を例に説明した。しかしながらこれに限定されることはない。個体間境界は、少なくとも、隣り合う２つの作物個体の面積に応じて決定されればよい。例えば、個体間境界決定部１５は、隣り合う２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を通る線（直線又は曲線）を、個体間境界とするようにしてもよい。これにより、実施形態の推定装置１では、隣り合う２つの作物個体の面積に応じた個体間境界を決定することができる。したがって、上述した効果と同様の効果を奏することができる。

また、上述した実施形態では、空撮画像を用いて圃場に植生する作物個体の生長領域を推定する場合を例に説明した。しかしながら空撮画像に限定されることはない。推定装置１は、圃場に植生する作物における位置と大きさとが少なくとも判る程度に撮像されたデジタル画像であれば、任意の画像を用いて作物個体の生長領域を推定してよい。

また、上述した実施形態では、圃場に植生する作物個体の生長領域を推定する場合を例に説明した。しかしながら圃場に限定されることはない。推定装置１は、少なくとも作物が植生する場所であれば、任意の場所における作物個体の生長領域を推定してよい。

なお、推定装置１の機能部が行う処理（演算及び制御に関する処理）の全部または一部の処理を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することで各部の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。

以上、本発明の実施形態を、図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計変更等も含まれる。

１ 推定装置
１０ 画像取得部
１１ 画像合成部
１２ 植生部分抽出部
１３ 個体数推定部
１４ 個体位置決定部
１５ 個体間境界決定部
１６ 個体生長領域算出部
１７ 記憶部

Claims (12)

1. 対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定する個体間境界決定部と、
   前記個体間境界決定部によって決定された前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出する生長領域算出部と、
   を備える推定装置。
2. 前記対象画像を取得する画像取得部と、
   前記画像取得部によって取得された前記対象画像に基づいて、前記対象場所において作物が植生している領域である植生部分を抽出する植生部分抽出部と、
   前記植生部分抽出部によって抽出された前記植生部分に含まれる作物個体の個体数を推定する個体数推定部と、
   前記個体数推定部によって推定された前記植生部分に含まれる作物個体の個体数に基づいて、前記植生部分に含まれる作物個体の位置である個体位置を決定する個体位置決定部と、
   を更に備え、
   前記個体間境界決定部は、前記個体位置決定部によって決定された前記個体位置を、前記対象場所において隣り合う作物のそれぞれの位置とする、
   請求項１に記載の推定装置。
3. 前記個体数推定部は、前記植生部分における面積、及び形状に関する特徴に基づいた分類を行うことにより、前記植生部分における個体数を推定する、
   請求項２に記載の推定装置。
4. 前記個体数推定部は、前記植生部分における形状に関する特徴として、当該形状における円形度、及び最小外接矩形の長辺長と短辺長とを用いる、
   請求項３に記載の推定装置。
5. 前記個体数推定部は、サポートベクタマシンに、前記植生部分における形状に関する特徴を入力して得られる分類結果を用いて、前記植生部分における個体数を推定する、
   請求項２から請求項４のいずれか一項に記載の推定装置。
6. 前記個体位置決定部は、前記個体数推定部により前記植生部分に複数の個体が含まれると推定された場合、当該植生部分における主軸に垂直な方向に存在する植生部分の画素数の分布に基づいて前記個体位置を決定する、
   請求項２から請求項５のいずれか一項に記載の推定装置。
7. 前記個体位置決定部は、前記植生部分に含まれる画素の分布を主成分分析して得られる第１主成分軸を、前記植生部分における主軸とする、
   請求項６に記載の推定装置。
8. 前記個体位置決定部は、前記植生部分を貫通する直線のうち、前記植生部分を貫通する距離が最も長い直線を、前記植生部分における主軸とする、
   請求項６に記載の推定装置。
9. 前記個体間境界決定部は、隣り合う２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を通る線を、前記個体間境界とする、
   請求項１から請求項８の何れか一項に記載の推定装置。
10. 前記個体間境界決定部は、隣り合う２つの作物個体の個体位置を焦点とし、当該２つの作物個体の個体位置を結ぶ線分を、当該２つの作物個体の面積の比率に応じて分割する点を頂点とする双曲線を、前記個体間境界とする、
    請求項１から請求項９の何れか一項に記載の推定装置。
11. 個体間境界決定部が、対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定し、
    生長領域算出部が、前記個体間境界決定部によって決定された前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出する、
    推定方法。
12. コンピュータに、
    対象場所が撮像された対象画像に基づいて、前記対象場所において隣り合う作物のそれぞれの位置及び被覆面積に基づいて、前記隣り合う作物における生長領域の境界である個体間境界を決定させ、
    前記個体間境界を用いて、前記対象場所に植生している作物のそれぞれの生長領域を算出させる、
    ためのプログラム。

Images