[发明专利]一种无效像元在轨检测方法

权利要求书

1.一种无效像元在轨检测方法，其特征在于，包括：

从自然目标中选择一个或多个作为类均匀源；

在不同太阳高度角和不同积分时间下，分别对选定的一个或多个类均匀源所在的区域进行多次成像，得到多张遥感图像；

对所述多张遥感图像中满足设定条件的遥感图像进行均一化处理，得到Num景均一化遥感图像；

将所述Num景均一化遥感图像转换为矩阵向量表达式，根据所述矩阵向量表达式，判定得到实测无效像元；

根据所述实测无效像元的位置和参考无效像元的位置，得到可变积分时间面阵CMOS在轨期间的无效像元位置集；

其中，所述将所述Num景均一化遥感图像转换为矩阵向量表达式，根据所述矩阵向量表达式，判定得到实测无效像元，包括：

将均一化处理后的Num＝θ×t×m×n景均一化遥感图像整理为矩阵向量表达式，通过Kolomogorov-Smirnov检验，判定得到无效像元；

其中，将均一化处理后的Num＝θ×t×m×n景均一化遥感图像整理为矩阵向量表达式，通过Kolomogorov-Smirnov检验，判定得到无效像元，具体流程如下：

将所述Num景均一化遥感图像的相同像元(i,j)，逐像元整理成Num维的矩阵

将所述Num景均一化遥感图像的均值，整理成Num维的矩阵

采用矩阵向量表达式，对以Num景均一化遥感图像逐像元为样本的数据和以Num景均一化遥感图像为样本的数据，分别按值x进行累计频率函数计算，分别记为和其中，为所有Num景均一化遥感图像的最大值；

将所述Num景均一化遥感图像逐像元样本数据的累计频率 函数分布与Num景均一化遥感图像均值样本数据的累计频率 函数分布按照如下检验公式进行比较：

若D大于临界值D(Num,α)，则判定该像元为无效像元；其中，D(Num,α)表示：显著水平为α且样本容量为Num对应的与的差距D。

2.根据权利要求1所述的方法，其特征在于，作为类均匀源的自然目标满足如下条件：

自然目标的面积大于面阵CMOS相机的幅宽；

自然目标的均匀性优于3％。

3.根据权利要求1所述的方法，其特征在于，所述在不同太阳高度角和不同积分时间下，分别对选定的一个或多个类均匀源所在的区域进行多次成像，得到多张遥感图像，包括：

根据遥感卫星的轨道预报，通过面阵CMOS相机对选定的一个或多个类均匀源所在的区域，在不同太阳高度角和不同积分时间下进行多次成像，得到多张遥感图像。

4.根据权利要求1所述的方法，其特征在于，所述对所述多张遥感图像中满足设定条件的遥感图像进行均一化处理，得到Num景均一化遥感图像，包括：

从所述多张遥感图像中筛选得到满足设定条件的待处理遥感图像；

对所述待处理遥感图像进行均一化处理，得到所述Num景均一化遥感图像。

5.根据权利要求4所述的方法，其特征在于，所述从所述多张遥感图像中筛选得到满足设定条件的待处理遥感图像，包括：

从所述多张遥感图像中选择云量覆盖小于10％或云量覆盖大于90％的遥感图像，作为所述待处理遥感图像。

6.根据权利要求4所述的方法，其特征在于，所述对所述待处理遥感图像进行均一化处理，得到所述Num景均一化遥感图像，包括：

根据如下公式，对所述待处理遥感图像进行均一化处理，得到所述Num景均一化遥感图像：

其中，为均一化处理后的Num景均一化遥感图像的数字值；DN(θ,t；m,n)为所述类均匀源中的第m个类均匀源、第θ个太阳高度角和第t个积分时间下的第n次遥感图像的数字值；μ_{DN(θ,t；m,n)}为遥感图像的均值；m、θ、t和n均大于1。