分享目標(biāo)跟蹤入門篇之相關(guān)濾波

前 言

目標(biāo)跟蹤是計(jì)算機(jī)視覺(jué)領(lǐng)域的一個(gè)重要問(wèn)題，目前廣泛應(yīng)用在體育賽事轉(zhuǎn)播、安防監(jiān)控和無(wú)人機(jī)、無(wú)人車、機(jī)器人等領(lǐng)域。簡(jiǎn)單來(lái)說(shuō)，目標(biāo)跟蹤就是在連續(xù)的視頻序列中，建立所要跟蹤物體的位置關(guān)系，得到物體完整的運(yùn)動(dòng)軌跡。

給定圖像第一幀的目標(biāo)坐標(biāo)位置，計(jì)算在下一幀圖像中目標(biāo)的確切位置。在運(yùn)動(dòng)的過(guò)程中，目標(biāo)可能會(huì)呈現(xiàn)一些圖像上的變化，比如姿態(tài)或形狀的變化、尺度的變化、背景遮擋或光線亮度的變化等。目標(biāo)跟蹤算法的研究也圍繞著解決這些變化和具體的應(yīng)用展開(kāi)。目前目標(biāo)跟蹤的難點(diǎn)主要包含：

形態(tài)變化 - 姿態(tài)變化是目標(biāo)跟蹤中常見(jiàn)的干擾問(wèn)題。運(yùn)動(dòng)目標(biāo)發(fā)生姿態(tài)變化時(shí)， 會(huì)導(dǎo)致它的特征以及外觀模型發(fā)生改變， 容易導(dǎo)致跟蹤失敗。例如：體育比賽中的運(yùn)動(dòng)員、馬路上的行人。

尺度變化 - 尺度的自適應(yīng)也是目標(biāo)跟蹤中的關(guān)鍵問(wèn)題。當(dāng)目標(biāo)尺度縮小時(shí)， 由于跟蹤框不能自適應(yīng)跟蹤， 會(huì)將很多背景信息包含在內(nèi)， 導(dǎo)致目標(biāo)模型的更新錯(cuò)誤；當(dāng)目標(biāo)尺度增大時(shí)， 由于跟蹤框不能將目標(biāo)完全包括在內(nèi)， 跟蹤框內(nèi)目標(biāo)信息不全， 也會(huì)導(dǎo)致目標(biāo)模型的更新錯(cuò)誤。因此， 實(shí)現(xiàn)尺度自適應(yīng)跟蹤是十分必要的。

遮擋與消失 - 目標(biāo)在運(yùn)動(dòng)過(guò)程中可能出現(xiàn)被遮擋或者短暫的消失情況。當(dāng)這種情況發(fā)生時(shí)， 跟蹤框容易將遮擋物以及背景信息包含在跟蹤框內(nèi)， 會(huì)導(dǎo)致后續(xù)幀中的跟蹤目標(biāo)漂移到遮擋物上面。若目標(biāo)被完全遮擋時(shí)， 由于找不到目標(biāo)的對(duì)應(yīng)模型， 會(huì)導(dǎo)致跟蹤失敗。

圖像模糊 - 光照強(qiáng)度變化， 目標(biāo)快速運(yùn)動(dòng)， 低分辨率等情況會(huì)導(dǎo)致圖像模型， 尤其是在運(yùn)動(dòng)目標(biāo)與背景相似的情況下更為明顯。因此， 選擇有效的特征對(duì)目標(biāo)和背景進(jìn)行區(qū)分非常必要。

01

目標(biāo)跟蹤算法發(fā)展

跟蹤算法主要從經(jīng)典算法到基于核相關(guān)濾波算法，再到基于深度學(xué)習(xí)的跟蹤算法。早期經(jīng)典的跟蹤方法比如 Meanshift、Particle Filter 和 Kalman Filter。Meanshift 方法是一種基于概率密度分布的跟蹤方法，使目標(biāo)的搜索一直沿著概率梯度上升的方向，迭代收斂到概率密度分布的局部峰值上。

首先 Meanshift 會(huì)對(duì)目標(biāo)進(jìn)行建模，比如利用目標(biāo)的顏色分布來(lái)描述目標(biāo)，然后計(jì)算目標(biāo)在下一幀圖像上的概率分布，從而迭代得到局部最密集的區(qū)域。Meanshift 適用于目標(biāo)的色彩模型和背景差異比較大的情形，早期也用于人臉跟蹤。

由于 Meanshift 方法的快速計(jì)算，它的很多改進(jìn)方法也一直適用至今。粒子濾波（Particle Filter）方法是一種基于粒子分布統(tǒng)計(jì)的方法。以跟蹤為例，首先對(duì)跟蹤目標(biāo)進(jìn)行建模，并定義一種相似度度量確定粒子與目標(biāo)的匹配程度。

在目標(biāo)搜索的過(guò)程中，它會(huì)按照一定的分布（比如均勻分布或高斯分布）撒一些粒子，統(tǒng)計(jì)這些粒子的相似度，確定目標(biāo)可能的位置。

在這些位置上，下一幀加入更多新的粒子，確保在更大概率上跟蹤上目標(biāo)。Kalman Filter 常被用于描述目標(biāo)的運(yùn)動(dòng)模型，它不對(duì)目標(biāo)的特征建模，而是對(duì)目標(biāo)的運(yùn)動(dòng)模型進(jìn)行了建模，常用于估計(jì)目標(biāo)在下一幀的位置。

另外，經(jīng)典的跟蹤方法還有基于特征點(diǎn)的光流跟蹤，在目標(biāo)上提取一些特征點(diǎn)，然后在下一幀計(jì)算這些特征點(diǎn)的光流匹配點(diǎn)，統(tǒng)計(jì)得到目標(biāo)的位置。在跟蹤的過(guò)程中，需要不斷補(bǔ)充新的特征點(diǎn)，刪除置信度不佳的特征點(diǎn)，以此來(lái)適應(yīng)目標(biāo)在運(yùn)動(dòng)中的形狀變化。本質(zhì)上可以認(rèn)為光流跟蹤屬于用特征點(diǎn)的集合來(lái)表征目標(biāo)模型的方法。

基于核相關(guān)濾波的跟蹤算法如MOSSE、CSK、KCF、BACF、SAMF將通信領(lǐng)域的相關(guān)濾波（衡量?jī)蓚€(gè)信號(hào)的相似程度）引入到了目標(biāo)跟蹤中，相關(guān)濾波的跟蹤算法始于 2012 年 P.Martins 提出的 CSK 方法，作者提出了一種基于循環(huán)矩陣的核跟蹤方法，并且從數(shù)學(xué)上完美解決了密集采樣（Dense Sampling）的問(wèn)題，利用傅立葉變換快速實(shí)現(xiàn)了檢測(cè)的過(guò)程。

在訓(xùn)練分類器時(shí)，一般認(rèn)為離目標(biāo)位置較近的是正樣本，而離目標(biāo)較遠(yuǎn)的認(rèn)為是負(fù)樣本。利用快速傅立葉變換，CSK 方法的跟蹤幀率能達(dá)到 100～400fps，奠定了相關(guān)濾波系列方法在實(shí)時(shí)性應(yīng)用中的基石。

利用深度學(xué)習(xí)訓(xùn)練網(wǎng)絡(luò)模型，得到的卷積特征輸出表達(dá)能力更強(qiáng)，在目標(biāo)跟蹤上，初期的應(yīng)用方式是把網(wǎng)絡(luò)學(xué)習(xí)到的特征，直接應(yīng)用到相關(guān)濾波或 Struck 的跟蹤框架里面，從而得到更好的跟蹤結(jié)果，網(wǎng)絡(luò)不同層的卷積輸出都可以作為跟蹤的特征。

總的來(lái)說(shuō)相比于光流法、Kalman、Meanshift等傳統(tǒng)算法，相關(guān)濾波類算法跟蹤速度更快，深度學(xué)習(xí)類方法精度高。

具有多特征融合以及深度特征的追蹤器在跟蹤精度方面的效果更好。

使用強(qiáng)大的分類器是實(shí)現(xiàn)良好跟蹤的基礎(chǔ)。

尺度的自適應(yīng)以及模型的更新機(jī)制也影響著跟蹤的精度。

02

相關(guān)濾波器思想

相關(guān)濾波跟蹤的基本思想就是，設(shè)計(jì)一個(gè)濾波模板，利用該模板與目標(biāo)候選區(qū)域做相關(guān)運(yùn)算，最大輸出響應(yīng)的位置即為當(dāng)前幀的目標(biāo)位置。

其中 y 表示響應(yīng)輸出， x 表示輸入圖像， w 表示濾波模板。利用相關(guān)定理，將相關(guān)轉(zhuǎn)換為計(jì)算量更小的點(diǎn)積。

分別是y，x，w 的傅里葉變換。相關(guān)濾波的任務(wù)，就是尋找最優(yōu)的濾波模板w。目前的難點(diǎn)：一般的相關(guān)濾波都是固定學(xué)習(xí)率的線性加權(quán)更新模型，不需要顯式保存訓(xùn)練樣本，每幀樣本訓(xùn)練的模型與已有目標(biāo)模型，以固定權(quán)值加權(quán)來(lái)更新目標(biāo)模型，這樣以往的樣本信息都會(huì)逐漸失效，而最近幾幀的樣本信息占模型的比重很大。

如果出現(xiàn)目標(biāo)定位不準(zhǔn)確、遮擋、背景擾動(dòng)等情況，固定學(xué)習(xí)率方式會(huì)平等對(duì)待這些“有問(wèn)題”的樣本，目標(biāo)模型就會(huì)被污染導(dǎo)致跟蹤失敗。另外相關(guān)濾波模板類特征（HOG）對(duì)快速變形和快速運(yùn)動(dòng)效果不好，但對(duì)運(yùn)動(dòng)模糊光照變化等情況比較好。

03

相關(guān)濾波器發(fā)展

MOSSE

相關(guān)濾波跟蹤的開(kāi)篇之作，利用目標(biāo)的多個(gè)樣本作為訓(xùn)練樣本，以生成更優(yōu)的濾波器。MOSSE 以最小化平方和誤差為目標(biāo)函數(shù)，用m個(gè)樣本求最小二乘解。

CSK

CSK針對(duì)MOSSE算法中采用稀疏采樣造成樣本冗余的問(wèn)題，擴(kuò)展了嶺回歸、基于循環(huán)移位的近似密集采樣方法、以及核方法。MOSSE與CSK處理的都是單通道灰度圖像，引入了循環(huán)移位和快速傅里葉變換，極大地提高了算法的計(jì)算效率。

但是離散傅里葉變換也帶來(lái)了一個(gè)副作用：邊界效應(yīng)。 針對(duì)邊界效應(yīng)，有2個(gè)典型處理方法：在圖像上疊加余弦窗調(diào)制；增加搜索區(qū)域的面積。加余弦窗的方法，使搜索區(qū)域邊界的像素值接近0，消除邊界的不連續(xù)性。余弦窗的引入也帶來(lái)了缺陷： 減小了有效搜索區(qū)域。

例如在檢測(cè)階段，如果目標(biāo)不在搜索區(qū)域中心，部分目標(biāo)像素會(huì)被過(guò)濾掉。如果目標(biāo)的一部分已經(jīng)移出了這個(gè)區(qū)域，很可能就過(guò)濾掉僅存的目標(biāo)像素。其作用表現(xiàn)為算法難以跟蹤快速運(yùn)動(dòng)的目標(biāo)。擴(kuò)大搜索區(qū)域能緩解邊界效應(yīng)，并提高跟蹤快速移動(dòng)目標(biāo)的能力。但缺陷是會(huì)引入更多的背景信息，可能造成跟蹤漂移。

CN

CN在CSK的基礎(chǔ)上擴(kuò)展了多通道顏色。將RGB的3通道圖像投影到11個(gè)顏色通道，分別對(duì)應(yīng)英語(yǔ)中常用的語(yǔ)言顏色分類，分別是black， blue， brown， grey， green， orange， pink， purple， red， white， yellow，并歸一化得到10通道顏色特征。也可以利用PCA方法，將CN降維到2維。

DCF KCF

從DCF到KCF多了Gaussian-kernel，performance上升0.21%，fps下降46.46%，kernel-trick雖然有用但影響較小，如果注重速度可以摒棄，如果追求極限性能可以用。KCF可以說(shuō)是對(duì)CSK的完善。

論文中對(duì)嶺回歸、循環(huán)矩陣、核技巧、快速檢測(cè)等做了完整的數(shù)學(xué)推導(dǎo)。KCF在CSK的基礎(chǔ)上擴(kuò)展了多通道特征。KCF采用的HoG特征，核函數(shù)有三種高斯核、線性核和多項(xiàng)式核，高斯核的精確度最高，線性核略低于高斯核，但速度上遠(yuǎn)快于高斯核。

SAMF

SAMF基于KCF，特征是HoG+CN。SAMF實(shí)現(xiàn)多尺度目標(biāo)跟蹤的方法比較直接，類似檢測(cè)算法里的多尺度檢測(cè)方法。由平移濾波器在多尺度縮放的圖像塊上進(jìn)行目標(biāo)檢測(cè)，取響應(yīng)最大的那個(gè)平移位置及所在尺度。因此這種方法可以同時(shí)檢測(cè)目標(biāo)中心變化和尺度變化。

DSST fDSST

從DSST到fDSST做了特征壓縮和scale filter加速即特征降維和插值，performance上升6.13%，fps上升83.37%。 DSST將目標(biāo)跟蹤看成目標(biāo)中心平移和目標(biāo)尺度變化兩個(gè)獨(dú)立問(wèn)題。首先用HoG特征的DCF訓(xùn)練平移相關(guān)濾波，負(fù)責(zé)檢測(cè)目標(biāo)中心平移。然后用HoG特征的MOSSE（這里與DCF的區(qū)別是不加padding）訓(xùn)練另一個(gè)尺度相關(guān)濾波，負(fù)責(zé)檢測(cè)目標(biāo)尺度變化。

2017年發(fā)表的文章又提出了加速版本fDSST。 尺度濾波器僅需要檢測(cè)出最佳匹配尺度而無(wú)須關(guān)心平移情況。DSST將尺度檢測(cè)圖像塊全部縮小到同一個(gè)尺寸計(jì)算特征（CN+HoG），再將特征表示成一維（沒(méi)有循環(huán)移位），尺度檢測(cè)的響應(yīng)圖也是一維的高斯函數(shù)。 DSST本來(lái)就是對(duì)尺度自適應(yīng)問(wèn)題的快速解決方案（支持33個(gè)尺度還比SAMF快很多），在fDSST中MD大神又對(duì)DSST進(jìn)行加速：

平移濾波器：PCA方法將平移濾波器的HOG特征從31通道降維到18通道，這一步驟與上面的CN特征類似，直接用PCA進(jìn)行降維，作者提到由于這里用了線性核，所以不需要CN中所用的平滑子空間約束，更加簡(jiǎn)單粗暴。

由于HOG特征天然會(huì)降低響應(yīng)分辨率（cell_size=4），這里也采用簡(jiǎn)單粗暴的方法，將響應(yīng)圖的分辨率上采樣到原始圖像分辨率，也就是響應(yīng)圖插值以提高檢測(cè)精度，方法是三角插值，等價(jià)于頻譜添0，方法更加簡(jiǎn)單粗暴，但這一步會(huì)增加算法復(fù)雜度，而且方法太簡(jiǎn)單也必然效果較差。

尺度濾波器：QR方法將尺度濾波器的HOG特征（二特征，沒(méi)有循環(huán)移位）~1000*17降維到17*17，由于自相關(guān)矩陣維度較大影響速度，為了效率這里沒(méi)有用PCA而是QR分解。多尺度數(shù)量是17（DSST中的一半），響應(yīng)圖是1*17，這里也通過(guò)插值方法將尺度數(shù)量從17插值到33以獲得更精確的尺度定位。

SRDCF

SRDCF與CFLB的思路都是擴(kuò)大搜索區(qū)域，同時(shí)約束濾波模板的有效作用域解決邊界效應(yīng)。給濾波模板增加一個(gè)約束，對(duì)接近邊界的區(qū)域懲罰更大，或者說(shuō)讓邊界附近濾波模板系數(shù)接近0，速度比較慢。

CFLB/BACF

使搜索區(qū)域內(nèi)，目標(biāo)區(qū)域以外的像素為0，CFLB僅使用單通道灰度特征，最新BACF將特征擴(kuò)展為多通道HOG特征。CFLB和BACF采用Alternating Direction Method of Multipliers（ADMM）快速求解。

DAT

不是一種相關(guān)濾波方法，而是一種基于顏色統(tǒng)計(jì)特征方法。DAT統(tǒng)計(jì)前景目標(biāo)和背景區(qū)域的顏色直方圖，這就是前景和背景的顏色概率模型，檢測(cè)階段，利用貝葉斯方法判別每個(gè)像素屬于前景的概率，得到像素級(jí)顏色概率圖

STAPLE STAPLE+CA

從Staple到STAPLE+CA加入Context-Aware約束項(xiàng)，performance上升3.28%，fps下降43.18%，說(shuō)明約束項(xiàng)有效，但犧牲了大量fps。STAPLE結(jié)合了模板特征方法DSST和顏色統(tǒng)計(jì)特征方法DAT。 相關(guān)濾波模板類特征（HOG）對(duì)快速變形和快速運(yùn)動(dòng)效果不好，但對(duì)運(yùn)動(dòng)模糊光照變化等情況比較好；而顏色統(tǒng)計(jì)特征（DAT）對(duì)變形不敏感，而且不屬于相關(guān)濾波框架沒(méi)有邊界效應(yīng)，但對(duì)光照變化和背景相似顏色不好。因此，這兩類方法可以互補(bǔ)。

C-COT

圖像特征的表達(dá)能力在目標(biāo)跟蹤中起著至關(guān)重要的作用。以HoG+CN為代表的圖像特征，性能優(yōu)秀而且速度優(yōu)勢(shì)非常突出，但也成為性能進(jìn)一步提升的瓶頸。以卷積神經(jīng)網(wǎng)絡(luò)（CNN）為代表的深度特征，具有更強(qiáng)大特征表達(dá)能力、泛化能力和遷移能力。將深度特征引入相關(guān)濾波也就水到渠成。

LMCF

LMCF提出了兩個(gè)方法，多峰目標(biāo)檢測(cè)和高置信度更新。多峰目標(biāo)檢測(cè)對(duì)平移檢測(cè)的響應(yīng)圖做多峰檢測(cè)，如果其他峰峰值與主峰峰值的比例大于某個(gè)閾值，說(shuō)明響應(yīng)圖是多峰模式，以這些多峰為中心重新檢測(cè)，并取這些響應(yīng)圖的最大值作為最終目標(biāo)位置。

高置信度更新：只有在跟蹤置信度比較高的時(shí)候才更新跟蹤模型，避免目標(biāo)模型被污染。一個(gè)置信度指標(biāo)是最大響應(yīng)。另一個(gè)置信度指標(biāo)是平均峰值相關(guān)能量（average peak-to correlation energy， APCE），反應(yīng)響應(yīng)圖的波動(dòng)程度和檢測(cè)目標(biāo)的置信水平。

CSR-DCF

提出了空域可靠性和通道可靠性方法?？沼蚩煽啃岳脠D像分割方法，通過(guò)前背景顏色直方圖概率和中心先驗(yàn)計(jì)算空域二值約束掩膜。這里的二值掩膜就類似于CFLB中的掩膜矩陣P。CSR-DCF利用圖像分割方法更準(zhǔn)確地選擇有效的跟蹤目標(biāo)區(qū)域。通道可靠性用于區(qū)分檢測(cè)時(shí)每個(gè)通道的權(quán)重。

ECO ECO-HC

ECO是C-COT的加速版，從模型大小、樣本集大小和更新策略三個(gè)方便加速，速度比C-COT提升了20倍，加量還減價(jià)，在VOT2016數(shù)據(jù)庫(kù)上EAO提升了13.3%，當(dāng)然最厲害的還是hand-crafted features版本的ECO-HC有60FPS，接下來(lái)分別看看這三步。

第一減少模型參數(shù)，既然CN特征和HOG特征都能降維，那卷積特征是不是也可以試試？這就是ECO中的加速第一步，也是最關(guān)鍵的一步，F(xiàn)actorized Convolution Operator分解卷積操作，效果類似PCA，但Conv. Feat.與前面的CN和HOG又不一樣：

CNN特征維度過(guò)于龐大，在C-COT中是96+512=608通道，需要降很多很多維才能保證速度，而無(wú)監(jiān)督降維如果太多會(huì)直接影響效果（對(duì)比通用方法 - 取特征值的95%以上的維度，保留信息量）；

雖然CNN特征遷移能力比較強(qiáng)，但這并不是針對(duì)跟蹤問(wèn)題專門訓(xùn)練的特征，對(duì)跟蹤問(wèn)題有用的信息隱藏在大量CNN激活值中，如果簡(jiǎn)單的無(wú)監(jiān)督降維，可能會(huì)過(guò)濾掉那些雖然不顯著，但對(duì)跟蹤問(wèn)題有效的特征信息。當(dāng)然HOG和CN特征也有同樣的問(wèn)題。

通過(guò)使用PCA， 有監(jiān)督降維：

P就是那個(gè)降維矩陣，放在目標(biāo)函數(shù)中優(yōu)化得到，具體求解比較復(fù)雜看論文吧，用PCA作為P的初始值去迭代優(yōu)化，采用Gauss-Newton和Conjugate Gradient方法。但每幀都迭代優(yōu)化降維矩陣速度反而會(huì)更慢，大神告訴我們僅在第一幀優(yōu)化這個(gè)降維矩陣就可以了，第一幀優(yōu)化完成后這個(gè)降維矩陣就是固定的，后續(xù)幀都直接用。

Factorized Convolution Operator減少了80%的卷積特征還能略微提升性能，HC版本從31+11降維到10+3速度提升非常明顯。至于為什么降維了還能提升效果，論文中說(shuō)參數(shù)太多容易過(guò)擬合，也可能是判別力較低或者無(wú)用的通道響應(yīng)圖會(huì)成為噪聲。

第二減少樣本數(shù)量，這個(gè)是針對(duì)Adaptive decontamination of the training set的加速，C-COT中要保存400個(gè)樣本，但視頻相鄰幀之間的相似性非常高，存在大量相似的冗余樣本，而且每次更新都要用樣本集中所有樣本做優(yōu)化，速度非常慢。

ECO中改為緊湊的生成樣本空間模型compact generative mode，采用Gaussian Mixture Model （GMM）合并相似樣本，建立更具代表性和多樣性的樣本集，需要保存和優(yōu)化的樣本集數(shù)量降到C-COT的1/8。用特征距離衡量?jī)蓚€(gè)樣本的相似程度，樣本合并方法是兩個(gè)樣本特征的權(quán)值相加，樣本特征按照權(quán)值加權(quán)合并。

第三改變更新策略，以前CF方法都是每幀更新，這種過(guò)更新不僅慢，而且會(huì)導(dǎo)致模型對(duì)最近幾幀嚴(yán)重過(guò)擬合，對(duì)遮擋、變形和平面外旋轉(zhuǎn)等突然變化過(guò)度敏感，但對(duì)大多數(shù)方法都是無(wú)可奈何的，因?yàn)槿鏚CF等方法不保存樣本，這一幀不更新就再也沒(méi)機(jī)會(huì)了。

但ECO保存了所有樣本的代表性樣本集，所以完全沒(méi)有必要每幀都更新，這里采用了sparser updating scheme（稀疏更新策略），即每隔5幀更新一次模型參數(shù)，這不但提高了算法速度，而且提高了對(duì)突變，遮擋等情況的穩(wěn)定性，三步優(yōu)化中稀疏更新對(duì)效果提升最大。

由于ECO的樣本集是每幀都更新的，稀疏更新并不會(huì)錯(cuò)過(guò)間隔期的樣本變化信息，但這種方法可能不適合沒(méi)有樣本集的方法，如KCF，因?yàn)闆](méi)有保存樣本集。

編輯：jq