(57)【要約】  情報処理装置２０は光路輝度補正の精度を高める構成を備える。情報処理装置２０の第１中間光路輝度算出部２１は、画像供給装置３１から供給される波長帯域毎の観測画像を構成する各画素に関連付けられている電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する。第２中間光路輝度算出部２２は、大気の状態を表す情報を含む第２中間光路輝度を算出する。最終光路輝度算出部２３は、波長帯域毎に、観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、第１中間光路輝度および第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する。光路輝度補正部２４は、波長帯域毎に、観測画像の全画素から最終光路輝度を減算する。
2017169946.tif 000002
【特許請求の範囲】
【請求項１】
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する第１中間光路輝度算出手段と、
前記大気の状態を表す情報を含む第２中間光路輝度を算出する第２中間光路輝度算出手段と、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する最終光路輝度算出手段と、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する光路輝度補正手段と、
を備える情報処理装置。
【請求項２】
前記第２中間光路輝度算出手段は、前記大気の状態を表す情報を利用して算出される光路輝度と前記第１中間光路輝度との乖離度を表すコストが最小となる場合の光路輝度を第２中間光路輝度として算出する請求項１に記載の情報処理装置。
【請求項３】
前記第２中間光路輝度算出手段は、前記観測画像が得られた場合の観測環境を利用したシミュレーションにより得られる前記電磁波の伝達状況に基づいて前記第２中間光路輝度を算出する請求項１に記載の情報処理装置。
【請求項４】
前記最終光路輝度算出手段は、前記観測画像と、前記観測画像が得られた場合の観測環境とのうちの少なくとも一方に基づいて前記暗画素の仮定に対する信頼度を算出し、算出した信頼度を重み係数として前記最終光路輝度を算出する請求項１又は請求項２又は請求項３に記載の情報処理装置。
【請求項５】
前記観測画像は地表の画像であり、
前記最終光路輝度算出手段は、前記観測画像の各画素に対応する地表領域が予め設定された特定の地表物を含む領域であるか否かを判断し、前記特定の地表物を含む領域であると判断された画素の数を利用して前記重み係数を算出し、当該算出した重み係数に基づいて前記最終光路輝度を算出する請求項４に記載の情報処理装置。
【請求項６】
請求項１乃至請求項５の何れか一つに記載の情報処理装置と、
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を前記情報処理装置に供給する画像供給装置と
を備える情報処理システム。
【請求項７】
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出し、
前記大気の状態を表す情報を含む第２中間光路輝度を算出し、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出し、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する
情報処理方法。
【請求項８】
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する処理と、
前記大気の状態を表す情報を含む第２中間光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する処理と
をコンピュータに実行させる処理手順が示されているコンピュータプログラムを記憶するプログラム記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、観測画像に含まれる情報を補正する技術に関する。
【背景技術】
【０００２】
人工衛星や航空機などの高所から観測装置によって地表を観測する技術は、リモートセンシングと呼ばれる。リモートセンシングでは、地表の所定の範囲の領域から放射される光などの電磁波の強さが観測されることが多い。リモートセンシングによって得られる観測結果は、画像に関連付けられた画素値として表されることが多い。つまり、画素値は、画像において、観測された領域の地表における位置に対応する画素に関連付けられるデータである。例えば、観測装置がイメージセンサである場合、観測結果は画像として得られる。その画像に含まれる画素値は、イメージセンサの受光素子に入射した光（観測光）の強さに応じた値である。【０００３】
なお、画素値が、観測された光の明るさを表す値である場合、画素値に代えて、その観測光の明るさを表す値を輝度値とも記載する。また、観測は、特定の範囲の波長帯域に含まれる波長の光が選択的に透過するフィルタを使用して行われることが多い。透過光の波長帯域が異なる複数のフィルタを使用することにより、複数の波長帯域毎の観測光の強さが観測結果として得られる。【０００４】
さらに、物体は、その表面の材質や状態によって、波長毎に異なる強度の光を反射することが知られている。物体における波長毎の光の反射率は、表面反射率と呼ばれる。当該物体の表面反射率は、リモートセンシングによって得られる画像の各画素に関連付けられる輝度値に基づいて算出することが可能である。この物体の表面反射率の情報に基づき、物体の状態や材質を取得するアプリケーション（応用技術）が期待されている。このようなアプリケーションが適用される分野として、営農支援や資源探査がある。【０００５】
営農支援や資源探査では、画像の各画素に関連付けられている波長帯域毎の輝度値（換言すれば、物体の表面反射率）を用いて、その農作物の生育具合の取得や鉱物の判別を行うことが期待される。的確な営農支援や資源探査を実現するためには、農作物や鉱石などの地表物の正確な情報を得る必要がある。【０００６】
なお、特許文献１には、スペクトルセンサにより計測された放射輝度を補正することにより、大気の影響を除く技術が開示されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１５−３２２０５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
リモートセンシングによって画像として得られる観測値（輝度値）は、観測対象による表面反射だけでなく、太陽光の明るさ、大気による光の吸収、そして大気による光の散乱に影響される。このことにより、観測値（輝度値）は、観測対象物体が持つ表面反射率に加えて、太陽光の明るさや大気の透過率に依存する成分と、大気で散乱されセンサに入射する光の成分を含む式（１）によって表される。
2017169946.tif 000003
【０００９】
式（１）におけるＬ（λ）は波長λにおける観測光の輝度値を表す。Ｒ（λ）は波長λにおける地表物の表面反射率を表す。α（λ）は波長λにおける太陽光の明るさと大気透過率に関連する成分（係数成分）を表す。β（λ）は、太陽光が大気で散乱され、地表物を経由せずにセンサに到達した波長λにおける光（散乱光）の輝度に関連する成分（加算成分）を表す。なお、太陽光が大気で散乱され、地表物を経由せずにセンサに到達した光（散乱光）の輝度を光路輝度とも記載する。【００１０】
式（１）によって表される観測光の輝度値Ｌ（λ）と地表物の表面反射率Ｒ（λ）との関係は、任意の波長帯域で観測された観測値に対して成り立つ。式（１）に表される係数成分α（λ）および加算成分β（λ）は、いわゆる環境雑音である。図１０は、高所撮影による観測画像における、観測光の輝度値Ｌ（λ）と、地表物の表面反射率Ｒ（λ）および環境雑音との関係をイメージで表す図である。【００１１】
ところで、環境雑音（つまり、太陽光の明るさ、大気の透過率、大気による散乱光の明るさ（光路輝度））は、日周運動による太陽高度の変化や大気の揺らぎ等の環境変化によって変動する。このため、地表物の正確な情報を得るためには、観測時における環境雑音を推定し、環境変化に依存する環境雑音の成分を補正する必要がある。つまり、加算成分β（λ）を推定し、この推定した加算成分β（λ）を式（１）から差し引くことによって、式（２）に表される補正された光路輝度Ｌ’（λ）が得られる。
2017169946.tif 000004
【００１２】
しかしながら、光路は複雑であるため、正確な光路輝度を得ることが難しいという問題が有る。【００１３】
本発明は上記課題を解決するために成されたものである。すなわち、本発明の主な目的は、光路輝度の補正精度を高める技術を提供することにある。
【課題を解決するための手段】
【００１４】
上記目的を達成するために、本発明の情報処理装置は、
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する第１中間光路輝度算出部と、
前記大気の状態を表す情報を含む第２中間光路輝度を算出する第２中間光路輝度算出部と、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する最終光路輝度算出部と、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する光路輝度補正部と、
を備える。【００１５】
また、本発明の情報処理システムは、
本発明の情報処理装置と、
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を前記情報処理装置に供給する画像供給装置と
を備える。【００１６】
さらに、本発明の情報処理方法は、
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出し、
前記大気の状態を表す情報を含む第２中間光路輝度を算出し、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出し、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する。【００１７】
さらにまた、本発明のプログラム記憶媒体は、
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する処理と、
前記大気の状態を表す情報を含む第２中間光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する処理と
をコンピュータに実行させるコンピュータプログラムを記憶する。【００１８】
なお、本発明の前記主な目的は、本発明の情報処理装置に対応する本発明の情報処理方法によっても達成される。また、本発明の前記主な目的は、本発明の情報処理装置、本発明の情報処理方法に対応する本発明のコンピュータプログラムおよびそれを格納する記憶媒体によっても達成される。
【発明の効果】
【００１９】
本発明によれば、光路輝度の補正精度を高めることができる。
【図面の簡単な説明】
【００２０】
【図１】本発明に係る第１実施形態の情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【図２】第１実施形態の情報処理装置の動作例を表すフローチャートである。
【図３】情報処理装置を実現するハードウェアの一構成例を説明する図である。
【図４】提案の情報処理装置および情報処理システムの構成を説明するブロック図である。
【図５】本発明に係る第２実施形態の情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【図６】本発明に係る第３実施形態の情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【図７】本発明に係る第４実施形態の情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【図８】本発明に係るその他の実施形態の情報処理装置の構成を簡略化して表すブロック図である。
【図９】図８の情報処理装置を用いる情報処理ステムの一例を説明する図である。
【図１０】観測光の輝度値と、地表物の表面反射率および環境雑音との関係をイメージで表す図である。
【発明を実施するための形態】
【００２１】
ここで、まず、本発明者が提案している情報処理装置の一構成例を説明する。
【００２２】
本発明者による提案の情報処理装置は、光路輝度の情報を含む観測画像に基づき光路輝度を算出し、この算出された光路輝度を補正する機能を備えている。すなわち、提案の情報処理装置は、波長帯域γにおける観測画像について、当該観測画像を構成する全ての画素にそれぞれ関連付けられている輝度値のうちの最も値が小さい最小輝度値を中間光路輝度βＷ（γ）と検知する。そして、提案の情報処理装置は、中間光路輝度βＷ（γ）を用いて大気状態を表す大気状態パラメタの真値Ｘｅを算出する。さらに、提案の情報処理装置は、大気状態パラメタＸを考慮した光路輝度βＭ（γ，Ｘ）を算出する数式におけるパラメタＸに、算出した大気状態パラメタの真値Ｘｅを代入することによって、最終的な光路輝度βＭ（γ，Ｘｅ）を算出する。
【００２３】
以降の説明において、大気状態パラメタＸを考慮して波長帯域γの光路輝度を算出する数式を光路輝度モデルと称し、βＭ（γ，Ｘ）と表す。
【００２４】
図４は、上述した提案の情報処理装置の一構成例を表すブロック図である。図４における情報処理装置５は、画像読み込み部５１と、暗画素推定部５２と、光路輝度推定部５３と、光路輝度補正部５４とを備える。画像読み込み部５１は画像供給装置６から観測画像を読み込む機能を備えている。暗画素特定部５２は、波長帯域毎に、読み込んだ観測画像から最小輝度値を検知する機能を備えている。ここで、波長帯域γにおける最小輝度値を中間光路輝度の算出値βＷ（γ）とする。
【００２５】
光路輝度推定部５３は、暗画素推定部５２により算出された中間光路輝度βＷ（γ）と、大気状態パラメタＸの値の変化に応じて値が変化する光路輝度モデルβＭ（γ，Ｘ）とを用いて、大気状態パラメタの真値Ｘｅを算出する機能を備えている。つまり、光路輝度推定部５３は、大気状態パラメタＸの値を変化させながら、各大気状態パラメタＸの値における光路輝度モデルβＭ（γ，Ｘ）および中間光路輝度βＷ（γ）に基づいてコストを算出する。ここでのコストとは、光路輝度モデルβＭ（γ，Ｘ）と中間光路輝度βＷ（γ）の乖離度を表す値である。
【００２６】
そして、光路輝度推定部５３は、コストが最も小さくなる場合の大気状態パラメタＸの値を大気状態パラメタの真値Ｘｅとする。また、光路輝度推定部５３は、光路輝度モデルβＭ（γ，Ｘ）の大気状態パラメタＸに算出値（真値）Ｘｅを代入して得られる値βＭ（γ，Ｘｅ）を最終光路輝度βｅ（γ）とする。
【００２７】
光路輝度補正部５４は、光路輝度推定部５３により算出された最終光路輝度βｅ（γ）を、観測画像の全ての画素における輝度値Ｌ（γ）から減算することにより、光路輝度を補正し、補正後の情報を出力装置７に出力する。
【００２８】
情報処理装置５は、上記のような構成を備えることにより、観測画像に含まれている情報の光路輝度を補正し、補正後の情報を出力できる。しかしながら、光路輝度モデルβＭ（γ，Ｘ）が内包する誤差の影響により、どのようなパラメタＸの値を選んでも光路輝度モデルβＭ（γ，Ｘ）と真の光路輝度との差が許容値を超えて大きくなる場合がある。この場合には、最終光路輝度βｅ（γ）が真の光路輝度と高精度に一致せず、情報処理装置５は、光路輝度を高精度に補正できないという問題が発生する。
【００２９】
そこで、本発明者は、光路輝度を算出する光路輝度モデルが誤差を内包する場合にも、光路輝度を高精度に補正すべく本発明を考え出した。
【００３０】
以下に、本発明に係る実施形態を図面を参照しつつ説明する。
【００３１】
＜第１実施形態＞
図１は、本発明に係る第１実施形態の情報処理装置およびそれを備えた情報処理システムの構成を簡略化して表すブロック図である。第１実施形態における情報処理システム１００は、大別して、情報処理装置１と、画像供給装置４と、出力装置３とを備えている。情報処理装置１は、画像供給装置４および出力装置３と通信可能に接続されている。
【００３２】
画像供給装置４は、例えば、互いに異なる複数の波長帯域について、観測対象にて反射した電磁波の強さを観測し、観測した結果を観測画像として出力する撮影装置である。あるいは、画像供給装置４は、そのような撮影装置によって観測された結果である観測画像を記憶するハードディスクなどの記憶装置やサーバ装置等であってもよい。
【００３３】
なお、以降の説明において、Ｎ（Ｎは２以上の整数）は、画像供給装置４が観測対象を観測する電磁波の波長帯域の数である。
【００３４】
画像供給装置４は、波長帯域に関連付けられたＮ個の観測画像を、情報処理装置１に供給する機能を備える。画像供給装置４は、Ｎ個の観測画像だけでなく、各観測画像に関連付けられた波長帯域を表す中心波長もしくは波長帯域の上限値および下限値をも情報処理装置１に供給してもよい。
【００３５】
例えば、観測対象が地表である場合には、画像供給装置４は、飛行機や人工衛星に搭載され、互いに異なる複数の波長帯域において、上空から、地表において反射した反射光を観測光の一部として観測する。具体例を挙げると、画像供給装置４は、カメラであり、互いに異なるＮ個の波長帯域のいずれかの光を選択的に透過する帯域通過フィルタを通ったＮ個の波長帯域の各観測光を撮影する。そして、画像供給装置４は、Ｎ個の撮影画像（観測結果）をそれぞれ観測画像として出力する。なお、画像供給装置４は、上空から地表を観測するのではなく、地表あるいは地表の近くから遠方の地表を観測してもよい。また、Ｎ個の各波長帯域の幅は、均一でなくてもよい。
【００３６】
画像供給装置４から出力されるＮ個の観測画像は、互いに異なるＮ個の波長帯域のいずれかについて観測した観測対象の明るさ分布を表す画像である。それぞれの観測画像において、観測画像を構成する各画素の輝度値は、当該画素に対応する方向から届いた観測光の強さを表す。なお、第１実施形態では、Ｎ個の観測画像は、例えば、同一の観測対象を観測することによって得られた画像であるが、異なる観測対象を観測して得られた画像であってもよい。ただ、この場合には、光路輝度の値が各画像の全観測領域について同一あるいはほぼ同一の場合である。異なる観測対象を観測して得られたＮ個の観測画像を利用する場合には、例えばシステム管理者がＮ個の画像を選択する。
【００３７】
出力装置３は、例えば、ディスプレイ装置である。また、出力装置３は、例えば、観測画像から撮影物体の材質や状態を示す情報を抽出する物体同定装置であってもよい。
【００３８】
情報処理装置１は、記憶部１１と、第１中間光路輝度算出部１２と、第２中間光路輝度算出部１３と、光路輝度補正部１４と、最終光路輝度算出部１５とを備える。
【００３９】
記憶部１１は、磁気ディスクあるいは半導体メモリ等の記憶デバイスである。記憶部１１は、画像供給装置４から供給されたＮ個の観測画像を、当該画像が観測された波長帯域の情報に関連付けて記憶する。観測画像に関連付けられる波長帯域の情報は、例えば、当該波長帯域を示す中心波長あるいは波長帯域の上限値および下限値の組み合わせでもよいし、各波長帯域に付与された識別子であってもよい。
【００４０】
第１中間光路輝度算出部１２と第２中間光路輝度算出部１３と光路輝度補正部１４と最終光路輝度算出部１５は、電子回路の場合もあれば、コンピュータプログラムとそのコンピュータプログラムに従って動作するプロセッサによって実現される場合もある。
【００４１】
第１中間光路輝度算出部１２は、記憶部１１に記憶されたＮ個の波長帯域のそれぞれに関連付けられた観測画像を構成する全画素の中から、輝度値が基準を満たす画素の輝度値を第１中間光路輝度として検知する機能を備える。例えば、第１中間光路輝度算出部１２は、基準を満たす輝度値として、全画素の輝度値の中から最小の輝度値を検知する。また、第１中間光路輝度算出部１２は、全画素の輝度値を小さい順に並べたときに、全画素数に所定の割合を掛け合わせた乗算値である順位の輝度値を、基準を満たす輝度値として検知してもよい。
【００４２】
以降の説明において、第１中間光路輝度算出部１２が検知する輝度値に対応する画素を暗画素と称する。暗画素は、観測される地表の暗い領域について表面反射率が零とみなせると仮定される画素である。また、Ｎ個の波長帯域におけるｉ（ｉは１乃至Ｎのいずれかの整数）番目の波長帯域をγｉと表す。さらに、Ｎ個の波長帯域γｉにおいて検知された暗画素の輝度値を第１中間光路輝度βＤ（γｉ）と表す。
【００４３】
第２中間光路輝度算出部１３は、第１中間光路輝度算出部１２によって算出された第１中間光路輝度βＤ（γｉ）を用いて、第２中間光路輝度βＥ（γｉ）を算出する機能を備える。
【００４４】
以降の説明において、第２中間光路輝度算出部１３が第２中間光路輝度βＥ（γｉ）を算出する過程で用いる光路輝度モデルβＭ（γｉ，Ｘ）に用いられるＸは、大気の状態を表すＱ個の大気状態パラメタＸｊ（ｊは１乃至Ｑのいずれかの整数）の集合を表す。大気状態パラメタＸｊは、大気に含まれる粒子（分子、エアロゾル、砂および煤等を含む粒子）に関する情報を表す。つまり、大気状態パラメタＸｊは、例えば、オングストローム指数、大気混濁度、大気透過率、大気分子の光学的厚さ、エアロゾルの光学的厚さ、視程、もしくはこれらを用いて算出される値である。
【００４５】
光路輝度モデルβＭ（γｉ，Ｘ）の具体例を、式（３）〜式（５）に表す。式（３）〜式（５）では、Ｎ個の波長帯域γｉにおける太陽光による放射照度をＩ（γｉ）と表す。
2017169946.tif 000005
【００４６】
2017169946.tif 000006
【００４７】
2017169946.tif 000007
【００４８】
例えば、第２中間光路輝度算出部１３は、Ｎ個の波長帯域γｉにおける太陽光による標準的な放射照度を放射照度Ｉ（γｉ）として用いる。また、第２中間光路輝度算出部１３は、観測環境を表す情報を外部の装置から取得し、下記の参考文献１に記載されている放射照度を表す物理モデルを用いて放射照度Ｉ（γｉ）を算出してもよい。観測環境を表す情報とは、例えば、観測時の太陽高度、観測年月日、観測時刻、観測対象領域の緯度および経度、観測対象領域の地形、観測時の水蒸気量およびエアロゾル量等を含む情報である。
［参考文献１］：R. E. Bird and C. Riordan， ”Simple Solar Spectral Model for Direct and Diffuse Irradiance on Horizontal and Tilted Planes at the Earth’s Surface for Cloudless Atmospheres”， Journal of Climatology and Applied Meteorology， Vol．25， pp．87-97
【００４９】
第２中間光路輝度算出部１３は、上述したような物理モデルを使用することにより、場所と時刻とから算出される太陽天頂角、および、大気の状態を表すパラメタを用いて直射日光と大気散乱光をシミュレーションにより算出する。これにより、第２中間光路輝度算出部１３は、例えば、晴天時の太陽光の分光放射照度を波長毎に算出できる。この物理モデルを用いる場合には、まず、シミュレーション結果に及ぼす影響の少ないパラメタの値が設定される。そして、第２中間光路輝度算出部１３は、設定されたパラメタの値と、観測環境を表す情報とのうち、太陽天頂角および大気の状態を表す値を用いて、晴天時における太陽光の単位波長あたりの放射照度である分光放射照度を算出する。そして、第２中間光路輝度算出部１３は、晴天時における太陽光の分光放射照度を、波長帯域γｉにおける複数の波長について足し合わせる。これにより、第２中間光路輝度算出部１３は、Ｎ個の波長帯域γｉ毎に放射照度Ｉ（γｉ）を算出する。
【００５０】
第２中間光路輝度算出部１３は、第１中間光路輝度算出部１２によって算出された第１中間光路輝度βＤ（γｉ）を用いて、例えば次のようにして第２中間光路輝度βＥ（γｉ）を算出する。すなわち、第２中間光路輝度算出部１３は、まず、第１中間光路輝度βＤ（γｉ）を用いて、光路輝度モデルβＭ（γｉ，Ｘ）の大気状態パラメタの真値Ｘｅを算出する。具体的には、例えば、第２中間光路輝度算出部１３は、大気状態パラメタＸの値を変えながら、各Ｘの値における光路輝度モデルβＭ（γｉ，Ｘ）と第１中間光路輝度βＤ（γｉ）を用いてコストＣを算出する。コストＣは、第１中間光路輝度βＤ（γｉ）と、観測時における実際の光路輝度との乖離度を表す値である。式（６）は、コストＣを算出する数式の一例である。
2017169946.tif 000008
【００５１】
第２中間光路輝度算出部１３は、コストＣが最も小さくなるときの大気状態パラメタＸの値を大気状態パラメタの真値Ｘｅとして算出する。
【００５２】
大気状態パラメタの真値Ｘｅを算出した後に、第２中間光路輝度算出部１３は、光路輝度モデルβＭ（γｉ，Ｘ）における大気状態パラメタＸに算出値（真値）Ｘｅを代入することによって、第２中間光路輝度βＥ（γｉ）を算出する。
【００５３】
最終光路輝度算出部１５は、第１中間光路輝度算出部１２によって算出された第１中間光路輝度βＤ（γｉ）と、第２中間光路輝度算出部１３によって算出された第２中間光路輝度βＥ（γｉ）を用いて、最終光路輝度βＰ（γｉ）を算出する機能を備えている。
【００５４】
最終光路輝度算出部１５は、式（７）に表される数式を用いて、第１中間光路輝度βＤ（γｉ）と第２中間光路輝度βＥ（γｉ）とを重み付けして足し合わせ、その結果を最終光路輝度βＰ（γｉ）として算出する。式（７）において、ａ（γｉ）は、暗画素の仮定に対する信頼度を波長帯域γｉ毎に表す重み係数である。以降の説明においてａ（γｉ）を暗画素仮定信頼係数とも称す。
2017169946.tif 000009
【００５５】
例えば、暗画素仮定信頼係数ａ（γｉ）には、０乃至１の任意の値が設定される。例えば、画像供給装置４から、８個の波長帯域を観測して得られた観測画像が情報処理装置１に供給されるとする。この場合における波長帯域γ１，γ２，γ３，γ４，γ５，γ６，γ７，γ８の各中心波長が、４２５ｎｍ，４８０ｎｍ，５４５ｎｍ，６０５ｎｍ，７２５ｎｍ，８３２．５ｎｍ，９５０ｎｍであるとする。この場合、例えば、各暗画素仮定信頼係数ａ（γ１），ａ（γ２），ａ（γ３），ａ（γ４），ａ（γ５），ａ（γ６），ａ（γ７），ａ（γ８）には、１，１，０，０，０，０，０，１が設定される。最終光路輝度算出部１５は、そのように設定された重み係数と、第１中間光路輝度βＤ（γｉ）と、第２中間光路輝度βＥ（γｉ）と、式（７）とを利用して、最終光路輝度βＰ（γｉ）を算出する。
【００５６】
光路輝度補正部１４は、記憶部１１が記憶した観測画像、および、最終光路輝度算出部１５が算出した最終光路輝度βＰ（γｉ）を用いて、波長帯域毎に、観測画像の全画素にそれぞれ関連付けられている輝度値から最終光路輝度を減算する機能を備えている。光路輝度補正部１４は、その減算により補正された画像を光路輝度補正画像とし、出力装置３に出力する。すなわち、減算された観測画像は、記憶部１１が記憶する観測画像から、最終光路輝度算出部１５によって算出された光路輝度の成分が除かれた画像（光路輝度補正画像）である。
【００５７】
なお、第１実施形態では、第２中間光路輝度と最終光路輝度と光路輝度モデルは、観測画像と関連付けられているＮ個の波長帯域γｉが考慮されているが、任意の個数の、任意の範囲の波長帯域における光路輝度や光路輝度モデルであってもよい。この場合、予め設定された個数の設定された範囲の波長帯域における光路輝度もしくは光路輝度モデルを用いて、観測画像と関連付けられているＮ個の波長帯域γｉ毎の光路輝度が算出される。
【００５８】
次に、図２のフローチャートを参照して、第１実施形態における情報処理システム１００の情報処理装置１の動作について説明する。
【００５９】
まず、記憶部１１が、画像供給装置４から供給されたＮ個の観測画像を当該観測画像が観測された波長帯域の情報に関連付けて記憶する（ステップＳ１０１）。
【００６０】
第１中間光路輝度算出部１２は、波長帯域毎に関連付けられた観測画像について、全画素の中から暗画素（輝度値が基準を満たす画素）の輝度値を第１中間光路輝度として検知する（ステップＳ１０２）。
【００６１】
第２中間光路輝度算出部１３は、光路輝度モデルにおける大気状態パラメタの値Ｘを変化させながら当該光路輝度モデルに基づいた光路輝度を算出し、算出した各光路輝度と第１中間光路輝度との乖離度を表すコストＣを算出する（ステップＳ１０３）。
【００６２】
第２中間光路輝度算出部１３は、コストＣが最小となるときの大気状態パラメタを大気状態パラメタの真値Ｘｅとする（ステップＳ１０４）。
【００６３】
そして、第２中間光路輝度算出部１３は、光路輝度モデルの大気状態パラメタＸに大気状態パラメタの真値Ｘｅを代入することにより第２中間光路輝度を算出する（ステップＳ１０５）。
【００６４】
最終光路輝度算出部１５は、第１中間光路輝度算出部１２による第１中間光路輝度と、第２中間光路輝度算出部１３による第２中間光路輝度とを、暗画素仮定の信頼度に応じた重み係数を用いて重み付けして足し合わせる。これにより、最終光路輝度算出部１５は、最終光路輝度を算出する（ステップＳ１０６）。
【００６５】
光路輝度補正部１４は、記憶部１１が記憶した観測画像の全画素から最終光路輝度を減算することにより光路輝度補正画像を生成する。光路輝度補正部１４は、その光路輝度補正画像を出力装置３に出力する（ステップＳ１０７）。
【００６６】
第１実施形態における情報処理装置１および情報処理システム１００は、光路輝度算出に用いられる光路輝度モデルが誤差を内包する場合にも、光路輝度を高精度に算出することができる。その理由は、最終光路輝度算出部１５は、第１中間光路輝度算出部１２による第１中間光路輝度と、第２中間光路輝度算出部１３による第２中間光路輝度とを、暗画素仮定の信頼度を表す重み係数を用いて重み付けして足し合わせる。これによって、最終光路輝度算出部１５が最終光路輝度を算出するからである。
【００６７】
すなわち、図４に表される発明者提案の情報処理装置５は、前述したように、各観測画像の最小輝度値を中間光路輝度としている。そして、情報処理装置５は、その中間光路輝度と光路輝度モデルを用いて大気状態パラメタの真値を算出し、算出した真値を光路輝度モデルに代入することにより最終光路輝度を算出している。しかしながら、その光路輝度算出に用いられる光路輝度モデルが誤差を内包する場合には、どのような大気状態パラメタの値を代入しても光路輝度モデルが光路輝度を高精度に算出できない。この場合には、最終光路輝度は光路輝度と高精度に一致しないため、情報処理装置５は、光路輝度を高精度に算出できない。
【００６８】
これに対して、第１実施形態における情報処理装置１および情報処理システム１００は、暗画素仮定の信頼度が高い波長帯域においては、暗画素仮定に基づき算出された第１中間光路輝度に近い値を最終光路輝度に設定できる。その結果、光路輝度算出に用いられる光路輝度モデルが誤差を内包する場合にも、光路輝度モデルの内包する誤差の影響を緩和し、光路輝度を高精度に算出することができる。
【００６９】
なお、第１実施形態の情報処理装置１を構成する各機能部は、専用のハードウェア（電子回路）によって実現することができる。また、少なくとも、第１中間光路輝度算出部１２と第２中間光路輝度算出部１３と光路輝度補正部１４と最終光路輝度算出部１５は、ソフトウェアプログラムの機能（処理）単位（ソフトウェアモジュール）と捉えることもできる。ただし、これら機能部は、説明の便宜上の構成であり、実装に際しては、様々な構成が想定され得る。第１実施形態の情報処理装置１を実現するハードウェア環境の一例を図３を参照して説明する。
【００７０】
図３は、第１実施形態の情報処理装置１を実行可能な情報処理装置（コンピュータ）の構成を例示的に表す図である。図３における情報処理装置１は、構成要素として下記を備えている。
・ＣＰＵ（Central Processing Unit）９０１
・ＲＯＭ（Read Only Memory）９０２
・ＲＡＭ（Random Access Memory）９０３
・ハードディスク（記憶装置）９０４
・外部装置との通信インタフェース９０５
・ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）等の記録媒体９０７に格納されたデータを読み書き可能なリーダライタ９０８
・入出力インタフェース９０９
情報処理装置１は、これらの構成がバス（通信線）９０６を介して接続されたコンピュータにより実現可能である。
【００７１】
情報処理装置１のＣＰＵ９０１が、前述した各機能部１２〜１５を実現するコンピュータプログラムを実行することにより、各機能部１２〜１５が実現される。なお、そのコンピュータプログラムは読み書き可能な揮発性のメモリ（ＲＡＭ９０３）又はハードディスク９０４等の不揮発性の記憶デバイスに格納される。
【００７２】
また、情報処理装置１へのコンピュータプログラムの供給手法は、例えば、ＣＤ−ＲＯＭ等の各種記録媒体９０７を介して当該装置内にインストールする手法や、インターネット等の通信回線を介して外部よりダウンロードする手法等がある。また、コンピュータプログラムや当該コンピュータプログラムを構成するコードは、格納された記録媒体９０７によって構成されると捉えることもできる。
【００７３】
＜第２実施形態＞
以下に、本発明に係る第２実施形態を説明する。なお、第２実施形態の説明において、第１実施形態の情報処理装置および情報処理システムを構成する構成部分と同一名称部分には同一符号を付し、その共通部分の重複説明は省略する。
【００７４】
図５は、第２実施形態に係る情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【００７５】
第２実施形態では、画像供給装置４は、第１実施形態で説明した機能に加えて、観測したときの観測環境を表す情報（観測環境情報）を情報処理装置１に提供する機能を備えている。観測環境情報は、例えば、観測時の太陽高度、観測年月日、観測時刻、観測対象領域の緯度および経度、観測対象領域の地形、観測対象領域における水蒸気量およびエアロゾル量等を含む。
【００７６】
記憶部１１は、画像供給装置４から供給された観測環境情報をも観測画像に関連付けて記憶する。
【００７７】
第２実施形態の情報処理装置１における最終光路輝度算出部１５以外の構成は、第１実施形態の情報処理装置１と同様である。
【００７８】
最終光路輝度算出部１５は、第１中間光路輝度算出部１２による第１中間光路輝度βＤ（γｉ）と、第２中間光路輝度算出部１３による第２の中間光路輝度βＥ（γｉ）とを用いる。さらに、最終光路輝度算出部１５は、記憶部１１に記憶されている観測画像と観測環境情報とのうちの少なくとも一方とをも用いる。最終光路輝度算出部１５は、それら情報に基づいて最終光路輝度βＰ（γｉ）を算出する。例えば、最終光路輝度算出部１５は、式（７）に表される数式を用いて、第１中間光路輝度βＤ（γｉ）と第２中間光路輝度βＥ（γｉ）を重み付けして足し合わせ、その結果を最終光路輝度βＰ（γｉ）として算出する。式（７）において、第１実施形態でも述べたように、ａ（γｉ）は、暗画素仮定に対する信頼の度合を観測波長帯域γｉ毎に表す重み係数である。
【００７９】
具体的には、最終光路輝度算出部１５は、記憶部１１に記憶されている観測環境情報と観測画像のうちの少なくとも一方を用いて暗画素仮定信頼係数を算出する。例えば、最終光路輝度算出部１５は、波長帯域γｉ毎に関連付けられた観測画像の各画素に対応する地表領域が、一つもしくは複数の特定地表物を含むか否かを判定する。そして、最終光路輝度算出部１５は、特定地表物が含まれる画素の数を算出し、算出した画素数、もしくは、画像の全画素数からその算出した画素数を差し引いた数に基づいて暗画素仮定信頼係数ａ（γｉ）を設定（算出）する。特定地表物は、例えば、水、人工物、植生、土壌、雪である。
【００８０】
より具体的には、最終光路輝度算出部１５は、Ｎ個の波長帯域に関連付けられた各観測画像におけるｊ番目の画素から得られる輝度値Ｌｊ（γ１）， …， Ｌｊ（γＮ）のうち、赤に相当する波長帯域に関連付けられた輝度値Ｌｊ（γＲ）と、近赤外に相当する波長帯域に関連付けられた輝度値Ｌｊ（γＮＩＲ）とを選択する。そして、最終光路輝度算出部１５は、それら選択した輝度値Ｌｊ（γＲ），Ｌｊ（γＮＩＲ）と、式（８）に示す数式とを用いて、正規化差分植生指数Ｉｊを算出する。
2017169946.tif 000010
【００８１】
さらに、最終光路輝度算出部１５は、正規化差分植生指標Ｉｊの値が閾値よりも大きい場合にｊ番目の画素を植生領域と判定する。最終光路輝度算出部１５は、そのような処理を、観測画像を構成する全ての画素に対し実行することによって、植生を含む地表領域の画素数を算出し、当該算出した画素数に基づいて暗画素仮定信頼係数を設定（算出）する。
【００８２】
なお、最終光路輝度算出部１５は、正規化差分植生指標の代わりに、下記の参考文献２に記載されている植生指標（例えば、ＳＡＶＩ（Soil Adjusted Vegetation Index）、ＲＶＩ（Rate Vegetation Index）、ＮＲＶＩ（Normalized Ratio Vegetation Index）、ＴＶＩ（Transformed Vegetation Index）、ＣＴＶＩ（Corrected Transformed Vegetation Index）、ＴＴＶＩ（Thiam’s Transformed Vegetation Index）、ＥＶＩ（Enhanced Vegetation Index）など）を用いてもよい。また、最終光路輝度算出部１５は、正規化差分植生指標の代わりに、参考文献３に記載されている水指数（例えば、ＮＤＷＩ（Normalized Difference Water Index））や土壌指数（例えば、ＮＤＳＩ（Normalized Difference Soil Index））やＮＨＦＤ（Non−Homogeneous Feature Difference）等を用いてもよい。
［参考文献２］：Nikolaos G． Silleos， Thomas K． Alexandridis， Ioannis Z． Gitas and Konstantinos Perakis， “Vegetation Indices： Advances Made in Biomass Estimation and Vegetation Monitoring in the Last 30 Years”， Geocarto International， Volume 21， Issue 4， 2006
［参考文献３］：A． F． Wolf， “Using WorldView−2 Vis−NIR Multispectral Imagery to Support Land Mapping and Feature Extraction Using Normalized Difference Index Ratios”， May 8， 2012
【００８３】
例えば、最終光路輝度算出部１５は、初期値として、全観測波長帯域に関連付けられる暗画素仮定信頼係数ａ（γｉ）（ｉ＝１，・・・，Ｎ）を１に設定する。そして、最終光路輝度算出部１５は、画像の全画素数から植生領域の画素数を差し引いた数が予め設定された値よりも小さい場合には、近赤外に相当する波長帯域γＮＩＲの暗画素仮定信頼係数ａ（γＮＩＲ）を０に設定する。また、水指数に基づいて水領域の画素が検知される場合であって、かつ、画像の全画素数から水領域の画素の数を差し引いた数が予め設定された値よりも小さい場合がある。この場合には、最終光路輝度算出部１５は、近紫外に相当する波長帯域γＮＵＶの暗画素仮定信頼係数ａ（γＮＵＶ）を０に設定する。さらに、最終光路輝度算出部１５は、人工物領域の画素が検知される場合であって、かつ、人工物領域の画素の数が予め設定された値よりも大きい場合には、全観測波長帯域に関連付けられる暗画素仮定信頼係数を１に設定する。
【００８４】
なお、最終光路輝度算出部１５は、記憶部１１に記憶されている観測対象の緯度・経度・観測日時の情報を用いて、観測画像の画素のうち陰の領域が含まれる画素の数を算出する。そして、最終光路輝度算出部１５は、この算出した画素数を用いて、全観測波長帯域に関連付けられる暗画素仮定信頼係数を設定してもよい。例えば、緯度・経度・観測日時と陰量を対応付けたルックアップテーブルが予め与えられている。最終光路輝度算出部１５は、記憶部１１に記憶されている観測対象の緯度・経度・観測日時およびルックアップテーブルに基づいて、陰の領域が含まれる画素の数を算出する。そして、最終光路輝度算出部１５は、その算出した画素の数に基づいて暗画素仮定信頼係数を設定する。また、最終光路輝度算出部１５は、記憶部１１に記憶されている観測環境情報および観測画像の両方を用いて、暗画素仮定信頼係数を設定してもよい。例えば、最終光路輝度算出部１５は、観測画像を用いて特定地表物が含まれる画素の数を算出し、また、観測環境情報を用いて陰の領域が含まれる画素の数を算出する。そして、最終光路輝度算出部１５は、特定地表物が含まれる画素の数および陰の領域が含まれる画素の数の両方を用いて、暗画素仮定信頼係数を設定する。
【００８５】
最終光路輝度算出部１５は、上記のように設定された暗画素仮定信頼係数と、第１中間光路輝度βＤ（γｉ）と、第２中間光路輝度βＥ（γｉ）と、式（７）とを利用して、最終光路輝度βＰ（γｉ）を算出する。
【００８６】
この第２実施形態の情報処理装置１および情報処理システム１００は、暗画素仮定信頼係数を用いて最終光路輝度を算出することにより、光路輝度モデルが誤差を内包する場合にも光路輝度モデルの内包する誤差の影響を緩和し、光路輝度を高精度に算出できる。
【００８７】
＜第３実施形態＞
以下に、本発明に係る第３実施形態を説明する。なお、第３実施形態の説明において、第１と第２実施形態の情報処理装置および情報処理システムを構成する構成部分と同一名称部分には同一符号を付し、その共通部分の重複説明は省略する。
【００８８】
図６は第３実施形態に係る情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【００８９】
この第３実施形態においても、第２実施形態と同様に、画像供給装置４は、情報処理装置１に向けて観測環境情報を提供し、記憶部１１は、その観測環境情報を観測画像に関連付けて記憶する。
【００９０】
第３実施形態の情報処理装置１における第２中間光路輝度算出部１３以外の構成は、第１実施形態の情報処理装置１と同様である。
【００９１】
第３実施形態では、第２中間光路輝度算出部１３は、観測環境情報を記憶部１１から取得し、第１中間光路輝度に代えて、取得した観測環境情報を用いて第２中間光路輝度をシミュレーションにより算出する。例えば、第２中間光路輝度算出部１３は、参考文献４に記載されているシミュレータを利用する。
［参考文献４］：Daniel Schlapfer， Dr． sc． Nat．， ReSe． “MODO User Manual Ver.5”， 2011
【００９２】
第３実施形態の情報処理装置１および情報処理システム１００は、観測環境情報を考慮して第２中間光路輝度を算出することにより、光路輝度を高精度に算出できる。
【００９３】
＜第４実施形態＞
以下に、本発明に係る第４実施形態を説明する。なお、第４実施形態の説明において、第１や第２の実施形態の情報処理装置および情報処理システムを構成する構成部分と同一名称部分には同一符号を付し、その共通部分の重複説明は省略する。
【００９４】
図７は第４実施形態に係る情報処理装置および情報処理システムの構成を簡略化して表すブロック図である。
【００９５】
第４実施形態の情報処理装置１における第２中間光路輝度算出部１３と最終光路輝度算出部１５以外の構成は、第１あるいは第２の実施形態の情報処理装置１と同様である。
【００９６】
第４実施形態では、第２中間光路輝度算出部１３は、観測環境情報を記憶部１１から取得し、第１中間光路輝度に代えて、取得した観測環境情報を用いて第２中間光路輝度を算出する。また、最終光路輝度算出部１５は、記憶部１１から取得した観測環境情報をも考慮して、最終光路輝度を算出する。
【００９７】
この第４実施形態の情報処理装置１および情報処理システム１００は、観測環境情報を考慮して最終光路輝度を算出するので、光路輝度をより高精度に算出できる。
【００９８】
＜その他の実施形態＞
なお、本発明は第１〜第４の実施形態に限定されず、様々な実施の形態を採り得る。例えば、本発明に係る情報処理装置のその他の形態として、図８に表されるような構成をも採り得る。図８に表される情報処理装置２０は、図９に表されるような画像供給装置３１と接続し、情報処理システム３０を構成する。画像供給装置３１は、大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる波長帯域毎の観測画像を情報処理装置２０に供給する機能を備えている。
【００９９】
図８に表される情報処理装置２０は、第１中間光路輝度算出部２１と、第２中間光路輝度算出部２２と、最終光路輝度算出部２３と、光路輝度補正部２４とを備えている。
【０１００】
第１中間光路輝度算出部２１は、画像供給装置３１から供給される波長帯域毎の観測画像を構成する各画素に関連付けられている電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出す機能を備えている。さらに、第１中間光路輝度算出部２１は、その輝度値に基づいた第１中間光路輝度を算出する機能を備えている。
【０１０１】
第２中間光路輝度算出部２２は、大気の状態を表す情報を含む第２中間光路輝度を算出する機能を備えている。
【０１０２】
最終光路輝度算出部２３は、波長帯域毎に、観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、第１中間光路輝度および第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する機能を備えている。
【０１０３】
光路輝度補正部２４は、波長帯域毎に、観測画像の全画素から最終光路輝度を減算する機能を備えている。
【０１０４】
この情報処理装置２０およびそれを備えた情報処理システム３０は、上記のように、第１中間光路輝度と第２中間光路輝度に、暗画素の仮定に対する信頼度を表す重み係数を用いて重み付けする。そして、情報処理装置２０および情報処理システム３０は、その重み付けされた第１中間光路輝度と第２中間光路輝度を利用して算出された最終光路輝度に基づいて、光路輝度を補正する。このため、情報処理装置２０および情報処理システム３０は、暗画素の仮定に対する信頼度が低い第１や第２の中間光路輝度が光路輝度の補正に関与する度合いを低くすることができる。これにおり、情報処理装置２０および情報処理システム３０は、光路輝度補正の精度を高めることができる。
【０１０５】
以上、上述した実施形態を模範的な例として本発明を説明した。しかしながら、本発明は、上述した実施形態には限定されない。即ち、本発明は、本発明のスコープ内において、当業者が理解し得る様々な態様を適用することができる。
【０１０６】
この出願は、２０１６年３月３０日に出願された日本出願特願２０１６−０６８５４１を基礎とする優先権を主張し、その開示の全てをここに取り込む。
【符号の説明】
【０１０７】
１  情報処理装置
３  出力装置
４  画像供給装置
１２  第１中間光路輝度算出部
１３  第２中間光路輝度算出部
１４  光路輝度補正部
１５  最終光路輝度算出部
【図１】
2017169946.tif 000011
【図２】
2017169946.tif 000012
【図３】
2017169946.tif 000013
【図４】
2017169946.tif 000014
【図５】
2017169946.tif 000015
【図６】
2017169946.tif 000016
【図７】
2017169946.tif 000017
【図８】
2017169946.tif 000018
【図９】
2017169946.tif 000019
【図１０】
2017169946.tif 000020
【手続補正書】
【提出日】平成30年9月10日(2018.9.10)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項８
【補正方法】変更
【補正の内容】
【請求項８】
大気を通った互いに異なる複数の波長帯域の電磁波を観測することにより得られる前記波長帯域毎の観測画像を構成する各画素に関連付けられている前記電磁波の観測値である輝度値のうち、予め設定された基準を満たす輝度値を取り出し、当該輝度値に基づいた第１中間光路輝度を算出する処理と、
前記大気の状態を表す情報を含む第２中間光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の画素のうち、物体により反射された電磁波の影響を受けていないと仮定される画素である暗画素の仮定に対する信頼度を表す重み係数を用いて、前記第１中間光路輝度および前記第２中間光路輝度を重み付けして足し合わせることにより最終光路輝度を算出する処理と、
前記波長帯域毎に、前記観測画像の全画素から前記最終光路輝度を減算する処理と
をコンピュータに実行させるコンピュータプログラム。
【国際調査報告】
2017169946.tif 000021
2017169946.tif 000022