CN114821685A - PSO (particle swarm optimization) based non-contact 3D fingerprint identification method for optimizing Delaunay triangulation - Google Patents

PSO (particle swarm optimization) based non-contact 3D fingerprint identification method for optimizing Delaunay triangulation Download PDF

Info

Publication number
CN114821685A
CN114821685A CN202210520446.1A CN202210520446A CN114821685A CN 114821685 A CN114821685 A CN 114821685A CN 202210520446 A CN202210520446 A CN 202210520446A CN 114821685 A CN114821685 A CN 114821685A
Authority
CN
China
Prior art keywords
dimensional
fingerprint
minutiae
tetrahedron
threshold parameter
Prior art date
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.)
Granted
Application number
CN202210520446.1A
Other languages
Chinese (zh)
Other versions
CN114821685B (en
Inventor
杨东亮
宋昌江
孙思文
刘彤军
高凤娇
王云龙
费磊
丛晓丹
李磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute Of Intelligent Manufacturing Heilongjiang Academy Of Sciences
Original Assignee
Institute Of Intelligent Manufacturing Heilongjiang Academy Of Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute Of Intelligent Manufacturing Heilongjiang Academy Of Sciences filed Critical Institute Of Intelligent Manufacturing Heilongjiang Academy Of Sciences
Priority to CN202210520446.1A priority Critical patent/CN114821685B/en
Publication of CN114821685A publication Critical patent/CN114821685A/en
Application granted granted Critical
Publication of CN114821685B publication Critical patent/CN114821685B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/12Fingerprints or palmprints
    • G06V40/13Sensors therefor
    • G06V40/1312Sensors therefor direct reading, e.g. contactless acquisition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2415Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/25Design optimisation, verification or simulation using particle-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/12Fingerprints or palmprints
    • G06V40/1347Preprocessing; Feature extraction
    • G06V40/1359Extracting features related to ridge properties; Determining the fingerprint type, e.g. whorl or loop
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/12Fingerprints or palmprints
    • G06V40/1365Matching; Classification
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/12Fingerprints or palmprints
    • G06V40/1382Detecting the live character of the finger, i.e. distinguishing from a fake or cadaver finger
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Abstract

A PSO-based optimized Delaunay triangularization non-contact 3D fingerprint identification method belongs to the field of fingerprint identification algorithms. In order to make the recognition result of the algorithm of non-contact three-dimensional space fingerprint acquisition and recognition more accurate. Providing a PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method, and estimating a three-dimensional azimuth angle of a three-dimensional minutia point based on the improvement of a principal component analysis technology; constructing an expanded tetrahedral dissection structure; extracting features of the expanded tetrahedral dissection structure; comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons meeting the matching conditions and calculating matching scores, and determining whether the tetrahedrons are matched or not according to the matching scores; optimizing Delaunay triangulated fingerprint comparison algorithm parameters based on an improved particle swarm optimization PSO; and obtaining a final model. According to the invention, by optimizing the identification comparison parameters, a more accurate comparison model can be established, so that the identification precision of the whole fingerprint identification method is improved.

Description

PSO (particle swarm optimization) based non-contact 3D fingerprint identification method for optimizing Delaunay triangulation
Technical Field
The invention relates to a fingerprint identification method, in particular to a PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method.
Background
The traditional identity authentication methods based on passwords, passwords and the like have the defects of easy forgetting, easy loss and easy counterfeiting, and can not meet the requirements of people on identity authentication and information security. The identification technology based on the fingerprint characteristics is a more convenient and reliable authentication mode.
The traditional contact type fingerprint authentication mode adopts a contact type acquisition mode. The collection equipment that sets up in public places such as all kinds of halls of handling affairs, station carries out the fingerprint collection in-process, and collection equipment can be contacted by many users, and whether the user can worry before the user uses can infect at fingerprint collection surface bacterium, whether has sanitary problem. Moreover, the pressing force of the finger is different every time, the deformation degree of the finger is different, the definition degree of the finger fingerprint image is different, and the identification failure can be caused. Multiple contacts with the acquisition device are required to re-identify. In order to make up for the defects of the contact type acquisition and identification method, and on the basis of obtaining fundation of a natural science fund (No: LH2020F049) project of Heilongjiang province, the invention carries out related research design, establishes a non-contact type three-dimensional space form fingerprint acquisition and identification algorithm, and designs a non-contact type 3D fingerprint identification algorithm based on a particle swarm algorithm and expansion Delaunay triangulation in order to enable identification to be more accurate.
Disclosure of Invention
The invention aims to provide a non-contact 3D fingerprint identification method based on improved PSO (particle swarm optimization) optimized Delaunay triangularization, so that the identification result of a non-contact three-dimensional space fingerprint acquisition and identification algorithm is more accurate.
A non-contact 3D fingerprint identification method based on PSO optimized Delaunay triangulation is realized by the following steps:
estimating a three-dimensional azimuth angle of the three-dimensional minutiae point improved based on a principal component analysis technology;
constructing an expanded tetrahedron subdivision structure;
extracting features based on the expanded tetrahedral dissection structure;
wherein the extracted features include: edges, internal angles, two-dimensional direction angles, three-dimensional azimuth angles, and tetrahedral types;
comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons which accord with the matching conditions, counting the matching scores by using the counting results, and determining whether the tetrahedrons are matched or not according to the matching scores;
optimizing Delaunay triangulated fingerprint comparison algorithm parameters based on an improved particle swarm optimization PSO; optimizing and determining comparison parameters to obtain a final model;
delaunay triangularized Chinese means a triangulation algorithm.
Preferably, the step of estimating the three-dimensional azimuth angle of the three-dimensional minutiae point improved based on the principal component analysis technology specifically comprises:
(1) acquiring a two-dimensional fingerprint image;
(2) preprocessing the two-dimensional fingerprint image; wherein the preprocessing operation comprises: extracting a fingerprint foreground area, fingerprint ridge line and fingerprint minutiae;
(3) extracting two-dimensional minutiae information on the preprocessed two-dimensional fingerprint image, wherein the two-dimensional minutiae information comprises a coordinate x, a coordinate y and a two-dimensional direction angle theta;
(4) adding a space coordinate z and a three-dimensional azimuth angle phi to the two-dimensional minutiae points to map the two-dimensional minutiae points into three-dimensional minutiae points of a three-dimensional space;
the acquisition process of the space coordinate z is as follows: simulating to form a three-dimensional curved surface through a three-dimensional fingerprint, obtaining the three-dimensional curved surface, and obtaining a position attribute z through xy of a two-dimensional minutia;
(5) estimating a three-dimensional azimuth angle phi of the three-dimensional minutiae added in the last step, specifically:
setting a point F to represent a three-dimensional minutia point;
using all vertices v in m neighbourhood around the three-dimensional minutiae point F i ∈R 3 With these vertices v i As a basis, calculating the normal direction of the three-dimensional minutiae point F by using a principal component analysis method; order to
Figure BDA0003641347740000021
Wherein the content of the first and second substances,
Figure BDA0003641347740000022
t represents transposition; r 3 A three-dimensional space is represented in which,
Figure BDA0003641347740000023
denotes the median value, t denotes v i The number of (2); by
Figure BDA0003641347740000024
Formed covariance matrix
Figure BDA0003641347740000025
Wherein T represents transpose; w is a i The weight is represented by a weight that is,
Figure BDA0003641347740000026
improving an estimation method of the normal direction of a three-dimensional minutiae F through a covariance matrix C, wherein F represents the three-dimensional minutiae, a represents a weight parameter, and m represents the number of neighborhoods;
(6) setting the number of the fingerprint minutiae as k, and then extracting the fingerprint minutiae; the set of fingerprint minutiae is denoted F k =(x k ,y k ,z kkk ) Wherein (x) k ,y k ,z k ) Is a three-dimensional minutiae space coordinate, (theta) kk ) The two-dimensional direction angle of the three-dimensional detail point and the three-dimensional azimuth angle of the three-dimensional detail point are obtained;
the feature vectors of the template fingerprint minutiae are represented as: p ═ F i p I 1.. M }, the feature vector of the fingerprint minutiae to be identified is expressed as:
Figure BDA0003641347740000027
the template minutiae feature set P comprises M minutiae, and the minutiae feature set Q to be identified comprises N minutiae.
Preferably, said step of constructing an extended tetrahedrally subdivided structure is performed by,
the method comprises the following steps of reconstructing an expanded Delaunay tetrahedron subdivision structure by using point cloud Dironey triangulation, and specifically comprises the following steps:
(1) applying a Delaunay triangulation algorithm, and obtaining a tetrahedral structure G ═ { F, E } by adopting the Delaunay triangulation algorithm in a three-dimensional space, wherein E is an edge set of the tetrahedral structure, and F ═ { F ═ F } i 1.. N } represents a three-dimensional minutiae point set of a certain fingerprint; the four vertices of any one tetrahedron are denoted as (F) 1 ,F 2 ,F 3 ,F 4 ) Each vertex is denoted as F k =(x k ,y k ,z k ),k=1,2,3,4;
(2)F i Is a minutia in the F set, N i ={F j |(F i ,F j ) E represents in Delaunay triangulated structure G and with F i A detail point set formed by connecting vertexes;
(3) is obtained by Delaunay triangulation algorithm in N i Group G of tetrahedrons of i ={N i ,E i In which G is i Is the resulting tetrahedral set, E i Is a set of edges of a tetrahedral structure, N i Representing a set of vertices;
(4) the tetrahedron sets are put together to form an expanded tetrahedron set S, S ═ G- 1 ∪…∪G n And n is the number of minutiae.
Preferably, the step of extracting features of the expanded tetrahedron subdivision structure specifically includes:
let the feature of the tetrahedron be expressed as:
Figure BDA0003641347740000031
where, t is the number of this tetrahedral type,
Figure BDA0003641347740000032
denotes the largest of all internal angles in the tetrahedron, l max And l min Respectively representing the longest and shortest edges of the tetrahedron,
Figure BDA0003641347740000033
indicates the difference between the direction angles of the two vertexes of the longest side,
Figure BDA0003641347740000034
indicating the difference in the azimuth angle between the two vertices of the shortest side,
Figure BDA0003641347740000035
representing the difference in azimuth of the two vertices of the longest side,
Figure BDA0003641347740000036
representing the difference of the azimuth angles of the two vertexes of the shortest side;
the numbering method of the tetrahedron types comprises the following steps:
first of all, the first step is to,
Figure BDA0003641347740000037
is the largest of all the inner angles in the four faces, the largest of the triangular faces
Figure BDA0003641347740000038
Is marked as V 1
The vertex of the second largest interior angle in the triangular face is then labeled V 2 The vertex of the smallest interior angle in the triangular surface is marked as V 3 The remaining one vertex of the tetrahedron is marked V 4
Finally, the four vertex types are represented using vectors (t) 1 ,t 2 ,t 3 ,t 4 ),t i Is V i The type number Class of the tetrahedron can be obtained corresponding to the following table.
Preferably, the step of comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons which meet the matching conditions, counting the matching scores by using the counting result, and determining whether the tetrahedrons are matched according to the matching scores; the method specifically comprises the following steps:
(1) comparing each tetrahedron, and judging whether the tetrahedrons are matched:
setting two fingerprints to be compared as a fingerprint A and a fingerprint B, wherein any tetrahedral feature on the fingerprint A is represented as
Figure BDA0003641347740000039
Any tetrahedral feature on the fingerprint B is represented as
Figure BDA00036413477400000310
Judging the matching condition:
Figure BDA00036413477400000311
ΔL max =|L max -L' max |,Δθ max =|θ-θ'|,Δφ min =|φ-φ'|
if t ≠ t', the two tetrahedrons for comparison are considered to be unmatched, and comparison of other characteristics is not carried out any more;
if t is equal to t', continuously judging whether the t is satisfied
Figure BDA0003641347740000041
ΔL min <T L ,Δθ max <T θ ,Δθ min <T θ ,Δφ max <T φ ,Δφ min <T φ If yes, the two tetrahedrons are considered to be matched and recorded; wherein the content of the first and second substances,
Figure BDA0003641347740000042
T L 、T θ 、T φ respectively representing a threshold parameter of a maximum internal angle, a side length threshold parameter of a tetrahedron, a direction angle threshold parameter of two vertexes of an edge and an azimuth angle threshold parameter of the two vertexes of the edge, wherein the values are set;
(2) counting the number of matched minutiae:
if the tetrahedron is matched, marking the four minutiae of the tetrahedron as matching, and counting the number r of the matched minutiae;
(3) and (3) counting the matching scores:
calculating the matching score S of the obtained data r of the matched minutiae by using the formula:
Figure BDA0003641347740000043
wherein N is A Number of minutiae, N, for fingerprint A B Number of minutiae for fingerprint B;
(4) and (3) recognition results:
matching score S and a set matching score threshold value T S Making a comparison when S>T S If so, matching the fingerprint A with the fingerprint B, otherwise, not matching the two fingerprints; wherein, T S Indicating a match score threshold to be achieved, the value of which is set.
Preferably, the improved particle swarm optimization PSO-based optimized Delaunay triangularization fingerprint comparison calculationA method parameter; optimizing and determining the comparison parameters to obtain the final model, wherein the fingerprint comparison algorithm parameters comprise a threshold parameter for optimizing the maximum internal angle
Figure BDA0003641347740000046
Edge length threshold parameter T of tetrahedron L Direction angle threshold parameter T of two vertexes of an edge θ The azimuth angle threshold parameter T of two vertexes of the edge φ And a matching score threshold T S (ii) a The method specifically comprises the following steps:
establishing a plurality of matched fingerprint pairs and unmatched fingerprint pairs as training sets, and acquiring optimal parameters by using a particle swarm algorithm PSO; utilizing particle swarm optimization PSO to carry out threshold parameter on maximum internal angle
Figure BDA0003641347740000047
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Optimizing the five parameters; wherein, the optimizing process is as follows:
(1) initializing particle swarm size to be 40, iterating times to be 200, searching threshold parameters of optimal maximum internal angle
Figure BDA0003641347740000048
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Defining the position and velocity of the particle swarm;
(2) determining fitness function
Figure BDA0003641347740000044
Wherein N represents the number of training samples, and ACC represents the identification accuracy of the training samples;
(3) updating individual optimal positions p im And the population optimum position p gm
(4) Updating the speed and position of the particle;
Figure BDA0003641347740000045
wherein, ω represents the inertial weight; m is an element of [1, M ]],i∈[1,40];v im And x im Respectively representing the speed and the position of the m-th particle; c. C 1 And c 2 Is a learning factor; r is 1 ,r 2 ∈[0,1]Is a random number; p is a radical of formula im And p gm Respectively representing the individual optimal position and the group optimal position searched by the i particle at present; updating particle positions
Figure BDA0003641347740000051
(5) If the error is less than 0.01 or the maximum iteration times are reached, ending the optimization of the particle swarm optimization algorithm to obtain the optimal particle position as the threshold parameter of the optimal maximum internal angle
Figure BDA0003641347740000052
Optimal length threshold parameter T L Optimal direction angle threshold parameter T θ Optimal azimuth angle threshold parameter T φ Optimal matching score threshold T S Therefore, an optimized three-dimensional fingerprint comparison algorithm based on the extended Delaunay triangulation is obtained.
The invention has the beneficial effects that:
the invention adopts non-contact type to collect fingerprints, thereby avoiding cross infection among users.
The fingerprint identification algorithm designed by the invention has the advantage of high operation speed, and the time required by identification is reduced.
Firstly, the invention carries out a novel tetrahedron type coding method, and particularly relates to a process of adding a space coordinate z and a three-dimensional azimuth angle phi to a two-dimensional minutia to map the two-dimensional minutia into a three-dimensional minutia of a three-dimensional space, then constructing an expanded tetrahedron subdivision structure and extracting characteristics, and designing a unique tetrahedron type number in the process. In the subsequent comparison process, once the numbers of the tetrahedrons are not matched, the comparison of other information is not carried out, so that the method makes the model matching more efficient. Therefore, the invention has higher calculation efficiency and reduces 60-70% of the time required for identification.
In addition, the invention adds the step of optimizing the Delaunay triangularization fingerprint comparison algorithm parameters based on the improved PSO, and is used for optimizing the identification comparison parameters to establish a more accurate comparison model, thereby improving the identification precision of the whole fingerprint identification method, and the identification precision reaches more than 98%. Specifically, the point cloud Dironi triangulation is introduced to obtain a tetrahedron of the three-dimensional fingerprint minutiae, and matching is performed on the obtained tetrahedron of the three-dimensional fingerprint minutiae. This can greatly reduce the amount of computation, but when false minutiae occur or minutiae are lost, it can cause tetrahedral structure errors, and thus tetrahedral matching errors. In order to overcome the problem, an expanded tetrahedron set is provided, and the robustness in tetrahedron matching is improved.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic diagram of an expanded tetrahedron subdivision scheme in accordance with the present invention;
FIG. 3 is a schematic illustration of a tetrahedron in the extended tetrahedron subdivision arrangement of the present invention;
FIG. 4 is a flow chart of a three-dimensional fingerprint matching algorithm according to the present invention;
FIG. 5 is an example of features involved in the step of extracting features involved in the present invention;
FIG. 6 is an example of a 2D fingerprint image to which the present invention relates;
FIG. 7 is a diagram of the extraction of two-dimensional minutiae points on a two-dimensional fingerprint image according to the present invention;
FIG. 8 is a ROC curve relating to a simulation experiment of the present invention.
Detailed Description
The first specific implementation way is as follows:
in the embodiment, as shown in fig. 1, the non-contact 3D fingerprint identification method based on PSO optimized Delaunay triangulation is implemented by the following steps:
step one, because the geometric coordinates of the vertex of the three-dimensional fingerprint surface are easily disturbed by noise, the normal direction of the vertex can be changed, and the estimation error of the three-dimensional azimuth angle of the three-dimensional minutiae point is larger. The invention is presented herein:
estimating a three-dimensional azimuth angle of the three-dimensional minutiae point improved based on a principal component analysis technology; the method is used for improving the precision of three-dimensional azimuth estimation;
step two, constructing an expanded tetrahedron subdivision structure;
step three, extracting characteristics based on the expanded tetrahedral subdivision structure;
wherein the extracted features include: edges, internal angles, two-dimensional direction angles, three-dimensional azimuth angles, and tetrahedral types;
as shown in FIG. 5, the azimuth angle is θ and the azimuth angle φ, the azimuth angle θ is the direction on the X-Y plane; the azimuth angle phi is a direction in space;
comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons which accord with the matching conditions, counting the matching scores by using the counting results, and determining whether the tetrahedrons are matched or not according to the matching scores;
fifthly, optimizing Delaunay triangularized fingerprint comparison algorithm parameters based on an improved Particle Swarm Optimization (PSO); optimizing and determining comparison parameters to obtain a final model;
delaunay triangularized Chinese means a triangulation algorithm.
The method provided by the invention verifies the identification accuracy through the acquired three-dimensional data. The simulation result shows that the method has the advantages of high operation speed and short identification time.
The invention provides a 3D fingerprint identification algorithm based on improved particle swarm optimization and expanded Delaunay triangulation, which has the following effects:
the expanded Delaunay triangulation is proposed to perform tetrahedrization on the three-dimensional fingerprint minutiae, so that the influence of wrong minutiae on the whole fingerprint triangulation is reduced, and the robustness of a tetrahedral structure is improved.
And 2, optimizing parameters in the three-dimensional fingerprint matching model by utilizing a particle swarm optimization algorithm to obtain the parameters of the optimal matching model, thereby greatly improving the identification precision of the model.
3, the novel tetrahedron type coding makes the model matching more efficient and has higher calculation efficiency.
The second embodiment is as follows:
different from the specific embodiment, in the non-contact 3D fingerprint identification method based on PSO optimized Delaunay triangulation of the present embodiment, the step of estimating the three-dimensional azimuth angle of the three-dimensional minutiae based on the principal component analysis technology improvement specifically includes:
(1) acquiring a two-dimensional fingerprint image; the two-dimensional fingerprint image is shown in FIG. 6;
(2) preprocessing the two-dimensional fingerprint image; wherein the preprocessing operation comprises: extracting a fingerprint foreground area, fingerprint ridge line and fingerprint minutiae;
(3) extracting two-dimensional minutiae information on the preprocessed two-dimensional fingerprint image, wherein the two-dimensional minutiae information comprises a coordinate x, a coordinate y and a two-dimensional direction angle theta; wherein the two-dimensional minutiae plot is shown in FIG. 7;
(4) adding a space coordinate z and a three-dimensional azimuth angle phi to the two-dimensional minutiae points to map the two-dimensional minutiae points into three-dimensional minutiae points of a three-dimensional space;
the acquisition process of the space coordinate z is as follows: simulating to form a three-dimensional curved surface through a three-dimensional fingerprint, obtaining the three-dimensional curved surface, and obtaining a position attribute z through xy of a two-dimensional minutia;
(5) estimating a three-dimensional azimuth angle phi of the three-dimensional minutiae added in the last step, specifically:
in order to improve the estimation accuracy of the three-dimensional azimuth angle of the three-dimensional minutiae, setting a point F to represent the three-dimensional minutiae;
using all vertices v in m neighbourhood around the three-dimensional minutiae point F i ∈R 3 With these vertices v i As a basis, calculating the normal direction of the three-dimensional minutiae point F by using a principal component analysis method; order to
Figure BDA0003641347740000071
Wherein the content of the first and second substances,
Figure BDA0003641347740000072
t represents transposition; r 3 To representThe three-dimensional space is a space with three-dimensional,
Figure BDA0003641347740000077
denotes the median value, t denotes v i The number of (2); by
Figure BDA0003641347740000073
Formed covariance matrix
Figure BDA0003641347740000074
Wherein T represents transpose; w is a i The weight is represented by a weight that is,
Figure BDA0003641347740000075
improving an estimation method of the normal direction of a three-dimensional minutiae F through a covariance matrix C, wherein F represents the three-dimensional minutiae, a represents a weight parameter, and m represents the number of neighborhoods;
(6) setting the number of the fingerprint minutiae as k, and then extracting the fingerprint minutiae by the conventional method for extracting the fingerprint minutiae; the set of fingerprint minutiae is denoted F k =(x k ,y k ,z kkk ) Wherein (x) k ,y k ,z k ) Is a three-dimensional minutiae space coordinate, (theta) kk ) The two-dimensional direction angle of the three-dimensional minutiae and the three-dimensional azimuth angle of the three-dimensional minutiae are obtained;
the feature vectors of the template fingerprint minutiae are represented as: p ═ F i p I 1.. M }, the feature vector of the fingerprint minutiae to be identified is expressed as:
Figure BDA0003641347740000076
the template minutiae feature set P comprises M minutiae, and the minutiae feature set Q to be identified comprises N minutiae.
The third concrete implementation mode:
different from the specific embodiment, in the non-contact 3D fingerprint identification method based on PSO optimized Delaunay triangulation according to the embodiment, the step of constructing the extended tetrahedron subdivision structure specifically,
the matching problem of minutiae can be understood as the best correspondence between the elements in the set P and the set Q, resulting in a matching score. The existing matching method usually needs to align elements in two sets, and the alignment operation has high calculation complexity. This reduces the computation by 60-70%, but when false minutiae occur or minutiae are missing, it leads to tetrahedral structure errors and thus tetrahedral matching errors. In order to overcome the problem, an expanded tetrahedron set is provided, and the robustness in tetrahedron matching is improved.
The method comprises the following steps of reconstructing an expanded Delaunay tetrahedron subdivision structure by using point cloud Dironey triangulation, and specifically comprises the following steps:
(1) applying the existing Delaunay triangulation algorithm, and obtaining a tetrahedral structure G (F, E) by adopting the Delaunay triangulation algorithm in a three-dimensional space, wherein E is an edge set of the tetrahedral structure, and F (F) i N represents a three-dimensional minutiae set of a certain fingerprint; the four vertices of any one tetrahedron are denoted as (F) 1 ,F 2 ,F 3 ,F 4 ) Each vertex is denoted as F k =(x k ,y k ,z k ),k=1,2,3,4;
(2)F i Is a minutia in the F set, N i ={F j |(F i ,F j ) E represents in Delaunay triangulated structure G and with F i A detail point set formed by connecting vertexes;
(3) is obtained by Delaunay triangulation algorithm in N i Group G of tetrahedrons of i ={N i ,E i In which G is i Is the resulting tetrahedral set, E i Is a set of edges of a tetrahedral structure, N i Representing a set of vertices;
(4) putting these tetrahedron sets together to form an expanded tetrahedron set S, S ═ G- 1 ∪…∪G n And n is the number of minutiae.
The extended tetrahedron set S includes not only the original tetrahedron join G, but also the set N of vertices i A generated tetrahedron set; higher matching scores may also be obtained when a loss of minutiae occurs.
The fourth concrete implementation mode:
different from the third specific embodiment, in the non-contact 3D fingerprint identification method based on PSO optimization Delaunay triangulation according to the third embodiment, the step of extracting features based on an extended tetrahedral subdivision structure specifically includes:
let the feature of the tetrahedron be expressed as:
Figure BDA0003641347740000081
where, t is the number of this tetrahedral type,
Figure BDA0003641347740000082
denotes the largest of all internal angles in the tetrahedron, l max And l min Respectively representing the longest and shortest edges of the tetrahedron,
Figure BDA0003641347740000083
indicates the difference between the direction angles of the two vertexes of the longest side,
Figure BDA0003641347740000084
indicating the difference in the azimuth angle between the two vertices of the shortest side,
Figure BDA0003641347740000085
representing the difference in azimuth of the two vertices of the longest side,
Figure BDA0003641347740000086
representing the difference of the azimuth angles of the two vertexes of the shortest side;
the tetrahedron type number t is a new type table and represents the combination type of the detail nodes, and each detail node can be divided into two types of end points or branch points; each tetrahedron can construct 16 types of numbers according to the detail point types of four vertexes of each tetrahedron;
the numbering method of the tetrahedron types comprises the following steps:
first of all, the first step is to,
Figure BDA0003641347740000091
is the largest of all the inner angles in the four faces, the largest of the triangular faces
Figure BDA0003641347740000092
Is marked as V 1
The vertex of the second largest interior angle in the triangular face is then labeled V 2 The vertex of the smallest interior angle in the triangular surface is marked as V 3 The remaining one vertex of the tetrahedron is marked V 4
Finally, the four vertex types are represented using vectors (t) 1 ,t 2 ,t 3 ,t 4 ),t i Is V i The type number Class of the tetrahedron can be obtained corresponding to the following table. For example, assuming that the type of four vertices of a tetrahedron is represented as (e, e, e, b), e represents that the point type is an end point, and b represents that the point type is a bifurcation point, the tetrahedron has a corresponding type number of 2.
The fifth concrete implementation mode:
different from the first, second, or fourth specific embodiments, in the PSO-optimized Delaunay triangulation-based non-contact 3D fingerprint identification method of the present embodiment, the steps of comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons that meet the matching conditions, counting the matching scores according to the counting results, and determining whether to match according to the matching scores are specifically:
(1) comparing each tetrahedron, and judging whether the tetrahedrons are matched:
setting two fingerprints to be compared as a fingerprint A and a fingerprint B, wherein any tetrahedral feature on the fingerprint A is represented as
Figure BDA0003641347740000093
On the fingerprint BAny one of the tetrahedral features is represented as
Figure BDA0003641347740000094
Judging the matching condition:
Figure BDA0003641347740000095
ΔL max =|L max -L' max |,Δθ max =|θ-θ'|,Δφ min =|φ-φ'|
if t ≠ t', the two tetrahedrons for comparison are considered to be unmatched, and comparison of other characteristics is not carried out any more;
if t is equal to t', continuously judging whether the t is satisfied
Figure BDA0003641347740000096
ΔL max <T L ,ΔL min <T L ,Δθ max <T θ ,Δθ min <T θ ,Δφ max <T φ ,Δφ min <T φ If yes, the two tetrahedrons are considered to be matched and recorded; wherein the content of the first and second substances,
Figure BDA0003641347740000098
T L 、T θ 、T φ respectively representing a threshold parameter of a maximum internal angle, a side length threshold parameter of a tetrahedron, a direction angle threshold parameter of two vertexes of an edge and an azimuth angle threshold parameter of the two vertexes of the edge, wherein the values are set;
(2) counting the number of matched minutiae:
if the tetrahedron is matched, marking the four minutiae of the tetrahedron as matching, and counting the number r of the matched minutiae;
(3) and (3) counting the matching scores:
calculating the matching score S of the obtained data r of the matched minutiae by using the formula:
Figure BDA0003641347740000097
wherein N is A Number of minutiae for fingerprint A,N B Number of minutiae for fingerprint B;
(4) and (3) recognition results:
matching score S and a set matching score threshold value T S The comparison is performed, as shown in fig. 4, in which the threshold is set in the process of comparing two fingerprints. When S is>T S If so, matching the fingerprint A with the fingerprint B, otherwise, not matching the two fingerprints; wherein, T S Indicating a match score threshold to be achieved, the value of which is set.
The sixth specific implementation mode:
different from the fifth embodiment, in the non-contact 3D fingerprint identification method based on PSO optimized Delaunay triangulation of the present embodiment, the improved particle swarm optimization PSO optimized Delaunay triangulation based fingerprint comparison algorithm parameters; in the step of optimizing and determining the comparison parameters to obtain the final model, the fingerprint comparison algorithm parameters include a threshold parameter for optimizing the maximum internal angle
Figure BDA0003641347740000106
Edge length threshold parameter T of tetrahedron L Direction angle threshold parameter T of two vertexes of an edge θ The azimuth angle threshold parameter T of two vertexes of the edge φ And a matching score threshold T S (ii) a The method specifically comprises the following steps:
threshold parameter due to maximum internal angle
Figure BDA0003641347740000101
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S The accuracy of the image recognition algorithm is set, so that a plurality of matched fingerprint pairs and unmatched fingerprint pairs are established as a training set, and the Particle Swarm Optimization (PSO) is used for obtaining the optimal parameters; utilizing particle swarm optimization PSO to carry out threshold parameter on maximum internal angle
Figure BDA0003641347740000107
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Optimizing the five parameters, and further improving the accuracy of fingerprint identification based on the extended Delaunay triangulation; wherein, the optimizing process is as follows:
(1) combining the requirement of a Delaunay triangularization fingerprint comparison algorithm, initializing the particle swarm size to be 40, the iteration number to be 200, and searching for the threshold parameter of the optimal maximum internal angle
Figure BDA0003641347740000102
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Defining the position and velocity of the particle swarm;
(2) determining fitness function
Figure BDA0003641347740000103
Wherein N represents the number of training samples, and ACC represents the identification accuracy of the training samples;
(3) updating individual optimal positions p im And the population optimum position p gm
(4) Updating the speed and position of the particle;
Figure BDA0003641347740000104
wherein, ω represents the inertial weight; m is an element of [1, M ]],i∈[1,40];v im And x im Respectively representing the speed and the position of the m-th particle; c. C 1 And c 2 Is a learning factor; r is 1 ,r 2 ∈[0,1]Is a random number; p is a radical of formula im And p gm Respectively representing the individual optimal position and the group optimal position searched by the i particle at present; updating particle positions
Figure BDA0003641347740000105
(5) If the error is less than 0.01 or the maximum iteration times are reached, ending the optimization of the particle swarm optimization algorithm to obtain the optimal particle position as the threshold parameter of the optimal maximum internal angle
Figure BDA0003641347740000111
Optimal length threshold parameter T L Optimal direction angle threshold parameter T θ Optimal azimuth angle threshold parameter T φ Optimal matching score threshold T S Therefore, an optimized three-dimensional fingerprint comparison algorithm based on the extended Delaunay triangulation is obtained.
Simulation verification:
collecting experimental data:
a structured light three-dimensional fingerprint acquisition system which is designed and built by taking a DLP4500 type three-dimensional scanner and a Grasshopper3 camera as main components is used for acquiring a large amount of three-dimensional fingerprint data through the system, and the identification accuracy of the method is verified through the acquired three-dimensional data. As shown in fig. 8, the change in receiver operating characteristic curve (ROC) is shown. The first method is a point matching method, and the second method is a fingerprint identification method related to the application. It can be found from experiments that in the collected data set, method one and method two, their corresponding Equal Error Rates (EER) are 0.16 and 0.13, respectively. The fingerprint identification method has the advantages of lower equal error rate and better identification performance.
The embodiments of the present invention are disclosed as the preferred embodiments, but not limited thereto, and those skilled in the art can easily understand the spirit of the present invention and make various extensions and changes without departing from the spirit of the present invention.

Claims (6)

1. A non-contact 3D fingerprint identification method based on PSO optimization Delaunay triangulation is characterized in that: the method is realized by the following steps:
estimating a three-dimensional azimuth angle of the three-dimensional minutiae point improved based on a principal component analysis technology;
constructing an expanded tetrahedron subdivision structure;
extracting features based on the expanded tetrahedral dissection structure;
wherein the extracted features include: edges, internal angles, two-dimensional direction angles, three-dimensional azimuth angles, and tetrahedral types;
comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons which accord with the matching conditions, counting the matching scores by using the counting results, and determining whether the tetrahedrons are matched or not according to the matching scores;
optimizing the fingerprint comparison algorithm parameters of Delaunay triangulation based on the improved particle swarm optimization PSO; optimizing and determining comparison parameters to obtain a final model;
the Delaunay triangularization Chinese meaning is a triangulation algorithm.
2. The PSO-based optimized Delaunay triangulation-based non-contact 3D fingerprint identification method as claimed in claim 1, wherein: the three-dimensional azimuth angle estimation method based on the improved three-dimensional minutiae of the principal component analysis technology specifically comprises the following steps:
(1) acquiring a two-dimensional fingerprint image;
(2) preprocessing the two-dimensional fingerprint image; wherein the preprocessing operation comprises: extracting a fingerprint foreground area, fingerprint ridge line extraction and fingerprint minutiae;
(3) extracting two-dimensional minutiae information on the preprocessed two-dimensional fingerprint image, wherein the two-dimensional minutiae information comprises a coordinate x, a coordinate y and a two-dimensional direction angle theta;
(4) adding a space coordinate z and a three-dimensional azimuth angle phi to the two-dimensional minutiae points to map the two-dimensional minutiae points into three-dimensional minutiae points of a three-dimensional space;
the acquisition process of the space coordinate z is as follows: simulating to form a three-dimensional curved surface through a three-dimensional fingerprint, obtaining the three-dimensional curved surface, and obtaining a position attribute z through xy of a two-dimensional minutia;
(5) estimating a three-dimensional azimuth angle phi of the three-dimensional minutiae added in the last step, specifically:
setting a point F to represent a three-dimensional minutia point;
using all vertices v in m neighbourhood around the three-dimensional minutiae point F i ∈R 3 With these vertices v i As a basis, calculating the normal direction of the three-dimensional minutiae point F by using a principal component analysis method; order to
Figure FDA0003641347730000011
Wherein the content of the first and second substances,
Figure FDA0003641347730000012
t represents transposition; r 3 A three-dimensional space is represented in which,
Figure FDA0003641347730000013
denotes the median value, t denotes v i The number of (2); by
Figure FDA0003641347730000014
Formed covariance matrix
Figure FDA0003641347730000015
Wherein T represents transpose; w is a i The weight is represented by a weight that is,
Figure FDA0003641347730000016
improving an estimation method of the normal direction of a three-dimensional minutiae F through a covariance matrix C, wherein F represents the three-dimensional minutiae, a represents a weight parameter, and m represents the number of neighborhoods;
(6) setting the number of the fingerprint minutiae as k, and then extracting the fingerprint minutiae; the set of fingerprint minutiae is denoted F k =(x k ,y k ,z kkk ) Wherein (x) k ,y k ,z k ) Is a three-dimensional minutiae space coordinate, (theta) kk ) The two-dimensional direction angle of the three-dimensional minutiae and the three-dimensional azimuth angle of the three-dimensional minutiae are obtained;
the feature vectors of the template fingerprint minutiae are represented as: p ═ F i p I 1.. M }, the feature vector of the fingerprint minutiae to be identified is expressed as:
Figure FDA0003641347730000021
the template minutiae feature set P comprises M minutiae, and the feature set Q of the minutiae to be identified comprises N minutiaeAnd (5) minutiae points.
3. The PSO-based optimized Delaunay triangulation-based non-contact 3D fingerprint identification method as claimed in claim 2, wherein: the process of the step of constructing an extended tetrahedrally subdivided structure is described,
the method comprises the following steps of reconstructing an expanded Delaunay tetrahedron subdivision structure by using point cloud Dironey triangulation, and specifically comprises the following steps:
(1) applying a Delaunay triangulation algorithm, and obtaining a tetrahedral structure G ═ { F, E } by adopting the Delaunay triangulation algorithm in a three-dimensional space, wherein E is an edge set of the tetrahedral structure, and F ═ { F ═ F } i N represents a three-dimensional minutiae set of a certain fingerprint; the four vertices of any one tetrahedron are denoted as (F) 1 ,F 2 ,F 3 ,F 4 ) Each vertex is denoted as F k =(x k ,y k ,z k ),k=1,2,3,4;
(2)F i Is a minutia in the F set, N i ={F j |(F i ,F j ) E represents in Delaunay triangulated structure G and with F i A detail point set formed by connecting vertexes;
(3) is obtained by Delaunay triangulation algorithm i Group G of tetrahedrons of i ={N i ,E i In which G is i Is the resulting tetrahedral set, E i Is a set of edges of a tetrahedral structure, N i Representing a set of vertices;
(4) the tetrahedron sets are put together to form an expanded tetrahedron set S, S ═ G- 1 ∪…∪G n And n is the number of minutiae.
4. The PSO-based optimized Delaunay triangulation-based non-contact 3D fingerprint identification method as claimed in claim 3, wherein: the step of extracting features of the expanded tetrahedron subdivision structure specifically comprises the following steps:
let the feature of the tetrahedron be expressed as:
Figure FDA0003641347730000022
where, t is the number of this tetrahedral type,
Figure FDA0003641347730000023
denotes the largest of all internal angles in the tetrahedron, l max And l min Respectively representing the longest and shortest edges of the tetrahedron,
Figure FDA0003641347730000024
indicates the difference between the direction angles of the two vertexes of the longest side,
Figure FDA0003641347730000025
indicating the difference in the azimuth angle between the two vertices of the shortest side,
Figure FDA0003641347730000026
representing the difference in azimuth of the two vertices of the longest side,
Figure FDA0003641347730000027
representing the difference of the azimuth angles of the two vertexes of the shortest side;
the numbering method of the tetrahedron types comprises the following steps:
first of all, when a user wants to use the apparatus,
Figure FDA0003641347730000031
is the largest of all the inner angles in the four faces, the largest of the triangular faces
Figure FDA0003641347730000032
Is marked as V 1
The vertex of the second largest interior angle in the triangular face is then labeled V 2 The vertex of the smallest interior angle in the triangular surface is marked as V 3 One vertex remaining in the tetrahedron is marked V 4
Finally, the four vertex types are represented using vectors (t) 1 ,t 2 ,t 3 ,t 4 ),t i Is V i The type number Class of the tetrahedron can be obtained corresponding to the following table.
5. The PSO-based optimized Delaunay triangulation-based non-contact 3D fingerprint identification method as claimed in claim 4, wherein: the method comprises the following steps of comparing each tetrahedron to be matched according to the extracted features, counting the tetrahedrons which accord with the matching conditions, counting the matching scores by using the counting results, and determining whether the tetrahedrons are matched according to the matching scores, wherein the steps are as follows:
(1) comparing each tetrahedron, and judging whether the tetrahedrons are matched:
setting two fingerprints to be compared as a fingerprint A and a fingerprint B, wherein any tetrahedral feature on the fingerprint A is represented as
Figure FDA0003641347730000033
Any tetrahedral feature on the fingerprint B is represented as
Figure FDA0003641347730000034
Judging the matching condition:
Figure FDA0003641347730000035
ΔL max =|L max -L' max |,Δθ max =|θ-θ'|,Δφ min =|φ-φ'|
if t ≠ t', the two tetrahedrons for comparison are considered to be unmatched, and comparison of other characteristics is not carried out any more;
if t is equal to t', continuously judging whether the t is satisfied
Figure FDA0003641347730000036
ΔL max <T L ,ΔL min <T L ,Δθ max <T θ ,Δθ min <T θ ,Δφ max <T φ ,Δφ min <T φ If so, the two are consideredThe tetrahedra are matched and recorded; wherein the content of the first and second substances,
Figure FDA0003641347730000037
T L 、T θ 、T φ respectively representing a threshold parameter of a maximum internal angle, a side length threshold parameter of a tetrahedron, a direction angle threshold parameter of two vertexes of an edge and an azimuth angle threshold parameter of the two vertexes of the edge, wherein the values are set;
(2) counting the number of matched minutiae:
if the tetrahedron is matched, marking the four minutiae of the tetrahedron as matching, and counting the number r of the matched minutiae;
(3) and (3) counting the matching scores:
calculating the matching score S of the obtained data r of the matched minutiae by using the formula:
Figure FDA0003641347730000038
wherein N is A Number of minutiae, N, for fingerprint A B Number of minutiae for fingerprint B;
(4) and (3) recognition results:
matching score S and a set matching score threshold value T S Making a comparison when S > T S If so, matching the fingerprint A with the fingerprint B, otherwise, not matching the two fingerprints; wherein, T S Indicating a match score threshold to be achieved, the value of which is set.
6. The PSO-based optimized Delaunay triangulation-based non-contact 3D fingerprint identification method as claimed in claim 5, wherein: the improved particle swarm optimization PSO is used for optimizing the Delaunay triangularization fingerprint comparison algorithm parameters; in the step of optimizing and determining the comparison parameters to obtain the final model, the fingerprint comparison algorithm parameters include a threshold parameter for optimizing the maximum internal angle
Figure FDA0003641347730000041
Edge length threshold parameter T of tetrahedron L The direction angle threshold value of two vertexes of the edgeNumber T θ The azimuth angle threshold parameter T of two vertexes of the edge φ And a matching score threshold T S (ii) a The method specifically comprises the following steps:
establishing a plurality of matched fingerprint pairs and unmatched fingerprint pairs as training sets, and acquiring optimal parameters by using a particle swarm algorithm PSO; utilizing particle swarm optimization PSO to carry out threshold parameter on maximum internal angle
Figure FDA0003641347730000042
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Optimizing the five parameters; wherein, the optimizing process is as follows:
(1) initializing particle swarm size to be 40, iterating times to be 200, searching threshold parameters of optimal maximum internal angle
Figure FDA0003641347730000043
Length threshold parameter T L Direction angle threshold parameter T θ Azimuth angle threshold parameter T φ Matching score threshold T S Defining the position and velocity of the particle swarm;
(2) determining fitness function
Figure FDA0003641347730000044
Wherein N represents the number of training samples, and ACC represents the identification accuracy of the training samples;
(3) updating individual optimal positions p im And the population optimum position p gm
(4) Updating the speed and position of the particle;
Figure FDA0003641347730000045
wherein, ω represents the inertial weight; m is an element of [1, M ]],i∈[1,40];v im And x im Respectively representing the speed and the position of the m-th particle; c. C 1 And c 2 Is a learning factor; r is 1 ,r 2 ∈[0,1]Is a random number; p is a radical of im And p gm Respectively represent i particlesThe current searched individual optimal position and the group optimal position; updating particle positions
Figure FDA0003641347730000046
(5) If the error is less than 0.01 or the maximum iteration times are reached, ending the optimization of the particle swarm optimization algorithm to obtain the optimal particle position as the threshold parameter of the optimal maximum internal angle
Figure FDA0003641347730000047
Optimal length threshold parameter T L Optimal direction angle threshold parameter T θ Optimal azimuth angle threshold parameter T φ Optimal matching score threshold T S Therefore, an optimized three-dimensional fingerprint comparison algorithm based on the extended Delaunay triangulation is obtained.
CN202210520446.1A 2022-05-12 2022-05-12 PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method Active CN114821685B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210520446.1A CN114821685B (en) 2022-05-12 2022-05-12 PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210520446.1A CN114821685B (en) 2022-05-12 2022-05-12 PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method

Publications (2)

Publication Number Publication Date
CN114821685A true CN114821685A (en) 2022-07-29
CN114821685B CN114821685B (en) 2023-05-02

Family

ID=82515967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210520446.1A Active CN114821685B (en) 2022-05-12 2022-05-12 PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method

Country Status (1)

Country Link
CN (1) CN114821685B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101620677A (en) * 2009-03-11 2010-01-06 刘鸣宇 Fingerprint identification method based on triangulation and LOD technology
US20140044322A1 (en) * 2012-08-08 2014-02-13 The Hong Kong Polytechnic University Contactless 3D Biometric Feature identification System and Method thereof
CN104009973A (en) * 2014-05-08 2014-08-27 电子科技大学 Fingerprint detail information hiding and recovering method based on set polynomial conversion and harmonics
CN110096634A (en) * 2019-04-29 2019-08-06 成都理工大学 A kind of house property data vector alignment schemes based on particle group optimizing
CN112487867A (en) * 2020-11-03 2021-03-12 杭州电子科技大学 Visual constraint fingerprint identification method based on enhanced triangulation
WO2021048174A1 (en) * 2019-09-10 2021-03-18 Thales Dis France Sa Method for determining a match between a candidate fingerprint and a reference fingerprint

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101620677A (en) * 2009-03-11 2010-01-06 刘鸣宇 Fingerprint identification method based on triangulation and LOD technology
US20140044322A1 (en) * 2012-08-08 2014-02-13 The Hong Kong Polytechnic University Contactless 3D Biometric Feature identification System and Method thereof
CN104009973A (en) * 2014-05-08 2014-08-27 电子科技大学 Fingerprint detail information hiding and recovering method based on set polynomial conversion and harmonics
CN110096634A (en) * 2019-04-29 2019-08-06 成都理工大学 A kind of house property data vector alignment schemes based on particle group optimizing
WO2021048174A1 (en) * 2019-09-10 2021-03-18 Thales Dis France Sa Method for determining a match between a candidate fingerprint and a reference fingerprint
CN112487867A (en) * 2020-11-03 2021-03-12 杭州电子科技大学 Visual constraint fingerprint identification method based on enhanced triangulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
CHENHAO LIN ET AL.: "Tetrahedron Based Fast 3D Fingerprint Identification Using Colored LEDs Illumination" *

Also Published As

Publication number Publication date
CN114821685B (en) 2023-05-02

Similar Documents

Publication Publication Date Title
CN100559398C (en) Automatic deepness image registration method
Lamdan et al. Object recognition by affine invariant matching
Gökberk et al. 3D shape-based face representation and feature extraction for face recognition
WO2017219391A1 (en) Face recognition system based on three-dimensional data
Si et al. Dense registration of fingerprints
CN101770566B (en) Quick three-dimensional human ear identification method
CN102779269B (en) Human face identification algorithm based on image sensor imaging system
CN109902585B (en) Finger three-mode fusion recognition method based on graph model
CN111428689B (en) Face image feature extraction method based on multi-pool information fusion
CN107958443A (en) A kind of fingerprint image joining method based on crestal line feature and TPS deformation models
Uz et al. Minutiae-based template synthesis and matching for fingerprint authentication
CN110516533A (en) A kind of pedestrian based on depth measure discrimination method again
CN108052912A (en) A kind of three-dimensional face image recognition methods based on square Fourier descriptor
CN101882150A (en) Three-dimensional model comparison and search method based on nuclear density estimation
CN112084895B (en) Pedestrian re-identification method based on deep learning
Yu et al. Sparse ICP with resampling and denoising for 3D face verification
Smeets et al. Fusion of an isometric deformation modeling approach using spectral decomposition and a region-based approach using ICP for expression-invariant 3D face recognition
CN109887009B (en) Point cloud local matching method
Lin et al. Local feature tensor based deep learning for 3d face recognition
JP2006285627A (en) Device and method for retrieving similarity of three-dimensional model
Rattani et al. Robust feature-level multibiometric classification
CN115497125B (en) Fingerprint identification method, system, computer equipment and computer readable storage medium
CN116363177A (en) Three-dimensional point cloud registration method based on geometric transformation and Gaussian mixture model
CN114821685B (en) PSO (particle swarm optimization) optimization Delaunay triangulation-based non-contact 3D fingerprint identification method
CN112767462B (en) Point cloud single-point alignment method based on ridge-valley characteristics and depth characteristic descriptors

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant