CN101615248B - Age estimation method, equipment and face recognition system - Google Patents

Age estimation method, equipment and face recognition system Download PDF

Info

Publication number
CN101615248B
CN101615248B CN2009101310593A CN200910131059A CN101615248B CN 101615248 B CN101615248 B CN 101615248B CN 2009101310593 A CN2009101310593 A CN 2009101310593A CN 200910131059 A CN200910131059 A CN 200910131059A CN 101615248 B CN101615248 B CN 101615248B
Authority
CN
China
Prior art keywords
data
age
dimensionality reduction
subclass
estimation
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.)
Expired - Fee Related
Application number
CN2009101310593A
Other languages
Chinese (zh)
Other versions
CN101615248A (en
Inventor
王蕴红
左坤隆
郝韬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huawei Technologies Co Ltd
Beihang University
Beijing University of Aeronautics and Astronautics
Original Assignee
Huawei Technologies Co Ltd
Beihang University
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 Huawei Technologies Co Ltd, Beihang University filed Critical Huawei Technologies Co Ltd
Priority to CN2009101310593A priority Critical patent/CN101615248B/en
Publication of CN101615248A publication Critical patent/CN101615248A/en
Application granted granted Critical
Publication of CN101615248B publication Critical patent/CN101615248B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

The embodiment of the invention provides an age estimation method and equipment. The method comprises the following steps: processing face images of users to extract initial characteristic data; selecting subsets of the initial characteristic data as low-dimension characteristic data with dimensionality less than that of the initial characteristic data; carrying out relevant component analysis on the low-dimension characteristic data to obtain training data; according to the obtained training data, training regression analysis parameters used for age estimation; adopting the regression analysis parameters to estimate the ages of the users. The embodiment of the invention can reduce the computation of the age estimation in the face recognition system.

Description

Age estimation method, equipment and face identification system
Technical field
The present invention relates to face identification system, more specifically, relate to age estimation method and equipment in the face identification system.
Background technology
Realize the face identification system of a robust, inevitably will solve the influence of factors such as for example illumination, expression, attitude and age.Estimation of Age is exactly the method according to pattern classification, utilizes the characteristic of facial image, and this individual age is estimated.Man-machine interaction based on the age has huge potential application in daily life.Yet the age recognition technology is but related to seldom under study for action automatically.A chief reason is that the characteristic of some uniquenesses of the aging effect of people face makes estimation of Age need use some off-gauge sorting techniques.This also makes estimation of Age become a challenging problem.
If can solve the age identification problem, then be expected to realize following target:
The man-machine interaction of given age: if computing machine can be judged user's age exactly, computing environment and type of interaction can be made adjustment according to user's age.Be different from the man-machine interaction of standard, such system can combine to guarantee that the teenager is away from unhealthy information with the control of Internet access security.
Facial image index based on the age: age identification automatically can be used for the facial image in the database is sorted; Also can it be used in the electron album, the user can come their photo is sorted, manages through setting the range of age easily like this.
Realization from moving face ageing system: age identification automatically depends on because the perception and the classification of the variation of the facial characteristics that change of age produced, and same method also is needed from the aging simulation system of moving face.
A kind of age recognition methods based on EHMM (built-in type hidden Markov model) is arranged in the prior art; At first; Through analysis, confirm the nonlinear relationship between the variation of age and people's face key feature, and set up nonlinear model based on this relation to great amount of samples; Utilize PARI (partial ageingratio image, local ageing ratio chart picture) to estimate the texture of aging people's face then; At last, utilize the characteristic of from the image of rebuilding, extracting to train EHMM to be used for recognition of face.
In realizing process of the present invention, the inventor finds to exist in the above-mentioned prior art following problem:
The realization cost of existing method is bigger, extract a large amount of profile informations and texture information respectively to the sample set of each age bracket, and operand is very big.
Summary of the invention
The embodiment of the invention provides age estimation method, equipment and the face identification system in a kind of face identification system, to reduce the operand of estimation of Age.
The embodiment of the invention provides a kind of age estimation method, comprising: processing face images of users is to extract initial characteristic data; The subclass of selecting said initial characteristic data is as low dimensional feature data, and the dimension of wherein said low dimensional feature data is less than the dimension of said initial characteristic data; To the constituent analysis of being correlated with of said low dimensional feature data, to obtain training data; According to the said training data that is obtained, training is used for the regretional analysis parameter of estimation of Age; And utilize said regretional analysis parameter to estimate said user's age.
The embodiment of the invention also provides a kind of estimation of Age equipment, comprising: the primitive character extraction module, and processing face images of users is to extract initial characteristic data; Feature selection module, the subclass of selecting said initial characteristic data is as low dimensional feature data, and the dimension of wherein said low dimensional feature data is less than the dimension of said initial characteristic data; Relevant component analysis module is to the constituent analysis of be correlated with of said low dimensional feature data, with the acquisition training data; Training module according to by the said training data that said relevant component analysis module obtained, is trained the regretional analysis parameter that is used for estimation of Age; And estimation module, utilize said regretional analysis parameter to estimate said user's age.
The embodiment of the invention also provides a kind of face identification system, comprises above-mentioned estimation of Age equipment.
The embodiment of the invention adopts the method for relevant constituent analysis to carry out dimensionality reduction for facial image, adopts the mode that returns to carry out estimation of Age, can reduce the operand in the age estimation method.
Description of drawings
In order to be illustrated more clearly in the technical scheme of the embodiment of the invention; To do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art below; Obviously, the accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills; Under the prerequisite of not paying creative work property, can also obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the process flow diagram that illustrates according to the age estimation method of the embodiment of the invention.
Fig. 2 is the process flow diagram that illustrates according to the relevant constituent analysis step of the embodiment of the invention.
Fig. 3 is the comparison synoptic diagram according to the average error result of the quadratic regression estimation of Age of the embodiment of the invention and prior art.
Fig. 4 is the test sample book exemplary plot.
Fig. 5 is the comparison synoptic diagram according to calculated performance on the FG-Net database of the embodiment of the invention and prior art.
Fig. 6 is the synoptic diagram that illustrates according to the estimation of Age equipment of the embodiment of the invention.
Fig. 7 is the synoptic diagram that illustrates according to the face identification system of the embodiment of the invention.
Embodiment
To combine the accompanying drawing in the embodiment of the invention below, the technical scheme in the embodiment of the invention is carried out clear, intactly description, obviously, described embodiment is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills are not making the every other embodiment that is obtained under the creative work prerequisite, all belong to the scope of the present invention's protection.
Fig. 1 is the process flow diagram that illustrates according to the age estimation method 100 of the embodiment of the invention.
In this method 100, at S105, processing face images of users is to extract initial characteristic data.By being processed data set (for example, facial image) through calculating or measure generation one stack features value.Because this is one group of data that data set itself gets, therefore be called as initial characteristic data.
In addition, in S105, can carry out pre-service to user's facial image earlier.For example, with the image size normalization.
At S110, the subclass of selecting initial characteristic data is as low dimensional feature data.The purpose of feature selecting is exactly in order to reduce the dimension of feature space, keep as far as possible the effectively information of classification that can be used for.For facial image; Each width of cloth facial image data all constitutes a point of higher dimensional space, but have in numerous dimensions of higher dimensional space can not effectively embody characteristic of division, so carry out linear transformation; The deletion redundant information obtains the most effectively characteristic of division information.
Because primitive character quantity maybe be a lot, in order to adapt to engineering time, must reduce number of features as far as possible.Feature selecting is exactly that the numerical value that from primitive character, obtains is selected the subclass that is used to classify according to certain criterion, as the characteristic of division of dimensionality reduction (dropping to the m dimension from the n dimension).In other words, feature selecting is exactly to utilize the feature space of low dimension of existing characteristic parameter structure, and the useful information that contains in the primitive character is mapped on a few key character, ignores less important information.Just to a n dimension original feature vector
X=[x 1,x 2,…,x n,] T
Through certain conversion, carry out dimensionality reduction, produce the proper vector of low dimension
Y=[y 1,y 2,…,y m,] T,m<n
Wherein Y contains the principal character of X.
For high dimensional data is mapped with many kinds of methods in low n-dimensional subspace n, principal component analysis (PCA) (PCA) is simple and effective a kind of.PCA is proposed by A.Turk and A.P.Pentland, and it is an expansion of Karhunen-Loeve transformation.Karhunen-Loeve transformation is the method for a kind of typical feature extraction and data representation, is the optimal mapping that under the meaning of square error minimum, obtains data compression, is applicable to probability distribution arbitrarily.Principal component analysis (PCA) is the best approach of carrying out feature extraction at present, has extensively been applied in the middle of a lot of fields such as pattern-recognition and computer vision.PCA is to locate a low n-dimensional subspace n that can satisfy the minimum variance criterion and expresses original high dimensional data.For given data set (initial characteristic data) X = { x i } i = 1 n , If the definition projection matrix is p, so just find to make following formula value reach maximum subspace:
p=arg?max?p TSp
Here, scatter matrix
S = Σ i = 1 n ( x - x ‾ ) ( x i - x ‾ ) T ,
Wherein x is { x i} I=1 nAverage vector.
At S115, to the low dimensional feature data that in S110, the obtain constituent analysis of being correlated with, to obtain training data.
Relevant constituent analysis (RCA; Relevant Component Analysis) can find and reduce the interior overall irrelevant variation of data set.The RCA method distributes bigger weights for relevant dimension through linear transformation, for uncorrelated dimension is distributed less weights, thereby changes the feature space of data representation.The RCA method is the class methods between principal component analysis (PCA) (PCA) and linear discriminant analysis (LDA); It has utilized the sample subclass (sample in each sample subclass belongs to classification identical but the possibility Unknown Label) with relation of equivalence; Estimate relevant dimension and irrelevant dimension through sample subclass, and change (White Transformation) through leucismus and reduce non-correlation with relation of equivalence.
The RCA algorithm is equivalent to utilize the equivalence information of part sample to do differentiation dimensionality reduction and decorrelation.
Fig. 2 is the process flow diagram that illustrates according to the relevant constituent analysis step S115 of the embodiment of the invention.
At S1151, at least one equivalent samples subclass of the low dimensional feature data of structure.Each equivalent samples subclass belongs to age-grade label.For a whole sample, can obtain the equivalence information of its part sample, promptly in these part samples, which sample belongs to identical classification.Need not know the label of classification; Can be expressed as the age label that need not know these facial images here, as long as but can confirm that some sample belongs to same age A1, some sample belongs to same age A2... (wherein A1, A2 the unknown).Can utilize this equivalence information, the original image collection is carried out dimensionality reduction and eliminates non-correlation.
At S1154, the equivalent samples subclass that is utilized among the S1151 to be constructed is differentiated dimensionality reduction to low dimensional feature data, to obtain by the further dimensionality reduction sample of dimensionality reduction.
In the process of dimensionality reduction; Covariance matrix as whole sample set (promptly between the class of the part sample with equivalence information capable of using; Estimation of covariance matrix value between class low dimensional feature data); Utilize cFDA (combined Fisher discriminant analysis, combination Fisher discriminatory analysis) to carry out dimensionality reduction then.
Particularly, for a given data set (low dimensional feature data) { X i} I=1 NAnd n the equivalent samples subclass that obtains from this data centralization C j = { x Ij } i = 1 n j C j = { x Ij } i = 1 n j , Calculate the interior covariance matrix of class of subclass of equal value:
C ^ = 1 N Σ j = 1 n Σ i = 1 n j ( x ji - m j ) ( x ji - m j ) T
Wherein N is the number of samples of data centralization, m jIt is subclass C j = { x Ij } i = 1 n j In the mean value of sample.
Use
Figure G2009101310593D00055
data set is differentiated dimensionality reduction, dimension reduction method is described below.
The dimension of representing low dimensional feature data with D.For one group of given equivalent samples subclass { C j} J=1 n, the order of covariance matrix in the class of calculating equivalent samples subclass R = Σ j = 1 n ( | C j | - 1 ) , Wherein | C j| represent the size of j equivalent samples subclass.
If D>R, then use PCA with the data dimension dimensionality reduction to α R, 0<α<1 wherein.
The covariance matrix S of computing whole t, and with the estimated value of covariance matrix in the class of equivalent samples subclass as covariance in the class of whole sample set, i.e. S w=C finds the solution following equality
A ~ = arg A max AS t A t AS w A t
And obtain the expression (i.e. " dimensionality reduction sample data " X ') of data in new lower dimensional space with
Figure G2009101310593D00058
.
At S1157, reduce the non-correlation of dimensionality reduction sample data, to obtain training data.
In the process that reduces non-correlation, definable is for the sample subclass of same age label, and the direction that its variance is big is irrelevant direction.Can adopt leucismus to bring covariance matrix is carried out conversion, make that the weights of irrelevant direction are lower, so just reduce the age non-correlation in the sample subclass.
Particularly, can recomputate covariance matrix in the class of subclass of equal value behind the dimensionality reduction
C ~ = 1 N Σ j = 1 n Σ i = 1 n j ( x ~ ji - m j ) ( x ~ ji - m j ) T
Where
Figure G2009101310593D00062
indicates an equivalent subset
Figure G2009101310593D00063
In the low-dimensional subspace representation.
The leucismus that calculates
Figure G2009101310593D00064
changes:
C ~ : W = C ~ - 1 / 2
Wherein W is that leucismus changes the result.
The sample that data are concentrated carries out following conversion
X new=WX′
Wherein X ' expression is the data point (dimensionality reduction sample data) behind the dimensionality reduction, X NewIt is the training data that is used for the regretional analysis parameter of sport career age estimation.
Return Fig. 1, after S115 obtained training data, at S120, according to the training data that in S115, obtains, training was used for the regretional analysis parameter of estimation of Age.
Regretional analysis (regression analysis) is a kind of statistical analysis technique of confirming complementary quantitative relationship between two or more parameter.
Represent that afterwards, estimation of Age has just become a multiple linear regression problem in the PCA space to sample set (low dimensional feature data) finding a low dimension PCA subspace.Following formula has provided concrete regression model:
age = f ( M ) : ⇔ L ^ = f ( y ) ^
Here
Figure G2009101310593D00067
is the age label that estimation obtains; F () is unknown regression equation, is an estimation of regression equation.For matrix form, following formula can be expressed as:
L = Y ^ B + e ,
Var(e)=σ 2I
The L here is an age label vector;
Figure G2009101310593D00069
is a known matrix, and this matrix is made up of the PCA proper vector of sample set.Var () representes variance (Variance), and I is a unit matrix.B needs our the coefficient vector parameter in training stage study.Error vector e is a unobservable random variable, and it is that 0 variance is σ that this stochastic variable should satisfy average 2And irrelevance.
B is tried to achieve by following formula:
B ^ = ( Y ^ T Y ) - 1 Y ^ T L = HL
Just can obtain regression equation then:
L ^ = Y ^ B ^
Residual vector e ^ = L - L ^ Satisfy:
E ( e ^ ) = 0 , Var ( e ^ ) = σ 2 ( I - H )
Take all factors into consideration result's accuracy and algorithm complex, can select quadratic equation to do regretional analysis.With the projection vector y of each sample image in the PCA subspace iImport as characteristic.
l ^ i = b ^ 0 + b ^ 1 T y i + b ^ 2 T y i 2
Here,
Figure G2009101310593D00072
Be image pattern x iThe estimation age,
Figure G2009101310593D00073
Be side-play amount,
Figure G2009101310593D00074
Be to estimate the regression equation coefficient vector.y iIt is the sample data behind the dimensionality reduction.
Therefore, can obtain:
L ^ = [ l ^ , . . . , l ^ n ] T
B ^ T = [ b ^ 0 b ^ 1 ( 1 ) . . . b ^ 1 ( d ) b ^ 2 ( 1 ) . . . b ^ 2 ( d ) ] T
Y ^ = [ 1 n × 1 [ y 1 . . . y n ] T [ y 1 2 . . . y 2 2 ] T ]
Wherein 1 N * 1Each element of representing a n * 1 dimension all is 1 matrix.
For sample set, above-mentioned formula becomes corresponding matrix form:
age = offset + w 1 T b + w 2 T ( b 2 )
Wherein age is the estimation age vector of corresponding sample set, and offset is the side-play amount vector, and b is a sample set, w 1, w 2It is the regression coefficient matrix.
Can obtain above-mentioned matrix of coefficients w through the training sample training 1, w 2And offset vector offset.We have just obtained the quadratic function value like this.
In S120, obtain after the regretional analysis parameter,, utilize the regretional analysis parameter that is obtained to come the age of estimating user at S125.After obtaining the quadratic regression analytical parameters, test sample book as function input substitution formula, just can be obtained the estimation of Age value of each sample.
Utilized the part classification information of sample set based on the feature dimension reduction method of RCA, i.e. equivalent samples subclass, utilization should equivalence sample subclass be found irrelevant direction, and utilized leucismus to bring the non-correlation that reduces sample set.The effect of RCA dimensionality reduction is the application that combines quadratic regression, reduces non-correlation and dimensionality reduction, not only can save the expense of computing, can also improve to return the precision of estimating.This is because the parameter that returns too much can need more training sample, and in reality, is difficult to satisfy the requirement of large sample amount, and RCA has solved this problem well.Each step of embodiment of the invention method can be adjusted order according to actual needs.
Fig. 3 is the comparison synoptic diagram according to the average error result of the quadratic regression estimation of Age of the embodiment of the invention and prior art.Fig. 4 is the test sample book exemplary plot.
Inventor of the present invention tests on the MORPH storehouse, and the data that will pass through after the PCA dimension-reduction treatment are divided into training sample and test sample book two parts according to 1: 1 ratio.Design for test is following:
Prior art adopts PCA+ to return.That is to say, use PCA that image set is carried out dimensionality reduction, carry out regression forecasting again, calculate the mean absolute error of tieing up the estimation of Age of 80 dimensions from 1 respectively.
According to one embodiment of present invention, adopt RCA+ to return: at first to use PCA with image dimensionality reduction to 80 dimension, adopt the cFDA dimensionality reduction again among the RCA then, carry out regression forecasting again, investigate the mean absolute error of the estimation of Age of from 1 to 80 dimension.(sample proportion that is used for subclass of equal value here is 30%, and when the cFDA in using RCA carried out dimensionality reduction, each age label was regarded as one type).
Particularly, the data that will pass through after the PCA dimension-reduction treatment are divided into training sample and test sample book two parts according to 1: 1 ratio.In test, at first original image is tieed up with PCA dimensionality reduction to 80.Wherein utilize PCA that original image is carried out dimension-reduction treatment.
Secondly, utilize RCA to reduce the non-correlation and the further dimensionality reduction of data.Here selected part has equivalent samples subclass of some images compositions of same age label, and constructs some equivalent samples subclass with all ages and classes label.Here we choose original sample and concentrate 30% sample composition equivalent samples subclass.Calculate covariance matrix in the class of this equivalence sample subclass (for fear of the linear discriminant dimensionality reduction of influence use to(for) regression result, with each age label as a type) then.
C ^ = 1 N Σ j = 1 n Σ i = 1 n j ( x ji - m j ) ( x ji - m j ) T
And using the interior covariance of class that obtains from the sample subclass, we carry out linear discriminant dimensionality reduction (as above S1154 is said) to original image set.For the sample behind the dimensionality reduction, we use leucismus to change further minimizing non-correlation
C ~ : W = C ~ - 1 / 2
At last, we adopt quadratic equation to carry out regression treatment to training sample.
l ^ i = b ^ 0 + b ^ 1 T y i + b ^ 2 T y i 2
Here,
Figure G2009101310593D00084
Be image pattern x iThe estimation age,
Figure G2009101310593D00085
Be side-play amount,
Figure G2009101310593D00086
Be to estimate the regression equation coefficient vector.y iIt is the sample data behind the dimensionality reduction.
For sample set, above-mentioned equality becomes corresponding matrix form:
age = offset + w 1 T b + w 2 T ( b 2 )
Age is the estimation age vector of corresponding sample set, and offset is the side-play amount vector, and b is a sample set, w 1, w 2It is the regression coefficient matrix.
Can obtain above-mentioned matrix of coefficients w through the training sample training 1, w 2And offset vector offset.So just obtained the quadratic function value.Test sample book as function input substitution formula, just can be obtained the estimation of Age value of each sample.
The result is as shown in Figure 3, and the X axle is represented the intrinsic dimensionality got, and the Y axle is represented mean absolute error (year).
According to Fig. 3 result, can see under the situation of the prior art of only using PCA that when the dimension of hanging down n-dimensional subspace n was lower, mean absolute error had tangible increase.And can keep lower stable mean absolute error according to the RCA method of the embodiment of the invention.And the RCA method has obtained the result who is superior to PCA when hanging down dimension.
In addition, for the method for the comparison embodiment of the invention and the difference of method on calculated performance and estimated efficiency of prior art, inventor of the present invention has carried out following test.
On FG-Net, test, compare with the method for PCA+SVR (Support Vector Regression, support vector regression), the method that in the test each algorithm was all adopted 1: 1 is divided into test set and training set to the FG-Net storehouse.Test findings and shown in Figure 5.Fig. 5 is the comparison synoptic diagram according to calculated performance on the FG-Net database of the embodiment of the invention and prior art, and wherein the X axle is the intrinsic dimensionality behind the dimensionality reduction, and the Y axle is second for the time of image training and test in the test, unit.By contrast; The method of the RCA+QR (Quadratic Regression, quadratic regression) that the embodiment of the invention adopts has apparent in view improvement on calculated performance, and with it more approaching result is arranged on estimated performance; Consider the requirement of practical application for the degree of accuracy of estimation of Age; Inventor of the present invention finds that the method that adopts RCA+QR significantly improves the performance of calculating on the basis that can keep the low mean absolute error of estimating.
Table 1 is the operation result on the FG-Net storehouse
Method SVM SVR LAR4 LAR8 LAR16 LAR32 RCA+QR
Mean absolute error 7.16 5.16 5.07 5.07 5.12 6.03 5.19
LAR4 wherein, LAR8, LAR16 and LAR32 are respectively the robustness methods of the part adjustment on the SVR basis, carried out, and the time complexity of its calculating will be higher than SVR.Only compared SVR and the RCA+QR complexity on computing time here.
Notice that the actual age of three samples shown in Figure 4 is respectively 16 years old, 39 years old and 62 years old.The estimation age of utilizing the method for the embodiment of the invention to obtain is 18.2,32.0 and 56.9.Accuracy of estimation for old-age group is lower, and its reason is mostly the training storehouse of being adopted is that living photo of the remote past is second-rate, and texture is fuzzy, and the non-interference of expression and attitude etc. is naturally arranged.
Employing can reduce the operand of estimation of Age according to the age estimation method of the embodiment of the invention.
Fig. 6 is the synoptic diagram that illustrates according to the estimation of Age equipment 200 of the embodiment of the invention.
This estimation of Age equipment 200 comprises primitive character extraction module 205, feature selection module 210, relevant component analysis module 215, training module 220 and estimation module 225.
Estimation of Age equipment 200 can be carried out above-described method 100.Particularly, primitive character extraction module 205 processing face images of users are to extract initial characteristic data.Feature selection module 210 selects the subclass of initial characteristic data as low dimensional feature data, and wherein the dimension of low dimensional feature data is less than the dimension of initial characteristic data.The constituent analysis of be correlated with of relevant 215 pairs of low dimensional feature data of component analysis module is with the acquisition training data.Training module 220, according to the training data that is obtained by relevant component analysis module 215, training is used for the regretional analysis parameter of estimation of Age.Estimation module 225 utilizes the regretional analysis parameter that is obtained by training module 220 training to come the age of estimating user.
Each module of estimation of Age equipment 200 with and the operation can repeat no more at this corresponding to above-mentioned age estimation method 100.
Employing can reduce the operand of estimation of Age according to the estimation of Age equipment of the embodiment of the invention.
Fig. 7 is the synoptic diagram that illustrates according to the face identification system 300 of the embodiment of the invention.This people's face device systems 300 can comprise above-mentioned estimation of Age equipment 200.It will be understood by those skilled in the art that according to actual needs people's face device systems 300 can also comprise other equipment (not shown), as be used to read facial image equipment, be used to equipment of storing facial image etc.
Employing can reduce the operand of estimation of Age according to the face identification system of the embodiment of the invention.
It will be appreciated by those skilled in the art that the module in the device among the embodiment can be distributed in the device of embodiment according to the embodiment description, also can carry out respective change and be arranged in the one or more devices that are different from present embodiment.The module of the foregoing description can be merged into a module, also can further split into a plurality of submodules.
The embodiment of the invention has proposed a kind of age estimation method based on relevant constituent analysis and quadratic regression, estimates for the sample that facial image is concentrated.The embodiment of the invention is used relevant constituent analysis when carrying out dimensionality reduction for facial image, reduced the non-correlation of facial image collection and in dimensionality reduction, considered the age discriminant information.For age-grade facial image, bring to changing the bigger lower weights of direction (irrelevant dimension) distribution through leucismus, thereby having the facial image of same age, assurance keeps less dispersion at lower dimensional space.On the other hand, the embodiment of the invention is applied to regretional analysis in the estimation of Age.Particularly, at first, utilize principal component analysis (PCA) that sample is carried out dimensionality reduction.Secondly, utilize RCA to reduce the non-correlation of sample set, and further select suitable differentiation dimensionality reduction space.In the dimensionality reduction subspace, use quadratic regression, set up recurrence and anticipation function.
The embodiment of the invention can be applicable to people's age is estimated, for making marks and put image library in order in the video image storage technical support is provided, retrieves to wait for suspect and also can dwindle range of search, the raising recall precision.In addition, the embodiment of the invention can provide technical foundation for recognition of face, character classification by age etc.
It should be noted that the embodiment of the invention described employing PCA and carried out dimensionality reduction, but it will be understood by those skilled in the art that other dimension reduction methods also can be applicable to the present invention, thereby also within the scope of the invention.
Those of ordinary skills can recognize; The unit and the algorithm steps of each example of describing in conjunction with embodiment disclosed herein; Can realize with electronic hardware, computer software or the combination of the two; For the interchangeability of hardware and software clearly is described, the composition and the step of each example described prevailingly according to function in above-mentioned explanation.These functions still are that software mode is carried out with hardware actually, depend on the application-specific and the design constraint of technical scheme.The professional and technical personnel can use distinct methods to realize described function to each certain applications, but this realization should not thought and exceeds scope of the present invention.
The software module that the method for describing in conjunction with embodiment disclosed herein or the step of algorithm can use hardware, processor to carry out, perhaps the combination of the two is implemented.Software module can place the storage medium of any other form known in random access memory (RAM), internal memory, ROM (read-only memory) (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or the technical field.
Although illustrated and described some embodiments of the present invention, it will be understood by those skilled in the art that under the situation that does not break away from principle of the present invention and spirit can carry out various modifications to these embodiment, such modification should fall in the scope of the present invention.

Claims (11)

1. the age estimation method in the face identification system is characterized in that, comprising:
Processing face images of users is to extract initial characteristic data;
The subclass of selecting said initial characteristic data is as low dimensional feature data, and the dimension of wherein said low dimensional feature data is less than the dimension of said initial characteristic data;
To the constituent analysis of being correlated with of said low dimensional feature data, to obtain training data;
According to the said training data that is obtained, training is used for the regretional analysis parameter of estimation of Age;
Utilize the said user's of said regretional analysis parameter estimation age,
It is wherein said that the constituent analysis of being correlated with comprises to obtain training data to said low dimensional feature data:
Construct at least one equivalent samples subclass of said low dimensional feature data, wherein each said equivalent samples subclass belongs to age-grade label;
Utilize the equivalent samples subclass of being constructed, said low dimensional feature data are differentiated dimensionality reduction, to obtain by the dimensionality reduction sample data of further dimensionality reduction;
Reduce the non-correlation of said dimensionality reduction sample data, to obtain said training data.
2. age estimation method as claimed in claim 1 is characterized in that, the subclass of the said initial characteristic data of said selection comprises as low dimensional feature data: the subclass that said initial characteristic data is selected in the employing principal component analysis (PCA) is as said low dimensional feature data.
3. age estimation method as claimed in claim 1 is characterized in that, the equivalent samples subclass that said utilization is constructed is differentiated dimensionality reduction to said low dimensional feature data and comprised:
Confirm the interior covariance matrix of class of said equivalent samples subclass;
As Estimation of covariance matrix value in the class of said low dimensional feature data, utilize combination Fisher discriminatory analysis that said low dimensional feature data are differentiated dimensionality reduction covariance matrix in the class of said equivalent samples subclass.
4. age estimation method as claimed in claim 1 is characterized in that, the non-correlation of the said dimensionality reduction sample data of said minimizing comprises:
For the equivalent samples subclass that belongs to age-grade label, defining the big direction of its variance is irrelevant direction.
5. age estimation method as claimed in claim 1 is characterized in that, the non-correlation of the said dimensionality reduction sample data of said minimizing comprises:
Confirm the interior covariance matrix of class of said dimensionality reduction sample data;
Covariance matrix in the class of said dimensionality reduction sample data is carried out leucismus change, obtain leucismus and change the result;
Change result and said dimensionality reduction sample data based on said leucismus, obtain said training data.
6. the estimation of Age equipment in the face identification system is characterized in that, comprising:
The primitive character extraction module is used for processing face images of users to extract initial characteristic data;
Feature selection module, the subclass that is used to select said initial characteristic data is as low dimensional feature data, and the dimension of wherein said low dimensional feature data is less than the dimension of said initial characteristic data;
Relevant component analysis module is used for to the constituent analysis of be correlated with of said low dimensional feature data, with the acquisition training data;
Training module is used for training the regretional analysis parameter that is used for estimation of Age according to by the said training data that said relevant component analysis module obtained;
Estimation module is used to utilize said regretional analysis parameter to estimate said user's age,
Wherein said relevant component analysis module comprises:
Tectonic element is used to construct at least one equivalent samples subclass of said low dimensional feature data, and wherein each said equivalent samples subclass belongs to age-grade label;
Differentiate the dimensionality reduction unit, be used to utilize the equivalent samples subclass of being constructed, said low dimensional feature data are differentiated dimensionality reduction, to obtain by the dimensionality reduction sample data of further dimensionality reduction;
The decorrelation unit is used to reduce the non-correlation of said dimensionality reduction sample data, to obtain said training data.
7. estimation of Age equipment as claimed in claim 6 is characterized in that, the subclass that said initial characteristic data is selected in said feature selection module employing principal component analysis (PCA) is as said low dimensional feature data.
8. estimation of Age equipment as claimed in claim 6; It is characterized in that; The interior covariance matrix of the class of said equivalent samples subclass is confirmed in said differentiation dimensionality reduction unit; As Estimation of covariance matrix value in the class of said low dimensional feature data, utilize combination Fisher discriminatory analysis that said low dimensional feature data are differentiated dimensionality reduction covariance matrix in the class of said equivalent samples subclass.
9. estimation of Age equipment as claimed in claim 6 is characterized in that, said decorrelation module is for the equivalent samples subclass that belongs to age-grade label, and defining the big direction of its variance is irrelevant direction.
10. estimation of Age equipment as claimed in claim 6; It is characterized in that; Said decorrelation module is confirmed the interior covariance matrix of the class of said dimensionality reduction sample data, said dimensionality reduction sample data is carried out leucismus change, to obtain said training data; Change result and said dimensionality reduction sample data based on said leucismus, obtain said training data.
11. a face identification system is characterized in that, comprises like each described estimation of Age equipment of claim 6 to 10.
CN2009101310593A 2009-04-21 2009-04-21 Age estimation method, equipment and face recognition system Expired - Fee Related CN101615248B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101310593A CN101615248B (en) 2009-04-21 2009-04-21 Age estimation method, equipment and face recognition system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101310593A CN101615248B (en) 2009-04-21 2009-04-21 Age estimation method, equipment and face recognition system

Publications (2)

Publication Number Publication Date
CN101615248A CN101615248A (en) 2009-12-30
CN101615248B true CN101615248B (en) 2012-03-21

Family

ID=41494880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101310593A Expired - Fee Related CN101615248B (en) 2009-04-21 2009-04-21 Age estimation method, equipment and face recognition system

Country Status (1)

Country Link
CN (1) CN101615248B (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5827225B2 (en) * 2010-06-21 2015-12-02 ポーラ化成工業株式会社 Age estimation method
CN102567719B (en) * 2011-12-26 2014-07-02 东南大学 Human age automatic estimation method based on posterior probability neural network
CN103188427B (en) * 2011-12-30 2016-08-10 华晶科技股份有限公司 Image capture unit and the control method thereof of image feature value group can be simplified
CN103824054B (en) * 2014-02-17 2018-08-07 北京旷视科技有限公司 A kind of face character recognition methods based on cascade deep neural network
TWI716344B (en) * 2014-02-24 2021-01-21 日商花王股份有限公司 Aging analyzing method, aging care counselling method using aging analyzing method, aging analyzing device and computer readable recording medium
CN104463190B (en) * 2014-10-30 2017-08-29 华为技术有限公司 Age estimation method and equipment
CN104573652B (en) 2015-01-04 2017-12-22 华为技术有限公司 Determine the method, apparatus and terminal of the identity of face in facial image
CN106529378B (en) * 2015-09-15 2019-04-02 中国科学院声学研究所 A kind of the age characteristics model generating method and age estimation method of asian ancestry's face
CN106529377A (en) * 2015-09-15 2017-03-22 北京文安智能技术股份有限公司 Age estimating method, age estimating device and age estimating system based on image
CN105550641B (en) * 2015-12-04 2020-03-31 康佳集团股份有限公司 Age estimation method and system based on multi-scale linear differential texture features
CN106056083B (en) * 2016-05-31 2019-08-13 腾讯科技(深圳)有限公司 A kind of information processing method and terminal
CN106200432A (en) * 2016-07-29 2016-12-07 深圳前海弘稼科技有限公司 Control method and device, intelligence plantation equipment for intelligence plantation equipment
CN108022116B (en) * 2016-11-01 2021-06-29 北京京东尚科信息技术有限公司 Method, system and terminal equipment for modeling user
CN107580262B (en) * 2017-09-25 2018-05-04 湖北工业大学 Smart television Password Input sets system
CN108460630B (en) * 2018-02-12 2021-11-02 广州虎牙信息科技有限公司 Method and device for carrying out classification analysis based on user data

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王蕴红等.融合全局与局部特征的子空间人连识别算法.《计算机学报》.2005,第28卷(第10期),1657-1663. *
郑成勇.基于小波变换和回归分析的人脸识别.《五邑大学学报(自然科学版)》.2007,第21卷(第1期),16-20. *

Also Published As

Publication number Publication date
CN101615248A (en) 2009-12-30

Similar Documents

Publication Publication Date Title
CN101615248B (en) Age estimation method, equipment and face recognition system
Wang et al. Manifold–manifold distance and its application to face recognition with image sets
CN109871860B (en) Daily load curve dimension reduction clustering method based on kernel principal component analysis
Räsänen et al. Data-based method for creating electricity use load profiles using large amount of customer-specific hourly measured electricity use data
He et al. A two-stage genetic algorithm for automatic clustering
Räsänen et al. Feature-based clustering for electricity use time series data
Marzban et al. Cluster analysis for verification of precipitation fields
CN103065160B (en) Based on the hyperspectral image classification method that the collaborative expression in local and neighborhood information retrain
D’Urso et al. Autoregressive model-based fuzzy clustering and its application for detecting information redundancy in air pollution monitoring networks
CN101807258B (en) SAR (Synthetic Aperture Radar) image target recognizing method based on nuclear scale tangent dimensionality reduction
CN111160401A (en) Abnormal electricity utilization judging method based on mean shift and XGboost
Gu et al. Local learning regularized nonnegative matrix factorization
CN101515328B (en) Local projection preserving method for identification of statistical noncorrelation
Cziráky et al. Regional development assessment: A structural equation approach
Chen et al. Learning spatial-temporal varying graphs with applications to climate data analysis
CN102411610A (en) Semi-supervised dimensionality reduction method for high dimensional data clustering
CN108596227B (en) Mining method for dominant influence factors of electricity consumption behaviors of users
CN102722578B (en) Unsupervised cluster characteristic selection method based on Laplace regularization
Li et al. Cluster-based data filtering for manufacturing big data systems
CN106778714A (en) LDA face identification methods based on nonlinear characteristic and model combination
Pérez et al. Interactive feature space extension for multidimensional data projection
Lu et al. A hybrid sales forecasting scheme by combining independent component analysis with K-means clustering and support vector regression
Orazi et al. A regional Kohonen map of financial inclusion and related macroeconomic variables
CN104537377A (en) Image data dimension reduction method based on two-dimensional kernel entropy component analysis
CN104573727A (en) Dimension reduction method of handwritten digital image

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120321

Termination date: 20180421