惠 飛,彭 娜,景首才,周 琪,賈 碩

(長安大學(xué) 信息工程學(xué)院,西安 710064)

0 概述

異常駕駛行為通常指行車速度超過速度限制、駕駛員突然改變速度或不斷地改變車輛側(cè)向位置等非平穩(wěn)行為,此類行為是導(dǎo)致交通事故的重要誘因。因此,為交通管理部門提供可靠的異常駕駛檢測方法具有非常重要的實(shí)際意義。

目前異常駕駛行為檢測方法主要基于視頻數(shù)據(jù)和車輛GPS軌跡數(shù)據(jù)進(jìn)行研究。文獻(xiàn)[1]利用視頻數(shù)據(jù)對車輛的軌跡和速度進(jìn)行分析,在背景差分法的基礎(chǔ)上通過建立跟蹤車輛信息鏈實(shí)現(xiàn)對車輛跟蹤和異常行為檢測;文獻(xiàn)[2]利用背景差分法和均值漂移算法獲取車輛位置、速度、方向等判別參數(shù),對3種判別參數(shù)的狀態(tài)函數(shù)加權(quán)融合檢測車輛異常行為;文獻(xiàn)[3]利用GPS提供的位置、時間、速度信息,以車輛正常行駛情況下的縱向和橫向加速度為閾值,檢測和識別駕駛員異常駕駛行為;文獻(xiàn)[4]以大量重特大交通事故數(shù)據(jù)為基礎(chǔ),通過聚類分析從危險駕駛行為角度研究外部影響因素與駕駛行為的關(guān)系。筆者通過對文獻(xiàn)的分析發(fā)現(xiàn),目前基于視頻的異常行為檢測,主要針對單車受限場景,難以對運(yùn)輸全過程進(jìn)行監(jiān)控[5]。而對GPS的軌跡分析,也以對單車進(jìn)行先驗閾值判斷為主,大多是設(shè)定檢測閾值而觸發(fā)報警,缺乏數(shù)據(jù)深層次分析與信息挖掘步驟,不具有智能化辨識與特性分析功能。

為此,本文提出一種基于GPS數(shù)據(jù)的駕駛行為異常檢測方法。首先根據(jù)時間、速度、加速度、方向、轉(zhuǎn)角等全局與局部特征及其對應(yīng)的13個統(tǒng)計量,構(gòu)建車輛駕駛行為的特征屬性,然后利用基于結(jié)構(gòu)相似度度量的凝聚層次聚類方法對駕駛行為進(jìn)行聚類建模,進(jìn)而實(shí)現(xiàn)異常行為檢測。本文方法將測試數(shù)據(jù)與訓(xùn)練數(shù)據(jù)集獲得的模型庫進(jìn)行對比檢測,無需使用經(jīng)驗指標(biāo)來設(shè)定各類異常的閾值,從而有效避免人為因素對檢測結(jié)果的影響。

1 駕駛行為聚類與異常檢測方法框架

本文對商用車司機(jī)異常駕駛行為進(jìn)行檢測,主要檢測超速、急加速、急減速、違章變道等行為。首先基于文獻(xiàn)[6]提出的全局和局部運(yùn)動特征分類方法,總結(jié)關(guān)于車輛運(yùn)行GPS數(shù)據(jù)所體現(xiàn)的司機(jī)駕駛行為特征,從GPS數(shù)據(jù)中提取運(yùn)動特征及其統(tǒng)計量,構(gòu)建車輛多特征駕駛行為模型,提出基于多特征的駕駛行為聚類與異常檢測方法,方法框架如圖1所示。其中數(shù)據(jù)準(zhǔn)備與預(yù)處理為大數(shù)據(jù)剔除無關(guān)噪聲數(shù)據(jù)基本環(huán)節(jié),本文不再贅述。

圖1 基于多特征駕駛行為聚類的異常檢測方法框架

1.1 駕駛行為多特征提取

1.1.1 全局特征提取

對于位置點(diǎn)數(shù)據(jù)全局特征的提取,首先根據(jù)軌跡數(shù)據(jù)點(diǎn)的時間信息,提取出數(shù)據(jù)點(diǎn)發(fā)生的17個全局時間特征,即:“月”(Month,如1月);“日”(Day,如2日);“星期”(Week,星期日到星期六分別為0～6);“周末/周內(nèi)”(Workday,其屬性為0/1屬性,周內(nèi)為1,周末為0);“節(jié)假日”(Holiday,其屬性為0/1屬性,法定節(jié)假日為1,法定工作日為0);“時段”(Time,本文將1 d 24 h等間隔劃分為12個one-hot[7-8]時段特征,值為0/1)。在提取全局時間特征后,根據(jù)數(shù)據(jù)點(diǎn)信息提取其他運(yùn)動參數(shù)。

軌跡是一系列連續(xù)GPS點(diǎn)的集合,記T=。其中,每個位置點(diǎn)pi包含位置信息(xi,yi)及時間信息(ti),記pi=。根據(jù)GPS數(shù)據(jù)信息提取兩位置點(diǎn)間的距離di、時間間隔Δti、速度vi、加速度ai、方向directioni、轉(zhuǎn)角θi以及相關(guān)統(tǒng)計量(均值、中值、標(biāo)準(zhǔn)差等)。然后將針對車輛駕駛行為中典型的異常駕駛行為(超速、急加/減速、違章變道)分別進(jìn)行特征提取。

1.1.2 超速駕駛行為特征提取

通過對兩數(shù)據(jù)點(diǎn)間的行駛距離、速度和時間間隔分析,能夠?qū)崿F(xiàn)對超速行為的識別。統(tǒng)計相關(guān)超速信息,通過對超速持續(xù)時間以及超速行為之間的間隔距離來判斷是否為違章超速行為,而非如超車等正常超速行為[9]。為判斷超速行為是否為違章超速行為,引入違章超速行駛行為距離判斷閾值ɑ與界定時長δ,按照行業(yè)判定經(jīng)驗取:

a≥500 m,δ≈30 s

(1)

根據(jù)超速行為行駛距離di引入“違章超速”特征值k,k=1表示違章超速,k=0表示未違章超速:

(2)

將市區(qū)內(nèi)行車速度按其大小值分為4個等級:低速(v≤20 km/h),中速(20 km/h60 km/h),本文根據(jù)司機(jī)駕駛速度等級提取“低速”“中速”“高速” “超速”4個特征(皆為one-hot特征,即值為0/1,0為否,1為是)。

1.1.3 急加/減速駕駛行為特征提取

在交通運(yùn)輸行業(yè)相關(guān)法規(guī)中,急加/減速行為并不屬于違法違規(guī)行為,但經(jīng)常性的急加/減速行為對制動系統(tǒng)及系統(tǒng)部件都有損害,并且急加/減速行為發(fā)生次數(shù)也是判斷司機(jī)開車是否平穩(wěn)的標(biāo)準(zhǔn)之一,對駕駛安全也有一定程度的影響。因此,在全局運(yùn)動特征提取的基礎(chǔ)上,根據(jù)點(diǎn)pi、pi+1間速度vi、vi+1和時間間隔Δti計算點(diǎn)pi+1的加速度值,Δti要求不大于2 s,如式(3)所示。若Δti過長,則無意義。

(3)

查閱司機(jī)駕駛行為“急加/減速”判定標(biāo)準(zhǔn),將加速度按其大小值分為4個標(biāo)準(zhǔn)等級:一級標(biāo)準(zhǔn)(ai≥2.78 m/s2,持續(xù)時間≥2 s),二級標(biāo)準(zhǔn)(2.22 m/s2≤a<2.78 m/s2,持續(xù)時間≥2 s),三級標(biāo)準(zhǔn)(1.67 m/s2≤a<2.22 m/s2，持續(xù)時間≥2 s),四級標(biāo)準(zhǔn)(ai<1.67 m/s2,持續(xù)時間≥2 s)。因此,本文根據(jù)加速度等級增加4個one-hot特征。

根據(jù)2個數(shù)據(jù)點(diǎn)pi與pi+1速度差值的正負(fù)判斷位置點(diǎn)處于加速或是減速,提取“加速/減速”特征,并根據(jù)2個數(shù)據(jù)點(diǎn)pi與pi+1之間的時間間隔Δti,增加“加速度持續(xù)時間”特征。

1.1.4 變道駕駛行為特征提取

在變道過程中,高速行駛狀態(tài)下的車輛轉(zhuǎn)向角度是很小的。因此,在提取出數(shù)據(jù)點(diǎn)全局特征的基礎(chǔ)上,還需計算出2個數(shù)據(jù)點(diǎn)之間行駛路線的彎曲程度和轉(zhuǎn)角。其中曲率si為兩位置點(diǎn)pi與pi+1之間的移動距離與其直線距離之比,表示兩點(diǎn)之間路徑的彎曲度,如式(4)所示。

si=(dist(pi-1,pi)+dist(pi,pi+1))/dist(pi-1,pi+1)

(4)

位置點(diǎn)的方向由數(shù)據(jù)集提供,表示每個位置點(diǎn)在當(dāng)前位置從正北方向順時針旋轉(zhuǎn)的角度Δθi。轉(zhuǎn)角Δθi指連續(xù)2個位置點(diǎn)間方向轉(zhuǎn)動角度,根據(jù)兩點(diǎn)各自的方向信息θi與θi+1進(jìn)行差值計算所得,如式(5)所示。

Δθi=|θi+1-θi|

(5)

兩位置點(diǎn)的方向與轉(zhuǎn)角示例圖如圖2所示。

圖2 位置點(diǎn)P1與P2的方向與轉(zhuǎn)角示例

根據(jù)一條軌跡為t=中所含軌跡點(diǎn)數(shù)n,統(tǒng)計其在一條行駛軌跡中車輛變道次數(shù)m并增加“變道次數(shù)”特征。對于是否未違章變道行為,將根據(jù)車輛連續(xù)變道次數(shù)m引入“違章變道”特征。根據(jù)行業(yè)經(jīng)驗,設(shè)定連續(xù)變道次數(shù)大于3次判定為違章變道。

1.1.5 各特征統(tǒng)計量計算

在已有的距離、速度、加速度、曲率、轉(zhuǎn)角5個全局運(yùn)動特征的基礎(chǔ)上,分別計算其對應(yīng)的13個數(shù)學(xué)統(tǒng)計量:均值,中值,標(biāo)準(zhǔn)差,變異系數(shù),最大3個數(shù),最小3個數(shù),自相關(guān)系數(shù),偏度系數(shù)和峰度系數(shù)。其中,用變異系數(shù)來衡量數(shù)據(jù)的離散程度。計算最大/最小的3個數(shù)[10]是為了考慮定位系統(tǒng)獲取的GPS數(shù)據(jù)存在定位精度和誤差的因素。偏度系數(shù)和峰度系數(shù)用來分別衡量統(tǒng)計分布的左右對稱性和形狀的陡峭程度[11]。自相關(guān)系數(shù)是用來度量2個位置點(diǎn)在時間上的相關(guān)程度[12]。一條軌跡的時間序列可表示為xt∈{x1,x2,…,xn},為計算其自相關(guān)系數(shù),可先用式(6)計算自協(xié)方差系數(shù){ck},其中,k為時滯長度,N為序列長度。再利用式(7)計算時滯為1時的自相關(guān)系數(shù)。

(6)

r1=c1/c0

(7)

偏度系數(shù)與峰度系數(shù)的計算可參考文獻(xiàn)[11]。最后,對5個全局運(yùn)動特征分別計算其對應(yīng)13個統(tǒng)計量。圖3顯示了部分運(yùn)動特征及其統(tǒng)計量結(jié)果。

圖3 部分運(yùn)動特征及其統(tǒng)計量

1.2 基于結(jié)構(gòu)距離的相似度度量方式

當(dāng)檢測目標(biāo)為運(yùn)動位置點(diǎn)表現(xiàn)出的駕駛行為的異常情況時,不能直接采用常用的軌跡之間的相似度度量方式[13-14]。本文借鑒圖像質(zhì)量評價中的結(jié)構(gòu)相似度(Structural Similarity,SSIM)[15]思想,采用SSIM方法來度量數(shù)據(jù)點(diǎn)駕駛行為特征之間的距離。

結(jié)構(gòu)距離(Structure Distance,SDist)指對軌跡結(jié)構(gòu)各個特征屬性分別采取不同的距離計算方式來計算,最后結(jié)合為一個距離計算標(biāo)準(zhǔn)。本文數(shù)據(jù)點(diǎn)包含6個全局特征,其軌跡結(jié)構(gòu)為pi=,各異常行為特征中還包含一些one-hot特征。因此,數(shù)據(jù)點(diǎn)結(jié)構(gòu)距離共包含7個部分:時間距離timeDist,位置距離locDist,速度距離speedDist,加速度距離accDist,曲率距離curDist,轉(zhuǎn)角距離angleDist以及one-hot特征距離Dist。

定義W={Wt,Wv,Wd,Wa,Ws,Wθ,Wo}為各部分特征權(quán)重向量,分別對應(yīng)軌跡的特征向量。通過這7個部分的距離構(gòu)成了軌跡相似度的計算,如式(8)所示。其中,時間距離采用余弦相似度度量方式計算。位置距離采用歐氏距離計算。速度和加速度均采取各統(tǒng)計量和的平均。曲率和轉(zhuǎn)角均采用差的絕對值衡量。

SDist(pi,pj)=Wa×timeDist+Wv×speedDist+

Wd×locDis+Wa×accDist+

Ws×curDist+Wθ×

angleDist+Wo×Dist

(8)

由于軌跡中每個特征的值域不同,需要對結(jié)構(gòu)距離進(jìn)行歸一化處理,如式(9)所示,其中Normalized()為距離的歸一化函數(shù)。因為兩點(diǎn)間的結(jié)構(gòu)相似度是對稱的,所以SSIM(pi,pj)=SSIM(pj,pi)。

SSIM(pi,pj)=1-Normalized(SDist(pi,pj))

(9)

下文將針對超速、急加/減速、急剎車、違章變道等典型異常駕駛行為,在已有的全局和局部特征以及相似度度量方式的基礎(chǔ)上,分別進(jìn)行凝聚層次聚類與異常檢測。

2 基于多特征的聚類與異常檢測

2.1 基于多特征的凝聚層次聚類算法

聚類是將抽象的數(shù)據(jù)集劃分成由相似數(shù)據(jù)構(gòu)成的若干個對象集或者簇,使數(shù)據(jù)對象在同一簇中相似度越大,類間相似度越小,聚類效果就越好。聚類分析是異常行為檢測的基礎(chǔ),在傳統(tǒng)的基于聚類的異常點(diǎn)檢測方法中,對象是否被檢測為離群點(diǎn)很大程度依賴于簇的個數(shù)[16]。聚類算法中有部分算法需設(shè)定初始聚類中心和聚類數(shù)目,不能自動確定簇的個數(shù),聚類結(jié)果不穩(wěn)定。

凝聚層次聚類算法具有聚類準(zhǔn)確率高、聚類結(jié)果穩(wěn)定的優(yōu)點(diǎn),無需選取初始聚類中心和設(shè)定任何參數(shù)。但其計算復(fù)雜度大,為了降低復(fù)雜度,自動確定聚類數(shù)目,本文在層次聚類方法中引入拉普拉斯特征映射思想自動確定聚類簇數(shù)。

2.1.1 拉普拉斯特征映射

步驟1通過特征相似矩陣Sn×n構(gòu)造拉普拉斯特征矩陣Ln×n=D-1/2Sn×nD-1/2,其中矩陣D為由矩陣S獲得的對角矩陣。

步驟2對Ln×n矩陣進(jìn)行特征值分解并按降序排列λ1≥λ2≥…≥λn;計算相鄰特征值之差,若第i個特征值和第i+1個特征值之間差異最大,則確定聚類個數(shù)k=argmax|λi+1-λi|。

步驟3構(gòu)造m×k矩陣Lm×k=[l1,l2,…,lk],其中l(wèi)k為矩陣Ln×n的前k個特征值對應(yīng)的特征向量;然后對Lm×k矩陣進(jìn)行歸一化處理,得到矩陣X=[x1,x2,…,xm]后再對低維數(shù)據(jù)xi進(jìn)行后續(xù)聚類。

2.1.2 凝聚層次聚類

在對特征矩陣進(jìn)行拉普拉斯變換后,根據(jù)降維后的特征矩陣及聚類數(shù)進(jìn)行最后的層次聚類。基于凝聚的層次聚類算法的偽代碼如下:

算法CURE層次聚類

輸入數(shù)據(jù)集D={p1,p2,…,pn}

輸出k個聚類簇及質(zhì)心

1.3.3 患兒治療依從性評價 由于哮喘患兒是特殊群體，其對于治療的依從性是針對于家屬，臨床藥師對實(shí)驗組患兒家屬在1年時間內(nèi)，每個月評估其治療的依從性，以能夠穩(wěn)定持續(xù)用藥為較好，以分值5分表示；以雖未能持續(xù)穩(wěn)定進(jìn)行藥物治療，但對于疾病控制仍存在相應(yīng)的重視及治療為一般，以分值3-4分表示；對于用藥治療頻率較低，完全不能按照療程進(jìn)行治療為較差，以分值1-2分表示。

1.對數(shù)據(jù)集建立點(diǎn)與點(diǎn)之間的距離矩陣Sn×n。

2.根據(jù)上述步驟1～步驟3,得到低維數(shù)據(jù)矩陣X及其聚類個數(shù)k。

3.loop

4.將每個對象歸為一類,得到n類,每類包含一個對象。

5.將距離矩陣Sn×n中的值升序排列,選取最小距離的2個類合并為一個類。

6.重新計算新類與舊類之間的距離。

7.重復(fù)前三步,直到最后合并成一個類或達(dá)到設(shè)定條件則退出循環(huán)。

8.輸出k個聚類簇及其質(zhì)心。

2.2 基于類標(biāo)記的異常檢測算法

正常駕駛行為有較高的重復(fù)性和相似性特點(diǎn),而異常駕駛行為相較正常駕駛行為可能因不同場景和交通狀況而表現(xiàn)出不同,與正常模式匹配存在很大的差距。在聚類結(jié)果中,有些簇的規(guī)模可能很小,在無法區(qū)分噪音和離群點(diǎn)下,本文將這些簇標(biāo)記為異常。由上一節(jié)的聚類簇結(jié)果,在設(shè)定參數(shù)α的條件下,根據(jù)閾值標(biāo)記簇為正?；虍惓４?再將待測點(diǎn)與簇進(jìn)行對比檢測。

2.2.1 簇結(jié)果標(biāo)記

文獻(xiàn)[18]定義了大小簇的概念,假設(shè)C={c1,c2,…,ck}為簇集,簇成員個數(shù)記為|ci|,將其按成員數(shù)目降序排序得到{|c1|>|c2|>…>|ck|};給定參數(shù)α,當(dāng)|c1|+|c2|+…+|cb|≥αn(n為原始數(shù)據(jù)集數(shù)目)時,則定義b(1≤b≤k-1)為大簇LC與小簇SC的邊界。那么,可以將簇標(biāo)記為大簇LC={ci|i≤b}與小簇SC={ci|i≥b},并將ci(i≤b)標(biāo)記為正常簇,ci(i≥b)標(biāo)記為異簇。計算大小簇邊界的算法偽代碼如下:

算法計算大小簇邊界

輸入聚類簇及簇內(nèi)數(shù)據(jù),參數(shù)α

輸出簇邊界b

1.計算各簇ci的成員個數(shù)|ci|(i=1,2,…,k),并對|ci|按降序排序,設(shè)置參數(shù)α

2.for(i=1;i≤k;i++)

5.記錄大小簇邊界b并跳出循環(huán)

6.end for

2.2.2 異常檢測算法

假設(shè)某車輛的一條軌跡數(shù)據(jù)為t={p1,p2,…,pm},m指該軌跡長度。計算每個點(diǎn)pi與所有簇之間的距離dist(pi,cj)并找到距離最小的dist(pi,cmin)。若dist(pi,cmin)大于簇半徑R,則將該位置點(diǎn)標(biāo)記為未知屬性點(diǎn),重新計算簇半徑,將R′記為新的簇半徑;否則找出簇cmin所屬簇的標(biāo)記,若標(biāo)記為正常則點(diǎn)pi為正常點(diǎn);若標(biāo)記為異常則點(diǎn)pi為異常點(diǎn)。

3 實(shí)驗與結(jié)果分析

3.1 實(shí)驗環(huán)境設(shè)置與數(shù)據(jù)準(zhǔn)備

實(shí)驗環(huán)境:Windows10操作系統(tǒng),Intel Xeon 2.50 GHz CPU,16.0 GB內(nèi)存。開發(fā)環(huán)境:Microsoft Visual Studio 2008,Matlab12,Anaconda2.0。采用Matlab統(tǒng)計學(xué)工具箱計算方向和轉(zhuǎn)角的統(tǒng)計量,其他特征采用python計算,并采用python語言對數(shù)據(jù)點(diǎn)進(jìn)行聚類分類及后續(xù)異常檢測。

實(shí)驗數(shù)據(jù)采用真實(shí)的西安市2016年的商用車運(yùn)行GPS數(shù)據(jù),包括時間、經(jīng)度、緯度、速度、方向等信息。實(shí)驗選取經(jīng)度在lon∈[108.7,109.2]、緯度在lat∈[34.1,34.5]范圍的西安市區(qū)內(nèi)數(shù)據(jù),共包含59 999個軌跡點(diǎn)即4 953條軌跡,實(shí)驗中將2/3的數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)集,其余數(shù)據(jù)為用以驗證與調(diào) 整聚類模型準(zhǔn)確率及其他性能的驗證數(shù)據(jù)集。

3.2 結(jié)果分析

3.2.1 聚類結(jié)果

圖4所示為本文基于多特征凝聚層次聚類結(jié)果與文獻(xiàn)[19]中基于因子分析的駕駛行駛聚類結(jié)果對比。其中用黑色實(shí)線框標(biāo)出兩處聚類細(xì)節(jié),可以看出圖4(a)中將一段道路上的GPS數(shù)據(jù)點(diǎn)聚類為同一類簇,而在圖4(b)中更細(xì)致地區(qū)分了軌跡點(diǎn)之間的差異。

圖4 文獻(xiàn)[19]方法與本文方法聚類結(jié)果對比

對于本文采用的聚類算法,針對超速駕駛行為從聚類準(zhǔn)確率與算法執(zhí)行時間2個性能上,與文獻(xiàn)[19]中的基于多特征的兩層聚類算法進(jìn)行比較,比較結(jié)果如圖5所示,其中聚類準(zhǔn)確率采取將聚類結(jié)果與手動標(biāo)記的245個數(shù)據(jù)點(diǎn)進(jìn)行比較的方式計算。

聚類準(zhǔn)確率評估結(jié)果如圖5(a)所示,從中可以看出,聚類準(zhǔn)確率隨著軌跡數(shù)據(jù)規(guī)模的增大而增大。當(dāng)軌跡數(shù)目達(dá)到800左右,軌跡點(diǎn)數(shù)為9 948時,2種方法的聚類準(zhǔn)確率趨于穩(wěn)定,本文方法的聚類準(zhǔn)確率在95%左右,文獻(xiàn)[19]的聚類準(zhǔn)確率在88%左右。由實(shí)驗結(jié)果可知,本文的聚類算法與文獻(xiàn)[19]采用的方法比較在聚類準(zhǔn)確率上有一定的提高。由于本文方法同時采用了數(shù)據(jù)點(diǎn)的時間全局特征、兩點(diǎn)間行駛距離、速度大小、速度行業(yè)設(shè)定等級、是否違章超速及其統(tǒng)計量等與速度相關(guān)的完備因素來衡量駕駛行為間的相似度,更好地提取數(shù)據(jù)點(diǎn)聚類的特征工程。

執(zhí)行效率評估結(jié)果如圖5(b)所示,從中可以看出,隨著軌跡數(shù)據(jù)量的增大,本文方法在聚類效率上的優(yōu)勢就越加明顯。當(dāng)軌跡數(shù)目達(dá)到1 000條時,本文方法的聚類時間減少了24%左右。這是因為本文方法在數(shù)據(jù)點(diǎn)聚類前將數(shù)據(jù)劃分到多個網(wǎng)格區(qū)域,后期同步執(zhí)行聚類算法,且在計算數(shù)據(jù)點(diǎn)間的相似度時引入譜聚類中的拉普拉斯映射思想,降低聚類數(shù)據(jù)集的維數(shù)并能自動確定聚類數(shù)目。因此,本文方法相比文獻(xiàn)[19]方法在聚類準(zhǔn)確率和效率上有一定程度的提高。

圖5 聚類算法性能對比

3.2.2 異常駕駛行為檢測結(jié)果

圖6所示為驗證數(shù)據(jù)集中針對超速駕駛行為異常檢測結(jié)果。其中,灰色點(diǎn)是表現(xiàn)為正常的數(shù)據(jù)點(diǎn),“×”為檢測出的異常點(diǎn)。實(shí)驗數(shù)據(jù)采用訓(xùn)練數(shù)據(jù)集與驗證數(shù)據(jù)集生成的聚類簇作為模型庫。對于測試數(shù)據(jù)集,由于不能采集到及時的商用車GPS數(shù)據(jù),因此采取在測試車輛上安裝GPS定位器的方法采集人為的異常駕駛數(shù)據(jù),包括超速、急加/減速、突然左/右變道等幾項駕駛異常行為。實(shí)驗中讓3位志愿者司機(jī)共做了20次超速駕駛、40次變速駕駛(包括加速/減速)、20次變道,其余100次為正常行駛實(shí)驗,共2 000個數(shù)據(jù)點(diǎn)。

圖6 驗證數(shù)據(jù)集檢測結(jié)果示例

將本文方法與文獻(xiàn)[20]基于各類異常指標(biāo)的駕駛行為異常檢測方法進(jìn)行對比,實(shí)驗結(jié)果如表1所示,其中,M1指檢測為異常駕行為的數(shù)目,M2為正確檢測為異常駕駛行為的數(shù)目。從中可以看出,本文方法的平均檢測準(zhǔn)確率高于文獻(xiàn)[20]中提出的基于各類異常指標(biāo)的異常檢測方法。因為文獻(xiàn)[20]采用的異常檢測方法是根據(jù)已有的一些經(jīng)驗指標(biāo)來設(shè)定各類異常的閾值,對數(shù)據(jù)的測試有很強(qiáng)的人為主觀因素。而本文方法的檢測結(jié)果是將待測數(shù)據(jù)集與基于訓(xùn)練數(shù)據(jù)集和驗證數(shù)據(jù)集獲得的模型庫進(jìn)行比較獲得,在對數(shù)據(jù)集進(jìn)行訓(xùn)練時是基于提取的各類全局與局部運(yùn)動特征以及相關(guān)統(tǒng)計量等因素,將 “多且相似”的正常的駕駛行為都將歸為一個簇,而訓(xùn)練中的異常數(shù)據(jù)也將歸為異常簇,因此,測試時不受人為客觀因素的影響,檢測準(zhǔn)確率更高。

表1 檢測準(zhǔn)確性對比

4 結(jié)束語

基于GPS數(shù)據(jù)挖掘,本文針對商用車司機(jī)5種典型的異常駕駛行為進(jìn)行檢測。在最大化提取每個位置點(diǎn)駕駛行為特征的基礎(chǔ)上,采用基于聚類的異常檢測方法使準(zhǔn)確率和效率得到提升,且異常檢測準(zhǔn)確率相比傳統(tǒng)設(shè)定閾值的方法更高。但本文對算法的檢測性能分析是建立在實(shí)驗采集的異常數(shù)據(jù)集上,不包含商用車司機(jī)駕駛車輛時的一些實(shí)際信息,使檢測準(zhǔn)確率受到影響,下一步將對此進(jìn)行改進(jìn)。