气溶胶反演输入转化错误_气溶胶反演

气溶胶反演方法

利用环境小卫星多光谱数据反演：

方法一：

1.

利用

SPSS

计算出大气光学厚度与大气参数(

ρ

0

(大气的路径辐射项等效反

射率)

、

S

(大气下界的半球反射率)

、

T(

μ

s

)T(

μ

v

)(

大气上行下行透过率

)

)的对

应关系，

据此建立查找表，

然后利用多波段数据进行地气解耦，

得到大气光学厚

度。

2.

构建查询表

利用

6S

模型构建气溶胶光学厚度查询表，输入参数为：太阳

天顶角，气溶胶模式，

550nm

波长处气溶胶光学厚度的等级，查找表计算的波段

(第一和第三波段)

，海拔高度。

3.

数据预处理

(

1

)重采样，为了加快运算速度和提高信噪比；

(

2

)辐射定

标，将图像的

DN

值转化为表观反射率。

4.

结果反演

根据获得的表观反射率计算出

NDVI

(用于识别暗目标)

，利用获

得的太阳高度角对查找表进行插值，

得到要计算波段的不同大气光学厚度下的大

气参数：

ρ

0

、

S

、

T(

μ

s

)T(

μ

v

)

。

5.

图像平滑与成图输出

在获得大气光学厚度后，对结果图像进行平滑处理，

达到内插部分非暗目标点的监测值并抑制异常点的目的，采用

9

×

9

像元的距离

加权平均的滤波方法进行；将结果导入

ArcMap

中，进行叠加矢量图，分等定级

以及添加图名图例等操作，制成专题图。

方法二：

1.

对要反演气溶胶光学厚度的卫星图像惊醒地理和辐射率校正

2.

然后用

MODTRAN

模型模拟生成

τ

(气溶胶光学厚度)和

ρ

(地表反射率)的

查找表

3.

接着判断卫星观测到的地表像元反射率

L

obs

与

MODTRAN

模拟的大气总辐射

L

total

是否相等。

4.

如果不等，就改变

ρ

，再用

MODTRAN

重新计算

L

total

，再判断是否相等。

5.

如果相等，则根据

ρ

和

τ

的关系曲线，由反演出的地表反射率

ρ

map

，计算到

气溶胶光学厚度分布

τ

map

。

利用环境小卫星高光谱数据反演：

方法：

1.

选择用于反演的波段

2.

假设利用某种气溶胶模式条件下，计算红波段和近红外波段表观反射率，不

考虑临近效应影响。大气散射在可见—近红外波段影响是比较大的，在可见

波段影响最大，其次是近红外波段，在中波红外接近于零，因此，在利用近

红外波段反演气溶胶光学特性之前，可以在清洁大气的假设下利用

6S

或者

MODTRAN

辐射传输模型进行大气校正。

3.

利用

NDVI

识别出浓密植被作为暗像元。

4.

根据红波段反射率与近红外波段反射率的对应关系，确定红波段的地表反射

率。