(57)【要約】【課題】  孤立した樹体または間隔をおいた樹体について、樹冠構造測定装置または魚眼レンズを装着した撮像装置を用いて非破壊的な手法により、多くの時間を要さずに樹体生産力を評価できるようにする。【解決手段】  樹体の外形を表すものとして回転半長円体モデルを決定し、主幹から所定の距離の位置に光学的樹体構造測定装置を配置して測定した樹葉を透過した光の強度と、樹葉を透過しない光の強度とを測定し、半長円体モデルにより樹体への入射した光が通過する光路長を計算し、測定された光の強度と、光路長とから樹体総葉面積を求める。光学的樹体構造測定装置の代わりに魚眼レンズを装着した撮像装置を配置し、撮像により取得された画像データについての演算処理によっても樹体総葉面積が求められる。【選択図】  図４
2007069736.tif 000002
【特許請求の範囲】
【請求項１】
孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について樹体の透過光を測定することにより樹体生産能力を評価する方法であって、
１本の樹体の主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体とする回転半長円体モデルを決定することと、
該主幹から所定の距離だけ離れた地面近辺の位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して該光学的樹体構造測定装置により測定対象となる複数の天頂角で入射し樹葉を通過した光の強度を測定することと、
天空からの光が遮られない位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して該光学的樹体構造測定装置により測定対象となる前記複数の天頂角での樹葉を通過しない光の強度を測定することと、
前記複数の天頂角で入射し樹葉を通過した光の強度と樹葉を通過しない光の強度とから該複数の天頂角での空隙率を求めることと、
測定対象とする前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積を回転半長円体モデルに基づいて計算することと、
前記複数の天頂角での空隙率と前記複数の天頂角に対応する光路長とから各天頂角に対応する光の減衰量を求めることと、
該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求めることと、
からなることを特徴とする樹体生産能力を評価する方法。
【請求項２】
前記測定された複数の天頂角での入射光の強度から、次式（１）
【数１】
2007069736.tif 000021
により該複数の天頂角での空隙率Ｔ（θｉ）を求め、次式（２）
【数２】
2007069736.tif 000022
により葉面積密度（ＬＡＤ）を求めるようにしたことを特徴とする請求項１に記載の樹体生産能力を評価する方法。
【請求項３】
孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について、樹体の透過光を測定することにより樹体生産能力を評価する方法であって、
１本の樹体の主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体とする回転半長円体モデルを決定することと、
該主幹から所定の距離だけ離れた地面近辺の位置に魚眼レンズを装着した撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置することと、
前記撮像装置により測定対象となる複数の天頂角からの入射光を含む画角の画像を撮像することと、
前記撮像装置により取得された画像データについて前記複数の天頂角での入射光の開空度を求める演算処理を行うことと、
測定対象とする前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積を回転半長円体モデルに基づいて計算することと、
前記複数の天頂角からの入射光の開空度と前記複数の天頂角に対応する光路長とから各天頂角に対応する光の減衰量を求めることと、
該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求めることと、
からなることを特徴とする樹体生産能力を評価する方法。
【請求項４】
前記撮像装置により撮像された画像について、次式（５）
【数５】
2007069736.tif 000023
により前記複数の天頂角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000024
により葉面積密度（ＬＡＤ）を求めるようにしたことを特徴とする請求項３に記載の樹体生産能力を評価する方法。
【請求項５】
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしたことを特徴とする請求項１〜４のいずれか１項に記載の樹体生産能力を評価する方法。
【請求項６】
魚眼レンズと、該魚眼レンズの結像位置に配設された撮像素子と、撮像動作制御部と、前記撮像素子により取得された画像データについて演算処理を行う演算回路と、該演算回路による演算結果を出力する出力部とを備えてなる撮像装置であって、前記演算回路は、孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について樹体の透過光を測定する際に、１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に前記撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置して撮像し取得された画像について測定対象となる複数の天頂角をそれぞれ含む複数の角度範囲の画像部分におる画像データから前記複数の角度範囲での開空度を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記撮像装置まで通過する光路長及び樹冠体積を求め、前記開空度と前記複数の天頂角に対する光路長とから各天頂角を含む角度範囲に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うものであることを特徴とする樹体生産能力を評価するための撮像装置。
【請求項７】
前記演算回路は取得された画像データについて、次式（５）
【数５】
2007069736.tif 000025
により前記複数の天長角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000026
により葉面積密度（ＬＡＤ）を求める演算処理を行うようにしたものであることを特徴とする請求項６に記載の樹体生産能力を評価するための撮像装置。
【請求項８】
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしたことを特徴とする請求項６または７のいずれか１項に記載の樹体生産能力を評価するための撮像装置。
【請求項９】
孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について該樹体の透過光を測定する際に、１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して測定された複数の天頂角で入射し樹葉を通過した光の強度と、天空からの光が遮られない位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して測定された前記複数の天頂角での樹葉を通過しない光の強度とから該複数の天頂角での空隙率を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて求められた前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積とから各天頂角に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うことにより樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
【請求項１０】
前記測定された複数の天頂角での入射光の強度から、次式（１）
【数１】
2007069736.tif 000027
により該複数の天頂角での空隙率Ｔ（θｉ）を求め、次式（２）
【数２】
2007069736.tif 000028
により葉面積密度（ＬＡＤ）を求めるようにしたことを特徴とする請求項９に記載の樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
【請求項１１】
孤立して存在する樹体または間隔をおいて複数存在する樹体について樹体の透過光を測定する際に１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に魚眼レンズを装着した撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置して撮像し取得された画像について測定対象となる複数の天頂角をそれぞれ含む複数の角度範囲の画像部分におる画像データから前記複数の角度範囲での開空度を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記撮像装置まで通過する光路長及び樹冠体積を求め、前記開空度と前記複数の天頂角に対応する光路長とから各天頂角を含む角度範囲に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うことにより樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
【請求項１２】
複数の天頂角からの入射光を含む画角について取得された画像データについて、次式（５）
【数５】
2007069736.tif 000029
により前記複数の天長角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000030
により葉面積密度（ＬＡＤ）を求めるようにしたことを特徴とする請求項１１に記載の樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
【請求項１３】
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしたことを特徴とする請求項９〜１２のいずれか１項に記載の樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、樹体生産能力を評価する方法、樹体生産能力を評価するための撮像装置及び樹体生産能力を評価するためのプログラムに関する。
【背景技術】
【０００２】
樹体生産能力の指標の一つである樹体の葉面積を非破壊で推定する方法について、これまでに種々検討されており、その１つとして光学的樹体構造測定装置を用いて樹葉の透過光を計測することにより樹体総葉面積を推定する方法がある。これは、孤立した樹体ないし間隔をおいて多数存在する樹体下の位置に光学的樹体構造測定装置を配置して、樹葉の透過光を計測して得られたデータを用いて樹体の総葉面積を推定し、樹体生産力を評価するものである。【０００３】
光学的樹体構造測定装置としてプラントキャノピーアナライザー（ＰＣＡ：商品名）を用いて樹体の葉面積推定を行うことについて、次のような文献に開示されている。
【非特許文献１】中野敬之「プラントキャノピーアナライザーの茶樹への適用」（日本作物学会記事69(3):419-423(2000)）
【非特許文献２】J.S.Broadhead et al.「Comparison of method for leaf area in three rows」 Agricultural and Forest Meteorology 115:151-161(2003))
【０００４】
非特許文献１では、孤立した群落である茶樹の樹列の葉面積診断にプラントキャノピーアナライザーを適用している。測定時にセンサーを樹体中心部に向け、方位方向（水平方向）の視野を規制するため開度９０°のビューキャップを装着し、また天頂角方向の視野を０°から６０°に規制する視野規制用フィルターを取り付けている。この方法では、葉面積推定の際に一様な個体群の仮定のまま計算し、また、茶樹の場合に葉群に比して枝の存在量が多いために葉面積指数を過大評価することになり、葉面積推定には使用し難いものとなっている。【０００５】
非特許文献２では、ケニアのサバンナにおける樹木（Croton megalocarpus、Melia volkensil）の樹列の葉面積密度（単位体積当たりの葉の面積の比率）の推定にプラントキャノピーアナライザーを適用し、樹列の横断面を楕円形と仮定して推定した光路長を用い、孤立樹の測定方法に基づいて葉面積密度を算出している。この方法では、算出結果と実測した葉面積密度との相関が低くなり、そのため新たに個々の葉の傾斜角の天頂角分布を実測しそれに基づいたモデルにより算出結果と実測した値との相関を高めている。しかしながら、新たなモデルの追加のために個々の葉の傾斜角を測定するのに多大の労力と時間を要するため、実用的な方法とは言えないものであった。【０００６】
代表的に用いられている光学的樹体構造測定装置として、ＰＣＡ（プラントキャノピーアナライザー：米国ＬＩ−ＣＯＲ社製のＬＡＩ−２０００）があるが、これは、図１に示すようなもので、先端側の入射角１４８°の範囲の光を取り込む魚眼レンズを含むレンズ系、反射鏡、フィルター、検出器を含む構成であり、検出器は受光素子のシリコンが同心円状に配設されていて、５つの異なる天頂角の光を検出するようにしてあり、検出器での測定により得られたデータを処理し、コンピュータへのデータの転送制御を行うものである。
ＰＣＡでの測定対象は、一般的には個体群の高さｚの３倍以上の半径方向の広がりをもち、方位による葉面積の分布及び高さがほぼ一様な植物個体群であり、図２に示すような断面形状になっている。この場合、地表面で個体群の透過率を測定した際の天頂角θにおける光路長Ｓ（θ）は以下の式により推定される。
Ｓ（θ）＝ｚ／ｃｏｓθ【０００７】
しかしながら、上記式の仮定の成立しない孤立した群落（孤立した１本の樹、すなわち孤立樹等）の場合、一般的には幹の根元で群落外側へセンサーを向けた透過光測定とともに、測定地点から天頂角θにおける光路長を実測する必要がある。孤立群落の形状及び葉面積は方位により一様ではないため、４方位程度の測定を行うことが必要になり、結局、合計２０箇所の長さ測定（５箇所の光路長×４方位）を行う必要がある。さらに、センサーを外側に向けるため、個体群中心部上方は測定対象にならない。また、この方法で求められるのは葉面積密度（樹体単位体積当たりの葉面積）であるため、総葉面積を算出するためには群落の体積について別途測定する必要がある。【０００８】
本発明者は、カンキツ樹の３３樹体について、ＰＣＡを用いて樹体の透過光を計測し、前述のような方位による葉面積の分布及び高さがほぼ一様な植物個体群を想定して樹体の葉面積の推定を行った。その結果では、ＰＣＡによる葉面積指数（単位面積当たりの葉の面積の比率）と樹冠葉面積指数（樹冠を土地に投影した面積当たりの葉の面積の比率）との関係は図３に示すようになり、葉面積密度の高低に関係なく相関が低く（ｒ＝０．４１８）、分散している。このように、ＰＣＡによる葉面積指数と樹冠葉面積指数との相関が低くなる要因として、測定対象となるカンキツ樹の樹体が本来ＰＣＡの対象とする高さの均一な樹体群ではないこと等が考えられる。
【発明の開示】
【発明が解決しようとする課題】
【０００９】
光学的樹冠構造測定装置を用いたこれまでの樹冠葉面積の推定方法では、葉面積密度の低い旧葉（冬を越した葉）のみの樹体には適用できず、計測することが困難になることがあり、新葉（冬を過ぎて新たに展開した葉）と旧葉の比率の評価が行えないという問題があった。また、孤立した、あるいは間隔をおいて存在し高さが均一な樹体群とならない樹体に対して樹冠葉面積の推定を実施する場合に、実施に要する時間が多くなるとともに、測定結果として得られる葉面積指数と樹冠葉面積指数との相関が低いものとならざるを得なかった。そのため、実施に多くの時間を要することなく、測定結果として得られる葉面積指数と樹冠葉面積指数との相関を高いものとし、また、必要な経費を少なくすることが求められていた。
【課題を解決するための手段】
【００１０】
本発明は、前述した課題を解決すべくなしたものであり、本発明による樹体生産能力を評価する方法は、孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について樹体の透過光を測定することにより樹体生産能力を評価する方法であって、１本の樹体の主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体とする回転半長円体モデルを決定することと、該主幹から所定の距離だけ離れた地面近辺の位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して該光学的樹体構造測定装置により測定対象となる複数の天頂角で入射し樹葉を通過した光の強度を測定することと、天空からの光が遮られない位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して該光学的樹体構造測定装置により測定対象となる前記複数の天頂角での樹葉を通過しない光の強度を測定することと、前記複数の天頂角で入射し樹葉を通過した光の強度と樹葉を通過しない光の強度とから該複数の天頂角での空隙率を求めることと、測定対象とする前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積を回転半長円体モデルに基づいて計算することと、前記複数の天頂角での空隙率と前記複数の天頂角に対応する光路長とから各天頂角に対応する光の減衰量を求めることと、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求めることと、からなるものである。【００１１】
前記測定された複数の天頂角での入射光の強度から、次式（１）
【数１】
2007069736.tif 000003
により該複数の天頂角での空隙率Ｔ（θｉ）を求め、次式（２）
【数２】
2007069736.tif 000004
により葉面積密度（ＬＡＤ）を求めるようにしてもよい。【００１２】
また、孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について、樹体の透過光を測定することにより樹体生産能力を評価する方法であって、１本の樹体の主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体とする回転半長円体モデルを決定することと、該主幹から所定の距離だけ離れた地面近辺の位置に魚眼レンズを装着した撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置することと、前記撮像装置により測定対象となる複数の天頂角からの入射光を含む画角の画像を撮像することと、前記撮像装置により取得された画像データについて前記複数の天頂角での入射光の開空度を求める演算処理を行うことと、測定対象とする前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積を回転半長円体モデルに基づいて計算することと、前記複数の天頂角からの入射光の開空度と前記複数の天頂角に対応する光路長とから各天頂角に対応する光の減衰量を求めることと、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求めることと、からなるものとしてもよい。【００１３】
前記撮像装置により撮像された画像について、次式（５）
【数５】
2007069736.tif 000005
により前記複数の天頂角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000006
により葉面積密度（ＬＡＤ）を求めるようにしてもよい。
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしてもよい。【００１４】
本発明による樹体生産能力を評価するための撮像装置は、魚眼レンズと、該魚眼レンズの結像位置に配設された撮像素子と、撮像動作制御部と、前記撮像素子により取得された画像データについて演算処理を行う演算回路と、該演算回路による演算結果を出力する出力部とを備えてなる撮像装置であって、前記演算回路は、孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について樹体の透過光を測定する際に、１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に前記撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置して撮像し取得された画像について測定対象となる複数の天頂角をそれぞれ含む複数の角度範囲の画像部分におる画像データから前記複数の角度範囲での開空度を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記撮像装置まで通過する光路長及び樹冠体積を求め、前記開空度と前記複数の天頂角に対する光路長とから各天頂角を含む角度範囲に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うものである。【００１５】
前記演算回路は取得された画像データについて、次式（５）
【数５】
2007069736.tif 000007
により前記複数の天長角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000008
により葉面積密度（ＬＡＤ）を求める演算処理を行うようにしてもよい。
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしてもよい。【００１６】
本発明による樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラムは、孤立して存在する樹体または間隔をおいて複数存在する個々の樹体について該樹体の透過光を測定する際に、１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して測定された複数の天頂角で入射し樹葉を通過した光の強度と、天空からの光が遮られない位置に光学的樹体構造測定装置を入射光軸が水平面に対して垂直上方に向くようにして配置して測定された前記複数の天頂角での樹葉を通過しない光の強度とから該複数の天頂角での空隙率を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて求められた前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記光学的樹体構造測定装置までに通過する光路長及び樹冠体積とから各天頂角に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うものである。【００１７】
前記測定された複数の天頂角での入射光の強度から、次式（１）
【数１】
2007069736.tif 000009
により該複数の天頂角での空隙率Ｔ（θｉ）を求め、次式（２）
【数２】
2007069736.tif 000010
により葉面積密度（ＬＡＤ）を求めるようにしてもよい。【００１８】
また、孤立して存在する樹体または間隔をおいて複数存在する樹体について樹体の透過光を測定する際に１本の樹体の主幹から所定の距離だけ離れた地面近辺の位置に魚眼レンズを装着した撮像装置を入射光軸が水平面に対して垂直上方に向くようにして配置して撮像し取得された画像について測定対象となる複数の天頂角をそれぞれ含む複数の角度範囲の画像部分におる画像データから前記複数の角度範囲での開空度を求め、前記樹体についてその主幹を軸とし樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半長円体として決定された回転半長円体モデルに基づいて前記複数の天頂角に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から前記撮像装置まで通過する光路長及び樹冠体積を求め、前記開空度と前記複数の天頂角に対応する光路長とから各天頂角を含む角度範囲に対応する光の減衰量を求め、該光の減衰量から葉面積密度を求めこれと前記樹冠体積とから樹体総葉面積を求める演算処理を行うものとしてもよい。【００１９】
複数の天頂角からの入射光を含む画角について取得された画像データについて、次式（５）
【数５】
2007069736.tif 000011
により前記複数の天長角での入射光に対する開空度をＴ（θｉ）求め、次式（２）
【数２】
2007069736.tif 000012
により葉面積密度（ＬＡＤ）を求めるようにしたことを特徴とする請求項１１に記載の樹体生産能力を評価するための演算処理をコンピュータ上で実行するためのプログラム。
前記回転半長円体モデルを、主幹を軸として樹高及び平均樹冠半径をそれぞれ２つの径とする上に凸の回転半楕円体モデルとしてもよい。
【発明の効果】
【００２０】
本発明においては、孤立した、または間隔をおいた樹体に対して半長円体モデルを適用し、光学的樹体構造測定装置を用いて樹体の透過光を測定することにより、あるいは魚眼レンズを装着した撮像装置を用いて測定対象となる複数の天頂角からの入射光を含む画角の画像を撮像することにより、葉面積密度の低い樹体や旧葉のみの樹体に対しても樹体総葉面積が適切に評価でき、かつ多くの時間、労力を要することなく葉面積指数と樹冠葉面積指数との相関が高い、実測の場合と同等の測定結果が得られるという効果がある。また、魚眼レンズを装着した撮像装置を用いることにより、測定を実行するための経費を軽減できる。
【発明を実施するための最良の形態】
【００２１】
本発明において、高さの均一な樹体群とならない、孤立した、あるいは間隔をおいて存在する樹体について樹体総葉面積を推定するにあたり、樹体について回転半長円体モデルを設定し、それに基づき光学的樹体構造測定手段あるいは魚眼レンズを装着した撮像装置を用いて測定して得られたデータから葉面積指数を求めることにより樹体生産力の評価を行う。そこで、まず回転半長円体モデルを用いた樹体生産力の評価ついて説明する。【００２２】
〔回転半長円体モデルを用いた樹体生産力の評価〕
図４に回転半長円体モデルの例として回転半楕円体モデルを示す。このモデルは、樹体の断面形状を、樹高と平均樹冠半径とをそれぞれ長径及び短径とする上に凸の半楕円とし、この半楕円の長径または短径の一方を軸とする回転体が樹体を表すものとして考えるものである。図４において、回転半楕円体モデルによる樹体の高さをｈ、半径をａとする。ｈ，ａは個々の樹体に応じて決まる量である。Ｓは入射光が回転半楕円体モデルでの回転半楕円体面における入射点から地表の測定点まで通過する光路長であり、天頂角θ（地表面に対して垂直の方向と測定方向とのなす角）によって異なる。【００２３】
主幹表面から距離ｃの位置で地面近辺に光学的樹体構造測定装置を配置して測定対象となる複数の天頂角で入射し樹葉を通過した光の強度と、天空からの光が遮られない位置に光学的樹体構造測定装置を配置して複数の天頂角での樹葉を通過しない光の強度とを測定して複数の天頂角での空隙率を求め、あるいは、主幹表面から距離ｃの位置で地面近辺に魚眼レンズを装着した撮像装置を配置して測定対象となる複数の天頂角を含む画角の撮像を行い取得された画像データについて演算処理を行って測定対象となる複数の天頂角で入射した光についての開空度を求める。【００２４】
〔光学的樹体構造測定装置による測定〕
回転半楕円体モデルに基づき代表的な光学的樹体構造測定装置であるプラントキャノピーアナライザー（ＰＣＡ）を用いて実際に樹体生産力を評価し、測定の有効性について、カンキツ樹の樹体を対象として検証を行った例について説明する。カンキツ樹の場合、主幹の径は１０ｃｍ程度で、主幹表面からの距離ｃは数１０ｃｍ（２０〜１００ｃｍ）程度になる。ｃの値を設定し、主幹表面からの距離がｃで地面近辺にＰＣＡを入射光軸が鉛直上方を向くように配置して樹葉を透過した光の強度（群落下部の光強度）を計測した。また、天空からの光が遮られない位置にＰＣＡを入射光軸が鉛直上方を向くように配置して複数の天頂角での入射光の強度（群落上部の光強度）を測定した。天空からの光が遮られない位置としては、樹体と同等の高さ、またはそれ異常の高さとするのがよい。ＰＣＡでは図４において異なる５つの天頂角θ（θ１，θ２，・・・，θ５）の方向について、樹葉を透過した光の強度を計測するものであり、用いたＰＣＡでは、天頂角θは、それぞれ７°，２２°，３８°，５２°，６８°としてある。【００２５】
ＰＣＡでの測定は、太陽の直達光が卓越する条件においては過小評価になるので、太陽光を直接受けないように散乱光が卓越する曇天日、もしくは晴天日の場合に太陽高度が低く直達光の影響が少ない早朝あるいは夕方に行うことが良好な測定結果を得るために望ましい。【００２６】
回転半楕円体モデルの場合、各天頂角θ１，θ２，・・・，θ５に対して入射光が回転半楕円体モデルでの回転半楕円体面における入射点から測定点まで通過する光路長Ｓ（θ１）〜Ｓ（θ５）がそれぞれ求められる。ただし、ＰＣＡを配置した測定位置が回転半楕円体の軸となる主幹から外れていることにより、光路長は方位によって異なる。そのため、どの方位による光路長を規定しておく必要があり、ここでは、ＰＣＡの位置から主幹の付け根に向かう方位での光路長をＳ（θ１）〜Ｓ（θ５）とする。各光路についてＰＣＡで計測した樹葉を透過した光の強度（群落下部の光強度）及び樹葉を透過しない光の強度（群落上部の光強度）と光路長（Ｓ（θｉ））とから光量の減衰量が求められる。回転半楕円体モデルを用いて得られた葉面積密度（ＬＡＤ）に、同様に半楕円体モデルによる寸法データから算出した樹体体積（Ｖ）を乗ずることで樹体総葉面積（ＬＡ）が算出される。【００２７】
まず、式（１）
【数１】
2007069736.tif 000013
により天頂角θｉでの空隙率（透過率）Ｔ（θｉ）を求め、このように求められたＴ（θｉ）と式（２）
【数２】
2007069736.tif 000014
（ｗ（θｉ）は天頂角ごとの重み付け係数である。）により葉面積密度（ＬＡＤ）を求める。次に、回転半楕円体モデルにより樹冠体積（Ｖ）を式（３）
【数３】
2007069736.tif 000015
（Ｒｍは平均樹冠半径（ｍ）、Ｈは樹高（ｍ）である。）により求め、式（２）と式（３）とから、総葉面積（ＬＡ）を式（４）
【数４】
2007069736.tif 000016
により求める。【００２８】
本発明において、回転半楕円体モデルによるＰＣＡでの測定をカンキツ樹の１９樹体について行った。主幹の根元から距離ｃが６０ｃｍの地表位置に樹体中心部に向けてＰＣＡを配置して測定した。測定項目としては、ＰＣＡから樹体中心部に向かう方位として４方位の透過光の測定と、樹体の寸法（４方位の樹冠半径と樹高の５箇所）となる。樹体形状を、４方位の樹冠半径を平均した平均樹冠半径をもつ上に凸の回転半楕円体と仮定しているため、上記５箇所の寸法測定により、樹冠体積（Ｖ）と測定地点での各天頂角θでの光路長を推定可能である。【００２９】
一方、対比のために実測により、樹体の葉面積を求めた。これは、樹体の葉の枚数を計数し、サンプル的に選んだ数十枚の葉の面積を計測し、それから総葉面積を求めるものである。上記ＰＣＡを用いた測定により得られた樹体の葉面積密度（樹体単位体積当たりの葉面積：ＬＡＤ）と実測による葉面積密度とを対比して、図５に示す結果が得られた。この場合回帰直線の傾きはほぼ１となり、また、かなり高い相関（ｒ＝０．６１７）が得られた。【００３０】
回転半楕円体モデルに基づいてＰＣＡを用いて得られた樹体総葉面積と実測による樹体総葉面積との関係を図６に示す。同図において、きわめて樹体が大きい１点だけ大きく外れているが、全体として非常に高い相関（ｒ＝０．８９７）が得られている。【００３１】
なお、本発明の実施例としては、回転半長円体モデルとして、上に凸の半楕円を中心軸の周りに回転した形状を考えたが、これはモデルにより光路長を計算する上での手法であり、回転体の断面形状は単に二次曲線としての楕円に限らず、光路の計算に利用できる他の曲線としてもよく、回転半長円体は回転半楕円体、あるいは特殊な場合として半球体を含む、より一般的な回転体の形状として考えることができ、樹体の形態に応じて適宜使い分けられるものである。【００３２】
従来のＰＣＡによる樹冠葉面積指数の推定方法では、葉面積密度の低い樹体には適用できず、また、旧葉のみの樹体を計測すること、新旧葉比の評価が行うことが困難であった。これは旧葉のみの樹体が新葉展開により葉面積を増加していく過程においては、葉面積密度と樹体体積との両方が増加するということにもよるのであるが、その点からも本発明による半長円体モデルとＰＣＡとを用いた総葉面積推定の方法は樹体生産力の評価指標として非常に有効であると言える。【００３３】
また、ＰＣＡを用いた従来の総葉面積の推定方法では、本発明のように回転半楕円体モデルを用いていないので、２０箇所の光路長（５箇所の光路長×４方位）を測定する必要があり、葉面積密度から総葉面積を算出するために、群落の体積を別途測定する必要があって、測定に多くの時間、労力を要しており、さらに、ＰＣＡを主幹の根元に配置して外方に向けるため、樹体中心部上方が測定対象とならないなどの要因があって、良好な測定結果が得難いものであった。それに対し、本発明による総葉面積の測定方法では、これらの難点が解消され、従来の方法に比べて実施に要する時間が１／３以下になると言える。【００３４】
〔魚眼レンズを装着した撮像装置による測定〕
回転半長円体モデルにより、光学的測定手段として魚眼レンズを装着した撮像装置（デジタルカメラ）を用いた測定によっても、樹体の総葉面積を推定できる。本発明において、回転半長円体モデルとして図４に示すような上に凸の回転半楕円体モデルを考え、主幹の表面からの距離ｃが３０ｃｍの位置で値面近辺の位置に視野角１８４度の全周魚眼レンズを装着した撮像装置を、レンズを鉛直上方に向けて配置し、撮像することにより測定を行った。撮像の際の天候、時間等の条件はＰＣＡを用いた場合と同様である。【００３５】
撮像装置では撮像素子（ＣＣＤ）により二次元的な輝度分布を示す画像が得られるが、その画像をパーソナルコンピュータに取り込み、画像解析を行って開空度を算出することにより、ＰＣＡの場合と同様に総葉面積を求める。ＰＣＡのシステムにおいて同心上に配置された測光素子により５つの天長角θ１，θ２，・・・，θ５（７°，２２°，３８°，５２°，６８°））の方向について樹葉を透過した光の強度を計測するのであるが、撮像装置による画像を用いる場合、画像において分布する輝度を幹と葉群を含む樹体の部分とそれ以外の部分（主として天空）とに分けて２値化した後に、それぞれ５つの天長角θ１，θ２，・・・，θ５を中心とした設定した角度範囲（０〜１４°，１５〜２９°，３０〜４４°，４５〜５９°，６０〜７４°）について樹体部分（黒）と天空部分（白）とのそれぞれのピクセル数を計数し、式（５）
【数５】
2007069736.tif 000017
により開空度Ｔ（θｉ）を求める。開空度は、光学的樹体構造測定装置を用いた場合の空隙率と同等のものとして規定される量である。開空度の計算は、撮影全方位範囲のうち主幹部を向いた中心角９０°の扇形部分について行った。さらに、この開空度Ｔ（θｉ）と、回転半楕円体モデルにより求めた光路長Ｓ（θｉ）とから、式（２）
【数２】
2007069736.tif 000018
（ｗ（θｉ）は天頂角ごとの重み付け係数であり、Ｔ（θｉ）は天頂角に対応する開空度である。）により葉面積密度ＬＡＤを求め、回転半楕円体モデルにより樹冠体積（Ｖ）を式（３）
【数３】
2007069736.tif 000019
（Ｒｍは平均樹冠半径（ｍ）、Ｈは樹高（ｍ）である。）により求め、式（６）と式（３）とから、総葉面積（ＬＡ）を式（４）
【数４】
2007069736.tif 000020
により求める。葉面積密度の算出において、主幹部を画像内に多く含む６８°（６０〜７４°）の値を除いている。
このようにして、魚眼レンズを装着した撮像装置での撮像により取得された画像データからの演算処理により樹体総葉面積が求められる。一方で、樹体の葉の枚数を計数し、サンプル的に選んだ数十枚の葉の面積を計測し、それから総葉面積を求めることにより実測値を得ている。【００３６】
魚眼レンズを装着した撮像装置による画像を用いた測定により得られた樹体の葉面積密度（ＬＡＤ）と実測による葉面積密度とを対比して、図７に示す結果が得られた。この場合、ＰＣＡを用いた場合よりも相関が高くなっており（それぞれＲ２＝０．５５３，０．４７２）、過小評価傾向は魚眼レンズを装着した撮像装置を用いた場合により大きくなっている。これは、葉群の重なりが存在する場合の測定機器への影響の程度が異なるためであると考えられ、魚眼レンズを装着した撮像装置による画像では葉群が何枚重なっても画像に写らなければ評価できないのに対して、ＰＣＡでは透過光を測定するために葉群の重なりをある程度まで評価できることによるものと考えられる。【００３７】
回転半楕円体モデルに基づいて魚眼レンズを装着した撮像装置を用いた測定により得られた樹体総葉面積と実測による樹体総葉面積とを対比して、図８に示す結果が得られた。他の値から大きく外れた１点を除外して、魚眼レンズを装着した撮像装置の場合にも、ＰＣＡの場合と同様に高い相関が得られている（それぞれＲ２＝０．７６４，０．７９０）。
魚眼レンズを装着した撮像装置を用いる場合、測定に要する経費をＰＣＡの場合よりも少なくすることができる。また、ＰＣＡの場合には受光素子が各天頂角に対応する同心円状のものであるために、測定結果が受光素子上の１点で受ける強い光の影響を受けることがあるが、デジタルカメラでの撮像素子は多数の画素ごとに受光した結果を用いるために、このような強い光の影響は少ないと言える。さらに、魚眼レンズを装着した撮像装置の場合にも、回転半長円体モデルとして、縦方向の断面形状を楕円以外のものとしてもよいことは、ＰＣＡの場合と同様である。【００３８】
樹体生産能力を評価するための魚眼レンズを装着した撮像装置は、図９に示すように、撮像光学系の魚眼レンズ１、撮像素子（ＣＣＤ）２、撮像駆動回路３、ディスプレイ４を含む撮像装置としての構成に加えて、複数の天頂角に対する開空度を求め、葉面積密度、樹体総葉面積を求める演算処理を行うための演算回路５を含むものとする。演算回路５は撮像装置のボディ内にＩＣチップを内蔵した形態としてもよく、あるいは撮像装置のＩ／Ｏ端子６を介して接続されるユニットの形態としてもよい。さらに、このような内蔵型あるいは外部接続ユニット型の演算回路に対し、演算に必要なパラメータ（樹体の高さ、平均樹冠半径等）の入力手段を備えるようにする。【００３９】
このように、魚眼レンズを装着した撮像装置に演算回路を備える形態でない場合には、撮像装置により取得された画像を記憶したメモリーカードを用いて、パーソナルコンピュータにより樹体生産能力を評価するための演算を行うか、撮像装置のＩ／Ｏ端子から画像データをケーブルを介してパーソナルコンピュータに転送してパーソナルコンピュータにより樹体生産能力を評価するための演算を行う。本発明は、このように撮像装置により取得された画像データから回転半長円体モデルに基づいて葉面積密度、樹体総葉面積を求める演算処理を行うためのプログラムとしても特徴づけられるものである。
【図面の簡単な説明】
【００４０】
【図１】樹冠構造測定装置の概略的構成を示す図である。
【図２】測定対象となる一般的な樹体群の断面形状を示す図である。
【図３】ＰＣＡを用いた従来の方法により得られた葉面積指数と樹冠面積指数との関係を示すグラフである。
【図４】回転半楕円体モデルについて説明する図である。
【図５】ＰＣＡを用いて回転半楕円体モデルに基づき得られた樹体の葉面積密度と実測による葉面積密度とを対比して示すグラフである。
【図６】ＰＣＡを用いて回転半楕円体モデルを用いて得られた樹体総葉面積と実測による樹体総葉面積との関係を示す図である。
【図７】魚眼レンズを装着した撮像装置を用いて回転半楕円体モデルに基づき得られた樹体の葉面積密度と実測による葉面積密度とを対比して示すグラフである。
【図８】魚眼レンズを装着した撮像装置を用いて回転半楕円体モデルを用いて得られた樹体総葉面積と実測による樹体総葉面積との関係を示す図である。
【図９】魚眼レンズを装着した撮像装置の概略的構成を示す図である。
【符号の説明】
【００４１】
１  魚眼レンズ
２  撮像素子
３  撮像駆動回路
４  ディスプレイ
５  演算回路
【図１】
2007069736.tif 000031
【図２】
2007069736.tif 000032
【図３】
2007069736.tif 000033
【図４】
2007069736.tif 000034
【図５】
2007069736.tif 000035
【図６】
2007069736.tif 000036
【図７】
2007069736.tif 000037
【図８】
2007069736.tif 000038
【図９】
2007069736.tif 000039
【国際調査報告】
2007069736.tif 000040
2007069736.tif 000041