CN102915527A - Face image super-resolution reconstruction method based on morphological component analysis - Google Patents

Face image super-resolution reconstruction method based on morphological component analysis Download PDF

Info

Publication number
CN102915527A
CN102915527A CN2012103897128A CN201210389712A CN102915527A CN 102915527 A CN102915527 A CN 102915527A CN 2012103897128 A CN2012103897128 A CN 2012103897128A CN 201210389712 A CN201210389712 A CN 201210389712A CN 102915527 A CN102915527 A CN 102915527A
Authority
CN
China
Prior art keywords
resolution
image
low
face
residual error
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
CN2012103897128A
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.)
Sun Yat Sen University
National Sun Yat Sen University
Original Assignee
National Sun Yat Sen 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 National Sun Yat Sen University filed Critical National Sun Yat Sen University
Priority to CN2012103897128A priority Critical patent/CN102915527A/en
Publication of CN102915527A publication Critical patent/CN102915527A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

The invention discloses a face image super-resolution reconstruction method based on morphological component analysis. The method comprises the following steps: upsampling a low-resolution input image to acquire an interpolated image; acquiring a whole high-resolution face image from the interpolated image by a morphological component analysis method; downsampling the whole high-resolution face image and subtracting the whole high-resolution face image from the low-resolution input image to acquire a low-resolution residual image, and performing face detail information compensation by performing neighboring construction on an image block at each face position to acquire a high-resolution residual image; and combining the high-resolution residual image with the whole high-resolution face image to acquire the face image super-resolution result and finish resolution reconstruction. The invention also provides a method for performing face image super-resolution reconstruction and expression normalization simultaneously on the basis of morphological component analysis. By the method, detailed face detail information can be acquired, the ringing phenomenon is eliminated and high image quality is achieved.

Description

Face image super-resolution reconstruction method based on the morphology constituent analysis
Technical field
The present invention relates to image processing field, particularly a kind of face image super-resolution reconstruction method based on the morphology constituent analysis.
Background technology
The task of image super-resolution is to infer high-definition picture from one or more low resolution input picture.Be widely used in real life, particularly the Super-resolution reconstruction of facial image is built the aspects such as long-distance video monitoring or Video processing in important application.
Traditional face image super-resolution method combines the peculiar classical priori of facial image, and a lot of face image super-resolution methods are appearring closely during the decade, as make up based on differential information radio-frequency component, based on principal component analysis (PCA) (PCA) set up block mold reconstruct high-resolution human face image, in conjunction with based on the block mold of PCA and the methods such as " whole faces+residual error face " that makes up based on the partial model of belief propagation.Although existing face image super-resolution reconstruction method has been obtained good effect, but still have greatly improved the space, a part that wherein needs to promote is whole face feature, existing algorithm often causes people's face of reconstruct and real people's face unlike same person, and generation ringing effect, this is because existing estimation model only to global face rather than all the components modeling of input picture, has therefore reduced the degree of accuracy of estimating.
Summary of the invention
Fundamental purpose of the present invention is to overcome the shortcoming of prior art with not enough, face image super-resolution reconstruction method based on the morphology constituent analysis is provided, people's face that the method is rebuild is more as real people's face, can obtain more detailed facial detail information, eliminate ringing, obtain better image quality.In addition, the present invention can also be extended to and solve simultaneously face image super-resolution and the normalized problem of expression.
Purpose of the present invention realizes by following technical scheme: the face image super-resolution reconstruction method based on the morphology constituent analysis may further comprise the steps:
(1) by interpolation method low resolution input picture up-sampling, obtain the image after the interpolation;
(2) image after using the morphology component analyzing method to the resulting interpolation of step (1) carries out picture breakdown, resolves into whole high-resolution human face image and anti-sharpening mask;
(3) the whole high-resolution human face image that step (2) is obtained carries out down-sampling, then reduced by the low resolution input picture, obtain the low-resolution residual error image, then come image is carried out the facial detail information compensation by the image block on the face location of low-resolution residual error image everywhere is carried out neighbour's reconstruct, obtain the high-resolution residual error image;
(4) the whole high-resolution human face image that the high-resolution residual error image in the step (3) and step (2) is obtained merges, and obtains final face image super-resolution result, finishes super-resolution rebuilding.
Up-sampling Interpolation Process formula in the described step (1) is: x Zoom=s ↑ x l=s ↑ s ↓ (h*x h) ≈ h*x h, s ↑ be up-sampling interpolation operation symbol wherein, s ↓ be the down-sampling operational character, x lThe low resolution input picture, x ZoomThe image after the interpolation, x hBe true high-definition picture, h is the some propagator of camera, and * is the convolution operation symbol; Image x after the interpolation ZoomWith equation expression be: x Zoom=x h+ x s, x wherein sBe comprise on the occasion of with the radio-frequency component of negative value, be called anti-sharpening mask, the image x after the interpolation ZoomX hVersion after fuzzy.
Preferably, the step of carrying out picture breakdown of the use morphology component analyzing method in the described step (2) specifically refers to:
(2-1) choose at random different people's face high-definition pictures, will obtain people's face low-resolution image behind its down-sampling, obtain image after the interpolation by interpolation method again, the image after the interpolation and corresponding high-definition picture are subtracted each other, obtain anti-sharpening mask; With high-definition picture and anti-sharpening mask as training sample, by the average of principal component analysis (PCA) (Principal Component Analysis, PCA) methodology acquistion to true high-definition picture
Figure BDA00002254016400021
Average with anti-sharpening mask
Figure BDA00002254016400022
Obtain all N eigenwert of corresponding true high-definition picture
Figure BDA00002254016400023
All N eigenwert with anti-sharpening mask
Figure BDA00002254016400024
And the eigenmatrix P of corresponding true high-definition picture hEigenmatrix P with anti-sharpening mask s, x then hAnd x sBe expressed as respectively:
Figure BDA00002254016400025
α wherein hAnd α sProjection coefficient, P hAnd P sBeing used as dictionary describes respectively
Figure BDA00002254016400027
With
Figure BDA00002254016400028
(2-2) based on the morphology component analyzing method, its core is to find the solution an optimized process, and the formula of finding the solution is min α h , α s | | α h | | 0 + | | α s | | 0 , Constraint condition is:
Figure BDA00002254016400031
Wherein, ε is the fault tolerant rate, and (m, n) is pixel coordinate, and S is enlargement factor, x lIt is the low resolution input picture; Finally, by optimum solution
Figure BDA00002254016400032
Produce the subsignal of separation
Figure BDA00002254016400033
With
Figure BDA00002254016400034
Above-mentioned solution formula is a NP Solve problems, and its sparse solution can approach with the solution of L1 Norm minimum, can find the solution by using alternately the method for iteration threshold, and each iteration all will be the iteration result super ellipsoid that stretches back
Figure BDA00002254016400035
Within the scope, and guarantee Reconstruction Constraints, namely input low-resolution image x lWith high-definition picture x hRelation need to satisfy And threshold value and iterations can obtain suitable value by experiment.Solution procedure is specific as follows:
(2-2-1) initialization: input maximum iteration time I MaxInput threshold value maxT h, maxT s, minTh and minT s, these threshold values are: max T h = | | P h T ( x zoom - x h ‾ - x s ‾ ) | | ∞ ,
Figure BDA00002254016400039
MinT hAnd minT sInput value is the square root of t eigenvalue of maximum corresponding to the corresponding true high-definition picture that uses the PCA method and obtain and t major component of anti-sharpening mask; Input back projection operational symbol p; Set x h, x sInitial value is 0, T h=maxT h, T s=maxT sSet iteration variable k, allow k get I from l Max, step (2-2-2) is the process of iteration;
(2-2-2) keep x sConstant calculating x hThe process of new value:
Calculate first residual error r = x zoom - x h ‾ - x s ‾ - x h - x s ,
Then calculate x h+ r is to P hProjection coefficient:
Passing threshold T hTo parameter alpha hCarrying out the hard-threshold processing obtains Convergent-divergent Obtain to satisfy constraint condition
Figure BDA000022540164000314
α sGet in the last iterative process and upgrade x sObtain
Figure BDA000022540164000315
Calculate, constraint condition is as follows:
| | x zoom - x h ‾ - x s ‾ - P h α h - P s α s | | 2 ≤ ϵ ,
Figure BDA00002254016400042
Figure BDA00002254016400043
i=1…N,
x l ( m , n ) = 1 S 2 Σ i = 1 S Σ j = 1 S x h ( Sm + i , Sn + j ) . ;
Reconstruct x h = x h ‾ + P h α h 2 ;
Back projection correction: with x hReplace with
Figure BDA00002254016400046
Then, keep x hConstant calculating x sNew value process is as follows:
Calculate first residual error
Then calculate x s+ r is to P sProjection coefficient:
Figure BDA00002254016400048
Passing threshold T sTo parameter alpha sCarrying out the hard-threshold processing obtains
Figure BDA00002254016400049
Convergent-divergent
Figure BDA000022540164000410
Obtain to satisfy constraint condition
Figure BDA000022540164000411
α hCalculate with front portion Constraint condition is as follows:
| | x zoom - x h ‾ - x s ‾ - P h α h - P s α s | | 2 ≤ ϵ ,
Figure BDA000022540164000415
i=1…N,
x l ( m , n ) = 1 S 2 Σ i = 1 S Σ j = 1 S x h ( Sm + i , Sn + j ) . ;
Reconstruct x s = x s ‾ + P s α s 2 ;
Upgrade threshold value T h = T h - max T h - min T h I max - 1 , T s = T s - max T s - min T s I max - 1 ;
When iterations then withdraws from iterative program and obtains optimum solution greater than the maximum iteration time that arranges
Figure BDA000022540164000420
Produce the subsignal of separation
Figure BDA000022540164000421
With
Figure BDA000022540164000422
Concrete, in the described step (3), the concrete steps that obtain the low-resolution residual error image are: the whole high-resolution human face image that is obtained by step (2)
Figure BDA00002254016400051
With low resolution input picture x l, definition low-resolution residual error image is:
Figure BDA00002254016400052
By based on neighbour's reconstruct of the image block on the face location of everywhere from
Figure BDA00002254016400053
Set up the high-resolution residual error image
Figure BDA00002254016400054
Estimation.
Concrete, described step (3) is to neighbour's reconstructing method of the image block on the face location of everywhere, and is specific as follows:
Obtain people's face low-resolution image as training sample behind the people's face high-definition picture that (3-1) will choose at random and its down-sampling.Low-resolution image applying step (2) is obtained whole high-resolution human face image, then by the subduction of people's face full resolution pricture, obtain high-resolution residual error image training sample.Obtained low-resolution residual error image training sample by the subduction of people's face low-resolution image behind the whole high-resolution human face image down sampling.
(3-2) with the low-resolution residual error image of inputting
Figure BDA00002254016400055
Be divided into image block, make between the image block to overlap.
(3-3) for the low-resolution residual error image block of each input, k arest neighbors low-resolution residual error image block training sample by same face location carries out portfolio restructuring, then the weights that keeping combination synthesize corresponding high-resolution residual error image block by substituting low-resolution residual error image block training sample for corresponding high-resolution residual error image block training sample.
(3-4) obtain complete high-resolution residual error image by the integrated high-resolution residual error image block of averaged overlapping portions
Figure BDA00002254016400056
As another kind of preferred version, based on the Multichannel Decomposition ability of morphology component analyzing method, to the expression shape change modeling in the morphology constituent analysis, expansion algorithm is to carrying out simultaneously face image super-resolution and expression normalization by the design dictionary.Face image super-resolution reconstruction method based on the morphology constituent analysis may further comprise the steps:
(1) study obtains the mutual incoherent expression dictionary of several correspondences to some classes expression by principal component analytical method;
(2) for the low resolution input picture of input, judge the expression that it is affiliated.Then, on the one hand carry out up-sampling, obtain the image of the espressiove after the interpolation; Use on the other hand the morphology component analyzing method, based on its corresponding expression dictionary of expressing one's feelings, the low resolution input picture is carried out picture breakdown, resolve into the amimia input picture of whole low resolution and low resolution espressiove mask;
(3) image of the espressiove after using the morphology component analyzing method to the resulting interpolation of step (2) carries out picture breakdown, based on its corresponding expression dictionary of expressing one's feelings, resolves into the amimia image of whole high resolving power, expression mask and anti-sharpening mask;
(4) the amimia image of the resulting whole high resolving power of step (3) is carried out down-sampling, then and the amimia input picture of whole low resolution that obtains of step (2) subtract each other, obtain the low-resolution residual error image; To this low-resolution residual error image, come image is carried out the facial detail information compensation by the image block on the face location of everywhere being carried out neighbour's reconstruct, obtain the high-resolution residual error image;
(5) the amimia image of whole high resolving power that the high-resolution residual error image in the step (4) and step (3) is obtained merges, and obtains final face image super-resolution result, finishes super-resolution rebuilding.
Further, in the described step (2), by adopting LBP(Local Binary Patterns, local binary patterns) sorter judges the expression in the low resolution input picture.
As preferably, face image super-resolution reconstruction method based on the morphology constituent analysis, all low resolution input pictures are all judged whether espressiove of its people's face by sorter before processing, if do not have, then adopt first method to process, if espressiove then adopts second method to process.
The present invention compared with prior art has following advantage and beneficial effect:
1, the present invention is interpreted as the picture breakdown problem to the face image super-resolution problem, and the framework that uses the morphology component analyzing method to set up one three step has been realized the decomposition of image effectively, experiment shows, method of the present invention is more effective than existing method, can obtain more detailed facial detail information, eliminated ringing, obtained better image quality, and under the scene of laboratory environment and reality, can both realize satisfactory results.
2, the invention allows for a kind of method, be extended to and solve simultaneously face image super-resolution and the normalized problem of expression, this is first successful method in the existing method.Thereby can further improve the accuracy of facial image identification.
Description of drawings
Fig. 1 is the algorithm flow chart of the embodiment of the invention 1;
Fig. 2 is the algorithm flow chart of the embodiment of the invention 2.
Embodiment
The present invention is described in further detail below in conjunction with embodiment and accompanying drawing, but embodiments of the present invention are not limited to this.
Embodiment 1
As shown in Figure 1, the face image super-resolution reconstruction method based on the morphology constituent analysis may further comprise the steps:
(1) by interpolation method low resolution input picture up-sampling, obtain the image after the interpolation;
(2) image after using the morphology component analyzing method to the resulting interpolation of step (1) carries out picture breakdown, resolves into whole high-resolution human face image and anti-sharpening mask;
(3) the whole high-resolution human face image that step (2) is obtained carries out down-sampling, then reduced by the low resolution input picture, obtain the low-resolution residual error image, then come image is carried out the facial detail information compensation by the image block on the face location of low-resolution residual error image everywhere is carried out neighbour's reconstruct, obtain the high-resolution residual error image;
(4) the whole high-resolution human face image that the high-resolution residual error image in the step (3) and step (2) is obtained merges, and obtains final face image super-resolution result, finishes super-resolution rebuilding.
Up-sampling Interpolation Process formula in the described step (1) is: x Zoom=s ↑ x l=s ↑ s ↓ (h*x h) ≈ h*x h, s ↑ be up-sampling interpolation operation symbol wherein, s ↓ be the down-sampling operational character, x lThe low resolution input picture, x ZoomThe image after the interpolation, x hBe true high-definition picture, h is the some propagator of camera, and * is the convolution operation symbol; Image x after the interpolation ZoomWith equation expression be: x Zoom=x h+ x s, x wherein sBe comprise on the occasion of with the radio-frequency component of negative value, be called anti-sharpening mask, the image x after the interpolation ZoomX hVersion after fuzzy.
The step that use morphology component analyzing method in the described step (2) carries out picture breakdown specifically refers to:
(2-1) choose at random different people's face high-definition pictures, will obtain people's face low-resolution image behind its down-sampling, obtain image after the interpolation by interpolation method again, the image after the interpolation and corresponding high-definition picture are subtracted each other, obtain anti-sharpening mask; High-definition picture and anti-sharpening mask as training sample, are obtained the average of true high-definition picture by principal component analytical method study
Figure BDA00002254016400071
Average with anti-sharpening mask Obtain all N eigenwert of corresponding true high-definition picture
Figure BDA00002254016400073
All N eigenwert with anti-sharpening mask
Figure BDA00002254016400074
And the eigenmatrix P of true high-definition picture hEigenmatrix P with anti-sharpening mask s, x then hAnd x sBe expressed as respectively:
Figure BDA00002254016400075
Figure BDA00002254016400076
α wherein hAnd α sProjection coefficient, P hAnd P sJust being used as dictionary describes respectively
Figure BDA00002254016400077
With
Figure BDA00002254016400078
(2-2) based on the morphology component analyzing method, solution formula is:
min α h , α s | | α h | | 0 + | | α s | | 0 ,
Constraint condition is as follows:
| | x zoom - x h ‾ - x s ‾ - P h α h - P s α s | | 2 ≤ ϵ ,
Figure BDA00002254016400082
Figure BDA00002254016400083
i=1…N,
x l ( m , n ) = 1 S 2 Σ i = 1 S Σ j = 1 S x h ( Sm + i , Sn + j ) .
Wherein, ε is the fault tolerant rate, and (m, n) is pixel coordinate, and S is enlargement factor, x lIt is the low resolution input picture; Finally, by optimum solution
Figure BDA00002254016400085
Produce the subsignal of separation
Figure BDA00002254016400086
With
Figure BDA00002254016400087
Above-mentioned solution formula is a NP Solve problems, and its sparse solution can approach with the solution of L1 Norm minimum, can find the solution by using alternately the method for iteration threshold, and each iteration all will be the iteration result super ellipsoid that stretches back
Figure BDA00002254016400088
Within the scope, and guarantee Reconstruction Constraints, namely input low-resolution image x lWith high-definition picture x hRelation need to satisfy
Figure BDA00002254016400089
And threshold value and iterations can obtain suitable value by experiment.
In the described step (3), the concrete steps that obtain the low-resolution residual error image are: the whole high-resolution human face image that is obtained by step (2)
Figure BDA000022540164000810
With low resolution input picture x l, definition low-resolution residual error image is:
Figure BDA000022540164000811
By based on neighbour's reconstruct of the image block on the face location of everywhere from
Figure BDA000022540164000812
Set up the high-resolution residual error image
Figure BDA000022540164000813
Estimation.
Described step (3) is to neighbour's reconstructing method of the image block on the face location of everywhere, and is specific as follows:
Obtain people's face low-resolution image as training sample behind the people's face high-definition picture that (3-1) will choose at random and its down-sampling.Low-resolution image applying step (2) is obtained whole high-resolution human face image, then by the subduction of people's face full resolution pricture, obtain high-resolution residual error image training sample.Obtained low-resolution residual error image training sample by the subduction of people's face low-resolution image behind the whole high-resolution human face image down sampling.
(3-2) with the low-resolution residual error image of inputting
Figure BDA000022540164000814
Be divided into image block, make between the image block to overlap.
(3-3) for the low-resolution residual error image block of each input, k arest neighbors low-resolution residual error image block training sample by same face location carries out portfolio restructuring, then the weights that keeping combination synthesize corresponding high-resolution residual error image block by substituting low-resolution residual error image block training sample for corresponding high-resolution residual error image block training sample.
(3-4) obtain complete high-resolution residual error image by the integrated high-resolution residual error image block of averaged overlapping portions
Figure BDA00002254016400091
Embodiment 2
Present embodiment except following characteristics other structures with embodiment 1: based on the Multichannel Decomposition ability of morphology component analyzing method, to the expression shape change modeling in the morphology constituent analysis, expansion algorithm is to carrying out simultaneously face image super-resolution and expression normalization by the design dictionary.Face image super-resolution reconstruction method based on the morphology constituent analysis may further comprise the steps:
(1) set five kinds of basic human face expressions: close one's eyes, frown, smile, surprised, open one's mouth.Respectively five kinds of expression study are obtained mutual incoherent expression dictionary by principal component analytical method, specific as follows:
(1-1) to each expression, choose at random different people face high-definition picture and the corresponding amimia high-definition picture of this expression.People's face low-resolution image and the corresponding amimia low-resolution image of this expression will be obtained behind its down-sampling.Obtain image after the interpolation by interpolation method again.With this facial expression image after the interpolation and corresponding amimia image subtraction, obtain high resolving power expression mask.With this facial expression image of low resolution and corresponding amimia image subtraction, obtain low resolution expression mask.
(1-2) with the expression mask of high-resolution and low-resolution as training sample, learn respectively to obtain expression average, expressive features value and the corresponding expression dictionary of high-resolution and low-resolution by principal component analytical method.
(2) for the low resolution input picture of input, judge the expression that it is affiliated.Then, on the one hand carry out up-sampling, obtain the image of the espressiove after the interpolation; Use on the other hand the morphology component analyzing method, based on its express one's feelings corresponding low resolution expression average, expressive features value and expression dictionary, the low resolution input picture is carried out picture breakdown, resolve into the amimia input picture of whole low resolution With low resolution espressiove mask x Le
(3) because the facial expression image after the interpolation
Figure BDA00002254016400093
Can be expressed as the amimia image x after the interpolation ZoomWith high resolving power expression mask x eAnd:
Figure BDA00002254016400094
And the amimia image after the interpolation can further be expressed as x Zoom=x h+ x s, then Therefore, the image of the espressiove after using the morphology component analyzing method to the resulting interpolation of step (2) carries out picture breakdown, based on its express one's feelings corresponding high resolving power expression average, expressive features value and expression dictionary, resolve into the amimia image of whole high resolving power The expression mask
Figure BDA00002254016400097
With anti-sharpening mask
(4) the amimia image of the resulting whole high resolving power of step (3) is carried out down-sampling, then and the amimia input picture of whole low resolution that obtains of step (2) subtract each other, obtain the low-resolution residual error image
Figure BDA00002254016400101
To this low-resolution residual error image, come image is carried out the facial detail information compensation by the image block on the face location of everywhere being carried out neighbour's reconstruct, obtain the high-resolution residual error image
Figure BDA00002254016400102
(5) the amimia image of whole high resolving power that the high-resolution residual error image in the step (4) and step (3) is obtained merges, and obtains final face image super-resolution result
Figure BDA00002254016400103
Finish super-resolution rebuilding.
In the described step (2), judge expression in the low resolution input picture by adopting the LBP sorter.
The same on the process of using the morphology component analyzing method in the described step (1), step (3) and the Method And Principle of embodiment 1, just embodiment 1 becomes 2 parts to picture breakdown, and embodiment 2 has added the expression mask, therefore need to become 3 parts to picture breakdown.
Embodiment 3
Present embodiment except following characteristics other structures with embodiment 1: in the present embodiment, face image super-resolution reconstruction method based on the morphology constituent analysis, all low resolution input pictures are all judged whether espressiove of its people's face by sorter before processing, if do not have, then adopt embodiment 1 described method to process, if espressiove then adopts embodiment 2 described methods to process.
The present invention describes effect of the present invention by following experiment: experimental selection CES-PEAL-R1 face data storehouse, this database comprises 99594 images of 1040 people, all image scalings on the storehouse are done simple calibration to the 128*96 pixel and according to the position of eyes and face, as shown in table 1, the subset of 5 positive face images is used in two experiments, and existing certain methods also has been implemented with as a comparison:
The content of form 1.CAS-PEAL-R1 database 5 subsets
First experiment is the visual quality evaluation experimental.The common collection of CAS-PEAL-R1 is used in this experiment, common concentrated each personage comprises 1040 images, with amplification factor 4 (or 8) original high-definition picture is carried out the level and smooth and down-sampled low-resolution image that obtains corresponding 32*24 (perhaps 16*12), and they are called high resolving power-low-resolution image pair.For each amplification factor, select 1000 high resolving power-low-resolution images to as training set, remaining 40 images are to as test set, in particular for the framework that uses " whole face+residual error face ", with 700 images to training block mold, with 300 images to training the residual compensation model.For morphology constituent analysis parameter, maximum threshold value is Lookup protocol as follows: max T h = | | P h T ( x zoom - x h ‾ - x s ‾ ) | | ∞ With
Figure BDA00002254016400111
Threshold value minT hAnd minT sBe the square root that uses t eigenvalue of maximum corresponding to corresponding true high-definition picture that the PCA method obtains and t major component of anti-sharpening mask, algorithm can obtain stable output when the observation number of times that can obtain working as iteration was above 50 times by experiment.For the quantized image vision quality, calculate the difference that 2 indexs are measured image and the true high-definition picture of reconstructing human face super resolution, one of them is square error (MSE), formula is
MSE = Σ i = 1 N | | I test , i h - J test , i h | | 2 2 / ( WHN ) Wherein
Figure BDA00002254016400113
True high-definition picture,
Figure BDA00002254016400114
Be the facial image of reconstruct, it is picture traverse and height that W and H differentiate, and N is the quantity of test pattern, and another one is structural similarity index (SSIM).Present embodiment institute employing method is that the MSE value that obtained in 4 o'clock is that 85, SSIM value is 0.9 in amplification factor, is that the MSE value that obtained in 8 o'clock is that 260, SSIM value is 0.79 in amplification factor, has good image quality assessment compared to more existing algorithms.
Second experiment is super-resolution reconstruction and expression normalization experiment: select 240 images as test set from the expression subset, remaining image is used as training set in the subset, and other parameters of experiment arrange consistent with experiment before.The result shows that the method for present embodiment can obtain expressionless high-definition picture from the low resolution input picture of multiple expression, and the recognition of face rate is higher by 15% than the method for bilinear interpolation, and is higher by 5% than the method for using original high-definition picture with multiple expression.
Above-described embodiment is the better embodiment of the present invention; but embodiments of the present invention are not restricted to the described embodiments; other any do not deviate from change, the modification done under Spirit Essence of the present invention and the principle, substitutes, combination, simplify; all should be the substitute mode of equivalence, be included within protection scope of the present invention.

Claims (7)

1. based on the face image super-resolution reconstruction method of morphology constituent analysis, it is characterized in that, may further comprise the steps:
(1) by interpolation method low resolution input picture up-sampling, obtain the image after the interpolation;
(2) image after using the morphology component analyzing method to the resulting interpolation of step (1) carries out picture breakdown, resolves into whole high-resolution human face image and anti-sharpening mask;
(3) the whole high-resolution human face image that step (2) is obtained carries out down-sampling, then reduced by the low resolution input picture, obtain the low-resolution residual error image, then come image is carried out the facial detail information compensation by the image block on the face location of low-resolution residual error image everywhere is carried out neighbour's reconstruct, obtain the high-resolution residual error image;
(4) the whole high-resolution human face image that the high-resolution residual error image in the step (3) and step (2) is obtained merges, and obtains final face image super-resolution result, finishes super-resolution rebuilding.
2. the face image super-resolution reconstruction method based on the morphology constituent analysis according to claim 1 is characterized in that, the up-sampling Interpolation Process formula in the described step (1) is: x Zoom=s ↑ x i=s ↑ s ↓ (h*x h) ≈ h*x h, s ↑ be up-sampling interpolation operation symbol wherein, s ↓ be the down-sampling operational character, x lThe low resolution input picture, x ZoomThe image after the interpolation, x hBe true high-definition picture, h is the some propagator of camera, and * is the convolution operation symbol; Image x after the interpolation ZoomWith equation expression be: x Zoom=x h+ x s, x wherein sBe comprise on the occasion of with the radio-frequency component of negative value, be called anti-sharpening mask, the image x after the interpolation ZoomX hVersion after fuzzy.
3. the face image super-resolution reconstruction method based on the morphology constituent analysis according to claim 1 is characterized in that, the step that the use morphology component analyzing method in the described step (2) carries out picture breakdown specifically refers to:
(2-1) choose at random different people's face high-definition pictures, will obtain people's face low-resolution image behind its down-sampling, obtain image after the interpolation by interpolation method again, the image after the interpolation and corresponding high-definition picture are subtracted each other, obtain anti-sharpening mask; High-definition picture and anti-sharpening mask as training sample, are obtained the average of true high-definition picture by principal component analytical method study Average with anti-sharpening mask
Figure FDA00002254016300012
Obtain all N eigenwert of corresponding true high-definition picture
Figure FDA00002254016300013
All N eigenwert with anti-sharpening mask
Figure FDA00002254016300014
And the eigenmatrix P of corresponding true high-definition picture hEigenmatrix P with anti-sharpening mask s, x then hAnd x sBe expressed as respectively:
Figure FDA00002254016300015
α wherein hAnd α sProjection coefficient, P hAnd P sBeing used as dictionary describes respectively
Figure FDA00002254016300017
With
Figure FDA00002254016300018
(2-2) based on the morphology component analyzing method, solution formula is:
min α h , α s | | α h | | 0 + | | α s | | 0 ,
Constraint condition is as follows:
| | x zoom - x h ‾ - x s ‾ - P h α h - P s α s | | 2 ≤ ϵ ,
Figure FDA00002254016300023
Figure FDA00002254016300024
i=1…N,
x l ( m , n ) = 1 S 2 Σ i = 1 S Σ j = 1 S x h ( Sm + i , Sn + j ) ;
Wherein, ε is the fault tolerant rate, and (m, n) is pixel coordinate, and S is enlargement factor, x lIt is the low resolution input picture; Finally, produce the subsignal of separation
Figure FDA00002254016300026
With
Figure FDA00002254016300027
Figure FDA00002254016300028
Be whole high-resolution human face image.
4. the face image super-resolution reconstruction method based on the morphology constituent analysis according to claim 1, it is characterized in that, in the described step (3), the concrete steps that obtain the low-resolution residual error image are: the whole high-resolution human face image that is obtained by step (2)
Figure FDA00002254016300029
With low resolution input picture x l, definition low-resolution residual error image is: x l l = x l - 1 S 2 Σ i = 1 S Σ j = 1 S x h g ( Sm + i , Sn + j ) , By based on neighbour's reconstruct of the image block on the face location of everywhere from
Figure FDA000022540163000211
Set up the high-resolution residual error image
Figure FDA000022540163000212
Estimation.
5. the face image super-resolution reconstruction method based on the morphology constituent analysis according to claim 4 is characterized in that, described step (3) is to neighbour's reconstructing method of the image block on the face location of everywhere, and is specific as follows:
Obtain people's face low-resolution image as training sample behind the people's face high-definition picture that (3-1) will choose at random and its down-sampling; Low-resolution image applying step (2) is obtained whole high-resolution human face image, then by the subduction of people's face full resolution pricture, obtain high-resolution residual error image training sample; Obtained low-resolution residual error image training sample by the subduction of people's face low-resolution image behind the whole high-resolution human face image down sampling;
(3-2) with the low-resolution residual error image of inputting
Figure FDA000022540163000213
Be divided into image block, make between the image block to overlap;
(3-3) for the low-resolution residual error image block of each input, k arest neighbors low-resolution residual error image block training sample by same face location carries out portfolio restructuring, then the weights that keeping combination synthesize corresponding high-resolution residual error image block by substituting low-resolution residual error image block training sample for corresponding high-resolution residual error image block training sample;
(3-4) obtain complete high-resolution residual error image by the integrated high-resolution residual error image block of averaged overlapping portions
Figure FDA00002254016300031
6. based on the face image super-resolution reconstruction method of morphology constituent analysis, it is characterized in that, may further comprise the steps:
(1) study obtains the mutual incoherent expression dictionary of several correspondences to some classes expression by principal component analytical method;
(2) for the low resolution input picture of input, judge the expression under it, then, carry out up-sampling on the one hand, obtain the image of the espressiove after the interpolation; Use on the other hand the morphology component analyzing method, based on its corresponding expression dictionary of expressing one's feelings, the low resolution input picture is carried out picture breakdown, resolve into the amimia input picture of whole low resolution and low resolution espressiove mask;
(3) image of the espressiove after using the morphology component analyzing method to the resulting interpolation of step (2) carries out picture breakdown, based on its corresponding expression dictionary of expressing one's feelings, resolves into the amimia image of whole high resolving power, expression mask and anti-sharpening mask;
(4) the amimia image of the resulting whole high resolving power of step (3) is carried out down-sampling, then and the amimia input picture of whole low resolution that obtains of step (2) subtract each other, obtain the low-resolution residual error image; To this low-resolution residual error image, come image is carried out the facial detail information compensation by the image block on the face location of everywhere being carried out neighbour's reconstruct, obtain the high-resolution residual error image;
(5) the amimia image of whole high resolving power that the high-resolution residual error image in the step (4) and step (3) is obtained merges, and obtains final face image super-resolution result, finishes super-resolution rebuilding.
7. the face image super-resolution reconstruction method based on the morphology constituent analysis according to claim 6 is characterized in that, in the described step (2), judges expression in the low resolution input picture by adopting the LBP sorter.
CN2012103897128A 2012-10-15 2012-10-15 Face image super-resolution reconstruction method based on morphological component analysis Pending CN102915527A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012103897128A CN102915527A (en) 2012-10-15 2012-10-15 Face image super-resolution reconstruction method based on morphological component analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103897128A CN102915527A (en) 2012-10-15 2012-10-15 Face image super-resolution reconstruction method based on morphological component analysis

Publications (1)

Publication Number Publication Date
CN102915527A true CN102915527A (en) 2013-02-06

Family

ID=47613882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103897128A Pending CN102915527A (en) 2012-10-15 2012-10-15 Face image super-resolution reconstruction method based on morphological component analysis

Country Status (1)

Country Link
CN (1) CN102915527A (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103325104A (en) * 2013-06-09 2013-09-25 中山大学 Facial image super-resolution reestablishing method based on iterative sparse representation
WO2015192316A1 (en) * 2014-06-17 2015-12-23 Beijing Kuangshi Technology Co., Ltd. Face hallucination using convolutional neural networks
CN106846286A (en) * 2017-01-05 2017-06-13 四川大学 Video super-resolution algorithm for reconstructing based on various complementary priori
CN107154023A (en) * 2017-05-17 2017-09-12 电子科技大学 Face super-resolution reconstruction method based on generation confrontation network and sub-pix convolution
CN107492070A (en) * 2017-07-10 2017-12-19 华北电力大学 A kind of single image super-resolution computational methods of binary channels convolutional neural networks
CN108765297A (en) * 2018-06-14 2018-11-06 厦门大学 Super resolution ratio reconstruction method based on circuit training
CN108830789A (en) * 2018-05-03 2018-11-16 哈尔滨理工大学 Automobile active safety image processing method based on rarefaction representation
CN110288677A (en) * 2019-05-21 2019-09-27 北京大学 It is a kind of based on can deformation structure pedestrian image generation method and device
WO2019184462A1 (en) * 2018-03-30 2019-10-03 Boe Technology Group Co., Ltd. Image processing apparatus, image processing method thereof, image processing system, and training method thereof
CN110310293A (en) * 2019-07-04 2019-10-08 北京字节跳动网络技术有限公司 Human body image dividing method and equipment
WO2020000878A1 (en) * 2018-06-26 2020-01-02 北京字节跳动网络技术有限公司 Method and apparatus for generating image
WO2020000877A1 (en) * 2018-06-26 2020-01-02 北京字节跳动网络技术有限公司 Method and device for generating image
WO2020187042A1 (en) * 2019-03-19 2020-09-24 京东方科技集团股份有限公司 Image processing method, device and apparatus, and computer readable medium
WO2021155675A1 (en) * 2020-02-07 2021-08-12 腾讯科技(深圳)有限公司 Image processing method and apparatus, computer-readable storage medium, and computer device
WO2021258530A1 (en) * 2020-06-22 2021-12-30 北京大学深圳研究生院 Image resolution processing method, device, apparatus, and readable storage medium
US11449751B2 (en) 2018-09-30 2022-09-20 Boe Technology Group Co., Ltd. Training method for generative adversarial network, image processing method, device and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011090789A1 (en) * 2010-01-22 2011-07-28 Thomson Licensing Method and apparatus for video object segmentation
CN102231204A (en) * 2011-06-20 2011-11-02 南京航空航天大学 Sequence image self-adaptive regular super resolution reconstruction method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011090789A1 (en) * 2010-01-22 2011-07-28 Thomson Licensing Method and apparatus for video object segmentation
CN102231204A (en) * 2011-06-20 2011-11-02 南京航空航天大学 Sequence image self-adaptive regular super resolution reconstruction method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YAN LIANG等: "Face hallucination based on morphological component analysis", 《HTTP://WWW.SCIENCEDIRECT.COM/SCIENCE/ARTICLE/PII/S0165168412002812》 *
YAN LIANG等: "Face Hallucination under an Image Decomposition Perspective", 《PATTERN RECOGNITION (ICPR), 2010 20TH INTERNATIONAL CONFERENCE ON》 *
邹耀贤: "基于KPCA的人脸超分辨率算法研究", 《万方学位论文数据库》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103325104A (en) * 2013-06-09 2013-09-25 中山大学 Facial image super-resolution reestablishing method based on iterative sparse representation
CN103325104B (en) * 2013-06-09 2016-05-11 中山大学 Based on the face image super-resolution reconstruction method of iteration sparse expression
WO2015192316A1 (en) * 2014-06-17 2015-12-23 Beijing Kuangshi Technology Co., Ltd. Face hallucination using convolutional neural networks
CN105960657A (en) * 2014-06-17 2016-09-21 北京旷视科技有限公司 Face hallucination using convolutional neural networks
CN106846286A (en) * 2017-01-05 2017-06-13 四川大学 Video super-resolution algorithm for reconstructing based on various complementary priori
CN106846286B (en) * 2017-01-05 2019-10-29 四川大学 Video super-resolution algorithm for reconstructing based on a variety of complementary priori
CN107154023A (en) * 2017-05-17 2017-09-12 电子科技大学 Face super-resolution reconstruction method based on generation confrontation network and sub-pix convolution
CN107154023B (en) * 2017-05-17 2019-11-05 电子科技大学 Based on the face super-resolution reconstruction method for generating confrontation network and sub-pix convolution
CN107492070A (en) * 2017-07-10 2017-12-19 华北电力大学 A kind of single image super-resolution computational methods of binary channels convolutional neural networks
CN107492070B (en) * 2017-07-10 2019-12-03 华北电力大学 A kind of single image super-resolution calculation method of binary channels convolutional neural networks
CN110322400A (en) * 2018-03-30 2019-10-11 京东方科技集团股份有限公司 Image processing method and device, image processing system and its training method
CN110322400B (en) * 2018-03-30 2021-04-27 京东方科技集团股份有限公司 Image processing method and device, image processing system and training method thereof
WO2019184462A1 (en) * 2018-03-30 2019-10-03 Boe Technology Group Co., Ltd. Image processing apparatus, image processing method thereof, image processing system, and training method thereof
US11189013B2 (en) 2018-03-30 2021-11-30 Boe Technology Group Co., Ltd. Image processing apparatus, image processing method thereof, image processing system, and training method thereof
CN108830789A (en) * 2018-05-03 2018-11-16 哈尔滨理工大学 Automobile active safety image processing method based on rarefaction representation
CN108765297A (en) * 2018-06-14 2018-11-06 厦门大学 Super resolution ratio reconstruction method based on circuit training
CN108765297B (en) * 2018-06-14 2020-07-17 厦门大学 Super-resolution reconstruction method based on cyclic training
WO2020000878A1 (en) * 2018-06-26 2020-01-02 北京字节跳动网络技术有限公司 Method and apparatus for generating image
WO2020000877A1 (en) * 2018-06-26 2020-01-02 北京字节跳动网络技术有限公司 Method and device for generating image
US11449751B2 (en) 2018-09-30 2022-09-20 Boe Technology Group Co., Ltd. Training method for generative adversarial network, image processing method, device and storage medium
WO2020187042A1 (en) * 2019-03-19 2020-09-24 京东方科技集团股份有限公司 Image processing method, device and apparatus, and computer readable medium
CN110288677A (en) * 2019-05-21 2019-09-27 北京大学 It is a kind of based on can deformation structure pedestrian image generation method and device
CN110310293A (en) * 2019-07-04 2019-10-08 北京字节跳动网络技术有限公司 Human body image dividing method and equipment
WO2021155675A1 (en) * 2020-02-07 2021-08-12 腾讯科技(深圳)有限公司 Image processing method and apparatus, computer-readable storage medium, and computer device
WO2021258530A1 (en) * 2020-06-22 2021-12-30 北京大学深圳研究生院 Image resolution processing method, device, apparatus, and readable storage medium

Similar Documents

Publication Publication Date Title
CN102915527A (en) Face image super-resolution reconstruction method based on morphological component analysis
CN112634137B (en) Hyperspectral and panchromatic image fusion method for extracting multiscale spatial spectrum features based on AE
Lim et al. DSLR: Deep stacked Laplacian restorer for low-light image enhancement
CN110119780B (en) Hyper-spectral image super-resolution reconstruction method based on generation countermeasure network
CN110136062B (en) Super-resolution reconstruction method combining semantic segmentation
CN108022213A (en) Video super-resolution algorithm for reconstructing based on generation confrontation network
CN105046672B (en) A kind of image super-resolution rebuilding method
CN110415199B (en) Multispectral remote sensing image fusion method and device based on residual learning
Huang et al. Deep hyperspectral image fusion network with iterative spatio-spectral regularization
CN112541864A (en) Image restoration method based on multi-scale generation type confrontation network model
Chen et al. Convolutional neural network based dem super resolution
CN110084773A (en) A kind of image interfusion method based on depth convolution autoencoder network
CN105488759B (en) A kind of image super-resolution rebuilding method based on local regression model
CN105631807A (en) Single-frame image super resolution reconstruction method based on sparse domain selection
CN111861886B (en) Image super-resolution reconstruction method based on multi-scale feedback network
Guan et al. Srdgan: learning the noise prior for super resolution with dual generative adversarial networks
CN116664397B (en) TransSR-Net structured image super-resolution reconstruction method
CN114170088A (en) Relational reinforcement learning system and method based on graph structure data
CN115578255A (en) Super-resolution reconstruction method based on inter-frame sub-pixel block matching
CN109559278B (en) Super resolution image reconstruction method and system based on multiple features study
CN112686830A (en) Super-resolution method of single depth map based on image decomposition
CN109886869A (en) A kind of unreal structure method of face of the non-linear expansion based on contextual information
CN113240581A (en) Real world image super-resolution method for unknown fuzzy kernel
Ren et al. A lightweight object detection network in low-light conditions based on depthwise separable pyramid network and attention mechanism on embedded platforms
Han et al. MPDNet: An underwater image deblurring framework with stepwise feature refinement module

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C05 Deemed withdrawal (patent law before 1993)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130206