DIP---基于MATLAB的人臉識別算法課程設(shè)計報告

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)字圖像處理實驗報告院 系： 計算機科學(xué)學(xué)院 班 級： 計科11303 小組成員： 張世柳、鄧偉養(yǎng)、蘭洋、馮威 成員學(xué)號： 實驗名稱： 基于MATLAB的人臉識別算法 實驗時間： 2015.10.01 - 2015.10.19 實驗地點： 東4教2號機房 專心-專注-專業(yè)目 錄一、緒論 隨著科技的發(fā)展，人類社會的進步，傳統(tǒng)身份識別由于容易遺失，容易被破解已不能起到身份識別作用。人們需要更加安全可靠的身份識別技術(shù)。而生物特征的獨一無二，不易丟失和被復(fù)制的特性很好滿足了身份識別的需要。同時隨著計算機科學(xué)技術(shù)和生物醫(yī)學(xué)的發(fā)展使得利用生物特征識別成為了可能。在生物特征識別領(lǐng)域

2、，由于人臉識別的操作快速簡單，結(jié)果直觀，準(zhǔn)確可靠，不需要人的配合等優(yōu)點已成為人們關(guān)注的焦點。主成分分析（PCA）通過提取高維度的人臉圖像的主元，使得圖像在低維度空間中被處理來降低了圖像處理的難度。由于其有效的解決了圖像空間維數(shù)過高的問題，已經(jīng)成為人臉識別領(lǐng)域非常重要的理論，本實驗研究的是基于MATLAB人臉識別算法的實現(xiàn)。二、實驗設(shè)計（一）實驗題目機器人視覺基于MATLAB的人臉識別算法（二）實驗?zāi)康?. 初步了解人臉識別的特征法；2. 學(xué)會使用主成分分析算法（PCA）；3. 通過功能模塊實現(xiàn)人臉識別系統(tǒng)；4. 完成數(shù)字圖像處理課程的作業(yè)要求。三、實驗準(zhǔn)備（一）環(huán)境準(zhǔn)備MATLAB 7.0（二

3、）知識準(zhǔn)備1. MATLAB的優(yōu)勢特點：（1) 高效的數(shù)值計算及符號計算功能，能使用戶從繁雜的數(shù)學(xué)運算分析中解脫出來；（2) 具有完備的圖形處理功能，實現(xiàn)計算結(jié)果和編程的可視化；（3) 友好的用戶界面及接近數(shù)學(xué)表達式的自然化語言，使學(xué)者易于學(xué)習(xí)和掌握；（4) 功能豐富的應(yīng)用工具箱，為用戶提供了大量方便實用的處理工具。2. 主成分分析（Principal Component Analysis，PCA），將多個變量通過線性變換以選出較少個數(shù)重要變量的一種多元統(tǒng)計分析方法。PCA是數(shù)字圖像處理中經(jīng)常用到的降維方法，在處理有關(guān)數(shù)字圖像處理方面的問題時，比如經(jīng)常用的圖像的查詢問題，在一個幾萬或者幾百萬甚

4、至更大的數(shù)據(jù)庫中查詢一幅相近的圖像。這時，我們通常的方法是對圖像庫中的圖片提取響應(yīng)的特征，如顏色，紋理，sift，surf，vlad等特征，然后將其保存，建立響應(yīng)的數(shù)據(jù)索引，再對要查詢的圖像提取相應(yīng)的特征，與數(shù)據(jù)庫中的圖像特征對比，找出與之最近的圖片。如果為了提高查詢的準(zhǔn)確率，通常會提取一些較為復(fù)雜的特征，如sift，surf等，一幅圖像有很多個這種特征點，每個特征點又有一個相應(yīng)的描述該特征點的128維的向量，設(shè)想如果一幅圖像有300個這種特征點，那么該幅圖像就有300*vector（128維）個，如果我們數(shù)據(jù)庫中有一百萬張圖片，這個存儲量是相當(dāng)大的，建立索引也很耗時，所以用PCA將其降維。四

5、、算法設(shè)計（一）問題描述1. 主成分的一般定義設(shè)有隨機變量X1，X2，Xp， 其樣本均數(shù)記為 ， ， ，樣本標(biāo)準(zhǔn)差記為S1，S2，Sp。首先作標(biāo)準(zhǔn)化變換，我們有如下的定義： (1) 若C1=a11x1+a12x2+ +a1pxp，且使 Var(C1)最大，則稱C1為第一主成分； (2) 若C2=a21x1+a22x2+a2pxp， ，(a21，a22，a2p)垂直于(a11，a12，a1p)，且使Var(C2)最大，則稱C2為第二主成分； (3) 類似地，可有第三、四、五主成分，至多有p個。2. 主成分的性質(zhì)主成分C1，C2，Cp具有如下幾個性質(zhì)： （1) 主成分間互不相關(guān)，即對任意i和j，C

6、i 和Cj的相關(guān)系數(shù) Corr(Ci，Cj)=0 i¹j(2) 組合系數(shù)(ai1，ai2，aip)構(gòu)成的向量為單位向量，(3) 各主成分的方差是依次遞減的， 即 Var(C1)Var(C2)Var(Cp) (4) 總方差不增不減， 即 Var(C1)+Var(C2)+ +Var(Cp) =Var(x1)+Var(x2)+ +Var(xp) =p 這一性質(zhì)說明，主成分是原變量的線性組合，是對原變量信息的一種改組，主成分不增加總信息量，也不減少總信息量。 (5) 主成分和原變量的相關(guān)系數(shù) Corr(Ci，xj)=aij =aij (6) 令X1，X2，Xp的相關(guān)矩陣為R, (ai1，ai

7、2，aip)則是相關(guān)矩陣R的第i個特征向量(eigenvector)。而且，特征值li就是第i主成分的方差， 即 Var(Ci)= li 其中l(wèi)i為相關(guān)矩陣R的第i個特征值(eigenvalue) ll1l2lp03. 主成分的數(shù)目的選取前面已指出，設(shè)有p個隨機變量，便有p個主成分。由于總方差不增不減，C1，C2等前幾個綜合變量的方差較大,而Cp，Cp-1等后幾個綜合變量的方差較小, 嚴(yán)格說來，只有前幾個綜合變量才稱得上主(要)成份，后幾個綜合變量實為“次”(要)成份。實踐中總是保留前幾個，忽略后幾個。 保留多少個主成分取決于保留部分的累積方差在方差總和中所占百分比(即累計貢獻率)，它標(biāo)志著前

8、幾個主成分概括信息之多寡。實踐中，粗略規(guī)定一個百分比便可決定保留幾個主成分；如果多留一個主成分，累積方差增加無幾，便不再多留。（二）PCA算法的功能實現(xiàn)1. 人臉空間的建立假設(shè)一幅人臉圖像包含N個像素點，它可以用一個N維向量表示。這樣，訓(xùn)練樣本庫就可以用i（i=1，.，M）表示。協(xié)方差矩陣C的正交特征向量就是組成人臉空間的基向量，即特征臉。將特征值由大到小排列：12.r,其對應(yīng)的特征向量為k。這樣每一幅人臉圖像都可以投影到由u1,u2,.,ur張成的子空間中。因此，每一幅人臉圖像對應(yīng)于子空間中的一點。同樣，子空間的任意一點也對應(yīng)于一幅圖像。2. 特征向量的選取雖然協(xié)方差矩陣最多有對應(yīng)于非零特征

9、值的k（k遠(yuǎn)小于M）個特征向量，但是通常情況下，k仍然很大，而事實上，根據(jù)應(yīng)用的要求，并非所有的特征向量都有需要保留，而特征空間投影的計算速度是直接與創(chuàng)建子空間所用的特征向量的數(shù)目相關(guān)，若考慮到計算時間的因素，可以適當(dāng)?shù)臏p去一些信息量少的特征向量，而且，去掉這些特征向量之后不一定不利于分類結(jié)果，有的情況下反而能夠提高識別性能。3. 人臉識別有了這樣一個由"特征臉"張成的降維子空間，任何一幅人臉圖像都可以向其投影得到一組坐標(biāo)系數(shù)，這組系數(shù)表明了該圖像在子空間中的位置，從而可以作為人臉識別的依據(jù)。換句話說，任何一幅人臉圖像都可以表示為這組"特征臉"的線性組合

10、，其加權(quán)系數(shù)即是K-L變換的展開系數(shù)，也可以稱為該圖像的代數(shù)特征。因此，在獲得特征臉之后，就可以對每一類別的典型樣本進行投影，由此得到每個人臉的投影特征從而構(gòu)成人臉特征向量，作為下一步識別匹配的搜索空間。輸入圖像及其在人臉空間上的投影（重構(gòu)圖像），人臉圖像在人臉空間中的投影變化不明顯而非人臉圖像的投影變化明顯。因此，檢測一幅圖像中是否存在人臉的基本思想是，計算該圖像中任意位置處的局部圖像與人臉空間之間的距離。其中，是局部圖像是否為人臉的度量。因此，計算給定圖像任意一點上的，就可以得到一映射圖（x,y）。4. 識別流程五、程序?qū)崿F(xiàn)（一）人臉識別程序1. 用戶界面2. 選擇圖片3. 圖片選擇后4.

11、 查找后（二）測試及結(jié)果分析實驗在兩個圖庫上測試，一個是自建人臉庫，該庫包含10個不同人物，每人有5張不同表情和姿態(tài)下的圖片，總共50幅。另一個是ORL人臉庫，該庫包含40個不同人物，每人有10張圖片，共400幅。用訓(xùn)練樣本進行測試，識別率為100%。而隨著訓(xùn)練樣本的增加，識別率會有所提升，由于標(biāo)準(zhǔn)人臉庫在采集時考慮了多種因素，人臉圖像比較標(biāo)準(zhǔn)，所以識別率較自建人臉庫識別率高，另外因為自建人臉庫的圖片太少，即訓(xùn)練樣本太少，也會對結(jié)果產(chǎn)生影響，效果不是很好。進行直方圖均衡化比灰度歸一化的識別率高，預(yù)處理對識別的效果起著至關(guān)重要的作用。而此次實驗的預(yù)處理還比較粗糙，PCA也只是起到了簡單的特征臉降

12、維的作用，要有更好的效果，還必須尋找更好的特征表達，使得可以盡量消除光照、表情、遮掩和姿勢的影響。（三）程序源碼1. “讀取圖片”按鈕% 讀取待查找圖片global im;%由于要在兩個按鈕函數(shù)中使用，故使用全局變量filename, pathname=. uigetfile('*.bmp','選擇圖片');str = pathname, filename;%合成路徑+文件名im = imread(str);%讀取圖片axes( handles.axes1);%使用第一個axesimshow(im);title('待查找')%顯示圖片2. “開始

13、查找”按鈕% PCA人臉識別global im;%使用全局變量imgdata=;%訓(xùn)練圖像矩陣for i=1:10 for j=1:5 a=imread(strcat('C:UsersdellDesktopORLs',num2str(i),'',num2str(j),'.bmp'); b=a(1:112*92); % b是行矢量 1×N，其中N10304 b=double(b); imgdata=imgdata; b; % imgdata 是一個M * N 矩陣，imgdata中每一行數(shù)據(jù)一張圖片，M50 end;end;imgdata

14、=imgdata' %每一列為一張圖片imgmean=mean(imgdata,2); % 平均圖片，N維列向量for i=1:50 minus(:,i) = imgdata(:,i)-imgmean; % minus是一個N*M矩陣，是訓(xùn)練圖和平均圖之間的差值end;covx=minus'* minus; % M * M 階協(xié)方差矩陣COEFF, latent,explained = pcacov(covx'); %PCA，用協(xié)方差矩陣的轉(zhuǎn)置來計算以減小計算量% 選擇構(gòu)成95%的能量的特征值i=1;proportion=0;while(proportion <

15、95) proportion=proportion+explained(i); i=i+1;end;p=i-1;% 訓(xùn)練得到特征臉坐標(biāo)系i=1;while (i<=p && latent(i)>0) base(:,i) = latent(i)(-1/2)*minus * COEFF(:,i); % base是N×p階矩陣，用來進行投影，除以latent(i)(1/2)是對人臉圖像的標(biāo)準(zhǔn)化 i = i + 1;end% 將訓(xùn)練樣本對坐標(biāo)系上進行投影,得到一個 p*M 階矩陣為參考reference = base'*minus;% 測試過程在測試圖片文

16、件夾中選擇圖片，進行查找測試a=im;b=a(1:10304);b=double(b);b=b'object = base'*(b-imgmean);distance=;%最小距離法，尋找和待識別圖片最為接近的訓(xùn)練圖片for k=1:50 temp= norm(object - reference(:,k); if (distance > temp) which = k; distance = temp; end;end;%找出距離最近的圖片所在的位置num1 = ceil(which/5);%第num1個文件夾num2 = mod(which,5);%第num2個圖片文

17、件if (num2 = 0) num2 = 5;end;I=imread(strcat('C:UsersdellDesktopFaces',num2str(num1),'',num2str(num2),'.bmp');%讀取該圖片axes( handles.axes2);%使用第2個axes%輸出判斷所最接近的訓(xùn)練樣本的人臉圖片并顯示其人名switch num1 case 1 imshow(I);title('此人為 Peter');%顯示圖片及人名 case 2 imshow(I);title('此人為 Ben'

18、); case 3 imshow(I);title('此人為 Nick'); case 4 imshow(I);title('此人為 Park'); case 5 imshow(I);title('此人為 Linkin'); case 6 imshow(I);title('此人為 Andy'); case 7 imshow(I);title('此人為 Arthur'); case 8 imshow(I);title('此人為 Lisa'); case 9 imshow(I);title('此人為 Dana'); case 10 imshow(I);title('此人為