WO2017221641A1

WO2017221641A1 - 植物生育指標測定装置、該方法および該プログラム

Info

Publication number: WO2017221641A1
Application number: PCT/JP2017/020022
Authority: WO
Inventors: 片桐　哲也; 片桐　禎人
Original assignee: コニカミノルタ株式会社
Priority date: 2016-06-22
Filing date: 2017-05-30
Publication date: 2017-12-28
Also published as: JP6981408B2; JPWO2017221641A1

Abstract

本発明は、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像を生成する植物生育指標測定装置、該方法および該プログラムである。本発明は、所定の時間間隔で画像を生成し、前記生育指標を求めるために有意な画像のみを画像記憶部に記憶するように、前記生成された画像を画像記憶部に記憶するか否かを判定し、その結果、前記画像記憶部に記憶すると判定された場合のみ、前記生成された画像を画像記憶部に記憶する。

Description

植物生育指標測定装置、該方法および該プログラム

　本発明は、植物における生育の度合いを表す生育指標を求める植物生育指標測定装置、植物生育指標測定方法および植物生育指標測定プログラムに関する。

　農業では、高品質および安定多収穫な農作物の植物を育てるために、例えば追肥時期や追肥量等の施肥管理を適切に実施する必要がある。そのために、現状の植物の状態が判定される。この判定には、例えばＮＤＶＩ（Ｎｏｒｍａｌｉｚｅｄ　Ｄｉｆｆｅｒｅｎｃｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、正規化植生指標）、ＤＶＩ（Ｄｉｆｆｅｒｅｎｃｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、差植生指標）およびＲＶＩ（Ｒａｔｉｏ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、比植生指標）等の、植物における生育の度合いを表す生育指標が用いられる。このような生育指標は、圃場を、互いに異なる波長帯域で撮像した複数の画像に基づいて求められる。ここで、圃場は、通常、前記画像を生成するカメラの画角に対して広いため、圃場全体を１回の撮像による１枚の画像で写すことが難しい。このため、例えば、特許文献１では、互いに隣り合う圃場をオーバーラップして複数の画像を取得し、同一位置の圃場で色を補正し、この補正した複数の画像を１枚に合成し、これによって圃場全体の１枚の画像を得て、圃場の画像のＲＧＢ値から窒素相対分布を求める作物の生育診断方法が提案されている（例えば特許文献１の図１１）。

　ところで、作物の生育状況を適切に把握するためには、適度な画像解像度の画像で、生育指標を求める必要がある。例えば、６４０ドット×４８０ドットであるＶＧＡクラスの画像解像度を持つ撮像装置で高さ２０ｍから圃場を撮像することが好ましい。この場合、１回の撮像による１枚の画像には、１６．６ｍ×１２．４ｍの範囲で圃場が写り込み、複数の画像を１枚の画像に合成する際に、互いに隣接する画像間で、圃場の同じ領域を写し込んだ重複部分が存在していることが好ましいため、例えば、３０ｍ×１００ｍの圃場全体の生育指標を求めるために、約４０枚の画像が必要となり、したがって、その画像の記録（保存、記憶）に比較的大容量の記憶装置が必要となる。

　一方、圃場の個々の領域それぞれを写した複数の画像を得るために、前記撮像装置を搭載した、例えばヘリコプター等の航空機を圃場上空で例えば九十九折りで往復移動させる必要がある。前記撮像装置が所定の時間間隔で自動的に撮像および記録する場合、折り返し場所の画像は、圃場の画像ではないため、生育指標の演算には、不必要な画像となり、したがって、無駄な記録が生じてしまう。また、前記航空機を無線で遠隔操縦する場合、その未熟な操縦技術によっても、無駄な記録が生じてしまう虞がある。

　そして、このような生育指標の演算には不必要な画像が複数の画像に含まれると、これら複数の画像を合成する際に、無駄な合成の処理が生じ、あるいは、不必要な画像を除去する処理が必要となり、したがって、合成の処理の効率が低下してしまう。

特開２００３－３３９２３８号公報

　本発明は、上述の事情に鑑みて為された発明であり、その目的は、記憶装置の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる植物生育指標測定装置および植物生育指標測定方法、植物生育指標測定プログラムを提供することである。

　上記並びにその他の本発明の目的、特徴及び利点は、以下の詳細な記載と添付図面から明らかになるであろう。

第１ないし第３実施形態にかかる植物生育指標測定装置の構成を示す図である。圃場に生育する作物の条または条間に応じた周波数に基づく画像記憶の判定処理を説明するための図である。圃場の生育指標に基づく画像記憶の判定処理を説明するための図である。太陽の写り込みの検出に基づく画像記憶の判定処理を説明するための図である。第１実施形態における植物生育指標測定装置の動作を示すフローチャートである。第１実施形態における植物生育指標測定装置において、その飛行軌跡と画像記憶位置との関係を説明するための図である。図６Ａに示す飛行軌跡において、植物生育指標測定装置によって撮像され生成された画像を説明するための図である。第２実施形態における植物生育指標測定装置の動作を示すフローチャートである。第３実施形態における植物生育指標測定装置の動作を示すフローチャートである。第４実施形態における植物生育指標測定システムの構成を示す図である。前記植物生育指標測定システムに用いられる画像生成装置の動作を示すフローチャートである。前記植物生育指標測定システムに用いられる植物生育指標測定装置の動作を示すフローチャートである。図１２に示すフローチャートにおける対応点探索の処理を説明するための図である。図１２に示すフローチャートにおける画像変換の処理を説明するための図である。ＮＤＶＩ合成画像の一例を示す図である。低コントラストによる合成精度の低下を説明するための図である。露出のばらつきによる合成精度の低下を説明するための図である。第１変形形態におけるパノラマ合成処理を説明するための図である。第２変形形態におけるパノラマ合成処理を説明するための図である。

　以下、本発明にかかる実施の一形態を図面に基づいて説明する。なお、各図において同一の符号を付した構成は、同一の構成であることを示し、適宜、その説明を省略する。本明細書において、総称する場合には添え字を省略した参照符号で示し、個別の構成を指す場合には添え字を付した参照符号で示す。

　実施形態における植物生育指標測定装置は、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成する植物生育指標測定装置であって、所定の時間間隔で画像を生成する画像生成部と、前記画像を記憶する画像記憶部と、前記生育指標を求めるために有意な画像のみを前記画像記憶部に記憶するように、前記画像生成部で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定部と、前記画像判定部で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成部で生成された画像を前記画像記憶部に記憶する画像記憶制御部とを備えるものである。以下、このような植物生育指標測定装置について、第１ないし第３実施形態で、より具体的に説明する。

　（第１実施形態）
　図１は、第１ないし第３実施形態にかかる植物生育指標測定装置の構成を示す図である。なお、図１は、第１ないし第３実施形態における各植物生育指標測定装置Ｍａ、Ｍｂ、Ｍｃの各構成を示し、図１において、添え字ａは、その構成が第１実施形態における植物生育指標測定装置Ｍａに関する構成であることを示し、添え字ｂは、その構成が第２実施形態における植物生育指標測定装置Ｍｂに関する構成であることを示し、そして、添え字ｃは、その構成が第３実施形態における植物生育指標測定装置Ｍｃに関する構成であることを示す。図２は、圃場に生育する作物の条または条間に応じた周波数に基づく画像記憶の判定処理を説明するための図である。図２Ａは、一例として、植物生育指標測定装置Ｍａが比較的低い高度で圃場を撮像した画像を示し、図２Ｂは、一例として、植物生育指標測定装置Ｍａが比較的高い高度で圃場を撮像した画像を示す。図３は、圃場の生育指標に基づく画像記憶の判定処理を説明するための図である。図３Ａは、一例として、圃場の可視画像を示し、図３Ｂは、そのＮＤＶＩ画像を示す。図３Ｃは、一例として、畦等の圃場外の可視画像を示し、図３Ｄは、そのＮＤＶＩ画像を示す。図４は、太陽の写り込みの検出に基づく画像記憶の判定処理を説明するための図である。

　第１実施形態における植物生育指標測定装置Ｍａは、例えば、図１に示すように、航空機１と、自機状態測定部２ａと、制御処理部３ａと、記憶部４ａと、画像生成部５と、インターフェース部（ＩＦ部）６とを備える。

　航空機１は、大気中を飛行する装置であり、例えば、気球、飛行船、飛行機、ヘリコプタおよびマルチコプタ等である。航空機１は、有人機であって良いが、好ましくは、無線操縦飛行（誘導飛行）または自律飛行による無人機（ドローン）である。航空機１は、本実施形態では、制御処理部３ａに接続され、誘導飛行または自律飛行による制御処理部３ａの制御に従って飛行する。

　自機状態測定部２ａは、制御処理部３ａに接続され、制御処理部３ａの制御に従って自機の飛行状態を測定する装置である。第１実施形態では、自機状態測定部２ａは、自機の位置（緯度Ｘ、経度Ｙ）および高度Ｚを測定する。このために、自機状態測定部２ａは、地球上における現在の自機の位置（緯度Ｘ、経度Ｙ）および高度Ｚを測定する位置高度測定部２１を備える。位置高度測定部２１は、例えば、衛星測位システムのＧＰＳ（Ｇｌｏｂａｌ　Ｐｏｓｉｔｉｏｎｉｎｇ　Ｓｙｓｔｅｍ）を備えて構成される。なお、ＧＰＳは、ＤＧＳＰ（Ｄｉｆｆｅｒｅｎｔｉａｌ　ＧＳＰ）等の誤差を補正する補正機能を持ったＧＰＳであっても良い。位置高度測定部２１は、その測定結果（自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を制御処理部３ａへ出力する。前記自機の測定位置および測定高度は、後述のように生成される複数の画像、に基づいて求められる複数の生育指標画像、を合成して１枚の生育指標合成画像を生成する際に用いられ、前記複数の生育指標画像それぞれの位置およびその撮影高度を示す。生育指標画像は、生育指標を表す画素値を持つ画素で構成された画像である。

　画像生成部５は、制御処理部３ａに接続され、制御処理部３ａの制御に従って、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像を、所定の時間間隔で生成する装置である。画像生成部５は、その生成した画像を制御処理部３ａへ出力する。前記生育指標は、例えば、ＮＤＶＩ（Ｎｏｒｍａｌｉｚｅｄ　Ｄｉｆｆｅｒｅｎｃｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、正規化植生指標、ＮＤＶＩ＝（Ｉ－Ｒ）／（Ｉ＋Ｒ））、ＲＶＩ（Ｒａｔｉｏ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、比植生指標、ＲＶＩ＝Ｒｉ’／Ｒｖ）、ＤＶＩ（Ｄｉｆｆｅｒｅｎｃｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、差植生指標、ＤＶＩ＝Ｒｉ’－Ｒｖ）、ＴＶＩ（Ｔｒａｎｓｆｏｒｍｅｄ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、ＴＶＩ＝ＮＤＶＩ＋０．５）^０．５）、および、ＩＰＶＩ（Ｉｎｆｒａｒｅｄ　Ｐｅｒｃｅｎｔａｇｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ、ＩＰＶＩ＝Ｒｉ’／（Ｒｉ’＋Ｒｖ）＝（ＮＤＶＩ＋１）／２）等である。画像生成部５は、これら各種の生育指標のうち、採用された生育指標に応じた適宜な波長帯域の画像を生成する。本実施形態では、一例として、前記生育指標は、ＮＤＶＩである。このため、画像生成部５は、前記画像として、可視画像および近赤外画像を生成する。より具体的には、画像生成部５は、可視撮像部５１と、赤外撮像部５２とを備える。

　可視撮像部５１は、可視光の画像（可視画像）を生成する装置である。より具体的には、本実施形態では、ＮＤＶＩを求めるために、可視撮像部５１は、例えば、波長６５０ｎｍを中心波長とする比較的狭帯域で光を透過する第１バンドパスフィルタ、前記第１バンドパスフィルタを透過した測定対象の可視光の光学像を所定の第１結像面上に結像する第１結像光学系、前記第１結像面に受光面を一致させて配置され、前記測定対象の可視光の光学像を電気的な信号に変換する第１イメージセンサ、前記第１イメージセンサの出力に対し公知の画像処理を施して可視画像を生成する第１デジタルシグナルプロセッサ（ＤＳＰ）等を備えて構成される、いわゆるカメラ等である。可視撮像部５１は、この撮像して生成した可視画像を制御処理部３ａへ出力する。

　赤外撮像部５２は、赤外光（例えば近赤外光）の画像（赤外画像（近赤外画像））を生成する装置である。より具体的には、本実施形態では、ＮＤＶＩを求めるために、赤外撮像部５２は、例えば、波長８００ｎｍを中心波長とする比較的狭帯域で光を透過する第２バンドパスフィルタ、前記第２バンドパスフィルタを透過した測定対象の赤外光の光学像を所定の第２結像面上に結像する第２結像光学系、前記第２結像面に受光面を一致させて配置され、前記測定対象の赤外光の光学像を電気的な信号に変換する第２イメージセンサ、前記第２イメージセンサの出力に対し公知の画像処理を施して近赤外画像を生成する第２ＤＳＰ等を備えて構成される、いわゆる赤外カメラ等である。赤外撮像部５２は、この撮像して生成した近赤外画像を制御処理部３ａへ出力する。

　ＩＦ部６は、制御処理部３ａに接続され、制御処理部３ａの制御に従って、例えば外部機器との間でデータの入出力を行う回路であり、例えば、シリアル通信方式であるＲＳ２３２Ｃのインターフェース回路、Ｂｌｕｅｔｏｏｔｈ（登録商標）規格を用いたインターフェース回路、ＩｒＤＡ（Ｉｎｆｒａｒｅｄ　Ｄａｔａ　Ａｓｓｃｏｉａｔｉｏｎ）規格等の赤外線通信を行うインターフェース回路、および、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ　Ｓｅｒｉａｌ　Ｂｕｓ）規格を用いたインターフェース回路等である。また、ＩＦ部６は、有線または無線によって通信する通信カード等であり、例えばイーサネット環境等の通信ネットワークを介して例えば外部機器との間で通信しても良い（イーサネットは登録商標）。

　記憶部４ａは、制御処理部３ａに接続され、制御処理部３ａの制御に従って、各種の所定のプログラムおよび各種の所定のデータを記憶する回路である。前記各種の所定のプログラムには、例えば、当該植物生育指標測定装置Ｍａの各部を当該各部の機能に応じて制御する制御プログラムや、画像生成部５で生成された画像に基づいて、前記圃場に生育する作物の条または条間に応じた周波数を求める画像周波数処理プログラムや、画像生成部５で生成された画像に基づいて、生育指標（本実施形態ではＮＤＶＩ）を求める生育指標処理プログラムや、画像生成部５で生成された画像に基づいて、太陽の写り込みを検出する太陽写込検出プログラムや、生育指標を求めるために有意な画像のみを後述の画像記憶部４１に記憶するように、画像生成部５で生成された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行う画像判定プログラムや、前記画像判定プログラムで画像記憶部４１に記憶すると判定された場合のみ、画像生成部５で生成された画像を画像記憶部４１に記憶する画像記憶制御プログラム等の制御処理プログラムが含まれる。前記各種の所定のデータには、画像生成部５で生成された画像（本実施形態では可視画像および近赤外画像）や、前記判定処理で用いる予め設定された各閾値等の、各プログラムを実行する上で必要なデータ等が含まれる。記憶部４ａは、例えば不揮発性の記憶素子であるＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）や書き換え可能な不揮発性の記憶素子であるＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ　Ｅｒａｓａｂｌｅ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）等を備える。そして、記憶部４ａは、前記所定のプログラムの実行中に生じるデータ等を記憶するいわゆる制御処理部３ａのワーキングメモリとなるＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）等を含む。

　記憶部４ａは、画像生成部５で生成され、後述の画像判定部３５ａで記憶すると判定された画像（本実施形態では可視画像および近赤外画像）を記憶する画像記憶部４１を機能的に備える。画像記憶部４１は、例えば、画像生成部５で画像を生成した回数を表すシリアル番号、画像生成部５で生成され、画像判定部３５ａで記憶すると判定された画像、および、当該画像を撮像する際に位置高度測定部２１で測定された自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を互いに対応付けて記憶する。そして、画像生成部５で生成された画像が、画像判定部３５ａで記憶しないと判定された場合、画像記憶部４１は、前記シリアル番号、および、画像生成部５で生成され、画像判定部３５ａで記憶しないと判定された画像を撮像する際に位置高度測定部２１で測定された自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を互いに対応付けて記憶する。この場合、前記自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）は、非記憶画像情報の一例である。前記非記憶画像情報は、画像判定部３５ａで画像記憶部４１に記憶しないと判定された場合に、画像生成部５で生成された画像であって画像判定部３５ａで画像記憶部４１に記憶しないと判定された前記画像を指し示す所定の情報である。なお、非記憶画像情報としての測定高度Ｚは、省略されても良い。

　制御処理部３ａは、植物生育指標測定装置Ｍａの各部を当該各部の機能に応じてそれぞれ制御し、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成するための回路である。制御処理部３ａは、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）およびその周辺回路を備えて構成される。制御処理部３ａは、機能的に、制御部３１と、画像周波数処理部３２と、生育指標処理部３３と、太陽写込検出部３４と、画像判定部３５ａと、画像記憶制御部３６とを備える。

　制御部３１は、当該植物生育指標測定装置Ｍａの各部を当該各部の機能に応じて制御し、当該植物生育指標測定装置Ｍａの全体制御を司るものである。

　画像周波数処理部３２は、画像生成部５で生成された画像に基づいて、圃場に生育する作物の条または条間に応じた周波数を求めるものである。より具体的には、画像周波数処理部３２は、エッジを検出するための画像フィルタを用いることによって、画像生成部５で生成された画像からエッジを検出し、このエッジを検出した画像から直線のハフ変換によって複数の直線を前記条または前記条間として求め、この求めた複数の直線の間隔から前記周波数を求める。この周波数を求める際に用いる画像は、可視画像であって良く、近赤外画像であって良く、これら可視画像および近赤外画像に基づいて、生育指標処理部３３によって求められた生育指標画像であって良い。画像周波数処理部３２は、この求めた周波数を画像判定部３５ａに通知する。

　例えば、前記圃場が水田である場合、前記作物は、稲であり、前記条は、稲の条である。植物生育指標測定装置Ｍａが比較的低い高度を飛行している場合、圃場の画像は、例えば、図２Ａに示す画像となり、その画像に基づいて求められた、作物の条または条間に応じた周波数は、比較的低い周波数となる。一方、植物生育指標測定装置Ｍａが比較的高い高度を飛行している場合、圃場の画像は、例えば、図２Ｂに示す画像となり、その画像に基づいて求められた、作物の条または条間に応じた周波数は、比較的高い周波数となる。このように、画像生成部５で生成された画像に基づいて求められた、圃場に生育する作物の条または条間に応じた周波数は、植物生育指標測定装置Ｍａの高度（飛行高度）と相関する。上述したように、作物の生育状況を適切に把握するためには、適度な画像解像度の画像で、生育指標を求める必要があり、このため、画像生成部５における可視撮像部５１および赤外撮像部５２の各画像解像度に応じた適度な高度（例えば１５ｍや２０ｍや２５ｍ等）で撮像して生成した画像である必要がある。このため、本実施形態では、画像生成部５で生成された画像に基づいて求められた、圃場に生育する作物の条または条間に応じた周波数から、植物生育指標測定装置Ｍａの高度（飛行高度）が推定され、その高度で生成された画像に対するその記憶の許否が、画像判定部３５ａによって後述のように判定されている。

　生育指標処理部３３は、画像生成部５で生成された画像に基づいて、生育指標を求めるものである。より具体的には、本実施形態では、生育指標処理部３３は、画像生成部５における可視撮像部５１および赤外撮像部５２それぞれで生成された可視画像および近赤外画像から、ＮＤＶＩを表した画素値を持つＮＤＶＩ画像を求める。ＮＤＶＩ画像の各画素値は、可視画像の画素値をＩとし、近赤外画像の画素値をＲとする場合に、ＮＤＶＩ画像の当該画素位置に対応する可視画像および近赤外画像の画素位置同士において、ＮＤＶＩ＝（Ｉ－Ｒ）／（Ｉ＋Ｒ）によって求められる。例えば、ＮＤＶＩ画像の画素位置（０、０）における画素値は、可視画像の画素位置（０、０）における画素値および近赤外画像の画素位置（０、０）における画素値から求められ、また例えば、ＮＤＶＩ画像の画素位置（１０、１５）における画素値は、可視画像の画素位置（１０、１５）における画素値および近赤外画像の画素位置（１０、１５）における画素値から求められる。なお、可視撮像部５１と、赤外撮像部５２との視差が考慮され、可視画像の画素位置および近赤外画像の画素位置の中の少なくともいずれか一方の画素位置が前記視差を修正するようにシフトされた上で、ＮＤＶＩが求められても良い。また、可視撮像部５１の前記第１結像光学系と赤外撮像部５２の前記第２結像光学系は、例えば画角や歪曲収差等の光学特性が互いに同等であるものとする。

　植物生育指標測定装置Ｍａは、生育指標を求めるための画像が必要であるため、上述したように、圃場ではない画像を記録（記憶）しても無駄になる。圃場ではない領域は、例えば畦等であり、圃場ではない前記領域には、畦道があったり、水路があったり、雑草が生育していたりする。これらの反射率は、圃場に生育する作物の反射率と異なり、通常、比較的低い。このため、例えば図３Ａに示すように、圃場を写した画像に基づくＮＤＶＩ画像は、圃場のＮＤＶＩが高いので、例えば図３Ｂに示すように、全面的に高い画素値の画像となる。一方、例えば図３Ｃに示すように、圃場ではない領域を多く含む画像に基づくＮＤＶＩ画像は、圃場ではない領域のＮＤＶＩが低いので、例えば図３Ｄに示すように、ＮＤＶＩが低い領域を多く含む画像となる。このため、本実施形態では、画像生成部５で生成された画像に基づいて求められた生育指標画像（ここではＮＤＶＩ画像）から、前記画像に圃場ではない領域が多く含まれるか否かが推定され、前記画像に対するその記憶の許否が、画像判定部３５ａによって後述のように判定されている。

　太陽写込検出部３４は、画像生成部５で生成された画像に基づいて、太陽の写り込みを検出するものである。より具体的には、太陽写込検出部３４は、例えば、画像生成部５で生成された可視画像から、画素値が飽和した画素の総数を求めることによって、画素値が飽和した画素領域のサイズ（広さ、画素数）を求め、この求めたサイズが、予め設定された所定の閾値（第３閾値）以上である画素領域を前記太陽の写り込みとして検出する。前記第３閾値は、太陽の写り込んだ複数のサンプルの画像から、適宜に設定される。

　植物生育指標測定装置Ｍａは、生育指標を求めるための画像が必要であるため、例えば、図４に示すように、太陽が写り込んだ画像を記録（記憶）しても、太陽が写り込んだ領域の画像では正しい生育指標（本実施形態ではＮＤＶＩ）が求められないから、無駄になる。このため、本実施形態では、画像生成部５で生成された画像に基づいて太陽の写り込みが検出され、その検出結果に基づいて前記画像に対するその記憶の許否が、画像判定部３５ａによって後述のように判定されている。

　画像判定部３５ａは、前記生育指標を求めるために有意な画像のみを画像記憶部４１に記憶するように、画像生成部５で生成された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行うものである。本実施形態では、画像判定部３５ａは、画像周波数処理部３２で求めた周波数に基づいて、前記判定処理を行う。より具体的には、画像判定部３５ａは、前記判定処理として、画像周波数処理部３２で求めた周波数が、予め設定された所定の撮像高度（飛行高度）で撮像した画像に写る前記条または前記条間に応じた周波数（撮像高度対応周波数）に、予め設定された所定の第１マージンの範囲内で一致するか否かを判定し、前記一致する場合を画像記憶部４１に記憶する場合と判定し、前記一致しない場合を画像記憶部４１に記憶しない場合と判定する。前記第１マージンは、前記撮像高度の許容範囲（撮像高度の下限値および上限値）に応じて設定される。そして、本実施形態では、画像判定部３５ａは、生育指標処理部３３で求めた生育指標に基づいて、前記判定処理を行う。より具体的には、画像判定部３５ａは、前記判定処理として、生育指標処理部３３で求めたＮＤＶＩ画像において、予め設定された所定の閾値（第１閾値）以下の画素値を持つ画素の総数が予め設定された所定の閾値（第２閾値）以上であるか否かを判定し、前記第２閾値以上ではない場合を画像記憶部４１に記憶する場合と判定し、前記第２閾値以上である場合を画像記憶部４１に記憶しない場合と判定する。前記第１閾値は、圃場のＮＤＶＩと圃場ではない領域のＮＤＶＩとを区別できるように、これらの複数のサンプルから適宜に設定される。前記第２閾値は、全画素数に対する、例えば５０％（画像全体に対し圃場ではない領域が５０％）や、６０％（画像全体に対し圃場ではない領域が６０％）や、７５％（画像全体に対し圃場ではない領域が７５％）等の画素数に設定される。さらに、本実施形態では、画像判定部３５ａは、太陽写込検出部３４の検出結果に基づいて、前記判定処理を行う。より具体的には、画像判定部３５ａは、太陽写込検出部３４が太陽の写り込みを検出しない場合を、画像記憶部４１に記憶する場合と判定し、太陽写込検出部３４が太陽の写り込みを検出した場合を、画像記憶部４１に記憶しない場合と判定する。より詳しくは、画像判定部３５ａは、画像周波数処理部３２で求められた周波数に基づく前記判定処理による第１判定結果、生育指標処理部３３で求められたＮＤＶＩ画像に基づく前記判定処理による第２判定結果、および、太陽写込検出部３４の検出結果に基づく前記判定処理による第３判定結果のＡＮＤが演算され、第１ないし第３判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、第１ないし第３判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。

　画像記憶制御部３６は、画像判定部３５ａで画像記憶部４１に記憶すると判定された場合のみ、画像生成部５で生成された画像を画像記憶部４１に記憶するものである。そして、画像記憶制御部３６は、画像判定部３５ａで画像記憶部４１に記憶しないと判定された場合に、画像生成部５で生成された画像であって画像判定部３５ａで画像記憶部４１に記憶しないと判定された前記画像を指し示す所定の情報を非記憶画像情報として画像記憶部４１に記憶する。より具体的には、画像記憶制御部３６は、画像判定部３５ａで画像記憶部４１に記憶すると判定された場合、画像生成部５で画像を生成した回数を表すシリアル番号、画像生成部５で生成され、画像判定部３５ａで記憶すると判定された画像、および、当該画像を撮像する際に位置高度測定部２１で測定された自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を互いに対応付けて画像記憶部４１に記憶し、画像判定部３５ａで画像記憶部４１に記憶しないと判定された場合、前記シリアル番号、および、画像生成部５で生成され、画像判定部３５ａで記憶しないと判定された画像を撮像する際に位置高度測定部２１で測定された自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を互いに対応付けて画像記憶部４１に記憶する。

　次に、本実施形態の動作について説明する。図５は、第１実施形態における植物生育指標測定装置の動作を示すフローチャートである。図６は、第１実施形態における植物生育指標測定装置において、その飛行軌跡と画像記憶位置との関係を説明するための図である。図６Ａは、飛行軌跡を示し、画像生成部５で画像を生成した画像生成位置を符号●で示している。図６Ｂは、図６Ａに示す飛行軌跡において、画像生成部５で生成した画像を記憶した位置を符号●で示している。図７は、図６Ａに示す飛行軌跡において、植物生育指標測定装置によって撮像され生成された画像を説明するための図である。図７Ａは、一例として、図６Ａに符号Ｔ１で示す各画像生成位置で生成された各画像のうちの１つの画像を示し、図７Ｄは、一例として、図６Ａに符号Ｔ２で示す各画像生成位置で生成された各画像のうちの１つの画像を示し、図７Ｂおよび図７Ｃは、一例として、図６Ａに符号Ｔ３で示す各画像生成位置で生成された各画像のうちの２つの画像を示す。

　このような構成の植物生育指標測定装置Ｍａは、その電源が投入されると、必要な各部の初期化を実行し、その稼働を始める。その制御処理プログラムの実行によって、制御処理部３ａには、制御部３１、画像周波数処理部３２、生育指標処理部３３、太陽写込検出部３４、画像判定部３５ａおよび画像記憶制御部３６が機能的に構成される。

　植物生育指標測定装置Ｍａは、誘導飛行または自律飛行による制御処理部３ａの制御に従って、測定対象の圃場の周辺領域から飛び立ち、前記圃場の上空に向かう。そして、植物生育指標測定装置Ｍａは、前記圃場の上空を、略一定の所定の高度で、圃場の一方端から他方端へ飛行し、前記他方端に到達すると、前方に所定距離だけ進むように折り返して前記他方端から前記一方端へ飛行し、前記一方端に到達すると、前方に所定距離だけ進むように折り返して前記一方端から前記他方端へ飛行する。植物生育指標測定装置Ｍａは、以下同様に、飛行することによって、測定対象の圃場の上空を九十九折りに飛行する。このような離陸から測定対象の圃場の上空を飛行している間に、植物生育指標測定装置Ｍａは、所定の時間間隔（サンプリング間隔）で、以下の各処理Ｓ１１～Ｓ１８を繰り返し実行する。

　すなわち、図５において、まず、植物生育指標測定装置Ｍａは、制御処理部３ａの制御部３１によって、互いに同期するように、画像生成部５で画像を生成して取得し、自機状態測定部２ａで位置および高度を取得する（Ｓ１１）。画像生成部５では、前記画像として、可視撮像部５１で可視画像が生成され、赤外撮像部５２で近赤外画像が生成される。自機状態測定部２ａでは、位置高度測定部２１で自機の位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）が測定される。これら可視画像、近赤外画像、ならびに、自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）は、記憶部４ａのワーキングメモリに一旦記憶される。

　次に、植物生育指標測定装置Ｍａは、制御処理部３ａの画像周波数処理部３２によって、処理Ｓ１１で取得された画像に基づいて、前記条または条間に応じた周波数を求め、この求めた周波数を記憶部４ａのワーキングメモリに一旦記憶する（Ｓ１２）。

　次に、植物生育指標測定装置Ｍａは、制御処理部３ａの生育指標処理部３３によって、処理Ｓ１１で取得された画像に基づいて、生育指標、本実施形態ではＮＤＶＩ画像を求め、この求めたＮＤＶＩ画像を記憶部４ａのワーキングメモリに一旦記憶する（Ｓ１３）。

　次に、植物生育指標測定装置Ｍａは、制御処理部３ａの太陽写込検出部３４によって、処理Ｓ１１で取得された画像に基づいて、太陽の写り込みを検出し、この検出結果を記憶部４ａのワーキングメモリに一旦記憶する（Ｓ１４）。

　次に、植物生育指標測定装置Ｍａは、制御処理部３ａの画像判定部３５ａによって、処理Ｓ１１で取得された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行う（Ｓ１５）。この判定の結果、前記画像を画像記憶部４１に記憶しないと判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍａは、次に、処理Ｓ１６を実行した後に、処理Ｓ１８を実行し、前記画像を画像記憶部４１に記憶すると判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍａは、次に、処理Ｓ１７を実行した後に、処理Ｓ１８を実行する。

　より具体的には、まず、画像判定部３５ａは、前記判定処理として、処理Ｓ１２で画像周波数処理部３２によって求めた周波数が、前記撮像高度対応周波数に、前記第１マージンの範囲内で一致するか否かを判定し、第１判定結果として、前記一致する場合を画像記憶部４１に記憶する場合と判定し、前記一致しない場合を画像記憶部４１に記憶しない場合と判定する。次に、画像判定部３５ａは、前記判定処理として、処理Ｓ１３で生育指標処理部３３によって求めたＮＤＶＩ画像において、前記第１閾値以下の画素値を持つ画素の総数が前記第２閾値以上であるか否かを判定し、第２判定結果として、前記第２閾値以上ではない場合を画像記憶部４１に記憶する場合と判定し、前記第２閾値以上である場合を画像記憶部４１に記憶しない場合と判定する。次に、画像判定部３５ａは、第３判定結果として、太陽写込検出部３４の検出結果に応じて、太陽写込検出部３４が太陽の写り込みを検出しない場合を、画像記憶部４１に記憶する場合と判定し、太陽写込検出部３４が太陽の写り込みを検出した場合を、画像記憶部４１に記憶しない場合と判定する。そして、画像判定部３５ａは、これら第１判定結果、第２判定結果および第３判定結果のＡＮＤが演算され、第１ないし第３判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、第１ないし第３判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。この最終的な判定において、画像記憶部４１に記憶しない場合と判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍａは、次に、処理Ｓ１６を実行した後に、処理Ｓ１８を実行し、画像記憶部４１に記憶する場合と判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍａは、次に、処理Ｓ１７を実行した後に、処理Ｓ１８を実行する。

　この処理Ｓ１６では、植物生育指標測定装置Ｍａは、制御処理部３ａの画像記憶制御部３６によって、処理Ｓ１１で取得した画像、すなわち、記憶部４ａのワーキングメモリに記憶されている画像を破棄（削除、消去）し、シリアル番号に、処理Ｓ１１で取得した自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）、すなわち、記憶部４ａのワーキングメモリに記憶されている自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を対応付けて画像記憶部４１に記憶する。

　前記処理Ｓ１７では、植物生育指標測定装置Ｍａは、制御処理部３ａの画像記憶制御部３６によって、シリアル番号に、処理Ｓ１１で取得した画像および自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）、すなわち、記憶部４ａのワーキングメモリに記憶されている画像および自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を対応付けて画像記憶部４１に記憶する。

　そして、前記処理Ｓ１８では、植物生育指標測定装置Ｍａは、制御処理部３ａによって、終了か否かを判定し、終了である場合には（Ｙｅｓ）、本処理を終了し、終了でない場合には（Ｎｏ）、処理を処理Ｓ１１に戻す。

　一例では、例えば、図６Ａに示す飛行軌跡において、植物生育指標測定装置Ｍａが離陸および上昇中である符号Ｔ１では、例えば、図７Ａに示す画像が画像生成部５によって生成され、この生成された画像には、圃場ではない領域が多く写り込んでいるため、生育指標処理部３３で求められたＮＤＶＩ画像に基づく前記判定処理による第２判定結果が、画像記憶部４１に記憶しない場合となる。したがって、符号Ｔ１では、処理Ｓ１１ないし処理Ｓ１５の各処理が実行された結果、処理Ｓ１６が実行され、前記画像が破棄され、画像記憶部４１には、非記憶画像情報が記憶される。なお、図示しないが、植物生育指標測定装置Ｍａが降下中および着陸である場合も同様である。

　また例えば、図６Ａに示す飛行軌跡において、植物生育指標測定装置Ｍａが圃場の上空を適度な高度で飛行中である符号Ｔ３では、例えば、図７Ｂや図７Ｃに示す画像が画像生成部５によって生成され、この生成された画像には、圃場が種に写り込んでいるため、前記第１ないし第３判定結果の全てが、画像記憶部４１に記憶する場合となる。したがって、符号Ｔ３では、処理Ｓ１１ないし処理Ｓ１５の各処理が実行された結果、処理Ｓ１７が実行され、前記画像が画像記憶部４１に記憶される。

　また例えば、図６Ａに示す飛行軌跡において、植物生育指標測定装置Ｍａが折り返しの飛行中である符号Ｔ２では、例えば、図７Ｄに示す画像が画像生成部５によって生成され、この生成された画像には、圃場ではない領域が多く写り込んでいるため、生育指標処理部３３で求められたＮＤＶＩ画像に基づく前記判定処理による第２判定結果が、画像記憶部４１に記憶しない場合となる。したがって、符号Ｔ２では、処理Ｓ１１ないし処理Ｓ１５の各処理が実行された結果、処理Ｓ１６が実行され、前記画像が破棄され、画像記憶部４１には、非記憶画像情報が記憶される。

　したがって、図６Ａに示す飛行軌跡において、画像生成部５は、図６Ａに示す各符号●の各画像生成位置それぞれで画像を生成するが、各画像生成位置それぞれで図５に示すフローチャートの各処理が実行される結果、画像記憶部４１に記憶される画像は、図６Ｂに示す各符号●の位置で画像生成部５によって生成された画像のみとなる。図６Ａと図６Ｂとを比較すると分かるように、植物生育指標測定装置Ｍａが、離陸中、上昇中および折り返しの飛行中である場合に、画像生成部５によって生成された画像は、破棄され、画像記憶部４１には、非記憶画像情報が記憶され、一方、植物生育指標測定装置Ｍａが、圃場の上空を適度な高度で飛行中である場合に、画像生成部５によって生成された画像は、画像記憶部４１に記憶される。この結果、画像記憶部４１には、生育指標を求めるために有意な画像のみが記憶され、生育指標を求めるために有意ではない画像は、破棄される。

　そして、このように画像記憶部４１に記憶された複数の可視画像および複数の近赤外画像に基づいて複数の生育指標画像が求められ、これら求められた複数の生育指標画像が、測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）を用いた公知の手法によって１枚の生育指標合成画像に合成され、圃場全体の生育指標が得られる。この際に、前記複数の生育指標画像に、不必要な画像が含まれると、無駄な合成の処理が生じる、あるいは、不必要な画像を除去する処理が必要となるが、本実施形態では、前記生育指標を求めるために有意な画像のみを画像記憶部４１に記憶するように、植物生育指標測定装置Ｍａが構成されているので、合成の処理における効率の低下をより抑制できる。なお、このような生育指標画像の演算やその画像合成は、当該植物生育指標測定装置Ｍａで実行されて良く、また、外部装置によって実行されて良い。

　以上説明したように、本実施形態における植物生育指標測定装置Ｍａおよびこれに実装された植物生育指標測定方法は、画像生成部５で生成された画像のうち、画像判定部３５ａで、前記生育指標を求めるために有意な画像として、画像記憶部４１に記憶すると判定された画像のみ、画像記憶部４１に記憶するので、画像記憶部４１の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　上記植物生育指標測定装置Ｍａおよび植物生育指標測定方法は、画像生成部５で生成された画像に基づいて求められた、圃場に生育する作物の条または条間に応じた周波数に基づいて、前記判定処理を行うので、画像生成部５で生成された画像で前記判定処理を実施できる。

　上記植物生育指標測定装置Ｍａおよび植物生育指標測定方法は、画像生成部５で生成された画像に基づいて求められた、前記生育指標に基づいて、前記判定処理を行うので、画像生成部５で生成された画像で前記判定処理を実施できる。

　上記植物生育指標測定装置Ｍａおよび植物生育指標測定方法は、画像生成部５で生成された画像に基づく太陽写込検出部３４の検出結果に基づいて、前記判定処理を行うので、画像生成部５で生成された画像で前記判定処理を実施できる。

　上記植物生育指標測定装置Ｍａおよび植物生育指標測定方法は、画像記憶部４１に記憶しないと判定された画像の非記憶画像情報を記憶するので、画像記憶部４１に記憶しないと判定された画像のみを、後に取り直すことができる。

　次に、別の実施形態について説明する。

（第２実施形態）
　上述の第１実施形態における植物生育指標測定装置Ｍａでは、自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）は、複数の生育指標画像を合成して１枚の生育指標合成画像を生成する際に、これら各複数の生育指標画像における相互の位置関係や実際の撮像高度を表すために、取得され、画像生成部５で生成した画像の記録（記憶）の許否は、前記自機の測定位置や測定高度を必要とせずに、画像生成部５で生成した前記画像を、圃場の特徴の観点から分析することで、決定されたが、第２実施形態における植物生育指標測定装置Ｍｂは、前記自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）に基づいて、画像生成部５で生成した画像の記録（記憶）の許否を決定するものである。

　このような第２実施形態における植物生育指標測定装置Ｍｂは、例えば、図１に示すように、航空機１と、自機状態測定部２ｂと、制御処理部３ｂと、記憶部４ｂと、画像生成部５と、ＩＦ部６とを備える。これら第２実施形態の植物生育指標測定装置Ｍｂにおける航空機１、位置高度測定部２１を備える自機状態測定部２ｂ、画像生成部５およびＩＦ部６は、それぞれ、第１実施形態の植物生育指標測定装置Ｍａにおける航空機１、位置高度測定部２１を備える自機状態測定部２ａ、画像生成部５およびＩＦ部６と同様であるからその説明を省略する。

　記憶部４ｂは、制御処理部３ｂに接続され、制御処理部３ｂの制御に従って、各種の所定のプログラムおよび各種の所定のデータを記憶する回路である。前記各種の所定のプログラムには、例えば、当該植物生育指標測定装置Ｍｂの各部を当該各部の機能に応じて制御する制御プログラムや、位置高度測定部２１で測定された自機の位置および後述の圃場領域情報記憶部４２に記憶された圃場領域情報に基づいて、前記判定処理を行う画像判定プログラムや、前記画像判定プログラムで画像記憶部４１に記憶すると判定された場合のみ、画像生成部５で生成された画像を画像記憶部４１に記憶する画像記憶制御プログラム等の制御処理プログラムが含まれる。前記各種の所定のデータには、画像生成部５で生成された画像（本実施形態では可視画像および近赤外画像）や、測定対象である圃場の所在領域を表す圃場領域情報や、前記判定処理で用いる予め設定された閾値等の、各プログラムを実行する上で必要なデータ等が含まれる。記憶部４ｂは、記憶部４ａと同様に、例えばＲＯＭ、ＥＥＰＲＯＭおよびＲＡＭ等を備える。

　そして、第２実施形態では、記憶部４ｂは、画像記憶部４１に加えてさらに、図１に破線で示すように、前記圃場領域情報を記憶する圃場領域情報記憶部４２を機能的に備える。記憶部４ｂの画像記憶部４１は、記憶部４ａの画像記憶部４１と同様であるので、その説明を省略する。圃場領域情報記憶部４２は、測定対象である圃場の所在領域を表す圃場領域情報を記憶するものである。例えば、本実施形態では、位置高度測定部２１が測定位置として緯度Ｘおよび経度Ｙを出力するので、これに応じて、圃場領域情報記憶部４２は、圃場の輪郭形状を表す複数の点の各位置の緯度Ｘおよび経度Ｙを前記圃場領域情報として記憶する。圃場が例えば矩形である場合には、前記矩形における４個の各頂点の各位置の緯度Ｘおよび経度Ｙが前記圃場領域情報として圃場領域情報記憶部４２に記憶される。

　制御処理部３ｂは、植物生育指標測定装置Ｍｂの各部を当該各部の機能に応じてそれぞれ制御し、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成するための回路である。制御処理部３ｂは、制御処理部３ａと同様に、例えば、ＣＰＵおよびその周辺回路を備えて構成される。制御処理部３ｂは、第２実施形態では、機能的に、制御部３１と、画像判定部３５ｂと、画像記憶制御部３６とを備える。

　制御部３１は、第１実施形態と同様に、当該植物生育指標測定装置Ｍｂの各部を当該各部の機能に応じて制御し、当該植物生育指標測定装置Ｍｂの全体制御を司るものである。

　画像判定部３５ｂは、前記生育指標を求めるために有意な画像のみを画像記憶部４１に記憶するように、画像生成部５で生成された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行うものである。本実施形態では、画像判定部３５ｂは、位置高度測定部２１で測定された自機の測定位置（緯度Ｘ、経度Ｙ）および圃場領域情報記憶部４２に記憶された圃場領域情報に基づいて、前記判定処理を行う。より具体的には、画像判定部３５ｂは、前記判定処理として、位置高度測定部２１で測定された自機の測定位置（緯度Ｘ、経度Ｙ）が、圃場領域情報記憶部４２に記憶された圃場領域情報で表される前記圃場の所在領域内であるか否かを判定し、前記圃場の所在領域内である場合を画像記憶部４１に記憶する場合と判定し、前記圃場の所在領域内ではない場合（前記圃場の所在領域外である場合）を画像記憶部４１に記憶しない場合と判定する。本実施形態では、画像判定部３５ｂは、さらに、位置高度測定部２１で測定された自機の測定高度Ｚに基づいて、前記判定処理を行う。より具体的には、画像判定部３５ｂは、前記判定処理として、位置高度測定部２１で測定された自機の測定高度Ｚが、予め設定された所定の閾値（第４閾値）以上であるか否かを判定し、前記第４閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第４閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。前記第４閾値は、例えば１５ｍや２０ｍや２５ｍ等の、予め設定された所定の撮像高度（飛行高度）に設定される。より詳しくは、画像判定部３５ｂは、自機の測定位置（緯度Ｘ、経度Ｙ）および圃場領域情報に基づく前記判定処理による第１１判定結果、および、自機の測定高度Ｚに基づく前記判定処理による第１２判定結果のＡＮＤが演算され、第１１および第１２判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、第１１および第１２判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。

　画像記憶制御部３６は、上述した第１実施形態の植物生育指標測定装置Ｍａにおける画像記憶制御部３６と同様であるので、その説明を省略する。

　次に、本実施形態の動作について説明する。図８は、第２実施形態における植物生育指標測定装置の動作を示すフローチャートである。

　このような構成の植物生育指標測定装置Ｍｂは、その電源が投入されると、必要な各部の初期化を実行し、その稼働を始める。その制御処理プログラムの実行によって、制御処理部３ｂには、制御部３１、画像判定部３５ｂおよび画像記憶制御部３６が機能的に構成される。そして、第２実施形態における植物生育指標測定装置Ｍｂも、第１実施形態における植物生育指標測定装置Ｍａと同様に、誘導飛行または自律飛行による制御処理部３ｂの制御に従って、測定対象の圃場の周辺領域から飛び立ち、前記圃場の上空に向い、圃場の上空を九十九折りに飛行する。このような離陸から測定対象の圃場の上空を飛行している間に、植物生育指標測定装置Ｍｂは、所定の時間間隔（サンプリング間隔）で、以下の各処理Ｓ２１～Ｓ２６を繰り返し実行する。

　すなわち、図８において、まず、植物生育指標測定装置Ｍｂは、制御処理部３ｂの制御部３１によって、画像生成部５で画像を生成して取得し（Ｓ２１）、これに同期して自機状態測定部２ｂで位置（緯度Ｘ、経度Ｙ）および高度Ｚを取得する（Ｓ２２）。これら可視画像、近赤外画像、ならびに、自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）は、記憶部４ｂのワーキングメモリに一旦記憶される。

　次に、植物生育指標測定装置Ｍｂは、制御処理部３ｂの画像判定部３５ｂによって、処理Ｓ２１で取得された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行う（Ｓ２３）。この判定の結果、前記画像を画像記憶部４１に記憶しないと判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍｂは、次に、処理Ｓ２４を実行した後に、処理Ｓ２６を実行し、前記画像を画像記憶部４１に記憶すると判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍｂは、次に、処理Ｓ２５を実行した後に、処理Ｓ２６を実行する。

　より具体的には、まず、画像判定部３５ｂは、前記判定処理として、位置高度測定部２１で測定された自機の測定位置（緯度Ｘ、経度Ｙ）が、圃場領域情報記憶部４２に記憶された圃場領域情報で表される前記圃場の所在領域内であるか否かを判定し、第１１判定結果として、前記圃場の所在領域内である場合を画像記憶部４１に記憶する場合と判定し、前記圃場の所在領域内ではない場合（前記圃場の所在領域外である場合）を画像記憶部４１に記憶しない場合と判定する。次に、画像判定部３５ｂは、前記判定処理として、位置高度測定部２１で測定された自機の測定高度Ｚが、予め設定された所定の第４閾値以上であるか否かを判定し、第１２判定結果として、前記第４閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第４閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。そして、画像判定部３５ｂは、これら第１１判定結果および第１２判定結果のＡＮＤが演算され、第１１および第１２判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、第１１および第１２判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。この最終的な判定において、画像記憶部４１に記憶しない場合と判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍｂは、次に、処理Ｓ２４を実行した後に、処理Ｓ２６を実行し、画像記憶部４１に記憶する場合と判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍｂは、次に、処理Ｓ２５を実行した後に、処理Ｓ２６を実行する。

　これら第２実施形態の植物生育指標測定装置Ｍｂにおける処理Ｓ２４、処理Ｓ２５および処理Ｓ２６は、それぞれ、第１実施形態の植物生育指標測定装置Ｍａにおける処理Ｓ１６、処理Ｓ１７および処理Ｓ１８と同様であるので、その説明を省略する。

　以上説明したように、第２実施形態における植物生育指標測定装置Ｍｂおよびこれに実装された植物生育指標測定方法は、第１実施形態と同様に、画像生成部５で生成された画像のうち、画像判定部３５ｂで、前記生育指標を求めるために有意な画像として、画像記憶部４１に記憶すると判定された画像のみ、画像記憶部４１に記憶するので、画像記憶部４１の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　また、上記植物生育指標測定装置Ｍｂおよび植物生育指標測定方法は、位置高度測定部２１で測定した自機の測定位置（緯度Ｘ、経度Ｙ）を用いるので、圃場の上空で撮像して生成した画像のみを確実に画像記憶部４１に記憶できる。

　次に、別の実施形態について説明する。

（第３実施形態）
　上述の第１実施形態における植物生育指標測定装置Ｍａでは、画像生成部５で生成した画像の記録（記憶）の許否は、前記自機の測定位置や測定高度を必要とせずに、画像生成部５で生成した前記画像を、圃場の特徴の観点から分析することで、決定されたが、第３実施形態における植物生育指標測定装置Ｍｃは、自機の高度Ｚ、方位φ、速度Ｓおよび加速度Ｋのうちの少なくとも１つに基づいて、画像生成部５で生成した画像の記録（記憶）の許否を決定するものである。

　このような第３実施形態における植物生育指標測定装置Ｍｃは、例えば、図１に示すように、航空機１と、自機状態測定部２ｃと、制御処理部３ｃと、記憶部４ｃと、画像生成部５と、ＩＦ部６とを備える。これら第３実施形態の植物生育指標測定装置Ｍｃにおける航空機１、位置高度測定部２１を備える自機状態測定部２ｃ、画像生成部５およびＩＦ部６は、それぞれ、第１実施形態の植物生育指標測定装置Ｍａにおける航空機１、位置高度測定部２１を備える自機状態測定部２ａ、画像生成部５およびＩＦ部６と同様であるからその説明を省略する。

　記憶部４ｃは、制御処理部３ｃに接続され、制御処理部３ｃの制御に従って、各種の所定のプログラムおよび各種の所定のデータを記憶する回路である。前記各種の所定のプログラムには、例えば、当該植物生育指標測定装置Ｍｃの各部を当該各部の機能に応じて制御する制御プログラムや、自機の高度、方位、速度および加速度のうちの少なくとも１つに基づいて、前記判定処理を行う画像判定プログラムや、前記画像判定プログラムで画像記憶部４１に記憶すると判定された場合のみ、画像生成部５で生成された画像を画像記憶部４１に記憶する画像記憶制御プログラム等の制御処理プログラムが含まれる。前記各種の所定のデータには、画像生成部５で生成された画像（本実施形態では可視画像および近赤外画像）や、前記判定処理で用いる予め設定された各閾値等の、各プログラムを実行する上で必要なデータ等が含まれる。記憶部４ｂは、記憶部４ａと同様に、例えばＲＯＭ、ＥＥＰＲＯＭおよびＲＡＭ等を備える。そして、第３実施形態では、記憶部４ｃは、記憶部４ａと同様に、画像記憶部４１を機能的に備える。

　制御処理部３ｃは、植物生育指標測定装置Ｍｃの各部を当該各部の機能に応じてそれぞれ制御し、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成するための回路である。制御処理部３ｃは、制御処理部３ａと同様に、例えば、ＣＰＵおよびその周辺回路を備えて構成される。制御処理部３ｃは、第３実施形態では、機能的に、制御部３１と、画像判定部３５ｃと、画像記憶制御部３６とを備える。

　制御部３１は、第１実施形態と同様に、当該植物生育指標測定装置Ｍｃの各部を当該各部の機能に応じて制御し、当該植物生育指標測定装置Ｍｂの全体制御を司るものである。

　画像判定部３５ｃは、前記生育指標を求めるために有意な画像のみを画像記憶部４１に記憶するように、画像生成部５で生成された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行うものである。本実施形態では、画像判定部３５ｃは、自機の高度、方位、速度および加速度のうちの少なくとも１つに基づいて、前記判定処理を行う。

　より具体的には、自機の高度Ｚに基づいて前記判定処理を行う場合では、画像判定部３５ｃは、第２実施形態の画像判定部３５ｂと同様に、前記判定処理として、位置高度測定部２１で測定された高度Ｚが、前記第４閾値以上であるか否かを判定し、第２１判定結果として、前記第４閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第４閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。

　自機の方位φに基づいて前記判定処理を行う場合では、制御処理部３ｃは、図１に破線で示すように、位置高度測定部２１で測定された測定位置（緯度Ｘ、経度Ｙ）に基づいて自機の方位φを求める方位処理部３７を機能的にさらに備え、画像判定部３５ｃは、方位処理部３７によって、位置高度測定部２１で測定された今回の測定位置（緯度Ｘ、経度Ｙ）および前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）に基づいて公知の演算手法によって方位φ（＝ｆ（Ｘ、Ｘ_－１、Ｙ、Ｙ_－１））を求め、この求めた方位φが、予め設定された所定の方位（例えば圃場の輪郭における長辺に沿った方向で九十九折りに往復しながら前記画像を生成する場合に前記長辺の方位等）Ｙ_ｒｅｆに、予め設定された所定の第２マージンの範囲内で一致するか否かを判定し、第２２判定結果として、前記一致する場合を画像記憶部４１に記憶する場合と判定し、前記一致しない場合を画像記憶部４１に記憶しない場合と判定する。

　なお、上述では、方位φは、位置高度測定部２１によって、互いに異なる各画像生成位置（互いに異なる各タイミング）で測定された２組の緯度Ｘ、Ｘ_－１および経度Ｙ、Ｙ_－１から方位処理部３７によって求められたが、図１に破線で示すように、自機状態測定部２ｃが、地球上における現在の方位φを測定する方位測定部２２をさらに備え、この方位測定部２２で測定された方位φが上述で用いられても良い。方位測定部２２は、例えば３軸方位計（３軸地磁気センサ）を備えて構成される。

　自機の速度Ｓに基づいて前記判定処理を行う場合では、制御処理部３ｃは、図１に破線で示すように、位置高度測定部２１で測定された測定位置（緯度Ｘ、経度Ｙ）に基づいて自機の速度Ｓを求める速度処理部３８を機能的にさらに備え、画像判定部３５ｃは、前記速度処理部３８によって、位置高度測定部２１で測定された今回の測定位置（緯度Ｘ、経度Ｙ）および前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）ならびにこれの測定時間間隔に基づいて公知の演算手法によって速度Ｓ（＝√（（Ｘ_－１－Ｘ）^２＋（Ｙ_－１－Ｙ）^２）／（前回の測定から今回の測定までの時間ｔ））を求め、前記判定処理として、この求めた速度Ｓが、予め設定された所定の第５閾値以上であるか否かを判定し、第２３判定結果として、前記第５閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第５閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。植物生育指標測定装置Ｍｃは、圃場の上空では、所定の速度（例えば１ｍ／フレームや、２ｍ／フレームや、３ｍ／フレーム等）で飛行するが、必要に応じてホバリング等を行うことがあり、ホバリング中の画像は、略同じであるので、１枚あれば充分である。このため、前記第５閾値は、前記圃場の上空における前記所定の速度と、ホバリングの飛行速度と区別できる値に適宜に設定される。

　なお、上述では、速度Ｓは、位置高度測定部２１によって、互いに異なる各画像生成位置（互いに異なる各タイミング）で測定された２組の緯度Ｘ、Ｘ_－１および経度Ｙ、Ｙ_－１ならびにそれらの測定時間間隔から速度処理部３８によって求められたが、図１に破線で示すように、自機状態測定部２ｃが、自機の速度Ｓを測定する、例えば速度計等の速度測定部２３をさらに備え、この速度測定部２３で測定された速度Ｓが上述で用いられても良い。

　自機の加速度Ｋに基づいて前記判定処理を行う場合では、制御処理部３ｃは、図１に破線で示すように、位置高度測定部２１で測定された測定位置（緯度Ｘ、経度Ｙ）に基づいて自機の加速度Ｋを求める加速度処理部３９を機能的にさらに備え、画像判定部３５ｃは、前記加速度処理部３９によって、位置高度測定部２１で測定された今回の測定位置（緯度Ｘ、経度Ｙ）、前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）および前々回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）ならびにこれらの各測定時間間隔に基づいて公知の演算手法によって加速度Ｋを求め、前記判定処理として、この求めた加速度Ｋが、予め設定された所定の第６閾値以上であるか否かを判定し、第２４判定結果として、前記第６閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第６閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。ホバリングでは、振動が生じることが多く、前記第６閾値は、複数のサンプルから、ホバリングによって生じる振動の加速度を識別できる適宜な値に設定される。

　なお、上述では、加速度Ｋは、位置高度測定部２１によって、互いに異なる各画像生成位置（互いに異なる各タイミング）で測定された３組の緯度Ｘ、Ｘ_－１、Ｘ_－２および経度Ｙ、Ｙ_－１、Ｙ_－２ならびにそれらの各測定時間間隔から加速度処理部３９によって求められたが、図１に破線で示すように、自機状態測定部２ｃが、自機の加速度Ｓを測定する、例えば加速度計等の加速度測定部２４をさらに備え、この加速度測定部２４で測定された加速度Ｓが上述で用いられても良い。

　そして、自機の高度Ｚ、方位φ、速度Ｓおよび加速度Ｋのうちの複数に基づいて、前記判定処理を行う場合では、画像判定部３５ｃは、これら複数それぞれの各判定結果のＡＮＤが演算され、前記複数それぞれの各判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、前記複数それぞれの各判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。

　ここでは、一例として、画像判定部３５ｃは、自機の方位φおよび速度Ｓに基づいて、前記判定処理を行う。

　次に、本実施形態の動作について説明する。図９は、第３実施形態における植物生育指標測定装置の動作を示すフローチャートである。

　このような構成の植物生育指標測定装置Ｍｂは、その電源が投入されると、必要な各部の初期化を実行し、その稼働を始める。その制御処理プログラムの実行によって、制御処理部３ｃには、制御部３１、画像判定部３５ｂ、画像記憶制御部３６、方位処理部３７および速度処理部３８が機能的に構成される。そして、第３実施形態における植物生育指標測定装置Ｍｃも、第１実施形態における植物生育指標測定装置Ｍａと同様に、誘導飛行または自律飛行による制御処理部３ｃの制御に従って、測定対象の圃場の周辺領域から飛び立ち、前記圃場の上空に向い、圃場の上空を九十九折りに飛行する。このような離陸から測定対象の圃場の上空を飛行している間に、植物生育指標測定装置Ｍｃは、所定の時間間隔（サンプリング間隔）で、以下の各処理Ｓ３１～Ｓ３９を繰り返し実行する。

　すなわち、図９において、まず、植物生育指標測定装置Ｍｃは、制御処理部３ｃの制御部３１によって、画像生成部５で画像を生成して取得し（Ｓ３１）、これに同期して自機状態測定部２ｃで位置（緯度Ｘ、経度Ｙ）および高度Ｚを取得する（Ｓ３２）。これら可視画像、近赤外画像、ならびに、自機の測定位置および測定高度（緯度Ｘ、経度Ｙ、高度Ｚ）は、記憶部４ｃのワーキングメモリに一旦記憶される。

　次に、植物生育指標測定装置Ｍｃは、制御処理部３ｃの速度処理部３８によって、処理Ｓ３２で取得された今回の測定位置（緯度Ｘ、経度Ｙ）および前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）ならびにこれの測定時間間隔に基づいて、自機の速度Ｓを求める（Ｓ３３）。この求められた自機の速度Ｓは、記憶部４ｃのワーキングメモリに一旦記憶される。

　次に、植物生育指標測定装置Ｍｃは、制御処理部３ｃの方位処理部３７によって、処理Ｓ３２で取得された今回の測定位置（緯度Ｘ、経度Ｙ）および前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）に基づいて、自機の方位φを求める（Ｓ３４）。この求められた自機の方位φは、記憶部４ｃのワーキングメモリに一旦記憶される。

　次に、植物生育指標測定装置Ｍｃは、制御処理部３ｃによって、今回の測定位置（緯度Ｘ、経度Ｙ）を、前回の測定位置（緯度Ｘ_－１、経度Ｙ_－１）とし、前回の測定位緯度（緯度Ｘ_－１、経度Ｙ_－１）を更新する（Ｘ_－１←Ｘ、Ｙ_－１←Ｙ）（Ｓ３５）。

　次に、植物生育指標測定装置Ｍｃは、制御処理部３ｃの画像判定部３５ｃによって、処理Ｓ３１で取得された画像を画像記憶部４１に記憶するか否かを判定する判定処理を行う（Ｓ３６）。この判定の結果、前記画像を画像記憶部４１に記憶しないと判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍｃは、次に、処理Ｓ３７を実行した後に、処理Ｓ３９を実行し、前記画像を画像記憶部４１に記憶すると判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍｃは、次に、処理Ｓ３８を実行した後に、処理Ｓ３９を実行する。

　より具体的には、まず、画像判定部３５ｃは、処理Ｓ３４で求めた方位φが、前記所定の方位Ｙ_ｒｅｆに、前記所定の第２マージンの範囲内で一致するか否かを判定し、第２２判定結果として、前記一致する場合を画像記憶部４１に記憶する場合と判定し、前記一致しない場合を画像記憶部４１に記憶しない場合と判定する。次に、画像判定部３５ｃは、前記判定処理として、処理Ｓ３３で求めた速度Ｓが、前記第５閾値以上であるか否かを判定し、第２３判定結果として、前記第５閾値以上である場合を画像記憶部４１に記憶する場合と判定し、前記第５閾値以上ではない場合を画像記憶部４１に記憶しない場合と判定する。そして、画像判定部３５ｃは、これら第２２判定結果および第２３判定結果のＡＮＤが演算され、第２２および第２３判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、第２２および第２３判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。この最終的な判定において、画像記憶部４１に記憶しない場合と判定した場合（Ｎｏ）には、植物生育指標測定装置Ｍｃは、次に、処理Ｓ３７を実行した後に、処理Ｓ３９を実行し、画像記憶部４１に記憶する場合と判定した場合（Ｙｅｓ）には、植物生育指標測定装置Ｍｃは、次に、処理Ｓ３８を実行した後に、処理Ｓ３９を実行する。

　これら第３実施形態の植物生育指標測定装置Ｍｃにおける処理Ｓ３７、処理Ｓ３８および処理Ｓ３９は、それぞれ、第１実施形態の植物生育指標測定装置Ｍａにおける処理Ｓ１６、処理Ｓ１７および処理Ｓ１８と同様であるので、その説明を省略する。

　以上説明したように、第３実施形態における植物生育指標測定装置Ｍｃおよびこれに実装された植物生育指標測定方法は、第１実施形態と同様に、画像生成部５で生成された画像のうち、画像判定部３５ｃで、前記生育指標を求めるために有意な画像として、画像記憶部４１に記憶すると判定された画像のみ、画像記憶部４１に記憶するので、画像記憶部４１の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　また、上記植物生育指標測定装置Ｍｃおよび植物生育指標測定方法は、自機の方位φに基づいて前記判定処理を行うので、自機の折り返し中に生成された画像の記憶が低減でき、自機の速度Ｓに基づいて前記判定処理を行うので、自機のホバリング中に生成された画像の記憶が低減できる。

　なお、自機の高度Ｚに基づいて前記判定処理を行う場合では、上記植物生育指標測定装置Ｍｃおよび植物生育指標測定方法は、自機の上昇中または降下中に生成された画像の記憶が低減できる。自機の加速度に基づいて前記判定処理を行う場合では、上記植物生育指標測定装置Ｍｃおよび植物生育指標測定方法は、自機のホバリング中の振動による加速度で前記ホバリングを検出できるから、前記ホバリング中に生成された画像の記憶が低減できる。

　なお、上述の実施形態において、第１ないし第３実施形態における植物生育指標測定装置Ｍａ、Ｍｂ、Ｍｃは、適宜に組み合わされても良い。例えば、第１および第２実施形態における植物生育指標測定装置Ｍａ、Ｍｂが組み合わされる。また例えば、第１および第３実施形態における植物生育指標測定装置Ｍａ、Ｍｃが組み合わされる。また例えば、第２および第３実施形態における植物生育指標測定装置Ｍｂ、Ｍｃが組み合わされる。第１ないし第３実施形態における植物生育指標測定装置Ｍａ、Ｍｂ、Ｍｃ全てが組み合わされる。これら各組合せにおいて、植物生育指標測定装置Ｍは、各判定結果のＡＮＤを演算し、各判定結果の全てが画像記憶部４１に記憶する場合と判定した場合（前記ＡＮＤの演算結果が１である場合）に、最終的に、画像記憶部４１に記憶する場合と判定し、各判定結果のうちの少なくともいずれか１つが画像記憶部４１に記憶しない場合と判定した場合（前記ＡＮＤの演算結果が０である場合）に、最終的に、画像記憶部４１に記憶しない場合と判定する。

　また、上述の実施形態では、画像生成部５は、可視撮像部５１と赤外撮像部５２とを備えて構成されたが、例えば特許第５１６８３５３号公報に開示されているように、可視画像および赤外画像を生成する１個の撮像部を備えて構成されても良い。このような可視画像および赤外画像を生成する１個の撮像部は、測定対象の光学像を所定の第３結像面上に結像する単一の第３結像光学系と、前記第３結像面に受光面を一致させて配置され、前記測定対象の光学像を電気的な信号に変換する第３イメージセンサ単一の光学系と、前記第３イメージセンサ上に配置されたＲ／Ｇ／Ｂ／Ｉｒ、あるいは、Ｗ／Ｙ／Ｒ／Ｉｒとを備えて構成される。前記Ｒは、主に赤色光を透過させるフィルタであり、前記Ｇは、主に緑色光を透過させるフィルタであり、前記Ｂは、主に青色光を透過させるフィルタであり、そして、前記Ｉｒフィルタは、主に赤外光を透過させるフィルタである。前記Ｗは、主に白色光を透過させるフィルタであり、前記Ｙは、主に黄色光を透過させるフィルタである。

　（第４実施形態）
　上述のように、圃場は、通常、前記画像を生成するカメラの画角に対して広いため、圃場全体を１回の撮像による１枚の画像で写すことが難しい。このため、例えば、特開２００３－３３９２３８号公報（特許文献１）では、互いに隣り合う圃場をオーバーラップして複数の画像を取得し、同一位置の圃場で色を補正し、この補正した複数の画像を１枚に合成し、これによって圃場全体の１枚の画像を得て、圃場の画像のＲＧＢ値から窒素相対分布を求める作物の生育診断方法が提案されている（例えば特許文献１の図１１）。

　ところで、前記特許文献１のように、複数の画像を繋ぎ合わせて１枚に合成する、いわゆるパノラマ合成処理を実施する場合、互いに隣接する画像間で連続的に繋ぎ合わせるために、繋ぎ合わせの基準となる繋ぎ合わせ基準点を抽出する必要がある。しかしながら、圃場は、土壌や作物が一面に広がっており、土や葉の反射率が近いため、圃場を写した可視画像は、そのコントラストが低くなる傾向にあることから、前記繋ぎ合わせ基準点の抽出が難しく、合成の精度（繋ぎ合わせの精度）が低下してしまう。また、圃場を写した近赤外画像も、作物が一面に広がっている場合、そのコントラストが低くなる傾向にあり、同様に、合成の精度が低下してしまう。さらに、可視画像や近赤外画像の撮像の際に、太陽が写り込んだり、圃場ではない被写体が写り込んだりすると、露出が大きく変化してしまうため、露出オーバーや露出アンダーが生じてしまうことから、繋ぎ合わせ基準点の抽出にばらつきが生じて合成の精度が低下してしまう。

　そこで、第４実施形態では、より精度良く合成した生育指標の画像を生成できる植物生育指標測定装置、植物生育指標測定方法、植物生育指標測定プログラムおよび植物生育指標測定システムについて開示する。

　この第４実施形態における植物生育指標測定装置は、植物における生育の度合いを表す生育指標を求める装置であり、特に、生育指標を求めるための画像を生成する撮像部（カメラ）によって１回で撮像される１枚の画像に写り込む圃場の領域よりも広い領域での圃場の生育指標画像を生成する。前記生育指標画像は、生育指標を表した画像値の画素で構成された画像である。すなわち、本実施形態では、圃場が、通常、前記撮像部（カメラ）の画角に対して広いため、圃場全体を１回の撮像による１枚の画像で写すことが難しいため、複数の画像が、いわゆるパノラマ合成処理によって、連続的に繋ぎ合わされて１枚の画像（合成画像）に合成される。このために、本実施形態における植物生育指標測定装置は、前記生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得部と、前記画像取得部で所得された複数の生育指標画像を、パノラマ合成処理することによって１枚の生育指標合成画像を生成する画像合成部を備える。

　前記画像取得部は、第１波長帯域での第１光学像を撮像することによって第１画像を生成する第１撮像部と、前記第１波長帯域と異なる第２波長帯域での第２光学像を撮像することによって第２画像を生成する第２撮像部と、前記第１および第２撮像部それぞれによって生成された、圃場の同一領域における第１および第２画像に基づいて前記生育指標画像を生成する生育指標画像生成部とを備え、前記第１および第２撮像部それぞれは、圃場の互いに異なる領域を含むように、複数の第１および第２画像を生成し、前記生育指標画像生成部は、前記複数の第１および第２画像に基づいて前記複数の生育指標画像を生成するように構成されて良いが、ここでは、前記画像取得部は、外部機器から、前記複数の生育指標画像を取得するインターフェース部を備えて構成される場合について以下に説明する。すなわち、前記複数の生育指標画像を生成する画像生成装置と、前記画像生成装置で生成された前記複数の生育指標画像を取得してこれら複数の生育指標画像を合成して１枚の画像（生育指標合成画像）を生成する植物生育指標測定装置とを備える植物生育指標測定システムについて以下に説明する。

　図１０は、第４実施形態における植物生育指標測定システムの構成を示す図である。このような実施形態の植物生育指標測定システムＳは、例えば、図１０に示すように、画像生成装置Ｍと、植物生育指標測定装置Ｐとを備える。

　画像生成装置Ｍは、測定対象の圃場における所定の画像を生成する装置である。前記所定の画像は、生育指標画像を求めるための元画像（例えば生育指標がＮＤＶＩである場合、ＮＤＶＩを求めるための可視画像および近赤外画像）であって良いが、本実施形態では、前記元画像に基づいて生成した生育指標画像である。画像生成装置Ｍは、例えば竿等の長尺なロッドの先端に取り付けられ、上方から圃場を俯瞰した前記圃場の生育指標画像を生成したり、前記圃場に隣接して比較的高い建物があれば、前記建物から前記圃場の生育指標画像を生成したり等しても良いが、本実施形態では、航空機を備え、上空から前記圃場の生育指標画像を生成するように構成されている。

　なお、本実施形態では、生育指標がＮＤＶＩである場合について以下に説明するが、生育指標は、これに限定されるものではなく、他の指標であって良い。例えば、生育指標は、ＲＶＩであって良い（ＲＶＩ＝Ｒｉ’／Ｒｖ、ＲＶＩ（ｘ、ｙ）＝ｒｉ’（ｘ、ｙ）／ｒｖ（ｘ、ｙ））。また例えば、生育指標は、ＤＶＩであって良い（ＤＶＩ＝Ｒｉ’－Ｒｖ、ＤＶＩ（ｘ、ｙ）＝ｒｉ’（ｘ、ｙ）－ｒｖ（ｘ、ｙ））。また例えば、生育指標は、ＴＶＩ（Ｔｒａｎｓｆｏｒｍｅｄ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ）であって良い（ＴＶＩ＝ＮＤＶＩ＋０．５）^０．５）。また例えば、生育指標は、ＩＰＶＩ（Ｉｎｆｒａｒｅｄ　Ｐｅｒｃｅｎｔａｇｅ　Ｖｅｇｅｔａｔｉｏｎ　Ｉｎｄｅｘ）であって良い（ＩＰＶＩ＝Ｒｉ’／（Ｒｉ’＋Ｒｖ）＝（ＮＤＶＩ＋１）／２）。生育指標がこれら他の指標である場合も以下と同様に、互いに異なる圃場の領域を含む複数の生育指標画像を合成した生育指標画像（合成生育指標画像、結合生育指標画像、パノラマ生育指標画像）を生成できる。

　より詳しくは、画像生成装置Ｍは、例えば、図１０に示すように、航空機６１と、位置高度測定部６２と、方位測定部６３と、制御処理部６４と、可視撮像部６５と、赤外撮像部６６と、記憶部６７と、インターフェース部（ＩＦ部）６８とを備える。

　航空機６１は、大気中を飛行する装置であり、例えば、気球、飛行船、飛行機、ヘリコプタおよびマルチコプタ等である。航空機６１は、有人機であって良いが、好ましくは、無線操縦飛行（誘導飛行）または自律飛行による無人機（ドローン）である。航空機６１は、本実施形態では、制御処理部６４に接続され、誘導飛行または自律飛行による制御処理部６４の制御に従って飛行する。

　位置高度測定部６２は、制御処理部６４に接続され、制御処理部６４の制御に従って、地球上における現在の位置および高度を測定する装置であり、例えば、衛星測位システムのＧＰＳ（Ｇｌｏｂａｌ　Ｐｏｓｉｔｉｏｎｉｎｇ　Ｓｙｓｔｅｍ）を備えて構成される。位置高度測定部６２は、その測定結果（測定位置（緯度Ｘ、経度Ｙ、高度Ｚ））を制御処理部６４へ出力する。なお、ＧＰＳは、ＤＧＳＰ（Ｄｉｆｆｅｒｅｎｔｉａｌ　ＧＳＰ）等の誤差を補正する補正機能を持ったＧＰＳであっても良い。

　方位測定部６３は、制御処理部６４に接続され、制御処理部６４の制御に従って、地球上における現在の方位を測定する装置であり、例えば３軸方位計（３軸地磁気センサ）を備えて構成される。方位測定部６３は、その測定結果（測定方位（方位角Ｙ））を制御処理部６４へ出力する。

　可視撮像部６５は、制御処理部６４に接続され、制御処理部６４の制御に従って可視光の画像（可視画像）を生成する装置である。可視撮像部６５は、第１波長帯域（ここでは可視波長帯域）での第１光学像を撮像することによって第１画像（ここでは可視画像）を生成する第１撮像部の一例に相当する。より具体的には、本実施形態では、ＮＤＶＩを求めるために、可視撮像部６５は、例えば、波長６５０ｎｍを中心波長とする比較的狭帯域で光を透過する第１バンドパスフィルタ、前記第１バンドパスフィルタを透過した測定対象の可視光の光学像を所定の結像面上に結像する第１結像光学系、前記第１結像面に受光面を一致させて配置され、前記測定対象の可視光の光学像を電気的な信号に変換する第１イメージセンサ、前記第１イメージセンサの出力に対し公知の画像処理を施して可視画像を生成する第１デジタルシグナルプロセッサ（ＤＳＰ）等を備えて構成される、いわゆるカメラ等である。可視撮像部６５は、この撮像して生成した可視画像を制御処理部６４へ出力する。

　赤外撮像部６６は、制御処理部６４に接続され、制御処理部６４の制御に従って赤外光（例えば近赤外光）の画像（赤外画像）を生成する装置である。赤外撮像部６６は、前記第１波長帯域と異なる第２波長帯域（ここでは赤外波長帯域）での第２光学像を撮像することによって第２画像（ここでは赤外画像）を生成する第２撮像部の一例に相当する。より具体的には、本実施形態では、ＮＤＶＩを求めるために、赤外撮像部６６は、例えば、波長８００ｎｍを中心波長とする比較的狭帯域で光を透過する第２バンドパスフィルタ、前記第２バンドパスフィルタを透過した測定対象の赤外光の光学像を所定の結像面上に結像する第２結像光学系、前記第２結像面に受光面を一致させて配置され、前記測定対象の赤外光の光学像を電気的な信号に変換する第２イメージセンサ、前記第２イメージセンサの出力に対し公知の画像処理を施して赤外画像を生成する第２ＤＳＰ等を備えて構成される、いわゆる赤外カメラ等である。赤外撮像部６６は、この撮像して生成した赤外画像を制御処理部６４へ出力する。

　ＩＦ部６８は、制御処理部６４に接続され、制御処理部６４の制御に従って、例えば植物生育指標測定装置Ｐ等の外部機器との間でデータの入出力を行う回路であり、例えば、シリアル通信方式であるＲＳ２３２Ｃのインターフェース回路、Ｂｌｕｅｔｏｏｔｈ（登録商標）規格を用いたインターフェース回路、ＩｒＤＡ（Ｉｎｆｒａｒｅｄ　Ｄａｔａ　Ａｓｓｃｏｉａｔｉｏｎ）規格等の赤外線通信を行うインターフェース回路、および、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ　Ｓｅｒｉａｌ　Ｂｕｓ）規格を用いたインターフェース回路等である。また、ＩＦ部６８は、有線または無線によって通信する通信カード等であり、例えばイーサネット環境等の通信ネットワークを介して例えば植物生育指標測定装置Ｐや生育指標画像を管理するサーバ装置等の外部装置との間で通信しても良い（イーサネットは登録商標）。

　記憶部６７は、制御処理部６４に接続され、制御処理部６４の制御に従って、各種の所定のプログラムおよび各種の所定のデータを記憶する回路である。前記各種の所定のプログラムには、例えば、当該画像生成装置Ｍの各部を当該各部の機能に応じて制御する制御プログラムや、測定対象の生育指標画像（本実施形態ではＮＤＶＩの画像（ＮＤＶＩ画像））を生成する生育指標画像生成プログラム等の制御処理プログラムが含まれる。前記各種の所定のデータには、可視撮像部６５で生成された可視画像や、赤外撮像部６６で生成された赤外画像や、これら可視画像および赤外画像に基づいて生成された生育指標画像（本実施形態ではＮＤＶＩ画像）等の、各プログラムを実行する上で必要なデータ等が含まれる。記憶部６７は、例えば不揮発性の記憶素子であるＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）や書き換え可能な不揮発性の記憶素子であるＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ　Ｅｒａｓａｂｌｅ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）等を備える。そして、記憶部６７は、前記所定のプログラムの実行中に生じるデータ等を記憶するいわゆる制御処理部６のワーキングメモリとなるＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）等を含む。

　そして、記憶部６７は、これら可視画像、赤外画像および生育指標画像（ＮＤＶＩ画像）を記憶する画像記憶部６７１を機能的に備える。画像記憶部６７１は、可視画像、赤外画像および生育指標画像（ＮＤＶＩ画像）それぞれに、位置高度測定部６２で測定した測定位置および方位測定部６３で測定した測定方位を、対応付けて、これらを記憶する。例えば、測定位置および測定方位は、可視画像、赤外画像および生育指標画像（ＮＤＶＩ画像）それぞれのヘッダに収容されることで、可視画像、赤外画像および生育指標画像（ＮＤＶＩ画像）それぞれに対応付けられる。

　なお、上述では、画像記憶部６７１は、これら可視画像、赤外画像および生育指標画像（ＮＤＶＩ画像）を記憶したが、記憶部６７は、可視画像および赤外画像を一時的に記憶し、これら可視画像および赤外画像に基づいて制御処理部６４によって生育指標画像（ＮＤＶＩ画像）を生成した後に、これら一時的に記憶した可視画像および赤外画像を消去（削除）し、画像記憶部６７１は、測定位置および測定方位を対応付けた生育指標画像（ＮＤＶＩ画像）のみを記憶しても良い。

　制御処理部６４は、画像生成装置Ｍの各部を当該各部の機能に応じてそれぞれ制御し、測定対象の生育指標画像（本実施形態ではＮＤＶＩ画像）を生成するための回路である。制御処理部６４は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）およびその周辺回路を備えて構成される。より詳しくは、制御処理部６４は、航空機６１の飛行を制御し、可視撮像部６５および赤外撮像部６６によって、互いに同期するように撮像して可視画像および赤外画像を生成し、これら生成した可視画像および赤外画像に基づいてＮＤＶＩ画像を求め、位置高度測定部６２および方位測定部６３によって、測定位置および測定方位を取得し、これら取得した測定位置および測定方位を可視画像、赤外画像およびＮＤＶＩ画像に対応付けて、これらを画像記憶部６７１に記憶する。なお、上述したように、画像記憶部６７１は、測定位置および測定方位を対応付けた生育指標画像（ＮＤＶＩ画像）のみを記憶しても良い。

　そして、これら位置高度測定部６２、方位測定部６３、制御処理部１３、可視撮像部６５、赤外撮像部６６、記憶部６７およびＩＦ部６８は、航空機６１に搭載され、適宜な位置に配置される。可視撮像部６５および赤外撮像部６６は、航空機６１によって当該画像生成装置Ｍが飛行している場合に、各光軸が互いに平行で、かつ、下方、好ましくは鉛直下方を向くように、航空機６１に配置される。

　一方、植物生育指標測定装置Ｐは、複数の生育指標画像を、いわゆるパノラマ合成処理によって、連続的に繋ぎ合わせて１枚の生育指標合成画像を生成する装置である。このような植物生育指標測定装置Ｐは、例えば、図１０に示すように、ＩＦ部７１と、制御処理部７２と、記憶部７３と、入力部７４と、出力部７５とを備える。

　入力部７４は、制御処理部７２に接続され、例えば、生育指標合成画像の生成開始を指示するコマンド等の各種コマンド、および、例えば測定対象の圃場の名称の入力等の生育指標を測定する上で必要な各種データを植物生育指標測定装置Ｐに入力する機器であり、例えば、所定の機能を割り付けられた複数の入力スイッチ、キーボードおよびマウス等である。出力部７５は、制御処理部７２に接続され、制御処理部７２の制御に従って、入力部７４から入力されたコマンドやデータ、および、植物生育指標測定装置Ｐによって測定された生育指標合成画像を出力する機器であり、例えばＣＲＴディスプレイ、ＬＣＤ（液晶ディスプレイ）および有機ＥＬディスプレイ等の表示装置やプリンタ等の印刷装置等である。

　なお、入力部７４および出力部７５からタッチパネルが構成されてもよい。このタッチパネルを構成する場合において、入力部７４は、例えば抵抗膜方式や静電容量方式等の操作位置を検出して入力する位置入力装置であり、出力部７５は、表示装置である。このタッチパネルでは、表示装置の表示面上に位置入力装置が設けられ、表示装置に入力可能な１または複数の入力内容の候補が表示され、ユーザが、入力したい入力内容を表示した表示位置を触れると、位置入力装置によってその位置が検出され、検出された位置に表示された表示内容がユーザの操作入力内容として植物生育指標測定装置Ｐに入力される。このようなタッチパネルでは、ユーザは、入力操作を直感的に理解し易いので、ユーザにとって取り扱い易い植物生育指標測定装置Ｐが提供される。

　ＩＦ部７１は、ＩＦ部６８と同様に、制御処理部７２に接続され、制御処理部７２の制御に従って、例えば画像生成装置Ｍ等の外部機器との間でデータの入出力を行う回路であり、例えば、シリアル通信方式であるＲＳ２３２Ｃのインターフェース回路、Ｂｌｕｅｔｏｏｔｈ（登録商標）規格を用いたインターフェース回路、ＩｒＤＡ（規格等の赤外線通信を行うインターフェース回路、および、ＵＳＢ規格を用いたインターフェース回路等である。また、ＩＦ部７１は、有線または無線によって通信する通信カード等であり、例えばイーサネット環境等の通信ネットワークを介して例えば画像生成装置Ｍや生育指標画像を管理するサーバ装置等の外部装置との間で通信しても良い（イーサネットは登録商標）。

　記憶部７３は、制御処理部７２に接続され、制御処理部７２の制御に従って、各種の所定のプログラムおよび各種の所定のデータを記憶する回路である。前記各種の所定のプログラムには、例えば、当該植物生育指標測定装置Ｐの各部を当該各部の機能に応じて制御する制御プログラムや、複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成プログラム等の制御処理プログラムが含まれる。前記画像合成プログラムには、画像の歪曲を補正する歪曲補正処理を行う歪曲補正プログラムや、前記複数の生育指標画像それぞれに対応付けられている各測定位置に基づいて前記複数の生育指標画像それぞれの配置位置を求める位置合せ処理を行う位置合せプログラムや、前記複数の生育指標画像における互いに隣接する２個の第Ａおよび第Ｂ画像を選定し、前記選定した第Ａおよび第Ｂ画像間で複数の対応点を抽出する対応点探索処理を行う対応点探索プログラムや、前記対応点探索プログラムで抽出された複数の対応点に基づいて、前記第Ａ画像に前記第Ｂ画像が連続するように前記第Ｂ画像を回転および並進する画像変換行列を求める変換行列演算処理を行う変換行列演算プログラムや、前記変換行列演算プログラムで求めた画像変換行列を用いて前記第Ｂ画像を画像変換する画像変換処理を行う画像変換プログラムや、前記第Ａ画像と前記画像変換プログラムで画像変換した前記第Ｂ画像とを合成する合成処理を行う合成プログラムが含まれる。前記各種の所定のデータには、前記複数の生育指標画像（本実施形態ではＮＤＶＩ画像）やこれらを１枚に合成した生育指標合成画像（本実施形態ではＮＤＶＩ合成画像）等の、各プログラムを実行する上で必要なデータ等が含まれる。記憶部７３は、例えば不揮発性の記憶素子であるＲＯＭや書き換え可能な不揮発性の記憶素子であるＥＥＰＲＯＭ等を備える。そして、記憶部７３は、前記所定のプログラムの実行中に生じるデータ等を記憶するいわゆる制御処理部６のワーキングメモリとなるＲＡＭ等を含む。なお、記憶部７３は、比較的大容量のデータを記憶できるハードディスク等を備えても良い。

　そして、記憶部７３は、これら複数の生育指標画像（ＮＤＶＩ画像）を記憶する画像記憶部７３１を機能的に備える。より具体的には、測定位置および測定方位を対応付けた複数の生育指標画像（ＮＤＶＩ画像）を記憶する。より詳しくは、画像記憶部７３１は、測定位置および測定方位をヘッダに収容した複数の生育指標画像（ＮＤＶＩ画像）を記憶する。

　なお、上述では、画像記憶部７３１は、測定位置および測定方位を対応付けた複数の生育指標画像（ＮＤＶＩ画像）のみを記憶したが、画像記憶部７３１は、その元画像である、測定位置および測定方位を対応付けた複数の可視画像および複数の赤外画像を記憶しても良い。

　制御処理部７２は、植物生育指標測定装置Ｐの各部を当該各部の機能に応じてそれぞれ制御し、複数の生育指標画像に基づいて１枚の生育指標合成画像を生成するための回路である。制御処理部７２は、例えば、ＣＰＵおよびその周辺回路を備えて構成される。制御処理部７２は、機能的に、制御部７２１と画像合成部７２２とを備える。

　制御部７２１は、当該植物生育指標測定装置Ｐの各部を当該各部の機能に応じて制御し、当該植物生育指標測定装置Ｐの全体制御を司るものである。制御部７２１は、ＩＦ部７１によって、画像生成装置Ｍから、画像生成装置Ｍで圃場の互いに異なる領域を含むように生成された、互いに異なる複数の生育指標画像を取得する。

　画像合成部７２２は、ＩＦ部７１で取得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成するものである。画像合成部７２２は、機能的に、歪曲補正部７２２１と、位置合せ部７２２２と、対応点探索部７２２３と、変換行列演算部７２２４、画像変換部７２２５と、合成部７２２６とを備える。

　歪曲補正部７２２１は、生育指標画像の歪曲を補正する歪曲補正処理を行うものである。

　位置合せ部７２２２は、複数の生育指標画像それぞれに対応付けられている各測定位置および各測定方位に基づいて、複数の生育指標画像それぞれの配置位置を求める位置合せ処理を行うものである。これによって複数の生育指標画像の中で、互いに隣接する育指標画像が求められる。

　対応点探索部７２２３は、前記複数の生育指標画像における互いに隣接する２個の第Ａおよび第Ｂ画像を選定し、前記選定した第Ａおよび第Ｂ画像間で複数の対応点を抽出する対応点探索処理を行うものである。後述の、図１２に示すフローチャートの説明では、前記第Ａ画像は、Ｎ番目の生育指標画像（ＮＤＶＩ画像）Ｉ_Ｎであり、前記第Ｂ画像は、Ｎ＋１番目の生育指標画像（ＮＤＶＩ画像）Ｉ_Ｎ＋１である。

　変換行列演算部７２２４は、対応点探索部７２２３で抽出された複数の対応点に基づいて、前記第Ａ画像に前記第Ｂ画像が連続するように前記第Ｂ画像を回転および並進する画像変換行列を求める変換行列演算処理を行うものである。

　画像変換部７２２５は、変換行列演算部７２２４で求めた画像変換行列を用いて前記第Ｂ画像を画像変換する画像変換処理を行うものである。

　合成部７２２６は、前記第Ａ画像と画像変換部７２２５で画像変換した前記第Ｂ画像とを合成する合成処理を行うものである。

　次に、本実施形態の動作について説明する。図１１は、前記植物生育指標測定システムに用いられる画像生成装置の動作を示すフローチャートである。図１２は、前記植物生育指標測定システムに用いられる植物生育指標測定装置の動作を示すフローチャートである。図１３は、図１２に示すフローチャートにおける対応点探索の処理を説明するための図である。図１３Ａは、互いに隣接するＮＤＶＩ画像のうちの一方を示し、図１３Ｂは、その他方を示す。図１４は、図１２に示すフローチャートにおける画像変換の処理を説明するための図である。図１４は、図１３Ａに示すＮＤＶＩ画像を画像変換したＮＤＶＩ画像である。図１５は、ＮＤＶＩ合成画像の一例を示す図である。図１５は、図１３Ａおよび図１３Ｂに示す各ＮＤＶＩ画像を、一方を画像変換後に合成したＮＤＶＩ合成画像である。

　まず、画像生成装置Ｍの動作について説明する。このような植物生育指標測定システムＳの画像生成装置Ｍは、その電源が投入されると、必要な各部の初期化を実行し、その稼働を始める。

　画像生成装置Ｍは、誘導飛行または自律飛行による制御処理部６４の制御に従って、測定対象の圃場を上空を、略一定の所定の高度で、圃場の一方端から他方端へ飛行し、前記他方端に到達すると、前方に所定距離だけ進むように折り返して前記他方端から前記一方端へ飛行し、前記一方端に到達すると、前方に所定距離だけ進むように折り返して前記一方端から前記他方端へ飛行する。画像生成装置Ｍは、以下同様に、飛行することによって、測定対象の圃場の上空を九十九折りに飛行する。このような測定対象の圃場の上空を飛行している間に、画像生成装置Ｍは、所定の時間間隔（サンプリング間隔）で、以下の各処理を繰り返し実行し、この繰り返しの各回で、可視撮像部６５および赤外撮像部６６の各画角の範囲に相当する前記圃場の一部の領域を可視および近赤外で撮像することで可視画像および近赤外画像を生成し、これらに基づいて生育指標画像を生成している。前記所定の時間間隔は、例えば飛行速度や、飛行高度や、可視撮像部６５および赤外撮像部６６の各画角等を考慮することによって、各可視画像および各近赤外画像が同種の互いに隣接する画像同士において一部で重複するように、設定される。

　より具体的には、図１１において、画像生成装置Ｍは、制御処理部６４の制御に従った可視撮像部６５および赤外撮像部６６によって、圃場を上空から撮像して可視画像および近赤外画像を生成し、取得する（Ｓ１１１）。

　次に、画像生成装置Ｍは、制御処理部６４によって、これら取得した可視画像および近赤外画像に基づいて生育指標画像、本実施形態では、ＮＤＶＩ画像を求める（Ｓ１１２）。ＮＤＶＩ画像の各画素値は、可視画像の画素値をＩとし、近赤外画像の画素値をＲとする場合に、ＮＤＶＩ画像の当該画素位置に対応する可視画像および近赤外画像の画素位置同士において、ＮＤＶＩ＝（Ｉ－Ｒ）／（Ｉ＋Ｒ）によって求められる。例えば、ＮＤＶＩ画像の画素位置（０、０）における画素値は、可視画像の画素位置（０、０）における画素値および近赤外画像の画素位置（０、０）における画素値から求められ、また例えば、ＮＤＶＩ画像の画素位置（１０、１５）における画素値は、可視画像の画素位置（１０、１５）における画素値および近赤外画像の画素位置（１０、１５）における画素値から求められる。なお、可視撮像部６５と、赤外撮像部６６との視差が考慮され、可視画像の画素位置および近赤外画像の画素位置の中の少なくともいずれか一方の画素位置が前記視差を修正するようにシフトされた上で、ＮＤＶＩが求められても良い。また、可視撮像部６５の前記第１結像光学系と赤外撮像部６６の前記第２結像光学系は、例えば画角や歪曲収差等の光学特性が互いに同等であるものとする。

　次に、画像生成装置Ｍは、制御処理部６４の制御に従った位置高度測定部６２および方位測定部６３によって、画像生成装置Ｍにおける位置（緯度Ｘ、経度Ｙ、高度Ｘ）および方位（方位角Ｙ）を測定し、取得する（Ｓ１１３）。

　次に、画像生成装置Ｍは、制御処理部６４によって、処理Ｓ１１２で求めた生育指標画像のヘッダに、処理Ｓ１１３で取得した測定位置（緯度Ｘ、経度Ｙ、高度Ｘ）および測定方位（方位角Ｙ）を収容し（Ｓ１１４）、処理Ｓ１１２で求めた生育指標画像を画像記憶部６７１に記憶して保存する（Ｓ１１５）。

　そして、画像生成装置Ｍは、制御処理部６４によって、終了か否かを判定し、終了である場合には（Ｙｅｓ）、本処理を終了し、終了でない場合には（Ｎｏ）、処理を処理Ｓ１１１に戻す（Ｓ１１６）。

　このように動作することで、画像生成装置Ｍは、前記圃場の一部の各領域それぞれに対応する複数の生育指標画像（ＮＤＶＩ画像）を生成し、これらを画像記憶部６７１に記憶する。

　次に、植物生育指標測定装置Ｐの動作について説明する。上記植物生育指標測定システムＳの植物生育指標測定装置Ｐは、その電源が投入されると、必要な各部の初期化を実行し、その稼働を始める。その制御処理プログラムの実行によって、制御処理部７２には、制御部７２１および画像合成部７２２が機能的に構成され、画像合成部７２２には、歪曲補正部７２２１、位置合せ部７２２２、対応点探索部７２２３、変換行列演算部７２２４、画像変換部７２２５および合成部７２２６が機能的に構成される。そして、画像生成装置Ｍによって生成された圃場の各領域それぞれに対応する複数の生育指標画像（本実施形態ではＮＤＶＩ画像）がＩＦ部７１から取得され、記憶部７３の画像記憶部７３１に記憶される。例えば、画像生成装置ＭのＩＦ部６８と植物生育指標測定装置ＰのＩＦ部７１とが互いに有線または無線によって接続され、ＩＦ部６８およびＩＦ部７１を介して画像生成装置Ｍから植物生育指標測定装置Ｐへ前記複数の生育指標画像（ＮＤＶＩ画像）が転送される。また例えば、ネットワーク上に、前記複数の生育指標画像（ＮＤＶＩ画像）を管理するサーバ装置が設けられ、ＩＦ部６８および前記ネットワークを介して画像生成装置Ｍから前記サーバ装置へ前記複数の生育指標画像（ＮＤＶＩ画像）が転送され、その後、植物生育指標測定装置Ｐの要求に応じて、前記ネットワークおよびＩＦ部７１を介して前記サーバ装置から植物生育指標測定装置Ｐへ前記複数の生育指標画像（ＮＤＶＩ画像）が転送される。また例えば、ＵＳＢメモリ等のメモリカードが用意され、画像生成装置Ｍから、ＩＦ部６８に接続された前記メモリカードへ前記複数の生育指標画像（ＮＤＶＩ画像）が転送され、その後、前記メモリカードをＩＦ部７１に接続して前記メモリカードから植物生育指標測定装置Ｐへ前記複数の生育指標画像（ＮＤＶＩ画像）が転送される。なお、前記メモリカードは、ＩＦ部６８に接続された場合には上述の画像記憶部６７１として使用され、ＩＦ部７１に接続された場合には、上述の画像記憶部７３１として使用されても良い。

　このように前記複数の生育指標画像（ＮＤＶＩ画像）が記憶部７３の画像記憶部７３１に記憶されると、図１２において、まず、植物生育指標測定装置Ｐは、画像合成部７２２の歪曲補正部７２２１によって、前記複数の生育指標画像（ＮＤＶＩ画像）の中の、Ｎ番目の生育指標画像（ＮＤＶＩ画像）Ｉ_Ｎと、Ｎ＋１番目の生育指標画像（ＮＤＶＩ画像）Ｉ_Ｎ＋１を、公知の常套手法によって、例えば可視撮像部６５の前記第１結像光学系や赤外撮像部６６の前記第２結像光学系等に起因する歪曲を補正する歪曲補正処理を行う（Ｓ１２１）。歪曲補正のためのデータは、例えば生育指標画像（ＮＤＶＩ画像）のヘッダに収容されて良く、また例えば入力部７４からユーザ（オペレータ）によって入力されて良い。なお、前記Ｎは、前記複数の生育指標画像（ＮＤＶＩ画像）それぞれに割り付けられたシリアル番号であり、前記複数の生育指標画像（ＮＤＶＩ画像）の総数がＫである場合に、１からＫまでの整数である。前記Ｎは、図１２に示すフローチャートが開始される際に、Ｎ＝１に初期化される。また、図１２に示すフローチャートの以下の説明では、生育指標画像は、その一例としてのＮＤＶＩ画像を指すものとし、特に、括弧書きで注釈しない。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２の位置合せ部７２２２によって、Ｎ番目の生育指標画像Ｉ_Ｎと、Ｎ＋１番目の生育指標画像Ｉ_Ｎ＋１とを、それらヘッダに収容された測定位置（緯度Ｘ、経度Ｙ、高度Ｘ）および測定方位（方位角Ｙ）に基づいて、前記複数の生育指標画像を１枚の画像に合成した生育指標合成画像Ｉｓ上での、Ｎ番目の生育指標画像Ｉ_ＮおよびＮ＋１番目の生育指標画像Ｉ_Ｎ＋１それぞれの配置位置を求める位置合せ処理を行う（Ｓ１２２）。より具体的には、例えば、生育指標合成画像Ｉｓには、基準として基準高度および基準方向が設定され、生育指標画像Ｉは、これら基準高度および基準方向で生成された画像となるように、拡大、縮小および回転等によって修正され、その緯度Ｘおよび経度Ｙが生育指標画像Ｉの中央位置（対角線の交点位置）とみなされ、前記緯度Ｘおよび経度Ｙに応じて生育指標合成画像Ｉｓ上の座標が生育指標画像Ｉに与えられる。これによって生育指標画像Ｉの生育指標合成画像Ｉｓにおける配置位置が決定される。なお、画像生成装置Ｍは、本実施形態では、上述のように、前記複数の可視画像および近赤外画像を生成するので、Ｎ番目の生育指標画像Ｉ_Ｎと、Ｎ＋１番目の生育指標画像Ｉ_Ｎ＋１とは、生育指標合成画像Ｉｓ上において、左右方向（水平方向）および上下方向（垂直方向）のうちのいずれかで隣り合う。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２の対応点探索部７２２３によって、Ｎ番目の生育指標画像Ｉ_Ｎと、Ｎ＋１番目の生育指標画像Ｉ_Ｎ＋１との間で、複数の対応点を抽出する対応点探索処理を行う（Ｓ１２３）。より具体的には、まず、対応点探索部７２２３は、Ｎ番目の生育指標画像Ｉ_ＮおよびＮ＋１番目の生育指標画像Ｉ_Ｎ＋１のうちの一方を、例えばＮ＋１番目の生育指標画像Ｉ_Ｎ＋１を基準画像とし、その他方を、例えばＮ番目の生育指標画像Ｉ_Ｎを参照画像とする。次に、対応点探索部７２２３は、基準画像から複数の特徴点を抽出する。より詳しくは、対応点探索部７２２３は、基準画像において、画素値が周辺画素の画素値より有意に異なる画素を特徴点として抽出する。例えばその画素値が周辺画素の画素値より予め設定された閾値以上で異なる画素が、特徴点として抽出される。互いに隣接する所定の画素数（例えば４画素や９画素等）以上の領域が特徴点として抽出されても良い。そして、対応点探索部７２２３は、基準画像から抽出した各特徴点に対応する各対応点を参照画像から、公知の常套手法によって、探索する。その一例が図１３に示されている。図１３Ａは、基準画像を示し、その●で示す各点は、特徴点Ｐ_Ｎ＋１であり、図１３Ｂは、参照画像を示し、その●で示す各点は、図１３Ａに示す各特徴点Ｐ_Ｎ＋１に対応する対応点Ｐ_Ｎである。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２の変換行列演算部７２２４によって、対応点探索部７２２３で抽出された複数の対応点に基づいて、Ｎ番目の生育指標画像Ｉ_ＮにＮ＋１番目の生育指標画像Ｉ_Ｎ＋１が連続するようにＮ＋１番目の生育指標画像Ｉ_Ｎ＋１を回転および並進する画像変換行列Ｍを求める変換行列演算処理を行う（Ｓ１２４）。より具体的には、画像変換行列Ｍは、例えばアフィン変換を表す行列であり、変換式は、次式１によって表される。ここで、式１における右辺の列ベクトル（ｘ、ｙ、１）は、Ｎ＋１番目の生育指標画像Ｉ_Ｎ＋１から抽出された特徴点の座標（画素位置）であり、式１における左辺の列ベクトル（ｘ’、ｙ’、１）は、Ｎ番目の生育指標画像Ｉ_Ｎから探索された前記特徴点に対応する対応点の座標（画素位置）である。列ベクトル（ｘ、ｙ、１）に作用する画像変換行列Ｍにおける１行１列、１行２列、２行１列および２行２列の各成分ａ、ｂ、ｃ、ｄは、回転の変換行列を表し、１行３列および２行３列の各成分ｔ_ｘ、ｔ_ｙは、並進（平行移動）の変換行列を表す。処理Ｓ１２３で求めた特徴点およびその対応点の各座標（画素位置）が式１に代入され、例えば最小二乗法によって、これら各成分ａ、ｂ、ｃ、ｄ、ｔ_ｘ、ｔ_ｙが求められ、これによって画像変換行列Ｍが求められる。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２の画像変換部７２２５によって、変換行列演算部７２２４で求めた画像変換行列Ｍを用いてＮ＋１番目の生育指標画像Ｉ_Ｎ＋１を画像変換する画像変換処理を行う（Ｓ１２５）。例えば、図１３Ａに示すＮ＋１番目の生育指標画像Ｉ_Ｎ＋１は、画像変換行列Ｍによって図１４に示す画像ＣＩ_Ｎ＋１に変換される。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２の合成部７２２６によって、Ｎ番目の生育指標画像Ｉ_Ｎと画像変換部７２２５で画像変換した変換後のＮ＋１番目の生育指標画像ＣＩ_Ｎ＋１とを合成する合成処理を行う（Ｓ１２６）。より具体的には、合成部７２２６は、生育指標合成画像Ｉｓ上の座標系において、処理Ｓ１２２で求めた配置位置で、Ｎ番目の生育指標画像Ｉ_Ｎと画像変換部７２２５で画像変換した変換後のＮ＋１番目の生育指標画像ＣＩ_Ｎ＋１と画像加算する（Ｉｓ＝Ｉ_Ｎ＋ＣＩ_Ｎ＋１）。例えば、図１３および図１４に示す例では、図１５に示すように、図１３Ｂに示すＮ番目の生育指標画像Ｉ_Ｎと、図１４に示す変換後のＮ＋１番目の生育指標画像ＣＩ_Ｎ＋１とが１枚に重ね合わされる。なお、図１５では、図１３Ｂに示すＮ番目の生育指標画像Ｉ_Ｎと、図１４に示す変換後のＮ＋１番目の生育指標画像ＣＩ_Ｎ＋１とが、視覚的に区別できるように、便宜上、それらの境界を白線で明示しているが、実際の生育指標合成画像Ｉｓでは、前記白線は、不要である。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２によって、次の生育指標画像Ｉを処理するために、Ｎを１だけインクリメントする（Ｎ←Ｎ＋１、Ｓ１２７）。

　次に、植物生育指標測定装置Ｐは、画像合成部７２２によって、全ての生育指標画像Ｉについて処理を終了したか否かを、すなわち、Ｎが最大値ＭＡＸ以上か否かを判定する（Ｓ１２８）。この判定の結果、全ての生育指標画像Ｉについて処理を終了していない場合、すなわち、Ｎが最大値ＭＡＸ以上ではない場合（Ｎｏ）、画像合成部７２２は、処理を処理Ｓ１２１に戻し、一方、前記判定の結果、全ての生育指標画像Ｉについて処理を終了している場合、すなわち、Ｎが最大値ＭＡＸ以上である場合（Ｙｅｓ）、画像合成部７２２は、本処理を終了し、生育指標合成画像Ｉｓを出力部７５に出力する。なお、必要に応じて、ＩＦ部２１から外部機器へ生育指標合成画像Ｉｓが出力されても良い。

　以上説明したように、本実施形態における植物生育指標測定システムＳ、植物生育指標測定装置Ｐならびにこれに実装された植物生育指標測定方法および植物生育指標測定プログラムは、可視画像や近赤外画像では無く、複数の生育指標画像Ｉをパノラマ合成処理するので、より精度良く合成した生育指標合成画像Ｉｓを生成できる。

　図１６は、低コントラストによる合成精度の低下を説明するための図である。図１６Ａは、ＮＤＶＩ画像の一例を示し、図１６Ｂは、可視画像の一例を示し、図１６Ｃは、近赤外画像の一例を示す。図１６Ａに示すＮＤＶＩ画像は、図１６Ｂに示す可視画像と図１６Ｃに示す近赤外画像に基づいて求められた画像である。図１７は、露出のばらつきによる合成精度の低下を説明するための図である。図１７Ａおよび図１７Ｂは、ＮＤＶＩ画像の一例を示し、図１７Ｃは、太陽が水田に写り込んだ場合の画像の一例を示し、図１７Ｄは、露出アンダーの場合の画像の一例を示し、図１７Ｅは、適正露出の場合の画像の一例を示し、図１７Ｆは、図１７Ｄに示す画像と図１７Ｅに示す画像とを合成した場合の画像を示す。

　圃場は、土壌や作物が一面に広がっており、土や葉の反射率が近いため、圃場を写した可視画像は、例えば、図１６Ｂに示すように、そのコントラストが低くなる傾向にある。このため、可視画像では、繋ぎ合わせの基準となる繋ぎ合わせ基準点の抽出が難しく、合成の精度（繋ぎ合わせの精度）が低下してしまう。また、圃場を写した近赤外画像も、作物が一面に広がっている場合、例えば図１６Ｃに示すように、そのコントラストが低くなる傾向にあり、同様に、合成の精度が低下してしまう。しかしながら、圃場における個々の作物が均一に生育することは少ないため、生育指標画像Ｉは、例えば図１６Ａに示すように、可視画像や近赤外画像と比較すると、コントラストが高く、繋ぎ合わせ基準点が抽出し易い。したがって、対応点探索する場合に、基準画像の特徴点が抽出し易く、その参照画像の対応点も探索し易い。

　また、可視画像や近赤外画像の撮像の際に、例えば図１７Ｃに示すように、太陽が写り込んだり、例えば図１７Ｄに示すように、圃場ではない被写体が写り込んだりすると、露出が大きく変化してしまうため、露出オーバーや露出アンダーが生じてしまう。このため、このような画像を合成すると、繋ぎ合わせ基準点の抽出にばらつきが生じて合成の精度が低下してしまう。例えば図１７Ｄに示す露出アンダーの画像と、図１７Ｅに示す、圃場のみが被写体で適正露出の画像とが合成されると、図１７Ｆに示すように、繋ぎ合わせを境界に、その前後（左右）で明暗や色調等に段差が生じ、合成の精度が低下してしまう。しかしながら、生育指標画像Ｉが波長帯域の異なる複数の画像に基づいて生成するので、生育指標画像Ｉは、例えば図１７Ａおよび図１７Ｂに示すように、可視画像や近赤外画像と比較すると、露出差が少ない。

　したがって、上記植物生育指標測定システムＳ、植物生育指標測定装置Ｐ、植物生育指標測定方法および植物生育指標測定プログラムは、パノラマ合成処理で複数の生育指標画像Ｉに基づくことで、より精度良く合成した生育指標合成画像Ｉｓを生成できる。

　なお、上述の第４実施形態では、特徴点と対応点とに基づいて画像変換行列Ｍが求められたが、前記複数の生育指標画像Ｉ_Ｎそれぞれにおける緯度Ｘ_Ｎ、経度Ｙ_Ｎ、高度Ｚ_Ｎおよび方位θ_Ｎに基づいて、アフィン変換を行う画像変換行列Ｍが求められても良い。すなわち、前記生育指標画像Ｉは、第１波長帯域での圃場の第１光学像を前記圃場の上方（例えば上空）から撮像する第１撮像部（例えば生育指標画像がＮＤＶＩ画像である場合では可視撮像部）によって生成された第１画像、および、前記第１波長帯域と異なる第２波長帯域での前記圃場の第２光学像を前記圃場の上方（前記上空）から撮像する第２撮像部（例えば生育指標画像がＮＤＶＩ画像である前記場合では赤外撮像部）によって生成された第２画像に基づいて生成され、前記複数の生育指標画像Ｉ_Ｎそれぞれは、当該生育指標画像Ｉ_Ｎを生成する際に用いられた第１および第２画像それぞれの撮像時における第１および第２撮像部の緯度Ｘ_Ｎ、経度Ｙ_Ｎ、高度Ｚ_Ｎおよび方位θ_Ｎを備え、前記画像合成部は、前記複数の生育指標画像Ｉ_Ｎそれぞれについて、当該生育指標画像Ｉ_Ｎの画素位置を、当該生育指標画像Ｉ_Ｎの備える緯度Ｘ_Ｎ、経度Ｙ_Ｎ、高度Ｚ_Ｎおよび方位θ_Ｎに基づくアフィン変換を行うことによって、前記生育指標合成画像Ｉｓにおける画素位置へ変換することによって前記パノラマ合成処理しても良い。

　図１８は、第１変形形態におけるパノラマ合成処理を説明するための図である。図１８の上段は、各生育指標画像Ｉ_Ｎを示し、図１８の下図は、生育指標合成画像Ｉｓの座標系を示す。

　アフィン変換は、公知のように、線形変換と平行移動（並進）とを組み合わせた変換であり、上述の式１もしくは次式２で表される（例えば、「アフィン変換とは」、［ｏｎｌｉｎｅ］、[２０１６年５月１９日検索]、インターネット（ＵＲＬ：ｈｔｔｐ／／ｄ．ｈａｔｅｎａ．ｎｅ．ｊｐ／Ｚｅｌｌｉｊ／２０１２０５２３／ｐ１））。

　ここで、式２の右辺の列ベクトル（ｘ、ｙ）は、図１８の上段に示すように、生育指標画像Ｉでの画素位置であり（Ｎ番目の生育指標画像Ｉ_Ｎではｘ＝ｘ_Ｎ、ｙ＝ｙ_Ｎ）、列ベクトル（ｘ’、ｙ’）は、図１８の下段に示すように、生育指標合成画像Ｉｓの座標系での画素位置を示す。式２の右辺第１項における、１行１列、１行２列、２行１列および２行２列の各成分ａ、ｂ、ｃ、ｄから成る２行２列の行列は、回転の変換行列Ｒ（θ）を表し、式２の右辺第２項の列ベクトル（ｔ_ｘ、ｔ_ｙ）は、並進（平行移動）の変換行列を表す。

　回転の変換行列Ｒ（θ）における各成分ａ、ｂ、ｃ、ｄは、方位測定部（３軸方位計）１３の値θから、次式３によって与えられる（Ｎ番目の生育指標画像Ｉ_Ｎではθ＝θ_Ｎ）。なお、撮像時の高度が基準高度と異なる場合には、回転行列Ｒ（θ）に、撮像時の高度が基準高度と同等となるように、撮像時の高度に応じたスケーリング係数が乗算される。例えば、成分ａは、ｃｏｓθに前記スケーリング係数を乗算した値である。

　そして、並進（平行移動）の変換行列における各成分ｔ_ｘ、ｔ_ｙは、生育指標合成画像Ｉｓの座標系における原点（０、０）からの距離ｄ［ｍ］と、水平方向（ｘ’軸）との角度φ［ｒａｄ］を、緯度Ｘおよび経度Ｙから、求め、画素数に変換することで、次式４（４－１、４－２）で与えられる。

　ここで、ｐｉｘ１ｍは、メートル単位を画素数の単位に変換するための係数であり、例えば、実測することによって求められる。また、前記距離ｄ［ｍ］および角度φ［ｒａｄ］の算出方法は、公知であり、次式５（５－１、５－２）によって与えられる（例えば「２地点間の距離と方位角」、［ｏｎｌｉｎｅ］、[２０１６年５月２２日検索]、インターネット（ＵＲＬ：ｈｔｔｐ／／ｋｅｉｓａｎ．ｃａｓｉｏ．ｊｐ／ｅｘｅｃ／ｓｙｓｔｅｍ／１２５７６７０７７９））。なお、Ｘ１およびＹ１は、それぞれ、Ａ地点の緯度および経度であり、Ｘ２およびＹ２は、それぞれ、Ｂ地点の緯度および経度であり、ｒは、地球を球体とみなした赤道半径（＝６３７８．１３７ｋｍ）である。ここでは、Ａ地点またはＢ地点が、生育指標合成画像Ｉｓの座標系における原点（０、０）となる。

ｄ＝ｒ×ａｒｃｃｏｓ［ｓｉｎ（Ｙ１）ｓｉｎ（Ｙ２）＋ｃｏｓ（Ｙ１）ｃｏｓ（Ｙ２）ｃｏｓ（Ｘ２－Ｘ１）］　　　・・・（５－１）
φ＝９０－ａｒｃｔａｎ［ｓｉｎ（Ｘ２－Ｘ１），ｃｏｓ（Ｙ１）ｔａｎ（Ｙ２）－ｓｉｎ（Ｙ１）ｃｏｓ（Ｘ２－Ｘ１）］　　　・・・（５－２）

　また、前記画像合成部は、耕作地または前記耕作地に生育する作物に由来する幾何的な所定の目印に基づいて前記パノラマ合成処理しても良い。すなわち、前記画像合成部は、前記生育指標画像Ｉから前記所定の目印を抽出する目印抽出処理を行う目印抽出部と、前記抽出した所定の目印に基づいて、予め設定された所定の基準に対して前記生育指標画像Ｉを回転および並進する画像変換行列を求める変換行列演算処理を行う変換行列演算部と、前記変換行列演算部で求めた画像変換行列を用いて前記生育指標画像Ｉを画像変換する画像変換処理を行う画像変換部と、前記画像変換部で画像変換した前記生育指標画像Ｉを合成する合成処理を行う合成部とを備え、前記複数の生育指標画像Ｉ_Ｎ全てを合成するまで、前記目印抽出部によって前記目印抽出処理を行い、前記変換行列演算部によって前記変換行列演算処理を行い、前記画像変換部によって前記画像変換処理を行い、前記合成部によって前記合成処理を行うことによって前記パノラマ合成処理しても良い。前記所定の基準は、前記複数の生育指標画像Ｉ_Ｎ全てに共通していれば良く、任意に適宜な値に設定される。前記所定の目印は、例えば、前記耕作地における畦である。畦は、水田と水田との間の土地であり、畦には、道や水路が含まれる。また例えば、前記所定の目印は、前記耕作地に生育する稲の条、また、前記稲の条間である。前記所定の目印に対応する生育指標の値は、例えば土壌や雑草等で反射率が作物の反射率と異なるため、作物に対応する生育指標の値と異なり、一般に小さい値（低い値）となる。このため、前記所定の目印は、生育指標画像Ｉの中で比較的目立ち、抽出し易い。したがって、このような植物生育指標測定装置Ｐは、前記所定の目印に基づいてパノラマ合成処理するので、パノラマ合成処理がし易くなる。

　図１９は、第２変形形態におけるパノラマ合成処理を説明するための図である。図１９Ａは、生育指標画像Ｉの一例を示し、図１９Ｂは、図１９Ａにおける中央に位置する領域の画像から、エッジを検出したエッジ検出結果を示し、図１９Ｃは、図１９Ｂに示す画像から、直線のハフ変換によって検出された直線を示す。

　より具体的には、例えば、図１７Ａや図１７Ｂに示すような、前記生育指標画像Ｉから、画素値が予め設定された閾値以下の画素値を持つ画素を抽出することによって、前記所定の目印として畦ＲＰが検出される。この畦ＲＰに基づいて、予め設定された所定の基準に対して前記生育指標画像Ｉを並進する並進の変換行列が求められる。

　また、直線のハフ変換によって、水田に生育する稲の条、また、前記稲の条間が検出され、これらの傾きの角度から、回転の変換行列が求められる。より詳しくは、前記生育指標画像Ｉに対し、エッジを検出するためのいわゆるエッジフィルタが作用され、エッジが検出される。例えば、図１９Ａに示す生育指標画像Ｉに対し、エッジを検出すると、図１９Ｂに示すエッジが検出される。次に、このエッジの各画像から直線のハフ変換によって直線が、稲の条または前記稲の条間として検出される。例えば、図１９Ｂに示すエッジの画像から、図１９Ｃに示すように、複数の直線が検出される。通常、複数の直線が検出されるので、各直線の傾き角度が平均され、平均値の傾き角度と、予め設定された所定の基準角度との差が回転量とされ、回転の変換行列が求められる。

　これら求められた生育指標画像Ｉに対する並進の変換行列および回転の変換行列から、前記生育指標画像Ｉを回転および並進する画像変換行列Ｍが求められる。

　そして、この求めた画像変換行列Ｍで前記生育指標画像Ｉが画像変換され、この画像変換後の前記生育指標画像Ｉが、生育指標合成画像Ｉｓ上の座標系における配置位置に配置され、これに隣接する他の前記生育指標画像Ｉと画像加算され、合成される。

　このような各処理が前記複数の生育指標画像Ｉ_Ｎ全てにつて実施され、生育指標合成画像Ｉｓが生成される。

　また、上述の第４実施形態では、植物生育指標測定装置Ｐは、複数の生育指標画像を画像生成装置Ｍから、ＩＦ部７１を介して取得したが、植物生育指標測定装置Ｐと画像生成装置Ｍとは、一体に構成されても良い。すなわち、上述の第４実施形態では、前記画像取得部は、その一例としてのＩＦ部７１であったが、前記画像取得部は、第１波長帯域での第１光学像を撮像することによって第１画像を生成する第１撮像部（例えば生育指標画像がＮＤＶＩ画像である場合では可視撮像部）と、前記第１波長帯域と異なる第２波長帯域での第２光学像を撮像することによって第２画像を生成する第２撮像部（例えば生育指標画像がＮＤＶＩ画像である前記場合では赤外撮像部）と、前記第１および第２撮像部それぞれによって生成された、圃場の同一領域における第１および第２画像に基づいて前記生育指標画像を生成する生育指標画像生成部とを備え、前記第１および第２撮像部それぞれは、圃場の互いに異なる領域を含むように、複数の第１および第２画像を生成し、前記生育指標画像生成部は、前記複数の第１および第２画像に基づいて前記複数の生育指標画像を生成しても良い。

　あるいは、生育指標画像Ｉは、植物生育指標測定装置Ｐで生成されても良い。すなわち、前記画像取得部は、外部機器から、圃場の互いに異なる領域を含むように生成され、互いに異なる第１および第２波長帯域での複数の第１および第２画像を取得するインターフェース部と、前記インターフェース部で取得した複数の第１および第２画像に基づいて前記複数の生育指標画像を生成する生育指標画像生成部とを備えても良い。

　また、上述の第４実施形態では、画像生成装置Ｍは、可視画像を生成する可視撮像部６５、および、赤外画像を生成する赤外撮像部６６を備えたが、例えば特許第５１６８３５３号公報に開示されているように、可視画像および赤外画像を生成する１個の撮像部を備えても良い。このような可視画像および赤外画像を生成する１個の撮像部は、測定対象の光学像を所定の第３結像面上に結像する単一の第３結像光学系と、前記第３結像面に受光面を一致させて配置され、前記測定対象の光学像を電気的な信号に変換する第３イメージセンサ単一の光学系と、前記第３イメージセンサ上に配置されたＲ／Ｇ／Ｂ／Ｉｒ、あるいは、Ｗ／Ｙ／Ｒ／Ｉｒとを備えて構成される。前記Ｒは、主に赤色光を透過させるフィルタであり、前記Ｇは、主に緑色光を透過させるフィルタであり、前記Ｂは、主に青色光を透過させるフィルタであり、そして、前記Ｉｒは、主に赤外光を透過させるフィルタである。前記Ｗは、主に白色光を透過させるフィルタであり、前記Ｙは、主に黄色光を透過させるフィルタである。

　また、上述の第４実施形態において、画像生成装置Ｍが等速直線運動している場合には、並進量の変動は、少ないので、並進量は、予め設定された所定の固定値とされ、画像から求めなくても良い。

　本明細書は、上記のように様々な態様の技術を開示しているが、そのうち主な技術を以下に纏める。

　一態様にかかる植物生育指標測定装置は、植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成する植物生育指標測定装置であって、所定の時間間隔で画像を生成する画像生成部と、前記画像を記憶する画像記憶部と、前記生育指標を求めるために有意な画像のみを前記画像記憶部に記憶するように、前記画像生成部で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定部と、前記画像判定部で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成部で生成された画像を前記画像記憶部に記憶する画像記憶制御部とを備える。

　このような植物生育指標測定装置は、前記画像生成部で生成された画像のうち、画像判定部で、前記生育指標を求めるために有意な画像として、前記画像記憶部に記憶すると判定された画像のみ、前記画像記憶部に記憶するので、前記画像記憶部の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　他の一態様では、上述の植物生育指標測定装置において、前記画像生成部で生成された画像に基づいて、前記圃場に生育する作物の条または条間に応じた周波数を求める画像周波数処理部をさらに備え、前記画像判定部は、前記画像周波数処理部で求めた周波数に基づいて、前記判定処理を行う。好ましくは、上述の植物生育指標測定装置において、前記画像判定部は、前記判定処理として、前記画像周波数処理部で求めた周波数が、所定の撮像高度で撮像した画像に写る前記条または前記条間に応じた周波数（撮像高度対応周波数）に、所定の第１マージンの範囲内で一致するか否かを判定し、前記一致する場合を前記画像記憶部に記憶する場合と判定する。好ましくは、上述の植物生育指標測定装置において、前記画像周波数処理部は、前記画像生成部で生成された画像からエッジを検出し、前記エッジを検出した画像から直線のハフ変換によって複数の直線を前記条または前記条間として求め、前記求めた複数の直線の間隔から前記周波数を求める。

　このような植物生育指標測定装置は、前記画像生成部で生成された画像に基づいて求められた、前記圃場に生育する作物の条または条間に応じた周波数に基づいて、前記判定処理を行うので、前記画像生成部で生成された画像で前記判定処理を実施できる。

　他の一態様では、上述の植物生育指標測定装置において、前記画像生成部で生成された画像に基づいて、前記生育指標を求める生育指標処理部をさらに備え、前記画像判定部は、前記生育指標処理部で求めた生育指標に基づいて、前記判定処理を行う。好ましくは、上述の植物生育指標測定装置において、前記画像生成部は、可視画像を生成する可視撮像部と、近赤外画像を生成する赤外撮像部とを備え、前記生育指標処理部は、前記画像生成部の可視撮像部および赤外撮像部それぞれで生成された可視画像および近赤外画像に基づいて、前記生育指標の一例に相当するＮＤＶＩを表した画素値を持つＮＤＶＩ画像を求め、前記画像判定部は、前記判定処理として、前記生育指標処理部で求めたＮＤＶＩ画像において、所定の第１閾値以下の画素値を持つ画素の総数が所定の第２閾値以上であるか否かを判定し、前記第２閾値以上ではない場合を前記画像記憶部に記憶する場合と判定する。

　このような植物生育指標測定装置は、前記画像生成部で生成された画像に基づいて求められた、前記生育指標に基づいて、前記判定処理を行うので、前記画像生成部で生成された画像で前記判定処理を実施できる。

　他の一態様では、上述の植物生育指標測定装置において、前記画像生成部で生成された画像に基づいて、太陽の写り込みを検出する太陽写込検出部をさらに備え、前記画像判定部は、前記太陽写込検出部の検出結果に基づいて、前記判定処理を行う。好ましくは、上述の植物生育指標測定装置において、前記太陽写込検出部は、前記画像生成部で生成された画像から、画素値が飽和した画素領域のサイズ（広さ、画素数）を求め、前記求めたサイズが所定の第３閾値以上である画素領域を前記太陽の写り込みとして検出し、前記画像判定部は、前記太陽写込検出部が前記太陽の写り込みを検出しない場合を前記画像記憶部に記憶する場合と判定する。

　このような植物生育指標測定装置は、前記画像生成部で生成された画像に基づく前記太陽写込検出部の検出結果に基づいて、前記判定処理を行うので、前記画像生成部で生成された画像で前記判定処理を実施できる。

　他の一態様では、上述の植物生育指標測定装置において、自機の高度、方位、速度および加速度のうちの少なくとも１つを測定する自機状態測定部と、前記画像生成部、前記画像記憶部、前記画像判定部、前記画像記憶制御部および前記自機状態測定部を搭載し、大気中を飛行する航空機とをさらに備え、前記画像判定部は、前記自機状態測定部の測定結果に基づいて、前記判定処理を行う。好ましくは、上述の植物生育指標測定装置において、前記自機状態測定部は、高度を測定する高度測定部を備え、前記画像判定部は、前記判定処理として、前記自機状態測定部で測定された高度が、所定の第４閾値以上であるか否かを判定し、前記第４閾値以上である場合を前記画像記憶部に記憶する場合と判定する。好ましくは、上述の植物生育指標測定装置において、前記自機状態測定部は、方位を測定する方位測定部を備え、前記画像判定部は、前記自機状態測定部で測定された方位が、所定の方位（例えば圃場の輪郭における長辺に沿った方向で九十九折りに往復しながら前記画像を生成する場合に前記長辺の方位等）に、所定の第２マージンの範囲内で一致するか否かを判定し、前記一致する場合を前記画像記憶部に記憶する場合と判定する。好ましくは、上述の植物生育指標測定装置において、前記自機状態測定部は、速度を測定する速度測定部を備え、前記画像判定部は、前記判定処理として、前記自機状態測定部で測定された速度が、所定の第５閾値以上であるか否かを判定し、前記第５閾値以上である場合を前記画像記憶部に記憶する場合と判定する。好ましくは、上述の植物生育指標測定装置において、前記自機状態測定部は、加速度を測定する加速度測定部を備え、前記画像判定部は、前記判定処理として、前記自機状態測定部で測定された加速度が、所定の第６閾値以上であるか否かを判定し、前記第６閾値以上ではない場合を前記画像記憶部に記憶する場合と判定する。好ましくは、上述の植物生育指標測定装置において、前記自機状態測定部は、自機の高度、方位、速度および加速度のうちの複数を測定し、前記画像判定部は、前記自機状態測定部における複数の測定結果それぞれついて、当該測定結果に基づいて前記判定処理を行い、前記複数の測定結果それぞれについての各判定処理の各判定処理結果のＡＮＤによって、前記画像記憶部に記憶するか否かを判定する。

　自機の高度に基づいて前記判定処理を行う場合、上記植物生育指標測定装置は、自機の上昇中または降下中に生成された画像の記憶が低減できる。自機の方位に基づいて前記判定処理を行う場合、上記植物生育指標測定装置は、自機の折り返し中に生成された画像の記憶が低減できる。自機の速度に基づいて前記判定処理を行う場合、上記植物生育指標測定装置は、自機のホバリング中に生成された画像の記憶が低減できる。自機の加速度に基づいて前記判定処理を行う場合、上記植物生育指標測定装置は、自機のホバリング中の振動による加速度で前記ホバリングを検出できるから、前記ホバリング中に生成された画像の記憶が低減できる。

　他の一態様では、上述の植物生育指標測定装置において、自機の位置を測定する位置測定部と、前記圃場の所在領域を表す圃場領域情報を記憶する圃場領域情報記憶部とをさらに備え、前記画像判定部は、前記位置測定部で測定された自機の位置および前記圃場領域情報記憶部に記憶された圃場領域情報に基づいて、前記判定処理を行う。好ましくは、上述の植物生育指標測定装置において、前記画像判定部は、前記判定処理として、前記位置測定部で測定された自機の位置が、前記圃場領域情報記憶部に記憶された圃場領域情報で表される前記圃場の所在領域内であるか否かを判定し、前記圃場の所在領域内である場合を前記画像記憶部に記憶する場合と判定する。

　このような植物生育指標測定装置は、圃場の上空で撮像して生成した画像のみを確実に画像記憶部に記憶できる。

　他の一態様では、上述の植物生育指標測定装置において、前記画像記憶制御部は、さらに、前記画像判定部で前記画像記憶部に記憶しないと判定された場合に、前記画像生成部で生成された画像であって前記画像判定部で前記画像記憶部に記憶しないと判定された前記画像を指し示す所定の情報を非記憶画像情報として前記画像記憶部に記憶する。好ましくは、上述の植物生育指標測定装置において、自機の位置を測定する位置測定部をさらに備え、前記画像記憶制御部は、前記画像生成部で生成された画像であって前記画像判定部で前記画像記憶部に記憶しないと判定された前記画像を前記画像生成部で生成した際の前記位置測定部で測定された自機の位置を、前記非記憶画像情報として前記画像記憶部に記憶する。

　このような植物生育指標測定装置は、前記画像記憶部に記憶しないと判定された画像の非記憶画像情報を記憶するので、前記画像記憶部に記憶しないと判定された画像のみを、後に取り直すことができる。

　他の一態様にかかる植物生育指標測定方法は、測定対象の圃場を上方から撮像した画像に基づいて、植物における生育の度合いを表す生育指標を求める植物生育指標測定方法であって、所定の時間間隔で画像を生成する画像生成工程と、前記生育指標を求めるために有意な画像のみを画像記憶部に記憶するように、前記画像生成工程で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定工程と、前記画像判定工程で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成工程で生成された画像を前記画像記憶部に記憶する画像記憶制御工程とを備える。

　このような植物生育指標測定方法は、前記画像生成工程で生成された画像のうち、画像判定工程で、前記生育指標を求めるために有意な画像として、前記画像記憶部に記憶すると判定された画像のみ、前記画像記憶部に記憶するので、前記画像記憶部の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　他の一態様にかかる植物生育指標測定プログラムは、測定対象の圃場を上方から所定の時間間隔で撮像した複数の画像に基づいて、植物における生育の度合いを表す生育指標を求める植物生育指標測定プログラムであって、コンピュータに、前記所定の時間間隔で撮像した複数の画像に対し、前記生育指標を求めるために有意な画像のみを画像記憶部に記憶するように、前記画像生成工程で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定工程と、前記画像判定工程で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成工程で生成された画像を前記画像記憶部に記憶する画像記憶制御工程とを実行させるための植物生育指標測定プログラムである。

　このような植物生育指標測定プログラムは、前記所定の時間間隔で撮像した複数の画像のうち、画像判定工程で、前記生育指標を求めるために有意な画像として、前記画像記憶部に記憶すると判定された画像のみ、前記画像記憶部に記憶するので、前記画像記憶部の記憶容量をより低減でき、合成の処理における効率の低下をより抑制できる。

　一態様にかかる植物生育指標測定装置は、植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得部と、前記画像取得部で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成部とを備える。好ましくは、上述の植物生育指標測定装置において、前記画像取得部は、第１波長帯域での第１光学像を撮像することによって第１画像を生成する第１撮像部と、前記第１波長帯域と異なる第２波長帯域での第２光学像を撮像することによって第２画像を生成する第２撮像部と、前記第１および第２撮像部それぞれによって生成された、圃場の同一領域における第１および第２画像に基づいて前記生育指標画像を生成する生育指標画像生成部とを備え、前記第１および第２撮像部それぞれは、圃場の互いに異なる領域を含むように、複数の第１および第２画像を生成し、前記生育指標画像生成部は、前記複数の第１および第２画像に基づいて前記複数の生育指標画像を生成する。好ましくは、上述の植物生育指標測定装置において、前記画像取得部は、外部機器から、前記複数の生育指標画像を取得するインターフェース部である。好ましくは、上述の植物生育指標測定装置において、前記画像取得部は、外部機器から、圃場の互いに異なる領域を含むように生成され、互いに異なる第１および第２波長帯域での複数の第１および第２画像を取得するインターフェース部と、前記インターフェース部で取得した複数の第１および第２画像に基づいて前記複数の生育指標画像を生成する生育指標画像生成部とを備える。好ましくは、上述の植物生育指標測定装置において、前記画像合成部は、前記複数の生育指標画像における互いに隣接する２個の第Ａおよび第Ｂ画像を選定し、前記選定した第Ａおよび第Ｂ画像間で複数の対応点を抽出する対応点探索処理を行う対応点探索部と、前記対応点探索部で抽出された複数の対応点に基づいて、前記第Ａ画像に前記第Ｂ画像が連続するように前記第Ｂ画像を回転および並進する画像変換行列を求める変換行列演算処理を行う変換行列演算部と、前記変換行列演算部で求めた画像変換行列を用いて前記第Ｂ画像を画像変換する画像変換処理を行う画像変換部と、前記第Ａ画像と前記画像変換部で画像変換した前記第Ｂ画像とを合成する合成処理を行う合成部とを備え、前記複数の生育指標画像全てを合成するまで、前記対応点探索部によって前記対応点探索処理を行い、前記変換行列演算部によって前記変換行列演算処理を行い、前記画像変換部によって前記画像変換処理を行い、前記合成部によって前記合成処理を行うことによって前記パノラマ合成処理する。好ましくは、上述の植物生育指標測定装置において、前記生育指標画像は、第１波長帯域での圃場の第１光学像を前記圃場の上方（例えば上空）から撮像する第１撮像部によって生成された第１画像、および、前記第１波長帯域と異なる第２波長帯域での前記圃場の第２光学像を前記圃場の上方（前記上空）から撮像する第２撮像部によって生成された第２画像に基づいて生成され、前記複数の生育指標画像それぞれは、当該生育指標画像を生成する際に用いられた第１および第２画像それぞれの撮像時における第１および第２撮像部の緯度、経度、高度および方位を備え、前記画像合成部は、前記複数の生育指標画像それぞれについて、当該生育指標画像の画素位置を、当該生育指標画像の備える緯度、経度、高度および方位に基づくアフィン変換を行うことによって、前記生育指標合成画像における画素位置へ変換することによって前記パノラマ合成処理する。

　本発明者は、圃場における個々の作物が均一に生育することは少ないため、生育指標画像は、可視画像や近赤外画像と比較すると、コントラストが高く、繋ぎ合わせ基準点の抽出がし易いことを見出した。そして、本発明者は、生育指標画像が波長帯域の異なる複数の画像に基づいて生成するので、生育指標画像は、可視画像や近赤外画像と比較すると、露出差が少ないことを見出した。したがって、上記植物生育指標測定装置は、可視画像や近赤外画像では無く、複数の生育指標画像をパノラマ合成処理するので、より精度良く合成した生育指標合成画像を生成できる。

　他の一態様では、上述の植物生育指標測定装置において、前記画像合成部は、耕作地または前記耕作地に生育する作物に由来する幾何的な所定の目印に基づいて前記パノラマ合成処理する。好ましくは、上述の植物生育指標測定装置において、前記所定の目印は、前記耕作地における畦である。畦は、水田と水田との間の土地であり、畦には、道や水路が含まれる。好ましくは、上述の植物生育指標測定装置において、前記所定の目印は、前記耕作地に生育する稲の条、また、前記稲の条間である。好ましくは、上述の植物生育指標測定装置において、前記画像合成部は、前記生育指標画像から前記所定の目印を抽出する目印抽出処理を行う目印抽出部と、前記抽出した所定の目印に基づいて、予め設定された所定の基準に対して前記生育指標画像を回転および並進する画像変換行列を求める変換行列演算処理を行う変換行列演算部と、前記変換行列演算部で求めた画像変換行列を用いて前記生育指標画像を画像変換する画像変換処理を行う画像変換部と、前記画像変換部で画像変換した前記生育指標画像を合成する合成処理を行う合成部とを備え、前記複数の生育指標画像全てを合成するまで、前記目印抽出部によって前記目印抽出処理を行い、前記変換行列演算部によって前記変換行列演算処理を行い、前記画像変換部によって前記画像変換処理を行い、前記合成部によって前記合成処理を行うことによって前記パノラマ合成処理する。

　前記所定の目印に対応する生育指標の値は、作物に対応する生育指標の値と異なり、一般に小さい値（低い値）となるため、前記所定の目印は、生育指標画像の中で比較的目立ち、抽出し易い。このため、上記植物生育指標測定装置は、前記所定の目印に基づいてパノラマ合成処理するので、パノラマ合成処理がし易くなる。

　他の一態様にかかる植物生育指標測定方法は、植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得工程と、前記画像取得工程で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成工程を備える。

　他の一態様にかかる植物生育指標測定プログラムは、コンピュータに、植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得工程と、前記画像取得工程で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成工程とを実行させるためのプログラムである。

　このような植物生育指標測定方法および植物生育指標測定プログラムは、可視画像や近赤外画像では無く、複数の生育指標画像をパノラマ合成処理するので、より精度良く合成した生育指標合成画像を生成できる。

　他の一態様にかかる植物生育指標測定システムは、画像生成装置と、植物生育指標測定装置とを備える植物生育指標測定システムであって、前記画像生成装置は、植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように、互いに異なる前記複数の生育指標画像を生成し、前記植物生育指標測定装置は、これら上述のいずれかの植物生育指標測定装置である。

　このような植物生育指標測定システムは、これら上述のいずれかの植物生育指標測定装置を用いるので、より精度良く合成した生育指標合成画像を生成できる。

　この出願は、２０１６年６月２２日に出願された日本国特許出願特願２０１６－１２３４３２と２０１６年６月２２日に出願された日本国特許出願特願２０１６－１２３４３３とを基礎とするものであり、その内容は、本願に含まれるものである。

　本発明を表現するために、上述において図面を参照しながら実施形態を通して本発明を適切且つ十分に説明したが、当業者であれば上述の実施形態を変更および／または改良することは容易に為し得ることであると認識すべきである。したがって、当業者が実施する変更形態または改良形態が、請求の範囲に記載された請求項の権利範囲を離脱するレベルのものでない限り、当該変更形態または当該改良形態は、当該請求項の権利範囲に包括されると解釈される。

　本発明によれば、植物における生育の度合いを表す生育指標を求める植物生育指標測定装置、植物生育指標測定方法および植物生育指標測定プログラムが提供できる。

Claims

　植物における生育の度合いを表す生育指標を求めるための、測定対象の圃場を上方から撮像した画像、を生成する植物生育指標測定装置であって、
　所定の時間間隔で画像を生成する画像生成部と、
　前記画像を記憶する画像記憶部と、
　前記生育指標を求めるために有意な画像のみを前記画像記憶部に記憶するように、前記画像生成部で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定部と、
　前記画像判定部で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成部で生成された画像を前記画像記憶部に記憶する画像記憶制御部とを備える、
　植物生育指標測定装置。
　前記画像生成部で生成された画像に基づいて、前記圃場に生育する作物の条または条間に応じた周波数を求める画像周波数処理部をさらに備え、
　前記画像判定部は、前記画像周波数処理部で求めた周波数に基づいて、前記判定処理を行う、
　請求項１に記載の植物生育指標測定装置。
　前記画像生成部で生成された画像に基づいて、前記生育指標を求める生育指標処理部をさらに備え、
　前記画像判定部は、前記生育指標処理部で求めた生育指標に基づいて、前記判定処理を行う、
　請求項１に記載の植物生育指標測定装置。
　前記画像生成部で生成された画像に基づいて、太陽の写り込みを検出する太陽写込検出部をさらに備え、
　前記画像判定部は、前記太陽写込検出部の検出結果に基づいて、前記判定処理を行う、
　請求項１に記載の植物生育指標測定装置。
　自機の高度、方位、速度および加速度のうちの少なくとも１つを測定する自機状態測定部と、
　前記画像生成部、前記画像記憶部、前記画像判定部、前記画像記憶制御部および前記自機状態測定部を搭載し、大気中を飛行する航空機とをさらに備え、
　前記画像判定部は、前記自機状態測定部の測定結果に基づいて、前記判定処理を行う、
　請求項１に記載の植物生育指標測定装置。
　自機の位置を測定する位置測定部と、
　前記圃場の所在領域を表す圃場領域情報を記憶する圃場領域情報記憶部とをさらに備え、
　前記画像判定部は、前記位置測定部で測定された自機の位置および前記圃場領域情報記憶部に記憶された圃場領域情報に基づいて、前記判定処理を行う、
　請求項１に記載の植物生育指標測定装置。
　前記画像記憶制御部は、さらに、前記画像判定部で前記画像記憶部に記憶しないと判定された場合に、前記画像生成部で生成された画像であって前記画像判定部で前記画像記憶部に記憶しないと判定された前記画像を指し示す所定の情報を非記憶画像情報として前記画像記憶部に記憶する、
　請求項１ないし請求項６のいずれか１項に記載の植物生育指標測定装置。
　測定対象の圃場を上方から撮像した画像に基づいて、植物における生育の度合いを表す生育指標を求める植物生育指標測定方法であって、
　所定の時間間隔で画像を生成する画像生成工程と、
　前記生育指標を求めるために有意な画像のみを画像記憶部に記憶するように、前記画像生成工程で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定工程と、
　前記画像判定工程で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成工程で生成された画像を前記画像記憶部に記憶する画像記憶制御工程とを備える、
　植物生育指標測定方法。
　測定対象の圃場を上方から所定の時間間隔で撮像した複数の画像に基づいて、植物における生育の度合いを表す生育指標を求める植物生育指標測定プログラムであって、
　コンピュータに、
　前記所定の時間間隔で撮像した複数の画像に対し、前記生育指標を求めるために有意な画像のみを画像記憶部に記憶するように、前記画像生成工程で生成された画像を前記画像記憶部に記憶するか否かを判定する判定処理を行う画像判定工程と、
　前記画像判定工程で前記画像記憶部に記憶すると判定された場合のみ、前記画像生成工程で生成された画像を前記画像記憶部に記憶する画像記憶制御工程と
　を実行させるための植物生育指標測定プログラム。
　植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得部と、
　前記画像取得部で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成部とを備える、
　植物生育指標測定装置。
　前記画像合成部は、耕作地または前記耕作地に生育する作物に由来する幾何的な所定の目印に基づいて前記パノラマ合成処理する、
　請求項１０に記載の植物生育指標測定装置。
　植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得工程と、
　前記画像取得工程で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成工程を備える、
　植物生育指標測定方法。
　コンピュータに、
　植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように生成された、互いに異なる前記複数の生育指標画像を取得する画像取得工程と、
　前記画像取得工程で所得された複数の生育指標画像を、パノラマ合成処理することによって１個の生育指標合成画像を生成する画像合成工程と
　を実行させるための植物生育指標測定プログラム。
　画像生成装置と、植物生育指標測定装置とを備える植物生育指標測定システムであって、
　前記画像生成装置は、植物における生育の度合いを表す生育指標を表した画像値の画素で構成された複数の生育指標画像であって、圃場の互いに異なる領域を含むように、互いに異なる前記複数の生育指標画像を生成し、
　前記植物生育指標測定装置は、請求項１０または請求項１１に記載の植物生育指標測定装置である、
　植物生育指標測定システム。