[发明专利]环境一号HJ-1A/B星CCD数据的云检测方法

权利要求书

1.一种环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，包括如下步骤：

S1：利用MOD09Q1数据建立地表反射率数据库；

S2：对环境一号HJ-1A/B星的CCD数据进行辐射定标的预处理，获取CCD数据各波段数据的表观反射率；

S3：将S2中环境一号HJ-1A/B星的CCD数据各波段数据的表观反射率数据与从S1中地表反射率数据库中选取成像日期内的地表反射率数据进行比较进行云像元识别，并生成第一层云掩膜；

S4：利用环境一号HJ-1A/B星的CCD数据的第3波段的表观反射率数据值域范围进行云像元识别，并生成第二层云掩膜；

S5：利用环境一号HJ-1A/B星的CCD数据的第4波段与第3波段表观反射率数据比值的值域范围进行云像元识别，并生成第三层云掩膜；

S6：根据步骤S3-S5的云像元识别结果，得到环境一号HJ-1A/B星的CCD数据云掩膜。

2.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，S1进一步包括如下步骤：

S1.1：获取一段时间内的MOD09Q1数据；

S1.2：将所述MOD09Q1数据中同时期的多个MOD09Q1影像的第1波段地表反射率数据和第2波段地表反射率数据分别进行拼接，形成一个第1波段地表反射率数据和一个第2波地表反射率段数据，并设置数据输出格式为GeoTiff，一个所述MOD09Q1影像包含多个第1波段地表反射率数据和多个第2波段地表反射率数据；

S1.3：重复S1.2，获取在S1.2中设定的一段时间内的第1波段和第2波段地表反射率数据，并将所述第1波段和所述第2波段地表反射率数据分别组成两个数据库，所述数据库中存储的数据均为GeoTiff格式数据。

3.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，S2进一步包括如下步骤：

S2.1：从环境一号HJ-1A/B星的CCD数据的辅助数据文件中获取定标系数，并将所述的CCD数据的第3波段和4波段的像元亮度DN值转换为表观辐亮度；

S2.2：从环境一号HJ-1A/B星的CCD数据的辅助数据文件中获取太阳天顶角，将S2.1中的表观辐亮度归一化为第3波段和4波段的表观反射率。

4.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，S3进一步包括如下步骤：

S3.1：从环境一号HJ-1A/B星的CCD数据的辅助数据文件中获取成像日期；

S3.2：根据S3.1中的成像日期，从S1中的地表反射率数据库中选择最接近CCD数据成像日期的地表反射率数据；

S3.3：将环境一号HJ-1A/B星CCD数据的第3波段和4波段表观反射率数据分别与S3.2中成像日期内的地表反射率数据的第1波段和2波段地表反射率数据做差，若CCD数据中像元所对应的这两个差值至少有一个大于0.2，则该像元被判定为云像元，赋值为1，其余无云像元赋值为0，生成第一层云掩膜。

5.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，在S4中，若环境一号HJ-1A/B星CCD数据中像元的第3波段表观反射率值大于0.22，则该像元被判定为云像元，赋值为1，其余无云像元赋值为0，生成第二层云掩膜。

6.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，在S5中，若环境一号HJ-1A/B星CCD数据中像元的第4波段表观反射率与第3波段表观反射率的比值大于0.87，且小于1.15，则该像元被判定为云像元，赋值为1，其余无云像元赋值为0，生成第三层云掩膜。

7.如权利要求1所述的环境一号HJ-1A/B星CCD数据的云检测方法，其特征在于，S6进一步包括如下步骤：

S6.1：将S3-S5中所述第一、二、三层云掩膜进行相加求和，生成新的云掩膜；

S6.2：将所述新的云掩膜中，值为1的像元赋值为0，表示无云像元；值为2的像元保留原值不变，表示可能为云的像元；值为3的像元赋值为1，表示云像元。