[發(fā)明專利]基于PAR數(shù)據(jù)的SIF異常檢測方法、裝置、系統(tǒng)、終端及介質

有效

• 申請?zhí)枺?/td>: 202211441181.2
• 申請日: 2022-11-17
• 公開（公告）號: CN115753712B
• 公開（公告）日: 2023-05-30
• 發(fā)明（設計）人: 姚永軍;李穎池;金磊
• 申請（專利權）人: 北京華益瑞科技有限公司
• 主分類號: G01N21/64
• 分類號: G01N21/64;G01N21/01;G06F17/18
• 代理公司: 北京市盛峰律師事務所 11337
• 代理人: 于國富

本發(fā)明公開一種基于PAR數(shù)據(jù)的SIF異常檢測方法、裝置、系統(tǒng)、終端及介質，涉及遙感數(shù)據(jù)處理領域。包括：獲取光譜數(shù)據(jù)；基于所述光譜數(shù)據(jù)，構建包括至少一個檢測時間單元的檢測時間組，所述檢測時間單元包括：一上行光譜數(shù)據(jù)采集時刻t、與所述采集時刻t相鄰的兩個下行光譜數(shù)據(jù)采集時刻t?1和t+1；基于預先獲取的PAR數(shù)據(jù)，獲取并計算任意一個檢測時間單元中的t時刻的PAR分別與兩個相鄰時刻的PAR的差值絕對值PARsubgt;1/subgt;和PARsubgt;2/subgt;，當PARsubgt;1/subgt;和PARsubgt;2/subgt;均大于等于預設PARsubgt;閾值/subgt;時，則判定基于該檢測時間單元采集時刻的對應光譜數(shù)據(jù)計算得到的SIF為異常值并標記。本發(fā)明通過對PAR監(jiān)測儀監(jiān)測的PAR數(shù)據(jù)結果的比較，校驗SIF異常與否，剔除SIF異常值，實現(xiàn)減少整體SIF結果誤差。

技術領域

本發(fā)明涉及遙感數(shù)據(jù)處理領域，尤其涉及一種基于PAR數(shù)據(jù)的SIF異常檢測方法、裝置、系統(tǒng)、終端及介質。

背景技術

日光誘導葉綠素熒光(Solar-Induced?Chlorophyll?Fluorescence/Sun-InducedChlorophyll?Fluorescence，SIF)是植物在自然光照射下產(chǎn)生的一種熒光，與光合作用高度相關，是探測光合作用的探針。目前，研究者對SIF的探測已廣泛應用于對光合作用的追蹤以及植被脅迫監(jiān)測。

在SIF的相關研究中，高光譜測定上、下行輻射數(shù)據(jù)質量是反演植被SIF的關鍵，植被SIF亮度只占植被反射輻射(即上行輻射)輻亮度的5％左右，因此，光譜數(shù)據(jù)的一些細微偏差將大幅度影響SIF計算準確率。在現(xiàn)有SIF觀測計算方法中，基于雙光譜儀的SIF觀測計算方法因整體不確定性大、成本較高，因此在大多數(shù)情況下并不適用。而采用單光譜儀交替的SIF觀測計算方法則是較為實用的一種觀測計算方式。但是在陰雨、多云等光線條件快速變化的情況下，采用單光譜儀獲取的某時次光譜數(shù)據(jù)會出現(xiàn)巨大偏差，在基于包含巨大偏差的光譜數(shù)據(jù)的基礎上所獲取SIF計算值準確率低。

綜上，目前針對單光譜儀交替測量的高光譜數(shù)據(jù)還沒有檢測SIF異常值的有效方法，因此，亟需找到一種剔除存有偏差光譜數(shù)據(jù)的方法以提高SIF計算的準確率。

發(fā)明內容

本發(fā)明的目的在于提供一種基于PAR數(shù)據(jù)的SIF異常檢測方法、裝置、系統(tǒng)、終端及介質，以解決因存有巨大偏差光譜數(shù)據(jù)而導致SIF計算準確率低的問題。

本發(fā)明第一方面提供一種基于PAR數(shù)據(jù)的SIF異常檢測方法，應用于中央控制單元，包括：

獲取光譜數(shù)據(jù)；

基于所述光譜數(shù)據(jù)，構建包括至少一個檢測時間單元的檢測時間組，所述檢測時間單元包括：一上行光譜數(shù)據(jù)采集時刻t、與所述采集時刻t相鄰的兩個下行光譜數(shù)據(jù)采集時刻t-1和t+1；

基于預先獲取的PAR數(shù)據(jù)和任意一個檢測時間單元記錄的采集時刻，獲取并計算所述任意一個檢測時間單元中各個采集時刻的PAR數(shù)據(jù)，并計算差值絕對值PAR₁和PAR₂，其中，PAR_t表示采集時刻t的PAR數(shù)據(jù)；PAR_t-1表示采集時刻t-1的PAR數(shù)據(jù)；PAR_t+1表示采集時刻t+1的PAR數(shù)據(jù)；

PAR₁＝|PAR_t-PAR_t-1|????????????????????????????(1)

PAR₂＝|PAR_t-PAR_t+1|????????????????????????????(2)