发明名称 虹膜纹理归一化处理方法
摘要 本发明为一种虹膜纹理归一化处理方法,主要包括采用Log算子对采样的虹膜图像进行边缘提取,采用Hough变换,定位虹膜内外圆的圆心和半径,并根据其坐标偏差情况,采用不同的校正模型分别对生理性偏差和采样偏差进行纠正;根据瞳孔缩放情况,对不均匀缩放进行校正。最后将采样纹理点展开成矩形虹膜纹理图像。本发明方法大大提高了虹膜图像的可靠性,也降低了对图像采集的要求,为图像的自动采集提供了可能,并进一步为虹膜识别技术实际应用创造了条件。
申请公布号 CN1202490C 申请公布日期 2005.05.18
申请号 CN03115895.1 申请日期 2003.03.19
申请人 上海交通大学 发明人 施鹏飞;宫雅卓;邢磊
分类号 G06K9/00;G06T3/00 主分类号 G06K9/00
代理机构 上海正旦专利代理有限公司 代理人 陆飞
主权项 1、一种眼睛虹膜纹理归一化处理方法,其特征在于具体步骤如下:(1)采用高斯-拉普拉斯变换对虹膜边缘进行提取,以平滑掉虹膜纹理图像采集中造成的噪声;(2)采用Hough变换,来定位虹膜内外圆的圆心与半径;(3)计算虹膜内外圆圆心的坐标偏差,根据该偏差的范围确定属于生理性偏差或采样偏差:当虹膜内外圆心偏差d>阈值do时,判定其为生理性偏差;反之,判定其为采样偏差;其中,阈值do取3-5个像素;(4)对圆心偏差进行纠正:采用不同的校正模型分别对生理性偏差和采样偏差进行纠正;其中,对于生理性偏差,采用如下校正模型进行归一化处理:<math> <mrow> <mi>x</mi> <mo>=</mo> <mo>{</mo> <mi>r</mi> <mo>&times;</mo> <mo>[</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>&times;</mo> <mfrac> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mi>sin</mi> <mi>c</mi> </mrow> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>]</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>}</mo> <mo>&times;</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> </mrow> </math> <math> <mrow> <mi>y</mi> <mo>=</mo> <mo>{</mo> <mi>r</mi> <mo>&times;</mo> <mo>[</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>d</mi> <mo>)</mo> </mrow> <mo>&times;</mo> <mfrac> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mi>sin</mi> <mi>c</mi> </mrow> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>]</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>}</mo> <mo>&times;</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math> 其中θ∈(0,2π),r∈[0,1],R1为虹膜内圆半径,R2为虹膜外圆半径,C=∠BO1O2,d=∠O2BO1,X1、y1为虹膜内圆圆心坐标点,x、y为虹膜纹理放射线上笛卞尔坐标系下的任一点坐标,r、θ为该点在极坐标系下的坐标;对于采样偏差,采用如下校正模型进行归一化处理:x=r×AB×cosa+x1 y=r×AB×sina+y1 (9)其中<math> <mrow> <mi>AB</mi> <mo>=</mo> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>A</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>B</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>A</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>B</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math> <math> <mrow> <mo>=</mo> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>&times;</mo> <mi>cos</mi> <mi>b</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>cos</mi> <mi>b</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>&times;</mo> <mi>sin</mi> <mi>b</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>&times;</mo> <mi>sin</mi> <mi>b</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> </math> <math> <mrow> <mi>a</mi> <mo>=</mo> <mi>arctg</mi> <mfrac> <mrow> <msub> <mi>y</mi> <mi>B</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>A</mi> </msub> </mrow> <mrow> <msub> <mi>x</mi> <mi>B</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>A</mi> </msub> </mrow> </mfrac> </mrow> </math> x1、y1为虹膜内圆圆心O1的坐标点,R1为内圆半径,x2、y2为虹膜外圆圆心O2坐标点,R2为外圆半径,r∈[0,1],A、B为径向采样线与虹膜的内、外圆的交点,b∈(0,2π),b为AO1与水平线的夹角;(5)计算虹膜内外边缘半径比,并以此作为眼睛瞳孔缩放程度的度量;(6)采用校正函数,以虹膜内外边缘半径比为参数,进行径向的非线性采样,对瞳孔不均匀缩放进行校正;进行校正所采用的校正函数为<math> <mrow> <msub> <mi>r</mi> <mi>i</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mo>+</mo> <mo>&infin;</mo> </mrow> </munderover> <msub> <mi>a</mi> <mi>j</mi> </msub> <mo>&times;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mi>j</mi> </msubsup> <mo>)</mo> </mrow> <mo>&times;</mo> <mi>r</mi> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>&times;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>&times;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>a</mi> <mn>3</mn> </msub> <mo>&times;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mn>3</mn> </msubsup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <msub> <mrow> <mo>+</mo> <mi>a</mi> </mrow> <mi>n</mi> </msub> <mo>&times;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mi>n</mi> </msubsup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>)</mo> </mrow> <mo>&times;</mo> <mi>r</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math> 于是,每幅图像经畸变校正后的标记点半径Ri表示为Ri=F(Si)×r×(R1-R0)+R0 (11)其中,R0、R1为分别由Hough变换求得的序列中每幅虹膜图像的瞳孔半径和虹膜半径,Si为每幅图像对应的R0/R1比值,r为{Si}中最小值所对应图像的标记点的归一化半径,r∈[0,1];(7)将采样纹理点展开成矩形的虹膜纹理图像。
地址 200030上海市华山路1954号