SeaWinds散射計風矢量反演的象元篩選、濾波方法和系統(tǒng)的制作方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及海洋遙感領(lǐng)域，特別是涉及一種SeaWinds散射計風矢量反演的象元 篩選、濾波方法和系統(tǒng)。
【背景技術(shù)】
[0002] 海風是作用在海洋表層的重要動力因子，驅(qū)動著海洋從海面波到深層洋流系統(tǒng)的 各種尺度的海洋運動，小則形成海洋波浪，大則推動洋流發(fā)展。海風通過調(diào)節(jié)海水與大氣之 間的熱量、水汽和化學物質(zhì)的交換，來調(diào)整海-氣耦合作用，維持全球和區(qū)域氣候模式，該 耦合作用對全球和區(qū)域氣候有重要影響，甚至能引起全球環(huán)境的變化，如厄爾尼諾現(xiàn)象。而 且，海風可以納入到區(qū)域或全球的數(shù)值天氣預報系統(tǒng)，提高天氣預報的能力。因此，獲得高 精度、高時空分辨率的海風數(shù)據(jù)對于海洋動力學、氣象、氣候?qū)W研究和人類合理利用風能資 源具有十分重要的科研價值和現(xiàn)實意義。
[0003] 海面風場的常規(guī)觀測系統(tǒng)主要通過船舶、海上浮標及沿岸站等。對于覆蓋全球面 積70%的海洋而言，常規(guī)手段獲得的風場資料十分有限而且花費巨大。衛(wèi)星遙感技術(shù)為海 面風場的測量提供了一個全新的手段，它具有覆蓋面廣、觀測連續(xù)、受天氣影響小等優(yōu)點， 已經(jīng)進入業(yè)務化運行。能夠提供全球海面風場信息的星載傳感器包括散射計、微波輻射計、 合成孔徑雷達（SAR)、衛(wèi)星高度計，但只有散射計能大范圍同時提供包括風速大小和風向在 內(nèi)的完整的風矢量信息。
[0004] 散射計反演風場分為三個個步驟：第一步，建立地球物理模型函數(shù)（geophysical modelfunction,GMF)，第二步，由地球物理模型函數(shù)反演得到2-4個模糊解，第三步，采用 圓中數(shù)濾波方法去除模糊解，得到唯一的真解。圓中數(shù)濾波算法有效運行的前提條件是濾 波窗口中的風矢量誤差比較大（風向反演誤差大于45度）的象元數(shù)目不能超過一半。由 于高風速時SeaWinds散射計觀測的反演誤差整體較大，因而在使用圓中數(shù)濾波算法時，濾 波窗口中反演誤差較大的象元數(shù)目往往超過一半，使得圓中數(shù)濾波效果不佳，使得最終反 演誤差偏高。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的一個目的在于提供一種SeaWinds散射計風矢量反演的象元篩選方法和 系統(tǒng)，可以有助于提高高風速情形下的風矢量反演精度。
[0006] 本發(fā)明的這一目的通過如下技術(shù)方案實現(xiàn)：
[0007] -種SeaWinds散射計風矢量反演的象元篩選方法，包括如下步驟：
[0008] 確定SeaWinds散射計風矢量反演中的象元的各模糊解以及各所述模糊解的總灰 度值；
[0009] 將各所述模糊解與數(shù)值風場數(shù)據(jù)進行對比，找出與數(shù)值風場數(shù)據(jù)最接近的目標模 糊解；
[0010] 判斷所述目標模糊解的總灰度值是否大于預設(shè)的灰度閾值；
[0011] 若是，對所述象元作一標記。
[0012] -種SeaWinds散射計風矢量反演的象元篩選系統(tǒng),包括：
[0013] 處理模塊，用于確定SeaWinds散射計風矢量反演中的象元的各模糊解以及各所 述模糊解的總灰度值；
[0014] 對比模塊，用于將各所述模糊解與數(shù)值風場數(shù)據(jù)進行對比，找出與數(shù)值風場數(shù)據(jù) 最接近的目標模糊解；
[0015] 第一判斷模塊，用于判斷所述目標模糊解的總灰度值是否大于預設(shè)的灰度閾值；
[0016] 標記模塊，用于在所述第一判斷模塊的判斷結(jié)果為是時，對所述象元作一標記。
[0017] 根據(jù)上述本發(fā)明的方案，其是首先確定象元的各模糊解及其相應的總灰度值，再 通過將各模糊解與數(shù)值風場數(shù)據(jù)進行對比獲得目標模糊解，判斷該目標模糊解的總灰度值 是否大于預設(shè)的灰度閾值，若是，對象元作一標記，SeaWinds散射計風矢量反演中的各象元 都可以按照這種方式進行篩選，最后得到的有標記的象元即為篩選出的反演誤差較小的象 元，可以基于這些象元進行圓中數(shù)濾波，由于篩選過程可有效剔除風矢量反演誤差較大的 象元，從而較大程度地提高了圓中數(shù)濾波算法的精度，即提高了最終的風矢量反演精度。
[0018] 本發(fā)明的另一目的在于提供一種SeaWinds散射計風矢量反演的濾波方法和系 統(tǒng)，可以有助于提高高風速情形下的風矢量反演精度。
[0019] 本發(fā)明的這一目的通過如下技術(shù)方案實現(xiàn)：
[0020] -種SeaWinds散射計風矢量反演的濾波方法，包括如下步驟：
[0021] 采用如上所述的SeaWinds散射計風矢量反演的象元篩選方法對SeaWinds散射計 風矢量反演中的象元進行篩選；
[0022] 選取濾波窗口，獲取該濾波窗口中的有標記的象元數(shù)目；
[0023] 判斷所述象元數(shù)目是否大于零；
[0024] 若是，則根據(jù)所述象元數(shù)目對應的預設(shè)計算方式計算所述濾波窗口的風矢量；
[0025] 若否，則根據(jù)L2B數(shù)據(jù)的獲得所述濾波窗口的風矢量。
[0026] -種SeaWinds散射計風矢量反演的濾波系統(tǒng)，包括如上所述的SeaWinds散射計 風矢量反演的象元篩選系統(tǒng)，還包括：
[0027] 選取模塊，用于選取濾波窗口，獲取該濾波窗口中的有標記的象元數(shù)目；
[0028] 第二判斷模塊，用于判斷所述象元數(shù)目是否大于零；
[0029] 濾波模塊，用于在所述第二判斷模塊的判斷結(jié)果為是時，根據(jù)所述象元數(shù)目對應 的預設(shè)計算方式計算所述濾波窗口的風矢量，在所述第二判斷模塊的判斷結(jié)果為否時，根 據(jù)L2B數(shù)據(jù)的獲得所述濾波窗口的風矢量。
[0030]根據(jù)上述本發(fā)明的方案，其是首先采用如上所述的SeaWinds散射計風矢量反演 的象元篩選方法對高風速情形下SeaWinds散射計風矢量反演中的象元進行篩選，再選取 濾波窗口，獲取該濾波窗口中的有標記的象元數(shù)目，判斷所述象元數(shù)目是否大于零，若是， 則根據(jù)所述象元數(shù)目對應的預設(shè)計算方式計算所述濾波窗口的風矢量，若否，則根據(jù)L2B 數(shù)據(jù)的獲得所述濾波窗口的風矢量，由于篩選過程可有效剔除濾波窗口中的風矢量反演誤 差較大的象元，從而較大程度地提高了圓中數(shù)濾波算法的精度；同時，由于高風速情形下反 演誤差很大的象元出現(xiàn)幾率較大，在整個濾波窗口中總灰度值大于設(shè)定的灰度閾值的象元 的數(shù)目并不穩(wěn)定，可能為〇,可能為1，也可能為2,根據(jù)不同的象元數(shù)目對應的計算方式計 算得到濾波窗口的風矢量，也提高了最終的風矢量反演精度。
【附圖說明】
[0031] 圖1為SeaWinds散射計地球物理模型函數(shù)的以雷達后向散射系數(shù)和相對方位角 為坐標的模型表；
[0032] 圖2為同一象元觀測值在風矢量解空間中的圖譜示意圖；
[0033] 圖3為本發(fā)明的SeaWinds散射計風矢量反演的象元篩選方法實施例的流程示意 圖；
[0034] 圖4為模糊解區(qū)域交點分布示意圖；
[0035] 圖5為確定模糊解的總灰度值的方式在其中一個實施例中的流程示意圖；
[0036] 圖6為圖譜算法風矢量反演誤差與總灰度關(guān)系示意圖；
[0037] 圖7為本發(fā)明的SeaWinds散射計風矢量反演的濾波方法實施例的流程示意圖；
[0038] 圖8為本發(fā)明的SeaWinds散射計風矢量反演的象元篩選系統(tǒng)實施例的結(jié)構(gòu)示意 圖；
[0039] 圖9為圖8中的處理模塊在其中一個實施例中的結(jié)構(gòu)示意圖；
[0040] 圖10為本發(fā)明的SeaWinds散射計風矢量反演的濾波系統(tǒng)實施例的流程示意圖。
【具體實施方式】
[0041] 為使本發(fā)明的目的、技術(shù)方案及優(yōu)點更加清楚明白，以下結(jié)合附圖及實施例，對本 發(fā)明進行進一步的詳細說明。應當理解，此處所描述的【具體實施方式】僅僅用以解釋本發(fā)明， 并不限定本發(fā)明的保護范圍。
[0042] 為了便于理解本發(fā)明的方案，以下首先對傳統(tǒng)的散射計風場反演方式進行說明。
[0043] 1)傳統(tǒng)風場反散射計風場反演方式的一般步驟
[0044] 傳統(tǒng)風場反散射計風場反演方式可分為三個步驟：建立地球物理模型函數(shù)，風矢 量解搜索（反演算法）以及風矢量模糊解去除。從1978年第一顆專門用于海洋遙感的衛(wèi)星 SEASAT-A發(fā)射到現(xiàn)在，國外學者提出了多種用于散射計海面風場遙感的模型函數(shù)及反演算 法。由于地球物理模型函數(shù)的雙調(diào)和特性以及觀測中存在的誤差，使得風矢量搜索的結(jié)果 不唯一，而是存在1-4個極小值，稱為模糊解。由于從單一象元很難獲得足夠的信息將其排 除，因而需要一個獨立的模糊解去除程序，將其排除以便獲得唯一的真解。
[0045] 2)散射計工作的基本原理
[0046] 散射計是一種特殊的主動微波傳感器，主要是用來獲取海面雷達截面積等多個測 量數(shù)據(jù)。通過向地表發(fā)射一定波長的微波，將那些在同一地面單元具有不同幾何觀測參數(shù) 的測量數(shù)據(jù)結(jié)合起來，經(jīng)過一定處理而獲得該單元海面風速和風向。散射計對海面風矢量 的測量是一種間接關(guān)系。
[0047] 海面微波散射主要有兩種機制：當入射角小于20°時，微波回向散射主要是鏡面 反射；當入射角大于20°時，微波回向散射主要是布拉格散射。而散射計的入射角一般大 于20°，如SeaWinds散射計內(nèi)外兩波束的入射角分別為46°和54°，所以散射計工作時， 海面微波散射以布拉格散射為主。在入射波長為微波的條件下，能引起布拉格散射的只能 是疊加于海浪上的厘米波或毛細波，而且布拉格散射與厘米波波峰與雷達視線之間的夾角 有很大的相關(guān)性。
[0048] 風從大氣向海面?zhèn)鬟f動量時，會使海面變得粗糙，形成小尺度的厘米波，厘米波數(shù) 量與海面上的磨擦風速成正比相關(guān)，而厘米波波峰方向?qū)︼L向敏感，所以能通過對微米波 的觀測來反推出海面的風矢量與風向。散射計就是測量海面的后向散射能量大小，并估算 標準化的海面雷達截面積NRCS(NormalizedRadarCrossSection)。
[0049] 3)海面風場反演的地球物理模型函數(shù)
[0050] 利用散射計在不同的參數(shù)條件下（入射角和極化方式）測得的海面后徑雷達散射 截面積進行風場反演，需要知道海面雷達后向散射截面積（NRCS)和風速、相對方位角（風 向與雷達觀測方向之間的夾角）、入射角以及極化方式等幾個變量之間的關(guān)系。模型函數(shù)的 目的就是準確描述雷達后向散射截面積（NRCS)與海面環(huán)境參數(shù)和儀器參數(shù)的關(guān)系。模型 函數(shù)對于散射計的設(shè)計及反演算法具有關(guān)鍵的作用，它的精度可大大地影響到風場反演的 精度以及散射計設(shè)計的可行性等重要問題。
[0051] 從第一顆專門用于海洋遙感的衛(wèi)星SEASAT發(fā)射到現(xiàn)在，人們不斷提出用于散射 計海面風場遙感的模型函數(shù)和反演方法，模型函數(shù)的非線性決定了風場反演方法