CN101493409A - 一种自动提取水体污染信息的方法 - Google Patents

Abstract

本发明公开了一种自动提取水体污染信息的方法，包括以下步骤：(1)计算获得被测水域的出水反射率；(2)建立水质遥感模型；(3)利用水体光谱曲线，根据上述水质遥感模型反演得到水质参数。通过对该方法进行大量应用试验，结果表明该方法具有较高的精度，基本达到实用化的要求。试验结果图像可以很好地反演区域水质的空间分布情况，水质污染的分布特征与所掌握的情况相符。可以很好的反映跨境污染的情况。可以较好地反映大、中型污染源情况，包括污染源位置、污染的规模及其影响范围，对污染的控制与治理可提供极有用的信息。

Description

一种自动提取水体污染信息的方法

技术领域

本发明属于水质遥感技术领域，特别涉及一种通过遥感反演自动提取水体污染信息的方法。

背景技术

遥感技术可以很好反映大范围内的水体水质状况，在区域水环境监测方面具有极大的潜力，而遥感反演的定量化程度直接决定了水环境遥感监测的效能。故水质遥感的定量极为重要。早期水污染遥感是根据污染水域色调变化的程度来对污染情况进行定性调查。

现阶段主要采用两类方法。一类是采用经验模型方法，即测量各种水体的光谱特征，并用回归分析等方法建立遥感数据的某个可见光波段或波段的比值与污染浓度之间经验公式，以此来对水污染信息进行定量提取。这种方法简单易行，物理意义不清楚，对水面实测数据及其与遥感数据的同步性依赖较大。此外，该方法假定海水其它性质基本不变。该假设在河口地区显然很难得到满足，当海水含沙量变化时，该方法的精度会受到显著影响。另一类是采用Dekker A G.1993提出的半经验公式 $R\left(0\right)\left(\mathrm{\λ}\right)=f*\frac{b_b\left(\mathrm{\λ}\right)}{a\left(\mathrm{\λ}\right)+b_b\left(\mathrm{\λ}\right)}$ 来进行计算，该方法虽在一定程度上反映了物理学规律，但物理意义并不会完清晰，公式内的f参数需要用经验公式确定，故实适性仍受局限，实际应用中通常效果并不好。解决此问题的最好方法是建立反射光谱与污染浓度之间的物理模型并将混浊度列为变量。

本申请的发明人曾经根据水体反射光谱产生的物理过程，建立污染浓度与水体反射光谱之间的物理模型，并应用于珠江中水域的污染信息提取，取得了较好的效果，但是该模型没有考虑黄色物质的作用，对悬浮的散射与反射特征的考虑也不够全面，在水混浊度较高的情况下，水质参数计算的精度仍受较大影响。

发明内容

本发明的目的在于提供一种具有更高的精度和稳定性的水体污染信息提取方法，用于水体混浊度、COD和叶绿素浓度和水深的定量遥感反演。

本发明所采用的技术方案：一种自动提取水体污染信息的方法，包括以下步骤：

(1)计算获得被测水域的出水反射率；

(2)建立水质遥感模型；

(3)利用水体光谱曲线，根据上述水质遥感模型反演得到水质参数。

上述水质参数包括清洁水、悬浮泥沙、综合污染物和叶绿素的散射率和吸收系数。为便于反演计算，在测量时分别选取深度较大的清洁水、混浊水、重污染水和高叶绿素含量水体进行测量，然后消除水面镜面反射的影响，计算出水反射率，再进行散射率和吸收系数的反演。

上述步骤(1)中的出水反射率R_wo(λ)的计算公式为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{R\left(\mathrm{\λ}\right)-{\mathrm{\ρ}}_{\mathrm{ra}}}{[T_a(\mathrm{\θ},\mathrm{\λ})+{\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})/\cos \mathrm{\θ}]T_a(0,\mathrm{\λ})T_{\mathrm{ws}}}-R_{\mathrm{TM}5}$

式中：R(λ)为遥感数据转换得以到的行星反射率，ρ_ra为大气上行散射反射率，T_a(θ，λ)为大气在入射方向上的透过率，ω′_ad(θ，λ)为下行散射光的半球散射系数，θ为太阳天顶角，T_ws为水面透过率，R_TM5为TM₅计算的反射率。

上述步骤(2)中的水质遥感模型为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})}(1-e^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H})+\frac{1}{\cos \mathrm{\θ}}{R}_b{e}^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H}$

式中：R_wo(λ)为出水反射率，ω(Θ)为水层垂向上的散射比，Θ为水层垂向上的散射角，θ为水体中的光入射角，α为水体吸收系数，β为水体散射系数，R_b为水底底质反射率，H为水深， $\mathrm{\μ}=1+\frac{1}{\cos \mathrm{\θ}}.$

上述步骤(3)中的水体光谱曲线是用ASD连续光谱仪在水面测量获得。

本发明所述水体污染信息的提取方法，通过对该方法进行大量应用试验，包括对珠江三角洲进行多时相水污染遥感监测和广东省全省大中型水库的水质遥感监测。结果表明该方法具有较高的精度，基本达到实用化的要求。试验结果图像可以很好地反演区域水质的空间分布情况，水质污染的分布特征与所掌握的情况相符。可以很好的反映跨境污染的情况。可以较好地反映大、中型污染源情况，包括污染源位置、污染的规模及其影响范围，对污染的控制与治理可提供极有用的信息。

附图说明

图1是清洁水体散射系数与反射率的曲线图；

图2是清洁水体散射系数的曲线图；

图3是黄色物质吸收、散射系数的曲线图；

图4是叶绿素、悬浮泥沙吸收、散射系数的曲线图；

图5是耗氧性有机物吸收系数的曲线图；

图6是1998年12月22日珠江三角洲水污染(COD)分布遥感图(卫星遥感数据：TM)；

图7是2002年11月7日珠江三角洲水网悬浮物混浊度分布遥感图(卫星遥感数据：TM)；

图8是2001年11月20日珠江三角洲水污染(COD)遥感分布图(卫星遥感数据：TM)；

图9是2006年11月10日珠江三角洲水污染(COD)遥感分布图(卫星遥感数据：TM)；

图10是2001年11月20日珠江三角洲主要水污染参数遥感反演结果验证——COD遥感反演值与测量值相关性分析；

图11是2001年11月20日珠江三角洲主要水污染参数遥感反演结果验证——综合污染浓度遥感结果与测量值相关分析；

图12是2006年11月10日珠江三角洲主要水污染参数遥感反演结果验证——COD计算值与测量值对比；

图13是2006年11月10日珠江三角洲主要水污染参数遥感反演结果验证——遥感反演综合浓度与测量浓度对比。

具体实施方式

下面结合附图对本发明的具体结构作进一步的描述。

1.技术原理

1.1水域出水反射率的遥感计算

首先对卫星遥感数据进行辐射定标。水域之上卫星遥感仪器接收到的短波辐射由水面辐射和大气散射辐射组成，可表示为：

$L\left(\mathrm{\λ}\right)=\frac{1}{\mathrm{\π}}[F_s\left(\mathrm{\λ}\right)+E_d\left(\mathrm{\λ}\right)]R_w\left(\mathrm{\λ}\right)T_{\mathrm{ar}}(0,\mathrm{\λ})T_{\mathrm{ws}}+L_{\mathrm{ar}}\left(\mathrm{\λ}\right)---\left(1\right)$

其中，λ为波长，F_s为水面上太阳直射光的辐照度、E_d为由水面上天空散射光辐照度，R_w水综合反射(包括水面反射和水体散射)率，T_ar(0)为垂向上的大气透过率，T_ws为水面透过率，是表征水面镜面反射对入射光的削弱的一个参数。L_ar为整层大气上行散射辐亮度。

水面太阳直射光辐照度为：

$E_s\left(\mathrm{\λ}\right)=T_{\mathrm{ar}}(\mathrm{\θ},\mathrm{\λ})F_0\left(\mathrm{\λ}\right)\frac{{r_0}^2}{r}\cos \mathrm{\θ}---\left(2\right)$

T_a(θ，λ)大气在入射方向上的透过率，θ为太阳天顶角；E₀为大气层外太阳平均辐照度；r₀、r分别及平均及成像时的日地距离，二者之比可从天文年历中查到。一般情况下r₀/r≈1。

整层大气上行散射辐亮度为：

L_ar(λ)＝L_r(λ)+L_a(λ)

L_r、L_a大气层瑞利散射和气溶胶散射上行辐射，由下式计算(Gordon，1989)：

L_x＝F₀τ_xω_xP_x/(4πcosθ_v)

下标x为r时表示瑞散射，为a时表示气溶胶散射；τ_x，ω_x，P_x分别为大气分子或气溶胶散射的光学厚度、单次散射比和散射相函数。整层大气上行散射为：

$L_{\mathrm{ar}}\left(\mathrm{\λ}\right)=F_0\left(\mathrm{\λ}\right)\frac{{\mathrm{\ω}}_a\left(\mathrm{\λ}\right){\mathrm{\τ}}_a{P}_a\left(\mathrm{\Θ}\right)+{\mathrm{\ω}}_r\left(\mathrm{\λ}\right){\mathrm{\τ}}_r{P}_r\left(\mathrm{\Θ}\right)}{4\mathrm{\π}}$

Θ为散射角。令： ${\mathrm{\ρ}}_{\mathrm{ra}}=\frac{{\mathrm{\ω}}_a\left(\mathrm{\λ}\right){\mathrm{\τ}}_a{P}_a\left(\mathrm{\Θ}\right)+{\mathrm{\ω}}_r\left(\mathrm{\λ}\right){\mathrm{\τ}}_r{P}_r\left(\mathrm{\Θ}\right)}{4},$ 为大气上行散射反射率。

入射天空光是天空散射光的半球积分，可表示为：

ω_a为整层大气的散射系数，P_a(Θ)为大气散射相函数；Θ为散射角，是θ，α的函数。令：

为下行散射光的半球散射系数，有：

$E_d\left(\mathrm{\λ}\right)=\frac{F_0}{\mathrm{\π}}\left(\mathrm{\λ}\right){\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})---\left(3\right)$

T_ws为水面透过率，表征水面镜面反射对入射光的削减，由下式计算：

T_ws＝(1-ρ_ws)cosθ             (4)

式中：ρ_ws为水面各方向太阳镜面反射率的平均值，通常取ρ_ws＝0.06。

水面向上的辐射由两部分组成：一部分为出自水中的辐射，包括水体散射和水底反射，相应的反射率为出水反射率，记为R_wo；另一部分为水面的镜面反射，相应的反射率记为R_m，故水综合反射率为：

R_w＝R_wo+R_m            (5)

将(2)～(5)式代入(1)式得：

$L\left(\mathrm{\λ}\right)=\frac{F_0\left(\mathrm{\λ}\right)}{\mathrm{\π}}\left\{\right[T_a(\mathrm{\θ},\mathrm{\λ})\cos \mathrm{\θ}+{\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})](R_{\mathrm{wo}}\left(\mathrm{\λ}\right)+R_m)T_a(0,\mathrm{\λ})T_{\mathrm{ws}}+{\mathrm{\ρ}}_{\mathrm{ra}}\left(\mathrm{\λ}\right)\}---\left(6\right)$

行星反射率： $R\left(\mathrm{\λ}\right)=\frac{\mathrm{\πL}\left(\mathrm{\λ}\right)}{F_0\left(\mathrm{\λ}\right)\cos \mathrm{\θ}},$ 将(6)代入得：

入射方向、观测方向和波长一定时，T_a、ω′_ad和ρ_ra取决于大气分子和大气中的气溶胶的性质及浓度。大气分子浓度和性质恒定且已知，气溶胶的性质(粒度和吸收率等)则变化较大。但其类型与所在地区的地理位置、地貌类型和季节密切相关。目前国际上已将全球气溶胶按地域分成海洋型、大陆型、城市型、郊区型等十多种类型，对每种气溶胶的性质均作了测定，并有相应的程序进行计算。进行大气纠正后得到的反射率图像象元值为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{R\left(\mathrm{\λ}\right)-{\mathrm{\ρ}}_{\mathrm{ra}}}{[T_a(\mathrm{\θ},\mathrm{\λ})\cos +{\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})/\cos \mathrm{\θ}]T_a(0,\mathrm{\λ})T_{\mathrm{ws}}}-R_m$

R(λ)为遥感数据转换得以到的行星反射率。水面镜面反射比R_m与波长无关。对红外波段的TM₅和TM₇，水为强吸收，R_wo基本为零，故 $R_{T{M}_5}\≈R_m,$ 因此每个波段减去TM₅计算的反射率R_TM5，可得到出水反射率R_wo(λ)，即：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{R\left(\mathrm{\λ}\right)-{\mathrm{\ρ}}_{\mathrm{ra}}}{[T_a(\mathrm{\θ},\mathrm{\λ})+{\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})/\cos \mathrm{\θ}]T_a(0,\mathrm{\λ})T_{\mathrm{ws}}}-R_{\mathrm{TM}5}---\left(8\right)$

1.2水质遥感模型

水的反射光谱来自水分子和水中悬浮颗粒物的散射光R_ws(λ)和水底反射光R_g(λ)。即：

R_wo(λ)＝R_ws(λ)+R_g(λ)          (9)

二者传播到水面过程中受水的吸收、散射而削减。纯水的反射率和吸收率是恒定的，受污染或含其它杂质后水体的散射和吸收性质发生变化，且不同的波段变化不同。例如，珠江河口水域对水反射率产生重要影响的主要是水中悬浮泥沙和污染物，叶绿素等其它物质因含量有限而影响不大。污染物主要来自工业污染和生活污水，污染源众多，并经充分混合，导致污染物光谱特征较为均一。如只考虑一次散射和污染物与悬浮泥沙的影响，水深为h处的薄层水的一散射对水面出射光的贡献为：

$\mathrm{dL}\left(\mathrm{\λ}\right)=\frac{F_0\cos \mathrm{\θ\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\π}}{e}^{-(\mathrm{\α}+\mathrm{\β})(1+\frac{1}{\cos \mathrm{\θ}})h}$

其中：F₀-入射光强；h-水深；β＝β_w+D_yβ_y+D_sβ_s+D_uβ_u，为水体散射系数，是水分子散射、悬浮泥沙散射和污染物散射率之和。β_w、β_y、β_s、β_u分别为水分子、黄色物质、中以泥沙为主的悬浮颗粒和水中污染物的散射系数，D_y、D_u和D_s分别为黄色物质、污染物综合浓度和泥沙浑浊度；

α＝α_w+D_yα_y+D_sα_s+D_uα_u，为水体吸收系数。α_w、α_y、α_s、α_u分别为水、黄色物质、水中泥沙和污染物质的吸收系数。θ为水体中的光入射角；Θ为水层垂向上的散射角，对TM数据其值为ω(Θ)为水层垂向上的散射比。

ω(Θ)＝β_wP_w(Θ)+D_yβ_yP_y(Θ)+D_sβ_sP_s(Θ)+D_uβ_uP_u(Θ)         (10)

P_w(θ)、P_u(θ)、P_s(θ)和P_u(θ)分别为水分子、悬浮泥沙和污染物质的散射相函数。整层水体散射在垂向上对应的反射率R(λ)为：

$R_{\mathrm{ws}}\left(\mathrm{\λ}\right)=\frac{\mathrm{\π}}{F_0\left(\mathrm{\λ}\right)\cos \mathrm{\θ}}\underset{0}{\overset{h}{\∫}}\mathrm{dL}\left(\mathrm{\λ}\right)\mathrm{dh}=\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4(1+\frac{1}{\cos \mathrm{\θ}})(\mathrm{\α}+\mathrm{\β})}(1-e^{(1+\frac{1}{\cos \mathrm{\θ}})(\mathrm{\α}+\mathrm{\β})H})---\left(11\right)$

根据水底反射的物理机制可求出水底反射光的辐亮度方程为：

$L_g=\frac{E_0\cos \mathrm{\θ}\·R_b}{\mathrm{\π}}{e}^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H}$

R_b为水底反射率。由 $R_g\left(\mathrm{\λ}\right)=\frac{\mathrm{\π}{L}_g}{E_0\cos \mathrm{\θ}}$ 得到水底反射光的出水反射率：

$R_g\left(\mathrm{\λ}\right)=\frac{1}{\cos \mathrm{\θ}}{R}_b{e}^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H}---\left(12\right)$

R_b为水底底质反射率。将(11)式代入(9)式得：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})}(1-e^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H})+\frac{1}{\cos \mathrm{\θ}}{R}_b{e}^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H}---\left(13\right)$

其中 $\mathrm{\μ}=1+\frac{1}{\cos \mathrm{\θ}}.$ 该式即为水质遥感模型。

2.水质参数的测量

2.1测量方法

水质参数包括清洁水、悬浮泥沙、综合污染物和叶绿素的散射率和吸收系数。这些参数是用水体光谱曲线，根据上述水质遥感模型反演得到。水体光谱曲线是用ASD连续光谱仪在水面测量获得。为便于反演计算，测量时分别选取深度较大的清洁水、混浊水、重污染水和高叶绿素含量水体进行测量，然后消除水面镜面反射的影响，计算出水反射率，再进行散射率和吸收系数的反演。

3.2测量得到的典型水体光谱测量与散射、吸收系数

测量得到的清洁水体及主要杂质的散射、吸收系数如图1至图5所示。

3.基于多波段遥感数据的模型求解

3.1主要污染物浓渡的计算

对于有一定水浑浊度或一定程度污染的水体，当水深＞3m时，到达水面的水底反射光极其微弱，可以忽略，其效果等同过水深极大。即对式(13)，H→∞，上式右边指数部分趋于0，将(5)式带入，上式简化为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)\≈\frac{{\mathrm{\β}}_w{P}_w\left(\mathrm{\Θ}\right)+D_y{\mathrm{\β}}_y{P}_y\left(\mathrm{\Θ}\right)+D_s{\mathrm{\β}}_s{P}_s\left(\mathrm{\Θ}\right)+D_u{\mathrm{\β}}_u{P}_u\left(\mathrm{\Θ}\right)+D_c{\mathrm{\β}}_c{P}_c\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}({\mathrm{\α}}_w+D_y{\mathrm{\α}}_y+D_s{\mathrm{\α}}_s+D_u{\mathrm{\α}}_u+D_c{\mathrm{\α}}_c+{\mathrm{\β}}_w+D_y{\mathrm{\β}}_y+D_s{\mathrm{\β}}_s+D_c{\mathrm{\β}}_c)}---\left(13\right)$

通常水中污染物以化学污染物质为主，对水光谱的影响以吸收为主，对水体散射性质影响不大；忽略污染物的散射作用，则上式可以进一步简化为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)\≈\frac{{\mathrm{\β}}_w{P}_w\left(\mathrm{\Θ}\right)+D_y{\mathrm{\β}}_y{P}_y\left(\mathrm{\Θ}\right)+D_s{\mathrm{\β}}_s{P}_s\left(\mathrm{\Θ}\right)+D_c{\mathrm{\β}}_c{P}_c\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}({\mathrm{\α}}_w+D_y{\mathrm{\α}}_y+D_s{\mathrm{\α}}_s+D_u{\mathrm{\α}}_u+D_c{\mathrm{\α}}_c+{\mathrm{\β}}_w+D_y{\mathrm{\β}}_y+D_s{\mathrm{\β}}_s+D_c{\mathrm{\β}}_c)}$

式中：R(λ)为遥感数据计算出的出水反射率。天顶角θ已知，所有的吸收率和散射率均可通过测试获得，黄色物质浓度D_y具有区域稳定性可通过测量得到；未知数只有污染浓度\泥沙浑浊度和叶绿素浓度D_u、D_s和D_c，将三个波段的遥感数据代入，即可得到一个三元一次方程组，解之可求污染物综合浓度D_u、D_s和D_c。

3.2水深的计算

对于清洁水体，或和杂质含量已知，可根据(13)式计算水深H：

$H=\frac{1}{\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})}\ln \left[\frac{R_{\mathrm{wo}}\left(\mathrm{\λ}\right)-\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})}}{(\frac{1}{\cos \mathrm{\θ}}{R}_b-\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})})}\right]$

4.应用实例

通过对该方法进行大量应用试验，包括对珠江三角洲进行多时相水污染遥感监测和广东省全省大中型水库的水质遥感监测。结果表明该方法具有较高的精度，基本达到实用化的要求。以下是利用该方法在珠江三角洲进行多时相水污染遥感监测应用示范的部分结果。

采用1998、2000、2002、2005、2006和2007六个时相的TM(或中巴)数据进行主要水质参数的遥感监测。用地面同步实地水样分析数据对反演结果进行对比，检验模型的精度；对于污染源的遥感，进行了实地考查和取样分析。

遥感结果图如图6至图9所示，卫星同步监测结果验证结果如图10至图13所示。应用结果表明，混浊度和COD的遥感反演结果与水面同步测量数据的相关系数均在0.7以上，具有较高的精度，基本达到实用化的要求。结果图像可以很好地反演区域水质的空间分布情况，水质污染的分布特征与所掌握的情况相符。可以很好的反映跨境污染的情况。可以较好地反映大、中型污染源情况，包括污染源位置、污染的规模及其影响范围。对污染的控制与治理可提供极有用的信息。

Claims (6)

1.一种自动提取水体污染信息的方法，其特征在于，包括以下步骤：

(1)计算获得被测水域的出水反射率；

(2)建立水质遥感模型；

(3)利用水体光谱曲线，根据上述水质遥感模型反演得到水质参数。

2.根据权利要求1所述自动提取水体污染信息的方法，其特征在于，所述水质参数包括清洁水、悬浮泥沙、综合污染物和叶绿素的散射率和吸收系数。

3.根据权利要求2所述自动提取水体污染信息的方法，其特征在于，在测量时分别选取深度较大的清洁水、混浊水、重污染水和高叶绿素含量水体进行测量，然后消除水面镜面反射的影响，计算出水反射率，再进行散射率和吸收系数的反演。

4.根据权利要求1所述自动提取水体污染信息的方法，其特征在于，所述步骤(1)中的出水反射率R_wo(λ)的计算公式为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{R\left(\mathrm{\λ}\right)-{\mathrm{\ρ}}_{\mathrm{ra}}}{[T_a(\mathrm{\θ},\mathrm{\λ})+{\mathrm{\ω}}_{\mathrm{ad}}^{\′}(\mathrm{\θ},\mathrm{\λ})/\cos \mathrm{\θ}]T_a(0,\mathrm{\λ})T_{\mathrm{ws}}}-R_{\mathrm{TM}5}$

式中：R(λ)为遥感数据转换得以到的行星反射率，ρ_ra为大气上行散射反射率，T_a(θ，λ)为大气在入射方向上的透过率，ω′_ad(θ，λ)为下行散射光的半球散射系数，θ为太阳天顶角，T_ws为水面透过率，R_TM5为TM₅计算的反射率。

5.根据权利要求1所述自动提取水体污染信息的方法，其特征在于，所述步骤(2)中的水质遥感模型为：

$R_{\mathrm{wo}}\left(\mathrm{\λ}\right)=\frac{\mathrm{\ω}\left(\mathrm{\Θ}\right)}{4\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})}(1-e^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H})+\frac{1}{\cos \mathrm{\θ}}{R}_b{e}^{-\mathrm{\μ}(\mathrm{\α}+\mathrm{\β})H}$

式中：R_wo(λ)为出水反射率，ω(Θ)为水层垂向上的散射比，Θ为水层垂向上的散射角，θ为水体中的光入射角，α为水体吸收系数，β为水体散射系数，R_b为水底底质反射率，H为水深， $\mathrm{\μ}=1+\frac{1}{\cos \mathrm{\θ}}.$

6.根据权利要求1所述自动提取水体污染信息的方法，其特征在于，所述步骤(3)中的水体光谱曲线是用ASD连续光谱仪在水面测量获得。

Images