JP2015026369A - 林相解析装置、林相解析方法及びプログラム - Google Patents

Abstract

【課題】航空写真等から森林変化の検出を高精度で自動的に行うことを可能とする。【解決手段】林相区画生成部２０は、森林を含む同一の対象地域を互いに異なる時期に上空から撮影した２時期のオルソ画像データ３０から、各画素での画素値が当該画素に対応する地点についての２時期のオルソ画像データ３０の画素値の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて対象地域を複数領域に分割して森林の林相区画を生成する。特徴情報抽出部２２は２時期の各オルソ画像データ３０から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する。林相変化判別部２４は各林相区画について、２時期のオルソ画像データ３０の特徴情報に基づいて林相の変化の有無を判別する。【選択図】図１

Description

本発明は、上空から撮影した２時期の画像から森林の林相の変化を検出する林相解析装置、林相解析方法及びプログラムに関する。

航空機などによる上空からのリモートセンシングは、広範囲での地上の状況把握が可能であり、伐採や風倒被害等による森林変化の検出や、森林の林相区分図の作成に利用されている。ちなみに、我が国においては林相区分図の更新は数年程度ごとに行われている。

航空写真等の上空からの高分解能画像データを利用した森林変化の検出や林相区分図の更新は従来、基本的に人間が目視判読によって行っている。しかしながら、広範囲にわたる目視判読は膨大な時間を要するため、より効率的な森林変化の検出、ならびに、林相区分図更新の手法が求められている。

そこで、森林変化の自動判読についてこれまで数多くの研究が行われてきた。当該研究の多くは、まず、変化前後の２時期の画像それぞれから林相区分図を自動的に作成し、それらを比較して森林変化を検出しようとするものである。

また、従来、衛星等のリモートセンシング画像のスペクトル情報を利用し、２時期の画像の差分や比演算等によって森林変化を検出する方法が知られている。

特開２００８−０４６８３７号公報 特開２０１１−０２４４７１号公報

例えば、林相区分図の自動生成に関し、森林を撮影した画像を領域分割して、類似した特徴を持つ複数の小領域である林相区画を生成し、当該林相区画ごとに画像分類を行って林相区分・樹種分類図を作成する手法が研究されている。ここで、画像領域分割手法を利用した林相区画の生成は、画像のスペクトル情報をベースとして行われる。しかしながら、森林の画像における異なる林分間の特徴の違いはそれほど顕著でなく、一方、同じ林分内でも日向領域と日陰領域の違いにより画像の特徴のばらつきが大きい。また、使用する画像データの撮影状況（撮影時期、バンド数など）も領域分割処理の結果に影響を与える。そのため、画像領域分割により林相区画を適切に生成することは容易ではない。

このような事情から現状の自動生成される林相区分図における誤分類は比較的高く、自動生成した２時期の林相区分図を比較して変化検出を行う手法は誤検出の割合が高くなるという問題を有する。

また、２時期の画像のスペクトル情報を利用して変化検出を自動的に行う手法に関しては、森林の季節的な変動に応じて画像のスペクトル情報が大きく変わり得る。しかし、航空写真等の高分解能画像を利用する場合は、飛行機の調達や天候等の影響により、撮影を行う２時期にて季節変動に関する森林の状態を揃えることができるとは限らない。そのため、伐採等による森林変化の正確な検出が難しくなることがあった。

本発明は上記問題点を解決するためになされたものであり、航空写真等の高分解能の画像に基づいて森林変化の検出を高精度で自動的に行うことを可能とし、ひいては林相区分図の自動更新を可能とする林相解析装置、林相解析方法及びプログラムを提供することを目的とする。

（１）本発明に係る林相解析装置は、森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成手段と、前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出手段と、前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別手段と、を有する。

（２）上記（１）に記載する林相解析装置において、さらに、前記第１の原画像の撮影時における既知の林相区分図に基づいて前記第２の原画像の撮影時における目的とする林相区分図を生成する林相区分図生成手段を備え、前記林相区分図生成手段が、前記林相の変化が検出された前記林相区画について、前記第２の原画像の林相を前記テクスチャ情報とスペクトル情報との少なくとも一方に基づいて判別し、当該林相により前記既知の林相区分図を修正して前記目的の林相区分図を生成するように構成することができる。

（３）上記（１）又は（２）に記載する林相解析装置において、さらに、前記複合画像の隣接する複数画素からなる領域を画素として分解能を低下させた劣化複合画像を生成する画像劣化手段を備え、前記林相区画生成手段が、前記劣化複合画像の１又は複数種類の前記区画用特徴量の類似性に基づいて、前記対象地域を複数領域に分割し林相区画を生成する初期区画手段と、隣接する林相区画を１又は複数種類の前記区画用特徴量に基づいて結合して新たな林相区画を生成する逐次区画処理を１回又は複数回行い、前記対象地域における前記森林の林相区画を出力する後続区画手段と、を備え、前記後続区画手段が少なくとも１回の前記逐次区画処理を前記複合画像の画素における前記区画用特徴量、又は前記複合画像及び前記劣化複合画像それぞれの前記区画用特徴量に基づいて行うように構成することができる。

（４）本発明に係る林相解析方法は、森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成ステップと、前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出ステップと、前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別ステップと、を有する。

（５）本発明に係るプログラムは、コンピュータに、森林を含む対象地域における林相解析を行わせるためのプログラムであって、当該コンピュータを、森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成手段、前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出手段、及び、前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別手段、として機能させる。

本発明によれば、航空写真等の高分解能の画像に基づいて森林変化の検出を高精度で自動的に行うことが可能となり、ひいては林相区分図の自動更新が可能となる。

本発明の第１の実施形態である林相解析システムの概略の構成を示すブロック図である。 日向／日陰分布パターンの例を示す説明図である。 日向／日陰境界パターンの例を示す説明図である。 ＬＢＰ画像の例を示す説明図である。 正規化処理前のオルソ画像の例を示す画像である。 図５の画像を正規化処理して得られた画像である。 本発明の実施形態である林相解析システムにおける林相区画生成の概略の処理フロー図である。 階層的な領域分割処理の最初の階層にて低分解能画像から生成された林相区画の例を示す画像である。 階層的な領域分割処理の途中階層にて生成された林相区画の例を示す画像である。 図９に示す階層より後の階層にて生成された林相区画の例を示す画像である。 本発明の実施形態である林相解析システムにおける林相変化箇所の検出及び林相区分図の生成の概略の処理フロー図である。 本発明の第２の実施形態である林相解析システムの概略の構成を示すブロック図である。

以下、本発明の実施の形態（以下実施形態という）について、図面に基づいて説明する。

［第１の実施形態］
図１は、第１の実施形態である林相解析システム２の概略の構成を示すブロック図である。本システムは、演算処理装置４、記憶装置６、入力装置８及び出力装置１０を含んで構成される。演算処理装置４として、本システムの処理を行う専用のハードウェアを作ることも可能であるが、本実施形態では演算処理装置４は、コンピュータ及び、当該コンピュータ上で実行されるプログラムを用いて構築される。

演算処理装置４は、コンピュータのＣＰＵ（Central Processing Unit）からなり、後述する林相区画生成部２０、特徴情報抽出部２２、林相変化判別部２４、林相区分図生成部２６として機能する。

記憶装置６はＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）等のメモリ装置である。記憶装置６は演算処理装置４にて実行される各種のプログラムや、本システムの処理に必要な各種データを記憶し、演算処理装置４との間でこれらの情報を入出力する。例えば、記憶装置６には、オルソ画像データ３０、林相区画に関するスペクトル条件３２及び幾何条件３４、林相を特徴付けるテクスチャ情報３６及びスペクトル情報３８が予め格納される。また、記憶装置６には、既知の林相区分図データ５０を予め格納しておくことができる。

オルソ画像データ３０は航空機等から撮影された空中写真画像に基づいて生成される。本システムではオルソ画像データ３０は、森林を含む対象地域を上空から撮影した原画像として用いられる高分解能画像データであり、例えば、赤（Ｒ）、緑（Ｇ）、青（Ｂ）の３成分からなる３バンドのマルチスペクトル画像、又はこれらに近赤外（ＮＩＲ）を加えた４成分からなる４バンドのマルチスペクトル画像である。本システムは森林の変化を検出することができ、当該変化を検出する２時期Ｔ１，Ｔ２それぞれにおけるオルソ画像データ３０が記憶装置６に格納される。

スペクトル条件３２及び幾何条件３４は、対象地域を領域分割して生成される林相区画が画像にて満たすべき条件である。スペクトル条件３２は林相区画を構成する画素群の色（Color）や濃度・明るさから抽出されるスペクトル特徴量に関する条件を規定する。一方、幾何条件３４は林相区画の幾何特徴量である形状（Shape）に関する条件を規定する。例えば、形状は、コンパクトネス（Compactness）やスムースネス（Smoothness）といったパラメータで表される。ここで、スペクトル条件３２及び幾何条件３４は、各時期の原画像に関してではなく、それらをまとめて１画像とした複合画像に関する条件として規定される。この点については後述する。

テクスチャ情報３６及びスペクトル情報３８は、林相それぞれを特徴付ける情報である。テクスチャ情報３６は高分解能画像における画像テクスチャを表現する情報である。本実施形態では、画像テクスチャとして日向／日陰分布パターン、日向／日陰境界パターン及びＬＢＰ（Local Binary Pattern）画像を用いることに対応して、テクスチャ情報３６として日向／日陰分布パターン情報４０、日向／日陰境界パターン情報４２及びＬＢＰ（Local Binary Pattern）画像情報４４が記憶装置６に格納される。例えば、各テクスチャ情報３６は画像テクスチャの平均値や標準偏差などの統計量で表現することができる。

スペクトル情報３８は高分解能画像におけるスペクトルを表現する情報であり、画素値から抽出されるスペクトル特徴量で表現することができる。例えば、各スペクトル特徴量として各色の平均値、標準偏差などの統計量を用いることができる。

入力装置８は、キーボード、マウスなどであり、ユーザが本システムへの操作を行うために用いる。

出力装置１０は、ディスプレイ、プリンタなどであり、林相区画や林相区画での林相変化の有無や林相区分図（林分図）を画面表示、印刷等によりユーザに示す等に用いられる。また、林相区画や林相区分に関するデータを他のシステムで利用できるよう、データとして出力してもよい。

林相区画生成部２０は、２時期のオルソ画像データ３０に基づく複合画像に対して階層的に領域分割処理を行って林相区画を生成する。林相区画生成部２０は、複合画像生成部６０、低分解能画像生成部６２、初期区画部６４及び後続区画部６６を有する。

複合画像生成部６０は、同一の対象地域を撮影した２時期Ｔ１，Ｔ２の高分解能の原画像（オルソ画像データ３０）から複合画像を定義する。複合画像は、その各画素での画像情報が当該画素に対応する地点についての２時期の原画像の画像情報の組からなる画像であり、原画像と同様の高分解能の画像である。具体的には、複合画像の各画素での画像情報、つまり画素値は、当該画素と同一地点に対応する原画像の画素値の組で与えられる。例えば、原画像がＲ，Ｇ，Ｂからなる３バンドのマルチスペクトル画像である場合には、複合画像は時期Ｔ１の原画像の３バンドと時期Ｔ２の原画像の３バンドとを合わせた６バンドのマルチスペクトル画像を構成する。

なお、本実施形態では幾何補正後のオルソ画像を用いた場合を例に説明するが、オルソ画像の代わりに上空から撮影した画像を用いる場合には、時期Ｔ１の原画像と時期Ｔ２の原画像との間にて画素の配置にずれがある場合がある。この場合には事前に幾何補正を行い、画素配置のずれを吸収し、両原画像の画像情報を対応付けて複合画像の画素値を定義する。

低分解能画像生成部６２は、高分解能の複合画像をディグラデーション処理して分解能を低下させた低分解能画像（劣化複合画像）を生成する画像劣化手段であり、低分解能画像の画素は複合画像における隣接する複数画素からなる領域に対応する。ここでは区別を容易にするために、高分解能画像である複合画像（及び原画像）の画素を原画素、低分解能画像の画素を拡大画素と表現する。例えば、低分解能画像は複合画像を一様なサンプリング密度でダウンサンプリングして生成することができ、拡大画素の画素値は複合画像にて対応する領域内の複数の原画素を代表する値に設定することができる。例えば、低分解能画像生成部６２は、複合画像に対し所定のフィルタを用いた畳み込みにより平滑化処理を行い、拡大画素の画素値を当該拡大画素内の所定位置の原画素での平滑化処理後の画素値に設定することができる。一般的な画像では、注目画素に距離が近い画素の輝度値は注目画素の輝度値と近くなる場合が多く、注目画素から遠ざかるほど、注目画素の輝度値との差が大きくなる場合が多い。このような性質を考慮し、注目画素に近いほど平均値を計算するときの重みを大きくし、遠くなるほど重みを小さくするように２次元ガウス分布を用いて画像の平滑化を行うガウシアンフィルタが知られている。拡大画素の画素値を求める際の平滑化処理に、このガウシアンフィルタを採用することができる。また、より簡単な手法として例えば、拡大画素の画素値はそれに対応する複数の原画素の画素値の平均値、中央値などで定義することもできる。

初期区画部６４は階層的な領域分割処理の最初の階層の処理を行う。具体的には複合画像のスペクトル特徴量及び幾何特徴量を区画用特徴量とし、初期区画部６４は低分解能画像生成部６２により生成された低分解能画像を区画用特徴量の類似性に基づいて複数領域に分割し１次の林相区画を生成する。初期区画部６４は、互いに隣接する複数の拡大画素からなる領域を１次の林相区画とするか否かを、それら拡大画素の画素値のスペクトル特徴量の類似性（スペクトル条件）と、それら拡大画素を結合して得られる領域についての幾何特徴量の類似性（幾何条件）とに基づいて決定する。

後続区画部６６は階層的な領域分割処理の第２の階層以降の処理を行う。具体的には後続区画部６６は、隣接する林相区画を結合して新たな林相区画を生成する逐次区画処理を１回又は複数回行う。逐次区画処理は、互いに隣接する複数の低次の林相区画について、それらのスペクトル特徴量の類似性と、それらを結合して得られる高次の林相区画についての幾何特徴量の類似性とに基づいて、それら低次の林相区画を結合するか否かを決める。

本実施形態では後続区画部６６は逐次区画処理として、初期区画部６４にて生成された低分解能画像（劣化複合画像）の領域分割処理結果である１次の林相区画内、あるいは先行する領域分割処理結果である低次の林相区画内における低分解能画像と高分解能画像（複合画像）との両方、あるいは高分解能画像のみを用いて領域分割処理を行う。低分解能画像からは、林相区画の結合可否の判断に用いる画素値のスペクトル特徴量として、林相区画内に包含される拡大画素群についてのスペクトル特徴量を求める。一方、高分解能画像からは、林相区画の結合可否の判断に用いる画素値のスペクトル特徴量として、林相区画内に包含される原画素群についてのスペクトル特徴量を求める。なお、後続区画部６６は少なくとも１回の逐次区画処理が、高分解能画像の画素値のスペクトル特徴量、又は高分解能画像及び低分解能画像それぞれの画素値のスペクトル特徴量に基づいて行う領域分割処理である構成とすることもできる。すなわち、後続区画部６６が複数回の逐次区画処理を行う場合に、そのうち一部の回数は低分解能画像の画素値だけからスペクトル特徴量を抽出し、それに基づいて領域分割処理を行ってもよい。

画像領域分割には、例えば領域併合（region merging）に基づく手法を利用することができる。この手法において画像領域分割処理の対象とされる画像に存在する隣接する２つのオブジェクトを統合するかどうかは、統合後に生成される新しいオブジェクトの異質性（heterogeneity）と統合前のオブジェクトの異質性との間の変化を評価することによって決定される。ちなみに、初期区画部６４での処理ではオブジェクトは低分解能画像の画素であり、後続区画部６６での処理ではオブジェクトは、既に生成されている林相区画である。

具体的には、領域併合によるオブジェクトの異質性の変化Δｈは、併合前後におけるスペクトル異質性の変化Δｈ_ｐと形状異質性の変化Δｈ_ｔとから次式によって算出される。

ここで、ｗ_ｐはスペクトル異質性の重み、ｗ_ｔは形状異質性の重みである。

併合前後のスペクトル異質性の変化Δｈ_ｐは、対象画像の各バンドにおける併合前後のオブジェクト内の画素値の標準偏差を用いて、次式によって計算される。

ここで、Ｎは画像のバンド数、ｗ_ｉはバンドｉの重み、ｎ_ａｂは併合後の新しいオブジェクトの画素数、ｎ_ａ，ｎ_ｂは併合前の２つのオブジェクトの画素数、σ_ｉ，ａｂは併合後のオブジェクトのバンドｉにおける標準偏差、σ_ｉ，ａ，σ_ｉ，ｂは併合前の２つのオブジェクトのバンドｉにおける標準偏差である。

また、併合前後の形状異質性の変化Δｈ_ｔは、コンパクトネスとスムースネスという２つの基準で次式により定義される。

ここで、Δｈ_ｃは併合前後のコンパクトネスの変化、Δｈ_ｓは併合前後のスムースネスの変化、ｗ_ｃはコンパクトネスの重み、ｗ_ｓはスムースネスの重みである。

オブジェクトのコンパクトネス基準はオブジェクトの周囲長と面積から計算され、一方、スムースネス基準はオブジェクトの周囲長と境界ボックスの直径（長軸）から計算される。具体的にはΔｈ_ｃ，Δｈ_ｓは次式で定義される。

ここで、ｌ_ａｂは併合後のオブジェクトの周囲長、ｌ_ａ，ｌ_ｂは併合前の２つのオブジェクトの周囲長、ｓ_ａｂは併合後のオブジェクトの面積、ｓ_ａ，ｓ_ｂは併合前の２つのオブジェクトの面積、ｂ_ａｂは併合後のオブジェクトの境界ボックスの直径、ｂ_ａ，ｂ_ｂは併合前の２つのオブジェクトの境界ボックスの直径、ｎ_ａｂは併合後の新しいオブジェクトの画素数、ｎ_ａ，ｎ_ｂは併合前の２つのオブジェクトの画素数である。

併合前後のオブジェクトの異質性の変化Δｈが設定されたしきい値を超えない場合、領域の併合処理が実施され、しきい値を上回る場合、領域の併合処理を停止する。設定されたしきい値はスケールパラメータ（scale parameter）と呼ばれ、画像の分割処理によって生成されるオブジェクトの大きさを表している。スケールパラメータが大きいほど、より多くのオブジェクトが併合され、領域分割によって最終的に生成されるオブジェクトのサイズが大きくなる。

（１）式に示すように、初期区画部６４における拡大画素の結合の判断、又は後続区画部６６における低次の林相区画の結合の判断に対するスペクトル特徴量及び幾何特徴量それぞれの寄与比率は重みｗ_ｐ，ｗ_ｃにより調節することができる。ここで、初期区画部６４及び後続区画部６６の一方又は両方は、林相区画の生成判断において、幾何条件を用いずスペクトル特徴量の類似性だけに基づいて行う構成にすることもできる。

特徴情報抽出部２２は各原画像から特徴情報としてテクスチャ情報又はスペクトル情報あるいはこれら双方を抽出する特徴情報抽出手段であり、テクスチャ解析部７０及びスペクトル解析部７２を有する。

テクスチャ解析部７０は、オルソ画像データ３０からテクスチャ情報を抽出するテクスチャ情報抽出手段である。テクスチャ解析部７０は、日向／日陰分布パターン情報抽出部８０、日向／日陰境界パターン情報抽出部８２及びＬＢＰ画像情報抽出部８４を有する。

日向／日陰分布パターン情報抽出部８０はオルソ画像データ３０を二値化処理して日向と日陰とが区別された日向／日陰分布パターンをテクスチャ情報として抽出する。図２は日向／日陰分布パターンの例を示す説明図である。図２（ａ）〜（ｄ）はそれぞれスギ林、ヒノキ林、広葉樹林、非森林の例であり、左右に並ぶ２つの画像のうち左側がオルソ画像データ３０、右側が日向／日陰分布パターンである。当該日向／日陰分布パターンにおいて白領域が日向、黒領域が日陰である。ちなみに、ここでの二値化のしきい値は大津の手法により決定している。

日向／日陰境界パターン情報抽出部８２はオルソ画像データ３０における輝度勾配を二値化処理して日向と日陰との境界領域が現れた日向／日陰境界パターンをテクスチャ情報として抽出する。図３は日向／日陰境界パターンの例を示す説明図である。図３（ａ）〜（ｄ）は図２と同様、それぞれスギ林、ヒノキ林、広葉樹林、非森林の例であり、左右に並ぶ２つの画像のうち左側がオルソ画像データ３０、右側が日向／日陰境界パターンである。当該日向／日陰境界パターンにおいて白領域が輝度勾配がしきい値以上の領域であり、黒領域がしきい値未満の領域である。ここでも、しきい値は大津の手法により決定している。日向／日陰境界パターンからは日向と日陰とが切り替わる空間的な頻度の多寡が読み取れる。

ＬＢＰ画像情報抽出部８４はオルソ画像データ３０から局所二値パターン演算子を用いて得られるＬＢＰ画像をテクスチャ情報として抽出する。図４はＬＢＰ画像の例を示す説明図である。図４（ａ）〜（ｄ）は図２及び図３と同様、それぞれスギ林、ヒノキ林、広葉樹林、非森林の例であり、左右に並ぶ２つの画像のうち左側がオルソ画像データ３０、右側がＬＢＰ画像である。ＬＢＰ画像は原画像の詳細な模様構造パターンを反映しており、しかも画像のコントラストの影響を受けにくいという特性を有する。

テクスチャ解析部７０は抽出された日向／日陰分布パターン、日向／日陰境界パターン及びＬＢＰ画像の特徴量として、それぞれの平均値や標準偏差などの統計量を算出する。ちなみに算出される統計量の種類は、記憶装置６に日向／日陰分布パターン情報４０、日向／日陰境界パターン情報４２及びＬＢＰ画像情報４４として格納される統計量の種類に対応している。

スペクトル解析部７２はオルソ画像データ３０からスペクトル情報を抽出するスペクトル情報抽出手段である。本実施形態ではスペクトル解析部７２はＲ，Ｇ，Ｂ，ＮＩＲの４バンドのマルチスペクトル画像に対して次式で示される正規化処理を行い、正規化後の成分Ｒ’，Ｇ’，Ｂ’からなる画像からスペクトル特徴量を抽出する。またスペクトル解析部７２は日向／日陰分布パターン情報に基づいて各林相区画における日向領域のみからスペクトル特徴量を抽出する。

図５及び図６はオルソ画像データ３０に対する正規化処理の例を示す画像であり、図５は原画像、すなわち正規化処理前のオルソ画像、図６は図５に対応するオルソ画像データ３０を正規化処理して得られた画像である。図５及び図６はグレースケールで表示しているが本来、図５はＲＧＢ合成のカラー画像であり、図６はＲ’Ｇ’Ｂ’合成のカラー画像である。例えば、樹種等の違いに応じた近赤外成分の大きさの相違により、原画像と正規化処理画像との間には色合いなどに違いが生じる。図５と図６とにおける原画像と正規化処理画像との相違は、画像内での濃淡の変化が両画像間にて相違していることから読み取れる。

なお、スペクトル解析部７２は正規化処理を行わずにスペクトル特徴量を抽出する構成や、日陰を含んだ各林相区画の全体におけるスペクトル特徴量を抽出する構成とすることもできる。

林相変化判別部２４は、各林相区画について、時期Ｔ１，Ｔ２における原画像の特徴情報に基づいて林相の変化の有無を判別する。ここで、特徴情報は特徴情報抽出部２２により抽出されるテクスチャ特徴量又はスペクトル特徴量あるいはこれら双方である。例えば、各林相区画における変化有無は、時期Ｔ１，Ｔ２における原画像の特徴情報を対比することで検出することができる。具体的には、時期Ｔ１，Ｔ２間での特徴量の差分や比を求め、得られた結果をしきい値と比較して変化の有無を判別する。なお、スペクトル特徴量はバンド別の複数成分からなる。このスペクトル特徴量を変化検出に用いる場合には、それら各バンドの特徴量の２時期での差分値の平方和に応じた指標を算出し、これをしきい値と比較することができる。また、変化の有無として統計的な手法を用いることもできる。例えば、２時期の特徴量それぞれを変量とする多変量空間での林相区画の分布を求め、当該分布における位置に基づいて変化の有無を判別することができる。例えば、基本的には、林相が変化していない区画の分布はスギ林、ヒノキ林、広葉樹林、非森林といった林相の種類ごとに、テクスチャ情報３６やスペクトル情報３８として記憶装置６に格納される平均値や標準偏差に応じた位置、拡がりのピークを形成し、それらピークからはずれた位置となる林相区画は２時期間で林相が変化していると判断することができる。当該判断は具体的にはテクスチャ情報３６及びスペクトル情報３８に格納される各林相の平均値、標準偏差を指標として行うことができる。

林相区分図生成部２６は、２時期Ｔ１，Ｔ２のうち一方時期における既知の林相区分図に基づいて他方時期における林相区分図を生成する。例えば、時期Ｔ１を時期Ｔ２に先行する時期とし、時期Ｔ１における林相区分図が林相区分図データ５０として与えられている場合に、林相区分図生成部２６は時期Ｔ２における林相区分図を生成する更新処理を行う。具体的には、林相区分図生成部２６は、林相変化判別部２４により時期Ｔ１，Ｔ２間にて林相の変化が検出された林相区画について、時期Ｔ２の原画像での林相をテクスチャ情報又はスペクトル情報あるいはこれら双方に基づいて判別する。そして、当該林相により時期Ｔ１における既知の林相区分図を修正して時期Ｔ２における目的の林相区分図を生成する。

以下、林相解析システム２の処理の流れを説明する。図７は林相解析システム２における林相区画生成の概略の処理フロー図である。

林相解析システム２は複合画像生成部６０により、時期Ｔ１に撮影されたオルソ画像データ３０と時期Ｔ２に撮影されたオルソ画像データ３０に基づいて複合画像を生成する（Ｓ５）。なお、ここでは説明の便宜上、オルソ画像データ３０とは別途に複合画像を定義しているが、複合画像は本質的には２時期のオルソ画像データ３０の組であり、例えば記憶装置６上に１枚の画像として記憶されるような形で複合画像を生成することは必ずしも行わなくてもよい。

林相解析システム２は低分解能画像生成部６２により複合画像から低分解能画像（劣化複合画像）を生成し（Ｓ１０）、次に、階層的な画像領域分割処理Ｓ２０を行って、内部が一様な林相からなる領域に対応した林相区画を生成する（Ｓ３０）。領域分割処理Ｓ２０は初期区画部６４及び後続区画部６６により行われる。領域分割処理は低分解能画像から開始され、低分解能画像を構成する拡大画素の画素値を用いた林相区画生成処理（以下、区画処理とする）が最初に少なくとも１回行われ、その後、複合画像を構成する原画素及び低分解能画像を構成する拡大画素の画素値、あるいは原画素のみの画素値を用いた区画処理が少なくとも１回行われる。具体的には、初期区画部６４は低分解能画像を領域分割してスケールが小さい林相区画を生成する（Ｓ２２）。後続区画部６６は低分解能画像と高分解能画像とを併用し、林相区画同士を結合することによりスケールが大きくなった林相区画を生成する（Ｓ２４）。対象地域から最終的に抽出された林相区画のデータ９０は林相変化判別部２４にて利用される。

領域分割処理はスペクトル条件３２及び幾何条件３４を満たすように行われる。その際、スペクトル条件３２と幾何条件３４との比重、つまりそれぞれを領域分割に寄与させる度合は調節することができる。

ここで、領域分割処理の階層ごとに、林相区画を生成する際の条件は異なり得る。すなわち、初回の拡大画素を結合して初期の林相区画を生成する際や、低次の林相区画同士を結合して高次の林相区画を生成する際に、スケールパラメータであるΔｈは生成される林相区画のスケールが徐々に大きくなることを可能とするように設定される。Δｈの増加により、スペクトル条件３２に関しては、結合対象となる複数の拡大画素又は林相区画についてのスペクトル特徴量の類似性の判断基準（Δｈ_ｐ）が領域分割処理が高階になるにつれ緩和される。また、幾何条件３４に関しても幾何特徴量の類似性の判断基準（Δｈ_ｔ）が緩和され、領域分割処理が高階になるにつれより複雑な形状あるいは大きな面積を有する林相区画が許容されるようになり結合が促進される。一方、スペクトル条件３２と幾何条件３４とで総合的に結合が緩和されればよいので、スペクトル条件３２により結合を緩和しつつ、幾何条件３４は結合を或る程度抑制する方向に変化させより単純な形状の林相区画の生成を促す条件に設定することもできる。例えば、そのような調整は（１）式の重みｗ_ｐ，ｗ_ｔを変えることで可能であり、具体的にはｗ_ｐを低下させ、一方、ｗ_ｔを増加させて、Δｈ_ｐがΔｈに寄与しにくくし、Δｈ_ｔがΔｈに寄与しやすくすることで実現できる。

さて、既に述べたように、高分解能画像では森林の日向と日陰の違いにより同じ林相内でも画素値のばらつきが大きいので、従来、隣接している林相を区別して正確に区画するのが困難であるという問題があった。これに対し、林相解析システム２は原画像に基づく複合画像を低分解能化し、階層的な領域分割処理のうち少なくとも最初の１回を低分解能画像にて行う。画像を低分解能化すると同じ林相内での画素値のばらつきが平滑化・緩和されるので、低分解能画像に基づく区画処理は日向と日陰とでの画素値の違いの影響を受けにくくなり、隣接する異なる林相を区別し易い。なお、低分解能での区画処理を階層的な複数回の区画処理のうちの早い段階で行う方が、当該効果は高くなることが期待できる。

一方、階層的な区画処理の中で、低分解能画像と高分解能画像との両方、あるいは高分解能画像のみを用いて区画処理を行うことで、最終的に得られる林相区画に高分解能画像の情報を反映させることができる。

図８〜図１０は林相区画生成部２０の処理の例を示す画像であり、２時期のうちの一方にて撮影した対象地域の同一の空中写真画像に林相区画の境界を黒線で描いている。図８は初期区画部６４による処理結果であり、低分解能画像を利用して生成された初期の林相区画を示している。図９及び図１０は後続区画部６６による処理結果であり、図１０の画像は図９の画像より高階の領域分割処理の結果を示している。すなわち、図８、図９、図１０の順に階層的に領域分割が行われる。ここで、上述したように林相区画の生成は低分解能画像と高分解能画像とを併用して行われる。具体的には、図８〜図１０に示す３階層が林相区画を生成する領域分割処理の全階層である場合には、図８の林相区画は低分解能画像を利用して生成され、図９及び図１０の林相区画は、先行する領域分割の処理結果である林相区画内における低分解能画像と高分解能画像との両方、あるいは高分解能画像のみを利用して生成される。そして、図１０の林相区画が最終的に抽出された林相区画となる。

本システムでは、２時期に撮影された画像の両方に基づいて領域分割処理を行うので、これにより自動生成される林相区画では、２時期の間に森林変化が生じている領域と変化していない領域とが別々の区画となる。ちなみに、この段階では各林相区画における森林変化の有無は不明であり、後続の処理にて各林相区画における森林変化の有無が判別される。

ここで、各時期の画像から林相区分図を作成し、それらを比較して林相分類が変化した領域を抽出する従来手法では、各時期における林相の誤分類が２時期間での森林変化の検出精度を低下させる。これに対し、本システムでは２時期の画像それぞれを領域分割し林相分類するのではなく、まず２時期の画像を総合して領域分割を行い林相区画を生成する。この領域分割処理は従来手法で生じる精度低下を回避できる。つまり、領域分割処理では２時期間で森林が変化した領域と変化していない領域とが好適な精度で別々の区画に分けられる。そして、当該林相区画をベースにした後続の処理での森林変化の検出精度が向上する。

図１１は林相解析システム２における林相変化箇所の検出及び林相区分図の生成の概略の処理フロー図である。

林相解析システム２は林相区画生成部２０により生成される各林相区画での林相変化の有無の判別などに用いる特徴情報を特徴情報抽出部２２により２時期のオルソ画像データ３０それぞれから抽出する。具体的には、特徴情報抽出部２２は特徴情報として、テクスチャ解析部７０によりオルソ画像データ３０からテクスチャ情報を抽出する（Ｓ４０，Ｓ５０，Ｓ６０）。テクスチャ解析部７０は日向／日陰分布パターン情報抽出部８０により日向／日陰分布パターンを抽出し（Ｓ４０）、日向／日陰境界パターン情報抽出部８２により日向／日陰境界パターンを抽出し（Ｓ５０）、またＬＢＰ画像情報抽出部８４により、原画像における詳細な模様構造パターンが現れるＬＢＰ画像を抽出し（Ｓ６０）、テクスチャ情報を求める。

２時期のオルソ画像データ３０が季節変動における森林の状態が同じである場合には、スペクトル情報は林相変化の有無の判別に用いることができる。また、スペクトル情報は林相区分図生成部２６による林相区分図の更新処理にて利用され得る。そこで、林相解析システム２は特徴情報として、２時期のオルソ画像データ３０それぞれからスペクトル情報も抽出してもよい。具体的には、特徴情報抽出部２２はスペクトル解析部７２によりオルソ画像データ３０を正規化処理し、正規化後の画像からスペクトル特徴情報を抽出する（Ｓ７０）。さらにスペクトル解析部７２は日向／日陰分布パターン情報抽出部８０により抽出した日向／日陰分布パターンを用いて、森林のスペクトルの特徴が良好に現れ得る日向部分のみからスペクトル特徴情報を抽出する（Ｓ８０）。

林相解析システム２は林相変化判別部２４により、林相区画データ９０に規定される各林相区画について、２時期のオルソ画像データ３０それぞれから抽出した特徴情報に基づいて林相の変化の有無を判別し、林相区画単位で対象地域における林相の変化箇所を抽出する（Ｓ９０）。

林相区分図生成部２６は、記憶装置６から先行時期Ｔ１における林相区分図データ５０を読み出し、林相区分図データ５０が表す既知の林相区分図のうち林相変化が検出された林相区画における林相を時期Ｔ２における林相に更新して、時期Ｔ２における林相区分図のデータ１００を生成する（Ｓ１００）。

林相区分図生成部２６は、林相の変化が検知された林相区画について時期Ｔ２のオルソ画像データ３０から日向／日陰分布パターン情報抽出部８０、日向／日陰境界パターン情報抽出部８２及びＬＢＰ画像情報抽出部８４により抽出された日向／日陰分布パターン、日向／日陰境界パターン及びＬＢＰ画像がそれぞれ、記憶装置６にテクスチャ情報３６として格納されている日向／日陰分布パターン情報４０、日向／日陰境界パターン情報４２及びＬＢＰ画像情報４４それぞれのどの林相に対応する特徴を有しているかを判定する。また、林相区分図生成部２６は林相の変化が検知された林相区画について時期Ｔ２のオルソ画像データ３０からスペクトル解析部７２により抽出されたスペクトル特徴量が記憶装置６に格納されているスペクトル情報３８のどの林相に対応する特徴を有しているかを判定する。そして林相区分図生成部２６はこれら３種類のテクスチャ情報についての林相の判定結果とスペクトル情報についての林相の判定結果とから例えば、最近隣法などを用いて総合的に林相区画の林相を決定する。

なお、本実施形態では、３種類のテクスチャ情報とスペクトル情報とでまず別々に林相を判定しているが、それらの特徴の組み合わせに林相を対応付けた情報を記憶装置６に予め格納し、オルソ画像データ３０から抽出された特徴の組み合わせに基づいて林相を決定してもよい。また、３種類のテクスチャ情報のうちの一部だけを用いる構成としたり、テクスチャ情報とスペクトル情報との一方だけを用いる構成としたりすることもできる。

また、テクスチャ情報は上述の３種類以外のものを用いることもできる。例えば、ウェーブレット解析により抽出した画像テクスチャを用いることができる。具体的には、ウェーブレット解析で抽出した高周波成分の情報をテクスチャ情報として画像のスペクトル特徴情報と併用して林相を判定することができる。

テクスチャ解析部７０が抽出する日向／日陰分布パターン、日向／日陰境界パターン及びＬＢＰ画像は高分解能画像に現れる画像テクスチャ情報である。林相解析システム２はこれを林相変化検出の特徴情報として利用する。当該特徴情報は従来の画像テクスチャを表現する統計情報に比べ、林相の違いをより好適に表現することができるので、当該特徴情報を用いることで林相変化検出をより正確に行うことができる。また、画像テクスチャ情報を利用することで、オルソ画像データ３０における季節変動の影響が軽減されるので、画像スペクトル情報を利用する従来手法と比較して、森林伐採や風倒被害などの森林変化をより安定かつ正確に検出することが可能となり、林相区分図の更新や森林の精密管理をより効率的に行うことができる。

本実施形態ではオルソ画像データ３０は航空写真としたが、本発明は、例えば衛星写真などの他の高分解能画像を用いた林相解析にも適用可能である。また、本実施形態ではオルソ画像を例に説明したが、必ずしもオルソ画像を用いる必要は無く、上空から撮影した画像を適宜用いることができる。

［第２の実施形態］
第２の実施形態の林相解析システムにおいて上記第１の実施形態と同じ構成要素は同一の符号を付して説明を省略する。以下、第２の実施形態の林相解析システム２００について第１の実施形態の林相解析システム２との相違点を説明する。

図１２は第２の実施形態である林相解析システム２００の概略の構成を示すブロック図である。林相解析システム２００と林相解析システム２との相違点は林相区画をどのような条件に基づいて行うかという点にある。第１の実施形態では林相区画生成部２０はスペクトル特徴量についてのスペクトル条件３２及び幾何特徴量についての幾何条件３４を満たすように領域分割処理を行って林相区画を生成している。これに対し、第２の実施形態の林相解析システム２００では、領域併合前後のオブジェクト、つまり、劣化画像の画素又は林相区画の色、形状、テクスチャ及び面積を区画用特徴量とし、これらについての色条件２０２、形状条件２０４、テクスチャ条件２０６及び面積条件２０８を満たすように領域分割処理を行って林相区画を生成する。

色条件２０２、形状条件２０４、テクスチャ条件２０６及び面積条件２０８は、オルソ画像（オルソ画像データ３０）を領域分割して生成されるオブジェクト領域である林相区画が満たすべき条件である。色条件２０２は林相区画におけるオルソ画像の画素値に関する条件を規定する。ちなみに「色条件」における「色」はオルソ画像の画素値を意味し、ＲＧＢの情報からなるスペクトル情報よりも広義で用いている。具体的には、ここでの「色」は第１の実施形態のスペクトル特徴量に相当し、例えば、ＮＩＲの情報を含むマルチスペクトル画像の画素値もここでの色の概念に含まれ得る。形状条件２０４は林相区画の形状に関する条件を規定する。テクスチャ条件２０６は林相区画の画像テクスチャに関する条件を規定する。画像テクスチャはテクスチャ特徴量で評価される。面積条件２０８は林相区画の面積に関する条件を規定する。なお、林相区画の形状、面積は林相区画の幾何特徴量である。

初期区画部６４は、互いに隣接する複数の画素からなる領域を１次の林相区画とするか否かを、それら画素の画素値の類似性、領域の画像テクスチャの類似性、及び領域の形状や面積などの幾何特徴量の類似性に基づいて決定する。ここで、結合される領域間や結合前後の領域間での画素値の類似性は、画素値の平均値などの統計量などの特徴量に基づいて判断することができる。

後続区画部６６は第１の実施形態と同様、逐次区画処理を行う。本実施形態では逐次区画処理は、互いに隣接する複数の低次の林相区画について、それらにおける画素値の類似性に関する色条件２０２と、それらを結合して得られる高次の林相区画についての形状条件２０４と、テクスチャ条件２０６と、面積条件２０８とに基づいて、それら低次の林相区画を結合するか否かを決める。

画像領域分割には第１の実施形態と同様、領域併合に基づく手法を利用する。本実施形態では、領域併合によるオブジェクトの異質性の変化Δｈは、併合前後における色の異質性の変化Δｈ_colorと、形状の異質性の変化Δｈ_shapeと、テクスチャの異質性の変化Δｈ_textureと、面積の異質性の変化Δｈ_aeraとから次式によって算出される。

ここで、ｗ_colorは色の異質性の重み、ｗ_shapeは形状の異質性の重み、ｗ_textureはテクスチャの異質性の重み、ｗ_areaは面積の異質性の重みである。

併合前後の色の異質性の変化Δｈ_color、形状の異質性の変化Δｈ_shape、テクスチャの異質性の変化Δｈ_texture、面積の異質性の変化Δｈ_aeraは、例えば、それぞれ次式によって計算される。

ここで、ｎ_ａｂは併合後の新しいオブジェクトの画素数、ｎ_ａ，ｎ_ｂは併合前の２つのオブジェクトの画素数である。また、ＣＩ_ａｂは併合後のオブジェクトの色情報指数（例えば、画素値）、ＣＩ_ａ，ＣＩ_ｂは併合前の２つのオブジェクトの色情報指数、ＳＩ_ａｂは併合後のオブジェクトの形状情報指数（例えば、スムースネス）、ＳＩ_ａ，ＳＩ_ｂは併合前の２つのオブジェクトの形状情報指数、ＴＩ_ａｂは併合後のオブジェクトのテクスチャ情報指数（例えば、画素値の標準偏差値）、ＴＩ_ａ，ＴＩ_ｂは併合前の２つのオブジェクトのテクスチャ情報指数、ＡＩ_ａｂは併合後のオブジェクトの面積情報指数（例えば、面積値）、ＡＩ_ａ，ＡＩ_ｂは併合前の２つのオブジェクトの面積情報指数である。

（６）式に示すように、初期区画部６４における画素の結合の判断、又は後続区画部６６における低次の林相区画の結合の判断に対する色条件、形状条件、テクスチャ条件及び面積条件それぞれの寄与比率は重みｗ_color，ｗ_shape，ｗ_texture，ｗ_areaにより調節することができる。ここで、初期区画部６４及び後続区画部６６の一方又は両方は、林相区画の生成判断において、色条件、形状条件、テクスチャ条件及び面積条件の何れか１つもしくは複数の組み合わせに基づいて行う構成にすることもできる。

２，２００ 林相解析システム、４ 演算処理装置、６ 記憶装置、８ 入力装置、１０ 出力装置、２０ 林相区画生成部、２２ 特徴情報抽出部、２４ 林相変化判別部、２６ 林相区分図生成部、３０ オルソ画像データ、３２ スペクトル条件、３４ 幾何条件、３６ テクスチャ情報、３８ スペクトル情報、４０ 日向／日陰分布パターン情報、４２ 日向／日陰境界パターン情報、４４ ＬＢＰ画像情報、５０ 林相区分図データ、６０ 複合画像生成部、６２ 低分解能画像生成部、６４ 初期区画部、６６ 後続区画部、７０ テクスチャ解析部、７２ スペクトル解析部、８０ 日向／日陰分布パターン情報抽出部、８２ 日向／日陰境界パターン情報抽出部、８４ ＬＢＰ画像情報抽出部、２０２ 色条件、２０４ 形状条件、２０６ テクスチャ条件、２０８ 面積条件。

Claims (5)

1. 森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成手段と、
   前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出手段と、
   前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別手段と、
   を有することを特徴とする林相解析装置。
2. 請求項１に記載の林相解析装置において、
   前記第１の原画像の撮影時における既知の林相区分図に基づいて前記第２の原画像の撮影時における目的とする林相区分図を生成する林相区分図生成手段を有し、
   前記林相区分図生成手段は、前記林相の変化が検出された前記林相区画について、前記第２の原画像の林相を前記テクスチャ情報と前記スペクトル情報との少なくとも一方に基づいて判別し、当該林相により前記既知の林相区分図を修正して前記目的の林相区分図を生成すること、
   を特徴とする林相解析装置。
3. 請求項１又は請求項２に記載の林相解析装置において、
   前記複合画像の隣接する複数画素からなる領域を画素として分解能を低下させた劣化複合画像を生成する画像劣化手段を有し、
   前記林相区画生成手段は、
   前記劣化複合画像の１又は複数種類の前記区画用特徴量の類似性に基づいて、前記対象地域を複数領域に分割し林相区画を生成する初期区画手段と、
   隣接する林相区画を１又は複数種類の前記区画用特徴量に基づいて結合して新たな林相区画を生成する逐次区画処理を１回又は複数回行い、前記対象地域における前記森林の林相区画を出力する後続区画手段と、
   を有し、
   前記後続区画手段は少なくとも１回の前記逐次区画処理を前記複合画像の画素における前記区画用特徴量、又は前記複合画像及び前記劣化複合画像それぞれの前記区画用特徴量に基づいて行うこと、
   を特徴とする林相解析装置。
4. 森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成ステップと、
   前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出ステップと、
   前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別ステップと、
   を有することを特徴とする林相解析方法。
5. コンピュータに、森林を含む対象地域における林相解析を行わせるためのプログラムであって、当該コンピュータを、
   森林を含む同一の対象地域を互いに異なる時期に上空から撮影した第１の原画像及び第２の原画像から、各画素での画像情報が当該画素に対応する地点についての前記両原画像の画像情報の組からなる複合画像を定義し、当該複合画像の１又は複数種類の区画用特徴量の類似性に基づいて前記対象地域を複数領域に分割して前記森林の林相区画を生成し出力する林相区画生成手段、
   前記各原画像から特徴情報としてテクスチャ情報とスペクトル情報との少なくとも一方を抽出する特徴情報抽出手段、及び、
   前記各林相区画について、前記両原画像の前記特徴情報に基づいて林相の変化の有無を判別する林相変化判別手段、
   として機能させることを特徴とするプログラム。

Images