朱曉雯，侯 宇，劉玉海，吳克儉**

(1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東 青島 266100；2.中國(guó)海洋大學(xué)數(shù)學(xué)科學(xué)學(xué)院，山東 青島 266100；3.中國(guó)科學(xué)院海洋研究所，山東 青島 266100；4.中科曙光國(guó)際信息產(chǎn)業(yè)有限公司,山東 青島 266100)

有效波高(Significant Wave Height, SWH)是海浪研究中一項(xiàng)重要的海洋要素。在物理海洋學(xué)中，有效波高在傳統(tǒng)上被定義為最高三分之一波(H1/3)的平均波高[1]，即前三分之一大波波高。有效波高是重要的海浪要素，也是描述海況主要參量之一，被廣泛應(yīng)用于海洋數(shù)值模擬、海岸工程和航海運(yùn)輸?shù)阮I(lǐng)域。高分辨率的SWH數(shù)據(jù)能為多尺度物理海洋學(xué)研究以及海浪精細(xì)化預(yù)報(bào)、防災(zāi)減災(zāi)、海氣耦合模式的發(fā)展等工作提供更加有效的數(shù)據(jù)支撐。在推動(dòng)海-氣耦合模式的過程中，需要對(duì)每個(gè)海洋網(wǎng)格增強(qiáng)像素化，這樣才能更好地解析中小尺度的氣候變化過程，然而隨著空間分辨率的提高，所需要的計(jì)算量級(jí)也急速增長(zhǎng)。

隨著人工智能技術(shù)的迅速發(fā)展和顯卡(GPU)在加速運(yùn)算中的廣泛使用，深度學(xué)習(xí)方法被廣泛運(yùn)用于海洋和氣象等自然科學(xué)研究中[2-8]。Wang等[9]歸納總結(jié)了適用于不同類型時(shí)空數(shù)據(jù)的深度學(xué)習(xí)方法，并分類統(tǒng)計(jì)了近年來深度學(xué)習(xí)方法在氣候和天氣學(xué)中的研究結(jié)果。除了傳統(tǒng)的插值方法之外，近年來有很多深度學(xué)習(xí)的方法也被用來解決海洋數(shù)據(jù)的降尺度問題。降尺度是從低分辨率(Low-Resolution, LR)重建高分辨率(High-Resolution, HR)數(shù)據(jù)的過程。二維數(shù)據(jù)的降尺度處理可以對(duì)應(yīng)計(jì)算機(jī)視覺和模式識(shí)別領(lǐng)域中圖像的超分辨率問題[10]。基于神經(jīng)網(wǎng)絡(luò)的超分辨率算法的基本思想是建立從低分辨率圖像到高分辨率圖像的一個(gè)映射[11]，Aurélien等[12]使用卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)對(duì)海表面溫度(Sea Surface Temperature, SST)圖像進(jìn)行高分辨率處理，利用深度神經(jīng)網(wǎng)絡(luò)解決海洋遙感數(shù)據(jù)的降尺度問題，取得了優(yōu)于傳統(tǒng)方法的效果。在進(jìn)行大批量海洋數(shù)據(jù)的高分辨率處理時(shí)，傳統(tǒng)的簡(jiǎn)單插值方法速度快，但精確度稍低(如雙三次插值)，復(fù)雜一點(diǎn)的插值方法雖能取得較好的差之效果，但是計(jì)算復(fù)雜，處理速度較慢(如克里金插值)。基于以上事實(shí)，作者思考能否改進(jìn)深度學(xué)習(xí)的方法，找到一種既能快速計(jì)算高分辨率處理效果，又能媲美甚至超越復(fù)雜插值方法的方案。

本文運(yùn)用了多種方法對(duì)SWH數(shù)據(jù)進(jìn)行處理，包括傳統(tǒng)的雙三次插值(Bicubic interpolation)、克里金插值(Kriging interpolation)以及改進(jìn)的深度學(xué)習(xí)方法：超分辨率卷積神經(jīng)網(wǎng)絡(luò)(Super-Resolution CNN，簡(jiǎn)稱SRCNN)，用峰值信噪比(Peak Signal-to-Noise Ratio，PSNR)等多種評(píng)價(jià)指標(biāo)比較了各種方法的高分辨率處理效果。

1 數(shù)據(jù)和方法介紹

1.1 數(shù)據(jù)來源及預(yù)處理

本文的SWH數(shù)據(jù)來源于歐洲中期天氣預(yù)報(bào)中心(European Center for Medium Range Weather Forecasting, ECMWF)發(fā)布的ERA Interim, Daily數(shù)據(jù)集，與之前發(fā)布的ERA-15和ERA-40產(chǎn)品相比，該數(shù)據(jù)集使用了最新的四維變分同化技術(shù)，其中海浪數(shù)據(jù)運(yùn)用耦合海浪模式數(shù)據(jù)和大量的衛(wèi)星數(shù)據(jù)、浮標(biāo)觀測(cè)數(shù)據(jù)進(jìn)行同化[13]。選取東太平洋海區(qū)的有效波高數(shù)據(jù)，空間范圍為:109.625°W～132°W, 0°～22.375°S，空間分辨率為0.125°×0.125°，時(shí)間范圍從2011年1月1—31日，時(shí)間分辨率6 h，即每天四個(gè)時(shí)刻數(shù)據(jù)(0時(shí)、6時(shí)、12時(shí)、18時(shí))，共124個(gè)時(shí)刻的SWH場(chǎng)，數(shù)據(jù)大小為180×180。

Dee等[14]指出，在模式開發(fā)和數(shù)據(jù)同化方面，ECMWF數(shù)據(jù)具有良好的效果，而且從多個(gè)角度驗(yàn)證了ERA-Interim數(shù)據(jù)集的優(yōu)越性。Mooney等[15]在愛爾蘭冬季氣溫的模擬結(jié)果中發(fā)現(xiàn)，運(yùn)用ERA-Interim數(shù)據(jù)模擬的結(jié)果更符合實(shí)際的觀測(cè)氣溫，兩者相關(guān)性系數(shù)更高，而且比使用ERA-40和NNRP-1數(shù)據(jù)模擬的結(jié)果稍好。以上事實(shí)有力證明了本文所選數(shù)據(jù)的可靠性，基于這些數(shù)據(jù)的結(jié)果也是科學(xué)可信的。

將高分辨率數(shù)據(jù)處理成縮小六倍的低分辨率數(shù)據(jù)用于測(cè)試，數(shù)據(jù)大小變?yōu)?0×30，空間分辨率從0.125°×0.125°降為0.75°×0.75°。每種方法展示了兩次實(shí)驗(yàn)測(cè)試，測(cè)試數(shù)據(jù)選擇標(biāo)準(zhǔn)是：在測(cè)試集中選擇形態(tài)分布差異較大的兩個(gè)時(shí)刻，以保證模型對(duì)于處理不同分布狀態(tài)的低分辨率海浪數(shù)據(jù)具有普適性效果。

1.2 評(píng)價(jià)指標(biāo)介紹

1.2.1 RMSE和PSNR 均方根誤差(Root-Mean-Square Error，RMSE)是一種常用的測(cè)量數(shù)值之間差異的量度，是估計(jì)量與真實(shí)量之差期望值的平方根，對(duì)于兩個(gè)m×n的矩陣I1和I2，它們的RMSE定義如式(1)，本文中RMSE單位為米(m)。

(1)

峰值信噪比(Peak Signal-to-Noise Ratio，PSNR)經(jīng)常用作圖像處理等領(lǐng)域中信號(hào)重建質(zhì)量的測(cè)量方法，是一種評(píng)價(jià)圖像質(zhì)量的客觀標(biāo)準(zhǔn)。通常圖像經(jīng)過處理后，在某種程度上輸出的圖像會(huì)與原始圖像有所差別。為了衡量經(jīng)過處理后的圖像質(zhì)量，使用PSNR作為評(píng)價(jià)圖像質(zhì)量的指標(biāo)，判斷圖像處理程序結(jié)果是否令人滿意。PSNR通過均方誤差(Mean Square Error，MSE)進(jìn)行定義，對(duì)于兩張m×n的灰度圖像I1和I2，它們的均方誤差(MSE)定義如式(2)：

(2)

峰值信噪比PSNR定義如式(3)：

(3)

其中：MAXI表示圖像像素點(diǎn)顏色的最大數(shù)值，對(duì)于灰度圖像，圖像每個(gè)采樣點(diǎn)用8位表示，則MAXI=255。為了統(tǒng)一比較標(biāo)準(zhǔn)，在計(jì)算PSNR時(shí)使用標(biāo)準(zhǔn)化到0～255范圍的SWH數(shù)據(jù)。

1.2.2 SSIM 結(jié)構(gòu)相似度(Structural Similarity,簡(jiǎn)稱SSIM)[16]被用來測(cè)量?jī)蓚€(gè)圖像之間的相似性。SSIM的設(shè)計(jì)是為了改進(jìn)傳統(tǒng)的圖像評(píng)價(jià)方法，如峰值信噪比(PSNR)和均方誤差(MSE)。結(jié)構(gòu)相似性的基本思想假定自然影像是高度結(jié)構(gòu)化的[17]，即在自然影像中相鄰像素之間有很強(qiáng)的關(guān)聯(lián)性，而海浪有效波高數(shù)據(jù)也具有這種關(guān)聯(lián)性，而這樣的關(guān)聯(lián)性承載了海洋要素的結(jié)構(gòu)信息。相較于傳統(tǒng)圖像衡量指標(biāo)(如PSNR)，結(jié)構(gòu)相似性(SSIM)在圖像品質(zhì)的衡量上更能符合人眼對(duì)圖像品質(zhì)的判斷。

給定兩張圖像x和y，兩者的結(jié)構(gòu)相似性定義如公式(4)所示：

(4)

式中：l(x,y)是亮度項(xiàng)；c(x,y)是對(duì)比度項(xiàng)，s(x,y)是結(jié)構(gòu)項(xiàng)；α>0,β>0,γ>0，這三項(xiàng)是調(diào)整l(x,y)、c(x,y)、s(x,y)相對(duì)重要性的參數(shù)；μx及μy、σx及σy分別為x和y中像素點(diǎn)的平均值和標(biāo)準(zhǔn)差，σxy為x和y的協(xié)方差；C1,C2,C3都為常數(shù)，用來保持l(x,y)、c(x,y)、s(x,y)的穩(wěn)定。結(jié)構(gòu)相似性指標(biāo)的值越大，代表兩張圖像相似度越高。本文在具體計(jì)算時(shí)，令α=β=λ=1及C3=C2/2，得到式(5)：

(5)

計(jì)算SSIM時(shí)同樣使用標(biāo)準(zhǔn)化到0～255范圍的SWH數(shù)據(jù)。

1.2.3 余弦相似性(Cosine similarity) 對(duì)于兩個(gè)圖像(或矩陣)A和B，其余弦相似性(Cosine similarity)由點(diǎn)積和向量長(zhǎng)度給出，如式(6)所示，式中Aij和Bij分別代表A和B對(duì)應(yīng)點(diǎn)的值。

對(duì)于圖像來說，像素值都為非負(fù)值，余弦相似性的范圍從0到1，1表示兩張圖像是完全相同的，0通常表示它們之間是獨(dú)立的，而在這之間的值則表示中間的相似性或相異性。為了統(tǒng)一比較標(biāo)準(zhǔn)，在計(jì)算余弦相似度時(shí)用標(biāo)準(zhǔn)化到0～255范圍的SWH數(shù)據(jù)。

(6)

1.3 方法介紹

1.3.1 雙三次插值 在數(shù)學(xué)上，雙三次插值(Bicubic Interpolation)是常用的二維數(shù)據(jù)的插值方法。插值后得到的曲面比雙線性插值或最近鄰插值方法得到的相應(yīng)曲面更加光滑。雙三次插值可以使用拉格朗日多項(xiàng)式、三次樣條或三次卷積算法。在本文中使用由Keys[18]提出的三次卷積算法。

本文中將30×30的低分辨率數(shù)據(jù)用雙三次插值的方法高分辨率處理為180×180的規(guī)格，即分辨率提高了六倍，分辨率從0.75°×0.75°變?yōu)?.125°×0.125°。

1.3.2 Kriging 在最初在地質(zhì)統(tǒng)計(jì)學(xué)中，克里金法(Kriging)或稱高斯過程回歸法(Gaussian process regression)是一種插值方法，插值后的值由先驗(yàn)協(xié)方差控制的高斯過程建模得出。在適當(dāng)?shù)南闰?yàn)假設(shè)下，克里金法給出了對(duì)中間值的最佳線性無偏估計(jì)[19]。該方法廣泛應(yīng)用于空間分析和計(jì)算機(jī)實(shí)驗(yàn)領(lǐng)域。Jessica等[20]運(yùn)用Kriging方法確定鯨類動(dòng)物保護(hù)的潛在區(qū)域；Lin等[21]將普通克里金法應(yīng)用于氣溫和光合有效輻射量等氣候要素的插值。本文運(yùn)用了普通的克里金方法(Ordinary kriging)進(jìn)行計(jì)算。

同雙三次插值方法(Bicubic interpolation)一樣，將30×30的低分辨率數(shù)據(jù)用Kriging方法高分辨率處理為180×180的規(guī)格，分辨率從0.75°×0.75°變?yōu)?.125°×0.125°。

1.3.3 SRCNN 超分辨率卷積神經(jīng)網(wǎng)絡(luò)(Super-Resolution CNN，簡(jiǎn)稱SRCNN)是由Dong等[22]提出的解決超分辨問題的深度學(xué)習(xí)方法，計(jì)算過程如圖1所示。對(duì)于給定的低分辨率圖像Y, SRCNN的第一卷積層提取一組特征圖，第二層將這些特征映射到高分辨率的特征圖上，最后一層結(jié)合空間鄰域內(nèi)的預(yù)測(cè)，生成最終的高分辨率圖像F(Y)[22]。該方法的實(shí)質(zhì)是運(yùn)用卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Network，CNN)直接學(xué)習(xí)低、高分辨率圖像之間的映射。本次研究中，低分辨率的數(shù)據(jù)是先由高分辨率數(shù)據(jù)進(jìn)行切割，得到33×33的高分辨率切片，再將切片進(jìn)行六倍下采樣，然后對(duì)縮小的切片用臨近插值法處理得到與高分辨率數(shù)據(jù)尺寸相同的，但“模糊”的低分辨率數(shù)據(jù)。這樣就得到了許多一一對(duì)應(yīng)的低、高分辨率數(shù)據(jù)，建模學(xué)習(xí)這些低、高分辨率數(shù)據(jù)的映射關(guān)系。

圖1 SRCNN過程展示

運(yùn)用SRCNN方法時(shí)，我們需要?jiǎng)澐钟?xùn)練集和測(cè)試集，為了避免臨近時(shí)間造成訓(xùn)練集和測(cè)試集的高度相似，本文選取前100個(gè)時(shí)刻的高分辨率數(shù)據(jù)作為訓(xùn)練集，剩余的數(shù)據(jù)用作測(cè)試。

傳統(tǒng)的SRCNN超分辨率方法用于處理圖像，即數(shù)據(jù)類型為整型，范圍在0～255之間，而這并不能完全滿足實(shí)際海洋問題的應(yīng)用。Aurélien等[12]解決了SST圖像數(shù)據(jù)的高分辨率問題，但并沒有將其結(jié)果轉(zhuǎn)化為真實(shí)的SST數(shù)據(jù)，因此也無法為進(jìn)一步的海洋要素分析提供真實(shí)數(shù)據(jù)，且由于圖像數(shù)據(jù)均為整數(shù)，即使轉(zhuǎn)化為真實(shí)海洋數(shù)據(jù)，數(shù)據(jù)取值范圍也僅有256個(gè)值，數(shù)據(jù)精度有限，這就大大降低了深度學(xué)習(xí)方法的實(shí)用性。為了解決該問題，并且與傳統(tǒng)的SRCNN方法作比較，本文設(shè)計(jì)了兩種實(shí)驗(yàn)方案(SRCNN的方案一和方案二，分別用SRCNN_1、SRCNN_2表示，下同)：

SRCNN_1：將SWH數(shù)據(jù)標(biāo)準(zhǔn)化到0～255范圍，存為三通道灰度圖，此時(shí)數(shù)據(jù)均為圖片格式，數(shù)據(jù)類型為整數(shù)型，用搭建好的SRCNN架構(gòu)訓(xùn)練一個(gè)六倍超分辨率重建的模型，此模型的輸入輸出格式均為圖片格式，想要得到真實(shí)的SWH數(shù)據(jù)，需要逆標(biāo)準(zhǔn)化的過程，將0～255范圍的圖片數(shù)據(jù)轉(zhuǎn)化為真實(shí)SWH數(shù)據(jù)。

SRCNN_2：對(duì)傳統(tǒng)的SRCNN框架進(jìn)行改進(jìn)，在模型的數(shù)據(jù)讀取部分和數(shù)據(jù)預(yù)處理部分進(jìn)行改進(jìn)，使模型能做到大批量讀取真實(shí)SWH數(shù)據(jù)，并在模型內(nèi)部進(jìn)行標(biāo)準(zhǔn)化過程。重新訓(xùn)練浮點(diǎn)數(shù)矩陣數(shù)據(jù)的超分辨率模型，訓(xùn)練好的模型可以實(shí)現(xiàn)輸入低分辨率SWH數(shù)據(jù)，輸出高分辨率SWH數(shù)據(jù)，無需額外進(jìn)行數(shù)據(jù)類型的轉(zhuǎn)化，使用效果更加便捷，處理數(shù)據(jù)效率更快。

將兩種方案分別在服務(wù)器上進(jìn)行訓(xùn)練，除了訓(xùn)練集的數(shù)據(jù)類型不同外，其他參數(shù)均相同，網(wǎng)絡(luò)結(jié)構(gòu)的參數(shù)設(shè)置如表1所示。

表1 SRCNN網(wǎng)絡(luò)結(jié)構(gòu)

SRCNN方法使用搭建在LINUX服務(wù)器上的TensorFlow框架進(jìn)行計(jì)算。TensorFlow是谷歌大腦(Google Brain)團(tuán)隊(duì)開發(fā)的第二代機(jī)器學(xué)習(xí)系統(tǒng)[23]。該方法用兩塊型號(hào)為NVIDIA GeForce GTX 1080Ti的顯卡進(jìn)行計(jì)算，極大提高了計(jì)算速度。具體參數(shù)的設(shè)置見表2。

表2 SRCNN方法參數(shù)設(shè)置

3 高分辨率處理效果

對(duì)于每種方法，本文選擇了兩個(gè)SWH分布差異較大的時(shí)刻進(jìn)行實(shí)驗(yàn)，分別為2011年1月29日18時(shí)和31日6時(shí)的SWH場(chǎng)。圖2和3分別展示了兩次實(shí)驗(yàn)的低分辨率、高分辨率SWH場(chǎng)和各種方法的處理效果。

與原始高分辨率SWH相比，四種方法的高分辨率處理效果在空間分布趨勢(shì)方面均沒有很大的出入。因此在分析SWH分布這種定性問題時(shí)，雙三次插值、Kriging插值和SRCNN兩個(gè)方案的結(jié)果均能滿足要求。

從細(xì)節(jié)分析，SRCNN_1的結(jié)果如圖2(e)、3(e)所示，具有明顯的片狀結(jié)構(gòu)，SWH場(chǎng)分布連續(xù)性不如其他的方法。這是因?yàn)樵摲椒ㄊ褂玫挠?xùn)練集是圖像格式，轉(zhuǎn)化真實(shí)的SWH數(shù)據(jù)時(shí)，數(shù)值范圍為0～255的整數(shù)，而SRCNN_2的設(shè)計(jì)很好地解決了這一問題，數(shù)據(jù)連續(xù)性有了明顯提高，更加適合處理SWH的高分辨率問題。

((a)低分辨率SWH；(b)原始高分辨率SWH；(c)雙三次插值；(d)kriging插值；(e)SRCNN_1；(f)SRCNN_2。(a)Low resolution SWH；(b)Original high resolution SWH；(c)Bicubic interpolation；(d)Kriging interpolation；(e)SRCNN_1；(f)SRCNN_2.)

((a)低分辨率SWH；(b)原始高分辨率SWH；(c)雙三次插值；(d)kriging插值；(e)SRCNN_1;(f)SRCNN_2。(a)Low resolution SWH；(b)Original high resolution SWH；(c)Bicubic interpolation(d)Kriging interpolation；(e)SRCNN_1;(f)SRCNN_2.)

4 分析與討論

4.1 指標(biāo)分析

分別計(jì)算四種方法所得結(jié)果的評(píng)價(jià)指標(biāo)(SSIM，PSNR，RMSE，COSINE)(見表3)。四種方法的SSIM 和COSINE 兩項(xiàng)指標(biāo)相差不大，結(jié)果都很接近1，說明從視覺上來看，四種方法的插值效果與原始高分辨率的數(shù)據(jù)非常相似，在定性分析SWH分布特征時(shí)均能滿足要求。從PSNR指標(biāo)來看，Kriging方法最優(yōu)，其次是SRCNN_2和SRCNN_1，而這三種方法與Bicubic方法相比，在PSNR上有了較大幅度的提高。RMSE指標(biāo)方面，Kriging和SRCNN_2最優(yōu)，明顯優(yōu)于其他兩種算法，說明Kriging和SRCNN_2方法的整體誤差比較小，誤差范圍在±0.01 m左右。但是由于Kriging方法對(duì)于不同時(shí)刻的數(shù)據(jù)都要重新進(jìn)行大量計(jì)算，因此處理大批量數(shù)據(jù)時(shí)面臨計(jì)算效率低、計(jì)算量大的問題。而SRCNN的方法在訓(xùn)練好一個(gè)區(qū)域的模型之后，可以快速、大量地處理該區(qū)域不同狀態(tài)的SWH數(shù)據(jù)，每次處理過程僅僅是多層矩陣運(yùn)算，大大提高了計(jì)算效率，具有快速、大批量處理數(shù)據(jù)的能力。

表3 海表面有效波高場(chǎng)高分辨率處理效果的不同指標(biāo)

從各項(xiàng)指標(biāo)結(jié)果上看，Kriging算法與兩種SRCNN算法都可以得到很好的高分辨率結(jié)果，且明顯優(yōu)于傳統(tǒng)的Bicubic算法。而SRCNN_2的各個(gè)指標(biāo)都要優(yōu)于SRCNN_1，表明本文對(duì)于SRCNN方法的改造有明顯的提升效果。SRCNN_2在與Kriging方法效果相當(dāng)?shù)那疤嵯?，大大縮短了計(jì)算時(shí)間和計(jì)算量，處理大量的低分辨率數(shù)據(jù)時(shí)，深度學(xué)習(xí)的方法具有更高的計(jì)算效率。

4.2 誤差分布分析

本文計(jì)算了每一點(diǎn)的估計(jì)誤差，并求得各種方法所造成誤差的最大、最小值(見表4)。分析表中數(shù)據(jù)可得，在最大、最小估計(jì)誤差方面，Bicubic算法、Kriging算法和SRCNN_1效果相近，SRCNN_2效果優(yōu)于這三種算法，不容易出現(xiàn)很大的偏差。雖然上文的結(jié)果顯示Kriging整體均方根誤差比較小，但其在某些點(diǎn)的處理效果不佳，會(huì)造成較大的誤差。SRCNN_2相比SRCNN_1有了進(jìn)一步的優(yōu)化，與其他三種方法相比，SRCNN_2的結(jié)果在極端海浪狀況的分析問題上更具有優(yōu)勢(shì)。

表4 海表面有效波高場(chǎng)高分辨率處理誤差的最大、最小值

進(jìn)一步分析誤差的空間分布規(guī)律，圖4和5分別展示了兩次實(shí)驗(yàn)的估計(jì)誤差分布場(chǎng)，分析誤差分布特點(diǎn)可以看出：Kriging方法的整體估計(jì)(見圖5(b))偏高在邊緣處的估計(jì)偏低且誤差大，兩種SRCNN算法的誤差分布場(chǎng)整體趨勢(shì)與Bicubic算法一致，但誤差范圍明顯小于Bicubic算法，SRCNN_2比SRCNN_1的誤差分布更加均勻，Bicubic、SRCNN_1和 SRCNN_2這三種方法的誤差場(chǎng)分布與原始海浪SWH場(chǎng)的分布有關(guān)。

((a)雙三次插值；(b)kriging插值；(c)SRCNN方案一；(d)SRCNN方案二。(a)Bicubic interpolation；(b)Kriging interpolation；(c)SRCNN_1；(d)SRCNN_2.)

((a)雙三次插值；(b)kriging插值；(c)SRCNN_1；(d)SRCNN_2。(a)bicubic interpolation；(b)Kriging interpolation；(c)SRCNN_1；(d)SRCNN_2.)

4.3 近岸海域的測(cè)試

為了進(jìn)一步測(cè)試SRCNN_2方法的效果，選取0.75°S～23.125°S，66.875°W～89.25°W的近岸海域進(jìn)行測(cè)試，測(cè)試效果如圖6所示。相比于雙三次插值，SRCNN_2方法在海岸線的還原上效果較好，但是靠近海岸線部分的值偏小。誤差原因可能如下：在卷積過程中，對(duì)陸地的處理是將其設(shè)置為0值，因此在海陸過渡處的插值偏?。挥?xùn)練集中的數(shù)據(jù)均不包含近岸處的浪場(chǎng)，沒有很好地將近岸的高、低分辨率數(shù)據(jù)的映射關(guān)系學(xué)到網(wǎng)絡(luò)中。在今后的實(shí)驗(yàn)中，在訓(xùn)練集中加入近岸的數(shù)據(jù)，會(huì)得到較好的改善。

((a)低分辨率SWH；(b)原始高分辨率SWH；(c)雙三次插值；(d)SRCNN_2。(a)Low resolution SWH；(b)Original high；resolution SWH；(c)Double cubic interpolation；(d)SRCNN_2.)

4.4 總結(jié)與展望

總結(jié)以上研究得出結(jié)論：整體指標(biāo)排名：Kriging>SRCNN_2>SRCNN_1>Bicubic,計(jì)算效率排名：SRCNN_2>SRCNN_1>Bicubic> Kriging，誤差極值效果排名：SRCNN_2>SRCNN_1>Bicubic> Kriging。

綜上所述，改進(jìn)的SRCNN方法(SRCNN_2)是一種在整體效果、局部細(xì)節(jié)和計(jì)算效率方面均比較優(yōu)秀的高分辨率處理方法，運(yùn)用顯卡進(jìn)行并行運(yùn)算大大節(jié)約了計(jì)算時(shí)間，實(shí)現(xiàn)了用深度學(xué)習(xí)方法直接處理原始海浪數(shù)據(jù)，提高了深度學(xué)習(xí)方法在海洋數(shù)據(jù)處理中的實(shí)際應(yīng)用價(jià)值。相比于傳統(tǒng)的Bicubic和Kriging方法，SRCNN_2的結(jié)果為SWH定性分析、定量分析和異常情況分析等提供了更加可靠的數(shù)據(jù)支持。

在之后的研究中，作者將使用不同的神經(jīng)網(wǎng)絡(luò)模型進(jìn)行高分辨率處理，擴(kuò)大研究的空間范圍和時(shí)間范圍，討論不同方法誤差分布的規(guī)律及產(chǎn)生原因，增加測(cè)試量，在地形復(fù)雜的近岸海域區(qū)域進(jìn)行計(jì)算和實(shí)驗(yàn)，進(jìn)一步提高結(jié)果的說服力。