基于AIS數(shù)據(jù)的船舶航跡多維預(yù)測(cè)方法

高大為， 朱永生， 張金奮， 鄢博冉， 何延康， 閆 柯

(1.西安交通大學(xué) 現(xiàn)代設(shè)計(jì)與轉(zhuǎn)子軸承教育部重點(diǎn)實(shí)驗(yàn)室， 陜西 西安 710049；2.武漢理工大學(xué) 智能交通系統(tǒng)研究中心， 湖北 武漢 430063)

船舶實(shí)時(shí)航行軌跡的預(yù)測(cè)精度不高是造成船舶碰撞、船-橋碰撞等典型水上交通事故的重要原因之一。[1]船舶的軌跡可描述船舶狀態(tài)，也可用于船舶的導(dǎo)航、避碰和海事監(jiān)管，從而提升船舶航行安全性。隨著近年來(lái)數(shù)據(jù)挖掘(Data Mining, DM)技術(shù)的興起，越來(lái)越多的技術(shù)手段被用于船舶軌跡的挖掘與預(yù)測(cè)，試圖基于歷史船舶航跡數(shù)據(jù)來(lái)預(yù)測(cè)未來(lái)一段時(shí)間內(nèi)的船舶航跡，從而感知潛在的航行風(fēng)險(xiǎn)，以便及早采取有效措施避免船舶的碰撞問(wèn)題。

1)在軌跡預(yù)測(cè)問(wèn)題上，DM技術(shù)已在多個(gè)領(lǐng)域得到了應(yīng)用。ZHANG等[2]提出使用矩陣神經(jīng)網(wǎng)絡(luò)來(lái)預(yù)測(cè)氣旋軌跡，避免對(duì)數(shù)據(jù)矢量化所造成的空間信息關(guān)聯(lián)性的缺失，并通過(guò)試驗(yàn)證明其相對(duì)于循環(huán)神經(jīng)網(wǎng)絡(luò)具有更強(qiáng)的優(yōu)勢(shì)；HAO等[3]提出一種基于長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)(Long Short-Term Memory,LSTM)的行人軌跡預(yù)測(cè)方法，不僅考慮社會(huì)鄰域的影響，還考慮了場(chǎng)景布局對(duì)預(yù)測(cè)結(jié)果的影響；LI等[4]基于測(cè)距傳感器和相機(jī)圖像中所學(xué)習(xí)到的三維姿態(tài)數(shù)據(jù)，提出一種三自由度行人姿態(tài)軌跡預(yù)測(cè)方法;ALTCHé等[5]針對(duì)目前的先進(jìn)駕駛輔助系統(tǒng)(Advanced Driving Assistance Systems, ADAS)無(wú)法進(jìn)行中期預(yù)測(cè)的問(wèn)題，提出一種預(yù)測(cè)未來(lái)公路上車輛橫向和縱向軌跡的方法，從而將車輛所能預(yù)測(cè)到的軌跡延長(zhǎng)。

2)在船舶的軌跡跟蹤方面，相關(guān)技術(shù)也日漸成熟。基于船舶自動(dòng)識(shí)別系統(tǒng)(Automatic Identification System, AIS)進(jìn)行船舶軌跡預(yù)測(cè)的研究也越來(lái)越多。然而，由于AIS系統(tǒng)自身的問(wèn)題、不適當(dāng)?shù)氖褂靡约靶盘?hào)傳輸問(wèn)題等原因，數(shù)據(jù)在時(shí)間和空間兩個(gè)維度常存在較大的誤差。針對(duì)該問(wèn)題，ZHAO等[6]對(duì)數(shù)據(jù)集進(jìn)行物理完整性、空間邏輯完整性和時(shí)間準(zhǔn)確度等多重完整性維度的系統(tǒng)性分析，提出一種提高數(shù)據(jù)質(zhì)量的算法;LIU等[7]使用張量分解工具提取了AIS數(shù)據(jù)中固有的運(yùn)動(dòng)模式，并基于鏈路預(yù)測(cè)技術(shù)來(lái)恢復(fù)全數(shù)據(jù)的張量，以解決船舶數(shù)據(jù)的缺失問(wèn)題；徐鐵等[8]提出使用kalman濾波算法來(lái)解決由AIS數(shù)據(jù)不可靠而導(dǎo)致的船舶軌跡不準(zhǔn)確或者誤差較大的問(wèn)題;鄧胡濱等[9]使用kalman濾波器來(lái)消除接收信號(hào)強(qiáng)度指示(Received Signal Strength Indicator, RSSI) 技術(shù)中環(huán)境因素的干擾。然而，基于預(yù)測(cè)算法的AIS數(shù)據(jù)需要保證具有等時(shí)間間隔的位置信息，而這些方法可提高數(shù)據(jù)質(zhì)量，卻難以直接進(jìn)行預(yù)測(cè)。

3)在船舶的航跡預(yù)測(cè)問(wèn)題上， SIMSIR等[10]基于船舶交通管理系統(tǒng)(Vessel Traffic Services, VTS)，使用船舶的位置和速度數(shù)據(jù)訓(xùn)練人工神經(jīng)網(wǎng)絡(luò)，依靠所預(yù)測(cè)的船舶航跡數(shù)據(jù)來(lái)解決狹窄水域的船舶預(yù)警問(wèn)題;GAN等[11]基于大量的船舶軌跡數(shù)據(jù)，利用模糊C均值(Fuzzy C-Means, FCM)聚類算法對(duì)船舶歷史軌跡進(jìn)行聚類，建立船舶速度、載重、自重、最大功率、船長(zhǎng)、船寬、船型和水位等已知因素的關(guān)系模型，解決了較長(zhǎng)距離的船舶軌跡的預(yù)測(cè)問(wèn)題，是一種典型的區(qū)域預(yù)測(cè)模型；權(quán)波等[12]提出使用基于LSTM模型對(duì)AIS數(shù)據(jù)進(jìn)行建模和學(xué)習(xí)，實(shí)現(xiàn)船舶航行軌跡多維度特征表達(dá)，滿足對(duì)船舶軌跡預(yù)測(cè)的精確度和實(shí)時(shí)性的需求；徐婷婷等[13]提出基于BP神經(jīng)網(wǎng)絡(luò)的船舶航跡預(yù)測(cè)方法，將航向和航速作為輸入，經(jīng)度、緯度差作為輸出來(lái)建立航跡預(yù)測(cè)的模型實(shí)現(xiàn)實(shí)時(shí)的航跡預(yù)測(cè)。文獻(xiàn)[12]和文獻(xiàn)[13]屬于典型的當(dāng)前航跡預(yù)測(cè)模型。然而，這些方法一方面需以信噪比較低的AIS數(shù)據(jù)進(jìn)行預(yù)測(cè)，當(dāng)噪聲較大時(shí)，難以實(shí)現(xiàn)理想的效果；另一方面，只使用當(dāng)前航跡數(shù)據(jù)預(yù)測(cè)難以達(dá)到理想的預(yù)測(cè)精度。

針對(duì)以上問(wèn)題，本文提出一種消噪和多維預(yù)測(cè)相結(jié)合的方法來(lái)預(yù)測(cè)船舶航跡。針對(duì)AIS數(shù)據(jù)中經(jīng)、緯度值波動(dòng)性大，有局部增大或減小趨勢(shì)的特點(diǎn)，首先使用三次樣條插值、中值濾波方法對(duì)經(jīng)度、緯度差值進(jìn)行消噪處理；然后對(duì)數(shù)據(jù)進(jìn)行時(shí)間和空間兩個(gè)維度上的建模，利用LSTM進(jìn)行訓(xùn)練并綜合兩方面的信息，最后利用訓(xùn)練后的模型預(yù)測(cè)船舶的航跡信息。

1 LSTM的理論模型

LSTM作為一種捕捉序列化數(shù)據(jù)中動(dòng)態(tài)信息的神經(jīng)網(wǎng)絡(luò)，通過(guò)記憶存儲(chǔ)單元實(shí)現(xiàn)對(duì)序列數(shù)據(jù)的準(zhǔn)確預(yù)測(cè)，已在機(jī)器翻譯等領(lǐng)域有了較為成熟的應(yīng)用。[14-15]LSTM在數(shù)據(jù)預(yù)測(cè)方面有強(qiáng)大的功能，相比于許多預(yù)測(cè)方法可獲得更高的預(yù)測(cè)精度。其體系結(jié)構(gòu)的核心思想是能夠存儲(chǔ)時(shí)變狀態(tài)的存儲(chǔ)單元，其非線性門單元可控制信息的流入和流出。LSTM通過(guò)遺忘門、輸入門和輸出門等3種結(jié)構(gòu)實(shí)現(xiàn)信息保護(hù)和控制，其結(jié)構(gòu)圖見圖1。圖1中：將各個(gè)時(shí)刻的信息輸入到不同的循環(huán)單元中，循環(huán)單元間可傳遞和輸出信息，并由末尾的循環(huán)單元形成最后的輸出。

圖1 LSTM結(jié)構(gòu)圖

圖1中遺忘門、輸入門和輸出門分別為

Ct=σ(Wf×[ht-1,xt]+bf)×Ct-1

(1)

Ct=Ct+σ(Wi×[ht-1,xt]+bi)×

tanh(WC×[ht-1,xt]+bC)

(2)

ht=σ(Wo×[ht-1,xt]+bo)×tanh(Ct)

(3)

式(1)～式(3)中：Ct為當(dāng)前時(shí)刻t的狀態(tài)胞體；ht為胞體的輸出；xt為當(dāng)前胞體的輸入；Wf為各個(gè)門的權(quán)重；bf為所對(duì)應(yīng)的偏差；σ為sigmod激活函數(shù)。

2 基于LSTM的船舶多維預(yù)測(cè)方法

多維預(yù)測(cè)方法包括AIS數(shù)據(jù)預(yù)處理、多維預(yù)測(cè)模型建立和預(yù)測(cè)。其中：AIS數(shù)據(jù)預(yù)處理是對(duì)數(shù)據(jù)進(jìn)行消噪，使其更接近于真實(shí)值；多維預(yù)測(cè)模型則考慮了當(dāng)前航跡數(shù)據(jù)和歷史航跡數(shù)據(jù)，使預(yù)測(cè)結(jié)果更加準(zhǔn)確。

2.1 AIS數(shù)據(jù)預(yù)處理

由于AIS數(shù)據(jù)在時(shí)間和空間兩個(gè)維度存在較大的誤差，甚至有時(shí)存在較明顯的數(shù)據(jù)缺失問(wèn)題，先對(duì)AIS數(shù)據(jù)進(jìn)行預(yù)處理。受到AIS數(shù)據(jù)發(fā)送頻率不一致等因素的影響，接收到的數(shù)據(jù)的時(shí)間間隔在大多數(shù)情況下是不統(tǒng)一的。因此，采用三次樣條插值法來(lái)擬合數(shù)據(jù)。假設(shè)時(shí)間序列區(qū)間[a,b]可分成n個(gè)區(qū)間，即[(x0,x1),(x1,x2),…,(xn-1,xn)]，共有n+1個(gè)數(shù)據(jù)點(diǎn)且對(duì)應(yīng)的函數(shù)值分別為y0,y1,…,yn。若插值函數(shù)S(x)滿足

S(xi)=yi(i=1,2,…,n)

(4)

式(4)在各個(gè)小區(qū)間上都是不高于三次的多項(xiàng)式，假設(shè)為

Si(x)=aix3+bix2+cix+di

(5)

式(5)中：ai、bi、ci和di為待定系數(shù)。若三次樣條方程滿足曲線光滑的條件，即S(x)、S′(x)和S″(x)連續(xù)，則除了根據(jù)每個(gè)節(jié)點(diǎn)處的值建立方程組外，還可根據(jù)該連續(xù)條件建立方程組

S(xi-0)=S(xi+0)

(6)

S′(xi-0)=S′(xi+0)

(7)

S″(xi-0)=S″(xi+0)

(8)

在此基礎(chǔ)上，再加上邊界條件即可求得所有三次樣條方程的系數(shù)，從而構(gòu)造出整個(gè)插值函數(shù)。當(dāng)計(jì)算某個(gè)點(diǎn)的插值時(shí)，則可由其所處的區(qū)間代入方程求得。

由于船舶的速度和航向經(jīng)常發(fā)生變化，經(jīng)緯度值具有較大的波動(dòng)性，使用原始的AIS數(shù)據(jù)進(jìn)行建模預(yù)測(cè)容易造成較大的誤差。因此，在三次樣條插值之后，進(jìn)一步使用中值濾波的方法，分析該時(shí)刻鄰域內(nèi)的多個(gè)數(shù)據(jù)來(lái)計(jì)算在該時(shí)刻的信息。

在對(duì)原始AIS數(shù)據(jù)進(jìn)行插值處理后，接著使用中值濾波來(lái)對(duì)其進(jìn)行去噪。設(shè)AIS數(shù)據(jù)的第i個(gè)點(diǎn)為x(i)，由其周圍2N個(gè)數(shù)據(jù)點(diǎn)組成的數(shù)據(jù)集合為

[x(i-N),…,x(i),…,x(i+N)]

(9)

則進(jìn)行中值濾波后該點(diǎn)的值為

x′(i)=Med[x(i-N),…,x(i),…,x(i+N)]

(10)

式(10)中：Med為該數(shù)據(jù)集合的中間值。

中值濾波的優(yōu)點(diǎn)在于對(duì)如AIS數(shù)據(jù)這種單位時(shí)間內(nèi)的位移波動(dòng)較大的數(shù)據(jù)，可起到較好的消噪效果。

2.2 多維預(yù)測(cè)模型

本文將船舶利用當(dāng)前運(yùn)動(dòng)航跡進(jìn)行預(yù)測(cè)時(shí)的模型稱為“當(dāng)前航跡模型(Current Trajectory Prediction Model, CTPM)”。CTPM已在多篇文獻(xiàn)有所介紹。[12]該模型由船舶當(dāng)前航行的數(shù)據(jù)建模，可反映船舶的實(shí)時(shí)狀態(tài)信息、自然環(huán)境狀態(tài)信息等。然而，該模型未能考慮船舶的空間位置特點(diǎn)。而船舶的航行狀態(tài)與所在水域的空間位置具有較大的關(guān)聯(lián)性。例如:在較窄的水域內(nèi)，船舶的航行速度會(huì)較慢;在較寬的水域內(nèi)，船舶的航行速度則會(huì)較快；在船舶接近港口時(shí)，速度會(huì)較慢；在彎曲航道，船舶一般會(huì)采取轉(zhuǎn)向操作。因此，在CTPM的基礎(chǔ)上，建立“區(qū)域預(yù)測(cè)模型(Regional Prediction Model, RPM)”。

RPM是指船舶利用預(yù)測(cè)時(shí)刻所處區(qū)域的位置信息建立模型，通過(guò)該區(qū)域的歷史軌跡來(lái)分析所預(yù)測(cè)時(shí)刻的航行軌跡。RPM與CTPM是一種互補(bǔ)關(guān)系，通過(guò)將兩者結(jié)合來(lái)提高對(duì)船舶航跡的預(yù)測(cè)能力。CTPM提取當(dāng)前船舶的航行狀態(tài)特征；RPM提取空間位置的水域特點(diǎn)對(duì)船舶航行狀態(tài)所影響的特征。

由于在短時(shí)間內(nèi)，相鄰AIS數(shù)據(jù)間的經(jīng)緯度差值一般較小(10-5到10-4甚至更小的范圍內(nèi))，而該數(shù)值的量級(jí)會(huì)給模型的訓(xùn)練帶來(lái)較多的困難。因此，選用經(jīng)緯度差值來(lái)作為預(yù)測(cè)模型在單個(gè)時(shí)刻的輸入和輸出，同時(shí)，將航速、航向信息作為模型的考慮因素，見圖2。圖2中：對(duì)于CTPM來(lái)說(shuō)，使用相鄰數(shù)據(jù)點(diǎn)的差值作為輸入；而對(duì)于RPM而言，則將其中一個(gè)航跡設(shè)為參考數(shù)據(jù)，利用與其差值作為輸入。

圖2 LSTM模型在t時(shí)刻的輸入與輸出

對(duì)RPM來(lái)說(shuō)，設(shè){Xi,j|i∈m,j∈n} 為船舶的某區(qū)域歷史航跡，其共有m段航跡，n個(gè)位置。

Xi,j=[xi,j,yi,j,vi,j,dij]

(11)

式(11)中：xi,j、yi,j、vi,j和dij分別為船舶在某個(gè)位置的經(jīng)度、緯度、速度和航向。為避免網(wǎng)絡(luò)訓(xùn)練中的數(shù)值問(wèn)題，并使網(wǎng)絡(luò)快速收斂，先對(duì)各個(gè)航跡數(shù)據(jù)進(jìn)行歸一化操作為

(12)

(13)

i∈{1,2,…,m})

(14)

同時(shí)，以當(dāng)前航跡數(shù)據(jù)來(lái)訓(xùn)練CTPM為

(15)

式(15)中：

(16)

再用全連接層將二者融合來(lái)進(jìn)行訓(xùn)練，最終的預(yù)測(cè)結(jié)果為

(17)

使用均方誤差(Mean Squared Error, MSE)作為損失函數(shù)訓(xùn)練模型為

(18)

模型使用自適應(yīng)矩估計(jì)優(yōu)化器(Adam)對(duì)其進(jìn)行優(yōu)化。Adam是一種對(duì)隨機(jī)目標(biāo)函數(shù)執(zhí)行一階梯度優(yōu)化的算法。模型隱藏層中的單元數(shù)通過(guò)在{85, 170, 255, 340, 425} 中搜尋最優(yōu)的值來(lái)作為此處的參數(shù)，用于記憶和存儲(chǔ)過(guò)去的狀態(tài)。最大訓(xùn)練周期數(shù)的選取需要保證模型可以收斂，在此根據(jù)模型多次試驗(yàn)結(jié)果選擇為1 500，循環(huán)層個(gè)數(shù)設(shè)置為3，預(yù)測(cè)步數(shù)設(shè)置為1。使用訓(xùn)練好的模型預(yù)測(cè)數(shù)據(jù)后，再將數(shù)據(jù)還原成經(jīng)緯度值。船舶航跡預(yù)測(cè)模型見圖3。

圖3 船舶航跡預(yù)測(cè)模型

2.3 船舶航跡預(yù)測(cè)流程

對(duì)于船舶的航跡預(yù)測(cè)來(lái)說(shuō)，一方面，前一時(shí)刻的經(jīng)緯度差值位置信息決定著待預(yù)測(cè)時(shí)刻的經(jīng)緯度值位置信息；另一方面，船舶的航速、航向也對(duì)經(jīng)緯度差值存在著較大的影響。因此，以經(jīng)緯度差值、航速、航向等4個(gè)參數(shù)作為網(wǎng)絡(luò)的輸入，以經(jīng)、緯度差值作為輸出。方法流程見圖4。

圖4 所提出方法流程

1)對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，對(duì)其進(jìn)行初步判定、插值處理等，從而獲得等時(shí)間間隔的經(jīng)緯度、航向和速度信息。

2)選取指定間隔的相鄰數(shù)據(jù)點(diǎn)的經(jīng)緯度差值。

3)通過(guò)中值濾波方法進(jìn)行濾波。

4)基于船舶空間歷史數(shù)據(jù)，構(gòu)造區(qū)域范圍內(nèi)的船舶軌跡預(yù)測(cè)模型RPM。

5)基于船舶的實(shí)時(shí)數(shù)據(jù)，構(gòu)造船舶當(dāng)前航跡預(yù)測(cè)模型CTPM。

6)聯(lián)合對(duì)以上模型進(jìn)行訓(xùn)練。

7)獲得實(shí)時(shí)預(yù)測(cè)的經(jīng)緯度差值。

8)將其還原成經(jīng)緯度值。

3 試驗(yàn)數(shù)據(jù)分析

為測(cè)試提出的軌跡預(yù)測(cè)算法的性能，分別使用海上AIS數(shù)據(jù)和內(nèi)河AIS數(shù)據(jù)來(lái)訓(xùn)練模型并預(yù)測(cè)。

3.1 案例1：海上航跡預(yù)測(cè)

海上AIS數(shù)據(jù)為大連、煙臺(tái)2個(gè)港口之間的數(shù)據(jù)，該水域的經(jīng)緯度值見圖5。由圖5可知：在相同任務(wù)的航行中每艘船舶的航行軌跡在大體趨勢(shì)上存在較大的相關(guān)性，但在各個(gè)位置處仍然存在著一定的差別。本文選取其中一段數(shù)據(jù)來(lái)驗(yàn)證所提出方法的效果。

圖5 某個(gè)范圍內(nèi)的船舶經(jīng)緯度值

首先對(duì)AIS數(shù)據(jù)進(jìn)行預(yù)處理，根據(jù)初步判斷法則去除明顯的異常數(shù)據(jù)。然后使用三次樣條插值，使其在時(shí)間上等間隔。對(duì)預(yù)處理后的經(jīng)緯度數(shù)據(jù)進(jìn)行差值處理，即獲得等間隔時(shí)間段內(nèi)船舶的航行里程信息，然后對(duì)經(jīng)緯度差值進(jìn)行消噪處理，結(jié)果見圖6。由圖6可知：消噪前的經(jīng)緯度差值數(shù)據(jù)(圖6中灰色部分)存在以下特點(diǎn)：在經(jīng)緯度差值較為穩(wěn)定的階段，其以某個(gè)數(shù)據(jù)值為中心，存在“正弦式”波動(dòng)的規(guī)律；在某些階段其經(jīng)緯度差會(huì)出現(xiàn)局部的增大或減小。對(duì)于前者，對(duì)該數(shù)據(jù)的消噪可使用均值濾波法，然而考慮到后者中出現(xiàn)的局部增大或減小，使用中值濾波可同時(shí)滿足以上兩方面消噪的需求。消噪后的數(shù)據(jù)如圖6中黑色曲線所示，經(jīng)中值濾波消噪后數(shù)據(jù)的波動(dòng)性得到較好的解決，且能保留局部出現(xiàn)的增大或減小的趨勢(shì)。

a) 經(jīng)度差

選取經(jīng)緯度從(121.504 9°E, 37.771 0°N)到(121.585 3°E, 37.983 1°N)的一段數(shù)據(jù)進(jìn)行分析，通過(guò)尋優(yōu)設(shè)置隱含層數(shù)為170，預(yù)測(cè)結(jié)果見圖7。

a) 整體預(yù)測(cè)結(jié)果

圖7中：實(shí)際值為實(shí)際的航行軌跡；訓(xùn)練值為訓(xùn)練網(wǎng)絡(luò)所使用的的經(jīng)緯度差值部分；測(cè)試值為進(jìn)行預(yù)測(cè)的部分。由圖7可知：船舶的航行軌跡預(yù)測(cè)精度較高，與實(shí)際值的趨勢(shì)一致。圖7b結(jié)果為在圖7a的基礎(chǔ)上進(jìn)行局部放大的某段訓(xùn)練數(shù)據(jù)，圖7c和圖7d為測(cè)試數(shù)據(jù)上船舶某兩段輕微轉(zhuǎn)彎時(shí)的軌跡，預(yù)測(cè)結(jié)果可看出，雖然測(cè)試值與實(shí)際值出現(xiàn)了一定的偏差，但是誤差較小，預(yù)測(cè)效果良好。

預(yù)測(cè)誤差可更為直觀的反映模型的性能，預(yù)測(cè)得到的經(jīng)緯度誤差見圖8。再結(jié)合該案例中的時(shí)刻差(6 s)可知:二者的均值分別為6.23×10-5°和8.89×10-5°。同時(shí)，本文將該方法與傳統(tǒng)的LSTM方法進(jìn)行對(duì)比見表1。由表1可知：在經(jīng)度、緯度兩個(gè)方向上，所提出的方法均能獲得較小的誤差，證實(shí)所提出方法的有效性。

a) 經(jīng)度誤差

表1 兩種預(yù)測(cè)方法的平均誤差的對(duì)比

3.2 案例2：內(nèi)河航跡預(yù)測(cè)

為進(jìn)一步測(cè)試提出的軌跡預(yù)測(cè)算法的性能，使用某段內(nèi)河的AIS數(shù)據(jù)進(jìn)行分析。所研究水域的經(jīng)緯度數(shù)據(jù)見圖9。由圖9可知：在相同任務(wù)的航行中各個(gè)船舶的航行軌跡在大體趨勢(shì)上存在較大的相關(guān)性，但在各個(gè)位置仍存在著一定的差別。本文選取船舶直線行駛、船舶轉(zhuǎn)彎、船舶加速和船舶減速時(shí)的AIS數(shù)據(jù)等4段典型數(shù)據(jù)，來(lái)驗(yàn)證所提出方法的效果。該4段數(shù)據(jù)以不同的標(biāo)記展示在圖9中。

圖9 某個(gè)范圍內(nèi)的船舶經(jīng)緯度值

在進(jìn)行模型預(yù)測(cè)之前，首先使用上述數(shù)據(jù)預(yù)處理方法對(duì)航跡數(shù)據(jù)進(jìn)行三次樣條插值和中值濾波，隨后再進(jìn)行航跡預(yù)測(cè)。通過(guò)尋優(yōu)設(shè)置隱含層數(shù)為170，預(yù)測(cè)結(jié)果見圖10～圖13。由圖10～圖13可知：該預(yù)測(cè)方法在船舶的各個(gè)階段都實(shí)現(xiàn)了較為準(zhǔn)確的預(yù)測(cè)。

在直線行駛階段如圖12所示，船船航跡的預(yù)測(cè)誤差最小，其最大預(yù)測(cè)誤差低至7.0×10-7°。在圖13所示的轉(zhuǎn)彎階段，船舶的預(yù)測(cè)誤差相對(duì)較大些，這是由于訓(xùn)練數(shù)據(jù)中沒(méi)有學(xué)習(xí)到船舶轉(zhuǎn)彎時(shí)的狀態(tài)信息。同時(shí)，在船舶加減速階段如圖10和圖11所示，其預(yù)測(cè)誤差相比直線階段也要略大一些。然而，相比于所預(yù)測(cè)時(shí)間內(nèi)的船舶位移，該誤差在可接受的范圍之內(nèi)。為驗(yàn)證方法的有效性，同時(shí)，將該方法與傳統(tǒng)的LSTM方法進(jìn)行對(duì)比，結(jié)果見表2。由表2可知：雖然直線階段二者差別不大，但是在其他的航行狀態(tài)下，所提出的方法的預(yù)測(cè)誤差明顯較小。通過(guò)以上分析，說(shuō)明所提出的預(yù)測(cè)模型在內(nèi)河AIS數(shù)據(jù)中也可達(dá)到較好的效果。

a) 預(yù)測(cè)結(jié)果

a) 預(yù)測(cè)結(jié)果

a) 預(yù)測(cè)結(jié)果

a) 預(yù)測(cè)結(jié)果

表2 兩種預(yù)測(cè)方法的平均均方根誤差的對(duì)比

4 結(jié)束語(yǔ)

本文提出一種基于船舶AIS數(shù)據(jù)的航跡多維預(yù)測(cè)方法，首先通過(guò)3次樣條插值、中值濾波等技術(shù)對(duì)數(shù)據(jù)進(jìn)行預(yù)處理，然后在此基礎(chǔ)上基于歷史軌跡數(shù)據(jù)和當(dāng)前航跡數(shù)據(jù)構(gòu)建多維預(yù)測(cè)算法來(lái)提高預(yù)測(cè)的精度。在數(shù)據(jù)預(yù)處理環(huán)節(jié)，通過(guò)多種方法保證數(shù)據(jù)的可靠性；在航跡多維預(yù)測(cè)環(huán)節(jié)，分別建立區(qū)域預(yù)測(cè)模型和當(dāng)前航跡預(yù)測(cè)模型，從時(shí)間和空間兩個(gè)維度同時(shí)對(duì)航跡進(jìn)行預(yù)測(cè)并綜合其結(jié)果。通過(guò)海上AIS數(shù)據(jù)和內(nèi)河AIS數(shù)據(jù)的試驗(yàn)驗(yàn)證，結(jié)果表明所提出的方法可實(shí)現(xiàn)對(duì)船舶航跡較為準(zhǔn)確的預(yù)測(cè)。