沈澤君,丁飛飛,楊文元,b

(閩南師范大學 a.福建省粒計算及其應(yīng)用重點實驗室; b.計算機學院,福建 漳州 363000)

0 概述

視頻跟蹤技術(shù)是計算機視覺領(lǐng)域的一個重要研究方向,其目的是在連續(xù)的圖像序列中建立一個可更新跟蹤目標的動態(tài)表觀模型[1]。跟蹤技術(shù)廣泛應(yīng)用于自動駕駛、醫(yī)學圖像跟蹤分析、安防監(jiān)控、軍事偵察等領(lǐng)域[2]。雖然在過去幾十年中視頻跟蹤技術(shù)得到了快速的發(fā)展[3-5],但其還面臨一定挑戰(zhàn),目標外觀變化或目標模糊等不確定因素均給跟蹤過程增加了難度。

近年來,基于生成學習類算法和判別學習類算法的跟蹤技術(shù)相繼被提出,如文獻[6]提出的自適應(yīng)局部結(jié)構(gòu)稀疏表觀模型和文獻[7]提出的利用稀疏協(xié)同模型引入判別信息的稀疏表觀模型。雖然上述模型提高了目標跟蹤的精準度,但未解決如何降低復雜度、提高魯棒性的問題。判別式相關(guān)濾波(DCF)跟蹤算法在平衡跟蹤精度和速度之間取得了較好效果。DCF跟蹤算法通過傅里葉變換將卷積計算轉(zhuǎn)換為點乘計算,使其具有較高的計算效率,有效地降低了算法的復雜度。但是,該算法對樣本的數(shù)量有較高要求。為解決該問題,文獻[8]提出一種最小均方誤差輸出總和濾波器(MOSSE)。在MOSSE的基礎(chǔ)上,研究人員提出了大量的改進算法[1],如文獻[9]提出的基于循環(huán)矩陣結(jié)構(gòu)核相關(guān)濾波跟蹤算法(CSK)。為了進一步提高對尺度變化的處理效果,文獻[10]提出具有特征融合的尺度自適應(yīng)核相關(guān)濾波跟蹤器(SAMF)。近幾年,基于結(jié)構(gòu)相關(guān)濾波器[11]、基于長時域跟蹤相關(guān)濾波器[12-14]、基于端到端相關(guān)濾波器[15-16]、基于核相關(guān)濾波器[17-18]和基于過濾式相關(guān)濾波器[19-21]相繼被提出。其中,將DCF和深度卷積神經(jīng)網(wǎng)絡(luò)(CNN)特征相結(jié)合,可以得到較好的效果[22]。但是,在已有的相關(guān)濾波器跟蹤算法中,通過分析各特征之間的粒度關(guān)系來提高算法精準度和時間效率的研究較少[23-24],并且多數(shù)跟蹤算法只考慮單一特征的提取或者只對不同特征進行簡單結(jié)合。

粒計算思維通過劃分不同的粒度來求解問題[25]。粒子是粒計算模型的最小單位,一個較小的特征可以構(gòu)成一個粒子,一個覆蓋面較廣的特征可以構(gòu)成一個較大的粒子。不同粒子可能位于不同的粒層,不同的?；瘶藴蕦ν粋€跟蹤算法生成的粒度也不同[26]。粒計算主要考慮3個問題,一是如何選取合適的?；瘶藴蕘泶_定粒層,二是如何在不同粒層上求解問題,三是如何在各個粒層間相互置換和融合[25]。

本文提出一種基于粒計算思維的多粒度相關(guān)濾波視頻跟蹤算法(MGCF)。MGCF算法對輸入的視頻圖像中每一幀進行基礎(chǔ)特征提取,組成一個基礎(chǔ)特征池,通過選取特征池中不同的基礎(chǔ)特征構(gòu)造出不同粒度的特征序列,根據(jù)特征序列構(gòu)成不同粒度的相關(guān)濾波器。對該幀進行各粒度跟蹤,依據(jù)對比評估和自我評估來計算各粒度下輸出結(jié)果的穩(wěn)健性評估得分,根據(jù)評估得分的高低選擇最優(yōu)的粒度,該粒度的跟蹤結(jié)果即為該幀的輸出,系統(tǒng)累計各幀的輸出并作為最終的跟蹤結(jié)果。

1 相關(guān)工作

1.1 濾波視頻跟蹤算法

DCF跟蹤算法的計算效率較高,其利用傅里葉變換將傳統(tǒng)的卷積運算轉(zhuǎn)化成點乘運算。設(shè)輸入的特征為x,g為相關(guān)濾波器,根據(jù)卷積理論可得:

(1)

(2)

訓練后的卷積核為:

(3)

根據(jù)上述方法可以在對圖像進行循環(huán)卷積操作時大幅降低算法的復雜度。對一個n×n的圖像分別進行傳統(tǒng)的卷積操作和利用傅里葉變換進行點乘運算,其算法復雜度分別為O(n4)和O(n2lgn)。

利用卷積的性質(zhì),可以通過第i個輸入快速找到其對應(yīng)的標簽yi[27]。將相關(guān)濾波器當成一個分類器,假設(shè)存在關(guān)系f(xi)=yi,則可以將訓練過程等價于優(yōu)化以下的目標函數(shù):

(4)

其中,δ為模型參數(shù),λ為防止過擬合的正則項系數(shù),G為損失函數(shù),具體如下:

G(f(δ,xi),yi)=(yi-f(δ,xi))2

(5)

函數(shù)f(xi)的表現(xiàn)形式多樣,在式(5)中可表示為線性操作,如f(xi)=〈δ,xi〉+b。其中,〈〉為點積,b為常數(shù)偏值。在此線性定義下,可以利用二次誤差的正則最小二乘法進行優(yōu)化,得到如下閉合解:

δ=(XTX+λI)-1XTy

(6)

其中,X矩陣的行向量表示訓練樣本,y向量為對應(yīng)的標簽信息,I為單位矩陣。

可以通過引入核函數(shù)來提高算法性能。例如,利用β(x)將輸入值x映射到一個非線性的特征空間并匯總,使得參數(shù)δ可以用一個線性組合類來表示:

(7)

根據(jù)文獻[28]的研究結(jié)果,最終的函數(shù)f(xi)如式(8)所示:

(8)

其中,κ(xi,xj)=〈β(xi),β(xj)〉,核矩陣k中對應(yīng)的元素kij=κ(xi,xj),則求解的優(yōu)化問題可以轉(zhuǎn)變?yōu)?

α=(K+λI)-1y

(9)

跟蹤的過程就是利用一個已訓練好的參數(shù)和一個基于樣本所構(gòu)成的模型來檢測待測區(qū)域。當待測區(qū)域出現(xiàn)時,最終的響應(yīng)可以由式(10)得到:

(10)

1.2 粒計算的粒度理論

(11)

B的粒度反映其分辨能力[29]。GD(B)表示在U中隨機挑選2個對象的B-不可分辨的可能性的大小,其值越大表示分辨能力越弱。

2 多粒度相關(guān)濾波器

如圖1所示,MGCF算法對輸入的圖像提取基礎(chǔ)特征構(gòu)成一個特征池,通過對特征池中的基礎(chǔ)特征進行組合,構(gòu)成不同粒度下的相關(guān)濾波器,然后根據(jù)其穩(wěn)健性評估得分的高低選擇最優(yōu)結(jié)果,粒度的跟蹤輸出就作為該幀的結(jié)果。

圖1 MGCF算法的系統(tǒng)流程

2.1 基礎(chǔ)特征池與粒度池

MGCF算法有一個粒度池用以收納各粒度下的相關(guān)濾波器,其相關(guān)濾波器是由基礎(chǔ)特征池中的4組基礎(chǔ)特征組合形成。第一組基礎(chǔ)特征為文獻[30]中使用的方向梯度直方圖(HOG)特征,記為C1,余下3組使用文獻[31]提出的HOG特征版本,第二組使用其1維～18維,第三組使用其19維～27維,第四組使用其28維～31維,分別記為C2、C3、C4。綜上,基礎(chǔ)特征池為{C1,C2,C3,C4},并組合形成14個相關(guān)濾波器,記為Gi。各粒度下的相關(guān)濾波器對應(yīng)的特征組合如表1所示。

表1 各粒度相關(guān)濾波器對應(yīng)的特征組合

Table 1 Characteristic combination corresponding to each granularity of correlation filter

相關(guān)濾波器特征組合相關(guān)濾波器特征組合G1{C1}G8{C2,C3}G2{C2}G9{C2,C4}G3{C3}G10{C3,C4}G4{C4}G11{C1,C2,C3}G5{C1,C2}G12{C1,C2,C4}G6{C1,C3}G13{C2,C3,C4}G7{C1,C4}G14{C1,C2,C3,C4}

雖然與G11～G14相比,G1～G4的特征較少,但是它們?yōu)楦櫧Y(jié)果提供了多樣性,這對于整個跟蹤系統(tǒng)的判斷是必不可少的。

2.2 多粒度相關(guān)跟蹤器

本文MGCF算法是可多線程并行的跟蹤算法,并且可以基于各粒度下的跟蹤判斷綜合評價出一個最優(yōu)解給每一幀。如圖2所示,節(jié)點為各粒度下的預(yù)測結(jié)果,通過對比評估和自我評估可得出每個粒度的穩(wěn)健性評估得分。

圖2 評估過程示意圖

2.2.1 對比評估

(12)

(13)

(14)

(15)

(16)

2.2.2 自我評估

(17)

2.2.3 各粒度的穩(wěn)健性評估與篩選

(18)

其中,μ為參數(shù)。在各幀下,系統(tǒng)根據(jù)各粒度的穩(wěn)健性評估得分來篩選出最優(yōu)粒度,以此作為當前的跟蹤結(jié)果。

2.3 MGCF算法模型

如圖3所示(彩色效果見《計算機工程》官網(wǎng)電子版),多粒度相關(guān)濾波跟蹤算法共有14個粒度濾波器對圖像進行跟蹤,在每一幀中,選擇穩(wěn)健性評估得分最高的粒度作為該幀的結(jié)果,如MGCF算法在第32幀中選擇粒度13、在第142幀中選擇粒度7、在第206幀中選擇粒度6。MGCF算法模型如圖4所示。

圖3 MGCF算法的粒度選擇過程

圖4 MGCF算法模型

MGCF算法步驟具體如下:

算法1MGCF算法

輸入視頻序列X={x1,x2,…,xn},第一幀的初始位置

1.對第i幀的圖像數(shù)據(jù)進行基礎(chǔ)特征提取

2.由基礎(chǔ)特征組合各粒度的相關(guān)濾波器

4.進入下一幀,返回步驟1,直到n

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

在OTB-2013和OTB-2015 2個公開數(shù)據(jù)集上,將本文MGCF算法與SRDCF、C-COT、SAMF、Staple、KCF和DCFNet 6個視頻跟蹤算法進行比較與分析。

3.1 實驗設(shè)置

實驗使用CPU為Intel(R)Core(TM)i7-6500U、GPU為GeForce 920MX、內(nèi)存為8 GB的電腦,MATLAB使用2018b版本。式(18)中的μ設(shè)置為0.1,實驗過程中的其他參數(shù)如表2所示。

表2 實驗參數(shù)設(shè)置

3.2 結(jié)果分析

在實驗中,使用魯棒性作為定量分析標準,將文本算法與SRDCF、C-COT、SAMF、Staple、KCF和DCFNet 6個算法進行像素誤差和重疊率比較。像素誤差是預(yù)測位置與實際位置的像素位移差,重疊率是重疊區(qū)域除以2個矩形總面積的比值。魯棒性評估分為空間魯棒性(SRE)和時間魯棒性(TRE)2種,MGCF算法在空間魯棒性中能獲得較高的性能,但是在時間魯棒性上表現(xiàn)不足。

3.2.1 像素誤差性能分析

如圖5(a)～圖5(e)所示(中括號內(nèi)數(shù)字為精準度平均值),MGCF算法在快速運動、平面內(nèi)外旋轉(zhuǎn)和尺度變化這4種情況下能取得最佳效果,精準度分別為91.7%、93.1%、91.3%和90.9%,比性能第2的算法分別高出1.3%、2.3%、2.7%和1.7%。但是,在遮擋情況下MGCF算法表現(xiàn)不佳,精準度只有93.3%,比SAMF算法低1.4%。

圖5 7種算法的像素誤差性能分析

3.2.2 重疊率性能分析

如圖6(a)～圖6(e)所示(中括號內(nèi)數(shù)字為成功率平均值),MGCF算法在快速運動、平面內(nèi)外旋轉(zhuǎn)和尺度變化4種情況下均能取得最佳性能,成功率分別為86.1%、85.4%、82.2%和83.9%,比性能第2的算法分別高出2.4%、4.3%、5.8%和1.0%。但是,在遮擋情況下MGCF算法表現(xiàn)不佳,成功率只有74.6%,比Staple算法低5.3%。

圖6 7種算法的重疊率性能分析

3.2.3 綜合比較

如圖7所示,MGCF算法在空間魯棒性下的像素誤差和重疊率分別為94.2%和87.0%,C-COT、Staple、SRDCF的精準度分別為87.8%、78.9%、74.8%,成功率分別為75.6%、70.4%、67%。MGCF相對C-COT、Staple、SRDCF的精準度分別提高6.4%、15.3%、19.4%,成功率分別提高11.4%、16.6%、20%。

圖7 7種算法的精準度和成功率對比

3.2.4 運算時效性分析

本實驗采用CPU為Intel(R)Core(TM)i7-6500U的計算機進行運算,各算法運行速度如表3所示,MGCF算法速度為7.8 FPS,該算法使用4組基礎(chǔ)特征構(gòu)造出14個粒度相關(guān)濾波器,對每一幀的圖像進行自我評估和對比評估,并據(jù)此得出最終的穩(wěn)健性評估值。由于MGCF算法根據(jù)穩(wěn)健性評估值的高低排序來篩選最后的跟蹤結(jié)果,需要分別用各個粒度特征來更新濾波模板,然后對每一幀進行跟蹤評估,其運算量較大,導致運算時效性不高。

表3 各算法運行速度比較

4 結(jié)束語

本文結(jié)合粒計算思維提出一種多粒度相關(guān)濾波跟蹤算法MGCF,該算法在多個相關(guān)濾波器中提出了粒度的概念,從粒度角度將分散的多線程跟蹤進行有效地結(jié)合,解決了傳統(tǒng)視頻跟蹤算法的特征單一性問題。實驗結(jié)果表明,MGCF算法具有較好的視頻跟蹤性能,成功率較高,但在時間魯棒性實驗中,存在多個干擾因素時其性能不佳,精準度較低。下一步將結(jié)合深度特征和遷移學習等技術(shù)解決MGCF算法在時間魯棒性實驗中精準度不高的問題。