JP2000131218A - 水中物質の濃度推定方法及び装置、記録媒体 - Google Patents
水中物質の濃度推定方法及び装置、記録媒体Info
- Publication number
- JP2000131218A JP2000131218A JP10304040A JP30404098A JP2000131218A JP 2000131218 A JP2000131218 A JP 2000131218A JP 10304040 A JP10304040 A JP 10304040A JP 30404098 A JP30404098 A JP 30404098A JP 2000131218 A JP2000131218 A JP 2000131218A
- Authority
- JP
- Japan
- Prior art keywords
- concentration
- light
- water
- radiation amount
- substance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 165
- 239000000126 substance Substances 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000003287 optical effect Effects 0.000 claims abstract description 30
- 238000010521 absorption reaction Methods 0.000 claims abstract description 17
- 230000005855 radiation Effects 0.000 claims description 113
- 238000004364 calculation method Methods 0.000 claims description 32
- 238000012545 processing Methods 0.000 claims description 19
- 238000012546 transfer Methods 0.000 claims description 16
- 230000007613 environmental effect Effects 0.000 claims description 15
- 238000011156 evaluation Methods 0.000 claims description 12
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 description 28
- 238000006243 chemical reaction Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 9
- 239000000243 solution Substances 0.000 description 8
- 239000013535 sea water Substances 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 230000008033 biological extinction Effects 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000012086 standard solution Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
画像データを利用して水中物質の濃度を高精度で推定す
る装置を提供する。 【解決手段】 水面下向き放射量算出部16では、撮像
環境情報に基づいて特定水域の水面に入射する放射量を
算出する。水面上向き放射量算出部(1)15は、変換
された衛星画像データの放射輝度データに基づいて該特
定水域における光の観測放射量を算出する。水面上向き
放射量算出部(2)17は、当該放射輝度と水中物質の
光学特性情報とに基づいて該特定水域における光の理論
放射量を算出する。水中物質濃度推定部18は、光の理
論放射量における吸収及び散乱係数に含まれる水中物質
の濃度を、非線形最小2乗問題に即した濃度推定関数に
基づいて導出する。
Description
グを利用した画像解析技術に係り、より詳しくは、人工
衛星等のリモートセンシングによるマルチスペクトル画
像データから特定水域における水中物質(溶存物質、懸
濁物質等、以下同じ)の濃度を推定する手法に関する。
ここでマルチスペクトル画像データとは、光を用いて所
定の波長帯域毎に撮影した複数の画像データから成るも
のである。
ングを利用して、国土や海域等の環境を広域且つ周期的
に観測する技術が知られている。また、この技術の適用
例として、マルチスペクトル画像データを用いた、特定
水域(海域)における水中物質の濃度推定手法について
の検討がなされている。従来のこの種の濃度推定方法
は、例えば、対象となる特定水域の複数の観測地点にお
いて実測された水中物質の濃度データ(以下、トゥルー
スデータ:TruthData)と、リモートセンシングにより
撮影された当該特定水域のマルチスペクトル画像データ
に含まれる輝度値との相関関係を利用して、その画像全
体における水中物質の濃度を推定するものである。この
ような手法は、経験的な回帰関係式や統計的手法に基づ
いた、所謂、Empiricalな枠組みによるもので
ある。
な経験的な回帰関係式や統計的手法に即した従来の枠組
みによる濃度推定手法では、例えば、時空間的変動が大
きく且つ複数の汚染物質が混在するような沿岸域を濃度
推定の対象とする場合に、特定水域毎に解析精度がバラ
つくために、現地観測データである膨大なトゥルースデ
ータを再度取り直さなければならないという問題があっ
た。また、膨大なトゥルースデータを取り直した場合に
は、その都度、回帰関係式等の計算をやり直す必要があ
り、演算処理が煩雑になるばかりでなく、濃度推定に係
る演算サイクルが長くなり、トゥルースデータと濃度推
定値との濃度誤差を低減させることが難しいという問題
があった。
画像データを利用して特定水域における水中物質の濃度
を推定する場合に、当該画像データの撮像領域や、撮像
時における大気状態等の撮像環境が変化した場合であっ
ても濃度の推定精度を一定水準に維持することを可能と
なる、改良された濃度推定方法を提供することにある。
本発明の他の課題は、上記濃度推定方法の実施に適した
濃度推定装置を提供することにある。本発明の他の課題
は、上記濃度推定方法及び装置等をコンピュータ装置上
で実現するための記録媒体を提供することにある。
は、特定水域をその上空から複数の波長帯域で撮影して
得たマルチスペクトル画像データと撮影時の環境情報、
例えば画像の解像度、画像センサの検出波長、大気の状
態、太陽の位置等の情報とを取得して前記特定水域にお
ける水中物質の濃度を導出する装置における濃度推定方
法であって、前記マルチスペクトル画像データから前記
特定水域における第1光放射量を導出するとともに、前
記環境情報及び既定の水中物質群の光学特性情報から前
記特定水域における第2光放射量を導出し、導出した前
記第1光放射量および第2光放射量に基づいて光学特性
が表す水中物質の濃度を推定する処理を、前記装置上で
実行することを特徴とする。
程式を水中における光の多流束モデルを用いて近似化す
ることにより第1光放射量及び第2光放射量を要素に含
む濃度推定関数を作成し、この濃度推定関数に前記導出
した第1光放射量及び第2光放射量を適用することで、
前記水中物質の濃度を推定するものである。前記濃度推
定関数は、例えば第1および第2光放射量が前記水中物
質の濃度に非線形に依存する形態のものであり、前記光
学特性情報は、例えば前記水中物質の光に対する粒子レ
ベルのミクロな反射特性及び吸収特性を含む情報であ
る。
装置は、前記マルチスペクトル画像データ及び前記環境
情報を取得する取得手段と、取得したマルチスペクトル
画像データを放射輝度データに換算するとともに換算さ
れた放射輝度データに基づいて前記特定水域における第
1光放射量を導出する第1光放射量導出手段と、前記環
境情報及び既定の水中物質群の光学特性情報から前記特
定水域における第2光放射量を導出する第2光放射量導
出手段と、光の放射伝達方程式を水中における光の多流
束モデルを用いて近似化することにより前記観測放射量
及び前記理論放射量を要素に含む濃度推定関数を作成
し、前記導出した観測放射量及び理論放射量を前記作成
した濃度推定関数に適用して前記濃度を推定する濃度推
定手段と、を備えたことを特徴とする。
記環境情報から前記特定水域の水面に入射する光の下向
き放射量を算出するとともに、前記算出した下向き放射
量及び既定の水中物質群の光学特性情報に基づいて前記
特定水域の水面から放射される光の放射量を算出し、こ
れらの算出結果に基づいて前記第2光放射量を導出する
ように構成される。また、前記特定水域の水面から放射
される、前記波長帯域毎の最大波長及び最小波長の放射
輝度の積分値を算出し、この積分値に基づいて前記第2
光放射量を導出するように構成される。
第1光放射量及び第2光放射量が水中物質の濃度に対し
て非線形に依存する、所定の非線形最小2乗問題に基づ
いて前記濃度推定関数を作成し、前記第2光放射量にお
ける吸収係数及び散乱係数に含まれる水中物質の濃度要
素を前記濃度推定関数に適用して濃度を推定するととも
に、この推定結果に即して前記濃度要素を更新するよう
に構成される。この更新は、前記第2光放射量と前記第
1放射量との差分に基づく評価値が最小になるまで反復
される。これにより、最適な濃度が得られる。
は、前記マルチスペクトル画像データと前記環境情報と
を取得する機能を有するコンピュータ装置に読み取られ
たときに、前記コンピュータ装置に少なくとも下記の処
理を実行させるためのプログラムが記録された記録媒体
である。 (1)前記マルチスペクトル画像データから前記特定水
域における第1光放射量を導出するとともに、前記環境
情報及び既定の水中物質群の光学特性情報から前記特定
水域における第2光放射量を導出する処理、(2)光の
放射伝達方程式を水中における光の多流束モデルを用い
て近似化することにより前記観測放射量及び前記理論放
射量を要素に含む濃度推定関数を作成し、前記導出した
観測放射量及び前記理論放射量を前記作成した濃度推定
関数に適用して前記水中物質の濃度を推定する処理。
定水域のマルチスペクトル画像データ(以下、衛星画像
データ)を用いて、特定水域の水中物質の濃度推定を行
う場合の実施の形態を説明する。濃度推定は、例えばコ
ンピュータを用いた濃度推定装置によって行う。ここで
は、この濃度推定装置の主制御手段(CPU)が、以下
のようにして衛星画像データと光(電磁波を含む)の放
射伝達理論とに基づく水中物質の濃度を算出すること
で、濃度推定を行う場合の例を挙げる。
について説明する。なお、ここでは、以下、海域を対象
とした例を用いて説明する。光の放射伝達現象を表す放
射伝達方程式は、水深をz、光の伝達方向(天頂角)を
θ、光の伝達方向(方位角)をφ、光の波長をλ、消散
係数をc、入射する放射輝度をL、散乱により生じる放
射輝度をLE *、源泉関数をS、cosθをμとすると、
以下の式(1)のようになる。
算出することは非常に困難であり、解を得るための様々
な数値的解法が提案されている。そこで本実施形態で
は、まず、水中における光の多流束モデルに基づく解法
を用いて、上記式(1)の放射伝達過程を近似化する。
この場合の多流束モデル、例えば四流束モデルの概念図
を図6に示す。この図6に示すように、水深zにおいて
微少区間dzに下向きに入射する波長λの直達放射照度
をEcd(z,λ)、拡散放射照度をEd(z,λ)と
し、上向きに入射する直達放射照度をEcu(z,λ)、
拡散放射照度をEu(z,λ)とすれば、波長λの光の
放射伝達過程は、以下に示す式(2)〜(5)の連立微
分方程式により表される。
海水の吸収係数、s(z,λ)は波長λの光に対する海
水の散乱係数、α(z,λ)は波長λの光に対する前方
散乱確率、β(z,λ)は波長λの光に対する後方散乱
確率、η(z,λ)は波長λの光の直達放射照度と拡散
放射照度の比である。上向きの光の拡散放射照度E
u(z,λ)は、上記式(2)〜(5)に示す連立微分
方程式を、海面と海底面との状態を境界条件として解く
ことにより得られるが、水深zの関数として解析的に算
出することは、非常に困難である。そこで本実施形態で
は、物質の鉛直濃度を一定と仮定することにより、拡散
放射輝度Eu(z,λ)を求める。但し、水深に比べて
光の透過距離が十分短いものとし、海面からの反射光は
考慮しないものとする。この場合の拡散放射照度E
u(z,λ)は、以下の式(6)〜(10)で表され
る。
直下の直達放射照度、Ei(λ)は海面から入射する海
面直下の拡散放射照度、ri(λ)は海面での海中拡散
光の反射率、bf(λ)は前方散乱係数、bb(λ)は後
方散乱係数である。また、海面直下の上向きの放射照度
Euo−(λ)は、上記式(6)に、「z=0」を代入す
ることにより、以下に示す式(11)で算出される。
星画像センサは、バンド(チャネル)毎に決められた波
長帯の光を検出する。そのため、衛星では、海面から放
射される光の放射輝度(Luo+(θv,φv,λ))が、
バンド毎の検出波長帯積分値として観測される。即ち、
バンドiで検出される波長帯の光における海面での放射
輝度Qi(sea)は、以下の式(12)で表される。
値、λ2は、バンドiの検出波長の最大値、θv,φvは
それぞれ光の伝達方向の天頂角の方位角と方位角であ
る。また、海面から放射される放射輝度(Luo+
(θv,φv,λ))は、以下に示す式(13)により算
出される。
−(λ)/Luo+(0,φv,λ)」、Esun(θsun,
λ)は海面直上での下向きの直達放射照度、E
sky(λ)は海面直上での下向きの拡散放射照度、T
surfは海中から大気中への透過率、f1は放射輝度変換
を含む下向きの拡散放射照度が海面で反射する確率、f
2は放射輝度変換を含む下向き直達放射照度が海面で反
射する確率である。海面から放射された光Qi(se
a)は、大気による吸収及び散乱の影響を受けながら大
気中を進み、衛星に到達する。衛星では、太陽放射光に
おける大気での散乱だけを受けて衛星に到達する光成分
も検出するため、衛星で観測される放射輝度Qi(sa
t)は、海面での放射輝度Qi(sea)のほかに、大
気や海面での吸収及び散乱によるノイズを含むものとな
る。そこで、このノイズによる影響を衛星で観測される
放射輝度Qi(sat)から取り除くことによって得ら
れる海面からの放射輝度Q’i(sat)は、海面での
放射輝度Qi(sea)と等価となり、バンド数をNと
すれば、以下に示す式(14)で表されるようになる。
ように、海面直下の上向きの放射照度(Euo−(λ))
は、波長λの光に対する海水の吸収係数a(z,λ)、
及び波長λの光に対する海水の散乱係数s(z,λ)に
含まれる各物質の濃度要素mkに非線形に依存するもの
であり、上述の連立方程式群を、濃度要素mkに関して
解析的に解くことは困難である。そこで水中物質の濃度
の推定に際し、まず、衛星画像データからバンド毎に大
気や海面での吸収及び散乱の影響を取り除いた海面から
の放射輝度Q’i(sat)を算出する。この場合、上
記式(12)の水面における光の上向き放射輝度Q
i(sea)を、「Q’i(sat)=Qi(sea)」
とする、すなわち衛星画像データによる算出値(第1光
放射量)と理論値(第2光放射量)とを等価にすること
により、吸収係数及び散乱係数に含まれる各水中物質の
濃度要素mkを逆推定することが可能となる。具体的に
は、以下の式(15)で示される評価関数Pを定義する
とともに、当該評価関数を最小化することにより最適な
濃度要素mkを導出する。
m’kに対する水面における光の上向き放射輝度Qi(s
ea)の値、σiはQ’i(sat)に含まれる誤差の標
準偏差である。この非線形最適化手法による上記式(1
5)に基づいて濃度要素mkの反復的な更新を行う。具
体的には、初期値として与えられた試行解から当該試行
解を改良していく過程を定義し、該過程の閾値と評価関
数Pの算出値とを比較して、評価関数Pが減少しなくな
るまで、つまり最小化するまでこの過程を繰り返す。こ
のような反復的な演算アルゴリズムにより最適な濃度要
素mkが導出され、正しい濃度の推定が可能となる。
伝達方程式を水中における光の多流束モデルを用いて近
似化し、上記式(12)の水面における光の上向き放射
輝度Qi(sea)について、画像データによる算出値
(第1光放射量)と理論値(第2光放射量)とを等価に
して、上記式(15)のような濃度の評価関数、すなわ
ち濃度推定関数を作成する。そして、この濃度推定関数
に基づいて、吸収及び散乱係数に含まれる各水中物質の
濃度についての解、すなわち濃度の逆推定を行う。
上記式(1)〜(5)と比較して解析が容易になるよう
に定式化されており、衛星画像データや大気状態等の撮
像環境情報を取り直した場合であっても迅速に最適解が
得られるため、撮像領域や撮像環境情報の変化にフレキ
シブルに対応することができ、トゥルースデータとの濃
度誤差を一定値以下に低減させることが可能となる。
例を説明する。図1は、本発明の濃度推定装置が具備す
る機能ブロック構成図である。図中、実線は処理の流
れ、破線はデータの流れを表す。この濃度推定装置1
は、オペレーティングシステム(OS)を搭載したコン
ピュータ装置によって実現される。すなわち、当該コン
ピュータ装置が所定のプログラムを読み込んで実行する
ことにより形成される、画像データ取得部11、撮像環
境情報取得部12、水中物質光学特性情報取得部13、
データ変換処理部14、水面上向き放射量算出部(1)
15、水面下向き放射量算出部16、水面上向き放射量
算出部(2)17、及び水中物質濃度推定部18を備え
て濃度推定装置1を構成している。
部または外部記憶装置に格納されて随時読み取られて実
行されるようになっており、当該コンピュータ装置上に
おいて上記各機能ブロック11〜18を形成可能な形態
で記録されたものである。但し、使用時にコンピュータ
装置上で各機能ブロック11〜18が形成されれば本発
明を実施することができるので、例えばインストール前
は、コンピュータ装置と分離可能なCD−ROM(コン
パクトディスク型ROM)やFD(フレキシブルディス
ク)等の可搬性記録媒体、或いは構内ネットワークに接
続されたプログラムサーバ等に記録されているものであ
ってもよい。
データ入力用インタフェースを介して入力される、水中
物質の濃度の推定対象となる特定水域の衛星画像データ
を読み込んでデータ変換処理部14へ送出するものであ
る。上記画像データ入力用インタフェースは、例えば、
パラボラアンテナを含む受信装置のように、人工衛星か
らの衛星画像データを入力する手段であるが、予め人工
衛星から取り込んだ衛星画像データを保持する外部記憶
装置等を用いても良い。
よる撮像の際に、緑バンド、赤バンド、近赤外バンド、
赤外バンド、熱赤外バンド等を含む光の観測波長帯域毎
に分割して観測された複数の画像データである。撮像環
境情報取得部12は、撮像環境情報入力用インタフェー
スを介して入力される撮像環境情報を読み込んで水面上
向き放射量算出部(1)15及び水面下向き放射量算出
部16に送出するものである。この撮像環境情報は、画
像データの解像度、図示しない画像センサによる検出波
長、大気状態、太陽位置等を含む、画像データ撮像時の
環境に関する観測情報である。また、上記撮像環境情報
入力用インタフェースは、例えば、撮像環境情報を予め
保持した外部記憶装置あるいは当該情報を個別的に入力
するためのキーボード等である。
用インタフェースを介して入力される水中物質の光学特
性情報を読み込んで水中物質濃度算出部17へ送出する
ものである。光学特性情報は、具体的には、水中に存在
する溶存物質や懸濁物質の光に対する粒子レベルのミク
ロな反射特性及び吸収特性を含む情報である。また、水
中物質光学特性情報入力用インタフェースは、例えば、
水中物質群の光学特性情報を予めデータベースの形態で
保持した外部記憶装置あるいは当該情報を個別的に入力
するためのキーボード等を用いて適宜構成される。
を画素毎の輝度値に変換するとともに、各輝度値を図示
しない記憶手段に格納されている所定の換算式により衛
星画像データ毎の放射輝度データに変換するものであ
る。この換算式は、個々の衛星画像データにおける光の
波長帯域の種類に応じて定まるものであり、例えば、N
ASA(米国航空宇宙局)等が提案している公知の換算
式を用いることができる。変換された放射輝度データ
は、水面上向き放射量算出部(1)15に入力される。
星画像データにおける放射輝度データ及び撮像環境情報
に基づいて、特定水域の水面から実際に放射されてい
る、光の観測波長帯域毎の放射量(以下、便宜上、「観
測放射量」と称する)を算出するものである。ここで
は、水面から上向きに放射されている光の放射輝度(以
下、「上向き放射輝度」)を算出する。この上向き放射
輝度は、例えば、光子に関する公知のモンテカルロ・シ
ュミレーションモデルを用いて算出することができる。
この算出値は、水中物質濃度推定部18に入力されるよ
うになっている。
情報取得部12における撮像環境情報に基づいて、水面
に入射する下向きの光の波長帯域毎の放射量(以下、
「水面下向き放射輝度」)、すなわち、水面直上での下
向き太陽直達光及び拡散光の強さを算出するものであ
る。この水面下向き放射輝度は、例えば、上述の光子に
関する所定のモンテカルロ・シュミレーションモデルを
用いて算出することができる。この算出値は、水面上向
き放射量算出部(2)17に入力されるようになってい
る。
面下向き放射輝度と水中物質群の光学特性情報とに基づ
いて、特定水域の水面から上向きに放射される、光の観
測波長帯域毎の理論上の放射量(以下、便宜上、「理論
放射量」と称する)を算出するものである。具体的に
は、光の放射伝達現象を、水中における光の多流束モデ
ルに基づいて定式化し、上記式(12)のような、上向
きの放射輝度算出式から放射輝度を算出するように構成
される。
5、水面上向き放射量算出部(2)17、水面下向き放
射量算出部16における各放射量の算出に関する技術に
ついては、例えば、「画像解析ハンドブック、高木幹雄
他 著、東京大学出版会」等に紹介されている。
して算出された観測放射量と理論放射量とに基づいて、
水中に存在する物質の濃度を推定するものである。実際
には、上記式(15)のような濃度の評価関数Pを作成
し、この評価関数Pを最小化(最適化)させることで、
吸収及び散乱係数に含まれる各水中物質の濃度mkの逆
推定を行う。この逆推定は、例えば、非線形最小2乗問
題の標準的解法である、Levenberg-Marquardt法を適用
することによって実現が可能である。推定された濃度結
果は、図示しない出力用インタフェースを介して適宜出
力される。この出力用インタフェースは、例えば、CR
T表示装置やプリンタのような可視化手段、あるいは外
部記憶装置等を用いて構成される。
水域における水中物質の濃度を推定する場合の手順につ
いて、図2〜5を参照して説明する。
説明する。濃度推定装置1では、画像データ取得部11
から特定水域における複数の衛星画像データ、撮像環境
情報取得部12から当該衛星画像データの撮像時におけ
る撮像環境情報、水中物質光学特性取得部13から水中
物質光学特性情報を、それぞれ対応する入力用インタフ
ェースを介して取得する(ステップS101)。取得し
た各画像データは、データ変換処理部14により、画素
毎の輝度値に変換され、さらに、所定の換算式を用いて
各々放射輝度データに変換される(ステップS10
2)。水面上向き放射量算出部15では、放射輝度デー
タ及び撮像環境情報に基づいて、光の観測放射量を算出
する(ステップS103)。また、水面下向き放射量算
出部16で、撮像環境情報に基づく水面下向き放射輝度
を算出するとともに、水面上向き放射量算出部(2)1
7で、当該水面下向き放射輝度と水中物質光学特性情報
と基づく光の理論放射量を算出する(ステップS10
4)。その後、上記理論放射量と観測放射量とに基づい
て水中物質濃度推定部18で、水中物質の濃度を推定す
る(ステップS105)。
プS104で用いる下向き放射輝度は、モンテカルロ・
シュミレーションモデルに基づいて算出することができ
る。以下、この場合の手法を、図3及び4を参照して詳
しく説明する。
入力用インタフェースを介して撮像環境情報を取得する
(ステップS201)。水面下向き放射量算出部16
は、取得した撮像環境情報に基づいて光の一次反射位置
を決定する(ステップS202)。次に、水面下向き放
射量算出部16は、撮像環境情報における太陽入射角に
従って、また、水面上向き放射量算出部(1)15の場
合は任意の方向へ光子を入射する(ステップS20
3)。この光子は、所定の乱数発生アルゴリズムに基づ
いて発生されるものである。
の場合には(ステップS204:No)、該光子数に基づ
いて自由行程距離を算出するとともに(ステップS20
5)、当該自由行程距離から、光子の終点座標を決定す
る(ステップS206)。決定された終点座標が大気上
面ならば、光子数をカウントしてステップS203に戻
り(ステップS207:大気上面、S208)、光子を
入射する。また、終点座標が水面以下ならば、光子数を
カウントしてステップS203に戻り、光子を入射する
(ステップS207:水面、S208)。
ップS207:大気中)、大気状態を判定する。大気状
態が吸収ならば(ステップS209:吸収)、ステップ
S203に戻り、光子を入射するとともに、大気状態が
散乱ならば(ステップS209:散乱)、散乱角を算出
して(ステップS210)ステップS203に戻り、光
子を入射する。ステップS204において、光子数が予
め設定された数以上ならば(ステップS204:Ye
s)、図4に移り、発生光子数及びカウント光子数から
光子数の割合、すなわち光子の到達確率を算出する(ス
テップS211)。
ステップS201で取得された撮像環境情報における太
陽放射照度に基づいて、拡散放射照度及び直達放射照度
を算出する(ステップS212:下向き、S213)。
では、画像データ入力用インタフェースを介して画像デ
ータ取得部11により衛星画像データを取得する((ス
テップS212:上向き、S214)。そして、データ
変換処理部14で変換された光の観測波長帯域毎の衛星
画像データに関する放射輝度データに基づいて、水面上
向き放射輝度Q’i(sat)、すなわち光の観測放射
量を算出する(ステップS215)。
質の濃度推定方法を、図5を参照して説明する。濃度推
定装置1は、まず、水面上向き放射量算出部(2)17
において、水面下向き放射量算出部16により海面直上
での下向き太陽直達光の放射照度Esun及び拡散光の放
射照度Eskyとを取得する(ステップS301〜30
2)。また、水中物質光学特性情報取得部13におい
て、水中物質光学特性情報入力用インタフェースを介し
て水中物質の光学特性情報を取得する(ステップS30
3)。
また、導出すべき水中物質の濃度に対する初期値を設定
するとともに(ステップS304)、取得した上記太陽
直達光の放射照度Esun、拡散光の放射照度Esky、水中
物質の光学特性情報、及び設定された水中物質の濃度初
期値に基づいて、水面上向き放射輝度Q’i(se
a)、すなわち光の理論放射量を算出する(ステップS
305)。
プS305により算出された光の理論放射量Q’i(s
ea)と衛星画像データにおける観測放射量Q’i(s
at)との差分に基づいた評価値Pを算出する(ステッ
プS306)。この評価値は、例えば非線形最小2乗問
題に即した上述の濃度推定関数に基づいて算出されるも
のである。そして、算出された評価値Pが最小となるま
で濃度を更新する(ステップS307:No〜308)。
さらに、ステップS305に戻って、光の理論放射量
Q’i(sea)を算出して処理を繰り返す。このステ
ップS305〜308での処理により、光の理論放射量
における吸収係数及び散乱係数に含まれる水中物質の最
適な濃度要素が導出される。なお、本処理手順図では、
処理順序を便宜的に設定して説明したものであり、実際
の処理では、例えば、ステップS301〜304は、並
列的に行われる。
は、簡易な手法によって最適な濃度が導出されるので、
経験的な回帰関係式や統計的手法に即した従来手法と比
較して、水域毎に膨大なトゥルースデータを再度取得す
ることなく、推定対象となる水域全体における水中物質
の濃度推定精度を一定値以上維持することが可能とな
る。
数の汚染物質が混在するような沿岸域を濃度推定特定水
域とする場合であっても、上述の濃度推定関数を用いる
ことにより、複数の水中物質の濃度を同時且つ高精度に
推定することが可能となる。
となく、濃度の推定精度を一定値以上維持することが可
能となり、水域全体に対するロバスト性が向上する。こ
のことから、従来手法と比較してトゥルースデータの収
集時間が短縮され、迅速な濃度推定が可能となるととも
に、トゥルースデータと濃度推定値との濃度誤差を低減
させることが可能になる。従って、撮影水域や撮影時の
大気状態が頻繁に変化する場合であっても、濃度推定値
を迅速に算出できるようになり、現地観測データに対す
る濃度誤差を一定値以下に維持することができる。
によれば、衛星画像データに係る撮像領域や撮像環境が
変化した場合であっても、当該特定水域における濃度の
推定精度を一定以上に維持することが可能となる。従っ
て、光学理論的な枠組みに基づいた合理的且つ汎用性の
高い水中物質の濃度推定が容易に実現可能となる特有の
効果がある。
ロック図。
図。
理手順図。
理手順図。
する処理手順図。
Claims (11)
- 【請求項1】 特定水域をその上空から複数の波長帯域
で撮影して得たマルチスペクトル画像データと撮影時の
環境情報とを取得して前記特定水域における水中物質の
濃度を導出する装置における濃度推定方法であって、 前記マルチスペクトル画像データから前記特定水域にお
ける第1光放射量を導出するとともに、前記環境情報及
び既定の水中物質群の光学特性情報から前記特定水域に
おける第2光放射量を導出し、導出した前記第1光放射
量および第2光放射量に基づいて光学特性が表す水中物
質の濃度を推定する処理を、 前記装置上で実行することを特徴とする、水中物質の濃
度推定方法。 - 【請求項2】 特定水域をその上空から複数の波長帯域
で撮影して得たマルチスペクトル画像データと撮影時の
環境情報とを取得して前記特定水域における水中物質の
濃度を導出する装置における濃度推定方法であって、 前記マルチスペクトル画像データから前記特定水域にお
ける第1光放射量を導出するとともに、前記環境情報及
び既定の水中物質群の光学特性情報から前記特定水域に
おける第2光放射量を導出する処理と、 光の放射伝達方程式を水中における光の多流束モデルを
用いて近似化することにより前記第1光放射量及び前記
第2光放射量を要素に含む濃度推定関数を作成し、この
濃度推定関数に前記導出した第1光放射量及び前記第2
光放射量を適用して前記水中物質の濃度を推定する処理
と、 を前記装置上で実行することを特徴とする、水中物質の
濃度推定方法。 - 【請求項3】 前記濃度推定関数は、前記第1光放射量
及び第2光放射量が前記水中物質の濃度に非線形に依存
する形態のものであることを特徴とする、 請求項1または2記載の濃度推定方法。 - 【請求項4】 前記光学特性情報は、前記水中物質の光
に対する粒子レベルの反射特性及び吸収特性を含む情報
であることを特徴とする、 請求項1ないし3のいずれかの項記載の濃度推定方法。 - 【請求項5】 特定水域をその上空から複数の波長帯域
で撮影して得たマルチスペクトル画像データ及び撮影時
の環境情報を取得する取得手段と、 取得した前記マルチスペクトル画像データを放射輝度デ
ータに換算するとともに、換算された放射輝度データに
基づいて前記特定水域における第1光放射量を導出する
第1光放射量導出手段と、 前記環境情報及び既定の水中物質群の光学特性情報から
前記特定水域における第2光放射量を導出する第2光放
射量導出手段と、 光の放射伝達方程式を水中における光の多流束モデルを
用いて近似化することにより前記第1光放射量及び前記
第2光放射量を要素に含む濃度推定関数を作成し、この
濃度推定関数に前記導出した観測放射量及び理論放射量
を適用して前記濃度を推定する濃度推定手段と、 を備えたことを特徴とする、水中物質の濃度推定装置。 - 【請求項6】 前記第2光放射量導出手段は、前記環境
情報から前記特定水域の水面に入射する光の下向き放射
量を算出するとともに、前記算出した下向き放射量及び
既定の水中物質群の光学特性情報に基づいて前記特定水
域の水面から放射される光の放射量を算出し、これらの
算出結果に基づいて前記第2光放射量を導出することを
特徴とする、 請求項5記載の濃度推定装置。 - 【請求項7】 前記第2光放射量導出手段は、前記特定
水域の水面から放射される、前記波長帯域毎の最大波長
及び最小波長の放射輝度の積分値を算出し、この積分値
に基づいて前記第2光放射量を導出することを特徴とす
る、 請求項5記載の濃度推定装置。 - 【請求項8】 前記第1光放射量導出手段及び前記第2
光放射量導出手段が、光子に関する所定のモンテカルロ
・シミュレーションモデルに基づいて構成されているこ
とを特徴とする、 請求項5乃至7のいずれかの項記載の濃度推定装置。 - 【請求項9】 前記濃度推定手段は、前記第1光放射量
及び第2光放射量が水中物質の濃度に対して非線形に依
存する、所定の非線形最小2乗問題に基づいて前記濃度
推定関数を作成し、前記第2光放射量における吸収係数
及び散乱係数に含まれる水中物質の濃度要素を前記濃度
推定関数に適用して濃度を推定するとともに、この推定
結果に即して前記濃度要素を更新するように構成されて
いることを特徴とする、 請求項5記載の濃度推定装置。 - 【請求項10】 前記濃度推定手段は、前記第2光放射
量と前記第1放射量との差分に基づく評価値が最小にな
るまで前記濃度要素を反復的に更新することを特徴とす
る、 請求項8記載の濃度推定装置。 - 【請求項11】 特定水域をその上空から複数の波長帯
域で撮影して得たマルチスペクトル画像データと撮影時
の環境情報とを取得する機能を有するコンピュータ装置
に読み取られたときに、前記コンピュータ装置に少なく
とも下記の処理を実行させるためのプログラムが記録さ
れた記録媒体。 (1)前記マルチスペクトル画像データから前記特定水
域における第1光放射量を導出するとともに、前記環境
情報及び既定の水中物質群の光学特性情報から前記特定
水域における第2光放射量を導出する処理、(2)光の
放射伝達方程式を水中における光の多流束モデルを用い
て近似化することにより前記観測放射量及び前記理論放
射量を要素に含む濃度推定関数を作成し、前記導出した
観測放射量及び前記理論放射量を前記作成した濃度推定
関数に適用して前記水中物質の濃度を推定する処理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP30404098A JP3348834B2 (ja) | 1998-10-26 | 1998-10-26 | 水中物質の濃度推定方法及び装置、記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP30404098A JP3348834B2 (ja) | 1998-10-26 | 1998-10-26 | 水中物質の濃度推定方法及び装置、記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000131218A true JP2000131218A (ja) | 2000-05-12 |
JP3348834B2 JP3348834B2 (ja) | 2002-11-20 |
Family
ID=17928337
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP30404098A Expired - Lifetime JP3348834B2 (ja) | 1998-10-26 | 1998-10-26 | 水中物質の濃度推定方法及び装置、記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3348834B2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007225358A (ja) * | 2006-02-22 | 2007-09-06 | National Institute For Agro-Environmental Science | 陸水中の土壌由来懸濁物濃度の推定方法 |
JP2010251958A (ja) * | 2009-04-14 | 2010-11-04 | Mitsubishi Electric Corp | 画像処理装置 |
CN105203466A (zh) * | 2015-09-17 | 2015-12-30 | 中国科学院南京地理与湖泊研究所 | 一种富营养化湖泊非藻华条件下藻类总存量遥感估算方法 |
CN105300862A (zh) * | 2015-11-13 | 2016-02-03 | 金陵科技学院 | 基于云处理车载移动大气颗粒污染物的环境检测方法及*** |
CN110108609A (zh) * | 2019-05-07 | 2019-08-09 | 国网湖南省电力有限公司 | 一种基于卫星多通道的pm2.5浓度值分布情况获取方法及*** |
CN112378828A (zh) * | 2020-12-11 | 2021-02-19 | 中科三清科技有限公司 | 基于卫星遥感数据反演大气细颗粒物浓度的方法及装置 |
-
1998
- 1998-10-26 JP JP30404098A patent/JP3348834B2/ja not_active Expired - Lifetime
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007225358A (ja) * | 2006-02-22 | 2007-09-06 | National Institute For Agro-Environmental Science | 陸水中の土壌由来懸濁物濃度の推定方法 |
JP4719915B2 (ja) * | 2006-02-22 | 2011-07-06 | 独立行政法人農業環境技術研究所 | 陸水中の土壌由来懸濁物濃度の推定方法 |
JP2010251958A (ja) * | 2009-04-14 | 2010-11-04 | Mitsubishi Electric Corp | 画像処理装置 |
CN105203466A (zh) * | 2015-09-17 | 2015-12-30 | 中国科学院南京地理与湖泊研究所 | 一种富营养化湖泊非藻华条件下藻类总存量遥感估算方法 |
CN105300862A (zh) * | 2015-11-13 | 2016-02-03 | 金陵科技学院 | 基于云处理车载移动大气颗粒污染物的环境检测方法及*** |
CN110108609A (zh) * | 2019-05-07 | 2019-08-09 | 国网湖南省电力有限公司 | 一种基于卫星多通道的pm2.5浓度值分布情况获取方法及*** |
CN110108609B (zh) * | 2019-05-07 | 2021-11-02 | 国网湖南省电力有限公司 | 一种基于卫星多通道的pm2.5浓度值分布情况获取方法及*** |
CN112378828A (zh) * | 2020-12-11 | 2021-02-19 | 中科三清科技有限公司 | 基于卫星遥感数据反演大气细颗粒物浓度的方法及装置 |
CN112378828B (zh) * | 2020-12-11 | 2021-09-17 | 中科三清科技有限公司 | 基于卫星遥感数据反演大气细颗粒物浓度的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
JP3348834B2 (ja) | 2002-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5233541A (en) | Automatic target detection process | |
Fournier et al. | Computer-based underwater imaging analysis | |
WO2018037920A1 (ja) | 画像処理装置、画像処理方法及びコンピュータ読み取り可能記録媒体 | |
Frederickx et al. | A forward scattering dipole model from a functional integral approximation | |
Polonsky et al. | Wide-angle imaging lidar deployment at the ARM Southern Great Plains site: Intercomparison of cloud property retrievals | |
JP3348834B2 (ja) | 水中物質の濃度推定方法及び装置、記録媒体 | |
Montes et al. | New algorithm for atmospheric correction of hyperspectral remote sensing data | |
CN117853950B (zh) | 基于光谱库和机器学习的多光谱卫星图像处理方法及*** | |
Nevis et al. | Advantages of three-dimensional electro-optic imaging sensors | |
Stephan et al. | Model based image restoration for underwater images | |
CN113888420A (zh) | 一种基于修正模型的水下图像复原方法、设备及存储介质 | |
CN111640117A (zh) | 一种寻找建筑物渗漏源位置的方法 | |
Meyers et al. | Incorporation of polarization into the DIRSIG synthetic image generation model | |
JPH06300845A (ja) | 画像情報処理装置 | |
Sundberg | Forward modeling of cloud shadows and the impact of cloud shadows on remote sensing data products | |
Mariano et al. | Hyperspectral material identification on radiance data using single-atmosphere or multiple-atmosphere modeling | |
US20240125935A1 (en) | Arithmetic operation system, training method, and non-transitory computer readable medium storing training program | |
US20240126953A1 (en) | Arithmetic operation system, training method, and non-transitory computer readable medium storing training program | |
US20240126952A1 (en) | Arithmetic operation system, training method, and non-transitory computer readable medium storing training program | |
Raqueño et al. | A model-based approach to hyperspectral analysis of water constituents | |
Caimi et al. | Real-time 3D underwater imaging and mapping using a laser line scan technique | |
Adler-Golden et al. | Cloud effects in hyperspectral imagery from first-principles scene simulations | |
Richtsmeier et al. | Full-spectrum scene simulation | |
Zege et al. | Recent advancement in computer simulating performance of the sounding and detecting ocean lidars in very shallow waters | |
de Paula et al. | Bathymetric DEM derived from Lyzenga’s algorithm (multispectral Bathymetry) using Landsat 8 images: a case study of Fort Lauderdale coast |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313115 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070913 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080913 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080913 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090913 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090913 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100913 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100913 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110913 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110913 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120913 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120913 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130913 Year of fee payment: 11 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |