[发明专利]卫星遥感图像的晴空数据集处理方法

权利要求书

1.一种卫星遥感图像的晴空数据集处理方法，其特征在于，包括：

获取遥感卫星的观测数据；

对所述观测数据依次进行大气反射率校正处理、图像分块投影处理及当日数据合成处理，生成彩色卫星遥感中间图像；以及

对所述彩色卫星遥感中间图像进行去除云污染处理，生成对应下垫面的晴空数据集。

2.根据权利要求1所述的方法，其特征在于，生成对应下垫面的晴空数据集包括：基于数字高程模型对遥感栅格数据进行渲染，生成模拟山体阴影的立体图像。

3.根据权利要求1所述的方法，其特征在于，所述去除云污染处理包括：基于最小反射率法和最大归一化植被指数合成法，对所述彩色卫星遥感中间图像进行晴空合成去云处理；其中，所述最小反射率法确定所述彩色卫星遥感中间图像各像元反射率的最小值为合成后的晴空反射率，所述最大归一化植被指数合成法确定所述彩色卫星遥感中间图像各像元归一化植被指数的最大值为合成后的晴空反射率。

4.根据权利要求3所述的方法，其特征在于，所述去除云污染处理还包括：根据不同下垫面的光谱特性，基于西格玛合成法去除所述彩色卫星遥感中间图像的云阴影；所述西格玛合成法根据所述彩色卫星遥感中间图像各像元反射率值的标准差及预设的反射率阈值，确定并删除所述彩色卫星遥感中间图像的云阴影像元，其中所述反射率阈值根据下垫面类型进行预先设置。

5.根据权利要求1-4任一项所述的方法，其特征在于，所述大气反射率校正处理包括根据下述公式确定校正后的大气反射率：

ρ_TOA(μ_s，μ_v，Φ)＝ρ₀(μ_s，μ_v，Φ)+T(μ_s)T(μ_v)ρ_s(μ_s，μ_v，Φ)/[1-ρ_s(μ_s，μ_v，Φ)S]

其中，μ_s、μ_v和Φ为球坐标系下的观测几何参数，T(μ_s)和T(μ_v)为与所述观测几何参数相关的大气通过率，ρ_TOA(μ_s，μ_v，Φ)为校正后的大气反射率，ρ₀(μ_s，μ_v，Φ)为路径辐射反射率，ρ_s(μ_s，μ_v，Φ)为大气表观反射率，S为大气反照率。

6.根据权利要求1-4任一项所述的方法，其特征在于，所述图像分块投影处理包括：全球等经纬度投影处理和/或全球极射赤面投影处理。

7.根据权利要求1-4任一项所述的方法，其特征在于，所述当日数据合成处理包括：根据最优太阳角度，对经所述图像分块投影处理生成的彩色卫星遥感分块图像进行当日多轨道合成。

8.一种卫星遥感图像的晴空数据集处理装置，其特征在于，包括：

数据获取模块，用于获取遥感卫星的观测数据；

图像生成模块，用于对所述观测数据依次进行大气反射率校正处理、图像分块投影处理及当日数据合成处理，生成彩色卫星遥感中间图像；以及

图像处理模块，用于对所述彩色卫星遥感中间图像进行去除云污染处理，生成对应下垫面的晴空数据集。

9.一种电子设备，包括：存储器、处理器及存储在所述存储器中并可在所述处理器中运行的可执行指令，其特征在于，所述处理器执行所述可执行指令时实现如权利要求1-7任一项所述的方法。

10.一种计算机可读存储介质，其上存储有计算机可执行指令，其特征在于，所述可执行指令被处理器执行时实现如权利要求1-7任一项所述的方法。