本發(fā)明涉及衛(wèi)星應(yīng)用領(lǐng)域，具體而言，涉及海冰標(biāo)識(shí)確定方法及裝置。

背景技術(shù)：

海冰是指直接由海水在海面上凍結(jié)而成的咸水冰，亦包括進(jìn)入海洋中的大陸冰川(冰山和冰島)、河冰及湖冰。海冰對高緯度地區(qū)以至極地地區(qū)的水文、熱力循環(huán)、洋流和生態(tài)系統(tǒng)都有重大影響。海冰過多時(shí)可能會(huì)導(dǎo)致海港封港，堵塞航道，擠壓船舶等問題，它是海洋5種主要災(zāi)害之一(另外4種是海嘯、風(fēng)暴潮、災(zāi)害海浪和赤潮)，素有白色殺手之稱。

準(zhǔn)確地檢測地表類型對反演風(fēng)矢量非常重要，可避免海洋sigma0測量結(jié)果被陸地或海冰的sigma0測量結(jié)果污染，因此有必要對受到海冰影響的觀測結(jié)果進(jìn)行海冰標(biāo)識(shí)。

但，傳統(tǒng)的海冰標(biāo)識(shí)方法過于復(fù)雜，計(jì)算周期過長，不適用于實(shí)效性要求過高的情況。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的在于提供海冰標(biāo)識(shí)確定方法，以提高確定海冰標(biāo)識(shí)的過程的便捷程度。

第一方面，本發(fā)明實(shí)施例提供了海冰標(biāo)識(shí)確定方法，包括：

建立地球表面的等經(jīng)緯度網(wǎng)格；

獲取預(yù)定時(shí)間內(nèi)的所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值；

獲取所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海面溫度數(shù)值；

根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格的海冰標(biāo)識(shí)情況。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第一種可能的實(shí)施方式，其中，所述預(yù)定時(shí)間為一周。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第二種可能的實(shí)施方式，其中，步驟所述根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第一待優(yōu)化海冰網(wǎng)格點(diǎn)；所述第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值，和/或所述第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值；

將所述第一待優(yōu)化海冰網(wǎng)格點(diǎn)和第二待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，所述第二待優(yōu)化海冰網(wǎng)格點(diǎn)與所述第一待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第三種可能的實(shí)施方式，其中，步驟所述根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第三待優(yōu)化海冰網(wǎng)格點(diǎn)；所述第三待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值；

將所述第三待優(yōu)化海冰網(wǎng)格點(diǎn)和第四待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，所述第四待優(yōu)化海冰網(wǎng)格點(diǎn)與所述第三待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離，且所述第四待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第四種可能的實(shí)施方式，其中，步驟所述根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第五待優(yōu)化海冰網(wǎng)格點(diǎn)；所述第五待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值；

將所述第五待優(yōu)化海冰網(wǎng)格點(diǎn)和第六待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，所述第六待優(yōu)化海冰網(wǎng)格點(diǎn)與所述第五待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離，且所述第六待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第五種可能的實(shí)施方式，其中，步驟所述獲取預(yù)定時(shí)間內(nèi)的所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值包括：

獲取一周內(nèi)每個(gè)網(wǎng)格點(diǎn)的海冰密度數(shù)值；

將獲取到的每個(gè)網(wǎng)格點(diǎn)海冰密度數(shù)值的最高值保留，并將其余值刪除。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第六種可能的實(shí)施方式，其中，還包括：

建立海冰標(biāo)識(shí)統(tǒng)計(jì)表；

將全部所述實(shí)際海冰網(wǎng)格點(diǎn)的行列號(hào)均寫入所述海冰標(biāo)識(shí)統(tǒng)計(jì)表中。

結(jié)合第一方面，本發(fā)明實(shí)施例提供了第一方面的第七種可能的實(shí)施方式，其中，還包括：

若獲取到的后向散射系數(shù)觀測結(jié)果，則提取所述后向散射系數(shù)觀測結(jié)果所對應(yīng)的后向散射系數(shù)測量脈沖的第一經(jīng)緯度坐標(biāo)；

計(jì)算所述第一經(jīng)緯度坐標(biāo)在所述等經(jīng)緯度網(wǎng)格中所對應(yīng)的目標(biāo)行列號(hào)；

在所述海冰標(biāo)識(shí)統(tǒng)計(jì)表中查找所述目標(biāo)行列號(hào)所對應(yīng)的網(wǎng)格是否為實(shí)際海冰網(wǎng)格點(diǎn)，若是，則將所述后向散射系數(shù)觀測結(jié)果設(shè)置海冰標(biāo)識(shí)。

第二方面，本發(fā)明實(shí)施例還提供了一種海冰標(biāo)識(shí)確定裝置，包括：

建立模塊，用于建立地球表面的等經(jīng)緯度網(wǎng)格；

第一獲取模塊，用于獲取預(yù)定時(shí)間內(nèi)的所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值；

第二獲取模塊，用于獲取所述等經(jīng)緯度網(wǎng)格的各個(gè)網(wǎng)格點(diǎn)的海面溫度數(shù)值；

確定模塊，用于根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格的海冰標(biāo)識(shí)情況。

結(jié)合第二方面，本發(fā)明實(shí)施例提供了第二方面的第一種可能的實(shí)施方式，其中，所述預(yù)定時(shí)間為一周。

本發(fā)明實(shí)施例提供的海冰標(biāo)識(shí)確定方法，與現(xiàn)有技術(shù)中需要通過相對復(fù)雜的算法來確定海冰標(biāo)識(shí)相比，其通過先建立地球表面的等經(jīng)緯度網(wǎng)格；再分別獲取預(yù)定時(shí)間內(nèi)的所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值；和獲取所述等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海面溫度數(shù)值；最后再根據(jù)所述海冰密集度數(shù)值和所述海面溫度數(shù)值確定所述等經(jīng)緯度網(wǎng)格的海冰標(biāo)識(shí)情況。即，本申請所提供的方案通過這兩種能夠方便獲取到數(shù)據(jù)(海冰密集度和海面溫度數(shù)值)來進(jìn)行海冰標(biāo)識(shí)，雖然降低了一定的準(zhǔn)確度，但保證了效率。

為使本發(fā)明的上述目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉較佳實(shí)施例，并配合所附附圖，作詳細(xì)說明如下。

附圖說明

為了更清楚地說明本發(fā)明實(shí)施例的技術(shù)方案，下面將對實(shí)施例中所需要使用的附圖作簡單地介紹，應(yīng)當(dāng)理解，以下附圖僅示出了本發(fā)明的某些實(shí)施例，因此不應(yīng)被看作是對范圍的限定，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他相關(guān)的附圖。

圖1示出了相關(guān)技術(shù)中，衛(wèi)星進(jìn)行地球掃描的第一種示意圖；

圖2示出了相關(guān)技術(shù)中，衛(wèi)星進(jìn)行地球掃描的第二種示意圖；

圖3示出了本發(fā)明實(shí)施例所提供的海冰標(biāo)識(shí)確定方法的基本流程圖；

圖4示出了本發(fā)明實(shí)施例所提供的海冰標(biāo)識(shí)確定裝置的基本方框圖。

具體實(shí)施方式

下面將結(jié)合本發(fā)明實(shí)施例中附圖，對本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。通常在此處附圖中描述和示出的本發(fā)明實(shí)施例的組件可以以各種不同的配置來布置和設(shè)計(jì)。因此，以下對在附圖中提供的本發(fā)明的實(shí)施例的詳細(xì)描述并非旨在限制要求保護(hù)的本發(fā)明的范圍，而是僅僅表示本發(fā)明的選定實(shí)施例?；诒景l(fā)明的實(shí)施例，本領(lǐng)域技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)的前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。

本申請所提供的方法涉及衛(wèi)星遙感測量，下面現(xiàn)對部分技術(shù)進(jìn)行說明，具體如下。

1，2012年8月發(fā)射的海洋二號(hào)，搭載有我國第一個(gè)可業(yè)務(wù)化運(yùn)行的微波散射計(jì)HY2-SCAT。HY2-SCAT主要用于全球海面風(fēng)場觀測，測風(fēng)風(fēng)速范圍為4～24m/s，風(fēng)速精度為2m/s或10％；風(fēng)向測量范圍為0～360°，風(fēng)向精度為±20°。HY2-SCAT工作頻率為13.256GHz，采用筆形波束圓錐掃描方式，通過筆形波束以固定仰角圍繞天底方向旋轉(zhuǎn)，在衛(wèi)星平臺(tái)順軌方向的運(yùn)動(dòng)中形成一定的地面覆蓋刈幅(如圖1所示)；散射計(jì)系統(tǒng)包括VV和HH兩個(gè)極化方式，分別以不同入射角進(jìn)行觀測，在平臺(tái)的運(yùn)動(dòng)過程中對同一分辨單元可獲取不同極化方式，不同入射角度的多次測量結(jié)果(如圖2所示)，以克服海面風(fēng)場方向反演的多值模糊問題。其中內(nèi)波束采用HH極化方式，入射角為41°，對應(yīng)地面足印大小約為23km×31km，刈幅寬度為1400km。外波束采用VV極化方式，入射角為48°，對應(yīng)地面足印大小約為25km×38km，刈幅寬度為1700km。

2，關(guān)于HY-2散射計(jì)數(shù)據(jù)產(chǎn)品簡介

海洋二號(hào)衛(wèi)星散射計(jì)目前可提供的數(shù)據(jù)產(chǎn)品分為L1B級(jí)產(chǎn)品數(shù)據(jù)產(chǎn)品、L2A級(jí)數(shù)據(jù)產(chǎn)品、L2B級(jí)數(shù)據(jù)產(chǎn)品和L3級(jí)數(shù)據(jù)產(chǎn)品。與本發(fā)明所提供的技術(shù)方案相關(guān)的為L1B級(jí)數(shù)據(jù)和L2A級(jí)數(shù)據(jù)。本發(fā)明所提供的技術(shù)方案所解決的問題是從L1B級(jí)數(shù)據(jù)文件到L2A級(jí)數(shù)據(jù)產(chǎn)品生產(chǎn)過程中的海冰標(biāo)識(shí)問題。

其中L1B數(shù)據(jù)為以遙測幀的時(shí)間為順序進(jìn)行存儲(chǔ)的散射計(jì)觀測數(shù)據(jù)。每個(gè)遙測幀包括96個(gè)散射計(jì)測量脈沖，每個(gè)測量脈沖包括后向散射系數(shù)，每個(gè)脈沖足跡的地理位置、以及用來描述測量數(shù)據(jù)的質(zhì)量和不確定性等信息的參數(shù)。

L2A產(chǎn)品文件包括衛(wèi)星平臺(tái)在一個(gè)空間軌道內(nèi)獲得的每個(gè)雷達(dá)后向散射sigma0測量值。此外，L2A產(chǎn)品也包含一些與每個(gè)sigma0測量值相對應(yīng)的輔助數(shù)據(jù)元素。這些輔助數(shù)據(jù)元素列出了每個(gè)sigma0測量值的位置、質(zhì)量、以及不確定性等相關(guān)信息。L2A產(chǎn)品中的sigma0以風(fēng)矢量單元進(jìn)行分組。每一個(gè)風(fēng)矢量單元行對應(yīng)地面測量刈幅的一個(gè)交軌切割。每一個(gè)L2A風(fēng)矢量單元都是一個(gè)25km的正方形。因此，需要1624個(gè)風(fēng)矢量單元行完成對地球的一次完整覆蓋。

3、SSMI海冰數(shù)據(jù)文件

SSMI海冰輔助數(shù)據(jù)為EUMETSAT Ocean&Sea Ice Satellite Application Facility(OSI)提供的數(shù)據(jù)產(chǎn)品，通過FTP方式發(fā)布，時(shí)效性為2天(也就是說可以獲得兩天之前的數(shù)據(jù))，采用極方位立體投影(Polar stereographic)方式組織數(shù)據(jù)。

4、ECMWF海面溫度預(yù)報(bào)數(shù)據(jù)

ECMWF海面溫度預(yù)報(bào)數(shù)據(jù)為歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts簡稱ECMWF)提供的海面溫度預(yù)報(bào)產(chǎn)品。該產(chǎn)品為網(wǎng)格化的數(shù)值產(chǎn)品，本發(fā)明用到的海溫預(yù)報(bào)數(shù)據(jù)產(chǎn)品空間分辨率為0.25°×0.25°，時(shí)間間隔為6小時(shí)。

5、海冰標(biāo)識(shí)的目的和意義

準(zhǔn)確地檢測地表類型對反演風(fēng)矢量非常重要，可避免海洋sigma0測量結(jié)果被陸地或海冰的sigma0測量結(jié)果污染，因此有必要對受到海冰影響的觀測結(jié)果進(jìn)行海冰標(biāo)識(shí)。

下面對本申請所提供的方法進(jìn)行說明：

如圖3所示，本申請?zhí)峁┝撕１鶚?biāo)識(shí)確定方法，包括：

S101，建立地球表面的等經(jīng)緯度網(wǎng)格；

S102，獲取預(yù)定時(shí)間內(nèi)的等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值；

S103，獲取等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海面溫度數(shù)值；

S104，根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格的海冰標(biāo)識(shí)情況。

步驟S101中，等經(jīng)緯度網(wǎng)格指的是覆蓋在地球表面的分辨率為0.25°×0.25°等經(jīng)緯度網(wǎng)格，第一維度對應(yīng)緯度，第二維度對應(yīng)經(jīng)度。實(shí)際操作中，一般會(huì)選擇將等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)均確定其海冰標(biāo)識(shí)情況，也可以是只針對某一個(gè)區(qū)域進(jìn)行海冰標(biāo)識(shí)的確定。

步驟S102中，需要獲取一定時(shí)間內(nèi)的海冰密集度數(shù)值，海冰密集度數(shù)值通?？梢酝ㄟ^SSMI海冰數(shù)據(jù)產(chǎn)品來獲得，當(dāng)然，還可以通過其他途徑來獲取海冰密集度數(shù)值，如通過AMSE海冰數(shù)據(jù)產(chǎn)品和AMSR海冰數(shù)據(jù)產(chǎn)品等來獲得。SSMI海冰數(shù)據(jù)產(chǎn)品一般是獲取到2天以內(nèi)的數(shù)據(jù)。但海冰數(shù)據(jù)產(chǎn)品在生成的過程中并不必然能夠覆蓋到地表的每個(gè)位置，可能會(huì)造成某一個(gè)產(chǎn)品中，地表的某些點(diǎn)區(qū)域是沒有密集度數(shù)據(jù)的，進(jìn)而，為了保證能夠獲取到每個(gè)位置的海冰密集度數(shù)值，應(yīng)當(dāng)將預(yù)定時(shí)間設(shè)置為一周，也就是連續(xù)獲取一周內(nèi)的SSMI海冰數(shù)據(jù)產(chǎn)品，這可能導(dǎo)致某一個(gè)點(diǎn)位置會(huì)有多個(gè)海冰密集度數(shù)值，為了保證海冰判斷沒有遺漏(保證海冰均能夠被識(shí)別出來)，此處應(yīng)當(dāng)保留各個(gè)密集度數(shù)值的最大值，并將除最大值外的其他值刪除。

步驟S103中，需要獲取海綿溫度數(shù)據(jù)，一般可以選擇ECMWF和NCEP海面溫度數(shù)據(jù)信息。此處獲取的海溫?cái)?shù)據(jù)優(yōu)選為最近時(shí)刻的。

步驟S104，在確定了海冰密集度數(shù)值和海面溫度數(shù)值后，便可以依據(jù)這兩個(gè)數(shù)值來確定海冰標(biāo)識(shí)了，或者說確定步驟S101中的等經(jīng)緯度網(wǎng)格中每個(gè)網(wǎng)格點(diǎn)是否有海冰。

步驟S104在具體執(zhí)行的時(shí)候有4種方式，可以依據(jù)不同的使用情況(如避免遺漏、避免誤判、避免遺漏和避免誤判的均衡)的不同，來選擇適當(dāng)?shù)姆绞健?/p>

具體的，下面先將各種不同的方式進(jìn)行說明：

第一種方式，步驟S104，即步驟根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第一待優(yōu)化海冰網(wǎng)格點(diǎn)；第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值，或第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值；

將第一待優(yōu)化海冰網(wǎng)格點(diǎn)和第二待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，第二待優(yōu)化海冰網(wǎng)格點(diǎn)與第一待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離。

第二種方式，步驟S104，即步驟根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第一待優(yōu)化海冰網(wǎng)格點(diǎn)；第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值，以及第一待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值；

將第一待優(yōu)化海冰網(wǎng)格點(diǎn)和第二待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，第二待優(yōu)化海冰網(wǎng)格點(diǎn)與第一待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離。

第三種方式，步驟S104，即步驟根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第三待優(yōu)化海冰網(wǎng)格點(diǎn)；第三待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值；

將第三待優(yōu)化海冰網(wǎng)格點(diǎn)和第四待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，第四待優(yōu)化海冰網(wǎng)格點(diǎn)與第三待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離，且第四待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值。

第四種方式，步驟S104，即步驟根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格上的海冰標(biāo)識(shí)情況包括：

確定第五待優(yōu)化海冰網(wǎng)格點(diǎn)；第五待優(yōu)化海冰網(wǎng)格點(diǎn)的海面溫度數(shù)值小于預(yù)設(shè)的第二數(shù)值；

將第五待優(yōu)化海冰網(wǎng)格點(diǎn)和第六待優(yōu)化海冰網(wǎng)格點(diǎn)作為實(shí)際海冰網(wǎng)格點(diǎn)，第六待優(yōu)化海冰網(wǎng)格點(diǎn)與第五待優(yōu)化海冰網(wǎng)格點(diǎn)的距離小于預(yù)設(shè)的參考距離，且第六待優(yōu)化海冰網(wǎng)格點(diǎn)的海冰密集度數(shù)值大于預(yù)設(shè)的第一數(shù)值。

上述四種方式相比，均是采用先確定核心點(diǎn)(如第一、三、五待優(yōu)化海冰網(wǎng)格點(diǎn))，再確定外圍點(diǎn)(如第二、四、六待優(yōu)化海冰網(wǎng)格點(diǎn))的方式來確定，也就是，這幾種方式均是采用確定核心點(diǎn)后，采用外延、外擴(kuò)的方式來增加被判斷為實(shí)際海冰網(wǎng)格點(diǎn)的網(wǎng)格點(diǎn)數(shù)量。但，上述四種方式相比，也有相當(dāng)?shù)膮^(qū)別，也就是核心點(diǎn)的確定方式各不相同，并且，外圍點(diǎn)的確定方式也各不相同。

相比之下，第一種方式更側(cè)重于避免遺漏，主要是采用第一種方式的情況下，確定的核心點(diǎn)的數(shù)量是很多的，即只要海冰密集度符合要求，或者是海面溫度符合要求就可以作為核心點(diǎn)，進(jìn)而在核心點(diǎn)數(shù)量較多的情況下，其所確定出的外圍點(diǎn)的數(shù)量也是較多的。

和第一種方式相比，第二種方式比較嚴(yán)謹(jǐn)，主要是其確定出的核心點(diǎn)數(shù)量較少，進(jìn)而最終確定出的實(shí)際海冰網(wǎng)格點(diǎn)較少。

第三種方式和第四種方式均是在判斷策略上進(jìn)行了一定程度的調(diào)整，以適用于不同的使用環(huán)境。上述這四種方式中，預(yù)設(shè)的參考距離為40-60KM，優(yōu)選為50KM，第一數(shù)值可以設(shè)置為0，第二數(shù)值可以設(shè)置為2°。

在以效率為優(yōu)先考慮條件的情況下，應(yīng)當(dāng)使用上述第一種方式。通過該種方式的執(zhí)行，可以確保沒有被標(biāo)識(shí)為海冰的網(wǎng)格單元一定沒有海冰標(biāo)識(shí)，但是會(huì)有部分不是海冰的區(qū)域被標(biāo)記為海冰，相對于該算法提供的海冰標(biāo)記的便利程度(計(jì)算效率提高)，同時(shí)權(quán)衡多剔除一部分?jǐn)?shù)據(jù)付出的代價(jià)，遠(yuǎn)小于由于海冰未能被正確標(biāo)記而產(chǎn)生錯(cuò)誤數(shù)據(jù)的代價(jià)，這種通過多剔除一小部分?jǐn)?shù)據(jù)來保證所有的海冰區(qū)域都能被標(biāo)識(shí)的策略是可接受，并且高效的。

為了便于使用，可以在確定了海冰標(biāo)識(shí)之后，再建立一個(gè)數(shù)據(jù)表進(jìn)行海冰標(biāo)識(shí)的記錄。具體的，本申請所提供的方法還包括：

建立海冰標(biāo)識(shí)統(tǒng)計(jì)表；

將全部實(shí)際海冰網(wǎng)格點(diǎn)的行列號(hào)均寫入海冰標(biāo)識(shí)統(tǒng)計(jì)表中。

需要說明的是，海冰標(biāo)識(shí)統(tǒng)計(jì)表有兩種形式，第一種是該表中只記錄被標(biāo)識(shí)為實(shí)際海冰網(wǎng)格點(diǎn)的行列號(hào)，進(jìn)而使用者通過查找指定行列號(hào)的網(wǎng)格點(diǎn)是否存在于該數(shù)據(jù)表中，就可以知曉指定的網(wǎng)格點(diǎn)是否是實(shí)際海冰網(wǎng)格點(diǎn)。第二種是該表中記錄每一個(gè)網(wǎng)格點(diǎn)的行列號(hào)，并且在記錄行列號(hào)的同時(shí)，為每個(gè)網(wǎng)格點(diǎn)設(shè)置標(biāo)記，用以標(biāo)識(shí)該網(wǎng)格點(diǎn)是否是實(shí)際海冰網(wǎng)格點(diǎn)。具體如實(shí)際海冰網(wǎng)格點(diǎn)標(biāo)記為1，非實(shí)際海冰網(wǎng)格點(diǎn)則標(biāo)記為0。

在確定了數(shù)據(jù)表之后，可以按照如下方式使用該數(shù)據(jù)表，即本申請所提供的海冰標(biāo)識(shí)確定方法，還包括：

若獲取到的后向散射系數(shù)觀測結(jié)果，則提取后向散射系數(shù)觀測結(jié)果所對應(yīng)的后向散射系數(shù)測量脈沖的第一經(jīng)緯度坐標(biāo)；

計(jì)算第一經(jīng)緯度坐標(biāo)在等經(jīng)緯度網(wǎng)格中所對應(yīng)的目標(biāo)行列號(hào)；

在海冰標(biāo)識(shí)統(tǒng)計(jì)表中查找目標(biāo)行列號(hào)所對應(yīng)的網(wǎng)格是否為實(shí)際海冰網(wǎng)格點(diǎn)，若是，則將后向散射系數(shù)觀測結(jié)果設(shè)置海冰標(biāo)識(shí)。

也就是，在確定了海冰標(biāo)識(shí)統(tǒng)計(jì)表后，再獲取到后向散射系數(shù)觀測結(jié)果時(shí)，可以先將后向散射系數(shù)測量脈沖的第一經(jīng)緯度坐標(biāo)進(jìn)行坐標(biāo)轉(zhuǎn)換，以得到對應(yīng)的目標(biāo)行列號(hào)，之后再依據(jù)該目標(biāo)行列號(hào)在海冰標(biāo)識(shí)統(tǒng)計(jì)表中進(jìn)行查找，最后如果該目標(biāo)行列號(hào)對應(yīng)的網(wǎng)格點(diǎn)是實(shí)際海冰網(wǎng)格點(diǎn)，則對后向散射系數(shù)觀測結(jié)果設(shè)置海冰標(biāo)識(shí)，反之，則不設(shè)置海冰標(biāo)識(shí)，或者是設(shè)置無海冰標(biāo)識(shí)。

下面，以一個(gè)具體的實(shí)例來說明本申請所提供的海冰標(biāo)識(shí)確定方法。

步驟1、根據(jù)待海冰標(biāo)識(shí)的L1B數(shù)據(jù)文件名(該數(shù)據(jù)文件名是使用時(shí)間進(jìn)行命名的)，確定觀測時(shí)間(如年、月、日)。

步驟2、海冰標(biāo)識(shí)查找表初始化

設(shè)置大小為720×1440的二維數(shù)組作為海冰標(biāo)識(shí)查找表ice_mask(720,1440)(對應(yīng)分辨率為0.25°×0.25°等經(jīng)緯度網(wǎng)格，第一維度對應(yīng)緯度，第二維度對應(yīng)經(jīng)度)，并對每個(gè)網(wǎng)格賦初值為0。

步驟3、根據(jù)步驟1中確定的觀測時(shí)間，提取距當(dāng)前觀測時(shí)間一星期以內(nèi)的SSMI海冰數(shù)據(jù)產(chǎn)品，形成SSMI產(chǎn)品文件列表。

步驟4、根據(jù)步驟2中生成的產(chǎn)品文件列表，依次讀取每個(gè)SSMI海冰數(shù)據(jù)文件，并將其轉(zhuǎn)換到步驟2中得到的分辨率為0.25°×0.25°等經(jīng)緯度網(wǎng)格(與海冰標(biāo)識(shí)查找表對應(yīng))中。并對海冰密集度大于0的海冰標(biāo)識(shí)查找表網(wǎng)格單元賦值為1；(即某個(gè)網(wǎng)格單元的位置上只要在一個(gè)星期之內(nèi)有過海冰，就把這個(gè)網(wǎng)格節(jié)點(diǎn)標(biāo)記為海冰標(biāo)識(shí)，即賦值為1)。

步驟5、根據(jù)觀測時(shí)間，提取包含觀測時(shí)間當(dāng)天的ECMWF預(yù)報(bào)海溫?cái)?shù)據(jù)文件名。

步驟6、根據(jù)步驟5生成的ECMWF預(yù)報(bào)海溫文件名，讀取海溫預(yù)報(bào)數(shù)據(jù)，并將其轉(zhuǎn)換到步驟1中得到的分辨率為0.25°×0.25°等經(jīng)緯度網(wǎng)格上，并對海面溫度小于2℃的海冰標(biāo)識(shí)查找表網(wǎng)格單元賦值為1；

步驟7、分別以每個(gè)賦值為1的網(wǎng)格單元為中心，向外延拓50km。具體操作方法為，對海冰標(biāo)識(shí)查找表ice_mask的每個(gè)單元進(jìn)行循環(huán)，若該單元值為1，則以該單元為中心的5×5窗口大小的單元均賦值為1。進(jìn)而，建立了信息完整的海冰標(biāo)識(shí)查找表ice_mask。

步驟8、根據(jù)海冰標(biāo)識(shí)查找表ice_mask，對L1B中的每個(gè)后向散射系數(shù)觀測結(jié)果進(jìn)行海冰標(biāo)識(shí)。具體過程如下：

1)提取后向散射系數(shù)測量脈沖對應(yīng)的經(jīng)緯度，并分別標(biāo)記為lat,lon。

2)利用下式，計(jì)算后向散射系數(shù)測量脈沖所對應(yīng)的風(fēng)矢量單元在海冰標(biāo)識(shí)查找表ice_mask中對應(yīng)的行列號(hào)；

I_index＝ROUND((lat+90.)*IGRID/180-0.5)；

J_index＝ROUND(lon*YGRID/360-0.5)；

其中，I_index為海冰查找表二維數(shù)組數(shù)組元素的行號(hào)，J_index為海冰查找表二維數(shù)組數(shù)組元素的列號(hào)，ROUND為取整符號(hào),lat為后向散射測量脈沖中心對應(yīng)的緯度，lon為后向散射系數(shù)測量脈沖中心對應(yīng)的經(jīng)度。IGRID一般取4，對應(yīng)0.25°的緯向分辨率。JGRID一般取4，對應(yīng)0.25°的經(jīng)向分辨率。

3)海冰標(biāo)識(shí)賦值

根據(jù)上一個(gè)步驟中獲得的行列號(hào)，查找在海冰標(biāo)識(shí)查找表ice_mask的對應(yīng)網(wǎng)格單元，若該網(wǎng)格單元值為1，則對該后向散射系數(shù)觀測結(jié)果置海冰標(biāo)識(shí)，否則不置海冰標(biāo)識(shí)。

由于散射計(jì)在地面觀測足印為橢圓形，采用傳統(tǒng)的方法，需要判斷該橢圓是否全部或者部分落在海冰區(qū)域，這就需要首先計(jì)算觀測組足印的中心點(diǎn)(對應(yīng)天線發(fā)射的波束在地面的交點(diǎn))，橢圓的朝向(對應(yīng)天線觀測的方位角)、橢圓的大小(對應(yīng)天線的方向圖)，并在此基礎(chǔ)上，確定該橢圓是否有部分或者全部區(qū)域落在海冰上，計(jì)算量巨大。采用本申請所提供的方法，由于波束中心點(diǎn)距橢圓的邊緣不會(huì)超過50km，通過將查找表外延50km的方法，只需要確保波束中心點(diǎn)所對應(yīng)的查找表單元沒有海冰，即可保證該觀測結(jié)果沒有被海冰影響。

與前述的方法相對應(yīng)的，本申請還提供了海冰標(biāo)識(shí)確定裝置，如圖4所示，包括：

建立模塊201，用于建立地球表面的等經(jīng)緯度網(wǎng)格；

第一獲取模塊202，用于獲取預(yù)定時(shí)間內(nèi)的等經(jīng)緯度網(wǎng)格的每個(gè)網(wǎng)格點(diǎn)的海冰密集度數(shù)值；

第二獲取模塊203，用于獲取等經(jīng)緯度網(wǎng)格的各個(gè)網(wǎng)格點(diǎn)的海面溫度數(shù)值；

確定模塊204，用于根據(jù)海冰密集度數(shù)值和海面溫度數(shù)值確定等經(jīng)緯度網(wǎng)格的海冰標(biāo)識(shí)情況。

優(yōu)選的，預(yù)定時(shí)間為一周。

所屬領(lǐng)域的技術(shù)人員可以清楚地了解到，為描述的方便和簡潔，上述描述的系統(tǒng)、裝置和單元的具體工作過程，可以參考前述方法實(shí)施例中的對應(yīng)過程，在此不再贅述。

所述作為分離部件說明的單元可以是或者也可以不是物理上分開的，作為單元顯示的部件可以是或者也可以不是物理單元，即可以位于一個(gè)地方，或者也可以分布到多個(gè)網(wǎng)絡(luò)單元上?？梢愿鶕?jù)實(shí)際的需要選擇其中的部分或者全部單元來實(shí)現(xiàn)本實(shí)施例方案的目的。

所述功能如果以軟件功能單元的形式實(shí)現(xiàn)并作為獨(dú)立的產(chǎn)品銷售或使用時(shí)，可以存儲(chǔ)在一個(gè)計(jì)算機(jī)可讀取存儲(chǔ)介質(zhì)中?；谶@樣的理解，本發(fā)明的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻(xiàn)的部分或者該技術(shù)方案的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，該計(jì)算機(jī)軟件產(chǎn)品存儲(chǔ)在一個(gè)存儲(chǔ)介質(zhì)中，包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本發(fā)明各個(gè)實(shí)施例所述方法的全部或部分步驟。而前述的存儲(chǔ)介質(zhì)包括：U盤、移動(dòng)硬盤、只讀存儲(chǔ)器(ROM，Read-Only Memory)、隨機(jī)存取存儲(chǔ)器(RAM，Random Access Memory)、磁碟或者光盤等各種可以存儲(chǔ)程序代碼的介質(zhì)。

以上所述，僅為本發(fā)明的具體實(shí)施方式，但本發(fā)明的保護(hù)范圍并不局限于此，任何熟悉本技術(shù)領(lǐng)域的技術(shù)人員在本發(fā)明揭露的技術(shù)范圍內(nèi)，可輕易想到變化或替換，都應(yīng)涵蓋在本發(fā)明的保護(hù)范圍之內(nèi)。因此，本發(fā)明的保護(hù)范圍應(yīng)所述以權(quán)利要求的保護(hù)范圍為準(zhǔn)。