"計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成装置であって、前記航空レーザ測量データから生成された前記格子のXY座標値と地盤標高値からなるDEMデータと、前記XY座標値と表層標高値からなるDSMデータと、前記XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力部と、前記記憶手段に保存された前記DEMデータと前記DSMデータから、同一XY座標値に関する前記地盤標高値と前記表層標高値を読み出し、前記表層標高値から前記地盤標高値を差し引いた樹高値を、前記XY座標値とともにDHMデータとして前記記憶手段に保存するDHMデータ作成部と、前記DHMデータの任意の格子を標本格子と設定し、前記標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、前記最大傾斜角度を平均化した地上開度値を、前記標本格子の前記XY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算部と、前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度用カラーテーブルを前記反射強度値で検索し、前記反射強度値に対して設定された第1の表示色成分を取得して、前記XY座標値とともに反射強度画像データとして前記記憶手段に保存する反射強度段彩部と、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高用カラーテーブルを前記樹高値で検索し、前記樹高値に対して設定された第2の表示色成分を取得して、前記XY座標値とともに樹高画像データとして前記記憶手段に保存する樹高段彩部と、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度用カラーテーブルを前記地上開度値で検索し、前記地上開度値に対して設定された第3の表示色成分を取得して、前記XY座標値とともに樹冠立体画像データとして前記記憶手段に保存する地上開度段彩部と、前記記憶手段に保存された各画像データからそれぞれ画像を生成し、重ね合わせて植生図として端末に表示する地図合成部とを含むことを特徴とする植生図作成装置。前記反射強度用カラーテーブル、前記樹高用カラーテーブルおよび前記地上開度用カラーテーブルを、前記端末からの入力により変更する合成画像調整部をさらに含む請求項1記載の植生図作成装置。前記合成画像調整部は、ユーザが表示したい画像の種類を前記端末から設定した画像選択表を記憶手段に保存し、前記地図合成部は、前記画像選択表を参照し、選択された画像のみを前記端末に表示する請求項2記載の植生図作成装置。前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度値変換テーブルを前記反射強度値で検索し、前記反射強度値に対して設定されたa成分値に変換し、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高値変換テーブルを前記樹高値で検索し、前記樹高値に対して設定されたb成分値に変換し、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度値変換テーブルを前記地上開度値で検索し、前記地上開度値に対して設定されたL成分値に変換し、前記a成分値、前記b成分値、前記L成分値を前記XY座標値とともにLabカラー植生画像データとして前記記憶手段に保存するLabカラー段彩部をさらに含む請求項1乃至3のいずれか1項に記載の植生図作成装置。計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成装置であって、前記航空レーザ測量データから生成された前記格子のXY座標値と地盤標高値からなるDEMデータと、前記XY座標値と表層標高値からなるDSMデータと、前記XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力部と、前記記憶手段に保存された前記DEMデータと前記DSMデータから、同一XY座標値に関する前記地盤標高値と前記表層標高値を読み出し、前記表層標高値から前記地盤標高値を差し引いた樹高値を、前記XY座標値とともにDHMデータとして前記記憶手段に保存するDHMデータ作成部と、前記DHMデータの任意の格子を標本格子と設定し、前記標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、前記最大傾斜角度を平均化した地上開度値を、前記標本格子の前記XY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算部と、前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度値変換テーブルを前記反射強度値で検索し、前記反射強度値に対して設定されたa成分値に変換し、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高値変換テーブルを前記樹高値で検索し、前記樹高値に対して設定されたb成分値に変換し、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度値変換テーブルを前記地上開度値で検索し、前記地上開度値に対して設定されたL成分値に変換し、前記a成分値、前記b成分値、前記L成分値を前記XY座標値とともにLabカラー植生画像データとして前記記憶手段に保存するLabカラー段彩部と、前記Labカラー植生画像データから画像を生成し、植生図として端末に表示する地図合成部とを含むことを特徴とする植生図作成装置。計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成プログラムであって、コンピュータに、前記航空レーザ測量データから生成された前記格子のXY座標値と地盤標高値からなるDEMデータと、前記XY座標値と表層標高値からなるDSMデータと、前記XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力手段、前記記憶手段に保存された前記DEMデータと前記DSMデータから、同一XY座標値に関する前記地盤標高値と前記表層標高値を読み出し、前記表層標高値から前記地盤標高値を差し引いた樹高値を、前記XY座標値とともにDHMデータとして前記記憶手段に保存するDHMデータ作成手段、前記DHMデータの任意の格子を標本格子と設定し、前記標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、前記最大傾斜角度を平均化した地上開度値を、前記標本格子の前記XY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算手段、前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度用カラーテーブルを前記反射強度値で検索し、前記反射強度値に対して設定された第1の表示色成分を取得して、前記XY座標値とともに反射強度画像データとして前記記憶手段に保存する反射強度段彩手段、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高用カラーテーブルを前記樹高値で検索し、前記樹高値に対して設定された第2の表示色成分を取得して、前記XY座標値とともに樹高画像データとして前記記憶手段に保存する樹高段彩手段、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度用カラーテーブルを前記地上開度値で検索し、前記地上開度値に対して設定された第3の表示色成分を取得して、前記XY座標値とともに樹冠立体画像データとして前記記憶手段に保存する地上開度段彩手段、前記記憶手段に保存された各画像データからそれぞれ画像を生成し、重ね合わせて植生図として端末に表示する地図合成手段としての機能を実行させるための植生図作成プログラム。前記コンピュータに、前記反射強度用カラーテーブル、前記樹高用カラーテーブルおよび前記地上開度用カラーテーブルを、前記端末からの入力により変更する合成画像調整手段としての機能をさらに実行させるための請求項6記載の植生図作成プログラム。前記合成画像調整手段に、ユーザが表示したい画像の種類を前記端末から設定した画像選択表を記憶手段に保存し、前記地図合成手段に、前記画像選択表を参照し、選択された画像のみを前記端末に表示することを実行させるための請求項7記載の植生図作成プログラム。前記コンピュータに、前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度値変換テーブルを前記反射強度値で検索し、前記反射強度値に対して設定されたa成分値に変換し、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高値変換テーブルを前記樹高値で検索し、前記樹高値に対して設定されたb成分値に変換し、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度値変換テーブルを前記地上開度値で検索し、前記地上開度値に対して設定されたL成分値に変換し、前記a成分値、前記b成分値、前記L成分値を前記XY座標値とともにLabカラー植生画像データとして前記記憶手段に保存するLabカラー段彩手段部としての機能をさらに実行させるための請求項6乃至8のいずれか1項に記載の植生図作成プログラム。計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成プログラムであって、コンピュータに、前記航空レーザ測量データから生成された前記格子のXY座標値と地盤標高値からなるDEMデータと、前記XY座標値と表層標高値からなるDSMデータと、前記XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力手段、前記記憶手段に保存された前記DEMデータと前記DSMデータから、同一XY座標値に関する前記地盤標高値と前記表層標高値を読み出し、前記表層標高値から前記地盤標高値を差し引いた樹高値を、前記XY座標値とともにDHMデータとして前記記憶手段に保存するDHMデータ作成手段、前記DHMデータの任意の格子を標本格子と設定し、前記標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、前記最大傾斜角度を平均化した地上開度値を、前記標本格子の前記XY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算手段、前記記憶手段に保存された前記反射強度データから前記反射強度値を読み出し、予め前記記憶手段に保存しておいた反射強度値変換テーブルを前記反射強度値で検索し、前記反射強度値に対して設定されたa成分値に変換し、前記DHMデータから前記樹高値を読み出し、予め前記記憶手段に保存しておいた樹高値変換テーブルを前記樹高値で検索し、前記樹高値に対して設定されたb成分値に変換し、前記DHM地上開度データから前記地上開度値を読み出し、予め前記記憶手段に保存しておいた地上開度値変換テーブルを前記地上開度値で検索し、前記地上開度値に対して設定されたL成分値に変換し、前記a成分値、前記b成分値、前記L成分値を前記XY座標値とともにLabカラー植生画像データとして前記記憶手段に保存するLabカラー段彩手段、前記Labカラー植生画像データから画像を生成し、植生図として端末に表示する地図合成手段としての機能を実行させるための植生図作成プログラム。オルソフォト画像は視覚的効果が優れている一方で、次のような欠点がある。1オルソフォト画像は上空から撮影するために、地形や地物などにより影が生じるので、見えない個所がある。2オルソフォト画像は対象地域内を一定の面積で順番に撮影した複数の画像であるため、隣接する画像を接合する際に接合部にずれが生じたり、見えない場所が発生したりする。また、各画像を撮影する際の明るさ等の条件が時々刻々変化する中で撮影されるため、各画像に色調の違いがあり、色による解析などでは誤差が発生し得る。3オルソフォト画像は地物に対する撮影方向によっては、地物の倒れ込み等により、位置ずれや見えない場所が発生する。4オルソフォト画像は平面的なカラー写真であるために、山岳部の河川や谷が樹木に覆っている場合には、河川や谷を容易に判断できない。また、山岳部の尾根も、尾根に生えている樹木の殆どが同系色例えば緑である場合は、尾根であることを容易に判断できない。すなわち、平面的なカラー写真画像では立体感に欠けると共に、表層下を容易に把握することができない。一方、航空レーザ測量データは、なんらかの方法で色情報を付加しないと、非特許文献1に掲載されている画像のように、作成される画像は白黒の画像となり、視覚的な分析効率が悪い。本発明は、以上の課題を解決するためになされたもので、上述した欠点を持つオルソフォト画像を用いずに、航空レーザ測量データのみから樹木の分布を一目で把握することができる植生図を作成する植生図作成装置および植生図作成プログラムを得ることを目的とする。本発明の植生図作成装置は、計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成装置である。航空レーザ測量データから生成された格子のXY座標値と地盤標高値からなるDEMデータと、XY座標値と表層標高値からなるDSMデータと、XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力部と、記憶手段に保存されたDEMデータとDSMデータから、同一XY座標値に関する地盤標高値と表層標高値を読み出し、表層標高値から地盤標高値を差し引いた樹高値を、XY座標値とともにDHMデータとして記憶手段に保存するDHMデータ作成部と、DHMデータの任意の格子を標本格子と設定し、標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、最大傾斜角度を平均化した地上開度値を、標本格子のXY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算部と、記憶手段に保存された反射強度データから反射強度値を読み出し、予め記憶手段に保存しておいた反射強度用カラーテーブルを反射強度値で検索し、反射強度値に対して設定された第1の表示色成分を取得して、XY座標値とともに反射強度画像データとして記憶手段に保存する反射強度段彩部と、DHMデータから樹高値を読み出し、予め記憶手段に保存しておいた樹高用カラーテーブルを樹高値で検索し、樹高値に対して設定された第2の表示色成分を取得して、XY座標値とともに樹高画像データとして記憶手段に保存する樹高段彩部と、DHM地上開度データから地上開度値を読み出し、予め記憶手段に保存しておいた地上開度用カラーテーブルを地上開度値で検索し、地上開度値に対して設定された第3の表示色成分を取得して、XY座標値とともに樹冠立体画像データとして記憶手段に保存する地上開度段彩部と、記憶手段に保存された各画像データからそれぞれ画像を生成し、重ね合わせて植生図として端末に表示する地図合成部とを含むことを要旨とする。さらに、本発明の植生図作成装置は、記憶手段に保存された反射強度データから反射強度値を読み出し、予め記憶手段に保存しておいた反射強度値変換テーブルを反射強度値で検索し、反射強度値に対して設定されたa成分値に変換し、DHMデータから樹高値を読み出し、予め記憶手段に保存しておいた樹高値変換テーブルを樹高値で検索し、樹高値に対して設定されたb成分値に変換し、DHM地上開度データから地上開度値を読み出し、予め記憶手段に保存しておいた地上開度値変換テーブルを地上開度値で検索し、地上開度値に対して設定されたL成分値に変換し、a成分値、b成分値、L成分値をXY座標値とともにLabカラー植生画像データとして記憶手段に保存するLabカラー段彩部を含むことを要旨とする。また、本発明の植生図作成プログラムは、計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成プログラムである。コンピュータに、航空レーザ測量データから生成された格子のXY座標値と地盤標高値からなるDEMデータと、XY座標値と表層標高値からなるDSMデータと、XY座標値と反射強度値からなる反射強度データを、外部記憶媒体から記憶手段に取り込む入力手段、記憶手段に保存されたDEMデータとDSMデータから、同一XY座標値に関する地盤標高値と表層標高値を読み出し、表層標高値から地盤標高値を差し引いた樹高値を、XY座標値とともにDHMデータとして記憶手段に保存するDHMデータ作成手段、DHMデータの任意の格子を標本格子と設定し、標本格子から複数方向毎に検索範囲内に存在する対象格子までの格子間距離と樹高差から最大傾斜角度をそれぞれ求め、最大傾斜角度を平均化した地上開度値を、標本格子のXY座標値とともにDHM地上開度データとして記憶手段に保存する地上開度計算手段、記憶手段に保存された反射強度データから反射強度値を読み出し、予め記憶手段に保存しておいた反射強度用カラーテーブルを反射強度値で検索し、反射強度値に対して設定された第1の表示色成分を取得して、XY座標値とともに反射強度画像データとして記憶手段に保存する反射強度段彩手段、DHMデータから樹高値を読み出し、予め記憶手段に保存しておいた樹高用カラーテーブルを樹高値で検索し、樹高値に対して設定された第2の表示色成分を取得して、XY座標値とともに樹高画像データとして記憶手段に保存する樹高段彩手段、DHM地上開度データから地上開度値を読み出し、予め記憶手段に保存しておいた地上開度用カラーテーブルを地上開度値で検索し、地上開度値に対して設定された第3の表示色成分を取得して、XY座標値とともに樹冠立体画像データとして記憶手段に保存する地上開度段彩手段、記憶手段に保存された各画像データからそれぞれ画像を生成し、重ね合わせて植生図として端末に表示する地図合成手段としての機能を実行させることを要旨とする。さらに、本発明の植生図作成プログラムは、計測地域全体に所望の間隔の格子を設定した航空レーザ測量データのみを使用して植生図を作成する植生図作成プログラムである。コンピュータに、記憶手段に保存された反射強度データから反射強度値を読み出し、予め記憶手段に保存しておいた反射強度値変換テーブルを反射強度値で検索し、反射強度値に対して設定されたa成分値に変換し、DHMデータから樹高値を読み出し、予め記憶手段に保存しておいた樹高値変換テーブルを樹高値で検索し、樹高値に対して設定されたb成分値に変換し、DHM地上開度データから地上開度値を読み出し、予め記憶手段に保存しておいた地上開度値変換テーブルを地上開度値で検索し、地上開度値に対して設定されたL成分値に変換し、a成分値、b成分値、L成分値をXY座標値とともにLabカラー植生画像データとして記憶手段に保存するLabカラー段彩手段部としての機能をさらに実行させることを要旨とする。本発明では、上述した欠点を持つオルソフォト画像を用いないため、影の影響、接合部のずれ、区域ごとの色調の違い、地物の倒れ込みを生じることなく、鮮明な植生図を作成することができる。さらに、航空レーザ測量データを元に段彩を施し、樹木の分布が把握できる立体的な植生図を作成することで、植生の色調、樹冠形状、そして樹高の情報を一目で読み取ることができる。また、樹種別に植生域を色分けすることで、機械的に林相区分線を引いた植生図を作成することができ、従来のオルソフォト画像による林相区分の目視判読作業を効率的かつより正確に行うことができる。本発明は、図1に例示したような、オルソフォト画像を用いずに、航空レーザ測量データ20のみから樹木の分布が把握できる立体的な植生図を作成する装置を提供することを目的としている。そこで、航空レーザ測量データ20から得られる樹木の特徴、すなわち、樹種を識別する指標として、樹高、樹冠形状、および葉の組成を選定した。樹種によって、低木であったり高木であったりと樹高が異なり、枝や葉の茂り方で樹冠形状も異なる。さらに、樹種によって、葉の形状や色が異なる。したがって、これらの特徴を航空レーザ測量データ20から数値化して地図に表し、それらを重ね合わせることで、樹種を識別できる植生図を作成することができる。まず、本発明で用いる入力データを、図2にデータ構成を例示しながら定義する。航空レーザ測量データ20とは、航空レーザ測量システムを用いて、航空機に搭載したレーザ測距儀から地上に向けてレーザを照射し、地上からのレーザ反射パルスとの時間差より求めた地上までの距離やレーザ反射パルスの強度を、計測地点の位置とともに計測したデータである。数値標高モデルDigitalElevationModelデータ以下、DEMデータ21というとは、計測地域全体に所望の格子間隔dたとえば0.2mや0.5mなどの格子構造を設定し、航空レーザ測量データ20のうち、レーザ反射パルスのうち主に最後に返ってきたパルスラストパルスによって計測された標高データから、地表面以外の建物や樹木などを取り除くフィルタリングを行い、標高値内挿補間法によって得た地盤の格子状の標高データであり、図2aにデータ構成を例示したように、格子番号ii=1,2,,nを付与した各格子の中心点のX座標経度Xi、Y座標緯度Yi、Z座標地盤標高値Zgiによって構成される。標高値内挿補間法の例としては、航空レーザ測量データ20の同じ標高値を結んだ等高線図を作成し、この等高線図に対して不整三角形網TINを作成して地盤を復元し、TINと各格子の交点の高さを求める方法がある。数値表層モデルDigitalSurfaceModelデータ以下、DSMデータ22というとは、航空レーザ測量データ20のうち、レーザ反射パルスのうち主に最初に返ってきたパルスファーストパルスによって計測された建物や樹木などを含む表層の標高データであり、図2bにデータ構成を例示したように、計測地域全体に設定したDEMデータと同じ格子構造の格子番号iを付与した各格子の中心点のX座標経度Xi、Y座標緯度Yi、Z座標表層標高値Zsiによって構成される。反射強度データ23とは、照射したレーザパルスに対するレーザ反射パルスの強さを計測したもので、図2cにデータ構成を例示したように、計測地域全体に設定したDEMデータと同じ格子構造の格子番号iを付与した各格子の中心点のX座標経度Xi、Y座標緯度Yi、反射強度値Fiによって構成される。以上のDEMデータ21、DSMデータ22、および反射強度データ23により、計測地域全体の各格子のXY座標値と地盤標高値Zgi、表層標高値Zsi、および反射強度値Fiを得ることができる。得られたデータからどのように樹種を識別するための指標を数値化し、数値化された指標から植生図を作成するかを、具体的な実施例を示しながら以下に説明する。なお、格子とは1辺の長さが格子間隔dの区画であるが、以下の説明において、格子の各種の値とは、格子のデータとして設定された格子の中心点の座標値や各種の値を示すものとする。<実施の形態1>本発明における実施の形態1の第1の植生図作成装置10には、図1に示すように、データ処理部として、入力部11と、DHMデータ作成部12と、地上開度計算部13と、反射強度段彩部14と、樹高段彩部15と、地上開度段彩部16と、地図合成部17と、合成画像調整部18が含まれる。また、第1の植生図作成装置10の記憶手段図示せずには、DEMデータ21と、DSMデータ22と、反射強度データ23と、数値高度モデルDigitalHeightModelデータ以下、DHMデータ24というと、DHM地上開度データ25と、反射強度用カラーテーブル26と、樹高用カラーテーブル27と、地上開度用カラーテーブル28と、画像選択表29が記憶されている。さらに、第1の植生図作成装置10の記憶手段には、端末19に植生図34を表示するために生成される反射強度画像データ31と、樹高画像データ32と、樹冠立体画像データ33が保存される。入力部11は、定義を上述した航空レーザ測量データ20から得られる入力データであるDEMデータ21と、DSMデータ22と、反射強度データ23を、外部記憶媒体から第1の植生図作成装置10の記憶手段に取り込む処理を行う。入力部11は、CD-ROMやUSBメモリなどの外部記憶媒体からデータを読込む手段であってもよいし、第1の植生図作成装置10に通信ネットワークを介して接続された外部装置の記憶媒体から、通信ネットワーク経由でダウンロードする手段であってもよい。DHMデータ作成部12は、入力部11によって記憶手段に保存されたDEMデータ21とDSMデータ22から、それぞれ各格子の地盤標高値Zgiと表層標高値Zsiを読み出し、表層標高値Zsiから地盤標高値Zgiを差し引いた高度値を、図3aに例示したように、各格子のXY座標値とともにDHMデータ24として記憶手段に保存する。本発明では、樹木を対象としているので、表層標高値Zsiから地盤標高値Zgiを差し引いた高度値は各樹木の樹高値Zhiとなり、DHMデータ24を樹高データとも呼ぶ。地上開度計算部13は、DHMデータ作成部12によって記憶手段に保存されたDHMデータ24から各格子の地上開度値θiを計算して、図3bに例示したように、各格子のXY座標値と地上開度値θiをDHM地上開度データ25として記憶手段に保存する。ここで、開度とは地形特徴線抽出の手法のひとつで、ある地点が周囲に比べて地上に突き出ている程度及び地下に食い込んでいる程度を数値化し、地形を3次元表現する際に利用される。開度には地上開度と地下開度の2種類がある。地上開度とは、ある地点から見たときに、空がどれだけ広く見えるかを数値化したもので、地下開度とは、ある地点から見たときに、地下空間がどれだけ広がっているかを数値化したものである。本発明では、地上開度を各樹木の樹冠の傾斜を示す指標として利用する。たとえば、針葉樹のように枝の横への広がりが少ない樹種は、樹冠が急傾斜であり、地上開度が小さくなる。一方、枝を横に広げる広葉樹のような樹種は、樹冠の傾斜がなだらかであり、地上開度が大きくなる。地上開度計算の詳細は後述する。反射強度段彩部14は、入力部11によって記憶手段に保存された反射強度データ23から各格子の反射強度値Fiを読み出し、予め記憶手段に保存しておいた反射強度用カラーテーブル26図4a参照に設定されている反射強度値Fiと表示色成分Ciの対応関係に従って、各格子の段彩を行い、計測地区の反射強度画像データ31を生成し、記憶手段に保存する。樹高段彩部15は、DHMデータ作成部12によって記憶手段に保存されたDHMデータ24から各格子の樹高値Zhiを読み出し、予め記憶手段に保存しておいた樹高用カラーテーブル27図4b参照に設定されている樹高値Zhiと表示色成分Ciの対応関係に従って、各格子の段彩を行い、計測地区の樹高画像データ32を生成し、記憶手段に保存する。地上開度段彩部16は、地上開度計算部13によって記憶手段に保存されたDHM地上開度データ25から各格子の地上開度値θiを読み出し、予め記憶手段に保存しておいた地上開度用カラーテーブル28図4c参照に設定されている地上開度値θiと表示色成分Ciの対応関係に従って、各格子の段彩を行い、計測地区の樹冠立体画像データ33を生成し、記憶手段に保存する。記憶手段に保存された反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33は、図5に例示したように、各格子のXY座標値と表示色成分Ciから構成される。表示色成分Ciのデータは、各カラーテーブルで定義されているデータであり、RGB値のように3要素であったり、グレースケール値のみの1要素であったりする。地図合成部17は、反射強度段彩部14、樹高段彩部15、および地上開度段彩部16によってそれぞれ記憶手段に保存された反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33から、各格子のXY座標値と表示色成分Ciを読み出し、各画像を生成し、端末19に重ね合わせて表示して、植生図34を作成する。また、図示していないが、地図合成部17は、接続されたプリンタに反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33を重ね合わせた植生図34のデータを出力して印刷することもできる。合成画像調整部18は、反射強度段彩部14、樹高段彩部15、および地上開度段彩部16が各格子の段彩を行う際に使用する反射強度用カラーテーブル26、樹高用カラーテーブル27、および地上開度用カラーテーブル28を、端末19からの入力により記憶手段に保存し、必要に応じて変更する。さらに、合成画像調整部18は、植生図34を作成するための反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33のうち、ユーザが表示したい画像データのみを端末19から選択させ、画像選択表29として記憶手段に保存する。地図合成部17は、画像選択表29を参照し、選択された画像のみを端末19に表示するようにする。次に、図4a~cに例示したカラーテーブルの構成について説明する。図4aに示した反射強度用カラーテーブル26の例では、反射強度値Fiを最弱値から最強値までを一定間隔ごとに分割し、表示色成分Ciとして、どのような成分構成かを示す成分モードM本例では、RGBと、それぞれの値に対してどの色を割当てるかをRGBの各値によって設定している。たとえば、実施の形態1では、反射強度値Fiが強いほど濃い青となるように、反射強度値Fiの最弱値から最強値までの範囲に、反射強度用カラーテーブル26の横のカラースケールに表色した白から青のグラデーションが割当たるよう、テーブル内の各反射強度値Fiに対してRGBの各値をそれぞれ設定しておく。図4bに示した樹高用カラーテーブル27の例では、樹高値Zhiの最低値から最高値までを一定間隔ごとに分割し、表示色成分Ciとして、どのような成分構成かを示す成分モードM本例では、RGBと、それぞれの値に対してどの色を割当てるかをRGBの各値によって設定している。たとえば、実施の形態1では、樹高が高いほど濃い緑となるように、樹高値Zhiの最低値から最高値までの範囲に、樹高用カラーテーブル27の横のカラースケールに表色した白から緑のグラデーションが割当たるよう、テーブル内の各樹高値Zhiに対してRGBの各値をそれぞれ設定しておく。図4cに示した地上開度用カラーテーブル28の例では、樹冠立体画像データ33をグレースケールで表示するために、地上開度値θiの最小値から最大値までを一定間隔ごとに分割し、表示色成分Ciとして、どのような成分構成かを示す成分モードM本例では、グレースケールと、それぞれの値に対してどの程度の灰色とするかをグレースケール値によって設定している。たとえば、実施の形態1では、地上開度が大きいほど濃い灰色となるように、地上開度値θiの最小値から最大値までの範囲に、地上開度用カラーテーブル28の横のカラースケールに表色した白から黒のグラデーションが割当たるよう、テーブル内の各地上開度値θiに対してグレースケール値をそれぞれ設定しておく。実施の形態1では、RGBやグレースケールによって青系、緑系、およびグレースケールの表示色成分Ciを指定しているが、使用目的や航空レーザ計測した地域の状態に応じて、表示色成分Ciを青系、緑系、およびグレースケール以外の組み合わせに変更できる。次に、本発明における実施の形態1の第1の植生図作成装置10の処理の流れについて、図6~図11を参照しながら説明する。DHMデータ作成部12の処理入力部11による入力データの記憶手段への保存が完了すると、DHMデータ作成部12は、図6に示す手順によりDHMデータ24を作成する。まず、記憶手段に保存されたDEMデータ21を読み込むステップS12-1。また、記憶手段に保存されたDSMデータ22を読み込むステップS12-2。次に、DEMデータ21から第1の格子格子番号i=1のXY座標値X1,Y1と地盤標高値Zg1を取得するステップS12-3。取得したXY座標値X1,Y1をキー情報として、DSMデータ22から同一のXY座標値X1,Y1をもつ格子を検索し、表層標高値Zs1を取得するステップS12-4。つづいて、取得した表層標高値Zs1から地盤標高値Zg1を差し引き、樹高値Zh1を算出するステップS12-5。最後に、算出した樹高値Zh1を、第1の格子のXY座標値X1,Y1とともに、DHMデータ24として記憶手段に保存するステップS12-6。ここで、処理を行った格子がDEMデータ21内の最後の格子であるかを判定する。たとえば、図2aのようにDEMデータ21に各格子のデータとともに全格子数最大格子番号i=nを設定しておき、処理を行った格子の格子番号iがnになったか判定し、nと等しい場合は、全格子について処理が完了したことになるので、処理を終了するステップS12-7。一方、処理を行った格子の格子番号iがnより小さい場合は、格子番号iに1を加算することで次の格子を指定し、ステップS12-3の処理から繰り返すステップS12-8。ステップS12-3からステップS12-6の処理を繰り返すことで、格子番号i=2、3、、nの樹高値Zh2、Zh3、、Zhnを算出し、図3aに図示したDHMデータ24が完成する。地上開度計算部13の処理本発明では、各樹木の樹冠の傾斜を示す指標として利用する地上開度値θiは、図8に図示したように、格子番号iの格子以下、地上開度計算部13の処理では標本格子Piというと、標本格子Piの地点から見通し距離以内の範囲に存在する各格子以下、地上開度計算部13の処理では対象格子Qmというとの間に形成される傾斜角から算出する。ここで、見通し距離は標本格子Piを中心とする円の半径Rで表し、その円内が検索範囲となる。半径Rが樹木間隔より長過ぎると、隣木の頂点の格子が検索範囲に含まれてしまい、隣木の頂点に対して最大傾斜角をなしてしまう。また、半径Rが樹木間隔に対して短過ぎると、樹冠内の山谷形状を検出し、一つの樹冠に対して複数の樹木頂点が抽出される場合がある。そのため、半径Rは計測地域の代表的な樹木間隔未満に設定するのが良い。標本格子Piの地上開度値θiは、標本格子Piから8方向に存在する対象格子Qmを検索し、各対象格子Qmとの間で形成される天頂からの角度θmの中から各方向について最大傾斜角度θmaxを抽出し、その合計を方向数で割って算出した平均傾斜角度である。図8aのXY平面上の格子位置の図に示したように、8方向をXY平面上でX軸となす角度によって、D0、D45、D90、D135、D180、D225、D270、D315とする。図8aに示した◆印が標本格子Pi、●印が各方向線上に位置する対象格子Qmであり、○印は対象外の格子を示している。たとえ各方向線上に位置していても検索範囲外に存在する格子は対象外となっている。地上開度計算部13は、DHMデータ作成部12によるDHMデータ24の作成が完了すると、DHM地上開度データ25を作成するが、図8bのXYZ空間の格子位置の図を参照しながら、図7に示したフローチャートにしたがって処理を説明する。まず、地上開度計算部13は、記憶手段に保存されたDHMデータ24を読み込むステップS13-1。次に、上述したように、検索範囲を指定するために、計測地域の代表的な樹木間隔未満となる半径Rを設定するステップS13-2。読み込んだDHMデータ24から、地上開度値θiを計算する格子番号iの標本格子Piを、たとえば格子番号iの最小値によって設定するステップS13-3。標本格子PiのXY座標値Xi,Yiと格子間隔dから、1方向の最初の対象格子QmのXY座標値Xm,Ymを、検索する方向によって以下のように設定するステップS13-4。方向D0の場合:Xm,Ym=Xid,Yi方向D45の場合:Xm,Ym=Xid,Yid方向D90の場合:Xm,Ym=Xi,Yid方向D135の場合:Xm,Ym=Xi-d,Yid方向D180の場合:Xm,Ym=Xi-d,Yi方向D225の場合:Xm,Ym=Xi-d,Yi-d方向D270の場合:Xm,Ym=Xi,Yi-d方向D315の場合:Xm,Ym=Xid,Yi-dステップS13-4で設定したXY座標値Xm,YmをキーにDHMデータ24を検索し、該当の対象格子QmがDHMデータ24に存在するかを判定するステップS13-5。該当の対象格子Qmが存在した場合、標本格子Piから対象格子QmまでのXY平面上の格子間距離Lmを、式1によって算出するステップS13-6。Lm={Xm-Xi2Ym-Yi2}1/2式1次に、対象格子Qmが検索範囲内かを判定するために、算出した格子間距離Lmが半径R以下か比較するステップS13-7。対象格子Qmが検索範囲内であれば、標本格子Piと対象格子Qmの樹高差△Zhmを式2によって算出するステップS13-8。△Zhm=Zhm-Zhi式2算出した格子間距離Lmと樹高差△Zhmから標本格子Piと対象格子Qmのアスペクト比Kmを式3によって算出し、記憶手段に一時的に記憶しておくステップS13-9。Km=△Zhm/Lm式31か所の対象格子Qmについて処理を完了すると、検索する方向によって、以下のように対象格子QmのXY座標値Xm,Ymを更新してステップS13-5に戻り、次の対象格子Qmについて処理を繰り返すステップS13-10。方向D0の場合:Xm,Ym=Xmd,Ym方向D45の場合:Xm,Ym=Xmd,Ymd方向D90の場合:Xm,Ym=Xm,Ymd方向D135の場合:Xm,Ym=Xm-d,Ymd方向D180の場合:Xm,Ym=Xm-d,Ym方向D225の場合:Xm,Ym=Xm-d,Ym-d方向D270の場合:Xm,Ym=Xm,Ym-d方向D315の場合:Xm,Ym=Xmd,Ym-d処理の途中、ステップS13-5とステップS13-7の判定で、検索を行っている方向に該当する対象格子Qmが存在しなくなった場合は、標本格子Piに関して全方向について処理が完了したかを、たとえば方向がD315となったか比較して判定するステップS13-11。まだ検索すべき方向がある場合は、次の検索方向を設定してステップS13-4に戻り、処理を繰り返すステップS13-12。標本格子Piに関して全方向について処理が完了した場合は、ステップS13-9によって記憶されていた全アスペクト比Kmから各方向の最大値Kmaxを抽出し、式4、式5によって各方向の最大傾斜角度θmaxを算出し、さらに各方向の最大傾斜角度θmaxの合計を対象格子Qmが存在した方向数で割った平均値を地上開度値θiとし、標本格子PiのXY座標値とともにDHM地上開度データ25として記憶手段に保存するステップS13-13。Max_angle=tan-1Kmax式4最大傾斜角度θmax=90-Max_angle式5ここで、対象格子Qmが存在した方向数は、標本格子Piが計測地域のどの位置に存在するかにより異なり、計測地域の内側に存在する場合は、図8a示した全8方向となる。一方、標本格子Piが計測地域の縁に存在する場合、計測地域の縁の向こう側には他の格子が存在しないために、対象格子Qmが存在した方向数は、計測地域外を指し示す3方向を除いた5方向のみとなり、標本格子Piが計測地域の4角に存在する場合は、計測地域外を指し示す5方向を除いた3方向のみとなる。対象格子Qmが存在した方向数は、たとえばステップS13-5で該当の対象格子QmがDHMデータ24に存在するかを判定し、存在した場合にカウントしておく。最後に、DHMデータ24内の全格子について地上開度値θiの算出が済んだか、たとえば標本格子Piの格子番号iが最大値か判定し、全格子について完了した場合は、処理を終了するステップS13-14。さらに処理すべき格子がある場合は、格子番号iに1を加算して次の標本格子Piに更新して、ステップS13-4から処理を繰り返すステップS13-15。次に、入力部11により記憶手段へ保存された反射強度データ23、DHMデータ作成部12により作成されたDHMデータ24、および地上開度計算部13により作成されたDHM地上開度データ25のそれぞれから、各種画像データを作成する手順について説明する。各種画像データ作成手順の流れを説明するにあたり、図9に反射強度段彩部14の処理のフローチャートを示した。樹高段彩部15と地上開度段彩部16の処理も、図1に示すように参照するデータが異なるのみで、基本的に図9のフローチャートと同様である。反射強度段彩部14の処理図9に示すように、反射強度段彩部14は、入力部11により記憶手段へ保存された反射強度データ23を読み込むステップS14-1。読み込んだ反射強度データ23から、たとえば格子番号iが最小の格子を指定して、該当する格子のXY座標値と反射強度値Fiを取得するステップS14-2。次に、取得した反射強度値Fiをキーにして反射強度用カラーテーブル26を検索し、取得した反射強度値Fiが該当する表示色成分Ciを取得するステップS14-3。取得した表示色成分Ciを、ステップS14-2で取得した格子のXY座標値とともに、反射強度画像データ31として記憶手段に保存するステップS14-4。最後に、処理を行った格子が反射強度データ23に含まれる最後の格子であるか、たとえば格子番号iが最大値未満であるかを判定するステップS14-5。最後の格子でない場合は、格子番号iを1加算して、反射強度データ23の次の格子を指定し、ステップS14-2に戻り、ステップS14-5までの処理を繰り返すステップS14-6。最後の格子である場合は、図4aに例示した反射強度用カラーテーブル26に設定されている表示色成分Ciがどのような成分構成かを示す成分モードMを、図5に例示したように、反射強度データ23に付加して処理を終了するステップS14-7。樹高段彩部15と地上開度段彩部16の処理前述したように、樹高段彩部15と地上開度段彩部16の処理の流れは、反射強度段彩部14の処理の流れと同様である。ステップS14-1に対応する処理で、反射強度データ23に代えて、樹高段彩部15はDHMデータ作成部12により作成されたDHMデータ24を、地上開度段彩部16は地上開度計算部13により作成されたDHM地上開度データ25を読み込む。ステップS14-2に対応する処理で、格子のXY座標値とともに、反射強度値Fiに代えて、樹高段彩部15は樹高値Zhiを、地上開度段彩部16は地上開度値θiを取得する。ステップS14-3に対応する処理で、樹高段彩部15は、取得した樹高値Zhiをキーにして樹高用カラーテーブル27を検索し、取得した樹高値Zhiが該当する表示色成分Ciを取得する。また、地上開度段彩部16は、取得した地上開度値θiをキーにして地上開度用カラーテーブル28を検索し、取得した地上開度値θiが該当する表示色成分Ciを取得する。ステップS14-4に対応する処理で、取得した格子のXY座標値とともに、取得した表示色成分Ciを、樹高段彩部15は樹高画像データ32として、地上開度段彩部16は樹冠立体画像データ33として記憶手段に保存する。樹高段彩部15と地上開度段彩部16は、反射強度段彩部14と同様に、それぞれ以上の処理を最後の格子まで繰り返す。樹高段彩部15と地上開度段彩部16は、図4bと図4cに例示した表示色成分Ciがどのような成分構成かを示す成分モードMを、図5に例示したように、それぞれ樹高画像データ32と樹冠立体画像データ33に付加する。地図合成部17の処理次に、各段彩部によって作成された画像データをユーザの端末19に表示する地図合成部17の処理について、図10のフローチャートを参照しながら説明する。地図合成部17は、各段彩部によってそれぞれ作成された反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33を記憶手段から読み込むステップS17-1。さらに、地図合成部17は、記憶手段に設定された画像選択表29を読み込むステップS17-2。地図合成部17は、反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33のうち、画像選択表29によってユーザが表示するように選択した画像データから画像を生成し、植生図34として端末19に表示するステップS17-3。地図合成部17は、画像を端末19に表示する際に、画像データに設定されている成分モードMによって表示色成分Ciの各値の意味を判断し、必要に応じて成分モードMの示す表色系から端末19の表色系に変換を行う。地図合成部17は、選択された画像データが複数ある場合は、それらの画像を重ね合わせて端末19に表示する。また、ユーザの要求により、選択された画像データから生成した画像をプリンタに出力するステップS17-4。合成画像調整部18の処理合成画像調整部18は、ユーザが端末19に表示される画像を変更したり、画像の色調を調整したりできるように、画像選択表29と各カラーテーブルの設定を変更する手段であり、図11に処理のフローチャートを示す。合成画像調整部18は、記憶手段に保存された各カラーテーブルと画像選択表29を読み込むステップS18-1。合成画像調整部18は、地図合成部17によって端末19に表示されている植生図34の左右または上下のいずれかに、読み込んだデータを表示するステップS18-2。合成画像調整部18は、ユーザが端末19に表示されている植生図34を元に、画像選択表29の画像選択の変更や、各カラーテーブルの色調の変更を設定すると、その設定データを受信するステップS18-3。各カラーテーブルの色調の変更は、カラーテーブルの各表示色成分Ciの値を1つ1つ入力する他に、図4で示したカラースケールをカラーテーブルの横に表示し、その両端の色を指定するだけで、カラーテーブルの範囲内で均等に分割され、各表示色成分Ciの値をカラーテーブルに自動設定するようにしてもよい。合成画像調整部18は、記憶手段に保存された各カラーテーブルと画像選択表29を、受信した設定データよって更新するステップS18-4。合成画像調整部18によって各カラーテーブルと画像選択表29が更新されると、各段彩部が起動され、更新されたカラーテーブルにしたがって画像データが再作成され、地図合成部17は、画像選択表29で新たな選択された画像を端末19に表示することになる。本実施の形態1にて、端末19には地図合成部17によって、図12に示したように、反射強度画像データ31からは、図12aに示す反射強度画像31aが青色系グラデーションで表示され、樹高画像データ32からは、図12bに示す樹高画像32aが緑色系グラデーションで表示され、樹冠立体画像データ33からは、図12cに示す樹冠立体画像33aがグレースケールで表示される。地図合成部17は、画像選択表29で全ての画像を表示するように選択されていれば、全ての画像を重ね合わせて、図12dに示す3画像を合成した植生図34を表示する。図12dに示す植生図34と比較するため、図12eに同じ地域を撮影したオルソフォト画像を示した。全くオルソフォト画像を用いずに、同様の森林の画像を生成できている。さらに、図12dに示す植生図34では、オルソフォト画像のように影になって色が濃くなっている部分がなく、各樹木が区別できる程度に明確に表示されている。また、オルソフォト画像では同じ緑色に撮影された樹木も、図12dに示す植生図34では、樹種が反映されて違った色調で表示されている。<実施の形態2>上述した実施の形態1では、反射強度、樹高、および地上開度の3指標のそれぞれに対して表示色成分Ciを割当てて各画像を生成したのち、3画像を合成することで植生図34を表示している。一方、本発明の実施の形態2では、反射強度、樹高、および地上開度の3指標によって1つの表示色成分Ciを決定し、1画像の植生図34を直接作成する。そこで、実施の形態2では、国際照明委員会CIECommissionInternationaldeI’Eclariageで規定されるLab表色系を利用して植生図34を作成する。図13aに図示したように、Lab表色系以下、Labカラーというの色空間は、L軸、a軸、およびb軸によって定義され、Lで明度を表し、色相と彩度を示す色度をaとbで表す。aとbは色の方向を示し、aは赤、-aは緑、bは黄、-bは青のそれぞれの方向を示している。実施の形態2では、図13bに図示したようなLabカラー色空間を定義する。つまり、L=0黒からL=255白までのL軸と、a=0緑からa=255赤までのa軸と、b=0青からb=255黄までのb軸によってLabカラー色空間を定義する。ここで、地上開度をL軸、反射強度をa軸、樹高をb軸に割当てることで、3指標により1つの表示色が決定することになる。図14には、本発明における実施の形態2の第2の植生図作成装置10aの構成のうち、図1に示した第1の植生図作成装置10と異なり、以下の説明に必要な構成を示し、第1の植生図作成装置10と同じその他の構成要素は省略している。第2の植生図作成装置10aには、Labカラー段彩部40が含まれ、記憶手段には、反射強度値変換テーブル41と、樹高値変換テーブル42と、地上開度値変換テーブル43と、Labカラー段彩部40により生成されるLabカラー植生画像データ44が保存される。Labカラー段彩部40は、入力部11によって記憶手段に保存された反射強度データ23から各格子の反射強度値Fiを読み出し、予め記憶手段に保存しておいた反射強度値変換テーブル41図15a参照を用いて反射強度値Fiをa成分値に変換する。また、Labカラー段彩部40は、DHMデータ作成部12によって記憶手段に保存されたDHMデータ24から各格子の樹高値Zhiを読み出し、予め記憶手段に保存しておいた樹高値変換テーブル42図15b参照を用いて樹高値Zhiをb成分値に変換する。さらに、Labカラー段彩部40は、地上開度計算部13によって記憶手段に保存されたDHM地上開度データ25から各格子の地上開度値θiを読み出し、予め記憶手段に保存しておいた地上開度値変換テーブル43図15c参照を用いて地上開度値θiをL成分値に変換する。Labカラー段彩部40は、変換して得られたa成分値、b成分値、L成分値を、各格子のXY座標値とともにLabカラー植生画像データ44として記憶手段に保存する。次に、図15a~cに例示した変換テーブルの構成について説明する。図15aに示した反射強度値変換テーブル41の例では、反射強度値Fiの最弱値をa成分値の0、最強値をa成分値の255とし、反射強度値Fiの範囲をa=0~255の範囲に均等に割り当てている。図15bに示した樹高値変換テーブル42の例では、樹高値Zhiの最低値をb成分値の0、最高値をb成分値の255とし、樹高値Zhiの範囲をa=0~255の範囲に均等に割り当てている。図15cに示した地上開度値変換テーブル43の例では、地上開度値θiの最小値をL成分値の0、最大値をL成分値の255とし、地上開度値θiの範囲をL=0~255の範囲に均等に割り当てている。Labカラー段彩部40によって生成されるLabカラー植生画像データ44の構成は、図16に示すように、成分モードMとしてLabが設定され、各格子のXY座標値と、L成分値、a成分値、およびb成分値からなる。次に、Labカラー段彩部40の処理の流れについて、図17と図18を参照しながら説明する。Labカラー段彩部40の処理図17に示すように、Labカラー段彩部40は、入力部11により記憶手段へ保存された反射強度データ23を読み込み、反射強度値変換処理を開始する。ステップS40-1。読み込んだ反射強度データ23から、たとえば格子番号iが最小の格子を指定して、該当する格子のXY座標値と反射強度値Fiを取得するステップS40-2。次に、取得した反射強度値Fiをキーにして反射強度値変換テーブル41を検索し、取得した反射強度値Fiが該当するa成分値を取得するステップS40-3。取得したa成分値を、ステップS40-2で取得した格子のXY座標値とともに、Labカラー植生画像データ44として記憶手段に保存するステップS40-4。最後に、処理を行った格子が反射強度データ23に含まれる最後の格子であるか、たとえば格子番号iが最大値未満であるかを判定するステップS40-5。最後の格子でない場合は、格子番号iを1加算して、反射強度データ23の次の格子を指定し、ステップS40-2に戻り、ステップS40-5までの処理を繰り返すステップS40-6。最後の格子である場合は、樹高値変換処理に移行する。図17に示すように、Labカラー段彩部40は、DHMデータ作成部12により作成されたDHMデータ24を読み込むステップS40-7。読み込んだDHMデータ24から、たとえば格子番号iが最小の格子を指定して、該当する格子のXY座標値と樹高値Zhiを取得するステップS40-8。次に、取得した樹高値Zhiをキーにして樹高値変換テーブル42を検索し、取得した樹高値Zhiが該当するb成分値を取得するステップS40-9。取得したb成分値を、既に記憶手段に保存してあるLabカラー植生画像データ44の、同一XY座標の格子のデータに追加するステップS40-10。最後に、処理を行った格子がDHMデータ24に含まれる最後の格子であるか、たとえば格子番号iが最大値未満であるかを判定するステップS40-11。最後の格子でない場合は、格子番号iを1加算して、DHMデータ24の次の格子を指定し、ステップS40-8に戻り、ステップS40-11までの処理を繰り返すステップS40-12。最後の格子である場合は、地上開度値変換処理に移行する。図18に示すように、Labカラー段彩部40は、地上開度計算部13により作成されたDHM地上開度データ25を読み込むステップS40-13。読み込んだDHM地上開度データ25から、たとえば格子番号iが最小の格子を指定して、該当する格子のXY座標値と地上開度値θiを取得するステップS40-14。次に、取得した地上開度値θiをキーにして地上開度値変換テーブル43を検索し、取得した地上開度値θiが該当するL成分値を取得するステップS40-15。取得したL成分値を、既に記憶手段に保存してあるLabカラー植生画像データ44の、同一XY座標の格子のデータに追加するステップS40-16。最後に、処理を行った格子がDHM地上開度データ25に含まれる最後の格子であるか、たとえば格子番号iが最大値未満であるかを判定するステップS40-17。最後の格子でない場合は、格子番号iを1加算して、DHMデータ24の次の格子を指定し、ステップS40-14に戻り、ステップS40-17までの処理を繰り返すステップS40-18。最後の格子である場合は、処理を終了する。地図合成部17は、実施の形態1で示した反射強度画像データ31、樹高画像データ32、および樹冠立体画像データ33とともに、Labカラー段彩部40が作成したLabカラー植生画像データ44を記憶手段から読み込み、画像選択表29によってユーザが表示するように選択した画像データから画像を生成し、植生図34として端末19に表示する。地図合成部17によって、図19aに示したように、Labカラー植生画像データ44からLabカラー植生画像44aが、植生図34として端末19に表示される。実施の形態1の図12dに示した植生図34と比較すると、各樹木の表示色の差がさらに明確になっている。実際の現地調査では、図19aのLabカラー植生画像44aに円で示した、樹高が高く反射強度が強いために赤色系で表示された箇所には、図19bのような葉が茂ったミズナラが群生し、樹高は同様に高いが反射強度が弱いために緑色系で表示された箇所には、図19cのような葉がまばらなサワグルミが群生していた。さらに、青色系で表示された箇所は、樹高が極めて低い、すなわち地盤が現れている場所であり、樹木と明確に区別できる。たとえば、図20aに示したように、Labカラー植生画像44aと立体地形画像を重ねて植生図34として端末19に表示することで、さらに立体的に樹種の植生を視覚的に認識しやすくなる。立体地形画像として、特許第4272146号公報に開示された立体画像作成装置による方法によって作成した画像を利用できる。すなわち、第1および第2の植生図作成装置10、10aに立体地形画像作成部を設け、地上開度計算部13と同様の処理によって、記憶手段に保存されたDEMデータ21からDEM地上開度値を求め、さらにDEMデータ21上に空気層を押し当てた立体を裏返した反転DEMデータ21からDEM地下開度値を求め、両方の差分をグレースケールで表示し、さらに傾斜を赤色系グラデーションで表示した擬似カラー画像を作成する。作成された擬似カラー画像は、尾根や山頂部分が白っぽく、また谷や窪地が黒っぽく表現され、傾斜が急な部分ほど赤く表現されることで1画像でも立体感があり、赤色立体画像ともよばれる。また、外部の立体画像作成装置の作成された立体地形画像を第1および第2の植生図作成装置10、10aに直接読み込む手段を設けても構わない。本発明では、画像選択表29の選択項目に立体地形画像を追加し、地図合成部17は、立体地形画像が選択されていれば、立体地形画像を他の選択画像に重ね合わせて端末19に表示する。この際、他の画像と重ね合わせることを考慮して、立体地形画像の彩度と明度を弱くして表示するなどの調整を適宜行うものとする。図20aはグレースケールとなっているためわかりにくいが、実際の画像では、円で示した左上部は各樹木が緑色系で表示され、中央部は赤色系で表示され、右下部は赤と緑が混在して表示されていて、実際に現地調査を行うと、それぞれの地区の樹種が異なっている。図20aを元に、同系色の範囲の輪郭を描画していくと、図20bに示したような林相区分図を得ることができる。このように、実施の形態1では、オルソフォト画像のように自然の緑に似た色合いで植生図34を表示することで、オルソフォト画像の代替として利用できる一方、実施の形態2では、樹種がより明確に区別できるように、Lab表色系によって、樹種の指標により自然の樹木の色主に緑以外の赤、青、黄など様々な色で樹木などを表示している。このように、本発明により、樹種の分布が視覚的に認識しやすくなり、現地の植生調査を効率よく行うことができる。以上に示した本実施の形態は、本発明の技術的思想を具体化するための装置や処理手順を例示するものであって、構成部品の配置や処理の順番等を限定するものではない。本発明の技術的思想は、特許請求の範囲に記載された技術的範囲内において、種々の変更を加えることができる。図面は模式的なものであり、装置の構成等は現実のものとは異なることに留意すべきである。"