王海霞，徐 進(jìn)，王慶名，陳 麗，龐璽斌

(1.大連海事大學(xué) 航海學(xué)院，遼寧 大連 116026；2.大連海事大學(xué) 校友事務(wù)與合作處，遼寧 大連 116026)

0 引言

海洋風(fēng)場是海洋表層運(yùn)動(dòng)的主要推力，它管控能源和動(dòng)力的轉(zhuǎn)移，能夠?qū)⒖諝鈳У胶Ｑ?，在氣候變化中占有非常重要的作用。海面風(fēng)場常被應(yīng)用于各類海洋科學(xué)研究，如海洋環(huán)流、熱帶風(fēng)暴、大氣與海洋的耦合效應(yīng)等[1]。因此，監(jiān)測海面風(fēng)場變化有著極其重要的意義。

目前，國外在航海雷達(dá)海上風(fēng)場的反演方面已經(jīng)存在一些研究成果。Ziemer等人[2]首次應(yīng)用二維傅里葉變換，分析2張帶有海雜波的雷達(dá)圖像，計(jì)算的海浪譜與傳統(tǒng)浮標(biāo)獲取的海浪譜值較為接近。Alpers等[3]和Vesecky等[4]對(duì)雷達(dá)海浪成像作了詳細(xì)描述，但他們均未考慮流的影響。Lee等人[5]基于電磁波海面散射機(jī)理，指出海面的粗糙度與風(fēng)速之間，呈現(xiàn)正相關(guān)特性。Hatten等人[6]發(fā)現(xiàn)航海雷達(dá)散射截面風(fēng)向及電磁波照射方向存在一定關(guān)系，運(yùn)用這種關(guān)系可以判斷出風(fēng)向，但此方法存在一定的局限性，往往不能滿足所需的全方位航海雷達(dá)散射截面。Dankert等人[7]提出了一種基于光流運(yùn)動(dòng)估計(jì)的技術(shù)，從X波段航海雷達(dá)圖像中獲得風(fēng)速和風(fēng)向，該方法相對(duì)簡單，不需對(duì)雷達(dá)系統(tǒng)的相位進(jìn)行校正，也不存在角度模糊問題，所獲得的海面風(fēng)場具有較高的時(shí)空分辨率。Liu等人[8]針對(duì)未被遮擋的低散射強(qiáng)度雷達(dá)圖像，提出了一種散射強(qiáng)度值選取法，并用對(duì)偶曲線擬合法減小了風(fēng)向反演誤差。

國內(nèi)應(yīng)用雷達(dá)監(jiān)測海面風(fēng)場研究的起步較晚，楊勁松[9]用合成孔徑雷達(dá)(SAR)，采用二維傅里葉變換波數(shù)譜的方法，反演出了風(fēng)向，并用COMD4模式函數(shù)算出了風(fēng)速。金麗潔[10]利用浪高反演算法[11]，采用多重信號(hào)分類算法完成定向，實(shí)現(xiàn)了浪高場的反演。張璇等人[12]分析目前交叉極化對(duì)風(fēng)場反演的優(yōu)勢，總結(jié)了多極化SAR圖像反演海面風(fēng)場的最新研究成果，并對(duì)多極化SAR圖像反演海面風(fēng)場的業(yè)務(wù)應(yīng)用發(fā)展進(jìn)行了初步展望。

本文基于X波段航海雷達(dá)反演海面風(fēng)場的基本原理，以采集到的雷達(dá)原始圖像為數(shù)據(jù)基礎(chǔ)，基于梯度算法求得主風(fēng)向，并與大連氣象臺(tái)的實(shí)測值進(jìn)行對(duì)比，以期為海面風(fēng)場的實(shí)時(shí)監(jiān)測提供有效的技術(shù)手段。

1 海面風(fēng)場反演方法

1.1 數(shù)據(jù)來源及設(shè)備參數(shù)

X波段船載雷達(dá)信號(hào)采集到的原始灰度圖像(設(shè)備參數(shù)詳見表1)監(jiān)測范圍為0.75海里，圖像大小為1 024×1 024(相同地點(diǎn)連續(xù)獲取的3幅圖像如圖1所示)。監(jiān)測數(shù)據(jù)采集位置：大連港附近海域(121.830°E，38.913°N)。船載平臺(tái)：大連海事大學(xué)“育鯤”輪號(hào)，采集時(shí)間為2015年8月14日20時(shí)4分。數(shù)據(jù)中含有大量的含雜波信息，用于反演海面海浪與風(fēng)向等信息的研究。

表1 育鯤輪航海雷達(dá)設(shè)備參數(shù)

參數(shù)項(xiàng)參數(shù)值工作極化方式水平工作頻段X波段雷達(dá)天線長度/ft8天線旋轉(zhuǎn)速度/RPM28脈沖工作頻率/Hz3 000/1 800 /785雷達(dá)控制臺(tái)23Inch液晶工控機(jī)工作量程/海里1～12

圖1 雷達(dá)圖像原圖

1.2 實(shí)驗(yàn)流程

海面風(fēng)場反演過程包括雷達(dá)信號(hào)的實(shí)時(shí)采集和海面風(fēng)場信息提取2部分。航海雷達(dá)風(fēng)場反演流程如圖2所示。

圖2 航海雷達(dá)風(fēng)場反演流程

1.3 圖像預(yù)處理

1.3.1 圖像求和

圖像求和即將圖像對(duì)應(yīng)像素點(diǎn)的值進(jìn)行求和，從而增強(qiáng)海浪信息所在位置的灰度強(qiáng)度。n幅圖像中，每個(gè)像素點(diǎn)f(x，y)求和公式為

(1)

1.3.2 中值濾波法

中值濾波基于排序統(tǒng)計(jì)理論，將每一像素的灰度值設(shè)置為該點(diǎn)某鄰域窗口內(nèi)的所有像素點(diǎn)灰度值的中值，公式為

g(x，y)=median({f(x-i，y-j)，(i，j)∈W)}，

(2)

式中，g(x，y)為輸出像素值；f(x-i，y-i)表示輸入像素值；W是含有基數(shù)點(diǎn)的模板窗口，一般為3×3窗口或5×5窗口。

1.3.3 梯度法反演風(fēng)向

風(fēng)向反演主要是基于海面風(fēng)引起的紋理可在雷達(dá)回波圖像中成像，Dankertd等人[13]研究發(fā)現(xiàn)雷達(dá)圖像中存在與風(fēng)向平行的紋理信息，分析雷達(dá)圖像中的海浪紋理可以獲取風(fēng)向[14]。梯度是用于描述圖像灰度變化情況的二維列向量，能反映變化的劇烈程度[15]。通過將圖像進(jìn)行分割并統(tǒng)計(jì)梯度方向與大小得到該區(qū)域的梯度特征分布的統(tǒng)計(jì)結(jié)果[16]，梯度的角度為圖像像素變化的方向，梯度的模為梯度方向的變化率。若求出海浪紋理變化最大的梯度方向，即可得到與之垂直的主風(fēng)向。算法實(shí)施步驟如下。

(1)圖像重采樣

對(duì)雷達(dá)原始圖像采用含有一定濾波作用的算子進(jìn)行計(jì)算，平滑圖像的同時(shí)，壓縮圖像像素?cái)?shù)量?？捎盟阕庸饺缦拢?/p>

R=B2S2B4，

(3)

(4)

(5)

式中，S2表示2×2線性插值平均，先用四階算子平滑，然后縮減圖像，最后用二階算子再次平滑圖像。

(2)局部梯度的計(jì)算

用最優(yōu)化的Sobel算子進(jìn)行局部梯度計(jì)算：

(6)

Dy=DxT，

(7)

式中，T表示轉(zhuǎn)置矩陣。此時(shí)圖像的梯度G為

G=(Dx+iDy)*(A)，

(8)

式中，A表示圖像；*表示卷積，G′與G″的計(jì)算方式如下所示：

G′=R*(G2)，

(9)

G″=R*(|G2|)。

(10)

由文獻(xiàn)[17]可知，一致性參數(shù)C為：

(11)

(3)確定風(fēng)向

應(yīng)用一致性參數(shù)和梯度的模，做梯度方向角度的加權(quán)直方圖。因?yàn)轱L(fēng)向與梯度變化最大的方向垂直，在直方圖中最大尖峰值加上90°，即為主風(fēng)向。

2 實(shí)驗(yàn)結(jié)果與分析

2.1 雷達(dá)原始圖像求和處理結(jié)果

在獲取的雷達(dá)圖像中包含許多信息，但在獲取過程中由于各種原因會(huì)存在不同形式的干擾。因此在反演前應(yīng)當(dāng)先進(jìn)行圖像處理。首先將連續(xù)獲取的雷達(dá)原始圖像進(jìn)行圖像求和，增強(qiáng)海浪信息所在位置的灰度強(qiáng)度，如圖3所示。

圖3 雷達(dá)原始圖像求和

2.2 中值濾波處理結(jié)果

采用5×5窗口，對(duì)雷達(dá)原始圖像求和處理后的圖像進(jìn)行中值濾波處理，降低同頻干擾，如圖4所示。在進(jìn)行了預(yù)處理后，為消除圖像存在的尾跡效應(yīng)，截取了沒有尾跡效應(yīng)的部分截取了沒有尾跡效應(yīng)的部分，即(275：402，445：572)中的圖像窗口進(jìn)行風(fēng)向反演。

圖4 中值濾波處理

2.3 局部梯度法反演風(fēng)向

采用Matlab編寫程序?qū)D像進(jìn)行梯度計(jì)算，并做出相應(yīng)的梯度加權(quán)直方圖，如圖5所示。由圖5可以看出，在梯度角度為88°時(shí)出現(xiàn)的頻率最大。因?yàn)橹黠L(fēng)向與梯度角度直方圖的峰值方向垂直，應(yīng)當(dāng)在獲得的梯度角度加上90°，即為所求的主風(fēng)向178°。

圖5 梯度加權(quán)直方圖

3 驗(yàn)證

根據(jù)大連氣象臺(tái)提供的2015年8月14日平均風(fēng)向?yàn)?83°，與實(shí)測結(jié)果相比風(fēng)向偏差5°，反演的風(fēng)向與預(yù)報(bào)結(jié)果是比較吻合的。因本文所用的雷達(dá)數(shù)據(jù)極化方式為水平極化，可以很好地消除180°模糊問題，因此，反演結(jié)果不存在180°模糊。

4 結(jié)束語

以大連海事大學(xué)教學(xué)實(shí)習(xí)船“育鯤”輪采集的航海雷達(dá)原始數(shù)據(jù)為基礎(chǔ)，基于梯度算法進(jìn)行風(fēng)場反演。在實(shí)驗(yàn)過程中，存在多種因素可能會(huì)影響計(jì)算的結(jié)果。首先，在雷達(dá)成像過程中，獲取到的風(fēng)場信息包含了雷達(dá)重力波和毛細(xì)波的影響。此外，在處理同頻干擾和尾跡效應(yīng)問題時(shí)，雖然采用了中值濾波，但干擾并未完全消除。本文提出的方法還未考慮其他因素，如沒有考慮到船舶運(yùn)動(dòng)對(duì)于雷達(dá)回波的影響；遇到遮擋物時(shí)，在遮擋物后方向形成盲區(qū)，不能準(zhǔn)確顯示出該區(qū)域的真實(shí)海面狀況。這些對(duì)于風(fēng)場反演的結(jié)果都將存在影響，使結(jié)果產(chǎn)生誤差。而且，計(jì)算結(jié)果為風(fēng)場的主風(fēng)向，并不是該區(qū)域內(nèi)的所有風(fēng)向，各區(qū)域點(diǎn)之間也會(huì)存在差異。上述難點(diǎn)都需要在未來的工作中，繼續(xù)研究。