(57)【要約】      （修正有）【課題】森林内と森林外との日射量を同時測定しなくても、さらには森林外との日射量の測定自体をなくしても、高精度に葉面積指数を求め得るようにする。【解決手段】葉面積指数を間接的に測定するに際し、広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影し、前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求め、前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求め、前記輝度値比を用いて相対日射量を推定し、前記相対日射量から葉面積指数を求めるものである。【選択図】図２
2007171033.tif 000002
【特許請求の範囲】
【請求項１】
葉面積指数を間接的に測定するに際し、
広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影し、
前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求め、
前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求め、
前記輝度値比を用いて相対日射量を推定し、
前記相対日射量から葉面積指数を求めることを特徴とする、葉面積指数の間接測定方法。
【請求項２】
前記細分領域は、前記電子式撮像素子の一または複数の画素に対応するものである、請求項１記載の葉面積指数の間接測定方法。
【請求項３】
前記所定領域の、天頂から所定天頂角までを天頂角によって三つの大領域に分割して、前記大領域毎に前記輝度値比の平均値を求め、
前記大領域毎の輝度値比の平均値から前記相対日射量としてのギャップの割合を推定し、
前記ギャップの割合から葉面積指数を求めることを特徴とする、請求項１または２記載の葉面積指数の間接測定方法。
【請求項４】
広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影する撮影手段と、
前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求める画像選択・輝度値演算手段と、
前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求める輝度値比演算手段と、
前記輝度値比を用いて相対日射量を推定する相対日射量推定手段と、
前記相対日射量から葉面積指数を求める葉面積指数演算手段と、
を具えてなる、葉面積指数の間接測定システム。
【請求項５】
前記細分領域は、前記電子式撮像素子の一または複数の画素に対応するものである、請求項４記載の葉面積指数の間接測定システム。
【請求項６】
前記輝度値比演算手段は、前記所定領域の、天頂から所定天頂角までを天頂角によって三つの大領域に分割して、前記大領域毎に前記輝度値比の平均値を求め、
前記相対日射量推定手段は、前記大領域毎の輝度値比の平均値から前記相対日射量としてのギャップの割合を推定し、
前記葉面積指数演算手段は、前記ギャップの割合から葉面積指数を求めることを特徴とする、請求項４または５記載の葉面積指数の間接測定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、森林のバイオマス量の代表的な指標値である葉面積指数を間接的に測定する方法およびシステムに関するものである。
【背景技術】
【０００２】
葉面積指数（LAI：Leaf Area Index）とは、単位水平面積当りの地上部の全ての葉の片面の面積の合計値のことであり、例えばLAI=２なら、１０ｍ２の面積内の葉の片面の面積の合計が２０ｍ２となって、水平な地面に葉を隙間なく並べれば２枚重なると解釈できる。
【０００３】
森林のバイオマス量の代表的な指標値である上記葉面積指数は、森林の熱収支モデル、蒸発散モデル、光合成モデル等の主要パラメータであり、地球温暖化研究の全球モデルの植生表現の主要パラメータにも採用されている。また森林に関わる様々な研究の試験地の説明項目としても使用されている。
【０００４】
かかる葉面積指数は、全木伐採や標準木伐採によって直接測定することができるが、全木伐採は、伐採に対する抵抗感および、季節変化やプロセスの重視のため近年は少なくなっており、標準木伐採は、種構成の多様な森林では困難であった。
【０００５】
一方、葉面積指数は相対日射量の関数となることが知られており（非特許文献１参照）、これを利用した葉面積指数の間接測定も実用化され、普及している（例えば、非特許文献２参照）。
【０００６】
また、魚眼レンズを用い、天頂角・方位角で分割した領域毎の相対照度を測定して高精度に葉面積指数を求めることも行われている（非特許文献３参照）。
【非特許文献１】Monsi & Saeki, 1953, Jpn. J. Bot. 14, 22-52.
【非特許文献２】J. M. Welles and J. M. Norman 1991, Instrument for indirect measurement of canopy architecture, Agronomy Journal, 83.818-825.
【非特許文献３】J. M. Norman & G. S. Campbell, Canopy structure, p.301-325, in R. W. Pearcy, J. Ehleringer, H. A. Mooney and P. W. Rundel (Eds.) Plant Physiological Ecology, 1989, Chapman & Hall, London, pp.457.
【発明の開示】
【発明が解決しようとする課題】
【０００７】
ところで、相対日射量を得るためには、森林内と森林外との日射量を同時測定することが必要であるが、樹高も空間的な広がりも大きな森林では、森林内外の日射量の同時測定は容易でないことが多いため、実際には時間差を置いて測定されることが多く、大きな誤差要因となっていた。
【０００８】
すなわち、時刻や天候の変化で森林外の日射量は数オーダー変動することもあり、例えば２オーダー異なれば、葉面積指数の値が±４以上も異なってしまい、これは到底無視できない程度であって、葉面積指数の測定値の信頼性を大きく損ねていた。
【０００９】
そこでこの状況を改善するため従来、日射が森林の林冠部を通過するとき近赤外と赤色光とで減衰の程度が異なり、森林内の日射は森林外の日射に比べ近赤外光／赤色光の光量比が大きくなるという点に着目し、この比を用いて相対日射量を推定して葉面積指数を求めることも提案されていたが、結局この手法は失敗して放棄されてしまった（非特許文献３参照）。
【００１０】
その失敗の理由は、従来は画像から光量を測定していたわけではないが、その測定は結局（近赤外光の積算値）／（赤色光の積算値）で計算していたことに相当するため、葉群の少ない部位（すなわち輝度の高い部位）の情報の重みが大きくなり、葉群の多い部位（すなわち輝度の低い部位）の情報の重みが低くなったためと考えられる。
【課題を解決するための手段】
【００１１】
この発明は、上記従来の課題を有利に解決することを目的とするものであり、この発明の葉面積指数の間接測定方法は、葉面積指数を間接的に測定するに際し、広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影し、前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求め、前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求め、前記輝度値比を用いて相対日射量を推定し、前記相対日射量から葉面積指数を求めることを特徴とするものである。
【００１２】
また、この発明の葉面積指数の間接測定システムは、広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影する撮影手段と、前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求める画像選択・輝度値演算手段と、前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求める輝度値比演算手段と、前記輝度値比を用いて相対日射量を推定する相対日射量推定手段と、前記相対日射量から葉面積指数を求める葉面積指数演算手段と、を具えてなるものである。
【発明の効果】
【００１３】
上記間接測定方法によれば、細分領域毎に輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択し、その選択した画像を所定露出時間で撮像したものとした場合のその画像の輝度値を求めるので、各細分領域の輝度値が極めて多段階の階調値で高精度に表現でき、しかも所定領域としての全天等を細分した細分領域毎にその高精度の輝度値を求めるので、葉群の少ない部位（すなわち輝度の高い部位）の情報と葉群の多い部位（すなわち輝度の低い部位）の情報との重みの差をほとんど若しくは全くなくすことができ、その高精度の輝度値を用いて近赤外光と赤色光との輝度値比を求めて相対日射量を推定し、そこから葉面積指数を求めるから、森林内と森林外との日射量を同時測定しなくても、さらには森林外との日射量の測定自体をなくしても、高精度に葉面積指数を求めることができる。
【００１４】
なお、この発明の葉面積指数の間接測定方法においては、前記細分領域は、前記電子式撮像素子の一または複数の画素に対応するものであっても良く、このようにすれば、各細分領域の輝度値を容易に求めることができる。
【００１５】
さらに、この発明の葉面積指数の間接測定方法においては、前記所定領域の、天頂から所定天頂角までを天頂角によって三つの大領域に分割して、前記大領域毎に前記輝度値比の平均値を求め、前記大領域毎の輝度値比の平均値から前記相対日射量としてのギャップの割合を推定し、前記ギャップの割合から葉面積指数を求めることとしても良く、このようにすれば、前記非特許文献２によれば天頂から所定天頂角までを天頂角によって三つに分割した領域毎の近赤外光と赤色光との光量比と市販のLAI測定器である後述のLAI-2000で測定したギャップの割合との間に一価関係があることが判明しているので、この非特許文献２に倣って、天頂から所定天頂角までを天頂角によって三つに分割した大領域毎の輝度値比の平均値と市販のLAI測定器のギャップの割合との間の一価関係を用いて前記輝度値比の平均値をギャップの割合に変換し、そのギャップの割合から非特許文献２に倣って計算してLAI値を求めることができる。
【００１６】
一方、上記間接測定システムによれば、撮影手段が、広角レンズおよび電子式撮像素子を用いて、近赤外光と赤色光とのそれぞれについて、露出時間を複数種類に変えて所定領域の画像を撮影し、画像選択・輝度値演算手段が、前記所定領域を細分した細分領域毎に、近赤外光と赤色光とのそれぞれについて、前記複数種類の露出時間の画像中で輝度値が露出不足から露出過多までの間の階調の中央値に最も近い画像を選択して、その選択した画像を所定露出時間に正規化した場合の輝度値を求め、輝度値比演算手段が、前記細分領域毎に求めた近赤外光と赤色光との輝度値から、前記細分領域毎の近赤外光と赤色光との輝度値比を求め、そして葉面積指数演算手段が、前記輝度値比を用いて相対日射量を推定する相対日射量推定手段と、前記相対日射量から葉面積指数を求めることから、上記間接測定方法を実施し得て、各細分領域の輝度値が極めて多段階の階調値で高精度に表現でき、しかも所定領域としての全天等を細分した細分領域毎にその高精度の輝度値を求めるので、葉群の少ない部位（すなわち輝度の高い部位）の情報と葉群の多い部位（すなわち輝度の低い部位）の情報との重みの差をほとんど若しくは全くなくすことができ、その高精度の輝度値を用いて近赤外光と赤色光との輝度値比を求めて相対日射量を推定し、そこから葉面積指数を求めるから、森林内と森林外との日射量を同時測定しなくても、さらには森林外との日射量の測定自体をなくしても、高精度に葉面積指数を求めることができる。
【００１７】
なお、この発明の葉面積指数の間接測定システムにおいては、前記細分領域は、前記電子式撮像素子の一または複数の画素に対応するものであっても良く、このようにすれば、各細分領域の輝度値を容易に求めることができる。
【００１８】
さらに、この発明の葉面積指数の間接測定システムにおいては、前記輝度値比演算手段は、前記所定領域の、天頂から所定天頂角までを天頂角によって三つの大領域に分割して、前記大領域毎に前記輝度値比の平均値を求め、前記相対日射量推定手段は、前記大領域毎の輝度値比の平均値から前記相対日射量としてのギャップの割合を推定し、前記葉面積指数演算手段は、前記ギャップの割合から葉面積指数を求めても良く、このようにすれば、前記非特許文献２によれば天頂から所定天頂角までを天頂角によって三つに分割した領域毎の近赤外光と赤色光との光量比と市販のLAI測定器である後述のLAI-2000で測定したギャップの割合との間に一価関係があることが判明しているので、この非特許文献２に倣って、天頂から所定天頂角までを天頂角によって三つに分割した大領域毎の輝度値比の平均値と市販のLAI測定器のギャップの割合との間の一価関係を用いて前記輝度値比の平均値をギャップの割合に変換し、そのギャップの割合から非特許文献２に倣って計算してLAI値を求めることができる。
【発明を実施するための最良の形態】
【００１９】
以下、本発明の実施の形態を実施例によって、図面に基づき詳細に説明する。ここに、図１は、この発明の葉面積指数の間接測定方法の一実施例の実施に用いる、この発明の葉面積指数の間接測定システムの一実施例の構成を示す説明図、図２は、上記実施例の葉面積指数の間接測定方法の手順を示すフローチャート、そして図３（ａ）および（ｂ）は、上記実施例の方法で求めた画素毎のIR／赤色光の平均値とギャップの割合との関係および従来のIR累積値／赤色光累積値の平均値とギャップの割合との関係をそれぞれ示す関係図である。
【００２０】
図１に示す上記実施例の間接測定システムは、三脚１上に固定されたデジタルカメラ２と、そのデジタルカメラ２に装着された広角レンズとしての魚眼レンズ３と、その魚眼レンズ３に装着された光学フィルタ４と、デジタルカメラ２に制御命令を送るとともにデジタルカメラ２からその撮影した画像データを取り込んで記録するパーソナルコンピュータ（パソコン）５と、そのパソコン５に電源を供給するバッテリ６とを具えてなる。
【００２１】
ここで、デジタルカメラ２は可視光と近赤外光との両方に感度をもつ必要があり、このようなデジタルカメラ２としては、例えばソニー株式会社製IEEE1394カメラを使用できる。また魚眼レンズ３としては、例えば株式会社タムロン製1/2インチ用単焦点、画角104.0×77.0を使用できる。さらに、光学フィルタ４としては、赤色光撮影時は例えばシグマ光機株式会社製近赤外吸収フィルタおよび同社製シャープカットフィルタを二枚重ねで使用でき、近赤外光撮影時は同社製赤外透過フィルタを使用できる。そしてデジタルカメラ２とパソコン５との間のインターフェースとしては、例えばラトックシステム株式会社製IEEE1394PC-Cardを使用することができる。
【００２２】
かかる実施例の間接測定システムを用いて葉面積指数の間接測定を行うに際しては、図２に示すように、先ずステップＳ１〜Ｓ３で撮影を行う。すなわち、ステップＳ１で、測定対象とする林冠下に三脚１を据え、水準器を用いてデジタルカメラ２の視軸が鉛直上方を向くようにデジタルカメラ２を固定する。デジタルカメラ２の絞りは、続くステップＳ２，Ｓ３でシャッタ速度を例えば1/100,000秒から0.4秒まで順次に変えて撮影したときに露出状態が露出不足から露出過多まで画像が変化するように適宜調整しておく。
【００２３】
次にステップＳ２で、上記赤外透過フィルタ等の近赤外光のみが透過する光学フィルタ４を魚眼レンズ３の前方に装着して、パソコン５からの制御命令でシャッタ速度（露出時間）を例えば1/100,000秒、1/50,000秒、1/20,000秒、0.000117秒、0.000284秒、0.000535秒、0.001037秒、0.00205秒、0.00413秒、0.00831秒、0.016669秒、0.033305秒、0.066995秒、0.133333秒、0.2秒、0.4秒の１６種類に順次変えながらデジタルカメラ２で全天画像を１６枚撮影し、それらの画像データを全てパソコン５に転送して保存する。なお、この撮影の際、シャッタ速度（露出時間）の変更、撮影の指令、画像データのデジタルカメラ２からパソコン５への転送保存等は全て、パソコン５上にあらかじめインストールしたプログラムにより全自動で処理される。
【００２４】
次にステップＳ３で、上記近赤外吸収フィルタおよびシャープカットフィルタの二枚重ね等の赤色光のみが透過する光学フィルタ４を魚眼レンズ３の前方に装着して、ステップＳ２と同様に、パソコン５からの制御命令でシャッタ速度（露出時間）を例えば1/100,000秒、1/50,000秒、1/20,000秒、0.000117秒、0.000284秒、0.000535秒、0.001037秒、0.00205秒、0.00413秒、0.00831秒、0.016669秒、0.033305秒、0.066995秒、0.133333秒、0.2秒、0.4秒の１６種類に順次変えながらデジタルカメラ２で全天画像を１６枚撮影し、それらの画像データを全てパソコン５に転送して保存する。なお、この撮影の際も、シャッタ速度（露出時間）の変更、撮影の指令、画像データのデジタルカメラ２からパソコン５への転送保存等は全て、パソコン５上にあらかじめインストールしたプログラムにより全自動で処理される。
【００２５】
従って、上記デジタルカメラ２、魚眼レンズ３および光学フィルタ４と、上記ステップＳ２，Ｓ３を実行する上記パソコン５とは、撮影手段に相当し、上記のようにして撮影した画像は、例えば1024×768個の画素を持ち、各画素は、青・緑・赤の三バンドについてそれぞれ例えば256階調の輝度値を持つ。そして、近赤外光と赤色光とは共に、赤バンドの輝度値として捉えられる。
【００２６】
次いでステップＳ４〜Ｓ８で計算を行う。なお、この計算も、パソコン５上にあらかじめインストールしたプログラムにより全自動で処理される。
【００２７】
すなわち、ステップＳ４で、ステップＳ２で撮影した例えば上記１６枚の近赤外光画像から、その画像の全画素（例えば1024×768個の画素）について細分領域としての一画素毎に、画素の輝度値が階調の中央値（例えば256階調の場合は128）に最も近い画像を選び出し、選んだ画像の露出時間を用い、その画像を例えば1/100秒で正規化した場合の輝度値を計算する。この処理により近赤外光画像の輝度値は７オーダー程度の階調値を表現できる。
【００２８】
次いでステップＳ５で、ステップＳ４と同様に、ステップＳ３で撮影した例えば上記１６枚の赤色光画像から、その画像の全画素（例えば1024×768個の画素）について細分領域としての一画素毎に、画素の輝度値が階調の中央値（例えば256階調の場合は128）に最も近い画像を選び出し、選んだ画像の露出時間を用い、その画像を例えば1/100秒で正規化した場合の輝度値を計算する。この処理により赤色光画像の輝度値も７オーダー程度の階調値を表現できる。従って、上記ステップＳ４，Ｓ５を実行する上記パソコン５は、画像選択・輝度値演算手段に相当する。
【００２９】
次いでステップＳ６で、上記近赤外光画像および赤色光画像の全画素（例えば1024×768個の画素）について一画素毎に、近赤外光の正規化輝度値と赤色光の正規化輝度値との比（近赤外光正規化輝度値／赤色光正規化輝度値）を計算して求め、天頂角（zenith angle：鉛直真上を０度とする）で分けた三つの大領域（領域１：０〜14度，領域２：16〜30度，領域３：31〜45度）について大領域毎の上記輝度値比の平均値を計算する。従って、上記ステップＳ６を実行する上記パソコン５は、輝度値比演算手段に相当する。
【００３０】
ところで、前述のように非特許文献２によれば、天頂から所定天頂角までを天頂角によって三つに分割した領域毎の近赤外光と赤色光との光量比と市販のLAI測定器であるLI-COR社製LAI測定器LAI-2000で測定したギャップの割合との間には一価関係があることが判明しているので、上記輝度値比の平均値と市販のLAI測定器である例えば上記LAI-2000で測定したギャップの割合との一価的関係をあらかじめ求めておき、この関係を利用して次のステップＳ７で、領域毎の上記輝度値比の平均値をギャップの割合に変換し、続くステップＳ８で、それらのギャップの割合からLAIの値を非特許文献２に倣って計算して求める。従って、上記ステップＳ７を実行する上記パソコン５は、相対日射量推定手段に相当し、また上記ステップＳ８を実行する上記パソコン５は、葉面積指数演算手段に相当する。
【００３１】
すなわち、本願発明者が上記の如くして正規化輝度値を求めた近赤外光画像および赤色光画像とそれらの画素毎の輝度値比を輝度で表した画像とを対比してみると、赤色光画像も近赤外光画像も何れも、葉の少ない領域が明るく、幹は暗くなっているが、両者の輝度値の比をとった画像では、葉の少ない領域も幹も近赤外光正規化輝度値／赤色光正規化輝度値比が低い（画像では黒っぽい）のに対し、葉の多い領域は近赤外光正規化輝度値／赤色光正規化輝度値比が高い値（画像では明るい）を示すことがわかった。これは、近赤外光正規化輝度値／赤色光正規化輝度値比が高い領域では、葉群によって日射の透過が妨げられる割合が大きいからであると推定される。
【００３２】
図３（ａ）は、上記実施例の方法で求めた1024×768個の画素毎の近赤外光正規化輝度値／赤色光正規化輝度値比の上記三つの大領域についての大領域毎の平均値（横軸）と市販のLAI測定器（LAI-2000）のギャップの割合（縦軸）との関係を示す図であり、例えば図中直線で示すように、近赤外光／赤色光の輝度値比の大領域毎の平均値と市販のLAI測定器である上記LAI-2000で測定したギャップの割合との間には一価的な関係が見出せる。従って、この関係を利用して上記ステップＳ７，Ｓ８のようにギャップの割合を求め、そのギャップの割合からLAIの値を、非特許文献２に倣って、以下の如くして計算することができる。
【００３３】
すなわち、天頂角で区切った領域（以下、リングと呼ぶ）のLAI値は、そのリングのギャップの割合から、次式（１）で求められることが知られている（非特許文献２）。
「あるリングのLAI値」=（−１／ｋ）×ln(「そのリングのギャップの割合」)×cos(「そのリングの天頂角代表値」)・・・（１）
ここで、ギャップの割合は、相対照度として測定される。ｋは、生態学の分野において吸光係数と呼ばれ、葉の傾斜角に対応する指数であり、0.5が採用されることが多い。
【００３４】
なお、非特許文献２に倣って、リングの区切りを、リング１：0.0〜12.3度、リング２：16.7〜28.6度、リング３：32.4〜43.4度とし、リングの天頂角代表値を、7度、23度、38度としてもよい。非特許文献２に倣えば、この３つのリングのLAI値を平均する場合は、0.034：0.104：0.160の比で重みをつければよい。すなわち、吸光係数ｋを0.5とするとともに重みの値を0.114、0.349、0.537として平均LAI値を求めることができる。
【００３５】
これに対し図３（ｂ）は、画素毎の処理をしていない近赤外光平均値／赤色光平均値の比（横軸）とギャップの割合（縦軸）との間の関係を示しているが、この図では、図示のようにそれらの関係はあいまいである。この図における横軸は、カメラを用いていない先人達の研究に対応すると考えられる。
【００３６】
従って、この実施例の葉面積指数の間接測定方法および間接測定システムによれば、森林内と森林外との日射量を同時測定しなくても、さらには森林外との日射量の測定自体をなくしても、高精度に葉面積指数を求めることができる。しかもこの実施例の方法によれば、従来は測定困難だった針葉樹の葉や斑入り広葉樹の葉も、十分良好に測定することができる。
【００３７】
以上、図示例に基づき説明したが、この発明は上述の例に限定されるものでなく、特許請求の範囲の記載範囲内で適宜変更することができ、例えば、上記実施例の方法およびシステムでは細分領域を一画素として正規化輝度値を求めているが、その代わりに纏まった複数画素としても良い。
【００３８】
また、例えば上記実施例の間接測定システムのデジタルカメラの鉛直上方に葉を挟み込むチャンバ部を配置して、撮影する所定領域内にその葉が入るようにすることで、葉の形態や大きさの測定の標準化に用いることができ、さらに、そのチャンバ部の上方に人工光源を配置することで、葉の赤色光透過特性／近赤外光透過特性（SPAD値）ひいては葉内クロロフィル濃度の間接測定による樹種や樹齢による光合成特性の比較等に用いることができ、この場合に対象物を従来は測定困難だった針葉樹の葉や斑入り広葉樹の葉としても、十分良好に測定することができる。
【００３９】
さらに、上記チャンバ部の上方でなく下方にデジタルカメラと並置して人工光源を配置することで、リモートセンシングによる植物調査におけるグランドトュルース（地上検証測定）にも用いることができ、この場合に対象物を従来は測定困難だった針葉樹の葉や斑入り広葉樹の葉としても、十分良好に測定することができる。
【００４０】
そして上記のような室内で用いる間接測定システムでは、パーソナルコンピュータ（パソコン）はバッテリ駆動でなく商用電源駆動としても良い。また、デジタルカメラの制御および画像データの保存を行うコンピュータと、計算処理を行うコンピュータとは別個にしても良い。
【産業上の利用可能性】
【００４１】
かくしてこの発明の葉面積指数の間接測定方法および間接測定システムによれば、森林内と森林外との日射量を同時測定しなくても、さらには森林外との日射量の測定自体をなくしても、高精度に葉面積指数を求めることができる。
【図面の簡単な説明】
【００４２】
【図１】この発明の葉面積指数の間接測定方法の一実施例の実施に用いる、この発明の葉面積指数の間接測定システムの一実施例の構成を示す説明図である。
【図２】上記実施例の葉面積指数の間接測定方法の手順を示すフローチャートである。
【図３】（ａ）および（ｂ）は、上記実施例の方法で求めた画素毎のIR／赤色光の平均値とギャップの割合との関係および従来のIR累積値／赤色光累積値の平均値とギャップの割合との関係をそれぞれ示す関係図である。
【符号の説明】
【００４３】
１  三脚
２  デジタルカメラ
３  魚眼レンズ
４  光学フィルタ
５  パーソナルコンピュータ（パソコン）
６  バッテリ
【図１】
2007171033.tif 000003
【図２】
2007171033.tif 000004
【図３】
2007171033.tif 000005