发明名称 一种基于线性规划模型的生物特征图像分析与识别方法
摘要 本发明公开了一种基于线性规划模型的生物特征图像识别方法,该方法包括以下步骤:生物特征训练样本图像归一化;得到基于局部子区域特征的生物特征图像表达;若特征个数过大,基于并行随机采样处理和线性规划模型选择最优特征,若否基于线性规划模型选择最优特征;将最优特征按权重大小降序排列,选择排名靠前的特征通过SVM得到分类器模型;待测试生物特征图像归一化;提取其最优局部子区域特征并输入分类器模型中,得到该生物特征图像的识别结果。本发明得到的生物特征图像特征鲁棒性好,识别精度高,因此本发明方法效率高,能够提高大规模生物特征数据比对的速度和精度,可用于生物特征识别的身份认证系统和其他需要安全性防范的诸应用系统中。
申请公布号 CN102902980B 申请公布日期 2015.12.02
申请号 CN201210339394.4 申请日期 2012.09.13
申请人 中国科学院自动化研究所 发明人 谭铁牛;孙哲南;王立彬
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 宋焰琴
主权项 一种基于线性规划模型的生物特征图像识别方法,其特征在于,该方法包括以下步骤:步骤S1,对生物特征训练样本图像进行归一化处理;步骤S2,对于归一化后的生物特征训练样本图像,提取其局部子区域特征,得到基于局部子区域特征的生物特征图像表达;步骤S3,判断所述生物特征图像表达中所包含的局部子区域特征个数是否过大,如是则转向步骤S5,如否则进入步骤S4;步骤S4,基于线性规划模型从局部子区域特征中选择最优特征;步骤S5,基于并行随机采样处理和线性规划模型从局部子区域特征中选择最优特征;步骤S6,将所述步骤S4或者S5得到的最优特征按其权重的取值大小进行降序排列,选择排名靠前的一定数量的特征,通过支持向量机SVM方法训练学习得到分类器模型;步骤S7,输入待测试生物特征图像;步骤S8,按照所述步骤S1对该待测试生物特征图像进行归一化处理;步骤S9,对于归一化后的待测试生物特征图像,提取其所述步骤S4或S5得到的最优局部子区域特征;步骤S10,将提取出的最优局部子区域特征输入到所述步骤S6得到的分类器模型中,从而得到对于该生物特征图像的识别结果;所述步骤S4进一步包括以下子步骤:步骤S41,建立具有最大间隔性质的损失函数模型;步骤S42,建立对所述损失函数模型中的特征权重的非负约束,即第i维特征的权重w<sub>i</sub>≥0i=1...D,其中,D为特征的总维数;步骤S43,计算经过所述步骤S42约束的每个局部子区域特征的可区分性参数信息f;步骤S44,对于每个样本的损失函数模型,加入相应非负的松弛变量ξ=(ξ<sub>1</sub>,ξ<sub>2</sub>,...ξ<sub>N</sub>)的限制,其中,N为样本数目;步骤S45,基于所述步骤S43得到的可区分性参数信息f和所述步骤S44得到的损失函数模型,建立标准线性规划模型,并采用单纯形算法求解出非负的特征权重向量,从而得到对于分类识别起到关键作用的特征;所述具有最大间隔性质的损失函数模型为Hinge函数:Loss=max(0,1‑yf(x)),其中,y代表类标记,f(x)代表分类器函数;上式等价于下面的线性不等式约束表达式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>+</mo></msubsup><mo>&le;</mo><mi>&alpha;</mi></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>+</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>-</mo></msubsup><mo>&GreaterEqual;</mo><mi>&beta;</mi></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>-</mo></msup></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0000763938040000021.GIF" wi="495" he="154" /></maths>其中,<img file="FDA0000763938040000023.GIF" wi="58" he="82" />和<img file="FDA0000763938040000024.GIF" wi="63" he="74" />分别表示第j个正、负样本,N<sup>+</sup>和N<sup>‑</sup>分别表示正负样本的个数,α和β分别为正负样本特征均值,w为特征权重向量;所述标准线性规划模型可表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>f</mi><mo>&CenterDot;</mo><mi>w</mi><mo>+</mo><mi>&lambda;</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&xi;</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000763938040000025.GIF" wi="369" he="137" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>+</mo></msubsup><mo>&le;</mo><mi>&alpha;</mi><mo>+</mo><msub><mi>&xi;</mi><mi>j</mi></msub></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>+</mo></msup><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>-</mo></msubsup><mo>&GreaterEqual;</mo><mi>&beta;</mi><mo>-</mo><msub><mi>&xi;</mi><mi>j</mi></msub></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>-</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>&xi;</mi><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd><mtd><mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000763938040000026.GIF" wi="662" he="245" /></maths>其中,<img file="FDA0000763938040000027.GIF" wi="362" he="143" />为目标函数,表示对特征权重向量作l<sub>1</sub>范数优化;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>+</mo></msubsup><mo>&le;</mo><mi>&alpha;</mi><mo>+</mo><msub><mi>&xi;</mi><mi>j</mi></msub></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>+</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mo>&CenterDot;</mo><msubsup><mi>x</mi><mi>j</mi><mo>-</mo></msubsup><mo>&GreaterEqual;</mo><mi>&beta;</mi><mo>-</mo><msub><mi>&xi;</mi><mi>j</mi></msub></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>...</mo><msup><mi>N</mi><mo>-</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>w</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>&xi;</mi><mo>&GreaterEqual;</mo><mn>0</mn></mrow></mtd><mtd><mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000763938040000028.GIF" wi="624" he="245" /></maths>为对于该目标函数的约束;w代表特征权重,为待优化的目标变量;f为先验的可区分性参数信息;<img file="FDA0000763938040000029.GIF" wi="71" he="98" /><img file="FDA00007639380400000210.GIF" wi="50" he="74" />分别代表正、负样本的特征;N<sup>+</sup>,N<sup>‑</sup>分别代表正、负样本的个数;α,β分别为训练集合中两类样本匹配分数的平均值;λ表示待定的参数。
地址 100190 北京市海淀区中关村东路95号