(57)【要約】  画像計測システムは、計測対象である植物３１の近傍に設置された反射体１０と、反射体１０を撮影し、撮影画像を出力する撮影手段２１と、撮影画像内の反射体１０の領域を検出し、領域の色情報を取得する画像分析手段４１と、色情報に基づいて、領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段４２と、絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段４４と、光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段４５とを備える。
2014103181.tif 000002
【特許請求の範囲】
【請求項１】
計測対象である植物の近傍に設置された反射体を撮影し、撮影画像を出力し、
前記撮影画像内の前記反射体の領域を検出し、
前記領域の色情報を取得し、
前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出し、
前記絶対分光放射照度分布に基づいて、光量子束密度を算出し、
前記光量子束密度に基づいて、葉面積指数を算出する
ことを特徴とする画像計測方法。
【請求項２】
反射体が、色が既知である完全拡散反射面を有する
請求項１記載の画像計測方法。
【請求項３】
光量子束密度と葉面積指数とを対応付けて記憶し、
算出された光量子束密度に対応する前記葉面積指数を抽出することによって葉面積指数を算出する
請求項１または請求項２記載の画像計測方法。
【請求項４】
色情報に基づいて、領域の輝度値を含む三刺激値を算出し、
前記輝度値に基づいて、前記領域の放射照度の絶対値である絶対放射照度を算出し、
前記三刺激値に基づいて、前記領域の色度を算出し、
前記色度に基づいて、前記領域における分光放射照度の相対分布である相対分光放射照度分布を算出し、
前記絶対放射照度と前記相対分光放射照度分布とに基づいて、絶対分光放射照度分布を算出する
請求項１から請求項３のうちのいずれか１項に記載の画像計測方法。
【請求項５】
計測対象である植物の近傍に設置された反射体と、
前記反射体を撮影し、撮影画像を出力する撮影手段と、
前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析手段と、
前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段と、
前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段と、
前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段とを備えた
ことを特徴とする画像計測システム。
【請求項６】
反射体が、色が既知である完全拡散反射面を有する
請求項５記載の画像計測システム。
【請求項７】
光量子束密度と葉面積指数とを対応付けて記憶する記憶手段を備え、
葉面積指数算出手段は、
算出された光量子束密度に対応する葉面積指数を前記記憶手段から抽出することによって葉面積指数を算出する
請求項５または請求項６記載の画像計測システム。
【請求項８】
絶対分光放射照度分布算出手段は、
色情報に基づいて、領域の輝度値を含む三刺激値を算出し、
前記輝度値に基づいて、前記領域の放射照度の絶対値である絶対放射照度を算出し、
前記三刺激値に基づいて、前記領域の色度を算出し、
前記色度に基づいて、前記領域における分光放射照度の相対分布である相対分光放射照度分布を算出し、
前記絶対放射照度と前記相対分光放射照度分布とに基づいて、絶対分光放射照度分布を算出する
請求項５から請求項７のうちのいずれか１項に記載の画像計測システム。
【請求項９】
計測対象である植物の近傍に設置された反射体と、前記反射体を撮影し、撮影画像を出力する撮影手段とを備えた画像計測システムに用いられる画像計測装置であって、
前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析手段と、
前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段と、
前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段と、
前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段とを備えた
ことを特徴とする画像計測装置。
【請求項１０】
計測対象である植物の近傍に設置された反射体と、前記反射体を撮影し、撮影画像を出力する撮影手段とを備えた画像計測システムに用いられる画像計測プログラムであって、
コンピュータに、
前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析処理と、
前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出処理と、
前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出処理と、
前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出処理と
を実行させるための画像計測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、葉面積指数を計測するための画像計測方法、画像計測システム、画像計測装置および画像計測プログラムに関する。
【背景技術】
【０００２】
葉面積指数（以下、ＬＡＩ：Ｌｅａｆ  Ａｒｅａ  Ｉｎｄｅｘ）は、農耕地や森林内の植物群落等における単位面積当たりのある方向（例えば、垂直方向）に重なる葉の面積の総和をあらわす。ＬＡＩは、施設栽培等において、植物生育や栽培状況を把握するための一指標として用いられている。【０００３】
ＬＡＩを計測する方法として、例えば、照度センサを用いて間接的に計測する方法が提案されている。この方法では、例えば、植物群落内の上方と下方とで照度を計測し、計測したそれらの照度に基づいてＬＡＩを推定する。【０００４】
また、これに関連する技術として、特許文献１には、植物群落内に設置した反射体を撮影し、その色情報から照度を算出し、ＬＡＩを推定する方法が記載されている。【０００５】
さらに、特許文献２には、ＬＡＩの間接計測方法が記載されている。特許文献２に記載された方法では、間接計測システムは、広角レンズおよび電子式素子を用いて、近赤外光と赤色光とのそれぞれについて、所定領域の画像を撮影する。次いで、間接計測システムは、所定領域を細分した細分領域ごとに、近赤外光と赤色光とのそれぞれについて、輝度値を求める。そして、間接計測システムは、詳細領域ごとの近赤外光と赤色光との輝度値比を求め、輝度値比に基づいて相対日射量を推定し、相対日射量から葉面積指数を求める。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】ＷＯ２０１２／０７３５１９号公報
【特許文献２】特開２００７−１７１０３３号公報
【非特許文献】
【０００７】
【非特許文献１】H. W. Jensen著、「フォトンマッピング」、オーム社出版局、p29,32
【非特許文献２】D. Judd, D. L. MacAdam and G. Wyszecki, "Spectral Distribution of Typical Daylight as a Function of Correlated Color Temperature,"Journal of the Optical Society of America, vol. 54, No. 8, August, 1964, pp.1031-1040
【非特許文献３】E. Kaneko, M. Toda, H. Aoki and M. Tsukada, ''Daylight Spectrum Model under Weather Conditions from Clear Sky to Cloudy",Proc. of IEEE 21st International Conference on Pattern Recognition, 2012.
【非特許文献４】Monsi and Saeki, "On the Factor Light in Plant Communities and its Importance for Matter Production,"Annals of Botany 95, 2005, pp.549-567, Originally published as : Japanese Journal of Botany 14, March,1953, pp.22-52
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１の手法では、光の強度として照度を利用し、ＬＡＩを算出している。しかし、照度とＬＡＩの相関は低いため、算出されたＬＡＩの精度は低い。【０００９】
また、特許文献２の手法では、方向が制限される間接計測システムにより近赤外光と赤色光との輝度値の比からＬＡＩを算出する。しかし、輝度値を用いても、ＬＡＩ算出の精度は低く、さらに光源である太陽光に対して方向が制限されるという課題もある。【００１０】
そこで、本発明は、ＬＡＩを精度良く算出することができる画像計測方法、画像計測システム、画像計測装置および画像計測プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明による画像計測方法は、計測対象である植物の近傍に設置された反射体を撮影し、撮影画像を出力し、前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得し、前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出し、前記絶対分光放射照度分布に基づいて、光量子束密度を算出し、前記光量子束密度に基づいて、葉面積指数を算出することを特徴とする。【００１２】
本発明による画像計測システムは、計測対象である植物の近傍に設置された反射体と、前記反射体を撮影し、撮影画像を出力する撮影手段と、前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析手段と、前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段と、前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段と、前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段とを備えたことを特徴とする。【００１３】
本発明による画像計測装置は、計測対象である植物の近傍に設置された反射体と、前記反射体を撮影し、撮影画像を出力する撮影手段とを備えた画像計測システムに用いられる画像計測装置であって、前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析手段と、前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段と、前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段と、前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段とを備えたことを特徴とする。【００１４】
本発明による画像計測プログラムは、計測対象である植物の近傍に設置された反射体と、前記反射体を撮影し、撮影画像を出力する撮影手段とを備えた画像計測システムに用いられる画像計測プログラムであって、コンピュータに、前記撮影画像内の前記反射体の領域を検出し、前記領域の色情報を取得する画像分析処理と、前記色情報に基づいて、前記領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出処理と、前記絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出処理と、前記光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出処理とを実行させることを特徴とする。
【発明の効果】
【００１５】
本発明によれば、ＬＡＩを精度良く算出することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明によるＬＡＩ画像計測システムを用いた計測系を植物群落に対して正面側から見た正面図である。
【図２】ＬＡＩ画像計測システムを用いた計測系を植物群落に対して上面側から見た上面図である。
【図３】本発明によるＬＡＩ画像計測システムの構成の一例を示すブロック図である。
【図４】輝度値Ｙと放射照度の関係の一例をプロットしたグラフである。
【図５】本発明によるＬＡＩ画像計測システムを用いてＬＡＩを算出する動作の一例を示すフローチャートである。
【図６】本発明によるＬＡＩ画像計測システムの主要な構成を示すブロック図である。
【発明を実施するための形態】
【００１７】
以下、本発明の実施形態を、図面を参照して説明する。
【００１８】
図１は、本発明によるＬＡＩ画像計測システム（葉面積指数画像計測システム）を用いた計測系を植物群落に対して正面側から見た正面図である。図１には、植物群落３０がある断面で切断された断面図が示されている。図２は、ＬＡＩ画像計測システムを用いた計測系を植物群落に対して上面側から見た上面図である。なお、図１に示すように、本実施形態で用いる「垂直方向」は、地面に対して垂直な方向を指し、「水平方向」は、地面に対して水平な方向を指す。
【００１９】
図１に示すように、本実施形態における植物群落３０は、説明を簡単にするために、農作物用の植物が一列に配置された一群の植物群であるものとする。また、植物群落３０は、図２に示すように奥行き方向（図１において正面から奥に向かっていく方向。以下、長手方向ともいう。）に植物が二列に配列された一群の植物群であるが、植物が一列、又は二列よりも多く配置された一群の植物群であってもよい。
【００２０】
なお、本実施形態では、田畑やビニールハウス等の農耕地内の農作物（例えば、トマトやキュウリ）用の植物群落のＬＡＩを計測する用途にＬＡＩ画像計測システムが適用される場合の例を説明する。しかし、本実施形態に示す例に限らず、例えば、森林内の樹木の群落のＬＡＩを計測する用途にＬＡＩ画像計測システムが適用されてもよい。
【００２１】
図１および図２に示すように、植物群落３０の近傍には、垂直方向に複数の反射体１０が配列された配列群が設置されている。また、図２に示すように、カメラ２０は、植物群落３０に遮られることなく反射体１０を撮影可能な位置に設置されている。本例では、図２に示すように、カメラ２０は、反射体１０の正面側に配置される。カメラ２０は、各反射体１０とカメラ２０との間に反射体１０を遮る葉等の障害物がないように配置される。
【００２２】
図１に示す本実施形態の例では、１つの植物群落３０あたり垂直方向に５つの反射体１０が配置され、垂直方向に計１０個の反射体１０が配置されている。しかし、配置可能な反射体１０の数は、本実施形態の例に限られない。例えば、反射体１０が１２個以上配置（１つの植物群落３０あたり垂直方向に反射体１０を６個以上配置）されてもよい。
【００２３】
反射体１０は、具体的には、表面が完全拡散（ランバート）反射の特性を有し、色情報（反射率）が既知である球体によって実現される。反射体１０は、例えば、各種多面体構造、または平らな反射板であってもよい。
【００２４】
本実施形態における反射体１０は、完全白色である（反射率ρ=１０である）。また、反射体１０の色（反射率ρ）は、既知であれば、完全白色に限らず、赤色やグレイ、黄色、橙色であってもよい。ただし、反射体１０の色は、少なくとも、植物群落３０の植物と区別可能にするため、緑色以外の色であることが好ましい。
【００２５】
また、各反射体１０は、植物群落３０の近傍に配置されることによって、植物群落３０によって生じる影等により、上方に配置されたものが最も反射する光の強度が強く、下方に配置されたもの程、反射する光の強度が弱くなる。そのため、垂直方向に複数配置された反射体１０から反射される光の強度を測定し相互に比較することによって、垂直方向のＬＡＩを計測することができ、垂直方向で見た場合の植物群落３０内での葉の重なり具合を把握することができる。なお、植物群落３０の近傍とは、植物群落３０から所定距離以内（例えば、５０ｃｍ以内）のことである。各反射体１０は、垂直方向に配列される際に、各反射体１０から反射される反射光の間に適度な光の強度の強弱が生じるように配置される。
【００２６】
カメラ２０は、具体的には、カラー画像を撮影可能なデジタルカメラ等の撮影装置によって実現される。なお、カメラ２０は、カラー画像を撮影可能なものであれば、静止画像を撮影するものに限らず、例えば、動画像を撮影可能なビデオカメラによって実現されてもよい。また、カメラ２０は、反射体１０を撮影し、撮影画像を後述するＬＡＩ画像計測装置４０（図１および図２において図示せず）に出力する。
【００２７】
なお、本実施形態では、太陽光が照射された反射体１０を撮影し垂直方向のＬＡＩを計測するため、カメラ２０は、球体である反射体１０の上半球を撮影可能な位置（例えば、全ての反射体１０よりも垂直方向に高い位置）に配置されていることが望ましい。
【００２８】
また、図１では、植物群落３０の近傍の１箇所に垂直方向に向かって５つの反射体１０を配置して計測を行う例を示しているが、植物群落３０の長手方向に所定間隔（例えば、５０ｃｍ間隔）で同様に垂直方向に５つずつ反射体１０を配置（すなわち、長手方向の側から見たときに反射体１０を格子状に配置）して計測を行ってもよい。そのようにすれば、植物群落３０の長手方向におけるＬＡＩの分布状況も計測することができる。
【００２９】
図３は、ＬＡＩ画像計測システムの構成の一例を示すブロック図である。図３に示すように、ＬＡＩ画像計測システムは、図２に示した反射体１０を撮影するカメラ２０と、ＬＡＩ画像計測装置４０とを含む。
【００３０】
図３に示すように、ＬＡＩ画像計測装置４０は、画像分析手段４１、絶対分光放射照度分布算出手段４２、データベース４３、ＰＰＦＤ（Ｐｈｏｔｏｓｙｎｔｈｅｔｉｃ  ｐｈｏｔｏｎ  ｆｌｕｘ  ｄｅｎｓｉｔｙ、光量子束密度）算出手段４４、ＬＡＩ（葉面積指数）算出手段４５、データベース４６、及びＬＡＩ出力手段４７を含む。また、図３に示すように、カメラ２０が出力する撮影画像は、ＬＡＩ画像計測装置４０に出力される。なお、ＬＡＩ画像計測装置４０は、具体的には、プログラムに従って動作するパーソナルコンピュータ等の情報処理装置によって実現される。
【００３１】
画像分析手段４１は、具体的には、プログラムに従って動作する情報処理装置のＣＰＵによって実現される。画像分析手段４１は、カメラ２０から入力した撮影画像の三刺激値ＸＹＺ（以下、ＸＹＺ値）を算出する。
【００３２】
画像分析手段４１は、既知の反射体１０の色情報や形状情報に基づいて、撮影画像中から反射体１０が写っている領域を特定（検出）し、特定した領域の輝度値を求める。画像分析手段４１は、具体的には、撮影画像中の白色の物体が写っている領域を特定し、特定した領域の部分画像を抽出する。なお、本実施形態では、白色の反射体１０が複数配置されているので、画像分析手段４１は、複数の領域を特定し、それぞれ特定した領域の部分画像を抽出する。そして、画像分析手段４１は、それぞれ、抽出した部分画像のＸＹＺ値を算出する。
【００３３】
以下では、カメラ２０で撮影したカラー画像中の対象物（反射体１０）のＸＹＺ値を算出する方法の一例を説明する。まず、画像分析手段４１は、入力画像の色情報および形状情報を用いて対象物の領域を検出する。次いで、画像分析手段４１は、対象物の領域の色情報ＲＧＢを取得し、取得した色情報ＲＧＢに基づいて、三刺激値ＸＹＺを計算する。
【００３４】
以下、カラー画像入力装置で撮影されたＲＧＢに関して、このＲＧＢのＲＧＢ蛍光体の色度および白色の色度は、カラー画像入力装置の色特性としてあらかじめ指定されており、ＲＧＢデータは三刺激値ＸＹＺに一意に変換できるものとする。以下では、その変換方法の一例を示す。入力画像のＲＧＢと三刺激値ＸＹＺとの関係は、次式（１）のように定められる。
【００３５】
【数１】
2014103181.tif 000003
【００３６】
式（１）において、Ｍ３×３は、３×３変換行列である。また、ＲＧＢは、ガンマ補正がかかっていないリニアな特性であるものとする。変換式は、ＲＧＢからＸＹＺに一意に変換できればよく、例えば、次の式（２）のように定められてもよい。
【００３７】
【数２】
2014103181.tif 000004
【００３８】
式（２）において、Ｍ３×９は、３×９変換行列である。
【００３９】
なお、Ｍ３×３やＭ３×９は、既知のカラーパッチを用いたカラー画像入力装置のカラーキャリブレーションを実施することで予め算出しておくことが可能である。
【００４０】
絶対分光放射照度分布算出手段４２は、撮影画像内の反射体１０の領域の分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する。以下に示す例では、絶対分光放射照度分布算出手段４２は、まず、反射体１０がランバート反射の特性を有し、反射率ρが既知である物体であることを利用して絶対放射照度を算出する。
【００４１】
Ｉ，Ｒをそれぞれランバート反射の特性を有する物体面上でのある点での放射照度、放射輝度、ρを物体の反射率とする。双方向反射率分布関数（ＢＲＤＦ）をｆＢＲＤＦと表記すると、非特許文献１に記載されているように、ｆＢＲＤＦは、次式（３）のように定義される。
【００４２】
【数３】
2014103181.tif 000005
【００４３】
また、非特許文献１に記載されているように、ランバート反射の特性を有する物体面のｆＢＲＤＦは、次式（４）のように表される。
【００４４】
【数４】
2014103181.tif 000006
【００４５】
本実施形態の例では、反射体１０は完全白色であるため、式（４）においてρ=１．０とすることにより、次式（５）が得られる。
【００４６】
【数５】
2014103181.tif 000007
【００４７】
絶対分光放射照度分布算出手段４２は、撮影時におけるカメラレンズの絞り値Ｆ、シャッタースピード、ゲインの情報を得ることで、式（１）または式（２）から算出される三刺激値Ｙ（輝度値Ｙ）を、画像中の反射体１０の各画素の位置からカメラ方向に放射される絶対放射輝度Ｒ（単位：Ｗ／ｍ２／ｓｒ）に変換する。
【００４８】
この場合、例えば、予め絞り値Ｆ、シャッタースピード、輝度値Ｙ、絶対放射輝度を対応付けたテーブルを用意しておく。絶対分光放射照度分布算出手段４２は、絞り値Ｆ、シャッタースピード、求めた輝度値Ｙを入力し、上述のテーブルを利用することによって、対応する絶対放射輝度Ｒをテーブルから抽出する。また、絶対分光放射照度分布算出手段４２は、抽出した値にさらに補間計算などを用いた計算をしてもよい。
【００４９】
さらに、反射体１０の表面がランバート反射する特性を有するため、式（３）と式（５）を用いると絶対放射輝度Ｒは絶対放射照度Ｉに変換される。
【００５０】
または、絶対分光放射照度分布算出手段４２は、絶対放射照度をテーブルから直接算出してもよい。その場合、予め絞り値Ｆ、シャッタースピード、輝度値Ｙ、絶対放射照度を対応付けたテーブルを用意しておく。絶対分光放射照度分布算出手段４２は、絞り値Ｆ、シャッタースピード、求めた輝度値Ｙに対応する絶対放射照度Ｉをテーブルから抽出する。また、絶対分光放射照度分布算出手段４２は、抽出した値にさらに補間計算などを行ってもよい。
【００５１】
または、絶対分光放射照度分布算出手段４２は、テーブルを用いずに、変換式を用いて絶対放射照度を算出してもよい。この場合、予め絞り値Ｆ、シャッタースピード、輝度値Ｙ、絶対放射輝度または絶対放射照度を対応付けた変換式を定めておく。そして、絶対分光放射照度分布算出手段４２は、この変換式を利用して、絞り値Ｆ、シャッタースピード、求めた輝度値Ｙを絶対放射輝度または絶対放射照度に変換すればよい。
【００５２】
図４は、絞り値Ｆ、シャッタースピード、ゲインを固定した場合の、輝度値Ｙと放射照度の関係を示したグラフである。図からわかるように、これらの値を用いることにより近似の変換式が容易に算出できる。
【００５３】
次に、絶対分光放射照度分布算出手段４２は、式（６）を用いて三刺激値ＸＹＺから色度ｘ，ｙを算出する。
【００５４】
【数６】
2014103181.tif 000008
【００５５】
次に、絶対分光放射照度分布算出手段４２は、色度ｘ，ｙを用いて、物体面からカメラ方向に向かって放射される放射輝度の分光分布（相対分光放射輝度分布）を求める。
【００５６】
以下、色度ｘ，ｙから相対分光放射輝度分布を算出する方法の一例として、非特許文献２に示されたＣＩＥ昼光モデルを利用した例を説明する。絶対分光放射照度分布算出手段４２は、次式（７）に示すように色度ｘ，ｙを用いてＭ１，Ｍ２を算出する。
【００５７】
【数７】
2014103181.tif 000009
【００５８】
λを放射される光の波長とした場合、非特許文献２で与えられるＣＩＥ昼光の平均ベクトルＳ０（λ）、第一主成分ベクトルＳ１（λ）、第二主成分ベクトルＳ２（λ）を用いて、式（８）に示すように相対分光放射輝度分布Ｒ（λ）（式（８）におけるベクトルＲ（チルダ付き）（λ））が得られる。
【００５９】
【数８】
2014103181.tif 000010
【００６０】
この相対分光放射輝度分布Ｒ（λ）は、λ=５６０での値を１とした場合の相対分布である。
【００６１】
反射体１０の表面がランバート反射する特性を有するため、ここで得られた相対分光放射輝度分布Ｒ（λ）は、相対分光放射照度分布Ｉｒ（λ）と等しい。
【００６２】
本例では、絶対分光放射照度分布算出手段４２は、ＣＩＥ昼光モデルを用いて相対分光放射照度分布を算出したが、この例に限らず、例えば非特許文献３に示されたＢｉｒｄモデルなどを用いてもよい。
【００６３】
次に、絶対分光放射照度分布算出手段４２は、得られた絶対放射照度Ｉと相対分光放射照度分布Ｉｒ（λ）から、絶対分光放射照度分布Ｉ（λ）を算出する。具体的には、絶対放射照度Ｉを用いて相対分光放射照度分布Ｉｒ（λ）の各波長での値に定数を乗算し、Ｉ（λ）を算出する。例えば、相対分光放射照度分布Ｉｒ（λ）の４００ｎｍ〜７００ｎｍの積分値が絶対放射照度Ｉと一致するように定数が定められる。そして、絶対分光放射照度分布算出手段４２は、相対分光放射照度分布Ｉｒ（λ）の各波長での値にその定数を乗算し、Ｉ（λ）を算出する。
【００６４】
データベース４３は、具体的には、磁気ディスク装置や光ディスク装置等の記憶装置によって実現される。本実施形態では、データベース４３は、絞り値Ｆ、シャッタースピード、輝度値Ｙを絶対放射輝度または絶対放射照度に変換するための変換テーブルを記憶する。データベース４３が記憶する変換テーブルは、具体的には、絞り値Ｆ、シャッタースピード、輝度値Ｙと絶対放射輝度または絶対放射照度を対応付けて含む。なお、データベース４３が記憶する変換テーブルは、例えば、予め、いくつかのサンプルとなる条件下において、それぞれ、専用装置で計測した絶対放射輝度または絶対放射照度と、絞り値Ｆ、シャッタースピード、撮影画像から求めた輝度値Ｙとを設定することによって構築される。
【００６５】
ＰＰＦＤ算出手段４４は、得られた絶対分光放射照度分布Ｉ（λ）に基づいて、光合成に利用される光の強度（光量子束密度）を算出する。
【００６６】
以下、ＰＰＦＤ算出手段４４が光の強度（光量子束密度）を算出する方法の例として、光合成有効光量子束密度の算出方法を示す。光量子束密度は、単位時間、単位面積を通過する４００ｎｍ〜７００ｎｍの波長の光量子の物質量である（単位：μｍｏｌ／ｍ２／ｓ）。光量子束密度の定義より、放射照度との間に以下の関係が成り立つ。
【００６７】
【数９】
2014103181.tif 000011
【００６８】
光量子１個当たりのエネルギー（Ｊ／個）は、ｈｃ／λと表される。λは、波長、ｈは、プランク定数、ｃは、光速度である。ｈｃ／λと式（９）より、次式（１０）が成り立つ。
【００６９】
【数１０】
2014103181.tif 000012
【００７０】
但し、各記号の意味はそれぞれ以下の通りである。
Ｉ（λ）：波長λの光の放射照度
Ｐ（λ）：波長λの光の光量子束密度
Ｐ：光量子束密度（４００ｎｍ〜７００ｎｍ）
【００７１】
式（１０）を変形することにより、次式（１１）が得られる。
【００７２】
【数１１】
2014103181.tif 000013
【００７３】
また、式（１１）で得られたＰ（λ）を４００ｎｍ〜７００ｎｍまで積分することにより、次式（１２）に示すように光量子束密度Ｐが得られる。
【００７４】
【数１２】
2014103181.tif 000014
【００７５】
ＬＡＩ算出手段４５は、得られた光量子束密度（ＰＰＦＤ）に基づいてＬＡＩ（葉面積指数）を算出する。高さｉにおける反射体１０の位置における光量子束密度Ｐｉと、ある基準となる反射体１０の位置における光量子束密度Ｐ０とを比較し、高さｉの上方に存在する葉のＬＡＩ（これをＦｉと表記する。）を算出する。なお、基準となる反射体１０の位置とは、本例では、最も上部に取り付けられた反射体１０の位置、すなわち、反射体１０と太陽との間に葉による影等の影響が殆どなく、基準の光量子束密度として用いることができる位置とする。
【００７６】
本実施形態では、ＬＡＩ算出手段４５は、データベース４６が記憶するＬＡＩ変換テーブルから、ＰＰＦＤ算出手段４４が算出した高さｉに設置された反射体１０の位置におけるＰＰＦＤと、ある基準となる反射体１０の位置のＰＰＦＤの差分値に対応するＬＡＩを抽出することによって、高さｉでのＬＡＩを求める。
【００７７】
または、ＬＡＩ算出手段４５は、非特許文献４に示された式（１３）を用いた演算処理を行うことによって、ＬＡＩを求めてもよい。
【００７８】
【数１３】
2014103181.tif 000015
【００７９】
式（１３）において、Ｋは、吸収係数を示しており、植物毎に値が異なるとともに、同じ植物であっても天候や時間帯等の外部要因によっても値が異なる。また、Ｆは、積算葉面積指数、すなわちＬＡＩである。
【００８０】
データベース４６は、光量子束密度をＬＡＩに変換するためのＬＡＩ変換テーブルを記憶する。具体的には、データベース４６が記憶するＬＡＩ変換テーブルでは、光量子束密度とＬＡＩとが対応付けられている。データベース４６が記憶するＬＡＩ変換テーブルは、例えば、予め、いくつかのサンプルとなる条件下において、一般に用いられる光量子センサを用いて計測した光量子束密度と、その時に求めたＬＡＩとを設定することによって構築される。
【００８１】
ＬＡＩ出力手段４７は、ＬＡＩ算出手段４５が算出したＬＡＩをコンピュータディスプレイなどに表示する。
【００８２】
次に、ＬＡＩ画像計測システムの動作を説明する。図５は、ＬＡＩ画像計測システムを用いてＬＡＩを計測する動作の一例を示すフローチャートである。本実施形態では、図１および図２に示すように、カメラ２０は、植物群落３０に遮られることなく反射体１０を撮影可能な位置に設置されている。
【００８３】
まず、カメラ２０は、反射体１０が配置された方向を撮影する（ステップＳ１０）。ＬＡＩ画像計測装置４０は、カメラ２０から撮影画像を入力する。
【００８４】
ＬＡＩ画像計測装置４０は、例えば、記憶装置（図示せず）に予め記憶されている反射体１０の色情報（本実施形態では、白色）を取得し、カメラ２０から入力した撮影画像中の白色の物体が写っている領域を特定し、特定した領域の部分画像を抽出する。なお、本実施形態では、白色の反射体１０が複数配置されているので、画像分析手段４１は、複数の領域を特定し、それぞれ特定した領域の部分画像を抽出する。そして、ＬＡＩ画像計測装置４０は、それぞれ、抽出した部分画像の色情報を算出する（ステップＳ１１）。具体的には、画像分析手段４１は、対象物の領域の色情報ＲＧＢを取得し、取得した色情報ＲＧＢに基づいて、三刺激値ＸＹＺを算出する。
【００８５】
また、ＬＡＩ画像計測装置４０は、算出した輝度値Ｙから、放射照度の絶対値である絶対放射照度を算出する（ステップＳ１２）。本実施形態では、ＬＡＩ画像計測装置４０は、データベースが記憶する放射照度変換テーブルから、算出した輝度値Ｙと撮影時の絞り値Ｆ、シャッタースピードに対応する放射照度を算出する。なお、ＬＡＩ画像計測装置４０は、予め定められた輝度値Ｙ、絞り値Ｆ、シャッタースピードを用いた変換式から、算出したＹ値、撮影時の絞りとシャッタースピードに対応する放射照度を算出してもよい。
【００８６】
また、ＬＡＩ画像計測装置４０は、算出したＸＹＺ値から、色度ｘｙ値を算出し、相対分光放射照度分布を算出する（ステップＳ１３）。本実施形態では、ＬＡＩ画像計測装置４０は、ＣＩＥ昼光モデルを用いて色度ｘｙ値から相対分光放射照度分布を算出する。
【００８７】
次に、ＬＡＩ画像計測装置４０は、算出した放射照度の絶対値と相対分光放射照度分布に基づいて、絶対分光放射照度分布を算出する（ステップＳ１４）。本実施形態では、ＬＡＩ画像計測装置４０は、分光放射照度分布の４００〜７００ｎｍの積分値が、算出した絶対放射照度と一致するように定数を定め、相対分光放射照度分布にその定数を乗じて絶対分光放射照度分布を得る。
【００８８】
次に、ＬＡＩ画像計測装置４０は、算出した絶対分光放射照度分布から光量子束密度（ＰＰＦＤ）を算出する（ステップＳ１５）。
【００８９】
次に、ＬＡＩ画像計測装置４０は、算出したＰＰＦＤから葉面積指数（ＬＡＩ）を算出する（ステップＳ１６）。本実施形態では、ＬＡＩ画像計測装置４０は、算出した上下のＰＰＦＤの差分値に対応するＬＡＩを、ＬＡＩ変換テーブルから抽出することによってＬＡＩを算出する。なお、ＬＡＩ画像計測装置４０は、算出したＰＰＦＤから式（１３）に従ってＬＡＩを算出してもよい。
【００９０】
ＬＡＩの算出が終了したら、ＬＡＩ画像計測装置４０は、算出したＬＡＩをディスプレイ等に表示する（ステップＳ１７）。
【００９１】
本実施形態のＬＡＩ画像計測システムは、ＬＡＩとの相関が高い光量子束密度を算出し、光量子束密度に基づいてＬＡＩを算出するので、より精度良くＬＡＩを算出することができる。
【００９２】
図６は、本発明による画像計測システムの主要な構成を示すブロック図である。図６に示すように本発明による画像計測システムは、主要な構成として、計測対象である植物３１の近傍に設置された反射体１０と、反射体１０を撮影し、撮影画像を出力する撮影手段２１と、撮影画像内の反射体１０の領域を検出し、領域の色情報を取得する画像分析手段４１と、色情報に基づいて、領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段４２と、絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段４４と、光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段４５とを備える。
【００９３】
また、上記の実施形態では、以下の（１）〜（５）に示すような画像計測システムおよび画像計測装置も開示されている。
【００９４】
（１）計測対象である植物（例えば、植物群落３０のうちの一部の植物）の近傍に設置された反射体（例えば、反射体１０）と、反射体を撮影し、撮影画像を出力する撮影手段（例えば、カメラ２０）と、撮影画像内の反射体の領域を検出し、領域の色情報を取得する画像分析手段（例えば、画像分析手段４１）と、色情報に基づいて、領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段（例えば、絶対分光放射照度分布算出手段４２）と、絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段（例えば、ＰＰＦＤ算出手段４４）と、光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段（例えば、ＬＡＩ算出手段４５）とを備えたことを特徴とする画像計測システム（例えば、ＬＡＩ画像計測システム）。
【００９５】
（２）画像計測システムは、反射体が、色が既知である完全拡散反射面を有するように構成されていてもよい。
【００９６】
（３）画像計測システムは、光量子束密度と葉面積指数とを対応付けて記憶する記憶手段（例えば、データベース４６）を備え、葉面積指数算出手段が、算出された光量子束密度に対応する葉面積指数を記憶手段から抽出することによって葉面積指数を算出するように構成されていてもよい。
【００９７】
（４）画像計測システムは、絶対分光放射照度分布算出手段が、色情報に基づいて、領域の輝度値を含む三刺激値を算出し、輝度値に基づいて、領域の放射照度の絶対値である絶対放射照度を算出し、三刺激値に基づいて、領域の色度を算出し、色度に基づいて、領域における分光放射照度の相対分布である相対分光放射照度分布を算出し、絶対放射照度と相対分光放射照度分布とに基づいて、絶対分光放射照度分布を算出するように構成されていてもよい。
【００９８】
（５）計測対象である植物（例えば、植物群落３０のうちの一部の植物）の近傍に設置された反射体（例えば、反射体１０）と、反射体を撮影し、撮影画像を出力する撮影手段（例えば、カメラ２０）とを備えた画像計測システムに用いられる画像計測装置（例えば、ＬＡＩ画像計測装置４０）であって、撮影画像内の反射体が撮影された領域を検出し、領域の色情報を取得する画像分析手段（例えば、画像分析手段４１）と、色情報に基づいて、領域における分光放射照度の絶対値の分布である絶対分光放射照度分布を算出する絶対分光放射照度分布算出手段（例えば、絶対分光放射照度分布算出手段４２）と、絶対分光放射照度分布に基づいて、光量子束密度を算出する光量子束密度算出手段（例えば、ＰＰＦＤ算出手段４４）と、光量子束密度に基づいて、葉面積指数を算出する葉面積指数算出手段（例えば、ＬＡＩ算出手段４５）とを備えたことを特徴とする画像計測装置。
【００９９】
この出願は、２０１２年１２月２６日に出願された日本出願特願２０１２−２８２０５６を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【０１００】
以上、実施形態を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
【産業上の利用の可能性】
【０１０１】
本発明は、施設栽培等において、植物生育や栽培状況を把握する用途に適用される。
【符号の説明】
【０１０２】
１０  反射体
２０  カメラ
２１  撮影手段
３０  植物群落
３１  植物
４０  ＬＡＩ画像計測装置
４１  画像分析手段
４２  絶対分光放射照度分布算出手段
４３，４６  データベース
４４  光量子束密度（ＰＰＦＤ）算出手段
４５  葉面積指数（ＬＡＩ）算出手段
【図１】
2014103181.tif 000016
【図２】
2014103181.tif 000017
【図３】
2014103181.tif 000018
【図４】
2014103181.tif 000019
【図５】
2014103181.tif 000020
【図６】
2014103181.tif 000021
【国際調査報告】
2014103181.tif 000022
2014103181.tif 000023
2014103181.tif 000024
2014103181.tif 000025