C++指纹识别系统论文

11 绪 纹自动识别技术的历史、现状和未来自动指纹识别技术是一种有着广泛应用背景的身份鉴定技术。本论文主要介绍了自动指纹识别系统的指纹图像处理和指纹匹配两部分。指纹识别的基本原理为：采用细节点坐标模型来做细节匹配，即对指纹的脊线末梢和指纹的脊线分支点提取特征点来鉴定指纹。具体的自动指纹识别系统框图见图11元前7000年到6000 年以前，古叙利亚和中国，指纹作为身份鉴别己经开始应用。考古发现，在这个时代，一些粘土陶器上留有陶艺匠人的指纹，中国的一些文件上印有起草者的大拇指指纹，在 古城市的房屋留有砖匠一对大拇指指纹的印记等。虽然指纹的一些特征己经被人们认识和接受而不能证明，但指纹己广泛应用社会的各个方面。19世纪初，科学研究发现了至今仍然承认的指纹的两个重要特征：一是两个不同手指的指纹脊线的式样(同，另外一个是指纹脊线的式样终生不改变。这个研究成果使得指纹在犯罪鉴别中得以正式应用。(主要代表性的事件有：1896年阿根廷首次应用，然后是1901年的苏格兰，20世纪初其他国家也相继应用到犯罪鉴别中)20世纪 60年代，由于计算机可以有效的处理图形，人们开始着手研究利用计算机来处理指纹。从那时起，自动指纹识别系统 0世纪80年代，个人电脑、光学扫描这两项技术的革新，使得它们作为指纹取像的工具成为现实，从而使指纹识别可以在其他领域中得以应用，比如代替在(90年代后期)，低价位取像设备的引入及其飞速发展，可靠的比对算法的发现为个人身份2识别应用的增长提供了舞台。对指纹识别技术来说，被广泛接受意味着在影响千万人的日常生活的各个地方使用。通过取代个人识别码和口令，指纹识别技术可以阻止非授权的访问，可以防止盗用窝电话、智能卡、桌面 作站及其计算机网络；在通过电话、网络进行的金融交易时进行身份认证；在建筑物或工作场所，指纹识别技术可以取代钥匙、证件、图章和卡阅读器。然而，必须指出的是，指纹识别技术到目前在国内乃至世界上仍未普及，其主要原因是：(1)速度慢，目前一般系统辨别时需要数秒钟甚至更长；(2)价格高，系统过于昂贵；(3)性能差，错误拒识率（错误接受率（高。因此，对快速而精确的指纹识别方法的研究有着重要的科学意义和应用价值。一般来说，一个指纹自动识别系统(要包括指纹图像提取子系统、指纹识别子系统和指纹压缩存储子系统几个部分组成。其中，指纹识别子系统是整个系统的核心部分，包括指纹图象处理、特征提取、指纹分类及匹配等四部分。本文将主要对指纹图象处理和特征提取的部分算法加以讨论和研究。文的主要研究工作本文主要的研究内容是在环境下实现指纹图像的小波变换和预处理，同时还研究指纹特征提取及匹配算法问题。全文对小波的基本理论进行了简要介绍，并对指纹图像小波变换的实现、基于小波变换频率特性的应用指纹图像滤波处理、指纹图像增强处理的理论以及编程实现进行了比较深入地研究，并列出了实现程序。此外，本文还在对指纹图像的预处理方法进行分析和研究的基础上，提出了一套指纹预处理算法，采用该算法可以得到清晰的指纹点线图，效果较为突出。3第2章 波分析图像处理小波变换是最近20多年来发展起来的用于信号分析和信号处理的一种新的域变换技术。由于小波变换是把信号在不同尺度上进行小波展开，它更适合于处理突变信号和非平稳信号。小波变换技术已受到人们极大的重视。一个信号可表示成一系列正弦和余弦函数之和。但就意味我们可以确定信号中包含的所有频率，但不能确定具有这些频率的信号出现在什么时候，只适宜处理平稳信号。在非平稳信号的分析中，人们希望存在一种变换函数，能够满足在高频信号中，有相对小的时间间隔以便给出较高的精度，而在低频信号中能够以相对较宽的时间间隔给出完全的信息。小波是有限宽度的基函数，这些基函数不仅在频率上而且在位置上是变换的，它们是有限宽度的波。基于它们的变换称为小波变换。小波变换具有时间一频率自动伸缩能力，这种能力可以在任何希望的频率范围上产生频谱信息。小波理论的提出可追溯到1910年975年981年明了小波函数的存在性。1984年法国地理学家除了人们对于小波函数是否存在的疑虑，从而真正掀起了小波研究的热潮。同年，将小波函数的构造统一于多分辨率分析的框架之下，同时，988年989年，990年，崔锦泰和王建忠构造了基于样条函数的单正交小波函数。1992年，此，小波理论系统的构架得以建立，它为原来信号处理领域里各自独立的方法，如多尺度分析、拉普拉斯金字塔、精确重建滤波器组等提供了一个统一的理论框架。随着理论研究的不断深入和应用领域的不断扩展，小波分析越来越显示出它的独特魅力。41、小波概念小波变换的定义是把某一被称为基本小波(也叫母小波函数 t做位移后，再在不同尺度下与待分析的信号 x t 做内积: 1, , 0z x t (效的频域表示是: *, 2 e d ( 中 X 和 分别是 x t 和 t 的傅里叶变换。小波变换有以下特点:1) 有多分辨率(也叫多尺度(特点，可以由粗及细地逐步观察信号。2) 可以看成用基本频率特性为 的带通滤波器在不同尺度下对信号做滤波。由于傅里叶变换的尺度特性可知，这组滤波器具有品质因数恒定，即相对带宽(带宽与中心频率之比)恒定的特点。注意，越大相对频率越低。3) 适当地选择基小波，使 t 在时域上为有限支撑， 在频域上也比较集中，就可以使在时、频域都具有表征信号局部特征的能力，因此有利于检测信号的瞬态或奇异点。正如上所述，小波分析的一个主要优点就是能够分析信号的局部特征。比如说，采用小波分析可以发现叠加在一个非常规范的正弦信号上的一个非常小的畸变信号的出现时间。传统的傅里叶变换只能得到平坦的频谱上的两个尖峰。利用小波分析可以非常准确地分析出信号在什么时刻发生畸变。小波分析可以检测出许多其他分析方法忽略的信号特性，例如，信号的趋势、信号的高阶不连续点、自相似特性。小波分析还能以非常小的失真度实现对信号的压缩与消噪，它在图像数据压缩方面的潜力已经得到确认。在二维情况下，小波分析除了“显微”能力外还具有“极化”能力(即方向选择性)，因而引人注意。2、小波在图像处理中的应用图像处理是小波分析应用的重要领域，近年来小波分析已被证明是进行图像处理强有力的工具之一，由于小波分析技术可以将信号或图像分层次按小波基展开，并且可以根据图像的性质及事先给定的图像处理要求确定到底要展开到哪一级为止，从而不仅能5有效地控制计算量，满足实时处理的需要，而且可以方便地实现通常由子带编码技术实现的累进编码(即采取逐步浮现的方式传送多媒体图像)。同时，小波变换具有放大、缩小和平移的功能，能够很方便地产生各种分辨率的图像，从而适合于不同分辨率图像的处理。因此，近十年来，基于小波变换的图像压缩算法得到了很大发展，取得了许多重要的成果，而且越来越成为该领域研究和应用的热点。续小波变换1、定义对任意信号 2( ) ( )f t L R ，f(t)的连续小波变换定义是: ,1, , 0a bR t bW f a b f t dt (2.1)果 t 满足容许条件; 2( )R d ，则称 t 为容许小波，可由 ( , )W f a b 恢复原信号1 , 2( ) ( , ) ( )a bR t C W f a b t a (中 2( )RC d ，根据 21, , , ,2f g f g f g L R 12, 2a e f a dW f a b (见f(t)的连续小波变换就相当于f(t)通过传递函数为 ( )a 的系统的输出 。2、特性在介绍小波特性之前先给出窗函数的概念。（1）定义:非平凡函数 )， )x x 也是属于)的，一个窗函数的的中心 ，定义分别是:6 22 221 Rt x x ( 122 22221 R x t x (定与都是窗函数，窗函数的中心与半径分别用 给出，那么函数,a b 是中心在b+半径等于a 的一个窗函数，因此由(可知连续小波积分给出了信号f(t)具有时间窗a ,b+a 的局部信息;窗函数的中心与半径分别用w*与 表示，则由(可知除了一个倍数 12 /2a a 与一个线性相位移。 外，连续小波变换还给出了 * */ / , / /a a a a ,a b 对应一系列带通系统，并且可以得出以下结论:（1）信号f(t)的连续小波变换就是一系列带通滤波器对f(t)滤波后的输出，2）通滤波器的带宽和中心频率也变化。心频率变大，带宽变宽:心频率变小，带宽变窄f(t)通过带通滤波器滤波，对分析信号的局部特性很有价值。信号变化缓慢的地方，主要为低频成份，频率范围也较窄，此时小波变换带通滤波器应相当于之，信号发生突变的地方，主要是高频成份，频率范围也较宽，小波变换的带通滤波器相当于伸缩因子波的范围也从低频到高频变化。散小波变换将连续小波变换的尺度时间位移得到离散小波变换。通常： 2 ja j z ,2 。所以，信号f(t)的离散小波变换(义为:,( , ) ( ) ( )j f j k f t t (, 0( ) 2 (2 )j jj k t t (中 , ( )j k t 取离散正交小波基， ( , )W f j k 是尺度t)的离散小波变换，离散小波变换具有以下的特性;71、 , ( )j k t 是小波函数以 t 在尺度上的伸缩和时域上的平移得到的。随 ( )j k t 在频域上处于不同的频段，随, ( )j k t 在时域上处于不同的时段，所以离散小波变换是一种信号的时间一频率分析;2 、尺度, ( )j k t 在时域上伸展，在频域上收缩，中心频率降低变换的时域分辨率降低，频域分辨率提高;尺度, ( )j k t 在时域上收缩，在频域上伸展，中心频率升高，变换的时域分辨率提高，分辨率的时频域分析。频域分辨率降低。所以，离散小波变换是一种多分辨率的时频域分析。纹图像二维小波变换概念及算法1、概念：一维信号的离散小波变换很容易推广到二维的情况。假设 ( )x 是一个一维的尺度函数， ( )y 是相应的小波函数，那么，可以得到一个二维小波变换的基础函数： 1 ,x y x y ， 2 ,x y x y ， 3 ,x y x y 2、算法图像可以看作是二维的矩阵，一般假设图像矩阵的大小为NN，且有N2n（那么每次小波变换后，图像便分解为4个大小为原来尺寸14的子块频带区域，如图2别包含了相应频带的小波系数，相当于在水平方向和竖直方向上进行隔点采样。进行下一层小波变换时，变换数据集中在2式（等式（明了图像小波变换的数学原型。频带保持了原始图像内容信息，图像的能量集中于此频带：102 2( , ) ( . ), ( 2 , 2 )j jf m n f x y x m y n (频带保持了图像水平方向上的高频边缘信息：11 12 2( , ) ( . ), ( 2 , 2 )j jf m n f x y x m y n （频带保持了图像竖直方向上的高频边缘信息：12 22 2( , ) ( . ), ( 2 , 2 )j jf m n f x y x m y n （频带保持了图像在对角线方向上的高频信息：13 32 2( , ) ( . ), ( 2 , 2 )j jf m n f x y x m y n (中表示内积运算。8图2次离散小波变换后的频率分布 图2维小波变换编程实现及运行结果下面将详细地按编程步骤对图像小波变换地实现进行分析。（1）函数描述）函数参数*维指针，指向原始的图像数据*波变换系数，存放一次水平变换后的小波系数*波变换系数，存放一次数值变换后的小波系数像属性参数，数值为原始图像的高度值像属性参数，数值为原始图像高度值的一半像属性参数，数值为原始图像的宽度值像属性参数，数值为原始图像宽度值的一半波变换的层数，数值为1层波变换因子，）函数代码 ：横坐标9：纵坐标; /图像矩阵的纵坐标;,第一次（时）时; ); /逻辑非操作后数值加1; /系数预测 /完成一个偶系数的边界处理(1);()2);0=0+，整数到整数的变换; ); /水平方向的变换结束/竖直方向的变换开始，数据源未水平变换后的小波系数;);(1);(1);(1); /边界处理(2);(2);/边界处理(; );(2);/存放小波系数，高频频带的系数则削弱其幅值; +3)/4*4;4if(if(,y,*i=0;,1, /允许图像复原操作标志150;0;y=0;yx)yx;yx=yx; y=0;y1; 1;在先给菜单选项定义消息图22息完成小波系数的显示。 /读取数字图像的文件头，获取图像的属性参数4);18+3)/4*4;if(if(,y,*;,3, /允许图像复员操作标志0;0; /得到小波系数的极大值和极小值y=0;yx)yx;yx=yx; /计算出小波系数的极值差(y=0; ,i=0;+3)/4*4;if(if(;,x,y,*;,i=0;+3)/4*4;if(if(;,x,y,*29;yx)yx;/yx=yx; =y=0;, /小波系数的正则化处理/正则化处理后，将其存放在设定显示设备的颜色灰度范围是0255y=0;+3)/4*4;if(if(,y,，存放极小值，存放极大值; /原始图像数据存储，数据用于图像对比度增强*46i=0;(y=0;);/找到: /判断是否是8里为了方便，只处理8它的可以类推）48:=256) /提示用户目前只支持256色位图的平滑！,系统提示, /解除锁定:); /返回 /创建对话框（为平均模板）=1.0;=1.0;=1.0;=0.0;=0.0;=1.0;=1.0;=1.0;=0.0;=0.0;0=1.0;1=1.0;2=1.0;3=0.0;4=0.0;5=0.0;6=0.0;7=0.0;8=0.0;9=0.0;0=0.0;1=9