[發(fā)明專利]基于海洋二號衛(wèi)星掃描微波輻射計的海冰密集度獲取方法

有效

• 申請?zhí)枺?/td>: 201410498638.2
• 申請日: 2014-09-25
• 公開（公告）號: CN104198052B
• 公開（公告）日: 2017-07-14
• 發(fā)明（設計）人: 石立堅;王其茂;林明森;鄒斌;黃磊
• 申請（專利權）人: 國家衛(wèi)星海洋應用中心
• 主分類號: G01J5/60
• 分類號: G01J5/60
• 代理公司: 北京超凡志成知識產(chǎn)權代理事務所(普通合伙)11371
• 代理人: 吳開磊

1.基于海洋二號衛(wèi)星掃描微波輻射計的海冰密集度獲取方法，其特征在于，包括：

根據(jù)典型特征區(qū)域的第一亮溫數(shù)據(jù)，得到亮溫特征值，其中所述典型特征區(qū)域包括一年冰特征區(qū)域、多年冰特征區(qū)域及海水特征區(qū)域，所述第一亮溫數(shù)據(jù)通過海洋二號衛(wèi)星掃描微波輻射計獲取，其中，根據(jù)所述一年冰特征區(qū)域的18.7GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,FY,18.7V；根據(jù)所述一年冰特征區(qū)域的18.7GHz頻段水平極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,FY,18.7H；根據(jù)所述一年冰特征區(qū)域的37GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,FY,37V；根據(jù)所述多年冰特征區(qū)域的18.7GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,MY,18.7V；根據(jù)所述多年冰特征區(qū)域的18.7GHz頻段水平極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,MY,18.7H；根據(jù)所述多年冰特征區(qū)域的37GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,MY,37V；根據(jù)所述海水特征區(qū)域的18.7GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,OW,18.7V；根據(jù)所述海水特征區(qū)域的18.7GHz頻段水平極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,OW,18.7H；根據(jù)所述海水特征區(qū)域的37GHz頻段垂直極化方向的亮溫數(shù)據(jù)，得到亮溫特征值T_b,OW,37V；

按預先劃分的像元，獲取監(jiān)測區(qū)域中每個所述像元的第二亮溫數(shù)據(jù)，并利用典型特征區(qū)域的亮溫特征值及每個像元的亮溫數(shù)據(jù)分別計算每個所述像元所對應的一年冰密集度及多年冰密集度；

根據(jù)每個像元的所述一年冰密集度及多年冰密集度的計算結果，將所述監(jiān)測區(qū)域劃分為一年冰區(qū)域、多年冰區(qū)域及海水區(qū)域；

利用所述一年冰區(qū)域、所述多年冰區(qū)域及所述海水區(qū)域的預設監(jiān)測時間的亮溫數(shù)據(jù)，得到每個所述監(jiān)測時間所對應的亮溫特征值，作為所述監(jiān)測區(qū)域的動態(tài)亮溫特征值，其中所述第二亮溫數(shù)據(jù)通過海洋二號衛(wèi)星掃描微波輻射計獲取；

根據(jù)所述監(jiān)測區(qū)域的動態(tài)亮溫特征值，獲取所述監(jiān)測區(qū)域的動態(tài)海冰密集度；

其中，利用公式計算得到一年冰密集度C_FY，利用公式計算得到多年冰密集度C_MY，利用公式C_T＝C_FY+C_MY計算得到海冰密集度C_T，其中，D＝D₀+D₁*PR+D₂*GR+D₃*PR*GR，PR為監(jiān)測區(qū)域的極化梯度率，GR為監(jiān)測區(qū)域的光譜梯度率，M_i、F_i、D_i為亮溫系數(shù)，i＝0～3；

其中，利用T_b,FY,18.7V、T_b,FY,18.7H、T_b,FY,37V、T_b,MY,18.7V、T_b,MY,18.7H、T_b,MY,37V、T_b,OW,18.7V、T_b,OW,18.7H及T_b,OW,37V得到亮溫系數(shù)M_i、F_i及D_i，i＝0～3；其中M_i、F_i及D_i的計算方法為：

M₀＝A₄B₀-A₀B₄ F₀＝A₀B₂-A₂B₀ D₀＝A₄B₂-A₂B₄

M₁＝A₅B₀-A₁B₄ F₁＝A₁B₂-A₃B₀ D₁＝A₅B₂-A₃B₄

M₂＝A₄B₁-A₀B₅ F₂＝A₀B₃-A₂B₁ D₂＝A₄B₃-A₂B₅

M₃＝A₅B₁-A₁B₅，F(xiàn)₃＝A₁B₃-A₃B₁，D₃＝A₅B₃-A₃B₅；

A₀＝-T_b,OW,18.7V+T_b,OW,18.7H B₀＝-T_b,OW,37V+T_b,OW,18.7V

A₁＝T_b,OW,18.7V+T_b,OW,18.7H B₁＝T_b,OW,37V+T_b,OW,18.7V

A₂＝T_b,MY,18.7V-T_b,MY,18.7H+A₀ B₂＝T_b,MY,37V-T_b,MY,18.7V+B₀

A₃＝-T_b,MY,18.7V-T_b,MY,18.7H+A₁ B₃＝-T_b,MY,37V-T_b,MY,18.7V+B₁

A₄＝T_b,FY,18.7V-T_b,FY,18.7H+A₀ B₄＝T_b,FY,37V-T_b,FY,18.7V+B₀

A₅＝-T_b,FY,18.7V-T_b,FY,18.7H+A₁，B₅＝-T_b,FY,37V-T_b,FY,18.7V+B₁；

其中，18.7GHz垂直極化方向可簡化表示為18.7V；18.7GHz水平極化方向可簡化表示為18.7H；37GHz垂直極化方向可簡化表示為37V。