吳子揚(yáng), 周大可, 楊 欣, 胡陽(yáng)明
(南京航空航天大學(xué) 自動(dòng)化學(xué)院, 南京 210016)
近年來(lái), 由于計(jì)算機(jī)圖形學(xué)和模式識(shí)別理論的不斷發(fā)展, 基于人臉面部的應(yīng)用越來(lái)越受到關(guān)注。隨著人臉檢測(cè)、 人臉識(shí)別等應(yīng)用日益商業(yè)化, 人們開(kāi)始認(rèn)識(shí)到以圖片形式表達(dá)的人臉信息在實(shí)際處理中普遍受到光照、 姿態(tài)等多方面的制約。因此, 對(duì)人臉計(jì)算機(jī)表達(dá)的相關(guān)研究也由傳統(tǒng)的二維空間轉(zhuǎn)向三維空間。由于直接獲取人臉的三維數(shù)據(jù)需要昂貴的硬件, 故由人臉的二維圖像恢復(fù)其三維信息普遍成為當(dāng)前研究的熱點(diǎn)[1]。
由面部圖像恢復(fù)三維信息的經(jīng)典方法是Blanz等[2,3]和Patel等[4]提出的形變模型算法, 通過(guò)對(duì)一系列的三維人臉樣本進(jìn)行統(tǒng)計(jì)計(jì)算參考人臉, 求解全局優(yōu)化方程實(shí)現(xiàn)參考人臉到待重建樣本的匹配, 其本質(zhì)是線性組合思想, 即通過(guò)某一子空間的基表達(dá)該空間的具體實(shí)例。由于形變模型方法需要通過(guò)迭代求得全局最優(yōu)解, 故需要在局部細(xì)節(jié)和重建速度上做出一定的犧牲, 王成章等[5]提出了利用改進(jìn)的遺傳算法求解形變因子, 在形變模型的基礎(chǔ)上節(jié)省了11 min以上。Blanz等[6]和胡陽(yáng)明等[7]提出了通過(guò)少量特征點(diǎn)重建三維模型的方法, 即稀疏形變模型法。稀疏形變模型同樣利用線性組合思想, 通過(guò)對(duì)面部少量特征點(diǎn)近似求解優(yōu)化方程, 得到形變系數(shù), 從而完成參考模型的形變。其優(yōu)勢(shì)在于重建速度快, 且重建精度不會(huì)損失太多。后續(xù)的工作主要建立在統(tǒng)計(jì)得到的參考模型[8-10]之上, 通過(guò)不同方法對(duì)參考模型進(jìn)行形變或修改以得到最終重建結(jié)果。
Liao等[11]將局部形變模型與明暗恢復(fù)(SFS: Shadow From Shading)方法結(jié)合, 對(duì)人臉的局部區(qū)域分別重建, 通過(guò)RBF(Radical Basis Function)徑向基插值得到最終結(jié)果。但由于局部形變模型較整體形變模型而言更易受到訓(xùn)練樣本的影響, 且該方法需要首先在二維圖像上檢測(cè)出關(guān)鍵區(qū)域, 故統(tǒng)計(jì)得到每個(gè)區(qū)域的形變模型較困難。文獻(xiàn)[12,13]提出通過(guò)參考模型和待重建人臉面部灰度值計(jì)算光照反射率, 通過(guò)迭代尋優(yōu)的方法恢復(fù)整體的深度值, 但由于需要比較參考模型和待重建模型的灰度值, 該方法在很大程度上會(huì)受到光照條件的影響。為了解決單純的稀疏形變模型對(duì)人臉模型細(xì)節(jié)描述能力不夠[14]這一缺陷, 龔勛[15]提出了雙形變模型方法, 即將兩種弱形變模型結(jié)果通過(guò)加權(quán)方式結(jié)合得到最終三維重建結(jié)果。該方法對(duì)細(xì)節(jié)的表達(dá)較為豐富, 但不足在于權(quán)值的計(jì)算方法以及兩種弱形變模型參數(shù)的確定。此外, 還有部分研究集中于對(duì)原始形變模型優(yōu)化方程的改進(jìn)[16,17]。
筆者在前述工作的基礎(chǔ)上, 通過(guò)對(duì)稀疏形變模型的研究, 提出了一種由粗到精的重建方案, 即分步重建。首先通過(guò)對(duì)部分點(diǎn)的稀疏重建得到一個(gè)外形相似的粗略模型, 替代原有參考模型, 并將剩余點(diǎn)根據(jù)新的參考模型進(jìn)行第2步重建。實(shí)驗(yàn)表明, 該方法在不降低重建速度的基礎(chǔ)上能有效提高重建精度, 為三維人臉重建的研究提供一定的借鑒, 也為后續(xù)的三維人臉識(shí)別打下基礎(chǔ)。
將三維人臉子空間Rn視為線性子空間, 其中任意一個(gè)元素都可由其余元素線性表示, 即
(1)
通常情況下, 采集到的三維人臉數(shù)據(jù)由點(diǎn)云組成, 而每個(gè)三維人臉的形狀信息也可用點(diǎn)云的坐標(biāo)表示, 即Si=(x1,y1,z1,x2,y2,z2,…,xk,yk,zk)T, 其中k為點(diǎn)云的數(shù)目。如果直接使用所有三維人臉數(shù)據(jù)重構(gòu)新的三維人臉的計(jì)算量巨大, 一般采用主成分分析(PCA: Principal Component Analysis)理論對(duì)樣本數(shù)據(jù)庫(kù)進(jìn)行降維, 得到
(2)
根據(jù)式(2)可知, 求取β需要知道測(cè)試人臉三維模型上點(diǎn)云的坐標(biāo)信息, 但實(shí)際測(cè)量對(duì)象是二維圖片, 僅能獲取二維人臉上若干點(diǎn)的位置信息。文獻(xiàn)[6]給出了由少量二維點(diǎn)數(shù)據(jù)恢復(fù)整體三維人臉信息方法, 即稀疏形變模型理論。
設(shè)測(cè)量得到的二維稀疏坐標(biāo)信息為d=(x1,y1,x2,y2,…,xl,yl)T∈R2l, 其中l(wèi)為二維取點(diǎn)的數(shù)目。設(shè)該人臉三維模型全體點(diǎn)坐標(biāo)數(shù)據(jù)D到二維稀疏點(diǎn)坐標(biāo)的映射關(guān)系為L(zhǎng), 則d=LD, 根據(jù)式(2), 可做出以下推導(dǎo)
(3)
LΔS=Δd=LPβ=Pfβ
(4)
Δd=Pfβ
(5)
其中Δd為關(guān)鍵特征點(diǎn)的二維坐標(biāo)與對(duì)應(yīng)平均模型上三維點(diǎn)坐標(biāo)投影的距離,Pf代表從主成分中選取對(duì)應(yīng)特征點(diǎn)的二維坐標(biāo)集合。通常, 直接采用式(5)計(jì)算β會(huì)帶來(lái)過(guò)擬合現(xiàn)象, 故一般采用文獻(xiàn)[6]所述基于貝葉斯最大后驗(yàn)概率的最優(yōu)解法。令目標(biāo)函數(shù)為
E(β)=‖Pfβ-Δd‖2+η‖β‖2
(6)
其中η為權(quán)重系數(shù)。根據(jù)奇異值分解(SVD: Sigular Value Decomposition), 有Pf=UΛVT, 其中Λ為特征值構(gòu)成的對(duì)角矩陣,U由Pf(Pf)T的特征向量構(gòu)成,V由(Pf)TPf的特征向量構(gòu)成。這樣可得其最優(yōu)解為
(7)
將式(7)的計(jì)算結(jié)果帶入式(2), 可得到三維數(shù)據(jù)的重建結(jié)果
(8)
圖1給出擾動(dòng)系數(shù)η變化對(duì)重建結(jié)果的影響。由圖1可見(jiàn),η越大, 重建結(jié)果越接近平均人臉, 誤差越大;η過(guò)小會(huì)導(dǎo)致過(guò)匹配, 重建結(jié)果嚴(yán)重變形。
圖1 擾動(dòng)系數(shù)η對(duì)重建結(jié)果的影響
a 第1步重建特征點(diǎn) b 所有特征點(diǎn)
用于測(cè)試的二維人臉首先需要進(jìn)行特征點(diǎn)的定位, 并將定位后的二維坐標(biāo)與三維平均人臉坐標(biāo)建立對(duì)應(yīng)關(guān)系。這種對(duì)應(yīng)關(guān)系通常通過(guò)空間坐標(biāo)系到平面坐標(biāo)系的投影建立, 將二維坐標(biāo)轉(zhuǎn)換為圖像坐標(biāo), 其原點(diǎn)一般取鼻尖點(diǎn); 而將三維坐標(biāo)轉(zhuǎn)換為柱面坐標(biāo), 鼻尖點(diǎn)對(duì)應(yīng)于三維坐標(biāo)的(0,Ymax,0)。在實(shí)際操作中, 由于測(cè)試人臉受到光照不均或姿態(tài)不完全校正等影響, 無(wú)法準(zhǔn)確定位出所有特征點(diǎn), 特別是面部區(qū)域非顯著點(diǎn)。根據(jù)式(7), 所有點(diǎn)的信息均需要帶入運(yùn)算才能求得最優(yōu)形變因子, 如果存在定位不準(zhǔn)確的點(diǎn), 將導(dǎo)致求得的形變因子不準(zhǔn)確。因此, 在傳統(tǒng)的稀疏重建方法中, 二維面部每個(gè)特征點(diǎn)定位的精度都會(huì)直接影響重建結(jié)果的精度。
為克服一次重建存在的問(wèn)題, 筆者采用三維人臉稀疏分步重建的方法。利用由粗到精的思路, 先根據(jù)面部定位較準(zhǔn)確的關(guān)鍵特征點(diǎn)進(jìn)行重建, 得到重建的粗略模型, 并根據(jù)第1步重建結(jié)果進(jìn)行第2步重建。兩步中選取的特征點(diǎn)如圖2所示。
由式(8)可得到由上述特征點(diǎn)重建后的第1步重建結(jié)果
(9)
其中η1為第1步重建的擾動(dòng)因子。
將S1作為第2步重建的參考人臉, 可得到其余特征點(diǎn)的二維坐標(biāo)與新的參考人臉三維坐標(biāo)投影的距離Δd′。在此基礎(chǔ)上將Δd′帶入式(6), 可得
(10)
其中η2為第2步重建的擾動(dòng)因子??梢缘玫椒植街亟ǖ淖罱K結(jié)果
(11)
在第1步重建過(guò)程中, 選擇每個(gè)區(qū)域的顯著點(diǎn), 如眼角、 鼻翼和嘴角點(diǎn)作為重建特征點(diǎn), 以確定三維人臉的大體形狀和各個(gè)區(qū)域的位置。由于選取的點(diǎn)較少, 一方面雖然重建出的模型細(xì)節(jié)表達(dá)不完善, 但能大體確定出待重建人臉的形狀和區(qū)域的位置; 另一方面由于每個(gè)區(qū)域只選取顯著的特征點(diǎn), 避免了同一區(qū)域非顯著點(diǎn)對(duì)顯著點(diǎn)的干擾。在第2步重建過(guò)程中, 由于已得到了一個(gè)粗略的形狀模型, 故需要選取較多的特征點(diǎn)對(duì)該模型進(jìn)行形變, 以達(dá)到對(duì)模型細(xì)節(jié)的重建。
在兩步重建的過(guò)程中, 擾動(dòng)系數(shù)的選擇也應(yīng)采用不同的方式。第1步重建是為得到大體形狀以及區(qū)域位置, 應(yīng)選取較小的擾動(dòng)系數(shù)以減小加入的干擾; 第2步重建是為得到較精確的模型, 應(yīng)選取使重建結(jié)果誤差較小的擾動(dòng)系數(shù)。
筆者提出的分步重建算法流程如圖3所示。
圖3 分步重建算法流程
實(shí)驗(yàn)以BJUT-3D Face Database中的部分三維人臉數(shù)據(jù)為訓(xùn)練樣本, 測(cè)試樣本為該庫(kù)其余樣本、 部分CAS-PEAL二維人臉數(shù)據(jù)庫(kù)以及部分現(xiàn)實(shí)照片。實(shí)驗(yàn)采用的硬件環(huán)境為主頻1.80 GHz的雙核Intel CPU和2.0 GByte內(nèi)存, 軟件環(huán)境為Matlab 2012a和VC++6.0。
實(shí)驗(yàn)分為兩部分。
1) 庫(kù)內(nèi)測(cè)試。訓(xùn)練樣本和測(cè)試樣本均來(lái)自BJUT-3D數(shù)據(jù)庫(kù), 選取已規(guī)范化(每個(gè)人臉由23 676個(gè)點(diǎn)構(gòu)成)的70個(gè)三維人臉數(shù)據(jù)作為訓(xùn)練集, 30個(gè)三維人臉數(shù)據(jù)作為測(cè)試集。將測(cè)試集樣本通過(guò)正交投影得到二維數(shù)據(jù)庫(kù), 直接在該二維庫(kù)上重建, 并與原三維信息進(jìn)行比較。
2) 庫(kù)外測(cè)試。即選取已規(guī)范化的100個(gè)三維人臉數(shù)據(jù)作為訓(xùn)練集, 直接對(duì)實(shí)際照片進(jìn)行重建。實(shí)際照片來(lái)自CAS-PEAL二維人臉數(shù)據(jù)庫(kù)。
實(shí)驗(yàn)共選取人面部47個(gè)特征點(diǎn), 對(duì)筆者提出的三維人臉稀疏分步重建與基于稀疏形變模型的三維重建分別進(jìn)行實(shí)驗(yàn), 并與現(xiàn)有相關(guān)三維重建方法做比較。
對(duì)于三維人臉數(shù)據(jù)庫(kù)內(nèi)樣本, 由于已知其真實(shí)三維信息, 故實(shí)驗(yàn)采用平均歐氏距離[6]作為重建精度評(píng)價(jià)標(biāo)準(zhǔn), 將重建得到的三維數(shù)據(jù)與原始數(shù)據(jù)直接比較。重建精度(重建誤差)計(jì)算如下
(12)
其中So,Sr分別是兩個(gè)形狀的曲面,V是曲面上的頂點(diǎn)集合,Vo,Vr分別是So,Sr上的頂點(diǎn)。eeuc越小, 則代表重建結(jié)果與真實(shí)結(jié)果越吻合, 即重建精度越高, 重建誤差越小。
實(shí)驗(yàn)第1步共選取人面部12個(gè)關(guān)鍵特征點(diǎn)以及3對(duì)位于邊緣輪廓上的點(diǎn), 初步確定待重建人臉的臉型與關(guān)鍵區(qū)域的大體位置; 第2步采用面部其余特征點(diǎn), 以確定關(guān)鍵區(qū)域的具體細(xì)節(jié)信息。對(duì)兩步重建中的擾動(dòng)因子分別取值, 計(jì)算其最終重建誤差。經(jīng)過(guò)多次試驗(yàn),η1=0.85,η2=3.4時(shí), 平均重建誤差最小。故對(duì)于測(cè)試集所有樣本, 均取η1=0.85,η2=3.4。選用稀疏形變模型(SRSD: Surface Reconstruction from Sparse Data)重建算法作比較, 部分實(shí)驗(yàn)結(jié)果如圖4所示。
a 二維投影人臉 b 第1步結(jié)果 c 第2步結(jié)果 d 真實(shí)三維信息 e SRSD重建結(jié)果
由圖4可看出, 分步形變的結(jié)果更加接近真實(shí)三維人臉, 且在面部細(xì)節(jié)上與稀疏形變模型相比表現(xiàn)得更加細(xì)致。而稀疏形變模型的結(jié)果普遍與第1次形變的結(jié)果較為相似, 這是由于稀疏形變模型所選用的部分特征點(diǎn)對(duì)整個(gè)重建結(jié)果的影響較小, 僅采用人臉部分關(guān)鍵特征點(diǎn)即可形變得到初步的三維人臉重建模型。
表1 相關(guān)算法的重建誤差比較
筆者選擇與筆者相關(guān)的基于稀疏形變模型的三維人臉重建算法和文獻(xiàn)[17]提出的SDM(Sparse Deforming Model)算法作比較, 分別計(jì)算其歐氏距離, 最終結(jié)果為測(cè)試集上30個(gè)樣本的誤差均值。精度對(duì)比如表1所示。其中SRSD與SDM的保留主成分?jǐn)?shù)目為69, 其余參數(shù)設(shè)置均參照文獻(xiàn)[5]和文獻(xiàn)[15]。
結(jié)合圖3與重建誤差實(shí)驗(yàn)比較可以得出, 筆者算法與SRSD相比, 精度提高了約10%, 與SDM相比, 精度提高了1.85 mm。
實(shí)驗(yàn)同時(shí)比較了不同算法的速度, 將常見(jiàn)三維人臉重建算法所需時(shí)間與筆者所述算法進(jìn)行對(duì)比。其中SRSD和筆者算法均用Matlab 2012a實(shí)現(xiàn), 包含用于主成分分析的時(shí)間, 其余算法時(shí)間分別引自文獻(xiàn)[2]、 文獻(xiàn)[6]、 文獻(xiàn)[12]和文獻(xiàn)[8]。實(shí)驗(yàn)結(jié)果如表2所示。
表2 相關(guān)算法的重建速度比較
從相關(guān)算法的速度比較可看出, 筆者提出的分步三維重建算法速度與傳統(tǒng)的稀疏形變模相比, 盡管在第1輪重建結(jié)果的基礎(chǔ)上引入第2輪重建, 但Matlab實(shí)現(xiàn)的算法速度與之相當(dāng)(稍慢0.02 s)。與經(jīng)典的形變模型算法相比, 筆者提出的算法速度具有很大的提高。綜合重建精度和重建速度比較, 筆者提出的基于分步重建的三維人臉稀疏重建算法具有一定的優(yōu)勢(shì)。
該實(shí)驗(yàn)采用CAS-PEAL二維人臉數(shù)據(jù)庫(kù)中POSE子集中隨機(jī)抽取的個(gè)人照片作為測(cè)試庫(kù)。對(duì)于庫(kù)中的圖片, 實(shí)驗(yàn)對(duì)正臉照片進(jìn)行重建, 重建的結(jié)果以不同角度的側(cè)面與數(shù)據(jù)庫(kù)中對(duì)應(yīng)角度的照片進(jìn)行比較, 部分實(shí)驗(yàn)結(jié)果如圖5所示。
a 輸入2D圖像 b 三維重建結(jié)果 c 輸入圖像側(cè)面 d 重建結(jié)果側(cè)面 e 輸入圖像側(cè)面 f 重建結(jié)果側(cè)面
對(duì)真實(shí)照片的重建結(jié)果可看出, 筆者提出的方法能一定程度上還原二維人臉的三維信息, 相比于傳統(tǒng)的重建方法, 筆者方法能在重建精度上得到一定的提高。但由于訓(xùn)練樣本庫(kù)較小, 且方法的根本是由少量點(diǎn)控制整體形變等原因, 三維重建的結(jié)果仍然在細(xì)節(jié)上存在一定誤差。
筆者在傳統(tǒng)稀疏形變模型三維人臉重建方法的基礎(chǔ)上, 提出了一種基于分步重建的稀疏三維重建方法, 并在BJUT-3D Face Database和CAS-PEAL數(shù)據(jù)庫(kù)中做了測(cè)試。實(shí)驗(yàn)表明, 該算法能在重建速度較快的基礎(chǔ)上, 與經(jīng)典算法相比提高一定的重建精度。由于采用兩步重建三維人臉, 每步的特征點(diǎn)數(shù)量的選擇以及參數(shù)設(shè)置都會(huì)影響最終結(jié)果和重建精度, 因此, 相關(guān)參數(shù)的選擇和優(yōu)化將是進(jìn)一步研究的重點(diǎn)。
參考文獻(xiàn):
[1]徐琳, 袁保宗, 高文. 真實(shí)感人臉建模研究的進(jìn)展與展望 [J]. 軟件學(xué)報(bào), 2003, 14(4): 804-810.
XU Lin, YUAN Bao-zong, GAO Wen. Development and Prospect on Realistic Facial Modeling [J]. Journal of Software, 2003, 14(4): 804-810.
[2]BLANZ V, VETTER T. A Morphable Model for the Synthesis of 3D Faces [C]∥ Computer Graphics Proceedings, SIGGRAPH’99. New York, USA: [s.n.], 1999: 187-194.
[3]BLANZ V, VETTER T. Face Recognition Based on Fitting a 3D Morphable Model [J]. IEEE Trans Pattern Analysis and Machine Intelligence, 2003, 25(9): 1063-1074.
[4]PATEL A, SMITH W A P. 3D Morphable Face Models Revisited [C]∥Proc IEEE Conf Computer Vision and Pattern Recognition. Miami, FL, USA: IEEE, 2009: 1327-1334.
[5]王成章, 尹寶才, 孫艷豐, 等. 改進(jìn)的基于形變模型的三維人臉建模方法 [J]. 自動(dòng)化學(xué)報(bào), 2007, 33(3): 232-239.
WANG Cheng-zhang, YIN Bao-cai, SUN Yan-feng, et al. An Improved 3D Face Modeling Method Based on Morphable Model [J]. Acta Automatica Sinica, 2007, 33(3): 232-239.
[6]BLANZ V, VETTER T. Reconstructing the Complete 3D Shape of Faces from Partial Information [J]. It-Information Technology, 2002, 44(6): 295-302.
[7]胡陽(yáng)明, 周大可, 鹿樂(lè), 等. 基于形變模型的三維人臉快速重建改進(jìn)算法 [J]. 吉林大學(xué)學(xué)報(bào): 信息科學(xué)版, 2012, 30(6): 605-609.
HU Yang-ming, ZHOU Da-ke, LU Le, et al. Improved 3D Fast Face Reconstruction Algorithm Based on Morphable Model [J]. Journal of Jilin University: Information Science Edition, 2012, 30(6): 605-609.
[8]胡永利. 真實(shí)感三維人臉建模及應(yīng)用研究 [D]. 北京: 北京工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院, 2004.
HU Yong-li. Realistic 3D Face Modeling and Its Application [D]. Beijing: School of Computer Science and Technology, Beijing University of Technology, 2004.
[9]PATEL A, SMITH W A P. Driving 3D Morphable Models Using Shading Cues [J]. Pattern Recognition, 2012, 45(5): 1993-2004.
[10]尹寶才, 何晏晏, 孫艷豐, 等. 三維人臉的非均勻重采樣對(duì)齊算法 [J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2007, 33(2): 213-218.
YIN Bao-cai, HE Yan-yan, SUN Yan-feng, et al. Nonuniform Resampling Based on Method for Pixel-Wise Correspondence between 3D Faces [J]. Journal of Beijing University of Technology, 2007, 33(2): 213-218.
[11]LIAO H B, CHEN Q H, ZHOU Q J, et al. Rapid 3D Face Reconstruction by Fusion of SFS and Local Morphable Model [J]. Journal of Visual Communication and Image Representation, 2012, 23(6): 924-931.
[12]KEMELMACHER-SHLIZERMAN I, BASRI R. 3D Face Reconstruction from a Single Image Using a Single Reference Face Shape [J]. IEEE Trans Pattern Analysis and Machine Intelligence, 2011, 33(2): 394-405.
[13]KEMELMACHER-SHLIZERMAN I, BASRI R. Molding Face Shapes by Example [C]∥Proc European Conference on Computer Vision. Graz, Austn: [s.n.], 2006: 277-288.
[14]龔勛, 王國(guó)胤. 基于特征點(diǎn)的三維人臉形變模型 [J]. 軟件學(xué)報(bào), 2009, 20(3): 724-733.
GONG Xun, WANG Guo-yin. 3D Face Deformable Model Based on Feature Points [J]. Journal of Software, 2009, 20(3): 724-733.
[15]龔勛. 基于單張二維圖片的三維人臉建模 [D]. 重慶: 西南交通大學(xué)信息科學(xué)與技術(shù)學(xué)院, 2008.
GONG Xun. 3D Face Modeling Based on a 2D Image [D]. Chongqing: School of Information Science and Technology, Southwest Jiaotong University, 2008.
[16]徐成華, 王蘊(yùn)紅, 譚鐵牛. 三維人臉建模與應(yīng)用 [J]. 中國(guó)圖像圖形學(xué)報(bào), 2004, 9(8): 893-903.
XU Cheng-hua, WANG Yun-hong, TAN Tie-niu. Overview of Research on 3D Face Modeling [J]. Journal of Image and Graphics, 2004, 9(8): 893-903.
[17]柴秀娟, 山世光, 卿來(lái)云, 等. 基于3D人臉重建的光照、 姿態(tài)不變?nèi)四樧R(shí)別 [J]. 軟件學(xué)報(bào), 2006, 17(3): 525-534.
CHAI Xiu-juan, SHAN Shi-guang, QING Lai-yun, et al. Pose and Illumination Invariant Face Recognition Based on 3D Face Reconstruction [J]. Journal of Software, 2006, 17(3): 525-534.