スペクトル解析装置
【課題】 スペクトル解析により対象物を抽出する場合、従来使用されていた2バンド(波長)では誤差が大きく、3バンド以上利用する場合はバンドの選択が困難で、かつ従来の波長幅が狭いセンサの利用では、スペクトルの微小な変化が抽出できず、利用することができなかった。
【解決手段】 RE波長域よりも短波長の波長域からRE波長域にかけてのスペクトル変化分(反射率を検出する場合は、下に凸のスペクトル)を検出するか、あるいは/及び、RE波長域からRE波長域よりも長波長の波長域にかけてのスペクトル変化分(反射率を検出する場合は、上昇するスペクトル)を検出し、植生領域を解析する
【解決手段】 RE波長域よりも短波長の波長域からRE波長域にかけてのスペクトル変化分(反射率を検出する場合は、下に凸のスペクトル)を検出するか、あるいは/及び、RE波長域からRE波長域よりも長波長の波長域にかけてのスペクトル変化分(反射率を検出する場合は、上昇するスペクトル)を検出し、植生領域を解析する
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、所定の物体の物理的特性を解析するためのスペクトルデータを解析する技術に関する。
【背景技術】
【0002】
衛星画像が広域かつ反復的に撮影できる特徴を生かし、衛星画像から植生領域を抽出し、モニタリングすることが行われてきた。植生は森林や畑領域だけでなく、都市域にも存在しており、定量的な緑覆率の算出が求められている。特に、詳細な植生領域の抽出は、都市域の緑覆率把握、防災、環境等といった観点から重要である。
【0003】
衛星や航空機に搭載される多バンドセンサで撮影すると、植生の特徴である赤領域におけるクロロフィル吸収を利用して、植生の領域を抽出することができる。従来の衛星センサは、可視から近赤外域にかけて4バンドの波長帯を有することが多く、この4バンドの波長に該当する赤領域と近赤外領域の反射率の差を利用していた。例えば、赤と近赤外の2バンドを利用する、簡便なNDVI(『森林リモートセンシング、加藤正人著、(株)日本林業調査会発行、2004年』参照)と呼ばれる指標がよく利用されている。また他の手法として、複数のバンドを利用して、例えば土地被覆分類を作ることで植生領域を抽出する場合もある。
【0004】
今後、新たにバンド数の多い衛星が打ちあがることもあり、より多くのバンド(情報)が利用できるようになってきた。これまでも、航空機等に搭載されるハイパースペクトルセンサ(超多バンドセンサ)は、同様に多くのバンドを持ち、詳細なスペクトル情報を取得してきた。図5を参照してハイパースペクトルセンサとマルチセンサのバンド構成を説明する。従来の衛星画像はマルチスペクトルと呼ばれる波長幅の広い4バンドで構成されている。それに対し、ハイパースペクトルセンサは波長幅の狭い、20バンド以上のバンドで構成されている。ハイパースペクトルセンサで撮られた画像から植生領域を抽出する場合、上記NDVIのように2バンドを利用した植生指標や、全バンドを利用して植生を分類する手法が利用されてきた。なお、今後打ち上げられる衛星には、図5のスーパースペクトルセンサ、またはこれに近いバンド幅、バンド数のセンサである。
【先行技術文献】
【非特許文献】
【0005】
【非特許文献1】(『森林リモートセンシング、加藤正人著、(株)日本林業調査会発行、2004年』
【発明の概要】
【発明が解決しようとする課題】
【0006】
NDVIのように従来の植生指標は2バンドの差分値のみに着目しており、特に人工物(建物の屋根)を誤認識することがよくあった。また、土地被覆分類図を作成して植生領域を抽出する場合、各処理対象画像に対して教師データを与える必要があり、大量の画像を処理することができない。
【0007】
また、ハイパースペクトルセンサのデータを利用する場合、特徴的なバンドを選択することが難しく、またバンドに対する波長幅が狭いため、特定のバンドが対象物体の特性変化を表現することができなかった。
【0008】
そこで、本発明は、前記問題に鑑みてなされたものであり、対象とする領域をスペクトル画像から精度よく、かつ簡便に抽出することを課題とする。
【課題を解決するための手段】
【0009】
本発明は、RE波長域よりも短波長の波長域からRE波長域にかけてのスペクトル変化分(反射率を検出する場合は、下に凸のスペクトル)を検出するか、あるいは/及び、RE波長域からRE波長域よりも長波長の波長域にかけてのスペクトル変化分(反射率を検出する場合は、上昇するスペクトル)を検出し、植生領域を解析するものである。換言すれば、RE波長域を含む、少なくとも3つの波長域のスペクトルを検出して,植生領域を解析するものである。
【0010】
実際には、R,RE,NIRの波長域,または/及びY,R,REの波長域の組み合わせのスペクトルを検出すると、精度良く検出できる。
【0011】
このように、従来用いられることのなかった、RE波長域のスペクトルを用いることで、後に述べるように、植生領域の検出を、エラーが少ない状態で行える。
【発明の効果】
【0012】
本発明によれば、対象領域を高精度かつ簡便に抽出できる。また、従来の植生指標であるNDVI等に比較して、誤差が少なく植生などの対象領域を抽出可能である。さらに、植生の活性度を抽出できる。
【図面の簡単な説明】
【0013】
【図1】本実施形態のスペクトル画像解析装置の機能構成図である。
【図2】本実施形態のスペクトル画像解析装置の物理構成図である。
【図3】本実施形態における処理フローを示す図である。
【図4】代表的な植生スペクトルの図である。
【図5】代表的なスペクトルセンサの波長幅および波長位置を示した図である。
【図6】本実施形態におけるカモフラージュされた物体を抽出する図の説明図である。
【図7】本実施形態におけるカモフラージュされた物体の抽出処理フローを示す図である。
【図8】本実施形態において、地図更新を行う説明図である。
【図9】本実施形態において、差分作成部を追加した、スペクトル画像解析装置101の機能構成図である。
【図10】本実施形態において、地図更新を行う処理のフロー図である。
【図11】代表的な植生のコニファーと緑色のペンキのスペクトルを示した図である。
【図12】本実施形態において、植生の生育状態を抽出する図である。
【発明を実施するための形態】
【0014】
以下、本発明を実施するための最良の形態(以下、「実施形態」という。)について、図面を参照(言及図以外の図面も適宜参照)しながら詳細に説明する。なお、本実施形態において、各蓄積部はハードディスクや計算機上のメモリによって実現される。その代案として、ネットワーク上のハードディスクを用いてもよい。また、各処理部は、コンピュータにおいてプログラムを読み込むことで実行される。その代案として、ハードウェアとソフトウェアの協調によってなされるものであってもよい。
【実施例1】
【0015】
図1、図2、および図4を参照して、本実施形態のスペクトル画像解析装置の構成について説明する。図1は、本実施形態のスペクトル画像解析装置101の機能構成図である。図2は、本実施形態のスペクトル画像解析装置101の物理構成図である。図4は代表的な植生のスペクトル図である。
【0016】
図1に示すように、スペクトル画像解析装置101は、処理部としての波長データ抽出部104、植生指標演算部105、スペクトル画像解析部106、画像表示部108と、記憶部(蓄積部)としてのデータ蓄積部101、演算結果蓄積部102とを備えて構成される。
【0017】
データ蓄積部102は、対象とする画像を蓄積(記憶)する。対象画像は、多バンドセンサで撮影された画像であり、また位置情報が付属されている。対象画像は衛星画像に限らず、航空画像、ヘリコプター画像、無人偵察機画像等でもよい。また、画像だけでなく、地上観測等で得られるスペクトルデータでもよい。また,スペクトルデータとは,反射率だけでなく,放射率,吸光度,放射照度,輝度等も含む。さらに放射率は,反射率+放射率=1の関係からも導くことができる。
【0018】
波長データ抽出部104は、データ蓄積部102の画像から、所望の波長帯に相当するバンドデータを抽出する。後記の植生指標演算部で使用する、少なくとも3バンドのデータを抽出する。ここで、所望の波長と近似した波長帯のデータを抽出してもよいし、複数の波長帯を積分し、所望の波長データとして積分値を用いてもよい。
【0019】
PNVI演算部105は、物体の物理的スペクトル特性(以下,スペクトル特性と記す)に着目して、例えば植生指標を演算する。以下では植生を例にとって説明する。まず、植生のスペクトル特性について、図4を使って説明する。図4の横軸は波長、縦軸は反射率である。バンド名401は、対応する代表波長の名前を表し、スペクトル402、403は代表的な植生のスペクトル(反射特性)を示している。植生スペクトルの形状から、R波長域でスペクトルが下に凸となり、RからRE、NIRにかけてスペクトルが上昇する傾向がわかる。これは、Rの波長域に存在するクロロフィル吸収特徴のためで、クロロフィルを持つ植物特有の形状である。植物の種類、成長状態によってスペクトルの絶対値および形状は異なるが、スペクトル402、403を比較しても上記クロロフィル吸収の特徴は変化しないことが分かる。
【0020】
ここで、3バンドのデータを使用する場合の植生指標をPair−wise Normalized Vegetation Index (以下、PNVIと記載)と名づけ、以下の式により計算する。PNVIは使用する波長帯によって、例えば以下の式で表現する。
【0021】
【数1】
【0022】
ただし、(NIR−RE)<0または(RE−R)<0の時、PNVI1はマイナス符号を付ける。
また,以下の式でも表現できる。
【0023】
【数2】
【0024】
ただし、(Y−RE)<0または(RE−R)<0の時、PNVI2はマイナス符号を付ける。
ここで、NIRは近赤外波長のバンド、REはレッドエッジと呼ばれる730nm付近のバンド、Rは赤領域の660nm付近のバンド、Yは黄色領域の615nm付近のバンドを示す。なお、図4に示した通り、NIRはNIR1,NIR2の2バンドあるが、いずれを用いても良い。ここで、図4の主要バンド領域は、CBが400-450nm、Bが450-520nm、Gが520-600nm、Yが580-650nm、Rが630-690nm、REが710-750nm、NIR1が760-900nm、NIR2が860-1050nmである。
【0025】
PNVI1は植生に特徴的なクロロフィル吸収帯に着目し、赤領域の吸収から近赤外領域にかけてスペクトルが立ちあがることを表現している。また、PNVI2は同様にクロロフィル吸収帯に着目し、黄色とレッドエッジ領域に挟まれた赤領域で、スペクトルが下に凸の特徴を持つことを表現している。
【0026】
従来のNDVI等の植生指標では、緑色の屋根を植生として誤分類することがよくあったが、PNVIでは緑色の屋根を植生以外として区別できる。これについて図11を用いて説明する。
【0027】
図11は、代表的な植生であるコニファーと、緑色のペンキのスペクトルを示している。通常のNDVIはRとNIRの差分をとるため、コニファーと緑色ペンキはどちらも高い値を示す。それに対し、例えばPNVI1では、コニファーのスペクトルがRよりREが高く、REよりNIRが高いため、高い値を示す。緑色ペンキは、RよりREが低く、REよりNIRが高いため、マイナスの値を示す。そのため、NDVIでは区別し難い緑色の屋根を、PNVIを使うことで植生と区別することが容易となる。
【0028】
上記では、PNVI1またはPNVI2をそれぞれ単独で使用する場合を説明した。ただし、植生指標に使用する波長帯はこれに限らない。また、3バンドだけでなく、4バンドまたはそれ以上でもよい。また、PNVI1とPNVI2両方を使用して、更に精度を高めても良い。
【0029】
演算結果蓄積部103は、植生指標演算部105で算出された植生指標結果を蓄積する。
【0030】
対象領域抽出部106は、植生指標演算部105で算出された多値の植生指標算出値から、所望の閾値を設定して植生領域とそれ以外の領域を分けることで、植生領域を抽出する。設定する閾値は手動により任意に定めてもよいし、画像内の情報を利用して、自動的に設定してもよい。抽出した植生領域とそれ以外の領域は2値化し、2値化画像として演算結果蓄積部103に格納する。
【0031】
画像重ね合わせ部107は、対象領域抽出部106で作成された植生の2値化画像と、データ蓄積部102に格納されている原画像を重ね合わせ、植生領域を強調する。植生領域部分のみ赤色で塗り、植生領域をわかりやすくする。
【0032】
画像表示部108は、画像重ね合わせ部107で作成した重ね合わせ画像を表示する。また、データ蓄積部102から読み出した原画像を並べて表示し、植生領域を確認し易くする。さらに、画像の拡大、縮小などを並べた画像同士で連動させ、視認し易くする。
【0033】
なお、上記は、スペクトル特性として、反射率に基づいて記載したが、放射率を検出して演算しても良い。
【0034】
次に、スペクトル画像解析装置101の物理構成について説明する。図2に示すように、スペクトル画像解析装置101は、ディスプレイ201、プリンタ202、演算部203、メモリ204、HDD(Hard Disk Drive)205、キーボード206およびマウス207を備えて構成され、それらがバス208によって接続されている。
【0035】
ディスプレイ201は、入力された画像等のデータを表示する。
【0036】
プリンタ202は、入力されたデータを印刷機等へ出力する。
【0037】
演算部203は、各種の処理エンジン2031を有し、処理命令が入力されると対応するエンジンソフトをメモリ204から読み出し、また、必要があればHDD205から蓄積されているデータを読み出して所定の演算を行う。また、計算結果をメモリ204に出力する。
【0038】
メモリ204は、処理エンジン管理部2041、演算情報管理部2042および設定条件管理部2043を有する。処理エンジン管理部2041は、演算部203での計算に使用する各エンジンのソフトを管理する。演算情報管理部2042は、演算部203で計算された植生指標情報を管理する。設定条件管理部2043は、キーボード206やマウス207から入力された条件を保存し、管理する。
【0039】
HDD205は、データ蓄積部102、演算結果蓄積部103を有し、それぞれが該当データを格納している。各データはID番号「1305」のようにIDで管理されており、効率的にデータにアクセスすることができる。HDD205に対してデータの読み出し命令が入力されると、対象データがアクセスされ読み出される。なお、HDD205はネットワーク上に設置されていてもよい。スペクトル画像解析装置101はネットワークインターフェイスを備え,HDD205など各構成はネットワーク上に設置されてもよい。
【0040】
キーボード206及びマウス207は、各種設定条件を入力するユーザによる操作のための手段である。
【0041】
図3は、本実施形態における処理フローを示す図である。
【0042】
ステップ301で、処理対象の画像入力を受け付け、データ蓄積部102に画像データを蓄積する。
【0043】
ステップ302で、PNVI演算部105は、波長データ抽出部104から抽出された所望のバンドデータを用いて、画素毎にPNVIを算出する。算出した結果は演算結果蓄積部103に格納する。ここで、対象としている画像は,2m×2mの画素サイズであるが,これに限らず1m〜1kmサイズの画素でも対応可能である。
【0044】
ステップ303で、ユーザはキーボード206から閾値を入力し、対象領域抽出部106は、入力された閾値を受けつける。閾値の入力はマウス207から入力してもよい。
【0045】
ステップ304で、対象領域抽出部106は、演算結果蓄積部103に格納されているPNVI算出結果画像を2値化する。
【0046】
ステップ305で、画像重ね合わせ部107でPNVI算出結果画像と元画像の重ね合わせ画像を作成し、画像表示部108でディスプレイ201に表示する。
【実施例2】
【0047】
次に、図6、図7を参照して、スペクトル画像解析作成装置の利用例として、カモフラージュされた物体を抽出する例について説明する。
【0048】
図6は、カモフラージュされた物体が撮影されている画像から、カモフラージュされた物体を抽出する例である。スペクトル画像601には、人工物6011とカモフラージュされた人工物6012、植生6013、道路6014が撮影されている。ステップ701で入力された画像は、画像表示部108でディスプレイ201に表示される。カモフラージュされた人工物6012は、植生6013と同様の色で表現されており、通常のRGB画像では区別がつかない。ステップ702で、このスペクトル画像601に対し、PNVI演算部105により算出した結果がPNVI画像602である。ステップ703で、画像表示部108により表示されたPNVI602画像では、植生と植生以外が区別して表現される。人工物6021とカモフラージュされた人工物6022、道路6024が植生以外と表現され、植生6023は植生として抽出される。ステップ704にて、表示されたスペクトル画像601とPNVI画像602を比較すると、スペクトル画像601では植生と同様の色で表現されていた領域6012が、PNVI画像602では、植生以外の領域6022として表現されている。この違いから、比較結果画像603では、カモフラージュされた物体6032として認識できる。
【実施例3】
【0049】
次に、図8、図9、図10を参照して、スペクトル画像解析作成装置の利用例として、地図更新を行う例について説明する。
【0050】
図8は時系列のスペクトル画像から、建物の建て替えを抽出する例である。スペクトル画像801とスペクトル画像802は同一地域を異なる時間で撮影した画像であり、あらかじめ位置合わせ等の処理が済んでいる画像である。PNVI画像803はスペクトル画像801の画像から算出したPNVI結果画像、PNVI画像804はスペクトル画像802の画像から算出した結果画像である。
【0051】
図9は、スペクトル画像解析装置101に差分画像作成部901を追加した装置図である。また、図10は、地図更新を行う処理ステップである。
【0052】
ステップ1001で時系列画像が入力されると、データ蓄積部102に格納される。
【0053】
ステップ1002で、PNVI演算部105によって対応するスペクトル画像の画素毎にPNVIが算出される。算出した画像はPNVI画像803、804で示される。ステップ1003では、差分画像作成部901で、算出されたPNVI画像803、804の差画像805が作成される。ここで、PNVI画像803、804は対象領域抽出部106によって2値化画像となっているため、単純な引き算処理によって、差画像805を作成することができる。また、対象領域抽出部106の2値化処理を行わず、多値画像のまま差分画像を作成してもよい。その場合は、PNVI値が比較的低い領域を、あいまいな植生領域として表現できる。
【0054】
ステップ1004では、作成した差画像805を表示する。PNVI画像803では植生領域以外であった領域8041、8042は、PNVI画像804では植生として表現されており、建物が取り壊され植生領域になったことがわかる。差画像805では、建物から植生への変化領域8051、8052として図示される。また、植生領域に建物等が建設された場合でも、同様の処理により抽出することができる。
【実施例4】
【0055】
次に、図12を参照して、スペクトル画像解析作成装置の利用例として、植生の生育状態を抽出する例について説明する。
【0056】
図12は、農場における生育状態の異なる圃場の例である。スペクトル画像1201には、生育状態が良い圃場1203と、生育状態が普通の圃場1204と、生育状態が不良の1205が撮影されている。通常のRGB画像で表示すると、全ての圃場は同様の色で表現される。このスペクトル画像1201に対して、PNVI演算部105でPNVI画像1202を算出すると、生育状態に対応したPNVI値が算出される。圃場1206はPNVI値が高く、圃場1207はPNVI値が平均的値、圃場1208ではPNVI値が低くなる。これは、図4を参照して説明すると、植物の生育状態(活性状態)は、REの位置およびNIRの反射の高さ、Rの吸収の深さで表される。生育状態が良いと、NIRの値が高くなり、Rが低くなり、また生育段階によってはREが低くなる。逆に生育状態が不良になると、NIRの値が低くなり、Rの値が高くなり、REが高くなる。これらの傾向は、PNVIの式で表現でき、PNVI値が高いと、すなわち生育状態が良いことを示す。このように、PNVIを利用することで、植物の領域のみならず、植物の生育状態を抽出することが可能となる。
【0057】
以上で本実施形態の説明を終えるが、本発明の態様はこれらに限定されるものではない。
【0058】
例えば、PNVIは植生の特徴を利用したバンドの組み合わせを定式化しているが、植生に限らず、例えば地質、水質、人工物等の特徴を鑑みて、所望の波長を利用することで、対象物体を拡大させることができる。また、PNVIの式は一例であって、スペクトル特徴を抽出する3バンド以上の組み合わせを、例えば直線近似、2次微分などの式で表してもよい。例えば,以下の式で表現できる。
【0059】
【数3】
【0060】
ここで,ρは定数である。
さらに、ハイパースペクトルセンサのように波長幅が狭く、スペクトル形状の微小な差を1バンドでは検知できない場合、複数のバンドを積分する,平均値を取る,などの組み合わせによって疑似的に波長幅の広いデータを作成することで、同様にPNVIを適用することができる。また、スペクトル形状の微小な差を必要としない場合は、適当なバンドの組み合わせを用いてもよい。
【0061】
また応用例として、例えば植生と建物を区別する場合を例にとったが、これに限らず土壌、水域、道路等の区別に利用してもよい。
【0062】
その他、ハードウェア、プログラム等の具体的な構成について、本発明の主旨を逸脱しない範囲で適宜変更が可能である。
【符号の説明】
【0063】
101…スペクトル画像解析装置、102…データ蓄積部、103…演算結果蓄積部、104…波長データ抽出部、105…PNVI演算部、106…対象領域抽出部、107…画像重ね合わせ部、108…画像表示部、201…ディスプレイ、202…プリンタ、203…演算部、204…メモリ、205…HDD、206…キーボード、207…マウス、401…バンド名、402…植生スペクトル1、403…植生スペクトル2、601…スペクトル画像、602…PNVI算出結果図、603…比較結果画像、801…スペクトル画像1、802…スペクトル画像2、803…PNVI算出結果図1、804…PNVI算出結果図2、805…差画像、901…差分画像作成部、1201…スペクトル画像、1202…PNVI算出結果図
【技術分野】
【0001】
本発明は、所定の物体の物理的特性を解析するためのスペクトルデータを解析する技術に関する。
【背景技術】
【0002】
衛星画像が広域かつ反復的に撮影できる特徴を生かし、衛星画像から植生領域を抽出し、モニタリングすることが行われてきた。植生は森林や畑領域だけでなく、都市域にも存在しており、定量的な緑覆率の算出が求められている。特に、詳細な植生領域の抽出は、都市域の緑覆率把握、防災、環境等といった観点から重要である。
【0003】
衛星や航空機に搭載される多バンドセンサで撮影すると、植生の特徴である赤領域におけるクロロフィル吸収を利用して、植生の領域を抽出することができる。従来の衛星センサは、可視から近赤外域にかけて4バンドの波長帯を有することが多く、この4バンドの波長に該当する赤領域と近赤外領域の反射率の差を利用していた。例えば、赤と近赤外の2バンドを利用する、簡便なNDVI(『森林リモートセンシング、加藤正人著、(株)日本林業調査会発行、2004年』参照)と呼ばれる指標がよく利用されている。また他の手法として、複数のバンドを利用して、例えば土地被覆分類を作ることで植生領域を抽出する場合もある。
【0004】
今後、新たにバンド数の多い衛星が打ちあがることもあり、より多くのバンド(情報)が利用できるようになってきた。これまでも、航空機等に搭載されるハイパースペクトルセンサ(超多バンドセンサ)は、同様に多くのバンドを持ち、詳細なスペクトル情報を取得してきた。図5を参照してハイパースペクトルセンサとマルチセンサのバンド構成を説明する。従来の衛星画像はマルチスペクトルと呼ばれる波長幅の広い4バンドで構成されている。それに対し、ハイパースペクトルセンサは波長幅の狭い、20バンド以上のバンドで構成されている。ハイパースペクトルセンサで撮られた画像から植生領域を抽出する場合、上記NDVIのように2バンドを利用した植生指標や、全バンドを利用して植生を分類する手法が利用されてきた。なお、今後打ち上げられる衛星には、図5のスーパースペクトルセンサ、またはこれに近いバンド幅、バンド数のセンサである。
【先行技術文献】
【非特許文献】
【0005】
【非特許文献1】(『森林リモートセンシング、加藤正人著、(株)日本林業調査会発行、2004年』
【発明の概要】
【発明が解決しようとする課題】
【0006】
NDVIのように従来の植生指標は2バンドの差分値のみに着目しており、特に人工物(建物の屋根)を誤認識することがよくあった。また、土地被覆分類図を作成して植生領域を抽出する場合、各処理対象画像に対して教師データを与える必要があり、大量の画像を処理することができない。
【0007】
また、ハイパースペクトルセンサのデータを利用する場合、特徴的なバンドを選択することが難しく、またバンドに対する波長幅が狭いため、特定のバンドが対象物体の特性変化を表現することができなかった。
【0008】
そこで、本発明は、前記問題に鑑みてなされたものであり、対象とする領域をスペクトル画像から精度よく、かつ簡便に抽出することを課題とする。
【課題を解決するための手段】
【0009】
本発明は、RE波長域よりも短波長の波長域からRE波長域にかけてのスペクトル変化分(反射率を検出する場合は、下に凸のスペクトル)を検出するか、あるいは/及び、RE波長域からRE波長域よりも長波長の波長域にかけてのスペクトル変化分(反射率を検出する場合は、上昇するスペクトル)を検出し、植生領域を解析するものである。換言すれば、RE波長域を含む、少なくとも3つの波長域のスペクトルを検出して,植生領域を解析するものである。
【0010】
実際には、R,RE,NIRの波長域,または/及びY,R,REの波長域の組み合わせのスペクトルを検出すると、精度良く検出できる。
【0011】
このように、従来用いられることのなかった、RE波長域のスペクトルを用いることで、後に述べるように、植生領域の検出を、エラーが少ない状態で行える。
【発明の効果】
【0012】
本発明によれば、対象領域を高精度かつ簡便に抽出できる。また、従来の植生指標であるNDVI等に比較して、誤差が少なく植生などの対象領域を抽出可能である。さらに、植生の活性度を抽出できる。
【図面の簡単な説明】
【0013】
【図1】本実施形態のスペクトル画像解析装置の機能構成図である。
【図2】本実施形態のスペクトル画像解析装置の物理構成図である。
【図3】本実施形態における処理フローを示す図である。
【図4】代表的な植生スペクトルの図である。
【図5】代表的なスペクトルセンサの波長幅および波長位置を示した図である。
【図6】本実施形態におけるカモフラージュされた物体を抽出する図の説明図である。
【図7】本実施形態におけるカモフラージュされた物体の抽出処理フローを示す図である。
【図8】本実施形態において、地図更新を行う説明図である。
【図9】本実施形態において、差分作成部を追加した、スペクトル画像解析装置101の機能構成図である。
【図10】本実施形態において、地図更新を行う処理のフロー図である。
【図11】代表的な植生のコニファーと緑色のペンキのスペクトルを示した図である。
【図12】本実施形態において、植生の生育状態を抽出する図である。
【発明を実施するための形態】
【0014】
以下、本発明を実施するための最良の形態(以下、「実施形態」という。)について、図面を参照(言及図以外の図面も適宜参照)しながら詳細に説明する。なお、本実施形態において、各蓄積部はハードディスクや計算機上のメモリによって実現される。その代案として、ネットワーク上のハードディスクを用いてもよい。また、各処理部は、コンピュータにおいてプログラムを読み込むことで実行される。その代案として、ハードウェアとソフトウェアの協調によってなされるものであってもよい。
【実施例1】
【0015】
図1、図2、および図4を参照して、本実施形態のスペクトル画像解析装置の構成について説明する。図1は、本実施形態のスペクトル画像解析装置101の機能構成図である。図2は、本実施形態のスペクトル画像解析装置101の物理構成図である。図4は代表的な植生のスペクトル図である。
【0016】
図1に示すように、スペクトル画像解析装置101は、処理部としての波長データ抽出部104、植生指標演算部105、スペクトル画像解析部106、画像表示部108と、記憶部(蓄積部)としてのデータ蓄積部101、演算結果蓄積部102とを備えて構成される。
【0017】
データ蓄積部102は、対象とする画像を蓄積(記憶)する。対象画像は、多バンドセンサで撮影された画像であり、また位置情報が付属されている。対象画像は衛星画像に限らず、航空画像、ヘリコプター画像、無人偵察機画像等でもよい。また、画像だけでなく、地上観測等で得られるスペクトルデータでもよい。また,スペクトルデータとは,反射率だけでなく,放射率,吸光度,放射照度,輝度等も含む。さらに放射率は,反射率+放射率=1の関係からも導くことができる。
【0018】
波長データ抽出部104は、データ蓄積部102の画像から、所望の波長帯に相当するバンドデータを抽出する。後記の植生指標演算部で使用する、少なくとも3バンドのデータを抽出する。ここで、所望の波長と近似した波長帯のデータを抽出してもよいし、複数の波長帯を積分し、所望の波長データとして積分値を用いてもよい。
【0019】
PNVI演算部105は、物体の物理的スペクトル特性(以下,スペクトル特性と記す)に着目して、例えば植生指標を演算する。以下では植生を例にとって説明する。まず、植生のスペクトル特性について、図4を使って説明する。図4の横軸は波長、縦軸は反射率である。バンド名401は、対応する代表波長の名前を表し、スペクトル402、403は代表的な植生のスペクトル(反射特性)を示している。植生スペクトルの形状から、R波長域でスペクトルが下に凸となり、RからRE、NIRにかけてスペクトルが上昇する傾向がわかる。これは、Rの波長域に存在するクロロフィル吸収特徴のためで、クロロフィルを持つ植物特有の形状である。植物の種類、成長状態によってスペクトルの絶対値および形状は異なるが、スペクトル402、403を比較しても上記クロロフィル吸収の特徴は変化しないことが分かる。
【0020】
ここで、3バンドのデータを使用する場合の植生指標をPair−wise Normalized Vegetation Index (以下、PNVIと記載)と名づけ、以下の式により計算する。PNVIは使用する波長帯によって、例えば以下の式で表現する。
【0021】
【数1】
【0022】
ただし、(NIR−RE)<0または(RE−R)<0の時、PNVI1はマイナス符号を付ける。
また,以下の式でも表現できる。
【0023】
【数2】
【0024】
ただし、(Y−RE)<0または(RE−R)<0の時、PNVI2はマイナス符号を付ける。
ここで、NIRは近赤外波長のバンド、REはレッドエッジと呼ばれる730nm付近のバンド、Rは赤領域の660nm付近のバンド、Yは黄色領域の615nm付近のバンドを示す。なお、図4に示した通り、NIRはNIR1,NIR2の2バンドあるが、いずれを用いても良い。ここで、図4の主要バンド領域は、CBが400-450nm、Bが450-520nm、Gが520-600nm、Yが580-650nm、Rが630-690nm、REが710-750nm、NIR1が760-900nm、NIR2が860-1050nmである。
【0025】
PNVI1は植生に特徴的なクロロフィル吸収帯に着目し、赤領域の吸収から近赤外領域にかけてスペクトルが立ちあがることを表現している。また、PNVI2は同様にクロロフィル吸収帯に着目し、黄色とレッドエッジ領域に挟まれた赤領域で、スペクトルが下に凸の特徴を持つことを表現している。
【0026】
従来のNDVI等の植生指標では、緑色の屋根を植生として誤分類することがよくあったが、PNVIでは緑色の屋根を植生以外として区別できる。これについて図11を用いて説明する。
【0027】
図11は、代表的な植生であるコニファーと、緑色のペンキのスペクトルを示している。通常のNDVIはRとNIRの差分をとるため、コニファーと緑色ペンキはどちらも高い値を示す。それに対し、例えばPNVI1では、コニファーのスペクトルがRよりREが高く、REよりNIRが高いため、高い値を示す。緑色ペンキは、RよりREが低く、REよりNIRが高いため、マイナスの値を示す。そのため、NDVIでは区別し難い緑色の屋根を、PNVIを使うことで植生と区別することが容易となる。
【0028】
上記では、PNVI1またはPNVI2をそれぞれ単独で使用する場合を説明した。ただし、植生指標に使用する波長帯はこれに限らない。また、3バンドだけでなく、4バンドまたはそれ以上でもよい。また、PNVI1とPNVI2両方を使用して、更に精度を高めても良い。
【0029】
演算結果蓄積部103は、植生指標演算部105で算出された植生指標結果を蓄積する。
【0030】
対象領域抽出部106は、植生指標演算部105で算出された多値の植生指標算出値から、所望の閾値を設定して植生領域とそれ以外の領域を分けることで、植生領域を抽出する。設定する閾値は手動により任意に定めてもよいし、画像内の情報を利用して、自動的に設定してもよい。抽出した植生領域とそれ以外の領域は2値化し、2値化画像として演算結果蓄積部103に格納する。
【0031】
画像重ね合わせ部107は、対象領域抽出部106で作成された植生の2値化画像と、データ蓄積部102に格納されている原画像を重ね合わせ、植生領域を強調する。植生領域部分のみ赤色で塗り、植生領域をわかりやすくする。
【0032】
画像表示部108は、画像重ね合わせ部107で作成した重ね合わせ画像を表示する。また、データ蓄積部102から読み出した原画像を並べて表示し、植生領域を確認し易くする。さらに、画像の拡大、縮小などを並べた画像同士で連動させ、視認し易くする。
【0033】
なお、上記は、スペクトル特性として、反射率に基づいて記載したが、放射率を検出して演算しても良い。
【0034】
次に、スペクトル画像解析装置101の物理構成について説明する。図2に示すように、スペクトル画像解析装置101は、ディスプレイ201、プリンタ202、演算部203、メモリ204、HDD(Hard Disk Drive)205、キーボード206およびマウス207を備えて構成され、それらがバス208によって接続されている。
【0035】
ディスプレイ201は、入力された画像等のデータを表示する。
【0036】
プリンタ202は、入力されたデータを印刷機等へ出力する。
【0037】
演算部203は、各種の処理エンジン2031を有し、処理命令が入力されると対応するエンジンソフトをメモリ204から読み出し、また、必要があればHDD205から蓄積されているデータを読み出して所定の演算を行う。また、計算結果をメモリ204に出力する。
【0038】
メモリ204は、処理エンジン管理部2041、演算情報管理部2042および設定条件管理部2043を有する。処理エンジン管理部2041は、演算部203での計算に使用する各エンジンのソフトを管理する。演算情報管理部2042は、演算部203で計算された植生指標情報を管理する。設定条件管理部2043は、キーボード206やマウス207から入力された条件を保存し、管理する。
【0039】
HDD205は、データ蓄積部102、演算結果蓄積部103を有し、それぞれが該当データを格納している。各データはID番号「1305」のようにIDで管理されており、効率的にデータにアクセスすることができる。HDD205に対してデータの読み出し命令が入力されると、対象データがアクセスされ読み出される。なお、HDD205はネットワーク上に設置されていてもよい。スペクトル画像解析装置101はネットワークインターフェイスを備え,HDD205など各構成はネットワーク上に設置されてもよい。
【0040】
キーボード206及びマウス207は、各種設定条件を入力するユーザによる操作のための手段である。
【0041】
図3は、本実施形態における処理フローを示す図である。
【0042】
ステップ301で、処理対象の画像入力を受け付け、データ蓄積部102に画像データを蓄積する。
【0043】
ステップ302で、PNVI演算部105は、波長データ抽出部104から抽出された所望のバンドデータを用いて、画素毎にPNVIを算出する。算出した結果は演算結果蓄積部103に格納する。ここで、対象としている画像は,2m×2mの画素サイズであるが,これに限らず1m〜1kmサイズの画素でも対応可能である。
【0044】
ステップ303で、ユーザはキーボード206から閾値を入力し、対象領域抽出部106は、入力された閾値を受けつける。閾値の入力はマウス207から入力してもよい。
【0045】
ステップ304で、対象領域抽出部106は、演算結果蓄積部103に格納されているPNVI算出結果画像を2値化する。
【0046】
ステップ305で、画像重ね合わせ部107でPNVI算出結果画像と元画像の重ね合わせ画像を作成し、画像表示部108でディスプレイ201に表示する。
【実施例2】
【0047】
次に、図6、図7を参照して、スペクトル画像解析作成装置の利用例として、カモフラージュされた物体を抽出する例について説明する。
【0048】
図6は、カモフラージュされた物体が撮影されている画像から、カモフラージュされた物体を抽出する例である。スペクトル画像601には、人工物6011とカモフラージュされた人工物6012、植生6013、道路6014が撮影されている。ステップ701で入力された画像は、画像表示部108でディスプレイ201に表示される。カモフラージュされた人工物6012は、植生6013と同様の色で表現されており、通常のRGB画像では区別がつかない。ステップ702で、このスペクトル画像601に対し、PNVI演算部105により算出した結果がPNVI画像602である。ステップ703で、画像表示部108により表示されたPNVI602画像では、植生と植生以外が区別して表現される。人工物6021とカモフラージュされた人工物6022、道路6024が植生以外と表現され、植生6023は植生として抽出される。ステップ704にて、表示されたスペクトル画像601とPNVI画像602を比較すると、スペクトル画像601では植生と同様の色で表現されていた領域6012が、PNVI画像602では、植生以外の領域6022として表現されている。この違いから、比較結果画像603では、カモフラージュされた物体6032として認識できる。
【実施例3】
【0049】
次に、図8、図9、図10を参照して、スペクトル画像解析作成装置の利用例として、地図更新を行う例について説明する。
【0050】
図8は時系列のスペクトル画像から、建物の建て替えを抽出する例である。スペクトル画像801とスペクトル画像802は同一地域を異なる時間で撮影した画像であり、あらかじめ位置合わせ等の処理が済んでいる画像である。PNVI画像803はスペクトル画像801の画像から算出したPNVI結果画像、PNVI画像804はスペクトル画像802の画像から算出した結果画像である。
【0051】
図9は、スペクトル画像解析装置101に差分画像作成部901を追加した装置図である。また、図10は、地図更新を行う処理ステップである。
【0052】
ステップ1001で時系列画像が入力されると、データ蓄積部102に格納される。
【0053】
ステップ1002で、PNVI演算部105によって対応するスペクトル画像の画素毎にPNVIが算出される。算出した画像はPNVI画像803、804で示される。ステップ1003では、差分画像作成部901で、算出されたPNVI画像803、804の差画像805が作成される。ここで、PNVI画像803、804は対象領域抽出部106によって2値化画像となっているため、単純な引き算処理によって、差画像805を作成することができる。また、対象領域抽出部106の2値化処理を行わず、多値画像のまま差分画像を作成してもよい。その場合は、PNVI値が比較的低い領域を、あいまいな植生領域として表現できる。
【0054】
ステップ1004では、作成した差画像805を表示する。PNVI画像803では植生領域以外であった領域8041、8042は、PNVI画像804では植生として表現されており、建物が取り壊され植生領域になったことがわかる。差画像805では、建物から植生への変化領域8051、8052として図示される。また、植生領域に建物等が建設された場合でも、同様の処理により抽出することができる。
【実施例4】
【0055】
次に、図12を参照して、スペクトル画像解析作成装置の利用例として、植生の生育状態を抽出する例について説明する。
【0056】
図12は、農場における生育状態の異なる圃場の例である。スペクトル画像1201には、生育状態が良い圃場1203と、生育状態が普通の圃場1204と、生育状態が不良の1205が撮影されている。通常のRGB画像で表示すると、全ての圃場は同様の色で表現される。このスペクトル画像1201に対して、PNVI演算部105でPNVI画像1202を算出すると、生育状態に対応したPNVI値が算出される。圃場1206はPNVI値が高く、圃場1207はPNVI値が平均的値、圃場1208ではPNVI値が低くなる。これは、図4を参照して説明すると、植物の生育状態(活性状態)は、REの位置およびNIRの反射の高さ、Rの吸収の深さで表される。生育状態が良いと、NIRの値が高くなり、Rが低くなり、また生育段階によってはREが低くなる。逆に生育状態が不良になると、NIRの値が低くなり、Rの値が高くなり、REが高くなる。これらの傾向は、PNVIの式で表現でき、PNVI値が高いと、すなわち生育状態が良いことを示す。このように、PNVIを利用することで、植物の領域のみならず、植物の生育状態を抽出することが可能となる。
【0057】
以上で本実施形態の説明を終えるが、本発明の態様はこれらに限定されるものではない。
【0058】
例えば、PNVIは植生の特徴を利用したバンドの組み合わせを定式化しているが、植生に限らず、例えば地質、水質、人工物等の特徴を鑑みて、所望の波長を利用することで、対象物体を拡大させることができる。また、PNVIの式は一例であって、スペクトル特徴を抽出する3バンド以上の組み合わせを、例えば直線近似、2次微分などの式で表してもよい。例えば,以下の式で表現できる。
【0059】
【数3】
【0060】
ここで,ρは定数である。
さらに、ハイパースペクトルセンサのように波長幅が狭く、スペクトル形状の微小な差を1バンドでは検知できない場合、複数のバンドを積分する,平均値を取る,などの組み合わせによって疑似的に波長幅の広いデータを作成することで、同様にPNVIを適用することができる。また、スペクトル形状の微小な差を必要としない場合は、適当なバンドの組み合わせを用いてもよい。
【0061】
また応用例として、例えば植生と建物を区別する場合を例にとったが、これに限らず土壌、水域、道路等の区別に利用してもよい。
【0062】
その他、ハードウェア、プログラム等の具体的な構成について、本発明の主旨を逸脱しない範囲で適宜変更が可能である。
【符号の説明】
【0063】
101…スペクトル画像解析装置、102…データ蓄積部、103…演算結果蓄積部、104…波長データ抽出部、105…PNVI演算部、106…対象領域抽出部、107…画像重ね合わせ部、108…画像表示部、201…ディスプレイ、202…プリンタ、203…演算部、204…メモリ、205…HDD、206…キーボード、207…マウス、401…バンド名、402…植生スペクトル1、403…植生スペクトル2、601…スペクトル画像、602…PNVI算出結果図、603…比較結果画像、801…スペクトル画像1、802…スペクトル画像2、803…PNVI算出結果図1、804…PNVI算出結果図2、805…差画像、901…差分画像作成部、1201…スペクトル画像、1202…PNVI算出結果図
【特許請求の範囲】
【請求項1】
解析対象画像を記憶するデータ蓄積部と、
前記解析対象画像から、所定の波長域に相当するバンドデータを抽出する波長データ抽出部と、
前記抽出されたバンドデータから、植生領域を演算するための演算部とを有し、
前記演算部は、RE波長域を含む少なくとも3つの波長域のスペクトルに基づいて、演算することを特徴とするスペクトル画像解析装置。
【請求項2】
前記少なくとも3つの波長域は、前記RE波長域、R波長域、NIR波長域であることを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項3】
前記演算部は、前記R波長域から前記RE波長域及び前記NIR波長域にかけて、上昇する反射率のスペクトルを用いて演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項4】
前記演算部は、
【数1】
(但し、NIRは近赤外波長のバンド、REはレッドエッジ波長のバンド、Rは赤領域波長のバンドを示す)
を演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項5】
前記少なくとも3つの波長域は、前記RE波長域、Y波長域、R波長域であることを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項6】
前記演算部は、前記Y波長域から前記R波長域及び前記RE波長域にかけて、下に凸の反射率のスペクトルを用いて演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項7】
前記演算部は、
【数2】
(但し、Yは黄色波長のバンド、Rは赤領域波長のバンド、REはレッドエッジ波長のバンドを示す)
を演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項8】
前記演算部は、更に、前記植生領域と植生以外の人工物とを区別し、植生領域を抽出する手段を有することを特徴とするスペクトル画像解析装置。
【請求項9】
更に、データ蓄積部と画像表示部を有し、
前記データ蓄積部は、同一地域を異なる時間で撮影した画像を蓄積しており、
前記演算部は、更に、前記異なる時間で撮影した画像同士の前記植生領域の差分を演算し、
前記画像表示部は、前記差分の画像を表示することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項10】
更に、データ蓄積部と画像表示部を有し、
前記データ蓄積部は、同一地域を異なる時間で撮影した画像を蓄積しており、
前記演算部での演算結果に基づいて、前記画像表示部に、前記植生の生育状態を表示させることを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項1】
解析対象画像を記憶するデータ蓄積部と、
前記解析対象画像から、所定の波長域に相当するバンドデータを抽出する波長データ抽出部と、
前記抽出されたバンドデータから、植生領域を演算するための演算部とを有し、
前記演算部は、RE波長域を含む少なくとも3つの波長域のスペクトルに基づいて、演算することを特徴とするスペクトル画像解析装置。
【請求項2】
前記少なくとも3つの波長域は、前記RE波長域、R波長域、NIR波長域であることを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項3】
前記演算部は、前記R波長域から前記RE波長域及び前記NIR波長域にかけて、上昇する反射率のスペクトルを用いて演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項4】
前記演算部は、
【数1】
(但し、NIRは近赤外波長のバンド、REはレッドエッジ波長のバンド、Rは赤領域波長のバンドを示す)
を演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項5】
前記少なくとも3つの波長域は、前記RE波長域、Y波長域、R波長域であることを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項6】
前記演算部は、前記Y波長域から前記R波長域及び前記RE波長域にかけて、下に凸の反射率のスペクトルを用いて演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項7】
前記演算部は、
【数2】
(但し、Yは黄色波長のバンド、Rは赤領域波長のバンド、REはレッドエッジ波長のバンドを示す)
を演算することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項8】
前記演算部は、更に、前記植生領域と植生以外の人工物とを区別し、植生領域を抽出する手段を有することを特徴とするスペクトル画像解析装置。
【請求項9】
更に、データ蓄積部と画像表示部を有し、
前記データ蓄積部は、同一地域を異なる時間で撮影した画像を蓄積しており、
前記演算部は、更に、前記異なる時間で撮影した画像同士の前記植生領域の差分を演算し、
前記画像表示部は、前記差分の画像を表示することを特徴とする請求項1記載のスペクトル画像解析装置。
【請求項10】
更に、データ蓄積部と画像表示部を有し、
前記データ蓄積部は、同一地域を異なる時間で撮影した画像を蓄積しており、
前記演算部での演算結果に基づいて、前記画像表示部に、前記植生の生育状態を表示させることを特徴とする請求項1記載のスペクトル画像解析装置。
【図1】
【図2】
【図3】
【図6】
【図7】
【図8】
【図9】
【図10】
【図12】
【図4】
【図5】
【図11】
【図2】
【図3】
【図6】
【図7】
【図8】
【図9】
【図10】
【図12】
【図4】
【図5】
【図11】
【公開番号】特開2011−69757(P2011−69757A)
【公開日】平成23年4月7日(2011.4.7)
【国際特許分類】
【出願番号】特願2009−221909(P2009−221909)
【出願日】平成21年9月28日(2009.9.28)
【出願人】(000233055)株式会社日立ソリューションズ (1,610)
【Fターム(参考)】
【公開日】平成23年4月7日(2011.4.7)
【国際特許分類】
【出願日】平成21年9月28日(2009.9.28)
【出願人】(000233055)株式会社日立ソリューションズ (1,610)
【Fターム(参考)】
[ Back to top ]