CN102332156A - Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering - Google Patents

Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering Download PDF

Info

Publication number
CN102332156A
CN102332156A CN201110264780A CN201110264780A CN102332156A CN 102332156 A CN102332156 A CN 102332156A CN 201110264780 A CN201110264780 A CN 201110264780A CN 201110264780 A CN201110264780 A CN 201110264780A CN 102332156 A CN102332156 A CN 102332156A
Authority
CN
China
Prior art keywords
image
filtering
domain filtering
fingerprint
frequency domain
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.)
Pending
Application number
CN201110264780A
Other languages
Chinese (zh)
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.)
CHANGZHOU LENCITY INFORMATION TECHNOLOGY Co Ltd
Jiangxi University of Finance and Economics
Original Assignee
CHANGZHOU LENCITY INFORMATION TECHNOLOGY Co Ltd
Jiangxi University of Finance and Economics
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 CHANGZHOU LENCITY INFORMATION TECHNOLOGY Co Ltd, Jiangxi University of Finance and Economics filed Critical CHANGZHOU LENCITY INFORMATION TECHNOLOGY Co Ltd
Priority to CN201110264780A priority Critical patent/CN102332156A/en
Publication of CN102332156A publication Critical patent/CN102332156A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Collating Specific Patterns (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to the technical field of fingerprint image enhancement treatment and particularly relates to a fingerprint enhancement method based on a time domain and a frequency domain simultaneously for filtering, and a final enhanced image can be obtained by performing frequency domain filtering enhancement treatment while performing time domain filtering enhancement treatment and further performing fusion treatment on the enhanced image. By adopting the fingerprint enhancement method, the enhanced filtering can be respectively performed on a fingerprint image at the time domain and the frequency domain, thereby overcoming the deficiencies of the traditional method and greatly enhancing the treatment effect of the fingerprint image.

Description

A kind of while is based on the fingerprint Enhancement Method of time-domain and frequency-domain filtering
Technical field
The present invention relates to fingerprint image enhancement process technical field, particularly a kind of while is based on the fingerprint Enhancement Method of time-domain and frequency-domain filtering.
Background technology
Traditional fingerprint Enhancement Method is mainly two kinds of time-domain filtering and frequency filterings.Time-domain filtering is mainly through the structure suitable filters; In time domain fingerprint image is carried out convolution algorithm; Thereby reach the purpose of enhancing, like methods such as wave filter, nonlinear diffusion wave filter, context wave filter and scale-space filtering based on Gabor filter.The frequency domain Fourier transform can convert in frequency and does product calculation do convolution algorithm in time domain; Thereby fast and effeciently image is handled; Reach the purpose that strengthens image, as based on direction frequency domain filter, anisotropic wave filter, directly or methods such as Short Time Fourier Transform analysis and " Log-Gabor " wave filter.But; These fingerprint Enhancement Method all are to use single wave filter; And the parameter of selective filter is mostly based on the streakline structure of fingerprint; Because fingerprint image has the characteristic of non-stationary, the streakline structure of inferior quality fingerprint image is complicated especially, and the effect of such filtering reinforcement method often can not be satisfactory.
Summary of the invention
Technical matters to be solved by this invention is: the enhancement method of fingerprint image of a kind of while based on time-domain filtering and frequency domain filtering is provided; In time domain the streakline of fingerprint image is repaired effectively; Select fingerprint image is strengthened filtering from direction and frequency domain respectively at frequency domain, thereby strengthen fingerprint image.
The technical solution adopted for the present invention to solve the technical problems is: a kind of while is based on the enhancement method of fingerprint image of time-domain filtering and frequency domain filtering; After in the time-domain filtering enhancement process, carrying out the frequency domain filtering enhancement process, again the enhancing image is carried out fusion treatment and finally strengthened image.
Described time-domain filtering enhancement process comprises former schematic partial plan picture normalization, asks local streakline direction, compensation filter and smoothing processing.
Described frequency domain filtering enhancement process comprises former figure global image normalization, asks the adjustment of directional diagram, Gabor Filtering Processing and gradation of image.
Described enhancing image co-registration is handled and is comprised that corresponding pixel points adds up and the gray scale adjustment.
The invention has the beneficial effects as follows: this fingerprint Enhancement Method can strengthen filtering to fingerprint image respectively in time domain and frequency domain, has overcome the deficiency of classic method, has greatly strengthened the treatment effect of fingerprint image.
Description of drawings
Below in conjunction with accompanying drawing and embodiment the present invention is further specified;
Fig. 1 is the theoretical block diagram of implementation of the present invention;
Fig. 2 is the theoretical block diagram of time-domain filtering of the present invention;
Fig. 3 is the former figure of fingerprint of time domain fingerprint image enhancing step 1 and the image after the local normalization;
Fig. 4 is the theoretical block diagram of frequency domain filtering of the present invention;
Fig. 5 is that time-domain and frequency-domain of the present invention strengthens the image co-registration block diagram;
Fig. 6 is that fingerprint image 1 strengthens experimental result;
Fig. 7 is that fingerprint image 2 strengthens experimental result.
Among Fig. 3, (a) be the former figure of fingerprint of time domain fingerprint image enhancing step 1; (b) be image after time domain fingerprint image enhancing step 1 local normalization; Among Fig. 6, (a) be the former figure of fingerprint that fingerprint image 1 strengthens experiment; (b) be that the fingerprint that fingerprint image 1 strengthens behind the time-domain filtering of experiment strengthens image; (c) be that the fingerprint that fingerprint image 1 strengthens behind the frequency domain filtering of experiment strengthens image; (d) the fingerprint enhancing image after the time-domain filtering that strengthens experiment for fingerprint image 1 strengthens with frequency domain filtering simultaneously; Among Fig. 7, (a) be the former figure of fingerprint that fingerprint image 2 strengthens experiment; (b) be that the fingerprint that fingerprint image 2 strengthens behind the time-domain filtering of experiment strengthens image; (c) be that the fingerprint that fingerprint image 2 strengthens behind the experiment frequency domain filterings strengthens image; (d) be the fingerprint image 2 fingerprint enhancing image after strengthening the experiment time-domain filterings and strengthening simultaneously with frequency domain filtering.
Embodiment
The time domain fingerprint image strengthens
Streakline fracture and streakline fuzzy region to fingerprint can be extracted the fact that exerts an influence to fingerprint characteristic; Construct the convolution template of suitable time domain filtering and choose reasonable wave filter; Combing and repairing are carried out in these streakline zones; Can reach the purpose that strengthens fingerprint image, the selection of filter parameter mainly is a structure convolution template parameter, and concrete steps as shown in Figure 2 are following:
Step 1: topography's normalization
Consider that each local intensity profile of image is unbalanced, thus adopt local normalization to handle to image, definition expectation average M among this paper 0=128, expectation variance V 0=128 * 128, local window size W * W is (8 * 8).Topography normalization be to each pixel img of image (i j) averages and variance for the part at center, again to central point img (i j) standardizes, thereby obtains better effect, and key step is following:
(1) obtains that (i j) is the average M and the variance V of the topography at center with img;
M ( i , j ) = 1 R * L Σ j - w / 2 j + w / 2 Σ i - w / 2 i + w / 2 img ( i , j ) , ( i , j ) ∈ img
V ( i , j ) = Σ j - w / 2 j + w / 2 Σ i - w / 2 i + w / 2 ( img ( i , j ) - M ( i , j ) ) 2 , ( i , j ) ∈ img - - - ( 1 )
(2) by following formula to central point img (i, j) processing of standardizing;
norimg ( i , j ) = M 0 + V 0 ( img ( i , j ) - M ( i , j ) ) V ( i , j ) , ( i , j ) ∈ img ∩ V ( i , j ) ≠ 0
norimg(i,j)=M0+100*(img(i,j)-M(i,j),(i,j)∈img?∩V(i,j)=0(2)
(3) (i, j) edge value is handled to norimg.
if?norimg(i,j)<0,norimg(i,j)=0,
esle?if?norimg(i,j)>255,norimg(i,j)=255(3)
Result after normalization is handled as stated above is shown in Fig. 3 (b).
Step 2: ask local streakline direction
Fingerprint image adopts gradient algorithm that the fingerprint local direction is estimated after normalization is handled.Fingerprint image after the normalization is divided into non-overlapping, the size of piece is W * W, and (i j) for center and mask operator carry out convolution, like the sobel operator, obtains the horizontal gradient component G of each pixel in the piece with the pixel of piece central point X(u, v), VG (vertical gradient) component G Y(u, v) and whole horizontal gradient component G XX(i, j), VG (vertical gradient) component G XY(i, j), through the upwards horizontal gradient component and (4)~(6) the formula acquisition by formula of VG (vertical gradient) component of point of counterparty.
Gxx = Σ u = i - w / 2 i + w / 2 Σ v = j - w / 2 j + w / 2 ( Gx 2 ( u , v ) - Gy 2 ( u , v ) ) - - - ( 4 )
Gxy = Σ u = i - w / 2 i + w / 2 Σ v = j - w / 2 j + w / 2 2 Gx ( u , v ) Gy ( u , v ) - - - ( 5 )
O ( x , y ) = 1 2 arctan ( Gxy Gxx ) - - - ( 6 )
In order accurately to estimate the streakline direction; After to the horizontal gradient component and VG (vertical gradient) component of obtaining pixel; It is carried out gaussian filtering handle, then directional diagram is handled, then directional diagram is carried out smoothly with Gaussian filter; Carry out gaussian filtering at last again and handle, for next step compensation filter is got ready.
Step 3: compensation filter and smoothing processing
Repair to cut off crestal line and suppress non-crestal line information according to the half-tone information of the direction of crestal line and neighborhood and compensate filtering, smothing filtering is that low frequency strengthens the spatial domain filtering technique, and its purpose is a Fuzzy Processing and abate the noise.The smothing filtering of spatial domain adopts simple average method, promptly asks the average gray value of this pixel neighborhood pixels point.The size of neighborhood is directly related with level and smooth effect, and neighborhood is big more, and smooth effect is good more, but neighborhood is excessive, marginal information is lost in a large number, thereby output image is thickened, and therefore needs choose reasonable neighborhood size.By formula (7~9) compensate directional diagram:
enhimg ( i , j ) = Σ m = - w / 2 w / 2 Σ n = - k / 2 k / 2 I ( i ′ , j ′ ) / [ ( w + 1 ) * ( h + 1 ) ] - - - ( 7 )
i′=i+mcos(oimg(i,j))+n*sin(oimg(i,j)) (8)
j′=j-msin(oimg(i,j))+n*cos(oimg(i,j)) (9)
Wherein oimg (i is the local direction obtained in second step j), and wherein window size is w * h, w=4, h=14, the selection of these values is to come according to the empirical value of repeatedly experiment.
The frequency domain fingerprint image strengthens
After time domain compensation filtering enhancing, the fingerprint image lines is comparatively clear, but in order to make fingerprint image more effectively extract characteristic, can take to strengthen based on the frequency domain of Gabor wave filter, and promptly the frequency domain fingerprint image strengthens.The frequency domain Fourier transform can convert in frequency and does product calculation do convolution algorithm in time domain, thereby fast and effeciently image is handled, and reaches the purpose that strengthens image.The frequency filter based on the Gabor wave filter of structure has good direction and frequency selectivity, can on direction and frequency domain selectivity, strengthen filtering to fingerprint image respectively at frequency domain, and its concrete performing step as shown in Figure 4 is following:
Step 1: global image normalization
After the process enhancement process first time; Each local intensity profile of fingerprint image is relatively more balanced; Can adopt the global image normalization that each gray values of pixel points on the image is adjusted; Adjust to one to the contrast of different original images and average gray fixedly on the rank, reduce the influence that different fingerprint image contrast differences are brought, definition expectation gray average mean 0=180, expectation gray variance var 0=3600, key step is following:
(1) gray-scale value of calculated fingerprint image and variance
M = Σ 1 R Σ 1 L enhimg ( i , j ) R * L , ( i , j ) ∈ enhimg - - - ( 10 )
V = Σ 1 R Σ 1 L ( enhimg ( i , j ) - M ) 2 R * L , ( i , j ) ∈ enhimg - - - ( 11 )
R, L be the row value and the train value of presentation video respectively, and enhimg is the image information after strengthening for the first time.
(2) by following formula image is standardized
G ( i , j ) = mean 0 + var 0 * ( enhimg ( i , j ) - M ) 2 / V , ( i , j ) ∈ enhimg > M mean 0 - var 0 * ( enhimg ( i , j ) - M ) 2 / V , ( i , j ) ∈ enhimg , otherwise - - - ( 12 )
Step 2: ask directional diagram
In the image on the different directions point corresponding gray scale be vicissitudinous, the gray scale in crestal line differs very little, and is maximum but the gray scale on the vertical direction differs, and can obtain directional diagram through compute gradient, key step is following:
(1) utilizes sobel to calculate each pixel and do two convolution, obtain the VG (vertical gradient) component G of each pixel Y(i is j) with horizontal gradient component G X(i, j);
Figure BDA0000089712700000064
(2) the direction expansion with each pixel gradient vector is twice;
Gsx(i,j)=Gx(i,j) 2-Gy(i,j) 2
Gsy(i,j)=2*Gx(i,j)*Gy(i,j) (13)
(3) ask the average gradient at each pixel place vectorial, local field size is W 1* W 1, (W 1=15), formula is:
A _ Gsx ( i , j ) = Σ u = - W 1 / 2 + i u = W 1 / 2 + i Σ v = - W 1 / 2 + j v = - W 1 / 2 + j Gsx ( u , v )
A _ Gsy ( i , j ) = Σ u = - W 1 / 2 + i u = W 1 / 2 + i Σ v = - W 1 / 2 + j v = - W 1 / 2 + j Gsyy ( u , v ) - - - ( 14 )
Then
Q ( i , j ) = 0.5 * arctan ( A _ Gsy ( i , j ) / A _ Gsx ( i , j ) ) , A _ Gsx ( i , j ) ≠ 0 0 , otherwise - - - ( 15 )
Figure BDA0000089712700000074
that obtains can not express all directions of fingerprint, and the angle of direction by formula (16) is adjusted:
O ( i , j ) = Q ( i , j ) + &pi; 2 , A _ Gsy > 0 Q ( i , j ) + &pi; , A _ Gsx ( i , j ) < 0 &cap; A _ Gsy ( i , j ) > 0 Q ( i , j ) , otherwise - - - ( 16 )
Wherein (i j) is adjusted angle to O.
Step 3:Gabor Filtering Processing
The Gabor wave filter has good direction and frequency selectivity, can when preserving correct crestal line and valley line structure, remove noise with it as frequency filter, and this paper adopts the Gabor wave filter to realize the enhancing of fingerprint image.Definition W 2=10 are filtering mask size, O X=3 are the logical size of frequency band, O YBe the logical size of direction band, f=0.1 is the streakline average frequency, and key step is following:
(1) utilizes the gabor wave filter to be used in reference to print image and handle, strengthen template antithetical phrase piece with formula (17) as sub-piece and strengthen;
E ( i , j ) = exp ( - 1 2 ( x 2 Ox 2 + y 2 Oy 2 ) ) * cos ( 2 * &pi; * f * x )
x=v*sin(O(i,j))+u*cos(O(i,j))
y=v*cos(O(i,j))-u*sin(O(i,j)) (17)
(2) the sub-piece after will strengthening by formula (18) be merged into the complete finger print image;
E ( i , j ) = &Sigma; u = - W 2 / 2 W 2 / 2 &Sigma; v = - W 2 / 2 W 2 E ( u , v ) + G ( i - u , j - v ) - - - ( 18 )
Wherein E is the image after strengthening, and the Gabor wave filter can reasonably carry out filtering to fingerprint image at each topography's piece, improves signal noise ratio (snr) of image, thereby can preserve comparatively complete sum fingerprint ridge information clearly.
Strengthen image co-registration
In conjunction with like Fig. 5, strengthen image co-registration and may further comprise the steps:
Step 1: the fusion of pixel;
Time domain enhancing fingerprint image and frequency domain enhancing fingerprint image are added up according to each corresponding pixel,
Realize the fusion of pixel.
Step 2: gradation of image adjustment;
The enhancing image that merges is at last carried out the gray scale adjustment, and key step is following:
(1) obtains enhancing back fingerprint image maximum gradation value ma and minimum gradation value mi, get difference val=ma-mi then;
(2) each pixel with E deducts minimum value;
E(i,j)=E(i,j)-mi (19)
(3) utilize formula (20) that each gray values of pixel points is adjusted.
E2(i,j)=E(i,j)*256/val (20)
Specific embodiment
Present embodiment adopts the Matlab programming language on the platform of matlab7.0, to realize; Chosen two fingerprint images and done experiment test, experimental result such as Fig. 6 and shown in Figure 7 can find out from experimental result; For different images; The effect that different filtering modes are obtained is different, Fig. 6 and Fig. 7 former figure before the effect after out-of-date, frequency domain strengthen is superior to strengthening, but two pieces of fingerprint images are best passing through the reinforced effects that obtains after time-domain filtering and the frequency domain filtering enhancement process.Reason is: time-domain filtering is handled the contrast of the crestal line and the valley line that can effectively improve fingerprint, and along the combing of fingerprint local direction, the crestal line of compensation fingerprint image; Particularly streakline ruptures and the streakline fuzzy region; Frequency domain filtering utilizes frequency bandpass filter on frequency domain, to carry out the fingerprint image enhancing, has good direction and selects and frequency selective characteristic, can make full use of the frequency and the directional information of streakline in the regional area; Outstanding streakline inherent structure; Remove picture noise, can fingerprint ridge line and the valley line structure is distortionless remain, reach the purpose of enhancing fingerprint image.Fig. 6 (d) and Fig. 7 (d) are the images of handling through time-domain and frequency-domain, combine the advantage of time-domain filtering and frequency domain filtering, and therefore, the reinforced effects of two images is best.

Claims (4)

1. while is characterized in that based on the enhancement method of fingerprint image of time-domain filtering and frequency domain filtering: in the time-domain filtering enhancement process, carry out the frequency domain filtering enhancement process, carry out fusion treatment and finally strengthened image strengthening image again.
2. a kind of while according to claim 1 is characterized in that based on the enhancement method of fingerprint image of time-domain filtering and frequency domain filtering: described time-domain filtering enhancement process comprises former schematic partial plan picture normalization, asks local streakline direction, compensation filter and smoothing processing.
3. a kind of while according to claim 1 is characterized in that based on the enhancement method of fingerprint image of time-domain filtering and frequency domain filtering: described frequency domain filtering enhancement process comprises former figure global image normalization, asks the adjustment of directional diagram, Gabor Filtering Processing and gradation of image.
4. a kind of while according to claim 1 is characterized in that based on the enhancement method of fingerprint image of time-domain filtering and frequency domain filtering: described enhancing image co-registration is handled and is comprised that corresponding pixel points adds up and the gray scale adjustment.
CN201110264780A 2011-09-07 2011-09-07 Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering Pending CN102332156A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110264780A CN102332156A (en) 2011-09-07 2011-09-07 Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110264780A CN102332156A (en) 2011-09-07 2011-09-07 Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering

Publications (1)

Publication Number Publication Date
CN102332156A true CN102332156A (en) 2012-01-25

Family

ID=45483921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110264780A Pending CN102332156A (en) 2011-09-07 2011-09-07 Fingerprint enhancement method based on time domain and frequency domain simultaneously for filtering

Country Status (1)

Country Link
CN (1) CN102332156A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109214160A (en) * 2018-09-14 2019-01-15 温州科技职业学院 A kind of computer network authentication system and method, computer program
CN109858418A (en) * 2019-01-23 2019-06-07 上海思立微电子科技有限公司 The treating method and apparatus of fingerprint image

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101329726A (en) * 2008-07-30 2008-12-24 电子科技大学 Method for reinforcing fingerprint image based on one-dimensional filtering

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101329726A (en) * 2008-07-30 2008-12-24 电子科技大学 Method for reinforcing fingerprint image based on one-dimensional filtering

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JUCHENG YANG ET AL: "Effective Algorithm for Fingerprint Enhancement in Poor Images", 《IPCV 2008》 *
廖开阳 等: "基于Gabor滤波器的指纹图像快速增强", 《计算机工程与应用》 *
杨巨成 等: "指纹图像增强新方法", 《第十五届全国图象图形学学术会议》 *
王莹 等: "指纹图像增强算法研究", 《科学技术与工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109214160A (en) * 2018-09-14 2019-01-15 温州科技职业学院 A kind of computer network authentication system and method, computer program
CN109858418A (en) * 2019-01-23 2019-06-07 上海思立微电子科技有限公司 The treating method and apparatus of fingerprint image
CN109858418B (en) * 2019-01-23 2021-10-15 上海思立微电子科技有限公司 Fingerprint image processing method and device

Similar Documents

Publication Publication Date Title
CN103093433B (en) Natural image denoising method based on regionalism and dictionary learning
CN103295191B (en) Multiple scale vision method for adaptive image enhancement and evaluation method
CN105654433A (en) Color image enhancement method based on improved multi-scale Retinex
CN104574337B (en) Based on the image enchancing method that bilateral gamma correction and multi-scale image merge
CN102289792A (en) Method and system for enhancing low-illumination video image
CN102096909B (en) Improved unsharp masking image reinforcing method based on logarithm image processing model
CN102890820B (en) Based on shearlet conversion and the image de-noising method of Wiener filtering
CN103778613B (en) A kind of adaptive polarization SAR images filter method of window
CN101615292B (en) Accurate positioning method for human eye on the basis of gray gradation information
CN101930604B (en) Infusion method of full-color image and multi-spectral image based on low-frequency correlation analysis
CN102298773B (en) Shape-adaptive non-local mean denoising method
CN104200437A (en) Image defogging method
CN106556940A (en) A kind of background suppression method in TFT LCD screens automatic optics inspection
CN104182947A (en) Low-illumination image enhancement method and system
Pandey et al. Enhancing the quality of satellite images by preprocessing and contrast enhancement
CN104580829A (en) Terahertz image enhancing method and system
CN103020918A (en) Shape-adaptive neighborhood mean value based non-local mean value denoising method
CN101877124A (en) Method and system for filtering medical image
CN107909112A (en) The fusion method that a kind of infrared light intensity is combined with polarization image multiclass argument
CN105913400A (en) Device for obtaining high-quality and real-time beautiful image
CN107784639A (en) A kind of polygon filtering and noise reduction method of unmanned aerial vehicle remote sensing image improvement
CN104616259B (en) A kind of adaptive non-local mean image de-noising method of noise intensity
CN107169946A (en) Image interfusion method based on non-negative sparse matrix Yu hypersphere color transformation
CN107203769A (en) Image characteristic extracting method based on DCT and LBP Fusion Features
CN101957984B (en) Image de-noising method based on parametric estimation of non-local shrinkage factor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20120125