[发明专利]一种多光谱遥感卫星影像的去雾方法

权利要求书

1.一种多光谱遥感卫星影像的去雾方法，其特征在于，包括如下步骤：

步骤1、对多光谱遥感影像数据进行大气校正；

步骤2、计算基准雾透射率图t(x)；

步骤3、根据计算的所述基准雾透射率图t(x)，计算遥感影像每个位置处对应的γ(x)；

步骤4、计算全局大气光值A；

步骤5、求取不同波段对应的波长值λj；

步骤6、对多光谱每个通道波段进行去雾处理；

步骤7、对去雾后的每个所述波段进行气溶胶补偿处理；

其中，步骤3中，计算所述遥感影像每个所述位置处对应的γ(x)的公式为：

。

2.如权利要求1所述的一种多光谱遥感卫星影像的去雾方法，其特征在于，所述步骤2中，所述基准雾透射率图的计算方法，包括：

步骤a1、采用大气校正后多光谱数据波段1(band1)和波段2(band2)数据，外推出一个新的波段数据，计算公式为：

band_extrapol(x)＝2*band1(x)-0.95*band2(x)

步骤a2、根据外推出的所述波段band_extrapol，估计雾厚度图HTM，对外推出的波段以每个像元为中心构建一个局部窗口(w×w)，并通过搜索这个局部块内的最小值像元，即黑暗像元来构建整张影像的所述雾厚度图HTM(x)，计算公式为：

其中Ω(x)表示w×w的局部块，x为该局部块的中心位置，局部块尺寸w为常数，

步骤a3、根据获得的所述雾厚度图HTM(x)计算基准雾透射率图t(x)，计算公式为：

t(x)＝1-ωHTM(x)

其中，ω为一个调节系数，ω的取值在1.38和1.8两个数值之间。

3.如权利要求1所述的一种多光谱遥感卫星影像的去雾方法，其特征在于，所述步骤4中，计算所述全局大气光值A的方法，包括：

估计影像上前0.01％最亮像素值得平均值作为所述全局大气光值A；

为了防止存在较多高亮度像元和整体像素值较低的情况对所述全局大气光值A影响，设置了如下式的约束：

使得所述全局大气光值A取值范围在0.45～0.75之间。

4.如权利要求1所述的一种多光谱遥感卫星影像的去雾方法，其特征在于，所述步骤5中，求取不同所述波段对应的所述波长值λj，包括：

取每个所述通道对应的所述波长值范围的中间值。

5.如权利要求1所述的一种多光谱遥感卫星影像的去雾方法，其特征在于，所述步骤6中，各个所述通道波段去雾计算公式为：

其中I_j表示多光谱影像的第j波段影像，J_j表示该波段去雾后的结果，A是全局大气光值，λ_j为对应第j波段波长值，λ₁为第1个波段波长值，t(x)为基准雾透射率图，γ(x)表征大气散射模型，是雾厚度条件相关的物理量。

6.如权利要求1所述的一种多光谱遥感卫星影像的去雾方法，其特征在于，所述步骤7中，每个所述波段进行气溶胶补偿处理的公式为：

其中，j表示多光谱卫星影像的波段序号，DN_j表示去雾后的波段反射率值，DN’_j表示气溶胶补偿后的波段反射率值，表示原始未去雾波段中无雾区域的平均反射率值，表示无雾区域经去雾处理后平均反射率值，abs()表示取绝对值运算。