CA2378345C

CA2378345C - Method and system for fingerprint template matching

Info

Publication number: CA2378345C
Application number: CA2378345A
Authority: CA
Inventors: Anthony P. Russo
Original assignee: Precise Biometrics AB
Current assignee: Precise Biometrics AB
Priority date: 1999-07-15
Filing date: 1999-11-19
Publication date: 2010-07-06
Anticipated expiration: 2019-11-19
Also published as: CN1359504A; US7299360B2; AU1630800A; EP1203344A1; AU768531B2; CA2378345A1; ATE309581T1; DE69928297D1; WO2001006445A1; CN1167031C; EP1203344B1; US6681034B1; US20040128521A1; DE69928297T2

Abstract

A method for matching templates and a system thereof. Each template includes a plurality of data chunks, each data chunk representing a minutia and comprising a location, a minutia angle and a neighborhood. The location is represented by two coordinates. In one embodiment, each coordinate and the minutia angle are quantized. The neighborhood includes positional parameters with respect to a selected minutia for a predetermined number of neighbor minutiae. In one embodiment, a neighborhood boundary is drawn around the selected minutia and neighbor minutiae are selected from the enclosed area. A reference template is compared to a measured template on a chunk-by-chunk basis. A chunk from each of the template is loaded into a random access memory (RAM). The location, minutia angle and neighborhood of the reference data chunk are compared with the location, minutia angle and neighborhood of the measured data chunk, respectively. In one embodiment, the comparison uses straight subtractions. If the differences for all the parameters meet their respective predetermined tolerances, the measured data chunk matches the reference data chunk. If the number of data chunk matches is equal to or is greater than a predetermined data chunk match rate, the measured template matches the reference template. The neighborhoods are compared by comparing each positional parameter. If all the positional parameters match, the neighbors match. If a predetermined number of neighbor matches is met, the neighborhoods match.

Description

~.'e~'6a, ~.s~,.: h' a ...av, APUY~' -t ~~~~ ~~ ~
CVO 01/06445 PCT/tIS99I2T5~
-ru~k~xern lar t~
METHOD AND SYSTEM FOR FINGERPRINT TEMPLATE MATCHING
5 Appendix A, which is:a part of the present disc ure, is ' a listing of software code for embodiments~.of ponents~of ' this'invention, which are described more mpletely below.
A portion of the disclosure o is patent document contains material which is s ect to copyright protection.
The copyright owner ha - o objection to the facsimile reproduction by one of the patent document or the patent disclosure it appears in the Patent and Trademark Office pate files or records, but otherwise reserves all copyright FIELD OF TF~IE INVENTION
This invention relates to fingerprint verification. More particularly, this invention relates to methods for matching fingerprint templates and-structures thereo=.
3ACKGROLJND OF THE INVENTION v ~ ' ' Biometric identification is used to verify the identity of a person by digitally measuring selected features of some physical characteristic and comparing those measurements with those filed for the person in a reference database, or sometimes on a smart card carried by the person. Physical characteristics that are being used include fingerprints, voiceprints, the geometry of the hand, the pattern of blood vessels on the wrist and on the retina of the eye, the topography of the iris of the eye, facial patterns and the dynamics of writing a signature and typing on 'a keyboard.
The fingerprint is one of the most widely used physical characteristics in biometric identification. Fingerprints are utilized as the most reliable means to identify individuals.
because of two outstanding characteristics; namely that they remain unchanged all through life~and differ from individual to individual. Fingerprints consist of raised friction ridges of skin separated by recessed valleys of skin. Fingerprint ~~~ .
CA 02378345 2002-O1-04 ~ -=."..a:.._..._ ..~~.:,..

WO 01!06a45 CA 02378345 2002 O1-04 p~/Ug99I27561 "minutiae" are conventionally defined as ridge endings or ridge bifurcations where a ridge splits into two ridges.
Since the matching of a fingerprint as an image can require large file capacity to store the measured and referenced features and complex computation to match measured features with'reference features, fingerprint identification is carried out using the positional relationship of features or minutiae that are extracted from the fingerprint image. The minutiae representing a fingerprint image are digitized and stored in a digital memory, which may be a read-only memory, magnetic tape, or the Like. A fingerprint digitized in this manner may be compared with reference fingerprints stored in the memol~y.'. , For the comparison to work, the reference ~~ fingerprint and the measured fingerprint must be extracted, '.
characterized, and digitized so that the fingerprint templates contain the same information and have the same format.
One application of fingerprint identification is to use the fingerprint to unlock a smart card which often contains encrypted information. Presently, a PIN (personal 20. identification number) is required to be.entered by the user .
'- before the encrypted information can be extracted f-ram the smart card and used. The use of a PIN has many drawbacks. For example, the accounting or card issuing organization faces significant administrative cost in handling the secret codes and the card alders have t memorize the secret codes.
~Snaert (~0.S a ~Q j Some privacy laws, ~or example, the privacy laws in Europe, require sensitive =nformatior_ such as a reference fingerprint template to be stored on the smart card in a way that it cannot leave the card without first unlocking the card with a PIN. Therefore, for a fingerprint matching scheme to be practical in Europe, the fingerprint template matching algorithm must be executed by a microprocessor in the smart card. Otherwise, the smart card must be unlocked by some other mechanism, such as a PIN, before the reference fingerprint template can be read. The difziculty of executing conventional fingerprint template matching algorithm on a smart card is mainly due to the limited computational capabilities and memory of a conventional smart card. Fo= example, a conventional _.._ -. ~1J7~BCli2~-Q .; t~~'~~' ~~ ~
r ,'~~~' '~. "~w~fi. ~' ~ t '.::,3chsher.~~aku~e~
2a EP-A-0 159 037 describes a fingerprint identification system comprising an input device, a control unit, a work memory, a card reader and a memory card. The memory card stores reference fingerprint information of the user. The ~5 ~ reference .fingerprint .information comprises a plurality' of minutiae. For each minutia the memory card stores its position, its direction and its relationship to adjacent minutiae. When an identification of a user is to be carried out, the input device produces an image of the user's fingerprint. The card reader reads the reference fingerprint information from the user's memory card and loads in into the work memory. The control unit then compares the reference fingerprint information with the fingerprint image produced by the input device.
. , . .. _ ., . . . . _ . : . .
25=1~i 08 2~t)1.
CA 02378345 2002-O1-04 ...

WO 01/06445 ~~~~ PGTILiS99I:751S~
smart card typically has less than 512 bytes of RAM (with 256 bytes being typical) and between 1 Kilobyte and I6 Kilobytes of memory. An 8-bit RISC (reduced instruction set computer) microprocessor has a speed between 1 MegaFiertz to 10 MegaHertz .5 which a quite slow given the magnitude of the computations required to complete a comparison between a measured ..' .
fingerprint aad a~reference fingerprint. ~In effect, the.
hardware constraints prevent the use of fingerprint to unlock data from smart cards.
. In addition to hardware constraints, another important .
design criterion is cost. U.S. Patent No. 4,582,985 (hereinafter, the '985 patent) entitled "Data Carrier", issued April 15;:?986, to Bo Lofberg and hereby incorporated by '. reference in'its entirety, describes a data carrier of a card ~ .
type. According to the '985 patent, the data carrier includes a fingerprint verificaf;ion device for carrying out the verification process. The verification device includes.a sensor device for sensing a finger tip of the owner and obtaining the corresponding finger print line information. A
specialized smart.card must be used to accommodate the.sensor ' device since a conventional smart.card does not provide such_~ .~
accommodations. In addition, since the fingerprint template generation, storage, and comparison are done at the data carrier, the microprocessor of a conventional smart card is not adequate. Hence, a specialized smart card having a more capable processor must be used, thereby increasing cost.
Therefore, what is needed are systems and fingerprint template matching algorithms that can be executed by a microprocessor with low memory and low computational capacities, trereby keeping the cost of the smart card at an acceptable level.
SUNIMARY~ OF THS INVENTION
verification system and a fingerprint temp ching algorithm are provided. Th _ rprint template matching algorithm is _ a of being executed by a microprocessor with <~~se'~t p4~es30.,3b~
_;_ 2fi-CA 02378345 2002-O1-04 -.

~1~2~, ~.~~t~ DE~~s~ f~ 9~g~~
~. <3.~~:,~x~u.
3a According to a first aspect, the present invention relates to a method for matching templates, comprising:
providing a reference template, said reference template comprising a plurality o reference data chunks;
characterised by providing a measured template_to be compared, said measured~template comprising a plurality of measured data chunks; loading one of said reference data chunks and one of said measured data chunks into a random access memory (RAM) for comparison, said RAM
storing only one reference data chunk and one measured data chunk during said comparison; and comparing said one of said reference data chunk with said one of said measured data chunks.
The fingerprint template matching algorithm is capable of being executed by a microprocessor with low memory and low computational capacities.
According to a second aspect the present invention relates to a fingerprint verification system, characterised by a smart card reader, said smart card reader comprising a fingerprint sensor; and a first microprocessor, said first microprocessor generating a measured template having a plurality of measured data 2S chunks from data read by said fingerprint sensor; a smart card, said smart card comprising: a static memory, said static memory storing a reference template having a plurality of reference data chunks; a second~micro- ' processor,~said second microprocessor executing a matching algorithm for determining whether said measured .
template matches said reference template; and a random access memory, said random access memory (RAM) storing at least one reference data chunk and at least one measured data chunk during execution of said matching algorithm;
3S and a communication channel between said smart card and said smart card reader.
'27.. _ . . . . . _~; _ CA 02378345 2002-O1-04 _ _ ...., .r~"~~.."~.

Phi ~e_ :2Qi~.8 2aQ~~~ DESK ;~w~~ 99~5~0~ ~~~ , 3b According to a third aspect, the present invention relates to a fingerprint verification system, comprising a fingerprint sensor; a static memory, said static memory storing a reference template having a plurality of reference data chunks; a~microprocessor, said micro- ' processor generating a measured template from data read by said fingerprint sensor and executing a matching algorithm for determining whether said measured template matches said reference template, characterised in that said measured template having a plurality of measured data chunks and by a random access memory, said random access memory (RAM) storing one of said reference data chunks and one of said measured data chunks during execution of said matching algorithm.
. . _ . _ . ._. ~ . _ ..._ .. :;~,.. . ~ .i .: ~ w. :~
28 .. . ,1 ~-08 2' _.~....~ .,~....M~.k»~

In one embodiment of this invention, a reference fingerprint template and a measured fingerprint template are generated from a reference fingerprint image and a fingerprint image to be verified, respectively. Each template comprises a plurality of data chunks, each data chunk representing a minutia and comprising a location, a minutia angle and a neighborhood. In one embodiment, the location is represented by two coordinates (x~, y~), the coordinates having a center point (0,0) at the upper left hand corner. The minutia angle 6~ is the angle between the x-axis and a line tangential to the ridge line at the minutia. In one embodiment, each coordinate and the minutia angle are quantized to a selected number of bits. In general, the amount of quantization is a function of the available memory and the degree of accuracy desired.
The neighborhood is made up of a predetermined number of neighbor minutiae which are selected for every minutia extracted from a fingerprint image. Each neighbor minutia is characterized by three positional parameters with respect to the selected minutia. The positional parameters include a distance and two angles.
In one embodiment, an optional neighborhood boundary is drawn around a randomly-selected minutia. In one embodiment, if the number of neighbor minutiae within the neighborhood boundary is less than the predetermined number, all neighbor minutiae within the neighborhood boundary are selected. In another embodiment, a predetermined number of the closest neighboring minutiae are selected. In yet another embodiment, a predetermined number of neighbor minutiae giving the best spread around the randomly selected minutia are selected. In one embodiment, the minutiae having the farthest distance from each other are selected. In another embodiment, minutiae that are very close, e.g., less than approximately 10 pixels, to each other are not selected. In another embodiment, an arbitrary one of the very close minutiae is selected. In one embodiment, a quadrant is arbitrarily drawn using the randomly-selected minutia as the center point and a predetermined number of minutiae are selected from each quadrant.

An x-axis is drawn in the direction of the line tangential to the ridge at the randomly-selected minutia. A y-axis, perpendicular to the x-axis, is drawn. The x-axis and the y-axis intersect at the randomly-selected minutia. Hence, the randomly-selected minutia now has a position of (0,0) and is designated as the "center minutia." A first line is drawn between the center minutia and one of its neighbor minutiae.
The first line has a distance di which is one of the positional parameters representing the neighbor minutia. An angle ~i between the first line and the x-axis is the second positional parameter representing the neighbor minutia. A second line is then drawn from the neighbor minutia to intersect the x-axis, in the direction of the line tangential to the ridge of neighbor minutia. The angle cpi between the x-axis and the second line is the third positional parameter representing the neighbor minutia.
In one embodiment, each positional parameter is quantized to a selected number of bits. For example, each angle ~i and cpi is quantized to six bits and the distance di is quantized to five bits.
A data chunk from each fingerprint template is loaded into a random access memory (RAM) to be compared. In one embodiment, the data chunks are sorted, but such sorting is optional. In one embodiment, the data chunks are sorted in accordance to their x-coordinates. In another embodiment, the data chunks are sorted in accordance to their y-coordinates.
In one embodiment, the data chunks in the reference template are sorted. In another embodiment, the data chunks in the measured template are sorted. In yet another embodiment, both the measured template and the reference template are sorted.
Each of the characterization parameters, i.e., location, minutia angle, and neighborhood, of the measured data chunk to be compared is compared with its respective counterpart in the reference data chunk. The order of the comparison may be varied. For example, the location parameter in the x-coordinates of the measured minutiae are compared first to the corresponding x-coordinates of the reference minutiae. In another embodiment, location parameter in the y-coordinates of the measured minutiae are compared first to the corresponding y-coordinates of the reference minutiae. In one embodiment, the comparison is straight subtraction. If the difference for a parameter pair is equal to or is less than a respective predetermined tolerance, the parameters match. If all of the parameters match, the data chunks match. On the other hand, if one of the parameters fail to match, the data chunks do not match and the next set of data chunks are compared.
The data chunk comparison begins with the first reference data chunk and the first measured data chunk and terminates if the predetermined number of data chunk matches have been reached (i.e., the templates match) or if all of the data chunk combinations have been compared. In one embodiment, if the first reference data chunk and the first measured data chunk do not match, the next measured data chunk is loaded into the RAM
to be compared with the first reference data chunk. The process continues until all the measured data chunks in the measured fingerprint template have been compared. The next reference data chunk and the first measured data chunk are then loaded into the RAM to be compared. The process continues until all the reference data chunks have been compared. In another embodiment, if the first reference data chunk does not match the first measured data chunk, the second reference data chunk is loaded into the RAM to be compared with the first measured data chunk. The process continues until all of the reference chunks have been compared. The second measured data chunk and the first reference data chunk are then loaded into the RAM to be compared. The process continues until all of the measured data chunks have been compared.
In one embodiment, where the data chunks are sorted, the comparison between a reference data chunk and the remaining uncompared measured data chunks is terminated if the difference between the reference data chunk and the measured data chunk exceeds the tolerance. In one embodiment, the next comparison starts with the next reference data chunk and the first matching measured data chunk.
In the alternative embodiment, the comparison between a measured data chunk and the remaining uncompared reference data chunks is terminated if the difference between the measured data chunk and the reference data chunk exceeds the tolerance.
In one embodiment, the next comparison starts with the next reference data chunk and the first matching measured data chunk.
In one embodiment, the neighbor minutiae in a neighborhood are sorted by distance di. The neighborhood comparison starts with the first neighbor in the reference neighborhood and the first neighbor in the measured neighborhood and terminates when the predetermined neighbor matches are reached (i.e., the neighborhoods match) or when all of the neighbor combinations have been compared. Each reference's positional parameter is compared with its corresponding measured positional parameter.
In one embodiment, the comparison is done with straight subtraction. If the difference meets a predetermined tolerance, the parameters match. If all of the positional parameters in the neighbor minutia match, the neighbor minutiae match. If one of the positional parameters does not match, the neighbor minutiae do not match.
In one embodiment, the reference fingerprint template is stored in a static memory in a smart card and the fingerprint template matching algorithm is executed in a microprocessor in the smart card. In one embodiment, the reference fingerprint template is generated and stored at a smart card reader.
The above-described fingerprint template matching algorithm may be executed by a microprocessor with low memory and low computational capacities because the matching algorithm does not require computational intensive operations and does not require comparisons that exceeds 8 bits.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a smart card system having a smart card reader for reading a smart card.
FIG. 2A shows a reference fingerprint template divided into data chunks.
FIG. 2B shows a measured fingerprint template to be compared divided into data chunks.

FIG. 3 shows how a minutia is represented.
FIG. 4A shows a circular neighborhood boundary.
FIG. 4B shows a square neighborhood boundary.
FIG. 4C shows a method of neighbor selection.
S FIG. 4D shows an angle to be represented in complex domain.
FIG. 4E shows an embodiment of a method of neighbor selection.
FIG. S shows construction of positional parameters for a neighbor minutia.
FIG. 6 shows two neighborhoods with errors between them.
FIG. 7A shows a reference data chunk comprising location, minutia angle and neighborhood parameters.
FIG. 7B shows a measured data chunk comprising location, minutia angle and neighborhood parameters.
FIG. 7C shows a reference neighborhood having various positional parameters.
FIG. 7D shows a measured neighborhood having various positional parameters.
FIG. 8A shows a flowchart of a method of template comparison.
FIG. 8B shows a flowchart of another method of template comparison.
FIG. 9A shows a flowchart of a method of detailed data chunk comparison.
FIG. 9B shows a flowchart of a method of detailed data chunk comparison.
FIG. 9C shows a flowchart of a method of detailed data chunk comparison.
FIG. 10 shows a flowchart of a method of detailed neighborhood comparison.
FIG. 11A shows a flowchart of a method of detailed neighbor comparison.
FIG. 11B shows a flowchart of a method of detailed neighbor comparison, with the neighbors sorted.
FIG. 12 shows a flowchart of a method of fingerprint template comparison in a multiple user situation.
_g_ . r ~ :~r~tc~2Q~ ~:~QQ1'~ flES~~'s~
I
~l'O Ql/~6445 CA 02378345 2002-O1-04 p~~'599127561 Use of the~same reference numbers in different figure' indicates similar or Like elements.
DETAILED DESCRIPTION
p ~ . The following description is mean. to~be illustrative only . ' _ ' .
and not limiting. Other embodiments of this invention will be ' obvious in view of the following description to those skilled in the semiconductor processing arts.
FIG. 1 shows one implementation of the present invention.
FIG. 1 shows a smart card reader system having a smart card 100 and a smart card reader 110. Smart ca=d 100 comprises a static memory 102 forwstoring a reference fingerprint.template, a microprocessor 104 for execL_ing a fingerprint template . .
matching algorithm, and RAM (random access msmory~ 106. Smart j, 1~ card reader 110 has a RAM 112 for storing a measured 'j.
fingerprint template, a microprocessor 114, and a fingerprint sensor 116. In one embodiment, fingerprint sensor 116 is a r . capacitive fingerprint sensor, ~ ~~° y'-~a ~eQ~r~.~o~ ;" ~T a ._ ' ~ _ v,~ .. ' .
~ re~ernnno i n i t CTf1'; Yoh Smart card reader 110 can be a stand-alone unit at a point-of-sale or ' #.
can be embedded in an ATM, for example. .. !
I:
TAhen a user wishes to gain access to smart card 100, 2~ he/she inserts smart card 100 into smart card reader 110.
Smart card reader 110 then prompts the use. to place his/her finger on fingerprint sensor 116. The prompt can be either visual or audible. Microprocessor 1.4 on smart card reader 110 then. digitizes the electrical signal measured at fingerprint sensor 116, creates a fingerprint image, extracts minutiae from the fingerprint image, characterizes each minutia, creates a neigaborhood for ,each minutia, and generates a measured fingerprint template prom the characterization parameters and the neiczborhood. The meassred ~inge=-p=in~

3~ template is then stored in RAM 112. Template construction is discussed in detail below.
Minutiae can be extracted from a fingerprint using any known methods. For example, ridge skeletonization technique 29 ,..~ , .. . - . ~ 6 ~D~ X1701 . . .. . _F.a: t k __ r can be used. This technique uses an iterative process to thin a ridge line by discarding pixels that do not affect the continuity of the ridge line until a single pixel line is remaining, the end point of which is the minutia. Other types of minutia extraction technique are described in, for example, U.S. Patent No. 4,646,352 entitled "Method and Device for Matching Fingerprints with Precise Minutia Pairs Selected From Coarse Pairs" issued February 24, 1987 to Asai et al. and U.S.
Patent No. 4,135,147 entitled "Minutiae Pattern Matcher" issued January 16, 1979 to Riganati et al. The minutia extraction technique is not described in detail in this application because it is not part of the invention.
Smart card reader 110 then sends the measured fingerprint template to smart card 100 through a communication channel 120 at a selected speed, typically 9.6 Kilobits/second, or 9600 baud. Microprocessor 104 in smart card 100 receives the measured fingerprint template from smart card reader 110, reads the reference fingerprint template from static memory 102 into RAM 106 chunk by chunk as needed, and executes a fingerprint template matching algorithm to compare the two fingerprint templates.
A collection of data chunks representing a fingerprint, either reference or measured is referred to herein as a fingerprint template. A data chunk consists information of an extracted minutia it corresponds to, including location, minutia angle and neighborhood, which will be discussed below.
FIGS. 2A, 2B, 3, 4A-4E, 5, 6 and 7A-7D show in detail how a fingerprint template is constructed. FIG. 2A shows a reference fingerprint template 200 having a plurality of data chunks 201, 202, 203, etc. Typically, reference fingerprint template 200 comprises 50 chunks, i.e., typically, 50 minutiae are extracted from a fingerprint image. In one embodiment, the minimum number of minutiae required is 10. The minimum number of minutiae is generally determined by the required matching accuracy. In one embodiment, the maximum number of minutiae allowed is 128. Number 128 is chosen as the maximum number of minutiae allowed because with a typical number of 50, 128 gives a sufficient range of the number of minutiae that may be extracted from a fingerprint image. Of course, the minimum and the maximum numbers of minutiae extracted may be modified to suit a particular application. Reference fingerprint template 200 is stored, for example, in static memory 102 on smart card 100.
FIG. 2B shows a measured fingerprint template 210 having a plurality of data chunks 211, 212, 213, etc. Similarly, measured fingerprint template 210 typically comprises 50 data chunks. Measured fingerprint template 210 is generated, for example, by smart card reader 110 and is stored in RAM 112 in smart card reader 110. The number of data chunks extracted from a measured fingerprint may be different from the number of data chunks extracted from a reference fingerprint due to possible false minutiae. In addition, the number of data chunks extracted from a measured fingerprint may be different from the number of data chunks extracted from a reference fingerprint because only a partial fingerprint is read.
However, a match can nevertheless be had if the number of fingerprint data chunk matches satisfies a predetermined number of fingerprint data chunk matches. Therefore, the number of fingerprint data chunks in the measured fingerprint template does not need to be the same as the number of fingerprint data chunks in the reference fingerprint template.
In one embodiment, all the data chunks in both reference fingerprint template 200 and measured fingerprint template 210 contain the same number of bits. A data chunk is constructed for each and every minutia extracted from a fingerprint image.
Therefore, the number of data chunks making up each fingerprint template is equal to the number of minutiae extracted from a fingerprint image. For example, if eighty minutiae are extracted from a fingerprint image, the fingerprint template will contain eighty data chunks.
FIG. 3 shows how each minutia is represented. An x-axis and a y-axis are constructed so that the point (0,0) is at the upper left hand corner of a fingerprint image. The x-axis and the y-axis are drawn regardless of the rotation or the translation of the fingerprint, as will be discussed later.
Each minutia in a fingerprint image is characterized by location (x~,y~) and a minutia angle 6~, where j is the number of minutiae extracted from a fingerprint image. Location (x~,y~) is the location of minutia j with respect to x- and y-axis. Minutia angle 0~ is the angle between the x-axis and a line tangential to a ridge 301 where minutia j is extracted.
Of importance, the minutiae of both the reference fingerprint and the measured fingerprint must contain the same characterization parameters.
A neighborhood NHS is constructed for each minutia j, as described below. FIGS. 4A-4E illustrate various ways of selecting neighbor minutiae in constructing a neighborhood NHS
for minutia j. A neighborhood is constructed for each and every minutia extracted from a fingerprint image. In one embodiment, a neighborhood boundary is drawn around minutia j.
For example, a circle 300 of a fixed radius r, e.g., 90 pixels, is drawn using minutia j as the center point, as shown in FIG.
4A, to provide a neighborhood boundary. In an alternative embodiment, a box 310, representing the neighborhood boundary, is drawn using minutia j as the center point, as shown in FIG.
4B. Of course, any predetermined neighborhood boundary of any shape or size can be drawn.
A predetermined number of neighbor minutiae are selected from the enclosed area. In one embodiment, fifteen (15) minutiae (excluding minutia j) within the neighborhood boundary are selected as the neighbor minutiae for neighborhood NHS. If the enclosed area contains less than the predetermined number of neighbor minutiae, all of the neighbor minutiae contained in the enclosed area are selected. It is noted that neighborhood NHS contains multiple neighbor minutiae N1-Ni where i is a predetermined number, so that there is enough information to make the neighborhood NHS unique enough to meet the accuracy requirement. The number of the neighbors, in general, depends on the uniqueness of the neighborhood NHS desired. For example, if only four neighbors N1-N4 are selected for a neighborhood, the probability of having equivalent or matching neighborhoods is much higher than the corresponding probability associated with a neighborhood of fifteen neighbors.
Similarly, a neighborhood having thirty neighbors would have a lower probability of being equivalent to other neighborhoods.
The uniqueness of the neighborhood determines how often a good fingerprint is rejected (i.e., false rejection rate) and how often a bad fingerprint is accepted (i.e., false acceptance rate). The number of the neighbor minutiae also depends on the acceptable computational complexity or speed that can be satisfactorily performed by a microprocessor.
In one embodiment, the neighbor minutiae are selected such as to maximize the spread around the center point, i.e., minutia j. Referring to FIG. 4C, an x-axis is drawn in the direction of the line tangential to a ridge line at minutia j.
A line is drawn from minutia j to a minutia N1, forming an angle al between the x-axis and the line. Similarly, a line is drawn from minutia j to each of the remaining minutiae within the enclosed area, each forming an angle ak with the x-axis, where k is the number of the minutiae within the enclosed area.
A mean is calculated from all the angles 6k.
In one embodiment, the mean is calculated in a complex domain since mean calculation in a degree domain would be more complicated due to the cyclical nature of the angles. A
complex number C can be described as a sum of its real and imaginary parts a and b, respectively, i.e., C = a + ib where i is the square root of -1.
Referring to FIG. 4D, angles can be expressed in terms of complex numbers by the following equality:
Ck = cos (6k) - i sin (ak) where 6k = tan-1 (bk/ax) Thus, the real part of angle 6k is ak = cos (6k) and the imaginary part of angle 6k is bk = sin(ak). To get an average of all the angles 6k, each angle ak is converted to its complex equivalent. The real part ak and the imaginary part bk are then averaged separately. The averaged real part aaVg and the averaged imaginary part bag are then converted back to a real ~k(ak) number representation 6a~g, for example, aaVg - , k bang - k~k) , and 6aVg - tan ~~avg / aavg) ' k In one embodiment, each angle 6k is examined to determine how close it is to the average angle aaVg by subtracting angle ak from average angle aaVg to get a difference akaiff. A
predetermined number of minutiae that have the largest difference akaiff are selected so that the entire neighborhood has maximum standard angular deviation. It is noted that if there is a cluster of minutiae such as cluster 301 shown in FIG. 4C, the average would be closer to cluster 301. Thus the minutiae that are within cluster 301 would be closer to the mean. In one embodiment, if the number of minutiae within the boundary is more than the predetermined number of neighbor minutiae, the minutiae within cluster 301 would not be selected. If one or more minutiae must be selected from cluster 301 to meet the predetermined number of neighbor minutiae, the minutiae are selected randomly. In other words, the minutiae within cluster 301 are less likely to be selected than the minutiae outside of cluster 301.
In one embodiment, if multiple minutiae, such as minutiae NS and N6, have the same angles, the distances between minutia NS and minutia j and between minutia N6 and minutia j are used for a tie breaker, e.g., the shorter the distance, the more likely the minutia would be selected. In another embodiment, if a minutia, i.e., minutia N3, is too close to minutia j, e.g., 10 pixels, the minutia is not selected. In yet another embodiment, the neighbor minutiae that are furthest away (in distance) from each other are selected. In one embodiment, if two or more neighbor minutiae are very close to each other, e.g., 10 pixels, such as neighbor minutia NS and neighbor minutia N~, either both neighbor minutiae NS and N~ are eliminated or one of neighbor minutiae NS and N~ is arbitrarily selected.
The neighborhood boundary is optional in selecting neighbor minutiae for a neighborhood. For example, the neighbor minutiae are selected such as to maximize the spread around the center point, i.e., minutia j, as described above.
However, in this embodiment, all of the minutiae extracted from the fingerprint are used to calculate the mean. All the other embodiments described above are also applicable to neighbor minutiae selection without a boundary. In an additional embodiment, a quadrant is arbitrarily drawn with minutia j as the center point and a predetermined number of neighbor minutiae are then selected from each quadrant, as shown in Fig.
4E. In general, any neighbor minutiae extraction method may be used to extract the predetermined number i of neighbor minutiae for constructing a neighborhood NHS.
After the neighbor minutiae are selected, a neighborhood NHS is constructed by deriving three positional parameters for each neighbor minutia. The positional parameters are derived in reference to minutia j. These positional parameters include a distance di and two angles ~i and cpi, where i is the number of neighbor minutiae.
FIG. 5 shows how the positional parameters are derived.
An x-axis is drawn in the direction of the line tangential to a ridge line at minutia j. A y-axis perpendicular to the x-axis is then drawn, intersecting the x-axis at minutia j. Hence, minutia j has a position of (0,0) and will be designated as the "center minutia No." A line 500 is drawn from center minutia No to a neighbor minutia Ni. Line 500 has a distance di. An angle is created between the x-axis and line 500. Another line 501 is drawn from neighbor Ni in the direction of a line tangential to a ridge line for neighbor minutia Ni to intersect and extend through the x-axis. The angle between the extended line 501 and the x-axis is angle cpi. Hence, neighbor minutia Ni may be represented by (di, Vii, cpi) . Similarly, neighbor minutia N1 is at a distance dl from center minutia No . An angle ~1 is created between line 502 and the x-axis. Angle cpl is created between line 503 and the x-axis. Neighbor minutia N2 is at a distance d2 from center minutia No. An angle ~2 is created between line 504 and the x-axis. An angle cp2 is created between line 505 and the x-axis. The process is repeated for all i neighbor minutiae of center minutia No. Hence, for a neighborhood having 15 neighbors (i.e., i=15) each of the 15 neighbors N1, N2, N3, . . . , N15 is represented by (di, Vii, cpi) .
Importantly, distance di and angles ~i and cpi are independent of any rotation or translation of a fingerprint image because these are relative positions with respect to center minutia No.

Hence, when a fingerprint image is rotated or translated, all the minutiae associated with that center minutia No are translated and rotated in the same fashion and amount and their relative positions with respect to the center minutia No remain the same.
In one embodiment, each of the parameters di, Vii, and are quantized to a selected number of bits. Quantization allows the reduction of the number of bits to represent a value. More specifically, during quantization, all the numbers that fall between two boundary values are discretized to a quantization value between the boundary numbers.
There are several approaches to derive the quantization value. In linear quantization, a quantization value can be derived by, e.g., a practical approach or a centroid approach.
In the practical approach, the quantization value can be calculated by rounding the quotient of an equation, e.g., x/N, where x is the unquantized value and N is the quantization interval, to the closest integer number. For example, if the unquantized value is 5 and the quantization interval is 2, the quantization value is obtained by rounding 5/2 = 2.5, or 3.
Similarly, if the unquantized value is 5 and the quantization interval is 4, then quantization value is obtained by rounding 5/4 - 1.25 or 1. In a centroid approach, the quantization value is computed as jx p(x) dx C= i , f p (x) dx where C is the quantization value; i and j are the boundary numbers of the interval to be quantized; and p is the probability of occurrence of each number x.
In non-linear quantization, a look-up table may be created to specify the quantization value for each unquantized value.
For example, the lookup table may specify for unquantized values from 0 to 1, the quantization value is 1; for unquantized valued from 1 to 3, the quantization value is 2;
for unquantized values from 3 to 6, the quantization value is 5; and so on. Non-linear quantization has the advantage of setting the quantization value in accordance to the errors associated with each minutia. For example, for a given minutia, the closer the neighbor minutia is to the given minutia, the smaller the error. Similarly, the further the neighbor minutia, the larger the error. Therefore, for neighbor minutiae closer to the given minutia, the parameters associated with the neighbor minutia needs to be quantized more finely than the neighbor minutiae that are further away.
Although any amount of quantization may be used, the amount of quantization depends on the preciseness of the match desired and the memory space available. As is well known in the art, the more the quantization, the less the preciseness of the representation and the less the storage space is required.
For example, for a set of unquantized values 1 through 10, 4 bits are required to represent each number. If the quantization interval is 2, then the largest quantization value needs to be represented is 10/2 or 5. Hence, all 10 numbers can be represented by 3 bits, requiring less storage space.
Similarly, if the quantization interval is 5, the largest quantization value needs to be represented is 10/5 or 2.
Hence, unquantized values 1 through 10 can now be represented by 2 bits, requiring even less storage space. However, the largest error in the first example would be 1 and the largest error in the second example would be 2.5. Hence, the more the quantization, the less the storage requirement and the larger the error.
In one embodiment, distance di is quantized to five bits.
In another embodiment, angles ~i and cpi are quantized to six bits each.
FIG. 6 illustrates two neighborhoods with slight errors, both rotational and translational. Reference neighborhood 510 is constructed with a center minutia No and two neighbor minutiae N1 and N2. Measured neighborhood 512 is constructed with a center minutia No' and two neighbor minutiae N1' and NZ'.
As can be seen, reference neighborhood 510 and measured neighborhood 512 show some resemblance, i.e., neighbor minutiae N1 and NZ and N1' and NZ', respectively, have similar positional relationships with respect to their center minutiae No and No', respectively. However, there are errors between the two neighborhoods. For example, there are slight differences in the distances as well as the angles. The slight differences may be due to the fact that a finger has a fluid surface, thus the surface of the finger stretches and varies under different conditions. Hence, these two slightly varied neighborhoods may be the fingerprint from the same finger, but created at different times. The fingerprint matching algorithm, therefore, must take these slight variations into account because if a perfect match is required, the numerous erroneous rejections will render fingerprint identification impractical.
Because matching tolerances must be accepted, the errors caused by a certain amount of quantization becomes unimportant in the resulting accuracy of the comparison. For example, if the tolerance of distance di is set to 1 and the maximum error caused by quantization of distance di is also 1, if both the reference distance and the measured distance falls into the same quanta, the distances match even though the actual distances may be different. The acceptable tolerance, therefore, determines how much quantization can be done for each parameter.
Another advantage for quantization is to increase the data transmission speed because the size of the data packets becomes smaller after quantization. In one embodiment, the quantized values are compressed using well-known compression method such as, but not limited to, run-length encoding or L-Z-W (Lemtel-Ziv-Welch) compression. In the run-length encoding, for example, if there are five number zero in a row, instead of storing the numbers as "0 0 0 0 0," the numbers are stored as "5 0." In the L-Z-W compression, a dictionary of frequently repeated groups of (8-bit) characters is built on a per-file basis in order to represents frequent character-strings as shorter bit patterns. Hence, compression allows the data to be represented by fewer bits, thus allowing the data transmission to be faster and data storage to be smaller. In one embodiment, the compressed values are transmitted from smart card reader 110 into smart card 100 and are decompressed at smart card 100.

In one embodiment, the location information, i.e., (x~,y~) and the minutia angle 0~ are quantized to a selected number of bits. For example, each of the position measurements x~ and y~
is quantized to four bits and the minutia angle 0~ is quantized to six bits.
The size of a data chunk is determined by the size of the parameters representing the location of the minutia, the minutia angle, and its neighborhood. In one embodiment, each data chunk in the fingerprint template is 34 bytes in size. In this embodiment, the fingerprint template is made of 1-byte location information (i.e., two 4-bit measurements x~ and y~), 6-bit minutia angle 0~ information, and fifteen 17-bit (5 bits for di, 6 bits each for angles ~i and cpi) neighborhood inf ormat ion .
As will be seen later, the fingerprint template matching algorithm compares the fingerprint templates on a chunk-by-chunk basis. Hence, the fingerprint template matching algorithm requires RAM 106 to be of a size large enough to store at least two data chunks, e.g., 68 bytes of data. In a conventional smart card where the typical RAM is 256 bytes in size, no additional memory is required because the required memory in accordance with this invention is within, in fact, well under, the memory confinement of a conventional smart card.
FIG. 7A shows reference data chunk 201 of reference fingerprint template 200 in detail. Reference data chunk 201 includes three parameters: location 702, minutia angle 704, and neighborhood 706. Location 702 includes x- and y- coordinates.
Neighborhood 706, in one embodiment, includes 15 neighbor minutiae N1 through N15. Each neighbor N1 through N15 further includes positional parameters such as distance di between the center minutia No and the neighbor minutia Ni and two related angles Vii, cpi, where i is 1 through 15, as shown in FIG. 7C.
Similarly, FIG. 7B shows data chunk 211 of measured fingerprint template 210 in detail. Measured data chunk 211 includes three parameters: location 752, minutia angle 754, and neighborhood 756. Location 752 includes x- and y-coordinates x~' and y~'. Neighborhood 756, in one embodiment, includes 15 neighbors Nl' through Nls' . Each neighbor Nl' through Nls' further includes positional parameters such as distance di' between the center minutia No' and a neighbor Ni' and two related angles Vii', cpi', where i is 1 through 15, as shown in FIG. 7D. Although FIGs. 7A and 7B show the same number of neighbors in reference fingerprint data chunk 201 and measured fingerprint data chunk 211, it is not necessary because a neighborhood match may be found if the number of neighbor minutia matches satisfies the predetermined number of neighbor minutia matches. For example, for a certain neighborhood boundary drawn, the enclosed area at the reference fingerprint contains enough minutiae to select the predetermined number of neighbor minutiae, e.g., 15 minutiae. However, the enclosed area at the measured fingerprint only contains less than the predetermined number of neighbor minutiae, e.g., 12 minutiae.
If the predetermined neighbor minutia match is 7 and there are 7 matches between the measured neighbor minutiae and the reference neighbor minutiae, there is a match between the reference neighborhood and the measured neighborhood even though neighborhood 706 contains different number of neighbors from neighborhood 756.
FIG. 8A illustrates the fingerprint template matching process in a smart card application. It is noted that although a fingerprint template matching algorithm is described with respect to a smart card and a smart card reader, the fingerprint template matching algorithm described hereinafter can also be executed by, for example, a personal computer, a microprocessor system, or a palm computer.
The template matching process starts with step 600. In one embodiment, the reference data chunks and the measured data chunks are sorted (step 601), but the sorting is not required.
In one embodiment, the data chunks are sorted according to its location information which is represented by the x- and y-coordinates. For example, the data chunks are sorted according to increasing x-coordinate x~ followed by increasing y-coordinate y~. Of course, any sorting technique can be used to sort the data chunks. In general, minutia angle 0~ is not used for sorting due to the cyclical nature of angle measurements.

In one embodiment, either the reference data chunks or the measured data chunks are sorted. In another embodiment, both the reference data chunks and the measured data chunks are sorted. The data chunks are sorted for efficiency purpose, as will be explained below.
Because the fingerprint templates are compared on a chunk-by-chunk basis and the fingerprint templates are serialized, no random access to the fingerprint templates is needed. In other words, once a data chunk has been used, that data chunk is discarded from RAM 106 and a new data chunk is read in. This allows the storage of only one data chunk from each template in RAM 106 during the comparison process.
The first data chunk of measured fingerprint template 210, e.g., measured chunk 211, is loaded from, for example, RAM 112 in smart card reader 110 into RAM 106 (step 602) through communication channel 120 (also refer to FIG. 1). Also in step 602, the first data chunk of reference fingerprint template 200, e.g., reference chunk 201, is loaded from static memory 102 into RAM 106 of smart card 100. Microprocessor 104 compares the two data chunks, i.e. measured chunk 211 and reference chunk 201 in RAM 106 (step 604). How the data chunks are compared in step 604 is discussed in detail later with reference to FIGS. 9A, 9B and 9C.
If the two fingerprint template data chunks, e.g., reference chunk 201 and measured chunk 211, match (step 606), a counter is incremented in step 608. If the number of fingerprint data chunk matches is equal to or greater than a predetermined number (step 610), measured fingerprint template 210 is considered to match reference fingerprint template 200 (step 612) and the fingerprint template matching process terminates in step 614. In one embodiment, more than four data chunk matches is considered a match.
If the two fingerprint template data chunks, e.g., reference chunk 201 and measured chunk 211, do not match (step 606) or if the number of data chunk matches is less than the predetermined number (step 610), it is determined whether the measured chunk used for the comparison, e.g., measured chunk 211, is the last measured chunk in measured fingerprint template 200 (step 616).
If the measured chunk, e.g. measured chunk 211, is not the last measured chunk in measured fingerprint template 210, the next measured chunk, e.g., measured chunk 212, in measured fingerprint template 210 is loaded into RAM 106 from RAM 112 via communication channel 120 (step 618). The newly loaded measured chunk, e.g., measured chunk 212, is compared with the reference chunk in RAM 106, e.g., reference chunk 201 (step 604). The process continues until all measured chunks in measured fingerprint template 210 have been compared or until a predetermined number of data chunk matches have been found.
If the measured chunk, e.g., measured chunk 201, is the last reference chunk in measured fingerprint template 210, then it is determined whether the reference chunk in RAM 106 used for comparison, e.g., reference chunk 201, is the last data chunk in reference fingerprint template 210 (step 620).
If the reference chunk, e.g., reference chunk 201, is not the last reference chunk in reference fingerprint template 200, the next reference chunk in reference fingerprint template 200, e.g., reference chunk 202, is loaded, for example, from static memory 102 into RAM 106 (step 622). Also in step 622, the first measured chunk in fingerprint template 210, e.g., measured chunk 211, is loaded into RAM 106. The newly loaded reference chunk, e.g., reference chunk 202, and first measured chunk, e.g., measured chunk 211, are then compared in step 604.
The process is continued until all reference chunks in reference fingerprint template 200 are compared or until the number of data chunk matches exceeds the predetermined number.
If the reference chunk, e.g., reference chunk 201, is the last reference chunk in reference fingerprint template 200, reference fingerprint template 200 and measured fingerprint template 210 do not match (step 624) and the fingerprint template matching process terminates in step 614.
Alternatively, the fingerprint template matching process can be carried out in accordance with FIG. 8B. In FIG. 8B, the optional sorting is eliminated; the reference chunks are rotated first; the measured chunks are then rotated. In particular, the fingerprint template matching process starts at step 650. The first data chunk of measured fingerprint template 210, e.g., measured chunk 211, is loaded from smart card reader 110 into RAM 106 (step 652) through communication channel 120. Also in step 652, the first data chunk of reference fingerprint template 200, e.g., reference chunk 201, is loaded from static memory 102 into RAM 106 of smart card 100. Microprocessor 104 compares the two fingerprint template data chunks, e.g., measured chunk 211 and reference chunk 201 in RAM 106 (step 654). How the data chunks are compared in step 654 is discussed in detail below with reference to FIG.
9C.
If the two data chunks, e.g., chunks 201 and 211, match (step 656), a counter is incremented in step 658. If the number of data chunk matches is equal to or is greater than a predetermined number (step 660), measured fingerprint template 210 matches reference fingerprint template 200 (step 662) and the process ends in step 664.
On the other hand, if the two data chunks, e.g., data chunks 201 and 211 do not match (step 656) or if the number of matches of the data chunks is less than the predetermined number (step 660), it is determined whether the measured chunk used for comparison, e.g., reference chunk 201, is the last reference chunk in reference fingerprint template 200 (step 666) .
If the reference chunk, e.g., reference chunk 201, is not the last data chunk in reference fingerprint template 200, the next reference chunk, e.g., reference chunk 202, in reference fingerprint template 200 is loaded from static memory 102 into RAM 106 (step 668). The newly loaded reference chunk, e.g., reference chunk 202, is compared with the measured chunk in RAM
106, e.g., measured chunk 211, in step 654. The process continues.
If the reference chunk, e.g., reference chunk 201, is the last reference chunk in reference fingerprint template 200, it is determined whether the measured chunk in RAM 106 used for comparison, e.g., measured chunk 211, is the last measured chunk in measured fingerprint template 210 (step 670).

If the measured chunk, e.g., measured chunk 211, is not the last data chunk in measured fingerprint template 210, the next measured chunk in measured fingerprint template 210, e.g., measured chunk 212, is loaded from RAM 112 into RAM 106 of smart card 100 via communication channel 120 (step 672). The first reference chunk in fingerprint template 200, e.g., reference chunk 201, is also loaded into RAM 106 (step 672).
The newly loaded measured chunk, e.g., measured chunk 212, and the first reference chunk, e.g., reference chunk 201, are then compared (step 654).
If the measured chunk, e.g., measured chunk 211, is the last measured chunk in measured fingerprint template 210, measured fingerprint template 210 and reference fingerprint template 200 do not match (step 674) and the fingerprint template matching process terminates at step 664.
FIG. 9A illustrates a flowchart of step 604 of how the fingerprint template data chunks are compared in detail. A
reference chunk includes a positional parameter (x~,y~) where x~
is the x-coordinate and y~ is the y-coordinate of a reference minutia. A measured chunk includes a positional parameter of (x~' ,y~' ) where x~' is the x-coordinate and y~' is the y-coordinate of a measured minutia. In this embodiment, the measured data chunks are sorted in accordance with x-coordinates x~ followed by y-coordinates y~. In one embodiment, the x-coordinates are compared by subtracting measured x-coordinate x~' from reference x-coordinate x~, i.e. x~-x~'=Xjdiff (step 850) . If the absolute value of difference x~diff is greater than the predetermined tolerance x~tolerance (step 852) , then it is determined whether the raw difference x~diff is less than the predetermined tolerance x~toleran~e (step 868) . If difference x~diff is less than x~tolerance, the rest of the measured data chunks will not be compared against the reference data chunk and it is determined whether the reference chunk used is the last reference chunk (step 620 in FIG. 8A) and the process continues. If the difference x~diff is greater than Xjtolerance~
the chunks do not match (step 872) and the process continues at step 606 (FIG. 8A).

If the absolute value of the difference xjdiff is less than or equal to a predetermined tOleranCe Xjtolerance ( step 852 ) , i.e., the x-coordinates match, the y-coordinates are compared.
Measured y-coordinate yj' is subtracted from reference y-coordinate yj , i . a . yj-yj' =yjdiff (step 854 ) . If the absolute value of the difference yjdiff is greater than the predetermined tolerance yjtolerance (step 856) , it is determined whether the raw difference yjdiff is less than yjtolerance. If difference yjdiff 1S
less than yjtolerance. the remaining measured data chunks will not be compared against the reference data chunk and it is determined whether the reference chunk used for comparison is the last reference chunk (step 620 in FIG. 8A). If difference yjdiff 1S greater than yjtolerance~ the chunks do not match (step 872) and the process continues at step 606 (FIG. 8A). If the absolute value of the difference yjdiff 1S less than or equal to a predetermined tolerance yjtoierance ( step 856 ) , i . a . , the y-coordinates match, the minutia angles 0j and 8j' are compared.
It is noted that a straight subtraction in degree domain does not work well with angles due to their cyclical nature.
Measured minutia angle 0j' is subtracted from reference minutia angle 0j , i . a . 8j -6j' =0jdiff ( step 858 ) . If the absolute value of the difference 0jdiff is less than or equal to a predetermined tolerance 0j tolerance (step 860 ) , i . a . , the minutia angles match.
If the absolute value of the difference 0jdiff is greater than the predetermined tolerance 0jtolerance~ it is determined whether the absolute value of the difference 0jdiff is greater than or equal to the difference between the maximum allowable quantized minutia angle ejmax and the predetermined tolerance 0j tolerance (step 861). If value l0jdiffl is greater than or equal to the difference between values 0jm~ arid 0jtolerance~ the minutia angles match.
An example for angle comparison is described below. In the embodiment where the angles are quantized to 6 bits, thus having values 0 to 63, 0 is the minimum value for 0j (0jmin) and 63 is the maximum value for 0j (OjmaX) . The quantized values are then subtracted, e.g. , ~0j-0j' ~=Ajdiff, if minutia angle 0j has a quantized value of 61 while minutia angle 8j' has a quantized value of 63, and the match tolerance is 2. Then ~61-63~=2, which is within the tolerance range. Thus, there is a match between minutia angles 0j and 0j'. In a situation where angle 0j has a quantized value of 0 while minutia angle 0j' has a quantized value of 63, and the match tolerance is 2. Then ~0-63~=63, which is out of the tolerance range. However, due to the cyclical nature, the difference between 63 and 0 is actually only 1, which is within the tolerance. Therefore, in addition to checking whether ~ 0j-0j' ~ ~0jtolerance~ whether ejdiff~ejmax-ejtolerance is also checked. In the above example, ejdiff~ejmax-ejtolerance ( 63>_63 -2 ) is true . Therefore, minutia angles 0j and 0j' match.
If the minutia angles 0j and 8j' match, neighborhood NHj is compared with neighborhood NHj' in step 862, which is described in detail in reference to FIG. 10.
If the neighborhood NHj from the reference chunk matches neighborhood NHj' from the measured chunk (step 864), the data chunks match (step 866). Conversely, if any of the parameters, i.e. x-coordinates, y-coordinates, minutia angles and neighborhoods, fail to match, the data chunks do not match (step 872) . The typical value for xjtolerance and is 4 yjtolerance and the typical value for 0j tolerance is 7. These typical values are derived from experimentation.
Since the measured data chunks are sorted in accordance with the x-coordinates xj, the comparisons for the remaining measured data chunks can be eliminated because the remaining measured data chunks will also have x-coordinates xj' exceeding the tolerance. Hence, sorting cuts down the number of comparisons. The number of comparisons can be further reduced by sorting the data chunks by y-coordinates yj.
A flag can be set so that only a subset of the measured data chunks is used for further comparisons. FIG. 9B
illustrates a flowchart of step 604 of how the fingerprint template data chunks are compared in detail. The x-coordinates are compared by subtracting measured x-coordinate xj' from reference x-coordinate xj , i . a . xj-xj' =Xjdiff (step 880 ) . If the absolute value of the difference x~diff is less than or equal to a predetermined tOleranCe Xjtolerance (step 881) , it is determined whether this is the first match between reference x-coordinate x~ and measured x-coordinate x~' (step 885). If the match is the first match, a flag is set for the measured chunk containing the first matching x-coordinate x~' (step 886).
If the absolute value of the difference x~diff is greater than the predetermined tOleranCe X~tolerance (step 881) , then it is determined whether the difference x~diff is less than the predetermined tolerance x~toleran~e (step 882) . If difference xjdiff 1S less than xjtolerance~ the rest of the measured data chunks will not be compared against the reference data chunk and it is determined whether the reference chunk used is the last reference chunk (step 883). If the reference chunk is not the last chunk in the reference fingerprint template, the next reference chunk and the start chunk is loaded into the RAM for comparison (step 884) and the process continues with step 880.
If the reference chunk is the last chunk in the reference fingerprint template, the fingerprint templates do not match (step 624 in FIG. 8A).
If the match is not the first match (step 885) or after the start chunk has been flagged (step 886), the y-coordinates are compared. Measured y-coordinate y~' is subtracted from reference y-coordinate y~, i.e. y~-y~'=yjdiff (step 887) . If the absolute value of the difference yjdiff is greater than the predetermined tolerance yjtoierance (step 888) , it is determined whether the difference yjdiff 1S less than y~toierance (step 889) .
If difference yjdiff is less than y~toieran~e~ the remaining measured data chunks will not be compared against the reference data chunk and it is determined whether the reference chunk used for comparison is the last reference chunk (step 883). If difference yjdiff 1S greater than yjtolerance~ the chunks do not match (step 896) and the process continues at step 606 (FIG.
8A) .
If the absolute value of the difference yjdiff 1S less than or equal to a predetermined tolerance y~toierance (step 888) , the remaining steps 890 through 896 are similar to steps 858-866 described in reference to FIG: 9A.
As can be seen from FIGs. 9A and 9B, x-coordinates and y-coordinates comparison filters the number of comparisons.
Although the fingerprint template matching algorithm is capable of comparing a fingerprint with a rotated and/or translated same fingerprint, the algorithm can execute faster if the finger is constrained in a certain orientation. By allowing translation and rotation, none, or very few values can be eliminated because the tolerance must be sufficiently large so as not to eliminate everything during the comparison and thus generating a false rejection.
FIG. 9C illustrates a flowchart of step 654 of how the fingerprint template data chunks are compared in detail, where every chunk is compared. Parameters xj and xj' are compared (step 806) , e.g. , by straight subtractions. If value xjdiff, i . a . , ~ xj-xj' ~ =xjdiff, is greater than a tunable, predetermined tolerance x arameter x does not match arameter x '.
jtolerance~ P j p j Therefore, the data chunks do not match (step 824) and the next data chunk is read (refer to FIG. 8B) . When value x~diff is equal to or is less than the tolerance Xjtolerance (step 812) , parameter xj matches parameter xj' (step 808) and parameter yj is compared with parameter yj' in step 810 with a straight subtraction.
If value yjdiff. i . a . ~ yj' -yj I =yjdiff, is greater than a tunable, predetermined tolerance yjtolerance~ parameter yj does not match parameter yj' (step 812) and the data chunks do not match (step 824). The next data chunk is read (see FIG. 8B). When value yjaiff is equal to or is less than the tolerance yjtolerance.
parameter yj matches parameter yj' (step 812) and parameter 0j is compared with parameter 0j' in step 814.
If minutia angle difference 0jdiff 1S greater than a tunable, predetermined minutia angle tolerance Ajtolerance~
parameter 0j does not match parameter 0j' (step 816) and the data chunks do not match (step 824) and the next data chunk is read (see FIGS. 8A and 8B). When minutia angle difference 0jdiff is equal to or is less than the minutia angle tolerance ejtolerance or when difference 0jdiff 1S greater than or equal to the maximum allowable quantized minutia angle minus the predetermined tolerance 0jtolerance~ parameter 0j matches parameter 0j' (step 816) and neighborhood NHj is compared with neighborhood NHj' in step 818, which is described in detail in reference to FIG. 10. If the neighborhood NHj from the reference chunk matches neighborhood NHS' from the measured chunk (step 820), the data chunks match (step 822).
Conversely, if the neighborhood NHS and neighborhood NHS' do not match, the data chunks do not match. The sequence of comparison may vary. For example, parameter y~ maybe compared before comparing parameter x~; angle 0~ maybe compared before comparing parameter y~; and so on.
FIG. 10 shows a flowchart of the neighborhood matching algorithm of steps 862, 893 and 818. In step 900, the first neighbor in an reference chunk (hereinafter reference neighbor), e.g., N1 (FIG. 7A), and the first neighbor in a measured chunk (hereinafter measured neighbor), e.g., N1' (FIG.
7B), are compared. Step 900 is described in detail below.
If the neighbors match (step 902), a neighbor counter is incremented in step 904. If the counter value is equal to or is greater than a predetermined number of neighbor matches (step 906), the neighborhoods match (step 918). For example, 7 neighbor matches for a 15-neighbor neighborhood can be considered a match of the neighborhoods.
If the neighbors do not match (step 902) or if the number of neighbor matches is less than the predetermined number of neighbor matches (step 906), it is determined whether the measured neighbor used in comparison is the last neighbor in the measured chunk (step 908).
If the measured neighbor is not the last measured neighbor in the measured chunk, the next measured neighbor is read (step 910) and the new measured neighbor is compared with the reference neighbor already in use. If the measured neighbor is the last measured neighbor in the measured chunk, it is determined whether the reference neighbor is the last reference neighbor in the reference chunk (step 912). If the reference neighbor is not the last reference neighbor in the reference chunk, the next reference neighbor and the first measured neighbor are read (step 914) and compared (step 900). If, however, the reference neighbor is the last reference neighbor in the reference chunk, the neighborhoods do not match (step 916 ) .

FIG. 11A illustrates step 900 in FIG. 10 in detail and shows a neighbor to neighbor comparison. The distances di of an reference neighbor is compared with a distance di' of a measured neighbor in step 1000. This distance comparison is done with straight subtractions, i . a . ~ di-di' ~ =didiff . If the distance difference diaiff meets the distance tolerance ditolerance (step 1002), angle ~i is compared with angle Vii' (step 1004).
Angle comparison is again done in a similar manner as minutia angle comparison described above. Specifically, angle ~i matches angle ~i' 1f I (~j-~1j' I <~jtolerance arid ~j-~)j' >~jmax-~jtolerance(Step 1006) .
Angles cpi and cpi' are then compared in step 1008.
Similarly, angle cpi matches angle cpi' if ~ cpj -cpj' ~ <(pjtolerance and ' (ste 1010). The nei hbors match (ste ~j ~~jmax-~jtolerance p g P
1012) if all parameters for a neighbor, e.g., distance di and angles ~i and cpj, match. The typical value for ditolerance is 1 and the typical value for ~itolerance arid ~itolerance is 6. These typical values are derived from experimentation.
If any of the parameters, i.e., distance di or angle ~i or cpi, do not match, the neighbors do not match (step 1014).
Again, the sequence of comparison, i . a . , di, ~i and cpi, may be modified.
FIG. 11B shows an alternative embodiment of neighbor comparison process in step 900. In this embodiment, the measured chunks are sorted in accordance with distance di'.
Measured distances di' is subtracted from reference distance di to derive a difference didiff (step 1050) . If the distance difference d;,aiff meets a predetermined distance tolerance ditolerance (step 1052) , angle ~i is compared with angle Vii' (step 1054) . If angles ~i and Vii' match (step 1056) , cpi and cpi' are compared in step 1058. If angles cpi and cpi' match (step 1060) , the neighbors match (step 1062).
If measured distance di does not match reference distance di', it is determined whether the difference didiff 1S less than the predetermined tolerance ditolerance ( step 1053 ) . If the difference dydiff 1S less than the predetermined tolerance d_ itolerance~ the neighbors do not match and the remaining measured neighbors are not compared. The process continues with step 912 in FIG. 10. If any of the parameters, i.e., distance di or angle ~i or cpi, do not match, the neighbors do not match (Step 1064) .
In one embodiment, for multiple users where multiple reference fingerprints are stored in a database, if the measured fingerprint template does not match the first reference fingerprint template, the next reference fingerprint template is compared with the measured fingerprint template, as shown in FIG. 12. Multiple reference fingerprint template comparison starts with step 680. The first reference fingerprint template and the measured fingerprint template are loaded into the RAM (one chunk at a time)(step 682). The fingerprint templates are compared in step 684, in accordance with the template comparison algorithm described above, with reference to FIGs. 8A and 8B. Whether there is a match between the fingerprint templates is determined in step 686. If there is a match (step 688), the process ends in step 696. If there is no match, it is determined whether the last reference fingerprint template has been used (step 690). If the last reference fingerprint template has been used, there is no match (step 692) and the process ends in step 696. However, if the last reference fingerprint template has not been used (step 690), the next reference fingerprint template is loaded to be compared with the measured fingerprint template (step 694). In general, the number of reference templates can be stored is limited by the memory capacity. For example, in a conventional smart card, the typical static memory available is 1 Kilobytes to 16 Kilobytes. Therefore, the number of reference templates can be stored is limited by that memory capacity.
In another embodiment for a multi-user system, the measured fingerprint template is compared with all the reference templates in a similar manner as described above.
However, in this embodiment, for each template comparison, all data chunks in both the measured fingerprint template and the reference template are compared. For example, even when a match is found for a pair of fingerprint templates, e.g., more than four data chunk matches, the remaining combinations of data chunks are compared. The number of data chunk matches for each reference fingerprint template is stored. If the measured fingerprint template matches more than one reference fingerprint templates, then the user producing the reference fingerprint template that gives the highest number of data chunk matches is identified as the user.
The algorithm in accordance with the present invention may be executed in a microprocessor having minimal capabilities, such as an 8-bit microprocessor on a conventional smart card because every parameter used in the algorithm is less than 8 bits. For example, location parameters x~ and y~ are quantized to four bits each; minutia angle 8~ and angles ~i and cpi are quantized to six bits each; and distance di is quantized to five bits. In addition, the fingerprint template matching algorithm in accordance with the present invention does not involve computational intensive processes involving for example, floating points, multiplications, and divisions.
Instead, the algorithm in accordance with the present invention only uses subtractions. Hence, the algorithm can be executed by a weak microprocessor having very low speed, e.g., between 1 and 10 MegaHertz, such as a microprocessor in a conventional smart card.
Although the invention has been described with reference to particular embodiments, the description is only an example of the invention's application and should not be taken as a limitation. Various other adaptations and combinations of features of the embodiments disclosed are within the scope of the invention as defined by the following claims.

Claims

I claim:

1. A method for matching templates, comprising:
providing a reference template, said reference template comprising a plurality of reference data chunks;
providing a template to be compared, said measured template comprising a plurality of measured data chunks;
loading one of said reference data chunks and one of said measured data chunks into a random access memory (RAM) for comparison, said RAM storing only one reference data chunk and one measured data chunk during said comparison; and comparing said one of said reference data chunks with said one of said measured data chunks.

2. The method of Claim 1, wherein providing a reference template comprises creating a reference data chunk for each minutia making up a reference fingerprint image.

3. The method of Claim 1, wherein providing a measured template comprises creating a measured data chunk for each minutia making up a fingerprint image to be verified.

4. The method of Claim 1, wherein each of said reference data chunks and each of said measured data chunks comprise:
a location of a minutia;
a minutia angle of said minutia; and neighborhood information of said minutia.

5. The method of Claim 4, further comprising comparing a location of said reference data chunk to a location of said measured data chunk.

6. The method of Claim 4, further comprising representing said location with a first coordinate and a second coordinate.

7. The method of Claim 6, further comprising quantizing said first and second coordinates.

8. The method of Claim 7, wherein said quantizing said first and second coordinates comprises quantizing each of said first and second coordinates to equal to or less than eight bits.

9. The method of Claim 7 wherein said quantizing is a function of a predetermined size of a memory.

10. The method of Claim 7 wherein said quantizing is a function of a predetermined accuracy.

11. The method of Claim 6, further comprising comparing a first coordinate of said reference data chunk to said a coordinate of said measured data chunk.

12. The method of Claim 11, wherein said comparing comprises subtracting said first coordinate of said measured data chunk from said first coordinate of said reference data chunk to obtain a first coordinate difference.

13. The method of Claim 12, wherein said first coordinate of said measured data chunk matches said first coordinate of said reference data chunk if said first coordinate difference is less than a predetermined first coordinate tolerance.

14. The method of Claim 6, further comprising comparing a second coordinate of said reference data chunk to a second coordinate of said measured data chunk.

15. The method of Claim 14, wherein said comparing comprises subtracting said second coordinate of said measured data chunk from said second coordinate of said reference data chunk to obtain a second coordinate difference.

16. The method of Claim 15, wherein said second coordinate of said measured data chunk matches said second coordinate of said reference data chunk if said second coordinate difference is less than a predetermined second coordinate tolerance.

17. The method of Claim 4, further comprising quantizing said minutia angle.

18. The method of Claim 17, further comprising quantizing said minutia angle to equal to or less than eight bits.

19. The method of Claim 4, further comprising comparing a minutia angle of one of said reference data chunks to a minutia angle of one of said measured data chunks.

20. The method of Claim 19, wherein said comparing comprises subtracting said minutia angle of said measured data chunk with said minutia angle of said reference data chunk to obtain a minutia angle difference.

21. The method of Claim 20, wherein said minutia angle of said measured data chunk matches said minutia angle of said reference data chunk if said minutia angle difference is less than a predetermined minutia angle tolerance.

22. The method of Claim 19, wherein said comparing comprises comparing said minutia angle of said measured data chunk with said minutia angle of said reference data chunk in a complex domain.

23. The method of Claim 4, further comprising selecting a predetermined number of neighbors.

24. The method of Claim 23, further comprising constructing a neighborhood boundary surrounding said minutia.

25. The method of Claim 24, wherein constructing a neighborhood boundary comprises constructing a neighborhood boundary of a predetermined shape and size.

26. The method of Claim 24, wherein said neighborhood boundary is circular having said minutia as a center for said neighborhood boundary.

27. The method of Claim 24, wherein said neighborhood boundary is square having said minutia as a center for said neighborhood boundary.

28. The method of Claim 23, further comprising selecting the closest minutiae to said minutia as said neighbors.

29. The method of Claim 23, further comprising arbitrarily selecting a neighbor from a plurality of closely spaced minutiae.

30. The method of Claim 23, further comprising disposing all closely spaced minutiae.

31. The method of Claim 23, further comprising:
constructing an x-axis having said minutia as a center point and in a direction tangential to a ridge where said minutia is extracted;
constructing a y-axis perpendicular to said x-axis;
constructing a line between said minutia and said neighbor, said line having a distance, said line being at a first angle from said x-axis;
constructing a line from said neighbor in a direction tangential to a ridge where said neighbor is extracted, said line intersecting said x-axis and forming a second angle from said x-axis; and representing said neighbor with said distance, said first angle, and said second angle.

32. The method of Claim 31, further comprising quantizing said distance.

33. The method of Claim 32, further comprising quantizing said distance to equal to or less than eight bits.

34. The method of Claim 31, further comprising quantizing said first angle.

35. The method of Claim 34, further comprising quantizing said first angle to equal to or less than eight bits.

36. The method of Claim 31, further comprising quantizing said second angle.

37. The method of Claim 36, further comprising quantizing said second angle to equal to or less than eight bits.

38. The method of Claim 31, further comprising comparing a distance of a reference neighbor in said reference data chunk to a distance of a neighbor in said measured data chunk.

39. The method of Claim 38, wherein said comparing comprises subtracting said distance of said measured data chunk with said distance of said reference data chunk to obtain a distance difference.

40. The method of Claim 39, wherein said distance of said measured data chunk matches said distance of said reference data chunk if said distance difference is less than a predetermined distance tolerance.

41. The method of Claim 31, further comprising comparing a first angle of a reference neighbor in said reference data chunk to a first angle of a neighbor in said measured data chunk.

42. The method of Claim 41, wherein said comparing comprises subtracting said first angle of said measured data chunk with said first angle of said reference data chunk to obtain a first angle difference.

43. The method of Claim 42, wherein said first angle of said measured data chunk matches said first angle of said reference data chunk if said first angle difference is less than a predetermined first angle tolerance.

44. The method of Claim 41, wherein said comparing comprises comparing said first angle of said measured data chunk with said first angle of said reference data chunk in a complex domain.

45. The method of Claim 31, further comprising comparing a second angle of a reference neighbor in said reference data chunk to a second angle of a neighbor in said measured data chunk.

46. The method of Claim 45, wherein said comparing comprises subtracting said second angle of said measured data chunk with said second angle of said reference data chunk to obtain a second angle difference.

47. The method of Claim 46, wherein said second angle of said measured data chunk matches said second angle of said reference data chunk if said second angle difference is less than a predetermined second angle tolerance.

48. The method of Claim 45, wherein said comparing comprises comparing said second angle of said neighbor in said reference data chunk to said second angle of said neighbor in said measured data chunk in a complex domain.

49. The method of Claim 31, wherein a neighbor of said reference data chunk matches a neighbor of said measured data chunk if a distance, a first angle and a second angle of a reference neighbor in said reference data chunk match a distance, a first angle and a second angle of a neighbor in said measured data chunk.

50. The method of Claim 49, wherein a neighborhood of said reference data chunk matches a neighborhood of said measured data chunk if a number of neighbor matches is greater than a predetermined neighbor matching rate.

51. The method of Claim 23, further comprising calculating a standard deviation for each of said neighbors.

52. The method of Claim 23, further comprising maximizing a spread of said neighbors around said minutia.

53. The method of Claim 4, further comprising comparing a neighborhood of said reference data chunk to a neighborhood of said measured data chunk.

54. The method of Claim 53, wherein said reference data chunk matches said measured data chunk if a location, a minutia angle and a neighborhood of said reference data chunk match a location, a minutia angle and a neighborhood of said measured data chunk.

55. The method of Claim 54, wherein said reference template matches said measured template if a number of data chunk matches is equal to or is greater than a predetermined data chunk match rate.

56. The method of Claim 1, further comprising sorting said reference data chunks.

57. The method of Claim 1, further comprising sorting said measured data chunks.

58. The method of Claim 56, further comprising setting a flag in said measured data chunk indicating a first match between said first coordinate of said reference data chunk and said first coordinate of said measured data chunk.

59. The method of Claim 58, further comprising comparing said second coordinate of said reference data chunk to said second coordinate of said measured data chunk starting from said flagged measured data chunk.

60. The method of Claim 57, further comprising setting a flag in said reference data chunk indicating a first match between said first coordinate of said reference data chunk and said first coordinate of said measured data chunk.

61. The method of Claim 60, further comprising comparing said second coordinate of said reference data chunk to said second coordinate of said measured data chunk starting from said flagged reference data chunk.

62. The method of Claim 20, wherein said minutia angle of said measured data chunk matches said minutia angle of said reference data chunk if said minutia angle difference is greater than or equal to a difference between a maximum allowable minutia angle and a predetermined minutia angle tolerance.

63. A fingerprint verification system, comprising, a smart card reader, said smart card reader comprising:
a fingerprint sensor; and a first microprocessor, said first microprocessor generating a measured template having a plurality of measured data chunks from data read by said fingerprint sensor;

a smart card, said smart card comprising:
a static memory, said static memory storing a reference template having a plurality of reference data chunks;

a second microprocessor, said second microprocessor executing a matching algorithm for determining whether said measured template matches said reference template; and a random access memory;~said random access memory (RAM) storing at least one reference data chunk and at least one measured data chunk during execution of said matching algorithm; and a communication channel between said smart card and said smart card reader.

64. The system of Claim 63, wherein a measured data chunk is loaded into said RAM of said smart card through said communication channel.

65. The system of Claim 63, wherein each of said measured data chunks and each of said reference data chunks comprise:
a location of a minutia;
a minutia angle of said minutia; and, a neighborhood of said minutia.

66. The system of Claim 65, wherein said location comprises:
a first coordinate; and a second coordinate.

67. The system of Claim 66, wherein said first coordinate is quantized.

68. The system of Claim 67, wherein said first coordinate is quantized to equal to or less than eight bits.

69. The system of Claim 66, wherein said second coordinate is quantized.

70. The system of Claim 69, wherein said second coordinate is quantized to equal to or less than eight bits.

71. The system of Claim 65, wherein said minutia angle is quantized.

72. The system of Claim 71, wherein said minutia angle is quantized to equal to or less than eight bits.

73. The system of Claim 65, wherein said neighborhood comprises positional parameters of a plurality of a predetermined number of neighbors.

74. The system of Claim 73, wherein said positional parameters comprises:
a distance between said minutia and a neighbor minutia;
a first angle between a first coordinate in a direction tangential to a ridge where said minutia is extracted and a line drawn between said minutia and said neighbor minutia; and a second angle between said first coordinate and a second coordinate in a direction tangential to a ridge where said neighbor minutia is extracted.

75. The system of Claim 74, wherein said distance is quantized.

76. The system of Claim 75, wherein said distance is quantized to equal to or less than eight bits.

77. The system of Claim 74, wherein said first angle is quantized.

78. The system of Claim 77, wherein said first angle is quantized to equal to or less than eight bits.

79. The system of Claim 74, wherein said second angle is quantized.

80. The system of Claim 79, wherein said second angle is quantized to equal to or less than eight bits.

81. The system of Claim 73, wherein said neighbors are sorted.

82. The system of Claim 63, wherein said second microprocessor is an 8-bit microprocessor.

83. The system of Claim 63, wherein said second microprocessor has a speed between approximately 1 MegaKertz and approximately 10 MegaHertz.

84. The system of Claim 63, wherein said RAM is approximately 256 Kilobyte in size.

85. The system of Claim 63, wherein said reference data chunks are sorted.

86. The system of Claim 63, wherein said measured data chunks are sorted.

87. The system of Claim 63, wherein said RAM stores only one reference data chunk and only one measured data chunk during execution of said matching algorithm.

88. A fingerprint verification system, comprising:
a fingerprint sensor;
a static memory, said static memory storing a reference template having a plurality of reference data chunks;
a microprocessor, said microprocessor generating a measured template from data read by said fingerprint sensor and executing a matching algorithm for determining whether said measured (template matches said reference template, wherein said measured template having a plurality of measured data chunks; and by a random access memory, said random access memory (RAM) storing one of said reference data chunks and one of said measured data chunks during execution of said matching algorithm.

89. The system of Claim 88, wherein each of said measured data chunks and each of said reference data chunks comprise:
a location of a minutia;
a minutia angle of said minutia; and a neighborhood of said minutia.

90. The system of Claim 89, wherein said location comprises:
a first coordinate; and a second coordinate.

91. The system of Claim 90, wherein said first coordinate is quantized.

92. The system of Claim 91, wherein said first coordinate is quantized to equal to or less than eight bits.

93. The system of Claim 90, wherein said second coordinate is quantized.

94. The system of Claim 93, wherein said second coordinate is quantized to equal to or less than eight bits.

95. The system of Claim 89, wherein said minutia angle is quantized.

96. The system of Claim 95, wherein said minutia angle is quantized to equal to or less than eight bits.

97. The system of Claim 89, wherein said neighborhood comprises positional parameters of a plurality of a predetermined number of neighbors.

98. The system of Claim 97, wherein said positional parameters comprise:

a distance between said minutia and a neighbor minutia;
a first angle between a first coordinate in a direction tangential to a ridge where said minutia is extracted and a line drawn between said minutia and said neighbor minutia; and a second angle between said first coordinate and a second coordinate in a direction tangential to a ridge where said neighbor minutia is extracted.

99. The system of Claim 98, wherein said distance is quantized.

100. The system of Claim 99, wherein said distance is quantized to equal to or less than eight bits.

101. The system of Claim 98, wherein said first angle is quantized.

102. The system of Claim 101, wherein said first angle is quantized to equal to or less than eight bits.

103. The system of Claim 98, wherein said second angle is quantized.

104. The system of Claim 103, wherein said second angle is quantized to equal to or less than eight bits.

105. The system of Claim 97, wherein said neighbors are sorted.

106. The system of Claim 88, wherein said reference data chunks are sorted.

107. The system of Claim 88, wherein said measured data chunks are sorted.