(57)【要約】【課題】  樹木を撮像して得られる画像を解析することにより、その画像に含まれる樹幹や葉や空を識別することが可能な画像解析システムを提供する。【解決手段】  樹木が撮影された全天画像Ｄ１を解析して対象を識別する画像解析システムＳであり、全天画像Ｄ１に含まれる略扇形状領域を識別する樹幹部識別手段ｍ２を備える。さらに、所定の方法で全天画像Ｄ１を二値化画像に変換する前処理手段ｍ１を備え、前記樹幹部識別手段ｍ２は、前処理手段ｍ１により生成された二値化画像Ｄ５について、中心から円周方向に黒ピクセルが所定の割合で含まれる第１ラインを探索し、当該第１ラインの中心側端部を起点として円周方向に黒ピクセルが所定の割合で含まれる二つの第２ラインを探索し、二つの第２ラインと円弧とに囲まれた扇形パターンＰを記憶し、隣接する第１ラインについてその扇形パターンＰを合成することが好ましい。【選択図】  図１
2005253407.tif 000002
【特許請求の範囲】
【請求項１】
樹木が撮影された全天画像を解析して対象を識別する画像解析システムであり、
全天画像に含まれる略扇形状領域を識別する樹幹部識別手段を備えることを特徴とする画像解析システム。
【請求項２】
所定の方法で全天画像を二値化画像に変換する前処理手段を備え、前記樹幹部識別手段は、前処理手段により生成された二値化画像について、中心から円周方向に黒ピクセルが所定の割合で含まれる第１ラインを探索し、当該第１ラインの中心側端部を起点として円周方向に黒ピクセルが所定の割合で含まれる二つの第２ラインを探索し、二つの第２ラインと円弧とに囲まれた扇形パターンを記憶し、隣接する第１ラインについてその扇形パターンを合成することを特徴とする請求項１記載の画像解析システム。
【請求項３】
前記各扇形パターンの第１ラインのライン長を微分し、微分値が０より大であり且つ第１ラインが互いに隣接する複数の扇形パターンを合成し、その合成範囲の幅を樹幹の幅とすることを特徴とする請求項２記載の画像解析システム。
【請求項４】
前記各扇形パターンの第１ラインのライン長を微分し、微分値が０より大であり且つ第１ラインが互いに隣接する複数の扇形パターンのなかから、最大微分値を有する扇形パターンの第１ラインのライン長を樹幹の高さとすることを特徴とする請求項２記載の画像解析システム。
【請求項５】
樹木が撮影された画像を解析して対象を識別する画像解析システムであり、
彩度が所定範囲であり、且つ、色成分のうち緑成分の輝度が最も高いピクセルの集合を識別する葉部識別手段を備えることを特徴とする画像解析システム。
【請求項６】
樹木が撮影された画像を解析して対象を識別する画像解析システムであり、
彩度が所定範囲であり、且つ、総ての色成分の平均輝度が所定範囲のピクセルの集合を識別する空部識別手段を備えることを特徴とする画像解析システム。
【請求項７】
樹木が撮影された画像を分析して対象を識別する画像解析システムであり、
各色成分について、輝度値とピクセル数とを軸とするヒストグラムを作成し、閾値と閾値との間のピクセル数が所定数となるように閾値を設定し、閾値と閾値の間隔が略均一となるように間隔を狭広させるとともに、その狭広に連動させてピクセルの位置を移動させるコントラスト調整手段を備えることを特徴とする画像解析システム。
【請求項８】
前記前処理手段が、前記コントラスト調整手段によりコントラスト調整を行い、且つ、前記葉部識別手段により識別された領域と、前記空部識別手段により識別された領域を白ピクセルとし、これらの領域を除く部分を黒ピクセルとして二値化画像に変換することを特徴とする請求項２記載の画像解析システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、樹木が撮影された画像を用いて、その画像から樹幹や葉や空を識別する画像解析システムに関する。
【背景技術】
【０００２】
従来から、林内環境調査などにおいては、林内の樹木の本数や葉の量などが評価要素として利用されている。たとえば、林内環境を評価するにあたって、樹木の本数は植生評価の重要なパラメータとなり、樹冠付近における葉の占める割合は、地表に対する日照状態や降雨遮断量や蒸発散量を評価する上で重要なパラメータとなる。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
ところで、近年はデジタル画像情報の収集効率が飛躍的に向上したことを踏まえ、各種林内調査等においても、デジタル画像の実用的な活用が期待できるようになってきた。とくに、一部の機種では魚眼レンズや曲面レンズ等による全天画像（広角画像）の撮影も可能であり、効率的な全天画像の画像解析処理が行えるようになった。
【０００４】
しかしながら、全天画像は林内環境調査等において一部使用されるに留まり、その利用例は非常に少ない。一方で、林内で撮影した全天画像には様々な情報が含まれており、デジタル画像処理技術を効果的に活用することにより、多様な情報を得ることができるものと考えられる。とくに、樹木の葉や樹幹、背景である空は、林内環境調査等において利用価値が高く、全天画像から樹木の樹幹や葉、背景である空を識別することにより、その情報を様々に利用することが可能となる。たとえば、林内環境調査においては、全天画像に含まれる樹幹を識別することにより植生状況を把握することができ、葉や空を識別することにより地表における日照状態や降雨遮断量などを評価することが可能となる。
【０００５】
そこで、本発明の目的は、樹木を撮像して得られる画像を解析することにより、その画像に含まれる樹幹や葉や空を識別することが可能な画像解析システムを提供することにある。
【課題を解決するための手段】
【０００６】
本発明の画像解析システムは、樹木が撮影された全天画像を解析して対象を識別する画像解析システムであり、全天画像に含まれる略扇形状領域を識別する樹幹部識別手段を備えることを特徴とする。全天画像においては、樹幹の形状が略扇形状となる一定のパターンを示す傾向にある。この発明によれば、この傾向を利用して略扇形状領域を識別することにより、この略扇形状領域を全天画像に含まれる樹幹として識別・利用することができる。
【０００７】
本発明の画像解析システムは、所定の方法で全天画像を二値化画像に変換する前処理手段を備え、前記樹幹部識別手段は、前処理手段により生成された二値化画像について、中心から円周方向に黒ピクセルが所定の割合で含まれる第１ラインを探索し、当該第１ラインの中心側端部を起点として円周方向に黒ピクセルが所定の割合で含まれる二つの第２ラインを探索し、二つの第２ラインと円弧とに囲まれた扇形パターンを記憶し、隣接する第１ラインについてその扇形パターンを合成することを特徴とする。前処理手段により例えばピクセルの彩度や輝度等を閾値とする方法で全天画像を二値化画像に変換し、その二値化画像について黒ピクセルを手がかりとしながら上記探索を行い、隣接する第１ラインについて扇形パターンを合成することにより、この合成範囲が前記扇形状領域とほぼ一致することとなり、この合成範囲を樹幹部分として識別することが可能となる。
【０００８】
本発明の画像解析システムは、前記各扇形パターンの第１ラインのライン長を微分し、微分値が０より大であり且つ第１ラインが互いに隣接する複数の扇形パターンを合成し、その合成範囲の幅を樹幹の幅とすることが好ましい。この発明によれば、全天画像において互いに一部が重なり合う樹幹が存在しても、微分によって樹幹の重なりを排除し、各々独立した樹幹としてその幅を求めることができる。
【０００９】
本発明の画像解析システムは、前記各扇形パターンの第１ラインのライン長を微分し、微分値が０より大であり且つ第１ラインが互いに隣接する複数の扇形パターンのなかから、最大微分値を有する扇形パターンの第１ラインのライン長を樹幹の高さとすることが好ましい。この発明によれば、ライン長を微分することにより、上述のように重なり合って撮影された樹幹を各々独立した樹幹として識別することができ、その中で最大微分値を有する第１ラインの長さを樹幹の高さとするため、たとえ樹木が一部重なり合って撮影されていたとしても、独立した樹幹としてその高さを求めることができる。
【００１０】
本発明の画像解析システムは、樹木が撮影された画像を解析して対象を識別する画像解析システムであり、彩度が所定範囲であり、且つ、色成分のうち緑成分の輝度が最も高いピクセルの集合を識別する葉部識別手段を備えることを特徴とする。樹木を撮影して得られる画像において、樹木の葉の部分は、ピクセルの彩度が一定範囲に集中する傾向にあり、さらに、色成分のうち緑成分の輝度が最も高い傾向にある。本発明によれば、この傾向を利用することにより、葉の部分を識別することができる。
【００１１】
本発明の画像解析システムは、樹木が撮影された画像を解析して対象を識別する画像解析システムであり、彩度が所定範囲であり、且つ、総ての色成分の平均輝度が所定範囲のピクセルの集合を識別する空部識別手段を備えることを特徴とする。樹木を撮影して得られる画像において、空の部分は、ピクセルの彩度が一定範囲に集中する傾向にあり、さらに、総ての色成分の平均輝度が一定範囲に集中する傾向にある。本発明によれば、この傾向を利用することにより、空の部分を識別することができる。
【００１２】
本発明の画像解析システムは、樹木が撮影された画像を分析して対象を識別する画像解析システムであり、各色成分について、輝度値とピクセル数とを軸とするヒストグラムを作成し、閾値と閾値との間のピクセル数が所定数となるように閾値を設定し、閾値と閾値の間隔が略均一となるように間隔を狭広させるとともに、その狭広に連動させてピクセルの位置を移動させるコントラスト調整手段を備えることを特徴とする。撮影された画像は、撮影条件やカメラの機種などの違いを原因として、明るさの分布に偏りがある。この発明によれば、ヒストグラムを分割する閾値の間隔が略均等となるように各間隔を狭広させることにより、各輝度値に対するピクセル数をほぼ均一化する。これにより、高すぎる輝度を抑え、低過ぎる輝度を上昇させて明るさの分布の偏りをなくし、コントラストを調整することができる。
【００１３】
本発明の画像解析システムは、前記前処理手段が、前記コントラスト調整手段によりコントラスト調整を行い、且つ、前記葉部識別手段により識別された領域と、前記空部識別手段により識別された領域を白ピクセルとし、これらの領域を除く部分を黒ピクセルとして二値化画像に変換することが好ましい。この発明によれば、コントラスト調整を行い、葉と空の部分を除いた二値化画像を処理するため、樹幹の識別精度を高めることができる。
【発明の効果】
【００１４】
本発明の画像解析システムによれば、全天画像に含まれる略扇形状領域を識別することにより、この略扇形状領域を樹幹として識別・利用することができる。この略扇形状領域は、例えば、全天画像を所定の方法で二値化した二値化画像について黒ピクセルを探索することにより精度良く識別することができる。黒ピクセルが所定割合以上含まれるラインのライン長を微分することにより、重なって撮影された複数の樹幹を各々独立した樹幹として、その幅や高さを求めることも可能となる。さらに、ピクセルの彩度や緑成分の輝度を解析することにより葉の部分を識別することが可能となり、ピクセルの彩度や平均輝度を解析することにより空の部分を識別することが可能となる。コントラスト調整することにより、さらに識別精度を高めることができる。本発明の画像解析システムで識別された樹幹や葉や空の情報は、林内環境調査などにおいて高い利用価値を有し、例えば、本システムにより識別された樹幹をカウントすることにより立木の本数を算出したり、葉や空の割合を地表における日照状態や降雨遮断量の評価などに利用したりすることが可能である。
【発明を実施するための最良の形態】
【００１５】
以下、本発明の画像解析システムＳについて、図面を参照しながら説明する。
【００１６】
図１は、本発明の画像解析システムＳのブロック図である。画像解析システムＳは、中央演算処理装置（ＣＰＵ）１と、メモリやハードディスクなどの記憶手段２とを備えるパーソナルコンピュータ等のコンピュータシステムである。画像解析システムＳは、デジタルカメラ等の撮影手段Ａ、及び、液晶ディスプレイなどの表示手段Ｂと接続可能となっている。本実施の形態では、撮影手段Ａと画像解析システムＳとは、インターネット等のネットワークを介して接続されている。画像解析システムＳには、画像解析プログラムが記憶されており、この画像解析プログラムが、後述の前処理手段ｍ１（コントラスト調整手段ｍ１１、葉部識別手段ｍ１２、空部識別手段ｍ１３、輪郭抽出手段ｍ１４）、樹幹部識別手段ｍ２（特徴抽出手段ｍ２１、識別処理手段２２）とを備える画像解析システムＳとしてコンピュータシステムを機能させている。
【００１７】
図２は、画像解析システムＳで使用される全天画像Ｄ１を示す図である。全天画像Ｄ１は、魚眼レンズを取り付けた撮影手段Ａを、調査対象の林分で光軸を林床の近似平面にほぼ直交するようにして、魚眼レンズを天方向に合わせ、所定の撮影高（例えば約１．２ｍ）で全天写真を撮影して得られる広角画像である。露出は、樹幹の色が明瞭に識別できることを目安に決定されている。全天画像Ｄ１は、外周が円形であり、中心から放射線状に樹木が並んでいる。全天画像Ｄ１に含まれる樹木は、中心側が先端であり、円周側が根元となっており、樹幹部分は略扇形状を呈する傾向を示す。撮影手段Ａは、レンズを天方向に向けて撮影されるため、背景の大部分は空となっている。全天画像Ｄ１は撮像手段ＡのメモリにJPEG形式で保存され、ネットワークを介して画像解析システムＳに転送される。
【００１８】
図３は、画像解析システムＳにより実現される画像解析方法を説明するフローチャートである。画像解析方法は、前処理ステップとしてコントラスト調整ステップとフィルタリングステップが行われ、その後に樹幹部識別ステップとして特徴抽出ステップと識別ステップとが順次行われる。
【００１９】
＜前処理ステップ＞
（コントラスト調整ステップ）
コントラスト調整ステップは、撮影条件やカメラの違いを原因としたコントラストのばらつきを解消するために行われるステップであり、コントラスト調整手段ｍ１１により実現される。図４は、コントラスト調整ステップの詳細を説明するフローチャートである。コントラスト調整手段ｍ１１は、全天画像Ｄ１の色要素ＲＧＢ各々について、輝度値をＸ軸、ピクセル数をＹ軸とする輝度の分布（ヒストグラム）を作成する（図４（ａ））。つぎに、ヒストグラム上のピクセル数を合計して総ピクセル数を算出し、閾値と閾値との間のピクセル数が所定数（望ましくは総ピクセル数／３２）よりも大であり且つその数に最も近い値となるように、ヒストグラムをＸ軸方向に分割する閾値を設定する（図４（ｂ））。閾値と閾値との間にプロットされているピクセル数はほぼ同一となる。このヒストグラムについて、Ｘ軸とＹ軸のメモリ（座標）を固定した状態で、閾値と閾値との間隔がほぼ均一になるように狭広させ、その狭広に連動させてプロットされているピクセルを移動させることによりコントラスト調整を行う（図４（ｃ））。以上の処理を色要素ＲＧＢそれぞれについて行う。図４（ｄ）にコントラスト調整済の全天画像Ｄ２を示す。この処理により、高すぎる輝度が抑えられ、低過ぎる輝度が上昇させられて、明るさの偏りが解消され、コントラストが調整された全天画像Ｄ２が得られる。
【００２０】
（フィルタリングステップ）
フィルタリングステップは、葉や空や下草等の樹幹幹以外の部分や、その他ノイズを除去し、後述の特徴抽出がスムーズに行えるようにフィルタリングを行うステップであり、画像解析システムＳの葉部識別手段ｍ１２と空部識別手段ｍ１３と輪郭抽出手段ｍ１４とにより実現される。図５は、フィルタリングステップを説明するフローチャートである。このフィルタリングステップにおいて、全天画像Ｄ２を所定の方法で二値化した全天画像Ｄ５に変換する。
【００２１】
二値化画像作成  
葉部識別手段ｍ１２は、彩度が所定範囲内であり、且つ、色成分ＲＧＢのうち緑成分Ｇの輝度が最も高いピクセルを白ピクセル（１）とする。葉部識別手段ｍ１２により白（１）とされたピクセルの集合が葉であると識別できる。また、空部識別手段ｍ１３は、彩度が所定範囲であり、且つ、総ての色成分の平均輝度が所定範囲のピクセルを白ピクセル（１）とする。空部識別手段ｍ１３により白（１）とされたピクセルの集合が空であると識別できる。そして、白ピクセル（１）以外のピクセルを黒ピクセル（０）とする。なお、識別精度の観点から、葉部識別手段ｍ１２と空部識別手段ｍ１３で使用される彩度の範囲は２５５階調のうち５〜３０が好ましく、さらには２０が好ましい。また、空部識別手段ｍ１３で使用される平均輝度の範囲は２３０以上２５５以下が好ましく、さらには２４０であることが好ましい。以上の処理により全天画像Ｄ２が二値化画像Ｄ３に変換される。この二値化された全天画像（二値化画像）Ｄ３は、樹幹が粗く抽出された状態となっている。なお、彩度とは各ピクセルのＲＧＢの中で輝度が最も高い色から輝度が最も低い色を引いた値である。
【００２２】
輪郭抽出  
つぎに、輪郭抽出手段ｍ１４は、コントラスト調整後の全天画像Ｄ２に対してＸ方向にハイパスフィルタ（ラプラシアンフィルタ）処理を行った画像と、Ｙ方向に同様の処理を行った画像とをピクセル単位で合成して全天画像Ｄ４を生成する。全天画像Ｄ４は、輪郭が抽出された状態となっている。輪郭抽出精度の観点から、フィルタ係数は｛−１／６, −１／３, １,−１／３, −１／６｝が好ましいが、これに限られるわけではない。なお、ハイパスフィルタとは、画像の明るさが急激に変化するエッジの検出や高周波成分の強調による画像の先鋭化に用いられるフィルタであり、ラプラシアンフィルタとは、一般的にエッジ検出や先鋭化に用いられる二次微分フィルタの一種であり、同様の機能を果たすものであれば他のフィルタを使用しても良い。
【００２３】
画像の合成  
さらに、二値化後の全天画像Ｄ３と、輪郭抽出後の全天画像Ｄ４とを各ピクセルで合成して全天画像Ｄ５を生成する。全天画像Ｄ５は、樹幹部分と地面が黒ピクセル（０）の状態となる。なお、フィルタリングステップでは、二値化後の全天画像Ｄ３のみを生成するか、又は、輪郭抽出後の全天画像Ｄ４のみを生成して、いずれか一方の全天画像Ｄ３，Ｄ４を後述の樹幹識別ステップで使用することも可能であるが、樹幹識別精度の観点からは、本実施の形態のように二つの全天画像Ｄ３,Ｄ４を合成した全天画像Ｄ５を生成して使用することが望ましい。
【００２４】
＜樹幹識別ステップ＞
樹幹部識別ステップは、全天画像に含まれる略扇形状領域を識別するステップである。全天画像においては、樹幹の形状が略扇形状となる一定のパターンを示す傾向にある。そこで、略扇形状領域を識別することにより、その識別した領域を樹幹の情報として利用することが可能となる。具体的には、図３に示すように、下記の特徴抽出ステップと識別ステップとを含み、各ステップは、樹幹抽出手段ｍ２の一部として機能する特徴抽出手段ｍ２１と識別処理手段ｍ２２により実現される。
【００２５】
（特徴抽出ステップ）
特徴抽出ステップは、特徴抽出手段ｍ２１により実現される。特徴抽出手段ｍ２１は、フィルタリング処理後の全天画像Ｄ５から、樹幹の特徴を持つ部分である略扇形状領域を抽出する機能を備える。図６は特徴抽出ステップを説明するフローチャートであり、図７はその処理の説明図である。
【００２６】
黒ピクセル探索１  
まず、特徴抽出手段ｍ２１は、全天画像Ｄ５について、中心から円周方向に向かって黒ピクセルを探索する（図７（ａ））。黒ピクセルが見つかった場合は、黒ピクセル探索２に進む。見つからない場合は、円周上の隣のピクセルに向かって繰り返し黒ピクセルを探索する（図７（ｂ））。なお、ここで円周上とは、全天画像の円周上でも良いし、全天画像の同心円上でも良い。本実施の形態では、全天画像よりも一回り小さい同心円の円周上までの探索を行っており、後述する位置ａ，ｂ，ｃ１，ｄ１，，，ｃｘ，ｄｘは全天画像よりも一回り小さい同心円の円周上に位置している。
【００２７】
黒ピクセル探索２  
黒ピクセルが見つかった場合、最初に黒ピクセルを発見した位置（起点）aから円周位置（終点）ｂまでの第１ラインＬ１上において黒ピクセルの割合が所定値以上であるかを確認する。すなわち、（第１ラインＬ１上の黒ピクセル数）／（ラインＬ１上の総ピクセル数）が所定値以上であるかを確認する（図７（ｃ）（ｄ）（ｅ））。黒ピクセルの割合は９０％以上とすると、ノイズの影響を受けずに精度の高い特徴抽出を行うことができるため好ましい。黒ピクセルの割合が所定値以上（９０％以上）である場合は、扇形パターン作成１に進む。所定値未満（９０％未満）の場合は、第１ラインＬ１上の次の黒ピクセルを探索する（図７（ｆ）。）なお、ここで「位置」は、ピクセルに対応する。
【００２８】
扇形パターン作成１
黒ピクセルが所定の割合（９０％以上）で含まれている第１ラインＬ１があった場合、円周上における位置（終点）ｂの両隣の位置（隣点）ｃ１,位置（隣点）ｄ１に対して、位置（起点）ａから位置（隣点）ｃ１の第２ラインＬｃ１上、及び、位置（起点）ａから位置（隣点）ｄ１までの第２ラインＬｄ１上に含まれる黒ピクセルの割合が所定値以上（好ましくは９０％以上）であるかどうかを確認する（図７（ｇ））。すなわち、第１ラインＬ１を挟むように二つの第２ラインＬｃ１、Ｌｄ１を設定して探索する。黒ピクセルの割合は、（第２ラインＬｃ１と第２ラインＬｄ１に含まれる黒ピクセル総数）／（第２ラインＬｃ１と第２ラインＬｄ１のピクセル総数）で算出する。黒ピクセルの割合が所定値以上（９０％以上）である場合は扇形パターン作成２に進む。黒ピクセル数の割合が所定値未満（９０未満）である場合は、ここで位置（起点）ａからの探索は終了し、位置（終点）ｂを円周上の一方方向にずらして、同様に黒ピクセル探索１，２を行う。
【００２９】
扇形パターン作成２  
位置（起点）ａから位置（隣点）ｃ１まで及び位置（隣点）ｄ１までのラインＬｃ１，Ｌｄ１上に黒ピクセルが所定の割合で含まれている場合、さらに円周上において位置（終点）ｂから離れる方向にある両隣（位置ｃ１の左隣の位置（隣点）ｃ２、位置（隣点）ｄ１の右隣の位置（隣点）ｄ２）に対して、同様に位置（起点）ａから位置（隣点）ｃ２,位置（隣点）ｄ２までの第２ラインＬｃ２，Ｌｄ２上に黒ピクセル数の割合が所定以上（９０％以上）であるかどうかを確認する（図７（ｈ））。所定以上（９０％以上）である場合、円周上の位置をｃ３，ｄ３，，，ｃｘ，ｄｘと左右に１ピクセルずつずらしながら第２ラインＬｃ３，Ｌｄ３，，，Ｌｃｘ，Ｌｄｘに対して同一処理を繰り返す。位置（隣点）ｃ(ｘ＋１)，ｄ(ｘ＋１)において黒ピクセルの割合が所定値未満（９０％未満）だった場合は、そこで位置（起点）ａからの探索は終了する。この繰り返しによって、位置（起点）ａを中心とし、第２ラインＬｃｘ，Ｌｄｘと、探索が完了した両隣の位置（隣点）ｃｘ，ｄｘを結ぶ円弧とにより扇形パターンＰが作成される（図７（ｉ））。ここで、位置（終点）ｂから、探索が完了した位置（隣点）ｃｘ、ｄｘまでの円周上にあるピクセルを、その扇形パターンＰの幅ピクセルと定義する。また、位置（終点）ｂから位置（隣点）ｃｘ又は位置（隣点）ｄｘまでの距離（ピクセルの数）を、その扇形パターンＰの幅ｗと定義する。
【００３０】
繰り返し処理  
位置（終点）ｂを円周に沿って１ピクセルずつずらしながら、位置（終点）ｂごとに、黒ピクセル探索１、黒ピクセル探索２、扇形パターン作成１、扇形パターン作成２を繰り返し実行する。扇形パターン作成１や扇形パターン作成２において黒ピクセルの割合が所定値未満（９０％未満）となった時点でその作成した扇形パターンＰを記憶手段に格納する。繰り返し処理により、位置（終点）ｂごとに扇形パターンＰが作成され、扇形パターンＰの集合としてデータが格納される。扇形パターンＰのデータの格納形式は、円周上のピクセル数分のメモリ領域に、それぞれ中心からのピクセル数、及び、幅ピクセル数を格納し、扇形パターン数値データとする（図６参照）。黒ピクセルの割合が９０％未満であった場合のデータは、「０，０」となる。
【００３１】
（識別ステップ）
識別ステップは、識別処理手段ｍ２２により実現される。識別処理手段ｍ２２は、全天画像において互いに一部が重なり合う樹幹が存在しても、微分によって樹幹の重なりを排除し、各々独立した樹幹として識別し、各樹幹の幅や樹幹の高さを求める機能を備える。図８は、識別ステップを説明するフローチャートであり、図９はその説明図である。
【００３２】
第１ラインのライン長算出  
扇形パターン数値データをもとに、各扇形パターンＰについて、位置（起点）ａから位置（終点）ｂまでの第１ラインＬ１のライン長（ピクセル数）ｈを計算する。ここで、扇形パターン数値データには、中心から位置（起点）ａまでの距離が格納されているため、ライン長ｈは（中心から円周までの半径距離）−（中心から位置（起点）ａまでの距離）により算出する（図９（ａ）（ｂ））。図９（ｃ）（ｄ）は、各扇形パターンＰについて求められたライン長ｈをグラフ化した説明図である。
【００３３】
第１ラインのライン長微分  
識別処理手段ｍ２２は、各扇形パターンＰのライン長ｈを微分し、扇形パターン微分値データを作成する。図９（ｅ）は扇形パターン微分値データをグラフ化した説明図である。
【００３４】
樹幹の幅・高さ識別  
扇形パターン微分値データにおいて、０より大の微分値を持つ扇形パターンが一つ、又は、０より大の微分値を有し且つ第１ラインが互いに隣接する複数の扇形パターンＰが出現している場合は、当該一つの扇形パターンＰ、又は、複数の扇形パターンＰの合成を一つの独立した樹幹とみなし、その合成範囲の幅Ｗを算出する。幅Ｗは、（０より大の微分値を有し且つ第１ラインが互いに隣接する複数の扇形パターンＰの数）＋（その複数の扇形パターンＰのなかで、両端に位置する扇形パターンＰ２，Ｐ３の幅ｗ）により求められる。なお、微分値が０より大の扇形パターンＰが一つの場合は、幅Ｗ＝（１＋幅ｗ＋（幅ｗ−１））により求められる。上記複数の扇形パターンＰの集合のうち、最大微分値を有する扇形パターンＰ１の第１ラインのライン長ｈがその樹幹の高さＨとなる。これによれば、全天画像Ｄ１に樹幹が重なり合って撮影されている場合でも、微分によりその重なりを排除して、独立した１本の樹幹とし、その高さＨや幅Ｗを求めることができる。具体的には、扇形パターン微分数値データを一方方向（図９中の識別方向）から探索し、微分値が０より大の扇形パターンＰあり、一つ前の扇形データＰが微分値０以下の扇形パターンＰ３から、０以下になる直前の扇形パターンＰ２を探し、それぞれの幅を加算して木の幹の幅とする（図９（ｆ）（ｇ））。
【００３５】
描画  
扇形パターン微分値データを探索して、連続して出現する微分値０より大の扇形パターンＰの集合のうち、最大の微分値を有する扇形パターンＰ１の位置（起点）ａを頂点として、扇形ＰＰを描写する（図９（ｈ）（ｉ）（ｊ））。図１の全天画像Ｄ６は、この描画の結果得られる全天画像である。この扇形ＰＰは樹幹部分である略扇形領域とほぼ一致するものであり、この扇形を樹幹として識別することができる。
【００３６】
以上、本実施の形態の画像解析システムＳを使用すれば、樹幹や葉や空を識別することができる。これらの情報は、森林環境調査などにおいて、立木の本数をカウントしたり、地表に対する日照状態の把握に利用することが可能である。
【００３７】
なお、本実施の形態では、全天画像を例に説明したが、葉部識別手段による葉の識別や、空部識別手段による空の識別に関しては、全天画像に限らず、パノラマ画像や一般的な画像などの様々な画像を対象として空や葉を識別することが可能である。
【図面の簡単な説明】
【００３８】
【図１】本発明の画像解析システムのブロック図。
【図２】画像解析システムで使用される全天画像を示す図。
【図３】画像解析システムにより実現される画像解析方法を説明するフローチャート。
【図４】コントラスト調整ステップの詳細を説明するフローチャート。
【図５】フィルタリングステップを説明するフローチャート。
【図６】特徴抽出ステップを説明するフローチャート。
【図７】特徴抽出ステップの処理を説明する説明図。
【図８】識別ステップを説明するフローチャート。
【図９】識別ステップを説明する説明図。
【符号の説明】
【００３９】
Ｓ      画像解析システム
Ａ      撮影手段
Ｂ      表示装置
ｍ１    前処理手段
ｍ１１  コントラスト調整手段
ｍ１２  葉部識別手段
ｍ１３  空部識別手段
ｍ１４  輪郭抽出手段
ｍ２    樹幹部識別手段
ｍ２１  特徴抽出手段
ｍ２２  識別処理手段
１      ＣＰＵ
２      記憶手段
Ｄ１，Ｄ２，Ｄ３，Ｄ４，Ｄ５，Ｄ６    全天画像
Ｐ      扇形パターン
ＰＰ    扇形
【図１】
2005253407.tif 000003
【図２】
2005253407.tif 000004
【図３】
2005253407.tif 000005
【図４】
2005253407.tif 000006
【図５】
2005253407.tif 000007
【図６】
2005253407.tif 000008
【図７】
2005253407.tif 000009
【図８】
2005253407.tif 000010
【図９】
2005253407.tif 000011