[發(fā)明專(zhuān)利]一種河口渾濁帶識(shí)別方法和識(shí)別系統(tǒng)

有效

• 申請(qǐng)?zhí)枺?/td>: 202011466246.X
• 申請(qǐng)日: 2020-12-14
• 公開(kāi)（公告）號(hào): CN112763425B
• 公開(kāi)（公告）日: 2021-12-28
• 發(fā)明（設(shè)計(jì)）人: 王重洋;李丹;鄭瓊;姜浩;陳金月;彭詠石;王力;劉楊曉月;陳水森;周霞
• 申請(qǐng)（專(zhuān)利權(quán)）人: 廣東省科學(xué)院廣州地理研究所
• 主分類(lèi)號(hào): G01N21/25
• 分類(lèi)號(hào): G01N21/25;G06K9/00
• 代理公司: 廣州駿思知識(shí)產(chǎn)權(quán)代理有限公司 44425
• 代理人: 潘桂生

1.一種河口渾濁帶識(shí)別方法，其特征在于，包括以下步驟：

通過(guò)衛(wèi)星獲取河口的遙感影像；

對(duì)遙感影像進(jìn)行預(yù)處理；

根據(jù)預(yù)處理后的所述遙感影像計(jì)算河口的總懸浮物含量；

根據(jù)預(yù)處理后的所述遙感影像計(jì)算河口的葉綠素含量；

根據(jù)所述河口的總懸浮物含量和河口的葉綠素含量通過(guò)以下方式計(jì)算河口渾濁帶系數(shù)：

TMZI＝[Log(TSS)-Log(Chla)]/[Log(TSS)+Log(Chla)]；

其中，TMZI表示河口渾濁帶系數(shù)，TSS表示河口的總懸浮物含量，Chla表示河口的葉綠素含量；

將所述河口渾濁帶系數(shù)與預(yù)設(shè)的渾濁帶閾值比對(duì)，將所述遙感影像中所述河口渾濁帶系數(shù)大于所述渾濁帶閾值的范圍確定為河口最大渾濁帶范圍；

其中，所述根據(jù)預(yù)處理后的所述遙感影像計(jì)算河口的總懸浮物含量，包括以下步驟：

從地表反射率提取紅光波段反射率和近紅外波段反射率，根據(jù)所述紅光波段反射率和近紅外波段反射率通過(guò)以下方式計(jì)算河口的總懸浮物含量：

若R₁≥M，

若R₁＜M，

其中，R₁為紅光波段反射率，R₂為近紅外波段反射率，M為預(yù)設(shè)的紅外光波段反射率閾值，a、b和c為預(yù)設(shè)的總懸浮物含量計(jì)算的參量；

所述根據(jù)預(yù)處理后的所述遙感影像計(jì)算河口的葉綠素含量，包括以下步驟：

從所述地表反射率提取藍(lán)光波段反射率、綠光波段反射率和紅光波段反射率，根據(jù)所述藍(lán)光波段反射率、綠光波段反射率和紅光波段反射率通過(guò)以下方式計(jì)算河口的葉綠素含量：

其中，B為藍(lán)光波段反射率，G為綠光波段反射率，R₁為紅光波段反射率；x和y為預(yù)設(shè)的葉綠素含量計(jì)算的參量；e為自然常數(shù)。

2.根據(jù)權(quán)利要求1所述的河口渾濁帶識(shí)別方法，其特征在于，所述將所述河口渾濁帶系數(shù)與預(yù)設(shè)的渾濁帶閾值比對(duì)，將所述遙感影像中所述河口渾濁帶系數(shù)大于所述渾濁帶閾值的范圍確定為河口最大渾濁帶范圍后，還包括以下步驟：

顯示所述遙感影像，并在所述遙感影像中標(biāo)示所述河口最大渾濁帶范圍。

3.根據(jù)權(quán)利要求1所述的河口渾濁帶識(shí)別方法，其特征在于，所述對(duì)遙感影像進(jìn)行預(yù)處理包括以下步驟：

對(duì)遙感影像進(jìn)行輻射定標(biāo)，獲取輻射亮度；

對(duì)輻射定標(biāo)后的遙感影像進(jìn)行大氣校正，獲取包含藍(lán)光波段反射率、綠光波段反射率、紅光波段反射率和近紅外波段反射率的地表反射率。

4.根據(jù)權(quán)利要求3所述的河口渾濁帶識(shí)別方法，其特征在于，所述對(duì)遙感影像進(jìn)行輻射定標(biāo)，獲取輻射亮度，包括以下步驟：

獲取衛(wèi)星的載荷觀測(cè)值、定標(biāo)斜率和定標(biāo)截距；

將衛(wèi)星的載荷觀測(cè)值通過(guò)以下方式轉(zhuǎn)出成輻射亮度：L＝Gain*DN+Bias；其中，L為輻射亮度，DN為衛(wèi)星的載荷觀測(cè)值，Gain為定標(biāo)斜率，Bias為定標(biāo)截距。

5.根據(jù)權(quán)利要求4所述的河口渾濁帶識(shí)別方法，其特征在于，所述對(duì)輻射定標(biāo)后的遙感影像進(jìn)行大氣校正，獲取地表反射率，包括以下步驟：

獲取所述輻射亮度，通過(guò)以下方式將所述輻射亮度轉(zhuǎn)換為地表反射率：

ρs＝Y(jié)/(1+Xc*Y)；

Y＝Xa*L-Xb；

其中，ρs是地表反射率，L是輻射亮度；Xa、Xb、Xc由大氣校正模型計(jì)算得到的用于將所述輻射亮度轉(zhuǎn)換為地表反射率的轉(zhuǎn)換參數(shù)。

6.一種計(jì)算機(jī)設(shè)備，其特征在于：包括儲(chǔ)存器、處理器以及儲(chǔ)存在所述儲(chǔ)存器中并可被所述處理器執(zhí)行的計(jì)算機(jī)程序，所述處理器執(zhí)行所述計(jì)算機(jī)程序時(shí)實(shí)現(xiàn)如權(quán)利要求1-5任意一項(xiàng)所述的方法的步驟。