# Patent

## Title
Method and apparatus for judging direction of blur and computer-readable recording medium storing a program therefor

## Abstract

Without use of special hardware such as an angular rate sensor, a direction of blur can be judged with high accuracy. A parameter acquisition unit obtains weighting parameters for principal components representing directions of blur in a predetermined structure in an image by fitting to the structure a mathematical model generated by a method of Active Appearance Models (AAM) using a plurality of sample images representing the structure in different conditions of blur. A blur direction judgment unit judges the direction of blur based on the weighting parameters, and a blur width acquisition unit finds a width of blur based on an edge component found by an edge detection unit in the direction perpendicular to the direction of blur. A blur correction unit corrects the blur in the image based on the direction and the width of blur.

## Background
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to a method of judging a direction of blur represented in an input image and an apparatus therefor. The present invention also relates to a computer-readable recording medium storing a program that causes a computer to execute the method.
2. Description of the Related Art
As a result of reduction in size and weight of imaging devices such as digital cameras and camera phones, camera shake tends to occur due to imaging device movement at the time of shuttering. In addition, progress is being made in increase in the number of pixels and in installation of a zooming function as standard equipment. Therefore, camera shake tends to become conspicuous in photographed images.
As a method of detecting such camera shake has been known a method using an angular rate sensor installed in an imaging device (see U.S. Patent Application Publication No. 20030002746).
However, the method described in U.S. Patent Application Publication No. 20030002746 needs special hardware such as an angular rate sensor, and installation of such hardware in a small imaging device such as a camera phone is difficult.
Furthermore, no method has been disclosed for detecting camera shake based only on a photographed image in the case where no information is available on camera shake detected by such hardware.
SUMMARY OF THE INVENTION
The present invention has been conceived based on consideration of the above circumstances. An object of the present invention is therefore to provide a method enabling judgment of a direction of blur such as camera shake with high accuracy and without special hardware such as an angular rate sensor, an apparatus for carrying out the judgment, and a computer-readable recording medium storing a program therefor.
A blur direction judgment method of the present invention comprises the steps of:
obtaining one or more weighting parameters for one or more statistical characteristic quantities representing a direction of blur in a predetermined structure in an input image by fitting a model representing the structure to the structure in the input image, the model having been obtained by carrying out predetermined statistical processing on a plurality of images representing the predetermined structure in different conditions of blur, and the model representing the structure by one or more statistical characteristic quantities including the statistical characteristic quantity or quantities representing the direction of blur and by one or more weighting parameters for weighting the statistical characteristic quantity or quantities according to an individual characteristic of the structure; and
judging the direction of blur represented in the input image according to a value or values of the weighting parameter or parameters having been obtained.
A blur direction judgment apparatus of the present invention is an apparatus for carrying out the blur direction judgment method described above. More specifically, the blur direction judgment apparatus of the present invention comprises:
parameter acquisition means for obtaining one or more weighting parameters for one or more statistical characteristic quantities representing a direction of blur in a predetermined structure in an input image by fitting a model representing the structure to the structure in the input image, the model having been obtained by carrying out predetermined statistical processing on a plurality of images representing the predetermined structure in different conditions of blur, and the model representing the structure by one or more statistical characteristic quantities including the statistical characteristic quantity or quantities representing the direction of blur and by one or more weighting parameters for weighting the statistical characteristic quantity or quantities according to an individual characteristic of the structure; and
blur direction judgment means for judging the direction of blur represented in the input image according to a value or values of the weighting parameter or parameters having been obtained by the parameter acquisition means.
A computer-readable recording medium of the present invention stores a program that causes a computer to execute the blur direction judgment method described above (that is, a program that causes a computer to function as the means described above).
Hereinafter, the blur direction judgment apparatus, the blur direction judgment method, and the blur direction judgment program of the present invention are described in detail.
As a method of generating the model representing the predetermined structure in the present invention, a method of AAM (Active Appearance Model) can be used. An AAM is one of approaches in interpretation of the content of an image by using a model. For example, in the case where a human face is a target of interpretation, a mathematical model of human face is generated by carrying out principal component analysis on face shapes in a plurality of images to be learned and on information of luminance after normalization of the shapes. A face in a new input image is then represented by principal components in the mathematical model and corresponding weighting parameters, for face image reconstruction (T. F. Cootes et al., “Active Appearance Models”, Proc. European Conference on Computer Vision, vol. 2, pp. 484-498, Springer, 1998, Germany; hereinafter referred to as Document 1).
The blur may be camera shake or movement of a subject.
The conditions of blur represent at least the direction of blur. In addition, the conditions may include a condition representing the magnitude of blur.
It is preferable for the predetermined structure to be suitable for modeling. In other words, variations in shape and luminance of the predetermined structure in images thereof preferably fall within a predetermined range. Especially, it is preferable for the predetermined structure to generate the statistical characteristic quantity or quantities contributing more to the shape and luminance thereof through the statistical processing thereon. Furthermore, it is preferable for the predetermined structure to have a high probability of being a main part of image. More specifically, the predetermined structure can be a human face.
The plurality of images representing the predetermined structure in different conditions of blur may be images obtained by actually photographing the predetermined structure in the different conditions of blur. Alternatively, the images may be generated through simulation based on an image of the structure photographed in a specific condition of blur.
It is preferable for the predetermined statistical processing to be dimension reduction processing that can represent the predetermined structure by the statistical characteristic quantity or quantities of fewer dimensions than the number of pixels representing the predetermined structure. More specifically, the predetermined statistical processing may be multivariate analysis such as principal component analysis. In the case where principal component analysis is carried out as the predetermined statistical processing, the statistical characteristic quantity or quantities refers/refer to principal components obtained through the principal component analysis.
In the case where the predetermined statistical processing is principal component analysis, principal components of higher orders contribute more to the shape and luminance than principal components of lower orders.
In the statistical characteristic quantity or quantities, at least information based on luminance in the structure needs to be represented, since the blur is represented by distribution of luminance in the image.
The statistical characteristic quantity or quantities representing the direction of blur may be a single statistical characteristic quantity or a plurality of statistical characteristic quantities. For example, the statistical characteristic quantity or quantities may be statistical characteristic quantities that vary in respective directions of blur.
The (predetermined) structure in the input image may be detected automatically or manually. In addition, the present invention may further comprise the step (or means) for detecting the structure in the input image. Alternatively, the structure may have been detected in the input image in the present invention.
A plurality of models in the present invention may be prepared for respective properties of the predetermined structure. In this case, the steps (or means) may be added to the present invention for obtaining at least one property of the structure in the input image and for selecting one of the models according to the at least one obtained property. The weighting parameter or parameters can be obtained by fitting the selected model to the structure in the input image.
The properties refer to gender, age, and race in the case where the predetermined structure is human face. The property may be information for identifying an individual. In this case, the models for the respective properties refer to models for respective individuals.
As a specific method of obtaining the property may be listed image recognition processing having been known (such as image recognition processing described in Japanese Unexamined Patent Publication No. 11 (1999) -175724). Alternatively, the property may be inferred or obtained based on information such as GPS information accompanying the input image.
Fitting the model representing the structure to the structure in the input image refers to calculation or the like for representing the structure in the input image by the model. More specifically, in the case where the method of AAM described above is used, fitting the model refers to finding values of weighting parameters for the respective principal components in the mathematical model.
As a specific method of judging the direction of blur represented in the input image according to the weighting parameter or parameters having been obtained for the statistical characteristic quantity or quantities representing the direction of blur, a relationship is experimentally and statistically found in advance between the direction of blur and the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur. Based on the relationship is then found the direction of blur corresponding to the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities obtained by fitting the model to the structure in the input image.
A width of blur may also be found in the direction of blur having been judged. More specifically, the width may be found based on the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur having been judged. In this case, a relationship is experimentally and statistically found in advance between the width of blur in the direction thereof and the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur, and the width of blur is found based on the relationship. Alternatively, a width of an edge perpendicular to the direction of blur having been judged may be found and used as the width of blur. In this case, edges perpendicular to the direction of blur are found through direction-dependent edge detection processing having been known. A mean width of the detected edges is found as the width of blur.
Furthermore, the blur represented in the input image may be corrected based on the direction of blur having been judged and the width of blur having been found. At this time, the blur in the structure in the input image may be corrected by changing the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur. In addition, sharpness correction may be carried out on the entire input image for enhancing the edges in the direction substantially perpendicular to the direction of blur having been judged, at processing strength in accordance with the width of blur having been found (see Japanese Unexamined Patent Publication No. 2004-070421). Furthermore, blur correction may be carried out on the structure in the input image by changing the value or value of the weighting parameter or parameters and on the remaining part of the input image according to the sharpness correction described above.
According to the blur direction judgment method, the blur direction judgment apparatus, and the computer-readable recording medium storing the blur direction judgment program of the present invention, the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur in the structure in the input image is/are obtained by fitting to the structure in the input image the model representing the predetermined structure by use of the statistical characteristic quantity or quantities including the statistical characteristic quantity or quantities representing the direction of blur and the weighting parameter or parameters therefor. Based on the value or values of the weighting parameter or parameters having been obtained, the direction of blur represented in the input image is judged. Therefore, without use of special hardware such as an angular rate sensor, blur such as camera shake can be detected with high accuracy, and the direction of blur can be judged in the present invention in a photographed image even in the case where no information on blur is detected at the time of photography.
In the case where the step (or the means) for detecting the structure in the input image is added, automatic detection of the structure can be carried out, leading to improvement in operability.
In the case where the plurality of models are prepared for the respective properties of the predetermined structure in the present invention while the steps (or the means) are added for obtaining the property of the structure in the input image and for selecting one of the models in accordance with the at least one obtained property, if the weighting parameter or parameters is/are obtained by fitting the selected model to the structure in the input image, the structure in the input image can be fit by the model that is more suitable. Therefore, processing accuracy is improved.
In the case where the steps (means) are added for finding the width of blur in the direction of blur having been judged and for correcting the blur in the input image based on the direction of blur and the width of blur having been found, the blur can be corrected based on a result of high-accuracy blur detection, which leads to generation of a preferable image.
In the case where the width of blur is found based on the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur, the width of blur can be found without image analysis processing such as edge detection, which enables simpler configuration and improvement in processing efficiency.
Likewise, in the case where the blur in the structure in the input image is corrected by changing the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur, the blur can be corrected without other image processing such as sharpness correction, which enables simpler configuration and improvement in processing efficiency.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows hardware configuration of a digital photograph printer in an embodiment of the present invention;
FIG. 2 is a block diagram showing functions and a flow of processing in the digital photograph printer in the embodiment and in a digital camera in another embodiment of the present invention;
FIGS. 3A and 3B show examples of screens displayed on a display of the digital photograph printer and the digital camera in the embodiments;
FIG. 4 is a block diagram showing details of blur direction judging process and blur correcting process in one aspect of the present invention;
FIG. 5 is a flow chart showing a procedure for generating a mathematical model of face image in the present invention;
FIG. 6 shows an example of how feature points are set in a face;
FIG. 7 shows how a face shape changes with change in values of weight coefficients for eigenvectors of principal components obtained through principal component analysis on the face shape;
FIG. 8 shows luminance in mean face shapes converted from face shapes in sample images;
FIG. 9 shows how luminance in a face change with change in values of weight coefficients for eigenvectors of principal components obtained by principal component analysis on the luminance in the face;
FIG. 10 is a block diagram showing details of blur direction judging process and blur correcting process in another aspect of the present invention;
FIG. 11 shows a structure of a reference table used in the blur correcting process and an example of values therein;
FIG. 12 is a block diagram showing an advanced aspect of the blur correcting process in the present invention; and
FIG. 13 shows the configuration of the digital camera in the embodiment of the present invention.



DESCRIPTION OF THE PREFERRED EMBODIMENTS
Hereinafter, embodiments of the present invention are described with reference to the accompanying drawings.
FIG. 1 shows the hardware configuration of a digital photograph printer in an embodiment of the present invention. As shown in FIG. 1, the digital photograph printer comprises a film scanner 51, a flat head scanner 52, a media drive 53, a network adapter 54, a display 55, a keyboard 56, a mouse 57, a hard disc 58, and a photographic print output machine 59, all of which are connected to an arithmetic and control unit 50.
In cooperation with a CPU, a main storage, and various input/output interfaces, the arithmetic and control unit 50 controls a processing flow regarding an image, such as input, correction, manipulation, and output thereof, by executing a program installed from a recording medium such as a CD-ROM. In addition, the arithmetic and control unit 50 carries out image processing calculation for image correction and manipulation. A blur direction judging process and a blur correcting process of the present invention is also carried out by the arithmetic and control unit 50.
The film scanner 51 photoelectrically reads an APS negative film or a 135-mm negative film developed by a film developer (not shown) for obtaining digital image data P0 representing a photograph image recorded on the negative film.
The flat head scanner 52 photoelectrically reads a photograph image represented in the form of hard copy such as an L-size print, for obtaining digital image data P0.
The media drive 53 obtains digital image data P0 representing a photograph image recorded in a recording medium such as a memory card, a CD, or a DVD. The media drive 53 can also write image data P2 to be output therein. The memory card stores image data representing an image photographed by a digital camera, for example. The CD or the DVD stores data of an image read by the film scanner regarding a previous print order, for example.
The network adapter 54 obtains image data P0 from an order reception machine (not shown) in a network photograph service system having been known. The image data P0 are image data used for a photograph print order placed by a user, and sent from a personal computer of the user via the Internet or via a photograph order reception machine installed in a photo laboratory.
The display 55 displays an operation screen for input, correction, manipulation, and output of an image carried out by the digital photograph printer. A menu for selecting the content of operation and an image to be processed are also displayed thereon, for example. The keyboard 56 and the mouse 57 are used for inputting a processing instruction.
The hard disc 58 stores a program for controlling the digital photograph printer. In the hard disc 58 are also stored temporarily the image data P0 obtained by the film scanner 51, the flat head scanner 52, the media drive 53, and the network adapter 54, in addition to image data P1 having been subjected to image correction (hereinafter referred to as the corrected image data P1) and the image data P2 having been subjected to image manipulation (the image data to be output).
The photograph print output machine 59 carries out laser scanning exposure of photographic printing paper, image development thereon, and drying thereof, based on the image data P2 representing the image to be output. The photograph print output machine 59 also prints printing information on the backside of the paper, cuts the paper for each print, and sorts the paper for each order. The manner of printing may be a laser exposure thermal development dye transfer process or the like.
FIG. 2 is a block diagram showing functions of the digital photograph printer and the flow of processing carried out therein. As shown in FIG. 2, the digital photograph printer comprises image input means 1, image correction means 2, image manipulation means 3, and image output means 4 in terms of the functions. The image input means 1 inputs the image data P0 of an image to be printed. The image correction means 2 uses the image data P0 as input, and carries out automatic image quality correction of the image represented by the image data P0 (hereinafter, image data and an image represented by the image data are represented by the same reference code) through image processing according to a predetermined image processing condition. The image manipulation means 3 uses the corrected image data P1 having been subjected to the automatic correction as input, and carries out image processing according to an instruction from an operator. The image output means 4 uses the processed image data P2 as input, and outputs a photographic print or outputs the processed image data P2 in a recording medium.
The image correction means 2 carries out processing such as white balance adjustment, contrast correction, sharpness correction, and noise reduction and removal. The image manipulation means 3 carries out manual correction on a result of the processing carried out by the image correction means 2. In addition, the image manipulation means 3 carries out image manipulation such as trimming, scaling, conversion to sepia image, conversion to monochrome image, and compositing with an ornamental frame, in addition to the blur direction judging process and the blur correcting process of the present invention.
Operation of the digital photograph printer and the flow of the processing therein are described next.
The image input means 1 firstly inputs the image data P0. In the case where an image recorded on a developed film is printed, the operator sets the film on the film scanner 51. In the case where image data stored in a recording medium such as a memory card are printed, the operator sets the recording medium in the media drive 53. A screen for selecting a source of input of the image data is displayed on the display 55, and the operator carries out the selection by using the keyboard 56 or the mouse 57. In the case where “film” has been selected as the source of input, the film scanner 51 photoelectrically reads the film set thereon, and carries out digital conversion. The image data P0 generated in this manner are then sent to the arithmetic and control unit 50. In the case where “hard copy” such as a photographic print has been selected, the flat head scanner 52 photoelectrically reads the hard copy set thereon, and carries out digital conversion. The image data P0 generated in this manner are then sent to the arithmetic and control unit 50. In the case where “recording medium” such as a memory card has been selected, the arithmetic and control unit 50 reads the image data P0 stored in the recording medium such as a memory card set in the media drive 53. In the case where an order has been placed in a network photograph service system or by a photograph order reception machine in a store, the arithmetic and control unit 50 receives the image data P0 via the network adapter 54. The image data P0 obtained in this manner are temporarily stored in the hard disc 58.
The image correction means 2 then carries out the automatic image quality correction on the image represented by the image data P0. More specifically, publicly known processing such as white balance adjustment, contrast correction, sharpness correction, and noise reduction and removal is carried out based on a setup condition set on the printer in advance, according to an image processing program executed by the arithmetic and control unit 50. The corrected image data P1 are output to be stored in a memory of the arithmetic and control unit 50. Alternatively, the corrected image data P1 may be stored temporarily in the hard disc 58.
The image manipulation means 3 thereafter generates a thumbnail image of the corrected image P1, and causes the display 55 to display the thumbnail image. FIG. 3A shows an example of a screen displayed on the display 55. The operator confirms displayed thumbnail images, and selects any one of the thumbnail images that needs manual image-quality correction or order processing for image manipulation while using the keyboard 56 or the mouse 57. In FIG. 3A, the image in the upper left corner (DSCF0001) is selected. As shown in FIG. 3B as an example, the selected thumbnail image is enlarged and displayed on the display 55, and buttons are displayed for selecting the content of manual correction and manipulation on the image. The operator selects a desired one of the buttons by using the keyboard 56 or the mouse 57, and carries out detailed setting of the selected content if necessary. The image manipulation means 3 carries out the image processing according to the selected content, and outputs the processed image data P2. The image data P2 are stored in the memory of the arithmetic and control unit 50 or stored temporarily in the hard disc 58. The program executed by the arithmetic and control unit 50 controls image display on the display 55, reception of input from the keyboard 56 or the mouse 57, and image processing such as manual correction and manipulation carried out by the image manipulation means 3.
The image output means 4 finally outputs the image P2. The arithmetic and control unit 50 causes the display 55 to display a screen for image destination selection, and the operator selects a desired one of destinations by using the keyboard 56 or the mouse 57. The arithmetic and control unit 50 sends the image data P2 to the selected destination. In the case where a photographic print is generated, the image data P2 are sent to the photographic print output machine 59 by which the image data P2 are output as a photographic print. In the case where the image data P2 are recorded in a recording medium such as a CD, the image data P2 are written in the CD or the like set in the media drive 53.
Below is described the blur direction judging process of the present invention followed by the blur correcting process carried out in the case where “Blur Correction” is selected in the screen shown in FIG. 3B. FIG. 4 is a block diagram showing details of the blur direction judging process and the blur correcting process. As shown in FIG. 4, the blur direction judging process and the blur correcting process is carried out by a face detection unit 31, a parameter acquisition unit 32, a blur direction judgment unit 33, an edge detection unit 34, a blur width acquisition unit 35, and a blur correction unit 36. The face detection unit 31 detects a face region P1f in the image P1. The parameter acquisition unit 32 fits to the detected face region P1f a mathematical model M generated by a method of AAM (see aforementioned Document 1) based on a plurality of sample images representing human faces in different conditions of blur, and obtains weight coefficients λd1, λd2, λd3, . . . , λdJ, for principal components representing directions of blur in the face region P1f. The blur direction judgment unit 33 judges a direction D of blur based on the weight coefficients having been obtained. The edge detection unit 34 detects edge components E perpendicular to the direction D in the image P1. The blur width acquisition unit 35 finds a width W of blur based on the edge component E. The blur correction unit 36 corrects the blur in the input image P1 based on the direction D and the width W of blur. The processing described above is controlled by the program installed in the arithmetic and control unit 50.
The mathematical model M is generated according to a flow chart shown in FIG. 5, and installed in advance together with the programs described above. Hereinafter, how the mathematical model M is generated is described.
For each of the sample images representing human faces in different directions and magnitudes of blur, feature points are set therein as shown in FIG. 6 for representing face shape (Step #1). In this case, the number of the feature points is 122. However, only 60 points are shown in FIG. 6 for simplification. Which part of face is represented by which of the feature points is predetermined, such as the left corner of the left eye represented by the first feature point and the center between the eyebrows represented by the 38th feature point. Each of the feature points may be set manually or automatically according to recognition processing. Alternatively, the feature points may be set automatically and later corrected manually upon necessity.
Based on the feature points set in each of the sample images, mean face shape is calculated (Step #2). More specifically, mean values of coordinates of the feature points representing the same part are found among the sample images.
Principal component analysis is then carried out based on the coordinates of the mean face shape and the feature points representing the face shape in each of the sample images (Step #3). As a result, any face shape can be approximated by Equation (1) below:




  
    
      
        S
        =
        
          
            S
            0
          
          +
          
            
              ∑
              
                i
                =
                1
              
              n
            
            ⁢
            
              
                p
                i
              
              ⁢
              
                b
                i
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




S and S0 are shape vectors represented respectively by simply listing the coordinates of the feature points (x1, Y1, . . . , X122, Y122) in the face shape and in the mean face shape, while pi and bi are an eigenvector representing the ith principal component for the face shape obtained by the principal component analysis and a weight coefficient therefor, respectively. FIG. 7 shows how face shape changes with change in values of the weight coefficients b1 and b2 for the eigenvectors pi and P2 as the highest and second-highest order principal components obtained by the principal component analysis. The change ranges from −3sd to +3sd where sd refers to standard deviation of each of the weight coefficients b1 and b2 in the case where the face shape in each of the sample images is represented by Equation (1). The face shape in the middle of 3 faces for each of the components represents the face shape in the case where the value of the corresponding weight coefficient is the mean value. In this example, a component contributing to face outline has been extracted as the first principal component through the principal component analysis. By changing the weight coefficient b1, the face shape changes from an elongated shape (corresponding to −3sd) to a round shape (corresponding to +3sd). Likewise, a component contributing to how much the mouth is open and to length of chin has been extracted as the second principal component. By changing the weight coefficient b2, the face changes from a state of open mouth and long chin (corresponding to −3sd) to a state of closed mouth and short chin (corresponding to +3sd). The smaller the value of i, the better the component explains the shape. In other words, the ith component contributes more to the face shape as the value of i becomes smaller.
Each of the sample images is then subjected to conversion (warping) into the mean face shape obtained at Step #2 (Step #4). More specifically, shift values are found between each of the sample images and the mean face shape, for the respective feature points. In order to warp pixels in each of the sample images to the mean face shape, shift values to the mean face shape are calculated for the respective pixels in each of the sample images according to 2-dimensional 5-degree polynomials (2) to (5) using the shift values having been found:




  
    
      
        
          x
          ′
        
        =
        
          x
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            x
          
        
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          y
          ′
        
        =
        
          y
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            y
          
        
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          x
        
        =
        
          
            ∑
            
              i
              =
              0
            
            n
          
          ⁢
          
            
              ∑
              
                j
                =
                0
              
              
                n
                -
                i
              
            
            ⁢
            
              
                a
                ij
              
              ·
              
                x
                i
              
              ·
              
                y
                j
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          x
        
        =
        
          
            ∑
            
              i
              =
              0
            
            n
          
          ⁢
          
            
              ∑
              
                j
                =
                0
              
              
                n
                -
                i
              
            
            ⁢
            
              
                b
                ij
              
              ·
              
                x
                i
              
              ·
              
                y
                j
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




In Equations (2) to (5) above, x and y denote the coordinates of each of the feature points in each of the sample images while x′ and y′ are coordinates in the mean face shape to which x and y are warped. The shift values to the mean shape are represented by Δx and Δy with n being the number of dimensions while aij and bij are coefficients. The coefficients for polynomial approximation can be found by using a least square method. In the case that the coordinates of a feature point become non-integer values (that is, values including decimals), pixel values of four pixels (having integer coordinates) that surround the coordinates after warping are found through linear approximation of the pixel values. More specifically, for 4 pixels surrounding the coordinates of the non-integer values generated by warping, the pixel values for each of the 4 pixels are determined in proportion to a distance thereto from the coordinates generated by warping. FIG. 8 shows how the face shape of each of 3 sample images is changed to the mean face shape.
Thereafter, principal component analysis is carried out, using as variables the luminance of each of the pixels in each of the sample images after the change to the mean face shape (Step #5). As a result, the luminance in the mean face shape converted from any arbitrary face image can be approximated by Equation (6) below:




  
    
      
        A
        =
        
          
            A
            0
          
          +
          
            
              ∑
              
                i
                =
                1
              
              m
            
            ⁢
            
              
                q
                i
              
              ⁢
              
                λ
                i
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




In Equation (6), A denotes a luminance vector (a1, . . . , am) represented by listing the luminance at each of the pixels in the mean face shape (where a represents the luminance while 1 to m refer to subscripts for identifying the respective pixels with m being the total number of pixels in the mean face shape). A0 is a mean face luminance vector represented by listing mean values of the luminance at each of the pixels in the mean face shape while qi and λi refer to an eigenvector representing the ith principal component for the luminance in the face obtained by the principal component analysis and a weight coefficient therefor, respectively. The smaller the value of i is, the better the component explains the luminance. In other words, the component contributes more to the luminance as the value of i becomes smaller.
FIG. 9 shows how the luminance values of faces change with change in values of the weight coefficients λi1 and λi2 for the eigenvectors qi1 and qi2 representing the i1th and i2th principal components obtained through the principal component analysis. The change in the weight coefficients ranges from −3sd to +3sd where sd refers to standard deviation of each of the values of the weight coefficients λi1, and λi2 in the case where the luminance in each of the face sample images are represented by Equation (6) above. For each of the principal components, the face in the middle of the 3 images corresponds to the case where the corresponding weight coefficient λi1, or λi2 is the mean value. In the examples shown in FIG. 9, a component contributing to presence or absence of beard has been extracted as the i1th principal component through the principal component analysis. By changing the weight coefficient λi1, the face changes from the face with dense beard (corresponding to −3sd) to the face with no beard (corresponding to +3sd). Since the face images having different directions and magnitudes of blur are used as the sample images in this embodiment, a component contributing to blur has been extracted as the i2th principal component through the principal component analysis. By changing the weight coefficient λi2, the magnitude of blur in the right-left direction changes. Likewise, components contributing to blur in other directions are extracted as other principal components. How each of the principal components contributes to what factor is determined through interpretation. The directions of blur can be represented by a plurality of principal components as in the above example or by a single principal component.
Through the processing from Step #1 to #5 described above, the mathematical model M is generated. In other words, the mathematical model M is represented by the eigenvectors pi representing the face shape and the eigenvectors qi representing the face luminance in the mean face shape, and the number of the eigenvectors is far smaller for pi and for qi than the number of pixels forming the face image. In other words, the mathematical model M has been compressed in terms of dimension. In the example described in Document 1, 122 feature points are set for a face image of approximately 10,000 pixels, and a mathematical model of face image represented by 23 eigenvectors for face shape and 114 eigenvectors for face luminance has been generated through the processing described above. By changing the weight coefficients for the respective eigenvectors, 90% of variations or more in face shape and luminance can be expressed.
A flow of the blur direction judging process and the blur correcting process based on the method of AAM using the mathematical model M is described below, with reference to FIG. 4.
The face detection unit 31 reads the image data P1, and detects the face region P1f in the image P1. More specifically, a first characteristic quantity representing a direction of a gradient vector showing a direction and a magnitude of an edge at each of the pixels in the image P1 is firstly input to a plurality of first detectors (which will be described later) for judgment as to whether a face candidate region exists in the image P1, as has been described in Japanese Unexamined Patent Publication No. 2005-108195 (hereinafter referred to as Reference 2). In the case where a face candidate region exists, the region is extracted, and the magnitude of the gradient vector is normalized at each of the pixels in the face candidate region. A second characteristic quantity representing the direction and the magnitude of the normalized gradient vector is then input to second detectors (which will be described later) for judgment as to whether the extracted face candidate region is a true face region. In the case where the region is a true face region, the region is detected as the face region P1f. The first/second detectors have been obtained through training using a method of machine learning such as AdaBoost that uses as input the first/second characteristic quantity calculated for face sample images and non-face sample images.
As a method of detection of the face region P1f may be used various known methods such as a method using a correlation score between an eigen-face representation and an image as has been described in Published Japanese Translation of a PCT Application No. 2004-527863 (hereinafter referred to as Reference 3). Alternatively, the face region can be detected by using a knowledge base, characteristics extraction, skin-color detection, template matching, graph matching, and a statistical method (such as a method using a neural network, SVM, and HMM), for example. Furthermore, the face region P1f may be specified manually by use of the keyboard 56 and the mouse 57 when the image P1 is displayed on the display 55. Alternatively, a result of automatic detection of the face region may be corrected manually.
The parameter acquisition unit 32 fits the mathematical model M to the face region P1f. More specifically, the parameter acquisition unit 32 reconstructs the image according to Equations (1) and (6) described above while sequentially changing the values of the weight coefficients bi and λi for the eigenvectors pi and qi corresponding to the principal components in order of higher order in Equations (1) and (6). The parameter acquisition unit 32 then finds the values of the weight coefficients bi and λi causing a difference between the reconstructed image and the face region P1f to become minimal (see Reference 3 for details). Among the weight coefficients having been found, the weight coefficients corresponding to the principal components representing the directions of blur are the weight coefficients λd1, λd2, λd3, . . . , λdJ. It is preferable for the values of the weight coefficients bi and λi to range only from −3sd to +3sd, for example, where sd refers to the standard deviation in each of distributions of bi and λi when the sample images used at the time of generation of the model are represented by Equations (1) and (6). In the case where the values do not fall within the range, it is preferable for the weight coefficients to take the mean values in the distributions. In this manner, erroneous application of the model can be avoided.
The blur direction judgment unit 33 judges the direction D of blur based on the weight coefficients λd1, λd2, λd3, . . . , λdJ having been obtained. More specifically, a direction vector representing the direction of blur is found in the face region P1f by adding unit direction vectors representing the respective directions of blur weighted by the absolute values of the weight coefficients corresponding thereto, and the direction represented by the direction vector is used as the direction D of the blur. Alternatively, a reference table may be found experimentally and statistically in advance for defining a relationship between a combination of the values of the weight coefficients λd1, λd2, λd3, . . . , λdJ and the direction D of the blur. The direction D of blur may be found with reference to the reference table according to the values of the weight coefficients λd1, λd2, λd3, . . . , λdJ having been obtained.
The edge detection unit 34 detects the edge components E in the direction perpendicular to the direction D through direction-dependent edge detection processing according to an operator Ei (where i=1, 2, . . . denotes an index for identifying each of the directions of edge) that detects an edge perpendicular to the direction D by using spatial first derivative or the like.
The blur width acquisition unit 35 finds a mean value of widths of the edge components E as the width W.
The blur correction unit 36 refers to a reference table Tw based on the width W of blur, and obtains an enhancement factor αj (where j (=1, 2, . . . ) refers to an index for identifying the enhancement factor) representing a degree of edge enhancement according to the magnitude of the width W. The reference table Tw, defines a relationship between the width W of blur and the enhancement factor αj, based on experimental and statistical data.
The blur correction unit 36 further carries out selective sharpening for enhancing the edges in the direction perpendicular to the direction of blur in the image data P1, based on the enhancement factor αj and the operator Ei once used by the edge detection unit 34. The blur correction unit 36 then outputs the image data P2. More specifically, the blur correction unit 36 carries out the conversion processing represented by the following equation:

P2(x,y)=P1(x,y)+αj·Ei·L(x,y)  (7)

where P2 (x,y) and P1 (x,y) respectively denote the output image and the input image while L (x,y) refers to a Laplacian image of the input image P1.

As has been described above, according to the blur direction judging process and the blur correcting process in the embodiment of the present invention, the parameter acquisition unit 32 obtains the weight coefficients λd1, λd2, λd3, . . . , λdJ corresponding to the principal components representing the directions of blur in the face region P1f detected by the face detection unit 31 in the image P1 by fitting to the face region P1f the mathematical model M generated according to the method of AAM using the sample images representing human faces in different conditions of blur. The blur direction judgment unit 33 then judges the direction D of blur based on the weight coefficients λd1, λd2, λd3, . . . , λdJ and the blur width acquisition unit 35 finds the width W of blur based on the edge components E detected by the edge detection unit 34 in the direction perpendicular to the direction D. The blur correction unit 36 then corrects the blur in the input image P1 based on the direction D and the width W of blur. Therefore, without special hardware such as an angular rate sensor, the detection and correction of blur can be carried out with high accuracy, and the image can be obtained in a more preferable state.
In the embodiment described above, the edge detection unit 34 may find the operator Ei based on the weight coefficients λd1, λd2, λd3, . . . , λdJ by referring to a reference table found experimentally and statistically in advance for defining a relationship between the operator Ei and the combination of the values of the weight coefficients λd1, λd2, λd3, . . . , λdJ without explicit judgment of the direction D of blur by the blur direction judgment unit 33. In this case, since the operator Ei that detects the edge perpendicular to the direction of blur is found with reference to the reference table, the direction D of blur is implicitly judged.
In the embodiment described above, the edge detection unit 34 may extract the edges in the direction perpendicular to the direction D of blur, among the edges in the directions and the magnitudes calculated by the face detection unit 31 at the respective pixels in the image P1 for finding the first characteristic quantity. In this manner, redundant edge detection is avoided, which leads to improvement in processing efficiency.
In the embodiment described above, the weight coefficients λd1, λd2, λd3, . . . , λdJ contribute only to the direction D of blur. However, in the case where the weight coefficients λd1, λd2, λd3, . . . , λdJ also contribute to the width of blur in the respective directions corresponding thereto (such as the case where the width becomes larger as the absolute value of each of the weight coefficients λd1, λd2, λd3, . . . , λdJ becomes larger), the width W of blur can be found directly from the weight coefficients λd1, λd2, λd3, . . . , λdJ without the processing by the edge detection unit 34. In this manner, the configuration becomes simpler and the processing becomes more efficient.
More specifically, a block diagram in FIG. 10 shows the configuration in this case. As shown in FIG. 10, a blur correction unit 33′ obtains the operator Ei for detecting the edge in the direction perpendicular to the direction of blur and the enhancement factor αj representing the degree of edge enhancement according to the magnitude of blur, with reference to a reference table T according to the weight coefficients λd1, λd2, λd3, . . . , λdJ obtained by the parameter acquisition unit 32. FIG. 11 shows a structure of the reference table T and an example of values therein. As shown in FIG. 11, the enhancement factor αj is related to the operator Ei for each of the values of the weight coefficients λd1, λd2, λd3, . . . , λdJ. This relationship is found experimentally and statistically in advance. Since the blur correction unit 33′ finds the operator Ei for detecting the edge perpendicular to the direction of blur with reference to the reference table, the blur correction unit 33′ implicitly detects the direction D of blur.
Before the reference table is referred, a combined parameter λD may be found as a linear combination of the weight coefficients as shown by Equation (8) below wherein λi is a coefficient representing a rate of contribution of the ith principal component corresponding to the weight coefficient λi to the direction and the width of blur:




  
    
      
        
          λ
          D
        
        =
        
          
            ∑
            
              i
              =
              1
            
            J
          
          ⁢
          
            
              β
              i
            
            ⁢
            
              λ
              i
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




In this case, the reference table T needs to relate only the combined parameter λD the operator Ei, and the enhancement factor αi. Therefore, the content of setting is simplified.

In addition, a function using the weight coefficients λd1, λd2, λd3, . . . , λdJ as input and the operator Ei and the enhancement factor αi as output may be defined so that the operator Ei and the enhancement factor αi can be found based on the function.
In the embodiment described above, the known selective sharpening is used as the blur correcting process to be carried out on the entire image P1. However, the blur correction may be carried out on the face region P1f by changing the values of the weight coefficients λd1, λd2, λd3, . . . , λdJ to values corresponding to a state with no blur (such as 0).
In the embodiment described above, the mathematical model M is unique. However, a plurality of mathematical models Mi (i=1, 2, . . . ) may be generated for respective properties such as race, age, and gender, for example. FIG. 12 is a block diagram showing details of blur direction judging process and the blur correcting process in this case. As shown in FIG. 12, a property acquisition unit 37 and a model selection unit 38 are added, which is different from the embodiment shown in FIG. 4. The property acquisition unit 37 obtains property information AK of a subject in the image P1. The model selection unit 38 selects a mathematical model MK generated only from sample images representing subjects having a property represented by the property information AK.
The mathematical models Mi have been generated based on the method described above (see FIG. 5), only from sample images representing subjects of the same race, age, and gender, for example. The mathematical models Mi are stored by being related to property information Ai representing each of the properties that is common among the samples used for the model generation.
The property acquisition unit 37 may obtain the property information AK by judging the property of the subject through execution of known recognition processing (such as processing described in Japanese Unexamined Patent Publication No. 11 (1999)-175724) on the image P1. Alternatively, the property of the subject may be recorded at the time of photography as accompanying information of the image P1 in a header or the like so that the recorded information can be obtained. The property of the subject may be inferred from accompanying information. In the case where GPS information representing a photography location is available, the country or a region corresponding to the GPS information can be identified, for example. Therefore, the race of the subject can be inferred to some degree. By paying attention to this fact, a reference table relating GPS information to information on race may be generated in advance. By inputting the image P1 obtained by a digital camera that obtains the GPS information at the time of photography and records the GPS information in a header of the image P1 (such as a digital camera described in Japanese Unexamined Patent Publication No. 2004-153428), the GPS information recorded in the header of the image data P1 is obtained. The information on race related to the GPS information may be inferred as the race of the subject when the reference table is referred to according to the GPS information.
The model selection unit 38 obtains the mathematical model MK related to the property information AK obtained by the property acquisition unit 37, and the parameter acquisition unit 32 fits the mathematical model MK to the face region P1f in the image P1.
As has been described above, in the case where the mathematical models Mi corresponding to the properties have been prepared, if the model selection unit 38 selects the mathematical model MK related to the property information AK obtained by the property acquisition unit 37 and if the parameter acquisition unit 32 fits the selected mathematical model MK to the face region P1f, the mathematical model MK does not have eigenvectors contributing to variations in face shape and luminance caused by difference in the property represented by the property information AK. Therefore, the face region P1f can be represented only by eigenvectors representing factors determining the face shape and luminance other than the factor representing the property. Consequently, processing accuracy improves, and the image can be obtained in higher quality.
From a viewpoint of improvement in processing accuracy, it is preferable for the mathematical models for respective properties to be specified further so that a mathematical model for each individual as a subject can be generated. In this case, the image P1 needs to be related to information identifying each individual.
In the embodiment described above, the mathematical models are installed in the digital photograph printer in advance. However, from a viewpoint of processing accuracy improvement, it is preferable for mathematical models for different human races to be prepared so that which of the mathematical models is to be installed can be changed according to a country or a region to which the digital photograph printer is going to be shipped.
The function for generating the mathematical model may be installed in the digital photograph printer. More specifically, a program for causing the arithmetic and control unit 50 to execute the processing described by the flow chart in FIG. 5 is installed therein. In addition, a default mathematical model may be installed at the time of shipment thereof. In this case, the mathematical model may be customized based on images input to the digital photograph printer. Alternatively, a new model different from the default model may be generated. This is especially effective in the case where the mathematical models for respective individuals are generated.
In the embodiment described above, the individual face image is represented by the weight coefficients bi and λi for the face shape and the luminance. However, variation in the face shape is correlated to variation in the luminance. Therefore, a new appearance parameter c can be obtained for controlling both the face shape and the luminance as shown by Equations (9) and (10) below, through further execution of principal component analysis on a vector (bi, b2, . . . , bi, . . . , λ1, λ2, . . . , λi, . . . ) combining the weight coefficients bi and λi:

S=S0+Qsc  (9)

A=A0+QAc  (10)

A difference from the mean face shape can be represented by the appearance parameter c and a vector Qs, and a difference from the mean luminance can be represented by the appearance parameter c and a vector QA.
In the case where this model is used, the parameter acquisition unit 32 finds the face luminance in the mean face shape based on Equation (10) above while changing a value of the appearance parameter c. Thereafter, the face image is reconstructed by conversion from the mean face shape according to Equation (9) above, and the value of the appearance parameter c causing a difference between the reconstructed face image and the face region P1f to be minimal is found.
Another embodiment of the present invention can be installation of the blur direction judging process and the blur correcting process in a digital camera. FIG. 13 shows the configuration of such a digital camera. As shown in FIG. 13, the digital camera has an imaging unit 71, an A/D conversion unit 72, an image processing unit 73, a compression/decompression unit 74, a flash unit 75, an operation unit 76, a media recording unit 77, a display unit 78, a control unit 70, and an internal memory 79. The imaging unit 71 comprises a lens, an iris, a shutter, a CCD, and the like, and photographs a subject. The A/D conversion unit 72 obtains digital image data P0 by digitizing an analog signal represented by charges stored in the CCD of the imaging unit 71. The image processing unit 73 carries out various kinds of image processing on image data such as the image data P0. The compression/decompression unit 74 carries out compression processing on image data to be stored in a memory card, and carries out decompression processing on image data read from a memory card in a compressed form. The flash unit 75 comprises a flash and the like, and carries out flash emission. The operation unit 76 comprises various kinds of operation buttons, and is used for setting a photography condition, an image processing condition, and the like. The media recording unit 77 is used as an interface with a memory card in which image data are stored. The display unit 78 comprises a liquid crystal display (hereinafter referred to as the LCD) and the like, and is used for displaying a through image, a photographed image, various setting menus, and the like. The control unit 70 controls processing carried out by each of the units. The internal memory 79 stores a control program, image data, and the like.
The functions of the image input means 1 in FIG. 2 are realized by the imaging unit 71 and the A/D conversion unit 72. Likewise, the functions of the image correction means 2 are realized by the image processing unit 73 while the functions of the image manipulation means 3 are realized by the image processing unit 73, the operation unit 76, and the display unit 78. The functions of the image output means 4 are realized by the media recording unit 77. All of the functions described above are realized under control of the control unit 70, by using the internal memory 79 in addition.
Operation of the digital camera and a flow of processing therein are described next.
The imaging unit 71 causes light entering the lens from a subject to form an image on a photoelectric surface of the CCD when a photographer fully presses a shutter button. After photoelectric conversion, the imaging unit 71 outputs an analog image signal, and the A/D conversion unit 72 converts the analog image signal output from the imaging unit 71 to a digital image signal. The A/D conversion unit 72 then outputs the digital image signal as the digital image data P0. In this manner, the imaging unit 71 and the A/D conversion unit 72 function as the image input means 1.
Thereafter, the image processing unit 73 carries out white balance adjustment processing, gradation correction processing, density correction processing, color correction processing, and sharpness processing, and outputs corrected image data P1. In this manner, the image processing unit 73 functions as the image correction means 2.
The image P1 is displayed on the LCD of the display unit 78. As a manner of this display can be used display of thumbnail images as shown in FIG. 3A. While operating the operation buttons of the operation unit 76, the photographer selects and enlarges one of the images to be processed, and carries out selection from a menu for further manual image correction or manipulation. In the case where “Blur Correction” is selected at this stage, the control unit 70 starts a program for blur direction judgment and blur correction stored in the internal memory 79, and causes the image processing unit 73 to carry out the blur direction judging process and the blur correcting process (see FIG. 4) using the mathematical model M stored in advance in the internal memory 79, as has been described above. Processed image data P2 are then output. In this manner, the functions of the image manipulation means 3 are realized.
The compression/decompression unit 74 carries out compression processing on the image data P2 according to a compression format such as JPEG, and the compressed image data are written via the media recording unit 77 in a memory card inserted in the digital camera. In this manner, the functions of the image output means 4 are realized.
By installing the blur direction judging process and the blur correcting process of the present invention as the image processing function of the digital camera, the same effect as in the case of the digital photograph printer can be obtained.
The manual correction and manipulation may be carried out on the image having been stored in the memory card. More specifically, the compression/decompression unit 74 decompresses the image data stored in the memory card, and the image after the decompression is displayed on the LCD of the display unit 78. The photographer selects desired image processing as has been described above, and the image processing unit 73 carries out the selected image processing.
Furthermore, the mathematical models for respective properties of subjects described by FIG. 12 may be installed in the digital camera. In addition, the processing for generating the mathematical model described by FIG. 5 may be installed therein. A person as a subject of photography is often fixed to some degree for each digital camera. Therefore, if the mathematical model is generated for the face of each individual as a frequent subject of photography with the digital camera, the model can be generated without variation of individual difference in face. Consequently, the blur direction judging process and the blur correcting process can be carried out with extremely high accuracy for an image including the face of the person.
A program for causing a personal computer or the like to carry out the blur direction judging process and the blur correcting process of the present invention may be incorporated with image editing software. In this manner, a user can use the blur direction judging process and the blur correcting process of the present invention as an option of image editing and manipulation on his/her personal computer, by installation of the software from a recording medium such as a CD-ROM to the personal computer, or by installation of the software through downloading of the software from a predetermined Web site on the Internet.

## Claims

1. A blur direction judgment method comprising the steps of:
obtaining one or more weighting parameters for one or more statistical characteristic quantities representing a direction of blur in a predetermined structure in an input image by fitting a model representing the structure to the structure in the input image, the model having been obtained by carrying out predetermined statistical processing on a plurality of images representing the predetermined structure in different conditions of blur, and the model representing the structure by one or more statistical characteristic quantities including the statistical characteristic quantity or quantities representing the direction of blur and by one or more weighting parameters for weighting the statistical characteristic quantity or quantities according to an individual characteristic of the structure; and
judging the direction of blur represented in the input image according to a value or values of the weighting parameter or parameters having been obtained.

2. A blur direction judgment apparatus comprising:
parameter acquisition means for obtaining one or more weighting parameters for one or more statistical characteristic quantities representing a direction of blur in a predetermined structure in an input image by fitting a model representing the structure to the structure in the input image, the model having been obtained by carrying out predetermined statistical processing on a plurality of images representing the predetermined structure in different conditions of blur, and the model representing the structure by one or more statistical characteristic quantities including the statistical characteristic quantity or quantities representing the direction of blur and by one or more weighting parameters for weighting the statistical characteristic quantity or quantities according to an individual characteristic of the structure; and
blur direction judgment means for judging the direction of blur represented in the input image according to a value or values of the weighting parameter or parameters having been obtained by the parameter acquisition means.

3. The blur direction judgment apparatus according to claim 2 further comprising blur width acquisition means for finding a width of blur in the direction of blur having been judged by the blur direction judgment means.

4. The blur direction judgment apparatus according to claim 3, wherein the blur width acquisition means finds the width of blur based on the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur.

5. The blur direction judgment apparatus according to claim 4 further comprising blur correction means for correcting the blur represented in the input image based on the direction of blur having been judged by the blur direction judgment means and the width of blur obtained by the blur width acquisition means.

6. The blur direction judgment apparatus according to claim 5, wherein the blur correction means corrects the blur in the structure in the input image by changing the value or values of the weighting parameter or parameters for the statistical characteristic quantity or quantities representing the direction of blur.

7. The blur direction judgment apparatus according to claim 2 further comprising detection means for detecting the structure in the input image, wherein
the parameter acquisition means obtains the weighting parameter or parameters by fitting the model to the structure having been detected.

8. The blur direction judgment apparatus according to claim 2 further comprising selection means for obtaining at least one property of the structure in the input image and for selecting the model corresponding to the obtained property from a plurality of the models representing the structure for respective properties of the predetermined structure, wherein
the parameter acquisition means obtains the weighting parameter or parameters by fitting the selected model to the structure in the input image.

9. The blur direction judgment apparatus according to claim 2, wherein the predetermined structure is a human face.

10. The blur direction judgment apparatus according to claim 2, wherein the model and the fitting of the model are realized by a method of Active Appearance Model.

11. A non-transitory computer-readable recording medium storing a blur direction judgment program causing a computer to function as:
parameter acquisition means for obtaining one or more weighting parameters for one or more statistical characteristic quantities representing a direction of blur in a predetermined structure in an input image by fitting a model representing the structure to the structure in the input image, the model having been obtained by carrying out predetermined statistical processing on a plurality of images representing the predetermined structure in different conditions of blur, and the model representing the structure by one or more statistical characteristic quantities including the statistical characteristic quantity or quantities representing the direction of blur and by one or more weighting parameters for weighting the statistical characteristic quantity or quantities according to an individual characteristic of the structure; and
blur direction judgment means for judging the direction of blur represented in the input image according to a value or values of the weighting parameter or parameters having been obtained by the parameter acquisition means.