[發明專利]一種可見光遙感圖像反射率獲取方法

有效

• 申請號: 201910370509.8
• 申請日: 2019-05-06
• 公開（公告）號: CN110096668B
• 公開（公告）日: 2023-06-02
• 發明（設計）人: 吳榮華;孫凌;蔣慶;吳鵬;張鵬
• 申請（專利權）人: 國家衛星氣象中心;吳榮華
• 主分類號: G06F17/10
• 分類號: G06F17/10
• 代理公司: 北京力量專利代理事務所(特殊普通合伙) 11504
• 代理人: 王鴻遠

1.一種可見光遙感圖像反射率獲取方法，其特征在于，包括：

接收可見光譜段對地觀測通道的預定天數內全部冷空視場觀測計數值；

利用中位數絕對離差法計算最小冷空視場觀測計數值，該最小冷空視場觀測計數值為暗電流計數值；

根據暗電流計數值計算可見光遙感圖像反射率；

利用中位數絕對離差法計算最小冷空視場觀測計數值，具體包括：

按照時序對預定天數內全部冷空視場觀測計數值進行排序，構建冷空視場觀測計數值時間序列數據集；

選取冷空視場觀測計數值時間序列數據集的中位數；

根據冷空視場觀測計數值時間序列數據集的中位數獲取冷空視場觀測計數值時間序列數據集的絕對離差數據集；

選取絕對離差數據集的中位數；

根據冷空視場觀測計數值時間序列數據集的中位數和絕對離差數據集的中位數獲取最小冷空視場觀測計數值。

2.根據權利要求1所述的可見光遙感圖像反射率獲取方法，其特征在于，通過以下公式獲取最小冷空視場觀測計數值：

DC_est＝DC_median-β*MAD

其中，DC_est為最小冷空視場觀測計數值，DC_median為冷空視場觀測計數值時間序列數據集的中位數，MAD為絕對離差數據集的中位數，β為常數。

3.根據權利要求2所述的可見光遙感圖像反射率獲取方法，其特征在于，通過以下公式根據暗電流計數值計算可見光遙感圖像反射率：

ρ＝k(DN-DC_Est)

其中，ρ為可見光遙感圖像反射率，k為輻射定標系數，DN為觀測計數值，DC_Est為暗電流計數值。

4.根據權利要求3所述的可見光遙感圖像反射率獲取方法，其特征在于，還包括：當計算歷史數據的可見光遙感圖像反射率時，根據同一歷史時間的暗電流計數值計算可見光遙感圖像反射率。

5.根據權利要求4所述的可見光遙感圖像反射率獲取方法，其特征在于，當計算歷史數據的可見光遙感圖像反射率時，通過以下公式根據同一歷史時間的暗電流計數值計算可見光遙感圖像反射率：

ρ_歷史時間＝k(DN_歷史時間-DC_{Est歷史時間})

其中，DC_{Est歷史時間}＝DC_{median歷史時間}-β*MAD_歷史時間，ρ_歷史時間為歷史時間的可見光遙感圖像反射率，k為輻射定標系數，DN_歷史時間為歷史時間的觀測計數值，DC_{Est歷史時間}為歷史時間的暗電流計數值。

6.根據權利要求3所述的可見光遙感圖像反射率獲取方法，其特征在于，還包括：當計算實時的可見光遙感圖像反射率時，根據提前預定天數內的暗電流計數值計算可見光遙感圖像反射率。

7.根據權利要求6所述的可見光遙感圖像反射率獲取方法，其特征在于，還包括：

ρ_實時＝k(DN_實時-DC_Est提前)

其中，DCEst_提前＝DCmedian_提前-β*MAD_提前，ρ_實時為實時的可見光遙感圖像反射率，k為輻射定標系數，DN_實時為實時的觀測計數值，DC_Est提前為提前預定天數內的暗電流計數值。

8.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，該程序被處理器執行時實現如權利要求1-7任一項所述方法的步驟。

9.一種計算機設備，包括存儲器、處理器及存儲在存儲器上并在處理器上運行的計算機程序，其特征在于，所述處理器執行所述程序時實現如權利要求1-7任一項所述方法的步驟。