(57)【要約】【課題】短時間且つ非破壊で、植物体の広範囲の成熟度を得ることができる植物体の検査装置を提供する。【解決手段】植物体の検査装置１００は、検査対象の植物体Ｐ１を撮影し、撮影された画像の予め定められた範囲内の各画素において、植物体が成熟するに従ってスペクトル強度が減少する第１波長帯Ａと、植物体が成熟するに従ってスペクトル強度が増加する第２波長帯Ｂと、における分光反射スペクトルを取得するスペクトル取得部２０と、分光反射スペクトルの第１波長帯及び第２波長帯のスペクトル強度を用いて、植物体の成熟度を算出する演算部３０と、を備える。【選択図】図１
2015184069.tif 000002
【特許請求の範囲】
【請求項１】
検査対象の植物体を撮影し、撮影された画像の予め定められた範囲内の各画素において、前記植物体が成熟するに従ってスペクトル強度が減少する第１波長帯と、前記植物体が成熟するに従ってスペクトル強度が増加する第２波長帯と、における分光反射スペクトルを取得するスペクトル取得部と、
前記分光反射スペクトルの前記第１波長帯及び前記第２波長帯のスペクトル強度を用いて、前記植物体の成熟度を算出する演算部と、
を備えることを特徴とする植物体の検査装置。
【請求項２】
前記演算部は、
複数の前記分光反射スペクトルを平均して平均スペクトルを得て、
前記平均スペクトルの前記第１波長帯の積分値Ａｓｕｍと、前記平均スペクトルの前記第２波長帯の積分値Ｂｓｕｍとを算出し、
前記積分値Ａｓｕｍ，Ｂｓｕｍと前記植物体の前記成熟度との関係を表す予め定められた関係式を用いて、前記成熟度を算出する
ことを特徴とする請求項１に記載の植物体の検査装置。
【請求項３】
前記関係式は、前記成熟度をＸとして、Ｘ＝（Ｂｓｕｍ−Ａｓｕｍ）／（Ｂｓｕｍ＋Ａｓｕｍ）であることを特徴とする請求項２に記載の植物体の検査装置。
【請求項４】
前記第１波長帯のスペクトル強度は、前記植物体に含まれるポルフィリン系色素の量の増加に応じて増加することを特徴とする請求項１から請求項３の何れかに記載の植物体の検査装置。
【請求項５】
前記ポルフィリン系色素は、クロロフィルを含むことを特徴とする請求項４に記載の植物体の検査装置。
【請求項６】
前記第２波長帯のスペクトル強度は、前記植物体に含まれる、カロテノイド系色素と、フラボノイド系色素と、キノイド系色素と、ベタレイン系色素との何れかの量の増加に応じて増加することを特徴とする請求項１から請求項５の何れかに記載の植物体の検査装置。
【請求項７】
前記第１波長帯は、５０５〜５１５ｎｍの波長帯であり、
前記第２波長帯は、６４５〜６６５ｎｍの波長帯であることを特徴とする請求項１から請求項５の何れかに記載の植物体の検査装置。
【請求項８】
前記予め定められた範囲は、前記植物体の画像の全体であることを特徴とする請求項１から請求項７の何れかに記載の植物体の検査装置。
【請求項９】
前記スペクトル取得部は、分光イメージングカメラであることを特徴とする請求項１から請求項８の何れかに記載の植物体の検査装置。
【請求項１０】
前記演算部は、
前記各画素の前記分光反射スペクトルについて、前記第１波長帯の積分値ＬＡｓｕｍと、前記第２波長帯の積分値ＬＢｓｕｍとを算出し、
前記各画素について、前記積分値ＬＡｓｕｍ，ＬＢｓｕｍと前記植物体の局所的な成熟度との関係を表す予め定められた関係式を用いて、前記局所的な成熟度を算出し、
前記各画素について前記局所的な成熟度を階調表示して、前記局所的な成熟度の分布を表す画像を表示する
ことを特徴とする請求項１から請求項９の何れかに記載の植物体の検査装置。
【請求項１１】
検査対象の植物体を撮影し、撮影された画像の予め定められた範囲内の各画素において、前記植物体が成熟するに従ってスペクトル強度が減少する第１波長帯と、前記植物体が成熟するに従ってスペクトル強度が増加する第２波長帯と、における分光反射スペクトルを取得するステップと、
前記分光反射スペクトルの前記第１波長帯及び前記第２波長帯のスペクトル強度を用いて、前記植物体の成熟度を算出するステップと、
を備えることを特徴とする植物体の検査方法。
【請求項１２】
算出された前記成熟度が予め定められたしきい値以上である前記植物体を出荷または収穫するステップを備えることを特徴とする請求項１１に記載の植物体の検査方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、分光反射スペクトルを用いた植物体の検査装置及び検査方法に関する。
【背景技術】
【０００２】
一般的に、果物や野菜などの植物体の成熟度は、見た目の色に基づいて判断している。そして、この成熟度に応じて、植物体の収穫時期や出荷時期が判断されている。つまり、収穫時期や出荷時期の判断は、人間の勘に頼っている部分が多く、新規就農者にとっては適切な判断が容易ではない。そのため、誰もが分かりやすい植物体の成熟度の判断方法が望まれている。
【０００３】
成熟度を栄養成分から判断する場合、植物体の破壊検査が一般的である。しかし、破壊検査では、例えば、検査対象の植物体をフリーズドライして粉末化したものを検査するため、検査対象の植物体は商品にできない。また、検査には数日を要するため、検査結果に応じて他の植物体の収穫時期や出荷時期を判断するには遅すぎる。
【０００４】
また、非破壊検査として、赤外線などの特殊な光を植物体に照射し、その反射光を検出し、反射光の強度から栄養成分を把握する検査が知られている。しかし、この検査は、植物体の一部のスポット検査であるため、植物体の全体の栄養成分を把握することはできない。
【０００５】
なお、他の非破壊検査として、植物体のデジタル画像を撮影し、このデジタル画像から算出した色彩情報と植物体の色素量との相関に関する情報に基づいて、植物体の色彩情報から色素含有量を算出する技術も知られている（特許文献１参照）。しかし、この技術では、精度が十分に高くはない。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００５−３１５８７７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、このような点を考慮してなされたものであり、短時間且つ非破壊で、植物体の広範囲の成熟度を得ることができる植物体の検査装置及び検査方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明による植物体の検査装置は、
検査対象の植物体を撮影し、撮影された画像の予め定められた範囲内の各画素において、前記植物体が成熟するに従ってスペクトル強度が減少する第１波長帯と、前記植物体が成熟するに従ってスペクトル強度が増加する第２波長帯と、における分光反射スペクトルを取得するスペクトル取得部と、
前記分光反射スペクトルの前記第１波長帯及び前記第２波長帯のスペクトル強度を用いて、前記植物体の成熟度を算出する演算部と、
を備えることを特徴とする。
【０００９】
本発明による植物体の検査装置において、
前記演算部は、
複数の前記分光反射スペクトルを平均して平均スペクトルを得て、
前記平均スペクトルの前記第１波長帯の積分値Ａｓｕｍと、前記平均スペクトルの前記第２波長帯の積分値Ｂｓｕｍとを算出し、
前記積分値Ａｓｕｍ，Ｂｓｕｍと前記植物体の前記成熟度との関係を表す予め定められた関係式を用いて、前記成熟度を算出してもよい。
【００１０】
本発明による植物体の検査装置において、
前記関係式は、前記成熟度をＸとして、Ｘ＝（Ｂｓｕｍ−Ａｓｕｍ）／（Ｂｓｕｍ＋Ａｓｕｍ）であってもよい。
【００１１】
本発明による植物体の検査装置において、
前記第１波長帯のスペクトル強度は、前記植物体に含まれるポルフィリン系色素の量の増加に応じて増加してもよい。
【００１２】
本発明による植物体の検査装置において、
前記ポルフィリン系色素は、クロロフィルを含んでもよい。
【００１３】
本発明による植物体の検査装置において、
前記第２波長帯のスペクトル強度は、前記植物体に含まれる、カロテノイド系色素と、フラボノイド系色素と、キノイド系色素と、ベタレイン系色素との何れかの量の増加に応じて増加してもよい。
【００１４】
本発明による植物体の検査装置において、
前記第１波長帯は、５０５〜５１５ｎｍの波長帯であり、
前記第２波長帯は、６４５〜６６５ｎｍの波長帯であってもよい。
【００１５】
本発明による植物体の検査装置において、
前記予め定められた範囲は、前記植物体の画像の全体であってもよい。
【００１６】
本発明による植物体の検査装置において、
前記スペクトル取得部は、分光イメージングカメラであってもよい。
【００１７】
本発明による植物体の検査装置において、
前記演算部は、
前記各画素の前記分光反射スペクトルについて、前記第１波長帯の積分値ＬＡｓｕｍと、前記第２波長帯の積分値ＬＢｓｕｍとを算出し、
前記各画素について、前記積分値ＬＡｓｕｍ，ＬＢｓｕｍと前記植物体の局所的な成熟度との関係を表す予め定められた関係式を用いて、前記局所的な成熟度を算出し、
前記各画素について前記局所的な成熟度を階調表示して、前記局所的な成熟度の分布を表す画像を表示してもよい。
【００１８】
本発明による植物体の検査方法は、
検査対象の植物体を撮影し、撮影された画像の予め定められた範囲内の各画素において、前記植物体が成熟するに従ってスペクトル強度が減少する第１波長帯と、前記植物体が成熟するに従ってスペクトル強度が増加する第２波長帯と、における分光反射スペクトルを取得するステップと、
前記分光反射スペクトルの前記第１波長帯及び前記第２波長帯のスペクトル強度を用いて、前記植物体の成熟度を算出するステップと、
を備えることを特徴とする。
【００１９】
本発明による植物体の検査方法において、
算出された前記成熟度が予め定められたしきい値以上である前記植物体を出荷または収穫するステップを備えてもよい。
【発明の効果】
【００２０】
本発明によれば、短時間且つ非破壊で、植物体の広範囲の成熟度を得ることができる。
【図面の簡単な説明】
【００２１】
【図１】一実施形態に係る植物体の検査装置の概略的な構成を示す図である。
【図２】一実施形態に係る成熟度の算出処理を示すフローチャートである。
【図３】一実施形態に係る局所的な成熟度の分布の画像化処理を示すフローチャートである。
【図４】４つの検体の分光反射スペクトルの平均スペクトルを示す図である。
【図５】（ａ）は、図４の４つの検体の画像であり、（ｂ）は、図４の４つの検体の局所的な成熟度の分布を表す画像であり、（ｃ）は、図４の４つの検体の成熟度を示す表である。
【発明を実施するための形態】
【００２２】
以下、図面を参照して本発明の一実施形態について説明する。なお、本明細書に添付する図面においては、図示と理解のしやすさの便宜上、適宜縮尺および縦横の寸法比等を、実物のそれらから変更し誇張してある。
【００２３】
図１は、一実施形態に係る植物体の検査装置１００の概略的な構成を示す図である。図１に示すように、植物体の検査装置１００は、光源１０と、分光イメージングカメラ（スペクトル取得部）２０と、演算部３０と、を備える。
【００２４】
検査対象の植物体Ｐ１は、成熟に伴い色が変化するものであり、例えば、野菜、果物及び穀物等である。また、植物体Ｐ１は、成熟に伴い増加する色素と、成熟に伴い減少する色素とを含むものである。
【００２５】
光源１０は、検査対象の植物体Ｐ１に可視光を含む照明光Ｌ１を照射する。例えば、照明光Ｌ１の波長は、３５０〜１０５０ｎｍである。照明光Ｌ１の強度の波長依存性は、低い方が好ましい。
【００２６】
分光イメージングカメラ２０は、照明された植物体Ｐ１からの反射光Ｌ２を検出することで、植物体Ｐ１を撮影する。この時、分光イメージングカメラ２０は、撮影された画像の各画素において、分光反射スペクトル（分光イメージングデータ）を取得する。分光イメージングカメラ２０は、例えば、３５０〜１０５０ｎｍの波長の範囲で、５ｎｍ間隔で分光して分光反射スペクトルを取得する。言い換えれば、分光イメージングカメラ２０は、各画素について分光反射スペクトルを合成することにより、カラー画像を取得できる。分光イメージングカメラ２０は、例えば、ハイパースペクトルカメラ  ＮＨシリーズ（エバ・ジャパン株式会社製）を用いてもよい。
【００２７】
なお、分光イメージングカメラ２０は、撮影された画像の全画素において、上記波長の全範囲で分光反射スペクトルを測定できなくてもよい。つまり、分光イメージングカメラ２０は、撮影された画像の予め定められた範囲内の各画素において、植物体Ｐ１が成熟するに従ってスペクトル強度が減少する第１波長帯Ａと、植物体Ｐ１が成熟するに従ってスペクトル強度が増加する第２波長帯Ｂと、における分光反射スペクトルを取得できればよい。
【００２８】
演算部３０は、分光イメージングカメラ２０で取得された複数の分光反射スペクトルの第１波長帯Ａ及び第２波長帯Ｂのスペクトル強度を用いて、植物体Ｐ１の成熟度を算出する。演算部３０は、例えば、コンピュータなどである。具体的な演算部３０の機能は、フローチャートを参照して後述する。
【００２９】
分光反射スペクトルの第１波長帯Ａのスペクトル強度は、植物体Ｐ１に含まれるポルフィリン系色素の量の増加に応じて増加し、その量の減少に応じて減少する。ポルフィリン系色素は、クロロフィルを含む。この色素は、緑色を有し、植物体Ｐ１の成熟に伴い減少する。
【００３０】
分光反射スペクトルの第２波長帯Ｂのスペクトル強度は、植物体Ｐ１に含まれる、カロテノイド系色素と、フラボノイド系色素と、キノイド系色素と、ベタレイン系色素との何れかの量の増加に応じて増加し、その量の減少に応じて減少する。これらの色素は、植物体Ｐ１の成熟に伴い増加する。第２波長帯Ｂは、第１波長帯Ａより長波長側の帯域である。
【００３１】
表１は、カロテノイド系色素、フラボノイド系色素及びベタレイン系色素のそれぞれの色素の一例と、その色素の色と、その色素を多く含む植物体Ｐ１の一例とを示す。
【表１】
2015184069.tif 000003
【００３２】
図２は、一実施形態に係る成熟度Ｘの算出処理を示すフローチャートである。
【００３３】
まず、分光イメージングカメラ２０により、照明光Ｌ１で照明された検体（植物体Ｐ１）を撮影する（ステップＳ１）。この時、前述のように、撮影された画像の各画素において分光反射スペクトルが取得される。
【００３４】
この処理に先立ち、分光イメージングカメラ２０により、検体が存在しない状態で、照明光Ｌ１で照明された白板等を撮影して、補正用の分光反射スペクトルを取得してもよい。このような補正用の分光反射スペクトルを用いて、ステップＳ１で取得された分光反射スペクトルを補正することにより、照明光Ｌ１の強度の波長依存性の影響を除去して、検体の分光反射スペクトルをより正確に取得できる。
【００３５】
次に、演算部３０により、撮影された画像から算出エリアを抽出する（ステップＳ２）。この処理は、例えば、周知の画像処理を用いて検体の画像の全体を算出エリアとして自動的に抽出してもよい。検体の画像の全体を算出エリアとすることで、検体の全体の成熟度Ｘを算出できる。また、このような算出エリアの自動的な抽出は、選別ラインなどにおいて複数の検体の成熟度を算出する場合に好ましい。なお、操作者が検体の画像の任意の領域を指定することで、この領域を算出エリアとして抽出してもよい。
【００３６】
次に、演算部３０により、抽出した算出エリアにおいて、取得された複数の分光反射スペクトルを平均して平均スペクトルを得る（ステップＳ３）。
【００３７】
次に、演算部３０により、平均スペクトルの第１波長帯Ａの積分値Ａｓｕｍと、平均スペクトルの第２波長帯Ｂの積分値Ｂｓｕｍと、を算出する（ステップＳ４）。
【００３８】
最後に、演算部３０により、積分値Ａｓｕｍ，Ｂｓｕｍと植物体Ｐ１の成熟度Ｘとの関係を表す予め定められた関係式を用いて、成熟度Ｘを算出する。本実施形態では、この関係式は、成熟度をＸとして、Ｘ＝（Ｂｓｕｍ−Ａｓｕｍ）／（Ｂｓｕｍ＋Ａｓｕｍ）である。この関係式は、文献：低温科学，vol.67，2009年，503ページに記載されたＮＤＶＩと称される指標の計算式を参考にして導出している。つまり、演算部３０により、この関係式に平均スペクトルの積分値Ａｓｕｍ，Ｂｓｕｍを代入し、成熟度Ｘを算出する（ステップＳ５）。例えば、成熟度Ｘは、図示しない表示部に表示される。成熟度Ｘは、−１から＋１の範囲の値であり、１に近い程、成熟していることを表す。
【００３９】
ステップＳ５の後、算出された成熟度Ｘが予め定められたしきい値以上である植物体Ｐ１を出荷または収穫するステップを備えてもよい。
【００４０】
図３は、一実施形態に係る局所的な成熟度ＬＸの分布の画像化処理を示すフローチャートである。この画像化処理は、図２を参照して説明した成熟度Ｘの算出処理と共に行われる。ここで、「局所的な成熟度」とは、各画素に対応する成熟度のことをいう。
【００４１】
まず、図２のステップＳ１と同じ処理を行う。
次に、演算部３０により、撮影された画像内の各画素の分光反射スペクトルについて、第１波長帯Ａの積分値ＬＡｓｕｍと、第２波長帯Ｂの積分値ＬＢｓｕｍとを算出する（ステップＳ１２）。
【００４２】
次に、各画素について、積分値ＬＡｓｕｍ，ＬＢｓｕｍと植物体Ｐ１の局所的な成熟度ＬＸとの関係を表す予め定められた関係式を用いて、局所的な成熟度ＬＸを算出する。本実施形態では、この関係式は、局所的な成熟度をＬＸとして、ＬＸ＝（ＬＢｓｕｍ−ＬＡｓｕｍ）／（ＬＢｓｕｍ＋ＬＡｓｕｍ）である。つまり、演算部３０により、各画素について、この関係式にスペクトルの積分値ＬＡｓｕｍ，ＬＢｓｕｍを代入し、局所的な成熟度ＬＸを算出する（ステップＳ１３）。
【００４３】
最後に、演算部３０により、各画素について、局所的な成熟度ＬＸを２０段階で階調表示して、局所的な成熟度ＬＸの分布を表す画像を表示する（ステップＳ１４）。画像は、例えば、図示しない表示部に表示される。なお、局所的な成熟度ＬＸは、−１から＋１の範囲の値である。階調表示の階調数は、２０段階に限らない。
【００４４】
次に、上述した植物体の検査装置１００を用いて、植物体Ｐ１であるトマトの成熟度Ｘを算出する一例について、図４，５を参照して説明する。以下の例では、成熟度Ｘが互いに異なる検体１から検体４の４つのトマトについて説明する。
【００４５】
トマトの場合、成熟に伴い、ポルフィリン系色素のクロロフィルが減少して、カロテノイド系色素のリコピンが増加することを検出する。この場合、第１波長帯Ａは、５０５〜５１５ｎｍの波長帯である。第２波長帯Ｂは、６４５〜６６５ｎｍの波長帯である。
【００４６】
図４は、４つの検体の分光反射スペクトルの平均スペクトルを示す図である。図４の横軸は、波長を単位ｎｍで表し、縦軸は、スペクトル強度を任意の単位で表している。また、図４は、波長３５０ｎｍから９５０ｎｍの範囲を示している。ここでは、検体１から検体４のそれぞれについて、トマト全体の画像の分光反射スペクトルを平均している。
【００４７】
図５（ａ）は、図４の４つの検体の画像であり、図５（ｂ）は、図４の４つの検体の局所的な成熟度ＬＸの分布を表す画像であり、図５（ｃ）は、図４の４つの検体の成熟度Ｘを示す表である。
【００４８】
図５（ａ）は、各画素について得られた分光反射スペクトルを合成して、カラー表示した画像である。つまり、この画像は、人間がトマトを肉眼で観察した場合に認識できる像と同等である。
【００４９】
検体１は、図５（ａ）に示すように、肉眼で観察した場合に緑色に観察される。図４に示すように、検体１の分光反射スペクトルは、第１波長帯Ａにおいて４つの検体中で最も強度が強く、第２波長帯Ｂにおいて４つの検体中で最も強度が弱い。よって、検体１の積分値Ａｓｕｍは、４つの検体中で最も大きく、積分値Ｂｓｕｍは最も小さい。従って、検体１の成熟度Ｘは、図５（ｃ）に示すように０．０２７であり、４つの検体中で最低になっている。また、図５（ｂ）に示すように、局所的な成熟度ＬＸは、検体１の全体において低いことが分かる。
【００５０】
検体２は、図５（ａ）に示すように、肉眼で観察した場合に薄い橙色に観察される。図４に示すように、検体２の分光反射スペクトルは、第１波長帯Ａにおいて検体１より強度が弱く、第２波長帯Ｂにおいて検体１より強度が強い。よって、検体２の積分値Ａｓｕｍは、検体１より小さく、積分値Ｂｓｕｍは、検体１より大きい。従って、成熟度Ｘは、図５（ｃ）に示すように０．８１６であり、検体１より高くなっている。また、図５（ｂ）に示すように、局所的な成熟度ＬＸは、検体２の大部分において低いが、一部では高いことが分かる。
【００５１】
検体３は、図５（ａ）に示すように、肉眼で観察した場合に検体２より濃い橙色に観察される。しかし、肉眼では、検体２と検体３の色の違いを明確に判断することは容易ではない。図４に示すように、検体３の分光反射スペクトルは、第１波長帯Ａにおいて検体２より強度が弱く、第２波長帯Ｂにおいて検体２より強度が強い。よって、検体３の積分値Ａｓｕｍは、検体２より小さく、積分値Ｂｓｕｍは、検体２より大きい。従って、成熟度Ｘは、図５（ｃ）に示すように０．９３３であり、検体２より高くなっている。また、図５（ｂ）に示すように、局所的な成熟度ＬＸは、検体３の大部分において高いが、一部では低いことが分かる。
【００５２】
検体４は、図５（ａ）に示すように、肉眼で観察した場合に赤色に観察される。図４に示すように、検体４の分光反射スペクトルは、第１波長帯Ａにおいて検体３より強度が弱く、第２波長帯Ｂにおいて検体３より強度が強い。よって、検体４の積分値Ａｓｕｍは、検体３より小さく、積分値Ｂｓｕｍは、検体３より大きい。従って、成熟度Ｘは、図５（ｃ）に示すように０．９３３であり、検体３より高くなっている。また、図５（ｂ）に示すように、局所的な成熟度ＬＸは、検体４の全体において高いことが分かる。
【００５３】
このように、肉眼では検体２と検体３との色の違いを明確に判断することは容易ではないため、複数のトマトの中から検体３と同程度又はそれ以上の成熟度のものだけを肉眼で正確に選別することは容易ではない。しかし、算出された成熟度Ｘは、検体２と検体３とにおいて明確に異なっている。そのため、成熟度Ｘを用いれば、例えば、選別のしきい値を０．９以上に設定することにより、複数のトマトの中から検体３と同程度又はそれ以上の成熟度のものを正確且つ容易に選別することができる。従って、新規就農者等であっても、植物体Ｐ１の適切な収穫時期や出荷時期を容易に判断できる。
【００５４】
以上で説明したように、本実施形態によれば、分光イメージングカメラ２０で植物体Ｐ１を撮影して取得された分光反射スペクトルの第１波長帯Ａ及び第２波長帯Ｂのスペクトル強度を用いて植物体Ｐ１の成熟度Ｘを算出するようにしているので、短時間且つ非破壊で植物体Ｐ１の成熟度Ｘを得ることができる。
【００５５】
また、撮影された画像の予め定められた範囲内の各画素において得られた分光反射スペクトルを用いているので、植物体Ｐ１の広範囲の成熟度Ｘを算出できる。そして、予め定められた範囲を植物体Ｐ１の画像の全体に設定することで、植物体Ｐ１の全体の成熟度Ｘを算出できる。
【００５６】
さらに、１つの波長のみのスペクトル強度ではなく、ある幅を有する波長帯のスペクトル強度を用いているので、植物体Ｐ１の育った環境に依存せず、高精度に成熟度Ｘを算出できる。即ち、植物体Ｐ１の育った環境（土壌、日射量等）に応じて、植物体Ｐ１が成熟するに従ってスペクトル強度が増加又は減少する波長が若干異なる可能性があるが、このような影響を受けないようにできる。
【００５７】
さらに、成熟に伴って減少する積分値Ａｓｕｍと、成熟に伴って増加する積分値Ｂｓｕｍとを用いるようにしているので、誤差等の影響を抑えて、高精度に成熟度Ｘを算出できる。
【００５８】
また、光源１０は、可視光を含む照明光Ｌ１を照射できればよいため、特殊な光源を必要としない。従って、植物体の検査装置１００を簡単な構成で実現できる。
【００５９】
（変形例）
以上の実施形態では、第１波長帯Ａ及び第２波長帯Ｂは連続した波長帯である一例について説明したが、これに限らない。例えば、第１波長帯Ａ及び第２波長帯Ｂの少なくとも何れかは、離散した複数の波長帯であってもよい。この場合、成熟に伴って複数の波長帯のスペクトル強度が変化するような色素を含む植物体Ｐ１の成熟度Ｘを、正確に算出できる。
【００６０】
また、以上の実施形態では、植物体Ｐ１の表面を観察する一例について説明したが、これに限らない。例えば、表面を観察するより内部を観察した方が成熟度Ｘをより正確に算出できる場合には、植物体Ｐ１の切断面を観察してもよい。この場合、切断した状態で出荷される植物体Ｐ１に適用することが好ましい。また、例えば、皮が除去された植物体Ｐ１の表面を観察してもよい。このような植物体Ｐ１の例として、玉ねぎなどが挙げられる。
【００６１】
また、成熟度Ｘを算出するための関係式として、以上の実施形態で説明した関係式とは異なる式を用いてもよい。
【００６２】
さらに、成熟度Ｘの算出処理は、図２のフローチャートの例に限らない。例えば、図３のフローチャートのステップＳ１３で各画素について局所的な成熟度ＬＸを算出した後に、複数の局所的な成熟度ＬＸを平均して植物体Ｐ１全体の成熟度Ｘを得てもよい。
【００６３】
また、図３のフローチャートを参照して説明した画像化処理は、省略してもよい。
【００６４】
また、以上の説明では、分光イメージングカメラ（スペクトル取得部）２０としてハイパースペクトルカメラを用いる一例について説明したが、これに限らない。分光イメージングカメラ２０として、撮影された画像の予め定められた範囲内の各画素において、第１波長帯Ａと第２波長帯Ｂとにおける分光反射スペクトルを取得できる装置であれば、どのような呼称の装置を用いてもよい。例えば、ハイパースペクトルカメラに代えて、マルチスペクトルカメラを用いても良い。マルチスペクトルカメラとは、スペクトル分解能がハイパースペクトルカメラより大きい装置である。つまり、分光イメージングカメラ２０は、反射画像や分光反射スペクトルが得られるCCD、CMOSなどの受光センサを有する装置であれば良い。
また、光源１０として、光の波長が少なくとも第１波長帯Ａと第２波長帯Ｂとを含むＬＥＤ（Light Emitting Diode）等の光源を用いてもよい。
【００６５】
上述した実施形態で説明した植物体の検査装置１００の演算部３０の少なくとも一部は、ハードウェアで構成してもよいし、ソフトウェアで構成してもよい。ソフトウェアで構成する場合には、演算部３０の少なくとも一部の機能を実現するプログラムをフレキシブルディスクやＣＤ−ＲＯＭ等の記録媒体に収納し、コンピュータに読み込ませて実行させてもよい。記録媒体は、磁気ディスクや光ディスク等の着脱可能なものに限定されず、ハードディスク装置やメモリなどの固定型の記録媒体でもよい。
【００６６】
また、演算部３０の少なくとも一部の機能を実現するプログラムを、インターネット等の通信回線（無線通信も含む）を介して頒布してもよい。さらに、同プログラムを暗号化したり、変調をかけたり、圧縮した状態で、インターネット等の有線回線や無線回線を介して、あるいは記録媒体に収納して頒布してもよい。
【００６７】
以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００６８】
１０  光源
２０  分光イメージングカメラ（スペクトル取得部）
３０  演算部
１００  植物体の検査装置
Ｐ１  植物体
Ａ  第１波長帯
Ｂ  第２波長帯
Ａｓｕｍ  第１波長帯の積分値
Ｂｓｕｍ  第２波長帯の積分値
Ｘ  成熟度
【図１】
2015184069.tif 000004
【図２】
2015184069.tif 000005
【図３】
2015184069.tif 000006
【図４】
2015184069.tif 000007
【図５】
2015184069.tif 000008