PDF《D数据大气校正》

2 5

0 0野射率可表示为pT o A( 肛. , 肛. , 妒) =p. ( 肛. , p. , 妒) +丁( 肛. ) 丁( , . . ) p. ( ps , p. , 9) /[

1 一p. ( 肛. , 肛. , 妒) 5] ( 1) 式中,p,,

..是大气层顶表观反射率;

p . 是大气的为大气下界的半球反射率;

丁( Ⅳ. ) 和r(Ⅳ. ) 分别为路径辐射项等效反射率;

p .为地表二向反射率;

Ⅳ. 向下 和向上整层大气透过率( 直射和漫射) . = c o s 日, , 肛. =c o s p . , 臼 .和 臼 .分 别为观测天顶角与太.S、r(, ..)和r(Ⅳ..)取决于单次散射比甜. 、 气溶胶光阳天顶角;

9 为太阳与传感器之间的相对方位角;

S 学厚度丁和 气 溶胶 散射相函数 尸.将式(1)转化为p.(p.,肛. ,

9 ) = [ p T oA( , ∞. , 肛. , 妒) 一Po( 肛.,p, , 妒)] / [ 丁( 肛.) 丁(肛,) + pT o A( . L '

. , 肛,,

9)S ―po( 肛. , 肛.,妒) S ] (2) l S 、 p . 、 丁( ∥. ) 、 丁( ∥. ) 为进行大气校正时所需要个关于变量S、 p . , 和r(Ⅳ. ) r ( 肛. ) 的三元方程组;

解 的 4个 参数.由于丁(弘. ) 和r( 肛. ) 在式中总是以乘 方程组,求出一组大气参数值.重复上述过程,就可积形式出现,因此将丁( 肛. ) r( 肛. ) 作为一个参数考建立大气参数查找表,进而完成遥感影像的大气校虑.故进行大气校正可利用式( 2) , 输入jS、p. 、 丁正( 图1) . ( 肛.) 及丁 ( ∥ . . ) 等参数以及需要进行大气校正的大气I―1 层顶表观反射率pT,.., 经计算可得经过大气校正的l大气层顶表观反射率 p , o A l l 全I.一一.. 一r――2. l

6 tl 地表反射率ps . 上l.;

l根据研究区所 处的地理位置以及图像 的获取时l各种观测几何条件、大气、气溶胶参数 l l I I i I 间, 可以确定应该选择的 大 气气 溶胶模式分别为中F―――――――― ― 参I纬度夏季大 气模 式和大陆型 气溶 胶模 式. 确定了 气1大气榭E 计算I I 翟I 溶胶模式也就确定了 单次散射比∞ . , 和气溶胶散射]:

1 羹l相函 数P . 气溶胶总 含量值则由5

5 0 n m 波长的气l地表反射率l l 表l溶胶光学厚度描述, 分别考虑1

0 个气溶胶光学厚度'

――1值( 丁=0.

1、0.

2、0.

3、0.

4、0.

5、0.

6、0.

7、0.

8、0.

9、 图1 大气校正方法 流程 图1.0) , 7个太阳天顶角( 日.:10.、

200、

300、40.、

500、F'

g.1F10wcha'

'

o '

a t m0

8 p h

8 7

1 0

0 0 '

'

e c '

'

o n me t h o d

6 0

0 、

7 0

0 ) , 7个观测天顶角( 口,,

=10. 、

2 0 . 、

3 00 、

4 0 . 、

5 0

0 、

6 0

0 、

7 0

0 ) ,

1 9个太阳入射与传感器观测方位角3结果与分析差(妒=00 、

1 0

0 、

2 0

0 、

3 0 .、

4 0

0、

5 0

0 、

6 0

0 、

7 0 . 、

8 0 . 、

9 0

0 、

1 0

0 0 、

1 1

00 、

1 2 0. 、

1 3

00 、

1 4

0 0 、

1 5

0 . 、

1 6

0 0 、

1 7

0 0 、

1 8

0 . ) 为直观说明大气校正后的效果,利用上述模型和3个地面反射率(p. =

0 .

0 1 、

0 .

5

0、 0.

9 9) 来模拟对禹城 地区5月份获取的HJ 一1 影像进行大气校计算大气层顶表观反射率p,...将对应的地表反射正.从 图 2中 可以看出,校正后的图像较原图像更率和大气层顶表观反射率代入式(2) , 即可建立一加清晰,对比度 明显增强.图2大气校正前 ( 左) 后( 右) ¨ J 一1 C C I ) 彩色合成 图像 Fi g .