[发明专利]一种多角度偏振卫星的数据预处理方法

权利要求书

1.一种多角度偏振卫星的数据预处理方法，其特征在于，包括以下步骤：

1）读取整轨数据中的观测角度个数信息，判断有效的总行数和总列数，将获取的整轨卫星数据进行分割得到初始分块；

2）判断所述初始分块中是否存在至少一列完整的数据列；如果存在，作为待处理数据块，执行步骤4）

3）将初始分块按行进行分割，分割为若干存在至少一列完整的数据列的子块；所述子块作为待处理数据块；

4）提取所述待处理数据块中的观测数据；

5）从所述完整的数据列中一个端点为起点向另一端点依次进行遍历，当前点作为参考像元，下一个点作为待处理像元，构成参考像元-待处理像元对；所有参考像元-待处理像元对组成参考列邻近像元对数据集；

6）以所述完整的数据列中每个点为起点，向左依次进行遍历，当前点作为参考像元，下一个点作为待处理像元，构成参考像元-待处理像元对；所有参考像元-待处理像元对组成左邻近像元对数据集；向右依次遍历，采用相同方法组成右邻近像元对数据集；

7）将参考像元所有观测角度的数据与待处理像元当前观测角度的数据采用公式（1）和公式（2）计算，判断构成的每个所述参考像元-待处理像元对是否匹配；

δz=VZA_ref,i-VZA_test （1）；

δa=VAA_ref,i-VAA_test （2）；

其中，δz为观测天顶角差值，i为观测角度序号，VZA_ref,i为参考像元的第i个观测天顶角，VZA_test为待处理像元当前观测角度的观测天顶角，δa为观测方位角差值，VAA_ref,i为参考像元的第i个观测方位角，VAA_test为待处理像元当前观测角度的观测方位角；

δz小于观测天顶角阈值Tz且δa小于观测方位角Ta；则判断为匹配，进入步骤9）；

否则，判断为不匹配；进入步骤8）

8）选取待处理像元的下一个观测角度的数据与参考像元通过步骤7）进行匹配，直至匹配；

9）该待处理像元作为参考像元，匹配的观测角度作为参考像元的起始观测角度；

10）采用步骤7）相同的方法处理所述参考列邻近像元对数据集、所述左邻近像元对数据集和所述右邻近像元对数据集，直至数据集中所有的参考像元-待处理像元对完成匹配；

11）将经匹配处理的所有参考像元的各观测角度对应的观测数据代替原数据，完成多角度定序重组。

2.如权利要求1所述的多角度偏振卫星的数据预处理方法，其特征在于，所述Tz的范围为0.5-1.0°；Ta的范围为1.0-2.0°。

3.如权利要求1所述的多角度偏振卫星的数据预处理方法，其特征在于，所述观测数据包括地理坐标域信息和各波段、各角度的数据域信息。

4.如权利要求3所述的多角度偏振卫星的数据预处理方法，其特征在于，所述地理坐标域信息包括地表高程、海陆掩膜、经纬度、格网行列号。

5.如权利要求3所述的多角度偏振卫星的数据预处理方法，其特征在于，所述数据域信息包括辐射亮度、偏振Stokes分量、观测几何。

6.如权利要求1所述的多角度偏振卫星的数据预处理方法，其特征在于，还包括以下步骤：

12) 根据待处理像元的起始观测角度序号AnS和有效观测角度数N，计算观测角度移动次序步长：

将起始观测角度序号AnS及以后的数据与起始观测角度序号AnS以前的数据换位,对数据进行重组，计算参考像元和待处理像元各观测角度下VZA和VAA的总体差值，判断是否所有观测角度均已满足匹配：

其中，i表示观测角度序号；所述观测角度共N个，下标ref和test分别表示参考像元和待处理像元，Tz为观测天顶角阈值和Ta为观测方位角阈值；

如果满足，则数据校核成功；如果不满足，则对新的数据进行步骤1）至步骤12）的处理；发现并剔除异常数据所在的观测角度。