[發明專利]一種無效像元在軌檢測方法

有效

• 申請號: 201611154200.8
• 申請日: 2016-12-14
• 公開（公告）號: CN106846292B
• 公開（公告）日: 2019-08-27
• 發明（設計）人: 王愛春;傅俏燕;閔祥軍;陸書寧;潘志強;韓啟金;張學文;劉李;邵俊;李照洲;趙航
• 申請（專利權）人: 中國資源衛星應用中心
• 主分類號: G06T7/00
• 分類號: G06T7/00;G06T5/00;G06T5/50
• 代理公司: 中國航天科技專利中心 11009
• 代理人: 范曉毅

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。