基于FPGA的混沌加密虹膜識別系統設計(二)

　　7.1.3 虹膜外邊緣的確定

　　(1) 虹膜外邊緣的特征分析

　　由圖1中所示的虹膜圖像可以看出，虹膜外邊緣的主要特點是：較相對與虹膜內邊緣而言，邊緣處灰度變化不是特別明顯，有一小段漸變的區域。也就是說，虹膜內部灰度趨近于一致這個事實，在參考文獻[8]中，介紹的環量積分算子應該式是一種有效的方法。

　　即：

?

　　(7-10)

　　(2) 采用環量積分算子實現虹膜外邊緣的檢測

　　如上分析，虹膜環量積分算子是檢測虹膜外邊緣的一種有效手段，為了克服虹膜紋理對環量線積分的影響，本文對式(7-1)作了如下改進，將環量線積分

　　改變為求

?

　　的圓環狀區域的面積分。即：

?

　　(7-11)

　　各符號的意義與(7-1)相同，為了便于計算，將其離散化可得:

?

　　(7-12)

　　式中

?

　　分別為t，

?

　　增長的步長，n，k，l，m分別為求卷積時高斯函數的中心、圓環中心、圓環的寬度、旋轉角度等參數的增量。顯而易見，式(7-12)和式(7-11)并不完全等價，式(7-11)中積分號內的部分的意義為積分區域內各個點的灰度的平均值，式(7-12)計算的是積分區域各個點灰度的總和，但由于(7-12)中角度

?

　　的步長

?

　　和圓環寬度t 的步長

?

　　以及圓環的寬度t 都是固定的，也就是說，對應不同的的圓環，從其中提取的計算環量積分的點的個數都是固定的，因而兩者只相差一個比例常數，并不影響判斷。注意，由式(7-1)到式(7-11)的改進過程中用到了卷積的性質：

?

　　(7-13)

　　如果使用式(7-12)在整個圖像空間中搜索，則系統開銷過大，本文將充分利用已經求得的瞳孔中心的位置參量，設定虹膜外邊緣的圓心與瞳孔中心相差5 各像素，從而將在整個區域內的搜索簡化為在5×5 的矩形區域內的搜索，大大減少了算法的時間復雜度。考慮到虹膜圖像的上部和下部易于受到眼皮和睫毛的干擾，在計算環量積分時，

?

　　取值限定在

?

　　的范圍內。

　　為了進一步減少系統搜索的開銷，本文采用一種由粗到精的取點與計算方法，設點

?

　　為搜索點，以t 為半徑增量，依次計算式(7-11)所示的環量面積分算子，在搜索空間內求得通過環量面積分的的最大值初步確定圓心和半徑

?

　　之后，再以

?

　　為圓心，在

?

　　區域內使用式(7-1)精確搜索，以確定圓的精確大小。

　　7.2 虹膜圖像的展開

　　為了便于對虹膜圖片分析，一般的系統中都要將其展開成矩形。

　　我們采用內圓圓心為中心，以虹膜的寬度為半徑建立極坐標系，將虹膜在極坐標系(ρ，θ)下展開成為橫坐標為θ，縱坐標為ρ的720*50 的矩形區域，展開的過程中，必然會出現新的圖像中某些點無法與原圖像中的點進行匹配的情況，通常情況下應進行插值處理，一般情況下，插值有以下幾種方法：(1)0級內插法，即將該點周圍四個鄰點中離它最近的一個點的像素的灰度級做為它的灰度級。(2)1 級內插法，亦稱雙線性內插法，是根據周圍四個點的灰度在兩個方向上進行線性內插，從而對原圖像中不存在的點計算出其近似值而不是用其鄰近點的像素來代替。(3)三次卷積法，是利用多項式來逼近理論上的最佳插值函數

?

　　的方法。由于0 級插值法缺乏一定的精度，而三次卷積法又計算量過大，本文中采取雙線性內插法。使用極坐標的優點是：眼睛旋轉的變化，可以轉變為θ方向的平移。

　　7.3 虹膜圖像的二值化

　　如上文所述，虹膜表面有許多斑點、凹陷區和皺紋組成，這些特征形成與遺傳和胚胎發育過程，含有豐富的信息。而且終生不變，從這些信息中可以用不同角度用不同的方法提取出用于區分不同虹膜的特征，進而進行身份識別。

　　本文采用虹膜圖像的結構特征分析方法進行虹膜識別。結構特征通常包括控制點，角，線段等等，結構特征具有直觀性好、穩定性高、抗噪聲能力強、編碼效率高等優點。利用結構特征進行編碼可以方便地解決虹膜圖像殘缺問題和局部編碼的區域劃分誤差。由于結構特征的直觀性，便于形成統一的數據格式，有利于應用的推廣。在灰度圖中，虹膜區域的紋理特征表現為虹膜區域內的灰度變化，記錄這些灰度變化，對圖像進行二值化處理，將圖像背景和灰度急劇變化的區域分開，可以作為虹膜識別的依據，下面就探討虹膜圖像二值化的方法。

　　7.3.1 虹膜圖像二值化方法探討

　　從直觀來看，提取圖像灰度變化的方法只需設定一定的閾值就提取其變換的信息，但是這種方法在提取灰度變化信息時卻有一定的局限性。如圖9所示為一圖像的截面圖。橫坐標表示截面的伸展方向，縱坐標表示對應點的灰度，若取閾值為A，則BC的之間的灰度變化體現不出來，同理，若取閾值為C，則無法體現AB之間的閾值變化。

　　7.3.2 邊緣檢測與濾波器的選擇

　　對于灰度圖像中的各點，其灰度值的一階或二階導數能夠很好地體現圖像邊界點，本文將探討用二階導數來求邊緣點，以反映圖像灰度的變化，從而確定特征點，但是由于噪聲信號的影響，一般應先對信號進行平滑濾波，設信號g(x)，如平滑濾波器的沖擊響應函數用h(x)表示，則濾波后的信號為g(x)=f(x)*h(x),然后再對g(x)求二階導數以檢測邊緣點。

　　由于微分運算與卷積運算次序有以下互換關系：

?

　　因此可以先平滑，后微分的兩步運算合并，并將平滑平滑微分濾波器的導數稱為一階微分濾波器，將

?

　　稱為二階微分濾波器，平滑濾波器應滿足以下條件：

　　(1)當

?

　　為偶函數;

　　(2)

?

　　;

　　(3)h(x)一階及兩階可微;

　　上述第二個條件保證了信號經平滑濾波器h(x)濾波之后，其均值不變。

　　Marr 提出用下述的高斯函數作為平滑濾波器：

?