CN101303729B - A New Fingerprint Singularity Detection Method - Google Patents
- ️Wed Jun 02 2010
CN101303729B - A New Fingerprint Singularity Detection Method - Google Patents
A New Fingerprint Singularity Detection Method Download PDFInfo
-
Publication number
- CN101303729B CN101303729B CN200810138119XA CN200810138119A CN101303729B CN 101303729 B CN101303729 B CN 101303729B CN 200810138119X A CN200810138119X A CN 200810138119XA CN 200810138119 A CN200810138119 A CN 200810138119A CN 101303729 B CN101303729 B CN 101303729B Authority
- CN
- China Prior art keywords
- mrow
- msub
- point
- msup
- munder Prior art date
- 2008-07-01 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 43
- 230000035772 mutation Effects 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000012847 principal component analysis method Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims 1
- 238000000605 extraction Methods 0.000 abstract description 10
- 238000005516 engineering process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Landscapes
- Collating Specific Patterns (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种新的指纹奇异点检测方法。它解决了目前奇异点提取方法严重依赖指纹方向场、不能有效处理低质量指纹图像、难以适于工程应用等问题，该方法分为前后两个阶段，前一阶段对现有的Poincare index庞加莱指数方法进行了有效的改进，并利用其提取指纹图像中的候选奇异点，后一阶段有效利用Gaussian-Hermite高斯-埃尔米特矩对候选奇异点进行去伪处理，该方法有效结合了奇异点周围邻域的纹线方向变化信息和奇异点周围局部区域的纹线一致性变化趋势信息，能准确可靠地提取指纹图像中的奇异点，具有抗噪性强、准确可靠性高、工程应用价值大等特点。
The invention discloses a new fingerprint singular point detection method. It solves the problems that the current singular point extraction method relies heavily on the fingerprint direction field, cannot effectively deal with low-quality fingerprint images, and is difficult to be suitable for engineering applications. The method is divided into two stages. The Lay index method has been effectively improved and used to extract the candidate singular points in the fingerprint image. In the latter stage, the Gaussian-Hermite Gauss-Hermite moment is effectively used to remove the false processing of the candidate singular points. This method effectively combines the The ridge direction change information in the neighborhood around the singular point and the ridge consistency change trend information in the local area around the singular point can accurately and reliably extract the singular point in the fingerprint image, which has strong noise resistance, high accuracy and reliability, and engineering Features such as great application value.
Description
技术领域technical field
本发明涉及一种指纹图像检测方法，具体地说一种实用自动指纹识别系统(AFIS)中的低质量指纹图像的新的指纹奇异点检测方法。The invention relates to a fingerprint image detection method, in particular to a novel fingerprint singular point detection method for low-quality fingerprint images in an automatic fingerprint identification system (AFIS).
背景技术Background technique
目前在实用自动指纹识别系统中，指纹分类技术是加快系统识别速度的关键技术之一，而现今主流的分类技术多是依据奇异点的数目、类型和位置等信息来实现的，而且在处理低质量指纹图像时所采用的纹理匹配算法也需要准确可靠的奇异点信息。现今自动指纹识别系统中采用的主流的奇异点提取方法，绝大多数依赖于指纹方向场的准确提取，但在处理低质量指纹图像时，由于可靠的计算纹线方向本身就是一个难题，因而这些方法提取的奇异点不仅定位不够准确，在纹线方向计算有误的地方以及一些噪声污染的地方，还往往容易检测到许多虚假的奇异点。这使得这些方法难以有效满足工程应用。在实际的应用中需要一种准确可靠的奇异点提取算法。At present, in the practical automatic fingerprint identification system, the fingerprint classification technology is one of the key technologies to speed up the identification speed of the system, but the current mainstream classification technology is mostly realized based on the number, type and position of singular points, and is low in processing. The texture matching algorithm used in quality fingerprint images also requires accurate and reliable singular point information. Most of the mainstream singular point extraction methods used in today's automatic fingerprint identification systems rely on the accurate extraction of the fingerprint direction field. The singular points extracted by the method are not only inaccurately positioned, but also often detect many false singular points easily in places where the calculation of the ridge direction is wrong and where some noise is polluted. This makes these methods difficult to effectively meet engineering applications. In practical applications, an accurate and reliable singular point extraction algorithm is needed.
发明内容Contents of the invention
本发明的目的就是为了解决现有的奇异点提取方法严重依赖指纹方向场、不能有效处理低质量指纹图像、难以适于工程应用等问题，提供一种准确可靠的新的指纹奇异点检测方法，该方法分为前后两个阶段，前一阶段对现有的Poincare index庞加莱指数方法进行了有效的改进，并利用其提取指纹图像中的候选奇异点，后一阶段效利用Gaussian-Hermite高斯-埃尔米特矩矩对候选奇异点进行去伪处理，该方法有效结合了奇异点周围邻域的纹线方向变化信息和奇异点周围局部区域的纹线一致性变化趋势信息，能准确可靠地提取指纹图像中的奇异点，具有抗噪性强、准确可靠性高、工程应用价值大等特点。The purpose of the present invention is to provide a new accurate and reliable fingerprint singular point detection method to solve the problems that the existing singular point extraction method relies heavily on the fingerprint direction field, cannot effectively process low-quality fingerprint images, and is difficult to be suitable for engineering applications. The method is divided into two stages. In the former stage, the existing Poincare index method is effectively improved and used to extract candidate singular points in the fingerprint image. In the latter stage, the Gaussian-Hermite Gaussian -The Hermitian moment de-aliases the candidate singular points. This method effectively combines the ridge direction change information of the neighborhood around the singular point and the ridge consistency change trend information of the local area around the singular point, which can be accurate and reliable. It has the characteristics of strong anti-noise, high accuracy and reliability, and great engineering application value.
为实现上述目的，本发明采用如下技术方案：To achieve the above object, the present invention adopts the following technical solutions:
一种准确可靠的奇异点提取方法，它的方法为，An accurate and reliable singular point extraction method, its method is,
(1)预备阶段：包括指纹图像的背景分离和方向场信息计算。(1) Preparatory stage: including the background separation of the fingerprint image and the calculation of the direction field information.
(2)第一阶段：候选奇异点的提取，包括以下几个步骤：(2) The first stage: the extraction of candidate singular points, including the following steps:
(a)方向域O′中的每一个点(i，j)按公式(1)计算其相应的Poincare index庞加莱指数PG，C(i，j)。本算法中采用了两条长度分别为5×5和7×7的封闭曲线来计算每个点的Poincare index庞加莱指数值，只要其中一条的计算结果符合奇异点条件，则认定该点为候选奇异点，若两条曲线检测出的奇异点类型不一致，则直接认定该点为伪点。(a) Each point (i, j) in the direction domain O' calculates its corresponding Poincare index PG, C (i, j) according to the formula (1). In this algorithm, two closed curves with lengths of 5×5 and 7×7 are used to calculate the Poincare index value of each point. As long as the calculation result of one of them meets the singular point condition, the point is considered to be Candidate singular points, if the types of singular points detected by the two curves are inconsistent, the point is directly identified as a false point.
(b)对上一步检测到的core核心点和delta三角点运用聚类算法按欧氏距离分别进行聚类，并统计各个聚类包含的奇异点的个数Nm。Nm＜N，(N表示指纹模式区应该能够检测到的奇异点的个数，这里N取为25)的聚类被删除。(b) Use the clustering algorithm to cluster the core points and delta triangle points detected in the previous step according to the Euclidean distance, and count the number N m of singular points contained in each cluster. N m <N, (N represents the number of singular points that should be detectable in the fingerprint pattern area, where N is taken as 25) clusters are deleted.
(c)求取剩余的各个聚类的平均突变程度，并按奇异点类型、突变程度大小分别对这些聚类进行排序，选取前M(本方法M＝3，core核心点点聚类和delta三角点点聚类都取3个，当实际聚类个数少于3个时取实际聚类个数)个突变程度较大的聚类作为候选奇异点聚类，候选奇异点最终就定位在这些聚类的质心上。(c) Find the average mutation degree of each remaining cluster, and sort these clusters according to the type of singular point and the degree of mutation, and select the top M (this method M=3, core core point clustering and delta triangle Take 3 point clusters, and when the actual number of clusters is less than 3, take the actual number of clusters) clusters with a greater degree of mutation as candidate singular point clusters, and the candidate singular points are finally located in these clusters on the centroid of the class.
(3)第二阶段：候选奇异点的去伪，包含以下几个步骤：(3) The second stage: de-falsification of candidate singular points, including the following steps:
(a)对每个候选奇异点，按公式(2)计算其周围半径为6τ的圆形邻域内每个像素的分布一致性coherence值(τ为平均纹线距离)。(a) For each candidate singular point, calculate the distribution consistency coherence value of each pixel in a circular neighborhood with a radius of 6τ (τ is the average ridge distance) according to formula (2).
(b)然后将该圆形区域划分为32个扇形和一个中心小圆形区域并计算各个部分的平均coherence值，圆形区域的半径和两个圆带的宽度均为2τ。(b) Then divide the circular area into 32 sectors and a central small circular area and calculate the average coherence value of each part. The radius of the circular area and the width of the two circular bands are both 2τ.
(c)最后，比较16个方向每个方向上的三个区域的平均分布一致性coherence值的大小(中心圆形区域为16个方向公共的区域)，若自里向外分布一致性coherence值越来越大，则标记该方向为有效方向。对某一候选奇异点而言，若有效方向数大于等于10，则提取该点为奇异点，否则，认为该点为伪点。若奇异点的某些方向出了指纹边界则将剩余的方向作为参考方向，在这些方向上计算有效方向。此外，经实验分析，计算每个象素分布一致性coherence值的窗口大小取为(4τ+1)×(4τ+1)较为合适。(c) Finally, compare the average coherence value of the three areas in each direction of the 16 directions (the central circular area is the area common to the 16 directions), if the coherence value is distributed from the inside to the outside If it becomes larger and larger, this direction is marked as a valid direction. For a candidate singular point, if the number of effective directions is greater than or equal to 10, the point is extracted as a singular point, otherwise, the point is considered as a false point. If some directions of the singular point are out of the fingerprint boundary, the remaining directions are taken as reference directions, and effective directions are calculated on these directions. In addition, through experimental analysis, it is more appropriate to take (4τ+1)×(4τ+1) as the window size for calculating the coherence value of each pixel distribution consistency.
第一阶段，所述步骤(a)中，计算每一个奇异点(i，j)的Poincare index庞加莱指数值时按如下的公式进行：In the first stage, in the step (a), the Poincare index Poincare index value of each singular point (i, j) is calculated according to the following formula:
PoincarePoincare (( ii .. jj )) == 11 22 ππ ΣΣ kk == 00 NψNψ -- 11 ΔΔ (( kk )) ,,
ΔΔ (( kk )) == δδ (( kk )) ,, ifif || δδ (( kk )) || << ππ 22 ππ ++ δδ (( kk )) ,, ifδifδ (( kk )) ≤≤ -- ππ 22 ππ -- δδ (( kk )) ,, otherwiseotherwise -- -- -- (( 11 ))
δ(k)＝O′(ψx(i′)，ψy(i′))-O′(ψx(i)，ψy(i))，δ(k)=O'(ψ x (i'), ψ y (i'))-O'(ψ x (i), ψ y (i)),
i′＝(i+1)mod Nψ，i'=(i+1)mod N ψ ,
其中ψx(i)和ψy(i)分别是以给定点为中心的具有Nψ个像素的封闭曲线上第k个点的x和y坐标。δ(k)表示两相邻方向角的差值，Δ(k)表示对差值调整以后的结果，i′表示第i个点之后的下一个点。where ψ x (i) and ψ y (i) are the x and y coordinates, respectively, of the kth point on a closed curve with Nψ pixels centered at a given point. δ(k) represents the difference between two adjacent direction angles, Δ(k) represents the result after adjusting the difference, and i' represents the next point after the i-th point.
本算法中采用了5×5和7×7的正方形封闭曲线，该公式的形成基于对原有Poincare Index庞加莱指数方法的改进，用于计算某点对应封闭曲线上方向角变化的累积。In this algorithm, 5×5 and 7×7 square closed curves are used. The formation of this formula is based on the improvement of the original Poincare Index method, which is used to calculate the accumulation of direction angle changes on the closed curve corresponding to a certain point.
所述步骤(a)中，奇异点条件是指本方法提出的奇异点应满足的条件，该条件是对公式(1)的计算过程进行约束，具体是指，在沿封闭曲线做方向角变化的累加时附加以下限制条件：In the step (a), the singular point condition refers to the condition that the singular point proposed by this method should satisfy, and this condition is to constrain the calculation process of the formula (1), specifically refers to the direction angle change along the closed curve The following restrictions are added to the accumulation of :
(1)统计这Nψ个方向的符号变化次数，若方向由正到负和由负到正各发生一次，且仅有一次，则继续计算其Poincare index庞加莱指数值，否则，认为该点是普通点。(1) Count the number of sign changes in these Nψ directions. If the direction changes from positive to negative and from negative to positive, and only once, continue to calculate its Poincare index value, otherwise, consider this point It's ordinary.
(2)统计这Nψ个差中绝对值 
(3)如果最终Poincare index庞加莱指数的值为1/2，那么该给定点(i，j)就被确定为core核心点，如果Poincare index庞加莱指数值为-1/2，那么该给定点(i，j)就被确定为delta三角点。(3) If the value of the final Poincare index Poincare index is 1/2, then the given point (i, j) is determined as the core core point, if the value of the Poincare index Poincare index is -1/2, then The given point (i, j) is determined as a delta triangle point.
第一阶段，所述步骤(b)中，N取值25是基于对试验结果的分析，在指纹的奇异点区能够检测到的奇异点数目是相对稳定的且集中在奇异点区的最内部，对于低质量指纹图像，只要奇异点区不受大噪声的污染，那么在奇异点区内能够检测到的奇异点的数目也是稳定的，本方法N取25。In the first stage, in the step (b), the N value of 25 is based on the analysis of the test results, and the number of singular points that can be detected in the singular point area of the fingerprint is relatively stable and concentrated in the innermost part of the singular point area , for low-quality fingerprint images, as long as the singular point area is not polluted by large noise, the number of singular points that can be detected in the singular point area is also stable, and N is set to 25 in this method.
第一阶段，所述步骤(c)中，突变程度是指，封闭曲线上相邻两方向角作差过程中，绝对值 
第二阶段，所述步骤(a)中，本方法利用分布一致性coherence值来表达纹线的一致性信息，而分布一致性coherence值的计算采用了Gaussian-Hermite高斯-埃尔米特矩。本方法用四个矩来描述指纹的纹线一致性信息：M0，1，M1，0，M0，3，M3，0。并对这四个不同阶的Gaussian-Hermite矩做如下定义：In the second stage, in the step (a), the method uses the coherence value of the distribution consistency to express the consistency information of the ridges, and the calculation of the coherence value of the distribution consistency adopts the Gaussian-Hermite Gaussian-Hermite moment. This method uses four moments to describe the ridge consistency information of the fingerprint: M 0,1 , M 1,0 , M 0,3 , M 3,0 . And the Gaussian-Hermite moments of these four different orders are defined as follows:
Mm uu (( xx ,, ythe y )) == λMλ M 1,01,0 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) ++ (( 11 -- λλ )) Mm 3,03,0 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) Mm vv (( xx ,, ythe y )) == λMλ M 0,10,1 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) ++ (( 11 -- λλ )) Mm 0.30.3 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) -- -- -- (( 22 ))
这里λ(0＜λ＜1)是不同阶的Gaussian-Hermite高斯-埃尔米特矩的联合权重系数(本方法中λ取值0.5)。利用定义(2)，指纹图像中的每一个象素(i，j)都可以计算得到一个特征向量[Mu，Mv]T。在指纹奇异点区和非奇异点区[Mu，Mv]T的分布具有不同的特征，在非奇异点区[Mu，Mv]T沿着垂直于脊线的方向分布，在奇异点区[Mu，Mv]T则均匀地分布在各个方向上，本方法采用主分量分析方法提取[Mu，Mv]T的分布特性，[Mu，Mv]T协方差矩阵Cm由下式定义：Here λ (0<λ<1) is the joint weight coefficient of Gaussian-Hermite moments of different orders (λ takes the value of 0.5 in this method). Using definition (2), each pixel (i, j) in the fingerprint image can be calculated to obtain a feature vector [M u , M v ] T . The distribution of [M u , M v ] T in the fingerprint singular point area and non-singular point area has different characteristics. In the non-singular point area [M u , M v ] T distributes along the direction perpendicular to the ridge line. The point area [M u , M v ] T is evenly distributed in all directions. This method uses the principal component analysis method to extract the distribution characteristics of [M u , M v ] T , and the covariance matrix of [M u , M v ] T Cm is defined by:
CC Mm == ΣΣ ww (( Mm uu -- mm uu )) 22 ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) ΣΣ ww (( Mm vv -- mm vv )) 22 -- -- -- (( 33 ))
其中， 
设λ1和λ2是协方差矩阵CM的两个特征值，则当λ1＞＞λ2时，[Mu，Mv]T的分布主要是沿着长轴方向分布，也即沿着垂直于脊线的方向分布，而在噪声区或者奇异点区λ1和λ2的值是很接近的，因此，定义[Mu，Mv]T的分布一致性coherence特征如下：Let λ1 and λ2 be the two eigenvalues of the covariance matrix C M , then when λ1>>λ2, the distribution of [M u , M v ] T is mainly distributed along the long axis, that is, along the direction perpendicular to the ridge The direction distribution of the line, and the values of λ1 and λ2 in the noise area or singular point area are very close. Therefore, the distribution consistency coherence characteristics of [M u , M v ] T are defined as follows:
cohdrencecohdrence == λλ 11 22 -- λλ 22 22 λλ 11 22 ++ λλ 22 22 -- -- -- (( 44 ))
== (( ΣΣ ww (( Mm uu -- mm uu )) 22 ++ ΣΣ ww (( Mm vv -- Mm vv )) 22 )) (( ΣΣ ww (( Mm uu -- mm uu )) 22 -- ΣΣ ww (( Mm vv -- mm vv )) 22 )) 22 ++ 44 (( ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) )) 22 (( ΣΣ ww (( Mm uu -- mm uu )) 22 )) 22 ++ (( ΣΣ ww (( Mm vv -- mm vv )) 22 )) 22 ++ 22 (( ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) )) 22
由此，纹线一致性越好，分布一致性coherence越大。Therefore, the better the consistency of the ridges, the greater the coherence of the distribution consistency.
本发明的有益效果：由于采用了两阶段的处理方法，奇异点提取的可靠性有了较大程度的提高，特别是在处理低质量指纹图像时，奇异点的漏检、误检现象有了较大程度的减少。此外，对原有的Poincare index庞加莱指数方法进行了改进，提高其抗噪能力，避免了对方向场的多次平滑处理，使最终提取的奇异点准确性得到有效提高。满足了实用自动指纹识别系统(AFIS)的应用需求。Beneficial effects of the present invention: due to the adoption of the two-stage processing method, the reliability of singular point extraction has been greatly improved, especially when processing low-quality fingerprint images, the phenomenon of missed detection and false detection of singular points has been eliminated. a greater degree of reduction. In addition, the original Poincare index method is improved to improve its anti-noise ability, avoid multiple smoothing of the direction field, and effectively improve the accuracy of the final extracted singular points. It meets the application requirements of practical automatic fingerprint identification system (AFIS).
附图说明Description of drawings
图1为去伪模板图像；Fig. 1 is to remove false template image;
图2为算法流程图。Figure 2 is the flow chart of the algorithm.
具体实施方式Detailed ways
下面结合附图与实施例对本发明做进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and embodiments.
一种准确可靠的奇异点提取方法，它的方法为，An accurate and reliable singular point extraction method, its method is,
(1)预备阶段：包括指纹图像的背景分离和方向场信息计算。(1) Preparatory stage: including the background separation of the fingerprint image and the calculation of the direction field information.
(2)第一阶段：候选奇异点的提取，包括以下几个步骤：(2) The first stage: the extraction of candidate singular points, including the following steps:
(a)方向域O′中的每一个点(i，j)按公式(1)计算其相应的Poincare index庞加莱指数PG，C(i，j)。本算法中采用了两条长度分别为5×5和7×7的封闭曲线来计算每个点的Poincareindex庞加莱指数值，只要其中一条的计算结果符合奇异点条件，则认定该点为候选奇异点，若两条曲线检测出的奇异点类型不一致，则直接认定该点为伪点。(a) Each point (i, j) in the direction domain O' calculates its corresponding Poincare index PG, C (i, j) according to the formula (1). In this algorithm, two closed curves with lengths of 5×5 and 7×7 are used to calculate the Poincare index value of each point. As long as the calculation result of one of them meets the singular point condition, the point is considered as a candidate Singular point, if the types of singular points detected by the two curves are inconsistent, the point is directly identified as a false point.
(b)对上一步检测到的core核心点和delta三角点运用聚类算法按欧氏距离分别进行聚类，并统计各个聚类包含的奇异点的个数Nm，Nm＜N(这里N取为25)的聚类被删除。(b) Use the clustering algorithm to cluster the core points and delta triangle points detected in the previous step according to the Euclidean distance, and count the number N m of singular points contained in each cluster, N m <N (here The clusters where N is taken as 25) are deleted.
(c)求取剩余的各个聚类的平均突变程度，并按奇异点类型、突变程度大小分别对这些聚类进行排序，选取前M(本方法M＝3，core点聚类和delta点聚类都取3个，当实际聚类个数少于3个时取实际聚类个数)个突变程度较大的聚类作为候选奇异点聚类，候选奇异点最终就定位在这些聚类的质心上。(c) Calculate the average mutation degree of the remaining clusters, and sort these clusters according to the type of singular point and the degree of mutation, and select the top M (this method M=3, core point clustering and delta point clustering If the number of clusters is less than 3, the actual number of clusters is taken as the number of clusters with a large mutation degree as candidate singular point clusters, and the candidate singular points are finally located in the clusters of these clusters. on the centroid.
(3)第二阶段：候选奇异点的去伪，如图2所示，包含以下几个步骤：(3) The second stage: the de-falsification of candidate singular points, as shown in Figure 2, includes the following steps:
(a)对每个候选奇异点，按公式(2)计算其周围半径为6τ的圆形邻域内每个像素的分布一致性coherence值(τ为平均纹线距离)。(a) For each candidate singular point, calculate the distribution consistency coherence value of each pixel in a circular neighborhood with a radius of 6τ (τ is the average ridge distance) according to formula (2).
(b)然后按图1所示的模板将该圆形区域划分为32个扇形和一个中心小圆形区域并计算各个部分的平均分布一致性coherence值，圆形区域的半径和两个圆带的宽度均为2τ。(b) Then divide the circular area into 32 sectors and a central small circular area according to the template shown in Figure 1, and calculate the average distribution consistency coherence value of each part, the radius of the circular area and two circular bands The width of is 2τ.
(c)最后，比较16个方向每个方向上的三个区域的平均分布一致性coherence值的大小(中心圆形区域为16个方向公共的区域)，若自里向外分布一致性coherence值越来越大，则标记该方向为有效方向。对某一候选奇异点而言，若有效方向数大于等于10，则提取该点为奇异点，否则，认为该点为伪点。若奇异点的某些方向出了指纹边界则将剩余的方向作为参考方向，在这些方向上计算有效方向。此外，经实验分析，计算每个象素分布一致性coherence值的窗口大小取为(4τ+1)×(4τ+1)较为合适。具体的算法流程图如图2所示。(c) Finally, compare the average coherence value of the three areas in each direction of the 16 directions (the central circular area is the area common to the 16 directions), if the coherence value is distributed from the inside to the outside If it becomes larger and larger, this direction is marked as a valid direction. For a candidate singular point, if the number of effective directions is greater than or equal to 10, the point is extracted as a singular point, otherwise, the point is considered as a false point. If some directions of the singular point are out of the fingerprint boundary, the remaining directions are taken as reference directions, and effective directions are calculated on these directions. In addition, through experimental analysis, it is more appropriate to take (4τ+1)×(4τ+1) as the window size for calculating the coherence value of each pixel distribution consistency. The specific algorithm flow chart is shown in Figure 2.
第一阶段，所述步骤(a)中，计算每一个点(i，j)的Poincare index值时按如下的公式进行：In the first stage, in the step (a), the Poincare index value of each point (i, j) is calculated according to the following formula:
PoincarePoincare (( ii .. jj )) == 11 22 ππ ΣΣ kk == 00 NψNψ -- 11 ΔΔ (( kk )) ,,
ΔΔ (( kk )) == δδ (( kk )) ,, ifif || δδ (( kk )) || << ππ 22 ππ ++ δδ (( kk )) ,, ifδifδ (( kk )) ≤≤ -- ππ 22 ππ -- δδ (( kk )) ,, otherwiseotherwise -- -- -- (( 11 ))
δ(k)＝O′(ψx(i′)，ψy(i′))-O′(ψx(i)，ψy(i))，δ(k)=O'(ψ x (i'), ψ y (i'))-O'(ψ x (i), ψ y (i)),
i′＝(i+1)mod Nψ，i'=(i+1)mod N ψ ,
其中ψx(i)和ψy(i)分别是以给定点为中心的具有Nψ个像素的封闭曲线上第k个点的x和y坐标。δ(k)表示两相邻方向角的差值，Δ(k)表示对差值调整以后的结果，i′表示第i个点之后的下一个点。where ψ x (i) and ψ y (i) are the x and y coordinates, respectively, of the kth point on a closed curve with Nψ pixels centered at a given point. δ(k) represents the difference between two adjacent direction angles, Δ(k) represents the result after adjusting the difference, and i' represents the next point after the i-th point.
本算法中采用了5×5和7×7的正方形封闭曲线，该公式的形成基于对原有Poincare Index庞加莱指数方法的改进，用于计算某点对应封闭曲线上方向角变化的累积。In this algorithm, 5×5 and 7×7 square closed curves are used. The formation of this formula is based on the improvement of the original Poincare Index method, which is used to calculate the accumulation of direction angle changes on the closed curve corresponding to a certain point.
所述步骤(a)中，奇异点条件是指本方法提出的奇异点应满足的条件，该条件是对公式(1)的计算过程进行约束，具体是指，在沿封闭曲线做方向角变化的累加时附加以下限制条件：In the step (a), the singular point condition refers to the condition that the singular point proposed by this method should satisfy, and this condition is to constrain the calculation process of the formula (1), specifically refers to the direction angle change along the closed curve The following restrictions are added to the accumulation of :
(1)统计这Nψ个方向的符号变化次数，若方向由正到负和由负到正各发生一次，且仅有一次，则继续计算其Poincare index庞加莱指数值，否则，认为该点是普通点。(1) Count the number of sign changes in these Nψ directions. If the direction changes from positive to negative and from negative to positive, and only once, continue to calculate its Poincare index value, otherwise, consider this point It's ordinary.
(2)统计这Nψ个差中绝对值 
(3)如果最终Poincare index庞加莱指数的值为1/2，那么该给定点(i，j)就被确定为core核心点，如果Poincare index庞加莱指数值为-1/2，那么该给定点(i，j)就被确定为delta三角点。(3) If the value of the final Poincare index Poincare index is 1/2, then the given point (i, j) is determined as the core core point, if the value of the Poincare index Poincare index is -1/2, then The given point (i, j) is determined as a delta triangle point.
第一阶段，所述步骤(b)中，N取值25是基于对试验结果的分析，在指纹的奇异点区能够检测到的奇异点数目是相对稳定的且集中在奇异点区的最内部，对于低质量指纹图像，只要奇异点区不受大噪声的污染，那么在奇异点区内能够检测到的奇异点的数目也是稳定的，本方法N取25。In the first stage, in the step (b), the N value of 25 is based on the analysis of the test results, and the number of singular points that can be detected in the singular point area of the fingerprint is relatively stable and concentrated in the innermost part of the singular point area , for low-quality fingerprint images, as long as the singular point area is not polluted by large noise, the number of singular points that can be detected in the singular point area is also stable, and N is set to 25 in this method.
第一阶段，所述步骤(c)中，突变程度是指，封闭曲线上相邻两方向角作差过程中，绝对值的那个差值，平均突变程度是指每一聚类内所有像素的突变程度的平均。In the first stage, in the step (c), the degree of mutation refers to the absolute value of the difference between two adjacent direction angles on the closed curve. The difference, the average mutation degree refers to the average of the mutation degree of all pixels in each cluster.
第二阶段，所述步骤(a)中，本方法利用分布一致性coherence值来表达纹线的一致性信息，而分布一致性coherence值的计算采用了Gaussian-Hermite矩。本方法用四个矩来描述指纹的纹线一致性信息：M0，1，M1，0，M0，3，M3，0。并对这四个不同阶的Gaussian-Hermite高斯-埃尔米特矩做如下定义：In the second stage, in the step (a), the method uses the coherence value of the distribution consistency to express the consistency information of the ridges, and the calculation of the coherence value of the distribution consistency adopts Gaussian-Hermite moments. This method uses four moments to describe the ridge consistency information of the fingerprint: M 0,1 , M 1,0 , M 0,3 , M 3,0 . And the Gaussian-Hermite Gaussian-Hermite moments of these four different orders are defined as follows:
Mm uu (( xx ,, ythe y )) == λMλ M 1,01,0 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) ++ (( 11 -- λλ )) Mm 3,03,0 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) Mm vv (( xx ,, ythe y )) == λMλ M 0,10,1 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) ++ (( 11 -- λλ )) Mm 0.30.3 (( xx ,, ythe y ,, II (( xx ,, ythe y )) )) -- -- -- (( 22 ))
这里λ(0＜λ＜1)是不同阶的Gaussian-Hermite矩的联合权重系数(本方法中λ取值0.5)。利用定义(2)，指纹图像中的每一个象素(i，j)都可以计算得到一个特征向量[Mu，Mv]T。在指纹奇异点区和非奇异点区[Mu，Mv]T的分布具有不同的特征，在非奇异点区[Mu，Mv]T沿着垂直于脊线的方向分布，在奇异点区[Mu，Mv]T则均匀地分布在各个方向上，本方法采用主分量分析方法提取[Mu，Mv]T的分布特性，[Mu，Mv]T协方差矩阵CM由下式定义：Here λ (0<λ<1) is the joint weight coefficient of Gaussian-Hermite moments of different orders (in this method, λ takes a value of 0.5). Using definition (2), each pixel (i, j) in the fingerprint image can be calculated to obtain a feature vector [M u , M v ] T . The distribution of [M u , M v ] T in the fingerprint singular point area and non-singular point area has different characteristics. In the non-singular point area [M u , M v ] T distributes along the direction perpendicular to the ridge line. The point area [M u , M v ] T is evenly distributed in all directions. This method uses the principal component analysis method to extract the distribution characteristics of [M u , M v ] T , and the covariance matrix of [M u , M v ] T C M is defined by:
CC Mm == ΣΣ ww (( Mm uu -- mm uu )) 22 ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) ΣΣ ww (( Mm vv -- mm vv )) 22 -- -- -- (( 33 ))
其中， 
设λ1和λ2是协方差矩阵CM的两个特征值，则当λ1＞＞λ2时，[Mu，Mv]T的分布主要是沿着长轴方向分布，也即沿着垂直于脊线的方向分布，而在噪声区或者奇异点区λ1和λ2的值是很接近的，因此，定义[Mu，Mv]T的分布一致性coherence特征如下：Let λ1 and λ2 be the two eigenvalues of the covariance matrix C M , then when λ1>>λ2, the distribution of [M u , M v ] T is mainly distributed along the long axis, that is, along the direction perpendicular to the ridge The direction distribution of the line, and the values of λ1 and λ2 in the noise area or singular point area are very close. Therefore, the distribution consistency coherence characteristics of [M u , M v ] T are defined as follows:
cohdrencecohdrence == λλ 11 22 -- λλ 22 22 λλ 11 22 ++ λλ 22 22 -- -- -- (( 44 ))
== (( ΣΣ ww (( Mm uu -- mm uu )) 22 ++ ΣΣ ww (( Mm vv -- Mm vv )) 22 )) (( ΣΣ ww (( Mm uu -- mm uu )) 22 -- ΣΣ ww (( Mm vv -- mm vv )) 22 )) 22 ++ 44 (( ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) )) 22 (( ΣΣ ww (( Mm uu -- mm uu )) 22 )) 22 ++ (( ΣΣ ww (( Mm vv -- mm vv )) 22 )) 22 ++ 22 (( ΣΣ ww (( Mm uu -- mm uu )) (( Mm vv -- mm vv )) )) 22
由此，纹线一致性越好，分布一致性coherence越大。Therefore, the better the consistency of the ridges, the greater the coherence of the distribution consistency.
Claims (3)
1. A new fingerprint singular point detection method is characterized in that the method comprises the following steps,
(1) a preparation stage: carrying out background separation and direction field information calculation on the fingerprint image;
(2) the first stage is as follows: the method for extracting the candidate singular points comprises the following steps:
(a) calculating for each point (i, j) in the orientation field O' its corresponding Poincare index value PG，C(i, j) calculating Poincare i at each point using two closed curves of 5 × 5 and 7 × 7 lengths, respectivelyDetermining a point as a candidate singular point as long as a calculation result of one of the ndex poincare index values meets the singular point condition, and directly determining the point as a pseudo point if the types of the singular points detected by the two curves are inconsistent; the calculation of the Poincare index value for each singular point (i, j) is performed according to the following formula:
<math><mrow><mi>Poincare</mi><mrow><mo>(</mo><mi>i</mi><mo>.</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>N</mi><mi>ψ</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><mi>Δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>
<math><mrow><mi>Δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi><mo>|</mo><mi>δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo><</mo><mfrac><mi>π</mi><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mi>π</mi><mo>+</mo><mi>δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mi>ifδ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>≤</mo><mo>-</mo><mfrac><mi>π</mi><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mi>π</mi><mo>-</mo><mi>δ</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mi>otherwise</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>
δ(k)＝O′(ψx(i′)，ψy(i′))-O′(ψx(i)，ψy(i))，
i′＝(i+1)modNψ，
wherein psix(i) And psiy(i) Each centered on a given point having NψX and y coordinates of a kth point on a closed curve of the pixels; δ (k) represents a difference between two adjacent direction angles, Δ (k) represents a result after adjustment of the difference, and i' represents a point next to the i-th point; the algorithm adopts 5 multiplied by 5 and 7 multiplied by 7 square closed curves, and the formula is formed based on the improvement of the original Poincare index method and is used for calculating the accumulation of the change of the upper azimuth angle of the closed curve corresponding to a certain point;
the singular point condition is a condition which should be satisfied by the singular point, the condition is to constrain the calculation process of the formula (1), and specifically, the following limiting conditions are added when the direction angle change is accumulated along the closed curve:
count the NψIf the direction is changed from positive to negative and from negative to positive once and only once, the Poincare index Poincare index value is continuously calculated, otherwise, the point is regarded as a common point;
count the NψIn a difference
If the number of the differences is more than one, the point is considered to be a common point;
if the final Poincare index value is 1/2, the detected point (i, j) is determined as core point, and if the Poincare index value is-1/2, the detected point (i, j) is determined as delta two-corner point;
(b) clustering the candidate singular points detected in the last step by using a clustering algorithm according to Euclidean distances, and counting the number N of the singular points contained in each clusterm，NmClusters < N are deleted; wherein, N represents the number of singular points that should be detected in the fingerprint mode area, and the number is stable, where N is 25;
(c) calculating the average mutation degree of the rest clusters, sorting the clusters according to the types of singular points and the mutation degrees, and selecting the first M clusters with larger mutation degreesThe clusters are used as candidate singular point clusters, and the candidate singular points are finally positioned on the mass centers of the clusters; the mutation degree refers to that in the process of making difference between two adjacent direction angles on the closed curve,
the average mutation degree refers to the average mutation degree of all pixels in each cluster;
(3) and a second stage: the false removing of the candidate singular point comprises the following steps:
for each candidate singular point, calculating the distribution consistency coherence value of each pixel in a circular neighborhood with the peripheral radius of 6 tau, wherein tau is the average streak line distance;
then dividing the circular area into 32 sectors and a central small circular area and calculating the average coherence value of each part, wherein the radius of the circular area and the width of two circular bands are both 2 tau;
finally, comparing the sizes of the average distribution consistency coherence values of the three regions in each direction in 16 directions, wherein the central circular region is a region common to the 16 directions, and if the distribution consistency coherence values from inside to outside are larger and larger, marking the direction as an effective direction; for a candidate singular point, if the number of effective directions is more than or equal to 10, extracting the point as the singular point, otherwise, considering the point as a pseudo point; if some directions of the singular points form the fingerprint boundary, the rest directions are used as reference directions, and effective directions are calculated in the directions; the specific process is as follows:
expressing the consistency information of the lines by the distribution consistency coherence value, wherein the distribution consistency coherence value is calculated by adopting Gaussian-Hermite special moment; four moments are used to describe the ridge consistency information of a fingerprint: m0，1，M1，0，M0，3，M3，0(ii) a And the four different orders of Gaussian-Hermite moments are defined as follows:
<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>M</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>λ</mi><msub><mi>M</mi><mn>1,0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>3,0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>λ</mi><msub><mi>M</mi><mn>0,1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>0,3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>
wherein lambda is the combined weight coefficient of Gaussian-Hermite special moments of different orders, wherein lambda is more than 0 and less than 1, and lambda takes 0.5; using equation (2), a feature vector [ M ] can be calculated for each pixel (i, j) in the fingerprint imageu，Mv]T(ii) a In fingerprint singular point region and nonsingular point region [ M ]u，Mv]THas different characteristics in the non-singular point region [ Mu，Mv]TDistributed along the direction perpendicular to the ridge line and in the singular point region [ Mu，Mv]TThen uniformly distributed in all directions, and adopting principal component analysis method to extract [ M [, M [ ]u，Mv]TDistribution characteristics of [ M ]u，Mv]TCovariance matrix CMIs defined by the formula:
<math><mrow><msub><mi>C</mi><mi>M</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><munder><mi>Σ</mi><mi>w</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mtd><mtd><munder><mi>Σ</mi><mi>u</mi></munder><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><munder><mi>Σ</mi><mi>w</mi></munder><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow></mtd><mtd><munder><mi>Σ</mi><mi>w</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>
wherein,
mum representing the average within the window WuValue, mvM representing the average within the window WvThe value n x n is the size of the window W, the value of n in the method is 4 tau +1, and tau is the average fringe distance;
let λ 1 and λ 2 be the covariance matrix CMWhen λ 1 > λ 2, [ M ] isu，Mv]TIs mainly distributed along the long axis direction, i.e. along the direction perpendicular to the ridge lineToward the distribution, while the values of λ 1 and λ 2 are very close in the noise region or singular point region, therefore, define [ Mu，Mv]TThe distribution consistency coherence is characterized as follows:
<math><mrow><mi>coherence</mi><mo>=</mo><mfrac><mrow><msup><msub><mi>λ</mi><mn>1</mn></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>λ</mi><mn>2</mn></msub><mn>2</mn></msup></mrow><mrow><msup><msub><mi>λ</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>λ</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>
<math><mrow><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><munder><mi>Σ</mi><mi>w</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munder><mi>Σ</mi><mi>w</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><msqrt><msup><mrow><mo>(</mo><munder><mi>Σ</mi><mi>u</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><munder><mi>Σ</mi><mi>u</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mrow><mo>(</mo><munder><mi>Σ</mi><mi>u</mi></munder><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow><mrow><msup><mrow><mo>(</mo><munder><mi>Σ</mi><mi>u</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><munder><mi>Σ</mi><mi>u</mi></munder><msup><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mrow><mo>(</mo><munder><mi>Σ</mi><mi>u</mi></munder><mrow><mo>(</mo><msub><mi>M</mi><mi>u</mi></msub><mo>-</mo><msub><mi>m</mi><mi>u</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>M</mi><mi>v</mi></msub><mo>-</mo><msub><mi>m</mi><mi>v</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>
thus, the better the line uniformity, the greater the distribution uniformity coherence.
2. The method as claimed in claim 1, wherein in the first stage (b), N is the number of singular points that should be detected in the fingerprint mode region, and the value 25 is based on the analysis of the test result, the number of singular points that can be detected in the singular point region of the fingerprint is relatively stable and is concentrated in the innermost part of the singular point region, and for the low-quality fingerprint image, as long as the singular point region is not polluted by the large noise, the number of singular points that can be detected in the singular point region is also stable, and N is 25 in the method.
3. The method as claimed in claim 1, wherein the window size of the coherence value of each pixel distribution in the second stage is (4 τ +1) × (4 τ + 1).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810138119XA CN101303729B (en) | 2008-07-01 | 2008-07-01 | A New Fingerprint Singularity Detection Method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810138119XA CN101303729B (en) | 2008-07-01 | 2008-07-01 | A New Fingerprint Singularity Detection Method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101303729A CN101303729A (en) | 2008-11-12 |
| CN101303729B true CN101303729B (en) | 2010-06-02 |
Family
ID=40113628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200810138119XA Expired - Fee Related CN101303729B (en) | 2008-07-01 | 2008-07-01 | A New Fingerprint Singularity Detection Method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101303729B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102460470B (en) * | 2009-06-11 | 2014-12-03 | 杜比实验室特许公司 | Trend analysis in content identification based on fingerprinting |
| CN102592276B (en) * | 2011-12-20 | 2014-07-09 | 中山大学 | Image area description method and copy image detection method based on same |
| CN103218616B (en) * | 2013-05-05 | 2016-04-13 | 西安电子科技大学 | Based on the image outline characteristic extraction method of Gauss-Hermite special moment |
| CN104866815B (en) * | 2015-04-23 | 2018-11-30 | 杭州电子科技大学 | Fingerprint epipole accurate positioning method based on image spatial feature |
| CN107024623B (en) * | 2016-02-02 | 2020-05-22 | 深圳市汇顶科技股份有限公司 | Test system and device of data acquisition chip and control method thereof |
| TWI606405B (en) * | 2016-05-30 | 2017-11-21 | 友達光電股份有限公司 | Image processing method and image processing system |
| CN106127753B (en) * | 2016-06-20 | 2019-07-30 | 中国科学院深圳先进技术研究院 | CT images body surface handmarking's extraction method in a kind of surgical operation |
| CN106447651A (en) * | 2016-09-07 | 2017-02-22 | 遵义师范学院 | Traffic sign detection method based on orthogonal Gauss-Hermite moment |
| CN107368780B (en) * | 2017-06-07 | 2020-08-28 | 西安电子科技大学 | Fingerprint registration point extraction method based on central singular point |
| CN116188024B (en) * | 2023-04-24 | 2023-07-04 | 山东蓝客信息科技有限公司 | Medical safety payment system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005038020A (en) * | 2003-07-15 | 2005-02-10 | Rikogaku Shinkokai | Fingerprint authentication device, computer system and network system |
| CN1327387C (en) * | 2004-07-13 | 2007-07-18 | 清华大学 | Method for identifying multi-characteristic of fingerprint |
| CN100342391C (en) * | 2003-12-24 | 2007-10-10 | 中国科学院自动化研究所 | Automatic fingerprint classification system and method |
| CN101114335A (en) * | 2007-07-19 | 2008-01-30 | 南京大学 | All-angle fast fingerprint recognition method |
| CN101145196A (en) * | 2006-09-13 | 2008-03-19 | 中国科学院自动化研究所 | A Fast Fingerprint Recognition Method Based on Singular Topological Structure |
-
2008
- 2008-07-01 CN CN200810138119XA patent/CN101303729B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005038020A (en) * | 2003-07-15 | 2005-02-10 | Rikogaku Shinkokai | Fingerprint authentication device, computer system and network system |
| CN100342391C (en) * | 2003-12-24 | 2007-10-10 | 中国科学院自动化研究所 | Automatic fingerprint classification system and method |
| CN1327387C (en) * | 2004-07-13 | 2007-07-18 | 清华大学 | Method for identifying multi-characteristic of fingerprint |
| CN101145196A (en) * | 2006-09-13 | 2008-03-19 | 中国科学院自动化研究所 | A Fast Fingerprint Recognition Method Based on Singular Topological Structure |
| CN101114335A (en) * | 2007-07-19 | 2008-01-30 | 南京大学 | All-angle fast fingerprint recognition method |
Non-Patent Citations (4)
| Title |
|---|
| 丁晋俊等.基于方向的指纹奇异点提取.计算机技术与发展17 2.2007,17(2),109-110,114. |
| 丁晋俊等.基于方向的指纹奇异点提取.计算机技术与发展17 2.2007,17(2),109-110,114. * |
| 王林等.基于Gaussian-Hermite矩的指纹奇异点定位.软件学报17 2.2006,17(2),242-249. |
| 王林等.基于Gaussian-Hermite矩的指纹奇异点定位.软件学报17 2.2006,17(2),242-249. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101303729A (en) | 2008-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101303729B (en) | 2010-06-02 | A New Fingerprint Singularity Detection Method |
| CN104318548B (en) | 2017-02-15 | Rapid image registration implementation method based on space sparsity and SIFT feature extraction |
| CN110189344B (en) | 2022-09-13 | Method for extracting straight line segment in electric power disconnecting link image |
| CN105335973B (en) | 2018-08-03 | Apply to the visual processing method of strip machining production line |
| CN104851095B (en) | 2017-12-26 | The sparse solid matching method of workpiece image based on modified Shape context |
| CN102073873B (en) | 2013-02-13 | Method for selecting SAR (spaceborne synthetic aperture radar) scene matching area on basis of SVM (support vector machine) |
| CN102629380B (en) | 2014-06-18 | Remote sensing image change detection method based on multi-group filtering and dimension reduction |
| CN108460780A (en) | 2018-08-28 | A kind of adhesion grain of rice image partition method based on background framework characteristic |
| CN104867139B (en) | 2018-02-09 | A kind of remote sensing image clouds and shadow detection method based on radiation field |
| CN101526994B (en) | 2011-02-02 | Fingerprint image segmentation method irrelevant to collecting device |
| CN104036480A (en) | 2014-09-10 | Surf algorithm based quick mismatching point eliminating method |
| CN104036516A (en) | 2014-09-10 | Camera calibration checkerboard image corner detection method based on symmetry analysis |
| CN104680509A (en) | 2015-06-03 | Real-time circular printing image defect detecting method |
| CN114169222A (en) | 2022-03-11 | Road surface condition discrimination method and road surface condition discrimination device |
| CN105513094A (en) | 2016-04-20 | Stereo vision tracking method and stereo vision tracking system based on 3D Delaunay triangulation |
| CN103218809A (en) | 2013-07-24 | Image measuring method of pearl length parameter |
| CN109781003B (en) | 2020-10-09 | A next best measurement pose determination method for structured light vision system |
| CN101840513A (en) | 2010-09-22 | Method for extracting image shape characteristics |
| CN103413116A (en) | 2013-11-27 | Effective fingerprint direction field calculating method |
| CN110210511A (en) | 2019-09-06 | A kind of improvement PCA-SIFT method for registering images based on cosine measure |
| CN101794442A (en) | 2010-08-04 | Calibration method for extracting illumination-insensitive information from visible images |
| CN106127677A (en) | 2016-11-16 | Surface in kind sampling point set boundary characteristic recognition methods based on fractional sample projected outline shape |
| CN109508674B (en) | 2021-08-13 | Airborne Down-View Heterogeneous Image Matching Method Based on Region Division |
| CN104036232B (en) | 2017-05-24 | Image edge feature analysis-based necktie pattern retrieval method |
| CN101430789A (en) | 2009-05-13 | Image edge detection method based on Fast Slant Stack transformation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 2008-11-12 | C06 | Publication | |
| 2008-11-12 | PB01 | Publication | |
| 2009-01-07 | C10 | Entry into substantive examination | |
| 2009-01-07 | SE01 | Entry into force of request for substantive examination | |
| 2010-06-02 | C14 | Grant of patent or utility model | |
| 2010-06-02 | GR01 | Patent grant | |
| 2014-08-20 | C17 | Cessation of patent right | |
| 2014-08-20 | CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100602 Termination date: 20130701 |