(57)【要約】【課題】マルチスペクトルセンサ、ハイパースペクトルセンサ等により得られた画像を、より適切に補正できる画像処理装置を得る。【解決手段】処理対象画像を、大気影響変化物による影響の空間分布に基づき画像分割部５が分割する。分割した分割領域ごとに、大気散乱放射輝度計算部７が大気散乱放射輝度を、大気透過率計算部８が大気透過率を算出する。算出した大気散乱放射輝度を基に大気散乱放射輝度補正部１１が、更に、算出した大気透過率を基に大気透過率補正部１２が、処理対象画像の放射輝度を補正する。【選択図】図１
2016126566.tif 000002
【特許請求の範囲】
【請求項１】
入力された画像において、大気散乱及び大気減衰を変化させる物質による影響を検出し、前記物質による影響の空間分布を解析する検出解析部と、
前記検出解析部による解析結果を用いて、前記画像を複数の領域に分割する分割部と、
前記分割部により分割された分割領域ごとに、大気散乱放射輝度の環境依存成分を算出するパラメータを、前記分割領域内の放射輝度の中から選択する解析部と、
前記解析部が選択した前記パラメータと、前記大気散乱放射輝度のバンド比とを用いて前記環境依存成分を算出し、前記環境依存成分と前記大気散乱放射輝度のバンド比とを用いて前記分割領域ごとに前記大気散乱放射輝度を算出する大気散乱放射輝度計算部と、
前記分割領域の前記大気散乱放射輝度と大気透過率とを用いて、前記画像内で当該分割領域に対応する領域内の放射輝度を補正する補正部とを備えることを特徴とする画像処理装置。
【請求項２】
前記大気散乱放射輝度計算部により算出された前記分割領域ごとの前記大気散乱放射輝度を用いて、前記分割領域ごとに前記大気透過率の修正値を算出する大気透過率計算部を備え、
前記補正部は、前記分割領域の前記大気散乱放射輝度と前記大気透過率の修正値とを用いて、前記画像内で当該分割領域に対応する領域内の放射輝度を補正することを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記大気散乱放射輝度計算部により算出された前記分割領域ごとの前記大気散乱放射輝度及び前記大気透過率計算部により算出された前記大気透過率の前記分割領域ごとの前記修正値から、画像内の各画素に対応させる前記大気散乱放射輝度及び前記修正値を補間によって算出する補間部を備えることを特徴とする請求項２記載の画像処理装置。
【請求項４】
前記検出解析部は、前記物質による影響が発生する範囲を示す位置情報を用いて、前記画像において前記物質による影響を検出することを特徴とする請求項１から請求項３のうちのいずれか１項記載の画像処理装置。
【請求項５】
検出解析部が、入力された画像において、大気散乱及び大気減衰を変化させる物質による影響を検出し、前記物質による影響の空間分布を解析する検出解析ステップと、
分割部が、前記検出解析ステップによる解析結果を用いて、前記画像を複数の領域に分割する分割ステップと、
解析部が、前記分割ステップにより分割された分割領域ごとに、大気散乱放射輝度の環境依存成分を算出するパラメータを、前記分割領域内の放射輝度の中から選択する解析ステップと、
大気散乱放射輝度計算部が、前記解析ステップにより選択された前記パラメータと、前記大気散乱放射輝度のバンド比とを用いて前記環境依存成分を算出し、前記環境依存成分と前記大気散乱放射輝度のバンド比とを用いて前記分割領域ごとに前記大気散乱放射輝度を算出する大気散乱放射輝度計算ステップと、
補正部が、前記分割領域の前記大気散乱放射輝度と大気透過率とを用いて、前記画像内で当該分割領域に対応する領域内の放射輝度を補正する補正ステップとを備えることを特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、マルチスペクトルセンサ、ハイパースペクトルセンサ等を用いて得られた画像に対し処理を施すための画像処理装置及び画像処理方法に関するものである。
【背景技術】
【０００２】
衛星、航空機等に搭載されるマルチスペクトルセンサ、ハイパースペクトルセンサ等を用いて地球を観測すると、地表面の分光反射特性の違いにより種々の分類を行うことができる。
しかし、衛星、航空機等に搭載されるこれらのセンサには、大気散乱、大気減衰等の大気伝搬による影響を受けた地表面からの反射光が伝達されるため、地表面の被写体が本来有する分光反射特性とは異なった特性を示す。さらに、大気伝搬による地表面からの反射光の変化量は、大気中に含まれる分子、エアロゾル等の大気の状態によって変化することから、センサが取得する計測値には、ばらつきが生じる。
このため、一般的に、大気伝搬による影響を補正することで、精度を高めた地表面の計測が行われる。例えば、特許文献１の画像補正装置では、複数の異なる波長帯域で撮影して得た画像から操作者が目視によりサンプル部位を選択し、当該サンプル部位の輝度値を利用した画像補正が行われる。しかしながら、特許文献１の画像補正装置は、操作者による目視が必要であり、自動的に画像補正を行うものではない。そこで、特許文献２では、センサにより観測された画像のヒストグラム解析を行うことで算出した最小の放射輝度値を利用することにより、自動的に画像補正を行うことが可能なように、画像処理装置が構成されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００５―１５７５６１号公報
【特許文献２】特開２０１３―２２５２４３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
もや、薄雲等により大気中の粒子分布にむらが生じている場合は、もや、薄雲等の濃淡に応じて大気伝搬による影響が部分的に変化する。
しかしながら、上記特許文献２に記載の画像処理装置では、画像中の各画素の放射輝度値を、一様な補正値を用いて補正しているため、大気伝搬による影響の部分的な変化を考慮した補正は行われず、地表面の計測精度が低下していた。つまり、適切な補正が行われていないという課題があった。
【０００５】
この発明は、上記のような課題を解決するためになされたもので、マルチスペクトルセンサ、ハイパースペクトルセンサ等により観測された画像に対して、より適切な補正を行うことのできる画像処理装置を得ることを目的とする。
【課題を解決するための手段】
【０００６】
この発明に係る画像処理装置は、入力された画像において、大気散乱及び大気減衰を変化させる物質による影響を検出し、物質による影響の空間分布を解析する検出解析部と、検出解析部による解析結果を用いて、画像を複数の領域に分割する分割部と、分割部により分割された分割領域ごとに、大気散乱放射輝度の環境依存成分を算出するパラメータを、分割領域内の放射輝度の中から選択する解析部と、解析部が選択したパラメータと、大気散乱放射輝度のバンド比とを用いて環境依存成分を算出し、環境依存成分と大気散乱放射輝度のバンド比とを用いて分割領域ごとに大気散乱放射輝度を算出する大気散乱放射輝度計算部と、分割領域の大気散乱放射輝度と大気透過率とを用いて、画像内で当該分割領域に対応する領域内の放射輝度を補正する補正部とを備えることを特徴とするものである。
【発明の効果】
【０００７】
この発明によれば、画像に対してより適切な補正が行われる。
【図面の簡単な説明】
【０００８】
【図１】この発明の実施の形態１に係る画像処理装置の構成図である。
【図２】この発明の実施の形態１に係る画像処理装置の構成図である。
【図３】この発明の実施の形態１に係る画像処理装置の処理を示すフローチャートである。
【図４】処理対象画像の彩度分布を示す図である。
【図５】大気影響変化物の分布を示す図である。
【図６】補正量補間部による補間処理の様子を説明する図である。
【図７】この発明の実施の形態２に係る画像処理装置の構成図である。
【発明を実施するための形態】
【０００９】
実施の形態１．
図１に、この発明の実施の形態１に係る画像処理装置の構成図を示す。
画像処理装置は、画像入力部１と、対象検出部２と、対象分布解析部３と、分割方式決定部４と、画像分割部５と、ヒストグラム解析部６と、大気散乱放射輝度計算部７と、大気透過率計算部８と、画像結合部９と、補正量補間部１０と、大気散乱放射輝度補正部１１と、大気透過率補正部１２と、スペクトルデータ解析部１３と、解析結果出力部１４と、センサパラメータ記憶部１５とを有する。
【００１０】
画像入力部１は、バンド数ｍ（ｍは正の整数）のマルチスペクトルセンサ、ハイパースペクトルセンサ等とデータ通信する機能（例えば、通信機器、ネットワークインタフェース）を有し、バンド数ｍのマルチスペクトルセンサ、ハイパースペクトルセンサ等により観測された処理対象画像が、入力される。なお、以下では、画像入力部１に入力される処理対象画像が、バンド数ｍのマルチスペクトルセンサによって観測されたものである場合を例に説明する。
【００１１】
対象検出部２は、画像入力部１により入力された処理対象画像において、大気中の粒子分布にむらを生じさせ大気散乱及び大気減衰が変化する原因となる物質（以下、大気影響変化物と称す）による影響を、検出する。大気影響変化物としては、例えば、もや、薄雲、光化学スモッグ等が挙げられる。
対象分布解析部３は、対象検出部２が検出した大気影響変化物による影響の空間分布を、解析する。
なお、対象検出部２と対象分布解析部３とが、検出解析部を構成する。
【００１２】
分割方式決定部４は、対象分布解析部３が解析した大気影響変化物による影響の空間分布に基づき、画像入力部１により入力された処理対象画像の分割方式を決定する。分割方式とは、分割サイズ、分割数等である。
画像分割部５は、分割方式決定部４が決定した分割方式に基づき、画像入力部１により入力された処理対象画像を、ｎ（ｎは正の整数）個の分割領域に分割する。
なお、分割方式決定部４と画像分割部５とが、分割部を構成する。
【００１３】
ヒストグラム解析部６は、画像分割部５により分割された各分割領域に対して、個別にヒストグラム解析を行う。そして、ヒストグラム解析部６は、大気散乱放射輝度計算部７が大気散乱放射輝度を算出するためのパラメータを分割領域ごとに算出する。
大気散乱放射輝度計算部７は、ヒストグラム解析部６が算出したパラメータと、センサパラメータ記憶部１５に記憶されているパラメータとを用いて、画像分割部５により分割された分割領域ごとに大気散乱放射輝度を算出する。
大気透過率計算部８は、大気散乱放射輝度計算部７が算出した大気散乱放射輝度と、センサパラメータ記憶部１５に記憶されているパラメータとを用いて、画像分割部５により分割された分割領域ごとに大気透過率を算出する。
【００１４】
画像結合部９は、大気散乱放射輝度計算部７が算出した分割領域ごとの大気散乱放射輝度を、画像入力部１により入力された処理対象画像での各分割領域の配置関係と同様に並べて結合する。さらに画像結合部９は、大気透過率計算部８が算出した分割領域ごとの大気透過率を、画像入力部１により入力された処理対象画像での各分割領域の配置関係と同様に並べて結合する。
補正量補間部１０は、画像結合部９により並べられた大気散乱放射輝度及び大気透過率が連続的に変化するように、補間処理を施す。
【００１５】
大気散乱放射輝度補正部１１は、補正量補間部１０が補間した大気散乱放射輝度を用いて、画像入力部１により入力された処理対象画像を補正する。
大気透過率補正部１２は、補正量補間部１０が補間した大気透過率を用いて、大気散乱放射輝度補正部１１により補正された処理対象画像を、更に補正する。
なお、大気散乱放射輝度補正部１１と大気透過率補正部１２とが、補正部を構成する。
【００１６】
スペクトルデータ解析部１３は、大気散乱放射輝度補正部１１及び大気透過率補正部１２により補正された処理対象画像をスペクトル解析し、その解析結果であるスペクトル解析画像を出力する。
解析結果出力部１４は、スペクトルデータ解析部１３が出力したスペクトル解析画像を、不図示のディスプレイに表示する。
センサパラメータ記憶部１５は、マルチスペクトルセンサにより観測される各バンドに依存する固有値である既知の大気透過率、大気散乱放射輝度のバンド比である既知の大気散乱放射輝度比等をセンサパラメータとして記憶している。
【００１７】
対象検出部２、対象分布解析部３、分割方式決定部４、画像分割部５、ヒストグラム解析部６、大気散乱放射輝度計算部７、大気透過率計算部８、画像結合部９、補正量補間部１０、大気散乱放射輝度補正部１１、大気透過率補正部１２、スペクトルデータ解析部１３及び解析結果出力部１４は、それぞれが例えばＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）を実装している半導体集積回路又はワンチップマイコン等で構成される。センサパラメータ記憶部１５は、ハードディスク等の各種メモリで構成される。
【００１８】
なお、上記のように画像処理装置の各部を専用のハードウェアで構成するのではなく、画像処理装置を例えばコンピュータで構成してもよい。この場合、図２に示すように、センサパラメータ記憶部１５をコンピュータの内部メモリ又は外部メモリであるメモリ１００上に構成するとともに、画像入力部１、対象検出部２、対象分布解析部３、分割方式決定部４、画像分割部５、ヒストグラム解析部６、大気散乱放射輝度計算部７、大気透過率計算部８、画像結合部９、補正量補間部１０、大気散乱放射輝度補正部１１、大気透過率補正部１２、スペクトルデータ解析部１３及び解析結果出力部１４の処理内容を記述したプログラムをメモリ１００に格納し、当該コンピュータのプロセッサ２００が当該プログラムを実行するようにすればよい。プロセッサ２００が当該プログラムを実行することで、センサパラメータ記憶部１５以外の画像処理装置の各部が実現される。また、プロセッサ２００及びメモリ１００がそれぞれ複数のプロセッサ及び複数のメモリで構成され、当該複数のプロセッサ及び当該複数のメモリが連携して画像処理装置の各部が実現されるようにしてもよい。
【００１９】
次に、画像処理装置による処理の一例について、図３に示すフローチャートを用いて説明する。
マルチスペクトルセンサは、太陽光が地表面で反射されて生じた反射光を観測する。なお、このときマルチスペクトルセンサは、各バンドの波長帯に対する地表面の被写体本来の放射輝度ではなく、当該本来の放射輝度に大気散乱、大気減衰等の大気伝搬による影響が加わったものを観測する。つまり、本来の放射輝度をＬ’’ｂａｎｄ、大気中の粒子によって太陽光が散乱された光がマルチスペクトルセンサに入射することで発生する放射輝度である大気散乱放射輝度をＬｓｃａｔｔ，ｂａｎｄ、各バンドに依存する大気透過率をτｂａｎｄとすると、マルチスペクトルセンサが観測する各バンドの放射輝度Ｌｉｎｐｕｔ，ｂａｎｄは、以下の式（１）で示される。ここで、ｂａｎｄ＝１〜ｍであり、Ｌｉｎｐｕｔ，ｂａｎｄ、τｂａｎｄ、Ｌ’’ｂａｎｄ、Ｌｓｃａｔｔ，ｂａｎｄは、バンドごとの値をとる。
Ｌｉｎｐｕｔ，ｂａｎｄ＝τｂａｎｄ×Ｌ’’ｂａｎｄ＋Ｌｓｃａｔｔ，ｂａｎｄ
・・・・（１）
【００２０】
マルチスペクトルセンサは、観測した放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを用いてバンドごとに画像を形成する。バンドごとに形成された計ｍ個の各画像は、画像内の各画素に割り当てられる画素データの集合体から成り、各画素データには、対応する画素位置での放射輝度Ｌｉｎｐｕｔ，ｂａｎｄ等が含まれている。このようにして成る画像が、画像処理装置に出力される。
【００２１】
画像入力部１は、マルチスペクトルセンサが観測した計ｍ個の画像を、処理対象画像として取り入れる（ステップＳＴ１）。
続いて、対象検出部２が、画像入力部１により入力された処理対象画像において、大気影響変化物による影響を検出する（ステップＳＴ２）。
このとき、対象検出部２は、計ｍ個のバンドごとの画像のうちの少なくともいずれか１つにおいて、例えば海洋、影等の空間的にほぼ一様の反射率特性を持つ被写体を用いて、大気影響変化物が原因となって生じている変動を検出することにより、大気影響変化物による影響を検出する。処理対象画像がＲＧＢ色空間を用いて表現された画像である場合は、まず下記の式（２）〜（６）によりＨＳＶ色空間に変換する。海洋、影等の領域は、彩度Ｓが低い領域として現れるため、対象検出部２は、閾値以下の彩度を示す部分の占める割合が多い領域を、海洋、影等の領域として検出する。図４には、このようにして検出される海洋の領域Ａ１及び影の領域Ａ２を、点線で囲って示している。図４は、処理対象画像の彩度分布のイメージ図であり、彩度が低いほど濃色で示している。そして、対象検出部２は、海洋、影等の領域内の各画素に割り当てられてほぼ同様の数値を示すはずの各画素データの中で、数値が変動している画素データを、大気影響変化物による影響として検出する。つまり、海洋、影等の領域内で画素データの数値が変動している部分が、大気中に大気影響変化物が存在している部分に相当する。
なお、人工物、土壌、植生等は、そもそもそれらに由来した画素データの数値変動がそれらの領域内で現れる。このため、上記のように海洋、影等の領域を対象に画素データの数値変動を検出することで、人工物、土壌、植生等が原因となって生じている画素データの数値変動は検出せずに、大気影響変化物が原因となって生じている画素データの数値変動のみを検出することができる。図４において、海洋の領域Ａ１及び影の領域Ａ２以外の領域が、人工物、土壌、植生等の領域である。
2016126566.tif 000003
2016126566.tif 000004
2016126566.tif 000005
2016126566.tif 000006
2016126566.tif 000007
【００２２】
続いて、対象分布解析部３が、対象検出部２により検出された大気影響変化物による影響の空間分布を解析する（ステップＳＴ３）。図５は、大気影響変化物による影響を抜き出してその分布を分かりやすく示したものであるが、大気影響変化物による影響の分布は、領域Ａ１，Ａ２だけでなく領域Ａ１，Ａ２以外の領域Ｂにも及ぶ。
例えば、対象分布解析部３は、大気影響変化物の影響によるテクスチャが一様であると仮定して、十分に広い海洋の領域Ａ１を２次元フーリエ変換することで、大気影響変化物の影響によるテクスチャの周波数解析を行う。また、十分な広さを持たず、影の領域Ａ２のように離散的に配置された領域が対象検出部２により検出された場合は、領域Ａ２の周囲に位置する領域Ｂ内の各画素に割り当てられた画素データを、領域Ａ２内の各画素に割り当てられた画素データを用いて補間した後、領域Ａ２，Ｂを同様に２次元フーリエ変換することで、大気影響変化物の影響によるテクスチャの周波数解析を行ってもよい。
あるいは、対象分布解析部３は、対象検出部２により検出された大気影響変化物による影響の強弱を２値化して領域分割し、その中での最小領域のサイズを解析してもよい。これは、画素データの数値変動が大きく、大気影響変化物による影響を強く受けている画素を黒、画素データの数値変動が小さく、大気影響変化物による影響をさほど受けていない画素を白とし、同じ色が連続して隣り合う範囲を１つの領域とした際に、これら全ての領域の中で最もサイズが小さい領域を求める処理に相当する。
【００２３】
続いて、分割方式決定部４が、対象分布解析部３による大気影響変化物による影響の空間分布の解析結果に基づき、画像入力部１により入力された処理対象画像の分割方式を決定する（ステップＳＴ４）。
例えば２次元フーリエ変換による周波数解析結果を用いる場合、対象分布解析部３が算出した大気影響変化物による影響の空間周波数を基に当該影響の分布を解析し、分割サイズ及び分割数を決定する。
あるいは、大気影響変化物による影響の強弱による２値化結果を用いる場合、対象分布解析部３が算出した最小領域のサイズを基に分割サイズ及び分割数を決定する。
【００２４】
続いて、画像分割部５が、分割方式決定部４により決定された分割方式を用いて処理対象画像をｎ個に分割する（ステップＳＴ５）。これにより、画像入力部１が処理対象画像として取り入れた計ｍ個の画像それぞれが、同じ位置でｎ個に分割される。なお、分割された各分割領域には、処理対象画像内においてどの位置の領域であるのかを示す画像内位置情報が付加される。
以降のステップＳＴ７〜ＳＴ１１は、分割領域ｉ（ｉ＝１〜ｎ）ごとに行われ、まず１番目の分割領域ｉ（この場合、ｉ＝１）を用いて処理が行われる（ステップＳＴ６）。なお、上記の通り計ｍ個の画像は同じ位置でそれぞれｎ個に分割されており、１番目の分割領域ｉ（ｉ＝１）は、当該ｍ個の画像ごと、つまり、バンドごとに１つずつ存在していることになる。２番目〜ｎ番目の分割領域ｉについても同様である。
ヒストグラム解析部６は、１〜ｎ番目の全ての分割領域ｉに対して処理が行われたか否かを判定する（ステップＳＴ７）。
全ての分割領域ｉに対して処理が行われた場合（ステップＳＴ７；ＹＥＳ）、後述するステップＳＴ１２の処理に移る。
【００２５】
一方、全ての分割領域ｉに対して処理が行われていない場合（ステップＳＴ７；ＮＯ）、ステップＳＴ８の処理に移る。
ヒストグラム解析部６は、分割領域ｉのヒストグラム解析を行い、分割領域ｉ内の各画素に割り当てられている放射輝度Ｌｓｅｎｓｏｒ,ｂａｎｄの中で、バンドごとに最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）を算出する（ステップＳＴ８）。これは、分割領域ｉ内の低反射物体を参照する処理に相当する。
【００２６】
続いて、大気散乱放射輝度計算部７が、ヒストグラム解析部６により算出された各バンドの最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）と、センサパラメータ記憶部１５に記憶されている大気散乱放射輝度比αｂａｎｄとを用いて、大気中の粒子による太陽光の散乱光がマルチスペクトルセンサに入射することで発生する大気散乱放射輝度のうちのマルチスペクトルセンサの観測環境に依存する成分Ｌｓｃａｔｔ０を、以下の式（７）により算出する。観測環境とは、大気中の粒子等である。
2016126566.tif 000008
そして、大気散乱放射輝度計算部７は、観測環境依存成分Ｌｓｃａｔｔ０と大気散乱放射輝度比αｂａｎｄとを用いて、以下の式（８）により大気中の粒子による太陽光の散乱光がマルチスペクトルセンサに入射することで発生する大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを算出する（ステップＳＴ９）。なお、算出した大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄには、ステップＳＴ５で付加された分割領域ｉの画像内位置情報が紐付けられる。
2016126566.tif 000009
【００２７】
続いて、大気透過率計算部８が、大気散乱放射輝度計算部７により算出された大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを用いて、以下の式（９）によりセンサパラメータ記憶部１５に記憶されている各バンドに依存する固有値である大気透過率τｂａｎｄ０を補正し、大気透過率τｂａｎｄを算出する（ステップＳＴ１０）。式（９）は、大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが増加すると、大気透過率が低下するという特徴を利用したものである。式中のβｂａｎｄは、各バンドに依存する係数であり、センサパラメータ記憶部１５に記憶されている。なお、算出した大気透過率τｂａｎｄには、ステップＳＴ５で付加された分割領域ｉの画像内位置情報が紐付けられる。
2016126566.tif 000010
続いて、大気透過率計算部８は、処理対象を次の分割領域に指定する（ステップＳＴ１１）。
【００２８】
ステップＳＴ７〜ＳＴ１１の処理を繰り返し、全ての分割領域ｉに対して処理が行われた場合（ステップＳＴ７；ＹＥＳ）は、続いて画像結合部９が、大気散乱放射輝度計算部７及び大気透過率計算部８により分割領域ｉごとに算出された大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄ及び大気透過率τｂａｎｄが割り当てられた各分割領域ｉを、画像入力部１により入力された処理対象画像での配置関係に戻して再結合する（ステップＳＴ１２）。つまり例えば、画像結合部９は、大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄに紐付けられている画像内位置情報と同じ画像内位置情報が付加されている分割領域ｉを特定することで、各分割領域ｉに、対応する大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを割り当てていく。そして、対応する大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが割り当てられた各分割領域ｉを、付加されている画像内位置情報を基に、元々の処理対象画像での配置関係と同じに配置する。大気透過率τｂａｎｄに関しても同様である。ステップＳＴ１２での画像結合部９の処理は、大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄで構成された大気散乱放射輝度画像と、大気透過率τｂａｎｄで構成された大気透過率画像を作成する処理に相当する。
【００２９】
続いて、補正量補間部１０が、画像結合部９により配置された大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄ及び大気透過率τｂａｎｄの分布が連続的になるように、補間処理を施す（ステップＳＴ１３）。
例えば、分割領域ｉの中心に位置する画素に当該分割領域ｉの大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを割り当て、隣り合う分割領域の中心間の画素には、線形補間関数等で算出した値を割り当てることとする。図６は、補正量補間部１０による大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄの補間処理を視覚的なイメージ図で示したものである。図６（ａ）は、各分割領域ｉに、対応する大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが割り当てられた大気散乱放射輝度画像である。図６では、区切られた１マスが１つの分割領域ｉに相当し、大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが大きいほど濃色で示している。１つの分割領域ｉ内の各画素には、同一の大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが割り当てられており、分割領域ｉの境界を境にして大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄは不連続に変化する。一方、図６（ｂ）は、図６（ａ）に対して補正量補間部１０が補間処理を施した状態の大気散乱放射輝度画像である。大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄが連続的に変化していることが分かる。
なお、大気透過率τｂａｎｄの補間処理についても同様である。
【００３０】
続いて、大気散乱放射輝度補正部１１が、補正量補間部１０により補間処理が施された大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを用いて、画像入力部１により入力された処理対象画像内の各画素に割り当てられている放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを補正する（ステップＳＴ１４）。大気散乱放射輝度補正部１１は、以下の式（１０）に示すように、処理対象画像内の各画素に割り当てられている放射輝度Ｌｉｎｐｕｔ，ｂａｎｄから、補間処理が施された大気散乱放射輝度画像内の各画素に割り当てられている大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを減算することで、放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを補正した放射輝度Ｌ’ｂａｎｄを算出する。
L’band=Linput,band - Lscatt,band・・・・（１０）  
【００３１】
続いて、大気透過率補正部１２が、補正量補間部１０により補間処理が施された大気透過率τｂａｎｄを用いて、大気散乱放射輝度補正部１１により算出されて処理対象画像内の各画素に割り当てられた放射輝度Ｌ’ｂａｎｄを補正する（ステップＳＴ１５）。大気透過率補正部１２は、以下の式（１１）に示すように、処理対象画像内の各画素に割り当てられた放射輝度Ｌ’ｂａｎｄを、補間処理が施された大気透過率画像内の各画素に割り当てられている大気透過率τｂａｎｄで除算することで、放射輝度Ｌ’ｂａｎｄを補正した放射輝度Ｌ’’ｂａｎｄを算出する。放射輝度Ｌ’’ｂａｎｄは、式（１）で示した本来の放射輝度に相当する。
L’’band = L’band / τband・・・・（１１）
【００３２】
続いて、スペクトルデータ解析部１３が、大気散乱放射輝度補正部１１及び大気透過率補正部１２により算出された放射輝度Ｌ’’ｂａｎｄが各画素に割り当てられた処理対象画像をスペクトル解析して、その解析結果であるスペクトル解析画像を出力する（ステップＳＴ１６）。例えば、バンド間の反射率の違いを利用した植生指標ＮＤＶＩ(Ｎｏｒｍａｌｉｚｅｄ Ｄｉｆｆｅｒｅｎｃｅ Ｖｅｇｅｔａｔｉｏｎ Ｉｎｄｅｘ)、水指標ＮＤＷＩ(Ｎｏｒｍａｌｉｚｅｄ Ｄｉｆｆｅｒｅｎｃｅ Ｗａｔｅｒ Ｉｎｄｅｘ)等を算出して、植生の識別、水域・非水域の識別が可能なスペクトル解析画像として出力する。
【００３３】
続いて、解析結果出力部１４が、スペクトルデータ解析部１３が出力したスペクトル解析画像を、不図示のディスプレイに表示する（ステップＳＴ１７）。
【００３４】
以上のように、この実施の形態１に係る画像処理装置によれば、処理対象画像を大気影響変化物による影響の空間分布に基づき分割し、分割した分割領域ごとに大気散乱放射輝度及び大気透過率を算出して、処理対象画像を補正することとした。従って、もや、薄雲等が起因となる大気伝搬による影響の部分的な変化を考慮した補正が可能であり、マルチスペクトルセンサ、ハイパースペクトルセンサ等が観測した処理対象画像を、適切に補正することができる。
【００３５】
なお、上記では、ヒストグラム解析部６は、分割領域ｉのヒストグラム解析を行い、当該分割領域ｉの中で最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）を算出するものとした。しかしながら、ヒストグラムの中で下位数％を除き、残りのヒストグラムの中から、最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）を算出してもよい。これは、マルチスペクトルセンサ、ハイパースペクトルセンサ等では、受光素子の欠陥、内部回路の不具合等により、実際の入射輝度とは異なる値の画素データが割り当てられる不良画素が画像内に発生するからである。このとき、不良画素に割り当てられる画素データは、輝度がゼロ又は飽和値を示すため、このような不良値を除いたものの中から最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）を算出することで、補正の精度を高めることが可能となる。
例えば、処理対象画像内の不良画素が全体の１％である場合、ヒストグラムの中で下位１％を除き、残りの上位９９％のヒストグラムの中から、最小の放射輝度ｍｉｎ（Ｌｓｅｎｓｏｒ,ｂａｎｄ）を算出する。
【００３６】
また、補正量補間部１０を設けずに、大気散乱放射輝度補正部１１及び大気透過率補正部１２が、大気散乱放射輝度計算部７が算出した各分割領域ｉの大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄ及び大気透過率計算部８が算出した各分割領域ｉの大気透過率τｂａｎｄを用いて、処理対象画像における放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを補正してもよい。この場合でも、処理対象画像は、従来の補正処理を用いる場合よりも適切に補正される。もちろん、補正量補間部１０を設ければ、更に適切な補正が可能となる。
【００３７】
また、大気透過率計算部８を設けずに、大気透過率補正部１２が、センサパラメータ記憶部１５に記憶されている大気透過率τｂａｎｄ０を用いて、処理対象画像における放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを補正してもよい。この場合でも、大気散乱放射輝度補正部１１は各分割領域ｉごとの大気散乱放射輝度Ｌｓｃａｔｔ，ｂａｎｄを用いて放射輝度Ｌｉｎｐｕｔ，ｂａｎｄを補正するので、処理対象画像は、従来の補正処理を用いる場合よりも適切に補正される。もちろん、大気透過率計算部８を設ければ、更に適切な補正が可能となる。
【００３８】
実施の形態２．
図７に、この発明の実施の形態２に係る画像処理装置の構成図を示す。
実施の形態２に係る画像処理装置は、画像入力部１と、対象検出部２と、対象分布解析部３と、分割方式決定部４と、画像分割部５と、ヒストグラム解析部６と、大気散乱放射輝度計算部７と、大気透過率計算部８と、画像結合部９と、補正量補間部１０と、大気散乱放射輝度補正部１１と、大気透過率補正部１２と、スペクトルデータ解析部１３と、解析結果出力部１４と、センサパラメータ記憶部１５と、位置情報記憶部１６とを有する。つまり、実施の形態２に係る画像処理装置は、実施の形態１に係る画像処理装置に、位置情報記憶部１６を追加したものである。
なお、図７において、図１と同一又は相当の部分については同一の符号を付し、その説明を省略又は簡略化する。
【００３９】
位置情報記憶部１６は、大気影響変化物が発生する範囲の位置情報を予め記憶している。例えば、都市部等で見られる光化学スモッグの発生地域、山間部で見られる濃霧発生地域等の位置情報を記憶している。この先見的な位置情報は、例えば緯度及び経度のデータとして表されている。
位置情報記憶部１６は、センサパラメータ記憶部１５と同様に、ハードディスク等の各種メモリで構成される。または、センサパラメータ記憶部１５と同様に、コンピュータの内部メモリ又は外部メモリ上に構成されてもよい。
【００４０】
対象検出部２は、位置情報記憶部１６に記憶された位置情報を用いて、画像入力部１により入力された処理対象画像の中で位置情報に示された範囲のみを対象に、大気影響変化物による影響を検出する。
つまり、図３を用いて説明したステップＳＴ２では、画像入力部１により入力された処理対象画像の全範囲を対象に、海洋、影等の空間的にほぼ一様の反射率特性を持つ領域を検出し、当該領域内での画素データの数値変動を大気影響変化物による影響として検出した。しかしながら、実施の形態２に係る画像処理装置の対象検出部２は、画像入力部１により入力された処理対象画像の中で位置情報に示された範囲のみを対象に、海洋、影等の領域を検出し、当該領域内での画素データの数値変動を大気影響変化物による影響として検出する。このように、処理対象画像に対し、大気影響変化物による影響を検出する必要のない範囲については対象検出部２の処理対象外とすることで、処理にかかる計算量が軽減される。
【００４１】
実施の形態２に係る画像処理装置は、対象検出部２による処理が行われるステップＳＴ２以外は、図３に示した実施の形態１に係る画像処理装置と同様の処理を行うので、その説明については省略する。
【００４２】
以上のように、この実施の形態２に係る画像処理装置によれば、位置情報記憶部１６を備えることで、実施の形態１に示す効果に加えて、処理にかかる計算量を軽減することができる。
【００４３】
なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
【符号の説明】
【００４４】
１  画像入力部、２  対象検出部、３  対象分布解析部、４  分割方式決定部、５  画像分割部、６  ヒストグラム解析部、７  大気散乱放射輝度計算部、８  大気透過率計算部、９  画像結合部、１０  補正量補間部、１１  大気散乱放射輝度補正部、１２  大気透過率補正部、１３  スペクトルデータ解析部、１４  解析結果出力部、１５  センサパラメータ記憶部、１６  位置情報記憶部、１００  メモリ、２００  プロセッサ。
【図１】
2016126566.tif 000011
【図２】
2016126566.tif 000012
【図３】
2016126566.tif 000013
【図４】
2016126566.tif 000014
【図５】
2016126566.tif 000015
【図６】
2016126566.tif 000016
【図７】
2016126566.tif 000017