劉 嵩,譚建軍

(湖北民族學(xué)院 信息工程學(xué)院,湖北 恩施445000)

1 引 言

人臉識(shí)別是計(jì)算機(jī)視覺中的一個(gè)熱點(diǎn),其較為關(guān)鍵的技術(shù)就是特征提取。主元分析法(Principle Component Analysis,PCA)是經(jīng)典的特征提取方法,它提取的特征是代數(shù)特征,滿足投影方差最大的原則,且是相互正交的[1]。但PCA 方法沒有利用樣本數(shù)據(jù)的高階信息,而高階信息更有利于人臉分類?；诤说姆蔷€性變換因其在支持向量機(jī)(SVM)中出色應(yīng)用而引起眾多研究人員關(guān)注,為了利用樣本的高階信息,研究人員將核函數(shù)的思想引入到人臉識(shí)別中。Scholkopf 提出了核主元分析(Kernel Principle Component Analysis,KPCA)的方法來提取特征,并應(yīng)用于人臉識(shí)別[2]。KPCA 基本思想是:首先通過非線性映射將原始樣本數(shù)據(jù)從低維空間變換到高維特征空間,然后在這個(gè)新的空間中使用PCA 方法。為避免“維數(shù)災(zāi)難”問題,引用支持向量機(jī)(Support Vector Machine,SVM)中的核函數(shù),即用滿足Mercer 條件的核函數(shù)來計(jì)算特征空間中向量的內(nèi)積[3-5],避免了顯式的求解映射函數(shù)。實(shí)驗(yàn)結(jié)果表明,KPCA不僅能夠抽取非線性特征,而且有著比主成分分析更好的識(shí)別效果。

然而,基于KPCA 的人臉識(shí)別方法需要多次計(jì)算內(nèi)積來求解核矩陣,其耗時(shí)占總識(shí)別時(shí)間的80%以上,特別是在樣本維數(shù)較高的時(shí)候更為明顯。另外考慮到在工程實(shí)踐中對(duì)成本、體積、后期升級(jí)等方面的要求,傳統(tǒng)的基于PC 機(jī)的人臉識(shí)別算法在實(shí)際應(yīng)用中受到很大限制。隨著數(shù)字信號(hào)處理(Digital Signal Processing,DSP)技術(shù)的快速發(fā)展,其運(yùn)算能力顯著提高,能實(shí)時(shí)處理海量數(shù)據(jù),同時(shí)擁有較高的性價(jià)比[6]。鑒于此,本文提出了一種利用DSP 實(shí)現(xiàn)核主元分析的方法,實(shí)驗(yàn)結(jié)果證明了該方法的有效性。

2 KPCA 算法

假設(shè)空間映射為Υ,x1, x2, …, xM為訓(xùn)練樣本,則映射可以定義如下:

假設(shè)特征空間數(shù)據(jù)滿足中心化的條件,即

則特征空間中的協(xié)方差矩陣為

現(xiàn)求Q 的特征值λ≥0 和特征向量ν滿足下式:

即有

考慮到所有的特征向量可表示為 Υ(x1),Υ(x2), …, Υ(xM)的線性張成,即

則有

定義M×M 維矩陣K:

則式(7)可以簡(jiǎn)化為

顯然滿足

求解式(10)就可以得到特征值和特征向量,對(duì)于測(cè)試樣本在特征向量空間Vk的投影為

將內(nèi)積用核函數(shù)替換則有

當(dāng)式(2)不成立時(shí),需進(jìn)行調(diào)整:

則核矩陣可修正為

常見的核函數(shù)有線性核函數(shù)、d 階多項(xiàng)式核函數(shù)和高斯徑向基函數(shù)核函數(shù)3 種,分別為

3 算法的實(shí)現(xiàn)與優(yōu)化

3.1 DSP 硬件平臺(tái)

TI 公司的TMS320C6000 系列DSP 具有較強(qiáng)的數(shù)據(jù)處理能力,特別適合做圖像處理?？紤]到后期的升級(jí)和總的開發(fā)成本,本系統(tǒng)選擇TMS320C6201芯片作為核心處理器搭建硬件平臺(tái)。TMS320C6201是一種高性能的定點(diǎn)數(shù)字信號(hào)處理器[7]。

系統(tǒng)實(shí)現(xiàn)的硬件平臺(tái)總體框圖如圖1 所示,系統(tǒng)以DSP 為核心, 包括DSP 處理模塊、電源模塊、LED 顯示模塊、存儲(chǔ)模塊和擴(kuò)展接口等部分,其中電源模塊為系統(tǒng)提供電能供應(yīng),顯示模塊顯示識(shí)別結(jié)果,存儲(chǔ)模塊包括SDRAM 模塊和FLASH 模塊,分別用來運(yùn)行程序和固化程序,接口模塊是為了后期系統(tǒng)升級(jí)預(yù)留的接口。

圖1 系統(tǒng)總體框圖Fig.1 Diagram of system hardware

3.2 算法實(shí)現(xiàn)

移植到TMS320C6201 平臺(tái)上的人臉識(shí)別算法實(shí)現(xiàn)流程如圖2 所示。系統(tǒng)上電后,首先進(jìn)行初始化,然后按預(yù)先設(shè)定的程序讀取訓(xùn)練圖像,在預(yù)處理階段運(yùn)行幾何歸一化和灰度歸一化程序。幾何歸一化是指將圖像中人臉變換到同一位置和同樣大小,得到大小、位置相當(dāng)?shù)膱D像?；叶葰w一化將圖像的灰度值和方差歸一化到一個(gè)特定的區(qū)間,從而將減少光照變化對(duì)人臉識(shí)別率的影響。訓(xùn)練圖像和測(cè)試圖像經(jīng)過預(yù)處理后再通過核主元分析方法提取人臉圖像的高維非線性特征,最后經(jīng)過分類器進(jìn)行分類識(shí)別并顯示識(shí)別結(jié)果。

圖2 算法流程圖Fig.2 Algorithm flow chart

3.3 算法優(yōu)化

優(yōu)化是指通過提高硬件資源的并行利用程度,提高代碼運(yùn)行速度,減少運(yùn)行時(shí)間[8]。DSP 的算法優(yōu)化一般經(jīng)過3 個(gè)階段,流程如圖3 所示。

圖3 算法優(yōu)化流程圖Fig.3 Algorithm op timization flow chart

第一階段是直接使用C 語言實(shí)現(xiàn)算法功能。TI公司提供了片上支持庫(Chip Support Library,CSL)來簡(jiǎn)化程序開發(fā)。CSL 是一個(gè)C 語言程序庫,給出了配置和控制片上外設(shè)的C 語言接口。雖然也可以用寫匯編程序的方法來操作片上外設(shè),但使用片上支持庫可以帶來以下好處:使外設(shè)的使用變得容易、縮短開發(fā)時(shí)間、提供了所有外設(shè)寄存器和寄存器字段的字符表示、方便將代碼移植到其他DSP。

第二階段是C 語言程序的優(yōu)化。對(duì)于耗時(shí)較長(zhǎng)的程序段,采取從軟件入手, 兼顧硬件的優(yōu)化原則,重新組織程序結(jié)構(gòu)。

第三階段是匯編語言優(yōu)化。找出不能滿足效率要求的代碼段,用線性匯編語言改寫,利用匯編優(yōu)化器進(jìn)行優(yōu)化。

在本實(shí)驗(yàn)中主要采取了如下優(yōu)化手段。

(1)存儲(chǔ)器優(yōu)化

DSP 內(nèi)部存儲(chǔ)空間不大,在優(yōu)化過程中需要合理配置片內(nèi)存儲(chǔ)器,避免反復(fù)的搬移數(shù)據(jù)增大程序運(yùn)行時(shí)間。因?yàn)閳D像數(shù)據(jù)都比較大,可以存儲(chǔ)在外部存儲(chǔ)器,而求得的特征矢量因?yàn)橐磸?fù)使用,可以存儲(chǔ)在內(nèi)部存儲(chǔ)器,同時(shí)開辟一部分片內(nèi)程序作為cache,提高程序訪問數(shù)據(jù)的速度。

(2)編譯器優(yōu)化

編譯器選項(xiàng)用來控制編譯器的行為,選擇不同的編譯器選項(xiàng),通過設(shè)置參數(shù)來優(yōu)化編譯過程, 可以提高編譯后匯編代碼的效率。

(3)EDMA 技術(shù)

高效的數(shù)據(jù)處理能力也需要高效的數(shù)據(jù)供給,DMA 數(shù)據(jù)搬移技術(shù)可在沒有CPU 參與的情況下進(jìn)入后臺(tái)操作,達(dá)到高效搬移的功能。

(4)軟件流水

在人臉識(shí)別中有大量的循環(huán)運(yùn)算,軟件流水是用來安排循環(huán)指令,使這個(gè)循環(huán)的多次迭代并行執(zhí)行的一種技術(shù),是提高代碼性能的最關(guān)鍵優(yōu)化方式。可以采用人工的干預(yù),使編譯器生成效率更高的軟件流水,這些干預(yù)包括消除冗余循環(huán)、向編譯器傳遞循環(huán)次數(shù)信息、循環(huán)展開等。

(5)線性匯編

在識(shí)別程序中,內(nèi)積運(yùn)算的C 語言程序效率較低,使用線性匯編重新改寫,可以大大提高計(jì)算速度。

(6)定點(diǎn)運(yùn)算

由于TMS320C6201 是定點(diǎn)系列的DSP,使用浮點(diǎn)數(shù)運(yùn)算勢(shì)必會(huì)加大DSP 的運(yùn)行時(shí)間,不利于DSP算法的運(yùn)行效率,所以在程序中將運(yùn)算轉(zhuǎn)為定點(diǎn)運(yùn)算,可以節(jié)約運(yùn)算時(shí)間。

4 實(shí)驗(yàn)測(cè)試

實(shí)驗(yàn)使用的人臉圖像集來自O(shè)RL 人臉數(shù)據(jù)庫。該數(shù)據(jù)庫的人臉圖像分別是在不同時(shí)期、不同光照、不同頭部角度和不同表情條件下攝制而得,由40 個(gè)人每人10 幅共400 幅構(gòu)成,是目前使用最廣泛的標(biāo)準(zhǔn)圖像。在實(shí)驗(yàn)中隨機(jī)抽取數(shù)據(jù)庫中的部分圖像作為訓(xùn)練樣本,在余下部分中選擇相同數(shù)目的圖像作為測(cè)試樣本,訓(xùn)練樣本數(shù)目與測(cè)試樣本數(shù)目相等,而且統(tǒng)一選擇高斯徑向基函數(shù)作為核函數(shù),核函數(shù)的參數(shù)取值為10,選擇的分類器是最小距離分類器。每一個(gè)實(shí)驗(yàn)都重復(fù)30 次,對(duì)實(shí)驗(yàn)結(jié)果取平均值。

實(shí)驗(yàn)中首先對(duì)比了KPCA 算法在PC 機(jī)運(yùn)行和本硬件平臺(tái)運(yùn)行的人臉識(shí)別率和識(shí)別時(shí)間,并且也對(duì)比了程序優(yōu)化對(duì)算法的改善。樣本數(shù)目和測(cè)試數(shù)目都為200,實(shí)驗(yàn)結(jié)果如表1 所示。從表1 可以看出,在兩種平臺(tái)之上運(yùn)行KPCA 算法,其識(shí)別率沒有變化,均為93%,但是DSP 平臺(tái)在識(shí)別時(shí)間方面大為減小,PC 平臺(tái)平均消耗時(shí)間為7.8 s,而DSP 平臺(tái)的平均消耗時(shí)間在優(yōu)化前是5.0 s,而優(yōu)化后僅僅需要1.6 s,識(shí)別速度提高了68%,效率顯著提高。

本文也比較了在DSP 平臺(tái)上不同的樣本數(shù)目對(duì)算法移植的影響。表2 對(duì)比了樣本數(shù)目分別為100 和200 時(shí)算法的優(yōu)化效果。從表2 可以看出,樣本數(shù)目越大, 優(yōu)化效果越明顯,主要原因在于DSP在處理大量數(shù)據(jù)的時(shí)候優(yōu)勢(shì)更為明顯。

表2 DSP 平臺(tái)不同樣本數(shù)目的識(shí)別時(shí)間對(duì)比Tab le 2 Recognition rate comparison between different sample numbers

5 結(jié)束語

本文將基于核主元分析的人臉識(shí)別算法移植到TI 公司的TMS320C6201 DSP 平臺(tái),并進(jìn)行了相關(guān)優(yōu)化。經(jīng)過優(yōu)化后的移植算法在保證識(shí)別率的前提下,節(jié)約了大量的識(shí)別時(shí)間,效率優(yōu)勢(shì)明顯,而且硬件系統(tǒng)運(yùn)行穩(wěn)定,可擴(kuò)展性強(qiáng),有一定的參考價(jià)值。如何在DSP 平臺(tái)上進(jìn)一步提高人臉識(shí)別率和運(yùn)行效率是下一步研究工作的重點(diǎn)。

[ 1] 俞寧.一種用于人臉識(shí)別的新PCA 算法[J] .電訊技術(shù),2009,49(12):65-69.

YU Ning.A Novel PCA Algorithm For Face Recognition[ J] .Telecommunication Engineering, 2009, 49(12):65 -69.(in Chinese)

[2] Scholkopf B, Smola A.NonlinearComponent Analysis as a Kernel Eigenvalue Problem[ J] .Neural Computer, 1998, 10(2):1299-1319.

[3] Gu Y F, Liu Y, Zhang Y.A Selective KPCA Algorithm based on H igh -order Statistics for Anomaly Detection in Hyperspectral Imagery[ J] .IEEE Geoscience and Remote Sensing Letters,2008, 5(1):43-47.

[4] LIU Chang-ping.Gabor based kernel PCA with fractional power polynom ial models for face recognition[ J] .IEEE Transactions on Pattern Analysis and M achine Intelligence, 2004,26(5):572-581.

[5] 杜卓明, 屠宏,耿國(guó)華.KPCA 方法過程研究與應(yīng)用[ J] .計(jì)算機(jī)工程與應(yīng)用, 2010,46(7):8-10.

DU Zhuo-ming,TU Hong,GENG Guo-hua.KPCA Method Research and Application Process[ J] .Computer Engineering and Applications,2010,46(7):8-10.(in Chinese)

[6] 舒紅霞, 王繼紅.分布式DSP 實(shí)時(shí)操作系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J] .小型微型計(jì)算機(jī)系統(tǒng),2010,31(7):1407-1410.

SHU Hong-xia,WANG Ji-hong.Design and Implementation of a Distributed DSP Real Time Operating System[ J] .Mini-Micro Systems,2010,31(7):1407-1410.(in Chinese)

[7] Eyre J, Bier J.The Evolution of DSP Processors From Early Architectures to The Latest Developments [ J] .IEEE Signal Processing Magazine, 2000, 17(2):43-51.

[8] 鄒土土土, 張超.基于DSP 的人臉識(shí)別算法實(shí)現(xiàn)與優(yōu)化[ J] .計(jì)算機(jī)應(yīng)用,2010,30(3):854-856.

ZOU Yao, ZHANG Chao.Implementation and Optimization of Face Recognition Algorithm based on DSP[ J] .Journal of Computer Applications,2010, 30(3):854-856.(in Chinese)